非均质材料中的初始渐变应力计算问题，每种材料分界面上的应力须手工计算，本人不理解怎么计算的，求助大家。
例子为彭文斌的Flac3D实用教程，P139.
new
gen zone brick 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 gravity 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 -10,0
ini syy 4e4 grad 0,2.0e4,0 range y -15,-10
ini syy 7e4 grad 0,2.2e4,0 range y -15,-25
step 1
请教：syy 0 4e4 和7e4是怎么计算出来的，一直困扰我，谢谢！！
