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

第二章 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 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(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.25deriveprint 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 xxxsum = 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_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) 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 = 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 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 fric 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_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_commandend_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 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 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 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.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;给定边界条件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 cont 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 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 afill ;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)endloopendafillAshow4-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 = 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-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 Poissons 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 bulkpauseprint 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复杂形状网格生成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 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 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 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;对称得出完整模型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 .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 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 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 cont zdisp outl on4-8立方体洞穴网格生成newdef parmrad=4.0 len=10.0 in_size=6 rad_size=10 endparmgen 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 0 30plot block group4-9球体洞穴网格生成newdef parmrad=4.0 len=10.0 in_size=6 rad_size=10 endparmdef 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*kbz=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 radmake_sphere;对称得出完整模型gen zon reflect dip 0 dd 0 ori 0 0 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 abovestep 1;网格显示plot set rotation 30 0 30plot block group plot 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 100end_commandend_loopendset n_steps=100 stress_inc=-1e3hist zone sxx 6,0,0plot create sxx_histplot add hist 1plot showSuperstep;网格显示plot set rotation 30 0 30plot block group plot 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_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.9,6.1 y 0,6 z 0,2hist zone sxx 6,0,0step 100;网格显示plot set rotation 30 0 30plot grid red fap greennewgen 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 y range y 4.9 5.1ini szz -5e4 grad 0 0 -1e4app szz -5e4 range z -.1 .1solvemodel null range x 1,4 y 1 4 z 3 5step 100plot set plane dip 90 dd 180 origin 0,2.5,0plot add surf plane behind yellplot add vel plane behind blackplot set rotation 30 0 30plot block groupplot show4-12位移边界newgen zone brick size 4,4,4 p3 2,0,3.464model elasprop bulk 1e8 shear .3e8macro left_boun plane dip 60 dd 270 origin 0.1,0,0 abovemacro right_boun plane dip 60 dd 270 origin 3.9,0,0 belowapply nvel 0.1 plane dip 60 dd 270 range left_bounapply nvel 0.1 plane dip 120 dd 90 range right_bounstep 1;网格显示plot set rotation 30 0 30plot block group plot grid velnewnewgen zone brick size 4 4 4model elasticprop bulk 1e8 shear .3e8apply nstress -1e6 plane dip 0 dd 0 range y 3.9 4.1apply nvel 0.0 plane dip 0 dd 0 range y -.1 .1apply nvel 0.0 plane dip 0 dd 0 range x -.1 .1step 10;网格显示plot set rotation 30 0 30plot block groupplot grid velnewgen zone brick size 10 5 5mod elprop shear 1e8 bulk 2e8fix x y z range x -.1 .1 y 0 5 z 0 5fix x y z range x 0 10 y 0 5 z -.1 .1fix x y z range x 9.9 10.1 y 0 5 z 0 5table 1def find_addhead = nullp_gp = gp_headloop while p_gp # nullx_pos = gp_xpos(p_gp)if x_pos = width thennew = get_mem(2)mem(new) = headmem(new+1) = p_gphead = newendifp_gp = gp_next(p_gp)endloopendset width=10.0find_adddef apply_velwhile_steppingad = headloop while ad # nullp_gp = mem(ad+1)gp_xvel(p_gp) = vel_max * gp_zpos(p_gp) / heightgp_zvel(p_gp) = -vel_max * (gp_xpos(p_gp) - width) / heightad = mem(ad)endloopendset large vel_max=1e-2 height=5.0step 100;网格显示plot set rotation 30 0 30plot block group plot grid vel4-13不考虑重力影响的均匀应力new ;set grav 0 0 0gen zone brick size 6 6 6model elasini sxx=-5e6 syy=-1e7 szz=-2e7apply sxx=-5e6 range x -0.1 0.1apply sxx=-5e6 range x 5.9 6.1apply syy=-1e7 range y -0.1 0.1apply syy=-1e7 range y 5.9 6.1apply szz=-2e7 range z -0.1 0.1apply szz=-2e7 range z 5.9 6.1prop bulk 8e9 shear 5e9ini dens 2000fix x range z -.1 .1step 10plot set rotation 30 0 30 ;网格三维显示plot cont zdisp outl on4-14考虑应力梯度的均匀材料new gen zone brick size 10 10 10 p1 20,0,0 p2 0,20,0 p3 0,0,20model mohrprop bulk 5e9 shear 3e9 fric 35ini density 2500set gravity 0,0,-10fix x y z range z -0.1 0.1ini szz -5.0e6 grad 0,0,2.5e4ini syy -2.5e6 grad 0,0,1.25e4ini sxx -2.5e6 grad 0,0,1.25e4apply szz -4.5e6 range z 19.9 20.1apply szz -5.0e6 range z -0.1 0.1apply sxx -2.5e6 grad 0,0,1.25e4 range x -0.1 0.1apply sxx -2.5e6 grad 0,0,1.25e4 range x 19.9 20.1apply syy -2.5e6 grad 0,0,1.25e4 range y -0.1 0.1apply syy -2.5e6 grad 0,0,1.25e4 range y 19.9 20.1step 10;网格三维显示plot set rotation 30 0 30 plot cont zdisp outl on4-15考虑应力梯度的非均匀材料new gen zone brick size 10 10 10 &p0 0,-25,0 p1 20,-25,0 p2 0,0,0 p3 0,-25,20model elasprop bulk 5e9 shear 3e9ini density 1600 range y -10,0ini density 2000 range y -15,-10ini density 2200 range y -25,-15set gravity 0,-10,0fix x range x -.1 .1fix x range x 19.9 20.1fix z range z -.1 .1fix z range z 19.9 20.1fix y range y -25.1 -24.9ini syy 0.0 grad 0,1.6e4,0 range y -10,0ini syy 4e4 grad 0,2.0e4,0 range y -15,-10ini syy 7e4 grad 0,2.2e4,0 range y -15,-25step 10;网格三维显示plot set rotation 30 0 30 plot cont zdisp outl on4-16非均匀网格应力初始化new gen zone radcyl size 3 8 4 5 fill p1 10,0,0 p2 0,10,0 p3 0,0,10mode elasticprop shear 3e8 bulk 5e8fix x range x -.1 .1fix x range x 9.9 10.1fix y range y -.1 .1fix y range y 9.9 10.1fix z range z -.1 .1ini szz = -2.5e5 grad 0,0,2.5e4ini density 2500set grav 0,0,-10step 10;网格三维显示plot set rotation 30 0 30 plot cont zdisp outl on4-17不规则自由面应力初始化gen zone brick size 15 15 10 p0 0,0,0 edge=100.0model elasticprop shear 3e8 bulk 5e8def mountaingp = gp_headloop while gp # nullzz = sqrt(gp_xpos(gp)2 + gp_ypos(gp)2)dz = 0.06 * sin(0.2 * zz + 100.0) ; Sum Fourier terms fordz = dz + 0.06 * sin(0.22 * zz - 20.3) ; quasi-random surfacedz = dz - 0.04 * sin(0.33 * zz + 33.3) ; topology.gp_zpos(gp) = 0.5 * gp_zpos(gp) * (1.0 + dz)gp = gp_next(gp)end_loopendmountainfix x range x -.1 .1fix x range x 99.9 100.1fix y range y -.1 .1fix y range y 99.9 100.1fix z range z -.1 .1set grav 0,0,-10ini density=2000ini szz=-2.0e6 (grad 0,0,2.0e4) sxx=-4.0e6 (grad 0,0,4.0e4) syy=-4.0e6 (grad 0,0,4.0e4)step 100plot set rotation 30 0 30 plot cont zdisp outl on;网格三维显示4-18非均网格内部压实new gen zone brick size 8 8 10 ratio 1.2 1 1model mohrini dens 2000prop bulk 2e8 shear 1e8prop fric 30fix x range x -.1