岩石试件端面摩擦效应数值模拟研究

试件端面摩擦效应直接影响试件内的塑性等效应变、侧向位移的分布和单元应力应变曲线。本文运用ANSYS中的接触单元模拟了平面应变状态下端面摩擦效应对塑性等效应变、侧向位移和单元应力应变曲线的影响,得到了不同摩擦系数时塑性等效应变及侧向位移的渐进变化形式。当接触面摩擦较小时,塑性等效应变图案为上下两个X形网络,侧向位移上下分布均匀;当接触面摩擦增大时,塑性等效应变网络向中部靠拢并且明显增大,侧向位移上下分布不均匀,中部较上下端面位移大;当试件端面侧向位移被限制,即摩擦力很大时,塑性等效应变网络变为一个X形局部化带,侧向位移分布更加不均匀,中部明显隆起。     

轴对称单轴压缩岩石材料应变局部化剪切带数值模拟研究

利用FLAC3D数值软件中的应变硬化—软化本构模型,模拟轴对称单轴压缩条件下岩石试件应变局部化剪切带的形成过程,同时对不同端面效应下长径比效应对应变局部化剪切带的影响进行分析。结果表明,岩石试件的应变局部化剪切带的启动、发展到最终形成是一个渐进的演化过程。端面效应和长径比效应均影响应变局部化剪切带的形状和岩石强度。光滑边界应变局部化剪切带形成于试件上下端面,摩擦边界形成于试件中部区域,且摩擦边界时试件的峰值强度高于光滑边界,峰值强度后的软化阶段陡峭于光滑边界。随着岩石试件长径比的增加,岩石的峰值强度在逐步降低,达到峰值强度所需的时步却呈增加的趋势。这一研究成果有助于加深理解岩石破坏过程。     

孔隙流体对岩体变形局部化的影响及数值模拟研究

本文首先叙述了孔隙流体对固体变形局部化影响的若干研究进展,并利用三维连续体快速拉格朗日分析,模拟了孔隙压力对岩体变形局部化应变场(剪切带)的影响。数值模拟结果表明,孔隙压力对变形局部化应变场有较大影响,随着孔隙压力的增加剪切带趋于不明显,岩体试件倾向于发生拉伸破坏。      

平面应变状态下岩石剪切带网络数值模拟研究

研究了围压、端面约束、试件高度及界面特性对剪切带网络形成的影响。首先对影响剪切带网络形成的因素（包括试件高度、围压及温度、试件端部轴向应变大小的计算方法的差异）进行了分析。之后采用FLAC 3D对平面剪切带网络进行了数值模拟,其中摩擦角及内聚力为应变软化。探讨了倾向于出现剪切带网络的若干条件。若端面缚束较强且存在一定的侧压力,试件中部出现多重剪切带。若端面缚束较弱,试件端部易于形成剪切带网络,且随试件高度的增加剪切带条数增加。界面法向和切向刚度越大,剪切带穿越断层并按其固有方向延伸能力越强,剪切带网络格局越明显。获得的数值结果可在实验中找到佐证,并且可以用来解释地震中的一些剪切应变局部化现象。     

基于无网格法的含缺陷岩石变形局部化数值模拟研究

利用自主开发的基于无网格法的岩石力学计算分析程序,采用Mohr-Coulomb破坏准则,考虑材料塑性屈服后的软化特性,研究单轴压缩条件下含一个缺陷的岩石试样的变形特性,探讨含缺陷材料的岩石变形局部化的演化规律。计算结果表明,岩石试件的破坏是一个渐进破坏的过程;加载过程中试件中的缺陷弱粒子首先出现粒子破坏并产生声发射,随加载进行试件产生微裂隙并扩张形成一个变形局部化带,局部化带扩张、发展,最终形成一个明显的剪切破坏带,缺陷粒子对试件的破坏起到了控制作用。岩石的声发射记录了岩石材料的塑性破坏过程,岩石试件的剪切破坏带与塑性区的发展具有非常类似的规律,但塑性区面积要大于剪切破坏带,从所处位置来看剪切破坏带位于塑性区之内,通过研究声发射特征信息发现变形局部化发生在应力-加载步（应变）达到峰值强度前。     

饱和砂土局部变形带模拟的有限元数值实现

基于有限变形理论,推导了Newton-Raphson 迭代算法在k+1步增量表达的矩阵形式,实现了饱和砂土变形局部化的有限元数值计算,得到了饱和砂土发生局部化变形的准则。基于Galerkin 方法,得到了位移场和应力场的空间离散化矩阵方程;由土体局部变形带的连续性条件,引入第1切线算子,推导出了砂土等颗粒状媒介发生局部化变形的必要条件。基于此核心算法,编制了有限元计算程序,模拟了饱和砂土在不排水条件下平面压缩过程中剪切带的形成与发展;通过比较分析,研究了有限元网格粗细对于土体局部变形带的影响,结果表明,网格粗细的病态依赖只是微小的,它只与变形条带的宽度有关,对于土体所表现出来的其他力学特性没有影响。     

岩体假三轴压缩及变形局部化剪切带数值模拟

岩体试件应力－应变全过程曲线，剪切带图案及其演化规律不仅在理论上，而且在工程实践中有着重要的意义，利用快速拉格朗日分析法模拟了尺寸效应，加载速率和围压等因素对假三轴试件抗压应力－应变全过程曲线，变形局部化启动，稳定，剪切带图案及演化规律的影响。     

上海黏性土剪切带形成的局部化变形试验研究

上海第五层灰色原状黏性土是上海市轨道交通车站换乘枢纽以及地下立交工程和超深基坑开挖施工所在土层,对这层土的局部化变形和剪切带形成规律的试验研究有助于深入了解它的基本力学性状,并为工程实践提供有意义的指导性结论。利用同济大学自制的平面应变仪,对该层黏土进行大量固结不排水平面应变试验,通过对试验结果的分析研究,定量得到上海黏性土变形局部化和剪切带形成的一些规律。     

不同弹模时圆形巷道围岩的应变局部化数值模拟

利用基于运动方程求解的FLAC程序模拟了不同弹性模量时圆形巷道的应变局部化过程。在计算中, 采用了"先加载, 后挖洞"方式。模拟结果表明, 在弹性模量中等或较小的条件下, 由于开挖卸荷而产生的围岩向空腔内部的"涌入"(运动)现象比较明显, 进而围岩中诱发了较多的破坏、V形坑或短剪切带式破坏。在上述条件下, 基于静力平衡的解析及数值方法不再适用, 而且是偏于不安全的。弹性模量越小, 开挖之后, 围岩维持其均匀周向(或环向)变形的能力越差, 这种轴对称变形越容易被打破, 也被打破得越早。尽管随着弹性模量的增加, 破坏单元数的变化并不是单调的, 但是在总体上, 随着弹性模量的增加, 破坏单元数目降低, 最大剪切应变增量急剧单调下降。      

不同围压下岩石抗压强度与变形参数尺寸效应的数值模拟

岩石抗压强度和变形参数是岩石工程设计的重要指标。由于岩石是典型的非均质材料,其强度和变形特性随样品尺 寸的变化而不同。本文采用PFC2D程序模拟了不同围压下不同尺寸岩样的压缩试验。结果表明（1） 岩样具有明显的尺寸效 应。同一围压下,尺寸越大,岩石强度、峰值应变和压缩模量越小,尺寸的变化对岩样的破坏模式影响较小;（2） 岩样具 有明显的围压效应。同一尺寸的岩样,随着围压的增大,岩石强度、峰值应变和压缩模量均增加,其中强度和峰值应变随 围压的增加呈线性增加。同时,随着围压的增大,岩石破裂模式由轴向劈裂破坏向剪切破坏变化;（3） 围压的存在会影响 岩样的尺寸效应。不同尺寸岩样的强度和峰值应变在相同围压区间内的增量基本相同,同时随着围压的增大,其强度和峰 值应变增加,进而使岩石强度和峰值应变的尺寸效应弱化;而不同尺寸岩样的压缩模量在相同围压区间内的增长率大致相 同,因而造成围压对压缩模量尺寸效应的影响较小     

工程岩体地质力学作用过程的数值模拟

工程岩体地质力学作用过程的数值模拟, 是研究岩石工程设计、施工问题的有效途径。本文提出了工程岩体计算数据智能化准备与确定子系统、地质力学作用过程的仿真数值模拟及其快速图形处理子系统, 并将其集成为一体。系统采用人机交互界面及多系统的自动转换, 使整个分析系统能充分发挥其模拟分析及图形处理功能, 方便实用。     

巷道围岩松动圈影响因素的敏感性数值计算

数值计算是预测围岩松动圈的一种有效方法,而各种影响因素的参数输入又是数值计算方法可靠性的关键所在。本文结合任楼煤矿中三轨道大巷已有实测位移和松动圈厚度,分析了各影响因素对围岩的影响性,研究了各影响因素对松动圈厚度的敏感性。计算结果表明,围压正向变化到一定程度后,增加的围压有一部分消耗在了已破碎围岩的二次破碎; 在实际岩层中开挖,只要巷道开挖就必然产生松动圈; 无构造应力情况下,围岩松动圈厚度主要取决于围岩的岩性和结构,各影响因素对松动圈厚度的敏感程度从高到低依次为围岩强度、巷道跨度、巷道断面尺寸、巷道围压和支护条件。     

采动诱发岩体移动破坏过程数值模拟研究

本文基于岩石介质的宏观非线性主要是由非均质性和各向异性造成的, 应用新的数值计算软件RFPA(2D), 对采动引起岩体失稳破坏的全过程进了数值模拟研究。     

裂隙岩体水力特征数值模拟试验的初步结果

本文利用文献中的数据进行数值模拟试验,结果表明:流量等效和渗透作用力变形效果等效互不包容;不同级别结构面体系难以用连续介质方法处理,同级别结构面体系在工程作用力变化范围内,可由其静态水力特征推演动态特征;工程岩体结构面间充分连通的假设过于粗糙,常规试验手段获取的量值倾向于局部特征的反映。     

泊松比对岩样破坏模式及全部变形特征的影响

利用编写的计算岩样全部变形特征的FISH函数, 采用FLAC模拟了泊松比不同时单缺陷岩石试样的破坏及全部变形特征。在峰前及峰后, 本构模型分别取为线弹性模型及莫尔库仑剪破坏与拉破坏复合的应变软化模型。高泊松比使岩样发生由单一剪切破坏向复杂破坏转变、破坏区域的面积增加、剪切带倾角降低, Coulomb、Roscoe及Arthur理论对此无法解释。不同泊松比时计算得到的峰前应力-轴向应变曲线、应力-侧向应变曲线、侧向应变-轴向应变曲线、体积应变-轴向应变曲线的线性阶段与平面应变压缩条件下的线弹性解吻合。若泊松比超过1/3, 通过计算得到的平面应变压缩泊松比可大于0.5, 这被数值模拟确认。泊松比的增加使峰后的侧向应变-轴向应变曲线、体积应变-轴向应变曲线、计算得到的泊松比-轴向应变曲线变得不陡峭, 使峰后的应力-侧向应变曲线变得陡峭, 使破坏的前兆变得不明显。      

岩体裂隙扩展过程的数值模拟

本文根据断裂力学理论及节理岩体的等效连续模型,探讨了岩体裂隙扩展过程的数值模拟方法,文末算例说明了方法的有效性。     

颗粒配比对岩石力学特征影响的数值模拟研究

岩石是矿物颗粒的集合体同时也是一种重要的非均质材料,了解它的力学特征对岩土工程及矿产开采都具有重要的指导作用。作为典型的颗粒材料,颗粒单元体的粒径分布配比必然影响着岩石的宏观力学表现。通过设置不同体积配比下的颗粒材料单元体,利用PCF2D软件模拟了相同颗粒材料单元体不同配比下岩石模型的力学特征。模拟结果表明颗粒单元体配比对岩石的力学特征有明显的影响。在模拟过程中大颗粒的配比显著影响着岩石的抗压强度,大颗粒含量相对越高,抗压强度越大。而细颗粒的配比影响着岩石的抗拉强度,细颗粒含量相对越高,抗拉强度越大,但是过多的细颗粒会降低岩石的抗拉强度。考虑岩石压缩过程中裂缝形态的影响。结果表明均匀分布、5:2:3、7:2:1的颗粒配比形成了贯穿裂缝,而1:2:7和3:2:5的颗粒配比未能形成贯穿裂缝,且细颗粒配比越高,裂缝数目出现高值的概率也越大。