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

各位大侠，我用FLAC计算地震时边坡的稳定性，按照静力计算，可是计算结果很不符合实际，麻烦各位给我看看 ;2Q~0a|  
h;3cd0  
rest 初始应力.sav O#Wh TDF"  
ini xdisp 0 ydisp 0 zdisp 0 "fq{Y~F%`  
ini xvel 0 yvel 0 zvel 0 i*CZV|t US  
def SSR C!7>1I~5  
ait1=0.02       -(K9s!C!.  
k11=0.0       )Wb0u0)_  
k12=1.2     mS0udHod  
ks=(k11+k12)/2 5Enotp[  
loop while (k12-k11)>ait1 }`+B=h-dW  
  E1=1.573e9 Pb*5eXk  
  poi1=0.21 r_;9'#&'  
  coh1=18000/ks   GKcv<G208  
  fri1=(atan((tan(25*pi/180))/ks))*180/pi /rSH"$  
dila1=0.0 U4G`ZKv(!  
  ten1=1e6 Hty0qr3  
  grav0=-10 qY[xpm  
  dens1=2000 A/`%/0e  
  K1=E1/(3*(1-2*poi1)) tnLAJ+-M  
  G1=E1/(2*(1+poi1)) %\i9p]=  
  command F`9]=T0  
set grav 1.5 0 grav0 IUtx!.]4  
model mohr &IlU|4`R%  
pro bulk 9.1e8 she 6.5e8 co 1.90e4 fric 26 & "--t e
 
range group soil `Qeg   
ini dens 2000 range group soil ;NRF=d>  
pro bulk K1 she G1 dens dens1 coh coh1 & 0?>dCu\  
friction fri1 dil dila1 tens ten1 range group stone vB+ '  
pro bulk 32.5e8 she 35.1e8 co 1.19e6 fric 45 & c&L"N!4z  
range group rock Zdn~`Q{  
ini dens 2650 range group rock     sN5B7)Vc  
set mech ratio 9.8e-6 "1,pHR-+R  
    solve step 10000 CW<N: F.9  
  endcommand Y
tO|D  
  if mech_ratio<1.0e-5 wb~@7,D  
    k11=ks H*9~yT'Q  
    k12=k12 %w7]@VZ  
  else ]5%/3P,/  
    k12=ks M`S0u~#tI  
    k11=k11 `mN4_\]  
  endif %Z*sU/^  
    ks=(k11+k12)/2 \rPbK+G.  
endloop 8zMu7,E  
fosfile0='_fos0.15g'+'.sav' 6d+p7x  
command IT$25ZF  
  save fosfile0 Afk$?wkL  
endcommand #h/Mbj~S  
end yV^s,P1  
SSR )XWP\ h  
pr ks