4模型边界条件及初始应力场的合理确定.pptVIP

第四节 初始地应力场的确定 (二)岩体构造应力场 2．弹塑性状态 由式(4-18)，Hcr＝0的条件为： 即 时岩体都为弹性，岩层中应力按式(4-13)～式(4-16)计算。 岩体处于由弹性进入塑性的临界状态，塑性区深度为零(在地表)。 浅部岩体要进入塑性，此时先按式(4-18)计算塑性区深度Hcr。 第四节 初始地应力场的确定 (二)岩体构造应力场 2．弹塑性状态 若岩层处于Hcr深度以下，按式(4-13)～式(4-16)计算岩层中应力；若岩层处于Hcr以上，则按下列式子计算岩层中地应力值σ1、σ2、σ3、σT(考虑岩体为理想塑性体)： 第四节 初始地应力场的确定 (二)岩体构造应力场 2．弹塑性状态 (2)在σ1、σ3水平方向，σ2在垂直方向。此时在Hcr处有： 代人式(4-6)，整理得Hcr的表达式为： 由式(4-23)可讨论Hcr＝0的条件为： 第四节 初始地应力场的确定 (二)岩体构造应力场 2．弹塑性状态 考虑一般只有软岩进入天然塑性状态，软岩的岩体结构特征常数s值极小, 可忽略，于是有： 岩体都为弹性，岩层中应力按式(4-13)～式(4-16)计算； 岩体由弹性进入塑性的临界状态； 浅部岩体要进入塑性，此时先按式(4-7)计算塑性区深度Hcr。 第四节 初始地应力场的确定 (二)岩体构造应力场 2．弹塑性状态 若岩层处于Her深度以下，则按式(4-13)~式(4-16)计算岩层中应力；若岩层处于Hcr以上，则按下列式子计算岩层中的地应力值σl、σ3和σT(考虑岩体为理想弹塑性体)： 第四节 初始地应力场的确定 三、地应力值对模拟结果的影响分析 为了研究地应力值对模型计算结果的影响，根据表4-2。模拟中重点考虑不同水平应力对模拟结果的影响，改变模型的侧压系数，计算至初始平衡，对比分析不同水平应力下的巷道围岩变形规律。 2.0 10 1.0 5 1.8 9 0.8 4 1.6 8 0.6 3 1.4 7 0.4 2 1.2 6 0.2 l 侧压系数 模型号 侧压系数 模型号 表4-2 侧压系数变化表 第四节 初始地应力场的确定 三、地应力值对模拟结果的影响分析 图4-9 巷道围岩位移状况 第四节 初始地应力场的确定 三、地应力值对模拟结果的影响分析 图4-10 不同侧压下巷道顶板位移图 * 第四章 模型边界条件及初始应力场的合理确定 第一节 模型边界条件的类型及其选择 数值模型的边界条件由场变量(如应力、位移)组成，它是场变量在模型边界的指定。力学边界条件主要有两种类型：应力边界条件和位移边界条件。 第一节 模型边界条件的类型及其选择 一、应力边界 自由面是指定应力边界的一个特例。缺省时模型边界不受应力和任何限制。 Itasca软件允许在指定的范围内加线性变化的应力或力。加应力梯度的典型情况是用来再现由重力引起的随深度增加的应力影响，但要保证所加梯度与初始化应力场指定的梯度和重力加速度的值相一致。 第一节 模型边界条件的类型及其选择 二、位移边界 Itasca软件中不能直接控制位移。实际上，它们在计算中不会发生作用。为了对一条边界加给定的位移，必须在给定的时步对那条边界指定边界速度。如果要求的位移是D，在时间增量T内施加速度u，即D＝uT，这里T＝△tN，△t是时间步长，N是时步数(或循环数)。实际上，u应保持很小而N很大以便减小对模拟系统的振动。 位移边界特别是零位移边界(约束边界)是计算模型中不可缺少的，没有约束的计算模型在不平衡力的作用下将产生平动或转动。 第一节 模型边界条件的类型及其选择 三、边界元边界 边界元边界条件是一种人工边界，用来模拟各向同性、线弹性介质无限(或半无限)区域的效应，直接以边界单元求解。下面以离散单元法(UDEC)为例介绍这种方法。 该方法沿UDEC模型外边界耦合到离散单元区域，离散单元区域嵌于弹性区域中。这个区域首先受到初始地应力场。如果在离散单元区域作一些改变，如开挖，那么开挖引起的位移和最终的应力状态在离散单元与边界单元界面上必须满足连续性。 三、边界元边界 边界单元节点与离散单元区域周边上的块体角点一致。确定每个边界单元节点上产生的位移，通过节点位移矢量乘以边界单元区域的刚度矩阵直接确定产生的节点力： f = K u (4-1) 式中 K——边界单元刚度矩阵； f——诱发的节点力矢量；u——节点位移矢量。 诱发节点力与总节点力的关系为： ft = f0 + f (4-2) 式中 fo——初始节点力矢量。 这样，方程(4-2)可写成： ft = Ku + f0