Injury properties of porcine lung under blast load
-
摘要: 为了深入研究爆炸冲击波作用下生物体肺部的力学响应和损伤特性,首先建立了猪胸部有限元模型,借助新研制的PVDF (polyvinylidene fluoride)柔性冲击波压力传感器测试了激波管试验中动物的体表压力,验证了有限元模型的准确性。然后,使用已验证的模型开展了不同比例距离下猪肺部损伤特性研究,分析了在不同强度冲击波作用下肺部的损伤程度和损伤区域,并建立了胸肺部表皮压力峰值与肺损伤的关系。最后,通过开展爆炸试验,获得了不同比例距离下猪的肺部损伤情况和胸部表皮压力曲线,验证了所建立的胸肺部表皮压力峰值与肺损伤关系的正确性。Abstract: In order to study the mechanical response and injury characteristics of living lungs under the action of blast shock waves, this study first established a finite element model of the chest of a small pig and used a newly developed PVDF (polyvinylidene fluoride) flexible pressure sensor to test the surface pressure of animals in shock tube experiments, verifying the accuracy of the finite element model. Secondly, the validated model was used to study the characteristics of lung injury of pigs at different blast distances, analyzing the lung injury location and severity under different intensities of shock waves. Moreover, the relation between the peak pressure on the surface of the chest and the grade of lung injury was established. Finally, blast tests were performed to obtain the lung injury of pigs as well as the chest surface pressure histories at different blast distances, which can be used to verify the correctness of the established relation.
-
Key words:
- animal finite element model /
- blast shock wave /
- lung injury /
- peak pressure
-
高超声速飞行器热防护系统(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. Bone material parameters of chest model
骨骼结构 密度/(kg∙m−3) 杨氏模量/MPa 泊松比 屈服强度/MPa 胸椎骨 2500 11000 0.40 — 胸椎间盘 1040 300 0.40 — 胸骨皮质骨 2 000 14000 0.30 90 肋间软骨 1500 12.5 0.40 4.9 肋骨松质骨 1000 40 0.45 1.8 肋骨皮质骨 2 000 11500 0.30 90 胸骨松质骨 1000 40 0.45 1.8 表 2 胸部模型软组织材料参数
Table 2. Soft tissue material parameters of chest model
组织 密度/(kg∙m−3) 线性体积模量/MPa 阻尼系数 剪切模量/MPa 心脏 1000 100 0.1 0 肺 288 150 0.4 0.04 胃 1000 4.59 0.1 0 肌肉 1000 2 000 0.4 0.04 表 3 不同比例距离下冲击波的等效参数
Table 3. The equivalent parameters of the shock wave at each proportional distance
序号 WTNT/kg R/m Z/(m∙kg−1/3) Δpm/kPa τ+/ms 1 1.08 2.30 2.27 150 2 2 1.30 2.17 2.00 200 2 3 1.70 2.00 1.66 300 2 4 2.00 1.86 1.48 400 2 表 4 胸廓运动速度、猪胸肺部表皮压力峰值与肺部损伤的对应关系
Table 4. Relations between thoracic motility velocity, the peak pressure on the chest and the lung injury in pigs
v/(m∙s−1) [30] p/kPa 肺损伤 <3.6 <232 无伤 3.6~7.5 232~609 微伤至轻伤 4.3~9.8 294~865 轻伤至中伤 7.5~16.9 609~ 1766 中伤至重伤 >12.8 > 1227 50%死亡 表 5 爆炸试验工况
Table 5. Blast test record
试验编号 猪编号 猪质量/kg TNT质量/kg 爆距/m 比例距离/(m∙kg−1/3) 1 1# 24.0 10 5.0 2.32 2# 25.5 10 6.0 2.78 3# 27.0 10 7.0 3.25 2 4# 27.5 10 3.5 1.62 5# 24.5 10 4.0 1.86 6# 24.0 10 4.0 1.86 3 7# 24.0 10 2.5 1.16 8# 25.5 10 3.0 1.39 9# 26.0 10 4.0 1.86 表 6 各个工况下猪肺部解剖伤情
Table 6. Anatomical injuries of pig lungs at various blast distances
爆距/m 肺部解剖伤情 解剖结果(AIS评估) 7.0 没有明显损伤 无伤 6.0 轻微损伤 轻微肺损伤 5.0 轻微损伤 轻微肺损伤 4.0 气管没有异常,左肺出血,右肺出血较严重 轻度肺损伤 3.0 气管没有异常,左肺出血较严重,右肺出血严重 中度到重度肺损伤 2.5 气管中有大量血沫,左肺出血较严重,右肺出血非常严重 重度肺损伤 表 7 猪肺伤情汇总与损伤判据判断伤情
Table 7. Summary of pig lung injury and the injury criterion determines the severity of the injury
爆距/m 体表压力峰值/kPa 解剖结果(AIS评估) 预测结果(本文损伤关系) 7.0 128 无伤 无伤 6.0 270 轻微肺损伤 微伤至轻伤 5.0 403 轻微肺损伤 微伤至轻伤 4.0 559 轻度肺损伤 轻伤至中伤 3.0 1253 中度到重度肺损伤 中伤至重伤 2.5 1600 重度肺损伤 50%死亡 -
[1] HUGHES S M, BORDERS III C W, ADEN J K, et al. Long-term outcomes of thoracic trauma in U. S. service members involved in combat operations [J]. Military Medicine, 2020, 185(11/12): 131–136. DOI: 10.1093/milmed/usaa165. [2] SALEM M H, BUX G M K, HUDEEL A. War thoracic wounds among civilians casualties in Aden during the 2015 [J]. Electronic Journal of University of Aden for Basic and Applied Sciences, 2020, 1(3): 159–166. DOI: 10.47372/ejua-ba.2020.3.39. [3] ZHANG B. Blast lung injury [M]//WANG Z G, JIANG J X. Explosive Blast Injuries: Principles and Practices. Singapore: Springer, 2023: 295–300. DOI: 10.1007/978-981-19-2856-7_19. [4] CERNAK I, SAVIC J, IGNJATOVIC D, et al. Blast injury from explosive munitions [J]. The Journal of Trauma: Injury, Infection, and Critical Care, 1999, 47(1): 96–103. DOI: 10.1097/00005373-199907000-00021. [5] SCOTT T E, KIRKMAN E, HAQUE M, et al. Primary blast lung injury: a review [J]. British Journal of Anaesthesia, 2017, 118(3): 311–316. DOI: 10.1093/bja/aew385. [6] SMITH J E. The epidemiology of blast lung injury during recent military conflicts: a retrospective database review of cases presenting to deployed military hospitals, 2003–2009 [J]. Philosophical Transactions of the Royal Society B: Biological Sciences, 2011, 366(1562): 291–294. DOI: 10.1098/rstb.2010.0251. [7] TELAND J A. Review of blast injury prediction models: FFI-rapport 2012/00539 [R]. Norway: Norwegian Defence Research Establishment (FFI), 2012. [8] COOPER G J, TOWNEND D J, CATER S R, et al. The role of stress waves in thoracic visceral injury from blast loading: modification of stress transmission by foams and high-density materials [J]. Journal of Biomechanics, 1991, 24(5): 273–285. DOI: 10.1016/0021-9290(91)90346-O. [9] 唐献述, 王树民, 龙源, 等. 爆炸空气冲击波对动物伤害效应试验研究 [J]. 工程爆破, 2012, 18(2): 104–106, 96. DOI: 10.3969/j.issn.1006-7051.2012.02.028.TANG X S, WANG S M, LONG Y, et al. Experimental study on the effect of explosion air shockwave on the animal injury [J]. Engineering Blasting, 2012, 18(2): 104–106, 96. DOI: 10.3969/j.issn.1006-7051.2012.02.028. [10] 王海宾, 赵英虎, 高莉, 等. 甲烷爆炸冲击波作用下密闭管道内动物损伤效应试验研究 [J]. 兵工学报, 2018, 39(8): 1639–1647. DOI: 10.3969/j.issn.1000-1093.2018.08.022.WANG H B, ZHAO Y H, GAO L, et al. Experimental investigation into the damage effect of methane explosion shock wave on animals in enclosed pipeline [J]. Acta Armamentarii, 2018, 39(8): 1639–1647. DOI: 10.3969/j.issn.1000-1093.2018.08.022. [11] VASSOUT P, FRANKE R, PARMENTIER G, et al. Mesures de pression et d’accélérations intracorporelles chez le porc exposé à des ondes de choc fortes en champ libre [R]. France: French-German Research Institute of Saint Louis, 1986: 86. [12] 陈海斌, 王正国, 杨志焕, 等. 冲击波传播的三个时段模拟实验中动物肺的损伤 [J]. 爆炸与冲击, 2000, 20(3): 264–269. DOI: 10.11883/1001-1455(2000)03-0264-6.CHEN H B, WANG Z G, YANG Z H, et al. Injury of animal lungs in the experiments to simulate the three phases of shock wave propagation [J]. Explosion and Shock Waves, 2000, 20(3): 264–269. DOI: 10.11883/1001-1455(2000)03-0264-6. [13] 陈海斌, 王正国. 爆炸性减压对兔肺的损伤作用 [J]. 中华创伤杂志, 2000, 16(2): 109–111. DOI: 10.3760/j:issn:1001-8050.2000.02.016.CHEN H B, WANG Z G. Injury action on rabbit lung by explosive decompression [J]. Chinese Journal of Trauma, 2000, 16(2): 109–111. DOI: 10.3760/j:issn:1001-8050.2000.02.016. [14] 段维勋. 胸部爆炸伤动物模型的建立及伤情特点分析的实验研究 [D]. 西安: 第四军医大学, 2002: 1–64.DUAN W X. Experimental researches on establishment of animal model of thoracic explosive injury and analysis of wound characteristics [D]. Xi’an: Air Force Medical University, 2002: 1–64. [15] 王峰, 杨志焕, 朱佩芳, 等. 高原冲击伤伤情特点的实验研究 [J]. 创伤外科杂志, 2008, 10(6): 549–551. DOI: 10.3969/j.issn.1009-4237.2008.06.026.WANG F, YANG Z H, ZHU P F, et al. Experimental study on characteristics of blast injury at high altitude [J]. Journal of Traumatic Surgery, 2008, 10(6): 549–551. DOI: 10.3969/j.issn.1009-4237.2008.06.026. [16] 袁丹凤, 杨傲, 麻超, 等. 冲击波强度与幼年大鼠肺冲击伤程度的量效关系 [J]. 中国医学物理学杂志, 2021, 38(6): 780–784. DOI: 10.3969/j.issn.1005-202X.2021.06.022.YUAN D F, YANG A, MA C, et al. Dose-effect relationship between shock wave intensity and blast lung injury in juvenile rats [J]. Chinese Journal of Medical Physics, 2021, 38(6): 780–784. DOI: 10.3969/j.issn.1005-202X.2021.06.022. [17] 张良. 准静态撞击下猪胸腔内应力分布及其三维重现与肺损伤 [D]. 重庆: 中国人民解放军陆军军医大学, 2005: 1–56.ZHANG L. Distribution and 3-D reconstruction of intrathoracic stress in the swine chest subjected to impact [D]. Chongqing: Army Medical University, 2005: 1–56. [18] 杨春霞. 羊肺脏有限元模型的建立及其在冲击波作用下的仿真分析 [D]. 重庆: 重庆大学, 2010: 1–60.YANG C X. Establishmen of finite element model of sheep lung and finite-element simulation of blast wave [D]. Chongqing: Chongqing University, 2010: 1–60. [19] MORDAKA J, MEIJER R, VAN ROOIJ L, et al. Validation of a finite element human model for prediction of rib fractures: SAE technical paper 2007-01-1161 [R]. TNO Science and Industry: SAE, 2007. [20] KIMPARA H, LEE J B, YANG K H, et al. Development of a three-dimensional finite element chest model for the 5th percentile female [R]. Wayne State University, Toyota Central R & D Labs., Inc.: SAE, 2005. [21] SHIN J, UNTAROIU C, LESSLEY D, et al. Thoracic response to shoulder belt loading: investigation of chest stiffness and longitudinal strain pattern of ribs: SAE technical paper 2009-01-0384 [R]. Center for Applied Biomechanics, University of Virginia: SAE, 2009. DOI: 10.4271/2009-01-0384. [22] IWAMOTO M, NAKAHIRA Y, TAMURA A, et al. Development of advanced human models in THUMS [C]//Proceedings of the 6th European LS-DYNA Users’ Conference. Nagakute City: Toyota Central R&D Labs., Inc., 2007. [23] YANG K H, HU J W, WHITE N A, et al. Development of numerical models for injury biomechanics research: a review of 50 years of publications in the Stapp car crash conference [C]//Proceedings of the 50th Stapp Car Crash Conference. Dearborn, Michigan: SAE, 2006: 429–490. [24] RUAN J, EL-JAWAHRI R, CHAI L, et al. Prediction and analysis of human thoracic impact responses and injuries in cadaver impacts using a full human body finite element model [C]//Proceedings of the 47th Stapp Car Crash Conference. San Diego: SAE, 2003. [25] ROBIN S. HUMOS: human model for safety: a joint effort towards the development of refined human-like car occupant models [C]//Proceedings of the 17th International Technical Conference on the Enhanced Safety of Vehicles. Amsterdam: SAE, 2001. [26] RUAN J S, EL-JAWAHRI R, BARBAT S, et al. Biomechanical analysis of human abdominal impact responses and injuries through finite element simulations of a full human body model [C]//Proceedings of the 49th Stapp Car Crash Conference. National Library of Medicine: SAE, 2005. [27] 范志强, 常瀚林, 何天明, 等. 基于PVDF复合压电效应的低强度冲击波柔性测量 [J]. 爆炸与冲击, 2023, 43(1): 013102. DOI: 10.11883/bzycj-2022-0152.FAN Z Q, CHANG H L, HE T M, et al. Flexible measurement of low-intensity shock wave based on coupling piezoelectric effect of PVDF [J]. Explosion and Shock Waves, 2023, 43(1): 013102. DOI: 10.11883/bzycj-2022-0152. [28] GREER A. Numerical modeling for the prediction of primary blast injury to the lung [D]. Waterloo: University of Waterloo, 2007: 117–141. [29] 卢芳云, 蒋邦海, 李翔宇, 等. 武器战斗部投射与毁伤 [M]. 北京: 科学出版社, 2013: 1–316. [30] AXELSSON H, YELVERTON J T. Chest wall velocity as a predictor of nonauditory blast injury in a complex wave environment [J]. The Journal of Trauma: Injury, Infection, and Critical Care, 1996, 40(3S): 31S–37S. DOI: 10.1097/00005373-199603001-00006. -