全国大学生数学建模大赛国家一等奖优秀论文系泊系统的设计.docx
《全国大学生数学建模大赛国家一等奖优秀论文系泊系统的设计.docx》由会员分享,可在线阅读,更多相关《全国大学生数学建模大赛国家一等奖优秀论文系泊系统的设计.docx(59页珍藏版)》请在冰点文库上搜索。
全国大学生数学建模大赛国家一等奖优秀论文系泊系统的设计
系泊系统的设计
摘要
本文详细对系泊系统的各个机构进行了力学分析,针对系泊系统的要求,建立优化模型,求解系泊系统在多种环境下的最优解,使得浮标游动范围,吃水程度和钢桶倾斜角度尽可能的小。
针对问题一,本文对系泊系统的受力及力矩进行了分析,基于浮标倾斜的考虑,得到了平衡状态下关于受力平衡及力矩平衡的方程组。
由于方程组数量较多及相互影响的特点,直接求解十分困难。
因此我们考虑以浮标两边的浸水长度hi,h2为变量,利用搜索算法对方程组进行求解,并得到相应的结果。
如当风速为12m/s时,钢桶的倾斜角度°,从上到下钢管的倾斜角度分别为°、°、°、°浮标吃水深度,浮标游动区域半径。
针对问题二,首先将风速为36m/s的情况代入问题一建立的模型中,但是得到的结果不满足题目所给定的要求。
则考虑在重物球质量一定的条件下,以浮标的吃水深度和游动区域及钢桶的倾斜角为目标,建立了一个单决策变量的多目标最优系泊模型,相比于问题一,此问的变量更多,更加难于求解,故考虑将多目标转化成单目标的问题进行求解,并继续使用搜索法对问题进行求解。
最后找到了三组可行解,其中最优解是重力球的质量为2102kg.
针对问题三,本文中有三个决策变量以及三个变系数,相比于前两问,无论是计算量还是计算维数,难度更大。
为了求解该问,建立了一个多决策变量的多目标变系数的最优系泊系统模型,为了简便运算,我们建立了变步长的搜索算法,并最终求解得到结果,得到的一组解为:
选用了III型号的锚链,重物球质量为2800kg,锚链长度为。
针对论文的实际情况,对论文的优缺点做了评价,文章最后还给出了其他的改进方向,以用于指导实际应用。
关键词:
系泊系统设计;力的平移定理;多目标;优化模型;搜索算法
1.问题的重述
一个由浮标系统、系泊系统和水声通讯系统组成的近浅海观测网的传输节点。
可以简化看作是一个浮于海平面的圆柱体浮标通过四节钢管链接装有通讯设备的钢桶,钢桶再通过锚链链接一个可以移动的锚,锚沉在海床上。
为了不让锚被拖行要求锚链末端与锚的链接处的切线方向与海床的夹角不超过16度,并且要保证通讯设备的工作效果,钢桶与电焊锚链链接出可悬挂重物体,使得钢桶的倾斜角度(钢桶与竖直线的夹角)不超过5度。
需要建立模型讨论并解决下列的几个问题:
1.若海水静止,在选用II型电焊锚链,重物球的质量为1200kg,布放水
深为18m海床平坦,海水密度为x103kg/m3的条件下,分别计算海面风速为12m/s
和24m/s时钢桶和各节钢管的倾斜角度、锚链形状、浮标的吃水深度和游动区域。
2.在问题1的假设下,计算海面风速为36m/s时钢桶和各节钢管的倾斜角度、
锚
链形状和浮标的游动区域。
并调节重物球的质量,使得钢桶的倾斜角度不超过5
度,锚链在锚点与海床的夹角不超过16度。
3.若布放海域的实测水深介于16m-20m之间。
布放点的海水速度最大可达到
s、风速最大可达到36m/s0给出考虑风力、水流力和水深情况下的系泊系统设计,
分析不同情况下钢桶、钢管的倾斜角度、锚链形状、浮标的吃水深度和游动区域
2.模型的假设
(1)假设两钢管用铰链链接在一起,可以自由转动。
(2)假设风是沿平行海平面的方向吹来。
(3)假设忽略锚链和重物球的浮力影响。
(4)假设忽略由于其它原因(如温度,湿度,受力等)而对各个系统产生的形变。
(5)假设浮标不会出现因为风力过大而被吹翻的情况。
(6)假设各个系统之间连接处的长度忽略不计。
3.符号说明
Fwind:
浮标所受水平方向的风力大小;
FbuoyantA:
浮标所受浮力大小;
FbuoyantB:
每节钢管所受浮力大小;
FbuoyantC:
钢桶所受浮力大小;
:
第一节钢管对浮标拉力的大小(或浮标对第一节钢管拉力的大小);
Fb:
:
第i节钢管下端所受拉力的大小(i1,2,3,4);
Fc(Fd°):
钢桶所受锚链的拉力大小(或锚链上方第一段所受钢桶的拉力大小);
Fq:
锚链第j段(自上而下)下端所受拉力大小(j12|[|,N,其中N为张紧状态锚链的段数);
Ff:
锚与海底间的摩擦力;
Fmax:
锚与海底间的最大静摩擦力;
FwaterA:
浮标所受的水流力大小;
FwaterBi:
第:
节钢管所受水流力大小(i1,2,3,4);
FwaterC:
钢桶所受水流力大小;
:
锚与海底间的最大静摩擦系数;
Fn:
海底对锚的支持力;
Ga:
浮标的重力;
Gb:
每节钢管的重力;
Gc:
表示钢桶的重力;
Gd:
锚链每段链环的重力;
Ge:
锚的重力;
Bo:
第一节钢管对浮标拉力与水平方向的夹角;
Bi:
第i节钢管下端所受拉力与水平方向的夹角(i1,2,3,4);
c(Do):
钢桶所受锚链的拉力与水平方向的夹角;
Dj.锚链第j段(自上而下)下端所受拉力与水平方向的夹角(j1,2,|||,N);
A:
浮标的水平倾斜角度;
b:
第i节钢管的水平倾斜角度(i123,4);
c:
钢桶的水平倾斜角度;
hi,h2:
圆柱体浮标两端的浸水长度,其中h,h2.
h:
浮标的吃水深度;
Vwind:
水平方向的风速;
Swind:
浮标在风向法平面的投影面积;
VoverflowA:
浮标的排水体积;
Vwater:
水流速度大小;
mA:
浮标的质量;
mB:
每节钢管的质量;
mC:
钢桶系统的质量;
m:
重物球的质量;
:
海水密度;
D:
锚链的线密度;
CDN:
法向阻力系数;
CDN:
切向阻力系数;
rA:
浮标的底面半径;
rB:
每节钢管的底面半径;
LB:
每节钢管的长度;
rC:
钢桶的底面半径;
LC:
钢桶的长度;
LD:
锚链每段链环的长度;
H:
海水深度;
s:
锚链某点至锚链上端顶点的弧长;
SA:
浮标浸入水下部分在水流力法平面的投影面积;
SBi:
第i节钢管在水流力法平面的投影面积(i1,2,3,4);
SC:
钢桶在水流力法平面的投影面积;
T:
锚链某点处所受拉力的大小;
:
锚链某点处所受拉力与水平方向的夹角;
XA:
平衡状态下浮标的水平距离;
XB:
平衡状态下钢管的水平距离;
XC:
平衡状态下钢桶的水平距离;
XD:
平衡状态下锚链的水平距离;
X:
平衡状态下整体的水平距离;
yD:
锚链竖直方向上的高度;
max:
锚不滑动时,锚链在锚点与海床的最大夹角;
D:
单位长度上锚链所受的法向水流力大小;
D:
单位长度上锚链所受的切向水流力大小;
4.模型的建立与求解
问题的分析
系统平衡的两种状态一一张紧状态和松弛状态及游动区域的理解
钢管、
本文考虑的在风力及水流力作用下的系泊系统的设计问题,首先我们考虑一种简单情形,即海水静止且风向水平情况,此时整个系泊系统可能会因风力移动达到一个平衡的状态,由于风向恒定我们可以认为最终的平衡状态中浮标、钢桶、锚链和锚大致位于同一个平面内,如图1所示。
图1传输节点示意图
在平衡状态下,各系统所受力及力矩分别达到平衡,此时浮标会与锚形成一定的水平距离。
此时的平衡状态可能分为两种:
一种是张紧状态,即锚链的每一段都受到了前后端的拉力从而使整个系统平衡;另一种是松弛状态,即锚链上方的一部分存在拉力,而下方部分散落在水底呈无拉力状态达到平衡。
在张紧状态下达到平衡时,浮标到达某个定点,此时该点与锚得水平距离为
半径以锚和浮标的水平距离为半径、锚所在位置为圆心做一个圆面,该区域即浮
标的游动区域,如图2所示。
但在松弛状态下由于松弛部分锚链的可伸缩性,浮标的位置会处于一定的区域,当松弛部分锚链长度最大时浮标位于最远位置,
以此时浮标与锚的水平距离为半径、锚所在位置为圆心的圆形区域即其游动区域。
图2浮标游动区域俯视示意图
浮标的倾斜问题及对吃水深度的理解
在对组成各个系统的零件受力分析的过程中,会出现无法将物体看作是均匀的情况,这样一来,物体的各力的受力点将不再是质点,为此应尽可能将受力点移动到质心上。
这时,根据力的平移定理,在将各力平移到质心上去的过程中会产生力矩,且对于一个平衡的物体其力矩和应为零。
在力矩平衡的作用下,我们发现平衡状态下浮标会产生一定的倾斜角度。
3所示。
假设在浮标的底面与海平面平行,即没有倾角时,其受力分析如图
图3浮标垂直海平面的受力分析图
且方
由图可知除浮力FbuoyantA外,风力卩⑷血和钢管拉力Fb°都会产生非零力矩,向相同,导致力矩和不为零,从而会产生旋转,不会达到平衡状态,我们认为在平衡状态下,浮标底面相对海平面存在一个倾斜角度a。
h2,
吃水深度为浮标底部到海面的距离。
由于浮标倾斜,导致浮标各点侵入海水的长度不同,如图4所示。
假设浮标两边侵入水中的长度分别为h1,h2,其中h则吃水深度即两边浸水长度的平均值。
2
图4倾斜浮标的浸水示意图
问题1模型的建立
根据对问题的分析发现,问题旨在讨论在不同情形达到平衡状态下,分析各
钢管的倾斜角度、锚链形状、浮标的方程,通过已知的信息,求得系统各部分所
受拉力大小及其角度情况和各部分本身的倾斜角度。
浮标的力平衡方程
通过浮标的受力分析可得,浮标受到竖直向上的浮力FbuoyantA,水平向右的风
力Fwind,重力Ga以及钢管对其的拉力Fb°,其受力分析图如图5所示:
图5浮标的受力分析图
根据平面力系平衡原理,可得浮标的力平衡方程如下:
在水平方向上,有
FB0COSb0Fwind,
(1)
其中,Bo为钢管对浮标拉力与水平方向的夹角,Fwind为水平风力的大小,有
如下的近似公式(S为物体在风法平面的投影面积)
2
Fwind0.625S(Vwind);
(2)
在竖直方向上,有
FBosinBo
mAg
FbuoyantA,
(3)
其中,有
而
FbuoyantA
gVoverflowA,
(4)
VoverflowA
r2,1
rAh—
A2
h)2rA2(hh2),
(5)
在这儿,mA为浮标质量,hnb分别表示浮标两边的浸水长度(见图4),rA为浮标的底面半径,VoverflowA为浮标的排水体积,为海水密度,g为重力加速度。
浮标在风向法平面的投影面积公式
在计算风力Fwind时,需用到浮标在风向法平面的投影面积S,由于浮标倾斜,
该投影区域并非一标准的矩形,而是由两个部分构成,其中的下方部分为矩形,而上方部分为倾斜圆在水平面上投影的一半,即半椭圆的一半,最终得到浮标在风向法平面的投影面积如图6所示。
图6浮标在风向法平面的投影面积示意图
通过平面几何关系可得,第一部分矩形的面积S满足:
S=2rA[LA0.5(0h2)]麻hTW
第二部分半椭圆的面积S2满足:
钢管的平衡方程
通过分析可得,每节钢管的受力分析图如图7所示
图7第i节钢管的受力分析图
根据平面力系平衡原理及力的平移定理(力矩平衡公式),有如下方程:
水平方向力平衡方程
FBi
cosBi
Fbi
cosBi(i1,2,3,4),
(8)
竖直方向力平衡方程
FBisinBi
mBg
FBii
sinBi1FbuoyantB(i1,2,3,4),
(9)
力矩平衡方程
1
FBi2Sin(Bi
Bi)Lb
FBii
1
2si门(耳1bJLbQ123,4).
(10)
在这儿,FbuoyantB为每节钢管所受浮力大小,满足
(11)
而FBi为第i节钢管下端所受拉力大小,mB为每节钢管的质量,g为重力加速度,B是第i节钢管与水平方向的夹角(钢管的倾斜角度),Bi是拉力FBi与水平方向的夹角,rB为钢管的底面半径,Lb为钢管的长度。
钢桶系统的平衡方程
钢桶的受力分析如图8所示
图8钢桶的受力分析图
根据平面力系平衡原理及力的平移定理(力矩平衡公式),有如下方程:
水平方向力平衡方程
竖直方向力平衡方程
力矩平衡方程
在这儿,FbuoyantC为钢桶系统的浮力,有
FbuoyantC
grC2LC,
(15)
而Fb4表示钢管对钢桶拉力的大小,b4表示钢管对钢桶拉力与水平方向的夹
角,Fc表示钢桶受下方锚链拉力的大小,C表示Fc与水平方向的夹角,me表示
钢桶的质量,m表示重物球的质量,re表示钢桶的底面半径,Le表示钢桶的长度。
锚链系统的平衡方程
由于锚链无档普通链环,我们可以将其看作无弹性悬垂线,这样可以将锚链
视作柔性的。
因此可用微元法的思想分析其受力平衡状态。
考虑锚链其中自上而下一小段弧长(S到sds)上的受力情况,如图9所示。
图9锚链微元的受力分析图
图中ds表示取的任意一段小弧长,T表示上段锚链对所取小弧长的拉力,而
TdT表示下段锚链对小弧长的拉力;是拉力T与水平方向的夹角;d是拉
力TdT与水平方向的夹角。
根据这一弧长微元上的受力平衡,可得出以下方程:
水平方向上,有
竖直方向上,有
Tsin(TdT)sin(
其中D表示单位长度锚链的质量。
将其展幵可得
d)
Ddsg
(17)
TcosTcoscosd
dTcoscosd
Tsin
sind
dTsin
sind(18)sindDdsg
TsinTsincosd
dTsincosd
Tcos
sind
dTcos
由泰勒公式sind
d(d)3
3!
III可知,
当角度变化d
很小时,忽略高阶无
穷小的影响,故sind近似等于d,cosd近似等于1,从而可得如下的近似方
程:
0dTcosTsinddTsind(伯)
DdsgdTsinTcosddTcosd
忽略高阶无穷小dTd项,(其中dTd都为无穷小,故它们的乘积为高阶的
无穷小)对上式化简可得如下一阶非线性微分方程组:
ddg
cos
dsT
dT.
Dgsinds
且有初值条件(0)c,T(0)Fc
锚的平衡方程及平衡状态条件
10所示:
锚的受力直接关系着最终是否达到平衡状态,其受力分析如图
图10锚的受力分析图
(21)
水平方向上,有
FDnCOSDnFf,
在竖直方向上,有
由于判断整个系统平衡的关键即对锚的水平拉力是否超过锚与海底之间的
最大静摩擦力Fmax,即得到平衡状态的判断条件为
又由海水深度为18m,可知
hh24
HsinaLbsin(bJLcsincyD
2i1
其中,H表示海水深度;d^sina表示浮标底面轴心到海平面的距离,h.,0
2
4
分别表示浮标两边浸入水中的长度,A浮标的倾斜角度,LBsin(B)表示4个
i1
钢管在水中的高度和,LB表示每节钢管的长度,Bi表示第i节钢管的倾斜角度;
Lc表示钢桶的长度,c表示钢桶的倾斜角度,yd表示锚链竖直方向的高度。
游动区域的分析
经过分析和求解可以知道在平衡状态下的各个系统水平距离分别为:
浮标的水平距离
钢管的水平距离
(26)
钢桶的水平距离
(27)
XcLccosc,
锚链的水平距离
则整体的水平距离(游动半径)为
XXAXBXCXD(29)
问题一模型的求解
模型的求解思路
通过对系泊系统各部分进行单独受力分析和力矩分析,得到了系统各处拉力大小及其倾角等未知量在平衡状态下的关系,由于方程数量较多,且相互影响,使得直接计算非常复杂。
因此我们首先考虑各方程之间的关系,得到未知量满足一定规律的递推公式,对问题进行系列简化,然后再采用搜索算法进行求解。
方程的简化及递推公式的推导
锚链系统微分方程的转化
由于锚链是无档普通链环,实际中是环环相扣的离散系统,为了计算方便,我
们首
先考虑对锚链系统对应的一阶非线性微分方程组(20式)进行转换[3],由ds与水平
方向及竖直方向上的位移之间的关系
dxcosds,dysinds,(30)
可得链环连接处的受力大小及角度的递推公式
其中,FD0表示链环上方第一段所受钢桶带来拉力的大小,D0表示链环上方
第一段所受钢桶带来拉力与水平方向的角度,由于其与链环对钢桶拉力FC互为反作用力,有Fd0Fc,d0c;Fq表示链环第j段(自上而下)下端所受拉力大小,d表示链环第j段下端所受拉力与水平方向的夹角(j12|||,N),N为链环的段数。
水平力平衡公式的推导及简化
整个系泊系统的每个部分均在两个水平力的作用下平衡,有
FwindFBicosBiFccosCFDjcosDj(i1,2,3,4;j1,2,川,N),(32)
其中风力Fwind由浮标两边的入水深度0山2决定。
竖直方向力平衡公式的推导及简化
由竖直方向上受力平衡可得到如下递推公式
其中浮标的浮力FbuoyantA也由浮标两边的入水深度^,人2决定。
系统各段所受拉力大小及其角度的计算公式
将各拉力竖直方向分力与水平分力相比,可得相应角度正切值tan的计算公
式如
各个部件在竖直方向的高度的计算
(1)找出浮标的竖直方向高度,这里记为Ha(HaQR)。
如图所示:
图11•浮标的竖直方向高度示意图
图中,h,MN,hOP,HaQR。
在直角RMQ和直角NKP中,因为
RMQNPK(同位角相等),所以这两个三角形另外两个角相等,即:
RQMKNP
arctan
KP
i
(35)
NK
所以
HaRQQMcosRQM
h1h2
2
cos(arcta门(_^))
KN
(36)
即:
HrA(h1
h2)
(37)
A7(h1h2)2
(2rA)2
其中,
5表示浮标的半径,为1m。
(2)钢管的竖直方向高度,这里记为Hb:
。
根据浮标的平衡方程、水平方向的方程(Fb,cosB0Fwind)和竖直方向的方
程(FbosinBoGaFbuoyantA)相除得:
即可以求出B),这样再带入浮标水平方向平衡方程则可求出Fb0。
再把B0和
Fb0的值带入钢管的平衡方程,我们可以求出钢管的FBi拉力及其倾斜角度Bi和钢
管的倾斜角Bi,又可以根据钢管的几何关系(如图6所示):
HBHBilBisin(Bi)
其中HB表示四根钢管的总高度,
ii
HBi表示第i节钢管的高度,Bi表示第i节
钢管的倾斜角度
3)钢桶的竖直方向高度,这里记为HC
把上个环节的计算出来的结果带入钢桶系统的平衡方程,同样也可以计算出
钢桶倾斜角度C,锚链对钢桶的拉力Fc及其方向(与X轴的夹角C)0则可以得
出钢桶高度的公式:
HClCsinC(40)
其中lC为钢桶的长度,它是一个常数,大小为1m0
(4)钢桶的竖直方向高度,这里记为HD0
根据已求得的钢桶对锚链的拉力Fc及其方向,把这对初始值带入锚链平衡微
分方程,可解出函数关系0就可以表示出锚链的高度:
D22.05
HD0dssin(d)(41)
c0
其中,ds表示锚链一微段的长度,d表示那一微段的倾斜角度。
故整个系泊系统总高度为
HHaHbHeHd(42)
吃水深度h,和g的值范围的确定
根据题目可知浮标的高度为2m,贝y有:
0h2m,0h,2m,h2为较长
的一条高。
只有当浮标在海水中时,浮标的重力与浮标的浮力二力平衡,即:
Ga
gVoverflowA
(43)
可以计算出此时h1d
0.31mo那么可以进一步确定范围:
0.31
h,h22,
(44)
算法步骤
(1)给定一组浮标吃水深度值[Uh?
],,(得到0.31h,h22),可以依次计算得到H,Dn;
(2)判断如果同时满足|H18,(取较小值,比如)和0Dn16两个条件,则输出结果,并停止计算;如果不满足,则取步长有h.h.»h2h2h2,并转到步骤
(1)o
计算结果及分析
(1)当风速为12ms时,得到的结果如下表所示:
表1.风速为12ms时的求解结果
h1
h2
海水深度
锚链水平倾
角
钢桶的倾角
所用无档链环个
数
O
O
144
钢管1的倾
角
钢管2的倾
角
钢管3的倾
角
钢管4的倾
角
游动区域半
径
浮标的吃水深度
O
O
O
O
得到锚链的形状如图12所示:
图12.风速为12ms时锚链的形状
以上求解得到的结果满足题目的要求,且所用的链环个数为144个,小于链
的总个数,说明达到平衡时锚链并未完全绷紧
(2)当风速为24ms时,得到的结果如下表所示:
表2.风速为24ms时的求解结果
hi
h2
海水深度
锚链水平倾
角
钢桶的倾角
所用无档链环个
数
O
O
185
钢管1的倾
角
钢管2的倾
角
钢管3的倾
角
钢管4的倾
角
游动区域半
径
浮标的吃水深度
O
O
O
O
得到锚链的形状如图13所示:
图13.风速为24ms时锚链的形状
以上求解得到的结果满足题目的要求,且所用的链环个数为185个,小于链
的总个数,说明达到平衡时锚链并未完全绷紧。
问题2的模型分析与建立
风速为36ms时系泊状态及其分析
当风速达到36ms时,若不考虑角度的限制,可通过问题一的递推式及搜索算法求得如下结果:
表3.风速为36ms时的求解结果
hi
h2
海水深度
锚链水平倾
角
钢桶的倾角
所用无档链环个
数
O
O
210
钢管1的倾
角
钢管2的倾
角
钢管3的倾
角
钢管4的倾
角
游动区域半
径
浮标的吃水深度
得到锚链的形状如图14所示:
图14.风速为36ms时锚链的形状
此时,锚链末端与锚的链接处的切线方向与海床的夹角为度,超过了16度,
钢桶的倾斜角度(钢桶与竖直线的夹角)为度,超过5度,已经不能够达到平衡状态,因此需考虑调节重物球质量,使其重新达到平衡状态。
基于重物球质量的定常系统最优系泊模型
当风速达到36ms时,系泊系统已经不再满足角度的要求,从而不再处于平
衡状态,因