The graded density impactor(GDI)dynamic loading technique is crucial for acquiring the dynamic physical property parameters of materials used in weapons.The accuracy and timeliness of GDI structural design are key to ...The graded density impactor(GDI)dynamic loading technique is crucial for acquiring the dynamic physical property parameters of materials used in weapons.The accuracy and timeliness of GDI structural design are key to achieving controllable stress-strain rate loading.In this study,we have,for the first time,combined one-dimensional fluid computational software with machine learning methods.We first elucidated the mechanisms by which GDI structures control stress and strain rates.Subsequently,we constructed a machine learning model to create a structure-property response surface.The results show that altering the loading velocity and interlayer thickness has a pronounced regulatory effect on stress and strain rates.In contrast,the impedance distribution index and target thickness have less significant effects on stress regulation,although there is a matching relationship between target thickness and interlayer thickness.Compared with traditional design methods,the machine learning approach offers a10^(4)—10^(5)times increase in efficiency and the potential to achieve a global optimum,holding promise for guiding the design of GDI.展开更多
The internal and external flow fields during vented explosions of methane were characterized through numerical simulation,and the capability of numerical simulation thereof was validated by previous experimental data ...The internal and external flow fields during vented explosions of methane were characterized through numerical simulation,and the capability of numerical simulation thereof was validated by previous experimental data at three ignition positions.The venting mechanism was revealed by the simulated concentration distribution,temperature profile,and airflow velocity.The results show rear ignition results in the external methane mass distribution taking the form of"mushroom"and columnar flames in the external space,which can be expressed as a third-order polynomial relationship with distance;central ignition forms a relationship of the form y=AxB.Front ignition causes the temperature to show a tendency to repeated oscillations(rising,falling,and rising).Central ignition generates the maximum vented airflow velocity(V_(max)=320 m/s)upon vent opening.The results indicate that it is acceptable to apply numerical simulation of methane explosions in practice.展开更多
The detonation of fuel-rich explosives yields combustible products that persistently burn upon mixing with ambient oxygen,releasing additional energy through a phenomenon known as the afterburning effect.This process ...The detonation of fuel-rich explosives yields combustible products that persistently burn upon mixing with ambient oxygen,releasing additional energy through a phenomenon known as the afterburning effect.This process greatly influences the evolution of confined blast loading and the subsequent structural response,which is crucial in confined blast scenarios.Given the complex nature of the reaction process,accurate analysis of the afterburning effect remains challenging.Previous studies have either overlooked the mechanisms of detonation product combustion or failed to provide experimental validation.This study introduces a three-dimensional model to effectively characterize the combustion of detonation products.The model integrates chemical reaction source terms into the governing equations to consider the combustion processes.Numerical simulations and experimental tests were conducted to analyze the combustion and energy release from the detonation products of fuel-rich explosives in confined spaces.Approximately 50%of the energy was released during the combustion of detonation products in a confined TNT explosion.Although the combustion of these products was much slower than the detonation process,it aligned with the dynamic response of the structure,which enhanced the explosive yield.Excluding afterburning from the analysis reduced the center-point deformation of the structure by 30%.Following the inclusion of afterburning,the simulated quasistatic pressure increased by approximately 45%.Subsequent comparisons highlighted the merits of the proposed approach over conventional methods.This approach eliminates the reliance on empirical parameters,such as the amount and rate of energy release during afterburning,thereby laying the foundation for understanding load evolution in more complex environments,such as ships,buildings,and underground tunnels.展开更多
To address the problems of low accuracy by the CONWEP model and poor efficiency by the Coupled Eulerian-Lagrangian(CEL)method in predicting close-range air blast loads of cylindrical charges,a neural network-based sim...To address the problems of low accuracy by the CONWEP model and poor efficiency by the Coupled Eulerian-Lagrangian(CEL)method in predicting close-range air blast loads of cylindrical charges,a neural network-based simulation(NNS)method with higher accuracy and better efficiency was proposed.The NNS method consisted of three main steps.First,the parameters of blast loads,including the peak pressures and impulses of cylindrical charges with different aspect ratios(L/D)at different stand-off distances and incident angles were obtained by two-dimensional numerical simulations.Subsequently,incident shape factors of cylindrical charges with arbitrary aspect ratios were predicted by a neural network.Finally,reflected shape factors were derived and implemented into the subroutine of the ABAQUS code to modify the CONWEP model,including modifications of impulse and overpressure.The reliability of the proposed NNS method was verified by related experimental results.Remarkable accuracy improvement was acquired by the proposed NNS method compared with the unmodified CONWEP model.Moreover,huge efficiency superiority was obtained by the proposed NNS method compared with the CEL method.The proposed NNS method showed good accuracy when the scaled distance was greater than 0.2 m/kg^(1/3).It should be noted that there is no need to generate a new dataset again since the blast loads satisfy the similarity law,and the proposed NNS method can be directly used to simulate the blast loads generated by different cylindrical charges.The proposed NNS method with high efficiency and accuracy can be used as an effective method to analyze the dynamic response of structures under blast loads,and it has significant application prospects in designing protective structures.展开更多
With different structural forms of ventilation pipes have various attenuation effects on incident shock waves while meeting ventilation requirements.The attenuation mechanism and the propagation law of shock waves in ...With different structural forms of ventilation pipes have various attenuation effects on incident shock waves while meeting ventilation requirements.The attenuation mechanism and the propagation law of shock waves in ventilation pipes of different structures are investigated by experiments and numerical simulations.Furthermore,for the same structure,the effects of peak pressure and positive pressure time on the attenuation rate are discussed.It is found that the attenuation rate increases with the incident shock wave pressure,and the shock wave attenuation rate tends to reach its limiting value k for the same structure and reasonably short positive pressure time.Under the same conditions,the attenuation rate is calculated using the pressure of the shock wave as follows:diffusion chamber pipe,branch pipe and selfconsumption pipe;the attenuation rate per unit volume is calculated as follows:self-consumption pipe,branch pipe and diffusion chamber pipe.In addition,an easy method is provided to calculate the attenuation rate of the shock wave in single and multi-stage ventilation pipes.Corresponding parameters are provided for various structures,and the margin of error between the formulae and experimental results is within 10%,which is significant for engineering applications.展开更多
Most of the existing studies on tunnel blast wave are based on spherical or grouped charges, however,conventional weapons are mostly cylindrical rather than spherical. In order to analyze the impact of cylindrical cha...Most of the existing studies on tunnel blast wave are based on spherical or grouped charges, however,conventional weapons are mostly cylindrical rather than spherical. In order to analyze the impact of cylindrical charges on the tunnel blast wave loads and to develop a quantitative calculation method, this study carried out experimental and numerical research. Initially, external explosion experiments were conducted using both 35 kg spherical charges and cylindrical charges with aspect ratio of 4.8 at two different distances from the tunnel entrance. Comparative analysis of the blast wave parameters in the tunnel revealed that the explosive equivalent of the cylindrical charges was significantly higher than that of the spherical charges. To address this, an equivalent coefficient κ based on the spherical charges was proposed for the cylindrical charges. Subsequently, numerical simulations were conducted for the experimental conditions, and the numerical simulation results match the experiments well. Through numerical calculations, the reliability of the equivalent coefficient κ under the experimental conditions was verified, and comparison analysis indicated that the explosion energy of cylindrical charges spreads more radially, resulting in more explosion energy entering the tunnel, which is the fundamental reason for the increase in tunnel blast wave loads. Additionally, analyzing the explosion energy ratio entering the tunnel is an effective method for calculating the equivalent coefficient κ. Finally, through more than one hundred sets of numerical calculation results, the impact of the proportional distance λ and the ratio of charge mass to the tunnel cross-section dimension φ on the equivalence coefficients κ was investigated. An empirical formula for the equivalence coefficient κ was derived through fitting, and the accuracy of the formula was validated through literature experimental results. The research findings of this paper will provide valuable guidance for the calculation of blast wave loads in tunnel.展开更多
In the context of deep rock engineering,the in-situ stress state is of major importance as it plays an important role in rock dynamic response behavior.Thus,stress initialization becomes crucial and is the first step ...In the context of deep rock engineering,the in-situ stress state is of major importance as it plays an important role in rock dynamic response behavior.Thus,stress initialization becomes crucial and is the first step for the dynamic response simulation of rock mass in a high in-situ stress field.In this paper,stress initialization methods,including their principles and operating procedures for reproducing steady in-situ stress state in LS-DYNA,are first introduced.Then the most popular four methods,i.e.,explicit dynamic relaxation(DR)method,implicit-explicit sequence method,Dynain file method and quasi-static method,are exemplified through a case analysis by using the RHT and plastic hardening rock material models to simulate rock blasting under in-situ stress condition.Based on the simulations,it is concluded that the stress initialization results obtained by implicit-explicit sequence method and dynain file method are closely related to the rock material model,and the explicit DR method has an obvious advantage in solution time when compared to other methods.Besides that,it is recommended to adopt two separate analyses for the whole numerical simulation of rock mass under the combined action of in-situ stress and dynamic disturbance.展开更多
The numerical simulation of a blast wave of a multilayer composite charge is investigated.A calculation model of the near-field explosion and far-field propagation of the shock wave of a composite charge is establishe...The numerical simulation of a blast wave of a multilayer composite charge is investigated.A calculation model of the near-field explosion and far-field propagation of the shock wave of a composite charge is established using the AUTODYN finite element program.Results of the near-field and far-field calculations of the shock wave respectively converge at cell sizes of 0.25-0.5 cm and 1-3 cm.The Euler--fluxcorrected transport solver is found to be suitable for the far-field calculation after mapping.A numerical simulation is conducted to study the formation,propagation,and interaction of the shock wave of the composite charge for different initiation modes.It is found that the initiation mode obviously affects the shock-wave waveform and pressure distribution of the composite charge.Additionally,it is found that the area of the overpressure distribution is greatest for internal and external simultaneous initiation,and the peak pressure of the shock wave exponentially decays,fitting the calculation formula of the peak overpressure attenuation under different initiation modes,which is obtained and verified by experiment.The difference between numerical and experimental results is less than 10%,and the peak overpressure of both internal and external initiation is 56.12% higher than that of central single-point initiation.展开更多
The design of solid armature of railgun should take full account of its operating conditions and material properties because the armature is subjected to dynamic loading conditions and experiences a complicated electr...The design of solid armature of railgun should take full account of its operating conditions and material properties because the armature is subjected to dynamic loading conditions and experiences a complicated electrical,thermal and mechanical process in the interior ballistic cycle.In this paper present,we first introduced a multi-physical field model of railgun,followed by several examples to investigate the launching process.Especially,we used the explicit finite element method,in which material nonlinearity and geometric nonlinearity were accounted,to investigate the deform behaviors of solid armature.The results show that the dynamic mechanical process of armature is dependent on the armature geometry,material and exciting electric current.By the numerical simulation,the understanding of the fracture mechanism of solid armature was deepened.展开更多
In actual production,deep coal mine roadways are often under typical static-dynamic coupling stress(SDCS)conditions with high ground stress and strong dynamic disturbances.With the increasing number of disasters and a...In actual production,deep coal mine roadways are often under typical static-dynamic coupling stress(SDCS)conditions with high ground stress and strong dynamic disturbances.With the increasing number of disasters and accidents induced by SDCS conditions,the safe and efficient production of coal mines is seriously threatened.Therefore,it is of great practical significance to study the deformation and failure characteristics of the roadway surrounding rock under SDCS.In this paper,the effects of different in-situ stress fields and dynamic load conditions on the surrounding rock are studied by numerical simulations,and the deformation and failure characteristics are obtained.According to the simulation results,the horizontal stress,vertical stress and dynamic disturbance have a positive correlation with the plastic failure of the surrounding rock.Among these factors,the influence of the dynamic disturbance is the most substantial.Under the same stress conditions,the extents of deformation and plastic failure of the roof and ribs are always greater than those of the floor.The effect of horizontal stresses on the roadway deformation is more notable than that of vertical stresses.The results indicate that for the roadway under high-stress conditions,the in-situ stress test must be strengthened first.After determining the magnitude of the in-situ stress,the location of the roadway should be reasonably arranged in the design to optimize the mining sequence.For roadways that are strongly disturbed by dynamic loads,rock supports(rebar/cable bolts,steel set etc.)that are capable of maintaining their effectiveness without failure after certain dynamic loads are required.The results of this study contribute to understanding the characteristics of the roadway deformation and failure under SDCS,and can be used to provide a basis for the support design and optimization under similar geological and geotechnical circumstances.展开更多
Based on RNG k-ε turbulence model and sliding grid technique, solid-liquid two-phase three-dimensional(3-D) unsteady turbulence of full passage in slurry pump was simulated by means of Fluent software. The effects of...Based on RNG k-ε turbulence model and sliding grid technique, solid-liquid two-phase three-dimensional(3-D) unsteady turbulence of full passage in slurry pump was simulated by means of Fluent software. The effects of unsteady flow characteristics on solid-liquid two-phase flow and pump performance were researched under design condition. The results show that clocking effect has a significant influence on the flow in pump, and the fluctuation of flow velocity and pressure is obvious, particularly near the volute tongue, at the position of small sections of volute and within diffuser. Clocking effect has a more influence on liquid-phase than on solid-phase, and the wake-jet structure of relative velocity of solid-phase is less obvious than liquid-phase near the volute tongue and the impeller passage outlet. The fluctuation of relative velocity of solid-phase flow is 7.6% smaller than liquid-phase flow at the impeller outlet on circular path. Head and radial forces of the impeller are 8.1% and 85.7% of fluctuation, respectively. The results provide a theoretical basis for further research for turbulence, improving efficient, reducing the hydraulic losses and wear. Finally, field tests were carried out to verify the operation and wear of slurry pump.展开更多
In order to reveal the dynamic process of hard-thick roof inducing rock burst, one of the most common and strongest dynamic disasters in coal mine, the numerical simulation is conducted to study the dynamic loading ef...In order to reveal the dynamic process of hard-thick roof inducing rock burst, one of the most common and strongest dynamic disasters in coal mine, the numerical simulation is conducted to study the dynamic loading effect of roof vibration on roadway surrounding rocks as well as the impact on stability. The results show that, on one hand, hard-thick roof will result in high stress concentration on mining surrounding rocks; on the other hand, the breaking of hard-thick roof will lead to mining seismicity, causing dynamic loading effect on coal and rock mass. High stress concentration and dynamic loading combination reaches to the mechanical conditions for the occurrence of rock burst, which will induce rock burst. The mining induced seismic events occurring in the roof breaking act on the mining surrounding rocks in the form of stress wave. The stress wave then has a reflection on the free surface of roadway and the tensile stress will be generated around the free surface. Horizontal vibration of roadway surrounding particles will cause instant changes of horizontal stress of roadway surrounding rocks; the horizontal displacement is directly related to the horizontal stress but is not significantly correlated with the vertical stress; the increase of horizontal stress of roadway near surface surrounding rocks and the release of elastic deformation energy of deep surrounding coal and rock mass are immanent causes that lead to the impact instability of roadway surrounding rocks. The most significant measures for rock burst prevention are controlling of horizontal stress and vibration strength.Key words展开更多
Rockbursts were frequently encountered in the construction of deeply buried tunnels at the Jinping-II hydropower station, Southwest China. In those cases, the existence of large structural planes, such as faults, was ...Rockbursts were frequently encountered in the construction of deeply buried tunnels at the Jinping-II hydropower station, Southwest China. In those cases, the existence of large structural planes, such as faults, was usually observed near the excavation boundaries. The formation mechanism of the “11·28” rockburst, which was a typical rockburst and occurred in a drainage tunnel under a deep burial depth, high in-situ stress state and complex geological conditions, has been difficult to explain. Realistic failure process analysis(RFPA3D) software was adopted to numerically simulate the whole failure process of the surrounding rock mass around the tunnel subjected to excavation. The spatial distribution of acoustic emission derived from numerical simulation contributed to explaining the mechanical responses of the process. Analyses of the stress, safety reserve coefficient and damage degree were performed to reveal the effect of faults on the formation of rockbursts in the deep tunnel. The existence of faults results in the formation of stress anomaly areas between the tunnel and the fault. The surrounding rock mass failure propagates toward the fault from the initial failure, to different degrees. The relative positions and angles of faults play significant roles in the extent and development of surrounding rock mass failure, respectively. The increase in the lateral stress coefficient leads to the aggravation of the surrounding rock mass damage, especially in the roof and floor of the tunnel. Moreover, as the rock strength-stress ratio increases, the failure mode of the near-fault tunnel gradually changes from the stress-controlled type to the compound-controlled type. These findings were consistent with the microseismic monitoring results and field observations, which was helpful to understand the mechanical behavior of tunnel excavation affected by faults. The achievements of this study can provide some references for analysis of the failure mechanisms of similar deep tunnels.展开更多
Fetr6 is an underground mine using the stope-and-pillar mining method. As there was some evidence regarding pillar failure in this mine, improving works such as roof support and replacing existing pillars with concret...Fetr6 is an underground mine using the stope-and-pillar mining method. As there was some evidence regarding pillar failure in this mine, improving works such as roof support and replacing existing pillars with concrete pillars (CP) were carried out. During the construction of the second CP, in the space between the remaining pillars, one of the pillars failed leading to the progressive failure of other pillars until 4 000 m 2 of mine had collapsed within a few minutes. In this work, this phenomenon is described by applying both numerical and empirical methods and the respective results are compared. The results of numerical modelling are found to be closer to the actual condition than those of the empirical method. Also, a width-to-height (W/H) ratio less than 1, an inadequate support system and the absence of a detailed program for pillar recovery are shown to be the most important causes of the Domino failure in this mine.展开更多
Pile foundation is widely used in the offshore engineering. The pile can be seriously destroyed by the soil liquefaction during strong earthquakes. The potentials of liquefaction and damages of pile foundation due to ...Pile foundation is widely used in the offshore engineering. The pile can be seriously destroyed by the soil liquefaction during strong earthquakes. The potentials of liquefaction and damages of pile foundation due to the liquefaction can be reduced by the implementation of the drainage in the liquefiable foundation. A patented pile technology, named rigidity-drain pile, was introduced. The partial section of the pile body was filled by materials with higher penetrability which forms some effective drainage channels in the pile. The principles and construction methods were presented. 3D models for both rigidity-drain pile and ordinary pile were built in FLAC3D code. The dynamic loadings were applied on the bottom of the model. According to the numerical results, in the case of the rigidity-drain pile, the water in the relevant distance range around the pile flows toward the pile drainage, the contour of the pore pressure shows a funnel form. Contrast to the ordinary pile, the rigidity-drain pile can dissipate the accumulated excess pore water, maintain effective stress and obviously reduce the possibility of surrounding soil liquefaction.展开更多
Foundation settlement is of great significance for high-fill engineering in collapsible loess areas.To predict the construction settlement of Lüliang Airport located in Shanxi Province,China,a plane strain finite...Foundation settlement is of great significance for high-fill engineering in collapsible loess areas.To predict the construction settlement of Lüliang Airport located in Shanxi Province,China,a plane strain finite element method considering the linear variation in the modulus,was carried out in this paper based on the results of geotechnical tests.The stress and deformation of four typical sections caused by layered fill are simulated,and then the settlement of the high-fill airport is calculated and analyzed by inputting three sets of parameters.The relative soft parameters of loess geomaterials produce more settlement than the relatively hard parameters.The thicker the filling body is,the greater the settlement is.The filling body constrained by mountains on both sides produces less settlement than the filling body constrained by a mountain on only one side even the filling thickness is almost the same.The settlement caused by the original subbase accounts for 56%−77%of the total settlement,while the fill soils themselves accounts for 23%−44%of the total settlement,which is approximately consistent with the field monitoring results.It provides a good reference for predicting the settlement of similar high-fill engineering.展开更多
In this study,the effect of varied loading ratio(mass of the explosive/mass of flyer plate)on the nature of interface,temperature and pressure developed in aluminum-steel explosive cladding is presented.Increase in th...In this study,the effect of varied loading ratio(mass of the explosive/mass of flyer plate)on the nature of interface,temperature and pressure developed in aluminum-steel explosive cladding is presented.Increase in the loading ratio,R,enhances the pressure developed,kinetic energy utilization and deformation work performed.Interfacial microstructures exhibit the formation of molten layer at few spots,owing to the increase in temperature beyond the melting point of parent alloy.The increase in temperature and the quantum of pressure developed were determined by numerical simulation performed in Ansys AUTODYN by employing smoothed particle hydrodynamics(SPH)method.The positioning of the experimental conditions on the weldability window is presented as well.展开更多
In order to control the dust pollution produced by air leg rock drill in the trolley area during the excavation of long-distance single ended tunnel,the full-scale physical model of working face was established by usi...In order to control the dust pollution produced by air leg rock drill in the trolley area during the excavation of long-distance single ended tunnel,the full-scale physical model of working face was established by using FLUENT software,and the numerical simulation analysis of tunnel drilling ventilation and dust removal parameters was carried out.The results show that it is difficult to control the dust pollution of the face by conventional ventilation,and the drilling dust is distributed in the range of 10 m from the face;after the introduction of the long pressure and short suction ventilation scheme,when the ratio of compressed air volume to exhaust air volume is 0.72,the height of the pressure fan is 2.5 m,the distance between the pressure fan and the palm face is 20 m,and the exhaust fan is 12 m away from the palm,the dust concentration control efficiency of the working face is increased by about 60%.Therefore,in the similar long-distance single head tunnel construction,it is appropriate to adopt the dust removal method of long-distance short suction and exhaust fan to ensure the working environment.展开更多
To predict three-dimensional temperature distribution of molten aluminum and its influencing factors inside an industrial aluminum holding furnace,a fluid-solid coupled method was presented.The fluid-solid coupled mat...To predict three-dimensional temperature distribution of molten aluminum and its influencing factors inside an industrial aluminum holding furnace,a fluid-solid coupled method was presented.The fluid-solid coupled mathematics models of aluminum holding furnace in the premixed combustion processing were established based on mass conservation,moment conservation,momentum conservation,energy conservation and chemistry species conservation.Computational results agree well with the test data of the typical condition.The maximum combustion temperature is 1 850 K.The average temperature of the molten aluminum is 1 158 K,and the maximum temperature difference is about 240 K.The average temperature increases 0.3 ℃ while the temperature of combustion air increases 1 ℃.The optimal excess air ratio is 1.25-1.30.展开更多
A stratum grouting-soil-structure interaction model which simplified the grouted zone into a series of spherical grout bulbs was established using FLAC3D program. The hypothetical non-uniform expansion process to reac...A stratum grouting-soil-structure interaction model which simplified the grouted zone into a series of spherical grout bulbs was established using FLAC3D program. The hypothetical non-uniform expansion process to reach an assigned volume strain due to soil compression by grouting was achieved by imposing radial velocity on outer mesh nodes of these spheres. This new method avoids the repeated trial calculation needed in the traditional method which applied a fictitious expanding pressure in the grouting element. The deformation and additional internal forces of structure were investigated during each grouting strategy and the influences of various stiffness of grouting proof curtain and bearing capacity of pile tip were discussed simultaneously. The numerical model is proved to be effective to replicate general behavior expected in the field and is capable of modeling the uplifting effect for the surface structure by grouting.展开更多
基金supported by the Guangdong Major Project of Basic and Applied Basic Research(Grant No.2021B0301030001)the National Key Research and Development Program of China(Grant No.2021YFB3802300)the Foundation of National Key Laboratory of Shock Wave and Detonation Physics(Grant No.JCKYS2022212004)。
文摘The graded density impactor(GDI)dynamic loading technique is crucial for acquiring the dynamic physical property parameters of materials used in weapons.The accuracy and timeliness of GDI structural design are key to achieving controllable stress-strain rate loading.In this study,we have,for the first time,combined one-dimensional fluid computational software with machine learning methods.We first elucidated the mechanisms by which GDI structures control stress and strain rates.Subsequently,we constructed a machine learning model to create a structure-property response surface.The results show that altering the loading velocity and interlayer thickness has a pronounced regulatory effect on stress and strain rates.In contrast,the impedance distribution index and target thickness have less significant effects on stress regulation,although there is a matching relationship between target thickness and interlayer thickness.Compared with traditional design methods,the machine learning approach offers a10^(4)—10^(5)times increase in efficiency and the potential to achieve a global optimum,holding promise for guiding the design of GDI.
基金supported by the Young Scientists Fund of National Natural Science Foundation of China(Grant Nos.12202202 and 12202494)the National Key Research and Development Program of China(Grant No.2021YFC3100700)。
文摘The internal and external flow fields during vented explosions of methane were characterized through numerical simulation,and the capability of numerical simulation thereof was validated by previous experimental data at three ignition positions.The venting mechanism was revealed by the simulated concentration distribution,temperature profile,and airflow velocity.The results show rear ignition results in the external methane mass distribution taking the form of"mushroom"and columnar flames in the external space,which can be expressed as a third-order polynomial relationship with distance;central ignition forms a relationship of the form y=AxB.Front ignition causes the temperature to show a tendency to repeated oscillations(rising,falling,and rising).Central ignition generates the maximum vented airflow velocity(V_(max)=320 m/s)upon vent opening.The results indicate that it is acceptable to apply numerical simulation of methane explosions in practice.
基金supported by the National Natural Science Foundation of China(Grant Nos.52171318 and 12202329)Joint Foundation of the Ministry of Education(Grant No.8091B022105)。
文摘The detonation of fuel-rich explosives yields combustible products that persistently burn upon mixing with ambient oxygen,releasing additional energy through a phenomenon known as the afterburning effect.This process greatly influences the evolution of confined blast loading and the subsequent structural response,which is crucial in confined blast scenarios.Given the complex nature of the reaction process,accurate analysis of the afterburning effect remains challenging.Previous studies have either overlooked the mechanisms of detonation product combustion or failed to provide experimental validation.This study introduces a three-dimensional model to effectively characterize the combustion of detonation products.The model integrates chemical reaction source terms into the governing equations to consider the combustion processes.Numerical simulations and experimental tests were conducted to analyze the combustion and energy release from the detonation products of fuel-rich explosives in confined spaces.Approximately 50%of the energy was released during the combustion of detonation products in a confined TNT explosion.Although the combustion of these products was much slower than the detonation process,it aligned with the dynamic response of the structure,which enhanced the explosive yield.Excluding afterburning from the analysis reduced the center-point deformation of the structure by 30%.Following the inclusion of afterburning,the simulated quasistatic pressure increased by approximately 45%.Subsequent comparisons highlighted the merits of the proposed approach over conventional methods.This approach eliminates the reliance on empirical parameters,such as the amount and rate of energy release during afterburning,thereby laying the foundation for understanding load evolution in more complex environments,such as ships,buildings,and underground tunnels.
基金financially supported by the National Natural Science Foundation of China(Grant Nos.52271317 and 52071149)the Fundamental Research Funds for the Central Universities(HUST:2019kfy XJJS007)。
文摘To address the problems of low accuracy by the CONWEP model and poor efficiency by the Coupled Eulerian-Lagrangian(CEL)method in predicting close-range air blast loads of cylindrical charges,a neural network-based simulation(NNS)method with higher accuracy and better efficiency was proposed.The NNS method consisted of three main steps.First,the parameters of blast loads,including the peak pressures and impulses of cylindrical charges with different aspect ratios(L/D)at different stand-off distances and incident angles were obtained by two-dimensional numerical simulations.Subsequently,incident shape factors of cylindrical charges with arbitrary aspect ratios were predicted by a neural network.Finally,reflected shape factors were derived and implemented into the subroutine of the ABAQUS code to modify the CONWEP model,including modifications of impulse and overpressure.The reliability of the proposed NNS method was verified by related experimental results.Remarkable accuracy improvement was acquired by the proposed NNS method compared with the unmodified CONWEP model.Moreover,huge efficiency superiority was obtained by the proposed NNS method compared with the CEL method.The proposed NNS method showed good accuracy when the scaled distance was greater than 0.2 m/kg^(1/3).It should be noted that there is no need to generate a new dataset again since the blast loads satisfy the similarity law,and the proposed NNS method can be directly used to simulate the blast loads generated by different cylindrical charges.The proposed NNS method with high efficiency and accuracy can be used as an effective method to analyze the dynamic response of structures under blast loads,and it has significant application prospects in designing protective structures.
文摘With different structural forms of ventilation pipes have various attenuation effects on incident shock waves while meeting ventilation requirements.The attenuation mechanism and the propagation law of shock waves in ventilation pipes of different structures are investigated by experiments and numerical simulations.Furthermore,for the same structure,the effects of peak pressure and positive pressure time on the attenuation rate are discussed.It is found that the attenuation rate increases with the incident shock wave pressure,and the shock wave attenuation rate tends to reach its limiting value k for the same structure and reasonably short positive pressure time.Under the same conditions,the attenuation rate is calculated using the pressure of the shock wave as follows:diffusion chamber pipe,branch pipe and selfconsumption pipe;the attenuation rate per unit volume is calculated as follows:self-consumption pipe,branch pipe and diffusion chamber pipe.In addition,an easy method is provided to calculate the attenuation rate of the shock wave in single and multi-stage ventilation pipes.Corresponding parameters are provided for various structures,and the margin of error between the formulae and experimental results is within 10%,which is significant for engineering applications.
文摘Most of the existing studies on tunnel blast wave are based on spherical or grouped charges, however,conventional weapons are mostly cylindrical rather than spherical. In order to analyze the impact of cylindrical charges on the tunnel blast wave loads and to develop a quantitative calculation method, this study carried out experimental and numerical research. Initially, external explosion experiments were conducted using both 35 kg spherical charges and cylindrical charges with aspect ratio of 4.8 at two different distances from the tunnel entrance. Comparative analysis of the blast wave parameters in the tunnel revealed that the explosive equivalent of the cylindrical charges was significantly higher than that of the spherical charges. To address this, an equivalent coefficient κ based on the spherical charges was proposed for the cylindrical charges. Subsequently, numerical simulations were conducted for the experimental conditions, and the numerical simulation results match the experiments well. Through numerical calculations, the reliability of the equivalent coefficient κ under the experimental conditions was verified, and comparison analysis indicated that the explosion energy of cylindrical charges spreads more radially, resulting in more explosion energy entering the tunnel, which is the fundamental reason for the increase in tunnel blast wave loads. Additionally, analyzing the explosion energy ratio entering the tunnel is an effective method for calculating the equivalent coefficient κ. Finally, through more than one hundred sets of numerical calculation results, the impact of the proportional distance λ and the ratio of charge mass to the tunnel cross-section dimension φ on the equivalence coefficients κ was investigated. An empirical formula for the equivalence coefficient κ was derived through fitting, and the accuracy of the formula was validated through literature experimental results. The research findings of this paper will provide valuable guidance for the calculation of blast wave loads in tunnel.
基金Project(41630642)supported by the Key Project of National Natural Science Foundation of ChinaProject(51974360)supported by the National Natural Science Foundation of ChinaProject(2018JJ3656)supported by the Natural Science Foundation of Hunan Province,China。
文摘In the context of deep rock engineering,the in-situ stress state is of major importance as it plays an important role in rock dynamic response behavior.Thus,stress initialization becomes crucial and is the first step for the dynamic response simulation of rock mass in a high in-situ stress field.In this paper,stress initialization methods,including their principles and operating procedures for reproducing steady in-situ stress state in LS-DYNA,are first introduced.Then the most popular four methods,i.e.,explicit dynamic relaxation(DR)method,implicit-explicit sequence method,Dynain file method and quasi-static method,are exemplified through a case analysis by using the RHT and plastic hardening rock material models to simulate rock blasting under in-situ stress condition.Based on the simulations,it is concluded that the stress initialization results obtained by implicit-explicit sequence method and dynain file method are closely related to the rock material model,and the explicit DR method has an obvious advantage in solution time when compared to other methods.Besides that,it is recommended to adopt two separate analyses for the whole numerical simulation of rock mass under the combined action of in-situ stress and dynamic disturbance.
基金funded by the National Natural Science Foundation of China under NO.11202103Qing-lan Project of Jiangsu Province。
文摘The numerical simulation of a blast wave of a multilayer composite charge is investigated.A calculation model of the near-field explosion and far-field propagation of the shock wave of a composite charge is established using the AUTODYN finite element program.Results of the near-field and far-field calculations of the shock wave respectively converge at cell sizes of 0.25-0.5 cm and 1-3 cm.The Euler--fluxcorrected transport solver is found to be suitable for the far-field calculation after mapping.A numerical simulation is conducted to study the formation,propagation,and interaction of the shock wave of the composite charge for different initiation modes.It is found that the initiation mode obviously affects the shock-wave waveform and pressure distribution of the composite charge.Additionally,it is found that the area of the overpressure distribution is greatest for internal and external simultaneous initiation,and the peak pressure of the shock wave exponentially decays,fitting the calculation formula of the peak overpressure attenuation under different initiation modes,which is obtained and verified by experiment.The difference between numerical and experimental results is less than 10%,and the peak overpressure of both internal and external initiation is 56.12% higher than that of central single-point initiation.
文摘The design of solid armature of railgun should take full account of its operating conditions and material properties because the armature is subjected to dynamic loading conditions and experiences a complicated electrical,thermal and mechanical process in the interior ballistic cycle.In this paper present,we first introduced a multi-physical field model of railgun,followed by several examples to investigate the launching process.Especially,we used the explicit finite element method,in which material nonlinearity and geometric nonlinearity were accounted,to investigate the deform behaviors of solid armature.The results show that the dynamic mechanical process of armature is dependent on the armature geometry,material and exciting electric current.By the numerical simulation,the understanding of the fracture mechanism of solid armature was deepened.
基金Projects(52074166,51774195,51704185)supported by the National Natural Science Foundation of ChinaProject(2019M652436)supported by the China Postdoctoral Science Foundation。
文摘In actual production,deep coal mine roadways are often under typical static-dynamic coupling stress(SDCS)conditions with high ground stress and strong dynamic disturbances.With the increasing number of disasters and accidents induced by SDCS conditions,the safe and efficient production of coal mines is seriously threatened.Therefore,it is of great practical significance to study the deformation and failure characteristics of the roadway surrounding rock under SDCS.In this paper,the effects of different in-situ stress fields and dynamic load conditions on the surrounding rock are studied by numerical simulations,and the deformation and failure characteristics are obtained.According to the simulation results,the horizontal stress,vertical stress and dynamic disturbance have a positive correlation with the plastic failure of the surrounding rock.Among these factors,the influence of the dynamic disturbance is the most substantial.Under the same stress conditions,the extents of deformation and plastic failure of the roof and ribs are always greater than those of the floor.The effect of horizontal stresses on the roadway deformation is more notable than that of vertical stresses.The results indicate that for the roadway under high-stress conditions,the in-situ stress test must be strengthened first.After determining the magnitude of the in-situ stress,the location of the roadway should be reasonably arranged in the design to optimize the mining sequence.For roadways that are strongly disturbed by dynamic loads,rock supports(rebar/cable bolts,steel set etc.)that are capable of maintaining their effectiveness without failure after certain dynamic loads are required.The results of this study contribute to understanding the characteristics of the roadway deformation and failure under SDCS,and can be used to provide a basis for the support design and optimization under similar geological and geotechnical circumstances.
基金Project(51375498)supported by the National Natural Science Foundation of China
文摘Based on RNG k-ε turbulence model and sliding grid technique, solid-liquid two-phase three-dimensional(3-D) unsteady turbulence of full passage in slurry pump was simulated by means of Fluent software. The effects of unsteady flow characteristics on solid-liquid two-phase flow and pump performance were researched under design condition. The results show that clocking effect has a significant influence on the flow in pump, and the fluctuation of flow velocity and pressure is obvious, particularly near the volute tongue, at the position of small sections of volute and within diffuser. Clocking effect has a more influence on liquid-phase than on solid-phase, and the wake-jet structure of relative velocity of solid-phase is less obvious than liquid-phase near the volute tongue and the impeller passage outlet. The fluctuation of relative velocity of solid-phase flow is 7.6% smaller than liquid-phase flow at the impeller outlet on circular path. Head and radial forces of the impeller are 8.1% and 85.7% of fluctuation, respectively. The results provide a theoretical basis for further research for turbulence, improving efficient, reducing the hydraulic losses and wear. Finally, field tests were carried out to verify the operation and wear of slurry pump.
基金Project(51404243)supported by the National Natural Science Foundation of ChinaProject(2014QNB26)supported by the Fundamental Research Funds for the Central Universities,China
文摘In order to reveal the dynamic process of hard-thick roof inducing rock burst, one of the most common and strongest dynamic disasters in coal mine, the numerical simulation is conducted to study the dynamic loading effect of roof vibration on roadway surrounding rocks as well as the impact on stability. The results show that, on one hand, hard-thick roof will result in high stress concentration on mining surrounding rocks; on the other hand, the breaking of hard-thick roof will lead to mining seismicity, causing dynamic loading effect on coal and rock mass. High stress concentration and dynamic loading combination reaches to the mechanical conditions for the occurrence of rock burst, which will induce rock burst. The mining induced seismic events occurring in the roof breaking act on the mining surrounding rocks in the form of stress wave. The stress wave then has a reflection on the free surface of roadway and the tensile stress will be generated around the free surface. Horizontal vibration of roadway surrounding particles will cause instant changes of horizontal stress of roadway surrounding rocks; the horizontal displacement is directly related to the horizontal stress but is not significantly correlated with the vertical stress; the increase of horizontal stress of roadway near surface surrounding rocks and the release of elastic deformation energy of deep surrounding coal and rock mass are immanent causes that lead to the impact instability of roadway surrounding rocks. The most significant measures for rock burst prevention are controlling of horizontal stress and vibration strength.Key words
基金Project(42177143) supported by the National Natural Science Foundation of ChinaProject(2020JDJQ0011) supported by the Science Foundation for Distinguished Young Scholars of Sichuan Province,China。
文摘Rockbursts were frequently encountered in the construction of deeply buried tunnels at the Jinping-II hydropower station, Southwest China. In those cases, the existence of large structural planes, such as faults, was usually observed near the excavation boundaries. The formation mechanism of the “11·28” rockburst, which was a typical rockburst and occurred in a drainage tunnel under a deep burial depth, high in-situ stress state and complex geological conditions, has been difficult to explain. Realistic failure process analysis(RFPA3D) software was adopted to numerically simulate the whole failure process of the surrounding rock mass around the tunnel subjected to excavation. The spatial distribution of acoustic emission derived from numerical simulation contributed to explaining the mechanical responses of the process. Analyses of the stress, safety reserve coefficient and damage degree were performed to reveal the effect of faults on the formation of rockbursts in the deep tunnel. The existence of faults results in the formation of stress anomaly areas between the tunnel and the fault. The surrounding rock mass failure propagates toward the fault from the initial failure, to different degrees. The relative positions and angles of faults play significant roles in the extent and development of surrounding rock mass failure, respectively. The increase in the lateral stress coefficient leads to the aggravation of the surrounding rock mass damage, especially in the roof and floor of the tunnel. Moreover, as the rock strength-stress ratio increases, the failure mode of the near-fault tunnel gradually changes from the stress-controlled type to the compound-controlled type. These findings were consistent with the microseismic monitoring results and field observations, which was helpful to understand the mechanical behavior of tunnel excavation affected by faults. The achievements of this study can provide some references for analysis of the failure mechanisms of similar deep tunnels.
文摘Fetr6 is an underground mine using the stope-and-pillar mining method. As there was some evidence regarding pillar failure in this mine, improving works such as roof support and replacing existing pillars with concrete pillars (CP) were carried out. During the construction of the second CP, in the space between the remaining pillars, one of the pillars failed leading to the progressive failure of other pillars until 4 000 m 2 of mine had collapsed within a few minutes. In this work, this phenomenon is described by applying both numerical and empirical methods and the respective results are compared. The results of numerical modelling are found to be closer to the actual condition than those of the empirical method. Also, a width-to-height (W/H) ratio less than 1, an inadequate support system and the absence of a detailed program for pillar recovery are shown to be the most important causes of the Domino failure in this mine.
基金Project (50639010) supported by the National Natural Science Foundation of China
文摘Pile foundation is widely used in the offshore engineering. The pile can be seriously destroyed by the soil liquefaction during strong earthquakes. The potentials of liquefaction and damages of pile foundation due to the liquefaction can be reduced by the implementation of the drainage in the liquefiable foundation. A patented pile technology, named rigidity-drain pile, was introduced. The partial section of the pile body was filled by materials with higher penetrability which forms some effective drainage channels in the pile. The principles and construction methods were presented. 3D models for both rigidity-drain pile and ordinary pile were built in FLAC3D code. The dynamic loadings were applied on the bottom of the model. According to the numerical results, in the case of the rigidity-drain pile, the water in the relevant distance range around the pile flows toward the pile drainage, the contour of the pore pressure shows a funnel form. Contrast to the ordinary pile, the rigidity-drain pile can dissipate the accumulated excess pore water, maintain effective stress and obviously reduce the possibility of surrounding soil liquefaction.
基金Project(2020M670604)supported by the China Postdoctoral Science FoundationProject(41790434)supported by the National Natural Science Foundation of China。
文摘Foundation settlement is of great significance for high-fill engineering in collapsible loess areas.To predict the construction settlement of Lüliang Airport located in Shanxi Province,China,a plane strain finite element method considering the linear variation in the modulus,was carried out in this paper based on the results of geotechnical tests.The stress and deformation of four typical sections caused by layered fill are simulated,and then the settlement of the high-fill airport is calculated and analyzed by inputting three sets of parameters.The relative soft parameters of loess geomaterials produce more settlement than the relatively hard parameters.The thicker the filling body is,the greater the settlement is.The filling body constrained by mountains on both sides produces less settlement than the filling body constrained by a mountain on only one side even the filling thickness is almost the same.The settlement caused by the original subbase accounts for 56%−77%of the total settlement,while the fill soils themselves accounts for 23%−44%of the total settlement,which is approximately consistent with the field monitoring results.It provides a good reference for predicting the settlement of similar high-fill engineering.
文摘In this study,the effect of varied loading ratio(mass of the explosive/mass of flyer plate)on the nature of interface,temperature and pressure developed in aluminum-steel explosive cladding is presented.Increase in the loading ratio,R,enhances the pressure developed,kinetic energy utilization and deformation work performed.Interfacial microstructures exhibit the formation of molten layer at few spots,owing to the increase in temperature beyond the melting point of parent alloy.The increase in temperature and the quantum of pressure developed were determined by numerical simulation performed in Ansys AUTODYN by employing smoothed particle hydrodynamics(SPH)method.The positioning of the experimental conditions on the weldability window is presented as well.
基金Project(51874016)supported by the National Natural Science Foundation of China。
文摘In order to control the dust pollution produced by air leg rock drill in the trolley area during the excavation of long-distance single ended tunnel,the full-scale physical model of working face was established by using FLUENT software,and the numerical simulation analysis of tunnel drilling ventilation and dust removal parameters was carried out.The results show that it is difficult to control the dust pollution of the face by conventional ventilation,and the drilling dust is distributed in the range of 10 m from the face;after the introduction of the long pressure and short suction ventilation scheme,when the ratio of compressed air volume to exhaust air volume is 0.72,the height of the pressure fan is 2.5 m,the distance between the pressure fan and the palm face is 20 m,and the exhaust fan is 12 m away from the palm,the dust concentration control efficiency of the working face is increased by about 60%.Therefore,in the similar long-distance single head tunnel construction,it is appropriate to adopt the dust removal method of long-distance short suction and exhaust fan to ensure the working environment.
基金Project(2006AA03Z523) supported by the National High-Tech Research and Development Program of ChinaProject(08C26224302178) supported by the Innovation Foundation of Central South University,China
文摘To predict three-dimensional temperature distribution of molten aluminum and its influencing factors inside an industrial aluminum holding furnace,a fluid-solid coupled method was presented.The fluid-solid coupled mathematics models of aluminum holding furnace in the premixed combustion processing were established based on mass conservation,moment conservation,momentum conservation,energy conservation and chemistry species conservation.Computational results agree well with the test data of the typical condition.The maximum combustion temperature is 1 850 K.The average temperature of the molten aluminum is 1 158 K,and the maximum temperature difference is about 240 K.The average temperature increases 0.3 ℃ while the temperature of combustion air increases 1 ℃.The optimal excess air ratio is 1.25-1.30.
基金Project(2007AA11Z134)supported by the National High Technology Research and Development Program of ChinaProject(10JJ4035)supported by the Natural Science Foundation of Hunan Province,ChinaProject(2010ybfz046)supported by the Fund of Excellent Doctoral Dissertation of Central South University,China
文摘A stratum grouting-soil-structure interaction model which simplified the grouted zone into a series of spherical grout bulbs was established using FLAC3D program. The hypothetical non-uniform expansion process to reach an assigned volume strain due to soil compression by grouting was achieved by imposing radial velocity on outer mesh nodes of these spheres. This new method avoids the repeated trial calculation needed in the traditional method which applied a fictitious expanding pressure in the grouting element. The deformation and additional internal forces of structure were investigated during each grouting strategy and the influences of various stiffness of grouting proof curtain and bearing capacity of pile tip were discussed simultaneously. The numerical model is proved to be effective to replicate general behavior expected in the field and is capable of modeling the uplifting effect for the surface structure by grouting.
