《《隧道与地下工程FLAC解析方法》全部命令流汇总.docx》由会员分享，可在线阅读，更多相关《《隧道与地下工程FLAC解析方法》全部命令流汇总.docx（44页珍藏版）》请在三一文库上搜索。
1、.第二章 FLAC3D 原理及入门指南2-1 定义一个 FISH 函数newdef abcabc = 25 * 3 + 5Endprint abc2-2 使用一个变量newdef abchh = 25abc = hh * 3 + 5EndPrint hhPrint abc2-3 对变量和函数的理解newdef abchh = 25abc = hh * 3 + 5Endset abc=0 hh=0print hhprint abcprint hhnewdef abcabc = hh * 3 + 5endset hh=25print abcset abc=0 hh=0print hhprint a
2、bcprint hh.2-4 获取变量的历史记录newgen zone brick size 1 2 1model mohrprop shear=1e8 bulk=2e8 cohes=1e5 tens=1e10fix x y z range y -0.1 0.1apply yvel -1e-5 range y 1.9 2.1plot set rotation 0 0 45plot block groupdef get_adad1 = gp_near(0,2,0)ad2 = gp_near(1,2,0)ad3 = gp_near(0,2,1)ad4 = gp_near(1,2,1)endget_
3、addef loadload=gp_yfunbal(ad1)+gp_yfunbal(ad2)+gp_yfunbal(ad3)+gp_yfunbal(ad4)endhist loadhist gp ydis 0,2,0step 1000plot his 1 vs -22-5 用 FISH 函数计算体积模量和剪砌模量newdef derives_mod = y_mod / (2.0 * (1.0 + p_ratio)b_mod = y_mod / (3.0 * (1.0 - 2.0 * p_ratio)endset y_mod = 5e8 p_ratio = 0.25deriveprint b_m
4、odprint s_mod2-6 在 FLAC输入中使用符号变量Newdef derives_mod = y_mod / (2.0 * (1.0 + p_ratio)b_mod = y_mod / (3.0 * (1.0 - 2.0 * p_ratio)endset y_mod = 5e8 p_ratio = 0.25.derivegen zone brick size 2,2,2model elasticprop bulk=b_mod shear=s_modprint zone prop bulkprint zone prop shear2-7 控制循环Newdef xxxsum = 0pr
5、od = 1loop n (1,10)sum = sum + nprod = prod * nend_loopendxxxprint sum, prodnewgen zone brick p0 (0,0,0) p1 (-10,0,0) p2 (0,10,0) p3 (0,0,-10)model elasplot set rotation 0 0 45plot block groupdef installpnt = zone_headloop while pnt #nullz_depth = -z_zcen(pnt)y_mod = y_zero + cc * sqrt(z_depth)z_pro
6、p (pnt,shear ) = y_mod / (2.0*(1.0+p_ratio)z_prop(pnt,bulk ) = y_mod / (3.0*(1.0-2.0*p_ratio)pnt = zone_next(pnt)end_loopendset p_ratio=0.25 y_zero=1e7 cc=1e8install2-8 拆分命令行new ;example of a sum of many thingsdef long_sumtemp = v1 + v2 + v3 + v4 + v5 + v6 + v7 + v8 + v9 + v10.long_sum = temp + v11
7、+ v12 + v13 + v14 + v15end2-9 变量类型newdef haveoneaa = 2bb = 3.4cc = Have a nice day dd = aa * bbee = cc +, old chapendhaveoneprint fish2-10 IF条件语句newdef abcif xx 0 thenabc = 1000elseabc = -1000end_ifendset xx = 10print abcset xx = 0print abc2-11 索单元自动生成newgen zone brick size 10 3 5plot set rotation 0
8、 0 45plot block groupdef place_cablesloop n (1,5)z_d = float(n) - 0.5commandsel cable beg 0.0,1.5,z_d end 7.0,1.5,z_d nseg 7end_commandend_loop.endplace_cablesplot grid sel geom rednewgen zone brick size 10 3 5plot set rotation 15 0 60plot block groupmod mohrprop bulk 1e8 shear .3e8 fric 35prop coh
9、1e3 tens 1e3ini dens 1000set grav 0,0,-10fix x y z range z -.1 .1fix y range y -.1 .1fix y range y 2.9 3.1fix x range x -.1 .1fix x range x 9.9 10.1set largehist unbalsolvesave cab_str.savini xdis 0 ydis 0 zdis 0hist gp xdisp 0,1,5def place_cablesloop n (1,5)z_d = 5.5 - float(n)z_t = z_d + 0.5z_b =
10、z_d - 0.5commandfree x range x -.1,.1 z z_b z_tsolvesel cable beg 0.0,0.5,z_d end 7.0,0.5,z_d nseg 7sel cable beg 0.0,1.5,z_d end 7.0,1.5,z_d nseg 7sel cable beg 0.0,2.5,z_d end 7.0,2.5,z_d nseg 7sel cable prop emod 2e10 ytension 1e8 xcarea 1.0 &gr_k 2e10 gr_coh 1e10 gr_per 1.0end_commandend_loopend
11、place_cablessave cab_end.savplot sketch sel cable force red.2-12 圆形隧道开挖模拟计算; 建立模型gen zon radcyl p0 0 0 0 p1 6 0 0 p2 0 1 0 p3 00 6 &size 4 28 4 dim 3 3 3 3 rat 1 1 1 1.2group outsiderockgen zone cshell p0 0 0 0 p1 3 0 0 p2 0 1 0 p3 0 0 3 &size 1 2 84 dim 2.7 2.7 2.7 2.7 rat 1 1 1 1 group concretline
12、r fillgroup insiderockgen zon reflect dip 90 dd 90 orig 0 0 0gen zon reflect dip 0 dd 0 ori 0 0 0gen zon brickp00 0 6 p1 6 0 6p2 0 1 6p3 0 0 13 size 4 2 6 groupoutsiderock1gen zon brickp00 0 -12 p1 6 0 -12 p2 0 1 -12 p3 0 0 -6 size 4 2 5 groupoutsiderock2gen zon brickp06 0 0 p1 21 0 0p2 6 1 0p3 6 0
13、6 size 10 2 4 groupoutsiderock3gen zon reflect dip 0 dd 0 orig 0 0 0 range group outsiderock3gen zon brick p0 6 0 6 p1 21 0 6 p2 6 1 6 p3 6 0 13 size 10 2 6 group outsiderock4 gen zon brick p0 6 0 -12 p1 21 0 -12 p2 6 1 -12 p3 6 0 -6 size 10 2 5 group outsiderock5gen zon reflect dip 90 dd 90 orig 0
14、0 0 range x -0.1 6.1 z 6.1 13.1 gen zon reflect dip 90 dd 90 orig 0 0 0 range x -0.1 6.1 z -6.1 -12.1 gen zon reflect dip 90 dd 90 orig 0 0 0 range x 6.1 21.1 z -12.1 13.1; 绘制模型图plot block groupplot add axes red;plot set rotation 0 0 45用于显示三维模型; 设置重力set gravity 0 0 -10; 给定边界条件fix z range z -12.01,-1
15、1.99fix x range x -21.01,-20.99fix x range x 20.99,21.01fix y range y -0.01 0.01fix y range y 0.99,1.01; 求解自重应力场model mohrini density 1800 ;围岩的密度prop bulk=1.47e8shear=5.6e7fric=20coh=5.0e4 tension=1.0e4; 体积、剪切、摩擦角、凝聚力、抗拉强度set mech ratio=1e-4solvesave Gravsol.savplot cont zdisp outl onplot cont szz.;
16、 毛洞开挖计算initial xdisp=0 ydisp=0 zdisp=0model null range group insiderock any group concretliner anyplot block groupplot add axes redset mech ratio=5e-4solvesave Kaiwsol.savplot cont zdispplot cont sdispplot cont szzplot cont xzz; 模筑衬砌计算model elas range group concretliner anyplot block groupplot add a
17、xes redini density 2500 range group concretliner any;衬砌混凝土的密度prop bulk=16.67e9,shear=12.5e9 range group concretliner any ;衬砌混凝土的体积弹模、剪切弹模set mech ratio=1e-4solvesave zhihusol.savplot cont zdispplot cont sdispplot cont szzplot cont xzz; 完成计算分析第四章 FISH 语言及建模技术4-1 数组newdef afill ;fill matrix with rando
18、m numbersarray var(4,3)loop m (1,4)loop n (1,3)var(m,n) = urandendloopendloopenddef ashow ;display contents of matrixloop m (1,4).hed = msg = +string(m)loop n (1,3)hed = hed+ +string(n)msg = msg + +string(var(m,n)endloopif m = 1dum = out(hed)endifdum = out(msg)endloopendafillAshow4-2 函数操作newdef xxxa
19、a =2*3 xxx = aa + bbendnewdef stress_sumstress_sum = 0.0pnt = zone headloop while pnt # nullStress_sum = stress_sum + z_sxx(pnt)pnt=z_next(pnt)end_loopendnewdef stress_sumsum = 0.0pnt = zone headloop while pnt # nullsum = sum + z_sxx(pnt)pnt=z next(pnt)end loopstress_sum = sumend.newdef h_var_1ipz =
20、 z_near(1,2,3)H_var_1 = sxx(ipz) + facH_var_2 = syy(ipz) + facH_var_3 = szz(ipz) + facH_var_4 = sxy(ipz) + facH_var_5 = syz(ipz) + facH_var_6 = sxz(ipz) + facend4-3 函数删除与重定义newdef joeii=out( This is A Function )enddef fredjoeendfreddef joeii=( This is A New Function )endfred4-4 字符串newdef in_defxx =
21、in(msg+ ( +default:+string(default)+): )if type(xx) = 3in_def = defaultelsein_def = xxendifend;def moduli_datadefault = 1.0e9msg=Input Young s modulusY_mod = in_defdefault = 0.25msg=Input Poissons ratio.p_ratio = in_defif p_ratio = 0.5 thenii = out( Bulk mod is undefined at Poissons ratio = 0.5)ii =
22、 out( Select a different value - )p_ratio = in_defendifs_mod = y_mod / (2.0 * (1.0 + p_ratio)b_mod = y_mod / (3.0 * (1.0 - 2.0 * p_ratio)end;moduli_datagen zone brick size 2,2,2model elasticprop bulk = b_mod shear = s_modprint p_ratio y_mod b_mod s_modpauseprint zone prop bulkpauseprint zone prop sh
23、ear4-5 马蹄形隧道网格newgen zone radcyl size 5 10 6 12 rat 1 1 1 1.2 &p0 0,0,0 p1 100,0,0 p2 0,200,0 p3 0,0,100gen zone radtun size 5 10 5 12 rat 1 1 1 1.2 &p0 0,0,0 p1 0,0,-100 p2 0,200,0 p3 100,0,0gen zone reflect dip 90 dd 270 origin 0,0,0plot set rotation 30 0 30plot block group4-6 复杂形状网格生成new; 主隧道网格ge
24、n zon radcyl p0 15 0 0 p1 23 0 0 p2 15 50 0 p3 15 0 8 &size 4 10 6 4 dim 4 4 4 4 rat 1 1 1 1 fillgen zon reflect dip 90 dd 90 orig 15 0 0gen zon reflect dip 0 ori 0 0 0; 辅助隧道网格gen point id 1 (2.969848,0.0,-0.575736)gen point id 2 (2.969848,50.0,-0.575736)gen zon radcyl p0 0 0 -1 p1 7 0 0 p2 0 50 -1
25、p3 0 0 8 p4 7 50 0 &p5 0 50 8 p6 7 0 8 p7 7 50 8 p8 point 1 p10 point 2 &.size 3 10 6 4 dim 3 3 3 3 rat 1 1 1 1gen zon radcyl p0 0 0 -1 p1 0 0 -8 p2 0 50 -1 p3 7 0 0 p4 0 50 -8 & p5 7 50 0 p6 7 0 -8 p7 7 50 -8 p9 point 1 p11 point 2 &size 3 10 6 4 dim 3 3 3 3 rat 1 1 1 1; 衬砌网格sel shell range cyl end
26、1 0 0 -1 end2 0 50 -1 rad 3; 隧道外围边界网格gen zone radtun p0 7 0 0 p1 50 0 0 p2 7 50 0 p3 15 0 50 p4 50 50 0 & p5 15 50 50 p6 50 0 50 p7 50 50 50 &p8 23 0 0 p9 7 0 8 p10 23 50 0 p11 7 50 8 & size 6 10 3 10 rat 1 1 1 1.1gen zone brick p0 0 0 8 p1 7 0 8 p2 0 50 8 p3 0 0 50 &p4 7 50 8 p5 0 50 50 p6 15 0 50
27、p7 15 50 50 & size 3 10 10 rat 1 1 1.1; 对称得出 1/2 模型gen zon reflect dip 0 ori 0 0 0 range x 0 23 y 0 50 z 8 50gen zon reflect dip 0 ori 0 0 0 range x 23 50 y 0 50 z 0 50; 建立主隧道和辅助隧道块名称group service range cyl end1 0 0 -1 end2 0 50 -1 rad 3group main range cyl end1 15 0 0 end2 15 50 0 rad 4; 对称得出完整模型ge
28、n zon reflect dip 90 dd 270 ori 0 0 0; 网格显示plot set rotation 30 0 30plot block group4-7 网格连接newgen zone brick size 4 4 4 p0 0,0,0 p1 4,0,0 p2 0,4,0 p3 0,0,2gen zone brick size 8 8 4 p0 0,0,2 p1 4,0,2 p2 0,4,2 p3 0,0,4attach face range z 1.9 2.1model elasprop bulk 8e9 shear 5e9fix z range z -.1 .1fix
29、 x range x -.1 .1fix x range x 3.9 4.1fix y range y -.1 .1fix y range y 3.9 4.1apply szz -1e6 range z 3.9 4.1 x 0,2 y 0,2hist unbal; 网格显示plot set rotation 30 0 30plot block group.solvesave att.savplot cont zdisp outl on单一网格（小）程序newgen zone brick size 8 8 8 p0 0,0,0 p1 4,0,0 p2 0,4,0 p3 0,0,4model el
30、asprop bulk 8e9 shear 5e9fix z range z -.1 .1fix x range x -.1 .1fix x range x 3.9 4.1fix y range y -.1 .1fix y range y 3.9 4.1apply szz -1e6 range z 3.9 4.1 x 0,2 y 0,2hist unbal; 网格显示plot set rotation 30 0 30plot block groupsolvesave noattx.savplot cont zdisp outl on单一网格（大）程序newgen zone brick size
31、 4 4 4 p0 0,0,0 p1 4,0,0 p2 0,4,0 p3 0,0,4model elasprop bulk 8e9 shear 5e9fix z range z -.1 .1fix x range x -.1 .1fix x range x 3.9 4.1fix y range y -.1 .1fix y range y 3.9 4.1apply szz -1e6 range z 3.9 4.1 x 0,2 y 0,2hist unbal; 网格显示plot set rotation 30 0 30plot block groupsolvesave noattx.savplot
32、 cont zdisp outl on.4-8 立方体洞穴网格生成newdef parmrad=4.0len=10.0in_size=6rad_size=10endparmgen zone radbrick edge len size in_size in_size in_size rad_size & rat 1.0 1.0 1.0 1.2 dim rad rad rad; 对称得出完整模型gen zon reflect dip 0 dd 0 ori 0 0 0gen zon reflect dip 90 dd 270 ori 0 0 0; 网格显示plot set rotation 30
33、0 30plot block group4-9 球体洞穴网格生成newdef parmrad=4.0len=10.0in_size=6rad_size=10endparmdef make_spherep_gp=gp_headloop while p_gp#nullpx=gp_xpos(p_gp)py=gp_ypos(p_gp)pz=gp_zpos(p_gp)dist=sqrt(px*px+py*py+pz*pz)if dist0 thenk=rad/distax=px*kay=py*kaz=pz*kmaxp=max(px,max(py,pz)k=len/maxpbx=px*kby=py*k.b
34、z=pz*ku=(maxp-rad)/(len-rad)gp_xpos(p_gp)=ax+u*(bx-ax)gp_ypos(p_gp)=ay+u*(by-ay)gp_zpos(p_gp)=az+u*(bz-az)end_ifp_gp=gp_next(p_gp)end_loopendgen zone radbrick edge len size in_size in_size in_size rad_size & rat 1.0 1.0 1.0 1.2 dim rad rad rad make_sphere; 对称得出完整模型gen zon reflect dip 0 dd 0 ori 0 0
35、0gen zon reflect dip 90 dd 270 ori 0 0 0; 网格显示plot set rotation 15 0 30plot block group4-10 应力边界newgen zone brick size (4,4,4) p0 (0,0,0) p1 (4,0,0) &p2 (0,4,0) p3 (2,0,3.464)model elasticprop bulk 1e8 shear .3e8apply nstress -1e6 range plane dip 60 dd 270 origin 0.1,0,0 above step 1; 网格显示plot set r
36、otation 30 0 30plot block groupplot grid red fap green4-11 改变应力边界条件newgen zone brick size 6 6 6model elasprop bulk 1e8 shear 7e7fix x range x -0.1 0.1def supersteploop ns (1,n_steps)x_stress = stress_inccommandapply sxx add x_stress range x 5.9,6.1 y 0,6 z 0,2step 100.end_commandend_loopendset n_ste
37、ps=100 stress_inc=-1e3hist zone sxx 6,0,0plot create sxx_histplot add hist 1plot showSuperstep; 网格显示plot set rotation 30 0 30plot block groupplot grid red fap greennewgen zone brick size 6 6 6model elasprop bulk 1e8 shear 7e7fix x range x -0.1 0.1def x_stressx_stress = stress_inc * stependset stress
38、_inc = -1e3apply sxx 1.0 hist x_stress range x 5.9,6.1 y 0,6 z 0,2hist zone sxx 6,0,0hist x_stressstep 100; 网格显示plot set rotation 30 0 30plot grid red fap greennewgen zone brick size 6 6 6model elasprop bulk 1e8 shear 7e7fix x range x -0.1 0.1table 1 0,0 100,-1e5apply sxx 1.0 hist table 1 range x 5.
39、9,6.1 y 0,6 z 0,2hist zone sxx 6,0,0step 100;网格显示plot set rotation 30 0 30plot grid red fap greennew.gen zone brick size 6,6,6 p1 6,0,-1model elasprop bulk 8e9 shear 5e9apply sxx -2e6 range x -0.1 0.1apply sxx -2e6 range x 5.9 6.1step 500plot bound fapplot grid dispnewgen zone brick size 5,5,5model elasprop bulk 8e9 shear 5e9set grav 0 0 -10ini dens 1000fix x range x -.1 .1fix x range x 4.9 5.1fix y range y -.1 .1fix