巷道开挖，塑性区求解释

cumtmbqq1

此矩形巷道宽8m，高3.2m，简化为平面应变模型，用摩尔库伦模型计算，
此分析遇到如下几个问题：
（1）若不考虑剪涨，即把剪胀角设为0，则计算不收敛
（2）考虑剪涨，打开大变形即几何非线性，计算时间步算到0.978时即不收敛
（3）摩尔库伦破坏准则考虑的只有剪切破坏，然而在巷道顶板处显然为拉伸破坏，用摩尔库伦模型不能体现出来，这样算出来的结果与现实差距较大，因为顶板处肯定有塑性区，此问题解决方法应该是要编写摩尔库伦模型与最大拉应力准则即第一强度理论共同作用的umat，不知哪位大侠手中有子程序。
材料参数如下：
*Material, name=Material-1
*Density
2000.,
*Elastic
4e+09, 0.25
*Mohr Coulomb
30.,30.
*Mohr Coulomb Hardening
2e+06,0.
模型

网格

边界条件

塑性区

最大拉应力准则与摩尔库伦准则umat中的公式推导参考文献：

bigcrown

怎么实现我不懂，不过你要做的东西和我要做的，有相同之处，我也想把这个拉伸破裂区搞出来，我看过国外的文章，他们是计算一个最小主应力还是平均主应力，当其值小于-c cot摩擦角，时，就默认发生了拉伸破裂，如果你能把这个力表示出来，这个拉伸破裂就出来了。

cumtmbqq1

2# bigcrown
看来做巷道用abaqus还真是不太好了。还是用flac吧。