弹体侵彻混凝土的优化模型及可视化仿真研究

苏永超; 宁建国; 许香照

doi:10.11858/gywlxb.20240811

弹体侵彻混凝土的优化模型及可视化仿真研究

doi: 10.11858/gywlxb.20240811

北京理工大学爆炸科学与安全防护全国重点实验室, 北京　100081

基金项目: 国家自然科学基金（12372350）

详细信息

作者简介:
苏永超（1999－），男，硕士研究生，主要从事侵彻力学及可视化研究. E-mail：s15802440892@163.com

通讯作者:
许香照（1989－），男，博士，特别副研究员，主要从事冲击动力学研究. E-mail：xzxu@bit.edu.cn

中图分类号: O346.5; O521.9
计量
- 文章访问数: 326
- HTML全文浏览量: 81
- PDF下载量: 42
出版历程
- 收稿日期: 2024-05-13
- 修回日期: 2024-05-31
- 录用日期: 2024-09-18
- 网络出版日期: 2024-12-02
- 刊出日期: 2025-04-03

Optimization Model and Visualization Simulation of Projectile Penetration into Concrete

State Key Laboratory of Explosion Science and Safety Protection, Beijing Institute of Technology, Beijing 100081, China

摘要

摘要: 利用可视化仿真技术研究弹体侵彻混凝土的毁伤机理和靶板响应是爆炸冲击领域的重要课题。混凝土作为常见的建筑材料，在遭受爆炸冲击或高速弹体侵彻时，其毁伤行为复杂多变。介绍了一种理论研究与可视化技术相结合的可视化仿真方法。基于空腔膨胀理论建立了优化的侵彻计算模型，可以预测弹体侵彻混凝土的侵彻深度特征。利用可视化物理引擎，对弹体运动轨迹、开坑孔径、靶板损伤、碎石飞溅等进行了细致的表征处理和仿真，增强了场景的真实性和可靠性。开发的可视化仿真系统不仅能够从多角度观察弹体侵彻混凝土的过程，而且能够高效、准确地分析和预测弹体侵彻混凝土靶的损伤行为和动力响应，在建筑工程设计和安全评估中具有重要的应用前景，为理解和探索混凝土侵彻机理提供了新的视角。
- 弹体侵彻 /
- 混凝土 /
- 可视化仿真 /
- 物理引擎
Abstract: Using visual simulation technology to investigate the damage mechanism and target response of projectile penetration into concrete is an important research topic in the field of explosive impact. Concrete, as a common building material, has complex and varied damage behavior when subjected to explosive impact or high-speed projectile penetration. Herein, a visual simulation method is introduced, which is based on the combination of theoretical research and visualization technology. An optimized model of penetration calculation is established based on the theory of cavity expansion, which can predict the characteristics of the penetration depth of concrete penetrated by the projectile. Using a visualization physics engine, the trajectory of the projectile, the aperture of the open pit, the damage of the target slab, and the debris splash are carefully characterized and simulated, which enhances to the realism and reliability of the scene. The developed visual simulation system can not only observe the process of projectile penetration into concrete from multiple perspectives, but also efficiently and accurately analyze and predict the damage behavior and dynamic response of projectile penetration into concrete targets. It has important application prospects in the design and safety assessment of construction projects, providing a novel perspectives for understanding and exploring the mechanism of concrete penetration.
- projectile penetration /
- concrete /
- visual simulation /
- physics engine

HTML全文

镁铝合金（MB2）作为一类特殊的合金材料，具有低密度、高强度、易机械加工、耐腐蚀等特点，广泛应用于车辆工程、航空航天等领域，其动态加载下的力学和物理特性对相关结构设计等具有重要意义。国内外对MB2合金的早期研究主要集中于力学特性方面，通过开展低压动态响应特性实验研究，获得了材料的弹塑性响应^[1]、动态损伤^[2]及层裂行为^[3]的初步认识，近年来Millett等^[4]通过一维冲击加载研究了MB2合金在早期变形和位错条件下随加载应力、加载脉宽而变化的弹塑性和剪切强度行为。另外，对动态加载下MB2合金的物态方程及相变研究已开展了一系列工作，主要集中于冲击Hugoniot数据测量^[5]，然而与动态加载下物理、力学特性紧密相关的相变研究尚处于起步阶段。声速作为应力扰动在材料中传播的定量表征参量，是获知材料动态响应特性（如相变、屈服强度、剪切模量）的主要途径之一^[6–7]，然而相关的研究工作却鲜见报道。

本研究拟采用反向碰撞实验技术^[8–9]，结合具有高时空分辨率的全光纤激光干涉测速技术DPS（Doppler Pin System）^[10]，对MB2合金开展30～73 GPa压力范围内的冲击Hugoniot及声速测量实验，并与早期实验数据进行对比验证，分析MB2合金的冲击熔化行为。

1. 实　验

本研究涉及低压力区的声速测量，为此采用反向碰撞实验设计，即将样品材料制作成飞片，撞击LiF透明窗口，其原理如图1所示。在拉氏坐标下，飞片以速度W直接撞击窗口（t=t₁），在飞片和窗口中分别产生左行冲击波和右行冲击波，引起飞片/窗口界面粒子速度的突跃。当飞片中的左行冲击波到达后界面时，将反射中心稀疏波，该稀疏波的波速就是材料在冲击压缩下的声速。如果样品材料发生冲击熔化，则中心稀疏波将以单一塑性波的形式在样品内传播，并在飞片/窗口界面处（t=t₂）发生卸载，引起粒子速度的下降。如果样品材料没有发生冲击熔化，则中心稀疏波包括传播速度较快的弹性波和传播速度相对较慢的塑性波。当传播速度较快的弹性卸载波到达飞片/窗口界面时（t=t₂），界面粒子速度下降，在速度剖面上形成第1个拐点；当塑性卸载波到达飞片/窗口界面时（t=t₃），界面粒子速度再次突变，在速度剖面上形成第2个拐点；当后续稀疏波陆续到达飞片/窗口界面时，会导致界面粒子速度的连续下降；如果在卸载过程中样品材料发生相变，则也会在速度剖面上形成拐点。

图 1 反向碰撞实验示意图

Figure 1. Schematic of backward-impact experimental configuration

下载: 全尺寸图片幻灯片

根据图2所示的界面连续性条件（其中p为冲击压力，u为粒子速度），结合Rankine-Hugoniot关系^[11]，可以得到样品内的冲击波速度D_s为

图 2 反向碰撞实验p-u图

Figure 2. p-u relation for backward-impact experiment

下载: 全尺寸图片幻灯片

${D_{\rm{s}}} = \frac{{{\rho _{0{\rm{w}}}}{D_{\rm{w}}}{u_{\rm{w}}}}}{{{\rho _{0{\rm{s}}}}(W - {u_{\rm{w}}})}}\quad\quad\quad\quad\quad\quad\quad(1)\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad$

式中：D、 $\rho$ 、W分别为冲击波速度、密度和飞片速度，下标s和w分别对应样品和窗口。如果窗口材料的D-u曲线满足线性关系

${D_{\rm{w}}} = {C_{0{\rm{w}}}} + {\lambda _{\rm{w}}}{u_{\rm{w}}}\quad\quad\quad\quad\quad\quad\quad(2)\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad$

式中：C_0w和 ${\lambda _{\rm{w}}}$ 为窗口材料的Hugoniot参数，则（1）式可表示为

${D_{\rm{s}}} = \frac{{{\rho _{0{\rm{w}}}}({C_{0{\rm{w}}}} + {\lambda _{\rm{w}}}{u_{\rm{w}}}){u_{\rm{w}}}}}{{{\rho _{0{\rm{s}}}}(W - {u_{\rm{w}}})}}$

(3)

因此，只需要测得飞片（样品）击靶速度W和窗口的波后粒子速度u_w，就可以获得样品的粒子速度u_s（u_s=W–u_w）和对应的冲击波速度D_s。随后，由波系作用（见图1（b））的几何关系可知，样品材料Hugoniot状态的拉格朗日纵波声速C_L为

${C_{\rm{L}}} = \frac{{{D_{\rm{s}}}{h_{\rm{s}}}}}{{{D_{\rm{s}}}{t_{12}} - {h_{\rm{s}}}}}$

(4)

相应的欧拉纵波声速C_l为

${C_{\rm{l}}} = \frac{{{D_{\rm{s}}}{h_{\rm{s}}}}}{{{D_{\rm{s}}}{t_{12}} - {h_{\rm{s}}}}}\frac{{{D_{\rm{s}}} - {u_{\rm{s}}}}}{{{D_{\rm{s}}}}}$

(5)

式中：h_s为样品厚度，下标1和2对应波剖面上不同的时间点。

可以看到，反向碰撞法以波剖面测量为基础，波系作用简单，通过波剖面的粒子速度及时间信息得到高压声速，实验数据具有较高的精度，但是由于可供选择的窗口材料种类较少，目前使用的LiF窗口的阻抗较低，导致实验压力范围有限。

2. 实验结果与分析

实验在中国工程物理研究院流体物理研究所 $\varnothing$ 30 mm二级轻气炮上进行。将MB2合金飞片安装在弹丸上，将弹丸发射至稳定的弹道速度W，并撞击LiF单晶窗口，通过DPS测量飞片击靶速度以及飞片/窗口界面粒子速度，获得待测材料的冲击波速度和声速。为了提高测试界面对DPS入射光的反射效率，避免靶室残留气体对测试的干扰，窗口击靶面上镀1 ${\text{μ}}{\rm{m}}$ 厚的铝膜并贴8 ${\text{μ}}{\rm{m}}$ 厚铝箔。LiF窗口折射率修正采用Rigg等^[12]的公式

${u_{\rm{w}}} = 0.7895 \times {u^{0.9918}}$

(6)

式中：u_w代表修正后的窗口界面粒子速度，u代表实测界面粒子速度，二者单位均为km/s。其中LiF密度为2.638 g/cm³，冲击波速度D和粒子速度u的关系为D=5.150+1.352u（单位km/s）^[5]。

图3给出了6发实验测得的MB2/LiF界面粒子速度剖面（平台高度随加载压力的增大而升高）。通过界面粒子速度剖面，由（1）式～（6）式可得MB2样品在30～73 GPa冲击压力范围内的冲击波速度-粒子速度和声速-压力数据，实验结果列于表1，其中密度采用排水法测量，实测值为(1.775±0.004) g/cm³。冲击波速度-粒子速度数据如图4所示，图中还显示了美国洛斯阿拉莫斯国家实验室（LASL）发表的实验数据^[5]。从图4中可以看出，本研究获得的实验数据与已有实验数据具有较好的一致性。

表 1 MB2样品冲击实验参数及结果

Table 1. Shock experiment parameters and results of MB2

Exp.No.	h_s/mm	W/(km·s^–1)	u_w/(km·s^–1)	u_s/(km·s^–1)	D_s/(km·s^–1)	p/GPa	C_l/(km·s^–1)
1	1.980±0.004	3.949±0.020	1.589±0.016	2.360±0.026	7.303±0.166	30.6±0.4	7.983±0.297
2	1.997±0.004	4.358±0.020	1.763±0.018	2.595±0.027	7.607±0.175	35.0±0.5	9.101±0.394
3	1.978±0.004	5.379±0.027	2.195±0.022	3.184±0.030	8.317±0.191	47.0±0.7	9.167±0.402
4	1.981±0.004	5.928±0.030	2.435±0.024	3.493±0.039	8.746±0.214	54.2±0.9	8.912±0.380
5	1.981±0.004	6.100±0.030	2.514±0.025	3.586±0.039	8.907±0.213	56.7±0.9	8.601±0.348
6	1.987±0.004	7.220±0.036	3.003±0.030	4.217±0.050	9.747±0.245	73.0±1.2	9.723±0.463

下载: 导出CSV

| 显示表格

图 3 MB2/LiF界面粒子速度剖面

Figure 3. Particle velocity profile of MB2/LiF interface

下载: 全尺寸图片幻灯片

图 4 冲击波速度与粒子速度的关系

Figure 4. Shock velocity vs. particle velocity

下载: 全尺寸图片幻灯片

从图3可以看到：6发实验测得的界面粒子速度剖面质量良好，粒子速度剖面对应的冲击波、卸载波到达样品/窗口界面的特征信号清晰；加载压力为30.6和35.0 GPa时卸载剖面的弹-塑性特征明显，表明MB2合金在该冲击压力下尚未完全熔化；随着加载压力的升高，卸载剖面的弹-塑性特征逐渐消失，当加载压力达到73.0 GPa时，弹-塑性卸载特征完全消失，表明MB2已完全进入熔化相区，与图5所示的不同加载压力下声速转变特征一致。

图 5 声速与冲击压力的关系

Figure 5. Sound velocity vs. shock pressure

下载: 全尺寸图片幻灯片

如图5所示，当加载压力由30.6 GPa增加到35.0 GPa时，纵波声速逐渐增大，由7.983 km/s增大至9.101 km/s；但是，当加载压力增大至47.0 GPa时，纵波声速逐渐向体波声速偏转，跃变为9.167 km/s，预示着随着加载压力的升高，材料内部的剪切效应减小，冲击熔化发生；直至压力达到56.7 GPa时，纵波声速转变为体波声速，由此进一步确认了冲击加载下MB2合金发生熔化。Urtiew等^[13]通过理论预测MB2合金在57 GPa附近开始熔化，该结果与本研究根据声速判定的冲击熔化区域基本一致，从而进一步证实了MB2合金在该压力范围内发生冲击熔化，只是Urtiew等预估的理论压力略高。图5中的实线是根据以下公式^[11]计算得到的体波声速曲线

$C_{\rm{b}}^2 = \frac{1}{2}{V^2}{\left(\frac{\gamma }{V}\right)_{\rm{H}}}{p_{\rm{H}}} - {V^2}\frac{{{\rm{d}}{p_{\rm{H}}}}}{{{\rm{d}}V}}\left[1 - \frac{1}{2}{\left(\frac{\gamma }{V}\right)_{\rm{H}}}({V_0} - {V_{\rm{H}}})\right]$

(7)

式中：C_b为体波声速； $\rho$ 、 $\gamma$ 分别为密度和Grüneisen系数， $\rho \gamma = {\rho _0}{\gamma _0}$ ， $\rho _0$ 为材料初始密度， $\gamma _0$ =1.43^[14]。图5中的虚线是基于吴-经方程^[15]计算得到的体波声速曲线。可以看出，理论计算的体波声速曲线较冲击加载实验值偏低约10%。这主要是由于合金材料自身成键及结合能的影响，难以采用传统的混合法则或物理模型准确计算物态方程的基本参数，如Grüneisen系数等，由此导致理论预测结果与实验结果出现差异（通常理论值偏低）。

根据不确定度传递关系^[16]，当全部直接测量量（输入量）Y_i彼此独立不相关时，由其确定的间接测量量z的合成不确定度 ${\text{Δ}} {\textit{z}}$ 由下式确定

${\text{Δ}} {{\textit{z}}^2} = \sum\limits_{i = 1}^N {{{\left( {\frac{{\partial f}}{{\partial {y_{{i}}}}}} \right)}^2}{\text{Δ}} y_{{i}}^2}$

(8)

式中：y_i为输入量Y_i的直接测量值，z为被测量的测量值，f为z和y_i的函数关系，N为输入量的总个数， ${\text{Δ}} {y_{{i}}}$ 为直接测量值的不确定度。因而，当密度为 ${\rho _{0{\rm{s}}}}$ 、厚度为 ${h_{\rm{s}}}$ 的飞片撞击窗口时，样品/窗口界面处粒子速度跳跃，实验测得飞片速度 $W$ 、界面粒子速度 ${u_{\rm{w}}}$ 、时间间隔（剖面平台） ${t_{12}}$ ，相应的测量不确定度为 ${\text{Δ}} {\rho _{0{\rm{s}}}}$ 、 ${\text{Δ}} {h_{\rm{s}}}$ 、 ${\text{Δ}} W$ 、 ${\text{Δ}} {u_{\rm{w}}}$ 、 ${\text{Δ}} {t_{12}}$ ，而窗口Hugoniot参数（ ${\rho _{0{\rm{w}}}}$ 、 ${C_{0{\rm{w}}}}$ 、 ${\lambda _{\rm{w}}}$ ）的不确定度（ ${\text{Δ}} {\rho _{0{\rm{w}}}}$ 、 ${\text{Δ}} {C_{0{\rm{w}}}}$ 、 ${\text{Δ}} {\lambda _{\rm{w}}}$ ）已知。由于各测量量是独立测量的，（ ${C_{0{\rm{w}}}}$ ， ${\lambda _{\rm{w}}}$ ）相互作用项较小，可忽略不计，据此可根据不确定度传递律确定声速测量的不确定度。

如图6所示，就该例反碰撞实验（No.2）而言，影响声速测量不确定度的因素很多，包括样品初始密度、厚度、飞片速度、界面粒子速度、追赶时间、窗口材料冲击Hugoniot参数等。样品内部冲击压缩状态（如粒子速度、冲击波速度、冲击压力等）均通过飞片速度、界面粒子速度及窗口Hugoniot参数计算获得，影响声速测量不确定度的主要因素在于飞片速度、界面粒子速度及稀疏波追赶时间（平台时间），所占比例（即对（8）式中各平方项求和，下同）约为声速测量不确定度的99%，其余如初始密度、厚度等参量在当前诊断水平下对声速测量不确定度的贡献较小，约占总体的1%。在现有诊断条件下，飞片速度采用DPS直接测量，测量扩展不确定度不大于0.5%；界面粒子速度剖面的测量不确定度主要受平台区速度及稀疏波追赶时间的影响，影响因素包含干涉信号数据转换精度、窗口折射率修正、起跳及卸载时刻的判断等，综合而言，平台区速度测量的扩展不确定度不大于1%，追赶时间测量的扩展不确定度约6 ns。总体而言，传递至声速的测量扩展不确定度不超过5%。

图 6 反碰撞法测量声速实验中不确定度

Figure 6. Uncertainties for backward-impact experiment based on law of propagation

下载: 全尺寸图片幻灯片

3. 结　论

采用反向碰撞实验技术，结合具有高时空分辨率的DPS，获得了MB2合金在30～73 GPa压力范围内的冲击Hugoniot及声速数据。随着加载压力的升高，MB2合金纵波声速呈现出明显的向体波声速转变的趋势，预示着材料内部的剪切效应逐渐减小，冲击熔化发生，其相变压力区间为40～57 GPa。该实验结果与不同加载压力下卸载波剖面对应的弹-塑性转变特征完全一致，由此进一步确认了冲击加载下MB2合金熔化相变的发生。

感谢中国工程物理研究院流体物理研究所黄金、康强、叶素华、方茂林、向曜民、陈志云等在实验过程中给予帮助。

图 1 13.5～36.2 MPa评价指标

Figure 1. 13.5–36.2 MPa regression evaluation indicator

下载: 全尺寸图片幻灯片

图 2 51.0～62.8 MPa评价指标

Figure 2. 51.0–62.8 MPa regression evaluation indicator

下载: 全尺寸图片幻灯片

图 3 优化模型预测结果的误差对比

Figure 3. Comparison of errors in the prediction results of the optimisation model

下载: 全尺寸图片幻灯片

图 4 孔洞与裂缝纹理

Figure 4. Hole and crack texture

下载: 全尺寸图片幻灯片

图 5 孔洞计算函数

Figure 5. Hole calculation function

下载: 全尺寸图片幻灯片

图 6 碎石材质纹理

Figure 6. Gravel material texture

下载: 全尺寸图片幻灯片

图 7 烟雾材质纹理

Figure 7. Smoke material texture

下载: 全尺寸图片幻灯片

图 8 毁伤效果

Figure 8. Destruction effect

下载: 全尺寸图片幻灯片

图 9 可视化系统结构

Figure 9. Structure of the visualization system

下载: 全尺寸图片幻灯片

图 10 仿真模块之间的关系

Figure 10. Relationship of simulation modules

下载: 全尺寸图片幻灯片

图 11 弹体参数输入

Figure 11. Input of projectile parameters

下载: 全尺寸图片幻灯片

图 12 靶体参数输入

Figure 12. Input of target parameters

下载: 全尺寸图片幻灯片

图 13 数据传递接收

Figure 13. Data transmission reception

下载: 全尺寸图片幻灯片

图 14 数据对比

Figure 14. Comparison of data

下载: 全尺寸图片幻灯片

图 15 侵彻过程对比

Figure 15. Comparison of penetration processes

下载: 全尺寸图片幻灯片

图 16 毁伤对比

Figure 16. Comparison of destruction

下载: 全尺寸图片幻灯片

表 1 30.5 mm口径弹体侵彻51.0 MPa靶板（ρ=2300 kg/m³）深度实验数据

Table 1. Experimental data on the depth of penetration of a 30.5 mm projectile into a 51.0 MPa target slab (ρ=2300 kg/m³)

m/kg	v/(m·s⁻¹)	H_exp/m	H_LMC/m	$\delta_{\rm{LMC}}$ /%	H_MCT/m	$\delta_{\rm{MCT}}$ /%	H_LNC/m	$\delta_{\rm{LNC}}$ /%
1.60	405	0.37	0.28	25.02	0.28	23.87	0.23	37.79
	446	0.42	0.33	20.53	0.34	19.19	0.28	33.66
	545	0.56	0.49	12.69	0.50	10.90	0.41	26.04
	651	0.78	0.69	12.14	0.70	9.99	0.59	24.42
	804	1.05	1.04	1.10	1.07	1.94	0.91	12.93
	821	1.23	1.08	11.96	1.12	9.18	0.96	22.28
	900	1.41	1.30	7.53	1.35	4.30	1.17	17.37
	1 009	1.75	1.65	5.91	1.71	2.15	1.50	14.50
	1 069	1.96	1.85	5.39	1.93	1.36	1.70	13.24
	1 201	2.03	2.36	16.20	2.47	21.87	2.21	8.70

m/kg	v/(m·s⁻¹)	H_exp/m	H_HMC/m	$\delta_{\rm{HMC}}$ /%	H_HMCT/m	$\delta_{\rm{HMCT}}$ /%	H_HNC/m	$\delta_{\rm{HNC}}$ /%
1.60	405	0.37	0.29	21.88	0.32	13.51	0.24	34.55
	446	0.42	0.35	16.85	0.39	8.14	0.29	30.19
	545	0.56	0.52	7.79	0.57	1.40	0.44	22.17
	651	0.78	0.73	6.41	0.80	2.51	0.62	20.47
	804	1.05	1.12	6.48	1.22	16.13	0.96	8.40
	821	1.23	1.17	5.11	1.27	3.45	1.01	18.25
	900	1.41	1.41	0.10	1.54	8.99	1.23	13.11
	1 009	1.75	1.79	2.37	1.95	11.34	1.57	10.14
	1 069	1.96	2.02	3.15	2.20	12.18	1.79	8.85
	1 201	2.03	2.58	27.13	2.81	38.34	2.32	14.10

下载: 导出CSV

表 2 实验验证数据

Table 2. Experimental validation data

Case	Projectile			Target slab		v/(m·s⁻¹)	H/m
Case	m/kg	d/mm	$\delta_{\rm{CRH}}$	f_c/MPa	ρ/(kg·m⁻³)	v/(m·s⁻¹)	H/m
1	1.600	30.5	3	51.0	2300	545	0.560
2	1.600	30.5	3	51.0	2300	1 201	2.030
3	0.478	20.3	3	58.4	2320	610	0.491
4	0.478	20.3	3	58.4	2320	1 009	1.300
5	1.620	30.5	3	58.4	2320	445	0.460
6	1.620	30.5	3	58.4	2320	980	1.950
7	0.480	20.3	3	62.8	2300	821	0.760

下载: 导出CSV

表 3 系统测试数据对比

Table 3. Comparison of system test data

Case	H/m		Error/%
Case	Experiment	Systematic prediction	Error/%
1	0.560	0.642	14.6
2	2.030	2.434	19.9
3	0.491	0.500	1.8
4	1.300	1.185	8.8
5	0.460	0.428	6.9
6	1.950	1.697	12.9
7	0.760	0.819	7.7

下载: 导出CSV

参考文献(38)

[1]	金丰年, 刘黎, 张丽萍, 等. 深钻地武器的发展及其侵彻 [J]. 解放军理工大学学报(自然科学版), 2002, 3(2): 34–40. JIN F N, LIU L, ZHANG L P, et al. Development of projectiles and their penetration [J]. Journal of PLA University of Science and Technology (Natural Science Edition), 2002, 3(2): 34–40.
[2]	王涛, 余文力, 王少龙, 等. 国外钻地武器的现状与发展趋势 [J]. 导弹与航天运载技术, 2005(5): 51–56. WANG T, YU W L, WANG S L, et al. Development of deep-drilling weapons and their penetration [J]. Missiles and Space Vechicles, 2005(5): 51–56.
[3]	NING J G, REN H L, LI Z, et al. A mass abrasion model with the melting and cutting mechanisms during high-speed projectile penetration into concrete slabs [J]. Acta Mechanica Sinica, 2022, 38(10): 121597. doi: 10.1007/s10409-022-21597-x
[4]	YANKELEVSKY D, FELDGUN V. The embedment of a high velocity rigid ogive nose projectile into a concrete target [J]. International Journal of Impact Engineering, 2020, 144: 103631. doi: 10.1016/j.ijimpeng.2020.103631
[5]	NING J G, MENG F L, MA T B, et al. Failure analysis of reinforced concrete slab under impact loading using a novel numerical method [J]. International Journal of Impact Engineering, 2020, 144: 103647. doi: 10.1016/j.ijimpeng.2020.103647
[6]	CHEN X G, LU F Y, ZHANG D. Penetration trajectory of concrete targets by ogived steel projectiles experiments and simulations [J]. International Journal of Impact Engineering, 2018, 120: 202–213. doi: 10.1016/j.ijimpeng.2018.06.004
[7]	焦登伟. 混凝土/钢筋混凝土高速侵彻贯穿问题的数值模拟与实验研究 [D]. 北京: 北京理工大学, 2018. JIAO D W. Numerical simulation and experimental research on high-speed penetration of concrete/reinforced concrete [D]. Beijing: Beijing Institute of Technology, 2018.
[8]	CHEN X W, LI Q M. Deep penetration of a non-deformable projectile with different geometrical characteristics [J]. International Journal of Impact Engineering, 2002, 27(6): 619–637. doi: 10.1016/S0734-743X(02)00005-2
[9]	ROSENBERG Z, DEKEL E. Analytical solution of the spherical cavity expansion process [J]. International Journal of Impact Engineering, 2009, 36(2): 193–198. doi: 10.1016/j.ijimpeng.2007.12.014
[10]	WU H, HEN X W, HE L L, et al. Stability analyses of the mass abrasive projectile high-speed penetrating into concrete target. Part Ⅰ: engineering model for the mass loss and nose-blunting of ogive nosed projectiles [J]. Acta Mechanica Sinica, 2014, 30(6): 933–942. doi: 10.1007/s10409-014-0090-1
[11]	XU X Z, SU Y C, FENG Y B. Optimization analysis of state equation and failure criterion for concrete slab subjected to impact loading [J]. International Journal of Impact Engineering, 2024, 186: 104872. doi: 10.1016/j.ijimpeng.2023.104872
[12]	WANG M, WANG K L, ZHAO Q C et al. LQR control and optimization for trajectory tracking of biomimetic robotic fish based on unreal engine [J]. Biomimetics, 2023, 8(2): 236. doi: 10.3390/biomimetics8020236
[13]	丁钱. 基于模型驱动的复杂系统装备互操作可视化仿真与验证平台设计与实现 [D]. 南京: 南京邮电大学, 2018. DING Q. Design and implementation of a visual simulation and verification platform for equipment interoperability of complex system based on model driven [D]. Nanjing: Nanjing University of Posts and Telecommunications, 2018.
[14]	CHEN R C, AKKUS I E, FRANCIS P. SplitX: high-performance private analytics [C]//Proceedings of the ACM SIGCOMM 2013 Conference on SIGCOMM. Hong Kong, China: Association for Computing Machinery, 2013: 315–326.
[15]	SHEWCHENKO N, FONRNIER E, WONNACOTT M, et al. A vulnerability/lethality model for the combat soldier, a new paradigm—basis and initial development [C]//Personal Armour Systems Symposium. Nuremberg-Fürth, Germany, 2012.
[16]	贺歆. 弹目遭遇动态可视化仿真研究 [D]. 南京: 南京理工大学, 2004. HE X. Simulation study on dynamic visualization of projectile encounters [D]. Nanjing: Nanjing University of Science and Technology, 2004.
[17]	刘云. 弹目遭遇可视化仿真关键技术研究 [D]. 南京: 南京理工大学, 2004. LIU Y. Research on the key technology of visual simulation of projectile encounter [D]. Nanjing: Nanjing University of Science and Technology, 2004.
[18]	RUSYAK I, SUFIYANOV V, KOLOLEV S, et al. Software complex for simulation of internal and external ballistics of artillery shot [C]//International Conference on Military Technologies (ICMT). Brno, Czech Republic: IEEE, 2015.
[19]	BUŻANTOWICZ W. Matlab script for 3D visualization of missile and air target trajectories [J]. International Journal of Computer and Information Technology, 2016, 5(5): 419–422.
[20]	SCHOEDER W J, MARTIN, AVILA , et al. The visualization toolkit user’s guide [M]. New York, USA: Kitware Inc., 2001.
[21]	刘建斌, 徐豫新, 高鹏, 等. 火箭杀爆弹毁伤幅员仿真 [J]. 兵工学报, 2016, 37(Suppl 2): 159–164. LIU J B, XU Y X, GAO P, et al. Simulation of the extent of damage caused by rocket explosive projectiles [J]. Acta Armamentarii, 2016, 37(Suppl 2): 159–164.
[22]	高鹏, 徐豫新, 李可, 等. 火箭杀爆弹动爆毁伤幅员计算方法与可视化仿真实现 [J]. 兵工学报, 2016, 37(Suppl 2): 216–220. GAO P, XU Y X, LI K, et al. Calculation method and visual simulation implementation of kinetic explosion damage amplitude of rocket explosive projectiles [J]. Acta Armamentarii, 2016, 37(Suppl 2): 216–220.
[23]	YANG Y B, QIAN L X, MING Q U. Lethality test simulation system for conventional warhead based on numerical results [C]//Asia Simulation Conference/6th International Conference on System Simulation and Scientific Computing. Beijing, 2005: 267–271.
[24]	杨云斌, 屈明, 钱立新. 破片战斗部威力仿真方法与仿真软件研究 [J]. 计算机仿真, 2007, 24(10): 14–19. YANG Y B, QU M, QIAN L X. Lethality simulation method & software for fragmentation warhead [J]. Computer Simulation, 2007, 24(10): 14–19.
[25]	WANG Y, ZHANG W Y, NING J G. Streamline-based visualization of 3D explosion fields [C]//International Conference on Computational Intelligence and Security. Sanya, Hainan: IEEE, 2011.
[26]	LI B X, LIANG Z G, PENG S, et al. Bullet external ballistic visualization simulation software design [C]//International Conference on Machine Learning and Computer Application. Shenyang, Liaoning: IEEE, 2021.
[27]	张进强, 蒋夏军. 外弹道可视化仿真研究与实现 [J]. 计算机工程与科学, 2015, 37(2): 372–378. ZHANG J Q, JIANG X J. Research and implementation of visualization simulation for external trajectory [J]. Computer Engineering & Science, 2015, 37(2): 372–378.
[28]	程翔, 李苑青, 王丽华. 基于OpenGL的六自由度三维弹道仿真技术研究 [J]. 电子科技, 2017, 30(6): 15–20. CHENG X, LI W Q, WANG L H. Research on six degree of freedom 3D trajectory simulation based on OpenGL [J]. Electronic Science and Technology, 2017, 30(6): 15–20.
[29]	王彪, 程鹏举, 韩卓茜, 等. 基于C#与STK的多雷达跟踪弹道导弹系统设计与实现 [J]. 计算机科学与应用, 2020, 10(6): 1185–1193. doi: 10.12677/CSA.2020.106123 WANG B, CHENG P J, HAN Z X, et al. Design and implementation of multi-radar tracking ballistic missile system based on C# and STK [J]. Computer Science and Application, 2020, 10(6): 1185–1193. doi: 10.12677/CSA.2020.106123
[30]	纪录, 吴国东, 王志军, 等. 基于STK的弹箭半实物飞行实时可视化仿真 [J]. 火力与指挥控制, 2020, 45(2): 170–174, 179. JI L, WU G D, WANG Z J, et al. Missile-rocket semi-physical flight simulation based on STK true time visualization technology [J]. Fire Control & Command Control, 2020, 45(2): 170–174, 179.
[31]	TIAN X J, TAO T J, LOU Q X, et al. Modification and application of limestone HJC constitutive model under the impact load [J]. Lithosphere, 2022, 2021: 6443087. doi: 10.2113/2022/6443087
[32]	FELDGUN V R, YANKELEVSKY D Z. Constitutive equations for reliable projectile penetration analysis into a concrete medium [J]. International Journal of Protective Structures, 2019, 11(2): 159–184.
[33]	DONG H, WU H J, LI J Z, et al. Dynamic spherical cavity expansion analysis of concrete/rock based on hoek-brown criterion [J]. Journal of Physics: Conference Series, 2020, 1507: 032012. doi: 10.1088/1742-6596/1507/3/032012
[34]	YANG H W, ZHANG J, WANG Z Y, et al. Numerical study on the resistance of rigid projectiles penetrating into semi-infinite concrete targets [J]. Acta Mechanica Sinica, 2021, 37(3): 482–493. doi: 10.1007/s10409-021-01054-6
[35]	LI Z, XU X Z. Theoretical investigation on failure behavior of ogive-nose projectile subjected to impact loading [J]. Materials, 2020, 13(23): 5372. doi: 10.3390/ma13235372
[36]	JADON A, PATIL A, JADON S. A comprehensive survey of regression based loss functions for time series forecasting [J/OL]. Arxiv, 2022: 2211.02989.
[37]	FORRESTAL M J, ALTMAN B S, CARGILE J D, et al. An empirical equation for penetration depth of ogive-nose projectiles into concrete targets [J]. International Journal of Impact Engineering, 1994, 15(4): 395–405. doi: 10.1016/0734-743X(94)80024-4
[38]	FORRESTAL M J, FREW D J, HANCHAK S J, et al. Penetration of grout and concrete targets with ogive-nose steel projectiles [J]. International Journal of Impact Engineering, 1996, 18(5): 465–476. doi: 10.1016/0734-743X(95)00048-F