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

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

图片:00.png

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

命令流： @0z0m;8  
nuSN)}b<Q  
new Ug7`ez4vw  
def _parm TF=k(@9J?  
3qiJwo
>  
        ;-----------------几何参数---------------------------------------------------- q9^Y?`  
c6,s+^^  
        p_d1=0.15                      ;-----------桩径------------------------------- l Io9,K
e  
        tslope=7.5                     ;-----------坡顶------------------------------- A<SOT>m]  
        dslope=-7.5                    ;-----------坡脚------------------------------- 1l+kO,X]  
        ngrad01=1.5                    ;-----------坡率------------------------------- 5L-lpT8P  
        s1=3.0*p_d1                    ;-----------桩排距S1=3D------------------------- ACigeK^C}E  
        s2=3.0*p_d1                    ;-----------桩间距S2=3D---------------------- d&|z=%9xl  
        s=2.0*s2                       ;-----------模型厚度S=2S2---------------------- 6F*-qb3  
        sl1=s-0.1 heL$2dZ5H  
        sl2=s+0.1 Tr8AG>  
y9C;T(oi;  
        ;-----------------土体材料参数---------------------------------------------------- 1E
5a(  
         "x(>Sj\%I  
        s_Emod=1.0e7 _[OF"X2  
        s_nu=0.3 ~h)@e\Kc  
        s_fri=20.0 f]tv`<Q7  
        s_coh=10e3 lt{lpH  
        s_tens=0 R9Sf!LR  
        soil_dens=2000 B;Xoa,  
        soil_bulk=s_Emod/3/(1-2*s_nu) ><NI'q*cQ  
        soil_shear=s_Emod/2/(1+s_nu) <0u\dU  
A%Bgp?B  
x&C%4Y_]  
        ;-----------------桩位---------------------------------------------------- +0XL5('2  
         e1 x^PT  
        ;-----------------1# 桩----------------------     `^7:7Wr]=  
        x_beg01=-s1/2 wMb)6YZs  
        y_beg01=0.5*s2     -t8hi+NK  
        z_beg01=10.0+5.0*(7.5+x_beg01)/7.5 ,3Wb4so  
zL:&Q<  
    x_end01=-s1/2 jR{-  
        y_end01=0.5*s2 w9RF2J  
    z_end01=0 .
dx 4,|6  
         bTLMd$  
        ;-----------------2# 桩---------------------- Jd5\&ma  
        x_beg02=-s1/2 k"xGA*B|  
        y_beg02=1.5*s2     pOC% oj  
        z_beg02=10.0+5.0*(7.5+x_beg01)/7.5 \  Md 3  
}p~OCW!  
    x_end02=-s1/2 '|l%rv  
        y_end02=1.5*s2 3k#[(phk  
    z_end02=0 sl/=g   
#D+7TWDwNt  
        ;-----------------3# 桩---------------------- Q [rj  
        x_beg03=s1/2 ?y~TCqV  
        y_beg03=0.5*s2     @#RuSc  
        z_beg03=10.0+10.0*(7.5+x_beg03)/15 oZTgN .q  
&v^!y=Bt  
    x_end03=s1/2 I eG=J4:*  
        y_end03=0.5*s2 vQ:wW',i  
    z_end03=0 G' Blp  
E7Ibp79}N  
        ;-----------------4# 桩---------------------- !z11" c  
        x_beg04=s1/2 j-7u>s-l  
        y_beg04=1.5*s2     +I t#Z3  
        z_beg04=10.0+10.0*(7.5+x_beg03)/15 v" }WP34  
G&q'#3ieC  
    x_end04=s1/2 1/B]TT  
        y_end04=1.5*s2 wfgqgPo!v  
    z_end04=0 ?4XnEDAm  
pb!V|#u"  
     aaDP9FW9e  
end dj'm, k b  
_parm GCDwWCxh  
pSHSgd~&  
;-----------------建立模型---------------------------------------------------- wV(AT$  
z{9=1XY  
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 M)&Io6>  
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 w|IjQ1{  
NXpmT4  
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 veeI==]  
>F1G!#$0  
group slope range group 12 any group 21 any group 22 any *G9sy_  
attach face LL&ud_Y  
qO-9 x0v#  
model mohr X) V7bVW  
prop shear soil_shear bulk soil_bulk fric s_fri coh s_coh tens s_tens [4sEVu}  
ini dens soil_dens 9Ycn0  
7R}9oK_I  
def pileproperty R}8XRe  
   d_pile=p_d1 
tz26=8  
   A_pile=pi*(d_pile)^2/4    ;截面积 |EKu2We*  
   Iy_pile=pi*(d_pile)^4/64  ;Y轴惯性矩 ,57$N&w  
   Iz_pile=pi*(d_pile)^4/64  ;Z轴惯性矩 kb71q:[  
   Ij_pile=Iy_pile+Iz_pile   ;惯性积 >M]6uf  
   C_pile=pi*d_pile          ;周长 `vd= ec  
   emod_pile=2.5e10          ;微型桩弹性模量 {( #zcK  
   nu_pile=0.2               ;微型桩泊松比 o*">KqU`b  
end pileproperty x2co>.
i  
pileproperty j~;;l!({i  
hoa7   
         sel pile id=1 begin=(x_beg01,y_beg01,z_beg01) end=(x_end01,y_end01,z_end01) nseg=15 zN#*G i'  
         sel pile id=1 begin=(x_beg02,y_beg02,z_beg02) end=(x_end02,y_end02,z_end02) nseg=15 Mi+H#xx16  
         sel pile id=1 begin=(x_beg03,y_beg03,z_beg03) end=(x_end03,y_end03,z_end03) nseg=15 +#2)kg 9_  
         sel pile id=1 begin=(x_beg04,y_beg04,z_beg04) end=(x_end04,y_end04,z_end04) nseg=15 1B|8ZmFJj  
         e,>%Z@92(  
         sel pile id=1 prop emod=emod_pile nu=nu_pile  &  v,=v  
                       xcarea=A_pile per=C_pile  & jjRUL.  
                       xcj=Ij_pile xciy=Iy_pile xciz=Iz_pile  & pY@Y?Jj  
                       cs_sk=1.3e11 cs_scoh=1.5e10 cs_sfric=20  & _A98  
                       cs_nk=1.3e9 cs_ncoh=1.0e4 cs_nfric=20 cs_ngap=off ~vHk
&r]|  
   7p.>\YtoR}  
;-------------------------边界条件-------------------------------- "13"`!m  
> _ <'D  
fix x     range x -17.6 -17.4 {:Aw_z:'  
fix x     range x 17.4 17.6 ;}qhc l+  
fix y     range y -0.1 0.1 L*h{'<Bz  
fix y     range y sl1 sl2 [}OgSP9i  
fix x y z range z -0.1 0.1 ndink$  
aTd D`h  
;----------------------求解安全系数FOS----------------------- "g>.{E5  
~e `Bq>  
set grav 0,0,-10 {`(>O"_[Q  
5c5oSy+  
solve fos file 0.5l.sav

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