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

各位大侠，我用FLAC计算地震时边坡的稳定性，按照静力计算，可是计算结果很不符合实际，麻烦各位给我看看 Im' :sJ31  
+s*l#'Q  
rest 初始应力.sav _($-dJ{  
ini xdisp 0 ydisp 0 zdisp 0 CL<KBmW7  
ini xvel 0 yvel 0 zvel 0 >2/wzs
W  
def SSR 1K*f4BnDr~  
ait1=0.02       DO1N`7@o  
k11=0.0       "M5ro$qZ}  
k12=1.2     Ls$g-k%c@Q  
ks=(k11+k12)/2 6ljRV)  
loop while (k12-k11)>ait1 P98X[0&  
  E1=1.573e9 \0^rJ1*  
  poi1=0.21 M|Lw`?T  
  coh1=18000/ks   G\,A> mT/P  
  fri1=(atan((tan(25*pi/180))/ks))*180/pi E^EU+})Ujr  
dila1=0.0 TW-zh~|F  
  ten1=1e6 iEJQ#5))0  
  grav0=-10 \{+nXn  
  dens1=2000 iuU3*yyn  
  K1=E1/(3*(1-2*poi1)) kb}]sj  
  G1=E1/(2*(1+poi1)) dPyBY]`  
  command dx?njR  
set grav 1.5 0 grav0 >m_p\$_  
model mohr n[MIa]dK  
pro bulk 9.1e8 she 6.5e8 co 1.90e4 fric 26 & b0z{"  
range group soil Ae{4AZ  
ini dens 2000 range group soil tE9_dR^K  
pro bulk K1 she G1 dens dens1 coh coh1 & x9o^9QJh  
friction fri1 dil dila1 tens ten1 range group stone @L>NN>?SGQ  
pro bulk 32.5e8 she 35.1e8 co 1.19e6 fric 45 & Vf,~MG  
range group rock Edn$0D68u_  
ini dens 2650 range group rock     G(|(y=ck  
set mech ratio 9.8e-6 G4iLCcjY  
    solve step 10000 !ovZ>,1  
  endcommand 1A/c/iC  
  if mech_ratio<1.0e-5 2t`9_z
qLw  
    k11=ks meM.?kk(  
    k12=k12 +v[$lh+  
  else M5ySs\O4  
    k12=ks M\.T 0M_  
    k11=k11 HBga'x
J  
  endif ,d [b"]Zy  
    ks=(k11+k12)/2 :86luLFm  
endloop 7p>-oR"  
fosfile0='_fos0.15g'+'.sav' A*26'  
command GZhfA ;O,  
  save fosfile0 3EW f|6RI  
endcommand XpAq=p0;  
end "TH-A6v1  
SSR @|UIV  
pr ks