梯度蜂窝夹芯板在爆炸荷载作用下的动力响应*

李世强，李 鑫，吴桂英，王志华，赵隆茂

(太原理工大学应用力学与生物医学工程研究所，山西 太原 030024)

梯度蜂窝夹芯板在爆炸荷载作用下的动力响应*

李世强，李 鑫，吴桂英，王志华，赵隆茂

(太原理工大学应用力学与生物医学工程研究所，山西 太原 030024)

利用弹道冲击摆锤系统对分层梯度蜂窝夹芯板在爆炸荷载下的动力响应进行了实验研究，分析了梯度蜂窝夹芯板在爆炸荷载作用下的变形失效模式，并与传统非梯度蜂窝夹芯板的抗爆性能做了对比。通过一维应力波理论，分析了应力波在梯度芯层中的传播规律。应力波透射系数在梯度试件中比非梯度芯层中小，而且相对密度递减的芯层组合有最小的应力波透射系数。综合考虑结构变形失效模式，后面板挠度，芯层压缩量以及应力波传播特点得到：分层梯度蜂窝夹芯板的抗爆性能明显优于传统的非梯度夹芯板，在所研究的荷载范围内，芯层相对密度从大到小排列试件的抗爆性能相对较好。

爆炸力学；动力响应；爆炸荷载；蜂窝夹芯板；功能梯度材料

随着汽车工业及航空航天技术的不断发展，夹芯复合结构由于比强度和比刚度高而被广泛的应用于防撞吸能构件的设计中。研究人员做了许多有关夹芯结构准静态[1]、动态[2]力学行为的工作，涉及的芯层结构包括蜂窝[3]、金属泡沫[4]、波纹板[5]及点阵结构[6]等多种形式，但大多数研究限于非梯度夹芯结构的力学行为。

分层梯度夹芯结构被认为是提高夹芯结构强度，减轻重量更为有效的途径。E.Etemadi等[7]利用LS-DYNA有限元程序分析了功能梯度芯梁在圆柱形子弹冲击下的力学行为，他们发现与传统均质梁相比，最大接触力和最大应变都有明显的降低。L.Cui等[8]提出了一种关于功能梯度聚合物泡沫的理论模型，在此模型的基础上，利用有限元方法研究了功能梯度泡沫的吸能特性，他们指出：增加梯度层之间的密度差异会降低芯层的吸能特性，相反减小层间密度差异能够提高芯层的吸能效率。N.Gardner等[9]分析了功能梯度夹芯梁在冲击波作用下的力学响应，发现增加芯层数目能够有效的提高结构的抗冲击性能。X.Liu等[10-11]通过数值模拟计算，分析了梯度金属泡沫夹芯圆柱壳和夹芯板在爆炸荷载下的动力响应，并与传统非梯度夹芯结构做了对比研究，他们发现在相同的荷载作用下梯度结构的抗爆炸性能明显优于非梯度结构，并且密度递减排列的梯度结构具有最佳的抗爆性能。Y.Li等[12]和N.A.Apetre等[13]对梯度夹芯结构的动力响应做了数值模拟计算，指出通过合理的芯层设计能够提高芯层的吸能效率，减轻或者完全避免结构的破坏。L.Zhang等[14]对梯度波纹夹芯板在冲击波作用下的动力响应进行了实验研究，他们发现平缓递减梯度芯层排列方式能够有效的减小后面板的挠度，提高结构的抗冲击性能。

本文中利用弹道冲击摆锤系统对分层梯度蜂窝夹芯板在爆炸荷载作用下的动力响应做了实验研究，分析梯度蜂窝夹芯板在爆炸荷载作用下的变形失效模式，并与传统非梯度蜂窝夹芯板的抗爆性能做了对比，并通过一维应力波理论，分析应力波在梯度芯层中的传播规律。

1 试样制备与实验方法

试件采用5052铝合金作为前后面板，厚度为1.0 mm，蜂窝芯层基体材料采用1200铝合金，蜂窝孔为标准的正六边形，壁厚τ=0.04 mm，芯层根据蜂窝孔边长a的不同分为取3种：a=1.5 mm (记为S)，a=2.0 mm(记为M)，a=2.5 mm(记为L)。各芯层的初始相对密度为：

(1)

式中：τ为蜂窝孔壁厚，θ=30°，ρs为基体材料密度。

通过不同的排列方式，可得到6种密度梯度组合形式不同的试件，记为G1～G6，每种试件芯层从前到后的排列方式见表1；另外，制作3种芯层密度不同的非梯度试件，记为UG1、UG2和UG3，芯层分别为S、M、L型。试件的长和宽均为300 mm，厚度为23.0 mm。为了防止蜂窝板在变形过程中层间相互嵌入破坏，采用厚度为0.1 mm的5052铝合金箔膜作为隔层，芯层与面板及隔层板之间不粘结，为自由接触。夹芯板试件如图1所示，其中C1、C2、C3分别代表从前到后的3个芯层。各种材料力学性能参数如表2所示，其中σ0.2为材料的屈服强度，E为弹性模量，ν为泊松比。

图1 夹芯板试件Fig.1 Specimens of sandwich plates

表1 试件分组与实验结果Table 1 Specimen configurations and blast loading results

表2 材料力学性能参数Table 2 Mechanical properties of aluminum alloys

图2为冲击摆锤系统，冲击摆锤由一根长约1.5 m的工字钢、配重以及夹具3部分组成，通过4根钢绞线悬挂在空中，试件夹具通过螺栓固定在摆锤的夹具中。TNT炸药在试件正前方起爆，起爆后整个系统在爆炸冲击作用下作前后摆动，通过测定摆锤系统摆动的位移来确定试件所受的冲量，摆锤总质量M=151.3 kg，摆动半径为2.69 m，摆动最大角度约为2°，摆动周期T=3.14 s，摆动最大位移为100 mm。TNT装药采用圆柱形(直径为D，长度为L)，如图3所示。

图2 冲击摆锤系统Fig.2 Ballistic pendulum system

图3 TNT装药Fig.3 TNT charge

2 实验结果与分析

2.1 实验结果

表1给出了17个试件的芯层组合方式及实验结果，其中：W为TNT药量，R为TNT炸药距试件前面板的距离，I为爆炸作用的冲量，γ为前后面板的挠度，δ为芯层的压缩量。冲量I是通过冲击摆锤系统测量得到，前后面板的挠度γ及芯层的压缩量δ均为试件最终的残余值。冲量的表达式为[15]：

I=Mx1(2π/T)eβT/4

(2)

β=ln(x1/x2)/T

(3)

式中：x1和x2分别为t=T/4和t=3T/4时摆锤的位移值。

2.2 失效模式分析

夹芯结构的失效模式可以分为3类[16]：模式Ⅰ，非弹性整体大变形；模式Ⅱ，非弹性局部变形；模式Ⅲ，局部贯穿撕裂。对于非梯度试件，3种失效模式如图4所示。图5和图6分别给出了W=20 g，R=200 mm时，梯度试件的剖面图和G1试件的芯层压缩区域分化情况。

图4 非梯度试件变形模式Fig.4 Deformation modes of ungraded sandwich plates

图5 梯度试件变形模式Fig.5 Deformation modes of graded sandwich plates

图6 芯层压缩区域划分Fig.6 Failure pattern of the honeycomb core

从表1以及图4～6可以看出，大多试件的前面板均呈现出局部变形或贯穿失效模式。分别从变形失效模式和塑性区域大小两方面分析夹芯板局部变形特征。从图4中可以看出随着装药距离的减小，非梯度夹芯结构的变形模式从非弹性大变形逐渐向局部压缩贯穿失效模式发展，当装药距离较近时，前面板呈花瓣状撕裂破坏，芯层贯穿，后面板也出现了较明显的局部变形(图4(c))，在后面板上出现局部的鼻状凸起。从图5可以看出，芯层的排列组合对结构的变形有明显的影响，当第1层蜂窝芯层的胞孔较小时，前面板主要为局部变形，前两层芯层变形以胞孔屈曲密实化和整体弯曲变形为主，第3层蜂窝芯层与后面板为整体弯曲大变形；当第1层蜂窝板的胞孔较大时，前后面板呈现整体大变形，芯层主要呈现部分压缩与整体弯曲变形。并且从图5中还可以发现，G6试件各部分变形相对一致，表现出了较好的整体变形性能。

通过对芯层的压缩程度将芯层从中心到边界划分了3个不同的区域，分别为：完全密实化区域；部分密实化区域，非压缩区域[17]。在非梯度试件中也出现了相似的芯层变形区域，如图4所示，如果药量较大，或炸药距离较前面板较近时，在芯层中部完全密实化区域会发生部分贯穿或完全贯穿失效。从图6可以看出，在W=20 g，R=200 mm时，G1试件(芯层相对密度从前到后递增排列)第1层蜂窝中心区域完全被压实，从中心向外芯层的压缩量逐渐减小，第3层蜂窝芯层中部没有出现压实区域，较小的中心区域被部分地压缩，主要呈现出整体的弯曲变形。芯层主要压缩区域和密实化区域的近似呈圆形，表3给出了梯度试件在W=20 g，R=200 mm时变形区域面积对比，其中S1和S2分别为部分压实区域和完全压实区域面积。从表3中可以看出，相对密度较小的芯层变形区域和密实化区域较大，当相对密度最小的芯层(L型)在结构的不同位置时，压缩区域面积在78.5～254.3 mm2范围内。不同试件中芯层相对密度最大的蜂窝芯层(S型)压缩区域相差较大，当它作为C3芯层时压缩区域面积在7.1～12.7 mm2范围内，但当其作为C1芯层时压缩区域面积在103.8～132.7 mm2范围内。由于增加相对密度加大的芯层的压缩量和压缩区域，可以有效提高结构的吸能效果，所以相对密度较大的芯层靠近冲击端时有利于结构抗爆性能的提高。

2.3 抗爆性能对比

结构抗爆性能可以通过对前后面板挠度、变形失效模式、芯层压缩量和能量吸收比例以及应力波传播特性等方面来分析。本小节综合夹芯板后面板挠度和芯层压缩情况与失效模式对比分析不同夹芯板的抗爆性能。图7和图8分别给出了梯度试件与非梯度试件后面板挠度残余值及芯层压缩量的对比。从图7可以看出，在3种工况下，梯度试件后面板挠度明显小于非梯度试件。当W=20g，R=200mm时，与非梯度试件UG3(a=2.5 mm)相比，梯度试件G1、G2、G3、G4、G5和G6后面板挠度分别减小了7.9%、27.5%、20.4%、22.5%、26.7%和 24.6%。在6种梯度组合中，“L-S-M”、“S-L-M”、“S-M-L”3中排列方式表现出了较好的抗爆炸冲击性能。从图8可以看出，不同的芯层排列对各芯层的压缩量有明显的影响，但孔径较大的芯层(L型)在结构的任何位置都有5～6 mm的压缩量，当孔径较小的芯层(S型)靠近前面板时，其压缩量明显增加，由于能量耗散主要通过芯层逐渐压缩屈曲，所以增加芯层的压缩量是提高夹芯结构能量吸收的有效途径。在G6试件中3层芯层都有较大的压缩量，且总压缩量是所有试件中最大的，其呈现出整体变形模式，所以综合以上变形模式，后面板挠度以及芯层的压缩情况，G6试件(芯层相对密度从前到后逐渐减小)的抗爆炸冲击性能较好。

表3 梯度试件芯层变形区域面积对比(W=20 g，R=200 mm)Table 3 Comparison of the deformation area of the core layers (W=20 g，R=200 mm)

图7 不同工况下后面板残余挠度Fig.7 Permanent mid-point deflections of the back-face-sheet under differente conditions

图8 W=20 g，R=200 mm时的芯层压缩量Fig.8 Core compressions while W=20 g, R=200 mm

2.4 应力波在梯度芯层中的传播

梯度芯层中应力波的传播特性对芯层的能量吸收和耗散有明显的影响，而影响应力波传播特性的主要因素为芯层材料的波阻抗[18-20]。由于应力波传播特性较为复杂，目前有关这方面的研究都集中在一维应力波理论的基础上，应力波从一种介质传播到另一种介质中时，在界面上发生反射与透射2种现象，在考虑夹芯结构作为防护层时，主要关注应力波的透射，减小透射波是提高夹芯结构防护性能的主要方法，在一维状态下，透射应力波为：

(4)

式中：ΔσT是透射应力波的应力增量，ΔσI是入射应力波的应力增量，K为透射系数，(ρ0c0)1和(ρ0c0)2分别是前方材料和后方材料的波阻抗，c0为材料的声速。如果(ρ0c0)1<(ρ0c0)2，则ΔσT>ΔσI如果(ρ0c0)1>(ρ0c0)2，则ΔσT<ΔσI。如果有多层材料，则透射系数可表示为[21]:

(5)

式中：(ρ0c0)i为第i层材料的波阻抗。本文中，(ρ0c0)L>(ρ0c0)M>(ρ0c0)S, 所以芯层的组合次序对应力的传播有明显的影响，当芯层相对密度从前到后逐渐减小时，透射系数Kn是最大的，相反，当芯层相对密度从前到后逐渐增加时，透射系数Kn最小。所以G6试件芯层的组合方式有利于减小透射应力波。然而，式(4)和(5)均未考虑材料在变形过程中波阻抗的变化，而随着蜂窝的压缩，蜂窝板的相对密度ρ在不断的改变，满足：

ρ=ρ0/(1-εt)

(6)

式中：ρ0为蜂窝板的初始相对密度，εt为不同时刻蜂窝芯层的压缩应变。

同时蜂窝芯层中的波速也不断变化，而且波速的变化更为复杂，在冲击波在不断压缩的芯层中传播时，冲击波的传播特性变得更加复杂，式(5)已可能不再适用，有关这方面的结论还需要进一步做数值模拟和理论研究。

3 结 论

通过爆炸荷载作用下蜂窝夹芯板的实验，对比分析了梯度夹芯板与传统非梯度结构变形失效模式以及抗爆性能，实验结果表明：(1)夹芯板变形模态主要呈现两种变形模式：非弹性整体大变形和局部贯穿失效；梯度夹芯板芯层的组合方式对结构的变形模式有明显的影响，当第1层蜂窝板的胞孔较小时，前面板呈现局部变形；当第1层蜂窝板的胞孔较大时，前面板呈现整体大变形。(2)分层梯度试件的抗爆性能明显优于非梯度试件，在本文研究的荷载范围内，相对密度从前到后逐渐减小的芯层组合具有最佳的抗爆炸冲击性能。(3)在一维应力波理论下，不考虑波阻抗在芯层压缩过程中的变化，相对密度从前到后逐渐减小的芯层组合有利于应力波的耗散和结构抗冲击性能的提高。但是当冲击波在压缩后的芯层中不断的发生反射和透射时，应当考虑波阻抗变化的影响。

[1] Gibson L J, Ashby M F. Cellular solids:structure and properties[M]. 2nd ed. UK: Cambridge University Press, 1997.

[2] Xu S, Beynon J H, Ruan D, et al. Strength enhancement of aluminium honeycombs caused by entrapped air under dynamic out-of-plane compression[J]. International Journal of Impact Engineering, 2012,47(4):1-13.

[3] Zhao H, Gary G. Crushing behavior of aluminum honeycombs under impact loading[J]. International Journal of Impact Engineering, 1998,21(10):827-836.

[4] Zhang X, Zhang H, Wen Z. Experimental and numerical studies on the crush resistance of aluminum honeycombs with various cell configurations[J]. International Journal of Impact Engineering, 2014,66:48-59.

[5] Liang C-C, Yang M-F, Wu P-W. Optimum desing of metallic corrugated core sandwich panals subjected to blast loads[J]. Ocean Engineering, 2001,28(7):825-861.

[6] McShane G J, Radford D D, Deshhpand V S, et al. The response of clamped sandwich plates subjected to shock loading[J]. European Journal of Mechanics：A: Solids, 2006,25:215-129.

[7] Etemadi E, Khatibi A A, Takaffoli M. 3D finite element simulation of sandwich panels with a functionally graded core subjected to low velocity impact[J].Composite Structures, 2009,89(1):28-34.

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

[9] Gardner N, Wang E, Shukla A. Performance of functionally graded sandwich composite beams under shock wave loading[J]. Composite Structures, 2012,94(5):1755-1770.

[10] Liu X, Tian X, Lu T J, et al. Blast resistance of sandwich-walled hollow cylinders with graded metallic foam cores[J]. Composite Structures, 2012,94(8):2485-2493.

[11] Liu X, Tian X, Lu T, et al. Sandwich plates with functionally graded metallic foam cores subjected to air blast loading[J].International Journal of Mechanical Sciences,2014,84:61-72.

[12] Li Y, Ramesh K T, Chin E S C. Dynamic characterization of layered and graded structures under impulsive loading[J].International Journal of Solids and Structures, 2001,38(34/35):6045-6061.

[13] Apetre N A, Sankar B V, Ambur D R. Low-velocity impact response of sandwich beams with functionally graded core[J]. International Journal of Solids and Structures, 2006,43(9):2479-2496.

[14] Zhang L, Hebert R, Wright J T, et al. Dynamic response of corrugated sandwich steel plates with graded cores[J]. International Journal of Impact Engineering, 2014,65:185-194.

[15] 敬霖,王志华,赵隆茂.爆炸荷载作用下结构冲量的测量[J].实验力学,2009,24(2):151-156. Jing Lin, Wang Zhihua, Zhao Longmao. Measurement of impulse acted on a structure subjected to blast loading[J]. Journal of Experimental Mechanics, 2009,24(2):151-156.

[16] Nurick G N, Langdon G S, Chi Y, et al. Behaviour of sandwich panels subjected to intense air blast：Part 1: Experiments[J]. Composite Structures, 2009,91:433-441.

[17] Zhu F, Zhao L, Lu G, et al. Deformation and failure of blast-loaded metallic sandwich panels- Experimental investigations[J]. International Journal of Impact Engineering, 2008,35(8):937-951.

[18] Makris A, Frost D L, Nerenberg J, et al. Attenuation of a blast wave with a cellular material[C]∥Proceedings of the 20th International Symposium on Shock Waves (ISSW/20). Pasadena, CA, USA, 1996,2:1387-1392.

[19] Bruck H A. A one-dimensional model for designing functionally graded materials to manage stress waves[J]. International Journal of Solids and Structures, 2000,37(44):6383-6395.

[20] Samadhiya R, Mukherjee A, Schmauder S. Characterization of discretely graded materials using acoustic wave propagation[J]. Computation Materials Science, 2006,37(1/2):20-28.

[21] 宋博,胡时胜,王礼立.分层材料的不同排列次序对透射冲击波强度的影响[J].兵工学报,2000,21(3):272-274.Song Buo, Hu Shisheng, Wang Lili. Influnence on the transmitted intensity of shock wave through different tactic orders of layered materials[J]. Acta Armamentarii, 2000,21(3):272-274.

(责任编辑 王小飞)

Dynamic response of functionally graded honeycomb sandwich plates under blast loading

Li Shiqiang, Li Xin, Wu Guiying, Wang Zhihua, Zhao Longmao

(InstituteofAppliedMechanicsandBiomedicalEngineering,TaiyuanUniversityofTechnology,Taiyuan030024,Shanxi,China)

In this paper we report on the tests that investigate the blast resistance of graded sandwich plates. The deformation model, the back-face-sheet deflections and the core compressions have been compared with the test results obtained from tests done on structures with ungraded core layers. The stress transfer characteristics are analyzed based on the one dimensional stress wave theory, indicating that the stress wave transferred factor is smaller in the graded core layers and it is smallest in the relative density-tapered core arrangement specimen. By considering the deformation model, back-face-sheet deflections, core compressions and stress transfer characteristics, the blast resistance of the graded sandwich plates is found to be better than that of the ungraded ones, and in the present loading conditions, the relative density-tapered core arrangement from the front sheet to the back sheet is found to have the best blast resistance.

mechanics of explosion; dynamic response; blast loading; honeycomb sandwich panels; functionally graded materials

10.11883/1001-1455(2016)03-0333-07

2014-10-13；

2015-02-10

国家自然科学基金项目(11172196)；山西省自然科学基金项目(2014011009-1)

李世强(1986- )，男，博士；

吴桂英，wgy2005112@163.com。

O381国标学科代码：13035

A