MC模型中弹性模量的设置

总所周知，弹塑性模型包括弹性部分和塑性部分。ABAQUS中这两部分是分开设置的。具体为在Property模块中单击，如图1设置弹性部分，如图2选择塑性部分。

图1

图2

ABAQUS中内置了很多弹性模型，包括线弹性、孔隙弹性等。孔隙弹性主要配合修正剑桥模型使用。简单模型中常用的是Elastic下的Isotropic各向同性，此时只需给出弹性模量和泊松比即可。（图3）

图3

由于MC模型相对简单，参数少，很多学生喜欢用它，但是在设置模量时容易犯以下几个错误：

（1）直接将弹性模量E取为岩土勘探报告中的压缩模量，Es。这个错误是很明显的，压缩模量是侧向变形收到限制的条件下，竖向应力与应变的比值。如果土体是弹性体，压缩模量与弹性模量之间有相应的转化公式

，由于泊松比小于0.5，所以E总是小于Es的。

（2）按理论公式将压缩模量转换后的值设为弹性模量。这里存在两个问题：一是压缩模量是根据土体单元体的压缩试验得到，通常取e-p曲线上两点连线确定，其与压力范围有关系，而且加载过程实际上是个弹塑性过程，根据Es转化后的模量也只能称为变形模量，与弹塑性理论中的弹性模量还是有差异的。二是大量现场测试到的变形模量与压缩模量之间的关系与理论公式大都相反，对多数土来说，变形模量反而大于压缩模量。高大钊老师在他的《实用土力学》页上有比较详细的说明，他将其归结于土的结构性，钻土取样扰动了土的结构，室内压缩试验结果不能反映原状土的压缩特性。

（3）设置模量时不考虑实际排水条件。照理来说，土的力学行为由有效应力控制，分析时应进行包含孔压单元的有效应力分析，并设置实际的排水条件和渗透参数。但方便起见，在两种情况下可进行相应简化，即完全排水条件和完全不排水条件。这两种情况都不变采用含孔压单元，但是后者的强度参数取总应力参数，泊松比取接近0.5（模拟无体积变形），弹性模量自然要取不排水条件下的数值。对于黏土，国内外有一些根据不排水强度估算模量的经验公式。

实际上，同学们之所以选用MC模型，主要是为了分析极限承载力（强度问题），不会有人用MC模型去算固结沉降的（按照MC模型，等压加载不会引起屈服，与土体实际性状不一致），那我们在选弹性模量时自然就要模拟土体在偏压力作用之下的特点。如图4所示，如不考虑MC模型的硬化，其偏应力应变关系为一折线，实际试验做出来是个曲线。因此，我们所选择的模量要尽可能使两条曲线接近，一般的做法是选E50，即1/2破坏应力对应的割线模量。如果要考虑弹性模量随围压的变化，可以将土层划分为多个区域，每个土层采用不同的值来近似模拟。

图4