蔡正宇,丁圓圓,王士龍,鄭志軍,虞吉林

(1.中國科學(xué)技術(shù)大學(xué)近代力學(xué)系中國科學(xué)院材料力學(xué)行為和設(shè)計重點實驗室，安徽 合肥 230026;2.寧波大學(xué)力學(xué)與材料科學(xué)研究中心，浙江 寧波 315211)

梯度多胞犧牲層的抗爆炸分析*

蔡正宇1,丁圓圓2,王士龍1,鄭志軍1,虞吉林1

(1.中國科學(xué)技術(shù)大學(xué)近代力學(xué)系中國科學(xué)院材料力學(xué)行為和設(shè)計重點實驗室，安徽 合肥 230026;2.寧波大學(xué)力學(xué)與材料科學(xué)研究中心，浙江 寧波 315211)

運用一維非線性塑性沖擊波模型和細(xì)觀有限元模型對密度梯度多胞犧牲層的抗爆炸性能進(jìn)行了分析?；诼薀o關(guān)的剛性-塑性硬化模型，建立了描述沖擊波在多胞犧牲層中傳播的控制方程，分別給出了正、負(fù)密度梯度多胞材料在指數(shù)型爆炸載荷作用下的響應(yīng)特性。研究了可正好吸收爆炸能量的梯度多胞犧牲層的臨界厚度與載荷強度、覆蓋層質(zhì)量、多胞材料的密度梯度等參數(shù)之間的關(guān)系，給出了以臨界厚度和支撐端應(yīng)力峰值為指標(biāo)的密度梯度設(shè)計圖。運用二維細(xì)觀有限元模型驗證了基于非線性塑性沖擊波模型的抗爆炸分析的有效性。

多胞材料；犧牲層結(jié)構(gòu)；爆炸載荷；沖擊波；有限元方法

多胞材料具有優(yōu)異的能量吸收和抗沖擊性能，被廣泛用于汽車、列車、航空航天和軍事防護(hù)等領(lǐng)域[1]。梯度多胞金屬材料由于其可設(shè)計性受到廣泛關(guān)注，通過改變胞元尺寸、胞壁厚度等參數(shù)的分布情況引入密度梯度，使得材料呈現(xiàn)出漸進(jìn)的力學(xué)性能，可望獲得更優(yōu)異的能量吸收和抗沖擊性能[2]。多胞犧牲層是由覆蓋層和多胞芯層組成的夾芯結(jié)構(gòu)，通常放置在被保護(hù)結(jié)構(gòu)的表面。在爆炸載荷作用下，覆蓋層獲得較大的沖擊速度，芯層壓縮坍塌吸收大量的能量從而保護(hù)主體結(jié)構(gòu)不發(fā)生破壞[3]。在有限的空間下，芯層厚度需要設(shè)計得盡可能小，同時被保護(hù)結(jié)構(gòu)所承受的沖擊應(yīng)力要維持在許可應(yīng)力之下。最小芯層厚度和支撐端最小應(yīng)力是犧牲層作為防護(hù)結(jié)構(gòu)的兩個主要設(shè)計目標(biāo)。

已有研究表明，多胞材料的密度分布對結(jié)構(gòu)的沖擊響應(yīng)和能量吸收影響顯著[4]。密度分布會導(dǎo)致應(yīng)力不均勻性，在強動載荷作用下可能會在不同位置產(chǎn)生多個局部化變形帶。多數(shù)實驗研究局限于由均勻泡沫粘結(jié)而成的階梯式、不連續(xù)的梯度材料[5-6]。目前，對連續(xù)密度分布的多胞材料表現(xiàn)出的力學(xué)性能已有一定認(rèn)識。在具有線性密度分布的多胞材料抗沖擊行為方面，Wang Xiaokai等[7]認(rèn)為正梯度泡沫鋁可有效降低初始應(yīng)力峰值，Zheng Jie等[8]認(rèn)為在瞬態(tài)沖擊下，選用正、負(fù)梯度的泡沫金屬桿不同組合形式，能夠?qū)崿F(xiàn)保護(hù)沖擊端或者支撐端的功能。D.Karagiozova等[9]對比高速沖擊下多胞材料中沖擊波傳播的距離，認(rèn)為密度梯度多胞材料與均勻密度多胞材料相比在能量吸收方面沒有明顯優(yōu)勢。

多胞金屬是理想的吸能材料，常用于抗爆炸夾芯結(jié)構(gòu)的芯層。如果多胞犧牲層的芯層厚度設(shè)計得過大，吸收完爆炸能量后仍有未利用的區(qū)域；如果太小，材料完全壓實，沖擊波在支撐端反射后，出現(xiàn)應(yīng)力增強現(xiàn)象。因此，對芯層厚度進(jìn)行優(yōu)化設(shè)計以充分發(fā)揮其吸能能力顯得極為重要。在工程應(yīng)用中，最常見的是單層單一密度犧牲層，A.G.Hanssen等[10]基于剛性-理想塑性-鎖定(rigid,perfectly-plastic,locking,R-PP-L)模型[1]，給出了近距離爆炸載荷作用下芯層最小厚度的理論解。G.W.Ma等[11]引入雙層密度犧牲層到抗爆結(jié)構(gòu)中，發(fā)現(xiàn)低密度泡沫靠近爆炸端、高密度泡沫靠近支撐端的構(gòu)型可以吸收更多的能量，但是主體結(jié)構(gòu)承受應(yīng)力較大。采用與之相反的構(gòu)型，Liao Shenfei等[12]達(dá)到了降低支撐端應(yīng)力的目標(biāo)，給出了臨界厚度的設(shè)計方法。丁圓圓等[13]基于更精確的剛性-塑性硬化(rigid-plastic hardening, R-PH)模型[14]揭示了均勻密度多胞犧牲層中沖擊波的傳播規(guī)律。因此，改變密度的分布可以滿足更多的設(shè)計需求。

本文中，選用密度線性分布的梯度多胞金屬作為芯層材料，在相同爆炸載荷下，以芯層材料所需要的臨界厚度和支撐端應(yīng)力峰值為指標(biāo)，根據(jù)沖擊波在密度梯度多胞金屬中的傳播規(guī)律和犧牲層的響應(yīng)特性，考察不同梯度多胞金屬吸能的優(yōu)劣，以期為工程防護(hù)設(shè)計提供參考。

1 理論模型

1.1 問題描述

當(dāng)爆炸源到犧牲層的距離遠(yuǎn)大于被保護(hù)結(jié)構(gòu)的尺寸時，爆炸載荷產(chǎn)生的球面波可以近似為平面波[15]。最早由Taylor發(fā)展了爆炸脈沖入射到自由面上的一維沖擊波理論，單位面積上的爆炸載荷可以假設(shè)為一個指數(shù)衰減波，形如：

(1)

式中：t為時間，p0為初始峰值載荷，τ為衰減時間常數(shù)，易見爆炸載荷的沖量為p0τ。

以密度梯度多胞金屬材料為芯層，在爆炸載荷作用方向上密度漸進(jìn)增加的稱為正梯度犧牲層，反之則稱為負(fù)梯度犧牲層，見圖1。本文中選用密度沿加載方向線性分布的多胞金屬，其相對密度分布：

(2)

式中：γ為密度梯度參數(shù)，ρ0為平均相對密度，L0為芯層X方向上的厚度。γ取值的正、負(fù)分別對應(yīng)于正、負(fù)梯度材料。多胞犧牲層的覆蓋層假設(shè)為剛性質(zhì)量塊，面積質(zhì)量記為m，支撐端視為固定的。

多胞金屬的壓縮行為包含3個階段：彈性、平臺和壓實段。R-PP-L模型廣泛應(yīng)用于描述多胞材料的應(yīng)力應(yīng)變關(guān)系，以此發(fā)展了一維沖擊波理論[16]以研究多胞金屬的動態(tài)響應(yīng)。然而，R-PP-L模型高估了材料的壓潰應(yīng)力和能量吸收能力。Zheng Zhijun等[14]提出了率無關(guān)的R-PH模型以表征均勻泡沫金屬材料在準(zhǔn)靜態(tài)壓縮下的力學(xué)行為，其應(yīng)力應(yīng)變關(guān)系可以寫作:

(3)

式中:初始壓潰應(yīng)力σ0、應(yīng)變硬化參數(shù)C都是與相對密度有關(guān)的量[17]。

圖1 梯度多胞犧牲層的抗爆炸示意圖Fig.1 Schematic diagrams for the anti-blast analysis of graded cellular sacrificial claddings

1.2 正梯度犧牲層的抗爆炸分析模型

在爆炸載荷作用下，正梯度多胞金屬中沖擊波傳播的理論分析模型如圖1(a)所示。在初始構(gòu)形中建立坐標(biāo)系，沖擊波陣面的Lagrange坐標(biāo)為Φ(t)。根據(jù)R-PH模型假定，當(dāng)犧牲層承受到爆炸載荷時，彈性前驅(qū)波在沖擊端產(chǎn)生并以無窮大的速度向支撐端傳播，同時沖擊波也在沖擊端產(chǎn)生并以其特定的速度向前傳播，當(dāng)塑性沖擊波掃過某個位置時，質(zhì)點速度、應(yīng)變和應(yīng)力從波前的{0,0,σ0(ρ(Φ))}突變?yōu)椴ê蟮膡v(t),εb(t),σb(t)}，其中v為速度，σ0(ρ(Φ))為當(dāng)?shù)氐某跏級簼?yīng)力。

由應(yīng)力波理論，跨過沖擊波陣面的質(zhì)量守恒和動量守恒分別給出:

(4)

(5)

式中:ρs為基體材料密度。聯(lián)立可得波后應(yīng)力為:

(6)

另外，聯(lián)立式(3)和(6)可以求出波后應(yīng)變:

(7)

進(jìn)一步考察覆蓋層和沖擊波波后區(qū)域組成的部分，由牛頓運動定律，有:

(8)

結(jié)合上述方程，可以得到控制方程:

(9)

初始條件為Φ(0)=0和v(0)=0。這是非線性微分方程組，沒有顯式解，因此本文中采用四階Runge-Kutta法進(jìn)行數(shù)值求解。為了獲得犧牲層的臨界厚度Lc，先給定一個足夠大的厚度L0，若沖擊波停止的位置不在支撐端，則調(diào)整L0的值，直至沖擊波停止的位置正好在支撐端，此時的L0即為臨界厚度。

1.3 負(fù)梯度犧牲層的抗爆炸分析模型

多胞金屬在壓潰過程中，一旦局部應(yīng)力達(dá)到當(dāng)?shù)氐那姸染蜁a(chǎn)生局部變形。對于負(fù)梯度犧牲層，沖擊端的應(yīng)力急劇上升而發(fā)生逐層壓潰[18]，而不同于正梯度情形，支撐端也會出現(xiàn)沖擊波且朝著相反的方向傳播，如圖1(b)所示。這是由于多胞金屬密度低的位置強度低，負(fù)梯度犧牲層強度在沖擊方向上遞減。所以，支撐端達(dá)到當(dāng)?shù)爻跏級簼?yīng)力后出現(xiàn)坍塌并逐漸向強度高的地方擴展，可以看成一個向沖擊端傳播的塑性壓縮波。在負(fù)梯度犧牲層中形成了雙波結(jié)構(gòu)，分別記為Shock 1和Shock 2。這2個波陣面把犧牲層分成3部分：第1部分是質(zhì)量塊到Shock 1緊后方區(qū)域，在爆炸加載下速度急劇上升，之后隨著波后應(yīng)力增大而不斷減速，速度記為v1；第2部分是中間未變形區(qū)域，在其左右兩端的波前應(yīng)力差作用下，v2由0開始增加到與第1部分速度相同，最后減速直至犧牲層運動停止；第3部分為支撐端到Shock 2波后的材料，速度為0。

第1階段：0

(10)

而初始條件為Φ1(0)=0,v1(0)=0,Φ2(0)=L0以及v2(0)=0。

第2階段：t1

(11)

而定解條件由第1階段結(jié)束時給出。同樣采用四階Runge-Kutta法進(jìn)行數(shù)值求解。

根據(jù)該理論模型，首先給定L0一個足夠大的初值進(jìn)行響應(yīng)分析，若運動停止時犧牲層仍然有未變形區(qū)域，則調(diào)整L0的值重新計算，直至運動結(jié)束時Shock 2的Lagrange坐標(biāo)與運動第1階段Shock 1停止的位置重合，犧牲層得到了充分利用，此時選定的厚度L0即為臨界厚度Lc。

2 細(xì)觀有限元模型

考察具有不同相對密度的二維均勻Voronoi蜂窩[19]在準(zhǔn)靜態(tài)壓縮下的力學(xué)行為，如圖2(a)所示?；w材料密度ρs為2 700 kg/m3，材料滿足彈性-理想塑性假定，楊氏模量為69 GPa，泊松比為0.3，屈服強度為175 MPa。通過R-PH模型對名義應(yīng)力應(yīng)變曲線進(jìn)行擬合，可以確定其參數(shù)隨相對密度的變化特征。結(jié)果表明，初始壓潰應(yīng)力σ0和硬化參數(shù)C與相對密度近似呈二次方關(guān)系，如圖2(b)所示，擬合關(guān)系為：

(12)

式中：σys為基體材料的屈服應(yīng)力。

圖2 不同相對密度Voronoi蜂窩的準(zhǔn)靜態(tài)應(yīng)力應(yīng)變曲線和初始壓潰應(yīng)力和硬化參數(shù)的擬合結(jié)果Fig.2 Quasi-static stress-strain curves and fitting results of initial crushing stress parameter and strain hardening for Voronoi honeycombs with different relative densities

利用二維隨機Voronoi技術(shù)[19]和變胞元尺寸分布法[5]構(gòu)建密度梯度多胞金屬的細(xì)觀有限元模型。固定胞壁厚度h0為0.26 mm，樣本不規(guī)則度k取0.4，生成梯度Voronoi蜂窩的平均相對密度ρ0= 0.1。運用ABAQUS/Explicit進(jìn)行數(shù)值模擬，胞壁選用S4R殼單元，基體材料性質(zhì)同上所述。數(shù)值模擬情形為平面應(yīng)變狀態(tài)，即所有節(jié)點在面外方向上的位移限制為零，并且通過施加幾何約束使得在前后2個面上的對應(yīng)節(jié)點在面內(nèi)的運動保持一致，以避免殼單元在面外方向上發(fā)生扭曲。通過網(wǎng)格收斂性分析，單元的平均長度約為0.5 mm。所有接觸面定義為通用接觸，摩擦因數(shù)取為0.2。正、負(fù)梯度模型寬度均為80 mm，厚度方向取1 mm，長度取為理論模型得到的臨界厚度，分別為257.9和320.0 mm?？贡芏忍荻榷喟麪奚鼘拥挠邢拊Ｐ腿鐖D3所示。右端為靜止的剛性壁面，左端覆蓋剛性質(zhì)量塊，外表面作用指數(shù)型爆炸載荷，上下邊界為自由邊界。

圖3 抗爆炸密度梯度多胞犧牲層的細(xì)觀有限元模型Fig.3 Cell-based finite element (FE) models of graded cellular sacrificial claddings for blast mitigation

3 結(jié)果與討論

3.1 理論結(jié)果與有限元結(jié)果的比較

本節(jié)中爆炸載荷參數(shù)p0=20 MPa，τ=0.15 ms，覆蓋層m=5.4 kg/m2，多胞金屬平均相對密度為ρ0=0.1。

正、負(fù)梯度犧牲層波后質(zhì)點速度歷史曲線如圖4所示。正梯度犧牲層沖擊端速度在很短的時間內(nèi)達(dá)到最大值，隨后開始迅速下降，最后較緩慢地減小。臨界厚度為爆炸能量恰好被犧牲層吸收時多胞芯層的最小厚度，即為沖擊波掃過的距離。經(jīng)歷響應(yīng)時間tm=2.8 ms，犧牲層速度降為零。由圖可見，沖擊波理論模型預(yù)測結(jié)果與有限元計算結(jié)果吻合很好。沖擊波傳到的位置恰好為支撐端，以此設(shè)計出的臨界厚度Lc=257.9 mm。負(fù)梯度犧牲層的響應(yīng)比較復(fù)雜，中間未變形區(qū)域從靜止開始加速運動，與沖擊端獲得共同速度后緩慢地下降為零。在第2階段Shock 2到達(dá)Shock 1停止的位置，若恰好爆炸能量完全被吸收，設(shè)計出的臨界厚度Lc=320.0 mm，最終時刻tm=4.5 ms。由圖可見，理論預(yù)測與有限元計算結(jié)果吻合較好，但有限元模型得到的響應(yīng)時間要比理論模型得到的響應(yīng)時間長。爆炸載荷作用的初始階段，沖擊端處速度高，兩者都能觀察到連續(xù)的塑性壓潰面。正梯度犧牲層很規(guī)整地逐層坍塌，而負(fù)梯度犧牲層的壓潰面密實程度相對較低，后續(xù)的壓潰區(qū)域逐漸變得松散，表現(xiàn)出了變形局部化的特點。特別地，負(fù)梯度靠近支撐端附近可以明顯地觀察到胞元變形，起初并未完全坍塌閉合，隨著中間區(qū)域沖擊速度由零開始不斷增大，支撐端附近產(chǎn)生隨機剪切帶[14]，說明了壓縮波產(chǎn)生于支撐端附近。后期，沖擊速度不斷衰減，正梯度犧牲層靠近沖擊端也出現(xiàn)類似的變形模式，如圖5所示。

圖4 速度歷史曲線Fig.4 Time history of velocity

隨著爆炸載荷的作用，正梯度犧牲層由于其高密度端靠近被保護(hù)結(jié)構(gòu)，支撐端應(yīng)力隨著沖擊波陣面上透射波的增強而增強；而負(fù)梯度犧牲層支撐端應(yīng)力由Shock 2波后應(yīng)力確定，如圖6所示。正梯度吻合較好，負(fù)梯度理論預(yù)測的支撐端應(yīng)力變化趨勢與有限元計算結(jié)果一致，但是應(yīng)力值偏大。主要是由于理論模型采用連續(xù)體假設(shè)，認(rèn)為在響應(yīng)的起始階段支撐端立刻出現(xiàn)塑性沖擊波，而Voronoi模型表現(xiàn)為胞棱彎曲變形且胞元未完全坍塌壓實，因此理論解高估了支撐端應(yīng)力。如果支撐端應(yīng)力不大于許可應(yīng)力，這種基于臨界厚度的抗爆炸設(shè)計便是有效的。

圖5 變形圖Fig.5 Deformation patterns

圖6 支撐端應(yīng)力的歷史曲線Fig.6 Time history of stress at the support end

圖7 臨界厚度隨初始加載峰值的變化Fig.7 Critical thickness versus peak pressure of blast load

圖8 臨界厚度隨覆蓋層面積質(zhì)量的變化關(guān)系Fig.8 Critical thickness versus area mass of cover plate

3.2 臨界厚度的影響參數(shù)分析

本節(jié)運用理論模型分析爆炸載荷強度、覆蓋層質(zhì)量以及密度梯度變化對犧牲層臨界厚度的影響，為工程設(shè)計提供參考。

首先對于給定的密度梯度犧牲層，考察爆炸載荷對臨界厚度的影響，如圖7所示。結(jié)果表明在爆炸沖量p0τ保持不變時，臨界厚度Lc幾乎不變；而當(dāng)τ保持不變時，Lc隨著爆炸初始峰值近似地呈線性增加。同樣地，保持p0不變，Lc隨τ也會有類似的變化關(guān)系。由此可知，爆炸載荷是以沖量p0τ的形式影響Lc的大小。

其次，考慮梯度多胞犧牲層吸收相同爆炸載荷時，臨界厚度隨覆蓋層質(zhì)量的變化規(guī)律。結(jié)果表明，臨界厚度隨著覆蓋層質(zhì)量增大而減小，如圖8所示。這是由于爆炸載荷是以沖量形式加載的，覆蓋層質(zhì)量的大小決定了爆炸載荷轉(zhuǎn)化為動能的能力，即質(zhì)量越大，由爆炸載荷所做的功越小。在設(shè)計允許范圍內(nèi)，增大覆蓋層質(zhì)量是有利的。

最后，考察密度梯度參數(shù)γ對多胞犧牲層臨界厚度的影響。由于使用了相同的平均相對密度的芯層，因此吸收同等強度的爆炸能量下，犧牲層臨界厚度的大小直接反映需要材料的多少,如圖9所示。

在圖9(a)中，虛線對應(yīng)于均勻密度多胞芯層的臨界厚度。顯然，正梯度芯層需要的厚度小于均勻密度多胞芯層的臨界厚度，而負(fù)梯度芯層臨界厚度相對較大。在工程上選用正梯度可以設(shè)計出更薄、更輕的犧牲層。另外，γ表征了密度分布的陡峭程度，隨著其絕對值的增加，材料局部密度的差異增大。正、負(fù)梯度犧牲層的臨界厚度都隨著|γ|的增大而減小，因此選取更大的密度梯度參數(shù)絕對值，對于減小芯層厚度也是有利的。

圖9 臨界厚度和支撐端應(yīng)力峰值隨密度梯度參數(shù)的變化Fig.9 Variations of critical thickness and peak stress at the support end with the gradient parameter

支撐端應(yīng)力峰值σm隨γ的變化關(guān)系如圖9(b)所示。對于相同的|γ|，選用正梯度芯層帶來的應(yīng)力峰值更大。正梯度芯層應(yīng)力峰值出現(xiàn)在響應(yīng)結(jié)束時刻，即沖擊波運動至支撐端、應(yīng)力達(dá)到密度最大位置對應(yīng)的初始壓潰應(yīng)力，它和覆蓋層質(zhì)量的選取無關(guān)并且與γ近似呈二次方關(guān)系。負(fù)梯度多胞芯層的應(yīng)力峰值為沖擊波波后應(yīng)力的最大值，與不同運動階段的響應(yīng)有關(guān)。隨著|γ|增大，應(yīng)力峰值起初維持在很低的水平，當(dāng)|γ|大于某個特定值之后，應(yīng)力峰值逐漸上升，例如m取5.4 kg/m2時，應(yīng)力峰值為0.53 MPa。同時，覆蓋層對負(fù)梯度芯層的支撐端應(yīng)力也有影響，質(zhì)量越大，應(yīng)力峰值越大。

進(jìn)一步分析可知，支撐端應(yīng)力峰值和臨界厚度這2個設(shè)計指標(biāo)之間呈負(fù)相關(guān)，即在工程設(shè)計中要獲得較低的應(yīng)力峰值必然要增加芯層厚度。如果希望被保護(hù)結(jié)構(gòu)承受應(yīng)力最小，選用|γ|較小的負(fù)梯度或者均勻密度效果更好，在受限的空間下，均勻密度多胞芯層最優(yōu)。當(dāng)給定覆蓋層質(zhì)量并限定支撐端應(yīng)力峰值時，選擇適當(dāng)?shù)恼荻葏?shù)可得到最小的犧牲層厚度。因此，圖9可為密度梯度多胞犧牲層的抗爆炸設(shè)計提供依據(jù)。

4 結(jié) 論

研究了密度梯度多胞犧牲層的抗爆炸性能，采用率無關(guān)的剛性-塑性硬化模型描述多胞金屬材料的應(yīng)力應(yīng)變關(guān)系，對于正梯度和負(fù)梯度情形分別建立了一維非線性塑性沖擊波模型，揭示了塑性沖擊波在梯度多胞金屬中的傳播規(guī)律，得到了梯度多胞犧牲層的響應(yīng)特性，進(jìn)而給出了芯層臨界厚度的設(shè)計方法?？紤]爆炸載荷強度、覆蓋層質(zhì)量和密度梯度參數(shù)對臨界厚度的影響，發(fā)現(xiàn)可以通過增大覆蓋層質(zhì)量和密度梯度參數(shù)設(shè)計出更小的臨界厚度。以臨界厚度和支撐端應(yīng)力峰值為指標(biāo)，對比不同梯度多胞金屬吸收相同爆炸載荷能力的優(yōu)劣。結(jié)果表明，通過增大梯度參數(shù)值，正梯度能夠有效地降低芯層厚度，但是會帶來支撐端應(yīng)力峰值的增大；如果要實現(xiàn)應(yīng)力峰值最小，同時考慮盡可能小的厚度，均勻密度是較理想的選擇。采用細(xì)觀有限元模型檢驗了犧牲層抗爆炸分析的可靠性，理論預(yù)測結(jié)果與有限元結(jié)果吻合較好，說明基于非線性塑性硬化沖擊波模型設(shè)計梯度多胞犧牲層是有效的。

[1] Reid S R, Peng C. Dynamic uniaxial crushing of wood[J]. International Journal of Impact Engineering, 1997,19(5/6):531-570.

[2] Liu Jiagui, Hou Bing, Lu Fangyun, et al. A theoretical study of shock front propagation in the density graded cellular rods[J]. International Journal of Impact Engineering, 2015,80:133-142.

[3] Guruprasad S, Mukherjee A. Layered sacrificial claddings under blast loading: Part I: Analytical studies[J]. International Journal of Impact Engineering, 2000,24(9):957-973.

[4] Cui L, Kiernan S, Gilchrist M D. Designing the energy absorption capacity of functionally graded foam materials[J]. Materials Science and Engineering A, 2009,507(1):215-225.

[5] Zhang Jianjun, Wang Zhihua, Zhao Longmao. Dynamic response of functionally graded cellular materials based on the Voronoi model[J]. Composites Part B: Engineering, 2016,85:176-187.

[6] 吳鶴翔,劉穎.梯度變化對密度梯度蜂窩材料力學(xué)性能的影響[J].爆炸與沖擊,2013,33(2):163-168. Wu Hexiang, Liu Ying. Influences of density gradient variation on mechanical performances of density graded honeycomb materials[J]. Explosion and Shock Waves, 2013,33(2):163-168.

[7] Wang Xiaohai, Zheng Zhijun, Yu Jilin. Crashworthiness design of density-graded cellular metals[J]. Theoretical and Applied Mechanics Letters, 2013,3(3):9-13.

[8] Zheng Jie, Qin Qinghua, Wang T J. Impact plastic crushing and design of density-graded cellular materials[J]. Mechanics of Material, 2016,94:66-78.

[9] Karagiozova D, Alves M. Propagation of compaction waves in cellular materials with continuously varying density[J]. International Journal of Solids and Structures, 2015,71:323-337.

[10] Hanssen A G, Enstock L, Langseth M. Close-range blast loading of aluminum foam panels[J]. International Journal of Impact Engineering, 2002,27(6):593-618.

[11] Ma G W, Ye Z Q. Energy absorption of double-layer foam cladding for blast alleviation[J]. International Journal of Impact Engineering, 2007,34(2):329-347.

[12] Liao Shenfei, Zheng Zhijun, Yu Jilin, et al. A design guide of double-layer cellular claddings for blast alleviation[J]. International Journal of Aerospace and Lightweight Structures, 2013,3(1):109-133.

[13] Sawle D R. Hypervelocity impact in thin sheets, semi-infinite targets at 15 km/s[J]. AIAA Journal, 1970,8(7):1240-1244.

[14] 丁圓圓,王士龍,鄭志軍,等.多胞犧牲層的抗爆炸分析[J].力學(xué)學(xué)報,2014,46(6):825-833. Ding Yuanyuan, Wang Shilong, Zheng Zhijun, et al. Anti-blast analysis of cellular sacrificial cladding[J]. Chinese Journal of Theoretical and Applied Mechanics, 2014,46(6):825-833.

[15] Fleck N A, Deshpande V S. The resistance of clamped sandwich beams to shock loading[J]. Journal of Applied Mechanics, 2004,71(3):386-401.

[16] Tan P J, Reid S R, Harrigan J J, et al. Dynamic compressive strength properties of aluminium foams: Part II: Shock theory and comparison with experimental data and numerical models[J]. Journal of the Mechanics and Physics of Solids, 2005,53(10):2206-2230.

[17] Wang Shilong, Ding Yuanyuan, Wang Changfeng, et al. Dynamic material parameters of closed-cell foams under high-velocity impact[J]. International Journal of Impact Engineering, 2017,99:111-121.

[18] 王長峰,鄭志軍,虞吉林.泡沫桿撞擊剛性壁的動態(tài)壓潰模型[J].爆炸與沖擊,2013,33(6):587-593. Wang Changfeng, Zheng Zhijun, Yu Jilin. Dynamic crushing models for a foam rod striking a rigid wall[J]. Explosion and Shock Waves, 2013,33(6):587-593.

[19] Zheng Zhijun, Yu Jilin, Li Jianrong. Dynamic crushing of 2D cellular structures: A finite element study[J]. International Journal of Impact Engineering, 2005,32(1/2/3/4):650-664.

(責(zé)任編輯 張凌云)

Anti-blast analysis of graded cellular sacrificial cladding

Cai Zhengyu1, Ding Yuanyuan2, Wang Shilong1, Zheng Zhijun1, Yu Jilin1

(1.CASKeyLaboratoryofMechanicalBehaviorandDesignofMaterials,DepartmentofModernMechanics,UniversityofScienceandTechnologyofChina,Hefei230026,Anhui,China;2.MechanicsandMaterialsScienceResearchCenter,NingboUniversity,Ningbo315211,Zhejiang,China)

The blast mitigation behavior of a density-graded cellular sacrificial cladding is investigated by using a nonlinear plastic shock model and a cell-based finite element model. Based on a rate-independent, rigid-plastic hardening idealization, a theoretical approach is applied to analyze the propagation of shock wave in density-graded cellular rods subjected to blast loading. The influences of the intensity of blast load, the cover mass and the density gradient parameter of the cellular material on the critical thickness, which is the minimum thickness of the core layer when the energy of explosion is fully absorbed, are investigated. A design guide of density-gradient is provided which considers the critical thickness of the cellular core as well as the peak stress at the support end. The validity of the anti-blast analysis of the graded cellular sacrificial cladding based on the nonlinear plastic shock model is verified by using cell-based finite element models.

cellular materials; sacrificial cladding; blast load; shock wave; finite element method

10.11883/1001-1455(2017)03-0396-09

2017-02-20；

2017-04-28

國家自然科學(xué)基金項目(11372307,11372308)

蔡正宇(1990— )，男，碩士研究生; 通信作者: 鄭志軍，zjzheng@ustc.edu.cn。

O382 國標(biāo)學(xué)科代碼： 13035

A