2016年全国大学生数学建模大赛国家一等奖优秀论文,系泊系统的设计

整理文档很辛苦,赏杯茶钱您下走!

免费阅读已结束,点击下载阅读编辑剩下 ...

阅读已结束,您可以下载文档离线阅读编辑

资源描述

1系泊系统的设计摘要本文详细对系泊系统的各个机构进行了力学分析,针对系泊系统的要求,建立优化模型,求解系泊系统在多种环境下的最优解,使得浮标游动范围,吃水程度和钢桶倾斜角度尽可能的小。针对问题一,本文对系泊系统的受力及力矩进行了分析,基于浮标倾斜的考虑,得到了平衡状态下关于受力平衡及力矩平衡的方程组。由于方程组数量较多及相互影响的特点,直接求解十分困难。因此我们考虑以浮标两边的浸水长度12,hh为变量,利用搜索算法对方程组进行求解,并得到相应的结果。如当风速为12m/s时,钢桶的倾斜角度1.0405°,从上到下钢管的倾斜角度分别为1.0086°、1.0146°、1.0206°、1.0267°,浮标吃水深度0.735m,浮标游动区域半径14.4429m。针对问题二,首先将风速为36m/s的情况代入问题一建立的模型中,但是得到的结果不满足题目所给定的要求。则考虑在重物球质量一定的条件下,以浮标的吃水深度和游动区域及钢桶的倾斜角为目标,建立了一个单决策变量的多目标最优系泊模型,相比于问题一,此问的变量更多,更加难于求解,故考虑将多目标转化成单目标的问题进行求解,并继续使用搜索法对问题进行求解。最后找到了三组可行解,其中最优解是重力球的质量为2102kg.针对问题三,本文中有三个决策变量以及三个变系数,相比于前两问,无论是计算量还是计算维数,难度更大。为了求解该问,建立了一个多决策变量的多目标变系数的最优系泊系统模型,为了简便运算,我们建立了变步长的搜索算法,并最终求解得到结果,得到的一组解为:选用了III型号的锚链,重物球质量为2800kg,锚链长度为23.4m。针对论文的实际情况,对论文的优缺点做了评价,文章最后还给出了其他的改进方向,以用于指导实际应用。关键词:系泊系统设计;力的平移定理;多目标;优化模型;搜索算法21.问题的重述一个由浮标系统、系泊系统和水声通讯系统组成的近浅海观测网的传输节点。可以简化看作是一个浮于海平面的圆柱体浮标通过四节钢管链接装有通讯设备的钢桶,钢桶再通过锚链链接一个可以移动的锚,锚沉在海床上。为了不让锚被拖行要求锚链末端与锚的链接处的切线方向与海床的夹角不超过16度,并且要保证通讯设备的工作效果,钢桶与电焊锚链链接出可悬挂重物体,使得钢桶的倾斜角度(钢桶与竖直线的夹角)不超过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。给出考虑风力、水流力和水深情况下的系泊系统设计,分析不同情况下钢桶、钢管的倾斜角度、锚链形状、浮标的吃水深度和游动区域。2.模型的假设(1)假设两钢管用铰链链接在一起,可以自由转动。(2)假设风是沿平行海平面的方向吹来。(3)假设忽略锚链和重物球的浮力影响。(4)假设忽略由于其它原因(如温度,湿度,受力等)而对各个系统产生的形变。(5)假设浮标不会出现因为风力过大而被吹翻的情况。(6)假设各个系统之间连接处的长度忽略不计。3.符号说明windF:浮标所受水平方向的风力大小;buoyantAF:浮标所受浮力大小;buoyantBF:每节钢管所受浮力大小;buoyantCF:钢桶所受浮力大小;0BF:第一节钢管对浮标拉力的大小(或浮标对第一节钢管拉力的大小);iBF:第i节钢管下端所受拉力的大小(1,2,3,4i);0()CDFF:钢桶所受锚链的拉力大小(或锚链上方第一段所受钢桶的拉力大小);3jDF:锚链第j段(自上而下)下端所受拉力大小(1,2,,jN,其中N为张紧状态锚链的段数);fF:锚与海底间的摩擦力;maxF:锚与海底间的最大静摩擦力;waterAF:浮标所受的水流力大小;iwaterBF:第i节钢管所受水流力大小(1,2,3,4i);waterCF:钢桶所受水流力大小;:锚与海底间的最大静摩擦系数;NF:海底对锚的支持力;AG:浮标的重力;BG:每节钢管的重力;CG:表示钢桶的重力;DG:锚链每段链环的重力;EG:锚的重力;0B:第一节钢管对浮标拉力与水平方向的夹角;iB:第i节钢管下端所受拉力与水平方向的夹角(1,2,3,4i);0()CD:钢桶所受锚链的拉力与水平方向的夹角;jD:锚链第j段(自上而下)下端所受拉力与水平方向的夹角(1,2,,jN);A:浮标的水平倾斜角度;iB:第i节钢管的水平倾斜角度(1,2,3,4i);c:钢桶的水平倾斜角度;12,hh:圆柱体浮标两端的浸水长度,其中12hh;h:浮标的吃水深度;4windv:水平方向的风速;windS:浮标在风向法平面的投影面积;overflowAV:浮标的排水体积;waterV:水流速度大小;Am:浮标的质量;Bm:每节钢管的质量;Cm:钢桶系统的质量;m:重物球的质量;:海水密度;D:锚链的线密度;DNC:法向阻力系数;DNC:切向阻力系数;Ar:浮标的底面半径;Br:每节钢管的底面半径;BL:每节钢管的长度;Cr:钢桶的底面半径;CL:钢桶的长度;DL:锚链每段链环的长度;H:海水深度;s:锚链某点至锚链上端顶点的弧长;AS:浮标浸入水下部分在水流力法平面的投影面积;iBS:第i节钢管在水流力法平面的投影面积(1,2,3,4i);CS:钢桶在水流力法平面的投影面积;5T:锚链某点处所受拉力的大小;:锚链某点处所受拉力与水平方向的夹角;AX:平衡状态下浮标的水平距离;BX:平衡状态下钢管的水平距离;CX:平衡状态下钢桶的水平距离;DX:平衡状态下锚链的水平距离;X:平衡状态下整体的水平距离;Dy:锚链竖直方向上的高度;max:锚不滑动时,锚链在锚点与海床的最大夹角;D:单位长度上锚链所受的法向水流力大小;D:单位长度上锚链所受的切向水流力大小;4.模型的建立与求解4.1问题的分析4.1.1系统平衡的两种状态——张紧状态和松弛状态及游动区域的理解本文考虑的在风力及水流力作用下的系泊系统的设计问题,首先我们考虑一种简单情形,即海水静止且风向水平情况,此时整个系泊系统可能会因风力移动达到一个平衡的状态,由于风向恒定我们可以认为最终的平衡状态中浮标、钢管、钢桶、锚链和锚大致位于同一个平面内,如图1所示。6图1传输节点示意图在平衡状态下,各系统所受力及力矩分别达到平衡,此时浮标会与锚形成一定的水平距离。此时的平衡状态可能分为两种:一种是张紧状态,即锚链的每一段都受到了前后端的拉力从而使整个系统平衡;另一种是松弛状态,即锚链上方的一部分存在拉力,而下方部分散落在水底呈无拉力状态达到平衡。在张紧状态下达到平衡时,浮标到达某个定点,此时该点与锚得水平距离为半径以锚和浮标的水平距离为半径、锚所在位置为圆心做一个圆面,该区域即浮标的游动区域,如图2所示。但在松弛状态下由于松弛部分锚链的可伸缩性,浮标的位置会处于一定的区域,当松弛部分锚链长度最大时浮标位于最远位置,以此时浮标与锚的水平距离为半径、锚所在位置为圆心的圆形区域即其游动区域。图2浮标游动区域俯视示意图4.1.2浮标的倾斜问题及对吃水深度的理解在对组成各个系统的零件受力分析的过程中,会出现无法将物体看作是均匀的情况,这样一来,物体的各力的受力点将不再是质点,为此应尽可能将受力点移动到质心上。这时,根据力的平移定理,在将各力平移到质心上去的过程中会产生力矩,且对于一个平衡的物体其力矩和应为零。在力矩平衡的作用下,我们发现平衡状态下浮标会产生一定的倾斜角度。假设在浮标的底面与海平面平行,即没有倾角时,其受力分析如图3所示。7图3浮标垂直海平面的受力分析图由图可知除浮力buoyantAF外,风力windF和钢管拉力0BF都会产生非零力矩,且方向相同,导致力矩和不为零,从而会产生旋转,不会达到平衡状态,我们认为在平衡状态下,浮标底面相对海平面存在一个倾斜角度A。吃水深度为浮标底部到海面的距离。由于浮标倾斜,导致浮标各点侵入海水的长度不同,如图4所示。假设浮标两边侵入水中的长度分别为12,hh,其中12hh,则吃水深度即两边浸水长度的平均值122hh。图4倾斜浮标的浸水示意图4.2问题1模型的建立根据对问题的分析发现,问题旨在讨论在不同情形达到平衡状态下,分析各钢管的倾斜角度、锚链形状、浮标的方程,通过已知的信息,求得系统各部分所受拉力大小及其角度情况和各部分本身的倾斜角度。4.2.1浮标的力平衡方程通过浮标的受力分析可得,浮标受到竖直向上的浮力buoyantAF,水平向右的风力windF,重力AG以及钢管对其的拉力0BF,其受力分析图如图5所示:8图5浮标的受力分析图根据平面力系平衡原理,可得浮标的力平衡方程如下:在水平方向上,有00cos,BBwindFF(1)其中,0B为钢管对浮标拉力与水平方向的夹角,windF为水平风力的大小,有如下的近似公式(S为物体在风法平面的投影面积)20.625();windwindFSv(2)在竖直方向上,有00sin,buoyantABBAFmgF(3)其中,有,buoyantAoverflowAFgV(4)而2221211211()(),22overflowAAAAVrhrhhrhh(5)在这儿,Am为浮标质量,12,hh分别表示浮标两边的浸水长度(见图4),Ar为浮标的底面半径,overflowAV为浮标的排水体积,为海水密度,g为重力加速度。4.2.2浮标在风向法平面的投影面积公式在计算风力windF时,需用到浮标在风向法平面的投影面积S,由于浮标倾斜,该投影区域并非一标准的矩形,而是由两个部分构成,其中的下方部分为矩形,而上方部分为倾斜圆在水平面上投影的一半,即半椭圆的一半,最终得到浮标在风向法平面的投影面积如图6所示。图6浮标在风向法平面的投影面积示意图9通过平面几何关系可得,第一部分矩形的面积1S满足:11222122=2[0.5()],()(2)AAAArSrLhhhhr(6)第二部分半椭圆的面积2S满足:212(),4AhhrS(7)其中,Ar为浮标的底面半径,12,hh分别表示浮标两边的浸水长度,AL表示浮标的长度。4.2.3钢管的平衡方程通过分析可得,每节钢管的受力分析图如图7所示。图7第i节钢管的受力分析图根据平面力系平衡原理及力的平移定理(力矩平衡公式),有如下方程:水平方向力平衡方程11coscos(1,2,3,4),iiiiBBBBFFi(8)竖直方向力平衡方程11sinsin(1,2,3,4),iiiiBBBBBbuoyantBFmgFFi(9)力矩平衡方程1111sin()sin()(1,2,3,4).22iiiiiiBBBBBBBBFLFLi(10)在这儿,buoyantBF为每节钢管所受浮力大小,满足2,buoyantBBBFgrL(11)而iBF为第i节钢管下端所受拉力大小,Bm为每节钢管的质量,g为重力加速度,iB是第i节钢管与水平方向的夹角(钢管的倾斜角度),iB是拉力iBF与水平方向的夹角,10Br为钢管的底面半径,BL为钢管的长度。4.2.4钢桶系统的平衡方程钢桶的受力分析如图8所示。图8钢桶的受力分析图根据平面力系平衡原理及力的平移定理(力矩平衡公式),有如下方程:水平方向力平衡方程44coscos,CCBBFF(12)竖直方向力平衡方程44sin()sin,CCCBBbuoyantCFmmgFF(13)力矩平衡方程44111sin()sin()cos

1 / 38
下载文档,编辑使用

©2015-2020 m.777doc.com 三七文档.

备案号:鲁ICP备2024069028号-1 客服联系 QQ:2149211541

×
保存成功