Rheological properties of Composition B in slow cook-off process
-
摘要: 为进一步探究熔铸炸药在烤燃过程中内部各物理场的变化情况,以B炸药为研究对象,完整地建立了基于Bingham流体模型的B炸药黏度计算模型并应用于慢速烤燃的数值模拟。通过数值模拟得到了B炸药在整个升温过程中上中下3个内部测点处的温度变化曲线并以烤燃试验加以验证,观察了弹体内部温度场与对流场的变化特点。结果表明:升温速率为1 ℃/min时,B炸药相变后逐渐开始流动,内部的温度场分布也随之改变,炸药出现自热反应与最终响应的区域都在弹体上部;升温速率为0.055 ℃/min时,炸药相变后内部很长时间内仍表现出类固相温度场的分布特点,当炸药出现自热反应后,才逐渐开始流动,温度场也逐渐转变为典型的液相温度场,炸药最终响应点在弹体上部,但最早出现自热反应的区域在弹体中心。Abstract: In order to investigate the changes of the internal physical fields of melt-castable explosives in cook-off, Composition B was chosen as the object. A complete viscosity model of Composition B based on the Bingham flow model was first established, and then applied in the numerical simulations of slow cook-off. In this way, the temperature curves of three inner measuring pointsthat located in the upper, middle and lower respectively were obtained and further testified with cook-off experimental measurements. Moreover, the variation of the inner temperature field in the whole process was observed as well. The results showed that when the heating rate was 1 ℃/min, the viscosity flow of Composition B appeared soon after the phase change, and the inner temperature field changed with that. The self-heating and ignition occurred in the upside area of the shell. But when the heating rate was 0.055 ℃/min, the inner temperature field was still like a solid phase after the phase change was completely done for a long time, and the viscosity flow appeared after the self-heating started, the inner temperature field gradually began to change like a liquid phase just at that time. The ignition area was in the upside of the shell too, but the self-heating area was in the middle of the shell. The contradictory points of view in previous studies can be preliminarily explained by this model.
-
Key words:
- slow cook-off test /
- melt-castable explosives /
- viscosity /
- Bingham fluid /
- responses
-
烤燃试验是分析与评估含能材料热刺激响应等级的重要手段,近年来一直是国内外钝感弹药研究的热门话题。按照装药方式,炸药可大致分为压装炸药、熔铸炸药、浇注炸药与塑态挤注炸药四大类[1],其中相比于其他类型的炸药,熔铸炸药的烤燃过程最为复杂[2]。在达到炸药的相变点之前,熔铸炸药表现为固态,其内部的温度分布特征与其他固体炸药类似。但当炸药开始相变后,具有了一定的流动性,就不再是固态,物质与能量的输运方式发生了改变,这是熔铸炸药与固态炸药烤燃特征出现差异的根本原因。
B炸药(TNT与RDX的质量分数为40%和60%)作为典型的熔铸炸药,相变前后物理场的主要区别是自然对流的出现。但对于相变后B炸药内部的流动情况,出现了两种截然不同的观点。一种观点认为B炸药熔化后,内部并未形成对流,如Maienschein等[3]在烤燃试验中通过内部温度测点观察了B炸药在1 ℃/h的升温速率下内部的温度变化,认为炸药内部并未形成对流,自热反应最早出现在炸药的几何中心。Hobbs等[2]也认为B炸药中的对流情况不能用简单的Boussinesq近似假设解释,炸药内悬浮的大量黑索金(RDX)颗粒将很大程度上阻碍对流的发生。但是另一种观点认为炸药发生相变后内部存在对流,如:Glascoe等[4]首先通过高速摄影观察到B炸药响应时刻点火位置出现在弹体上端,认为对流传热对B炸药的烤燃响应过程有重要影响;Nichols等[5]也认为B炸药相变后内部存在对流,对流传热会导致炸药熔化后内部温度最高点出现在上部,响应点也在上部。
B炸药发生相变后,除了温度场外,影响对流形成的最大因素是炸药自身的黏度。在早期的研究工作中,大多将B炸药的黏度近似认为是一个常数或是简单的随温度变化的关系式[5-8]。Nunez等[9]发现液态B炸药存在剪切稀化现象,是一种非牛顿流体。Zerkle等[10]通过落球黏度计分析了液态B炸药的黏度,验证了其是剪切稀化流体,黏度会随应变率变化而变化,Davis等[11-12]在此基础上发现液态B炸药与TNT/HMX熔铸炸药的流动特征符合Bingham流体特性。
本文在前人研究的基础上,分析归纳液态B炸药的特性,建立在热刺激全过程中B炸药的黏度模型,在此基础上,利用慢烤数值模拟对B炸药烤燃过程特别是熔化后的温度场分布、流动特性以及点火位置进行分析,并通过试验加以验证,以期为进一步认识熔铸炸药的热刺激响应机制和该类弹箭的热安全使用提供理论和试验依据。
1. B炸药的黏度模型
Bingham流体是一类特殊的非牛顿流体,在低应变率下表现为类固性,但当所受应力超过某一临界值时,开始屈服并出现黏性流动[13],其所受应力与黏度的数学表达式为[14]:
μ={∞τ<τyμapp+τy/|γ|τ≥τy (1) 式中:μ为黏度;τ为应力;τy是屈服应力阈值;μapp是其出现黏性流动时的表观黏度;γ是应变率。当τ超过τy时,B炸药屈服并出现黏性流动。
对于Bingham流体的屈服应力阈值,可用经典逾渗理论模型表达[11]:
τy=C(φ−φc)n (2) 式中:φ为流体中的固相粒子质量分数,在B炸药中则对应扩散在溶液体系中的RDX颗粒的质量分数,本文对于RDX颗粒,只考虑颗粒浓度对黏度的影响,忽略颗粒度大小的差异;φc为逾渗阈值;C和n都是通过试验拟合出的常数,此处各参数取值见表1[11]。
表 1 屈服应力阈值计算参数Table 1. Parameters to calculate the yield stress thresholdφc C/Pa n 0.336 6 874 2.37 关于B炸药悬浮液的表观黏度,可以Quemada方程表述[15]:
μapp=μc(1−φφmax (3) 式中:μc为连续组分中无悬浮粒子溶剂的黏度,B炸药中则对应TNT的黏度;φmax为在不同应变率下悬浮液体系中固相粒子的极限质量分数。
根据Zerkle等[10]的落球黏度试验,μc可用以下经验公式表达:
{\mu _{\rm{c}}} = 0.012\;7{\left( {T - {T_{{\rm{m}}.{\rm{TNT}}}}} \right)^{0.304}} (4) 式中:T为当前温度,Tm,TNT为TNT的熔化温度,取值80 ℃。
B炸药悬浮液中最大固相粒子浓度φmax是一个与应变率相关的函数[16],其表达式为:
{\varphi _{\max }} = {\left[ {\varphi _0^{ - 1} - \left( {\varphi _0^{ - 1} - \varphi _\infty ^{ - 1}} \right){{\left( {1 + Z{\gamma ^{ - m}}} \right)}^{ - 1}}} \right]^{ - 1}} (5) 式中:φ0和φ∞分别为液相与类固相悬浮液体系中的极限固相粒子质量分数;Z与m为根据试验拟合出的常数。具体计算取值见表2[10]。
表 2 φmax计算参数Table 2. Calculating parameters of φmaxφ0 φ∞ Z m 0.55 0.713 8.73 0.398 对于相变后B炸药中的RDX,存在3种状态:(1)溶解在液态TNT中;(2)以固态悬浮于溶液体系中;(3)因受热发生分解以气体产物的形式逃逸出溶液体系。RDX的分解速率会随着温度的升高而加快,溶解度也随温度的升高而增大。所以,悬浮于液态B炸药中的固相RDX粒子浓度分数φ是一个与温度相关的函数,可以用以下经验公式表达[16]:
{\varphi = {\varphi _{\rm{a}}} - \left( {A + B{T^*} + C{T^{*2}} + D{T^*}} \right)} (6) {{T^*} = \frac{{T - {T_{{\rm{m}},{\rm{TNT}}}}}}{{{T_{{\rm{m}},{\rm{RDX}}}} - {{{T}}_{{\rm{m}},{\rm{NT}}}}}}} (7) 式中:φa为RDX的初始质量分数;A、B、C、D都是根据测量结果拟合出的常数;Tm,RDX是RDX的理论熔点。各参数具体取值见表3。
表 3 φ计算参数Table 3. Calculating parameters of φφa A B C D Tm,RDX/℃ 0.6 0.042 6 0.075 3 0.134 0.747 4 204 至此,一个基于Bingham流体,包含温度、剪切应变率与固相粒子浓度的完整B炸药黏度计算模型已建立完成。
2. 慢速烤燃数值模拟
使用多物理场分析软件FLUENT对B炸药的慢速烤燃进行数值模拟分析,将上述所建立的黏度模型编写成UDF代码嵌入FLUENT中[17]。
2.1 有限元模型
建立B炸药烤燃弹的网格模型,网格尺寸控制在1 mm左右。弹体尺寸为
\varnothing 35 mm×112 mm,壳体厚度为3.5 mm,材料为45#钢;药柱尺寸为\varnothing 28 mm×105 mm。在弹体的几何中心与上下端距药柱边界17 mm处设置3个温度监测点,如图1所示。2.2 参数设置
壳体材料的具体参数见表4[8],除黏度外,B炸药的其他参数见表5与表6[7, 18]。
表 4 壳体材料参数Table 4. Material parameters for the shell材料 密度/(kg∙m−3) 比热容/(J∙kg−1) 导热系数/(W∙m−1∙K−1) 45#钢 7 850 475 15 表 5 B炸药物性参数Table 5. Physical parameters of Comp B密度/(kg∙m−3) 比热容/(J∙kg−1) 导热系数/(W∙m−1∙K−1) 熔化热/(kJ∙kg−1) 相变起始温度Ts/℃ 相变结束温度T1/℃ \begin{array}{*{20}{c}} {1\;690}&{T {\simfont\text{≤}} {T_s}} \\ {1\;690 - 0.675\left( {T - {T_{\rm{s}}}} \right)}&{T {\simfont\text{>}} {T_{\rm{s}}}} \end{array} 1 126 \begin{array}{*{20}{c}} {0.17}&{T {\simfont\text{<}} {T_{\rm{s}}}} \\ {0.17 - 0.025\left( {T - {T_{\rm{m}}}} \right)}&{{T_{\rm{s}}} {\simfont\text{≤}} T {\simfont\text{≤}} {T_{\rm{l}}}} \\ {0.15}&{T {\simfont\text{>}} {T_{\rm{l}}}} \end{array} 128 80 82 表 6 B炸药化学反应动力学参数Table 6. Chemical kinetic parameters of Comp B活化能/(kJ∙mol−1) 反应热/(J∙kg−1) 指前因子 22 5.82×106 2.01×1018 设置外壁升温速率分别为1 ℃/min与0.055 ℃/min,进行数值模拟计算。其中升温速率为0.055 ℃/min时,先以10 ℃/h的升温速率升至50 ℃并保温1 h,再以0.055 ℃/h的升温速率升温直至响应。
2.3 数值模拟
2.3.1 温度场变化特性分析
各升温速率下弹体内部温度测点曲线分别见图2与图3。从图2与图3可知,对于同一尺寸的烤燃弹,升温速率越快,炸药内部与壳体外壁面的温差越大。1 ℃/min升温速率时,炸药相变完成后中心温度仍低于底部温度,随着升温的进一步进行,才逐渐表现出“顶部温度高,底部温度低”的液相温度场分布特点。两种升温速率下的数值模拟结果都表明,当炸药临近响应时,温度最高点都出现在测点1附近,这与Glascoe等[4]和Nichols等[5]的研究相符。
通过数值模拟能观察到各升温速率下弹体内部温度场的变化情况,分别如图4与图5所示。
从图4可知,当升温较快时,B炸药完全熔化后,在较短时间内炸药内部仍表现出类似固相温度场的同心椭圆分布特点,说明此时炸药虽已完全熔化为液态,但并未发生明显的黏性流动,热传导仍是主要的传热方式,对流传热对温度场的影响非常微弱。但是,随着温度的继续升高,炸药的黏度降低,流动性增强,黏性流动发生后,微观上极薄的两流层之间的剪切应变率增大,对流强度逐渐增大,对流传热对整个系统的温度场的影响也逐渐增强。在黏性流动的影响下,温度场的分布也逐渐转变为典型的液相温度场层状分布的特点。当温度继续升高时,烤燃弹上部开始出现明显的自热反应,最终发生响应。
从图5可知,当升温较慢时,炸药完全熔化后,内部温度场在很长时间内都表现为固相温度场的特征,没有发生明显的黏性流动,直到炸药中心发生自热反应后,内部温度场才逐渐开始发生变化,这也与Maienschein等[3]在试验中观察到的现象相符。升温速率较慢时,B炸药完全熔化后,很长时间内其内部的流场强度都很小,对流传热对炸药内部温度场的影响很弱,只有当发生明显的自热反应后,炸药内部的温度差增大,才逐渐出现明显的黏性流动,温度场分布特征也随之发生变化,最终烤燃弹上部发生响应。值得注意的是,两种升温速率下,响应时刻点火位置都出现在弹体上部。但是,最早出现自热反应的区域却有明显差异。
2.3.2 流场变化特征及分析
两种升温速率下烤燃弹内部的流场平均速度变化曲线见图6。
从图6可知,1 ℃/min的升温速率下,烤燃弹内的流场平均流速曲线的变化趋势可大致分为4个阶段:首先是短暂的稳定期,随后是上升期,然后急速下降,最后急速上升直至响应。当升温速率为0.055 ℃/min时,烤燃弹内部的流场平均流速变化只表现出较长时间的稳定期与急速上升期两个特征阶段。
在1 ℃/min的升温速率下,当外壁温度约110 ℃时,弹体内部的炸药开始屈服并出现流动,流动速度逐渐增大,此时炸药还未发生明显的自热反应,这说明外壁温度快速升高导致的药柱中心位置与靠近壁面位置之间较大的温度差是引起炸药屈服并流动的原因。随着黏性流动的出现,炸药逐渐体现出剪切稀化的特点,流场强度逐渐增强,但随后出现了急速的下跌,又迅速升高直至相应。出现下跌的原因是炸药发生了自热反应,内部流场的流动方向出现了变化,其速度矢量图见图7。
在自热反应发生前,靠近壁面的炸药温度较高,密度较低,在浮升力的驱动下沿着壁面向弹体顶部流动,驱使顶部的炸药从弹体中心向下流动,平衡后炸药内部总体的流动趋势是从靠近壁面处向上流动至顶部,再从顶部中心向底部流动。而当炸药发生自热反应后,顶部的温度高于外壁,顶部靠近外壁处的炸药温度低于顶部中心,在浮升力的驱动下开始沿壁面向下流动,而弹体中心的炸药则向上流动,总体流动趋势与发生自热反应前刚好相反。由于自热反应引起的流动方向的改变造成了流场平均速度的骤跌与骤升,且流速曲线出现骤跌处的时间点可以用来准确判断炸药出现自热反应的时刻,即图6中的A点。
当升温速率为0.055 ℃/min时,在炸药出现自热反应之前,内部温度差很小,流速几乎为零,在很长时间内液态炸药在宏观上都处于几乎静止的状态,所以温度云图表现出类似于固相温度场的特点。当出现明显的自热反应后,内部温度差开始增大,炸药逐渐屈服并开始出现黏性流动,且流动速率呈类指数形式上升,最终响应。炸药自热反应前后烤燃弹内部的流场速度矢量见图8。这说明在升温速率较慢时,炸药自热反应导致的内部温度差升高是使其屈服并流动的主要原因,平均流速曲线出现突变的时刻即是炸药出现自热反应的准确时间点,即图6中的B点。
升温速率较慢时,炸药发生明显自热反应的区域在弹体中心,这也初步解释了前人的研究工作中出现的矛盾观点。无论升温速率快慢,自热反应出现后,炸药均表现出较强的流动性,响应点均在靠近药柱顶部的位置,且都可用流场流速的突变点准确判断炸药出现自热反应的时刻。
3. 试验验证
慢速烤燃试验采用与图1等尺寸的烤燃弹,并在内部设置3个温度测点,且测温点位置与数值模拟中完全一致。使用加热套筒对烤燃弹均匀加热,加热套筒尺寸为
\varnothing 86 mm×194 mm,材料为45钢,厚度2 mm,烤燃弹与加热装置照片见图9。将电加热带均匀缠绕在加热套筒外壁并裹上硅酸铝保温层,用铁丝将烤燃弹体竖直吊放在加热套筒中心,通过MR13温控仪控制加热功率并记录温度数据,采样周期为1 s。响应后收集弹体残骸照片见图10。升温速率为1 ℃/min时,B炸药在184 ℃时发生响应,初步判断响应等级为部分爆轰,全过程中外壁与内部各温度测点的温度曲线见图11。升温速率为0.055 ℃/min时,烤燃弹在181 ℃时发生响应,初步判断响应等级为爆轰,此升温速率下各测点的温度曲线见图12,升温速率越低,炸药响应越剧烈。
两种升温速率下测点1的计算值与试验值曲线比较分别见图13与图14,两种升温速率下测点1温度的计算值与试验值的变化情况基本吻合,且由图11与图12也能看出两种升温速率下烤燃试验中最早出现自热反应的区域和最终点火点的位置分布与数值模拟的结果相同,由此可间接印证前文所建立的黏度模型能很好地表征烤燃过程中B炸药熔化后的内部真实流动情况。
4. 结 论
(1)在慢速烤燃过程中,B炸药熔化后虽变为液态,但仍会表现出一定的“类固性”,只有当浮升力产生的内部应力超过某一极限时,才会屈服并出现黏性流动,影响其内部温度分布特点。升温速率较快时,炸药熔化后内部各处的温度差较大,各部分逐渐屈服并流动,对流传热导致其温度场逐渐转变为典型的液相温度场,炸药最终在上部出现自热反应并响应。而升温速率较慢时,炸药熔化后内部的温度差较小,在很长时间内炸药均保持宏观静止,温度场也体现出典型固相温度场的分布特点,当中心发生自热反应后,内部温差增大,炸药才逐渐开始流动并在上部发生响应。
(2)不同升温速率下B炸药熔化后内部的流场流动特点也有很大的差异。升温速率较快时,炸药在熔化后出现自然对流,流场平均流速逐渐增大,自热反应发生后流场平均流速先骤跌再急剧上升直至响应,骤跌的出现是由于自热反应发生前后炸药内部流场的流动方向发生了变化。升温速率较慢时,炸药在熔化后内部流场的平均流速仍接近于0,直到炸药发生自热反应后才开始急剧上升直至响应。内部流场平均流速出现突变的时间点可用于准确判断炸药出现自热反应的时刻。
(3)通过试验结合数值模拟验证了基于Bingham流体的B炸药黏度模型能很好地描述B炸药在慢速烤燃过程中的流动情况,可用于预测与分析B炸药烤燃过程中自热反应区域与点火位置的分布规律,此模型对于其他种类熔铸炸药或含铝炸药是否适用值得进一步研究。
-
表 1 屈服应力阈值计算参数
Table 1. Parameters to calculate the yield stress threshold
φc C/Pa n 0.336 6 874 2.37 表 2 φmax计算参数
Table 2. Calculating parameters of φmax
φ0 φ∞ Z m 0.55 0.713 8.73 0.398 表 3 φ计算参数
Table 3. Calculating parameters of φ
φa A B C D Tm,RDX/℃ 0.6 0.042 6 0.075 3 0.134 0.747 4 204 表 4 壳体材料参数
Table 4. Material parameters for the shell
材料 密度/(kg∙m−3) 比热容/(J∙kg−1) 导热系数/(W∙m−1∙K−1) 45#钢 7 850 475 15 表 5 B炸药物性参数
Table 5. Physical parameters of Comp B
密度/(kg∙m−3) 比热容/(J∙kg−1) 导热系数/(W∙m−1∙K−1) 熔化热/(kJ∙kg−1) 相变起始温度Ts/℃ 相变结束温度T1/℃ \begin{array}{*{20}{c}} {1\;690}&{T {\simfont\text{≤}} {T_s}} \\ {1\;690 - 0.675\left( {T - {T_{\rm{s}}}} \right)}&{T {\simfont\text{>}} {T_{\rm{s}}}} \end{array} 1 126 \begin{array}{*{20}{c}} {0.17}&{T {\simfont\text{<}} {T_{\rm{s}}}} \\ {0.17 - 0.025\left( {T - {T_{\rm{m}}}} \right)}&{{T_{\rm{s}}} {\simfont\text{≤}} T {\simfont\text{≤}} {T_{\rm{l}}}} \\ {0.15}&{T {\simfont\text{>}} {T_{\rm{l}}}} \end{array} 128 80 82 表 6 B炸药化学反应动力学参数
Table 6. Chemical kinetic parameters of Comp B
活化能/(kJ∙mol−1) 反应热/(J∙kg−1) 指前因子 22 5.82×106 2.01×1018 -
[1] 智小琦. 弹箭炸药装药技术 [M]. 北京: 兵器工业出版社, 2012: 7−8.ZHI X Q. Charge technology of explosive on projectile and rocket [M]. Beijing: Ordnance Industry Press, 2012: 7−8. [2] HOBBS M L, KANESHIGE M J, ERIKSON W W. Predicting large-scale effects during cookoff of PBXs and melt-castable explosives [C] // Proceedings of the 26th International Colloquium on the Dynamics of Explosions and Reactive Systems. Boston: MA: OSTI, 2017: 152−158. [3] MAIENSCHEIN J L, McCLELLAND M A, WARDELL J F, et al. ALE3D model predictions and experimental analysis of the cookoff response of Comp B [C] // Proceedings of Joint Army Navy NASA Air Force (JANNAF) Meeting. Colorado Springs, CO, USA: LLNL, 2003: 51−64. [4] GLASCOE E A, DEHAVEN M R, MCCLELLAND M, et al. Mechanisms of Comp-B thermal explosions [C] // Proceedings of the 15th International Detonation Symposium. San Francisco: LLNL, 2014: 376−385. [5] NICHOLS A L, SCHOFIELD S. Modeling the response of fluid/melt explosives to slow cook-off [C] // Proceedings of the 15th International Detonation Symposium. San Francisco: LLNL, 2014: 1128−1136. [6] FEDOROFF B T, SHEFFIELD O E, KAYE S M. Encyclopedia of explosives and related items [M]. Dover, NJ: Picatinny Arsenal, 1962. [7] HOBBS M L, KANESHIGE M J, ANDERSON M U. Cookoff of a melt-castable explosive (COMP-B) [C] // Proceedings of the 27th Propulsion Systems Hazards Joint Subcommittee Meeting. Monterrey: SNL-NM, 2012: 1020−1034. [8] SANHYE W, DUBOIS C, LAROCHE I, et al. Numerical modeling of the cooling cycle and associated thermal stresses in a melt explosive charge [J]. AIChE Journal, 2016, 62(10): 3797–3811. DOI: 10.1002/aic.15288. [9] NUNEZ M P, ZERKLE D K, ZUCKER J M. The rheology of molten Composition B [R]. NM: Los Alamos, 2012. [10] ZERKLE D K, NUNEZ M P, ZUCKER J M. Molten composition B viscosity at elevated temperature [J]. Journal of Energetic Materials, 2016, 34(4): 368–383. DOI: 10.1080/07370652.2015.1102179. [11] DAVIS S M, ZERKLE D K, SMILOWITZ L B, et al. Integrated rheology model: explosive composition B-3 [J]. Journal of Energetic Materials, 2018, 36(4): 398–411. DOI: 10.1080/07370652.2018.1451573. [12] DAVIS S M, ZERKLE D K. Short communication: estimation of yield stress/viscosity of molten Octol [J]. AIP Advances, 2018, 8(5): 055202. [13] MACOSKO C W. Rheology principles, measurements and applications [M]. New York: VCH Publishers, 1994: 92−98. [14] MORRISON F A. Understanding rheology [M]. New York: Oxford University Press, 2001: 232. [15] QUEMADA D. Rheology of concentrated disperse systems and minimum energy dissipation principle: I. viscosity-concentration relationship [J]. Rheologica Acta, 1977, 16(1): 82–94. DOI: 10.1007/BF01516932. [16] ZHOU J Z Q, UHLHERR P H, LUO F T. Yield stress and maximum packing fraction of concentrated suspensions [J]. Rheologica Acta, 1995, 34(6): 544–561. DOI: 10.1007/BF00712315. [17] Fluent Inc. FLUENT user’s guide [M]. US: Fluent Inc, 2006. [18] MCCLELLAND M A, GLASCOE E A, NICHOLS A L, et al. ALE3D simulation of incompressible flow, heat transfer, and chemical decomposition of Comp B in slow cookoff experiments [C] // Proceedings of International Detonation Symposium. San Francisco: LLNL, 2014: 517−528. 期刊类型引用(5)
1. 徐跃跃,应武江,关通,王泽钰,高家乐,苗飞超,张向荣,周霖. DNP基熔铸炸药烤燃试验与数值模拟研究. 北京理工大学学报. 2025(02): 126-136 . 百度学术
2. 王琦,智小琦,肖游,郝春杰. 基于UCM模型的B炸药慢烤泄压结构的作用分析. 爆炸与冲击. 2022(04): 36-44 . 本站查看
3. 张海军,聂建新,王领,王栋,胡峰,郭学永. 主控提前点火对复合推进剂慢速烤燃响应的影响. 爆炸与冲击. 2022(10): 3-13 . 本站查看
4. 徐瑞,智小琦,王帅. 缓释结构对B炸药烤燃响应烈度的影响. 高压物理学报. 2021(03): 149-157 . 百度学术
5. 周捷,智小琦,王帅,范兴华. 升温速率与流变特性对B炸药慢烤响应的影响. 爆炸与冲击. 2020(12): 70-78 . 本站查看
其他类型引用(3)
-