周 雷 朱 凌 李應(yīng)剛* 余同希 郭開嶺 蔡 偉

(武漢理工大學(xué)高性能船艦技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室1) 武漢 430063) (武漢理工大學(xué)交通學(xué)院2) 武漢 430063) (高新船舶與深海開發(fā)裝備協(xié)同創(chuàng)新中心3) 武漢 430063) (香港科技大學(xué)力學(xué)與航空工程系4) 香港 999077)

0 引 言

全球氣候變暖，北冰洋海冰加速融化，極地蘊(yùn)藏的豐富資源和北極航道戰(zhàn)略價(jià)值都在日益突顯.隨著北極航道的開通，將會(huì)有越來(lái)越多的船舶進(jìn)入北極航道進(jìn)行能源開發(fā)與極地科學(xué)考察，然而北極航道上覆蓋著冰層，隨機(jī)分布大量冰排及浮冰，給船舶航行帶來(lái)了極大的困難和潛在的威脅.因此，極地船舶與冰相互作用及結(jié)構(gòu)抗冰碰性能問(wèn)題愈來(lái)愈引起人們的關(guān)注.

目前，船體結(jié)構(gòu)抗冰體碰撞研究主要集中在試驗(yàn)方法和數(shù)值計(jì)算方法.Gagnon等[1-3]進(jìn)行了冰體與Terry Fox破冰船碰撞試驗(yàn)，還進(jìn)行了冰體與船體結(jié)構(gòu)碰撞模型試驗(yàn)，分析了碰撞過(guò)程中結(jié)構(gòu)的接觸壓力和碰撞力隨時(shí)間變化歷程，并與有限元計(jì)算結(jié)果進(jìn)行了對(duì)比.Ritch等[4]也進(jìn)行了Terry Fox破冰船與小型冰山碰撞試驗(yàn)，測(cè)出了船體結(jié)構(gòu)局部碰撞過(guò)程中的局部壓力.Manuel等[5]進(jìn)行了一系列淡水冰與船舶結(jié)構(gòu)的準(zhǔn)靜態(tài)壓縮實(shí)驗(yàn)，船舶結(jié)構(gòu)產(chǎn)生了較大的塑性變形.Kim等[6-7]進(jìn)行了不同加載速率和不同壓頭的冰體壓縮試驗(yàn).Combescur等[8]進(jìn)行了冰體的高速?zèng)_擊試驗(yàn)，研究了冰體破碎失效和板的變形.張健等[9]通過(guò)模型試驗(yàn)的方法對(duì)船體舷側(cè)板架與冰體碰撞進(jìn)行了研究.Song等[10-11]進(jìn)行了實(shí)驗(yàn)室制取淡水冰體與加筋板結(jié)構(gòu)沖擊試驗(yàn)，冰體和結(jié)構(gòu)一起發(fā)生破壞.數(shù)值計(jì)算中冰體材料的數(shù)值模型是一個(gè)難題，目前主要采用的是各向同性彈塑性模型、彈性斷裂模型、可壓碎泡沫型模型[12-13]，以及一些其它自定義材料數(shù)值模型.Kim等[14-15]采用各向同性彈性斷裂模型(*MAT_ISOTROPIC_ELASTIC_FAILURE)來(lái)對(duì)冰體材料進(jìn)行數(shù)值模擬.Zhu等[16-17]利用該材料進(jìn)行了冰與船體板碰撞數(shù)值仿真，提出了一種可以保守估計(jì)局部冰載荷的計(jì)算方法，在此基礎(chǔ)上還進(jìn)行了船體板在楔形冰重復(fù)碰撞下的動(dòng)態(tài)響應(yīng)研究，研究了板的塑性累積損傷和能量吸收情況.

在上述研究中，與冰撞擊的對(duì)象幾乎都是船體鋼板，相對(duì)于船體鋼板而言，夾層結(jié)構(gòu)具有更好的抗沖擊性能.Guo等[18]對(duì)泡沫鋁夾層板進(jìn)行了重復(fù)沖擊實(shí)驗(yàn)，研究了泡沫鋁夾層板抗重復(fù)沖擊性能，并分析了結(jié)構(gòu)響應(yīng)隨著沖擊次數(shù)的變化規(guī)律.Zhu等[19]利用實(shí)驗(yàn)方法研究了低溫對(duì)泡沫鋁夾層板單次和重復(fù)沖擊結(jié)構(gòu)響應(yīng)的影響規(guī)律.盡管已有很多人對(duì)泡沫鋁夾層板抗沖擊性能進(jìn)行了研究，然而，對(duì)泡沫鋁夾層板抗冰體沖擊性能的研究相對(duì)較少，冰體破碎對(duì)泡沫鋁夾層板動(dòng)態(tài)沖擊響應(yīng)的影響尚不清楚.因此，本文利用ANSYS/LS-DYNA有限元軟件, 研究了冰碰載荷作用下船用泡沫鋁夾層板的動(dòng)態(tài)沖擊響應(yīng).在此基礎(chǔ)上，對(duì)比分析了冰體撞擊和剛體撞擊對(duì)泡沫鋁夾層板動(dòng)態(tài)沖擊響應(yīng)的差異，為極地船舶結(jié)構(gòu)抗沖擊防護(hù)設(shè)計(jì)提供了理論依據(jù).

1 楔形冰撞擊泡沫鋁夾層板數(shù)值仿真

1.1 仿真模型

泡沫鋁夾層板由上、下面板及泡沫鋁芯層構(gòu)成，結(jié)構(gòu)簡(jiǎn)圖見圖1.本文采用ANSYS/LS-DYNA對(duì)泡沫鋁夾層板受楔形冰撞擊的動(dòng)態(tài)沖擊過(guò)程進(jìn)行了模擬，建立了圖2所示的有限元模型.夾層板面板的尺寸為400 mm×400 mm，上、下面板厚度為1 mm，芯層厚度為20 mm，楔形冰對(duì)中碰撞泡沫鋁夾層板，其幾何尺寸見圖3.泡沫鋁芯層和楔形冰體采用實(shí)體單元(solid164)，面板采用殼單元(shell163).面板與芯層之間采用共節(jié)點(diǎn)形式，楔形冰與面板之間定義侵蝕面面接觸(eroding surface to surface).泡沫鋁夾層板四周剛性約束，楔形冰體保留碰撞方向的自由度，其他自由度全部約束，賦予楔形冰體初速度3.3 m/s.為了提高計(jì)算精度，將楔形冰體前段碰撞區(qū)的網(wǎng)格進(jìn)行局部加密.

圖2 楔形冰-泡沫鋁夾層板有限元模型

圖3 碰撞模型幾何尺寸

圖1 泡沫鋁夾層板結(jié)構(gòu)簡(jiǎn)圖

泡沫鋁采用可壓碎泡沫模型(crushable foam)，彈性部分定義彈性模量和泊松比，塑性部分參考Zhu等[20]的實(shí)驗(yàn)數(shù)據(jù)輸入真實(shí)應(yīng)力應(yīng)變曲線轉(zhuǎn)化得到的塑性應(yīng)力應(yīng)變曲線，見圖4a)，鋁面板采用分段線性塑性模型(piecewise linear plasticity)，彈性部分定義彈性模量和泊松比，塑性部分參考Xi等[21]的實(shí)驗(yàn)數(shù)據(jù)輸入真實(shí)應(yīng)力應(yīng)變曲線轉(zhuǎn)化得到的塑性應(yīng)力應(yīng)變曲線，見圖4b).

圖4 塑性應(yīng)力與塑性應(yīng)變間的關(guān)系曲線、

1.2 冰材料模型及驗(yàn)證

就目前的研究來(lái)說(shuō)，由于冰的復(fù)雜性，冰的材料特性的模擬仍處于不斷探索當(dāng)中.本文采用ANSYS/LS-DYNA材料庫(kù)中彈性斷裂模型*MAT_013(MAT_ISOTROPIC_ELASTIC_FAILTURE)來(lái)模擬冰的材料模型，彈性斷裂材料模型采用von-Mises屈服準(zhǔn)則作為冰的破壞準(zhǔn)則.當(dāng)有效塑性應(yīng)變大于失效應(yīng)變或壓力小于截?cái)鄩毫r(shí)冰材料就會(huì)發(fā)生破碎失效.結(jié)合文獻(xiàn)[14]的實(shí)驗(yàn)數(shù)據(jù)，給出該模型的基本材料參數(shù)，見表1.

表1 冰的材料參數(shù)(*MAT_013)

為了驗(yàn)證冰的本構(gòu)模型及其參數(shù)的合理性, 開展了楔形冰撞擊船體板實(shí)驗(yàn)，見圖5.實(shí)驗(yàn)采用自行設(shè)計(jì)的水平?jīng)_擊試驗(yàn)裝置[22].實(shí)驗(yàn)中，楔形冰以3.3 m/s的速度正撞船體板中心，低碳鋼板的尺寸為1 200 mm×400 mm×5.85 mm.采用LS-DYNA中MAT-24定義船體板材料模型，材料的密度7 850 kg/m3、彈性模量206 GPa，輸入真實(shí)應(yīng)力應(yīng)變曲線見圖6.圖7將實(shí)驗(yàn)和仿真的碰撞加速度進(jìn)行了比較，由圖7可知，數(shù)值仿真結(jié)果與實(shí)驗(yàn)結(jié)果吻合得較好，驗(yàn)證了此種冰材料模型和數(shù)值仿真方法的可靠性.

圖5 楔形冰-船體板碰撞實(shí)驗(yàn)

圖6 船體板真實(shí)應(yīng)力-應(yīng)變曲線

圖7 實(shí)驗(yàn)和仿真碰撞加速度對(duì)比

2 冰體碰撞和剛體碰撞泡沫鋁夾層板結(jié)構(gòu)動(dòng)態(tài)響應(yīng)

為了比較冰體和剛體撞擊泡沫鋁夾層板的區(qū)別，按照?qǐng)D3的幾何尺寸建立剛體撞擊泡沫鋁夾層板的數(shù)值模型.保證沖擊能量和沖頭形狀和不變，研究泡沫鋁夾層板在剛性頭和冰頭沖擊下的沖擊動(dòng)態(tài)響應(yīng).

2.1 泡沫鋁夾層板沖擊響應(yīng)過(guò)程分析

為了研究泡沫鋁夾層板變形模態(tài)在冰碰和剛碰下有何不同，截取不同時(shí)刻不同中剖面下泡沫鋁夾層板上、下面板的變形模態(tài)進(jìn)行比較，結(jié)果見圖8.由圖8可知，在楔形冰體的撞擊下，泡沫鋁夾層板的上、下面板有著相似的變形模態(tài)；在楔形剛體的撞擊下，泡沫鋁夾層板的上、下面板變形模態(tài)存在差異，上面板在楔形剛體的撞擊下有著明顯的塑性鉸.兩種撞擊工況下，泡沫鋁夾層板的上面板變形模態(tài)存在差異的原因是：冰頭相對(duì)于夾層板的強(qiáng)度要弱，當(dāng)作用力或者變形超過(guò)了楔形冰單元的承受范圍，楔形冰單元網(wǎng)格就會(huì)失效，從而不會(huì)在接觸區(qū)產(chǎn)生較大的作用力，而剛性頭相對(duì)與夾層板的強(qiáng)度要強(qiáng)，并且剛性沖頭不會(huì)變形，因此當(dāng)剛頭與上面板發(fā)生碰撞時(shí)，沖頭兩尖端的接觸力會(huì)較大，導(dǎo)致上面板與尖端相接觸區(qū)域局部變形嚴(yán)重，從而產(chǎn)生明顯的塑性鉸.兩種撞擊情況下，碰撞力隨時(shí)間的變化曲線見圖9.圖10為泡沫鋁夾層板的上、下面板在冰體撞擊和剛體撞擊下的變形時(shí)程曲線，由圖10可知，在相同沖擊能量和沖頭形狀下，剛體撞擊產(chǎn)生的變形明顯比冰體大.從材料特性上分析，冰體破碎會(huì)使得碰撞力減小，因此相較于剛體撞擊，冰體撞擊工況下泡沫鋁夾層板上、下面板的變形較小.

圖8 兩種撞擊工況下泡沫鋁夾層板變形模態(tài)對(duì)比圖

圖9 碰撞力對(duì)比圖

圖10 上、下面板變形對(duì)比圖

2.2 泡沫鋁夾層板能量吸收分析

為了研究泡沫鋁夾層板能量吸收情況在冰碰和剛碰下有何不同，分別給出了泡沫鋁夾層板各個(gè)部分的能量時(shí)程曲線和撞擊過(guò)程能量轉(zhuǎn)化圖，見圖11～12.由圖11可知，在剛體碰撞下，上面板吸收的能量最多，其次是下面板，泡沫鋁芯層吸收的能量最少；在冰體碰撞下，上面板和泡沫鋁芯層吸收的能量相差無(wú)幾，下面板吸收的能量最少.由圖12可知，在初始碰撞能量相同的情況下，剛體碰撞中泡沫鋁夾層板吸收的能量明顯大于冰體碰撞中泡沫鋁夾層板吸收的能量.在剛體碰撞中，總的碰撞能分成了兩個(gè)部分，絕大部分能量被泡沫鋁夾層板所吸收，成為了變形能，小部分能量轉(zhuǎn)化為剛體的反彈能.在冰體碰撞中，總的碰撞能分成了三個(gè)部分，一部分被泡沫鋁夾層板吸收變成變形能，一部分能量隨著冰體破碎被耗散，剩下的小部分能量轉(zhuǎn)化為楔形冰的反彈能.從能量的角度分析，在總的沖擊能量一定的情況下，冰體破碎耗散大部分沖擊能量，從而導(dǎo)致泡沫鋁夾層板吸收的能量比剛體碰撞少，產(chǎn)生的塑性變形也較剛體碰撞要小.由此可見，由于冰體碰撞中存在冰體破碎這一現(xiàn)象，它影響了泡沫鋁夾層板各個(gè)組分吸收能量的大小和比例，從而導(dǎo)致了泡沫鋁夾層板相較于剛體撞擊表現(xiàn)出不同的沖擊動(dòng)態(tài)響應(yīng).

圖11 泡沫鋁夾層板各個(gè)部分能量時(shí)程曲線

圖12 撞擊過(guò)程能量轉(zhuǎn)化曲線圖

3 抗沖擊性能參數(shù)影響研究

3.1 冰碰速度的影響

為了研究楔形冰碰撞速度對(duì)泡沫鋁夾層板的結(jié)構(gòu)動(dòng)態(tài)響應(yīng)的影響，針對(duì)圖3幾何尺寸，選取速度分別為2.0，2.5，3和3.5 m/s 進(jìn)行研究.不同速度下，泡沫鋁夾層板上、下面板中點(diǎn)處位移隨時(shí)間變化的曲線見圖13a)～b).由圖13a)～b)可知，上、下面板中點(diǎn)的位移隨著時(shí)間的變化規(guī)律基本相同，先隨著時(shí)間逐漸增大，直到達(dá)到峰值，然后夾層板開始釋放彈性能，位移值減小，直到楔形冰與夾層板發(fā)生分離，位移在很小的范圍內(nèi)上下震蕩.泡沫鋁夾層板上、下面板中點(diǎn)的最大變形和最終變形隨著碰撞速度變化曲線見圖13c).由圖13c)可知，上、下面板中點(diǎn)的最大變形和最終變形隨著碰撞速度的變化規(guī)律基本相同，隨著碰撞速度的增大而增大.同一速度下，下面板的變形比上面板的變形要小.這是因?yàn)樯厦姘搴托緦游樟舜蟛糠帜芰?，傳遞到下面板的能量較小，下面板受到的沖擊較小，所產(chǎn)生的變形也相應(yīng)較小.不同速度下，上、下面板的變形隨著速度的增大而增大.這是因?yàn)楸鏊俣仍黾訉?dǎo)致沖擊能量增加，夾層板各個(gè)部分吸收的能量隨之增加，從而導(dǎo)致變形增大.

圖13 結(jié)構(gòu)響應(yīng)與碰撞速度的關(guān)系

3.2 芯層厚度的影響

為了研究泡沫鋁夾層板芯層厚度對(duì)泡沫鋁夾層板的結(jié)構(gòu)動(dòng)態(tài)響應(yīng)的影響，針對(duì)圖3幾何尺寸，選取芯層厚度分別為20，25和30 mm進(jìn)行研究.不同芯層厚度下，泡沫鋁夾層板上、下面板中點(diǎn)處位移隨時(shí)間變化的曲線見圖14a)～b).由圖14a)～b)可知，上、下面板中點(diǎn)的位移隨著時(shí)間的變化規(guī)律基本相同，先隨著時(shí)間逐漸增大，直到達(dá)到峰值，然后夾層板開始釋放彈性能，位移值減小，直到楔形冰與夾層板發(fā)生分離，位移在很小的范圍內(nèi)上下震蕩.泡沫鋁夾層板上、下面板中點(diǎn)的最大變形和最終變形隨著芯層厚度變化曲線見圖14c).由圖14c)可知，上、下面板中點(diǎn)的最大變形和最終變形隨著芯層厚度的變化規(guī)律基本相同，隨著芯層厚度的增加而減小.當(dāng)芯層厚度逐漸增厚，可以看出泡沫鋁夾層板下面板變形減少量較上面板更為明顯，根據(jù)泡沫鋁準(zhǔn)靜態(tài)壓縮應(yīng)力-應(yīng)變曲線可知泡沫鋁具有平臺(tái)應(yīng)力階段，芯層越厚，可壓縮量越多，進(jìn)而吸收的能量越多，從而傳遞到下面板的能量越小，變形越小.由此可知，在面板厚度不變、沖擊能量不變等情況下，可以通過(guò)增加泡沫鋁夾層板芯層厚度來(lái)減小上、下面板的沖擊變形，同時(shí)也使其具備更優(yōu)秀的防護(hù)效果.

圖14 結(jié)構(gòu)響應(yīng)與芯層厚度的關(guān)系

4 結(jié) 論

1) 利用ANSYS/LS-DYNA材料庫(kù)中彈性斷裂模型可以較好的模擬冰在碰撞過(guò)程中的破碎現(xiàn)象，并與實(shí)驗(yàn)結(jié)果取得了較好的一致性.

2) 冰碰載荷作用下，泡沫鋁夾層板的上、下面板和泡沫鋁芯層有著不同的形變量，泡沫鋁芯層吸收能量的特性導(dǎo)致上、下面板的變形存在差別.

3) 冰體撞擊和剛體撞擊對(duì)泡沫鋁夾層板的上面板變形模態(tài)有差異，剛體撞擊泡沫鋁夾層板上面板有明顯的塑性鉸，而冰體撞擊泡沫鋁夾層板上面板無(wú)明顯塑性鉸.與剛體撞擊相比，冰體撞擊泡沫鋁夾層板由于冰體破碎會(huì)導(dǎo)致能量耗散，從而使得泡沫鋁夾層板吸收能量和產(chǎn)生的塑性變形都較前者低.