岩土力学中的塑性流动仿真与分析

为了保证岩土工程建设的安全性,在施工时要求具有特定地基和结构加固。实地测试的成本极高,因此仿真就显得非常实用,甚至必不可少。人们开发了很多数值模型来深入研究土壤行为。在这里,我们将向您介绍 COMSOL Multiphysics 中用于研究土壤的运用最为普遍的模型,及对隧道开挖实例进行分析。

岩土工程快速入门

建筑界普遍存在这样一个趋势:海上结构物建造的水域越来越深;建筑物之间的距离越来越近;海上风力发电机建造在离海岸很远的深海中,这使其可能面临着极其严苛的负载条件。因此,近几十年来,岩土工程师开发了多种数值仿真来应对这种建筑趋势以确保建筑的安全性。

岩土力学中的塑性流动仿真与分析的图1
“Paris Metro construction 03300288-3″。已获 Public domain 许可,通过 Wikimedia Commons 共享。

塑性与岩土材料

塑性是指材料能稳定地发生永久变形而不破坏其完整性的能力,金属、土壤、岩石、混凝土等材料便具有这样的特性。当造成弹性形变的应力上升到达一个特定的应力级别——屈服应力时,材料开始产生塑性形变。
弹/塑性行为是与路径相关的,应力取决于材料的之前的变形行为。因此,塑性模型通常与应力变化速率直接关联,而非应力和塑性应变。整个行业中应用最为广泛、最著名的塑性模型是以 von Mises 屈服面为基础的,该模型中塑性流动不因压力的大小而改变。因此,屈服条件及塑性流动只以偏应力张量为基础。
然而,因为分析土壤物质时需考虑摩擦和膨胀的影响,所以该模型对此类材料无效。让我们来看看该如何解决这个问题,并简单介绍一下 COMSOL Multiphysics® 仿真软件中不同的土壤塑性模型。

土壤及岩石的塑性

对于土壤和岩石等材料,摩擦和膨胀的影响是不可忽略的。众所周知,这类材料对压力非常敏感,当施加压力时会产生不同的拉伸和压缩行为。因此,上文所述的 von Mises 模型不适用于这类材料,人们转而通过求解屈服函数来研究摩擦材料的行为。
让我们通过如下图所示的滑块来阐释这类材料的摩擦行为和塑性流动。
岩土力学中的塑性流动仿真与分析的图2
如图,向滑块施加一个垂直载荷  岩土力学中的塑性流动仿真与分析的图3  及一个切向载荷  岩土力学中的塑性流动仿真与分析的图4 。假设滑块被放置在一个带有静摩擦系数为  岩土力学中的塑性流动仿真与分析的图5  的表面上,根据库仑定律,滑块在滑动前可以承受的最大力为  岩土力学中的塑性流动仿真与分析的图6  。因此,当满足以下条件时,滑块开始滑动:
(1)
岩土力学中的塑性流动仿真与分析的图7
滑动方向为水平方向。对于切向载荷来说,当  岩土力学中的塑性流动仿真与分析的图8  时,滑块不会滑动;当  岩土力学中的塑性流动仿真与分析的图9  时,滑块将会向外加载荷  岩土力学中的塑性流动仿真与分析的图10  的方向滑动。第一个土壤塑性模型是在 Mohr-Coulomb 准则的基础上开发的,它是对获得连续材料和多轴应力状态方法的一种概括。它被定义为:当任何平面内的剪应力与平均正应力的组合达到临界条件时,材料开始屈服甚至断裂。该临界条件如下所示:
(2)
岩土力学中的塑性流动仿真与分析的图11
在这里, 岩土力学中的塑性流动仿真与分析的图12  是剪应力, 岩土力学中的塑性流动仿真与分析的图13  是垂直应力, 岩土力学中的塑性流动仿真与分析的图14  是内聚力,代表垂直应力为 0 时的抗剪强度, 岩土力学中的塑性流动仿真与分析的图15  是出自著名的库仑摩擦模型的内摩擦系数。该方程代表 Mohr 平面内的两条直线。当三个 Mohr 圆都位于这两条直线之间时,应力状态为安全;而当这三个圆中的一个与这两条直线相切时,即为临界状态(材料开始屈服)。
岩土力学中的塑性流动仿真与分析的图16
Mohr-Coulomb 屈服行为。Mohr 圆以主应力  岩土力学中的塑性流动仿真与分析的图17  、 岩土力学中的塑性流动仿真与分析的图18  及  岩土力学中的塑性流动仿真与分析的图19  为基础。正如你看到的,当其中一个圆与屈服面相切时,材料便会发生屈服。
如上图所示,应力状态按照 
岩土力学中的塑性流动仿真与分析的图20  
  岩土力学中的塑性流动仿真与分析的图21  
这两个公式进行计算。因此,屈服准则及方程 2 可以被广义表达为以下形式:
(3)
岩土力学中的塑性流动仿真与分析的图22
它甚至可以被看作是以库仑摩擦为基础的更通用系列准则的一个特例,可以表达为以 应力张量的不变量 为基础的方程式:
(4)
岩土力学中的塑性流动仿真与分析的图23

岩土力学中的塑性流动仿真与分析的图24
Mohr-Coulomb 屈服函数表示法。

Mohr-Coulomb 准则对主应力空间内的六棱锥进行了定义,这为直接对该准则进行分析提供了很大的便利。但是,由于尖角(例如屈服面的法线在尖角处是未被定义的)的存在,很难对本构方程从数值计算的角度进行处理。
为了避免与尖角有关的问题, 另一个屈服准则——Drucker-Prager 屈服准则被开发了出来,它是通过对 von Mises 屈服准则进行修改而得到的,这一准则不仅将库仑摩擦纳入考虑范畴,还建立了与流体静压的依赖关系:
(5)
岩土力学中的塑性流动仿真与分析的图25
这个方程表示主应力平面内的平滑圆锥,而不是六棱锥。如果系数  岩土力学中的塑性流动仿真与分析的图26  和  岩土力学中的塑性流动仿真与分析的图27  与Mohr-Coulomb 准则中的系数互相匹配,即如下所示:
(6)
岩土力学中的塑性流动仿真与分析的图28
Drucker-Prager 屈服面会穿过 Mohr-Coulomb 六棱锥的内部或外部顶点,这取决于符号  岩土力学中的塑性流动仿真与分析的图29  为正还是为负。 塑性流动 方向来源于所谓的“塑性势”,它既可以是与塑性势的方向相同的,即与塑性相关;也可以不同,即与塑性无关,发生屈服现象(屈服函数)。我们可以据此开发很多不同的非关联流动法则。
在 Drucker-Prager 模型中,可以利用关联法则使体积塑性流动不为零。因此,在施加压力的情况下体积会发生变化。然而,这与许多土壤材料,特别是颗粒材料的行为是相矛盾的。相反,可以在塑性行为是等体积的(体积守恒)情况中使用非关联流动法则,这能更好地反映颗粒材料的塑性行为。

岩土力学中的塑性流动仿真与分析的图30
Drucker-Prager 屈服函数表示法。

COMSOL Multiphysics 中土壤塑性的非关联流动法则

接下来,我将向您展示如何借助 COMSOL Multiphysics 对土壤塑性使用非关联法则。针对软件的任意一种塑性模型,都可以使用非关联塑性法则。
如果您选择 Mohr-Coulomb 模型,有两种不同的基本方法可解决非关联塑性。由 Drucker-Prager 模型和Mohr-Coulomb 屈服函数可得到相同的塑性势,但是其静压轴斜率不同,即膨胀角取代了摩擦角(如以下截图所示)。
此外,当使用 Drucker-Prager 模型与 Mohr-Coulomb 准则进行匹配时,可以轻易地将膨胀角与您所需的非关联法则进行匹配。例如,可通过膨胀角为零来求解上述非关联法则。
最后要说明的一项要点是,我们开发了一种名为椭圆端盖的实用特征来避免材料超出特定压力范围的非物理行为。事实上,真实材料维持弹性变形的压力范围并不是无限的。因此,我们可以使用 COMSOL Multiphysics 中的椭圆端盖特征解决这个问题。

岩土力学中的塑性流动仿真与分析的图31
土壤塑性特征设置窗口。

让我们利用目前学到的所有知识来分析隧道开挖实例。同时,我们也将借此机会找出上述不同特征所产生的效果。

隧道开挖实例

为避免建设过程中的塌方事故,工人需对隧道进行必要的加固。隧道开挖过程的模拟,对所需加固量的预测非常重要。
以下模型用于对隧道开挖过程中的土壤行为进行模拟。对不同土壤模型中的地表沉降(即沿自由地表的垂直位移)及塑性区域进行计算和对比。我们将要用到的几何结构如下图所示。为了让我们的模型更加逼真,我们可以使用无限元来扩大土壤区域,同时保持足够小的计算区域,以便在相对较短的时间内找出解决方案。

岩土力学中的塑性流动仿真与分析的图32
该几何结构由 100 米深和 100 米宽的土壤层及另加 20 米的无限元构成。一条直径为 10 米的隧道位于离对称轴 10 米远、离地表 20 米深的位置。

首先,需要在隧道开挖前在土壤中加入地应力。然后,我们可以计算在移除与隧道对应的土壤后的塑性行为。地应力须被并入第二步。这些操作可以直接在 COMSOL Multiphysics 中进行设置。
我们可以从添加一个计算地应力的固定步骤开始。然后,在第二步中添加土壤可塑性特征,并进行相同的计算步骤。最后,我们可以计算出结果。为了在第二步中加入预应力,我们可以在固体力学界面中加入预应力和应变特征,如下图所示。
岩土力学中的塑性流动仿真与分析的图33
预应力和应变特征用于在隧道开挖时将第一步中的地应力作为预应力加入第二步中。变量 solid.sx 及 solid.sxy为应力张量的 x-分量及 xy-分量,其他参数也是如此。
第一幅图显示了由第一步计算得出的地应力。这些应力来自重力荷载。

岩土力学中的塑性流动仿真与分析的图34
隧道开挖前土壤中的 von Mises 应力。

第二幅图显示了隧道开挖后的应力分布。地应力由第一步计算得出。请注意,如所预见的,隧道周围的 von Mises 应力及隧道形状的变形将会增加。

岩土力学中的塑性流动仿真与分析的图35
隧道开挖后土壤中的 von Mises 应力。

如上所述,在移除隧道区域时会增加一个塑性特征,且土壤会经历一个塑性行为。如下图中带有关联塑性流动的 Drucker-Prager 模型所示。塑性区域集中在隧道周围。对此区域的分析有助于深入研究导致土壤更容易发生变形的情况。因此,这有助于我们得出避免坍塌所需的加固量和所需的隧道形状。
岩土力学中的塑性流动仿真与分析的图36
隧道开挖后的塑性区域。
我们对四种不同情况下的隧道开挖进行了仿真,以比较上文所述的不同土壤模型,以及研究内聚力对土壤行为的影响。其结果以地表沉降作为标准。
接下来,我们用一个一维绘图来展示以下结论:内聚力越低,变形越大。需要注意的是,Mohr-Coulomb 模型中的土壤比 Drucker-Prager 模型中的土壤更坚硬。带有零膨胀角的非关联准则防止了土壤在压力下膨胀,因此地表沉降愈发严重。

岩土力学中的塑性流动仿真与分析的图37
不同塑性模型与材料性质下的地表沉降比较。


本文来自 :COMSOL 博客

COMSOL岩土力学流动仿真

岩土力学中的塑性流动仿真与分析的评论3条

岩土力学中的塑性流动仿真与分析的相关案例教程

最近将J2弹塑性用户材料子程序UMAT的编写的基本理论和编程实践进行了一些总结,记录一下留存。弹塑性材料主要包含屈服条件,流动法则,硬化准则。 屈服函数主要是表征屈服条件,一般用F表示,表明应力满足某种关系时材料到达屈服,进入塑性。常见的有Mises屈服,tresca屈服,Drucker-prager屈服,Mohr—Coulomb屈服等。如果以主应力分量建立笛卡尔坐标系,则这些屈服条件在坐标系中可
abaqus混凝土坍落度试验 此示例说明了在 Abaqus 中使用扩展的 Drucker-Prager 塑性模型解决涉及有限变形的问题。 Abaqus 提供了 Drucker-Prager 类的三种不同的屈服标准。在所有三个中,屈服函数都取决于材料中的围压和偏应力。最简单的是子午线 (p-q) 平面中的一条直线。其他屈服标准是子午平面中的双曲曲面和一般指数曲面。 “扩展的 Drucke
1 引言 在过去的一个笔记中,曾经提及过Comsol Multiphysics多物理场数值模拟软件(最新发布的工业软件(采矿,结构和数值模拟),但这个软件在岩土工程中很少用到。根据调查,GeotechSet数据集内仅包括了不到10篇相关文献。本笔记简述了Comsol的岩土力学模块。 2 Comsol岩土力学模块 Comsol的岩土力学模块(Geomechanics Module)是结构力学模块的附
微信公众号:【OpenSEES抗震笔记】 为STKO for OpenSEES软件官方公众号。公众号致力于成为结构抗震专业学习平台,分享抗震相关的动力学理论,实用编程,以及常用软件技巧。特别针对基于OpenSEES的非线性分析和相应理论以及二次开发。 00 小编寄语 Angus' opnion 非常开心能够收到冯博士使用STKO 助力OpenSEES在边坡稳定研究的心得和感悟。边坡稳定是Angus
1. 引言 岩土材料不同于钢和混凝土材料,在选择本构模型时更困难。岩土材料有三个特点,这三个特点导致在数值计算时会产生特殊的问题。 第一个问题是物理不稳定性。如果岩土材料破坏时出现软化行为,就会出现物理不稳定性。当物理不稳定性发生时,部分材料储存的能量会以动能形式释放,这会导致解答可能无法收敛。 第二个特点是非线性材料的路径依赖性。大多数地质力学系统有无数的解决方案可以满足描述系统的平衡、兼容性和
影响力
粉丝
内容
获赞
收藏
    项目客服
    培训客服
    3 10