Shearing Process Research Articles

Dynamic shear behaviours of granite under coupled static and high-rate loadings

Understanding dynamic shear behaviours of rock under quasi-static confinement and dynamic loads is not only of great importance in underground engineering structures, but also beneficial to reveal earthquake mechanics regarding pre-stressed rocks fracturing from initial intact to rapid shear sliding. Here we propose a dynamic shear test setup and present full-field observations of stressed rock failure at high shear rate (up to γ˙ = 220 s−1) and sliding at high slip velocity over 10 m/s in the order of coseismic slip, with the integration of a novel triaxial Hopkinson bar (Tri-HB) and high-speed three-dimensional Digital Image Correlation (3D-DIC) techniques. Dynamic shear strength and friction coefficient are determined during the dynamic shear test under confinement, with different shear strain rates and pre-stresses. The dynamic shearing process can be divided into three stages, i.e. shear deformation, fracture degradation and dynamic frictional sliding. The high energy density of shear fracture band is observed, the rate and stress state dependences of rock shear strength and dynamic frictional coefficient are revealed. The dynamic shear strength of granite follows the Mohr-Coulomb criteria under varying pre-stresses and strain rates. The gouge materials with size of 1∼10 μm are generated within the high-rate shearing zone, and this appears to control the frictional resistance. The results could help determination of rock shear strength and post-shear failure, as well as explain the shear rupture and the resistance to slip mechanism of rock faults and strong near-surface earthquake events.

Study on strain localization of frozen sand based on uniaxial compression test and discrete element simulation

Strain localization has always been an important subject in frozen soil mechanics and engineering. To evaluate the development of local strain and the formation of shear bands in frozen soil, uniaxial compression tests have been conducted on frozen sand at various temperatures and particle grades. The strain localization evolution law of the frozen soil is analyzed utilizing the digital image correlation (DIC) method. The test results reveal that the entire process of shear band generation, development, and formation in frozen soil can be well captured. Within the testing temperature and particle grade range in this study, in comparison to particle grade, temperature exerts a more pronounced influence on the shear band angle which increases as temperature decreases. It is discovered that the width of the shear band increases with the decrease in temperature and the increase in mean particle diameter d50. Subsequently, a discrete element method (DEM) model is developed to examine the microscopic mechanical characteristics of frozen sand in uniaxial compression tests. The reliability of the DEM model is verified through comparative analysis with test results. Besides, the development law of the strain localization of the frozen soil obtained by the numerical simulation is revealed based on the quantitative analysis of the rotational angle, bonding state, and displacement of the soil particle during the shearing process.

Stiffness Anisotropy and Micro-Mechanism of Calcareous Sand with Different Particle Breakage Ratios Subjected to Shearing Based on DEM Simulations

Stress-induced anisotropy in calcareous sand can cause an uneven displacement in island reef engineering. In this study, stiffness, as a quantitative indicator, is explored to reveal the stress-induced anisotropy in calcareous sand. Based on the discrete element method, the stiffness anisotropic characteristics of calcareous sand during shearing, as well as the impact of particle breakage, are investigated by numerical simulations. Both the macro and micro responses, i.e., the maximum shear modulus, contact normal, strong and weak contact normal force, and the direction of particle breakage, are explored for calcareous sand with different particle breakage ratios. The results show that calcareous sand exhibits notable anisotropy during shearing, with the maximum shear modulus in the vertical direction (deviatoric stress direction) being significantly greater than that in the horizontal direction. Moreover, the higher the particle breakage rate, the lower the stiffness and its anisotropy. The micro-mechanism results indicate that the primary particle breakage during the shearing process occurs in the vertical direction. That is, the particle breakage weakens the strong contact force in the vertical direction, leading to a redistribution of the strong contact forces from the vertical direction to other directions. This redistribution mainly manifests in a decrease in the anisotropy of contact normal and contact vector within the sample, as well as a decrease in the proportion of strong contact forces in the overall contacts. This, in turn, reduces the shear strength and stiffness of calcareous sand, particularly in the vertical direction, and results in a decrease in the maximum shear modulus and its anisotropy. The maximum reduction can be up to 50% of the original value. These insights can provide a certain theoretical support for the uneven displacement and long-term stability of calcareous sand for islands and reefs.

Experimental study of the shear behavior of concrete-rock interfaces under static and dynamic loading in the context of low confinement stress

Modified static and dynamic shear tests are introduced to study the shear behavior of concrete-rock interfaces under static and dynamic loading in the context of low confinement stresses. Static and dynamic shear tests of concrete-sandstone and concrete-granite interfaces are performed using these techniques. Three levels of interface roughness are considered: smooth, bush-hammered, and rough rock surfaces. The results of these tests show that in both static and dynamic regimes, the shear evolution of concrete-rock interfaces can be described according to three successive stages: the shear stress accumulation, the shear slip, and the residual shear stress stage. The main parameters driving the shear process are the concrete-rock bonds, the interface roughness, and the residual friction. However, unlike in the static shear evolution, in the dynamic shear evolution, the concrete-rock bonds and the roughness seem active in the shear stress accumulation stage. Furthermore, the correlation between the shear strength and the normal stress is stronger in static than dynamic conditions. The significance of the normal stress on the dynamic shear strength appears more important in rough concrete-granite interfaces than in the other two interfaces. Lastly, the dynamic peak shear strengths of all the interfaces tested are three to four times higher than their static counterparts.

Continuous tow shearing for the automated manufacture of defect-free complex 3D geometry composite parts

It has always been challenging to manufacture a composite structure with complex geometry using automated fibre placement (AFP) process. The tool surface is difficult to be tessellated using fibre tapes with a finite width, without producing gaps and overlaps, and fibre steering along non-geodesic fibre paths produces various defects such as tape buckling and pull-up. In this work, a defect-free fibre-steering process on a complex surface was demonstrated by realising the continuous tow shearing (CTS) process in three dimensions. A new head control algorithm was developed, which defines the head orientations and trajectories based on the pin-jointed net model to manipulate the fibre tow using both in-plane shear and out-of-plane twisting deformations. Fibre steering process utilising this new algorithm was tested on a doubly-curved surface, using an industrial robot arm equipped with a CTS prototype head, and its layup quality and accuracy were assessed by using a three-dimensional laser profile scanner. An experimental comparison with the conventional AFP process showed that the new control algorithm enables defect-free fibre steering on complex 3D surfaces allowing for simplifying the design and analysis of novel composite structures.

Tribological properties of polyaromatic organics as lubricating grease

The tribological properties of polyaromatic organics (xPyE) as lubricating grease additives are studied by simulations and experiments. Molecular Dynamics (MD) simulations are used to simulate the mechanical properties of four kinds of polyaromatic organics (2P1E, 3P2E, 4P3E, and 5P4E) with different aromatic ring numbers and a mass fraction under 300 MPa and 5 m/s relative shear velocity of the two solid walls. The friction force, positive pressure, stress, and adsorption energy are calculated quantitatively. Meantime, the friction coefficient of lubricating grease and the wear morphology of the sample are measured by Multifunctional Friction and Wear Testing Machine. The results show that the organic additives adsorbed on the wall form a protective film during the shearing process, improving the bearing capacity and reducing friction and wear. Among them, the 3P2E organic additive has an appropriate number of aromatic rings. Its single molecular adsorption energy and unit adsorption energy are higher than other organic additives and can uniformly lay on the wall to form a stable, protective film. This demonstrates the best friction performance and improves the load-bearing capacity of the lubricating film, reducing maximum wall stress by 53.4 %.

Shear damage mechanisms of jointed rock mass: a macroscopic and mesoscopic study.

The joints are existing throughout the underground rock mass. It is of great significance to investigate the shear performance of the rock mass to maintain the stability of the underground structure. In this study, we conducted orthogonal tests to determine the proportion of rock-like materials, and used JRC curves to make specimen molds and then prepare the specimens. We conducted straight shear tests and uniaxial compression tests to determine the various mechanical parameters of the rock-like materials. Next, we carried out the compression and shear tests to investigate the shear characteristics of the specimens, and study the damage pattern and shear strength of the jointed rock mass under different confining pressures and roughness levels. The mesoscopic displacements in the shear process of joints were analyzed by using ABAQUS. The test results show that the effect of the confining pressure on the shear strength of the joint plane is relatively obvious, and a larger confining pressure indicates a larger shear strength. The effects of different joint plane roughness and shear rated on the shear characteristics of the joint plane are also significant. The mesoscopic displacement difference inside the joint plane with higher roughness is relatively large, and the stress concentration phenomenon is obvious and lasts longer, which leads to the faster destruction of the specimen with higher roughness and the higher destruction degree. Therefore, we suggest that the priority should be given to the reinforcement of jointed rock mass with high roughness during the construction to prevent sudden destabilization and failure.

Cyclic shear behavior and BoBiLSTM-based model for soil-rock mixture-concrete interfaces

In order to study the impact of surface roughness on the cyclic shear characteristics of the Soil-Rock Mixture and concrete interface, a series of cyclic shear tests were conducted using a large indoor direct shear apparatus. The effects of three concrete surface roughness coefficients JRC (0.4, 9.5, 16.7), five rock content levels (0%, 25%, 50%, 75%, 100%), and three cyclic shear displacement amplitudes (1, 3, 6 mm) on interface cyclic shear stress and Soil-Rock Mixture shear deformation were analyzed. A Bidirectional Long Short-Term Memory (BoBiLSTM) model was proposed, utilizing Bayesian optimization and k-fold cross-validation for hyperparameter tuning to streamline the model parameter selection process and enhance the prediction accuracy of the stress-strain relationship under cyclic loading. The experimental results show that, under five rock content levels, as the concrete surface roughness coefficient and cyclic shear displacement amplitude increase, the interface average peak shear stress increases accordingly. The interface average peak shear stress of the sample with 75% rock content is the highest; in terms of vertical displacement, the sample with 50% rock content has the maximum displacement, while the sample with 25% rock content has the minimum. The two types of samples show different soil deformation patterns in the two shear directions during the cyclic shearing process; as the shear displacement amplitude increases from 1 mm to 3 mm and 6 mm, the greater the concrete surface roughness, the smaller the change in shear stiffness and damping ratio. Compared to traditional Long Short-Term Memory (LSTM) models, the BoBiLSTM model demonstrated improvements in the average metrics of R2, RMSE, and MAPE by 0.32%, 57.25%, and 72.32%, respectively.

Exploring the relevance between load-bearing capacity and surface friction behavior based on a layered hydrogel cartilage prototype

Cartilage is well lubricated over a lifetime and this phenomenon is attributed to both of the surface hydration lubrication and the matrix load-bearing capacity. Lubricious hydrogels with a layered structure are designed to mimic cartilage as potential replacements. While many studies have concentrated on improving surface hydration to reduce friction, few have experimentally detected the relationship between load-bearing capacity of hydrogels and their interface friction behavior. In this work, a bilayer hydrogel, serving as a cartilage prototype consisted of a top thick hydrated polymer brush layer and a bottom hydrogel matrix with tunable modulus was designed to investigate this relationship. The coefficient of friction (COF, μ) is defined as the sum of interfacial component (μInt) and deformation/hysteresis component (μHyst). The presence of the top hydration layer effectively dissipates contact stress and reduces the interface interaction (μInt), leading to a stable and low COF. The contribution of mechanical deformation (μHyst) during the sliding shearing process to COF can be significantly reduced by increasing the local mechanical modulus, thereby enhancing the load-bearing capacity. These results show that the strategy of coupling surface hydration layer with a high load-bearing matrix can indeed enhance the lubrication performance of hydrogel cartilage prototypes, and implies a promising routine for designing robust soft matter lubrication system and friction-control devices.

Micromechanical modeling of hollow cylinder torsional shear test on sand using discrete element method

Previous studies on the hollow cylinder torsional shear test (HCTST) have mainly focused on the macroscopic behavior, while the micromechanical responses in soil specimens with shaped particles have rarely been investigated. This paper develops a numerical model of the HCTST using the discrete element method (DEM). The method of bonded spheres in a hexagonal arrangement is proposed to generate flexible boundaries that can achieve real-time adjustment of the internal and external cell pressures and capture the inhomogeneous deformation in the radial direction during shearing. Representative angular particles are selected from Toyoura sand and reproduced in this model to approximate real sand particles. The model is then validated by comparing numerical and experimental results of HCTSTs on Toyoura sand with different major principal stress directions. Next, a series of HCTSTs with different combinations of major principal stress direction (α) and intermediate principal stress ratio (b) is simulated to quantitatively characterize the sand behavior under different shear conditions. The results show that the shaped particles are horizontally distributed before shearing, and the initial anisotropic packing structure further results in different stress–strain curves in cases with different α and b values. The distribution of force chains is affected by both α and b during the shear process, together with the formation of the shear bands in different patterns. The contact normal anisotropy and contact force anisotropy show different evolution patterns when either α or b varies, resulting in the differences in the non-coaxiality and other macroscopic responses. This study improves the understanding of the macroscopic response of sand from a microscopic perspective and provides valuable insights for the constitutive modeling of sand.

Impact of rise and fall phases of shear on platelet activation and aggregation using microfluidics.

Blood flow disorders are often the result of the non-physiological narrowing of blood arteries caused by atherosclerosis and thrombus. The blood then proceeds through rising-peak-decreasing phases as it passes through the narrow area. Although abnormally high shear is known to activate platelets, the shear process that platelets undergo in small arteries is complex. Thus, understanding how each shear phase affects platelet activation can be used to improve antiplatelet therapy and decrease the risk of side effects like bleeding. Blood samples were sheared (68.8 ms,5200s-1) in vitro by the microfluidic technique, and platelet activation levels (P-selectin and integrin αIIbβ3) and von Willebrand factor (vWF) binding to platelets were analyzed by flow cytometry. Post-stenosis platelet aggregation was dynamically detected using microfluidic technology. We studied TXA2, P2Y12-ADP, and integrin αIIbβ3-fibrinogen receptor pathways by adding antiplatelet drugs, such as acetylsalicylic acid (ASA, an active ingredient of aspirin that inhibits platelet metabolism), ticagrelor (hinders platelet activation), and tirofiban (blocks integrin αIIbβ3 receptor) in vitro, respectively, to determine platelet activation function mediated by transient non-physiological high shear rates. We demonstrated that platelets can be activated under transient pathological high shear rates. The shear rise and fall phases influenced shear-induced platelet activation by regulating the binding of vWF to platelets. The degree of platelet activation and aggregation increased with multiple shear rise and fall phases. ASA did not inhibit shear-mediated platelet activation, but ticagrelor and tirofiban effectively inhibited shear-mediated platelet activation. Our data demonstrated that the shear rise and fall phases play an important role in shear-mediated platelet activation and promote platelet activation and aggregation in a vWF-dependent manner. Blocking integrin αIIbβ3 receptor and hindering P2Y12-ADP were beneficial to reducing shear-mediated platelet activation.

Predicting the magnitude of technogenic earthquakes during underground mining of the Zhezkazgan ore field

Purpose. Determining a relationship between the shear/failure of the rock mass and the technogenic earthquakes caused by them during underground mining of ore deposits with the derivation of an equation for the dependence of the technogenic earthquake indicators for the conditions of the Zhezkazgan copper-ore field. Methods. The research methodology consists of studying and analyzing early research on technogenic earthquakes by the method of statistical data processing. Calculation of rock mass deformation distribution in the study area is based on numerical modeling in the Comsol Multiphysics 5.6 and MATLAB 2020 software package environment. Calculations are performed based on solving a plane problem using the Finite Element Method (FEM). Findings. Based on the transition from the failure area to subsidence and shear values of the overlying rock mass stratum, an equation for the dependence of the earthquake magnitude on the numerical values of the mass subsidence or shear has been obtained. Originality. For the first time, based on the physics and geomechanics of rock mass shear processes, empirical-analytical formulas have been obtained that make it possible to predict the technogenic earthquake magnitude during underground mining of ore deposits in the conditions of the Zhezkazgan copper-ore field. Practical implications. Preliminary predictive calculations made by the obtained formulas for the conditions of the active mine No. 31 of the East-Zhezkazgan mine, ТОО Kazakhmys Smelting, show acceptable results of magnitude value, comparable to in-situ measurements during field mining. This prediction makes it possible to pre-calculate the technogenic earthquake magnitude at the stage of designing mining operations and make appropriate scientifically sound decisions during further mining of the field.

HARDWARE AND SOFTWARE COMPLEX FOR STABILITY CONTROL OF QUARRY SIDES USING OPTICAL FIBER

This article explores solutions to address the issue of implementing a shift monitoring system in quarries, which is the primary cause of collapsing sides and edges. For guaranteed safety and continuous monitoring of the quarry's activities, advanced fiber-optic sensors have been integrated. Their design is based on the use of single-mode optical fibers, which provides the most accurate assessment of displacements in mountainous terrain up to a distance of 30 km. Extensive laboratory studies have been carried out to deeply analyze the deformation and shear processes occurring in this geographical area. The results obtained from these analyses bode well for a significant improvement in quarry performance and a reduction in the likelihood of emergency situations.

The Thermophysical Aspects of the Transformation of Porous Structures in Versatile Nanostructured Materials

The technology of obtaining porous nanostructures is based on ecological organosilicon materials and their uses in some spheres of human life, for example, for medical preparations, for thermal insulation of building structures and industrial equipment, and for cleaning. The purpose of this study was to establish correlations between various experimental parameters (shear stress, speed pulsations, temperature, viscosity, and processing time) and the rheological characteristics of suspensions obtained by the method of liquid-phase dispersion; it was a study of hydrodynamic effects and the processes of heat and mass exchange in liquid systems during the liquid-phase dispersion of hydrogel monoliths by means of discrete-pulse activation in a special rotary apparatus. The dehydration of hydrogels was carried out by two methods: convective drying in a layer and spraying in the coolant flow. Experiments have shown that the key parameters for obtaining stable homogeneous suspensions are a synergistic combination of concentration factors and processing time. To obtain adsorbents in the form of pastes with specified adsorption properties and a monolith size of up to 300 μm, the optimal parameters were a hydrogel concentration of 70% and a processing time in the double-recirculation mode. Xerogels obtained by convective drying are a polydisperse mixture of strong monoliths and fragile aggregates. In contrast, xerogel monoliths obtained by spray drying show great homogeneity in terms of dispersion and strength characteristics. The rheological parameters of the hydrogel dispersions, which depend on the concentration and hydrodynamic treatment modes, are the dominant factors affecting the moisture extraction during drying. This study marks the first investigation into the resilience of porous organosilicon structures against the influence of intense turbulence fields and mechanical stresses experienced within the rotor apparatus during suspension production.

Stress-free single-cell transcriptomic profiling and functional genomics of murine eosinophils.

Eosinophils are a class of granulocytes with pleiotropic functions in homeostasis and various human diseases. Nevertheless, they are absent from conventional single-cell RNA sequencing atlases owing to technical difficulties preventing their transcriptomic interrogation. Consequently, eosinophil heterogeneity and the gene regulatory networks underpinning their diverse functions remain poorly understood. We have developed a stress-free protocol for single-cell RNA capture from murine tissue-resident eosinophils, which revealed distinct intestinal subsets and their roles in colitis. Here we describe in detail how to enrich eosinophils from multiple tissues of residence and how to capture high-quality single-cell transcriptomes by preventing transcript degradation. By combining magnetic eosinophil enrichment with microwell-based single-cell RNA capture (BD Rhapsody), our approach minimizes shear stress and processing time. Moreover, we report how to perform genome-wide CRISPR pooled genetic screening in ex vivo-conditioned bone marrow-derived eosinophils to functionally probe pathways required for their differentiation and intestinal maturation. These protocols can be performed by any researcher with basic skills in molecular biology and flow cytometry, and can be adapted to investigate other granulocytes, such as neutrophils and mast cells, thereby offering potential insights into their roles in both homeostasis and disease pathogenesis. Single-cell transcriptomics of eosinophils can be performed in 2-3 d, while functional genomics assays may require up to 1 month.

Failure mechanism of hot dry rock under the coupling effect of thermal cycling and direct shear loading path

Inducing shear failure in reservoir rocks is a critical strategy for enhancing hot dry rock resource exploitation, as it promotes the development of an extensive fracture network. In this study, granite was subjected to direct shear tests under three distinct temperatures and five thermal cycling treatments. Utilizing acoustic emission (AE) and digital image correlation (DIC) techniques, the shear failure traits and the evolutionary patterns of cracks in granite under diverse thermal cycling treatments were probed. Experimental results reveal that the direct shear strength of granite initially decreases, subsequently rebounds with increasing thermal cycles, and exhibits a steady decrease with temperature increase. Effects of thermal shock provoke a surge in AE counts and energy during the initial half of the shearing process, steering the failure mode of the samples from brittle to plastic. Furthermore, the macro crack morphology changes from a singular shear crack to a pair of parallel shear cracks, becoming intricately complex at elevated temperatures (600 °C). Additionally, the initiation of the main shear crack in granite consistently manifests a mixed tensile-shear failure, with the proportion of shear action displaying a positive correlation with temperature. In contrast, the proximal secondary cracks are purely tension-initiated. Finally, a competitive mechanism between rock strengthening and weakening induced by thermal cycles treatment was proposed to explain the remarkable rebound in direct shear strength following thermal cycling treatment. The current work provides insights for the mechanism of shear failure in rocks under different thermal cycling conditions in geothermal reservoir stimulation.

Spectral energy transfer analysis of a forced homogeneous isotropic turbulence using triple decomposition of velocity gradient tensor

Triple decomposition of the velocity gradient tensor is exploited to investigate the role of the simple-shear, normal-strain, and rigid-body rotation, in the turbulent kinetic energy transfer of a forced, high-Reynolds-number homogeneous isotropic turbulent flow. Splitting the total energy flux into three elementary partial fluxes shows that the rigid-body rotation has a secondary contribution to turbulent kinetic energy transfer, while within 80% of the total energy flux is produced through the simple-shear and normal-strain processes. Energy is dominantly extracted through the shearing process from larger motion scales and injected into smaller scales by the normal-straining process. In addition, each partial energy flux contains a strong conservative part and a weak non-conservative part. Besides the total kinetic energy transfer, all partial energy fluxes also show scale-invariance behaviour over the inertial subrange. Statistically, total energy flux shows a strong correlation with the simple-shear and normal-strain processes. However, the relatively poor correlation between these two dominant processes, implies a fairly low efficiency for turbulent kinetic energy transfer. These results can provide new insights into sub-grid scale turbulence modelling based on the physical mechanism of simple-shear and normal-strain processes.

Turbulence calculation based on the extended Navier–Stokes equations

In this study, phenomenological observations and the Kreuer interpretation of the origin of viscosity were used to develop a computational method for solving the turbulence problem of incompressible viscous Newtonian fluids based on extended Navier–Stokes (N–S) equations. The shear process in fluid flow was hypothesized to be accompanied by eddy formation, and the effects of eddies on the convection and diffusion were considered. The classical N–S equations were improved to obtain extended N–S equations. The extended equations are closed, and the sources of the velocity fluctuations are explicitly considered to be additional convection and diffusion. The extended equations are compatible with the classical N–S equations; thus, they can describe laminar and turbulent flows in a unified manner. In fluid flow simulations, the equations describing the mean flow quantities could be directly obtained from the extended N–S equations without any additional turbulence models. A numerical investigation was carried out to verify the extended equations by exploring the flow over a cube placed in a channel. The simulation results were compared with both the large eddy simulation and experimental results.

Numerical investigation of fluid flowing through rough fractures subject to shear

Comprehending fluid flow in rock masses is essential for modern underground engineering, including chemical energy extraction, nuclear pollutant remediation, and hydrocarbon utilization, complicated by shear-induced and surface roughness effects in fractures. This study employed numerical simulations to investigate the fluid flow behavior in fractures with different surface roughness under shear, where the shear direction is perpendicular to the flow direction. The nonlinear flow of the fluid is observed to have a strong correlation with the confining pressure (Pz), roughness (JRC), and shear displacement (u). The generation of eddy currents is frequently linked to the presence of flow channel intricacies and the velocity of flow at a microscopic scale. The Forchheimer equation could describe the process of nonlinear phenomena accentuation very well. The fracture under Pz caused a reduction in hydraulic transmissivity (T) due to compression. Furthermore, the T changes dramatically as the shear process progresses. Based on the analysis of the Forchheimer coefficient (β) and critical Reynolds number (Rec) parameters used to determine the response of nonlinear flow, it appears that an increase in Pz facilitates the transition of the fluid into a nonlinear flow state. Conversely, shearing has the opposite effect and reduces the tendency toward nonlinear flow.

Shear Resistance Evolution of Geogrid-Aggregate Interfaces under Direct Shear: Insights from 3D DEM Simulations

This paper presents a mesoscopic evaluation of the shear resistance evolution of geogrid-aggregate interfaces subjected to direct shear loading. A three-dimensional discrete element method (DEM) model was developed based on experimental data. The tensile response of geogrid were simulated through a series of calibration tests. Aggregate with complex particle shapes were simulated to accurately capture the interlocking effect among aggregates based on the real particle surface. The individual shear resistance components were quantified based on particle displacement field and contact distribution characteristics. The influences of aperture-aggregate size ratio and geogrid stiffness on the shear resistance components are discussed. The results indicate that the peak value of shear resistance component follows a descending order from frictional resistance of aggregate, to passive resistance of transverse rib, and to geogrid-aggregate interface frictional resistance. During the shear process, the frictional resistance of aggregate becomes active first, followed by the geogrid-aggregate interface frictional resistance, and then the development of passive resistance of transverse ribs starts with a certain lag. Optimizing the geogrid-aggregate size ratio and utilizing geogrids with higher rib stiffness could enhance the passive resistance of transverse ribs but would not significantly affect the geogrid-aggregate interface frictional resistance and frictional resistance of aggregate.