Analysis of mass abrasion of high-speed penetrator influenced by aggregate in concrete target
-
摘要: 高速侵彻时,弹靶之间发生强烈的局部作用,引起弹体头部发生质量侵蚀,从而影响弹体的侵彻性能。在弹体侵彻过程中,混凝土中的骨料对弹体的质量侵蚀有显著影响。本文通过对高速侵彻混凝土弹体的质量侵蚀实验数据进行分析,进一步分析讨论了混凝土骨料对弹体质量侵蚀的影响。将混凝土靶体视为骨料和砂浆基质混合的二相复合材料,引入混凝土骨料的体积分数χ和骨料的剪切强度τ1代替骨料的莫氏硬度H,给出无量纲骨料修正因子β,建立了修正的弹体质量损失工程模型。模型预测结果与现有的实验数据符合得很好,更准确地表征了混凝土骨料对弹体质量损失的影响。Abstract: The extreme local interaction between the projectile and target will cause mass erosion of the projectile during high-speed penetrating, and then decreases penetration performance of the penetrator. Aggregates in the concrete target will affect mass loss of the projectile obviously. Analysis about the experimental data is conducted to further discuss the effect of concrete aggregate on mass loss of the residual projectile after high speed penetration into concrete target. By assuming the concrete as a two-phase composite composed of mortar and aggregate and introducing the volume fraction and shear strength of aggregate instead of the aggregate Moh’s hardness, a modified engineering model is presented to predict the mass loss of projectile by giving a dimensionless modified factor β affected by aggregate. The modified model is in good agreement with available experimental data and can better characterize the effect of aggregate on the mass abrasion of penetrator into concrete target.
-
Key words:
- high speed penetration /
- mass abrasion /
- aggregate /
- volume fraction /
- shear strength
-
动能弹高速侵彻混凝土时,弹靶界面的高温、高压以及高速摩擦作用使得弹体发生明显的质量侵蚀现象[1]。质量侵蚀引起弹体头形钝化,随着侵彻速度的提高,侵蚀效应变得愈加显著,甚至可能导致弹体发生动态屈曲、弯曲和破裂等结构破坏失效以及诱发弹道倾斜等不稳定性,严重影响弹体的侵彻性能[2-6]。因此,研究高速侵彻混凝土弹体的侵蚀效应对弹体的结构设计和优化具有重要意义。
Forrestal等[2]和Frew等[7]开展了不同弹材的尖卵形弹体高速侵彻混凝土的实验,首次关注到弹体的质量侵蚀现象,完整地记录了弹体的质量损失,并发现实验弹体头部表面有熔化和剧烈刮擦的迹象。另外,何翔等[5]、杨建超等[6]和Mu等[8]也开展了弹体高速侵彻混凝土的实验,发现类似的实验现象。对文献[2]获得的实验数据进行分析,Silling等[9]发现,当撞击速度低于约1 km/s时,侵彻后弹体的质量损失与初始动能存在线性关系,当撞击速度高于约1 km/s时,侵彻后弹体质量损失基本保持不变。Chen等[10-11]进一步综合分析文献[2,7]的实验结果,发现侵蚀后的剩余弹头仍接近尖卵形,同时指出混凝土骨料硬度对弹体的质量损失有显著影响。Jones等[12]假设弹体的质量损失全部来源于侵彻过程中弹靶间剧烈摩擦造成的弹头表面材料的熔化脱落,建立了相应的弹体质量损失的理论模型。He等[13]基于Jones等[12]的工作在其模型中引入了骨料硬度的影响,得到了考虑骨料影响的弹体质量损失预测公式。
本文进一步深入分析混凝土骨料对弹体侵蚀效应的影响,将混凝土视为骨料和砂浆混合的二相复合材料,采用混凝土骨料体积分数χ和骨料剪切强度
τ1 代替骨料莫氏硬度H,引入无量纲骨料修正因子β,建立β修正的弹体质量损失工程预测模型,并基于泰勒展开,推导Silling关系中弹体质量损失与弹体初始动能的线性系数的解析表达式。1. 骨料修正的弹体质量损失模型
通过拟合实验数据,Silling等[9]得出当初始撞击速度
v0≤1km/s 时弹体的质量损失与弹体的初始动能满足线性关系,弹体的相对质量损失δ可表示为δ=ΔMM0=Cv20/2 (1) 式中:线性系数C通过实验数据拟合得到,单位为s2/km2。
Jones等[12]假设弹体的质量损失全部来源于侵彻过程中弹靶间剧烈摩擦造成的弹体头部表面材料熔化脱落,并认为摩擦功全部转化为热,且全部用于熔化弹体表面材料,得到弹体质量损失表达式为:
ΔM=πd2τ0Na0Z4κQ (2) 式中:κ为热功当量;Q为单位质量弹体材料的熔化热,在模型中,Q取纯铁的熔化热,则κQ =1 264.8×103 J/kg;d为弹体直径;Z为计及质量损失的弹体侵彻深度;
τ0 为靶体的剪切强度,根据相应的侵彻深度实验数据拟合得到;对于任意弹体,弹头的无量纲纵截面面积Na 为Na=8d2∫b0ydx (3) 式中:b为弹头长度,x为弹体轴向坐标,y为弹头轮廓线的函数;特殊地,对于尖卵形弹体,式(3a)可简化为:
Na=4ψ2cos−1(1−12ψ)−(2ψ−1)√4ψ−1 (3) 式中:
ψ=s/d 为尖卵形弹头的曲径比,s为尖卵形母线的曲率半径;当ψ=ψ0 (ψ0 初始弹头的曲径比)时,Na=Na0 为初始弹头的无量纲纵截面面积。He等[13]在Jones模型的基础上,考虑骨料的影响,引入无量纲骨料硬度η,得到η修正的弹体质量损失表达式:
ΔM=ηπd2τ0Na0Z4kQ (4) 式中:η=H/H0,H为骨料的莫氏硬度,H0=7为参考硬度,η为骨料莫式硬度的相对值。当骨料为石英石时,η=1;骨料为石灰石时,η=3/7≈0.43。He等[13]在其修正模型中采用Tresca准则估计混凝土的剪切强度为
τ0=fc/√3 [1, 14, 15],并取Z为刚性弹的侵彻深度:Z=2M0πd2N∗0ρtln(1+ρtN∗0v20Sfc) (5) 式中:
ρt 为靶体密度;fc为靶体无约束抗压强度;S是与fc相关的无量纲经验常数,取S=82.6fc−0.544 或S=72.0fc−0.5 [1, 15];弹头形状因子N∗ 为:N∗=8d2∫b0yy′31+y′2dx (6a) 对于尖卵形弹体,上式可简化为:
N∗=13ψ−124ψ2 (6) 当
ψ=ψ0 时,N∗=N∗0 为初始弹头形状因子。将式(5)代入式(4),He等[13]得到η修正的弹体相对质量损失的预测公式:
δ=ητ0Na02ρtN∗0κQln(1+ρtN∗0v20Sfc) (7) 式中:N*0为初始弹头形状因子。
对比式(4)和式(2)可以发现,He等[13]模型仅仅是通过引入骨料莫式硬度的相对值η来修正弹体的质量损失。而莫式硬度H本身只是按照十种矿物之间相对硬度顺序表示的值,并非绝对硬度值。因此,He等[13]考虑骨料影响的弹体质量损失模型较为粗糙。事实上,弹体高速侵彻混凝土时,弹靶之间发生强烈的局部作用,因此考虑局部效应和混凝土靶体的非均匀性,式(2)中混凝土的剪切强度
τ0 受到骨料的影响。式(4)可以视为将混凝土的剪切强度τ0 修正为ητ0 。本文将混凝土靶体视为骨料和砂浆的二相混合材料,根据Hill等人[16]提出的混合物与各相材料之间的关系,引入骨料的剪切强度τ1 和体积分数χ,对混凝土靶体的剪切强度由τ0 修正为τR :1τR=1−χτ0+χτ1 (8) 即
τR=1(1−χ)+χτ0τ1τ0=βτ0 (8) 式中:χ为骨料的体积分数,
τ1 为骨料的剪切强度。由于骨料对于混凝土整体的无约束抗压强度fc影响较小,因此将τ0=fc/√3 视为砂浆的剪切强度。由式(8b)可得到无量纲骨料修正因子β,β与骨料的体积分数χ以及砂浆和骨料的相对强度τ0/τ1 相关。β=1(1−χ)+χτ0τ1=1(1−χ)+χfc√3τ1 (9) 则修正的弹体相对质量损失可表示为:
δR=βτ0Na02ρtN∗0κQln(1+ρtN∗0v20Sfc) (10) 式(10)给出了本文考虑骨料影响的弹体相对质量损失的解析表达式,其中
βτ0 是考虑局部效应的混凝土剪切强度,将混凝土视为骨料和砂浆的混合材料,从而引入骨料对弹体质量损失的影响。图1给出了不同骨料含量χ和骨料剪切强度
τ1 下的β−fc 曲线。其中图1(a)表示强度较高的石英石骨料(取骨料体积分数为40%、50%、60%,剪切强度为20 MPa、30 MPa);图1(b)表示强度较低的石灰石骨料(取骨料体积分数为40%、50%、60%,剪切强度为8 MPa、10 MPa)。从图1可以看出,对于相同体积分数和强度的骨料,骨料修正因子β随着靶体无约束抗压强度fc的增大而减小。当骨料强度大于靶体强度时,β>1;当骨料强度小于靶体强度时,β<1。而He等[13]的骨料修正因子η只与骨料自身的硬度相关,与靶体强度无关,这会导致:当高强度骨料、低强度靶体时,其模型会低估骨料的影响,从而低估弹体的质量损失;而当低强度骨料、高强度靶体时,其模型会高估骨料的影响,从而高估弹体的质量损失。式(1)表明,当弹体撞击速度低于某一临界值时,弹体的质量损失与初始动能成正比。类似Wu等[17]的分析,可给出弹体的上限撞击速度
vc0 以及比例系数C的解析式。将式(10)进行泰勒级数展开,并和式(1)对比即可得到比例系数C为:C=βNa0√3κQS(1−12N∗0ρtSfcv20+...) (11) 若式(11)中的一阶项
12N∗0ρtSfcv20<15 ,即v0<vc0=√2Sfc/(5N∗0ρt) 时,仅保留其常数项,同时考虑忽略了高阶项会使得C值偏大,可取比例系数C的解析式为:C=0.8βNa0√3κQS (12) 2. 实验验证
表1给出了尖卵形弹体高速侵彻混凝土实验的相关参数[2, 6-7]。本节基于各组实验数据对弹体高速侵彻混凝土靶的质量损失模型进行验证。
表 1 实验弹靶参数Table 1. Parameters of targets and projectiles工况 fc/MPa ρt/(kg·m−3) 骨料 H 弹体材料 Yp/MPa ρp/(kg·m−3) m0/kg d/mm L/d ψ0 1-1 13.5 2 000 石英石 7 4340 钢 1 481 7 850 0.064 12.9 6.88 3 1-2 13.5 2 000 0.064 12.9 6.88 4.25 2-1 21.6 2 000 0.064 12.9 6.88 3 2-2 21.6 2 000 0.064 12.9 6.88 4.25 3 62.8 2 300 0.478 20.3 10 3 4 51.0 2 300 1.6 30.5 10 10 5 58.4 2 320 石灰石 3 4340 钢/AerMet100 1 481/1 820 7 850 0.478 20.3 10 3 6 58.4 2 320 4340 钢/AerMet100 1 481/1 820 1.62 30.5 3 7 34.8 2 300 60Si2 Mn/Tc4 1 300/1 030 0.155 14 4.25 8 48.6 2 300 60Si2 Mn/20#钢 1 300/450 0.155 14 4.25 9 61.3 2 300 60Si2 Mn/45#钢 1 300/680 0.155 14 4.25 10 76.4 2 300 60Si2 Mn/35CrMnSi 1 300/1 540 0.155 14 4.25 在文献[2, 6-7]中未直接给出骨料体积分数和剪切强度相关实验参数。根据通常混凝土骨料的体积分数在40%~60%之间,本文取骨料体积分数χ=50%。同时,利用工况1-1和工况6的弹体质量损失实验数据拟合得到石英石和石灰石骨料的剪切强度分别为20和8 MPa,进而应用于其他工况。由式(9)计算得到所有工况的不同fc值下的骨料修正因子η和β,见表2。
表 2 无量纲骨料修正因子η和βTable 2. Dimensionless modified factors η and βfc/MPa 骨料类型 η β fc/MPa 骨料类型 η β fc/MPa 骨料类型 η β 13.5 石英石 1 1.44 62.8 石英石 1 0.71 58.4 石灰石 0.43 0.38 21.6 石英石 1 1.23 34.8 石灰石 0.43 0.57 61.3 石灰石 0.43 0.37 51 石英石 1 0.81 48.6 石灰石 0.43 0.44 76.4 石灰石 0.43 0.31 图2~13给出了本文β修正的弹体质量损失模型与He等[13]的η修正模型的预测结果和实验数据的对比。正如前文分析,在工况1-1和1-2,即高强度骨料、低强度靶体时,He等[13]的
η 修正模型低估了弹体质量损失,如图2和图3;而在工况10,即低强度骨料、高强度靶体时,He等[13]的η 修正模型高估了弹体质量损失,如图13。表3~4给出两组高强度骨料、低强度靶体和低强度骨料、高强度靶体工况(工况1-2和工况10)的模型预测结果误差数据对比。从预测结果来看,本文β修正的弹体质量损失模型与现有的实验数据较符合,预测结果优于He等[13]的η修正模型,更为准确地表征了骨料对弹体质量损失的影响。表 3 工况1-2的模拟结果与实验对比Table 3. Comparison of experimental and simulation results at case 1-2工况1-2弹速/(m·s−1) 345 585 722 900 1 063 He等[13]模型误差 −33% −41% −35% −27% −13% 本文模型误差 −5% −15% −7% 6% 25% 表 4 工况10的模拟结果与实验对比Table 4. Comparison of experimental and simulation results at case 10工况10弹速/(m·s−1) 847 975 1 124 1 165 1 250 1 315 1 376 1 382 He等[13]模型误差 46% 23% 21% 102% 88% 49% 98% 57% 本文模型误差 5% −11% −12% 45% 35% 7% 42% 13% 图14~25给出了在较低速度
v0<vc0=√2Sfc/5N∗0ρt 下,弹体相对质量损失的实验值以及Wu等[17]的线性近似解和本文的β修正的线性近似解的对比。从图中可以看出,通过β修正的线性近似解也更接近实验值。值得说明的是,本文为获得弹体质量损失的解析式忽略了弹体头部形状的钝化,而初始形状的无量纲纵截面面积
Na 相对偏大,会高估弹体的质量损失,因此通过拟合实验数据得到的骨料剪切强度会相对偏低。后续可通过进一步考虑弹体的头形钝化,建立相应的弹体质量损失数值模型。3. 结 语
本文考虑混凝土骨料影响的弹体损失工程模型预测结果与实验数据较吻合,模型的有效性得到验证。将混凝土靶体视为骨料和砂浆的二相复合材料,引入混凝土骨料的体积分数χ和骨料的剪切强度
τ1 代替骨料的莫氏硬度H,给出无量纲骨料修正因子β,更好地表征了混凝土骨料对弹体质量损失的影响,体现了混凝土骨料和砂浆对弹体质量侵蚀的耦合作用,物理意义更加明确。根据β修正的弹体相对质量损失表达式,基于泰勒级数展开,给出了Silling关系中弹体质量损失与弹体初始动能的线性系数的解析式。 -
表 1 实验弹靶参数
Table 1. Parameters of targets and projectiles
工况 fc/MPa ρt/(kg·m−3) 骨料 H 弹体材料 Yp/MPa ρp/(kg·m−3) m0/kg d/mm L/d ψ0 1-1 13.5 2 000 石英石 7 4340 钢 1 481 7 850 0.064 12.9 6.88 3 1-2 13.5 2 000 0.064 12.9 6.88 4.25 2-1 21.6 2 000 0.064 12.9 6.88 3 2-2 21.6 2 000 0.064 12.9 6.88 4.25 3 62.8 2 300 0.478 20.3 10 3 4 51.0 2 300 1.6 30.5 10 10 5 58.4 2 320 石灰石 3 4340 钢/AerMet100 1 481/1 820 7 850 0.478 20.3 10 3 6 58.4 2 320 4340 钢/AerMet100 1 481/1 820 1.62 30.5 3 7 34.8 2 300 60Si2 Mn/Tc4 1 300/1 030 0.155 14 4.25 8 48.6 2 300 60Si2 Mn/20#钢 1 300/450 0.155 14 4.25 9 61.3 2 300 60Si2 Mn/45#钢 1 300/680 0.155 14 4.25 10 76.4 2 300 60Si2 Mn/35CrMnSi 1 300/1 540 0.155 14 4.25 表 2 无量纲骨料修正因子η和β
Table 2. Dimensionless modified factors η and β
fc/MPa 骨料类型 η β fc/MPa 骨料类型 η β fc/MPa 骨料类型 η β 13.5 石英石 1 1.44 62.8 石英石 1 0.71 58.4 石灰石 0.43 0.38 21.6 石英石 1 1.23 34.8 石灰石 0.43 0.57 61.3 石灰石 0.43 0.37 51 石英石 1 0.81 48.6 石灰石 0.43 0.44 76.4 石灰石 0.43 0.31 表 3 工况1-2的模拟结果与实验对比
Table 3. Comparison of experimental and simulation results at case 1-2
工况1-2弹速/(m·s−1) 345 585 722 900 1 063 He等[13]模型误差 −33% −41% −35% −27% −13% 本文模型误差 −5% −15% −7% 6% 25% 表 4 工况10的模拟结果与实验对比
Table 4. Comparison of experimental and simulation results at case 10
工况10弹速/(m·s−1) 847 975 1 124 1 165 1 250 1 315 1 376 1 382 He等[13]模型误差 46% 23% 21% 102% 88% 49% 98% 57% 本文模型误差 5% −11% −12% 45% 35% 7% 42% 13% -
[1] 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. [2] 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. [3] 孔祥振, 方秦, 吴昊, 等. 弹体高速侵彻混凝土靶质量损失和头形钝化的数值模拟研究 [J]. 固体力学学报, 2015, 36(S1): 65–73.KONG Xiangzhen, FANG qin, WU Hao, et al. Numerical study on mass loss and nose-blunting for the projectile during the high-speed penetration on concrete target [J]. Chinse Journal of Solid Mechanics, 2015, 36(S1): 65–73. [4] 刘志林, 王晓鸣, 李文彬, 等. 考虑侵蚀效应的卵形弹丸侵彻混凝土介质模型研究 [J]. 弹道学报, 2017, 29(2): 19–25. doi: 10.3969/j.issn.1004-499X.2017.02.004LIU Zhilin, WANG Xiaoming, LI Wenbin, et al. Model of ogive-nose projectile penetrating concrete target considering effect of mass loss and nose blunting [J]. Journal of Ballistics, 2017, 29(2): 19–25. doi: 10.3969/j.issn.1004-499X.2017.02.004 [5] 何翔, 徐翔云, 孙桂娟, 等. 弹体高速侵彻混凝土效应的实验研究 [J]. 爆炸与冲击, 2010, 30(1): 1–6. DOI: 10.11883/1001-1455(2010)01-0001-06.HE Xiang, XU Xiangyun, SUN Guijuan, et al. Experimental investigation on projectiles’ high-velocity penetration into concrete targets [J]. Explosion and Shock Waves, 2010, 30(1): 1–6. DOI: 10.11883/1001-1455(2010)01-0001-06. [6] 杨建超, 何翔, 金栋梁. 弹体高速侵彻混凝土质量侵蚀特性试验研究 [J]. 防护工程, 2010, 32(1): 6–10.YANG Jianchao, HE Xiang, JIN Dongliang. Experimental investigation into mass loss of projectiles that penetrate concrete targets at high velocity [J]. Protective Engineering, 2010, 32(1): 6–10. [7] FREW D J, HANCHAK S J, GREEN M L, et al. Penetration of concrete targets with ogive-nose steel rods [J]. International Journal of Impact Engineering, 1998, 21(6): 489–497. DOI: 10.1016/S0734-743X(98)00008-6. [8] MU Z C, ZHANG W. An investigation on mass loss of ogival projectiles penetrating concrete targets [J]. International Journal of Impact Engineering, 2011, 38(8−9): 770–778. DOI: 10.1016/j.ijimpeng.2011.04.002. [9] SILLING S A, FORRESTAL M J. Mass loss from abrasion on ogive-nose steel projectiles that penetrate concrete targets [J]. International Journal of Impact Engineering, 2007, 34(11): 1814–1820. DOI: 10.1016/j.ijimpeng.2006.10.008. [10] CHEN X W, HE L L, YANG S Q. Modeling on mass abrasion of kinetic energy penetrator [J]. European Journal of Mechanics A/Solids, 2010, 29(1): 7–17. DOI: 10.1016/j.euromechsol.2009.07.006. [11] 陈小伟, 杨世全, 何丽灵. 动能侵彻弹体的质量侵蚀模型分析 [J]. 力学学报, 2009, 41(5): 739–747. DOI: 10.3321/j.issn:0459-1879.2009.05.017.CHEN Xiaowei, YANG Shiquan, HE Liling. Modeling on mass abrasion of kinetic energy penetrator [J]. Chinese Journal of Theoretical and Applied Mechanics, 2009, 41(5): 739–747. DOI: 10.3321/j.issn:0459-1879.2009.05.017. [12] JONES S E, FOSTER J C, TONESS O A, et al. An estimate for mass loss from high velocity steel penetrators [C] // Pressure Vessels and Piping Conference. New York: ASME, 2002: 227−237. DOI: 10.1115/PVP2002-1149. [13] HE L L, CHEN X W, HE X. Parametric study on mass loss of penetrators [J]. Acta Mechanica Sinica, 2010, 26(4): 585–597. DOI: 10.1007/s10409-010-0341-8. [14] CHEN X W. Dynamics of metallic and reinforced concrete targets subjected to projectile impact [D]. Singapore: Nanyang Technological University, 2003: 1−62. [15] 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. [16] HILL R. Elastic properties of reinforced solids: Some theoretical principles [J]. Journal of the Mechanics and Physics of Solids, 1963, 11(5): 357–372. DOI: 10.1016/0022-5096(63)90036-X. [17] WU H, CHEN X W, HE L L, et al. Stability analyses of the mass abrasive projectile high-speed penetrating into concrete target, Part I: engineering model for the mass loss and nose-blunting of ogive-nosed projectiles [J]. Acta Mechanica Sinica, 2014, 30(6): 933–942. DOI: 10.1007/s10409-014-0090-1. -