Crack propagation and damage evolution of metallic cylindrical shells under internal explosive loading

This paper investigates the three-dimensional crack propagation and damage evolution process of metallic column shells under internal explosive loading.The calibration of four typical failure parameters for 40CrMnSiB steel was conducted through experiments and subsequently applied to simulations.The numerical simulation results employing the four failure criteria were compared with the differences and similarities observed in freeze-recovery tests and ultra-high-speed tests.This analysis addressed the critical issue of determining failure criteria for the fracture of a metal shell under internal explosive loads.Building upon this foundation,the damage parameter D_(c),linked to the cumulative crack density,was defined based on the evolution characteristics of a substantial number of cracks.The relationship between the damage parameter and crack velocity over time was established,and the influence of the internal central pressure on the damage parameter and crack velocity was investigated.Variations in the fracture modes were found under different failure criteria,with the principal strain failure criterion proving to be the most effective for simulating 3D crack propagation in a pure shear fracture mode.Through statistical analysis of the shell penetration fracture radius data,it was determined that the fracture radius remained essentially constant during the crack evolution process and could be considered a constant.The propagation velocity of axial cracks ranged between 5300 m/s and 12600 m/s,surpassing the Rayleigh wave velocity of the shell material and decreasing linearly with time.The increase in shell damage exhibited an initial rapid phase,followed by deceleration,demonstrating accelerated damage during the propagation stage of the blast wave and decelerated damage after the arrival of the rarefaction wave.This study provides an effective approach for investigating crack propagation and damage evolution.The derived crack propagation and damage evolution law serves as a valuable reference for the development of crack velocity theory and the construction of shell damage evolution modes.

Numerical Simulation of the Flat Ribbon Wound Explosion Containment Vessels Subjected to Internal Blast Loading

In order to constitute engineering design methods of the flat ribbon wound explosion containment vessels, the dynamic response of such vessels subjected to internal explosion loading is simulated using LS-DYNA3D. Three winding angles, 10°, 15°and 20°, are considered. It is shown that among ribbon vessels investigated, the center displacement of outermost ribbons of the vessel with 10°winding angle is the smallest under the same blast loading. The response of vessels loaded in inner core is local. From the center of the cylindrical shell to the bottom cover, the maximum strain gradually decreases. The ribbons are subjected to tension in the length direction and compression in the width direction. Blasting shock energy concentrates on where is close to center section of blasting. For comparison, numerical simulation of a monobloc thick-walled explosion containment vessel is also investigated. It can be found that the biggest deformation of the flat ribbon wound explosion containment vessels is bigger than that of the monobloc thick-walled explosion containment vessel in the center section of blasting under the same TNT. Numerical results are approximately in agreement with experimental ones. It is proved that the ribbon vessels have the valuable properties of ' leak before burst at worst' compared with the monobloc vessels through numerical simulation.

Failure mechanisms and destruction characteristics of cemented coal gangue backfill under compression effect of non-uniform load

Backfill mining is one of the most important technical means for controlling strata movement and reducing surface subsidence and environmental damage during exploitation of underground coal resources. Ensuring the stability of the backfill bodies is the primary prerequisite for maintaining the safety of the backfilling working face, and the loading characteristics of backfill are closely related to the deformation and subsidence of the roof. Elastic thin plate model was used to explore the non-uniform subsidence law of the roof, and then the non-uniform distribution characteristics of backfill bodies’ load were revealed. Through a self-developed non-uniform loading device combined with acoustic emission (AE) and digital image correlation (DIC) monitoring technology, the synergistic dynamic evolution law of the bearing capacity, apparent crack, and internal fracture of cemented coal gangue backfills (CCGBs) under loads with different degrees of non-uniformity was deeply explored. The results showed that: 1) The uniaxial compressive strength (UCS) of CCGB increased and then decreased with an increase in the degree of non-uniformity of load (DNL). About 40% of DNL was the inflection point of DNL-UCS curve and when DNL exceeded 40%, the strength decreased in a cliff-like manner;2) A positive correlation was observed between the AE ringing count and UCS during the loading process of the specimen, which was manifested by a higher AE ringing count of the high-strength specimen. 3) Shear cracks gradually increased and failure mode of specimens gradually changed from “X” type dominated by tension cracks to inverted “Y” type dominated by shear cracks with an increase in DNL, and the crack opening displacement at the peak stress decreased and then increased. The crack opening displacement at 40% of the DNL was the smallest. This was consistent with the judgment of crack size based on the AE b-value, i. e., it showed the typical characteristics of “small b-value-large crack and large b-value-small crack”. The research results are of significance for preventing the instability and failure of backfill.

Study on energy release characteristics of reactive material casings under explosive loading

Reactive Materials(RMs),a new material with structural and energy release characteristics under shockinduced chemical reactions,are promising in extensive applications in national defense and military fields.They can increase the lethality of warheads due to their dual functionality.This paper focuses on the energy release characteristics of RM casings prepared by alloy melting and casting process under explosive loading.Explosion experiments of RM and conventional 2A12 aluminum alloy casings were conducted in free field to capture the explosive fireballs,temperature distribution,peak overpressure of the air shock wave and the fracture morphology of fragments of reactive material(RM)warhead casings by using high-speed camera,infrared thermal imager temperature and peak overpressure testing and scanning electron microscope.Results showed that an increase of both the fireball temperature and air shock wave were observed in all RM casings compared to conventional 2A12 aluminum ally casings.The RM casings can improve the peak overpressure of the air shock wave under explosion loading,though the results are different with different charge ratios.According to the energy release characteristics of the RM,increasing the thickness of RM casings will increase the peak overpressure of the near-field air shock wave,while reducing the thickness will increase the peak overpressure of the far-field air shock wave.

Explosion resistance performance of reinforced concrete box girder coated with polyurea:Model test and numerical simulation

To study the anti-explosion protection effect of polyurea coating on reinforced concrete box girder,two segmental girder specimens were made at a scale of 1:3,numbered as G(without polyurea coating)and PCG(with polyurea coating).The failure characteristics and dynamic responses of the specimens were compared through conducting explosion tests.The reliability of the numerical simulation using LS-DYNA software was verified by the test results.The effects of different scaled distances,reinforcement ratios,concrete strengths,coating thicknesses and ranges of polyurea were studied.The results show that the polyurea coating can effectively enhance the anti-explosion performance of the girder.The top plate of middle chamber in specimen G forms an elliptical penetrating hole,while that in specimen PCG only shows a very slight local dent.The peak vertical displacement and residual displacement of PCG decrease by 74.8% and 73.7%,respectively,compared with those of specimen G.For the TNT explosion with small equivalent,the polyurea coating has a more significant protective effect on reducing the size of fracture.With the increase of TNT equivalent,the protective effect of polyurea on reducing girder displacement becomes more significant.The optimal reinforcement ratio,concrete strength,thickness and range of polyurea coating were also drawn.

Fragment spatial distribution of prismatic casing under internal explosive loading

Non-cylindrical casings filled with explosives have undergone rapid development in warhead design and explosion control.The fragment spatial distribution of prismatic casings is more complex than that of traditional cylindrical casings.In this study,numerical and experimental investigations into the fragment spatial distribution of a prismatic casing were conducted.A new numerical method,which adds the Lagrangian marker points to the Eulerian grid,was proposed to track the multi-material interfaces and material dynamic fractures.Physical quantity mappings between the Lagrangian marker points and Eulerian grid were achieved by their topological relationship.Thereafter,the fragment spatial distributions of the prismatic casing with different fragment sizes,fragment shapes,and casing geometries were obtained using the numerical method.Moreover,fragment spatial distribution experiments were conducted on the prismatic casing with different fragment sizes and shapes,and the experimental data were compared with the numerical results.The effects of the fragment and casing geometry on the fragment spatial distributions were determined by analyzing the numerical results and experimental data.Finally,a formula including the casing geometry parameters was fitted to predict the fragment spatial distribution of the prismatic casing under internal explosive loading.

Micro defects formation and dynamic response analysis of steel plate of quasi-cracking area subjected to explosive load

As the protective component,steel plate had attracted extensive attention because of frequently threats of explosive loads.In this paper,the evolution of microstructure and the mechanism of damage in the quasi-cracking area of steel plate subjected to explosive load were discussed and the relationships between micro defects and dynamic mechanical response were revealed.After the explosion experiment,five observation points were selected equidistant from the quasi-cracking area of the section of the steel plate along the thickness direction,and the characteristics of micro defects at the observation points were analyzed by optical microscope(OM),scanning electron microscope(SEM) and electron backscattered diffraction(EBSD).The observation result shows that many slip bands(SBs) appeared,and the grain orientation changed obviously in the steel plate,the two were the main damage types of micro defects.In addition,cracks,peeling pits,grooves and other lager micro defects were appeared in the lower area of the plate.The stress parameters of the observation points were obtained through an effective numerical model.The mechanism of damage generation and crack propagation in the quasicracking area were clarified by comparing the specific impulse of each observation point with the corresponding micro defects.The result shows that the generation and expansion of micro defects are related to the stress area(i.e.the upper compression area,the neutral plane area,and the lower tension area).The micro defects gather and expand at the grain boundary,and will become macroscopic damage under the continuous action of tensile stress.Besides,the micro defects at the midpoint of the section of the steel plate in the direction away from the explosion center(i.e.the horizontal direction) were also studied.It was found that the specific impulse at these positions were much smaller than that in the thickness direction,the micro defects were only SBs and a few micro cracks,and the those decreased with the increase of the distance from the explosion center.

Theoretical and experimental analyses of casing collapsing strength under non-uniform loading

Failure data from oilfield showed that casings which were designed according to API standards were deformed and collapsed in salt formations. The main reason for decrease in strength may be caused by non-uniform loading(NUL) that was not considered in traditional casing collapsing strength design or that the designing method should be improved and developed. Obviously, the calculation of casing collapse strength is one of the key factors in casing design. However, the effect of NUL on casing collapse strength was generally neglected in the present computational methods. Therefore, a mechanical model which can calculate casing collapse strength under NUL was established based on the curved beam theory of the elasticity and was solved using displacement method. Simultaneously, three anti-collapse experiments were performed on C110 casing under NUL, and the strain and deformation laws of three casings in the process of collapse were obtained by the electrical method. Yield limit of every casing was obtained by analyzing those data. Experimental results are consistent with the results of calculation of new model. It indicates that the model can be used to calculate yield limit loading of casings under NUL.

Microscopic defects formation and dynamic mechanical response analysis of Q345 steel plate subjected to explosive load

As the basic protective element, steel plate had attracted world-wide attention because of frequent threats of explosive loads. This paper reports the relationships between microscopic defects of Q345 steel plate under the explosive load and its macroscopic dynamics simulation. Firstly, the defect characteristics of the steel plate were investigated by stereoscopic microscope(SM) and scanning electron microscope(SEM). At the macroscopic level, the defect was the formation of cave which was concentrated in the range of 0-3.0 cm from the explosion center, while at the microscopic level, the cavity and void formation were the typical damage characteristics. It also explains that the difference in defect morphology at different positions was the combining results of high temperature and high pressure. Secondly, the variation rules of mechanical properties of steel plate under explosive load were studied. The Arbitrary Lagrange-Euler(ALE) algorithm and multi-material fluid-structure coupling method were used to simulate the explosion process of steel plate. The accuracy of the method was verified by comparing the deformation of the simulation results with the experimental results, the pressure and stress at different positions on the surface of the steel plate were obtained. The simulation results indicated that the critical pressure causing the plate defects may be approximately 2.01 GPa. On this basis, it was found that the variation rules of surface pressure and microscopic defect area of the Q345 steel plate were strikingly similar, and the corresponding mathematical relationship between them was established. Compared with Monomolecular growth fitting models(MGFM) and Logistic fitting models(LFM), the relationship can be better expressed by cubic polynomial fitting model(CPFM). This paper illustrated that the explosive defect characteristics of metal plate at the microscopic level can be explored by analyzing its macroscopic dynamic mechanical response.

Study on dynamic response of multi-degree-of-freedom explosion vessel system under impact load

In order to study the dynamic response and calculate the axial dynamic coefficient of the monolayer cylindrical explosion vessel,the wall of vessel is simplified as a multi-degree-of-freedom(MDoF) undamped elastic foundation beam.Decoupling the coupled motion equation and using Duhamel's integrals,the solutions in generalized coordinates of the equations under exponentially decaying loads,square wave loads and triangular wave loads are calculated.These solutions are consistent in form with the solutions of single-degree-of-freedom(SDoF) undamped forced vibration simplified model.Based on the model,equivalent MDoF design method(also called MDoF dynamic coefficient method) of cylindrical explosion vessel is proposed.The traditional method can only predict the dynamic coefficient of torus portion around the explosion center,but this method can predict that of the vessel wall at any axial n dividing point position.It is verified that the prediction accuracy of this model is greatly improved compared with the SDoF model by comparing the results of this model with SDoF model and numerical simulation in different working conditions.However,the prediction accuracy decreases as the scaled distance decreases and approaches the end of the vessel,which is related to the accuracy of the empirical formula of the implosion load,the simplification of the explosion load direction,the boundary conditions,and the loading time difference.

Generation of electromotive force in igneous rocks subjected to non-uniform loading

When one end of an air-dry igneous rock block was uniaxially loaded in laboratory, there appeared an electromotive force that made electric currents flow from the stressed volume to the unstressed volume. Quartz-free rocks such as gabbro also generated this force, stronger than quartz-bearing rocks such as granite. This indicates that the piezoelectric effect of quartz and the electrokinetic effect of pore water do not make a large contribution toward generating the electromotive force. We focus on peroxy bond that is one of the abundant lattice defects in igneous rock-forming minerals. When mechanical loading deforms the lattice structure around this defect and breaks its bond, its energy levels change and act like an accepter. As an electron is trapped at this defect from a neighbor 02- site, a positive hole is activated there. They attempt to diffuse toward the unstressed volume through the valence band and are simultaneously affected by the attractive electric force with the electrons trapped in peroxy bonds. This leads to a polarization in the stressed volume and the generation of electromotive force between the stressed and unstressed volumes. Similar electromotive force may be generated in the Earth＇s crust where inhomogeneous stress/strain is changing.

Numerical Simulation of Concrete Plate Damaged Under Explosive Loading

The grisliness after-effects can be induced by explosion accident with the collapsing of the structures, the demolishing of the equipments and the casualty of the human beings. Isolation belt constructed between the blast point and the construction is one of the useful design schemes for blast resistance. The nonlinear procedure ANSYS/LSDYNA970 is used to simulate the contact detonation and the isolation belt of blast resistance filled with the air or water respectively. The results indicate that the maximal damage can be caused by the contact detonation, and the isolation belt of blast resistent filled with water can reduce the damage greatly.

Strength and Toughness of Steel Fibre Reinforced Reactive Powder Concrete Under Blast Loading

The blast resistance of structures used in buildings needs to be investigated due to the increased threat of a terrorist attack. The damage done by Composition B or Powergel to steel fibre reinforced reactive powder concrete （SFRPC） panels and ordinary reinforced concrete （RC） panels of equivalent static flexural strength is compared. A 0. 5 kg charge was detonated at a distance of 0. 1 m from the 1. 3 m × 1. 0 m × 0. 1 m （thick） panels, which were simply supported and spaning 1.3 m. Dynamic displacement measurements, high-speed video recording and visual examination of the panels for spall and breach were undertaken. The SFRPC panels withstood the bare charge blast better than the reinforced ordinary concrete panels. Neither type of panel was breached using a O. 5 kg charge, The RC panel exhibited more spalling when Composition B was used. Under successive Composition B loading conditions, the RC panel was breached. In comparison the SFRPC panel was not breached. Exposure to fragmenting charge loading conditions confirmed these performance differences between the SFRPC panel and the reinforced ordinary concrete panel.

Experimental investigation on dynamic response and damage models of circular RC columns subjected to underwater explosions

Reinforced concrete(RC) columns are widely used as supporting structures for high-piled wharfs.The study of damage model of a RC column due to underwater explosion is a critical issue to assess the wharfs antiknock security.In this study,the dynamic response and damage model of circular RC columns subjected to underwater explosions were investigated by means of scaled-down experiment models.Experiments were carried out in a 10.0 m diameter tank with the water depth of 2.25 m,under different explosive quantities(0.025 kg-1.6 kg),stand-off distances(0.0 m-7.0 m),and detonation depths(0.25 m-2.0 m).The shock wave load and dynamic response of experiment models were measured by configuring sensors of pressure,acceleration,strain,and displacement.Then,the load distribution characteristics,time history of test data,and damage models related to present conditions were obtained and discussed.Three damage models,including bending failure,bending-shear failure and punching failure,were identified.In addition,the experie nce model of shock wave loads on the surface of a RC column was proposed for engineering application.

Response of Box-Type Structures Under Internal-Blast Loading

The tests of box-type structures under internal-blast loading are carried out. Then a numerical analysis of the test structures is done using a fully coupled numerical finite element model. The break-up process of the structure is simulated. The failure modes of the simulated structure agree well with the experimental results. The effects of the size of the reinforcing bars and the detailing of connections among the rebars in the concrete on the throw velocity of the fragments are discussed.

Radial explosion strain and its fracture effect from confined explosion with charge of cyclonite

Instrumented experiments were conducted in concrete models to study the explosion-induced radial strain and fracture effect of rock-like media under confined explosion with a charge of cyclonite. As a charge was exploded, two different radial strain waves were sequentially recorded by a strain gage at a distance of 80 mm from the center of charge. Through the attenuation formula of the maximum compressive strain(εrmax), the distribution of εrmax and its strain rate（ ） between the charge and gage were obtained. The effect of the two waves propagating outwards on the radial fracture of surrounding media was discussed. The results show that the two waves are pertinent to the loading of shock energy （Es） and bubble energy （Eb） against concrete surrounding charge, respectively. The former wave lasts for much shorter time than the latter. The peak values of εrmax and of the former are higher than those of the latter, respectively.

Non-capacity transport of non-uniform bed load sediment in alluvial rivers

Rivers often witness non-uniform bed load sedim ent transport. For a long tim e, non-uniform bed load transport has been assum ed to be at capacity regime determined exclusively by local flow. Yet whether the capacity assumption for non-uniform bed load transport is justified remains poorly understood. Here, the relative time scale of non-uniform bed load transport is evaluated and non-capacity and capacity models are compared for both aggradation and degradation cases with observed data. As characterized by its relative time scale, the adaptation of non-uniform bed load to capacity regime should be fulfilled quickly. However, changes in the flow and sedim ent inputs from upstream or tributaries hinder the adaptation. Also, the adaptation to capacity regime is size dependent, the finer the sediment size the slower the adaptation is, and vice versa. It is shown that the capacity model may entail considerable errors compared to the non-capacity model. For modelling of non-uniform bed load, non-capacity modelling is recommended, in which the temporal and spatial scales required for adaptation are explicitly appreciated.

A laboratory method to simulate seismic waves induced by underground explosions

The seismic waves induced by underground explosions generate geological hazards affecting deep buried tunnels such as rockbursts and engineering-induced earthquakes. This issue is difficult to study through full-scale testing due to the expense and unpredictable danger. To solve this problem, the authors developed experimental apparatus and presented a laboratory method to simulate seismic waves induced by underground explosions. In this apparatus, a combined structure of a diffusive-shaped water capsule and a special-shaped oil capsule was designed. This structure can provide an applied confining stress and freely transmit the stress wave generated by external impact. Therefore, the coupled loading of in situ stress and seismic waves induced by underground explosions in the deep rock mass was simulated. The positive pressure time and peak value of the stress wave could be adjusted by changing the pulse-shaper and the initial impact energy. The obtained stress waves in the experiments correspond to that generated by 0.15-120 kt of TNT equivalent explosion at a scaled distance of 89.9-207.44 m/kt.

Finding the buckling load of non-uniform columns using the iteration perturbation method

The aim of this study is to calculate the critical load of variable inertia columns. The example studied in this paper can be used as a paradigm for other non-uniform columns. The wavelength of equivalent vibratory system is used to calculate the critical load of the trigonometrically varied inertia column. In doing so, the equilibrium equation of the column is theoretically studied using the perturbation method. Accuracy of the calculated results is evaluated by comparing the solution with numerical results. Effect of improving the initial guess on the solution accuracy is investigated. Effects of varying parameters of the trigonometrically varied inertia and the uniformly tapered columns on their stability behavior are studied. Finally, using the so-called ＂perfectibility＂ parameter, two design goals, i.e., being lightweight and being strong, are studied for the discussed columns.

EFFECT OF ORIGINAL EXPLOSION STRESS WAVE ON PROPAGTION OF EXPL

ＴｈｅｄｅｔｏｎａｔｉｏｎｂｅｈａｖｉｏｕｒｓｏｆｔｈｅｃｙｌｉｎｄｒｉｃｃａｒｔｒｉｄｇｅｓｗｉｔｈａｎｄｗｉｔｈｏｕｔａｎａｘｉａｌｃａｖｉｔｙａｎｄｍａｄｅｏｆＨｅｘｏｇｅｎ（ＲＤＸ）ｉｎｌｏｗｅｒｄｅｎｓｉｔｙ，ａｎｄｔｈｅｐａｒａｍｅｔｅｒｓｏｆｅｘｐｌｏｓｉｏｎｓｔｒｅｓｓｗａｖｅｓｇｅｎｅｒａｔｅｄｉｎｃｏｎｃｒｅｔｅｍｏｄｅｌｓｆｒｏｍｔｈｅｅｘｐｌｏｓｉｏｎｏｆｔｈｏｓｅｃａｒｔｒｉｄｇｅｓｗｅｒｅａｎａ－ｌｙｚｅｄ．Ｆｒａｇｍｅｎｔａｔｉｏｎｅｘｐｅｒｉｍｅｎｔｓ，ｕｓｉｎｇｃａｒｔｒｉｄｇｅｓｗｉｔｈａｎｄｗｉｔｈｏｕｔｃａｖｉｔｙｗｅｒｅｃｏｎｄｕｃｔｅｄｉｎｐｌｅｘｉｇｌａｓｓｍｏｄｅｌｓｉｎｏｒｄｅｒｔｏｅｘ－ｐｌｏｒｅｔｈｅｉｎｆｌｕｅｎｃｅｏｆｖａｒｉａｔｉｏｎｓｉｎｔｈｅｐａｒａｍｅｔｅｒｓｏｆｏｒｉｇｉｎａｌｅｘｐｌｏｓｉｏｎｓｔｒｅｓｓｗａｖｅ，ｄｕｅｔｏｃｈａｎｇｅｉｎｔｈｅｒａｔｉｏｏｆｓｈｏｃｋｅｎｅｒｇｙ（Ｅｓ）ｔｏｂｕｂｂｌｅｅｎｅｒｇｙ（Ｅｂ）ｕｎｄｅｒｃｏｕｐｌｅｃｈａｒｇｅａｓｔｏｔａｌｅｘｐｌｏｓｉｏｎ（Ｅ）ｒｅｌｅａｓｅｄｂｙｔｈｅｃａｒｔｒｉｄｇｅ?