Numerical simulation of a novel blast wave mitigation device

This paper proposes a novel blast wave mitigation device, consisting of a piston–cylinder assembly. A shock wave is induced inside the device when it is subject to a blast wave. The shock wave propagates inside the device and is reflected repeatedly. The physical processes within the blast wave mitigation device are simulated numerically. Numerical predictions are in excellent agreement with analytical solutions for special cases of the investigated problem that are available in the literature. The peak pressure on the base of the device caused by the blast wave is studied using a number of design parameters. The numerical simulation shows that, although the transmitted impulse remains practically unchanged, the peak pressure of the blast wave can be reduced by as much as 98%, or even higher, depending on the design parameters chosen.     

Prediction of blast wave effects on a developed site

The guidelines for protecting against and mitigating explosion hazards require knowledge and either the experimental or theoretical evaluation of blast wave parameters. To this end, this article proposes a numerical method for simulating blast wave propagation in complex geometries. This method permits an on-the-ground TNT-like explosion and the subsequent blast wave to be simulated, with the possibility of modifying the ground topology by adding a number of obstacles. The numerical model is explored from both a qualitative and quantitative point of view by comparing it to experimental data, with a correct determination of the wave amplitude and phase signature at different key-positions. The pressure and Mach number distributions deduced from simulations in complex congested areas highlight various fluid–solid interactions, such as regular reflections and diffractions. The iso-damage diagrams of different obstacle layouts are also presented and discussed.     

Two regimes of blast wave propagation in a bubbly liquid

The evolution of a spherically symmetric layer in a two-phase bubbly liquid under the action of a blast wave is numerically studied. It is demonstrated that such a medium can feature two characteristic regimes of blast wave propagation.     

Pressure wave propagation in a short tunnel caused by passing trains

The aerodynamic behaviour of pressure transients in railway tunnels due to the passage of train are accompanied by the wave reflection between the entrance and exit of the tunnel as well as the tunnel wall and train surface. This study presents a numerical method using the Harten-Yee upwind TVD scheme to simulate compression wave transients and reflection in a railway tunnel. The strength of the compression wave generated by the entrance of a train into a tunnel emerged within the applied calculation procedure, without any special consideration of the boundary condition for the entrance or exit of the tunnel as well as the train nose and tail. Good agreements were obtained when the initial pressure rise at tunnel entrance and the process of pressure variations in a short tunnel, which were calculated by the presented numerical method, were compared with experimental data. The pressure wave reflection at the entrance and exit, the unsteady flow around the train and the transmitted wave at tunnel exit portal were simulated. The transmitted wave at the tunnel exit portal is also discussed. In addition a method to determine the non-reflective boundary condition for multi-dimensional Euler equation is given.

     

Pressure wave propagation in a short tunnel caused by passing trains

The aerodynamic behaviour of pressure transients in railway tunnels due to the passage of train are accompanied by the wave reflection between the entrance and exit of the tunnel as well as the tunnel wall and train surface. This study presents a numerical method using the Harten-Yee upwind TVD scheme to simulate compression wave transients and reflection in a railway tunnel. The strength of the compression wave generated by the entrance of a train into a tunnel emerged within the applied calculation procedure, without any special consideration of the boundary condition for the entrance or exit of the tunnel as well as the train nose and tail. Good agreements were obtained when the initial pressure rise at tunnel entrance and the process of pressure variations in a short tunnel, which were calculated by the presented numerical method, were compared with experimental data. The pressure wave reflection at the entrance and exit, the unsteady flow around the train and the transmitted wave at tunnel exit portal were simulated. The transmitted wave at the tunnel exit portal is also discussed. In addition a method to determine the non-reflective boundary condition for multi-dimensional Euler equation is given.     

坑道内爆炸冲击波传播规律的研究

摘 要：炸药在坑道内爆炸将产生沿坑道传播的空气冲击波,冲击波与坑道壁相互作用,使得坑道内冲击波的传播规律明显不同于自由大气中冲击波的传播规律。使用有限元计算软件LS-DYNA对已有的坑道试验进行了数值模拟,通过数值模拟结果与试验结果数据比较,验证了有限元计算中所使用的计算模型和参数取值的合理性。利用数值模拟的计算结果,使用量纲分析理论,拟合了距离爆炸中心点一定距离处爆炸冲击波超压峰值的计算公式,并与试验结果进行对比,分析了该公式的适用性,为坑道中爆炸冲击波的传播规律研究提供了参考。
     

Calculation of the process of blast wave diffraction in a cylindrical channel

We performed a numerical investigation of the process of transition of a spherical leading front into a plane one in a cylindrical channel. The processes of the collision of reflected shock waves, formation of a nonstationary grating-like structure of flow, and of the overtaking interaction of shock waves are investigated. We found that in the presence of hot gas layers on the walls of the channel a plane head front is not formed.Notation r,z radial and axial coordinates - t time - p pressure - T temperature - u,v radial and axial velocity components - a distance between the planes - b half-distance between the planes - d channel diameter (in experimental data) - R ₀ inital radius of the shock wave front - p ^*(r, z),u ^*(r, z),v ^*(r, z),T ^*(r, z) distributions of pressure, radial, and axial velocity, and temperature at the initial time moment - L length of the zone of formation of a plane front - T ₁,T ₂ temperature of gas in the lower and upper thermal layers - h ₁,h ₂ thicknesses of the lower and upper thermal layers Department of Theoretical Problems of the Russian Academy of Sciences, Moscow Physicotechnical Institute, Moscow. Translated from Inzhernerno-Fizicheskii Zhurnal, Vol. 68, No. 4, pp. 583–588, July–August, 1995.     

Temperature dependence of fracture toughness in cleavage fracture of iron and iron alloys

In order to examine the temperature dependence of fracture toughness in cleavage fracture, exploratory work was carried out. Then effects of alloying elements and micro-structure on the low temperature fracture toughness were studied quantitatively in iron and iron alloys.The results indicate that (1) the relationship between fracture toughness G_ic and testing temperature T at low temperatures is G_ic = G_o exp (T/T_e), where G_e and T_o are the material constants; (2) G_o exhibits a stro dependence on solute carbon and nitrogen contents but is independent of micro-structure and other elements; (3) G_o values increases with increasing solute carbon and nitrogen contents; (4) T_o depends on the structure; (5) $\text{1}\text{T}{}_{\mathrm{o}}$ values increase with increasing nickel and manganese contents, to the contrary, decrease with increasing carbon, silicon and phosphorus contents; and (6) $\text{1}\text{T}{}_{\mathrm{o}}$ values increase with decreasing grain size.     

Ground iron blast furnace slag as a matrix for cellulose-cement materials

The use of ground iron blast furnace slag (BFS) as a low-cost alternative to ordinary Portland cement (OPC) binders in fibre-cement products was examined. Both high quality softwood fibres and residual sisal from agricultural waste were chemically pulped and used as reinforcement. Composites based on several different binder formulations consisting of slag chemically activated by mixtures of gypsum and hydrated lime displayed their optimum strength and fracture toughness properties at fibre contents between 8% and 12%, with values in the ranges of 14.7–24.5 MPa and 1.13–2.36 kJ/m², respectively. Corresponding flexural moduli lay in the range 4.3–7.8 GPa and, at 12% fibre content, the composites possessed water absorption values up to 34% by mass and densities in the region of 1.3 g/cm³. A formulation of BFS activated by 10% gypsum and 2% lime presented a good compromise between strength and energy absorption combined with a reasonable price.     

Brittle fracture of solids caused by intense electron beams

In the mid 1960s, powerful pulse electron beam accelerators having a voltage of some millions of volts were invented and later used to fracture various materials. Experimental data analysis allowed discovery of a new mode of fracture in several ductile crystals caused by a specific energy supply to the crack tip. The mode differs from well known thermomechanical modes of fracture caused by the “ heat-thermostress-crack ” mechanism. This new mode is called the electron fracture mode (EFM). It is characterized by the following three special features, (i) Initial macrocracks in a specimen do not affect the threshold of fracture; that is, the value of the beam intensity at which the specimen breaks, (ii) The fracture of different materials, which can be very ductile at usual mechanical loads, occurs in a brittle manner; that is, the specimen usually splits by a crack without any residual deformation, (iii) The splitting cracks propagate with supersonic velocities. These data are controversial from the point of view of common fracture mechanics and, hence, they cannot be understood or explained from the traditional position.The purpose of the present study is to create a simple practical model of the EFM. The basic viewpoint can be briefly summarized as follows: during irradiation of a solid by a high intensity electron beam, some solid plasma clots are formed and act as “ blades ” or “ wedges,” cutting the crystalline specimen.In the Introduction, experimental data on the EFM are analyzed and discussed, while the peculiarities of the EFM are specified. As a result, it is concluded that the processes caused by the EFM are unusual for the common concepts of fracture mechanics. In Section 2 the invariant Γ-integrals of an electromagnetic deformable medium are modified for supersonic singularities. The basic model and some problems serving to explain and describe the EFM are formulated. In Section 3, the relativistic electron interactions in beams are considered. Using Γ-integrals, we derive the law of the interaction of two moving relativistic charges; that is, the generalized Coulomb's law for relativistic charges. In particular, when two relativistic electrons, e, move with the same velocity, v, one behind the other along a rectilinear trajectory, the force, F, acting upon the rear electron is equal to: where R is the distance between the electrons, c is the speed of light in the vacuum, and a is the phase-speed of light in a medium having electromagnetic constants, μ, , and '. It appears that two electrons moving faster than the phase-speed of light attract one to the other, as distinct from the common Coulomb law. Hence, the beams of such relativistic electrons tend to self-pack and self-compress. The latter problem is studied using a periodic chain model of the electron beam. In Section 4, the dynamic elastic problem of supersonic cutting by a thin wedge is formulated and solved, and the drag force is calculated. In Section 5, the problem of deceleration of the moving wedge is solved in quasi-steady approximation. The length of a resulting cut, that is, the final crack, is determined. Some applications of the analytical solutions are given. In Section 6, the theoretical results are analyzed and compared with experimental results. The role of relativistic electrons is estimated and some parameters of solid-state electron plasma clots are defined. In the Conclusion, the necessity of further study of this mysterious phenomenon is emphasized.     

Predicting fracture of layered composites caused by internal instability

The present paper estimates the critical strain at which ply instability occurs in compressible layered composites under uniaxial or biaxial compression. This is achieved by using two methods: the continuum approach and the model of piecewise-homogeneous medium in conjunction with the three-dimensional stability theory. The accuracy of the continuum theory is examined by taking into account the influence of layer thickness, stiffness properties and biaxiality of loads. Also, upper and lower bounds for the critical buckling strain of laminates with interfacial defects (cracks with connected edges) are determined using the results for perfectly bonded and sliding layers.     

Creep fracture phenomena in iron embrittled by liquid copper

Tensile creep studies of the embrittlement of notched iron by liquid copper in the range 1100 to 1130° C have shown the embrittlement to be of the delayed failure type, from diffusion-controlled, grain boundary penetration by copper with accompanying extended surface notching. Linking of surface cracks and internal cavities along copper-penetrated grain boundaries was the final failure mode. Crack growth was determined by the rate and amount of copper diffusion ahead of a crack and not by the feed of liquid to the crack tip.     

Surface alloying of pure iron by molybdenum and tin by a continuous wave CO2 laser

It is now well established that considerable improvement in the mechanical/chemical properties of near surface regions of materials can be achieved by the process of laser surface alloying. The change in chemistry at the surface is attained through the process of melting and mixing of a predeposited coating and a thin layer of the substrate. Rapid solidification of this molten region at the surface then results in the development of very interesting microstructural features. In the present work, an attempt was made to surface alloy pure iron by molybdenum and/or tin by using a continuous wave CO₂ laser. Morphology of the resulting microstructural features was examined by optical as well as scanning electron microscopy and the compositional details were determined using an electron probe microanalyser. Microhardness measurements were carried out as a function of depth from the laser-treated surface. This paper discusses the results of these investigations and delineates the roles of the various parameters on the chemistry and microstructure of the surface alloys formed as a result of laser treatment.     

温压炸药爆炸冲击波在坑道内的传播规律研究

温压炸药爆炸冲击波对坑道内的人员和设备构成严重威胁。在TNT爆炸试验数据验证的基础上,利用AUTODYN软件建立了炸药堵口爆炸的数值计算模型。基于JWL-Miller能量释放模型计算原理,通过与TNT冲击波的对比,研究了某型温压炸药爆炸冲击波在坑道内的传播特性;通过与空旷地面爆炸冲击波的对比,研究了坑道对温压炸药爆炸冲击波的约束作用。研究结果表明:温压炸药具有更大的破坏威力,温压炸药的超压峰值和正相冲量平均为TNT的1.91倍和1.82倍,其超压和冲量等效TNT当量系数分别为1.94和2.21;坑道对温压炸药冲击波的约束作用明显,其超压峰值和正相冲量平均值分别为空旷地面上的13.55倍和15.21倍。     

Study of an interaction of a blast wave with a destructible screen made of a granular material

The results of an experimental study of three cases of the interaction of a shock wave with screens made of an easily destructible material are given. The dependence of the coefficient of attenuation on the thickness and location of the screen is obtained.     

Deformation and fracture of sintered iron

Conclusion Metallographic-fractographic evaluation of the distribution of microdeformation processes in the volume of porous iron at the static tensile test demonstrates the influence of porosity on the concentration of deformation flows into the microvolume of connections. The degradation of deformation strengthened active volumes with the subsequent degradation of macrodeformation characteristics can occur in geometrically and structurally slightly developed connections.Published in Fiziko-Khimicheskaya Mekhanika Materialov, No. 3, pp. 50–54, May–June, 1992.     

Simulation of contact-induced fracture

The boundary element method is applied to the analysis of crack growth induced by the action of an indenter on monolithic ceramic materials and on thin ceramic coatings. The path followed by the crack has non-trivial geometry. The direction of crack propagation is determined using the maximum mechanical energy release rate criterion. The chief advantage of the method is the ease with which problems having complex geometry or loading (giving rise to complex crack trajectories) can be treated. In this paper we present models of contact damage induced by cylindrical and spherical indenters and show that both the crack geometry and the crack energetics can be accurately represented. In particular, the method allows the study of the early stages of crack growth, which are difficult to observe experimentally.     

刚性地面对爆炸冲击波传播规律的影响

为研究刚性地面对爆炸冲击波传播特性的影响,采用ANSYS/LS-DYNA对TNT炸药在刚性地面工况下起爆后的整个过程进行了数值分析,得到了刚性地面工况下爆炸冲击波的传播规律以及特征参数,并将其与自由空气工况下爆炸冲击波的传播规律以及特征参数进行了对比分析。建立了TNT炸药在刚性地面工况下起爆后,比例距离在0.5 m/kg^(1/3)≤Z≤3 m/kg(1/3)≤Z≤3 m/kg^(1/3)之间时的冲击波超压峰值简便计算公式,分析了起爆高度对刚性地面附近冲击波的影响规律,并得到了刚性地面反射增强效应作用范围的判别公式。数值分析所得结果可为结构或构件在刚性地面上发生爆炸时荷载的初步确定提供一定的参考。