求助 [复制链接]

new ^y KkWB*  
gen zone brick p0 0 0 0 p1 60 0 0 p2 0 60 0 p3 0 0 90 & q\|RI;W  
p4 60 60 0 p5 0 60 90 p6 60 0 150 p7 60 60 150 & F",TP,X  
size 6 6 10 a]xGzv5  
model elas *Uj;a.  
pro bulk 10e10 she 10e10 `b]wyP  
ini den 2500 LME&qKe5  
apply sxx -1e9 grad 0 0 1.1111111e7 range x -.1 .1 h]+;"v6 /  
apply sxx -1e9 grad 0 0 6.6666666e6 range x 59.9 60.1 6H5o/)Q~  
apply syy -1e9 grad 0 0 8.3333333e6 range y -.1 .1 (Y8LyY  
apply syy -1e9 grad 0 0 8.3333333e6 range y 59.9 60.1 G2?#MO  
apply szz -1e8 grad 0 0 8.3333333e5 ran z 0 120 6vR6=@(`>  
set grav 0 0 -10 
C"B'Dj  
step 30000 
f1+  
ini xdisp 0 ydisp 0 zdisp 0 sQ,xTWdj  
ini xvel 0 yvel 0 zvel 0 hITYBPqRO  
plo cont szz -2>s#/%  
grad部分1.1111111e7 、6.6666666e6 等是怎么得出来的，而且在不同的边界梯度值怎么还不一样？ 8iOHav4  
请高手指点！ EP>Lh7E9n  

1.1111111e7  = sxx ( -1e9 ) /  z方向的高度 ( 因为 1.1111111e7  在 z 方向位置上) h6^|f%\w*i  
a:`<=^:4,  
a0gg<Ml  
至于不同边界应力梯度不一样,只能说明原位地应力梯度本来就不一样