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

各位大侠，我用FLAC计算地震时边坡的稳定性，按照静力计算，可是计算结果很不符合实际，麻烦各位给我看看 Ryl:a\  
g!FuY/%+
 
rest 初始应力.sav $By<$  
ini xdisp 0 ydisp 0 zdisp 0 0%32=k7O[  
ini xvel 0 yvel 0 zvel 0 IY_iB*T3jt  
def SSR h8x MI  
ait1=0.02       OsTc5K.U~  
k11=0.0       kVuUjP6(c  
k12=1.2     ,cXD.y  
ks=(k11+k12)/2 5 )A1\  
loop while (k12-k11)>ait1 ,L|%"K]yM  
  E1=1.573e9 Ja|5 @  
  poi1=0.21 HSw;^E)1  
  coh1=18000/ks   5^}\4.eXo  
  fri1=(atan((tan(25*pi/180))/ks))*180/pi I )~GZ  
dila1=0.0 w0[6t#$F  
  ten1=1e6 N,<uf@LQ  
  grav0=-10 ({ +!`}GY  
  dens1=2000 `:ArT}F  
  K1=E1/(3*(1-2*poi1)) kS%Ydy#:'  
  G1=E1/(2*(1+poi1)) +xRK5+}9  
  command +UCG0D  
set grav 1.5 0 grav0 Hf%@3X  
model mohr u^^vB\"^  
pro bulk 9.1e8 she 6.5e8 co 1.90e4 fric 26 & p99]  
range group soil pD8+ 4;A  
ini dens 2000 range group soil bYcV$KJk  
pro bulk K1 she G1 dens dens1 coh coh1 & H/"-Z;0{  
friction fri1 dil dila1 tens ten1 range group stone t DO=P c  
pro bulk 32.5e8 she 35.1e8 co 1.19e6 fric 45 & G?9"Y%  
range group rock &r[`>B{tP  
ini dens 2650 range group rock     99]R$eT8  
set mech ratio 9.8e-6 gK+4C  
    solve step 10000 d}OTO10  
  endcommand Lt2u,9  
  if mech_ratio<1.0e-5 d|jNf</`  
    k11=ks ;RH;OE,A  
    k12=k12 9-sw!tKx  
  else M5i%jZk  
    k12=ks XP`Nf)3{Yd  
    k11=k11 fPU`/6  
  endif af6M,{F  
    ks=(k11+k12)/2 #9=Vg  
endloop pXtl 6K%  
fosfile0='_fos0.15g'+'.sav' ['b}QW@Fx  
command {WJm  
  save fosfile0 l\f*
d6o  
endcommand 3t.l5m Rg5  
end m\;R2"H%  
SSR bes<qy  
pr ks