相关FLAC3D程序命令流(--单桩承载力)

相关FLAC程序命令流单桩竖向受力模拟newtitle竖向加载桩;桩长11m,直径0.6m;建立土体网格gen zo radcyl p0 (0 0 0) p1 (8 0 0 ) p2 (0 0 -1.5) p3 (0 8 0) & P4 (8 0 -1.5) P5 (0 8 -1.5) P6 (8 8 0) P7 (8 8 -1.5) & P8 (0.3 0 0) P9 (0 0.3 0) P10 (0.3 0 -1.5) P11 (0 0.3 -1.5) &size 3 3 6 15 ratio 1 1 1 1.15 ;粉质粘土2group fznt2 range z -1.5 0gen zo radcyl p0 (0 0 -1.5) p1 (8 0 -1.5) p2 (0 0 -3.0) p3 (0 8 -1.5) & p4 (8 0 -3.0) p5 (0 8 -3.0) p6 (8 8 -1.5) p7 (8 8 -3.0) & p8 (0.3 0 -1.5) p9 (0 0.3 -1.5) p10 (0.3 0 -3.0) p11 (0 0.3 -3.0) &size 3 3 6 15 ratio 1 1 1 1.15 ;粉细砂group fxs range z -3.0 -1.5gen zo radcyl p0 (0 0 -3.0) p1 (8 0 -3.0) p2 (0 0 -5.5) p3 (0 8 -3.0) & p4 (8 0 -5.5) p5 (0 8 -5.5) p6 (8 8 -3.0) p7 (8 8 -5.5) & p8 (0.3 0 -3.0) p9 (0 0.3 -3.0) p10 (0.3 0 -5.5) p11 (0 0.3 -5.5) &size 3 5 6 15 ratio 1 1 1 1.15 ;粉质粘土4group fznt4 range z -5.5 -3.0gen zo radcyl p0 (0 0 -5.5) p1 (8 0 -5.5) p2 (0 0 -7.0) p3 (0 8 -5.5) & p4 (8 0 -7.0) p5 (0 8 -7.0) p6 (8 8 -5.5) p7 (8 8 -7.0) & p8 (0.3 0 -5.5) p9 (0 0.3 -5.5) p10 (0.3 0 -7.0) p11 (0 0.3 -7.0) &size 3 3 6 15 ratio 1 1 1 1.15 ;粉质粘土5group fznt5 range z -7.0 -5.5gen zo radcyl p0 (0 0 -7.0) p1 (8 0 -7.0) p2 (0 0 -8.5) p3 (0 8 -7.0) & p4 (8 0 -8.5) p5 (0 8 -8.5) p6 (8 8 -7.0) p7 (8 8 -8.5) & p8 (0.3 0 -7.0) p9 (0 0.3 -7.0) p10 (0.3 0 -8.5) p11 (0 0.3 -8.5) &size 3 3 6 15 ratio 1 1 1 1.15 ;粉土group ft range z -8.5 -7.0gen zo radcyl p0 (0 0 -8.5) p1 (8 0 -8.5) p2 (0 0 -11) p3 (0 8 -8.5) & p4 (8 0 -11) p5 (0 8 -11) p6 (8 8 -8.5) p7 (8 8 -11) & p8 (0.3 0 -8.5) p9 (0 0.3 -8.5) p10 (0.3 0 -11) p11 (0 0.3 -11) &size 3 5 6 15 ratio 1 1 1 1.15 ;细砂1group xs1 range z -11 -8.5;桩体以下网格gen zo radcyl p0 (0 0 -11.0) p1 (8 0 -11.0) p2 (0 0 -14.5) p3 (0 8 -11.0) & p4 (8 0 -14.5) p5 (0 8 -14.5) p6 (8 8 -11.0) p7 (8 8 -14.5) & p8 (0.3 0 -11.0) p9 (0 0.3 -11.0) p10 (0.3 0 -14.5) p11 (0 0.3 -14.5) &size 3 7 6 15 ratio 1 1 1 1.15 fill ;细砂2group xs2 range z -14.5 -11.0gen zo radcyl p0 (0 0 -14.5) p1 (8 0 -14.5) p2 (0 0 -15.0) p3 (0 8 -14.5) & p4 (8 0 -15.0) p5 (0 8 -15.0) p6 (8 8 -14.5) p7 (8 8 -15.0) & p8 (0.3 0 -14.5) p9 (0 0.3 -14.5) p10 (0.3 0 -15.0) p11 (0 0.3 -15.0) &size 3 1 6 15 ratio 1 1 1 1.15 fill ;粉质粘土6group fznt6 range z -15.0 -14.5gen zo reflect dd 270 dip 90;设置分界面interface 1 face range cylin end1 0 0 0 end2 0 0 -1.5 radius 0.31 & cylin end1 0 0 0 end2 0 0 -1.5 radius 0.29 notinterface 2 face range cylin end1 0 0 -1.5 end2 0 0 -3.0 radius 0.31 & cylin end1 0 0 -1.5 end2 0 0 -3.0 radius 0.29 notinterface 3 face range cylin end1 0 0 -3.0 end2 0 0 -5.5 radius 0.31 & cylin end1 0 0 -3.0 end2 0 0 -5.5 radius 0.29 notinterface 4 face range cylin end1 0 0 -5.5 end2 0 0 -7.0 radius 0.31 & cylin end1 0 0 -5.5 end2 0 0 -7.0 radius 0.29 notinterface 5 face range cylin end1 0 0 -7.0 end2 0 0 -8.50 radius 0.31 & cylin end1 0 0 -7.0 end2 0 0 -8.50 radius 0.29 notinterface 6 face range cylin end1 0 0 -8.50 end2 0 0 -11 radius 0.31 & cylin end1 0 0 -8.50 end2 0 0 -11 radius 0.29 notinterface 7 face range cylin end1 0 0 -10.90 end2 0 0 -11.10 radius 0.31 ;建立桩体网格gen zo cyl p0 (0 0 12) p1 (0.3 0 12) p2 (0 0 1) p3 (0 0.3 12) & p4(0.3 0 1) p5(0 0.3 1) &size 3 22 6gen zo reflect dd 270 dip 90 range z 1 12.0group pile range z 1 12.1ini z add -12.0 range group pileattach face;设置材料属性mo elas range group pile mo mo range group fznt2prop bulk 6.349e6 shear 3.101e6 range z -1.5 0prop coh 14.7e3 fric 21.7 range group fznt2mo mo range group fxsprop bulk 6.818e6 shear 3.516e6 range z -3.0 -1.5prop coh 0 fric 25 range group fxsmo mo range group fznt4prop bulk 6.944e6 she 3.968e6 range z -5.5 -3.0prop coh 14.0e3 fric 20 range group fznt4mo mo range group fznt5prop bulk 7.639e6 she 4.365e6 range z -7.0 -5.5prop coh 15.2e3 fric 23 range group fznt5mo mo range group ftprop bulk 8.472e6 she 5.984e6 range z -8.5 -7.0prop coh 12.0e3 fric 24.5 range group ftmo mo range group xs1prop bulk 20.556e6 she 6.175e6 range z -11 -8.5prop coh 0.0 fric 18 range group xs1mo mo range group xs2prop bulk 20.556e6 she 6.175e6 range z -14.5 -11prop coh 0.0 fric 18 range group xs2mo mo range group fznt6prop bulk 25.510e6 she 9.398e6 range z -15.0 -14.5prop coh 25e3 fric 25 range group fznt6; 分界面性质interface 1 prop kn 2.097e13 ks 2.097e13 fric 20 coh 30000interface 2 prop kn 2.301e13 ks 2.301e13 fric 20 coh 30000interface 3 prop kn 2.447e13 ks 2.447e13 fric 20 coh 30000interface 4 prop kn 2.692e13 ks 2.692e13 fric 20 coh 30000interface 5 prop kn 1.223e13 ks 1.223e13 fric 20 coh 30000interface 6 prop kn 1.958e13 ks 1.958e13 fric 20 coh 30000interface 7 prop kn 7.608e13 ks 7.608e13 fric 20 coh 30000ini dens 1800 range group fznt2ini dens 1800 range group fxsini dens 1800 range group fznt4ini dens 1800 range group fznt5ini dens 1800 range group ftini dens 1800 range group xs1ini dens 1800 range group xs2ini dens 1800 range group fznt6ini dens 1800 range group pile ;边界条件fix z range z -15.1 -14.9fix x range x -8.1 -7.9fix x range x 7.9 8.1fix y range y -0.1 0.1fix y range y 7.9 8.1set grav 0 0 -9.80;初始条件ini szz 0.0 grad 0 0 18000 range z -8.5 0 ini szz 10200 grad 0 0 19200 range z -15.0 -8.5ini sxx 0.0 grad 0 0 6000 range z -8.5 0ini sxx 3400 grad 0 0 6400 range z -15.0 -8.5ini syy 0.0 grad 0 0 6000 range z -8.5 0ini syy 3400 grad 0 0 6400 range z -15.0 -8.5water density 1000water table origin 0 0 -8.5 normal 0 0 -1ini dens 1920 range z -15.0 -8.5hist unbalsolve save pile0.sav; 设桩以后初始平衡mo elas range group pileprop bulk 15.56e9 she 11.67e9 range group pileini dens 2500 range group pilesolve rat 1.0e-6save pile1.sav; 分级加载ini state 0ini xdis 0 ydis 0 zdis 0apply szz -0.7077e6 range z -0.1 0.1 group pilehist gp zdis 0 0 0hist gp zdis 0 0 -11hist zone szz 0 0 0hist zone szz 0 0 -1hist zone szz 0 0 -2hist zone szz 0 0 -3hist zone szz 0 0 -4hist zone szz 0 0 -5hist zone szz 0 0 -6hist zone szz 0 0 -7hist zone szz 0 0 -8hist zone szz 0 0 -9hist zone szz 0 0 -10hist zone szz 0 0 -11solve save pile200.savapply szz -1.0616e6 range z -0.1 0.1 group pilesolvesave pile300.savapply szz -1.4154e6 range z -0.1 0.1 group pilesolvesave pile400.savapply szz -1.7693e6 range z -0.1 0.1 group pilesolvesave pile500.savapply szz -2.1231e6 range z -0.1 0.1 group pilesolvesave pile600.savapply szz -2.4770e6 range z -0.1 0.1 group pilesolvesave pile700.savapply szz -2.8309e6 range z -0.1 0.1 group pilesolvesave pile800.savapply szz -3.1847e6 range z -0.1 0.1 group pilesolvesave pile900.savapply szz -3.5386e6 range z -0.1 0.1 group pilesolvesave pile1000.sav单桩水平受力模拟newtitle侧向加载桩gen zo b size 16 2 16 p0 0 0 0 p1 16 0 0 p2 0 0 -2 p3 0 16 0group fznt2 range z -2 0gen zo b size 16 2 16 p0 0 0 -2 p1 16 0 -2 p2 0 0 -4 p3 0 16 -2group fxs range z -4.0 -2.0gen zo b size 16 3 16 p0 0 0 -4 p1 16 0 -4 p2 0 0 -7 p3 0 16 -4group fznt4 range z -7.0 -4.0gen zo b size 16 2 16 p0 0 0 -7 p1 16 0 -7 p2 0 0 -9 p3 0 16 -7group fznt5 range z -9.0 -7.0gen zo b size 16 2 16 p0 0 0 -9 p1 16 0 -9 p2 0 0 -11 p3 0 16 -9group ft range z -11.0 -9.0gen zo b size 16 4 16 p0 0 0 -11 p1 16 0 -11 p2 0 0 -15 p3 0 16 -11group xs range z -15.0 -11.0;土体设置为摩尔库伦模型mo mo range group fznt2prop bulk 6.349e6 she 3.101e6 coh 14.7e3 fric 21.7 range group fznt2mo mo range group fxsprop bulk 6.818e6 shear 3.516e6 coh 0 fric 25 range group fxsmo mo range group fznt4prop bulk 6.944e6 she 3.968e6 coh 14.0e3 fric 20 range group fznt4mo mo range group fznt5prop bulk 7.639e6 she 4.365e6 coh 15.2e3 fric 23 range group fznt5mo mo range group ftprop bulk 3.472e6 she 1.984e6 coh 12.0e3 fric 24.5 range group ftmo mo range group xsprop bulk 5.556e6 she 3.175e6 coh 0.0 fric 18 range group xs;初始条件ini dens 1800 range z -8.5 0.0;边界条件fix z range z -14.9 -15.1fix x range x -0.1 0.1fix x range x 15.9 16.1fix y range y -0.1 0.1fix y range y 15.9 16.1set grav 0 0 -10.0 ini szz 0.0 grad 0 0 18000 range z -8.5 0ini szz 10200 grad 0 0 19200 range z -15.0 -8.5ini sxx 0.0 grad 0 0 6000 range z -8.5 0ini sxx 3400 grad 0 0 6400 range z -15.0 -8.5ini syy 0.0 grad 0 0 6000 range z -8.5 0 ini syy 3400 grad 0 0 6400 range z -15.0 -8.5;地下水位water density 1000water table origin 0 0 -8.5 norm 0 0 -1ini dens 1920 range z -15.0 -8.5hist id=1 unbalsolvesave cepile0.sav;在土体中心创建桩sel pile id=1 begin=(8 8 0 ) end=(8 8 -2) nseg=2sel pile prop emod=28e9 nu=0.2 xcarea=0.2826 xcj=1.27e-2 xciy=6.3585e-3 &xciz=6.3585e-3 per=1.884 cs_sk=2.097e10 cs_scoh=2.769e4 cs_sfric=13.02 &cs_nk=2.097e8 cs_ncoh=0.0 cs_nfric=0.0 cs_ngap=off ;sel pile id=1 begin=(8 8 -2) end=(8 8 -4) nseg=2sel pile prop emod=28e9 nu=0.2 xcarea=0.2826 xcj=1.27e-2 xciy=6.3585e-3 &xciz=6.3585e-3 per=1.884 cs_sk=2.301e10 cs_scoh=0.0 cs_sfric=15 &cs_nk=2.301e8 cs_ncoh=0.0 cs_nfric=0.0 cs_ngap=off ;sel pile id=1 begin=(8 8 -4) end=(8 8 -7) nseg=3sel pile prop emod=28e9 nu=0.2 xcarea=0.2826 xcj=1.27e-2 xciy=6.3585e-3 &xciz=6.3585e-3 per=1.884 cs_sk=2.447e10 cs_scoh=2.638e4 cs_sfric=12 &cs_nk=2.447e8 cs_ncoh=0.0 cs_nfric=0.0 cs_ngap=off ;sel pile id=1 begin=(8 8 -7) end=(8 8 -9) nseg=2sel pile prop emod=28e9 nu=0.2 xcarea=0.2826 xcj=1.27e-2 xciy=6.3585e-3 &xciz=6.3585e-3 per=1.884 cs_sk=2.692e10 cs_scoh=2.864e4 cs_sfric=13.8 &cs_nk=2.692e8 cs_ncoh=0.0 cs_nfric=0.0 cs_ngap=off ;sel pile id=1 begin=(8 8 -9) end=(8 8 -11) nseg=2sel pile prop emod=28e9 nu=0.2 xcarea=0.2826 xcj=1.27e-2 xciy=6.3585e-3 &xciz=6.3585e-3 per=1.884 cs_sk=1.223e10 cs_scoh=2.261e4 cs_sfric=14.7 &cs_nk=1.223e8 cs_ncoh=0.0 cs_nfric=0.0 cs_ngap=off ;hist id=21 sel pilesel force Fy end1 cid 1;顶部剪力hist id=22 sel pilesel force Fy end2 cid 1hist id=23 sel pilesel force Fy end2 cid 2hist id=24 sel pilesel force Fy end2 cid 3hist id=25 sel pilesel force Fy end2 cid 4hist id=26 sel pilesel force Fy end2 cid 5hist id=27 sel pilesel force Fy end2 cid 6hist id=28 sel pilesel force Fy end2 cid 7hist id=29 sel pilesel force Fy end2 cid 8hist id=30 sel pilesel force Fy end2 cid 9hist id=31 sel pilesel force Fy end2 cid 10hist id=32 sel pilesel force Fy end2 cid 11;hist id=33 sel pilesel moment Mz end1 cid 1;弯矩hist id=34 sel pilesel moment Mz end2 cid 1hist id=35 sel pilesel moment Mz end2 cid 2hist id=36 sel pilesel moment