惰性气体和水蒸气对长直空间燃气爆炸超压及其振荡的抑制作用

惰性气体和水蒸气对长直空间燃气爆炸超压及其振荡的抑制作用

Inert Gas and Water Vapor Suppressing Overpressure and Its Oscillation of Gas Explosion in Long Straight Space

惰性气体和水蒸气对长直空间燃气爆炸超压及其振荡的抑制作用

1. 有限元模型

2. 结果与讨论

3. 结　论

1.1 计算模型

1.2 计算工况

1.3 材料参数

2.1 管道应力分析

2.2 管道位移分析

2.3 管道等效应变分析

2.4 管道振速分析

doi: 10.11858/gywlxb.20200654

计量

doi: 10.11858/gywlxb.20200654

计量

出版历程

作者简介:
刘　洋（1996－），男，硕士研究生，主要从事燃气爆炸荷载研究. E-mail：20145035@cqu.edu.cn

通讯作者:
李　展（1990－），男，博士，讲师，主要从事燃气爆炸灾害效应研究. E-mail：lz.9008@163.com

作者简介: 刘 洋（1996－），男，硕士研究生，主要从事燃气爆炸荷载研究. E-mail：20145035@cqu.edu.cn

通讯作者: 李 展（1990－），男，博士，讲师，主要从事燃气爆炸灾害效应研究. E-mail：lz.9008@163.com

刘洋; 李展; 方秦; 王森佩; 陈力

doi:10.11858/gywlxb.20200654

刘洋^1,,
李展^1,*, ,,
方秦¹,
王森佩¹,
陈力²

1.
陆军工程大学爆炸冲击防灾减灾国家重点实验室，江苏南京　210007
2.
东南大学爆炸安全防护教育部工程研究中心，江苏南京　211189

基金项目: 江苏省自然科学基金（BK20190571）；国家自然科学基金（52008392）

详细信息

中图分类号: O381; TD712
计量
- 文章访问数: 2004
- HTML全文浏览量: 1118
- PDF下载量: 30
出版历程
- 收稿日期: 2020-12-11
- 修回日期: 2021-01-13

LIU Yang^1
,,
LI Zhan^{1,*
, ,},
FANG Qin¹,
WANG Senpei¹,
CHEN Li²

1.
State Key Laboratory of Disaster Prevention and Mitigation of Explosion and Impact, Army Engineering University of PLA, Nanjing 210007, Jiangsu, China
2.
Engineering Research Center of Safety and Protection of Explosion & Impact of Ministry of Education, Southeast University, Nanjing 211189, Jiangsu, China

摘要

摘要: 长直空间燃气爆炸超压及其振荡将对人员和结构安全产生不利影响。为减轻燃气爆炸危害，基于CFD软件FLACS建立了长直管道空间燃气爆炸数值模型，并对模型进行了验证。利用已验证的数值模型，研究了添加不同体积分数CO₂、N₂和水蒸气的化学当量比CH₄/空气混合气体的爆炸，讨论了惰性气体和水蒸气的体积分数对爆炸超压及其振荡的影响，并对比了3种气体的抑爆效果。结果表明：CO₂、水蒸气和N₂的体积分数每增加10%，密闭管道气体爆炸的最终超压将分别下降81、47、65 kPa，尾端泄爆管道分别下降24、25、20 kPa，3种气体的体积分数分别为25%、26%、30%时，爆炸被完全抑制；CO₂、水蒸气和N₂均能有效抑制爆炸超压的振荡，压力振幅和压力振荡频率均随添加气体体积分数的增加而减小；CO₂对爆炸超压及其振荡的抑制效果最好，水蒸气次之，N₂最弱，这与3种气体的物理特性及其抑爆机理的差异有关。
- 水蒸气 /
- 惰性气体 /
- 燃气爆炸抑制 /
- FLACS /
- 压力振荡
Abstract: The overpressure and its oscillation of gas explosion in long straight space will endanger personnel and structures. In order to reduce the potential hazards, a numerical model of gas explosions in long straight tube was established based on the CFD software FLACS and was verified by using the existing experimental data. Based on the validated numerical models, the suppression of CO₂, N₂ and water vapor on the CH₄/air mixture explosions with stoichiometric concentration was investigated. The influence of the volume fraction of inert gas and water vapor on the explosion overpressure and its oscillation was considered and discussed. It is shown that for every 10% increase in the volume fraction of CO₂, water vapor and N₂, the final overpressure of the gas explosion in the closed tube drops by 81, 47 and 65 kPa, the final overpressure in the end-vented tube decreases by 81, 47 and 65 kPa. When the volume fractions of CO₂, water vapor and N₂ are 25%, 26%, and 30%, respectively, the explosion is completely suppressed. CO₂, water vapor and N₂ can effectively suppress the oscillation of explosive overpressure, both the pressure amplitude and the pressure oscillation frequency decrease with the increase of the volume fraction of added gas. CO₂ has the best suppression effect on explosion overpressure and its oscillation, followed by water vapor, and N₂ is the weakest. This phenomenon is related to the differences in the physical properties of the three gases and the suppression mechanisms.
- water vapor /
- inert gas /
- gas explosion suppression /
- FLACS /
- pressure oscillation

HTML全文

作为液体和气体长距离运输的一种重要方式，管道运输在国家经济发展和国民生活中发挥着重要作用。然而，随着城镇化进程的加速，城市管网系统密集分布，爆炸作用引起的管道安全问题受到国内外广泛关注^[1-3]。都的箭等^[4]通过实验研究发现，正对爆心管段背面受到很大的轴向拉应力作用，且管道受爆炸载荷的影响主要与爆心距有关。Ji等^[5]研究了X70钢管在局部爆炸载荷下的动力响应，发现管道的挠度和损伤程度随炸药量和接触面积的增大而增大，且壁厚对管道损伤和失效后的运动有重要作用。数值模拟是研究爆炸问题的一种重要方法，只要方法得当，模拟效果可与实际情况相吻合^[6-7]。为此，梁政等^[8]利用数值模拟方法研究了管道埋深、药量和管道壁厚因素对爆炸载荷下的埋地管道动力响应的影响。房冲^[9]通过模拟研究发现，在爆炸载荷下充水管道的变形量、位移和峰值压强都比内空管道小。余洋等^[10]采用野外实验与数值计算相结合的方法研究了初始条件对钢质方管在侧向局部爆炸载荷作用下损伤破坏效应的影响。

迄今为止，对爆炸载荷作用下焊缝区附近埋地钢管的动力响应的相关研究鲜有报道。基于此，以两种含Y型焊缝（坡口有2 mm余高焊缝和坡口无余高焊缝）的埋地X70钢管为例，采用有限元软件ANSYS/LS-DYNA，数值模拟研究爆炸载荷作用下焊缝区附近埋地X70钢管的动力响应规律，以期为埋地管线附近的爆破施工设计和埋地管线的安全防护提供一定的理论参考。

采用cm-g-μs单位制，建立由TNT炸药、黄土和焊接管道组成的计算模型，如图1所示。模型纵向长38.4 cm，管道中心到模型侧面的宽度为130.0 cm，模型整体高271.6 cm，其中：TNT炸药为边长14.0 cm的正方体，采用中心起爆方式；焊接管道为外径1 016.0 mm、壁厚14.6 mm的X70钢管。焊缝选取两种尺寸^[11]，分别为Y型坡口有余高(H = 2.0 mm)焊缝和Y型坡口无余高(H = 0)焊缝，如图2所示。为了提高计算收敛速度，将焊缝尺寸设计图进行适当的简化，简化模型如图3所示。两种焊缝均不考虑分层焊接工艺的影响，焊缝与管道采用共节点方式连接。

图 1 计算模型

Figure 1. Calculation model

下载: 全尺寸图片幻灯片

图 2 焊缝的设计尺寸

Figure 2. Design of weld size

下载: 全尺寸图片幻灯片

图 3 焊缝的简化模型

Figure 3. Simplified diagram of weld model

下载: 全尺寸图片幻灯片

考虑到计算模型的对称性，取1/2模型建模。炸药、黄土、管道及焊缝选用SOLID164六面体实体单元，用扫掠方式划分网格，并对焊缝位置进行网格细化处理。炸药和黄土采用欧拉网格，焊接管道和焊缝采用拉格朗日网格，运用任意拉格朗日-欧拉算法及管土间流固耦合算法模拟爆炸载荷作用下埋地焊接管道的动力响应。在土体外侧和底面设置透射边界条件，模型对称面施加对称约束。

为了初步揭示爆炸载荷作用下两种焊缝形式的埋地焊接管道的动力响应规律，选取药包尺寸为14.0 cm × 14.0 cm × 14.0 cm的TNT炸药，对埋深为1.5 m的焊缝有余高（H = 2.0 mm）管道（管道A）和焊缝无余高（H = 0）管道（管道B），在炸高分别为60.0、85.0和110.0 cm的3种条件下的6种工况进行模拟计算，如表1所示，其中，h_B为炸高。

表 1 计算工况

Table 1. Calculation conditions

Weld type	Buried depth of pipeline/m	Size of TNT/(cm × cm × cm)	h_B/cm
No weld reinforcement (H = 0)	1.5	14.0 × 14.0 × 14.0	60.0, 85.0, 110.0
Weld reinforcement (H = 2.0 mm)	1.5	14.0 × 14.0 × 14.0	60.0, 85.0, 110.0

下载: 导出CSV

| 显示表格

TNT炸药选用高能炸药模型（Mat_High_Explosive_Burn）和JWL状态方程定义。JWL状态方程表达式为

${p_{\rm{z}}}= A\left( {1 - \frac{\omega }{{{R_1}\nu }}} \right){{\rm{e}}^{ - {R_1}\nu }} + B\left( {1 - \frac{\omega }{{{R_2}\nu }}} \right){{\rm{e}}^{ - {R_2}\nu }} + \frac{{\omega E}}{\nu }$

(1)

式中：p_z为爆炸产物的压力，A、B、R₁、R₂、ω为TNT材料常数，v为爆炸产物的相对比容，E为炸药初始内能。炸药密度ρ_z、爆速D以及JWL状态方程参数见表2^[12]。

表 2 炸药材料参数^[12]

Table 2. Material parameters of explosive^[12]

ρ_z/(g·cm^–3)	D/(m·s^–1)	p/GPa	A/GPa	B/GPa	R₁	R₂	ω	E/(J·cm^–3)
1.58	6 930	21	373.77	3.75	4.15	0.90	0.35	6 000

下载: 导出CSV

| 显示表格

黄土选用泡沫模型（Mat_Soil_and_Foam）描述。该材料模型的应力屈服常数f为

$f = {S\!_{ij}}{\delta _{ij}}/2 - \left( {{a_0} + {a_1}{p_{\rm{t}}} + {a_2}{{p_{\rm{t}}}^2}} \right)$

(2)

式中：S_ij为土体材料的Cauchy偏应力张量，δ_ij为土体材料的Kronecker系数，a₀、a₁、a₂分别为土体摩擦角、土体黏聚力和土体爆炸动载效应的影响系数，p_t为土体压力。a₀、a₁、a₂由土工实验测得的内摩擦角和土壤黏聚力参数确定，土体密度ρ_t、剪切模量G、体积模量K等参数见表3^[13-14]。

表 3 土体材料参数^[13–14]

Table 3. Material parameters of soil^[13–14]

ρ_t/(g·cm^–3)	G/MPa	K/MPa	a₀/Pa²	a₁/Pa	a₂
1.8	41.14	87.87	2.12 × 10⁸	5.23 × 10³	3.22 × 10^–2

下载: 导出CSV

| 显示表格

X70钢管道和焊缝均采用双线性随动材料模型（Mat_Plastic_Kinematic）描述，遵循von Mises屈服准则，其表达式为

$\sigma =\begin{cases} {{E_{\rm{s}}}\varepsilon }&{\varepsilon \leqslant {\varepsilon _{\rm{e}}}}\\ {{\sigma _{\rm{y}}} + {E_{\rm{t}}}(\varepsilon - {\varepsilon _{\rm{e}}})}&{\varepsilon > {\varepsilon _{\rm{e}}}} \end{cases}$

(3)

式中：σ为应力；σ_y为屈服应力；E_s为弹性模量；E_t为切线模量，0 < E_t < E_s；ε为应变；ε_e为弹性极限应变。管道和焊缝的具体材料参数见表4^[15-18]，其中，μ为泊松比。

表 4 管道及焊缝材料参数^[15–18]

Table 4. Material parameters of pipe and weld bead^[15–18]

Material	ρ/(g·cm^–3)	μ	E_s/GPa	E_t/GPa	σ_y/GPa
X70-pipeline^[15–16]	7.90	0.3	210	13.5	0.48
Weld bead^[17–18]	7.25	0.3	220	15.3	0.55

下载: 导出CSV

| 显示表格

图4为边长14.0 cm的正方体TNT炸药爆炸时，炸高h_B为60.0 cm，埋深为1.5 m的两种X70管道焊缝附近的von-Mises应力云图。由图4可以看出：当传播时间为1 440 μs时，爆炸应力波阵面已经接触管道；当传播时间为1 600 μs时，焊缝有2.0 mm余高的管道A和焊缝无余高的管道B的最大应力增幅分别为81.4 MPa和43.0 MPa；当传播时间为1 920 μs时，管道A和管道B的最大应力均大于焊缝与管道的材料屈服应力，且应力沿迎爆面正对爆心位置向外扩展，其中管道A的应力呈“十”字形扩展，而管道B的应力以椭圆形向四周扩展；在3 520 μs时，管道应力集中主要沿裂缝位置发展，管道A和管道B的应力最大值分别为601.2 MPa和591.0 MPa；在6 080 μs时，管道A和管道B继续变形但应力减小，应力最大值分别减小到581.8 MPa和565.8 MPa；在9 120 μs时，管道A和管道B的应力集中基本消失。

图 4 不同时刻X70钢管道的von Mises应力

Figure 4. von Mises stress of X70 steel pipe at different moments

下载: 全尺寸图片幻灯片

图5和图6分别为两种管道外表面上正对爆心位置的焊缝与管道分界面处焊缝单元与管道单元的应力时程曲线。在管道受爆炸载荷作用阶段，两种管道的应力在大约480 μs内呈跳跃式上升。其主要原因是管道为瞬时受力，一部分爆炸能量使管道变形并向管道四周传递，导致焊缝与管道分界面处两个典型单元的应力呈降低趋势，此现象与图4的应力云图吻合。根据应力集中系数和余高关系的经验公式^[19]可得：管道A和管道B的应力集中系数分别为1.016和1.008，即随着余高增大，应力集中系数逐渐增大。对比图5和图6可知，管道A的焊缝单元应力峰值较高，应力下降趋势相对较陡。这也说明焊缝余高的存在使得焊缝与管道分界面的截面尺寸突变增大，从而导致焊缝有余高的焊接管道受应力集中的影响较大。在1 912 μs时，图5（管道A）和图6（管道B）的焊缝单元应力最大，分别约为560.0 MPa（焊缝的屈服强度为550 MPa）和545.6 MPa。同时，管道A的焊缝处首先达到管道屈服强度（480 MPa），按照von Mises屈服准则，管道A开始进入局部塑性变形阶段，此时管道B的应力尚未达到材料的屈服强度。

图 5 管道A(H = 2.0 mm)典型单元的应力时程曲线

Figure 5. Stress-time curves of classic element of A pipe (H = 2.0 mm)

下载: 全尺寸图片幻灯片

图 6 管道B(H = 0)典型单元的应力时程曲线

Figure 6. Stress-time curves of classic element of B pipe (H = 0)

下载: 全尺寸图片幻灯片

埋地X70管道的迎爆面和背爆面的最大位移如表5所示。从表5可知，由于爆炸冲击波的一部分能量在土中被耗散，且随着爆炸冲击波在土中传播距离的增大，两种焊缝形式管道的迎爆面和背爆面的最大位移均呈现减小的趋势。当炸高h_B从60.0 cm增加到85.0 cm以及从85.0 cm增加到110.0 cm时，管道A和管道B迎爆面的最大位移减小量分别为2.303 cm、0.715 cm和2.300 cm、0.572 cm，而管道A和管道B背爆面的最大位移减小量分别为0.391 cm、0.235 cm和0.373 cm、0.280 cm。两种焊缝形式管道迎爆面的最大位移减小量大于背爆面，这是由于爆炸冲击载荷在土中传播后直接作用于管道迎爆面，对管道迎爆面产生的影响较大，土体对管道背爆面具有一定的支撑作用，从而减小了管道背爆面位移。在相同炸高下管道A比管道B的最大位移大，且在炸高为60.0、85.0和110.0 cm时，两种焊缝形式管道迎爆面的最大位移差值分别为0.270、0.267和0.124 cm，即随着炸高的增大，两种焊缝形式的埋地焊接管道最大位移的差值逐渐减小，也说明当炸高较小时，管道A整体抵抗变形的能力弱于管道B。然而，随着炸高的增大，作用于管道的能量减小^[20]，管道塑性变形较小，使得这种现象逐渐模糊。

表 5 埋地X70管道的迎爆面和背爆面的最大位移

Table 5. Maximum displacement of explosion-front and explosion-back surfaces of buried X70 pipeline

Types of weld	h_B/cm	Maximum displacement/cm
Types of weld	h_B/cm	Explosion-front surface	Explosion-back surface
Weld reinforcement (H = 2.0 mm)	60.0	5.482	0.846
	85.0	3.179	0.455
	110.0	2.464	0.220
No weld reinforcement (H = 0)	60.0	5.212	0.943
	85.0	2.912	0.570
	110.0	2.340	0.290

下载: 导出CSV

| 显示表格

表6为两种不同类型焊缝的埋地焊接管道在不同炸高下的最大等效应变统计。从表6可知，管道A和管道B的最大等效应变均随炸高的增大而减小。当炸高从60.0 cm增大到85.0 cm时，管道A和管道B的最大等效应变分别减小约58.12%和61.13%；当炸高从85.0 cm增大到110.0 cm时，管道A和管道B的最大等效应变分别减小约45.92%和38.05%，在炸高相同时，管道A的最大等效应变大于管道B，且管道A的最大等效应变位于焊缝余高表面，而管道B的最大等效应变在焊缝与管道处一定范围内沿纵向分布。这在一定程度上说明管道B能更好地协调焊缝与管道分界处的应变，有利于保障焊缝与管道的局部协同变形性能。

表 6 不同炸高时埋地X70管道的最大等效应变

Table 6. Maximum effective strain of buried X70 pipeline with different blasting heights

Types of weld	h_B/cm	Peak effective strain/10^–3
Weld reinforcement (H = 2.0 mm)	60.0	9.937
	85.0	4.162
	110.0	2.251
No weld reinforcement (H = 0)	60.0	6.877
	85.0	2.673
	110.0	1.656

下载: 导出CSV

| 显示表格

表7为不同炸高下两种焊缝形式管道的迎爆面和背爆面处焊缝位置中心单元X方向的最大振动速度。从表7可以看出，两种焊缝形式管道的迎爆面和背爆面的最大振动速度均随着炸高增大而减小，且迎爆面的最大振速均大于背爆面。这说明迎爆面受爆炸地震波的影响较大。当炸高h_B为60.0、85.0和110.0 cm时，管道B的迎爆面的最大振动速度较管道A大，迎爆面差值分别为1.600、0.539和0.329 m/s，而背爆面差值在0.200 m/s以内。管道峰值速度随着管壁厚度的增大而减小^[12]，由于管道A增加了管道在焊缝位置的径向厚度，可将其视为管道焊缝位置的壁厚增大导致管道A的峰值振速减小。这说明管道A抵抗爆炸振动的性能优于管道B，且在炸高为60.0 cm时，管道A抵抗振动性能的优势较为明显。

表 7 埋地X70管道的迎爆面和背爆面最大振速

Table 7. Maximum vibration velocity of explosion-front and explosion-back surfaces of buried X70 pipeline

Types of weld	h_B/cm	Maximum vibration velocity/(m·s^–1)
Types of weld	h_B/cm	Explosion-front surface	Explosion-back surface
Weld reinforcement (H = 2.0 mm)	60.0	22.748	4.431
	85.0	9.316	2.817
	110.0	4.503	1.693
No weld reinforcement (H = 0)	60.0	24.348	4.294
	85.0	9.855	2.867
	110.0	4.832	1.746

下载: 导出CSV

| 显示表格

图7为不同炸高时两种焊缝形式的管道典型单元的速度时程曲线。当炸高h_B分别为60.0、85.0和110.0 cm时，管道A和管道B达到最大振速的时间分别为2 560 μs和2 560 μs、4 500 μs和4 600 μs、7 200 μs和7 200 μs，两种焊缝形式的管道达到最大振速的时间差值均在100 μs以内。这说明两种焊缝形式的管道达到最大振速的时间主要受炸高的影响，受焊缝形式的影响较小。

图 7 管道典型单元的速度时程曲线

Figure 7. Velocity-time curve of typical pipeline elements

下载: 全尺寸图片幻灯片

（1）当炸高为60.0 cm时，两种焊缝形式的埋地X70焊接管道在爆炸载荷作用下焊缝位置均出现应力集中，但焊缝有余高的管道受应力集中影响较大，且会先于焊缝无余高管道进入屈服阶段。

（2）当炸高为60.0～110.0 cm时，由于爆炸载荷直接作用于迎爆面，且管土间的相互作用对管道背爆面具有一定的支撑作用，两种焊缝形式管道迎爆面的最大位移均大于背爆面的最大位移。当炸高为60.0、85.0 cm时，焊缝有余高的管道整体抵抗变形的能力明显弱于焊缝无余高的管道。

（3）焊缝无余高管道较焊缝有余高管道在焊缝与管道分界处的应变更为协调，能更好地保障焊缝与管道的局部协同变形性能。

（4）在相同的爆炸载荷下，焊缝有余高管道抵抗振动的性能优于焊缝无余高管道。药量相同条件下，相对于焊缝形式，炸高对含焊缝区管道的最大振速起主要作用。

图 1 管道模型及其网格划分

Figure 1. Tube model and its meshing

下载: 全尺寸图片幻灯片

图 2 模型验证工况^[38]

Figure 2. Cases of model verification^[38]

下载: 全尺寸图片幻灯片

图 3 模型验证（测点1）

Figure 3. Numerical model verification (pressure sensor 1)

下载: 全尺寸图片幻灯片

图 4 混合气体体积示意图

Figure 4. Diagram of premixed gas volume

下载: 全尺寸图片幻灯片

图 5 不同CO₂体积分数条件下的超压时程曲线

Figure 5. Overpressure-time history curves of CO₂ with different volume fractions

下载: 全尺寸图片幻灯片

图 6 最终超压及升压时间随CO₂体积分数的变化

Figure 6. Change of final overpressure and pressure rising time with CO₂ volume fraction

下载: 全尺寸图片幻灯片

图 7 不同CO₂体积分数下超压振幅随时间的变化

Figure 7. Variation of overpressure amplitude of different volume fractions of CO₂ with time

下载: 全尺寸图片幻灯片

图 8 不同CO₂体积分数下超压振荡频率随时间的变化

Figure 8. Variation of overpressure oscillation frequency of different volume fractions of CO₂ with time

下载: 全尺寸图片幻灯片

图 9 不同N₂体积分数条件下的超压时程曲线

Figure 9. Overpressure-time history curves of N₂ with different volume fractions

下载: 全尺寸图片幻灯片

图 10 最终超压及升压时间随N₂体积分数的变化

Figure 10. Change of final overpressure and pressure rising time with N₂ volume fraction

下载: 全尺寸图片幻灯片

图 11 不同N₂体积分数下超压振幅随时间的变化

Figure 11. Variation of overpressure amplitude of different volume fractions of N₂ with time

下载: 全尺寸图片幻灯片

图 12 不同N₂体积分数下超压振荡频率随时间的变化

Figure 12. Variation of overpressure oscillation frequency of different volume fractions of N₂ with time

下载: 全尺寸图片幻灯片

图 13 不同水蒸气体积分数条件下的超压时程曲线

Figure 13. Overpressure-time history curves of N₂ with different volume fractions

下载: 全尺寸图片幻灯片

图 14 最终超压及升压时间随水蒸气体积分数的变化

Figure 14. Change of final overpressure and pressure rising time with water vapor volume fraction

下载: 全尺寸图片幻灯片

图 15 不同水蒸气体积分数下超压振幅随时间的变化

Figure 15. Variation of overpressure amplitude of different volume fractions of water vapor with time

下载: 全尺寸图片幻灯片

图 16 不同水蒸气体积分数下超压振荡频率随时间的变化

Figure 16. Variation of overpressure oscillation frequency of different volume fractions of water vapor with time

下载: 全尺寸图片幻灯片

图 17 不同CO₂体积分数条件下的超压时程曲线

Figure 17. Overpressure-time history curves of CO₂ with different volume fractions

下载: 全尺寸图片幻灯片

图 18 最终超压和升压时间随CO₂体积分数的变化

Figure 18. Change of final overpressure and pressure rising time with CO₂ volume fraction

下载: 全尺寸图片幻灯片

图 19 不同CO₂体积分数下超压振幅随时间的变化

Figure 19. Variation of overpressure amplitude of different volume fractions of CO₂ with time

下载: 全尺寸图片幻灯片

图 20 不同CO₂体积分数下超压振荡频率随时间的变化

Figure 20. Variation of overpressure oscillation frequency of different volume fractions of CO₂ with time

下载: 全尺寸图片幻灯片

图 21 不同N₂体积分数条件下的超压时程曲线

Figure 21. Overpressure-time history curves of N₂ with different volume fractions

下载: 全尺寸图片幻灯片

图 22 最终超压及升压时间随N₂体积分数的变化

Figure 22. Change of final overpressure and pressure rising time with N₂ volume fraction

下载: 全尺寸图片幻灯片

图 23 不同N₂体积分数下超压振幅随时间的变化

Figure 23. Variation of overpressure amplitude of different volume fractions of N₂ with time

下载: 全尺寸图片幻灯片

图 24 不同N₂体积分数下超压振荡频率随时间的变化

Figure 24. Variation of overpressure oscillation frequency of different volume fractions of N₂ with time

下载: 全尺寸图片幻灯片

图 25 不同水蒸气体积分数条件下的超压时程曲线

Figure 25. Overpressure-time history curves of water vapor with different volume fractions

下载: 全尺寸图片幻灯片

图 26 最终超压及升压时间随水蒸气体积分数的变化

Figure 26. Change of final overpressure and pressure rising time with water vapor volume fraction

下载: 全尺寸图片幻灯片

图 27 不同水蒸气体积分数下超压振幅随时间的变化

Figure 27. Variation of overpressure amplitude of different volume fractions of water vapor with time

下载: 全尺寸图片幻灯片

图 28 不同水蒸气体积分数下超压振荡频率随时间的变化

Figure 28. Variation of overpressure oscillation frequency of different volume fractions of water vapor with time

下载: 全尺寸图片幻灯片

图 29 密闭管道内CO₂、N₂和水蒸气对爆炸的抑制作用对比

Figure 29. Comparison of suppression of CO₂, N₂ and water vapor on the explosion in the closed tube

下载: 全尺寸图片幻灯片

图 30 尾端泄爆管道内CO₂、N₂和水蒸气对爆炸的抑制作用对比

Figure 30. Comparison of suppression of CO₂, N₂ and water vapor on the explosion in the end-vented tube

下载: 全尺寸图片幻灯片

表 1 添加气体的体积分数（密闭/泄爆）

Table 1. Volume fraction of added gas (closed/end-vented) %　

CO₂ N₂ Water vapor (H₂O) CO₂ N₂ Water vapor (H₂O)

0 0 0 20 25 20
4 5 4 24 30 24
8 10 8 25 31 26
12 15 12 27
16 20 16

下载: 导出CSV

参考文献(41)

[1]	李建国, 郭昭胜, 张永生, 等. 室内燃气爆炸对混凝土结构的破坏及有限元分析 [J]. 消防科学与技术, 2018, 37(4): 563–566. doi: 10.3969/j.issn.1009-0029.2018.04.041 LI J G, GUO Z S, ZHANG Y S, et al. The concrete structural failure caused by internal gas explosion and its finite element analysis [J]. Fire Science and Technology, 2018, 37(4): 563–566. doi: 10.3969/j.issn.1009-0029.2018.04.041
[2]	闫秋实, 刘晶波, 伍俊. 典型地铁车站内爆炸致人员伤亡区域的预测研究 [J]. 工程力学, 2012, 29(2): 81–88. YAN Q S, LIU J B, WU J. Estimation of casualty areas in subway station subjected to terrorist bomb [J]. Engineering Mechanics, 2012, 29(2): 81–88.
[3]	王波, 杜扬, 齐圣, 等. 油气爆炸在细长密闭管道内的振荡传播特性 [J]. 振动与冲击, 2017, 36(17): 97–103, 126. WANG B, DU Y, QI S, et al. Oscillation propagation characteristics of gasoline-air mixture explosion in elongated closed tubes [J]. Journal of Vibration and Shock, 2017, 36(17): 97–103, 126.
[4]	XIAO H H, MAKAROV D, SUN J H, et al. Experimental and numerical investigation of premixed flame propagation with distorted tulip shape in a closed duct [J]. Combustion and Flame, 2012, 159(4): 1523–1538.
[5]	LI H W, GUO J, YANG F Q, et al. Explosion venting of hydrogen-air mixtures from a duct to a vented vessel [J]. International Journal of Hydrogen Energy, 2018, 43(24): 11307–11313.
[6]	YAO Z F, DENG H X, ZHAO W L, et al. Experimental study on explosion characteristics of premixed syngas/air mixture with different ignition positions and opening ratios [J]. Fuel, 2020, 279: 118426.
[7]	LV P F, ZHANG J X, JU M H, et al. Influence of branch pipes on the deflagration characteristics of methane in confined space [J]. Process Safety Progress, 2020, 40(1): e12152.
[8]	李国庆, 杜扬, 王波, 等. 点火位置对管道内油气泄压爆炸超压特性影响 [J]. 振动与冲击, 2017, 36(24): 204–212. LI G Q, DU Y, WANG B, et al. Effects of ignition position on overpressure characteristics of vented gasoline-air mixture explosion in a pipe [J]. Journal of Vibration and Shock, 2017, 36(24): 204–212.
[9]	PHYLAKTOU H, ANDREWS G E. Gas explosions in long closed vessels [J]. Combustion Science and Technology, 1991, 77(1/2/3): 27–39.
[10]	郑立刚, 朱小超, 于水军, 等. 浓度和点火位置对氢气-空气预混气爆燃特性影响 [J]. 化工学报, 2019, 70(1): 408–416. ZHENG L G, ZHU X C, YU S J, et al. Effect of concentration and ignition position on characteristics of premixed hydrogen-air deflagration [J]. CIESC Journal, 2019, 70(1): 408–416.
[11]	SEARBY G. Acoustic instability in premixed flames [J]. Combustion Science and Technology, 1992, 81(4/5/6): 221–231.
[12]	XING H D, XU Q M, SONG X Z, et al. The effects of vent area and ignition position on pressure oscillations in a large L/D ratio duct [J]. Process Safety and Environmental Protection, 2020, 135: 166–170.
[13]	朱传杰, 林柏泉, 江丙友, 等. 瓦斯爆炸在封闭管道内冲击振荡特征的数值模拟 [J]. 振动与冲击, 2012, 31(16): 8–12, 17. doi: 10.3969/j.issn.1000-3835.2012.16.002 ZHU C J, LIN B Q, JIANG B Y, et al. Numerical simulation on oscillation and shock of gas explosion in a closed end pipe [J]. Journal of Vibration and Shock, 2012, 31(16): 8–12, 17. doi: 10.3969/j.issn.1000-3835.2012.16.002
[14]	ZHU C J, LIN B Q, JIANG B Y, et al. Numerical simulation of blast wave oscillation effects on a premixed methane/air explosion in closed-end ducts [J]. Journal of Loss Prevention in the Process Industries, 2013, 26(4): 851–861.
[15]	韦世豪, 杜扬, 王世茂, 等. 不同形状受限空间内油气爆燃特性的实验研究 [J]. 中国安全生产科学技术, 2017, 13(5): 41–47. WEI S H, DU Y, WANG S M, et al. Experimental study on deflagration characteristics of gasoline-air mixture in confined space with different shapes [J]. Journal of Safety Science and Technology, 2017, 13(5): 41–47.
[16]	韦世豪, 杜扬, 王世茂, 等. 非均匀容积式受限空间油气爆炸超压与火焰特征 [J]. 后勤工程学院学报, 2017, 33(4): 30–36. doi: 10.3969/j.issn.1672-7843.2017.04.006 WEI S H, DU Y, WANG S M, et al. Overpressure and flame characteristics of gasoline-vapor explosion in the non-uniform volumetric confined space [J]. Journal of Logistical Engineering University, 2017, 33(4): 30–36. doi: 10.3969/j.issn.1672-7843.2017.04.006
[17]	李毅. 管道中氢-空和甲烷-空预混火焰传播与压力振荡研究 [D]. 焦作: 河南理工大学, 2015.
[18]	路长, 李毅, 潘荣锟. 管道氢气-空气预混气体爆炸特征的试验研究 [J]. 安全与环境学报, 2016, 16(3): 38–42. LU C, LI Y, PAN R K, et al. Experimental study on explosion tendency of hydrogen-air premixed gases in the duct [J]. Journal of Safety and Environment, 2016, 16(3): 38–42.
[19]	王亚磊, 郑立刚, 于水军, 等. 约束端面对管内甲烷爆炸特性的影响 [J]. 爆炸与冲击, 2019, 39(9): 139–148. WANG Y L, ZHENG L G, YU S J, et al. Effect of vented end faces on characteristics of methane explosion in duct [J]. Explosion and Shock Waves, 2019, 39(9): 139–148.
[20]	WANG Z R, NI L, LIU X, et al. Effects of N₂/CO₂ on explosion characteristics of methane and air mixture [J]. Journal of Loss Prevention in the Process Industries, 2014, 31: 10–15. doi: 10.1016/j.jlp.2014.06.004
[21]	LI M, XIAO Y, DENG J, et al. Effect of CO₂ on explosion limits of flammable gases in goafs [J]. Mining Science and Technology, 2010, 20(2): 193–197.
[22]	LI M H, XU J C, WANG C J, et al. Thermal and kinetics mechanism of explosion mitigation of methane-air mixture by N₂/CO₂ in a closed compartment [J]. Fuel, 2019, 255(1): 115747.
[23]	张迎新, 吴强, 刘传海, 等. 惰性气体N₂/CO₂抑制瓦斯爆炸实验研究 [J]. 爆炸与冲击, 2017, 37(5): 906–912. ZHANG Y X, WU Q, LIU C H, et al. Experimental study on coal mine gas explosion suppression with inert gas N₂/CO₂ [J]. Explosion and Shock Waves, 2017, 37(5): 906–912.
[24]	李凌飞. 甲烷爆炸特性及其抑爆技术研究[D]. 太原: 中北大学, 2012.
[25]	任韶然. 惰性及特种可燃气体对甲烷爆炸特性的影响实验及分析 [J]. 天然气工业, 2013, 33(10): 110–115. doi: 10.3787/j.issn.1000-0976.2013.10.019 REN S R. An experimental study of effects of insert and special flammable gases on methane’s explosion characteristics [J]. Natural Gas Industry, 2013, 33(10): 110–115. doi: 10.3787/j.issn.1000-0976.2013.10.019
[26]	王华, 葛岭梅, 邓军. 惰性气体抑制矿井瓦斯爆炸的实验研究 [J]. 矿业安全与环保, 2008, 35(1): 4–7. doi: 10.3969/j.issn.1008-4495.2008.01.002 WANG H, GE L M, DENG J. Experimental study of using inert gas to suppress mine gas explosion [J]. Mining Safety & Environmental Protection, 2008, 35(1): 4–7. doi: 10.3969/j.issn.1008-4495.2008.01.002
[27]	李成兵. N₂/CO₂/H₂O抑制甲烷爆炸化学动力学机理分析 [J]. 中国安全科学学报, 2010, 20(8): 88–92. doi: 10.3969/j.issn.1003-3033.2010.08.014 LI C B. Chemical kinetics mechanism analysis of N₂/CO₂/H₂O suppressing methane explosion [J]. China Safety Science Journal, 2010, 20(8): 88–92. doi: 10.3969/j.issn.1003-3033.2010.08.014
[28]	XIE Y L, WANG J H, ZHANG M, et al. Experimental and numerical study on laminar flame characteristics of methane oxy-fuel mixtures highly diluted with CO₂ [J]. Energy and Fuels, 2013, 27: 6231–6237. doi: 10.1021/ef401220h
[29]	秦文茜, 王喜世, 谷睿, 等. 超细水雾作用下瓦斯的爆炸压力及升压速率 [J]. 燃烧科学与技术, 2012, 18(1): 90–95. QIN W Q, WANG X S, GU R, et al. Methane explosion overpressure and overpressure rise rate with suppression by ultra-fine water mist [J]. Journal of Combustion Science and Technology, 2012, 18(1): 90–95.
[30]	康泉胜, 李振明, 王睿, 等. 超细水雾对管内丙烷爆炸火焰抑制效果的试验研究 [J]. 安全与环境学报, 2015, 15(6): 111–114. KANG Q S, LI Z M, WANG R, et al. Experimental study on the inhibition effects of the ultra-fine water mist on the propane explosion flame in the tube [J]. Journal of Safety and Environment, 2015, 15(6): 111–114.
[31]	高旭亮. 超细水雾抑制甲烷爆炸实验与数值模拟 [D]. 大连: 大连理工大学, 2014.
[32]	刘丹, 司荣军, 李润之. 环境湿度对瓦斯爆炸特性的影响 [J]. 高压物理学报, 2015, 29(4): 307–312. doi: 10.11858/gywlxb.2015.04.011 LIU D, SI R J, LI R Z. Ambient humidity on explosion characteristics of methane/air mixture [J]. Chinese Journal of High Pressure Physics, 2015, 29(4): 307–312. doi: 10.11858/gywlxb.2015.04.011
[33]	GEXCON A S. FLACS V10. 8 user’s manual [Z]. Bergen, Norway: GEXCON AS, 2017.
[34]	LAUNDER B E, SPALDING D B. The numerical computation of turbulent flows [J]. Computer Methods in Applied Mechanics and Engineering, 1974, 3(2): 269–289.
[35]	陈晓坤, 丁园月, 程方明, 等. CO₂对矿井多组分可燃性气体抑爆特性的影响 [J]. 煤炭科学技术, 2015, 43(3): 43–47. CHEN X K, DING Y Y, CHENG F M, et al. Influence of CO₂ on explosion suppression characteristics of multicomponent flammable gas in coal mine [J]. Coal Science and Technology, 2015, 43(3): 43–47.
[36]	罗振敏, 王涛, 程方明, 等. 小尺寸管道内二氧化碳抑制甲烷爆炸效果的实验及数值模拟 [J]. 爆炸与冲击, 2015, 35(3): 393–400. doi: 10.11883/1001-1455-(2015)03-0393-08 LUO Z M, WANG T, CHENG F M, et al. Experimental and numerical studies on the suppression of methane explosion using CO₂ in a mini vessel [J]. Explosion and Shock Waves, 2015, 35(3): 393–400. doi: 10.11883/1001-1455-(2015)03-0393-08
[37]	张印, 赵东风, 刘义. 基于FLACS的CH₄/CO₂/air混合气爆炸参数分析 [J]. 中国安全生产科学技术, 2016, 12(9): 36–40. ZHANG Y, ZHAO D F, LIU Y. Analysis on explosion parameters of CH₄/CO₂/air mixed gas based on FLACS [J]. Journal of Safety Science and Technology, 2016, 12(9): 36–40.
[38]	LI Z, CHEN L, YAN H C, et al. Gas explosions of methane-air mixtures in a large-scale tube [J]. Fuel, 2021, 285: 119239.
[39]	杨春丽, 刘艳, 胡玢, 等. 氮气和水蒸气对瓦斯爆炸基元反应的影响及抑爆机理分析 [J]. 高压物理学报, 2017, 31(3): 301–308. doi: 10.11858/gywlxb.2017.03.012 YANG C L, LIU Y, HU F, et al. Effect of nitrogen an water vapor on methane-air mixture explosion elementary reaction suppression mechanism [J]. Chinese Journal of High Pressure Physics, 2017, 31(3): 301–308. doi: 10.11858/gywlxb.2017.03.012
[40]	袁渭兰, 吕卫民, 贾忠湖. 气体动力学[M]. 北京: 科学出版社, 2013: 55.
[41]	卓士创, 董慎行. 弹簧振子在定常干摩擦阻尼作用下的振动 [J]. 大学物理, 2001(6): 17–21, 33. doi: 10.3969/j.issn.1000-0712.2001.06.005 ZHUO S C, DONG S X. The oscillation of a spring oscillator under constant dry frictional damping [J]. College Physics, 2001(6): 17–21, 33. doi: 10.3969/j.issn.1000-0712.2001.06.005

施引文献

资源附件(0)

访问统计

点击查看大图

图(30) / 表(1)

文章访问数: 2004
HTML全文浏览量: 1118
PDF下载量: 30

1. 有限元模型
1.1 计算模型
1.2 计算工况
1.3 材料参数
2. 结果与讨论
2.1 管道应力分析
2.2 管道位移分析
2.3 管道等效应变分析
2.4 管道振速分析
3. 结　论

Water vapor (H₂O)

Water vapor (H₂O)