-
根据中国船级社(CCS)《海上高速船入级及建造规范》[1]中的定义,气垫船是指船舶不论在静止还是运动时,其全部重量或大部分重量能被连续产生的气垫所支承的船舶。气垫船的航行阻力很小,可使航速高达60~80 km/h,可以完成多种使命任务,在军用和民用领域均有着十分广阔的发展潜力和应用前景。
从上世纪80~90年代开始,各研究机构开始意识到在船舶设计初期对船体结构进行疲劳强度评估的重要性。用于评估的方法主要有简化计算方法和直接计算方法2种,其中简化计算方法因简单易行并能大致估算出典型节点的疲劳寿命,一般用于初步估算。气垫船一般为铝合金制造,尺度较小,其刚度较一般船舶相对较弱;此外,该船特有的航行方式、垫升风机等因素也使得气垫船的振动疲劳问题较常规船舶更为严重[2]。但是,由于气垫船的运行原理不同于一般船型,其受力较特殊,故船体弯矩、加速度、损伤度的计算等势必与普通船舶相比有很大的差异。目前,国际上各船级社提出的疲劳强度简化计算方法并不完全适用于气垫船,并且学者们的注意力主要集中在局部振动强度、振动噪声和总体强度等问题上[2-5],针对主艇体疲劳强度的研究较少,而各船级社给出的疲劳强度规范又难以直接用于气垫船,故本文将参考CCS的《船体结构疲劳强度指南》[6]和《海上高速船入级及建造规范》[7]、国家军用标准《水面舰艇结构设计计算方法》[8],以及及国外船级社的相关疲劳规范[9-10],在现有的计算常规船型疲劳强度规范的基础上,提出一种适用于气垫船疲劳强度计算的简化方法,并以某气垫船为例进行验证分析。
-
谱分析法及简化计算方法是疲劳强度分析最常用的2种方法,其主要是基于S-N曲线(应力寿命曲线,通常由试验测得)和线性疲劳累积损伤理论。疲劳损伤度D及S-N曲线的表达和计算如图1、图2及式(1)和式(2)所示。图2中,σr为材料的疲劳极限,在该应力水平下,材料可以承受无限次应力循环。
$$N{S^m}{\rm{ = }}A$$ (1) $$D{\rm{ = }}\frac{{{n_1}}}{{{N_1}}} + \frac{{{n_2}}}{{{N_2}}} + \frac{{{n_3}}}{{{N_3}}} + ... = \sum\limits_i^K {\frac{{{n_i}}}{{{N_i}}}} $$ (2) 式中:S为应力范围;N为循环次数;Ni为应力范围为Si时结构发生疲劳破坏的循环次数;ni为应力循环实际发生的次数;m,A为试验确定的S-N曲线参数。
如果船体在设计寿命内遭受的交变应力的长期分布符合f (S)分布,根据积分思想并结合S-N曲线表达式,损伤度的计算可按式(3)计算。
$$D{\rm{ = }}\int {\frac{{{\rm{d}}n}}{N}} {\rm{ = }}\int {\frac{{{N_{\rm{L}}}f(S){\rm{d}}S}}{N} = } {N_{\rm{L}}}\int {\frac{{f(S){\rm{d}}S}}{N} = } \frac{{{N_{\rm{L}}}}}{A}\int {\frac{{f(S){\rm{d}}S}}{N}} $$ (3) 式中:NL为相应时间内应力范围的循环次数;dn为应力范围微元dS中的载荷循环次数。
各规范大多是以此为基础进行相关规定。结合上式可以看出,若能得到所考虑的设计寿命周期内船体结构所受外界交变载荷的长期分布f (S),要计算出船体结构的疲劳累积损伤程度以及寿命并不困难。然而,由于海况的不确定性以及船体的复杂性,难以直接获得这一长期分布,但结合大量的实船数据,发现该长期分布近似符合双参数的Weibull分布,故而在简化计算方法中采用Weibull分布对f (S)进行拟合。Weibull分布包含2个参数,形状参数与尺度参数,其中形状参数一般回归为与船体主尺度相关的表达式,尺度参数可根据超越概率的定义计算求得。
对于疲劳强度简化计算,现行各规范虽然在具体参数的选取及所用简化公式上有所不同,但大体流程是一致的。可根据船体梁弯矩计算梁应力,根据局部载荷计算局部应力,据此得到应力范围,然后再根据线性累积损伤理论并选取合适的S-N曲线计算损伤度及疲劳寿命[11]。
-
气垫船的工作状态有多种:垫升高速航行状态、排水低速航行状态、登陆垫着陆工作状态等,根据各工作状态对疲劳损伤度的贡献,本文将排水低速航行状态及垫升高速航行状态纳入了校核当中。气垫船的排水航行状态与一般排水型船舶没有本质的差别。但采用垫升状态航行时,其载荷情况有着很大的不同:首先,引起船体梁应力的波浪弯矩将由垫升气压产生的垫升弯矩所代替;其次,船底板受到的垫升压力、气道因为垫升而受到的压力等气垫船垫升状态所独有的载荷也需要予以考虑[12]。根据上述分析,结合气垫船的特点,确定排水航行及垫升航行状态下的疲劳强度校核流程分别如图3和图4所示。
根据气垫船实际工作中各状态的时间占比,确定总在航率为0.7,其中垫升状态为0.63,排水状态为0.07。
航行状态确定之后,共选取迎浪、横浪和斜浪等6个工况。在每个工况中,各总体载荷成分按比例组成,且有一个载荷成分达到最大值,以作为该工况的主要载荷。参考规范并结合气垫船的航行状态,得到各载荷分量的分配系数如表1所示。
表 1 各工况下各载荷分量的分配系数
Table 1. Distribution coefficients of each load component under different working conditions
浪向 工况 静水(垫升)弯矩 中垂弯矩 中拱弯矩 水平弯矩 迎浪 1 1 1 0 0 2 1 0 −1 0 斜浪 3 1 0.2 0 −0.2 4 1 0 −0.2 0.2 横浪 5 1 0.1 0 −1 6 1 0 −0.1 1 -
根据实船数据确定的Weibull分布的形状参数对计算结果以及疲劳评估的准确性影响很大,目前,各规范对此的规定也各不相同,例如,CCS疲劳指南和共同规范(HCSR规范)均采用形状参数为1、超越概率为10−2的做法进行疲劳强度计算,挪威船级社(DNV)规范的形状参数则根据不同的情况进行计算,并对应10−4的超越概率。根据文献[13],如果形状参数选取得足够准确,那么任何超越概率都可以选择,但当形状参数h取为1、对应的超越概率为10−2时,所得结果能在最大程度上减小因形状参数不准确所带来的误差。由于缺少大量的气垫船实船数据,故选取形状参数h=1、超越概率为10−2进行后续的疲劳强度简化计算,这与CCS规范以及共同规范是一致的。根据波浪载荷相关知识,我们假设在船舶寿命周期内会遭遇108次应力循环,即超越概率为10−8,同时,引入概率水平转换因子fp,计算式如下:
$${f_{\rm{p}}} = \frac{{\ln {{10}^2}}}{{\ln {{10}^8}}} = 0.25$$ (4) -
疲劳强度是一个一段时间内的概念。随着使用年限的增长,骨材、板材等势必会发生腐蚀,从而达不到最开始的设计要求,所以评估疲劳强度时应考虑船体结构工作运营时产生的磨损腐蚀的影响,并结合CCS规范,在计算热点应力时引入腐蚀修正系数。
1) 行疲劳强度简化计算时,对于船体梁弯曲正应力的热点应力,腐蚀修正系数
$ {f}_{\mathrm{c}\mathrm{h}}=1.05 $ 。2) 进行疲劳强度简化计算时,对于侧向载荷作用下的弯曲正应力的热点应力,腐蚀修正系数
$ {f}_{\mathrm{c}\mathrm{h}}=1.1 $ 。需要说明的是,CCS规范给出的腐蚀修正系数主要用于设计年限为20年的钢质船,对于铝合金制气垫船,目前还没有规范给出有关其腐蚀余量或修正的计算方法。铝制气垫船的设计年限一般约为15年,且其材质的抗腐蚀性相对较好,所以在对材料的具体特性了解不多的情况下,按照上述说明选取系数对材料进行腐蚀修正不失为一种偏于安全的做法。
-
1) 排水航行加速度。
当气垫船以排水状态航行时,其加速度可以借鉴一般常规排水型船舶的计算。垂向加速度的计算如式(5)和式(6)[6]所示。
$$ {a_0} = \left( {\frac{{3C}}{L} + \frac{{{C_{\rm{V}}}V}}{{\sqrt L }}} \right){f_{\rm{r}}}{f_{\rm{p}}} $$ (5) $$ {a_z} = 7{a_0}{\sqrt {{C_{\rm{b}}}}} $$ (6) 式中:
$ {a}_{0} $ 为加速度常数;$ {C}_{\mathrm{V}} $ 为航速参数,可按$ \sqrt{L}/50 $ 计算求得且不大于0.2;L为船长;V为服务航速(需特别注意,此处为排水低速航行时的航速),kn;C为波浪系数;fp为概率水平转换因子,取为0.25;fr为航区系数,结合气垫船工作的海域,近海取0.9,沿海取0.85;$ {a}_{z} $ 为垂向加速度;$ {C}_{\mathrm{b}} $ 为方形系数。结合相关试验测得加速度数据,可知中垂状态加速度值应为中拱状态的1.5倍左右。
2) 垫升航行加速度。
当气垫船以垫升状态航行时,可按式(7)计算重心位置的加速度acg[1]:
$${a_{{\rm{cg}}}} = \left( {0.44\sqrt {{H_{1/3}}} + 0.022} \right)\frac{{{V_H}}}{{\sqrt[3]{\varDelta }}}g{f_{\rm{r}}}{f_{\rm{p}}}$$ (7) 式中:
$ {H}_{1/3} $ 为气垫船航行海域波浪的有义波高;$ {V}_{H} $ 为在该有义波高波浪中航行时的速度,kn;$\varDelta$ 为排水量,t;g为重力加速度,m/s2。任意位置的垂向加速度az按式(8)[1]计算:
$${a_{{z}}} = {k_{\rm{a}}}{a_{{\rm{cg}}}}$$ (8) 式中,
$ {k}_{\mathrm{a}} $ 为垂向加速度分布系数,舯后取1.0,艏部取2.0,船艏与船舯之间线性插值取值。 -
气垫船排水状态下的载荷计算可借鉴一般排水型船舶的计算方法。该船与一般排水型船相比最大的区别在于船体上环绕有一圈围裙。在气垫船的实际运动过程中,围裙的运动或是变形都会带动周围的水运动,其对疲劳载荷的影响不可忽略,因此针对排水航行状态下的气垫船,结合试验提出了一种简化模型方案,即采用一定长度的刚性附体,使其与围裙的作用效果相近。简化示意图如图5和图6所示。
之后,利用三维频域势流理论方法(COMPASS-WALCS),对带有附体的简化模型的运动响应及波浪载荷进行预报。类似地,也可以采用相同的方法对2.2.3节所述排水状态的加速度(不限于垂向加速度)进行预报计算。通过上述“简化”方法建立模型并提交计算,得到某剖面的垂向载荷,将其无因次化之后与试验所得数据进行对比,具体如图7所示。图中,λ/L为波长船长比,
${\bar M_{\rm{V}}}$ 为垂向波浪弯矩的无因次量,试验1、试验2为重复性验证试验。通过对比结果,可知与没有附体(模型2)的结果相比,考虑附体(模型1)的结果更接近于试验值,证明了该方法的合理性。 -
气垫船垫升状态的载荷计算相比排水状态要复杂,船体弯矩可以通过以下2种方式计算:
1) 根据目标实船,按照一定的缩尺比制作模型,并在水池中设定与实际海况相当的试验波浪,通过模型试验获得船体弯矩。全垫升气垫船对模型制作水平、风机垫升能力等的要求较高,虽然模型试验是获得数据最直接的方法,但是难度和成本相对较高。
2) 根据文献[13-14],基于弗劳德–克雷洛夫假设,结合垫升原理频域数值计算方法,对垫升航行状态下气垫船的载荷进行计算。
首先,通过气垫船的主尺度等数据信息建立目标船的几何模型,并按照质量分布建立质量模型;然后,应用全垫升理论,适当考虑风机特性、围裙触水及动态响应,以及气垫压缩性的影响,将船体-围裙-波面三者间的气室气体看成一个具有一定体积与压强的气垫,在单位规则波的波动作用下,其将形成气室内体积与压强的波动,波动的气垫压强会诱导船体形成各自由度的加速度响应和运动响应,据此得到船体在波浪中运动的数学模型,在求解运动的基础上对载荷进行计算。
-
按照前文,得到垫升工况及排水工况的梁载荷后,再根据式(9),将静水(静垫升)弯矩与波浪(垫升)弯矩相结合,得到船体梁载荷引起的名义应力σnh[6]为
$$ {\sigma _{{\rm{nh}},ij(k)}} = {\sigma _{{\rm{sw}},(k)}} + {C_{{\rm{WV}},ij}}{\sigma _{{\rm{WV}},ij(k)}} + {C_{{\rm{WH}},ij}}{\sigma _{{\rm{WH}},ij(k)}} $$ (9) 式中:CWV,CWH为载荷组合因子,其具体的取值见表1;σsw,(k)为静水弯矩引起的船体梁弯曲正应力,其中k的取值为1或2分别代表垫升及排水工况;σWV,ij(k),σWH,ij(k)分别为垂向弯矩和水平弯矩引起的船体梁弯曲应力;i=1,2,3,分别代表表1中的迎浪、斜浪和横浪;j=1,2,分别代表中垂或中拱。
-
装载甲板装载车辆所产生的局部装载压力
${P_{\rm{c}}}$ 按式(10)计算:$${P_{\rm{c}}} = {P_{{\rm{cs}}}} + {P_{{\rm{cw}}}}$$ (10) 式中:
$ {P}_{\mathrm{c}\mathrm{s}} $ 为静水中的装载压力,根据所装载武器装备的重量和轮印面积等进行计算;$ {P}_{\mathrm{c}\mathrm{w}} $ 为由气垫船在高速航行过程中所产生的垂向运动所引起的惯性力。若装载履带式车辆,则其载荷面积取为履带着地面积;若装载轮式车辆,结合文献[15],如果板格上2个轮印之间的距离很小,可将其视为一个轮印,则其单侧载荷面积为a×(2b+e),其中e为同一侧两轮的间距,a,b分别为轮印的长度及宽度。
$ {P}_{\mathrm{c}\mathrm{w}} $ 按式(11)[6]计算:$${p_{{\rm{cw}}}} = 0.65M{a_{\textit{z}}}$$ (11) 式中:0.65为结合CCS疲劳指南所确定的、升沉加速度对垂向加速度的贡献度;
$ {M}$ 为所装载武器装备的重量;$ {a}_{\textit{z}} $ 为该处的垂向加速度(考虑方向)。同时,考虑携载系数
$ f $ ,具体计算如式(12)[6]所示:$$f = 1 + 0.5\frac{{{a_{\textit{z}}}}}{g}$$ (12) -
液体装载的计算与常规排水型船舶基本类似。由液体装载引起的总压力PL按式(13)计算,取值不小于0。
$${p_{\rm{L}}} = {p_{{\rm{Ls}}}} + {p_{{\rm{Lw}}}}$$ (13) 式中:
$ {P}_{\mathrm{L}\mathrm{s}} $ 为由液体货物引起的静压力,按式(14)[6]计算;$ {P}_{\mathrm{L}\mathrm{w}} $ 为由液体货物引起的惯性压力,按式(15)[6]计算。$${p_{{\rm{Ls}}}} = {\rho _{\rm{L}}}g({{\textit{z}}_{{\rm{TOP}}}} - {\textit{z}})$$ (14) $$ {p_{{\rm{Lw}}}} = {\rho _{\rm{L}}}{\left[ {{f_{{\rm{ull}} - x}}{a_x}({x_B} - x) + {f_{{\rm{ull}} - y}}{a_y}({y_B} - y) + {a_{\textit{z}}}({{\textit{z}}_B} - {\textit{z}})} \right]^{\left[ 6 \right]}} $$ (15) 式中:
$ {\rho }_{\mathrm{L}} $ 为液体货物密度,t/m3,对于货油、LNG和压载水,分别取为0.9,0.5和1.025;${\textit{z}}_{\mathrm{T}\mathrm{O}\mathrm{P}}$ 为船舶正浮状态下液舱顶的垂向坐标;z为载荷点的垂向坐标;$ {f}_{\mathrm{u}\mathrm{l}\mathrm{l}-x} $ 为纵向充装系数;$ {f}_{\mathrm{u}\mathrm{l}\mathrm{l}-y} $ 为横向充装系数;其他参数以及其计算见中国船级社2018版疲劳指南。 -
进行强度校核时,气垫船在工作状态下气道中产生的局部压力不可忽略,其数值大体上可以由气垫压力通过囊压比计算求得。而气垫船在正常航行状态中垫升风机的效率一般保持在90%以上,故可近似将气垫所提供的较为持续稳定的支持力视为均布压力
$ {P}_{\mathrm{f}} $ 。气道压力$ {P}_{\mathrm{a}} $ 的计算如式(16)所示:$${p_{\rm{a}}} = \frac{{{p_{\rm{f}}}}}{\lambda } = \frac{{Tg}}{{s\lambda }}$$ (16) 式中:T为气垫船主船体重量,t;
$ s $ 为气垫接触面积,$ {\mathrm{m}}^{2} $ ;$ \lambda $ 为囊压比,通常情况下取1.39。 -
1) 在排水航行状态下,船底板受外部水动压力,其计算可参考低速航行的排水型船舶,依照CCS疲劳指南进行计算
2) 在垫升航行状态下,船底板受气垫中气体产生的向上的局部压力,该压力应与船体总重量相当,根据牛顿第二定律,由垂向加速度产生的垂向惯性力、重力以及向上的垫升力三者平衡,据此得到垫升力即可获得外部垫升压力。
由局部车辆搭载、液体装载引起的载荷,以及对疲劳的影响也可根据有限元方法进行计算,然后,可通过2.3及2.4节的计算得到船体梁应力及局部应力,考虑应力集中系数修正、腐蚀系数修正等所得到的热点应力,据此得到热点应力范围以及设计应力范围[2]。
$$ {S_{\rm{D}}}_{(k)}{\rm{ = }}\max ({f_{\rm{m}}}{S_{{\rm{h}},{i_{(k)}}}}) $$ (17) 式中:SD(k)为设计应力范围;Sh,i(k)为热点应力范围,即载荷工况“成对”工况的差值;fm为平均应力修正系数,参考CCS疲劳指南,可按式(18)和式(19)[2]进行计算。
焊接节点:
$$ {f_{{\rm{m}},i(k)}} = \left\{ \begin{aligned} & \min \left[1.0,\;0.85 + 0.3\frac{{{\sigma _{_{{\rm{m}},i(k)}}}}}{{{C_{\rm{s}}}{S_{_{{\rm{h}},i(k)}}}}}\right]{\rm{ , }}\;\;{\sigma _{_{{\rm{m}},i(k)}}} > 0 \\& \max \left[0.7,\;0.85 + 0.3\frac{{{\sigma _{_{{\rm{m}},i(k)}}}}}{{{C_{\rm{s}}}{S_{_{{\rm{h}},i(k)}}}}}\right]{\rm{ , }}\;\;{\sigma _{_{{\rm{m}},i(k)}}} < 0 \end{aligned} \right. $$ (18) 母材自由边:
$${f_{{\rm{m}},i(k)}} = \left\{ \begin{aligned} & \min \left[1.0,\; 0.85 + 0.4\frac{{{\sigma _{_{{\rm{m}},i(k)}}}}}{{{C_{\rm{s}}}{S_{_{{\rm{h}},i(k)}}}}}\right]{\rm{ , }}\;\;{\sigma _{_{{\rm{m}},i(k)}}} > 0 \\& \max \left[0.6,\; 0.8 + 0.4\frac{{{\sigma _{_{{\rm{m}},i(k)}}}}}{{{C_{\rm{s}}}{S_{_{{\rm{h}},i(k)}}}}}\right]{\rm{ , }}\;\;{\sigma _{_{{\rm{m}},i(k)}}} < 0 \end{aligned} \right.$$ (19) 式中:σm,i(k)为热点平均应力;Sh,j(k)为热点应力范围;CS为系数,按1.6+0.002 5L计算。
-
高速气垫船的材质一般为铝合金,结合规范及试验所得的铝合金试件S-N曲线,确定S-N曲线的参数如表2所示。表中:m为S-N曲线反斜率,取为3;K为S-N曲线参数;Sq为S-N曲线两线段交点处的应力幅值。
表 2 铝合金S-N曲线参数
Table 2. Parameters of S-N curve of aluminum alloy
名称 m K Sq 母材自由边 3 1.28×1011 29.5 焊接节点 3 2.5×1011 36.8 -
结构节点在载荷工况为k时的累积损伤度
$ {D}_{k} $ 按式(20)计算。由线性累积损伤理论,将所有工况下的损伤度相加即为最后的损伤度D。$$D = \sum {{D_k} = \sum {\frac{{{N_{\rm{D}}}{\alpha _k}}}{K}\frac{{{S_{{\rm{D}}(k)}^m}}}{{{{(\ln {N_{\rm{L}}})}^{m/\xi }}}}} } \Gamma \left(1 + \frac{m}{\xi }\right)$$ (20) 式中:
$ {N}_{\mathrm{D}}$ 为运营期内(15年)的载荷循环总次数,这里取为0.75×108;$ {N}_{\mathrm{L}}$ 为载荷谱回复周期的循环次数,取为102;$ {\alpha }_{k}$ 为工况为k的时间分配系数,根据气垫船的实际工作状态及各自所占时间比例确定,具体见第2节的时间分配;ξ为Weibull分布形状参数,取为1;Γ为完全Gamma函数值,可查表获得。船舶结构在设计寿命期间,根据线性累积原理,其总累积损伤度可由各工况下的损伤度相加求得[16]。符合疲劳强度要求的船舶结构应满足D≤1.0,且疲劳寿命按
$ {T}_{\mathrm{d}} $ /D即可计算求得。此处,$ {T}_{\mathrm{d}} $ 疲劳设计寿命,根据本船特点取为15年。 -
根据载荷特点及实际情况,可知船肿剖面附近为全船比较危险的地方。因21#剖面位置位于船肿,且船体底部为燃油舱,并搭载有车辆,会产生局部的装载压力,故选取21#横剖面的节点进行实船的全船疲劳强度简化计算分析。剖面图及节点选取、编号如图8。
按照2.3节所述,船体弯矩计算如表3和表4所示。将波浪弯矩乘以概率水平修正因子(0.25)及航区系数(0.9)即可得到修正后的数值。
表 3 排水航行状态载荷
Table 3. Load of drained navigation state
载荷类型 数值 修正后数值 备注 静水弯矩/(kN·m) 4 400 4 400 中垂状态 垂向波浪弯矩/(kN·m) 17 488 3 934.8 中垂状态 垂向波浪弯矩/(kN·m) 20 629 4 641.5 中拱状态 水平波浪弯矩/(kN·m) 1 430 321.75 – 表 4 垫升航行状态载荷
Table 4. Load of cushion navigation state
载荷类型 数值 修正后数值 备注 静垫升弯矩/(kN·m) 5 430 5 430 中垂状态 垂向垫升弯矩/(kN·m) 19 760 4 446 中垂状态 垂向垫升弯矩/(kN·m) 17 662 3 973.95 中拱状态 根据前文所述,21#剖面的疲劳强度简化计算如表5~表7所示。
表 5 船体梁载荷引起的名义应力
Table 5. Nominal stress caused by hull girder load
纵骨节点
编号垫升状态 排水状态 迎浪1 迎浪2 迎浪1 迎浪2 随浪1 随浪2 横浪1 横浪2 1 −12.329 −1.818 −10.405 0.302 −6.492 −4.318 −5.903 −4.995 2 −12.329 −1.818 −10.405 0.302 −6.513 −4.297 −5.797 −4.914 3 −12.329 −1.818 −10.405 0.302 −6.524 −4.285 −5.740 −4.914 4 −11.743 −1.731 −9.910 0.287 −6.228 −4.068 −5.398 −4.680 5 −8.268 −1.219 −6.978 0.202 −4.399 −2.849 −3.728 −3.295 6 −8.268 −1.219 −6.978 0.202 −4.383 −2.866 −3.809 −3.295 7 −3.033 −0.447 −2.560 0.074 −1.642 −1.017 −1.228 −1.209 8 −1.709 −0.252 −1.442 0.042 −0.981 −0.517 −0.411 −0.681 9 3.544 0.523 2.991 −0.087 1.791 1.316 2.072 1.412 10 13.206 1.947 11.145 −0.323 6.871 4.707 6.735 5.263 11 10.999 1.622 9.282 −0.269 5.745 3.898 5.497 4.383 12 15.413 2.272 13.008 −0.377 8.063 5.450 7.640 6.143 13 19.096 2.815 16.116 −0.467 9.976 6.767 9.537 7.610 14 20.364 3.002 17.186 −0.498 10.656 7.197 10.079 8.115 15 20.458 3.016 17.266 −0.500 10.745 7.192 9.930 8.153 16 −2.245 −0.331 −1.895 0.055 −1.204 −0.765 −0.967 −0.895 17 22.035 3.249 18.596 −0.539 11.558 7.761 10.768 8.782 18 21.997 3.243 18.564 −0.538 11.525 7.760 10.815 8.766 19 −2.245 −0.331 −1.895 0.055 −1.192 −0.776 −1.024 −0.895 20 −2.245 −0.331 −1.895 0.055 −1.179 −0.789 −1.090 −0.895 表 6 不同载荷工况下热点应力
Table 6. Hot spot stress under different load cases
纵骨节点
编号垫升状态 排水状态 迎浪1 迎浪2 迎浪1 迎浪2 随浪1 随浪2 横浪1 横浪2 1 −54.755 −38.243 −22.458 −6.629 −16.672 −13.458 −15.801 −14.459 2 −54.755 −38.243 −22.458 −6.629 −16.703 −13.427 −15.645 −14.339 3 −54.755 −38.243 −22.458 −6.629 −16.720 −13.410 −15.560 −14.339 4 −53.834 −38.108 −21.726 −6.650 −16.282 −13.088 −15.055 −13.993 5 −24.200 −13.128 −10.316 0.299 −6.504 −4.212 −5.511 −4.871 6 −24.200 −13.128 −10.316 0.299 −6.480 −4.237 −5.631 −4.871 7 5.322 9.325 −3.371 0.098 −2.162 −1.340 −1.617 −1.592 8 8.464 10.616 −1.813 0.053 −1.233 −0.650 −0.517 −0.856 9 15.253 10.789 3.759 −0.109 2.251 1.654 2.604 1.775 10 19.801 2.919 16.711 −0.484 10.302 7.058 10.098 7.891 11 16.491 2.431 13.918 −0.403 8.614 5.845 8.243 6.572 12 23.110 3.407 19.504 −0.565 12.090 8.172 11.455 9.210 13 27.269 4.020 19.291 −0.559 11.941 8.100 11.416 9.110 14 29.079 4.287 20.571 −0.596 12.755 8.615 12.064 9.714 15 29.214 4.307 20.667 −0.599 12.862 8.608 11.886 9.759 16 −800.714 −752.320 −789.639 −761.085 −788.553 −762.372 −788.182 −762.577 17 31.466 4.639 22.260 −0.645 13.835 9.289 12.890 10.511 18 31.412 4.631 22.222 −0.644 13.796 9.289 12.946 10.493 19 −800.714 −752.320 −789.639 −761.085 −788.535 −762.390 −788.272 −762.577 20 −800.714 −752.320 −789.639 −761.085 −788.515 −762.411 −788.374 −762.577 表 7 设计应力范围、损伤度及疲劳寿命
Table 7. Design stress range,damage degree and fatigue life
纵骨节点编号 属性 设计应力范围
(排水航行状态)/(kN·m)损伤度
(排水航行状态)设计应力范围
(垫升航行状态)/(kN·m)损伤度
(垫升航行状态)疲劳寿命/年 1 船底纵桁(燃油舱) 15.83 0.009 99 11.56 0.035 02 333.3 2 船底纵桁(燃油舱) 15.83 0.009 99 11.56 0.035 02 333.3 3 船底纵桁(燃油舱) 15.83 0.009 99 11.56 0.035 02 333.3 4 船底纵桁(燃油舱) 15.08 0.008 63 11.01 0.030 25 385.7 5 内底板纵桁 10.61 0.003 01 7.75 0.010 56 >1 000 6 内底板纵桁 10.61 0.003 01 7.75 0.010 56 >1 000 7 浮箱甲板纵桁 3.47 0.000 11 4.00 0.001 45 >1 000 8 气道纵桁 1.87 0.000 02 2.15 0.000 23 >1 000 9 气道甲板纵桁 3.87 0.000 15 4.46 0.002 02 >1 000 10 舷侧纵骨 17.20 0.012 81 16.31 0.098 31 135.0 11 纵舱壁纵桁 14.32 0.007 40 13.58 0.056 79 233.7 12 纵舱壁纵骨 20.07 0.020 37 19.03 0.156 30 84.9 13 甲板纵桁与横梁相交节点 19.85 0.019 71 22.46 0.256 77 54.3 14 甲板纵桁与横梁相交节点 21.17 0.023 90 23.95 0.311 38 44.7 15 甲板纵桁与横梁相交节点 21.27 0.024 23 24.06 0.315 74 44.1 16 装载甲板纵骨与横梁相交节点 28.55 0.058 66 33.88 0.881 59 16.0 17 甲板纵骨与横梁交接点 22.90 0.030 28 25.91 0.394 50 35.3 18 甲板纵骨与横梁交接点 22.87 0.030 12 25.87 0.392 47 35.5 19 装载甲板纵骨与横梁相交节点 28.55 0.058 66 33.88 0.881 59 16.0 20 装载甲板纵桁与横梁相交节点 28.55 0.058 66 33.88 0.881 59 16.0 由上述计算可以看出,目标气垫船21#剖面所选纵骨节点的疲劳寿命均符合疲劳设计寿命的要求,即疲劳寿命均大于15年。其中,装载甲板纵骨和甲板纵骨的寿命与其他纵骨节点相比偏低,其原因在于:对装载甲板纵骨而言,因搭载着较重的货物,会引起较大的惯性力变化,从而使得局部应力引起的应力范围变大,故而寿命较低;对甲板纵骨而言,因其位置远离中和轴,受总应力影响较大,故而寿命偏低。而同样远离中和轴的船底板纵骨,由于其尺寸较大,自身结构强度高于甲板纵骨,故而其疲劳寿命相对较高。
-
在国际上现行的主流规范中,关于船体疲劳强度评估方法,尤其是简化计算方法的主要研究对象是排水型常规船舶,针对气垫船尚无明确的疲劳强度评估方法。本文从参数选取、载荷计算、局部载荷、损伤度计算等方面讨论了气垫船疲劳强度评估的简化方法,并采用该方法,以实船为例进行了疲劳强度简化计算,主要得到以下结论:
1) 气垫船的运行方式与普通排水型船舶相比极其特殊,且其典型航行工况(如垫升航行状态和排水航行状态)的航行状态存在较大差异,导致其疲劳载荷的确定较为复杂,因此,气垫船疲劳载荷的输入是疲劳计算的重点和难点。
2) 气垫船尺寸较小,与普通排水型船舶相比,其在同等海况下发生的振动,即载荷循环的次数会更多,由搭载货物引起的局部应力对疲劳的影响也更严重,这些对于疲劳强度来说都是十分不利的,且计算结果也表明,对于车辆搭载甲板支撑纵骨及远离中和轴位置的纵骨节点,应对其疲劳问题给予足够的重视。这也证明本文所讨论的疲劳强度简化计算方法对于气垫船主艇体的初步疲劳评估来说是适用的。
3) 结合实际发生疲劳破坏的位置(例如发动机短舱、舷侧气道、气道甲板螺栓孔等)可以看出,简化计算对于局部位置的计算效果一般,考虑气垫实际状态、气道内气体分布等因素的疲劳强度直接计算方法还有待进一步的研究。
A simplified method for structural fatigue analysis of hovercraft
-
摘要:
目的 气垫船的运行原理不同于一般船型,其受力特殊,速度较一般舰船高出很多,目前国际上各船级社所提出的疲劳强度简化计算方法并不完全适用于气垫船,需要提出一种适用于气垫船疲劳强度校核的简化计算方法 方法 针对气垫船结构的疲劳强度问题,根据气垫船特殊的结构及运营模式等,结合主流规范、试验及相关研究,从工况确定、载荷计算、累计损伤推导等方面进行分析,提出一种适用于气垫船结构疲劳分析的简化方法,并以某气垫船为例进行算例分析。 结果 简化计算结果显示典型纵骨节点的疲劳寿命均符合疲劳强度要求,算例分析结果表明该方法具有一定的适用性。 结论 所提简化方法对气垫船的研发设计及校核具有参考意义。 Abstract:Objective The operation principle of a hovercraft is different from those of general ship types, with special stresses and much higher speeds than general ships. At present, the simplified fatigue strength calculation method proposed by international classification societies is not fully applicable to hovercraft, so it is necessary to put forward a suitable simplified calculation method. Methods According to the special structure and operation mode of hovercraft, and combined with mainstream specifications, experiments and related research, this paper analyzes the working conditions, load calculations, cumulative damage derivation and other aspects, and puts forward a simplified method suitable for the fatigue analysis of hovercraft structures using a specific hovercraft as an example for calculation. Results The simplified calculation results show that the fatigue life of typical longitudinal nodes meets the requirements of fatigue strength, and the analysis results of the example show that this method has certain applicability. Conclusions The simplified method proposed in this paper can be used as a reference for the development, design and verification of hovercraft. -
Key words:
- hovercraft /
- fatigue strength /
- simplified calculation /
- specifications
-
表 1 各工况下各载荷分量的分配系数
Table 1. Distribution coefficients of each load component under different working conditions
浪向 工况 静水(垫升)弯矩 中垂弯矩 中拱弯矩 水平弯矩 迎浪 1 1 1 0 0 2 1 0 −1 0 斜浪 3 1 0.2 0 −0.2 4 1 0 −0.2 0.2 横浪 5 1 0.1 0 −1 6 1 0 −0.1 1 表 2 铝合金S-N曲线参数
Table 2. Parameters of S-N curve of aluminum alloy
名称 m K Sq 母材自由边 3 1.28×1011 29.5 焊接节点 3 2.5×1011 36.8 表 3 排水航行状态载荷
Table 3. Load of drained navigation state
载荷类型 数值 修正后数值 备注 静水弯矩/(kN·m) 4 400 4 400 中垂状态 垂向波浪弯矩/(kN·m) 17 488 3 934.8 中垂状态 垂向波浪弯矩/(kN·m) 20 629 4 641.5 中拱状态 水平波浪弯矩/(kN·m) 1 430 321.75 – 表 4 垫升航行状态载荷
Table 4. Load of cushion navigation state
载荷类型 数值 修正后数值 备注 静垫升弯矩/(kN·m) 5 430 5 430 中垂状态 垂向垫升弯矩/(kN·m) 19 760 4 446 中垂状态 垂向垫升弯矩/(kN·m) 17 662 3 973.95 中拱状态 表 5 船体梁载荷引起的名义应力
Table 5. Nominal stress caused by hull girder load
纵骨节点
编号垫升状态 排水状态 迎浪1 迎浪2 迎浪1 迎浪2 随浪1 随浪2 横浪1 横浪2 1 −12.329 −1.818 −10.405 0.302 −6.492 −4.318 −5.903 −4.995 2 −12.329 −1.818 −10.405 0.302 −6.513 −4.297 −5.797 −4.914 3 −12.329 −1.818 −10.405 0.302 −6.524 −4.285 −5.740 −4.914 4 −11.743 −1.731 −9.910 0.287 −6.228 −4.068 −5.398 −4.680 5 −8.268 −1.219 −6.978 0.202 −4.399 −2.849 −3.728 −3.295 6 −8.268 −1.219 −6.978 0.202 −4.383 −2.866 −3.809 −3.295 7 −3.033 −0.447 −2.560 0.074 −1.642 −1.017 −1.228 −1.209 8 −1.709 −0.252 −1.442 0.042 −0.981 −0.517 −0.411 −0.681 9 3.544 0.523 2.991 −0.087 1.791 1.316 2.072 1.412 10 13.206 1.947 11.145 −0.323 6.871 4.707 6.735 5.263 11 10.999 1.622 9.282 −0.269 5.745 3.898 5.497 4.383 12 15.413 2.272 13.008 −0.377 8.063 5.450 7.640 6.143 13 19.096 2.815 16.116 −0.467 9.976 6.767 9.537 7.610 14 20.364 3.002 17.186 −0.498 10.656 7.197 10.079 8.115 15 20.458 3.016 17.266 −0.500 10.745 7.192 9.930 8.153 16 −2.245 −0.331 −1.895 0.055 −1.204 −0.765 −0.967 −0.895 17 22.035 3.249 18.596 −0.539 11.558 7.761 10.768 8.782 18 21.997 3.243 18.564 −0.538 11.525 7.760 10.815 8.766 19 −2.245 −0.331 −1.895 0.055 −1.192 −0.776 −1.024 −0.895 20 −2.245 −0.331 −1.895 0.055 −1.179 −0.789 −1.090 −0.895 表 6 不同载荷工况下热点应力
Table 6. Hot spot stress under different load cases
纵骨节点
编号垫升状态 排水状态 迎浪1 迎浪2 迎浪1 迎浪2 随浪1 随浪2 横浪1 横浪2 1 −54.755 −38.243 −22.458 −6.629 −16.672 −13.458 −15.801 −14.459 2 −54.755 −38.243 −22.458 −6.629 −16.703 −13.427 −15.645 −14.339 3 −54.755 −38.243 −22.458 −6.629 −16.720 −13.410 −15.560 −14.339 4 −53.834 −38.108 −21.726 −6.650 −16.282 −13.088 −15.055 −13.993 5 −24.200 −13.128 −10.316 0.299 −6.504 −4.212 −5.511 −4.871 6 −24.200 −13.128 −10.316 0.299 −6.480 −4.237 −5.631 −4.871 7 5.322 9.325 −3.371 0.098 −2.162 −1.340 −1.617 −1.592 8 8.464 10.616 −1.813 0.053 −1.233 −0.650 −0.517 −0.856 9 15.253 10.789 3.759 −0.109 2.251 1.654 2.604 1.775 10 19.801 2.919 16.711 −0.484 10.302 7.058 10.098 7.891 11 16.491 2.431 13.918 −0.403 8.614 5.845 8.243 6.572 12 23.110 3.407 19.504 −0.565 12.090 8.172 11.455 9.210 13 27.269 4.020 19.291 −0.559 11.941 8.100 11.416 9.110 14 29.079 4.287 20.571 −0.596 12.755 8.615 12.064 9.714 15 29.214 4.307 20.667 −0.599 12.862 8.608 11.886 9.759 16 −800.714 −752.320 −789.639 −761.085 −788.553 −762.372 −788.182 −762.577 17 31.466 4.639 22.260 −0.645 13.835 9.289 12.890 10.511 18 31.412 4.631 22.222 −0.644 13.796 9.289 12.946 10.493 19 −800.714 −752.320 −789.639 −761.085 −788.535 −762.390 −788.272 −762.577 20 −800.714 −752.320 −789.639 −761.085 −788.515 −762.411 −788.374 −762.577 表 7 设计应力范围、损伤度及疲劳寿命
Table 7. Design stress range,damage degree and fatigue life
纵骨节点编号 属性 设计应力范围
(排水航行状态)/(kN·m)损伤度
(排水航行状态)设计应力范围
(垫升航行状态)/(kN·m)损伤度
(垫升航行状态)疲劳寿命/年 1 船底纵桁(燃油舱) 15.83 0.009 99 11.56 0.035 02 333.3 2 船底纵桁(燃油舱) 15.83 0.009 99 11.56 0.035 02 333.3 3 船底纵桁(燃油舱) 15.83 0.009 99 11.56 0.035 02 333.3 4 船底纵桁(燃油舱) 15.08 0.008 63 11.01 0.030 25 385.7 5 内底板纵桁 10.61 0.003 01 7.75 0.010 56 >1 000 6 内底板纵桁 10.61 0.003 01 7.75 0.010 56 >1 000 7 浮箱甲板纵桁 3.47 0.000 11 4.00 0.001 45 >1 000 8 气道纵桁 1.87 0.000 02 2.15 0.000 23 >1 000 9 气道甲板纵桁 3.87 0.000 15 4.46 0.002 02 >1 000 10 舷侧纵骨 17.20 0.012 81 16.31 0.098 31 135.0 11 纵舱壁纵桁 14.32 0.007 40 13.58 0.056 79 233.7 12 纵舱壁纵骨 20.07 0.020 37 19.03 0.156 30 84.9 13 甲板纵桁与横梁相交节点 19.85 0.019 71 22.46 0.256 77 54.3 14 甲板纵桁与横梁相交节点 21.17 0.023 90 23.95 0.311 38 44.7 15 甲板纵桁与横梁相交节点 21.27 0.024 23 24.06 0.315 74 44.1 16 装载甲板纵骨与横梁相交节点 28.55 0.058 66 33.88 0.881 59 16.0 17 甲板纵骨与横梁交接点 22.90 0.030 28 25.91 0.394 50 35.3 18 甲板纵骨与横梁交接点 22.87 0.030 12 25.87 0.392 47 35.5 19 装载甲板纵骨与横梁相交节点 28.55 0.058 66 33.88 0.881 59 16.0 20 装载甲板纵桁与横梁相交节点 28.55 0.058 66 33.88 0.881 59 16.0 -
[1] 中国船级社. 海上高速船入级及建造规范[S]. 北京: 人民交通出版社, 2015. China Classification Society. Rules for construction and classification of sea-going high speed craft[S]. Beijing: People's Communications Press, 2015 (in Chinese). [2] 李建彰. 全垫升气垫船结构强度评估及局部振动分析[D]. 哈尔滨: 哈尔滨工程大学, 2014. LI J Z. Structure strength assessment and local vibration analysis of air cushion vehicle[D]. Harbin: Harbin Engineering University, 2014 (in Chinese). [3] 徐伟. 全垫升式气垫船总振动计算方法研究[D]. 大连: 大连理工大学, 2018. XU W. The hull vibration calculation method research of the air cushion craft[D]. Dalian: Dalian University of Technology, 2018 (in Chinese). [4] 李建彰, 任慧龙, 刘宁, 等. 砰击载荷作用下气垫船船底结构瞬态响应分析[J]. 舰船科学技术, 2014, 36(3): 24–28. LI J Z, REN H L, LIU N, et al. Transient response analysis on bottom structure of the air cushion vehicle under slamming load[J]. Ship Science and Technology, 2014, 36(3): 24–28 (in Chinese). [5] 曲鹏. 小型气垫船结构强度有限元分析研究[D]. 哈尔滨: 哈尔滨工程大学, 2013. QU P. Finite element analysis of structural strength of small hovercraft[D]. Harbin: Harbin Engineering University, 2013 (in Chinese). [6] 中国船级社. 船体结构疲劳强度指南[S]. 北京: 人民交通出版社, 2018. China Classification Society. Guidelines for fatigue strength of hull structures[S]. Beijing: People's Communications Press, 2018 (in Chinese). [7] ABS. Guide for building and classing high speed naval craft[S]. Houston: American Bureau of Shipping, 2016. [8] DNV. Fatigue assessment of ship structure[S]. Hovik: DNV, 2003. [9] 中国人民解放军总装备部. 水面舰艇结构设计计算方法: GJB/Z 119-99[S]. [S.l.: s.n.], 1999. General Armament Department of PLAN. Calculation method of surface ship structure design: GJB/Z 119-99[S]. [S.l.: s.n.], 1999 (in Chinese). [10] 刘宁, 唐首祺, 任慧龙. 基于国内外部分规范的疲劳强度简化计算对比研究[C]//2019年船舶结构力学学术会议论文集. 武汉: 船舶力学学术委员会, 2019: 153-158. LIU N, TANG S Q, REN H L. Comparative study on simplified calculation of fatigue strength based on some domestic and foreign codes[C]//Proceedings of the Conference on Ship Structural Mechanics in 2019. Wuhan: Academic Committee on Ship Mechanics, 2019: 153-158 (in Chinese). [11] 李桂春, 李铮捷. 气垫船气囊与气垫[M]. 北京: 国防工业出版社, 2011: 48-64. LI G C, LI Z J. Air bag and air cushion of hovercraft[M]. Beijing: National Defense Industry Press, 2011: 48-64 (in Chinese). [12] YANG P, GU X K. Discussion of some important parameters in fatigue loading calculation for ship structural[J]. Journal of Ship Mechanics, 2012, 16(6): 658–667. [13] 刘宁, 王晓强, 任慧龙, 等. 全垫升气垫船在波浪中运动的非线性理论研究[J]. 华中科技大学学报(自然科学版), 2014, 42(4): 91–95. LIU N, WANG X Q, REN H L, et al. Nonlinear theory research on motion response of air cushion vehicle in waves[J]. Journal of Huazhong University of Science and Technology (Nature Science Edition), 2014, 42(4): 91–95 (in Chinese). [14] 唐尧. 全垫升气垫船载荷及运动响应计算方法研究[D]. 哈尔滨: 哈尔滨工程大学, 2011. TANG Y. Computational method of motions and wave loads on air cushion vehicle[D]. Harbin: Harbin Engineering University, 2011 (in Chinese). [15] 邱吉廷, 吴嘉蒙. 汽车运输船车辆甲板结构构件的设计研究[C]//纪念徐秉汉院士船舶与海洋结构力学学术会议论文集. 无锡: 中国造船工程学会, 2011: 410-422. QIU J T, WU J M. Design and research on structural members of vehicle deck of motor transport ships[C]//Proceedings of the Academic Conference on Ship And Ocean Structural Mechanics In Memory of Academician Xu Binghan. Wuxi: China Society of Shipbuilding Engineering, 2011: 410-422 (in Chinese). [16] 甄春博, 冯亮, 任慧龙, 等. 一种三体船结构疲劳分析的简化方法[J]. 舰船科学技术, 2015, 37(11): 27–30. ZHEN C B, FENG L, REN H L, et al. One simplified method for fatigue assessment of trimaran structure[J]. Ship Science and Technology, 2015, 37(11): 27–30 (in Chinese). -