Research on the Ballistic Performance of Cement Mortar
-
摘要: 对于水泥砂浆的抗弹性能研究,目前很少考虑靶体所处的应力状态,为此基于自研的真三轴静载混凝土侵彻实验装置和水泥砂浆抗弹性能实验结果,讨论了水泥砂浆在不同应力状态下的开坑深度和开坑阻力。应用侵彻深度的经验公式和基于HJC模型的有限元数值计算方法,对比分析了水泥砂浆侵彻实验,结果表明,对于低速冲击过程,采用UMIST公式和HJC模型的数值分析对开坑深度的预测较为有效。应力状态对开坑深度有明显的影响,即随着侧限增加,水泥砂浆的三轴强度提高,弹丸的开坑深度减小。应用基于HJC模型的数值分析方法,研究了弹丸开坑过程中弹体内的加速度波形和y轴支撑杆上的波形,结果表明:弹丸开坑过程对两种波形都有影响,其中y轴支撑杆上的波形可以更好地反映开坑过程。虽然数值模拟结果与实验波形的趋势基本一致,但是应力幅值有一定的差异,说明基于HJC模型的数值分析对开坑阻力的计算能力尚待提高。Abstract: The stress state of the target was seldom taken into account in the investigation of the ballistic performance of cement mortar. Based on the self-developed penetration experimental system of concrete under true-triaxial confinement and the experimental results of the anti-bullet performance of cement mortar, the depth and resistance of opening pit under different stress states were discussed in the present paper. The empirical formula of penetration depth and finite element method (FEM) based on HJC model were used to analyze the penetration behaviors of cement mortar results. The results showed that under the lower velocity impact, the UMIST formula and the HJC model were both effective in the prediction of pit depth. At the same time, the stress state had an obvious influence on pit depth. With the increase of the lateral limit, the cement mortar strength increases and the pit depth of the projectile decreases. The acceleration wave in bullet and the wave in the y-axis support rod were calculated by FEM based on HJC model. The results showed that the process of projectile opening pit would be recorded by these two waveforms, and the wave structure in the y-axis support rod would be more significantly. Although the tendency of the simulation results was basically consistent with the experimental waves, there was difference in the stress amplitude to some degree, which also indicated that the calculation method of the pit opening resistance based on HJC model needed to be improved.
-
Key words:
- ballistic performance /
- pit depth /
- pit resistance /
- cement mortar /
- true-triaxial stress state
-
负泊松比蜂窝结构又称拉胀结构,因具有许多常规结构不具备的力学特性[1],而成为研究热点。蜂窝材料具有较高的相对刚度、强度和高效的能量吸收能力,在抗剪切、抗屈曲、提高硬度以及抗疲劳等方面拥有独特的优越性[2-3],在一些应用领域中发挥着关键作用,如汽车、航空、军事、医学领域[4]。多孔结构的力学性能主要取决于细观上的拓扑结构。近年来,通过改变细观结构,人们发现负泊松比结构具有很多特殊优势,因而被广泛应用[5]。
马芳武等[6]研究了一种内凹三角形负泊松比结构,通过改变内凹角度,分析了冲击端和固定端的平台应力和能量吸收能力,并与内凹六边形进行了对比。Zhang等[7]分析了内凹六边形蜂窝在两个正交方向上的后继屈服拉伸行为,同时考虑结构的塑性影响以及孔壁的非线性行为分析模型,提出了单胞结构的塑性铰变形机制,得到了单胞结构的应力-应变曲线。Li等[8-9]对内凹蜂窝结构进行分级、强化,并将正弦曲线引入内凹蜂窝结构,得到了新的改进模型,进而分析了结构的泊松比和能量吸收变化。邓小林等[10]研究了全参数化的正弦曲线蜂窝结构,以不同振幅、不同厚度建立模型,研究了蜂窝结构在不同冲击速度下的动力响应,发现正弦曲线蜂窝较常规六边形蜂窝有更好的能量吸收效果。崔世堂等[11]利用有限元模拟方法研究了负泊松比蜂窝结构面内冲击动力学特性,发现平台应力和结构的比吸能随冲击速度的增大而增高,随胞元扩展角的增大而降低。陈鹏等[12]研究了具有零泊松比特征的半凹角蜂窝结构,并将其与正六边形蜂窝和内凹负泊松比蜂窝在面内冲击荷载作用下的抗冲击性能进行对比分析,数值结果表明,半凹角蜂窝的抗冲击性能介于正六边形蜂窝和内凹蜂窝之间。Hu等[13]通过理论分析和数值模拟,研究了内凹角度和壁长对内凹负泊松比蜂窝在大变形下的单轴动态冲击性能的影响,推导出冲击过程中平均冲击应力的经验公式。Zhang等[14]通过有限元模拟,研究了内凹蜂窝x方向的平面内动态冲击行为,发现内凹蜂窝的面内动态性能不仅与冲击速度和边缘厚度有关,还受蜂窝壁角的影响。Li等[15]通过单轴和双轴压缩模拟以及理论分析,研究了正六边形蜂窝结构的面内压缩动态力学性能,分析了双轴压缩的变形模式,结果表明:相比单轴冲击,双轴冲击下在x和y方向的真实应力增强,能量吸收能力也得到了提高,且完全致密化阶段比单轴压缩阶段更平滑。此外,Li等[16]研究了六边形、内凹、混合3种蜂窝模型在单、双轴冲击下的面内动态力学性能,结果表明:正交双轴冲击下,六边形蜂窝表现出3种变形模式,内凹和混合型蜂窝没有明显的过渡模式,由于负泊松比效应的影响,内凹蜂窝具有较差的耗能能力。
值得注意的是,自然界中的蜂窝结构和人造蜂窝结构在细观上总存在一定的缺陷,从而引起结构的不规则性,力学性能也会发生一定的变化。Ajdari等[17]通过数值模拟研究了正六边形和不规则二维蜂窝的平面内动态冲击问题,分析了孔壁缺失和空间扰动形成的结构微观不规则性对力学性能的影响。Alkhader等[18]用函数定义六边形蜂窝、随机Voronoi泡沫以及正方形和三角形拓扑结构等多种二维拓扑结构的不规则程度,以研究其单轴压缩响应,结果表明,相对于以弯曲为主的结构,以拉伸为主的结构有表现出灾难性屈服后软化反应的趋势,而不规则性则会导致更多的弯曲现象。Liu等[19]对内凹蜂窝材料的面内动态冲击过程进行了数值模拟,并在此基础上定义了内凹蜂窝结构的不规则性,分析发现,在准静态下不规则的内凹蜂窝比规则的正六边形蜂窝能吸收更多的能量,但这种情况在高速撞击下逆转。Zheng等[20]通过数值模拟研究了坐标扰动和Voronoi随机模型两种不规则模型与正六边形蜂窝在不同冲击速度下的变形模式和平台冲击力,得到不规则性结构更具复杂性的结论。Zhu等[21]研究了孔的不规则性对二维随机泡沫弹性性能的影响,构造了不规则度不同的周期性随机结构,并通过数值模拟确定了其有效弹性性能,结果表明,二维随机泡沫体形状越不规则,有效弹性模量和剪切模量越大,在一定的压比相对密度下,体积模量越小。
综上所述,实际中蜂窝结构往往是不规则的,且易受双轴冲击载荷作用。而关于不规则结构在双轴冲击下的研究较少,为此本工作将针对不规则内凹负泊松比蜂窝结构在双轴冲击下的面内冲击响应,分析规则度和冲击速度对结构变形影响的规律。
1. 计算模型
1.1 不规则内凹蜂窝结构有限元模型的建立
采用如图1所示的节点扰动方法来建立不规则内凹蜂窝的有限元模型。
如图1所示,将规则的内凹六边形蜂窝结构的每个节点按照式(1)中的方法进行随机扰动
{yi=y0+Δysinθmxi=x0+Δxcosθm (1) 式中:
$\theta _{\rm m}$ 为角度随机值(${0^\circ } \leqslant \theta_{\rm m} \leqslant {360^\circ }$ ),${x_0}$ 、${y_0}$ 为节点的原始坐标,${x_i}$ 、${y_i}$ 为扰动后的节点坐标。为了保证随机扰动之后模型的棱壁不会重叠,需要对$\Delta x$ 、$\Delta y$ 进行如式(2)的限制{μ=Δx2+Δy20⩽μ⩽μm (2) 式中:
$\;\mu $ 为节点扰动的随机长度;$\;\mu_{\rm{m}}$ 为扰动的最大长度,$0 \leqslant {\mu _{\rm{m}}} \leqslant {l_1}/2$ 。内凹负泊松比蜂窝结构的不规则度可以定义为K=2μml1 (3) 式中:
${l_1}$ 为规则蜂窝结构的最短棱壁长度。假设蜂窝结构所有棱壁的厚度均相同,则可通过改变棱壁的厚度来调节蜂窝结构的相对密度。本研究采用15%的相对密度进行分析,图2显示了部分模型。
图2中内凹蜂窝结构的相对密度
$\Delta \rho $ 可以表示为Δρ=ρ∗ρs=N∑i=1li×tL1×L2 (4) 此外,对于规则的内凹负泊松比蜂窝,其相对密度
$\Delta\rho_{\rm r} $ 也可以表示为Δρr=12tl1(l2/l1+2)cosα(l2/l1+sinα) (5) 式中:
${\;\rho ^*}$ 为模型的密度,${\;\rho \rm_s}$ 为基体材料的密度,${l_i}$ 为各个孔壁的长度,$t$ 为孔壁的厚度,$N$ 为孔壁的数量,${L_1}$ 、${L_2}$ 为整个蜂窝结构的长度和宽度,l2为规则蜂窝结构的最长棱壁长度。1.2 有限元模型
采用ABAQUS/EXPLICIT软件进行分析。模型的边界条件设置:在两个正交方向上,将模型置于两块刚性板之间、底部刚性板之上,底部和左端的刚性板作为固定端, 顶部和右端作为冲击端,冲击速度为3~100 m/s,同时约束内凹蜂窝结构的面内自由度,如图3所示。建立的内凹蜂窝结构的主要参数为L1 = 129.9 mm, L2 = 120.0 mm, l1 = 5 mm, l2 = 10 mm,
$\theta $ = 60°。由于蜂窝铝具有高强度和高刚度的良好力学性能,本研究采用铝合金作为基体材料,主要参数为:密度$\;\rho$ = 2700 kg/m3,弹性模量E = 72 GPa,泊松比为0.33,屈服强度${\sigma _y}$ = 103 MPa,并采用线性强化模型,图4为结构基体材料的本构关系,其中Et为切线模量,$\sigma\rm{_s} $ 为线性强化模型的屈服强度。蜂窝细胞数量为15 × 15,可保证材料不受尺寸效应的影响。所有单元均采用4节点壳单元进行网格划分,网格单元尺寸为0.5 mm,节点数为28660,网格数为 19 540,建立无摩擦和通用接触。2. 结果与讨论
2.1 双轴冲击下的变形模式
为了对双轴冲击条件进行分类,采用与双轴冲击有关的参数
$\lambda $ ,表达式为$\lambda = {v_x}/{v_y}$ ,其中${v_x}$ 、${v_y}$ 分别为$x$ 、$y$ 方向的冲击速度。这里只讨论$\lambda = 1$ 的情况,显然$\lambda = 1$ 时为等双轴冲击。首先研究了不同规则度(K = 0, 0.6, 1.0)的内凹负泊松比在不同冲击速度(6、50和100 m/s)下的变形模式。图5、图6和图7给出了内凹蜂窝结构的结构变形情况。需要说明的是,为更好地展示变形结果,每隔约10%的应变截取一张变形模态图,同时为了清晰、规律地显示图像,所有图形都设置了相同的大小。从图5可以看出,对于规则的内凹蜂窝,在等双轴低速冲击过程中,结构首先在交叉处棱壁堆积,从而使内部先形成四边形,结构整体的变形在近端和远端都较均匀。这与文献[15]中内凹蜂窝的变形是一致的,也验证了本模型的有效性。持续的压缩使孔壁进一步堆积形成局部致密化,结构的致密过程主要是局部致密。由于负泊松比效应的影响,材料在一个方向受压时,其另一个正交方向会出现颈缩。因此,在双向冲击受压的情况下,结构会更早进入完全密实阶段。随着冲击速度的增大,结构从冲击端(上部和右端)开始密实,而固定端几乎没有变形。随着应变增加,致密向固定端传递,直至完全进入密实化。从图5中第2行和第3行图像还可以看出,随着冲击速度的增大,蜂窝结构的下端会产生部分“翘起”现象,这是由于负泊松比效应的影响会导致结构颈缩,且结构与固定端端部是无绑定约束,从而造成这类现象。
与规则蜂窝不同的是,不规则蜂窝结构在低速冲击下,其内部不会形成较为规则的四边形。这是由于不规则度的存在使结构棱壁处的堆叠也变得不规则。此外,从图5~图7中
$\varepsilon = 0.5$ 列可以看出,由于不规则度的引入,结构的变形模式由局部密实转变为整体密实,从而使内凹蜂窝结构在相同压缩程度下,密实化程度明显降低。在高速冲击下(v = 100 m/s),从图5~图7中可以看出,不规则程度越高,冲击端的致密程度越大。这是因为高速冲击下,结构在冲击端的密实主要是棱壁的弯曲折叠过程,随着不规则度的增加,棱壁的弯曲折叠受到的约束增大,向固定端传递的速度也会降低,所以不规则蜂窝结构的密实过程会更长,而在冲击端密实程度也会更高。此外,从图6、图7中$\varepsilon = 0.6$ 、v = 100 m/s对应的变形情况可以看出,固定端还有尚未变形进入密实的孔,说明不规则蜂窝结构具有较长的平台阶段,能够承受更大的压缩变形。2.2 不同冲击速度下的应力-应变曲线
图8和图9给出了蜂窝结构在双轴冲击下两个正交方向冲击端的名义应力-应变曲线,其中名义应力
$\sigma $ 通过冲击端的反力除以对应截面的原始面积获得,名义应变$\varepsilon$ 通过冲击位移除以对应的原长获得。从图中可以看出,内凹蜂窝结构在不同方向上的$\sigma \text{-}\, \varepsilon$ 曲线均表现出典型多孔材料在受压时所具有的弹性阶段、平台阶段和密实阶段3部分。从图8和图9中v = 6 m/s时的曲线可以看出,对于K = 0时的应力-应变曲线,在应变接近0.4处,结构变形的平台阶段均出现一个上升的阶梯,并且x方向最明显。结合2.1节关于变形模态的分析,认为这主要是由于在等低速双轴冲击下内凹蜂窝结构变形主要经历两种棱壁堆叠过程,即堆叠形成四边形以及四边形的进一步弯曲堆叠。由于第1步的堆叠,棱壁基本不会屈曲,主要是旋转折叠,因此这一平台阶段的应力水平较低;第2步的堆叠主要是棱壁的屈曲折叠,所以此阶段的应力水平较高。从图8和图9中也可以看出,K = 0时,结构会更早进入密实化阶段,而不规则度的引入使结构拥有较长的平台阶段,密实化阶段出现滞后现象,此现象与2.1节中变形模态的分析结果是一致的。随着冲击速度的增大,平台阶段的应力升高,说明结构的能量吸收能力随着冲击速度的增大而增强。
蜂窝结构的平台应力一般表示为
σp=1εd−ε0∫εdε0σ(ε)dε (6) 式中:
${\sigma \rm_p}$ 为平台应力;${\varepsilon _0}$ 为对应初始应力峰值的名义应变;${\varepsilon \rm_d}$ 为锁定应变,为蜂窝结构密实化阶段所对应的应变;$\sigma \left( \varepsilon \right)$ 为名义应力-应变曲线。图10给出了不规则度不同的内凹蜂窝结构在两个正交方向上不同冲击速度下的平台应力变化趋势。从图10可以看到:随着冲击速度的增大,两个方向上的平台应力值都会上升;对于K = 0的规则蜂窝结构,其在两个方向上的平台应力相差较大,这是结构的各向异性所导致的。引入不规则度时,在高速冲击下两个方向上的平台应力变化大小及趋势都较接近,说明结构的各向异性降低,这一点从2.1节内凹蜂窝结构的变形模态中也可以看出。
2.3 内凹蜂窝结构的能量吸收能力
在动态冲击过程中,能量主要由材料的塑性变形消耗。采用比塑性耗散能表征单位质量的能量吸收能力,表达式为
W=EPEDM (7) 式中:
${E\rm_{PED}}$ 为塑性耗散能,可以从有限元分析软件中直接获得;M为结构的质量。图11给出了内凹蜂窝结构在6、50和100 m/s 3种不同冲击速度下的比塑性能量耗散与
$y$ 方向冲击应变的关系。从图11中可以看出,当应变较低时,比塑性耗散能$W$ 上升较缓慢,且所有曲线基本重合。这表明在早期,不规则度对内凹蜂窝结构的影响较小。随着压缩程度的增加,$W $ 增加的速率变大,且K = 0时,$W $ 增加得最快,表明结构开始进入密实阶段,这是由结构的负泊松比效应引起的。对于不规则蜂窝结构,曲线上升得较缓慢,表明不规则度的引入使结构的平台阶段延长,结构具有更强的能量吸收能力。3. 结 论
采用有限元方法研究了具有不同不规则度内凹负泊松比结构的面内双轴冲击响应,得到了以下结论。
(1)内凹蜂窝结构的变形受冲击速度的影响。随着冲击速度的提高,蜂窝结构的变形逐渐转向逐层致密,受结构负泊松比效应的影响,在等高速双轴压缩时,结构的固定端会有局部“翘起”现象。此外,由于不规则度的引入,在低速冲击下,结构的密实化过程从局部致密转变为整体致密,从而导致在相同的压缩程度下,结构的密实化程度降低。
(2)随着冲击速度的增大,平台阶段的应力上升,能量吸收能力更强,比塑性耗散能也上升。不规则度的引入延长了平台阶段,降低了结构的各向异性程度,从而提高了结构的能量吸收能力。
-
表 1 水泥砂浆的HJC本构模型参数
Table 1. Parameters of HJC model for cement mortar
${\;\rho {_0} }$/(kg·m−3) $G$/GPa AHJC/GPa BHJC/GPa CHJC NHJC $f{'} $/MPa 1844 1.32 0.66 1.335 0.0018 0.845 14.4 T/MPa ${\dot \varepsilon{_0} }$/s−1 $\varepsilon $f,min Smax ${p{\rm{_c}} }$/MPa ${\;\mu {\rm{_c}} }$ ${p{\rm{_l} } }$/GPa 2.0 1 0.01 80.24 13.8 0.0075 1.096 ${\;\mu {\rm{_l} } }$ k1/GPa k2/GPa k3/GPa D1 D2 0.15 85 −171 208 0.006629 1.0 表 2 弹丸JC本构模型参数
Table 2. Parameters of JC model of projectile
${\rho{_0} }$/(kg·m−3) G/GPa T0/K c/(J·kg−1·K−1) AJC/MPa BJC/MPa nJC 7830 77 293 477 792 510 0.26 CJC Tm/K d1 d2 d3 d4 d5 0.014 1793 0.05 3.44 −2.12 0.002 0.61 表 3 无量纲侵彻深度公式参数
Table 3. Formula parameters of dimensionless penetration depth
Stress state k0 m1 m2 m3 No confinement 0.90 0.70 1.21 0.60 Unilateral confinement 1.02 0.55 1.25 0.65 Bilateral confinement 1.05 0.40 1.30 0.68 -
[1] 徐松林, 王鹏飞, 赵坚, 等. 基于三维Hopkinson杆的混凝土动态力学性能研究 [J]. 爆炸与冲击, 2017, 37(2): 180–185. doi: 10.11883/1001-1455(2017)02-0180-06XU S L, WANG P F, ZHAO J, et al. Dynamic behavior of concrete under static triaxial loading using 3D-Hopkinson bar [J]. Explosion and Shock Waves, 2017, 37(2): 180–185. doi: 10.11883/1001-1455(2017)02-0180-06 [2] 徐松林, 王鹏飞, 单俊芳, 等. 真三轴静载作用下混凝土的动态力学性能研究 [J]. 振动与冲击, 2018, 37(15): 59–67. doi: 10.13465/j.cnki.jvs.2018.15.008XU S L, WANG P F, SHAN J F, et al. Dynamic behavior of concrete under static tri-axial loadings [J]. Journal of Vibration and Shock, 2018, 37(15): 59–67. doi: 10.13465/j.cnki.jvs.2018.15.008 [3] XU S L, SHAN J F, ZHANG L, et al. Dynamic compression behaviors of concrete under true triaxial confinement: an experimental technique [J]. Mechanics of Materials, 2020, 140: 103220. doi: 10.1016/j.mechmat.2019.103220 [4] 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 [5] 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 [6] CHEN X W, LI J C. Analysis on the resistive force in penetration of a rigid projectile [J]. Defence Technology, 2014, 10(3): 285–293. doi: 10.1016/j.dt.2014.06.007 [7] 沈河涛. 弹丸侵彻混凝土介质效应的研究[D]. 北京: 北京理工大学, 1996.SHEN H T. Study on the effect of projectile penetrating concrete medium [D]. Beijing: Beijing Institute of Technology, 1996. [8] BACKMAN M E, GOLDSMITH W. The mechanics of penetration of projectiles into targets [J]. International Journal of Engineering Science, 1978, 16(1): 1–99. doi: 10.1016/0020-7225(78)90002-2 [9] 薛建锋, 沈培辉, 王晓鸣. 弹体侵彻混凝土开坑阶段阻力的计算 [J]. 高压物理学报, 2016, 30(6): 499–504. doi: 10.11858/gywlxb.2016.06.010XUE J F, SHEN P H, WANG X M. Resistance during cratering for projectile penetrating into concrete target [J]. Chinese Journal of High Pressure Physics, 2016, 30(6): 499–504. doi: 10.11858/gywlxb.2016.06.010 [10] 蒋志刚, 甄明, 刘飞, 等. 钢管约束混凝土抗侵彻机理的数值模拟 [J]. 振动与冲击, 2015, 34(11): 1–6. doi: 10.13465/j.cnki.jvs.2015.11.001JIANG Z G, ZHEN M, LIU F, et al. Simulation of anti-penetration mechanism of steel tube confined concrete [J]. Journal of Vibration and Shock, 2015, 34(11): 1–6. doi: 10.13465/j.cnki.jvs.2015.11.001 [11] 朱翔, 陆新征, 杜永峰, 等. 外包钢管加固RC柱抗冲击试验研究 [J]. 工程力学, 2016, 33(6): 23–33. doi: 10.6052/j.issn.1000-4750.2014.11.0991ZHU X, LU X Z, DU Y F, et al. Experimental study on impact resistance of reinforced conceret columns strengthened with steel jackets [J]. Engineering Mechanics, 2016, 33(6): 23–33. doi: 10.6052/j.issn.1000-4750.2014.11.0991 [12] 甄明, 蒋志刚, 万帆, 等. 钢管约束混凝土抗侵彻性能试验 [J]. 国防科技大学学报, 2015, 37(3): 121–127. doi: 10.11887/j.cn.201503020ZHEN M, JIANG Z G, WAN F, et al. Steeltube confined concrete targets penetration experiments [J]. Journal of National University of Defense Technology, 2015, 37(3): 121–127. doi: 10.11887/j.cn.201503020 [13] 蒙朝美, 宋殿义, 蒋志刚, 等. 多边形钢管约束混凝土靶抗侵彻性能试验研究 [J]. 振动与冲击, 2018, 37(13): 14–19. doi: 10.13465/j.cnki.jvs.2018.13.003MENG C M, SONG D Y, JIANG Z G, et al. Tests for anti-penetration performance of polygonal steel tube-confined concrete targets [J]. Journal of Vibration and Shock, 2018, 37(13): 14–19. doi: 10.13465/j.cnki.jvs.2018.13.003 [14] 徐松林, 单俊芳, 王鹏飞, 等. 三轴应力状态下混凝土的侵彻性能研究 [J]. 爆炸与冲击, 2019, 39(7): 071101. doi: 10.11883/bzycj-2019-0034XU S L, SHAN J F, WANG P F, et al. Penetration performance of concrete under triaxial stress [J]. Explosion and Shock Waves, 2019, 39(7): 071101. doi: 10.11883/bzycj-2019-0034 [15] 陈丽娜, 单俊芳, 周李姜, 等. 应力状态对水泥砂浆侵彻性能的影响 [J]. 振动与冲击, 2020, 39(15): 32–40. doi: 10.13465/j.cnki.jvs.2020.15.005CHEN L N, SHAN J F, ZHOU L J, et al. Effects of stress state on penetration performance of cement mortar [J]. Journal of Vibration and Shock, 2020, 39(15): 32–40. doi: 10.13465/j.cnki.jvs.2020.15.005 [16] MEYER C S. Development of geomaterial parameters for numerical simulations using the Holmquist-Johnson-Cook constitutive model for concrete: ARL-TR-5556 [R]. Orlando: Army Research Laboratory, 2011. [17] JOHNSON G R, COOK W H. Fracture characteristics of three metals subjected to various strains, strain rates, temperatures and pressures [J]. Engineering Fracture Mechanics, 1985, 21(1): 31–48. doi: 10.1016/0013-7944(85)90052-9 [18] Army Corps of Engineers. Fundamentals of protective design: AT1207821 [R]. Army Corps of Engineers, 1946. [19] National Defense Research Committee. Effects of impact and explosion: summery technical report of division 2 [R]. Washington DC: National Defense Research Committee, 1946. [20] KENNEDY R P. A review of procedures for the analysis and design of concrete structures to resist missile impact effects [J]. Nuclear Engineering and Design, 1976, 37(2): 183–203. doi: 10.1016/0029-5493(76)90015-7 [21] BARR P. Guidelines for the design and assessment of concrete structures subjected to impact [R]. London, UK: UK Atomic Energy Authority, Safety and Reliability Directorate, 1990. [22] YOUNG C W. Penetration equations: SAND 97-2426 [R]. Albuquerque, NM, US: Sandia National Laboratories, 1997. [23] REID S R, WEN H M. Predicting penetration, cone cracking, scabbing and perforation of reinforced concrete targets struck by flat-faced projectiles: UMIST Report ME/AM/02.01/TE/G/018507/Z [R]. Manchester: University of Manchester Institute of Science and Technology, 2001. [24] LI Q M, CHEN X W. Dimensionless formulae for penetration depth of concrete target impacted by a non-deformable projectile [J]. International Journal of Impact Engineering, 2003, 28(1): 93–116. doi: 10.1016/S0734-743X(02)00037-4 [25] FORRESTAL M J, FREW D J, HICKERSON J P, et al. Penetration of concrete targets with deceleration-time measurements [J]. International Journal of Impact Engineering, 2003, 28(5): 479–497. doi: 10.1016/S0734-743X(02)00108-2 [26] 王琳, 王富耻, 王鲁, 等. 空心弹体垂直侵彻混凝土靶板的应变测试研究 [J]. 北京理工大学学报, 2002, 22(4): 453–456. doi: 10.3969/j.issn.1001-0645.2002.04.014WANG L, WANG F C, WANG L, et al. Strain measurement in hollow projectiles impacting concrete targets [J]. Journal of Beijing Institute of Technology, 2002, 22(4): 453–456. doi: 10.3969/j.issn.1001-0645.2002.04.014 [27] 张磊, 任新见, 孔德锋. 钢筋混凝土HJC模型的研究和改进[C]//第四届全国工程安全与防护学术会议. 洛阳, 2014: 134−138.ZHANG L, REN X J, KONG D F. Research and improvement of HJC model of steel reinforced concrete [C]//Proceedings of the 4th National Conference of Engineering Safety and Protection. Luoyang, 2014: 134−138. -