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.

The overconsolidation ratio (OCR) is a critical factor in determining the mechanical behaviour of overconsolidated clays. On the basis of the three requirements for the peak strength line, a continuous and smooth peak strength line is constructed from the perspective of the peak stress ratio, and then a new yield function for overconsolidated clays is developed. The developed yield function in the stress space is characterized by an elliptical curve. The evolution of the developed yield function in the stress space is captured by a new hardening parameter, which is constructed by integrating the proposed peak strength surface with the subloading surface concept. By combining the developed yield function with the non-orthogonal plastic flow rule, a non-orthogonal elastoplastic constitutive model of overconsolidated clays is established to consider the influence of the OCR on strength and deformation. The proposed model requires seven material parameters, all of which have a clear physical meaning and can be easily determined via conventional laboratory tests. Three typical stress paths are employed to demonstrate the essential features of the proposed model. The effectiveness of the proposed model is confirmed by comparing the experimental data with corresponding model predictions.

For many solids, irreversible deformation is often accompanied by changes in the internal structure, impacting the reversible responses, a phenomenon termed elasto-plastic coupling. This coupling has been observed experimentally in various geomaterials, including clayey and sandy soils, as well as hard and soft rocks. Fabric anisotropy, which characterizes the internal structure, is a distinct feature of soils and significantly influences both reversible and irreversible behaviors. In this study, we adopted a coupling formulation based on the framework of anisotropic critical state theory (ACST) to describe the anisotropic elasto-plastic coupling response of soils. The formulation incorporates a deviatoric fabric tensor F, which consistently quantifies the internal structure of soils in both reversible and irreversible range, into a hyperelastic formulation and a plastic model, respectively. A novel evolution rule of F, defined based on the current stress ratio and plastic strain, is proposed, where the direction gradually aligns with the loading direction and the norm achieves different asymptotic values depending on the applied loading paths. This allows for the representation of evolved anisotropy effects on elasticity, dilatancy and strength simultaneously, providing a natural description of elasto-plastic coupling. Within this coupling framework, any anisotropic model within ACST can serve as the plastic platform for developing the elasto-plastic coupling models with anisotropic hyperelasticity. Herein, a bounding surface plastic model is utilized for illustration. The proposed model's performance is demonstrated by especially comparing simulated results to test data on evolving elastic stiffness ratios and overall elastoplastic responses under varying monotonic and cyclic loading conditions.

The time-dependent behaviour of soft and clayey soils treated with Deep Cement Mixing (DCM) columns is important for analyzing the long-term performance of civil engineering infrastructures. Previous studies on DCMinstalled composite soil (CS) have primarily focused on examining the soil strength and stiffness characteristics. The limited focus on the time-dependent settlement and stress-strain distribution of CS underscores the need for a more comprehensive understanding of this complex phenomenon. In this study, a lab-scale physical ground model is designed and developed to investigate the time-dependent settlement profile of the composite Montmorillonitic Clay soil (MMC). The settlement behaviour of the ground model is assessed using Creep Hypothesis B and the results are further validated with the Power Law Model. Additionally, a FEM-based numerical simulation is performed to examine the time-dependent settlement and the stress distribution between the column and surrounding clay soil at different depths. The results from the physical model test show that the time-dependent parameter of the ground model (i.e., DCM column installed in MMC) is proportionate to the loading rate until the failure of the DCM column is reached. However, the time-dependent parameter was found to be decreased by 59.04 % in the post-failure phase of the DCM column. This reduction indicates that the DCM column was the primary load-bearing component before its failure. The numerical study shows that the pore water pressure dissipation in the clay soil and DCM column interface was similar at various depths. The top and bottom sections of the DCM column possess higher stress levels, which demonstrates its susceptibility for failure in the DCM column.

This paper presents a novel micropolar-based hypoplastic model to reproduce the stress-strain relationship of face mask chips-sand mixtures (MSMs) and their localized deformation. Based on a critical state hypoplastic model, a non-polar hypoplastic model for MSMs is first developed with modifications and new features: (1) the cohesion induced by face mask chips is considered by introducing an additional stress tensor into the Cauchy stress tensor; (2) the initial stiffness variation in MSMs is described with a modified tangential modulus; and (3) the effective skeleton void ratio concept is introduced to capture the initial and critical void ratio variations in MSMs. The model is then extended to its micropolar terms by incorporating the micropolar theory, which includes an internal length parameter and a couple stress induced by particle rotation, with the advantage of overcoming the mesh dependency problem in the conventional finite element method (FEM) based simulations. Moreover, the new micropolar hypoplastic formulations are implemented into a FEM code. The onset and evolution of shear bands in MSMs are investigated by simulating a series of biaxial tests on both pure sand and MSMs. Numerical results are also compared to experimental observations, demonstrating that the developed micropolar hypoplastic model can adeptly capture the shear band propagation in MSMs and their mechanical responses.

This article presents the findings of a comprehensive assessment of the predictive capabilities and limitations of advanced geotechnical numerical tools utilizing two sophisticated constitutive models for sands: the hardening soil model with small strains and hypoplasticity with intergranular strain. The evaluation is based on simulations of laboratory and centrifuge tests under monotonic and cyclic loading conditions. Initially, these models were calibrated and assessed using an experimental database on Fontainebleau sand. This database encompasses a range of laboratory results, including isotropic compression, drained monotonic triaxial, and undrained cyclic triaxial tests with varying initial conditions. The models, in general, provided good representation for monotonic experiments while some discrepancies were observed in undrained cyclic experiments. Subsequently, the calibrated models were employed to replicate a series of centrifuge tests involving a pile embedded in the same sand. The pile was subjected to various episodes of monotonic and cyclic lateral loading. In general, the models accurately replicated the experimental observations from tests conducted under monotonic loading conditions. Some small discrepancies were found in pile tests subjected to cyclic loading, these were however minor when compared to issues in predicting cyclic element tests at undrained conditions.

Accurately describing the solid-like and fluid-like behaviors of granular media is crucial in geotechnical engineering. While the unified frictional-collisional model, integrating rate-independent frictional and ratedependent collisional stresses, is widely used for solid-fluid phase transitions, an effective model is still under investigation, and comprehensive analyses are lacking. This study addresses these gaps by developing an enhanced elastoplasticity-based frictional-collisional model. The frictional stress is modeled using a critical-statebased elastoplasticity approach, and the collisional stress is formulated through an enhanced kinetic theory incorporating particle stiffness. Subsequently, comprehensive element simulations are conducted to explore the effects of concentration, particle stiffness, and strain rate paths on the model. The proposed model's effectiveness is also validated against experimental data. Finally, a detailed comparison with the typical mu(I) rheology model and a state-equation-based phase transition model is conducted. Our analyses show that the developed model effectively captures strain rate path and particle stiffness through the collisional stress component, while concentration-dependent characteristics are captured through both frictional and collisional stress components. Through comparative analyses, we also found that both the state-equation-based and elastoplasticity-based models depict solid-like behavior and replicate the rheology of granular media in a fluid-like state, similar to the mu(I) model. However, they differ in implementing critical state theory: the state-equation-based model acts as a partial-range phase transition model, describing stress evolution from the critical state to the fluid-like state, while the proposed elastoplasticity-based model serves as a full-range phase transition model, covering stress evolution from the initial to the fluid-like state.

Accurate continuum modelling of granular flows is essential for predicting geohazards such as flow-like landslides and debris flows. Achieving such precision necessitates both a robust constitutive model for granular media and a numerical solver capable of handling large deformations. In this work, a novel unified phase transition constitutive model for granular media is proposed that follows a generalized Maxwell framework. The stress is divided into an elastoplastic part and a viscous part. The former utilizes a critical-state-based elastoplasticity model, while the latter employs a strain acceleration-based mu(I) rheology model. Key characteristics such as nonlinear elasticity, nonlinear plastic hardening, stress dilatancy, and critical state concept are incorporated into the elastoplasticity model, and the non-Newtonian mu(I) rheology model considers strain rate and strain acceleration (i.e., a higher-order derivative of strain) to capture changes in accelerated and decelerated flow conditions. A series of element tests is simulated using the proposed unified phase transition model, demonstrating that the novel theory effectively describes the transition of granular media from solid-like to fluid-like states in a unified manner. The proposed unified model is then implemented within the material point method (MPM) framework to simulate 2D and 3D granular flows. The results show remarkable consistency with results from experiments and other numerical methods, demonstrating the model's accuracy in capturing solid-like behaviour during inception and deposition, as well as liquid-like behaviour during propagation.

Conventional plasticity assumes that a yield surface exists and the direction of plastic strain increment (DPSI) is uniquely dependent on the current stress state. Triaxial stress probing tests of yield and plastic flow of sand have been conducted using discrete-element modelling with polyhedral particles resembling the shapes of Toyoura sand. It is found that a yield surface does not exist, but a memory surface (MS) separating two types of distinct sand behaviour can be established. Within the MS, the DPSI is primarily controlled by the stress increment, and the magnitude of plastic strain increment is insensitive to the stress increment direction. When the stress state is on or outside the MS, a much larger plastic strain increment is observed if the stress increment points outside the MS, and the DPSI is dependent on both the current stress state and stress increment. The shape and size of the MS, which can be modelled by the SANISAND yield function, are dependent on the soil density and evolve with plastic strain.

The majority of existing effective stress-based constitutive models approach thermal effects through the temperature dependency of surface tension and its effects on the soil-water retention curve (SWRC) and effective stress. Experimental tests and theoretical studies, however, suggest that the temperature effect on surface tension alone is not sufficient to properly explain thermal-induced changes in the effective stress and SWRC. This study focuses on the temperature-dependent elastoplastic behavior of low plasticity unsaturated soils by developing a set of constitutive-level relations that incorporate temperature-dependent SWRC and effective stress models. These models account for the effect of temperature on the enthalpy, contact angle, and surface tension. The application of the presented constitutive relations was demonstrated and validated for low plasticity soils, specifically incorporating temperature effects into the hardening modulus, specific volume change, yield stress of the modified Cam-Clay model, and stress-strain relationships. The proposed relationships are incorporated in any effective stress-based constitutive model for modeling temperature dependency of elastoplastic response in low plasticity unsaturated soils. Employing these relationships can enhance the numerical simulation of low plasticity unsaturated soils under thermo-mechanical or other coupled processes involving temperature-dependent conditions.