不同围压条件下的圆形巷道岩爆过程模拟

利用FLAC模拟了不同围压条件下圆形巷道的岩爆过程。为了模拟巷道开挖,利用编写的FISH函数删除巷道内部的单元。岩石服从摩尔库仑剪破坏与拉破坏复合的破坏准则,破坏之后呈现应变软化-理想塑性行为。模拟结果表明:当围压较低时,剪切应变集中区域呈圆环状,围岩能保持稳定,不出现剪切带;当围压增加到一定程度时,围岩中出现了"狗耳"形的V形坑,发生岩爆,但围岩也还能保持稳定;当围压进一步增加时,围岩中出现了多条狭长的剪切带,巷道的整个断面均遭到了破坏,发生强烈的岩爆。随着围压的增加,V形岩爆坑变大、变深,剪切带花样的对称性变差;在高围压时,剪切带花样与塑性力学中的滑移线网有类似之处。     

三维盾构隧道开挖面极限支护压力的上限解

为了研究掘进过程中盾构隧道开挖面附近土体的稳定性,假设土体为纯粘性土,且为理想弹塑性材料,服从摩尔-库仑屈服准则,参照矩形基础承载力问题的三维Hill机构,建立了盾构隧道开挖面土体处于主动、被动极限平衡状态下的破坏模式。从塑性极限分析上限法的基本原理出发,推导出三维盾构隧道开挖面极限支护压力上限解的计算公式。结合某算例,讨论了三维盾构隧道开挖面主动、被动极限支护压力上限解与平面应变上限解的关系,以及开挖面极限支护力与土体粘聚力、h/D的相互关系。研究结果表明:三维盾构隧道开挖面主动极限支护力要比平面应变上限解小,而被动极限支护力要比平面应变上限解大;随着粘聚力的增加,三维盾构隧道主动、被动极限支护力与平面应变上限解的差别也逐渐增大;盾构隧道开挖面主动、被动极限支护力的上限解均随h/D的增大而增大,且h/D越大,盾构隧道开挖面主动、被动极限支护力随土体粘聚力的变化速率也越快。     

弹性模量依赖围压的圆形洞室围岩破坏行为数值模拟

在考虑拉伸截断及非线性屈服准则的自定义本构模型中,又考虑了弹性模量(弹模)受围压的影响,模拟了两种网格、不同本构模型及加载方式条件下的洞室围岩的破坏行为。在第1种网格中,洞室围岩被剖分为辐射网格,采用"先挖洞,后加载"的方式及理想弹-塑性本构模型进行计算;在第2种网格中,洞室围岩被剖分为矩形网格,采用"先加载,后挖洞"的方式及弹-脆-塑性本构模型进行计算。计算结果表明,对于第1种网格,随着弹模的围压影响系数的增加,洞室围岩的应变、位移及塑性区半径均有所减小,被监测单元的环向应力的峰值有所降低,不同本构模型的计算结果的差别逐渐显现。对于第2种网格,当弹模的围压影响系数较低时,围岩中形成了4个对称的V形坑,随着弹模的围压影响系数的增加,V形坑的个数及尺寸明显减少,这可由洞室围岩的应力、位移及剪切应变增量均有所减小得到解释。     

不同轴压/侧压下逐步开挖马蹄形巷道围岩分区破裂化数值模拟

基于空间局部化视角研究分区破裂化的发生机制。采用三维非均质应变软化模型,在逐步开挖和不同轴压与侧压之比条件下研究了马蹄形（U形）巷道围岩不同剖面和测线上剪切应变增量的分布特征。结果表明：U形巷道围岩的分区破裂化呈圆环形;分区破裂化可出在现掌子面的前、后方;在掌子面前方,随着远离掌子面,分区破裂化有所减弱,直至消失;轴压与侧压之比越大,分区破裂化越明显;在垂直于巷道轴线的不同平面上,分区破裂化有所不同,这与不同剪切带的特征不同及发育不同有关;共轭剪切带的出现,会使分区破裂化变得复杂。启动于巷道表面的空间圆锥剪切面的充分发展产生了分区破裂化。     

黏土半巷道围岩模型过巷道轴线平面破坏过程实验观测∗

为了观测过巷道轴线平面上应变局部化带及裂纹的发展演化规律,开展了黏土半巷道围岩模型三轴压缩实验研究。半巷道围岩模型的观测表面处于平面应变状态,且分步开挖巷道。在每次巷道开挖之前,对半巷道围岩模型卸荷。每次巷道开挖之后,再对半巷道围岩模型加荷。采用数字图像相关方法观测了观测表面的最大剪切应变场,并统计了裂纹面积百分比随纵向应变的演化规律。得到以下结论:当纵向应变较高时,在半巷道围岩模型的观测表面上,出现了由巷道两帮或两帮与掌子面交点处发展出的应变局部化带;另外,还出现了距离巷道表面较远的应变局部化带。由此可以推出,垂直于巷道轴线平面上将会出现分区破裂化。裂纹面积百分比随纵向应变的增加而增加;随着阈值(可识别的裂纹面积的最小值)的增加,裂纹面积百分比降低,但降低的速度越来越慢。剪裂纹主要位于巷道的两帮;拉裂纹遍布各处,呈现一定的等间距性。     

发电机定子冷却水管路堵塞的原因分析及处理

利用FLAC模拟了不同岩石峰后脆性条件下圆形巷道围岩的破坏区分布及能量释放规律。在计算中,采用了应变软化本构关系及"先加载,后挖洞"方式。模拟结果表明,随着峰后脆性的增加,围岩中发生破坏的单元变多,释放的弹性应变能总量增加,巷道围岩越来越难于再次达到平衡状态;尽管在一些算例中出现了沿环向发展的破坏区,但是不认为这与分区破裂化现象有关;难于采用二维的滑移线或破坏特征线解释分区破裂化现象。围岩分区破裂化现象的研究可从三维连续介质模型的空间局部化视角出发,亦可从二维颗粒体模型的颗粒相互作用机制视角出发,这两种思路都引入了二维连续介质模型中不具备的因素。     

卸荷时间对圆形巷道围岩开裂及径向应力波传播影响的数值模拟

采用提出的连续-非连续方法,研究了卸荷时间(卸荷快慢)对圆形巷道围岩开裂、径向应力波传播及围岩径向应力随时间演变规律的影响。该方法借鉴了拉格朗日元法、变形体离散元法及虚拟裂纹模型的思想。首先,研究了卸荷时间对围岩开裂的影响;然后,研究了卸荷时间对径向应力波传播的影响;最后,研究了卸荷时间对围岩径向应力随时间演变规律的影响。结果表明:卸荷越快,围岩中产生的裂纹越多,最大不平衡力越大,波动越剧烈。增加卸荷时间,开挖边界附近的开裂位置由连续分布向间隔分布转变:当卸荷时间较短时,开挖边界附近围岩呈现连续分布的开裂形态,与这些位置的最大主应力均超过抗拉强度引起拉裂有关;当卸荷时间适中时,围岩呈现出与V形坑类似的开裂形态,可能与这些位置存在剪应力有关。卸荷越快,处于扰动区(径向压应力下降区)内的应力环越多,越清晰,围岩径向应力波动越剧烈,所达到的最大径向拉应力越大或最小径向压应力越小。     

先动后静荷载联合作用下结构性软粘土力学响应

为研究结构性软粘土在先动后静荷载作用下的力学响应特性,利用动静扭剪仪对天津滨海结构性软粘土进行了不同动应力幅值和波形的循环荷载作用,当动力累积塑性应变达到不同控制值后,再对土样进行静力三轴剪切试验。试验结果表明,在循环荷载作用下,累积塑性变形不断增大,土体结构破损加剧,随着先期动应力幅值和累积塑性应变的增大;之后,静力三轴应力应变关系曲线由应变软化型逐渐变为硬化型,其区分界限(1.25%)即为结构屈服应变;以结构屈服应变为分界点,相同围压下静力破坏应力随累积塑性应变的增大呈先快后缓的下降趋势。     

地质体材料剪切带内部的常剪切应变点及速度分布分析

对于峰后线性应变软化的地质体材料,剪切带内部的总剪切应变等于弹性剪切应变(由经典弹性理论描述)及由微结构效应而引起的局部塑性剪切应变(由非局部理论或梯度塑性理论描述)之和。若剪切应力—塑性剪切应变曲线的斜率的绝对值(称之为软化模量)小于剪切弹性模量的两倍,则在剪切带的任一剖面内存在两个总剪切应变不依赖于剪切应力的点,称之为常剪切应变点。在这两个点上,弹性剪切应变的降低和局部塑性剪切应变的增加处于平衡状态,总剪切速度达到它的最大值或最小值。在两个常剪切应变点之间,局部总速度随剪切应力率的降低而增加。剪切带内部的局部总速度分布是非线性的,这与通常采用的剪切带内部速度的线性分布假定(忽略微结构效应)不同。     

基于非线性屈服准则及主应力判据的圆形巷道围岩岩爆过程的数值模拟

鉴于考虑拉伸截断的线性莫尔-库仑屈服准则不能考虑岩石在高应力条件下的非线性屈服特征,根据虎克-布朗本构模型,将上述屈服准则在受压区的线性屈服函数修改为非线性形式。采用3种岩爆的主应力判据,判别圆形巷道开挖之后围岩中各种级别岩爆的分布及演变规律。计算采用"先加载,后挖洞"的方式,岩石服从弹—脆—塑性本构模型。研究发现,发生高级别岩爆的单元数少于发生低级别岩爆的单元数;轻微岩爆区的形态更接近于塑性区及剪切应变增量的高值区。根据巴顿判据,发生重岩爆的单元数较多,而且主要发生的是拉伸岩爆;根据陶振宇判据,发生轻微及中等岩爆的单元数均多于根据谷明成-陶振宇判据判别的结果。上述两种判据判别的结果均表明,在不考虑拉伸岩爆的条件下,只有位于围岩内部的单元才有可能发生高级别的岩爆,而位于巷道表面的单元一般仅发生轻微岩爆。     

孔隙压力对含随机缺陷岩石破坏过程及全部变形特征的影响

对于平面应变压缩条件下含有随机缺陷的岩样,利用FLAC研究了孔隙压力对岩样破坏过程、全部变形特征及前兆的影响。以前编写的若干FISH函数,被用于生成缺陷和计算轴向、侧向、体积应变及侧向应变与轴向应变比值的负值(计算得到的泊松比)。在峰值应力之后,密实岩石单元服从线性应变软化行为及随后的理想塑性行为,而材料缺陷呈现理想塑性行为。当孔隙压力较高时,应力—侧向应变曲线具有一平台;破坏的前兆更明显;变形后岩样的体积总是大于原始体积;在初始加载阶段、均匀变形阶段及峰后变形阶段,由于明显的侧向膨胀,计算得到的泊松比远大于0.5。当孔隙压力较低时,在峰值应力之前,变形后岩样的体积小于原始体积,体积扩容出现于峰值应力之后,引起了负的体积应变。利用广义虎克定律,解释了平面应变弹性状态下数值结果的合理性。对岩样进行带状区域扫描后,确认随机缺陷的初始分布与岩样的最终破坏形态紧密相关。     

地震块体模型的共轭剪切破裂带数值模拟

对地震块体模型的应变率场的数值模拟结果表明 ,地震块体模型 ,尤其是当界面法向和切向刚度较大时 ,能产生明显的共轭剪切破裂带。地震块体模型存在老断层 ,呈现各向异性 ,因而 ,促进了某一方向上的主条带充分发展 ,与主条带共轭的次条带被抑制。存在复杂地质构造可能是地震带的间隔距离并不严格相等的原因。地震发生在剪切破裂带交叉部位的原因是这一位置剪切应变率较高。考察应变率场能揭示出地震块体模型的应变局部化特征。在构造运动作用下 ,新生构造与老构造并不重叠是可能的 ,因而 ,新生构造的交汇点就不一定出现在老断层上。地震块体模型的中心块体的尺寸越大 ,中心块体越稳定     

基于滑移线网络法的黄土隧道坍塌拱分析及其承载评价∗

针对隧道稳定分析计算方法复杂、计算量大等问题,基于滑移线场提出一种简洁方便的图解法具有重要的理论意义和工程意义。首先发展完善了无支护黄土隧道全断面围岩滑移线网络图分析方法,然后,通过MATLAB实现了滑移线网络法的可视化分析程序,结合试验性黄土隧道塌落监测成果验证了无支护黄土隧道全断面围岩滑移线网络图分析方法的分析程序,最后,评价分析了黄土抗剪强度内摩擦角和黏聚力对黄土隧洞围岩坍塌范围及土拱极限承载力的影响。结果表明:黄土土拱极限承载力随黏聚力的增大而增大且敏感程度升高,黄土土拱极限承载力随摩擦角的增大而减小且敏感程度降低,坍塌范围随摩擦角的增大而减小。     

圆形巷道围岩的环向应力及发生岩爆单元分布的数值模拟

利用FLAC模拟了不同水平方向压力(小于竖直方向压力)及岩石峰后不同脆性条件下的圆形巷道破坏过程.岩石服从莫尔-库仑剪破坏与拉破坏复合的破坏准则,破坏之后呈现应变软化-理想塑性行为.监测了模型中第1象限对角线上的单元环向应力分布及演化规律.根据徐林生和王兰生提出的环向应力岩爆判据,判断模型中各单元是否发生岩爆.模拟结果...