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

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

图片:00.png

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

命令流： )vB,eZq  
}| BnG"8  
new 1>$fLbmkI  
def _parm @<5?q:9.8  
=IH z@CU  
        ;-----------------几何参数---------------------------------------------------- "=0#pH1o  
&<Bx1\ ~V  
        p_d1=0.15                      ;-----------桩径------------------------------- 8S_i;  
        tslope=7.5                     ;-----------坡顶------------------------------- L
`%v#R  
        dslope=-7.5                    ;-----------坡脚------------------------------- #\kYGr-G)  
        ngrad01=1.5                    ;-----------坡率------------------------------- 9|Cu2  
        s1=3.0*p_d1                    ;-----------桩排距S1=3D------------------------- %Y"@VcN  
        s2=3.0*p_d1                    ;-----------桩间距S2=3D---------------------- kDM\IyM<\  
        s=2.0*s2                       ;-----------模型厚度S=2S2---------------------- J;pn5k~3  
        sl1=s-0.1 ULq#2l  
        sl2=s+0.1 K4Mv\!Q<8  
d>z?JDt  
        ;-----------------土体材料参数---------------------------------------------------- ~l~ai>/  
         =6Dz<Lq  
        s_Emod=1.0e7 L3^WI( 8m  
        s_nu=0.3 \Lxsg!wtJ  
        s_fri=20.0  hfpSxL  
        s_coh=10e3 vDgf}  
        s_tens=0 %a;N)1/  
        soil_dens=2000 ] 8sVXZ  
        soil_bulk=s_Emod/3/(1-2*s_nu) _C3l2v'I$  
        soil_shear=s_Emod/2/(1+s_nu) #WD}XOA  
P>/n!1c  
fHek!Jv.  
        ;-----------------桩位---------------------------------------------------- \e%H5Wx
  
         uUXvBA?l  
        ;-----------------1# 桩----------------------     \vVGfG?6  
        x_beg01=-s1/2 K~p\B  
        y_beg01=0.5*s2     zK`z*\
 
        z_beg01=10.0+5.0*(7.5+x_beg01)/7.5 ENwDW#U9  
\K+LKa)  
    x_end01=-s1/2 A ^YHtJ  
        y_end01=0.5*s2 }v[*V 
 
    z_end01=0 i?uJ<BdU[  
         ,}i`1E1=  
        ;-----------------2# 桩---------------------- SG1fu<Q6J  
        x_beg02=-s1/2 Z}(,OZh  
        y_beg02=1.5*s2     eJ+V!K'H2  
        z_beg02=10.0+5.0*(7.5+x_beg01)/7.5 `}/&}Sp  
3+gp_7L  
    x_end02=-s1/2 VY)!bjW.  
        y_end02=1.5*s2 F.)b`:g  
    z_end02=0 FWNO/)~t  
PT7L65  
        ;-----------------3# 桩---------------------- #!d@;=[\  
        x_beg03=s1/2 \8%64ZL`  
        y_beg03=0.5*s2     u?[dy n  
        z_beg03=10.0+10.0*(7.5+x_beg03)/15 )K2,h5zU  
+5Yf9  
    x_end03=s1/2 F0O"rN{  
        y_end03=0.5*s2 a $pxt!6  
    z_end03=0 Qo>VN`v  
<4,n6$E  
        ;-----------------4# 桩---------------------- +;7Rz_.6f  
        x_beg04=s1/2 H tIl;E  
        y_beg04=1.5*s2     4-@D`,3L  
        z_beg04=10.0+10.0*(7.5+x_beg03)/15 Fv \yhR  
iV5x-G`  
    x_end04=s1/2 9ZjSM,+  
        y_end04=1.5*s2 H-GlCVq~  
    z_end04=0 XkZ82w#b  
irSdqa/  
     K4]c   
end : [o0Va2 d  
_parm 9/[3xhB4  
k23*F0Dv  
;-----------------建立模型---------------------------------------------------- lg_X|yhL  
Vk/CV2  
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 0*S2_&Q)  
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 b&BSigrvou  
gbOd(ugH  
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 +@),
Fk_  
Np>[mNmga  
group slope range group 12 any group 21 any group 22 any [ay~l%x  
attach face RkVU^N"  
iX0s4  
model mohr P+!j[X^  
prop shear soil_shear bulk soil_bulk fric s_fri coh s_coh tens s_tens : E`N0UA  
ini dens soil_dens )jbYWR*&  
ZOGH.`  
def pileproperty N5u.V\F!z\  
   d_pile=p_d1 [m7^Euury  
   A_pile=pi*(d_pile)^2/4    ;截面积 l?:!G7ie  
   Iy_pile=pi*(d_pile)^4/64  ;Y轴惯性矩 8<}f:9/  
   Iz_pile=pi*(d_pile)^4/64  ;Z轴惯性矩 XJl2_
#  
   Ij_pile=Iy_pile+Iz_pile   ;惯性积 |7Z7_YWs  
   C_pile=pi*d_pile          ;周长 *rPUVhD_  
   emod_pile=2.5e10          ;微型桩弹性模量 @[M5$,"  
   nu_pile=0.2               ;微型桩泊松比 ?[>+'6  
end pileproperty &]gw[ `  
pileproperty wykk</eQ.i  
1df}gG  
         sel pile id=1 begin=(x_beg01,y_beg01,z_beg01) end=(x_end01,y_end01,z_end01) nseg=15 -=aI!7*"$  
         sel pile id=1 begin=(x_beg02,y_beg02,z_beg02) end=(x_end02,y_end02,z_end02) nseg=15 +$Q33@F5l  
         sel pile id=1 begin=(x_beg03,y_beg03,z_beg03) end=(x_end03,y_end03,z_end03) nseg=15 V_JM@VN}Kk  
         sel pile id=1 begin=(x_beg04,y_beg04,z_beg04) end=(x_end04,y_end04,z_end04) nseg=15 J,ZvaF  
         t0XM#9L  
         sel pile id=1 prop emod=emod_pile nu=nu_pile  & 9V*h:[6a(  
                       xcarea=A_pile per=C_pile  & Xk[;MZ[  
                       xcj=Ij_pile xciy=Iy_pile xciz=Iz_pile  & ZSj^\JU  
                       cs_sk=1.3e11 cs_scoh=1.5e10 cs_sfric=20  & \M>}-j`v  
                       cs_nk=1.3e9 cs_ncoh=1.0e4 cs_nfric=20 cs_ngap=off y#HDJ=2  
   f.ku v"  
;-------------------------边界条件-------------------------------- wehZ7eqm  
FCv3ZF?K  
fix x     range x -17.6 -17.4 "Gx(-NH+  
fix x     range x 17.4 17.6 5bZf$$b  
fix y     range y -0.1 0.1 5#+G7 'k  
fix y     range y sl1 sl2 #gbJ$1s  
fix x y z range z -0.1 0.1 
eIjn~2^  
`z<k7ig  
;----------------------求解安全系数FOS----------------------- b_xn80O  
/PbN!r<1  
set grav 0,0,-10 pj.}VF!d  
iOI8'`mk  
solve fos file 0.5l.sav

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