船只补给悬链线系统运动特性分析

吴卫国 ,谭振东,王振涛,冷永刚,杜尊峰

(1.天津大学 机械学院, 天津 300072； 2.军事交通运输研究所, 天津 300161； 3.天津市科委, 天津 300051 ； 4.天津大学 水利工程仿真与安全国家重点实验室,天津 300072)

船只补给悬链线系统运动特性分析

吴卫国1,2,谭振东3,王振涛2,冷永刚1,杜尊峰4

(1.天津大学 机械学院, 天津 300072； 2.军事交通运输研究所, 天津 300161； 3.天津市科委, 天津 300051 ； 4.天津大学 水利工程仿真与安全国家重点实验室,天津 300072)

船只索道补给越来越广泛用于两船或船岸之间物资传送，本文综合考虑船只与锚泊系统的耦合以及船只运动对索道的影响，通过计算波浪中船只补给悬链线系统的动力学响应，得到了船只补给悬链线系统的运动特性，为系统设计提供依据。运用索道悬链线理论和有限元方法，对模型进行了计算，并对数值模拟结果和理论结果进行了比较。结果表明：数值模拟结果和理论结果整体吻合较好；三级海况下，连接设备时应将船只侧面对着固定物，连接点应在船只的中心线上。高海况条件下，应考虑垂荡的影响，理论模型简化因素较多，索道跨距80 m之后，应以仿真分析模型为主。

船只补给；悬链线系统；运动特性；动力学响应；有限元方法

随着海洋开采技术的发展，许多大的石油或天然气平台逐步出现，使得海上大型浮体结构的设计研究成为热点问题，其中浮体和浮体之间的连接技术又是一个关键性问题。在此方面，许多研究者研究了相连浮体的水作用响应[1-2]。Fang等[3]采用二维切片法计算了两平行布置的条状浮体间的水动力相互作用。Duncan等[4]采用Oortmerssen[5]的方法计算了两船体的运动响应，并在计算中考虑了锚泊缆和防碰垫的作用。Chen and Mahrenholtz用边界元的方法，讨论了两漂浮圆柱之间的连接刚度、间距及水深等参数的变化对圆柱运动响应的影响[6]。Newman[7]采用边界积分的方法，对两个相同的均质方箱在固接和铰接情况下的运动响应进行了数值计算。沈庆等[8]对于铰联结系泊多浮体系统，将其作为铰联结无根树系统，用多刚体力学方法进行了运动学分析。勾莹等[9]应用边界积分方程方法研究了波浪与2个相连三维浮体的相互作用问题，计算中考虑了浮体相互间的水动力干涉, 通过联立各浮体的运动方程以及运动连续条件求得浮体运动量。Hong等[10]通过数值计算研究系泊多浮体运动响应特性。Sun等[11]基于线性绕射理论，并采用拉格朗日乘子方法，探讨了两浮体在铰接、固接、弹性连接情况下的运动响应以及不同连接形式下连接点的受力情况。刘文玺等[12]在浮体和系泊系统的时域耦合分析中，将浮体的分析方法和系泊缆索的分析方法结合，建立了浮体和系泊系统非线性运动响应和缆索张力的耦合计算方法。许鑫等[13]对目前海洋工程中较为常见的多浮体系统进行了介绍，并就国内外在多浮体系统的水动力方面的研究进行了系统的阐述。

目前，浮体与浮体在水作用力下相互作用的研究已达到一定深度，但是关于浮体和固定物体连接作用的研究还不多见，尤其是船只与固定物体连接的动力学响应研究则更是鲜见报道。

1 船只补给悬链线系统理论分析

船只补给悬链线系统是由船只、悬链线与固定物体组成。悬链线连接一般为索道连接，本文中悬链线统一称为索道。船只补给悬链线连接主要用于海上浮体与相对静止设施之间的物资和人员传送，它区别于浮体与浮体连接，表现为连接的功用不同，连接的方式不同，连接的距离也有差别。此外，船只补给悬链线系统的连接方式也不同于山地丛林的索道连接方式[14]，山地丛林索道主要是在固定的物体上展开。舰船补给索道连接[15]在舰船航行补给时，两船处于航行状态，保持适当的航速, 并利用操舵装置产生横向力, 抵消索道张力, 维持一定的船距。

船只对固定目标补给时要求船只在锚泊状态下仅通过索道连接固定目标，索道只能依靠补偿装置自身来控制。具体连接工作方式是：带有补给品的船只当接近固定目标时，应在开阔的水域内锚泊确定其位置，如图1中的浮体。而后人员在船只上抛送绳索一端到固定物目标上，将绳索与固定物的某处位置连接，然后绳索的另一端与船只连接，并通过恒张力补偿装置将绳索收紧，实现船只和固定物的索道连接。船只水线以上部分会受到风载作用，而船只水线以下部分会受到海流作用，船只受到的这两种作用载荷由船只抛锚的锚链力进行平衡，使船只在平衡位置附近漂浮运动。船只在水中漂浮运动时，索道会出现拉紧或松弛，系统中恒张力补偿装置自动收放索道使得索道处于张紧状态，从而保持恒定张力。这里恒张力补偿装置如何得到索道跨距以及补偿量等参数并严格控制索道张力，相关研究可见文献[16]，本文只考虑船与固定物之间相对运动时(靠近或者离开)缆的张力基本保持恒定，以确保物资能够在船只和固定物之间保持平稳运送的情况。

图1 船只和固定物的连接方式示意图Fig.1 The connecting way of vessel and fixation object

船只和固定物之间有一定的距离空间，二者通过悬链线相互作用。船只和固定物的距离变化与两个船只或浮体间的距离变化是不同的，两个船只间或浮体间的距离是通过航行或运动来保持。而固定物和船只的距离是通过船只抛锚固定的受力平衡方式保持。因此有必要研究船只锚泊情况下的波浪运动响应，为船只补给悬链线系统分析奠定基础。

1.1 船只运动分析

1.1.1 船只运动平衡位置及稳定性分析

船只在波浪中受到风、浪、流引起的力以及力矩的作用如图2所示。

图2 船只受力分析示意图Fig.2 The analysis of vessel forces

作用在船只的纵向总载荷FXT、横向总载荷FYT及总摇艏力矩MXYT分别为

(1)

式中：FXW是纵向风力，FYW是横向风力，FXC是纵向流力，FYC是横向流力，XD是纵向二阶波浪漂移力，YD是横向二阶波浪漂移力，MXYC是流的摇艏力矩，MXYW是风的摇艏力矩，ND是漂移力的摇艏力矩。

对于单点系泊船只，平衡位置的计算方程为

(2)

式中：∑M为合力矩，ARM为船只中心和船艏系索点处的距离(一般取ARM=0.48LOA)，LOA为总船长。

1.1.2 船只三维运动响应

船只在水中的三维运动响应主要是受到不规则波的作用，这种不规则波可被分解成为多个规则子波进行处理，通过三维格林函数法计算规则波中船只的运动响应，然后进行叠加就可以得到船只总体三维运动响应。在不规则波的多规则波分解中，频率为n规则波能够被获得，命名为ωi，i=1,2，…,n。

(3)

(4)

式中：MS为质量矩阵，A为附加质量矩阵，B为阻尼系数矩阵，KW为恢复力系数矩阵，Fe为激励力，CS为粘性阻尼，从试验或经验公式中获得，KL是缆索的刚度矩阵。通过三维格林函数计算船只在规则子波的运动，然后再进行叠加，进而得到船只在不规则波的运动响应。

1.1.3 船只与锚泊链耦合动力分析

在对船只与锚泊链耦合动力分析中，船只的水动力学特性采用三维势流理论进行求解。船只的拖拽力采用莫里森公式进行计算分析，锚泊线的动力分析采用杆理论和有限元模型，在总体坐标系中建立控制方程进行计算，锚泊线和补给船舶之间的连接采用平动、转动弹簧以及阻尼器进行模拟，通过调整弹簧和阻尼器的控制参数可以模拟各种形式的连接，船只和锚泊线耦合系统的方程，如下所示

(5)

式中:KL是锚泊线和弹簧(阻尼)连接器的联合刚度矩阵，KB是刚度矩阵，KC和KCT分别是包含船只和锚泊线耦合项的耦合刚度矩阵及其转换矩阵，UL和UB分别是锚泊线和船只的位移列阵，FL和FB分别是作用在锚泊线和船只上的力列阵和力矩列阵。

1.2 悬链线(索道)计算分析

索道受力的计算理论主要有悬链线理论和抛物线理论，研究高架索道特性的方法主要有近似法和精确法；抛物线理论假定垂向载荷沿水平方向均匀分布，悬链线理论假定垂向载荷沿着索道长度弧度方向均匀分布。当索道只受到重力时，线的真实形状应该是悬链线[17-18]。根据文献[18]可得到索道悬链线的力学分析模型。

1.2.1 无载荷索道计算模型

钢索在跨度内沿长度方向均匀分布，以钢索最低点为原点建立直角坐标系，见图3。

图3 钢索曲线受力分析图Fig.3 The analysis of ropeway forces

悬链线的方程是一个双曲余弦函数，其方程形式为

(6)

式中：R为无载荷索道钢索受到的水平分力，w是悬链线即钢索的单位长度重量，悬链线的斜率角θ。

(7)

钢索各点受到的拉力T为

(8)

1.2.2 考虑集中载荷的作用

集中载荷是通过货物的重力引起的，货物重力Q作用在钢索时，其受力可表示为

(9)

(10)

式中：R为有荷载索道受到的水平拉力，Ta为有荷载索道AD段的平均拉力，Tb为有荷载索道DB段的平均拉力，k为AD的距离与DB距离的比值；l0为A、B支点的水平跨距，TB为有荷载索道AD段最低点的水平坐标，TL为有荷载索道DB段最低点的水平坐标，TR为有荷载索道AD段的悬链线长度，L2为有荷载索道DB段的悬链线长度。

1.3 综合运动分析

船只在波浪中的运动受到风、浪、流等影响较大，运动响应也与船只的类型有关，在低海况条件下，船只偏离平衡位置的距离通过索道的张紧装置得到补充，索道连接点的位置也变化较小，因此，低海况条件下，计算其综合运动时，可考虑索道长度的变化，忽略船只运动对索道的其他影响。连接索道受力分析模型如图4所示，根据文献[19]可以得到货物的运动方程。

图4 船只补给悬链线系统的计算模型Fig.4 The calculation model of vessel supply catenary system

货物沿索道运动的轨迹为[19]

(11)

式中：q为钢索的单位长度质量，h为A、B两点的高度差，Q为货物(包括吊钩等)重力，β为AB与水平方向的夹角，H为钢索的水平方向的分力，L为索道跨距，C0为常数。

货物运动方程为[19]

(12)

式中：m为货物(包括吊钩等)质量；c为阻尼系数，c=2ζω0m；ζ为无量纲阻尼比；αL和αR为如图4所示的PA、PB与水平方向的夹角，L0为A点和Q点的水平距离。

当x为L/2时，挠度最大，从式(11)可得最大挠度fmax与水平力H的关系为

(13)

由分析知：随着索道跨距的增加，最大挠度也随着增大，货物有可能落水；高度差h越大，最大挠度减小，最低点纵坐标值越小，即离水面越远；索道水平力与最大挠度成反比关系，索道张力越大，最大挠度越小，但是对索道的要求越高。因此，必须控制索道的跨距大小，保持两端一定的高度差，确保张力在安全范围内，这样才能安全运输物资。下面将对船只及悬链线系统的运动进行仿真分析。

2 数值仿真模拟

船只与固定物的连接，可以以运输船和海上石油平台为例。大型石油平台相对较大，可考虑为固定物。考虑到运输船和平台连接的安全性，连接设备在三级海况下展开，运输船的型深为4.7 m，非满载下，轻载状态时的高度为4.3 m，船舶吃水为1.955 m，悬挂点的高度为7.05 m，坐标系的原点建立在吃水线平面，船舶被分成20个部分，在每个部位均用面单元建立网格，运输船的模型见图5。三级海况下的环境载荷见表1，通过仿真模拟得到三级海况下运输船的运动位移、速度和加速度，结果见表2～4。三级海况下的运输船耦合运动分析结果见图6。

图5 运输船网格Fig.5 The grid of carrier

三级海况的环境载荷见表1。

表1 三级海况环境载荷

表2三级海况下的运动位移

表3 三级海况下的速度

表4 三级海况下的加速度

(a)运输船的纵荡

(b)运输船的横荡

(c)运输船的垂荡

(d)运输船的横摇

(e)运输船的纵摇

(f)运输船的艏摇图6 三级海况下的运输船耦合运动分析结果Fig.6 The analysis result of coupling motion at the level 3 sea states

由分析可知，三级海况下，对船只和固定物连接影响较大的因素为纵荡和纵摇，纵荡的最大值可达到1.2 m，纵摇的速度达到2.5 m/s，纵摇加速度达到4.4 m/s2，这对索道连接是极其不利的。因此，在连接设备时应该使船只侧面对着固定物，连接点应该在船只的中心线上，降低纵荡和纵摇对索道的影响。垂荡会引起连接点位置的变化，在三级海况下垂荡的最大值为0.15 m，对连接点位置影响较小，可忽略其影响，在高海况条件下，应考虑该因素。

在考虑了锚泊系统耦合运动后，船只的运动较为复杂，在做周期运动的同时，也会随着锚链的影响做小范围运动，运动的幅度较小，对于索道的连接影响较小，但要注意锚的位置变化，高海况下，会发生走锚的现象，对于船只与固定物的连接较为不利。

连接绳索的跨距分别取40、60、80 m，三个索道的参数为：16、12、12 mm，绳索的长度为100、200、100 m。16 mm的索道单位质量为0.99 kg/m，其他类型的绳索单位质量为0.6 kg/m，货物的重量为50 kg，通过理论计算得到不同跨距下拉力及其挠度，见表5。

表5 不同跨距下索道拉力及其挠度

Table 5 The ropeway tension and deflection in the different ropeway spans

索道跨距最小值状态实际工作状态最小拉力/N挠度/m挠度/m水平拉力/N最大拉力/N索道长度/m801725.524.472488.82840.881.6601404.217.871410.51826.062.13401085.511.347639.71149.1243.1

从表5来看，张力变化幅度随着挠度增大而减小。在相同的挠度时，张力随着距离的增大而减小。

下面在ansys中模拟索道的连接，索道采用Beam 188单元建模，等效截面为9.5 mm，钢索密度为7.85 t/m3，重力加速度为9.8 m/s2，两端采用绞支，如图7所示。仿真时接收索必须同放索保持相同速度，考虑到货物的重力为490 N，ANSYS仿真结果在表6中，ANSYS仿真结果与悬链线计算的结果比较见表7。

图7 索道模拟单元划分示意Fig.7 The simulation of ropeway

有限元仿真结果在表6中，有限元仿真结果与悬链线计算的结果比较见表7。

表6 索道的有限元计算结果

表7 索道的理论结果和仿真结果对比

Table 7 The comparison between theoretical and simulation results

索道跨距7/mANSYS结果理论结果差值拉力/N索道重量/N拉力/N索道重量/N拉力/N索道重量/N802868.981778.72840.81759.2928.1819.41601842.51135418261339.5216.5114.48401152.21939.31149.12929.233.0910.07

图8 理论与仿真得到索道拉力差值随跨距变化Fig.8 The ropeway tension difference between theoretical results and simulation results in different ropeway spans

图9 牵引拉力随离岸距离变化关系(索道跨距80 m)Fig.9 The ropeway tension in different ropeway spans. (The ropeway span is 80 m)

从表7和图8来看，仿真结果和理论结果整体吻合，随着跨距的增加，拉力差值在增大，在跨距为80 m时，拉力差值为28 N，与拉力的比值约为1%，理论分析值相对仿真结果值偏小，表明理论模型简化因素较多，随着跨距增加，索道跨距80 m之后，应以仿真分析模型为主。

从图9来看，在索道跨距80 m时，在船上牵引岸上的索具过程中，牵引拉力随距岸距离是变化的，牵引拉力在40 m的范围内变化不大，超越这一范围后，随着距离的增加呈线性增加，也会越来越费力。

3 结论

本文主要分析了船只补给悬链线系统的运动特性，并考虑了船只与锚泊系统的耦合以及船只运动对索道的影响，通过理论分析和仿真模拟，得到了船只悬链线系统连接动力学响应，并对结果进行了分析，得到以下结论：

1)通过对船只与悬链线运动分析，得到了系统连接的运动特性，为船只悬链线系统设计提供理论依据。

2)三级海况下，对船只和固定物连接影响较大的因素为纵荡和纵摇，在连接设备时应将船只侧面对着固定物，连接点应在船只的中心线上，降低纵荡，纵摇对索道的影响。

3)垂荡会引起连接点位置的变化，在三级海况下垂荡的最大值为0.15 m，对连接点位置影响较小，可忽略其影响，在高海况条件下，应考虑该因素。

4)考虑了锚泊系统耦合运动后，船只的运动较为复杂，在做周期运动的同时，也会做小范围运动，对船只和固定物连接影响较小，但应注意防止走锚现象的发生。

5)仿真结果和理论结果基本保持一致，随着跨距的增加，拉力差值在增大，跨距小于80 m时，拉力差值与拉力的比值为1%；理论分析值相对仿真结果值偏小，表明理论模型简化因素较多，随着跨距增加，索道跨距80 m之后，应以仿真分析模型为主。

6)在索道跨距80 m时，在船上牵引岸上的索具过程中，牵引拉力随距岸距离是变化的，牵引拉力在40 m的范围内变化不大，超越这一范围后，随着距离的增加呈线性增加，也会越来越费力。

[1]CHEN S, MAHRENHOLTZ O. Interaction of water waves and floating twin cylinders in beam waves[J]. Applied ocean research, 1994, 14(6): 371-379.

[2]MCTAGGART K, CUMMING D, HSIUNG C C, et al. Seakeeping of two ships in close proximity[J]. Ocean engineering, 2003, 30(8): 1051-1063.

[3]FANG M C, KIM C H. An analysis of water shipping between two floating platforms in the beam wave[C]//Proceedings of the Fifth International Symposium on Offshore Mechanics and Arctic Engineering. Tokyo, 1986, 1: 499-505.

[4]DUNCAN J H, BARR R A, LIU Y Z. Computations of the coupled response of two bodies in a seaway[C]//International Workshop on Ship and Platform Motions. Berkeley, CA, 1983.

[5]VAN OORTMERSSEN G. Hydrodynamic interaction between two structures floating in waves[C]//Proceedings of 2nd International Conference on Behaviour of Off-shore Structures. London, 1979: 339-356.

[6]CHEN S, MAHRENHOLTZ O. Interaction of water waves and floating twin cylinders in beam waves[J]. Applied ocean research, 1992, 14(6): 371-379.

[7]NEWMAN J N. Wave effects on deformable bodies[J]. Applied ocean research, 1994, 16(1): 47-59.

[8]沈庆, 陈徐均. 系泊多浮体系统波浪运动响应的动力学分析[J]. 解放军理工大学学报, 2002, 1(4): 31-36. SHEN Qing, CHEN Xujun. Dynamic analysis of moored multi-body floating system’s responses to regular waves[J]. Journal of PLA university of science and technology, 2002, 1(4): 31-36.

[9]勾莹, 滕斌, 宁德志. 波浪与两相连浮体的相互作用[J]. 中国工程科学, 2004, 6(7): 75-80, 93. GOU Ying, TENG Bin, NING Dezhi. Interaction effects between wave and two connected floating bodies[J]. Engineering science, 2004, 6(7): 75-80, 93.

[10]HONG S Y, KIM J H, CHO S K, et al. Numerical and experimental study on hydrodynamic interaction of side-by-side moored multiple vessels[J]. Ocean engineering, 2005, 32(7): 783-801.

[11]SUN L, TAYLOR R E, CHOO Y S. Responses of interconnected floating bodies[J]. The IES journal part a: civil & structural engineering, 2011: 4(3): 143-156.

[12]刘文玺, 周其斗, 张纬康, 等. 浮体与系泊系统的耦合动力分析[J]. 船舶力学, 2014, 18(8): 940-957. LIU Wenxi, ZHOU Qidou, ZHANG Weikang, et al. Dynamic analysis of the coupled floating body/mooring system[J]. Journal of ship mechanics, 2014, 18(8): 940-957.

[13]许鑫, 杨建民, 李欣, 等. 海洋工程中多浮体系统的水动力研究综述[J]. 中国海洋平台, 2014, 29(4): 1-8, 13. XU Xin, YANG Jianmin, LI Xin, et al. Review of the research on multi-body system hydrodynamics[J]. China offshore platform, 2014, 29(4): 1-8, 13.

[14]张立忠. 我国矿山载人索道现状及发展前景[J]. 中国索道, 2002, 2(5): 1-4. ZHANG Lizhong. Development and application of carrying people ropeway in the mine[J]. China ropeway, 2002, 2(5): 1-4.

[15]任爱娣, 张得伟, 何学军. 航行横向补给舰船运动对高架索张力影响[J]. 舰船科学技术, 2014, 36(7): 94-96, 102. REN Aidi, ZHANG Dewei, HE Xuejun. Influence of tension in highline cable of alongside replenishment by the motion of ships[J]. Ship science and technology, 2014, 36(7): 94-96, 102.

[16]余建星, 顾鹏. 海上干货补给技术[J]. 海洋技术, 2005, 24(3): 105-110. YU Jianxing, GU Peng. The technique of solid cargo supplying on the sea[J]. Ocean technology, 2005, 24(3): 105-110.

[17]郑丽凤, 周新年, 巫志龙. 悬索的理论计算与实测误差分析[J]. 北华大学学报: 自然科学版, 2010, 11(2): 162-168. ZHENG Lifeng, ZHOU Xinnian, WU Zhilong. Error analysis among theoretical calculation and measure data of flexible cable[J]. Journal of Beihua university: natural science, 2010, 11(2): 162-168.

[18]周新年. 架空索道理论与实践[M]. 北京: 中国林业出版社, 1996: 140-200. ZHOU Xinnian. Theory and practice of aerial ropeway[M]. Beijing: China Forestry Publishing House, 1996: 140-200.

[19]余建星, 顾鹏. 海上分离式轻型补给系统研究[J]. 天津大学学报, 2006, 39(11): 1270-1271. YU Jianxing, GU Peng. Offshore minitype detached replenishment system[J]. Journal of Tianjin university, 2006, 39(11): 1270-1271.

Analysis of motion characteristics of supply vessel catenary systems

WU Weiguo1,2，TAN Zhendong3，WANG Zhentao2，LENG Yonggang1，DU Zunfeng4

(1.School of Mechanical Engineering, Tianjin University, Tianjin 300072, China; 2.Military Transportation Institute, Tianjin 300161, China; 3.Tianjin Science and Technology Commission, Tianjin 300051, China; 4.State Key Laboratory of Hydraulic Engineering Simulation and Safety, Tianjin University, Tianjin 300072, China)

The cableway has been widely used by supply vessels in shore-to-ship or ship-to-ship transfers of goods and materials. By analyzing the motion of the catenary system in supply vessels, and considering the effects of the coupling of vessels and mooring systems and the shipmotion on cableways, we calculated the dynamic response of the catenary systems of supply vessels, from which obtained the motion characteristics of vessel catenary systems. The result provides the basis for the cableway design. Based on the cableway catenary theory and using the finite element method (FEM), we established a model and performed simulations. We then compared the numerical simulation and theoretical results, which showed that the results were mostly consistent. In a sea state code of 3, vessels should be connected to the device side facing the fixture and the connection point should be on the centerline of the vessel. In a very rough ocean environment, the impact of heave should be considered. As there are many simplified factors in the theoretical model, when the cableway span exceeds 80 m, the simulation model should be considered first.

vessel supply; catenary system; motion characteristics; dynamic response; FEM

2015-09-14.

时间：2016-12-21.

国家自然科学基金项目(51109158)；天津市应用基础与前沿技术研究计划(15JCZDJC32200).

吴卫国(1983-), 男, 博士； 冷永刚(1964-), 男, 教授，博士生导师.

冷永刚，E-mail：leng_yg@tju.edu.cn.

10.11990/jheu.201509045

U644

A

1006-7043(2017)01-126-07

吴卫国,谭振东,王振涛,等. 船只补给悬链线系统运动特性分析[J]. 哈尔滨工程大学学报， 2017, 38(1): 126-132. WU Weiguo，TAN Zhendong，WANG Zhentao,et al. Analysis of motion characteristics of supply vessel catenary systems[J]. Journal of Harbin Engineering University, 2017, 38(1): 126-132.

网络出版地址：http://www.cnki.net/kcms/detail/23.1390.u.20161221.1524.002.html