系泊系统的设计和探究

1、2016年高教社杯全国大学生数学建模竞赛承诺书我们仔细阅读了全国大学生数学建模竞赛章程和全国大学生数学建模竞赛参 赛规则（以下简称为“竞赛章程和参赛规则”，可从全国大学生数学建模竞赛网站下 载）。我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网 上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛章程和参赛规则的，如果引用别人的成果或 资料（包括网上资料），必须按照规定的参考文献的表述方式列出，并在正文引用处予 以标注。在网上交流和下载他人的论文是严重违规违纪行为。我们以中国大学生名誉和诚信郑重承诺，严格遵守竞赛章程和参

2、赛规则，以保证竞 赛的公正、公平性。如有违反竞赛章程和参赛规则的行为，我们将受到严肃处理。我们授权全国大学生数学建模竞赛组委会，可将我们的论文以任何形式进行公开展 示（包括进行网上公示，在书籍、期刊和其他媒体进行正式或非正式发表等）。我们参赛选择的题号（从A/B/C/D中选择一项填写）：我们的报名参赛队号（12位数字全国统一编号）：参赛学校（完整的学校全称，不含院系名）：参赛队员（打印并签名）：1.指导教师或指导教师组负责人（打印并签名）：（指导教师签名意味着对参赛队的行为和论文的真实性负责）日期：年一月一日（请勿改动此页内容和格式。此承诺书打印签名后作为纸质论文的封面，注意电子版论 文中不得

3、出现此页。以上内容请仔细核对，如填写错误，论文可能被取消评奖资格。）2016年高教社杯全国大学生数学建模竞赛送全国评阅统一编号（赛区组委会填写）：全国评阅随机编号（全国组委会填写）：（请勿改动此页内容和格式。此编号专用页仅供赛区和全国评阅使用，参赛队打印后装订到纸质论文的第二页上。注意电子版论文中不得出现此页。）系泊系统的设计和探究摘要 本文利用牛顿力学定律，力矩平衡原理、非线性规划、循环遍历法等方法对系泊系 统进行了设计与探究。通过对系泊系统各组件和浮标运用牛顿经典力学体系进行分析，得 到了各个情况下的钢桶倾斜角度、锚链状态、浮标吃水深度和游动区域。 第一问是求解在风速为12m/s和24m/

4、s时，浮标的吃水深度和游动区域、钢 桶以及钢管的倾斜角度和锚链形态。对于此，首先，我们对浮标、钢管、钢桶、链环进行了基于静力平衡的力学分析，并得到了一系列的方程组；接着，由于钢管、钢桶、链环还满足力矩平衡状态，故得到系泊系统平衡时的刚体力学方程组；然后， 根据系泊系统各组成部件的倾斜角度可以得出其在翌直方向的投影，且它们的投影之和与海水深度存在几何约束。最后，以这个几何约束条件为前提，运用MATLAB中的循环遍历法对该复杂的非线性方程组进行求解，得到在风速为12m/s 时，1至4号钢管的倾斜角度依次为1.1459, 1.1516, 1.1573,1.1688，钢桶的倾斜角 度为1.1860，浮

5、标的吃水深度为0.6870m,浮标的游动区域以锚为圆心，12.5362m为 半径的范围圆。锚链形状为拖地状态，拖地长度6.7m。当风速为24m/s时，从上至下各钢管的倾斜角度为4.408, 4.437, 4.465, 1.1688，钢桶 的倾斜角度为4.566,浮标的吃水深度为0.6974m,游动范围为以锚在海平面度的投影为 圆心，半径为17.8592m的圆。 对于第二问，求解当海面风速为36m/s时，浮标的吃水深度和游动区域、钢桶 以及钢管的倾斜角度和锚链形态。利用第一问中的力学方程和程序，求得钢桶的倾角为 19.5951和四节钢管的倾斜角度依次为19.756、19.755、19.916、2

6、0.076。浮标的游 动区域为以锚在海面上的投影为圆心，半径为18.8828m的圆。由于部分数据与问题二 中钢桶的倾斜角度不超过5，锚链在锚点与海床的夹角不超过16的要求不符，所以通过调节重物球的质量使钢桶的倾斜角度和锚链在锚点与海 床的夹角处在要求的范围之内。借助MATLA程序中的循环遍历法，可以求得重物球的质量3770kg。 对于问题三，在设计系泊系统时，必须要使其能在最恶劣的情况下也能使用，故这 里只讨论风速为36m/s,海水速度为1.5m/s时的系泊系统各构件的状态。假设在问题二 的情况下，重物球质量为3770kg，海水深度分别为16m和20m的情况下，浮标的吃水深度和游动区域、钢桶以

7、及钢管的倾斜角度和锚链形态。以钢桶倾斜角度和浮标游动范围为评价指标，运用TOPSIS的方法得出在海水深度为16m，在风速为36m/s，水流速度为1.5m/s时，五种锚链的钢桶倾斜角都为 3.8291 ,比较浮标浮动范围得到，V号锚链浮动半径最小为18.9175m 。在海水深度为20m,其他条件不变时，得出钢桶倾斜角度和浮标游动范围数据不易直接比较，故在此 采用TOPSIS方法进行比较，得到五种锚链与最有方案的接近程度大小关系：IV III II I V，因此可以得出在海水深度为20m时IV号锚链最为合适。关键词：牛顿第二定律力矩平衡原理非线性规划循环遍历法、问题重述近浅海观测网的传输节点由浮标

8、系统、系泊系统和水声通讯系统组成（如图一所示）。某 型传输节点的浮标系统可简化为低面直径2m高2m的圆柱体，浮标的质量为1000kg。系泊系统由钢管、钢桶、重物球、电焊锚链和特制的抗拖移锚组 成。锚的质量为600kg，锚链选用无档普通链环，常用型号及其参数已在附表中列出。 钢管共4节，每节长度1m直径为50mm每节钢管的质量为10kg水声通讯系统安装在 一个长1m外径30cm的密封圆柱形钢桶内，设备和钢桶总质量为100kg。钢桶上接第4 节钢管，下接电焊锚链第1节。现由题可知，如果要使该系统正常工作，那么就要使锚 链末端和锚链接处的切线方向与海床的夹角不超过16度，否则锚会被拖行，致使节点移

9、位；若钢桶的倾斜角度（钢桶与翌直线的夹角）超过5度，则设备的工作效果较差。因 此为了控制钢桶的倾斜角度，可在钢桶与电焊锚链连接处悬挂重物球。图1传输节点示意图（仅为结构模块示意图，未考虑尺寸比例）系泊系统的设计问题就是确定锚链的型号、长度和重物球的质量，使得浮标的吃水 深度和游动区域及钢桶的倾斜角度尽可能小。问题1某型传输节点选用II型电焊锚链22.05m,选用的重物球的质量为1200kg。现将该节点布放在水深为18m海水密度为1.025 1 03kg/m3的海床平坦的海 域。请计算当海面风速分别为12m/s和24m/s且海水静止时，钢桶和各节钢管的倾斜角 度、锚链形状、浮标的吃水深度和游动区

10、域。问题2请在问题1的假设下，计算海面风速为36m/s时钢桶和各节钢管的倾斜角 度、锚链形状和浮标的游动区域。并试调节重物球的质量，使得钢桶的倾斜角度不超过5 度，锚链在锚点与海床的夹角不超过16度。问题3而由于潮汐等因素的影响，布放海域的水深实际介于16m20r之间。布放 点的海水速度最大可达到1.5m/s、风速最大可达到36m/s。请给出考虑风力、水流力和 水深情况下的系泊系统设计，分析不同情况下钢桶、钢管的倾斜角度、锚链形状、浮标的 吃水深度和游动区域。二、问题分析本题是以系泊系统设计为背景的力学分析问题。2.1问题一分析首先，建立以锚为原点，其余各部件都在第一象限的平面直角坐标系。然后

11、，根据牛顿经典力学理论，对浮标、钢管和钢桶分别进行受力分析，可得出三 个方程组。接着，由于钢桶与钢管之间存在着力矩平衡，因此可以得出钢桶、钢管之间的关系方程；对于锚链，可对其每一个小段进行受力分析，得出其有关方程 式；因为系泊系统各组成部分在翌直方向上的投影之和为海平面与海床之间的距离，即 18米，为有效限定条件，所以可列出浮标、钢管、钢桶、锚链在翌直方向上投影之和与 18米之间的关系方程式。最后，使用迭代算法在MATLABA件上对所有方程进行求解；分别求出风速在 12m/s和24m/s时时钢桶和各节钢管的倾斜角度、锚链形状、浮标的吃水深度和游动区 域。2.2问题二分析根据问题一中的求解思路，

12、将风速值36m/s代入，求解出数值，再与题目中的条件 钢桶的倾斜角小于等于5度，锚在锚点与海床的夹角小于等于16度进行对比，判断其 是否合理。如不合理，运用循环遍历法，逐渐增加重物球的质量，直至得到一个合理的质 量。2.3问题三分析由于潮汐等因素的影响，布放海域的实测水深介于16m20n之间。布放点的海水 速度最大可达到1.5m/s、风速最大可达到36m/s。请给出考虑风力、水流力和水深情况 下的系泊系统设计，分析不同情况下钢桶、钢管的倾斜角度、锚链形状、浮标的吃水深度 和游动区域。三、模型假设1、海平面与海底平面水平。2、浮标始终海平面保持垂直状态。3、在一定时间内，风速和风向恒定。4、风向

13、为水平风向。5、该地区重力加速度为9.8m/s2。&假设重力球、锚链、钢管为同一种材料，均为密度为7.9 103kg/m3的钢。7、锚所受浮力忽略不计。四、符号说明符号单位符号意义F风N浮标所受风力F浮标N浮标所受浮力G厂标N浮标的重力hm浮标吃水深度F管N钢管所受浮力G E管N钢管的重力管1N钢管所受拉力(i=1, 2, 3, 4, 5)度钢管拉力与翌直线夹角(i=1, 2, 3, 4, 5)片度钢管倾斜角度(j=1,2,3,4)Lgm钢管长度FtN钢桶所受浮力GtN钢桶的重力TtN钢桶所受拉力e度钢桶倾斜角度Ltm钢桶长度F链N链环所受浮力=连N链环的重力T链iN链环所受拉力(i=1,2,

14、3.211 )半i度链环拉力与翌直线夹角V i度链环倾斜角度(i=1,2,3.210)G锚N锚的重力N锚所受拉力FfN锚所受摩擦力FnN锚的支持力GqN重物球的重力Rm浮标的游动半径五、模型的建立与求解5.1问题一的模型建立与求解5.1.1模型一的建立已知重力球的质量1200kg,海水深度18米，海水密度1.025 xi03kg/mao求当海面 风速分别为12m/s和24m/s且海水静止时，钢桶和各节钢管的倾斜角度、锚链形状、浮 标的吃水深度和游动区域。我们可以在系统处于风力恒定的平衡状态下，对系泊系统进行 局部受力分析。图2对浮标进行受力分析浮标受力分析图由图2可知，浮标在水平方向上受到翌直

15、向上的浮力F浮，翌直向下的重力G标，水平向右的风力F风和钢管产生的斜向下的拉力T其中：F浮g(d/ 2)2F 风=0.625 ( H-h)根据牛顿第一正律列出卜列方程组F 浮=G 标 Ti cos -F 风= si化简得:(5)G管图3钢管受力分析图i=atan F 风 F 浮-G 标 HTi= .( )对钢管进行受力分析F管为钢管所受浮力，G管为钢管自身重力，T管i表示第i-1个钢管对第i个钢管的拉 力，冷表示第i-1个钢管对第i个钢管的拉力与其浮力所成的夹角。其中，1表示第i个钢 管，i= 1, 2, 3, 4, 5。图3表示的为第一个钢管的受力分析图，由题知，四个钢管材质大小都相同，因

16、此，所受浮力和自身重力也相同。其中，当i=1时，T管i表示浮标对第一个钢管的拉力， T管2表示第一个钢管对第二个钢管的拉力，表示浮标对第一个钢管的拉力与其浮力所成的夹角，:2表示第一个钢管对第二个钢管的拉力与其浮力所成的力，浮标对其有一个斜向上的拉力，第二根钢管对其有一个斜向下的拉力。根据受力分析平衡列出力学方程组。翌直方向：F 管 +T 管 iSi n_： 二 G 管 T 管 i+1 CS ：i+1水平方向：丁管C0S: i =丁管问 Sin , i+1化简得：Si+1 二 arctanT 管i Sin :打：T管i COS：、 F 管G 管夹角。由图可知，第一根钢管在翌直方向上受到翌直向上

17、的浮力和翌直向下的重对钢管倾斜角度进行分析因为钢管存在固定的转轴，所以由力矩平衡可画出下列图4的受力分析(10)T 管j+1 二 T 管 iSin : Sin+1 图。其中表示第j个钢管的倾斜角度，j= 1, 2, 3, 4。和图4钢管倾斜角度示意图图4为钢管倾斜角度示意图，（j=1,2,3,4 ）表示第j个钢管的倾斜角度，根据力矩平衡原理列出力学方程组。(11)1 L T管i d sin (:）L（G管-F管）siL T管i sin ( : jT.-j)21L 2 管F 管)sin 1 j(12化简得:(13):etan并孔心卫2 T| cos :（4）钢桶进行受力分析Ft为钢桶所受浮力，T

18、管5为第四根钢管对钢桶的拉力，T链1为第一根锚链对钢桶的拉力，Gt为钢桶自身重力，Gq为重物球的重力。：_5为T管5与Ft所成夹 角，5为T链1与Gt所成夹角。图5钢桶受力分析图(14)通过图5对钢桶进行受力分析得:T 管 5 C0S：5 Ft = 丁链1 C0S 冷 Gt TOC o 1-5 h z T 管 5 部二 5 叮链 i sin (15)化简得：1 = arctanMT 管 5 sin 二 5： T 管 5 cos 二 5 Ft-Gt-Gq 门(16)T 链 1=Tt5Sin： sin 1(17)（5）对钢桶倾斜角度进行分析因为钢桶存也在固定的转轴，所以由力矩平衡可画出下列图6的受

19、力分析图。其中二表示钢桶的倾斜角度。图6钢桶力矩平衡分析示意图由于钢桶还符合刚体的力矩平衡，可以列出以下方程：(18)1LtT 管 5驴（花：5）=Lt（GFt）s(19)(20)v - arctan.210其中，1表示第i个锚链管，i= 1 , 2, 3T G连图6链环受力分析图图6为链环受力分析图，通过对链环进行受力分析得:（21）化简得T链isin：i二 T链.Sin 帚i+i= arctan 连jn “T链fos (22)(23)(24)T 链 icos口 链=丁连 i+icos -i+iG 链T链i+i二T链isin；isin 帚（7）对连环倾斜角度进行分析个钢管的倾斜角度，i= 1

20、,2, 3,-.210.LTi链isin(:-11 -2L （G链-F链）COS(25)T链i 1sin(i - 1+1) =1 Li (G链i- F链i) COSi(26)由于连环符合刚体的力矩平衡,可以列出以下方程,其中i表示第图7链环力矩平衡分析化简得：=arcta n 2 T I5.1.2模型一的求解1风速为12m/s（1）计算各节钢管和钢管的倾斜角度。用牛顿经典力学理论对整个系统进行力学分析，可以得到海底平面对锚的摩擦力 Ff，浮标所受浮力F浮（F浮=gV,V二二R2h）和浮标吃水深度h的值（浮标所受浮力F浮和浮标吃水深度h的值上面已求出）。将锚链、锚和重物球看做一 个整体用牛顿经典

21、力学理论对其进行力学分析，可得钢桶所受锚链拉力，根据力矩平衡可 得钢桶与水平线夹角：5 ；同理，根据牛顿经典力学理论，依次对钢管1,2,3,4进行力 学分析，可求出：j， : 2，3，: 4。1）以H为几何约束条件，运用MATLABA的循环遍历法对以钢管的受力分析的方程组（9）（ 10）以及力矩平衡方程（13）进行求解，得到在风速为12m/s时，钢管的倾斜角度见表2。表2第一个钢管的倾斜角度1.1459第二个钢管的倾斜角度1.1516第二个钢管的倾斜角度1.1573第四个钢管的倾斜角度1.1688（2）计算钢桶的倾斜角度运用MATLABA的循环遍历法对以钢管的受力分析的方程组（9）（ 10）进

22、行 求解得：T 管 5 =11607N=84.3:5将数据代入力矩平衡方程（13）进行化简得钢桶的倾斜角度为：-arcta n 2 T 一、讷：5_（G5F5）=1.18602，凡 cos%（3）求解水深度h由于系泊系统各组成部分在翌直方向上的投影之和等于面与海床之间的距离，即18 米，以此为限定条件，列出浮标、钢管、钢桶、锚链在翌直方向上投影之和与18米之间 的关系方程式，再根据浮标体积进而求出浮标的吃水深度。设H为所有部件在翌直方向上的投影之和，即H=h+LG + Lt Lh其中，h为浮标吃水深度，lg为钢管在翌直方向的投影之和，lt钢桶翌直方向的投影，lh为锚 链在翌直方向的投影之和。4

23、210因为 LG=V LgCOS： ，Lti =04H=h+、Lg cos: “-1i -0二 Lt -5，Lh=vLh COs ；所以21057 LhCOs i -0由题可知，H III II I V，很容易得出在海水深度为20m,风速为36m/s，水流速度为1.5m/s时，IV号锚链最为合适。综合1、2两种情况可知当海水深度较浅时，选用V号锚链，系泊系统工作情况最佳，当海水深度较深时选用IV号锚链系泊系统工作情况最佳。六模型评价模型优点：使用MATLA清晰地表现出来锚链的形态，浮标的吃水深度，以及系 统各部件的倾斜角度，本论文运用图像和表格使整体模型更加清晰 明了。2）分析了在不同海水深度

24、状态下，使用不同锚链时，浮标吃水深度、 钢桶倾斜角度以及锚链在锚点与海床的夹角的变化情况。3）选用锚链时运用较为客观的评选方案，尽可能的避免了主观因素带 来的影响。模型缺点：1）假设的风速方向较为固定，在实际情形中，风速方向不是固定不变 的。2）未能给出系泊系统随各个变量的具体关系。clear;clc重物球的质量 链环的个数mq=2200; n=210;%mi n=i nf;%i nffor h=0:0.001:2%thital=zeros(1,4);%Ft=zeros(1,5);%a=zeros(1,5);%Ft2=zeros(1, n+1);%gama=zeros(1, n+1);%thi

25、ta2=zeros(1, n)+pi/2;%s=4-2*h;%v=12;%m=1000;%p=1025;%g=9.8;%v0=pi*h;%Ffen g=0.625*s*vA2;%Ffu=p*g*vO;Gfu=m*g;if Ffu-Gfu0为无穷大浮标的吃水深度范围钢管与翌直平面的夹角浮标与钢管，钢管与钢桶之间的力浮标，钢管，钢桶之间的力与翌直平面的夹角链环所受的力链环之间的力与翌直平面的角度链环与翌直平面的夹浮标受风面积风速浮标质量海水的密度标准重力加速度浮标所排海水体积风力大小enda(1)=ata n(Ffe ng/(Ffu-Gfu);Ft(1)=sqrt(Ffe ngA2+(Ffu-Gf

26、u)A2); % 钢管Vgua n=pi*0.025A2;Ggua n=10*g;Fgua nfu=p*g*Vgua n; for i=1:4a(i+1)=ata n( (Ft(i)*si n( a(i)/(Ft(i)*cos(a(i)+Fgua nfu-Ggua n); %钢管所受下一节钢管拉力与翌直平面的夹角Ft(i+1)=Ft(i)*s in (a(i)/si n( a(i+1);%钢管所受下一节钢管的拉力thital(i)=ata n(Ft(i)*si n( a(i)*1/(Fgua nfu-Ggua n)*0.5+Ft(i)*cos(a(i);%钢管与翌直平面的夹角end%钢桶的受力

27、分析Vt=0.15A2*pi;%钢桶的体积Vq=mq/7900;%重物球的体积Gt=100*g;Gq=mq*g;Ftfu=p*g*Vt;Fqfu=p*g*Vq;gama(1)=ata n(Ft (5) *si n(a (5) )/(Ftfu+Ft (5) *cos(a(5)-Gt-Gq+Fqfu);%钢桶所受锚链拉力与翌直平面的夹角Ft2(1)=Ft (5) *si n(a (5) )/si n( gama (1)；钢桶所受锚链拉力beta=ata n(Ft( 5)*si n(a(5)*1/(Ftfu-Gt)*0.5+Ft(5)*cos(a(5)*1); %钢桶与翌直平面的夹角mm=0.735

28、; roum=7900;度 7900kg/m*3Vm=mm/roum;Fmfu=p*g*Vm;Gm=mm*g;Lm=0.105;for1: n%锚链每节链环的质量 锚链所用钢的密每节链环的体积每节链环所受的浮力每节链环所受的重力每节链环的长度gama(i+1)=ata n(Ft2(i)*si n( gama(i)/(Ft2(i)*cos(gama(i)+Fmfu-链环所Gm); %受下一节链环拉力与翌直平面的夹角if gama(i+1)0gama(i+1)=gama(i+1)+pi;endFt2(i+1)=Ft2(i)*si n( gama(i)/si n( gama(i+1);*才而毛可链环

29、所受到的下一节链环对其的拉力thita2(i)=ata n(Ft2(i)*si n( gama(i)*Lm/(Fmfu-Gm)*Lm/2+Ft2(i)*cos(gama(i)*Lm);%链环与翌直线的夹角if thita2(i)0thita2(i)=thita2(i)+pi;end总高度浮标的吃水深度 逐次逼近end H=h+sum(cos(thital)+Lm*sum(cos(thita2)+cos(beta); % if abs(H-18)mi n min h=h;%min=abs(H-18);%mi nH=H; mi nthial=thital; min thita2=thita2; m

30、in beta=beta; mi nF t2=Ft2; end end t=1:210; xx1(t)=mi nthita2; xx2=xx1(211-t); y=cumsum(Lm*cos(xx2); x=cumsum(Lm*si n( xx2);for t=1:1:210;if y(t)0；x(t)=x(t);y(t)=0;endendplot(x,y)youdo ng=sum(si n(min thial)+Lm*sum(si n( mi nthita2)+si n(mi nbeta)问题二的程序:clc,clear n=210; lm%=;重物球的质量lh司;lbeta司;mi n=%

31、f;链环的个数for mq=1200:10:4000%inf为无穷大min h=0; minH=0;min teba=0;mi nthital=zeros(1,4);mi nthita2=zeros(1, n)+pi/2; minF t2=zeros(1, n+1); for h=0:0.001:2;thital=zeros(1,4);Ft=zeros(1,5);a=zeros(1,5);夹角Ft2=zeros(1, n+1);gama=zeros(1, n+1);thita2=zeros(1, n)+pi/2;beta=0; s=4-2*h;v=36;m=1000;p=1025; g=9.8

32、;v0=pi*h;Ffen g=0.625*s*VA2;%Ffu=p*g*v0; Gfu=m*g; ifFfu-Gfu0%浮标的吃水深度范围%钢管与翌直平面的夹角浮标与钢 管，钢管与钢桶之间的力浮标，钢管，钢桶 之间的力与翌直平面的链环之间的力链环之间的力与翌直平面的角度链环与翌直平面的夹角钢桶与翌直平面的夹角浮标受风面积浮标质量海水的密度标准重力加速度浮标所排海水体积风力大小enda(1)=ata n(Ffe ng/(Ffu-Gfu);Ft(1)=sqrt(Ffe ngA2+(Ffu-GfuF2);% 钢管Vgua n=pi*0.025A2;Ggua n=10*g;Fgua nfu=p*g*

33、Vgua n;for i=1:4a(i+1)=ata n( (Ft(i)*si n(a(i)/(Ft(i)*cos(a(i)+Fgua nfu-Ggua n);%钢管所受下一节钢管拉力与翌直平面的夹Ft(i+1)=Ft(i)*si n(a(i)/si n( a(i+1); %钢管所受下一节钢管的拉力钢管与翌直平面的夹角thital(i)=ata n(Ft(i)*si n( a(i)*1/(Fgua nfu-Ggua n)*0.5+Ft(i)*cos(a (i);end钢桶的体积重物球的体积Vt=0.15A2*pi;Vq=mq/7900;Gt=100*g;Gq=mq*g;Ftfu=p*g*Vt;

34、Fqfu=p*g*Vq;gama(1)=ata n(Ft (5)*si n(a(5)/(Ftfu+Ft(5)*cos(a(5)-Gt-Gq+Fqfu);%钢桶所受锚链拉力与翌直平面的夹角Ft2(1)=Ft (5)*si n(a (5) )/si n( gama)；%钢桶所受锚链拉力beta=ata n(Ft( 5)*si n(a(5)*1/(Ftfu-Gt)*0.5+Ft(5)*cos(a(5)*1);%钢桶与翌直平面的夹角每节链环的质量锚链所用钢的密度7900kg/m*3每节链环的体积每节链环所受的浮力每节链环所受的重力每节链环的长度mm=0.735;roum=7900;Vm=mm/roum

35、;Fmfu=p*g*Vm;Gm=mm*g;Lm=0.105;for i=1: ngama(i+1)=ata n(Ft2(i)*si n( gama(i)/(Ft2(i)*cos(gama(i)+Fmfu-Gm);%链环所受下一节链环拉力与翌直平面的夹角 if gama(i+1)0 gama(i+1)=gama(i+1)+pi;endFt2(i+1)=Ft2(i)*si n( gama(i)/si n(gama(i+1); %链环所受到的下一节链环对其的拉力thita2(i)=ata n(Ft2(i)*si n( gama(i)*Lm/(Fmfu-Gm)*Lm/2+Ft2(i)*cos(gama

36、(i)*Lm);%链环与翌直线的夹角if thita2(i)0thita2(i)=thita2(i)+pi;endendH=h+sum(cos(thital)+Lm*sum(cos(thita2)+cos(beta); %总高度if abs(H-18)5con ti nue;endif 90-mi nthita2( n)*180/pi16 con ti nue;endif abs(mi nH-18)0.2;con ti nueendlmq=lmq mq;lh=lh min h; lbeta=lbeta min beta;end问题二clear;clcmq=3770;%重物球的质量链环的个数n=

37、283;%为无穷大mi n=i nf;%inf浮标的吃水深度范围for h=0:0.001:2%钢管与翌直平面的夹角thital=zeros(1,4);%浮标与钢管，钢管与钢桶之间的力浮Ft=zeros(1,5);%标，钢管，钢桶之间的力与翌直平面a=zeros(1,5);%的夹角链环所受的力Ft2=zeros(1, n+1);%链环之间的力与翌直平面的角度链环与gama=zeros(1, n+1);%翌直平面的夹浮标受风面积thita2=zeros(1, n)+pi/2;%风速s=4-2*h;%浮标质量海水的密度标准重力加速度v=36;%浮标所排海水体积m=1000;%p=1025;%g=9

38、.8;%v0=pi*h;%Ffen g=0.625*s*VA2;%风力大小Ffu=p*g*vO;Gfu=m*g;sw 仁 2*h;vw=1.5;Fw 仁 374*sw1*vwA2if Ffu-Gfu0enda(1)=ata n( (Ffe ng+Fw1)/(Ffu-Gfu);Ft(1)=sqrt(Ffe ng+Fw1)A2+(Ffu-GfuF2); % 钢管Vgua n=pi*0.025A2;Ggua n=10*g;Fgua nfu=p*g*Vgua n;sw2=0.05*1;Fw2=374*sw2*vwA2for i=1:4a(i+1)=ata n(Ft(i)*si n(a(i)+Fw2)/

39、(Ft(i)*cos(a(i)+Fgua nfu-Ggua n);%钢管所受下一节钢管拉力与翌直平面的夹角Ft(i+1)=Ft(i)*si n(a(i)/si n(a(i+1);%钢管所受下一节钢管的拉力thital(i)=ata n( Ft(i)*si n(a(i)*1/(Fgua nfu-Ggua n)*0.5+Ft(i)*cos( a(i); % 钢 管与翌直平面的夹角end嗨冈桶的受力分析Vt=0.15A2*pi;%钢桶的体积Vq=mq/7900;%重物球的体积Gt=100*g;Gq=mq*g;Ftfu=p*g*Vt;Fqfu=p*g*Vq;sw3=0.3*1;fw3=374*sw3*vwA2gama(1)=ata n(Ft (5) *si n(a (5) )+fw3)/(Ftfu+Ft(5)*cos(a(5)-Gt-Gq+Fqfu); %钢桶所受锚链拉力与翌直平面的夹角Ft2(1)=(Ft (5) *si n(a( 5)+fw3)/si n(gama(1);% 钢桶所受锚链拉力beta=ata n( Ft(5)*si n( a(5)*1/(Ftfu-Gt)*0.5+Ft(5)*cos(a(