炸药爆轰产物Jones-Wilkins-Lee状态方程不确定参数

Jones-Wilkins-Lee (JWL)状态方程是一种不显含化学反应、由实验方法确定参数的半经验状态方程, 能比较精确地描述爆轰产物的膨胀驱动做功过程. 在JWL状态方程中有多个未知(不确定)参数需要确定. 传统的确定JWL状态方程参数的方法是“调参数”, 人为因素影响较大, 无法给出参数的不确定性信息. 本文利用贝叶斯分析方法研究了炸药的不确定参数, 该方法能够基于以往的认识、实验和模拟数据标定(calibration)不确定参数. 在本文结果中, 不确定参数的后验分布均值与文献结果相符合, 基于参数标定结果的数值模拟90%置信区间完全包含实验数据. 数值标定结果说明贝叶斯参数标定适用于确定样品炸药的JWL状态方程参数. 特别是, 在本文JWL状态方程参数标定过程中极大地减少了人为因素的影响.     

PETN炸药爆轰产物状态方程的理论研究

 采用van der Waals等效单组分流体模型和Ross硬球微扰理论软球修正模型,计算爆轰气相产物的状态方程;用石墨相、金刚石相、类石墨液相和类金刚石液相4种相态描述凝聚成分,由Gibbs自由能最小原理确定不同状态下的凝聚产物相态。对爆轰产物混合系统采用自由能最小原理,通过化学平衡方程组求解炸药爆轰产物系统的平衡组分。使用该理论计算PETN炸药Chapman-Jouguet（CJ）点的爆轰参数,其值与实验值符合得很好;同时计算了以CJ点为起始点的等熵卸载线,并与传统的Jones-Wilkins-Lee（JWL）状态方程的计算结果进行比较,发现计算的γ值是单调递减的,而JWL状态方程计算的γ值却出现了“双峰”现象。分析认为,传统的JWL状态方程给出的“双峰”变化,是由其函数形式自身决定的,并不对应实际物理过程。     

爆轰产物JWL状态方程应用研究

 通过对凝聚态炸药爆轰产物JWL状态方程的研究,提出了一种确定JWL状态方程参数的γ拟合法,即利用γ律状态方程来拟合JWL状态方程参数。这种新方法不需进行圆筒试验,与现有方法比较,具有经济、安全、方便、准确的特点。利用这种方法拟合确定了4种常用炸药TNT、C-4、PETN、HMX的JWL状态方程参数,通过与圆筒试验法得出的JWL状态方程p-V曲线进行对比,证明γ拟合法具有较高的精度,完全可以满足计算爆炸力学的应用。     

钝感炸药圆筒试验与爆轰产物JWL状态方程研究

 对以TATB为主的钝感炸药JB-9014进行了25 mm和50 mm两种装药直径的圆筒试验，测试了TU1圆筒在爆轰产物驱动下的膨胀过程R(t)关系。对圆筒试验进行了综合分析和二维数值模拟计算，综合评估了JB-9014的作功能力和圆筒试验的相似性。通过二维流体动力学数值模拟，确定得到了JB-9014炸药爆轰产物JWL状态方程参数。经对JB-9014二维平面滑移爆轰驱动试验的数值模拟检验，证明确定的JWL状态方程参数是可靠的，具有较高的精度和普适性。     

多元混合炸药爆轰驱动圆筒膨胀规律的理论确定方法

为快速预估任意配比的多元混合炸药爆轰产物的JWL(Jones-Wilkins-Lee)参数,提出了快速确定多元混合炸药爆轰驱动圆筒膨胀规律的理论方法,即在给定各组分爆轰产物JWL参数的前提下,根据能量守恒定律,采用Gurney模型,确定圆筒试验中多元混合炸药爆轰驱动圆筒膨胀距离随时间变化的曲线。同时,利用能量守恒原理以及经典爆轰理论中通过常γ状态方程得到的爆速、爆压和爆热之间的关系式,提出了确定多元混合炸药爆速和爆压的方法。采用该理论方法,分别计算了多元混合炸药PBXC03和PBXC10爆轰驱动圆筒膨胀规律及爆速和爆压,计算结果与前人的实验结果符合较好,验证了该理论方法的可行性和有效性。     

论VLW状态方程

 论述了VLW状态方程中爆轰产物分子势参数的确定；根据不同类型炸药爆轰性能参数计算结果对炸药爆轰产物碳在CJ态的相进行了分析讨论，认为不同的炸药，产物碳有石墨和金刚石两种相存在；将VLW状态方程和BKW状态方程进行了比较分析，指出了BKW状态方程与VLW状态方程的区别及其产生的原因，并给出了VLW状态方程计算含能材料（包括高能炸药、燃料空气炸药、民用炸药）的爆轰性能参数计算结果。结果表明，计算值与实验值符合较好。     

战斗部爆轰产物温度效应数值模拟

 利用BKW代码、JWL方程和LS-DYNA程序，研究了无约束、弱约束和强约束等三种战斗部约束条件下炸药爆轰产物密度随时间的变化关系，再结合爆轰产物HOM状态方程获得爆轰产物温度随时间的衰变函数。给出了战斗部爆轰产物温度的理论计算方法。研究表明，爆轰区内爆轰产物温度及其衰减速率并不均匀，战斗部炸药的约束条件对爆轰产物密度和温度有重要影响，强约束可以延缓爆轰产物密度和温度的衰减，相同情况下仅战斗部壳体厚度加倍并不会导致爆轰产物密度和温度的衰减情况同比例变化。     

水下爆炸气泡脉动周期的简便计算方法

将爆轰产物状态方程与气泡运动方程相结合,提出了一种水下爆炸气泡脉动周期的数值计算方法,能够快速计算不同炸药水下爆炸的气泡脉动周期。结果表明:采用JWL状态方程与采用γ律状态方程的气泡周期计算方法相比,前者计算TNT炸药水下爆炸气泡脉动周期的误差更小,计算误差小于3%;通过对比RS211炸药水下爆炸实验结果,进一步验证了该数值计算方法对含铝炸药水下爆炸气泡脉动周期的计算同样具有良好的适用性。     

炸药爆轰的连续介质本构模型和数值计算方法

假定炸药和爆轰产物处于局部热力学平衡状态, 即它们的压力和温度相同, 利用热力学基本关系建立炸药爆轰过程的连续介质本构模型的一般理论框架. 在此框架下, 炸药爆轰本构模型由一组常微分方程构成, 包括炸药和爆轰产物的状态方程、简单混合法则、化学反应速率方程和能量守恒方程, 易于由成熟的计算方法如梯形法等进行求解. 一组广义Maxwell型非线性固体本构形式的微分方程描述了压力和温度随时间的演化速率与应变率和化学反应速率的关系, 借助简单混合物理论, 其中的系数由炸药和爆轰产物的材料参数确定. 未反应的炸药和爆轰产物采用JWL状态方程, 化学反应率方程采用Lee-Tarver点火-燃烧二项式模型, 模拟PBX-9404炸药的一维冲击波起爆过程和爆轰波传播过程. 计算结果表明了本文给出的本构模型和相应计算方法的有效性. 关键词： 炸药爆轰 本构模型 化学反应率方程 数值模拟     

散心爆轰的数值模拟

本文在瞬时稳态爆轰的假定下,简要叙述了散心爆轰的不定常性质。在散心爆轰的数值模拟方法中,我们改进了人造粘性并用Cochran反应率和JWL状态方程来描述炸药的引爆过程。计了几个典型的散心爆轰问题,得到了比较满意的结果。     

镉镍蓄电池及其在大功率脉冲恒流源中的应用

 介绍了碱性镉镍蓄电池用作脉冲大电流恒流源的可行性，叙述了该类电池放电电流的规律和消除记忆效应的办法，并列出了用镉镍蓄电池组装有关恒流源样机的测试数据和实际考验结果。实验证明：在脉冲式大电流场合下，使用蓄电池作能源，可以消除常用的高功率脉冲源对市电电网的冲击，节省大量设备投资，获得具有高稳定度性能的恒流源。     

基于统计物理的爆轰产物物态方程研究

用van der Waals等效单组分流体模型和Ross硬球微扰理论软球修正模型，计算爆轰气相产物的状态方程；用石墨相、金刚石相、类石墨液相和类金刚石液相４种相态描述凝聚成分，由Gibbs自由能最小确定不同状态下的凝聚产物相态.对爆轰产物混合系统采用Gibbs自由能最小原理，通过化学平衡方程组求解炸药爆轰产物系统的平衡组分，计算结果与Becker-Kistiakowsky-Wilson (BKW)和Lennard-Jones-Devonshire结果相近.使用该理论对炸药的爆轰参数做了预言，与BKW，Jo 关键词： 爆轰产物 物态方程 化学平衡方程组     

CTVD格式数值计算非均质炸药爆轰问题

将高分辨率激波捕捉格式CTVD格式拓展应用到非均质炸药爆轰的数值模拟问题.增加了化学反应率控制方程,引入Lee-Tarver点火成长模型,未反应的固体炸药和化学反应气体产物都使用JWL形式状态方程.数值模拟了非均质固体炸药PBX-9404和TATB的冲击起爆问题.获得了较高的爆轰波分辨率和光滑解区的数值精度,对具有复杂物态方程形式的固体炸药爆轰问题,CTVD格式具有简单实用、高效和高分辨的特点.     

炸药条加载圆柱壳的数值模拟（Ⅰ）：流固耦合模拟

 对炸药条加载圆柱壳模拟X射线力学效应的试验进行了流固耦合数值模拟。首先建立了包括圆柱壳、两层缓冲橡胶、炸药条和空气在内的二维数值模型,其中,两层橡胶采用Ogden超弹性本构模型描述,炸药条及其爆轰产物采用高能炸药燃烧（High Explosive Burn）模型和JWL状态方程描述,空气采用多项式状态方程描述。然后采用多物质ALE流固耦合算法进行计算,获得了爆炸加载的物理图像、载荷传递与结构响应。结果表明,采用19条炸药条加载265 mm的圆柱壳模拟X射线余弦载荷作用,特征点应力响应基本等效;试验附加的两层橡胶对结构响应特征具有影响,在壳体质量较小和（或）刚度较低的情况下,紧密粘接会引起结构响应失真。     

Diffraction of sound by an elastic cylinder near the surface of an elastic halfspace

Diffraction of an acoustic wave by an elastic cylinder near the surface of an elastic halfspace is considered. The solution relies on a Helmholtz-type integral equation and uses the Green function of an elastic halfspace. The latter function is represented in the form of an integral over the Sommerfeld contour on the plane of a complex variable that has the meaning of the angle of the wave incidence on the halfspace boundary. An integral equation for the sound pressure distribution over the cylinder surface is derived. This equation is reduced to an infinite system of equations for the Fourier-series expansion coefficients of this distribution. The results obtained are valid for the diffraction of a cylindrical wave and a plane wave. They also describe the diffraction of a spherical wave when the transmitter and receiver are far from the cylinder and lie in one plane that is orthogonal to the cylinder axis.