Probabilistic analyses of soil consolidation by prefabricated vertical drains for single‐drain and multi‐drain systems

Natural soils are one of the most inherently variables in the ground. Although the significance of inherent soil variability in relation to reliable predictions of consolidation rates of soil deposits has long been realized, there have been few studies that addressed the issue of soil variability for the problem of ground improvement by prefabricated vertical drains. Despite showing valuable insights into the impact of soil spatial variability on soil consolidation by prefabricated vertical drains, available stochastic works on this subject are based on a single‐drain (or unit cell) analyses. However, how the idealized unit cell solution can be a supplement to the complex multi‐drain systems for spatially variable soils has never been addressed in the literature. In this study, a rigorous stochastic finite elements modeling approach that allows the true nature of soil spatial variability to be considered in a reliable and quantifiable manner, both for the single‐drain and multi‐drain systems, is presented. The feasibility of performing an analysis based on the unit cell concept as compared with the multi‐drain analysis is assessed in a probabilistic context. It is shown that with proper input statistics representative of a particular domain of interest, both the single‐drain and multi‐drain analyses yield almost identical results. Copyright © 2016 John Wiley & Sons, Ltd.     

Numerical stability for modelling of dynamic two‐phase interaction

Dynamic two‐phase interaction of soil can be modelled by a displacement‐based, two‐phase formulation. The finite element method together with a semi‐implicit Euler–Cromer time‐stepping scheme renders a discrete equation that can be solved by recursion. By experience, it is found that the CFL stability condition for undrained wave propagation is not sufficient for the considered two‐phase formulation to be numerically stable at low values of permeability. Because the stability analysis of the two‐phase formulation is onerous, an analysis is performed on a simplified two‐phase formulation that is derived by assuming an incompressible pore fluid. The deformation of saturated porous media is now captured in a single, second‐order partial differential equation, where the energy dissipation associated with the flow of the fluid relative to the soil skeleton is represented by a damping term. The paper focuses on the different options to discretize the damping term and its effect on the stability criterion. Based on the eigenvalue analyses of a single element, it is observed that in addition to the CFL stability condition, the influence of the permeability must be included. This paper introduces a permeability‐dependent stability criterion. The findings are illustrated and validated with an example for the dynamic response of a sand deposit. Copyright © 2015 John Wiley & Sons, Ltd.     

An improved implicit numerical integration of a non‐associated,three‐invariant cap plasticity model with mixed isotropic–kinematic hardening for geomaterials

The Sandia GeoModel is a continuum elastoplastic constitutive model that captures many features of the mechanical response for geological materials over a wide range of porosities and strain rates. Among the specific features incorporated into the formulation are a smooth compression cap, isotropic/kinematic hardening, nonlinear pressure dependence, strength differential effect, and rate sensitivity. This study attempts to provide enhancements regarding computational tractability, domain of applicability, and robustness of the model. A new functional form is presented for the yield and plastic potential functions. This reformulation renders a more accurate, robust, and efficient model as it eliminates spurious solutions attributed to the original form. In addition, we achieve a high‐performance implementation, because the local iterative method is allowed to recast residual vectors with a uniform dimensionality. The model is also furnished with a smooth, elliptical tension cap to account for the tensile yielding. Moreover, an efficient algorithm is introduced, which decreases the computational cost by differentiating the updated shear yield surface from the cap surface based on the trial relative stress state. Finally, various numerical examples including a large‐scale boundary value problem are presented to demonstrate the fidelity of the modified model and to analyze its numerical performance. Copyright © 2015 John Wiley & Sons, Ltd.