【正文】 plot block groupplot add axes red gen zon brick p0 0 0 p1 0 p2 0 1 p3 0 0 amp。 主隧道開挖25mmodel null range group main1step 1000save plot set rotation 0 0 60 plot cont zdisp outl on第六章 雙線鐵路隧道施工過程分析62 V級圍巖施工過程模擬（1）建立隧道1/4圓周模型gen zon radcyl p0 0 0 0 p1 0 0 p2 0 1 0 p3 0 0 amp。 輔助隧道開挖25 m計(jì)算rest ini xdis ydis zdis model null range group service y 0,25step 1000save 。p4 7 50 8 p5 0 50 50 p6 15 0 50 p7 15 50 50 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。p5 0 50 8 p6 7 0 8 p7 7 50 8 p8 point 1 p10 point 2 amp。網(wǎng)格三維顯示plot set rotation 30 0 30 plot cont zdisp outl on421加載順序new 。網(wǎng)格三維顯示plot set rotation 30 0 30 plot cont zdisp outl on。step 750。ini sxx grad 0,0,。 quasirandom surfacedz = dz * sin( * zz + ) 。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 .1fix x range x fix z range z .1fix z range z fix y range y ini syy grad 0,0 range y 10,0ini syy 4e4 grad 0,0 range y 15,10ini syy 7e4 grad 0,0 range y 15,25step 10。網(wǎng)格顯示plot set rotation 30 0 30plot block group plot grid vel413不考慮重力影響的均勻應(yīng)力new 。網(wǎng)格顯示plot set rotation 30 0 30plot grid red fap greennewgen zone brick size 6 6 6model elasprop bulk 1e8 shear 7e7fix x range x table 1 0,0 100,1e5apply sxx hist table 1 range x , y 0,6 z 0,2hist zone sxx 6,0,0step 100。網(wǎng)格顯示plot set rotation 15 0 30plot block group410應(yīng)力邊界newgen zone brick size (4,4,4) p0 (0,0,0) p1 (4,0,0) amp。對稱得出完整模型gen zon reflect dip 0 dd 0 ori 0 0 0gen zon reflect dip 90 dd 270 ori 0 0 0。網(wǎng)格顯示plot set rotation 30 0 30plot block group solvesave plot cont zdisp outl on單一網(wǎng)格（?。┏绦騨ewgen 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 .1fix x range x .1fix x range x fix y range y .1fix y range y apply szz 1e6 range z x 0,2 y 0,2hist unbal。對稱得出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。p8 23 0 0 p9 7 0 8 p10 23 50 0 p11 7 50 8 amp。size 3 10 6 4 dim 3 3 3 3 rat 1 1 1 1。 輔助隧道網(wǎng)格gen point id 1 (,)gen point id 2 (,)gen zon radcyl p0 0 0 1 p1 7 0 0 p2 0 50 1 p3 0 0 8 p4 7 50 0 amp。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 amp。fill matrix with random numbersarray var(4,3)loop m (1,4)loop n (1,3)var(m,n) = urandendloopendloopenddef ashow 。模筑襯砌計(jì)算model elas range group concretliner anyplot block groupplot add axes redini density 2500 range group concretliner any 。求解自重應(yīng)力場model mohr ini density 1800 。繪制模型圖plot block groupplot add axes red。gr_k 2e10 gr_coh 1e10 gr_per end_mandend_loopendplace_cablessave plot sketch sel cable force red212圓形隧道開挖模擬計(jì)算。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 索單元自動(dòng)生成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。 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 。給定邊界條件fix z range z ,fix x range x ,fix x range x ,fix y range y fix y range y ,。毛洞開挖計(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ì)算分析第四章 FISH語言及建模技術(shù)41數(shù)組newdef afill 。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 shear45馬蹄形隧道網(wǎng)格newgen zone radcyl size 5 10 6 12 rat 1 1 1 amp。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。p5 7 50 0 p6 7 0 8 p7 7