边坡稳定性计算 [复制链接]

Flac计算安全系数，结果velocity却显示错误。这是什么原因？

图片:00.png

发上命令流或者其他的图，就这一个图，看不明白呀！

命令流： !RX\">z  
zN}1Qh  
new DmDsn  
def _parm Uedzt  
hM}rf6B  
        ;-----------------几何参数---------------------------------------------------- Yq~$Q4  
_$s ;QI]x  
        p_d1=0.15                      ;-----------桩径------------------------------- 7ux0|l  
        tslope=7.5                     ;-----------坡顶------------------------------- )1 ]P4  
        dslope=-7.5                    ;-----------坡脚------------------------------- {OFbU  
        ngrad01=1.5                    ;-----------坡率------------------------------- 4n6EkTa  
        s1=3.0*p_d1                    ;-----------桩排距S1=3D------------------------- #E35%7*  
        s2=3.0*p_d1                    ;-----------桩间距S2=3D---------------------- 6)<g%bH!
  
        s=2.0*s2                       ;-----------模型厚度S=2S2---------------------- 5G5P#<Vv  
        sl1=s-0.1 * V;L|c  
        sl2=s+0.1 zTA+s 2  
oU/CXz?H  
        ;-----------------土体材料参数---------------------------------------------------- >6fc`3*!  
         .|>zQ(7YC  
        s_Emod=1.0e7 }:JE*D|  
        s_nu=0.3 q\+khy,k  
        s_fri=20.0 A4;EtW+F  
        s_coh=10e3
CjtBQ5  
        s_tens=0 z&fXxp  
        soil_dens=2000 <1")JDW  
        soil_bulk=s_Emod/3/(1-2*s_nu) 6z;C~_BV  
        soil_shear=s_Emod/2/(1+s_nu) },r30`)Q  

<dzfD;  
5f5bhBZ<  
        ;-----------------桩位---------------------------------------------------- kU uDA><1  
         Dn{ hU$*  
        ;-----------------1# 桩----------------------     +/!kL0[v  
        x_beg01=-s1/2 )qXl8HI  
        y_beg01=0.5*s2     LxD >eA  
        z_beg01=10.0+5.0*(7.5+x_beg01)/7.5 @9-/p^n1  
wHneVqI/U  
    x_end01=-s1/2 2.''Nt6|  
        y_end01=0.5*s2 b{T". @b  
    z_end01=0 yX%q7ex
 
         b4TZnO  
        ;-----------------2# 桩---------------------- )_[eqr  
        x_beg02=-s1/2 pkWJb!  
        y_beg02=1.5*s2     [C1.*Q+l  
        z_beg02=10.0+5.0*(7.5+x_beg01)/7.5 l!r2[T]I@7  
50MdZ;R-3  
    x_end02=-s1/2 5]C}044  
        y_end02=1.5*s2 <ezv  
    z_end02=0 `re9-HM  
$|J16tW  
        ;-----------------3# 桩---------------------- *Uq1q  
        x_beg03=s1/2 {+D 6o  
        y_beg03=0.5*s2     0 #*M'C#  
        z_beg03=10.0+10.0*(7.5+x_beg03)/15 E?$|`<o{|`  
\#1*r'V8  
    x_end03=s1/2 uu08q<B5b)  
        y_end03=0.5*s2 ]/byz_7]  
    z_end03=0 TL^af
-  
%V r vu5  
        ;-----------------4# 桩---------------------- nR%ASUx:Y  
        x_beg04=s1/2 :|j,x7&/{  
        y_beg04=1.5*s2     Qsv3`c  
        z_beg04=10.0+10.0*(7.5+x_beg03)/15 6T3uv,2  
%N((p[\H  
    x_end04=s1/2 fL3Px  
        y_end04=1.5*s2 O>8|Lc  
    z_end04=0 &8kc0Z@y  
ly,3,ok  
     61qs`N=k  
end UO3Q
wZ4j;  
_parm 5Q_T=TL  
>%qk2h>  
;-----------------建立模型---------------------------------------------------- QGv$~A[h  
-P I$SA,  
gen zone brick p0 -17.5 0 0 p1 -7.5 0 0 p2 -17.5 s 0 p3 -17.5 0 10 size 20 3 20 ratio 1 1 1 group 21 j?!BHNs  
gen zone brick p0 7.5 0 0 p1 17.5 0 0 p2 7.5 s 0 p3 7.5 0 20 size 20 3 20 ratio 1 1 1 group 22 _<;;CI3w  
\}CQo0v  
gen zone brick p0 dslope 0 0 p1 tslope 0 0 p2 dslope s 0 p3 dslope 0 10 p4 tslope s 0  p5 dslope s 10 p6 tslope 0 20.0 p7 tslope s 20.0 size 20 3 20 ratio 1 1 1 group 12 eN*=wOh  
#jY\l&E  
group slope range group 12 any group 21 any group 22 any 05$;7xnf(  
attach face 9 Vn  
^]nnvvp  
model mohr ZUDdLJ  
prop shear soil_shear bulk soil_bulk fric s_fri coh s_coh tens s_tens #&Xr2?E@  
ini dens soil_dens +rT%C&ze  
Y&vn`#   
def pileproperty &yu3nA:7D  
   d_pile=p_d1 g&z)y  
   A_pile=pi*(d_pile)^2/4    ;截面积 c eH8  
   Iy_pile=pi*(d_pile)^4/64  ;Y轴惯性矩 Z0o+&3a6  
   Iz_pile=pi*(d_pile)^4/64  ;Z轴惯性矩 )-2sk@y  
   Ij_pile=Iy_pile+Iz_pile   ;惯性积 7Jm&z/  
   C_pile=pi*d_pile          ;周长 9\2<#,R1q  
   emod_pile=2.5e10          ;微型桩弹性模量 /-Saz29f^Q  
   nu_pile=0.2               ;微型桩泊松比 2oY.MQD7iW  
end pileproperty %}< e;t-O  
pileproperty 4J#F;#iA  
VD=}GY33=  
         sel pile id=1 begin=(x_beg01,y_beg01,z_beg01) end=(x_end01,y_end01,z_end01) nseg=15 t(<^
of:  
         sel pile id=1 begin=(x_beg02,y_beg02,z_beg02) end=(x_end02,y_end02,z_end02) nseg=15 jg2UX  
         sel pile id=1 begin=(x_beg03,y_beg03,z_beg03) end=(x_end03,y_end03,z_end03) nseg=15 K})=&<M0  
         sel pile id=1 begin=(x_beg04,y_beg04,z_beg04) end=(x_end04,y_end04,z_end04) nseg=15 cvoE4&m!  
         93Z/|7  
         sel pile id=1 prop emod=emod_pile nu=nu_pile  & +TA'P$j  
                       xcarea=A_pile per=C_pile  & f?KHp|  
                       xcj=Ij_pile xciy=Iy_pile xciz=Iz_pile  & \BIa:}9O  
                       cs_sk=1.3e11 cs_scoh=1.5e10 cs_sfric=20  & !m7`E  
                       cs_nk=1.3e9 cs_ncoh=1.0e4 cs_nfric=20 cs_ngap=off aQ@9(j> F  
   ].E89_|O  
;-------------------------边界条件-------------------------------- l/=2P_8
+Z  
jZRf{  
fix x     range x -17.6 -17.4 x2-i1#j`;  
fix x     range x 17.4 17.6 FG-v71!h#  
fix y     range y -0.1 0.1 G8]DK3#  
fix y     range y sl1 sl2 "EA6RFRD  
fix x y z range z -0.1 0.1 j$2rU'  
N?Wx-pK  
;----------------------求解安全系数FOS----------------------- z ^e99dz  
X<pg^Y0  
set grav 0,0,-10 `2}Frw+?  
>[,ywRJ#_}  
solve fos file 0.5l.sav

我也在做抗滑桩的模拟 但是根据文献的FLAC的算例进行验算 发现加抗滑桩后得不出和文献相同的结果，不知道怎么回事。有时间了和lz探讨下。