边坡渗流案例.pdfVIP

在 COMSOL Multiphysics 5.6 版本中创建 堤坝边坡稳定性 此模型基于 《COMSOL 软件许可协议》5.6 版本授权。 所有商标均为其各自所有者的财产。请参见 /trademarks 。 简介 边坡稳定性分析是预测由于各种载荷和环境条件引起的地层沉降、变形和滑移的重要 方法。在堤坝中，边坡稳定性分析对于确定大坝的安全性非常重要。 本例模型受参考资料 1 中一个示例的启发。土壤中的孔隙压力通过达西定律建模，而 堤坝的弹塑性分析采用莫尔- 库仑准则。 用于研究边坡稳定性的方法称为强度折减法，其中莫尔- 库仑材料参数是安全系数 (FOS) 的函数。降低材料参数会导致土壤剪切强度降低，最终导致堤坝不稳定。这种现 象会在一定的荷载组合下导致路堤失稳。有关该方法的更多详细信息请参见参考资 料 1。 模型定义 图 1 显示堤坝的横截面。长度 L 、L 和 L 分别为 24 m 、5 m 和 24 m ，堤坝高度 L 为 1 2 3 4 12 m。水位为 10 m ，可能渗流高度为 4 m 。堤坝总宽度为 L 1 + L2 + L3 。 . L1 L2 L3 水库水位 L4 渗流高度 H w H s 图 1 ：堤坝示意图。 本例中使用平面应变近似在二维模式下对堤坝进行建模。重力和静水压力的影响的影 响也包含在内。莫尔- 库仑模型的材料属性根据安全系数参数 FOS 进行参数演化。参 数化研究逐渐增加 FOS 参数，从而降低了每个参数步骤中的土壤强度。对于 1.915 以 上的参数值，模型不收敛，这表明堤坝失稳。 莫尔- 库仑屈服函数及相关塑性势为 F Q m J2 + I 1 – k (1) 2 | 堤坝边坡稳定性 其中 sin   k C cos  (2) 3 参数化内聚力 C 和内摩擦角  为 c tan u tan s C  atan p 0+ atanp 0 (3) FOS FOS FOS   其中 c 是材料内聚力， 和 分别是非饱和图及饱和土的内摩擦角，p 为达西定律给 u s 出的孔隙压力。 结果与讨论 堤坝的压力水头如图 2 所示；在浸没壁上，它从 0 m 到 10 m 不等，而在渗流面上为 0 m 。正压力水头表示正孔隙压力，表示饱和土，而非饱和土用零压力水头表示。图中 的零压力水头线是将饱和土与非饱和土分开的潜水面位置。 对于大于 1.915 的 FOS 值，弹塑性分析不收敛；因此，进行仿真计算直到 FOS 值变为 1.915，这是由于塑性应变增加以及随后的剪切强度降低而导致路堤失稳的值。 . 图 2 ：堤坝中的压力水头。零压力水头线表示潜水面的位置。