A continuum finite element approach for rock failure and rubble formation

A continuum finite element formulation that models the failure of a rock mass, the collapse of this rubble into a void space and the associated bulking of the rubble is presented. This formulation is incorporated into a finite element computer program that treats the large deformation non-linear motion of axisymmetric and two-dimensional continua. A dynamic relaxation procedure is used for the solution algorithm. Example calculations are given that demonstrate the formulation.     

基于节理网络有限元的柱状节理岩体等效变形模量结构效应研究

因其特殊的地质成因,柱状节理成为一种具有强烈非连续性及各向异性的特殊岩体结构。以白鹤滩水电工程设计、建设中所遇到的柱状节理玄武岩体为背景,针对柱状节理岩体的等效变形模量问题,在探讨了其结构效应表征参数的基础上,利用节理网络有限元为工具,研究各结构效应表征参数对柱状节理岩体等效变形模量的影响。研究结果表明：在柱轴线横向平面上,对于柱体不规则程度,当柱体完全不规则相对于完全规则情况,等效模量提高大约10%左右;在其他条件不变情况下,当柱体平均边长从0.1 m增加至0.5 m时,等效模量从5.36 GPa增加至23.4 GPa;对不同加载方向的研究结果表明,柱状节理在柱轴横向平面上可以视作各向同性;当节理刚度线性递增时,相应的岩体等效模量也基本符合线性递增的规律。在平行柱体轴线平面,节理组2间距越大,柱状节理岩体的等效变形模量越大;随着错距比的增加,柱状节理岩体的等效变形模量呈现先增加后减少的趋势,当错距比为50%时,节理岩体的等效变形模量取最大值;对于节理刚度的影响,与柱轴横向平面类似的,随柱体轴线平面上节理刚度的线性增加,相应的岩体模量也基本符合线性递增的规律。该研究成果与已有的现场研究成果及三维数值模型成果相比,吻合良好。说明了文中提供的研究结论具有较好的可信度及参考价值。     

软弱围岩隧道施工全过程非线性有限元分析

结合赣州-龙岩铁路某单线隧道软弱围岩地质情况,采用非线性有限元法对软弱围岩条件下的铁路隧道湿喷纤维混凝土支护结构施工过程进行了数值模拟,分析了围岩和支护结构体的非线性力学行为的应力场分布、位移场,围岩塑性区分布特征,指出台阶步开挖时拱顶下沉、底板上鼓、墙腰收缩的主要控制因素.     

A finite element formulation for brine transport in rock salt

A set of coupled field equations is developed for transport of liquid brine in natural rock salt. The natural rock salt consists of individual crystals brought together so that only a portion of the crystal faces or grain boundaries contribute to the hydraulically connected pore space. Transport of brine inclusions within individual crystals is considered to be thermally driven; whereas transport along crystal interfaces or grain boundaries is considered to be pressure driven. The field equations for both transport mechanisms are developed and incorporated in a finite element program. An analytical solution to a one-dimensional boundary value problem is derived and compared to the finite element solution. An application of the finite element code to radioactive waste emplacement is briefly discussed.     

一种新的节理裂隙岩体弹塑性模型

提出一种节理裂隙岩体弹塑性模型,将结构划分为块体单元和缝单元,其中缝单元可以是实际节理裂隙,也可以是人为缝单元。以块体单元形心的刚体位移和块体的平均应变作为基本未知量。该模型能充分考虑节理裂隙材料和块体材料的本构关系,计算量不大,是一种位移不协调单元。该模型特别适合节理裂隙岩体的数值分析。     

危岩主控结构面变形破坏的扩展有限元法模拟分析

危岩是三峡库区典型的地质灾害类型之一,主要分为坠落式危岩和倾倒式危岩两种失稳形式,危岩主控结构面受荷断裂扩展是其发育成灾的关键核心。将危岩主控结构面类比为宏观裂纹,运用扩展有限元方法（XFEM）,求解结构面扩展过程的移动非连续问题,探索在荷载作用下危岩主控结构面的断裂扩展行为,以重庆万州太白岩危岩为例,考虑岩石的抗拉强度、主控结构面的几何位置与倾角对危岩的变形破坏模式与稳定性的影响。结果表明,随着岩石抗拉强度的降低,危岩的稳定性程度降低;随着主控结构面倾角的减小,危岩的稳定性程度提高,同时结构面开裂路径也明显改变;危岩体两种失稳模式的破坏机制都为拉剪机制,相对于倾倒式危岩,同等条件下坠落式危岩的稳定性更差。     

Distinct element simulation of ultimate bearing capacity in jointed rock foundations

In this paper, distinct element method numerical modeling is applied to evaluate bearing capacity of strip footing rested on anisotropic discontinuous rock mass. As yet, a little work has been carried out to investigate the effect of joint set orientation on the bearing capacity of rock mass. Generally, the overall behavior of rock mass under loading is very complicated and such analysis should include deformation determination, sliding along discontinuities and failure of rock material. Failure mechanism of rock mass depended on both geometrical parameters of joint sets and strength parameters of rock mass. In this research, it is assumed that rock mass contains one joint set, and therefore the anisotropy in bearing capacity and rock behavior is only due to the existence and orientation of the joint set. In this study, by assuming constant strength parameters and using Mohr–Coulomb failure criterion for the single joint set and nonlinear Hoek–Brown failure criterion for rock material, variation of the bearing capacity values and the type of failure mechanism of rock mass with different joint set dips is investigated. The obtained results indicate that the ultimate bearing capacity of rock mass containing one joint set varies between 27 and 86 % of intact rock.     

考虑剪胀角的节理化岩体隧道开挖进尺优化研究

合理进尺的选取是影响隧道开挖施工安全和效率的关键问题之一。以杨家坪隧道现场监测数据为根据,通过ABAQUS软件开展数值试算确定了最接近实际隧道围岩的剪胀角,建立了考虑剪胀角影响的Mohr-Coulomb模型的隧道开挖数值模型,并按1.5、2.0、2.5、3.0、3.5 m开挖进尺开展了计算,对开挖过程中隧道围岩的变形和力学响应规律进行了对比研究。结果显示,考虑剪胀角影响的数值模型计算结果的总体规律与实测结果一致性较好;有锚杆支护的情况下开挖进尺选择为3.0 m能较好地兼顾施工安全与效率,应注意控制距开挖面5.0 m范围内的拱顶变形。     

Micromechanics approach to poroelastic behavior of a jointed rock

The formulation of macroscopic poroelastic behavior of a jointed rock is investigated within the framework of a micro–macro approach. The joints are modeled as interfaces, and their behavior is modeled by means of generalized poroelastic state equations. Starting from Hill's lemma extended for a jointed medium and extending the concept of strain concentration to relate the joint displacement jump to macroscopic strain, the overall poroelastic constitutive equations for the jointed rock are formulated. The analysis emphasizes the main differences and similarities of the resulting behavior with respect to that characterizing ordinary porous media. It is shown that, unlike ordinary porous media, conditions on the poroelastic parameters of joints are required for the macroscopic drained stiffness to entirely define the poroelastic behavior. This is achieved, for instance, if the joint network is characterized by a unique Biot coefficient. Extension of the analysis to non‐linear poroelasticity is also outlined. Finally, the theoretical formulation is applied to two particular cases of jointed rock for which explicit expressions of the overall poroelastic parameters are derived. Copyright © 2011 John Wiley & Sons, Ltd.     

A mixed finite element approach for viscoelastic wave propagation

In this paper, we are interested in the modeling of wave propagation in viscoelastic media. We present a family of models which generalize the Zeners model. We achieve its mathematical analysis: existence and uniqueness of solutions, energy decay and propagation with finite speed. For the numerical resolution, we extend a mixed finite element method proposed in [8]. This method combines mass lumping with a centered explicit scheme for time discretization. For the resulting scheme, we prove a discrete energy decay result and provide a sufficient stability condition. For the numerical simulation in open domains we adapt the perfectly matched layers techniques to viscoelastic waves [23]. Various numerical results are presented.     

基于离散元法的节理岩体边坡稳定性分析

节理岩体边坡的稳定性在很大程度上取决于节理的强度及其分布形式。由于节理岩体边坡的失稳破坏具有大变形和非连续的特点,因此,离散单元法成为研究节理岩体边坡破坏机制的最有效方法之一。通过采用离散元软件PFC2D进行数值模拟,对完整岩石及节理的力学性能进行研究,并建立含密集节理的岩体边坡模型,讨论了节理连通率对边坡破坏形式的影响。结果表明,节理岩体边坡的失稳破坏是一个渐进的过程;在多组节理密集分布的岩体边坡中,连通率越大,其稳定性越差;随着连通率的减小,边坡的破坏形式由大范围的滑坡转变为局部崩塌的形式     

节理岩体的三维阶谱复合单元法初步研究

阐述了节理岩体复合单元法基本原理,引入了阶谱有限单元法中的虚节点和广义节点概念,提出了覆盖虚棱（面）和广义棱（面）的概念,并进一步将虚节点和广义节点分类。在此基础上构造了节理岩体的三维复合单元不同阶次时的位移函数并推导了对应的刚度矩阵。算例表明此法是合理有效的。提出的概念有助于复合单元法和传统有限单元法二者升阶谱的理解与统一实现,并为实现节理岩体三维复合单元法的P型自适应奠定了基础。     

Discrete element modelling of an underground explosion in a jointed rock mass

The purpose of this paper is to investigate the application of the discrete element program UDEC for modelling underground explosions in a jointed rock mass. A field underground explosion test has been conducted and UDEC and AUTODYN were utilised in a coupled manner to predict the rock mass response due to the underground explosion. The dynamic equation for the state of the rock material obtained from impact tests was incorporated into the calibration modelling. The comparison of modelling results with test results and empirical formulas shows that UDEC is capable of modelling explosion wave propagation in jointed rock mass with high reliability.     

隧道围岩结构地震动稳定性分析的动力有限元强度折减法

为了得到隧道围岩结构的地震动安全系数,借助通用有限元软件ANSYS,首先对水平地震作用下的模型进行模态分析,得到质量阻尼系数和刚度阻尼系数;其次由静力分析模型得到竖向边界上的水平向支座反力,然后将结构自重转化为温度边界条件,通过热分析得到模型各节点的温度,从而实现在动力分析中考虑重力的影响;最后采用悬臂梁动力分析模型,导入热分析获得的模型各节点的温度,并在竖向边界上施加水平向支座反力,通过不断折减围岩塑性区的凝聚力c和内摩擦角?,直到计算不收敛为止,从而得到隧道围岩结构的地震动安全系数。数值算例结果表明：采用的方法是可行的,将围岩结构自重转化为节点温度的措施解决了以往动力分析不能考虑结构自重的难点,进而为以后地震作用下隧道动力安全系数的计算及其工程应用提供了理论依据。     

Strain parameters of discontinuities in rock for finite element calculation

The bulk strain induced by excavation of rooms for storage of highly radioactive waste is usually unimportant, but the change in aperture and shear displacement of permeable fractures may be of practical importance. Hence, the normal and shear of such discontinuities have to be considered and are commonly predicted by using special deformation moduli termed “normal” and “shear” stiffnesses in numerical calculations. For the use of some finite element methods, these have to be converted to compression and shear moduli ,which, in turn, requires that the thickness of the discontinuity is known. The normal and shear strain can be expressed in terms of soil mechanical parameters, yielding the compression modulus M and the shear modulus G, with the required form for finite element calculations. By definition they are functions of the normal and shear stiffness K_n and K_s. It is concluded that calculation of the normal and shear strain of discontinuities with clastic fillings by the use of finite element methods and deformation moduli derived from stiffness numbers is very uncertain except when the geometry of the weaknesses can be accurately defined.     

黏弹性节理岩体的位移不连续模型研究

选用三参数标准线性固体作为岩石本构,提出了一种考虑岩体黏弹性的位移不连续模型;根据一维黏弹性波的特征线法,推导了节理处质点速度、应力和应变递推公式。首先,基于分离式霍普金森压杆（SHPB）对砂层进行试验,得到其应力-应变关系,并换算出砂层节理的法向刚度;接着,通过一维强间断黏弹性波的波速公式、高频波衰减系数以及任一频率下的衰减系数,确定出数值算法中的三参数。最后,基于自制的摆锤装置,探讨了一维应力波在节理岩体中的传播规律,试验中以两根长1 000 mm、直径为68.50 m的岩杆作为入射和透射杆,以3 mm砂层模拟节理。试验和数值结果吻合度良好,进一步验证了该方法的可靠性。     

裂隙岩体流固耦合双重介质模型的有限元计算

基于孔隙-裂隙岩体的双重孔隙介质流固耦合计算的微分方程,利用伽辽金有限元法提出的相应有限元公式,并基于岩体分类指标（RQD,RMR）提出了与岩体应力状态相关的渗透系数计算公式。编制了相应的有限元程序并给出了应用算例,将计算结果与相关文献作了比较,得出相关结论。     

Elastoplastic analysis of jointed rocks using a coupled finite element and boundary element method

This paper presents a coupled, elastoplastic, finite element and boundary element method for the two-dimensional, non-linear analysis of anisotropic jointed rock. The non-linear and anisotropic behaviour of a jointed rock mass is simulated by representing the mass as an equivalent anisotropic, elastoplastic continuum, so that the influence of the jointing system is ‘smeared’ across the continuum, i.e. the individual joints are not modelled as discrete entities. Numerical examples have been solved to verify the capability, accuracy and efficiency of the present technique. The proposed technique has also been applied to the analysis of tunnel excavation problems in plane strain. The effects of anisotropy and non-linearity of the jointed rock mass during excavation have been investigated in some detail.     

充水节理岩体中大断面隧洞开挖的数值模型研究

采用非饱和的渗流应力耦合模型分析了糯扎渡水电站2#导流洞的开挖过程,导流洞部分洞身穿过节理带,洞顶水头较高。详细地讨论了高水头和节理带对施工过程的影响。节理带采用描述高密度平行节理组的各向异性节理本构来逼近;考虑开挖引起的介质变形对渗透系数的影响;考虑排水引起的饱和度变化对渗透系数的影响。所用的非饱和瞬态耦合模型可以模拟出开挖引起的EDZ区域孔隙水压力急剧升高、有效应力减小、渗透系数动态的变化以及排水对洞室稳定性的提高。数值模拟的计算结果与国外类似试验的一般性观测结论相吻合,因此,可以用来评价水位以下隧洞施工方法和施工速度的合理性和经济性。     

A 2D coupled finite element and boundary element scheme to simulate the elastic behaviour of jointed rocks

In this paper a coupled finite and boundary element formulation is developed for the analysis of excavation in jointed rock. The presence of joints in the rock mass has been included implicitly by treating it as an appropriate anisotropic elastic continuum. The boundary element formulation for an anisotropic medium is briefly discussed. Good agreement has been found between numerical and analytical solutions for several example problems, demonstrating the accuracy of the present formulation. Numerical solutions are also presented for the problems of a deep circular tunnel and a basement excavated in a variety of jointed rock masses.