Semi-Analytical Solutions for Consolidation in Multi-Layered Unsaturated Silt with Depth-Dependent Initial Condition

This paper presents an analytical model for one-dimensional consolidation analysis of multi-layered unsaturated soils under depth-dependent initial conditions. The general solutions are derived explicitly using the Laplace transform. By combining these general solutions with interfacial continuity conditions between layers and the boundary conditions, the reduced-order system is solved via the Euler method to obtain analytical solutions in the Laplace domain. Numerical inversion of the Laplace transform is then performed using Crump’s method to yield the final analytical solutions in the time domain. The model incorporates initial conditions that account for both uniform and linear distributions of initial excess pore pressure within the soil stratum. The proposed solution is verified by reducing it to degenerated cases (e.g., uniform initial pressure) and comparing it with existing analytical solutions, showing excellent agreement. This confirms the model’s correctness and demonstrates its generalization to multi-layered systems with depth-dependent initial conditions. Focusing on a double-layered unsaturated soil system, the one-dimensional consolidation characteristics under depth-dependent initial conditions are investigated by varying the physical parameters of individual layers. The proposed solution can serve as a theoretical reference for the consolidation analysis of multi-layered unsaturated soils with depth-dependent initial conditions.