崔虎威 楊 平

（武漢理工大學(xué)交通學(xué)院 武漢 430063）

0 引 言

在船舶設(shè)計(jì)與強(qiáng)度評(píng)估中，如果不能明確結(jié)構(gòu)的極限承載能力，依舊單純采用基于線彈性理論的許用工作應(yīng)力方法，由于不能從根本上決定結(jié)構(gòu)的真實(shí)安全極限已經(jīng)不再適合當(dāng)前的設(shè)計(jì)理念和發(fā)展要求.將基于極限強(qiáng)度的理論引入船舶結(jié)構(gòu)的設(shè)計(jì)與強(qiáng)度評(píng)估比傳統(tǒng)方法更加合理.IACS《雙殼油船共同規(guī)范》也已將基于“極限強(qiáng)度”概念的非線性屈曲理論引入到技術(shù)標(biāo)準(zhǔn)之中，并于2006年4月1日起開始施行.

在設(shè)計(jì)工作中，準(zhǔn)確地計(jì)算和預(yù)測(cè)結(jié)構(gòu)的極限強(qiáng)度對(duì)于確定結(jié)構(gòu)的安全極限十分重要.雖然多年來對(duì)結(jié)構(gòu)極限強(qiáng)度已有很多理論研究、數(shù)值計(jì)算、模型試驗(yàn)方面的研究工作和成果，但在實(shí)際工程運(yùn)用上仍有相當(dāng)差距.國際船舶與海洋工程領(lǐng)域內(nèi)的學(xué)術(shù)組織ISSC（國際船舶結(jié)構(gòu)大會(huì)）近年來組織了國際上多所研究機(jī)構(gòu)和大學(xué)的知名學(xué)者對(duì)船體極限強(qiáng)度作了較為深入的研究，但研究的結(jié)果表明，盡管采用相同的結(jié)構(gòu)參數(shù)、材料特性、初始缺陷和載荷工況，相互間的結(jié)果仍非常分散，而且與試驗(yàn)也不甚吻合［1］.因此ISSC近年來一直在開展關(guān)于結(jié)構(gòu)極限狀態(tài)的研究工作，其中包括對(duì)板格，加筋板以及船體梁分別采用非線性FEM 進(jìn)行標(biāo)定計(jì)算［2－4］.

影響板格和加筋板極限強(qiáng)度的因素有很多，其中主要有焊接導(dǎo)致的初始缺陷，即焊接變形及焊接殘余應(yīng)力，海水腐蝕等等.初始變形在很大程度上影響板格和加筋板極限強(qiáng)度的大小，因此在各種計(jì)算方法中，初始變形是首要考慮的影響因素.在使用非線性FEM進(jìn)行實(shí)際計(jì)算分析時(shí)，如何有效合理地施加初始變形，選取合適的單元尺寸去獲得計(jì)算時(shí)間與計(jì)算精度的統(tǒng)一，以及邊界條件的確定都將對(duì)極限強(qiáng)度的計(jì)算結(jié)果產(chǎn)生影響.

本文采用非線性有限元分別對(duì)一塊板格和加筋板進(jìn)行具體極限強(qiáng)度計(jì)算，來探究上述問題，并為結(jié)構(gòu)非線性有限元計(jì)算提供參考.

1 結(jié)構(gòu)建模與非線性有限元分析

1.1 板格和加筋板模型

在實(shí)際船舶的甲板板架或船底板架中，板格由縱骨和強(qiáng)橫梁支持，加筋板由縱桁和橫框架支持，見圖1.

圖1 船體板架結(jié)構(gòu)

1.2 板格和加筋板模型的范圍與邊界條件

計(jì)算模型的范圍和邊界條件是相關(guān)聯(lián)的，一般存在著單跨（one－bay）與雙跨（1／2＋1＋1／2－bay）的選擇.

單跨模型不考慮其對(duì)模型轉(zhuǎn)動(dòng)自由度的約束，四邊均采取簡(jiǎn)支并在面內(nèi)保持直線.為了在實(shí)際的設(shè)計(jì)中獲得數(shù)學(xué)上的簡(jiǎn)化，此種模型和邊界經(jīng)常被采用.

雙跨模型由于擴(kuò)展了模型的尺寸范圍，使得模型獲得更加真實(shí)和客觀的轉(zhuǎn)動(dòng)自由度約束.但對(duì)雙跨模型來說，若加筋板的屈曲模態(tài)為奇數(shù)波型，則邊界條件與單跨模型相同；若加筋板的屈曲模態(tài)為偶數(shù)波型，則應(yīng)施加對(duì)稱邊界條件.

實(shí)際建模中為簡(jiǎn)化分析模型，一般不對(duì)支撐構(gòu)件進(jìn)行建模，其影響由施加自由度約束代替.本文對(duì)板格采取單跨模型，對(duì)于加筋板采取帶四根加強(qiáng)筋的雙跨模型為研究對(duì)象.

1.3 模型幾何尺寸與材料屬性及單元?jiǎng)澐?/h3>

板格 板格尺寸L×B×t＝4 300mm×815mm×14mm.

加筋板 板材L×B×t＝4 300mm×815mm×17.8mm；縱骨腹板hw×tw＝436mm×8mm；縱骨面板bf×tf＝172mm×17mm.

構(gòu)件材料選用常見的船用鋼材：E＝205.8 GPa，υ＝0.3，屈服強(qiáng)度бy＝315MPa.

單元尺寸的大小決定了模型網(wǎng)格的疏密程度.一般來講，網(wǎng)格越密則模型越精細(xì)，能使得計(jì)算結(jié)果到更精確.但是在有限元非線性計(jì)算中，網(wǎng)格的疏密也很大程度上決定著計(jì)算時(shí)間，尤其是在帶多根加強(qiáng)筋的加筋板極限強(qiáng)度計(jì)算中，更為明顯.為了在計(jì)算時(shí)間與計(jì)算結(jié)果精確性之中達(dá)到平衡，在確定單元尺寸時(shí)有一個(gè)原則，即在單元的密度增加對(duì)計(jì)算結(jié)果影響微小，可認(rèn)為上一步的單元尺寸已經(jīng)合適.在本文中，經(jīng)過多次計(jì)算比較和參考其他研究者的觀點(diǎn)，采取以下單元尺寸方案：沿板格寬度方向取10個(gè)單元，沿腹板高度方向取6個(gè)單元，沿面板寬度方向取2個(gè)單元即可，沿板和加強(qiáng)筋縱向，根據(jù)單元形狀盡量為正方形要求決定縱向單元尺寸.

1.4 初始變形的施加

現(xiàn)有的研究表明，光板和加筋板結(jié)構(gòu)極限強(qiáng)度對(duì)初始撓度很敏感，所以在分析其極限強(qiáng)度的時(shí)候須將初始變形施加到結(jié)構(gòu)中去.但是由于初始變形的形式相當(dāng)復(fù)雜，如果手動(dòng)進(jìn)行加載初始變形是件非常困難的事情.解決的方法之一是編制擴(kuò)展程序，將初始變形通過改變單元節(jié)點(diǎn)坐標(biāo)來實(shí)現(xiàn)［5］.目前一般的處理方法是先計(jì)算結(jié)構(gòu)的屈曲，將低階屈曲模態(tài)以指定的變形幅值施加到計(jì)算模型中.雖然這種方法在計(jì)算結(jié)構(gòu)的屈曲或極限強(qiáng)度時(shí)能大致計(jì)入初始變形的影響，但與實(shí)際結(jié)構(gòu)的初始變形一般情況下并不一致，因此得到的結(jié)構(gòu)的屈曲或極限強(qiáng)度結(jié)果將會(huì)出現(xiàn)偏差.

初始變形能夠明顯影響板和加筋板的極限強(qiáng)度特性，所以采取合適的方法去施加初始變形就顯得很關(guān)鍵.

對(duì)于板格而言，只需考慮板自身的局部屈曲，本文中將板的初始缺陷與不同的加載形式相關(guān)聯(lián)，通過提取第一階模態(tài)，擴(kuò)大到幅值Wopl＝B／200即可完成初始變形的施加.

對(duì)于加筋板，則復(fù)雜得多.本文考慮雙向受壓的加載情況.計(jì)入3種與屈曲相關(guān)聯(lián)的初始缺陷，同時(shí)考慮加載方式對(duì)屈曲模態(tài)以及初始變形的影響.

1）與雙向受壓屈曲模態(tài)相關(guān)聯(lián)的板的初始缺陷，幅值為B／200.

2）與雙向受壓屈曲模態(tài)相關(guān)聯(lián)，梁柱型的筋的初始缺陷，幅值為L(zhǎng)／1 000.

3）與雙向受壓屈曲模態(tài)相關(guān)聯(lián)，筋的側(cè)傾型的初始缺陷，幅值為L(zhǎng)／1000.

在計(jì)算中，首先要進(jìn)行特征值屈曲分析，而后據(jù)所得屈曲模態(tài)分別挑選出與以上3種初始變形形式相似的模態(tài)，使用Ansys中的Upgeom功能，分別挑選出與板和筋相關(guān)的節(jié)點(diǎn)進(jìn)行單獨(dú)的變形施加操作，最后將各部分?jǐn)U大到所需幅值后，疊加在一起就完成了加筋板初始變形的施加.

顯然，如何去挑選模態(tài)和計(jì)算者的判斷有關(guān)系.本文在處理時(shí)遵循了以下原則：板和加強(qiáng)筋的模態(tài)首先要均勻分布，例如板上“峰”、“谷”分布要均勻，筋的模態(tài)也要均勻協(xié)調(diào).在滿足均勻性要求時(shí)，以傾向低階模態(tài)為原則.選擇低階是因?yàn)槟B(tài)較低的屈曲優(yōu)先發(fā)生，進(jìn)而達(dá)到屈曲所需的壓力要小，從而使得所計(jì)算的極限強(qiáng)度值偏于安全.顯然根據(jù)均勻性的要求，加筋板的初始缺陷并不是簡(jiǎn)單地鎖定第一階.

2 非線性有限元計(jì)算

2.1 板格和加筋板的復(fù)合受壓載荷

選取以下4種載荷模式對(duì)板格進(jìn)行計(jì)算，考慮雙向軸壓與側(cè)壓的復(fù)合作用，見表1.

表1 板格上雙向軸壓與側(cè)壓載荷

選取以下4種載荷模式進(jìn)行計(jì)算，重點(diǎn)考慮不同軸向壓力加載下的加筋板極限強(qiáng)度特性，見表2.

表2 加筋板上雙向軸壓與側(cè)壓載荷

其中板格的計(jì)算由于考慮到了側(cè)壓的作用，將導(dǎo)致在非線性計(jì)算中產(chǎn)生多載荷步的計(jì)算問題.本文在處理時(shí)先施加側(cè)壓，以得到在前一步進(jìn)行特征值屈曲分析所獲初始缺陷基礎(chǔ)上的再次變形結(jié)果，而后施加軸向壓力.

2.2 非線性有限元計(jì)算

板格及加筋板的有限元模型見圖2～3，板格和加筋板的計(jì)算極限狀態(tài)時(shí)的von Mises應(yīng)力及應(yīng)力與應(yīng)變關(guān)系見圖4～13.

圖2 板格有限元模型

圖3 加筋板有限元模型

圖4 板格極限狀態(tài)時(shí)的應(yīng)力分布（σx∶σy＝1.0∶0.0）

圖5 板格極限狀態(tài)時(shí)的應(yīng)力分布（σx∶σy＝0.8∶0.2）

圖6 板格極限狀態(tài)時(shí)的應(yīng)力分布（σx∶σy＝0.7∶0.3）

圖7 板格極限狀態(tài)時(shí)的應(yīng)力分布（σx∶σy＝0.5∶0.5）

圖8 加筋板極限狀態(tài)時(shí)的應(yīng)力分布（σx∶σy＝1.0∶0.0）

圖9 加筋板極限狀態(tài)時(shí)的應(yīng)力分布（σx∶σy＝0.79∶0.21）

圖10 加筋板極限狀態(tài)時(shí)的應(yīng)力分布（σx∶σy＝0.4∶0.6）

圖11 加筋板極限狀態(tài)時(shí)的應(yīng)力分布（σx∶σy＝0.0∶1.0）

圖12 板格的應(yīng)力與應(yīng)變關(guān)系曲線

圖13 加筋板的應(yīng)力與應(yīng)變關(guān)系曲線

3 計(jì)算結(jié)果分析

現(xiàn)將本文所得計(jì)算值與文獻(xiàn)［6－7］進(jìn)行對(duì)比，結(jié)果見表3～4.其中計(jì)算值與參考值均為極限狀態(tài)時(shí)刻工作應(yīng)力與屈服強(qiáng)度的比值.

表3 板格計(jì)算值與文獻(xiàn)值的比較

表4 加筋板計(jì)算值與文獻(xiàn)值的比較

由表3可見，就板格而言，計(jì)算值與參考值之間誤差很微小.由表4可見，對(duì)于加筋板來說，兩者最大誤差在5%左右，以加筋板的復(fù)雜程度來看，計(jì)算結(jié)果是可以接受的.

板格和加筋板的誤差有此區(qū)別，是由于加筋板的初始缺陷的施加比板格復(fù)雜得多，還與計(jì)算者自己的判斷有很大關(guān)系的原因.

4 結(jié)束語

本文分別對(duì)2塊分別來源于實(shí)船結(jié)構(gòu)的板格和加筋板進(jìn)行了非線性有限元的極限強(qiáng)度計(jì)算，重點(diǎn)就兩者計(jì)算模型的選取，邊界條件的施加，以及初始變形的實(shí)現(xiàn)等方面進(jìn)行了論述.

從計(jì)算所得的結(jié)果和已有的研究結(jié)果相比較可知，兩者符合的比較好，證明了本文所采取的計(jì)算方法和處理方法是正確的，同時(shí)也佐證了Ansys作為一款功能強(qiáng)大的通用軟件，其在進(jìn)行極限強(qiáng)度非線性分析時(shí)的有效性.

本文極限強(qiáng)度的計(jì)算及處理方法可為結(jié)構(gòu)非線性計(jì)算分析提供一定參考.

［1］Yao T，Astrup O C，Caridis P，et al.Ultimate hull girder strength：proc.of 14th international ship and offshore structures congress［C］.Nagasaki，Japan，2000，2：321－391.

［2］Simonsen B C，Estefen S F.Ultimate strength：proc.of 15th international ship and offshore structures congress［C］.San Diego，USA，2003，1：265－328.

［3］Yao T，Brunner E.Ultimate strength：proc.of 16th international ship and offshore structures congress［C］.Southampton，UK，2006（1）：353－438.

［4］Paik J K，Brunner E.Ultimate strength：proc.of 17th international ship and offshore structures con－gress［C］.Seoul，Korea，2009，1：375－474.

［5］Zhang S，Khan I.Buckling and ultimate capability of the plates and stiffened panels in anxial compression［J］.Marine Structures，2009（22）：791－808.

［6］Paik J K，Seo J K.Nonlinear finite element method models for ultimate strenth analysis of steel stiffenedplate structures under combined biaxial compression and lateral pressure actions—Part I：Plate elements［J］.Thin－walled Structures，2009（47）：1 008－1 017.

［7］Paik J K，Seo J K.Nonlinear finite element method models for ultimate strenth analysis of steel stiffenedplate structures under combined biaxial compression and lateral pressure actions—Part II：Stiffened panels［J］.Thin－walled Structures，2009（47）：998－1 017.