基于MATLAB平台的凝聚相爆炸模型的分析与程序实现_王颖

5 7
#
(
1+
$Ps /
6 7
$P
P s/
0
P
0
)
1/
#2C0Fra bibliotek( 10)
Qs= Q0# ( 1+
6 7
$P s/
P 0) /
( 1+
1 7
$Ps / P0)
( 11)
Ts=
T 0( 1+
$P P0
s
)#
(
1+
1 7
$Ps / P0) / ( 1+
6 7
$Ps / P 0)
( 12)
C= C0( 1+
摘 要 主要利用 MATLA B 程序, 实现凝聚相爆炸事故( CPE) 的模型分析与计算, 通过对凝聚相爆炸模型的分析和后果 模拟, 得出死亡、重伤的区域半径, 为安全评价及分析的定量评价提供技术支撑。
关键词 MATLAB 凝聚相爆炸模型 死亡半径 重伤半径
Analysis of Condensed Phase Detonation Model Based on MATLAB Software and the Realization of the Procedure WANG Ying WU Ying CHEN N ing- ning CAI Xiang- ping
随着计算机技术的迅猛发展, 大量的工程软件可 以应用 于安全分析的相应环节之 中, MATLAB 就 是其中 之一[ 1] 。本 文主要论述 MATLAB 程序在凝聚相爆炸事故分析中的应 用, 通过对凝聚相 爆 炸( CPE) 模 型的 分 析 和后 果 模拟, 得出 死
亡、重伤的区域半径, 进而对伤害区域进行划分。 1 凝聚相爆炸( CPE) 模型分析计算
图 2 无输入值时冲击波超压与距离的三维等高线 当输入 c= 340, T = 288 时, p = 10 300 Pa, 其冲击波超 压 与距离的三维等高线如图 3。 2. 2 数据分析 安全分析工作中仅 仅依靠图 形来表 现结 果缺乏 相应 的 精确性, 难以表示具体的伤害 区域。这里 可以依 靠 MATLAB 的强大计算功能来准确地表达爆炸的伤害范围, 计算的精 度 可以根据实际需要灵 活调节。 根据冲 击波超 压对建 筑物 的 破坏作用和对人体的伤害作用表, 选取以 爆炸中心 300 m 以 内的范围进行区域 划分, 步 长初步 定为 1 m, 此 精确 度在 伤 害区域划分中已足够, 现截取部分数据如表 1。
受严重伤害或死亡, 其内径为 0, 外径记为 R 0. 5 , 表 示外圆 周 外人员因冲击波 作用 导致 肺 出血 死亡 的概 率为 50% , 它 与
爆炸当量间的关系由下式确定:
R0. 5= 13. 6( WTNT/ 1 000) 0. 37
( 2)
WTNT= E/ QTNT
( 3)
k
E E= QB, IWI
本文重点是研究爆炸的预防控制措施, 因此将凝 聚相爆
炸( CPE) 模型用 于模 拟评 价与 分析 烟花 爆竹 爆炸 事故 的后 果严重度、危险等级、灾害影响和破坏范围。
爆炸具有多种破坏形式, 如冲击波破坏, 热辐射危害, 一 次、二次碎片破坏, 爆炸产物毒气危害, 以及由爆炸引 起的火 灾等。其中 , 最危险、破坏力最强、破坏区域最大的是 冲击波
$Ps ) 1/ 2#
1+ 1+
1 7
$Ps / P0
6 7
$Ps / P0
1/ 2
( 13)
在标准状态下, P 0= 0. 1 MPa, Q0 = 0. 125 kg# s2/ m4 , T0 =
288 K, C0 = 340 m/ s。
2 MATLAB 实现凝聚相爆炸( CPE) 后果分析
MATLAB 强大的计算和绘图 能力, 使得对 事故后果 的计
基于以上基本条件, 用 MATLAB 可以得出冲击波超压和距 离的关系( 见图 1) , 以此来划分伤害后果区域,部分程序如下:
图 1 冲击波超压与距离函数 E= 16272; P0= 101300;
# 39 #
E0= P0/ E; Q1= ( E0) . ^( 1/ 3) ; R= 0: 10: 400; Q= R* Q1; Pa= P0. ^( 0. 137. * ( Q. ^( - 3) ) + 0. 199. * ( Q. ^( - 2) ) + 0. 269. * ( Q. ^( - 1) ) - 0. 119) ; plot( R, Pa) ; 为了更加直观地标志冲击波波阵面上的超压 , 本文采 用 了三维等高线来表示 爆炸产 生的冲 击波超压 和距离 之间 的 关系, 通过程序 得出了冲击 波超压 和距离 之间 三维等 高线。 其源程序代码如下: E= 16272; P0= 101300; E0= P0/ E; Q1= ( E0) . ^( 1/ 3) ; [ x , y] = meshgrid( - 300: 3: 300, - 300: 3: 300) ; R= ( x. ^2+ y. ^2) . ^1/ 2; Q= R* Q1; Pa= P0. ^( 0. 137. * ( Q. ^( - 3) ) + 0. 199. * ( Q. ^( - 2) ) + 0. 269. * ( Q. ^( - 1) ) - 0. 119) ; z= Pa/ P0; plot3( x, y, z) ; c= input(c c= c) ; T = input(cT = c) ; p= input(cp= c) ; contour3( x, y, Pa, c) ; contour3( x, y, Pa, T) ; contour3( x, y, Pa, p) ; xlabel(cR/ mc) ; ylabel(cR/ mc) ; zlabel(cP/ Pac) ; 当未输参数时, 其图形如图 2。
# 38 #
工业安全与环保 Industrial Saf ety and Environment al Protect ion
2009 年第 35 卷第 11 期 November 2009
基于 MATLAB 平台的凝聚相爆炸 模型的分析与程序实现
王颖 伍颖 陈宁宁 蔡香萍
( 中国地质大学( 武汉) 工程学院 武汉 430074)
于轻伤区半径计算, 可取 $Ps = 17 kPa。
波阵面上的其他参数( 如波 阵面的 传播速 度 U, 粒 子速
度 Vs , 密度 Qs , 温度 T s , 音速 C 等) , 均 可以由 理想气体 的激
波关系, 通过超压 $Ps 求出[ 5] 。
U= C0( 1+
6 7
$P s)
1/2
(9)
Vs =
( 5)
综上可以得出
WTNT= 1. 8 W f
( 6)
式中, WTNT 为 爆 源的 TNT 当 量, kg ; E 为 爆 源 的 总能 量, J; QTNT为 TNT 爆热, 可取为 4. 52 MJ/ kg; Wf 为爆源的质量, kg。
如果认为该圆周内 没有死亡 的人数 正好 等于圆 周外 死
依据上文 1. 2 节, 对 于重伤 区半径, 取 $Ps = 44 kPa; 对 于轻伤区半径计算, 取 $Ps = 17 kPa。由 表 1 可得, 轻 伤半径 R= 180 m, 死亡半径依据式( 3) 计算, 可得 R= 21. 8 m。 3 结论
( Engineering College , China University of Geosciences Wuhan 430074) Abstract This art icle mainly use MATLAB procedure to realize the condensed phase explosion( CPE) model analysis and calculation and obt ain the death region radius and the severely wounded region radius by CPE model analysis and the consequence simulat ion. This result can provide technical support for the safety evaluat ion and quant it ative evaluat ion analysis. Keywords MATLAB CPE model death radius severely wounded radius
# 40 #
图 3 冲击波超压与距离的三维等高线 > > R= 0: 10: 400;
Pa= 101300. ^( 0. 137. * ( ( 0. 137. * R) . ^( - 3) ) + 0. 199. * ( ( 0. 137. * R) . ^( - 2) ) + 0. 269. * ( ( 0. 137. * R) . ^( - 1) ) - 0. 119)
Q1/ 3/ R= Q01/ 3/ R0
(1)
则在不同距离上的超压是相同的。式中, R 为爆炸测 量点到
爆炸点的距离; R0 为 1 000 kg 的 TNT 爆 炸测 量点 到爆 炸点
的距离; Q 为某次炸药的 TNT 当量; Q 0 为 1 000 kg TNT。 爆炸相似律为我们 从小当量 爆炸所 获得 的结果 换算 到
$Ps / P0= 0. 137 Z- 3+ 0. 119 Z- 2 + 0. 269 Z- 1- 0. 119
( 1< $Ps / P0< 10)
(7)
Z= R# ( P0/ E) 1/3
(8)