浅水域生存工况下系泊浮筒系统三维运动.docxVIP

浅水域生存工况下系泊浮筒系统三维运动3张挺孙明光(中山大学, 广州 510275)摘要本文提出一个浅水域生存工况下, 多锚链系泊浮筒系统三维运动的工程计算模型, 对链系采用射的法计算, 使其与浮筒运动计算相匹配。对不同的风波流方向及链刚度, 计算了浮筒 的运动及链力, 结果显示本文的模型是合理的。关键词生存工况 浅水域 浮筒 链系 三维运动前言1悬链腿单点系泊系统, 由于其适应性广, 技术成熟, 安装容易, 对水深和油轮大小不是非常敏感, 因而在世界上被大量使用。 在我国, 也开始使用这类系统。在一般情况下, 认为系统处于作业工况 (即有油轮系泊) 时, 其系泊链张力最大。 这方面 已有不少的工作。 但是, 我们对一些系统的初步估算说明, 在浅水情形, 系统处于生存工况下, 虽然油轮撒离了, 但浮筒却有更激烈的运动, 从而导致系泊链会发生更大的张力。实践证实了这种情况存在的可能性。目前, 有关这方面的问题, 研究结果很少。文献1 仅研究多锚 链系泊浮筒的纵荡运动, 文献 2 则考虑了纵荡与升沉的非耦合运动, 文献 3 虽考虑了纵 荡、升沉及纵摇的耦合运动, 但属于弱非线性问题, 不适合浅水及生存工况的计算, 另外, 它难以 推广到三维运动情形, 文献4 研究了船舶在波中大振幅运动问题, 计算时湿表面采用瞬 时值, 文献 5研究了波流共存情形下的 CD 及 CM值计算方法, 文献 6用时域模拟了二阶波与系泊柱之间的相互作用, 文献 7 研究了物体靠近海底时, CD 、CM 值的变化。 我们注意到浮筒系泊系统, 在浅水生存工况条件下, 会发生大振幅三维运动, 前后链的运动十分复 杂, 前链处于极度张紧状态, 后链则处于松弛状态; 同时由于锚链系统受力大, 刚变曲线会出现高度非线性的状态, 甚至会超过平方关系很多; 另外, 由于运动是强非线性, 各运动分量间有强烈的耦合, 而且上浪等现象也会影响运动。 在浅水情况下, 附加质量及阻尼除了依赖于 频率外, 还与位置有关。 本文针对上述问题, 提出一个工程计算模型, 一方面要进行简化, 另 一方面也保持一定的准确度。计算模型及流体作用力2对于油轮撤离后, 在浅水域中处于生存工况下的多锚链的浮筒系泊系统 ( 图 1) , 现建立其 三维六自由度的非线性计算模型, 应用时域模拟来获得系统的运动情况。211系泊链三维动力学方程海洋工程第 15 卷2系统中系泊链有 6 条 ( 或 8 条) , 均有足够的长度, 在任何时刻都不出现拔桩现 象。因此, 系泊链可分为悬链部分和触底的 平直部分, 海底的摩擦相对影响较小, 可忽 略不计。 另外, 假设海流是定常的, 不随深度变化, 且没有垂直向分量。如图 2 所示, 以悬链的触地点为原点,建立直角坐标系 O X Y Z , O Z 轴垂直向上,_ _ _O X Y 平面位于海底, i , j , k 为单位矢量。_ _ _设, u , v , w 为悬链任意点 P处的一组点处的悬_直角坐标单位矢量, 其中 u 为 P__链切向, Υ为水平面与 u 之间的夹角, ?为 u_在水平面上投影与矢量 i 之间的夹角, 于图 1悬链腿系泊系统图 2 地面及悬链上的坐标系是悬链单位向量组与地面单位向量组有下述关系__co sΥ sin ?j +__sin Υ k_u = co sΥco s?i +__v = -sin ? i + co s?j_(1)___w = -sin Υco s?i -sin Υ sin ?j + co sΥ k_设 ds u 是位于 P 点处长度为 ds 的悬链元, 它的长度在 X , Y , Z 方向上的投影为:/? 1994-2015 China Academic Journal Electronic Publishing House. All rights reserved. 第 1 期浅水域生存工况下系泊浮筒系统三维运动3__dsco sΥco s?dx =dsu r i =__(2)dy =dsu r j =_dsco sΥ sin ?_dz =dsu r k =ds sin Υ悬链拉伸后的长度 ds’可表为 ds’= ds (1+ Ε) , 其中 Ε= T ?EA , 这里 T为链张力, E 为弹性模量, A为截面积, 根据牛顿定律, 悬链的三维动力学方程可写为:52 x5 T5 Υ5 t2 =5 s co sΥ- T sin Υ 5 s +F u co sΥ(1 +Ε) -F w sin Υ(1 +Ε)M52 y5 ?M 5 t2 = - T co s? 5 s -F v (1 +Ε)(3)52 z5 T5 Υ5 t2 =5 s sin Υ+ T