基于时间增量步的边坡稳定性计算精度分析.pdfVIP

112 公 路 与 汽 运 2014年9月 基于时间增量步的边坡稳定性计算精度分析 胡勇军，郑智雄 ，曹志伟 (长沙理工大学 土木与建筑学院，湖南 长沙 410004) 摘要：在有限元强度折减法的基础上，利用ABAQUS中提供的时间增量步，提 出一种计算边 坡稳定性的方法——基于时间增量步的有限元强度折减法；对该方法应用于边坡稳定性分析 的可 靠性进行验证，利用ABAQus对时间增量步控制的优势，对不同坡角下边坡稳定安全系数计算的 可靠性及不 同网格划分形式对其计算精度 的影响进行分析 ，结果表 明该方法求得 的安全 系数有 效、可靠，采用坡面局部加密网格计算安全系数具有很高的计算精度，并能提高计算效率。 关键词 ：公路 ；时间增量步；有限元强度折减法 ；边坡稳定性分析 ；安全系数 ；临界滑移面 中图分类号：U416．1 文献标志码 ：A 文章编号：1671—2668(2014)O5一O112一O3 随着数值计算 的不断发展 ，有 限元强度折减系 材料参数粘聚力 、内摩擦角随增量步时间变化的 数法在边坡稳定性分析中得到广泛应用 。与传统极 函数关系可通过ABAQUS中的场变量建立起来。 限平衡法相比，它不需要任何假设便能 自动地求得 (2)折减系数与增量步时间的关系。有限元强 任意形状的临界滑移面及对应 的安全系数，还可 以 度折减法中的折减系数F定义如下： 反映边坡失稳及塑性 区的开展过程 。但有 限元强度 Co 一 i ci一F (33)) 折减法在操作过程中要将折减系数手动输入有限元 软件 inp文件 中，并不断进行试算，其操作过程较繁 tango tan 一 — an~p；一—百一 L(44) 琐。为了简化计算过程，在有限元强度折减法基础 由式 (1)～(4)可得出： 上通过 ABAQus设置岩土体粘聚力与内摩擦角为 1 场变量中增量步时间 t的函数 ，从而通过控制增量 F一 · (5) n— bt 步时间t来实现对强度参数 的自动折减 ，经一次计 由有限元强度折减法安全系数的定义 ，按照不 算即可求得边坡的安全系数，提高计算效率。 收敛判据得到相应的折减系数即为安全系数。 1 基于场变量 中时间增量步的有限元强度 (3)可调参数 a 、b数值的确定 。ABAQUS软 折减法 件中增量步时问取值范围为 0≤t≤1，假设取 a、b 均为正值，由于折减系数要满足 F≥0，故 只需满足 (1)建立材料参数与增量步时间t的关系：