Phase Transition and Mechanical Property Modulation of Silicon Nitride at High Temperature and High Pressure
CUI Siwen, LI Shixin, MA Shuailing, LIAN Min, ZHAO Xingbin, TAO Qiang, ZHU Pinwen
2025, 39(12): 120101.
Review on Stress-Strain Rate Controllable Loading of Functionally Graded Materials
GUO Chengcheng, ZHANG Ruizhi, LIU Zhiqiang, HUANG Zihao, ZHANG Jian, LUO Guoqiang, SHEN Qiang
2025, 39(11): 110101.
Dynamical Mechanical Behaviors and Enhanced Ductility Mechanisms of Boron Carbide Based on Deep Potential Molecular Dynamics Simulations
LI Jun, SONG Jiahe, JI Wei, LIU Lisheng
2025, 39(11): 110102.
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
 doi: 10.11858/gywlxb.20251181
PDF(0)
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
 doi: 10.11858/gywlxb.20251176
PDF(0)
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
 doi: 10.11858/gywlxb.20251222
PDF(0)
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
 doi: 10.11858/gywlxb.20251227
PDF(1)
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
 doi: 10.11858/gywlxb.20251248
PDF(6)
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
 doi: 10.11858/gywlxb.20251221
PDF(2)
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
 doi: 10.11858/gywlxb.20251174
PDF(4)
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
 doi: 10.11858/gywlxb.20251224
PDF(2)
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
 doi: 10.11858/gywlxb.20251213
PDF(1)
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
 doi: 10.11858/gywlxb.20251056
HTML PDF(21)
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
 doi: 10.11858/gywlxb.20251054
HTML PDF(38)
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
 doi: 10.11858/gywlxb.20251193
PDF(1)
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
 doi: 10.11858/gywlxb.20251223
PDF(1)
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
 doi: 10.11858/gywlxb.20251198
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
 doi: 10.11858/gywlxb.20251191
PDF(3)
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
 doi: 10.11858/gywlxb.20251186
PDF(3)
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
 doi: 10.11858/gywlxb.20251173
PDF(5)
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
 doi: 10.11858/gywlxb.20251151
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
 doi: 10.11858/gywlxb.20251172
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
 doi: 10.11858/gywlxb.20251199
PDF(3)
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
 doi: 10.11858/gywlxb.20251164
PDF(4)
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
 doi: 10.11858/gywlxb.20251216
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
 doi: 10.11858/gywlxb.20251189
PDF(1)
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
 doi: 10.11858/gywlxb.20251155
PDF(1)
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
 doi: 10.11858/gywlxb.20251218
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
 doi: 10.11858/gywlxb.20251123
PDF(2)
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
 doi: 10.11858/gywlxb.20251169
PDF(6)
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
 doi: 10.11858/gywlxb.20251075
HTML PDF(22)
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
 doi: 10.11858/gywlxb.20251084
HTML PDF(27)
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
 doi: 10.11858/gywlxb.20251114
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
 doi: 10.11858/gywlxb.20251142
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
 doi: 10.11858/gywlxb.20251081
PDF(1)
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
 doi: 10.11858/gywlxb.20251140
PDF(4)
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
 doi: 10.11858/gywlxb.20251106
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
 doi: 10.11858/gywlxb.20251119
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
 doi: 10.11858/gywlxb.20251111
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
 doi: 10.11858/gywlxb.20251116
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
 doi: 10.11858/gywlxb.20251148
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
 doi: 10.11858/gywlxb.20251124
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
 doi: 10.11858/gywlxb.20251122
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
 doi: 10.11858/gywlxb.20251105
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
 doi: 10.11858/gywlxb.20251103
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
 doi: 10.11858/gywlxb.20251127
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
 doi: 10.11858/gywlxb.20251108
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
 doi: 10.11858/gywlxb.20251099
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
 doi: 10.11858/gywlxb.20251092
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
 doi: 10.11858/gywlxb.20251102
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
 doi: 10.11858/gywlxb.20251104
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
 doi: 10.11858/gywlxb.20251083
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
 doi: 10.11858/gywlxb.20251073
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
 doi: 10.11858/gywlxb.20251057
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
 doi: 10.11858/gywlxb.20251059
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
 doi: 10.11858/gywlxb.20251046
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.
Cover
2025, 39(12)  
PDF(0)
Abstract:
2025, 39(12): 1-2.  
HTML PDF(0)
Abstract:
Cover Feature
Phase Transition and Mechanical Property Modulation of Silicon Nitride at High Temperature and High Pressure
CUI Siwen, LI Shixin, MA Shuailing, LIAN Min, ZHAO Xingbin, TAO Qiang, ZHU Pinwen
2025, 39(12): 120101.   doi: 10.11858/gywlxb.20251060
HTML PDF(49)
Abstract:

Silicon nitride (Si3N4) ceramics are a novel category of structural ceramics, exhibiting both high reliability and economy value due to their unique physical and chemical properties. However, their strong covalent bonds make densification challenging, and phase transition control remains an issue. In this study, high temperature and high pressure (HTHP) sintering was utilized in combination with MgO-Y2O3 binary additives (mass ratio Si3N4∶MgO∶Y2O3 = 94∶3∶3) to achieve a synergistic sintering process involving both high pressure and liquid-phase formation. The effects of binary sintering additives on the sintering process, phase transition behavior, microscopic morphology and mechanical properties of Si3N4 under high pressure were systematically investigated by designing a two-layer comparative experimental assembly to ensure the same sintering temperature. The results show that the MgO-Y2O3 promotes the αβ phase transformation of Si3N4 through liquid phase formation, which reduces the phase transition onset temperature from 1800 ℃ to 1650 ℃. Furthermore, high pressure enhances grain rearrangement and densification. As a result, highly dense Si3N4 ceramics with an optimal Vickers hardness up to (24.5±1.9) GPa were successfully synthesized. These findings provide an effective strategy for preparing high-performance Si3N4 ceramics and hold significant implications for advancements in physics and material science.

Physical Property and Structure
Investigation of Mechanical Behavior in Nanocrystalline Palladium under High Pressure
LIU Boyu, LIU Jingyi, ZHUANG Yukai, WANG Qiming, ZHANG Youjun
2025, 39(12): 121101.   doi: 10.11858/gywlxb.20251133
HTML PDF(12)
Abstract:

The investigation of mechanical response characteristics of nanocrystalline metallic materials under extreme high-pressure conditions possesses significant scientific importance and engineering value. Using a diamond anvil cell combined with synchrotron radiation X-ray diffraction techniques, the mechanical behavior of palladium (Pd) with an average grain size of approximately 10 nm under static high pressure was studied. Within the investigated pressure range (0−111 GPa), the crystal structure of palladium remained stable. Analysis of diffraction peak positions and full width at half maximum (FWHM) at each pressure point enables determination of unit cell volume, grain size, and microscopic strain under high-pressure conditions. Fitting with the third-order Birch-Murnaghan equation yields bulk modulus of 288 GPa (hydrostatic) and 290 GPa (non-hydrostatic), and the yield strength is approximately 20 GPa. In addition, by integrating existing literatures, this study systematically explored the influence of size effects on mechanical properties. The yield strength of Pd progressively increases with decreasing grain size, exhibiting a 300% enhancement compared to Pd nanofibers. These findings provide crucial data for the structural design and application of nanocrystalline Pd under extreme conditions.

Material Synthesis
Sintering and Characterizing of WC-5Co Cemented Carbide under High Temperature and High Pressure
CHEN Chunhua, GUO Ruiang, HE Duanwei, WEI Xiuyan, HU Zuguang, YANG Jianyun
2025, 39(12): 123101.   doi: 10.11858/gywlxb.20251067
HTML PDF(9)
Abstract:

As an important type of powder metallurgy products, tungsten carbide-cobalt (WC-Co) cemented carbides are widely applied in numerous industrial fields. Traditional sintering methods suffer from high sintering temperatures, long sintering times, and grain coarsening, which make it difficult to meet the requirements for the preparation of high-strength fine-grained cemented carbides. In this study, the sintering behavior of the sub-micron WC-5Co (with a mass fraction of cobalt of 5%) cemented carbides was investigated using the high-temperature and high-pressure (HTHP) method. The experimental results show that the sintered body can achieve a relative density of 99.4% and a Vickers hardness value of (24.0±0.3) GPa under the treatment conditions of 5.0 GPa and 1 250 ℃, which is significantly better than that of the same composition of cemented carbides prepared by spark plasma sintering (SPS) and vacuum sintering methods. In comparison with the commercial cemented carbides with the same initial powder, the average grain size of WC in the sintered samples under high-temperature and high-pressure conditions is reduced by about 36%, and the Vickers hardness is increased by about 10%. In comparison with traditional sintering methods, the HTHP sintering method effectively reduces the sintering temperature, shortens the sintering time, and enhances sintering efficiency. Moreover, high pressure can effectively inhibit the abnormal growth of the WC grains during solid-state sintering. During liquid-phase sintering (at 1 450 ℃), an increase in sintering pressure results in significant refinement of the WC grains. When the sintering pressure is increased from 3.0 GPa to 5.0 GPa, the Vickers hardness of the samples is increased by about 7%.

Theoretical Study on the Polymerization Mechanism of Hydrogen-Doped Carbon Monoxide under High Pressure
LI Xiuyuan, PENG Zihuan, JIANG Chongwen, HUANG Zhihong, LI Nan
2025, 39(12): 123102.   doi: 10.11858/gywlxb.20251088
HTML PDF(12)
Abstract:

Carbon monoxide (CO), as a prototypical low-Z molecular system, can polymerize under high pressure to form polymeric carbon monoxide (p-CO). The polymerization mechanisms and structures are of fundamental importance for understanding pressure-induced bonding and exploring novel functional materials. However, progress in this field has been hindered by two major challenges: the high-pressure requirements for CO and the metastable property of p-CO at ambient pressure. Recent studies have shown that hydrogen (H2) doping can facilitate the polymerization of CO, but the polymerization mechanisms and structures are still poorly understood. In this work, molecular dynamics simulations were performed to investigate the influence of H2 on the polymerization progress of CO. The results reveal that a doping ratio of 10% can optimally reduce the polymerization pressure of CO. At 3–4 GPa, H2 physically induces the dimerization reaction of CO. At 5 GPa, the chemical inertness of H2 inhibits further polymerization of CO. When the pressure reaches 10 GPa, H2 participates in the polymerization reaction, forming C―H and O―H bonds. Finally, the polymerization produces a disordered three-dimensional network structure (p-CO/H) dominated by C―C and C=O bonds.

Dynamic Response of Matter
Dynamic Mechanical Properties and Damage Evolution in Granite under Coupled High-Temperature-Impact Cyclic Loading
GUO Hao, ZUO Xuchao, ZHONG Kangjie, WU Jun, LIU Lei
2025, 39(12): 124101.   doi: 10.11858/gywlxb.20251091
HTML PDF(14)
Abstract:

To study the influence of high-temperature-impact dual-cycle cumulative damage on the dynamic mechanical characteristics of granite, specimens with a height-to-diameter ratio of 0.8 were subjected to high-temperature cycling treatments at 100, 300, and 500 ℃ for 2, 4, and 6 cycles, respectively. The P-wave velocities before and after treatments were measured. Constant-amplitude cyclic impact tests under gas pressures of 0.25, 0.30, and 0.35 MPa were conducted using a split Hopkinson pressure bar (SHPB) system. The effects of temperature, high-temperature cycles, impact gas pressure, and impact cycles on the dynamic behavior of granite were systematically investigated. Based on Lemaitre’s continuum damage constitutive model and the strain equivalence principle, cumulative damage factors for high-temperature and impact cycles were defined, and the critical dual-cycle cumulative damage factor was analyzed. Results indicate that with increasing initial high-temperature cumulative damage and impact gas pressure, the crack morphology of granite during the first impact evolved from single cracks to complex crack networks with increased connectivity, while peak stress decreased and peak strain increased. Significant variations in peak stress and strain were observed between the first and final impacts during cyclic loading. The influence of factors on the cumulative impact damage under high-temperature cycling followed the order: impact gas pressure>temperature>number of high-temperature cycles. The critical high-temperature-impact cumulative damage factor was determined to range between 0.625 and 0.676. These findings provide theoretical support for safety assessment in deep underground resource extraction.

Impact Dynamic Characteristics and Damage Features of Early-Age Concrete-Mudstone Composite
PAN Chong, XIE Quanmin, SUN Jinshan, ZHOU Hui, PANG Chaolai, MA Jun
2025, 39(12): 124102.   doi: 10.11858/gywlxb.20251094
HTML PDF(9)
Abstract:

To investigate the dynamic mechanical properties of early-age concrete-mudstone composites under impact loading, split Hopkinson pressure bar (SHPB) tests integrated with a high-speed camera were conducted on composite specimens with curing ages of 1, 3, and 7 d. Digital image correlation (DIC) technology was employed to analyze the evolution of displacement and strain fields, systematically revealing the dynamic damage and failure characteristics of the composites. The test results indicate that as the strain rate increases, the composite specimens exhibit significant strain rate dependence across all curing ages, and their dynamic strength growth follows a logarithmic function model. The energy dissipation density increases linearly with incident energy. DIC measurements show that the maximum surface displacements of the 1, 3, and 7 d specimens are 1.564, 1.196, and 0.924 mm, respectively, with corresponding maximum strains of 1.886%, 1.352%, and 1.184%. This study elucidates the dynamic mechanical behavior and damage-failure mechanisms of early-age concrete-mudstone composites under impact loading, providing a theoretical foundation for damage prevention and control in surrounding rock-initial support structures during tunnel blasting construction.

High Pressure Applications
Discrete Element Analysis on the Influence of Block Shape and Spatial Arrangement on Shielding Performance
LUO Yuting, ZHAO Tingting, JU Kaixuan, WANG Zhiyong
2025, 39(12): 125101.   doi: 10.11858/gywlxb.20251087
HTML PDF(12)
Abstract:

Shelter layers, recognized as a critical component in modern military defense systems, are designed to protect strategically important rear assets. Given the prevalent application of block stones as primary filling materials in these protective structures, comprehensive investigations into their fundamental ballistic resistance mechanisms and practical performance optimization have become increasingly crucial for defense engineering. Therefore, the dynamic fragmentation behavior of block stones under high-velocity projectile impact was systematically simulated using discrete element spherical particles and bonded particle model (BPM) in this paper. The process of rigid projectile vertical penetration into densely packed block stone structure was numerically simulated, with the influences of block stone particle size, shape, and spatial arrangement characteristics on its anti-penetration performance being explored. The results revealed two primary energy dissipation mechanisms governing the shelter layer performance: over 90% of the kinetic energy of the projectile was dissipated through the combined effects of inter-block collisions and sliding friction during the penetration process. Meanwhile, the number of break block stones is found to be negatively correlated with the block stone particle size and positively correlated with the aspect ratio (long-to-short axis) of the block stones. Moreover, when single-sized block stones are used in a multi-layer staggered arrangement, the penetration depth of the projectile is primarily determined by the magnitude of the peak penetration resistance. Specifically, the maximum penetration resistance and the minimum penetration depth are exhibited by the scenario with circular block stones of 120 mm particle size. However, when a layered arrangement with block stone particle size gradiently decreasing along the direction of the impact face is employed, the overall blast-resistant effect of the structure is not effectively enhanced. Finally, the results of the study can provide a reference for understanding the underlying mechanism of the shielding effect of block stone structure.

Experimental Study on Ignition Time of Industrial Electronic Detonator Ignition Head under Different Ignition Voltages
GUAN Jiajia, LI Hongwei, YANG Lin, LIANG Hao, LIU Yanhong, LU Shaobao, ZHU Fuyuan
2025, 39(12): 125201.   doi: 10.11858/gywlxb.20251101
HTML PDF(10)
Abstract:

In response to the practical engineering problem that the firing capacitor of electronic detonators may lose power after being impacted, which subsequently lead to insufficient firing energy and unreliable detonation, the firing time of metal bridge and bridge-wire ignition heads under different firing voltages were studied through the high-speed schlieren testing system. The influence of the initial energy input to the ignition head on the thermal decomposition stage, flame growth period, and flame duration was clarified, and the relationship between the total ignition time and the ignition voltage and ignition energy for two types ignition heads was obtained. The results show that within the range of 13−21 V, the rate of change of the total firing time of the ignition heads with voltage decreases first and then increases. As the firing voltage increases, the time reduction rates of the thermal decomposition, flame growth period, and flame duration of the metal bridge model ignition heads are 66.2%, 76.6%, and 15.0% respectively. The time reduction rates of the three stages of the bridge-wire ignition heads are 28.0%, 39.2%, and 30.0% respectively, and the firing time of each stage of the metal bridge model ignition head is shorter than that of the bridge-wire. When the ignition energy is between 1.9 and 4.9 mJ, the firing consistency and accuracy of the metal bridge model and bridge-wire ignition head will be affected. When the ignition energy is less than 1.9 mJ, the metal bridge mode ignition heads will fail to fire due to insufficient firing energy. This study provides a basis for the design of the firing margin of electronic detonators, which will reduce the misfire rate of electronic detonators in small-hole blasting.

Effect of Circumferential Stress on Energy Evolution Mechanism of Rockburst under True Triaxial Unilateral Unloading Conditions
ZHANG Huaxu, HUANG Xin, GUO Jiaqi, SUN Feiyue, ZHU Zihui
2025, 39(12): 125301.   doi: 10.11858/gywlxb.20251066
HTML PDF(15)
Abstract:

In order to study the influence of hoop stress on the energy evolution in the process of rockburst catastrophe, a new true triaxial rockburst test system was used to carry out the rockburst simulation test of single-side rapid unloading, three-five-side force and vertical continuous loading under different circumferential stresses. The rockburst failure modes of granite samples under different hoop stresses were analyzed. Combined with the principle of energy conservation, the evolution law of individual energy in the process of rockburst catastrophe was revealed. The results show that there is an obvious energy competition evolution mechanism between dissipated energy and elastic strain energy under different circumferential stresses. The circumferential stress will significantly affect the damage degree and distribution range of the rock sample. The unloading surface of the rock sample with a circumferential stress of 178.992 MPa has the deepest damage degree. Under the action of high circumferential stress, the elastic strain energy in the rock sample releases faster after the peak point, and the development of rockburst has short-term characteristics. The conversion rate of dissipated energy is proportional to the circumferential stress, and the conversion rate of elastic strain energy is inversely proportional to the circumferential stress. From the absolute value of energy, the increase of circumferential stress will obviously increase the accumulation of elastic strain energy and the release of dissipated energy. The total energy conversion rate of the rock sample is greater than the elastic energy conversion rate, which in turn exceeds the dissipated energy conversion rate. All three parameters exhibit positive correlations with circumferential stress. The increase of the circumferential stress will obviously accelerate the conversion rates of total energy, elastic energy and dissipated energy.

Blasting Crater Test Law under Different Resistance Lines Based on DEM-PBM Method
HE Defu, WANG Kai, LI Yucheng, LI Zhuo, ZHANG Zhiyu, HUANG Yonghui, LUO Yi
2025, 39(12): 125302.   doi: 10.11858/gywlxb.20251086
HTML PDF(14)
Abstract:

In open-pit blasting project, the minimum resistance line (W) determines the dominant direction of blasting energy release and the maximum direction of rock velocity. This parameter is the core factor controlling the blasting action range, crushing effect and throwing path of the rock. Based on discrete element method-population balance model (DEM-PBM) coupling numerical simulation method of blasting crater, field blasting crater test was carried out through multi-scale modeling and test data collaborative analysis. The rock breaking law of blasting crater test under different minimum resistance lines was studied in order to optimize blasting parameters and improve blasting efficiency and safety. Based on the comparison of numerical simulation results, the blasting expansion area at 20 ms increases by 70.33% compared with that at 5 ms when W=1.05 m, and the blasting expansion area increases by 11.42% when W=1.85 m. By comparing the throwing effect of fragmentation under different minimum resistance lines, the throwing effect is the best when W=1.05 m. According to the results of field blasting crater test, when the length of the resistance line increases, the volume of the blasting crater becomes smaller, and the volume of the blasting crater is the largest when W=1.05 m, which is similar to the blasting cavity expansion law obtained through numerical simulation. The fractal law of fragmentation size in field blasting was analyzed. When W=1.05 m and the depth of blast hole is 1.2 m, more rock blocks are thrown, the fragmentation distribution is moderate, and the throwing effect of broken blocks is the best. This study provides reference for the optimization of field drilling and blasting parameters.

Damage and Fractal Analysis of Double-Hole Blasting
HUO Fei, LI Xianglong, WANG Jianguo, CUI Guangjiu, HOU Meng, KONG Xiangjie, SUN Fangyixiong
2025, 39(12): 125303.   doi: 10.11858/gywlxb.20251093
HTML PDF(13)
Abstract:

To investigate the damage and failure characteristics of the double-hole blasting medium, theoretical analysis and model experiments are conducted in this study. Based on the analytical solution of the elastic plane strain problem in the double-hole blasting model, a theoretical model is established to study the evolution of the dynamic stress field during double-hole blasting. The damage of double-hole blasting with different blast hole spacing is studied by the model experiment; the damage in different regions around the blast holes is analyzed by partitioned research and fractal dimension quantification. The study shows that as the blast hole spacing increases, the super position effect of stress waves decreases, the radius of the crushing zone and the average length of the main crack gradually increase, the number of cracks decreases first and then increases. Additioanlly, the increase of blast hole spacing lead to gradual decrease of the damage along the direction of the blast hole connection, while the damage perpendicular to the direction of the blast hole connection gradually increases. Small hole spacing is conducive to the penetration of cracks between holes and promotes cracks expansion along the direction of the blast hole connection. With increasing blast hole spacing, the damage variable in the left and right regions of the specimen gradually increases, and reaches a minimum value in the central region when the blast hole spacing is 50 mm. The damage variables in area Ⅰ and area Ⅱ first decrease and then increase, while the damage variable in area Ⅲ gradually decreases. The damage variable distribution in the polar coordinate system shows that area Ⅰ exhibits a uniform damage pattern, and area Ⅱ gradually transitions from an elliptical to a circular distribution. Based on the relationship between the fractal dimension of the damage area and the damage variable, a fractal damage model for double-hole blasting of poly(methyl methacrylate) (PMMA) material is constructed.

Experimental Study of the Effect of Shear Stress on Phase Transition in c-Axis CdS Single Crystal under Dynamic Loading
TANG Zhi-Ping, Gupta Y M
1989, 3(4): 290-297 .   doi: 10.11858/gywlxb.1989.04.005
[Abstract](14215) [PDF 8643KB](2276)
Abstract:
For a long time, the problem whether shear stress affects the phase transition initial pressure is not well solved. Duvall and Graham suggested that cadmium sulfide (CdS) crystal could be used to study the effect of shear stress on the initial pressure of phase transition in c-axis CdS single crystal specimens under high velocity impact systematically. The axial stress of initial phase transition measured is T=(3.250.1) GPa, corresponding to a mean pressure pT=(2.290.07) GPa, which agrees the value 2.3 GPa of static results quite well within the experimental error. The shear stress in this case, T=0.72 GPa, is as high as 31.5% of the mean pressure. This result shows that the mechanism of phase transition may be assumed only to relate to a critical mean pressure or critical thermodynamic state, and the effect of shear stress can be ignored.
Flattening of Cylindrical Shells under External Uniform Pressure at Creep
Shesterikov S A, Lokochtchenko A M
1992, 6(4): 247-253 .   doi: 10.11858/gywlxb.1992.04.002
[Abstract](10053) [PDF 2836KB](2017)
Abstract:
Experimental studies of the deformation of cylindrical shells under creep to fracture conditions are described in this paper. Analyses of three series of test shells are given and experimental and theoretical results are compared with each other.
The Generation of 90 GPa Quasi-Hydrostatic Pressures and the Measurements of Pressure Distribution
LIU Zhen-Xian, CUI Qi-Liang, ZOU Guang-Tian
1989, 3(4): 284-289 .   doi: 10.11858/gywlxb.1989.04.004
[Abstract](16114) [PDF 6073KB](2682)
Abstract:
Quasi-hydrostatic pressures up to 90 GPa were obtained at room temperature in the diamond cell by using solid argon as pressure medium. The pressure distribution was determined by measuring the special shift of the R1 line of ruby at different positions within the sample chamber. Experimental results showed that the pressure differences (p) between the pressures at each point within the chamber and the mean pressure (p) were very small, ratios of p/p were less than 1.5% when below 80 GPa. The shape of ruby R lines at 90 GPa is similar to that at ambient pressure. Thus, quasi-hydrostatic pressure near 100 GPa can be obtained by using solid argon as pressure medium. Moreover, the red shifts with pressures of the peak positions at 14 938 and 14 431 cm-1 in ruby emission spectra, were also examined. It concluded that the line, 14 938 cm-1, can be adopted in the pressure calibration.
A Study on Calculation of the Linear Thermal Expansion Coefficients of Metals
ZHENG Wei-Tao, DING Tao, ZHONG Feng-Lan, ZHANG Jian-Min, ZHANG Rui-Lin
1994, 8(4): 302-305 .   doi: 10.11858/gywlxb.1994.04.010
[Abstract](16883) [PDF 1350KB](1189)
Abstract:
Based on the expression of pressure at temperature T and in terms of the universal equation of state Debye model and the thermodynamic relations, a general expression for the calculation of the linear thermal expansion coefficients of metals is obtained. This formula applied to the calculation of Al, Cu, Pb. Calculated results are in good agreement with the experiments.
Development of Large Volume-High Static Pressure Techniques Based on the Hinge-Type Cubic Presses
WANG Hai-Kuo, HE Duan-Wei, XU Chao, GUAN Jun-Wei, WANG Wen-Dan, KOU Zi-Li, PENG Fang
2013, 27(5): 633-661.   doi: 10.11858/gywlxb.2013.05.001
[Abstract](15003) [PDF 12118KB](1133)
Abstract:
The large volume press (LVP) becomes more and more popular with the scientific and technological workers in the high pressure area, because it could generate relatively higher pressure, provide better hydrostatic pressure and could be utilized in conjunction with in situ X-ray diffraction, neutron diffraction and ultrasonic measurement. There have been generally two LVP techniques to generate high-pressure: the double-anvil apparatus and the multi-anvil apparatus (MAA). Hinge-type cubic presses, as the main apparatus in china, have been widely used in the fields of both scientific research and diamond industry. However, for a long time past, the maximum pressure using the conventional one-stage anvil system for hinge-type cubic press is about 6 GPa, and the techniques about two-stage apparatus (octahedral press) that could generate pressure exceed 20 GPa is blank in our country. To a certain extent, the backwardness of the LVP technology in china restricts the development of high pressure science and related subjects. In recent years, we designed two kinds of one-stage high pressure apparatus and the two-stage apparatus based on hinge-type cubic-anvil press, the one-stage high pressure apparatus and the two-stage apparatus using cemented carbide as anvils could generate pressures up to about 9 GPa and 20 GPa respectively. This article mainly reviews the mechanics structure, design of cell assembly, pressure and temperature calibration, design and preparation of the sintered diamond anvils and pressure calibration to 35 GPa using sintered diamond as two-stage anvils about the one-stage high pressure apparatus and the two-stage apparatus designed in our laboratory.
The Failure Strength Parameters of HJC and RHT Concrete Constitutive Models
ZHANG Ruo-Qi, DING Yu-Qing, TANG Wen-Hui, RAN Xian-Wen
2011, 25(1): 15-22 .   doi: 10.11858/gywlxb.2011.01.003
[Abstract](18268) [PDF 689KB](1101)
Abstract:
The analyzed and calculated results indicate that the concrete failure strength will decrease under higher hydrostatic pressure, when the original failure parameters of HJC and RHT models implemented in LS-DYNA and AUTODYN are adopted. A new method is introduced which using the characteristic strength of concrete to confirm the modified failure parameters of HJC and RHT models. The same physical experiment of concrete penetration was simulated using the modified HJC and RHT failure parameters respectively, and the numerical results demonstrated that the RHT model matched the experiments much better. But the numerical results with the HJC modified failure parameters were not enough satisfied, because the third invariant of the deviated stress tensor was not considered in the HJC model.
Research on Deformation Shape of Deformable Warhead
GONG Bai-Lin, LU Fang-Yun, LI Xiang-Yu
2010, 24(2): 102-106 .   doi: 10.11858/gywlxb.2010.02.004
[Abstract](9019) [PDF 1765KB](395)
Abstract:
Basing on the detonation theory, the structure of the deformable warhead was simplified to be double layer cylindrical shells under the detonation. Plastic hinges were introduced into the loading section of the shell, which contacted with the deforming charge, and the deforming charge was divided into small segments accordingly. Loading and movement of these segments were analyzed. Deforming shape of the cylindrical shell under the loading with equal distribution was bulgy, and the displacement of shell segments was obtained. Deforming charge with different thickness, according to the displacement of the segment, was set up to realize the same displacement of the shell segments on the loading direction. The D-shape was achieved theoretically, and the shape of deforming charge was designed accordingly. Numerical simulation validated the feasibility of the designed plan. The results indicate that the deformable warhead with the new-designed deforming charge can realize the D-shape.
Recent Progresses in Some Fields of High-Pressure Physics Relevant to Earth Sciences Achieved by Chinese Scientists
LIU Xi, DAI Li-Dong, DENG Li-Wei, FAN Da-Wei, LIU Qiong, NI Huai-Wei, SUN Qiang, WU Xiang, YANG Xiao-Zhi, ZHAI Shuang-Meng, ZHANG Bao-Hua, ZHANG Li, LI He-Ping
2017, 31(6): 657-681.   doi: 10.11858/gywlxb.2017.06.001
[Abstract](11988) [FullText HTML](5035) [PDF 2527KB](5035)
Abstract:

In the last 10 years or so, nearly all major Chinese universities, schools and research institutes with strong Earth science programs showed strong interest in developing a new research branch of High-Pressure Earth Sciences.As a result, many young Chinese scientists with good training from the universities in the west countries were recruited.This directly led to a fast growing period of about 10 years for the Chinese high-pressure mineral physics research field.Here we take the advantage of celebrating the 30th anniversary of launching the Chinese Journal of High Pressure Physics, and present a brief summary of the new accomplishments made by the Chinese scientists in the fields of high-pressure mineral physics relevant to Earth sciences.The research fields include:(1) phase transitions in the lower mantle; (2) high spin-low spin transitions of iron in lower mantle minerals; (3) physical properties of the Earth core; (4) electrical measurements of rocks; (5) electrical measurements of minerals; (6) elasticity of minerals (especially equation of states); (7) high-pressure spectroscopic studies; (8) chemical diffusions in minerals; (9) ultrasonic measurements under high pressure; (10) physical properties of silicate melts; (11) geological fluids.In sum, the last 10 years have seen a rapid development of the Chinese high-pressure mineral physics, with the number of scientific papers increasing enormously and the impact of the scientific findings enhancing significantly.With this good start, the next 10 years will be critical and require all Chinese scientists in the research field to play active roles in their scientific activities, if a higher and advanced level is the goal for the Chinese mineral physics community.

Modification of Tuler-Butcher Model with Damage Influence
JIANG Dong, LI Yong-Chi, GUO Yang
2009, 23(4): 271-276 .   doi: 10.11858/gywlxb.2009.04.006
[Abstract](11393) [PDF 402KB](871)
Abstract:
A modificatin of Tuler-Butcher model including damage influence was presented, which was incorporated into a hydrodynamic one-dimensional finite difference computer code, to simulate the process of spall fracture of 45 steel and Al-Li alloy. The calculated results are in good agreement with experimental data, and shows the correctness of the model.
Experiment and Numerical Simulation of Cylindrical Explosive Isostatic Pressing
CHEN Lang, LU Jian-Ying, ZHANG Ming, HAN Chao, FENG Chang-Gen
2008, 22(2): 113-117 .   doi: 10.11858/gywlxb.2008.02.001
[Abstract](14290) [PDF 1180KB](897)
Abstract:
The experiments of cylindrical explosive isostatic pressing were carried out. The internal temperatures in pressed explosives were measured by thermocouples. A thermal/structural coupled model of the explosive isostatic pressing was set up. The numerical simulations of cylindrical explosive were conducted. The calculated pressures and temperatures in explosives were given. The deformations,pressures and temperatures distribution were analyzed. The calculated results indicated that each surface center of the cylindrical explosive was sunken by isostatic pressing. During the isostatic pressing of cylindrical explosive, the internal temperature of the explosives increases, and the temperature and pressure are not uniform.
Perimeter-Area Relation of Fractal Island
LONG Qi-Wei
1990, 4(4): 259-262 .   doi: 10.11858/gywlxb.1990.04.004
[Abstract](16347) [PDF 1508KB](2413)
Abstract:
The relationship of perimeter with area (P/A relation) of fractal island is discussed. It is shown that Mandelbrot's fractal relation between Koch perimeter and area does not hold in the island with finite self-similar generations. This might be the reason why the fractal dimension measured with P/A relation varied with the length of yardstick in previous work.
Long-Distance Flight Performances of Spherical Fragments
TAN Duo-Wang, WEN Dian-Ying, ZHANG Zhong-Bin, YU Chuan, XIE Pan-Hai
2002, 16(4): 271-275 .   doi: 10.11858/gywlxb.2002.04.006
[Abstract](14916) [PDF 2450KB](943)
Abstract:
Using two-stage light gas gun and laser technique for velocity easurement, we studied the long-distance flight performances of spherical fragments with different materials and different diameters. The flight distance is 60~120 m, and the initial velocity is 1.2~2.2 km/s. The experimental results show that: (1) the velocity attenuation coefficient of spherical fragment is constant, and (2) the air drag coefficient is slightly affected by the initial velocity of spherical fragment, the air drag coefficient is a linear function of initial velocity.
Application Research on JWL Equation of State of Detonation Products
ZHAO Zheng, TAO Gang, DU Chang-Xing
2009, 23(4): 277-282 .   doi: 10.11858/gywlxb.2009.04.007
[Abstract](16428) [PDF 365KB](992)
Abstract:
By investigating the JWL equation of state of detonation products of condensed explosive, we present a method to determine JWL parameters by fitting. This approach does not require cylinder test and is more economical, secure, convenient and accurate than existing methods. Using this method, four kinds of common explosive, e.g., TNT, C-4, PETN and HMX have been studied. By comparing to the p-V curve of JWL equation of state given by cylinder test, we showed that the fitting has a high precision and meets the need of explosion mechanics application.
Design and Temperature Calibration for Heater Cell of Split-Sphere High Pressure Apparatus Based on the Hinge-Type Cubic-Anvil Press
CHEN Xiao-Fang, HE Duan-Wei, WANG Fu-Long, ZHANG Jian, LI Yong-Jun, FANG Lei-Ming, LEI Li, KOU Zi-Li
2009, 23(2): 98-104 .   doi: 10.11858/gywlxb.2009.02.004
[Abstract](15616) [PDF 4054KB](937)
Abstract:
A new type of heater cell for the split-sphere high pressure apparatus based on the hinge-type cubic-anvil press was reported. This heating apparatus has the advantages of being simple, low cost, fast temperature rising, good heat insulation, and the temperature signal can be easily extracted. Carbon tube was used as a heating element for side-heating in our experiments. The size of the sample in the cell can reach 3 mm in diameter, and 7 mm in height. The relationship between the heating electric power and cell temperature was calibrated with Pt6%Rh-Pt30%Rt thermocouples under different pressures. The experimental results indicate that the temperature can reach 1 700 ℃ under the oil hydraulic pressure of 40 MPa (cell pressure is about 10 GPa).The temperature can keep stable for more than 2 h under a fixed power.
Factors Analysis of Debris Cloud's Shape of Hypervelocity Impact
TANG Mi, BAI Jing-Song, LI Ping, ZHANG Zhan-Ji
2007, 21(4): 425-432 .   doi: 10.11858/gywlxb.2007.04.016
[Abstract](14051) [PDF 1599KB](796)
Abstract:
The numerical simulations of hypervelocity impact of Al-spheres on bumper at normal are carried out using the smoothed particle hydrodynamics (SPH) technique. The simulation results are compared with experimental results, and the simulated hole diameters of bumper and debris cloud are well consistent with experimental results. The effect of impact velocity, bumper thickness, projectile diameter, materials, shape of projectile, interval on produced debris cloud are further analyzed. Regarding the length and diameter as index, orthogonal design method is applied to analyze the primary and secondary relations on the debris cloud's index of the three factors, that is impact velocity, bumper thickness and projectile diameter. The results indicate that bumper thickness is the main influence factor of debris cloud's length while projectile diameter is the main influence factor of debris cloud's diameter.
Design of the Sample Assembly for Ultrasonic Measurement at High Pressure and 300 K in Six-Side Anvil Cell
WANG Qing-Song, WANG Zhi-Gang, BI Yan
2006, 20(3): 331-336 .   doi: 10.11858/gywlxb.2006.03.019
[Abstract](11645) [PDF 411KB](722)
Abstract:
We introduced briefly the principle of design of sample assembly for ultrasonic measurements at high pressure, and designed a new kind of sample assembly to measure the isothermal compression of Al and Cu at 300 K. Ideal quasi-hydrostatic loading was realized, and high-quality ultrasonic signals were obtained under high pressure. It was indicated that the design of sample assembly was reasonable. We analyzed in brief main uncertainty of ultrasonic measurement in six-side anvil cell at 300 K.
Detonation Shock Dynamics Calibration of JB-9014 Explosive at Ambient Temperature
TAN Duo-Wang, FANG Qing, ZHANG Guang-Sheng, HE Zhi
2009, 23(3): 161-166 .   doi: 10.11858/gywlxb.2009.03.001
[Abstract](14831) [PDF 794KB](917)
Abstract:
Detonation shock dynamics (DSD) is an approximation to the reactive Euler equations that allows numerically efficient tracking of curved detonation waves. The DSD parameters are the velocity curvature relation and the boundary angle. A computer code was developed to facilitate the calibration of these parameters for JB-9014 insensitive high explosive using the generalized optics model of DSD. Calibration data were obtained from measurements of the detonation velocities and fronts in JB-9014 rate sticks at ambient temperature, with diameters of 10~30 mm. The steady state detonation velocities and fronts predicted by these DSD parameters are in very good agreement with experiment.
Shock Wave Physics: The Coming Challenges and Exciting Opportunities in the New Century-Introduction of the 12th International Conference of Shock Compression of Condensed Matter (SCCM-2001)
GONG Zi-Zheng
2002, 16(2): 152-160 .   doi: 10.11858/gywlxb.2002.02.012
[Abstract](14854) [PDF 500KB](944)
Abstract:
The 12th Biennial International Conference of the APS Topical Group on Shock Compression of Condensed Matter (SCCM-2001) was introduced. Papers presented in SCCM-2001 were surveyed and the recent progresses on shock compression of condensed matter were retrospected. The basic paradigms and the great achievements of the physics and mechanics of condensed matter at high dynamic pressure and stress were surveyed and revaluated. The coming challenges and exciting opportunities of shock wave physics in the 21 century were prospected.
The Constitutive Relationship between High Pressure-High Strain Rate and Low Pressure-High Strain Rate Experiment
CHEN Da-Nian, LIU Guo-Qing, YU Yu-Ying, WANG Huan-Ran, XIE Shu-Gang
2005, 19(3): 193-200 .   doi: 10.11858/gywlxb.2005.03.001
[Abstract](11494) [PDF 416KB](864)
Abstract:
It is indicated that the constitutive equations at high strain rates proposed by Johnson-Cook(J-C), Zerilli-Armstrong (Z-A) and Bodner-Parton (B-P) collapse the data of flow stress in compression, tension, torsion, and shear into simple curve with the scalar quatities 'effective' stress and 'effective' strain, however, the collapsed data of flow stress did not include the data in the planar shock wave tests. The SCG constitutive equation proposed by Steinberg et al for the planar shock wave tests is discussed, which describes the coupled high pressure and high strain rate effects on the plastic deformation of materials. Basing on the recent experiments at elevated temperatures and high strain rates and the shear strength measurements during shock loading, the flow stress for tungsten at high pressure and high strain rates is estimated with J-C and SCG constitutive equations, respectively. It is concluded that the J-C, Z-A and B-P constitutive equations may not be appropriate to describe the plastic behavior of materials at high pressure and high strain rates, comparing with SCG constitutive equation. It is emphasized that the physical background of the constitutive equation at high pressure and high strain rates is different from that at low pressure and high strain rates.
Experimental Study on the Damage Effect of Compound Reactive Fragment Penetrating Diesel Oil Tank
XIE Chang-You, JIANG Jian-Wei, SHUAI Jun-Feng, MEN Jian-Bing, WANG Shu-You
2009, 23(6): 447-452 .   doi: 10.11858/gywlxb.2009.06.008
[Abstract](13305) [PDF 5649KB](1009)
Abstract:
Two new kinds of compound reactive fragments were designed and prepared, and the penetration tests of the compound reactive fragments against oil tank with diesel oil were performed. The compound reactive fragment is composed of shell, bare reactive fragment and coping. Bare reactive fragments prepared by high-temperature sintering in a vacuum container have two kinds of formulations, one is mixed aluminum powder with PTFE, another is mixed titanium powder with PTFE. Fragments were fired using 12.7 mm ballistic gun, and the penetration process against oil tank were recorded by high-speed camera. The experimental results show that two kinds of compound reactive fragment can penetrate through 6 mm thick oil tank and have obvious ignition effects. Comparing with inert fragments, compound reactive fragments have better capability of penetration and ignition.