A novel framework for nonlinear thermal elastic-viscoplastic (TEVP) constitutive relationships was proposed in this study, incorporating three distinct thermoplasticity mechanisms. These four TEVP formulations, combined with an existing TEVP constitutive equation presented in the companion paper, were integrated into a coupled consolidation and heat transfer (CHT) numerical model. The CHT model accounts for large strain, soil selfweight, creep strains, thermal-induced strains, the relative velocity of fluid and solid phases, varying hydraulic conductivity and compressibility during consolidation process, time-dependent loading, and heat transfer, including thermal conduction, thermo-mechanical dispersion, and advection. The performance of CHT model, incorporating different TEVP constitutive equations, was evaluated through comparing the simulation results with measurements from laboratory oedometer tests. Simulation results, including settlement, excess pore pressure and temperature profiles, showed good agreement with the experimental data. All four TEVP constitutive relationships produced identical results for the consolidation behavior of soil that in the oedometer tests. The TEVP constitutive equations may not have a significant effect on the heat transfer in soil layers because of the identical performance on simulating soil compression. The CHT model, incorporating the four TEVP constitutive equations, was then used to investigate the long-term consolidation and heat transfer behavior of a four layer soil stratum under seasonally cyclic thermal loading in a field test, with excellent agreement observed between simulated results and measured data.
A numerical model that accounts for fully coupled long-term large strain consolidation and heat transfer provides a more realistic analysis for various applications, including geothermal energy storage and extraction, buried power cables, waste disposal, groundwater tracers, and landfills. Despite its importance, existing models rarely simulate fully coupled large-strain long-term consolidation and heat transfer effectively. To address this research gap, this study presents a numerical model, called Consolidation and Heat Transfer (i.e., CHT), designed for one-dimensional (1D) coupled large-strain consolidation and heat transfer in layered soils, with the added capability to account for thermal creep. The model employs a piecewise-linear approach for the coupled long-term finite strain consolidation and heat transfer processes. The consolidation algorithm extends the functionality of the CS-EVP code by incorporating thermally induced strains. The heat transfer algorithm accounts for conduction, thermomechanical dispersion, and advection, assuming local thermal equilibrium between fluid and solid phases. Heat transfer is consistent with the spatial and temporal variation of void ratio and seepage velocity in the long-term consolidating layer. This paper details the development of the CHT model, presents verification checks against existing numerical solutions, and demonstrates its performance through several simulations. These simulations illustrate the effects of seepage velocity, thermal boundary conditions, and layered soil configurations on the coupled heat transfer and consolidation behavior of saturated compressible soils.
The use of horizontal drains assisted by vacuum loading is an effective method for speeding up the consolidation of dredged soil slurry. However, few studies developed models for the large strain consolidation of clayey slurry with prefabricated horizontal drains (PHDs) under self-weight and vacuum loading considering the effects of nonlinear compression and creep. This study introduces a PHD-assisted finite strain consolidation model considering nonlinear compression and limited creep by incorporating an improved elasto-viscoplastic constitutive equation. Firstly, the governing equations for the consolidation of very soft soil with PHDs were derived and solved by the finite-difference method. Subsequently, the proposed consolidation model was verified by comparing the calculations with the finite element solutions, a laboratory model test, and a field trial performed in Hong Kong. Good agreement with the numerical solutions and measured results indicates that the proposed model can capture the consolidation features with PHD combining staged filling and time-dependent vacuum loading. Then, the proposed model was used to estimate a self-weight consolidation test and field test in Japan to show the performance of the proposed model. Finally, parametric studies were conducted to explore the influence of nonlinear compression and creep on the consolidation of soft soil with PHDs.
Snow, characterized as a unique granular and low-density material, exhibits intricate behavior influenced by the proximity to its melting point and its three-phase composition. This composition entails a structured ice skeleton surrounded by voids filled with air and spread with liquid water. Mechanically, snow experiences dynamic transformations, including bonding/degradation between its grains, significant inelastic deformations, and a distinct rate sensitivity. Given snow's varied structures and mechanical strengths in natural settings, a comprehensive constitutive model is necessary. Our study introduces a pioneering formulation grounded on the modified Cam-Clay model, extended to finite strains. This formulation is further enriched by an implicit gradient damage modeling, creating a synergistic blend that offers a detailed representation of snow behavior. The versatility of the framework is emphasized through the careful calibration of damage parameters. Such calibration allows the model to adeptly capture the effects of diverse strain rates, particularly at high magnitudes, highlighting its adaptability in replicating snow's unique mechanical responses across various conditions. Upon calibration against established experimental benchmarks, the model demonstrates a suitable alignment with observed behavior, underscoring its potential as a comprehensive tool for understanding and modeling snow behavior with precision and depth.
Vertical drain assisted by vacuum and/or surcharge preloading is an effective method for improvement of soft ground with high water content. A large settlement will occur, and the water flow may deviate from the Darcy's law. The creep is also non-negligible in estimating the long-term settlement of such soft ground. To accurately predict the consolidation process, this study develops an axisymmetric finite strain consolidation model based on Barron's free-strain theory incorporating the creep, radial and vertical flows, non-Darcian flow law, and void ratio-dependent hydraulic conductivity during the consolidation process. First, to mathematically validate the model and highlight the new model's features, the existing model not considering the creep and the non-Darcy flow is also adopted as a reference for comparison based on a benchmark simulation. Then, Rowe cell tests involving non-Darcian flow are simulated by the new model to experimentally validate the predictive performance. Furthermore, the model is applied to simulate the consolidation process of a long-term monitoring embankment to examine the applicability of the model for engineering practice. All results demonstrate that the model is capable of accurately describing the consolidation of soft soils with vertical drains under combined loading with features of creep and non-Darcy flow.
Prefabricated vertical drains (PVDs) combined with vacuum and/or surcharge loading have been widely adopted to improve the strength of soft soils. Precise consolidation analysis is the theoretical basis for the design of preloading method with PVD. Current consolidation theories for layered soils with PVD seldom consider the influence of large strain, nonlinear creep, and self-weight loading simultaneously. This paper, thus, presents a finite strain elastic visco-plastic consolidation model, called RCS-EVP, for radial consolidation of layered soils with PVD. RCS-EVP is developed based on the piecewise-linear method. It takes into account nonlinear creep with limit creep strain, variable boundary conditions, anisotropy of soil hydraulic conductivity, and variable compressibility and hydraulic conductivity during the consolidation under self-weight, time-dependent surcharge and/or vacuum loading. The performance of RCS-EVP is evaluated by comparing with the results from finite element simulations and a laboratory physical model test. The variations of settlement and pore pressure of a soft soil ground improved by vacuum preloading with PVD are estimated using RCS-EVP. The results indicate that RCS-EVP provides good estimates of long-term consolidation of layered soils with PVD under both laboratory and in-situ conditions.