Contact Force Chains Research Articles

3D DEM investigation on macro-meso mechanical responses of gas hydrate-bearing sediments in process of unloading confining pressure

Well drilling, gas production, and reservoir stimulation can disrupt stress states of gas hydrate-bearing sediments (GHBSs) surrounding wellbores, potentially inducing wellbore instability due to confining pressure unloading. Therefore, it is very important to know about the macro-meso mechanical responses of GHBSs under an unloading confining pressure path (UCPP). Herein, a discrete element method (DEM) was used to perform triaxial tests on three-dimensional numerical GHBSs under the UCPP and a conventional drained path (CDP) to bridge the knowledge gap. Results show that GHBSs under the UCPP are more prone to failure and dilatation than under the CDP, due to the reduced bearing capacity of contact force chains, relatively unstable particle assembly, and rapid cementation breakage development. Additionally, the stress path has little effect on the cohesion and friction angle of GHBSs. Furthermore, under the UCPP, the failure and dilatation risk of GHBSs are positively correlated with unloading confining pressure rate, inversely related to hydrate saturation and initial effective confining pressure. The grey relational analysis illustrates that compared to hydrate saturation and initial effective confining pressure, the unloading confining pressure rate minimally affects the unloading confining pressure effect. These findings can serve for experimental studies and field applications (e.g., optimizing drilling mud density, sand production rate, and fracturing fluid flowback rate) concerning the unloading confining pressure responses of GHBSs.

Does double pre-notches have a greater impact than single pre-notch on the mechanical and fracture behavior of rock? Insights from three-point bending tests using the numerical approach of grain-based model

The mechanical performance of rock masses is significantly impacted by the existence of cracks. Distribution of the cracks may result in different mechanical or fracture behaviors of rock. To gain a preliminary understanding of such problems, a series of numerical tests of three-point bending (TPB) are designed and carried out with the consideration of single pre-notch and double pre-notches. The numerical method of grain-based model (GBM) approach is employed in this paper, with a validation implemented by comparing the results with previous experiments to demonstrate the reliability of this model. Afterward, a comprehensive investigation, including the mechanical behavior, strength characteristics, crack evolution and progressive failure mechanism of the above mentioned TPB tests are conducted. The conclusions are as follows: (1) During TPB tests, as the offset distance increases, the fracture mode of the single pre-notched specimens gradually transitions from Mode I to Mode II, with the peak fracture toughness occurring at a mixed Mode I-II fracture condition (2p/s = 0.4). Moreover, the failure angle, notch open degree and expansion area all tend to rise regardless of specimens of single or double pre-notches. (2) A thorough exploration of crack evolution shows that the specimens primarily fracture during the crack propagation and failure stages. Shear failure predominates among these stages, with microcracks mostly occurring inside the rock mineral grains. (3) When offset distance increases to a certain degree (2p/s = 0.8), the fracture of specimen ceases to be governed by the pre-notch; instead, it occurs directly from the middle of the beam specimen. This phenomenon can be attributed to the reduction of stress concentration at the notch tip, as inferred from analysis of the contact force chain. (4) In double pre-notched specimen, microcracks initiate at both pre-notches, but the final failure results from crack propagation at one notch tip. Compared to single pre-notched specimen, the presence of double pre-notches exhibits superior fracture resistance performance.

Penetration mechanism of open-ended pipe piles under jacking and driving methods

Different pile penetration methods will cause different pile penetration results. Therefore, it is very necessary to study the pile penetration mechanism under different pile penetration methods. In this paper, numerical simulations using the discrete element method (DEM) are carried out to investigate the meso-scopic mechanism of pipe pile penetration under jacking and driving methods. The mechanism is comprehensively compared and analyzed in terms of particle movement, changes in the force chain, displacement changes, resistance of penetration pile, soil lateral stress and soil vertical stress. Results show that the height of the soil plug in the jacked pile is lower than that in the driven pile. The degree of contact force chain at the bottom of the soil plug in the driven pile is slightly higher than that at the bottom of the soil plug in the pile of the jacked pile. The total resistance of the driven pile greater than that of the jacked pile. Within a certain distance below the pile tip, the jacked pile affects the vertical stresses of the soil in this area to a greater extent than the driven pile. Then, parametric sensitivity analysis is performed, which captures the effects of the penetration velocity in the jacking method and the hammering frequency in the driving method. This study can provide a reference for the selection of the penetration method of pipe piles in actual projects.

Inner structures of rapid free-surface granular avalanche over a small bump obstacle: Expansion fan, oblique shock wave, and contact anisotropy

This study investigates the inner flow characteristics of a rapid granular avalanche passing over a small bump obstacle fixed on an inclined chute using the discrete element method. Both the cross-sectional mean flow properties, such as free-surface height, mean flow velocity, and mean stresses, and the inner local flow properties, including granular temperature, coordination number, pressure, contact force orientation, and granular fabrics, were comprehensively investigated. Upstream of the obstacle, a wide compression region where mean stresses strengthen and exhibit anisotropy was observed. Employing the kinetic theory of granular gas, we revealed a smooth supersonic-to-subsonic transition near the obstacle, a phenomenon distinct from typical gas dynamics. These upstream flow phenomena are attributed to the generation of stream-wise-oriented contact force chains as the flow impacts the obstacle. Downstream of the obstacle, a complex non-monotonic expansion–compression–expansion process was observed. We demonstrated that this non-monotonic flow process reflects an inner gasdynamic-like phenomenon characterized by an expansion fan followed by an oblique shock wave. Moreover, the force chains and the inner shock structure were found to significantly influence the evolution of stream-wise velocity profiles. These findings underscore the significance of inner flow structures in shaping the dynamics of granular avalanche flow interacting with obstacles.

DEM simulation of effects of particle size distribution on meso-mechanical properties of slip zone soil

Particle Size Distribution (PSD) exerts a substantial influence on the mechanical properties of geological materials such as rocks and soils, which can be viewed at a microscale as an assembly of discrete particles. An exploration into the effects of particle gradation on the properties of these materials provides valuable insights into their nature. In the study, the Discrete Element Method (DEM) was used to conduct numerical shear tests on eight distinct groups of slip zone soil, each characterized by a different particle gradation. The aim was to examine the meso-mechanical properties and shear evolution laws of slip zone soil numerical samples with both optimal and sub-optimal PSDs. Findings underscore the pivotal role that PSD plays in various aspects, including dilatancy, the evolution of the displacement field, the network of contact force chains, the principal stress, and the distribution of normal and tangential contact forces within the slip zone soil. It was observed that the network of contact force chains in the numerical samples with an optimal PSD was more complex than in those samples with a sub-optimal PSD. Additionally, the distribution of principal stresses before and after shear was more uniformly balanced. This particle size-based study offers significant reference value for future investigations into the impact of PSD on the macroscopic and meso-mechanical properties of slip zone soil. By augmenting this knowledge, a more comprehensive understanding of the fundamental behavior of these materials can be attained, leading to improved prediction and management of geological risks.

Micromechanical analysis of contact erosion under cyclic loads using the coupled CFD‒DEM method

Different layers of soil often have distinct particle sizes. When exposed to the natural environment, soil is easily affected by natural rainfall, rising groundwater levels, and human activities, leading to particle contact erosion, which reduces the safety and service performance of the soil structure. In this paper, a coupled computational fluid dynamics–discrete element method (CFD–DEM) model was employed to investigate the particle migration phenomena, mechanical response of contact interfaces, variations in flow fields, and macroscopic deformation during the contact erosion process under cyclic loads at different frequencies and amplitudes. The conclusions are presented as follows: (1) Within one cycle of cyclic loading, both compression during loading and stress relaxation during unloading are the main factors triggering the migration of fine particles. (2) The migration and loss of fine particles mainly occur in the early stages of cyclic loading, where strong contact force chains are formed within the fine particle layer, leading to significant plastic deformation of the soil at the macroscopic level. (3) Under cyclic loading, changes in the soil pore structure cause an upwards hydraulic gradient in the initial quiescent water flow field. This hydraulic gradient can rupture weak contact force chains and cause particle pumping. (4) Increasing the frequency and amplitude of cyclic loading intensifies the erosion of fine particles, causing greater axial deformation of the soil. Compared to cyclic loading frequency, the amplitude of cyclic loading has a greater impact on contact erosion.

Observed Characterization of Multi‑level Retaining Structure for Deep Excavation of Subway Station

Traditional support structures cannot meet the complex conditions of different excavation depths and areas in underground transportation hubs. On the basis of fully considering the spatial position relationship of foundation pit groups, this article proposes a multilevel retaining system that meets the requirements of multilevel foundation pit excavation. The evolution law of the support structure during the excavation process of the inner pit was explored using on-site monitoring and numerical simulation methods. The results indicate that the excavation of the inner pit reduces the passive earth pressure, and the deformation of the outer support structure can be effectively suppressed by setting a retaining structure or a bottom slab in the bench zone. The excavation of the inner pit causes significant vertical deformation of the support structure adjacent to the foundation pit, while the impact on the structure far away from the foundation pit is relatively small. According to the contact force chain and soil pressure between the two rows of support structure, the soil in this area is divided into a “relaxation zone” and a “compression zone.” The evolution mechanism of earth pressure in the case of mutual-effect failure between two rows of piles is revealed. This paper addresses the deformation properties of multilevel support structures as well as the mechanism of earth pressure evolution between structures.

Evaluation of large axial strain of granular sand in post-liquefaction stage using DEM simulation

Field investigations in recent studies have revealed that significant damage to foundations often occurs in the post-liquefaction stage, during which sand undergoes substantial deformation. As a result, numerous researchers have dedicated considerable effort to studying the evolution of the axial strain at this stage. However, the accurate measurement of the exact value of this strain is challenging mainly owing to the limitations of the triaxial apparatus. In this study, the numerical method of discrete element method (DEM) is used to simulate undrained cyclic triaxial loading tests according to a large number of loading cycles. The mechanical behaviour of granular sand under different cell pressures was reproduced using the DEM method, wherein the liquefaction resistance increased with the cell pressure. The post-liquefaction axial strain exhibited a similar pattern of variation among the calculated cases at different cell pressures. The value of the double amplitude of the axial strain stabilised at a maximum value after a certain loading cycle. The microscopic characteristics derived from the DEM calculations, including the contact number and force chain network, were analysed to provide reasonable explanations. The sand fabric was re-established after a large number of loading cycles and controlled by cyclic loading, which helped stabilise the accumulated axial strain.

Multi-Scale Research on the Mechanisms of Soil Arching Development and Degradation in Granular Materials with Different Relative Density

Soil arching is significantly influenced by relative density, while its mechanisms have barely been analyzed. A series of DEM numerical simulations of the classical trapdoor test were carried out to investigate the multi-scale mechanisms of arching development and degradation in granular materials with different relative density. For analysis, the granular assembly was divided into three zones according to the particle vertical displacement normalized by the trapdoor displacement δ. The results show that before the maximum arching state (corresponding to the minimum arching ratio), contact forces between particles in a specific zone (where the vertical displacement of particles is larger than 0.1δ but less than 0.9δ) increase rapidly and robust arched force chains with large particle contact forces are generated. The variation in contact forces and force chains becomes more obvious as the sample porosity decreases. As a result, soil arching generated in a denser particle assembly is stronger, and the minimum value of the arching ratio is increased with the sample porosity. After the maximum arching state, the force chains in this zone are degenerated gradually, leading to a decrease in particle contact forces in microscale and an increase in the arching ratio in macroscale. The recovery of the arching ratio after the minimum value is also more significant in simulations with a larger relative density, as the degeneration of contact force chains is more obvious in denser samples. These results indicate the importance of contact force chain stabilities in specific zones for improving soil arching in engineering practice.

Effect of Confining Pressure on the Macro- and Microscopic Mechanisms of Diorite under Triaxial Unloading Conditions

In this study, the response mechanism between macro- and microscales of deep hard-rock diorite is investigated under loading and unloading conditions. Moreover, the statistical theory is combined with particle flow code simulations to establish a correlation between unloading rates observed in laboratory experiments and numerical simulations. Subsequent numerical tests under varying confining pressures are conducted to examine the macroscopic mechanical properties and the evolution of particle velocity, displacement, contact force chain failures, and microcracks in both axial and radial directions of the numerical rock samples during the loading and unloading phases. The findings indicate that the confining pressure strength curve displays an instantaneous fluctuation response during unloading, which intensifies with higher initial confining pressures. This suggests that rock sample damage progresses in multiple stages of expansion and penetration. The study also reveals that with increased initial confining pressure, there is a decrease in particle velocity along the unloading direction and an increase in particle displacement and the number of contact force chain failures, indicating more severe radial expansion of the rock sample. Furthermore, microcracks predominantly accumulate near the unloading surface, and their total number escalates with rising confining pressure, suggesting that higher confining pressures promote the development and expansion of internal microcracks.

Thermal contraction coordination behavior between unbound aggregate layer and asphalt mixture overlay based on the finite difference and discrete element coupling method

AbstractThe constraint action of the unbound aggregate layer underneath plays an important role in affecting the temperature strains in the top asphalt layer. The focus of the present paper is to investigate the interactive thermal contraction mechanisms between the asphalt mixture and granular base layers to offer a new perspective in promoting the understanding of the thermal cracking disease. In this paper, both experiments and simulations were conducted for the thermal contraction tests. A novel numerical modeling of virtual composite structures was proposed using wall‐zone coupling operation based on the finite difference and discrete element coupling method. The experimental composite structures containing three types of asphalt mixture overlays and five types of unbound aggregate base layers were developed for verification. The thermal contraction and restraint mechanism were revealed from both macro and microscopic scales. The proposed numerical modeling shows a 94% high accuracy. The particle displacement vector and contact force chain could explain the thermal contraction behavior of composite structures. Smaller particle motion displacement or stronger contact force chain result in a higher restraint strain of the asphalt overlay. The thermal contraction behavior can be coordinated through the compaction and loosening of unbound aggregates. The material parameters and cooling temperature differences in the asphalt overlay have slight effects on the constraint action of a base layer, while the gradation and mechanical parameters of the unbound aggregate layer show significant impaction. The parameters, cohesion C and friction angle φ, show a quadratic function with the restraint coefficient. This work has significant guidance on the selection of pavement structure and materials to improve the thermal cracking problem and lay the basis for the mechanical theoretical calculation to predict thermal cracking.

A numerical investigation on the effect of rotation on the cone penetration test

The cone penetration test (CPT) is one of the most popular in situ soil characterization tools. However, the test is often difficult to conduct in soils with high penetration resistance. To resolve the problem, a rotary CPT device has recently been adopted in practice by rotating the rod to increase the penetrability, particularly in deep dense sand. This study investigates the underlying mechanism of the rotation effects from a micromechanical perspective using models based on the discrete element method. With rotation, the cone penetration resistance ( qc) decreases by up to 50%, while the cone torque resistance ( tc) increases gradually. These results are also used to successfully assess existing theoretical solutions. The mechanical work required during penetration is observed to keep rising as the rotational velocity increases. Microscopic variables including particle displacement and velocity field show that rotation reduces the volume of disturbed soil during penetration and drives particles to rotate horizontally, while contact force chain and contact fabric indicate that rotation increases the number of radial and tangential contacts and the corresponding contact forces, forming a lateral stable structure around the shaft, which can reduce the force transmitted to the particles below the cone, thus decreasing the vertical penetration resistance.

Discrete Element Simulation of Macro and Micro Mechanical Properties of Round Gravel Material under Triaxial Stress

In this study, the effects of relative density and confining pressure on the shear characteristics of round gravel are investigated using a large-scale triaxial apparatus and the discrete element method. A simple and efficient numerical method for simulating flexible membranes is introduced. The results show that the stress-strain curves develop from hardened to softened type with increasing relative density, while the stress-strain curves develop from softened to hardened type with increasing confining pressure. As the axial strain increases, the strong contact force chains are vertically distributed, and the larger the relative density and confining pressure, the greater the number and thickness of the strong contact force chains. In the shear process, the distribution of average normal and tangential contact forces show “peanut-shaped” and “petal-shaped”, respectively. The increase in relative density increases the anisotropy of the specimen, while the increase in confining pressure results in a decrease. A linear relationship exists between the macroscopic stress ratio and the anisotropy coefficient. The anisotropy coefficient of the normal contact force provides the greatest contribution to the macroscopic shear strength (about 55%), followed by the anisotropy coefficient of the contact normal (about 26%) and that of the tangential contact force (about 19%).

Three dimensional discrete element modelling of the conventional compression behavior of gas hydrate bearing coal

To analyze the relationship between macro and meso parameters of the gas hydrate bearing coal (GHBC) and to calibrate the meso-parameters, the numerical tests were conducted to simulate the laboratory triaxial compression tests by PFC3D, with the parallel bond model employed as the particle contact constitutive model. First, twenty simulation tests were conducted to quantify the relationship between the macro–meso parameters. Then, nine orthogonal simulation tests were performed using four meso-mechanical parameters in a three-level to evaluate the sensitivity of the meso-mechanical parameters. Furthermore, the calibration method of the meso-parameters were then proposed. Finally, the contact force chain, the contact force and the contact number were examined to investigate the saturation effect on the meso-mechanical behavior of GHBC. The results show that: (1) The elastic modulus linearly increases with the bonding stiffness ratio and the friction coefficient while exponentially increasing with the normal bonding strength and the bonding radius coefficient. The failure strength increases exponentially with the increase of the friction coefficient, the normal bonding strength and the bonding radius coefficient, and remains constant with the increase of bond stiffness ratio; (2) The friction coefficient and the bond radius coefficient are most sensitive to the elastic modulus and the failure strength; (3) The number of the force chains, the contact force, and the bond strength between particles will increase with the increase of the hydrate saturation, which leads to the larger failure strength.

Discrete modelling of the mechanical response of Cuxhaven sand under shear and oedometric conditions using the rolling resistance contact model

The Distinct element method (DEM) is a promising approach to model the microscopic behaviour of granular materials, but the capability of the simulations to reproduce the mechanical response of real granular materials depends strongly on the contact model parameters utilized. The present study focuses on the calibration and validation of the rolling resistance linear model parameters for Cuxhaven sand based on the experimental results from triaxial and oedometer tests. A sensitivity analysis was conducted to explore the influence of contact model parameters on the sample preparation and the shear stage of triaxial tests. The influence of parameters like rolling resistance friction coefficient, inter-particle friction coefficient, effective modulus, and normal-to-shear stiffness ratio were investigated. For calibration of the DEM model, the input parameters were selected such that the simulations reproduce the macro mechanical characteristics like dilation angle, peak stress, and stiffness. The calibrated parameters were then validated by simulating a one-dimensional compression test. The results are in good agreement with the experiments, which proves the suitability of the calibrated parameters. In addition, the validated parameters were applied to investigate the mechanical behaviour including the evolution of contact force chains, non-coaxiality of principal stress and strain rate, and sample inhomogeneity during a simple shear test.

Experimental and numerical studies on mechanical behavior of pre-stressed infilled rock joints subjected to horizontal vibration loading

Underground jointed rock masses are usually subjected to a coupled effect of static pre-stress and vibration disturbance from rock-cutting in Tunnel Boring Machines (TBMs) tunneling. To investigate the mechanical behavior of pre-stressed infilled rock joints disturbed by horizontal vibration loading, a series of laboratory and numerical experiments were conducted, based on the in-situ vibration characteristics during TBM rock-cutting. The P-wave velocity and direct shear test techniques are combined to quantitatively characterize the damage of infilled rock joints. The experimental results reveal that the damage degree of rock joints exhibits a nonlinear increase with the increase of vibration frequency and amplitude, and the decrease of static pre-stress. Low static pre-stress (below 0.5 MPa), high vibration frequency (above 70 Hz) and high vibration amplitude (above 0.3 mm/s) are adverse factors to the stability of rock joints. From the PFC2D simulation, the contact force chain and contact force field were analyzed to elucidate the failure mechanism of infilled rock joints. The failure mode of infilled rock joints is closely related to the level of parameters, consisting with the results in laboratory experiments. In addition, the comparison of vibration results for weakly and strongly filling material demonstrates that grouting reinforcement enhanced the vibration resistance of in-situ jointed rock masses.

Dynamic simulation of powder spreading processes toward the fabrication of metal-matrix diamond composites in selective laser melting

A novel method for the integrated production of structural and functional metal-matrix diamond composites is provided by selective laser melting (SLM) in the field of diamond tools. However, the uneven distribution of diamond particles and loose packing density in the powder bed can easily result in poor performance of diamond tools. In this work, a discrete element method model of diamond/CuSn composite powders is developed to simulate blade-spreading processes. It investigates how powder bed quality and particle dynamics are affected by diamond powder physical properties and powder spreading process parameters. The findings reveal that coarse diamond particles exhibit strong force chains at low velocities, whereas contact force chains are weak for fine CuSn particles at high velocities. It shows that diamond particles with irregular shapes have poor diffusivity and spreadability due to the large friction and mechanical locking forces, which causes diamond particles segregation. Therefore, the packing density and uniformity of the powder bed are both improved by reducing the volume fraction and particle size of diamond particles. In addition, the powder bed quality is improved by increasing the powder layer thickness and decreasing the translational speed of the blade, which weakens the shear expansion and jamming of diamond particles. The results have some significance for optimizing powder spreading processes toward the production of metal-matrix diamond composites in SLM.

Particle gradations optimization for powder spreading in additive manufacturing

Optimizing particle gradations to improve the powder bed quality is of practical engineering interest for powder spreading in additive manufacturing (AM). The discrete element model of ceramic powder is introduced to simulate blade-spreading and roller-spreading processes. Based on the packing theory, the effect of particle gradations on the powder bed quality and particle microscopic behavior is analyzed. The results show that fine particle gradations weaken the wall effect. The powder shear dilation in blade-spreading strengthens the loosening effect, however, the compaction effect of roller-spreading weakens the loosening effect. As coarse particle gradations increase, contact force chains become loose, however, strong force arches lead to particle jamming, uneven distribution, and voids in the powder bed. Particle segregation occurs because fine particle gradations pass through the gap between the substrate and spreader, and are deposited at the front part of the powder bed, while coarse particle gradations due to jamming are deposited at the end part of the powder bed. When the particle gradation coefficient is in the range of 0.1 to 0.5, the powder bed quality and particle segregation phenomenon are significantly improved compared to Gaussian distribution. The results can provide valuable references for the selection of particle gradations in AM.

Research on the damage characteristics of rock masses based on double guide-hole blasting under high in-situ stress

In complex high in-situ stress conditions, how to achieve the ideal rock blasting effect through effective methods is often a difficult point of blasting operations. This paper analyzes the influence of guide holes on blasting effect by adding guide holes to the rock pre-treatment method. Based on the particle expansion method to carry out double-hole blasting experiments, the influence of the blasthole spacing and ground stress on the blasting effect is investigated from the levels of macroscopic cracking effect, microscopic particle contact and so on. The study shows that: (i) the setting of empty holes between the gun holes can enhance the crack penetration effect, and the penetration effect is more obvious when the distance between the gun holes and the empty holes is less than 2.5 times the radius of the crushed zone. (ii) At the level of contact force chain, when the distance between blastholes and empty holes is less than 2.5 times the radius of the crushing zone. The compressive stress in the direction perpendicular to the direction of the blasthole line inhibits crack development, and the tensile stress in the direction parallel to the blasthole line promotes crack development. The main stress direction is perpendicular to the direction of the blasthole line. (iii) As the distance between the blastholes increases, the effect of crack suppression by stresses in the vertical direction decreases, and the main force direction is parallel to the direction of the blasthole line.

Macro- and micro- deterioration mechanism of high-speed railway graded gravel filler during vibratory compaction

This study aims to explore the impact mechanism of the vibratory compaction deterioration in high-speed railway graded gravel (HRGG) fillers, which can contribute to the control compaction quality and enhance the service performance of subgrade. Firstly, vibratory compaction experiments were conducted with HRGG fillers using the intelligent vibratory compactor to reveal the vibratory compaction deterioration from the evolution of the dry density ρd, dynamic stiffness Krd. Secondly, X-ray computed tomography (X-CT) scan tests were conducted with the fillers of different compaction stages. Then, the evolution of coarse particle shape characteristics was obtained to reveal the key controlling factor of compaction deterioration. Finally, the high-precision 3D vibratory compaction discrete element method (DEM) models were established for different deterioration degrees. The micro-indicators, such as the coordination number, contact force chains, and fabric anisotropy were investigated to explore the inherent relationship between the compaction deterioration and the key controlling factor. The results showed that the Krd gradually decreased at the compaction deterioration. Additionally, it was found that the key controlling factor for compaction deterioration was the abrasion of coarse particles. From the DEM simulations with different compaction deterioration degrees (CDDS), the dynamic stiffness Krd decreased with increasing CDDS, which was consistent with the results of the vibratory compaction experiments. Moreover, it was observed that the micro-indicators exhibited a decreasing trend with increasing CDDS, indicating a decrease in the ability of coarse particles to wrap fine particles and particles interlocking, which disrupts the ability of the particle skeleton to withstand external loads. This study not only provides a novel approach to investigate the deterioration mechanism during vibratory compaction but also establishes a new theoretical basis for the controlling compaction quality of HRGG fillers.