考虑剪应力作用时横观各向同性岩体中圆形巷道位移反分析的惟一性

推导了考虑剪应力作用时横观各向同性岩体中圆形巷道的位移解析解,利用参数可辨识条件对考虑剪应力作用时横观各向同性岩体中圆形巷道位移反分析的惟一性进行了探讨。结果表明,无论布置多少个测点也不能惟一地反分析出所有7个参数,必须至少已知3个参数,才有可能惟一地反分析出其他参数;水平与竖直地应力分量是否为 0 对反分析惟一性有一定影响;3个地应力分量的可辨识性最好,各向同性面上的弹性模量和泊松比次之,垂直于各向同性面方向的弹性模量和泊松比最差。与不考虑剪应力时位移反分析惟一性结果相比,考虑剪应力时参数的可辨识性顺序不变,但是条件惟一的比例降低,不惟一和绝对惟一的比例增加,说明考虑剪应力能够提高位移反分析的惟一性。工程实例计算表明,推导的位移解析解和反分析惟一性探讨结果能够很好地指导岩体参数反演计算。     

横观各向同性岩体中马蹄形隧洞的位移反分析方法

将各向异性岩体中任意形状隧洞的位移解析解与智能算法--差分进化法（DE）结合,研究横观各向同性岩体中在各向异性面上开挖马蹄形隧洞的位移反分析问题,获得关于参数灵敏度、参数可惟一辨识性和测线布置方案的有关结论。计算发现,通过DE算法可获得很高精度的参数辨识结果;参数的可惟一辨识性不仅与灵敏度有关,还与具体的参数组合形式有关,在给定竖向地应力和各向同性面角度的情况下最多可以高精度惟一辨识5个参数;当待辨识参数的个数大于等于4时会产生不可完全高精度惟一辨识的参数组合,当待辨识参数的个数小于等于3时,任意的参数组合形式均可以被高精度惟一辨识。     

一种横观各向同性岩体蠕变模型

层状岩体由于层理面或者结构面的存在,在力学上具有横观各向同性的特点,现有的各向同性蠕变模型难以全面反映横观各向同性岩体的蠕变力学特性。为构建能够反映横观各向同性岩体的三维蠕变模型,以能反映瞬时应变、减速蠕变和稳态蠕变特征的黏弹性Burgers模型为基础,基于常泊松比假定,在三维各向同性蠕变本构方程的基础上,按照算子替换的方法,将横观各向同性柔度矩阵代替各向同性柔度矩阵,并考虑了平行和垂直层理方向岩体蠕变力学行为的差异性,推导了横观各向同性岩体的三维蠕变本构方程。根据本构方程的特点,提出了根据平行和垂直方向岩体蠕变试验结果进行三维蠕变本构模型中蠕变参数的辨识方法。将提出的模型应用于三轴蠕变试验参数辨识,从而获得了一套完整的三维蠕变参数,并与试验结果进行对比分析,从而验证了所提模型与方法的合理性和有效性。进一步,指出了传统蠕变试验设计方案的局限性,给出了横观各向同性材料蠕变试验设计建议。研究成果为研究岩体三维蠕变机制提供了新思路,可对岩体蠕变试验设计提供相应的科研支撑。     

层理性页岩力学参数对水力压裂裂缝形态的影响分析

岩石力学性质及地应力等地质力学参数是进行储层压裂改造的基础数据,影响着整个压裂设计的准确性。由于页岩层理发育及矿物组成的复杂性、页岩力学性质的横观各向同性,所以基于传统力学模型获取的压裂参数设计将丧失意义。以横向各向同性模型体现层理发育页岩储层的力学性质,通过简化的横纵弹性模量和泊松比可以较好地表达出页岩的横观各向同性,并总结了横观各向同性地应力解释方法。基于上述理论,应用商业化水力压裂数模软件对四川盆地某层理页岩储层进行了压裂设计,比较了基于横观各向同性模型与传统模型解释力学参数差异,非常规储层动静态力学参数转换,同时通过改变射孔位置研究了起裂点对裂缝形态的影响。研究发现,横观各向同性解释的杨氏模量、地应力要略高于传统方法,同时横观各向同性解释的地层地应力差较小,所以受不同射孔位置的变化对裂缝几何的形态影响要远远高于各向同性模型模拟的结果,特别是缝高变化更为明显。同时数据表明,在高应力区射孔将形成更为狭窄缝宽的水力裂缝,这会严重影响后续支撑剂的进入,以上结论将为提高页岩储层的压裂设计精确性给予重要帮助。     

数值计算图谱及其在位移反分析中的应用

本文采用了一种简便经济而又具有较强的实用价值的方法——数值计算图谱方法,并将其应用于军都山双线铁路隧道工程的岩体力学特性参数(弹性模量E和水平向地应力分量P)的位移反分析,获得了可靠的结果。     

风化岩体中压力型锚索锚固段的剪应力分析

岩体在风化过程中,其物理力学性质会随着风化程度的不同有所差异。结合均质岩体中压力型锚索受力特点,以及内部对称荷载作用下横观各向同性弹性体内的剪应力分布特点,推导出了压力型锚索在横观各向同性弹性岩体内的剪应力分布方程。假定风化岩体为横观各向同性弹性体,其弹性模量及剪切模量均按指数函数变化,分析结果表明：①由于风化岩体的弹性模量随着深度的变化,锚索锚固段的剪应力也随着岩体的弹性模量的增加（减小）而增加（减小）;②在承压板附近,锚索锚固段的剪应力随锚索孔半径的增加而减小,由于风化岩体剪切模量的变化,因而导致剪应力的变化量很大;③在承压板附近,锚索锚固段接触面上的剪应力随着垂直于岩体横向等向面的剪切模量的改变而相应的改变。研究成果可为实际工程中压力型锚索设计提供参考。     

水平条形均布荷载下横观各向同性弹性地基的位移分析

分析了水平条形均布荷载下横观各向同性地基位移规律,归纳了位移随弹性模量比n之不同而变化的规律,提出了水平条形均布荷载下横观各向同性地基位移的近似简化解,以供岩土工程设计人员参考。     

随机分布贯穿裂隙岩体变形特性研究

洞室开挖后,其周边通常会产生许多随机分布的贯穿裂隙,直接影响洞室围岩稳定,研究随机分布贯穿裂隙岩体的变形及变形特性具有重要意义。基于线弹性理论和线性刚度理论计算岩石和裂隙的位移,用概率的方法建立了随机分布贯穿裂隙岩体变形的计算模型,给出了随机分布贯穿裂隙岩体的等效弹性模量和等效泊松比,研究了岩石和裂隙的材料参数和几何参数对岩体等效弹性模量和等效泊松比的影响。可得如下结论：等效弹性模量和等效泊松比随着岩石弹性模量的增大而增大;等效泊松比随着岩石泊松比的增大而增大;等效弹性模量和等效泊松比随着裂隙法向刚度的增大而增大;随着剪切刚度的增大,等效弹性模量逐渐增大,而等效泊松比则逐渐减少;随着裂隙平均间距的增大,等效弹性模量逐渐减小,等效泊松比在平均倾斜角较小时逐渐增大,在平均倾斜角较大时逐渐减小;随着裂隙平均倾斜角的增大,等效弹性模量先减小后增大,而等效泊松比先增大后减小。模型能较全面地考虑构成岩体的岩石和裂隙的材料参数与几何参数对岩体变形的影响,其结果对研究洞室围岩的变形和工程设计有一定的参考价值。     

电爆冲击波对储层致裂效果影响因素的数值模拟研究

电爆冲击波在油气储层增渗解堵中具有巨大应用潜力,为了解各因素对电爆冲击波致裂效果的影响程度,本文采用ANSYS/LS[CD*2]DYNA软件,通过改变地应力、弹性模量、抗拉强度参数进行数值模拟。结果显示,地应力对致裂效果具有强烈的抑制作用,在双向等压地应力作用下,裂纹区形状为圆形,并且随着地应力的增大,岩体破裂度、破裂半径均逐渐减小。在双向不等压地应力作用下,随水平主应力差值增大,裂纹扩展方向由圆形逐渐转变为椭圆形,具有沿最大水平主应力方向扩展的特点;弹性模量对致裂效果具有重要影响,随着弹性模量的增大,岩体破裂度、破裂半径均呈增大趋势;随抗拉强度增大,岩体破裂度、破裂半径均减小。     

宜兴抽水蓄能电站试验洞的反分析研究

结合江苏宜兴抽水蓄能电站地下厂房试验洞量测资料的分析,研究了成层各向异性岩体中洞室工程问题的反分析方法。采用分层均质模型进行反分析计算。计算区域的各层岩层均被分别视为均质各向同性体,初始地应力假设呈线性分布。有限元分析中,对倾斜成层的岩层、软弱夹层及断层节理等分别划分了单元,据以模拟成层岩体的各向异性特征。为使计算过程简化,对离洞室较远的岩层拟按地质资料对材料特性参数赋值,反分析计算的目标未知数为初始地应力和试验洞相邻围岩的弹性模量。结果表明,与工程实践符合较好。     

锚杆支护圆形隧洞的等效强度参数及可靠性分析

通过采用均匀化方法,研究了圆形隧洞的锚杆支护特性,将高密度支护模式下的岩石和锚杆复合体考虑成均匀、连续、强度参数增强的等效材料,简化了岩石和锚杆间复杂的力学耦合问题。通过定义锚杆密度参数来反映不同支护模式的特性,建立锚杆密度参数与Mohr-Coulomb屈服准则中主要参数之间的关系,推导出等效弹性模量、等效黏聚力和等效内摩擦角的表达式,并分析比较了隧洞在支护前后的位移情况。结合可靠性理论,采用容许极限位移量作为失稳判据,分析了隧洞在支护前后的可靠性指标与破坏概率,结果表明,文中提出的方法简单可靠,锚杆支护对隧洞的位移限制效果明显,可显著提高隧洞的可靠性。     

层状岩体地下洞室破坏模式数值模型及验证

基于层面法向局部坐标系中的横观各向同性弹性本构模型,根据层面产状方位角,推导出整体坐标系中弹性应力的转换方法。针对层状岩体的各向异性特征,将岩体的破坏模式细分为5种,分别对每一种破坏模式建立各自的判别函数和屈服准则,根据塑性流动正交法则推导出相应的应力修正迭代计算方法。通过单轴、三轴压缩数值试验,对层状岩体的强度和变形性质进行研究,分析层面倾角和围压对层状岩体力学行为的影响。将所建模型与FLAC3D自带的横观各向同性弹性模型以及遍布节理模型进行对比,验证模型的正确性及精度。同时,将模型应用于三向等压理想圆形隧洞的开挖计算中,分析层状岩体破坏方式受层面方位角影响的变化规律,分析结果表明层面对围岩破坏模式和塑性区扩展方向起到了控制作用。     

Response of a tunnel deeply embedded in a viscoelastic medium

This paper presents a three‐dimensional energy‐based solution for the time‐dependent response of a deeply embedded and unsupported semi‐infinite tunnel of circular cross‐section. The tunnel is taken to be excavated quasi‐instantaneously from an infinite rock body that initially exhibits an isotropic stress state and that is made up of a homogeneous, isotropic and viscoelastic material. The viscoelastic behaviour is modelled by means of Burger's model, and the rock is taken to behave volumetrically linear elastic and to exhibit exclusively deviatoric creep. This viscoelastic problem is transformed into the Laplace domain, where it represents a quasi‐elastic problem. The displacement fields in the new solution are taken to be the products of independent functions that vary in the radial and longitudinal directions. The differential equations governing the displacements of the system and appropriate boundary conditions are obtained using the principle of minimum potential energy. The solutions for these governing equations in the Laplace domain are then obtained analytically and numerically using a one‐dimensional finite difference technique. The results are then transformed back into the time domain using an efficient numerical scheme. The accuracy of the new solution is comparable with that of a finite element analysis but requires much less computation effort. Copyright © 2011 John Wiley & Sons, Ltd.     

Updating the models and uncertainty of mechanical parameters for rock tunnels using Bayesian inference

Rock mechanical parameters and their uncertainties are critical to rock stability analysis, engineering design, and safe construction in rock mechanics and engineering. The back analysis is widely adopted in rock engineering to determine the mechanical parameters of the surrounding rock mass, but this does not consider the uncertainty. This problem is addressed here by the proposed approach by developing a system of Bayesian inferences for updating mechanical parameters and their statistical properties using monitored field data, then integrating the monitored data, prior knowledge of geotechnical parameters,and a mechanical model of a rock tunnel using Markov chain Monte Carlo(MCMC) simulation. The proposed approach is illustrated by a circular tunnel with an analytical solution, which was then applied to an experimental tunnel in Goupitan Hydropower Station, China. The mechanical properties and strength parameters of the surrounding rock mass were modeled as random variables. The displacement was predicted with the aid of the parameters updated by Bayesian inferences and agreed closely with monitored displacements. It indicates that Bayesian inferences combined the monitored data into the tunnel model to update its parameters dynamically. Further study indicated that the performance of Bayesian inferences is improved greatly by regularly supplementing field monitoring data. Bayesian inference is a significant and new approach for determining the mechanical parameters of the surrounding rock mass in a tunnel model and contributes to safe construction in rock engineering.     

考虑轴向力和渗透力时圆形隧道广义Hoek-Brown解

为了研究渗透水压力和轴向应力共同作用时隧道围岩的应力和位移变化趋势,将圆形隧道简化为轴对称模型,假定渗透场以渗透体积力作用在原应力场,以围岩开挖断面为假定平面,引入垂直于该平面的轴向应力。基于广义Hoek-Brown强度准则和非关联流动法则,推导出考虑轴向应力和渗透场共同作用时弹-脆-塑性围岩的应力和位移非线性解,采用数值算例分析了轴向应力和渗透力共同作用时隧道围岩塑性区应力场和位移场的分布规律。计算结果表明：与无渗透水压力作用下的模型相比,渗透力作用使得围岩塑性区各点位移增大,并且内外水头差越大,位移增大越明显。当轴向应力为中主应力时,围岩塑性区半径和塑性区各点应力增大,轴向应力为大主应力和小主应力时,围岩塑性区半径和塑性区应力变化较小。因此,渗透力的存在不利于隧道围岩的稳定性,并且轴向应力的大小对于富水地区隧道围岩的应力和位移分布具有较大影响。在施工设计时考虑渗透力以及轴向应力的共同影响对于保证隧道围岩稳定性具有重要意义。     

Elastic solutions for a transversely isotropic half-space subjected to buried asymmetric-loads

Elastic closed-form solutions for the displacements and stresses in a transversely isotropic half-space subjected to various buried loading types are presented. The loading types include finite line loads and asymmetric loads (such as uniform and linearly varying rectangular loads, or trapezoidal loads). The planes of transverse isotropy are assumed to be parallel to its horizontal surface. These solutions are directly obtained from integrating the point load solutions in a transversely isotropic half-space, which were derived using the principle of superposition, Fourier and Hankel transformation techniques. The solutions for the displacements and stresses in transversely isotropic half-spaces subjected to linearly variable loads on a rectangular region are never mentioned in literature. These exact solutions indicate that the displacements and stresses are influenced by several factors, such as the buried depth, the loading types, and the degree and type of rock anisotropy. Two illustrative examples, a vertical uniform and a vertical linearly varying rectangular load acting on the surface of transversely isotropic rock masses, are presented to show the effect of various parameters on the vertical surface displacement and vertical stress. The results indicate that the displacement and stress distributions accounted for rock anisotropy are quite different for those calculated from isotropic solutions. Copyright © 1999 John Wiley & Sons, Ltd.     

Generalized Hoek-Brown solution of circular tunnel considering effects of axial stress and seepage forceEI北大核心CSCD

This paper studies the joint effect of seepage force and axial stress on the stress and displacement of circular tunnel. The circular tunnel is simplified as an axisymmetric model and the seepage field is simplified as volumetric force in the stress field. The excavation cross-section of surrounding rock is assumed as a plane as well, and an axial stress perpendicular to the plane is further introduced. Nonlinear solutions for the stress and displacement of circular tunnel are deduced considering the joint effect of axial stress and seepage force, based on the generalized Hoek-Brown failure criterion and the non-associated flow rule in elastic-brittle-plastic rock mass. Numerical simulations are also employed to analyze the distribution of stress field and displacement field in plastic zone of a circular tunnel under the joint effect of axial stress and seepage force. The calculated results show that the displacement in plastic zone increases significantly with the gradient increment of the seepage pressure, compared with the situation without seepage force. The radius and stress of surrounding rock in plastic zone increase when axial stress is the intermediate principal stress, while the radius and stress have less change when axial stress is the major or minor principal stress. It can be concluded that the seepage force has negative effects on the stability of circular tunnel, and the axial stress significantly influences the stress and displacement of the circular tunnel, especially in water-rich areas. Therefore, it is necessary to consider the joint effects of axial stress and seepage force to ensure the stability of circular tunnel in water-rich area.