【正文】 loop 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。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。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。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。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復(fù)雜形狀網(wǎng)格生成new。 主隧道網(wǎng)格gen zon radcyl p0 15 0 0 p1 23 0 0 p2 15 50 0 p3 15 0 8 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。 輔助隧道網(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。p5 0 50 8 p6 7 0 8 p7 7 50 8 p8 point 1 p10 point 2 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。p5 7 50 0 p6 7 0 8 p7 7 50 8 p9 point 1 p11 point 2 amp。size 3 10 6 4 dim 3 3 3 3 rat 1 1 1 1。襯砌網(wǎng)格sel shell range cyl end1 0 0 1 end2 0 50 1 rad 3。 隧道外圍邊界網(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 15 50 50 p6 50 0 50 p7 50 50 50 amp。p8 23 0 0 p9 7 0 8 p10 23 50 0 p11 7 50 8 amp。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。p4 7 50 8 p5 0 50 50 p6 15 0 50 p7 15 50 50 amp。size 3 10 10 rat 1 1 。對(duì)稱得出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。對(duì)稱得出完整模型gen zon reflect dip 90 dd 270 ori 0 0 0。網(wǎng)格顯示plot set rotation 30 0 30plot block group47網(wǎng)格連接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 model 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。網(wǎng)格顯示plot set rotation 30 0 30plot block group solvesave plot cont zdisp outl on單一網(wǎng)格（小）程序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 .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。網(wǎng)格顯示plot set rotation 30 0 30plot block groupsolvesave plot cont zdisp outl on單一網(wǎng)格（大）程序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 .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。網(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 rad。對(duì)稱得出完整模型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 group49球體洞穴網(wǎng)格生成newdef parmrad= len= in_size=6 rad_size=10 endparmdef make_spherep_gp=gp_headloop while p_gpnullpx=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=(maxprad)/(lenrad)gp_xpos(p_gp)=ax+u*(bxax)gp_ypos(p_gp)=ay+u*(byay)gp_zpos(p_gp)=az+u*(bzaz)end_ifp_gp=gp_next(p_gp)end_loopendgen zone radbrick edge len size in_size in_size in_size rad_size amp。rat dim rad rad radmake_sphere。對(duì)稱得出完整模型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 15 0 30plot block group410應(yīng)力邊界newgen zone brick size (4,4,4) p0 (0,0,0) p1 (4,0,0) amp。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。網(wǎng)格顯示plot set rotation 30 0 30plot block group plot grid red fap green411改變應(yīng)力邊界條件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 red fap greennewgen zone brick size 6 6 6model elasprop bulk 1e8 shear 7e7fix x range x def x_stressx_stress = stress_inc * stependset stress_inc = 1e3apply sxx hist x_stress range x , y 0,6 z 0,2hist zone sxx 6,0,0hist x_stressstep 100。網(wǎng)格顯示plot set rotation 30 0 30plot grid red fap greennewgen zone brick size