Study on the detonation properties of extrusion cast explosives
-
摘要: 设计了4种以HMX为主炸药的挤注炸药(ECX)炸药配方,采用25 mm和50 mm 2种装药直径的圆筒试验分别对添加含能增塑剂和惰性增塑剂的两种挤注炸药进行做功能力研究。研究结果表明:含能增塑剂FM-1对ECX有显著的能量贡献;在相同主炸药含量(含HMX为88%)下,含FM-1为9.5%的ECX-05较含惰性增塑剂的ECX-02比动能提高22%~23%,优于国外同类型ECX(EX-08-EL)的水平,其加速金属能力与含HMX95.5%的压装PBX炸药(LX-14)相近。Abstract: Extrusion-cast explosive (ECX) is a new type of composite explosive with high energy and low vulnerability. It can be expected that for ECX wide applications will be found in warhead of high performance weapons. In this work, four ECX with HMX were designed. Cylinder expansion tests of explosives with diameter 25 mm and 50 mm were used to identify the power of ECX formulations with energetic and inert plasticizer. It was found that energetic plasticizer FM-1 had great contribution to the energy of ECX. That is, with the same percentage of HMX (88%), the specific energy of ECX-05 with energetic plasticizer was 22%~23% higher than that of ECX-02 with inert plasticizer, and was better than that of the same type ECX (EX-08-EL) produced abroad. The ability of accelerating metal was close to that of LX-14, i. e., pressed PBX with 95.5% HMX.
-
Key words:
- mechanics of explosion /
- detonation energy /
- cylinder expansion test /
- ECX /
- HMX /
- energetic plasticizer
-
在地震或地下爆炸激发的地震波传播研究中, 具有一定压力边界的有限球形空腔可以视为震源, 因此固体介质中的球面应力波传播理论对于研究地震或地下爆炸的震源特征、应力波的传播演化规律等具有重要作用。F.G.Blake[1]和H.L.Selberg[2]基于折合位移势(reduced displacement potential, RDP)给出了理想弹性球面波传播的处理方法。RDP是球面波波动方程的一个解, 球面波传播中涉及到的粒子速度、粒子位移、应力、应变等物理量均可通过RDP给出。H.C.Rodean[3]研究了理想弹性介质中球形发散压缩波的传播, 并讨论了球形发散应力波的特征频率(周期、波长)及其衰减规律。H.C.Rodean的工作对理想弹性球面应力波的传播规律进行了完整的总结。由于波传播和材料的动态力学性能相关, 研究者们基于不同的本构方程(弹塑性[4-8]、弹黏塑性[9-11]、Maxwell体[12]和Kelvin-Voigt体[13]黏弹性、线性和非线性ZWT[14-20]等)对球面应力波的传播特征进行了分析, 进一步完善了球面波传播理论。
在考虑介质的黏弹性条件下, 从震源辐射出的波频率特征和理想弹性条件下波的传播特征不同。为深入研究黏弹性球面应力波的频率响应特性随传播距离的变化特征, 本文中基于线黏弹性球面波Laplace域的理论解[19], 得到不同传播距离处粒子速度、应力、应变等物理量的传递函数, 并以标准线性固体模型为例, 重点讨论粒子速度频率响应函数的传播特征。以弹性半径为0.025 m的空腔爆炸为例, 采用Laplace数值逆变换方法[21]对粒子速度波形的演化进行分析, 给出粒子速度强间断幅值及粒子速度峰值随传播距离变化的衰减曲线。
1. 线黏弹性球面应力波的频率响应函数
无限介质中线黏弹性球面波的Laplace解可按如下方程计算[19]:
˜vr(r,s)=−[1rβ(s)+1r2]s˜φ(r,s) (1) ˜ur(r,s)=−[1rβ(s)+1r2]˜φ(r,s) (2) ˜σr(r,s)=1r[(˜K(s)+43˜G(s))β2(s)+4˜G(s)(1rβ(s)+1r2)]˜φ(r,s) (3) ˜σθ(r,s)=1r[(˜K(s)+23˜G(s))β2(s)−2˜G(s)(1rβ(s)+1r2)]˜φ(r,s) (4) ˜εr(r,s)=[1rβ2(s)+2r2β(s)+2r3]˜φ(r,s) (5) ˜εθ(r,s)=−[1r2β(s)+1r3]˜φ(r,s) (6) 式中:˜K(s)和˜G(s)分别为复体积模量和复剪切模量, ˜vr和˜ur分别为径向粒子速度和粒子位移, ˜σr和˜σθ分别为径向应力和切向应力, ˜εr和˜εθ分别为径向应变和切向应变, s为Laplace变量, β(s)=√ρ0s2˜K(s)+43˜G(s)为Laplace域的波传播系数。基于公式(1)~(6), 能够得到不同物理量的传递函数。假设弹性半径为r0, 则不同半径r(r≥r0)处的传递函数Hg(r, r0, s)(下标g代表公式(1)~(6)左端的各物理量)可以写为:
Hg(r,r0,s)=g(r,s)g(r0,s) (7) 令s=ωi(ω=2πf为圆频率, f为频率, i为虚数单位), 则公式(7)表征不同物理量传递函数对应的频率响应函数。需注意的是, 上面的分析中并没有基于特定的线黏弹性本构模型, 即上述公式具有分析线黏弹性球面发散应力波的普适性。
2. 基于标准线性固体模型的粒子速度频率响应函数传播特征
在地下爆炸或地震研究中, 地运动信号(粒子加速度、速度、位移)是常见的分析对象, 因此本文中仅对粒子速度的频率响应函数进行讨论(理论上粒子加速度、速度、位移的频率响应函数相同)。粒子速度的频率响应函数可以写为:
Hvr(r,r0,ωi)=[1rβ(ωi)+1r2][1r0β(ωi)+1r20]e−β(ωi)(r−r0) (8) 线黏弹性球面应力波的传播系数β(ω i)可以写为:
β(ωi)=α(ω)+k(ω)i (9) 式中:α(ω)为线黏弹性球面应力波衰减因子, k(ω)为线黏弹性球面应力波的波数(k(ω)=ω/c(ω), c(ω)为线黏弹性球面应力波的相速度)。把公式(9)代入公式(8), 则有:
Hvr(r,r0,ωi)=r0r=[α(ω)+k(ω)i+1r][α(ω)+k(ω)i+1r0]e−α(ω)(r−r0)e−k(ω)i(r−r0)=|Hvr(r,r0,ωi)|eiφvr(r,r0,ω) (10) 式中:|Hvr(r, r0, ωi)和φvr(r, r0, ω)分别为粒子速度响应函数的幅频特性和相频特性, 其表达式可以写为:
|Hvr(r,r0,ωi)|=r0r=√(α(ω)+r−1)2+k2(ω)√(α(ω)+r−10)2+k2(ω)e−α(ω)(r−r0) (11) φvr(r,r0,ω)=−k(ω)(r−r0)+Φ(r,ω)−Φ(r0,ω) (12) 公式(12)中Φ(r, ω)满足下式:
tan(Φ(r,ω))=k(ω)α(ω)+r−1 (13) 本文中基于标准线性固体模型展开对粒子速度频率响应函数传播特征的讨论。标准线性固体模型由1个线性弹簧和1个Maxwell体并联而成, 如图 1所示。此模型的复体积模量˜K(s)和复剪切模量˜G(s)可以写为:
˜Ks=13(1−2μ)(E0+E1θ1s1+θ1s)=K0+K1θ1s1+θ1s (14) ˜Gs=12(1+μ)(E0+E1θ1s1+θ1s)=G0+G1θ1s1+θ1s (15) 式中:μ为泊松比, E0和E1分别为线性弹簧和Maxwell体中弹簧的弹性模量, K0和K1分别为对应的体积模量, G0和G1分别为对应的剪切模量, θ1为Maxwell体的松弛因子。令s=ωi=2πfi, 把公式(14)~(15)代入公式(9), 可以给出衰减因子α(f)和波数k(f)的解析式[20]。由文献[20], 当f→∞时(对应高频状态):
α∞=limf→∞α(f)=E12(E0+E1)c∞θ1 (16) k∞=limf→∞k(f)=2πfc∞ (17) 式中:c∞代表高频体波波速, 可以写为:
c∞=√(1−μ)(E0+E1)(1−2μ)(1+μ)ρ0=√(K0+K1)+43(G0+G1)ρ0 (18) 当f→0时(对应低频状态):
α0=limf→0α(f)=2π2 c0 E1θ1E0f2 (19) k0=limf→0k(f)=2πfc0 (20) 式中:c0代表低频体波波速, 可以写为:
c0=√(1−μ)E0(1−2μ)(1+μ)ρ0=√K0+43G0ρ0 (21) 结合公式(11)、(16)~(21), 当f→∞时, 有:
|Hvr(r,r0,2πfi)|f→∞=r0re−α∞(r−r0) (22) 当f→0时, 有:
|Hvr(r,r0,2πfi)|f→0=r20r2 (23) 当公式(22)~(23)满足下式时:
|Hvr(r,r0,2πfi)|f→∞=r0re−α∞(r−r0)=|Hvr(r,r0,2πfi)|f→0=r20r2 (24) 速度响应函数对高频和低频的响应相同, 通过整理公式(24), 有:
lnr−α∞r=lnr0−α∞r0 (25) 定义函数χ(r), 有:
χ(r)=lnr−α∞r (26) 如图 2所示, 函数χ(r)在r∈(0, α−1∞]区间内是递增函数, 在r∈[α−1∞, ∞)区间内是递减函数, 在r=α−1∞处取得极大值。存在下述两种情况:(1)当r0≤α−1∞时, 在r∈[α−1∞, ∞)区间内存在r1满足式(25), 使得高频幅频响应和低频幅频响应相同。在r∈[r0, r1]区间内高频幅频响应高于低频幅频响应, 在r∈(r1, ∞)区间内高频幅频响应低于低频幅度响应;(2)当r0>α−1∞时, 在r∈[r0, ∞)区间内当且仅当r1=r0时满足公式(25), 由此可知在r∈(r0, ∞)区间内高频幅频响应均低于低频幅频响应。
以黄土标准线性固体模型参数(如表 1中所示)为例, 此时α∞≈3.55 m-1(α−1∞≈0.282m)。图 3给出了弹性半径r0=0.025 m(r0≤α−1∞)时粒子速度频率响应函数随传播距离的变化, 可以看出,图 3中曲线符合上述第一种情况的描述。当r>r1=r0>α−1∞时, 其粒子速度频率响应函数和图 3中r=2.5 m时的曲线类似, 高频幅频响应均低于低频幅频响应, 在频率f0(定义为卓越频率)处使得幅频响应取得最大值, 频率f1(定义为上限频率, 高于此频率的幅频响应均低于低频幅频响应, 即信号幅频响应主要集中在上限频率以下)处的幅频响应等于低频幅频响应。
表 1 黄土材料标准线性固体模型参数Table 1. Parameters of the standard linear solid model for loess密度ρ/(kg·m3) 弹性模量E0/GPa 弹性模量E1/GPa 松弛时间θ1/μs 泊松比μ 1 800 1.60 0.33 21.0 0.25 图 4给出了弹性半径r0分别为0.025、0.25、2.5、25和250 m时卓越频率f0和上限频率f1随传播距离的变化。对于特定场地的地下填实爆炸来说, 弹性半径的大小和爆炸当量相关, 爆炸当量越大则弹性半径越大。从图 4中可以看出, f0和f1随着传播距离的增加而减小, 爆炸当量越大, 卓越频率f0和上限频率f1越小。为对比近区和远区粒子速度响应函数的频率特征, 表 2给出了不同弹性半径时0.5和1 000 km处的卓越频率f0和上限频率f1。弹性半径为0.025~25 m时, 在接近1 000 km传播距离处, 此4条曲线在图 4中基本重合, 表 2中的数据也证实了这一点, 其卓越频率f0在2.33~2.48 Hz之间, 上限频率f1在11.32~11.71 Hz之间。可以看出不同弹性半径下, 上限频率f1的变化率比卓越频率f0要小。弹性半径为25 m时对应的爆炸当量是弹性半径为0.025 m时爆炸当量的109倍, 而在1 000 km处粒子速度频率响应函数的卓越频率f0和上限频率f1基本一致, 这说明当波的传播距离据爆心足够远时爆炸当量对地震信号的频谱特性没有显著影响。图 4中弹性半径为250 m时的频率曲线在1 000 km处卓越频率f0为1.19 Hz、上限频率f1为10.09 Hz。从变化趋势上看, 若传播距离继续增加, 弹性半径为250 m时的频率曲线会和弹性半径分别为0.025、0.25、2.5、25 m的4条频率曲线更加接近。
表 2 卓越频率f0和上限频率f1的变化Table 2. Variations of predominant frequency f0 and upper limit frequency f1弹性半径r0/m 0.5 km处的卓越频率f0/Hz 0.5 km处的上限频率f1/Hz 1 000 km处的卓越频率f0/Hz 1 000 km处的上限频率f1/Hz 0.025 110.76 419.65 2.48 11.71 0.25 109.28 414.41 2.48 11.71 2.5 73.82 361.39 2.47 11.70 25.0 26.92 278.37 2.33 11.32 250.0 9.43 184.49 1.19 10.09 3. 粒子速度波形的演化
前面对不同位置粒子速度的频率响应函数特征进行了讨论, 下面就粒子速度的波形演化特征进行数值模拟分析, 主要是研究强间断及粒子速度幅值随传播距离的变化规律。模拟中采用黄土作为示例材料(材料参数如表 1所示), 弹性半径取为0.025 m, 假设空腔边界的压力形式为:
σr0=σr(r0,t)=−(σde−t/t0+σs) (27) 式中:σd+σs为空腔边界上施加压力的强间断幅值, σs为稳态空腔压力, t0为压力的延迟时间。取σd=σs=1 MPa, t0=21 μs。
图 5给出了0.025~4.775 m区域内不同位置处的粒子速度波形。由公式(18)和表 1, 可以给出c∞=1 134.3 m/s。由于空腔边界处施加了强间断压力边界条件, 在空腔边界上粒子速度的强间断幅值vviscor_sd(r0)可以写为:
vviscor_sd(r0)=σd+σsρ0c∞=0.980 8 m/s (28) 粒子速度强间断幅值vviscor_sd(r)可以写为[19]:
vviscor_sd(r)=vviscor_sd(r0)r0re−α∞(r−r0) (29) 通过图 5可以看出, 粒子速度强间断幅值衰减很快, 在图 5(c)中, r≥2.275 m处粒子速度强间断幅值已经几乎观察不到了, 这主要是高频波衰减快造成的。在本文中给出的计算区域内, 含有强间断的粒子速度波形逐渐演化为具有一定上升时间的波形。这里需要指出, 即使r≥2.275 m处的粒子速度强间断从视觉上观察不到, 但理论上是存在的。另外, 通过图 5还可以看出, 当r≥0.775 m时, 粒子速度强间断幅值不再是粒子速度波形的峰值, 这主要是因为不同位置处粒子速度的幅频响应函数对不同频率的波响应不同, 这些位置处高频波的幅频响应不是整个幅频响应曲线的峰值。
图 6给出了黏弹性黄土介质中粒子速度峰值、强间断幅值随传播距离的变化规律。作为对比, 图 6还给出了理想弹性材料中粒子速度峰值随传播距离的变化规律(理想弹性材料中粒子速度强间断幅值和粒子速度峰值相等)。从图 6可以看出, 黏弹性介质中粒子速度幅值的衰减曲线介于理想弹性介质中粒子速度幅值衰减曲线和黏弹性介质中粒子速度强间断幅值衰减曲线之间。在近区, 理想弹性材料中粒子速度强间断幅值衰减曲线(仅有几何扩散效应)和黏弹性介质中粒子速度峰值(或强间断幅值)的衰减曲线的偏离程度较小, 而随着波传播距离的增加, 两者的偏离程度加大。这说明:(1)在近区, 球面波的几何扩散效应对粒子速度峰值(或强间断幅值)的衰减起主导作用;(2)在远区, 介质的黏性对粒子速度峰值(或强间断幅值)的衰减起主导作用。
4. 结论
基于球面发散应力波在Laplace域的理论解, 得到了不同传播距离处粒子速度、应力、应变等物理量的传递函数, 并以标准线性固体模型为例, 重点讨论了粒子速度的频率响应函数的传播特征, 给出了卓越频率和上限频率随传播距离变化的曲线。以弹性半径为0.025 mm的空腔爆炸为例, 采用Laplace数值逆变换方法对粒子速度波形的演化进行了模拟。基于上述分析, 得出如下结论:
(1) 黏弹性球面波的幅频响应函数除了和弹性半径r0及传播距离r有关之外, 还和黏弹性本构引入的本构衰减有关。理想弹性球面波的幅频响应函数仅和弹性半径r0及传播距离r有关;
(2) 在线黏弹性条件下, 当弹性半径r0≤α−1∞时, 随着球面波传播距离的增加, 粒子速度幅频响应函数的高频响应会由高于低频响应状态过渡到低于低频响应状态。当弹性半径r0>α−1∞时, 无论球面波传播距离如何变化, 粒子速度幅频响应函数的高频响应永远低于低频响应;
(3) 黏弹性介质中粒子速度幅值的衰减曲线介于理想弹性介质中粒子速度幅值衰减曲线和黏弹性介质中粒子速度强间断幅值衰减曲线之间。在近区, 球面波的几何扩散效应对粒子速度峰值(或强间断幅值)的衰减起主导作用。在远区, 介质的黏性对粒子速度峰值(或强间断幅值)的衰减起主导作用。
-
计量
- 文章访问数: 2660
- HTML全文浏览量: 141
- PDF下载量: 302
- 被引次数: 0