隧道及地下工程FLAC解析方法 命令流Word格式文档下载.docx

1、def loadload=gp_yfunbal（ad1）+gp_yfunbal（ad2）+gp_yfunbal（ad3）+gp_yfunbal（ad4）hist loadhist gp ydis 0,2,0step 1000plot his 1 vs -22-5用FISH函数计算体积模量和剪砌模量def derives_mod = y_mod / （2.0 * （1.0 + p_ratio）b_mod = y_mod / （3.0 * （1.0 - 2.0 * p_ratio）set y_mod = 5e8 p_ratio = 0.25deriveprint b_mod print s_mod

2、 2-6 在FLAC输入中使用符号变量Newgen zone brick size 2,2,2model elasticprop bulk=b_mod shear=s_modprint zone prop bulkprint zone prop shear2-7 控制循环def xxxsum = 0prod = 1loop n （1,10）sum = sum + nprod = prod * nend_loopxxxprint sum, prodgen zone brick p0 （0,0,0） p1 （-10,0,0） p2 （0,10,0） p3 （0,0,-10）model elasde

3、f installpnt = zone_headloop while pnt #nullz_depth = -z_zcen（pnt）y_mod = y_zero + cc * sqrt（z_depth） 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） set p_ratio=0.25 y_zero=1e7 cc=1e8 install2-8 拆分命令行new ;example of a sum of ma

4、ny thingsdef long_sum temp = v1 + v2 + v3 + v4 + v5 + v6 + v7 + v8 + v9 + v10long_sum = temp + v11 + v12 + v13 + v14 + v152-9 变量类型def haveoneaa = 2bb = 3.4cc = Have a nice daydd = aa * bbee = cc + , old chaphaveoneprint fish2-10 IF条件语句if xx 0 thenabc = 1000elseabc = -1000end_ifset xx = 10set xx = 02

5、-11 索单元自动生成gen zone brick size 10 3 5def 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_commandplace_cablesplot grid sel geom redplot set rotation 15 0 60mod mohrprop bulk 1e8 shear .3e8 fric 35prop coh 1e3 tens 1e3ini dens 1000set grav 0,0,-1

6、0fix 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_cablesz_d = 5.5 - float（n）z_t = z_d + 0.5z_b = z_d - 0.5free x range x -.1,.1 z z_b z_tsel cable be

7、g 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.0save cab_end.savplot sketch sel cable force red2-12圆形隧道开挖模拟计算;建立模型gen zon radcyl p0 0

8、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 group concretliner fill group insiderockgen zon reflect dip 90 dd 90 orig 0 0 0gen zon reflect dip 0 dd 0 ori 0

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

10、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 zon 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.

11、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 add axes redplot set rotation 0 0 45 用于显示三维模型设置重力set gravity 0 0 -10给定边界条件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求解自

12、重应力场model mohr ini density 1800 ;围岩的密度prop bulk=1.47e8 shear=5.6e7 fric=20 coh=5.0e4 tension=1.0e4 ;体积、剪切、摩擦角、凝聚力、抗拉强度set mech ratio=1e-4save Gravsol.sav plot cont zdisp outl onplot cont szz毛洞开挖计算initial xdisp=0 ydisp=0 zdisp=0 model null range group insiderock any group concretliner anyset mech rat

13、io=5e-4save Kaiwsol.savplot cont zdispplot cont sdisp plot cont xzz模筑衬砌计算model elas range group concretliner anyini density 2500 range group concretliner any ;衬砌混凝土的密度prop bulk=16.67e9,shear=12.5e9 range group concretliner any ;衬砌混凝土的体积弹模、剪切弹模save zhihusol.sav完成计算分析第四章 命令流按照顺序进行4-1数组def afill ;fill

14、matrix with random numbersarray var（4,3）loop m （1,4）loop n （1,3）var（m,n） = urandendloopdef ashow ;display contents of matrixhed = msg = +string（m）hed = hed + +string（n）msg = msg + +string（var（m,n）if m = 1dum = out（hed）endifdum = out（msg）afillAshow4-2函数操作aa = 2 * 3xxx = aa + bbdef stress_sumstress_su

15、m = 0.0pnt = zone headloop while pnt # nullStress_sum = stress_sum + z_sxx（pnt）pnt=z_next（pnt）sum = 0.0sum = sum + z_sxx（pnt）pnt=z next（pnt）end loopstress_sum = sumdef 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 =

16、syz（ipz） + facH_var_6 = sxz（ipz） + fac4-3函数删除与重定义def joeii=out（This is A Function）def fred joefredii=（This is A New Function）4-4字符串def in_defxx = in（msg+（+default:+string（default）+）:）if type（xx） = 3in_def = defaultin_def = xxdef moduli_datadefault = 1.0e9msg=Input Youngs modulus Y_mod = in_defdefaul

17、t = 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 = out（ Select a different value -）moduli_dataprop bulk = b_mod shear = s_modprint p_ratio y_mod b_mod s_modpause4-5马蹄形隧道网格gen zone radcyl size 5 10 6 12 rat 1 1 1 1.2 &p0

18、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 304-6复杂形状网格生成 主隧道网格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

19、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 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 d

20、im 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 &衬砌网格sel shell 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

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

22、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对称得出完整模型gen zon reflect dip 90 dd 270 ori 0 0 0网格显示4-7网格连接gen 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

23、,2 p1 4,0,2 p2 0,4,2 p3 0,0,4attach face range z 1.9 2.1prop bulk 8e9 shear 5e9fix z range z -.1 .1fix x range x 3.9 4.1fix y range y 3.9 4.1apply szz -1e6 range z 3.9 4.1 x 0,2 y 0,2plot block group save att.sav单一网格（小）程序gen zone brick size 8 8 8 p0 0,0,0 p1 4,0,0 p2 0,4,0 p3 0,0,4save noattx.sav单一网格（大）程序gen zone brick size 4 4 4 p0 0,0,0 p1 4,0,0 p2 0,4,0 p3 0,0