求助 [复制链接]

new +d2+w1o^V  
gen zone brick p0 0 0 0 p1 60 0 0 p2 0 60 0 p3 0 0 90 &  N_=7  
p4 60 60 0 p5 0 60 90 p6 60 0 150 p7 60 60 150 & ?A
8Uf=  
size 6 6 10 m Wsegq4  
model elas Z>Rd6o'  
pro bulk 10e10 she 10e10 e-]k{_wm  
ini den 2500 [64K?l0&  
apply sxx -1e9 grad 0 0 1.1111111e7 range x -.1 .1 wGBQ.Ve[  
apply sxx -1e9 grad 0 0 6.6666666e6 range x 59.9 60.1 o%+KS5v!  
apply syy -1e9 grad 0 0 8.3333333e6 range y -.1 .1 d9R0P2  
apply syy -1e9 grad 0 0 8.3333333e6 range y 59.9 60.1 Km'd=B>Jy  
apply szz -1e8 grad 0 0 8.3333333e5 ran z 0 120 O (tcu@vfl  
set grav 0 0 -10  bXQ(6P  
step 30000
4O[5,  
ini xdisp 0 ydisp 0 zdisp 0 lF8B+  
ini xvel 0 yvel 0 zvel 0 v7\~OOoH]  
plo cont szz +`yDWN?7  
grad部分1.1111111e7 、6.6666666e6 等是怎么得出来的，而且在不同的边界梯度值怎么还不一样？ <(>v|5K0]  
请高手指点！ st:[|`  

1.1111111e7  = sxx ( -1e9 ) /  z方向的高度 ( 因为 1.1111111e7  在 z 方向位置上) jO3Q@N0_  
rQ30)5^V|  
zblh_6  
至于不同边界应力梯度不一样,只能说明原位地应力梯度本来就不一样