爆炸冲击波在建筑群中传播规律的数值模拟研究

为了研究建筑群的布局特性对爆炸冲击波传播的干扰作用及爆炸对周围环境的影响,笔者基于任意拉格朗日多物质流固耦合算法,对空气采用ALE网格,爆轰产物采用JWL状态方程,利用LS-DYNA程序对TNT理想爆炸源的爆炸进行了数值模拟,并与经验公式的计算结果进行了对比,得出了数值模拟方法可行的结论。在该基础上,以一个建筑小区为例,对其遭受爆炸后的响应作了进一步模拟研究。结果表明,建筑物的密度和布局方式对爆炸场的数值与形态都有敏感作用,因此,在城市规划建筑布局阶段,应考虑建筑密度和布局对爆炸冲击波的影响。     

A numerical study of the evolution of the blast wave shape in rectangular tunnels

When the explosion of condensed materials occurs in square or circular cross-section tunnel, the subsequent blast wave reveals two patterns: three-dimensional close to the explosive charge and one-dimensional far from the explosion. Pressure decays for these two patterns have been thoroughly studied. However, when the explosion occurs in rectangular cross-section tunnel, which is the most regular geometry for underground networks, the blast wave exhibits a third, two-dimensional, patterns. In order to assess the range of these three patterns, several numerical simulation of blast waves were carried out varying the width and the height of the rectangular cross-section as well as the mass of the charge. Laws are presented to localize the transition zones between the 3D and the 2D patterns, and between the 2D and the 1D patterns, as functions of non-dimensional width and height. The numerical results of the overpressure are compared to existing 3D and 1D laws. An overpressure decay law is proposed to represent the 2D pattern. Knowing the two transition zones and the overpressure decays within these zones, an algorithm is presented to efficiently predict an overpressure map. This algorithm is validated by comparison with experimental data.     

燃料空气炸药一次引爆爆炸效应实验研究

为了提高燃料空气炸药（FAE）爆炸威力,设计制备了不同相态的FAE,并采用光测和电测方法,开展了开放空间下FAE一次引爆对比实验研究。结果表明:固态、液固混合及液态FAE被一次引爆后,均存在引爆中心装药、抛撒燃料、点火和爆炸4个阶段,云雾存续时间均长于等质量TNT装药的相应值,随着距离增加,固态FAE爆炸场超压变化规律为“减少-增加-减少”,其值高于等质量TNT装药的相应值;液态、液固混合FAE爆炸场超压变化规律为“增加-减少”,远场超压高于等质量TNT装药的相应值。     

建构筑物周边爆破施工安全控制距离及分析方法研究

为研究在居民区附近开展大面积爆破施工的安全可行性以及安全控制距离标准,采用数值模拟方法,分析了地铁车辆段爆破施工对周边高压线及居民楼的影响,根据分析结果,给出了安全控制距离及关键控制指标,说明采用数值模拟的方法适用于类似工况的可行性。另外,研究发现,在居民楼附近进行地铁车辆段爆破施工风险较高,建议在距建构筑物100 m范围内采用静态爆破或机械开挖的方式进行施工。     

Case study on the catastrophic explosion of a chemical plant for production of m-phenylenediamine

On March 21, in Xiangshui County, Yancheng city, China, a catastrophic explosion occurred at Tianjiayi Chemical Co., Ltd. that caused 78 deaths and injured more than 617 individuals, with a property loss of as much as US$ 100 million. An explosion crater with a diameter of 98 m was generated. An extensive consequence analysis based on the crater size and physical effects for the estimation of the quantity of exploded dinitrobenzene (DNB) was performed. By using the position recorded in the satellite photo and the layout of the plant, the crater site was determined to be located near a warehouse where dinitrobenzene was stored. From the scaling law associating the TNT equivalent with crater diameter, the TNT equivalent of exploded ditrobenzene was determined to be approximately 708 tons. Based on the results of the consequence analysis performed by three methodologies, it can be concluded that the Tianjiayi Chemical Co. warehouse not only stored a substantially large quantity of dinitrobenzene but was also located too close to adjacent structures without an adequate separation distance. Judging by the effect of explosive blast on personnel, a fatality radius was determined to be about 400 m from the explosion center.     

Numerical and reduced-scale experimental investigation of blast wave shape in underground transportation infrastructure

When an explosion occurs in a tunnel, the study of the blast wave quickly becomes complicated, owing to the multiple propagation patterns of the blast wave (incident wave, regular and Mach reflections) and to the geometrical conditions. Considering this problem, two patterns can be revealed. Near the explosive, the well-known free-field pressure wave can be observed. After multiple reflections on the tunnel's walls, this overpressure behaves like a one-dimensional (1D) wave. One aim of this paper is to determine the position of this transition spherical-to-planar wave propagation in a tunnel using both numerical and reduced-scale experiments, and thereby validate the dedicated law established in a previous work.For this purpose, a detonation of TNT in a tunnel with a cross-section of up to 55 m² is considered. Results show good agreement between the numerical simulations and experiments. The transition zone between the three-dimensional (3D) and the 1D wave is well detected. An application to a simplified subway station is also investigated which shows that significant planar waves can be transmitted to the neighboring stations via the junction tunnels.     

Small-scale physical explosions in shock tubes in comparison with condensed high explosive detonations

A series of small-scale experiments involving physical explosions in a 1.6 l pressure vessel was carried out. Explosions were initiated by spontaneous rupture of an aluminium membrane on one side of the vessel at a pressure in the range 1–1.2 MPa. The pressure waves released were measured at different distances along two separate shock tubes, one 10 m long and 200 mm in diameter (closed at one end by the high pressure vessel) and the other 15 m long and 100 mm in diameter.TNT equivalency was used for predicting the blast wave characteristics after vessel rupture. TNT equivalency was used because equations for prediction of peak pressure and impulse of the blast wave in 1-D geometry after detonations of condensed explosives are known. Some experiments with an equivalent amount of real explosive were carried out for comparison with the theoretical and experimental data obtained. The applicability of the TNT equivalency method presented for calculations of maximum pressure and shock wave impulse generated after rupture of the pressure vessel in 1-D geometry is discussed.     

One-step estimation method and nomogram to predict safety distances of pressurized gas pipelines from blast sources

Pressurized gas pipelines may be endangered by accidental or intentional activities causing excess ground vibrations leading to failure. Thus, many papers have been published on this issue, yet most of these efforts suffer from complexity and inability to ascertain that all engaged factors are taken into account. These factors are attributed, namely, to explosion source, explosion-ground interaction, energy transfer through the ground, ground-pipeline interaction, and response of the pipeline. The deterministic method developed previously by the author has derived useful prediction mathematical formulas that can be used to calculate, either the stand-off distance of purposefully used explosives, or the maximum quantity of explosives to be used at a given distance, or the equivalent stress of a pipeline induced by an explosive charge at a given standoff distance. The method can also be applied to accidental blast events in factories or other activities. Yet, the complex formulas developed are hard to be used in the field and a simpler approach is needed to facilitate and accelerate in-field estimations. Thus, a nomogram was prepared to make the safety estimations easy and fast and taking into account all engaged critical factors starting from the explosion source through the integrity of the structure.     

Evaluation of the power of the distributed blast type explosive

Many works have been done on power evaluation of explosives and some evaluation methods presented. However, because of the differences of the explosion characteristics between distributed blast and common condensed explosive, a more reasonable way for evaluating the power of distributed blast is needed. In the paper, for a given range in space, a TNT equivalency method is proposed. Both the physical and mathematical meaning of power are considered in the method. Instead of giving equivalency of some individual points, TNT equivalent value for a certain space range can be obtained based on the method.     

Measurement of gas detonation blast loads in semiconfined geometry

The ignition and explosion of combustible vapor clouds represents a significant hazard across a range of industries. In this work, a new set of gas detonations experiments were performed to provide benchmark blast loading data for non-trivial geometry and explosion cases. The experiments were designed to represent two different accident scenarios: one where ignition of the vapor cloud occurs shortly after release and another where ignition is delayed and a fuel concentration gradient is allowed to develop. The experiments focused on hydrogen-air and methane-oxygen detonations in a semiconfined enclosure with TNT equivalencies ranging from 9 g to 1.81 kg. High-rate pressure transducers were used to record the blast loads imparted on the interior walls of a 1.8 m × 1.8 m × 1.8 m test fixture. Measurements included detonation wave velocity, peak overpressure, impulse, and positive phase duration. A comparison of the pressure and impulse measurements with several VCE models is provided. Results show that even for the most simplified experimental configuration, the simplified VCE models fail to provide predictions of the blast loading on the internal walls of the test fixture. It is shown that the confinement geometry of the experiment resulted in multiple blast wave reflections during the initial positive phase duration portion of the blast loading, and thus, significantly larger blast impulse values were measured than those predicted by analytical models. For the pressure sensors that experienced normally-reflect blast waves for the initial blast impulse, the Baker-Strehlow and TNT equivalency models still struggled to accurately capture the peak overpressure and reflected impulse. The TNO multi-energy model, however, performed better for the case of simple normally-reflected blast waves. The results presented here may be used as validation data for future model or simulation development.     

城市危险源爆炸灾害应急损失评估研究

采用冲击波伤害-破坏超压准则,首先确定了不同当量TNT炸药爆炸引起的人员伤害区域和建筑物破坏区域;再通过确定危险物爆炸所产生热量值的大小,将各种危险物料转化为相应当量的TNT炸药,可用于评估各种不同危险物料爆炸所造成的损失。同时,结合GIS软件平台,编制了爆炸灾害评估模拟程序。在软件平台中以图形的形式给出了设定条件下人员伤害区域及建筑物破坏区域的范围,根据区域内的人员分布和建筑物造价、财产多少,给出了定量的人员伤亡和经济损失应急评估。     

爆炸应力波在层状节理岩体中的传播规律试验研究

为了研究爆炸应力波在岩体中的传播规律,基于相似理论,开展了爆炸应力波在层状节理岩体中衰减的实验研究。研究表明,节理阻碍 应力波的传播且存在各向异性;对于垂直或小角度入射节理的应力波,在垂直方向和纵向的衰减很明显并且幅度较大,在横向衰减不明显;对 于斜入射节理的应力波,横向和垂直方向没有明显的衰减,纵向的衰减很明显并且幅度较大;应力波的合成峰值速度与炮孔密集度之间的关系 密切,当炮孔间距超过一定距离时,实际测得的速度时程为最近药包爆破产生的应力波,并非总药量产生的应力波的叠加。     

TNT当量法预测某石化设备爆炸后果评价

以烃化反应器为模拟发生泄漏事故源,泄漏的物料形成爆炸性气体,发生爆炸事故,采用TNT当量法预测爆炸事故后果,通过公式,计算出因爆炸冲击波导致肺出血而引起死亡的概率为0.5的半径.     

Using computational fluid dynamics (CFD) for blast wave predictions

Explosions will, in most cases, generate blast waves. While simple models (e.g., Multi Energy Method) are useful for simple explosion geometries, most practical explosions are far from trivial and require detailed analyses. For a reliable estimate of the blast from a gas explosion it is necessary to know the explosion strength. The source explosion may not be symmetric; the pressure waves will be reflected or deflected when hitting objects, or even worse, the blast waves may propagate inside buildings or tunnels with a very low rate of decay. The use of computational fluid dynamics (CFD) explosion models for near and far field blast wave predictions has many advantages. These include more precise estimates of the energy and resulting pressure of the blast wave, as well as the ability to evaluate non-symmetrical effects caused by realistic geometries, gas cloud variations and ignition locations. This is essential when evaluating the likelihood of a given leak source as cause of an explosion or equally when evaluating the potential risk associated with a given leak source for a consequence analysis.In addition, unlike simple methods, CFD explosion models can also evaluate detailed dynamic effects in the near and far field, which include time dependent pressure loads as well as reflection and focusing of the blast waves. This is particularly valuable when assessing actual near-field blast damage during an explosion investigation or potential near-field damage during a risk analysis for a facility. One main challenge in applying CFD, however, is that these models require more information about the actual facility, including geometry details and process information. Collecting the necessary geometry and process data may be quite time consuming. This paper will show some blast prediction validation examples for the CFD model FLACS. It will also provide examples of how directional effects or interaction with objects can significantly influence the dynamics of the blast wave. Finally, the challenge of obtaining useful predictions with insufficient details regarding the geometry will also be addressed.     

Numerical simulation of far-field blast loads arising from large TNT equivalent explosives

Large TNT equivalent explosions usually arise from accidents occurring during the transportation, storage, and manufacturing of chemicals relevant to process industries. The blast wave generated by the explosion will spread and interact with the surrounding factories and storehouses, damaging the building structures within several kilometers and causing significant casualties and property losses. This study aims to develop an efficient numerical simulation method to predict blast loads to estimate the consequences of accidents involving far-field free air bursts or surface burst explosions. Before its interaction with the interested target, a blast wave is generated in the numerical model by specifying the initial and boundary conditions of the disturbed air. Based on empirical data of incident overpressure, an explicit formula to calculate the air particle velocity is derived from the governing equations of a perfect inviscid gas. A simplified path line method is proposed to calculate the air density. The proposed method is applied to the LS-DYNA CESE solver to simulate the blast loads on building structures in the far field. Validations against empirical data and experiments indicate that the proposed method is sufficiently accurate for engineering applications and, through a case study, presents a more efficient performance than the LOAD_BLAST_ENHANCED (LBE) and mapping methods.     

Blast wave clearing behavior for positive and negative phases

The purpose of the research was to improve prediction of response of buildings to blast waves by including the negative phase and considering clearing of both positive and negative phases. Commonly used structural design practices, which trace their origins to military design manuals, often ignore the negative phase as well as positive phase clearing. For high explosive threats, this approach is conservative in most circumstances. However, negative phase clearing had not previously been studied for blast waves, and the implications for structural response had not been evaluated. This paper presents results of modeling negative phase blast clearing behavior for a typical blast wave and discusses the differences from positive phase clearing. The implications of including positive and negative phase clearing in building blast damage analysis are also investigated through single-degree-of-freedom (SDOF) analyses.Blast waves from explosion sources like a vapor cloud explosion (VCE), pressure vessel burst or high explosive exhibit both positive and negative phases, and the relative magnitude of the positive and negative phases varies among explosion sources and the specific circumstances of each source. A fully reflected blast wave is produced if an incident blast wave were to strike an infinitely tall and wide wall in a normal orientation. Both the positive and negative phases of the blast wave are enhanced by the reflection process. However, when an incident blast wave strikes a wall of finite size in a normal orientation, rarefaction waves are created at the edges of the wall, and the rarefactions sweep down from the roof and inward from sides. The rarefaction waves result in a clearing effect for both the positive and negative phases.Clearing relieves some of the applied blast load on the reflected wall for the positive phase. However, this is not always the case for the negative phase. As shown by the results presented in this paper, clearing may either relieve or enhance the applied negative phase blast load, depending on the duration of the blast wave and the wall dimensions.The impact of negative phase clearing on structural response for generic building components was also investigated. Nonlinear SDOF methods were used to characterize response in terms of peak positive and negative displacements. It was found that the influence of the negative phase is significant and the peak structural response can occur during negative (outward) displacement.     

Computer simulation of shock waves transmission in obstructed terrains

Generation and transmission of blast waves in real terrains is of major importance for risk analysis procedures involving accidental explosion scenarios. The problem arises from the impact of overpressure wave on people and structures that may be lethal or catastrophic under certain conditions. In this paper, a CFD simulation of shock wave propagation in obstructed terrain is attempted. Overpressure histories as well as a series of critical parameters, namely the positive and negative peak overpressure, the arrival time, and the positive and negative phase duration at specific points within the domain were obtained during the simulation. Their comparison with experimental measurements from field-scale high explosive blast tests performed by HSE showed a reasonably good agreement indicating that CFD computer programs provide reliable tools for estimating explosive shocks in complex terrains.     

液体火箭共底破裂爆炸安全设防距离

针对航天发射场一旦发生低温推进剂泄漏而导致火箭爆炸,会对人员和财产造成重大损失的问题,采用TNT当量模型和TNO(The Netherlands Organization)多能模型计算不同摩尔百分比的氢氧推进剂混合反应时产生爆炸冲击波的危害性,并模拟爆炸冲击波造成的事故影响范围,然后对两种模型的仿真结果加以对比分析,根据最不利原则选取出最终需要的结果,最后划分出安全设防距离。由仿真结果可知,不同的氢氧混合摩尔百分比造成的爆炸后果不同,同时TNT当量模型在爆炸近场处高估了爆炸超压值,在爆炸远场处低估了爆炸超压值,而TNO多能模型在理论上有效地对这一缺陷进行了弥补。对航天发射场的安全布局起到了一定的参考价值。     

土介质炸点痕迹反演爆源特征的研究

论述了土介质炸点痕迹特征参量及其控制影响因素,通过对国内外现有研究成果的再分析和有效利用,得到爆源特征与炸点痕迹特征参量之间的关系。在此基础上,建立利用土介质上的炸点痕迹反演爆源特征的方法,并实现编程计算。通过实例分析说明了该方法的可行性。     

苯乙烯丙烯酸共聚物/碳黑混合体系粉尘爆炸特性研究

采用MIE-D1.2型最小点火能测试装置及20 L球型粉尘爆炸测试装置,对苯乙烯丙烯酸共聚物/碳黑混合体系粉尘的爆炸特性进行研究。结果表明,过74μm、58μm、47μm孔径筛的粉尘对静电火花敏感,其最小点火能表征值分别为610 mJ、361 mJ、201 mJ。随粉尘质量浓度增加,最小点火能呈现先减小后增加的规律。随粉尘粒径减小,最小点火能与粉尘质量浓度变化关系曲线向低粉尘质量浓度和低点火能量方向偏移,且对应的最敏感爆炸质量浓度从500 g/m~3降至200 g/m~3。随粉尘质量浓度增加,过147μm、74μm、47μm孔径筛的苯乙烯丙烯酸共聚物/碳黑混合体系粉尘爆炸压力及爆炸压力上升速率呈现先增加后减小趋势。在相同粉尘质量浓度下,中位径小于74μm的苯乙烯丙烯酸共聚物/碳黑混合体系粉尘,粉尘的爆炸压力增幅明显减小。苯乙烯丙烯酸共聚物/碳黑混合体系粉尘爆炸下限质量浓度为25 g/m~3,最大爆炸指数为14.636 MPa·m/s,爆炸危险等级划分为St1。