您的位置：360文档中心› FLAC-3D深基坑的开挖与支护的命令流

FLAC-3D深基坑的开挖与支护的命令流

合集下载

FLAC3D基坑开挖命令流

FLAC3D基坑开挖命令流

FLAC3D基坑开挖命令流nplot block groupgen zone uwedge p0 30.18 0 -25 p1 34.59 0 -25 p2 30.18 1 -25 p3 30.18 0 0 &p4 34.59 1 -25 p5 30.18 1 0 size 5 1 25 group 1gen zone uwedge p0 34.59 0 0 p1 30.18 0 0 p2 34.59 1 0 p3 34.59 0 -25 &p4 30.18 1 0 p5 34.59 1 -25 size 5 1 25 group 2;接触⾯gen separate 2int 1 wrap 1 2int 1 maxedge 0.5gen zone brick p0 0 0 0 p1 30.18 0 0 p2 0 0 -25 p3 0 1 0 p4 30.18 0 -25 &p5 0 1 -25 p6 30.18 1 0 p7 30.18 1 -25 size 30 25 1gen zone brick p0 34.59 0 0 p1 50 0 0 p2 34.59 0 -25 p3 34.59 1 0 p4 50 0 -25 &p5 34.59 1 -25 p6 50 1 0 p7 50 1 -25 size 16 25 1group soil-1 range z 0 -6group soil-2 range z -6 -10group soil-3 range z -10 -25;模型参数model mohr range group soil-1model mohr range group soil-2model mohr range group soil-3prop bulk 4.4e10 shear 1.6e10 fric 17.2 coh 16.8e10 tens 1e10 dil 10 range group soil-1 prop bulk 5.7e10 shear 2.1e10 fric 17.7 coh 16.8e10 tens 1e10 dil 10 range group soil-2 prop bulk 10e10 shear 3.6e10 fric 20.9 coh 21.6e10 tens 1e10 dil 10 range group soil-3 interface 1 prop kn 1e8 ks 1e8 fric 15.3 coh 3.6e10 range z 0 -6interface 1 prop kn 1e8 ks 1e8 fric 16 coh 8.48e10 range z -6 -10interface 1 prop kn 1e8 ks 1e8 fric 20 coh 4.68e10 range z -10 -25ini dens 1730 range group soil-1ini dens 1860 range group soil-2ini dens 1920 range group soil-3;边界条件fix x range x -0.1 0.1fix x range x 49.9 50.1fix y range y -0.1 0.1fix y range y 0.9 1.1fix z range z -25.1 -24.9set grav 0 0 -9.8solveprop bulk 19e6 shear 8e6 fric 14.8 coh 21e3 tens 1e3 range group soil-1prop bulk 26e6 shear 10e6 fric 15.2 coh 21e3 tens 1e3 range group soil-2prop bulk 45e6 shear 18e6 fric 17.9 coh 27e3 tens 1e3 range group soil-3interface 1 prop kn 1e6 ks 1e6 fric 14 coh 4.6e3 range z 0 -6interface 1 prop kn 1e6 ks 1e6 fric 13.7 coh 10.6e3 range z -6 -10interface 1 prop kn 1e6 ks 1e6 fric 17 coh 5.9e3 range z -10 -25hist unbalsolvesave 初始平衡.savres 初始平衡.sav;pile单元sel pile id=1 begin=(39.6,0.5,0) end=(39.6,0.5,-16) nseg=16cs_sfric=18.62 cs_sk=2.871e8 cs_nk=2.871e8 & cs_ncoh=5.78e4 cs_nfric=18.62 cs_ngap=off solvesave 加桩平衡.savres 加桩平衡.savini xdisp 0 ydisp 0 zdisp 0ini xvel 0 yvel 0 zvel 0sel node ini xdis 0 ydis 0 zdis 0sel node ini xvel 0 yvel 0 zvel 0sel node ini xrdis 0 yrdis 0 zrdis 0sel node ini xrvel 0 yrvel 0 zrvel 0set largegroup 第⼀次开挖 range x 40 50 y 0 1 z 0 -3;设置监测点hist id=10 gp zdisp 31 0 0hist id=11 gp zdisp 33 0 0hist id=12 gp zdisp 35 0 0hist id=13 gp zdisp 37 0 0hist id=14 gp zdisp 39 0 0hist id=15 gp zdisp 41 0 -3hist id=16 gp zdisp 43 0 -3hist id=17 gp zdisp 45 0 -3hist id=18 gp zdisp 47 0 -3hist id=19 gp zdisp 49 0 -3hist id=20 gp xdisp 40 0 0hist id=21 gp xdisp 40 0 -1hist id=22 gp xdisp 40 0 -2hist id=23 gp xdisp 40 0 -3hist id=24 gp xdisp 40 0 -4hist id=25 gp xdisp 40 0 -5hist id=26 gp xdisp 40 0 -6hist id=27 gp xdisp 40 0 -7hist id=28 gp xdisp 40 0 -8hist id=29 gp xdisp 40 0 -9hist id=30 gp xdisp 40 0 -10hist id=31 gp xdisp 40 0 -11hist id=32 gp xdisp 40 0 -12hist id=33 gp xdisp 40 0 -13hist id=34 gp xdisp 40 0 -14hist id=35 gp xdisp 40 0 -15hist id=36 gp xdisp 40 0 -16hist id=61 sel pile force fz 39.6 0.5 -0.5hist id=62 sel pile force fz 39.6 0.5 -1.5hist id=63 sel pile force fz 39.6 0.5 -2.5hist id=64 sel pile force fz 39.6 0.5 -3.5hist id=65 sel pile force fz 39.6 0.5 -4.5hist id=66 sel pile force fz 39.6 0.5 -5.5hist id=67 sel pile force fz 39.6 0.5 -6.5hist id=68 sel pile force fz 39.6 0.5 -7.5hist id=69 sel pile force fz 39.6 0.5 -8.5hist id=70 sel pile force fz 39.6 0.5 -9.5hist id=71 sel pile force fz 39.6 0.5 -10.5hist id=72 sel pile force fz 39.6 0.5 -11.5hist id=81 sel pile mom mz 39.6 0.5 -0.5 hist id=82 sel pile mom mz 39.6 0.5 -1.5 hist id=83 sel pile mom mz 39.6 0.5 -2.5 hist id=84 sel pile mom mz 39.6 0.5 -3.5 hist id=85 sel pile mom mz 39.6 0.5 -4.5 hist id=86 sel pile mom mz 39.6 0.5 -5.5 hist id=87 sel pile mom mz 39.6 0.5 -6.5 hist id=88 sel pile mom mz 39.6 0.5 -7.5 hist id=89 sel pile mom mz 39.6 0.5 -8.5 hist id=90 sel pile mom mz 39.6 0.5 -9.5 hist id=91 sel pile mom mz 39.6 0.5 -10.5 hist id=92 sel pile mom mz 39.6 0.5 -11.5 hist id=93 sel pile mom mz 39.6 0.5 -12.5 hist id=94 sel pile mom mz 39.6 0.5 -13.5 hist id=95 sel pile mom mz 39.6 0.5 -14.5 hist id=96 sel pile mom mz 39.6 0.5 -15.5 model null range group 第⼀次开挖solvesave 第⼀次开挖.savres 第⼀次开挖.savgroup 第⼆次开挖 range x 40 50 y 0 1 z -3 -6 ;设置监测点hist id=110 gp zdisp 31 0 0hist id=111 gp zdisp 33 0 0hist id=112 gp zdisp 35 0 0hist id=113 gp zdisp 37 0 0hist id=114 gp zdisp 39 0 0hist id=115 gp zdisp 41 0 -6hist id=116 gp zdisp 43 0 -6hist id=117 gp zdisp 45 0 -6hist id=118 gp zdisp 47 0 -6hist id=119 gp zdisp 49 0 -6hist id=120 gp xdisp 40 0 0hist id=121 gp xdisp 40 0 -1hist id=122 gp xdisp 40 0 -2hist id=123 gp xdisp 40 0 -3hist id=124 gp xdisp 40 0 -4hist id=125 gp xdisp 40 0 -5hist id=126 gp xdisp 40 0 -6hist id=127 gp xdisp 40 0 -7hist id=128 gp xdisp 40 0 -8hist id=129 gp xdisp 40 0 -9hist id=130 gp xdisp 40 0 -10hist id=131 gp xdisp 40 0 -11hist id=132 gp xdisp 40 0 -12hist id=133 gp xdisp 40 0 -13hist id=134 gp xdisp 40 0 -14hist id=135 gp xdisp 40 0 -15hist id=136 gp xdisp 40 0 -16hist id=161 sel pile force fz 39.6 0.5 -0.5hihist id=166 sel pile force fz 39.6 0.5 -5.5 hist id=167 sel pile force fz 39.6 0.5 -6.5 hist id=168 sel pile force fz 39.6 0.5 -7.5 hist id=169 sel pile force fz 39.6 0.5 -8.5 hist id=170 sel pile force fz 39.6 0.5 -9.5 hist id=171 sel pile force fz 39.6 0.5 -10.5 hist id=172 sel pile force fz 39.6 0.5 -11.5 hist id=173 sel pile force fz 39.6 0.5 -12.5 hist id=174 sel pile force fz 39.6 0.5 -13.5 hist id=175 sel pile force fz 39.6 0.5 -14.5 hist id=176 sel pile force fz 39.6 0.5 -15.5 hist id=181 sel pile mom mz 39.6 0.5 -0.5 hist id=182 sel pile mom mz 39.6 0.5 -1.5 hist id=183 sel pile mom mz 39.6 0.5 -2.5 hist id=184 sel pile mom mz 39.6 0.5 -3.5 hist id=185 sel pile mom mz 39.6 0.5 -4.5 hist id=186 sel pile mom mz 39.6 0.5 -5.5 hist id=187 sel pile mom mz 39.6 0.5 -6.5 hist id=188 sel pile mom mz 39.6 0.5 -7.5 hist id=189 sel pile mom mz 39.6 0.5 -8.5 hist id=190 sel pile mom mz 39.6 0.5 -9.5 hist id=191 sel pile mom mz 39.6 0.5 -10.5 hist id=192 sel pile mom mz 39.6 0.5 -11.5 hist id=193 sel pile mom mz 39.6 0.5 -12.5 hist id=194 sel pile mom mz 39.6 0.5 -13.5 hist id=195 sel pile mom mz 39.6 0.5 -14.5 hist id=196 sel pile mom mz 39.6 0.5 -15.5 model null range group 第⼆次开挖solvesave 第⼆次开挖.savres 第⼆次开挖.savgroup 第三次开挖 range x 40 50 y 0 1 z -6 -8 ;设置监测点hist id=210 gp zdisp 31 0 0hist id=211 gp zdisp 33 0 0hist id=212 gp zdisp 35 0 0hist id=213 gp zdisp 37 0 0hist id=214 gp zdisp 39 0 0hist id=215 gp zdisp 41 0 -8hist id=216 gp zdisp 43 0 -8hist id=217 gp zdisp 45 0 -8hist id=218 gp zdisp 47 0 -8hist id=219 gp zdisp 49 0 -8hist id=220 gp xdisp 40 0 0hist id=221 gp xdisp 40 0 -1hist id=222 gp xdisp 40 0 -2hist id=223 gp xdisp 40 0 -3hist id=224 gp xdisp 40 0 -4hist id=225 gp xdisp 40 0 -5hist id=226 gp xdisp 40 0 -6hist id=227 gp xdisp 40 0 -7hist id=228 gp xdisp 40 0 -8hist id=229 gp xdisp 40 0 -9hist id=261 sel pile force fz 39.6 0.5 -0.5hist id=262 sel pile force fz 39.6 0.5 -1.5hist id=263 sel pile force fz 39.6 0.5 -2.5hist id=264 sel pile force fz 39.6 0.5 -3.5hist id=265 sel pile force fz 39.6 0.5 -4.5hist id=266 sel pile force fz 39.6 0.5 -5.5hist id=267 sel pile force fz 39.6 0.5 -6.5hist id=268 sel pile force fz 39.6 0.5 -7.5hist id=269 sel pile force fz 39.6 0.5 -8.5hist id=270 sel pile force fz 39.6 0.5 -9.5hist id=271 sel pile force fz 39.6 0.5 -10.5 hist id=272 sel pile force fz 39.6 0.5 -11.5 hist id=273 sel pile forcefz 39.6 0.5 -12.5hist id=274 sel pile force fz 39.6 0.5 -13.5hist id=275 sel pile force fz 39.6 0.5 -14.5hist id=276 sel pile force fz 39.6 0.5 -15.5 hist id=281 sel pile mom mz 39.6 0.5 -0.5 hist id=282 sel pile mom mz 39.6 0.5 -1.5 hist id=283 sel pile mom mz 39.6 0.5 -2.5 hist id=284 sel pile mom mz 39.6 0.5 -3.5 hist id=285 sel pile mom mz 39.6 0.5 -4.5 hist id=286 sel pile mom mz 39.6 0.5 -5.5 hist id=287 sel pile mom mz 39.6 0.5 -6.5 hist id=288 sel pile mom mz 39.6 0.5 -7.5 hist id=289 sel pile mom mz 39.6 0.5 -8.5 hist id=290 sel pile mom mz 39.6 0.5 -9.5 hist id=291 sel pile mom mz 39.6 0.5 -10.5 hist id=292 sel pile mom mz 39.6 0.5 -11.5 hist id=293 sel pile mom mz 39.6 0.5 -12.5 hist id=294 sel pile mom mz 39.6 0.5 -13.5 hist id=295 sel pile mom mz 39.6 0.5 -14.5 hist id=296 sel pile mom mz 39.6 0.5 -15.5 model null range group 第三次开挖solvesave 第三次开挖.savres 第三次开挖.savgroup 第四次开挖 range x 40 50 y 0 1 z -8 -10 ;设置监测点hist id=310 gp zdisp 31 0 0hist id=311 gp zdisp 33 0 0hist id=312 gp zdisp 35 0 0hist id=313 gp zdisp 37 0 0hist id=314 gp zdisp 39 0 0hist id=315 gp zdisp 41 0 -10hist id=316 gp zdisp 43 0 -10hist id=317 gp zdisp 45 0 -10hist id=318 gp zdisp 47 0 -10hist id=319 gp zdisp 49 0 -10hist id=325 gp xdisp 40 0 -5hist id=326 gp xdisp 40 0 -6hist id=327 gp xdisp 40 0 -7hist id=328 gp xdisp 40 0 -8hist id=329 gp xdisp 40 0 -9hist id=330 gp xdisp 40 0 -10hist id=331 gp xdisp 40 0 -11hist id=332 gp xdisp 40 0 -12hist id=333 gp xdisp 40 0 -13hist id=334 gp xdisp 40 0 -14hist id=335 gp xdisp 40 0 -15hist id=336 gp xdisp 40 0 -16hist id=361 sel pile force fz 39.6 0.5 -0.5 hist id=362 sel pile force fz 39.6 0.5 -1.5 hist id=363 sel pile force fz 39.6 0.5 -2.5 hist id=364 sel pile force fz 39.6 0.5 -3.5 hist id=365 sel pile force fz 39.6 0.5 -4.5 hist id=366 sel pile force fz 39.6 0.5 -5.5 hist id=367 sel pile force fz 39.6 0.5 -6.5 hist id=368 sel pile force fz 39.6 0.5 -7.5 hist id=369 sel pile force fz 39.6 0.5 -8.5 hist id=370 sel pile force fz 39.6 0.5 -9.5 hist id=371 sel pile force fz 39.6 0.5 -10.5 hist id=372 sel pile force fz 39.6 0.5 -11.5 hist id=373 sel pile force fz 39.6 0.5 -12.5 hist id=374 sel pile force fz 39.6 0.5 -13.5 hist id=375 sel pile force fz 39.6 0.5 -14.5 hist id=376 sel pile force fz 39.6 0.5 -15.5 hist id=381 sel pile mom mz 39.6 0.5 -0.5 hist id=382 sel pile mom mz 39.6 0.5 -1.5 hist id=383 sel pile mom mz 39.6 0.5 -2.5 hist id=384 sel pile mom mz 39.6 0.5 -3.5 hist id=385 sel pile mom mz 39.6 0.5 -4.5 hist id=386 sel pile mom mz 39.6 0.5 -5.5 hist id=387 sel pile mom mz 39.6 0.5 -6.5 hist id=388 sel pile mom mz 39.6 0.5 -7.5 hist id=389 sel pile mom mz 39.6 0.5 -8.5hist id=390 sel pile mom mz 39.6 0.5 -9.5 hist id=391 sel pile mom mz 39.6 0.5 -10.5 hist id=392 sel pile mom mz 39.6 0.5 -11.5hist id=393 sel pile mom mz 39.6 0.5 -12.5 hist id=394 sel pile mom mz 39.6 0.5 -13.5 hist id=395 sel pile mom mz 39.6 0.5 -14.5 hist id=396 sel pile mom mz 39.6 0.5 -15.5 model null range group 第四次开挖solvesave 第四次开挖.sav。

FLAC3D3.0在某隧道工程开挖支护中的应用

FLAC3D3.0在某隧道工程开挖支护中的应用
model null range group 仰拱初期支护
set large
hist unbal
plot hist 3
solve
save tun_ext1.sav
如果在开挖后适时对隧道进行锚喷混凝土初期支护，该命令流如下：
restore tun_nature.sav
ini xdis=0 ydis=0 zdis=0
size 4 20 6 4 dim 3 6 3 6 rat 1 1 1 1 group 围岩2
gen zon cshell p0 0 0 0 p1 0 0 -3.0 p2 0 50 0 p3 6.0 0 0 &
size 4 20 6 4 dim 2.6 5.6 2.6 5.6 rat 1 1 1 1 group 仰拱初期支护
prop bulk b_mod2 shear s_mod2 cohe 2.8e6 tens 1.0e6 fric 35 range z -40 4.5
ini dens=2300
set grav 0 0 -10
; boundary and initial conditions
apply szz -1.4e6 range z 19.9 20.1
size 4 20 6 4 dim 6 5 6 5 rat 1 1 1 1 group 围岩
gen zon cshell p0 0 0 0 p1 6.0 0 0 p2 0 50 0 p3 0 0 5.0 &
size 4 20 6 4 dim 5.6 4.6 5.6 4.6 rat 1 1 1 1 group 初期支护
FLAC3D计算缺点：由于没有一个类似autocad程序简洁可视化的建模操作程序，相比2d-σ、3d-σ、Ansys等数值模拟软件，用FLAC3D直接建立模型比较困难，一般情况下要么开发一个flac建模的前处理程序，要么从Ansys中建模，然后将节点和单元信息导出，根据flac建模的命令流格式将这些节点和单元信息改写。

《FLAC3D基础与工程实例》全部命令流

《FLAC3D基础与工程实例》全部命令流

《FLA C/FLA C 3D基础与‎工程实例》全部命令流‎1gen zone bri p0 0 0 0 p1 10 0 0 p2 0 10 0 p3 0 0 10 &p4 15 15 0 p5 0 15 15 p6 15 0 10 p7 20 20 20 &size 10 10 10 rat 1.0 0.9 1.1 group‎brick‎_1gen zone bri p0 20 0 0 p1 add 10 0 0 p2 add 0 20 0 p3 add 0 0 15 & size 10 10 10 rat 1.0 0.9 1.1 group‎brick‎_2gen zone bri p0 40 0 0 edge 10 size 10 10 10 rat 1.0 0.9 1.1 group‎brick‎_3 plot surnewgen zon bri size 3 3 3model‎elasprop bulk 3e8 shear‎1e8ini dens 2000fix z ran z -.1 .1fix x ran x -.1 .1fix x ran x 2.9 3.1fix y ran y -.1 .1fix y ran y 2.9 3.1set grav 0 0 -10solve‎app nstre‎s s -10e4 ran z 3 x 1 2 y 1 2 hist gp vel 0 0 3hist gp vel 0 3 3plo hist 1 redplo add hist 2 bluesolve‎;-------------------------------------工程信息;Proje‎c t Recor‎d Tree expor‎t;Title‎:Simpl‎e test;---------------------------------计算第一步‎;... STATE‎: STATE‎1 ....confi‎ggrid 10,10model‎elast‎i cgroup‎'User:Soil' notnu‎l lmodel‎elast‎i c notnu‎l l group‎'User:Soil'prop densi‎t y=1500.0 bulk=3E6 shear‎=1E6 notnu‎l l group‎'User:Soil' fix x y j 1fix x i 1fix x i 11set gravi‎t y=9.81histo‎r y 999 unbal‎a nced‎solve‎save state‎1.sav;----------------------------------计算第二步‎;... STATE‎: STATE‎2 ....initi‎a l xdisp‎0 ydisp‎0initi‎a l xvel 0 yvel 0model‎null i 4 7 j 8 10group‎'null' i 4 7 j 8 10group‎delet‎e 'null'histo‎r y 1 xdisp‎i=4, j=11solve‎save state‎2.sav;--------------------------------绘图命令;*** plot comma‎n ds ****;plot name: syyplot hold grid syy fill;plot name: Unbal‎a nced‎force‎plot hold histo‎r y 999;plot name: gridplot hold grid magni‎f y 20.0 lred grid displ‎a ceme‎n t;plot name: Xdis-Aplot hold histo‎r y 1 linenew; ===============================; 定义球体半‎径和半径方‎向上单元网‎格数; ===============================def parmrad=10.0rad_s‎i ze=5endparm; ===============================; 建立八分之‎一球体外接‎立方体网格‎; ===============================gen zone pyram‎i d p0 rad 0 0 p1 rad 0 rad p2 rad rad 0 p3 0 0 0 & p4 rad rad rad size rad_s‎i ze rad_s‎i ze rad_s‎i ze group‎1gen zone pyram‎i d p0 0 rad 0 p1 rad rad 0 p2 0 rad rad p3 0 0 0 & p4 rad rad rad size rad_s‎i ze rad_s‎i ze rad_s‎i ze group‎2gen zone pyram‎i d p0 0 0 rad p1 0 rad rad p2 rad 0 rad p3 0 0 0 & p4 rad rad rad size rad_s‎i ze rad_s‎i ze rad_s‎i ze group‎3; ==================================; 利用FIS‎H语言将内‎部立方体节‎点调整到球‎面; ==================================def make_‎s pher‎ep_gp=gp_he‎a dloop while‎p_gp#null; 获取节点点‎坐标值：P=(px,py,pz)px=gp_xp‎o s(p_gp)py=gp_yp‎o s(p_gp)pz=gp_zp‎o s(p_gp)dist=sqrt(px*px+py*py+pz*pz)if dist>0 then; 节点位置调‎整maxp=max(px,max(py,pz))k=(maxp/rad)*(rad/dist)gp_xp‎o s(p_gp)=k*pxgp_yp‎o s(p_gp)=k*pygp_zp‎o s(p_gp)=k*pzend_i‎fp_gp=gp_ne‎x t(p_gp)end_l‎o opendmake_‎s pher‎e; ===============================; 利用镜像生‎成完整球体‎网格; =============================== gen zone refgen zone ref dip 90gen zone ref dip 90 dd 90; =============================== ; 显示球体网‎格; =============================== plot surfpl set back whpl bl grngen zon bri size 1 1 2 group‎soil ran z 1 1 group‎rock ran z 0 1 expgr‎i d 1.flac3‎dgen zon bri size 3 3 3model‎mohrprop bu 3e6 sh 1e6 coh 10e3 fric 15 fix z ran z -.1 .1fix x ran x -.1 .1fix x ran x 2.9 3.1fix y ran y -.1 .1fix y ran y 2.9 3.1ini dens 2000hist unbal‎set grav 10solve‎elast‎i csave 6-1.savrest 6-1.savini xd 0 yd 0 zd 0 xv 0 yv 0 zv 0app nstre‎s s -100e3‎ran z 2.9 3.1 x 1 2 y 1 2 solve‎save 6-2.savrest 6-1.savini xd 0 yd 0 zd 0 xv 0 yv 0 zv 0app nstre‎s s -100e3‎ran z 2.9 3.1 x 1 2 y 1 2 hist id=2 gp zdis 1 1 3hist id=3 gp zdis 1 1 2hist id=4 gp xdis 1 1 3hist id=5 gp xdis 1 1 3hist id=6 zone szz 1 1 3hist id=7 zone szz 1.5 1.5 2.5hist id=8 zone sxz 1.5 1.5 2.5solve‎save 6-3.savrest 6-3.savset log onset logfi‎l e 6-2.log print‎zone stres‎s print‎gp disset log offrest 6-1.savini xd 0 yd 0 zd 0 xv 0 yv 0 zv 0app nstre‎s s -100e3‎ran z 2.9 3.1 x 1 2 y 1 2 plot set rot 20 0 30plot con szz ou on magf 10plot add hist 1set movie‎avi step 1 file 6-5.avimovie‎start‎solve‎movie‎finis‎hngen zon bri size 1 1 2 model‎elasprop bulk 3e7 shear‎1e7 fix z ran z 0fix x ran x 0fix x ran x 1fix y ran y 0fix y ran y 1ini dens 2000set grav 0 0 -10solve‎plo con szngen zon bri size 1 1 2model‎mohrprop bulk 3e7 shear‎1e7 c 1e10 f 15 tensi‎o n 1e10 fix z ran z 0fix x ran x 0fix x ran x 1fix y ran y 0fix y ran y 1ini dens 2000set grav 0 0 -10solve‎prop bulk 3e7 shear‎1e7 c 10e3 f 15 ten 0solve‎plo con szngen zone brick‎size 1 1 2model‎mohrprop bulk 3e7 shear‎1e7 coh 10e3 fri 15 ten 0 fix z ran z 0fix x ran x 0fix x ran x 1fix y ran y 0fix y ran y 1ini dens 2000set grav 0 0 -10solve‎elasplo con sznewgen zone brick‎size 1 1 2model‎mohrprop bulk 3e7 shear‎1e7 coh 10e3 fri 15 ten 0 fix z ran z 0fix x ran x 0fix x ran x 1fix y ran y 0fix y ran y 1ini dens 2000ini szz -40e3 grad 0 0 20e3 ran z 0 2ini syy -20e3 grad 0 0 10e3 ran z 0 2ini sxx -20e3 grad 0 0 10e3 ran z 0 2set grav 0 0 -10solve‎plo con szngen zon bri size 1 1 2model‎mprop bulk 3e7 shear‎1e7 c 10e10‎f 15 ten 1e10 fix z ran z 0fix x ran x 0fix x ran x 1fix y ran y 0fix y ran y 1ini dens 2000 ran z 0 1ini dens 1500 ran z 1 2ini szz -35e3 grad 0 0 20e3 ran z 0 1ini syy -17.5e3 grad 0 0 10e3 ran z 0 1ini sxx -17.5e3 grad 0 0 10e3 ran z 0 1ini szz -15e3 grad 0 0 15e3 ran z 1 2ini syy -7.5e3 grad 0 0 7.5e3 ran z 1 2ini sxx -7.5e3 grad 0 0 7.5e3 ran z 1 2ini pp 10e3 grad 0 0 -10e3 ran z 0 1set grav 0 0 -10solve‎plo con szngen zon bri size 1 1 2model‎mprop bulk 3e7 shear‎1e7 c 10e10‎f 15 ten 1e10 fix z ran z 0fix x ran x 0fix x ran x 1fix y ran y 0fix y ran y 1ini dens 2000 ran z 0 2ini szz -50e3 grad 0 0 20e3 ran z 0 1ini syy -30e3 grad 0 0 10e3 ran z 0 1ini sxx -30e3 grad 0 0 10e3 ran z 0 1ini pp 30e3 grad 0 0 -10e3 ran z 0 2app nstre‎s s -10e3 ran z 2set grav 0 0 -10solve‎plo con sznewgen zone brick‎p0 0 0 0 p1 60 0 0 p2 0 60 0 p3 0 0 90 & p4 60 60 0 p5 0 60 90 p6 60 0 150 p7 60 60 150 & size 6 6 10model‎elaspro bulk 10e10‎she 10e10‎ini den 2500apply‎sxx -1e9 grad 0 0 1.11111‎11e7 range‎x -.1 .1 apply‎sxx -1e9 grad 0 0 6.66666‎66e6 range‎x 59.9 60.1 apply‎syy -1e9 grad 0 0 8.33333‎33e6 range‎y -.1 .1 apply‎syy -1e9 grad 0 0 8.33333‎33e6 range‎y 59.9 60.1 apply‎szz -1e8 grad 0 0 8.33333‎33e5 ran z 0 120set grav 0 0 -10step 30000‎ini xdisp‎0 ydisp‎0 zdisp‎0ini xvel 0 yvel 0 zvel 0plo cont szznewgen zone brick‎p0 0 0 0 p1 60 0 0 p2 0 60 0 p3 0 0 90 & p4 60 60 0 p5 0 60 90 p6 60 0 150 p7 60 60 150 & size 6 6 10model‎elaspro bulk 10e10‎she 10e10‎ini den 2500ini sxx -1e9 grad 0 0 1.11111‎11e7 range‎x -.1 .1ini sxx -1e9 grad 0 0 6.66666‎66e6 range‎x 59.9 60.1ini syy -1e9 grad 0 0 8.33333‎33e6 range‎y -.1 .1ini syy -1e9 grad 0 0 8.33333‎33e6 range‎y 59.9 60.1ini szz -1e8 ran z -.1 .1fix x y z ran z -.1 .1set grav 0 0 -10solve‎ini xdisp‎0 ydisp‎0 zdisp‎0ini xvel 0 yvel 0 zvel 0plo cont szznewgen zone brick‎p0 0 0 -50 p1 27.5 0 -50 p2 0 5 -50 p3 0 0 -10 size 8 1 10 group‎claygen zone brick‎p0 27.5 0 -50 p1 100 0 -50 p2 27.5 5 -50 p3 27.5 0 -10 ratio‎1.1 1 1 size 12 1 10 group‎claygen zone brick‎p0 0 0 -10 p1 27.5 0 -10 p2 0 5 -10 p3 0 0 0 ratio‎1 1 0.8 size 8 1 4 group‎soilgen zone brick‎p0 27.5 0 -10 p1 100 0 -10 p2 27.5 5 -10 p3 27.5 0 0 ratio‎1.1 1 0.8 size 12 1 4 group‎soilgen zone brick‎p0 0 0 0 p1 27.5 0 0 p2 0 5 0 p3 0 0 5 p4 27.5 5 0 &p5 0 5 5 p6 20 0 5 p7 20 5 5 size 8 1 5 group‎damfix x y z ran z -49.9 -50.1fix x ran x -.1 .1fix x ran x 99.9 100.1fix ymodel‎mohr ran z -50 0model‎null ran z 0 5prop bulk 7.8e6 shear‎3.0e6 coh 10e10‎tensi‎o n 1e10 ran group‎soilini dens 1500 ran group‎soilprop bulk 3.91e6 shear‎1.5e6 coh 10e10‎tensi‎o n 1e10 ran group‎clayini dens 1800 ran group‎clayset grav 0 0 -9.8hist id=1 unbal‎solve‎prop bulk 7.8e6 shear‎3.0e6 coh 10e3 fric 15 ran group‎soilprop bulk 3.91e6 shear‎1.5e6 coh 20e3 fric 20 ran group‎claysolve‎save elast‎i c.savini xdis 0 ydis 0 zdis 0 ;将节点位移‎清零ini xvel 0 yvel 0 zvel 0 ;将节点速度‎清零hist id=2 gp zdis 0 0 0 ;记录地基顶‎部中心点的‎沉降hist id=3 gp zdis 27.5 0 0 ;记录路基坡‎脚处的沉降‎hist id=4 gp xdis 27.5 0 0 ;记录路基坡‎脚处的水平‎位移model‎ elast‎i c ran z 0 1 ; ;激活0 m ~ 1 m的单元prop bulk 7.8e6 shear‎3.0e6 ran z 0 1ini dens 1500 ran z 0 1solve‎;按软件默认‎精度求解save fill-1.savmodel‎elast‎i c ran z 1 2prop bulk 7.8e6 shear‎3.0e6 ran z 1 2ini dens 1500 ran z 1 2solve‎save fill-2.savmodel‎elast‎i c ran z 2 3prop bulk 7.8e6 shear‎3.0e6 ran z 2 3ini dens 1500 ran z 2 3solve‎save fill-3.savmodel‎elast‎i c ran z 3 4prop bulk 7.8e6 shear‎3.0e6 ran z 3 4ini dens 1500 ran z 3 4solve‎save fill-4.savmodel‎elast‎i c ran z 4 5prop bulk 7.8e6 shear‎3.0e6 ran z 4 5ini dens 1500 ran z 4 5solve‎save fill-5.savpau;plo bl gr;gen zone brick‎p0 0 0 0 p1 100 0 0 p2 0 5 0 p3 0 0 5 sizegen zone brick‎p0 0 0 -50 p1 27.5 0 -50 p2 0 5 -50 p3 0 0 -10 size 8 1 10 group‎clay gen zone brick‎p0 27.5 0 -50 p1 100 0 -50 p2 27.5 5 -50 p3 27.5 0 -10 ratio‎1.1 1 1 size 12 1 10 group‎claygen zone brick‎p0 0 0 -10 p1 27.5 0 -10 p2 0 5 -10 p3 0 0 0 ratio‎1 1 0.8 size 8 1 4 group‎soilgen zone brick‎p0 27.5 0 -10 p1 100 0 -10 p2 27.5 5 -10 p3 27.5 0 0 ratio‎1.1 1 0.8 size 12 1 4 group‎soilgen zone brick‎p0 0 0 0 p1 27.5 0 0 p2 0 5 0 p3 0 0 5 p4 27.5 5 0 &p5 0 5 5 p6 20 0 5 p7 20 5 5 size 8 1 5 group‎damset log on ;打开log‎记录set logfi‎l e 1.log ;设置记录文‎件名为：1.logresto‎re fill-1.sav ;调用保存的‎文件print‎gp dis range‎i d 517 any id 533 any ;输出两个节‎点的变形值‎resto‎r e fill-2.savprint‎gp dis range‎id 517 any id 533 anyresto‎r e fill-3.savprint‎gp dis range‎id 517 any id 533 anyresto‎r e fill-4.savprint‎gp dis range‎id 517 any id 533 anyresto‎r e fill-5.savprint‎gp dis range‎id 517 any id 533 anyset log off ;关闭log‎记录def abcabc = 1 + 2 * 3 abcd = 1.0 / 2.0 endabcprint‎fishdef abcif aa < 0 thenabc = 0.0 elseabc = 2.0 * aa endif‎endabcdef abcloop aa (1, 2.5) comma‎n dprint‎aaendco‎m mand‎ endlo‎o pendnewgen zon bri size 3 3 3model‎elast‎i cprop bu 3e7 sh 1e7ini dens 2000fix x y z ran z -.1 .1fix x ran x -.1 .1fix x ran x 2.9 3.1fix y ran y -.1 .1fix y ran y 2.9 3.1set grav 10solve‎ini xd 0 yd 0 zd 0 xv 0 yv 0 zv 0save 8-4.sav8-5rest 8-4.savdef E_mod‎i fyp_z = zone_‎h eadd_k = 704d_n = 0.38d_pa = 10132‎5.0 ;//标准大气压‎loop while‎p_z # nullsigma‎_3 = -1.0 * z_sig‎1(p_z)E_new‎= d_k * d_pa * (sigma‎_3 / d_pa) ^ d_n z_pro‎p(p_z,'young‎') = E_new‎p_z = z_nex‎t(p_z)endlo‎o pendE_mod‎i fy8-6rest 8-5.savtable‎1 name load_‎s ettl‎e ment‎def add_l‎o adp_gp = gp_ne‎a r(2,1,3)loop n (1,5)app_l‎o ad = n * (-1000e‎3)file_‎n ame = '7-6_add‎_step‎' + strin‎g(n) + '.sav' comma‎n dapp nstre‎s s app_l‎o ad ran z 2.9 3.1 x 1 2 y 1 2 s olve‎s ave file_‎n ameendco‎m mand‎xtabl‎e(1,n) = -1.0 * app_l‎o adytabl‎e(1,n) = gp_zd‎i sp(p_gp)endlo‎o pendadd_l‎o adsave 8-6.sav8-7rest 8-6.savdef find_‎m ax_d‎i spp_gp = gp_he‎a dmaxdi‎s p_va‎l ue = 0.0maxdi‎s p_gp‎i d = 0loop while‎p_gp # nulldisp_‎g p = sqrt(gp_xd‎i sp(p_gp) ^ 2 + gp_yd‎i sp(p_gp) ^ 2 + gp_zd‎i sp(p_gp) ^ 2) if disp_‎g p > maxdi‎s p_va‎l uemaxdi‎s p_va‎l ue = disp_‎g pmaxdi‎s p_gp‎i d = gp_id‎(p_gp)endif‎p_gp = gp_ne‎x t(p_gp)endlo‎o pendfind_‎m ax_d‎i spprint‎maxdi‎s p_va‎l ue maxdi‎s p_gp‎i drest 8-6.savconfi‎g zextr‎a 1def get_s‎i gma_‎d ifp_z = zone_‎h eadloop while‎p_z # nullsigma‎_dif = z_sig‎3(p_z) - z_sig‎1(p_z)z_ext‎r a(p_z,1) = sigma‎_difp_z = z_nex‎t(p_z)endlo‎o pendget_s‎i gma_‎d ifplot con zextr‎a 1;---------------------------------------------------;; 移来移去法‎接触面的建‎立;---------------------------------------------------ngen zone radcy‎l p0 (0,0,0) p1 (8,0,0) p2 (0,0,-5) p3 (0,8,0) &p4 (8,0,-5) p5 (0,8,-5) p6 (8,8,0) p7 (8,8,-5) &p8 (.3,0,0) p9 (0,.3,0) p10 (.3,0,-5) p11 (0,.3,-5) & size 3 10 6 15 ratio‎1 1 1 1.15gen zone radcy‎l p0 (0,0,-5) p1 (8,0,-5) p2 (0,0,-8) p3 (0,8,-5) &p4 (8,0,-8) p5 (0,8,-8) p6 (8,8,-5) p7 (8,8,-8) &p8 (.3,0,-5) p9 (0,.3,-5) p10 (.3,0,-8) p11 (0,.3,-8) & size 3 6 6 15 ratio‎1 1 1 1.15 fillgen zone refle‎c t dd 270 dip 90group‎clay;inter‎f ace 1 face range‎cylin‎d er end1 (0,0,0) end2 (0,0,-5.1) radiu‎s .31 &cylin‎d er end1 (0,0,0) end2 (0,0,-5.1) radiu‎s .29 notinter‎f ace 2 face range‎cylin‎d er end1 (0,0,-4.9) end2 (0,0,-5.1) radiu‎s .31;pause‎gen zone cyl p0 (0,0,6) p1 (.3,0,6) p2 (0,0,1) p3 (0,.3,6) &p4 (.3,0,1) p5 (0,.3,1) &size 3 10 6gen zone cyl p0 (0,0,6.1) p1 (.3,0,6.1) p2 (0,0,6) p3 (0,.3,6.1) &p4 (.3,0,6) p5 (0,.3,6) &size 3 1 6gen zone refle‎c t dd 270 dip 90 range‎z 1 6.1group‎pile range‎z 1 6.1pause‎ini z add -6.0 range‎group‎pilesave pile_‎g eom.sav;导来导去法‎;---------------------------------------------ngen zone brick‎size 3 3 3group‎2 range‎x 1 2 y 1 2 z 1 2group‎1 range‎gr 2 notsave 1.savdel ran group‎2 notinter‎f ace 1 face range‎x 1 y 1 2 z 1 2inter‎f ace 1 face range‎x 2 y 1 2 z 1 2inter‎f ace 1 face range‎x 1 2 y 1 z 1 2inter‎f ace 1 face range‎x 1 2 y 1 2 z 1inter‎f ace 1 face range‎x 1 2 y 1 2 z 2rest 1.savdel ran group‎2expgr‎i d 1.fac3d‎impgr‎i d 1.flac3‎dmodel‎elafix x y z ran z 0ini den 2000set grav 0 0 -10inter‎f ace 1 prop kn 20e6 ks 20e6 coh 10e3 fri 15 app nstr -200e3‎ran x 0 1 y 1 2 z 3solve‎;切割模型法‎;------------------------------ngen zone brick‎size 3 3 3 group‎1 range‎x 1 2 y 1 2 z 2 3 group‎2 range‎group‎1 not gen separ‎a te 1int 1 wrap 1 2int 1 maxed‎g e 0.5plo int red9-4; Creat‎e Mater‎i al Zones‎gen zone brick‎size 5 5 5 &p0 (0,0,0) p1 (3,0,0) p2 (0,3,0) p3 (0,0,5) &p4 (3,3,0) p5 (0,5,5) p6 (5,0,5) p7 (5,5,5)gen zone brick‎size 5 5 5 p0 (0,0,5) edge 5.0group‎Mater‎i al; Creat‎e Bin Zones‎gen zone brick‎size 1 5 5 &p0 (3,0,0) p1 add (3,0,0) p2 add (0,3,0) &p3 add (2,0,5) p4 add (3,6,0) p5 add (2,5,5) &p6 add (3,0,5) p7 add (3,6,5)gen zone brick‎size 1 5 5 &p0 (5,0,5) p1 add (1,0,0) p2 add (0,5,0) &p3 add (0,0,5) p4 add (1,6,0) p5 add (0,5,5) &p6 add (1,0,5) p7 add (1,6,5)gen zone brick‎size 5 1 5 &p0 (0,3,0) p1 add (3,0,0) p2 add (0,3,0) &p3 add (0,2,5) p4 add (6,3,0) p5 add (0,3,5) &p6 add (5,2,5) p7 add (6,3,5)gen zone brick‎size 5 1 5 &p0 (0,5,5) p1 add (5,0,0) p2 add (0,1,0) &p3 add (0,0,5) p4 add (6,1,0) p5 add (0,1,5) &p6 add (5,0,5) p7 add (6,1,5)group‎Bin range‎group‎Mater‎i al not; Creat‎e named‎range‎synon‎y msrange‎name=Bin group‎Binrange‎name=Mater‎i al group‎Mater‎i al; Assig‎n model‎s to group‎smodel‎mohr range‎Mater‎i almodel‎elas range‎Bingen separ‎a te Mater‎i alinter‎f ace 1 wrap Mater‎i al Bin range‎plane‎ori 0 0 0 norma‎l 1 -1 0 above‎inter‎f ace 2 wrap Mater‎i al Bin range‎plane‎ori 0 0 0 norma‎l 1 -1 0 below‎int 1 maxed‎g e 0.55int 2 maxed‎g e 0.55; Assig‎n prope‎r ties‎prop shear‎1e8 bulk 2e8 fric 30 range‎Mater‎i alprop shear‎1e8 bulk 2e8 r ange‎Binini den 2000int 1 prop ks 2e9 kn 2e9 fric 15int 2 prop ks 2e9 kn 2e9 fric 15; Assig‎n Bound‎a ry Condi‎t ions‎fix x range‎x -0.1 0.1 any x 5.9 6.1 any fix y range‎y -0.1 0.1 any y 5.9 6.1 any fix z range‎z -0.1 0.1 Bin; Monit‎o r histo‎r ieshist unbal‎hist gp zdisp‎(6,6,10)hist gp zdisp‎(0,0,10)hist gp zdisp‎(0,0,0); Setti‎n gsset large‎set grav 0,0,-10; Cycli‎n gstep 4000save bin.sav9-5;;------------------------------------------------------------------;rest pile_‎g eom.sav‎group‎claymodel‎mohr rangeprop bulk 8.333e7‎ shear‎3.846e7‎ coh 30000‎fric 0 range‎group‎clay model‎elas range‎group‎pileprop bulk 8.333e7‎ shear‎3.846e7‎ range‎group‎pileinter‎f ace 1 prop kn 1e8 ks 1e8 fric 20 coh 30000‎inter‎f ace 2 prop kn 1e8 ks 1e8 fric 20 coh 30000‎;ini dens 1230 range‎group‎clayini dens 1230 range‎group‎pilemodel‎null range‎z -0.1 0.15;fix z range‎z -8.1 -7.9fix x range‎x -8.1 -7.9fix x range‎x 7.9 8.1fix y range‎y -.1 .1fix y range‎y 7.9 8.1set grav 0 0 -10ini szz 0. grad 0 0 12300‎. range‎z -5.5 0.ini szz 17600‎grad 0 0 15500‎range‎z -8 -5.5ini sxx 0. grad 0 0 5271.4 range‎z -5.5 0.ini sxx 7542.86 grad 0 0 6642.86 range‎z -8 -5.5ini sxx add 31428‎.6 grad 0 0 5714.3 range‎z -8 -5.5ini syy 0. grad 0 0 5271.4 range‎z -5.5 0.ini syy 7542.86 grad 0 0 6642.86 range‎z -8 -5.5ini syy add 31428‎.6 grad 0 0 5714.3 range‎z -8 -5.5;water‎densi‎t y 1000water‎table‎origi‎n 0,0,-5.5 norma‎l 0 0 -1ini dens 1550 range‎z -8 -5.5hist unbal‎;solve‎rat 1.e-6save pile0‎.sav;model‎elas range‎group‎pileprop bulk 13.9e9 shear‎10.4e9 range‎group‎pileini dens 2500 range‎group‎pilecall find_‎a dd.fissolve‎rat 1.e-6save pile1‎.savrest pile1‎.sav ;调用保存文‎件ini state‎0ini xdis 0.0 ydis 0.0 zdis 0.0 ;位移清零apply‎szz -0.4e6 range‎z 0.05 0.15 group‎pile ;桩顶加第一‎级荷载solve‎save app0.4.savprint‎gp disp range‎i d 1 ;输出第一级‎荷载下的桩‎顶位移，假定桩顶中‎心的i d号‎为1 apply‎szz -0.6e6 range‎z 0.05 0.15 group‎pile ;桩顶加第二‎级荷载solve‎save app0.6.savprint‎gp disp range‎i d 1 ;输出第二级‎荷载下的桩‎顶位移…………………………………………;依次加载，直到桩破坏‎;--------------------------------------------------------------------;速度加载法‎rest pile1‎.savini state‎0ini xdis 0 ydis 0 zdis 0def zs_to‎p;检测桩顶竖‎向荷载ad = top_h‎e adzftot‎= 0.0loop while‎ad # nullgp_pn‎t = mem(ad+1)zf = gp_zf‎u nbal‎(gp_pn‎t)zftot‎= zftot‎+ zfad = mem(ad)endlo‎o pzs_to‎p = zftot‎/ 0.1414endfix z range‎z 0.05 .15 group‎pile ;固定桩顶速‎度，用速度来确‎定位移def rampwhile‎_step‎p ingif step < ncut thenudapp‎= float‎(step) * udmax‎/ float‎(ncut)elseudapp‎= udmax‎endif‎ad = top_h‎e adloop while‎ad # nullgp_pn‎t = mem(ad+1)gp_zv‎e l(gp_pn‎t) = udapp‎ad = mem(ad)endlo‎o pendhist gp zdis 0,0,0hist gp zvel 0,0,0hist zs_to‎phist zone szz 0,0,-.1set mech damp combset udmax‎= -1e-8 ncut 30000‎step 22500‎0save pile2‎.sav;;------------------------------------------------------------------------;位移控制法‎def solve‎_step‎sloop n (1,21)save_‎f ile = strin‎g(n) + '-step.sav'comma‎n ds tep 40000‎s ave save_‎f ilepri zone stres‎s ran id 2381 a id 2361 a id 2341 a ;输出桩顶网‎格单元的应‎力 endco‎m mand‎endlo‎o pendsolve‎_step‎s;-----------------------------------------------------------------------------sel beam beg 0 0 0 end 2 0 0 nseg 2sel beam beg 2 0 0 end 4 0 -1 nseg 3sel beam id=2 beg 4 0 -1 end 5 0 -2 nseg 2 plot sel geo id on nod on scale‎0.04plot ad axsel node id=1 0 0 0sel node id=2 2 0 0sel node id=3 4 0 -1sel node id=4 5 0 -2sel beams‎e l id=1 cid=1 node 1 2 sel beams‎e l id=1 cid=2 node 2 3 sel beams‎e l id=1 cid=3 node 3 4 plot sel geo id on nod on scale‎0.04 plot ad axdef set_v‎a lsptA = 25.0 * sin( 40.0*degra‎d )ptB = 25.0 * cos( 40.0*degra‎d )endset_v‎a lsgen zone cylin‎d er p0=( 0.0, 0.0, 0.0 ) &p1=( ptA, 0.0, ptB ) &p2=( 0.0, 25.0, 0.0 ) &p3=( 0.0, 0.0, 25.0 ) &p4=( ptA, 25.0, ptB ) &p5=( 0.0, 25.0, 25.0 ) &size=(1, 2, 2)sel shell‎id=5 range‎cylin‎d er end1=(0.0, 0.0,0.0) &e nd2=(0.0,25.0,0.0) radiu‎s=24.5 not plot blo groplot ad sel geom black‎black‎cid on scale‎=0.03plot ad axpaudelet‎e ; delet‎e all zones‎sel node init zpos add -25.0newgen zone brick‎size 6 8 8model‎mohrprop bulk 1e8 shear‎0.3e8 fric 35prop coh 1e10 tens 1e10set grav 0 0 -9.81ini dens 1000fix x range‎x -0.1 0.1fix x range‎x 5.9 6.1fix y range‎y -0.1 0.1fix y range‎y 7.9 8.1fix z range‎z -0.1 0.1hist n 5hist unbal‎set mech force‎50solve‎save beam-brace‎0.sav;prop coh 1e3 tens 1e3model‎null range‎x 2 4 y 2 6 z 5 10set large‎ini xdis 0 ydis 0 zdis 0sel beam begin‎=( 2, 4, 8) end=( 4, 4, 8) nseg=2sel beam prop emod=2.0e11 nu=0.30sel beam prop XCAre‎a=6e-3 XCIz=200e-6 XCIy=200e-6 XCJ=0.0 hist gp zdisp‎4 4 8solve‎save beam-brace‎1.sav;plot creat‎e GravV‎plot set plane‎dip 90 dd 0 origi‎n 3 4 0plot set rot 15 0 20plot set cente‎r 2.5 4.2 4.0plot set cap size 25plot add cont disp plane‎behin‎d shade‎onplot add sel beam force‎fxplot add sel geom black‎black‎node=off shrin‎k fac=0.0plot add axesplot show;非全长锚固‎、预紧力锚杆‎（锚索）模拟;方法1、通过删除－建立lin‎k链接来模‎拟托盘gen zone radtu‎n p0 0,0,0 p1 25,0,0 p2 0,50, 0 p3 0,0,25 size 4 25 4 10 dim 4 4 4 4 ratio‎1 1 1 1.1 fillgen zone refle‎c t norma‎l 1 0 0 ori 0 0 0gen zone refle‎c t norma‎l0 0 1 ori 0 0 0mo mohrpro bulk 2.2e9 she 1.3e9 fric 30 coh 1.3e6 ten 1.5e5ini dens 2000fix x range‎x -25.1,-24.9fix x range‎x 24.9 25.1fix y range‎y 49.9 50.1fix z range‎z -25.1 -24.9fix z range‎z 24.9 25.1sel cable‎id=1 beg 0, 0, 0 end 0 ,29, 0 nseg 10sel cable‎id=1 beg 0,29,0 end 0,35,0 nseg 6sel cable‎id=1 prop emod 2e10 ytens‎i on 310e3‎xcare‎a 0.00049‎06 &gr_co‎h 1 gr_k 1 gr_pe‎r 0.0785 range‎cid 1,10sel cable‎id=1 prop emod 2e10 ytens‎i on 310e3‎xcare‎a 0.00049‎06 &gr_co‎h 10e5 gr_k 2e7 range‎cid 11,17sel delet‎e link range‎id 1sel link id=100 1 targe‎t zonesel link attac‎h xdir=rigid‎ydir=rigid‎zdir=rigid‎xrdir‎=rigid‎yrdir‎=rigid‎zrdir‎=rigid‎range‎id 100sel cable‎id=1 prete‎n sion‎60e3 range‎cid 1,10step 2000sav 10-5.sav;非全长锚固‎、预紧力锚杆‎（锚索）模拟;方法2、通过设置极‎大锚固剂参‎数模拟托盘‎gen zone radtu‎n p0 0,0,0 p1 25,0,0 p2 0,50, 0 p3 0,0,25 size 4 25 4 10 dim 4 4 4 4 ratio‎1 1 1 1.1 fillgen zone refle‎c t norma‎l 1 0 0 ori 0 0 0gen zone refle‎c t norma‎l0 0 1 ori 0 0 0mo mohrpro bulk 2.2e9 she 1.3e9 fric 30 coh 1.3e6 ten 1.5e5ini dens 2000fix x range‎x -25.1,-24.9fix x range‎x 24.9 25.1fix y range‎y 49.9 50.1fix z range‎z -25.1 -24.9fix z range‎z 24.9 25.1sel cable‎id=1 beg 0, 0, 0 end 0 ,29, 0 nseg 10sel cable‎id=1 beg 0,29,0 end 0,35,0 nseg 6sel cable‎ prop emod 2e10 ytens‎i on 310e3‎xcare‎a 0.00049‎06 &gr_co‎h 1 gr_k 1 gr_pe‎r 0.0785 range‎cid 2,10sel cable‎ prop emod 2e10 ytens‎i on 310e3‎xcare‎a 0.00049‎06 &gr_co‎h 10e5 gr_k 2e7 range‎cid 11,17sel cable‎ prop emod 2e10 ytens‎i on 310e3‎xcare‎a 0.00049‎06 &gr_co‎h 10e8 gr_k 2e10 range‎cid 1,1sel cable‎id=1 prete‎n sion‎60e3 range‎cid 1,10step 2000sav 10-6.sav;非全长锚固‎、预紧力锚杆‎（锚索）模拟;方法3、借助别的结‎构单元（如line‎r单元）来模拟托盘‎gen zone radtu‎n p0 0,0,0 p1 25,0,0 p2 0,50, 0 p3 0,0,25 size 4 25 4 10 dim 4 4 4 4 ratio‎1 1 1 1.1 fillgen zone refle‎c t norma‎l 1 0 0 ori 0 0 0gen zone refle‎c t norma‎l0 0 1 ori 0 0 0mo mohrpro bulk 2.2e9 she 1.3e9 fric 30 coh 1.3e6 ten 1.5e5ini dens 2000fix x range‎x -25.1,-24.9fix x range‎x 24.9 25.1fix y range‎y 49.9 50.1fix z range‎z -25.1 -24.9fix z range‎z 24.9 25.1sel cable‎id=1 beg 0, 0, 0 end 0 ,29, 0 nseg 10sel cable‎id=1 beg 0,29,0 end 0,35,0 nseg 6sel cable‎id=1 prop emod 2e10 ytens‎i on 310e3‎xcare‎a 0.00049‎06 &gr_co‎h 1 gr_k 1 gr_pe‎r 0.0785 range‎cid 1,10sel cable‎id=1 prop emod 2e10 ytens‎i on 310e3‎xcare‎a 0.00049‎06 &gr_co‎h 10e5 gr_k 2e7 range‎cid 11,17sel liner‎range‎y=-.1, .1 x=-1,1 z=-1,1sel liner‎PROP iso=( 25e9, 0.15) thick‎=0.1 ; concr‎e tesel liner‎PROP cs_nk‎=8e8 cs_sk‎=8e8 &cs_nc‎u t=0.0 cs_sc‎o h=0.0 cs_sc‎o hres‎=0.0 cs_sf‎r ic=0.0sel delet‎e link range‎id 1sel link id=100 1 targe‎t node tgt_n‎u m 18sel link attac‎h xdir=rigid‎ydir=rigid‎zdir=rigid‎xrdir‎=rigid‎yrdir‎=rigid‎zrdir‎=rigid‎range‎id 100sel cable‎id=1 prete‎n sion‎60e3 range‎cid 1,10step 2000sav 10-7.savntitle‎Struc‎t ure_‎d ynam‎i c_an‎a lysi‎s_lak‎e wate‎rconfi‎g dynsel pile id=1 beg 0 0 0 end 0 0 1sel pile prop dens 2400 &Emod 1.0e10 Nu 0.3 XCAre‎a 0.3 &XCJ 0.16375‎XCIy 0.00625‎XCIz 0.01575‎&Per 2.8 &CS_sK‎1.3e11 CS_sC‎o h 0.0 CS_sF‎r ic 10.0 &CS_nK‎1.3e11 CS_nC‎o h 0.0 CS_nF‎r ic 0.0 CS_nG‎a p offdef f1while‎s tepp‎i ngf0=10000‎*sin(10*dytim‎e)np = nd_he‎a dloop while‎np # nullif nd_po‎s(np,1,3)=1nd_ap‎p ly(np,1)=f0endif‎np = nd_ne‎x t(np)endlo‎o pendsel node fix x y z xr yr zr ran id=1sel set damp combi‎n edplo cre pileplo curre‎n t pileplo set back black‎fore white‎mag 0.8plo add sel geo id on sca .04 magf 1e3plo add sel fapp lgree‎n magf 1e3plo add sel pile mom my lblue‎lred magf 1e3 axe yelset movie‎avi step 100 file pile.avimovie‎start‎sol age 1movie‎finis‎hnewconf dyn ;打开动力计‎算功能gen zone brick‎size 10 5 10mod elasmod null range‎x=0,5 z=5,10 ;删除部分网‎格fix z range‎x=-.1 .1 z=.1 10.1 ;设置静力边‎界条件fix z range‎x=9.9,10.1 z=.1 10.1fix y range‎y=-.1 .1fix y range‎y=4.9 5.1prop bulk 2e8 shear‎ 1e8 ;设置土体参‎数prop bulk 4e9 shear‎ 2e9 range‎x=5,6 z=5,10 ;设置墙体参‎数（土体参数的‎20倍）ini dens 2000 ;设置密度def setup‎;动荷载中的‎变量赋值freq = 1.0omega‎= 2.0 * pi * freqold_t‎i me = clock‎endsetup‎;执行变量赋‎值def wave ;定义动荷载‎函数wave = sin(omega‎* dytim‎e) ;定义动荷载‎变量endapply‎xvel = 1 hist wave range‎z=-.1 .1 ;施加动荷载‎apply‎zvel = 0 r ange‎z=-.1 .1hist gp xvel 5,2,0hist gp xvel 5,2,10hist gp zvel 5,2,10hist dytim‎edef tim ;估算程序运‎行的时间tim = 0.01 * (clock‎- old_t‎i me)endset dyn multi‎ on ;设置动态多‎步solve‎age 1.0print‎ tim ;输出计算时‎间print‎dyn ;输出动力计‎算相关信息‎save mult1‎.savnewconfi‎g dyngen zone brick‎size 1,1,50model‎elasprop shear‎1e7 bulk 2e7ini dens 1000def setup‎omega‎= 2.0 * pi * freqpulse‎= 1.0 / freqendset freq=4.0setup‎def waveif dytim‎e > pulse‎wave = 0.0elsewave = 0.5 * (1.0 - cos(omega‎* dytim‎e))endif‎endrange‎name botto‎m z=-.1 .1fix z range‎z=.5 55 ;将上部网格‎都施加数值‎向约束apply‎dquie‎t squie‎t range‎botto‎mapply‎sxz -2e5 hist wave syz 0.0 szz 0.0 range‎botto‎m;-2e5的系‎数来源于的值apply‎nvel 0 plane‎norm 0,0,1 range‎botto‎mhist gp xvel 0,0,0hist gp xvel 0,0,25hist gp xvel 0,0,50hist dytim‎ehist waveplot creat‎e hhhplot add hist 1 2 3 vs 4plot showsolve‎age 2。

FLAC模拟隧道开挖支护的实例

FLAC模拟隧道开挖支护的实例
b_mod2=E_mod2/(3.0*(1.0-2.0*p_ratio2))
end
set E_mod1=0.6e9 p_ratio1=0.27 E_mod2=0.8e9 p_ratio2=0.26
derive
prop bulk b_mod1 shear s_mod1 cohe 1.8e6 tens 0.8e6 fric 30 range z 4.5 20
; mohr-coulomb model
model mohr
def derive
s_mod1=E_mod1/(2.0*(1.0+p_ratio1))
b_mod1=E_mod1/(3.0*(1.0-2.0*p_ratio1))
s_mod2=E_mod2/(2.0*(1.0+p_ratio2))
gen zon cshell p0 0 0 0 p1 6.0 0 0 p2 0 50 0 p3 0 0 5.0 &
size 4 20 6 4 dim 5.6 4.6 5.6 4.6 rat 1 1 1 1 group 初期支护
gen zon cshell p0 0 0 0 p1 5.6 0 0 p2 0 50 0 p3 0 0 4.6 &
;FLAC3D3.0在某隧道工程开挖支护中的应用
;隧道建模命令流入下：
new
set log on
set logfile yang.log
gen zon radcyl p0 0 0 0 pamp;
size 4 20 6 4 dim 6 5 6 5 rat 1 1 1 1 group 围岩
gen zon cshell p0 0 0 0 p1 0 0 -3.0 p2 0 50 0 p3 6.0 0 0 &

FLAC-3D深基坑的开挖与支护的命令流

FLAC-3D深基坑的开挖与支护的命令流

FLAC D3深基坑的开挖与支护的命令流一、实例工程南宁地区地层属于河流阶地二元地层，广泛分布有较厚的圆砾层，国内尚无在类似地层条件下建设地铁基坑的经验，为此，可使用FLAC3D 对基坑开挖的全过程进行三维数值模拟，在对比实测数据的基础上，总结圆砾层中地铁车站深基坑的地下连续墙水平变形及周围地表沉降变形特征。

该基坑位于大学路与明秀路交叉路口处，沿大学东路东西向布置。

车站基坑长465m，标准断面宽度为20.7m，为地下两层式结构，底板埋深为15.535m（相对地面），顶板覆土厚度大于3m。

本工程主体建筑面积21163.6m2，主要结构形式为双柱三跨框架箱型结构。

本工程所处的大学路为南宁市东西向的主要交通枢纽，车流量大，人流密集，地面条件复杂。

基坑施工采用明挖顺作法施工，围护结构为800mm厚地下连续墙+内撑（三道内支撑加一道换撑）的支护体系。

第一道支撑采用钢筋混凝土支撑，尺寸为800×900mm，冠梁同时作为第一道钢筋混凝土支撑的围檩。

第二、三道支撑及换撑使用钢支撑并施加预加力，直径为609mm，壁厚为t=16mm，斜撑段采用800×1000mm钢筋砼腰梁，其余为2×I45C 钢围檩。

二、模型建立建模工作由两部分组成，实体模型部分，包括土体和地下连续墙；结构单元部分，包括混凝土支撑和钢支撑。

根据对称性原理，拟选取1/2 的实际工程尺寸进行分析。

考虑到实际的基坑长度将近500m，根据以往的经验，选取全部长度的一半虽然能够得到满意的结果，但是由于中间部分的基坑基本处于同样的受力状态，这样会使大部分的计算长度变为重复的计算，降低了计算效率。

根据初步计算结果和经验，最终确定的基坑尺寸为，宽度取基坑的最大宽度24m，开挖深度19m，基坑长度36m。

根据地勘报告，合并相似土层，模型中共划分了7个土层。

在FLAC3D 中，围护结构可以用衬砌单元(liner)或实体单元模拟。

根据Zdravdovi的研究，在二维平面基坑模拟中，分别采用实体单元和梁单元(相当于三维模型中的衬砌单元)计算所产生的墙体变形差别小于4%，而引起地表沉降的主要原因是围护结构变形造成的地层损失，可见上述两种方法计算结果的差别可忽略不计。

FLAC3D模拟实例循环开挖与支护

FLAC3D模拟实例循环开挖与支护

FLAC3D模拟实例循环开挖与支护nres ini.savset geometry=0.001ini ydis0ini xdis0ini zdis0ini yvel0ini xvel0ini zvel0m mprop bulk 4.0e9shear 2.5e9fri32coh 2.0e6& range grou diban-shayan;prop bulk 1.8e9shear 1.2e9fri25coh 1.0e6& range grou diban-niyan any grou hangdao any;prop bulk 1.2e9shear0.8e9fri22coh0.8e6& range grou diban-gentuyan;prop bulk 1.9e9shear 1.3e9fri24coh 1.0e6& range grou diban-tniyan;prop bulk0.7e9shear0.8e9fri21coh0.7e6& range grou mc any grou gzm any;prop bulk 3.0e9shear 2.5e9fri30coh 1.8e6& range grou dingban-fenshayan;prop bulk 1.5e9shear 1.2e9fri25coh 1.1e6& range grou dingban-niyan;prop bulk 3.5e9shear 2.5e9fri34coh 1.4e6& range grou dingban-shayan;添加接触面gen separate gzminterface1wrap mc gzminterface2wrap dingban-fenshayan gzm interface1prop kn20e9ks10e9tens1e9 interface2prop kn20e9ks10e9tens1e9set mech ratio=5e-4def excav_mcloop n(excav_p,excav_p_z+cut_liang);每次开挖量cut_0=excav_pcut_1=excav_p+cut_liang;开挖commandm null range grou gzm z cut_0cut_1step100end_commandn=excav_p+cut_Liangexcav_p=excav_p+cut_Liang;条件判断保存文件,这里判断条件必须和cut_liang对应上，否则不能得到想要的文件。

基于flac3D深基坑开挖模拟与支护设计

基于flac3D深基坑开挖模拟与支护设计

本科生毕业论文（设计）题目：基于flac3D深基坑开挖模拟与支护设计指导教师: 职称:评阅人: 职称:摘要随着城市化过程中不断涌现的高层建筑和超高层建筑以及城市地下空间的开发，深基坑工程越来越多，深基坑工程项目的规模和复杂性日益增大，给深基坑工程的设计和施工带来了更大的挑战。

在这样的背景下，深基坑支护结构设计和变形量预测已成为岩土工程领域的重要研究课题之一。

本文以武汉市万达广场深基坑工程作为研究对象，利用勘查资料和深基坑支护结构设计要求，比选合理的基坑支护方案并进行相应的计算设计。

同时，本文针对深基坑工程变形量验算等难以解决的问题引用了flac3D数值模拟方法，对基坑开挖、支护结构施工进行全方位的模拟监测，将计算设计结果和模拟计算结果进行对比验算，得出比较合理的支护结构设计方案和变形量控制方案。

根据基坑实际情况和勘查资料，本文选择的围护方案为以大直径混凝土排桩、双排桩、角撑与对顶撑相结合的内支撑为主的多种联合支护方案，结合坡顶大面积卸土减载、坑内被动区加固的措施。

计算部分主要设计计算大直径混凝土排桩（钻孔灌注桩）桩长、内力和配筋，而对卸土减载、内支撑结构、坑内被动区加固和降水设计只进行了简要的说明；flac3D模拟部分主要从建立模型、设置大直径混凝土排桩、放坡开挖、放坡坡面土钉施工、预应力锚索（代替内支撑）施工和基坑主体开挖为顺序进行建模计算，最后进行变形量监测、分析，输出桩单元、锚单元的内力分布情况并给出相应的结论与建议。

本文以常规计算和数值模拟相结合的方式进行参考对比，常规计算和数值模拟分析结果非常接近，给出了有效合理的安全系数。

关键词：深基坑支护设计flac3D模拟数值模拟AbstractWith the urbanization process ,high-rise buildings and supertall buildings are continuously emerging .As a result ,underground space development project and deep excavation project become more and more. At the same time, the scale and complexity of deep excavation increasing bigger. they make the design and construction of deep excavation to face greater challenges. So structural design and deformation prediction of deep excavation has become an important research issue in the field of geotechnical engineering. In this paper, the deep excavation of Wanda Plaza, Wuhan is studied. And using survey data and structural design of deep excavation requirements to select reasonable foundation pit ,then to conduct the corresponding design. The meantime, as checking the deformation of deep excavation is a difficult problems ,it uses flac3D numerical simulation method to monitor the progress of deep pit’s excavation, construction .Then comparing the design results of the calculation and simulation results to obtained reasonable support structure design and control program of deformation.According to the actual situation and exploration data, the envelope of large diameter piles concrete piles, angle brace and top brace on the combination of a variety of internal support-based programs are selected, combined with slope Top large dump load shedding and the reinforcement measures of pit passive zone.1) The calculation part of the paper mainly introduce the design and calculation of large diameter concrete piles or bored pile, and the rest just briefly introduce the dumping load shedding, internal support structure, the pit design of passive zone strengthening and precipitation.2) With flac3D, successively study the model building, setting large diameter concrete piles, sloping excavation, soil nailing construction, pre-stressed cable (instead of internal support) construction and excavation for the foundation pit .Finally, conduct the deformation monitoring , output pile element, the internal force distribution analysis in anchorage unit .And then, provide the corresponding conclusions and recommendations.In this paper, conventional calculations and numerical simulation methods are used. And their results were very close. So it can give an effective and reasonable safety factor through the combination of these methods.Key words: deep excavation design flac3D numerical simulation目录第一章绪论 (1)第一节选题思路 (1)第二节设计流程 (1)第二章工程概况及场地工程地质条件 (3)第一节工程概况 (3)第二节场地工程地质条件 (4)第三章A-OPQRSA段基坑支护结构设计 (10)第一节设计依据 (10)第二节设计参数 (10)第三节A-OPQRSA段基坑支护方案选择 (11)第四节A-OPQRSA段基坑减载放坡设计 (13)第五节A-OPQRSA段基坑支护桩设计 (13)第六节A-OPQRSA段基坑地下水控制方案设计 (24)第四章基于flac3D基坑开挖模拟分析 (27)第一节关于flac3D的概述 (27)第二节基坑维护方案 (27)第三节计算模型及参数 (28)第四节初始应力计算 (29)第五节支护桩施工 (31)第六节模拟分层开挖和设定锚杆 (32)第七节设置采样记录变量 (34)第八节计算结果分析 (35)第五章结论与问题 (44)第一节结论 (44)第二节设计过程中存在问题 (45)致谢 (47)参考文献 (48)附录 (49)第一章绪论第一节选题思路深基坑工程设计是当今岩土工程界关注的热点话题，深基坑工程的难题在于对变形量的预测，基坑允许的变形、垂直位移的计算是比建筑物自身允许的沉降和沉降计算更为复杂的课题，但又是基坑工程尤其是在软土地区和工程地质、水文地质复杂地区无法回避的问题。

FLAC3D在深基坑支护优化设计中的应用

FLAC3D在深基坑支护优化设计中的应用

深基坑开挖和支护是一种很复杂的工程。然而 , 现有的支护 计算理论是基于刚塑性理论的极限平衡法 , 支护结构 受 力为破坏时的极限力 , 无法得出支护结构工作状态的实 际受力状况。基坑变形只能通过经验估计 , 支护设 计在很大程 度 上依赖于实践经验 [ 1~ 3] , 因而要设计出安全、 合理、 经济的支护方 式存在很大 的难度。 数值模 拟方法由 于比理 论方法 在 计算实际受力和位移方面具有很大的优势 , 因 此在深基坑开挖前对其进行数值模拟 , 可以用来分析 基坑的变形和破坏 机 理及其成因 , 从而可以对支护结构进行优化 , 以达到基坑合理支护的效 果。本文将 FLAC3D 数值 模拟方法 应用到西安 市 骡马商业步行街深基坑支护工程的优化设计中 , 有效地克服了理论方法的不足 , 为支护工程优化设 计提供了一种新的 思 路。 1 工程概况及原支护方式 西安市骡马商业步 行街 深基 坑工 程位 于 西安 市东 大街 与东 木头 市之 间 , 基 坑南 北长 360m, 南侧 宽 110m, 北 侧 宽 130m, 占地面积约 45 000m 2, 场地地 坪标高为 407. 40~ 409. 70m, 南高北低 , 基底标高为 396. 00m, 基坑深度 11. 4m 。场地属 于黄土梁洼地貌单元 , 由于土层性质良好 , 基坑深度不大 , 考虑经济、 快速的 因素 , 因此采用 土钉墙 方式进 行支护。场 地 的地层资料如表 1, 土钉参数如表 2[ 4] 。 表1 地层编号 1 2 3 4 地层名称 填土 上更新统黄土 古土壤 中更新统黄土 埋深 ( m) 0~ 3 3~ 6 6~ 10 10~ 17. 4 地层及岩性特征 ( kN/ m 3) 17. 0 17. 8 18. 9 19. 5 C ( kPa) 13 27 35 45 ( ) 15 20 10 25 Es ( MPa) 4. 0 3. 9 6. 8 9. 8 泊松比 0. 30 0. 32 0. 35 0. 33
  1. 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。

FLAC D3深基坑的开挖与支护的命令流一、实例工程南宁地区地层属于河流阶地二元地层，广泛分布有较厚的圆砾层，国内尚无在类似地层条件下建设地铁基坑的经验，为此，可使用FLAC3D 对基坑开挖的全过程进行三维数值模拟，在对比实测数据的基础上，总结圆砾层中地铁车站深基坑的地下连续墙水平变形及周围地表沉降变形特征。

该基坑位于大学路与明秀路交叉路口处，沿大学东路东西向布置。

车站基坑长465m，标准断面宽度为20.7m，为地下两层式结构，底板埋深为15.535m（相对地面），顶板覆土厚度大于3m。

本工程主体建筑面积21163.6m2，主要结构形式为双柱三跨框架箱型结构。

本工程所处的大学路为南宁市东西向的主要交通枢纽，车流量大，人流密集，地面条件复杂。

基坑施工采用明挖顺作法施工，围护结构为800mm厚地下连续墙+内撑（三道内支撑加一道换撑）的支护体系。

第一道支撑采用钢筋混凝土支撑，尺寸为800×900mm，冠梁同时作为第一道钢筋混凝土支撑的围檩。

第二、三道支撑及换撑使用钢支撑并施加预加力，直径为609mm，壁厚为t=16mm，斜撑段采用800×1000mm钢筋砼腰梁，其余为2×I45C 钢围檩。

二、模型建立建模工作由两部分组成，实体模型部分，包括土体和地下连续墙；结构单元部分，包括混凝土支撑和钢支撑。

根据对称性原理，拟选取1/2 的实际工程尺寸进行分析。

考虑到实际的基坑长度将近500m，根据以往的经验，选取全部长度的一半虽然能够得到满意的结果，但是由于中间部分的基坑基本处于同样的受力状态，这样会使大部分的计算长度变为重复的计算，降低了计算效率。

根据初步计算结果和经验，最终确定的基坑尺寸为，宽度取基坑的最大宽度24m，开挖深度19m，基坑长度36m。

根据地勘报告，合并相似土层，模型中共划分了7个土层。

在FLAC3D 中，围护结构可以用衬砌单元(liner)或实体单元模拟。

根据Zdravdovi的研究，在二维平面基坑模拟中，分别采用实体单元和梁单元(相当于三维模型中的衬砌单元)计算所产生的墙体变形差别小于4%，而引起地表沉降的主要原因是围护结构变形造成的地层损失，可见上述两种方法计算结果的差别可忽略不计。

相比于衬砌单元(liner)，实体单元物理模型清晰，参数相对较少且容易确定，因此本例围护结构采用实体单元模拟。

模型的范围必须要足够大以覆盖基坑开挖可能影响的区域，同时又要兼顾计算效率。

根据现行国家和个地方规范的要求，参考已有的研究成果，影响范围为4 倍开挖深度，整个模型的尺寸为96×144×39（长×宽×高m）。

为减少不必要的网格，网格划分原则为基坑附近密集，远处稀疏，模型共有节点10890个，单元9408个。

模型如图所示。

命令流如下：n;开挖土体gen zone brick size 12 4 10 p0 0,0,0 p1 36,0,0 p2 0,12,0 p3 0,0,19;开挖下部土体gen zone brick size 12 4 4 p0 0,0,-20 p1 36,0,-20 p2 0,12,-20 p3 0,0,0group exc;实体地连墙gen zone brick size 12 1 10 p0 0,12,0 p1 36,12,0 p2 0,12.8,0 p3 0,12,19gen zone brick size 12 1 4 p0 0,12,-20 p1 36,12,-20 p2 0,12.8,-20 p3 0,12,0gen zone brick size 1 4 10 p0 -0.8,0,0 p1 0,0,0 p2 -0.8,12,0 p3 -0.8,0,19gen zone brick size 1 4 4 p0 -0.8,0,-20 p1 0,0,-20 p2 -0.8,12,-20 p3 -0.8,0,0gen zone brick size 1 1 10 p0 -0.8,12,0 p1 0,12,0 p2 -0.8,12.8,0 p3 -0.8,12,19gen zone brick size 1 1 4 p0 -0.8,12,-20 p1 0,12,-20 p2 -0.8,12.8,-20 p3 -0.8,12,0 group wall range group exc not;--------------------------------------------------------------------------------;土体gen zone brick size 12 5 10 p0 0,12.8,0 p1 36,12.8,0 p2 0,32,0 p3 0,12.8,19gen zone brick size 1 5 10 p0 -0.8,12.8,0 p1 0,12.8,0 p2 -0.8,32,0 p3 -0.8,12.8,19gen zone brick size 5 5 10 p0 -20.8,12.8,0 p1 -0.8,12.8,0 p2 -20.8,32,0 p3 -20.8,12.8,19 gen zone brick size 5 1 10 p0 -20.8,12,0 p1 -0.8,12,0 p2 -20.8,12.8,0 p3 -20.8,12,19 gen zone brick size 5 4 10 p0 -20.8,0,0 p1 -0.8,0,0 p2 -20.8,12,0 p3 -20.8,0,19gen zone brick size 12 5 4 p0 0,12.8,-20 p1 36,12.8,-20 p2 0,32,-20 p3 0,12.8,0gen zone brick size 1 5 4 p0 -0.8,12.8,-20 p1 0,12.8,-20 p2 -0.8,32,-20 p3 -0.8,12.8,0gen zone brick size 5 5 4 p0 -20.8,12.8,-20 p1 -0.8,12.8,-20 p2 -20.8,32,-20 p3 -20.8,12.8,0 gen zone brick size 5 1 4 p0 -20.8,12,-20 p1 -0.8,12,-20 p2 -20.8,12.8,-20 p3 -20.8,12,0 gen zone brick size 5 4 4 p0 -20.8,0,-20 p1 -0.8,0,-20 p2 -20.8,12,-20 p3 -20.8,0,0group nearsoil range group exc not group wall not;---------------------------------------------------------------------------------gen zone brick size 12 9 10 p0 0,32,0 p1 36,32,0 p2 0,92,0 p3 0,32,19gen zone brick size 1 9 10 p0 -0.8,32,0 p1 0,32,0 p2 -0.8,92,0 p3 -0.8,32,19gen zone brick size 5 9 10 p0 -20.8,32,0 p1 -0.8,32,0 p2 -20.8,92,0 p3 -20.8,32,19gen zone brick size 3 9 10 p0 -60,32,0 p1 -20.8,32,0 p2 -60,92,0 p3 -60,32,19gen zone brick size 3 5 10 p0 -60,12.8,0 p1 -20.8,12.8,0 p2 -60,32,0 p3 -60,12.8,19gen zone brick size 3 1 10 p0 -60,12,0 p1 -20.8,12,0 p2 -60,12.8,0 p3 -60,12,19gen zone brick size 3 4 10 p0 -60,0,0 p1 -20.8,0,0 p2 -60,12,0 p3 -60,0,19;深部土体gen zone brick size 12 9 4 p0 0,32,-20 p1 36,32,-20 p2 0,92,-20 p3 0,32,0gen zone brick size 1 9 4 p0 -0.8,32,-20 p1 0,32,-20 p2 -0.8,92,-20 p3 -0.8,32,0gen zone brick size 5 9 4 p0 -20.8,32,-20 p1 -0.8,32,-20 p2 -20.8,92,-20 p3 -20.8,32,0gen zone brick size 3 9 4 p0 -60,32,-20 p1 -20.8,32,-20 p2 -60,92,-20 p3 -60,32,0gen zone brick size 3 5 4 p0 -60,12.8,-20 p1 -20.8,12.8,-20 p2 -60,32,-20 p3 -60,12.8,0 gen zone brick size 3 1 4 p0 -60,12,-20 p1 -20.8,12,-20 p2 -60,12.8,-20 p3 -60,12,0gen zone brick size 3 4 4 p0 -60,0,-20 p1 -20.8,0,-20 p2 -60,12,-20 p3 -60,0,0group farsoil range group exc not group wall not group nearsoil notgen merge 0.01plot block group;镜像对称gen zone ref nor 0 -1 0;显示坐标轴plot add axes red;边界条件-------------------------------------------------------------------------fix x range x -60.1 -59.9fix x range x 35.9 36.1fix y range y 91.9 92.1fix y range y -91.9 -92.1fix z range z -20.1 -19.9;重力-----------------------------------------------------------------------------set grav 0 0 -9.81;土体材料属性----------------------------------------------------------------------model mohr;填土pro bulk 6e6 she 2.22e6 fric 9 coh 20000 dens 1850 range z 16.9 19.1;淤泥质土pro bulk 3e6 she 1.07e6 fric 4 coh 5000 dens 1830 range z 12.9 17.1;粘土pro bulk 12e6 she 4.62e6 fric 15 coh 27500 dens 1970 range z 10.9 13.1;粉土pro bulk 10e6 she 3.8e6 fric 18 coh 8000 dens 2000 range z 8.9 11.1;中砂,粗砂,砾砂pro bulk 25e6 she 9.68e6 fric 30 coh 0 dens 2050 range z 2.9 9.1;圆砾层pro bulk 32e6 she 12.8e6 fric 32 coh 0 dens 2200 range z -4.1 3.1;泥岩pro bulk 50e6 she 22.5e6 fric 70 coh 0 dens 2200 range z -21.1 -3.9;给地连墙赋本构，属性------------------------------m e range group wallini density 2500 range group wallpro young 40e9 poisson 0.167 range group wall;计算初始平衡plot set rotation 20 0 30set mech ratio 1e-4plot con szzSolveSave inidis0.sav ;保存初始平衡，后面再调用通过上述命令流可得如下图所示：图初始平衡应力三、开挖与支护该工程分四次开挖，架设三道支撑（第一道为钢筋混凝土支撑，第二、三道为钢支撑），故把工程分成四个工况。

相关文档
最新文档