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

各位大侠，我用FLAC计算地震时边坡的稳定性，按照静力计算，可是计算结果很不符合实际，麻烦各位给我看看 9ec0^T  
\!m!ibr  
rest 初始应力.sav BZ2nDW*%  
ini xdisp 0 ydisp 0 zdisp 0 (}sDm~;s  
ini xvel 0 yvel 0 zvel 0 u1wg C#  
def SSR hzW{_Q.|?  
ait1=0.02       S\}?zlV  
k11=0.0       9N<*S'Z  
k12=1.2     t2Ip\>;9f  
ks=(k11+k12)/2 p|b+I"M  
loop while (k12-k11)>ait1 vP-3j  
  E1=1.573e9 F
ZM2   
  poi1=0.21 SHytyd  
  coh1=18000/ks   jGJ.P
vc>i  
  fri1=(atan((tan(25*pi/180))/ks))*180/pi B<(Pd  
dila1=0.0 N\c&PS  
  ten1=1e6 uAc@ Z-  
  grav0=-10 :EYu 4Y  
  dens1=2000 sbs[=LW4  
  K1=E1/(3*(1-2*poi1)) poTl|y @  
  G1=E1/(2*(1+poi1)) RfD{g"]y  
  command ]Jja  
set grav 1.5 0 grav0 0`V3s]%iu  
model mohr OiYNH~hv  
pro bulk 9.1e8 she 6.5e8 co 1.90e4 fric 26 & z|Hy>|+  
range group soil nMTLD  
ini dens 2000 range group soil '" ^ B&W  
pro bulk K1 she G1 dens dens1 coh coh1 & #4Dn@Gqh.Y  
friction fri1 dil dila1 tens ten1 range group stone ]jZiW1C*a  
pro bulk 32.5e8 she 35.1e8 co 1.19e6 fric 45 & KuIBYaK, g  
range group rock 71JM [2  
ini dens 2650 range group rock     lb-S0plw  
set mech ratio 9.8e-6 ,Le&I9*%  
    solve step 10000 R5m`;hF  
  endcommand ;r gH}r  
  if mech_ratio<1.0e-5 N>Vacc_[  
    k11=ks eo>/  
    k12=k12 ^fFtI?.6jI  
  else :D<:N*9i  
    k12=ks x:!C(Ep)  
    k11=k11 (~q#\  
  endif T@%;0Ro~  
    ks=(k11+k12)/2 k&MlQ2'!<  
endloop aQl?d<|+lk  
fosfile0='_fos0.15g'+'.sav' D?iy.Dg  
command 6U;Jg
_zS  
  save fosfile0 Cf.pTYSl  
endcommand "Cz
z,;0  
end pe\Txg
6  
SSR  R4&|t  
pr ks