敏感环境下基坑数值分析中土体本构模型的选择

数值分析已成为敏感环境下基坑工程分析的最重要手段,其关键是选择合适的土体本构模型和计算参数。在分析了岩土数值分析中常用土体本构模型特点的基础上,通过算例较系统地对比了各类模型在基坑开挖数值分析中的适用性。敏感环境下的基坑工程需重点关注墙后土体的变形,从满足工程需要和方便实用的角度出发,建议采用能考虑黏土的塑性和应变硬化特征、能区分加荷和卸荷且刚度依赖于应力水平的硬化类弹塑性模型,如MCC模型和HS模型进行分析。具体工程实例的分析,表明了硬化类弹塑性模型在敏感环境下基坑开挖数值分析中的适用性。     

A framework for automatic modeling of underground excavations and optimizing three‐dimensional boundary and finite element meshes derived from them—framework

Many problems in mining and civil engineering require using numerical stress analysis methods to repeatedly solve large models. Widespread acceptance of tunneling methods, such as New Austrian Tunneling Method, which depend heavily on numerical stress analysis tools and the fact that the effects of excavation at the face of a tunnel are distinctively three–dimensional (3D), necessitates the use of 3D numerical analysis for these problems. Stress analysis of a practical mining problem can be very lengthy, and the processing time can be measured in days or weeks at times. A framework is developed to facilitate efficient modeling of underground excavations and to create an optimal 3D mesh by reducing the number of surface and volume elements while keeping the result of stress analysis accurate enough at the region of interest, where a solution is sought. Fewer surface and volume elements mean fewer degrees of freedom in the numerical model, which directly translates into savings in computational time and resources. The mesh refinement algorithm is driven by a set of criteria that are functions of distance and visibility of points from the region of interest, and the framework can be easily extended by adding new types of criteria. This paper defines the framework, whereas a second companion paper will investigate its efficiency, accuracy and application to a number of practical mining problems. Copyright © 2012 John Wiley & Sons, Ltd.     

Prediction of dilation and permeability changes in rock salt

A numerical procedure has been developed for predicting dilation (porosity) and gas permeability changes in rock salt. The hierarchical single-surface constitutive model of Desai and co-workers is used a finite element program to calculate the state of stress and strain surrounding excavations in rock salt. The elastoplastic constitutive model accounts for strain hardening, a non-associative volumetric response and stress-path-dependent behaviour. The calculated stress and strain fields are used in a flow model based on the equivalent channel concept to predict permeability. Parameters for both the mechanical and permeability models are developed from laboratory test results. Two field experiments adjacent to underground excavations are modelled. The extent of the dilated rock zone around the excavation is predicted well, but the magnitude of the porosity and gas permeability is underpredicted very near the excavations. This discrepancy is attributed to model parameters derived from loading-only laboratory tests, whereas significant unloading occurs in the field. The shape of the yield surface was found to be an important factor in dilation and permeability predictions. Similar stress, strain and permeability fields were obtained with different model types (plane strain or axisymmetric) and initial stress states, and with instantaneous and progressive excavation.     

Developing constitutive models from EPR‐based self‐learning finite element analysis

A constitutive model that captures the material behavior under a wide range of loading conditions is essential for simulating complex boundary value problems. In recent years, some attempts have been made to develop constitutive models for finite element analysis using self‐learning simulation (SelfSim). Self‐learning simulation is an inverse analysis technique that extracts material behavior from some boundary measurements (eg, load and displacement). In the heart of the self‐learning framework is a neural network which is used to train and develop a constitutive model that represents the material behavior. It is generally known that neural networks suffer from a number of drawbacks. This paper utilizes evolutionary polynomial regression (EPR) in the framework of SelfSim within an automation process which is coded in Matlab environment. EPR is a hybrid data mining technique that uses a combination of a genetic algorithm and the least square method to search for mathematical equations to represent the behavior of a system. Two strategies of material modeling have been considered in the SelfSim‐based finite element analysis. These include a total stress‐strain strategy applied to analysis of a truss structure using synthetic measurement data and an incremental stress‐strain strategy applied to simulation of triaxial tests using experimental data. The results show that effective and accurate constitutive models can be developed from the proposed EPR‐based self‐learning finite element method. The EPR‐based self‐learning FEM can provide accurate predictions to engineering problems. The main advantages of using EPR over neural network are highlighted.     

Three‐dimensional inverse analyses of a deep excavation in Chicago clays

Numerical models are commonly used to estimate excavation‐induced ground movements. Two‐dimensional (2D) plain strain assumption is typically used for the simulation of deep excavations which might not be suitable for excavations where three‐dimensional (3D) effects dominate the ground response. This paper adapts an inverse analysis algorithm to learn soil behavior from field measurements using a 3D model representation of an excavation. The paper describes numerical issues related to this development including the generation of the 3D model mesh from laser scan images of the excavation. The inverse analysis to extract the soil behavior in 3D is presented. The model captures the measured wall deflections. Although settlements were not sufficiently measured, the predicted settlements around the excavation site reflected strong 3D effects and were consistent with empirical correlations. Copyright © 2010 John Wiley & Sons, Ltd.     

Assessment of strut forces for braced excavation in clays from numerical analysis and field measurements

One important consideration in the design of a braced excavation system is to ensure that the structural bracing system is designed both safely and economically. The forces acting on the struts are often determined using empirical methods such as the Apparent Pressure Diagram (APD) method developed by Peck (1969). Most of these empirical methods that were developed from either numerical analysis or field studies have been for excavations with flexible wall types such as sheetpile walls. There have been only limited studies on the excavation performance for stiffer wall systems such as diaphragm walls and bored piles. In this paper, both 2D and 3D finite element analyses were carried out to study the forces acting on the struts for braced excavations in clays, with focus on the performance for the stiffer wall systems. Subsequently, based on this numerical study as well as field measurements from a number of reported case histories, empirical charts have been proposed for determining strut loads for excavations in stiff wall systems.     

A Numerical Comparison of the Behavior of a Braced Excavation Using Two and Three-Dimensional Creep Plastic Analyses

Current study deals with investigating the effects of both time factor and the selection of a constitutive model type on predicting deformations of an excavation braced by nailing using two and three-dimensional finite element analyses. In addition, the effects of stress path and the type of defined initial conditions of the analytical model on deformations of the floor and walls of the excavation are also studied. Time factor, in the form of earth materials’ creep, can largely be entered into calculation of deformations of excavations by conducting viscoelastic and viscoplastic analyses. On the other hand, there hasn’t been done a comprehensive study regarding the creep behavior of excavations through comparing the results of two-dimensional and three-dimensional numerical analyses so far. The results showed that it’s largely possible to approach the actual deformation behavior of an excavation by considering the constitutive model of soft soil creep, SSC model, in the numerical plastic analyses. The effects of stress path on the deformation behavior of the excavation walls and excavation floor are investigated by using OCR stress ratio and POP stress difference; These two factors, both of which are also analogous, represent a boundary value for swelling behavior of the excavation floor and an increasing rate for the deformation behavior of the excavation walls since the increase in OCR or POP is equal to the increase in the soil lateral pressure coefficient at rest.     

A continuum mechanics based framework for optimizing boundary and finite element meshes associated with underground excavations—accuracy,efficiency and applications

A framework was developed to address the automatic optimization of the level of geometric detail required for stress analysis of underground excavations in mining, which was presented in the companion paper. The motivation for optimizing the mesh geometry stems from the over‐discretization of computational domain as the digital mine model is built while our knowledge of some of the input parameters is quite limited. Thus, the accuracy of the solution is not expected to be increased with a finely discretized mesh, only the computation time does. Therefore, it is acceptable if the results obtained from an optimized model have accuracy comparable to the uncertainty in input data (e.g. rock mass properties, geology, etc.). Although the mesh optimization framework automates the geometry optimization and reduces computation time, the accuracy of the solution from the resulting geometry must be evaluated to ensure the quality of the solution at the ‘region of interest’. Both a priori (mesh quality) and a posteriori (solution quality) measures are employed along with recording the mesh optimization time. Finally, the applicability of the mesh optimization framework is demonstrated by analysing a number of mining and civil engineering underground models. Copyright © 2005 John Wiley & Sons, Ltd.     

Modeling small-strain behavior of Taipei clays for finite element analysis of braced excavations

We examine the small-strain behavior of Taipei clays in braced excavation through a detailed analysis of a well-documented case history. Specifically, we analyze the case of the Taipei National Enterprise Center (TNEC) excavation using two soil models, the Modified Cam-clay model (MCC) and the three-Surface Kinematic Hardening model (3-SKH). Our finite element analysis includes a consideration of the over-consolidated stress state and the high initial shear modulus of the clay. Results show that the observed wall deflection and surface settlement can be satisfactorily predicted simultaneously using the 3-SKH model. This is an improvement on the MCC model, for which only wall deflection, not ground settlement, can be accurately predicted. This study re-confirms the importance of considering small-strain non-linear behavior for the over-consolidated stress state in finite element analyses of braced excavation responses.     

A multilaminate constitutive model accounting for anisotropic small strain stiffness

In this paper an extension of existing multilaminate soil models is presented, which can account for inherent and stress‐induced cross‐anisotropic elasticity in the small strain range and its dependency on the load history. In the multilaminate framework, material behaviour is formulated on a number of local planes in each stress point, and the macroscopic response of the material is obtained by integration of the local contributions. Strain‐induced anisotropy, which adds to the stiffness anisotropy inherently present in the material, is therefore intrinsically taken into account. Micro–macro relations between local parameters on plane level and global parameters on macroscopic level are obtained by the spectral decomposition of the global elastic compliance matrix. The model is implemented into a finite‐element code, and model predictions are compared with experimental data of triaxial tests on different soils involving small and large load cycles. The importance of cross‐anisotropic elasticity within the small strain range for predicting ground deformations in geotechnical boundary value problems is discussed at the example of an excavation problem. Copyright © 2012 John Wiley & Sons, Ltd.     

Variation with time of the factor of safety of slopes excavated in unsaturated soils

While the evolution with time of the factor of safety against slope failure is well understood for excavations performed under fully saturated conditions, this evolution has not been extensively studied for excavations carried out in unsaturated soils. The objective of this paper is to investigate the relationship between the factor of safety and time for an excavation performed in an unsaturated silty soil, employing the finite element method. A hypothetical boundary value problem is considered and two types of analyses performed; in the first one, unsaturated soil behaviour is modelled through appropriate constitutive and soil–water retention curve models, whereas in the second one full saturation is assumed. The effect that the saturated soil permeability, its variation with suction, the increase of apparent cohesion due to suction, the depth of the groundwater table and the hydraulic hysteresis have on the results of the unsaturated analysis is examined in a parametric study. The analyses results demonstrate that for unsaturated soils the factor of safety may increase with time, in contrast to what is commonly accepted to be the case in fully saturated soils. Furthermore, it is not possible to know in advance which one of the two types of analysis performed (i.e. the saturated or the unsaturated) will produce conservative results. It is, therefore, advisable to perform unsaturated analysis in geotechnical practice, when dealing with excavations in such soils.     

基坑垮塌的离散元模拟及冗余度分析

目前基于构件设计的基坑设计理论使得很多支护体系缺乏必要的冗余度,并导致了众多基坑事故。将离散元方法引入到基坑垮塌模拟和冗余度研究中,提出了基于局部破坏的冗余度分析模拟方法。进而以探究支撑端部连接对支护结构冗余度的影响为例,对两种连接情况的基坑的破坏过程进行了模拟,对其中的连续破坏现象进行了分析,并将模拟结果与实际工程进行了对比。结果表明,支护体系冗余度的提高可以有效地抵抗连续破坏的发生,防止基坑在局部损伤情况下,由一种形式的破坏引发其他形式的破坏,由局部破坏导致整体崩溃。另外,模拟结果也验证了离散元法在基坑破坏模拟中的适用性与可行性,为类似地下工程的破坏模拟和基坑冗余度的深入研究提供了参考。     

A computational framework to predict the water‐leakage pressure of segmental joints in underwater shield tunnels using an advanced finite element method

The construction of shield tunnels under riverbeds and seabeds has considerably increased over the past decades. Due to the ultra‐high water head, water leakage through tunnel joints is a major concern during a tunnel's service life. One practical solution to prevent groundwater penetration is to implement ethylene‐propylene‐diene‐monomer gaskets at the segmental joints. However, numerical simulation of fluid pressure penetration into rubber materials remains a challenging problem in computational mechanics. Severe mesh distortions can occur due to large deformation. Consequently, a convergent solution is difficult to achieve. This paper presents an Abaqus‐based numerical framework to solve the previously mentioned problem using the implicit finite element solver. The key aspects of this framework are twofold: (1) a remesh and re‐map algorithm to overcome the excessive mesh distortion, and (2) simulation of fluid penetration into the contact interface of the gaskets to reproduce the water‐leakage process at the tunnel joints. The proposed framework is first tested to simulate the gasket‐in‐groove mechanical behavior and is then validated using experimental data and the solution produced by an explicit finite element solver. The developed framework is then adopted to predict the water‐leakage pressure at gasketed tunnel joints to illustrate the practical applications. Finally, the numerical results are compared with experimental data to demonstrate the accuracy and robustness of the proposed method and confirm its superiority and effectiveness over existing methods. This novel method can be used by tunnel designers to analyze and estimate the waterproof behavior of gasketed joints in shield tunnels without performing extensive experimental testing works.     

基坑施工对临近运营地铁影响控制标准及关键技术

随着城市的快速发展,地铁车站周边开发的建筑越来越多,周边建筑施工对地铁车站将产生不可避免的不利影响。本文结合某大型建筑深基坑工程的设计、施工及监测数据,通过有限元计算,分析基坑开挖施工力学特征,研究基坑施工对临近地铁车站的变形影响。研究发现,基坑采用合理的地下连续墙、钻孔灌注桩等围护体系、适当的被动区加固方式、科学的分基坑开挖施工组织等措施,可以保证基坑开挖施工对临近车站的影响在可控制、可接受的范围内,为软土地区类似基坑的设计、施工提供参考经验。     

Assessment of basal stability for braced excavation systems using the finite element method

Conventional methods of predicting the basal stability of braced excavations are unable to take into consideration the stiffness of the retaining wall and the depth of penetration of the wall below the bottom of the excavation. A simple and improved procedure for predicting the stability of strutted excavations using the finite element method is presented. Detailed studies were carried out to assess the effects of the wall properties and soil geometry on the stability of the excavation.     

Modelling strain localization in granular materials using micropolar theory: mathematical formulations

It has been known that classical continuum mechanics laws fail to describe strain localization in granular materials due to the mathematical ill‐posedness and mesh dependency. Therefore, a non‐local theory with internal length scales is needed to overcome such problems. The micropolar and high‐order gradient theories can be considered as good examples to characterize the strain localization in granular materials. The fact that internal length scales are needed requires micromechanical models or laws; however, the classical constitutive models can be enhanced through the stress invariants to incorporate the Micropolar effects. In this paper, Lade's single hardening model is enhanced to account for the couple stress and Cosserat rotation and the internal length scales are incorporated accordingly. The enhanced Lade's model and its material properties are discussed in detail; then the finite element formulations in the Updated Lagrangian Frame (UL) are used. The finite element formulations were implemented into a user element subroutine for ABAQUS (UEL) and the solution method is discussed in the companion paper. The model was found to predict the strain localization in granular materials with low dependency on the finite element mesh size. The shear band was found to reflect on a certain angle when it hit a rigid boundary. Applications for the model on plane strain specimens tested in the laboratory are discussed in the companion paper. Copyright © 2006 John Wiley & Sons, Ltd.     

An approach to large-scale field stress determination

The construction of stable structures in rock masses requires knowledge of the in situ stresses at the scale of excavations. However, the measurements obtained by the conventional overcoring technique are related to a small scale (centimetres). To extrapolate them to the scales of interest to rock mechanics (from meters to kilometres) requires a large number of individual stress measurements, followed by statistical analysis to avoid a considerable scatter of the measured values. In this paper, a method is proposed based on (a) large-scale surface stress and modulus measurements using the cylindrical jack method complemented by a special measuring scheme and then (b) back analysis for a given excavation shape. The method allows the simultaneous reconstruction of the stress components at the scale of excavation. A numerical simulation for a cylindrical excavation in an isotropic rock mass demonstrates the high accuracy and robustness of the method. The presence of a fractured zone surrounding the excavation can hamper the stress reconstruction, hence special measures should be taken to conduct the measurements in competent rock.     

Comparison of two inverse analysis techniques for learning deep excavation response

Performance observation is a necessary part of the design and construction process in geotechnical engineering. For deep urban excavations, empirical and numerical methods are used to predict potential deformations and their impacts on surrounding structures. Two inverse analysis approaches are described and compared for an excavation project in downtown Chicago. The first approach is a parameter optimization approach based on genetic algorithm (GA). GA is a stochastic global search technique for optimizing an objective function with linear or non-linear constraints. The second approach, self-learning simulations (SelfSim), is an inverse analysis technique that combines finite element method, continuously evolving material models, and field measurements. The optimization based on genetic algorithm approach identifies material properties of an existing soil model, and SelfSim approach extracts the underlying soil behavior unconstrained by a specific assumption on soil constitutive behavior. The two inverse analysis approaches capture well lateral wall deflections and maximum surface settlements. The GA optimization approach tends to overpredict surface settlements at some distance from the excavation as it is constrained by a specific form of the material constitutive model (i.e. hardening soil model); while the surface settlements computed using SelfSim approach match the observed ones due to its ability to learn small strain non-linearity of soil implied in the measured settlements.     

软土隧道横向抗震分析的简化响应位移法

响应位移法能够反映软土隧道在地震荷载作用下的动力反应特性,是一种简便实用的隧道抗震设计简化分析方法。以方形截面隧道为例,采用地层－结构整体动力有限元方法,验证了响应位移法的准确性与合理性;为了简化响应位移法的计算过程,基于平面应变假定,采用弹性理论的复变函数方法推导了土弹簧刚度的解析表达式,并与有限元解进行了对比分析。将土弹簧刚度的解析表达式应用到响应位移法的计算中,同时采用地层剪应力及自由场地震反应位移的简化计算模式,将简化后的响应位移法和整体动力有限元法进行了对比分析,验证了简化方法的可行性。     

Accelerating strategies to the numerical simulation of large‐scale models for sequential excavation

In this paper, a novel combination of well‐established numerical procedures is explored in order to accelerate the simulation of sequential excavation. Usually, large‐scale models are used to represent these problems. Due to the high number of equations involved, the solver algorithm represents the critical aspect which makes the simulation very time consuming. The mutable nature of the excavation models makes this problem even more pronounced. To accomplish the representation of geometrical and mechanical aspects in an efficient and simple manner, the proposed solution employs the boundary element method with a multiple‐region strategy. Together with this representational system, a segmented storage scheme and a time‐ordered tracking of the changes form an adequate basis for the usage of fast updating methods instead of frontal solvers. The present development employs the Sherman–Morrison–Woodbury method to speed up the calculation due to sequential changes. The efficiency of the proposed framework is illustrated through the simulation of test examples of 2D and 3D models. Copyright © 2006 John Wiley & Sons, Ltd.