地层结构法命令流
地下结构的地层结构计算方法
模型建立要点
midas地层结构算例
第七章 地层结构法的适用性
位移清零
模型建立要点
midas地层结构算例
初始地 应力场
计算开挖边界 等效结点力
删除开挖网格 反向施加结点力
确定释放系数
第七章 地层结构法的适用性
荷载分步释放 与围岩特性
岩爆
模型建立要点
midas地层结构算例
高地应 力
0.7m
E砼＝23Gpa A砼＝0.28m2 I砼＝0.00183m4
E钢＝210Gpa A钢＝39.578×10-4m2 I钢＝2500×10-8m4
E A ＝ E砼 A砼+ E钢A钢/S E I ＝ E砼 I砼+ E钢I钢/S
取E ＝ E砼
A ＝ A砼+ E钢A钢/(SE砼) =0.3316 I ＝ I砼+ E钢I钢/(SE砼) =0.002155
岩土材料
• 根据岩土性质和计算目的选择适合的本构模型。 • 定量分析时应注意材料参数的确定，必要时采用反分析。
结构材料
• 弹性或弹塑性 • 初期支护内的钢拱架与喷射砼一般视为整体计算
加固地层材料
• 直接模拟 • 不模拟，作为安全储备 • 提高地层材料参数
第七章 地层结构法的适用性
边界条件
模型建立要点
576个四边形单元
35个梁单元
第七章 地层结构法的适用性
模型建立要点
midas地层结构算例
地层与结构连接
公共节点，变形协调
. . . 1 node . A. B.
不同节点，相互独立
. . .. . 2 nodes . A. B.
摩擦接触，接触单元
地层结构法
地层结构法
地层结构法是一种用于分析地层结构的研究方法，也被称为构造地质学。
它旨在研究地质构造，通过对岩石的观察和研究，推断出其在时空上的变化。
构造地质学的基本原理是，岩石不是以静态的形式存在的，而是在构造活动的作用下，在时空上发生变化的。
研究的重点是，通过解释地质构造的演化，来解释该地区的地质发展历史。
为了更准确地探索地层结构，地质学家们开发了许多研究方法，其中最重要的是地层结构法。
它是通过对岩石结构和外观形态的研究，来分析其历史演化，推断出其在时空维度上的变化。
地层结构法分为宏观和微观两个层次，宏观层次是指通过现场调查，描述地层构造的总体规律，微观层次则是指通过极细的观察，对细节进行研究，以深入了解地层构造的结构和细节。
在宏观观测上，地质学家可以通过对岩石的形态、分布特征、岩性及有机质等征象来解释地层构造的演化。
因为，其形态特征可用于推断地层构造的时空变化，分布特征可用于推断其发育过程，岩性及有机质特征可用于推断构造环境和特殊情况。
此外，微观观测也能够帮助我们更深入地研究地层构造。
通过研究岩石的粒度、颜色、纹理等特征，可以推断出其形成的构造环境，以及其历史发展过程。
通过对地层结构的观测和研究，可以推断出地层构造的演化、发育历史和形成环境，进而解释当地的地质构造。
而这一方法，也帮助了我们更好地理解地质历史，并进一步探索其地质结构。
地层结构法已经成为构造地质学研究中一种重要的理论和实践手段，在探寻地质构造的演化及其时空变化的过程中，得到了越来越多的应用。
它不但有助于我们更好地理解地质历史，而且也为更深入地探寻地质结构提供了可能。
地层、地层单位及其表示方法
地层、地层单位及其表示方法
地层是指地质起源的地下层次结构，是2500多万年来多次地壳构造活动的产物。
地层单位是指组成地层的基本构成，地层单位类型有岩性、沉积性和构造性。
地层的测量和表示形式有体积法、断面法和悬面法。
体积法是比较客观的测量地层的方法，即用地层体积表示地层厚度。
断面法是通过断裂面测量地层厚度，一般需要在地下取断层好的地方进行采集，以确定地层厚度。
悬面法是在露头单元中测量地层厚度，或者从野外边界到内部之间测量的地层厚度。
地层的表示形式有二维、三维和四维。
最常见的地层表示形式是二维，它是按照地层厚度在地面上放射或线性表示的。
三维表示形式是在断裂体积进行表示的，四维表示形式是加入了时间因素的三维表示形式，这是因为地层也受到环境和其他因素的影响。
总之，地层是地球显著特征之一，地层单位组成地层，而地层的表示形式又分为二维、三维和四维，以确定地层真实存在情况。
隧道及地下工程FLAC解析方法命令流
隧道及地下工程FLAC解析方法命令流第四章命令流按照顺序进行4-1数组newdef afill ;fill matrix with random numbers array 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 + ’ ’+s tring(n)msg = msg + ’ ’+string(var(m,n))endloopif m = 1dum = out(hed)endifdum = out(msg)endloopendafillAshow4-2函数操作newdef xxxaa = 2 * 3xxx = aa + bbendnewdef stress_sumstress_sum = 0.0pnt = zone headloop while pnt # nullStress_sum = stress_sum + z_sxx(pnt) pnt=z_next(pnt)end_loopendnewdef stress_sumsum = 0.0pnt = 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) + facend4-3函数删除与重定义newdef joeii=out(‘This is A Function’)enddef fredjoeendfreddef joeii=(‘This is A New Function’)endfred4-4字符串newdef in_defxx = in(msg+’(’+’default:’+string(default)+’):’) if type(xx) = 3in_def = defaultelsein_def = xxendifend;def moduli_datadefault = 1.0e9msg=’Input Young‘s modulus ’Y_mod = in_defdefault = 0.25msg=’Input Poisson‘s ratio ’p_ratio = in_defif p_ratio = 0.5 thenii = out(’ Bulk mod is undefined at Poisson‘s ratio = 0.5’) ii = out(’ Select a different val ue --’)p_ratio = in_defendifs_mod = y_mod / (2.0 * (1.0 + p_ratio))b_mod = y_mod / (3.0 * (1.0 - 2.0 * 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 shear4-5马蹄形隧道网格newgen 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,100gen 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,0gen zone reflect dip 90 dd 270 origin 0,0,0plot set rotation 30 0 30plot block group4-6复杂形状网格生成new; 主隧道网格gen zon radcyl p0 15 0 0 p1 23 0 0 p2 15 50 0 p3 15 0 8 & 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; 辅助隧道网格gen point id 1 (2.969848,0.0,-0.575736)gen point id 2 (2.969848,50.0,-0.575736)gen zon radcyl p0 0 0 -1 p1 7 0 0 p2 0 50 -1 p3 0 0 8 p4 7 50 0 &p5 0 50 8 p6 7 0 8 p7 7 50 8 p8 point 1 p10 point 2 & 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 &p5 7 50 0 p6 7 0 -8 p7 7 50 -8 p9 point 1 p11 point 2 & size 3 10 6 4 dim 3 3 3 3 rat 1 1 1 1;衬砌网格sel shell range cyl end1 0 0 -1 end2 0 50 -1 rad 3; 隧道外围边界网格gen zone radtun p0 7 0 0 p1 50 0 0 p2 7 50 0 p3 15 0 50 p4 50 50 0 & p5 15 50 50 p6 50 0 50 p7 50 50 50 &p8 23 0 0 p9 7 0 8 p10 23 50 0 p11 7 50 8 &size 6 10 3 10 rat 1 1 1 1.1gen zone brick p0 0 0 8 p1 7 0 8 p2 0 50 8 p3 0 0 50 &p4 7 50 8 p5 0 50 50 p6 15 0 50 p7 15 50 50 &size 3 10 10 rat 1 1 1.1;对称得出1/2模型gen zon reflect dip 0 ori 0 0 0 range x 0 23 y 0 50 z 8 50 gen 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 3 group main range cyl end1 15 0 0 end2 15 50 0 rad 4;对称得出完整模型gen zon reflect dip 90 dd 270 ori 0 0 0;网格显示plot set rotation 30 0 30plot block group4-7网格连接newgen zone brick size 4 4 4 p0 0,0,0 p1 4,0,0 p2 0,4,0 p3 0,0,2 gen zone brick size 8 8 4 p0 0,0,2 p1 4,0,2 p2 0,4,2 p3 0,0,4 attach face range z 1.9 2.1model elasprop bulk 8e9 shear 5e9fix z range z -.1 .1fix x range x -.1 .1fix x range x 3.9 4.1fix y range y -.1 .1fix y range y 3.9 4.1apply szz -1e6 range z 3.9 4.1 x 0,2 y 0,2hist unbal;网格显示plot set rotation 30 0 30plot block groupsolvesave att.savplot cont zdisp outl on单一网格（小）程序newgen zone brick size 8 8 8 p0 0,0,0 p1 4,0,0 p2 0,4,0 p3 0,0,4 model elasprop bulk 8e9 shear 5e9fix z range z -.1 .1fix x range x -.1 .1fix x range x 3.9 4.1fix y range y -.1 .1fix y range y 3.9 4.1apply szz -1e6 range z 3.9 4.1 x 0,2 y 0,2hist unbal;网格显示plot set rotation 30 0 30plot block groupsolvesave noattx.savplot cont zdisp outl on单一网格（大）程序newgen zone brick size 4 4 4 p0 0,0,0 p1 4,0,0 p2 0,4,0 p3 0,0,4 model elasprop bulk 8e9 shear 5e9fix z range z -.1 .1fix x range x -.1 .1fix x range x 3.9 4.1fix y range y -.1 .1fix y range y 3.9 4.1apply szz -1e6 range z 3.9 4.1 x 0,2 y 0,2hist unbal;网格显示plot set rotation 30 0 30plot block groupsolvesave noattx.savplot cont zdisp outl on4-8立方体洞穴网格生成newdef parmrad=4.0len=10.0in_size=6rad_size=10endparmgen zone radbrick edge len size in_size in_size in_size rad_size &rat 1.0 1.0 1.0 1.2 dim rad rad rad;对称得出完整模型gen zon reflect dip 0 dd 0 ori 0 0 0gen zon reflect dip 90 dd 270 ori 0 0 0;网格显示plot set rotation 30 0 30plot block group4-9球体洞穴网格生成newdef parmrad=4.0len=10.0in_size=6rad_size=10endparmdef make_spherep_gp=gp_headloop while p_gp#nullpx=gp_xpos(p_gp)py=gp_ypos(p_gp)pz=gp_zpos(p_gp)dist=sqrt(px*px+py*py+pz*pz)if dist>0 thenk=rad/distax=px*kay=py*kaz=pz*kmaxp=max(px,max(py,pz))k=len/maxpbx=px*kby=py*kbz=pz*ku=(maxp-rad)/(len-rad)gp_xpos(p_gp)=ax+u*(bx-ax)gp_ypos(p_gp)=ay+u*(by-ay)gp_zpos(p_gp)=az+u*(bz-az)end_ifp_gp=gp_next(p_gp)end_loopendgen zone radbrick edge len size in_size in_size in_sizerad_size & rat 1.0 1.0 1.0 1.2 dim rad rad radmake_sphere;对称得出完整模型gen zon reflect dip 0 dd 0 ori 0 0 0gen zon reflect dip 90 dd 270 ori 0 0 0;网格显示plot set rotation 15 0 30plot block group4-10应力边界newgen zone brick size (4,4,4) p0 (0,0,0) p1 (4,0,0) &p2 (0,4,0) p3 (2,0,3.464)model elasticprop bulk 1e8 shear .3e8apply nstress -1e6 range plane dip 60 dd 270 origin 0.1,0,0 above step 1;网格显示plot set rotation 30 0 30plot block groupplot grid red fap green4-11改变应力边界条件newgen zone brick size 6 6 6model elasprop bulk 1e8 shear 7e7fix x range x -0.1 0.1def supersteploop ns (1,n_steps)x_stress = stress_inccommandapply sxx add x_stress range x 5.9,6.1 y 0,6 z 0,2 step 100end_commandend_loopendset n_steps=100 stress_inc=-1e3hist zone sxx 6,0,0plot create sxx_histplot add hist 1plot showSuperstep;网格显示plot set rotation 30 0 30plot block groupplot grid red fap greennewgen zone brick size 6 6 6model elasprop bulk 1e8 shear 7e7fix x range x -0.1 0.1def x_stressx_stress = stress_inc * stependset stress_inc = -1e3apply sxx 1.0 hist x_stress range x 5.9,6.1 y 0,6 z 0,2 hist zone sxx 6,0,0hist x_stressstep 100;网格显示plot set rotation 30 0 30plot grid red fap greennewgen zone brick size 6 6 6model elasprop bulk 1e8 shear 7e7fix x range x -0.1 0.1table 1 0,0 100,-1e5apply sxx 1.0 hist table 1 range x 5.9,6.1 y 0,6 z 0,2 hist zone sxx 6,0,0step 100;网格显示plot set rotation 30 0 30plot grid red fap greennewgen zone brick size 6,6,6 p1 6,0,-1model elasprop bulk 8e9 shear 5e9apply sxx -2e6 range x -0.1 0.1apply sxx -2e6 range x 5.9 6.1step 500plot bound fapplot grid dispnewgen zone brick size 5,5,5model elasprop bulk 8e9 shear 5e9set grav 0 0 -10ini dens 1000fix x range x -.1 .1fix x range x 4.9 5.1fix y range y -.1 .1fix y range y 4.9 5.1ini szz -5e4 grad 0 0 -1e4app szz -5e4 range z -.1 .1solvemodel null range x 1,4 y 1 4 z 3 5step 100plot set plane dip 90 dd 180 origin 0,2.5,0plot add surf plane behind yellplot add vel plane behind blackplot set rotation 30 0 30plot block groupplot show4-12位移边界newgen zone brick size 4,4,4 p3 2,0,3.464model elasprop bulk 1e8 shear .3e8macro left_boun ’plane dip 60 dd 270 origin 0.1,0,0 above’macro ri ght_boun ’plane dip 60 dd 270 origin 3.9,0,0 below’apply nvel 0.1 plane dip 60 dd 270 range left_boun apply nvel 0.1 plane dip 120 dd 90 range right_bounstep 1;网格显示plot set rotation 30 0 30plot block groupplot grid velnewnewgen zone brick size 4 4 4model elasticprop bulk 1e8 shear .3e8apply nstress -1e6 plane dip 0 dd 0 range y 3.9 4.1 apply nvel 0.0 plane dip 0 dd 0 range y -.1 .1 apply nvel 0.0 plane dip 0 dd 0 range x -.1 .1 step 10;网格显示plot set rotation 30 0 30plot block groupplot grid velnewgen zone brick size 10 5 5mod elprop shear 1e8 bulk 2e8fix x y z range x -.1 .1 y 0 5 z 0 5fix x y z range x 0 10 y 0 5 z -.1 .1fix x y z range x 9.9 10.1 y 0 5 z 0 5table 1def find_addhead = nullp_gp = gp_headloop while p_gp # nullx_pos = gp_xpos(p_gp)if x_pos = width thennew = get_mem(2)mem(new) = headmem(new+1) = p_gphead = newendifp_gp = gp_next(p_gp)endloopendset width=10.0find_adddef apply_velwhile_steppingad = headloop while ad # nullp_gp = mem(ad+1)gp_xvel(p_gp) = vel_max * gp_zpos(p_gp) / heightgp_zvel(p_gp) = -vel_max * (gp_xpos(p_gp) - width) / height ad = mem(ad)endloopendset large vel_max=1e-2 height=5.0step 100;网格显示plot set rotation 30 0 30plot block groupplot grid vel4-13不考虑重力影响的均匀应力new ;set grav 0 0 0gen zone brick size 6 6 6model elasini sxx=-5e6 syy=-1e7 szz=-2e7apply sxx=-5e6 range x -0.1 0.1apply sxx=-5e6 range x 5.9 6.1apply syy=-1e7 range y -0.1 0.1apply syy=-1e7 range y 5.9 6.1apply szz=-2e7 range z -0.1 0.1apply szz=-2e7 range z 5.9 6.1prop bulk 8e9 shear 5e9ini dens 2000fix x range z -.1 .1step 10plot set rotation 30 0 30 ;网格三维显示plot cont zdisp outl on4-14考虑应力梯度的均匀材料newgen zone brick size 10 10 10 p1 20,0,0 p2 0,20,0 p3 0,0,20 model mohrprop bulk 5e9 shear 3e9 fric 35ini density 2500set gravity 0,0,-10fix x y z range z -0.1 0.1ini szz -5.0e6 grad 0,0,2.5e4ini syy -2.5e6 grad 0,0,1.25e4ini sxx -2.5e6 grad 0,0,1.25e4apply szz -4.5e6 range z 19.9 20.1apply szz -5.0e6 range z -0.1 0.1apply sxx -2.5e6 grad 0,0,1.25e4 range x -0.1 0.1apply sxx -2.5e6 grad 0,0,1.25e4 range x 19.9 20.1apply syy -2.5e6 grad 0,0,1.25e4 range y -0.1 0.1apply syy -2.5e6 grad 0,0,1.25e4 range y 19.9 20.1step 10;网格三维显示plot set rotation 30 0 30plot cont zdisp outl on4-15考虑应力梯度的非均匀材料newgen zone brick size 10 10 10 &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 -.1 .1fix x range x 19.9 20.1fix z range z -.1 .1fix z range z 19.9 20.1fix y range y -25.1 -24.9ini syy 0.0 grad 0,1.6e4,0 range y -10,0ini syy 4e4 grad 0,2.0e4,0 range y -15,-10ini syy 7e4 grad 0,2.2e4,0 range y -15,-25step 10;网格三维显示plot set rotation 30 0 30plot cont zdisp outl on4-16非均匀网格应力初始化newgen zone radcyl size 3 8 4 5 fill p1 10,0,0 p2 0,10,0 p3 0,0,10 mode elasticprop shear 3e8 bulk 5e8fix x range x -.1 .1fix x range x 9.9 10.1fix y range y -.1 .1fix y range y 9.9 10.1fix z range z -.1 .1ini szz = -2.5e5 grad 0,0,2.5e4ini density 2500set grav 0,0,-10step 10;网格三维显示plot set rotation 30 0 30plot cont zdisp outl on4-17不规则自由面应力初始化gen zone brick size 15 15 10 p0 0,0,0 edge=100.0model elasticprop shear 3e8 bulk 5e8def mountaingp = gp_headloop while gp # nullzz = sqrt(gp_xpos(gp)?2 + gp_ypos(gp)?2)dz = 0.06 * sin(0.2 * zz + 100.0) ; Sum Fourier terms fordz = dz + 0.06 * sin(0.22 * zz - 20.3) ; quasi-random surface dz = dz - 0.04 * sin(0.33 * zz + 33.3) ; topology.gp_zpos(gp) = 0.5 * gp_zpos(gp) * (1.0 + dz)gp = gp_next(gp)end_loopendmountainfix x range x -.1 .1fix x range x 99.9 100.1fix y range y -.1 .1fix y range y 99.9 100.1fix z range z -.1 .1set grav 0,0,-10ini density=2000ini szz=-2.0e6 (grad 0,0,2.0e4) sxx=-4.0e6 (grad 0,0,4.0e4) syy=-4.0e6 (grad 0,0,4.0e4) step 100plot set rotation 30 0 30plot cont zdisp outl on;网格三维显示4-18非均网格内部压实newgen zone brick size 8 8 10 ratio 1.2 1 1 model mohrini dens 2000prop bulk 2e8 shear 1e8fix x range x -.1 .1fix x range x 7.9 8.1fix y range y -.1 .1fix y range y 7.9 8.1fix z range z -.1 .1;ini szz -2.0e5 grad 0,0,2e4;ini sxx -1.5e5 grad 0,0,1.5e4;ini syy -1.5e5 grad 0,0,1.5e4set grav 10step 1000;pause;prop tens 1e10 coh 1e10;step 750;prop tens 0 coh 0;step 250;网格三维显示plot set rotation 30 0 30plot cont zdisp outl on4-19模型改变后初始应力变化newgen zone brick size 5 5 5model elasprop sh 2e8 bu 3e8fix x y z range z -.1 .1set grav 0 0 -10ini dens 2000Solve;网格三维显示plot set rotation 30 0 30plot cont zdisp outl on;模型改变model mohr range x 0 2 y 0 5 z 0 2prop sh 2e8 bu 3e8 fric 35 range x 0 2 y 0 5 z 0 2 Step 2000 ;网格三维显示plot set rotation 30 0 30plot cont zdisp outl on4-20应力与孔隙压力的初始化newconfig fluidgen zone brick size 8 5 10model elasmodel fl_isoprop bulk 1e9 shear 5e8prop poros 0.5 perm 1e-10ini fmod 2e9ini fdensity 1e3ini sat 0ini sat 1 range z -.1 5.1set grav 0 0 -10fix x range x -.1 .1fix x range x 7.9 8.1fix y range y -.1 .1fix y range y 4.9 5.1fix z range z -.1 .1ini pp 5.e4 grad 0,0,-1.e4 range z 0.0 5.ini szz -20e4 grad 0,0,20e3ini szz add -1.5e4 grad 0,0,.25e4 range z 5,6 ini szz add -2.5e4 grad 0,0,.5e4 range z 0,5 ini szz add -.25e4 range z 0,5 solve;网格三维显示plot set rotation 30 0 30plot cont zdisp outl onnewgen zone brick size 1 1 10model elasini dens 2500 range z 0 5ini dens 2250 range z 5 6ini dens 2000 range z 6 10prop bulk 1e9 shear 5e8set grav 10water dens 1000water table ori 0 0 5 normal 0 0 1fix x range x -.1 .1fix x range x 7.9 8.1fix y range y -.1 .1fix y range y 4.9 5.1fix z range z -.1 .1ini szz -20e4 grad 0,0,20e3ini szz add -1.5e4 grad 0,0,.25e4 range z 5,6 ini szz add -2.5e4 grad 0,0,.5e4 range z 0,5 ini szz add -.25e4 range z 0,5 solve;网格三维显示plot set rotation 30 0 30plot cont zdisp outl on4-21加载顺序new; 主隧道建模gen zon radcyl p0 15 0 0 p1 23 0 0 p2 15 50 0 p3 15 0 8 & 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;辅助隧道建模gen point id 1 (2.969848,0.0,-0.575736)gen point id 2 (2.969848,50.0,-0.575736)gen zon radcyl p0 0 0 -1 p1 7 0 0 p2 0 50 -1 p3 0 0 8 p4 7 50 0 &p5 0 50 8 p6 7 0 8 p7 7 50 8 p8 point 1 p10 point 2 &size 3 10 6 4 dim 3 3 3 3 rat 1 1 1 1 fillgen zon radcyl p0 0 0 -1 p1 0 0 -8 p2 0 50 -1 p3 7 0 0 p4 0 50 -8 &p5 7 50 0 p6 7 0 -8 p7 7 50 -8 p9 point 1 p11 point 2 &size 3 10 6 4 dim 3 3 3 3 rat 1 1 1 1 fill;隧道边界网格gen zone radtun p0 7 0 0 p1 50 0 0 p2 7 50 0 p3 15 0 50 p4 50 50 0 & p5 15 50 50 p6 50 0 50 p7 50 50 50 &p8 23 0 0 p9 7 0 8 p10 23 50 0 p11 7 50 8 &size 6 10 3 10 rat 1 1 1 1.1gen zone brick p0 0 0 8 p1 7 0 8 p2 0 50 8 p3 0 0 50 &p4 7 50 8 p5 0 50 50 p6 15 0 50 p7 15 50 50 &size 3 10 10 rat 1 1 1.1gen 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 50group main1 range cyl end1 15 0 0 end2 15 25 0 rad 4 group main2 range cyl end1 15 25 0 end2 15 50 0 rad 4 group service range cyl end1 0.0,0.0,-0.575736 &end2 0.0,50.0,-0.575736 rad 3.0save tun0.savrest tun0.sav; 初始应力场模量model mohrprop shear 0.36e9 bulk 0.6e9 coh 1e5 fric 20 tens 1e5 apply szz -1.4e6 range z 49.9 50.1fix z range z -50.1 -49.1fix x range x -.1 .1fix x range x 49.9 50.1fix y range y -.1 .1fix y range y 49.9 50.1ini sxx -1.4e6 syy -1.4e6 szz -1.4e6hist unbalhist gp xdis 3,0,-1hist gp zdis 0,0,2hist gp xdis 3,25,-1hist gp zdis 0,25,2step 1000save tun1.sav; 辅助隧道开挖25 m计算rest tun1.savini xdis 0.0 ydis 0.0 zdis 0.0model null range group service y 0,25step 1000save tun2.sav; 在辅助隧道上施加衬砌restore tun2.savsel shell id=1 range cyl end1 0 0 -1 end2 0 25 -1 rad 3 sel shell prop iso=(25.3e9, 0.266) thick = 0.5sel node fix y xr zr range y -0.1 0.1 ; symmetry cond. sel node fix x yr zr range x -0.1 0.1 ; symmetry cond. model mohr range group service y 25 50prop shear 0.36e9 bulk 0.6e9 coh 1e5 fric 20 tens 1e5 ; 主隧道开挖25mmodel null range group main1step 1000save tun3.savplot set rotation 0 0 60plot cont zdisp outl on。
地层学原理---第五章 地层层序 地层结构 地层
2） 二级层序：mesosequence（中层序）
该层序代表地球一个中等周期的变化，这种变化是 由太阳学穿越银道面周期发生的，其时限，据Innanen （1978）计算为67Ma，即太阳系两次穿越银道面的间 隔时间为33.5Ma，Stothers(1985)提出了（31±3）Ma和 （33±3）Ma两个值；Wang and shi(1996)认为最短时 间段为35—36Ma，史晓颖把35Ma称为地质历史上一个 重要的自然临界周期。二级层序不仅表现在全球海平 面变化、生物群绝灭、构造强动幕、岩浆活动、古气 侯、古地磁侧转等，而且还表现在古地理的演化上。
2. 李勇等（1994）提出3种岩石地层结构：进积型、加积型和 退积型，陆相地层均决于冲积体系进积、加积、退积的 演化过程。
3. 魏家庸等（1991）提出韵律结构（Rhythmic texture）、旋 回 结 构 （ Cyclic textare ） 、 均 一 结 构 （ Homogeneous texture）、加积结构（Aggradational textune）、退积结构 （ Retrogradational texture ） 和 进积 结 构 （ Progradational texture）
（二） 岩石地层结构的构成单元
（1） 岩层（bed）； （2） 岩层组（bed sets）； （3） 基本层序（essential sequence）； （4） 基本层序组（essential sequence set相当于
Parasequence)Байду номын сангаас
FLAC-3D深基坑的开挖与支护的命令流
FLAC D3深基坑的开挖与支护的命令流一、实例工程南宁地区地层属于河流阶地二元地层，广泛分布有较厚的圆砾层，国内尚无在类似地层条件下建设地铁基坑的经验，为此，可使用FLAC3D 对基坑开挖的全过程进行三维数值模拟，在对比实测数据的基础上，总结圆砾层中地铁车站深基坑的地下连续墙水平变形及周围地表沉降变形特征。
该基坑位于大学路与明秀路交叉路口处，沿大学东路东西向布置。
车站基坑长465m，标准断面宽度为20.7m，为地下两层式结构，底板埋深为15.535m（相对地面），顶板覆土厚度大于3m。
本工程主体建筑面积21163.6m2，主要结构形式为双柱三跨框架箱型结构。
本工程所处的大学路为南宁市东西向的主要交通枢纽，车流量大，人流密集，地面条件复杂。
基坑施工采用明挖顺作法施工，围护结构为800mm厚地下连续墙+内撑（三道内支撑加一道换撑）的支护体系。
第一道支撑采用钢筋混凝土支撑，尺寸为800×900mm，冠梁同时作为第一道钢筋混凝土支撑的围檩。
第二、三道支撑及换撑使用钢支撑并施加预加力，直径为609mm，壁厚为t=16mm，斜撑段采用800×1000mm钢筋砼腰梁，其余为2×I45C 钢围檩。
二、模型建立建模工作由两部分组成，实体模型部分，包括土体和地下连续墙；结构单元部分，包括混凝土支撑和钢支撑。
根据对称性原理，拟选取1/2 的实际工程尺寸进行分析。
考虑到实际的基坑长度将近500m，根据以往的经验，选取全部长度的一半虽然能够得到满意的结果，但是由于中间部分的基坑基本处于同样的受力状态，这样会使大部分的计算长度变为重复的计算，降低了计算效率。
根据初步计算结果和经验，最终确定的基坑尺寸为，宽度取基坑的最大宽度24m，开挖深度19m，基坑长度36m。
根据地勘报告，合并相似土层，模型中共划分了7个土层。
在FLAC3D 中，围护结构可以用衬砌单元(liner)或实体单元模拟。
根据Zdravdovi的研究，在二维平面基坑模拟中，分别采用实体单元和梁单元(相当于三维模型中的衬砌单元)计算所产生的墙体变形差别小于4%，而引起地表沉降的主要原因是围护结构变形造成的地层损失，可见上述两种方法计算结果的差别可忽略不计。
PLAXIS 建模方法汇总 - 筑信达
\图1 三维模型在2D AE 至2016版本和3D 2013至 2016版本中，建模主要包含两个方面，即地层建模和结构建模（结构单元以及荷载）。
1 地层建模地层建模常用方法主要有两种，即钻孔和外部文件导入（VIP）。
1.1钻孔一般情况下，特别是建立三维模型时，为了简化建模以及让网格划分更顺利，都要对实际地层进行简化。
此时，地层建模主要利用钻孔工具建模。
通过定义钻孔顶部和底部标高，可以生成当前钻孔土层厚度。
也可以利用多个钻孔创建非水平地层，如图2和图3所示。
钻孔之间利用内插的方法生成地层。
同时在每个钻孔中可以分别定义水头水力条件，生成初始地下水位。
图2 编辑钻孔图3 多钻孔生成的非水平地层1.2外部文件导入（VIP）如果想要创建勘察报告中提供的地质模型，一般通过导入外部文件的方法生成。
此时要将外部文件的几何信息进行处理，处理成PLAXIS 可识别的离散化三角形对象。
特别要注意不同软件版本对导入文件的几何信息要求不同，3D2016版本可以之间导入体，无需离散化三角形对象。
导入的文件格式支持3ds、dxf、txt（2D）、dwg 等，其中3D 2016不再支持dwg 格式。
导入形式包括导入土层表面和整个土层，如图4和图5所示。
图4 导入地层表面（3D）PLAXIS 建模方法汇总筑信达 孙立超总体来说，建模技术也是一种艺术。
不同的人对同样的工程项目会有不同的简化结果。
而简化的结果对网格划分和计算时间成本以及计算结果都有很大影响。
建模的一个基本原则就是简化的模型能够反映项目分析的主要问题。
在这个基本原则下，尤其是对于PLAXIS 3D 来说，对于以提交精确计算结果为目的时，模型不要贪大求全，能简化的可以适当简化，比如，不均匀分布地层的等效，排桩等效成板，建筑物等效成荷载等。
如果为了使得模型更加真实的反映工程项目三维模型兼顾计算时，建立的地层反映出最主要的不均匀分布情况，建立的建筑物或构筑物对象反映出主要传力体系，能用结构单元模拟的构件最好用结构单元模拟。
《隧道及地下工程FLAC解析方法》全部命令流汇总..
第二章FLAC3D原理及入门指南2-1定义一个FISH函数newdef abcabc = 25 * 3 + 5Endprint abc2-2使用一个变量newdef abchh = 25abc = hh * 3 + 5EndPrint hhPrint abc2-3对变量和函数的理解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 hh2-4获取变量的历史记录newgen zone brick size 1 2 1model mohrprop shear=1e8 bulk=2e8 cohes=1e5 tens=1e10fix x y z range y -0.1 0.1apply yvel -1e-5 range y 1.9 2.1plot 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 -22-5用FISH函数计算体积模量和剪砌模量newdef derives_mod = y_mod / (2.0 * (1.0 + p_ratio))b_mod = y_mod / (3.0 * (1.0 - 2.0 * p_ratio))endset y_mod = 5e8 p_ratio = 0.25deriveprint b_modprint s_mod2-6 在FLAC输入中使用符号变量Newdef derives_mod = y_mod / (2.0 * (1.0 + p_ratio))b_mod = y_mod / (3.0 * (1.0 - 2.0 * p_ratio))endset y_mod = 5e8 p_ratio = 0.25derivegen zone brick size 2,2,2model elasticprop bulk=b_mod shear=s_modprint zone prop bulkprint zone prop shear2-7 控制循环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 / (2.0*(1.0+p_ratio))z_prop(pnt, ’bulk’) = y_mod / (3.0*(1.0-2.0*p_ratio))pnt = zone_next(pnt)end_loopendset p_ratio=0.25 y_zero=1e7 cc=1e8install2-8 拆分命令行new ;example of a sum of many thingsdef long_sumtemp = v1 + v2 + v3 + v4 + v5 + v6 + v7 + v8 + v9 + v10long_sum = temp + v11 + v12 + v13 + v14 + v15end2-9 变量类型newdef haveoneaa = 2bb = 3.4cc = ’Have a nice day’dd = aa * bbee = cc + ’, old chap’endhaveoneprint fish2-10 IF条件语句newdef abcif xx > 0 thenabc = 1000elseabc = -1000end_ifendset xx = 10print abcset xx = 0print abc2-11 索单元自动生成newgen zone brick size 10 3 5plot set rotation 0 0 45plot block groupdef place_cablesloop n (1,5)z_d = float(n) - 0.5commandsel cable beg 0.0,1.5,z_d end 7.0,1.5,z_d nseg 7 end_commandend_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 -.1 .1fix y range y -.1 .1fix y range y 2.9 3.1fix x range x -.1 .1fix x range x 9.9 10.1set largehist unbalsolvesave cab_str.savini xdis 0 ydis 0 zdis 0hist gp xdisp 0,1,5def place_cablesloop n (1,5)z_d = 5.5 - float(n)z_t = z_d + 0.5z_b = z_d - 0.5commandfree x range x -.1,.1 z z_b z_tsolvesel cable beg 0.0,0.5,z_d end 7.0,0.5,z_d nseg 7sel cable beg 0.0,1.5,z_d end 7.0,1.5,z_d nseg 7sel cable beg 0.0,2.5,z_d end 7.0,2.5,z_d nseg 7sel cable prop emod 2e10 ytension 1e8 xcarea 1.0 & gr_k 2e10 gr_coh 1e10 gr_per 1.0end_commandend_loopendplace_cablessave cab_end.savplot sketch sel cable force red2-12圆形隧道开挖模拟计算;建立模型gen zon radcyl p0 0 0 0 p1 6 0 0 p2 0 1 0 p3 0 0 6 &size 4 2 8 4 dim 3 3 3 3 rat 1 1 1 1.2 group outsiderockgen zone cshell p0 0 0 0 p1 3 0 0 p2 0 1 0 p3 0 0 3 &size 1 2 8 4 dim 2.7 2.7 2.7 2.7 rat 1 1 1 1 group concretliner fill group insiderock gen 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 outsiderock1 gen zon brick p0 0 0 -12 p1 6 0 -12 p2 0 1 -12 p3 0 0 -6 size 4 2 5 group outsiderock2 gen zon brick p0 6 0 0 p1 21 0 0 p2 6 1 0 p3 6 0 6 size 10 2 4 group outsiderock3 gen 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 outsiderock5 gen zon reflect dip 90 dd 90 orig 0 0 0 range x -0.1 6.1 z 6.1 13.1gen zon reflect dip 90 dd 90 orig 0 0 0 range x -0.1 6.1 z -6.1 -12.1gen zon reflect dip 90 dd 90 orig 0 0 0 range x 6.1 21.1 z -12.1 13.1;绘制模型图plot block groupplot add axes red;plot set rotation 0 0 45 用于显示三维模型;设置重力set gravity 0 0 -10;给定边界条件fix z range z -12.01,-11.99fix x range x -21.01,-20.99fix x range x 20.99,21.01fix y range y -0.01 0.01fix y range y 0.99,1.01;求解自重应力场model mohrini density 1800 ;围岩的密度prop bulk=1.47e8 shear=5.6e7 fric=20 coh=5.0e4 tension=1.0e4 ;体积、剪切、摩擦角、凝聚力、抗拉强度set mech ratio=1e-4solvesave Gravsol.savplot cont zdisp outl onplot cont szz;毛洞开挖计算initial xdisp=0 ydisp=0 zdisp=0model null range group insiderock any group concretliner anyplot block groupplot add axes redset mech ratio=5e-4solvesave Kaiwsol.savplot cont zdispplot cont sdispplot cont szzplot cont xzz;模筑衬砌计算model elas range group concretliner anyplot block groupplot add axes redini density 2500 range group concretliner any ;衬砌混凝土的密度prop bulk=16.67e9,shear=12.5e9 range group concretliner any ;衬砌混凝土的体积弹模、剪切弹模set mech ratio=1e-4solvesave zhihusol.savplot cont zdispplot cont sdispplot cont szzplot cont xzz;完成计算分析第四章FISH语言及建模技术4-1数组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)endloopendafillAshow4-2函数操作newdef xxxaa = 2 * 3xxx = aa + bbendnewdef stress_sumstress_sum = 0.0pnt = zone headloop while pnt # nullStress_sum = stress_sum + z_sxx(pnt)pnt=z_next(pnt)end_loopendnewdef stress_sumsum = 0.0pnt = 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) + facend4-3函数删除与重定义newdef joeii=out(‘This is A Function’)enddef fredjoeendfreddef joeii=(‘This is A New Function’)endfred4-4字符串newdef in_defxx = in(msg+’(’+’default:’+string(default)+’):’)if type(xx) = 3in_def = defaultelsein_def = xxendifend;def moduli_datadefault = 1.0e9msg=’Input Young‘s modulus ’Y_mod = in_defdefault = 0.25msg=’Input Poisson‘s ratio ’p_ratio = in_defif p_ratio = 0.5 thenii = out(’ Bulk mod is undefined at Poisson‘s ratio = 0.5’)ii = out(’ Select a different val ue --’)p_ratio = in_defendifs_mod = y_mod / (2.0 * (1.0 + p_ratio))b_mod = y_mod / (3.0 * (1.0 - 2.0 * 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 shear4-5马蹄形隧道网格newgen 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,100gen 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,0gen zone reflect dip 90 dd 270 origin 0,0,0plot set rotation 30 0 30plot block group4-6复杂形状网格生成new; 主隧道网格gen zon radcyl p0 15 0 0 p1 23 0 0 p2 15 50 0 p3 15 0 8 &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; 辅助隧道网格gen point id 1 (2.969848,0.0,-0.575736)gen point id 2 (2.969848,50.0,-0.575736)gen zon radcyl p0 0 0 -1 p1 7 0 0 p2 0 50 -1 p3 0 0 8 p4 7 50 0 &p5 0 50 8 p6 7 0 8 p7 7 50 8 p8 point 1 p10 point 2 & 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 & p5 7 50 0 p6 7 0 -8 p7 7 50 -8 p9 point 1 p11 point 2 & size 3 10 6 4 dim 3 3 3 3 rat 1 1 1 1;衬砌网格sel shell range cyl end1 0 0 -1 end2 0 50 -1 rad 3; 隧道外围边界网格gen zone radtun p0 7 0 0 p1 50 0 0 p2 7 50 0 p3 15 0 50 p4 50 50 0 & p5 15 50 50 p6 50 0 50 p7 50 50 50 &p8 23 0 0 p9 7 0 8 p10 23 50 0 p11 7 50 8 &size 6 10 3 10 rat 1 1 1 1.1gen zone brick p0 0 0 8 p1 7 0 8 p2 0 50 8 p3 0 0 50 &p4 7 50 8 p5 0 50 50 p6 15 0 50 p7 15 50 50 &size 3 10 10 rat 1 1 1.1;对称得出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;对称得出完整模型gen zon reflect dip 90 dd 270 ori 0 0 0;网格显示plot set rotation 30 0 30plot block group4-7网格连接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 1.9 2.1model elasprop bulk 8e9 shear 5e9fix z range z -.1 .1fix x range x -.1 .1fix x range x 3.9 4.1fix y range y -.1 .1fix y range y 3.9 4.1apply szz -1e6 range z 3.9 4.1 x 0,2 y 0,2hist unbal;网格显示plot set rotation 30 0 30plot block groupsolvesave att.savplot cont zdisp outl on单一网格（小）程序newgen zone brick size 8 8 8 p0 0,0,0 p1 4,0,0 p2 0,4,0 p3 0,0,4 model elasprop bulk 8e9 shear 5e9fix z range z -.1 .1fix x range x -.1 .1fix x range x 3.9 4.1fix y range y -.1 .1fix y range y 3.9 4.1apply szz -1e6 range z 3.9 4.1 x 0,2 y 0,2hist unbal;网格显示plot set rotation 30 0 30plot block groupsolvesave noattx.savplot cont zdisp outl on单一网格（大）程序newgen zone brick size 4 4 4 p0 0,0,0 p1 4,0,0 p2 0,4,0 p3 0,0,4 model elasprop bulk 8e9 shear 5e9fix z range z -.1 .1fix x range x -.1 .1fix x range x 3.9 4.1fix y range y -.1 .1fix y range y 3.9 4.1apply szz -1e6 range z 3.9 4.1 x 0,2 y 0,2hist unbal;网格显示plot set rotation 30 0 30plot block groupsolvesave noattx.savplot cont zdisp outl on4-8立方体洞穴网格生成newdef parmrad=4.0len=10.0in_size=6rad_size=10endparmgen zone radbrick edge len size in_size in_size in_size rad_size &rat 1.0 1.0 1.0 1.2 dim rad rad rad;对称得出完整模型gen zon reflect dip 0 dd 0 ori 0 0 0gen zon reflect dip 90 dd 270 ori 0 0 0;网格显示plot set rotation 30 0 30plot block group4-9球体洞穴网格生成newdef parmrad=4.0len=10.0in_size=6rad_size=10endparmdef make_spherep_gp=gp_headloop while p_gp#nullpx=gp_xpos(p_gp)py=gp_ypos(p_gp)pz=gp_zpos(p_gp)dist=sqrt(px*px+py*py+pz*pz)if dist>0 thenk=rad/distax=px*kay=py*kaz=pz*kmaxp=max(px,max(py,pz))k=len/maxpbx=px*kby=py*kbz=pz*ku=(maxp-rad)/(len-rad)gp_xpos(p_gp)=ax+u*(bx-ax)gp_ypos(p_gp)=ay+u*(by-ay)gp_zpos(p_gp)=az+u*(bz-az)end_ifp_gp=gp_next(p_gp)end_loopendgen zone radbrick edge len size in_size in_size in_size rad_size & rat 1.0 1.0 1.0 1.2 dim rad rad radmake_sphere;对称得出完整模型gen zon reflect dip 0 dd 0 ori 0 0 0gen zon reflect dip 90 dd 270 ori 0 0 0;网格显示plot set rotation 15 0 30plot block group4-10应力边界newgen zone brick size (4,4,4) p0 (0,0,0) p1 (4,0,0) &p2 (0,4,0) p3 (2,0,3.464)model elasticprop bulk 1e8 shear .3e8apply nstress -1e6 range plane dip 60 dd 270 origin 0.1,0,0 above step 1;网格显示plot set rotation 30 0 30plot block groupplot grid red fap green4-11改变应力边界条件newgen zone brick size 6 6 6model elasprop bulk 1e8 shear 7e7fix x range x -0.1 0.1def supersteploop ns (1,n_steps)x_stress = stress_inccommandapply sxx add x_stress range x 5.9,6.1 y 0,6 z 0,2step 100end_commandend_loopendset n_steps=100 stress_inc=-1e3hist zone sxx 6,0,0plot create sxx_histplot showSuperstep;网格显示plot set rotation 30 0 30plot block groupplot grid red fap greennewgen zone brick size 6 6 6model elasprop bulk 1e8 shear 7e7fix x range x -0.1 0.1def x_stressx_stress = stress_inc * stependset stress_inc = -1e3apply sxx 1.0 hist x_stress range x 5.9,6.1 y 0,6 z 0,2 hist zone sxx 6,0,0hist x_stressstep 100;网格显示plot set rotation 30 0 30plot grid red fap greennewgen zone brick size 6 6 6model elasprop bulk 1e8 shear 7e7fix x range x -0.1 0.1table 1 0,0 100,-1e5apply sxx 1.0 hist table 1 range x 5.9,6.1 y 0,6 z 0,2 hist zone sxx 6,0,0step 100;网格显示plot set rotation 30 0 30plot grid red fap greennewgen zone brick size 6,6,6 p1 6,0,-1model elasprop bulk 8e9 shear 5e9apply sxx -2e6 range x -0.1 0.1apply sxx -2e6 range x 5.9 6.1step 500plot grid dispnewgen zone brick size 5,5,5model elasprop bulk 8e9 shear 5e9set grav 0 0 -10ini dens 1000fix x range x -.1 .1fix x range x 4.9 5.1fix y range y -.1 .1fix y range y 4.9 5.1ini szz -5e4 grad 0 0 -1e4app szz -5e4 range z -.1 .1solvemodel null range x 1,4 y 1 4 z 3 5step 100plot set plane dip 90 dd 180 origin 0,2.5,0plot add surf plane behind yellplot add vel plane behind blackplot set rotation 30 0 30plot block groupplot show4-12位移边界newgen zone brick size 4,4,4 p3 2,0,3.464model elasprop bulk 1e8 shear .3e8macro left_boun ’plane dip 60 dd 270 origin 0.1,0,0 above’macro right_boun ’plane dip 60 dd 270 origin 3.9,0,0 below’apply nvel 0.1 plane dip 60 dd 270 range left_bounapply nvel 0.1 plane dip 120 dd 90 range right_bounstep 1;网格显示plot set rotation 30 0 30plot block groupplot grid velnewnewgen zone brick size 4 4 4model elasticprop bulk 1e8 shear .3e8apply nstress -1e6 plane dip 0 dd 0 range y 3.9 4.1apply nvel 0.0 plane dip 0 dd 0 range y -.1 .1apply nvel 0.0 plane dip 0 dd 0 range x -.1 .1step 10;网格显示plot set rotation 30 0 30plot block groupplot grid velnewgen zone brick size 10 5 5mod elprop shear 1e8 bulk 2e8fix x y z range x -.1 .1 y 0 5 z 0 5fix x y z range x 0 10 y 0 5 z -.1 .1fix x y z range x 9.9 10.1 y 0 5 z 0 5table 1def find_addhead = nullp_gp = gp_headloop while p_gp # nullx_pos = gp_xpos(p_gp)if x_pos = width thennew = get_mem(2)mem(new) = headmem(new+1) = p_gphead = newendifp_gp = gp_next(p_gp)endloopendset width=10.0find_adddef apply_velwhile_steppingad = headloop while ad # nullp_gp = mem(ad+1)gp_xvel(p_gp) = vel_max * gp_zpos(p_gp) / heightgp_zvel(p_gp) = -vel_max * (gp_xpos(p_gp) - width) / height ad = mem(ad)endloopendset large vel_max=1e-2 height=5.0step 100;网格显示plot set rotation 30 0 30plot block groupplot grid vel4-13不考虑重力影响的均匀应力new ;set grav 0 0 0gen zone brick size 6 6 6model elasini sxx=-5e6 syy=-1e7 szz=-2e7apply sxx=-5e6 range x -0.1 0.1apply sxx=-5e6 range x 5.9 6.1apply syy=-1e7 range y -0.1 0.1apply syy=-1e7 range y 5.9 6.1apply szz=-2e7 range z -0.1 0.1apply szz=-2e7 range z 5.9 6.1prop bulk 8e9 shear 5e9ini dens 2000fix x range z -.1 .1step 10plot set rotation 30 0 30 ;网格三维显示plot cont zdisp outl on4-14考虑应力梯度的均匀材料newgen zone brick size 10 10 10 p1 20,0,0 p2 0,20,0 p3 0,0,20 model mohrprop bulk 5e9 shear 3e9 fric 35ini density 2500set gravity 0,0,-10fix x y z range z -0.1 0.1ini szz -5.0e6 grad 0,0,2.5e4ini syy -2.5e6 grad 0,0,1.25e4ini sxx -2.5e6 grad 0,0,1.25e4apply szz -4.5e6 range z 19.9 20.1apply szz -5.0e6 range z -0.1 0.1apply sxx -2.5e6 grad 0,0,1.25e4 range x -0.1 0.1apply sxx -2.5e6 grad 0,0,1.25e4 range x 19.9 20.1apply syy -2.5e6 grad 0,0,1.25e4 range y -0.1 0.1apply syy -2.5e6 grad 0,0,1.25e4 range y 19.9 20.1step 10;网格三维显示plot set rotation 30 0 30plot cont zdisp outl on4-15考虑应力梯度的非均匀材料newgen zone brick size 10 10 10 &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 -.1 .1fix x range x 19.9 20.1fix z range z -.1 .1fix z range z 19.9 20.1fix y range y -25.1 -24.9ini syy 0.0 grad 0,1.6e4,0 range y -10,0ini syy 4e4 grad 0,2.0e4,0 range y -15,-10ini syy 7e4 grad 0,2.2e4,0 range y -15,-25step 10;网格三维显示plot set rotation 30 0 30plot cont zdisp outl on4-16非均匀网格应力初始化newgen zone radcyl size 3 8 4 5 fill p1 10,0,0 p2 0,10,0 p3 0,0,10 mode elasticprop shear 3e8 bulk 5e8fix x range x -.1 .1fix x range x 9.9 10.1fix y range y -.1 .1fix y range y 9.9 10.1fix z range z -.1 .1ini szz = -2.5e5 grad 0,0,2.5e4ini density 2500set grav 0,0,-10step 10;网格三维显示plot set rotation 30 0 30plot cont zdisp outl on4-17不规则自由面应力初始化gen zone brick size 15 15 10 p0 0,0,0 edge=100.0model elasticprop shear 3e8 bulk 5e8def mountaingp = gp_headloop while gp # nullzz = sqrt(gp_xpos(gp)ˆ2 + gp_ypos(gp)ˆ2)dz = 0.06 * sin(0.2 * zz + 100.0) ; Sum Fourier terms fordz = dz + 0.06 * sin(0.22 * zz - 20.3) ; quasi-random surfacedz = dz - 0.04 * sin(0.33 * zz + 33.3) ; topology.gp_zpos(gp) = 0.5 * gp_zpos(gp) * (1.0 + dz)gp = gp_next(gp)end_loopendmountainfix x range x -.1 .1fix x range x 99.9 100.1fix y range y -.1 .1fix y range y 99.9 100.1fix z range z -.1 .1set grav 0,0,-10ini density=2000ini szz=-2.0e6 (grad 0,0,2.0e4) sxx=-4.0e6 (grad 0,0,4.0e4) syy=-4.0e6 (grad 0,0,4.0e4) step 100plot set rotation 30 0 30plot cont zdisp outl on;网格三维显示4-18非均网格内部压实newgen zone brick size 8 8 10 ratio 1.2 1 1model mohrini dens 2000prop bulk 2e8 shear 1e8prop fric 30fix x range x -.1 .1fix x range x 7.9 8.1fix y range y -.1 .1fix y range y 7.9 8.1fix z range z -.1 .1;ini szz -2.0e5 grad 0,0,2e4;ini sxx -1.5e5 grad 0,0,1.5e4;ini syy -1.5e5 grad 0,0,1.5e4set grav 10step 1000;pause;prop tens 1e10 coh 1e10;step 750;prop tens 0 coh 0;step 250;网格三维显示plot set rotation 30 0 30plot cont zdisp outl on4-19模型改变后初始应力变化newgen zone brick size 5 5 5model elasprop sh 2e8 bu 3e8fix x y z range z -.1 .1set grav 0 0 -10ini dens 2000Solve;网格三维显示plot set rotation 30 0 30plot cont zdisp outl on;模型改变model mohr range x 0 2 y 0 5 z 0 2prop sh 2e8 bu 3e8 fric 35 range x 0 2 y 0 5 z 0 2 Step 2000;网格三维显示plot set rotation 30 0 30plot cont zdisp outl on4-20应力与孔隙压力的初始化newconfig fluidgen zone brick size 8 5 10model elasmodel fl_isoini dens 2000prop bulk 1e9 shear 5e8prop poros 0.5 perm 1e-10ini fmod 2e9ini fdensity 1e3ini sat 0ini sat 1 range z -.1 5.1set grav 0 0 -10fix x range x -.1 .1fix x range x 7.9 8.1fix y range y -.1 .1fix y range y 4.9 5.1fix z range z -.1 .1ini pp 5.e4 grad 0,0,-1.e4 range z 0.0 5.ini szz -20e4 grad 0,0,20e3ini szz add -1.5e4 grad 0,0,.25e4 range z 5,6 ini szz add -2.5e4 grad 0,0,.5e4 range z 0,5 ini szz add -.25e4 range z 0,5solve;网格三维显示plot set rotation 30 0 30plot cont zdisp outl onnewgen zone brick size 1 1 10model elasini dens 2500 range z 0 5ini dens 2250 range z 5 6ini dens 2000 range z 6 10prop bulk 1e9 shear 5e8set grav 10water dens 1000water table ori 0 0 5 normal 0 0 1fix x range x -.1 .1fix x range x 7.9 8.1fix y range y -.1 .1fix y range y 4.9 5.1fix z range z -.1 .1ini szz -20e4 grad 0,0,20e3ini szz add -1.5e4 grad 0,0,.25e4 range z 5,6 ini szz add -2.5e4 grad 0,0,.5e4 range z 0,5 ini szz add -.25e4 range z 0,5solve;网格三维显示plot set rotation 30 0 30plot cont zdisp outl on4-21加载顺序new; 主隧道建模gen zon radcyl p0 15 0 0 p1 23 0 0 p2 15 50 0 p3 15 0 8 &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;辅助隧道建模gen point id 1 (2.969848,0.0,-0.575736)gen point id 2 (2.969848,50.0,-0.575736)gen zon radcyl p0 0 0 -1 p1 7 0 0 p2 0 50 -1 p3 0 0 8 p4 7 50 0 &p5 0 50 8 p6 7 0 8 p7 7 50 8 p8 point 1 p10 point 2 &size 3 10 6 4 dim 3 3 3 3 rat 1 1 1 1 fillgen zon radcyl p0 0 0 -1 p1 0 0 -8 p2 0 50 -1 p3 7 0 0 p4 0 50 -8 &p5 7 50 0 p6 7 0 -8 p7 7 50 -8 p9 point 1 p11 point 2 &size 3 10 6 4 dim 3 3 3 3 rat 1 1 1 1 fill;隧道边界网格gen zone radtun p0 7 0 0 p1 50 0 0 p2 7 50 0 p3 15 0 50 p4 50 50 0 & p5 15 50 50 p6 50 0 50 p7 50 50 50 &p8 23 0 0 p9 7 0 8 p10 23 50 0 p11 7 50 8 &size 6 10 3 10 rat 1 1 1 1.1gen zone brick p0 0 0 8 p1 7 0 8 p2 0 50 8 p3 0 0 50 &p4 7 50 8 p5 0 50 50 p6 15 0 50 p7 15 50 50 &size 3 10 10 rat 1 1 1.1gen 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 50group main1 range cyl end1 15 0 0 end2 15 25 0 rad 4group main2 range cyl end1 15 25 0 end2 15 50 0 rad 4group service range cyl end1 0.0,0.0,-0.575736 &end2 0.0,50.0,-0.575736 rad 3.0save tun0.savrest tun0.sav; 初始应力场模量model mohrprop shear 0.36e9 bulk 0.6e9 coh 1e5 fric 20 tens 1e5apply szz -1.4e6 range z 49.9 50.1fix z range z -50.1 -49.1fix x range x -.1 .1fix x range x 49.9 50.1fix y range y -.1 .1fix y range y 49.9 50.1ini sxx -1.4e6 syy -1.4e6 szz -1.4e6hist unbalhist gp xdis 3,0,-1hist gp zdis 0,0,2hist gp xdis 3,25,-1hist gp zdis 0,25,2step 1000save tun1.sav; 辅助隧道开挖25 m计算rest tun1.savini xdis 0.0 ydis 0.0 zdis 0.0model null range group service y 0,25step 1000save tun2.sav; 在辅助隧道上施加衬砌restore tun2.savsel shell id=1 range cyl end1 0 0 -1 end2 0 25 -1 rad 3sel shell prop iso=(25.3e9, 0.266) thick = 0.5sel node fix y xr zr range y -0.1 0.1 ; symmetry cond.sel node fix x yr zr range x -0.1 0.1 ; symmetry cond.model mohr range group service y 25 50prop shear 0.36e9 bulk 0.6e9 coh 1e5 fric 20 tens 1e5; 主隧道开挖25mmodel null range group main1step 1000save tun3.savplot set rotation 0 0 60plot cont zdisp outl on第六章双线铁路隧道施工过程分析6-2 V级围岩施工过程模拟（1）建立隧道1/4圆周模型gen zon radcyl p0 0 0 0 p1 10.55 0 0 p2 0 1 0 p3 0 0 10.55 &size 5 2 10 4 dim 5.55 5.55 5.55 5.55 rat 1 1 1 1.2 group outsiderocksyplot block groupplot add axes redgen zone cshell p0 0 0 0 p1 5.55 0 0 p2 0 1 0 p3 0 0 5.55 &size 1 2 10 4 dim 5.05 5.05 5.05 5.05 rat 1 1 1 1 group concretlinersy fill group insiderocksy gen zon radcyl p0 0 0 0 p1 0 0 -9 p2 0 1 0 p3 10.55 0 0 &size 5 2 10 4 dim 4 5.55 4 5.55 rat 1 1 1 1.2 group outsiderockxygen zone cshell p0 0 0 0 p1 0 0 -4 p2 0 1 0 p3 5.55 0 0 &size 1 2 10 4 dim 3.5 5.05 3.5 5.05 rat 1 1 1 1 group concretlinerxy fill group insiderockxy plot block groupplot add axes redgen zon brick p0 0 0 10.55 p1 10.55 0 10.55 p2 0 1 10.55 p3 0 0 25.55 &size 5 2 8 group outsiderock1gen zon brick p0 0 0 -34.55 p1 10.55 0 -34.55 p2 0 1 -34.55 p3 0 0 -9 &size 5 2 12 group outsiderock2gen zon brick p0 10.55 0 0 p1 50 0 0 p2 10.55 1 0 p3 10.55 0 10.55 &size 20 2 5 group outsiderock3gen zon brick p0 10.55 0 -9 p1 50 0 -9 p2 10.55 1 -9 p3 10.55 0 0 &size 20 2 5 group outsiderock3gen zon brick p0 10.55 0 10.55 p1 50 0 10.55 p2 10.55 1 10.55 p3 10.55 0 25.55 & size 20 2 8 group outsiderock4gen zon brick p0 10.55 0 -34.55 p1 50 0 -34.55 p2 10.55 1 -34.55 p3 10.55 0 -9 & size 20 2 12 group outsiderock5gen zon reflect dip 90 dd 90 orig 0 0 0plot set rotation 30 0 30 ;显示三维图，绕x轴负方向转30度set gravity 0 0 -10 ;设置重力加速度为z方向-10fix z range z -34.56,-34.54 ;设置底边界fix x range x -50.01,-49.99 ;设置左边界fix x range x 49.99,50.01 ;设置右边界fix y range y -0.01 0.01 ;设置前边界fix y range y 0.99,1.01 ;设置后边界model mohr ;莫尔～库仑模型ini density 2000 ;围岩的密度prop bulk=7.14e8 shear=3.333e9 fric=25 coh=2e5 tension=1.0e5Setp 6000 ;求解6000次Solve ;求解计算save Gsol.sav ;计算结果保存在Gsol.sav文件中plot cont zdisp ;绘制竖向位移场，如图6-14所示plot cont xdisp ;绘制竖向位移场，如图6-15所示plot cont szz ;绘制竖向应力场，如图6-16所示plot cont sxx ;绘制竖向应力场，如图6-17所示;1-左上半断面开挖rest Gsol.savplot block groupplot add axes redini xdis 0.0 ydis 0.0 zdis 0.0ini density 2200 range group outsiderocksy any group outsiderocksz anyprop bulk=1.923e9, shear=4.167e9 fric=35, coh=0.5e6 &tension=5e5, range group outsiderocksy any group outsiderocksz anymodel null range group insiderocksz any group concretlinersz anystep 2000save stepp1.savplot cont zdisprest stepp1.savsel shell id=1 range x -0.01 0.01 z -0.01 5.56 cyl end1 0 0 0 end2 0 1 0 rad 5.55sel shell id=1 range x -0.01 0.01 z -0.01 5.56sel shell id=1 prop iso=(21.0e9,0.2) thick=0.25step 2000plot cont zdispsave stepp2.sav;3-左下半断面开挖rest stepp2.savplot block groupplot add axes redmodel null range group insiderockxz any group concretlinerxz anystep 2000save stepp3.savplot cont zdisprest stepp3.savsel shell id=1 range x -0.01 0.01 z -4.01 0.01sel shell id=1 prop iso=(21.0e9,0.2) thick=0.25step 2000plot cont zdispsave stepp4.sav;5-右上半断面开挖rest stepp4.savplot block groupplot add axes redmodel null range group insiderocksy any group concretlinersy anystep 2000save stepp5.savplot cont zdisprest stepp5.savsel shell id=1 range x -0.01 5.56 z -0.01 5.56 cyl end1 0 0 0 end2 0 1 0 rad 5.55 sel shell id=1 prop iso=(21.0e9,0.2) thick=0.25step 2000plot cont zdispsave stepp6.sav;7-右下半断面开挖rest stepp6.savplot block groupplot add axes redmodel null range group insiderockxy any group concretlinerxy anystep 2000save stepp7.savplot cont zdisprest stepp7.savsel delete shell id=1 range x -0.01 0.01sel shell id=1 prop iso=(21.0e9,0.2) thick=0.25step 2000save stepp8.savrest stepp8.savmodel 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=12.5e9, shear=16.667e9 range group concretlinerxyprop bulk=12.5e9, shear=16.667e9 range group concretlinerxzprop bulk=12.5e9, shear=16.667e9 range group concretlinersyprop bulk=12.5e9, shear=16.667e9 range group concretlinerszstep 2000plot cont zdispsave stepp9.savplot block groupplot add axes red6.3 IV级围岩施工过程模拟new;1/4 上半部分gen zon radcyl p0 0 0 0 p1 10.55 0 0 p2 0 1 0 p3 0 0 10.55 size 5 2 10 4&dim 5.55 5.55 5.55 5.55 rat 1 1 1 1.2 group outsiderocksygen zone cshell p0 0 0 0 p1 5.55 0 0 p2 0 1 0 p3 0 0 5.55 size 1 2 10 4&dim 5.15 5.15 5.15 5.15 rat 1 1 1 1 group concretlinersy fill group insiderocksy;1/4 下半部分gen zon radcyl p0 0 0 0 p1 0 0 -9 p2 0 1 0 p3 10.55 0 0 &size 5 2 10 4 dim 4 5.55 4 5.55 rat 1 1 1 1.2 group outsiderockxygen zone cshell p0 0 0 0 p1 0 0 -4 p2 0 1 0 p3 5.55 0 0 size 1 2 10 4&dim 3.6 5.15 3.6 5.15 rat 1 1 1 1 group concretlinerxy fill group insiderockxy;1/2绘制上下地层网格gen zon brick p0 0 0 10.55 p1 10.55 0 10.55 p2 0 1 10.55 p3 0 0 25.55 size 5 2 8 group outsiderock1gen zon brick p0 0 0 -34.55 p1 10.55 0 -34.55 p2 0 1 -34.55 p3 0 0 -9 size 5 2 12 group outsiderock2;1/2右侧中间部分土体网格gen zon brick p0 10.55 0 0 p1 50 0 0 p2 10.55 1 0 p3 10.55 0 10.55 size 20 2 5 group outsiderock3gen zon brick p0 10.55 0 -9 p1 50 0 -9 p2 10.55 1 -9 p3 10.55 0 0 size 20 2 5 group outsiderock3。
隧道及地下工程FLAC解析方法 命令流
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
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
对全部高中资料试卷电气设备，在安装过程中以及安装结束后进行高中资料试卷调整试验；通电检查所有设备高中资料电试力卷保相护互装作置用调与试相技互术关，系电，力根保通据护过生高管产中线工资敷艺料设高试技中卷术资配0料不置试仅技卷可术要以是求解指，决机对吊组电顶在气层进设配行备置继进不电行规保空范护载高高与中中带资资负料料荷试试下卷卷高问总中题体资，配料而置试且时卷可，调保需控障要试各在验类最；管大对路限设习度备题内进到来行位确调。保整在机使管组其路高在敷中正设资常过料工程试况中卷下，安与要全过加，度强并工看且作护尽下关可都于能可管地以路缩正高小常中故工资障作料高；试中对卷资于连料继接试电管卷保口破护处坏进理范行高围整中，核资或对料者定试对值卷某，弯些审扁异核度常与固高校定中对盒资图位料纸置试，.卷保编工护写况层复进防杂行腐设自跨备动接与处地装理线置，弯高尤曲中其半资要径料避标试免高卷错等调误，试高要方中求案资技，料术编试交写5、卷底重电保。要气护管设设装线备备置敷4高、调动设中电试作技资气高，术料课中并中3试、件资且包卷管中料拒含试路调试绝线验敷试卷动槽方设技作、案技术，管以术来架及避等系免多统不项启必方动要式方高，案中为；资解对料决整试高套卷中启突语动然文过停电程机气中。课高因件中此中资，管料电壁试力薄卷高、电中接气资口设料不备试严进卷等行保问调护题试装，工置合作调理并试利且技用进术管行，线过要敷关求设运电技行力术高保。中护线资装缆料置敷试做设卷到原技准则术确：指灵在导活分。。线对对盒于于处调差，试动当过保不程护同中装电高置压中高回资中路料资交试料叉卷试时技卷，术调应问试采题技用，术金作是属为指隔调发板试电进人机行员一隔，变开需压处要器理在组；事在同前发一掌生线握内槽图部内纸故，资障强料时电、，回设需路备要须制进同造行时厂外切家部断出电习具源题高高电中中源资资，料料线试试缆卷卷敷试切设验除完报从毕告而，与采要相用进关高行技中检术资查资料和料试检，卷测并主处且要理了保。解护现装场置设。备高中资料试卷布置情况与有关高中资料试卷电气系统接线等情况，然后根据规范与规程规定，制定设备调试高中资料试卷方案。
地层学原理---第五章 地层层序 地层结构 地层
18
（五） 岩石地层结构的研究方法
1. 基本层序调查内容： ① 单层成分； ② 鉴别特殊成分和成因的夹层（生物层、粘土层、含矿层、风化壳等）； ③ 单层厚度、形态、岩石结构、沉积构造、遗迹化石、古生态、古流向等； ④ 基本层序内单层间的叠覆特征和接触关系，是否具旋回性和趋向性； ⑤ 基本层序的纵横向变化； ⑥ 与理想的相模式进行对比，对形成环境和沉积作用进行分析。
H
10
H
11
评述
根据基准面变化的相对振幅来划分沉积层序级别，能够 减少主观性、易于在实践中应用。主要依据：层序边界范围， 不整合面分布的范围，层序界面下伏地层的构造变形程度与 上覆地层的海侵范围，沉积系统与构造背景的改变程度。
H
12
二、地层结构（stratigraphic structure）
2. 李勇等（1994）提出3种岩石地层结构：进积型、加积型和 退积型，陆相地层均决于冲积体系进积、加积、退积的 演化过程。
3. 魏家庸等（1991）提出韵律结构（Rhythmic texture）、旋 回 结 构 （ Cyclic textare ） 、 均 一 结 构 （ Homogeneous texture）、加积结构（Aggradational textune）、退积结 构 （ Retrogradational texture ） 和 进 积 结 构 （Progradational texture）
H
19
1. 基本层序调查方法：在实测主干剖面是研究基本层序的主 要方法，辅之以主干地质路线和辅助地质路线的调查。
2. 应注意的问题：
① 相模式与基本层序的区别：相模式是根据研究古代岩石、现 代沉积物和现代沉积环境及作用的观测与实验所获的大量 资料，对特定沉积环境和沉积作用或其产物普遍特征的概 括和简化的表达形式，常见的类型有直观的环境模式和理 想化的作用——产物模式（自旋回性的垂向相模式）。因 此，相模式一般性的、代表环境和作用本质的理想化标准， 基本层序是具体的，包含偶然性、地方性特色的描述实体。
工程地质分析原理-名词解释
工程地质条件：与工程建筑有关的地质条件的总称。
包括岩土体工程性质、地质构造、地形地貌、水文地质特征、自然地质现象、工程动力地质作用和天然建筑材料等七方面。
工程地质问题：工程地质条件与工程建筑之间所存在的矛盾（稳定和安全）。
因工程建筑结构类型、受力特点和工作方式不同，存在各种工程地质问题。
工程地质分析的基本研究方法：1.自然历史分析法2.力学分析法3.模型模拟试验法4.工程地质类比法断层分类：走滑断层、逆断层、正断层，按断裂的主次关系划分为主断层、分支断层和次断层野外识别：走滑断层：1地表出露线平直，断层倾角较陡。
2断层面两侧相对的水平运动，相对的垂直升降很小。
3河流最易于沿这种断层发育，水工建筑物也就最易于受到这种活断层的威胁。
走滑断层：位移方向与断层走向相同或相近的大型平移断层活断层：一般理解为目前还在持续活动的断层，或在历史时期或近期地质时期活动过、极可能在不远的将来重新活动的断层。
活动断层的活动特性:蠕滑:持续不断缓慢蠕动的称为蠕滑或稳滑 ;粘滑:间断地、周期性突然错断的为粘滑，粘滑常伴有地震，是活断层的主要活动方式。
各类活断层在我国的分布1．分布于我国西部的多个NWW—NW向的长大弧形断裂，主要为逆走滑—走滑断层，其走滑速率由南而北有所减小。
2．活动逆冲断裂也主要分布于我国西部分、3．现代活动正断层主要分布于我国东部。
活断层区规划设计建筑物的原则(1)有低级别的活断层的场地优于有高级别的; 有活动时期老的断层的场地优于有活动时期新的; 有全新世(10000 a)内无活动的断层的场地优于有全新世内有活动的断层的场地等。
(2)尽可能避开主断层带；(3)如为逆断层或正断层类型，尽可能避开有强烈地表变形和分支、次生断裂发育的断层上盘(逆断层的上升盘、正断层下降盘)。
如有较大的正、逆断层，场地往往需要选在距主断面数千米之外。
地震：接近地球表面的岩层中弹性波传播所引起的震动称为地震。
震级：衡量地震本身大小的尺度烈度：衡量震动强烈程度的尺度基本烈度：一个地区今后一定时期内在一般场地条件下可能普遍遭受到的最大地震烈度。
命令流隧道及地下工程解析方法李围
第六章命令流按照顺序进行6-1 V级围岩施工过程模拟(1)建立隧道1/4圆周模型gen zon radcyl p0 0 0 0 p1 10.55 0 0 p2 0 1 0 p3 0 0 10.55 &size 5 2 10 4 dim 5.55 5.55 5.55 5.55 rat 1 1 1 1.2 group outsiderocksyplot block groupplot add axes redgen zone cshell p0 0 0 0 p1 5.55 0 0 p2 0 1 0 p3 0 0 5.55 &size 1 2 10 4 dim 5.05 5.05 5.05 5.05 rat 1 1 1 1 group concretlinersy fill group insiderocksy gen zon radcyl p0 0 0 0 p1 0 0 -9 p2 0 1 0 p3 10.55 0 0 &size 5 2 10 4 dim 4 5.55 4 5.55 rat 1 1 1 1.2 group outsiderockxygen zone cshell p0 0 0 0 p1 0 0 -4 p2 0 1 0 p3 5.55 0 0 &size 1 2 10 4 dim 3.5 5.05 3.5 5.05 rat 1 1 1 1 group concretlinerxy fill group insiderockxy plot block groupplot add axes redgen zon brick p0 0 0 10.55 p1 10.55 0 10.55 p2 0 1 10.55 p3 0 0 25.55 &size 5 2 8 group outsiderock1gen zon brick p0 0 0 -34.55 p1 10.55 0 -34.55 p2 0 1 -34.55 p3 0 0 -9 &size 5 2 12 group outsiderock2gen zon brick p0 10.55 0 0 p1 50 0 0 p2 10.55 1 0 p3 10.55 0 10.55 &size 20 2 5 group outsiderock3gen zon brick p0 10.55 0 -9 p1 50 0 -9 p2 10.55 1 -9 p3 10.55 0 0 &size 20 2 5 group outsiderock3gen zon brick p0 10.55 0 10.55 p1 50 0 10.55 p2 10.55 1 10.55 p3 10.55 0 25.55 &size 20 2 8 group outsiderock4gen zon brick p0 10.55 0 -34.55 p1 50 0 -34.55 p2 10.55 1 -34.55 p3 10.55 0 -9 &size 20 2 12 group outsiderock5gen zon reflect dip 90 dd 90 orig 0 0 0plot set rotation 30 0 30 。
地层结构法命令流
finish/clear/prep7!-----------------------------------------------------------------------------------------------------------------------------!隧道轮廓k,1001,,,k,1002,0,2.217k,1,0,10.258k,2,-5.832,6.629larc,1,2,1001,6.5 !创建圆弧，第一端点，第二端点，圆心所在一侧任一关键点，半径k,3,-6.5,3.758larc,2,3,1001,6.5k,4,-6.429,2.629larc,3,4,1001,9k,5,-6.367,2.217larc,4,5,1001,9k,6,-6.179,1.134l,5,6k,7,-5.617,1.134l,6,7k,8,-5.174,0.874larc,7,8,1001,1.802larc,8,1001,1002,15.75lcomb,2,3 !合并线! LCOMB, NL1, NL2, KEEP--------其中：NL1,NL2---拟合并的两条线号。
NL1可为ALL，或元件名。
! KEEP---是否保留输入的线及其公共关键点控制参数。
! KEEP=0则删除NL1和NL2及其公共关键点，如果已经划分网格则不能删除，或者依附于其它图素也不能删除! KEEP=1则保留线及其公共关键点，但公共关键点不依附于新创建的线。
!该命令可以合并独立线或依附于同面上的线，合并后便于网格划分。
!可合并的线可为直线或曲线，以及直线与曲线，可共线或不共线。
!当为多条时，应为多条首尾相连的线。
!无论在何种坐标系下执行合并，合并后的线不改变合并前的空间位置。
lcomb,4,5lcomb,7,8lsymm,x,all !对称，镜像!LSYMM, Ncomp, NL1, NL2, NINC, KINC, NOELEM, IMOVE!Ncomp---对称控制选项，可选X（缺省），Y，Z值。
08 地层结构模式计算方法
山岭隧道
19
洞室围岩应力分析
散元法、颗粒流法等，其中有限元和边界元建立在连续介 质力学的基础上，适合于小变形分析，是发展较早、较成 熟的方法，尤以有限元在地下工程中应用更为广泛。有限 元法不但可以模拟各种施工过程和各种支护效果，同时可 以考虑复杂地层情况（如断层、节理等地质构造以及地下 水等），对材料的各种非线性也可考虑。
山岭隧道
s
o
u
t
h
w
e
s
t
j
I
a
o
t
o
n
g
u
n
I
v
e
r
s
I
t
y
《山岭隧道》课程
第五章
支护结构设计
郑余朝
西南交通大学
Southwest Jiaotong University
1
主要内容
荷载结构模式计算方法 地层结构模式计算方法 复合式衬砌结构设计 单层衬砌结构设计 TBM管片衬砌结构设计
（6.1.9）
t
由边界条件
r
r
r r1
m1 Ec mEC C1 B1 1 c 1 2 B1 2 m1 1 (m1 1)r r
（6.1.10）
C1 B1 2 p w r1
C1 B1 2 p a r0
可得：
pa p w 2 2 C1 2 r0 r1 2 r0 r1
荷载-结构模型和地层结构模型计算实例：北京地铁十号线八达岭高速站明挖、暗挖结构设计
1.3 结构统一技术要求1.3.1.地下铁道结构中主要构件的设计使用年限为100年，主要构件是指结构的主体结构（梁、板、墙、柱）及基础结构。
相应结构可靠度理论的设计基准期均采用50年。
取1.1；1.3.2.地下铁道结构中永久构件的安全等级为一级，相应的结构构件重要性系数γ取0.9；在人防荷载或地震荷载临时构件的安全等级为三级，相应的结构构件重要性系数γ组合下，相应的结构构件重要性系数γ取1.0。
1.3.3. 地下铁道结构的地震作用按8度设防，地下结构框架柱的抗震等级暂按三级，梁、板、墙等构件暂按四级。
待“北京地铁十号线场地地震安全评估报告”完成后，依其为设计依据。
1.3.4.地下结构中露天或与无侵蚀性的水或土壤直接接触的迎土面混凝土构件的环境类别为二类，非迎土面及内部混凝土构件的环境类别为一类，两者均视为一般环境条件。
1.3.5. 结构构件在永久荷载和基本荷载组合作用下，应按荷载短期效应组合并考虑长期效应组合的影响进行结构构件裂缝验算。
二类环境混凝土构件的裂缝宽度（迎土面）应不大于0.2mm，一类环境（非迎土面及内部混凝土构件）混凝土构件的裂缝宽度均应不大于0.3mm，混凝土管片内外侧的裂缝宽度应不大于0.2mm。
当计及地震、人防或其他偶然荷载作用时，可不验算结构的裂缝宽度。
1.3.6. 地下铁道结构中主要构件的耐火等级为一级。
1.3.7. 在规定的设防部位，地下结构按5级人防的抗力标准进行验算。
1.3.8.当地下结构处于有侵蚀地段时，应采取抗侵蚀措施，混凝土抗侵蚀系数不得低于0.8。
2.工程概况八达岭高速站结构设计形式为两端明挖、中间暗挖车站，其中两端明挖部分为地下双层三跨岛式车站，车站中部过健德桥部分采用分离暗挖单洞通过，车站主体、风道基坑围护结构均采用钻孔灌注桩，基坑内设置横向钢支撑，出入口明挖部分采用型钢围檩，内设钢支撑，车站主体、风道、出入口的明挖部分采用现浇钢筋混凝土箱形框架结构，东北、西北两出入口跨路部分采用暗挖法施工。
8-隧道地层结构计算解读
第八章隧道支护地层—结构分析方法（杨林德、朱合华、丁文其）第一节一般规定8.1.1 设计模型1、与地层—结构分析方法相应的隧道设计模型是地层—结构模型。
其设计理念，是认为围岩具有自支承能力，支护（含衬砌结构，下同）的作用是加固围岩，并与围岩联合组成共同受力的整体，共同承受荷载的作用。
2、采用地层—结构分析方法进行设计计算时，计算区范围应同时包含支护和围岩地层。
3、采用地层—结构分析方法进行设计计算时，应同时考虑开挖施工步骤的影响。
初步设计阶段可按常规施工方法选定开挖施工步骤，施工图设计阶段应改按施工组织设计制定的技术方案确定。
4、采用地层—结构分析方法设计隧道时，应同时检验围岩的稳定性和支护结构的受力变形状态是否满足强度条件及按使用要求确定的变形量限制条件。
5、采用地层—结构分析方法设计隧道时，内衬结构的工作状态应为弹性受力状态，或经论证认为仍可保持稳定的弹塑性受力状态；初期支护（含开挖阶段增设的喷射混凝土层）和围岩的工作状态可为弹塑性受力状态。
8.1.2 地层—结构分析方法的适用地质条件1、地层—结构分析方法的基础理念，是认为围岩具有自支承能力，并可由其与支护结构共同组成承载体系。
因而这类方法适用于在具有一定自支承能力的围岩中建造的隧道支护结构的计算。
2、 V 级及 V 级以上的围岩都具有一定的自支承能力，因而都可采用地层—结构分析方法进行设计计算。
3、Ⅲ级及Ⅲ级以上的围岩自支承能力强，对在这些级别的围岩中建造的隧道，经验表明对支护结构根据经验选定设计参数时已可使围岩保持稳定，因而规范规定一般不要求进行计算。
4、 IV 级、 V 级围岩中建造的隧道一般采用复合式支护，对其宜采用地层—结构分析方法进行设计计算。
但对在V 级围岩中建造的浅埋隧道，围岩承载能力较低时仍宜采用荷载- 结构分析方法计算。
5、 IV 级围岩的自支承能力优于V 级围岩，采用地层—结构分析方法进行设计计算时，宜通过控制荷载释放过程，使IV 级围岩中隧道内衬结构经受的荷载相对较小，围岩的自支承能力可适度充分发挥。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
finish/clear/prep7!-----------------------------------------------------------------------------------------------------------------------------!隧道轮廓k,1001,,,k,1002,0,2.217k,1,0,10.258k,2,-5.832,6.629larc,1,2,1001,6.5 !创建圆弧，第一端点，第二端点，圆心所在一侧任一关键点，半径k,3,-6.5,3.758larc,2,3,1001,6.5k,4,-6.429,2.629larc,3,4,1001,9k,5,-6.367,2.217larc,4,5,1001,9k,6,-6.179,1.134l,5,6k,7,-5.617,1.134l,6,7k,8,-5.174,0.874larc,7,8,1001,1.802larc,8,1001,1002,15.75lcomb,2,3 !合并线! LCOMB, NL1, NL2, KEEP--------其中：NL1,NL2---拟合并的两条线号。
NL1可为ALL，或元件名。
! KEEP---是否保留输入的线及其公共关键点控制参数。
! KEEP=0则删除NL1和NL2及其公共关键点，如果已经划分网格则不能删除，或者依附于其它图素也不能删除! KEEP=1则保留线及其公共关键点，但公共关键点不依附于新创建的线。
!该命令可以合并独立线或依附于同面上的线，合并后便于网格划分。
!可合并的线可为直线或曲线，以及直线与曲线，可共线或不共线。
!当为多条时，应为多条首尾相连的线。
!无论在何种坐标系下执行合并，合并后的线不改变合并前的空间位置。
lcomb,4,5lcomb,7,8lsymm,x,all !对称，镜像!LSYMM, Ncomp, NL1, NL2, NINC, KINC, NOELEM, IMOVE!Ncomp---对称控制选项，可选X（缺省），Y，Z值。
!NL1,NL2,NINC---按增量NINC从NL1到NL2定义关键点的范围（缺省为NL1），NINC缺省为1。
NL1也可为ALL或元件名，此时NP2和NINC将被忽略。
!KINC---要创建的关键点编号增量，缺省时由系统自动指定（不会覆盖）。
!NOELEM---是否创建单元和节点控制参数。
NOELEM=0（缺省）如果存在单元和节点则生成；NOELEM=1不生成单元和节点。
!IMOVE---线是否被移动或重新创建。
IMOVE=0（缺省）原来的线不动，重新创建新线；当IMOVE=1不创建新线，原来的线移动到新位!置，此时编号不变（即ITIME、KINC和NOELEM均无效），且单元和节点一并移动。
!该命令要求当前坐标系为直角坐标系，线可以在任意象限。
同KSYMM相同，可通过设定当前坐标系为工作平面或局部坐标系而改变镜像位置。
nummrg,all !粘接!nummrg是一个使用率极高的命令，使用起来非常的方便。
它通常用来实现“粘接”的效果，即，合并坐标完全一致的两个图素，具体做法就是!删除编号大的图素，并用编号小的图素占据原来的位置。
这样操作以后，会实现在不同图素之间形成共用部分，从而在一定程度上取代Vglue命令。
!但是，对于复杂模型，在进行图素合并命令之前，对操作结果的预期不明朗之时，还是尽量使用vglue为好，因为这样能最大限度地保障有限元模型的连续性。
numcmp,all !压缩编号alls!-------------------------------------------------------------------------------------------------------------------------------------! 建立地层范围blc4,-40,-40,80,90!命令：BLC4, XCORNER, YCORNER, WIDTH, HEIGHT, DEPTH 通过一角点坐标和尺寸创建矩形面!XCORNER,YCORNER---矩形面或块体第1个角点在工作平面上的X和Y坐标。
!WIDTH---平行于工作平面X轴方向离XCORNER的距离。
!HEIGHT---平行于工作平面Y轴方向离YCORNER的距离!DEPTH---离工作平面的垂直距离，即平行于Z轴。
!DEPTH=0（缺省）则生成面。
如WIDTH或HEIGHT或DEPTH为负值，则为反方向距离。
asbl,all,all !从面中减去线nummrg,allnumcmp,allalls!-----------------------------------------------------------------------------------------------------------------------------!细分隧道分析模型wprota,,,90! 工作平面的旋转!命令：WPROTA, THXY, THYZ, THZX!其中THXY, THYZ, THZX为绕工作平面坐标系Z轴、X轴和Y轴的旋转角度。
asbw,all !从面中减去与工作面相交的部分wpoffs,,,20!将工作平面沿其自身坐标轴移动!命令：WPOFFS, XOFF, YOFF, ZOFF!其中XOFF, YOFF, ZOFF为工作平面坐标系内沿其X轴、Y轴和Z轴的偏移增量。
asbw,all !用工作面切分面wpoffs,,,-40asbw,allwpoffs,,2.139wprota,,90asbw,allwpoffs,,,-4asbw,allwpoffs,,,-20asbw,allwpoffs,,,46asbw,allwpoffs,,,-20.995wprota,,90wpoffs,,,-13.821asel,s,,,24,29,5asbw,allwpoffs,,,-0.562asel,s,,,10,17,7asel,a,,,13asbw,allwpoffs,,,-11.234asel,s,,,23,30,7asel,a,,,14asbw,allwpoffs,,,-0.562asel,s,,,35,37,2asbw,allwprota,,90asel,s,,,18,38,20asel,a,,,27,28,1asbw,alllcomb,1,8lcomb,3,14lcomb,54,37lcomb,57,39wpcsys,-1 !使工作平面的位置与当前坐标系的X-Y 平面或R-θ平面一致，但不改变视图方向。
! WPCSYS，WN，KCN!WN 是窗口号，该窗口中的视图方向修改为与工作平面垂直。
默认值是1。
如果WN 为负值，则视图方向不变。
!KCN 是坐标系号，默认为当前激活的坐标系。
如果KCN 是直角坐标系，工作平面与其X-Y 平面一致；如果是圆柱或球坐标系，则与R-θ平面一致。
allsnummrg,allnumcmp,all!------------------------------------------------------------------------------------------------------------------------------------!设置实体单元类型ET,1,PLANE42KEYOPT,1,1,0 !参见Elements Reference部分的Element Library中的plane42 KEYOPT,1,2,0KEYOPT,1,3,2 !设置为平面应变模式KEYOPT,1,5,0KEYOPT,1,6,0mp,ex,1,1.5e9 !!地层!!加固地层mp--定义线性材料属性mp,prxy,1,0.4mp,dens,1,2000tb,dp,1 !激活非线性材料的数据表tbdata,1,0.15e6,25 !定义tb数据表中的数据et,2,beam3*set,th1,0.25r,2,th1,th1*th1*th1/12,th1mp,ex,2,2.555e10 !!初期支护mp,prxy,2,0.2mp,dens,2,2500mp,ex,3,1.5e9 !!上台阶内部土体mp,prxy,3,0.4mp,dens,3,2000tb,dp,3tbdata,1,0.15e6,25mp,ex,4,1.5e9 !!下台阶内部土体mp,prxy,4,0.4mp,dens,4,2000tb,dp,4tbdata,1,0.15e6,25mp,ex,5,1.5e9 !!仰拱内部土体mp,prxy,5,0.4mp,dens,5,2000tb,dp,5tbdata,1,0.15e6,25!----------------------------------------------------------------------------------------------------------------------------------------------!!!初次划分单元网格lsel,s,,,1,3,2lsel,a,,,11,13,2lsel,a,,,4lesize,all,,,12 !确定划分指定线的网格数和网格尺寸alls!LESIZE, NL1, SIZE, ANGSIZ, NDIV, SPACE, KFORC, LAYER1, LAYER2, KYNDIV!NL1---Number of the line to be modified. If ALL, modify all selected lines [LSEL]. If NL1 = P, graphical! picking is enabled and all remaining command fields are ignored (valid only in the GUI).A component name may also be substituted for NL1.!SIZE---If NDIV is blank, SIZE is the division (element edge) length. The number of divisions is automatically calculated from the line length (rounded upward to! next integer). If SIZE is zero (or blank), use ANGSIZ or NDIV.!ANGSIZ---The division arc (in degrees) spanned by the element edge (except for straight lines, which always result in one division). The number of divisions is! automatically calculated from the line length (rounded upward to next integer).!NDIV---If positive, NDIV is the number of element divisions per line. If -1 (and KFORC = 1), NDIV is assumed to be zero element divisions per line. TARGE169 with! a rigid specification ignores NDIV and will always mesh with one element division.!SPACE---Spacing ratio. If positive, nominal ratio of last division size to first division size (if > 1.0, sizes increase, if < 1.0, sizes decrease). If negative,! |SPACE| is nominal ratio of center division(s) size to end divisions size. Ratio defaults to 1.0 (uniform spacing). For layer-meshing, a value of 1.0! normally is used. If SPACE = FREE, ratio is determined by other considerations!KFORC---KFORC 0-3 are used only with NL1 = ALL. Specifies which selected lines are to be modified.! 0 —Modify only selected lines having undefined (zero) divisions.! 1 —Modify all selected lines.! 2 —Modify only selected lines having fewer divisions (including zero) than specified with this command.! 3 —Modify only selected lines having more divisions than specified with this command.! 4 —Modify only nonzero settings for SIZE, ANGSIZ, NDIV, SPACE, LAYER1, and LAYER2. If KFORC = 4, blank or 0 settings remain unchanged.!LAYER1---Layer-meshing control parameter. Distance which defines the thickness of the inner mesh layer, LAYER1. Elements in this layer are uniformly-sized with! edge lengths equal to the specified element size for the line (either through SIZE orline-length/NDIV). A positive value for LAYER1 is interpreted as! an absolute length, while a negative value in interpreted as a multiplier on the specified element size for the line. As a general rule, the resulting! thickness of the inner mesh layer should be greater than or equal to the specified element size for the line. If LAYER1 = OFF, layer-meshing control! settings are cleared for the selected lines. The default value is 0.0!LAYER2---Layer-meshing control parameter. Distance which defines the thickness of the outer mesh layer, LAYER2. Elements in this layer transition in size from! those in LAYER1 to the global element size. A positive value of LAYER2 is interpreted as an absolute length, while a negative value is interpreted as a! mesh transition factor. A value of LAYER2 = -2 would indicate that elements should approximately double in size as the mesh progresses normal to LAYER1.! The default value is 0.0.!KYNDIV---! 0, No, and Off means that SmartSizing cannot override specified divisions and spacing ratios. Mapped mesh fails if divisions do not match. This defines! the specification as "hard".! 1, Yes, and On means that SmartSizing can override specified divisions and spacing ratios for curvature or proximity. Mapped meshing can override! divisions to obtain required matching divisions. This defines the specification as" soft".lsel,s,,,35,37,2lsel,a,,,48,51,1lsel,a,,,52lesize,all,,,8allslsel,s,,,36,70,34lsel,a,,,7,9,2lsel,a,,,71,76,5lsel,a,,,74,80,6lesize,all,,,11allslsel,s,,,33,65,32lsel,a,,,6,8,2lsel,a,,,38,69,31lsel,a,,,72,83,11lesize,all,,,1allslsel,s,,,39,41,1 lsel,a,,,25,32,7 lsel,a,,,87,91,4 lsel,a,,,86,90,4 lesize,all,,,2 allslsel,s,,,55,56,1 lsel,a,,,53,54,1 lsel,a,,,58,60,2 lsel,a,,,61 lsel,a,,,42,43,1 lsel,a,,,77,88,11 lsel,a,,,27,81,54 lsel,a,,,15,64,49 lsel,a,,,21,24,3 lesize,all,,,6 allslsel,s,,,12,19,7 lsel,a,,,17,30,13 lsel,a,,,2,5,3 lsel,a,,,34 lsel,a,,,16,31,15 lsel,a,,,57,59,2 lsel,a,,,44,47,3 lsel,a,,,26,29,3 lsel,a,,,85,92,7 lsel,a,,,89,93,4 lsel,a,,,66,68,1 lsel,a,,,22,75,53 lsel,a,,,78,84,6 lsel,a,,,23,62,39 lsel,a,,,10,18,8 lsel,a,,,45,46,1 lsel,a,,,28,73,45 lsel,a,,,79,82,3 lsel,a,,,14,20,6 lsel,a,,,63 lesize,all,,,6 alls!----------------------------------------------------------------------------------------------------------------------------------------!定义单元的材料属性及单元范围，并最终划分单元网格type,1 !选取单元类型mat,1 !单元材料属性amap,24,31,36,1,26 !简化面映射网格划分，amap,area,kp1,kp2,kp3,kp4----其中，area---拟划分的面号，kp1、kp2、kp3、kp4---指定的角点，3或4个都可以，顺序可任意。