• ISSN 1001-1455  CN 51-1148/O3
  • EI、Scopus、CA、JST收录
  • 力学类中文核心期刊
  • 中国科技核心期刊、CSCD统计源期刊

编织Kevlar/Epoxy复合材料层合板在冲击荷载下的动态响应

马小敏 李世强 李鑫 王志华 吴桂英

贾东, 黄西成, 胡文军, 张方举. 基于J-C模型的镁合金MB2动静态拉伸破坏行为[J]. 爆炸与冲击, 2017, 37(6): 1010-1016. doi: 10.11883/1001-1455(2017)06-1010-07
引用本文: 马小敏, 李世强, 李鑫, 王志华, 吴桂英. 编织Kevlar/Epoxy复合材料层合板在冲击荷载下的动态响应[J]. 爆炸与冲击, 2016, 36(2): 170-176. doi: 10.11883/1001-1455(2016)02-0170-07
Jia Dong, Huang Xicheng, Hu Wenjun, Zhang Fangju. Fracture behavior of magnesium alloy MB2 under quasi-static and dynamic tension loading based on Johnson-Cook model[J]. Explosion And Shock Waves, 2017, 37(6): 1010-1016. doi: 10.11883/1001-1455(2017)06-1010-07
Citation: Ma Xiaomin, Li Shiqiang, Li Xin, Wang Zhihua, Wu Guiying. Dynamic response of woven Kevlar/Epoxy composite laminatesunder impact loading[J]. Explosion And Shock Waves, 2016, 36(2): 170-176. doi: 10.11883/1001-1455(2016)02-0170-07

编织Kevlar/Epoxy复合材料层合板在冲击荷载下的动态响应

doi: 10.11883/1001-1455(2016)02-0170-07
基金项目: 

国家自然科学基金项目 11402163

山西省自然科学基金项目 2014011009-1

详细信息
    作者简介:

    马小敏(1988—),男,博士研究生

    通讯作者:

    吴桂英,wgy2005112@163.com

  • 中图分类号: O347.3

Dynamic response of woven Kevlar/Epoxy composite laminatesunder impact loading

  • 摘要: 通过编织Kevlar/Epoxy复合材料层合板的平头弹冲击实验,分析了结构在不同冲量下的变形失效模式以及结构的抗冲击性能。实验表明复合板的变形失效模式主要表现为:(1)弹性变形;(2)复合板表面嵌入失效及整体塑性大变形;(3)背面纤维拉伸断裂及分层失效。基于实验研究,运用LS-DYNA 971有限元程序对铺层数不同的复合板在冲击载荷作用下的动态响应过程进行了数值模拟,模拟结果与实验吻合较好,子弹作用区域边缘处首先发生近似圆形的嵌入失效,而在板背面发生近似正方形的破坏区域;计算中重点分析了铺层数对结构动力响应的影响,在一定冲量范围内,通过对铺层数的优化,能够有效地减小后面板挠度,提高结构的能量吸收效率,增强结构的抗冲击性能。
  • 镁合金是一种潜力巨大的轻质工程材料,有着广阔的应用前景,其中MB2是镁合金的典型代表。材料在复杂应力状态下的失效破坏行为是工程材料和结构设计的理论基础。目前,镁合金在动静态加载下的力学行为研究以单向应力状态为主[1-3],且能实现复杂应力动态加载的实验设备与技术尚不成熟。J.W.Hancock等[4]研究发现金属的延性明显依赖于应力的三轴状态。近年来,O.S.Hopperstad等[5-6]对结构钢进行了拉伸加载实验和数值分析,结果表明材料的等效破坏应变随应力三轴状态的变化趋势与J.W.Hancock等的研究结论一致;Y.B.Bao[7]在单轴加载条件下对铝合金2024进行了不同应力状态的实验研究;D.Anderson等[8]研究了应力三轴度和应变率对DP780钢破坏行为的影响。汤安民等[9]、李智慧等[10]研究了不同金属材料宏观断裂形式及断裂机理与应力状态的关系, 陈刚等[11]建立了应力三轴度和应变率相关的45钢损伤失效模型, 朱浩等[12]分析了应力三轴度和应变率对铝合金6063力学性能的影响, 张伟等[13]以应力三轴度和应变率效应为基础给出了铝合金7A04的本构关系和失效模型。

    所以,对于钢材、铝合金等常用金属材料,其失效破坏行为的研究已经很多,而对镁合金复杂应力状态下的失效破坏行为研究还很少。本文中将通过复杂应力下的动静态拉伸实验及数值模拟,分析应力状态和应变率对镁合金MB2失效破坏的影响,基于Johnson-Cook本构及失效模型建立适用于镁合金MB2的破坏准则,分析其宏观破坏模式及其微观失效机理,为镁合金MB2在实际工程中的应用提供强度设计依据和理论支撑。

    材料所受应力状态不同时,材料内产生的塑性变形与应力集中程度将不同,为了表征材料的复杂应力状态,引用应力三轴度η为应力状态参数,即平均应力与等效应力的比值σH/σ。在延性金属复杂应力状态的失效破坏研究中,依据Bridgman原理[14]设计的缺口试件拉伸实验是最主要的研究方法,根据圆弧缺口颈部应力方程解,可以得到应力三轴度的计算公式:

    η=13+ln(1+a2r22aR)
    (1)

    式中:aR分别为最小横截面的半径及缺口半径,r为到横截面中心的距离,a0为最小横截面初始半径。基于塑性不可压的假设,可以得到缺口处的等效应变定义公式如下:

    ˉε=2ln(a0/a)
    (2)

    a取试件断裂时的横截面直径af时,式(2)计算得到的即为试件的等效破坏应变εf。静态拉伸实验试件采用直径为5 mm的光滑圆柱以及最小横截面直径为6 mm,缺口曲率半径分别为9、6、3 mm的缺口圆柱试件。动态拉伸实验试件采用直径为3 mm的光滑圆柱以及最小横截面直径为2 mm,缺口曲率半径分别为1.0、1.5、2.0、3.0 mm的缺口圆柱试件。静态拉伸加载由材料试验机MTS 810完成,而动态拉伸加载则通过分离式Hopkinson拉杆(SHTB)实验装置完成,如图 1所示。

    图  1  SHTB实验系统
    Figure  1.  SHTB testing system

    对于光滑圆柱试件的动态拉伸实验,按照一维应力波理论进行数据处理分析可知,试件在不同应力幅值的入射波作用下表现出了明显的应变率效应,不同应变率下的真实应力应变曲线如图 2所示。

    图  2  镁合金MB2真实应力应变曲线
    Figure  2.  True stress-strain curves of magnesium alloy MB2

    在高速冲击载荷作用下,加载过程为绝热过程,材料在该过程中的温度变化ΔT可由塑性变形进行求解:

    ΔT=βρcpεp0σdεp
    (3)

    式中:β为摩擦能量转换系数,ρ为材料密度,cp为比定压热容,εp为等效塑性应变。对于镁合金MB2,摩擦能量转换系数通常取为0.925,密度为1 800 kg/m3,比定压热容为1 040~1 148 J/(kg·K)。

    鉴于镁合金MB2的应变率效应,这里引用经典的J-C本构模型对实验结果进行拟合。在室温、准静态(0.001 s-1)拉伸条件下,J-C本构模型的应变率项和温度项均为1,此时材料的本构模型退化为:

    σ=A+Bεnp
    (4)

    采用最小二乘法对参数ABn进行拟合可以得到:A=192 MPa、B=218.3 MPa、n=0.370 56。由于在300 s-1下试件未被拉断,所以这里以800和1 800 s-1应变率下的抗拉强度及其对应的等效塑性应变作为比较点与参考应变率0.001 s-1进行比较分析,获得参数C的算术平均值为0.015。对于反映温度效应的参数m,由于本研究未涉及高温下的动态拉伸实验,所以此处引用文献[15]中的研究结论,给出参数m的值为0.95。

    对于光滑圆柱试件和缺口圆柱试件,由于其加载过程中应力状态的差异,导致其动静态加载条件下的破坏形式也有所不同,如图 3所示。

    图  3  不同试件破坏模式
    Figure  3.  Macro fracture patterns of different specimens

    光滑圆柱试件在整体上发生剪切破坏,且在破坏前有一定的颈缩,试件断口处有一定的杯锥特征。缺口试件在准静态条件下整体上发生正拉断,断口边缘有杯锥特征;在动态加载下整体表现为杯锥形断裂。通过测量断裂后缺口处的最小横截面直径,得到了动态拉伸试件不同加载条件下的等效破坏应变,如表 1所示。

    表  1  不同类型试件的等效破坏应变
    Table  1.  Equivalent fracture strain of different specimens
    入射波幅值/MPa 等效破坏应变
    光滑 R=3.0 mm R=2.0 mm R=1.5 mm R=1.0 mm
    50 未断裂 0.333 8 0.250 2 0.254 1 0.205 3
    120 0.311 7 0.307 1 0.225 7 0.239 7 0.197 7
    350 0.289 3 0.249 2 0.180 0 0.175 8 0.177 6
    下载: 导出CSV 
    | 显示表格

    J-C失效模型类似于J-C本构模型,采用多项乘积的形式对应力状态参数、应变率参数以及温度效应参数进行解耦。由于在本研究中未涉及等效破坏应变的温度效应,所以此处给出简化后的表达形式如下:

    ˉεf=[D1+D2exp(D3η)][1+D4ln(˙ε/˙ε0)]
    (5)

    式中:D1D2D3D4为失效模型系数。研究发现,试件最小横截面处的应力三轴度是随等效应变的增加而不断变化的,整个加载过程的应力三轴度计算需要考虑到应变的累积效应[16],所以D.M.Goto等[17]、Y.Bao等[18]在研究不同延性金属材料时,对整个加载过程的平均应力三轴度ηav进行了定义:

    ηav=1ˉεfˉεf0η(ˉε)dˉε
    (6)

    采用ABAQUS/Explicit分析软件对镁合金MB2不同类型试件的动态拉伸实验过程进行模拟计算,鉴于SHTB和试件的对称性,以轴对称单元进行建模,以拉伸本构关系和实验中采集的入射波为输入条件,模拟得到了试件在加载过程中应力三轴度的变化情况,如图 4所示。

    图  4  不同类型试件应力三轴度变化
    Figure  4.  Stress triaxiality curves of different specimens

    同时,由于缺口试件动态加载过程中,其塑性变形集中于缺口区域,并不满足一维应力波原理,无法通过理论公式对应变率进行求解。所以这里根据缺口试件等效应变的定义,借助数值模拟得到最小横截面直径的变化和应力波作用时程,通过时间积分得到试件的平均应变率:

    ˙εav=1Δttft0˙ε(t)dt=2Δtln(a0af)
    (7)

    式中:Δt为试件发生拉伸变形的时间,t0为载荷作用开始的时刻,tf为载荷作用结束的时刻。通过数值计算得到不同类型试件入射波在试件上的作用时间,进而计算得到其平均应变率,如表 2所示。

    表  2  缺口试件的平均应变率
    Table  2.  Average strain rates of different specimens
    入射波幅值/MPa 平均应变率/s-1
    R=3 mm R=2 mm R=1.5 mm R=1 mm
    50 828 968 1 047 1 239
    120 2 315 2 573 2 569 2 694
    350 3 905 4 129 4 914 4 981
    下载: 导出CSV 
    | 显示表格

    在准静态(˙ε=˙ε0)拉伸条件下,J-C失效模型的应变率项为1,且镁合金MB2的等效破坏应变存在明显的转折点[19],通过拟合得到,应力三轴度小于0.66时,参数D1D2D3的值分别为1.04×10-5、0.073 7以及0.323;当应力三轴度大于0.66时,其值分别为105.7、-0.103以及0.236。对于镁合金MB2,J-C失效模型中的(εf/ε0-1)与ln(˙ε/˙ε0)并不满足线性关系,所以这里对模型中的应变率项进行了修正,修正后的失效模型为:

    ˉεf=[D1+D2exp(D3η)]{1+D5 [ 1-exp(D6ln(˙ε/˙ε0))]}
    (8)

    式中:D5D6为应变率项参数。通过拟合分析发现,修正模型中(εf/ε0-1)与ln(˙ε/˙ε0)的关系与实验数据吻合,D5D6的值分别为8.87×10-20、2.801,如图 5所示。

    图  5  应变率项参数拟合曲线
    Figure  5.  Fitted curve for the item of strain rates

    结合参数拟合的结果,可以得到镁合金MB2最终的破坏准则为以应力三轴度0.66为转折点的分段函数形式,如图 6所示。图中实心圆点为实验结果,曲面为破坏准则的拟合结果。

    图  6  镁合金MB2的破坏准则曲面
    Figure  6.  Fracture criterion surface of magnesium alloy MB2

    通过失效破坏曲面可以发现,镁合金MB2的断裂延性(等效破坏应变)随着应力三轴度的增加先增大后减小,随着应变率的增大而不断减小。从破坏模式上来看,随着应力三轴度的增大,材料先发生剪切破坏,随后又发生正拉断破坏;随着应变率的增大,光滑圆柱及缺口圆柱试件在动静态载荷作用下的破坏形式基本一致。所以,应力状态对断裂延性和破坏模式的影响存在明显的转折点,而应变率的影响则不存在转折点。

    为了对该失效模型进行验证,这里通过ABAQUS/Explicit分析软件和单元失效法对不同类型试件的动态失效破坏形式进行了数值模拟分析,其模拟结果如图 7所示。

    图  7  不同试件失效破坏模式的数值模拟结果
    Figure  7.  Fracture patterns of different specimens from numerical simulation

    从数值模拟结果可以看出,无论是光滑圆柱试件还是缺口圆柱试件,其最终破坏形式均呈现出杯锥形特征,这与实验现象是基本吻合的;但对于光滑圆柱试件,该模型不能体现试件的整体剪切断裂特征,这与J-C失效模型的物理涵义以及镁合金MB2在不同应力状态下的微观损伤机理有关。

    宏观破坏行为与材料的微观损伤变形机制密切相关,为进一步认识应力状态对镁合金MB2破坏行为的影响,对准静态条件下光滑圆柱及缺口试件的断口形貌进行微观扫描观察,如图 8所示。

    图  8  试件微观断口形貌
    Figure  8.  Micro fracture morphology of specimens

    对于光滑拉伸试件,微观断口上能看到河流花样特征即解理断裂,同时也呈现出了一定的微孔洞特征,即有韧性断裂的特点,如图 8(a)所示,属韧脆性混合断裂;对于缺口试件,不同缺口程度均表现出较明显的微孔洞断裂的特征,在图 8(b)中能清晰地看到韧窝聚合、连接的特点。这些特征与试件在宏观破坏形式上的表现是一致的,即光滑圆柱试件破坏是由混合损伤机制引起的,而缺口试件破坏则是由微孔洞损伤演化机制决定的。J-C失效模型在微观上反映的是材料的微孔洞损伤演化机理,在宏观上体现为该微观损伤机理导致的材料断裂延性变化,不能体现其他变形机制的作用,这就是模拟结果可以描述光滑圆柱试件的杯锥形破坏特点而无法描述其整体剪切特征的原因。

    (1) 通过缺口圆柱试件设计能够实现镁合金MB2不同程度的多向应力状态,结合SHTB装置可以实现其在复杂应力状态下的动态拉伸加载研究。

    (2) 在研究的实验条件范围内,应力状态对镁合金MB2破坏行为的影响存在明显的转折点,而应变率(5 000 s-1以内)则不存在转折点。通过模型修正和参数拟合得到了镁合金MB2的等效破坏应变准则,且其数值模拟结果可以较好地体现拉伸试件的杯锥形破坏特征。

    (3) 镁合金MB2不同应力状态下的破坏行为与其微观损伤机理密切相关,随着应力三轴度的增加,其微观损伤机制由混合断裂转变为韧窝断裂,这是导致其宏观延性和破坏模式发生变化的重要原因。

  • 图  1  实验装置示意图

    Figure  1.  Schematic diagram of experimental devices

    图  2  复合板的变形失效模式

    Figure  2.  The deformation and failure modes of the plates

    图  3  复合板背面中心点位移时程曲线

    Figure  3.  Displacement histories ofthe mid-span in the back face

    图  4  Kevlar纤维复合板有限元模型及冲击实验照片

    Figure  4.  Finite element model of the structure and its photo in the experiment

    图  5  Kevlar纤维复合板受撞击实验与数值模拟对比

    Figure  5.  Comparison of the experimental and simulated final deformation modes under impact

    图  6  子弹和Kevlar纤维复合板作用的过程

    Figure  6.  Process of projectile impacting the Kevlar laminates

    图  7  纤维复合板中应力分布时程曲线

    Figure  7.  Histories of stress distribution in the Kevlar laminates

    图  8  纤维复合板背面中心点量纲一残余挠度随量纲一冲量变化的规律

    Figure  8.  Relation between normalized residual deflectionand normalized impulse at the mid-span on the back face

    图  9  纤维复合板吸能效率随量纲一冲量变化的规律

    Figure  9.  Relation between energy absorption efficiencyand normalized impulse

    表  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 Ⅱ Ⅲ
    下载: 导出CSV

    表  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
    下载: 导出CSV

    表  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
    下载: 导出CSV
  • [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/njlgdxxb200705021

    Gu 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
  • 期刊类型引用(4)

    1. 祁武超,刘炳宏,田素梅. 镁合金带解锁装置热冲击断裂特性及拓扑优化. 兵器装备工程学报. 2023(11): 47-56 . 百度学术
    2. 周伦,苏兴亚,敬霖,邓贵德,赵隆茂. 6061-T6铝合金动态拉伸本构关系及失效行为. 爆炸与冲击. 2022(09): 113-124 . 本站查看
    3. 衣海娇,甄莹,曹宇光,张士华,史永晋. 6061-T6铝合金断裂应变与应力三轴度关系研究. 机械强度. 2020(03): 551-558 . 百度学术
    4. Dewang Zhao,Daxin Ren,Kunmin Zhao,Pan Sun,Xinglin Guo,Liming Liu. Ultrasonic Welding of Magnesium–Titanium Dissimilar Metals:A Study on Thermo-mechanical Analyses of Welding Process by Experimentation and Finite Element Method. Chinese Journal of Mechanical Engineering. 2019(06): 191-201 . 必应学术

    其他类型引用(11)

  • 加载中
图(9) / 表(3)
计量
  • 文章访问数:  6031
  • HTML全文浏览量:  2440
  • PDF下载量:  1127
  • 被引次数: 15
出版历程
  • 收稿日期:  2014-08-13
  • 修回日期:  2014-11-10
  • 刊出日期:  2016-03-25

目录

/

返回文章
返回