JournalofSolidRocketTechnology
固体火箭技术
Vol.42No.42019
①
固体推进剂装药枪击试验数值模拟与试验
李军强1,申依欣2,刘 春1,张佳钰3,张崇民1
付小龙4,樊学忠1,张国防4
3.西安北方秦川集团有限公司,西安 710032;4.陕西师范大学化学化工学院,西安 710062)
摘要:固体推进剂的枪击试验是评价固体推进剂及装药在受到枪击的情况下,是否发生燃烧、爆炸及爆轰等剧烈反应Dyna的计算结果。结果表明,以点火增长模型、Johnson⁃Cook本构模型和Gruneisen状态方程分别赋予推进剂、子弹及推进剂壳体,计算结果表明,推进剂发生明显的燃烧反应,且推进剂燃烧在壳体内部产生的压力显著增大。采用12.7mm枪击试验验证的试验结果表明,推进剂装药在12.7mm子弹以约850m/s速度撞击下,推进剂发生燃烧,此结果与模拟结果相一致。
关键词:应用化学;固体推进剂;装药;数值模拟;枪击
中图分类号:V435+.21 文献标识码:A 文章编号:1006⁃2793(2019)04⁃0451⁃05DOI:10.7673/j.issn.1006⁃2793.2019.04.005
(1.西安近代化学研究所,西安 710065;2.陕西省西安市阎良区95960部队,西安 710089;
的重要试验方法。采用LS⁃Dyna有限元仿真软件研究了枪击试验推进剂的反应过程,并采用12.7mm枪击试验验证了LS⁃
Numericalsimulationandexperimentonbulletimpact
characteristicsofasolidpropellantcharge
LIJunqiang1,SHENYixin2,LIUChun1,ZHANGJiayu3,ZHANGChongmin1,
2.Unit95960inYanLiangDistrict,Xi’ancity,ShaanxiProvince,Xi’an 710089,China;4.SchoolofChemistry&ChemicalEngineering,ShaanxiNormalUniversity,Xi’an 710062,China)
lookingatthelevelofresponses,suchascombustion,explosion,anddetonationduringtheimpactprocess.Inthispaper,thereactionprocessofthebulletimpacttestofpropellanthasbeenstudiedbyLs⁃Dynafiniteelementsimulationsoftware,andthecalculationre⁃sultswereverifiedbythe12.7mmbulletimpacttest.Thepropellant,bulletandpropellantshellaregivenbytheignitiongrowthmodel,Johnson⁃CookconstitutivemodelandGruneisenstateequation,respectively.Thecalculationresultsshowthatthepropellantbulletimpacttestshowthatthepropellantchargeundergoexplosionatavelocityofabout850m/swith12.7mmbullet,andthere⁃sultisconsistentwiththesimulationresults.
Keywords:appliedchemistry;solidpropellants;charge;numericalsimulation;bulletimpact
hasobviousignitingreactionandthedetonationpressureissignificantlyincreased.Theexperimentresultsverifiedbythe12.7mm
Abstract:Thebulletimpacttestofthesolidpropellantsisanimportantexperimentalevaluatingmethodoftheirvulnerabilityby
3.Xi’anBeifangQinchuanGroupCo.Ltd,Xi’an 710032,China;(1.Xi’anModernChemistryResearchInstitute,Xi’an 710065,China;
FUXiaolong4,FANXuezhong1,ZHANGGuofang4
0 引言
穿透推进剂及装药的过程[1-2],可评价固体推进剂及
收稿日期:2019⁃01⁃21;修回日期:2019⁃04⁃01。
装药在受到不同类型枪击的情况下,是否发生燃烧、爆炸及爆轰等剧烈反应,从而判断固体推进剂是否满足钝感需求的重要试验方法[3-5]。然而,由于枪击试验
固体推进剂的枪击试验是模拟在战场环境中子弹
①
作者简介:李军强(1979—),男,硕士,从事固体推进剂研究。E⁃mail:llijq@sohu.com通讯作者:付小龙(1982—),男,博士,从事固体推进剂研究。E⁃mail:fuxiaolong204@163.com
—451—
2019年8月固体火箭技术第42卷
现象较为复杂、反应剧烈,难以用试验过程进行详细描述,仅可用试验结果进行简单判断,且试验准备的周期较长、样品制备复杂、危险性大
[6-8]
2 计算方法
2.1 算法
仿真计算采用Lagrange算法。
模拟的方法,在枪击试验前进行计算,可有效预测反应趋势,揭示试验燃烧、爆炸等反应过程,大幅减少试验数量,节省人力物力成本。
枪击敏感性的理论模拟与试验研究[9-11]。国防科技大学庄建华等
[12]
。因此,采用数值
近年来,国内外学者更加重视炸药、推进剂等对
12.7mm根据高速变形条件下功⁃热转化理
论,建立了瞬态热传导的理论模型,并进行了数值模拟。研究结果表明,子弹射击固体火箭发动机,能使固体推进剂内形成高温热点,并进一步分析了该过程中推进剂高温热点的生成和发展过程。李小柱等[13]采用壁厚模拟发动机及压力传感器测试枪击试验过程中发动机的内压曲线,研究了钢、铝和玻璃钢壳体及3种推进剂对发动机枪击安全性的影响。研究结果表明,枪击子弹及壳体碎片具有一定的动能和温度,是引燃发动机装药的主因。子弹及壳体碎片的温度越高,发动机装药越容易被击燃。中国工程物理研究院代晓淦等
[14]
采用枪击试验对4种不同尺寸PBX
炸药进行了试验,并用有限元模拟方法对炸药样品的枪击试验进行了数值模拟。结果表明,枪击试验中随着PBX炸药长度的增加,炸药的反应程度显著增强,其模拟计算中撞击后的子弹速度与试验测试结果基本相符。
mm枪击试验推进剂在不同直径枪击下的反应过程本文采用LS⁃Dyna有限元仿真软件研究了12.7
同时,验证了相同试验条件下推进剂12.7mm枪击试,验。本文的研究可为固体推进剂枪击感度的研究提供理论指导与试验方法借鉴。
1 试验方法
烧弹射击带壳体的固体推进剂试样12.7mm枪击试验的原理:使用,在子弹的高速撞12.7mm穿甲燃
击及摩擦等因素作用下,子弹的部分机械能可转化为热能,固体推进剂在撞击及热量作用下可能发生热分解、燃烧或者爆炸反应。经综合分析,确定固体推进剂装药的响应等级法603.2。12.7mm枪击试验中所使用的固体推进枪击试验方法采用。
GJB772A—1997剂尺寸方
为
ϕ45mm×180mm,试验壳体采用45#钢,内径为ϕ50mm,壁厚4mm,子弹速度为(850±20)m/s,发射装置口径为12.7mm,射击距离为20~30m,弹丸质量为46g,枪击试验示意图见图1。—452—
图1 12.7mm枪击试验示意图
Fig.1 Diagramof12.7mmbulletimpacttest
2.2 2.2.1 本构模型推进剂冲击起爆的本构关系包含推进剂本构模型
2个方面:反应速率方程和材料本构模型。
推进剂的冲击起爆的反应速率是数学模型的核心
部分。Lee⁃Tarver提出的点火增长模型计算精确,在冲击起爆、破片撞击等领域应用广泛,包含1个三项式反应速率方程和 ddtF2个JWL状态方程。三项式反应速率方程:=I(1-F)b(ρρ
-1-a)x+G0
(Ignition)1(1-F)cFdpy(1)
+G (Growth)式中 F为反应份数2(1-F)eFgpz;t为时间 (Completion)
;ρp为压力;a为临界压缩度0为初始密度;ρ为当前密度;;y为燃烧项压强指数;b和c为点火和燃烧项的燃耗阶数;I和x控制热点数量,是冲击强度和作用时间的函数;d和G;g和G1为热点早期反应的控制参数制参数。
2为高压反应速率的控当F≥FMXIG时,点火速率为零;当F≥FMXGR时,增
长速率为零;当F≤FMNGR时,完成速率为零。其中,FMXIG为初始反应份数,FMXGR为第一阶段F的最大值,FMNGR为第二阶段F最小值。
Tarver反应速率方程需要与状态方程联合使用点火增长模型联合使用的状态方程是JWL,与状态
Lee⁃
2019年8月李军强,等:固体推进剂装药枪击试验数值模拟与试验第4期
方程,可表示为p=A(1-
ωωωE)eR1V+B(1-)eR2V+) (2)R1VR2VV
3 结果与讨论
3.1 物理模型及网格划分
为研究子弹撞击过程对固体推进剂的影响,了解子弹穿透样品的作用过程,采用LS⁃DYNA程序对枪击试验进行三维计算。计算过程中,子弹与壳体采用刚性体模型,固体推进剂弹塑性模型,枪击试验物理模型见图2,计算材料的主要性能参数见表1。对子弹、推进剂3)。模型采用mm⁃mg⁃μs的单位制,压力单位为Mbar。药柱实体模型统一采用映射方式进行网格划分(见图
式中 p为压力;V为相对体积;E为内能;A、B、R1、R2、ω分别为表征推进剂金属加速特性的常数,由圆筒试验标定。
推进剂的材料模型采用弹塑性流体模型,其本构2.2.2 金属本构模型
关系采用剪切模量、密度和屈服应力等参数来描述。
子弹采用刚体模型。壳体采用Gruneisen状态方程和计算和断裂准则Johnson⁃CookJohnson⁃Cook。本构模型包含两部分本构模型描述,材料为其中,流动应力的计算方程为
:45#流动应力的钢。σy=(A+B(εp
)n
)(1+C1nε
·p
·ε式中 A为静态屈服应力;B为应变硬化模量0
)(1-T∗m)(3)
;C为应变率相关系数;n为应变硬化指数;m为温度相关系数;σy为流体屈服应力;εp为有效塑性应变;·
εp/·
ε0为等效塑性应变率;T为温度;T∗为特征温度;T∗T=room)/(Tmelt-Troom),Tmelt为材料的熔点,Troom为室温(。
T-断裂处的应变定义为 εf=[D1+D2exp(D3σ
p-
)]×
式中 σ=[1.[15s+D41n(ε
)][1σ
∗
+D5T∗]
(4)
ijsij为常数]0.5,s为偏应力;D1、D2、D3、D4、D5的行为特性Gruneisen。
,其压缩材料压力为
状态方程用于模拟金属材料在高压下
ρ0C2μ[1+(1-
γ0
=)μ-
α p[1-(S22
μ2]
+1-1)μ-S2μμ2+1-S3μμ
3
+](5)
对于膨胀材料(γ1
0+αμ)Ep=,ρ其状态方程为
0C2μ+(γ0+αμ)E
(6)
式(5)、式(6)中,ρ线斜率的截距;γ0为材料密度;C为式(7)μs0为Gruneisen系数;α为对γ0⁃μ的一p曲阶体积修正量,μ=ρ/ρ能。
0-1;E为材料单位初始体积内μ表述为
s⁃μp曲线是材料冲击绝热线,即Hugoniot线,其μs=C+S1μp+S2(
μμp2 μs
)μp+S3(
μμp3s
)μp
(7)
式中μs为冲击波速度;μp为质点速度;S1s⁃μp曲线斜率的系数。
、S2、S3为图2 12.7mm枪击试验物理模型
Fig.2 Physicalmodelof12.7mmbulletimpacttest
表1 子弹、壳体与固体推进剂材料参数
Table1 Thematerialparametersofbullet,shell
andpropellant
材料密度/强度/杨氏模子弹(g7.83/cm3)MPa量/GPa泊松比2352940.26壳体7.852352120.26推进剂
1.70
0.5
10
0.48
图3 12.7mm枪击试验物理模型网格
Fig.3 Physicalmodelgridof12.7mmbulletimpacttest
3.2 枪击试验仿真结果
为考察带壳体的推进剂在枪击作用下的安全性能,分别研究了不同时刻推进剂的反应度和压力云图,
—453—
2019年8月固体火箭技术第42卷
见图4、图5。反应度达到1时代表推进剂起爆,此时对应的压力值便为推进剂的爆轰压力值。
燃烧反应程度加大。由图5可知,子弹以850m/s的速度冲击带壳体的推进剂,子弹打穿壳体后速度逐渐FringeLevelsFringeLevels1.000e+006.878e-018.883e-016.152e-017.765e-015.426e-016.648e-014.700e-015.531e-013.974e-014.413e-013.248e-013.296e-012.521e-012.179e-011.795e-011.061e-011.069e-01-5.609e-033.432e-02-1.173e-01-3.829e-02(a)403μs (b)407μs
FringeLevels1.000e+00Fringe8.812e-011.000e+00Levels7.625e-018.859e-016.437e-017.718e-015.249e-016.577e-014.062e-015.437e-012.874e-014.296e-011.686e-013.155e-014.987e-022.014e-01-6.889e-02-1.877e-01-2.676e-028.733e-02-1.408e-01(c)411图μs4 不同时刻反应度云图
(d)415μs
Fig.4 Responsivecloudmapatdifferenttimes
FringeLevelsFringeLevels8.181e-023.170e-016.581e-022.770e-014.981e-022.370e-013.381e-021.970e-011.781e-021.569e-011.815e-031.169e-01-1.418e-027.689e-02-3.018e-023.687e-02-4.618e-02-3.147e-03-6.218e-02-4.317e-02-7.818e-02-8.319e-02(a)403μs (b)407μs
FringeLevelsFringeLevels3.556e-014.438e-013.118e-013.908e-012.681e-013.379e-012.244e-012.849e-011.806e-012.320e-011.369e-011.790e-019.316e-021.261e-014.942e-027.310e-025.690e-032.015e-02-3.804e-02-3.280e-02-8.178e-02-8.575e-02(c)411图5 μs 不同时刻压力云图
(d)415μs
Fig.5 Pressurecloudmapatdifferenttimes
由图4可知,推进剂在407μs时刻反应度达到1,此时推进剂发生燃烧,随着时间的增长,压力值增大,—454—
减低至596.8m/s,但随着推进剂的燃烧,在燃烧压力的作用下,子弹速度增加至639.4m/s,在子弹自身速度和燃烧压力的共同作用下,子弹从壳体穿出。子弹速度⁃时间关系见图6。
900850800
)s/m750(/v700650600550
350
375
400
425450
475
500
t/µs
图6 子弹速度⁃时间关系图Fig.6 Relationshipbetweenvelocityand
timeofthebullet
3.3 枪击试验结果
使用12.7mm穿甲燃烧弹射击带壳体的固体推进剂试样,试验装置与试验结果见图7~图10。
图7 12.7mm枪击试验装置
Fig.7 Testdeviceof12.7mmbulletimpacttest
图8 12.7mm枪击试验样品
Fig.8 Testsampleof12.7mmbulletimpacttest
2019年8月李军强,等:固体推进剂装药枪击试验数值模拟与试验
Solids,2001,49(1):27⁃68.
第4期
dynamicloading[J].JournaloftheMechanicsandPhysicsof[3] SimhaCHM,BlessSJ,BedfordA.Computationalmodelingof
thepenetrationresponseofahigh⁃purityceramic[J].Interna⁃tionalJournalofImpactEngineering,2002,27(1):65⁃86.ofImpactEngineering,2003,28(1):13⁃31.
[4] LimCT,ShimVPW,NgYH.Finite⁃elementmodelingof
theballisticimpactoffabricarmor[J].InternationalJournal[5] GrujicicM,PanduranganB,ZhaoCL,etal.Hypervelocity
图9 枪击后样品抛出结果
Fig.9 Throwingresultsofsamplesaftershooting
图10 枪击后样品燃烧结果
Fig.10 Explosionresultsofsamplesaftershooting
由试验结果可知,推进剂装药在12.7mm子弹以约850m/s速度撞击下,推进剂发生燃烧,被抛出2m左右,壳体穿孔,见证板完好。试验结果表明,在子弹的高速撞击及摩擦等因素作用下,子弹的机械能迅速转化为热能,作用于固体推进剂时,使推进剂受热发生分解直至燃烧。
4 结论
试验的条件下推进剂装药的响应结果(1)采用LS⁃Dyna计算了推进剂在。由模拟结果可12.7mm枪击
知Gruneisen,分别以较能反应实际的结果状态方程赋予了推进剂点火增长模型、Johnson⁃Cook本构模型和,在推进剂反应度达到、子弹以及推进剂壳体1的条件下,推进剂发生明显的燃烧反应(2)采用12.7,mm推进剂装药在枪击试验验证了,且压力显著增大12.7mmLS⁃Dyna。
子弹以约
的计
850算结果m/s。速度撞击下结果表明结果相一致。
,推进剂发生燃烧,此结果与模拟参考文献:
[1] Ben⁃dorvancesinG,DubinskianalyticalA,ElperinofT.Ballisticplatepenetrationimpact:recentdynamics
ad⁃
[2] [J].AppliedZavattieriMechanicsmodelingReviews,2005,58(6):355⁃371.almodelofPD,RaghuramceramicmicrostructuresPV,EspinosasubjectedHD.Atocomputation⁃
multi⁃axial
impactthermalresistanceprotectionofsystemsreinforced[J].carbon⁃carbonAppliedSurface/carbon⁃foam[6] 2006,252:5035⁃5050.
Science,formanceGrujicicM,PanduranganB,ZecevicU,etal.Ballisticper⁃
trixing(AP)compositeofaluminaandhybrid/S⁃2lightweightglassfiber⁃reinforcedarmoragainstpolymer⁃ma⁃armorpierc⁃[7] inIannucciMaterialsandnon⁃APStructures,2007,3:287⁃312.projectiles[J].MultidisciplineModel
Engineering,2006,32(6):1013⁃1043.
compositeL.underProgressiveimpact[failureJ].InternationalmodellingJournalofwovenofcarbon
Impact
[8] ZaeramixedR,Sánchez⁃Sáez[9] posites:PartceramicParkH,KimA,2000,31:823⁃833.
/metalarmoursS.Modelingsubjectedofthetoadhesiveimpact[J].Com⁃layerin
compositeH.Damageresistanceofsinglelapadhesive[10] nalKimofstrainH,KeuneImpactjointsEngineering,2010,37(2):177⁃184.bytransverseimpact[J].InternationalJour⁃rates[J].JN.CompressiveJournalofMaterialsstrengthScience,oficeand2007,impact42:
[11] 2802⁃2806.
ShazlyiceunderM,PrakashuniaxialV,LerchB.Highstrain⁃ratebehaviorof
[12] 庄建华Solidsand值模拟,[J].毛佳Structures,2009,46:1499⁃1515.
compression[J].InternationalJournalof固体火箭技术,张为华,等.,2009,32(4):422⁃426.固体火箭发动机枪击过程数
ZHUANGJianhua,MAOgunshotJia,ZHANGprocessofWeihua,etsolidrocketal.Numerical
motor[J].
[13] 李小柱JournalsimulationofSolidoftheRocketTechnology,2009,32(4):422⁃426.究[J].弹箭与制导学报,裴养卫.固体火箭发动机枪击低易损性试验研
,2000,20(2):39⁃42.
LIstudyXiaozhu,ofofsolidPEIProjectiles,rocketYangwei.engineLowRockets,undervulnerabilityMissilespoppingandconditionexperimentGuidance,[J].[14] 2000,20(2):39⁃42.
Journal代晓淦炸药响应规律研究,申春迎,吕子剑[J].含能材料,等.枪击试验中不同尺寸,2008,16(4):432⁃435.PBX⁃2
DAIpropertiesXiaogan,fordifferentSHENChunying,sizePBX⁃2LVexplosivesZijian,etinal.bulletReactionim⁃16(4):432⁃435.
pacttest[J].ChineseJournalofEnergeticMaterials,2008,
(编辑:卓飞)
—455—
因篇幅问题不能全部显示,请点此查看更多更全内容