【正文】 k=step 2000plot cont zdisp save 。設置底邊界fix x range x , 。 在輔助隧道上施加襯砌restore sel shell id=1 range cyl end1 0 0 1 end2 0 25 1 rad 3sel shell prop iso=(, ) thick = sel node fix y xr zr range y 。 主隧道建模gen zon radcyl p0 15 0 0 p1 23 0 0 p2 15 50 0 p3 15 0 8 amp。 topology.gp_zpos(gp) = * gp_zpos(gp) * ( + dz)gp = gp_next(gp)end_loopendmountainfix x range x .1fix x range x fix y range y .1fix y range y fix z range z .1set grav 0,0,10ini density=2000ini szz= (grad 0,0,) sxx= (grad 0,0,) syy= (grad 0,0,)step 100plot set rotation 30 0 30 plot cont zdisp outl on。p2 (0,4,0) p3 (2,0,)model elasticprop bulk 1e8 shear .3e8apply nstress 1e6 range plane dip 60 dd 270 origin ,0,0 abovestep 1。size 6 10 3 10 rat 1 1 1 gen zone brick p0 0 0 8 p1 7 0 8 p2 0 50 8 p3 0 0 50 amp。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。建立模型gen zon radcyl p0 0 0 0 p1 6 0 0 p2 0 1 0 p3 0 0 6 amp。體積、剪切、摩擦角、凝聚力、抗拉強度set mech ratio=1e4solvesave plot cont zdisp outl onplot cont szz。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)格顯示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。網(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。step p8 23 0 0 p9 7 0 8 p10 23 50 0 p11 7 50 8 amp。size 5 2 12 group outsiderock2gen zon brick p0 0 0 p1 50 0 0 p2 1 0 p3 0 amp。plot cont zdisp 。1/2繪制上下地層網(wǎng)格gen zon brick p0 0 0 p1 0 p2 0 1 p3 0 0 size 5 2 8 group outsiderock1
gen zon brick p0 0 0 p1 0 p2 0 1 p3 0 0 9 size 5 2 12 group outsiderock2。4二次襯砌rest model elas range group concretlinerxymodel elas range group concretlinerxzmodel elas range group concretlinersymodel elas range group concretlinerszini density 2600 range group concretlinerxyini density 2600 range group concretlinerxzini density 2600 range group concretlinersyini density 2600 range group concretlinerszprop bulk=, shear= range group concretlinerxyprop bulk=, shear= range group concretlinerxzprop bulk=, shear= range group concretlinersyprop bulk=, shear= range group concretlinerszstep 2000save III級圍巖施工過程模擬new。1全斷面開挖rest ini xdis ydis zdis model null range group insiderocksz any group concretlinersz any amp。 group insiderockxy any group concretlinerxy anystep 2000save 。size 5 2 10 4 dim rat 1 1 1 group outsiderocksygen zone cshell p0 0 0 0 p1 0 0 p2 0 1 0 p3 0 0 size 1 2 10 4 amp。1/2右上側和下側土體網(wǎng)格gen zon brick p0 0 p1 50 0 p2 1 p3 0 size 20 2 8 group outsiderock4
gen zon brick p0 0 p1 50 0 p2 1 p3 0 9 size 20 2 12 group outsiderock5
。繪制豎向位移場，如圖615所示plot cont szz 。size 20 2 5 group outsiderock3gen zon brick p0 0 p1 50 0 p2 1 p3 0 amp。p4 7 50 8 p5 0 50 50 p6 15 0 50 p7 15 50 50 amp。網(wǎng)格三維顯示plot set rotation 30 0 30 plot cont zdisp outl on。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。對稱得出完整模型gen zon reflect dip 0 dd 0 ori 0 0 0gen zon reflect dip 90 dd 270 ori 0 0 0。size 3 10 6 4 dim 3 3 3 3 rat 1 1 1 1。模筑襯砌計算model elas range group concretliner anyplot block groupplot add axes redini density 2500 range group concretliner any 。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 60pl