FLAC上机题目2012模型及参数.docVIP

FLAC3D数值模拟上机题 计算模型分别如图1、2、3所示，边坡倾角分别为30°、45°、60°，岩土体参数为： 密度ρ＝2500 kg/m3, 弹性模量E＝1×108 Pa，泊松比μ＝0.3， 抗拉强度σt＝0.8×106 Pa，内聚力C＝4.2×104 Pa，摩擦角φ＝17°，膨胀角Δ＝20° 试用FLAC3D软件建立单位厚度的计算模型，并进行网格剖分，参数赋值，设定合理的边界条件，利用FLAC3D软件分别计算不同坡角情况下边坡的稳定性，并进行结果分析。 附 换算公式： 1 kN/m3= 100 kg/m3 剪切弹性模量： 体积弹性模量： 图1 倾角为30°的边坡(单位：m) 图2 倾角为45°的边坡(单位：m) 图3 倾角为60°的边坡(单位：m) 目前，边坡稳定性分析在实际工程中主要应用的是极限平衡法，而该方法中，基于经典极限平衡理论的条分法用得最为广泛，但是它存在很多不足，不能考虑岩土内部的应力应变关系，无法分析边坡破坏发展的过程，并需要事先知道滑移面。随着计算技术的发展，有限差分法逐渐在为边坡分析中非常重要的方法之一。强度折减法则是有限元分析中比较常用的方法，它通过不断折减材料状态，从而计算出边坡的安全因数。 本文利用FLAC3D计算模型，采用理想弹塑性材料，基于强度折减法，用数值计算中常见的Mohr—Coulomb对边坡进行非线性有限元分析，得到边坡的安全系数，从而判断边稳定性。 1.有限元强度折减法的原理 建立在强度缩小有限元分析基础上的边坡稳定分析的基本原理是将边坡强度参数粘聚力和内摩擦角同时除以一个折减系数F，得到一组新的和值然后作为一组新的材料参数输入,再进行试算。当计算不收敛时,对应的F被称为边坡的最小安全系数,此时,边坡达到极限状态,发生剪切破坏,同时可得到临界滑动面所在的塑性区。 2.安全系数的物理意义 Bishop等将土坡稳定安全系数F几定义为沿整个滑动面的抗剪强度与实际抗剪强度之比，即有: （1） 上式可以进一步写为如下的形式: （2） 强度折减法分析的实质是一种弹塑性有限元分析，与通常意义的弹塑性分析的差别在于，在分析过程中不是按荷载步调整加于计算模型上的荷载，而是在分析中，不断地按下式调整岩土材料的、值。 (3) (4) 式中: 、为岩土材料的粘聚力和内摩擦角; 、为经过折减后的粘聚力和内摩擦角; F为折减系数。 通常对于折减系数F，在开始计算时可以凭计算者对实际问题的判定而确定一个适当的初值，再进行弹塑性分析。若分析过程在指定的迭代步中收敛则认为折减系数的取值太小，可适当的加大折减系数并按式(3)和式((4)对岩土材料的、值进行调整，再按前述方法进行分析。当折减系数F达到某一值时计算就可能发散，此时所对应的折减系数F即认为是边坡的安全系数。 3.屈服准则 影响土坡稳定的关键因素是土体的抗剪强度，当剪切面上的最大剪应力达到了土体能够承受的破坏极限时，土坡将发生失稳破坏，此处岩土体的本构关系采用理想弹塑性模型，屈服准则为Mohr—Coulomb屈服准则 （5） 式中：为内摩擦角，为内聚力，为应力张量的第一个不变量； 为应力偏张量的第二不变量； 为应力罗德角。 4.强度折减法的极限状的判剧 强度折减法的原理和操作简单，其关键问题是临界状态的判定，目前主要有三类判定准则。 以数值计算的收敛性作为失稳判据 （2）以特征部位位移的突变性作为失稳判据。绘制研究对象多个特征点的位移与折减系数的关系曲线，以位移曲线出现较为明显的转折点对应的强度折减系数作为该研究对象的安全系数。 （3）以塑性区的贯通性作为失稳判据。当计算结果显示结构面塑性开展区相互贯通，此时的折减系数即为坡体稳定安全系数。 5.针对本次上机的边坡稳定性求解 5.1 边坡倾角为30°时的情形 图1 x方向上位移云图 图2 x方向上应力云图 图3 位移云图 5.2边坡倾角为45°时的情形 图 x方向上位移云图 位移云图 安全系数1.14 5.3边坡倾角为60°时的情形 安全系数1.007 4