请教隧道开挖模拟的初始地应力问题

dreamsoil-ep

我用下述材料来模拟隧道周围岩土性质。

*MATERIAL,NAME=ROCK1
*DRUCKER PRAGER,SHEAR CRITERION=LINEAR
38.0, 0.8, 22.
*DRUCKER PRAGER HARDENING
12.6,0.0
24.3,6.33E-5
30.6,1.9E-4
36.4,5.7E-4
37.9, 8.23E-04
36.4, 1.33E-03
31.1, 1.96E-03
23.8, 2.66E-03
13.1, 3.48E-03
7.28, 4.30E-03
4.37, 5.00E-03
3.40, 5.76E-03
2.91, 6.39E-03
1.94, 6.96E-03
1.94, 7.53E-03
1.46, 8.16E-03
1.46, 8.67E-03
1.46, 9.43E-03
由于是深埋隧道，所以，在所建模型的上端施加竖向荷载，代表其上的上覆土重
*DLOAD
SURFACE, P3, 19.6
初始地应力分布为
*INITIAL CONDITIONS,TYPE=STRESS,GEOSTATIC
EALL,-19.6,22700,-20.96465,-33000,0.33333
约束条件是隧道对称轴上的水平约束和隧道轮廓外表的水平和竖向的约束，上端面施加荷载，其他面即右端面和下端面用的是无限单元，所以不加约束。
分析步是
*STEP, UNSYMM=YES
*GEOSTATIC
可结果是不能收敛
我看abaqus的例子中用的也是一些复杂的本构模型，这样的分析应该是没问题的，但不知是为什么，还请各位指点

aba_aba

请上传模型

dreamsoil-ep

模型中网格有两个节点，在*MPC下的TIE和LINEAR中同时用到（为了refine，经验不足），出现警告。但总的计算应该没问题。因为，用弹性算没有问题。感谢aba_aba 的关注。

aba_aba

msg文件中没有警告信息，计算也收敛了，没什么问题。

dreamsoil-ep

十分谢谢！

chenfanjian

想问一下楼主,硬化曲线的数据是怎么得到的?