隧道及地下工程FLAC解析方法命令流.docVIP

第四章 命令流按照顺序进行 4-1数组 new def afill ;fill matrix with random numbers array var(4,3) loop m (1,4) loop n (1,3) var(m,n) = urand endloop endloop end def ashow ;display contents of matrix loop m (1,4) hed = ’ ’ msg = ’ ’+string(m) loop n (1,3) hed = hed + ’ ’+string(n) msg = msg + ’ ’+string(var(m,n)) endloop if m = 1 dum = out(hed) endif dum = out(msg) endloop end afill Ashow 4-2函数操作 new def xxx aa = 2 * 3 xxx = aa + bb end new def stress_sum stress_sum = 0.0 pnt = zone head loop while pnt # null Stress_sum = stress_sum + z_sxx(pnt) pnt=z_next(pnt) end_loop end new def stress_sum sum = 0.0 pnt = zone head loop while pnt # null sum = sum + z_sxx(pnt) pnt=z next(pnt) end loop stress_sum = sum end new def h_var_1 ipz = z_near(1,2,3) H_var_1 = sxx(ipz) + fac H_var_2 = syy(ipz) + fac H_var_3 = szz(ipz) + fac H_var_4 = sxy(ipz) + fac H_var_5 = syz(ipz) + fac H_var_6 = sxz(ipz) + fac end 4-3函数删除与重定义 new def joe ii=out(‘This is A Function’) end def fred joe end fred def joe ii=(‘This is A New Function’) end fred 4-4字符串 new def in_def xx = in(msg+’(’+’default:’+string(default)+’):’) if type(xx) = 3 in_def = default else in_def = xx endif end ; def moduli_data default = 1.0e9 msg=’Input Young‘s modulus ’ Y_mod = in_def default = 0.25 msg=’Input Poisson‘s ratio ’ p_ratio = in_def if p_ratio = 0.5 then ii = out(’ Bulk mod is undefined at Poisson‘s ratio = 0.5’) ii = out(’ Select a different value --’) p_ratio = in_def endif s_mod = y_mod / (2.0 * (1.0 + p_ratio)) b_mod = y_mod / (3.0 * (1.0 - 2.0 * p_ratio)) end ; moduli_data gen zone brick size 2,2,2 model elastic prop bulk = b_mod shear = s_mod print p_ratio y_mod b_mod s_mod pause print zone prop bulk pause print zone prop shear 4-5马蹄形隧道网格 new gen zone radcyl size 5 10 6 12 rat 1 1 1 1.2 p0 0,0,0 p1 100,0,0 p2 0,200,0 p3 0,0,100 gen zone radtun size 5 10 5 12 rat 1 1 1 1.2 p0 0,0,0 p1 0,0,-100 p2 0,200,0 p3 100,0,0 gen zone reflect dip 90 dd 270 origin 0,0,0 plot set rotation 30 0 30 plot block group 4-6复杂形状