許 俊 吳劍國 葉 帆 王凡超

（1.浙江工業(yè)大學(xué) 建筑工程學(xué)院，杭州 310014；2.中國船舶及海洋工程設(shè)計研究院 上海 200011）

引 言

船舶是以加筋板格為主要承力構(gòu)件的板梁組合結(jié)構(gòu)，對其開展循環(huán)載荷下的極限強度分析，研究加筋板格在循環(huán)載荷作用下的應(yīng)力應(yīng)變關(guān)系，對于進一步理解和研究船體梁在波浪載荷作用下的破壞機理十分重要。

20世紀(jì)80年代中期，日本學(xué)者Y.Fukumoto和H. Kusama［1-2］對板材和薄壁箱型梁等結(jié)構(gòu)在循環(huán)加載下的性能進行一系列實驗研究，獲得結(jié)構(gòu)平均應(yīng)力-平均應(yīng)變的滯回曲線。90年代初期，T. Yao［3］研究了構(gòu)件在加載-卸載-反向再加載這種循環(huán)載荷作用下結(jié)構(gòu)發(fā)生的屈曲、屈服崩潰行為。近年來，Shengming Zhang［4-5］運用非線性有限元方法對加筋板格在殘余應(yīng)力下的極限強度進行研究。

武漢理工大學(xué)單成巍、劉格暢等人［6-7］也對這一問題進行相關(guān)研究，他們對傳統(tǒng)的逐步崩潰法進行改進，并采用非線性有限元方法進行循環(huán)載荷作用下的加筋板逐步崩潰行為的研究，但他們都是對單一因素進行簡單分析。本文將多種因素進行綜合考慮，采用非線性有限元方法，分析殘余應(yīng)力、隨動強化和側(cè)向荷載混合作用對極限強度的影響。

1 有限元模型

加筋板格單元指的是由一根加強筋和與其相連的帶板所組成的結(jié)構(gòu)，它是船體結(jié)構(gòu)的重要承載構(gòu)件。當(dāng)加筋板格單元受到軸向壓力作用時，會發(fā)生以下5種失效模式：

（1）加強筋之間板的屈曲；

（2）加筋板格單元的面板受壓屈曲/屈服失效；

（3）加筋板格單元的帶板受壓屈曲/屈服失效；

（4）加筋板格單元的帶板受壓屈曲/屈服和面板受拉屈服組合破壞導(dǎo)致的失效；

（5）加強筋的側(cè)傾或筋的腹板局部屈曲導(dǎo)致的失效。

本文采用的加筋板格計算模型為單筋單跨模型，如圖1所示，b為帶板寬度，a為跨度（即強橫梁間距）。循環(huán)加載方式為軸向加載。

圖1 加筋板格與加筋板格單元示意圖

1.1 尺寸與邊界

加筋板格的尺寸參數(shù)見表1，應(yīng)力應(yīng)變關(guān)系為理想彈塑性。

表1 加筋板格模型尺寸

為在實際設(shè)計中得以簡化，單跨模型不考慮其對模型轉(zhuǎn)動自由度的約束，計算模型取1/2+1/2個橫梁間距，C1-D1為橫梁處，橫梁未建模，輔以邊界條件約束。采用參考文獻［5］中提出的約束方式，橫梁與縱筋相交處的筋腹板的Y向位移相同，以限制筋在橫梁處的側(cè)向變形。如下頁圖2所示：與橫梁平行的A1-A2邊和B1-B2邊施加對稱邊界，B1-B2邊上的所有節(jié)點的X方向位移被約束以限制剛體運動。A1-B1邊上所有節(jié)點的Y方向位移被約束。為利于收斂，循環(huán)加載方式采用位移控制，在A1-A2邊施加X方向位移， A1-A2的所有節(jié)點在X方向均有相同的位移。本文建模采用ABAQUS非線性有限元軟件。

圖2 單跨模型

圖3 模型初始?xì)堄鄳?yīng)力分布

1.2 幾何初始缺陷

加筋板均由焊接工藝制造而成，因此，板格中必然存在初始缺陷。初始撓度的有限元模擬一般分為局部缺陷和整體缺陷。如圖2所示，模型在G1點建立基準(zhǔn)坐標(biāo)系，a為橫梁的跨度， 為腹板高度，b為筋的間距，B為板橫向邊界的總長；板格屈曲半波數(shù)m滿足：（1）局部缺陷：

（2）整體缺陷：

式中：Wg、Wy分別為帶板和筋腹板的初始撓曲方程；a/1 000為帶板和筋腹板的最大初始撓曲值。

1.3 殘余應(yīng)力

鋼材焊接是局部受熱。在焊接階段，受熱的材料將會膨脹，但由于相鄰低溫度材料的限制，會產(chǎn)生壓應(yīng)力。在冷卻階段，受熱后的材料將會迅速恢復(fù)其原有形狀，從而在熱影響區(qū)產(chǎn)生拉應(yīng)力。在加筋板結(jié)構(gòu)中，焊縫沿著加強筋腹板與板相連區(qū)域分布（殘余拉應(yīng)力區(qū)） ；板的中部則分布著與之相平衡的殘余壓應(yīng)力。這些應(yīng)力統(tǒng)稱為殘余應(yīng)力。本文有限元模型殘余應(yīng)力分布如圖3所示。加筋板中殘余應(yīng)力的分布一般如下頁圖4（a）所示，為簡化，采用如圖4（b）所示的殘余應(yīng)力分布。其中，拉應(yīng)力σrtx、σrty分布在板的四周與縱骨、橫梁相交處，壓應(yīng)力σrcx、σrcy分布在板的中間。為達到平衡，需滿足如下公式：

ISSC（2009）［8］認(rèn)為帶板上Y方向殘余應(yīng)力及筋上殘余應(yīng)力對極限強度影響較小，僅考慮帶板上X方向殘余應(yīng)力分布，本文亦僅考慮帶板上X方向殘余應(yīng)力分布。

1.4 包辛格效應(yīng)

金屬在一個方向發(fā)生塑性變形后再反向變形時，其屈服強度下降的現(xiàn)象稱為包辛格效應(yīng)。加筋板中亦存在著包辛格效應(yīng)，采用ABAQUS軟件進行線性隨動強化模擬，下頁圖5為包辛格效應(yīng)的應(yīng)力應(yīng)變示意圖，強化模量參數(shù)KBC（KBC為BC線的斜率）為1 000 N/mm2。

圖4 殘余應(yīng)力分布

圖5 包辛格效應(yīng)示意圖

2 殘余應(yīng)力的影響分析

2.1 單向循環(huán)

單向循環(huán)方式為：0—1.4倍屈服應(yīng)變—2.1倍屈服應(yīng)變—3.2倍屈服應(yīng)變，考慮10%殘余壓應(yīng)力（即0.1倍屈服應(yīng)力）和15%殘余壓應(yīng)力兩種情況。圖6為模型在無殘余應(yīng)力和10%殘余壓應(yīng)力下，單向受壓循環(huán)的應(yīng)力應(yīng)變對比曲線。

2.2 雙向循環(huán)

應(yīng)變幅為1.4倍屈服應(yīng)變，循環(huán)4次；考慮10%殘余壓應(yīng)力和15%殘余壓應(yīng)力兩種情況。

圖6 模型不同殘余應(yīng)力的應(yīng)力應(yīng)變曲線

由圖6可知，殘余應(yīng)力會使加筋板達到極限強度時處于較大的應(yīng)變值，在后屈曲階段，殘余應(yīng)力存在使加筋板的結(jié)構(gòu)能力更大，無殘余應(yīng)力加筋板的極限強度下降更快。

模型在1.4倍屈服應(yīng)變幅時，不同殘余壓應(yīng)力下，隨循環(huán)次數(shù)變化下的曲線對比如下頁圖7所示。由該圖可知，有殘余應(yīng)力時，第二次循環(huán)時加筋板的極限強度比第一次大，隨著循環(huán)次數(shù)增加，加筋板的極限強度開始下降，下降幅度先大后小。

圖7 模型在1.4倍屈服應(yīng)變幅時，不同殘余壓應(yīng)力下隨循環(huán)次數(shù)變化下的對比曲線

模型在15%殘余壓應(yīng)力和3.2倍屈服應(yīng)變幅時，第n次循環(huán)達極限強度的應(yīng)力云圖見圖8。由該圖可知，第1次循環(huán)達到極限強度時，板和筋的絕大部分都達到了屈服應(yīng)力；隨著循環(huán)次數(shù)增加，加筋板達到極限狀態(tài)時，高應(yīng)力區(qū)越來越向橫梁之間的區(qū)域集中；第4次達到循環(huán)極限強度時，只有靠近橫梁之間區(qū)域的一小部分板和筋的應(yīng)力達到屈服強度。經(jīng)過多次循環(huán)后，加筋板極限強度下降變得緩慢。

圖8 模型在15%殘余壓應(yīng)力和3.2倍屈服應(yīng)變幅時，第n次循環(huán)達極限強度的應(yīng)力云圖

3 殘余應(yīng)力和側(cè)向荷載的混合作用

3.1 單向循環(huán)

下頁圖9（a）為在側(cè)壓0.16 MPa時，無殘余應(yīng)力和10%殘余壓應(yīng)力下，模型單向受壓循環(huán)時的應(yīng)力應(yīng)變對比曲線。

3.2 雙向循環(huán)

雙向循環(huán)方式為先壓后拉，循環(huán)4次。殘余應(yīng)力考慮3種工況：無殘余應(yīng)力、10%殘余壓應(yīng)力、15%殘余壓應(yīng)力；側(cè)壓考慮兩種工況：0 MPa、0.16 MPa；應(yīng)變幅考慮3種情況：1.4倍屈服應(yīng)變、2.1倍屈服應(yīng)變、3.2倍屈服應(yīng)變。

圖9 模型在0.16 MPa側(cè)壓時，不同殘余應(yīng)力的應(yīng)力應(yīng)變曲線

3.2.1 有側(cè)壓時，考慮不同殘余應(yīng)力

圖9（b）為0.16 MPa側(cè)壓、不同殘余壓應(yīng)力且等應(yīng)變幅時，模型在循環(huán)四次過程中應(yīng)力應(yīng)變對比曲線，圖10為模型0.16 MPa側(cè)壓、不同殘余壓應(yīng)力下，每次循環(huán)時的極限強度對比曲線。

圖10 模型0.16 MPa側(cè)壓時，不同殘余壓應(yīng)力隨著循環(huán)次數(shù)的應(yīng)力對比曲線

3.2.2 有側(cè)壓和相同殘余壓應(yīng)力時，不同的等幅應(yīng)變

下頁圖11為0.16 MPa側(cè)壓、15%殘余壓應(yīng)力時，模型在不同的等幅應(yīng)變下，4次循環(huán)過程中的應(yīng)力應(yīng)變曲線。

由圖9（a）可知，在一定側(cè)壓下，單向循環(huán)應(yīng)變幅較大時，不同殘余壓應(yīng)力的極限強度趨向同一值。由圖9（b）和圖10可知，在一定側(cè)壓下，殘余壓應(yīng)力較大時，不同的殘余壓應(yīng)力對極限強度的影響不明顯。由圖11可知，循環(huán)應(yīng)變幅很大時，隨著循環(huán)次數(shù)增加，不同殘余壓應(yīng)力下的極限強度趨向同一值，側(cè)壓使極限強度下降趨勢變大。等幅應(yīng)變幅越大，循環(huán)時的極限強度下降越多。殘余壓應(yīng)力、側(cè)壓和循環(huán)應(yīng)變幅均較大時，為最不利組合。

圖11 在0.16MPa側(cè)壓和15%殘余壓應(yīng)力時，不同等幅應(yīng)變下的對比曲線

4 殘余應(yīng)力和隨動強化的混合作用

僅考慮雙向循環(huán)，循環(huán)方式均為先壓后拉，每次均為等應(yīng)變幅。應(yīng)變幅為2.1倍屈服應(yīng)變，循環(huán)4次。殘余壓應(yīng)力考慮3種工況：不考慮殘余壓應(yīng)力、只考慮15%殘余壓應(yīng)力、同時考慮15%殘余壓應(yīng)力和隨動強化效應(yīng)。

由圖12可知，多次循環(huán)后，隨動強化效應(yīng)對加筋板的極限強度幾乎毫無影響。

圖12 在2.1倍屈服應(yīng)變幅時，不同殘余壓應(yīng)力下的應(yīng)力應(yīng)變曲線

5 結(jié) 論

循環(huán)應(yīng)變幅度、殘余壓應(yīng)力和側(cè)向荷載均能顯著降低加筋板的極限承載力，0.16 MPa側(cè)壓、15%殘余壓應(yīng)力、3.2倍屈服應(yīng)變幅時，循環(huán)4次后的極限強度下降39%（相比無殘余、無側(cè)壓、單次加載的極限強度）。在后屈曲階段，殘余應(yīng)力的存在使加筋板的結(jié)構(gòu)能力更大，無殘余應(yīng)力的加筋板的極限強度下降更快。多次循環(huán)后，隨動強化效應(yīng)對加筋板的極限強度幾乎無影響。應(yīng)力狀態(tài)重分布所造成的高應(yīng)力向中部轉(zhuǎn)移是導(dǎo)致加筋板循環(huán)時極限強度下降的主要原因。

本文后續(xù)還將對大跨及破損模型進行循環(huán)載荷研究，并對滯回曲線的耗能進行分析，多角度揭示加筋板極限強度的影響規(guī)律。

［參考文獻］

［1］ FUKUMTO Y， KUSAMA H. Local instability tests of plate elements under cyclic uniaxial loading. Journal of Structural Engineering Division［J］. ASCE， 1985， 111（ST5）：1051-1067.

［2］ FUKUMTO Y， KUSAMA H. Cyclic behavior of Plates under in-plane loading ［J］. Engineering Structures， 1985（7）：56-63.

［3］ YAO T， NIKOLOV P I. Buckling/plastic collapse of plates under cyclic loading［J］. Journal of the Society of Naval Architects of Japan， 1990， 168 ： 449-462.

［4］ IMTAZK，SHENGMING Z. Effects of welding-induced residual stress on ultimate strength of plates and stiffened panels［J］. Ships and Offshore Structures， 2011， 6（4）：297-309.

［5］ IMTAZ K，SHENGMING Z. Buckling and ultimate capability of plates and stiffened panels in axial compression［J］. Marine Structures， 2009， 22 ： 791-808.

［6］ 劉恪暢.循環(huán)彎曲下船體梁極限強度［D］. 武漢： 武漢理工大學(xué)， 2012.

［7］ 單成巍.循環(huán)荷載作用下船體結(jié)構(gòu)的極限強度非線性的有限元分析［D］. 武漢：武漢理工大學(xué)， 2013.

［8］ BRANNER K，CZUJKO J. Committee III.1-Ultimate Strength ［C］//Proceedings of the 17thInternational Ship and Offshore Structures Congress （ISSC）. Seoul，korea，2009： 375-474.