[求助]FLAC地震计算 [复制链接]

各位大侠，我用FLAC计算地震时边坡的稳定性，按照静力计算，可是计算结果很不符合实际，麻烦各位给我看看 wQP^WzNE  
%k-3?%&8  
rest 初始应力.sav ?Z Rs\+{vG  
ini xdisp 0 ydisp 0 zdisp 0 Rot@x r7Hc  
ini xvel 0 yvel 0 zvel 0 s2'yY(u/  
def SSR ewa wL"  
ait1=0.02       8+a
4>8[M  
k11=0.0       M)'HCnvs'  
k12=1.2     0B@SN)<kH  
ks=(k11+k12)/2 kzn5M&f>  
loop while (k12-k11)>ait1 +BVym~*^  
  E1=1.573e9 +$C4\$t  
  poi1=0.21 k7?N ?7w  
  coh1=18000/ks   [ oL.+  
  fri1=(atan((tan(25*pi/180))/ks))*180/pi @ezH'y-v  
dila1=0.0 M m[4yP%  
  ten1=1e6 4L ;% h  
  grav0=-10 !c}O5TI|#  
  dens1=2000 pm*xb
]8y  
  K1=E1/(3*(1-2*poi1)) z.$4!$q  
  G1=E1/(2*(1+poi1)) ORyE`h
 
  command P3cRl']  
set grav 1.5 0 grav0 u [m  
model mohr 6+.uU[x@  
pro bulk 9.1e8 she 6.5e8 co 1.90e4 fric 26 & ??PC k1X  
range group soil 9o;^[Ql-  
ini dens 2000 range group soil :=/85\P0SU  
pro bulk K1 she G1 dens dens1 coh coh1 & #33fGmd[  
friction fri1 dil dila1 tens ten1 range group stone ?ooe'V@  
pro bulk 32.5e8 she 35.1e8 co 1.19e6 fric 45 & yaah*1ip[  
range group rock 7ePqmB<.  
ini dens 2650 range group rock     @ zs'Y8  
set mech ratio 9.8e-6 l)Pu2!Ic  
    solve step 10000 *AoR==:ya  
  endcommand X%Z{K-  
  if mech_ratio<1.0e-5 P|.]DJ  
    k11=ks LQtj~c>X-|  
    k12=k12 !lf|7  
  else 0~ZFv Wv  
    k12=ks J?RabYd ~  
    k11=k11 OK{quM5  
  endif 3a?o3=  
    ks=(k11+k12)/2 *+j*{>E  
endloop gZLP\_CL  
fosfile0='_fos0.15g'+'.sav' "OQ^U_  
command 2_n7=&  
  save fosfile0 Fs{x(_LOr  
endcommand z.GMqW%B  
end :uWw8`  
SSR 85n1eE  
pr ks