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

各位大侠，我用FLAC计算地震时边坡的稳定性，按照静力计算，可是计算结果很不符合实际，麻烦各位给我看看 RC7]'4o  
! |<Fo'U  
rest 初始应力.sav svcK?^ HTe  
ini xdisp 0 ydisp 0 zdisp 0 61k"p2?+  
ini xvel 0 yvel 0 zvel 0 5YeM%%-S  
def SSR }HFN3cq;C  
ait1=0.02       Bg x'9p/  
k11=0.0       R6;=n"Ueb  
k12=1.2     \J
e0CD=e`  
ks=(k11+k12)/2 >4TaP*_  
loop while (k12-k11)>ait1 "Q@ronP(~  
  E1=1.573e9 r\'A i6  
  poi1=0.21 -g*4(w  
  coh1=18000/ks   ) l:[^$=,  
  fri1=(atan((tan(25*pi/180))/ks))*180/pi 1mOh{:1u  
dila1=0.0 iJ1"at  
  ten1=1e6 dMv=gdY  
  grav0=-10 Z`ID+  
  dens1=2000 :VRNs  
  K1=E1/(3*(1-2*poi1)) z[L8
$7L  
  G1=E1/(2*(1+poi1)) o,AAC
  
  command !Prg_6 `  
set grav 1.5 0 grav0 mwHB(7YS,  
model mohr 7q?, ?  
pro bulk 9.1e8 she 6.5e8 co 1.90e4 fric 26 & oAODp!_c  
range group soil $P/~rZ@M@  
ini dens 2000 range group soil OEA&~4&{7  
pro bulk K1 she G1 dens dens1 coh coh1 & ^ *k?pJ5  
friction fri1 dil dila1 tens ten1 range group stone 'vbsvT  
pro bulk 32.5e8 she 35.1e8 co 1.19e6 fric 45 & jFL #s&ft  
range group rock cPyE 6\lN  
ini dens 2650 range group rock     P}n_IV*@  
set mech ratio 9.8e-6 X86O lP)eX  
    solve step 10000 a|t$l=|DD  
  endcommand Jh,]r?Bd  
  if mech_ratio<1.0e-5 XDOY`N^L  
    k11=ks sBvzAVBL  
    k12=k12 9
6( v  
  else ;-~B)M_S`  
    k12=ks `yrB->|vG  
    k11=k11 )rxX+k+b/  
  endif xr4*{v  
    ks=(k11+k12)/2 I9_RlAd  
endloop 8lvV4
yb  
fosfile0='_fos0.15g'+'.sav' ;g+N&)n  
command g+vv
a"  
  save fosfile0 /Kql>$I  
endcommand RO+GK`J  
end gY/"cq  
SSR cQCSe,$ W  
pr ks