Finite-element modeling of a complex deep excavation in Shanghai

The excavation of the north square underground shopping center of Shanghai South Railway Station is a complex deep excavation using the top-down construction method. The excavation has a considerable size and is close to the operating Metro Lines. In order to predict the performance of the excavation more accurately, 3D finite-element analyses are conducted to simulate the construction of this complex excavation. The effects of the anisotropic soil stiffness, the adjacent excavation, and zone excavation on the wall deformation are investigated. It is shown that the numerical simulation with anisotropic soil stiffness yields a more reasonable prediction of the wall deflection than the case with isotropic soil stiffness. The deformation of the shared diaphragm wall between two excavations is influenced by the construction sequence of the two excavations. The zoned excavation can greatly reduce the diaphragm wall deformation. However, only the zoned excavation at the first excavation stage affects the deformation of the walls significantly. When the depth of the excavation increases, the zoned excavation has minor effect on the deformation of diaphragm walls.     

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

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

抛物线型屈服面人工冻土蠕变本构模型研究

利用自行研制的W3Z-200型冻土三轴压缩蠕变仪,对淮南地区深表土冻土进行了三轴蠕变试验。根据试验结果,提出以标准黏弹塑性模型为基础,并耦合温度自由度,应用抛物线型屈服准则对模型中的黏塑性项进行改进,建立了新的本构模型来描述高围压复杂应力状态下冻土的蠕变变形特征。基于ADINA有限元软件平台,开发了冻土本构子程序。对淮南某冻结井筒开挖过程进行了数值模拟,计算结果验证了冻土抛物线型黏弹塑本构模型的正确性和合理性。     

Rock Mass Stability Investigation Around Tunnels in an Underground Mine in USA

Stability and deformation of rock masses around tunnels in underground mines play significant roles on the safety and efficient exploitation of the ore body. Therefore, understanding of geomechanical behavior around underground excavations is important and necessary. In this study, a three-dimensional numerical model was built and stress analyses were performed by using 3DEC software for an underground mine in USA using the available information on stratigraphy, geological structures and mechanical properties of rock masses and discontinuities. Investigations were conducted to study the effect of the lateral stress ratio (K₀), material constitutive models, boundary conditions and rock support system on the stability of rock masses around the tunnels. Results of the stress, displacement, failure zone, accumulated plastic shear strain and post-failure cohesion distributions were obtained for these cases. Finally, comparisons of the deformation were made between the field deformation measurements and numerical simulations.     

进路分步开挖背景应力场演化规律研究

鲁中小官庄铁矿由于埋藏深、地压大、矿体倾角比较缓、矿岩破碎,是国内有名的难采矿山之一。进路开挖以后支护巷道破坏严重,其围岩变形为无收敛变形,也是国内支护较难的矿山之一。针对小官庄铁矿情况,为掌握其地压活动规律,应用东北大学岩石破裂过程分析系统（RFPA）,模拟其采用无底柱分段崩落法进路开挖过程,对进路开挖过程巷道围岩应力变化进行数值分析。发现随着进路开挖,导致进路出现片帮、底鼓、顶板下沉等现象,并且在矿岩接触带出现高应力集中,导致两进路之间的间柱破坏严重,并随着进路开挖应力逐步向新开挖两进路之间的间柱转移,开挖顺序造成边界矿体出现高应力集中,导致边界矿体难采。数值分析结果可以看出：进路的开挖顺序是导致小官庄铁矿进路围岩破坏的一个主要原因。     

人工冻土黏弹塑性本构参数反分析研究

通过大量的人工冻结黏土蠕变试验研究,获得在高偏应力下人工冻结黏土具有显著的蠕变特性：当应力水平较低时,人工冻土发生黏弹性蠕变变性;当应力水平较高时,人工冻土发生黏塑性蠕变变性。因此,以西原模型为基础,增加一个黏弹性元件项和用抛物线屈服函数代替塑性屈服项,构建了高应力下抛物线型黏弹塑性蠕变本构模型。改进的西原模型可以较好描述人工冻土剪切蠕变特性。基于小生境技术的遗传算法和人工冻土黏弹塑性模型有限元理论,在FEPG（有限元程序生成器）自动生成平台上研发了位移反分析程序。利用现场实测数据对深井冻结井筒开挖过程进行了位移反分析,获得了与试验相一致的冻土本构力学模型参数,数值分析结果与实测数据之间的误差在6%以内,反演结果验证了模型的可行性。冻结井筒开挖过程的位移反分析获得的人工冻土本构模型参数,能反映冻结壁施工过程和有效冻结范围内冻土本构参数平均值,因而具有广泛的代表性,对冻结壁应力场、位移场和稳定性预测具有重要工程意义。     

The use of two-part Hooke’s model (TPHM) to model the mine-by test at Mont Terri Site,Switzerland

The full-scale mine-by (MB) test conducted in 2008 in the Mont Terri underground rock laboratory (Switzerland) investigated the deformation and the coupled hydro-mechanical behavior of the Opalinus Clay in response to tunnel excavations. The Opalinus Clay is currently under investigation in Switzerland as a potential host rock for geologic disposal of high-level radioactive waste. To further improve the understanding and modeling of the coupled processes and their impact on the performance of a geologic repository in Opalinus Clay, a newly developed two-part Hooke’s model (TPHM) was implemented into a geomechanical simulator. A three-dimensional simulation model based on the TPHM was then developed to predict the deformation and pore pressure responses in the near field of the MB Niche 2 test at the Mont Terri Site. The usefulness and validity of the TPHM are demonstrated by the consistency between simulation results and field observations. Simulation results show that the pore pressure disturbance becomes visible at about 11 m ahead of the mine-by excavation advancing face (along the longitudinal direction of the MB Niche). The results also demonstrate that there exists a good correlation between the excavation damage zone (EDZ) and the pore-pressure evolution, which may have important practical implications for monitoring EDZ evolution with pore-pressure sensors. The simulation results, which are sensitive to the constitutive relationships used in the model, capture both the observed displacements and the size of the damage zone, whereas the approach based on the conventional Hooke’s law underestimates both. The comparison between simulated and observed results also indicates that laboratory-measured mechanical properties can be used to accurately predict field-scale mechanical deformations, as long as valid constitutive relationships are employed.     

Deep Cement Mixed Walls with Steel Inclusions for Excavation Support

Deep cement mixed (DCM) walls are widely used in supporting excavations in many parts of the world. In this paper, a case study of an excavation supported by a DCM wall with steel inclusions is analysed using a three-dimensional finite element model and based on the coupled theory of nonlinear porous media. The DCM wall is constructed with wide flange steel inclusions. The stress–strain behaviour of the DCM wall section is simulated using an extended version of the Mohr–Coulomb model, which considers the strain-softening behaviour of DCM columns beyond yield. The computed lateral deformations are compared with the field measurements to validate the numerical modelling procedure. Using the same case study, the internal stability of the wall against bending and shear failure modes is investigated. In addition, the lateral pressure distribution along the wall length is investigated because in practice design is carried out considering a uniform pressure distribution assuming rigid wall movements. A parametric study was carried out to investigate the viability of DCM walls in supporting excavations by varying the spacing between steel inclusions, wall thickness and initial lateral earth pressure. Based on the results of the parametric study, guidelines are proposed to select the most efficient geometric arrangement of steel inclusions within DCM walls.     

Coupled THM processes in EDZ of crystalline rocks using an elasto-plastic cellular automaton

This paper aims at a numerical study of coupled thermal, hydrological and mechanical processes in the excavation disturbed zones (EDZ) around nuclear waste emplacement drifts in fractured crystalline rocks. The study was conducted for two model domains close to an emplacement tunnel; (1) a near-field domain and (2) a smaller wall-block domain. Goodman element and weak element were used to represent the fractures in the rock mass and the rock matrix was represented as elasto-visco-plastic material. Mohr–Coulomb criterion and a non-associated plastic flow rule were adopted to consider the viscoplastic deformation in the EDZ. A relation between volumetric strain and permeability was established. Using a self-developed EPCA^2D code, the elastic, elasto-plastic and creep analyses to study the evolution of stress and deformations, as well as failure and permeability evolution in the EDZ were conducted. Results indicate a strong impact of fractures, plastic deformation and time effects on the behavior of EDZ especially the evolution of permeability around the drift.     

考虑表面形貌特征的岩体结构面蠕变特性

工程岩体开挖会遇到各种复杂应力条件,岩体在荷载作用下随时间会产生流变现象,这种“累进性破坏”导致工程岩体常出现滑坡、塌方、大变形和支护困难等问题.为解释天然岩体结构面的流变现象,以天然红砂岩结构面为研究对象,基于结构面三维形貌扫描试验、结构面分级加载剪切蠕变试验,对天然岩体结构面的剪切蠕变特性进行系统研究.研究发现,岩体结构面蠕变变形量与蠕变应力和结构面三维形貌特征指标正相关;岩体结构面蠕变曲线可分为3个阶段,即过渡蠕变、稳态蠕变和加速蠕变阶段,当法向应力一定时,结构面三维形貌特征指标越大,发生的蠕变破坏越剧烈;基于结构面剪切蠕变曲线与剪切蠕变速率曲线特征,建立了参数物理意义明确的岩体结构面剪切蠕变经验模型.     

重塑超固结上海软土力学特性及弹塑性模拟

对典型上海软土重塑样进行了围压不变和平均主应力不变的三轴排水剪切试验,得到重塑上海软土在不同初始超固结比和围压条件下的应力-应变关系,弄清了超固结比、围压以及应力路径对重塑上海软土的变形和强度特性的影响;根据土体的应力-应变曲线得到重塑上海软土的临界状态应力比及内摩擦角。采用姚仰平等建议的基于伏斯列夫面的超固结土本构模型,并根据等向压缩及三轴排水剪切试验确定其模型参数,对保持围压和平均主应力不变的三轴压缩试验进行了模型预测。预测结果表明,此超固结土本构模型能较好地反映重塑超固结上海软土的变形和强度特性。     

卸荷蠕变条件下软硬相接岩层非协调变形规律研究

软硬相接岩层在长时间跨度内的蠕变变形可能极为不协调,这对开挖工程安全稳定造成严重威胁。对某水电站引水隧洞同一开挖断面内砂质泥岩和泥质砂岩两种软硬相差较大的岩石进行卸荷条件下的蠕变试验和数值模拟,结果表明：在卸围压条件下,试样侧向蠕变变形较轴向发展更为迅速,呈现出显著的侧向扩容;初始围压水平较高时,试样卸荷蠕变破坏的侧向扩容效应更加显著;泥质砂岩在稳态蠕变阶段变形速率比砂质泥岩更慢,蠕变破坏时的应变值均比砂质泥岩蠕变破坏时的应变值要小,征兆性更弱。基于Burgers模型,引入损伤变量,建立带损伤的蠕变本构模型,模型曲线较好地描述了岩石破坏偏应力水平之前的线性蠕变特征和破坏偏应力水平下的非线性加速蠕变特征;软硬相接岩层蠕变数值模拟表明,蠕变变形初期,软硬相接点变形差别幅度不大,随着时间增长,在交接面处发生层间错动,对围岩整体稳定十分不利,工程实践中应予以足够关注。     

Efficiency of excavations with buttress walls in reducing the deflection of the diaphragm wall

Installation of buttress walls against diaphragm walls has been used as an alternative measure for the protection of adjacent buildings during excavation, but their mechanism in reducing movements has not yet been fully understood. This study performs three-dimensional finite element analyses of two excavation case histories, one in clay with T-shape buttress walls and another in dominant sand with rectangular buttress walls, to establish analysis model. Then, a series of parametric study were performed by varying soil types, types and length of buttress walls based on the above-mentioned excavations. Results show that the mechanism of buttress walls in reducing wall deflections mainly came from the frictional resistance between the side surface of buttress wall and adjacent soil rather than from the combined bending stiffness from diaphragm and buttress walls. The buttress wall with a length <2.0 m had a poor effect in reducing the wall deflection because the soil adjacent to the buttress wall had almost the same amount of movement as the buttress wall, causing the frictional resistance little mobilized. Since the frictional resistance of buttress walls in a deep excavation has fully been mobilized prior to the final excavation depth, the efficiency of buttress walls in reducing the wall deflection in a deep excavation was much less than that in a shallow excavation. Rectangular shape of buttress walls was of a better effect than T-shape in the shallow excavation because frictional resistance between buttress walls and adjacent soil played a major role in reducing the wall deflection rather than bearing resistance of the flange. When the excavation went deeper, the difference in reducing the wall deflection between the R-shape and T-shape became small.     

Investigation of the integrated retaining system to limit deformations induced by deep excavation

A series of three-dimensional finite element analyses of deep excavations with the integrated system between buttress walls and diaphragm walls was conducted to investigate the effect of the buttress wall intervals, treatments, locations, height, and thickness on limiting deformations induced by deep excavation. The integrated retaining system was formed by maintaining buttress walls when soil was excavated. The wall deflection control mechanism of the integrated retaining system mainly came from the combined stiffness between the buttress wall and the diaphragm wall. In addition, the ground settlement control mechanism came from the combined stiffness between the buttress wall and the diaphragm wall, and the frictional resistance between the buttress wall and the surrounding soil. For achieving 50% reduction in the wall deflection and the ground surface settlement, the length and intervals of buttress walls that were applied to the integrated retaining system were at least 4 and 8 m, respectively. When the deflection at the diaphragm wall head was well restrained, for example, by the floor slab, the position of the buttress wall head could be located at a depth the diaphragm wall starts to bulge out. In such a case, the performance between the full height and limited height of buttress walls was quite close. Furthermore, a new well-documented excavation project was analyzed to verify the performance of the integrated retaining system. Results showed that the integrated retaining system worked excellently if the joints between buttress walls and diaphragm walls were constructed properly.     

基于路基土蠕变效应的路基预应力损失模型研究

新型预应力路基技术可用于整治路基病害或强化既有线路基,研究其预应力损失规律对保障预应力路基工程的长期安全具有重要意义。基于弹性理论提出了预应力路基侧压力板与路基边坡接触面中心点后方水平路径上附加应力的计算方法,分析表明该特征路径上水平附加应力随距板−土接触面中心点距离的增加呈现良好的指数函数衰减关系。基于水平附加应力的指数函数扩散规律及路基土的蠕变行为,推导了路基土的蠕变变形计算公式。进一步基于蠕变变形与预应力钢筋回缩变形的协调性,建立了路基预应力的损失模型和计算方法,并结合FLAC3D数值仿真开展了对比分析。研究表明：钢筋预拉力损失的理论预测曲线与数值仿真曲线吻合良好,二者偏差不足5%,论证了路基预应力损失模型的有效性;钢筋预拉力可在锚固后的60 d内达到稳定,且预拉力的稳定值可达其初始值的85%～90%,说明新型预应力路基技术可为路基土提供相当可观的稳定附加围压,从而通过改善路基土的应力状态达到强化路基的目的。     

应力路径对地表变形特性影响的数值模拟

针对不同环境条件及工程因素对地表变形影响的不同特点,研究了不同应力状态下土的变形特性,得出其影响地表变形的基本规律。通过典型三轴试验对比分析了变形规律,由试验数据空间反演了各应力路径下的本构模型空间,编制相应的模块嵌入有限元程序,模拟了几种不同路径下地表沉降算例,得到了导致地表变形的几个基本规律：①不同的地表变形问题对应的应力路径千差万别,需要用不同的应力路径条件下的本构模型进行分析;②应力路径的影响是不可忽略的,甚至在诸多因素中起主导作用;③数值建模方法建立起来的模型通过权值系列的不同反映应力路径的影响,能较好应用于此类问题的分析;④CTC与TC路径极限体变相差近20 %;⑤CTC路径的剪缩区范围比CTE路径大,两种路径下剪缩体应变值最大相差近17.3 %;CTE路径较CTC路径的剪胀更为明显,剪胀体应变最大相对差值近98 %。     

A simple prediction model for wall deflection caused by braced excavation in clays

Deep excavations particularly in deep deposits of soft clay can cause excessive ground movements and result in damage to adjacent buildings. Extensive plane strain finite element analyses considering the small strain effect have been carried out to examine the wall deflections for excavations in soft clay deposits supported by retaining walls and bracing. The excavation geometry, soil strength and stiffness properties, and the wall stiffness were varied to study the wall deflection behavior. Based on these results, a simple Polynomial Regression (PR) model was developed for estimating the maximum wall deflection. Wall deflections computed by this method compare favorably with a number of field and published records.     

Site scale modeling of slow-moving landslides,a 3D viscoplastic finite element modeling approach

This paper presents an advanced 3D numerical methodology to reproduce the kinematics of slow active landslides, more precisely, to reproduce the nearly constant strain rate (secondary creep) and the acceleration/deceleration of the moving mass due to hydrological changes. For this purpose, finite element analyses are performed in a large area covering a long time-span (12 years), in order to exhibit different interacting slope movements. First, we perform a stability analysis using the shear strength reduction (SSR) technique with a Mohr-Coulomb failure criteria. It is done in order to compute factors of safety (FS) and to identify two different scenarios, the first one being stable (FS > 1) and the second one being unstable (FS < 1). In the studied test case, the Portalet landslide (Central Spanish Pyrenees), the first scenario corresponds to an initial stable configuration of the slope and the second one to an unstable excavated configuration. Second, taking the first scenario as an initial condition, a time-dependent analysis is performed using a coupled formulation to model solid skeleton and pore fluids interaction, and a simplified ground water model that takes into account daily rainfall intensity. In this case, a viscoplastic constitutive model based on Perzyna’s theory is applied to reproduce soil viscous behavior and the delayed creep deformation due to the excavation. The fluidity parameter is calibrated to reproduce displacements measured by the monitoring systems. Our results demonstrate that 3D analyses are preferable to 2D ones for reproducing in a more realistic way the slide behavior. After calibration, the proposed model is able to simulate successfully short- and medium-term predictions during stages of primary and secondary creep.     

南宁膨胀土非线性流变模型研究

对南宁非饱和膨胀土进行一系列的加载-卸载固结蠕变试验,得到了不同含水率的膨胀土在不同应力水平下应变-时间曲线和不同时刻的应力-应变等时曲线簇。通过对膨胀土各阶段蠕变变形的分析,得到了膨胀土的固结蠕变是包含有黏弹性变形和黏塑性变形成分的非线性蠕变。采用模型流变与经验流变理论相结合的方法,将膨胀土的蠕变变形分成线性和非线性两部分进行分析,蠕变柔量和线性黏弹性模量随时间减小,而线性黏塑性模量随时间增大。结合膨胀土的线性黏塑性模型和非线性黏塑性经验模型,建立能够描述膨胀土非线性流变的黏塑性本构方程。根据试验所得数据,通过非线性拟合方法得到各模型参数,用该模型进行数值分析得到的理论与蠕变试验结果对比非常吻合。研究成果为实际工程提供了流变变形计算分析的可靠依据。     

加筋砂土挡墙面板刚度效果的数值分析

为了对加筋砂土挡墙面板的刚度效果有一个合理、定量的把握,利用非线性硬化-软化弹塑性有限元方法对相关的一系列室内模型试验结果进行了系统全面的数值分析。有限元解析中所采用的砂土本构模型是以修正塑性功为基本状态量的弹塑性硬软化模型,该模型可以较为精确地模拟砂土的应力路径效应。结果表明,利用这种较高精度的有限元方法对加筋砂土挡墙的变形破坏进行解析,不仅能较好地模拟加筋砂土挡墙基础底面的平均压力与沉降之间的关系,同时也能较好地再现由于面板刚度变化对加筋砂土挡墙基础的承载力以及渐进性变形破坏的影响。随着面板刚度的增加,面板对剪切带的抑制作用将随之增加,具体表现在砂土围压σ3增大所带来的砂土强度σ1的提高,进而使得加筋砂土挡墙的峰值承载力也随之增大。