This study explores the effectiveness of soft viscoelastic biopolymer inclusions in mitigating cyclic liquefaction in loosely packed sands. This examination employs cyclic direct simple shear testing (CDSS) on loose sand treated with gelatin while varying the gelatin concentration and the cyclic stress ratio (CSR). The test results reveal that the inclusion of soft, viscoelastic gelatin significantly reduces shear strain and excess pore pressure during cyclic shear. Liquefaction potential, defined as the number of cycles to liquefaction (NL) at an excess pore pressure ratio (ru = Delta u/sigma ' vo) of 0.7, is substantially improved in gelatin-treated sands compared to gelatin-free sands. This improvement in liquefaction resistance is more pronounced as the inclusion stiffness increases. Furthermore, the viscoelastic pore-filling inclusion helps maintain skeletal stiffness during cyclic shearing, resulting in a higher shear modulus in gelatin-treated sand in both small and large-strain regimes. At a grain scale, pore-filling viscoelastic biopolymers provide structural support to the skeletal frame of a loosely packed sand. This pore filler mitigates volume contraction and helps maintain the effective stress of the soil structure, thereby reducing liquefaction potential under cyclic shearing. These findings underscore the potential of viscoelastic biopolymers as bio-grout agents to reduce liquefaction risk in loose sands.

Underground structures are subject to deterioration conditions in which water leakage occurs through cracks due to the long-term influence of soil and groundwater. Therefore, composite waterproofing sheets can play an important role in securing the leakage stability of structures by combining them with concrete structures. In this study, a total of eight composite waterproofing sheets were used according to the thickness of the compound and the properties of the material attached to the concrete, and the deformation characteristics at the bonding surface were identified through repeated tensile tests. Types A, B, and C, with a compound thickness of 1.35 to 1.85 mm and a single layer, had strong bonding performance, with a deformation rate of 0.5 to 2 x 10-4 and a DE/RE ratio of 0.3 to 1.3; tensile deformation progressed while maintaining integrity with the concrete at the bonding surface. Types D and E were viscoelastic and non-hardening compounds with a compound thickness of 1.35 to 3.5 mm, where the strain rate due to tensile deformation was the lowest, at 0.1 x 10-4 or less, and the DE/RE ratio was -5 to 3; therefore, when internal stress occurs, the high-viscosity compound absorbs it, and the material is judged to have low deformation characteristics. Types F, G, and H, which were 2 to 2.9 mm thick and had two layers using a core material, were found to have characteristics corresponding to tensile deformation, as the strain rate increased continuously from 0.2 to 0.5 x 10-4, and the DE/RE ratio increased up to 8 mm of tensile deformation.

This paper proposes a semi-analytical solution for one-dimensional consolidation of viscoelastic unsaturated soil considering a variable permeability coefficient under exponential loading. The governing equations of excess pore air pressure (EPAP) and excess pore water pressure (EPWP) were acquired by introducing the Merchant viscoelastic model. By employing Lee's correspondence principle and the Laplace transform, the solutions for EPAP and EPWP were derived under the boundary conditions of the permeable top surface and impermeable bottom surface. Crump's method was then used to execute the inverse Laplace transform, yielding a semi-analytical solution in the time domain. Through typical examples, the dissipation of EPAP and EPWP and the change of the average degree of consolidation over time under the influence of different elastic moduli, viscoelastic coefficients, and air-to-water permeability ratios were studied. The variation of the permeability coefficient and its influence on consolidation were also analyzed. The findings of this research show that the consolidation rate of viscoelastic unsaturated soil is slower than that of elastic unsaturated soil; however, an acceleration in the consolidation of the soil is observed when changes in the permeability coefficient are considered. These discoveries enhance our comprehension of the consolidation behaviors exhibited by viscoelastic unsaturated soil, thereby enriching the knowledge base on its consolidation traits.

Wave propagation in an ocean site is an essential research topic in various scientific fields, such as offshore geotechnical engineering, ocean seismology, and underwater acoustics. Previous studies have considered the seabed soil as elastic or poroelastic, ignoring the viscoelastic characteristics of its solid skeleton. Based on the fractional-derivative viscoelastic theory and the modified Biot theory, considering the flow-independent viscosity related to solid skeleton, this paper proposes a generalized viscoelastic wave equation for a fluid-saturated porous medium. The equation has a flexible mathematical form to describe soil rheological properties more accurately through fractional order. On this basis, the total wave field equation of an ocean site, modeled as the fluid-poroviscoelastic-solid media, is established. Then an analytical solution for wave propagation in an ocean site subjected to obliquely incident P and SV waves is obtained, and its degeneration and extension are studied. The proposed method is comprehensively validated through experiment, analytical, and numerical methods. Finally, a parameter analysis is performed to investigate the effects of water depth, seabed properties (including viscoelastic parameters, fractional order and permeability), and incident angle on the seismic response of a poroviscoelastic seabed.

A novel theoretical model is proposed to investigate the torsional response of a pile in fractional-order viscoelastic unsaturated transversely isotropic soil with imperfect contact. This model employs Biot's framework for three-phase porous media along with the theory of fractional derivatives. Unlike previous models that assume continuous displacement at the pile-soil interface, this study uses the Kelvin model to simulate relative slippage between pile-soil contact surfaces (imperfect contact). Incorporating fractional-order viscoelastic and transversely isotropic models to describe the stress-strain relationship, comprehensive dynamic governing equations are derived. Using the separation of variables method, inverse Fourier transform, and convolution theory, analytical solutions for the frequency domain response and semi-analytical solutions for the time domain response of the pile head under semi-sine pulse excitation are obtained. Using numerical examples, the effects of model parameters in the fractional-order viscoelastic constitutive model, pile-soil relative slip and continuity model, and soil anisotropy on the torsional complex impedance, twist angle, and torque are presented.

The pile foundation construction adjacent to an operational subway tunnel can induce the creep effects of the surrounding soil of the tunnel, resulting in the deformation of the existing tunnel lining and potentially compromising the safe operation of the tunnel. Therefore, the Mindlin solution and the generalized Kelvin viscoelasticity constitutive model were employed to establish the theoretical calculation model for the deformation of the adjacent subway tunnel caused by the pile construction. Then, the effect of pile construction on the deformation of adjacent tunnels under different pile-tunnel spacing was analyzed via three-dimensional numerical simulation and theoretical calculation methods and compared with the field monitoring data. The results showed that the theoretical and numerical data are in agreement with the field monitoring data. The theoretical model provides closer predictions to the field-measured values than the numerical simulation. As the distance between the pile and the tunnel increases, both the vertical settlement and the horizontal displacement of the subway tunnel lining exhibit a gradual reduction. In the hard plastic clay region of Hefei City (China), pile foundation construction near an operational subway tunnel can be classified into three distinct zones based on proximity to the tunnel: the high-impact zone (3.0 D). The pile foundation in high-, moderate-, and low-impact zones should be monitored for 7 days, 3 days, and 1 day, respectively, to ensure the stable deformation of the lining.

An analysis for the torsional dynamic response of end-bearing pile foundations embedded in a layered transversely isotropic geomaterial (soil/rock) is presented. The deformation of the transversely isotropic soil or rock is described by the method of separation of variables. The elasticity theory for a viscoelastic medium with frequency independent hysteretic material damping, and the Extended Hamilton's Principle are utilised to derive the differential equations describing pile and soil motions. The differential equations are solved analytically in an iterative algorithm. The accuracy of the analysis is verified with existing studies reported in the literature for pile foundations embedded in a homogeneous and layered soil deposit. The effect of the degree of anisotropy on the pile-soil response - dynamic pile-head stiffness, distribution of pile rotation and torque with depth, dimensionless soil displacement function for various values of pile slenderness and pile-soil stiffness ratios in a homogenous soil deposit is investigated. Design charts of static pile-head stiffness in a homogeneous soil deposit for a wide range of pile-soil stiffness and pile slenderness ratios, and degree of anisotropy are also reported. The effect of soil layering for a pile embedded in a two-layered soil deposit is also studied.

The behavior of soft soils distributed in coastal areas usually exhibits obvious time-dependent behavior after loading. To reasonably describe the stress-strain relationship of soft soils, this paper establishes a viscoelastic-viscoplastic small-strain constitutive model based on the component model and the hardening soil model with small-strain stiffness (HSS model). First, the Perzyna's viscoplastic flow rule and the modified Hardin-Drnevich model are introduced to derive a one-dimensional incremental Nishihara constitutive equation. Next, the flexibility coefficient matrix is utilized to extend the one-dimensional model to three-dimensional conditions. Then, by combining the HSS elastoplastic theory with the component model, the viscoelastic-viscoplastic small-strain constitutive model is subsequently established. To implement the proposed model for numerical analysis, the corresponding UMAT subroutine is developed using Fortran. After comparing the results of numerical simulations with those of existing literature, the reliability of the constitutive model and the program written in this paper is verified. Finally, numerical examples are designed to further analyze the effects of small-strain parameters and viscoelastic-viscoplastic parameters on the time-dependent behavior of soft soils.

Earthquakes contribute to the failure of anti-dip bedding rock slopes (ABRSs) in seismically active regions. The pseudo-static method is commonly employed to assess the ABRSs stability. However, simplifying seismic effects as static loads often underestimates rock slope stability. The development of a practical stability analysis approach for ABRSs, particularly in slope engineering design, is imperative. This study proposes a stability evaluation model for ABRSs, incorporating the viscoelastic properties of rock, to quantitatively assess the safety factor and failure surface under seismic conditions. The mathematical description of the pseudo-dynamic method, derived in this study, accounts for the viscoelastic properties of ABRSs and integrates the Hoek-Brown failure criterion with the Kelvin-Voigt stress-strain relationship of rocks. Furthermore, to address concurrent translation-rotation failure in ABRSs, upper bound limit analysis is utilized to quantify the safety factor. Through a comparison with existing literature, the proposed method considers the effect of harmonic vibration on the stability of ABRSs. The obtained safety factor is lower than that of the quasi-static method, with the resulting percentage change exceeding 5%. The critical failure surface demonstrates superior positional accuracy compared to the Aydan and Adhikary basal planes, with minimal error observed between the physical model test and the numerical simulation test. The parameter sensitivity analysis reveals that the inclination of ABRSs exhibits the highest sensitivity (Sk) value across the three levels of horizontal seismic coefficient (kh). The study aims to devise an expeditious calculation approach for assessing the stability of ABRSs during seismic events, intending to offer theoretical guidance for their stability analysis. (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/).

Soft soils exhibit significant time-dependent effects during long-term deformation. To precisely describe the long-term behavior of soft soils, it is necessary to employ elastoplastic theory and rheology principles for investigating the stress-strain relationship of the soils. In this paper, a super-subloading modified Cam-clay model is initially derived. Subsequently, by introducing the Kelvin model to describe the creep behavior of soils, and combining it with the modified Cam-clay model, an overconsolidated structural viscoelastic-elastoplastic model is further presented. After converting the equation into matrix form and programming it in Fortran, the proposed model is implemented by ABAQUS. Then, the accuracy of the developed model and program is verified through comparison with existing literature and experimental results. Finally, parametric analysis is conducted to explore the impact of viscoelasticity, structure, and overconsolidation on the responses of soft soils.