Internal erosion induces alterations in the initial microstructure of soils, simultaneously affecting physical, hydraulic, and mechanical properties. The initial soil composition plays a crucial role in governing the initiation and progression of seepage-induced suffusion. This study employs the controlled variable method to develop granular soil models with varying particle size ratios, initial fine particle contents, and coarse particle shapes. Seepage suffusion simulations coupled with microstructural analyses are conducted using the CFD-DEM approach. Results demonstrate that particle size ratio, fine particle content, and coarse particle shape exert distinct influences on cumulative erosion mass, fine particle distribution, contact fabric, and mechanical redundancy at both macroscopic and microscopic scales. This numerical investigation advances the fundamental understanding of internal erosion mechanisms and informs the development of micro-mechanical constitutive models. Furthermore, for binary granular media composed of coarse and fine particles, careful control of the particle size ratio and fine content is recommended when utilizing gap-graded soils in embankment and dam construction to improve structural resilience and resistance to internal erosion.

The cyclic behavior of clay significantly influences the dynamic response of offshore wind turbines (OWTs). This study presents a practical bounding surface model capable of describing both cyclic shakedown and cyclic degradation. The model is characterized by a simple theoretical framework and a limited number of parameters, and it has been numerically implemented in ABAQUS through a user-defined material (UMAT) subroutine. The yield surface remains fixed at the origin with isotropic hardening, while a movable projection center is introduced to capture cyclic hysteresis behavior. Cumulative plastic deviatoric strain is integrated into the plastic modulus to represent cyclic accumulation. Validation against undrained cyclic tests on three types of clay demonstrates its capability in reproducing stress-strain hysteresis, cyclic shakedown, and cyclic degradation. Additionally, its effectiveness in solving finite element boundary value problems is verified through centrifuge tests on large-diameter monopiles. Furthermore, the model is adopted to analyze the dynamic response of monopile OWTs under seismic loading. The results indicate that, compared to cyclic shakedown, cyclic degradation leads to a progressive reduction in soil stiffness, which diminishes acceleration amplification, increases settlement accumulation, and results in higher residual excess pore pressure with greater fluctuation. Despite its advantages, this model requires a priori specification of the sign of the plastic modulus parameter cd to capture either cyclic degradation or shakedown behavior. Furthermore, under undrained conditions, the model leads pstabilization of the effective stress path, which subsequently results in underestimation of the excess pore pressure.

Accurately modeling soil-fluid coupling under large deformations is critical for understanding and predicting phenomena such as slope failures, embankment collapses, and other geotechnical hazards. This topic has been studied for decades and remains challenging due to the nonlinear responses of geotechnical structures, which typically result from plastic yielding and finite deformation of the soil skeleton. In this work, we comprehensively summarize the theory involved in the soil-fluid coupling problem. Within a finite strain framework, we employ an elasto-plastic constitutive model with linear hardening to represent the solid skeleton and a nearly incompressible model for water. The water content influences the behavior of the solid skeleton by affecting its cohesion. The governing equations are discretized by material point method and two sets of material points are employed to independently represent solid skeleton and fluid, respectively. The proposed method is validated by comparing simulation results with experimental results for the impact of water on dry soil and wet soil. The capability of the method is further demonstrated through two cases: (1) the impact of a rigid body on saturated soil, causing water seepage, and (2) the filling of a ditch, which considers the erosion of the foundation. This work may provide a versatile tool for analyzing the dynamic responses of fluid and solid interactions, considering both mixing and separation phenomena.

This paper presents a rigorous, semi-analytical solution for the drained cylindrical cavity expansion in transversely isotropic sand. The constitutive model used for the sand is the SANISAND-F model, which is developed within the anisotropic critical state theory framework that can account for the essential fabric anisotropy of soils. By introducing an auxiliary variable, the governing equations of the cylindrical expansion problem are transformed into a system of ten first-order ordinary differential equations. Three of these correspond to the stress components, three are associated with the kinematic hardening tensor, three describe the fabric tensor, and the last one represents the specific volume. The solution is validated through comparison with finite element analysis, using Toyoura sand as the reference material. Parametric analyses and discussion on the impact of initial void ratio, initial mean stress level, at-rest earth pressure coefficient and initial fabric anisotropy intensity are presented. The results demonstrate that the fabric anisotropy of sand significantly influences the distribution of stress components and void ratio around the cavity. When fabric anisotropy is considered, the solution predicts lower values of radial, circumferential and vertical stresses near the cavity wall compared to those obtained without considering fabric anisotropy. The proposed solution is expected to enhance the accuracy of cavity expansion predictions in sand, which will have significant practical applications, including interpreting pressuremeter tests, predicting effects of driven pile installation, and improving the understanding of sand mechanics under complex loading scenarios.

Expansive soil, characterized by significant swelling-shrinkage behavior, is prone to cracking under wet-dry cycles, severely compromising engineering stability. This study combines experimental and molecular dynamics (MD) simulation approaches to systematically investigate the improvement effects and micromechanisms of polyvinyl alcohol (PVA) on expansive soil. First, direct shear tests were conducted to analyze the effects of PVA content (0 %-4 %) and moisture content (30 %-50 %) on the shear strength, cohesive force, and internal friction angle of modified soil. Results show that PVA significantly enhances soil cohesive force, with optimal improvement achieved at 3 % PVA content. Second, wet-dry cycle experiments revealed that PVA effectively suppresses crack propagation by improving tensile strength and water retention. Finally, molecular dynamics simulations uncovered the distribution of PVA between montmorillonite (MMT) layers and its influence on interfacial friction behavior. The simulations demonstrated that PVA forms hydrogen bonding networks, enhancing interlayer interactions and frictional resistance. The improved mechanical performance of PVAmodified soil is attributed to both nanoscale bonding effects and macroscale structural reinforcement. This study provides theoretical insights and technical support for expansive soil stabilization.

In geotechnical engineering, the development of efficient and accurate constitutive models for granular soils is crucial. The micromechanical models have gained much attention for their capacity to account for particle-scale interactions and fabric anisotropy, while requiring far less computational resources compared to discrete element method. Various micromechanical models have been proposed in the literature, but none of them have been conclusively shown to agree with the critical state theory given theoretical proof, despite the authors described that their models approximately reach the critical state. This paper modifies the previous CHY micromechanical model that is compatible with the critical state theory based on the assumption that the microscopic force-dilatancy relationship should align with the macroscopic stress-dilatancy relationship. Moreover, under the framework of the CHY model, the fabric anisotropy can be easily considered and the anisotropic critical state can be achieved with the introduction of the fabric evolution law. The model is calibrated using drained and undrained triaxial experiments and the results show that the model reliably replicates the mechanical behaviors of granular materials under both drained and undrained conditions. The compatibility of the model with the critical state theory is verified at both macroscopic and microscopic scales.

Pile foundations are frequently used in the construction of bridges, offshore platforms, and offshore wind turbines, which are often subjected to complex lateral cyclic loading from wind, wave, or current. These lateral loads usually come from different directions or constantly change their direction, which is ignored by most existing calculation models. A two-dimensional p -y model is proposed in this study for the lateral response of the pile subjected to multi-directional cyclic loading in sand. Without introducing additional parameters, the p -y response in two dimensions is coupled by developing the model within the framework of the bounding surface p -y model. Combined with the collapse and recompression model, the effect of sand collapse around the pile during cyclic loading is considered to approach reality. The pile lateral displacement and soil resistance are obtained in incremental form using the finite difference method in the two-dimensional case. By comparing with the model test results, it is demonstrated that the proposed model is able to reasonably predict the lateral cyclic response of the pile as well as the effects of multi-directional cyclic loading. The distribution and variation characteristics of the soil resistance are further discussed by analyzing the results calculated by the proposed model.

This study presents a novel micromorphic continuum model for sand-gravel mixtures with low gravel contents, which explicitly accounts for the influences of the particle size distribution, gravel content, and fabric anisotropy. This model is rigorously formulated based on the principle of macro-microscopic energy conservation and Hamilton's variational principle, incorporating a systematic analysis of the kinematics of coarse and fine particles as well as macro-microscopic deformation differentials. Dispersion equations for plane waves are derived to elucidate wave propagation mechanisms. The results demonstrate that the model effectively captures normal dispersion characteristics and size-dependent effects on wave propagation in these mixtures. In long-wavelength regimes, wave velocities are governed by macroscopic properties, whereas decreasing wavelengths induce interparticle scattering and multiple reflections, attenuating velocities or inhibiting waves, especially when wavelengths approach interparticle spacing. The particle size, porosity, and stiffness ratio primarily influence the macroscopic average stiffness, exhibiting consistent effects on dispersion characteristics across all wavelength domains. In contrast, the particle size ratio and gravel content simultaneously influence both macroscopic mechanical properties and microstructural organization, leading to opposing trends across different wavelength ranges. Model validation against experiments confirms its exceptional predictive ability regarding wave propagation characteristics, including relationships between lowpass threshold frequency, porosity, wave velocity, and coarse particle content. This study provides a theoretical foundation for understanding wave propagation in sand-gravel mixtures and their engineering applications.

Buried pipelines are essential for the safe and efficient transportation of energy products such as oil, gas, and various chemical fluids. However, these pipelines are highly vulnerable to ground movements caused by geohazards such as seismic faults, landslide, liquefaction-induced lateral spreading, and soil creep, which can result in potential pipeline failures such as leaks or explosions. Response prediction of buried pipelines under such movements is critical for ensuring structural integrity, mitigating environmental risks, and avoiding costly disruptions. As such, this study adopts a Physics-Informed Neural Networks (PINNs) approach, integrated with a transfer learning technique, to predict structural response (e.g., strain) of both unreinforced and reinforced steel pipes subjected to Permanent Ground Displacement (PGD). The PINN method offers a meshless, simulation-free alternative to traditional numerical methods such as Finite Element Method (FEM) and Finite Difference Method (FDM), while eliminating the need for training data, unlike conventional machine learning approaches. The analyses can provide useful information for in-service pipe integrity assessment and reinforcement, if needed. The accuracy of the predicted results is verified against Finite Element (FE) and Finite Difference (FD) methods, showcasing the capability of PINNs in accurately predicting displacement and strain fields in pipelines under geohazard-induced ground movement.

Mesh-free methods, such as the Smooth Particle Hydrodynamics (SPH) method, have recently been successfully developed to model the entire wetting-induced slope collapse process, such as rainfall-induced landslides, from the onset to complete failure. However, the latest SPH developments still lack an advanced unsaturated constitutive model capable of capturing complex soil behaviour responses to wetting. This limitation reduces their ability to provide detailed insights into the failure processes and to correctly capture the complex behaviours of unsaturated soils. This paper addresses this research gap by incorporating an advanced unsaturated constitutive model for clay and sand (CASM-X) into a recently proposed fully coupled seepage flow-deformation SPH framework to simulate a field-scale wetting-induced slope collapse test. The CASM-X model is based on the unified critical state constitutive model for clay and sand (CASM) and incorporates a void-dependent water retention curve and a modified suction-dependent compression index law, enabling the accurate prediction various unsaturated soil behaviours. The integration of the proposed CASM-X model in the fully coupled flow deformation SPH framework enables the successful prediction of a field-scale wetting-induced slope collapse test, providing insights into slope failure mechanisms from initiation to post-failure responses.