第36卷第2期爆炸与冲击V o l.36,N o.2 2016年3月E X P L O S I O N A N DS HO C K WA V E S M a r.,2016
D O I:10.11883/1001-1455(2016)02-0248-05
叠氮化铜微装药爆轰驱动飞片的数值模拟*
简国祚,曾庆轩,郭俊峰,李兵,李明愉
(北京理工大学爆炸科学与技术国家重点实验室,北京100081)
摘要:为了优化叠氮化铜微装药器件的设计,探究叠氮化铜爆轰驱动飞片的作用原理,根据微装药器件
的实际设计和相关实验,采用A N S Y S/L S-D Y N A流固耦合算法对叠氮化铜爆轰驱动飞片的作用过程作了数
值模拟三具体研究了加速膛长度对飞片的平整性和完整性的影响,分析了微装药的尺寸与飞片速度之间的关
系三研究结果表明:加速膛的长度对飞片的完整性二平整性和速度具有重要影响,在过长的加速膛中飞片飞行
时易发生破碎,加速膛过短飞片的驱动速度不能达到最佳三装药尺寸与飞片速度之间关系密切,装药直径对
飞片速度的前期成长影响不大,但对飞片获得的最大速度却有较为明显的影响;装药的直径大于0.8mm时,增加装药直径并不能使飞片的最大速度明显增加三
关键词:爆炸力学;M E M S引信;流固耦合;L S-D Y N A;叠氮化铜;钛飞片
中图分类号:O383;T Q563国标学科代码:13035文献标志码:A
叠氮化铜异常敏感且威力大,曾经是一种严格禁用的起爆药[1]三但随着M E M S引信向灵巧化二智能化二微型化以及高可靠性二低成本等方向发展[2],将叠氮化铜应用于M E M S引信系统中成为可能三J. T o o n[3]制备了具有纳米多孔铜,将其作为军用雷管炸药的前驱体材料三L.G e r a l d[4]研究了一种薄膜炸药雷管,在基板上沉积铜金属层,然后利用金属层在装置中 原位 反应合成线性尺寸小于1mm的起爆药,进而引爆下一级装药三J.D a n i e l[5]提出可将直径为2mm二厚度为0.5mm的多孔铜与叠氮酸 原位 反应生成一种起爆药三通过测试得出,该起爆药50%发火的临界能量为3.8m J(1m F,87V),具有体积小二成本低二能耗低等优点,可用于子母弹和40mm手榴弹等常规性弹药中三
为了进一步降低叠氮化铜的装药量,提高其使用的安全性,利用飞片聚能原理,叠氮化铜微装药爆轰驱动飞片冲击起爆钝感装药成为一项新的满足M E M S引信爆炸起爆序列需求的起爆技术三其原理是将具有孔隙结构的纳米多孔铜作为一种炸药前驱体材料组装于M E M S引信机构中,通过化学方法将纳米多孔铜 原位 转化成叠氮化物炸药,叠氮化铜爆轰驱动飞片冲击起爆下一级钝感装药三该微装药冲击起爆技术避免了装药和压药潜在的危险性,与M E M S技术集成可实现引信元器件的大批量生产,生产制备过程和爆炸作用后都不产生有害物质,对人与环境友好三
驱动飞片的完整性二平整性和速度是爆轰驱动技术中最关键的技术指标,影响这些指标的因素有微装药的尺寸二密度二飞片的尺寸和材料性质二加速膛的尺寸等等三通过系统的实验考核微装药爆轰驱动飞片的众多影响因素很困难三为此,本文中根据实际设计的M E M S微装药器件和部分实验,运用A N-S Y S/L S-D Y N A数值模拟探究叠氮化铜爆轰驱动飞片的作用过程,分析微装药爆轰驱动飞片作用原理,以期为器件的优化设计提供相关技术支持和理论依据三
1数值模型
叠氮化铜爆轰驱动飞片运动属于接触爆炸作用问题三对于接触爆炸,虽然可以使用L a g r a n g e算法,在叠氮化铜炸药和飞片之间定义接触来考虑相互作用,但由于炸药单元很容易在爆炸过程中发生严重畸变,往往使得计算中断三另外一种算法就是使用流固耦合方法计算炸药对飞片的作用过程中,炸药
*收稿日期:2014-09-05;修回日期:2014-12-02
第一作者:简国祚(1989 ),男,硕士研究生;通讯作者:曾庆轩,z e n g q i n g x u a n@b i t.e d u.c n三
图1数值模型
F i g
.1N u m e r i c a lm o d e l 使用欧拉算法,飞片使用拉格朗日算法,炸药单元不
存在畸变问题[
6]
三本文中采用流固耦合方法计算纳米多孔叠氮化铜炸药对飞片的作用三
数值模型如图1所示,
由叠氮化铜炸药二空气二玻璃纤维聚碳酸酯约束壳体和金属钛飞片4个部分组成三炸药和空气采用欧拉算法,空气域包含整个模型区域,钛飞片和聚碳酸酯约束壳体采用拉格朗日算法,通过*C O N S T R A I N E D _L A G R A N G E _I N _S O L I D 关键字进行耦合三建立二维模型,模拟叠氮化铜驱动飞片过程,采用c m -g -μ
s 单位制建模三 叠氮化铜采用高能炸药燃烧模型和J W L 状态方程三J W L 状态方程能够比较准确地描述爆轰产物
的膨胀驱动过程,其基本形式为[
7
]:p =A 1-ωR 1?è?
??÷V e -R 1V +B 1-ωR 2?è???÷V e -R 2V +ωE V
(1)式中:p 为爆轰产物压力,
V 为爆轰产物相对比容,E 为单位体积内能,状态方程各参数为[8]
:A =410G P a ,B =4.5G P a ,R 1=4.90,R 2=1.3,ω=0.3,ρ=2.
29g /c m 3,爆速D =4.7k m /s ,C -J 爆压p C J =12.55G P a 三 空气采用N U L L 材料模型和L I N E A R _P O L Y N OM I A L 状态方程,
线性多项式状态方程:p ?=C 0+C 1μ+C 2μ2+C 3μ3+(C 4+C 5μ+C 6μ2
)E (2
)式中:μ=1V -1,
p ?为爆轰压力;E 为单位体积内能;V 为相对体积三通常情况下,为方便计算,将空气看作理想气体,设置:C 0=C 1=C 2=C 3=C 6=0,C 4=C 5=
γ-1,γ为气体的质量热容比[9]
三 金属钛飞片和玻璃纤维聚碳酸酯约束壳体采用适用于金属和塑料的各向同性塑料模型[1
0]
,其中金属钛飞片厚度为28μm ,
具体材料参数如表1所示三表1钛飞片和聚碳酸酯的材料模型参数
T a b l e 1P a r a m e t e r s o f t i t a n i u ma n d p o l y
c a r b o n a t e 材料ρ
/(g 四c m -3)E 0/G P a G /G P a K /G P a
ν
钛4.51113.7643.755
94.8030.30聚碳酸酯
1.19
2.34
0.849
3.2560.38
2 结果分析
2.1 飞片的平整性和完整性分析
图2所示为不同时刻叠氮化铜爆轰驱动飞片的压力场图,
数值模型中,选用加速膛长度为0.68mm ,叠氮化铜装药尺寸为?1.0mm?0.6mm ,约束壳体的厚度为0.6mm 三
图2(a )t =0.02μs 时叠氮化铜驱动飞片的压力分布F i g .2(a )P r e s s u r e d i s t r i b u t i o no f f l y e r d r i v e nb y C u (N 3)2a t t =0.02s
图2(b )t =0.06μs 时叠氮化铜驱动飞片的压力分布
F i g .2(b )P r e s s u r e d i s t r i b u t i o no f f l y e r d r i v e nb y C u (N 3)2
a t t =0.06s
9
42 第2期 简国祚,等:叠氮化铜微装药爆轰驱动飞片的数值模拟
图2(c )t =0.11μs 时叠氮化铜驱动飞片的压力分布F i g .2(c )P r e s s u r e d i s t r i b u t i o no f f l y e r d r i v e nb y C u (N 3)2t =0.11μ
s 图2(d )t =0.20μs 时叠氮化铜驱动飞片的压力分布
F i g .2(d )P r e s s u r e d i s t r i b u t i o no f f l y e r d r i v e nb y C u (N 3)2
t =0.20μs
图2(e )t =0.32μs 时叠氮化铜驱动飞片的压力分布F i g .2(e )P r e s s u r e d i s t r i b u t i o no f f l y e r d r i v e nb y C u (N 3)2t =0.32μs 图2(f )t =0.38μs 时叠氮化铜驱动飞片的压力分布
F i g .2(f )P r e s s u r e d i s t r i b u t i o no f f l y e r d r i v e nb y C u (N 3)2
t =0.38μs
在炸药爆轰驱动作用下,
飞片被加速膛剪切成圆弧面状,冲击起爆下一级装药三由于点起爆的方式使得爆轰波从初始阶段就不能形成平面波,加之装药直径很小,爆轰波阵面率先到达飞片的中心部
位,导致这部分飞片最先受到爆轰波的作用,从而形成圆弧形的飞片[
11]
三数值模拟中的飞片形状与实际的飞片形状较为相符,能够较为形象地反应叠氮化铜驱动飞片的作用过程三
图3所示为飞片受到爆轰作用后破损情况的示意图三在建立的多组不同尺寸装药爆轰驱动钛飞片的数值模拟中,运用L S -P R E P O S T 后处理查看运算结果动画时,
剪切的钛飞片在加速膛中飞行到约700μm 时均出现不同程度的破裂,
随着爆轰作用时间的延长,钛飞片破损程度不断加剧直至成为碎片,与实验发现的加速膛超过700μm 后,叠氮化铜爆轰驱动的飞片不能全部起爆H N S -Ⅳ药柱现象相符合三实验设计叠氮化铜M E M S 微起爆器的加速膛长度分别为560二680二700二870和1070μm ,叠氮化铜药片尺寸为?0.8mm?0.6mm ,加速膛尺寸超过700μm 后,微起爆器均不能起爆下一级钝感H N S -Ⅳ炸药柱三因此,可认为由于加速膛过长和飞片破碎,碎片化的金属钛冲击下一级的钝感HN S -Ⅳ药柱时,不能形成完整且尺寸较大的应力作用面,进而不能起爆H N S -Ⅳ炸药柱三
图3爆轰作用后飞片破损示意图
F i g .3M o r p h o l o g i e s o f f l y e r b y d e t o n a t i o no f c o p p
e r a z i d e 052爆 炸 与 冲 击 第36卷
2.2 装药直径对飞片速度的影响
图4不同直径装药驱动飞片的速度时程曲线
F i g .4V e l o c i t y h i s t o r i e s o f f l y e r d r i v e nb y c h a r g
e w i t hd i
f f e r e n t d i a m e t e r s
图4所示为不同装药直径下叠氮化铜驱动飞片
的速度时程图,选用的装药直径分别为0.5二0.6二
0.8和1.0mm ,叠氮化铜装药的厚度均为0.6mm ,
采用相同的材料模型二状态方程二网格划分方式三从图4中可以看出,不同装药直径对飞片速度的前期成长影响不大,但对飞片运行的最大速度却有较为明显的影响;装药的直径大于0.8mm 时,增加装药直径并不能使飞片最大速度明显增加三
根据哈里顿原理可知,
在有限直径药柱内,爆轰波传播与化学反应的能量释放速度和侧向膨胀引起的能量耗散速度之比有关;能量损失随装药直径减小而增大,从而造成波阵面上参数的下降以及化学反应时间的增加
[12
]三通过纳米多孔铜 原位 合成
的叠氮化铜的颗粒度较小,爆轰反应进行的相对较
快,化学反应区内完成反应所经历的时间较短,爆轰波易传播,稳定爆轰的临界直径小三当微装药的厚度一定且直径较小时,微装药爆轰作用的能量有限,从而爆轰驱动飞片的速度不能达到最佳三当微装药的直径较大时,则不利于微装药器件的集成化设计三所以,微装药的直径存在一个最佳值三2.3 装药厚度对飞片速度的影响
图5不同厚度装药爆轰驱动飞片的速度时程曲线
F i g .5V e l o c i t y h i s t o r i e s o f f l y e r d r i v e nb y c h a r g
e w i t hd i
f f e r e n t t h i c k n e s s e s
根据微型起爆装置的实际情况,考察药厚0.3~0.6mm 情况三采用相同的材料模型二状态方程二网格划分方式,装药厚度分别为0.3二0.4二0.5和
0.6mm ,直径均为0.8mm ,
钛飞片的厚度为28μm 三图5所示为不同装药厚度下叠氮化铜驱动飞片的速度时程图三从图5中可以看出,飞片速度与装药厚度关系密切三当装药厚度为0.3mm 时,飞片获得的最大速度约为2247.3m /s
,平均速度为2008.6m /s
;随着装药厚度的增加,飞片的极值速度和平均速度逐渐增大三微装药的厚度也存在一个较为合适的区间值,这是因为叠氮化铜药片太薄时不能稳定爆轰驱动飞片,飞片的驱动速度小,达不到钝感装药的起爆速度阈值,无法冲击起爆下一级钝感装药;叠氮化铜药片太厚时,则存在炸药爆轰作用能量过剩问题,对飞片驱动速度的成长意义不大三
3 结 论
(1
)在数值模拟中飞片的剪切形状与飞片的实际形状较为相符,模型能够揭示叠氮化铜驱动飞片的作用过程三驱动飞片的完整性二平整性和速度与加速膛的长度关系密切,加速膛的长度应限制在
700μm 以下,
且长度不能过短三 (2)不同装药直径对飞片速度的前期成长影响不大,但对飞片达到的最大速度却有较为明显的影响;微装药直径设计为0.8mm 较为合适三
(3
)飞片获得的驱动速度与微装药厚度关系密切,采用越小的装药厚度,飞片获得的极值速度和平均速度越小;随着装药厚度的增加,飞片获得的极值速度和平均速度逐渐增大并趋于稳定三
1
52 第2期 简国祚,等:叠氮化铜微装药爆轰驱动飞片的数值模拟
252爆炸与冲击第36卷参考文献:
[1] P a t n a i kP.Ac o m p r e h e n s i v e g u i d e t o t h e h a z a r d o u s p r o p e r t i e s o f c h e m i c a l s u b s t a n c e s[M].J o h nW i l e y&S o n s I n c,
2007.
[2]牛兰杰,施坤林,赵旭,等.微机电技术在引信中的应用[J].探测与控制学报,2008,30(6):54-59.
N i uL a n j i e,S h i K u n l i n,Z h a oX u,e t a l.A p p l i c a t i o n o fM E M S o n f u z e s[J].J o u r n a l o f D e t e c t i o n&C o n t r o l,2008, 30(6):54-59.
[3] T o o nJ.E x p l o s i v e so nac h i p u n i q u e p o r o u sc o p p e r s t r u c t u r e se n a b l en e w g e n e r a t i o no fm i l i t a r y m i c r o-d e t o n a t o r s
[J].R e s e a r c hH o r i z o n s,2007,25(1):22-23.
[4] G e r a l dL.I n t e g r a t e d t h i n f i l me x p l o s i v em a c r o-d e t o n a t o r:U S7322294B1[P].2008.
[5] D a n i e l J.M E M Sm i c r o-d e t o n a t o r b a s e d f u z i n g[C]?T h e53r da n n u a l f u z e c o n f e r e n c e.2009.
[6]时党勇,李裕春,张胜民.基于A N S Y S/L S-D Y N A8.1进行显示动力分析[M].北京:清华大学出版社,2005.
[7]石少卿,康建功.A N S Y S/L S-D Y N A在爆炸与冲击领域内的工程应用[M].北京:中国建筑工业出版社,2011:2-31.
[8]曾庆轩,简国祚,李兵,等.叠氮化铜J W L状态方程参数拟合[J].火工品,2014(6):23-26.
Z e n g Q i n g x u a n,J i a nG u o z u o,L i B i n g,e t a l.T h e f i t t e d p a r a m e t e r s o f J W Le q u a t i o no f s t a t e f o r c o p p e r a z i d e[J].
I n i t i a t o r s&P y r o t e c h n i c s,2014(6):23-26.
[9]简国祚,曾庆轩,李明愉,等.基于A N S Y S/L S-D Y N A流固耦合算法的爆炸硬化仿真分析[C]?第十届全国爆炸力
学学术会议论文集.2014:208-212.
[10] L i v e r m o r eS o f t w a r eT e c h n o l o g y C o r p o r a t i o n.L S-D Y N Ak e y w o r du s e r sm a n u a l[M].V e r s i o n971.L i v e r m o r e,
C a l i f o r n i a,2011.
[11]覃文志,龙新平,何碧,等.B N C P驱动飞片的实验测试及数值模拟[J].火工品,2012,1(1):15-17.
Q i n W e n z h i,L o n g X i n p i n g,H e B i,e t a l.E x p e r i m e n t s a n d n u m e r i c a l s i m u l a t i o n o f f l y e r s d r i v e n b y B N C P[J].I n-
i t i a t o r s&P y r o t e c h n i c s,2012,1(1):15-17.
[12]张宝平,张庆明,黄风雷.爆轰物理学[M].北京:兵器工业出版社,2001:180-190.
S i m u l a t i o no f f l y e r s d r i v e nb y d e t o n a t i o no f c o p p e r a z i d e
J i a nG u o z u o,Z e n g Q i n g x u a n,G u o J u n f e n g,L i B i n g,L iM i n g y u
(S t a t eK e y L a b o r a t o r y o f E x p l o s i o nS c i e n c e a n dT e c h n o l o g y,
B e i j i n g I n s t i t u t e o f T e c h n o l o g y,B e i j i n g100081,
C h i n a)
A b s t r a c t:A i mt o o p t i m i z e t h e d e s i g n o fM E M S c o p p e r a z i d e f u z e a n d i n v e s t g a t e t h em e c h a n i s mu n d e r-l y i n g t h e p r o c e s s o f t h e c o p p e r a z i d e e x p l o s i v e-d r i v e n f l y e r p l a t e.A c c o r d i n g t o t h e a c t u a l d e s i g n o f t h e m i c r o-c h a r g e f u z e a n d i t s r e l a t e de x p e r i m e n t s,t h e p r o c e s so f t h ec o p p e ra z i d ee x p l o s i v e-d r i v e nf l y e r p l a t ew a s s i m u l a t e d a d o p t i n g t h e f l u i d-s o l i d c o u p l i n g a l g o r i t h mi nL S-D Y N A p r o g r a m.T h e i n f l u e n c e s o f t h eb a r r e l s l e n g t ho n t h e f l y e r sv e l o c i t y a n d i n t e g r i t y w e r e s t u d i e da n d t h e r e l a t i o n s h i p b e t w e e n t h em i c r o-c h a r g e s i z ea n dt h e f l y e r'sv e l o c i t y w e r ed i s c u s s e d.O u r r e s e a r c hr e s u l t s i n d i c a t e t h a t t h e
b a r r e l s l e n g t hh a s am a j o r i m p a
c t o n t h e f l y e r s v e l o c i t y a n
d i n t
e g r i t y.I t i s
f o u n d t h a t,w h e n i t i s a c-
c e l e r a t e
d i na l o n g b a r r
e l,t h e
f l y e r i s l i k e l y t o b em o r e f r a
g i l e a n d c a n n o t a c
h
i e v em a x i m a l d r i v i n g v e-l o c i t y.T h e s i z e o f t h e m i c r o-c h a r g e i su n i q u e l y r e l a t e d w i t ht h e f l y e r sv e l o c i t y i nt h a t t h e f l y e r s m a x i m u mv e l o c i t y i s s i g n i f i c a n t l y a f f e c t e db y t h e c h a r g e sd i a m e t e r.W i t ht h e i n c r e a s eo f t h e t h i c k-n e s s o f t h ec h a r g e,t h ea v e r a g ev e l o c i t y a n dt h e m a x i m a l v e l o c i t y w e r er a i s e d g r a d u a l l y.W h e nt h e c h a r g e d i a m e t e r i s a b o v e0.8mm,i t s i n f l u e n c e o n t h e f l y e r sm a x i m a l v e l o c i t y i s n o t r e m a r k a b l e. K e y w o r d s:m e c h a n i c so f e x p l o s i o n;M E M Sf u z e;f l u i d-s t r u c t u r e i n t e r a c t i o n;L S-D Y N A;c o p p e r a z-i d e;t i t a n i u mf l y e r
(责任编辑王易难)