Numerical study on dynamic response and spall damage of filter concrete under impact load
-
摘要: 借鉴局域共振材料的工作机制,通过在混凝土基体中嵌入滤波单元,设计出具有应力波衰减特性的滤波混凝土。通过将滤波混凝土结构简化为质量弹簧力学系统来分析滤波混凝土对应力波的衰减机制。采用数值模拟方法,对比研究了冲击荷载作用下普通混凝土模型和滤波混凝土模型中应力波的传播特性和层裂破坏模式。通过参数分析,研究了滤波单元的材料和几何属性对其储能效果的影响。研究结果表明:滤波单元有效降低了混凝土基体中应力波的传播速度和应力峰值;滤波单元的储能机制有效降低了混凝土基体中的能量;金属球的质量越大,滤波单元的储能效果越好,但弹性层的弹性模量和厚度需要通过适当分析进行设计以实现滤波单元的储能最大化;滤波混凝土基体的局部损伤耗散了荷载中的大量能量,有效降低了结构自由面附近的破坏程度。Abstract: Based on the working mechanism of local resonance materials, a filter concrete with stress wave attenuation characteristics is designed by embedding metal balls wrapped with elastic layer (filter units) in the concrete matrix. First, the stress wave attenuation mechanism of filter concrete is analyzed by simplifying the filter concrete structure into a mass-spring mechanical system. Then, the propagation velocity and peak stress of stress wave in normal concrete model and filter concrete model under impact load are compared by using numerical simulation approach. Through parameter analysis, the influence of the density of metal ball, elastic modulus and thickness of elastic layer on the energy storage of filter units are studied. Finally, the spalling damage patterns of normal concrete model and filter concrete model under impact load are compared. The results show that the filter units can effectively reduce the stress wave propagation velocity and magnitude of peak stress in the concrete matrix. The vibration of the metal balls and the deformation of the elastic layer form a good energy storage mechanism for filter units and effectively reduce the energy exerted by the impact load on the concrete matrix. The larger the mass of the metal balls, the better the energy storage effect of the filter units, while the elastic modulus and thickness of the elastic layer need to be designed through a proper analysis to maximize the energy storage of the filter units. The concrete matrix around the elastic layer has obvious stress concentration and local damage may occur, but the local damage of the filter concrete matrix dissipates a large amount of energy produced by the load, effectively reducing the degree of destruction near the free surface of the structure. Combined with the attenuation effect of the filter units on the stress wave, the filter concrete has achieved good impact resistance.
-
Key words:
- filter concrete /
- local resonance /
- spalling damage /
- stress wave /
- energy storage mechanism
-
电热化学发射技术使用等离子体发生器代替常规点火源,将电能转化为等离子体内能来点燃固体发射药。等离子体与固体发射药相互作用不仅能显著缩短发射药点火延迟时间[1-3],还能消除发射药的温度效应[4],提高燃速,从而有效提高炮口初速和动能[5]。常规发射通常使用指数燃速公式来描述发射药燃烧特性[6-7],C.R.Woodley等[8-9]假定电热化学发射时注入的电能只用于增强发射药的内能,引入了电功率增强因子,给出了与输入电功率相关的发射药燃速公式。M.J.Taylor等[10]进一步研究了辐射、发射药侵蚀和电热裂解对燃速的影响,考虑了引起燃速增强的电功率阈值的影响。李海元等[11-12]、A.Brik等[13]利用改进的密闭爆发器对等离子体增强发射药燃速的特性进行了实验研究; 实验表明等离子体作用下,发射药燃气压力变化迅速,燃气压力梯度增大,点火和燃烧特性发生显著变化。因此,等离子体作用下燃气压力梯度对发射药燃速的影响不可忽略。
本文中,拟通过密闭爆发器实验得到燃气压力曲线,分析等离子体作用下固体发射药的燃烧特性。基于指数燃速公式,结合燃气压力梯度和电功率增强对发射药燃速的影响,得到等离子体作用下固体发射药的瞬态燃速公式。
1. 基本原理
1.1 实验装置
图 1所示为等离子体点火密闭爆发器实验装置结构示意图,主要由脉冲成形网络、底喷式等离子体发生器和密闭爆发器等组成。密闭爆发器为耐高压厚壁圆筒,其中一端装有底喷式等离子体发生器,其主体由聚乙烯毛细管组成,一端采用杆状封闭电极,另一端装有敞开式环状电极,通过爆炸丝将两端电极导通。脉冲成形网络由多个模块并联组成,能通过时序设置产生幅值、脉宽可调的脉冲电流。固体发射药试样由硝化棉纸包裹放置于密闭爆发器内。
实验测试系统主要由传感器和数据采集设备组成。燃气压力采用Kistler6215压力传感器测量,发生器两端的电压和电流分别采用电阻分压器和Rogowski线圈测量。采用JV5200瞬态记录仪对脉冲电源和发生器的电参数以及密闭爆发器内的压力信号进行同步采集记录。
1.2 瞬态燃速理论模型
1.2.1 基本假设
忽略密闭爆发器内固体发射药试样燃烧过程中气相压力和温度随空间的变化,即假设任意时刻密闭爆发器内各位置的气相压力和温度均相同。
固体发射药试样燃烧过程满足几何燃烧定律假设:(1)装药的所有颗粒具有均一的理化性质以及完全相同的几何形状和尺寸;(2)所有药粒表面同时着火;(3)所有药粒具有相同的燃烧环境。
由于气相与密闭爆发器壁面存在换热,需考虑试样燃烧过程中的热损失。本文中采用热损失因数修正发射药火药力的方式来计算发射药试样燃烧过程中的热损失。
1.2.2 气相状态方程
考虑点火药的影响,得到密闭爆发器内的气体状态方程[14]:
p{V0−mρp[1−ψ(t)]−αmψ(t)−αigmig(t)}=T[Rmψ(t)+Rigmig(t)] (1) 式中:p为燃气压力,Pa;V0为密闭爆发器容积,m3;m为装药质量,kg;ρp为固体发射药密度,kg/m3;ψ (t)为t时刻发射药已燃体积分数;α为余容,m3/kg;αig为点火药余容,m3/kg;mig(t)为t时刻已燃点火药质量,kg;T为混合气体温度,K;R为发射药燃气气体常数,J/(mol·K);Rig为点火药燃气气体常数,J/(mol·K)。
1.2.3 气相能量守恒方程
气相能量守恒方程[14]为:
(1−c1)fmψ(t)+figmig(t)+(κ−1)cplEpl(t)=T[Rmψ(t)+Rigmig(t)] (2) 式中:cl为热损失因数;f为发射药火药力,J/kg;fig为点火药火药力,J/kg;κ为燃气比热比;cpl为密闭爆发器中输入等离子体电能的利用率;Epl(t)为t时刻输入等离子体发生器的电能,J。
式(2)通用于常规点火和等离子体点火。常规点火时,Epl(t)=0;等离子体点火时,mig(t)=0。
1.2.4 固体发射药瞬态燃速公式
根据Vieille定律,常规点火时固体发射药燃速公式通常采用指数形式[14]:
u=u1pn1 (3) 式中:u1为燃速系数,n1为燃速指数。
C.R.Woodley等[8]在Vieille定律的基础上引入电功率增强因子,得到等离子体增强发射药燃速公式:
u=u1pn1(1+βePe) (4) 式中:βe为电功率增强因子,MW-1; Pe为输入发生器的电功率,MW。
等离子体点火时,在Woodley燃速公式的基础上,加入燃气压力梯度对发射药燃速的影响[15],得到等离子体作用下固体发射药瞬态燃速公式:
u=u1pn1[1+A(t)n1u21p2n1+1dpdt](1+βePe) (5) 式中:A(t)为发射药燃烧过程中与压力及火焰结构有关的时间函数。
1.2.5 相关系数处理方法
由气体状态方程和能量守恒方程可知:
p{V0−mρp[1−ψ(t)]−αmψ(t)−αigmig(t)}=(1−c1)fmψ(t)+figmig(t)+(κ−1)cplEpl(t) (6) 首先,通过常规点火实验确定热损失因数。常规点火时, Epl(t)=0。发射药试样完全燃完时刻有:ψ(t) =1, p=pm, pm为实验测得的最高燃气压力。由此可求出发射药燃烧过程中的热损失因数:
c1=1+figmig(t)fm−pmf[V0m−α−mig(t)mαig] (7) 其次,通过等离子体点火实验确定电能利用率。等离子体点火时,mig(t)=0。同样发射药试样完全燃完时刻有:ψ(t) =1, p=pm。假定等离子体点火时热损失因数未发生变化,进一步求出等离子体点火过程中的电能利用率:
cpl=pm(V0−αm)−(1−c1)fm(κ−1)Epl(t) (8) 1.2.6 拟合精度
采用均方误差σ来衡量模拟压力曲线与实验压力曲线间的误差:
σ=√n∑i=1[ps(i)−pt(i)]2n (9) 式中:n为压力曲线上选取不同时刻的点数;ps(i)、pt(i)分别为同一时刻模拟和实验测得的压力。
2. 实验结果与分析
对4/7高固体发射药试样进行了2发常规点火和3发等离子体点火密闭爆发器实验。其中密闭爆发器容积为145 cm3,装药量为36.1 g; 实验中采用单个模块放电,电容约为1 300 μF,电感为40 μH。实验参数与结果如表 1所示,其中第3、4发实验发射药置于密闭爆发器中间,第5发实验发射药置于发生器出口处。表 1中, Uc为电容器放电电压;tig为发射药点火延迟时间(定义为爆发器内压力达到20 MPa的时刻);tend为发射药完全燃完时刻,即爆发器内压力达到峰值pm对应的时刻。
表 1 实验参数和结果Table 1. Experimental parameters and results编号 点火方式 Uc/kV tig/ms tend/ms pm/MPa 实验1 2号电点火 1.604 5.05 298 实验2 2号电点火 1.673 4.95 300 实验3 等离子体点火 8.3 0.261 2.93 318 实验4 等离子体点火 10.0 0.197 2.44 319 实验5 等离子体点火 10.1 0.167 2.16 329 实验测得的膛内压力曲线如图 2所示,由图可知等离子体点火大大缩短了发射药试样的点火延迟时间和燃烧过程时间;同时,随着等离子体的注入,密闭爆发器内燃气压力不断增大。
实验测得的膛内压力梯度曲线如图 3所示,由图 3可知在发射药燃烧初始时刻,等离子体作用下的燃气压力梯度远大于常规点火,而且等离子体明显影响了燃气的压力梯度,随着发射药的燃烧,影响程度逐渐降低。
图 4所示为常规点火时实验和数值模拟得到的压力曲线,其中数值模拟压力曲线分别由指数燃速公式和瞬态燃速公式推得。对比可知,指数燃速公式和瞬态燃速公式均有较高的精确度;但是在发射药燃烧后期瞬态燃速公式与实验结果更吻合。
通过分析常规点火和等离子体点火实验相关数据,得到4/7高发射药的电功率增强因子为0.005 MW-1。图 5~7分别给出了第3、4和5发等离子体点火实验时压力与放电功率随时间变化的曲线,其中压力曲线分别通过实验测量、Woodley燃速公式和瞬态燃速公式所得。当加入电功率增强因子后,Woodley燃速公式和瞬态燃速公式得到的压力曲线均符合实验曲线。因此,拟合得到的电功率增强因子适用于4/7高发射药试样。
结合图 3中的压力梯度曲线,分析等离子体输入过程中压力梯度与发射药位置和输入电功率的关系。对比第3、4发实验曲线可知,发射药处于相同位置时,早期燃气压力梯度峰值与输入电功率峰值成正比,且压力梯度变化趋势与输入电功率曲线相似。对比第4、5发实验曲线可知,在相似电功率下,缩短发生器与发射药距离提高了电能利用系数,从而增大了燃气压力梯度。因此,缩短等离子体与发射药距离,提高输入电功率均能获得更理想的等离子体增强效应。
对比图 5~7中各压力曲线可知,在燃烧中间段Woodley燃速公式所得压力曲线略低于实验压力曲线。其原因可能是Woodley燃速公式中仅考虑了等离子体注入期间的电增强效应,未考虑放电结束后的增强效应。因此,脉冲放电结束后Woodley燃速公式所得压力曲线与实验有偏差,而在瞬态燃速公式中通过压力梯度项考虑了放电结束后的燃速增强效应,提高了数值模拟精度。
进一步对比各燃速公式的拟合程度,计算得拟合的压力曲线的均方误差如表 2所示。由表 2可知,随着输入电功率的增加以及发生器与发射药间距的减小,Woodley燃速公式均方误差增大;瞬态燃速公式得到的压力均方误差小于Woodley燃速公式,且受到输入电能、电功率和发生器与发射药间距的影响更小。因此,瞬态燃速公式得到的膛压曲线与实验数据符合情况优于Woodley燃速公式,能更真实地反映等离子体作用下密闭爆发器内燃气压力的变化。
表 2 模拟压力曲线与实验压力曲线间的均方误差Table 2. Mean squared errors between simulated pressure curves and test ones编号 σ/MPa Woodley燃速公式 瞬态燃速公式 实验3 4.325 4.294 实验4 9.312 4.910 实验5 13.506 5.715 3. 结论
利用密闭爆发器实验,研究了等离子体增强4/7高固体发射药的燃速特性,实验表明等离子体作用下燃气压力梯度增大,燃速明显增强。根据实验数据,综合考虑燃气压力梯度和燃速增强因子的影响,拟合了4/7高固体发射药等离子体作用下的瞬态燃速公式,计算得到等离子体注入期间的燃速增强因子为0.005 MW-1。瞬态燃速公式与实验压力曲线符合程度比C.R.Woodley等[8]提出的燃速公式更高,能更精确描述固体发射药在等离子体作用下的燃烧过程。实验表明,在低装填密度时,缩短等离子体与发射药的距离和提高电功率均能获得更明显的等离子体增强效应。
-
表 1 混凝土的材料参数
Table 1. Material parameters of concrete
ρ/(kg·m−3) σc/MPa σt/MPa E/GPa μ a0y/MPa a1y 2 440 34 2.7 30 0.156 8.93 0.625 a2y/MPa−1 a0/MPa a1 a2/MPa−1 a1f a2f/MPa−1 6.437×10−3 11.82 0.446 2.02×10−3 0.442 2.957×10−3 表 2 滤波单元的材料参数
Table 2. Material parameters of a filter unit
材料 ρ/(kg·m−3) E/GPa μ 铅 11400 160 0.44 天然橡胶 900 0.047 0.42 表 3 滤波混凝土模型的几何参数
Table 3. Geometric parameters of the filter concrete model
L/mm D/mm l/mm r/mm T/mm 500 74 75 22 2 -
[1] WU J, ZHOU Y, 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. [2] 汪维. 钢筋混凝土构件在爆炸载荷作用下的毁伤效应及评估方法研究 [D]. 长沙: 国防科学技术大学, 2012.WANG W. Study on damage effects and assessments method of reinforced concrete structural members under blast loading [D]. Changsha, Hunan, China: National University of Defense Technology, 2012. [3] CHEN G, HAO Y F, HAO H. 3D meso-scale modelling of concrete material in spall tests [J]. Materials and Structures, 2015, 48(6): 1887–1899. DOI: 10.1617/s11527-014-0281-z. [4] 郭弦. 冲击作用下混凝土中应力波传播规律研究 [D]. 长沙: 国防科学技术大学, 2010.GUO X. Stress wave propagation in concrete structure under impact loading [D]. Changsha, Hunan, China: National University of Defense Technology, 2010. [5] 巫绪涛, 廖礼. 脆性材料中应力波衰减规律与层裂实验设计的数值模拟 [J]. 爆炸与冲击, 2017, 37(4): 705–711. DOI: 10.11883/1001-1455(2017)04-0705-07.WU X T, LIAO L. Numerical simulation of stress wave attenuation in brittle material and spalling experiment design [J]. Explosion and Shock Waves, 2017, 37(4): 705–711. DOI: 10.11883/1001-1455(2017)04-0705-07. [6] 俞鑫炉, 付应乾, 董新龙, 等. 混凝土一维应力层裂实验的全场DIC分析 [J]. 力学学报, 2019, 51(4): 1064–1072. DOI: 10.6052/0459-1879-19-008.YU X L, FU Y Q, DONG X L, et al. Full field DIC analysis of one-dimensional spall strength for concrete [J]. Chinese Journal of Theoretical and Applied Mechanics, 2019, 51(4): 1064–1072. DOI: 10.6052/0459-1879-19-008. [7] LIU Z Y, ZHANG X X, MAO Y W, et al. Locally resonant sonic materials [J]. Science, 2000, 289(5485): 1734–1736. DOI: 10.1126/science.289.5485.1734. [8] LIU Z Y, CHAN C T, SHENG P. Analytic model of phononic crystals with local resonances [J]. Physical Review B, 2005, 71(1): 014103. DOI: 10.1103/PhysRevB.71.014103. [9] MITCHELL S J, PANDOLFI A, ORTIZ M. Metaconcrete: designed aggregates to enhance dynamic performance [J]. Journal of the Mechanics and Physics of Solids, 2014, 65: 69–81. DOI: 10.1016/j.jmps.2014.01.003. [10] MITCHELL S J, PANDOLFI A, ORTIZ M. Investigation of elastic wave transmission in a metaconcrete slab [J]. Mechanics of Materials, 2015, 91: 295–303. DOI: 10.1016/j.mechmat.2015.08.004. [11] 张恩, 路国运, 杨会伟, 等. 超材料混凝土的带隙特征及对冲击波的衰减效应 [J]. 爆炸与冲击, 2020, 40(6): 063301. DOI: 10.11883/bzycj-2019-0252.ZHANG E, LU G Y, YANG H W, et al. Band gap features of metaconcrete and shock wave attenuation in it [J]. Explosion and Shock Waves, 2020, 40(6): 063301. DOI: 10.11883/bzycj-2019-0252. [12] JIN H X, HAO H, HAO Y F, et al. Predicting the response of locally resonant concrete structure under blast load [J]. Construction and Building Materials, 2020, 252: 118920. DOI: 10.1016/j.conbuildmat.2020.118920. [13] XU C, CHEN W, HAO H, et al. Static mechanical properties and stress wave attenuation of metaconcrete subjected to impulsive loading [J]. Engineering Structures, 2022, 263: 114382. DOI: 10.1016/j.engstruct.2022.114382. [14] OYELADE A, ABIODUN Y, SADIQ M O. Dynamic behaviour of concrete containing aggregate resonant frequency [J]. Journal of Computational Applied Mechanics, 2018, 49(2): 380–385. DOI: 10.22059/JCAMECH.2018.269048.339. [15] HUANG H H, SUN C T, HUANG G L. On the negative effective mass density in acoustic metamaterials [J]. International Journal of Engineering Science, 2009, 47(4): 610–617. DOI: 10.1016/j.ijengsci.2008.12.007. [16] LIU Z Y, CHAN C T, SHENG P. Three-component elastic wave band-gap material [J]. Physical Review B, 2002, 65(16): 165116. DOI: 10.1103/PhysRevB.65.165116. [17] 吴健, 白晓春, 肖勇, 等. 一种多频局域共振型声子晶体板的低频带隙与减振特性 [J]. 物理学报, 2016, 65(6): 064602. DOI: 10.7498/aps.65.064602.WU J, BAI X C, XIAO Y, et al. Low frequency band gaps and vibration reduction properties of a multi-frequency locally resonant phononic plate [J]. Acta Physica Sinica, 2016, 65(6): 064602. DOI: 10.7498/aps.65.064602. [18] EURO C. CEB-FIP model code 1990 [Z]. Lausanne, Switzerland: Thomas TelFord Sevices Ltd., 1993. DOI: 10.1680/ceb-fipmc1990.35430. [19] MALVAR L J, CRAWFORD J E. Dynamic increase factors for concrete [R]. Port Hueneme CA: Naval Facilities Engineering Service Center, 1998. [20] 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. [21] WU H J, ZHANG Q M, HUANG F L, et al. Experimental and numerical investigation on the dynamic tensile strength of concrete [J]. International Journal of Impact Engineering, 2005, 32(1): 605–617. DOI: 10.1016/j.ijimpeng.2005.05.008. 期刊类型引用(15)
1. 薛颖杰,陈智刚,杨芮,崔晋,康彦姝,付建平. 椭圆截面弹体侵彻多层混凝土靶的数值仿真. 中北大学学报(自然科学版). 2025(01): 35-44 . 百度学术
2. 张丁山,张博,付良,徐笑,李鹏飞. 弹体尾部斜锥面形状对侵彻偏转的影响. 高压物理学报. 2024(01): 145-154 . 百度学术
3. 刘泓甫,黄风雷,白志玲,段卓平. 刚性弹体带攻角斜侵彻贯穿混凝土靶板的理论模型. 兵工学报. 2023(08): 2381-2390 . 百度学术
4. 朱超,张晓伟,张庆明,张陶. 弹体斜侵彻双层钢板的结构响应和失效研究. 爆炸与冲击. 2023(09): 140-154 . 本站查看
5. 杨璞,李继承,陈建良,张斌,何丽灵,陈刚. 撞击姿态对构型弹体非正侵彻多层间隔钢靶弹道特性的影响规律. 爆炸与冲击. 2023(09): 123-139 . 本站查看
6. 张丁山,全嘉林,付良,张博,徐笑. 侵彻弹体尖卵形头部形状对偏转力矩的影响. 火炸药学报. 2023(09): 834-839 . 百度学术
7. 姜安邦,李典,李永清,侯海量. 偏转式抗侵彻防护技术研究现状. 材料导报. 2023(24): 99-107 . 百度学术
8. 陈柏翰,沈子楷,邹慧辉,王伟光,王可慧. 硬目标内弹体斜侵彻状态的基本演化特性. 兵工学报. 2023(S1): 59-66 . 百度学术
9. 陈柏翰,王笠镔,邹慧辉,王伟光,王可慧. 自旋对弹体侵彻效应的影响. 兵工学报. 2023(S1): 117-124 . 百度学术
10. 郭松林,高世桥,李泽章,李云彪. 弹引系统攻角侵彻混凝土仿真与试验研究. 兵器装备工程学报. 2022(01): 135-139+205 . 百度学术
11. 成丽蓉,汪德武,贺元吉,赵宏伟,杨杰,邓斌. 弹体斜侵彻多层间隔混凝土薄靶姿态偏转机理研究. 北京理工大学学报. 2022(10): 1009-1016 . 百度学术
12. 李鹏程,张先锋,刘闯,魏海洋,刘均伟,邓宇轩. 攻角和入射角对弹体侵彻混凝土薄靶弹道特性影响规律研究. 爆炸与冲击. 2022(11): 92-104 . 本站查看
13. 廉璞,牟东,青泽,刘军. 侵彻弹药姿态测量技术研究现状及发展. 探测与控制学报. 2021(02): 1-9 . 百度学术
14. 李东伟,刘俞平,王筱锋,王昭明,赵双双. 弹形参数对战斗部斜穿甲姿态偏转影响研究. 兵器装备工程学报. 2021(09): 40-44 . 百度学术
15. 张健东,武海军,李伟,李金柱,皮爱国. 头部带肋板的异形结构弹体斜贯穿混凝土薄靶实验和数值模拟. 高压物理学报. 2021(06): 186-194 . 百度学术
其他类型引用(5)
-