JOB-9003和JB-9014炸药平面爆轰驱动飞片的对比研究

采用双灵敏度激光速度干涉仪(VISAR)测量了JOB-9003和JB-9014两种炸药平面一维爆轰驱动飞片的速度历史,比较了二者的做功能力;并用DYNA2D程序进行了数值模拟。在有JO-9159炸药的条件下,实验和计算结果表明,JB-9014炸药的做功能力可以达到JOB-9003炸药的75.0%;无JO-9159炸药条件下的计算结果表明,JB-9014炸药的做功能力约为JOB-9003炸药的71.2%     

JB-9014炸药超压爆轰产物的状态方程

根据P.K.Tang等提出的对炸药爆轰产物超压状态方程建模时只对JWL状态方程CJ等熵线中 高压指数项做修正的研究思路,首先给定超压状态下内能等熵线的修正项,再根据热力学定律对内能等熵线 求微分而得到沿压力等熵线的修正项。对JB-9014炸药超压爆轰冲击Hugoniot实验数据和声速实验数据同 时进行拟合,得到了3个JWL状态方程在超压爆轰状态下的修正项,并进行了分析与比较。得到超压修正项 的方法简单,3组超压修正项与P.K.Tang的修正项一样,都能很好地拟合超压Hugoniot数据。在实验数据 范围外,对超压状态下的声速-压力实验数据的拟合精度有所差别。     

多级炸药爆轰高速驱动技术的Gurney模型优化分析

利用多级炸药 飞片装置可得到速度在10km/s以上、直径数厘米的末级平面飞片,在高压物态方程和高速碰撞等实验研究中有重要应用前景。改进了爆炸驱动的Gurney模型,推广用于典型的二级爆轰驱动(前级飞片 炸药层 后级飞片)装置的分析,得到前、后级飞片终速度的计算公式,与实验结果符合较好。根据对炸药层厚度优化计算结果的讨论可知,对于综合利用前级飞片动能和炸药层能量最有利的后级飞片速度处于2～3倍Gurney速度的范围。前级飞片速度很高时,后级飞片速度趋于两者直接弹性碰撞的结果。     

一类爆轰合成用乳化炸药的爆轰反应特征

结合热分解特性、爆轰产物状态及数值计算探讨所制备的含Fe、Mn(Zn)元素乳化炸药的爆轰反 应特征。对该系列乳化炸药及其爆轰产物分别进行DSC和DTG实验、XRD检测和TEM 表征,通过比较不 同类型炸药及其爆轰产物的热分解特征、爆轰产物的成分和形貌、晶型畸变度等研究该类炸药爆轰反应特征; 并通过数值计算研究和佐证了不同温度和压强状态下,爆轰固体产物可能存在状态。结果表明:硝酸铵有利 于炸药爆轰时形成相对比较均匀和平稳的爆轰反应结构,而这种结构有利于可重复性合成比较单一和均匀 的爆轰产物。高密度炸药的爆轰产物比低密度炸药的相对纯净。低密度炸药不完全爆轰,导致了爆轰产物杂 质较多;数值计算结果表明,爆轰固体产物分布相图可以辅助分析爆轰固体产物可能存在状态及Mn(Zn)的 爆轰合成机制。     

钝感炸药的超压爆轰与冲击起爆过程数值模拟

采用Hybrid反应率结合修正的JWL方程,研究了LX-17、超细TATB等钝感炸药的冲击起爆（SDT）过程,并计算了爆轰波的对碰现象。结果表明,该方法计算钝感炸药的冲击起爆过程与实验数据符合较好;计算爆轰波对碰区的峰值压力提高了10％。     

高密度B炸药的燃烧转爆轰实验研究

采用电探针及压力传感器测试技术对密度为1.597 g/cm3的固体B炸药（TNT/RDX=40/60）的燃烧转爆轰性能进行研究。实验结果表明,在较强的约束条件下（45号钢管,内径20 mm,外径64 mm,长500 mm）,B炸药形成了DDT现象,诱导爆轰距离为295~310 mm。     

PBX9502炸药爆轰产物的状态方程

采用vanderWaals等效单组分流体模型和Ross硬球微扰理论软球修正模型,计算了爆轰气相产 物的状态方程;用石墨相、金刚石相、类石墨液相和类金刚石液相4种相态描述凝聚成分,由Gibbs自由能最 小确定了不同状态下的凝聚产物相态。对爆轰产物混合系统采用自由能最小原理,通过化学平衡方程组求解 了炸药爆轰产物系统的平衡组分。使用该理论计算了高含碳炸药PBX9502Chapman-Jouguet(CJ)点的爆轰 参数,计算值与实验值符合很好;同时计算了3条等温线,并与Sesame库比较,发现温度超过1000 ℃时,计 算值与Sesame库的计算结果比较接近。在计算的5802K 等温线上发现了一个拐点,分析发现是由于在此 处游离态的碳发生了相变。     

含铝炸药水下滑移爆轰实验研究

为了研究含铝炸药水中爆炸近场范围内沿装药径向的冲击波特性及爆轰产物的膨胀规律,设计了水下滑移爆轰实验装置。采用高速扫描相机和阴影照相技术记录了一种RDX基含铝炸药的水下滑移爆轰过程,获得了清晰的水下滑移爆轰过程的扫描图像。通过对图像进行分析,得到了水中爆炸近场范围内沿装药径向冲击波的传播迹线及爆轰产物气泡的膨胀迹线,进而分析出冲击波传播速度、阵面压力及爆轰产物气泡的膨胀位移等参数的变化规律。此外,根据气泡的膨胀迹线,标定了该含铝炸药爆轰产物的JWL状态方程参数,并将其与Φ50mm圆筒实验确定的参数值在p-V图上进行了对比。结果显示,两者偏差较小,证实了该方法的可行性。     

利用液体炸药爆轰制备石墨烯薄片

强氧化性酸的环境下,石墨易形成石墨层间化合物(GICs),利用石墨这一特性,将天然石墨置于发烟硝酸中,并加入硝基甲烷,配制成液体炸药,使用塑料容器盛装后放入爆轰反应釜中引爆。收集爆轰产物,并利用XRD、EDX、SEM、TEM、Raman光谱、比表面积与孔隙度分析仪进行分析,结果显示:制备出的石墨烯具有良好的晶体特性并呈现极薄的片状结构,其比表面积达到天然石墨的9.16倍,平均厚度约为14.73 nm。     

炸药爆轰产物驱动不可压缩刚粘塑性柱壳的运动

从质量和动量守恒关系出发 ,结合炸药爆轰产物压力剖面的分析和假设 ,导出了不可压缩刚粘塑性柱壳在柱面散心和聚心爆轰产物推动作用下的运动方程。应用该方程式计算爆轰驱动柱壳的运动与实验结果及动力有限元程序DYNA2D的数值模拟结果十分接近 ,说明散心爆轰波场合波后产物流场可用柱面Taylor波的近似解合理描述 ,有限初始半径聚心爆轰波后产物流场可用平面Taylor波近似描述。     

一期水管冷却效应的数值模拟新方法

采用有限单元离散空间域，视冷却水管所在单元为子结构，其它单元为常规单元，然后，把混凝土结构瞬态温度场的有限元控制方程表示成状态方程，最后，利用精细积分法求解．结果表明，子结构、状态方程和精细积分技术的引用降低了对计算机贮存量的要求和提高了冷却水管附近温度梯度的精度，显示出极大的优越性．     

爆轰波拐角传播三维数值模拟

为了研究爆轰波在拐角中的传播特性以及拐角爆轰低压流场现象,用LSDYNA 3D程序对三种常见的拐角装药的爆轰波传播特性进行了数值模拟,清晰显示了爆轰波因拐角处爆轰波传播面积的变化而产生的衰减 增长过程,讨论了侧向三通中爆速的变化情况,并得到了拐角处低压爆轰区的尺寸。     

Plastic multipliers as driving variables of numerical simulation in elastoplasticity

Statics of elasto-plastic media is stated in terms of eigenstrains acting upon a background elastic problem with fixed (initial) stiffness. In order to minimize the number of unknowns and to provide computationally cheap algorithms plastic multipliers are used as main driving variables. Fixed-point type iterations are suggested for computing plastic multipliers within the load increment. Numerical experiments results are shown to be in good correlation with other software results. A range of problems is outlined where the usage of the proposed algorithms can be advantageous.     

弯管内爆轰波传播的流场显示和数值模拟

采用激光纹影系统拍摄了爆轰波在不同位置的流场照片. 用二阶附加半隐的龙格- 库塔法和五阶WENO格式 分别离散欧拉方程时间和空间导数项,用基元反应来描述爆轰化学反应过程,获得了压力、 温度、典型组元质量分数分布及数值胞格结构和爆轰波平均速度. 结果表明：受壁面稀疏波 和压缩波影响,爆轰波阵面发生畸变. 但由于弯管曲率半径较大,未出现爆轰波熄灭. 靠近 凹壁面的激波强度大于凸壁面侧,且凹壁面侧的反应区宽度较凸壁面侧要窄. 弯管出口处的 三波点数目较入口处减少,爆轰波衰减. 在出口直段,受扰动的爆轰波可恢复为自持爆轰波. 爆轰波流场、胞格结构、平均爆轰波速度的计算和实验结果定性一致.     

基于误差马尾图量化爆轰数值模拟结果的置信度

针对爆轰流体力学数值模拟过程中输入参数的不确定性, 通过抽样技术, 形成确定性爆轰流体力学程序的各种输入和数值求解, 建立输入参数与输出响应量的样本, 再通过概率框架下的误差累积分布函数与马尾图, 给出了爆轰数值模拟过程中输入参数不确定度对模拟结果影响的置信度量化方法。通过一维黎曼问题、平面爆轰问题计算了误差马尾图, 给出了二维爆轰拉氏自适应流体动力学LAD2D程序计算网格与模拟结果置信度的关系, 对多物理爆轰过程发展高置信度数值模拟软件有很好的借鉴作用。

     

A numerical simulation of reflection processes of a detonation wave on a wedge

Abstract. A two dimensional numerical simulation has been performed to study reflection processes of detonation waves on a wedge. The numerical scheme adopted is the flux corrected transport scheme and a two-step chemical reaction is assumed for a stoichiometric oxyhydrogen mixture diluted with argon. Transverse wave structures of the detonation are produced by artificial disturbances situated in front of a one-dimensional Chapman-Jouguet detonation wave. Numerical grids are generated by solving a Laplace equation. Results show that in the case where Mach reflection occurs, the cells in the Mach stem are smaller than those in the incident wave and are distorted in shape. There is also an initiating stage during which the cells in the Mach stem are created. The critical angle beyond which Mach reflection cannot occur is discussed. Received 15 October 1999 / Accepted 27 March 2000     

动态存储方法在气相爆轰波数值模拟中的应用

为减少反应流计算花费的时间,采用动态存储建立热力学数据表（In situ adaptive tabulation, ISAT）的方法来取代反应流计算中化学反应的直接积分（Direct integral, DI）,并针对气相H2/O2爆轰波的传播过程进行了数值模拟。采用两种不同的ISAT误差放大判据对一维爆轰波的压力、温度及组分变化进行了计算,并和DI的结果进行了比较以考察其计算精度。此外,还探讨了ISAT方法的效率。计算结果表明,和DI方法相比,在爆轰波传播计算中ISAT方法计算误差小于3％,化学反应计算速度提高了8倍。     

Meshless numerical simulation for fully nonlinear water waves

A meshless numerical model for nonlinear free surface water wave is presented in this paper. An approach of handling the moving free surface boundary is proposed. Using the fundamental solution of the Laplace equation as the radial basis functions and locating the source points outside the computational domain, the problem is solved by collocation of only a few boundary points. Present model is first applied to simulate the generation of periodic finite‐amplitude waves with high wave‐steepness and then is employed to simulate the modulation of monochromatic waves passing over a submerged obstacle. Good agreements are observed as compared with experimental data and other numerical models. Copyright © 2005 John Wiley & Sons, Ltd.     

圆柱形爆轰波的二维数值模拟

基于带化学反应的二维Euler方程,对圆柱形爆轰波的直接起爆和传播过程进行了二维数值模拟研究,拟分析起爆条件和初始压强对圆柱形爆轰波形成和传播的影响。研究发现,圆柱形爆轰波起爆成功向外传播的过程中,新的三波结构的生成标志着爆轰波进入稳定传播阶段。在起爆能量足够的情况下,起爆半径(曲率)的大小决定着三波结构初始形成时的数目和传播半径,起爆压强对其基本不产生影响;起爆半径大(曲率小)时,三波结构初始形成时的传播半径大、数目多,圆柱形爆轰波进入稳定传播阶段的传播距离长;数值模拟中,初始压强的提高,有助于圆柱形爆轰在较短的传播距离内进入稳定传播阶段。