为什么计算出来的初始地应力地表应力很大，达27KPa，不是应该是0吗？求大神解答。。
命令如下：
;边界条件
fix z range z -62.9 -63.1
fix x range x -99.9 -100.1
fix x range x 99.9 100.1
fix y range y -0.1 0.1
fix y range y 1.9 2.1
;本构及参数
model mohr
def modulus
local young=1.3e9
local v=0.3125
local f1=30
local c1=3.25e5
global jianqie=young/(2*(1+v))
global tiji = young/(3*(1-2*v))
global f=f1
global c=c1
global fz=1.2*f1
global cz=1.2*c1
end
@modulus
set grav 0 0 -9.81
pro bulk @tiji sh @jianqie fric @f coh @c dens 2000
set mech ratio 1e-5
;开始计算
so

建议你到技术邻问一下，有专门的专家回答相关问题的

楼主你的问题解决了吗

楼主，您解决了吗？我也遇到相同的问题了，可以解释下吗？拜托拜托

同问

你先 用弹性模型算一遍，再把模型用MC模型算一遍试试

应该是赋予密度参数不能用pro命令，要用ini命令，因为pro命令岩石密度不变，发生变形后质量会增加