全国大学生数学建模大赛国家一等奖优秀论文-系泊系统的设计.docx

系泊系统的设计摘 要本文详细对系泊系统的各个机构进行了力学分析,针对系泊系统的要求,建立优化模型,求解系泊系统在多种环境下的最优解,使得浮标游动范围,吃水程度和钢桶倾斜角度尽可能的小.针对问题一,本文对系泊系统的受力与力矩进行了分析,基于浮标倾斜的考虑,得到了平衡状态下关于受力平衡与力矩平衡的方程组.由于方程组数量较多与相互影响的特点,直接求解十分困难.因此我们考虑以浮标两边的浸水长度为变量,利用搜索算法对方程组进行求解,并得到相应的结果.如当风速为12m/s时,钢桶的倾斜角度1.0405°,从上到下钢管的倾斜角度分别为1.0086°、1.0146°、1.0206°、1.0267°,浮标吃水深度0.735m,浮标游动区域半径14.4429m.针对问题二,首先将风速为36m/s的情况代入问题一建立的模型中,但是得到的结果不满足题目所给定的要求.则考虑在重物球质量一定的条件下,以浮标的吃水深度和游动区域与钢桶的倾斜角为目标,建立了一个单决策变量的多目标最优系泊模型,相比于问题一,此问的变量更多,更加难于求解,故考虑将多目标转化成单目标的问题进行求解,并继续使用搜索法对问题进行求解.最后找到了三组可行解,其中最优解是重力球的质量为2102kg.针对问题三,本文中有三个决策变量以与三个变系数,相比于前两问,无论是计算量还是计算维数,难度更大.为了求解该问,建立了一个多决策变量的多目标变系数的最优系泊系统模型,为了简便运算,我们建立了变步长的搜索算法,并最终求解得到结果,得到的一组解为: 选用了III型号的锚链,重物球质量为2800kg,锚链长度为23.4m.针对论文的实际情况,对论文的优缺点做了评价,文章最后还给出了其他的改进方向,以用于指导实际应用.关键词：系泊系统设计；力的平移定理；多目标；优化模型；搜索算法1问题的重述一个由浮标系统、系泊系统和水声通讯系统组成的近浅海观测网的传输节点.可以简化看作是一个浮于海平面的圆柱体浮标通过四节钢管装有通讯设备的钢桶,钢桶再通过锚链一个可以移动的锚,锚沉在海床上.为了不让锚被拖行要求锚链末端与锚的处的切线方向与海床的夹角不超过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符号说明：浮标所受水平方向的风力大小；：浮标所受浮力大小；：每节钢管所受浮力大小；：钢桶所受浮力大小；：第一节钢管对浮标拉力的大小<或浮标对第一节钢管拉力的大小>；：第节钢管下端所受拉力的大小<>；：钢桶所受锚链的拉力大小<或锚链上方第一段所受钢桶的拉力大小>；：锚链第段<自上而下>下端所受拉力大小<,其中为张紧状态锚链的段数>；：锚与海底间的摩擦力；：锚与海底间的最大静摩擦力；：浮标所受的水流力大小；：第节钢管所受水流力大小；：钢桶所受水流力大小；：锚与海底间的最大静摩擦系数；：海底对锚的支持力；：浮标的重力；：每节钢管的重力；：表示钢桶的重力；：锚链每段链环的重力；：锚的重力；：第一节钢管对浮标拉力与水平方向的夹角；：第节钢管下端所受拉力与水平方向的夹角<>；：钢桶所受锚链的拉力与水平方向的夹角；：锚链第段<自上而下>下端所受拉力与水平方向的夹角<>；：浮标的水平倾斜角度；：第节钢管的水平倾斜角度<>；：钢桶的水平倾斜角度；：圆柱体浮标两端的浸水长度,其中；：浮标的吃水深度；：水平方向的风速；：浮标在风向法平面的投影面积；：浮标的排水体积；：水流速度大小；：浮标的质量；：每节钢管的质量；：钢桶系统的质量；：重物球的质量；：海水密度；：锚链的线密度；：法向阻力系数；：切向阻力系数；：浮标的底面半径；：每节钢管的底面半径；：每节钢管的长度；：钢桶的底面半径；：钢桶的长度；：锚链每段链环的长度；：海水深度；：锚链某点至锚链上端顶点的弧长；：浮标浸入水下部分在水流力法平面的投影面积；：第节钢管在水流力法平面的投影面积；：钢桶在水流力法平面的投影面积；：锚链某点处所受拉力的大小；：锚链某点处所受拉力与水平方向的夹角；：平衡状态下浮标的水平距离；：平衡状态下钢管的水平距离；：平衡状态下钢桶的水平距离；：平衡状态下锚链的水平距离；：平衡状态下整体的水平距离；：锚链竖直方向上的高度；：锚不滑动时,锚链在锚点与海床的最大夹角；：单位长度上锚链所受的法向水流力大小；：单位长度上锚链所受的切向水流力大小；4.模型的建立与求解4.1 问题的分析4.1.1 系统平衡的两种状态张紧状态和松弛状态与游动区域的理解本文考虑的在风力与水流力作用下的系泊系统的设计问题,首先我们考虑一种简单情形,即海水静止且风向水平情况,此时整个系泊系统可能会因风力移动达到一个平衡的状态,由于风向恒定我们可以认为最终的平衡状态中浮标、钢管、钢桶、锚链和锚大致位于同一个平面内,如图1所示.图1 传输节点示意图在平衡状态下,各系统所受力与力矩分别达到平衡,此时浮标会与锚形成一定的水平距离.此时的平衡状态可能分为两种：一种是张紧状态,即锚链的每一段都受到了前后端的拉力从而使整个系统平衡；另一种是松弛状态,即锚链上方的一部分存在拉力,而下方部分散落在水底呈无拉力状态达到平衡.在张紧状态下达到平衡时,浮标到达某个定点,此时该点与锚得水平距离为半径以锚和浮标的水平距离为半径、锚所在位置为圆心做一个圆面,该区域即浮标的游动区域, 如图2所示.但在松弛状态下由于松弛部分锚链的可伸缩性,浮标的位置会处于一定的区域,当松弛部分锚链长度最大时浮标位于最远位置,以此时浮标与锚的水平距离为半径、锚所在位置为圆心的圆形区域即其游动区域.图2 浮标游动区域俯视示意图4.1.2 浮标的倾斜问题与对吃水深度的理解在对组成各个系统的零件受力分析的过程中,会出现无法将物体看作是均匀的情况,这样一来,物体的各力的受力点将不再是质点,为此应尽可能将受力点移动到质心上.这时,根据力的平移定理,在将各力平移到质心上去的过程中会产生力矩,且对于一个平衡的物体其力矩和应为零.在力矩平衡的作用下,我们发现平衡状态下浮标会产生一定的倾斜角度.假设在浮标的底面与海平面平行,即没有倾角时,其受力分析如图3所示.图3 浮标垂直海平面的受力分析图由图可知除浮力外,风力和钢管拉力都会产生非零力矩,且方向相同,导致力矩和不为零,从而会产生旋转,不会达到平衡状态,我们认为在平衡状态下,浮标底面相对海平面存在一个倾斜角度.吃水深度为浮标底部到海面的距离.由于浮标倾斜,导致浮标各点侵入海水的长度不同,如图4所示.假设浮标两边侵入水中的长度分别为,其中,则吃水深度即两边浸水长度的平均值.图4 倾斜浮标的浸水示意图4.2 问题1模型的建立根据对问题的分析发现,问题旨在讨论在不同情形达到平衡状态下,分析各钢管的倾斜角度、锚链形状、浮标的方程,通过已知的信息,求得系统各部分所受拉力大小与其角度情况和各部分本身的倾斜角度.4.2.1 浮标的力平衡方程通过浮标的受力分析可得,浮标受到竖直向上的浮力,水平向右的风力,重力以与钢管对其的拉力,其受力分析图如图5所示：图5浮标的受力分析图根据平面力系平衡原理,可得浮标的力平衡方程如下：在水平方向上,有1其中,为钢管对浮标拉力与水平方向的夹角,为水平风力的大小,有如下的近似公式为物体在风法平面的投影面积2在竖直方向上,有3其中,有4而5在这儿,为浮标质量,分别表示浮标两边的浸水长度<见图4>,为浮标的底面半径,为浮标的排水体积,为海水密度,为重力加速度.4.2.2 浮标在风向法平面的投影面积公式在计算风力时,需用到浮标在风向法平面的投影面积,由于浮标倾斜,该投影区域并非一标准的矩形,而是由两个部分构成,其中的下方部分为矩形,而上方部分为倾斜圆在水平面上投影的一半,即半椭圆的一半,最终得到浮标在风向法平面的投影面积如图6所示.图6 浮标在风向法平面的投影面积示意图通过平面几何关系可得,第一部分矩形的面积满足：6第二部分半椭圆的面积满足：7其中,为浮标的底面半径,分别表示浮标两边的浸水长度,表示浮标的长度.4.2.3钢管的平衡方程通过分析可得,每节钢管的受力分析图如图7所示.图7第节钢管的受力分析图根据平面力系平衡原理与力的平移定理<力矩平衡公式>,有如下方程：水平方向力平衡方程8竖直方向力平衡方程9力矩平衡方程10在这儿,为每节钢管所受浮力大小,满足11而为第节钢管下端所受拉力大小,为每节钢管的质量,为重力加速度,是第节钢管与水平方向的夹角钢管的倾斜角度,是拉力与水平方向的夹角,为钢管的底面半径,为钢管的长度.4.2.4钢桶系统的平衡方程钢桶的受力分析如图8所示.图8钢桶的受力分析图根据平面力系平衡原理与力的平移定理<力矩平衡公式>,有如下方程：水平方向力平衡方程12竖直方向力平衡方程13力矩平衡方程14在这儿,为钢桶系统的浮力,有15而表示钢管对钢桶拉力的大小,表示钢管对钢桶拉力与水平方向的夹角,表示钢桶受下方锚链拉力的大小,表示与水平方向的夹角,表示钢桶的质量,表示重物球的质量,表示钢桶的底面半径,表示钢桶的长度.4.2.5锚链系统的平衡方程由于锚链无档普通链环,我们可以将其看作无弹性悬垂线,这样可以将锚链视作柔性的.因此可用微元法的思想分析其受力平衡状态.考虑锚链其中自上而下一小段弧长<到>上的受力情况,如图9所示.图9锚链微元的受力分析图图中表示取的任意一段小弧长,表示上段锚链对所取小弧长的拉力,而表示下段锚链对小弧长的拉力；是拉力与水平方向的夹角；是拉力与水平方向的夹角.根据这一弧长微元上的受力平衡,可得出以下方程：水平方向上,有16 竖直方向上,有17其中表示单位长度锚链的质量. 将其展开可得18由泰勒公式可知,当角度变化很小时,忽略高阶无穷小的影响,故近似等于,近似等于1,从而可得如下的近似方程：19忽略高阶无穷小项,其中都为无穷小,故它们的乘积为高阶的无穷小对上式化简可得如下一阶非线性微分方程组：20且有初值条件4.2.6锚的平衡方程与平衡状态条件锚的受力直接关系着最终是否达到平衡状态,其受力分析如图10所示：图10锚的受力分析图水平方向上,有21在竖直方向上,有22由于判断整个系统平衡的关键即对锚的水平拉力是否超过锚与海底之间的最大静摩擦力,即得到平衡状态的判断条件为23又由海水深度为18,可知24其中,表示海水深度；表示浮标底面轴心到海平面的距离,分别表示浮标两边浸入水中的长度,浮标的倾斜角度,表示4个钢管在水中的高度和,表示每节钢管的长度,表示第节钢管的倾斜角度；表示钢桶的长度,表示钢桶的倾斜角度,表示锚链竖直方向的高度.4.2.7 游动区域的分析经过分析和求解可以知道在平衡状态下的各个系统水平距离分别为：浮标的水平距离25钢管的水平距离26钢桶的水平距离27锚链的水平距离28则整体的水平距离<游动半径>为294.3问题一模型的求解4.3.1 模型的求解思路通过对系泊系统各部分进行单独受力分析和力矩分析,得到了系统各处拉力大小与其倾角等未知量在平衡状态下的关系,由于方程数量较多,且相互影响,使得直接计算非常复杂.因此我们首先考虑各方程之间的关系,得到未知量满足一定规律的递推公式,对问题进行系列简化,然后再采用搜索算法进行求解.4.3.2方程的简化与递推公式的推导4.3.2.1 锚链系统微分方程的转化由于锚链是无档普通链环,实际中是环环相扣的离散系统,为了计算方便,我们首先考虑对锚链系统对应的一阶非线性微分方程组<20式>进行转换3,由与水平方向与竖直方向上的位移之间的关系30可得链环连接处的受力大小与角度的递推公式31其中,表示链环上方第一段所受钢桶带来拉力的大小,表示链环上方第一段所受钢桶带来拉力与水平方向的角度,由于其与链环对钢桶拉力互为反作用力,有；表示链环第段<自上而下>下端所受拉力大小,表示链环第段下端所受拉力与水平方向的夹角,为链环的段数.4水平力平衡公式的推导与简化整个系泊系统的每个部分均在两个水平力的作用下平衡,有32其中风力由浮标两边的入水深度决定.4.3.2.3 竖直方向力平衡公式的推导与简化由竖直方向上受力平衡可得到如下递推公式其中浮标的浮力也由浮标两边的入水深度决定.4.3.2.4 系统各段所受拉力大小与其角度的计算公式将各拉力竖直方向分力与水平分力相比,可得相应角度正切值的计算公式如从而可得平衡状态下依赖于的各段拉力角度与大小的计算公式.4各个部件在竖直方向的高度的计算1找出浮标的竖直方向高度,这里记为<>.如图所示：图11. 浮标的竖直方向高度示意图图中,.在直角和直角中,因为同位角相等,所以这两个三角形另外两个角相等,即：35所以36即:37其中,表示浮标的半径,为1.2钢管的竖直方向高度,这里记为.根据浮标的平衡方程、水平方向的方程和竖直方向的方程相除得：38即可以求出,这样再带入浮标水平方向平衡方程则可求出.再把和的值带入钢管的平衡方程,我们可以求出钢管的拉力与其倾斜角度和钢管的倾斜角,又可以根据钢管的几何关系如图6所示：39其中表示四根钢管的总高度,表示第节钢管的高度,表示第节钢管的倾斜角度.3钢桶的竖直方向高度,这里记为.把上个环节的计算出来的结果带入钢桶系统的平衡方程,同样也可以计算出钢桶倾斜角度,锚链对钢桶的拉力与其方向与轴的夹角.则可以得出钢桶高度的公式：40其中为钢桶的长度,它是一个常数,大小为.4钢桶的竖直方向高度,这里记为.根据已求得的钢桶对锚链的拉力与其方向,把这对初始值带入锚链平衡微分方程,可解出函数关系.就可以表示出锚链的高度：41其中,表示锚链一微段的长度,表示那一微段的倾斜角度.故整个系泊系统总高度为424.3.2.6吃水深度和的值范围的确定 根据题目可知浮标的高度为2,则有：,为较长的一条高. 只有当浮标在海水中时,浮标的重力与浮标的浮力二力平衡,即：43可以计算出此时.那么可以进一步确定范围：,444.3.3 算法步骤1给定一组浮标吃水深度值,得到,可以依次计算得到,；2判断如果同时满足,取较小值,比如0.1和两个条件,则输出结果,并停止计算；如果不满足,则取步长有并转到步骤1.4.3.4 计算结果与分析1当风速为12时,得到的结果如下表所示：表1.风速为12时的求解结果h1h2海水深度锚链水平倾角钢桶的倾角所用无档链环个数0.68m0.79m17.9042m3.4838°1.0405°144钢管1的倾角钢管2的倾角钢管3的倾角钢管4的倾角游动区域半径浮标的吃水深度1.0086°1.0146°1.0206°1.0267°14.4429m0.735m得到锚链的形状如图12所示：图12. 风速为12时锚链的形状以上求解得到的结果满足题目的要求,且所用的链环个数为144个,小于链的总个数,说明达到平衡时锚链并未完全绷紧.2当风速为24时,得到的结果如下表所示：表2. 风速为24时的求解结果h1h2海水深度锚链水平倾角钢桶的倾角所用无档链环个数0.71m0.79m17.9117m11.7855°3.9299°185钢管1的倾角钢管2的倾角钢管3的倾角钢管4的倾角游动区域半径浮标的吃水深度3.814°3.8356°3.8575°3.8797°17.2252m0.75m得到锚链的形状如图13所示：图13. 风速为24时锚链的形状以上求解得到的结果满足题目的要求,且所用的链环个数为185个,小于链的总个数,说明达到平衡时锚链并未完全绷紧.4.4问题2的模型分析与建立4.4.1 风速为36时系泊状态与其分析当风速达到36时,若不考虑角度的限制,可通过问题一的递推式与搜索算法求得如下结果：表3. 风速为36时的求解结果h1h2海水深度锚链水平倾角钢桶的倾角所用无档链环个数0.759m0.7810m17.7796m17.7796°8.1268°210钢管1的倾角钢管2的倾角钢管3的倾角钢管4的倾角游动区域半径浮标的吃水深度7.8999°7.9424°7.9853°8.0287°18.3712m0.77m得到锚链的形状如图14所示：图14.风速为36时锚链的形状此时,锚链末端与锚的处的切线方向与海床的夹角为17.7796度,超过了16度,钢桶的倾斜角度<钢桶与竖直线的夹角>为8.1268度,超过5度,已经不能够达到平衡状态,因此需考虑调节重物球质量,使其重新达到平衡状态.4.4.2 基于重物球质量的定常系统最优系泊模型当风速达到36时,系泊系统已经不再满足角度的要求,从而不再处于平衡状态,因此我们考虑调整重物球质量.在此我们假设海水静止,考虑水平风速、水深、锚链型号、长度等均为定常数,基于重物球质量的选择,建立了单决策变量的多目标最优系泊模型.4.4.2.1目标函数一浮标的吃水深度尽可能小：45其中吃水深度,即浮标两边浸水长度的平均值； 二游动区域尽可能小,即游动半径尽可能最小：46其中游动半径<见29式>；三钢桶的倾斜角度尽可能小；47其中表示钢桶与水平面的夹角.4.4.2.2 约束条件<一>海水的深度约束,即平衡状态下整个系泊系统的竖直高度之和应等于水深,有 48其中表示浮标底面轴心到海平面的距离,分别表示浮标两边浸入水中的长度,浮标的倾斜角度,表示4个钢管在水中的高度和,表示每节钢管的长度,表示第节钢管的倾斜角度；表示钢桶的长度,表示钢桶的倾斜角度,表示锚链竖直方向的高度.二锚链底端拉力的夹角约束,即锚链在锚点与海床的夹角不超过16度,有49由问题1可知,满足如下微分方程组50其中表示锚链在锚点与海床的夹角,表示使得锚不滑动时,锚链在锚点与海床的最大夹角.三平衡约束条件,即达到系统各部分力和力矩平衡的各等式约束条件,包括式？综上：我们建立了基于重物球质量的多目标最优系泊模型如下51524.5问题2 的模型求解4.5.1多目标的转换对于问题2,这是一个多目标的优化问题,不能够直接进行求解,需要把多目标转化为单目标进行求解.具体方法为：先计算单目标5354得到倾角的最优值,再将作为约束条件加入到优化问题里面,这样问题转化为双目标优化问题.55<56> 再考虑吃水深度的优化,利用类似的方法得到其最优值作为约束条件,即：<57><58>4.5.2 可行解的范围 因为该问所给的风速变大,在初始的条件下不能再得到可行解了.所以这里除了给定、的值,还需要给定重物球的质量.在问题一种重力球质量为1200已经不能满足题目要求,所以重物球的重力应该满足：.为了找到重物球重力的上限,首先假设只有重物球的重力作用,且浮标完全淹没在水中.此时,对浮标进行受力分析,得到的受力分析图如图15所示：图15.浮标的受力分析图根据受力分析,当浮标处于平衡状态时有：<59>其中,；由于浮标完全被淹没,所以有：则：此时,.因此设定重物球重力的取值范围是.当重物球重力的值接近上限值时,浮标接近被完全淹没,此时浮标受风吹的面积小,即风力小,整个系统在水平方向只受风力和锚上的摩擦力的作用,且这两个力是一对平衡力,当风力较小时,系统将不会在水平方向上移动,当浮标完全沉没时,系统将完全不在水平方向移动,这时钢桶倾斜角度达到最小为0度,锚链处于松弛状态,且此时为平衡状态.4.5.3问题2 的具体算法步骤求解的基本思路同问题一是类似的,还是采用搜索算法.但是由于变量个数增加了一个,搜索的范围和计算量将大大增加,如果还采用问题一的定步长搜索的话,计算时间和结果可能难以保证.故采用变步长搜索算法,具体步骤为：1在搜索参数的有效范围内,给定一个较大的步长比如,与问题一的求解方法一样,得到对应的；2判断如果同时满足条件,取较小值,比如0.1和,则输出结果并结束计算；如果不满足,则记下最接近判断条件的一组；3以为新的初始点,选定一个较小的步长比如,按照上述方法再次计算,如果满足判断条件则停止；否则可以类似方法继续缩小步长,总可以找到可行解.4.5.4问题2 的计算结果根据上述的求解思路和方法,编写程序进行搜索求解,限定搜索的时间,得到以下结果：第一种结果：表4.问题2的第一种求解结果重物球的质量h1h2海水深度1968kg1m1m17.9081m锚链水平倾角钢桶的倾角所用链环数/总链环数游动区域半径13.7633°4.3475°207/21018.2717m在得到以上结果的情况下,将参数代入可以得到锚链的形状为：图16.第一种求解结果下锚链的形状其他结果：表5.问题2的第二种求解结果重物球的质量h1h2海水深度1969kg1m1m17.9038m锚链水平倾角钢桶的倾角所用链环数/总链环数游动区域半径12.7107°4.3475°210/21018.3094m表6.问题2的第三种求解结果重物球的质量h1h2海水深度2102kg1.04m1.04m17.9047m锚链水平倾角钢桶的倾角所用链环数/总链环数游动区域半径12.0327°3.9415°210/21018.1739m问题2的第二、第三种求解结果满足题目的要求,说明对于第二问的解题思路和方法是正确且可行的,如果按照此种方法就能找到最优的解,但由于算法的数据处理量过于庞大以与时间的限制,这里就不一一进行求解了.4.6问题3的模型建立与求解4.6.1最优系泊问题问题一和问题二的模型都是在水深、风速恒定、水流静止的特殊情况下建立的,而对于此问,题目给定水深、风速、水流的范围,考虑实际系泊系统中水流速度的影响,在水速、风速、水深不定的情况下,基于锚链型号、长度与重物球质量的选择,建立了多决策变量的多目标最优系泊模型.所以决策目标和第二题相同,考虑基于锚链型号、长度与重物球质量的变系数,那么模型的限制条件因不同情况而不同,下面以当水平风速和水流力同向的情况分析.问题3当水平风速和水流力同向的情况4.6.2.1浮标的力平衡方程根据理论力学的相关知识,对平衡时的浮标进行受力分析,得到浮标的受力分析图如图17所示：图17.浮标的受力分析图浮标受到竖直向上的浮力；水平向右的风力；浮标所受水流力；竖直向下的浮标重力；浮标受下方钢管的拉力.其中表示浮标的倾斜角度,表示拉力与轴方向的夹角.根据平面力系平衡原理,轴与轴方向上受力达到平衡,可得出浮标的平衡方程：水平方向：<60>其中,有近海风载荷的近似公式：为物体在风法平面的投影面积<61><62>近海水流力的近似公式：为物体在水流法平面的投影面积由问题一中浮标在风向法平面投影面积的分析可知,投影面积可分为矩形和椭圆两部分,可得：<63>竖直方向：<64>其中,有,<65><66><67>其中,、分别表示浮标左、右边进入水中桶壁的长度4.6.2.2钢管的平衡方程对第节钢管对其进行分析,得到受力分析图如图18所示：图18. 第节钢管的受力分析图为第节钢管的在水中的浮力；为第节钢管受到上方物体的拉力；：为第节钢管受到下方物体的拉力；钢管所受水流力；为第节钢管的重力；是第节钢管与轴的夹角钢管的倾斜角度；是拉力与轴的夹角；是拉力与轴的夹角.根据平面力系平衡原理和力的平移定理,与方向上受力平衡且力矩和为零,可得出钢管的平衡方程：力平衡：<水平方向><68>其中,有<69><70>力平衡：<竖直方向><71>其中,有,<72>力矩平衡：<73>4.6.2.3钢桶系统的平衡方程对于通信系统的钢桶,得到受力分析如图19所示：图19.钢桶的受力分析图表示上面第4节钢管对钢桶的拉力；表示拉力与轴的夹角；表示钢桶受下方锚链的拉力；钢管所受水流力；表示力与轴的夹角；表示重球对钢桶的重力；表示钢桶的自身所受的重力；表示钢桶在水中所受的浮力；为钢桶的倾斜角度钢桶与轴的夹角.根据平面力系平衡原理和力的平移定理,与方向上受力平衡且力矩和为零,可得出浮标的平衡方程：力平衡：<水平方向><74>其中,有<75><76>力平衡：<竖直方向><77>其中,有,<78>力矩平衡：<79>4.6.2.4锚链系统的平衡分析假设锚链是无弹性悬垂线,这样可以将锚链简化看作是柔软的,利于模型的简化.再考虑其中一小段弧长<到>上的受力情况,在有水流力的情况下,如图20所示：图20.锚链的受力分析图图中表示取的任意一段小弧长,表示下段锚链对所取小弧长的拉力,而表示上段锚链对小弧长的拉力；是拉力与水平方向的夹角；是拉力与水平方向的夹角.根据这一微段的平衡,可得出以下力系平衡方程：水平方向上,有<80> 竖直方向上,有<81>其中表示单位长度锚链的质量,和分别为单位长度上锚链所受的法向和切向水流力.对上式化简可得如下一阶非线性微分方程组：<82>通过查找资料4,和按下列公式来计算：<83><84>其中,和分别为法向和切向阻力系数,为水的密度,为水流速度综上：我们建立了基于重物球质量的多目标变系数最优系泊模型如下<85><86>模型求解在水流速度为1.5,风速为36这种情况下,相比问题二多了两个决策变量：锚链的型号与其长度,增加了三个变系数：水流速度、风速、水深.所以求解这个模型依然通过给定重物球重力,锚链长度,锚链型号,和.用问题二同样的思路去搜索解,最后搜索到结果为：表7.问题三的第一组解锚链型号III重物球的质量水流速h1h2海水深度锚链水平倾角钢桶的倾角2800kg1.5m/s1.31.316.0057m5.284°6.4744°钢管1的倾角钢管2的倾角钢管3的倾角钢管4的倾角所用链环数/总链环数锚链长度游动区域半径6.1184°6.172°6.1874°6.2029°195/19523.4m21.1156m表8.问题三的第二组解锚链型号IV重物球的质量水流速h1h2海水深度锚链水平倾角钢桶的倾角2960kg1.5m/s1.4m1.4m18.7723m0.0331°5.887°钢管1的倾角钢管2的倾角钢管3的倾角钢管4的倾角所用链环数/总链环数锚链长度游动区域半径5.5703°5.6178°5.6306°5.6434°131/16825.2m26.5694m得到的锚链形状如图所示：图21.第一组求解结果下锚链的形状图22.第二组求解结果下锚链的形状5模型的推广与改进方向本文建立了单一决策变量的多目标以与多决策变量的多目标优化模型,推导出系泊系统各个部分的平衡方程,巧妙的运用搜索算法,搜索算法极大的减少了运算量,提高了运算效率.根据建立的模型可以求解并设计一个最优的系泊系统.本文求解模型的方法可以推广应用到各种参数较多,方程数较多的方程组求解.以与本文建立的模型可以很好的运用于各种链状的物体受力分析的研究上去.但本文对更多优解的研究还不够深入,可以继续求解更有结进行研究.6模型的优缺点通过对系泊系统各个部分严谨的受力分析,得出了比较完美的受力平衡关系式,进而建立了比较完善的模型,思路严晰,结构严谨.只是由于运算量太大,没有找到更好的求解方法.如果时间足够充裕还可以对模型的求解进行更加深入的研究.7.参考文献1 姜启源. 数学模型第三版M. ：高等教育,1999.2 韩中庚. 数学建模方法与其应用第二版M. ：高等教育,2009.3 孙宁松,海上移动式平台锚泊定位系统锚索链受力分析J. 中国海洋平台,2008,23<2>4 赫春玲,滕斌,不均匀可拉伸单锚链系统的静力分析J. 中国海洋平台,2003,1845 滕斌,赫春玲,韩凌,Chebyshev多项式在锚链分析中的应用J. 中国工程科学,2005,718附录1、第一问和第二问程序%20160911%系泊系统的第一问和第二问%搜索方法clearclcclose allg=9.8;%重力加速度v=36;%风速rou=1025;%海水密度ma=1000;da=2;ha=2;%浮标lb=1;mb=10;db=0.05;%钢管mc=100;lc=1;dc=0.3;%钢桶ball=1200;%重物球ld=0.105;md=7;Ld=22.05;nd=Ld/ld;md=md*ld*g;%链环me=600;%锚Ffb=rou*g*pi*db2/4*lb;%钢管的浮力mb=mb*g;%钢管的重力t=1;for ball=1200%1200:2:3000for h1=0.759%0.9:0.01:1.1%0.3:0.05:1%0.31:0.01:1.5%0.68:0.0005:0.72 for h2=0.781%h1:0.01:1.5%0.75:0.0005:0.79 L=h1; R=0; %-浮标- S1=da*<ha-<h1+h2>/2>*da/sqrt<<h2-h1>2+da2> S2=pi*da2/8*<h2-h1>/da; S=S1+S2;%风载荷 Fw=0.625*S*v2;%风力 V=pi*da2/8*<h1+h2>%吃水体积 Ffa=rou*g*V;%浮标的浮力 Ta1=Ffa-ma*g;%竖直方向 Ta2=Fw;%水平方向 Ta=sqrt<Ta12+Ta22>%浮标的总拉力 thetaa=atand<Ta1/Ta2> %浮标拉力的角度 alphaa=atand<<h2-h1>/da>%浮标的倾斜角 %-浮标- %-钢管- for k=1:4 Tb2<k>=Ta2;%水平方向 if k=1 Tb1<k>=Ta1+Ffb-mb;%竖直方向 else Tb1<k>=Tb1<k-1>+Ffb-mb;%竖直方向 end thetab<k>=atand<Tb1<k>/Tb2<k>>%拉力角度 Tb<k>=sqrt<Tb1<k>2+Tb2<k>2>%拉力大小 if k=1%钢管的倾斜角 alphab<k>=atand<<Ta*sind<thetaa>+Tb<k>*sind<thetab<k>>>. /<Tb<k>*cosd<thetab<k>>+Ta*cosd<thetaa>>> else alphab<k>=atand<<Tb<k-1>*sind<thetab<k-1>>+Tb<k>*sind<thetab<k>>>. /<Tb<k>*cosd<thetab<k>>+Tb<k-1>*cosd<thetab<k-1>>>> end if k=1 xb<k>=ld*cosd<alphaa>%钢管k的x长度 else xb<k>=ld*cosd<alphab<k-1>>%钢管k的x长度 end R=R+xb<k>%目前游动总长度