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

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

图片:00.png

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

命令流： foE2rV/Y  
:ykZ7X&  
new n]coqJ  
def _parm sS(^7GARa  
8yFD2(#  
        ;-----------------几何参数---------------------------------------------------- _zm<[0(  
INzQ0z-z  
        p_d1=0.15                      ;-----------桩径------------------------------- }+#ag:M  
        tslope=7.5                     ;-----------坡顶------------------------------- Q:VD2<2  
        dslope=-7.5                    ;-----------坡脚------------------------------- qm]ljut  
        ngrad01=1.5                    ;-----------坡率------------------------------- ,bmTBZV  
        s1=3.0*p_d1                    ;-----------桩排距S1=3D------------------------- wQnr*kyza  
        s2=3.0*p_d1                    ;-----------桩间距S2=3D---------------------- 8gBqur{  
        s=2.0*s2                       ;-----------模型厚度S=2S2---------------------- K{>O.5  
        sl1=s-0.1 +I\bs.84  
        sl2=s+0.1 h/(9AO}t  
o(~JZik  
        ;-----------------土体材料参数---------------------------------------------------- 3[aJ=5  
         P!YT{}  
        s_Emod=1.0e7 cX.v^9kuX  
        s_nu=0.3 G';oM;~/|  
        s_fri=20.0 a/^YgrC\T  
        s_coh=10e3 9i|6  
        s_tens=0 x'JfRz  
        soil_dens=2000 +bf%]   
        soil_bulk=s_Emod/3/(1-2*s_nu) 2#W%--  
        soil_shear=s_Emod/2/(1+s_nu) |klL KX&  
a9jY^E'|n  
V|?  
        ;-----------------桩位---------------------------------------------------- p7H*Ff`  
         F<-Pbtw  
        ;-----------------1# 桩----------------------     05pCgI}F>  
        x_beg01=-s1/2 n7<<}wcV  
        y_beg01=0.5*s2     Z@C D1+G  
        z_beg01=10.0+5.0*(7.5+x_beg01)/7.5 L1C'V/g  
-R7f/a8  
    x_end01=-s1/2 [TO:-8$.  
        y_end01=0.5*s2 R?|_`@@A  
    z_end01=0 ;co{bk|rj  
         N}FG%a  
        ;-----------------2# 桩---------------------- D|-]"(2i  
        x_beg02=-s1/2 YN>#zr+~  
        y_beg02=1.5*s2     !8=uBS%  
        z_beg02=10.0+5.0*(7.5+x_beg01)/7.5 ?QVD)JI*k  
x|<|eRYK  
    x_end02=-s1/2 /e{Oqhf[n  
        y_end02=1.5*s2 &|E2L1  
    z_end02=0 ( v ~/glf  
N{p2@_fnB  
        ;-----------------3# 桩---------------------- Z^GriL  
        x_beg03=s1/2 <O\z`aA'q  
        y_beg03=0.5*s2     @1Zf&'/6  
        z_beg03=10.0+10.0*(7.5+x_beg03)/15 V(!b!i@  
x=au.@psBS  
    x_end03=s1/2 29Q5s$YD@  
        y_end03=0.5*s2 <"[}8  
    z_end03=0 'OX6eY5
 
Dh +^;dQ6  
        ;-----------------4# 桩---------------------- J?%D4AeS]v  
        x_beg04=s1/2 PL+fLCk,I  
        y_beg04=1.5*s2     ^<|If:|  
        z_beg04=10.0+10.0*(7.5+x_beg03)/15 )s=z i"  
bR&hI9`%F  
    x_end04=s1/2 tfv]AC7x  
        y_end04=1.5*s2 dnIBAe  
    z_end04=0 B4|%E$1+  
g\*gHHa  
     MB"?^~Sm  
end P<4jY?.  
_parm Va*Uwy?x/)  
BTd'bD~EA  
;-----------------建立模型---------------------------------------------------- s9[v_(W  
LK:|~UV?  
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 '{oe}].,
 
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 6o1.?t?  
(c[h,>`@:  
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 QdW%5lM+  
*.nqQhW  
group slope range group 12 any group 21 any group 22 any DD3J2J  
attach face ^*{xTB57  
w@%W{aUC  
model mohr @#Xzk?+  
prop shear soil_shear bulk soil_bulk fric s_fri coh s_coh tens s_tens ;:$Na=  
ini dens soil_dens Zp__  
o!\
O)  
def pileproperty acGmRP9g  
   d_pile=p_d1 ]B,S<*h  
   A_pile=pi*(d_pile)^2/4    ;截面积 $yFur[97C  
   Iy_pile=pi*(d_pile)^4/64  ;Y轴惯性矩 b0t];Gc%b  
   Iz_pile=pi*(d_pile)^4/64  ;Z轴惯性矩 MzG(+B  
   Ij_pile=Iy_pile+Iz_pile   ;惯性积 F~l3?3ZV  
   C_pile=pi*d_pile          ;周长 *iC t4J  
   emod_pile=2.5e10          ;微型桩弹性模量 B-&J]H  
   nu_pile=0.2               ;微型桩泊松比 -ZyFUGd%  
end pileproperty P^h2w%6'  
pileproperty ([9h.M6v  
7L-%5:1%  
         sel pile id=1 begin=(x_beg01,y_beg01,z_beg01) end=(x_end01,y_end01,z_end01) nseg=15 *10e)rzM  
         sel pile id=1 begin=(x_beg02,y_beg02,z_beg02) end=(x_end02,y_end02,z_end02) nseg=15 p6HZ2Q:a  
         sel pile id=1 begin=(x_beg03,y_beg03,z_beg03) end=(x_end03,y_end03,z_end03) nseg=15 SV\x2^Ea0  
         sel pile id=1 begin=(x_beg04,y_beg04,z_beg04) end=(x_end04,y_end04,z_end04) nseg=15 ?pF;{  
         10}Zoq|)n  
         sel pile id=1 prop emod=emod_pile nu=nu_pile  & \ I?;%  
                       xcarea=A_pile per=C_pile  & hCxL4LrF  
                       xcj=Ij_pile xciy=Iy_pile xciz=Iz_pile  & )
uX:f8  
                       cs_sk=1.3e11 cs_scoh=1.5e10 cs_sfric=20  & giu~"#0/F  
                       cs_nk=1.3e9 cs_ncoh=1.0e4 cs_nfric=20 cs_ngap=off U.^)|IHW  
   3.dUMJ$_  
;-------------------------边界条件-------------------------------- h;ShNU  
jZ{S{"j  
fix x     range x -17.6 -17.4 >$Fc=~;Ba  
fix x     range x 17.4 17.6 |
[{;*wtv  
fix y     range y -0.1 0.1 mML^kgy\N  
fix y     range y sl1 sl2 GO?-z0V  
fix x y z range z -0.1 0.1 U<6k!Y9ny  

r 2  
;----------------------求解安全系数FOS----------------------- dl":?D4H  
lP9I\Ge&  
set grav 0,0,-10 'g=yJ  
VhW;=y>}  
solve fos file 0.5l.sav

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