为什么计算出来的初始地应力地表应力很大，达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
命令如下：
;边界条件
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
