基于热点应力的超规范油船疲劳强度分析.pdf

新技术新工艺 2 0 2 3年 第9期 2 4 新技术新工艺 设计与计算基于热点应力的超规范油船疲劳强度分析沈海涛,郭哲璐(浙江浙能迈领环境科技有限公司,浙江 杭州 3 1 1 2 0 0)摘 要:为了方便清舱、减轻空船重量、提高货舱容积,对一艘7 5 0 0 DWT油船的舱壁结构做不设顶凳的超规范设计。基于J T P(双壳油船共同结构规范),运用热点应力法对该结构进行疲劳强度分析。应用有限元软件建立模型,并对热点位置精细化网格。采用数值计算模块计算不同工况下疲劳部位的热点应力。应用S-N曲线法计算出疲劳累积损伤并估算疲劳寿命,结果表明:该舱段结构疲劳寿命满足规范设计要求;在满足疲劳强度的前提下,舱壁无顶凳设计合理可行,并且具有一定的推广价值。关键词:油船;顶凳;热点应力;疲劳强度;有限元分析中图分类号:U 6 6 1.4 2 文献标志码:AA n a l y s i s o f F a t i g u e S t r e n g t h f o r S u p e r g a u g e T a n k e r s b a s e d o n H o t S p o t S t r e s sS HE N H a i t a o,GUO Z h e l u(Z h e j i a n g Z h e n e n g M a r i n e E n v i r o n m e n t a l T e c h n o l o g y C o.,L t d.,H a n g z h o u 3 1 1 2 0 0,C h i n a)A b s t r a c t:I n o r d e r t o c l e a n h o l d e a s i l y,l i g h t e n e m p t y s h i p a n d i n c r e a s e h o l d c a p a c i t y,a s u p e r g a u g e d e s i g n w a s p r o-p o s e d f o r a 7 5 0 0 DWT t a n k e r w i t h b u l k h e a d s o f n o n-u p p e r s t o o l.B a s e d o n J T P(J o i n t T a n k e r P r o j e t C o mm o n S t r u c t u r a l R u l e s),t h e f a t i g u e s t r e n g t h o f t h i s s t r u c t u r e w a s a n a l y z e d b y h o t s p o t s t r e s s m e t h o d.T h e m o d e l w a s s e t u p b y u s i n g f i n i t e e l e m e n t s o f t w a r e a n d t h e h o t s p o t p a r t m e s h w a s r e f i n e d.T h e h o t s p o t s t r e s s i n d i f f e r e n t w o r k i n g c o n d i t i o n s w a s c a l c u l a t e d t h r o u g h n u m e r i c a l c a l c u l a t i o n m o d u l e.T h e f a t i g u e c u m u l a t i v e d a m a g e w a s c a l c u l a t e d a n d t h e f a t i g u e l i f e w a s e s t i m a t e d b y t h e m e t h o d o f S-N c u r v e.T h e r e s u l t s h o w e d t h a t t h e f a t i g u e l i f e o f t h e c a r g o h o l d m e t t h e s t a n d a r d a n d d e s i g n r e q u i r e-m e n t s.T h e d e s i g n o f t h e b u l k h e a d s w i t h n o u p p e r s t o o l s w a s r e a s o n a b l e a n d f e a s i b l e,a n d i t h a d s o m e p r o m o t i o n a l v a l u e s.K e y w o r d s:t a n k e r s,u p p e r s t o o l,h o t s p o t s t r e s s,f a t i g u e s t r e n g t h,f i n i t e e l e m e n t a n a l y s i s 近几十年的造船资料1-3显示:典型的油船货油舱段在槽形舱壁上下端一般会设置顶凳和底凳以提高结构强度。然而顶凳和底凳的设置会有不少缺点,对此不少专家学者做了很多研究4-5。研究发现:槽形舱壁不设置顶凳不仅可以方便清舱,还可以减轻空船重量,提高货舱容积。这对于降低船舶建造和运营成本可以起到一定作用,因此舱壁无顶凳设计具有一定的现实意义。然而,槽形舱壁不设置顶凳会增加舱壁跨距,使应力水平提高,从而可能影响到结构的疲劳寿命6。因此对于无顶凳设计的舱段结构进行疲劳强度分析很有必要。本文以一艘7 5 0 0 DWT油船为例,根据 国内航行海船建造规范 进行结构设计。规范要求槽形舱壁必须设置顶凳。考虑增大舱容空间、方便清舱等因素,本船做无顶凳设计,因此属于超规范船舶。目前,针对此类超规范船舶的疲劳计算的研究较少,应用现有的国内规范无法直接校核。因此,为解决超规范油船疲劳计算问题,本文提出基于热点应力的疲劳强度评估方法。参照J T P(J o i n t T a n k e r P r o-j e t C o mm o n S t r u c t u r a l R u l e s)规范的相关要求,建立该船舶三舱段有限元模型,并对疲劳热点位置进行精细化网格划分,对船舶满载和压载等工况进行疲劳应力计算,对该船舱段结构进行疲劳强度评估。1 疲劳强度评估方法1.1 热点应力组成及确定热点应力法是针对结构局部热点的疲劳强度评估方法,该方法主要应用在焊接结构上,目前在船舶、海洋结构物、机械等领域都有一定的应用7。由于该方法考虑了构件的外形引起的热点应力集中,相比传统方法更适合工程实际。所谓热点应力就是指最大的结构应力或者说是结构最危险处的应力8。它由两部分组成:弯曲应力和膜应力,而这两者相加的最大值就是该结构的热点应力。有别于缺口应力,热点应力不包含因焊接缺口、裂纹等导致的局部应力集中产生的应力峰值,它只取决于焊缝接头的宏观尺寸和载荷大小9-1 0。由于焊缝表面与母材交界部位最容易产生疲劳裂纹,因此对于焊接结构,热点一般位于焊缝的焊趾处。对于如何确定结构的热点应力,不同规范有不同的规定。船体结构疲劳强度指南 规定:疲劳热点应力可由距离焊趾0.5t和1.5t(t为接头处板的净厚度)处的应力通过两点线性外推计算得出1 1;设计计算 新技术新工艺 设计与计算2 5 J T P规范规定:疲劳热点应力由距焊趾位置0.5t处的表面应力通过相邻两单元形心处的应力沿船宽方向线性内插计算得出1 2。本文参考J T P规范,使用第二种方法进行计算。根据J T P规范,边舱折角焊接接头部位的热应力就是垂直于焊接方向的单元应力。具体的热点应力位于边舱板和内底板中心线交线0.5t+xw t处(见图1)。图1 热点应力的位置示意图图1中,xw t为焊趾距离,内底板净厚度可由式1计算得到:t=tg r s-0.5tc o r r(1)式中,tg r s为总厚度;tc o r r为腐蚀增量。1.2 热应力计算根据J T P规范要求,须对油船满载和正常压载两种工况分别进行疲劳强度校核。当处于满载工况时,由货舱内部动水压力导致的应力范围Sw可由式2进行计算:Sw=0.4S1x+0.9S1y+0.9S1z(2)因此,对于油船满载工况,其热点应力范围S1可由式3进行计算:S1=fm0.8 5(Se1+0.2 5Se2)-0.3Sw(3)式中,fm为折减修正系数,由于强度计算时关键位置使用腐蚀增量0.5tc o r r,因此fm取0.9 5。对于油船正常压载工况,其热点应力范围可用式4求得:S2=fm0.8 5(Se1-0.2 5Se2)(4)为消除平均应力对计算结果的影响,应对两个工况下的热点应力范围进行修正,得到更为准确的总应力范围S R1和S R2:S R1=1.0S1(满载工况)(5)S R2=0.6S2(压载工况)(6)1.3 疲劳损伤及寿命评估为估算结构疲劳寿命,应先计算该结构的疲劳损伤。本文采用线性累积损伤模型进行计算,遵循P a l m g r e n-M i n e r规则,用DM值大小反映结构累积损伤程度,同时应满足DM1,否则可视为疲劳寿命无法满足设计要求。船舶结构的疲劳累积损伤DM可用式7进行计算:DM=2i=1DMi(7)式中,i=1或i=2分别代表满载工况和压载工况。假设长时间内结构应力范围遵循W e i b u l l概率分布,则对应不同工况下的疲劳损伤DMi可用式8进行计算:DMi=iNLK2SmR il nNR mi1+m (8)式中,NL为循环次数,与设计疲劳寿命对应,其大小可用式9进行计算:NL=f0U4 l o gL(9)式中,U为设计寿命,本船取2 5年;L为规范船长,为1 1 1.7 2 m;f0为考虑非航行状态的因子,取值0.8 5;i为两种载况下的船舶寿命系数,取值0.5;NR为对应概率为1 0-5的循环次数,取值1 0 0 0 0 0。为G a mm a函数;K2、m和i均为S-N曲线相关参数。为W e i b u l l概率分布参数,可用式1 0进行计算:=1.1-0.3 5L-1 0 03 0 0(1 0)参照J T P规范,舱段焊接结构的疲劳寿命可用疲劳累积损伤DM通过式1 1求得:疲劳寿命=设计寿命DM(1 1)为满足规范设计要求,通过计算求得的疲劳寿命应不小于设计寿命。本节中未提及的公式变量的具体含义详见参考文献1 3 中附录C所示。2 船体舱段模型2.1 船体结构及参数设计本文评估的是一艘7 5 0 0 DWT超规范油船,全船采用纵骨架式,货舱区域为双壳双底结构,出于方便清舱、提高舱容和减轻空船重量等因素考虑,本船槽形舱壁结构采用无顶凳设计(见图2和图3)。船舶具体的主尺度参数见表1。图2 纵舱壁结构图 新技术新工艺 2 0 2 3年 第9期 2 6 新技术新工艺 设计与计算图3 横舱壁结构图表1 船舶主尺度总长/m规范船长/m型深/m型宽/m吃水/m1 1 9.8 01 1 1.7 28.1 01 7.2 06.5 02.2 计算模型及网格设置本油船共设有5对货油舱,取中间2、3、4号3段货油舱作为计算对象,应用有限元软件建立模型(见图4)。模型纵向从F r 5 9肋位到F r 1 3 1+1 5 0肋位,横向采用全宽模式,垂向从船底到上甲板。图4 三舱段计算模型根据J T P规范要求,网格单元全部采用梁单元和板单元建模,按照骨架的实际排列规则划分网格。对于热点位置(底边舱和内底板交界处(见图5)进行精细化网格划分(见图6),网格尺寸为tt(t为内底板净厚度),并在热点所有方向延伸1 0个单元;此外,在热点附近5 0 0 mm范围内采用净厚度(即总厚度减去1/2的腐蚀增量0.5tc o r r)建模,其他区域建模厚度均采用总厚度扣除1/4的腐蚀增量(即0.2 5tc o r r)。图5 热点所处位置图6 精细化网格划分该油船舱段结构材料为“C C S-A”级钢,其材料属性参数如下:弹性模量E=2.0 61 05 N/mm2;泊松比=0.3 0;屈服强度s=2 3 5 MP a。2.3 边界条件本船舱段模型边界条件采用模型两端施加独立的弹簧单元方法设置,其中弹簧单元刚度c可用式1 2求得。具体边界条件设置见表2。c=E1+As-n e t5 0lt kn=0.7 7As-n e t5 0Elt kn(1 2)式中,变量含义详见参考文献1 1 中附录C。表2 端部边界约束位置平移转动xyzxyz后端后端(所有纵向单元)相关相关相关后端独立点约束Mv-e n dMh-e n d甲板,内底板和外板弹簧内壳、舷侧与纵舱壁弹簧前端前端(所有纵向单元)相关相关相关前端独立点Mv-e n dMh-e n d甲板,内底板和外板弹簧内壳、舷侧与纵舱壁弹簧 设计计算 新技术新工艺 设计与计算2 7 2.4 计算工况及载荷对于油船来说,在整个生命周期中主要存在满载航行、压载航行和停港或进坞修理3种状态,其中船舶停港或进坞修理时一般不会产生疲劳损伤,因此油船的疲劳累积损伤可认为是满载和压载这两种工况下的疲劳损伤之和。因为在这两种工况下船舶均处于航行状态,所以载荷计算只考虑动载荷。对于满载工况,需计算波浪载荷和舱内动水载荷;对于压载工况,由于船舶处于空舱状态,所以只需计算波浪载荷;同时舱段受到垂向和水平船体梁弯矩的影响,需对各工况下应力范围进行修正。具体工况设置见表3。其中各工况下的施加载荷的计算方法详见J T P规范。表3 计算工况工况应力分量施加载荷满载工况L1Se1 波浪动压力仅施加于被校核底边舱折角的船体一舷L2Se2 波浪动压力仅施加于未被校核底边舱折角的船体一舷L3Si x 纵向加速度引起的舱室内部动压力L4Si y 横向加速度引起的舱室内部动压力L5Si z 垂向加速度引起的舱室内部动压力压载工况L6Se1 波浪动压力仅施加于被校核底边舱折角的船体一舷L7Se2 波浪动压力仅施加于未被校核底边舱折角的船体一舷弯矩修正的载荷工况C1SV BM 单位垂向弯矩施加于舱段模型端部C2SHBM 单位水平弯矩施加于舱段模型端部3 计算结果及分析3.1 热点应力计算运用软件数值计算模块计算上述L1L7以及C1、C2共9种不同工况下的舱段结构的应力水平,进而确定疲劳热点位置处的应力大小。为使计算结果更加准确,考虑水平和垂向船体弯矩对热点应力的影响,可通过式1 3对初始数据进行修正(经过船体弯矩修正前后的计算结果汇总见表4):Sc_i=sc_i-MV_isV BM-MH_isHBM(1 3)式中,sc_i和Sc_i分别代表修正前后的热点应力分量;MV_i和MH_i分别代表垂向船体弯矩和水平船体弯矩;SV BM和SHBM分别代表C1和C2工况中计算得到的热点应力。表4 各工况下的热点应力范围工况应力分量修正前应力/MP a修正后应力/MP a满载工况L1Se11 1 6.31 1 4.7L2Se21 1.81 0.2L3Si x0.0 30.0 2L4Si y6 0.25 8.3L5Si z9 5.89 5.5压载工况L6Se18 1.28 0.1L7Se21 8.81 7.7弯矩修正的载荷工况C1SV BM0.41 0-90.41 0-9C2SHBM0.21 0-80.21 0-83.2 S-N曲线选取采用S-N曲线法对该油船结构进行疲劳强度分析,根据J T P规范要求,选取合适的S-N曲线(见图7)。其中S代表交变应力范围,N代表循环次数,这两者之间的函数关系可用式1 4表示:l o g(N)=l o g(K)-ml o g(S)(1 4)式中,m为曲线反斜率;K2为S-N曲线参数。图7 S-N曲线由于本油船疲劳热点部位用角焊缝焊接,焊缝周围均匀、致密、连续,而且无焊渣、焊瘤等表面缺陷,同时载荷施力方向与焊缝平行,因此选择D级S-N曲线进行疲劳损伤计算。该曲线的具体参数见表5。表5 S-N曲线参数S-N曲线类别mK2Sq/MP aD3.01.5 21 01 25 3.43.3 疲劳寿命评估疲劳强度计算遵循P a l m g r e n-M i n e r规则,根据上述方法进行油船结构疲劳寿命评估。根据式8计算求得满载工况和压载工况下的疲劳损伤:DM1=0.3 9 3;DM2=0.0 1 5。再根据式7计算求得该结构的疲劳累积损伤:DM=0.4 0 8。最后根据式1 1计算求得该油船结构的疲劳寿命:新技术新工艺 2 0 2 3年 第9期 2 8 新技术新工艺 设计与计算疲劳寿命=设计寿命DM=6 1.2 7 5(年)本船设计寿命为2 5年,因此该舱段结构设计符合规范要求。4 结语本文对一艘7 5 0 0 DWT油船做舱壁无顶凳的超规范设计。参照J T P规范,确定热点位置,运用有限元法计算出不同工况下的热点应力范围。应用S-N曲线法计算出热点位置的疲劳累积损伤并估算疲劳寿命,得到如下结论。1)底边舱和内底板交界处应力较为集中,当受到波浪载荷和舱内动水载荷的联合作用时,容易产生疲劳破坏,可对该部位做加强处理。2)对热点位置进行疲劳强度分析,计算得出该结构疲劳寿命约为6 1年,大于2 5年设计寿命,符合规范要求,从而验证舱壁无顶凳设计合理可行。3)槽形舱壁无顶凳设计具有方便清舱、减轻空船重量、提高货舱容积的优点,在满足疲劳强度前提下,该超规范设计具有一定的推广价值。本文结论可为同类船舶超规范设计和疲劳强度计算提供参考。参考文献1 马静敏,沈友徽,李华.基于i S I GHT的油船槽形横舱壁优化设计J.船海工程,2 0 0 5(1):3 5-3 7.2 许允,顾俊,胡斌.阿芙拉型油船节点优化设计J.船海工程,2 0 1 7,4 6(2):1 7-2 0.3 冯国庆,刘琳琳,任慧龙,等.HC S R疲劳评估中舱壁位移引起的纵骨端部附加应力研究J.船海工程,2 0 1 4,4(2):3 9-4 3.4 何皛磊,尹群.大型油船槽形舱壁结构形式设计及模糊综合评判J.船舶工程,2 0 1 0,3 2(1):1 5-1 9.5 甘水来,李国强,李勇.基于共同规范的散货船槽型横舱壁设计研究J.船舶与海洋工程,2 0 1 4(1):3 2-3 8.6 陈宝松,江晓俐.基于二次抛物线插值的船体疲劳强度评估J.船海工程,2 0 1 1,4 0(2):2 8-3 1,3 5.7 刘亮,李辉,周广喜,等.基于HC S R和C S R-O T的油船疲劳有限元法对比分析J.船海工程,2 0 1 3,4 2(6):3 8-4 1.8 刘新胜,李晓聪,杨丽娜,等.S MT焊点疲劳寿命的预估J.新技术新工艺,2 0 2 0(1 1):7 0-7 5.9 蒋文进,李良碧,罗广恩,等.F P S O典型节点疲劳寿命分析J.舰船科学技术,2 0 1 4,3 6(1 2):2 9-3 3,3 8.1 0 时宏森,杨涛,唐超,等.基于AN S Y S W o r k b e n c h的扭转弹簧疲劳寿命分析J.新技术新工艺,2 0 2 1(5):2 4-2 9.1 1 中国船级社.船体结构疲劳强度指南M.北京:人民交通出版社,2 0 1 5.1 2 I A C S,J o i n t t a n k e r p r o j e c t c o mm o n s t r u c t u r a l r u l e sS.1 3 I A C S,R u l e s J c o mm o n s t r u c t u r e r u l e s f o r o i l t a n k e r sS.作者简介:沈海涛(1 9 8 2-),男,高级工程师,硕士,主要从事船舶环保设备研发、有限元分析、大气污染治理等方面的工作。收稿日期:2 0 2 3-0 6-0 2责任编辑 郑练(上接第2 3页)送实验室办公室。3.开放基金课题的申报方式(1)申请者须按规定格式撰写 特种表面保护材料及应用技术国家重点实验室开放基金课题申请书(附件1),保证所有提交的申报材料的真实性,规定时间内向实验室提交签字盖章后的开放基金课题申请书电子文档两份(一份w o r d文档格式和一份加盖单位公章的P D F扫描件)、纸质申请书一式两份(加盖单位公章)。(2)基金课题申请书电子文档发送至实验室电子邮箱:q u f a n r i m p.c o m.c n;纸质原件邮寄至:湖北省武汉市硚口区宝丰二路1 2 6号,特种表面保护材料及应用技术国家重点实验室(邮编:4 3 0 0 3 0)。(3)本年度开放基金课题申请截止日期:2 0 2 3年1 0月1 3日。申请者如有较为突出的创新学术思想,或者需要开展与国家重大需求相关的创新研究工作,在与实验室联系后,可随时提出开放基金课题申请。(4)确定资助课题及资助金额后,将通知申请者填写任务书,未通过者不再另行通知。任务书由申请者所在单位签署审查意见并盖章后,将一式三份任务书邮寄至实验室联系人处,并按计划执行。4.联系方式联系人:瞿凡电子邮件:q u f a n r i m p.c o m.c n电话:[phone]通讯地址:湖北省武汉市硚口区宝丰二路1 2 6号(来源于特种表面保护材料及应用技术国家重点实验室)