The mechanism of shield-soil interaction has always been a significant issue in academia and industry. For the active articulated shield, the presence of the active articulation system has an inevitable impact on shield-soil interaction. Therefore, a shield-soil interaction model considering active articulation was proposed and numerically solved using the time-incremental method. The model was validated through a case study of the Binzhong interval in Fuzhou Metro Binhai Express. The influence of active articulation on shield heading, resultant moment of earth pressure on shield shell, and resultant propulsion moment was carefully studied. Key conclusions include: (1) The shield-soil interaction during continuous excavation is more accurately reflected by the model and its numerical solution method. (2) Increasing the pitch articulation angle significantly reduces shield heading. (3) With a smaller coefficient of subgrade reaction, the articulation angle is approximately linearly correlated with the resultant moment of earth pressure; this relationship transitions to nonlinearity as the coefficient increases. In upper-soft lower-hard strata, the resultant moment varies with pitch articulation direction and becomes more pronounced with larger articulation angles. (4) Within the small-angle attitude correction range, a certain articulation angle reduces the resultant propulsion moment, enabling efficient attitude control. These findings provide theoretical support for shield axis deviation calculation and shield attitude control strategy.

Evaluating the compaction quality of rock-filled subgrades rapidly and accurately poses a pressing challenge in highway engineering. To address this, this study establishes a discrete element-finite difference coupling model to simulate the response of rock-filled subgrades under impact loading. The primary parameters of the model are calibrated using indoor large-scale triaxial tests, and the model's accuracy is verified through comparisons between calculated and field data. Furthermore, this study conducts an in-depth analysis of the dynamic response results of five commonly used gradations of rock-filled subgrades under varying degrees of compaction, discussing the influence of gradation fractal dimension and porosity on subgrade deformation response. The findings are as follows: (1) A good exponential relationship between subgrade porosity and resilient modulus is identified, and the concept of settlement ratio is introduced, with a linear relationship between settlement ratio and subgrade porosity being verified. It is suggested that both resilient modulus and settlement ratio should be used as control indicators when evaluating subgrade compaction quality. (2) A prediction function for subgrade resilient modulus considering fill gradation and porosity is obtained, revealing that particle gradation has a significant impact on resilient modulus. Specifically, as the gradation fractal dimension approaches 2.31, the resilient modulus increases more rapidly with decreasing porosity. (3) A settlement ratio of zero corresponds to the ideal compaction state of the subgrade. This study establishes a prediction model for the critical resilient modulus of the subgrade in its ideal state, considering fill gradation, and finds that the critical modulus first increases and then decreases with increasing fractal dimension D, reaching a maximum when D=2.34. These findings aim to provide new methods and theories for evaluating the compaction quality of rock-filled subgrades in engineering.

Stratified rock slopes are prone to damage under strong earthquake, leading to geological disasters such as crumbling, landslides and debris flow, and their stability evaluation and support structure optimization are key issues for engineering construction and academic research. Based on field investigations, theoretical analyses, numerical simulations and physical model tests in strong earthquake regions, scholars at home and abroad have carried out a lot of fruitful researches on the damage mechanism and reinforcement measures of rock slopes in strong earthquake regions. Starting from four aspects, including destabilization and damage characteristics of laminated rock slopes, types of support structures, reinforcement mechanisms of support structures and new seismic support structures, the research status of rock slope support structures under strong earthquakes is systematically reviewed, the shortcomings in the current basic research and technical methods of support structures are indicated, and the future research and development directions of seismic support structures for slopes are prospected. This study provides theoretical support for revealing the instability mechanisms and reinforcement strategies of stratified rock slopes in strong earthquake regions, while establishing a scientific foundation for developing more reliable support structures.

The cohesion effect induced by liquid bridges between particles, which promotes particle aggregation, is widespread in both natural environments and engineering applications. Understanding the migration and clogging processes of particles in fractured media under the influence of capillary-cohesion is crucial for advancing particle transport knowledge. Through visualization experiments and seepage calculations, the processes of capillary-cohesive particle migration and clogging are studied. A phase diagram of clogging patterns in the space of capillary-cohesion and flowrate is proposed. Experimental results show that capillary-cohesion induces particle agglomeration, increasing effective particle diameter and significantly enhancing fracture clogging. Stripe-like clogging patterns occur at high flow rates, while complete clogging patterns or entrance sealing patterns occur at low flow rates. Hydrodynamic analysis reveals that fluid velocity distributions control the growth of clogging stripes and the change in residual flow channels in the complete clogging patterns. Furthermore, Smoluchowski theory effectively describes the linear growth behavior of clogging stripes over time. These findings elucidate the mechanism of capillary-cohesive particle migration and clogging in rock fractures, providing theoretical and technical guidance for evaluating and controlling particle transport in fractured media.

Combining bridge row piles with energy piles to create energy row piles can harness shallow geothermal energy for bridge deck deicing in winter and cooling in summer, respectively, while also supporting the mechanical loads of the bridge deck. This study investigates the thermo-mechanical response of energy row piles under heating-cooling cycles through field tests, and analyzes the interactions among energy row piles, slab, and unheated piles. An interface model considering the cyclic shear characteristics of the pile-soil interface is developed in a finite element software, and thermo-mechanical coupling numerical models of energy row piles are established to further explore the changes and mechanisms of long-term settlement of energy row pile under the combined effect of mechanical loads and heating-cooling cycles. The findings reveal that interactions among the energy row pile, slab, and unheated piles can result in load redistribution, leading to high thermally induced stresses of approximately 80% of the maximum thermally induced stress of the energy row pile (i.e. 1.1 MPa) at the top of the energy row piles due to strong restraining effects. Meanwhile, the slab experiences tensile stress exceeding the tensile strength of C30 concrete, reaching 3.75 MPa. Moreover, when the mechanical load is large, energy row piles progressively develop long-term settlement with an increasing number of thermal cycles, exhibiting a negative exponential growth pattern. This phenomenon is attributed to the mechanical load driving the pile-soil interface toward its limiting state, where cyclic shear at the interface readily induces plastic shear displacements, ultimately resulting in the long-term settlement of the energy row piles.

The problem of soil cutting widely exists in engineering fields such as tunnelling, port and waterway dredging, geological drilling, and civil construction. Accurately characterizing the three-dimensional soil failure surface in front of the cutting tool during the soil cutting process is of great significance for analyzing soil disturbance states, evaluating tool cutting performance, and understanding soil-tool interaction mechanisms. A nonlinear elastoplastic damage-based constitutive model is employed to describe the deformation and failure process of soil. Based on the characteristics of damage energy dissipation per unit area of the soil medium, a new numerical method is proposed to directly characterize the three-dimensional soil failure surface. Numerical simulations of flat-tool cutting processes under various operating conditions verify the effectiveness and robustness of the proposed characterization method. The influence of cutting angle and depth on the width, rupture distance, soil disturbance area, and shear failure angle of the three-dimensional soil failure surface is discussed in combination with theoretical calculations. Furthermore, the shape of the three-dimensional soil failure surface for complex-shaped tools obtained through this numerical method is consistent with experimental results, further validating the applicability of the proposed method for complex tool scenarios.

Vacuum preloading, as a widely adopted ground improvement method for saturated soft soils with high water content, is extensively applied in large-scale coastal reclamation projects. However, post-reinforcement bearing capacity remains insufficient in many engineering cases, particularly with limited strength improvement in deep soil layers. Numerous studies have demonstrated that the consolidation efficiency of vacuum preloading is constrained by two critical factors: depth-dependent attenuation of vacuum pressure and fine particle enrichment-induced clogging of drainage paths near prefabricated vertical drains. To address these challenges, this study integrates electro-osmosis with vacuum preloading (EVP) during the later stage of vacuum preloading in the dredger fill project of Yueqing Bay North Port Area. A large-scale model test pool was employed, where conventional vacuum preloading was conducted for 108 days until settlement stabilization, followed by a two-phase EVP intervention. The first phase lasted 11 days, after which electrode polarity was reversed for the second phase (6.5 days), totaling 17.5 days of EVP reinforcement. Post-EVP results revealed significant improvements: at depths of 20 cm, 60 cm, and 100 cm, soil water content decreased by 4.2%,4.84%, and 2.34%, respectively, while vane shear strength increased by 32%, 75%, and 61.1%. The test results indicate that superimposing the electro-osmosis method during the later stage of vacuum preloading can achieve a significant improvement in vane shear strength (with a water content reduction of less than 5%). Particularly for deep soil layers with low initial strength that are difficult to reinforce solely by vacuum preloading, the strength increased by 61%−75%, demonstrating effective reinforcement performance.

Tectonic fault cores are formed substantially of clay minerals. Even a slight change in mineral composition or in water saturation can result in a significant alteration of the sliding regime on the fault. We present results of laboratory experiments on a slider model set-up that was used to study the regularities of slip behavior in a model fault filled with gouge. The gouge consisted of quartz sand and clays of different types (bentonite, illite and kaolinite). The slip behavior essentially depended on gouge mineralogy. The accumulated stress could release via both fast and slow slips. The scaled kinetic energy for fast slips was 10⁻⁵–10⁻³, while that for the slowest slips was 10⁻⁹–10⁻⁷. Fast stick-slip is characteristic of model faults filled with quartz sand in dry and moistened conditions. A gradual transformation from stick-slip to stable sliding was observed for quartz sand/clay gouge as the clay content approached 20%. Under moistening clay, mineralogy played a key role. If the illite clay content was 5%, the moistening led to an increase in peak velocity by more than an order of magnitude; if the bentonite clay was 5%, it led to stabilization of sliding. While alteration in friction coefficient after moistening remained relatively small, the scaled kinetic energy could vary by several orders of magnitude.

In the constructural backfill mining, the composite bearing structure of 'backfill body-immediate roof' structure will be subjected to different loading rates depending on the mining speed and other conditions. According to the loading rate of 0.15−2.40 mm/min, the uniaxial compression test of five groups of rock-backfill composite were carried out, and digital image correlation technology and acoustic emission monitoring were carried out to analyze the evolutionary characteristics of its energy loss. It can be seen from the experiment that the strength of siltstone is significantly greater than the strength of the rock-backfill composite and the backfill body, and the strength of the combination is closer to the strength of the filling body than the siltstone. It can be seen that 0.60 mm/min is the critical load for this group of experiments. When the loading rate of the rock-backfill composite is 0.15−0.60 mm/min, the rock-backfill composite ultimately realizes the synergistic deformation of the siltstone and the backfill body in the rock-backfill composite and destruction of the rock-backfill composite in the process of loading, and when the loading rates are 1.20−2.40 mm/min, rock-backfill composite failed to achieve the collaborative deformation damage of the siltstone and the backfill body parts. When the loading rate is lower than 0.60 mm / min, due to the strength difference between the siltstone and the filling body and the non-uniform deformation of the contact interface between the two, a large crack penetrates the whole specimen. It can be seen that the final failure mode of each group of specimens is a tensile and shear mixed failure mode. By analyzing the dissipation energy changes of the rock-backfill composite and the backfill body, it can be seen that when the loading rate is greater than the critical loading rate, the pre-peak dissipation ratio of the rock-backfill composite is greater than that of the backfill body, and the composite can be destroyed in a coordinated manner. By calculating the energy storage coefficient and energy storage limit of the rock-backfill composite under different loading rates, it is found that when the loading rate is less than 0.60 mm/min, the higher the loading rate, the higher the energy storage limit of the combination specimen, and the speed of absorbing elastic energy is also rising synchronously. Finally, the backfill body part is destroyed first, and the energy released by the instantaneous damage is transmitted to the siltstone part of the rock-backfill composite, so that the elastic energy absorbed by the siltstone part can reach the energy storage limit. The crack in the backfill body part extends into the sandstone to achieve synergistic damage. The results of this study are intended to provide suggestions for ensuring the stability of the composite bearing structure of ' backfill body-immediate roof 'structure under different mining and filling rates.

Current machine learning models for recognizing geological conditions during shield tunneling heavily rely on precise geological data labelling, limiting their applicability in complex geological environments. To address this, we propose a continuous dynamic time warping (CDTW)-based agglomerative hierarchical clustering model (CDTW-Agglomerative), which integrates a linear interpolation framework to overcome DTW's discretization issues. An online learning mechanism is implemented for dynamic strata recognition. The model's accuracy and reliability are validated using Xiamen Metro Line 3 data, with generalization tested on Line 6 data. Results show recognition accuracies of 85% and 73% on the two datasets, demonstrating robust generalization. CDTW-Agglomerative outperforms DTW-Agglomerative, SoftDTW-Agglomerative, and CDTW-based models (K-means, K-medoids, Spectral clustering). Notably, it identifies cutterhead stratigraphy without requiring pre-labelled geological data, supporting intelligent decision-making for tunnelling parameters.