An elastoplastic analysis scheme for the cylindrical cavity expansion in offshore islands unsaturated soils considering anisotropy is established. The hydraulic properties and anisotropy caused by stress of unsaturated soils are coupled in an elastoplastic constitutive matrix for unsaturated soil to obtain the governing equations for the cylindrical cavity expansion problem, with an analytical solution that utilizes the original hydro-mechanical state of the soil as the initial conditions. Through a comparative analysis with other analytical solutions, the effectiveness of the new solution is verified. Moreover, the swelling response of the cylindrical cavity expansion in unsaturated soils is examined by systematically analyzing different parameters of the surrounding soil. The findings reveal that the development and rate of anisotropy in normal consolidated soil and over-consolidated soil exert a significant impact on the soil's mechanical characteristics. Nevertheless, the alteration in the model constant h has little effect on the soil's mechanical characteristics. The analytical solution introduces anisotropy and broadens the expansion theory of unsaturated soils to yield a more comprehensive theoretical framework for the comprehensive analysis of offshore islands' unsaturated soils.

Soils are known to be inherently anisotropic, resulting in complex responses to loading. This paper aims to develop an elastoplastic solution for the undrained expansion of a cylindrical cavity in sands adopting a non-associated and anisotropic model, SANISAND. The rigorous derivation of the stress-strain state of the soil element is provided following a standardized solving procedure. The dilatancy and crushing of the soil are invoked in the three-dimensional cavity expansion solution by adopting the critical state soil mechanics and limiting compression curve, respectively. By combining this with a governing equation that considers the undrained condition, the stress-strain state of the surrounding soil around the cavity can be determined. A subroutine is then implemented into the ABAQUS FEM simulation to verify the solution. The solutions are also validated against those based on an isotropic model, and anisotropic sand is used to investigate the effects of the initial effective mean stress, at-rest coefficient of earth pressure, and overconsolidation ratio on the stress distribution, stress path, and boundary surfaces.