摩尔-库仑模型与其在FLAC3D中的应用..doc

摩尔-库仑模型及其在FLAC3D中的应用 摘要: 本文首先阐述了塑性流动理论的增量方程，结合摩尔库仑破坏准则和拉伸破坏准则形成了FLAC3D中采用的摩尔库仑本构模型，并指出不同的应力计算值条件下的计算方法。最后通过模型试验与解析方法进行对比，发现FLAC3D计算结果与简单模型下精确的解析解吻合较好，但在变性较大时逐渐出现一定偏差。 关键词: 摩尔-库仑模型，增量方程，流动法则，FLAC3D 塑性流动理论的增量方程 一般情况下，破坏准则可表示为 （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） 莫尔库伦模型（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中，发生拉伸破坏，新的应力点在上，由流动法则推出。 塑性修正 首先考虑剪切破坏的情况。式（22）进行偏微分可得： （28） 将，，代换，，，式（17）变为：