路堤荷载下软土路基沉降的数值模拟.pdf

路堤荷载下软土路基沉降的数值模拟 程 志 刘劲松 摘 要 ：通过应用ABAQUS软件内部的有限元程序对路堤荷载下的沉降进行了模拟，初步预测了路堤的沉降量，从而对 指导路基施工以及后期的沉降观测具有一定的实际意义。 关键词：软土地基，沉降，数值模拟 中图分类号：U416．1 文献标识码：A 软土路基的沉降多采用现行公路软基规范中推荐的方法进 的莫尔圆提出的，破坏线是与这些莫尔圆相切的直线，如图 1所 行计算。近年来，随着有限单元法技术的不断成熟，越来越多的 示，Mohr—Coulomb的强度准则为：r=C—atan~。其中，r为剪切 本构模型以及计算方法可以应用于实际的工程，并且得到了很好 强度 ；为正应力；c为材料的粘聚力；≠为材料的内摩擦角。在 的实际效果 1【,2J。 平面上 ，Mohr—Coulomb模型为等边不等角的六边形，屈服面存在 在ABAQUS中进行岩土工程问题的弹塑性计算时，鉴于常 尖角。ABAOUS采用的本构模型是经典 Mohr—Coulomb屈服准则 规地质报告所提供的参数 ，我们常常选用 Mohr—Coulomb模型和 的扩展，该模型的流动势函数在子午面上的形状为双曲线，在 Drucker-Prager模型本构关系。本文分别应用两种模 型进行模 平面上没有尖角，因其势函数完全光滑，确保了塑性流动方向的 拟，从而考虑结果的合理性以及两者的差异性。 唯一性 。 1 弹塑性模型l3J 1．2 Drueker-Prager模 型 1．1 Mohr—Coulomb模 型 ABAOUS对经典的Drucker—Prager模型进行了扩展，扩展的 Mohr-Coulomb屈服准则假定 ：作用在某一点的剪应力等于该 Drucker—Prager模型的屈服面在 平面上不是圆形 的，屈服面在 点的抗剪强度时，该点发生破坏 ，剪切强度与作用在该面的正应 子午面上包括线性模型、双曲线模型和指数模型。本文采用线性 力成线性关系。Mohr—Coulomb模型是基于材料破坏时应力状态 模型。在子午面上 ，线性Drucker—Prag口模型的屈服轨迹如图2所 示 ，屈服准则的表达式为：F=t—ptanfl—d==0。 ／ ／ h ^／ 、＼、 图 3 路基有限元网格 d 1 1 J3 d3 d 1十 3 d — —_ d ～ … 一 !一～ 图 1 Mohr—coul0mb破坏模型 ／ j 户 图 4 Mo~ oulomb模型竖向沉降(单位 ：m) 图2 Drucker-Prager屈服轨迹 1．3 参数转换 在岩土工程中经常遇到平面应变问题，例如长隧道、堤坝、边 坡和路坡等等。因此，有必要建立适于平面应变流动及破坏的模 型参数，下面给出的转换关系是针对线性Drucker—Prager模型展 开的，其同样适用于高围压应力水平下的双曲线模型。平面应变 问题 中Molar—Coulomb模型和线性 Drucker-Prager模型参数的转 换关系如下式 ：相关联流动时，由 = 得： √3sin~ d 43cos~ ． 图5 Dru