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.

Gravel-bearing sandstone reservoirs represent a significant type of reservoir in oil and gas exploration. Due to the difference of the spatial random distribution the content and the shape of the gravel particles, these reservoirs exhibit complex mechanical properties and failure modes. In this study, a numerical model of gravel-bearing sandstone was developed by using the Finite Element Method (FEM) and were verified by the actual indoor experimental data. The effect of the gravel particle sizes, gravel content, and gravel types on the compressive peak strength and microcrack evolution processes are further analyzed. The results reveal that cracks initiate within the sandstone matrix surrounding the gravel and propagate through the gravel with continued loading. The primary factors governing the stability of gravel-bearing sandstone are the gravel radius and content. The variation in gravel penetration rate is synchronized with the changes in peak strength. By embedding gravel particles of different shapes into the model, it is observed that the peak compressive strength of round gravel is comparable to that of elliptical gravel, with both exhibiting higher peak strengths than angular gravel. Regression models demonstrate that the tensile strength difference between the gravel and the sandstone matrix is a critical parameter influencing gravel penetration. Confining pressure has a relatively minor effect on the elastic modulus, while its impact on peak compressive strength is significantly more pronounced.

Waves can cause significant accumulation of pore water pressure and liquefaction in seabed soils, leading to instability of foundations of marine hydrokinetic devices (MHKs). Geostatic shear stresses (existing around foundations, within slopes, etc.) can substantially alter the rate of pore pressure buildup, further complicating the liquefaction susceptibility assessments. In this study, the development of wave-induced residual pore water pressure and liquefaction within sandy seabed slopes supporting MHK structures is evaluated. Unlike most earlier studies that excluded the impact of shear stress ratios (SSR) on the residual pore pressure response of sloping seabeds, asymmetrical cyclic loadings are considered herein for a range of SSRs. To obtain wave-induced loading in the seabed (and cyclic shear stress ratios, CSRs), the poroelasticity equations governing the seabed response, coupled with those for fluid and structure domains, are solved simultaneously. Utilizing an experimental model based on anisotropic cyclic triaxial test data that includes CSR and SSR impacts, an equation for the rate of pore pressure buildup is developed and added as a source term to the 2D consolidation equation. Numerical investigations were performed by developing finite element models in time domain. The models were calibrated using particle swarm optimization method and validated against wave flume experimental data. The results indicate that the consideration of static shear stresses has led to sudden rise in residual pore pressures followed by fast dissipations at early and late time steps, respectively, beneath the structure. The exclusion of SSR is shown to cause significant overestimation of pore pressure accumulations at late cycles, potentially causing significant overdesign of MHK foundations. The impact of proximity to the free drainage boundary, CSR amplitude, and loading frequency on the accumulation of residual pore pressure is illustrated. The residual liquefaction susceptibility of the seabed is shown to decline by increase of the seabed slope angle.

Seepage problems in half-space domains are crucial in hydrology, environmental, and civil engineering, involving groundwater flow, pollutant transport, and structural stability. Typical examples include seepage through dam foundations, coastal aquifers, and levees under seepage forces, requiring accurate numerical modeling. However, existing methods face challenges in handling complex geometries, heterogeneous media, and anisotropic properties, particularly in multi-domain half-spaces. This study addresses these challenges by extending the modified scaled boundary finite element method (SBFEM) and using this method to explore steady seepage problems in complex half-space domain. In the modified SBFEM framework, segmented straight lines or curves, parallel to the far-field infinite boundary, are introduced as scaling lines, with a one-dimensional discretization applied to them, thereby reducing computational costs.Then the weighted residual method is applied to obtain the modified SBFEM governing equations and boundary conditions of steady-state seepage problem according to the Laplace diffusion equation and Darcy's law. Furthermore, the steady seepage matrix at infinity is obtained by solving the eigenvalue problem of Schur decomposition and then the 4th-order Runge-Kutta algorithm is used to iteratively solve until the seepage matrix at the boundary lines is reached. Comparisons between the present numerical results and solutions available in the published work have been conducted to demonstrate the efficiency and accuracy of this method. At the same time, the influences of the geometric parameters and complex half-space domain on the seepage flow characteristics in complex half-space domain are investigated in detail.

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.

The soil moisture content (SMC) of moist clay directly affects the traction performance of off-road tire. This study set up a high-fidelity interaction model between off-road tire and moist clay with various moisture content, developed by coupling the finite element method (FEM) and smoothed particle hydrodynamics (SPH) algorithm. The interaction behavior between pneumatic tire and moist clay is studied. Firstly, a finite element model of tire which can characterize the complex structure and nonlinear mechanical properties is established. The Drucker-Prager (D-P) constitutive model parameters of clay with various moisture levels are calibrated by soil mechanical test. The moist clay with various moisture content is modeled through the SPH algorithm. The hybrid FEM-SPH interaction model is used to define the tire-moist clay interaction. Moreover, a traction performance test device suitable for tire-moist clay is developed to verify the accuracy of the interaction model. The influence of soil moisture content and tire operating conditions include vertical load and inflation pressure on the longitudinal traction coefficient, rolling resistance coefficient and instantaneous sinkage of tire center are quantitatively analyzed. The purpose of this study is to provide accurate tire force information under moist clay for unmanned ground vehicle (UGV), which can improve the problem of wheel instantaneous sinkage of tire center and slip under moist clay, and effectively reduce the yaw phenomenon in the path tracking process.

This paper presents a method for analyzing slope stability in anisotropic and heterogeneous clay using a strength reduction finite element method (SRFEM) integrated with the level set method (LSM). Anisotropy refers to the inherent anisotropy in the clay's strength, while heterogeneity describes the spatial variability in strength parameters. The static LSM uses a zero level set function to model heterogeneous clay slopes. The method is validated through undrained slope stability analyses on different types of anisotropic clay and heterogeneous fields, showing its effectiveness in modeling anisotropic shear strength and capturing the characteristics of heterogeneous regions. The results indicate that the proposed method accurately predicts factors of safety and slip surfaces across various soil conditions, accounting for both anisotropic and heterogeneous characteristics.

This study investigates the underlying causes of pier displacement and cracking in a highway link bridge. The initial geological assessment ruled out slope instability as a contributing factor to pier movement. Subsequently, a comprehensive analysis, integrating in situ soil investigation and finite element modeling, was conducted to evaluate the influence of additional fill loads on the piers. The findings reveal that the additional filled soil loads were the primary driver of pier tilting and lateral displacement, leading to a significant risk of cracking, particularly in the mid- of the piers. Following the removal of the filled soil, visual inspection of the piers confirmed the development of circumferential cracks on the columns of Pier 7, with the crack distribution closely aligning with the high-risk zones predicted by the finite element analysis. To address the observed damage and residual displacement, a reinforcement strategy combining column strengthening and alignment correction was proposed and validated through load-bearing capacity calculations. This study not only provides a scientific basis for analyzing the causes of accidents and bridge reinforcement but, more importantly, it provides a systematic method for analyzing the impact of additional filled soil loads on bridge piers, offering guidance for accident analysis and risk assessment in similar engineering projects.

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.

Piled raft foundations are increasingly used in construction due to their cost efficiency, requiring fewer piles than traditional pile foundations. Their ability to withstand cyclic lateral loads, such as those from earthquakes and wind forces, is crucial for structural stability. Understanding their response under cyclic loading conditions is essential, and finite element modeling (FEM) is a valuable tool for analyzing these behaviors. A recent 3D FEM study examined the performance of piled raft foundations in clay soils, focusing on loading pattern, frequency, and number of cycles. Results showed that lateral load capacity decreased as cycle count and frequency increased, with full cyclic loading (FCL) having a more pronounced effect than half cyclic loading (HCL). The raft shared 20.57-39.07 % (HCL) and 27.68-55.13 % (FCL) of lateral loads at frequencies of 0.1-10 Hz over 20 cycles. Additionally, locked-in moments increased by 21 %, and the degradation factor ranged from 65 to 80 % for HCL and 70-90 % for FCL. These findings provide valuable insights into pile-soil interaction and foundation stability under cyclic lateral loading, ensuring more effective design strategies for structures exposed to dynamic forces. Future research should explore long-term cyclic effects to further optimize foundation performance.