求助 [复制链接]

new )e5 @  
gen zone brick p0 0 0 0 p1 60 0 0 p2 0 60 0 p3 0 0 90 & OnQdq^UB  
p4 60 60 0 p5 0 60 90 p6 60 0 150 p7 60 60 150 & #@<9S{F  
size 6 6 10 'AN>`\mR$  
model elas 3k#~yaoI  
pro bulk 10e10 she 10e10 BauU{:Sh  
ini den 2500 n]E?3UGD@W  
apply sxx -1e9 grad 0 0 1.1111111e7 range x -.1 .1 MxT-1&XL  
apply sxx -1e9 grad 0 0 6.6666666e6 range x 59.9 60.1 p w8 s8?  
apply syy -1e9 grad 0 0 8.3333333e6 range y -.1 .1 O T.*pk+<)  
apply syy -1e9 grad 0 0 8.3333333e6 range y 59.9 60.1 y=?)n\f  
apply szz -1e8 grad 0 0 8.3333333e5 ran z 0 120 /c uLc^(X  
set grav 0 0 -10 :oJ=iB'Zc  
step 30000 94B\5I}  
ini xdisp 0 ydisp 0 zdisp 0 0lhVqy}:}o  
ini xvel 0 yvel 0 zvel 0 8Z(Mvq]f&  
plo cont szz d3=KTTi\  
grad部分1.1111111e7 、6.6666666e6 等是怎么得出来的，而且在不同的边界梯度值怎么还不一样？ O,Cb"{qH8  
请高手指点！ 9eH$XYy
 

1.1111111e7  = sxx ( -1e9 ) /  z方向的高度 ( 因为 1.1111111e7  在 z 方向位置上) JU)dr4S?  
S[J}UpV  
JsODzw  
至于不同边界应力梯度不一样,只能说明原位地应力梯度本来就不一样