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

各位大侠，我用FLAC计算地震时边坡的稳定性，按照静力计算，可是计算结果很不符合实际，麻烦各位给我看看 i7-i!`<  
.`~?w+ ~  
rest 初始应力.sav bvJ@H Z$  
ini xdisp 0 ydisp 0 zdisp 0 =v^#MU{k?  
ini xvel 0 yvel 0 zvel 0 XYR q"{Id  
def SSR C-S>'\|8  
ait1=0.02       QxG^oxU}  
k11=0.0       xTuJ~$(  
k12=1.2     |pS]zD  
ks=(k11+k12)/2 m-$}'mEO  
loop while (k12-k11)>ait1 0?oL zw&  
  E1=1.573e9 Ep
O2%|@  
  poi1=0.21 p*5_+u  
  coh1=18000/ks  
S'^ q  
  fri1=(atan((tan(25*pi/180))/ks))*180/pi 1K#[Ef4  
dila1=0.0 ;o'r@4^&$R  
  ten1=1e6 lTu&
9)  
  grav0=-10 N$Ad9W?T  
  dens1=2000 
?\8  
  K1=E1/(3*(1-2*poi1)) 5.ab/uk;M  
  G1=E1/(2*(1+poi1)) =PY{Elf  
  command QY4;qA  
set grav 1.5 0 grav0 T16gq-h'  
model mohr UUDHknm"  
pro bulk 9.1e8 she 6.5e8 co 1.90e4 fric 26 & 7b08Lo7b  
range group soil kh#QT_y  
ini dens 2000 range group soil ZHjL8Iq  
pro bulk K1 she G1 dens dens1 coh coh1 & UapU:>!"`  
friction fri1 dil dila1 tens ten1 range group stone ,9d]-CuP;  
pro bulk 32.5e8 she 35.1e8 co 1.19e6 fric 45 & VqvjOeCbH  
range group rock 5?HwM[`  
ini dens 2650 range group rock     .'A1Eoo0d  
set mech ratio 9.8e-6 N@tKgx  
    solve step 10000 J*b Je"8  
  endcommand %\,9S`0  
  if mech_ratio<1.0e-5 '*L6@e#U  
    k11=ks _BA; H+M  
    k12=k12 n%29WF6Zf  
  else LI@BB:)[  
    k12=ks )V~=B]  
    k11=k11 EPn!6W5^  
  endif SF-E>s!XL  
    ks=(k11+k12)/2 ;<m
*ASM.3  
endloop T8Mqu`$r  
fosfile0='_fos0.15g'+'.sav' i$%Bo/Y   
command c*7|>7C$i  
  save fosfile0 lCUYE"o  
endcommand G=[<KtWa  
end !AJkd.  
SSR \,Ws=9f  
pr ks