Experimental and numerical simulation research on damage effect of jetting projectile charge (JPC) on reinforced concrete wall
-
摘要: 为了满足高侵深和大穿孔的要求，设计一种聚能杆式弹丸（jetting projectile charge, JPC），开展大尺寸钢筋混凝土墙的毁伤效应试验。在此基础上，基于修正参数的K&C（Karagozian & Case）模型进行数值模拟，研究JPC高速侵彻和爆炸冲击波对钢筋混凝土墙的联合破坏作用，分析墙体厚度对破坏效果的影响规律。结果表明，在1.67倍和2.50倍装药直径的炸高条件下，JPC均能够有效贯穿80 cm（6.67倍装药直径）厚的钢筋混凝土墙，形成直径大于6 cm（0.50倍装药直径）的柱状孔洞；聚能装药的多载荷毁伤特性决定了钢筋混凝土墙的破坏结果，爆炸冲击波能够加剧墙体正面开坑和背面崩落的破坏范围；墙体厚度对于墙体正面漏斗坑的直径与深度及内部侵彻孔洞直径均无显著影响；随着墙体厚度增大，背面漏斗坑直径逐渐减小，深度却逐渐增大。
-
关键词:
- 聚能杆式弹丸(JPC) /
- 钢筋混凝土墙 /
- 侵彻 /
- 冲击波 /
- 联合破坏
Abstract: To meet the requirements of a tandem penetrating warhead for high penetration depth and large perforation, a jetting projectile charge (JPC) was designed. The damage test of a large-scale reinforced concrete wall was carried out to analyze the impact of standoff distance on the damaging effect. By constructing a large air domain covering the whole reinforced concrete wall for the transmission of explosion shock wave and JPC, the coupling damage of JPC high-speed penetration and explosion shock wave to the reinforced concrete wall was considered. The damage evolution, strain rate and other parameters of the Karagozian & Case (K&C) model were modified, based on which a numerical model was established to simulate the whole process of the combined damage of JPC and explosion shock wave to the reinforced concrete wall. The reliability of the numerical model was fully verified by comparing the simulation and test results from the failure mode, crater depth and crater diameter of the reinforced concrete wall. On this basis, the combined damage effect of JPC and explosion shock wave on the reinforced concrete wall was further studied, and the influence of wall thickness on the damaging effect was analyzed. The results show that JPC can penetrate the reinforced concrete wall with a thickness of 80 cm (6.67 times of charge diameter) at the standoff distance of 1.67 times and 2.50 times of charge diameter, and form cylindrical holes with a diameter of more than 6 cm (0.50 times of charge diameter). The multi-load damage characteristic of shaped charge determines the damage result of the reinforced concrete wall, and the explosion shock wave can intensify the damage range of the front crater and back crater of the reinforced concrete wall. The wall thickness has no significant effect on the diameter and depth of the crater on the front of the wall and the diameter of the internal penetration hole. With the increase of the wall thickness, the crater diameter on the back gradually decreases and the crater depth on the back gradually increases.-
Key words:
- jetting projectile charge (JPC) /
- reinforced concrete wall /
- penetration /
- shock wave /
- combined damage
-
反舰导弹的战斗部在爆炸后形成空爆冲击波以及破片群，舰船上层建筑防护结构不仅会受到大量破片群以及冲击波的耦合伤害，也会出现破片群单独侵彻作用一定时间之后冲击波再开始作用[1]。破片群密集侵彻结构的破坏机理复杂多样，而纤维增强复合材料因自身力学性能复杂，其破坏机理的分析较金属靶板的破坏机理分析更困难[2]。由于复合材料比强度与比模量较高，应力波在复合材料结构中的传播速度高于在金属结构中的传播速度，将更快地产生应力波叠加效应，从而产生更广泛的破坏范围[3]。在应力波的叠加效应下，靶板的被撞击区域内能量密度与能量持续时间大幅增加，导致靶板在破片侵彻过程的抵抗能力大幅下降，破片的侵彻与穿透能力大幅提高，产生侵彻能力增强效应[4]。
目前开展的实验主要集中于单个破片的侵彻作用，对多破片的实验研究较少。主要是因为开展复合材料破片侵彻的实验研究比较困难，尤其是多破片的侵彻实验，难点在于无法精确控制破片弹道轨迹和冲击速度。本文中，将从数值研究角度对复合材料破片群侵彻作用进行定量研究。针对脆性材料毁伤问题，包括裂纹扩张、分层损伤和高速侵彻等不连续问题时，采用有限元方法等传统方法仍具有局限性，为此，Silling[5]于2000年提出了近场动力学理论（peridynamic, PD），作为一种基于非局部作用的新兴粒子法，近20年受到了广泛关注，尤其是在解决结构不连续性问题上有更好的应用前景。
本文中，基于键的复合材料近场动力学冲击模型，对破片群侵彻作用下纤维增强复合材料结构的损伤进行研究，分析复合材料结构破坏过程及损伤特性，探讨破片群的增强效应，并对破片速度、破片数量和破片群间距对侵彻能力增强效应的影响进行总结。
1. PD基本理论
在近场动力学方法中，物质点xi的作用范围由近场半径δ定义，称为近场范围。近场范围内的所有物质点称为物质点xi的近场域Hx，在近场域范围δ之外的粒子认为不与物质点xi产生相互作用力。物质点xi在其近场域内物质点xj的互相作用下产生变形，物质点xi与xj分别产生位移ui和uj，因此，物质点xi与xj之间的伸长率可以定义为：
sij=|xi+ui−xj−uj|−|xi−xj||xi−xj| (1) Silling[5]最早提出将物质点之间作用力以特殊的力矢量作为力密度函数的方法，以键表示这种成对相互作用的力矢量
f 。为简化表示，将物质点的位移向量与相对位置向量分别表示为η=u(x′,t)−u(x,t) 和ξ=x′−x ，则η+ξ 为变形后的两相对作用点的相对位置向量，得到的运动方程：ρ(x)¨u(x,t)=∫Hf(u′−u,x′−x,t)dH+b(x,t) (2) 式中：
f(u′−u,x′−x,t) 为力密度矢量，它定义为物质点x′ 施加在物质点x 上的单位体积平方的力矢；b(x,t) 为体密度矢量。最初的PD基于键模型以微弹性材料为例建立，这种键可以类比为一种机械弹簧，具有与之相似的属性：
(1)物质点相互作用近场力
f 与物质点之间键的伸长率s呈相关性，定义作用力f 的绝对值与伸长率s的比值为弹簧参数c。(2)当键的伸长率超过临界值
s0 时，认为键已经断裂，即物质点之间的作用力消失且不可恢复。(3)键在受压的情况下不会断裂，伸长率超过临界值后仍具有一定承载能力。
因此，由以上的假设可以将PD物质点本构力函数写为：
f(η,ξ)=(c1s−c2T)η+ξ|η+ξ| (3) 式中：T为成对作用物质点相对于环境温度的变化率平均值，伸长率s定义按照式(1)中的形式。对于各向同性材料而言，式(3)中的参数c1与c2可以通过计算各向同性材料均匀膨胀下的理论均质物体能量密度确定，物体同时也受到温度变化的影响。将PD和传统连续介质力学的能量密度等同c1=c=
18Kπδ4 ，c2=cα，可以推导出c1与c2的值，式中K为体积模量，α为材料的热膨胀系数。近场动力学中物质点的局部损伤定义为其近场域内失效的相互作用点数量与初始总数的加权比例，局部损伤的表达式可以写为：
D(x,t)=1−∫Hμ(t,ξ)dV′∫HdV′ (4) μ(t,ξ)={1s(t′,ξ)＜s0,0≤t′≤t0otherwise (5) 式(4)中局部损伤
D 的变化范围为0～1：损伤值D=1，代表物质点与其近场域内所有点的相互作用力全部终止，该物质点完全破坏失效；而损伤值D=0则表示该物质点近场域内的相互作用力保持完好。物质点的局部损伤可表征物体内可能产生的裂纹形式，物体中一系列微小裂纹随着外力作用而连接在一起，最终形成了宏观上的裂纹。2. 复合材料层合板损伤模式与破坏过程
2.1 破片群侵彻模型
层合板模型几何尺寸为200 mm×100 mm×16 mm，采用的铺层为
[45∘/0∘/−45∘/90∘]2s ，每一层厚度为1 mm，边界条件为四周固支，具体材料参数如表1所示。本文中，采用自行开发的近场动力学程序对复合材料层合板进行建模，模型的离散物质点的宽度Δx =1 mm，长度方向上设置200个物质点，宽度方向上设置100个物质点，厚度方向设置16个物质点，可得物质点体积V=1.0 mm3；近场半径δ =3.015Δx ，时间步大小取为Δt =2×10−8 s，其中近场半径和时间步长已通过验算满足模型计算稳定性要求。表 1 层合板材料性能参数Table 1. Material properties of the laminate参数 含义 数值 单位 E1 x方向弹性模量 125 GPa E2 y方向弹性模量 7.6 GPa E3 z方向弹性模量 7.6 GPa ν12 面内泊松比 0.344 ν13 面外泊松比 0.344 ν23 面外泊松比 0.46 G12 xy平面剪切模量 4.32 GPa G13 xz平面剪切模量 4.32 GPa G23 yz平面剪切模量 3.23 GPa ρ 密度 1678 kg/m3 Xt 纵向拉伸强度 2200 MPa Xc 纵向压缩强度 1100 MPa Yt 横向拉伸强度 50 MPa Yc 横向压缩强度 200 MPa 基于上述模型，本文中研究破片数量、破片间距和破片速度对于板的损伤模式的影响效应[6]。设定破片间距S=1～2 mm，破片速度v=200～1 500 m/s，破片的数量N=2～9，以研究破片数量在侵彻作用中的影响效应，破片的分布形式如图1～2所示。本文中采用球形破片，破片定义为刚性体，半径为8 mm，密度为7 850 kg/m3。
图 1 不同数量破片群排列示意图Figure 1. Schematic diagram of the arrangement of fragment groups with different numbers
图 2 破片群冲击工况示意图Figure 2. Schematic diagram of the impact condition of the fragment group2.2 损伤模式
根据近场动力学理论对于物质点损伤D的定义，物质点的损伤情况可以由一个取值[0,1]的区间表示，D=0表示物质点未损伤，D=1表示物质点与其近场域内的物质点对点力全部消失。为在后处理中可以直观表达模型的损伤情况，用色条表示材料在该损伤模式下的损伤情况，如图3所示。
破片群的侵彻破坏模式与破片间距、速度以及破片数量有关[7]。由表2可知，当破片群间距较小（S=1 mm）和破片速度较低（v=300 m/s）时，破片群侵彻较慢，板在初期形成鼓包变形，随后逐渐在迎弹面形成大变形和基体损伤，背弹面为纤维拉伸破坏。当速度较高时（v=800 m/s），破坏模式为集团冲塞纤维层拉伸破坏，迎弹面沿面内纤维方向发生基体损伤与剪切损伤，在一定范围内产生变形，背弹面主要损伤表现为剪切损伤，损伤范围相对迎弹面更广。当破片速度达到高速时（v=1 200 m/s），破坏模式依然表现为集团冲塞破坏，迎弹面及背弹面的主要损伤为剪切破坏，且背弹面的损伤范围较迎弹面大。
表 2 破片间距S=1mm 时，在破片群侵彻下层合板的损伤模式Table 2. Damage modes of the laminate subjected to fragment group penetration when the fragment spacingS=1mm v/(m·s−1) 基体损伤 剪切损伤 迎弹面 迎弹面 迎弹面 背弹面 300 
800 
1 200 
由表3可知，在破片间距较大的条件下（S=10 mm），当破片速度较低（v=300 m/s）时，由于间距较大，层合板未出现集团冲塞破坏，但在破片冲击位置之间出现了裂纹的连接产生附加效应，在迎弹面出现大范围剪切破坏，损伤范围沿水平以及纤维方向，背弹面的剪切损伤范围相对较小。当破片速度较高（v=800 m/s）时，损伤模式接近于独立穿孔的损伤模式，迎弹面冲击位置范围内产生损伤并在相互之间产生部分断裂情况，背弹面由于破片的高剩余速度产生大面积的拉伸断裂破坏。随着破片速度的提高，达到高速（v=1 200 m/s）时，层合板的损伤模式主要是迎弹面的剪切穿孔，以及在冲击位置附近的部分基体损伤，由于穿透速度过快，损伤范围较小，背弹面产生纤维拉伸破坏。
表 3 破片间距S=10mm 时，在破片群侵彻下层合板的损伤模式Table 3. Damage modes of the laminate subjected to fragment group penetration when the fragment spacingS=10mm v/(m·s−1) 基体损伤 剪切损伤 迎弹面 背弹面 迎弹面 背弹面 300 
800 
1 200 
由以上分析可知，单破片与破片群侵彻的复合材料层合板，其破坏模式有较大区别：单破片侵彻下，层合板的损伤模式以冲击部位的穿孔破坏以及小范围的剪切与基体损伤为主，而破片群侵彻下层合板的破坏模式主要是集团冲塞破坏、独立穿孔、大面积剪切拉伸断裂和基体损伤等。破坏程度与破片数量、间距和速度有关，其中破片间距主要影响集团冲塞破口损伤，随着间距的增大，逐渐转变为独立穿孔-附加效应。破片速度主要影响层合板中损伤的范围比例，随着速度的提高，层合板上下表面的损伤范围逐渐减小，层间损伤程度逐渐增大。
2.3 破片群侵彻过程
图4为破片间距S=1 mm、破片速度v=800 m/s的侵彻工况下，为了便于观察，以层合板的半剖模型表示层合板破坏的具体过程。由图4可知，当t=4 µs时，破片群开始侵彻层合板，破片撞击后层合板将产生层间的压缩波以及面内应力波，压缩波在冲击位置的厚度方向产生局部压缩变形直至破坏失效，而应力波开始逐渐在各破片着弹点之间产生剪切破坏，其宏观表征体现为各破片着弹点附近形成浅坑状凹陷；当t=8 µs时，破片群之间开始产生剪切的延伸连接损伤，随着时间的延长逐渐扩大。由于间距较小，多个破片的动能会汇聚叠加，冲击位置能量密度增加形成大量集团剪切破口。当t=18 µs时，随着破片群侵彻过程中速度的不断衰减，每一层单层板的剪切破坏逐渐转变为拉伸破坏，集团剪切破口进一步形成；t=38 µs时，集团冲塞破口形成，破片群与冲塞块依靠惯性继续飞出。
图 4 破片群侵彻下层合板的破坏过程（S=1 mm, v=800 m/s）Figure 4. Failure process of the laminate subjected to fragment group penetration (S=1 mm, v=800 m/s)图5为破片群S=10 mm、v=800 m/s的侵彻工况下，以板的半剖模型表示层合板破坏的具体过程。由图5可知，当t=4 µs时，破片群开始侵彻层合板，冲击位置由于破片的高速冲击作用产生局部压缩变形，板发生剪切破坏，开始形成弹坑。当t=8 µs时，破片群继续侵彻作用，由于间距较大，冲击位置之间传播的应力波更加均匀分散，破坏较小，迎弹面主要在冲击位置产生独立剪切穿孔。同时，压缩波传播到背弹面产生拉伸破坏。当t=16 µs时，破片群冲击位置之间逐渐有一定损伤出现，体现为破片冲击的附加效应。当t=42 µs时，背弹面进一步拉伸，产生剪切破口，破片群穿透。
图 5 破片群侵彻下层合板的破坏过程（S=10 mm, v=800 m/s）Figure 5. Failure process of the laminate subjected to fragment group penetration (S=10 mm, v=800 m/s)3. 破片群侵彻下复合材料层合板损伤关键因素的影响
层合板在破片群侵彻过程中，多个破片分别在冲击位置造成损伤变形，并在板内产生应力波，当破片穿透层合板的时间
t＞S/C ，其中C为应力波传播速度，使冲击位置的能量密度与能量持续时间大大增加，提高破片的侵彻能力，即应力波的相互叠加会对破片的侵彻产生增强效应[8-11]。当破片穿透层合板的时间t＜S/C ，破片穿透层合板的时间内，应力波未开始相互叠加，此时不会产生破片群的侵彻增强效应，侵彻能力接近于单破片侵彻能力。由此可知，破片间距S和破片速度v是破片群侵彻能力的主要影响因素，前者影响应力波的叠加发生时间，后者影响破片穿透的过程，两者共同决定了破片群的侵彻能力。此外，破片的数量也是影响破片群侵彻能力的重要因素，当破片数量过少，叠加效应不足以产生明显影响。采用控制单一变量方法，开展层合板破片群侵彻关键因素的影响规律研究。为将破片群侵彻能力量化分析，采用破片群穿透后剩余速度最大的破片表征该工况下破片群的侵彻能力。而根据图3中破片群的排列方式以及应力波叠加原理，中心破片将会受到主要的加强作用，其侵彻能力的增强效应最大，所以可以直接采用中心破片的侵彻状态表征相应工况下破片群的侵彻能力。
3.1 破片群数量的影响
保持破片群速度v=800 m/s和破片间距S=1 mm不变，仅改变破片数量获得不同工况下破片群侵彻的剩余速度，如图6所示。单破片侵彻时，破片剩余速度为472 m/s。破片数量为2枚时，两破片的应力波叠加效应开始产生一定影响，但由于破片数量较少，剩余速度为541 m/s。随着破片数量的增加，应力波的叠加效应对于中心破片的影响不断增强，破片剩余速度逐渐增高。
图 6 破片数量对破片群侵彻能力的影响Figure 6. Influence of the fragment number on the penetration ability of the fragment group当破片数量增加至5枚以上时，速度增高幅度相对较小，当破片数量到达9枚时，破片数量对于破片群侵彻能力的影响已经微乎其微。
3.2 破片间距的影响
根据3.1节的分析结果，固定破片数量N=9以及破片速度v=800 m/s，研究破片间距对破片群侵彻能力的影响。图7为不同破片间距下破片群侵彻剩余速度与单破片侵彻剩余速度对比。当S=0 mm时（设置的极端工况，每个破片都紧密地，无间隔地聚集在一起），破片群侵彻剩余速度比单破片侵彻剩余速度高50.4%，随着破片间距的增大，破片群剩余速度不断降低。当S=15 mm时，破片群与单破片侵彻剩余速度基本相等，可以认为破片群效应已退化为独立破片侵彻。这种破片群效应出现的主要原因在于，应力波的叠加所需时间随着破片冲击位置之间距离的增大而不断增长，直至出现破片已经开始穿透层合板而应力波叠加效应仍未作用，导致破片间距对破片群效应的不断减弱直至退化消失。
图 7 破片间距对破片群侵彻能力的影响Figure 7. Influence of the fragment spacing on the penetration ability of the fragment group图8为破片速度v=800 m/s和破片间距S=5 mm的破片群侵彻下，层合板的损伤过程。由图8可知，当t=4 µs时，破片群开始侵彻层合板，冲击位置开始具有法向速度并产生弹坑；当t=12 µs时，随着应力波的传播，迎弹面冲击位置之间开始产生变形损伤，背弹面开始产生局部断裂；当t=38 µs时，由于层间损伤导致背弹面纤维层的拉伸损伤断裂，破口形成，损伤失效材料与破片一起飞出。
图 8 破片群侵彻下层合板的破坏过程（S=5 mm, v=800 m/s）Figure 8. Failure process of the laminate subjected to fragment group penetration (S=5 mm, v=800 m/s)图9为破片速度v=800 m/s和破片间距S=20 mm的破片群侵彻下，层合板的损伤过程。由图9可知，当t=4 µs时，破片群着靶，层合板冲击位置获得法向速度产生弹坑，应力波开始向四周传播；当t=12 µs时，由于破片间距较大，应力波叠加效应已经可以忽略，破片在背弹面形成局部变形拉伸；当t=42 µs时，破片完全穿透层合板，在背弹面形成大范围剪切拉伸损伤，非冲击区域也产生一定的剪切拉伸损伤。
图 9 破片群侵彻下层合板的破坏过程（S=20 mm, v=800 m/s）Figure 9. Failure process of the laminate subjected to fragment group penetration (S=20 mm, v=800 m/s)综上所述，破片间距较小时应力波叠加效应明显，应力波在破片冲击位置以外造成损伤，尤其是层间损伤情况，使得连接处产生拉伸断裂破坏进而失去抗冲击能力。当破片间距较大时，应力波叠加效应影响较小，非冲击区域的位移较小，不足以让背弹面提前产生局部拉伸破坏，背弹面仍具有一定的抗冲击能力，类似于单破片侵彻情况。
3.3 破片冲击速度的影响
固定影响因素破片间距S=3 mm和破片数量N=9，研究破片速度v=300, 800, 1200, 1500 m/s 等4种工况下破片速度与破片群效应的关系。图10为该破片群与相同速度单破片工况下剩余速度对比情况。由图10可知，当v=300 m/s时，破片群剩余速度是单破片的1.78倍。随着破片速度的升高，破片群与单破片的剩余速度差距不断减小，即群效应不断减弱。当v=1 200 m/s时，破片群的侵彻增强效应退化消失为单破片侵彻情况。这是由于在破片间距不变条件下，穿透时间随着破片速度升高而缩短，应力波叠加效应影响时间缩短，即
t＜S/C ，因此破片群侵彻能力增强效应随破片速度的升高而不断退化直至消失。
图 10 不同初始速度下破片群与单破片侵彻能力对比Figure 10. Comparison of penetration ability between fragment group and single fragment at different initial velocities图11为在破片间距S=3 mm、破片速度v=1 200 m/s以及破片数量N=9枚的破片群对层合板的侵彻过程。由图11可知，当t=4 µs时，破片开始侵彻层合板；当t=8 µs时，面内应力波使冲击位置之间的部分产生较大的变形，层间产生剪切破坏；当t=12 µs时，冲击位置之间产生大量剪切损伤，背弹面受到层间压缩波的作用产生剪切损伤，破片最下端已到达层合板底部位置；当t=22 µs时，破片完全穿透层合板。由于破片的高速侵彻，应力波叠加效应不足以影响破片群侵彻能力，损伤模式已退化为单破片模式。
图 11 破片群侵彻下层合板的破坏过程（S=3 mm, v=1 200 m/s）Figure 11. Failure process of the laminate subjected to fragment group penetration (S=3 mm, v=1 200 m/s)综上所述，由于破片速度直接决定了穿透时间，对于破片群侵彻增强效应影响巨大。在相同破片间距工况下，随着破片冲击速度的提高，破片穿透时间缩短，从而应力波叠加效应对破片群侵彻作用有限，导致破片群增强效应退化至消失。
4. 结 论
利用破片冲击近场动力学模型，开展破片群对纤维增强复合材料层合板侵彻的损伤研究，分析其破坏过程，得到的主要结论如下。
(1)破片群侵彻损伤模式多样，包括集团冲塞、独立穿孔、大范围剪切拉伸破坏等，损伤模式与破片数量、破片速度和破片间距相关。
(2)在一定的破片间距与速度下，破片数量的增加对破片群侵彻能力的增强效应明显，5枚以下数量的破片侵彻能力增强效应明显。
(3)破片间距与破片群侵彻能力增强效应负相关，随着破片间距的减小，破片群增强效应提高。破片间距较小时，应力波传播所需时间减少，冲击动能迅速叠加造成层合板快速失去抗冲击能力。
(4)破片速度直接决定了破片穿透时间，在固定间距与数量的情况下，破片速度的提高使得破片更快穿透层合板，应力波叠加效应作用时间不足以影响破片群侵彻能力。
-
图 2 JPC成型过程(t=0, 15, 18, 24, 31, 62, 131 μs)
Figure 2. Forming process of the JPC (t=0, 15, 18, 24, 31, 62, 131 μs)
图 4 JPC成型参数随炸高变化
Figure 4. Forming parameters of the JPC varing with standoff distance
图 7 钢筋混凝土墙的破坏结果（炸高20 cm）
Figure 7. Failure of the reinforced concrete wall (standoff distance is 20 cm)
图 8 钢筋混凝土墙的破坏结果（炸高30 cm）
Figure 8. Failure of the reinforced concrete wall (standoff distance is 30 cm)
图 10 JPC聚能装药破坏钢筋混凝土墙的数值模型
Figure 10. Numerical model of reinforced concrete wall damaged by the JPC
图 11 钢筋混凝土墙破坏的数值模拟结果
Figure 11. Numerical simulation results of reinforced concrete wall failure
图 12 JPC聚能装药破坏钢筋混凝土墙的应力云图
Figure 12. Stress cloud diagram of reinforced concrete wall damaged by the JPC
图 13 钢筋混凝土墙背面破坏过程
Figure 13. Failure process of the back of the reinforced concrete wall
图 14 钢筋混凝土墙破坏的应力云图
Figure 14. Stress cloud diagram of failure of reinforced concrete wall
图 15 减小空气域后JPC聚能装药破坏钢筋混凝土墙的数值模型
Figure 15. Numerical model of reinforced concrete wall damaged by the JPC after reducing air domain
图 16 不同尺寸空气域的数值模拟结果对比（炸高20 cm）
Figure 16. Comparison of numerical simulation results of different air domains (standoff distance is 20 cm)
图 17 不同尺寸空气域的数值模拟结果对比（炸高30 cm）
Figure 17. Comparison of numerical simulation results of different air domains (standoff distance is 30 cm)
图 18 不同厚度钢筋混凝土墙的损伤分布
Figure 18. Damage distribution of the reinforced concrete wall with different thicknesses
图 19 不同厚度钢筋混凝土墙的破坏变化
Figure 19. Damage of reinforced concrete wall with different thicknesses
ρ/（g·cm−3） C4 C5 E0/MPa V 1.225×10−3 0.4 0.4 0.25 1 
下载: 导出CSV
ρ/（g·cm−3） A/GPa B/GPa R1 R2 ω E0/GPa 1.7 630 6.801 4.1 1.3 0.36 10
下载: 导出CSV
ρ/（g·cm−3） A/GPa B/GPa n c m Tm/K Tr/K 8.96 0.09 0.292 0.31 0.025 1.09 1356 293
下载: 导出CSV
表 4 钢筋混凝土墙毁伤试验结果
Table 4. Test results of the reinforced concrete wall
炸高/cm D1/cm D2/cm D3/cm D4/cm H1/cm H2/cm 20 57.5 40.5 6.5 6.2 9.8 10.4 30 39.3 54.8 6.3 6.5 8.8 12.6
下载: 导出CSV
表 5 混凝土材料参数
Table 5. Material parameters of concrete
ρ/（g·cm−3） A0/MPa 泊松比 b1 b2 长度单位换算系数 压力单位换算系数 2.3 28 0.2 0.82 1.03 0.3937 1.45×107
下载: 导出CSV
表 6 钢筋混凝土墙破坏的数值模拟与试验结果对比（炸高20 cm）
Table 6. Comparison between numerical simulation results and test results (standoff distance is 20 cm)
方法 D1/cm D2/cm D3/cm D4/cm H1/cm H2/cm 试验结果 57.5 40.5 6.5 6.2 9.8 10.4 数值模拟 62.4 44.4 5.9 5.6 9.4 11.5 相对误差/% 8.5 9.6 9.2 9.7 4.1 10.6
下载: 导出CSV
表 7 钢筋混凝土墙破坏的数值模拟与试验结果对比（炸高30 cm）
Table 7. Comparison between numerical simulation results and test results (standoff distance is 30 cm)
方法 D1/cm D2/cm D3/cm D4/cm H1/cm H2/cm 试验结果 39.3 54.8 6.3 6.5 8.8 12.6 数值模拟 42.2 55.6 5.8 6.8 8.9 13.1 相对误差/% 7.4 1.5 7.9 4.6 1.1 4.0
下载: 导出CSV
表 8 不同厚度钢筋混凝土墙破坏的数值模拟结果
Table 8. Simulation results of reinforced concrete wall with different thickness
δ/cm D1/cm D2/cm D3/cm D4/cm H1/cm H2/cm H3/cm H3/δ 60 63.2 64.2 5.9 4.6 9.2 8.6 42.2 0.703 70 61.8 52.6 5.9 5.5 9.2 10.6 50.2 0.717 80 62.4 44.4 5.9 5.6 9.4 11.5 59.1 0.738 90 63.4 38.8 5.8 5.7 9.6 13.8 66.6 0.740 100 62.8 27.6 5.8 5.6 9.5 16.2 74.3 0.743
下载: 导出CSV
-
[1] 张雷雷, 黄风雷, 段卓平. 爆炸成型杆式侵彻体对混凝土靶侵彻研究 [J]. 弹箭与制导学报, 2006(SA): 1140–1142. DOI: 10.3969/j.issn.1673-9728.2006.02.361.ZHANG L L, HUANG F L, DUAN Z P. Studies about jetting projectiles charge’s penetration to concrete target [J]. Journal of Projectiles, Rockets, Missiles and Guidance, 2006(SA): 1140–1142. DOI: 10.3969/j.issn.1673-9728.2006.02.361. [2] 李成兵, 沈兆武, 裴明敬. 聚能杆式弹丸侵彻混凝土实验研究 [J]. 振动与冲击, 2007(7): 85–88,183. DOI: 10.3969/j.issn.1000-3835.2007.07.020.LI C B, SHEN Z W, PEI M J. Experimental investigation of rod-shaped projectile penetrating concrete [J]. Journal of Vibration and Shock, 2007(7): 85–88,183. DOI: 10.3969/j.issn.1000-3835.2007.07.020. [3] 顾文彬, 胡亚峰, 刘建青, 等. 基于等质量变壁厚球缺罩聚能杆式射流成型特性研究 [J]. 爆破器材, 2013, 42(4): 14–19. DOI: 10.3969/j.issn.1001-8352.2013.04.004.GU W B, HU Y F, LIU J Q, et al. Characteristic study on jetting projectile charge formation of variable thickness sphere liner based on equi-mass [J]. Explosive Materials, 2013, 42(4): 14–19. DOI: 10.3969/j.issn.1001-8352.2013.04.004. [4] 付恒, 陈智刚. 等壁厚球缺形装药结构优化设计 [J]. 爆破器材, 2016, 45(5): 17–22. DOI: 10.3969/j.issn.1001-8352.2016.05.004.FU H, CHEN Z G. Optimization design of hemispherical charge structure with equal thickness [J]. Explosive Materials, 2016, 45(5): 17–22. DOI: 10.3969/j.issn.1001-8352.2016.05.004. [5] 张钧, 陈智刚, 李小军, 等. 变壁厚球缺罩杆式射流的形成与侵彻性能研究 [J]. 爆破器材, 2016, 45(1): 39–42. DOI: 10.3969/j.issn.1001-8352.2016.01.009.ZHANG J, CHEN Z G, LI X J, et al. Formation and penetration performances of jetting penetrator charge of hemispherical liners with variable thickness [J]. Explosive Materials, 2016, 45(1): 39–42. DOI: 10.3969/j.issn.1001-8352.2016.01.009. [6] 王维占, 冯顺山, 李小军, 等. 局部变壁厚球缺型药型罩侵彻威力性能的数值模拟及试验验证 [J]. 火炸药学报, 2018, 41(4): 420–424. DOI: 10.14077/j.issn.1007-7812.2018.04.018.WANG W Z, FENG S S, LI X J, et al. Numerical simulation and experimental verification of the penetrating power performance of hemispherical liner with locally variable-walled thickness [J]. Chinese Journal of Explosives & Propellants, 2018, 41(4): 420–424. DOI: 10.14077/j.issn.1007-7812.2018.04.018. [7] 张毅, 王志军, 崔斌, 等. 不同结构钛合金罩战斗部侵彻混凝土数值模拟 [J]. 兵器材料科学与工程, 2015, 38(2): 95–98. DOI: 10.14024/j.cnki.1004-244x.2015.02.008.ZHANG Y, WANG Z J, CUI B, et al. Numerical simulation on concrete target penetrated by warhead of different titanium alloy shaped charge [J]. Ordnance Material Science and Engineering, 2015, 38(2): 95–98. DOI: 10.14024/j.cnki.1004-244x.2015.02.008. [8] LI J, HAO H. Numerical study of concrete spall damage to blast loads [J]. International Journal of Impact Engineering, 2014, 68: 41–55. DOI: 10.1016/j.ijimpeng.2014.02.001. [9] WU J, ZHOU Y M, ZHANG R, et al. Numerical simulation of reinforced concrete slab subjected to blast loading and the structural damage assessment [J]. Engineering Failure Analysis, 2020, 118: 104926. DOI: 10.1016/j.engfailanal.2020.104926. [10] HU F, WU H, FANG Q, et al. Numerical simulations of shaped charge jet penetration into concrete-like targets [J]. International Journal of Protective Structures, 2017, 8(2): 237–259. DOI: 10.1177/2041419617706863. [11] LI W B, WANG X M, LI W B. The effect of annular multi-point initiation on the formation and penetration of an explosively formed penetrator [J]. International Journal of Impact Engineering, 2010, 37(4): 414–424. DOI: 10.1016/j.ijimpeng.2009.08.008. [12] GERAMI N D, LIAGHAT G H, MOGHADAS G H R S, et al. Analysis of liner effect on shaped charge penetration into thick concrete targets [J]. Journal of the Brazilian Society of Mechanical Sciences and Engineering, 2017, 39(8): 3189–3201. DOI: 10.1007/s40430-017-0797-6. [13] WANG C, MA T B, NING J G. Experimental investigation of penetration performance of shaped charge into concrete targets [J]. Acta Mechanica Sinica, 2008, 24(3): 345. DOI: 10.1007/s10409-008-0160-3. [14] WU Y C, CRAWFORD J E. Numerical modeling of concrete using a partially associative plasticity model [J]. Journal of Engineering Mechanics, 2015, 141(12): 04015051. DOI: 10.1061/(ASCE)EM.1943-7889.0000952. [15] HONG J, FANG Q, CHEN L, et al. Numerical predictions of concrete slabs under contact explosion by modified K&C material model [J]. Construction and Building Materials, 2017, 155(11): 1013–1024. DOI: 10.1016/j.conbuildmat.2017.08.060. [16] 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. [17] 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. [18] LIN X S, GRAVINA R J. An effective numerical model for reinforced concrete beams strengthened with high performance fibre reinforced cementitious composites [J]. Materials and Structures, 2017, 50(5): 212. DOI: 10.1617/s11527-017-1085-8. [19] LEE M J, PARK G K, KIM S W, et al. Damage characteristics of high-performance fiber-reinforced cement composites panels subjected to projectile impact [J]. International Journal of Mechanical Sciences, 2022, 214: 106919. DOI: 10.1016/j.ijmecsci.2021.106919. [20] CRAWFORD J E, WU Y C, MAGALLANES J M, et al. The importance of shear-dilatancy behaviors in RC columns [J]. International Journal of Protective Structures, 2013, 4(3): 341–377. DOI: 10.1260/2041-4196.4.3.341. [21] LI J, ZHANG Y X. Evaluation of constitutive models of hybrid-fibre engineered cementitious composites under dynamic loadings [J]. Construction and Building Materials, 2012, 30: 149–160. DOI: 10.1016/j.conbuildmat.2011.11.031. [22] 许香照, 马天宝, 郝莉. 大口径聚能装药侵彻厚混凝土靶板的数值模拟及实验研究 [J]. 中国科学:技术科学, 2016, 46(4): 339–348. DOI: 10.1360/N092016-00004.XU X Z, MA T B, HAO L. Experimental and numerical investigation on heavy-caliber shaped-charge penetration in thick concrete target [J]. Scientia Sinica: Technologica, 2016, 46(4): 339–348. DOI: 10.1360/N092016-00004. [23] ABDEL-KADER M. Numerical predictions of the behaviour of plain concrete targets subjected to impact [J]. International Journal of Protective Structures, 2018, 9(3): 313–346. DOI: 10.1177/2041419618759109. [24] 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. [25] MALVAR L J, ROSS C. A review of strain rate effects for concrete in tension [J]. ACI Materials Journal, 1998, 95: 735–739. [26] KONG X Z, FANG Q, ZHANG J H, et al. Numerical prediction of dynamic tensile failure in concrete by a corrected strain-rate dependent nonlocal material model [J]. International Journal of Impact Engineering, 2020, 137: 103445. DOI: 10.1016/j.ijimpeng.2019.103445. [27] Comite Euro-International Du Beton. CEB-FIP model code 1990: design code[M]. London:Thomas Telford Publishing, 1993. [28] 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. [29] LI Q M, TONG D J. Perforation thickness and ballistic limit of concrete target subjected to rigid projectile impact [J]. Journal of Engineering Mechanics, 2003, 129(9): 1083–1091. DOI: 10.1061/(ASCE)0733-9399(2003)129:9(1083. [30] LI P P, BROUWERS H J H, YU Q L. Influence of key design parameters of ultra-high performance fibre reinforced concrete on in-service bullet resistance [J]. International Journal of Impact Engineering, 2020, 136: 103434. DOI: 10.1016/j.ijimpeng.2019.103434. [31] MENG C M, TAN Q H, JIANG Z G, et al. Approximate solutions of finite dynamic spherical cavity-expansion models for penetration into elastically confined concrete targets [J]. International Journal of Impact Engineering, 2018, 114: 182–193. DOI: 10.1016/j.ijimpeng.2018.01.001. [32] CAO Y Y, YU Q L, TANG W H, et al. Numerical investigation on ballistic performance of coarse-aggregated layered UHPFRC [J]. Construction and Building Materials, 2020, 250: 118867. DOI: 10.1016/j.conbuildmat.2020.118867. 期刊类型引用(7)
1. 郝博，刘力维，姜琦. 基于混合算法计算炸药JWL状态方程参数的研究. 工程爆破. 2024(02): 42-48+97 . 
百度学术2. 杨峰，翟红波，苏健军，李尚青，肖洋，刘伟. TNT和温压炸药浅埋爆炸效应差异性研究. 兵器装备工程学报. 2024(05): 160-166 . 百度学术3. 张生，李述涛，王振清，陈叶青. 带壳装药毁伤元威力参数计算方法. 振动与冲击. 2024(13): 25-34 . 百度学术4. 崔浩，郭锐，宋浦，顾晓辉，周昊，邢柏阳. 基于遗传算法辨识炸药JWL状态方程参数的研究. 振动与冲击. 2022(09): 174-180+209 . 百度学术5. 崔浩，郭锐，顾晓辉，宋浦，杨永亮，江琳，俞旸晖. BP神经网络和圆筒能量模型标定炸药的JWL参数. 火炸药学报. 2021(05): 665-673 . 百度学术6. 刘子德，智小琦. DNAN熔铸混合炸药的慢烤试验与模拟. 兵器装备工程学报. 2018(12): 178-181 . 百度学术7. 谢东升，宋巍，王玮，孙滔，黄铁铮，蔡宣明. 高压输电线路覆冰清除爆破参数研究. 中北大学学报(自然科学版). 2018(06): 746-751 . 百度学术其他类型引用(4)
-
下载:
百度学术
计量
- 文章访问数: 529
- HTML全文浏览量: 216
- PDF下载量: 100
- 被引次数: 11