数学建模A题系泊系统设计

系泊系统的设计摘要此题要求观测近海观测网的组成，建立模型对其中系泊系统进展设计，在不同风速和水流的情况下确定锚链，重物球，钢管及浮标等的状态，从而使通讯设备的工作效果最正确。求解的具体流程如下：针对问题一，分别对系统中的受力物体在水平方向和竖直方向上的力进展分析，找出锚链对锚无拉力时的临界风速，运用力矩平衡求出钢管与钢桶的倾斜角度。对于锚链，将其等效为悬链线模型，根据风速不同判断锚链的状态，从而求出结果。针对问题二，需要调节重物球的质量，使通讯设备在QUOTE时能够正常工作。为了确定重物球的质量，首先将实际风速与临界风速进展比拟，判断此时系统中各物体的状态，与题目中数据进展比拟。在钢桶倾斜角度到达临界角度时，计算锚链与海床的夹角并于题中数据进展比拟，计算重物球的质量。在浮标完全没入海面时，计算相应条件下重物球的质量，从而确定满足条件的重物球的质量围。针对问题三，要求在不同条件下，求出系泊系统中各物体的状态。以型号I锚链为例，当水流方向与风速方向一样时，系统条件最差，分析在不同水深条件下的系泊系统设计。由题中条件确定系统设计的限制条件，对系统各物体进展受力分析，以使整体结果最小，即可得出最优的系泊系统设计。关键词：悬链线多目标非线性规划一、问题重述近浅海观测网的传输节点由浮标系统、系泊系统和水声通讯系统组成〔如图1所示〕。*型传输节点的浮标系统可简化为底面直径2m、高2m的圆柱体，浮标的质量为1000kg。系泊系统由钢管、钢桶、重物球、电焊锚链和特制的抗拖移锚组成。锚的质量为600kg，锚链选用无档普通链环，近浅海观测网的常用型号及其参数在附表中列出。钢管共4节，每节长度1m，直径为50mm，每节钢管的质量为10kg。要求锚链末端与锚的处的切线方向与海床的夹角不超过16度QUOTE，否则锚会被拖行，致使节点移位丧失。水声通讯系统安装在一个长1m、外径30cm的密封圆柱形钢桶，设备和钢桶总质量为100kg。钢桶上接第4节钢管，下接电焊锚链。钢桶竖直时，水声通讯设备的工作效果最正确。假设钢桶倾斜，则影响设备的工作效果。钢桶的倾斜角度〔钢桶与竖直线的夹角〕超过5度时，设备的工作效果较差。为了控制钢桶的倾斜角度，钢桶与电焊锚链处可悬挂重物球。系泊系统的设计问题就是确定锚链的型号、长度和重物球的质量，使得浮标的吃水深度和游动区域及钢桶的倾斜角度尽可能小。问题1*型传输节点选用II型电焊锚链22.05m，选用的重物球的质量为1200kg。现将该型传输节点布放在水深18m、海床平坦、海水密度为1.025×103kg/m3的海域。假设海水静止，分别计算海面风速为12m/s和24m/s时钢桶和各节钢管的倾斜角度、锚链形状、浮标的吃水深度和游动区域。问题2在问题1的假设下，计算海面风速为36m/s时钢桶和各节钢管的倾斜角度、锚链形状和浮标的游动区域。请调节重物球的质量，使得钢桶的倾斜角度不超过5度，锚链在锚点与海床的夹角不超过16度。问题3由于潮汐等因素的影响，布放海域的实测水深介于16m~20m之间。布放点的海水速度最大可到达1.5m/s、风速最大可到达36m/s。请给出考虑风力、水流力和水深情况下的系泊系统设计，分析不同情况下钢桶、钢管的倾斜角度、锚链形状、浮标的吃水深度和游动区域。二、模型假设1.不考虑流体对锚链的作用，忽略锚链本身的伸长，锚链沿长度均匀分布；2.假设风是二维的，只存在平行于水平面的风速，不存在垂直方向上的分量；三、符号说明符号含义浮标的浮力钢管的浮力钢桶的浮力重物球的浮力近海风荷载浮标近海水流力钢桶近海水流力重物球近海水流力浮标重力钢管重力钢桶重力重物球重力锚链重力浮标游动区域h浮标吃水深度四、问题分析4.1问题一分析问题一对于系泊系统的受力物体进展分析，分别对系统中的受力物体在水平方向和竖直方向上的力进展分析，找出锚链对锚无拉力时的临界风速，运用力矩平衡求出钢管与钢桶的倾斜角度。对于锚链，将其等效为悬链线模型，根据风速不同判断锚链的状态，从而求出结果。4.2问题二分析问题二需要调节重物球的质量，使通讯设备在QUOTE时能够正常工作。为了确定重物球的质量，首先将实际风速与临界风速进展比拟，判断此时系统中各物体的状态，并与数据进展比拟。在钢桶倾斜角度到达临界角度时，计算锚链与海床的夹角并与题目要求进展比拟，计算重物球的质量。在浮标完全没入海面时，计算相应条件下重物球的质量，综合确定满足条件的重物球的质量围。问题三分析问题三要求在不同条件下，求出系泊系统中各物体的状态。以型号I锚链为例，当水流方向与风速方向一样时，系统条件最差，分析在不同水深条件下的系泊系统设计。由题中条件确定系统设计的限制条件，对系统各物体进展受力分析，以使整体结果最小，即可得出最优的系泊系统设计。五、模型建立与求解5.1问题一模型问题分析问题一,系泊系统整体受力平衡，浮标受到恒定风力时，分别对系统中的受力物体在水平方向和竖直方向上的力进展分析，找出锚链对锚无拉力时的临界风速，运用力矩平衡求出钢管与钢桶的倾斜角度。对于锚链，将其等效为悬链线模型，根据风速不同判断锚链的状态，最后利用力矩平衡求出钢管与钢桶的倾斜角度,从而求出结果。模型建立Step1.系泊系统受力分析对于浮标，它受到水的浮力，自身重力，风载荷及第一根钢管对浮标的力。此题中海水静止，故所有近海水流力为0。将钢管对浮标的力在分别水平方向和竖直方向上进展分解，具体受力如下图，由浮标所受的力平衡得到：图浮标受力分析图式中，QUOTE、QUOTE分别为钢管对浮标的力再水平方向和竖直方向的分立。对于钢管，将第一节钢管对浮标的力设为QUOTE，后续第节钢管对第k节钢管的力为QUOTE，将钢管的力进展分解，钢管的受力情况如下图，因此：图钢管受力分析图式中，QUOTE为钢管收到的浮力，由于钢管体积较小，在水中钢管的浮力与重力相比很小，可忽略不计；QUOTE，QUOTE分别为第k跟钢管对前一根钢管在水平方向和竖直方向的分力。对于钢桶，它受到第四根钢管的力，水的浮力，自身重力以及锚链的拉力和重物球的拉力，具体受力如下图，则：图钢桶受力分析图式中，QUOTE和QUOTE分别表示锚链的拉力在水平方向和竖直方向上的分力。对于系泊系统，由浮标、钢管和钢桶得水平受力分析可得：其中，S为物体在风向法平面的投影面积，v为风速。由竖直方向受力分析得：设浮标的吃水深度为h，则：其中，QUOTE为海水密度，题中为QUOTE；QUOTE为重力加速度，此题取QUOTE；QUOTE为浮标底面积半径，此题中为；QUOTE为浮标质量，为；为每节钢管的质量，为；QUOTE为重物球的质量，此题中为，QUOTE为钢桶底面半径，QUOTE为钢桶的长度，钢管长度与钢桶相等。将重物球当作铁球处理，则重物球的密度为QUOTE，由题中条件可知，可以递推出：Step2.悬链线模型对于锚链，假设其质量分布均匀，可以将锚链作为悬链线处理，从而做出以下分析,由静力学平衡条件可知，在坐标系中，锚链水平分立和垂直分力的代数和为0，可以得到：图锚链受力分析图式中，QUOTE表示锚链对钢桶的拉力，QUOTE、QUOTE分别表示锚链的水平分力与竖直分力，QUOTE表示锚链的重力，QUOTE为锚链任一点与水平方向的夹角。因此，可推导出：QUOTE表示锚链单位长度的质量，即线密度；QUOTE为锚链在水中未触碰海底的长度。由曲线几何关系和力学关系及悬链线模型可求得，导线任一点的斜率为：式中,*、y分别表示锚链在水平方向和竖直方向的投影长度。对两式分别进展积分，并用临界点数据(QUOTE)来确定积分常数，则可得，进一步可求得方程：其中Step3.力矩分析平衡对于钢管，由于钢管的力矩平衡，可以得到：其中，QUOTE表示分别表示第k根钢管的倾斜角度，为水流对钢管的力，则：即为：式中，QUOTE分别是钢管的倾斜角度〔k取1，2，3，4〕，根据题中的数值，可以得到以下公式：代入数值后可求出钢管相应的倾斜角度。对于钢桶，由钢桶的力矩平衡，可以得到：其中，QUOTE表示分别表示钢桶的倾斜角度，于是：带入数值得：式中，QUOTE是钢桶的倾斜角度，由题中给出的数据可以计算出钢桶的倾斜角度。模型求解求解不同风速条件下系泊系统的状态，需要先求出系统的临界条件，即锚链与海床相切时的风速，然后将此临界风速与题目中所给风速进展比拟，判断锚链的状态，最终求出不同风速下系泊系统的状态。Step1.临界条件当系统处于临界状态时，锚链与海床底部相切，此时，锚链全部抬起且对锚无拉力，QUOTE。取锚链以上的局部进展受力分析，设锚链以上局部竖直高度为QUOTE，由受力分析模型得：由模型中的公式可以得到，此式是关于临界风速QUOTE的表达式。设锚链的竖直高度为y，可以得到：由题中条件可知，水深为QUOTE,则：取锚以上局部进展分析，由系统受力平衡可得，锚以上局部重力与浮力相等，由此可得：将题中各数值代入公式，编程可求得QUOTE。Step2.风速为QUOTE时由于风速QUOTE小于临界速度QUOTE的值，所以当风速为QUOTE时，锚链与海底接触。此时，竖直方向上的受力QUOTE，可以得到：由悬链线模型可以得到，解得：因此：由题中条件可知，水深为QUOTE,则：由锚链与海底深度的关系，编程求得浮标吃水深度h=0.688994。浮标的游动区域由锚链在水平方向的投影长度，以及钢管和钢桶在水方向的投影长度能共同决定，设浮标游动区域为QUOTE则可以得到：由锚链在水平方向的投影长度及钢管、钢桶的倾斜角度可以求得，浮标游动区域大小QUOTE。由钢管与钢桶的受力分析可以求得，四节钢管与钢桶的倾斜角度分别如表所示。表12m/s时钢管及钢桶倾斜角度物体倾斜角度第一节钢管1.136638第二节钢管1.146308第三节钢管1.155913第四节钢管1.165680钢桶1.109907Step3.风速为QUOTE时由于风速QUOTE大于临界速度QUOTE的值，所以当风速为QUOTE时，锚链完全倾斜。此时，锚链倾斜的长度QUOTE。由临界条件中锚以上局部重力与浮力相等，即：代入数值得：可求得QUOTE的值。由浮标，钢桶，钢管水平受力分析可得：由悬链线模型可得：由题中条件可知，水深为QUOTE,则：由锚链与海底深度的关系，编程求得浮标吃水深度h=0.703525。浮标的游动区域由锚链在水平方向的投影长度，以及钢管和钢桶在水方向的投影长度能共同决定，设浮标游动区域为QUOTE则可以得到：由锚链在水平方向的投影长度及钢管、钢桶的倾斜角度可以求得，浮标游动区域大小you=17.521641。钢管与钢桶的受力分析可以求得，四节钢管与钢桶的倾斜角度分别如表所示。表24m/s时钢管及钢桶倾斜角度物体倾斜角度第一节钢管4.321599第二节钢管4.368671第三节钢管4.403959第四节钢管4.439822钢桶4.2347445.2问题二模型问题分析问题二需要调节重物球的质量，使通讯设备在QUOTE时能够正常工作。为了确定重物球的质量，首先将实际风速与临界风速进展比拟，判断此时系统中各物体的状态，与题目中数据进展比拟。在钢桶倾斜角度到达临界角度时，计算锚链与海床的夹角并于题中数据进展比拟，计算重物球的质量。在浮标完全没入海面时，计算相应条件下重物球的质量，从而确定满足条件的重物球的质量围。模型建立Step1.受力分析与数据由于风速QUOTE大于临界速度QUOTE的值，所以当风速为QUOTE时，锚链完全倾斜，锚链倾斜的长度QUOTE。此时，钢管和钢桶的倾斜角度，锚链形状，浮标的吃水深度以及游动区域都可以通过问题一中的模型计算。此题中海水静止，故所有近海水流力为0。Step2.重物球质量求解模型〔1〕钢桶倾斜角度为临界角度时当浮筒倾斜角度QUOTE时，由悬链线模型可以得到：由问题一分析可得:代入数值得到：可解得浮标吃水深度。QUOTE表示重物球质量，此题中未知。由浮标，钢桶，钢管水平受力分析可得：由临界条件中锚以上局部重力与浮力相等，即：则：由悬链线模型可得：由题中条件可知，水深为QUOTE,则：根据浮标的吃水深度，由二分法可求出重物球的质量。〔2〕浮标没入海面时当浮标完全没入海面时，风荷载对浮标无作用力，系泊系统水平方向不受力，则钢管与钢桶，锚链均呈竖直状态，倾斜角度为0，即：此时浮标的吃水深度即为圆柱体的高度：锚链以上局部的竖直高度为：由题中条件可知，水深为QUOTE,则：可以求得，锚链在水中的局部长度为：取锚以上局部进展分析，由系统受力平衡可得，锚以上局部重力与浮力相等，由此可得：则重物球的重力为：代入数据得：5.2.3模型求解Step1.计算QUOTE时物体的状态当风速为QUOTE时，由锚链与海底深度的关系，编程求得浮标吃水深度h=0.703525。浮标的游动区域由锚链在水平方向的投影长度，以及钢管和钢桶在水方向的投影长度能共同决定，设浮标游动区域为QUOTE则可以得到：由锚链在水平方向的投影长度及钢管、钢桶的倾斜角度可以求得，浮标游动区域大小QUOTE。钢管与钢桶的受力分析可以求得，四节钢管与钢桶的倾斜角度分别如表所示。表36m/s时钢管及钢桶倾斜角度物体倾斜角度第一节钢管8.987921第二节钢管9.170037第三节钢管9.240551第四节钢管9.312162钢桶8.901957锚链在锚点与海床的夹角19.422643由题中数据，钢桶的倾斜角度不超过5度，锚链在锚点与海床的夹角不超过16度，计算结果均大于题目要求，故应调整重物球的质量。Step2.计算重物球的质量围当钢桶的倾斜角度为5度时，由模型可求得：QUOTEkg由题中条件可知，水深为QUOTE,则：由锚链与海底深度的关系，编程求得浮标吃水深度QUOTE。锚链在锚点与海床的夹角QUOTE。当浮标完全没入海面时因此，重物球的质量围为QUOTE。显然，由题意可知，随着重物球质量的增大，钢管与钢桶的倾斜角度和锚链在海床的夹角都不断减小。因此，假设锚链在锚点与海床的夹角QUOTE，重物球的质量小于钢桶倾斜角度为5度时重物球的最小质量，此时钢桶倾斜角度大于5度，假设不成立。可以得出，题中所求重物球的质量围是重物球的最大质量围。5.3问题三模型问题分析问题三要求在不同条件下，求出系泊系统中各物体的状态。以型号I锚链为例，当水流方向与风速方向一样时，系统条件最差，分析在不同水深条件下的系泊系统设计。由题中条件确定系统设计的限制条件，对系统各物体进展受力分析，以使整体结果最小，即可得出最优的系泊系统设计。5.3.2模型建立Step1.系统受力分析此题中，系统收到近海风荷载与近海水流力两个外力，且外力随风速大小与水流速大小的变化而变化。系统中各物体受力情况与模型一一样，结合条件，可判断设计系泊系统。Step2.确定目标函数根据题中要求，为了使通讯设备效果最好，浮标吃水深度，游动区域与钢桶倾斜角度应该最小，故本文中确定的优化目标为三个。〔1〕浮标吃水深度〔2〕浮标游动区域〔3〕钢桶倾斜角度Step3.确定限制条件5.3.3模型求解正文六、模型评价〔优点与缺点〕七、参考文献[1]悬链线方程，s://wenku.baidu./view/3343c4e6524de518964b7d9b.html"qq-pf-to=pcqq.c2c，2017年8月16日[2]夏运强,浮筒式防风单点系泊系统研究与应用,2007年6月1日[3][4]附录一三角函数值表角度正弦值〔sin〕正切值(cos)1

0.017452406

0.017455065

2

0.034899497

0.034921

3

0.052335956

0.052407779

4

0.069756474

0.069926812

5

0.087155743

0.087488664

6

0.104528463

0.105104235

7

0.121869343

0.122784561

8

0.139173101

0.140540835

9

0.156434465

0.15838444

10

0.173648178

0.176326981

11

0.190808995

0.194380309

12

0.207911691

0.212556562

13

0.224951054

0.230868191

14

0.241921896

0.249328003

15

0.258819045

0.267949192

附录二临界速度程序p=1025;g=9.8;s=22.05;

L=(1000+40+1300+22.05*7-120/786*p-0.15^2*p*pi)/(pi*p);

symsv;

f*=1.25*(2-L)*v*v;

a=pi*p*g*L;

f0=a-1000*g;f1=a-1010*g;f2=a-1020*g;f3=a-1030*g;f4=a-1040*g;f5=a-2340*g+120/786*p*g+0.15^2*pi*g*p;

%F=[f0,f1,f2,f3,f4,f5];

tan1=f*/(a-1005*g);cos1=sqrt(1/(1+tan1^2));sin1=sqrt(tan1^2/(1+tan1^2));

tan2=f*/(a-1015*g);cos2=sqrt(1/(1+tan2^2));sin2=sqrt(tan1^2/(1+tan1^2));

tan3=f*/(a-1025*g);cos3=sqrt(1/(1+tan3^2));sin3=sqrt(tan1^2/(1+tan1^2));

tan4=f*/(a-1035*g);cos4=sqrt(1/(1+tan4^2));sin4=sqrt(tan1^2/(1+tan1^2));

tan5=f*/(a-1090*g+60/786*p*g+0.5*0.15^2*pi*g*p);cos5=sqrt(1/(1+tan5^2));sin5=sqrt(tan1^2/(1+tan1^2));

fL=L+cos1+cos2+cos3+cos4+cos5;

*=f*/(7*g)*asinh(s*7*g/f*);

y=sqrt(s^2+(f*/(7*g))^2)-f*/(7*g);

ff=y+fL-18;

%fplot(ff,[23,25]);

g=matlabFunction(ff);

[v,k]=eff(g,23,25)附录二风速为QUOTE时p=1025;g=9.8;v=12;

symsL;

f*=1.25*(2-L)*v*v;

a=pi*p*g*L;

f0=a-1000*g;f1=a-1010*g;f2=a-1020*g;f3=a-1030*g;f4=a-1040*g;f5=a-2340*g+120/786*p*g+0.15^2*pi*g*p;

%F=[f0,f1,f2,f3,f4,f5];

tan1=f*/(a-1005*g);cos1=sqrt(1/(1+tan1^2));sin1=sqrt(tan1^2/(1+tan1^2));

tan2=f*/(a-1015*g);cos2=sqrt(1/(1+tan2^2));sin2=sqrt(tan1^2/(1+tan1^2));

tan3=f*/(a-1025*g);cos3=sqrt(1/(1+tan3^2));sin3=sqrt(tan1^2/(1+tan1^2));

tan4=f*/(a-1035*g);cos4=sqrt(1/(1+tan4^2));sin4=sqrt(tan1^2/(1+tan1^2));

tan5=f*/(a-1090*g+60/786*p*g+0.5*0.15^2*pi*g*p);cos5=sqrt(1/(1+tan5^2));sin5=sqrt(tan1^2/(1+tan1^2));

fL=L+cos1+cos2+cos3+cos4+cos5;

s=f5/(7*g);

*=f*/(7*g)*asinh(s*7*g/f*);

y=sqrt(s^2+(f*/(7*g))^2)-f*/(7*g);

ff=y+fL-18;

fplot(ff,[0.6,0.8]);

g=matlabFunction(ff);

%[h,k]=eff(g,0.6,0.64)%结果是0.62，带入求得s为负值，舍去[h,k]=eff(g,0.68,0.72);

*=vpa(subs(*,h));

s=vpa(subs(s,h));

y=vpa(subs(y,h));

fL=vpa(subs(fL,h));

%vpa(subs(F,h))

tan1=(subs(tan1,h));a1=vpa(atan(tan1));

tan2=subs(tan2,h);a2=vpa(asin(tan2));

tan3=subs(tan3,h);a3=vpa(asin(tan3));

tan4=subs(tan4,h);a4=vpa(asin(tan4));

tan5=subs(tan5,h);a5=vpa(asin(tan5));

h

a=vpa([a1,a2,a3,a4,a5]*180/pi)

you=22.05-s+*+vpa(subs(sin1+sin2+sin3+sin4+sin5,h))附录三风速为QUOTE时p=1025;g=9.8;v=36;s=22.05;

symsL;

f*=1.25*(2-L)*v*v;

a=pi*p*g*L;

f0=a-1000*g;f1=a-1010*g;f2=a-1020*g;f3=a-1030*g;f4=a-1040*g;f5=a-2340*g+120/786*p*g+0.15^2*pi*g*p;

tan1=f*/(a-1005*g);cos1=sqrt(1/(1+tan1^2));sin1=sqrt(tan1^2/(1+tan1^2));

tan2=f*/(a-1015*g);cos2=sqrt(1/(1+tan2^2));sin2=sqrt(tan1^2/(1+tan