Earthquake-induced soil liquefaction causes ground and foundation failures, and it challenges the scientific community to explore the liquefaction problem in deep deposit under strong shaking. Due to the capacity limitation of physical modelling facility, it is difficult to reproduce soil liquefaction response of deep sand ground by centrifuge shaking table test. To address this problem, a suite of centrifuge model tests were conducted with the aid of Iai's Type III generalized scaling law (i.e., GSL) to observe the liquefaction response of deep sand ground, where Models 1 and 2 were used to validate the GSL and Model 3 with the validated GSL stands for the deep sand ground with prototype depth of 80 m. The test results of Models 1 and 2 indicate that GSL generally performs well for small-strain shear modulus, nonlinear dynamic response of acceleration and the generation of excess pore water pressure, but leaves considerable errors for post-shaking dissipation process and ground settlement with large plastic strain. The test results of Model 3 indicate that liquefaction is also possible in depth of 30-40 m under shaking event of PBA = 0.4 g and Mw = 7.5. For deeper depth without triggering of liquefaction, a depthdependent power function relationship between the peak excess pore water pressure and Arias intensity has been established. The test results also revealed that consolidation and earthquake shaking history contribute to the development of soil anisotropy in a deep ground, leading to a continuous increase of anisotropy degree, which could be evaluated using the small-strain shear moduli in different stress planes under orthogonal stress conditions.

To better characterize the intricate coupled thermo-hydro-mechanical dynamic (THMD) response in twodimensional saturated soil and to enrich the research object of Green-Naghdi (G-N) generalized thermoelastic theory, this study innovatively combines the G-N generalized thermoelastic theory and Caputo's fractional order derivative, to obtain the new control equations, and to establish a new fractional order thermoelastic theoretical model. The article is solved by the normal mode analysis (NMA), which can eliminate the integration error and solve the complex fractional order partial differential control equations quickly at the same time. The effects of different boundary conditions of fractional order derivatives, porosity, frequency, and thermal conductivity coefficients on non-dimensional excess pore water pressure, temperature, vertical displacement, and vertical stress are also fully analyzed, and the distribution curves of high precision numerical solutions are given. The results show that the effect of frequency variation on each non-dimensional variable is obvious. The effects of fractional order derivatives, porosity and thermal conductivity coefficients on the non-dimensional variables vary depending on the boundary conditions. The results provide theoretical support for geotechnical and environmental engineering.

Rockfill, a coarse granular material commonly used in dam construction, exhibits complex mechanical behavior under generalized stress conditions. This paper investigates the mechanical properties of rockfill through a series of stress-path tests conducted on a self-developed, large-scale true triaxial apparatus with cubical specimens of 60 x 30 x 30cm. Three test series are carried out by varying the mean effective stress, the deviator stress and the Lode's angle, respectively. An elastoplastic constitutive model is presented to describe the behavior of rockfill. An improved dilatancy equation is introduced by considering the phase transformation stress ratio instead of the critical stress ratio.

Underground tunnels subjected to asymmetric load or ground conditions are susceptible to experiencing uneven longitudinal bending, shearing, and torsional deformations, which further induce cross sectional flattening and warping. The intrinsic damages caused by multiple deformation modes are critical for tunnel health and safety but have long been neglected in practice. In the paper, a three-dimensional analytical model for soil-tunnel interactions was proposed with multiple-mode deformations incorporated, where the tunnel is assumed as a thin-walled pipe resting on an elastic foundation with five deformation modes: bending, shearing, torsion, warping, and flattening. Besides, a three-dimensional variable soil spring model was adopted, accounting for the strata discontinuities in longitudinal and transverse directions. A finite element solution for the proposed model was derived under arbitrary external loads using the principle of minimum potential energy. The validity of the proposed model was substantiated through three case studies. Based on the model, the coupling relationship of tunnel structure in transverse and longitudinal directions was revealed. Furthermore, parametric analysis was conducted to reveal the impact of tunnel width-to-thickness ratio, soil resistance coefficient, and composite strata on tunnel behaviors. These results significantly contribute to a deeper understanding of the intricate behaviors of tunnels, offering potential advancements for improved tunnel design methodologies.

Stress-strain behavior of two different soil specimens subjected to cyclic compressive loading are studied herein, the goal being to present a simple dynamic uniaxial mem-modeling approach that aids physical insight and enables system identification. In this paper, mem stands for memory, i.e., hysteresis. Mem-models are hysteresis models transferred from electrical engineering using physical analogies. Connected in series, a mem-dashpot and mem-spring are employed to model inter-cycle strain ratcheting and intra-cycle gradual densification of the two soil specimens. Measured time histories of stress and strain are first decomposed so that the two modeling components, mem-dashpot and mem-spring, can be identified separately. This paper focuses on the mem-dashpot, a nonlinear generalization of a linear viscous damper. A mem-spring model is also devised built on an extended Masing model. Nonlinear dynamic simulations (with inertia) employing the identified mem-dashpot and mem-spring demonstrate how well the identified mem-models reproduce the measured early-time data (first 200 cycles of loading). Choices of state variables inherited from bond graph theory, the root of mem-models, are introduced, while MATLAB time integrators (i.e., ode solvers) are used throughout this study to explore a range of contrasting damper and spring models. Stiff solver and the state event location algorithm are employed to solve the equations of motion involving piecewise-defined restoring forces (when applicable). Computational details and results are relegated to the appendices. This is the first study to use single-degree-of-freedom (SDOF) system dynamic simulations to explore the consistency of mem-models identified from real-world data.

Accurate prediction of ground surface settlement (GSS) adjacent to an excavation is important to prevent potential damage to the surrounding environment. Previous studies have extensively delved into this topic but all under the limitations of either imprecise theories or insufficient data. In the present study, we proposed a physics-constrained neural network (PhyNN) for predicting excavation-induced GSS to fully integrate the theory of elasticity with observations and make full use of the strong fitting ability of neural networks (NNs). This model incorporates an analytical solution as an additional regularization term in the loss function to guide the training of NN. Moreover, we introduced three trainable parameters into the analytical solution so that it can be adaptively modified during the training process. The performance of the proposed PhyNN model is verified using data from a case study project. Results show that our PhyNN model achieves higher prediction accuracy, better generalization ability, and robustness than the purely data-driven NN model when confronted with data containing noise and outliers. Remarkably, by incorporating physical constraints, the admissible solution space of PhyNN is significantly narrowed, leading to a substantial reduction in the need for the amount of training data. The proposed PhyNN can be utilized as a general framework for integrating physical constraints into data-driven machine-learning models. (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/).

Cemented sandy gravel is often used to enhance the foundation soil of engineering projects. This paper presents results of triaxial tests on cemented sandy gravel specimens. We compared 8 cemented specimens and 4 uncemented specimens. The strength, dilatancy, and stiffness behavior of both cemented and uncemented specimens are compared. The strength of cemented specimens is significantly greater than that of uncemented specimens, and the cemented specimens demonstrate pronounced expansion characteristics. The peak friction angle of the cemented specimen shows a linear relationship with the confining pressure: psi = 68.1-18.2lg(sigma 3/pa). To quantify the structural strength of the cemented specimens, a structural damage parameter is introduced based on the differences in mechanical properties between the two materials. The structural damage parameter first increases and then decreases as shearing progresses, and a hump curve function is used to describe this behavior. In the frame of the generalized plasticity, a novel elastoplastic model is established, considering the structural parameter as a factor of the plastic modulus, loading vectors and plastic flow direction vectors. The calculated values fit well with the experimental results. The model can reflect the characteristics of cemented sandy gravel, in terms of stress softening, residual strength, and volumetric dilation. Finally, the model is used to evaluate the deformation of a sluice dam foundation after being enhanced with cemented sandy gravel. The results show that after treatment, both the settlement of the gate floor and the shear deformation of the waterstops can be reduced by more than 10%.

Loess is susceptible to loading effects such as significant changes in strength and volume variation caused by loading and wetting. In this study, considering the different connection states of pore water and gas in loess fabric, the gas phase closure case is incorporated into a unified form of the generalized effective stress framework, introducing a damage parameter considering the effects of closed pore gas. The loading effects of unsaturated loess under wide variations in saturation are described in a unified way, and the model performance is verified by corresponding stress and hydraulic path tests. The results indicated that the collapse response involves the initial void ratio of loess, and the coupled outwards motion of the loading-collapse (LC) yield surface under loading enhances its structural strength. Suction-enhanced yield stress requires a greater tensile stress to counteract its structural stability. The nucleation of bubbles at high saturation causes a decrease in yield stress. The loading effect exhibits a smaller collapse behavior when the influence of closed gas is considered, whereas the suction path does not cross the LC in the stress space under hydraulic action for the same parameters, which amplifies the influence of closed gas on loess deformation. (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/).

Seismic events and wave action can induce volumetric strain (ev) accumulation in saturated sandy soils, leading to damage to the ground surface and structures. A quantifiable relationship exists between the generation of ev in sandy soils under drained conditions and the development of pore water pressures under undrained conditions. In this study, the impact of relative density (Dr), cyclic stress path, and stress level on the characteristics of volumetric strain (ev) generation in saturated coral sands (SCS) was evaluated through drained tests employing various cyclic stress paths. The test findings demonstrate that the rate of ev accumulation in SCS is notably affected by the cyclic stress path. The rise in peak volumetric strain (evp) in SCS, as a function of the number of cycles, conforms to the arctangent function model. The unit cyclic stress ratio (USR) was employed as an indicator of complex cyclic loading levels. It was determined that coefficient (evp)u is positively correlated with USR at a specific Dr. At the same Dr, coefficient CN1 exhibits a positive correlation with USR, while coefficient CN2 displays a negative correlation with USR, following a power-law relationship. Irrespective of cyclic loading conditions, evp rises with an increase in generalized shear strain amplitude (yga). A power function model was established to represent the relationship between evp and yga. The coefficient 41 decreases as Dr increases. Comparisons were drawn between evp and yga for Ottawa sand and SCS. The results indicate that, as Dr of Ottawa sand increases from 30 % to 70 %, the coefficient 41 decreases from 1.54 to 0.73, representing a reduction of approximately 53 %. In contrast, under identical conditions, the coefficient 41 of SCS exhibits a less pronounced decrease, from 1.16 to 0.79, corresponding to a reduction of roughly 32 %. These observations suggest that variations in Dr have a more substantial impact on generating evp in Ottawa sand compared to SCS.

Soft clays exhibit significant challenges in geotechnical engineering due to their low permeability, high compressibility, and susceptibility to settlement under applied loads. These geological factors pose unique difficulties in predicting long-term settlement accurately and efficiently, particularly through Class C prediction methods that involve iterative processes with complex numerical models. To address these challenges, this study presents an efficient approach for Class C prediction of long-term settlement in soft clays. This approach integrates Bayesian updating with structural reliability methods (BUS) and the general simplified Hypothesis B method which is a semi-analytical method based on one-dimensional elastic visco-plastic (1D EVP) model. Unlike previous research that used Response Surface Model (RSM) with polynomial function for consolidation evaluation, the proposed approach enhances both accuracy and performance consistency under varying conditions. Additionally, by leveraging analytical solutions instead of iterative small-time steps required by Finite Element Method (FEM) or Finite Difference Method (FDM), the computational efficiency is also enhanced. The effectiveness of the proposed approach is demonstrated through its application to an embankment improved with prefabricated vertical drain (PVD) in Ballina, New South Wales, Australia. Comparative analyses demonstrate that the predicted settlements from this study, using only the monitoring settlement data collected prior to the 76th day of the project, align closely with the results from established RSM and FEM-based Bayesian back analysis approaches. The obtained results also indicate that the predicted settlements, based on 76 days of monitoring data, closely match field measurements at various depths, whether relying solely on settlement data or integrating additional pore water pressure data. For the Ballina embankment, over 40,000 consolidation analyses required for a single BUS simulation can be completed within 10 h using the general simplified Hypothesis B method, compared to months it might take with FEM or FDM approaches. This makes the proposed approach a practical tool for geotechnical engineers, enabling reliable settlement predictions early in the project timeline while maintaining low computational costs.