Mechanical properties of the mixed cellular material with soft matrix and its response to repeated impacts
-
摘要: 为了探索具有优异吸能性能的软基体混合胞孔材料的力学性能,研究该类材料在多次冲击下的冲击响应和材料的可恢复性,对一种软基体混合胞孔材料—人工软骨仿生超材料(artificial cartilage foam,ACF)进行了不同速度下的单轴拉伸和压缩实验,得到了ACF材料在不同应变率条件下的应力-应变曲线。并利用落锤冲击实验机对ACF材料和另一种软基体混合胞孔材料—发泡聚丙烯材料(expanded polypropylene,EPP)进行了多次冲击下的对比测试,得到了2种材料在单次和多次冲击下的动力学响应。实验结果表明:ACF材料是一种应变率敏感的材料,随着应变率的提升,材料的弹性模量、抗拉强度和抗压强度均逐渐提高;在50 J 冲击能量作用下,ACF材料能够吸收96%以上的冲击能量,远高于EPP材料的70%,ACF材料具有更加优异的吸能性能;5次冲击后ACF材料的最大峰值力、最大变形量和吸能能力几乎不变。相比于EPP材料,ACF材料有良好的可恢复性,且具有稳定的多次抗冲击能力。这些研究为软基体混合胞孔材料在多次冲击防护中的应用提供了实验依据。Abstract: The mixed cellular materials with soft matrix are a new type of cushioning and protective materials which have excellent energy absorption properties. In order to study the effect of strain rate on the mechanical behaviors of this kind of materials, uniaxial tensile and compression experiments were conducted on the artificial cartilage foam (ACF) material, at different velocities to obtain the stress-strain curves of the ACF material under different strain rate conditions. Based on the obtained stress-strain curves, the elastic moduli and material strengths of the ACF material were gained under different strain rate conditions. And the comparative tests of the ACF material and expanded polypropylene (EPP) material of the same size and thickness under multiple impacts were carried out by a drop-hammer impact test machine. By comparing the impact responses of the two materials under single and multiple impact loads, the energy absorption characteristics, and the stability of the energy absorption characteristics of the two materials were analyzed. The results reveal that the ACF material is a strain rate-sensitive material, and the stress-strain curves under different strain-rate conditions take on the same trend. The elastic modulus, tensile and compressive strengths of the material gradually increase with the increasing strain rate. Under the action of 50-J impact energy, the ACF material can absorb more than 96% of the impact energy, higher than the 70% of the impact energy absorbed by the EPP material. Moreover, the maximum displacement of the ACF material is only 40% of that of the EPP material. Therefore, the ACF material has more excellent energy absorption performance than the EPP material. The peak force, maximum displacement, and energy absorption ability of the ACF material were almost unchanged after five impacts. Compared with the EPP materials, the ACF material has more favorable recoverability and more stable repeated impact resistance. The research of the work can provide an experimental basis for the application of the mixed cellular materials with soft matrix in multiple impact protection.
-
高超声速飞行器热防护系统(thermal protection system, TPS)的结构材料必须具有轻质、耐高温、环境稳定等特点, 以满足长航程安全服役的要求。近年来, 陶瓷基复合材料在飞行器大面积热防护领域应用广泛, 这类材料可以比传统的金属TPS减重50%, 同时可以延长使用寿命, 降低制造成本[1]。
基于化学气相沉积技术(chemical vapor infiltration, CVI)[2]制造的碳纤维增韧碳化硅陶瓷基复合材料(C/SiCS)以其低密度、抗氧化[3]、耐高温[4]等优异特性逐渐得到重视, 目前已经成功应用于飞行器的鼻锥和机翼前缘的防热结构中。然而, 与金属材料不同, 陶瓷基材料不具备产生塑性变形的能力, 材料的内部分层、基体裂纹、纤维断裂等破坏是其吸收冲击能量的主要方式, 这些损伤会降低材料的力学性能, 从而严重威胁材料的服役安全。因此, 为了飞行器内部结构及设备在服役过程中的安全性考虑, 这类结构材料在冲击载荷作用下的破坏行为及其数值预测一直是众多学者关注的热点。
目前, 针对平纹编织C/SiC(2D-C/SiC)复合材料的研究并不多见, 已有的相关文献多数仅限于对材料简单力学性能的实验研究[5-6]。杨扬等[7]基于电炮加载装置驱动Mylar飞片完成了冲击速度在3 400~9 300 m/s范围内2D-C/SiC材料的超高速撞击实验, 并对材料损伤演变和破坏特征进行了分析研究。然而, 关于材料在低速冲击载荷作用下的实验及其数值模拟目前尚未发现。
材料的各向异性特征使得在对其进行模拟计算时面临着巨大的挑战, R.A.Clegg等[8]提出的正交各向异性模型可以较好地描述复合材料在冲击载荷作用下的力学响应, 已经在Kevlar-Epoxy复合材料的冲击模拟中得到了成功的应用[9]。然而, 该模型重点关注了材料的非线性应变硬化行为, 在模拟C/SiC这种典型的脆性材料时无法得到理想的数值结果。本文中在已完成的球形弹丸低速冲击2D-C/SiC薄板的实验基础上, 引入材料的各向异性本构模型, 基于Autodyn商用软件对2D-C/SiC复合材料在钢球弹丸冲击作用下的数值模拟问题进行研究。
1. 2D-C/SiC材料低速冲击实验
基于空气炮加载装置对2D-C/SiC复合材料薄板在钢球撞击作用下的力学响应问题进行了实验研究。在实验过程中, 使用高速摄影技术记录了2D-C/SiC薄板受撞击后碎片云团的发展过程。
1.1 实验材料与计划
实验中选用的2D-C/SiC复合材料由西北工业大学超高温复合材料重点实验室研制。首先, T300碳纤维正交编织成碳布预制体, 在低压炉内基于CVI技术沉积6次SiC基体, 得到孔隙率约15%、体积分数约45%的2D-C/SiC材料, 而后基于CVI沉积2次SiC涂层。最终加工成型后, 切割为大小115 mm×115 mm、厚度3 mm、密度约2.11 g/cm3的正方形薄板作为目标靶板。实验中选用5、6 mm两种直径的普通钢球作为冲击弹丸。具体实验计划及主要结果列于表 1中, 表中D表示弹丸直径, vp表示弹丸的冲击速度, E表示弹丸的动能。
表 1 低速冲击实验计划Table 1. Low-speed impact experimental scheme实验编号 D/mm vp/(m·s-1) E/J 实验结果 1 5 79.0 1.61 未穿透 2 5 92.7 2.22 未穿透 3 5 98.5 2.51 未穿透 4 5 130.0 4.37 穿透 5 6 144.0 9.26 穿透 6 6 147.0 9.65 穿透 7 6 211.0 19.89 穿透 8 6 218.0 21.23 穿透 9 6 219.0 21.43 穿透 1.2 碎片云团结构
限于实验条件, 高速摄像机置于靶板斜后方位置, 图 1给出了试样8的高速摄影视频截图。由图 1可见, 2D-C/SiC碎片云团的发展过程大致分为2个阶段:(1)簇状飞散(图 1(a))。弹丸接触靶板时, 在作用点局部区域迅速发生纤维断裂、基体粉碎等多种形式破坏, 大量固体颗粒成簇状形式堆积飞出; (2)两区域形成(图 1(b)~(d))。簇状颗粒群运动发展, 逐渐在碎片云团中形成了界限清晰的两区域结构, 即柱状区和分散区。在柱状区内, 颗粒分布密度较高、颗粒尺寸较小、运动速度较快, 并沿轴线向前运动; 在分散区内, 颗粒分布较为稀疏、颗粒尺寸较大、运动速度较慢, 且沿一定的飞散角辐射飞出。2D-C/SiC碎片云结构与以往的金属材料及树脂基复合材料截然不同, 是其脆性特征的典型表现。在数值模拟中, 建立合理的本构模型对材料特有的力学响应特征进行模拟计算是本文研究的重点。
2. 耦合的各向异性材料模型
2.1 正交各向异性弹性模型
复合材料在冲击载荷作用下的力学响应过程十分复杂, 通常是诸如体积响应、材料各向异性等多种因素的耦合。在数值计算中, 不仅球应变会影响偏应力的计算, 偏应变同样会在球应力的计算中起作用, 通常需要面临求解状态方程和本构关系这两个高度耦合的子模型[10]。传统的计算方法是将这两个子模型独立计算, 然而这样的过程较为繁琐, C.E.Anderson等[11]首先提出了一种适用于求解各向异性材料中多重响应耦合问题的计算方法。本文中也以此方法为基础, 建立如下的计算模型。
2D-C/SiC是典型的正交各向异性复合材料。基于塑性力学, 全应变可分解为与体积变化相关的球应变ε0和与形状变化相关的偏应变
之和。此时, 2D-C/SiC的本构关系可用如下全应力应变关系表述为:
[σ11σ22σ33σ23σ31σ12]=[C11C12C13000C12C22C23000C13C23C33000000C44000000C55000000C66][ε0+εd11−1ε0+εd22ε0+εd33ε23ε31ε12] 由于各向同性压力P是3个正应力的平均值的相反数, 因此由上式可以推导得出:
P=−13(C11+C22+C33+2C12+2C13+2C23)ε0−3∑i=113[C1i+C2i+C3i]εdii (1) 从(1)式中可以清晰地看出, 与各向同性材料中球应力仅依赖于体积应变不同, 各向异性材料中偏应变也会在一定程度上影响球应力的计算。其中ε0一项的系数可称为材料的等效体积模量, 记为A。
2.2 正交各向异性失效模型
失效模型主要包括失效初始化和后失效响应两个阶段。在失效初始化阶段, 采用最大应力判断准则, 即预先设置材料在3个拉伸方向和3个剪切方向的失效应力, 当任一方向的应力值达到失效应力时, 材料即在该方向失效, 进入后失效响应阶段。
对于材料的后失效响应, 需要在计算开始前按照如下准则进行相关设置:当材料在任一方向进入后失效阶段时, 材料在该方向的承载能力瞬间丧失, 但在其他方向的承载能力保持不变; 同时, 材料在3个剪切方向的承载能力下降[8]。例如, 假设1方向为薄板厚度方向, 2、3方向为面内的2个正交方向。那么, 当1方向拉伸应力大于拉伸失效应力或12(或13)方向剪切应力大于剪切失效应力时, 材料会发生分层破坏。此时, 材料在1方向应力立即置为0, 不再具有承载能力, 但在面内的2个方向承载能力不变; 同时由于分层破坏的发生, 材料在剪切方向的承载能力也将有所下降, 需要对刚度矩阵中的剪切刚度系数进行退化, 可对材料的刚度矩阵作如下修正(α为刚度退化因子):
[0σ22σ33σ23σ31σ12]=[0000000C22C230000C23C33000000αC44000000αC55000000αC66][ε11ε22ε33ε23ε31ε12] 表 2 正交弹性模型材料参数Table 2. Parameters of orthotropic elastic model刚度系数 体积模量 C11/GPa C22/GPa C33/GPa C12/GPa C23/GPa C31/GPa G12/GPa G23/GPa G31/GPa A/GPa 19 120.5 120.5 7.8 20 7.8 8.8 23 8.8 36.8 表 3 正交失效模型材料参数Table 3. Parameters of orthotropic failure model失效初始化 σfail-11/MPa σfail-22/MPa σfail-33/MPa σfail-12/MPa σfail-23/MPa σfail-31/MPa 50 455 455 40 165 40 后失效响应 刚度退化因子 失效模式 0.2 11失效
11方向拉伸
应力置022失效
22方向拉伸
应力置033失效
33方向拉伸
应力置012失效
11方向拉伸
应力置023失效
11方向拉伸
应力置031失效
11方向拉伸
应力置03. 数值模拟及分析
3.1 计算模型
由于正撞击的对称性, 基于Autodyn软件建立如图 2所示的二维半模型计算。2D-C/SiC靶板选用SPH求解器; 钢弹丸steel 4340基于FEM求解器计算; 弹丸与靶板之间设置接触, 接触间距0.01 mm; 在靶板边缘3 mm×4 mm范围内施加全局固定边界, 即限定此区域内的粒子在x和y方向速度恒为0。
下面通过对碎片云结构和典型破坏形式等实验结果的数值验证, 证明该模型对描述2D-C/SiC材料在冲击载荷作用下力学响应的合理性和可行性。
3.2 数值合理性的实验对比
3.2.1碎片云团结构的数值验证
图 3 (b)、(c)分别给出了试样8的数值模拟结果与对应的速度矢量图。通过对比可以看出, 实验结果和数值模拟结果吻合较好, 具有以下相似特征:(1)碎片云团均呈现明显的双区域结构, 即均存在沿冲击轴线方向的柱状区和轴线两侧的飞散区, 这是2D-C/SiC在冲击载荷作用下的特殊响应特征; (2)从速度矢量图可以看出, 柱状区粒子运动方向均为沿轴线方向; 而飞散区内的粒子速度方向与轴线有大小不同的夹角。这与实验结果是完全一致的, 表明该模型可以较好地表征2D-C/SiC在弹丸撞击下碎片云团的结构特征。
3.2.2 B扫描结果的数值验证
B扫描结果显示的实际上是试样内部损伤的二维截面图。图 4(a)给出试样8的B扫结果, 其中虚线区域的黄色部分表示试样; 图 4(b)~(c)给出工况8靶板扩展稳定后内部损伤的模拟计算结果。
从图 4看出:(1)数值模拟与B扫描显示的穿孔形状十分相似, 穿孔直径近似为弹丸直径(6 mm); (2)从B扫结果看, 靶材背面显示黄色较正面偏淡, 这是由于材料在靠近背面的区域损伤严重, 导致超声波信号无法完全被探头接收。相似的现象可以在数值模拟中显示得更加清晰, 图 4(c)给出了靶材穿孔区域的局部放大图。由图可见, 与正面相比, 靶材背面存在明显的材料分层与剥落损伤, 这与实验结果比较吻合。上述相似性也表明该模型可以较好地表征2D-C/SiC在弹丸撞击下的破坏特征。
3.2.3碎片云轴向速度的数值验证
碎片云的轴向速度是表征碎片云团发展状态的重要指标。实验中, 在靶板表面标记了刻度, 据此可以根据高速摄影记录计算出碎片云团柱状区前端的运动速度, 而将此速度视为碎片云轴向运动速度vaxis。同样, 在数值模拟中, 取靶板自由面与冲击轴线的交点为特征点近似柱状区端点。选取靶板穿透工况中冲击速度相差较大的4种情况进行验证, 图 5给出了不同冲击速度vaxis下随时间的变化历程。
从图 5可以看出, 在弹丸接触靶板后的短时间内, 特征点的轴向速度迅速增加并逐渐趋于稳定。将特征点最终的稳定速度视为对vaxis的近似, 图 6给出了实验数据与数值计算的对比结果。可以看出, 二者的吻合度较好, 并且随着冲击速度的增高, 相对误差逐渐减小, 其中冲击速度219 m/s时二者的相对误差仅为2.25%, 表明该模型可以较好地表征2D-C/SiC在冲击载荷作用下碎片云团的速度特性。
4. 极限侵彻深度预测模型
尽管碎片云团中柱状区内的粒子运动速度较高, 但由于颗粒质量较小, 对机体的损伤程度有限。碎片云团中对机体威胁较大的仍是穿透靶板的弹丸, 因此可以通过对弹丸剩余速度的分析研究, 建立合理有效的预测公式评估C/SiC靶板的抗冲击性能。以上通过3方面对比验证了文中使用的材料模型及相关参数的正确性和合理性, 下面将基于此模型对弹丸侵彻下C/SiC靶板的抗冲击特性进行预估。
考虑同一弹丸(直径6 mm)以不同冲击速度对不同厚度靶板进行撞击。设置靶板厚度δ为1、2、3、4、5 mm, 弹丸冲击速度vi为50、100、150、200、250、300、350、400 m/s, 计算得到弹丸剩余速度变化规律如图 7所示。可以看出:(1)对同一厚度的靶板, 剩余速度随着冲击速度的减小不断减小; 而对相同冲击速度, 剩余速度会随着靶板厚度的减小而不断增加; (2)剩余速度在速度较低靶板较厚时出现了负值, 这是由于弹丸在这些工况下未能穿透靶板, 出现回弹现象造成的。
评估材料抗冲击性能比较关心的重要参数是靶板的极限穿透速度, 即上述曲面与剩余速度为0的平面的交线。然而, 如果固定靶板厚度分别观察剩余速度随冲击速度的变化趋势(图 8), 可以看出在临界穿透状态附近(图中虚线), 剩余速度变化曲线斜率陡然增高, 说明在临界穿透状态附近剩余速度变化较剧烈; 而当冲击速度高于极限穿透速度之后, 弹丸剩余速度变化基本稳定, 呈准线性增长趋势。因此, 如果直接对剩余速度曲面进行拟合来求解极限穿透速度, 会导致在临界穿透状态附近的计算误差较大。
为了解决上述问题, 采用分步拟合方法:(1)假设靶板无限厚, 首先固定冲击速度, 分别拟合得到弹丸剩余速度随靶板厚度的关系, 进而求解在此冲击速度下弹丸的极限侵彻深度; (2)其次基于Poncelet型侵彻深度预测模型[12], 利用得到的极限侵彻深度T与弹丸冲击速度vi进行拟合, 得到如下关系:
T=0.0005362ln(1+0.756v2i) 此时, 便可以在已知靶板厚度的情况下, 通过反解上式预测弹丸的极限穿透速度。同时, 在工程设计中, 如果工程人员已经根据飞行器的服役环境预估出外来冲击物的撞击速度, 那么同样可以利用上式预估弹丸的极限侵彻深度, 也就可以为防护板厚度参数的确定提供一定的参考。
5. 结论
基于Autodyn软件选取正交各向异性材料模型, 推导了相关的材料参数, 并从碎片云团结构、B扫描内部损伤形貌和碎片云轴向速度三个方面对比了实验与模拟结果, 表明参数的合理性和正确性, 较好地描述了2D-C/SiC在钢球弹丸撞击下的主要损伤特征和力学特性。在此基础上, 基于数值模拟计算结果推导得到了钢球弹丸撞击2D-C/SiC材料的极限侵彻深度预测公式, 这可为TPS结构设计提供一定的工程参考价值。
-
表 1 不同应变率下人工软骨仿生超材料的拉伸强度和断裂应变
Table 1. Tensile strength and breaking strain of the artificial cartilage foam material in tensile experiments at different strain rates
应变率/s−1 拉伸强度/MPa 断裂应变 10−1 0.910 1.468 40 2.046 1.381 120 2.895 1.131 160 2.969 0.772 表 2 不同应变率下压缩实验的平台应力和弹性模量
Table 2. Plateau stresses and elastic moduli obtained in compression experiments at different strain rates
应变率/s−1 弹性模量/MPa 平台应力/MPa 10−1 1.006 0.30 50 23.896 2.22 100 34.069 3.75 150 74.522 4.50 表 3 ACF和EPP材料在不同冲击次数下的能量吸收
Table 3. Energy absorption of ACF and EPP materials under different impact times
No. ACF材料 EPP材料 吸能/J 吸能率/% 吸能率下降/% 吸能/J 吸能率/% 吸能率下降/% 1 44.7 97.1 31.5 71.6 2 45.2 96.6 0.5 26.2 60.0 11.6 3 45.0 96.4 0.7 25.6 58.2 13.4 4 45.0 96.7 0.4 23.8 52.5 19.1 5 45.1 96.7 0.4 25.1 56.3 15.3 -
[1] 王志华, 李世强, 李鑫, 等. 轻质多孔金属及其夹芯结构力学行为的研究进展 [J]. 太原理工大学学报, 2017, 48(3): 492–503. DOI: 10.16355/j.cnki.issn1007-9432tyut.2017.03.030.WANG Z H, LI S Q, LI X, et al. Advances in studies of the mechanical behavior of cellular metallic materials and structures [J]. Journal of Taiyuan University of Technology, 2017, 48(3): 492–503. DOI: 10.16355/j.cnki.issn1007-9432tyut.2017.03.030. [2] 刘培生, 崔光, 程伟. 多孔材料性能模型研究1: 数理关系 [J]. 材料工程, 2019, 47(6): 42–62. DOI: 10.11868/j.issn.1001-4381.2018.001407.LIU P S, CUI G, CHENG W. Study on property model for porous materials 1: mathematical relations [J]. Journal of Materials Engineering, 2019, 47(6): 42–62. DOI: 10.11868/j.issn.1001-4381.2018.001407. [3] GIBSON L J, ASHBY M F. Cellular solids: structure and properties [M]. Cambridge, UK: Cambridge University Press, 1997. [4] WANG N Z, CHEN X, AO L I, et al. Three-point bending performance of a new aluminum foam composite structure [J]. Transactions of Nonferrous Metals Society of China, 2016, 26(2): 359–368. DOI: 10.1016/s1003-6326(16)64088-8. [5] WAN T, LIU Y, ZHOU C X, et al. Fabrication, properties, and applications of open-cell aluminum foams: a review [J]. Journal of Materials Science & Technology, 2021, 62(3): 11–24. DOI: 10.1016/j.jmst.2020.05.039. [6] ZHANG J X, YANG Y, QIN Q H, et al. Dynamic compressive response of sinusoidal corrugated core sandwich plates [J]. International Journal of Applied Mechanics, 2018, 10(7): 1850075. DOI: 10.1142/S1758825118500758. [7] 朱凌, 郭开岭, 余同希, 等. 泡沫金属夹芯梁在重复冲击下的动态响应 [J]. 爆炸与冲击, 2021, 41(7): 073101. DOI: 10.11883/bzycj-2020-0198.ZHU L, GUO K L, YU T X, et al. Dynamic responses of metal foam sandwich beams to repeated impacts [J]. Explosion and Shock Waves, 2021, 41(7): 073101. DOI: 10.11883/bzycj-2020-0198. [8] 王飞, 李剑荣, 虞吉林. 铝蜂窝结构单向压缩、失稳和破坏机制研究 [J]. 力学学报, 2001, 33(6): 741–748. DOI: 10.3321/j.issn:0459-1879.2001.06.003.WANG F, LI J R, YU J L. A study of instability and collapse of aluminum honeycombs under uniaxial compression [J]. Chinese Journal of Theoretical and Applied Mechanics, 2001, 33(6): 741–748. DOI: 10.3321/j.issn:0459-1879.2001.06.003. [9] 张健, 赵桂平, 卢天健. 泡沫金属在冲击载荷下的动态压缩行为 [J]. 爆炸与冲击, 2014, 34(3): 278–284. DOI: 10.11883/1001-1455(2014)03-0278-07.ZHANG J, ZHAO G P, LU T J. High speed compression behaviour of metallic cellular materials under impact loading [J]. Explosion and Shock Waves, 2014, 34(3): 278–284. DOI: 10.11883/1001-1455(2014)03-0278-07. [10] 杨宝, 汤立群, 刘逸平, 等. 冲击条件下泡沫铝的细观变形特征分析 [J]. 爆炸与冲击, 2012, 32(4): 399–403. DOI: 10.11883/1001-1455(2012)04-0399-05.YANG B, TANG L Q, LIU Y P, et al. Meso deformation characteristics analysis of aluminum foam under impact [J]. Explosion and Shock Waves, 2012, 32(4): 399–403. DOI: 10.11883/1001-1455(2012)04-0399-05. [11] 李玲, 钱学登, 王科. 复合结构混凝土动力学特征及变形破坏机制 [J]. 煤矿安全, 2019, 50(4): 41–45. DOI: 10.13347/j.cnki.mkaq.2019.04.010.LI L, QIAN D K, WANG K. Dynamic characteristics and deformation failure mechanism of composite structure concrete [J]. Safety in Coal Mines, 2019, 50(4): 41–45. DOI: 10.13347/j.cnki.mkaq.2019.04.010. [12] 韩李斌, 杨黎明. 泡沫混凝土动态力学性能及破坏形式 [J]. 宁波大学学报(理工版), 2017, 30(1): 68–72. DOI: 10.3969/j.issn.1001-5132.2017.01.013.HANG L B, YANG L M. Dynamic properties and failure types of foamed concrete [J]. Journal of Ningbo University (Natural Science and Engineering), 2017, 30(1): 68–72. DOI: 10.3969/j.issn.1001-5132.2017.01.013. [13] 王必勤. EPDM发泡材料的制备及结构与性能研究[D]. 上海: 上海交通大学, 2006.WANG B Q. Preparation of EPDM and study on its structure and properties [D]. Shanghai, China: Shanghai Jiao Tong University, 2006. [14] 景鹏. 高g值冲击测试关键技术研究[D]. 太原: 中北大学, 2009.JING P. High-g impact test value of the research of key technologies [D]. Taiyuan, Shanxi, China: North University of China, 2009. [15] 鲁林, 李晓峰. 冲击环境作用下聚氨酯材料的应变率分布及吸能特性研究 [J]. 兵工学报, 2015, 36(S1): 213–219.LU L, LI X F. Research on strain rate distribution and energy absorption of polyurethane materials under shock environment [J]. Acta Armamentarii, 2015, 36(S1): 213–219. [16] 杨文叶, 姜子敬, 闫雪燕, 等. EPP材料在汽车座椅中的应用和性能研究 [J]. 汽车文摘, 2021(8): 58–62. DOI: 10.19822/j.cnki.1671-6329.20210062.TANG W Y, JINAG Z J, YAN X Y, et al. Study on the application and performance of EPP materials in car seats [J]. Automotive Digest, 2021(8): 58–62. DOI: 10.19822/j.cnki.1671-6329.20210062. [17] 金强维. EPP动态缓冲性能的研究[D]. 西安: 陕西科技大学, 2019.JIN Q W. Investigation on dynamic cushioning property of EPP [D]. Xi’an, Shaanxi, China: Shaanxi University of Science and Technology, 2019. [18] 陈润峰, 石庚辰, 张力, 等. 人工软骨仿生材料在引信缓冲中的应用 [J]. 探测与控制学报, 2020, 42(6): 10–14.CHEN R F, SHI G C, ZHANG L, et al. Application of artificial cartilage biomimetic material in fuze buffer [J]. Journal of Detection & Control, 2020, 42(6): 10–14. -