2026, 40(7): 070101.
Development of a 6.2 GPa Precompressed Static-Dynamic Compression Technique for Wide-Range Equation-of-State Investigations
MA Xuyang, TU Yuchun, HE Zhiyu, JIA Guo, FANG Zhiheng, WANG Peipei, HUANG Xiuguang
2026, 40(6): 060101.
Effect of Boron Nitride Content on the Explosion Performance of On-Site Mixed Emulsion Explosives
FU Jiakun, LIU Feng, ZHU Zhengde, CHEN Chuanbin
2026, 40(5): 050109.
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:
Dynamic response characteristics of soil-rock strata induced by transient waves
LIU Weiwei, JIA Lei, FU Yanqing, WANG Teng
 doi: 10.11858/gywlxb.20261089
PDF(0)
Abstract:
Clarifying the dynamic response characteristics of soil-rock strata under transient waves is of great significance for accurately evaluating and controlling the dangerous effects of dynamic effects such as blasting or earthquake. Aiming at the typical soil-rock stratum structure, the standard linear solid is introduced to describe the soil constitutive, and the physical analysis model of soil-rock stratum under transient wave action is constructed. Based on the characteristic line method and the displacement continuity method, the explicit mathematical description of transient wave propagation in soil-rock stratum is established. Finally, through parameter analysis, the influence of dimensionless thickness, relaxation time and equivalent dynamic modulus of soil on the propagation characteristics and site amplification effect of transient wave in soil-rock stratum is systematically discussed. The results show that: (1) when considering the viscoelasticity of soil, the amplification effect of surface site is small; (2) With the increase of dimensionless thickness, the particle velocity amplification effect of soil-rock interface decreases first and then stabilizes, the peak value of stress wave increases first and then decreases, and the surface velocity amplification effect increases first and then decreases. (3) With the increase of relaxation time, the amplification effect of interface vibration velocity is basically unchanged, the surface amplification effect increases first and then stabilizes, and the peak value of stress wave increases gradually. (4) With the increase of equivalent modulus, the amplification effect of surface and interface vibration velocity is almost unchanged, and the peak value of stress wave decreases gradually.
Study on dynamic mechanical properties and mesoscopic simulation of steel-polyoxymethylene hybrid fiber reinforced concrete
CAO Yibin, CHEN Ruihao, LANG Zhijun, ZHANG Lidan, ZHANG Wenyu, LIU Lei
 doi: 10.11858/gywlxb.20261096
PDF(0)
Abstract:
To reveal the dynamic mechanical properties and meso-damage characteristics of steel-polyoxymethylene hybrid fiber reinforced concrete (HFRC).In this study, the Split Hopkinson Pressure Bar (SHPB) tests were conducted to investigate the dynamic mechanical properties of plain concrete (PC), steel fiber reinforced concrete (SFRC), and HFRC. Scanning electron microscopy (SEM) was employed to observe the micromorphology and fracture characteristics of the fiber-matrix interface. A three-dimensional mesoscopic numerical model was established using LS-DYNA software, which incorporated polyhedral aggregates, mortar, interfacial transition zone (ITZ), and two types of fibers. The damage evolution, energy dissipation, and stress characteristics of HFRC were revealed from a mesoscopic perspective. The results show that polyoxymethylene (POM) fiber exerts a synergistic reinforcing effect with steel fiber, significantly improving the dynamic impact resistance of the material. Specifically, steel fiber exhibits strong bonding with the matrix and dissipates energy via interfacial slip, whereas POM fiber presents weak interfacial bonding and dissipates energy mainly through pull-out, fracture, and deformation. The two fibers function synergistically at different stages of crack propagation, effectively enhancing the crack resistance and energy dissipation capacity of the matrix. Meso-simulations indicate that mortar serves as the primary energy-dissipating component, and ITZ is the weakest region that fails first under impact loading. With increasing impact air pressure, the energy absorption and peak stress of each component increase significantly. Hybrid fibers can reduce energy reflection and kinetic energy loss, facilitating more impact energy to be absorbed and dissipated. The peak stress of ITZ is the most sensitive to loading rate, while aggregate shows the lowest sensitivity. The established mesoscopic model can well reproduce the dynamic mechanical behavior and damage characteristics of HFRC, providing theoretical and numerical references for the application of HFRC in impact-resistant engineering.
LI Rui, SUN Rui, LIU Keyang, HU Minghang, CHEN Yajing, WANG Quan, LI Xuejiao, XU Xiaomeng, CUI Xiaorong
 doi: 10.11858/gywlxb.20261101
PDF(6)
Abstract:
Energy Absorption and Cushioning Performance of Second-Order Hierarchical Corrugated Sandwich Structures under Static and Dynamic Loads
LIU Tao, LI Zhaokai, LIU Xiaoyong, XIE Jiamiao, HAO Wenqian
 doi: 10.11858/gywlxb.20261105
PDF(6)
Abstract:
Renowned for the nature of high specific stiffness and strength, lightweight, and excellent energy absorption capabilities, corrugated sandwich structures have found broad utility in fields such as aerospace and vehicle protection. This study introduces a hierarchical design concept into such structures by proposing a second-order hierarchical corrugated sandwich (SHCS) structure and exploring its mechanical and impact performance. An analytical expression for the peak load is derived, and the discrepancy between the theoretical calculations and the finite element analysis results is within 10%. Using numerical simulation methods, the deformation modes of structures is investigated with different numbers of minor supports (n). The deformation modes of structures are categorized into three modes: progressive buckling, transition buckling, and global buckling. Quasi-static compressive and dynamic impact responses of the structure are further analyzed under various impact velocities and geometric parameters of minor support structure, with particular focus on deformation mechanisms and energy absorption characteristics. The results indicate that under quasi-static compression, the thickness of minor support structure exerts a more significant influence on energy absorption. Increasing the thickness from 0.8 to 1.6 mm enhances the specific energy absorption by 67.3%. Under low-velocity impact, both the peak load and specific energy absorption follow a parabolic trend with increasing minor support panel thickness. During dynamic loading, the energy absorption performance of the structure improves as the impact velocity increases. Under high-velocity impact conditions, the peak load and specific energy absorption show an increasing trend with greater minor support panel thickness.
Nitrogen Aggregation Induced by High-Pressure High-Temperature Pretreatment and Its Effect on the Apparent Charge-State Response of NV Centers in Diamond
LI Weijian, CHEN Ning, WANG Hongwei, WANG Hao, ZHAO Shixiu, PAN Yilong
 doi: 10.11858/gywlxb.20261106
PDF(4)
Abstract:
This study investigated how the initial C-center nitrogen content, electron irradiation, vacuum annealing, and high pressure high temperature pretreatment affects NV center formation and charge-state regulation in HPHT type-Ib diamond. Two groups of HPHT diamond single crystals with different C-center nitrogen contents were used. The samples were treated by electron irradiation, vacuum annealing at 600 to 1000 ℃, and HPHT pretreatment at 5 GPa and 1100 to 1900 ℃. FTIR, Raman, and PL spectra were used to analyze nitrogen aggregation, lattice state, and NV-related photoluminescence. The results show that the high-nitrogen sample forms stronger NV-related emission, whereas irradiation-induced vacancies in the low-nitrogen sample tend to remain as GR1 centers. Increasing the vacuum annealing temperature enhances the total NV-related emission but reduces the relative response of negatively charged NV centers. HPHT pretreatment does not directly produce many NV centers. Instead, it promotes nitrogen aggregation and modifies the apparent charge-state ratio after subsequent irradiation and annealing. These results indicate that HPHT pretreatment can serve as a preceding processing parameter for tuning the initial defect state before irradiation and regulating the later charge-state response of NV centers.
Effects of impact velocity and pulse duration on spallation behavior and void evolution on a metastable Si4V5Mn5Cr10Co30Fe46 high-entropy alloy
LIU Zhengyuan, ZHANG Tuanwei, ZHAO Jiawei, DU Shiyu, ZHANG Lei, WANG Zhihua
 doi: 10.11858/gywlxb.20261102
PDF(14)
Abstract:
Metastable high-entropy alloys can exhibit coupled phase transformation, strain localization and damage evolution under shock loading, and their spall response is governed not only by loading intensity but also by the duration of the stress pulse. To clarify the distinct roles of impact velocity and pulse duration in spallation behavior and void evolution, plate-impact experiments were conducted on a Si4V5Mn5Cr10Co30Fe46 metastable high-entropy alloy using a single-stage light-gas gun. The free-surface velocity response, spall parameters, local microstructural evolution and three-dimensional void morphology under different loading conditions were systematically investigated. The results show that, at a fixed specimen thickness, increasing the impact velocity from 282 m/s to 553 m/s markedly raises the peak free-surface response and peak compressive stress, whereas the spall strength changes only slightly. Meanwhile, surface voids evolve from dispersed nucleation to localized clustering, accompanied by increased fractions of HCP/BCC phases and a pronounced rise in high-KAM regions near the voids, indicating that higher impact velocity promotes local phase transformation, lattice distortion and concentrated damage development. Micro-X-ray computed tomography further reveals that increasing impact velocity drives the voids towards larger volumes, stronger spatial concentration and more complex morphologies. In contrast, under nearly constant impact velocity, as the specimen thickness increases from 1 mm to 2 mm, the pulse duration is prolonged from 0.075 μs to 0.25 μs, and the spall strength correspondingly increases from 1.40 GPa to 1.83 GPa. The spatial distribution of voids gradually changes from dispersed to centrally concentrated, with an enhanced tendency for interconnection, indicating that a longer pulse duration is more favorable for damage accumulation and localized coalescence. By combining free-surface velocity histories, two-dimensional microstructural characterization and three-dimensional void statistics, it is shown that impact velocity mainly controls the instantaneous driving force for damage during spallation, whereas pulse duration primarily governs the time window for damage accumulation and the extent of localization. Together, these two factors determine the nucleation sites, growth paths and final failure mode of spall damage in this metastable high-entropy alloy.
Effect of ODEA Content on Micro-morphology and Rheological Properties of Mixed Emulsion Matrix
HE Zhiwei, LI Yuanlong, YUE Xing, HUANG Zhenyi, YUE Jiawei, HU Qianhao, GE Luyao
 doi: 10.11858/gywlxb.20261103
PDF(22)
Abstract:
To investigate the effect of ODEA on the microscopic morphology and rheological properties of on-site mixed emulsion matrix, four groups of samples were prepared by compounding Span-80 with different mass fractions of ODEA. The microscopic morphology and rheological properties of the samples were characterized using an optical microscope, a laser particle size analyzer and a rotational rheometer. The results show that with the increase of ODEA mass fraction, the average particle size and the dispersion degree of internal phase droplets of the on-site mixed emulsion matrix decrease first and then increase. When the ODEA content is 1 wt%, the sample has the smallest average particle size, the most uniform particle size distribution and the highest detonation velocity. The addition of ODEA improves the viscosity, storage modulus and cohesion of the composite system samples. Within the experimental temperature range, all four groups of samples can meet the pumping requirements. Compared with the single Span-80 system, the sample with 1 wt% ODEA has the strongest elastic recovery capacity under external force and the best shear stability. This study provides relevant experimental basis for the research on adding small-molecule alkanolamide co-emulsifiers into the formulation of on-site mixed emulsion explosives.
Advances in Impact Discharge Theory and Dynamic Phase Transition Mechanism in Ferroelectric Ceramics
FENG Qiu, GE Mingyue, XIONG Zhengwei, LI Jun, JIN Ke, GAO Zhipeng
 doi: 10.11858/gywlxb.20261083
PDF(18)
Abstract:
Ferroelectric ceramics have become a core medium for pulsed power devices owing to their high remanent polarization and shock-induced depolarization via structural phase transitions. This review systematically summarizes the research progress on the electrical response behaviors, theoretical models, and phase transition mechanisms of perovskite ferroelectric ceramics under shock wave loading from the 1950s to the present, with particular emphasis on the last two decades. The reviewed material systems include the lead-based lead zirconate titanate family and lead-free systems such as bismuth sodium titanate based, potassium sodium niobate based, bismuth ferrite based, and silver niobate based systems. Regarding material evolution, the transition from traditional high-performance lead-based dominance to novel lead-free systems featuring high energy density, high power density, and environmental compatibility is outlined. In terms of theoretical models, the universal evolution of discharge waveforms with shock pressure is analyzed, and the development from the constant-current source model, phase transition kinetic model, to the relaxation model is summarized, with a focused review on the "piezoelectric-ferroelectric" dual-mechanism discharge framework and the construction logic of the shock ferroelectric equation of state covering the full pressure range. Concerning physical mechanisms, the essential differences in phase transition behaviors between dynamic shock and static high-pressure loading are emphasized, and the microscopic mechanisms specifically the uniaxial stress lowering phase transition barriers and the polycrystalline orientation statistics broadening the transition window are elucidated. Finally, current limitations such as the phenomenological nature of piezoelectric models and the inadequacy of multiscale non-equilibrium simulations are identified, and future directions including machine-learning-potential molecular dynamics and texture engineering are discussed.
The Influence of Aluminum Powder Particle Size on the performance of HMX-Based Aluminized Explosives in Air Blast
GUO Liuwei, TAO Shixing, HUANG Longjie, LI Shengfu, MA Qingpeng
 doi: 10.11858/gywlxb.20261107
PDF(17)
Abstract:
To investigate the effect of aluminum powder particle size on the performance of HMX-based aluminized explosives, the shock overpressure and fireball performance in air blast of HMX-based aluminized explosives with aluminum powder particle sizes of 0.15 μm, 10 μm, and 50 μm were studied. The results indicate that both decreasing and increasing the aluminum powder particle size can enhance the peak value of the near-field shock overpressure. However, neither size change effectively influences the shock overpressure in farther field. The mechanisms by which small and large particle size aluminum powders enhance the shock overpressure differ: smaller particles accelerate the energy release process, while larger particles undergo a process of dispersion followed by energy release, enhancing the energy release at the edge of the explosion fireball. Changing the aluminum powder particle size can extend the positive pressure duration of the shock wave but results in a faster decay rate in the high-pressure segment, leading to a shorter pressure half-life and making it difficult to effectively improve the total impulse of the positive pressure zone or the impulse in the high-pressure segment. Additionally, altering the aluminum powder particle size affects the fireball performance. As the particle size increases, the size of the explosion fireball monotonically increases. However, both overly large or small aluminum powder particle sizes are detrimental to the surface temperature of the explosion fireball.
A prediction model for the remaining velocity of simultaneous multi-fragment impacts based on feedforward neural networks
LI Yuzhuo, KANG Yaowen, GAO Yueguang, FU Jianping, REN Kai, YANG Rui
 doi: 10.11858/gywlxb.20261074
PDF(15)
Abstract:
Numerical Simulation and Experimental Study of Deck-Charging Effects in Bench Blasting of Soft-Hard Interbedded Rock Masses
JIA Yingxin, HE Tiezhu, YANG Guoliang, LIU Qiang, HAN Zimo, MA Mingxue
 doi: 10.11858/gywlxb.20261100
PDF(17)
Abstract:
Due to the significant differences in mechanical properties among different rock strata in soft-hard interbedded rock masses, open-pit bench blasting is prone to problems such as excessive fragmentation of soft rock, insufficient breakage of hard rock, and large discreteness in muckpile fragment size, which adversely affect blasting quality and subsequent loading and hauling efficiency. To optimize blasting energy distribution and improve rock fragmentation in soft-hard interbedded rock masses, this study takes the Jilongde open-pit coal mine as the engineering background. A three-borehole fluid-solid coupling numerical model was established using ANSYS/LS-DYNA to comparatively analyze the blasting responses of continuous charging, air-deck charging, and water-deck charging under different deck ratios. Field blasting tests were also conducted to validate the numerical simulation results. The peak effective stress in the soft rock zone, the mean effective stress in the hard rock zone, and the box-counting dimension of the damage section were selected as evaluation indices to systematically reveal the influence of different charging structures on stress transmission, energy distribution, and damage range in soft and hard rock strata. The results show that deck charging can effectively regulate the release process of explosive energy along the borehole direction. As the deck ratio increases, the peak effective stress in the soft rock zone generally decreases, indicating that stress concentration in the soft rock is weakened. The mean effective stress in the hard rock zone first remains relatively stable and then decreases, suggesting that an appropriate deck ratio can maintain the stress level required for hard rock fragmentation, whereas an excessive deck ratio reduces the breakage effect of hard rock. Considering both stress response and damage fractal characteristics, deck ratios of 17.0 % and 18.5 % can effectively balance soft-rock fragmentation control and hard-rock breakage, among which the 17.0 % case exhibits a more favorable box-counting dimension of the damage section and a more reasonable damage distribution. Field test results further indicate that, compared with continuous charging, the 17.0 % deck charging scheme can significantly alleviate excessive pulverization of soft rock and improve muckpile fragment size uniformity. In particular, water-deck charging yields the lowest nonuniformity coefficient and the best overall blasting performance. The findings provide a reference for the optimization of charging structures and refined blasting design in open-pit mines with soft-hard interbedded rock masses.
WANG Nan, LI Xuhang, PAN Ruochen, SONG Guangjun, BIE Bixiong, FAN Duan, CAI Yang
 doi: 10.11858/gywlxb.20261093
PDF(17)
Abstract:
To obtain dynamic performance of materials at cryogenic temperature for development of polar science and engineering technology, a low-temperature Split Hopkinson Pressure Bar (SHPB) system is developed. Compared to traditional low-temperature SHPB system, this system places pressure bars and samples in the cryogenic environment, effectively mitigating temperature variance in the sample and local temperature gradient in the pressure bar due to heat conduction. This design reduces the complexity of data processing. Additionally, the vaporization of liquid nitrogen expels air from the low-temperature chamber, keeping the pressure bars and samples dry and preventing common freezing issues encountered in low-temperature SHPB experiments. As a result, this significantly enhances the accuracy and consistency of the experimental results. Moreover, compared to auto-assembling low-temperature SHPB devices, the new system is simpler and more convenient. The relevant parameters of the pressure bar and strain gauge under cryogenic conditions were calibrated. Finite element simulation indicate that the temperature gradient along the bars introduces a measurement deviation of up to 10% in stress in traditional SHPB system; in contrast, the specialized cryogenic SHPB system effectively reduces measurement errors and substantially enhances measurement accuracy. The dynamic performance of 2024Al alloy under cryogenic conditions was tested to validate the operational reliability of the developed SHPB system. Experimental results demonstrate that, compared with traditional SHPB devices, the novel system exhibits significantly higher repeatability and reliability in cryogenic experiments.
Numerical Simulation on Anti-Penetration Performance of SiC/Ti-Alloy Interpenetrating TPMS Structures
ZHOU Haotian, LI Yinan, TAN Bin, MENG Yuquan, SONG Weidong
 doi: 10.11858/gywlxb.20261068
PDF(36)
Abstract:
Ceramic/metal composite materials have been widely used in national defense, military industry, and aerospace fields as lightweight impact-resistant structures with high specific strength and high energy absorption efficiency. With the development of 3D printing technology, it has become possible to fabricate complex lattice structures based on Triply Periodic Minimal Surfaces (TPMS). In this paper, an interpenetrating TPMS ballistic composite structure composed of silicon carbide (SiC) ceramic and titanium alloy (TC4) is designed. A series of numerical simulations are carried out under single-projectile and double-projectile penetration conditions using ABAQUS. The damage modes, penetration depth, and ballistic limit velocity of the proposed structure and pure SiC target are compared and analyzed.The numerical results show that different interpenetrating TPMS structures exhibit distinct damage and failure modes. The three-dimensional topological configuration restrains crack propagation inside the ceramic, resulting in slighter overall damage than the pure SiC target. The damage caused by the second projectile further develops along the penetration region of the first projectile, accompanied by an increase in penetration depth. Compared with the pure SiC target, the three interpenetrating TPMS targets present smaller penetration depth and higher ballistic limit velocity. When the projectile can perforate the target, the P-type structure shows better ballistic performance against low-velocity projectiles, while the D-type structure is superior against high-velocity projectiles. It is demonstrated that, at the same areal density, interpenetrating TPMS targets possess better ballistic performance than pure SiC.This study can provide technical support and theoretical basis for the design of novel lightweight ceramic armor.
Effects of Joint Angle on Crack Propagation Behavior under Dynamic Impact and Quasi-static Loading
WANG Qisheng, LI Runpeng, WANG Xin, QIU Peng, LEI Jianyin, LIU Zhifang
 doi: 10.11858/gywlxb.20261097
PDF(33)
Abstract:
To reveal the influence of joint angle on the crack propagation behavior of PMMA under different loading conditions, a digital laser caustics system was employed in combination with Hopkinson pressure bar dynamic impact tests and quasi-static three-point bending tests. Experiments were carried out on PMMA three-point bending specimens containing joint defects with angles from 15° to 75°. The results show that the crack propagation process of jointed specimens can be divided into three phases: pre-crack initiation, crack-joint interface interaction, and crack re-initiation. With increasing joint angle, the distance between the end point of phase 1 and the re-initiation point of phase 3 gradually increases, indicating that the joint angle changes the interaction range near the crack-joint interface. Compared with static loading, the tensile-shear coupling at the crack tip and the stress-wave disturbance are more pronounced under dynamic impact loading. The peak crack propagation velocity of the 30° jointed specimen reaches 638.4 m/s during the re-initiation phase, which is higher than that of the 75° jointed specimen. In this phase, the peak mode-I stress intensity factor decreases from 2.31 MPa·m1/2 for the 15° specimen to 2.03 MPa·m1/2 for the 75° specimen. The peak mode-II stress intensity factor mainly appears in low- to medium-angle jointed specimens, indicating that a strong local shear disturbance exists near the crack tip at the instant of interface re-initiation. Compared with dynamic loading, the external load input and crack-tip stress accumulation under quasi-static loading are more gradual. Although the crack still exhibits dynamic propagation characteristics after instability, the shear component is weaker, and crack propagation tends to be dominated by mode-I opening. The crack-path complexity analysis based on the box-counting fractal dimension shows that the fractal dimensions of crack paths under dynamic impact loading are generally higher than those under quasi-static loading, indicating that dynamic loading is more likely to induce complex path evolution such as crack deflection, retention, and re-initiation. Among all specimens, the 30° jointed specimen has the largest fractal dimension and the highest crack-path complexity. The results indicate that the joint angle affects crack propagation behavior by changing the crack-interface interaction range, the tensile-shear stress ratio at the crack tip, and the crack-path complexity, providing a reference for fracture response analysis of defective engineering structures.
CFD-FEM Simulation on the Effect of Shell Perforation after Ablation on Combustion Heat Transfer and Explosive Ignition
ZHANG Jinxuan, ZHENG Songlin, YU Yin
 doi: 10.11858/gywlxb.20261078
PDF(34)
Abstract:
Fast cook-off assessment has always been a crucial topic in ammunition safety research. However, previous simulation studies typically did not consider the influence of factors such as outer casing ablation on ammunition ignition behavior. This paper pre-sets different casing ablation hole sizes and gap widths, employing the computational fluid dynamics methodology to obtain the flame temperature field under ablation conditions. The flame thermal load is then transferred as a boundary condition for finite element method calculations, thereby investigating the ignition behavior of the explosive after casing ablation. Simulation results indicate that when the casing gap is 5 mm, the ignition position of the explosive varies with the size of the ablation hole, consistently located at the edge of the hole. When the casing gap reaches 10 mm or more, the ignition position does not change with the ablation hole size and is consistently located directly above the explosive. On the other hand, the size of the casing ablation hole and the casing gap have no significant effect on the ignition time of the explosive. Mechanism analysis suggests that the ignition position is influenced by the reaction between the residual fuel in the high-temperature gas and the air within the casing gap. This study can provide a theoretical reference for accident emergency response, ammunition safety assessment, and improvement.
Numerical Simulation Study on Hypervelocity Impact of High-Entropy Alloy Protective Structure
YIN Yunfei, YANG Qiuzu, GUO Jia’ao, LI Zhiqiang
 doi: 10.11858/gywlxb.20251275
HTML PDF(120)
Abstract:

The space debris problem has become one of the most pressing challenges in the field of space environment protection. Currently, most spacecraft shielding systems adopt Whipple-type structures, in which aluminum alloys are commonly used as bumper materials. In this study, the smoothed particle hydrodynamics (SPH) method implemented in the AUTODYN software was employed to numerically investigate the hypervelocity impact of spherical projectiles on high-entropy alloy (HEA) protective structures. The characteristics of the resulting debris clouds, including fragment number, mass distribution, and momentum, were systematically analyzed under various impact conditions. In addition, the effects of impact velocity and the ratio of bumper thickness to projectile diameter (t/D) on the hypervelocity impact response of HEA shielding structures were examined. The results show that, under identical impact conditions, the debris cloud characteristics generated by HEA protective structures differ significantly from those of aluminum alloy structures, namely: the total number of fragments increases by approximately 51.86%; the number of small-mass fragments increases by about 79.56%, while the number of large-mass fragments decreases; and the maximum debris cloud momentum along the impact direction (z-direction) is less than 75% of that associated with aluminum alloy structures across multiple projectile diameters. Parametric analyses indicate that the expansion of the debris cloud is primarily governed by impact velocity, with higher velocities leading to faster expansion, while the influence of the t/D is relatively minor. In contrast, the mass and kinetic energy of hazardous fragments (large-mass/high-energy fragments) are mainly affected by t/D, with higher values resulting in greater impact severity. These findings provide theoretical support and reference for the application of high-entropy alloys in next-generation spacecraft shielding structures.

Experimental and numerical simulation calculation of spherical target explosion concrete throwing
DU Yong, ZHANG Yuping, DAI Xianghui, WANG Kehui, WANG Kaiqiang, ZHAO Shengwei
 doi: 10.11858/gywlxb.20261087
PDF(37)
Abstract:
In order to study the effect of blasting parameters on the throwing of crushed stone in engineering blasting, an experimental system for the throwing of spherical target explosive concrete based on digital high-speed image acquisition system was designed, and the experiments on the throwing of concrete under the condition of single-hole and three-hole charge were carried out, and the failure modes and throwing quality of spherical target after blasting under the condition of single-hole and three-hole charge were compared. The simulation model of spherical target explosion is established, and the accuracy of the model is verified by experiments. The concrete throwing characteristics of different explosive holes and different spherical target strength are simulated and calculated. The results show that when the quantity of explosive is 5g, the mass of the three holes is 89% more than that of the single hole explosive; when the quantity of explosive is 10g, the mass of the three holes is 82.1% more than that of the single hole explosive; when the quantity of explosive is constant, increasing the quantity of the holes can effectively improve the quality of the throwing concrete and the area of the blasting pit, which is helpful to improve the production efficiency and reduce the cost. The simulation results are in good agreement with the experimental data.
Study On The Effect Of Deep-Water Hydrostatic Pressure on the Performance of Ammonium-Amine Explosive
WANG Junhao, HUANG Wenyao, WANG Quan, TONG Kai, YUAN Wenjie, HUANG Daguang, WENG Qiuhong
 doi: 10.11858/gywlxb.20261082
PDF(33)
Abstract:
To investigate the explosion performance of ammonium-amine explosives under deep-water hydrostatic pressure conditions, an experimental simulation device for blasting charges in a deep-water hydrostatic pressure environment was designed. An optical microscope was used to observe the microscopic bubble changes in the explosive after exposure to atmospheric pressure and hydrostatic pressures of 0.1 MPa, 0.2 MPa, 0.3 MPa, and 0.4 MPa for 1 hour followed by pressure relief and recovery. Fiji image analysis technique was employed for microscopic characterization of the size distribution of sensitizing bubbles. The density of the explosive under hydrostatic pressure was measured. The detonation velocity of the explosive in a PVC charge tube with an outer diameter of 40 mm was tested under atmospheric pressure, after pressure relief, and under hydrostatic pressure. The results indicate that the deep-water hydrostatic pressure environment has a significant effect on the performance of ammonium-amine explosives. For the explosive after pressure relief and recovery, when the pressure ranges from 0 to 0.2 MPa, the number of microscopic bubbles increases with increasing pressure, leading to an increase in effective "hot spots", and the detonation velocity increases from 4313 m·s⁻¹ to 4621 m·s⁻¹. When the pressure exceeds 0.2 MPa, bubbles coalesce and merge, reducing the number of effective "hot spots", and at 0.4 MPa, the detonation velocity decreases to 4072 m·s⁻¹. For the explosive detonated under hydrostatic pressure, as the pressure increases from 0 to 0.4 MPa, the explosive density increases from 1.02 g·cm⁻³ to 1.34 g·cm⁻³, following the relationship ρ = 0.47P0.43 + 1.02. The detonation velocity increases from 4313 m·s⁻¹ under atmospheric pressure to 4448 m·s⁻¹ at 0.1 MPa, decreases to 3412 m·s⁻¹ at 0.3 MPa, and results in failure to detonate at 0.4 MPa. The ammonium-amine explosive exhibits a certain resistance to deep-water hydrostatic pressure and demonstrates good recovery after pressure relief.
Property Simulations of Ultrahigh-Pressure Hydrous Magnesium Silicate Phases at High Temperatures and Pressures: Exploring the Possibility of Deep Water in Terrestrial Super-Earths
REN Hang, ZHANG Wenqi, LIU Lei
 doi: 10.11858/gywlxb.20261061
PDF(37)
Abstract:
Hydrous magnesium silicate is considered as a potential water-rich reservoir in the early Earth's interior. Investigating its behavior under extreme high-temperature and high-pressure conditions is crucial for understanding the internal structure models and water presence potential of super-Earths. By using first-principles molecular dynamics simulations, this study systematically investigates the stability and elastic properties of β-Mg₂SiO₅H₂ within the pressure range of 500-900 GPa and temperature range of 2000-6000 K. The results indicate that the system remains thermodynamically stable across the entire studied pressure-temperature range, with no structural phase transitions observed. Calculations of the mean square displacement reveal the transition interval for the superionic state: at 500 GPa and 2000 K, 700 GPa and 3000 K, and 900 GPa and 3100 K, all atoms remain confined within the lattice, with the system in a normal state. When the temperature increases to 500 GPa and 4000 K, 700 GPa and 5000 K, and 900 GPa and 6000 K, the mean square displacement of H atoms exhibits a linear increase, while the framework atoms (Mg, Si, O) remain localized. The proton trajectories form a diffuse network, exhibiting the characteristics of a superionic state with a "solid framework + liquid-like ions". Simulation results show that the density of β-Mg₂SiO₅H₂ increases linearly with pressure. The shear modulus decreases nearly linearly with increasing temperature, while the bulk modulus shows a significant reduction at 900 GPa and 6000 K, likely in response to the transition to the superionic state. The variations in shear wave velocity and compressional wave velocity are primarily controlled by pressure and temperature, increasing with pressure and decreasing with temperature. Under high pressure (900 GPa), the superionic state leads to a notable decrease in the compressional wave velocity of β-Mg₂SiO₅H₂, suggesting that the superionic transition induces structural softening. This study confirms that β-Mg₂SiO₅H₂ can remain stable under the deep mantle pressure-temperature conditions of super-Earths with masses 5-8 times that of Earth and can transition to a superionic state under specific conditions. Its high water content of up to 11.4 wt% and efficient proton transport capability have significant implications for deep water cycles and the habitability of terrestrial planets, providing key theoretical insights for understanding planetary interior dynamics.
A Method for Simultaneous Measurement of Shock Temperature and Sound Velocity in Transparent Minerals
REN Lei, GAN Bo, HUANG Yuqian, HE Qing, ZHANG Youjun
 doi: 10.11858/gywlxb.20261079
PDF(40)
Abstract:
Method/t/nDynamic compression techniques serve as a critical approach to generate ultrahigh pressure and high temperature conditions, with wide applications in high-energy-density physics, geophysics, and defense-related research. Shock temperature is a key parameter for characterizing thermodynamic states and constructing equations of state, whereas sound velocity is highly sensitive to phase transitions and critically constrains elastic properties. Nevertheless, conventional shock experiments typically require separate measurements of sound velocity and temperature, which increases experimental complexity and may introduce state mismatch between different physical quantities. Here, we develop a method for the simultaneous measurement of shock temperature and longitudinal wave velocity in transparent minerals. This method determines shock temperature from thermal radiation emitted at the shock front that transmits through the uncompressed region, while the longitudinal sound velocity (<italic>V</italic>P) is derived from the temporal evolution of radiation signals during compression and release. Natural single-crystal calcite was selected to validate the method through shock experiments. At a shock pressure of 115.7 GPa, the shock temperature and <italic>V</italic>P were determined to be 3810 K and 9.41 km/s, respectively. The results indicate that calcite may undergo shock melting under these conditions, validating the feasibility of our method. Our study provides a novel approach for simultaneously measuring thermodynamic and elastic properties of transparent minerals under extreme conditions, providing important constraints for understanding deep-Earth processes and evaluating the dynamic response of transparent window materials under high-energy-density environments.
Fullerenes under High Temperature and High Pressure
SONG Jing, WANG Lin
 doi: 10.11858/gywlxb.20261069
PDF(36)
Abstract:
Effect of freezing temperature on mechanical properties of layered hail**
SHI Xiaopeng, FANG Jianglu, JIE Jiang
 doi: 10.11858/gywlxb.20261054
PDF(44)
Abstract:
This thesis aims to clarify the influence of layered structure and different freezing temperatures on the mechanical properties of hails. The transparent hails with a diameter of 50 mm and layered hails with distinct structures were prepared using a self-made mold. A series of quasi-static compression experiments were conducted for the simulated hails (both transparent hail and layered hail) through a universal testing machine at different freezing temperatures (-10 °C, -20 °C, -30 °C, -40 °C), and combined with microscopic structure analysis to determine the failure mechanism. Based on the experimental results, it is found that as the freezing temperature decreases, the compressive strength of both transparent hails and layered hails significantly increases. Layered hails are more sensitive to temperature changes. At -40 °C, the average compressive strength of layered hails is approximately 6 times higher than that at -10 °C, while that of transparent hails is approximately 4 times higher. When layered hails fail, they exhibit multiple cracks, while transparent hails have single main crack. The layered structure is rich in bubbles, which can change the crack direction and thereby enhance its compressive strength.
Pressure-Tuned Bandgap and Optoelectronic Properties of Molybdenum Diselenide
ZHANG Shenghan, QI Wenming, ZHU Zhikai, DONG Hongliang, DONG Lan, CHEN Zhiqiang
 doi: 10.11858/gywlxb.20261062
PDF(38)
Abstract:
MoSe2 is a promising candidate for tunable near-infrared photodetectors owing to its band gap alignment with optical fiber communication windows. We combine diamond anvil cell techniques with density functional theory calculations to investigate pressure-induced structural evolution and optoelectronic modulation in MoSe2. X-ray diffraction and Raman spectroscopy confirm that the hexagonal 2H phase remains stable up to 10 GPa, exhibiting pronounced anisotropic compression: the c-axis compressibility is approximately three times that of the a-axis. Infrared reflectivity spectra show a monotonic increase in reflectance with pressure, indicating a bandgap narrowing trend. First-principles calculations reveal a pressure-driven downward shift of the conduction band minimum, with the band gap narrowing linearly from 1.24 eV at ambient pressure to 0.77 eV at 4 GPa and optical absorption being enhanced. Photocurrent increases progressively from 0.4 to 3.9 GPa, peaking at approximately twice the ambient value before vanishing at 4.3 GPa due to overwhelming dark current. The consistency between theoretical predictions and experimental observations elucidates the electronic origin of the pressure-tuned optoelectronic response, providing a foundation for MoSe2-based pressure-modulated spectral devices.
Performance degradation and failure characteristics of lithium-ion batteries under impact loading
ZHOU Xuehui, HUANG Zixuan, ZHANG Xinchun, RAO Lixiang, YANG Shuai
 doi: 10.11858/gywlxb.20261077
PDF(37)
Abstract:
To address performance degradation and safety challenges in lithium-ion batteries under impact, this study aims to investigate the damage and failure characteristics of batteries under various collision forms. For NCR18650BD cylindrical lithium-ion batteries, electrochemical impedance spectroscopy (EIS), cyclic voltammetry (CV) methods, and the multi-physics in-situ monitoring platform were employed, the capacity decay, cycling ageing, and impedance degradation patterns of batteries under ball, cylindrical, and plane punch were analyzed, and the influence of collision form on failure and thermal behavior was revealed. Results show that at an impact energy of 5 J, only the plane punch causes capacity degradation due to lithium inventory loss in the battery at C/20 rate. At 3C rate, the smaller punch means the faster cycle ageing. Furthermore, the larger the punch size, the increased electrochemical impedance of the battery. At an impact energy of 20 J, the smaller punch results in the lower force peak, the more severe battery failure, and the higher probability of thermal runaway. As the impact velocity increases, the overall stiffness of the battery increases, leading to a worsening of force-electric failure. This research provides the theoretical foundation and technical support for the optimized design and safety assessment of lithium-ion batteries.
Study on the Effect of Vaseline Modification on the Microstructure and Detonation Performance of Beeswax-Based Emulsion Explosives
LI Ming, WU Hongbo, XIN Youli, WANG Xinqi, HU Pengfei, ZHAO Changxin, ZHANG Wei, ZHANG Chenxi
 doi: 10.11858/gywlxb.20261036
PDF(36)
Abstract:
To address engineering challenges associated with traditional beeswax-based emulsion explosives—such as high oil phase viscosity, difficulty in emulsification dispersion, and susceptibility to hardening and embrittlement at low temperatures—this study introduces Vaseline as a modifier. The effects of variations in oil phase molecular structure on the rheology, microstructure, thermal stability, storage stability, and detonation performance of the emulsion explosives were systematically investigated.The results indicate that the abundant branched alkanes in Vaseline reduce the viscosity of the oil phase system and effectively improve shearing efficiency during the emulsification process. As the mass fraction of Vaseline increases from 0% to 2.4%, the Sauter mean diameter D[3,2] of the emulsion matrix decreases from 8.38 μm to 5.85 μm. Although the onset decomposition temperature T0 and the temperature of maximum weight loss rate Tp decrease slightly, the thermal stability remains within industrial safety requirements. At a Vaseline mass fraction of 1.8%, the system achieves optimal microstructural uniformity with the Polydispersity Index (PDI) dropping to 1.73. Under this condition, the resistance to high-low temperature cycles increases from 10 to 18 cycles, and the detonation velocity reaches a peak of 5180.3 m/s. This study confirms the feasibility of enhancing the comprehensive performance of emulsion explosives by tuning the oil phase microstructure, providing an economical and effective formulation optimization strategy for improving the adaptability of industrial explosives to harsh, temperature-variable environments.
Numerical Study on the Blast Resistance of Cross-shaped Steel-reinforced Concrete Columns Subjected to Near-Field Explosion
REN Shuangchao, WANG Xueji, GUO Xuekang, LI Yi
 doi: 10.11858/gywlxb.20261035
PDF(40)
Abstract:
Steel-reinforced concrete (SRC) columns are widely adopted in critical and high-rise buildings due to their high load-bearing capacity. The failure of SRC columns may trigger progressive collapse of the entire structure when subjected to blast loading. However, research on the blast resistance of SRC columns under near-field explosions remains limited. To fill this gap, a refined finite element model of a cross-shaped SRC column in a real high-rise building was established in LS-DYNA. The effects of scaled distance and axial load ratio on the blast resistance and damage level of SRC columns were investigated under near-field explosion scenarios. The numerical results indicated that, at an axial compression ratio of 0.2, the proportion of shock load distributed within 1300 mm of the mid-height of the SRC column decreased from 61.3% to approximately 42.5% with increasing scaled distance. If interfacial debonding occurred between the confined core region and the unconfined concrete, the SRC column would exhibit flexural-shear failure. Conversely, if no such debonding occurred, flexural failure would be observed in the SRC column, and the damage level could not exceed a medium level. A higher axial load ratio was detrimental to blast resistance when the scaled distance <italic>Z</italic> ≤ 0.6 m/kg 1/3, whereas it became beneficial when <italic>Z</italic> ≥ 0.7 m/kg 1/3. Finally, empirical predictive models were established to estimate the average cumulative impulse and damage index of SRC columns, thereby providing a quantitative basis for rapid post-blast damage assessment.
Exploratory Study on the High-Temperature and High-Pressure Modification of Natural Jadeite Jade
CHEN Zhengjie, ZHOU Li, ZHENG Wei, ZHANG Jiawei, HE Duanwei
 doi: 10.11858/gywlxb.20261070
PDF(46)
Abstract:
This study used low-quality natural jadeite jade, whose principal mineral phase is jadeite, as the starting material, and carried out exploratory high-temperature and high-pressure modification experiments with a cubic-anvil press. Two processing routes were investigated: direct high-temperature/high-pressure treatment within the jadeite stability field, and a route involving melting-quenching in the low-pressure metastable region to prepare a glass precursor followed by crystallization in the high-pressure stability field. The effects of encapsulation material, pressure, temperature, and holding time on the formation and crystallization behavior of the NaAlSi2O6 glass precursor were systematically examined. The results show that the former route did not significantly improve the transparency or color uniformity of the samples, nor did it produce obvious densification of the microstructure. In contrast, in the latter route, high-quality transparent and bubble-free glass precursors were obtained without encapsulation and with Mo or Re encapsulation, among which the Re-encapsulated sample showed the best overall performance. For the Ta-encapsulated sample, the main phase after crystallization was NaTaO3. Crystallization experiments using the Re-encapsulated glass precursor further indicate that higher temperature is more favorable for overall crystallization, while prolonged holding time promotes crystallization from the rim toward the interior and gradually leads to the formation of a fibrous interwoven microstructure. These results demonstrate that, compared with direct high-temperature/high-pressure treatment within the jadeite stability field, the route involving structural reconstruction through a glass precursor followed by high-pressure crystallization exhibits greater potential for modification.
A Constitutive Model for Concrete under Low-Temperature Conditions and Its Application in Penetration Simulations
YU Baoxiang, NING Jianguo, XU Xiangzhao
 doi: 10.11858/gywlxb.20261058
PDF(47)
Abstract:
In polar resource development and high-latitude cold-region engineering, concrete structures are frequently subjected to the coupled effects of low temperatures and dynamic loading. To accurately investigate the mechanical response under these conditions, Split Hopkinson Pressure Bar (SHPB) tests were conducted within a temperature range of 20 °C to -20 °C to characterize the temperature- and strain-rate-dependent properties of concrete. Based on the experimental findings, an improved dynamic damage constitutive model was developed by incorporating a temperature coefficient into the damage evolution equation. This model was further integrated with the Equation of State (EOS) and the RHT yield criterion to fully account for low-temperature effects. The proposed model was numerically implemented via a Fortran-based Vectorized User Material (VUMAT) subroutine. The accuracy of this user-defined model was rigorously verified by comparing numerical stress waveforms and fracture morphologies with SHPB experimental results. Finally, the validated model was applied to simulate projectile penetration into low-temperature concrete targets. The results demonstrate that low temperatures significantly inhibit penetration damage and reduce penetration depth. This phenomenon is attributed to the combined filling and bonding effects of pore ice, which enhance the target’s resistance to impact and large deformations. These findings provide a solid theoretical basis and numerical framework for the impact-resistant design and safety assessment of structures in cold regions.
Study on Thermal Conductivity of AlN/diamond Composites Sintered under High-Pressure and High-Temperature
GONG Fa, LIANG Wenjia, WANG Qiming, LI Qian, LIU Hongwen, HE Peihong, HE Duanwei, PENG Fang
 doi: 10.11858/gywlxb.20261071
PDF(49)
Abstract:
Aluminum nitride ceramics are important heat-dissipation materials for high-power electronic devices. However, the high sintering temperatures required by conventional processing routes limit their practical application and increase fabrication costs. It is therefore necessary to develop preparation methods capable of achieving densification at relatively low temperatures. To address the difficulty of simultaneously obtaining high densification and high thermal conductivity in polycrystalline AlN ceramics under reduced-temperature sintering conditions, this work adopts a stepwise research strategy. First, the densification behavior and thermal conductivity of pure AlN under high-pressure assistance were investigated, and the optimal sintering conditions were identified. Under additive-free conditions, dense pure AlN ceramics with clean grain boundaries were fabricated at 5.0 GPa and 1400 °C, exhibiting a thermal conductivity of 101.6 W·m -1·K-1. Based on these optimized conditions, the AlN/diamond composite system was further studied, and the effect of diamond volume fraction on the structure and properties of the composites was systematically examined. The results show that the thermal conductivity of the composites first decreases and then increases with increasing diamond content, reaching 112.4 W·m-1·K-1 at 33.3 vol%. Mechanistic analysis indicates that interfacial thermal resistance dominates at low diamond contents, whereas at high diamond contents the enhancement of heat transport by thermally conductive diamond pathways becomes more significant. By taking full advantage of the processing benefits of high-temperature and high-pressure technology, this work achieves substantial improvement of AlN-based materials at temperatures lower than those required in conventional sintering, thereby providing a new route for the low-temperature fabrication of high-performance thermally conductive ceramics.
Effect of Nanocrystalline Grain Size on the Dynamic Structure and Damage of Iron
YU Jinmin, GUO Xiuxia, HE Zhiyu, SHAO Jianli
 doi: 10.11858/gywlxb.20251288
HTML PDF(128)
Abstract:

Grain size effect is one of the key factors governing the dynamic mechanical response of metallic materials. Phase transformation iron is selected as the model material, and a series of nanocrystalline polycrystals with identical topology and grain orientation distributions but different grain sizes are constructed to investigate size effects under a fixed grain configuration. Molecular dynamics simulations show that, under high strain rate uniaxial compression, all models undergo elastic deformation, αε phase transition, and high-pressure phase plastic deformation. During the elastic stage, grain boundaries act as a soft layer, leading to lower stresses in the fine grain models than in the coarse grain model. After the structural phase transition, grain boundaries hinder the plastic development of the new phase, so that the fine grain models exhibit higher stress than the coarse grain models. At the onset of phase transition, the threshold of phase transition of smaller grains is lower, and the transformed phase in fine grains mainly forms stacking fault structures, whereas twinning structures appear in relatively larger grains. With increasing strain, the disappearance of twinning and the reconstruction of stacking faults are observed in large grains. Under high strain rate tension, shear strain of grain boundary in the large grain models is highly localized, readily forming continuous shear bands that serve as preferred paths of crack propagation. After grain refinement, shear strain of grain boundary gradually evolves into a diffuse mode, and the effective paths of crack propagation are constrained by the network of grain boundaries. The change of grain boundary effects leads to a non-monotonic variation of fracture strength with the grain size.

Numerical Simulation of CoCrFeMnNi High Entropy Alloy Shaped Charge Jet and Its Penetration into a Target Plate
MENG Yuquan, LEI Rong, LIU Shanshan, WU Xiaobao, SONG Weidong
 doi: 10.11858/gywlxb.20251264
HTML PDF(126)
Abstract:

The continuous advancement of modern armor protection technology has imposed severe challenges on the damage power of shaped charge warheads. Traditional liner materials have become a major constraint in improving penetration depth due to their limited comprehensive performance. High entropy alloys (HEAs), owing to their unique multi-principal element design, exhibit core potentials, such as high strength, high hardness, and excellent dynamic fracture toughness, thereby making them highly promising candidate materials for new-generation liners. In this study, the CoCrFeMnNi high entropy alloy was prepared via selective laser melting (SLM) technology, and the mechanical properties of the alloy were tested under quasi-static and dynamic loadings. The Johnson-Cook (J-C) dynamic constitutive model and related parameters for the CoCrFeMnNi HEA were determined. Using LS-DYNA, shaped charge jet formation models for both copper and high entropy alloy were established, and numerical simulations were carried out on the processes of jet formation and target penetration for copper and HEA. The research results indicate that, when compared with copper, the CoCrFeMnNi HEA liner can form a more stable and continuous jet. Its unique formation and stretching-breakup mechanisms ultimately result in greater penetration depth, confirming the significant advantage of high entropy alloys in the enhanced damage effect.

Dynamic Tensile Properties of CuCrZr Alloy under Electro-Magnetic-Thermo-Mechanical Multifield Coupled Loading
SU Rina, ZHOU Zhongyu, CHEN Xuemiao, LUO Binqiang, WANG Guiji, TAN Fuli, ZHAO Jianheng
 doi: 10.11858/gywlxb.20261052
PDF(47)
Abstract:
CuCrZr alloy is one of the candidate materials for electromagnetic railgun rails. Obtaining its mechanical response under coupled electromagnetic-thermal-mechanical loading is of great significance for the engineering application of CuCrZr alloy. This paper proposes an electromagnetic ring expansion experimental technique assisted by external magnetic field, which stably achieves high strain rate loading above 10⁴/s without significantly increasing the induced current and Joule heating temperature rise in the metal sample ring. Based on this technique, the dynamic tensile properties of CuCrZr alloy under coupled electromagnetic-thermal-mechanical loading were investigated. The stress-strain curves and fracture strain of CuCrZr alloy under high current density, high strain rate, high temperature rise rate, and strong magnetic field environment were obtained. The relevant results provide important references for the application of CuCrZr alloy under multi-physics field coupling conditions.
Experimental Study on Dynamic Response and Damage Characteristics of Early-Age Steel Fiber Reinforced Shotcrete with High-Temperature Variable Curing
CENG Xiangge, XIE Quanmin, XU Yongkang, ZHOU Hui, ZHENG Zhibin, PAN Chong
 doi: 10.11858/gywlxb.20261020
PDF(49)
Abstract:
In the drill-and-blast construction of tunnels under high geotemperature conditions, the mechanism governing damage evolution in early-age Steel Fiber Reinforced Shotcrete (SFRS) subjected to combined thermal and mechanical loads has not been fully understood. A damage constitutive model was developed utilizing the Split Hopkinson Pressure Bar (SHPB) technique and the fractal dimension analysis of CT images. This model facilitated the study of the dynamic response and damage features of SFRS specimens (aged 1 to 3 days) subjected to impact loading, following a curing regime involving high and fluctuating temperatures. The results indicate that the high-temperature varying-temperature curing environment markedly degrades the dynamic mechanical properties of SFRS. The reductions in dynamic strength were measured at 16.1%, 38.1%, and 56.5% for specimens aged 1, 2, and 3 days, respectively, demonstrating a cumulative temperature-induced damage effect with increasing age. The 3-day-old SFRS demonstrated the highest energy absorption efficiency, with its proportion of dissipated energy increasing by 98.9% and 16.7% relative to the 1-day and 2-day-old specimens, respectively. The developed damage constitutive model yielded a goodness-of-fit exceeding 0.9, proving capable of representing the progression of the stress-strain curve for SFRS under impact loading across its elastic, yielding, and failure phases. Furthermore, the research reveals the dynamic damage evolution law of early-age SFRS under thermo-mechanical coupling, which can provide a theoretical basis for the design of support structures in high ground-temperature tunnels.
Dynamic Response of Circular Bridge Piers Subjected to Eccentric Vehicle Impact
LIU Yunting, ZHANG Jie, LIANG Shaomin, YAN Tao, WANG Zhiyong
 doi: 10.11858/gywlxb.20261024
PDF(55)
Abstract:
Existing studies have extensively investigated the global dynamic response of bridge piers under vehicle impact; however, most of them focus on frontal collisions, and the pier response under eccentric collisions remains insufficiently clarified. This study conducts a numerical investigation on the dynamic response of a circular bridge pier subjected to eccentric vehicle impact. A validated finite element model for vehicle-circular pier collision is established. Comparative analyses are performed between different eccentricity levels and the frontal-impact case in terms of failure-mode evolution, impact-force time histories, and key-section displacement responses. In addition, the effects of impact velocity and vehicle mass on the dynamic response under eccentric collision are examined. The results indicate that when the eccentricity is less than 60%, the damage in the impacted region and the bending cracks on the rear face are generally less severe than those under a central frontal collision. As the eccentricity increases such that the engine no longer directly participates in contact, the collision mechanism shifts, and the overall damage level becomes higher than that of the central frontal collision. The velocity effect is non-linear; when the speed increases to 120 km/h, damage is markedly aggravated, accompanied by a significant increase in peak structural displacement. Mass parameters exhibit stage-dependent dominance: the engine mass mainly governs the flexural-shear response during the engine-contact phase, whereas increasing cargo mass substantially amplifies the impact effect in the cargo-contact phase, promoting the pier damage from localized cracking to more extensive flexure-dominated failure.
JH-2 Constitutive Model Parameters and Blast Damage Criterion for Tempered Glass
DU Haoyuan, HAN Lei, REN Yunyan, JIANG Bonan
 doi: 10.11858/gywlxb.20261021
PDF(65)
Abstract:
Tempered glass is widely used in building structures, and its failure behavior and damage threshold under blast loading provide essential support for determining safety stand-off distances and explosion damage ranges. Investigating its failure mechanisms and damage criteria under blast shock waves therefore has important engineering significance. In this study, the mechanical behavior and failure mechanisms of tempered glass over a wide range of strain rates were investigated using a universal testing machine, a split Hopkinson pressure bar, and field emission scanning electron microscopy. Based on the mechanical test results, the material parameters of the Johnson Holmquist Ceramics (JH-2) constitutive model for tempered glass were determined through systematic analysis. Controlled blast loading tests were then conducted using a shock tube, and a finite element model was established to simulate the damage process. Comparisons between test observations and numerical results confirmed the validity of the calibrated constitutive parameters. Finally, a pressure-impulse damage curve was obtained to evaluate the damage threshold of tempered glass. The results provide a theoretical model and an analytical approach for the safety design and damage assessment of tempered glass structures.
Theoretical Modeling and Coupling Effects of Multiple Projectile Penetration-Explosion in Concrete Targets
ZHU Junlong, LI Yuan, WANG Zihao, ZHENG Zhijun
 doi: 10.11858/gywlxb.20261012
PDF(69)
Abstract:

Aiming at the complex multi-physics coupling process involving penetration and explosion in sequential strikes by multiple projectiles, there is currently a lack of corresponding theoretical prediction methods. To address this, this paper proposes a theoretical model based on a decoupling-modeling approach for predicting the entire process of multiple projectiles sequentially penetrating and exploding in concrete targets. The model separately describes the penetration and explosion stages: the penetration process employs dynamic cavity‑expansion theory to construct a resistance function, takes into account the influence of projectile inclination and prior damage, and solves the projectile trajectory by explicit finite‑difference integration of the equations of motion; crater morphology is estimated using empirical formulas, with the resulting crater morphology serving as the initial condition for the next strike, thereby enabling rapid prediction of cumulative damage throughout the sequential penetration‑explosion process. Furthermore, full‑process numerical simulations of a three‑projectile sequence were conducted to validate the theoretical predictions under several typical firing configurations. The results demonstrate that the theoretical model can accurately predict the sequential penetration depth. Although deviations in crater morphology exist for the first penetration-explosion event, the prediction accuracy improves significantly as the projectile sequence advances. The penetration depth of subsequent projectiles is notably enhanced due to concrete pre‑damage, with the greatest gain occurring inside the crater tunnel zone—yet being highly sensitive to impact‑point location—while inside the funnel zone the gain decays gradually with increasing inter‑projectile spacing. The proposed theoretical model is computationally efficient and can serve as a theoretical prediction tool for damage assessment in multi‑projectile sequential strikes.

Damage assessment and study of influencing factors in reinforced concrete beams under combined high temperature and impact
OUYANG Xin, WANG Wei, LIU Jiening, ZHOU Yongwang
 doi: 10.11858/gywlxb.20261018
PDF(64)
Abstract:
To study the impact resistance performance of reinforced concrete (RC) beams under high temperature and impact loading, a numerical model of RC beams under high temperature was established using LS-DYNA software based on existing experiments, and the validity of the numerical model was verified. Based on the numerical model, four failure modes of RC beams under high temperature and impact loading were summarized. The section damage factor D was introduced to evaluate the damage degree of RC beams under high temperature, and a predictive formula for the section damage factor in relation to impact height and temperature was fitted. The influence of temperature, beam span, impact location, and hammer shape on the impact resistance performance of RC beams under high-temperature conditions was analyzed. The results show that as the temperature increases, the mid span displacement of RC beams increases, while the impact force decreases. The larger the span to height ratio, the longer it takes for the mid span displacement to reach its peak. Under the same impact energy, the mid-span displacement caused by impact at the mid-span is greater than that caused by impact near the supports. Furthermore, hammerheads with flat contact surfaces cause greater damage to RC beams compared to those with curved contact surfaces.
Influence of Different Wire Diameter and Discharge Period on Copper Wire Electrical Explosion
ZHOU Zhangan, LU Yizhan, XIAO Bo
 doi: 10.11858/gywlxb.20261003
PDF(74)
Abstract:
Wire electrical explosion has important applications in the fields of Z-pinch, nano powder preparation and special processing. In this paper, the electrical explosion process of copper wire in vacuum was studied by numerical simulation, and the effects of wire diameter and discharge period on the electrical explosion process of copper wire were investigated. The simulation results show that the use of fine wire has higher energy efficiency than that of thick wire. In addition, the effect of electric explosion of filament is hardly affected by the speed of discharge. When the wire diameter increases, the faster the discharge, the worse the spatial uniformity of electric explosion, but on the other hand, it can also improve the energy utilization efficiency of electric explosion.
Prediction of Dynamic Mechanical Response of Materials Based on U-Net Model: Influence of Texture Representation Differences
GAO Xiang, ZHAO Dan, FANG Huiqing, WANG Jianjun, MA Shengguo, WANG Zhihua
 doi: 10.11858/gywlxb.20251280
PDF(67)
Abstract:
The utilization of neural networks for the prediction of alloy properties and the inverse design of alloy microstructures has emerged as a novel approach in the industry for understanding material performance and developing new alloys. Texture acts as a critical factor influencing microstructural evolution during alloy deformation. It is typically characterized by spatially uncorrelated discrete grain orientation Euler angles, spatial-orientation coupled Euler angles within a Representative Volume Element (RVE), or pole figures/inverse pole figures. However, identifying which texture representation method serves as the optimal input to maximize the performance of neural network models requires further investigation. Consequently, employing a modified U-Net model as the backbone architecture, this study evaluates and compares the impact of three texture representation methods - discrete Euler angles, spatial Euler angles, and pole figures-as model inputs on the overall performance of the neural network. The three trained neural network models were individually deployed to predict samples within the test set. The results demonstrate that employing pole figures as the texture representation method yields the optimal performance. Furthermore, the trained neural network models were utilized to predict the macroscopic stress-strain curves of the alloys by incorporating a one-dimensional (1D) convolutional layer at the output stage. Compared to traditional methods relying solely on fully connected layers, this modification significantly enhances the prediction accuracy of the stress-strain curves.
Impact Test and Crack Propagation Law of Rock-like Materials with Different Joint Geometric Parameters
FEI Honglu, CHEN Liangyu, YANG Pengliang, DING Wen, YANG Shitao, HU Gang, ZHOU Linli
 doi: 10.11858/gywlxb.20261025
PDF(59)
Abstract:
The crack-propagation behavior and dynamic response characteristics of artificial jointed rock-like materials under impact loading were investigated. Drop-weight impact tests were conducted using a self-developed apparatus, in which joint aperture, inclination angle, number of joints, and the distance between the joint and the loading surface were taken as geometric parameters. LS-DYNA was used to establish numerical models of drop-weight impact for vertical combined joints composed of multiple parallel vertical joints and for vertical-horizontal combined joints consisting of one set of vertical joints intersecting with one set of horizontal joints. The influence of joint geometric parameters on the failure mode, crack-propagation path, and stress-strain response of the jointed rock-like materials was analyzed. Based on the linear elastic fracture mechanics theory under the plane stress assumption, the relationship between the strain components at the joint end and the stress intensity factors was derived. The test results show that as the joint aperture increased from 0.5mm to 0.9mm, the maximum impact load decreased from 35.16kN to 22.07kN; as the joint inclination angle increased from 30° to 60°, the maximum impact load decreased from 47.17kN to 29.57kN; and as the number of joints increased from one to three, the maximum impact load decreased from 36.31kN to 23.69kN, all showing a decreasing trend. When the distance between the joint and the loading surface increased from 50mm to 90mm, the maximum impact load first increased from 38.19kN to 41.75kN and then decreased to 40.88kN, showing a nonlinear increase-decrease relationship. The peak strain at the measurement points was negatively correlated with the joint aperture, inclination angle, number of joints, and joint-loading surface distance. The numerical simulation results indicate that the overall stability of the model decreases significantly with an increasing number of vertical joints; in the vertical-horizontal combined joint model, enlarging the vertical joint aperture reduces the degree of bottom damage, and when the horizontal joint is located above the vertical joint, the model becomes more susceptible to instability and failure.
Lattice Dynamics of Mercury under High Pressure
LIU Peiyuan, ZHAO Bohao, WANG Lijuan, LIN Chuanlong, YANG Liuxiang, GOU Huiyang
 doi: 10.11858/gywlxb.20261011
PDF(67)
Abstract:
Mercury is a metal with special physical and chemical properties. High pressure can significantly alter its crystal structures and interatomic interactions. However, vibration-related experimental data of mercury under high pressure have been scarce for a long time, which limits the in-depth understanding of its high-pressure phase transition mechanisms. To fill this research gap, this study based on a low-wavenumber high-pressure Raman experimental platform, conducted in-situ high pressure Raman measurements using diamond anvil cell (DAC) technology, and combined with theoretical calculation methods to systematically investigate the laws of structural transformation and vibrational mode evolution of mercury under high pressure. The study successfully detected the vibrational signals of two solid phases of mercury under high pressure for the first time, obtained the Raman spectroscopy experimental results of mercury under high-pressure conditions, revealed the characteristics of pressure-dependent vibrational frequencies, and calculated the relevant thermodynamic parameters of two phases. The experimental data is highly consistent with the reported structural phase transitions of mercury, filling the long-standing experimental gap in this field and providing key experimental support for an in-depth understanding of the Phase transition behaviors and related mechanisms of mercury under high pressure.
Microscopic mechanisms of plastic deformation in high-entropy carbide (Zr0.2Hf0.2Ti0.2Nb0.2Ta0.2)C under quasi-isentropic and ramp compression
FENG Lanxi, LI Chuanying, LI Wanghui, ZHANG Xiaoqing, YAO Xiaohu
 doi: 10.11858/gywlxb.20261008
PDF(74)
Abstract:
High-entropy carbides (HECs), characterized by pronounced chemical disorder and lattice distortion, exhibit exceptional strength-toughness synergy and are promising candidates for impact protection and high-temperature structural applications. However, their microstructural evolution and stress response under extreme conditions, such as high stress and strain rates, remains poorly understood. In this work, a high-accuracy machine-learning interatomic potential is employed to investigate the representative multi-principal carbide (Zr0.2Hf0.2Ti0.2Nb0.2Ta0.2)C (HEC) through large-scale molecular dynamics (MD) simulations. To elucidate how multiple "high-entropy" effects govern atomic-scale plasticity in high-entropy carbides (HECs), large-scale molecular dynamics simulations are carried out to explore their response under quasi-isentropic compression and ramp-wave loading along three principal crystallographic orientations: [0 0 1], [0 1 ¯1], and [1 1 1]. The results demonstrate that the high-entropy effects profoundly reshape the initiation of plasticity, the competition among slip systems, and the localized deformation modes in HEC. Local stress fluctuations and lattice distortions enhance transient Bain-type stacking rearrangements within the sublattice and promote the synergistic activation of multiple slip systems. This leads to a transformation of shear band formation from isolated nucleation to a network-like propagation. Complementary first-principles calculations reveal that the carbon vacancy formation energy in high-entropy ceramics is significantly reduced compared to their single-component carbides. This reduction of vacancy formation energy facilitates preferential displacement of carbon atoms and their participation in shear band nucleation during compression. Furthermore, the comparison between different loading paths highlights the complexity of the high-entropy effects’ response. The quasi-isentropic loading path helps to unveil the intrinsic deformation mechanisms governed by the high-entropy effect itself, whereas the stress gradients inherent in ramp-wave loading couple with the high-entropy effect, leading to a premature triggering and intensification of plastic localization.
Pressure-Tuned Superconductivity in FePSe3 Thin Films
CHENG Yi, LI Meilun, XIAO Hong, LIN Chuanlong
 doi: 10.11858/gywlxb.20261001
PDF(80)
Abstract:
Two-dimensional layered materials constitute a unique class of compounds in which strongly covalent or ionic atomic planes are stacked via van-der-Waals forces. This weak interlayer coupling allows the thickness to be precisely tuned down to few-layer or monolayer, giving rise to a rich spectrum of dimensionality-dependent physical properties. In this work, we take the prototypical van der Waals layered compound FePSe3 as a model system and, by combining mechanical exfoliation with high-pressure techniques based on a diamond anvil cell (DAC), systematically investigate the electrical transport properties of both bulk FePSe3 and thin layers with different thicknesses under pressure. We focus on the combined effects of external pressure and reduced dimensionality on the normal-state transport behavior and superconductivity. Experimental results show that bulk FePSe3 exhibits pressure-tuned superconductivity, with the superconducting transition temperature <italic>Tc</italic> reaching a minimum around 15 GPa, accompanied by a concurrent minimum in the Hall coefficient <italic>R</italic>H. This behavior is consistent with previous reports on bulk materials, and suggests that pressure may induce a Fermi surface reconstruction. Compared to the bulk, the thin-layer FePSe3 samples show a suppressed superconducting state, characterized by a reduced <italic>Tc</italic>, and a monotonic decrease in <italic>R</italic>H with increasing pressure. This indicates that two-dimensional confinement in thin flakes suppresses the occurrence of Fermi surface reconstruction. These findings provide key experimental evidence for understanding the pressure-driven evolution of the electronic states in FePSe3.
Theoretical Calculation of Combustion Performance of Capsule-type Gun Propellant
LIANG Hao, LIANG Jinghao, DAI Zengjie, XU Bo, SONG Baolu, PANG Yu, HU Pin
 doi: 10.11858/gywlxb.20261009
PDF(65)
Abstract:
To address the critical requirement for energy release control under high-pressure conditions in gun bores, and to resolve the conflict between high loading density and combustion progressivity, a novel end-capped tubular combined charge structure containing spherical propellants based on the "capsule" structure concept is proposed. This design utilizes the tubular shell and end-capped discs to form a slow-burning "capsule shell", with high specific surface area discrete spherical propellants filled inside acting as fast-burning "capsule agents". Through physical isolation and geometric burn-through of the outer shell, a high-progressivity combustion mode characterized by "slow burning followed by fast burning" is established. A theoretical combustion model incorporating the form function and gas generation intensity (Γ) is developed to quantitatively analyze the influence of component geometric parameters on combustion performance. Calculation results indicate that the combined charge exhibits typical "step-wise" double-peak characteristics during the combustion process. The burning thickness ratio of the end-capped disc to the tubular shell (rp-k) is critical for regulating the phase of high-pressure gas intervention into the internal combustion, thereby determining the onset timing of the Γ surge. The burning thickness ratio of the spherical propellant to the tubular shell (rq-k) directly governs the magnitude of combustion progressivity. When rq-k<0.5, the spherical propellants with high specific surface area burn out before the tubular shell within the high-pressure field, triggering a drastic surge in Γ; moreover, a smaller rq-k corresponds to a larger initial burning surface area, leading to a more significant gain in the Γ peak. Conversely, when rq-k>0.5, the tubular shell burns out before the spherical propellants, resulting in a step-wise drop in intensity at the late stage of combustion (Ψ>0.95). Furthermore, an increase in loading density further amplifies this surface area augmentation effect induced by geometric phase transition. Theoretically, this combined charge with "capsule" characteristics possesses high combustion progressivity, offering a novel approach for the development of combined gun charges.
Elastic wave velocities of lawsonite and its implications for seismic velocity anomalies in subduction zones
MA Junsheng, LIU Chuanjiang, WANG Duojun, CAI Nao, ZHANG Rui
 doi: 10.11858/gywlxb.20261006
PDF(63)
Abstract:
Lawsonite is a key hydrous mineral in subducting oceanic crust, with an H₂O content of up to 11.5 wt.% and volume fractions that can reach ~50 vol.%, and it is stable over a wide pressure–temperature range. Constraining the high-pressure elastic properties of lawsonite is therefore essential for understanding seismic velocity structures in subduction zones and the deep water cycle. In this study, dense polycrystalline lawsonite was synthesized by hot-pressing natural lawsonite powder at 7 GPa and 1073 K for 2 hours. Ultrasonic interferometry was then used to measure P- and S-wave velocities of lawsonite up to 8 GPa at room temperature, from which the elastic moduli and their pressure derivatives were determined. The results show that both elastic wave velocities and moduli increase with pressure, yielding KS0=120(1)GPa, G0=51(1)GPa, Ks'=5.7(1), G'=0.9(1). Using these elastic properties, together with previous mineral-physics results, we constructed a lawsonite-bearing eclogite model, which indicates that the presence of lawsonite can provide a plausible explanation for the seismic velocity anomalies observed at depths of ~60–90 km within intermediate-temperature subducting slabs such as Nicaragua.
Structural Evolution, Dual Role of Hydrogen and Superconductivity in Lithium-Rich Li-N-H Compounds under High Pressure
LI Qiuyue, HAN Shuai, YANG Guochun
 doi: 10.11858/gywlxb.20261017
PDF(70)
Abstract:
Using first-principles calculations combined with CALYPSO structure prediction method, we systematically investigated the crystal structures, electronic properties, and superconducting behavior of the Li-N-H ternary system under the pressure ranging from 100 to 300 GPa. Six thermodynamically stable or metastable Li-rich compounds were identified: C2 Li2NH, P212121 Li2NH, I4/mmm Li3NH, Imm2 Li8NH, Cm Li9NH and P2/c Li10NH. The results reveal a distinct evolution in the chemical role of hydrogen with increasing lithium content. In phases with lower Li content, H atoms tend to form covalent bonds with N atoms, thereby achieving a stable closed-shell electronic configuration. As the Li content increases, H atoms progressively occupy lattice interstitial sites, acting as electron acceptors that trap excess interstitial anionic electrons (IAEs). This transformation effectively tunes the quantity, degree of localization, and spatial topology of IAEs. Correspondingly, C2 Li2NH, P212121 Li2NH and I4/mmm Li3NH are insulating non-electrides, while Imm2 Li8NH, Cm Li9NH and P2/c Li10NH are electrides, which exhibit metallic behavior. Notably, P2/c Li10NH demonstrates superconductivity with a predicted transition temperature of 4.8 K at 100 GPa, mainly originating from the strong electron-phonon coupling between H-p orbital electrons and low frequency phonon modes dominated by Li atoms. This work elucidates the dual functional role of hydrogen in high-pressure Li-N-H systems—from covalent N-H bonding coordination to IAEs capturing—and provides theoretical insights for the rational design of novel high-pressure electrides and superconducting materials.
Statistical Characteristics of Spallation Based on Stochastic Numerical Simulation
GUAN Youhao, ZHANG Hao, PEI Xiaoyang
 doi: 10.11858/gywlxb.20251290
PDF(70)
Abstract:
This work integrates stochastic theory with a phase-field model for spallation in ductile metals. By assigning the initial yield strength four distinct random distributions to characterize the random distribution of material defects and employing an explicit dynamic solver, the entire process of spall damage—from gradual evolution to instability and coalescence—was successfully simulated. The simulation results were validated through plate impact experiments and triangular wave loading experiments. These validations revealed the relationship between the heterogeneity of material yield strength and both the spall strength and the number/area of damage zones. The results indicate a negative correlation between the standard deviation of the initial yield strength and the spall strength, which holds for both single and multiple spall scenarios in ductile metals. For single spallation, regardless of the initial distribution of yield strength, the resulting spall strength follows a normal distribution. For multiple spallation, the number of initially nucleated damage zones increases linearly with the standard deviation, while the size of these zones follows a Weibull distribution. Under the same initial random distribution, the number of damage zones evolves over time, showing a trend of initial slow growth, subsequent acceleration until saturation, and a final decline after saturation. This trend corresponds to the typical process of damage evolution involving nucleation and coalescence during spallation.
Progress on Cross-Scale Design and Machine Learning Prediction of Penetration Resistance of Hybrid Fiber Reinforced Concrete
YU Xiaofeng, LUO Jianlin, WEN Yulei, ZHU Min, MA Minglei, LIU Chao, LIAN Chunming, CHEN Fengwei
 doi: 10.11858/gywlxb.20251258
PDF(95)
Abstract:
Hybrid Fiber-Reinforced Concrete (HFRC) significantly enhances penetration resistance through multi-scale fiber hybridization and multi-stage energy dissipation mechanisms. Compared to single-type fiber-reinforced concrete, HFRC exhibits superior dynamic strength, energy absorption capacity, and crack resistance, establishing it as a key structural material in military protective engineering. This paper systematically reviews recent advances in cross-scale design and machine learning (ML) prediction of the penetration resistance of HFRC. The analysis begins by examining how different fiber combinations generate cross-scale synergistic effects that collectively improve the dynamic strength and anti-penetration capacity of HFRC. Subsequently, the mechanisms and multi-scale transmission pathways through which nanomaterials enhance the crater resistance and anti-spalling capacity of HFRC by strengthening the matrix and interfaces are examined. Furthermore, this review elucidates how multi-scale structural characteristics such as fiber distribution, orientation, and interfacial bonding synergistically govern the evolution of penetration-induced damage and the resulting failure patterns. Finally, the predictive efficacy of ML models for penetration resistance of HFRC is evaluated, along with potential integration pathways between ML and traditional numerical simulation.
Research on parameter optimization of corrugated Whipple protective structure under hypervelocity impact
GUO Jiaao, YANG Qiuzu, LIU Xiaochuan, YIN Yunfei, LI Zhiqiang
 doi: 10.11858/gywlxb.20251276
PDF(104)
Abstract:
The geometric configuration of corrugated Whipple shields significantly influences their protective capability against hypervelocity impacts. To optimize the performance of corrugated Whipple shield structures under hypervelocity impact conditions, an integrated optimization method combining the finite element-smoothed particle hydrodynamics (FE-SPH) coupled algorithm with orthogonal experimental design was proposed. A reliable numerical simulation model was established, and the Z-axis momentum density was introduced as an evaluation index for protective performance. The effects of three geometric parameters-corrugation thickness, span, and angle-on the shielding effectiveness were systematically investigated. Orthogonal test results indicated that the order of influence of these factors was: thickness > angle > span. Further two-factor refined experiments were conducted, and a quadratic polynomial model was developed to identify the optimal combination of geometric parameters. The optimized configuration improved the protective performance by 33.72% compared to a flat panel. The study confirms that the optimized corrugated structure effectively promotes projectile fragmentation and debris cloud dispersion, facilitating three-dimensional redistribution of momentum, thereby significantly enhancing the shield's protective performance. This research provides a theoretical basis and a parameter optimization pathway for the design of spacecraft protective structures.
Synthesis and Characterization of P-Doped Diamond Crystals in the FeNiCo-C System under High Pressure and High Temperature
ZHANG Haobo, HU Meihua, LI Shangsheng, LIU Di, HE Shasha, LI Xiaoxiao, WANG Zhenyang
 doi: 10.11858/gywlxb.20251285
PDF(87)
Abstract:
To investigate the effects of phosphorus doping on diamond crystal growth, diamond single crystals doped with phosphorus were synthesized along the (111) plane using the temperature gradient method. The experiments were conducted under conditions of 5.5 GPa and 1300 ℃, with Fe₃P added into the FeNiCo-C system. The synthesized diamond samples were characterized by Fourier-transform infrared spectroscopy (FTIR), Raman spectroscopy, photoluminescence (PL) spectroscopy, and X-ray photoelectron spectroscopy (XPS). With increasing Fe₃P addition, the diamond color gradually became lighter. The crystal morphology changed from octahedral to hexoctahedral. The nitrogen impurity content in the diamonds showed a decreasing trend as more Fe₃P was added. This occurs because the addition of Fe₃P alters the catalyst properties, increasing the nitrogen solubility of the catalyst. Thus, fewer nitrogen atoms enter the diamond lattice. Phosphorus doping increases internal stress and induces lattice distortion in the diamond crystal, degrading its quality. This conclusion is supported by the shift and broadening of the Raman peak. The incorporation of phosphorus atoms inhibits the formation of NV⁻ centers in diamond crystals. XPS results confirm the successful incorporation of phosphorus into the diamond lattice. This study provides useful insights for understanding the synthesis mechanism of phosphorus-doped diamond crystals. It also supports potential applications of phosphorus-doped diamond crystals.
Optimization of Borehole Spacing and Decoupling Coefficient for Presplitting Blasting in Water-Bearing Borehole
SHEN Zewei, LIU Haoshan, ZHANG Zhiyu, HUANG Yonghui, HE Defu
 doi: 10.11858/gywlxb.20251292
PDF(88)
Abstract:
In water-filled borehole presplitting blasting,the incompressibility and high wave impedance of the water medium significantly alter the pathways of explosive energy transmission and the rock-breaking mechanisms.As a result, traditional parameter design methods developed for air-filled boreholes often lead to high overbreak ratios and excessive damage to the retained rock mass in water-bearing strata.Taking the water-rich slope of the Jianshan phosphate mine as the engineering background,this study establishes a coupled smoothed particle hydrodynamics–finite element (SPH–FEM) numerical model to systematically investigate the propagation characteristics of blast-induced stress waves,rock mass damage evolution,and crack propagation behavior under different borehole spacings and decoupling coefficients.The results indicate a strong correlation between the superposition of stress waves from adjacent boreholes and the coalescence of presplitting cracks.When the borehole spacing is 1.4 m and the decoupling coefficient is 2.34,the average crack propagation length reaches 49.48 cm,enabling the formation of regular and continuous through-going presplitting cracks while effectively suppressing excessive crushing around the borehole wall and the development of secondary cracks.Field tests further validate the reliability of the numerical simulations:under the optimized parameters,the half-hole rate of water-filled borehole presplitting blasting increases to 85%,and the acoustic reduction rate decreases by 18% compared with conventional blasting,demonstrating favorable damage control performance and crack-forming effectiveness under water-bearing conditions.The findings provide a useful reference for presplitting blasting parameter design in complex hydrogeological environments.
A WMA-SVM Model for Slope Stability Prediction
SUN Huafen, RAO Hui, HOU Kepeng, WANG Honglin, WANG Zeqi
 doi: 10.11858/gywlxb.20251241
PDF(87)
Abstract:
To enhance the prediction accuracy of data-driven models in slope stability classification, this study proposes a hybrid intelligent model (WMA-SVM) that integrates a novel Whale Migration Algorithm (WMA) with a Support Vector Machine (SVM). First, a heterogeneous dataset of slope cases from diverse engineering backgrounds was constructed. To address its significant class imbalance, a combined strategy using the Synthetic Minority Over-sampling Technique (SMOTE) and the Local Outlier Factor (LOF) algorithm was adopted to generate a high-quality balanced dataset. Subsequently, the WMA algorithm, which demonstrated superior optimization performance on eight benchmark test functions, was employed to optimize the hyperparameters of the SVM adaptively. Evaluation results show that the proposed WMA-SVM model significantly outperforms all benchmark models across various performance metrics. Moreover, based on the Permutation Feature Importance (PFI) method, the unit weight (γ), slope angle (β), and internal friction angle (φ) were identified as the most critical features influencing the classification outcomes for this dataset. Finally, the model's generalization capability was further validated through eight independent engineering case studies, revealing a high consistency between the predictions and the actual stability states. This research provides a modeling framework with considerable generalization potential for the intelligent analysis of slope stability.
Equivalent Bird-Strike Test Method and Fixture Design for the Trailing Edge of Aero-Engine Composite Fan Blades
SI Wulin, LI Wenhao, JIANG Xiaowei, LI You, ZHAO Zhenqiang, ZHANG Chao
 doi: 10.11858/gywlxb.20251271
PDF(97)
Abstract:
In order to investigate the response and damage behavior of aero-engine composite fan blades under bird-strike events, an equivalent bird-strike test method is proposed in which component-level flat plate specimens are used to replace full-scale fan blades, with the aim of reproducing the trailing-edge delamination observed in full-scale blades during bird strikes by means of component-level flat plate tests. By carrying out bird-strike tests and numerical simulations of flat plate specimens under different clamping schemes, the impact response characteristics of the specimens and the initiation and propagation of delamination damage under each scheme are systematically analyzed; on this basis, a component-level equivalent test method capable of effectively simulating trailing-edge delamination during the blade bird-strike process is proposed, and the baseline impact conditions that can induce one-sided trailing-edge delamination in typical composite laminates, including impact height, impact velocity and bird-cut ratio, are determined. Moreover, by comparing experimental and numerical results under different impact conditions, the accuracy of the numerical model is verified. On the basis of the validated numerical model, numerical simulations are further used to analyze the sensitivity of the established equivalent test method to the impact parameters, and to quantify the influence of impact height, impact velocity and bird-cut ratio on the impact response of composite flat plates. The results show that the equivalent test method proposed in this paper can reproduce, through composite flat plate tests, the local displacement response and delamination damage modes of full-scale blades under bird strikes, and that the experimental results exhibit good robustness.
YANG Xigui, LAI Shoulong
 doi: 10.11858/gywlxb.20251255
PDF(90)
Abstract:
Rhombohedral C60 holds significant potential for applications in fields such as two-dimensional materials and catalysis; however, the synthesis of high-purity, high-quality rhombohedral C60 remains challenging. In this study, rhombohedral C60 was successfully synthesized at 6 GPa and 650 ℃. Characterization by X-ray diffraction, Raman spectroscopy, X-ray photoelectron spectroscopy and spherical aberration-corrected transmission electron microscopy confirmed the obtained sample as a high-purity two-dimensional rhombohedral structure. The effects of pressure and temperature (6-10 GPa, 650-800 ℃) on the polymerization of C60 were investigated, clarifying the phase boundary between the rhombohedral phase and disordered amorphous carbon clusters. Temperature-dependent Raman spectroscopy revealed that the rhombohedral C60 polymer remains stable up to approximately 350 °C, beyond which it depolymerizes, reverting to the pristine face-centered cubic C60 molecules. This work defines a clear processing window for the synthesis of high-quality rhombohedral C60 and establishes an experimental foundation for its further application in functional materials.
First-Principles Study of the Effects of Phase Transitions and Decomposition on the Lattice Thermal Conductivity of Calcite
HONG Zheng, ZHU Yongqiang, XIONG Yuanmeng, LU Cheng, HE Kaihua
 doi: 10.11858/gywlxb.20251228
PDF(93)
Abstract:
The lattice thermal conductivity (κlatt) of minerals plays a critical role in controlling heat flow and temperature distribution in the Earth's interior. Calcite, primarily composed of calcium carbonate (CaCO₃), can be subducted into the deep Earth and serves as an important carbon source. As pressure and temperature conditions change with depth, CaCO₃ undergoes phase transitions and thermal decomposition, which significantly affect its physical properties. In this study, we investigate the effects on thermal conductivity of calcite induced phase transitions and thermal decomposition using first-principles calculations combined with lattice dynamics. Our results show that the calcite I → calcite II phase transition leads to a reduction in thermal conductivity, whereas subsequent phase transitions at higher pressures result in increase. The thermal conductivity of aragonite and post-aragonite increases nearly linearly with pressure increasing, and the latter exhibiting a stronger pressure dependence. Upon thermal decomposition, the CaO exhibits significantly higher thermal conductivity than that of calcite, which may enhance local heat transfer. Analysis of relevant thermodynamic parameters indicates that the changes in thermal conductivity induced by phase transitions and decomposition are collectively determined by phonon group velocity and anharmonic scattering rates.
Velocity Variation Law of Two Projectiles in Staggered Sequential Penetration into Concrete Targets with Limited Thickness
XU Baowen, ZHANG Dingshan
 doi: 10.11858/gywlxb.20251244
PDF(89)
Abstract:
To investigate the effects of dislocation distance and projectile diameter on the velocity variation of the second projectile during the sequential penetration of a concrete target, a theoretical model was developed to characterize the energy loss and velocity change during the displaced penetration process. Validation experiments were designed, and a comparative analysis was conducted among theoretical predictions, experimental data, and numerical simulations. The results indicate that displaced sequential penetration reduces the velocity decay of the second projectile, thereby enhancing its penetration depth. As the dislocation distance increases, the beneficial influence of the first projectile on the second projectile’s velocity retention diminishes. Beyond a critical dislocation distance, this effect becomes negligible. A larger diameter of the first projectile corresponds to a greater critical dislocation distance. Under test conditions involving penetration of a 1 m thick C40 concrete target at an initial velocity of 600 m/s, the critical dislocation distances for projectile diameters of 50 mm, 80 mm, and 100 mm were approximately 8d, 10d, and 14d, respectively. The maximum deviations between theoretical predictions and experimental results for the second projectile’s velocity were about 7.1%, while numerical simulations deviated by approximately 3.8% from the experimental data.
Wave-Cutting Efficiency and Mechanism of Single-Tube Multi-Row Hole Bubble Curtain
DU Mingran, LI Jirui, JIN Cong, CENG Huilian, TAN Caiyong
 doi: 10.11858/gywlxb.20251231
PDF(86)
Abstract:
To further optimize the wave-cutting efficiency of bubble curtain, field tests on underwater explosion shock wave attenuation using single-tube multi-row bubble holes and high-speed photography observations of bubble curtain morphology were designed. Additionally, the numerical calculation model for equivalent thickness of bubble curtain was studied using AUTODYN software. The results indicate that, under the same air flow rate, the number of bubble hole rows is a critical factor affecting wave-cutting efficiency. At a detonation center distance of 12m, the wave-cutting efficiencies for 1, 2, and 3 rows of holes are 89.92%, 97.25%, and 96.41%, respectively. At different detonation center distances, the cutting efficiency of the two-row hole bubble curtain is the best, and the cutting efficiency is all greater than 95.00%. Both the thickness and density of the bubble curtain are maximized with 2 rows of holes, and the thickness of the bubble curtain is the key factor determining wave-cutting efficiency. The equivalent thickness fitting formula established by combining experiments and simulations has high reliability, and the simulation model exhibits high accuracy. It is suggested that similar projects adopt single-tube 2-row hole bubble curtain to achieve convenient, efficient and low-cost wave-cutting efficiency.
Raman Scattering Study of Lattice Dynamics and Phase Transitions in Layered Perovskite Sr2Ta2O7 Ceramics under High Pressure
QIAN Chao, QI Wenming, Abliz Mattursun, HU Qingyang, WANG Yuanyuan, DONG Hongliang, CHEN Bin
 doi: 10.11858/gywlxb.20251269
PDF(98)
Abstract:
Strontium tantalate (Sr₂Ta₂O₇) is a ceramic with an orthogonal Cmcm space group phase. Due to its potential applications in the field of multiferroic materials, it has become a research hotspot in recent years. However, the regulatory mechanism and phase transition behavior of hydrostatic pressure on its complex lattice structure remain unclear, which, to some extent, limits the in-depth understanding of the "structure-property" relationship of this material. This study systematically investigates the lattice dynamic response characteristics of orthogonal Cmcm Sr₂Ta₂O₇ under high pressure up to 30 GPa using in-situ high-pressure Raman spectroscopy, marking the highest pressure study conducted on this system to date. The results indicate that when the pressure reaches 5 GPa, significant changes occur in the Raman vibrational modes of the material, a phenomenon attributed to a structural phase transition induced by symmetry breaking, corresponding to a transition from a commensurate phase to an incommensurate phase, consistent with previous research findings on Sr₂Ta₂O₇. As the pressure further increases to 20 GPa, a second phase transition may occur, which is identified as a first-order phase transition closely related to lattice disordering. However, the specific crystal structure of this high-pressure phase remains to be further confirmed in future studies. The Raman spectroscopy analysis suggests that the structural distortion of this high-pressure phase may follow a transformation pathway from orthogonal to monoclinic.
Coupling Mechanism of Wall Protection Blasting and Notched Blasting
GONG Yue, SU Hong, LIU Buqing, PAN Yan, WANG Cunguo, SUN Jinshan
 doi: 10.11858/gywlxb.20251238
PDF(80)
Abstract:
In order to improve the blasting forming accuracy and the protection of surrounding rock under deep and complex geological conditions, comparative experiments of protective blasting and protective notched-coupling blasting were conducted using a digital laser dynamic caustic experimental system and PMMA specimens. The crack propagation mechanisms and mechanical response characteristics of the two blasting modes were systematically investigated. The results show that, in the protective notched-coupling blasting, the main crack propagates stably along the preset notched direction, demonstrating an excellent directional control effect. Meanwhile, the secondary peak value of the crack-tip stress intensity factor is significantly higher than that of the single protective blasting, indicating a stronger dynamic stress concentration effect. The notched-coupling blasting also exhibits a higher overall crack propagation velocity with slower attenuation at later stages, reflecting enhanced persistence and stability. In addition, this blasting mode effectively reduces the length and number of cracks on the protected wall side, thereby providing better rock-mass protection. Overall, the protective notched-coupling blasting optimizes crack propagation behavior and improves the directionality and energy utilization efficiency of blasting through the combined mechanisms of notch guidance and energy re-concentration. These findings provide theoretical support and technical guidance for precision blasting design and engineering applications in deep rock masses.
A Dynamic Spherical Cavity Expansion Model for Ceramics Considering Shear-Dilatancy
LI Xiao, LIANG Xuan, WEN Heming
 doi: 10.11858/gywlxb.20251242
PDF(83)
Abstract:
The cavity expansion theory is often used to predict the penetration resistance of a target against a projectile. A dynamic spherical cavity expansion model for ceramic materials is suggested by considering shear-dilatancy effect through introducing a dilatancy-kinematic relation in comminuted region. The comminuted region is further divided into a linear comminuted region (satisfying the Mohr-Coulomb failure criterion) and a saturated comminuted region (satisfying the maximum shear strength), depending on whether the shear strength reaches its maximum (plateau). Firstly, equations for calculating radial stress at cavity surface are derived. Secondly, numerical simulations of cavity expansion process in ceramics at different expansion velocities are conducted. Finally, the effects of key parameters such as compressive strength and density on the cavity surface radial stress are discussed. It is shown that the model predictions for cavity radial stress and interface velocities of cracked and comminuted regions are in good agreement with numerical simulations. It is also shown that compressive strength plays a dominant role in enhancing cavity radial stress and the influence of density increases with increasing cavity expansion velocity.
Dynamic Plastic Deformation Mechanism of 301 Stainless Steel at Low Temperatures
WANG Pengfei, HUANG Tingting, CHEN Meiduo, ZHAN Junlan, TIAN Jie, XU Songlin
 doi: 10.11858/gywlxb.20251246
PDF(77)
Abstract:
Deep space exploration faces challenges from extreme temperatures and complex high-speed operating environments, placing higher demands on the low-temperature impact resistance of materials. In this study, a low-temperature Hopkinson bar impact experimental device under a vacuum liquid helium environment was developed to achieve dynamic loading of materials under ultra-low temperature conditions. The dynamic mechanical response of 301 stainless steel produced by two rolling processes was investigated under the combined effects of low temperature (30K-298K) and high strain rates (4000 s⁻¹-5000 s⁻¹). Experimental results show that the yield strength of both materials exhibits a significant negative correlation with temperature and a positive correlation with strain rate. The unidirectionally rolled samples displayed an anomalous increase in toughness at 77K. The study indicates that the unidirectional rolling process induces a higher content of martensitic phase, thereby endowing the material with greater strength. Microstructural characterization results reveal that the anomalies in macroscopic mechanical behavior stem from the competition of deformation mechanisms; at room temperature, the samples mainly exhibit a toughness fracture mechanism dominated by ductile dimples, whereas at low temperatures, they transition to a brittle fracture mode dominated by quasi-cleavage. On this basis, the Johnson-Cook constitutive model was used to fit the mechanical properties, showing good consistency with the experimental results. This research provides important experimental methods and theoretical support for the dynamic strength and toughness design of metallic materials under extreme low-temperature impact conditions.
The Role of Gradient Structure in the Integrated Performance of PDC Cutters
GAO Jun, ZHANG Zhicai, HOU Zhiqiang, WANG Chao, LI Hao, YANG Yikan, YANG Jiao, FANG Rui, TANG Yao, WANG Haikuo
 doi: 10.11858/gywlxb.20251233
PDF(81)
Abstract:
Addressing the urgent demand for high-performance polycrystalline diamond compact (PDC) cutters in deep/ultra-deep oil and gas exploration, this study optimized the PDC synthesis formulation through orthogonal experimental design. Under high pressure conditions (8.5 GPa and 1750 °C), we successfully fabricated both conventional homogeneous mixed PDC cutter (H-PDC) and gradient-structured PDC cutter (G-PDC) featuring a "fine-grained work layer/coarse-grained transition layer" structure. Microstructural characterization reveals that the gradient structure facilitates uniform distribution of cobalt binder, suppresses cobalt aggregation, enhances interlayer interfacial bonding, and generates higher residual compressive stress. The cobalt content in the G-PDC work layer is 9.16 wt.%. After acid leaching for cobalt removal, the cobalt content decreased to 2.49 wt.%. Performance evaluations demonstrate that G-PDC achieves a wear resistance lifespan of 920 passes, superior to H-PDC (800 passes). The average impact toughness of G-PDC reaches 740 J, representing approximately 107% improvement over H-PDC. Furthermore, the gradient structure alleviates thermal expansion mismatch, increasing the thermal stability temperature by about 30 °C. This research confirms that combining high pressure synthesis technology with gradient structural design can synergistically enhance the wear resistance, impact toughness, and thermal stability of PDC cutters, providing a viable pathway for developing next-generation superhard composites for extreme working conditions.
The Influence of Cottonseed Oil Content on the Rheological Properties and Anti-Vibration Performance of Site-Mixed Emulsion Matrix
YUE Xing, HE Zhiwei, HUANG Zhenyi, YUE Jiawei, HU Qianhao, LI Yuanlong
 doi: 10.11858/gywlxb.20251226
PDF(70)
Abstract:
To investigate the influence of cottonseed oil content on the rheological properties and anti-vibration performance of site-mixed emulsion explosive matrix, samples of site-mixed emulsion explosive matrix with different cottonseed oil contents were prepared. A rotational rheometer, HY-5A rotary speed-regulating vibrator, and water dissolution method were used to study the rheological properties and anti-vibration performance of the site-mixed emulsion matrix. The results show that when the mass fraction of cottonseed oil is not more than 2.5%, the viscosity of the site-mixed emulsion matrix increases gradually. In terms of the temperature environment of the pipeline of on-site mixing trucks, the viscosity of the emulsion matrix can meet the pumping requirements. The elastic modulus and cohesion increase gradually and then remain nearly stable. The anti-vibration performance first increases and then decreases, and when the cottonseed oil content is 2%, the anti-vibration performance is the best.
Research Progress on Two-Dimensional Diamond
MING Jiaxin, LI Jiayin, CHEN Yabin
 doi: 10.11858/gywlxb.20251248
PDF(88)
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.
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(76)
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>
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(70)
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(60)
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.
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(65)
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(89)
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(81)
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.

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(92)
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.
Cover
2026, 40(7)  
PDF(0)
Abstract:
2026, 40(7): 1-2.  
HTML PDF(0)
Abstract:
Mechanical Behavior of Lightweight Composite Protective Materials and Structures
2026, 40(7): 070101.   doi: 10.11858/gywlxb.20261119
HTML PDF(0)
Abstract:
Research Progress on Design Strategies and Impact Resistance of Heterogeneous Cellular Structures Material
LI Shiqiang, LI Zihao, WANG Zhihua, LU Guoxing
2026, 40(7): 070102.   doi: 10.11858/gywlxb.20261041
HTML PDF(66)
Abstract:

As lightweight and high-strength functional-structural integrated materials, cellular structural materials are widely applied in aerospace, automotive manufacturing, and biomedical fields. However, traditional single-configuration cellular materials (e.g., honeycomb structures and point-lattice lattices) gradually exhibit performance limitations under complex conditions such as impact shock waves, multi-directional impacts, or nonlinear deformations. Against this backdrop, heterogeneous cellular structure material (HCSM) have emerged as a research hot pot in impact protection. This paper systematically reviews recent design strategies and impact resistance performance of HCSM. HCSMs are primarily categorized into two types: topological configuration heterogeneity (including complementary and enhanced fusion) and material heterogeneity (e.g., filling with foam materials and shear-thickening materials). Through innovative “functional fusion” approaches, they overcome the performance bottlenecks of single-configuration cellular materials. The study further elucidates the synergistic reinforcement effects and deformation mechanisms of HCSM under impact loads, while analyzing their intrinsic mechanisms for improving energy absorption efficiency, stiffness, and stability. Despite significant progress in HCSM research, challenges remain in connectivity optimization, additive manufacturing process compatibility, complex condition validation, and multifunctional integration. Going forward, the integration of artificial intelligence and machine learning technologies holds promise for achieving end-to-end optimization of HCSMs from design to manufacturing, thereby providing new directions for developing next-generation high-performance impact-resistant structural materials.

Research Progress on Design Strategies and Mechanical Behaviors of Self-Locking Structure
XIONG Jian, YAN Chengrui, CHEN Zongbing
2026, 40(7): 070103.   doi: 10.11858/gywlxb.20261028
HTML PDF(64)
Abstract:

The self-locking structure achieves interlocking property through ingenious design of the connection mode between cells, which enables the cells to lock with each other without the need for any additional constraints. The self-locking structure possesses significant advantages, such as light weight, portability, rapid assembly, and disassembly. Therefore, this structure is widely applied in various fields, such as shock resistance and explosion prevention. Self-locking structures exist in many structures in nature. The design concepts of self-locking structures are introduced from three aspects: the inspiration from biomimetic self-locking structures, the energy absorption mechanism of periodic structures, and failure of shear bands in periodic structures. The research progress of two-dimensional unidirectional self-locking structures, three-dimensional multi-directional self-locking structures and curved self-locking structures are then respectively introduced based on the classification of self-locking direction. Among them, the multi-directional self-locking structure withstands more complex loading conditions. Therefore, research progress of three representative multi-directional self-locking structures based on dumbbell-type, bone stitching, and origami design are further introduced. Finally, the research on the self-locking structure is summarized, and its future research prospects are discussed.

Novel Buffering Metamaterials with a Long Load Plateau: Design and Mechanical Characterization
YOU Yixuan, YE Wenkang, ZHANG Tianpeng, HU Lingling
2026, 40(7): 070104.   doi: 10.11858/gywlxb.20261042
HTML PDF(9)
Abstract:

Metamaterials with both reversible deformation and a long load plateau meet the demand for cyclic buffering, offering great application prospects in protective engineering. However, current metamaterials generally suffer from low material utilization, which limits their load-bearing and energy absorption performance. To address these limitations, a novel buffering metamaterial with a long load plateau is proposed in this work. Composed of bilaterally symmetric double-arc structures and vertically symmetric curved plates, the metamaterial is capable of recoverable large deformation and overall cooperative load-bearing deformation, thereby improving material utilization and optimizing structural load capacity and energy absorption performance. Experimental tests and numerical simulations were conducted to validate the long load plateau and recoverable large deformation characteristics of the metamaterial. The influences of structural geometric parameters on its mechanical behavior were also systematically analyzed. The results indicate that the long load plateau can be effectively tuned by adjusting the thickness of lateral double arcs, the thickness of intermediate curved plates, and the central transverse span. Once the intermediate curved plates are removed, the long load plateau feature disappears, and the force-displacement curve presents an approximately linear variation. Finite element simulations at equal mass confirm that the proposed metamaterial possesses better buffering performance than similar structures without a long load plateau, and the underlying buffering mechanism is clarified. The findings provide a novel design strategy for improving the performance of metamaterials with a long load plateau, and facilitate their application in protective engineering.

Inverse Design and Boundary Effects in Impact Waveform Control of Graded Foam Metals
YIN Jiangnan, YANG Qiang, GAO Jinling, LIU Xiaochuan, LIU Jiagui, LU Tianjian
2026, 40(7): 070105.   doi: 10.11858/gywlxb.20261046
HTML PDF(47)
Abstract:

To achieve the precise generation of high-amplitude impact waveforms required for aviation safety testing and related fields, this study investigates the impact waveform regulation mechanisms of graded cellular metals under different boundary conditions. Based on the conservation laws of mass and momentum, theoretical models for impact waveform generation using graded cellular metals are established for both free and elastic boundary conditions. Furthermore, an inverse design method for density gradient is proposed, which incorporates an average relative density constraint combined with Gauss-Newton iteration, enabling the reverse solution from a prescribed acceleration waveform to the material density gradient distribution. Finite element results demonstrate that the proposed method can effectively generate required waveforms—such as triangular and half-sine waves—under both boundary conditions. The study also reveals that: free boundaries are more suitable for generating high amplitude and long duration waveforms, whereas elastic boundaries can improve the realizability of low-amplitude waveforms through stiffness regulation; boundary conditions do not alter the impact duration but exert a significant influence on waveform shape; and excessive impedance mismatch between adjacent layers will intensify waveform oscillations, thereby compromising the waveform generation accuracy. The proposed inverse design strategy for density gradients exhibits favorable versatility and provides both theoretical support and a practical design tool for the development of high-amplitude impact testing technologies.

Hierarchical Energy Absorption and Dynamic Response of Bionic Thin-Walled-Foam Composite Structures Based on Mechanical Matching Design
GAO Dandan, YAN Hao, ZHOU Ying, WANG Tao, HUANG Guangyan
2026, 40(7): 070106.   doi: 10.11858/gywlxb.20261043
HTML PDF(50)
Abstract:

To achieve the synergistic improvement of load-bearing stability and energy absorption in lightweight protective structures, a bio-inspired thin-walled-foam composite structure based on mechanical matching design was proposed. Three configurations of polylactic acid (PLA) bio-inspired shells were fabricated via additive manufacturing, and subsequently filled with polyurethane foam through an in-situ foaming process. Tensile tests, quasi-static compression tests, and drop-weight impact tests were conducted to investigate the foaming-induced thermal effects on the mechanical properties of the PLA shells and the structural response of the composites. Crashworthiness was evaluated using peak crushing force (PCF), plateau force, specific energy absorption (SEA), mean crushing force (MCF), and crushing force efficiency (CFE). Results show that the temperature rise during foaming reduces the elastic modulus and strength of PLA while improving its ductility, thereby enhancing the mechanical compatibility between the shell and foam. Consequently, the composite structures exhibit significantly increased plateau force and MCF, and their collapse mode transforms from local instability to progressive stacked crushing, leading to stable hierarchical energy absorption. Dynamic impact tests further demonstrate the superior load-bearing and energy absorption performance of the composite structures under high-energy impact. The results highlight the synergistic role of geometric configuration, material matching, and thermal-mechanical coupling in regulating the energy absorption behavior, providing guidance for the design of lightweight bio-inspired protective structures.

Effect of Polyurethane Colloid Content on the Ballistic Performance of Prepreg PBO Fiber Composite Laminates
SU Zihan, LI Xiangyu, LIANG Minzu, WANG Jie, LIN Yuliang, ZHANG Yuwu
2026, 40(7): 070107.   doi: 10.11858/gywlxb.20261034
HTML PDF(33)
Abstract:

To address the issue of weak interfacial bonding and insufficient ballistic performance in poly-p-phenylene ben-zobisoxazole (PBO) fiber composites caused by the chemically inert fiber surface, the waterborne polyurethane impregnation followed by a hot-pressing process was employed to fabricate PBO composite laminates with varying resin mass fractions (18.40% and 20.45%) and numbers of fiber layers (10, 20, and 30 layers). The tensile properties, ballistic limit, failure modes, and energy absorption mechanisms of the composites were investigated through quasi-static tensile testing, ballistic experiments, and X-ray computed tomography (CT). The results indicate that this process can effectively improve the formability and protective performance of the PBO fiber composite laminates. In comparison with the specimen exhibiting a resin mass fraction of 20.45%, the 30-layer specimen with a resin mass fraction of 18.40% demonstrated a 14.86% increase in tensile strength and a 28.33% increase in elastic modulus. The ballistic limit increases with resin content. The ballistic limit of the 10-layer and 20-layer specimens with a resin mass fraction of 20.45% increases by 9.9% and 5.3%, respectively. However, the increasing effect diminishes with a greater number of layers. For the 30-layer specimens, the difference in resin mass fraction between the two resins was less than 1%. The primary failure modes of the composite laminates include fiber shear fracture, matrix cracking and delamination, and fiber tensile fracture. Energy absorption is achieved through the synergistic mechanisms of fiber compressive deformation, shear, and tensile fracture. Efficiency of the energy absorption decreases with increasing impact velocity, while it increases with the number of layers under the same impact velocity. The effect of colloid content diminishes significantly for numerous layers. The study can provide a reference for the design of ballistic protection using PBO fiber composites.

Experimental Study on the Impact Dynamics Behavior of Ultrathin Carbon Fiber Composites
ZHAO Changfang, LIU Hao, ZHOU Caihua, ZHOU Zhitan, JI Liang
2026, 40(7): 070108.   doi: 10.11858/gywlxb.20251265
HTML PDF(93)
Abstract:

Carbon fiber reinforced polymer (CFRP), as an advanced composite material, are widely used in engineering applications. However, research on the dynamic mechanical behavior of ultrathin CFRP laminates remains relatively limited. In this study, unidirectional ultrathin prepreg and hot-pressing molding processes were employed to fabricate ultrathin CFRP laminates with a single ply thickness of only 0.1 mm. The strain rate effects on specimens with five different ply orientations—0°, 90°, 0°/90°, 45°, and ±45°—were systematically investigated. Quasi-static compression experiments indicated that the 45° ply orientation enhanced plastic behavior but reduced material strength and modulus, whereas the 90° ply orientation contributed to increased modulus and strength while reducing plastic deformation. Dynamic impact tests revealed that the 90° ply orientation improved both dynamic modulus and strength while decreasing yield strain. Although the 45° ply orientation reduced dynamic yield strength, it significantly increased the sensitivity of dynamic modulus and yield strain to strain rate. Compared with conventional CFRP laminates with a 0°/90° ply layup (ply thickness is 0.295 mm, and dynamic strength and modulus are 900 MPa and 10.12 GPa, respectively), the ultrathin CFRP composites developed in this study exhibited a 66% increase in fiber content per unit thickness; under the 0°/90° ply configuration, dynamic strength and modulus were enhanced by 123% and 926%, respectively. Based on the experimental data, a constitutive model for the ultrathin CFRP composites was established, and corresponding constitutive parameters were provided, offering a basis for predicting the mechanical behavior of CFRPs under different ply orientations and strain rates.

Formation Mechanism of Ceramic Cones in Boron Carbide Ceramic Composite Targets with Different Backplates under High-Velocity Impact
WANG Xinde, LI Mingshu, WANG Renjie, WANG Yonggang, JIANG Zhaoxiu
2026, 40(7): 070109.   doi: 10.11858/gywlxb.20261037
HTML PDF(36)
Abstract:

To investigate the influence of backplate mechanical properties on the formation mechanism of ceramic cones in boron carbide ceramic composite armor, four typical backplate materials, including 6061 aluminum alloy, 7075 aluminum alloy, T300 carbon fiber board, and ultra-high molecular weight polyethylene (UHMWPE) were selected. A combination of ballistic impact experiments conducted via a one-stage light gas gun and numerical simulations performed with LS-DYNA was adopted to systematically study the effects of backplate yield strength, stiffness, and wave impedance on the morphology and evolution of ceramic cones. The results indicate that: the load transfer from the ceramic cone to the backplate is not solely dependent on a single outer cone but is achieved through the synergistic action of multiple cracks, including the outer and inner cones; the yield strength of the backplate has no significant effect on the crack propagation of the main cone; regarding stiffness, the outer cone angle decreases linearly with increasing elastic modulus, while the inner cone angle increases exponentially; wave impedance alters the internal stress field of the ceramic by modulating stress wave reflection/transmission, resulting in a linear increase in the inner cone angle and an exponential decrease in the outer cone angle with increasing impedance.

Impact Response and Design Optimization of Triangular Corrugated Sandwich Beams: A Machine Learning Approach
LI Dong, ZHANG Xiaobin, LIU Zhifang, LEI Jianyin
2026, 40(7): 070111.   doi: 10.11858/gywlxb.20251287
HTML PDF(66)
Abstract:

To enhance the prediction accuracy of low-velocity impact performance and improve the structural design efficiency of triangular corrugated sandwich beams, this paper proposes a machine learning modeling and optimization process for the impact response of sandwich beams based on a hard-parameter-sharing multi-task learning (MTL) framework. A sample dataset is generated using finite element models, and the rationality of the models is validated against existing experimental results. Subsequently, an MTL model is trained to simultaneously predict the structural specific energy absorption (SEA), the maximum deflection of the top panel, and the initial peak load. The results show that the MTL model optimized via Bayesian optimization demonstrates strong predictive performance under a 50 J impact energy condition. The predictions align well with the finite element simulation results, with the coefficient of determination R2 for all output variables in the test set exceeding 0.989, thereby validating the effectiveness and reliability of the model in response prediction and engineering optimization analysis. Parameter sensitivity analysis reveals that the core cell count and core wall thickness have the most significant influence on structural stiffness, followed by the top panel thickness, while the bottom panel thickness has a relatively minor impact. Moreover, the core wall thickness exhibits a certain saturation threshold in terms of performance enhancement. In combination with the non-dominated sorting genetic algorithm Ⅱ (NSGA-Ⅱ), multi-objective optimization analysis are conducted focusing on deformation characteristics, energy absorption performance, and comprehensive performance, and yields optimal parameter configurations that meet different engineering design requirements for sandwich beams.

Crashworthiness of Bionic Fractal Multi-Cell Circular Tubes under Axial Load
GAO Jianming, ZHANG Xiaobin, LIU Zhifang, LEI Jianyin, LI Shiqiang
2026, 40(7): 070113.   doi: 10.11858/gywlxb.20251250
HTML PDF(86)
Abstract:

A bio-inspired fractal multi-cell circular (BFMC) tube with embedded regular polygons is proposed to address the gap between the need for high absorption and the limited performance of traditional thin-walled circular tubes. Inspired by biological structures and fractal hierarchy theory, geometric models of BFMC tubes embedded with square, pentagonal, and hexagonal cells are constructed. Numerical simulations are carried out to systematically investigate the effects of mass, fractal dimension, and the number of sides of the embedded polygons on the axial crushing performance, and the results are compared with those of typical multi-tubes. The results indicate that, under approximately equal mass conditions, the BFMC tube can significantly enhance the specific energy absorption and the load-bearing capacity owing to its fractal hierarchical and bio-inspired configurations. Its crashworthiness increases with mass, first decreases and then rises as the fractal dimension increases, and improves further as the number of polygon sides increases, while the peak force is only weakly affected. A theoretical model for predicting the mean crushing force of BFMC tubes is developed based on the super folding element theory and is validated through numerical simulations. This study provides theoretical support and structural design guidelines for developing high performance thin-walled energy absorption structures.

Effect of Interlayer Materials on the Interfacial Microstructure and Dynamic Mechanical Properties of 7B53 Aluminum Alloy Composite Plates
WAN Yu, CHEN Zejun, CAO Xianming, DANG Yuehui, CONG Fuguan, WANG Qiang
2026, 40(7): 070114.   doi: 10.11858/gywlxb.20261044
HTML PDF(39)
Abstract:

The choice of aluminum alloys used as the interlayers significantly influences the interfacial bonding properties and dynamic impact mechanical properties of 7B53 aluminum alloy composite plates (7A52/interlayer/7A63). In this study, the influence mechanism of different aluminum alloy interlayer materials (7A01, 6061, 2024 aluminum alloy) on the interfacial metallurgical bonding quality and the dynamic mechanical behavior at high strain rates (17003200 s−1) was systematically investigated using tensile-shear tests, Charpy impact tests, split Hopkinson pressure bar (SHPB) tests, and scanning electron microscopy (SEM). The results show that the composite plate with the 6061 interlayer exhibits the optimal interfacial bonding performance, achieving a maximum shear strength of 109.6 MPa, which is 36.5 MPa higher than that of the plate with the 7A01 interlayer (73.1 MPa). This improvement is attributed to the fact that the 6061 alloy promotes the formation of fine and uniform grains at the interface, thereby effectively strengthening the interfacial region. SHPB tests reveal that the inhomogeneous deformation of the interlayer interrupts the penetration of cracks into the 7A52 layer and promotes crack deflection along the interface. The composite plate with the 7A01 interlayer shows low strain rate sensitivity. The plate with the 6061 interlayer, while decreases in flow stress due to thermal softening within the strain rate range of 17002700 s−1, maintains stable deformation under high-velocity impact owing to its excellent ductility. Compared to the composite plate with the 2024 interlayer, the plate with the 6061 interlayer achieves higher plastic strain while retaining relatively high yield strength. The 6061 interlayer composite plate successfully achieves an effective integration of the high toughness of 7A52 and the high strength of 7A63, providing an important theoretical basis for the design of impact-resistant protective structures for armored vehicles.

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](15417) [PDF 8643KB](2590)
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](11085) [PDF 2836KB](2324)
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](17445) [PDF 6073KB](2856)
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](18138) [PDF 1350KB](1373)
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](16386) [PDF 12118KB](1527)
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](19571) [PDF 689KB](1323)
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](10039) [PDF 1765KB](601)
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](13028) [FullText HTML](5628) [PDF 2527KB](5628)
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](12508) [PDF 402KB](1146)
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](15480) [PDF 1180KB](1304)
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.
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](18039) [PDF 365KB](1418)
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.
Perimeter-Area Relation of Fractal Island
LONG Qi-Wei
1990, 4(4): 259-262 .   doi: 10.11858/gywlxb.1990.04.004
[Abstract](17477) [PDF 1508KB](2586)
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](16042) [PDF 2450KB](1143)
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.
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](17208) [PDF 4054KB](1348)
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.
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](12635) [PDF 411KB](946)
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.
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](15205) [PDF 1599KB](997)
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.
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](16103) [PDF 794KB](1158)
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.
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](12965) [PDF 416KB](1284)
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.
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](16108) [PDF 500KB](1308)
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.
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](14758) [PDF 5649KB](1447)
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.