求助 [复制链接]

new 9A9T'g)Du  
gen zone brick p0 0 0 0 p1 60 0 0 p2 0 60 0 p3 0 0 90 & 2{B ScI5K  
p4 60 60 0 p5 0 60 90 p6 60 0 150 p7 60 60 150 & *dsI>4%m  
size 6 6 10 F]s:`4  
model elas #IU^(W  
pro bulk 10e10 she 10e10 4AKPS&k;  
ini den 2500 uhB V)Qg  
apply sxx -1e9 grad 0 0 1.1111111e7 range x -.1 .1 kX+98?h-C  
apply sxx -1e9 grad 0 0 6.6666666e6 range x 59.9 60.1 as[! 9tB]
 
apply syy -1e9 grad 0 0 8.3333333e6 range y -.1 .1 d
*#.(C9^  
apply syy -1e9 grad 0 0 8.3333333e6 range y 59.9 60.1 /+m2|Ij(  
apply szz -1e8 grad 0 0 8.3333333e5 ran z 0 120 #WAX&<m  
set grav 0 0 -10 Sar1NkD#  
step 30000 j@{dsS:6  
ini xdisp 0 ydisp 0 zdisp 0 er3`ITp:dp  
ini xvel 0 yvel 0 zvel 0 6.6?Rp".  
plo cont szz 2)-4?uz~  
grad部分1.1111111e7 、6.6666666e6 等是怎么得出来的，而且在不同的边界梯度值怎么还不一样？ Q&/WVRD  
请高手指点！ m
!]J{OGG:  

1.1111111e7  = sxx ( -1e9 ) /  z方向的高度 ( 因为 1.1111111e7  在 z 方向位置上) A;o({9VH`Z  
k"%sdYkb!  
k;)mc+ ~+  
至于不同边界应力梯度不一样,只能说明原位地应力梯度本来就不一样