Numerical simulation on ground shock waves induced by hypervelocity penetration of a projectile into a limestone target
-
摘要: 为探究超高速动能武器的对地破坏效应及其影响因素，采用数值模拟方法对弹体超高速侵彻的地冲击规律进行了研究。首先，基于石灰岩静动态力学性能实验数据对材料模型参数进行了标定，并对已有弹体大范围着速侵彻石灰岩靶体进行了模拟，验证了所采用材料模型和数值模拟方法的合理性。随后，开展了钨合金长杆弹超高速侵彻石灰岩靶体的数值模拟，细致分析了地冲击传播的现象和机理：弹体超高速侵彻靶体时，弹靶交界面处会产生瞬时高压，并以应力波的形式在靶体中传播，对靶体内部造成破坏，且当弹体初速度高于3.0 km/s时，地冲击显著增强。最后，进一步研究了不同弹靶参数对地冲击的影响，发现从相对深度来看，弹体参数（弹体长径比、密度）对地冲击规律影响不大；而靶体特征特别是孔隙率对地冲击传播具有较大影响。Abstract: With the advancement of hypervelocity weapons such as the “Rods-from-God”, the damage and failure in targets induced by the hypervelocity penetrators have been a topic of current research, which is still not fully understood. To address this issue, numerical investigation was carried out on ground shock induced by hypervelocity penetration of projectiles into limestone targets. As the material model and corresponding parameters are crucial for the accurate numerical predictions, the parameters for the p-α equation of state and the Kong-Fang material model recently proposed to describe the limestone were firstly calibrated based on a large amount of dynamic tests. The smooth particle hydrodynamics (SPH) method was employed for simulating the target and the axisymmetric numerical model was used to improve the computational efficiency. The calibrated parameters and numerical algorithm were validated by numerically simulating a series of penetration tests on limestone targets with a broad range of striking velocities. Then, based on the validated numerical model, the penetration of long-rod tungsten projectile into a limestone target was simulated and the mechanism of the corresponding ground shock was discussed. It was found that a high pressure in the target was induced by the hypervelocity impact of the projectile, which then propagated into the target as a stress wave, leading to the damage and failure in the target. The induced ground shock wave increased with the increase of the initial projectile velocity, especially when the initial projectile velocity is over 3.0 km/s. Finally, parametric study was conducted to investigate the effects of the parameters related to the projectiles and targets on the ground shock wave. The parameters related to the projectiles, for example, the length-to-diameter ratio and density, which can influence on the damage area of the targets by influencing the depth of penetration, has limited influence on the ground shock wave from the view of the relative depth (the ratio of the depth to the penetration depth). While the target parameters, especially the porosity which can affect the wave propagation, have a great effect on the ground shock wave.
-
Key words:
- hypervelocity penetration /
- ground shock wave /
- limestone /
- material model
-
近年来，以“上帝之杖”为代表的超高速动能武器迅速发展，其终点毁伤效应成为武器研发及工程防护领域关注的热点问题。弹体对地打击的破坏模式主要表现为两个方面：(1) 弹体冲击靶体，侵入靶体直接造成破坏；(2) 侵彻冲击产生的瞬时高压以应力波的形式在靶体中传播，造成靶体内部的间接破坏效应（如开裂，震塌等破坏现象），即为侵彻的地冲击效应。已有实验表明，弹体侵彻混凝土和岩石等脆性靶体时会呈现刚性侵彻、变形不侵蚀侵彻和侵蚀侵彻3个典型阶段[1]，且会出现侵彻深度随初速度增加而减小的情况[2-3]。弹体超高速侵彻的直接侵彻深度相较于中低速侵彻没有明显增加，地冲击效应可能成为其对地破坏的主要途径。因此，研究弹体超高速侵彻的地冲击规律，对综合评价超高速动能武器的对地毁伤效应具有重要意义。
“上帝之杖”动能弹的入地侵彻速度可超过3 km/s[4]，但受限于实验条件，对弹体在该速度区间侵彻岩石类脆性材料的研究较少。Poelchau等[5]、李干等[6]、钱秉文等[2, 7]分别开展了不同弹体超高速侵彻岩石类脆性材料的实验研究，主要针对靶体开坑形态、侵彻深度的变化规律等进行研究。孔祥振等[8-9]对弹体超高速侵彻半无限金属靶的6种典型的理论模型进行了对比分析，发现其均不适用于混凝土等脆性靶体的侵彻深度预测。特别是针对弹体超高速侵彻岩石类材料的地冲击研究，目前还较匮乏且主要以数值模拟方法为主。基于GEODYN流体动力学计算软件和自定义的材料动态本构关系，Antoun等[10]开展了钨合金弹体超高速侵彻石灰岩靶体的数值模拟研究，重点探讨了初速度为4和6 km/s情况下的地冲击效应。邓国强等[11]以流体弹塑性模型为基础，采用无网格的光滑粒子流体动力学（smooth particle hydrodynamics, SPH）方法，对弹体超高速侵彻花岗岩靶体进行了数值模拟研究，分析了超高速打击的破坏现象及靶体中应力波的传播规律，并基于破坏特征和应力波传播特征提出了超高速侵彻和装药浅埋爆炸的等效方法[12]。基于超高速侵彻成坑与浅埋爆炸成坑的对比分析，王明洋等[13]提出了超高速侵彻地冲击的等效计算方法，并结合开展的合金钢弹体超高速侵彻花岗岩靶体实验对等效方法进行了验证。可以发现，不同学者选用的弹靶材料、本构模型参数以及数值模拟方法均不相同，且对地冲击的影响因素尚未有系统的分析。
本文中，以钨合金长杆弹和石灰岩靶体为研究对象，对超高速侵彻的地冲击进行数值模拟。由于超高速侵彻时石灰岩靶体处于高围压、高应变率状态，因此为保证数值模拟的可靠性，首先基于石灰岩材料的静动态力学实验数据对材料模型参数进行详细标定，并结合弹体大范围着速侵彻石灰岩靶体实验验证模型参数的准确性；然后基于上述验证的数值模型，对弹体超高速侵彻石灰岩靶体进行数值模拟，研究不同初速度下的地冲击规律，并进一步详细探讨弹靶参数对地冲击的影响规律。
1. 石灰岩材料模型参数的标定
石灰岩是一种多孔、中等强度的岩石材料，抗压强度约为60 MPa，初始密度和压实密度分别为2 280和2 700 kg/m³。已有石灰岩的静动态力学性能实验数据较多[14-15]，因此靶体选用石灰岩，数值模拟中材料模型选用前期自主开发的Kong-Fang模型。本节中首先简要介绍Kong-Fang材料模型，然后基于已有的静动态力学性能实验数据对模型参数进行标定。
1.1 Kong-Fang模型
在对已有脆性材料模型（HJC、RHT和K&C模型等）进行精细化评估的基础上，Kong等[16]提出了适用于爆炸冲击作用下的Kong-Fang材料模型，并将其植入LS-DYNA软件中。Zhang等[17]结合高速侵彻下脆性靶体材料的高围压屈服和动态拉伸特征，对Kong-Fang模型进行了改进，并将改进的Kong-Fang模型植入AUTODYN商业软件中，经验证能够较好地模拟脆性材料靶体的损伤破坏效应。本文中，模型介绍及数值模拟均基于改进后的Kong-Fang模型。
Kong-Fang模型的强度面包含最大强度面σm、屈服强度面σy和残余强度面σr，并通过损伤变量D和η 插值得到当前破坏强度面σ ：
σ={r[D(σr−σm)+σm]p≤0,tensiler[η(σm−σy)+σy]λ≤λm,p＞0,strainhardeningr[η(σm−σr)+σr]λ＞λm,p＞0,strainsoftening (1) 式中：r为当前子午面与压缩子午面的比值；p为压力，受压为正； λ 为修正的等效塑性应变；λm为峰值应力对应的λ。
Kong-Fang模型能够区分拉压损伤，总损伤D由拉压损伤耦合而成：
D=1−(1−Dc)[1−(Dt+ΔDt)] (2) 式中：Dc为压缩损伤，Dt为应力偏量造成的拉伸损伤，ΔDt为应力球量造成的拉伸损伤。强度面及损伤相关的参数标定详见1.2节。
Kong-Fang模型使用径向放大破坏强度面的方法考虑应变率效应，最终强度面按照下式确定：
σ=rfσ(p/rf) (3) 式中：rf为材料的动态强度增强因子，在Kong-Fang模型中采用Xu等[18]提出的经验公式，具体形式及参数标定详见1.3节。
1.2 强度面参数的标定
Kong-Fang模型中最大强度面σm需通过围压压缩实验数据标定，屈服强度面σy和残余强度面σr分别通过材料的应变硬化和应变软化曲线特征进行标定，各强度面具体形式为：
σm={3[p/(1−D)+T]p≤01.5/ψ(p+T)0＜p≤fc/3fc+p−fc/3a1+a2p/a2pfcfcp＞fc/3 (4) σr={0p＜00.7fc(3p/fc)a30≤p≤fc/30.7fc+p−fc/3a1+a2p/fcp＞fc/3 (5) σy={(3−9Tfyc)p+3T0≤p≤fyc/3fyc+p−fyc/3a1y+a2yp/fcp＞fyc/3 (6) 式中：T、fc和fyc分别为材料的单轴静态拉伸强度、压缩强度和屈服强度。a1、a2、a3、a1y和a2y为强度面中的经验参数，需根据实验数据进行标定；ψ 为拉压子午面之比，需要根据围压压缩（σ1＜σ2=σ3，σ1、σ2、σ3分别为第一、第二和第三主应力，以拉为正）和围压拉伸（σ1＞σ2=σ3）两种状态的强度面关系确定。
利用Fossum等[19]开展的三轴实验数据对原Kong-Fang模型的拉压子午面比ψ 标定为：
ψ(p)={0.5p≤00.5+1.5T/fcp=fc/31.15/[1+1.3/(3a1+2.3a2fc)]p=2.3fc/3p/(p+fc)p≥3fc (7) 根据Fossum等[19-20]和Walton等[21]开展的不同围压下的石灰岩围压压缩实验数据，拟合得到最大强度面曲线，如图1所示。同时考虑岩石材料的剪切饱和特性，限定强度面最大值为20fc。
图 1 最大强度面模型与实验数据的拟合Figure 1. Fitting of the maximum strength surface model to the experimental data根据石灰岩压缩实验的应力-应变曲线，取屈服强度fyc=0.8fc，并对屈服强度面、残余强度面及损伤相关的参数（见表1）进行拟合。基于标定的参数开展单轴压缩（UUC）及围压压缩（TXC）数值模拟，其应力-应变曲线与实验结果[21]的对比如图2所示。可以看出，Kong-Fang模型能够较好地模拟出石灰岩材料的应变硬化和应变软化特征，特别是石灰岩低围压下呈现脆性破坏而高围压下呈现塑性破坏这一特性。
表 1 石灰岩的强度模型参数Table 1. Parameters of the strength surface models for limestone最大强度面 残余强度面 屈服强度面 损伤参数 a1 a2 a3 a1y a2y λm 0.877 0.022 0.80 0.890 0.046 3×10−5 
图 2 石灰岩应力-应变曲线的数值模拟结果与实验数据的对比Figure 2. Comparison of stress-strain curves of limestone between experimental data and numerical simulation1.3 应变率效应的参数标定
Kong-Fang模型采用Xu等[18]提出的应变率增强因子计算公式：
rf=tanh[(lg˙ε−Wx)S][(Fm/Wy−1)+1]Wy (8) 式中：Fm、Wx、S和Wy为相关经验参数。利用已有的实验数据对石灰岩的压缩动态强度增强因子（dynamic increase factor of compression, DIFC）和拉伸动态强度增强因子（dynamic increase factor of tension, DIFT）分别进行标定，其参数如表2所示，应变率增强因子的曲线与相关实验结果[22-27]的对比如图3所示。
表 2 石灰岩的动态强度增强因子参数Table 2. Parameters for dynamic increase factors of limestone动态强度增强因子 Fm Wx S Wy DIFC 9 1.8 1.2 5.0 DIFT 10 1.6 1.6 5.5 
图 3 石灰岩动态强度增强因子随应变率的变化曲线Figure 3. Changes of the dynamic increase factors with strain rate for limestone1.4 状态方程参数的标定
p-α 状态方程是由Herrmann[28]提出的一种多孔材料的状态方程（equation of state, EOS），能够较好地描述孔隙压实过程中的材料性质变化，因此本文中的石灰岩采用该状态方程。p-α 状态方程根据材料孔隙压实过程分为3个阶段。第一阶段为压力小于孔隙坍塌压力pcrush时，采用线弹性状态方程：
p=Kμ (9) 式中：K为材料的体积模量，μ 为体积应变。
第二阶段为孔隙压实的过程，阶段的起始分别由孔隙坍塌压力pcrush和完全压实压力plock控制：
p=f(μ/α)/α (10) α=1 + (αp−1)(plock−pplock−pcrush)n (11) 式中：α 为材料的孔隙率，αp为对应于pcrush的材料孔隙率；f(μ)为材料在第三阶段密实段所采用的方程；n为压缩系数。第三阶段材料完全密实，可采用多项式状态方程描述：
p={A1μ+A2μ2+A3μ3μ≥0T1μμ＜0 (12) 式中：A1、A2、A3为压缩段参数，一般取A1 =K；T1为拉伸段参数，拉伸时按照线弹性考虑，即T1 =K。
对于线弹性段和孔隙压实段，可根据石灰岩的静水压力实验结果[19, 29]标定，而对于高压密实段主要通过平板撞击实验数据[30-31]进行参数标定，拟合得到的状态方程参数及曲线分别如表3和图4所示。
表 3 石灰岩的状态方程参数Table 3. Equation of state parameters for limestonepcrush/MPa plock/GPa n A1/GPa A2/GPa A3/GPa 100 1.44 3 22.53 –175.0 495.0 
图 4 石灰岩状态方程曲线与实验数据的拟合Figure 4. Equation of state of limestone fitted to experimental data1.5 模型验证
为了验证上述材料模型及参数的准确性，对Frew等[32]开展的4340钢弹体大范围打击速度侵彻石灰岩靶体的实验进行数值模拟，并对结果进行比较。实验中弹体直径分别为7.1和12.7 mm，长径比均为10，靶体尺寸根据弹体初始速度vs的不同而变化。基于前期经验[17]，数值模拟采用有限元与SPH相结合的方法，其中弹体采用有限元，靶体采用SPH粒子，并利用二维轴对称算法以提高计算效率。实验后弹体变形较小，因此视为刚性材料，石灰岩靶体由Kong-Fang模型描述且参数按1.4节中标定的选取。
数值模拟得到的不同直径弹体侵彻石灰岩靶体的侵彻深度（hp）与实验数据的对比如图5所示，图中展示了弹体直径分别为7.1和12.7 mm工况的结果对比。可以看出，数值模拟得到的侵彻深度与实验数据较接近，验证了石灰岩材料模型及其参数的可靠性。
图 5 数值模拟的侵彻深度与实验数据的对比Figure 5. Simulated depths of penetration at different initial projectile velocities compared with experimental data2. 弹体超高速侵彻石灰岩靶体的地冲击
当弹体撞击到靶体表面时，弹靶交界面处会产生瞬时高压力[33]。一方面，弹体继续向前运动，隧道区逐渐扩展，侵彻深度不断增大，且弹体的速度也逐渐降低，并可能发生变形、侵蚀等破坏现象；另一方面，在弹靶交界处产生的瞬时高压会在靶体中以应力波的形式传播，并对内部靶体产生不同程度的影响[9]，如图6所示。应力波在靶体中衰减、反射并发生相互作用，会造成石灰岩材料的压实、破碎或断裂，并可能进一步在靶体中形成开裂、开坑及震塌等破坏现象。
当弹体侵彻速度较低时，弹体对靶体的破坏主要表现为直接侵彻（隧道区）；而当弹体侵彻速度较高时，弹体会发生较严重的磨蚀甚至完全侵蚀，直接侵彻隧道区深度不会显著增大，而较高的侵彻速度会导致弹靶交界面处的压力升高，因此应力波对靶体的破坏效应凸显。基于上述验证的数值模拟方法及材料模型参数，本节对钨合金弹体超高速侵彻石灰岩靶体的地冲击进行研究。
2.1 数值模型
采用AUTODYN-2D轴对称算法，其中圆柱形钨合金弹体长2 m，直径为0.2 m，初速度区间考虑为1.5～5.0 km/s，采用有限元方法建模；石灰岩靶体长20 m，直径为10 m，采用SPH方法建模。弹体有限元网格大小按前期经验确定为10 mm×10 mm，同时对靶体的粒子间距大小进行了敏感性分析，分别开展了靶体粒子间距为50、40和30 mm的数值模拟，发现数值预测的侵彻深度随粒子间距减小趋于稳定。因此，同时考虑计算的准确性和计算效率，粒子间距确定为40 mm。在弹体着靶点正下方6、7、8、9、10、12、14、16、18和20 m处的靶体内部分别布置压力测点。
钨合金弹体采用Shock状态方程结合Steinberg强度模型，能够较好地体现钨合金材料在高速侵彻下的高应变率高压力学性能[2]，参数取AUTODYN中的默认值[34]。为了较好地描述弹体的变形及侵蚀，采用不同的单元删除算法对初速度为5.0 km/s的工况进行了收敛性论证。为了尽可能减少单元删除的影响，且综合考虑弹体的破坏模式、计算效率等因素，采用基于应变的单元删除算法，并确定单元删除阈值为0.8。石灰岩靶体采用1.1节中介绍的Kong-Fang模型及所标定的材料参数。
从模拟得到的弹体破坏程度来看，当弹体初速度低于2.5 km/s时，侵彻完成时弹体尚有残余长度；而当弹体速度高于3.0 km/s时，在侵彻过程中弹体被完全侵蚀，之后隧道区底部的靶体材料由于惯性会继续运动，造成侵彻深度进一步增大。不同撞击速度下模拟得到的最终侵彻深度如图7所示，其中黑色虚线为根据模拟结果拟合得到的趋势线。可以看出，在研究的速度区间内，侵彻可大致分为侵蚀侵彻阶段和流体侵彻阶段，侵彻深度随弹体初速度的变化整体呈先减小再增大的趋势，并最终趋于稳定，与前期的研究结论[3, 6]相符。并且最终的稳定侵彻深度大于流体动力学极限，这是由靶体材料的惯性运动造成的。
图 7 数值模拟得到不同撞击速度下的的侵彻深度Figure 7. Simulated depths of penetration at different initial projectile velocities2.2 超高速侵彻的应力波传播
以弹体初速度为5.0 km/s的工况为例，对弹体超高速侵彻石灰岩靶体的过程进行分析，其破坏现象和压力波传播过程如图8所示。弹体撞击靶体时，弹靶交界面产生瞬时压力（高达30 GPa）并开始向周围传播，侵彻过程中弹靶接触面的高压力状态会持续以应力波的形式传播，如图8(a)所示；1.4 ms时，弹体已被完全磨蚀，直接侵彻深度约为4.45 m，如图8(b)所示；弹体完全磨蚀后，坑底靶体材料在惯性作用下继续运动，弹坑继续扩张，在3.0 ms时，侵彻深度约为4.81 m；最终计算到100.0 ms时，侵彻深度为6.20 m且趋于稳定。从图8(c)～(d)可以看出，压力波基本呈半球形向靶体内传播，与地面浅埋爆炸具有一定的相似性；但由于侵彻过程具有方向性，沿深度方向的压力波峰值大于水平方向同等距离的压力波峰值。图8(d)中，当主要的地冲击已经传播到靶体内部较深处时，靶体靠近弹坑处仍然有压力分布，这是由于弹靶之间的相互作用会持续在弹靶交界面处产生新的高压状态，并持续向靶体中传播。
图 8 超高速侵彻过程中的破坏现象及压力波传播Figure 8. Damage and pressure wave propagation in the target during hypervelocity penetration其他弹体初速度下的压力波传播过程与上述算例基本一致。超高速侵彻造成的地冲击：一方面，会在隧道区周围造成高压区域，导致材料发生压实破碎甚至相变；另一方面，可能传入地下防护结构顶板，在顶板自由表面反射产生拉伸波，造成开裂、震塌等破坏效应。
2.3 超高速侵彻的地冲击分析
图9给出了不同侵彻速度下靶体中各测点沿深度方向的径向应力时程变化。可以看出不同侵彻速度下的应力波波形存在一定差异：当侵彻速度较低时，其升压时间较长，且可能出现压力进入平台段或开始下降、后再次升高的情况，这是由于侵彻完成时弹体尚有残留，会持续在弹靶接触面产生压力并向靶体内部传播；而当侵彻速度较高时，其应力迅速达到峰值并降低，这可能是由于弹体磨蚀较快，地冲击主要是由于初始侵彻阶段引起的瞬时高压造成的。当侵彻速度较高时，其升压时间较短，特别是弹体初速度为5.0 km/s时，其升压时间均在0.5 ms以内。需要说明的是，在弹体侵彻过程中，弹靶接触面上的应力峰值较大。对于布置深度小于或略大于侵彻深度的压力测点，会明显受到弹体的干扰，因此在总结地冲击规律时只选取了地冲击传播稳定后的测点数据，且不同初速度下的直接侵彻深度不同，因而地冲击传播稳定后的测点数量也不同。
图10给出了不同侵彻速度下应力峰值随深度变化的曲线，可以看出：随着深度的增大，应力峰值均不断降低，应力波在传播过程中不断衰减；随着侵彻速度的升高，同一深度处的应力峰值也会随之升高；当初速度高于3.0 km/s时，不同深度处的应力峰值升高较快。总体来看，弹体初速度较低时，侵彻产生的地冲击较小，对靶体内部的影响也较小；而当初速度逐渐升高时，地冲击的影响范围逐渐加深。
图 10 不同弹体初速度下应力峰值随深度的变化趋势Figure 10. Change of peak stress with depth at different initial projectile velocities3. 地冲击的影响因素
对于中、低速侵彻，弹靶材料、尺寸等因素对最终的破坏效应有显著的影响。对于超高速侵彻，已有研究表明：弹靶的强度不再起主导作用，破坏效应主要与弹体长径比（L/d）、弹靶密度比等因素有关[8]。而对应力波的传播，靶体材料的孔隙率会有较大影响[35]。因此，本节以弹体初速度5.0 km/s工况为基准，探究弹体长径比、弹靶密度比和靶体孔隙率等因素对地冲击的影响。
3.1 弹体长径比
弹体长径比，特别是弹体长度是影响侵彻深度的重要因素，本节采用改变弹体长度、保持弹体直径不变的方式探讨弹体长径比对侵彻深度的影响。数值模拟得到的侵彻深度如图11所示，可以看出，在弹径不变的情况下，随着弹体长度的增加，侵彻深度也不断增大。
图 11 侵彻深度与弹体长径比的关系（弹径不变）Figure 11. Depth of penetration versus length-to-diameter ratio at a constant projectile diameter图12(a)给出了不同长径比下应力峰值随深度的变化趋势，显然，随着侵彻深度的增大，地冲击影响的范围也会相应增大。因此引入相对深度（深度h与最终侵彻深度hp之比）作为评判地冲击影响范围的依据，如图12(b)所示。可以看出，当横坐标为相对深度时，不同长径比地冲击应力峰值变化趋势和数值基本一致。综合来看，在弹径不变的情况下，长径比对侵彻深度影响较大，而对地冲击规律影响不大。
3.2 弹体密度
在超高速侵彻阶段，弹靶密度比也是侵彻深度的重要影响因素之一。考虑到地质类材料的密度一般较为接近，因此只改变弹体的材料密度以考虑弹靶密度比对地冲击的影响。本节中保持弹体的尺寸不变，弹体材料选择密度差异较大的铝、钢作为对照，其计算条件及结果如表4所示，表中ρs为弹体密度，ρt为靶体材料密度。
表 4 不同弹体密度情况下的侵彻深度Table 4. Depths of penetration at different projectile densities材料 ρs/（kg∙m−3） （ρs/ρt）1/2 hp/m 铝 2785 1.10 2.58 钢 7830 1.85 4.69 钨 17000 2.73 6.20 3种弹体材料下靶体中不同深度的应力峰值变化趋势如图13所示。从图13(a)可以看出，弹体材料密度越高，地冲击影响的深度也越大，实际中超高速动能武器主要采用钨合金等高密度金属材料；而从图13(b)来看，不同弹体材料下地冲击影响的相对深度差别不大。总体来看，不同的弹体密度对侵彻深度有较大影响，但对地冲击规律的影响不大。
3.3 靶体孔隙率
地质材料性质差异大，对弹体的侵彻过程及应力波的传播都存在一定程度的影响。考虑到地质材料性质差异主要由孔隙率引起，本节中仅改变材料的孔隙率，保持靶体强度等参数不变，探讨孔隙率对地冲击的影响。
原石灰岩材料的孔隙率约为15%，在数值模拟时采用p-α 状态方程考虑其孔隙压实过程的性质变化。基于现有的实验数据，采用多项式状态方程，拟合得到无孔隙对照组岩石材料的状态方程曲线，如图14所示。
图 14 不同孔隙率岩石材料的状态方程Figure 14. Equations of state for limestones with different porosities数值模拟得到的不同孔隙率靶体内不同深度的应力峰值变化趋势如图15所示。可以看出，不同孔隙率的应力峰值差异较大，且孔隙率越高，其应力峰值越低。这主要是由于材料的孔隙率对应力波的传播有较大的影响，孔隙率越高，应力波衰减越快[35]，因此孔隙率对地冲击影响较大，这与Antoun[10]的研究结论一致。进一步分析可以发现，已有对地冲击的实验和数值模拟研究中采用的岩石不同[10-11, 13]，得到的地冲击数值差异较大，可能是由于孔隙率不同导致的。
基于上述讨论可以看出，弹体的参数变化（弹体的长径比、密度等）主要影响其侵彻深度，并会影响地冲击的破坏深度；但若从相对深度衡量，弹体参数变化对地冲击规律的影响不大。而靶体孔隙率影响靶体中应力波的传播，因此对地冲击有较大的影响。
4. 结 论
为探究超高速动能武器的对地破坏效应及其影响因素，采用数值模拟方法对弹体超高速侵彻的地冲击进行了研究。利用现有石灰岩静动态力学性能实验数据，对Kong-Fang材料模型进行了参数标定，并验证了材料模型及其参数的准确性。开展了钨合金弹体超高速侵彻石灰岩靶体的数值模拟，对不同速度侵彻下的靶体破坏过程和破坏机理进行了分析，得到以下结论。
（1）弹体超高速侵彻靶体时，弹靶接触时会产生瞬时高压，并以应力波的形式在靶体中传播，造成靶体内部不同程度的破坏。
（2）从靶体内部不同测点的地冲击应力时程曲线来看，当侵彻速度较低时，其升压和降压时间均较长；而当侵彻速度较高时，升压时间较短，可近似视为冲击波。
（3）从地冲击沿深度方向的应力峰值来看，随着深度的增大，应力峰值不断降低，地冲击不断衰减；随着侵彻速度的提高，同一深度处的应力峰值升高，地冲击显著增强；本文中的算例表明，当弹体侵彻速度高于3.0 km/s时，地冲击的影响较大。
（4）从相对深度来看，弹体参数（弹体长径比、密度）对地冲击规律的影响不大；靶体的特征特别是孔隙率对地冲击规律具有较大的影响。
-
图 1 最大强度面模型与实验数据的拟合
Figure 1. Fitting of the maximum strength surface model to the experimental data
图 2 石灰岩应力-应变曲线的数值模拟结果与实验数据的对比
Figure 2. Comparison of stress-strain curves of limestone between experimental data and numerical simulation
图 3 石灰岩动态强度增强因子随应变率的变化曲线
Figure 3. Changes of the dynamic increase factors with strain rate for limestone
图 4 石灰岩状态方程曲线与实验数据的拟合
Figure 4. Equation of state of limestone fitted to experimental data
图 5 数值模拟的侵彻深度与实验数据的对比
Figure 5. Simulated depths of penetration at different initial projectile velocities compared with experimental data
图 7 数值模拟得到不同撞击速度下的的侵彻深度
Figure 7. Simulated depths of penetration at different initial projectile velocities
图 8 超高速侵彻过程中的破坏现象及压力波传播
Figure 8. Damage and pressure wave propagation in the target during hypervelocity penetration
图 10 不同弹体初速度下应力峰值随深度的变化趋势
Figure 10. Change of peak stress with depth at different initial projectile velocities
图 11 侵彻深度与弹体长径比的关系（弹径不变）
Figure 11. Depth of penetration versus length-to-diameter ratio at a constant projectile diameter
图 14 不同孔隙率岩石材料的状态方程
Figure 14. Equations of state for limestones with different porosities
表 1 石灰岩的强度模型参数
Table 1. Parameters of the strength surface models for limestone
最大强度面 残余强度面 屈服强度面 损伤参数 a1 a2 a3 a1y a2y λm 0.877 0.022 0.80 0.890 0.046 3×10−5 
下载: 导出CSV
表 2 石灰岩的动态强度增强因子参数
Table 2. Parameters for dynamic increase factors of limestone
动态强度增强因子 Fm Wx S Wy DIFC 9 1.8 1.2 5.0 DIFT 10 1.6 1.6 5.5
下载: 导出CSV
表 3 石灰岩的状态方程参数
Table 3. Equation of state parameters for limestone
pcrush/MPa plock/GPa n A1/GPa A2/GPa A3/GPa 100 1.44 3 22.53 –175.0 495.0
下载: 导出CSV
表 4 不同弹体密度情况下的侵彻深度
Table 4. Depths of penetration at different projectile densities
材料 ρs/（kg∙m−3） （ρs/ρt）1/2 hp/m 铝 2785 1.10 2.58 钢 7830 1.85 4.69 钨 17000 2.73 6.20
下载: 导出CSV
-
[1] KONG X Z, WU H, FANG Q, et al. Projectile penetration into mortar targets with a broad range of striking velocities: test and analyses [J]. International Journal of Impact Engineering, 2017, 106: 18–29. DOI: 10.1016/j.ijimpeng.2017.02.022. [2] 钱秉文, 周刚, 李进, 等. 钨合金柱形弹超高速撞击水泥砂浆靶的侵彻深度研究 [J]. 爆炸与冲击, 2019, 39(8): 083301. DOI: 10.11883/bzycj-2019-0141.QIAN B W, ZHOU G, LI J, et al. Penetration depth of hypervelocity tungsten alloy projectile penetrating concrete target [J]. Explosion and Shock Waves, 2019, 39(8): 083301. DOI: 10.11883/bzycj-2019-0141. [3] 王杰, 武海军, 杨荷, 等. 高速/超高速侵彻半无限靶研究进展 [J]. 兵工学报, 2017, 38(S1): 73–88.WANG J, WU H J, YANG H, et al. Research progress in penetration of projectiles into semi-infinite targets at high-velocity/hypervelocity [J]. Acta Armamentarii, 2017, 38(S1): 73–88. [4] 李争, 刘元雪, 胡明, 等. “上帝之杖”天基动能武器毁伤效应评估 [J]. 振动与冲击, 2016, 35(18): 159–164. DOI: 10.13465/j.cnki.jvs.2016.14.026.LI Z, LIU Y X, HU M, et al. Damage effect evaluation of God stick space-based kinetic energy weapons [J]. Journal of Vibration and Shock, 2016, 35(18): 159–164. DOI: 10.13465/j.cnki.jvs.2016.14.026. [5] POELCHAU M H, KENKMANN T, THOMA K, et al. The MEMIN research unit: scaling impact cratering experiments in porous sandstones [J]. Meteoritics and Planetary Science, 2013, 48(1): 8–22. DOI: 10.1111/maps.12016. [6] 李干, 宋春明, 邱艳宇, 等. 超高速弹对花岗岩侵彻深度逆减现象的理论与实验研究 [J]. 岩石力学与工程学报, 2018, 37(1): 60–66. DOI: 10.13722/j.cnki.jrme.2017.0584.LI G, SONG C M, QIU Y Y, et al. Theoretical and experimental studies on the phenomenon of reduction in penetration depth of hyper-velocity projectiles into granite [J]. Chinese Journal of Rock Mechanics and Engineering, 2018, 37(1): 60–66. DOI: 10.13722/j.cnki.jrme.2017.0584. [7] 钱秉文, 周刚, 李进, 等. 钨合金弹体超高速撞击混凝土靶成坑特性研究 [J]. 北京理工大学学报, 2018, 38(10): 1012–1017. DOI: 10.15918/j.tbit1001-0645.2018.10.004.QIAN B W, ZHOU G, LI J, et al. Study of the crater produced by hypervelocity tungsten alloy projectile into concrete target [J]. Transactions of Beijing Institute of Technology, 2018, 38(10): 1012–1017. DOI: 10.15918/j.tbit1001-0645.2018.10.004. [8] 孔祥振, 方秦, 吴昊, 等. 长杆弹超高速侵彻半无限靶理论模型的对比分析与讨论 [J]. 振动与冲击, 2017, 36(20): 37–43. DOI: 10.13465/j.cnki.jvs.2017.20.007.KONG X Z, FANG Q, WU H, et al. Comparisons of long rod high velocity penetration models for semi-infinite targets [J]. Journal of Vibration and Shock, 2017, 36(20): 37–43. DOI: 10.13465/j.cnki.jvs.2017.20.007. [9] KONG X Z, WU H, FANG Q, et al. Rigid and eroding projectile penetration into concrete targets based on an extended dynamic cavity expansion model [J]. International Journal of Impact Engineering, 2017, 100: 13–22. DOI: 10.1016/ j.ijimpeng.2016.10.005. [10] ANTOUN T H, GLENN L A, WALTON O R, et al. Simulation of hypervelocity penetration in limestone [J]. International Journal of Impact Engineering, 2006, 33: 45–52. DOI: 10.1016/j.ijimpeng.2006.09.009. [11] 邓国强, 杨秀敏. 超高速武器对地打击效应数值仿真 [J]. 科技导报, 2015, 33(16): 65–71. DOI: 10.3981/j.issn.1000-7857.2015.16.010.DENG G Q, YANG X M. Numerical simulation of damage effect of hyper velocity weapon on ground target [J]. Science and Technology Review, 2015, 33(16): 65–71. DOI: 10.3981/j.issn.1000-7857.2015.16.010. [12] 邓国强, 杨秀敏. 超高速武器流体侵彻与装药浅埋爆炸效应的等效方法 [J]. 防护工程, 2015, 37(6): 27–32.DENG G Q, YANG X M. Effect equivalent method between fluid penetration of hypervelocity weapon and shallow detonation of explosive [J]. Protective Engineering, 2015, 37(6): 27–32. [13] 王明洋, 岳松林, 李海波, 等. 超高速弹撞击岩石的地冲击效应等效计算 [J]. 岩石力学与工程学报, 2018, 37(12): 2655–2663. DOI: 10.13722/j.cnki.jrme.2018.0473.WANG M Y, YUE S L, LI H B, et al. An equivalent calculation method of ground shock effects of hypervelocity projectile striking on rock [J]. Chinese Journal of Rock Mechanics and Engineering, 2018, 37(12): 2655–2663. DOI: 10.13722/j.cnki.jrme.2018.0473. [14] 方秦, 孔祥振, 吴昊, 等. 岩石Holmquist-Johnson-Cook模型参数的确定方法 [J]. 工程力学, 2014, 31(3): 197–204. DOI: 10.6052/j.issn.1000-4750.2012.10.0780.FANG Q, KONG X Z, WU H, et al. Determination of Holmquist-Johnson-Cook consitiutive model parameters of rock [J]. Engineering Mechanics, 2014, 31(3): 197–204. DOI: 10.6052/j.issn.1000-4750.2012.10.0780. [15] FOSSUM A F, SENSENY P E, PFEIFLE T W, et al. Experimental determination of probability distributions for parameters of a Salem limestone cap plasticity model [J]. Mechanics of Materials, 1995, 21(2): 119–137. DOI: 10.1016/0167-6636(95)00002-X. [16] KONG X Z, FANG Q, CHEN L, et al. A new material model for concrete subjected to intense dynamic loadings [J]. International Journal of Impact Engineering, 2018, 120: 60–78. DOI: 10.1016/j.ijimpeng.2018.05.006. [17] ZHANG S B, KONG X Z, FANG Q, et al. Numerical prediction of dynamic failure in concrete targets subjected to projectile impact by a modified Kong-Fang material model [J]. International Journal of Impact Engineering, 2020, 144: 103633. DOI: 10.1016/j.ijimpeng.2020.103633. [18] XU H, WEN H M. Semi-empirical equations for the dynamic strength enhancement of concrete-like materials [J]. International Journal of Impact Engineering, 2013, 60: 76–81. DOI: 10.1016/j.ijimpeng.2013.04.005. [19] FOSSUM A F, BRANNON R M. On a viscoplastic model for rocks with mechanism-dependent characteristic times [J]. Acta Geotechnica, 2006, 1(2): 89–106. DOI: 10.1007/s11440-006-0010-z. [20] FOSSUM A F. Rock penetration: finite element sensitivity and probabilistic modeling analyses [R]. Albuquerque, New Mexico: Sandia National Laboratories, 2004: 13–31. [21] WALTON G, HEDAYAT A, KIM E, et al. Post-yield strength and dilatancy evolution across the brittle-ductile transition in Indiana limestone [J]. Rock Mechanics and Rock Engineering, 2017, 50(7): 1691–1710. DOI: 10.1007/s00603-017-1195-1. [22] FREW D J, FORRESTAL M J, CHEN W. A split Hopkinson pressure bar technique to determine compressive stress-strain data for rock materials [J]. Experimental Mechanics, 2001, 41(1): 40–46. DOI: 10.1007/BF02323102. [23] CHAKRABORTY T, MISHRA S, LOUKUS J, et al. Characterization of three Himalayan rocks using a split Hopkinson pressure bar [J]. International Journal of Rock Mechanics and Mining Sciences, 2016, 85: 112–118. DOI: 10.1016/ j.ijrmms.2016.03.005. [24] KHAN A S, IRANI F K. An experimental study of stress wave transmission at a metallic-rock interface and dynamic tensile failure of sandstone, limestone, and granite [J]. Mechanics of Materials, 1987, 6(4): 285–292. DOI: 10.1016/0167-6636(87)90027-5. [25] KUBOTA S, OGATA Y, WADA Y, et al. Estimation of dynamic tensile strength of sandstone [J]. International Journal of Rock Mechanics and Mining Sciences, 2008, 45(3): 397–406. DOI: 10.1016/j.ijrmms.2007.07.003. [26] WANG Q Z, LI W, XIE H P. Dynamic split tensile test of flattened Brazilian disc of rock with SHPB setup [J]. Mechanics of Materials, 2009, 41(3): 252–260. DOI: 10.1016/j.mechmat.2008.10.004. [27] CHO S H, OGATA Y, KANEKO K. Strain-rate dependency of the dynamic tensile strength of rock [J]. International Journal of Rock Mechanics and Mining Sciences, 2003, 40(5): 763–777. DOI: 10.1016/S1365-1609(03)00072-8. [28] HERRMANN W. Constitutive equation for the dynamic compaction of ductile porous materials [J]. Journal of Applied Physics, 1969, 40(6): 2490–2499. DOI: 10.1063/1.1658021. [29] HEARD H C, ABEY A E, BONNER B P. High pressure mechanical properties of Indiana limestone: UCID-16501 [R]. USA: California University, Lawrence Livermore National Laboratory, 1974. [30] LARSON D B, ANDERSON G D. Plane shock wave studies of porous geologic media [J]. Journal of Geophysical Research: Solid Earth, 1979, 84(B9): 4592–4600. DOI: 10.1029/JB084iB09p04592. [31] MURRI W J, SMITH C W, MAHRER K D. Equation of state of rocks: PYU-1883 [R]. USA: Stanford Research Institute, 1974. [32] FREW D J, FORRESTAL M J, HANCHAK S J. Penetration experiments with limestone targets and ogive-nose steel projectiles [J]. Journal of Applied Mechanics, 2000, 67(4): 841–845. DOI: 10.1115/1.1331283. [33] MCFARLAND C, PAPADOS P, GILTRUD M. Hypervelocity impact penetration mechanics [J]. International Journal of Impact Engineering, 2008, 35(12): 1654–1660. DOI: 10.1016/j.ijimpeng.2008.07.080. [34] STEINBERG D J, COCHRAN S G, GUINAN M W. A constitutive model for metals applicable at high strain rate [J]. Journal of Applied Mechanics, 1980, 51(3): 1498–1504. DOI: 10.1063/1.327799. [35] 杜成成. 应力波在混凝土中传播特性及结构特征参数监测研究 [D]. 哈尔滨: 哈尔滨工业大学, 2018: 49–77.DU C C. Research on transmission properties of stress waves in concrete and monitoring of structural characteristic parameters [D]. Harbin: Harbin Institute of Technology, 2018: 49–77. 期刊类型引用(3)
1. 张哲伟，尚宏春，赵文，燕超，杨臻. 高应变率下金属动态行为表征和断裂研究. 兵器材料科学与工程. 2024(06): 40-46 . 
百度学术2. 潘鹏志，王兆丰，封雨捷，李雨芯. 深空探测任务中岩土力学若干问题及研究进展. 岩土力学. 2024(11): 3153-3172 . 百度学术3. 李国辉，刘基程，耿汉生，马林建，李增. 岩石高速-超高速侵彻效应研究进展. 防护工程. 2023(05): 67-78 . 百度学术其他类型引用(1)
-
下载:
百度学术
计量
- 文章访问数: 493
- HTML全文浏览量: 615
- PDF下载量: 226
- 被引次数: 4