教材例子
new
gen zone brick size 10 10 10 p0 0,-25,0 p1 20,-25,0 p2 0,0,0 p3 0,-25,20
plo blo gro
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
solve
非均质材料中渐变初应力，教材中说每种材料分界面上的应力必须经过手工计算，真实情况下弹性模量可能发生变化，将出现水平应
力，如果存在较高水平应力则需要用ini命令来设置。
我想知道 命令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
中syy 4e4和syy 7e4是如何计算得到的？
