For construction quality control, the compaction delay referred to as mellowing time (MT) is crucial for achieving the desired outcomes of the chemical soil stabilization process in the field. In the current study, fly ash-based geopolymer (GFA) is used as a chemical stabilizer for expansive clay because of its significance in resource utilization and waste repurposing for soil stabilization through an enhanced process. The MT-influenced macroscopic physicomechanical properties and microstructural and mineralogical properties of expansive clay treated with varying GFA and curing period (CP) were investigated. The significant amelioration of strength and compression properties is observed through the unconfined compression test, California bearing ratio test, and one-dimensional (1D) consolidation test with an increase in GFA content and CP. This improvement is caused by the formation of cementitious [(N, C)-A-S-H] compounds as confirmed by SEM, EDAX, and XRD analyses. Meanwhile, as the MT increases, a decline in both the strength and compression characteristics of the GFA-treated specimens is observed. However, these specimens exhibit a reversal in deformability and brittleness with an increase in MT, which can be attributed to the development of a porous aggregated soil structure resulting from initial hydration before densification. In addition, a generalized mathematical modeling framework was established based on three-dimensional (3D) response surface modeling to quantify the MT-influenced strength and brittleness-related characteristics using MT, GFA, and CP as predictors. The established mathematical framework showed generality and reasonable accuracy in the prediction based on the experimental data. This article outlines the implications for practitioners and researchers of using GFA for the stabilization of expansive clay considering MT-influenced mechanical characteristics in the field.

This article presents a micro-structure tensor enhanced elasto-plastic finite element (FE) method to address strength anisotropy in three-dimensional (3D) soil slope stability analysis. The gravity increase method (GIM) is employed to analyze the stability of 3D anisotropic soil slopes. The accuracy of the proposed method is first verified against the data in the literature. We then simulate the 3D soil slope with a straight slope surface and the convex and concave slope surfaces with a 90 degrees turning corner to study the 3D effect on slope stability and the failure mechanism under anisotropy conditions. Based on our numerical results, the end effect significantly impacts the failure mechanism and safety factor. Anisotropy degree notably affects the safety factor, with higher degrees leading to deeper landslides. For concave slopes, they can be approximated by straight slopes with suitable boundary conditions to assess their stability. Furthermore, a case study of the Saint-Alban test embankment A in Quebec, Canada, is provided to demonstrate the applicability of the proposed FE model. (c) 2025 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Published by Elsevier B.V. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/ 4.0/).

Microwave-assisted rock-breaking technology, as a novel hybrid approach, is anticipated to facilitate the efficient excavation of complex rock formations. It is therefore crucial to understand the damage and failure mechanisms of rocks that have been subjected to irradiation. In this study, uniaxial compression experiments were conducted on granite specimens after 1.4 kW microwave irradiation for varying durations. Furthermore, a numerical method was proposed to solve electromagnetic-thermal-mechanical coupling problems by integrating finite and discrete elements. The results demonstrated a differential temperature distribution (high temperature in the middle and low-temperature areas at the ends) in the granite specimens under microwave irradiation, which resulted in a notable reduction in their physical and mechanical properties. As the duration of irradiation increased, the rate of heating and the extent of strength reduction both diminished, while the morphology and distribution of cracks at ultimate failure became increasingly complex. The numerical method effectively addresses the simulation challenges associated with the electromagnetic selective heating of granite containing multiple polar minerals under microwave irradiation. This approach accounted for the non-uniform thermal expansion of the minerals and provided a comprehensive model of damage progression under compression. (c) 2025 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Published by Elsevier B.V. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/ 4.0/).

Weak structural plane deformation is responsible for the non-uniform large deformation disasters in layered rock tunnels, resulting in steel arch distortion and secondary lining cracking. In this study, a servo biaxial testing system was employed to conduct physical modeling tests on layered rock tunnels with bedding planes of varying dip angles. The influence of structural anisotropy in layered rocks on the micro displacement and strain field of surrounding rocks was analyzed using digital image correlation (DIC) technology. The spatiotemporal evolution of non-uniform deformation of surrounding rocks was investigated, and numerical simulation was performed to verify the experimental results. The findings indicate that the displacement and strain field of the surrounding layered rocks are all maximized at the horizontal bedding planes and decrease linearly with the increasing dip angle. The failure of the layered surrounding rock with different dip angles occurs and extends along the bedding planes. Compressive strain failure occurs after excavation under high horizontal stress. This study provides significant theoretical support for the analysis, prediction, and control of non-uniform deformation of tunnel surrounding rocks. (c) 2024 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Published by Elsevier B.V. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/ 4.0/).

To achieve the loading of the stress path of hard rock, the spherical discrete element model (DEM) and the new flexible membrane technology were utilized to realize the transient loading of three principal stresses with arbitrary magnitudes and orientations. Furthermore, based on the deep tunnel of China Jinping Underground Laboratory II (CJPL-II), the deformation and fracture evolution characteristics of deep hard rock induced by excavation stress path were analyzed, and the mechanisms of transient loading-unloading and stress rotation-induced fractures were revealed from a mesoscopic perspective. The results indicated that the stress-strain curve exhibits different trends and degrees of sudden changes when subjected to transient changes in principal stress, accompanied by sudden changes in strain rate. Stress rotation induces spatially directional deformation, resulting in fractures of different degrees and orientations, and increasing the degree of deformation anisotropy. The correlation between the degree of induced fracture and the unloading magnitude of minimum principal stress, as well as its initial level is significant and positive. The process of mechanical response during transient unloading exhibits clear nonlinearity and directivity. After transient unloading, both the minimum principal stress and minimum principal strain rate decrease sharply and then tend to stabilize. This occurs from the edge to the interior and from the direction of the minimum principal stress to the direction of the maximum principal stress on the epsilon 1-epsilon 3 1-epsilon 3 plane. Transient unloading will induce a tensile stress wave. The ability to induce fractures due to changes in principal stress magnitude, orientation and rotation paths gradually increases. The analysis indicates a positive correlation between the abrupt change amplitude of strain rate and the maximum unloading magnitude, which is determined by the magnitude and rotation of principal stress. A high tensile strain rate is more likely to induce fractures under low minimum principal stress. (c) 2024 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting by Elsevier B.V. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/ licenses/by-nc-nd/4.0/).

To investigate the long-term stability of deep rocks, a three-dimensional (3D) time-dependent model that accounts for excavation-induced damage and complex stress state is developed. This model comprises three main components: a 3D viscoplastic isotropic constitutive relation that considers excavation damage and complex stress state, a quantitative relationship between critical irreversible deformation and complex stress state, and evolution characteristics of strength parameters. The proposed model is implemented in a self-developed numerical code, i.e. CASRock. The reliability of the model is validated through experiments. It is indicated that the time-dependent fracturing potential index (xTFPI) at a given time during the attenuation creep stage shows a negative correlation with the extent of excavationinduced damage. The time-dependent fracturing process of rock demonstrates a distinct interval effect of the intermediate principal stress, thereby highlighting the 3D stress-dependent characteristic of the model. Finally, the influence of excavation-induced damage and intermediate principal stress on the time-dependent fracturing characteristics of the surrounding rocks around the tunnel is discussed. (c) 2024 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting by Elsevier B.V. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/ licenses/by-nc-nd/4.0/).

Understanding the pore water pressure distribution in unsaturated soil is crucial in predicting shallow landslides triggered by rainfall, mainly when dealing with different temporal patterns of rainfall intensity. However, the hydrological response of vegetated slopes, especially three-dimensional (3D) slopes covered with shrubs, under different rainfall patterns remains unclear and requires further investigation. To address this issue, this study adopts a novel 3D numerical model for simulating hydraulic interactions between the root system of the shrub and the surrounding soil. Three series of numerical parametric studies are conducted to investigate the influences of slope inclination, rainfall pattern and rainfall duration. Four rainfall patterns (advanced, bimodal, delayed, and uniform) and two rainfall durations (4h intense and 168-h mild rainfall) are considered to study the hydrological response of the slope. The computed results show that 17% higher transpiration-induced suction is found for a steeper slope, which remains even after a short, intense rainfall with a 100-year return period. The extreme rainfalls with advanced (PA), bimodal (PB) and uniform (PU) rainfall patterns need to be considered for the short rainfall duration (4 h), while the delayed (PD) and uniform (PU) rainfall patterns are highly recommended for long rainfall durations (168 h). The presence of plants can improve slope stability markedly under extreme rainfall with a short duration (4 h). For the long duration (168 h), the benefit of the plant in preserving pore-water pressure (PWP) and slope stability may not be sufficient.

The tangential displacement amplitude determines the mobilization of shear strength and the deforming and sliding displacements of soil-structure interfaces, and therefore plays a crucial role in the interface behavior. A series of three-dimensional (3D) simple-shear interface tests were conducted between gravel and steel to investigate the influence of tangential displacement amplitude on the tangential deformation, volumetric change, and shear strength. Test results show that deforming and sliding displacements are distinctly induced by shearing. The deforming displacement migrates toward the initial shear direction, caused by the shear orientation effect, and the migration becomes magnified and then stabilizes with cyclic shearing. The shear strength would not be mobilized when the tangential displacement amplitude is relatively small. It behaves in an anisotropic manner if mobilized and gradually degrades as cyclic shearing continues, attributed to the dominant particle crushing over the shear densification effect. Two critical tangential displacement amplitudes are found for the mobilization of shear strength, and determine whether the shear strength could be mobilized during cyclic shearing and whether it is immediately mobilized at initial shearing, respectively. The tangential displacement amplitude primarily affects deforming and sliding amplitudes and their migration, shear stiffness, irreversible and peak reversible normal displacements, peak and residual cyclic shear strength, and anisotropy extent, instead of their relationship patterns. An increased tangential displacement amplitude results in magnified sliding amplitude, decreased deforming weight, accelerated degradation of deforming amplitude, and magnified migration of deforming and sliding displacements. Additionally, large tangential displacement amplitude leads to large irreversible and peak reversible normal displacements, small shear stiffness, and small peak and residual cyclic shear strength. The peak reversible normal displacement is determined by and has a linear relationship with the deforming displacement, and the irreversible normal displacement presents perfect consistency behavior against shear work density, regardless of tangential displacement amplitude. The consistency behavior could be well described using a hyperbolic model, which significantly simplifies the 3D constitutive modeling of the interface.

Previous earthquake events indicate that pile foundations in liquefiable soils are vulnerable to damage due to the coupling of inertial and kinematic effects. Inclined piles are widely applied in structures located in liquefiable soils, but few investigations of the coupling of the superstructure-pile inertial and soil-pile kinematic effects have been conducted. To address this gap, this study adopted a three-dimensional (3D) numerical model to investigate the coupling of inertial and kinematic effects in pile foundations with different inclination angles. The pile head bending moment was employed to represent the pile response, while the soil surface displacement and structure acceleration were utilized to quantify the kinematic and inertial effects. The role of the inclination angle on the interactions between inertial and kinematic effects is herein considered for pile groups. In particular, the inertial effect significantly influences the behavior of pile groups with larger inclination angles, whereas the kinematic effect predominates the pile head moment in vertical pile groups. In this paper, the influence of the pile inclination angle, superstructure configuration, and earthquake intensity on the interactions was investigated. The principal findings revealed that the kinematic effect dominates in the vertical pile group irrespective of the properties of the superstructure, while the inertial effect plays a significant role in the response of the inclined pile groups, especially for superstructures with considerable heights. Inclined piles are vulnerable to damage due to the interaction of inertial and kinematic effects during earthquakes. This study conducted a series of three-dimensional (3D) finite-element simulations to investigate the interaction of inertial and kinematic effects in pile foundations with different inclination angles. The influence of pile inclination angle, superstructure height, and earthquake characteristics was investigated. In current practices, various codes and pseudostatic methods have been adopted to sum a percentage of the inertia-induced bending moment and another percentage of the kinematic-induced bending moment. This study indicates that under certain conditions, the simple summing of the bending moment induced by the inertial and kinematic effects could be inaccurate. The present study identified several factors that influence the interaction of inertial and kinematic effects on piles with different inclination angles. The inclined piles in liquefied soil, especially for supporting tall and heavy superstructure, attention should be given to the influence of inertial effect on the pile head bending moment.

In recent decades, research on renewable energy has been boosted by the emerging awareness of energy security and climate change and their consequences, such as the global cost of adapting to the climate impacts. Both onshore and offshore wind turbine farms have been considered as one of the main alternatives to fossil fuels. Their development currently involves seismic-prone areas, such as the Californian coastline and East Asia, where the risk of soil liquefaction is significant. Onshore wind turbines (OWTs) typically are founded on shallow rafts. Their operation can be affected strongly by the simultaneous presence of intense earthquakes and wind thrust, which may cause remarkable permanent tilting and loss of serviceability. In these conditions, accurate evaluation of the seismic performance of these structures requires the development of well-validated numerical tools capable of capturing the cyclic soil behavior and the build-up and contextual dissipation of seismic-induced pore-water pressures. In this paper, a numerical model developed in OpenSees, calibrated against the results of dynamic centrifuge tests, was used to evaluate the influence of some ground motion intensity Measures of the seismic behavior of OWTs included the amplitude, frequency content, strong-motion duration, and Arias intensity (energy content) of the earthquake, together with the effect of a coseismal wind thrust, which is not well explored in the literature. The seismic performance of an OWT was assessed in terms of peak and permanent settlement and tilting, the latter of which was compared with the threshold of 0.5 degrees typically adopted in practice.