非均质材料的初始渐变应力如何定义

cumt-李文峰

彭文斌版flac3D上面有一个例子：非均质材料给出初始应力，0-10m深度介质密度1600,10-15m深度介质密度2000,15-25m深度介质密度2200，下面是给出模型的初始渐变应力的命令：

new
gen zo b size 10 10 10 p0 0 -25 0 p1 20 -25 0 p2 0 0 0 p3 0 -25 20

model elas
prop bulk 5e9 shear 3e9
ini density 1600 range y -10 0
ini density 2000 range y -15 -10
ini density 2200 range y -25 -15

set grav 0 -10 0
fix x range x -.1 .1
fix x range x 19.9 20.1
fix z range z -.1 .1
fix z range z 19.9 20.1
fix y range y -25.1 -24.9

ini syy 0.0   grad 0 1.6e4 0 range y 0 -10
ini syy 4e4 grad 0 2.0e4 0 range y  -10 -15
ini syy 7e4 grad 0 2.2e4 0 range y -15 -25

step 1
pl contour syy



0-10m范围：syy应力初始值为0，应力梯度为1.6e4，
可为什么10--15m范围：应力初始值为4e4？该syy应力初始值加载在哪个边界面？
15-25m范围应力初始条件存在着同样的疑问.......
麻烦各位点拨.......