Water spewing and muck plugging often occur during earth pressure balance (EPB) shield machines tunnelling in water-rich sandy strata, even though the conventional foam has been employed to condition sandy soils. In this study, a novel thickened foaming agent suitable for EPB shield tunnelling in water-rich sandy strata is developed. In contrast to conventional foam-conditioned sands, the thickened foam-conditioned sand has a low permeability due to the consistent filling of soil pores with the thickened foam, and the initial permeability coefficient decreases by approximately two orders of magnitude. It also exhibits a suitable workability, which is attributed to the enhanced capability of the thickened foam to condition sandy soils. In addition, the effect of concentration on the stability of the foam is explained by the Gibbs-Marangoni effect, and conditioning mechanisms for the thickened foam on sands are discussed from the evolution of foam bubbles.

A hybrid material point/finite volume method for the numerical simulation of shallow water waves caused by large dynamic deformations in the bathymetry is presented. The proposed model consists of coupling the nonlinear shallow water equations for the water flow and a dynamic elastoplastic system for the seabed deformation. As a constitutive law, we consider a linear elastic-non-associative plastic model with the Drucker-Prager yield criterion allowing for large deformations under undrained cases. The transfer conditions between these models are achieved by using forces sampled from the hydraulic pressure and the friction terms along the interface between the seabed soil and shallow water. A detailed description regarding the coupled algorithm for the hybrid material point/finite volume method is presented. Several numerical examples are investigated to demonstrate the performance of the finite volume method for simulations of shallow water flow and the material point method for capturing the large deformation process of the solid phase. We also present numerical simulations of an undrained clay column collapse that induced shallow water waves and a dam-break problem to demonstrate the excellent performance of the proposed hybrid material point/finite volume method.

This paper presents a stability study on the collapse mechanisms of a plane-strain tunnel face in c-ϕ soils using the upper bound finite element method with rigid translatory moving elements (UBFELA-RTME) and nonlinear programming technique. Practical considerations are given to the unlined length influence behind the tunnel face. An advanced mesh adaptive updating strategy is adopted, aiming to improve the computational efficiency, the accuracy of upper-bound solutions, as well as the produced collapse mechanisms. The unlined length influence on the face stability and collapse mechanism of the tunnel face are determined with various combinations of tunnel depth ratios, soil friction angles, and dilatancy angles. Using the UBFELA-RTME with the Davis's approach and a mesh adapting strategy, the non-associated plasticity flow rule can be well approximated. The developed technique was validated against different numerical methods, and it is concluded that the tunnel face stability can be improved by increasing soil friction and dilatancy angles, and yet weakens as the unlined length increases where a mesh-liked collapse zone gradually appears on the tunnel vault top. It gradually evolves to a global collapse failure till the ground surface. The findings contribute to a better understanding of the ground surface failure under the unlined support length influence in tunnel construction.

This study presents the formulation and implementation of a fully implicit stabilised Material Point Method (MPM) for dynamic problems in two-phase porous media. In particular, the proposed method is built on a three-field formulation of the governing conservation laws, which uses solid displacement, pore pressure and fluid displacement as primary variables (u–p–U formulation). Stress oscillations associated with grid-crossing and pore pressure instabilities near the undrained/incompressible limit are mitigated by implementing enhanced shape functions according to the Generalised Interpolation Material Point (GIMP) method, as well as a patch recovery of pore pressures – from background nodes to material points – based on the same Moving Least Square Approximation (MLSA) approach investigated by Zheng et al. [1]. The accuracy and computational convenience of the proposed method are discussed with reference to several poroelastic verification examples, spanning different regimes of material deformation (small versus large) and dynamic motion (slow versus fast). The computational performance of the proposed method in combination with the PARDISO solver for the discrete linear system is also compared to explicit MPM modelling [1] in terms of accuracy, convergence rate, and computation time.

The Material Point Method (MPM) has been gaining increasing popularity as an appropriate approach to the solution of coupled hydro-mechanical problems involving large deformations. This study extends the implicit GIMP-patch method for coupled poroelastic problems recently proposed by Zheng et al. (2021b) to tackle large-deformation problems in (nearly) isochoric elastoplastic geomaterials, particularly by remedying the numerical inaccuracies caused by volumetric locking, such as spurious stress oscillations and an excessively stiff overall response of the system at hand. To overcome these difficulties in two-phase coupled analyses, the B¯ approach of Hughes (1980) is incorporated into an existing version of the implicit GIMP-patch method. Details regarding the formulation and implementation of the proposed method are provided, while several benchmark problems are numerically analysed to evaluate its performance in the presence of elastoplastic behaviour. Particular emphasis is placed on (i) mitigating effective stress oscillations and (ii) solving several two-phase, coupled, large deformation geotechnical problems. The numerical results confirm the suitability of the implicit B¯ GIMP-patch method for the solution of geotechnical problems spanning weak to strong hydro-mechanical coupling and small to large deformations.

This paper presents a single-point Material Point Method (MPM) for large deformation problems in two-phase porous media such as soils. Many MPM formulations are known to produce numerical oscillations and inaccuracies in the simulated results, largely due to numerical integration and stress recovery performed at non-ideal locations, cell crossing errors, and mass moving from one background grid cell to another. The same drawbacks lead to even worse consequences in the presence of an interstitial fluid phase, especially when undrained/incompressible conditions are approached. In this study, an explicit stabilised MPM, based on the Generalised Interpolation Material Point (GIMP) method with Selective Reduced Integration (SRI), is proposed to mitigate typical numerical oscillations in (nearly) incompressible coupled problems. It includes two additional features to improve stress and pore pressure recovery, namely (i) patch recovery of pore pressure increments based on a Moving Least Squares Approximation, and (ii) two-phase extension of the Composite Material Point Method for effective stress recovery. The combination of components leads to a new method named GC-SRI-patch. After a detailed description of the approach, its effectiveness is verified through analysing various consolidation problems, with emphasis on the representation of pore pressures in time and space.

The stability and collapse mechanism of tunnel faces are simplified conservatively to two-dimensional plane strain models along the longitudinal middle line of tunnel. Using the upper bound finite element method with rigid translatory moving element (UBFEM-RTME), a series of stability factors N_cr and collapse mechanisms displayed with active discontinuities are deduced. The influences of dimensionless buried depth ratio H/D, internal friction angle φ and dilatancy angle ψ on the variations of N_cr and mesh-like collapse mechanisms that are identical to the form of slip lines are discussed. A fitting formula of N_cr for the influence factors H/D and φ is deduced, and the effects of numbers and locations of active discontinuities are also investigated. This study illustrates that the UBFEM-RTME with combination of mesh adaptive updating strategies and reasonable and sufficient mesh density can improve the accuracy of the obtained N_cr values and the refinement of mesh-like collapse mechanism. The results reveal the main characteristics of the ultimate collapse mechanisms of tunnel faces, and they can provide theoretical supports for the stability evaluations of tunnel faces and pre-reinforcement scheme of soil strata.