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

各位大侠，我用FLAC计算地震时边坡的稳定性，按照静力计算，可是计算结果很不符合实际，麻烦各位给我看看 O!>#q4&]  
mfj%-)l9  
rest 初始应力.sav 7/M[T\c  
ini xdisp 0 ydisp 0 zdisp 0 @ddCVxd  
ini xvel 0 yvel 0 zvel 0 AxEdQRGk  
def SSR K91O$'J  
ait1=0.02       ?h1g$SBxk  
k11=0.0       w3i74C&0  
k12=1.2     Q|[^dju  
ks=(k11+k12)/2 U;kNo3=  
loop while (k12-k11)>ait1 u~,hTY(%  
  E1=1.573e9 2
v(Y'f.  
  poi1=0.21 4,@jSr|I3i  
  coh1=18000/ks   1OvoW Nx  
  fri1=(atan((tan(25*pi/180))/ks))*180/pi gs+nJ+b  
dila1=0.0 F,as>X#  
  ten1=1e6 BNk>D|D;  
  grav0=-10 S*n5d>;  
  dens1=2000 k@9hth2Q  
  K1=E1/(3*(1-2*poi1)) ){mqo
%{SO  
  G1=E1/(2*(1+poi1)) :CqR1_n%  
  command {li Q&AZ  
set grav 1.5 0 grav0 x _kT Wq  
model mohr N MkOx$  
pro bulk 9.1e8 she 6.5e8 co 1.90e4 fric 26 & uo?R;fX26  
range group soil pFs/ipZX^*  
ini dens 2000 range group soil ,2
U  
pro bulk K1 she G1 dens dens1 coh coh1 & (mbm',%-(  
friction fri1 dil dila1 tens ten1 range group stone k]AL\) &W  
pro bulk 32.5e8 she 35.1e8 co 1.19e6 fric 45 & +9]t]Vrw  
range group rock },X.a@:  
ini dens 2650 range group rock     V(;T{HW&  
set mech ratio 9.8e-6 `_.:O,^n^  
    solve step 10000 3rMi:*?  
  endcommand kbvF 9#  
  if mech_ratio<1.0e-5 QeT~s5 H  
    k11=ks Mk9J~'C_  
    k12=k12 cjtcEW  
  else ]w,
|WZm  
    k12=ks HoX={^aG%  
    k11=k11 ` >>]$ZJ  
  endif cxQ %tL+S&  
    ks=(k11+k12)/2 cSv;HN:  
endloop hVyeHbx  
fosfile0='_fos0.15g'+'.sav' Zqf ovG
 
command A|PZ<WAY  
  save fosfile0 3-AOB3](
 
endcommand H`k YDp  
end uzb|yV'B  
SSR \PZ;y=]p}  
pr ks