【正文】 第二章 FLAC3D原理及入門指南21定義一個(gè)FISH函數(shù)newdef abcabc = 25 * 3 + 5Endprint abc22使用一個(gè)變量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ù)計(jì)算體積模量和剪砌模量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 。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 + v15end29 變量類型newdef haveoneaa = 2bb = cc = ’Have a nice day’dd = aa * bbee = cc + ’, old chap’endhaveoneprint fish210 IF條件語句newdef abcif xx 0 thenabc = 1000elseabc = 1000end_ifendset xx = 10print abcset xx = 0print abc211 索單元自動生成newgen zone brick size 10 3 5plot set rotation 0 0 45plot block groupdef place_cables loop n (1,5)z_d = float(n) mandsel cable beg ,z_d end ,z_d nseg 7end_mandend_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 .1fix y range y .1fix y range y fix x range x .1fix x range x set largehist unbalsolvesave ini xdis 0 ydis 0 zdis 0hist gp xdisp 0,1,5def place_cablesloop n (1,5)z_d = float(n)z_t = z_d + z_b = z_d mandfree x range x ,.1 z z_b z_tsolvesel cable beg ,z_d end ,z_d nseg 7sel cable beg ,z_d end ,z_d nseg 7sel cable beg ,z_d end ,z_d nseg 7sel cable prop emod 2e10 ytension 1e8 xcarea amp。gr_k 2e10 gr_coh 1e10 gr_per end_mandend_loopendplace_cablessave plot sketch sel cable force red212圓形隧道開挖模擬計(jì)算。建立模型gen zon radcyl p0 0 0 0 p1 6 0 0 p2 0 1 0 p3 0 0 6 amp。 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。 size 1 2 8 4 dim 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 z gen zon reflect dip 90 dd 90 orig 0 0 0 range x z gen zon reflect dip 90 dd 90 orig 0 0 0 range x z 。繪制模型圖plot block groupplot add axes red。plot set rotation 0 0 45 用于顯示三維模型。設(shè)置重力set gravity 0 0 10。給定邊界條件fix z range z ,fix x range x ,fix x range x ,fix y range y fix y range y ,。求解自重應(yīng)力場model mohr ini density 1800 。圍巖的密度prop bulk= shear= fric=20 coh= tension= 。體積、剪切、摩擦角、凝聚力、抗拉強(qiáng)度set mech ratio=1e4solvesave plot cont zdisp outl onplot cont szz。毛洞開挖計(jì)算initial xdisp=0 ydisp=0 zdisp=0 model null range group insiderock any group concretliner anyplot block groupplot add axes redset mech ratio=5e4solvesave plot cont zdispplot cont sdisp plot cont szzplot cont xzz。模筑襯砌計(jì)算model elas range group concretliner anyplot block groupplot add axes redini density 2500 range group concretliner any 。襯砌混凝土的密度prop bulk=,shear= range group concretliner any 。襯砌混凝土的體積彈模、剪切彈模set mech ratio=1e4solvesave plot cont zdispplot cont sdisp plot cont szzplot cont xzz。完成計(jì)算分析第四章 FISH語言及建模技術(shù)41數(shù)組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)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 head