2. Guangzhou Marine Engineering Corporation, Guangzhou 510250, China
三体船作为水面舰船的一种新船型,在稳性、耐波性、快速性和甲板面积等方面相对于常规船型而言有许多优势。当前针对其开展的研究工作方兴未艾,但主要的研究工作是在船型开发[1, 2]、阻力[3]和耐波性[4, 5, 6, 7]方面。三体船的波浪载荷对结构设计与强度校核非常关键,但是由于波浪载荷这个问题比较复杂,国内外对其的研究工作还很有限,其中Bingham等[8]采用脉动源、移动-脉动源和频域水弹性法计算了三体船在迎浪规则波中的运动和载荷,英国劳氏船级社虽然也发布了三体船规范指导该种船型的建造和入级工作[9],但和具有一百多年历史的单体船设计规范相比,人们对高性能三体船的波浪载荷的研究还不充分,相应的规范设计法也很不成熟。国内,徐敏等[10]采用基于时域Rankine源法的软件计算了三体船的波浪载荷,并且对三体船布局进行了优化;汪雪良等[11]研究了三体船横向连接桥结构的波浪载荷,并将规范值和实验值做了比较;蒋清华等[12]采用三维频域直接计算法研究了某三体船的总纵强度;本课题组已经开发了基于刚体理论的三体船运动和波浪载荷计算方法,集成到了COMPASS-WACLS[13]软件系统中。由于三体船型对安全性和重量控制要求高,其独体的连接桥结构以及细长高速和首部外飘特性,在其船体结响应中包括颤振和弹振,因此在三体船载荷数值计算时考虑水动力的非线性特性和砰击效应的显得非常重要,开发一套非线性时域流固耦合分析方法以预报三体船波浪载荷是亟需解决的问题。
本文首先介绍波浪载荷计算的三维时域水弹性理论,采用有限元法对三体船进行模态分析,然后计算了三体船在迎浪当中的垂荡、纵摇以及船舯处的垂向波浪弯矩和剪力的时历;然后按照劳氏规范计算了该三体船的垂向设计波浪载荷,然后水弹性计算结果和长期预报值以及规范值进行比较分析。
1 波浪载荷三维时域水弹性理论 1.1 三体船模态分析
采用水弹性理论分析船舶的运动和波浪载荷时,首先需要对船体结构进行模态计算以获得结构自由振动时模态特征值。目前,对船体结构进行模态分析的方法有解析法(例如迁移矩阵法)和数值有限元法,有限元法又可以采用不同的计算模型(例如梁理论和三维结构模型)。本文选取船体干模态作为模态分析的主坐标,采用有限元法进行船体结构振动方程的求解。由于只是考虑了三体船迎浪状况下的响应,以下给出垂向振动的控制方程式:
$\left\{ \begin{matrix}
\mu \left( x \right)w\left( x,t \right)=V_{z}^{'}\left( x,t \right)+{{f}_{y}}\left( x,t \right) \\
{{J}_{y}}{{\theta }_{y}}\left( x,t \right)=V_{y}^{'}\left( x,t \right)-{{V}_{z}}\left( x,t \right) \\
\end{matrix} \right.$
(1)
$\left\{ \begin{matrix}
{{V}_{Z}}\left( x,t \right)=kG{{A}_{V}}\left( x \right)\left[ \gamma \left( x,t \right)+{{\alpha }_{V}}\left( x \right){{{\dot{\gamma }}}_{Z}}\left( x,t \right) \right] \\
{{M}_{y}}\left( x,t \right)=E{{I}_{V}}\left( x \right)\left[ {{\theta }^{'}}\left( x,t \right)+{{\beta }_{V}}\left( x \right){{{\dot{\theta }}}_{y}}\left( x,t \right) \right] \\
\end{matrix} \right.$
(2)
w′(x,t)=γZ(x,t)-θy(x,t)
(3)
式(1)~(3)构成了垂向运动控制方程。根据上述控制方程,以及两端自由梁的边界条件,采用适当的数值方法可以确定船体梁自由震动的固有频率和振型函数。图1给出了用于模态分析的三体船有限元模型,在图2中给出了前面几节点振型图。
1.2 非线性水弹性力学方程的建立与求解
在时域内,船体运动的非线性水弹性力学方程可以写作如下形式:
(a+A)${\ddot{p}}$r(t)+(b+B)
${\dot{p}}$r(t)+(c+C)
${\ddot{p}}$r(t)
=
FI(t)+FD(t)+Fsiam(t)
(4)
若将船体在波浪上的运动视作在波浪激励下的强迫振动,此时式(4)可以转化为
a${\ddot{p}}$r(t)+
b${\dot{p}}$r(t)+cpr(t)=F(t)
(5)
当船体在大幅波中摇荡运动时,某一时刻作用在船体上的流体载荷F(t)可以写作下式:
F(t)=Fs(t)+FI(t)+FR(t)+
FD(t)+Fsiam(t)
(6)
1)静水恢复力
船体瞬时平均湿表面上的流体静力载荷为
FS(t)=-ρg$\sum\limits_{k=1}^{m}{{{p}_{ka}}\iint\limits_{S\left( t \right)}{\bar{n}\cdot {{{\bar{u}}}_{r}}{{w}_{x}}ds=-}}\sum\limits_{k=1}^{m}{{{R}_{rk}}\left( t \right)}{{p}_{ka}}\left( t \right)$
(7)
2)入射波力
船体瞬时平均湿表面上的入射波力表达式为
FI(t)=-ρζa$\iint\limits_{S\left( t \right)}{\bar{n}\cdot {{{\bar{u}}}_{r}}}$(iω-U$\frac{\partial }{\partial x}$φIds
(8)
FI(t)=-ρζaiω0$\iint\limits_{S\left( t \right)}{\bar{n}\cdot {{{\bar{u}}}_{r}}}$φIds(r=1,2,…,m)
(9)
3)波浪砰击力
当船舶在波浪中航行时,会发生底部砰击、外飘砰击,由于三体船独特的构型,也会产生严重的连接甲板砰击。对底部砰击载荷进行预报,首先需要采用一定的方法判断砰击发生的条件。底部砰击发生时剖面的最大压力可表示为
${{p}_{\max }}=\frac{1}{2}\rho {{k}_{1}}\dot{Z}_{R0}^{2}$
(10)
Δ=$\frac{7.94}{{{10}^{3}}}\sqrt{L}$ | (11) |
把不同时刻各点压力的垂向分量沿整个剖面积分,即可得到砰击载荷Fb(xb,t)。
外张砰击载荷预报方法包括动量砰击理论和瓦格纳型砰击理论,分别适合于外张砰击力和外张砰击压力的计算,本文采用前者。按动量砰击理论有
Ff(xb,t)=$\left\{ \begin{matrix}
\frac{d{{m}_{\infty }}}{dz}\left( \frac{d{{Z}_{R}}}{dt} \right),\frac{d{{Z}_{R}}}{dt}\le 0 \\
0,\frac{d{{Z}_{R}}}{dt}>0 \\
\end{matrix} \right.$
(12)
Fslam(t)$\int\limits_{L}^{{}}{{}}$Ff(xb,t)w(x)dxb
(13)
对于连接桥砰击,本文作者在该方面做了大量的研究,具体方式可以参考文献[16, 17].
1.2.2 船体非线性运动方程的求解
通过以上的分析,船体在规则波中航行的运动微分方程为
a${\ddot{p}}$r(t)+b
${\dot{p}}$r(t)+cpr(t)=F(t)
(14)
同频域方法相同,由式(14)解出主坐标后,利用模态叠加原理,就可以得到船体结构的位移w(x,t)、弯矩M(x,t)、动剪切力V(x,t):
$\left\{ \begin{matrix}
w\left( x,t \right)=\sum\limits_{r=0}^{m}{{{p}_{r}}\left( t \right){{w}_{r}}\left( x \right)} \\
M\left( x,t \right)=\sum\limits_{r=0}^{m}{{{p}_{r}}\left( t \right){{M}_{r}}\left( x \right)} \\
V\left( x,t \right)=\sum\limits_{r=0}^{m}{{{p}_{r}}\left( t \right){{V}_{r}}\left( x \right)} \\
\end{matrix} \right.$
(15)
基于三维时域非线性水弹性程序,采用非线性设计波法进行极限、巡航工况中不同波高下规则波弯矩的响应曲线绘制,应用线性长期分析结果进行等效设计波波高的截取,采用得到的等效设计波,应用三维时域非线性水弹性程序在不同波长下进行计算,等到非线性设计波下的响应最大值,将响应最大值作为设计载荷用于强度评估。其中极限工况的参数:长期预报时航速取为巡航速度,垂向波浪弯矩超越概率值取为10-8,得到等效设计波波高h1,计算航速可取为较低航速;巡航工况的参数:长期预报时航速取为巡航速度,垂向波浪弯矩超越概率值取为10-5,得到等效设计波波高h2,计算航速可取为实际巡航速度。本计算实例的三体船主要参数如表1。
由图3可以得知,在巡航工况下,三体船的纵摇和横摇呈现规则简谐运动,而垂向波浪弯矩会遭遇高频载荷;在极限工况中,三体船的纵摇和横摇呈现明显的非线性,在垂向波浪弯矩中叠加明显的砰击弯矩。
3 波浪载荷的规范计算以英国国防部的资金支持和英国QinetiQ公司的技术支持为基础,英国劳氏船级社(LR)针对三体船这一特殊船型研发出劳氏三体船规范用于三体船的建造和入级工作。该规范对三体船的环境载荷等的若干方面进行了详细叙述,并规定采用直接计算程序作为新造船舶的设计依据。本节该规范对三体船船体梁纵向波浪载荷进行计算。将该计算结果与长期预报值、时域水弹性对比分析,对比结果可为理论计算方法的改进与规范计算公式的修正提供参考。规范中明确给出了三体船所受各项外载荷的计算公式。对于三体船纵向弯曲,分为中拱和中垂两种情况,但是三体船不仅存在主船体中拱和中垂的情况,还存在侧体中拱和中垂的情况,如图4、5所示。
所受的纵向波浪弯矩Mw的计算公式为
Mw=Ff·Df·M0
(16)
${{F}_{fs}}\left\{ \begin{matrix}
1.1R_{A}^{0.3}{{R}_{A}}\ge 1 \\
-1.1{{R}_{A}}<1 \\
\end{matrix} \right.;{{R}_{A}}=\frac{300\left( {{A}_{BF}}-0.5{{A}_{SF}} \right)}{{{L}_{R}}{{B}_{wl}}};$
${{D}_{f}}\left\{ \begin{matrix}
0,尾垂线及以后区域 \\
1,\left( 0.35-\frac{{{x}_{sh}}}{4{{L}_{R}}} \right){{L}_{R}}\tilde{\ }\left( 0.6-\frac{{{x}_{sh}}}{4{{L}_{R}}} \right) \\
0,首垂线及前区域 \\
\end{matrix} \right.$
${{L}_{f}}\left\{ \begin{matrix}
\begin{matrix}
0.041 & 2{{L}_{R}}-4.0,{{L}_{R}}\le 90m \\
\end{matrix} \\
0.75-{{\left( \frac{300-{{L}_{R}}}{100} \right)}^{1.5}},{{L}_{R}}>90m \\
\end{matrix} \right.$
把三体船的水弹性时域计算结果,与三维理论预报结果以及三体船规范结果进行比较分析,其中计算了巡航、极限两种工况最后把计算统计结果同三体船模型试验结果对比分析,对比结果列于表2。
kN·m | ||
数据源 | 数值 | |
船舯垂向弯矩中拱 | 船舯垂向 弯矩中垂 | |
劳氏规范值 | 250 990 | -403 452 |
水弹性计算值(巡航) | 302 022 | -444 860 |
水弹性计算值(极限) | 501 432 | -565 291 |
长期预报值(巡航) | 268 200 | -268 200 |
长期预报值(极限) | 365 900 | -365 900 |
从对比结果可以看出,按照10-8超越概率的长期预报值是基于线性理论的,不管是对于中拱还是中垂工况,其弯矩值值是相等的,无法考虑砰击以及其他非线性流体载荷影响;水弹性方法的计算出中拱幅值小于中垂幅值,充分体现实际三体船在海浪中的波浪载荷响应值;考虑了砰击效应和水动力非线性效应的水弹性方法在巡航工况时和劳氏三体船载荷的值比较接近,吻合较好;在水弹性计算结果中由表2可知长期预报值在巡航工况下的值比劳氏规范值略大,这个正好符合船体非线性的载荷特征。在极限工况中的水弹性预报的中拱和中垂弯矩均较巡航工况要大,并且比规范值更大;在极限工况中的水弹性预报的中拱和中垂弯矩相对误差为12.7%,说明在高海况中(极限工况)本文的时域水弹性理论也能给出较好的计算结果。
5 结论通过三维时域水弹性理论对三体船的垂向波浪设计载荷的研究,本文得出以下主要结论:
1)三体船的模态分析分析可以得知其结构比较细长,固有频率较低,在计算其纵向波浪载荷时计入水弹性效应非常重要;
2)本文三维时域水弹性方法可以考虑非线性的入射波力、静水恢复力和砰击力(底部砰击、外飘砰击和连接桥砰击,即使在高海况中也能给出相对合理的结果,其大小与三维水弹性计算时的巡航工况下大致相等;
3)英国劳氏三体船规范的计算公式能够反映出纵向的垂向波浪弯矩特性;
4)由于三体船规范计算值与巡航工况时的水弹性计算结果,但是远小于极限工况时的幅值,而极限工况也是三体船会遭遇的工况,建议在确定三体船设计载荷时,应该考虑极限波浪弯矩的计算。
[1] | 黄德波. 若干高速与高性能排水式舰船船型的近期发展[C]//第七届全国水动力学学术会议暨第十九届全国水动力学研讨会文集(上册). 哈尔滨, 中国,2005:1-8. |
[2] | 李培勇. 多体船型的性能及概念设计研究[D]. 上海:上海交通大学, 2002. LI Peiyong. Characteristic and concept design study of multihull[D]. Shanghai:Shanghai Jiao Tong University, 2002. |
[3] | 周广利. 三体船型阻力预报、优化与系列性试验分析研究[D]. 哈尔滨:哈尔滨工程大学, 2010. ZHOU Guangli. Resistance prediction, hull optimization and serial model tests of trimarans[D]. Harbin:Harbin Engineering University, 2010. |
[4] | 李培勇, 冯铁城, 裘泳铭. 三体船横摇运动[J]. 中国造船, 2003, 44(1):24-30. LI Peiyong, FENG Tiecheng, QIU Yongming. Investigation of trimarans roll motion characteristics[J]. Shipbuilding of China, 2003, 44(1):24-30. |
[5] | ZHANG J W. Roll damping characteristics of a trimaran displacement ship[J]. International Shipbuilding Progress, 1999, 46(448):445-472. |
[6] | FRANCESCUTO A. On the roll motion of a trimaran in beam waves[C]//Proceedings of the 11th International Offshore and Polar Engineering Conference (ISOPE 2001). Stavanger, Norway, 2001. |
[7] | DAVIS M R, Holloway D S. A comparison of the motions of trimarans, catamarans and monohulls[J]. Australian Journal of Mechanical Engineering, 2007, 4(2):183-196. |
[8] | BINGHAM A E, HAMPSHIRE J K, TEMAREL P, et al. Motions and loads of a trimaran travelling in regular waves[C]//Proceedings of the 6th International Conference on Fast Sea Transportation. Southampton, 2001. |
[9] | Lloyd's Register (LR). Incorporating Notice No.1, The rules for the classification of trimaran[S]. 2006. |
[10] | 徐敏, 张世联. 三体船波浪载荷与片体布局优化[J]. 上海交通大学学报, 2011, 45(4):470-474. XU Min, ZHANG Shilian. The wave loads and outrigger hulls' position optimization for trimaran traveling in waves[J]. Journal of Shanghai Jiao Tong University, 2011, 45(4):470-474. |
[11] | 汪雪良, 胡嘉骏, 顾学康, 等. 三体船横向结构波浪设计载荷试验与规范比较研究[J]. 船舶力学, 2011, 15(3):269-275.WANG Xueliang, HU Jiajun, GU Xuekang, et al. Comparative studies of the transverse structure design wave loads for a trimaran by model tests and rule calculations[J]. Journal of Ship Mechanics, 2011, 15(3):269-275. |
[12] | 周清华, 蒋学敏, 李祥宁. 基于波浪载荷直接计算的三体船总纵强度评估研究[J]. 船舶工程, 2012, 34(2):17-20. ZHOU Qinghua, JIANG Xuemin, LI Xiangning. Research on longitudinal strength assessment of trimaran based on wave loads direct calculation method[J]. Ship Engineering, 2012, 34(2):17-20. |
[13] | CCS. Compass wacls user manual[Z]. Beijing:CCS, 2013. |
[14] | KUMAI T. Damping factors in the higher modes of ship vibrations[J]. European Shipbuilding, 1958, 7(1):29-34. |
[15] | OCHI M K, MOTTER L E. Prediction of slamming characteristics and hull response for ship design[C]//Annual Meeting of Society of Naval Architects and Marine Engineers. Jersey City, USA, 1973:144-176. |
[16] | REN Huilong, CHEN Liangliang, FENG Guoqing, et al. Cross-deck slamming load calculation of trimaran based on modified plate slamming theory[J]. Applied Mechanics and Materials, 2013(477/478):325-329. |
[17] | 戴仰山, 沈进威, 宋竞正. 船舶波浪载荷[M]. 北京:国防工业出版社,2007. DAI Yangshan, SHEN Jinwei, SONG Jingzheng. Ship wave loads[M]. Beijing:National Defense Industry Press, 2007. |