插值方式對(duì)基于共同規(guī)范油船疲勞評(píng)估的影響

劉 俊 黃 銘

（上海交通大學(xué)船舶海洋與建筑工程學(xué)院 上海 200240）

0 引 言

2006年4月1日生效的《油船結(jié)構(gòu)共同規(guī)范（CSR）》是國(guó)際船級(jí)社協(xié)會(huì)（IACS）為了抵制低標(biāo)準(zhǔn)船舶而制定的.由于工作環(huán)境的原因，疲勞破壞是船舶結(jié)構(gòu)破壞的主要破壞形式之一.為了保證入級(jí)油船有足夠的疲勞強(qiáng)度，CSR對(duì)油船的疲勞評(píng)估在附錄C中作了詳細(xì)的的規(guī)定［1］.

大多數(shù)情況下船舶結(jié)構(gòu)節(jié)點(diǎn)的幾何形狀和受載情況比較復(fù)雜，常常不能在已有的S－N 曲線中找到可直接比擬的類型，此時(shí)需要采用熱點(diǎn)應(yīng)力法來(lái)進(jìn)行這些部位的疲勞評(píng)估.熱點(diǎn)處的應(yīng)力一般采用插值方法求取，由于熱點(diǎn)區(qū)域一般存在應(yīng)力集中現(xiàn)象，應(yīng)力梯度很高，而且梯度的變化也比較大，所以熱點(diǎn)應(yīng)力的大小與插值方式的選擇有直接關(guān)聯(lián)，不同的插值方式求得的熱點(diǎn)應(yīng)力值會(huì)不一樣，而且可能會(huì)有很大的差別［2］.共同規(guī)范中，規(guī)定在熱點(diǎn)處精細(xì)網(wǎng)格區(qū)內(nèi)使用尺寸為t×t（t為接頭處內(nèi)底板凈厚度）的4節(jié)點(diǎn)板殼單元進(jìn)行離散，熱點(diǎn)應(yīng)力值使用結(jié)構(gòu)交接點(diǎn)第一和第二個(gè)單元的應(yīng)力值內(nèi)插得到，這與以往船級(jí)社的通用做法有一定差別.

由于實(shí)施的時(shí)間還不長(zhǎng)，共同規(guī)范下船體結(jié)構(gòu)的熱點(diǎn)應(yīng)力法疲勞校核結(jié)果和以往各自船級(jí)社規(guī)范校核結(jié)果之間的差異，以及造成差異的影響因素等還沒有得到充分的研究和認(rèn)識(shí)，不利于新規(guī)范的應(yīng)用和實(shí)施.基于此，本文按照共同規(guī)范的規(guī)定對(duì)某一已投入營(yíng)運(yùn)的油船進(jìn)行了的疲勞分析，分析表明，該油船遠(yuǎn)不能滿足CSR的疲勞強(qiáng)度要求，這說(shuō)明CSR對(duì)于結(jié)構(gòu)強(qiáng)度的要求更高，基于CSR設(shè)計(jì)的油船更為安全.文中進(jìn)一步采用不同的插值方式，包括其他船級(jí)社曾經(jīng)采用的插值方式推求疲勞載荷下的熱點(diǎn)處應(yīng)力范圍，研究了不同插值方式對(duì)熱點(diǎn)應(yīng)力范圍計(jì)算值以及進(jìn)一步對(duì)疲勞損傷和疲勞壽命等估算值的影響，并探討了其影響程度，研究表明，正是CSR所采用的插值方式導(dǎo)致該油船的疲勞損傷和疲勞壽命的評(píng)估結(jié)果差強(qiáng)人意，CSR采用的求解熱點(diǎn)應(yīng)力范圍的插值方式是偏于保守的.

1 熱點(diǎn)應(yīng)力法

由于實(shí)際結(jié)構(gòu)往往比較復(fù)雜，名義應(yīng)力的求得以及對(duì)應(yīng)S－N 曲線的獲取不一定方便，所以熱點(diǎn)應(yīng)力法從20世紀(jì)60年代后期開始應(yīng)用于壓力容器和焊接管接頭的疲勞設(shè)計(jì)以來(lái)，應(yīng)用越來(lái)越廣泛，已經(jīng)擴(kuò)展到非管接頭的結(jié)構(gòu)［3］.所謂熱點(diǎn)應(yīng)力是指最大結(jié)構(gòu)應(yīng)力或“結(jié)構(gòu)中危險(xiǎn)截面上危險(xiǎn)點(diǎn)應(yīng)力”，一般多采用有限元計(jì)算或應(yīng)變片實(shí)際測(cè)量法求解.采用有限元法時(shí)，有限元網(wǎng)格的劃分只模擬結(jié)構(gòu)整體尺寸，不反映局部細(xì)微尺寸變化，計(jì)算求得結(jié)構(gòu)中的工作應(yīng)力，不包括焊縫形狀、裂紋、缺口等引起強(qiáng)烈局部應(yīng)力集中，只依賴于構(gòu)件接頭處的宏觀尺寸和載荷參量.實(shí)際應(yīng)用中，由于采用不同的熱點(diǎn)應(yīng)力插值方法和不同的有限元網(wǎng)格，不同研究者得到的疲勞評(píng)估結(jié)果相差較大.從以往研究者的工作和作者已做過的對(duì)比來(lái)看［4］，共同規(guī)范規(guī)定的網(wǎng)格精度是合適的，故僅在此討論插值方式的影響.

2 熱點(diǎn)應(yīng)力求解的插值方式

熱點(diǎn)應(yīng)力法求解時(shí)，比較普遍的是將熱點(diǎn)應(yīng)力定義為焊趾處的應(yīng)力值，此應(yīng)力值又由距焊趾一定距離的2個(gè)點(diǎn)的應(yīng)力值線性外插得到.由于熱點(diǎn)附近應(yīng)力梯度往往很大，求解出的熱點(diǎn)處應(yīng)力范圍大小受插值方式的影響非常明顯.常用比較典型的插值方式有：（1）針對(duì)焊接結(jié)構(gòu)中焊趾位于板的表面的a型熱點(diǎn)，國(guó)際焊接學(xué)會(huì)推薦采用距離焊趾0.4 t和1.0t處應(yīng)力值線性外推計(jì)算焊趾處的熱點(diǎn)應(yīng)力［5－6］，如圖1a）所示.（2）采用距離焊趾0.5 t和1.5 t處應(yīng)力值線性外推計(jì)算焊趾處的應(yīng)力作為熱點(diǎn)應(yīng)力也很常見，見圖1b），這也是很多船級(jí)社采用［7－9］的.（3）在管接頭進(jìn)行疲勞強(qiáng)度分析時(shí)，如C點(diǎn)為焊趾位置，A，B為外插點(diǎn)，如圖1c），則一般A與C之間的距離取，0.4～0.5tA與B之間的距離取為t，其中t為弦管壁厚，焊縫趾端處的熱點(diǎn)應(yīng)力也是通過外插法得到［10］.

圖1 插值示意圖

公同規(guī)范中規(guī)定，對(duì)焊接結(jié)構(gòu)節(jié)點(diǎn)，熱點(diǎn)應(yīng)力范圍指距焊趾位置 處、與焊接方向垂直的表面應(yīng)力（見圖2），而此值使用結(jié)構(gòu)交接點(diǎn)的第一和第二個(gè)單元各自的應(yīng)力通過線性內(nèi)插法得到，其中單元應(yīng)力取自高斯積分點(diǎn)［11］.可以看出，在CSR中計(jì)算熱點(diǎn)應(yīng)力的位置的表述以及插值方式都有所不同.前述3種方法都是外推（插），而共同規(guī)范規(guī)定為內(nèi)插，根據(jù)熱點(diǎn)應(yīng)力的分布特征可知外插方式得到的計(jì)算值會(huì)小于內(nèi)插方式，所以可以定性判斷共同規(guī)范采用的求解熱點(diǎn)應(yīng)力范圍的插值方式較之以前的常用方法嚴(yán)格，這有助于保證油船建造的高標(biāo)準(zhǔn).

圖2 CSR熱點(diǎn)應(yīng)力示意

3 有限元模型和不同插值方式

針對(duì)某油船采用板梁?jiǎn)卧⒘巳摱文Ｐ停ㄒ妶D3）進(jìn)行疲勞評(píng)估，熱點(diǎn)應(yīng)力范圍采用局部細(xì)化模型進(jìn)行計(jì)算（見圖4）.油船有限元計(jì)算模型的建立，包括腐蝕余量折減，工況選取以及約束和載荷的施加等等都按照共同規(guī)范的規(guī)定進(jìn)行，熱點(diǎn)的定義也按照共同規(guī)范采用，分析中惟一不同的只是熱點(diǎn)應(yīng)力范圍計(jì)算方法.考慮到上文（3）針對(duì)的是管接頭，不是板連接件，所以研究中采用了表1所列的4種熱點(diǎn)應(yīng)力范圍的求解方式進(jìn)行對(duì)比研究.

圖3 油船三艙段模型

圖4 局部細(xì)化模型（t×t）

表1 熱點(diǎn)應(yīng)力范圍的求解方式

4 插值方式的影響分析

采用表1所列的不同熱點(diǎn)應(yīng)力范圍計(jì)算方式進(jìn)行了疲勞分析，評(píng)估結(jié)果列于表2.從計(jì)算結(jié)果可以看出：

1）不同方法得到的評(píng)估結(jié)果相差很大，基于共同規(guī)范插值方式（方法1）計(jì)算得出的設(shè)計(jì)年限內(nèi)疲勞損傷（2.416）是船級(jí)社以往常用方式（方法2）得出的設(shè)計(jì)年限內(nèi)疲勞損傷計(jì)算值（1.067）的2.264倍，而方法2預(yù)期疲勞壽命估算值（23.44 a）是方法1預(yù)期疲勞壽命估算值（10.35a）的2.267倍.

2）基于共同規(guī)范插值方式（方法1）的熱點(diǎn)應(yīng)力范圍計(jì)算結(jié)果最大，相應(yīng)的設(shè)計(jì)年限內(nèi)（25a）疲勞損傷最大（2.416）而預(yù)期疲勞壽命最短（10.35a），遠(yuǎn)遠(yuǎn)低于共同規(guī)范關(guān)于油船疲勞強(qiáng)度要求，其計(jì)算結(jié)果和以直接計(jì)算結(jié)果作為熱點(diǎn)應(yīng)力值方式（方法4）的計(jì)算結(jié)果比較接近.

3）按照船級(jí)社以往常用方式（方法2）計(jì)算熱點(diǎn)應(yīng)力范圍，得到的設(shè)計(jì)年限（25a）內(nèi)疲勞損傷最?。?.067）和預(yù)期疲勞壽命最長(zhǎng)（23.44a），接近共同規(guī)范設(shè)計(jì)壽命25a的要求，同時(shí)可以看出該油船確實(shí)滿足當(dāng)時(shí)的入級(jí)標(biāo)準(zhǔn).

4）按照焊接學(xué)會(huì)推薦的熱點(diǎn)插值計(jì)算方法（方法3）得到的評(píng)估結(jié)果與方法2比較接近.

上面4種計(jì)算熱點(diǎn)應(yīng)力范圍的方法中，除方法4是直接采用計(jì)算值外，方法2和方法3采用距交接點(diǎn)遠(yuǎn)于熱點(diǎn)的兩點(diǎn)應(yīng)力外推（插）的方式，方法1采用與交接點(diǎn)距離分別大于和小于熱點(diǎn)的兩點(diǎn)處應(yīng)力范圍值內(nèi)插的方式，由于與交接點(diǎn)越近應(yīng)力和應(yīng)力梯度越大，所以方法2和3確定的熱點(diǎn)應(yīng)力范圍小于直接計(jì)算值；反之方法1確定的計(jì)算值則大于直接計(jì)算值，由于計(jì)算得到的熱點(diǎn)應(yīng)力范圍不同，疲勞評(píng)估結(jié)果也出現(xiàn)了差異.上述結(jié)果說(shuō)明，不同的插值方式得到的疲勞評(píng)估結(jié)果差異很大，插值方式的選取直接影響了疲勞損傷和疲勞壽命的計(jì)算值.

表2 不同插值方式下疲勞評(píng)估

5 結(jié) 論

研究表明，該油船盡管滿足入級(jí)時(shí)相應(yīng)的疲勞強(qiáng)度要求，卻不能滿足共同規(guī)范的疲勞強(qiáng)度規(guī)定，造成這種矛盾的直接原因是，依循不同規(guī)范計(jì)算熱點(diǎn)應(yīng)力范圍時(shí)采用的插值方式不同.共同規(guī)范規(guī)定的熱點(diǎn)應(yīng)力范圍插值方式有助于提高油船的建造標(biāo)準(zhǔn)，保障海上安全，但校核結(jié)果偏于保守.

［1］IACS.Common structural rules for oil tankers［R］.International Association of Classification Societies， 2006.

［2］Fricke W，Kahl A.Comparison of different structural stress approaches for fatigue assessment of welded ship structures［J］.Marine Structures.2005，18（7－8）：473－488.

［3］Zhi G X，Yamada K.A method of determining geometric stress for fatigue strength evaluation of steel welded joints［J］.International Journal of Fatigue，2004，26：1277－1293.

［4］劉 俊，薛鴻祥，陳 芳.某油輪基于共同規(guī)范疲勞評(píng)估及網(wǎng)格精度［J］.上海交通大學(xué)學(xué)報(bào)，2008，42（6）：910－913.

［5］賈法勇，霍立興，張玉鳳.熱點(diǎn)應(yīng)力有限元分析的主要影響因素［J］.焊接學(xué)報(bào)，2003，24（3）：27－30.

［6］Poutiainen I，Tanskanen P，Marquis G.Finite ele－ment methods for structural hot spot stress.International Journal of Fatigue［J］.2004（26）：1147－1157.

［7］ABS.Guide for building and classing membrane tank lng vessels［M］.New York：American Bureau of Shipping，2002.

［8］DNV.Fatigue assessment of ship structures［S］.Det Norske Veritas Classification Notes No30.2003.

［9］BV.Part B－h(huán)ull and stability［S］.Bureau Veritas Rules for the Classification of Steel Ships.2003.

［10］劉 剛，鄭云龍，趙德有.BINGO9000半潛式鉆井平臺(tái)疲勞強(qiáng)度分析［J］.船舶力學(xué)，2002，6（2）：54－63.

［11］羅廣恩，王自力，張延昌.基于JTP規(guī)范中熱點(diǎn)應(yīng)力的油船船體疲勞強(qiáng)度分析［J］.江蘇科技大學(xué)學(xué)報(bào)：自然科學(xué)版，2007，21（1）：17－20.