舰船科学技术  2019, Vol. 41 Issue (3): 8-12   PDF    
系泊多浮体耦合系统水动力特性数值分析
李鸿1, 张斐斐1, 李莉2, 秦超1     
1. 哈尔滨工程大学 航天与建筑工程学院,黑龙江 哈尔滨 150001;
2. 哈尔滨工程大学 材料科学与化学工程学院,黑龙江 哈尔滨 150001
摘要: 基于系泊浮筒链防撞系统,利用AQWA软件对不同海洋环境下的系泊多浮体耦合系统进行全耦合时域分析,得到浮体的运动响应和缆绳受力情况。结果表明,系泊缆绳最大张力出现在端部系泊缆上,中部系泊缆的迎浪侧缆绳比背浪侧缆绳张力大;单纯波浪环境下,浮体平均横向位移较小,在初始位置往复运动,波流共同作用时浮体发生前后摆动,在水流前进方向运动量较大,而且不能回到初始位置。最后将数值模拟结果与现有试验结果进行对比,验证了该全耦合时域分析方法的可行性,能够为系泊多浮体耦合系统的水动力特性研究提供参考依据。
关键词: 多浮体     运动响应     动力响应     全耦合时域分析     AQWA    
Numerical analysis of hydrodynamic characteristics of mooring multi-body coupled system
LI Hong1, ZHANG Fei-fei1, LI Li2, QIN Chao1     
1. Harbin Engineering University, College of Aerospace and Civil Engineering, Harbin 150001, China;
2. Harbin Engineering University, College of Materials Science and Chemical Engineering, Harbin 150001, China
Abstract: Based on the mooring buoy chain anti-ship collision system, AQWA software is used to conduct a fully coupled time-domain analysis of the mooring multi-body coupled system in different marine environments, and the motion response of the floating body and the force of the cable are obtained. The results show that the maximum tension of the mooring cable appears on the end mooring cable. The front side cable of the middle mooring cable has greater tension than the backside cable. In wave environment, the average lateral position of the floating body is relatively small and it is reciprocating in the initial position. The buoyant body oscillates back and forth when the wave-current interact, and the amount of movement in the forward direction of the water flow is large, and it cannot return to the initial position. Finally, the numerical simulation results are compared with the existing experimental results to verify the feasibility of the fully coupled time-domain analysis method and provide reference for the study of the hydrodynamic characteristics of the mooring multi-body coupled system.
Key words: floating multi-body     motion responses     dynamic responses     fully coupled time-domain analysis     AQWA    
0 引 言

各类码头的建成、海上运输的高速发展和我国领海海域的保护,使系泊多浮体系统得到广泛应用。针对桥梁的非通航孔设计了各类浮式防船撞设施,避免因桥墩受到船舶撞击而带来严重的人员伤亡和经济损失。在各类危险品码头也逐渐发展了浮体系泊式防撞系统,用于拦截由于恶劣环境或船舶机械故障引起的意外事故,保护码头的同时也避免了灾难性的海上污染。

沈庆[12]采用势流理论计算了规则波作用下的单系泊浮体受到的流体压力,使用Huston方法和Kane方程提出了用于求解多浮体系统中浮体质心位移和角位移的复数线性代数方程组。王翔[3]采用三维频域函数法得出波频和浪向对浮体的运动响应和锚缆张力产生很大的影响,且在低频范围内出现极值。王桂波[4]通过模态法研究了受到规则波作用的铰接多浮体结构频域运动响应,发现浮体的运动响应受到水深、铰接位置和波浪周期的影响。袁培银[5]对深海半潜式多浮体新型系泊系统进行了时域分析,新型设计可以在一定程度上改善平台的运动响应,满足系泊缆的安全要求,在钻井作业时保证平台的安全性。高峰[6]采用物理模型研究了多船过驳平台的系泊特性,试验得到各浮体在风浪流作用下的运动量,并测试了系泊缆张力的变化。于文太[7]使用Matlab-GUI开发的计算软件,分析了多浮筒悬链线系泊缆索姿态曲线的影响因素。王贝壳[8]提出了适用于非通航孔桥墩的防船撞拦截网装置,对比AQWA数值模拟结果与试验结果,验证了拦截网工程设计的有效性。夏雪[9]提出了小水线面浮式拦截系泊体系,结合物理模型试验和数值模拟,验证了该系统在深水海域具有更好的适应性。张周康[1011]使用计算流体力学软件和多刚体动力学理论,通过嵌套网格方法、VOF方法对斜浪作用下的浮体系统进行了运动响应的模拟分析,对比实验结果验证了数值模拟方法的正确性。在海洋动力环境要素作用下,王鑫[12]采用物理模型试验方法,对系泊浮筒的运动响应和系泊缆的动力响应进行研究。

系泊浮筒链防撞系统为多浮体耦合结构,涉及多浮体运动响应和和系泊缆绳张力问题,具有对各种海上建筑适应性强、成本低廉、易于施工修理等特点,但是由于各浮体之间存在结构和水动力的耦合作用,运动和动力响应与单浮体系统有明显不同,对其水动力特性缺乏系统研究,这在一定程度上限制了该系统的应用。因此,本文使用AQWA软件,在物理模型试验的基础上建立系泊多浮体耦合系统的水动力模型,对不同海洋环境下的多浮体系统,采用全耦合时域分析方法,进行浮体运动响应和系泊缆绳动力响应的数值计算,得到浮体横移量和系泊缆绳张力的数值解,并与物理模型试验的结果进行对比分析。

1 多浮体耦合系统

多个浮体相互连接组成系泊多浮体系统,多浮体系统的形式多种多样,海上浮式移动平台FPSO、桩群结构以及系泊防撞系统等,本文选取系泊浮筒链式防撞系统[12]进行系泊多浮体耦合系统的全耦合时域分析。

1.1 系泊浮筒链系统模型

浮筒直径4 m,长18 m,质量为80 kg,每个浮筒分布有4条系泊缆,浮筒间缆绳长5 m,浮筒链总长度为552 m,水深48 m,浮筒和缆绳的布置以及缆绳编号如图1所示。

图 1 系泊多浮体耦合系统布置图 Fig. 1 Arrangement of mooring multi-body coupled system
1.2 系泊缆绳参数

采用线性缆绳,王鑫[12]在物理模型试验中采用的缆绳受力变形曲线如图2所示,使用Matlab将试验中的缆绳受力变形曲线拟合成下列函数:

图 2 缆绳受力变形曲线 Fig. 2 Force-deformation curve of cable
$ y = 83.8{x^3} + 5\;251{x^2} + 90\;259.8x - 11\;701.8{\text{。}} $ (1)
1.3 系泊浮筒链系统水动力模型的建立

采用AQWA软件,在物理模型试验的基础上建立系泊系统的水动力模型,如图3所示。

图 3 系泊系统水动力模型 Fig. 3 The hydrodynamic model of mooring system
2 海洋环境动力要素 2.1 波浪谱

采用不规则波JONSWAP谱进行数值模拟,其表达式为[13]

$ {S_{JON}}\left( \omega \right) = A{S_{PM}}\left( \omega \right){\gamma ^{\exp \left( { - 0.5{{\left( {\frac{{\omega - {\omega _P}}}{{\sigma {\omega _P}}}} \right)}^2}} \right)}}{\text{。}} $ (2)

式中: $ \gamma $ 为谱峰升高因子; $ \sigma $ 为谱型参数; $ A = 1 - 0.287$ $ \ln \left( \gamma \right)$ 为无因次参数。

JONSWAP谱的谱峰升高因子 $ \gamma $ 取3.3、有义波高 $ {H_s}$ 和谱峰周期 $ {T_p}$ 均采用和物理模型试验波浪谱相同的参数,波浪要素如表1所示。

表 1 波浪要素 Tab.1 Wave parameters
2.2 海流参数

海流速度随时间变化缓慢,因此在工程设计中,为简化起见,海流被当作稳定的流动,垂线流速分布公式如下所示:

$ \frac{u}{{{U_ * }}} = \frac{1}{\kappa }\ln \left( {\frac{y}{{{y_0}}}} \right){\text{。}} $ (3)

式中: $ u$ 为离床面高度为z处的流速; $ {y_0}$ 为流速等于0处的床面以上距离; $ {U_ * }$ 为摩阻流速; $ \kappa $ 为卡门常数[14]。潮流方向和波浪传播方向一致,流速大小与物理模型试验一致。

3 结果与分析

考虑单纯波浪和波流共同作用的不同海洋环境,利用AQWA软件对系泊浮筒链系统采用全耦合时域分析方法,进行浮体运动响应和系泊缆动力响应的数值计算。浮体在波浪载荷作用下,随着波浪方向发生位移响应,当波浪或水流进一步在浮体周围形成绕射水流后,浮体随着绕流的不对称性会沿着系统轴线发生蛇形摆动,浮体之间会产生相互干扰,由于系泊缆的约束作用,浮体系统的运动主要表现为横向运动,为了验证数值模拟方法的可靠性和精确度,选取物理模型的实测波要素和水流速度,采用全耦合时域分析方法对系泊多浮体耦合系统进行了数值计算。

3.1 单纯波浪作用下系统的运动响应和动力响应

对系泊缆绳施加40 kN的预张力,浮筒之间的连接缆绳不施加预张力,在不同有义波高和谱峰周期的单纯波浪作用下,全耦合时域计算时间设置为600 s。在有义波高为0.85 m的单纯波浪作用下,系泊浮筒链系统的浮体横向运动时间历程曲线如图4所示,缆绳张力时间历程曲线如图5所示。

图 4 浮筒横向运动时域曲线 Fig. 4 Time domain curve of buoy lateral motion

图 5 缆绳张力时域曲线 Fig. 5 Time domain curve of cable tension

观察图4可知在0.85 m有义波高的单纯波浪作用下,浮筒沿着轴线做往复运动,最大横向位移出现在200 s附近,达到1.17 m,横向位移的平均值为0.17 m,表明浮筒基本可以回到初始位置。由图5(a)可知:浮筒前后往复运动的同时也在随着波峰和波谷上下运动,出现端部系泊缆张力小于预应力的情况,浮筒在前后往复运动时,端部系泊缆的张力随着浮筒横向运动分量的增大而增大,最大值为152.3 kN。由于浮筒链随着不对称绕流发生蛇形摆动,在端部系泊缆的约束和各浮筒之间的耦合作用下,中间浮筒的横向位移最大,导致中间浮筒两端的缆绳张力比其他浮筒间缆绳张力大,图5(b)图5(c)描述了2号和3号缆绳张力的时域曲线,最大值分别为141.6 kN和133.3 kN。观察图5(d)图5(g)发现:中间浮筒迎浪侧的4号和5号系泊缆张力的最大值为54 kN和56 kN,6号系泊缆的最大值为51 kN,由于浮筒链向后运动的最大位移为0.85,所以背浪侧缆绳的最大张力相对较小,为46.7 kN,在保证不走锚的情况下,要重点关注中部位置浮筒迎浪侧的系泊缆张力。结合图4图5发现:缆绳的最大张力与最大横向位移均出现在200~250 s之间。

不同波浪作用下浮筒横移量的试验和数值模拟结果见表2,试验和数值模拟得到的系泊缆张力见表3

表 2 单纯波浪作用下,中间浮筒横向运动分量最大值的试验和模拟对比 Tab.2 Experimental and simulation comparisons of the maximum lateral motion component of middle buoy (only wave)

表 3 单纯波浪作用下,缆绳张力最大值的试验和模拟对比 Tab.3 Experimental and simulation comparisons of the maximum tension of cable (only wave)

分析可知:采用全耦合时域分析模拟得到的浮筒横向运动分量的最大值与模型试验得到的结果较为接近,浮筒间缆绳张力的试验值和模拟值也普遍较为接近。物理试验中安装了21个浮筒,为了节省计算时间,水动力模型只选取了中间位置的7个浮筒,通过增大端部系泊缆的长度代替减少的浮筒数量,从而使水动力模型和物理试验模型的跨度保持一致,这也导致水动力模型系统的稳定性降低。全耦合时域分析得到的背浪方向浮筒横向运动分量的平均值普遍大于模型试验的平均值,这是由于水动力模型系统的稳定性没有物理模型试验的稳定性好,导致逃逸波浪和50年重现期波浪作用时背浪方向浮筒位移最大值时域分析与试验结果相差较大,迎浪侧5号缆绳张力时域分析最大值与试验最大值也相差较大。

3.2 波流共同作用下系统的运动响应和动力响应

考虑到允许工程船舶作业的海洋环境,且当地海流流速不会超过1.2 m/s[12],因此只考虑1.2 m/s的海流和2种作业波浪的共同作用下,对系泊浮筒链进行600 s的全耦合时域计算,系泊浮筒链系统的浮体横向运动时间历程曲线如图6所示,浮筒横向运动分量最大值的实验和模拟对比如表4所示,在50年重现期波浪的作用下,缆绳张力最大值的实验和模拟对比如表5所示。

图 6 波流共同作用下,浮筒横向运动时域曲线 Fig. 6 Time domain curve of buoy lateral motion (wave+1.2 m/s current)

表 4 波流共同作用下,浮筒横向运动分量最大值的的实验和模拟对比 Tab.4 Experimental and simulation comparisons of the maximum lateral motion component of middle buoy (wave+1.2 m/s current)

表 5 波流共同作用下,缆绳张力最大值的实验和模拟对比 Tab.5 Experimental and simulation comparisons of the maximum tension of cable (wave+1.2 m/s current)

观察图6可知:在波浪和水流的耦合作用下,浮筒发生前后摆动,在水流的冲击作用下沿水流前进方向运动量较大,而且不能回到初始位置;结合表4得到,浮筒横向位移的全耦合时域分析最大值与模型试验最大值较为接近,可以较为准确的模拟出浮筒在波流共同作用时横向位移的最大值,但是数值计算最小值和平均值大于物理模型试验平均值,这是由于水动力模型系统的稳定性没有物理模型试验的稳定性好,导致数值计算高估了系统横向运动分量的最小值和平均值。表5给出了系泊浮筒最大缆绳张力的汇总,出现在50年重现期波浪和1.2 m/s的海流同时作用下,依然是端部缆绳的张力最大,为755.6 kN,可以为锚固块体重量的设计提供参考依据。

4 结 语

本文采用全耦合时域分析方法,通过AQWA软件对浮筒链式系泊多浮体耦合系统进行了运动和动力响应数值模拟,并与现有物理模型试验进行了对比分析,得到下列结论:

1)在给定海洋动力要素的情况下,浮体横向位移的全耦合时域分析结果与模型试验结果较为接近,浮筒间缆绳张力的试验值和模拟值普遍较为接近。

2)由于水动力模型系统的稳定性没有物理模型试验的稳定性好,导致数值计算高估系统横向运动分量的最小值和平均值。

3)通过全耦合时域分析结果与物理模型试验结果的对比,分析了影响数值模拟结果与物理模型试验结果产生误差的原因,验证了全耦合时域分析方法的正确性,对系泊多浮体耦合系统的运动和动力响应的计算具有较好的精确度。

参考文献
[1]
沈庆, 陈徐均. 系泊多浮体系统波浪运动响应的动力学分析[J]. 解放军理工大学学报, 2000, 1(4): 31-36.
SHEN Qing, CHEN Xu-jun. Dynamic analysis of moored multi-body floating system’s responses to regular waves[J]. Journal of PLA University of Science and Technology, 2000, 1(4): 31-36.
[2]
沈庆, 陈徐均. 系泊多浮体系统流固耦合和浮体间耦合动力学分析[J]. 中国造船, 2002, 43(2): 81-84.
SHEN Qing, CHEN Xu-jun. Dynamic analysis of mooring multi-body system coupled with fluid and among bodies[J]. Ship Building of China, 2002, 43(2): 81-84. DOI:10.3969/j.issn.1000-4882.2002.02.012
[3]
王翔, 沈庆, 王广东, 等. 锚泊多浮体系统动态响应[J]. 解放军理工大学学报(自然科学版), 2007, 8(4): 368-374.
WANG Xiang, SHEN Qing, WANG Guang-dong, et al. Dynamic response of muti-floating-body system[J]. Journal of PLA University of Science and Technology, 2007, 8(4): 368-374.
[4]
王桂波, 勾莹, 滕斌, 等. 铰接多浮体系统在规则波作用下的运动响应[J]. 大连理工大学学报, 2014, 54(6): 618-625.
WANG Gui-bo, GOU Ying, TENG Bin, et al. Motion responses of hinged multiple floating bodies under regular wave action[J]. Journal of Dalian University of Technology, 2014, 54(6): 618-625.
[5]
袁培银, 赵宇, 郭建廷, 等. 多浮体式系泊系统设计及水动力性能分析[J]. 舰船科学技术, 2016, 38(2): 80-84.
YUAN Pei-yin, ZHAO Yu, GUO Jian-ting, et al. The design method of a new type of multi-body mooring system and the analysis of its hydrodynamic performance[J]. Ship Science and Technology, 2016, 38(2): 80-84. DOI:10.3404/j.issn.1672-7649.2016.2.016
[6]
高峰, 李焱, 沈文君, 等. 风浪流共同作用下的多浮体系泊性能试验研究[J]. 水道港口, 2017, 38(6): 556-560.
GAO Feng, LI Yan, SHEN Wen-jun, et al. The study on mooring characteristic for multi-body floating system under the action of wind, wave and flow[J]. Journal of Waterway and Harbor, 2017, 38(6): 556-560.
[7]
于文太, 辛凤娟, 梁学先, 等. 多浮筒悬链线系泊系统计算软件的开发及试验[J]. 船舶工程, 2017, 39(11): 31-36.
[8]
王贝壳, 陈涛, 杨黎明, 等. 非通航孔桥墩自适应拦截网防撞装置实船拦截试验与水动力计算[J]. 海洋工程, 2017, 35(1): 90-96.
[9]
夏雪, 张可成, 伏耀华, 等. 深水海域非通航孔桥梁防撞数值模拟[J]. 中国航海, 2018, 41(1): 84-87.
XIA Xue, ZHANG Ke-cheng, FU Yao-hua, et al. Numerical simulation of anti-collision design of non navigable spans of deep sea bridge[J]. Navigation of China, 2018, 41(1): 84-87. DOI:10.3969/j.issn.1000-4653.2018.01.017
[10]
张周康, 洪亮. 多浮体系统在波浪作用下运动响应的数值仿真研究[J]. 机械与电子, 2018, 36(1): 10-14.
ZHANG Zhou-kang, HONG Liang. Numerical simulation of multibody floating system’s motion response to waves[J]. Machinery & Electronics, 2018, 36(1): 10-14. DOI:10.3969/j.issn.1001-2257.2018.01.003
[11]
张周康, 洪亮. 斜浪下多浮体多浮体结构物运动响应的仿真分析[J]. 兵器装备工程学报, 2018, 39(3): 182-186. DOI:10.11809/bqzbgcxb2018.03.038
[12]
王鑫. 系泊浮体链式防撞系统水动力特性和防撞性能试验研究[D]. 大连: 大连理工大学, 2013.
[13]
高巍. ANSYS AQWA软件入门与提高[M]. 北京: 中国水利水电出版社, 2018.
[14]
倪志辉, 吴立春. 浅海潮流流速剖面参数的相位分析[J]. 人民黄河, 2009, 31(6): 53-57.
NI Zhi-hui, WU Li-chun. Phase analysis of profile parameters of tidal current velocity in shallow sea[J]. Yellow River, 2009, 31(6): 53-57. DOI:10.3969/j.issn.1000-1379.2009.06.021