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

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

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

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

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

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

地震作用下深埋隧道围岩形变的数值模拟研究

隧道围岩受地震荷载作用影响产生开裂变形,给交通、安全与经济带来极大影响,为提升隧道安全稳定性能,以某市隧道为对象,研究地震作用下深埋隧道围岩形变的数值模拟。综合考虑地震动作用与地震波扩散特性,利用FLAC软件构建三维动力模型;通过边界设置解决地震波反射问题,基于物理力学参数与地震荷载条件,通过相对变形法研究地震条件下深埋隧道围岩变形特性、能量聚集特性与安全性能。研究结果显示:深埋隧道围岩受地震荷载作用影响,形成挤让、椭圆化的变形走向;变形达到最高值时,应变能量密度高于2350 J/m^3的区域集中在隧道Ⅰ洞左侧围岩边墙和中夹岩柱上;各监测点安全系数呈非对称性,隧道右侧下角区域安全系数低于国家相关标准,计算结果与实际测量结果一致。     

油页岩地下开采巷道围岩的流变特性及工程应用

为研究油页岩巷道围岩的流变特性,采用岩土SJ-1B三轴蠕变仪对油页岩进行了三轴蠕变试验。根据试验结果,建立了油页岩非线性蠕变方程。根据上述非线性蠕变方程,以吉林桦甸油页岩矿一井运输大巷为工程背景,分析了围岩的应力场和位移场,并对不同支护强度和应力状态下的蠕变变形进行了系统的分析。研究结果表明,油页岩蠕变具有非线性,控制油页岩巷道围岩过量变形的根本途径是改变围岩的应力状态,以及适当提高锚杆、锚索的初始预应力。本文为有效控制油页岩开采时巷道围岩的有害变形提供了理论依据。     

基于岩体波速和Hoek-Brown准则的初始地应力场参数计算模型研究

在初始地应力场参数计算中引入岩体波速,通过基于岩体波速的Hoek-Brown准则建立初始地应力场中水平地应力σ_H的取值范围Δ_(σ_H)计算模型,并通过实测地应力数据验证了Δ_(σ_H)计算模型的合理性。在Δ_(σ_H)计算模型的基础上,推导出了σ_H/σ_V(水平地应力/垂直地应力)的取值范围Δ(σ_H/σ_V)计算模型,并通过Δ(σ_H/σ_V)计算模型从岩体稳定性的角度对σ_H/σ_V在地壳表面和深部的分布现象进行合理解释。通过Δ(σ_H/σ_V)计算模型给出的σ_H/σ_V的取值下限(σ_H/σ_V)_(min)和取值上限(σ_H/σ_V)_(max)算出侧压系数λ,并给出了侧压系数λ与深度h之间的函数关系λ=ah~b+1。     

不同场地条件下水下隧道周围土体位移变化的数值分析

以拟建的哈尔滨松花江隧道为背景,考虑土体—衬砌的材料非线性特性和相互作用,按平面应变假定,建立在不同场地条件下的若干接近实际的有限元计算模型,应用有限元计算软件ANSYS进行数值分析,得出不同场地条件下隧道周围土体位移的分布规律,遴选出了影响水下隧道周围土体变形的部分关键因素,为合理确定隧道覆土层安全厚度提供了一定的理论依据。     

半横向通风方式隧道内污染物分布的数值仿真

某城市公路隧道采用半横向通风方式进行火灾排烟,其排风孔位于隧道侧壁顶部。针对该城市公路隧道,运用CFD数值仿真方法建立1：1的三维数值仿真模型,模拟不同工况下,开启不同形状排风孔时,隧道内的风速、风压以及污染物浓度的分布情况。研究表明,排风孔所在位置的隧道横断面上,形成一个左侧区域负压较高、压强变化梯度大而右侧区域负压较小的＂排风横断面＂;随着排风孔开口高度的增加,排风横断面上负压分布向右侧扩大,右侧区域风速增加,有利于右侧区域内污染物的稀释,而左侧区域竖直方向风速减小,左下角区域风流方向紊乱,左侧区域污染物不能得到有效稀释,浓度较高。     

平面P波入射下深埋圆形复合式衬砌隧道抗减震机理研究

基于Fourier-Bessel级数展开法,研究深埋圆形三层复合式衬砌洞室在平面P波入射下的动应力集中问题,并给出三层衬砌洞室动应力集中系数级数解析解;依托某IX度地震区管道隧道实际工程,分析不同衬砌刚度组合和厚度组合对洞室动应力集中系数的影响。研究表明:注浆加固洞室围岩和设置减震层都可以降低二次衬砌动应力集中系数;增大围岩注浆区弹性模量和厚度,有利于减小衬砌动应力集中系数,最优围岩注浆区厚度为1倍洞室净空半径;减震层弹性模量降低,减震层厚度增大,二次衬砌动应力集中系数变小,减震层弹性模量宜低于围岩弹模1/20,最优减震层厚度宜取1/50的洞室内净空半径。最后针对实际管道隧道抗减震技术,考虑围岩稳定性,提出"围岩-加固圈-减震层-衬砌"新型减震结构,分析结果表明:对比其他三种抗减震措施,新型减震结构的减震效果最好。     

Characteristics of a Short-Term and Imminent Earthquake Precursor Field in the Northern Part of North China and Numerical Simulation of Source Process

According to studies of more than 20 earthquakes with MS≥5.0 in North China, seven features of the short-term and imminent earthquake precursors have been summarized in this paper. At the same time, taking the short-term and imminent earthquake stage as the physical process of a source's medium softening and fault creep, we calculated temporal variation of mean stress, maximum shear stress, body strain, and pore pressure in some certain points (supposed stations) in the source area and its adjacent area by using an anisotropic and nonlinear source model and a finite element method. According to an analysis of these theoretical curves, we conclude that the short-term and imminent earthquake precursors have such characteristics as complex shapes, exponential growth of the precursor number with tune, and precursors' migration from the outside area to the source area, which to a certain extent reveal the cause of the characteristics of the short-term and imminent earthquake precursor field.