FLAC3D在边坡稳定分析中的应用.docVIP

——————————————————————————————————————————————— FLAC3D在边坡稳定分析中的应用 第６期（总第３３７期） ［文章编号】１００９—２８４６（２０１０）０６－００３９—０３ 吉妹水利 ２０１０年６月 ＦＬＡＣ３Ｄ在边坡稳定分析中的应用 王学志１，翟诚１一，毕重１，朱朝艳１ （１．辽宁工业大学土木建筑工程学院，辽宁锦州 １２１０００； ２．辽宁工程技术大学力学与工程科学学院，辽宁阜新１２３０００） 【摘要】 针对极限平衡法对边坡进行稳定分析时不能解决边坡分析的变形问题。利用ＦＬＡＣ３Ｄ软件内嵌的强度折减法对 某一土质迫坡进行数值模拟。通过模拟可以得到边坡的安全系数以及逾坡土体的应力场、住移扬的变化规律，得出最大住移所在的位置。结果表明，ＦＬＡＣ３Ｄ在解决边坡稳定分析的变形方面有着较大的优势。 【关键词】 ＦＬＡＣ３Ｄ；强度折残法；安全系数；边坡稳定分析 ＂ＩＶ 【中圈分类号】 ８６１ ［文献标识码】 Ｂ 维土体、岩体或其他材料的力学行为，尤其是达 １ 引言 目前边坡稳定性分析的方法有很多，最常用 到屈服极限时的塑性流变特性。这里，在计算时 是利用ＦＬＡＣ３Ｄ软件内有一条自动查找安全系数的命令一ｓｏｌｖｅｆｏｓ［２】，该命令就是对内聚力ｃ和 内摩擦角‘ｐ进行折减，直至边坡处于临界破坏 的是极限平衡分析法，包括瑞典圆弧法、ｂｉｓｈｏｐ法、Ｊａｎｂｕ法等。极限平衡方法的特点是，只考 状态来确定安全系数，其实质就是强度折减法，该命令通过内插逼近的办法来确定安全系数。其 计算原理如下： 虑静力平衡条件和摩尔一库伦破坏准则，通过分 析土体在破坏时的力的平衡来求解问题，计算时将边坡土体当做刚体进行分析，这种方法并不能解决边坡分析的变形问题。ＦＬＡＣ３Ｄ软件的基本原理是三维快速拉格朗日法，它是基于三维显式有限差分法的数值分析方法。用来模拟岩土或其他材料的三维力学特性。将计算区域分为若干四面体单元，每个单元在给定的边界条件下遵循指定的线性或非线性本构关系。如果单元应力使得材料屈服或产生塑性流动，则单元网格可以随着材料的变形而变形，可以对连续介质进行大变形分析，相对局限于小变形假定的其他方法更适合岩土工程的破坏问题研究…。这里利用ＦＬＡＣ３Ｄ 滑坡的安全系数定义为实际抗剪强度指标与滑坡发生破坏临界状态时虚拟的折减强度指标的 比值： ， ，（素＋詈ｔａｎ９）以 ＼ｆｄＡ Ｌ］＿，其中ｃ，２素，妒７＝ ｌ（ｃ’＋ｔｒｔａｎｔｐ’）ｄＡ 』ｒｄＡ ，’ “ 欹件内嵌的强度折减法来进行数值模拟，对边坡 的稳定性给予评价。 ８ｒｃｔａｎ（警）。 通过逐步调整系数，得到不同的ｃ’、妒’，代人ＦＬＡＣ３Ｄ反复迭代计算，直至坡体达到临界状态，坡体达到临界状态的Ｋ值即为边坡安全系数‘”。在莫尔一库伦屈服准则下解得不收敛作 为破坏标准。２．２屈服准则 ２计算原理 ２．１ ＦＬＡＣ３Ｄ软件内嵌强度折减法 ＦＬＡＣ３ Ｄ是由美国ＩｔａｓｃａＣｏｎｓｕｌｔｉｎｇＧｒｏｕｐｌｎｅ 开发的的三维显式有限差分程序，它可以模拟三 ［收稿日期】２００９一０４一１４ 莫尔～库仑强度准则是目前工程中应用最为 【作者简介】王学志（１９７６一，。男．博士．剐教授，研究方向：混凝土断裂研究。 【注】国家自然科学基金贵助项目辽宁省教育厅科学技术研究项目。 一３９— 万方数据 吉林水利 广泛的破坏准则，其表达式为¨Ｊ： ＦＬＡＣ３Ｄ在边坡稳定分析中的应用王学志等２０１０年６月 ＝１８ｋＮ／ｍ３。地面以下为饱和粘土，容重＾ｙ＝２０ｋＮ／ｍ３。坡比是１：１，坡高是１０ｍ。Ｘ方向的 ／＝争ｌｓｉｎ妒＋（ｃｏｓＯ，一去ｓｉｎ驴ｉｎｃｐ）压一 ｃｃｏｓ‘，ｏ ２ 距离是４０ｍ，Ｙ方向的距离是３０ｍ，Ｚ轴负方向 的距离是１０ｍ。黏聚力ｃ＝１０ｋＮ／ｍｚ，摩擦角‘ｐ＝３０。，体积模量为２００ＭＰａ，剪切模量１００ＭＰａ。建立的力学模型如图１所示。模型共有１９ ６５０ ０ 式中：，２，以为应力张量第一不变量和应力偏 量第二不变量，０。为应力罗德角，Ｃ、９为黏聚力和内摩擦