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.展开更多
Ice load on underwater vehicles breaking through ice covers from underneath is a significant concern for researchers in polar exploration,and the research on this problem is still in its early stages.Both mechanical e...Ice load on underwater vehicles breaking through ice covers from underneath is a significant concern for researchers in polar exploration,and the research on this problem is still in its early stages.Both mechanical experimental measurement and numerical simulation pose research challenges.This study focuses on the ice load of a cylinder structure breaking upward through the ice sheet form underneath in the Small Ice Model Basin of China Ship Scientific Research Center(CSSRC SIMB).A high-speed camera system was employed to observe the ice sheet failure during the tests,in which,with the loading position as center,local radial cracks and circumferential cracks were generated.A load sensor was used to measure the overall ice load during this process.Meanwhile,a numerical model was developed using LS-DYNA for validation and comparison.With this model,numerical simulation was conducted under various ice thicknesses and upgoing speeds to analyze the instantaneous curves of ice load.The calculation results were statistically analyzed under different working conditions to determine the influence of the factors on the ice load of the cylinder.The study explores the measurement method about ice load of objects vertically breaking through model ice sheet and is expected to provide some fundamental insights into the safety design of underwater structures operating in ice waters.展开更多
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.”展开更多
[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.展开更多
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 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.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
The intersection is a widely used traffic line structure from the shallow tunnel to the deep roadway,and determining the subsidence hidden danger area of the roof is the key to its stability control.However,applying t...The intersection is a widely used traffic line structure from the shallow tunnel to the deep roadway,and determining the subsidence hidden danger area of the roof is the key to its stability control.However,applying traditional maximum equivalent span beam(MESB)theory to determine deformation range,peak point,and angle influence poses a challenge.Considering the overall structure of the intersection roof,the maximum equivalent triangular plate(METP)theory is proposed,and its geometric parameter calculation formula and deflection calculation formula are obtained.The application of the two theories in 18 models with different intersection angles,roadway types,and surrounding rock lithology is verified by numerical analysis.The results show that:1)The METP structure of the intersection roof established by the simulation results of each model successfully determined the location of the roof’s high displacement zone;2)The area comparison method of the METP theory can be reasonably explained:①The roof subsidence of the intersection decreases with the increase of the intersection angle;②The roof subsidence at the intersection of different roadway types has a rectangular type>arch type>circular type;③The roof subsidence of the intersection with weak surrounding rock is significantly larger than that of the intersection with hard surrounding rock.According to the application results of the two theories,the four advantages of the METP theory are compared and clarified in the basic assumptions,mechanical models,main viewpoints,and mechanism analysis.The large deformation inducement of the intersection roof is then explored.The J 2 peak area of the roof drives the large deformation of the area,the peak point of which is consistent with the center of gravity position of the METP.Furthermore,the change in the range of this peak is consistent with the change law of the METP’s area.Hence,this theory clarifies the large deformation area of the intersection roof,which provides a clear guiding basis for its initial support design,mid-term monitoring,and late local reinforcement.展开更多
To enhance the recuperation rate of the mine and comply with the stipulations of green mining technology, it is vital to expeditiously recuperate the coal pillar resources in the final stage, thus preventing the consi...To enhance the recuperation rate of the mine and comply with the stipulations of green mining technology, it is vital to expeditiously recuperate the coal pillar resources in the final stage, thus preventing the considerable squandering of resources. The coal pillar resource of the main roadway and its branch roadway constitutes a significant recovery subject. Its coal pillar shape is regular and possesses a considerable strike distance, facilitating the arrangement of the coal pillar recovery working face (CPRWF) for mining operations. However, for the remaining coal pillars with a thick and hard roof (THF) and multiple tectonic zones, CPRWF encounters challenges in selecting an appropriate layout, managing excessive roof pressure, and predicting mining stress. Aiming at the roadway coal pillar group with THF and multi-structural areas in specific projects, a method of constructing multi-stage CPRWF by one side gob-side entry driving (GSED) and one side roadway reusing is proposed. Through theoretical calculation of roof fracture and numerical simulation verification, combined with field engineering experience and economic analysis, the width of the narrow coal pillar (NCP) in the GSED is determined to be 10 m and the length of the CPRWF is 65 m. Concurrently, the potential safety hazard that the roof will fall asymmetrically and THF is difficult to break during CPRWF mining after GSED is analyzed and verified. Then, a control method involving the pre-cutting of the roof in the reused roadway before mining is proposed. This method has been shown to facilitate the complete collapse of THF, reduce the degree of mine pressure, and facilitate the symmetrical breaking of the roof. Accordingly, a roof-cutting scheme based on a directional drilling rig, bidirectional shaped polyvinyl chloride (PVC) pipe, and emulsion explosive was devised, and the pre-splitting of 8.2 m THF was accomplished. Field observations indicate that directional cracks are evident in the roof, the coal wall is flat during CPRWF mining, and the overall level of mining pressure is within the control range. Therefore, the combined application of GSED and roof-cutting technology for coal pillar recovery has been successfully implemented, thereby providing new insights and engineering references for the construction and pressure relief mining of CPRWF.展开更多
基金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.
文摘Ice load on underwater vehicles breaking through ice covers from underneath is a significant concern for researchers in polar exploration,and the research on this problem is still in its early stages.Both mechanical experimental measurement and numerical simulation pose research challenges.This study focuses on the ice load of a cylinder structure breaking upward through the ice sheet form underneath in the Small Ice Model Basin of China Ship Scientific Research Center(CSSRC SIMB).A high-speed camera system was employed to observe the ice sheet failure during the tests,in which,with the loading position as center,local radial cracks and circumferential cracks were generated.A load sensor was used to measure the overall ice load during this process.Meanwhile,a numerical model was developed using LS-DYNA for validation and comparison.With this model,numerical simulation was conducted under various ice thicknesses and upgoing speeds to analyze the instantaneous curves of ice load.The calculation results were statistically analyzed under different working conditions to determine the influence of the factors on the ice load of the cylinder.The study explores the measurement method about ice load of objects vertically breaking through model ice sheet and is expected to provide some fundamental insights into the safety design of underwater structures operating in ice waters.
文摘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.”
文摘[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.
文摘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 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.
基金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.
基金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.
文摘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.
基金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.
基金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.
文摘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.
基金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.
基金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.
文摘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(52204164)supported by the National Natural Science Foundation of ChinaProject(2021QNRC001)supported by the Young Elite Scientists Sponsorship Program by CAST,China。
文摘The intersection is a widely used traffic line structure from the shallow tunnel to the deep roadway,and determining the subsidence hidden danger area of the roof is the key to its stability control.However,applying traditional maximum equivalent span beam(MESB)theory to determine deformation range,peak point,and angle influence poses a challenge.Considering the overall structure of the intersection roof,the maximum equivalent triangular plate(METP)theory is proposed,and its geometric parameter calculation formula and deflection calculation formula are obtained.The application of the two theories in 18 models with different intersection angles,roadway types,and surrounding rock lithology is verified by numerical analysis.The results show that:1)The METP structure of the intersection roof established by the simulation results of each model successfully determined the location of the roof’s high displacement zone;2)The area comparison method of the METP theory can be reasonably explained:①The roof subsidence of the intersection decreases with the increase of the intersection angle;②The roof subsidence at the intersection of different roadway types has a rectangular type>arch type>circular type;③The roof subsidence of the intersection with weak surrounding rock is significantly larger than that of the intersection with hard surrounding rock.According to the application results of the two theories,the four advantages of the METP theory are compared and clarified in the basic assumptions,mechanical models,main viewpoints,and mechanism analysis.The large deformation inducement of the intersection roof is then explored.The J 2 peak area of the roof drives the large deformation of the area,the peak point of which is consistent with the center of gravity position of the METP.Furthermore,the change in the range of this peak is consistent with the change law of the METP’s area.Hence,this theory clarifies the large deformation area of the intersection roof,which provides a clear guiding basis for its initial support design,mid-term monitoring,and late local reinforcement.
基金Project(52204164) supported by the National Natural Science Foundation of ChinaProject(2023ZKPYSB01) supported by the Fundamental Research Funds for the Central Universities,China。
文摘To enhance the recuperation rate of the mine and comply with the stipulations of green mining technology, it is vital to expeditiously recuperate the coal pillar resources in the final stage, thus preventing the considerable squandering of resources. The coal pillar resource of the main roadway and its branch roadway constitutes a significant recovery subject. Its coal pillar shape is regular and possesses a considerable strike distance, facilitating the arrangement of the coal pillar recovery working face (CPRWF) for mining operations. However, for the remaining coal pillars with a thick and hard roof (THF) and multiple tectonic zones, CPRWF encounters challenges in selecting an appropriate layout, managing excessive roof pressure, and predicting mining stress. Aiming at the roadway coal pillar group with THF and multi-structural areas in specific projects, a method of constructing multi-stage CPRWF by one side gob-side entry driving (GSED) and one side roadway reusing is proposed. Through theoretical calculation of roof fracture and numerical simulation verification, combined with field engineering experience and economic analysis, the width of the narrow coal pillar (NCP) in the GSED is determined to be 10 m and the length of the CPRWF is 65 m. Concurrently, the potential safety hazard that the roof will fall asymmetrically and THF is difficult to break during CPRWF mining after GSED is analyzed and verified. Then, a control method involving the pre-cutting of the roof in the reused roadway before mining is proposed. This method has been shown to facilitate the complete collapse of THF, reduce the degree of mine pressure, and facilitate the symmetrical breaking of the roof. Accordingly, a roof-cutting scheme based on a directional drilling rig, bidirectional shaped polyvinyl chloride (PVC) pipe, and emulsion explosive was devised, and the pre-splitting of 8.2 m THF was accomplished. Field observations indicate that directional cracks are evident in the roof, the coal wall is flat during CPRWF mining, and the overall level of mining pressure is within the control range. Therefore, the combined application of GSED and roof-cutting technology for coal pillar recovery has been successfully implemented, thereby providing new insights and engineering references for the construction and pressure relief mining of CPRWF.
