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 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.展开更多
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.展开更多
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 this study, a three dimensional(3D) numerical model of six-degrees-of-freedom(6DOF) is applied to simulate the water entries of twin spheres side-by-side at different lateral distances and time intervals.The turbul...In this study, a three dimensional(3D) numerical model of six-degrees-of-freedom(6DOF) is applied to simulate the water entries of twin spheres side-by-side at different lateral distances and time intervals.The turbulence structure is described using the shear-stress transport k-ω(SST k-ω) model, and the volume of fluid(VOF) method is used to track the complex air-liquid interface. The motion of spheres during water entry is simulated using an independent overset grid. The numerical model is verified by comparing the cavity evolution results from simulations and experiments. Numerical results reveal that the time interval between the twin water entries evidently affects cavity expansion and contraction behaviors in the radial direction. However, this influence is significantly weakened by increasing the lateral distance between the two spheres. In synchronous water entries, pressure is reduced on the midline of two cavities during surface closure, which is directly related to the cavity volume. The evolution of vortexes inside the two cavities is analyzed using a velocity vector field, which is affected by the lateral distance and time interval of water entries.展开更多
The complex structure of the bottom of a high-speed train is an important source of train aerodynamic drag.Thus,improving the bottom structure is of great significance to reduce the aerodynamic drag of the train.In th...The complex structure of the bottom of a high-speed train is an important source of train aerodynamic drag.Thus,improving the bottom structure is of great significance to reduce the aerodynamic drag of the train.In this study,computational fluid dynamics(CFD)based on three-dimensional steady incompressible Reynolds-average Naiver-Stokes(RANS)equations and Realizable k-εturbulence model were utilized for numerical simulations.Inspired by the concept of streamlined design and the idea of bottom flow field control,this study iteratively designed the bogies in a streamlined shape and combined them with the bottom deflectors to investigate the joint drag reduction mechanism.Three models,i.e.,single-bogie model,simplified train model,and eight-car high-speed train model,were created and their aerodynamic characteristics were analyzed.The results show that the single-bogie model with streamlined design shows a noticeable drag reduction,whose power bogie and trailer bogie experience 13.92%and 7.63%drag reduction,respectively.The range of positive pressure area on the bogie is reduced.The aerodynamic drag can be further reduced to 15.01%by installing both the streamlined bogie and the deflector on the simplified train model.When the streamlined bogies and deflectors are used on the eight-car model together,the total drag reduction rate reaches 2.90%.Therefore,the proposed aerodynamic kit for the high-speed train bottom is capable to improve the flow structure around the bogie regions,reduce the bottom flow velocity,and narrow the scope of the train’s influence on the surrounding environment,achieving the appreciable reduction of aerodynamic drag.This paper can provide a new idea for the drag reduction of high-speed trains.展开更多
The utilization of prefabricated light modular radiant heating system has demonstrated significant increases in heat transfer efficiency and energy conservation capabilities.Within prefabricated building construction,...The utilization of prefabricated light modular radiant heating system has demonstrated significant increases in heat transfer efficiency and energy conservation capabilities.Within prefabricated building construction,this new heating method presents an opportunity for the development of comprehensive facilities.The parameters for evaluating the effectiveness of such a system are the upper surface layer’s heat flux and temperature.In this paper,thermal resistance analysis calculation based on a simplified model for this unique radiant heating system analysis is presented with the heat transfer mechanism’s evaluation.The results obtained from thermal resistance analysis calculation and numerical simulation indicate that the thermal resistance analysis method is highly accurate with temperature discrepancies ranging from 0.44℃ to−0.44℃ and a heat flux discrepancy of less than 7.54%,which can meet the requirements of practical engineering applications,suggesting a foundation for the prefabricated radiant heating system.展开更多
In order to improve the overall resilience of the urban infrastructures, it is required to conduct blast resistant design for important building structures in the city. For complex terrain in the city, it is recommend...In order to improve the overall resilience of the urban infrastructures, it is required to conduct blast resistant design for important building structures in the city. For complex terrain in the city, it is recommended to determine the blast load on the structures via numerical simulation. Since the mesh size of the numerical model highly depends on the explosion scenario, there is no generally applicable approach for the mesh size selection. An efficient method to determine the mesh size of the numerical model of near-ground detonation based on explosion scenarios is proposed in this study. The effect of mesh size on the propagation of blast wave under different explosive weights was studied, and the correlations between the mesh size effect and the charge weight or the scaled distance was described. Based on the principle of the finite element method and Hopkinson-Cranz scaling law, a mesh size measurement unit related to the explosive weight was proposed as the criterion for determining the mesh size in the numerical simulation. Finally, the applicability of the method proposed in this paper was verified by comparing the results from numerical simulation and the explosion tests and was verified in AUTODYN.展开更多
The reverse magnetohydrodynamic(MHD)energy bypass technology is a promising energy redis⁃tribution technology in the scramjet system,in augmented with a power generation equipment to supply the neces⁃sary long-distanc...The reverse magnetohydrodynamic(MHD)energy bypass technology is a promising energy redis⁃tribution technology in the scramjet system,in augmented with a power generation equipment to supply the neces⁃sary long-distance flight airframe power.In this paper,a computational model of the scramjet magnetohydrody⁃namic channel is developed and verified by using the commercial software Fluent.It is found that when the mag⁃netic induction intensity is 1,2,3,4 T,the power generation efficiency is 22.5%,22.3%,22.0%,21.5%,and decreases with the increase of the magnetic induction intensity,and the enthalpy extraction rate is 0.026%,0.1%,0.21%,0.34%,and increases with the increase of the magnetic induction intensity.The deceleration ef⁃fect of electromagnetic action on the airflow in the power channel increases with the increase of magnetic induc⁃tion intensity.The stronger the magnetic field intensity,the more obvious the decreasing effect of fluid Mach num⁃ber in the channel.The power generation efficiency decreases as the magnetic induction intensity increases and the enthalpy extraction rate is reversed.As the local currents gathering at inlet and outlet of the power generation area,total temperature and enthalpy along the flow direction do not vary linearly,and there are maximum and minimum values at inlet and outlet.Increasing the number of electrodes can effectively regulate the percentage of Joule heat dissipation,which can improve the power generation efficiency.展开更多
This paper explores the performances of a finite element simulation including four concrete models applied to a full-scale reinforced concrete beam subjected to blast loading. Field test data has been used to compare ...This paper explores the performances of a finite element simulation including four concrete models applied to a full-scale reinforced concrete beam subjected to blast loading. Field test data has been used to compare model results for each case. The numerical modelling has been, carried out using the suitable code LS-DYNA. This code integrates blast load routine(CONWEP) for the explosive description and four different material models for the concrete including: Karagozian & Case Concrete, Winfrith, Continuous Surface Cap Model and Riedel-Hiermaier-Thoma models, with concrete meshing based on 10, 15, and 20 mm. Six full-scale beams were tested: four of them used for the initial calibration of the numerical model and two more tests at lower scaled distances. For calibration, field data obtained employing pressure and accelerometers transducers were compared with the results derived from the numerical simulation. Damage surfaces and the shape of rupture in the beams have been used as references for comparison. Influence of the meshing on accelerations has been put in evidence and for some models the shape and size of the damage in the beams produced maximum differences around 15%. In all cases, the variations between material and mesh models are shown and discussed.展开更多
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 polyur...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.展开更多
During the propagation of high-power lasers within internal channels,the laser beam heats the propagation medium,causing the thermal blooming effect that degrades the beam quality at the output.The intricate configura...During the propagation of high-power lasers within internal channels,the laser beam heats the propagation medium,causing the thermal blooming effect that degrades the beam quality at the output.The intricate configuration of the optical path within the internal channel necessitates complex and time-consuming efforts to assess the impact of thermal blooming effect on the optical path.To meet the engineering need for rapid evaluation of thermal blooming effect in optical paths,this study proposed a rapid simulation method for the thermal blooming effect in internal optical paths based on the finite element method.This method discretized the fluid region into infinitesimal elements and employed finite element method for flow field analysis.A simplified analytical model of the flow field region in complex internal channels was established,and regions with similar thermal blooming effect were divided within this model.Based on the calculated optical path differences within these regions,numerical simulations of phase distortion caused by thermal blooming were conducted.The calculated result were compared with those obtained using the existing methods.The findings reveal that for complex optical paths,the discrepancy between the two approaches is less than 3.6%,with similar phase distortion patterns observed.For L-type units,this method and the existing methods identify the same primary factors influencing aberrations and exhibit consistent trends in their variation.This method was used to analyze the impact of thermal blooming effect in a straight channel under different gravity directions.The results show that phase distortion varies with changes in the direction of gravity,and the magnitude of the phase difference is strongly correlated with the component of gravity perpendicular to the optical axis.Compared to the existing methods,this approach offers greater flexibility,obviates the need for complex custom analysis programming.The analytical results of this method enable a rapid assessment of the thermal blooming effect in optical paths within the internal channel.This is especially useful during the engineering design.These results also provide crucial references for developing strategies to suppress thermal blooming effect.展开更多
[Objective]This study aims to develop a thermodynamically consistent phase-field framework for modeling the initiation and evolution of discontinuous structures in geomaterials.[Methods]Our model introduces crack driv...[Objective]This study aims to develop a thermodynamically consistent phase-field framework for modeling the initiation and evolution of discontinuous structures in geomaterials.[Methods]Our model introduces crack driving forces derived from the volumetric-deviatoric strain decomposition strategy,incorporating distinct tension,compression,and shear degradation mechanisms.Inertia effects capture compaction-band formation driven by wave-like disturbances,grain crushing,and frictional rearrangement.A monolithic algorithm ensures numerical stability and rapid convergence.[Results]The framework reproduces tensile,shear,mixed tensile-shear,and compressive-shear failures using the Benzeggagh-Kenane criterion.Validation against benchmark simulations-including uniaxial compression of rock-like and triaxial compression of V-notched sandstone specimens-demonstrates accurate predictions of crack initiation stress,localization orientation,and energy dissipation.[Conclusions]The framework provides a unified and robust numerical tool for analyzing the spatiotemporal evolution of strain localization and fracture in geomaterials.[Significance]By linking microscale fracture dynamics with macroscale failure within a thermodynamically consistent scheme,this study advances predictive modeling of rock stability,slope failure,and subsurface energy systems,contributing to safer and more sustainable geotechnical practice.展开更多
This study focuses on the thermal management of 4680-type cylindrical lithium-ion battery packs utilizing NCM811 chemistry.It establishes coupled multi-physics models for both immersion and serpentine cold plate cooli...This study focuses on the thermal management of 4680-type cylindrical lithium-ion battery packs utilizing NCM811 chemistry.It establishes coupled multi-physics models for both immersion and serpentine cold plate cooling systems.Through a combination of numerical simulation and experimental validation,the technical advantages and mechanisms of immersion cooling are systematically explored.Simulation results indicate that under a 3C fast-charging condition(inlet temperature 20℃,flow rate 36 L/min),the immersion cooling structure 3demonstrates a triple enhancement in thermal performance compared to the cold plate structure 1:a 13.06%reduction in peak temperature,a 31.67%decrease in overall maximum temperature difference,and a 47.62%decrease in single-cell temperature deviation,while also reducing flow resistance by 33.61%.Furthermore,based on the immersion cooling model,a small battery module comprising seven cylindrical cells was designed for thermal runaway testing via nail penetration.The results show that the peak temperature of the triggered cell was limited to 437.6℃,with a controllable temperature rise gradient of only 3.35℃/s and a rapid cooling rate of 0.6℃/s.The maximum temperature rise of adjacent cells was just 64.8℃,effectively inhibiting thermal propagation.Post-test disassembly revealed that the non-triggered cells retained>99.2%of their original voltage and>99%structural integrity,confirming the module’s ability to achieve“localized failure with global stability.”展开更多
The finite volume method was applied to numerically simulate the bottom pressure field induced by regular waves,vehicles in calm water and vehicles in regular waves.The solution of Navier-Stokes(N-S)equations in the v...The finite volume method was applied to numerically simulate the bottom pressure field induced by regular waves,vehicles in calm water and vehicles in regular waves.The solution of Navier-Stokes(N-S)equations in the vicinity of numerical wave tank's boundary was forced towards the wave theoretical solution by incorporating momentum source terms,thereby reducing adverse effects such as wave reflection.Simulations utilizing laminar flow,turbulent flow,and ideal fluid models were all found capable of effectively capturing the waveform and bottom pressure of regular waves,agreeing well with experimental data.In predicting the bottom pressure field of the submerged vehicle,turbulent simulations considering fluid viscosity and boundary layer development provided more accurate predictions for the stern region than inviscid simulations.Due to sphere's diffractive effect,the sphere's bottom pressure field in waves is not a linear superposition of the wave's and the sphere's bottom pressure field.However,a slender submerged vehicle exhibits a weaker diffractive effect on waves,thus the submerged vehicle's bottom pressure field in waves can be approximated as a linear superposition of the wave's and the submerged vehicle's bottom pressure field,which simplifies computation and analysis.展开更多
During upward horizontal stratified backfill mining,stable backfill is essential for cap and sill pillar recovery.Currently,the primary method for calculating the required strength of backfill is the generalized three...During upward horizontal stratified backfill mining,stable backfill is essential for cap and sill pillar recovery.Currently,the primary method for calculating the required strength of backfill is the generalized three-dimensional(3 D)vertical stress model,which ignores the effect of mine depth,failing to obtain the vertical stress at different positions along stope length.Therefore,this paper develops and validates an improved 3 D model solution through numerical simulation in Rhino-FLAC^(3D),and examines the stress state and stability of backfill under different conditions.The results show that the improved model can accurately calculate the vertical stress at different mine depths and positions along stope length.The error rates between the results of the improved model and numerical simulation are below 4%,indicating high reliability and applicability.The maximum vertical stress(σ_(zz,max))in backfill is positively correlated with the degree of rock-backfill closure,which is enhanced by mine depth and elastic modulus of backfill,while weakened by stope width and inclination,backfill friction angle,and elastic modulus of rock mass.Theσ_(zz,max)reaches its peak when the stope length is 150 m,whileσ_(zz,max)is insensitive to changes in rock-backfill interface parameters.In all cases,the backfill stability can be improved by reducingσ_(zz,max).The results provide theoretical guidance for the backfill strength design and the safe and efficient recovery of ore pillars in deep mining.展开更多
Whipple shields as sacrificial bumpers,safeguard the satellites against extremely fast,different-sized projectiles traveling through space in the low earth orbit.Typical Whipple shields comprise a front and rear plate...Whipple shields as sacrificial bumpers,safeguard the satellites against extremely fast,different-sized projectiles traveling through space in the low earth orbit.Typical Whipple shields comprise a front and rear plate,separated by a gap or space.Recent advancements have explored the use of foam,cellular cores,and alternative materials such as ceramics instead of aluminium for the plates.In the current work,the effect of including fluid cores(air/water)sandwiched between the front and rear plates,on the response to hypervelocity impact was explored through a numerical approach.The numerical simulation consisted of hypervelocity impact by a 2 mm diameter,stainless steel projectile,launched at speeds of 3 e9 km/s with a normal impact trajectory towards the Whipple shield.The front and rear bumpers,made of AA6061-T6,were each 1 mm thick.A space of 10 mm was taken between the plates(occupied by fluid).The key metrics analyzed were the perforation characteristics,stages of the debris cloud generation and propagation,energy variations(internal,kinetic and plastic work),temperature variations,and the fragmentation summary.From the computational analysis,employing water-core in Whipple shields could prevent the rear bumper perforation till 6 km/s,lower the peak temperatures at the front bumper perforation zones and debris tip,and generate fewer,larger fragments.展开更多
Investigating the blast effects and mechanisms on typical finite-sized obstacles is essential for optimizing defense strategies and designing more robust barriers to deter terrorists and protect critical locations.Thi...Investigating the blast effects and mechanisms on typical finite-sized obstacles is essential for optimizing defense strategies and designing more robust barriers to deter terrorists and protect critical locations.This study investigates the blasting effects and underlying mechanisms of concrete frustums subjected to contact explosions,employing both numerical simulations and field tests.It focuses on the effects of top and side blasting,with particular emphasis on fracture modes,damage patterns,and fragment sizes,as well as the causes of different failure modes and the propagation of stress waves.The study also explores the blasting effects of detonating explosives at varying positions along the side and with different charge amounts.The results show that side-blasting leads to complete fragmentation,with tensile waves playing a significant role in creating extensive damage zones that propagate parallel to the frustum's outer surface,concentrating damage near the surface.During top-blasting,the upper half of the frustum undergoes fragmentation,while the lower half experiences cracking.Tensile waves propagate from the top to the bottom surface,forming larger blocks in regions with lower wave intensity.Three distinct damage zones within the frustum were identified,and a series of mathematical formulas were derived to describe the relationship between the maximum fragment size and charge mass.As the charge mass increased from 1.0 kg to 4.0 kg,the maximum fragment size decreased.Detonation at the center of the frustum's side resulted in the most severe fragmentation,with a 51.8%reduction in fragment size compared to other detonation positions.Finally,four broken modes were classified,each influenced by charge mass and explosive location.This study provides valuable insights for optimizing civil blasting operations and designing protective engineering structures.展开更多
To explore the design criteria for composite charges and reveal the intrinsic relationship between the detonation wave propagation in composite charges and the overall energy distribution of shock waves,this study ana...To explore the design criteria for composite charges and reveal the intrinsic relationship between the detonation wave propagation in composite charges and the overall energy distribution of shock waves,this study analyzes the propagation and interaction processes of detonation waves in composite charges with different structural dimensions and explosive combinations. It also investigates the spatial distribution characteristics of the resulting shock wave loads. Based on dimensional analysis theory, a theoretical analysis of the shock wave overpressure distribution in free air fields is conducted. Utilizing the derived dimensionless function relationships, the hydrocode AUTODYN is employed to investigate the effects of charge structure parameters and explosive combinations on the internal overdriven detonation phenomena and the distribution of shock wave loads. It is found that the overdriven detonation phenomenon in the inner layer of composite charges increases the strength of the axial detonation wave,thereby enhancing the intensity of the primary end wave formed upon refraction into the air, which affects the distribution characteristics of the shock wave overpressure. Research has shown that increasing the thickness ratio and detonation velocity ratio of composite charges is beneficial for exacerbating the phenomenon of overdriven detonation, improving the primary end wave intensity and axial overpressure. This gain effect gradually weakens with the propagation of shock waves. When overdriven detonation occurs inside the composite charge, the detonation pressure first increases and then decreases. The Mach reflection pressure of the composite charge with a larger aspect ratio is attenuated to a greater extent. In addition, as the aspect ratio of the composite charge increases, the shock wave energy gradually flows from the axial direction to the radial direction. Therefore, as the aspect ratio of the composite charge increases, the primary end wave intensity and axial overpressure gradually decrease.展开更多
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.展开更多
基金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 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.
基金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.
文摘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.
基金China Academy of Launch Vehicle Technology(Grant No.CALT-2022-03)Science and Technology on Underwater Information and Control Laboratory(Grant No.2021-JCJQ-LB-030-05).
文摘In this study, a three dimensional(3D) numerical model of six-degrees-of-freedom(6DOF) is applied to simulate the water entries of twin spheres side-by-side at different lateral distances and time intervals.The turbulence structure is described using the shear-stress transport k-ω(SST k-ω) model, and the volume of fluid(VOF) method is used to track the complex air-liquid interface. The motion of spheres during water entry is simulated using an independent overset grid. The numerical model is verified by comparing the cavity evolution results from simulations and experiments. Numerical results reveal that the time interval between the twin water entries evidently affects cavity expansion and contraction behaviors in the radial direction. However, this influence is significantly weakened by increasing the lateral distance between the two spheres. In synchronous water entries, pressure is reduced on the midline of two cavities during surface closure, which is directly related to the cavity volume. The evolution of vortexes inside the two cavities is analyzed using a velocity vector field, which is affected by the lateral distance and time interval of water entries.
基金Project(2020YFA0710901)supported by the National Key Research and Development Program of ChinaProject(2023JJ30643)supported by the Natural Science Foundation of Hunan Province,China+1 种基金Project(12372204)supported by the National Natural Science Foundation of ChinaProject(2022ZZTS0725)supported by the Self-exploration and Innovation Project for Postgraduates of Central South University,China。
文摘The complex structure of the bottom of a high-speed train is an important source of train aerodynamic drag.Thus,improving the bottom structure is of great significance to reduce the aerodynamic drag of the train.In this study,computational fluid dynamics(CFD)based on three-dimensional steady incompressible Reynolds-average Naiver-Stokes(RANS)equations and Realizable k-εturbulence model were utilized for numerical simulations.Inspired by the concept of streamlined design and the idea of bottom flow field control,this study iteratively designed the bogies in a streamlined shape and combined them with the bottom deflectors to investigate the joint drag reduction mechanism.Three models,i.e.,single-bogie model,simplified train model,and eight-car high-speed train model,were created and their aerodynamic characteristics were analyzed.The results show that the single-bogie model with streamlined design shows a noticeable drag reduction,whose power bogie and trailer bogie experience 13.92%and 7.63%drag reduction,respectively.The range of positive pressure area on the bogie is reduced.The aerodynamic drag can be further reduced to 15.01%by installing both the streamlined bogie and the deflector on the simplified train model.When the streamlined bogies and deflectors are used on the eight-car model together,the total drag reduction rate reaches 2.90%.Therefore,the proposed aerodynamic kit for the high-speed train bottom is capable to improve the flow structure around the bogie regions,reduce the bottom flow velocity,and narrow the scope of the train’s influence on the surrounding environment,achieving the appreciable reduction of aerodynamic drag.This paper can provide a new idea for the drag reduction of high-speed trains.
基金Project(NB-2020-JG-07)supported by the Research and Engineering Application of Key Technologies for New Building Industrialization Project of China Northwest Architectural Design and Research Institute Co.,Ltd.Project(2023-CXTD-29)supported by the Key Scientific and Technological Innovation Team of Shaanxi Province,ChinaProject supported by the K.C.Wong Education Foundation。
文摘The utilization of prefabricated light modular radiant heating system has demonstrated significant increases in heat transfer efficiency and energy conservation capabilities.Within prefabricated building construction,this new heating method presents an opportunity for the development of comprehensive facilities.The parameters for evaluating the effectiveness of such a system are the upper surface layer’s heat flux and temperature.In this paper,thermal resistance analysis calculation based on a simplified model for this unique radiant heating system analysis is presented with the heat transfer mechanism’s evaluation.The results obtained from thermal resistance analysis calculation and numerical simulation indicate that the thermal resistance analysis method is highly accurate with temperature discrepancies ranging from 0.44℃ to−0.44℃ and a heat flux discrepancy of less than 7.54%,which can meet the requirements of practical engineering applications,suggesting a foundation for the prefabricated radiant heating system.
基金the funding supports of the National Key Research and Development Plan,China(Grant No.2022YFC3801800)National Natural Science Foundation of China(Grant Nos.52038010 and 52078368)。
文摘In order to improve the overall resilience of the urban infrastructures, it is required to conduct blast resistant design for important building structures in the city. For complex terrain in the city, it is recommended to determine the blast load on the structures via numerical simulation. Since the mesh size of the numerical model highly depends on the explosion scenario, there is no generally applicable approach for the mesh size selection. An efficient method to determine the mesh size of the numerical model of near-ground detonation based on explosion scenarios is proposed in this study. The effect of mesh size on the propagation of blast wave under different explosive weights was studied, and the correlations between the mesh size effect and the charge weight or the scaled distance was described. Based on the principle of the finite element method and Hopkinson-Cranz scaling law, a mesh size measurement unit related to the explosive weight was proposed as the criterion for determining the mesh size in the numerical simulation. Finally, the applicability of the method proposed in this paper was verified by comparing the results from numerical simulation and the explosion tests and was verified in AUTODYN.
文摘The reverse magnetohydrodynamic(MHD)energy bypass technology is a promising energy redis⁃tribution technology in the scramjet system,in augmented with a power generation equipment to supply the neces⁃sary long-distance flight airframe power.In this paper,a computational model of the scramjet magnetohydrody⁃namic channel is developed and verified by using the commercial software Fluent.It is found that when the mag⁃netic induction intensity is 1,2,3,4 T,the power generation efficiency is 22.5%,22.3%,22.0%,21.5%,and decreases with the increase of the magnetic induction intensity,and the enthalpy extraction rate is 0.026%,0.1%,0.21%,0.34%,and increases with the increase of the magnetic induction intensity.The deceleration ef⁃fect of electromagnetic action on the airflow in the power channel increases with the increase of magnetic induc⁃tion intensity.The stronger the magnetic field intensity,the more obvious the decreasing effect of fluid Mach num⁃ber in the channel.The power generation efficiency decreases as the magnetic induction intensity increases and the enthalpy extraction rate is reversed.As the local currents gathering at inlet and outlet of the power generation area,total temperature and enthalpy along the flow direction do not vary linearly,and there are maximum and minimum values at inlet and outlet.Increasing the number of electrodes can effectively regulate the percentage of Joule heat dissipation,which can improve the power generation efficiency.
基金This research has been conducted under SEGTRANS project,funded by the Centre for Industrial Technological Development(CDTI,Government of Spain).
文摘This paper explores the performances of a finite element simulation including four concrete models applied to a full-scale reinforced concrete beam subjected to blast loading. Field test data has been used to compare model results for each case. The numerical modelling has been, carried out using the suitable code LS-DYNA. This code integrates blast load routine(CONWEP) for the explosive description and four different material models for the concrete including: Karagozian & Case Concrete, Winfrith, Continuous Surface Cap Model and Riedel-Hiermaier-Thoma models, with concrete meshing based on 10, 15, and 20 mm. Six full-scale beams were tested: four of them used for the initial calibration of the numerical model and two more tests at lower scaled distances. For calibration, field data obtained employing pressure and accelerometers transducers were compared with the results derived from the numerical simulation. Damage surfaces and the shape of rupture in the beams have been used as references for comparison. Influence of the meshing on accelerations has been put in evidence and for some models the shape and size of the damage in the beams produced maximum differences around 15%. In all cases, the variations between material and mesh models are shown and discussed.
基金the Natural Science Foundation of Jiangsu Province(Grant No.BK20200494)China Postdoctoral Science Foundation(Grant No.2021M701725)+3 种基金Jiangsu Postdoctoral Research Funding Program(Grant No.2021K522C)Fundamental Research Funds for the Central Universities(Grant No.30919011246)National Natural Science Foundation of China(Grant No.52278188)Natural Science Foundation of Jiangsu Province(Grant No.BK20211196)。
文摘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.
文摘During the propagation of high-power lasers within internal channels,the laser beam heats the propagation medium,causing the thermal blooming effect that degrades the beam quality at the output.The intricate configuration of the optical path within the internal channel necessitates complex and time-consuming efforts to assess the impact of thermal blooming effect on the optical path.To meet the engineering need for rapid evaluation of thermal blooming effect in optical paths,this study proposed a rapid simulation method for the thermal blooming effect in internal optical paths based on the finite element method.This method discretized the fluid region into infinitesimal elements and employed finite element method for flow field analysis.A simplified analytical model of the flow field region in complex internal channels was established,and regions with similar thermal blooming effect were divided within this model.Based on the calculated optical path differences within these regions,numerical simulations of phase distortion caused by thermal blooming were conducted.The calculated result were compared with those obtained using the existing methods.The findings reveal that for complex optical paths,the discrepancy between the two approaches is less than 3.6%,with similar phase distortion patterns observed.For L-type units,this method and the existing methods identify the same primary factors influencing aberrations and exhibit consistent trends in their variation.This method was used to analyze the impact of thermal blooming effect in a straight channel under different gravity directions.The results show that phase distortion varies with changes in the direction of gravity,and the magnitude of the phase difference is strongly correlated with the component of gravity perpendicular to the optical axis.Compared to the existing methods,this approach offers greater flexibility,obviates the need for complex custom analysis programming.The analytical results of this method enable a rapid assessment of the thermal blooming effect in optical paths within the internal channel.This is especially useful during the engineering design.These results also provide crucial references for developing strategies to suppress thermal blooming effect.
文摘[Objective]This study aims to develop a thermodynamically consistent phase-field framework for modeling the initiation and evolution of discontinuous structures in geomaterials.[Methods]Our model introduces crack driving forces derived from the volumetric-deviatoric strain decomposition strategy,incorporating distinct tension,compression,and shear degradation mechanisms.Inertia effects capture compaction-band formation driven by wave-like disturbances,grain crushing,and frictional rearrangement.A monolithic algorithm ensures numerical stability and rapid convergence.[Results]The framework reproduces tensile,shear,mixed tensile-shear,and compressive-shear failures using the Benzeggagh-Kenane criterion.Validation against benchmark simulations-including uniaxial compression of rock-like and triaxial compression of V-notched sandstone specimens-demonstrates accurate predictions of crack initiation stress,localization orientation,and energy dissipation.[Conclusions]The framework provides a unified and robust numerical tool for analyzing the spatiotemporal evolution of strain localization and fracture in geomaterials.[Significance]By linking microscale fracture dynamics with macroscale failure within a thermodynamically consistent scheme,this study advances predictive modeling of rock stability,slope failure,and subsurface energy systems,contributing to safer and more sustainable geotechnical practice.
文摘This study focuses on the thermal management of 4680-type cylindrical lithium-ion battery packs utilizing NCM811 chemistry.It establishes coupled multi-physics models for both immersion and serpentine cold plate cooling systems.Through a combination of numerical simulation and experimental validation,the technical advantages and mechanisms of immersion cooling are systematically explored.Simulation results indicate that under a 3C fast-charging condition(inlet temperature 20℃,flow rate 36 L/min),the immersion cooling structure 3demonstrates a triple enhancement in thermal performance compared to the cold plate structure 1:a 13.06%reduction in peak temperature,a 31.67%decrease in overall maximum temperature difference,and a 47.62%decrease in single-cell temperature deviation,while also reducing flow resistance by 33.61%.Furthermore,based on the immersion cooling model,a small battery module comprising seven cylindrical cells was designed for thermal runaway testing via nail penetration.The results show that the peak temperature of the triggered cell was limited to 437.6℃,with a controllable temperature rise gradient of only 3.35℃/s and a rapid cooling rate of 0.6℃/s.The maximum temperature rise of adjacent cells was just 64.8℃,effectively inhibiting thermal propagation.Post-test disassembly revealed that the non-triggered cells retained>99.2%of their original voltage and>99%structural integrity,confirming the module’s ability to achieve“localized failure with global stability.”
文摘The finite volume method was applied to numerically simulate the bottom pressure field induced by regular waves,vehicles in calm water and vehicles in regular waves.The solution of Navier-Stokes(N-S)equations in the vicinity of numerical wave tank's boundary was forced towards the wave theoretical solution by incorporating momentum source terms,thereby reducing adverse effects such as wave reflection.Simulations utilizing laminar flow,turbulent flow,and ideal fluid models were all found capable of effectively capturing the waveform and bottom pressure of regular waves,agreeing well with experimental data.In predicting the bottom pressure field of the submerged vehicle,turbulent simulations considering fluid viscosity and boundary layer development provided more accurate predictions for the stern region than inviscid simulations.Due to sphere's diffractive effect,the sphere's bottom pressure field in waves is not a linear superposition of the wave's and the sphere's bottom pressure field.However,a slender submerged vehicle exhibits a weaker diffractive effect on waves,thus the submerged vehicle's bottom pressure field in waves can be approximated as a linear superposition of the wave's and the submerged vehicle's bottom pressure field,which simplifies computation and analysis.
基金Project(2024ZD1003704)supported by the Deep Earth Probe and Mineral Resources Exploration-National Science and Technology Major Project,ChinaProjects(51834001,52130404)supported by the National Natural Science Foundation of China。
文摘During upward horizontal stratified backfill mining,stable backfill is essential for cap and sill pillar recovery.Currently,the primary method for calculating the required strength of backfill is the generalized three-dimensional(3 D)vertical stress model,which ignores the effect of mine depth,failing to obtain the vertical stress at different positions along stope length.Therefore,this paper develops and validates an improved 3 D model solution through numerical simulation in Rhino-FLAC^(3D),and examines the stress state and stability of backfill under different conditions.The results show that the improved model can accurately calculate the vertical stress at different mine depths and positions along stope length.The error rates between the results of the improved model and numerical simulation are below 4%,indicating high reliability and applicability.The maximum vertical stress(σ_(zz,max))in backfill is positively correlated with the degree of rock-backfill closure,which is enhanced by mine depth and elastic modulus of backfill,while weakened by stope width and inclination,backfill friction angle,and elastic modulus of rock mass.Theσ_(zz,max)reaches its peak when the stope length is 150 m,whileσ_(zz,max)is insensitive to changes in rock-backfill interface parameters.In all cases,the backfill stability can be improved by reducingσ_(zz,max).The results provide theoretical guidance for the backfill strength design and the safe and efficient recovery of ore pillars in deep mining.
文摘Whipple shields as sacrificial bumpers,safeguard the satellites against extremely fast,different-sized projectiles traveling through space in the low earth orbit.Typical Whipple shields comprise a front and rear plate,separated by a gap or space.Recent advancements have explored the use of foam,cellular cores,and alternative materials such as ceramics instead of aluminium for the plates.In the current work,the effect of including fluid cores(air/water)sandwiched between the front and rear plates,on the response to hypervelocity impact was explored through a numerical approach.The numerical simulation consisted of hypervelocity impact by a 2 mm diameter,stainless steel projectile,launched at speeds of 3 e9 km/s with a normal impact trajectory towards the Whipple shield.The front and rear bumpers,made of AA6061-T6,were each 1 mm thick.A space of 10 mm was taken between the plates(occupied by fluid).The key metrics analyzed were the perforation characteristics,stages of the debris cloud generation and propagation,energy variations(internal,kinetic and plastic work),temperature variations,and the fragmentation summary.From the computational analysis,employing water-core in Whipple shields could prevent the rear bumper perforation till 6 km/s,lower the peak temperatures at the front bumper perforation zones and debris tip,and generate fewer,larger fragments.
基金the support provided by the Technology Innovation Project (Grant No. KYGYZB002201) for the research work
文摘Investigating the blast effects and mechanisms on typical finite-sized obstacles is essential for optimizing defense strategies and designing more robust barriers to deter terrorists and protect critical locations.This study investigates the blasting effects and underlying mechanisms of concrete frustums subjected to contact explosions,employing both numerical simulations and field tests.It focuses on the effects of top and side blasting,with particular emphasis on fracture modes,damage patterns,and fragment sizes,as well as the causes of different failure modes and the propagation of stress waves.The study also explores the blasting effects of detonating explosives at varying positions along the side and with different charge amounts.The results show that side-blasting leads to complete fragmentation,with tensile waves playing a significant role in creating extensive damage zones that propagate parallel to the frustum's outer surface,concentrating damage near the surface.During top-blasting,the upper half of the frustum undergoes fragmentation,while the lower half experiences cracking.Tensile waves propagate from the top to the bottom surface,forming larger blocks in regions with lower wave intensity.Three distinct damage zones within the frustum were identified,and a series of mathematical formulas were derived to describe the relationship between the maximum fragment size and charge mass.As the charge mass increased from 1.0 kg to 4.0 kg,the maximum fragment size decreased.Detonation at the center of the frustum's side resulted in the most severe fragmentation,with a 51.8%reduction in fragment size compared to other detonation positions.Finally,four broken modes were classified,each influenced by charge mass and explosive location.This study provides valuable insights for optimizing civil blasting operations and designing protective engineering structures.
基金funded by the National Natural Science Foundation of China(Grant No. 12302437)Jiangsu Provincial Natural Science Foundation (Grant No.SBK2023045424)。
文摘To explore the design criteria for composite charges and reveal the intrinsic relationship between the detonation wave propagation in composite charges and the overall energy distribution of shock waves,this study analyzes the propagation and interaction processes of detonation waves in composite charges with different structural dimensions and explosive combinations. It also investigates the spatial distribution characteristics of the resulting shock wave loads. Based on dimensional analysis theory, a theoretical analysis of the shock wave overpressure distribution in free air fields is conducted. Utilizing the derived dimensionless function relationships, the hydrocode AUTODYN is employed to investigate the effects of charge structure parameters and explosive combinations on the internal overdriven detonation phenomena and the distribution of shock wave loads. It is found that the overdriven detonation phenomenon in the inner layer of composite charges increases the strength of the axial detonation wave,thereby enhancing the intensity of the primary end wave formed upon refraction into the air, which affects the distribution characteristics of the shock wave overpressure. Research has shown that increasing the thickness ratio and detonation velocity ratio of composite charges is beneficial for exacerbating the phenomenon of overdriven detonation, improving the primary end wave intensity and axial overpressure. This gain effect gradually weakens with the propagation of shock waves. When overdriven detonation occurs inside the composite charge, the detonation pressure first increases and then decreases. The Mach reflection pressure of the composite charge with a larger aspect ratio is attenuated to a greater extent. In addition, as the aspect ratio of the composite charge increases, the shock wave energy gradually flows from the axial direction to the radial direction. Therefore, as the aspect ratio of the composite charge increases, the primary end wave intensity and axial overpressure gradually decrease.
基金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.
