Optimization Model and Visualization Simulation of Projectile Penetration into Concrete
-
摘要: 利用可视化仿真技术研究弹体侵彻混凝土的毁伤机理和靶板响应是爆炸冲击领域的重要课题。混凝土作为常见的建筑材料,在遭受爆炸冲击或高速弹体侵彻时,其毁伤行为复杂多变。介绍了一种理论研究与可视化技术相结合的可视化仿真方法。基于空腔膨胀理论建立了优化的侵彻计算模型,可以预测弹体侵彻混凝土的侵彻深度特征。利用可视化物理引擎,对弹体运动轨迹、开坑孔径、靶板损伤、碎石飞溅等进行了细致的表征处理和仿真,增强了场景的真实性和可靠性。开发的可视化仿真系统不仅能够从多角度观察弹体侵彻混凝土的过程,而且能够高效、准确地分析和预测弹体侵彻混凝土靶的损伤行为和动力响应,在建筑工程设计和安全评估中具有重要的应用前景,为理解和探索混凝土侵彻机理提供了新的视角。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.
-
Key words:
- projectile penetration /
- concrete /
- visual simulation /
- physics engine
-
由于爆炸焊接过程具有瞬时性、复杂性,人们在通过实验及相关设备来观察、研究爆炸焊接过程时,效果十分有限。随着计算机技术的发展,数值模拟技术的突破为爆炸焊接研究提供了帮助。迄今为止,国内外诸多学者利用数值模拟软件对爆炸焊接过程及参数开展了大量研究,取得了较理想的研究成果。Nassiri等[1]分别利用光滑粒子流体动力学(Smoothed particle hydrodynamics,SPH)法、任意拉格朗日-欧拉(Arbitrary Lagrange-Euler,ALE)法模拟了爆炸焊接的高速撞击过程,并研究了结合界面的剪切应力和速度分布,结果显示只有SPH法再现了碰撞时的射流现象。Abe[2]使用二维有限差分法研究了爆炸焊接波的形成机理,计算结果与模拟结果保持一致。Yuan等[3]采用SPH法研究了6061Al/AZ31B的爆炸焊接,模拟得到波形界面和射流现象,但是并没有详细指出波的形成机理。Tabbataee等[4]使用有限单元法模拟了射流。Mousavi等[5-6]采用欧拉法模拟了射流和类似流体状的波状界面,并成功预测了波形尺寸和射流速度。王宇新等[7]应用无网格物质点法(Material point method,MPM),对爆炸荷载作用下基板和复板的动态变形过程进行了三维数值模拟,并分析了碰撞点压力和速度的变化情况。刘江等[8]用SPH法进行二维建模,重现了实验中观察到的射流和界面波现象,并揭示了碰撞压力、有效塑性应变、温度等参数随时间的变化关系。周春华等[9-10]研究了均匀布药下结合界面波的变化规律后,提出了新的不等药厚布药工艺。
从众多学者的研究成果来看,数值模拟软件能够帮助人们全面、系统地了解爆炸焊接过程,降低实验成本,优化实验参数。在平板爆炸焊接中,通过炸药爆炸产生的爆轰能量推动复板与基板发生斜碰撞,碰撞界面在微秒级的时间内产生数千兆帕高压,发生塑性变形,从而实现基板和复板的固相焊接。由此可以看出,炸药是整个爆炸焊接过程的能量来源,因此对炸药药量和布药方式进行控制具有重要的研究意义。本研究应用LS-DYNA软件,采用SPH法进行二维建模,研究炸药与爆炸焊接界面波的关系及波的形成机理。
1. 计算模型
1.1 爆炸几何模型
SPH法是一种无网格法,适合处理爆炸、冲击等大变形问题,但是其搜索算法较为耗时,粒子越多,该问题越突出,计算效率越低[11]。为了提高计算结果的精确性,同时兼顾计算效率,炸药采用粗网格划分,基板和复板采用细网格划分。先采用LS-DYNA软件建立二维平面计算模型,再利用LS-prepost软件将网格全部生成SPH粒子(见图1)。复板和基板材料选取Cu和Q235钢,爆炸焊接选用硝铵炸药。模型尺寸:复板150 mm × 4 mm,基板150 mm × 15 mm,间隙6 mm,基板和复板网格尺寸为0.2 mm × 0.2 mm,炸药网格尺寸为0.4 mm × 0.4 mm。采用cm-g-μs单位制。
1.2 材料模型和状态方程
硝铵炸药选用Jones-Wilkins-Lee(JWL)[12]状态方程(EOS),计算公式为
p=AJ(1−ωR1V)e−R1V+BJ(1−ωR2V)e−R2V+ωE0V (1) 式中:AJ、BJ、R1、R2和ω为材料系数;p为爆轰产物压力,GPa;E0为初始比内能,kJ·cm–3;V为爆轰气体产物的相对比容。炸药的具体JWL状态参数见表1,其中ρ为密度,D为爆速。
表 1 硝铵炸药的JWL状态方程参数Table 1. JWL EOS parameters of ammonium nitrate explosiveρ/(kg·m−3) D/(m·s−1) AJ/GPa BJ/GPa R1 R2 ω 800 2 800 132.75 0.423 5.3 1.2 0.21 基板、复板均采用 Mie-Grüneisen状态方程[13]和Johnson-Cook材料模型[14]。Grüneisen状态方程被广泛应用于解决材料的大变形问题,它描述了粒子速度up和冲击速度D之间的基本关系,计算公式如下
ρ=ρH+Γρ(e−eH) (2) 其中
Γρ=Γ0ρ0=Constant (3) ρH=ρ0c0μ(1+μ)[1−(s−1)μ]2 (4) eH=12ρHρ0μ1+μ (5) μ=ρρ0−1 (6) 式中:Γ0为Grüneisen系数,ρ为当前密度,ρ0为初始密度,c0为体积声速,s为压缩比。
Johnson-Cook材料模型经常用来反映金属的力学行为,计算公式如下
σ=(A+Bεnp)(1+Cln˙ε∗p)(1−T∗m) (7) 其中无量纲温度T*表示为
T∗=(T−Tr)/(Tm−Tr) (8) 式中:εp为有效塑性应变;
˙ε∗p=˙εp/˙ε0p 为有效塑性应变率,˙ε0p 为参考应变率;A、B、C、m、n是与材料相关的常数;Tr为室温;Tm为熔点。Cu和Q235钢的Johnson-Cook材料模型参数列于表2,Grüneisen EOS参数列于表3,其中:c为声速,km/s; S1为拟合系数;a为Grüneisen系数Γ0的一阶体积修正,G为剪切模量。表 2 Cu和Q235钢的Johnson-Cook模型参数Table 2. Parameters of Johnson-Cook model of Cu and Q235 steelMaterial ρ/(g·cm−3) G/GPa A/GPa B/GPa n C m Tm/K Tr/K Cu 8.96 46 0.090 0.292 0.31 0.025 1.09 1 356 294 Q235 7.83 77 0.792 0.510 0.26 0.014 1.03 1 793 294 表 3 Cu和Q235钢的Grüneisen方程参数Table 3. Grüneisen EOS parameters of Cu and Q235 steelMaterial c/(km·s−1) S1 Γ0 a Cu 3.940 1.489 2.02 0.47 Q235 4.569 1.490 2.17 0.46 1.3 爆炸复合窗口
爆炸焊接质量的好坏极其依赖选取的爆炸焊接参数,如药量比、间距、碰撞速度、碰撞角等。参考文献[15]中的爆炸复合窗口,选取复板厚度δ1 = 4 mm时药量比(单位面积炸药与复板质量之比R)分别为R1 = 1.0、R2 = 1.5的两组参数开展研究。经计算得到对应的炸药厚度分别为44.8 mm和67.2 mm。由图2可知,R1靠近可焊窗口下限,R2位于可焊窗口内部,两组参数均位于爆炸复合窗口内部。
2. 模拟结果分析
2.1 均匀布药
为了研究炸药量与焊接界面波形的关系,首先进行均匀布药模拟。分别建立炸药厚度为44.8 mm和67.2 mm的计算模型,输出模拟图像,如图3所示。
为了更方便地了解爆炸焊接过程参数的变化规律,在结合面上等间距区域选取7个关键点,选取位置如图4所示。输出关键点的碰撞压力值,如表4所示。硝铵炸药的爆轰压力计算公式为
表 4 均匀布药方案关键点碰撞压力Table 4. Collision pressure of key points in uniform charge schemeKey point Pressure/GPa Key point Pressure/GPa R1 = 1.0 R2 = 1.5 R1 = 1.0 R2 = 1.5 A1 0.581 1.602 A5 4.999 7.086 A2 1.672 2.002 A6 5.754 8.576 A3 4.507 5.289 A7 0.526 1.031 A4 4.654 6.191 p=ρ0D21+γ (9) 式中:γ为硝铵炸药的多方指数,ρ0为炸药密度,D为炸药爆速。
由式(9)可以看出,模拟中使用的硝铵炸药的密度为常数,所以爆轰压力p只与炸药爆速D有关。由于炸药的爆速与单位面积的药量成正比[16],且R1 < R2,因此R2对应的爆轰压力要高于R1的爆轰压力。炸药的爆轰压力作用于复板,复板又传递到焊接界面,最终导致R2对应的基板、复板碰撞压力高于R1。根据图5的模拟结果可知,理论分析结果与模拟结果一致。
焊接界面波长的计算公式[17]为
λ=πvpd2c(cvc−1)2 (10) 式中:vp为碰撞速度;d为复板厚度;c为材料声速;vc为碰撞点移动速度,等于炸药的爆轰速度。
从式(10)可以看出,给定金属板材的材料参数后,复板厚度d和材料的声速c为常数,碰撞速度vp和碰撞点移动速度vc为变量,界面波的波长λ与vp成正比,与vc成反比。
从图3可以看出,复板上部均出现了打伤现象,这是由于前期炸药和板材划分的网格尺寸不同所致。R1和R2两组模拟均得到了波状焊接界面,这被认为是焊接良好的标志[18],模拟结果与文献[19]的金相实验结果吻合较好。从图3(a)和图3(b)看出,焊接界面在起爆端区域没有出现波状界面,随着距起爆端距离的不断增大,界面先后出现了无波到微波、微波到大波的转变。图3(a)中射流现象最先出现在4.299 6 μs时刻,而图3(b)中的射流最先出现在4.198 2 μs,早于前者,这是由于药量比越大,爆轰能量越大,导致射流出现较早。焊接界面的放大图像显示,波状界面先于射流出现,可以认为波的形成不一定需要射流,印证了Godunov提出的应力波机理[20]。从图3(c)、图3(d)可以清晰地看出,碰撞点前方在基板和复板碰撞作用下,板材表面喷射出大量射流,可以清洗掉金属表面的杂质和氧化层,提高焊接质量[21]。喷射出的射流中,基板的成分明显高于复板,基板材料是射流成分的主要来源,这是由于Cu的密度高于Q235钢,密度较小的材料更容易产生射流[22]。图3(e)和图3(f)为模拟图像,表明基板与复板已经完全复合。从复合板中央部分的取样结果来看,图3(f)的波状尺寸大于图3(e),由于小波增加了焊接面积,有助于提高焊接质量,波形越大,越容易出现裂缝和孔洞等缺陷,因此可以认为药量比R1优于R2的焊接质量,药量比尽量取在可焊窗口的下限。
图5是均匀布药下各个关键点的碰撞压力折线。考虑到边界效应的影响,去掉首尾两个关键点后,折线显示沿着爆轰方向焊接界面的碰撞压力逐渐增大,由于碰撞压力过大会产生熔化层,产生的稀疏波可能将已结合的复合板拉开,降低焊接质量,因此均匀布药不是最终选择。
2.2 梯形布药
根据均匀布药方式下界面波的分布特点,起爆端区域的结合呈平直界面,这是由于炸药起爆时能量不稳定[23],碰撞压力小。在爆炸焊接中,碰撞压力过小会导致界面未复合[16],即起爆端附近极易出现未复合现象,降低金属板材的利用率,这种现象被称为边界效应,在实验和工程实践中应尽量避免。碰撞压力要超过某个临界值,复合板才能获得良好的焊接质量,这也证明了爆炸焊接是一种压力焊[24]。沿着爆轰方向,结合界面逐渐由无波发展成微波、大波界面,这是爆炸产物的堆集以及基板和复板碰撞引起的振动能在待复合区叠加并共同作用的结果[25-26]。结合式(9)、式(10)可知,波状界面的尺寸与碰撞点的移动速度有关,碰撞压力与炸药爆速成正比,在平行法爆炸焊接装置中爆速等于碰撞点的移动速度[27],由此得出,炸药的爆轰能量赋予了复板巨大的压力,碰撞压力决定界面波形。因此可以通过调整炸药厚度来控制爆轰能量,使界面波形维持在微波状,进而提高复合板的焊接质量。
设计如图6所示的呈梯形装药结构,通过调整炸药起爆端和末端厚度a、b的值,来确定一个合适的倾斜装药角度,沿着爆轰方向控制碰撞点的压力基本保持不变,避免因大波状界面而产生涡旋和裂缝,提高复合板的结合强度,使整个焊接界面都保持在微波状结合状态。布药方案见表5。
表 5 梯形布药方案Table 5. Ladder charging schemeScheme a/mm b/mm Ⅰ 67.2 58.8 Ⅱ 67.2 50.4 Ⅲ 67.2 42.0 Ⅳ 67.2 33.6 由于网格划分较大,导致波形尺寸不易测量。当网格划分较小时,可以较清晰地观察界面的形貌特征。由前述分析得知,碰撞压力决定界面波形,因此将界面上的碰撞压力作为研究对象,使界面上各点处碰撞压力基本保持一致。
不同梯形布药方案下界面关键点的碰撞压力如图7所示,具体压力值列于表6。由图7中折线的波动情况可知,当采用梯形布药方式后,基板、复板之间的碰撞压力不再呈均匀布药下逐渐增长的趋势,而是呈较小幅度的上下波动,因此梯形布药方式可以消除不利因素,使焊接界面的碰撞压力基本保持不变,以获得波长较小、波幅较小的微波界面。观察图7发现,起爆端和末端的碰撞压力最小,这是由稀疏波引发的边界效应所致。在不考虑首尾两端的压力时,方案Ⅲ的碰撞压力波动最小,因此最终焊接界面的波形状态好于其他方案。而方案Ⅱ中界面上的碰撞压力起伏最大,最终可能导致基板与复板未复合。因此确定合适的首尾装药高度是梯形布药方案的关键。
表 6 梯形布药方案关键点碰撞压力Table 6. Collision pressure of key points of ladder charge schemeKey point Pressure/GPa Scheme Ⅰ Scheme Ⅱ Scheme Ⅲ Scheme Ⅳ A1 1.738 1.158 1.694 0.351 A2 2.549 2.548 2.141 1.936 A3 3.839 3.830 3.621 2.776 A4 5.230 7.547 5.439 5.665 A5 7.890 8.330 3.413 4.537 A6 6.423 3.505 3.731 2.064 A7 0.780 0.431 0.397 0.435 综上所述,采用梯形布药方式可以消除焊接界面波形不均匀现象,而且减少了炸药用量,同时也可以减小爆炸焊接实验现场引发的噪声。通过分析各个关键点的碰撞压力可知,方案Ⅲ的复合效果最好,即炸药起爆端和末端的厚度分别为67.2 mm和42.0 mm时装药倾斜角度最佳。在实际的工程实践中,梯形布药方式的应用并不广泛,原因是其影响因素较多,参数确定比较复杂。尽管减小炸药厚度会使碰撞压力减小,但碰撞压力的变化与炸药厚度并非成严格的正比关系,所以才会得到图7的模拟结果。如何优化梯形布药参数使界面波波形基本保持不变,还需要大量深入的研究。
3. 界面波形成机理
爆炸焊接过程具有瞬时性,因此很难用实验方法来研究界面波的形成与焊接参数之间的关系。虽然有关界面波形成机理的研究已经进行了数十年,但迄今为止仍然没有达成一个统一的认识[28]。期间众多学者提出了一些波形形成机制,主要有复板流侵彻机理[29]、Helm-holtz不稳定流机理[30-31]、涡脱落机理[32-33]和应力波机理[34]。随着数值模拟软件的发展,SPH法能够清晰地模拟出波形成过程,为解释界面波的形成机理提供了帮助。
图8展示了采用SPH法得到的爆炸焊接界面波形成过程。基板与复板在P点高速碰撞,此时P点可看作低黏塑性流体状态,巨大的碰撞压力使P点产生塑性变形。如图8(a)所示,复板与基板碰撞后形成初始射流和再入射流,再入射流在前方复板的阻碍下向斜上方喷出。沿着焊接方向向前发展,碰撞点P转移到波峰上,见图8(b),先前的再入射流被基板捕获,形成在波峰的前后方,此处极易形成漩涡、空洞等缺陷,从而降低复合板的焊接质量。碰撞点P转移到波峰上后,再入射流与波峰碰撞,沿斜下方喷出。图8(c)中碰撞点P越过波峰后形成新的碰撞点,在巨大的爆轰压力下基板再次发生塑性变形,形成新的波形,由此产生周期性的波状界面。由图8可见,SPH法模拟的波形成过程与复板流侵彻机理的一致性较好,证明了SPH法在解释爆炸焊接界面波形成过程中复板流侵彻机理的有效性。
4. 结 论
(1)采用LS-DYNA软件和SPH法,模拟了整个爆炸焊接过程,得到典型的波状界面和射流现象,理论分析与模拟结果具有一致性。
(2)在两种药量比、均匀布药方案中,忽略边界效应后的结果显示,沿着爆轰方向基板与复板之间的碰撞压力均逐渐增大,而且药量比越大,碰撞压力越大,界面波形越大。
(3)在梯形布药方案中,通过改变炸药起爆端和末端的高度设计了4种方案,结果显示梯形布药可以有效消除爆炸焊接界面波的不均匀现象,使界面波形尺寸基本保持一致,而且节省了炸药用量。当起爆端和末端的高度分别为67.2 mm和42.0 mm时效果最佳。
(4)介绍了几种波形成机理,展示并研究了界面波形成过程,模拟的波形成过程与复板流侵彻机理的一致性较好,证明了复板流侵彻机理解释界面波形成过程的有效性。
-
表 1 30.5 mm口径弹体侵彻51.0 MPa靶板(ρ=
2300 kg/m3)深度实验数据Table 1. Experimental data on the depth of penetration of a 30.5 mm projectile into a 51.0 MPa target slab (ρ=
2300 kg/m3)m/kg v/(m·s−1) Hexp/m HLMC/m δLMC/% HMCT/m δMCT/% HLNC/m δ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−1) Hexp/m HHMC/m δHMC/% HHMCT/m δHMCT/% HHNC/m δ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 表 2 实验验证数据
Table 2. Experimental validation data
Case Projectile Target slab v/(m·s−1) H/m m/kg d/mm δCRH fc/MPa ρ/(kg·m−3) 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 表 3 系统测试数据对比
Table 3. Comparison of system test data
Case H/m Error/% Experiment Systematic prediction 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 -
[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-Fü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.106123WANG 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 -