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

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

图片:00.png

  

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

命令流： >oqZ !V5[  
4vqNule  
new &O tAAE  
def _parm &jm[4'$ *z  
yXF|Sqv  
        ;-----------------几何参数---------------------------------------------------- KO,_6>8]U  
SE!L :  
        p_d1=0.15                      ;-----------桩径------------------------------- cU+>|'f&  
        tslope=7.5                     ;-----------坡顶------------------------------- [hU=mS8=^  
        dslope=-7.5                    ;-----------坡脚------------------------------- ]^@0
+!  
        ngrad01=1.5                    ;-----------坡率------------------------------- 3*WS"bt  
        s1=3.0*p_d1                    ;-----------桩排距S1=3D------------------------- 2Xgx*'t\  
        s2=3.0*p_d1                    ;-----------桩间距S2=3D---------------------- Jz$>k$!UD  
        s=2.0*s2                       ;-----------模型厚度S=2S2---------------------- q^Inb)FeN  
        sl1=s-0.1 3h&b
Z  
        sl2=s+0.1 a:u}d7T3e  
}W}G X(?P  
        ;-----------------土体材料参数---------------------------------------------------- H8i+'5x,?  
         G(\Ckf:  
        s_Emod=1.0e7 :qnokrGzB  
        s_nu=0.3 FpzP#;  
        s_fri=20.0 kCp)!hV
Q  
        s_coh=10e3 x,1=D~L}  
        s_tens=0 BXA]9eK  
        soil_dens=2000 1+^n!$  
        soil_bulk=s_Emod/3/(1-2*s_nu) 8c).8RLf  
        soil_shear=s_Emod/2/(1+s_nu) k4R4YI"jV  
MPvWCPB  
q6&67u0  
        ;-----------------桩位---------------------------------------------------- 
-f.R#J$2  
         };p~A-E=  
        ;-----------------1# 桩----------------------     )ddsyFGW  
        x_beg01=-s1/2 gj;@?o0  
        y_beg01=0.5*s2     \7 Mq $d  
        z_beg01=10.0+5.0*(7.5+x_beg01)/7.5 8IJ-]wHIb  
owMmCR  
    x_end01=-s1/2 d*%Mv[X:<  
        y_end01=0.5*s2 iGa}3pF  
    z_end01=0 $e1:Q#den2  
         :#SNpn=@  
        ;-----------------2# 桩---------------------- 1;`Fe":;vC  
        x_beg02=-s1/2 ]ymC3LV]  
        y_beg02=1.5*s2     cY~M4:
vgT  
        z_beg02=10.0+5.0*(7.5+x_beg01)/7.5 /e<5Np\X  
K;C_Z/<%  
    x_end02=-s1/2 dXO=ZU/N  
        y_end02=1.5*s2 NL-_#N$  
    z_end02=0 u _X}-U  
5)*6V&  
        ;-----------------3# 桩---------------------- !x_t`78T  
        x_beg03=s1/2 7@~QkTH~y  
        y_beg03=0.5*s2     YkI9d&ib+  
        z_beg03=10.0+10.0*(7.5+x_beg03)/15 h--45`cE  
-t>Z 9  
    x_end03=s1/2 os_WYQ4>j  
        y_end03=0.5*s2 R8uj3!3^  
    z_end03=0 Nm"<!a<F  
^$}/|d(  
        ;-----------------4# 桩---------------------- X~zRZ0  
        x_beg04=s1/2 qp@m&GH  
        y_beg04=1.5*s2     3W'FcE)|E  
        z_beg04=10.0+10.0*(7.5+x_beg03)/15 ifHU|0_=  
B~rU1Y)  
    x_end04=s1/2 a7=YG6[  
        y_end04=1.5*s2 e&ts\0  
    z_end04=0 %-r?=L  
V" }*"P-%  
     %yKKUZ~  
end Z+,CL/  
_parm

< GU  
_eh3qs:  
;-----------------建立模型---------------------------------------------------- + "zYn!0  
Oc~VHT  
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 \H[Yyp4  
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 S2bexbp0o  
~.H~XKw
 
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 ( mt*y]p?  
n5\}KZh  
group slope range group 12 any group 21 any group 22 any |C4o zl=O?  
attach face qZP:@r"  
G}q<{<+$  
model mohr %VS 2M #f  
prop shear soil_shear bulk soil_bulk fric s_fri coh s_coh tens s_tens o Fi) d[`  
ini dens soil_dens S<9gyW  
wE8]'o  
def pileproperty :^5>wD
u{  
   d_pile=p_d1 Qqh^E_O  
   A_pile=pi*(d_pile)^2/4    ;截面积 }Wqtip:L  
   Iy_pile=pi*(d_pile)^4/64  ;Y轴惯性矩 |lY`9-M`I  
   Iz_pile=pi*(d_pile)^4/64  ;Z轴惯性矩 w?i)/q  
   Ij_pile=Iy_pile+Iz_pile   ;惯性积 A`g.[7  
   C_pile=pi*d_pile          ;周长 oHk27U G  
   emod_pile=2.5e10          ;微型桩弹性模量 r;$r=Ufr  
   nu_pile=0.2               ;微型桩泊松比 wZ O@J|  
end pileproperty )u-ns5  
pileproperty 1Lc8fP$  
JC/d
:.  
         sel pile id=1 begin=(x_beg01,y_beg01,z_beg01) end=(x_end01,y_end01,z_end01) nseg=15 vv!Bo~L1,  
         sel pile id=1 begin=(x_beg02,y_beg02,z_beg02) end=(x_end02,y_end02,z_end02) nseg=15 w~p4S+k&  
         sel pile id=1 begin=(x_beg03,y_beg03,z_beg03) end=(x_end03,y_end03,z_end03) nseg=15 _=6vW^s  
         sel pile id=1 begin=(x_beg04,y_beg04,z_beg04) end=(x_end04,y_end04,z_end04) nseg=15 `LqnEutzc  
         n}f3Vrl  
         sel pile id=1 prop emod=emod_pile nu=nu_pile  & =8@RKG`>;  
                       xcarea=A_pile per=C_pile  & g(1"GKg3K  
                       xcj=Ij_pile xciy=Iy_pile xciz=Iz_pile  & d0TgqO{  
                       cs_sk=1.3e11 cs_scoh=1.5e10 cs_sfric=20  & i"vDRrDe  
                       cs_nk=1.3e9 cs_ncoh=1.0e4 cs_nfric=20 cs_ngap=off GG*BN<(>!  
   :)j7U3u  
;-------------------------边界条件-------------------------------- =m
7CJc
 
w gmWo8  
fix x     range x -17.6 -17.4 *XHj)DC;  
fix x     range x 17.4 17.6 :bI,rEW#_  
fix y     range y -0.1 0.1 G%Y*q(VrEu  
fix y     range y sl1 sl2 9Il'E6 J  
fix x y z range z -0.1 0.1 75<el.'H  
{({ R:!c  
;----------------------求解安全系数FOS----------------------- Nz}|%.GP"  
UC.8DaIPN  
set grav 0,0,-10 w~ijD ^g  
1]HHe*'Z  
solve fos file 0.5l.sav

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