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

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

图片:00.png

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

命令流： "> ]{t[Ib  
9o<5Z=  
new a\]glw\;  
def _parm -+)06BqF}  
XhEJF !  
        ;-----------------几何参数---------------------------------------------------- VxTrL}{(6  
4apy{W  
        p_d1=0.15                      ;-----------桩径------------------------------- $AyE6j_1gX  
        tslope=7.5                     ;-----------坡顶------------------------------- /Y_)dz^@  
        dslope=-7.5                    ;-----------坡脚------------------------------- 8Ts_;uId  
        ngrad01=1.5                    ;-----------坡率------------------------------- +B#+'  
        s1=3.0*p_d1                    ;-----------桩排距S1=3D------------------------- \:Vm7Zg  
        s2=3.0*p_d1                    ;-----------桩间距S2=3D---------------------- d:&=|k
Kw  
        s=2.0*s2                       ;-----------模型厚度S=2S2---------------------- aEvW<jHh  
        sl1=s-0.1 q:a-tdv2  
        sl2=s+0.1 d")TH3pG  
c80!Ub@  
        ;-----------------土体材料参数---------------------------------------------------- 2Q81#i'Cm  
         5N /NUs   
        s_Emod=1.0e7 b2vCr F;  
        s_nu=0.3 Y)*:'&~2e  
        s_fri=20.0 3<A$lG  
        s_coh=10e3 4mM?RGWv
 
        s_tens=0 =+ vl+h  
        soil_dens=2000 qC:QY6g$N  
        soil_bulk=s_Emod/3/(1-2*s_nu) SpJIEw  
        soil_shear=s_Emod/2/(1+s_nu) =b_/_b$q  
*i,@d&J y]  
=6H  
        ;-----------------桩位---------------------------------------------------- NR9=V  
         FhJtiw@  
        ;-----------------1# 桩----------------------     7N4)T'B  
        x_beg01=-s1/2 ?, m_q+  
        y_beg01=0.5*s2     /|LQ?n  
        z_beg01=10.0+5.0*(7.5+x_beg01)/7.5 ,vs#(d6G  
Y{D?
&x%yq  
    x_end01=-s1/2 (U([T-H  
        y_end01=0.5*s2 VxW>XxG0  
    z_end01=0 \IX|{]*D  
         34c+70x7  
        ;-----------------2# 桩---------------------- 2e^6Od!Y?  
        x_beg02=-s1/2 *6/OLAkyF  
        y_beg02=1.5*s2     j^6,V\;l  
        z_beg02=10.0+5.0*(7.5+x_beg01)/7.5 b"pN;v  
EV[ BB;eb  
    x_end02=-s1/2 9][A1+"  
        y_end02=1.5*s2 apwA  
    z_end02=0 bR*} s/  
M~l\rg8  
        ;-----------------3# 桩---------------------- -58Sb"f  
        x_beg03=s1/2 ?%ei+  
        y_beg03=0.5*s2     FM,o&0HSd  
        z_beg03=10.0+10.0*(7.5+x_beg03)/15 zT+ "Z(oz,  
e!N%  
    x_end03=s1/2 uFnq3m^u  
        y_end03=0.5*s2 <Gj]XAoe%  
    z_end03=0 "?S>}G\  
[Wn6d:  
        ;-----------------4# 桩---------------------- <tgfbY^nL  
        x_beg04=s1/2 Wxzh'c#\8  
        y_beg04=1.5*s2     OLE@35"v]  
        z_beg04=10.0+10.0*(7.5+x_beg03)/15 "sJ@_lp  
kJ"}JRA<  
    x_end04=s1/2 Jix;!("  
        y_end04=1.5*s2 _#y=T20'3  
    z_end04=0 Z+agS8e(  
qk=OodEMK  
     1xjw=  
end UHS{X~CS e  
_parm HE@P<  
gr\vC  
;-----------------建立模型---------------------------------------------------- <(_${zR  
{nH*Wu*^  
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 R_gON*9  
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 IeAUVRS)  
?
P"j5  
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 N**"u"CX  
C^ZDUj`  
group slope range group 12 any group 21 any group 22 any rqEP!S^  
attach face Dpkc9~z  
cv9-ZOxJ  
model mohr wpa^]l  
prop shear soil_shear bulk soil_bulk fric s_fri coh s_coh tens s_tens )NnkoCNeE  
ini dens soil_dens zTg&W7oz  
So]O`RJv  
def pileproperty J<-2dvq  
   d_pile=p_d1 \ mg  
   A_pile=pi*(d_pile)^2/4    ;截面积 9&Ny;oy#6  
   Iy_pile=pi*(d_pile)^4/64  ;Y轴惯性矩 (``|5;T\  
   Iz_pile=pi*(d_pile)^4/64  ;Z轴惯性矩 FB n . 4  
   Ij_pile=Iy_pile+Iz_pile   ;惯性积 Qk8YR5K   
   C_pile=pi*d_pile          ;周长 l-M .
C8N  
   emod_pile=2.5e10          ;微型桩弹性模量 S<>u  
   nu_pile=0.2               ;微型桩泊松比 tx]!|x" F  
end pileproperty Q#}c5TjVr  
pileproperty _F3= H]P  
M4K>/-9X+V  
         sel pile id=1 begin=(x_beg01,y_beg01,z_beg01) end=(x_end01,y_end01,z_end01) nseg=15 >/7[HhBT  
         sel pile id=1 begin=(x_beg02,y_beg02,z_beg02) end=(x_end02,y_end02,z_end02) nseg=15 :JCe,1!3@  
         sel pile id=1 begin=(x_beg03,y_beg03,z_beg03) end=(x_end03,y_end03,z_end03) nseg=15 X-_0wR  
         sel pile id=1 begin=(x_beg04,y_beg04,z_beg04) end=(x_end04,y_end04,z_end04) nseg=15 Jqqt@5Ni  
         5mg] su&#  
         sel pile id=1 prop emod=emod_pile nu=nu_pile  & g&d tOjM  
                       xcarea=A_pile per=C_pile  & XT\Q"=FD  
                       xcj=Ij_pile xciy=Iy_pile xciz=Iz_pile  & ,e+S7YX  
                       cs_sk=1.3e11 cs_scoh=1.5e10 cs_sfric=20  & V?4G~~F  
                       cs_nk=1.3e9 cs_ncoh=1.0e4 cs_nfric=20 cs_ngap=off s`o_ER  
   :.?%e{7  
;-------------------------边界条件-------------------------------- B^BbA-I  
U_AmRiy  
fix x     range x -17.6 -17.4 E7j9A`  
fix x     range x 17.4 17.6 T@j@IEGH  
fix y     range y -0.1 0.1 y-vQ4G5F|  
fix y     range y sl1 sl2 e'&<DE)  
fix x y z range z -0.1 0.1 ]`\~(*;[W9  
8~vE  
;----------------------求解安全系数FOS----------------------- 6adXE  
,{s
CI
/  
set grav 0,0,-10 =
R0#WMf$@  
"U iv[8B  
solve fos file 0.5l.sav

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