【正文】 end1 0 0 1 end2 0 50 1 rad 3。p5 0 50 8 p6 7 0 8 p7 7 50 8 p8 point 1 p10 point 2 amp。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 group46復雜形狀網(wǎng)格生成new。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)endloopendafillAshow42函數(shù)操作newdef xxxaa = 2 * 3xxx = aa + bbendnewdef stress_sumstress_sum = pnt = zone headloop while pnt nullStress_sum = stress_sum + z_sxx(pnt)pnt=z_next(pnt)end_loopendnewdef stress_sumsum = pnt = 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) + facend43函數(shù)刪除與重定義newdef joeii=out(‘This is A Function’)enddef fred joeendfreddef joeii=(‘This is A New Function’)endfred44字符串newdef in_defxx = in(msg+’(’+’default:’+string(default)+’):’)if type(xx) = 3in_def = defaultelsein_def = xxendifend。襯砌混凝土的密度prop bulk=,shear= range group concretliner any 。圍巖的密度prop bulk= shear= fric=20 coh= tension= 。plot set rotation 0 0 45 用于顯示三維模型。建立模型gen zon radcyl p0 0 0 0 p1 6 0 0 p2 0 1 0 p3 0 0 6 amp。第二章 FLAC3D原理及入門指南21定義一個FISH函數(shù)newdef abcabc = 25 * 3 + 5Endprint abc22使用一個變量newdef abchh = 25abc = hh * 3 + 5EndPrint hhPrint abc23對變量和函數(shù)的理解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 hh24獲取變量的歷史記錄newgen zone brick size 1 2 1model mohrprop shear=1e8 bulk=2e8 cohes=1e5 tens=1e10fix x y z range y apply yvel 1e5 range y plot 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 225用FISH函數(shù)計算體積模量和剪砌模量newdef derives_mod = y_mod / ( * ( + p_ratio))b_mod = y_mod / ( * ( * p_ratio))endset y_mod = 5e8 p_ratio = deriveprint b_mod print s_mod 26 在FLAC輸入中使用符號變量Newdef derives_mod = y_mod / ( * ( + p_ratio))b_mod = y_mod / ( * ( * p_ratio))endset y_mod = 5e8 p_ratio = derivegen zone brick size 2,2,2model elasticprop bulk=b_mod shear=s_modprint zone prop bulkprint zone prop shear27 控制循環(huán)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 / (*(+p_ratio)) z_prop(pnt, ’bulk’) = y_mod / (*(*p_ratio)) pnt = zone_next(pnt) end_loopendset p_ratio= y_zero=1e7 cc=1e8 install28 拆分命令行new 。 size 4 2 8 4 dim 3 3 3 3 rat 1 1 1 group outsiderockgen zone cshell p0 0 0 0 p1 3 0 0 p2 0 1 0 p3 0 0 3 amp。設(shè)置重力set gravity 0 0 10。體積、剪切、摩擦角、凝聚力、抗拉強度set mech ratio=1e4solvesave plot cont zdisp outl onplot cont szz。襯砌混凝土的體積彈模、剪切彈模set mech ratio=1e4solvesave plot cont zdispplot cont sdisp plot cont szzplot cont xzz。def moduli_datadefault = msg=’Input Young‘s modulus ’Y_mod = in_defdefault = msg=’Input Poisson‘s ratio ’p_ratio = in_defif p_ratio = thenii = out(’ Bulk mod is undefined at Poisson‘s ratio = ’)ii = out(’ Select a different value ’)p_ratio = in_defendifs_mod = y_mod / ( * ( + p_ratio))b_mod = y_mod / ( * ( * p_ratio))end。 主隧道網(wǎng)格gen zon radcyl p0 15 0 0 p1 23 0 0 p2 15 50 0 p3 15 0 8 amp。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 amp。 隧道外圍邊界網(wǎng)格gen zone radtun p0 7 0 0 p1 50 0 0 p2 7 50 0 p3 15 0 50 p4 50 50 0 amp。p4 7 50 8 p5 0 50 50 p6 15 0 50 p7 15 50 50 amp。對稱得出完整模型gen zon reflect dip 90 dd 270 ori 0 0 0。網(wǎng)格顯示plot set rotation 30 0 30plot block groupsolvesave plot cont zdisp outl on48立方體洞穴網(wǎng)格生成newdef parmrad= len= in_size=6 rad_size=10 endparmgen zone radbrick edge len size in_size in_size in_size rad_size amp。rat dim rad rad radmake_sphere。網(wǎng)格顯示plot set rotation 30 0 30plot block group plot grid red fap green411改變應力邊界條件newgen zone brick size 6 6 6model elasprop bulk 1e8 shear 7e7fix x range x def supersteploop ns (1,n_steps)x_stress = stress_incmandapply sxx add x_stress range x , y 0,6 z 0,2step 100end_mandend_loopendset n_steps=100 stress_inc=1e3hist zone sxx 6,0,0plot create sxx_histplot add hist 1plot showSuperstep。網(wǎng)格顯示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 apply nvel plane dip 0 dd 0 range y .1apply nvel plane dip 0 dd 0 range x .1step 10。網(wǎng)格三維顯示plot cont zdisp outl on414考慮應力梯度的均勻材料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 ini szz grad 0,0,ini syy grad 0,0,ini sxx grad 0,0,apply szz range z apply szz range z apply sxx grad 0,0, range x apply sxx grad 0,0, range x apply syy grad 0,0, range y apply syy grad 0,0, range y step 10。網(wǎng)格三維顯示plot set rotation 30 0