近年来，随着各种破坏性武器的不断发展以及恐怖主义袭击事件不断增多^[1-2]，提高现有工程结构的综合防护能力已经成为防护工作的当务之急。迄今为止，大部分工程结构中使用最广泛的材料仍然是普通混凝土。然而，普通混凝土在高速冲击荷载作用下会发生非常严重的脆性破坏，在结构的背面出现大面积的震塌、崩落以及碎片飞溅现象^[3-6]，对附近人员和设备的安全构成严重威胁。因此，找到合适的工程防护材料成为解决上述问题的关键。活性粉末混凝土（reactive powder concrete，RPC）^[7]是根据最紧密堆积原理制备成的具有高强度和良好韧性的水泥基材料。由于其本身的高强度以及钢纤维和基体之间的桥连作用，使得RPC相较于传统混凝土能够有效地减小迎弹面的开坑尺寸、侵彻深度，降低结构发生整体性破坏的机率^[8-9]。截至目前，已有很多学者针对RPC的抗侵彻性能开展了实验研究^[10-11]，但是关于RPC的实验数据仍然有限，无法系统地指导现有工程结构设计。进行大型抗侵彻实验存在实验成本高昂、周期长、动态数据采集困难等问题。

近些年，随着有限元技术的不断发展，通过数值模拟方法预测混凝土类材料在高速冲击荷载作用下的破坏形态日益受到学者们的青睐。充分发挥有限元软件的优势，可以有效地指导实验、工程顺利进行。然而，目前对RPC材料进行动态数值模拟所选用的本构模型仍然是传统的混凝土本构模型，例如：HJC模型^[12]、RHT模型^[13]、CSCM模型^[14]、K&C模型^[15]，这些模型都无法准确地反映RPC的拉伸延性特征以及应变率效应。因此，要想准确地描述RPC材料的基本力学特性，需要对上述模型的参数作出系统地调整。在以往的RPC高速冲击数值模拟研究中，运用较多的是K&C本构模型，该模型综合考虑了混凝土类材料的应变率效应、应变硬化软化现象、剪胀效应、围压效应，能够准确地描述混凝土类材料在大变形、高应变率和高静水压力下的力学特性。

本文中，首先对RPC靶体进行多次侵彻实验，得到RPC靶体破坏数据，并计算出RPC的Forrestal^[16-17]公式中与靶体材料相关的系数；然后基于K&C模型和现有RPC的基本力学性能实验数据，对K&C本构模型的强度面参数、损伤参数、损伤演化模型、应变率效应、状态方程参数进行修正，采用修正的K&C模型参数模拟上述多次侵彻实验并与实验结果进行对比，验证模拟结果的正确性；最后在质量为10 kg、直径为80 mm弹体正侵彻情况下，对长2 200 mm、宽2 200 mm、高1 260 mm的RPC靶体进行数值预测，得到弹体速度与侵彻深度的关系，并利用Forrestal公式验证数值预测结果的合理性，得到实验过程中弹体的最大加速度和弹体贯穿靶体时的极限速度，以期为实验的顺利进行提供指导。

1.   实验设计

1.1   靶体原材料及制备方法

实验所用原料为硅酸盐水泥（PC），二氧化硅质量分数大于95%、比表面积为18 200 m²/kg的硅灰（SF），减水率大于40%的聚羧酸高效减水剂，比表面积大于10 000 m²/kg的矿渣，细度模数为2.5的普通砂，长度为13 mm、直径为0.175 mm的钢纤维，普通自来水。将以上原料按表1的配合比配制钢纤维掺量为2%的RPC。RPC靶体的浇筑采用强制式卧轴搅拌机搅拌，先将混合好的水泥、粉煤灰、硅灰、精细砂与减水剂干拌2 min，随后加入水搅拌2 min，然后把钢纤维均匀加入继续搅拌5 min，并浇筑到预先制作好的钢桶当中。同时也浇筑了基本力学性能测试试件，常温养护28 d后的测试结果如表2所示，其中f_c和f_t分别为材料的无侧限抗压强度和无侧限抗拉强度，E为弹性模量，μ为泊松比，ρ为密度。

表  1  RPC材料配合比

Table  1.  Mixture proportions of RPC

kg/m3
材料	胶凝材料	砂	减水剂	水	钢纤维
RPC	1 238	928	17.7	234.3	159

下载: 导出CSV 

| 显示表格

表  2  RPC基本力学性能参数

Table  2.  Basic mechanical performance parameters of RPC

材料	fc/MPa	ft/MPa	E/GPa	μ	ρ/（g·cm−3）
RPC	120	9.27	46.2	0.22	2.44

下载: 导出CSV 

| 显示表格

1.2   弹体与靶体

实验所用弹体材质为30GrMnSiNi2A低合金超高强度钢，淬火后抗拉屈服强度可以达到1 650 MPa，由25 mm口径滑膛炮发射，在580 m/s弹速下弹体为刚性侵彻，弹体外形与尺寸如图1所示。弹体直径为25 mm，头部曲径比r_CRH=3，长径比为6，壁厚与弹径比为0.14，底托采用尼龙材质。弹体内部填入惰性材料以调整质量，平均发射质量为352.9 g。RPC靶体为圆柱体，尺寸为$\varnothing$600 mm×600 mm，侧面用4 mm厚钢板箍紧，直径为20倍弹径，可以忽略边界效应。

图  1  弹体尺寸

Figure  1.  Projectile sizes

下载: 全尺寸图片 幻灯片

1.3   侵彻实验方法

实验布置如图2所示。发射装置为25 mm口径滑膛炮，靶体安放在支架上，调整木楔让靶体表面与炮口轴线垂直，炮口距离靶体3 m。由滑膛炮分别发射3发弹丸对RPC靶体进行多次侵彻实验。着靶速度由测速板与电子计时仪测出。高速摄影机记录弹体着靶姿态，拍摄弹体飞行和弹靶初始作用过程。

图  2  侵彻实验设备布置

Figure  2.  Arrangement of penetration experiment equipments

下载: 全尺寸图片 幻灯片

2.   侵彻实验与结果分析

2.1   实验结果

对RPC靶体分别进行3次打击，打击位置均为靶体正中心，弹体着靶时的速度分别为511.5、552.5、560.0 m/s，图3给出了高速摄影机记录的511.5 m/s速度下典型的侵彻过程。

图  3  高速摄影机记录的弹体撞击靶体过程

Figure  3.  Process of the projectile impacting the targetrecorded by a high-speed camera

下载: 全尺寸图片 幻灯片

侵彻的实验结果如表3所示，其中v₀为弹体的初速度，h为弹体侵入靶体沿轴线方向的深度，S和H分别为迎弹面弹坑的面积和深度，N为迎弹面裂纹条数，W_max为最大裂纹宽度。

表  3  靶体多次侵彻实验结果

Table  3.  Experimental results of targets subjected to multiple penetrations

侵彻次数	v0/（m·s−1）	h/mm	S/cm2	H/mm	N	Wmax/mm
1	511.5	129.1	329.7	59.5	0	0
2	552.5	257.4	344.1	79.8	12	2
3	560.0	290.3	354.8	114.9	13	13

下载: 导出CSV 

| 显示表格

从表3可以看出：随着侵彻次数的增加，弹体的侵彻深度、迎弹面的弹坑深度、裂纹条数、最大裂纹宽度增加明显，但迎弹面弹坑面积基本没有变化。此外，为了准确预测混凝土类材料受弹体侵彻的侵彻深度，目前已有几十种经验公式可供参考，其中Forrestal公式较符合本次实验条件。选用第1次侵彻的实验结果来确定Forrestal公式中与靶体材料相关的系数，本实验侵彻深度大于2倍弹径，符合Forrestal公式的适用条件，侵彻深度L的计算公式为：

式中：M为弹体质量（kg），d为弹体直径（m），$S^{\!*}$为实验确定的常数，参数B和${v_{1}^2}$由下式给出：

联立式（1）～（2）可得：

为确定RPC Forrestal侵彻深度计算公式中的抗侵彻系数，将弹体第1次侵彻靶体的实验结果与Wu等^[11]、崔亚男^[18]、Feng等^[19]的实验结果代入式（3），分别得到不同强度的RPC Forrestal公式中的抗侵彻系数，如图4所示。可以看出，根据本次实验计算得到的Forrestal公式中的抗侵彻系数与Feng等、Wu等、崔亚男的实验结果较接近，可以认为本次实验结果有效，因此，本文中将120 MPa RPC的Forrestal公式中的抗侵彻系数定为S^*=7.454 3。

图  4  2%钢纤维掺量下不同强度RPC的Forrestal公式中的S^*

Figure  4.  S* in the Forrestal formula of RPC with differentcompressive strengths at 2% steel fiber content

下载: 全尺寸图片 幻灯片

2.2   靶体破坏形态

从图5(a)中可以看到，弹体第1次侵彻RPC靶体后，靶体迎弹面未出现任何裂缝，只是形成了一个较小的弹坑，并且实验过程中发现，第1发弹体侵入靶体一定深度后弹出，这主要是由于钢纤维与RPC基体之间的桥联作用以及RPC本身强度较高，加上弹速较低，导致弹体未能直接侵入到靶体当中。弹体第2次侵入靶体后，迎弹面出现了许多微小的裂缝，但是弹坑面积未出现明显的增加，弹体直接侵入到靶体中，未出现弹出现象，如图5(b)所示。弹体第3次侵入后，迎弹面可见的裂缝数量相较于第2次侵彻没有明显增多，但是裂缝的宽度都有一定程度的增加，出现了一条贯穿的径向裂缝，弹坑面积基本没有发生变化，如图5(c)所示。

图  5  RPC多次侵彻实验结果

Figure  5.  Experimental results of RPC multiple penetrations

下载: 全尺寸图片 幻灯片

为了观察弹体在靶体中的运动轨迹并回收弹体，对实验后的靶体进行了切割解剖，剖面如图6所示。由于第1发弹体在侵入靶体一定深度后直接弹出靶体，图6仅标出了第2发弹体和第3发弹体的运动轨迹，图6中的红色虚线为第2发和第3发弹体侵彻的弹道。可见，由于靶体中存在初始缺陷以及受到第1次侵彻的影响，第2次和第3次侵彻的弹道都发生了明显的偏转。

图  6  多次侵彻实验后的靶体剖面

Figure  6.  Profile of targets after multiple penetration experiments

下载: 全尺寸图片 幻灯片

3.   数值模拟

3.1   有限元模型

根据实验情况建立如图7所示的有限元模型，弹体、靶体、钢箍均采用SOLID164实体单元，并选用全尺寸模型进行计算。弹体与靶体之间设为面面侵蚀接触，靶体与钢箍设为面面自动接触。为了避免计算过程中单元发生畸变，引入MAT_ADD EROSION侵蚀准则，选用最大主应变控制单元失效^[20]，对于本次模拟侵蚀应变的取值，选取第1次侵彻模拟结果和实验数据吻合最佳时的侵蚀应变作为该组实验的侵蚀应变取值，即当单元的最大主应变超过0.15时，单元立即被删除。

图  7  有限元模型

Figure  7.  Finite element model

下载: 全尺寸图片 幻灯片

3.2   材料模型

材料模型选择和模型参数确定是决定模拟结果可靠的关键因素。弹体采用刚体材料模型，钢箍和钢筋均选用LS-DYNA中的P-K随动模型，弹体、钢箍以及钢筋的模型参数如表4所示，其中σ_y为屈服强度。

表  4  弹体、钢箍以及钢筋材料模型参数

Table  4.  Model parameters of projectile, steel culvert and steel bar material

材料	ρ/（kg·m−3）	E/GPa	μ	σy/MPa
弹体	7 850	210	0.3	1 650
钢箍/钢筋	7 800	210	0.3	300

下载: 导出CSV 

| 显示表格

活性粉末混凝土靶体的材料模型选用K&C模型，但K&C模型的参数是基于45.4 MPa混凝土实验数据由程序自动生成的，无法准确描述RPC的拉伸韧性、应变率关系和压力相关性，因此本文中将基于现有的RPC实验数据，通过调整K&C模型的参数，使之可以更好地描述RPC的基本力学特性。

3.2.1   强度面参数

K&C^[15]模型引入了3个失效面：初始屈服面、极限强度面、残余强度面，分别描述混凝土的初始强度、极限强度及残余强度的变化规律，其表达式为：

式中：a_i、a_iy和a_if（i=0，1，2）为由材料三轴围压实验数据确定的常数，p为静水压力，η(λ)为损伤演化关系，$\varphi$为拉压子午比。$\varphi$具体定义为：

当前失效面定义为：

式中：$r_{\rm{f}}$为应变率效应系数（$\psi$），${r'}$为当前子午与压缩子午的比值。

在K&C模型中，强度面的9个参数均参照45.4 MPa混凝土的三轴围压实验数据由LS-DYNA程序自动生成，它们与材料抗压强度的关系表示为：

搜集了近些年与RPC和混凝土相关的三轴围压实验数据^[21-31]，如图8所示，发现在低归一化静水压力下，混凝土与RPC归一化等效应力值基本一致。由于RPC高归一化静水压力的实验数据很难获得，因此本文中假定RPC在高归一化静水压力下归一化等效应力值与混凝土一致，以此获得RPC的极限强度面参数。从图8中还可以看出，采用K&C自动生成的强度面参数明显低于实验数据点，而修正后的RPC强度面参数与实验数据非常吻合。而RPC弹性强度面参数的确定依然按照K&C模型中的方法，即假定Δσ_y=0.45Δσ_m^[15]。Kong等^[32]在修改的K&C模型中认为混凝土类材料的残余强度面平行于极限强度面，并应用于混凝土靶体的贯穿数值模拟当中，取得了不错的模拟效果。本文中也采取相同的方法，这样就可以全面地确定RPC的强度面的参数，具体表示为：

图  8  混凝土和RPC三轴围压实验数据

Figure  8.  Triaxial confining pressure experimental data of concrete and RPC

下载: 全尺寸图片 幻灯片

3.2.2   损伤演化模型

K&C模型的当前失效面将随着损伤因子λ的增大先强化后软化，λ为等效塑性应变的函数，定义为：

式中：${\rm{d}}{\overline \varepsilon _{\rm{p}}} = \sqrt {\dfrac{2}{3}{\rm{d}}\varepsilon _{ij}^{\rm{p}}{\rm{d}}\varepsilon _{ij}^p}$为有效塑性应变增量，${\rm{d}}\varepsilon _{ij}^{\rm{p}}$为有效塑性应变增量张量，b₁和b₂分别为压缩和拉伸损伤参数。

当混凝土类材料三向受拉时，只产生体积应变而不产生形状应变，此时偏应变为零，因此损伤变量也为零。为考虑应力路径接近三向受拉时体积损伤应变的影响，K&C模型中引入了体积损伤变量Δλ：

式中：b₃为控制混凝土三向受拉时软化段的损伤参数；k_d为内变量因子；ε_V为体积应变，ε_{V, y}为屈服点附近的体积应变；f_d则表征应力路径与三向受拉状态靠近的程度。f_d由下式给出：

在K&C模型中，损伤演化模型由13对(λ, η)数据点来定义，这些相邻的数据点由插值函数来连接。其定义是针对单一强度（45.4 MPa）混凝土实验得到的数据，与之吻合较好，但无法应用到其他强度的RPC中，更不能描述RPC材料具有的拉伸延性特征，因此需要提出一个新的、可靠的损伤演化关系来模拟RPC损伤关系。通过不断地对单个单元进行拉伸、压缩模拟，最终确定损伤参数b₁和b₂分别为2.2和–7.25，将反映η和λ关系的13对点定义如表5所示。

表  5  损伤演化函数η(λ)

Table  5.  Damage evolution function η(λ)

λ	η		λ	η
0	0		4.0×10−6	0.51
2.7×10−5	0.62		6.7×10−4	0.37
6.8×10−5	0.92		1.2×10−3	0.27
8.0×10−5	0.99		2.0×10−3	0.20
1.0×10−4	1.00		5.5×10−3	0.10
1.4×10−4	0.96		1.6×10−2	0
2.6×10−4	0.66

下载: 导出CSV 

| 显示表格

图9展示了采用自动生成参数的K&C模型、新参数的K&C模型以及RPC单轴拉伸压缩实验数据，结果表明，新的K&C参数模拟的结果与实验数据更吻合。

图  9  拉伸和压缩应力应变曲线

Figure  9.  Tensile and compressive stress-strain curves

下载: 全尺寸图片 幻灯片

图10给出了体积损伤参数b₃对三轴拉伸应力应变曲线的影响，可以看出，三轴拉伸应力应变曲线的软化梯度随b₃的减小而减小。对三轴拉伸体积损伤参数b₃的取值，Weerhijm等^[33]假定：单轴拉伸的断裂应变和三轴拉伸的断裂应变相等。因此本文中b₃取值为0.018。

图  10  不同b₃值对应的三轴拉伸应力应变曲线

Figure  10.  Triaxial tension strain stress curvescorresponding to different b₃ values

下载: 全尺寸图片 幻灯片

3.2.3   应变率参数

材料的动态力学特性会随着加载速率的变化而变化。学者们针对普通混凝土的这种特性已经开展了深入的研究，发现在高加载速率下，混凝土的压缩强度会增加100%，拉伸强度会增加600%^[34]，并且提出了可靠的应变率公式加以描述，如K&C模型中的应变率公式被广泛应用于混凝土受到冲击、爆炸荷载作用下的数值模拟中。RPC材料的拉伸实验表明^[35]，RPC材料拥有较普通混凝土更明显的率依赖性。因此，在RPC材料动态数值模拟中考虑应变率效应就显得非常重要。此前，许多学者^[36-38]对RPC受冲击和爆炸荷载的数值模拟都采用基于普通混凝土提出的应变率公式，这显然无法较准确地模拟RPC在高应变率荷载下的动态力学响应。Park等^[38]对基体强度分别为56、81和180 MPa的RPC进行了拉伸应变率敏感性试验研究，并提出了拉伸应变率公式（式(14)）。Lin^[39]基于现有RPC的动态压缩实验数据，提出了动态压缩应变率公式（式(15)），并将其应用于RPC板抗爆数值模拟中，取得了不错的模拟效果。本文中也采用相同的计算公式（式(14)~(15)），并据此计算出RPC的应变率效应数据点，如表6所示。此外，最新版本的K&C模型中将混凝土动态拉伸增强因子和动态压缩增强因子上限分别设置为9.7和2.94，以避免高估混凝土材料在高应变率下的应变率效应。本文中也将这一设置引入到RPC的动态数值模拟：

表  6  活性粉末混凝土的K&C模型应变率效应特征点取值

Table  6.  K&C model strain rate characteristic points of reactive powder concrete

$\dot \varepsilon/{\rm{s}^{-1} }$	$\psi$		$\dot \varepsilon/{\rm{s}^{-1} }$	$\psi$		$\dot \varepsilon/{\rm{s}^{-1} }$	$\psi$
–30 000	9.97		–10	1.27		–1×10–4	1.07
–4 782	9.97		–3	1.24		–1×10–5	1.03
–1 000	5.41		–1	1.22		0	1.00
–300	3.45		–0.1	1.18		30	1.00
–100	2.29		–0.01	1.14		265	2.94
–25	1.28		–1×10–3	1.11		30 000	2.94

下载: 导出CSV 

| 显示表格

式中：${\psi _{\rm{t}}}$、${\psi _{\rm{c}}}$分别为动态拉伸增强因子和动态压缩增强因子，${\dot \varepsilon _0} = {10^{ - 6}}\;{{\rm{s}}^{ - 1}}$为静态应变率。

3.2.4   状态方程参数

K&C模型采用8号状态方程来定义混凝土压力与体积应变之间的关系：

式中：${E}_{0}$为初始体积内能，γ₀为温度特征系数，C(ε_V)和$T\left({\varepsilon }_{V}\right)$分别为压力和温度与体积应变之间的关系。在加载阶段，压力由式(4)确定，卸载沿卸载体积模量至压力截断点（图11(a)），重新加载首先沿卸载路径至卸载开始点，然后沿加载路径继续上升。

图  11  状态方程

Figure  11.  State equation

下载: 全尺寸图片 幻灯片

尽管K&C模型中的8号状态方程能够很好地描述混凝土压力与体积应变之间的关系，但是由于RPC的孔隙率远小于普通混凝土，导致原来的8号状态方程无法准确地描述RPC的压力与体积应变的关系。材料的状态方程所需的实验数据通常需要进行飞片撞击试验。因此本文中首先通过基本力学性能实验测试求得RPC的弹性模量和泊松比，并按照式（17）求得RPC材料的体积模量K_V，然后基于Marsh^[40]、高乐^[41]和严少华等^[42]的Hugoniot试验数据校核K&C模型中的状态方程，以获得适合于RPC的8号状态方程参数：

图11(b)展示了校核结果，可以看出，自动生成的8号状态方程明显低于试验数据，而校核后的8号状态方程与试验数据基本一致。校核后的状态方程参数如表7所示，其中ε_V为体积应变，σ_V为相应体积应变所对应的体积应力，K_av代表卸载模量。其他参数对计算结果基本没有影响，因此选用K&C模型自动生成的参数，至此关于RPC的K&C模型参数完全确定。

表  7  RPC的K&C模型8号状态方程参数

Table  7.  Parameters of No. 8 equation of state in the K&C model of RPC

εV
εV1	εV2	εV3	εV4	εV5	εV6	εV7	εV8	εV9	εV10
0	0.0015	0.0043	0.0101	0.0305	0.0513	0.0726	0.0943	0.174	0.208
σV/GPa
σV1	σV2	σV3	σV4	σV5	σV6	σV7	σV8	σV9	σV10
0	0.041	0.094	0.292	0.881	1.622	2.511	3.573	8.714	11.579
Kav/GPa
Kav1	Kav2	Kav3	Kav4	Kav5	Kav6	Kav7	Kav8	Kav9	Kav10
27.5	27.5	27.885	29.288	34.843	40.425	45.980	50.186	112.915	137.5

下载: 导出CSV 

| 显示表格

3.3   数值模拟结果与分析

为进一步观察靶体内部损伤状况并验证所选K&C模型参数的准确性，利用LS-DYNA中的重启动方法分别对3次侵彻进行数值模拟，该方法可以记录第1次弹体冲击RPC靶体后的应力、应变和损伤状态并输出重启动文件。在分析弹体第2次冲击RPC靶体时，LS-DYNA会读取弹体第1次冲击RPC靶体后的应力、应变和损伤状态等数据并继续进行计算，这样便可以实现对弹体多次冲击受损后RPC靶体的数值模拟，模拟结果如图12～14所示。利用数值模拟得到弹体第1次的侵彻深度和开坑面积分别为13.2 cm、315.2 cm²，而实验结果分别为12.9 cm、329.7 cm²。第2次侵彻深度和开坑面积的模拟结果分别为25.9 cm、336.4 cm²，而实验结果分别为25.7 cm、344.1 cm²。可知前两次模拟得到的侵彻深度与实验结果基本一致，误差不超过5%。图14给出了第3次侵彻的模拟结果，第3次模拟得到的侵彻深度为34.0 cm，与实验得到的侵彻深度（29.0 cm）存在较大误差。并且模拟得到的弹体偏转角度小于实验结果，这样就使得模拟得到的侵彻深度大于实验测量的结果。这主要是由于在进行有限元分析时将RPC靶体视为完全均质材料，因此在第1次和第2次模拟结果中可以看到，RPC靶体内部的损伤也是均匀分布的，弹体侵彻RPC靶体的弹道没有明显偏转。然而实验中RPC靶体并非完全均质，一部分区域存在初始缺陷，造成弹体第1次冲击后靶体内部的损伤更加不均匀。当弹体第2次冲击靶体时，弹体四周受到的阻力不一致，导致弹体沿损伤更严重的区域偏转。因此才会看到实验中第2发弹体和第3发弹体都发生了明显的偏转，而模拟结果无法再现实验中弹体发生的明显偏转。

图  12  第1次侵彻RPC靶体的模拟结果

Figure  12.  Simulation results of the first penetration of the RPC target

下载: 全尺寸图片 幻灯片

图  13  第2次侵彻RPC靶体的模拟结果

Figure  13.  Simulation results of the second penetration of the RPC target

下载: 全尺寸图片 幻灯片

图  14  第3次侵彻RPC靶体的模拟结果

Figure  14.  Simulation results of the third penetration of the RPC target

下载: 全尺寸图片 幻灯片

此外，K&C模型中将损伤定义为D=2λ/(λ+λ_m)，当D从0到1时，表示材料进入了强化阶段，当D从1到2时，表示材料进入了软化阶段，而当D=2时认为单元完全失效。从损伤云图12～14中可以看到，随着侵彻次数增多，靶体内部的损伤也不断扩大，同时迎弹面的裂缝数量也在增多。在迎弹面，D＞1.8的损伤面积没有发生太大变化，这与多次侵彻RPC靶体后迎弹面开坑面积基本相同的实验结果一致。因此可以认为RPC的K&C模型参数和所用的模拟方法具有可靠性。

4.   数值预测

基于上述RPC的K&C模型参数和模拟方法，采用LS-DYNA软件对RPC侵彻实验进行数值预测，为实验测试系统的选用和弹速设计提供指导。该实验的靶体尺寸为2 200 mm×2 200 mm ×1 260 mm，弹体直径为80 mm，长为361.35 mm，弹体质量为10 kg并为靶体配置钢筋，具体参数如图15所示。

图  15  弹体与钢筋分布示意图（单位：mm）

Figure  15.  Schematic diagram of the distribution of projectiles and reinforcement (unit: mm)

下载: 全尺寸图片 幻灯片

根据实验情况的对称性，为了提高计算效率和节省计算资源，建立如图16所示的1/4有限元模型，弹体与RPC靶体均采用SOLID164实体单元，钢筋采用BEAM161梁单元，靶体与钢筋之间选用CONSTRAINED_LAGRANGE_IN_SOLID耦合方式进行相互作用，弹体与靶体之间采用CONTACT_ERODING_SURFACE_TO_SURFACE面面侵蚀接触关键字控制，模型单位制m-s-kg，共包含节点2 395 628个、六面体单元1 592 608个、梁单元7 267个。弹体和钢筋均选用P-K随动强化模型，靶体采用修正的K&C模型参数。

图  16  有限元模型

Figure  16.  Finite element model

下载: 全尺寸图片 幻灯片

4.1   数值预测结果与分析

图17(a)～(c)为卵形弹以850 m/s的速度正侵彻RPC靶的数值模拟结果。可以看到，弹体高速侵彻对RPC产生了明显的局部破坏，形成了直径为0.356 m、深度为0.192 m的弹坑，同时形成了直径略大于弹径的圆柱形孔道，迎弹面产生了数十条主裂缝，弹体侵彻靶体所造成的损伤主要集中在柱形孔道。图17(d)～(e)分别显示了弹体的加速度、速度、位移时程曲线，可以看到弹体的最大加速度为63 000g，弹体的速度在侵入靶体后一直减小，直至零。

图  17  850 m/s弹速下的侵彻结果

Figure  17.  Penetration results at the projectile velocity of 850 m/s

下载: 全尺寸图片 幻灯片

图18显示了弹体以1 150 m/s的速度贯穿RPC靶体的模拟结果，同850 m/s的侵彻结果一样，靶体形成了直径为0.340 m、深度为0.213 m的漏斗坑。不同的是迎弹面产生了更多裂缝，且靶体内部损伤更严重，此外背弹面出现了弹体贯穿形成了直径为0.505 m、深度为0.264 m的贯穿坑和宏观主裂纹。从图18(d)～(e)中可以看出，弹体加速度在0.2 ms附近达到峰值，最大值为73 000g，相较于850 m/s的侵彻结果，加速度的峰值提高了16%，同时弹体在贯穿靶体后的剩余速度为100 m/s。

图  18  1 150 m/s弹速下的贯穿结果

Figure  18.  Penetration results at the projectile velocity of 1 150 m/s

下载: 全尺寸图片 幻灯片

通过上述数值模拟方法，可以得到不同速度弹体侵入RPC靶体后的侵彻深度，不同弹速下弹体侵入RPC靶体的深度与RPC-Forrestal公式的预测结果如图19所示，可以发现，采用本文中的数值模拟方法预测的侵彻深度与RPC-Forrestal公式预测的结果较吻合。因此可以认为本次数值模拟的结果有效。另外，根据数值模拟和RPC-Forrestal公式可以发现：当弹体速度超过1 100 m/s时，弹体会贯穿RPC靶体。

图  19  弹速与侵彻深度之间的关系

Figure  19.  Relationship between projectile velocity and penetration depth

下载: 全尺寸图片 幻灯片

5.   结　论

开展了直径25 mm弹体多次侵彻活性粉末混凝土靶实验，得到了靶体的破坏数据；确定了活性粉末混凝土的K&C模型参数；通过多次侵彻实验验证了改进后K&C模型参数的正确性，并对直径为80 mm弹体冲击2 200 mm×2 200 mm×1 260 mm活性粉末混凝土靶体的实验结果进行数值预测，得到以下结论：

（1）在相同弹速下沿同一点逐次侵彻RPC靶体，随着侵彻次数的增加，弹体的侵彻深度、迎弹面的弹坑深度、裂纹条数、最大裂纹宽度明显增大，但是迎弹面弹坑面积基本没有变化；

（2）通过侵彻实验确定了抗压强度为120 MPa RPC的Forrestal侵彻深度计算公式中抗侵彻系数S^*=7.454 3；

（3）基于RPC的准静态单轴压缩、拉伸实验、三轴围压实验、拉伸压缩应变率实验和飞片撞击实验数据，提出了一组适用于RPC的K&C本构模型参数；

（4）采用非线性有限元方法对侵彻实验进行了数值模拟，模拟结果与实验结果吻合较好，并采用同样的方法对实验进行了数值预测，得到了两种预设弹速下弹体的加速度、速度、位移时程曲线以及不同弹速下弹体侵入靶体的深度，并利用Forrestal经验公式验证了数值预测结果的准确性，最终确定了弹体贯穿靶体的速度至少要达到1 100 m/s。