《隧道及地下工程FLAC解析方法》全部命令流汇总.docx

1、隧道及地下工程FLAC解析方法全部命令流汇总第二章 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=0

2、print hhprint abcprint hh2-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(

3、1,2,1)endget_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.25deriv

4、eprint b_mod print s_mod 2-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.25derivegen zone brick size 2,2,2model elasticprop bulk=b_mod shear=s_modprint zone prop bulkprint zone prop shear2-7 控制循环Newdef xx

5、xsum = 0prod = 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_d

6、epth) z_prop(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=1e8 install2-8 拆分命令行new ;example of a sum of many thingsdef long_sum temp = v1 + v2 + v3 + v4 + v5 + v6 + v7 + v8 + v9 + v10long_sum

7、 = temp + v11 + v12 + v13 + v14 + v15end2-9 变量类型newdef haveoneaa = 2bb = 3.4cc = Have a nice daydd = 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

8、set rotation 0 0 45plot block groupdef place_cables loop 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_loopendplace_cablesplot grid sel geom rednewgen zone brick size 10 3 5plot set rotation 15 0 60plot block groupmod mohrprop bulk 1e8 shear .3e8 fr

9、ic 35prop coh 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

10、_d + 0.5z_b = 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_comm

11、andend_loopendplace_cablessave cab_end.savplot sketch sel cable force red2-12圆形隧道开挖模拟计算;建立模型gen zon radcyl p0 0 0 0 p1 6 0 0 p2 0 1 0 p3 0 0 6 & size 4 2 8 4 dim 3 3 3 3 rat 1 1 1 1.2 group outsiderockgen zone cshell p0 0 0 0 p1 3 0 0 p2 0 1 0 p3 0 0 3 & size 1 2 8 4 dim 2.7 2.7 2.7 2.7 rat 1 1 1 1

12、group concretliner fill group insiderockgen zon reflect dip 90 dd 90 orig 0 0 0gen zon reflect dip 0 dd 0 ori 0 0 0gen zon brick p0 0 0 6 p1 6 0 6 p2 0 1 6 p3 0 0 13 size 4 2 6 group outsiderock1gen zon brick p0 0 0 -12 p1 6 0 -12 p2 0 1 -12 p3 0 0 -6 size 4 2 5 group outsiderock2gen zon brick p0 6

13、0 0 p1 21 0 0 p2 6 1 0 p3 6 0 6 size 10 2 4 group outsiderock3gen 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 outsiderock4gen zon brick p0 6 0 -12 p1 21 0 -12 p2 6 1 -12 p3 6 0 -6 size 10 2 5 group outsiderock5gen zo

14、n reflect dip 90 dd 90 orig 0 0 0 range x -0.1 6.1 z 6.1 13.1gen zon reflect dip 90 dd 90 orig 0 0 0 range x -0.1 6.1 z -6.1 -12.1gen 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;给定边

15、界条件fix z range z -12.01,-11.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 mohr ini density 1800 ;围岩的密度prop bulk=1.47e8 shear=5.6e7 fric=20 coh=5.0e4 tension=1.0e4 ;体积、剪切、摩擦角、凝聚力、抗拉强度set mech ratio=1e-4solvesave Gravsol.sav plot con

16、t zdisp outl onplot cont szz;毛洞开挖计算initial xdisp=0 ydisp=0 zdisp=0 model null range group insiderock any group concretliner anyplot block groupplot add axes redset mech ratio=5e-4solvesave Kaiwsol.savplot cont zdispplot cont sdisp plot cont szzplot cont xzz;模筑衬砌计算model elas range group concretliner

17、anyplot block groupplot add axes 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 sdisp plot cont szzplot cont xzz;完成计算分析第四章 FISH语言及建模技术4-1数组newdef a

18、fill ;fill matrix with random 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)endloopenda

19、fillAshow4-2函数操作newdef xxxaa = 2 * 3xxx = 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

20、 = sumendnewdef h_var_1ipz = 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 fred joeendfreddef joeii=(This is A New Function)endfred4

21、-4字符串newdef in_defxx = in(msg+(+default:+string(default)+):)if type(xx) = 3in_def = defaultelsein_def = xxendifend;def moduli_datadefault = 1.0e9msg=Input Youngs modulus Y_mod = in_defdefault = 0.25msg=Input Poissons ratio p_ratio = in_defif p_ratio = 0.5 thenii = out( Bulk mod is undefined at Poiss

22、ons ratio = 0.5)ii = 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 bulkpaus

23、eprint zone prop shear4-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复杂形

24、状网格生成new; 主隧道网格gen 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

25、7 0 0 p2 0 50 -1 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

26、 range cyl end1 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

27、 p6 15 0 50 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;对称得出

28、完整模型gen 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 .

29、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 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

30、 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 cont zdisp outl on单一网格（大）程序newgen zone brick size 4 4 4 p0 0,0,0