求助 [复制链接]

new uXFI7vV6P  
gen zone brick p0 0 0 0 p1 60 0 0 p2 0 60 0 p3 0 0 90 & m0ra  
p4 60 60 0 p5 0 60 90 p6 60 0 150 p7 60 60 150 & ;y\/7E  
size 6 6 10 &nr{-][  
model elas |=YK2};  
pro bulk 10e10 she 10e10 U&])ow):  
ini den 2500 ,P}7e)3  
apply sxx -1e9 grad 0 0 1.1111111e7 range x -.1 .1 ohKoX$|p~  
apply sxx -1e9 grad 0 0 6.6666666e6 range x 59.9 60.1 D
s"
%=  
apply syy -1e9 grad 0 0 8.3333333e6 range y -.1 .1 B2]52Fg-"  
apply syy -1e9 grad 0 0 8.3333333e6 range y 59.9 60.1 hO$Gx*e$  
apply szz -1e8 grad 0 0 8.3333333e5 ran z 0 120 &v((tZ  
set grav 0 0 -10 +`Q]p"G  
step 30000 [q!]Ds" _  
ini xdisp 0 ydisp 0 zdisp 0 k-n`R)p:  
ini xvel 0 yvel 0 zvel 0 >Z_;ZMu)  
plo cont szz $%bd`d*S  
grad部分1.1111111e7 、6.6666666e6 等是怎么得出来的，而且在不同的边界梯度值怎么还不一样？ F*J1w|)F0  
请高手指点！ DVhBZ!u9  

1.1111111e7  = sxx ( -1e9 ) /  z方向的高度 ( 因为 1.1111111e7  在 z 方向位置上) )C?H m^#  
eC<RM Q4  
Ey:?!  
至于不同边界应力梯度不一样,只能说明原位地应力梯度本来就不一样