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

各位大侠，我用FLAC计算地震时边坡的稳定性，按照静力计算，可是计算结果很不符合实际，麻烦各位给我看看 ^2|G0d@.:  
&3o[^_Ti  
rest 初始应力.sav 7jf%-X  
ini xdisp 0 ydisp 0 zdisp 0 aDae0$lc.S  
ini xvel 0 yvel 0 zvel 0 tW=,o&C=  
def SSR No'^]r  
ait1=0.02       _w'N&#  
k11=0.0       ",pN.<F9O  
k12=1.2     E&RiEhuv  
ks=(k11+k12)/2 1z#0CX}Y/H  
loop while (k12-k11)>ait1 Q)H1
\  
  E1=1.573e9 PCx] >&  
  poi1=0.21 Xyf
7sHQ  
  coh1=18000/ks   ^EG@tB $<  
  fri1=(atan((tan(25*pi/180))/ks))*180/pi .t7mTpi  
dila1=0.0 |.qK69  
  ten1=1e6 q3~RK[OCq  
  grav0=-10 0VA$ Ige  
  dens1=2000 [^J2<\<0  
  K1=E1/(3*(1-2*poi1)) _i-
(`5  
  G1=E1/(2*(1+poi1)) V&j |St[  
  command ?11\@d  
set grav 1.5 0 grav0 gk*Md+  
model mohr ckZZ)lW`*  
pro bulk 9.1e8 she 6.5e8 co 1.90e4 fric 26 & ;`TSu5/  
range group soil $;ssW"7~Qn  
ini dens 2000 range group soil )(Z)yz  
pro bulk K1 she G1 dens dens1 coh coh1 & 2u> [[U1:  
friction fri1 dil dila1 tens ten1 range group stone Rzolue 8  
pro bulk 32.5e8 she 35.1e8 co 1.19e6 fric 45 & x;G~c5  
range group rock Gz>M`M`[
4  
ini dens 2650 range group rock     oI!L2  
set mech ratio 9.8e-6 RgPY,\_9+  
    solve step 10000 6*cY[R|q!  
  endcommand AIx,c1G]K  
  if mech_ratio<1.0e-5 K?4/x4p@  
    k11=ks 'A|OVyH  
    k12=k12 X+:>&&9  
  else B;ro(R  
    k12=ks @`</Z)  
    k11=k11 ~(|~Ze>  
  endif uMB|x,X I  
    ks=(k11+k12)/2 D>YbL0K>X~  
endloop +>F #{b  
fosfile0='_fos0.15g'+'.sav' <IX)D `mf  
command K]b_JDEk  
  save fosfile0 #*g.hL<  
endcommand =t+{)d.w  
end Hko(@z  
SSR L 1fK  
pr ks