基于爆炸损伤的头部有限元模型建立与验证

李涛; 常利军; 陈泰伟; 刘浚源; 肖淞铭; 蔡志华

doi:10.11883/bzycj-2024-0173

基于爆炸损伤的头部有限元模型建立与验证

doi: 10.11883/bzycj-2024-0173

湖南科技大学机械设备健康维护湖南省重点实验室，湖南湘潭411201

基金项目: 国家自然科学基金（12372356，11972158）；湖南省研究生科研创新项目（CX20221044）

详细信息

作者简介:
李　涛（2000－　），男，硕士研究生，litao200530@163.com

通讯作者:
蔡志华（1981－　），男，博士，教授，博士生导师，caizhihua003@163.com

中图分类号: O389
计量
- 文章访问数: 152
- HTML全文浏览量: 93
- PDF下载量: 78
- 被引次数: 4
出版历程
- 收稿日期: 2024-06-05
- 修回日期: 2024-10-23
- 网络出版日期: 2024-11-05
- 刊出日期: 2024-12-01

Establishment and verification of a head finite element model based on explosion injury

Hunan Provincial Key Laboratory of Health Maintenance for Mechanical Equipment, Hunan University of Science and Technology, Xiangtan 411201, Hunan, China

摘要

摘要: 为了更好地理解爆炸冲击波作用下头部的力学响应和损伤机制，利用计算机电子断层扫描与核磁共振医学图像获取了头部几何信息，开发了具有骨缝结构的精细化头部有限元模型。基于已有的激波管尸体实验，开展了正面、侧面与背面爆炸冲击数值模拟，通过对比颅内压-时间历程曲线与颅内压峰值，验证有限元模型的有效性。结果表明：在3种冲击方向下，颅内4个区域的压力峰值与文献实验仿真数据吻合较好；爆炸仿真中颅骨骨缝处有明显应力集中，骨缝线处头部有更大的损伤风险；同等爆炸冲击强度下，正面和背面冲击比侧面冲击对头部造成的损伤更严重。建立的头部模型可应用于爆炸载荷下的头部损伤研究，同时可探究骨缝对于头部生物力学响应的影响，对爆炸损伤研究具有重要意义。
- 爆炸 /
- 有限元模型 /
- 模型验证 /
- 生物力学 /
- 颅内压力
Abstract: In order to better understand the mechanical response and injury mechanism of the head under the action of explosive shock wave, the geometric information of the head was obtained through computerized tomography and magnetic resonance imaging, and a finite element model of the head with fine cranial bone and brain tissue was developed. Based on the existing blast tube cadaver experiments, forward, side, and backward explosive shock numerical simulations were conducted, and the cranial pressure-time history curves and peak cranial pressure were compared to validate the finite element model. The results show that the peak pressures of the four regions in the cranium under the three impact directions are in good agreement with the experimental and simulated data in the literature; there is obvious stress concentration at the suture line of the cranial bone in the simulated blast simulation; the head has a greater risk of injury at the suture line; and the front and back impacts cause more serious head injuries than the side impact under the same explosive shock intensity. The head model established can be used in the study of head injury under explosive loading, and the influence of suture on the biomechanical response of the head can be explored, which has important research significance for blast injury research.
- explosion /
- finite element model /
- model verification /
- biomechanics /
- intracranial pressure

HTML全文

近年来，随着高超声速飞行器的快速发展，与飞行平台共形的异构型战斗部对目标的侵彻行为引起了广泛关注。与传统回转体战斗部不同，异构型战斗部能更好地适应高超声速平台的气动外形，提高舱内空间的利用率，进而提升武器的毁伤效能。受弹体姿态和环境等因素的影响，弹体着靶时通常存在攻角和着角，使得弹体对目标的侵彻特性更加复杂^[1]。由于异构型战斗部结构的特殊性，在斜侵彻过程中更容易出现姿态偏转和弹道偏离现象。因此，在一定约束条件下，确定异构型战斗部（以下简称弹体）在侵彻过程中的结构响应和失效规律，对于战斗部的结构设计具有重要意义。

针对传统回转体弹体斜侵彻典型靶标的作用过程与机理，已有丰富的研究成果。对于半无限靶，以空腔膨胀理论为基础，提出了描述弹体斜侵彻行为的理论模型^[2-4]，能够对侵彻弹道、弹体姿态等进行预测；闪雨^[5]基于微分面力法、空腔膨胀理论以及自由面效应模型，对弹体侵彻过程进行简化，建立了弹体侵彻半无限靶的弹道预测模型。对于多层钢板，Goldsmith^[6]对非理想条件下弹靶作用的研究成果进行总结，提出可从弹体剩余速度、靶板破坏模式和弹体姿态偏转等方面开展研究。Gupta等^[7]、Iqbal等^[8-9]、杜华池等^[10]开展了弹体斜侵彻多层钢板靶的试验和数值模拟研究，得到了弹体材料、入射姿态、靶板材料及厚度等因素对弹道偏转、剩余速度和靶板破坏模式的影响规律。

在异构型弹体方面，王文杰等^[11]、Dong等^[12]开展了椭圆截面弹体侵彻混凝土半无限靶的研究，建立了刚性弹体的侵彻深度和侵彻阻力理论模型；刘子豪^[13]、Dai等^[14]结合椭圆截面弹体的结构参数，对混凝土半无限靶的侵彻过程进行分析，发现弹体侵彻能力与长短轴之比呈正相关。王浩等^[15]、田泽等^[16]结合靶板破坏模式和弹体的受力特征，建立了椭圆截面弹体斜侵彻双层钢板过程中弹体姿态的偏转模型。岳胜哲等^[17]对类椭圆截面弹体斜侵彻薄靶的姿态偏转进行了研究，发现在相同条件下，弹体姿态偏转幅度与截面不对称度呈正相关。王景琛等^[18]针对异构型弹体斜侵彻薄靶问题进行数值仿真，分析了着角、攻角和弹靶结构参数对弹体所受动态载荷的影响规律。

以上研究主要针对弹体的侵彻性能、弹道特性和靶板失效等问题，没有考虑弹体的结构响应和失效。大量实验结果^[19]表明，随着侵彻速度的提高，弹体结构将出现头部侵蚀、弯曲、断裂等失效行为。针对头部侵蚀现象，Silling等^[20]基于对试验结果的分析，认为着靶速度接近1 km/s时，弹体质量损失与初始动能成正比；何翔等^[21]、武海军等^[22]、He等^[23]、Zhao等^[24]提出了一系列与速度相关的质量损失（半）理论模型。针对弹体弯曲和断裂现象，陈小伟等^[25-26]给出了深侵彻弹体壁厚的设计依据，并结合实例进行了验证。皮爱国等^[27]开展了大长径比弹体侵彻混凝土靶的结构响应研究，给出了弹体剪力、弯矩的计算方法以及弹体弯曲的临界条件。王一楠等^[28]、张欣欣等^[29]基于自由梁理论，分析了弹体斜侵彻混凝土靶的结构弯曲现象，得到剪力、弯矩以及屈服函数的分布规律。刘坚成等^[30]通过反弹道试验研究了大长径比弹体侵彻厚靶的结构响应特性，发现弹体变形集中在头部，以弯曲和屈曲为主。

综上所述，目前对弹体结构响应的研究主要针对混凝土厚靶的侵彻过程，且仅考虑了弹体头部受载情况。事实上，在侵彻过程中，弹身同样会受到严重的冲击载荷作用。特别是弹体斜侵彻多层钢板时，由于弹靶接触位置随时间发生变化，弹体的结构响应更为复杂，现有基于头部受载的结构响应结论无法有效推广应用。因此，为了深入探究斜侵彻多层钢板过程中弹体的结构响应，本文中设计圆形、椭圆、非对称椭圆三种截面弹体，开展不同弹体斜侵彻双层钢板的试验研究，并进行数值仿真和弹体结构响应模型的建立，着重分析斜侵彻条件下弹体的动态载荷特性以及结构失效条件。

1. 弹体斜侵彻双层钢板的试验研究

1.1 试验准备

根据异构型战斗部的结构特点，设计了圆形截面（circular cross-section，CC）、椭圆截面（elliptical cross-section，EC）和非对称椭圆截面（asymmetric elliptical cross-section，AC）三种类型弹体。圆形截面弹体的头部形状系数η和非对称椭圆截面弹体的不对称度γ分别定义如下：

$\eta = \frac{\rho }{D}\text{，}\gamma = \frac{{{b_1}}}{{{b_2}}}$

(1)

式中：ρ为圆形截面弹体头部子午线曲率半径，D为弹体直径，b₁和b₂分别为上下非对称椭圆截面弹体的上、下两部分短半轴长度，a为弹体截面长半轴长度，如图1所示。

图 1 弹体结构参数示意图

Figure 1. Schematic diagram of projectile structural parameters

下载: 全尺寸图片幻灯片

基于圆形截面弹体，对长轴和短轴进行适当比例的缩放，可得到椭圆和非对称椭圆截面弹体，通过控制参数保证三种弹体的截面积一致。三种弹体的具体结构参数如表1所示，表中D为圆形截面弹体直径，A和B为椭圆截面弹体长轴、短轴长度，L为弹体总长度，h为弹体壁厚，m为弹体质量，其中非对称椭圆截面弹体不对称度为2。通过数控中心加工得到的弹体实物如图2所示。

表 1 三种截面弹体的结构参数

Table 1. Structural parameters of three projectiles with different cross-sections

截面形状	截面参数/mm			L/mm	h/mm	m/g
截面形状	D	A	B	L/mm	h/mm	m/g
圆形	30	−	−	180	4	506.2
椭圆形	−	33	27	180	4	509.2
非对称椭圆形	−	33	18/9	180	4	519.5

下载: 导出CSV

| 显示表格

图 2 加工后的三种截面弹体

Figure 2. Three types of projectiles after manufacture

下载: 全尺寸图片幻灯片

为了保证弹体的结构强度，弹体材料选用30CrMnSiNi2A高强度钢。粗加工后，对试件进行一次淬火，然后再精加工。图3给出了弹体材料的准静态拉伸试验结果，其弹性模量为210 GPa，屈服强度和极限强度分别为1707、2046 MPa，断裂应变为5%。靶板材料为45钢，尺寸为520 mm×520 mm，共有4种厚度，分别为12、8、6、4 mm。根据安排进行组合，设置双层靶板的垂直间距为360 mm。

图 3 弹体材料的准静态拉伸曲线

Figure 3. Quasi-static tensile curves of projectile material

下载: 全尺寸图片幻灯片

试验平台如图4所示，采用口径40 mm的一级轻气炮作为发射装置，在弹体出炮管后以气动方式实现弹托分离。为了模拟弹体斜侵彻工况，设置靶板的倾斜角度为30°。利用高速摄像机对弹体的侵彻过程进行记录，得到弹体的姿态和速度等参数。此外，在靶板后方布置沙箱，对弹体进行软回收。

图 4 试验系统示意图

Figure 4. Schematic diagram of experimental system

下载: 全尺寸图片幻灯片

1.2 试验结果

共进行了8发斜侵彻试验，其中椭圆截面和非对称椭圆截面弹体各2发、圆形截面弹体4发作为对照参考。图5给出了不同弹体侵彻双层钢板的典型时刻照片。

图 5 弹体侵彻过程的典型时刻

Figure 5. Typical moments for different projectiles during penetration process

下载: 全尺寸图片幻灯片

图6给出了侵彻过程中的弹体姿态的定义。如图6(a)所示，将弹体的姿态角θ定义为弹体轴线与水平方向的夹角，弹体头部方向在水平线以上时姿态角为正；攻角α为弹体轴线与速度的夹角，弹体头部方向在速度方向以上时攻角为正；着角β为靶板法线与弹体速度的夹角，当靶面向上时为正。在试验过程中，弹体初始速度保持水平，则弹体的初始姿态角θ与攻角α相等，初始着角β=30°。如图6(b)所示，试验过程中，椭圆和非对称椭圆截面弹体的长轴位于水平方向。

图 6 斜侵彻中弹体不同角度参数及截面布置示意图

Figure 6. Diagram for the altitude angles and cross section of the projectile in oblique penetration

下载: 全尺寸图片幻灯片

根据高速摄像机的图片结果，得到弹体速度、姿态角、长度等试验数据如表2所示。其中，v₀、θ₀、l₀为靶前数据，v₁、θ₁、l₁为靶后数据，弹体穿靶前后的剩余长度之比δ = l₁/l₀。

表 2 弹体侵彻不同靶板的试验结果

Table 2. Penetration experimental results

试验编号	靶板编号	靶板厚度/mm	速度			姿态角			长度
试验编号	靶板编号	靶板厚度/mm	v₀/(m·s⁻¹)	v₁/(m·s⁻¹)	ΔE/J	θ₀/(°)	θ₁/(°)	Δθ/(°)	l₀/mm	l₁/mm	δ/%
CC-1	1	8	474	425	11012	−2.64	−7.56	−4.92	180	145	81
CC-1	2	6	425	388	7520	−4.95	−9.88	−4.93	145	145	100
CC-2	1	4	445	413	6864	16.17	15.65	−0.52	180	136	75
CC-2	2	4	413	387	5200	13.68	11.31	−2.37	136	136	100
CC-3	1	8	608	544	18432	14.68	20.81	6.13	180	133	74
CC-3	2	8	544	453	22681	19.95	27.47	7.52	133	133	100
CC-4	1	12	499	409	20430	3.15	0	−3.15	180	180	100
CC-4	2	8	409	348	11544	−2.23	−13.95	−11.72	180	109	61
EC-1	1	8	486	442	10208	−2.12	−6.37	−4.25	180	143	79
EC-1	2	6	442	404	8037	−3.23	−4.85	−1.62	143	143	100
EC-2	1	8	493	453	9460	0	−2.23	−2.23	180	180	100
EC-2	2	6	453	416	8038	−3.28	−8.94	−5.66	180	138	77
AC-1	1	8	489	435	12474	−1.60	−8.65	−7.05	180	128	72
AC-1	2	6	435	389	9476	−5.44	−10.75	−5.31	128	128	100
AC-2	1	8	473	429	9922	−2.77	−9.54	−6.77	180	180	100
AC-2	2	6	429	395	7004	−7.24	−13.89	−6.65	180	133	74

下载: 导出CSV

| 显示表格

由表2可知，除CC-2、CC-3外，其余6发弹体以小攻角姿态着靶，在470～500 m/s入射速度下穿透双层靶板。从能量角度进行分析，忽略结构断裂时两部分弹体的速度差异，根据弹体总质量与靶前、靶后速度计算得到侵彻不同靶板的动能损失ΔE。由于第1层靶和第2层靶的厚度不同，弹体穿过首层靶板的动能损失更大。

由图5可知，在初始攻角较小（−3°～3°）时，穿靶过程中弹体呈逆时针方向偏转，出现弹体“低头”现象；随着侵彻层数的增加，弹体出靶后的姿态角逐渐增大，侵彻弹道向下偏转。对比靶板厚度为8 mm+6 mm的5组试验结果，可以看到非对称椭圆截面弹体穿过首层靶板的姿态角变化量Δθ大于圆截面弹体和椭圆截面弹体，这与文献[17]的结论相似。

对比CC-1和CC-4的结果，在初始攻角较小时，靶板厚度的小幅度增加并未改变弹体姿态与弹道的偏转方向。而从CC-1和CC-3来看，在靶板厚度相当时，当攻角从−2.6°增加到14.7°时，弹体偏转方向发生改变，出现了弹体“抬头”现象，侵彻弹道从向下偏转变为向上偏转。此外，根据CC-2和CC-3的试验结果，在大攻角条件下，当靶板厚度从8 mm减小至4 mm时，仍保持弹体“低头”与弹道向下偏转规律。可见，存在一个使弹体姿态和弹道轨迹发生改变的临界攻角，且该临界值随着靶板厚度的减小而增大。

由试验过程的高速摄影照片发现，8发试验中弹体均出现局部断裂，其中CC-1、CC-2、CC-3、EC-1、AC-1弹体在穿过首层靶板时出现断裂，而CC-4、EC-2和AC-2实验中弹体在穿过第2层靶板时出现断裂。由图5可知，弹体断裂主要发生在尾部出靶阶段，这是由于穿靶过程中弹体姿态偏转，导致出靶时尾部撞击靶板，引起弹身断裂。图7给出了5发入靶条件相近的弹体断裂情况，可以看到弹体断裂位置未产生明显的变形，主要破坏模式为脆性断裂。

图 7 弹体的破坏情况

Figure 7. Damages of projectiles

下载: 全尺寸图片幻灯片

根据表2，以上5发弹体断裂后的剩余长度约为初始长度的72%～81%。由图5可知，弹体因侵彻首层靶板而断裂后，未出现再次断裂现象。从弹体的结构强度来看，断裂导致弹体长径比减小，本身承载能力相对提高，而且穿靶后速度有所下降，弹体所受载荷减小，因此试验过程中弹体的结构破坏主要为“一次断裂”。由5发弹体的测量结果来看，圆截面弹体的剩余长度大于椭圆截面弹体，而非对称椭圆截面弹体最小。

2. 弹体斜侵彻双层钢板的仿真研究

2.1 数值仿真模型

为深入分析弹体侵彻过程中的动态载荷及结构响应细节，采用有限元程序Abaqus/Explicit对该过程进行数值仿真。基于试验中的弹体和靶板结构，采用C3D8R单元建立模型，弹体单元尺寸为1 mm，靶板中心区域单元尺寸为1 mm，外围网格尺寸为5 mm。为简化计算，采用对称模型，弹靶接触定义为通用接触，并对靶板的边界施加固定约束。

弹体和靶板材料均选用Johnson-Cook强度模型，其中弹体采用拉伸断裂应变ε_T作为损伤判据。结合材料准静态力学性能试验结果，将其设为0.05，靶板参数根据文献确定，具体参数如表3～4所示。表中，ρ为密度，E为弹性模量，v为泊松比，T_r为参考温度，T_m为熔点温度， $\dot\varepsilon$ ₀为参考应变率，A为准静态屈服强度，B、n为应变硬化常数，C为应变率常数，m为温度软化系数，ε_T为拉伸断裂应变；D₁、D₂、D₃、D₄、D₅为Johnson-Cook损伤参数。

表 3 弹体30CrMnSiNi2A材料参数^[31]

Table 3. Material parameters of 30CrMnSiNi2A^[31]

ρ/(g·cm⁻³)	E/GPa	v	T_r/K	T_m/K	A/MPa	B/MPa	n	m	C	$\dot\varepsilon$ ₀/s⁻¹	ε_T
7.85	210	0.3	294	1760	1600	810	0.479	1	0.04	2.1×10⁻³	0.05

下载: 导出CSV

| 显示表格

表 4 靶板45钢材料参数^[32]

Table 4. Material parameters of 45 steel^[32]

ρ/(g·cm⁻³)	E/GPa	v	A/MPa	B/MPa	n	m	C
7.8	210	0.33	507	320	0.28	1.06	0.064

T_r/K	T_m/K	$\dot\varepsilon$ ₀/s⁻¹	D₁	D₂	D₃	D₄	D₅
294	1760	1	0.1	0.76	1.57	0.005	−0.84

下载: 导出CSV

| 显示表格

2.2 试验有效性验证

根据试验结果，现有试验平台对500 g弹体的稳定发射速度约为500 m/s。考虑到与高超声速武器的末端速度存在差异，利用数值仿真方法，对不同速度下圆形弹体斜侵彻钢板的动态载荷进行分析，从而验证本文试验的有效性。仿真工况设置靶板倾角为30°，厚度为8 mm，弹体攻角为2°，速度为500、1000 m/s，两种速度下弹体的轴向和横向载荷结果如图8所示。

图 8 不同入射速度下弹体的载荷时程曲线

Figure 8. Time history curves of projectile load under different impact velocities

下载: 全尺寸图片幻灯片

结果表明，速度变化的影响主要体现在弹体头部穿靶阶段，对弹身穿靶阶段的载荷影响较小。在头部穿靶阶段，1000 m/s速度下弹体的载荷约为500 m/s速度下弹体载荷的2倍；在此过程中，弹体可能出现头部侵蚀、轴向屈曲等结构响应行为，可通过提高弹体材料强度和壳体厚度等方式加以避免。对于弹身穿靶阶段，两种速度下弹体载荷差异较小；由于横向载荷的作用，会出现弹体弯曲/断裂的结构响应行为。此外，高速条件下弹体的穿靶时间更短，姿态偏转幅度变小，这有利于降低弹身受到的横向载荷。因此，从以上分析来看，速度变化对弹体弯曲/断裂的结构响应行为影响较小，针对速度为500 m/s的弹体斜侵彻试验进行分析，对弹体弯曲/断裂的结构响应研究仍具有参考价值。

2.3 数值仿真和试验结果对比分析

现实中弹体以大攻角侵彻的现象较少，在此仅分析小攻角的情况。图9给出了速度为480 m/s、攻角为−2°、靶板着角为30°条件下三种弹体的弹道轨迹仿真结果，并与CC-1、EC-1、AC-1的试验结果进行了对比。由图9可知，数值仿真得到的三种弹体向下偏转的弹道轨迹以及弹体姿态与试验结果基本符合，穿靶后的弹体因姿态偏转而处于“低头”状态。

图 9 弹体侵彻轨迹的对比

Figure 9. Comparison of simulated and experimental results on penetration trajectories

下载: 全尺寸图片幻灯片

图10～13给出了三种弹体试验与数值仿真结果的对比。图10为弹体速度的对比结果，在一定误差范围内，椭圆截面弹体的剩余速度略高于其他两种弹体。图11给出了由数值仿真得到的不同截面弹体动能结果，可以看到，从侵彻首层靶板来看，在相同的入射条件下，椭圆弹体的动能损失比其他两种弹体小约13%。图12为弹体姿态角的对比结果，除EC-1弹体的试验结果以外，整体结果较为接近，同时表现出弹体在穿靶以及靶间飞行过程中的姿态变化。从失效情况来看，穿过首层靶板后，三种弹体尾部均发生断裂，且由于弹体所受载荷位置随弹体侵彻过程而不断变化，整体表现为“斜向”断裂模式，与试验结果吻合较好。图13为弹体剩余长度的对比结果，可以看到，数值仿真结果中弹体剩余长度比例约为74% ～ 78%，与试验结果较为接近。

图 10 弹体速度对比

Figure 10. Comparison of projectile velocities

下载: 全尺寸图片幻灯片

图 11 不同弹体的动能对比

Figure 11. Comparison of kinetic energies of different projectiles

下载: 全尺寸图片幻灯片

图 12 弹体姿态角对比

Figure 12. Comparison of projectile attitude angles

下载: 全尺寸图片幻灯片

图 13 弹体剩余长度对比

Figure 13. Comparison of projectile residual lengths

下载: 全尺寸图片幻灯片

基于上述分析结果，可见本文采用的数值仿真方法具有较好的可靠性，可用于弹体斜侵彻双层钢板的结构响应和失效分析。

2.4 弹体斜侵彻双层钢板的动态载荷

为了建立弹体结构动力响应的理论模型，需要得到弹体所受的动态载荷。为此，图14给出了速度为480 m/s、攻角为−2°、靶板着角为30°条件下三种弹体受到轴向和横向载荷的数值仿真结果。

图 14 弹体载荷时程曲线

Figure 14. Time history curves of projectile load

下载: 全尺寸图片幻灯片

如图14(a)所示，弹体穿过每层靶板的过程中经历了2次轴向载荷峰值。第1次载荷峰值出现在头部贯穿阶段，第2次载荷峰值出现在弹身后半段穿靶阶段。弹体的第1次轴向载荷峰值约为第2次载荷峰值的2～2.5倍。结合弹体断裂的分析结果，可认为轴向载荷并不是导致弹体破坏的主要原因。

从弹靶相互作用来看，横向载荷是造成斜侵彻弹体姿态偏转和结构断裂的主要原因。如图14(b)所示，依据横向载荷的变化情况，可将弹体侵彻过程分为3个阶段，分别对应载荷峰值①、②、③，图15给出了典型工况的等效应力云图。在头部压入阶段，如图15(a)所示，由于着角的存在，弹体头部下侧首先与靶板接触，导致弹体头部承受向上的作用力。在此期间达到载荷正向峰值①，由于作用时间较短，弹体头部向上偏转幅度较小。随后是头部贯穿阶段，随着弹体头部继续运动，弹靶接触面积增大，靶板变形更大并形成花瓣状破坏。如图15(b)所示，上侧靶板向背面弯曲的变形程度大于下侧靶板，导致弹体受到的横向载荷合力向下，且在头部完全穿过靶板时达到载荷峰值②。由于作用力较大且作用时间较长，弹体姿态向下偏转幅度较大。此外，在头部穿靶过程中，由于载荷的作用，距离头部约L/3处局部应力最大。最后是弹身穿靶阶段，如图15(c)所示，在先前载荷的作用下，弹体处于低头状态，弹身持续撞击上侧靶板而受到向下的作用力，在弹体尾部穿靶时达到载荷峰值③；在侵彻过程中，弹体在力矩作用下产生顺时针角速度，使得弹体出靶后出现姿态修正效果。

图 15 弹体侵彻的不同阶段

Figure 15. Penetration stages of projectile

下载: 全尺寸图片幻灯片

结合数值仿真结果以及轴向、横向载荷历程曲线，对弹体斜侵彻首层钢板的失效模式进行分析。图16给出了典型工况的塑性应变云图，从图中可见塑性应变主要分布于弹体头部和弹身中后段。如图16(a)所示，在头部穿靶过程中，靶板在弹体的挤压作用下变形、破坏，此时弹体头部受到以轴向载荷为主的作用力，使得弹体头部存在较严重的塑性应变。随着弹体继续运动，在横向载荷的作用下弹体姿态发生偏转，此时弹身受到以横向载荷为主的作用力。随着载荷逐渐增大，弹身达到屈服状态并产生塑性应变，引起弹体局部损伤，如图16(b)所示。由于载荷作用位置的移动，损伤部位受到载荷的弯曲作用而斜向扩展，最终导致弹体断裂，如图16(c)所示。因此，可将弹体的结构破坏总结为由“局部受压”和“整体弯曲”共同作用而引起的弹体断裂。

图 16 弹体的失效模式

Figure 16. Failure modes of projectile

下载: 全尺寸图片幻灯片

结合数值仿真得到的动态载荷结果，对三种弹体侵彻过程中的动态响应进行对比。对于首层靶板，相同入射条件下，三种弹体的载荷形式基本相同，载荷峰值的差异在9%以内。图17给出了三种弹体侵彻首层靶板过程中的姿态角变化过程，可见随着载荷作用位置的变化，弹体姿态角呈先增大后减小趋势，而且在三种弹体中，非对称椭圆截面弹体姿态偏转幅度最大。从弹体结构响应来看，当截面积相同时，由于截面形状的差异，三种弹体的截面屈服弯矩不同，从大到小依次为圆形、椭圆、非对称椭圆。说明在载荷相当的条件下，非对称椭圆截面弹体更容易达到屈服状态，断裂位置更加靠近头部。

图 17 三种弹体的姿态角对比

Figure 17. Comparison of attitude angle among three different projectiles

下载: 全尺寸图片幻灯片

3. 弹体斜侵彻双层钢板的结构响应模型

3.1 弹体响应的自由梁模型

由于高速侵彻过程中，弹体的边界条件与空间自由梁撞击过程一致，因此可借助自由梁模型对其动力响应进行分析。将弹体视为质量分布不均匀的自由梁，分为头部实心段和弹身空心段。由于实心段不存在结构强度问题，且相对较短，因此简化模型中不考虑截面形状的变化。如图18所示，设弹体总长度为L，其中实心段长为L₁，线密度为ρ₁，空心段线密度为ρ₂；横向载荷为F，作用位置距离左端部为e (e>L₁)，弹体左端的平动加速度为 $\ddot u$ ，相对质心的转动加速度为 $\ddot \theta$ 。

图 18 横向载荷作用下的自由梁模型

Figure 18. Free-free beam model for the projectile under lateral load

下载: 全尺寸图片幻灯片

对简化模型的主要参数无量纲化，可得自由梁线密度之比μ = ρ₂ / ρ₁，实心段长度λ = L₁ / L，载荷作用位置ξ = e / L，任意位置 $\bar x=x/L$ 。

设弹体的总质量M_t为：

$\frac{{{M_{\text{t}}}}}{{{\rho _1}L}} = \lambda (1 - \mu ) + \mu$

(2)

弹体质心的位置X_c为：

$\frac{{{X_{\text{c}}}}}{L} = \frac{1}{2}\frac{{{\lambda ^2}(1 - \mu ) + \mu }}{{\lambda (1 - \mu ) + \mu }}$

(3)

弹体对质心的转动惯量J_c为：

$\dfrac{{{J_{\text{c}}}}}{{{\rho _1}{L^3}}} = \dfrac{{{\lambda ^4}(1 - \mu ) + 4{\lambda ^3}\mu - 6{\lambda ^2}\mu + 4\lambda \mu + \dfrac{{{\mu ^2}}}{{1 - \mu }}}}{{12\left(\lambda + \dfrac{\mu }{{1 - \mu }}\right)}}$

(4)

由刚体动力学可得弹体的运动方程：

$\left\{ {\begin{array}{*{20}{c}} {{M_{\text{t}}}(\ddot u - \ddot \theta {X_{\text{c}}}) = F} \\ {{J_{\text{c}}}\ddot \theta = F({X_{\text{c}}} - e)} \end{array}} \right.$

(5)

从而解得弹体头部加速度和转动加速度为：

$\dfrac{{\ddot u}}{{F/({\rho _1}L)}} = \dfrac{{4{\lambda ^3} - 6{\lambda ^2}\xi + \dfrac{{2\mu }}{{1 - \mu }}(2 - 3\xi )}}{{{\lambda ^4}(1 - \mu ) + 4{\lambda ^3}\mu - 6{\lambda ^2}\mu + 4\lambda \mu + \dfrac{{{\mu ^2}}}{{1 - \mu }}}}$

(6)

$\dfrac{{\ddot \theta }}{{F/({\rho _1}{L^2})}} = \dfrac{{6{\lambda ^2} - 12\lambda \xi + \dfrac{{6\mu }}{{1 - \mu }}(1 - 2\xi )}}{{{\lambda ^4}(1 - \mu ) + 4{\lambda ^3}\mu - 6{\lambda ^2}\mu + 4\lambda \mu + \dfrac{{{\mu ^2}}}{{1 - \mu }}}}$

(7)

基于运动方程的求解和自由梁的受力分析，可得自由梁任意位置的无量纲剪力：

$\begin{split} &\dfrac{{Q(\bar x)}}{F} =\\ &\left\{ {\begin{array}{*{20}{l}} \dfrac{{3{{\bar x}^2}\left[ {{\lambda ^2} - 2\lambda \xi + \dfrac{\mu }{{1 - \mu }}(1 - 2\xi )} \right] - 2\bar x\left[ {2{\lambda ^3} - 3{\lambda ^2}\xi + \dfrac{\mu }{{1 - \mu }}(2 - 3\xi )} \right]}}{{{\lambda ^4}(1 - \mu ) + 4{\lambda ^3}\mu - 6{\lambda ^2}\mu + 4\lambda \mu + \dfrac{{{\mu ^2}}}{{1 - \mu }}}}& 0 {\text{≤}} \bar x {\text{＜}} \lambda \\ \dfrac{{{{\bar x}^2}\left[ {3{\lambda ^2} - 6\lambda \xi + \dfrac{{3\mu }}{{1 - \mu }}(1 - 2\xi )} \right] - \bar x\left[ {4{\lambda ^3} - 6{\lambda ^2}\xi + \dfrac{{2\mu }}{{1 - \mu }}(2 - 3\xi )} \right] - \left[ {4{\lambda ^3} - 3{\lambda ^2} + 6\lambda \xi (1 - \lambda ) + \dfrac{\mu }{{1 - \mu }}} \right]}}{{{\lambda ^4}\left(\dfrac{{1 - \mu }}{\mu }\right) + 4{\lambda ^3} - 6{\lambda ^2} + 4\lambda + \dfrac{\mu }{{1 - \mu }}}} - 1 & \lambda {\text{≤}} \bar x {\text{＜}} \xi \\ \dfrac{{(\bar x - 1)\left\{ {3\bar x\left[ {{\lambda ^2} - 2\lambda \xi + \dfrac{\mu }{{1 - \mu }}(1 - 2\xi )} \right] - \left[ {4{\lambda ^3} - 3{\lambda ^2} + \dfrac{\mu }{{1 - \mu }} + 6\lambda \xi (1 - \lambda )} \right]} \right\}}}{{{\lambda ^4}\left(\dfrac{{1 - \mu }}{\mu }\right) + 4{\lambda ^3} - 6{\lambda ^2} + 4\lambda + \dfrac{\mu }{{1 - \mu }}}} & \xi {\text{≤}} \bar x {\text{≤}} 1 \end{array}} \right.\end{split}$

(8)

对剪力积分，可得弹体任意位置的无量纲弯矩：

$\frac{{M(\bar x)}}{{FL}} = \int_0^{\bar x} {\frac{{Q(\bar x)}}{F}} {\rm{d}}\bar x$

(9)

皮爱国等^[24]考虑到弹体侵彻过程中轴向、横向载荷的耦合作用，提出弹体结构的塑性屈服准则。考虑到实验中弹体以脆性断裂为主，本文主要针对弹体结构的屈服强度进行分析。以材料弹性屈服强度σ_y作为屈服函数的条件，其表达式为：

$\psi (x) = \frac{{\left| {N(x)} \right|}}{{{N_{\text{y}}}(x)}} + \frac{{\left| {M(x)} \right|}}{{{M_{\text{y}}}(x)}} - 1$

(10)

弹体的轴向力N(x)可由下式求得：

$N(x) = \left\{ {\begin{split} & {N\dfrac{{m(x)}}{{m(L)}}}&{0 {\text{≤}} x {\text{≤}} e} \\ & {N(1 - \dfrac{{m(x)}}{{m(L)}})}&{e {\text{＜}} x {\text{≤}} L} \end{split}} \right.$

(11)

式中：N为弹靶作用位置处的轴向载荷，任意位置处的轴向载荷以弹靶作用位置e为分界，头部一侧处于受拉状态，尾部一侧处于受压状态，两侧载荷呈线性分布；m(x)为端部到任意位置的弹体质量。N_y(x)为轴向屈服载荷，M_y(x)为横向屈服弯矩，表达式为：

${N_{\text{y}}}(x) = {\sigma _{\text{y}}}S(x)$

(12)

${M_{\text{y}}}(x) = \frac{{{\sigma _{\text{y}}}I(x)}}{{{y_{\max }}(x)}}$

(13)

式中：σ_y为材料屈服强度，S(x)为任意位置的截面面积，I(x)为截面惯性矩，y_max(x)为截面各点距离中心轴的最大距离。考虑到弹体的侵彻姿态，本文以短轴方向上的截面抗弯能力为参考。

3.2 弹体结构响应的分析

对于本文中的弹体，实心段长度λ = 0.25，线密度之比为：

$\mu = \frac{{{\rho _2}}}{{{\rho _1}}} = 1 - {(1 - h/r)^2}$

(14)

式中：h为弹体壁厚，取为4 mm，r为与椭圆弹体截面积相等的圆形弹体截面半径，取为15 mm，则弹体结构简化模型中线密度之比μ = 0.46。代入式(9)～(10)，可得弹体不同位置受载时的无量纲弯矩分布，结果如图19所示。可见，随着载荷位置的移动，弹体无量纲弯矩分布随之变化。根据理论分析结果，当载荷作用位置位于λ<ξ<0.82时，弯矩最大位置即为载荷作用位置。结合轴向载荷的分布规律，由式(12)可知，轴力最大位置同样为载荷作用位置。因此，可通过校核载荷作用位置处的屈服函数，判断弹体受载时是否达到屈服状态，进而评估侵彻过程中某时刻的弹体结构失效情况。

图 19 无量纲弯矩分布

Figure 19. Distribution of dimensionless bending moment

下载: 全尺寸图片幻灯片

根据图14的轴向和横向载荷仿真结果，对三种截面弹体载荷作用位置的屈服函数进行计算，结果如图20所示。可以看到，由于弯矩在屈服函数中占比更大，因此弹身不同位置受载时的屈服函数分布与横向载荷曲线相似，出现2次峰值，其中第2次峰值对应位置已达到屈服条件。考虑到弹体材料具有应变率强化效应，且断裂前会产生塑性变形，因此弹体的实际断裂位置滞后于理论计算的屈服位置。表5给出了弹体剩余长度比例的试验、数值仿真以及理论模型结果的对比，可以看出理论模型结果与试验结果吻合较好，相对误差在8%以内。

图 20 移动载荷作用下弹体的屈服函数

Figure 20. Yield function of projectile under moving load

下载: 全尺寸图片幻灯片

表 5 弹体剩余长度的不同结果对比

Table 5. Comparison of results on projectile residual length

弹型	试验结果	数值仿真结果	理论模型结果	理论模型相对误差/%
圆形	0.81	0.78	0.75	7.41
椭圆形	0.78	0.75	0.73	6.42
非对称椭圆	0.72	0.74	0.71	1.38

下载: 导出CSV

| 显示表格

4. 结　论

开展了圆形、椭圆和非对称椭圆三种截面弹体斜侵彻双层钢板的试验研究，获得了不同弹体的弹道轨迹以及结构失效情况。在此基础上，针对弹体斜侵彻工况开展数值仿真工作，并结合试验结果，对弹靶作用过程以及弹体动态载荷进行分析。最终，基于空间自由梁理论，建立了弹体结构响应分析模型，得到载荷作用下弹体的剪力、弯矩以及屈服函数分布规律，并给出弹体结构强度分析方法。主要结论如下。

(1)根据试验结果，当弹体以正着角水平侵彻多层钢板时，存在一个临界攻角，当攻角小于该值时，弹体侵彻过程中会发生低头，随着穿靶层数的增加，姿态偏转幅度逐渐增大，弹道轨迹向下偏转。当攻角大于临界攻角时，弹体由低头转变为抬头，弹道轨迹向上偏转。该临界攻角随着靶板厚度的减小而增大。

(2)由于试验中，弹体材料韧性较差，弹体的失效模式为脆性断裂，且断裂位置为弹体中后段，距离头部0.72L～0.81L，其中非对称椭圆截面弹体的断裂位置最靠近头部。对双层钢板而言，在弹体侵彻首层靶板出现断裂后，随着弹体长度减小以及速度下降，后续穿靶过程中未出现再次断裂。

(3)通过数值仿真得到弹体的轴向和横向载荷，结果表明在斜侵彻过程中弹体会经历2次横向冲击载荷，第2次载荷峰值更大，是造成弹体结构断裂的主要原因。此外，根据三种弹体的数值仿真结果，在相同入射条件下，非对称椭圆截面弹体的载荷峰值及姿态偏转幅度较大。

(4)利用自由梁模型，得到了弹体任意位置受到横向载荷作用下的剪力和弯矩分布规律。结合弹体载荷的数值仿真结果，进一步建立了考虑弯矩和轴力耦合作用的弹体结构强度和失效分析方法。通过与数值仿真及实验结果的对比，该方法能够准确预测斜侵彻过程中弹体的断裂行为。

图 1 头部模型

Figure 1. Head models

下载: 全尺寸图片幻灯片

图 2 模型与实验示意图

Figure 2. Model and experimental diagram

下载: 全尺寸图片幻灯片

图 3 不同爆心距处压力随时间的变化

Figure 3. Variation of pressure with time at different distances away from explosion center

下载: 全尺寸图片幻灯片

图 4 不同强度冲击波正面冲击下额部颅内压对比

Figure 4. Comparison of frontal intracranial pressures under forward impact by different-strength shock waves

下载: 全尺寸图片幻灯片

图 5 不同强度冲击波正面冲击下脑室颅内压对比

Figure 5. Comparison of ventricle intracranial pressures under forward impact by different-strength shock waves

下载: 全尺寸图片幻灯片

图 6 不同强度冲击波正面冲击下顶部颅内压对比

Figure 6. Comparison of parietal intracranial pressures under forward impact by different-strength shock waves

下载: 全尺寸图片幻灯片

图 7 不同强度冲击波正面冲击下枕部颅内压对比

Figure 7. Comparison of occipital intracranial pressure under forward impact by different-strength shock waves

下载: 全尺寸图片幻灯片

图 8 颅骨最大应力云图

Figure 8. Maximum stress cloud map of skull

下载: 全尺寸图片幻灯片

图 9 不同强度不同方向冲击波冲击下的颅内压峰值对比

Figure 9. Comparison of intracranial pressure peak values under different intensity and direction impact wave shocks.

下载: 全尺寸图片幻灯片

表 1 颅骨组织材料参数^[15-16]

Table 1. Skull tissue material parameters^[15-16]

结构	ρ/(g∙cm⁻³)	E/GPa	ν	σ_y/MPa	$G$ /GPa	n	Cowper-Symonds模型		ε_p/%	来源
结构	ρ/(g∙cm⁻³)	E/GPa	ν	σ_y/MPa	$G$ /GPa	n	C	P	ε_p/%	来源
顶骨	2	11.5	0.3	90	1.15	0.1	2.5	7.0	0.02	文献[16]
颞骨	2	11.5	0.3	90	1.15	0.1	2.5	7.0	0.02	文献[16]
颧骨	2	11.5	0.3	90	1.15	0.1	2.5	7.0	0.02	文献[16]
枕骨	2	11.5	0.3	90	1.15	0.1	2.5	7.0	0.02	文献[16]
额骨	2	11.5	0.3	90	1.15	0.1	2.5	7.0	0.02	文献[16]
蝶骨	2	11.5	0.3	90	1.15	0.1	2.5	7.0	0.02	文献[16]
面骨	1.71	5.37	0.19	−	−	−	−	−	−	文献[15]
下颌骨	2	11.5	0.3	145	1.15	0.1	2.5	7.0	−	文献[16]

下载: 导出CSV

表 2 脑组织材料参数^{[15, 18-19]}

Table 2. Material parameters of brain tissue^{[15, 18-19]}

组织	ρ/(g∙cm⁻³)	E/MPa	ν	${G}_{0}$ /kPa	${G}_{\infty }$ /kPa	$\beta$ /s⁻¹	K/MPa	来源
头皮	1.2	16.7	0.42	−	−	−	−	文献[18]
脑脊液	1.04	−	−	100	20	100	1050	文献[19]
硬脑膜	1.13	31.5	0.45	−	−	−	−	文献[15]
软脑膜	1.13	31.5	0.45	−	−	−	−	文献[15]
大脑	1.06			1.66	0.928	16.95	557	文献[18]
小脑	1.06			1.16	0.928	16.95	557	文献[18]
脑干	1.04			1.66	0.928	16.95	557	文献[18]

下载: 导出CSV

表 3 不同爆心距处超压峰值的对比

Table 3. Comparison of overpressure peaks at different distances away from explosion center

d/m	Δp₊/kPa		相对误差/%
d/m	数值模拟	经验公式	相对误差/%
0.75	307	298.7	2.78
0.80	258	256.7	0.51
0.85	223	223.4	0.18
0.90	202	196.2	2.96
0.95	190	173.8	9.32

下载: 导出CSV

表 4 不同强度冲击波正面冲击下颅内压峰值对比

Table 4. Comparison of peak intracranial pressures under forward impact by different-strength shock waves

冲击波强度/kPa	颅内压峰值/kPa				备注
冲击波强度/kPa	额部	脑室	顶部	枕部	备注
75	102	32		−30	实验cad4^[17]
	162	46	53		实验cad5^[17]
	132	39			cad4与cad5的平均值^[17]
	103	55	55	−36	仿真^[17]
	155	52	62	−59	数值模拟
	4.51%	5.77%	11.29%	38.98%	最小误差
102	142	47		−40	实验cad4^[17]
	220	63	83		实验cad5^[17]
	181	55			cad4与cad5的平均值^[17]
	168	91	95	−55	仿真^[17]
	292	84	104	−112	数值模拟
	24.66%	8.33%	8.65%	50.89%	最小误差

下载: 导出CSV

表 5 不同强度冲击波侧面冲击下颅内压峰值对比

Table 5. Comparison of peak intracranial pressure under lateral impact by different-strength shock waves

冲击波强度/kPa	颅内压峰值/kPa				备注
冲击波强度/kPa	额部	脑室	顶部	枕部	备注
75	20	20		13	实验cad4^[17]
	38	39	54		实验cad5^[17]
	29	30			cad4与cad5的平均值^[17]
	69	42	49	44	仿真^[17]
	76	40	72	68	数值模拟
	9.21%	1.25%	25%	35.29%	最小误差
102		24		37	实验cad4^[17]
	94	45	62		实验cad5^[17]
		34.5			cad4与cad5的平均值^[17]
	94	50	55	70	仿真^[17]
	118	56	88	68	数值模拟
	20.34%	10.71%	37.5%	2.94%	最小误差

下载: 导出CSV

表 6 不同强度冲击波背面冲击下颅内压峰值对比

Table 6. Comparison of peak intracranial pressures under back impact by different-strength shock waves

冲击波强度/kPa	颅内压峰值/kPa				备注
冲击波强度/kPa	额部	脑室	顶部	枕部	备注
75	−55	19		53	实验cad4^[17]
	−133	19	88		实验cad5^[17]
	−94	19			cad4与cad5的平均值^[17]
	−62	40	53	112	仿真^[17]
	−106	42	97	104	数值模拟
	11.32%	4.76%	9.28%	7.69%	最小误差
102	−74	22		90	实验cad4^[17]
	−156	25	150		实验cad5^[17]
	−115	23.5			cad4与cad5的平均值^[17]
	−73	49	127	138	仿真^[17]
	−135	45	143	149	数值模拟
	14.81%	8.89%	4.89%	7.38%	最小误差

下载: 导出CSV

参考文献(24)

[1]	CONNELLY C, MARTIN K, ELTERMAN J, et al. Early traumatic brain injury screen in 6594 inpatient combat casualties [J]. Injury, 2017, 48(1): 64–69. DOI: 10.1016/j.injury.2016.08.025.
[2]	GALARNEAU M R, WOODRUFF S I, DYE J L, et al. Traumatic brain injury during Operation Iraqi Freedom: findings from the United States Navy-Marine Corps Combat Trauma Registry [J]. Journal of Neurosurgery, 2008, 108(5): 950–957. DOI: 10.3171/JNS/2008/108/5/0950.
[3]	BHATTACHARJEE Y. Shell shock revisited: solving the puzzle of blast trauma [J]. Science, 2008, 319(5862): 406–408. DOI: 10.1126/science.319.5862.406.
[4]	YU X C, AZOR A, SHARP D J, et al. Mechanisms of tensile failure of cerebrospinal fluid in blast traumatic brain injury [J]. Extreme Mechanics Letters, 2020, 38: 100739. DOI: 10.1016/j.eml.2020.100739.
[5]	KULKARNI S G, GAO X L, HORNER S E, et al. Ballistic helmets: their design, materials, and performance against traumatic brain injury [J]. Composite Structures, 2013, 101: 313–331. DOI: 10.1016/j.compstruct.2013.02.014.
[6]	MOSS W C, KING M J, BLACKMAN E G. Skull flexure from blast waves: a mechanism for brain injury with implications for helmet design [J]. Physical Review Letters, 2009, 103(10): 108702. DOI: 10.1103/PhysRevLett.103.108702.
[7]	GOELLER J, WARDLAW A, TREICHLER D, et al. Investigation of cavitation as a possible damage mechanism in blast-induced traumatic brain injury [J]. Journal of Neurotrauma, 2012, 29(10): 1970–1981. DOI: 10.1089/neu.2011.2224.
[8]	康越, 张仕忠, 张远平, 等. 基于激波管评价的单兵头面部装备冲击波防护性能研究 [J]. 爆炸与冲击, 2021, 41(8): 085901. DOI: 10.11883/bzycj-2020-0395. KANG Y, ZHANG S Z, ZHANG Y P, et al. Research on anti-shockwave performance of the protective equipment for the head of a soldier based on shock tube evaluation [J]. Explosion and Shock Waves, 2021, 41(8): 085901. DOI: 10.11883/bzycj-2020-0395.
[9]	MAO H J, ZHANG L Y, JIANG B H, et al. Development of a finite element human head model partially validated with thirty five experimental cases [J]. Journal of Biomechanical Engineering, 2013, 135(11): 111002. DOI: 10.1115/1.4025101.
[10]	COTTON R T, PEARCE C W, YOUNG P G, et al. Development of a geometrically accurate and adaptable finite element head model for impact simulation: the Naval Research Laboratory-Simpleware Head Model [J]. Computer Methods in Biomechanics and Biomedical Engineering, 2016, 19(1): 101–113. DOI: 10.1080/10255842.2014.994118.
[11]	GHAJARI M, HELLYER P J, SHARP D J. Computational modelling of traumatic brain injury predicts the location of chronic traumatic encephalopathy pathology [J]. Brain, 2017, 140(2): 333–343. DOI: 10.1093/brain/aww317.
[12]	CARMO G P, DYMEK M, PTAK M, et al. Development, validation and a case study: the female finite element head model (FeFEHM) [J]. Computer Methods and Programs in Biomedicine, 2023, 231: 107430. DOI: 10.1016/j.cmpb.2023.107430.
[13]	聂伟晓, 温垚珂, 董方栋, 等. 破片侵彻戴防弹头盔头部靶标钝击效应数值模拟 [J]. 兵工学报, 2022, 43(9): 2075–2085. DOI: 10.12382/bgxb.2022.0428. NIE W X, WEN Y K, DONG F D, et al. Numerical simulation of bludgeoning effect of fragments penetrating head target wearing bulletproof helmet [J]. Acta Armamentarii, 2022, 43(9): 2075–2085. DOI: 10.12382/bgxb.2022.0428.
[14]	栗志杰, 由小川, 柳占立, 等. 基于三维头部数值模型的颅脑碰撞损伤机理研究 [J]. 工程力学, 2019, 36(5): 246-56. DOI: 10.6052/j.issn.1000-4750.2018.04.0254. LI Z J, YOU X C, LIU Z L, et al. Study on the mechanism of brain injury during head impact based on the three-dimensional numerical head model [J]. Engineering Mechanics, 2019, 36(5): 246-256. DOI: 10.6052/j.issn.1000-4750.2018.04.0254.
[15]	张文超, 王舒, 梁增友, 等. 爆炸冲击波致颅脑冲击伤数值模拟研究 [J]. 北京理工大学学报, 2022, 42(9): 881-90. DOI: 10.15918/j.tbit1001-0645.2021.191. ZHANG W C, WANG S, LIANG Z Y, et al. Numerical simulation on traumatic brain injury induced by blast waves [J]. Transactions of Beijing institute of Technology, 2022, 42(9): 881-890. DOI: 10.15918/j.tbit1001-0645.2021.191.
[16]	毛征宇, 李泽民, 牛文鑫, 等. 不同载荷作用下头部生物力学响应仿真分析 [J]. 医用生物力学, 2016, 31(6): 532–539,547. DOI: 10.3871/j.1004-7220.2016.06.532. MAO Z Y, LI Z M, NIU W X, et al. The simulation analysis on biomechanical responses of human head under different loading conditions [J]. Journal of Medical Biomechanics, 2016, 31(6): 532–539,547. DOI: 10.3871/j.1004-7220.2016.06.532.
[17]	SHARMA S. Biomechanical analysis of blast induced traumatic brain injury: a finite element modeling and validation study of blast effects on human brain [M]. Detroit: Wayne State University, 2011.
[18]	李泽民. 子弹冲击防弹头盔动力学响应及防护性能仿真研究 [D]. 湘潭: 湖南科技大学, 2016. LI Z M. The bullet impact ballistic helmets simulation research on dynamic response and protective performance [D]. Xiangtan: Hunan University of Science and Technology, 2016.
[19]	PAVAN P G, NASIM M, BRASCO V, et al. Development of detailed finite element models for in silico analyses of brain impact dynamics [J]. Computer Methods and Programs in Biomedicine, 2022, 227: 107225. DOI: 10.1016/j.cmpb.2022.107225.
[20]	王智, 常利军, 黄星源, 等. 爆炸冲击波与破片联合作用下防弹衣复合结构防护效果的数值模拟 [J]. 爆炸与冲击, 2023, 43(6): 063202. DOI: 10.11883/bzycj-2022-0515. WANG Z, CHANG L J, HUANG X Y, et al. Simulation on the defending effect of composite structure of body armor under the combined action of blast wave and fragments [J]. Explosion and Shock Waves, 2023, 43(6): 063202. DOI: 10.11883/bzycj-2022-0515.
[21]	黄星源. 爆炸冲击波作用下颅脑损伤力学机制与头盔防护性能研究 [D]. 湘潭: 湖南科技大学, 2023. DOI: 10.27738/d.cnki.ghnkd.2023.000012. HUANG X Y. Research on mechanical mechanism of craniocerebral injury and protective performance of helmet under the blast wave [D]. Xiangtan: Hunan University of Science and Technology, 2023. DOI: 10.27738/d.cnki.ghnkd.2023.000012.
[22]	赵辉, 朱峰. 原发性颅脑冲击伤的生物力学机制 [J]. 创伤外科杂志, 2016, 18(6): 375–378. DOI: 10.3969/j.issn.1009-4237.2016.06.017. ZHAO H, ZHU F. The biomechanical mechanism of primary blast brain injury [J]. Journal of Traumatic Surgery, 2016, 18(6): 375–378. DOI: 10.3969/j.issn.1009-4237.2016.06.017.
[23]	康越, 马天, 黄献聪, 等. 颅脑爆炸伤数值模拟研究进展: 建模、力学机制及防护 [J]. 爆炸与冲击, 2023, 43(6): 061101. DOI: 10.11883/bzycj-2022-0521. KANG Y, MA T, HUANG X C, et al. Advances in numerical simulation of blast-induced traumatic brain injury: modeling, mechanical mechanism and protection [J]. Explosion and Shock Waves, 2023, 43(6): 061101. DOI: 10.11883/bzycj-2022-0521.
[24]	SAUNDERS R N, TAN X G, QIDWAI S M, et al. Towards identification of correspondence rules to relate traumatic brain injury in different species [J]. Annals of Biomedical Engineering, 2019, 47(9): 2005–2018. DOI: 10.1007/s10439-018-02.

施引文献

期刊类型引用(4)

1.	韩明海，刘闯，李鹏程，刘子涵，张先锋. 弹体高速侵彻花岗岩靶体的结构响应特性. 爆炸与冲击. 2025(01): 104-124 . 本站查看
2.	薛颖杰，陈智刚，杨芮，崔晋，康彦姝，付建平. 椭圆截面弹体侵彻多层混凝土靶的数值仿真. 中北大学学报(自然科学版). 2025(01): 35-44 . 百度学术
3.	蒋腾，武海军，邓希旻，全鑫，董恒，黄风雷. 变截面弹体斜侵彻两层间隔钢靶弹道特性. 兵工学报. 2025(02): 192-207 . 百度学术
4.	董恒，黄风雷，武海军，邓希旻，李萌，刘龙龙. 异形弹体高速侵彻/穿甲机理研究进展. 兵工学报. 2024(09): 2863-2887 . 百度学术