Consequences assessment of explosions in pipes using coupled FEM-SPH method

Explosions often lead to destruction of equipment, which is a difficult problem including complicated fluid-solid interactions. Most traditional CFD methods cannot synchronously solve the movements of fluids and large deformation and fracture of solids because such problem is usually accompanied with constantly moving-and-changing boundary conditions. In this paper, a coupled Finite Element Method-Smoothed Particle Hydrodynamics (FEM-SPH) method was proposed to simulate the dynamic processes of explosions in pipes. The propagation of blast wave and the fracture of pipe were captured in every timestep, where the energy dissipation caused by plastic deformation and crack propagation were fully considered. A rate-dependent failure criterion for high-strain-rate load conditions was employed in the numerical simulation, which was presented in our previous work and has been verified in the dynamic fracture behavior of steels for pressure vessels and pipes. In addition, a simpler formula was proposed to describe the attenuation of blast wave outside the pipe and the consequences caused by the explosions were assessed. Results revealed the interaction between blast wave and pipe, the leakage of detonation products, the attenuations of peak overpressures outside the pipe and the corresponding consequences at different distances. It is found that when considering the energy consumption during plastic deformation and crack propagation in coupled FEM-SPH method, the assessment results are more rational than that without considering such energy consumption.     

Propagation and intensity of blast wave from hydrogen pipe rupture due to internal detonation

The coupled fluid-structure-rupture model was developed to study the propagation and intensity of blast wave from hydrogen pipe rupture due to internal detonation. The dynamic rupture of pipe and propagation of blast wave were well coupled together in every timestep during the simulation. The numerical model was validated with experiments in terms of both typical rupture profiles and blast overpressures. Results reveal that crack branching of pipe can dramatically increase the rupture opening rate which controls the intensity and shape of the resultant blast wave. Due to the process of crack initiation and extension, the blast wave out of the pipe first forms and then is strengthened by the subsequent compression waves. This makes the maximum peak overpressure appears at a certain standoff distance above the rupture. Despite consuming some percentages of energy, the dynamic rupture of pipe generally presents positive effects (up to 2–3 times) on the blast wave intensity along the jetting direction due to the convergence effect of rupture opening on the release of internal high-pressure gas. Finally, through defining normalized overpressure and impulse based on the same hydrogen detonation in open spaces, the quantitative influences of pipe rupture on the blast wave intensity in cases of different detonation pressures and standoff distances are clarified.     

FMECA-模糊综合分析法在自旋滑车桥轴风险评价中的应用

FMECA法是一种定性的分析方法,它通过对目标系统下各单元的功能、潜在故障模式、故障原因、故障率及严重性进行分析,找出危险性较高的故障模式和设计中的薄弱环节。模糊综合分析法可以将定性的评价结果以定量的形式表达,但其评价结果具有模糊性。提出了一种将FMECA法与模糊综合分析法相结合的评价方法,既可以减少模糊综合分析评价结果的模糊性,又可以将FMECA法的定性描述以定量的形式表达,使评价结果更加客观、合理。并以大型游乐设施中自旋滑车的桥轴为例,运用FMECA-模糊综合分析法对其进行风险评价,计算结果表明其风险等级较低,与实际检验结果相符。     

管道中氢气-空气爆轰的多级泄爆数值模拟研究

为了研究管道内氢气的爆燃转爆轰及其抑制过程,对单个障碍物管道中氢气-空气混合物燃爆过程以及多级泄爆进行了二维数值模拟。基于氢气-空气19步详细化学反应动力学机理,以及k-ε湍流模型、概率密度函数输运方程和同位网格SIMPLE算法,采用计算流体软件Fluent进行模拟。结果表明:密闭管道无泄爆时,在距点火端1.5 m左右爆燃转为爆轰;泄爆口的位置对管道内氢气-空气预混气体的爆炸参数有重要影响,泄爆口位于管道中部时,能降低管道内爆轰超压,泄爆效果较好;位于管道中部单个泄爆口泄爆时,有效降低爆轰超压,管道中部设置2个泄爆口时,能通过压力和混合气体的泄放将管道中已经发生的爆轰衰减为爆燃;当有3个泄爆口泄爆时,管道中没有发生爆轰,达到良好的泄爆效果。     

Propagation Characteristics of Gas Explosion in Linked Vessels Based on DDT Criteria

To study the occurrence conditions and propagation characteristics of deflagration to detonation transition (DDT) in linked vessels, two typical linked vessels were investigated in this study. The DDT of the methane–air mixture under different pipe lengths and inner diameters was studied. Results showed that the CJ detonation pressure of the methane–air mixture was 1.86 MPa, and the CJ detonation velocity was 1987.4 m/s. Compared with a single pipe, the induced distance of DDT is relatively short in the linked vessels. With the increase in pipeline length, DDT is more likely to occur. Under the same pipe diameter, the DDT induction distance in the vessel–pipe–vessel structure is shorter than that in the vessel–pipe structure. With the increase in pipeline diameter, the length of the pipe required to form the DDT is reduced. For linked vessels in which detonation formed, four stages, namely, slow combustion, deflagration, deflagration to detonation, and stable detonation, occurred in the vessels. Moreover, for a pipe diameter of 60 mm and a length of 8 m, overdriven detonation occurred in the vessel–pipe–vessel structure.     

加卸载作用下裂隙岩芯渗透性变化规律研究

针对光滑平行板裂隙加载渗流规律难以准确反映实际工程裂隙岩体卸荷渗流特性的问题,利用水力压裂方法制作了天然岩芯裂隙,开展了裂隙试件受轴压、围压加载和卸载作用的渗流试验,并将裂隙试件的CT扫描图像导入ABAQUS软件,对粗糙裂隙受载时的应力和变形规律进行了数值模拟。结果表明:裂隙面的法向应力是影响裂隙渗流的主要因素;在加载阶段,粗糙裂隙面受载产生接触应力致使接触部分产生变形和破坏,裂隙接触刚度增大,裂隙宽度非线性减小,渗透率降低速率越来越小;在卸载阶段,弹性变形逐步恢复,裂隙宽度缓慢增加,渗透率先缓慢增大,当载荷卸到一定值后,在孔隙压力作用下裂隙张开、渗透率突然增大。研究结果表明了粗糙裂隙与光滑平行板裂隙对载荷敏感性不同的机理,反映了实际工程裂隙岩体卸荷渗流特性。     

3D Voronoi等效晶质模型在岩石破坏细观研究中的应用

基于3D Voronoi图提出一种能够反映岩石矿物晶粒结构的三维精细化建模方法,分析了晶粒排列均匀化程度及大小等微观结构的变化对模型宏观力学性质的影响程度,并进行了晶粒结构优化设计。通过模拟粉砂质板岩单轴压缩和巴西劈裂的破坏过程,以及与室内试验对比验证了其适用性,并从细观角度揭示了岩石的破坏机理。研究结果表明:通过裂纹生长速率大致可以表征室内试验中弹性变形、非弹性变形和峰值后宏观破坏阶段;岩样在单轴压缩和巴西劈裂破坏过程中以拉伸断裂为主,在宏观上表现为脆性断裂;岩样宏观断裂呈现沿晶和穿晶断裂的组合形态,弹性变形阶段裂纹的萌生和扩展以沿晶断裂为主,而穿晶断裂的扩展和贯通往往导致局部失稳,表现为非弹性变形阶段。     

深水钻井热交换作用下的井壁稳定性分析

为了研究深水钻井过程中热交换作用对井壁稳定性的影响规律,建立了钻井液、钻柱、海水和地层各介质热交换作用下温度控制方程,根据热弹性理论计算温度变化对井周应力的影响,结合地层强度准则与拉伸破坏准则确定深水井壁坍塌压力与破裂压力。结果表明,钻井液循环对地层产生的冷却作用能够降低井壁坍塌压力与破裂压力,但破裂压力降低幅度比坍塌压力的降低幅度稍大。因此,尽管钻井液的冷却效果能够抑制井壁坍塌,但应注意低温度的井壁更易形成诱导缝,甚至引发地层漏失。计算结果对深水钻井作业具有指导意义。     

含热熔孔洞缺陷的埋地聚乙烯管道应力分析及寿命预测

为了研究聚乙烯（PE）管道热熔孔洞缺陷与寿命之间的关系,基于Ansys软件对埋地状态下的含有热熔孔洞缺陷的PE管道进行应力分析,得到了含热熔孔洞缺陷的PE管道在不同内压下的最大Mises应力、环向应力和径向应力的变化数据;依据应力分析数据,使用Matlab编写程序,根据Suleiman双曲本构模型对PE管道进行了寿命预测。研究结果表明:含有热熔孔洞缺陷的PE管道,其最大应力随管道内压的增大而增大,寿命随着管道内压和缺陷体积的增大而减少,含缺陷管道寿命和内压关系可使用双对数函数来进行描述;使用应力分析和寿命预测相结合的方法,可以得到不同缺陷的PE管道寿命-内压关系式。     

初始温度和初始压力对气体爆轰参数影响的研究

利用已有的气体爆炸模型和包含初始压力、初始温度的气体爆轰参数的计算公式,从理论上研究初始压力和初始温度对气体爆轰参数的影响情况。使用VisualBasic语言编写计算程序,将计算值与文献值进行对比,具有较好的一致性。以甲烷-空气混合物为例,计算在98000Pa,280～400K及298K,0.1～0.5MPa的气体爆轰参数。计算结果表明,初始压力一定,混合物的爆轰压随初始温度的升高而减小,爆轰波速增大;初始温度一定,混合物的爆轰压随初始压力的增大而增大,爆轰波速基本不变;在初始温度和初始压力两个影响因素中,初始压力对混合物爆轰参数的影响明显大于初始温度。     

考虑岩石抗拉强度受围压影响的地层破裂压模型

预测地层破裂压力是地下流体注入工程中的关键问题,关系到工程的稳定性和安全性。根据地测资料,岩石的抗拉强度会随着围压产生变化,然而目前常见的地层破裂压预测模型均未考虑该因素。基于厚壁圆筒法,针对2种砂岩开展不同围压条件下的砂岩致裂试验,根据试验数据,计算各组试验中试样的抗拉强度,并拟合围压与岩石抗拉强度的关系函数。研究岩石抗拉强度随围压变化的规律,提出能在实际地层中表征该规律的参数,基于该参数对现有的地层破裂压模型进行修正,并结合实测数据验证修正模型的可靠性。研究结果表明:岩石的抗拉强度随着围压的增加呈现明显的线性正相关关系,低围压下,新模型的预测结果与实测数据较为接近。随着围压增加,二者的差异逐渐增大,说明选取的特征参数在低围压下能较好地描述抗拉强度受围压的影响效应,而高围压下参数的取值需进一步研究。     

An experimental study of flame acceleration and deflagration to detonation transition in representative process piping

The paper summarizes the results of experimental tests and accompanying analyses to investigate the factors that govern flame acceleration and potential transition to detonation in a relatively long unobstructed piping system. The overall aim of the work was to obtain sufficient experimental data so as to be able to develop and evaluate methodologies for classifying and predicting potential detonation flame acceleration and deflagration to detonation transition (DDT) hazard in industrial process pipes and mixtures. The present results show that the flame acceleration process in an unobstructed pipe exhibit three distinct phases: an initial establishment phase; a second rapid acceleration phase and a final transition to detonation phase. Test results with ethylene indicate that the acceleration process is not sensitive to initial pressure (all other parameters remaining constant) but can be sensitivity to initial pipe wall temperature or possibly mixture humidity. The presence of bends increases the local rate of turbulent combustion, an effect attributed to the additional turbulence generated downstream of the bend. For straight pipes, detonation was only observed to develop for hydrogen–air and ethylene–air mixtures. Detonation was not observed with methane, propane or acetone as fuel in the present piping apparatus.     

爆炸冲击波作用下圆形通风设施的动态破坏特性

为了研究瓦斯爆炸对通风设施造成的严重破坏及其导致的通风系统紊乱与灾情迅速扩展问题,基于LS-DYNA有限元软件,模拟研究巷道中瓦斯爆炸冲击波作用下圆形通风设施的动态破坏特性,分析应力、应变、速度、位移的动态特征及其破坏过程,并对圆形通风设施动态破坏机理进行初步探索。研究结果表明:爆炸冲击波作用下,圆形通风设施正、反面出现压应力与拉应力分布圆环,且压应力与拉应力峰值处于在外部受约束边界径向圆环附近区域与圆心位置,应变、速度和位移最大值均分布在次区域;爆炸冲击波作用下,通风设施表面形成多圈的正向和反向的应变相互交替,呈现W型的褶皱变形;爆炸冲击波作用下,环向裂纹集中在受约束的边界区域附近,发生脆性破坏,其他区域会出现较多的径向裂纹,可发生韧性破坏,整个断裂过程是脆性与韧性断裂的混合型断裂。     

Effect of ignition position on the run-up distance to DDT for hydrogen-air explosions

The method described in this paper enabled reliable and accurate positioning of an overdriven detonation by calculation of shock wave velocities (detonation and retonation) for hydrogen explosions in a closed 18 m long horizontal DN150 pipe. This enabled an empirical correlation between the ignition position and the run-up distance to DDT to be determined. It was shown that the initial ability of the flame to expand unobstructed and the piston-like effect of burnt gas expanding against the closed end of the tube contributed to initial flame acceleration and hence were able to affect the run-up distance to overdriven detonation. Flame speeds and rates of initial pressure rise were also used to explain how these two competing effects were able to produce a minimum in the run-up distance to DDT. The shortest run-up distance to DDT, relative to the ignition position, for this pipe and gas configuration was found when the ignition position was placed 5.6 pipe diameters (or 0.9 m) from the closed pipe end. The shortest run-up distance to DDT relative to the end of the pipe was recorded when the ignition source was placed 4.4 pipe diameters or 0.7 m from the pipe end.     

一起氨制冷系统管道断裂事故的原因分析

制冷系统的库房回气管道发生了断裂失效,断裂位置在环焊缝中心且断口平整,焊缝存在未焊透缺陷。管道采用管束整体发泡保温,增加了管道的刚性,不利于管道的位移补偿。通过对管道材质进行分析,化学成分和力学性能符合设计要求。对管道受力状态进行分析并进行应力计算,结果表明,管道焊接接头未焊透,削弱了焊接接头的承载能力。管道在"冷缩"时不能补偿而产生较大的拉应力,导致在管道强度薄弱的焊接接头处发生断裂失效。     

惰性物对爆炸的遏制

利用爆炸激波管技术研究了部分惰性物对爆炸特性的影响。用多通道辐射高温计测量了氮气对汽油爆轰温度的影响。研究表明，增加氮气含量可以较明显地降低汽油爆轰温度。采用化学当量配比氢氧混合物爆轰产生水蒸气的方法研究了水蒸气对汽油爆轰特性的影响。研究表明，水蒸气能明显降低爆炸压力，水蒸气压力增加到0．1MPa时可导致爆炸熄灭；硝基甲烷和氧气混合物中充人氮气后爆炸压力明显下降，采用光多通道分析系统(OMA谱仪)和多台单色谱仪的光谱测量结果表明，反应中间产物的CH3O、CH辐射强度迅速衰减，反应衰竭。     

Experimental study on detonation flame penetrating through flame arrester

In-line detonation flame arresters are important safety apparatus to prevent group tank fires caused by the spreading of fire through vapor connection lines. In this study, a DN50 experimental apparatus aimed at the detonation flame penetration characteristics and failure mechanisms in a flame arrester was set up, and a series of experiments were carried out with 6.6% C2H4 and air mixture. Pressure, and velocity of flame penetrating through flame arrester housing and filters were analyzed. Experimental results showed that the attenuation of pressure and velocity was proportional to the thickness of the filters. Two failure modes of the fire-extinguishing process in the flame arrester were captured directly with a high-speed camera. In Mode I, the detonation flame could go straight through the flame arrester filters when the filters were too thin. In Mode II, when the filters were not sufficiently thick, the remained shock wave pressure of detonation flame was still several times of the initial pressure and could rise sharply at the downstream contraction section, resulting in that the flammable gas at the downstream transition section could be compressed and reignited even the flame had been extinguished by filters. These conclusions are helpful to reveal the nature of failure modes of fire-extinguishing process and design flame arresters with high fire-resisting performance by structure improved.     

天然气管道局部等效应力及塑性变形评估

为了研究腐蚀及地面运动对埋地天然气管线安全性的协同影响,以X80管道为研 究对象,模拟腐蚀缺陷及土壤力作用于管道之上,利用有限元方法对有腐蚀缺陷与预应 变情况下的管道局部等效应力及塑性变形进行评估,结果表明腐蚀缺陷的深度对局部应 力和应力分布影响非常明显,在失效压力预测中起着决定性作用。随着腐蚀深度的增加 ,应力集中增强,导致内表面和外表面的等效应力大小进一步分化,腐蚀深度的增加对 管道内表面的等效应力的影响很大,但对有效塑性应变的影响却不大。模拟管道上施加 有纵向应变的土壤力,不论拉伸与压缩的情况下,都会降低管道的失效压力,在施加拉 伸预应变下的管道失效压力小于压缩预应变下的。塑性变形首先发生在外表面处,并扩 展到腐蚀缺陷相邻区域,管道内表面也具有一定的塑性变形,但强度低。     

储气井螺纹连接安全性有限元模拟分析及实验研究

储气井是一种新型压力容器,其螺纹连接的安全性目前既没有经过充分的研究论证,也没有经过实践的充分检验.通过实验和有限元模拟方法,对储气井螺纹连接的安全性进行了分析研究,得出了以下结论:(1)储气井模拟构件在极限内压实验过程中,在管体发生强度失效前,螺纹连接接头部位先发生泄漏失效,表明螺纹连接接头是储气井的薄弱环节;(2)储气井在极限内压条件下的破坏模式为“漏而不破”,这一特性有益于保障储气井的安全;(3)螺纹沟槽对于井管管体而言是一种“损伤”,但其并没有增高螺纹连接接头部位的结构应力,原因是因为接箍起到了加强作用;(4)螺纹连接接头部位螺纹牙底存在应力集中,可能会导致疲劳失效,建议设计时予以重点考虑,内压疲劳实验结果显示螺纹连接接头的抗疲劳性能满足储气井一般工况要求.     

土压平衡矩形顶管施工引起的地表沉降规律研究

为了分析与掌握矩形顶管施工引起的地表沉降变形规律,提出基于MIDAS-GTS有限元软件的分析方法,首先,通过工程实例验证MIDAS-GTS有限元软件在分析顶管施工引起地表沉降变形的有效性;其次,通过影响因素的敏感性分析确定不同影响因素作用下的地表沉降变形特性,这些影响因素主要包括开挖面支护压力、侧摩阻力、地层结构特性等;最后,通过地表沉降变形特性分析得出矩形顶管的适用范围。研究结果表明:基于MIDAS-GTS有限元软件的地表沉降规律分析方法和结果,可为矩形顶管施工过程提供相关参考和借鉴。