As the increasing demand for deep mineral resource extraction and the construction of deep vertical shafts by the artificial ground freezing method, the stability and safety of shaft that traverse thick alluvial depend significantly on their interaction with the surrounding deep frozen soil medium. Such interaction is directly conditioned by the mechanical properties of the deep frozen soil. To precisely capture these in-situ mechanical properties, the mechanical parameters tests using remodeled frozen specimens cannot ignore the disparities in consolidation history, stress environment and formation conditions between the deep and shallow soils. This study performs a series of long-term high-pressure K0 consolidation (where K0 represents the static earth pressure coefficient, describing the ratio of horizontal to vertical stress under zero lateral strain conditions), freezing under sustained load and unloading triaxial shear tests utilizing remodeled deep clay. This study presents the response of unloading strength and damage properties under varying consolidation stresses, durations, and freezing temperatures. The unloading strength increases sharply and then stabilizes with consolidation time. The unloading strength shows an approximate linear positive correlation with the consolidation stress, while a negative correlation with the freezing temperature. The strengthening rate of the unloading strength due to freezing temperature tends to decrease with increasing consolidation time. Additionally, an improved damage constitutive model was proposed and validated by incorporating the initial K0 stress state and a Weibull-based assumption for damage elements. Based on the back propagation (BP) neural network, a prediction method for the stress-strain curve was offered according to the consolidation stress level, initial stress state, and temperature. These results can provide references for improving the mechanical testing methods of deep frozen clay and revealing differences in mechanical properties between deep and shallow soils.

The soil strength of soft clay is influenced by strain rate effect. Models considering strain rate effect always ignore the impact of loading rate on pore pressure and have poor applicability to 3D engineering problems. Based on the classic inelastic core boundary surface model, a logarithmic rate function representing the strain rate effect of soft soil was introduced to the hardening law. A new parameter H was added to adjust the plastic modulus while another new parameter mu is introduced to account for the strain rate effect. A rate-effect boundary surface constitutive model suitable for saturated clay was subsequently proposed. Combined with the implicit integral numerical algorithm and stress-permeability coupling analysis, the innovative model was implemented in the finite element software and validated by comparing with the results of triaxial tests. By analysing the rate-effect of 11 types of soft soil, a formula to calculate the rate parameter was derived. The developed model was used to calculate the undrained vertical bearing capacity and sliding resistance of a movable subsea mudmat. The mudmat frictional coefficient from soil undrained to partial drained and finally undrained state was obtained and compared with those from the Modified Cam-Clay model. Identical results were obtained without considering the rate effect. When considering the strain rate effect on the improvement of soil strength, the friction resistance coefficient initially decreases and then increases with the decrease of the sliding speed, eventually stabilising after reaching the limit value. The rate-effect on the friction resistance coefficient is most prominent under undrained conditions with high sliding speeds. The soil strain rate effect is suggested to be considered in the design of the subsea mudmat avoid underestimating the friction resistance.

This study systematically investigated the pore structure response of kaolin and illite/smectite mixed-layer rich clay in a reconstituted state to one-dimensional (1D) compression by first performing oedometer tests on saturated clay slurries, followed by characterising their pore structure using multi-scale characterisation techniques, with the primary objective of advancing the current understanding of the microstructural mechanisms underlying the macroscopic deformation of such clays. Under 1D loading, the volume reduction observed at the macro level essentially represented the macroscopic manifestation of changes in inter-aggregate porosity at the pore scale. It was the inter-particle pores that were compressed, despite the interlayer pores remaining stable. Two distinct pore collapse mechanisms were identified: kaolin exhibited a progressive collapse of particular larger pore population in an ordered manner, whereas illite/smectite mixed-layer rich clay demonstrated overall compression of inter-aggregate pores. Accordingly, mathematical relationships between the porosity and compressibility parameters for these two soils were proposed, with the two exhibiting opposite trends arising from their distinct microstructural features. Approaching from the unique perspective of pore structure, quantitative analysis of pore orientation and morphology on the vertical and horizontal planes demonstrated some progressively increasing anisotropy during compression. These findings provide important insights into porescale mechanisms governing clay compression behaviour and enrich the limited microporosity database in soil mechanics.

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.

A group of earthquakes typically consists of a mainshock followed by multiple aftershocks. Exploration of the dynamic behaviors of soil subjected to sequential earthquake loading is crucial. In this paper, a series of cyclic simple shear tests were performed on the undisturbed soft clay under different cyclic stress amplitudes and reconsolidation degrees. The equivalent seismic shear stress was calculated based on the seismic intensity and soil buried depth. Furthermore, reconsolidation was conducted at the loading interval to investigate the influence of seismic history. An empirical model for predicting the variation of the accumulative dissipated energy with the number of cycles was established. The energy dissipation principle was employed to investigate the evolution of cyclic shear strain and equivalent pore pressure. The findings suggested that as the cyclic stress amplitude increased, incremental damage caused by the aftershock loading to the soil skeleton structure became more severe. This was manifested as the progressive increase in deformation and the rapid accumulation of dissipated energy. Concurrently, the reconsolidation process reduced the extent of the energy dissipation by inhibiting misalignment and slippage among soil particles, thereby enhancing the resistance of the soft clay to subsequent dynamic loading.

The foundation soil below the structure usually bears the combined action of initial static and cyclic shear loading. This experimental investigation focused on the cyclic properties of saturated soft clay in the initial static shear stress state. A range of constant volume cyclic simple shear tests were performed on Shanghai soft clay at different initial static shear stress ratios (SSR) and cyclic shear stress ratios (CSR). The cyclic behavior of soft clay with SSR was compared with that without SSR. An empirical model for predicting cyclic strength of soft clay under various SSR and CSR combinations was proposed and validated. Research results indicated that an increase of shear loading level, including SSR and CSR, results in a larger magnitude of shear strain. The response of pore water pressure is simultaneously dominated by the amplitude and the duration of shear loading. The maximum pore water pressure induced by smaller loading over a long duration may be greater than that under larger loading over a short duration. The initial static shear stress does not necessarily have a negative impact on cyclic strength. At least, compared to cases without SSR, the low-level SSR can improve the deformation resistance of soft clay under the cyclic loading. For the higher SSR level, the cyclic strength decreases with the increase of SSR.

The discrete element method (DEM) has demonstrated significant advantages in simulating soil-tool interaction and an appropriate contact model notable affected the simulation accuracy. The accuracy of numerical simulation is compromised due to the variations in soil properties when tillage implements are employed in clay-moist soil conditions. This study aims to establish a discrete element model of clay-moist soil based on the Edinburgh Elasto-Plastic Adhesion (EEPA) contact model. Calibration tests using a combination of direct shear tests and cone penetration tests were conducted to identify sensitive parameters that need to be calibrated in the model and analyze the effects of each parameter. The results indicated that contact plasticity ratio and surface energy had significant influence on representing the mechanical properties of clay-moist soil. Then, by utilizing scanning technology to acquire furrow shape data, soil bin test was conducted to validate the reliability of the calibration parameters. Using sensitive parameters as variables, the actual value of clay-moist soil with a moisture content of 33 % as the target value obtained from experimental tests. The optimal combination was: the coefficient of static friction of 0.45, the coefficient of rolling friction of 0.18, and the surface energy of 27.95 J.m-2, the contact plasticity ratio of 0.59. The relative error between the simulated draft force value and the actual measured value was 7.98 %, and the relative errors in the furrow type parameters did not exceed 5 %. The accuracy of the calibration results was verified through comparative analysis of simulation and empirical results. This study provides a scientific approach for employing DEM in modeling clay-moist soil-tool interaction.

Pile penetration in soft ground involves complex mechanisms, including significant alterations in the soil state surrounding the pile, which influence the pile negative skin friction (NSF) over time. However, the pile penetration process is often excluded from finite element analysis. This paper investigates the impact of pile penetration on the generation of NSF and dragload. A stable node-based smoothed particle finite element method (SNS-PFEM) framework is introduced for two-dimensional axisymmetric conditions and coupled consolidation, incorporating the ANICREEP model of soft soil with a modified cutting-plane algorithm. A field case study with penetration process is simulated to verify the numerical model's performance, followed by a parametric analysis on the effect of penetration rate on NSF during consolidation. Results indicate that without the pile penetration process in NSF analysis can result in an unsafely low estimation of NSF and dragload magnitudes. The penetration rate affects dragload only at the initial consolidation stage. As consolidation progresses, dragload converges to nearly the same magnitude across different rates. Additionally, current design methods inadequately predict the beta value (where beta is an empirical factor correlating vertical effective stress of soil with the pile skin friction) and its time dependency, for which a new empirical formula for the time-dependent beta value is proposed and successfully applied to other field cases.

Knowledge Gap: The aggregation of clay minerals-layered silicate nanoparticles-strongly impacts fluid flow, solute migration, and solid mechanics in soils, sediments, and sedimentary rocks. Experimental and computational characterization of clay aggregation is inhibited by the delicate water-mediated nature of clay colloidal interactions and by the range of spatial scales involved, from 1 nm thick platelets to flocs with dimensions up to micrometers or more. Simulations: Using a new coarse-grained molecular dynamics (CGMD) approach, we predicted the microstructure, dynamics, and rheology of hydrated smectite (more precisely, montmorillonite) clay gels containing up to 2,000 clay platelets on length scales up to 0.1 mu m. Simulations investigated the impact of simulation time, platelet diameters (6 to 25 nm), and the ratio of Na to Ca exchangeable cations on the assembly of tactoids (i.e., stacks of parallel clay platelets) and larger aggregates (i.e., assemblages of tactoids). We analyzed structural features including tactoid size and size distribution, basal spacing, counterion distribution in the electrical double layer, clay association modes, and the rheological properties of smectite gels. Findings: Our results demonstrate new potential to characterize and understand clay aggregation in dilute suspensions and gels on a scale of thousands of particles with explicit representation of counterion clouds and with accuracy approaching that of all-atom molecular dynamics (MD) simulations. For example, our simulations predict the strong impact of Na/Ca ratio on clay tactoid formation and the shear-thinning rheology of clay gels.

As a newly emerged solution for supporting the new generation of offshore wind turbines (OWTs), the pile-bucket foundation has received wide attention. However, little attention has been paid to the grouted connection that connects the monopile and bucket foundation. As the loadtransferring, yet vulnerable component, the fatigue mechanism of the grouted connection and its influence on the cyclic laterally-loaded response of OWT foundation are still not clear. In this study, a sophisticated three-dimensional (3D) finite element (FE) model of the pile-bucket foundation with grouted connection is constructed, which incorporates a hypoplastic clay model and the concrete damage plasticity (CDP) to consider the cyclic load effect on both soil and grout material. A modal analysis is first performed to verify the rationality of the proposed model. Then the influence of cyclic load frequency, load amplitude and stiffener arrangement on the accumulation of pile head displacement, stress distribution and crack development of the grouted connection is systematically analyzed. Results indicate that as load frequency approaches the eigen-frequency, the OWT structure tends to vibrate more intensively, leading to stress concentration and fatigue damage of the grouted material and rapid accumulation of the pile-head displacement. The influence of load amplitude on grout damage seems to be limited in the contact area in the simulated cases. Meanwhile, the installation of stiffeners slightly mitigates the pile head displacement accumulation, but also raises the risk of stress concentration and fatigue damage of the grouted connection. The numerical results reveal the load-transferring function and fatigue damage of the grouted connection, which could provide some reference for an optimized structure and dynamic design for the pile-bucket foundation under cyclic load.