搜索

x

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

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

周洪强 于明 孙海权 董贺飞 张凤国

引用本文:
Citation:

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

周洪强, 于明, 孙海权, 董贺飞, 张凤国

A continuum constitutive model and computational method of explosive detonation

Zhou Hong-Qiang, Yu Ming, Sun Hai-Quan, Dong He-Fei, Zhang Feng-Guo
PDF
导出引用
  • 假定炸药和爆轰产物处于局部热力学平衡状态, 即它们的压力和温度相同, 利用热力学基本关系建立炸药爆轰过程的连续介质本构模型的一般理论框架. 在此框架下, 炸药爆轰本构模型由一组常微分方程构成, 包括炸药和爆轰产物的状态方程、简单混合法则、化学反应速率方程和能量守恒方程, 易于由成熟的计算方法如梯形法等进行求解. 一组广义Maxwell型非线性固体本构形式的微分方程描述了压力和温度随时间的演化速率与应变率和化学反应速率的关系, 借助简单混合物理论, 其中的系数由炸药和爆轰产物的材料参数确定. 未反应的炸药和爆轰产物采用JWL状态方程, 化学反应率方程采用Lee-Tarver点火-燃烧二项式模型, 模拟PBX-9404炸药的一维冲击波起爆过程和爆轰波传播过程. 计算结果表明了本文给出的本构模型和相应计算方法的有效性.
    Some basic thermodynamic relationships are used to develop a theroretical framework for modeling the detonation in explosives, on the assumption that explosive and detonation product are in a local thermodynamic equilibrium state, i.e., their pressures and temperatures are identical. Using this framework, a continuum constitutive model for explosive detonation is composed of a group of ordinary differential equations including the state equations of explosive and its product, simple mixing law, reaction rate equation and energy conservation equation, which are easily solved by a mature computational method such as trapezoidal rule. A group of nonlinear constitutive equations in a generalized Maxwellian form describe the relationship among the time evolution rates of pressure and temperature, the strain rate, and the chemical reaction rate. Coefficients appearing in the constitutive equations are determined only by parameters of the explosive and the product through using simple mixing rule. The continuum constitutive model and the corresponding computational method are verified by simulating the detonation behaviour of PBX9404 impacted by high velocity Cu flyer, and in the simulation the JWL equation of state for unreacted explosive and detonation product and the two-term Lee-Tarver reaction rate equation are adopted.
    • 基金项目: 中国工程物理研究院科学技术发展基金(批准号:2013A0201010,2013B0101014)和国家自然科学基金(批准号:11272064)资助的课题.
    • Funds: Project supported by the Science and Technology Development Foundation of China Academy of Enginerring Physics (Grant Nos. 2013A0201010, 2013B0101014) and the National Natural Science Foundation of China (Grant No. 11272064).
    [1]

    Sun J S, Zhu J S 1995 Theory of Detonation Physics (Beijing: National Denfense Industry Press) (in Chinese) [孙锦山, 朱建士 1995 理论爆轰物理 (北京: 国工业出版社)]

    [2]

    Sun C W, Wei Y Z, Zhou Z K 2000 Application of Detonation (Beijing: National Denfense Industry Press) (in Chinese) [孙承纬, 卫玉章, 周之奎 2000 应用爆轰物理 (北京: 国工业出版社)]

    [3]

    Mader C L 2008 Numerical Modeling of Explosives and Propellants (New York: CRC Press) p373

    [4]

    Johnson J N, Tang P K, Forest C A 1985 J. Appl. Phys. 57 4323

    [5]

    Sun C W 1986 Chin. J. Comput. Phys. 3 142 (in Chinese) [孙承纬 1986 计算物理 3 142]

    [6]

    Zhang G R, Chen D N 1991 Detonation Dynamics of Agglomerate Detonation (Beijing: National Denfense Industry Press) p142 (in Chinese) [章冠人, 陈大年 1991 凝聚炸药起爆动力学 (北京: 国防工业出版社) 第142页]

    [7]

    Brenan K E, Campbell S L, Petzold L R 1996 Numerical Solution of Initial-Value Problems in Differential-Algebraic Equations (Philadelphia: SIAM) p4

    [8]

    Ascher Uri M, Petzold L R 1998 Computer Methods for Ordinary Differential Equations and Differential-Algebraic Equations (Philadelphia: SIAM)

    [9]

    Hayes D B 1975 J. Appl. Phys. 46 3438

    [10]

    Johnson J N 1981 J. Appl. Phys. 52 2812

    [11]

    Sun H Q, Zhang W H 2006 Chin. J. Energ. Mater. 14 16 (in Chinese) [孙海权, 张文宏2006 含能材料 14 16]

    [12]

    Lee E L, Tarver C M 1980 Phys. Fluids 23 2362

  • [1]

    Sun J S, Zhu J S 1995 Theory of Detonation Physics (Beijing: National Denfense Industry Press) (in Chinese) [孙锦山, 朱建士 1995 理论爆轰物理 (北京: 国工业出版社)]

    [2]

    Sun C W, Wei Y Z, Zhou Z K 2000 Application of Detonation (Beijing: National Denfense Industry Press) (in Chinese) [孙承纬, 卫玉章, 周之奎 2000 应用爆轰物理 (北京: 国工业出版社)]

    [3]

    Mader C L 2008 Numerical Modeling of Explosives and Propellants (New York: CRC Press) p373

    [4]

    Johnson J N, Tang P K, Forest C A 1985 J. Appl. Phys. 57 4323

    [5]

    Sun C W 1986 Chin. J. Comput. Phys. 3 142 (in Chinese) [孙承纬 1986 计算物理 3 142]

    [6]

    Zhang G R, Chen D N 1991 Detonation Dynamics of Agglomerate Detonation (Beijing: National Denfense Industry Press) p142 (in Chinese) [章冠人, 陈大年 1991 凝聚炸药起爆动力学 (北京: 国防工业出版社) 第142页]

    [7]

    Brenan K E, Campbell S L, Petzold L R 1996 Numerical Solution of Initial-Value Problems in Differential-Algebraic Equations (Philadelphia: SIAM) p4

    [8]

    Ascher Uri M, Petzold L R 1998 Computer Methods for Ordinary Differential Equations and Differential-Algebraic Equations (Philadelphia: SIAM)

    [9]

    Hayes D B 1975 J. Appl. Phys. 46 3438

    [10]

    Johnson J N 1981 J. Appl. Phys. 52 2812

    [11]

    Sun H Q, Zhang W H 2006 Chin. J. Energ. Mater. 14 16 (in Chinese) [孙海权, 张文宏2006 含能材料 14 16]

    [12]

    Lee E L, Tarver C M 1980 Phys. Fluids 23 2362

  • [1] 宋利伟, 石颖, 陈树民, 柯璇, 侯晓慧, 刘志奇. 地下黏弹性介质波动方程及波场数值模拟. 物理学报, 2021, 70(14): 149102. doi: 10.7498/aps.70.20210005
    [2] 刘富成, 刘雅慧, 周志向, 郭雪, 董梦菲. 双层耦合非对称反应扩散系统中的超点阵斑图. 物理学报, 2020, 69(2): 028201. doi: 10.7498/aps.69.20191353
    [3] 姜春华, 赵正予. 化学复合率对激发赤道等离子体泡影响的数值模拟. 物理学报, 2019, 68(19): 199401. doi: 10.7498/aps.68.20190173
    [4] 丁明松, 江涛, 董维中, 高铁锁, 刘庆宗, 傅杨奥骁. 热化学模型对高超声速磁流体控制数值模拟影响分析. 物理学报, 2019, 68(17): 174702. doi: 10.7498/aps.68.20190378
    [5] 李诗尧, 于明. 固体炸药爆轰的一种考虑热学非平衡的反应流动模型. 物理学报, 2018, 67(21): 214704. doi: 10.7498/aps.67.20172501
    [6] 华雪东, 王炜, 王昊. 考虑车与车互联通讯技术的交通流跟驰模型. 物理学报, 2016, 65(1): 010502. doi: 10.7498/aps.65.010502
    [7] 陈大伟, 王裴, 孙海权, 蔚喜军. 爆轰波对碰驱动平面锡飞层的动力学及动载行为特征研究. 物理学报, 2016, 65(2): 024701. doi: 10.7498/aps.65.024701
    [8] 王路, 徐江荣. 两相湍流统一色噪声法概率密度函数模型. 物理学报, 2015, 64(5): 054704. doi: 10.7498/aps.64.054704
    [9] 刘扬, 韩燕龙, 贾富国, 姚丽娜, 王会, 史宇菲. 椭球颗粒搅拌运动及混合特性的数值模拟研究. 物理学报, 2015, 64(11): 114501. doi: 10.7498/aps.64.114501
    [10] 陈石, 王辉, 沈胜强, 梁刚涛. 液滴振荡模型及与数值模拟的对比. 物理学报, 2013, 62(20): 204702. doi: 10.7498/aps.62.204702
    [11] 敬明, 邓卫, 王昊, 季彦婕. 基于跟车行为的双车道交通流元胞自动机模型. 物理学报, 2012, 61(24): 244502. doi: 10.7498/aps.61.244502
    [12] 白占国, 董丽芳, 李永辉, 范伟丽. 双层耦合Lengel-Epstein模型中的超点阵斑图. 物理学报, 2011, 60(11): 118201. doi: 10.7498/aps.60.118201
    [13] 尹增谦, 赵盼盼, 董丽芳, 房同珍. 有源等离子体在开放大气环境下的反应扩散过程研究. 物理学报, 2011, 60(2): 025206. doi: 10.7498/aps.60.025206
    [14] 赵啦啦, 刘初升, 闫俊霞, 蒋小伟, 朱艳. 不同振动模式下颗粒分离行为的数值模拟. 物理学报, 2010, 59(4): 2582-2588. doi: 10.7498/aps.59.2582
    [15] 卓士创, 闫长春. 可见光紫端748 THz处的高透过率负折射率材料模拟研究. 物理学报, 2010, 59(1): 360-364. doi: 10.7498/aps.59.360
    [16] 蔡利兵, 王建国. 介质表面高功率微波击穿的数值模拟. 物理学报, 2009, 58(5): 3268-3273. doi: 10.7498/aps.58.3268
    [17] 侯日立, 彭建祥, 经福谦. 一种计算金属剪切模量的本构模型:以Al为例. 物理学报, 2009, 58(9): 6413-6418. doi: 10.7498/aps.58.6413
    [18] 陈大年, 范春雷, 胡金伟, 谭华, 王焕然, 吴善幸, 俞宇颖. 高导无氧铜的高压与高应变率本构模型研究. 物理学报, 2009, 58(4): 2612-2618. doi: 10.7498/aps.58.2612
    [19] 尚华艳, 黄海军, 高自友. 基于元胞传输模型的可变信息标志选址问题研究. 物理学报, 2007, 56(8): 4342-4347. doi: 10.7498/aps.56.4342
    [20] 欧阳建明, 邵福球, 王 龙, 房同珍, 刘建全. 一维大气等离子体化学过程数值模拟. 物理学报, 2006, 55(9): 4974-4979. doi: 10.7498/aps.55.4974
计量
  • 文章访问数:  3122
  • PDF下载量:  666
  • 被引次数: 0
出版历程
  • 收稿日期:  2014-02-28
  • 修回日期:  2014-06-18
  • 刊出日期:  2014-11-05

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

  • 1. 北京应用物理与计算数学研究所, 北京 100094
    基金项目: 中国工程物理研究院科学技术发展基金(批准号:2013A0201010,2013B0101014)和国家自然科学基金(批准号:11272064)资助的课题.

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

English Abstract

参考文献 (12)

目录

    /

    返回文章
    返回