有限单元法在工程地质中的应用.ppt

(一)无厚度节理单元 古德曼(R.E.Goodman)提出的，又称“古德曼节理元”。 模拟闭合，不含充填物的节理裂隙。 单元是一段由两个直接接触的平面(模拟节理或裂隙的两个壁面)构成的四结点单元体(图4-1)。 最大主应力分布云图 最小主应力分布云图 二、支护的模拟 地下井巷围岩在重分布应力作用下将产生变形和破坏，为保证地下工程的正常使用，在施工过程中就要对围岩进行支护，支护的形式多种多样，支护和围岩之间相互力学作用过程也是十分复杂的。 利用有限单元法可以模拟这种相互作用，为选择合理的支护型式、合适的支护时间提供依据。 衬砌和锚喷是两种经常采用的支护结构形式。 锚杆的模拟可采用第二章第二节介绍的一维杆单元。 衬砌、支护的模拟过程如图4-8所示，每一步计算都以前一步形成的应力、位移作为出发点。 每一步开挖把要开挖部分的单元作为“空单元”，每一步衬砌把要衬砌的单元重新赋予衬砌材料的参数，此过程可简称为单元的“死”与“活”。 一、从地质模型到计算模型 二、岩土本构关系及参数选择 三、结果分析举例 本节主要内容： 第三节 开采沉陷及覆岩破坏的模拟 一、从地质模型到计算模型 （1）对预测区工程地质岩组柱状的分析 （2）对预测区工程地质剖面的分析 （3）对预测剖面垂直层面方向进行划分 （4）计算模型范围的划定 （5）计算边界条件的确定 从地质模型到计算模型的一般过渡步骤： 二、岩土本构关系及参数选择 岩层移动计算中应合理选择岩土的本构关系。这里的“合理”一方面是指所选择的本构关系要符合预测区岩土实验的结果和实际情况；另一方面是指所采用的本构关系要简便，参数容易取得且经济合理。 （1）横观各项同性线弹性模式 （2）各向同性线弹性模式 （3）弹塑性模式 岩土本构关系： 参数选择： 弹性模量E、泊松比μ、粘聚力c、内摩擦角 和抗拉强度σt。 结构面强度 岩体强度 岩石强度 岩体的抗剪强度介于岩块和软弱结构面的抗剪强度之间。岩体越完整，其c、 值越接近于岩块的。一般岩体的c、 值可在顺层理岩块指标和垂直层理岩块指标中选取。 岩体的抗拉强度是确定岩体是否产生拉破坏的指标。破碎岩体，取其抗拉强度为零，对于完整的岩体则采用岩块指标。由于成层岩体岩石力学性质各向异性比较明显，计算中采用采用各向异性材料分析时，可分别就平行层理和垂直层理方向选取不同的计算参数。 在计算参数的选取时应采用数理统计的方法，给参数一个概率的含义，计算结果也带有概率的含义。 同一矿区需要不断地积累经验数据，建立参数数据库。 三、结果分析举例 1-超前应力区；2-控顶区；3-冒落发展区；4-冒落岩块受压缩区；5-冒落岩块稳定区； W-垂直方向移动向量；u-水平方向移动向量 图 正常回采阶段采动影响的空间—时间关系 t(d) 初始期 活跃期 衰退期 前期变形过程（采动沉陷） 1.7mm/d V(t) W(t) V(mm/d) W(mm) 前期变形的阶段性特点示意图 采动岩体内平面应力场按照最大主应力和最小主应力的大小、方向和性质可归纳为3中类型（P92图4-13）： （1）双向拉应力区 （2）拉压应力区 （3）压应力区 采空区 煤柱 沉降中心区 沉降边缘区 沉降边缘区 沉降后地表面 沉降前地表面 图 采后地表变形分区示意图 覆岩破坏状况分为4各区域（P93图4-14）： （1）未破坏区 （2）塑性区或剪切破坏区 （3）拉张破坏区 （4）拉张裂隙区 下面举例说明其应用 30 2.20 0.80 7.35 0.32 2200 2000 充填层 42 6.02 3.54 22.43 0.22 38400 2620 40 砂质泥岩组 11 38 14.37 4.86 84.32 0.18 77300 2750 7 L7灰岩 10 40 7.83 2.33 19.97 0.24 3900 2600 3 泥岩 9 36 11.37 4.23 82.60 0.19 76100 2700 3 L8灰岩 8 42 7.62 0.93 12.97 0.30 1900 2600 10 底板泥岩组 7 25 1.20 0.64 4.86 0.38 1500 1300 7 二1煤 6 44 3.13 0.76 8.77 0.28 1700 2580 5 顶板泥岩 5 48 6.41 5.43 81.18 0.21 68400 2650 15 大占砂岩 4 45 6.62 1.93 21.60 0.22 17200 2600 40 砂质泥岩组 3 47 7.13 5.00 83.20 0.21 67500 2640 10 砂锅窑砂岩 2 32 0.71 0.91 5.50 0.37 1300 1