Experimental investigation and numerical prediction on resistance of reactive powder concrete to multiple penetration
-
摘要: 活性粉末混凝土（reactive powder concrete，RPC）具有超高的强度和优异的阻裂性能。为了研究RPC在多次冲击荷载下的损伤规律，采用25 mm口径滑膛炮对直径为600 mm、高600 mm的RPC圆柱形靶体进行了多次侵彻实验，得到了每次侵彻后靶体的破坏数据，并根据实验数据确定了Forrestal经验公式中的相关系数。基于K&C本构模型和现有RPC基本力学性能的实验数据，修正了K&C模型的强度面参数、损伤参数、状态方程参数、损伤演化模型以及应变率效应相关参数，系统地确定了RPC的K&C模型参数。采用LS-DYNA软件中的重启动功能模拟了弹体多次侵彻RPC靶体的破坏结果，模拟结果与实验结果基本一致，验证了模拟方法的有效性。对长2 200 mm、宽2 200 mm、高1 260 mm的RPC靶体抗侵彻实验进行了数值预测，得到了侵彻深度与弹速之间的关系、弹体贯穿靶体时的极限速度以及弹体侵彻过程中的峰值加速度。Abstract: Reactive powder concrete (RPC) has ultra-high strength and excellent crack resistance. To study the damage law of the RPC subjected to multiple impact loads, a 25 mm caliber smoothbore gun was used to penetrate the RPC cylindrical target with the diameter of 600 mm and the height of 600 mm. In addition, the experimental data of the target after each penetration was obtained, and the correlation coefficient in the Forrestal empirical formula was determined. Based on the K&C constitutive model and the existing experimental data of the RPC, the model parameters for the RPC were determined systematically by modifying the strength and surface parameters, damage parameters, equation-of-state parameters, damage evolution model, the strain rate effect. The restart function in the LS-DYNA software was used to simulate the damage results of the projectile repeatedly penetrating the RPC target. The simulation results are basically consistent with the experimental results, and the effectiveness of the simulation method is verified. Finally, the numerical prediction of the penetration resistance experiment of the RPC target with the length of 2 200 mm, the width of 2 200 mm, and the height of 1 260 mm was carried out. The relationship between the penetration depth and the projectile velocity, the minimum velocity of the projectile passing through the target and the peak acceleration during projectile penetration were obtained.
-
近年来，随着各种破坏性武器的不断发展以及恐怖主义袭击事件不断增多[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 m2/kg的硅灰（SF），减水率大于40%的聚羧酸高效减水剂，比表面积大于10 000 m2/kg的矿渣，细度模数为2.5的普通砂，长度为13 mm、直径为0.175 mm的钢纤维，普通自来水。将以上原料按表1的配合比配制钢纤维掺量为2%的RPC。RPC靶体的浇筑采用强制式卧轴搅拌机搅拌，先将混合好的水泥、粉煤灰、硅灰、精细砂与减水剂干拌2 min，随后加入水搅拌2 min，然后把钢纤维均匀加入继续搅拌5 min，并浇筑到预先制作好的钢桶当中。同时也浇筑了基本力学性能测试试件，常温养护28 d后的测试结果如表2所示，其中fc和ft分别为材料的无侧限抗压强度和无侧限抗拉强度，E为弹性模量，μ为泊松比，ρ为密度。
表 1 RPC材料配合比Table 1. Mixture proportions of RPCkg/m3 材料 胶凝材料 砂 减水剂 水 钢纤维 RPC 1 238 928 17.7 234.3 159 表 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 1.2 弹体与靶体
实验所用弹体材质为30GrMnSiNi2A低合金超高强度钢，淬火后抗拉屈服强度可以达到1 650 MPa，由25 mm口径滑膛炮发射，在580 m/s弹速下弹体为刚性侵彻，弹体外形与尺寸如图1所示。弹体直径为25 mm，头部曲径比rCRH=3，长径比为6，壁厚与弹径比为0.14，底托采用尼龙材质。弹体内部填入惰性材料以调整质量，平均发射质量为352.9 g。RPC靶体为圆柱体，尺寸为
∅ 600 mm×600 mm，侧面用4 mm厚钢板箍紧，直径为20倍弹径，可以忽略边界效应。1.3 侵彻实验方法
实验布置如图2所示。发射装置为25 mm口径滑膛炮，靶体安放在支架上，调整木楔让靶体表面与炮口轴线垂直，炮口距离靶体3 m。由滑膛炮分别发射3发弹丸对RPC靶体进行多次侵彻实验。着靶速度由测速板与电子计时仪测出。高速摄影机记录弹体着靶姿态，拍摄弹体飞行和弹靶初始作用过程。
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所示，其中v0为弹体的初速度，h为弹体侵入靶体沿轴线方向的深度，S和H分别为迎弹面弹坑的面积和深度，N为迎弹面裂纹条数，Wmax为最大裂纹宽度。
表 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 从表3可以看出：随着侵彻次数的增加，弹体的侵彻深度、迎弹面的弹坑深度、裂纹条数、最大裂纹宽度增加明显，但迎弹面弹坑面积基本没有变化。此外，为了准确预测混凝土类材料受弹体侵彻的侵彻深度，目前已有几十种经验公式可供参考，其中Forrestal公式较符合本次实验条件。选用第1次侵彻的实验结果来确定Forrestal公式中与靶体材料相关的系数，本实验侵彻深度大于2倍弹径，符合Forrestal公式的适用条件，侵彻深度L的计算公式为：
L=2Mπd2ρBln(1+ρBv21S∗fc)+2d (1) 式中：M为弹体质量（kg），d为弹体直径（m），
S∗ 为实验确定的常数，参数B和v21 由下式给出：B=8rCRH−124r2CRH,v21=2Mv20−πd3S∗fc2M+πd3Bρ (2) 联立式（1）～（2）可得：
S∗=2ρBMv20fc(2M+πd3Bρ)exp[(L−2d)πd2ρB/(2M)+ρBπd3/(2M+πd3Bρ)−1] (3) 为确定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 content2.2 靶体破坏形态
从图5(a)中可以看到，弹体第1次侵彻RPC靶体后，靶体迎弹面未出现任何裂缝，只是形成了一个较小的弹坑，并且实验过程中发现，第1发弹体侵入靶体一定深度后弹出，这主要是由于钢纤维与RPC基体之间的桥联作用以及RPC本身强度较高，加上弹速较低，导致弹体未能直接侵入到靶体当中。弹体第2次侵入靶体后，迎弹面出现了许多微小的裂缝，但是弹坑面积未出现明显的增加，弹体直接侵入到靶体中，未出现弹出现象，如图5(b)所示。弹体第3次侵入后，迎弹面可见的裂缝数量相较于第2次侵彻没有明显增多，但是裂缝的宽度都有一定程度的增加，出现了一条贯穿的径向裂缝，弹坑面积基本没有发生变化，如图5(c)所示。
为了观察弹体在靶体中的运动轨迹并回收弹体，对实验后的靶体进行了切割解剖，剖面如图6所示。由于第1发弹体在侵入靶体一定深度后直接弹出靶体，图6仅标出了第2发弹体和第3发弹体的运动轨迹，图6中的红色虚线为第2发和第3发弹体侵彻的弹道。可见，由于靶体中存在初始缺陷以及受到第1次侵彻的影响，第2次和第3次侵彻的弹道都发生了明显的偏转。
3. 数值模拟
3.1 有限元模型
根据实验情况建立如图7所示的有限元模型，弹体、靶体、钢箍均采用SOLID164实体单元，并选用全尺寸模型进行计算。弹体与靶体之间设为面面侵蚀接触，靶体与钢箍设为面面自动接触。为了避免计算过程中单元发生畸变，引入MAT_ADD EROSION侵蚀准则，选用最大主应变控制单元失效[20]，对于本次模拟侵蚀应变的取值，选取第1次侵彻模拟结果和实验数据吻合最佳时的侵蚀应变作为该组实验的侵蚀应变取值，即当单元的最大主应变超过0.15时，单元立即被删除。
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 活性粉末混凝土靶体的材料模型选用K&C模型，但K&C模型的参数是基于45.4 MPa混凝土实验数据由程序自动生成的，无法准确描述RPC的拉伸韧性、应变率关系和压力相关性，因此本文中将基于现有的RPC实验数据，通过调整K&C模型的参数，使之可以更好地描述RPC的基本力学特性。
3.2.1 强度面参数
K&C[15]模型引入了3个失效面：初始屈服面、极限强度面、残余强度面，分别描述混凝土的初始强度、极限强度及残余强度的变化规律，其表达式为：
Δσy={a0y+p/(a1y+a2yp)p≥0.15fc1.35ft+3p(1−3ft/fc)0≤p＜0.15fc1.35(p+ft)p＜0 (4) Δσm={a0+p/(a1+a2p)p≥fc/3(1.5/φ)(p+ft)0≤p＜fc/33(p/η+ft)p＜0 (5) Δσr=a0f+p/(a1f+a2fp) (6) 式中：ai、aiy和aif（i=0，1，2）为由材料三轴围压实验数据确定的常数，p为静水压力，η(λ)为损伤演化关系，
φ 为拉压子午比。φ 具体定义为：φ(p)={0.5p≤00.5+1.5ft/fcp=fc/31.15fc/[a0+2.3/(3a1+2.3a2fc)]p=2.3fc/30.753p=3fc1p≥8.45fc (7) 当前失效面定义为：
Δσ=rfΔσ(p/rf)=√3J2={[η(Δσm−Δσy)+Δσy]r′Hardening[η(Δσm−Δσr)+Δσr]r′Softing (8) 式中：
rf 为应变率效应系数（ψ ），r′ 为当前子午与压缩子午的比值。在K&C模型中，强度面的9个参数均参照45.4 MPa混凝土的三轴围压实验数据由LS-DYNA程序自动生成，它们与材料抗压强度的关系表示为：
{a0=0.2956fc,a1=0.4463,a2=0.0808/fca0y=0.2232fc,a1y=0.6250,a2y=0.2575/fca0f=0,a1f=0.4417,a2f=0.1183/fc (9) 搜集了近些年与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{a0=0.4426fc,a1=0.5698,a2=0.02516/fca0y=0.2797fc,a1y=0.8989,a2y=0.06850/fca0f=0,a1f=0.5698,a2f=0.02516/fc (10) 3.2.2 损伤演化模型
K&C模型的当前失效面将随着损伤因子λ的增大先强化后软化，λ为等效塑性应变的函数，定义为：
λ={∫t0d¯εp/[rf(1+p/ft)b1]p＞0∫t0d¯εp/[rf(1+p/ft)b2]p≤0 (11) 式中：
d¯εp=√23dεpijdεpij 为有效塑性应变增量，dεpij 为有效塑性应变增量张量，b1和b2分别为压缩和拉伸损伤参数。当混凝土类材料三向受拉时，只产生体积应变而不产生形状应变，此时偏应变为零，因此损伤变量也为零。为考虑应力路径接近三向受拉时体积损伤应变的影响，K&C模型中引入了体积损伤变量Δλ：
Δλ=b3fdkd(εV−εV,y) (12) 式中：b3为控制混凝土三向受拉时软化段的损伤参数；kd为内变量因子；εV为体积应变，εV, y为屈服点附近的体积应变；fd则表征应力路径与三向受拉状态靠近的程度。fd由下式给出：
fd={1−|√3J2/p|/0.10≤|√3J2/p|≤0.10|√3J2/p|＞0.1 (13) 在K&C模型中，损伤演化模型由13对(λ, η)数据点来定义，这些相邻的数据点由插值函数来连接。其定义是针对单一强度（45.4 MPa）混凝土实验得到的数据，与之吻合较好，但无法应用到其他强度的RPC中，更不能描述RPC材料具有的拉伸延性特征，因此需要提出一个新的、可靠的损伤演化关系来模拟RPC损伤关系。通过不断地对单个单元进行拉伸、压缩模拟，最终确定损伤参数b1和b2分别为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 图9展示了采用自动生成参数的K&C模型、新参数的K&C模型以及RPC单轴拉伸压缩实验数据，结果表明，新的K&C参数模拟的结果与实验数据更吻合。
图10给出了体积损伤参数b3对三轴拉伸应力应变曲线的影响，可以看出，三轴拉伸应力应变曲线的软化梯度随b3的减小而减小。对三轴拉伸体积损伤参数b3的取值，Weerhijm等[33]假定：单轴拉伸的断裂应变和三轴拉伸的断裂应变相等。因此本文中b3取值为0.018。
图 10 不同b3值对应的三轴拉伸应力应变曲线Figure 10. Triaxial tension strain stress curvescorresponding to different b3 values3.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˙ε/s−1 ψ ˙ε/s−1 ψ ˙ε/s−1 ψ –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 ψt={(˙ε/˙ε0)0.01465˙ε≤25s−10.002352(˙ε/˙ε0)0.3735˙ε＞25s−1 (14) ψc={1˙ε≤30s−10.728+0.008337(˙ε/˙ε0)˙ε＞30s−1 (15) 式中：
ψt 、ψc 分别为动态拉伸增强因子和动态压缩增强因子，˙ε0=10−6s−1 为静态应变率。3.2.4 状态方程参数
K&C模型采用8号状态方程来定义混凝土压力与体积应变之间的关系：
p=C(εV)+γ0T(εV)E0 (16) 式中：
E0 为初始体积内能，γ0为温度特征系数，C(εV)和T(εV) 分别为压力和温度与体积应变之间的关系。在加载阶段，压力由式(4)确定，卸载沿卸载体积模量至压力截断点（图11(a)），重新加载首先沿卸载路径至卸载开始点，然后沿加载路径继续上升。尽管K&C模型中的8号状态方程能够很好地描述混凝土压力与体积应变之间的关系，但是由于RPC的孔隙率远小于普通混凝土，导致原来的8号状态方程无法准确地描述RPC的压力与体积应变的关系。材料的状态方程所需的实验数据通常需要进行飞片撞击试验。因此本文中首先通过基本力学性能实验测试求得RPC的弹性模量和泊松比，并按照式（17）求得RPC材料的体积模量KV，然后基于Marsh[40]、高乐[41]和严少华等[42]的Hugoniot试验数据校核K&C模型中的状态方程，以获得适合于RPC的8号状态方程参数：
KV=E3(1−2μ) (17) 图11(b)展示了校核结果，可以看出，自动生成的8号状态方程明显低于试验数据，而校核后的8号状态方程与试验数据基本一致。校核后的状态方程参数如表7所示，其中εV为体积应变，σV为相应体积应变所对应的体积应力，Kav代表卸载模量。其他参数对计算结果基本没有影响，因此选用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 3.3 数值模拟结果与分析
为进一步观察靶体内部损伤状况并验证所选K&C模型参数的准确性，利用LS-DYNA中的重启动方法分别对3次侵彻进行数值模拟，该方法可以记录第1次弹体冲击RPC靶体后的应力、应变和损伤状态并输出重启动文件。在分析弹体第2次冲击RPC靶体时，LS-DYNA会读取弹体第1次冲击RPC靶体后的应力、应变和损伤状态等数据并继续进行计算，这样便可以实现对弹体多次冲击受损后RPC靶体的数值模拟，模拟结果如图12～14所示。利用数值模拟得到弹体第1次的侵彻深度和开坑面积分别为13.2 cm、315.2 cm2，而实验结果分别为12.9 cm、329.7 cm2。第2次侵彻深度和开坑面积的模拟结果分别为25.9 cm、336.4 cm2，而实验结果分别为25.7 cm、344.1 cm2。可知前两次模拟得到的侵彻深度与实验结果基本一致，误差不超过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模型参数。
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 depth5. 结 论
开展了直径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。
-
图 3 高速摄影机记录的弹体撞击靶体过程
Figure 3. Process of the projectile impacting the targetrecorded by a high-speed camera
图 4 2%钢纤维掺量下不同强度RPC的Forrestal公式中的S*
Figure 4. S* in the Forrestal formula of RPC with differentcompressive strengths at 2% steel fiber content
图 8 混凝土和RPC三轴围压实验数据
Figure 8. Triaxial confining pressure experimental data of concrete and RPC
图 10 不同b3值对应的三轴拉伸应力应变曲线
Figure 10. Triaxial tension strain stress curvescorresponding to different b3 values
图 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
图 15 弹体与钢筋分布示意图（单位：mm）
Figure 15. Schematic diagram of the distribution of projectiles and reinforcement (unit: mm)
图 17 850 m/s弹速下的侵彻结果
Figure 17. Penetration results at the projectile velocity of 850 m/s
图 18 1 150 m/s弹速下的贯穿结果
Figure 18. Penetration results at the projectile velocity of 1 150 m/s
图 19 弹速与侵彻深度之间的关系
Figure 19. Relationship between projectile velocity and penetration depth
表 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
表 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
表 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
表 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
表 6 活性粉末混凝土的K&C模型应变率效应特征点取值
Table 6. K&C model strain rate characteristic points of reactive powder concrete
˙ε/s−1 ψ ˙ε/s−1 ψ ˙ε/s−1 ψ –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
表 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
-
[1] 任辉启, 穆朝明, 刘瑞朝, 等. 精确制导武器侵彻效应与工程防护 [M]. 北京: 科学出版社, 2016.REN H Q, MU C M, LIU R C, et al. Penetration effects of precision guided weapons and engineering protection [M]. Beijing: Science Press, 2016. [2] 王振宇, 冯进技, 张殿臣. 国外小型钻地核武器的发展及防护建议 [C] // 中国土木工程学会防护工程分会第五届理事会暨第九次学术会议. 长春: 中国土木工程学会, 2004: 66–69.WANG Z Y, FENG J J, ZHANG D C. Development of foreign small earth-penetrating nuclear weapons and relevant preventing measures [C] // The 5th Council and 9th Academic Conference of Protection Engineering Branch of China Society of Civil Engineering, Changchun: China Civil Engineering Society, 2004: 66–69. [3] DANCYGIER A N, YANKELEVSKY D Z, JAEGERMANN C. Response of high performance concrete plates to impact of non-deforming projectiles [J]. International Journal of Impact Engineering, 2007, 34(11): 1768–1779. DOI: 10.1016/j.ijimpeng.2006.09.094. [4] GEBBEKEN N, GREULICH S, PIETZSCH A. Hugoniot properties for concrete determined by full-scale detonation experiments and flyer-plate-impact tests [J]. International Journal of Impact Engineering, 2006, 32(12): 2017–2031. DOI: 10.1016/j.ijimpeng.2005.08.003. [5] LEPPÄNEN J. Experiments and numerical analyses of blast and fragment impacts on concrete [J]. International Journal of Impact Engineering, 2005, 31(7): 843–860. DOI: 10.1016/j.ijimpeng.2004.04.012. [6] KENNEDY R P. A review of procedures for the analysis and design of concrete structures to resist missile impact effect [J]. Nuclear Engineering and Design, 1976, 37(2): 183–203. DOI: 10.1016/0029-5493(76)90015-7. [7] RICHARD P, CHEYREZY M. Composition of reactive powder concretes [J]. Cement and Concrete Research, 1995, 25(7): 1501–1511. DOI: 10.1016/0008-8846(95)00144-2. [8] 张文华, 张云升, 陈振宇. 超高性能混凝土抗缩比钻地弹侵彻试验及数值仿真 [J]. 工程力学, 2018, 35(7): 167–186. DOI: 10.6052/j.issn.1000-4750.2017.03.0237.ZHANG W H, ZHANG Y S, CHEN Z Y. Penetration test and numerical simulation of ultral-high performance concrete with a scaled earth penetrator [J]. Engineering Mechanics, 2018, 35(7): 167–186. DOI: 10.6052/j.issn.1000-4750.2017.03.0237. [9] FARNAM Y, MOHAMMADI S, SHEKARCHI M. Experimental and numerical investigations of low velocity impact behavior of high-performance fiber-reinforced cement based composite [J]. International Journal of Impact Engineering, 2010, 37(2): 220–229. DOI: 10.1016/j.ijimpeng.2009.08.006. [10] 赖建中, 朱耀勇, 徐升, 等. 超高性能水泥基复合材料抗多次侵彻性能研究 [J]. 爆炸与冲击, 2013, 33(6): 601–607. DOI: 10.11883/1001-1455(2013)06-0601-07.LAI J Z, ZHU Y Y, XU S, et al. Resistance of ultra-high-performance cementitious composites to multiple impact penetration [J]. Explosion and Shock Waves, 2013, 33(6): 601–607. DOI: 10.11883/1001-1455(2013)06-0601-07. [11] WU H, FANG Q, CHEN X W, et al. Projectile penetration of ultra-high performance cement based composites at 510–1320 m/s [J]. Construction and Building Materials, 2015, 74: 188–200. DOI: 10.1016/j.conbuildmat.2014.10.041. [12] HOLMQUIST T J, JOHNSON G R, COOK W H. A computational constitutive model for concrete subjected to large strains, high strain rates, and high pressures [C] // The 14th International Symposium on Ballistics. Quebec City: National Defense Research Establishment, Sweden, 1993. [13] RIEDEL W, THOMA K, HIERMAIWE S, et al. Penetration of reinforced concrete by BETA-B-500, numerical analysis using a new macroscopic concrete model for hydrocodes [C] // The 9th International Symposion on the Interaction of the Effects of Munitions with Structures, Berlin Strausberg, 1999. [14] MURRAY Y D. Users manual for LS-DYNA Concrete Material Model 159: FHWA-HRT-05-062 [R]. 2007. [15] MALVAR L J, CRAWFORD J E, WESEVICH J W, et al. A plasticity concrete material model for DYNA3D [J]. International Journal of Impact Engineering, 1997, 19(9–10): 847–873. DOI: 10.1016/S0734-743X(97)00023-7. [16] 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. [17] 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. [18] 崔亚男. 卵形头弹体撞击活性粉末混凝土失效特性实验研究[D]. 天津: 中国民航大学, 2019.CUI Y N. Experimental study on the failure characteristics of reactive powder concrete impacted by ogive-nose projectiles [D]. Tianjin: Civil Aviation University of China, 2019. [19] FENG J, GAO X D, LI J Z, et al. Influence of fiber mixture on impact response of ultra-high-performance hybrid fiber reinforced cementitious composite [J]. Composites Part B: Engineering, 2019, 163(1): 487–496. DOI: 10.1016/j.compositesb.2018.12.141. [20] 任根茂, 吴昊, 方秦, 等. 普通混凝土HJC本构模型参数确定 [J]. 振动与冲击, 2016, 35(18): 9–16. DOI: 10.13465/j.cnki.jvs.2016.18.002.REN G M, WU H, FANG Q, et al. Determinations of HJC constitutive model parameters for normal strength concrete [J]. Journal of Vibration and Shock, 2016, 35(18): 9–16. DOI: 10.13465/j.cnki.jvs.2016.18.002. [21] REN G M, WU H, FANG Q, et al. Triaxial compressive behavior of UHPCC and applications in the projectile impact analyses [J]. Construction and Building Materials, 2016, 113: 1–14. DOI: 10.1016/j.conbuildmat.2016.02.227. [22] AFSHIN N, MOHAMMAD S, MASOUD M, et al. Behavior of steel fiber-reinforced cementitious mortar and high performance concrete in triaxial loading [J]. ACI Materials Journal, 2015, 112(1): 95–103. DOI: 10.14359/51686837. [23] YU Z R, ZHAO H Z, AN M J, et al. Mechanical properties of reactive powder concrete under triaxial compression [J]. Journal of the China Railway Society, 2017(7): 121–126. DOI: 10.3969/j.issn.1001-8360.2017.07.017. [24] FARNAM Y, MOOSAVI M, SHEKARCHI M, et al. Behaviour of slurry infiltrated fiber concrete (SIFCON) under triaxial compression [J]. Cement and Concrete Research, 2010, 40: 1571–1581. DOI: 10.1016/j.cemconres.2010.06.009. [25] SIRIJARONCHAI K, EL-TAWIL S, PARRA-MONTESINOS G. Behavior of high performance fiber reinforced cement composites under multi-axial compressive loading [J]. Cement and Concrete Composites, 2010, 32: 62–72. DOI: 10.1016/j.cemconcomp.2009.09.003. [26] ZHANG K, ZHAO L Y, TAO N, et al. Experimental investigation and multiscale modeling of reactive powder cement pastes subject to triaxial compressive stresses [J]. Construction and Building Materials, 2019, 224: 242–254. DOI: 10.1016/j.conbuildmat.2019.07.049. [27] POLANCO-LORIA M. Improvements to the HJC concrete model in LS-DYNA: STF24 F01286 [R]. Trondheim, Norway, 2002. [28] 侯正纲, 三轴应力状态下混凝土强度研究[D]. 天津: 河北工业大学, 2006.HOU Z G. Research on concrete strength under triaxial stresses [D]. Tianjin: Hebei University of Technology, 2006. [29] 闫东明, 林皋, 徐平. 三向应力状态下混凝土动态强度和变形特性研究 [J]. 工程力学, 2007, 24(3): 58–64. DOI: 10.3969/j.issn.1000-4750.2007.03.010.YAN D M, LIN G, XU P. Dynamic strength and deformation of concrete in triaxial stress state [J]. Engineering Mechanics, 2007, 24(3): 58–64. DOI: 10.3969/j.issn.1000-4750.2007.03.010. [30] 熊益波, 胡永乐, 徐进, 等. 混凝土Johnson-Holmquist模型极限面参数确定 [J]. 兵工学报, 2010, 31(6): 746–751.XIONG Y B, HU Y L, XU J, et al. Determining failure surface parameters of the Johnson-Holmquist concrete constitutive model [J]. Acta Armamentarii, 2010, 31(6): 746–751. [31] 谢和平, 董毓利, 李世平. 不同围压下混凝土受压弹塑性损伤本构模型的研究 [J]. 煤炭学报, 1996, 21(3): 265–270. DOI: 10.3321/j.issn:0253-9993.1996.03.009.XIE H P, DONGY L, LI S P. Study of a constitutive model of elasto-plastic damage of concrete in axial compression test under different pressures [J]. Journal of China Coal Society, 1996, 21(3): 265–270. DOI: 10.3321/j.issn:0253-9993.1996.03.009. [32] KONG X Z, FANG Q, LI Q M, et al. Modified K&C model for cratering and scabbing of concrete slabs under projectile impact [J]. International Journal of Impact Engineering, 2017, 108: 217–228. DOI: 10.1016/j.ijimpeng.2017.02.016. [33] WEERHIJM J, DOORMAAL J C A M V. Tensile failure of concrete at high loading rates: new test data on strength and fracture energy from instrumented spalling tests [J]. International Journal of Impact Engineering, 2007, 34(3): 609–626. DOI: 10.1016/j.ijimpeng.2006.01.005. [34] MAALEJ M, QUEK S T, ZHANG J. Behavior of hybrid-fiber engineered cementitious composites subjected to dynamic tensile loading and projectile impact [J]. Journal of Materials in Civil Engineering, 2005, 17(2): 143–152. DOI: 10.1061/(ASCE)0899-1561(2005)17:2(143). [35] MALVAR L J. Review of static and dynamic properties of steel reinforcing bars [J]. ACI Materials Journal, 1998, 95(5): 609–616. DOI: 10.1016/S0886-7798(98)00088-1. [36] MAO L, BARNETT S, BEGG D, et al. Numerical simulation of ultra high performance fibre reinforced concrete panel subjected to blast loading [J]. International Journal of Impact Engineering, 2014, 64: 91–100. DOI: 10.1016/j.ijimpeng.2013.10.003. [37] MAO L, BARNETT S J. Investigation of toughness of ultra high performance fibre reinforced concrete (UHPFRC) beam under impact loading [J]. International Journal of Impact Engineering, 2017, 99: 26–38. DOI: 10.1016/j.ijimpeng.2016.09.014. [38] PARK J K, KIM S W, KIM D J. Matrix-strength-dependent strain-rate sensitivity of strain-hardening fiber-reinforced cementitious composites under tensile impact [J]. Composites Structure, 2017, 162: 313–324. DOI: 10.1016/j.compstruct.2016.12.022. [39] LIN X S. Numerical simulation of blast responses of ultra-high performance fibre reinforced concrete panels with strain-rate effect [J]. Construction and Building Materials, 2018, 176: 371–382. DOI: 10.1016/j.conbuildmat.2018.05.066. [40] MARSH S P. LASL shock Hugoniot data [M]. California: University of California Press, 1980. [41] 高乐. 活性粉末混凝土高压状态方程研究[D]. 广州: 广州大学, 2011.GAO L. Research on high pressure equation of RPC [D]. Guangzhou: Guangzhou University, 2011. [42] 严少华, 钱七虎, 周早生, 等. 高强混凝土及钢纤维高强混凝土高压状态方程的实验研究 [J]. 解放军理工大学学报(自然科学版), 2000, 1(6): 49–53. DOI: 10.3969/j.issn.1009-3443.2000.06.010.YAN S H, QIAN Q H, ZHOU Z S, et al. Experimental study of equation of state for high-strength concrete and high-strength fiber foncrete [J]. Journal of PLA University of Science and Technology, 2000, 1(6): 49–53. DOI: 10.3969/j.issn.1009-3443.2000.06.010. 期刊类型引用(4)
1. 康耕新，颜海春，张亚栋，刘明君，郝礼楷. 接触爆炸下混凝土墩破坏效应试验与数值模拟. 兵工学报. 2024(01): 144-155 . 
百度学术2. 常慧珠，张冬梅，麻旭东，李世中，贾保峰，赵紫良，何隆. 弹丸侵彻混凝土目标的热固耦合数值分析. 中北大学学报(自然科学版). 2024(02): 185-193 . 百度学术3. 崔冰，王景全，刘加平. UHPC桥梁研究进展与规模化应用技术路径分析. 中国公路学报. 2023(09): 1-19 . 百度学术4. 钟锐，张峰领. 超高性能混凝土、纤维增强高强及高延性混凝土抗侵彻性能比较. 硅酸盐学报. 2021(11): 2423-2434 . 百度学术其他类型引用(9)
-
下载:
百度学术
计量
- 文章访问数: 1023
- HTML全文浏览量: 574
- PDF下载量: 124
- 被引次数: 13