- 1、本文档共5页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
岩石中爆炸波传播的数值模拟.pdf
第 17 卷 第 1 、2 期 计 算 物 理 Vol. 17 ,No. 12
2000 年 3 月 CHIN ESE J OURNAL OF COMPU TA TIONAL PHYSICS Mar. , 2000
( )
[文章编号] 1001246X 2000
岩石中爆炸波传播的数值模拟
王铁良 , 张建鑫 , 韩学安
( 西北核技术研究所 ,陕西 西安 710024)
[摘 要] 采用一维球对称流体弹塑性模型模拟了一系列 TN T 填实爆炸引起的周围岩石动力学
参量的变化 ,根据计算结果拟合了岩石中峰值应力、粒子速度、加速度及位移的近似公式 ,并与一
些文献中提供的实测结果进行了比较 ,说明数值模拟可为爆破工程、安全评估提供参考。
[ 关键词] 化学爆炸 ;弹塑性流体 ;数值模拟
[ 中图分类号] O35 ;O241 [文献标识码] A
1 物理模型
考虑球形爆室 ,用一维球对称流体弹塑性模型模拟爆炸波在岩石中的传播。计算中所用
炸药为浇铸 TN T ,密度 1600 kg ·m - 3 ,岩石为某种花岗岩 ,所取的药量分别为 30 kg 、100 kg 、
500 kg 、1000 kg 。计算的对象由两部分组成:爆室部分和岩石部分 ,两部分通过爆室壁耦合。
爆室中的状态量由 B KW 状态方程求出 , 以此结果作为边界条件,求解岩石中的动力学参量 ,
再将算出的爆室壁的新位置作为爆室部分计算的条件。
11 爆室中状态参量的计算
爆室中的状态参量由B KW 方程给出 ,B KW 方程是 Los Alamos 科学实验室的固态爆轰产
物状态方程 ,对 TN T 炸药 ,其公式如下
( ) 2 ( ) 3 ( ) 4
ln p = A + B ln V + C ln V + D ln V + F ln V
其中, A = - 3 66524588562 , B = - 2 4671426161 , C = 2 284615733 ×10 - 1 ,
D = 6 0497692516 ×10 - 2 ,F = - 1 93025884483 ×10 - 2 . p 的单位为 102 GPa , V 的单
位为 10 - 3 m3 ·kg - 1 .
爆室的初始比容取炸药的初始比容
V 0 = 1/ ρ
0
其中, ρ 为炸药的初始密度。
0
将爆室压强作为岩石部分计算的边界条件, 计算岩石中的应力场和粒子运动, 经过 Δt 时
间后, 空腔比容 V 为
4 3
π
V = R / W
3
这里, R 是 t 时刻的爆室半径。
[ 收稿日期] ; [ 修订日期]
( )
[ 作者简介] 王铁良 1966~ ,男 ,陕西富平 ,助理研究员 ,从事计算物理方面的研究 ,西安 69 信箱 15 分箱 ,710024.
© 1994-2006 China Academic Journal Electronic Publishing House. All rights reserved.
第 1 、2 期 王铁良等 :岩石中爆炸波
文档评论(0)