求助 [复制链接]

new 2c@4<kyfP  
gen zone brick p0 0 0 0 p1 60 0 0 p2 0 60 0 p3 0 0 90 & ?&EPZqI  
p4 60 60 0 p5 0 60 90 p6 60 0 150 p7 60 60 150 & |VbF&*v`  
size 6 6 10 2#~5[PtP^  
model elas U4qp?g+:  
pro bulk 10e10 she 10e10 N(q%|h<Z/=  
ini den 2500 x!fvSoHp  
apply sxx -1e9 grad 0 0 1.1111111e7 range x -.1 .1 :$."x '  
apply sxx -1e9 grad 0 0 6.6666666e6 range x 59.9 60.1 EzqYHY+_r  
apply syy -1e9 grad 0 0 8.3333333e6 range y -.1 .1 %>`0hk88  
apply syy -1e9 grad 0 0 8.3333333e6 range y 59.9 60.1 0M(\xO  
apply szz -1e8 grad 0 0 8.3333333e5 ran z 0 120 H8.U#%  
set grav 0 0 -10 L3iYZ>]  
step 30000 mG@xehH  
ini xdisp 0 ydisp 0 zdisp 0 !q=ej^(S  
ini xvel 0 yvel 0 zvel 0 ~F~g$E2 }  
plo cont szz O&!>C7  
grad部分1.1111111e7 、6.6666666e6 等是怎么得出来的，而且在不同的边界梯度值怎么还不一样？ (.4mX t  
请高手指点！ D@*<p h=  

1.1111111e7  = sxx ( -1e9 ) /  z方向的高度 ( 因为 1.1111111e7  在 z 方向位置上) }D/0&<1  
pt<!b0G  
.P+om<~B  
至于不同边界应力梯度不一样,只能说明原位地应力梯度本来就不一样