吳 雄，李應(yīng)剛,3，肖 雯，蔡 偉，朱 凌

(1.高性能艦船技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室(武漢理工大學(xué))，武漢 430063；2.武漢理工大學(xué) 船海與能源動(dòng)力工程學(xué)院，武漢 430063；3.西安交通大學(xué) 機(jī)械結(jié)構(gòu)強(qiáng)度與振動(dòng)國(guó)家重點(diǎn)實(shí)驗(yàn)室，西安 710049)

隨著全球氣候變暖和極地海冰的快速消融，極地地區(qū)蘊(yùn)藏的豐富資源和極具戰(zhàn)略意義的航道價(jià)值日益凸顯。然而，極地洋面上漂浮著大量浮冰和小型冰山，極地船舶與海洋工程裝備在極地航行與工程作業(yè)過程中不可避免遭受浮冰碰撞，造成船體結(jié)構(gòu)損傷，威脅船舶航行安全及船員的生命。因此，浮冰碰撞下極地船舶與海洋工程裝備安全性顯得尤為重要[1-2]。

Liu等[3-5]較系統(tǒng)地研究了船舶與冰山碰撞的外部動(dòng)力學(xué)及內(nèi)部動(dòng)力學(xué)問題，提出了考慮到船體幾何形狀構(gòu)造和摩擦因數(shù)對(duì)能量損耗影響的計(jì)算公式。張健等[6]利用數(shù)值仿真的方法模擬了船-冰碰撞過程，研究并總結(jié)了冰體形狀、冰體質(zhì)量、碰撞速度等對(duì)船舶碰撞力、船舶損傷變形、能量吸收等因素的影響規(guī)律。Ritch等[7]進(jìn)行了Terry Fox破冰船與小型冰山碰撞試驗(yàn)，測(cè)量了船體結(jié)構(gòu)與小型冰山碰撞過程中的冰載荷。王林等[8]通過數(shù)值仿真建立了各向同性黏塑性冰層有限元模型，計(jì)算并比較了兩種典型破冰結(jié)構(gòu)的破冰機(jī)理和破冰效果,分析了船舶船首與冰層碰撞時(shí)的損傷特點(diǎn)。劉俊杰等[9]采用數(shù)值方法模擬了水面艦船在冰區(qū)航行遭受浮冰過程，得到了碰撞過程中艦船的動(dòng)態(tài)響應(yīng)及能量轉(zhuǎn)化特性。Gao等[10]采用各向同性彈塑性冰材料模型分析了船舶與冰山碰撞過程中的碰撞力和能量耗散。王健偉等[11]對(duì)船-冰碰撞進(jìn)行了數(shù)值模擬，研究了船舶在不同速度下與不同厚度冰層碰撞的動(dòng)態(tài)結(jié)構(gòu)響應(yīng)。Cai等[12-13]開發(fā)了混凝土冰材料數(shù)值模型，采用非線性有限元軟件ANSYS/LS-DYNA研究了浮冰碰撞下船體板結(jié)構(gòu)彈塑性動(dòng)態(tài)響應(yīng)特性，并開展了船-冰碰撞試驗(yàn)驗(yàn)證。

輕質(zhì)蜂窩夾芯結(jié)構(gòu)作為一種物理功能與結(jié)構(gòu)一體化結(jié)構(gòu)，具有高比強(qiáng)度、高比剛度、高韌性、高能量吸收等優(yōu)良的力學(xué)性能，在極地船舶與海洋工程結(jié)構(gòu)安全與碰撞沖擊防護(hù)領(lǐng)域具有廣泛的潛在應(yīng)用價(jià)值[14-15]。然而，目前對(duì)蜂窩夾芯結(jié)構(gòu)的動(dòng)力學(xué)行為研究主要集中在剛體沖擊，浮冰碰撞下蜂窩夾芯結(jié)構(gòu)動(dòng)態(tài)響應(yīng)特性尚不清楚。本文采用數(shù)值模擬與試驗(yàn)方法研究浮冰碰撞下蜂窩金屬夾芯板動(dòng)態(tài)響應(yīng)特性，揭示其塑性變形能量吸收機(jī)理與冰體能量耗散機(jī)制，分析浮冰碰撞沖擊位置及蜂窩芯層厚度對(duì)其動(dòng)態(tài)變形及能量分配影響規(guī)律。

1 數(shù)值仿真

1.1 數(shù)值仿真模型

由于極地浮冰的形狀一般呈不規(guī)則的多邊狀棱柱，且Liu等指出楔形浮冰對(duì)船舶損傷最為嚴(yán)重，因此本文選取楔形浮冰作為研究對(duì)象進(jìn)行試驗(yàn)和數(shù)值仿真。采用非線性有限元軟件ANSYS/LS-DYNA建立了楔形冰碰撞蜂窩金屬夾芯板有限元模型，如圖1(a)所示。蜂窩金屬夾芯板由上、下面板和正六邊形蜂窩芯層組成，上、下面板的邊長(zhǎng)為400 mm，厚度為1 mm，蜂窩芯層厚度為20 mm，正六邊形蜂窩邊長(zhǎng)為6 mm，單層壁厚為0.07 mm，如圖1(b)所示。蜂窩金屬夾芯板上、下面板采用A5052鋁合金材料，蜂窩芯層采用A3003-H19鋁合金材料。上、下面板和蜂窩芯層采用Shell 163單元建模，定義為分段線性彈塑性材料模型(*MAT_024)，通過萬(wàn)能試驗(yàn)機(jī)進(jìn)行準(zhǔn)靜態(tài)拉伸試驗(yàn)測(cè)試得到材料力學(xué)特性參數(shù)，如表1所示。楔形浮冰采用Solid 164單元建模，定義為混凝土冰材料模型(*MAT_078)，具體材料參數(shù)如表2所示。本文重點(diǎn)關(guān)注浮冰-蜂窩金屬夾芯板瞬態(tài)碰撞沖擊內(nèi)部動(dòng)力學(xué)行為過程，忽略流體作用及外部動(dòng)力學(xué)的影響，楔形冰碰撞沖擊過程中只保留z方向上自由度，約束其余的自由度，夾芯板上、下面板四周剛性固定約束，如圖1(a)所示。上、下面板與蜂窩芯層之間定義固-連接觸(*CONTACT_TIED_NODES_TO_SURFACE)，模擬面板與蜂窩芯層的粘接，上面板與楔形浮冰之間定義面-面侵蝕接觸(*CONTACT_ERODING_SURFACE_TO_SURFACE)，模擬碰撞過程中楔形冰破碎失效，由于浮冰碰撞沖擊過程中蜂窩芯層可能會(huì)出現(xiàn)褶皺，定義蜂窩芯層自接觸(*CONTACT_AUTOMATIC_SINGLE_SURFACE)。楔形冰和上、下面板的網(wǎng)格尺寸為5 mm，蜂窩芯層的網(wǎng)格尺寸為2 mm。為了保證計(jì)算精度，對(duì)沖擊區(qū)域進(jìn)行網(wǎng)格細(xì)化，細(xì)化網(wǎng)格為0.5 mm。初始時(shí)刻楔形冰與蜂窩夾芯板相距50 mm，賦予楔形冰3 m/s的初始沖擊速度，楔形冰的初始沖擊能量為360 J。采用ANSYS/LS-DYNA顯式動(dòng)力學(xué)中心差分方法對(duì)運(yùn)動(dòng)方程進(jìn)行顯式時(shí)間積分，即可求解得到浮冰碰撞下蜂窩金屬夾芯板結(jié)構(gòu)動(dòng)態(tài)響應(yīng)與能量吸收特性。

圖1 結(jié)構(gòu)示意圖Fig.1 Structure diagram

表1 面板和蜂窩芯層的材料參數(shù)(*MAT_024)Tab.1 Material parameters of face sheet and honeycomb core (*MAT_024)

表2 楔形冰材料模型參數(shù)(*MAT_078)Tab.2 Material parameters of ice wedge(*MAT_078)

1.2 數(shù)值仿真計(jì)算結(jié)果

圖2(a)為楔形浮冰碰撞下蜂窩金屬夾芯板動(dòng)態(tài)響應(yīng)過程，當(dāng)楔形浮冰與蜂窩金屬夾芯板碰撞后，楔形浮冰的破碎主要發(fā)生在楔形冰頭前端高壓應(yīng)力區(qū)，而其他區(qū)域基本沒有發(fā)生破壞。圖2(b)為上、下面板在碰撞過程中的變形模式，從圖2(b)可看出上面板在楔形浮冰碰撞下表現(xiàn)為局部凹陷與整體彎曲的耦合變形模式，而下面板主要表現(xiàn)為整體彎曲變形模式。由圖可知，有限元數(shù)值模型可以較好的模擬蜂窩夾芯結(jié)構(gòu)動(dòng)態(tài)變形與浮冰碰撞破碎失效過程。

圖2 楔形浮冰碰撞下蜂窩金屬夾芯板動(dòng)態(tài)響應(yīng)過程Fig.2 Dynamic impact process of AHSP under ice floe impact

為了進(jìn)一步深入分析蜂窩金屬夾芯板在楔形浮冰碰撞下動(dòng)態(tài)響應(yīng)與能量吸收特性，求解得到碰撞過程的沖擊力時(shí)程曲線、沖擊力-位移曲線、撓度時(shí)程曲線、能量分配曲線如圖3所示。由圖3可以看到，楔形浮冰碰撞沖擊力時(shí)程曲線呈半波正弦脈沖狀，沖擊力峰值達(dá)到27.5 kN。上、下面板在楔形浮冰碰撞下的變形逐漸增大，由于蜂窩芯層局部壓縮，上面板撓度大于下面板撓度，上面板中點(diǎn)的最大撓度為19 mm。當(dāng)楔形浮冰碰撞結(jié)束后，由于上、下面板材料彈性效應(yīng)出現(xiàn)回彈，上面板中點(diǎn)的最終撓度為15.5 mm。由圖3(d)可知，楔形浮冰碰撞過程的沖擊能量分配與剛體碰撞顯著不同，楔形浮冰碰撞沖擊能量主要轉(zhuǎn)化為蜂窩金屬夾芯板塑性變形能和楔形冰回彈動(dòng)能以及冰體破碎耗散能量，其中大部分能量被上面板和蜂窩芯層吸收以及冰體破碎耗散，下面板產(chǎn)生塑性變形較小，能量吸收少，說明蜂窩金屬夾芯板具有優(yōu)良的沖擊防護(hù)性能。

圖3 楔形浮冰碰撞下蜂窩金屬夾芯板動(dòng)態(tài)響應(yīng)特性Fig.3 Dynamic responses of honeycomb sandwich plates under ice floe impact

2 試驗(yàn)驗(yàn)證

為了驗(yàn)證浮冰碰撞下蜂窩金屬夾芯板動(dòng)態(tài)響應(yīng)數(shù)值仿真模型的可靠性，我們采用水平?jīng)_擊試驗(yàn)裝置進(jìn)行了楔形冰碰撞蜂窩金屬夾芯板動(dòng)態(tài)響應(yīng)試驗(yàn)，如圖4(a)所示。水平?jīng)_擊試驗(yàn)裝置包括安裝蜂窩金屬夾芯板的基座、夾具、沖擊小車、軌道。試驗(yàn)過程中制備蜂窩金屬夾芯板尺寸為560 mm×560 mm，蜂窩金屬夾芯板固定在基座上實(shí)現(xiàn)固支邊界條件，蜂窩金屬夾芯板的有效尺寸為400 mm×400 mm，如圖4(b)所示，其他材料參數(shù)和幾何參數(shù)與數(shù)值仿真模型保持一致。采用碎冰機(jī)處理得到5～10 mm尺寸的碎冰，利用楔形冰模具填充碎冰并加滿水，然后在放置在恒溫-25 ℃的工業(yè)冰箱中3天以上即可制備得到楔形冰沖頭如圖4(c)所示，楔形冰頭幾何參數(shù)與數(shù)值仿真模型保持一致。楔形冰放置在沖擊小車中，沖擊小車沿著軌道上升到指定高度后釋放即可獲得初始沖擊速度與沖擊能量，實(shí)現(xiàn)浮冰碰撞蜂窩金屬夾芯板試驗(yàn)。

圖4(d)為碰撞后的蜂窩夾芯板，從圖中可看出由于楔形冰頭碰撞沖擊作用，蜂窩金屬夾芯板上面板產(chǎn)生了局部凹陷。圖4(e)為碰撞后楔形冰損傷破碎圖，從圖中可以發(fā)現(xiàn)冰體破碎只發(fā)生在楔形冰頭前端區(qū)域。楔形浮冰碰撞蜂窩金屬夾芯板動(dòng)態(tài)響應(yīng)試驗(yàn)過程中，通過高速攝影測(cè)得楔形冰碰撞蜂窩金屬夾芯板的瞬時(shí)速度為3 m/s，通過M+P數(shù)據(jù)采集儀和B &K加速度傳感器測(cè)得碰撞過程的加速度-時(shí)間曲線，并通過公式轉(zhuǎn)化為沖擊力-時(shí)間曲線。試驗(yàn)中楔形冰質(zhì)量為18.02 kg，小車與配重質(zhì)量為61.98 kg，總質(zhì)量為80 kg，楔形冰與小車的初始沖擊動(dòng)能為360 J，與數(shù)值模型沖擊動(dòng)能保持一致。楔形冰碰撞蜂窩金屬夾芯板動(dòng)態(tài)響應(yīng)試驗(yàn)結(jié)束后，通過游標(biāo)卡尺測(cè)量上、下面板中線各個(gè)位置的最終撓度，并與仿真計(jì)算結(jié)果進(jìn)行對(duì)比。由圖5(a)可知，上、下面板中線最終撓度數(shù)值仿真結(jié)果與試驗(yàn)測(cè)試結(jié)果吻合較好，從圖5(b)可看出，數(shù)值計(jì)算得到的沖擊力-時(shí)間曲線與試驗(yàn)結(jié)果吻合較好，驗(yàn)證了數(shù)值仿真模型的可靠性。浮冰碰撞沖擊試驗(yàn)過程中沖擊時(shí)間相比數(shù)值仿真的沖擊時(shí)間長(zhǎng)，分析誤差產(chǎn)生原因主要來源于試驗(yàn)過程中邊界條件的彈性效應(yīng)。

圖4 楔形冰碰撞蜂窩金屬夾芯板試驗(yàn)(mm)Fig.4 Experiment test of honeycomb sandwich plate under ice floe impact(mm)

圖5 數(shù)值計(jì)算結(jié)果與試驗(yàn)結(jié)果對(duì)比Fig.5 Comparisons of numerical results with experimental results

3 參數(shù)分析

3.1 浮冰碰撞位置的影響

為了分析碰撞位置對(duì)楔形浮冰碰撞下蜂窩金屬夾芯板動(dòng)態(tài)響應(yīng)的影響，保持?jǐn)?shù)值仿真模型幾何參數(shù)和材料參數(shù)不變，計(jì)算楔形浮冰與蜂窩金屬夾芯板在3種碰撞位置下的碰撞力-時(shí)間曲線、碰撞力-位移曲線、上、下面板撓度和蜂窩芯層壓縮量曲線以及能量吸收率曲線，結(jié)果如圖6所示，其中3種浮冰碰撞沖擊位置在圖6(a)中標(biāo)出，由于楔形浮冰碰撞沖擊位置改變導(dǎo)致浮冰與蜂窩金屬夾芯板接觸剛度變化，蜂窩金屬夾芯板中間的位置1處接觸剛度最小，靠近剛性固定邊界的位置2和位置3處接觸剛度較大。由圖6(a)和圖6(b)可知，楔形浮冰碰撞沖擊位置對(duì)碰撞沖擊力時(shí)程曲線有顯著影響，3種碰撞沖擊位置下碰撞力峰值分別為29.9 kN，28.6 kN，25.9 kN，對(duì)應(yīng)的碰撞持續(xù)時(shí)間分別為18.0 ms，17.5 ms，16.0 ms。由圖6(c)可知，隨著楔形浮冰碰撞沖擊位置改變，上面板中心點(diǎn)撓度差異不大，下面板中心點(diǎn)撓度和蜂窩芯層壓縮量呈現(xiàn)相反的變化規(guī)律。由圖6(d)可知，隨著楔形浮冰碰撞沖擊位置改變，蜂窩金屬夾芯板上、下面板能量吸收率影響較小，蜂窩芯層能量吸收率與浮冰破碎耗散能量呈現(xiàn)相反的變化規(guī)律。綜上所述，當(dāng)碰撞沖擊位置接觸剛度較大時(shí)，楔形浮冰破碎失效越嚴(yán)重，碰撞力峰值越小，碰撞持續(xù)時(shí)間越短，冰體破碎耗散能量吸收率越大，冰體破碎能量耗散對(duì)于碰撞沖擊能量吸收起到越來越重要的作用。

圖6 浮冰碰撞沖擊位置對(duì)蜂窩金屬夾芯板動(dòng)態(tài)響應(yīng)的影響Fig.6 Effect of impact position on the dynamic response of honeycomb sandwich plate

3.2 蜂窩芯層厚度的影響

為了分析蜂窩芯層厚度對(duì)楔形浮冰碰撞下蜂窩金屬夾芯板動(dòng)態(tài)響應(yīng)的影響，保持?jǐn)?shù)值仿真模型其他幾何參數(shù)、材料參數(shù)以及沖擊參數(shù)不變，分別計(jì)算蜂窩芯層厚度分別為10 mm，20 mm，30 mm時(shí)蜂窩金屬夾芯板的碰撞力-時(shí)間曲線、碰撞力-位移曲線、上、下面板撓度和蜂窩芯層壓縮量曲線以及能量吸收率曲線，計(jì)算結(jié)果如圖7所示。由圖7(a)和圖7(b)中可知，蜂窩芯層厚度為10 mm，20 mm，30 mm時(shí)的碰撞力分別為27.4 kN，29.9 kN，29.4 kN，對(duì)應(yīng)的碰撞持續(xù)時(shí)間分別為20.0 ms，18.0 ms，16.5 ms，隨著蜂窩芯層厚度的增加，碰撞力的峰值基本保持不變，碰撞持續(xù)時(shí)間顯著縮短。由圖7(c)和圖7(d)可知，隨著蜂窩芯層厚度的增加，上、下面板的中心點(diǎn)撓度減小，蜂窩芯層壓縮量逐漸增大。上、下面板的能量吸收率逐漸下降，蜂窩芯層能量吸收率逐漸上升，楔形浮冰破碎失效增多，導(dǎo)致冰破碎耗散能量增加。綜上可知，蜂窩芯層厚度越大，上、下面板吸能作用越來越小，蜂窩芯層吸能效果越來越明顯，冰破碎耗散能量增加。

圖7 蜂窩芯層厚度對(duì)蜂窩金屬夾芯板動(dòng)態(tài)響應(yīng)的影響Fig.7 Effect of honeycomb core thickness on the dynamic response of honeycomb sandwich plate

4 結(jié) 論

本文通過非線性有限元軟件ANSYS/LS-DYNA建立了楔形浮冰碰撞下蜂窩金屬夾芯板動(dòng)態(tài)響應(yīng)有限元模型，并開展了楔形浮冰碰撞蜂窩金屬夾芯板動(dòng)態(tài)響應(yīng)試驗(yàn)驗(yàn)證，分析了碰撞沖擊位置和蜂窩芯層厚度對(duì)蜂窩金屬夾芯板動(dòng)態(tài)響應(yīng)的影響，可得出以下結(jié)論：

(1)楔形冰碰撞下蜂窩金屬夾芯板上面板表現(xiàn)為局部凹陷與整體彎曲的耦合變形模式，下面板表現(xiàn)為整體彎曲變形模式，沖擊能量轉(zhuǎn)化為蜂窩金屬夾芯板的變形能和楔形冰的回彈動(dòng)能以及冰體破碎耗散能量，數(shù)值仿真與試驗(yàn)結(jié)果吻合較好，驗(yàn)證了數(shù)值計(jì)算模型的準(zhǔn)確性。

(2)隨著楔形浮冰碰撞沖擊位置的接觸剛度逐漸增大，蜂窩金屬夾芯板塑性變形能量吸收率逐漸降低，冰體破碎耗散能量吸收率逐漸增大，冰體破碎能量耗散對(duì)于碰撞沖擊能量吸收起到越來越重要的作用。隨著蜂窩芯層厚度的增加，上、下面板的能量吸收率逐漸下降，蜂窩芯層能量吸收率逐漸上升，楔形浮冰破碎耗散能量逐漸增加。