精选优质文档-----倾情为你奉上摩尔-库仑模型及其在FLAC3D中的应用摘要: 本文首先阐述了塑性流动理论的增量方程，结合摩尔库仑破坏准则和拉伸破坏准则形成了FLAC3D中采用的摩尔库仑本构模型，并指出不同的应力计算值条件下的计算方法最后通过模型试验与解析方法进行对比，发现FLAC3D计算结果与简单模型下精确的解析解吻合较好，但在变性较大时逐渐出现一定偏差关键词: 摩尔-库仑模型，增量方程，流动法则，FLAC3D1. 塑性流动理论的增量方程一般情况下，破坏准则可表示为 （1） 式中，为已知屈服函数，用来判定塑性流动开始产生在主应力空间中，为一曲面，落在曲面内的应力点为弹性状态塑性状态下的应变增量可表示为弹性应变增量和塑性应变增量之和： （2）弹性应变增量和弹性应力增量的关系表示为： （3）式中，为弹性应变增量的线性方程流动法则规定了塑性应变增量向量的方向，即与塑性势面的方向垂直，表示为： （4）得到的新的应力矢量应满足屈服方程： （5）式（5）提供了一个估计塑性应变增量矢量的表达式。
将式（2）代入式（3），且考虑到为线性函数，得： （6）再将流动法则（4）代入得： （7）假定破坏函数为线性函数，式（5）可表示为： （8）式中，代表函数减去其常量值，对于位于屈服面上的应力点，，式（8）可转化为， （9）此时，定义新的应力分量为： （10） （11）根据式（11），可得： （12）综合式（9），（12），可得： （13）根据应力增量表达式（7），估算应力（11），新的应力（10）可表示为： （14）2. 莫尔库伦模型（IN FLAC3D）莫尔库伦模型的破坏包线包括两部分，一段剪切破坏包线和一段拉伸破坏包线与剪切破坏相对应的是相关联的流动法则，与拉伸破坏对应的是不相关联的流动法则 在FLAC3D中，莫尔库伦模型表示在，，主应力空间中，对应的应变分量为主应变，， 弹性增量方程主应力空间中，虎克定律的增量表达式可写为， （15）式中，和为由剪切模量和体积模量定义的材料常数。
(16)根据式（3），式（15）可改写为： （17）复合破坏准则莫尔库伦模型所采用的破坏准则为摩尔库仑准则和最大拉应力准则三个主应力为破坏准则在面表示如图1图1 FLAC3D莫尔库伦破坏准则破坏包线，在A到B上由莫尔库伦准则定义， （18）在B到C上由拉伸破坏准则定义， （19）式中，为摩擦角，为粘聚力，为抗拉强度 (20)由图1可见，材料的抗拉强度不能超过和交点对应的值，因此抗拉强度的最大值为 (21)2.3 流动法则塑性势面由两个方程来描述，和，分别用来定义剪切塑性流动和拉伸塑性流动函数为不相关联的流动法则， （22）式中，为膨胀角， （23）函数为相关联的流动法则， （24）统一的流动法可由函数定义，为和的对角线，如图2所示。
(25) (26) (27)图2 莫尔库伦模型－流动法则当由式（11）计算出的对应的应力点落在图2所示domain 1中，产生剪切破坏，应力点在相应曲线上，流动法则由塑性势方程推得如果点落在domain 2中，发生拉伸破坏，新的应力点在上，由流动法则推出2.4 塑性修正首先考虑剪切破坏的情况式（22）进行偏微分可得： （28）将，，代换，，，式（17）变为： （29）由式（14），（18）可得： （30） 且， （31）考虑拉伸破坏的情况，类似方法可得： （32）3. FLCA3D执行方法当在FLCA3D中运行莫尔库伦模型时，首先通过将由虎克定律计算出的应力增量叠加到原有应力上计算（），这时可计算出主应力如果主应力达到破坏准则，进入domain 1或者domain 2在第一种情况下，产生剪切破坏，,, 由式（30）求得；在第二种情况下，产生拉伸破坏，,, 由式（32）求得。
如果应力点在平面上落在包线内部，那么表明在这一步计算中没有出现塑性流动，新的主应力为4. 模型验证如图3所示，x，z方向主应力相等，边界条件为 （33）式中，是试样y方向的恒定变形速率，L为试样高度图3 模型实验边界条件在FLAC3D中，采用一单元进行模拟，计算坐标（0，1，0）在指定竖向位移下的竖向应力并与解析解对比，结果如图4，5经计算发现，解析解与数值解吻合较好，当点的竖向位移较大时，存在偏差5．结论将塑性理论的增量模型及摩尔库仑准则和拉伸破坏准则相结合，形成FLAC3D中采用的摩尔库仑模型针对不同的应力计算值，在FLAC3D中采用不同的处理方法最后通过模型试验与解析方法进行对比，发现FLAC3D计算结果与简单模型下精确的解析解吻合较好，但在变形较大时逐渐出现一定偏差图4 vs Sy (膨胀角100)图5 vs Sy (膨胀角00)6. 参考文献[1] Itasca Consulting Group, Inc. FLAC3D User Manuals, Version 2.1, Minneapolis, Minnesota, 2002.6[2] Wood，D.M. Soil Behavior and Criterion State Soil Mechanics. Cambridge: Cambridge University Press, 1990.[3]钱家欢,殷宗泽. 土工原理与计算. 北京:中国水利水电出版社,1996[4]郑颖人,沈珠江,龚晓南. 岩土塑性力学原理. 北京:中国建筑工业出版社,2002专心---专注---专业。