模型Z（+90m~-580m）+90m是地表，一开始我加的X,Y,Z三方向的应力，命令如下：
impgrid f:\Flacmodel_haitang.Flac3D
model mo
ini dens 2700 ran group 1
prop bulk 3.33e9 shear 2.5e9 fric 50 coh 20e6 ten 21e6 ran group 1
ini dens 2700 ran group 2
prop bulk 3.33e9 shear 2.5e9 fric 50 coh 20e6 ten 21e6 ran group 2
ini dens 2689 ran group 3 z -580 90
prop bulk 4.34e9 shear 2.37e9 fric 43.9 coh 15.1e6 ten 15.56e6 ran group 3
attach face all
fix x y z range z -580.1 -579.9
fix x y range x -450.1 -449.9
fix x y range x 2699.9 2700.1
fix x y range y -250.1 -249.9
fix x y range y 199.9 200.1
ini szz 0 grad 0,0,4.3e4
ini syy 0 grad 0,0,1.433e4
ini sxx 0 grad 0,0,1.433e4
set gravity 0,0,-9.81
solve
save e:\初始0.sav
X,Y方向的应力图不尽如人意：
后来只施加自重应力 还是不科学：
impgrid D:\Flacmodel_haitang.Flac3D
model elas
ini dens 2700 ran group 1
prop bulk 10e9 shear 7e9 fric 50 coh 20e6 ten 21e6 ran group 1
ini dens 2700 ran group 2
prop bulk 10e9 shear 7e9 fric 50 coh 20e6 ten 21e6 ran group 2
ini dens 2689 ran group 3
prop bulk 10e9 shear 7e9 fric 45 coh 15e6 ten 15e6 ran group 3
attach face all
apply nstress -6.6e6 range z -59.9 -60.1
fix x y z range z -580.1 -579.9
fix x range x -450.1 -449.9
fix x range x 2699.9 2700.1
fix y range y -250.1 -249.9
fix y range y 199.9 200.1
set gravity 0,0,-9.81
solve
save d:\初始elas.sav
模型X(-450,2700)Y(-250,200)Z(-580,90)网格尺寸是30*30*30的
请教各位大神 是不是我命令写错了 自重是按书上的步骤写的 ，构造的
σZ=容重*H，σx=σy=1/3σZ。请问哪里错了 ？谢谢各位了 
