Dynamic response of woven Kevlar/Epoxy composite laminatesunder impact loading
-
摘要: 通过编织Kevlar/Epoxy复合材料层合板的平头弹冲击实验,分析了结构在不同冲量下的变形失效模式以及结构的抗冲击性能。实验表明复合板的变形失效模式主要表现为:(1)弹性变形;(2)复合板表面嵌入失效及整体塑性大变形;(3)背面纤维拉伸断裂及分层失效。基于实验研究,运用LS-DYNA 971有限元程序对铺层数不同的复合板在冲击载荷作用下的动态响应过程进行了数值模拟,模拟结果与实验吻合较好,子弹作用区域边缘处首先发生近似圆形的嵌入失效,而在板背面发生近似正方形的破坏区域;计算中重点分析了铺层数对结构动力响应的影响,在一定冲量范围内,通过对铺层数的优化,能够有效地减小后面板挠度,提高结构的能量吸收效率,增强结构的抗冲击性能。
-
关键词:
- 固体力学 /
- 动态响应 /
- 冲击荷载 /
- Kevlar/Epoxy复合材料层合板
Abstract: Based on the blunt projectile impact test of woven Kevlar/Epoxy composite laminates, the deformation and failure modes of the composite laminates subjected to impact load were analyzed. Experimental results show that the deformation and failure behaviors were exhibited in the following ways: the global elastic deformation, the global plastic deformation with local embedded failure on the front surface, and the delaminated failure with fibers tension fracture on the back surface. The finite element software LS-DYNA 971 was employed to analyze the dynamic response of the woven Kevlar/Epoxy composite laminates subject to impact loading. Numerical simulation results show that there is a good agreement of the deformation/failure modes and the back face center-point deflection of the specimens, with those of the experimental results. The failure area on the front face is a circle embedded region, but a square failure region on the back face. The numerical simulation is focused on studying the effects of the number of layers on the dynamic response of the structure. Optimizing the number of the layers can effectively reduce the permanent deflection, increase the energy absorption efficiency and improve the impact resistance performance of the structure within a given range of impulses.-
Key words:
- solid mechanics /
- dynamic response /
- impact load /
- Kevlar/Epoxy composite laminates
-
金属及合金在强动载荷作用下将发生近表面破碎,产生大量高速运动碎片和粒子,对其他结构或器件造成损伤破坏。因此,金属及合金材料近表面的动态破碎是爆炸与冲击动力学领域高度关注的关键科学问题,对该问题的研究在装甲防护、航空航天、公共安全等方面具有重要应用价值[1-5]。Walsh等[6]和Asay等[7-9]早就开始了材料近表面的动态破碎研究。随着诊断技术的不断发展,研究人员通过实验研究发现,材料近表面动态破碎存在3种主要机制,即层裂、微层裂和微喷射,这3种机制强烈依赖于冲击压力、温度、相结构、材料微结构和加载波形等因素[10-12]。目前,相关研究基本都以简单一维应变加载为主,关于材料在复杂加载状态下动态破碎的研究还相对较少。因此,对复杂应变/应力条件下材料近表面动态破碎的认识严重不足。
马赫透镜作为一种高压加载的特殊方法,通过冲击波的马赫反射和对碰可以在材料内部形成马赫杆加载,且在马赫透镜的中心位置产生三角波,它已经被用于高压状态方程[13-14]和疏松材料压实[15]方面的研究。该三角波与平板撞击、爆轰加载和激光冲击下产生的一维应变三角波存在显著不同。前者的压力幅值在一定的时间内随传播距离的变化几乎不变,而后者的压力幅值则会因为追赶稀疏波的影响随时间的变化急剧衰减。因此,马赫杆加载虽然是一维应变加载,但其应力状态更复杂,为研究复杂加载条件下材料的动态破碎提供了新的机遇。
本文中,将数值模拟和实验相结合,利用大口径火炮在无氧铜样品中产生马赫杆加载,并利用激光粒子速度干涉仪测量样品自由面的粒子速度剖面,研究在马赫杆加载下无氧铜样品的近表面动态破碎行为,揭示在不同加载条件下无氧铜样品自由面的2种动态破碎机制。
1. 马赫杆加载模拟设计
马赫透镜构型如图1所示,内外圆柱由不同波阻抗的金属组成。飞片高速撞击后,外层圆柱材料中的冲击波速度比内圆柱中的高,使外侧冲击波能透射入内层圆柱中,形成锥形收敛,会聚到中心轴线上,并与内圆柱中的冲击波叠加产生再加载,使内圆柱的局部压强和粒子速度显著提升。边界处的外侧冲击波入射可看作是斜冲击波的反射状态,通过斜冲击波的入射可以在内圆柱材料中心线上产生稳定的马赫杆[14]。虽然马赫杆加载宏观上表现为一维应变状态,但整个内圆柱中的冲击波阵面随时间和空间剧烈变化,因此其实质的加载状态更复杂。
在内圆柱中,除了波阵面处的斜激波入射与反射,在界面处还有内外冲击波的相互作用,结果是向内侧圆柱中反射稀疏波。在入射波向中心汇聚时,反射的稀疏波也到达中心,跟随在马赫杆后。从而在马赫杆后形成低压卸载区,使冲击前沿压强到达峰值后又迅速下降,马赫杆波形就由方波转变成三角波。
基于马赫透镜的原理分析,对马赫杆加载进行了模拟设计。利用有限元分析软件AUTODYN,采用Lagrange方法和二维旋转轴对称结构建立模型。材料动力学模型选用最常见的Mie-Grüneisen状态方程[16]和Johnson-Cook(JC)本构模型[17],其模型参数分别见表1[18-20]和表2[17, 21-23]。模拟设计的两类马赫杆实验(Mach-1和Mach-2)参数如表3,分别对应在无氧铜中产生最高压力99.4和24.2 GPa。两类马赫杆实验的压力流场随时间的演化如图2所示。
表 3 设计参数Table 3. Parameters of experimental design实验编号 材料 飞片速度/(km·s−1) 厚度/mm 直径/mm 飞片 外圆 内圆 飞片 样靶 内圆 外圆 Mach-1 304不锈钢 LY12 铝 无氧铜 1.40 3.0 16.0 4.8 38.0 Mach-2 304不锈钢 TC4钛合金 无氧铜 0.50 12.0 26.0 14.0 45.0 图3为实验Mach-1中内圆柱中心线上等间距拉格朗日点的压力和粒子速度剖面模拟结果。0.7 μs前,冲击波为方波。0.7 μs后,峰值粒子速度开始提升,且波形开始变化。其变化来自两部分,分别是斜冲击的入射波和界面反射的稀疏波,导致在提升峰值粒子速度的同时,后续的卸载稀疏波也变得更陡峭,使冲击平台逐步缩短。2.2 μs时,粒子速度幅值基本保持不变,且能在较长时间内保持这种三角波的形式。
2. 马赫杆加载实验技术
内外圆柱之间采用间隙配合,用环氧树脂粘结。为确保冲击实验时飞片与靶装置正面碰撞,将样品安装在有机玻璃支架上,经校准后用样品环将样品固定在炮口上,如图4(a)所示。火炮驱动304钢飞片以设计速度撞击马赫透镜装置。采用多普勒光探针系统(Doppler pin system, DPS)测量内外圆柱自由面的粒子速度,DPS探头分布见图4(b)~(c)。在实验Mach-1中,在内圆柱中心处设置一根DPS探针,用于测量马赫杆中心的运动历程。在铝套筒后表面距圆心3.0、5.0、7.0、9.0 mm处,布置系列DPS探针,用以判断碰撞的倾斜角度以及比较模拟结果与实验结果,并确定模拟需要的黏性等参数。在Mach-2实验中,碰撞速度较低,马赫杆的范围更大,因此在内圆柱的半径0、1.0、1.4、1.8、2.6和3.4 mm处布置了DPS探针,对应图4(c)中的点E、F、G、H、I和J,用于分析马赫杆的径向行为。飞片的碰撞速度也通过DPS测速确定。
3. 实验结果与讨论
DPS测量获得的飞片碰撞速度分别为1325和605 m/s。在实验Mach-1中,铝套筒后表面不同位置处的粒子速度剖面如图5所示。由同心圆上4个DPS信号的起跳时间间隔可知,碰撞倾斜角度为0.376°,因此该实验中飞片与靶装置的碰撞可以近似为完全正碰。在实测弹速下,铝套筒后自由面粒子速度的数值模拟结果和实验结果具有较好的一致性,由此证实实验中实现了预期的马赫杆加载。马赫杆稳定加载后的峰值压力分别达到95.75和32.38 GPa。结合模拟进行进一步分析发现,在自由面半径为3.0 mm的位置处,由于靠近内外圆柱的界面,受到内边界反射稀疏波的影响,冲击上升沿的峰值强度有所降低;在自由面半径为5.0 mm的位置处,冲击波基本没有受到来自外侧和内测边界处卸载稀疏波的影响,是一个完整的矩形冲击波;而在自由面半径为7.0、9.0 mm的位置处,冲击波则是受到了外侧稀疏波的影响,粒子速度剖面呈现三角波的特征,冲击波在到达自由面后开始卸载,粒子速度明显下降。图5(b)~(d)的粒子速度剖面在1.0 μs左右出现回弹上升,应该是发生了层裂,而且越接近外侧,层裂出现越早。
实验Mach-1中内圆柱中心处的DPS测量结果见图6,可以看出无氧铜后自由面可能发生了微层裂,导致近表面物质碎化为细小碎片云并以接近自由面的速度脱离基体。这是由于稳定后的马赫杆中心处的压力由碰撞初始的28.92 GPa提升到了95.75 GPa,使无氧铜在卸载后进入近熔化区。
在实验Mach-2中,由于马赫杆的范围更大,因此可以研究马赫杆不同位置的动力学过程。模拟设计中,点E、F、G和H位于马赫杆平台范围内,而点I和J位于过渡区。如图7所示,粒子速度剖面都有明显的三角波层裂信号。为描述该层裂行为,数值模拟中采用最大拉应力断裂模型(设最大拉应力为2.6 GPa),可以获得实验Mach-2中不同位置处的粒子速度剖面。模拟结果与实验结果进行比较,可见两者的一致性较好,数值模拟较好复现了不同位置的层裂信号。该层裂行为是由于马赫杆产生的三角波加载到达自由面后经卸载产生拉伸应力所致。如图8所示,马赫透镜中内圆柱冲击波形成的层裂除了轴向应力的影响外还与其在径向分布的差异有关,不同位置的层裂规模和距离差异较大。内圆柱中的马赫杆冲击波阵面沿径向分布不均,冲击波阵面为锥状,点E~J的峰值速度和冲击波到达时间分别随半径增大而减小和提前(见图7)。这与平面加载的层裂行为存在显著的不同,平面加载的层裂在其径向上并无变化。中心线处的马赫杆虽然到达自由面的时间最迟,但强度最高,形成的三角波层裂位置也最靠前,层裂片速度最快,而两侧的冲击波强度相对较低,层裂位置也靠后。最终的层裂片沿径向各个位置的速度和厚度也不同,整体呈凸形(见图8)。
4. 结 论
以无氧铜为研究对象,模拟设计了2类不同冲击波压力的马赫杆加载实验。模拟结果与实验结果十分接近,马赫杆相关的主要特征也都比较符合,表明在实验上用马赫透镜实现了马赫杆加载。同时,观察到了无氧铜在2种不同马赫杆加载下产生的2种不同近表面破碎模式。在高压下,由于无氧铜卸载进入近熔化区,发生微层裂破坏。在低压下,无氧铜处于固态,三角波发生自由卸载形成层裂。与平面冲击产生的层裂不同,由于马赫杆冲击波状态径向分布不均匀,产生的层裂片呈凸形。
该工作对进一步了解金属材料的动态破碎行为具有一定价值,并为研究复杂加载状态下的材料破坏提供了新的实验方法。
-
表 1 实验数据
Table 1. Experimental data
No. n h/mm v/(m·s-1) I/(N·s) W/mm 变形失效模式 1 18 5.11 13.30 16.49 1.8 Ⅰ 2 18 4.98 24.70 30.63 8.3 Ⅱ 3 18 5.15 32.60 40.42 12.4 Ⅱ Ⅲ 4 18 5.13 36.47 45.22 13.5 Ⅱ Ⅲ 5 18 5.10 36.50 45.26 13.3 Ⅱ Ⅲ 表 2 Kevlar纤维平纹织物的材料参数
Table 2. Material properties of the Kevlar composite fabric
ρ/(kg·m-3) E/GPa Gab/GPa Gca/GPa ν Xt/GPa Xc/GPa Yt/GPa Yc/GPa Sc/GPa 1 400 59.5 5.18 5.18 0.34 1.20 0.23 1.20 0.23 0.12 表 3 Kevlar纤维复合板在不同速度冲击下实验与数值模拟对比
Table 3. Comparison of the experimental and simulated results at different impact velocities
v/(m·s-1) dmax/mm W/mm 实验 数值模拟 ε/% 实验 数值模拟 ε/% 13.30 10.3 9.3 -9.71 1.8 1.5 -16.67 24.70 14.2 14.8 4.23 8.7 7.9 -9.20 32.60 18.3 17.6 -3.83 12.4 11.2 -9.68 36.47 20.7 20.3 -1.93 13.5 13.6 0.71 36.50 19.9 20.3 2.01 13.3 13.6 2.21 -
[1] 王震鸣.复合材料力学和复合材料结构力学[M]. 2版.北京:机械工业出版社, 1991. [2] Cunniff P M. An analysis of the system effects in woven fabrics under ballistic impact[J]. Textile Research Journal, 1992, 62(9):495-509. doi: 10.1177/004051759206200902 [3] Cunniff P M. A semi-empirical model for the ballistic performance of textile-based personnel armor[J]. Textile Research Journal, 1996, 66(1):45-58. doi: 10.1177/004051759606600107 [4] 顾冰芳, 龚裂航, 徐国跃.Kevlar纤维叠层织物防弹机理和性能研究[J].南京理工大学学报:自然科学版, 2007, 31(5):638-642. http://d.old.wanfangdata.com.cn/Periodical/njlgdxxb200705021Gu Bingfang, Gong Liehang, Xu Guoyue. Bullet proof mechanism and performance of Kevlar multi-layers textile[J]. Journal of Nanjing University of Science and Technology: Natural Science Edition, 2007, 31(5):638-642. http://d.old.wanfangdata.com.cn/Periodical/njlgdxxb200705021 [5] Barauskas R, Abraitiene A. Computational analysis of impact of a bullet against the multilayer fabrics in LY-DYNA[J]. International Journal of Impact Engineer, 2007, 34(7):1286-1305. doi: 10.1016/j.ijimpeng.2006.06.002 [6] Kostopoulos V, Markopouls Y P, Giannopoulos G, et al. Finite element analysis of impact damage response of composite motorcycle safety helmets[J]. Composites Part B: Engineering, 2002, 33(2):99-107. doi: 10.1016/S1359-8368(01)00066-X [7] Taraghi I, Fereidoon A, Taheri-Behrooz F. Low-velocity impact response of woven Kevlar/Epoxy laminated composites reinforced with multi-walled carbon nanotubes at ambient and low temperature[J]. Materials and Design, 2014, 53(01):152-158. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=c73672e3c6e12f9da0a06afb6e7b5ab2 [8] Reis P N B, Ferreira J A M, Zhang Z Y, et al. Impact response of Kevlar composites with nanoclay enhanced epoxy matrix[J]. Composites Part B: Engineering, 2013, 46(3):7-14. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=65953763f38e4cd0ac1f3271593ec476 -