一类系泊系统的设计问题
Design of a Class of Mooring Systems
DOI: 10.12677/APP.2021.111002, PDF, HTML, XML, 下载: 557  浏览: 1,782 
作者: 刘春江:中国人民解放军陆军炮兵防空兵学院,安徽 合肥;李宏欣, 易绵竹:信息工程大学,河南 洛阳
关键词: 系泊系统受力分析最优状态静力平衡微元化Mooring System Stress Analysis Optimal State Static Balance Microelement
摘要: 系泊系统的设计主要确定系泊系统各个部件的参数问题。本文把这些物理部件集中看成一点进行受力分析,分析这些部件之间的关系,各个部件处于什么状态下使整个系泊系统处于最优状态,根据静力平衡条件列出平衡方程,给出各个节点的坐标以及各个钢管和钢桶与竖直方向的夹角之间的关系方程,对锚链进行微元化,获得系泊系统状态下各个参数及位置关系,使系泊系统处于最优状态。
Abstract: The design of mooring system mainly determines the parameter problem of each component of mooring system. In this paper, we take these physical components as a point for force analysis, analyze the relationship between these components, make the whole mooring system in the optimal state under which state each component is in, and list the equilibrium equation according to the static equilibrium condition.
文章引用:刘春江, 李宏欣, 易绵竹. 一类系泊系统的设计问题[J]. 应用物理, 2021, 11(1): 9-16. https://doi.org/10.12677/APP.2021.111002

1. 引言

当今世界对石油资源的需求增大,海上石油资源开发已逐渐成为新的能源发展方向,因此,为了满足日益增长的能源需求,各国都将油气资源的开发的重点转移到深海。海上平台广泛应用于钻井,采油,加工处理和储存。为了保证海上平台在大多数海况下都可以正常工作,系泊系统的作用至关重要。因此,需要相关的理论来更好的解决相关设备中遇到的问题。目前有多种用于浮体结构定位的系泊系统,如果按系泊方式划分,可以分为悬链线锚腿系泊 [1] [2] [3]、单锚腿系泊 [4] [5] [6]、悬链线系泊 [7] [8] [9] [10] [11]、张力腿系泊 [12] [13] [14] [15] [16]、系缆桩–缓冲系泊 [17] - [23] 等。一般来说,系泊方式是根据系泊力的大小、水深、系索长度、是否配置中间浮体、地形和海床等条件来确定的。为适应深海平台开发的需要,近几年来,深海系泊理论和技术的研究得到了蓬勃发展。

一般近浅海平台观测网的传输节点由浮标系统、系泊系统和水声通讯系统组成(如图1所示)。模型传输节点的浮标系统可简化为圆柱体。系泊系统由钢管、钢桶、重物球、电焊锚链和特制的抗拖移锚组成,锚链选用无档普通链环,钢管共n节组成,水声通讯系统安装在密封圆柱形钢桶内。钢桶上接第n节钢管,下接电焊锚链。钢桶竖直时,水声通讯设备的工作效果最佳。若钢桶倾斜,则影响设备的工作效果。钢桶的倾斜角度(钢桶与竖直线的夹角)超过一定度时,设备的工作效果较差。为了控制钢桶的倾斜角度,钢桶与电焊锚链链接处可悬挂重物球。

Figure 1. Transmission node diagram (structure module diagram only, size ratio not considered)

图1. 传输节点示意图(仅为结构模块示意图,未考虑尺寸比例)

在该系泊系统中,包括钢管、钢桶、重物球、电焊锚杆和特制的抗拖移锚组成。锚链末端与锚的连接处的切线方向与海床的夹角一般不超过β˚,否则,锚会被拖行。而钢桶的倾斜角度也被限制在0˚到α˚之间,0˚其工作效果最佳,倾斜角度愈大工作效果愈差,系泊系统的设计问题就是确定锚链的型号、长度和重物球的质量,使得浮标的吃水深度和游动区域及钢桶的倾斜角度尽可能小。关于系泊系统 [24] [25] [26] [27] [28] 有许多作者有研究,如潘斌 [14] 等研究了单腿浮标在外力作用下,其系泊系统的静力计算方法。系统包括:一个浮筒,一根锚泊线,一个锚。锚泊线可由多段不同重量和尺寸的索或链组成。计算中考虑了锚泊线的伸长变形以及流速沿水深的变化,孙宁松 [10] 对平台配备的定位锚泊系统进行分析和设计。

本文主要探讨了如何确定系泊系统浮标、钢管、钢桶、重物球、电焊锚链和特制的抗拖移锚等参数,锚链选用无档普通链环,钢管共n节组成,水声通讯系统安装在密封圆柱形钢桶内。钢桶上接第n节钢管,下接电焊锚链。钢桶竖直时,在上述情况下如何确定参数使水声通讯设备的工作效果最佳问题。

2. 预备知识

为了能更好分析、描述系统个部件之间的相互影响,个部件的受力情况,不妨对个部件做如下的假设。

1) 假设钢管之间用铰链链接可自由转动。

2) 假设风是沿平行海平面的方向吹来。

3) 假设锚链和重物球的浮力影响可以忽略。

4) 假设忽略其它原因(如温度,湿度等)对系统影响。

5) 假设浮标不会出现因为风力等因素而被吹翻的情况。

为了讨论方便,钢管由4节组成依次把浮标、钢管(每节)、钢桶、重物球、猫链、猫的质量记作: m i ( i = 1 , 2 , 3 , 4 , 5 , 6 ) ,管对浮标的拉力与竖轴的夹角记作 θ 10 ,钢管从上往下依次与竖轴的夹角记作 θ i + 1 i ( i = 1 , 2 , 3 ) ,钢桶与竖轴的夹角 θ 54 ,浮标高度为a;浮标的吃水深度为h;钢管长度为l;钢桶的高度为 l t ;锚链的长度为s。

3. 受力分析

浮标受的力由下几部分组成:在吃水深度为h的情况下浮标受到的浮力为 f = 2 R h ρ g ;浮标所受海风的应力为 f = 0.625 v 2 2 R ( a h ) ;浮标自身的重力为 m 1 g ,将浮标、钢管、钢桶、重物球看成一个整体,整体受力分析,以浮力 f 方向向上看作y轴正向,海风应力 f 与浮力垂直,看作x轴正向,钢管对浮标的作用力 f 1 看作是外力,分解为x轴上的力 f 1 x 与y轴上的力 f 1 y ,系统整体受力平衡,因此有如下关系:

f 1 x = f (3.1)

f 1 y = f m 1 g (3.2)

其中 m 1 g 为浮标的重力,当浮标的吃水深度为h,浮标的半径R,浮标的体积v及浮标的重力已知时,通过上方程计算得到钢管对浮标的作用力 f 1 ,当系泊系统处于平衡状态时,系统节点之间受到作用力与反作用力大小相等方向相反。从而钢管对浮标的作用力 f 1 与浮标对钢管的拉力大小相等,方向相反。

3.1. 对钢桶受力情况进行分析

下面分析钢桶受力情况(如图2),将各力直角系下的表示,钢桶受力有重物球的重力 G = ( 0 , m 4 g ) ;钢桶的重力 G = ( 0 , m 3 g ) ;钢桶的浮力 f = ( 0 , ρ V g ) V 是钢桶的体积,猫链的拉力 f = ( f x , f y ) 。力 G + f = ( f x , m 4 g f y ) 的作用点坐标为: B ( l sin φ , l cos φ ) ,l是桶的长度, φ 是桶的倾斜角;力 G + f 的作用点坐标为: A ( l 2 sin φ , l 2 cos φ ) ;应用虚位移原理得:

( G + f ) δ A + ( G + f ) δ B = 0 (3.3)

Figure 2. Steel barrel stress

图2. 钢桶受力

δ A = ( l 2 cos φ δ φ , l 2 sin φ δ φ ) δ B = ( l cos φ δ φ , l sin φ δ φ ) 代入上方程(2.3)化简得:

l 2 m 3 g sin φ + l 2 ρ V g sin φ + l f x cos φ l m 4 g sin φ l f y sin φ = 0 (3.4)

由此可解得 φ = arctan 2 f x m 3 g ρ V g + 2 m 4 g + 2 f y

3.2. 对钢管进行受力分析

与浮标链接的钢管必然存在一对作用力与反作用力,两力大小相等方向相反,钢管与钢管链接处必然存在一对作用力与反作用力,两力大小相等方向相反;我们这里主要考虑它的力的大小对系统影响;不妨假设钢管、浮标均在一个平面内。

第一节钢管受力情况讨论:

第一节钢管的受力有:钢管的重力 m 2 g ;钢管的浮力 ρ V g ;第二节对第一节钢管的拉力 f 21 = ( f 21 x , f 21 y ) ;受力状况如图2所示:

钢管的重力与浮力可以看作集中在 A ( k 2 sin θ 1 , k 2 cos θ 1 ) 点上,其中k为钢管的长度, φ 1 第一节钢管的倾斜角,A点受力 ( 0 , m 2 g + ρ V g ) B ( k sin θ 1 , k cos θ 1 ) 是第二与第一节钢管的接点,B点受力 f 21 = ( f 21 x , f 21 y ) V 是钢管的体积,应用虚位移原理得:

( 0 , m 2 g + ρ V g ) δ A + ( f 21 x , f 21 y ) δ B = 0 (3.5)

δ A = ( k 2 cos θ 1 δ θ 1 , k 2 sin θ 1 δ θ 1 ) δ B = ( k cos θ 1 δ θ 1 , k sin θ 1 δ θ 1 ) 代入上方程化简得:

k 2 m 2 g sin θ 1 δ θ 1 + k 2 ρ V g sin θ 1 δ θ 1 k f 21 x cos θ 1 δ θ 1 + k f 21 y sin θ 1 δ θ 1 = 0 (3.6)

由此可解得 θ 1 = arctan 2 k f 21 x k m 2 g + k ρ V g + 2 k f 21 y

第二节钢管受力情况讨论:

第二节钢管的受力有,钢管的重力 m 2 g ;钢管的浮力 ρ V g ;第三节对第二节钢管的拉力 f 32 = ( f 32 x , f 32 y ) 受力状况,钢管的重力与浮力可以看作集中在 A ( k 2 sin θ 2 , k 2 cos θ 2 ) 点上,其中k为钢管的长度, φ 2 第二节钢管的倾斜角,A点受力 ( 0 , m 2 g + ρ V g ) B ( k sin θ 2 , k cos θ 2 ) 是第三与第二节钢管的接点,B点受力 f 32 = ( f 32 x , f 32 y ) ;受力状况类似如图2所示,应用虚位移原理得:

( 0 , m 2 g + ρ V g ) δ A + ( f 32 x , f 32 y ) δ B = 0 (3.7)

δ A = ( k 2 cos θ 2 δ θ 2 , k 2 sin θ 2 δ θ 2 ) δ B = ( k cos θ 2 δ θ 2 , k sin θ 2 δ θ 2 ) 代入上方程化简得:

k 2 m 2 g sin θ 2 δ θ 2 + k 2 ρ V g sin θ 2 δ θ 2 k f 32 x cos θ 2 δ θ 2 + k f 32 y sin θ 2 δ θ 2 = 0 (3.8)

由此方程可解得 θ 2 = arctan 2 k f 32 x k m 2 g + k ρ V g + 2 k f 32 y

第三节钢管受力情况讨论:

第三节钢管的受力有,钢管的重力 m 2 g ;钢管的浮力 ρ V g ;第四节对第三节钢管的拉力 f 43 = ( f 43 x , f 43 y ) ,第三节受力:钢管的重力与浮力可以看作集中在 A ( k 2 sin θ 3 , k 2 cos θ 3 ) 点上,其中k为钢管的长度, φ 3 第三节钢管的倾斜角,A点受力 ( 0 , m 3 g + ρ V g ) B ( k sin θ 3 , k cos θ 3 ) 是第四与第三节钢管的接点,B点受力 f 43 = ( f 43 x , f 43 y ) ;受力状况类似如图2所示,应用虚位移原理得:

( 0 , m 2 g + ρ V g ) δ A + ( f 43 x , f 43 y ) δ B = 0 (3.9)

δ A = ( k 2 cos θ 3 δ θ 3 , k 2 sin θ 3 δ θ 3 ) δ B = ( k cos θ 3 δ θ 3 , k sin θ 3 δ θ 3 ) 代入上方程化简得:

k 2 m 2 g sin θ 3 δ θ 3 + k 2 ρ V g sin θ 3 δ θ 3 k f 43 x cos θ 3 δ θ 3 + k f 43 y sin θ 3 δ θ 3 = 0 (3.10)

由此方程可解得: θ 3 = arctan 2 k f 43 x k m 2 g + k ρ V g + 2 k f 43 y

第四节钢管受力情况讨论:

第四节钢管的受力有,钢管的重力 m 2 g ;钢管的浮力 ρ V g ;钢桶对第四节钢管的拉力 f 4 = ( f 4 x , f 4 y ) ,第四节受力:钢管的重力与浮力可以看作集中在 A ( k 2 sin θ 4 , k 2 cos θ 4 ) 点上,其中k为钢管的长度, φ 4 第四节钢管的倾斜角,A点受力 ( 0 , m 3 g + ρ V g ) B ( k sin θ 4 , k cos θ 4 ) 是钢桶与第三节钢管的接点,B点受力 f 4 = ( f 4 x , f 4 y ) ;受力状况类似如图2所示,应用虚位移原理得:

( 0 , m 2 g + ρ V g ) δ A + ( f 4 x , f 4 y ) δ B = 0 (3.11)

δ A = ( k 2 cos θ 4 δ θ 4 , k 2 sin θ 4 δ θ 4 ) δ B = ( k cos θ 4 δ θ 4 , k sin θ 4 δ θ 4 ) 代入上方程化简得:

k 2 m 2 g sin θ 4 δ θ 4 + k 2 ρ V g sin θ 4 δ θ 4 k f 4 x cos θ 4 δ θ 4 + k f 4 y sin θ 4 δ θ 4 = 0 (3.12)

由此方程可解得 θ 4 = arctan 2 k f 4 x k m 2 g + k ρ V g + 2 k f 4 y

3.3. 对锚链进行受力分析

锚链在水中必然是呈弯曲状态,在平衡状态下锚链的形状为一悬链线,是双曲余弦函数,最低点处受水平向左的拉力H,右悬挂点处受一个斜向上的拉力T,设T和水平方向夹角为 θ ,锚链的质量为 m 5 ,锚链的线密度为 ρ l ,s为锚链的长度,则猫链的重力为 G = ρ l s ,设猫链上任意一点上的切线的切角为 θ

根据力学原理可知,H、G和T三力平衡,可知 tan θ = G H ,假设水平张力在猫链上处处相等,对于任意一段猫链 l s ,该平衡均成立, tan θ = ρ l l s H ,而 tan θ = d y d x ,对该式 tan θ = ρ l l s H 取微分,则有

d ( tan θ ) = ρ l H d ( l s ) = ρ l H ( d x ) 2 + ( d y ) 2 = ρ l H 1 + ( tan θ ) 2 d x (3.13)

对方程(2.13)两边积分得

d ( tan θ ) 1 + ( t a n θ ) 2 = ρ l H d x (3.14)

解此方程得悬链线的一般方程为

y = ρ l H c h ( ρ l H x + c 1 ) + c 2 (3.15)

由初始条件:猫链拖地,并设拖地点为原点,则对于拖地点有 x = 0 , y = 0 , tan θ = 0 代入悬链线的一般方程得 c 1 = 0 , c 2 = ρ l H ,则猫链的方程

y = ρ l H c h ( ρ l H x + 1 ) (3.16)

根据题中给出的参考图,我们不妨把他们看成在一个平面内,以猫为坐标原点,海底为x轴,与海底垂直向上为y轴,我们假设浮标吃水深度为h,第i根钢管的倾斜角为 θ i ( i = 1 , 2 , 3 , 4 ) θ 4 为钢桶的倾斜角, θ 5 为猫链与钢桶连接处悬链线的切线的夹角, f i + 1 i ( i = 0 , 1 , 2 , 3 , 4 , 5 ) 为第 i + 1 根钢管对i根钢管的拉力,为了表达方便把浮标、钢桶、猫也看成钢管, i = 0 对应浮标。当系泊系统处于平衡状态是,水平方向与垂直方向的合力为零,我们把浮标、钢管、钢桶、重力球、猫链、猫所受力分解为水平方向与垂直方向得到平衡方程,求解方程得到系泊系统各个参数,得到问题的解。

1) 浮标受力情况

浮标受到的浮力为: f = 2 R h ρ g

浮标所受海风的应力为: f = 0.625 v 2 2 R ( a h )

浮标自身的重力为: m 1 g

与钢管连接处钢管的拉力 f 10 = ( f 10 x , f 10 y )

系泊系统处于平衡状态时,水平与垂直方向上合力为零。

即得平衡方程 { f m 1 g f 10 y = 0 f f 10 x = 0

2) 钢管受力情况

当系泊系统处于平衡状态时,根据物理力学知识,每个节点上作用力与反作用力大小相等,方向相反,每个节点上力矩总和为零。

第一节钢管受力分析,把第一节钢管与浮标连接处看成相对转动点,浮标对第一节钢管的拉力 f 10 = ( f 10 x , f 10 y ) ;钢管自身的重力为 m 2 g ;钢管的浮力为 ρ V g ;第二节钢管对第一节钢管的拉力 f 2 , 1 = ( f 2 , 1 x , f 2 , 1 y ) ;则这点受力方程式及力矩平衡方程

{ f 10 x f 21 x = 0 f 10 y f 21 y m 2 g + ρ V g = 0 f 21 sin ( θ 2 θ 1 ) + ( ρ V g m 2 g ) sin θ 1 = 0

第二节钢管受力情况类似

{ f 21 x f 32 x = 0 f 21 y f 32 y m 2 g + ρ V g = 0 f 32 sin ( θ 3 θ 2 ) + ( ρ V g m 2 g ) sin θ 2 = 0

第三节钢管受力情况类似

{ f 32 x f 43 x = 0 f 32 y f 43 y m 2 g + ρ V g = 0 f 43 sin ( θ 4 θ 3 ) + ( ρ V g m 2 g ) sin θ 3 = 0

第四节钢管受力情况类似

{ f 43 x f 54 x = 0 f 43 y f 54 y m 2 g + ρ V g = 0 f 54 sin ( θ 5 θ 4 ) + ( ρ V g m 2 g ) sin θ 4 = 0

钢桶受力情况类似

{ f 54 x f 65 x = 0 f 54 y f 65 y m 3 g m 4 g + ρ V g = 0

浮标的吃水深度为h,则当系泊系统处于平衡状态时,系泊系统位于水下部分的深度为

L = h + l ( cos θ 1 + cos θ 2 + cos θ 3 + cos θ 4 ) + l t cos θ 5 + s y (3.17)

其中 s y 锚链在水中的深度;当系泊系统处于平衡状态时,根据前面受力分析,利用每个节点之间作用力与反作用力之间的关系,从浮标出发考虑问题,当浮标的吃水深度为h等参数确定时,利用力学原理及平衡状态下课计算确定上述参数 θ i ( i = 1 , , 5 ) ;因此我们下面只要确定锚链在水中的深度 s y 即可,而猫链在水中

的方程是 y = ρ l H c h ( ρ l H x + 1 ) ,只要确定猫链与钢桶的连接点坐标, s y 的大小即可确定。因此只要确定海水

的深度,根据我们给出的模型,就可计算、模拟出系泊系统各项部件的参数要求,从而设计一套系泊系统。

本文着重分析了浮标、钢管、钢桶、重物球、猫链之间的受力关系及相互之间的影响,建立了它们之间受力方程。建立了猫链的曲线方程,得到只要确定海水深度,就可确定各部件在给定的条件下的参数,但各根钢管、猫链的张力分配问题没有考虑,需后续进一步讨论,其特点在工程设计中,类似的问题根据给出的设备情况可以参考该模型,其优势具有一定的通用参考价值,可以引用模型分析具体问题。特别是对每个部件的受力分析及部件之间的衔接分析对类似问题可引用。

4. 总结与展望

本文对系泊系统组成的各个部件参数进行了探究,通过悬链线理论 [11] 对部件及部件之间的受力进行分析,获得:1) 各个部件级部件之间的关系方程表达式,2) 部件之间受力方程,3) 锚链的方程表达式,通过方程优化在一定条件下得到这类系泊系统的各个物理参数。下一步需对这些参数及部件需要进一步优化、修改,同时要综合考虑风力、水流力等复杂情况下系泊系统设计,本文从数学角度进行分析,与实际情况有差异,望后续继续分析改进。

参考文献

参考文献

[1] 王丹, 刘家新. 一般状态下悬链线方程的应用[J]. 船海工程, 2007, 36(3): 26-28.
[2] 袁梦, 范菊, 朱仁传, 缪国平. 基于悬链线理论的系泊系统势能[J]. 上海交通大学学报, 2011, 45(4): 597-603.
[3] 康健, 陈国明, 朱本瑞. 四锚链工程船舶锚泊优化分析[J]. 石油矿场机械, 2013, 42(2): 63-66.
[4] 范会渠. 弹性悬链线方程参数变换法及其工程应用[J]. 力学与实践, 2010, 32(2): 32-35.
[5] 郭迎雪, 姚晴. 基于悬链线方程的系泊系统状态分析[J]. 应用数学进展, 2017, 6(3): 296-307.
[6] 刘涌国. 带悬锤的悬链线方程求解及其简化方法[J]. 水运工程, 2014(3): 103-106.
[7] 刘阳, 李杨, 闫琛. 悬链锚桩系泊系统锚桩结构设计分析[J]. 海洋工程装备与技术, 2019, 6(1): 199.
[8] 何英豪, 朱雨薇, 韩冰凌. 基于悬链线方程的系泊系统分析[J]. 工业设计, 2016(11): 136-137, 139.
[9] 王希春, 王静. 利用自悬链模型研究系泊系统的设计问题[J]. 数学学习与研究, 2018(13): 139.
[10] 孙宁松. 海上移动式平台锚泊定位系统锚索链受力分析[J]. 中国海洋平台, 1997, 23(2): 41-44.
[11] 张继明, 范秀涛. 海洋资料浮标锚泊系统的系泊力计算[J]. 山东科学, 2014, 27(2): 19-24.
[12] 王领, 刘生法. 深水大型浮台系泊系统水动力计算分析[J]. 广东造船, 2018(5): 16-18.
[13] 洪文渊. 绷紧式深海系泊缆与上部浮体的耦合动力分析[D]: [硕士学位论文]. 杭州: 中国计量学院, 2015.
[14] 潘斌, 高婕. 浮标系泊系统静力计算[J]. 重庆交通学院学报, 1997, 16(1): 68-73.
[15] 谭兵, 伍悦. 基于迭代搜索算法的系泊系统的优化设计[J]. 兰州文理学院学报, 2017, 31(5): 31-35, 45.
[16] 刘昌欣, 范兴奎. 基于静力学分析的系泊系统状态浅析[J]. 应用数学进展, 2017, 6(3): 357-366.
[17] 王宏伟. 深海系泊系统模型截断技术研究[D]: [博士学位论文]. 哈尔滨: 哈尔滨工程大学, 2011.
[18] 郑中强. 海洋系泊系统的非线性动力学研究[D]: [博士学位论文]. 青岛: 中国海洋大学, 2015.
[19] 张欣欣, 郭小天, 张亮. 非线弹性索对系泊系统性能的影响[J]. 应用科技, 2014(1): 84-88.
[20] 张隆辉. 系泊系统主动式截断混合模型试验技术仿真研究[D]: [硕士学位论文]. 北京: 中国舰船研究院, 2016.
[21] 王贝贝, 于蓉, 马灵灵, 周克元. 系泊系统的最优设计[J]. 保山学院学报, 2017(5): 41-43.
[22] 李斐, 吕文, 李琴. 水流影响下系泊系统设计的数学模型[J]. 烟台大学学报, 2017, 30(3): 173-179.
[23] 成海亮. 船舶T型系泊系统安全性研究[J]. 中国水运(下半月), 2017, 17(6): 22-23.
[24] 张路莹, 刘绍庆. 基于多目标规划的系泊系统的设计[J]. 曲阜师范大学学报, 2017, 43(3): 49-53.
[25] 吴婷婷. 开敞式码头系泊系统状态预测方法研究[D]: [硕士学位论文]. 大连: 大连交通大学, 2016.
[26] 徐剑峰. 浅水斜底地形下超大型浮体单模块系泊系统研究[D]: [硕士学位论文]. 上海: 上海交通大学, 2018.
[27] 尘颖, 许佰雁. 基于迭代法的系泊系统问题模型分析[J]. 洛阳师范学院学报, 2017, 36(2): 19-22.
[28] 袁培银, 赵宇. 多浮体式系泊系统设计及水动力性能分析[J]. 舰船科学技术, 2016, 38(2): 80-84.