Recently Accepted

Recently Accepted articles have been peer-reviewed and accepted, which are not yet assigned to volumes/issues, but are citable by Digital Object Identifier (DOI).
Display Method:
Reaction Evolution Characteristics of Ignited DNAN-Based Explosive Charges with Pre-cracks
YAO Xin, WANG Hui, SHEN Fei, QU Kepeng
, Available online  , doi: 10.11858/gywlxb.20251181
Abstract:
To investigate the influence of cracks and gaps on the reaction evolution characteristics of aluminum-containing DNAN-based explosives after the formation of mechanical induced hotspots, explosive charge samples with different initial cracks were fabricated. An explosive impact ignition device based on gun propellant combustion loading was designed. The evolution process following the ignition of explosives was simulated. Pressure changes and post-test morphological features of the explosives were recorded. Numerical simulations were conducted to analyze the stress field and reaction distribution of explosive charges with different initial cracks under the same loading conditions. The results indicate that the crack-free and single-line crack explosive charge with no gap debris remained intact, and pressure dropped rapidly after the peak with no reaction occurred, and the hot spot region was located at the bottom. While for the single-line crack explosive charge with 1mm gap, the explosive charge fractured and exhibited local low-order reactions, with a slow pressure decay process. Among these, the hot spot region of the single-line crack explosive shifted to the side surface, while the cross-line crack explosive formed dual hot spot regions on both the side surface and bottom, further enhancing the reaction intensity. This demonstrates that pre-cracks significantly influence the explosive reaction process by altering stress distribution and expanding hot spot regions.
Phase Transformation, Sintering Mechanism and Dynamics of Singlet-Doublet Al Nanosphere Collisions with Initial I-Shaped Configuration
JIANG Jun, SUN Weifu
, Available online  , doi: 10.11858/gywlxb.20251176
Abstract:
Molecular dynamics simulations are used to study the dynamics of a single Al nanosphere (singlet) colliding with an aggregate of two Al nanospheres (doublet) with initial I-shaped configuration. Depending on the initial impact velocity, there are four collision outcomes, namely bounce, adhesion, aggregation and melting. At a very low velocity, the repulsive force between the nanospheres will cause the nanospheres to rebound without contact, and the critical velocity of bounce decreases with the increase of the diameter of the nanosphere. As the velocity increases, the nanospheres are sintered together due to adhesion between them and the formation of new bonds. The phase transformation and atomic diffusion during singlet-doublet collisions are quantitatively characterized by Common Neighbor Analysis, Dislocation Analysis and mean square displacement to explore the underlying sintering mechanism. The critical impact velocity of singlet melting is obtained by monitoring the temperature of singlet with different diameters.
Translational-Rotational Decoupling Dynamics of High-Pressure Liquid Water under Quasi-Isentropic Compression
DENG Changhao, CHEN Bo, DAI Jiayu
, Available online  , doi: 10.11858/gywlxb.20251222
Abstract:
Water is widely present on Earth, as well as in ocean planets, ice giants, and interstellar space. Its phase states and properties across a broad thermodynamic range serve as the core foundation for addressing numerous key scientific issues, including biochemical reactions, environmental climate dynamics, and planetary internal structure evolution. Although liquid water exhibits abundant anomalous behaviors under extreme pressure conditions, research in this field has long been constrained by bottlenecks in experimental measurement techniques and the complexity of theoretical calculations. As a result, atomic-scale data available for analysis remain extremely scarce, which limits in-depth understanding of its underlying microscopic mechanisms. In this work, a deep learning interaction model was constructed based on high-precision ab initio data. Using molecular dynamics simulations, liquid water was compressed along an isentrope to pressures on the order of tens of thousands of atmospheres. The structural and dynamic properties of liquid water under different pressure conditions were systematically calculated and analyzed. The results show that the inherent tetrahedral local coordination environment of water molecules in liquid water is significantly disrupted under high pressure, leading to a marked enhancement in the rotational motion capability of water molecules as pressure increases. In sharp contrast, the translational motion of water molecules is strongly constrained in the highly condensed high-pressure environment. Furthermore, the mean squared displacement (MSD) of water molecules under high pressure exhibits a typical three-stage behavior characteristic of glassy systems, namely the ballistic transport region, the plateau region, and the diffusion region. From a macroscopic perspective, the significant reduction in translational motion capability is manifested as a substantial increase in shear viscosity. Of particular importance is that, unlike supercooled water under ambient pressure—where translational and rotational motions are strongly coupled—liquid water under dynamic high pressure exhibits an intrinsic decoupling phenomenon between translational and rotational motions. The findings of this work are expected to provide meaningful microscopic insights for investigating important scientific issues such as the response of materials under dynamic loading and the solidification of metastable liquids.
Energy Conversion Prediction Model of Expansion Tube under Near-Field Explosion Loading
QI Zizhen, LI Minghao, ZHANG Yuyan, LIANG Minzu, ZHANG Yuwu, LIN Yuliang
, Available online  , doi: 10.11858/gywlxb.20251227
Abstract:
The near-field region of an explosion is the core zone of munition-induced damage, characterized by the coupled loading of intense shock waves and detonation products. Currently, the mechanical response and energy conversion mechanisms of Expansion Tube Structures (ETS) under such extreme loading conditions remain unclear. In this study, ETS is adopted as a representative energy-absorbing structure to investigate its energy conversion behavior under the coupled action of near-field shock waves and detonation products. Based on experimental validation, numerical simulations are employed to analyze the characteristics of near-field blast loads and the dynamic response of ETS. Furthermore, a theoretical prediction formula for near-field blast loads is established, and a theoretical model for predicting energy conversion efficiency is developed under the strong-shock assumption. The results show that the energy conversion efficiency decreases significantly with increasing scaled distance, dropping below 10% when the scaled distance exceeds 0.8 m/kg¹/³. Moreover, the energy conversion efficiency exhibits a strong positive correlation with the specific impulse of the reflected wave, indicating that specific impulse is the key factor governing energy transfer. This work elucidates the intrinsic energy conversion mechanism of ETS under near-field coupled loading, and the proposed theoretical model provides a robust foundation for the design and performance evaluation of near-field protective structures.
Research Progress on Two-Dimensional Diamond
MING Jiaxin, LI Jiayin, CHEN Yabin
, Available online  , doi: 10.11858/gywlxb.20251248
Abstract:
Two-dimensional (2D) diamond, an atomically thin carbon-based material, not only inherits the exceptional properties of bulk diamond but is also expected to exhibit unique physical characteristics arising from nanoscale effects. Currently, research on 2D diamond remains in its infancy, being primarily driven by theoretical investigations, while experimental efforts have mainly focused on its controllable synthesis and structural characterization. Owing to pronounced interfacial effects, the direct application of conventional high-pressure synthesis methods to nanoscale systems is considerably limited, making it challenging to achieve a stable transition from sp2 to sp3 hybridization, thereby posing numerous critical scientific challenges for the study of 2D diamond. This review systematically summarizes recent theoretical and experimental advances in the structural features, synthesis strategies, and physicochemical properties of 2D diamond, and provides perspectives on future research directions and scientific opportunities in the field of 2D diamond.
Calculate the Viscosity of Iron-Sulfur Alloy under the Conditions of the Earth Inner Core Based on the Neural Network Potential
XU Yunfan, HE Yu, ZHANG Wei, LI Heping
, Available online  , doi: 10.11858/gywlxb.20251221
Abstract:
Viscosity is a key physical property that governs the dynamics and evolutionary history of the Earth inner core and plays an important role in the origin of seismic anisotropy. Previous studies have investigated the viscosity of pure iron in the hexagonal close-packed (HCP) and body-centered cubic (BCC) phases under inner-core conditions through computational simulations. However, the inner core also contains light elements such as carbon, hydrogen, oxygen, sulfur, and silicon, and the effect of these light elements on the viscosity of the inner core remains insufficiently understood. In this study, we constructed a neural-network potential (NNP) for Fe-S alloy under inner core conditions and employed it to perform large-scale molecular dynamics simulations. We systematically examined the effect of vacancy concentrations as low as 0.02% on the ionic transport properties of Fe-S alloy. Based on the self-diffusion coefficients of Fe atoms in the lattice, we further explored the creep mechanisms and viscosity of Fe-S alloys under core conditions. The results indicate that dislocation creep dominates the rheological behavior, yielding viscosities in the range of 1×1014-2×1016 Pa·s, which are consistent with constraints from free-core nutation and seismic observations.
Polyvinylidene Fluoride (PVDF) Phase Transitions under in-Situ High Pressure
ZHANG Shenghan, LI Ting, ZHANG Xiaojun, CHEN Zhiqiang
, Available online  , doi: 10.11858/gywlxb.20251174
Abstract:
Poly(vinylidene fluoride) (PVDF) is a versatile semi-crystalline polymer with outstanding piezoelectric, pyroelectric, and dielectric properties, widely used in sensors, energy devices, and biomedical applications. Its performance is governed by crystallinity and polymorphic structure (α, β, γ, δ, ε), among which the polar β phase possesses superior electromechanical characteristics compared with the non-polar α phase. Nevertheless, the α phase remains thermodynamically the most stable and the most readily obtained. In this work, the structural evolution of PVDF under high pressure was investigated by means of in situ X-ray diffraction and Fourier-transform infrared spectroscopy. At ambient conditions the powder consists primarily of the α phase with a minor fraction of α. Upon compression to 0–20 GPa, the α phase gradually diminishes; the emergence of new diffraction peaks and band shifts indicates sequential α → β and β → γ transformations, accompanied by a pronounced increase in β content and concomitant formation of γ. When the pressure exceeds 20 GPa, severe lattice distortion destroys long-range crystalline order, resulting in peak broadening and eventual amorphization. The study unveils the intricate interplay between pressure-induced chain rearrangement and polymorphic transitions, clarifies the high-pressure phase-transformation pathway and structural evolution of PVDF, and thereby deepens the structure–property understanding of this polymer. The findings also provide a theoretical basis for tailoring its performance under extreme conditions and for designing high-pressure technologies.
Theoretical Study on Structural Stability and Superionic Phase Transition of UH5 under High Pressure
DING Yuqing, JIA Xixi, ZHANG Wenhui, WANG Hui
, Available online  , doi: 10.11858/gywlxb.20251224
Abstract:
The thermodynamic, mechanical, and dynamical stability, along with the electronic properties of UH5 within 30 GPa, are systematically investigated using first-principles calculations. The experimentally synthesized orthorhombic, hexagonal, and cubic phases are all found to be magnetic materials, with spin polarizations of 82%, 100%, and 100%, respectively, and their thermodynamic stability decreases sequentially. Elastic constant and phonon calculations demonstrate that all three phases are mechanically and dynamically stable. Chemical bonding analysis indicates that this stability primarily originates from the prevalent covalent U-H interaction within the lattice. Furthermore, it is predicted that the orthorhombic phase, which has been experimentally quenched to 1 GPa, transforms into a superionic state at 1200 K, where hydrogen ions undergo rapid diffusion within the uranium sublattice interstices, achieving a diffusion coefficient of 1.2 × 10 -4 <italic>cm²/s.</italic>
FENG Yuheng, LIANG Anqi, LIU Xingyu, YIN Jianping, YI Jianya, ZHANG Xuepeng
, Available online  , doi: 10.11858/gywlxb.20251213
Abstract:
To investigate the influence of torpedo guidance nose configuration on the lethality of an underwater shaped charge warhead, a series of numerical simulations were performed using the AUTODYN finite element code. The damage performance of the shaped penetrator under different simulated nose structures was studied, analyzing the complete process including shock wave diffraction, behind-target load propagation, and target damage. The results indicate that both the penetrator head velocity and the behind-target hole diameter generally increase with the total length and number of layers of the simulated nose. Within a certain range, increasing the number of nose layers effectively optimizes the formation of the EFP, thereby enhancing its penetration capability. Furthermore, there exists an optimal total nose length that maximizes the head velocity while preventing necking and fracture of the penetrator.
Quantification of Uncertainty in Magnesium Oxide and Rhenium Pressure Standards Based on Bayesian Statistical Methods
DAI Feifan, XIANG Shikai, LI Weiwei, ZHANG Ruizhi, ZHANG Jian, LUO Guoqiang, WU Run, XIAN Yunting
, Available online  , doi: 10.11858/gywlxb.20251056
Abstract:

Accurate pressure measurement in static high-pressure experiments relies on the equation of state (EOS) of standard materials, where uncertainties in EOS parameters can significantly affect the accuracy of pressure predictions. This study focuses on magnesium oxide (MgO, B1 phase) and rhenium (Re, hexagonal close packed phase), employing Bayesian statistical methods and Markov Chain Monte Carlo (MCMC) simulation techniques to systematically quantify the uncertainty in pressure prediction during diamond anvil cell (DAC) experiments. By constructing a Bayesian framework with uniform prior distributions and normal likelihood functions, and integrating multiple sets of experimental data for parameter calibration, the results demonstrate that the Bayesian statistical approach successfully quantifies the posterior distribution of EOS parameters, revealing strong correlations between them, e.g., a negative correlation between Grüneisen parameter and initial volume for MgO, and a positive correlation between bulk modulus and Grüneisen parameter for Re. The uncertainty in pressure predictions for both MgO and Re increases significantly at higher pressures; for Re, this uncertainty also rises markedly with increasing temperature, whereas no clear trend is observed for MgO. This study provides pressure benchmarks with quantified uncertainties, contributing to improved accuracy in high-pressure experimental measurements. It holds significant reference value for ensuring the reliability of experimental data in materials science and geophysical research.

First-Principles Study on the Multiphase Equation of State of Tin
CHEN Kaile, WANG Yuechao, XU Yuanji, LIU Yu, XIAN Jiawei, WANG Lifang, JIAN Dan, LIU Haifeng, SONG Haifeng
, Available online  , doi: 10.11858/gywlxb.20251054
Abstract:

Metallic tin is a focal point in high-pressure physics research and a critical material of strategic importance in defense-related technologies. Due to the rich physical phases of tin, it is crucial to study the multiphase equation of state and phase boundaries of tin, whether in basic research or industrial applications. This work systematically investigates the high-temperature and high-pressure multiphase equation of state (EOS), phase boundaries, elastic modulus, sound velocities, and Hugoniot curves of tin using density functional theory (DFT) combined with the mean-field potential (MFP) method. The results not only provide the multiphase EOS of tin under extreme conditions but also demonstrate good agreement with experimental data for the β-γ phase boundary and ambient-pressure sound velocities of β-Sn. Furthermore, this study evaluates the effects of different density functionals (LDA, PBEsol, and SCAN) on the high-pressure EOS. The LDA and PBEsol functionals show superior consistency with experimental Hugoniot curves and ambient-pressure elastic moduli, while the SCAN functional exhibits larger deviations in phase boundary predictions but achieves closer agreement with experimental ambient-pressure sound velocities for β-Sn.

Simulation Study on Hard X-Ray Detection Efficiency for Microchannel Plate
YANG Jing, DAN Lianqiang
, Available online  , doi: 10.11858/gywlxb.20251193
Abstract:
An improved detection efficiency model for microchannel plate response to hard X-ray is described, which builds on previous models by incorporating a more detailed consideration of cross talking for photoelectron cross-section between atomic shells in the MCP bulk. An analytical investigation and numerical calculation of the detection efficiency were carried out as function of compositional parameters, channel diameter, wall thickness, MCP thickness, respectively. Furthermore, according to the calculation results and developed technologies, a group of optimized parameters were given out, and 45% detection efficiency was found response to 50 ~ 200 keV X-ray.
Effect of Boron Nitridecontent on the Explosion Performance of On-Site Mixed Emulsion Explosives
FU Jiakun, LIU Feng, ZHU Zhengde, CHEN Chuanbin
, Available online  , doi: 10.11858/gywlxb.20251223
Abstract:
To investigate the effect of boron nitride (BN) content on the explosive performance of on-site mixed bulk emulsion explosives, transmission electron microscopy (TEM) and optical microscopy were employed to characterize the microstructure of BN particles and the matrix, respectively. Air blast testing, the probe method, and the lead cylinder compression test were used to determine the shock wave parameters, detonation velocity, and brisance of the BN-containing on-site mixed bulk emulsion explosives. Combined with theoretical calculations, the influence of BN content on the microstructure and explosive properties was systematically studied. The test results indicate that the addition of BN did not significantly affect the stability of the internal phase droplets. As the BN content increased from 0% to 1.6%, the detonation velocity, brisance, and peak overpressure all exhibited a trend of initial increase followed by decrease: the detonation velocity increased from 3850.45m·s⁻¹ to 4724.89m·s⁻¹ and then decreased to 3903.20m·s⁻¹, with a maximum increase of 22.71%; the brisance increased from 13.86 mm to 19.87mm and then decreased to 17.18mm, with a maximum increase of 43.36%; the peak overpressure increased from 136.44kPa to 318.33kPa and then decreased to 285.41kPa, with a maximum increase of 133.31%; the specific impulse increased from 9.23Pa·s to 33.98Pa·s and then decreased to 31.99Pa·s, with a maximum increase of 268.15%. Based on the experimental results, introducing an appropriate amount of BN can significantly enhance the explosive performance of on-site mixed bulk emulsion explosives.
Impact-Induced Fracture Process and Energy Dissipation Characteristics of Copper-Bearing Albite Rock Based on FDEM
ZHANG Xiyuan, LI Xianglong, ZUO Ting, LIU Jinbao, WANG Jianguo, HU Tao, WANG Hao
, Available online  , doi: 10.11858/gywlxb.20251198
[Abstract](70) PDF (13)
Abstract:
In order to ensure the efficient recovery of copper resources, the copper-bearing albite rock samples were taken as the research object, and the impact loads with different strengths were applied by using the split Hopkinson pressure bar ( SHPB ). The crack propagation process was recorded by a high-speed camera system, and the energy dissipation law of the samples under different impact pressures was analyzed by combining the one-dimensional stress wave propagation theory and the law of conservation of energy. At the same time, a numerical model of the impact process of copper-bearing albite is established based on the finite-discrete element ( FDEM ) coupling algorithm. The results show that the incident energy and the peak stress increase with the increase of the impact pressure, and the degree of fragmentation of the sample also increases. When the incident energy is less than 140 J, the energy dissipation rate increases with the increase of the incident energy. When the incident energy is greater than 160 J, the energy dissipation rate decreases with the increase of the incident energy, and the energy dissipation rate reaches the maximum when the impact pressure is 0.35 MPa. The new crack area and the total impact energy increase with the increase of impact load. When the impact pressure is 0.30 MPa, the strain energy ratio is the smallest, indicating that the rock breaking efficiency of 0.30 MPa impact pressure is the highest. In the process of impact, tensile failure plays a dominant role and forms the main dominant area in the horizontal direction. The numerical model based on FDEM can effectively predict different shocks.
OUYANG Dehua, LIU Yuhan, PAN Jiazheng, LI Zhe, GUO Xiaoqiang, WANG Song, LIU Xingyu
, Available online  , doi: 10.11858/gywlxb.20251191
Abstract:
To enhance the safety of the non-lethal kinetic energy ammunition used in the current 38mm riot control guns in the country, the finite element - discrete element method was employed to numerically simulate the impact process of the 38mm spherical kinetic energy projectile filled with lead sand on a human body - like target. The modeling method and parameter selection were indirectly verified through a rigid - wall experiment, and data on the deformation process, kinetic energy, velocity, displacement, and energy transfer rate during the projectile impact on the target were obtained. Based on this, comparative analysis was conducted on different projectile velocities and wall thicknesses, and safety related shooting suggestions were proposed. The results show that the projectile undergoes significant deformation upon impact with the target, transforming into a disc like shape, while the target exhibits a circular indentation, with both deformations being partially recoverable to some extent. The wounding power of the projectile increases with velocity and decreases with wall thickness. The minimum safe shooting distances without causing abdominal skin penetration injuries for projectile wall thicknesses of 5mm, 7mm, and 9mm are 122.40m, 64.62m, and 31.26m, respectively.
Influence of Rock Mass Joints on Slot Blasting and Parameter Optimization Based on Discrete Element Method
SONG Yongkang, LIU Haoshan, ZHANG Zhiyu
, Available online  , doi: 10.11858/gywlxb.20251186
Abstract:
The grooving and blasting effect of the drilling and blasting method in roadway tunneling directly affects the blasting cycle efficiency, while the existing studies mostly ignore the influence of mesoscopic defects such as internal joints of rock mass. Based on the PFC 2D discrete element method, a discrete fracture network (DFN) is introduced to construct a rock mass model with different densities of joints, and the particle expansion method is used to simulate the groove blasting process, and the effects of joint density on crack propagation, energy dissipation and post-explosion block size are systematically analyzed. On this basis, the blast hole layout scheme is optimized, the original 6-hole layout is simplified to a 4-hole diamond-shaped layout, and the 15 ms differential detonation is used to improve the explosive energy utilization rate, and the post-detonation effect is similar to the original scheme. Field tests show that the optimization scheme effectively saves the actual production cost and reduces the drilling workload. The research results emphasize the importance of considering the joint defects of rock mass for the optimization of blasting parameters, and provide a theoretical basis and practical reference for efficient tunneling of rock roadways.
Effect of wire material on the energy deposition in electro-chemical coupling explosions
WANG Cheng, WANG Hangyu, LI Xinghan, WEI Ding, LIN Jiarui, CHEN Haodong, GAN Yundan
, Available online  , doi: 10.11858/gywlxb.20251173
Abstract:
To enhance the total energy output and power of energetic materials, this study employs plasma generated by electrically exploded metal wires to initiate the detonation of energetic materials, achieving coupled release of electrical and chemical energy. Using a self-built experimental system for electro-chemical coupling explosion, voltage and current curves during the explosion process were measured under ambient temperature and pressure in air. The electro-chemical coupling explosion was divided into four typical stages. The research indicates that the primary energy deposition of different metal materials occurs at distinct stages: nickel and copper wires, with their medium boiling points and high temperature coefficients of resistance, achieve efficient phase change energy deposition during the wire phase transition and current pause stages. During the plasma discharge stage, aluminum undergoes explosive vaporization due to oxide layer fracture and forms highly conductive plasma owing to its low ionization energy, leading to a significant leap in energy deposition. Tungsten, through latent heat accumulation in the liquid phase and a sharp increase in resistance, accounts for over 80% of its energy deposition during the plasma discharge stage. The study also reveals that the unique current pause phenomenon in electro-chemical coupling explosions is influenced by metal properties (such as temperature coefficient of resistance, boiling point, and latent heat of vaporization). Copper wires exhibit the longest current pause duration, while tungsten wires show no such phenomenon. This paper systematically investigates the power and energy deposition characteristics during electro-chemical coupling explosions, elucidates the influence mechanisms of metal materials on the energy release process, and provides experimental evidence and technical support for enhancing the total energy output and power of energetic materials.
A Machine Learning Potential Model for Simulating Dynamic Mechanical Response of Pb-Sn Alloy
WANG Xiaoyang, HOU Enze, WANG Han
, Available online  , doi: 10.11858/gywlxb.20251151
[Abstract](102) PDF (6)
Abstract:
Lead is a low-melting-point metal with a complex temperature-pressure phase diagram. Alloying with tin further reduces its melting temperature, making lead-tin alloys an important model material for studying dynamic mechanical responses and failure behavior. However, experimental characterization of atomic-scale dynamic failure mechanisms in Pb-Sn alloys remains challenging due to current technical limitations. Non-equilibrium molecular dynamics (NEMD) simulations can track atom trajectories and reveal key dynamic processes under dynamic loading-unloading. It thus serves as a critical alternative tool. Yet, the reliability of molecular dynamics relies on the accuracy of interatomic potentials, and currently, no high-accuracy potential exists for Pb-Sn alloys under dynamic conditions./t/nIn this work, we develop a machine-learning interatomic potential (DP-PbSn) for Pb-Sn alloys using a concurrent learning scheme. This potential achieves first-principles accuracy across a wide thermodynamic range (0-100 GPa, 0-5000 K), reliably predicting fundamental properties (e.g., lattice constants, elastic constants), defect energetics (e.g., surface energy, stacking fault energy, vacancy formation energy), as well as melting curves and shock Hugoniot curves, demonstrating its suitability for dynamic simulations. Leveraging this potential, we conduct preliminary NEMD simulations to investigate the dynamic mechanical responses of pure Pb and Pb-Sn alloys, elucidating the influence of Sn on phase transitions and plastic deformation under dynamic loading. The DP-PbSn serves as a robust theoretical tool for high-accuracy non-equilibrium molecular dynamics, providing essential insights for experimental studies on the dynamic damage behavior of Pb-Sn alloys.
A Review of Machine Learning Potentials in the Study of Materials Properties
LI Jinlong, WANG Hao, GENG Huayun
, Available online  , doi: 10.11858/gywlxb.20251172
[Abstract](102) PDF (16)
Abstract:
With the rapid advancement of artificial intelligence technologies and hardware capabilities, AI has gradually become a revolutionary tool driving transformative changes across multiple scientific research domains. In the field of materials science, machine learning methods have played a significant role in high-throughput materials design and property prediction. Over the past decade, machine learning-based approaches for constructing interatomic potentials have been widely applied in the study of material properties, providing crucial support for the theoretical design of novel materials and the in-depth understanding of their underlying microscopic mechanisms. This article reviews the development of machine learning potentials, introduces their fundamental workflows, outlines the principles of mainstream methods and their applications in materials property research, briefly discusses recent progress in emerging universal potential models, and concludes with an analysis of current challenges and future research directions.
NIU Xingang
, Available online  , doi: 10.11858/gywlxb.20251199
Abstract:
In order to explore the evolution law of damage to gas containing coal under impact, a gas containing coal separation Hopkinson pressure bar (SHPB) test system was used to conduct dynamic compression tests on coal with gas pressures of 0, 0.5, 1, 1.5, and 2 MPa. Based on energy theory, the deformation and failure process of gas containing coal under impact was analyzed, and the influence of gas pressure on coal energy parameters was explored. With the help of SMP strength criterion and Weibull distribution function, a dynamic damage constitutive model of gas containing coal considering gas impact coupling was established by combining energy consumption index. Research has shown that during the impact compression process, the time history curve of gas containing energy can be divided into a slow growth stage, an accelerated growth stage, and a stable stage; With the increase of gas pressure, the reflected energy of coal shows a linear increase trend, while the transmitted energy and dissipated energy show a linear decrease trend; The constitutive model curve of gas shock coupling damage has strong consistency with the experimental curve, which can accurately describe the damage evolution law of the entire stress-strain process of gas containing coal under impact.
Structural Optimization and Energy Absorption Characteristics of Double-Layer Variable-Diameter Energy-Absorbing Components for Anti-Impact Brackets
HAN Ruifu
, Available online  , doi: 10.11858/gywlxb.20251164
Abstract:
In order to effectively alleviate the destructive effect of ground impact on hydraulic supports, a double-layer variable energy absorbing member with higher energy absorption is proposed based on the research basis of single-layer variable energy absorbing components. Based on the energy method to analyze the energy dissipation theory of the expansion and reduction deformation of pipe fittings of different sections, the formula for calculating the bearing capacity of the components in the stable diameter reduction process under different combinations of corrugated tubes and round tubes is derived. Through numerical simulation analysis, the energy absorption curves, bearing capacity curves and deformation laws of eight different types of energy-absorbing components are obtained, and it is found that the double-layer reducer energy-absorbing member structure (SBY type) of the inner bellows and the outer round tube has better energy absorption performance. The influence of different structural parameters on the energy absorption effect was explored, and four structural parameters, namely inner tube wall thickness, outer tube wall thickness, corrugation radius and the number of inversion angles in the base, had significant effects on the energy absorption characteristic parameters. The experimental scheme was designed according to the Latin hypercube sampling method, and the structural parameters were optimized by using the Kriging surrogate model combined with the multi-objective particle swarm optimization algorithm, and finally the optimized structural parameters were selected as follows: inner tube wall thickness 6mm, outer tube wall thickness 2.9mm, corrugation radius 6.9mm, and base inversion angle 40°. The energy-absorbing components were made for axial quasi-static pressurization experiments to verify the accuracy and effectiveness of the simulation analysis and optimization results. The results show that the total energy absorption of the double-layer variable diameter energy-absorbing components is increased by 54.2%, the specific energy absorption energy is increased by 55.6%, the average bearing capacity is increased by 43.2%, and the standard deviation of the load is increased by 59.5%. The designed components have better energy absorption performance, which can make the yield anti-erosion process more reliable, and provide a theoretical basis and reference for the design of energy-absorbing components of deep roadway support hydraulic supports.
Metallic Hydrogen Ligand Compounds: A Potential Route to Superconducting Metallic Hydrogen at Ambient Pressure
ZHANG Zihan, DUAN Defang, CUI Tian
, Available online  , doi: 10.11858/gywlxb.20251216
[Abstract](123) PDF (16)
Abstract:
Metallic hydrogen, with its properties including room-temperature superconductivity and quantum fluidity, is known as the holy grail of high-pressure physics. However, since atomic metallic hydrogen requires pressures about 500 GPa, it has not been realized in experiments since its conception in 1935. To take advantage of properties the properties of metallic hydrogen in the future, it will be crucial to obtain it at ambient pressure. Current approaches to obtaining metallic hydrogen at low pressures rely on the "chemical precompression" in hydrides to induce metallization of hydrogen at low pressures, essentially identifying superconducting hydrides that can host the properties of metallic hydrogen. However, these superconducting hydrides currently lack distinct structural features, complicating the search for metallic hydrogen hosts. Here, we identify metallic hydrogen ligand compounds with hydrogen as the ligands as potential hosts for properties of metallic hydrogen at low pressures. The metallization of the non-bonding orbitals of the hydrogen ligands is a key criterion for determining whether a metallic hydrogen ligand compound can host metallic hydrogen properties. This article summarizes the main behaviors of hydrogen at ambient pressure, focusing on hydrogen ligand compounds at ambient pressure. Then, using a simple model of a one-dimensional hydrogen atom chain, we analyzed the causes of non-bonding orbital metallization and the physical picture of reduced stability pressure. The orbital characteristics of metallic hydrogen ligand compounds are then analyzed, highlighting their rules of superconductivity, topological properties, and the electronic structure that enable metallization. The analysis of metallic hydrogen ligand compounds presented in this article not only provides important structural information for future exploration of metal hydride superconductors but also provides an important theoretical foundation for realizing the properties of metallic hydrogen at ambient pressure.
Inhibition Mechanism of KHCO3-Containing Water Mist on Methane/Hydrogen Premixed Deflagration
HUANG Hui, LI Yuanbing, LI Xia, SHAO Peng
, Available online  , doi: 10.11858/gywlxb.20251189
Abstract:
Explosion prevention and mitigation technologies for hydrogen/methane gas mixtures represent a critical research area for ensuring the safe application of hydrogen energy. This study systematically investigates the inhibition mechanism of potassium bicarbonate (KHCO3)-containing fine water mist on methane/hydrogen premixed deflagration using a combined approach of experiment and numerical simulation. The results indicate that KHCO3-containing fine water mist exhibits a significant inhibitory effect on methane/hydrogen premixed deflagration, with its suppression performance positively correlated to the KHCO3 mass fraction. Taking the condition of X_(H_2 )=10% as an example, 11 wt% KHCO3 addition resulted in reductions of the maximum explosion pressure P_max and the average rate of pressure rise 〖(dp/dt)〗_avg by 34.64% and 44.57%, respectively. The laminar burning velocity was reduced by up to 66.43%. KHCO3 contributes to suppression through both physical and chemical mechanisms. Physically, droplet phase change (evaporation) absorbs heat and the generated steam dilutes the fuel mixture, thereby lowering the flame temperature and reducing reactant concentrations. Chemically, the decomposition of KHCO3 generates potassium compounds, which undergo the KOH → K → KOH recombination cycle to scavenge key radicals (•H, •O, •OH). This process competes with chain-branching reactions and interrupts the combustion chain reactions.Furthermore, the suppression process is governed by a competition between inhibitory and promotional effects. At high hydrogen blending ratios and high mass fractions, the physical evaporation efficiency becomes a bottleneck that constrains the chemical inhibition, leading to a saturation of the overall suppression efficiency. Nevertheless, a significant inhibitory effect is still maintained.
Preliminary XFEL Experimental Simulation Platform Based On HSWAP Engine
LIU Jin
, Available online  , doi: 10.11858/gywlxb.20251155
Abstract:
X-ray Free Electron Laser (XFEL) plays a critical role in diagnosing dynamic compression processes in micro- and meso-scale materials. To deepen our understanding of XFEL physics and optimize facility design, a preliminary XFEL experimental simulation platform was developed based on the High-Performance Computing (HPC) Simulation Workflow Application Platform (HSWAP). HSWAP provides workflow, component, and data linkage models for XFEL experiments, enabling flexible simulation of diverse processes through modular configurations. This platform was employed to investigate X-ray diffraction (XRD) of microscale materials and phase contrast imaging (PCI) of meso-scale explosive samples. Simulation results for XRD of a metallic sample under shock loading and PCI of voids in explosive materials demonstrate the platform's ability to accurately reproduce experimental dynamics. By integrating numerical models with data analysis, the platform enhances the design of XFEL experiments and provides a foundation for interpreting diagnostic capabilities in ultrafast processes. Future work will focus on refining simulation methods for meso-scale samples using phase-field approaches and high-Z materials under shock conditions.
Advances in the Application of Machine Learning to the Calculation of Mineral States and Properties in the Earth's Deep Interior
WANG Chuan, CENG Qiyu, CHEN Bo, YU Xiaoxiang, KANG Dongdong, DAI Jiayu
, Available online  , doi: 10.11858/gywlxb.20251218
[Abstract](119) PDF (14)
Abstract:
The deep interior of the Earth is under extreme high-temperature and high-pressure conditions. Its material composition, phase transition behavior, and physical properties are crucial for understanding the Earth's internal structure, dynamic processes, and evolution. Traditional experimental methods face challenges in maintaining thermodynamic states and diagnosing physical quantities under such extreme conditions. While first-principles calculations offer quantum-level precision, their computational efficiency limits their direct application to simulating deep-Earth minerals across large spatiotemporal scales. Machine learning methods present new opportunities. By constructing high-precision, efficient machine learning potential functions based on first-principles datasets, machine learning methods significantly extend the spatiotemporal scale of first-principles simulations, which provide revolutionary tools for studying the physical states, phase transitions, elasticity, and transport properties of deep-Earth minerals. This paper systematically reviews the progress of machine learning applications in studying major deep-Earth minerals, including those in the upper mantle, transition zone, lower mantle, subduction zone components, and core materials, and summarizes the representative achievements of machine learning methods in revealing phase transitions, thermal conductivity, diffusion, melting, and elastic properties, while also discussing current limitations and future research directions.
Bursting Performance Optimization of Reverse-Arched Bursting Discs Based on Variable Fidelity Surrogate Models
YU Yaowen, LIANG Hao, CHEN Zhanghai, PU Weiqiang
, Available online  , doi: 10.11858/gywlxb.20251123
Abstract:
To address the optimization design problem of the bursting performance of reverse-arched bursting discs (RABDs), a hierarchical Kriging (H-Kriging) surrogate model was constructed based on both high- and low-fidelity finite element analysis results. This model enables the rapid prediction of the burst pressure of RABDs, facilitating the development of a mathematical model for performance optimization and structural improvement. The results show that the H-Kriging surrogate model relating burst pressure to structural parameters based on high- and low-fidelity finite element models can significantly reduce computational cost while accurately predicting the burst pressure of RABDs. For the initial structural design scheme of RABDs, optimization was carried out using a genetic algorithm, with the optimized design accounting for manufacturing tolerance in disc thickness. This resulted in a 58.8% reduction in burst pressure fluctuation, significantly reducing the sensitivity of burst pressure to thickness manufacturing errors and providing valuable engineering reference.
Numerical Investigation on Cavity Evolution and Motion Characteristics of High Speed Water Entry Ogival Projectiles with Different Headforms
ZHENG Xiaobo, SONG Haisheng, ZOU Daoxun, YAO Weiguang, LI Teng, GUI Yulin, HE Yu, CHEN Yonglong
, Available online  , doi: 10.11858/gywlxb.20251169
Abstract:
At present, The trans-media weapon is one of research hotspot in the military field. Based on the reynolds time-averaged N-S equation, VOF multiphase flow model and modified Realizable k-ε turbulence model, a three-dimensional numerical simulation method is constructed to study the cavity evolution and motion characteristics of ogival nose projectiles with different head shapes during high-speed vertical water entry, and the influence of head shapes on the cavitation evolution and motion characteristics is analyzed. The results show that the numerical simulation and experimental data have good consistency in the evolution of cavity shape and projectiles velocity. The geometry of the projectile warhead significantly affects the formation mechanism of the cavity and the motion characteristics of the projectiles. The cavity of ogival nose projectiles and double-cone ogival nose projectiles initially appears in the shoulder area of the projectile body, while the cavity of cone-cylinder ogival nose projectiles starts in the head and quickly wraps the entire projectile body. Combined with the analysis of the fluid pressure field, it is shown that a low-pressure area appears on the double-cone ogival nose projectile, which leads to the slowdown of the projectile velocity attenuation. The head of the cone-cylinder ogival nose projectile forms a typical high-pressure area, which leads to the acceleration of the projectile velocity attenuation. In addition, the axial acceleration of the cone-cylinder ogival nose projectile is more than twice that of the other two projectiles.
Molecular Dynamics Simulation of Micro-Jetting and Spallation in Helium-Bubble Copper under Double Supported Shocks
WANG Xinxin, BAO Qiang, HE Anmin, SHAO Jianli, WANG Pei
, Available online  , doi: 10.11858/gywlxb.20251075
Abstract:

Micro-jetting and micro-spallation at metal interfaces under intense shock loading play pivotal roles in applications such as inertial confinement fusion (ICF). These phenomena exhibit inherent complexity due to their multi-scale dynamics, strong nonlinearity, and coupled multi-field interactions. Under extreme irradiation conditions, the formation of high-pressure nanoscale helium bubbles significantly alters interface failure mechanisms. Using molecular dynamics methods, we investigate micro-jet growth and damage evolution in helium-containing copper subjected to double supported shock loadings. Helium bubbles demonstrate lower critical activation stress thresholds for expansion compared to void nucleation, with these thresholds being dependent on bubble distribution and number density. Under low-pressure primary shocks, helium-containing metals produce more pronounced micro-jets than pure metals. During secondary shocks, helium bubbles promote jet fragmentation, resulting in higher maximum velocities at micro-jet tips while maintaining comparable velocity distributions in micro-jet bodies. Secondary shocks show negligible effects on bulk helium bubbles that were previously compressed by initial shocks and partially rebounded due to rarefaction waves without complete recovery. Near-surface ruptured bubble walls may reattach to bubble bases after secondary shocks, temporarily re-trapping helium atoms that are subsequently released during unloading-induced re-expansion and rupture. The collapse mechanism of helium bubbles under secondary shock is closely related to the helium bubbles size and the strength of secondary shock. This study establishes fundamental physical understanding and provides a theoretical foundation for future cross-scale investigations of coupled micro-jetting and micro-spallation evolution in irradiated helium-containing metals.

Crystal Structure and Physical Properties of Sr2He Compound under High Pressure
WANG Qingmu, ZHANG Pan, SHI Jingming, LI Yinwei
, Available online  , doi: 10.11858/gywlxb.20251084
Abstract:

By combining first-principles calculations under the framework of density functional theory (DFT) and the CALYPSO crystal structure prediction method, the structural stability of the inert element helium (He) and alkaline-earth metals under high-pressure conditions has been systematically investigated. The calculations reveal that among the alkaline-earth metals, strontium (Sr) forms compounds with He exhibiting relatively low energy values. Consequently, the crystal structure of Sr2He at 400 GPa was predicted. Electron localization function (ELF) and density of states (DOS) analyses show no tendency for covalent bond formation between Sr and He atoms. Furthermore, Bader charge analysis reveals ionic bonding between Sr and He atoms, with charge transfer occurring from He to Sr. These results provide key insights into the bonding mechanism of Sr2He. This study elucidates the crystal structure, bonding nature, and electronic properties of Sr2He, offering theoretical support for understanding the stability and physical properties of such metastable materials and providing important guidance for their experimental synthesis.

Experimental Study on Directional Rock Fracture by Energy-gathered Cutting under Dynamic Impact
WANG Wei, CHENG Mingfeng, LUO Xin, WANG Jinbao, ZOU Baoping, CAO Chunhui
, Available online  , doi: 10.11858/gywlxb.20251114
[Abstract](102) PDF (3)
Abstract:
This study addresses the challenge of excavating through heterogeneous tuffaceous sandstone formations in tunnel construction by proposing a novel energy-gathering slotting rock-breaking technique based on dynamic impact. Using self-developed geotechnical dynamic impact testing system, cylindrical tuffaceous sandstone specimens (Φ100 mm × 50 mm) were prepared with 10 mm thick polyurethane pads adhered to one end. Radially arranged holes of 3 mm, 6 mm, and 9 mm in diameter were drilled into the pads, each fitted with six corresponding energy-gathering nails. Seven groups of tests were conducted under impact air pressures ranging from 0.35 to 0.65 MPa to investigate the effects of varying impact energy and nail diameter on directional rock fracturing performance. The results show that as the impact pressure increases, the peak stress and energy absorption of the specimens rise significantly, with fracture patterns transitioning from primarily intergranular to transgranular cracking. The 3 mm nails were prone to local crushing and failed to produce effective through-cutting cracks, while the 9 mm nails caused blocky or pulverized failure under high pressure. In contrast, the 6 mm nails consistently induced stable, continuous, and directional fractures under various pressures, producing more transgranular cracks and demonstrating excellent energy utilization efficiency. Scanning electron microscopy confirmed the strain-rate effect of impact: cracks were predominantly intergranular under low strain rates (low impact forces) and became transgranular under high strain rates. This technique leverages the compressive-reflective-tensile stress chain mechanism inherent in dynamic fracture mechanics to achieve controlled, directional rock breaking without explosives or liquid media. By properly matching impact parameters and nail diameters, this method can efficiently guide crack propagation along predetermined paths in deep, heterogeneous rock masses, offering a promising strategy for controlling over- and under-excavation in complex geological tunneling conditions.
Discrete Element Simulation of Axially Compressed Energy Constitutive Relations in Defective Sandstone
chenbing
, Available online  , doi: 10.11858/gywlxb.20251142
[Abstract](107) PDF (4)
Abstract:
Dual Wavelength All-Fiber Laser Interferometric Velocity Retest Technique
WANG Wei, LIU Shenggang, GU Wei, TAO Tianjiong, MA Heli, WANG Xiang, WENG Jidong
, Available online  , doi: 10.11858/gywlxb.20251081
Abstract:
Based on the wavelength division multiplexing technology and the traditional optical path structure design of laser velocity interferometer, a dual wavelength laser velocity interferometer which has the retest function is proposed. The basic principle and optical path structure design of this dual wavelength laser velocity interferometer are briefly introduced. The dynamic examination experiment of low speed and high speed are carried out on gas guns, and the velocities of the same measuring point are simultaneous measured by this method with one single probe.The experimental results show that the deviation of the two wavelengths velocity measurement is within 1.5%.
Laminar Combustion and Explosion Characteristics of Ternary Premixed Fuels under High Pressure
CHEN Rui, JIANG Genzhu, TAO Juxiang
, Available online  , doi: 10.11858/gywlxb.20251140
Abstract:
The effect of ethanol addition on the combustion characteristics of hydrogen methane mixed fuel was studied. Based on the constant volume combustion system, The laminar combustion and explosion characteristics of the three-component fuel were studied at the initial temperature of 400K, different ethanol contents (20%, 50% and 80%), different pressures (1bar, 3bar and 4bar) and different equivalent ratios (φ=0.7-1.4). The results show that under all experimental conditions, the hydrodynamic instability is most obvious when the equivalent ratio is 1.1, which is positively correlated with pressure and ethanol concentration. The laminar combustion velocity of the pre-mixed fuel is linearly affected by the initial pressure and the concentration of ethanol. The fitting correlation is in good agreement with the experimental values, and the deviation is less than 7%. The kinetic analysis shows that R1 pair is the main element reaction to increase the flame velocity of laminar flow. The linear relationship between explosion characteristics and initial pressure was obtained. The relationship between linear correlation slope, intercept and equivalent ratio was accurately determined by polynomial fitting, and the quantitative relationship between explosion characteristics, equivalent ratio, initial pressure and ethanol concentration was obtained.
Mechanical Properties and Ignition Performance of Rare Earth Reactive Materials under Impact Loading
LI Shoujia, ZHANG Beichen, DOU Jihang, HAN Yuhang, ZHAO Hongwei, CHEN Xuefang, QIN Shuaiwei, LU Xiaoxia, BI Pengyu
, Available online  , doi: 10.11858/gywlxb.20251106
[Abstract](136) PDF (6)
Abstract:
Aluminum (Al), one of commonly used reactive metals, is widely appliedin reactive material systems. However, its relatively low reactivity restricts the energy release of systems. To improve the reactivity of Al, we introduced aluminum-cerium alloy (Al-Ce alloy) which include the rare-earth cerium with high reactivity into the system. This study investigated the mechanical properties and ignition performance of four reactive material systems under shock overload—Al2Ce/PTFE, Al/PTFE, Al2Ce/ammonium perchlorate (AP), and Al/AP. A split Hopkinson pressure bar (SHPB) system was used to study the dynamic stress-strain behavior, ignition delay, and combustion duration of the prepared samples. Thermal analysis was conducted to assess the influence of reactive metal content on the thermal decomposition of AP. Results showed there are three distinct shock-induced ignition modes: non-ignition, combustion, and combustion (deflagration). Both Al2Ce/PTFE and Al/PTFE exhibited poor ignition performance. The Al2Ce/AP system demonstrates higher ultimate strength and critical failure strain, achieving deflagration upon impact with significantly shorter ignition delay and combustion duration compared to Al/AP. The incorporation of cerium accelerates AP decomposition and substantially increased the enthalpy of the Al2Ce/AP system, resulting in more concentrated energy release. Ce effectively enhances the reactivity of aluminum, and its high reactivity accelerates the reaction kinetics of the reactive system. Besides, it significantly intensifies energy release under impact loading. In conclusion, rare earth aluminum alloy materials, due to their high reactivity advantage, are of great significance for the development of new aluminum-based impact reaction materials.
Numerical Simulation of Penetration Response in Salt Cavern Gas Storage Based on RHT Constitutive Model
HU Rongrong, FAN Jinyang, YANG Fan, LIANG Wuxing, JIANG Deyi
, Available online  , doi: 10.11858/gywlxb.20251119
[Abstract](174) PDF (10)
Abstract:
Underground salt cavern gas storage is an important energy infrastructure, and once it is damaged by impact, it will cause irreparable losses. Therefore, it is of great significance to propose key dynamic stability indexes to evaluate the safety of salt caverns under extreme impact loads. To investigate the dynamic response of salt cavern gas storage under high-speed penetration, the salt rock material was defined using the RHT constitutive model. A finite element model of the gas storage structure was established using ANSYS/LS-DYNA software to analyze the damage effects of a ground-penetrating weapon on the salt cavern structure. Numerical simulations were conducted for three scenarios with varying overburden thicknesses, focusing on four parameters: vertical displacement, vertical stress, effective plastic strain, and shear stress. These simulations revealed the failure mechanisms of the cavern roof and surrounding rock under dynamic impact, as well as the variation patterns of key stability indicators. The results demonstrated that reduced overburden thickness intensified the dynamic response of the surrounding rock and expands plastic deformation zones. Displacements of the roof and surrounding rock exhibited an initial increase followed by a decrease. Salt rock in low vertical stress regions experienced higher shear stresses, increasing susceptibility to failure. The surrounding rock accumulated greater plastic strain, showing heightened sensitivity to penetration-induced disturbances.
Characterization of Damage to Adjacent Fill Bodies by Blasting of Slit Packets
ZHU Benliu, LI Xianglong
, Available online  , doi: 10.11858/gywlxb.20251111
[Abstract](138) PDF (2)
Abstract:
In order to accurately regulate the damage effect of slit pack blasting on the filling body of the quarry in deep mines, this study focuses on the damage control mechanism of the peripheral hole spacing (500mm, 600mm, 700mm, 800mm). Based on the theory of elastic fluctuation and the dynamic propagation characteristics of shock waves in rocky media, the diffusion mechanism of the stress wave under the action of multi-media in the constrained orientation during slit packet blasting is established; and combined with the strong correlation between brittle concrete materials and the damage evolution of the filling body, the cross-media equivalence calibration framework of the RHT intrinsic model is established; and the “filling body-mining body” model is constructed on the basis of the numerical simulation software, ANSYS/LS-DYNA, and the “filling body-mining body” model is developed by using the numerical simulation software, ANSYS/LS-DYNA. Based on the numerical simulation software ANSYS/LS-DYNA, we constructed a multi-media dynamic coupling numerical model of “filling body-mineral body-cutting slit package”, arranged observation points at the junction of filling body-mineral body, and conducted a combined analysis of the peak stress change, the change of the blast vibration velocity, and the damage evolution of the filling body at the observation points. Then, based on the blasting test of the approach and return stage of the neighboring filling body in Jinchuan Three Mining Area, the blasting test of conventional packs, slit packs and different peripheral hole spacing was conducted. The test shows that: slit pack blasting triggers gas-phase jet and strain-energy convergence effects in the unconfined direction, synchronously suppresses the stress and vibration peaks in the confined direction, and achieves directional attenuation of the blasting load on the neighboring filling body; the field test shows that, compared with the conventional charge, the slit pack significantly reduces the degree of damage of the filling body by more than 36%; the degree of blasting damage and the peripheral hole spacing show a negative correlation, and the damage suppression efficiency is improved with the increase of the spacing. The damage suppression efficiency is improved when the spacing increases.
LI Tao, NI Yu, WANG Zhiliang
, Available online  , doi: 10.11858/gywlxb.20251116
[Abstract](114) PDF (3)
Abstract:
Keyword:Aiming at the problem that traditional peripheral blasting easily induces random cracks and cause surrounding rock failure, this study investigates the damage evolution mechanisms and dynamic response characteristics of presplitting hole directional blasting by integrating elastic mechanics theory with ANSYS/LS-DYNA numerical simulations. Initially, based on elastic mechanics theory, the mechanical mechanism by which presplitting holes guide crack propagation through reflected stress waves generating tensile stress concentration under explosive loading is elucidated. Subsequently, a numerical model of planar double-hole decoupled charge was established to systematically explore the influence of borehole spacing and in-situ stress field on the damage evolution. Finally, the dynamic variation patterns of peak pressure and peak particle vibration velocity near presplitting holes were analyzed. The results demonstrate that: (1) empty hole generates tensile stress concentration through reflected stress wave, significantly improving directional crack penetration while suppressing undesired crack bifurcation. (2) As the borehole spacing increases, the stress concentration effect of empty holes weakens, the peak pressure and particle peak velocity near the hole walls decrease, and the through-going cracks between holes also decrease. (3)Under the action of in-situ stress, emptyhole maintains its directional control capability within the stress field. With increasing in-situ stress, the horizontal peak pressure near empty holes increases, resulting in inhibited crack propagation between boreholes and weakened crack-guiding effect of empty holes. It is necessary to adjust the hole spacing to promote the formation of penetrating cracks.
High-Pressure Study on Structural Phase Transformation and Physical Properties of SrB2C2
GUO Hua, WANG Fan, ZHENG Baobing
, Available online  , doi: 10.11858/gywlxb.20251148
[Abstract](149) PDF (15)
Abstract:
Using the first-principles calculations based on density functional theory and the structure prediction method based on particle swarm optimization algorithm, this paper conducted a structural search for SrB2C2 in the pressure range of 0-350 GPa, and successfully determined that SrB2C2 belongs to the tetragonal system tI20-SrB2C2 under normal pressure, and belongs to the orthorhombic system oF40- SrB2C2 under high pressure. Based on the enthalpy difference curve, the phase transition pressure of SrB2C2 was determined to be 44.7 GPa. The stability of tI20-SrB2C2 and oF40- SrB2C2 at the corresponding pressure and the possibility of experimental synthesis were verified by calculating the phonon spectrum, elastic constants and formation enthalpy. It can be seen from the Young's modulus and shear modulus as a function of orientations that tI20-SrB2C2 has higher degree of mechanical anisotropy than oF40- SrB2C2, which can be ascribed to the fact that the sp2-hybridized boron-carbon bonds form the layered structure of tI20-SrB2C2, while the boron-carbon bonds of oF40- SrB2C2 are mainly sp3-hybridized covalent bonds, forming a more stable three-dimensional tetrahedral network structure. The calculation of the electronic structure shows that SrB2C2 is an indirect band gap semiconductor, and the calculation of the electronic localization function shows that the boron-carbon bonds in tI20-SrB2C2 and oF40-SrB2C2 are sp2 and sp3 covalent bonds, respectively.
Energy Conversion Characteristics and Temperature Field Simulation of Ignition Resistors for Industrial Electronic Detonators
LI Hongwei, ZHANG Liguo, ZHOU En, LIANG Hao, YANG Lin, ZHANG Wanlong
, Available online  , doi: 10.11858/gywlxb.20251124
[Abstract](149) PDF (3)
Abstract:
In order to study the dynamic changes of temperature during the ignition process of the ignition resistor and solve the matching problem of the energy storage capacitor-ignition resistor system, this paper measures the voltage-current changes and temperature changes of the ignition resistor under different capacitor discharge voltages through electrothermal experiments, infrared temperature measurement and numerical simulation methods. Combined with the surface conditions of the unmelted samples, Determine the critical fuse voltage of the ignition resistor, and obtain the law of the electrical characteristic curve and the temperature variation law of the ignition resistor. The results show that under the same voltage, the bridge membrane linear type has the shortest melting time. Under the same resistance value, the melting time and heating time of the bridge film ignition resistor are shorter than those of the bridge wire type, and the maximum temperature it can reach is also higher. During the power-on process of the bridge wire type and bridge film type S-shaped ignition resistors, heat is prone to accumulate at the corners and phase change occurs first.
The Disappearing Quartz-Coesite Path: the Phase Transition Mechanism of Silicon Dioxide from Machine Learning Simulations
ZHU Shengcai, DENG Pu, HOU Rui, ZHAO Yingliang
, Available online  , doi: 10.11858/gywlxb.20251122
[Abstract](184) PDF (15)
Abstract:
This study systematically investigates the structural phase transition mechanisms of silicon dioxide under high pressure using a high-dimensional neural network potential model combined with the stochastic surface walking algorithm. First, a global potential energy surface encompassing quartz, coesite, stishovite, and amorphous states was constructed, and the thermodynamic phase diagram was plotted, revealing the thermodynamic stability advantage of stishovite in high-pressure regions. Further analysis demonstrated that the energy barrier for the quartz-to-stishovite transition path significantly decreases under high pressure, indicating strong kinetic feasibility, whereas the coesite-to-stishovite pathway follows a single transition state mechanism with a slightly increasing energy barrier under pressure. Regarding the amorphization transition, the key role of the symmetry-deficient low-energy structure group in the high-pressure amorphization of quartz was clarified through sampling and identification, unveiling the "short-range order—medium-range disorder" structure as a defining characteristic of the amorphous state. Notably, no effective quartz-coesite transition path was observed during the study, suggesting that the amorphization transition inhibits this transformation pathway and revealing the kinetic principles underlying the absence of the quartz-coesite transition. This work comprehensively explores the mechanisms of crystalline and amorphous phase transitions in silicon dioxide under high pressure and provides theoretical foundations and methodological paradigms for high-pressure simulation studies of complex oxides.
Microstructure and Properties of the Energetic Structural Material of Ti1.5ZrNbMo0.5W0.5 High-Entropy Alloy
WU Xiaohan, HE Jinyan, ZHUANG Zhihua, ZHANG Xinggao, PENG Wenlian, CHEN Hao, XU Hanqing
, Available online  , doi: 10.11858/gywlxb.20251105
[Abstract](165) PDF (15)
Abstract:
With the increasing demand for enhanced mechanical properties and energy release capabilities in energetic structural materials, traditional materials struggle to concurrently achieve both high mechanical properties and energy release properties. In this study, a novel Ti1.5ZrNbMo0.5W0.5 high-entropy alloy was developed by powder metallurgy process, and its microstructure, mechanical properties, damage effectiveness and energy release mechanisms were comprehensively investigated. The result indicates that the sintered Ti1.5ZrNbMo0.5W0.5 alloy, characterized by high density and fine grain size, demonstrates superior quasi-static and dynamic compression properties. During the ballistic gun experiments, the Ti1.5ZrNbMo0.5W0.5 alloy projectile can penetrate the Q235 steel plate with thickness range of 6-10 mm at the speed range of about 600-1100 m/s. Meanwhile, after penetrating through the target, the fragment was broken into small-sized fragments and causing the severe energy release reaction. This energy release reaction is primarily driven by the substantial oxidation of Zr-rich regions, releasing significant thermal energy and successfully igniting the cotton and gasoline placed behind the steel target. This research provides a thorough characterization of the microstructure and mechanical properties of Ti1.5ZrNbMo0.5W0.5 alloy. Furthermore, it evaluates its overall performance in practical armor-piercing application and reveals its energy release mechanisms. The research results provides a theoretical foundation and experimental data for the further study and application of Ti-Zr-Nb-Mo-W system high-entropy alloy.
Prediction Model and Application of Rock Burst Tendency in Deep High Stress Areas
QI Yun, BAI Chenhao, DUAN Hongfei, DAI Lianpeng, LI Xuping, WANG Wei
, Available online  , doi: 10.11858/gywlxb.20251103
[Abstract](198) PDF (12)
Abstract:
To ensure the construction safety of geotechnical engineering in deep high stress areas, a combined rock burst intensity prediction model based on Whale Optimization Algorithm (WOA) and Extreme Gradient Boosting Tree (XGBoost) is proposed to address the suddenness and complexity of rock burst. Firstly, the main controlling factors that affect the intensity level of rock burst are analyzed, and the uniaxial compressive strength, maximum tangential stress, uniaxial tensile strength, brittleness coefficient, stress coefficient, and elastic energy index are selected to establish a prediction index system for rock burst intensity level. The original samples are processed using the Pearson correlation coefficient, multiple Imputation by Chained Equations (MICE), synthetic minority oversampling technique (SMOTE), and principal component analysis (PCA). Secondly, the maximum number of iterations, maximum depth of the tree, and learning rate of the XGBoost model were optimized through WOA, and the prediction results of the model were comprehensively evaluated using accuracy, precision, recall, F1 score, and Cohen Kappa coefficient. Finally, the model was applied to predict the rock burst intensity level of the Qinlingzhongnanshan highway Tunnel and the water diversion system for hydropower stations. Results show that the WOA-optimized XGBoost model achieves optimal performance when the maximum number of iterations, maximum tree depth, and learning rate are 51, 13, and 0.7325, respectively. Prediction results for rock burst intensity level using the WOA-XGBoost model outperform those of other intelligent algorithm models, verifying the model's high accuracy and reliability in predicting rock burst intensity level.
High-Pressure Preparation of High-Strength Wood Materials
YE Zi, ZHOU Xuefeng, XU Jianing, ZHOU Chenglin, YANG Yi, ZHENG Linpeng, CHEN Bin
, Available online  , doi: 10.11858/gywlxb.20251127
[Abstract](123) PDF (13)
Abstract:
Mineral resources on Earth are finite, but wood is renewable. Consequently, replacing limited industrial materials with modified wood remains a long-term human pursuit. This study processed samples of three wood types: (Balsa Ochroma lagopus, Basswood Tilia tuan, and African blackwood Dalbergia melanoxylon). The research team used a large volume cubic press to compress these samples at room temperature under high pressure. The effects of high-pressure treatment on the air-dry density, compressive strength, and elastic modulus of the three wood species were analyzed, and changes in their internal microstructures were observed using CT and SEM. The results showed that the physical and mechanical properties of all three wood species improved. After high-pressure processing at 5.5 GPa, the density of lightwood, linden, and blackwood increased by 239%, 112%, and 11%, respectively; the surface hardness increased by 79%, 46%, and 15%, respectively; and the compressive strength increased by 33%, 9%, and 28%, respectively. Notably, the specific strength of compressed African blackwood (101.55 kJ·kg-1) approaches that of aluminum alloys(109.23 kJ·kg-1). Results demonstrate that African Blackwood is lighter than ceramic materials. Furthermore, this wood offers superior electrical insulation and thermal insulation compared to aluminum alloy. Crucially, African Blackwood possesses high specific strength. This property gives it significant potential to replace aluminum alloy in numerous special environments. Such application supports sustainable development for future industries. Ultimately, this research opens new possibilities for high-value wood applications.
Static/Dynamic Mechanical Properties and Ballistic Behavior of 6061 Aluminum Alloy
FENG Zhijian, HU Menglei, ZHANG Xuefeng
, Available online  , doi: 10.11858/gywlxb.20251108
[Abstract](250) PDF (21)
Abstract:
Due to its excellent mechanical properties, aluminum alloy has significant engineering application value in aerospace, ship engineering, new energy, and electronic devices. However, during its service, it often needs to withstand dynamic impact loads. Studying its mechanical response behavior under high strain rate conditions has important theoretical and practical significance for engineering applications. This study takes 6061 aluminum alloy as the research object and conducts in-depth research on its static and dynamic mechanical properties and ballistic response characteristics through systematic experimental tests and numerical simulations. The experimental results show that within the strain rate range of 0.001 - 3800 s-1, 6061 aluminum alloy exhibits a significant strain rate strengthening effect. The flow stress increases significantly with the increase in strain rate, by 18.5%. However, its strain hardening behavior remains relatively stable under different strain rate conditions. The parameters of the Johnson-Cook constitutive model calibrated by the least square method can accurately describe the mechanical response of the material under different strain rates. The ballistic experiment results show that the ballistic limit of a spherical projectile penetrating a 6061 aluminum alloy target plate is 282.6 m/s, and the residual velocity has a good linear relationship with the incident velocity under the super-ballistic limit condition. The failure morphology analysis of the target plate reveals that its failure mode is related to the impact velocity: at low impact velocities, the main failure mode is overall deformation dominated by composite stress, while at high penetration velocities, it is mainly local shear failure. The finite element model established successfully reproduces the ballistic response and failure mode observed in the experiments, with an error of less than 5%, verifying the reliability of the fitted constitutive model parameters and numerical methods. Using an experimentally verified finite element model, the ballistic responses of spherical projectiles of different diameters penetrating a 6061 aluminum alloy target plate were studied. When the projectile diameters were 10mm, 8mm, and 6mm, the ballistic limit velocities of the target plate were 283m/s, 392m/s, and 443m/s respectively. Therefore, under the condition that the thickness of the target plate remains unchanged, the higher the projectile mass, the greater the ballistic limit velocity of the target plate. This study provides important theoretical basis and experimental data support for the engineering application of 6061 aluminum alloy under impact load conditions.
Shock Wave Propagation Law of Curved Tunnel and Curved Diffusion Tunnel
CHENG Hao, PENG Yong, XUE Xiaoguang, LU Qiu, LI Xiangyu, LI Zhibin
, Available online  , doi: 10.11858/gywlxb.20251099
[Abstract](108) PDF (16)
Abstract:
In view of the unclear attenuation law of shock wave in curved tunnel is unclear, the influence of radius and turning angle on shock wave propagation in curved tunnel is analyzed. It was found that its influence on the wave dissipation efficiency was limited, and the wave dissipation efficiency of curved tunnel is similar to that of direct turning tunnel with the same angle, which is basically less than 7.2%. In order to improve the attenuation efficiency of shock wave in curved tunnel, a new protective idea of constructing arc-shaped diffusion tunnels based on arc-shaped tunnels by setting up diffusion chambers was proposed. The influence laws of diffusion ratio and diffusion forms (inner diffusion, two side diffusion and outer diffusion) on the wave dissipation efficiency of curved diffusion tunnels were also discussed. The calculation shows that curved diffusion tunnel can greatly improve the attenuation efficiency of shock wave, and the attenuation rate can reach 55.9%. Among them, the outer diffusion curved tunnel has the highest wave dissipation efficiency, followed by the inner diffusion type and the double-sided diffusion type. Moreover, the wave attenuation efficiency increases continuously with the increase of the diffusion ratio. As the peak pressure of the shock wave increases, the wave attenuation efficiency of the curved diffusion tunnel also improves, reaching up to 64.4%. When the peak pressure continues to increase, the wave attenuation efficiency of the curved diffusion tunnel slightly decreases but remains basically unchanged. The wave attenuation efficiency of the curved diffusion tunnel decreases with the increase of the positive pressure duration of the shock wave. When the positive pressure duration is 100 ms, the wave attenuation efficiency drops to 25.4%. However, as the positive pressure duration further increases, the wave attenuation efficiency of the curved diffusion tunnel remains almost unchanged.
Physical Mechanisms of "Shock Cooling" at the Molecular Fluid/Window Interface under Shock Loading
LI Kewei, Muhammad Sabeeh Akram, YANG Lei, YUAN Wenshuo, LIU Fusheng
, Available online  , doi: 10.11858/gywlxb.20251092
[Abstract](168) PDF (13)
Abstract:
The physical mechanism of "shock cooling" at the molecular fluid/window interface has long puzzled the shock wave physics community. There are three distinct viewpoints explaining the cooling effect at the shock interface: (1) thermal equilibrium between the molecular fluid and the window; (2) extinction effect of the molten optical window; and (3) shock response characteristics of the molecular fluid. This paper comparatively investigates the shock radiative behavior and radiation temperature variation characteristics at the interfaces between the chemically active fluid CHBr3, the inert liquid argon (LAr), and the LiF optical window. Under the same shock pressure, the interface radiation characteristics of the two media exhibit distinct evolution features, indicating that the interface cooling effect is closely related to the fluid medium and its chemical activity. Therefore, the observational results of this paper strongly support that the interface cooling effect is caused by the shock response of the fluid itself, rather than the heat conduction mechanism or the window melting extinction mechanism.
Damage of a New Shaped Warhead to Water-Bearing Composite Structure
CUI He, FU Jianping, REN Kai, MIAO Chunzhuang, GAO Xiaotao, LI Taotao, FENG Xiongbo
, Available online  , doi: 10.11858/gywlxb.20251102
[Abstract](199) PDF (20)
Abstract:
In order to improve the penetration ability of the shaped charge warhead to the water-containing composite structure, a truncated cone-ball combined cover was designed, and its jet forming motion law in the water medium and the damage performance to the water-containing composite structure were explored by numerical simulation. It is found that in the process of penetrating the water-containing composite structure, compared with the sub-hemisphere-sphere combined cover and the U-shaped-sphere combined cover, the truncated cone-sphere combined cover has a larger jet length L and a higher jet head velocity Vh. The cavity channel formed in the water medium and the radial expansion velocity Vr of the water medium are the smallest. The residual kinetic energy EK and the residual velocity VK of the jet after the breakdown target plate are the largest. The influence of structural parameters such as cone angle α, height h, side wall thickness a1 and top wall thickness a2 on the jet shape and penetration performance of the truncated cone-spherical combined cover is investigated by simulation technology, and the orthogonal optimization experiment is designed. It can be found that the primary and secondary order of the influence of these structural parameters on the jet penetration performance is the cone angle α > height h > side wall thickness a1 > top wall thickness a2. When α = 26 °, h = 22mm, a1 = 4.0mm, a2 = 3.2mm, the penetration performance of the truncated cone-sphere combined cover is better, and the residual kinetic energy EK of the jet when penetrating the aftereffect target is 136.2KJ. This study has certain reference value for the design of shaped torpedo warhead and the improvement of torpedo warhead damage power.
Numerical Simulation Study on the Effect of the Wave Shaper on Shaped Charge Jet Formation in Centrally-Apertured Liners
LIANG Zhouguang, FU Jianping, REN Kai, YANG Rui, SHI Junqing, WANG Bo, GAO Yueguang, CHEN Zhigang
, Available online  , doi: 10.11858/gywlxb.20251104
[Abstract](240) PDF (30)
Abstract:
To investigate the application of wave shapers in small-caliber shaped charges, the effects of aligning wave shaper parameters with those of center-holed liners on detonation product leakage and jet penetration performance were studied. Based on the regular oblique reflection theory of detonation waves, quantitative relationships between wave shaper parameters and detonation wave initial incident angles/pressure distributions at various positions on the liner surface were derived. Systematic analysis using LS-DYNA software revealed the influence patterns of wave shaper diameter and height on jet formation and penetration performance. The results indicate that adding wave shapers to center-holed liners effectively increases the collapse pressure on the liner, suppresses detonation product leakage, enhances energy utilization efficiency, and improves jet penetration performance. Jet penetration capability initially increases and then decreases with increasing wave shaper diameter. Wave shaper height exhibits a multi-extremum response effect on jet performance. The largest penetration depth of 158.17 mm into 45# steel targets was achieved with a wave shaper diameter of 6 mm and a height of 4 mm, representing a 17.21% improvement compared to structures without wave shapers. These findings offer valuable insights for designing small-caliber shaped charge warheads.
Experimental and Simulation Study on Impact Response of Submarine Optoelectronic Composite Cables
LUO Longqi, LI Kebin, ZHAO Yuantao, ZHOU Fenghua, ZHENG Yuxuan
, Available online  , doi: 10.11858/gywlxb.20251083
[Abstract](169) PDF (17)
Abstract:
This thesis aims to explore the impact resistance mechanical properties of Submarine Optical-Electrical Composite Cable (SOCC) under different working conditions. Firstly, a hammer impact experiment was carried out on SOCC to reveal the structural deformation characteristics of the outer armour under different variables and to record the impact evolution process and the maximum degree of concave deformation; secondly, a finite element simulation analysis was carried out on SOCC, and a comparison analysis was made with the experimental results; lastly, the deformation characteristics of SOCC under the influence of different parameters were explored. The results show that both the inner and outer armour undergo depression deformation, while the copper armour, copper conductor, and optical cable armour mainly show bending deformation, coupled with local depression deformation. With the increase of impact energy, the time required for metal components to reach the maximum deformation decreases, and the rebound becomes faster. The impact angle does not have a significant effect on the depression deformation of the inner and outer armour and produces significant damage to other internal components, of which the upper component has the most serious deformation damage. This paper is conducive to the evaluation of the dynamic performance of SOCC and provides a reference for the design of SOCC protection measures in engineering.
Ballistic and Compression after Impact Behaviors of Carbon/Basalt Fiber Hybrid Laminate
ZHANG Longfei, LIN Gaojian, YUAN Ye, ZHOU Zhipeng
, Available online  , doi: 10.11858/gywlxb.20251073
[Abstract](241) PDF (12)
Abstract:
This study examines the ballistic and compression-after-impact (CAI) performance of carbon fiber/basalt fiber (CF/BF) hybrid laminates with different CF/BF ratios. The results show that the basalt fiber significantly enhances the energy absorption capacity of the hybrid laminates. C-scanning and electron microscopy were used to investigate the damage mechanisms, providing insights into the improved performance. Although the compression strength of the original hybrid laminates monotonically decreases with the increasing BF content ratio, the residual compression strength measured from CAI tests for the damaged laminates shows a locally wavy trend due to the competition between the increased energy absorption capacity and decreased original compression strength. This study offers guidance for designing lightweight, impact-resistant composite structures.
Influence of Temperature on Mechanical Properties and Spall Damage of Invar36 Alloy
TANG Zeming, HU Jianbo, HU Changming, CHEN Sen
, Available online  , doi: 10.11858/gywlxb.20251057
[Abstract](276) PDF (40)
Abstract:
This study systematically investigated the effects of temperature on the spall behavior of Invar36 alloy through plate impact experiments and microstructural characterization techniques. Utilizing a single-stage light gas gun loading platform combined with a high-temperature heating device, the experiments measured free surface velocity profiles and spall strength variations in samples with different segregation orientations within the temperature range of 20°C to 300°C. Results demonstrate that the spall strength of Invar36 alloy exhibits a linear decrease with increasing temperature, with elevated temperatures significantly weakening its dynamic tensile resistance. Microstructural damage analysis reveals that at room temperature, voids nucleate and propagate along element segregation bands, while high-temperature damage concentrates at grain boundaries. Elevated temperatures reduce the constraining effect of segregation and facilitate material softening through thermally activated dislocation motion. The research elucidates the central role of temperature in governing spall strength and damage mechanisms, providing a theoretical foundation for failure-resistant design of Invar alloys under high-temperature impact conditions.
Research Progress on Dynamic Mechanical Response Characteristics of High-Velocity Particle Flow Impacting Multilayer Sandwich Composite Structures
ZHENG Wei, WANG Kunxuan, WANG Dengwang, LI Jun, GAO Yubo
, Available online  , doi: 10.11858/gywlxb.20251059
[Abstract](311) PDF (15)
Abstract:
Multi-layer sandwich composite structures have significant applications in impact protection. In particular, they demonstrate superior protective performance when subjected to impacts from explosive fragment particle clusters. Based on an analysis of the impact resistance and failure mechanisms of single-layer materials, this paper reviews the research progress regarding the dynamic mechanical response characteristics of composite structures under both single-particle and multi-particle impacts. The results indicate that metallic materials predominantly exhibit features such as plastic deformation, crack propagation, and localized thermal softening. By contrast, ceramics rapidly disperse impact energy due to their high hardness and propensity for brittle fracture. Meanwhile, fiber-reinforced composites achieve hierarchical energy dissipation through their continuous fiber network. Studies on multi-layer sandwich structures show that high-speed particle impacts on the target plate have been found to induce phenomena such as localized stress wave propagation, micro-crack formation, and interfacial delamination. The mechanisms underlying impact resistance in these structures are complex. However, current research primarily focuses on the impact resistance of structures under single-impact conditions. The protective mechanisms under multi-particle impacts remain unclear, and the employed research methods are relatively limited. Experimentally, approaches such as the modified split Hopkinson pressure bar (SHPB) apparatus are predominantly utilized to achieve high-speed loading of particle clusters. Nevertheless, issues regarding secondary impacts and velocity limitations in these experiments have yet to be effectively resolved. In numerical simulations, the smoothed particle hydrodynamics–finite element method (SPH-FEM) coupling approach remains the mainstream method for investigating particle cluster impacts. However, concerns regarding the accuracy of these models still warrant further investigation.
Tensile Fracture Characteristics and Dynamic Crack Evolution Law of Concrete
LIU Jinhao, LI Jinzhu, YAO Zhiyan, ZHANG Liwei
, Available online  , doi: 10.11858/gywlxb.20251046
[Abstract](220) PDF (18)
Abstract:
To investigate the tensile fracture characteristics and crack evolution mechanisms of concrete, Brazilian disc quasi-static splitting tests and falling weight impact tests were conducted. The crack propagation and mechanical responses were analyzed using the finite element-cohesive element method (FCEM). Experimental results demonstrated that under quasi-static loading, concrete discs exhibited tensile fracture with a primary crack penetrating along the loading direction at the disc center, accompanied by minor parallel secondary cracks. Crack propagation primarily occurred within the mortar matrix and along aggregate-mortar interfaces. The tensile performance of three-dimensional concrete discs exhibited significant enhancement with increasing thickness-diameter ratio. Under dynamic impact loading, specimens maintained a center-initiated fracture pattern, where the main crack propagated along the loading diameter, while triangular crushing zones formed at the edges in contact with testing apparatus. With increasing drop height, the specimens sequentially exhibited four distinct failure modes: no crack initiation, crack initiation without penetration, complete crack penetration, and severe fragmentation. High-speed photography quantified time-dependent crack lengths, demonstrating prolonged crack propagation durations at reduced drop heights. Numerical simulations revealed a nonlinear decreasing trend in crack initiation time versus drop height, with an empirical formula established to describe their relationship.