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

各位大侠，我用FLAC计算地震时边坡的稳定性，按照静力计算，可是计算结果很不符合实际，麻烦各位给我看看 G5 WVr$  
hMO=#up&  
rest 初始应力.sav ^\=`ed
N0  
ini xdisp 0 ydisp 0 zdisp 0 R&k<AZ  
ini xvel 0 yvel 0 zvel 0 \~$#1D1f  
def SSR 8OU\V5i[,q  
ait1=0.02       :4/3q|cn  
k11=0.0       ;*Et[}3  
k12=1.2     &j"?\f?  
ks=(k11+k12)/2 ea 'D td  
loop while (k12-k11)>ait1 db7B^|Di  
  E1=1.573e9 ^}o2  
  poi1=0.21 Bj;'qB
>3  
  coh1=18000/ks   ",; H`V  
  fri1=(atan((tan(25*pi/180))/ks))*180/pi {4Cmu;u  
dila1=0.0 ##>H&,Dp[  
  ten1=1e6 FvjPdN/L?R  
  grav0=-10 qo bc<-  
  dens1=2000 dR,fXQm  
  K1=E1/(3*(1-2*poi1)) *.t7G  
  G1=E1/(2*(1+poi1)) l'_r:b  
  command Zb>?8  
set grav 1.5 0 grav0 $%#!bV  
model mohr <\^8fn  
pro bulk 9.1e8 she 6.5e8 co 1.90e4 fric 26 & (uE!+2C  
range group soil }Zn}  
ini dens 2000 range group soil V**~m9f  
pro bulk K1 she G1 dens dens1 coh coh1 & ]{@-HTt  
friction fri1 dil dila1 tens ten1 range group stone _]H&,</  
pro bulk 32.5e8 she 35.1e8 co 1.19e6 fric 45 & uy$e?{Jf  
range group rock p_%Rt"!  
ini dens 2650 range group rock     sUQ@7sTj  
set mech ratio 9.8e-6
@#l= l  
    solve step 10000 !_)[/q"  
  endcommand hHnYtq  
  if mech_ratio<1.0e-5 YNF k  
    k11=ks }19\.z&J  
    k12=k12 <PH#[dH  
  else \_f(M|  
    k12=ks 5U$0F$BBp  
    k11=k11 on`3&0,.  
  endif '\iCP1>+S  
    ks=(k11+k12)/2 NZ0;5xGR  
endloop HIZe0%WPw  
fosfile0='_fos0.15g'+'.sav' ;HO=  
command 2^nxoye  
  save fosfile0 H*CW1([  
endcommand 9rf)gU3{+L  
end lRFYx?y  
SSR 8<Av@9 *}  
pr ks