Application of Convergence Confinement Analysis to the study of preceding displacement of a squeezing rock tunnel

Summary This paper deals with a study on the geomechanical behavior of the rock mass around the tunnel. The Convergence Confinement Analysis is one of the most appropriate way to estimate this behavior. Elasto-plastic models proposed by Hoek and Brown, Egger and Kastner were adopted in this study.Measurements of the preceding displacement were performed in a road tunnel in Japan, the geology of which is Neogene squeezing mudstone. The measured data was compared with the results of Convergence Confinement Analysis.The preceding displacement was about 40% against the upper bench face converged displacement and 17% against the final converged displacement. Meanwhile, it increased drastically in lower bench excavation step.Back analysis was also executed to calculate the extension of the plastic zone of the rock mass, which was half of that estimated by Convergence Confinement Analysis in upper bench excavation step. But, in lower bench excavation step, the results of the two analyses were very similar. The support pressure and support stiffness were also calculated.     

Some numerical issues using element‐free Galerkin mesh‐less method for coupled hydro‐mechanical problems

A new formulation of the element‐free Galerkin (EFG) method is developed for solving coupled hydro‐mechanical problems. The numerical approach is based on solving the two governing partial differential equations of equilibrium and continuity of pore water simultaneously. Spatial variables in the weak form, i.e. displacement increment and pore water pressure increment, are discretized using the same EFG shape functions. An incremental constrained Galerkin weak form is used to create the discrete system equations and a fully implicit scheme is used for discretization in the time domain. Implementation of essential boundary conditions is based on a penalty method. Numerical stability of the developed formulation is examined in order to achieve appropriate accuracy of the EFG solution for coupled hydro‐mechanical problems. Examples are studied and compared with closed‐form or finite element method solutions to demonstrate the validity of the developed model and its capabilities. The results indicate that the EFG method is capable of handling coupled problems in saturated porous media and can predict well both the soil deformation and variation of pore water pressure over time. Some guidelines are proposed to guarantee the accuracy of the EFG solution for coupled hydro‐mechanical problems. Copyright © 2008 John Wiley & Sons, Ltd.     

Numerical modeling of combined effects of upward and downward propagation of shear bands on stability of slopes with sensitive clay

Modeling of progressive development of zones of large inelastic shear deformation (shear band) that results from strain‐softening behavior of sensitive clays could explain the failure mechanisms of large landslides. Because of toe erosion, a shear band can be initiated and propagated upward (inward) from the river bank. On the other hand, upslope surcharge loading could generate shear bands that might propagate down towards the river bank. In the present study, upward and downward propagation of shear bands and failure of sensitive clay slopes are modeled using the Coupled Eulerian Lagrangian approach in Abaqus finite element (FE) software. It is shown that the formation and propagation of shear bands are significantly influenced by kinematic constraints that change with displacements of the soil masses, and therefore the propagation of an existing shear band might be stopped and new shear bands could be formed. The main advantages of the present FE modeling are: (i) extremely large strains in the shear bands can be successfully simulated without numerical issues, (ii) a priori definition of shearing zones is not required to tackle severe strains; instead, the FE program automatically identifies the critical locations for shear band formation and propagation. Toe erosion could significantly increase the slope failure potential because of upslope surcharge loading. FE analyses with a thick and thin sensitive clay layers show that the global failure could occur at lower surcharge loads in the former as compared to the latter cases. Copyright © 2016 John Wiley & Sons, Ltd.     

Effects of soil fabric on behaviors of granular soils: Microscopic modeling

The effects of initial soil fabric on behaviors of granular soils are investigated by using Distinct Element Method (DEM) numerical simulation. Soil specimens are represented by an assembly of non-uniform sized spheres with different initial contact normal distributions. Isotropically consolidated triaxial compression loading and extension unloading in both undrained and drained conditions are simulated for vertically- and horizontally-sheared specimens. The numerical simulation results are compared qualitatively with the published experimental data and the effects of initial soil fabric on resulting soil behaviors are discussed, including the effects of specimen reconstitution methods, effects of large preshearing, and anisotropic characteristics in undrained and drained conditions. The effects of initial soil fabric and mode of shearing on the quasi-steady state line are also investigated. The numerical simulation results can systematically explain that the observed experimental behaviors of granular soils are due principally to their conditions of the initial soil fabric. This outcome provides insights into the observed phenomena in microscopic view.     

Use of microfluidic experiments to optimize MICP treatment protocols for effective strength enhancement of MICP-treated sandy soils

Microbially induced calcium carbonate (CaCO₃) precipitation (MICP) has been extensively studied for soil improvement in geotechnical engineering. The quantity and size of calcium carbonate crystals affect the strength of MICP-treated soil. In this study, microfluidic chip experiments and soil column experiments were conducted to optimize MICP treatment protocols for effective strength enhancement of MICP-treated sandy soils. The microscale experiments reveal that, due to Ostwald ripening, longer injection intervals allow crystals to dissolve and reprecipitate into larger crystals regardless of the concentration of cementation solution. Even though a cementation solution input rate of 0.042 mol/l/h is sufficient to maintain a high chemical transformation efficiency, a further reduction in the input rate by about four times resulted in an increase in the size of crystals produced by the end of treatment from about 40 to 60 μm. These findings were applied in soil column experiments. Results showed that significantly larger crystals and higher soil strength were achieved when the normalized rate of cementation solution injection was reduced from 0.042 to 0.021 mol/l/h. Crystal size and soil strength increased slightly more when the normalized input rate was further reduced from 0.021 to 0.010 mol/l/h. This study demonstrates how data from microscale microfluidic experiments that examine the effects of injection intervals and concentration of cementation solution on the properties of calcium carbonate crystals can be used to optimize MICP treatment in macroscale sand soil column experiments for effective strength enhancement.

     

The role of constitutive models in MPM simulations of granular column collapses

The granular column collapse is a well-established experiment which consists of having a vertical column of granular material on a flat surface and letting it collapse by gravity. Despite its simplicity in execution, the numerical modelling of a column collapse remains challenging. So far, much attention has been dedicated in assessing the ability of various numerical methods in modelling the large deformation and little to the role of the constitutive model on both the triggering mechanism and the flow behaviour. Furthermore, the influence of the initial density, and its associated dilatancy and strength characteristics, have never been included in the analyses. Most past numerical investigations had relied on simple constitutive relations which do not consider the softening behaviours. The aim of this study is to illustrate the influence of the constitutive model on the on-set of failure, the flow behaviour and the deposition profile using the material point method. Three constitutive models were used to simulate the collapse of two granular columns with different geometries and for two densities. The results of the simulations showed that the constitutive model had a twofold influence on the collapse behaviour. It defined the volume of the mobilised mass which spread along the flat surface and controlled the dissipation of its energy. The initial density was found to enhance the failure angle and flow behaviours and was more significant for small columns than for larger ones. The analysis of the potential energy of the mobilised mass explained the existence of two collapse regimes.     

包气带非水相污染物迁移离心模型研究中的相似比例关系

非水相污染物在包气带中的迁移是一个复杂的物理、化学及生物过程 ,离心模型是解决该问题的有效途径之一。采用理论分析和数值分析的方法 ,对两种湿润历史条件下 ,二相与三相体系中 ,毛细上升高度、达到毛细平衡所需要的时间以及迁移距离等与离心加速度之间的关系进行了研究 ,提出了存在于包气带非水相污染物迁移离心模型研究中的相似比例关系     

A two-scale constitutive framework for modelling localised deformation in saturated dilative hardening material

Undrained deformation of dilative sand generates negative excess pore pressure. It enhances the strength, which is called dilative hardening. This increased suction is not permanent. The heterogeneity at the grain scale triggers localisations causing local volume changes. The negative hydraulic gradient drives fluid into dilating shear zones. It loosens the soil and diminishes the shear strength. It is essential to understand the mechanism behind this internal drainage and to capture it numerically. The purpose of this paper is to develop a macroscopic constitutive relationship for the undrained deformation of saturated dense sand in the presence of a locally fully or partially drained shear band. Separate constitutive relations are generated for the band and intact material. Both time and scale dependence during pore fluid diffusion in saturated sand are captured, eliminating the mesh dependency for finite element implementations. The model is applied to the Gauss points that satisfy the bifurcation criterion. The proposed method is calibrated to recreate the undrained macroscopic response bestowed by an extra-small mesh. The microscopic behaviours inside and outside shear band predicted by this model are qualitatively in good agreement with individual material point behaviours inside and outside the shear band in the extra-small mesh. Depending on the loading rate and the shear band thickness, the response inside the band can be fully or partially drained, which governs the ultimate global strength. The calibrated model is exploited to simulate an upscaled biaxial compression test with semipermeable boundaries.