不同强度岩石中开挖圆形巷道的局部化过程模拟

利用FLAC模拟了不同粘聚力条件下圆形巷道的局部化过程。为了模拟巷道开挖,利用编写的F ISH函数删除巷道内部的单元。岩石服从莫尔库仑剪破坏与拉破坏复合的破坏准则,破坏之后呈现应变软化-理想塑性行为。文中模拟分为3步:首先,将静水压力施加在模型上,直到达到静力平衡状态;然后,利用编写的F ISH函数开挖巷道;最后,计算重新开始,直到达到静力平衡状态或者塑性流动状态。模拟结果表明,随着粘聚力的降低,巷道围岩的破坏模式首先由孔壁附近零星单元的破坏向4个对称的小V形坑式剪切破坏转变,然后由包含若干小V形坑的大V形坑式剪切破坏向巷道全断面的破坏转变。前三者破坏发生后,巷道围岩仍然能保持稳定。与最大塑性拉伸应变相比,最大剪切应变增量、最大塑性剪切应变要高得多;最大剪切应变增量、最大塑性剪切应变相差不大;随着粘聚力的增加,三者均越来越小。     

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

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

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

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

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

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

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

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

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

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

泥化夹层压缩试验的损伤特征及损伤方程研究

泥化夹层的损伤破坏是影响边坡稳定性的重要因素,为研究其损伤机理,进行了不同围压条件下的压缩试验,通过试验数据和试样破坏外观分析损伤情况;利用CT扫描仪在试验不同加载时刻对泥化夹层进行实时CT扫描,结合MATLAB图像二值化处理,得到不同围压下泥化夹层损伤识别图像,反映了内部损伤发展情况;最后根据推广的应变等价原理,通过变形模量确定损伤变量,获得泥化夹层损伤变量与轴向应变的关系。结果表明:①泥化夹层强度普遍表现出随围压增加而增加的规律,破坏应变取为10%,由于泥化夹层自身组成及结构存在差异性,导致个别试样的强度规律不明显;②泥化夹层内部存在薄弱层,加载过程中会出现阶段性的小幅度塌缩,强度曲线上反映为阶段性的应变突增现象;③无围压和一定围压情况下的损伤破坏方式不同,无围压时为脆性破坏,裂缝萌生、开展,一定围压时表现为塑性破坏,泥化夹层被压缩,孔洞裂隙闭合,破坏方式改变临界围压介于0～50kPa;④通过本文方法建立的损伤方程与不同围压下的室内试验结果拟合良好,该方法合理可行。     

深埋硬岩隧洞岩爆风险与结构面倾角关系探讨∗

针对深埋硬岩隧洞即时性应变-结构面滑移型岩爆预测预警问题。以埋深 2 375 m 的某深部地下实验室 7#、 8#实验室发生的一次极强岩爆作为现场测试对象,基于钻孔摄像技术实现隧洞围岩内结构面展布可视化,通过岩爆诱发洞段与安全通过洞段围岩内结构面展布差异的对比分析、围岩弹性应变能分布特征分析与岩爆统计规律分析,对岩爆孕育中结构面的控制作用进行了深入的探讨。研究结果表明,深埋硬岩隧洞边墙与围岩内发育的结构面平行或是小角度(小于 30°)相交时具有更高的即时性应变-结构面滑移型岩爆风险。7#、8#实验室钻孔摄像结果显示,岩爆诱发洞段围岩内结构面近乎平行于边墙发育,3 处安全通过洞段围岩内结构面与隧洞边墙均为大角度 (大于 30°)相交。围岩内不同的结构面展布形态为岩爆孕育提供了不同的能量累积环境,当围岩内结构面平行于边墙或是小角度相交于边墙时 ,结构面的存在将有利于岩爆孕育 。 所有发生于锦屏二级水电站排水洞 SK8+200~8+800 洞段的岩爆统计数据表明,隧洞边墙与围岩内结构面相交夹角较小(小于 30°)的区域具有更高的即时性应变-结构面滑移型岩爆风险,这与 7#、8#实验室现场测试及数值模拟结果相吻合。     

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

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

单轴压缩土样剪切带上主应变轴偏转角的分布及规律统计北大核心CSCD

根据土样纵向应变较高时最大剪切应变场中清晰剪切带所处位置沿剪切带方向布置平直测线,获得6条测线上从剪切带萌生至微裂纹刚出现过程中主应变轴偏转角的分布及演变规律,并对其进行了统计。研究发现：当微裂纹刚出现时,测线上主应变轴偏转角的高峰和低谷对应于2条共轭剪切带的交汇处附近或剪切带上最大剪切应变的局部高值区附近;在土样中部,剪切带上主应变轴偏转角主要分布在-5°-5°;随着纵向应变的增加,剪切带上主应变轴偏转角的标准差先降低,后趋于稳定;对于发生顺时针剪切的剪切带,主应变轴偏转角的均值主要为负值,而对于发生逆时针剪切的剪切带,主应变轴偏转角的均值以正值为主。     

不同变形失稳模式下隧道围岩稳定性判据研究的进展

计算机模拟隧道围岩失稳研究的困难在于失稳模型不清楚、失稳判据不明确。本文综合分析了目前国内外应用的隧道稳定性判别方法,总结了不同变形失稳机制下的计算模型及其量化的破坏判据,分析了各种方法的研究进展和存在局限,为围岩稳定性数值模拟提供快速、合理、可操作的判据。     

深埋岩溶隧道掌子面突水灾害可靠性上限分析∗

旨在研究高地应力、高压富水条件下深埋岩溶隧道掌子面突水灾害的可靠性问题。基于已有研究成果,考虑高地应力以及高压富水地质条件,采用极限分析上限法构建了以剪切破坏为主的深埋岩溶隧道掌子面三维防突机制。根据虚功率原理建立突水破坏过程中的能量方程,利用Hoek-Brown 强度准则求解了高地应力以及高压富水条件下防止掌子面突水所需要的支护力上限解。在极限破坏状态下,根据掌子面上施加的支护力与突水破坏时的围岩压力构建极限状态方程,建立了深埋岩溶隧道掌子面防突可靠度模型,并采用响应面法计算了掌子面发生突水灾害的失效概率。分析了水平地应力、溶腔水压力、岩体强度参数以及隧道洞径对支护力、潜在破坏长度的影响规律,给出了满足不同容许失效概率下深埋岩溶隧道预防突水灾害所需要的最小支护力以及在有限支护效应下能抵抗的最大破坏长度。将研究成果应用于实际工程案例中,与已有研究、现场结果相比较,验证了本文计算结果的有效性,可为今后类似深埋岩溶隧道的防突问题提供理论指导和参考。     

考虑岩体各向异性强度的隧洞围岩稳定性分析

岩体中普遍存在着大量的节理裂隙等不连续面,这些不连续面对岩体的强度有着重要的影响。以具有较密集成组节理的岩体为研究对象,以Drucker-Prager准则和裂隙应力张开准则为基础,分别考虑岩块和节理面各自的物理力学性质,研究岩体在各向异性强度准则下的理论方法,并建立了岩体各向异性本构模型,用以描述节理岩体在强度方面的各向异性及其破坏特征。将该模型用于隧洞围岩稳定性分析工程中,研究表明,该模型可以较好地解释节理岩体围岩的破坏特征,与工程实践相符。研究成果对我国地下石油、核废料储存以及深埋地下工程的设计与施工具有重要的指导意义。     

九华山隧道开挖爆破对古城墙影响的计算分析

计算分析了南京九华山隧道工程开挖爆破对明古城墙的影响。在对国内外流行的安全标准进行比较的基础上,指出了采用频率与震速结合的爆破震动安全标准的必要性,且采用德国爆破震动安全标准作为计算九华山隧道开挖爆破对古城墙的影响的判据;分别利用目前流行的震动速度和震动频率公式进行计算,判断开挖爆破对古城墙的影响程度及影响范围。结果表明,爆破震动对古城墙安全构成威胁,并计算得出开挖爆破最小安全距离。     

颗粒体材料在压缩位移控制加载条件下的力学行为及界面参数对它的影响

采用连续介质快速拉格朗日分析程序(FLAC),模拟了颗粒体材料在压缩位移控制加载条件下的力学行为及破坏过程.颗粒及它们之间的界面被离散为尺寸相同的单元.颗粒被视为各向同性的弹性材料,而界面被视为弹性-应变软化-理想塑性材料.界面完全破坏之后,其内聚力为零而内摩擦角不为零,因而适于模拟颗粒之间的摩擦行为.在双轴压缩条件下...     

隧道火灾时空温度场小波分析

提出了分布式光纤温度传感系统(DTS)时空温度场小波分析方法,完成了监测系统报警算法的设计;应用该算法,对隧道内弱风和强风条件下的模拟实验隧道火灾时空温度场进行了瞬态识别;基于隧道温度场监测数据时空矩阵的提取,对隧道的温度场变化规律做出了分析,并选取Ribo3.7小波基成功识别出隧道火灾时空温度场的奇异点.实验结果表明...