FLAC3D 5.0培训日程安排培训日程安排1.FLAC3DV5.0界面操作界面操作2.FLAC3D基本操作方法基本操作方法vs应用流程应用流程；；3.FLAC3D内置内置Fish语言的应用；语言的应用；4.FLAC3D结构单元结构单元vs接触单元；接触单元；5.FLAC3D渗流模块渗流模块6.其他其他1 StructuralElement2 FLAC3D结构单元1.结构单元的类型2.结构单元的建模方法3.结构单元的参数取值4.结构单元实例分析5.关于link3 FLAC3D中包含六种形式的结构单元，可以分成两类：中包含六种形式的结构单元，可以分成两类：Ø线型结构单元线型结构单元：：•梁梁单元单元(beam)•锚索单元锚索单元(cable)•桩单元桩单元(pile)Ø壳壳型型结构单元结构单元：：•壳单元壳单元(shell)•土工格栅土工格栅(geogrid)•衬砌单元衬砌单元(liner)ØFLAC3D中的结构单元是岩土工程中实际结构的一种“抽象”，即采用简单的单元形式来模拟复杂的结构体Ø结构单元由结构节点(node)和结构构件(SELs)构成Ø结构单元中的节点(node)可以与周围的实体网格(zone)或其它结构节点建立连接(link)，通过连接实现岩土体或结构与其它结构发生相互作用。
Ø注意：结构节点并不是简单地与实体网格的节点(gridpoint)建立联系，也不能建立node与gridpoint之间的link1、结构单元的类型4 Ø梁梁单元单元sel beam id 1 beg 4 0 -1 end 5 0 -2 nseg 4sel node id=1 0 0 0sel node id=2 2 0 0 sel node id=3 4 0 -1sel node id=4 5 0 -2sel beamsel cid=1 id=1 node 1 2 ;sel beamsel cid=2 id=1 node 2 3sel beamsel cid=3 id=1 node 3 4Ø桩桩单元单元sel pile id 1 beg 0 0 0 end 0 0 10 nseg 42、结构单元的建模方法两种建模方式各有各的优点，第二种方式适合建立复杂曲线结构单元（但但是要注意它不会自是要注意它不会自动建立动建立link！！若不！！若不手动手动link就无任何作就无任何作用用）5 Ø锚索锚索单元单元sel cable id 1 beg 4 0 -1 end 5 0 -2 nseg 42、结构单元的建模方法6 建立梁单元，并显示单元坐标系！2、结构单元的建模方法—线型结构单元起始点坐标并给定分段数目的方法；7 ID号相同，共用Node，ID不同，各个ID对应的结构单元有各自独立的node。
除非设置联系，否则即使节点位于同一位置也不会传递力结构单元的显示！GUI操作和命令操作（manual）！调整好显示效果后可以将显示的命令文件另存出来，以备下次使用最适用于几何模型相同，参数不同的，不同工况分析的比较）2、结构单元的建模方法—线型结构单元8 先建立节点再联接成单元的方法；2、结构单元的建模方法—线型结构单元9 Ø壳单元壳单元2、结构单元的建模方法—壳型结构单元10 def set_vals global ptA = 25.0 * sin( 40.0*degrad ) ; global ptB = 25.0 * cos( 40.0*degrad )end@set_valsgenerate zone cylinder p0=( 0.0, 0.0, 0.0 ) & p1=( @ptA, 0.0, @ptB ) & p2=( 0.0, 25.0, 0.0 ) & p3=( 0.0, 0.0, 25.0 ) & p4=( @ptA, 25.0, @ptB ) & p5=( 0.0, 25.0, 25.0 ) & size=(1, 2, 2)sel shell id=5 range cylinder end1=(0.0, 0.0,0.0) & end2=(0.0,25.0,0.0) radius=24.5 notplot add zg plot ad sel geom delete zones ; delete all zonessel node init zpos add -25.02、结构单元的建模方法—壳型结构单元11 通过附着在实体网格表面来生成shell单元。
The shells can then be repositioned if ecessary by using the SEL node init command2、结构单元的建模方法—壳型结构单元12 ØFLAC3D是岩土工程的专业软件，因此一般很少用来做专门的结构是岩土工程的专业软件，因此一般很少用来做专门的结构分析在涉及到结构单元的问题中，往往都要考虑结构与周围的实分析在涉及到结构单元的问题中，往往都要考虑结构与周围的实体单元的相互作用在结构单元的建模时要特别注意一个基本原则：体单元的相互作用在结构单元的建模时要特别注意一个基本原则：一个一个zone至多包含一个至多包含一个structure node!Ø因此在建立线型结构单元时，要特别注意因此在建立线型结构单元时，要特别注意nseg变量的大小变量的大小nseg太太小则会导致计算不精确，而太大就会违反结构单元建模的基本原则小则会导致计算不精确，而太大就会违反结构单元建模的基本原则2、结构单元的建模方法—注意事项13 Ø梁单元梁单元 •emod——弹性模量，弹性模量，E•nu——泊松比，泊松比，ν•xcarea——横截面积，横截面积，A•xciy——梁结构梁结构y轴惯性矩轴惯性矩, Iy•xciz——梁结构梁结构z轴惯性矩，轴惯性矩，Ix•xcij——极惯性矩，极惯性矩，J•density——密度，密度，ρ•pmoment——塑性矩，塑性矩，Mp•thexp——热膨胀系数，热膨胀系数，αt•ydirection——矢量矢量YØ锚索锚索单元单元 •emod——弹性模量，弹性模量， E•xcarea——横截面积，横截面积，A•gr_coh——单位长度上水泥浆粘结力单位长度上水泥浆粘结力cg•gr_fric——水泥浆的摩擦角水泥浆的摩擦角φg•gr_k——单位长度上单位长度上水泥浆刚度水泥浆刚度kg•gr_per——水泥浆外圈周长水泥浆外圈周长Pg•slide——大变形滑动标志大变形滑动标志•slide_tol——大变形滑动容差大变形滑动容差•ycomp——抗压强度抗压强度(力力)•density——密度密度•thexp——热膨胀系数热膨胀系数3、结构单元的参数取值14 3、结构单元的参数取值某些结构单元参数的取值要视具体情况而定，根据经验且必要时调整参数通过试算来确定。
15 4、结构单元实例分析4.1、简支梁（beam单元）承受两个相等集中载荷4.2、简支梁（shell单元）承受两个相等集中载荷16 Simple Beam – Two Equal Concentrated Loads4.1、简支梁（beam单元）承受两个相等集中载荷17 A simply supported beam is loaded by two equal concentrated loads, symmetrically placed as shown in Figure 1.9. The shear and moment diagrams for this configuration are also shown in the figure.The shear force magnitude,V, is equal to the applied concentrated load,P. The maximum moment,Mmax, occurs between the two loads and is equal to Pa. The maximum deflection of the beam,max, occurs at the center and is given by AISC (1980, p. 2-116) as4.1、简支梁（beam单元）承受两个相等集中载荷18 载荷(N)载荷距支座的距离(m)铰支座之间的距离(m)弹性模量(Pa)惯性矩(m^4)惯性矩(m^5)最大挠度(m)PaLEIyIzΔmax10000392.00E+11 2.00E-04 2.00E-04 0.006468750 根据理论公式计算得到：4.1、简支梁（beam单元）承受两个相等集中载荷19 newtitle "Simple Beam --- Two Equal Concentrated Loads Symmetrically Placed"; ======================================================; Create the grid, insure that nodes will exist at third points.sel beam id=1 begin=( 0, 0, 0) end=( 3, 0, 0) nseg=3sel beam id=1 begin=( 3, 0, 0) end=( 6, 0, 0) nseg=4sel beam id=1 begin=( 6, 0, 0) end=( 9, 0, 0) nseg=3; ======================================================; Assign beam propertiessel beam id=1 prop emod=2e11 nu=0.30 &xcarea=6e-3 xcj=0.0 xciy=200e-6 xciz=200e-6;========================================================; Specify model boundary conditions (including applied loads)sel node fix z xr yr ; restrict all non-beam modessel node fix y range id=1 ;sel node fix y range id=9 ; ; rollers at beam endssel node apply force=(0.0,-1e4,0.0) range id=2 ; apply point loadssel node apply force=(0.0,-1e4,0.0) range id=5 ;4.1、简支梁（beam单元）承受两个相等集中载荷20 ; ==========================================================; Setup histories for monitoring behavior.history add id=10 sel node ydisp id=7history add id=30 sel beamsel moment mz end2 cid=1 ; moment, right of SEL-1history add id=31 sel beamsel moment mz end1 cid=2 ; moment, left of SEL-2; ========================================================; Bring the problem to equilibriumsolve ratio=1e-7save equal-concent-loads; ========================================================; Print out beam responses.list sel beam forcelist sel beam momentlist sel node disp range id=7return4.1、简支梁（beam单元）承受两个相等集中载荷如何设置结构单元的跟踪变量！21 Beam_concent_loads_Example1.34.1、简支梁（beam单元）承受两个相等集中载荷22 挠度计算4.1、简支梁（beam单元）承受两个相等集中载荷23 剪力、弯矩计算24 剪力、弯矩计算这是节点力！！！25 26 梁单元局部坐标系：x轴从节点1到节点2，y轴在横截面中4.1、简支梁（beam单元）承受两个相等集中载荷27 弯矩矢量的指向，右手法则！弯矩矢量的指向，右手法则！4.1、简支梁（beam单元）承受两个相等集中载荷28 •List sel beamnodal forces: components are displayed in terms of the beam local coordinate systems. These are the forces exerted by the nodes on the beamSEL.小结：梁单元的常用命令4.1、简支梁（beam单元）承受两个相等集中载荷history sel beamselcid怎么找？（坐标or鼠标information？）29 selSelnode命令针对所有的结构单元30 Sel node fix keyword . . .31 newtitle "Simple Beam (modeled using shellSELs)"gen zone brick size 12,3,1 &p0 0,0,1.0 p1 9,0,1.0 p2 0 0 0 p3 0 1 1.0 ;shell 宽度为单位1（z方向上）sel shell id=1 crossdiag elemtype=dkt range y -0.1 0.1sel shell id=1 prop iso=(2e11, 0.0) thick=0.133887delete zonesel node fix x y xr yr range x=(-0.1, 0.1) ; support at left end – hinge 铰支座sel node fix y xr yr range x=( 8.9, 9.1) ; support at right end – roller 辊轴支座sel node fix z xr yr ; restrict non-beam deformation modessel node apply force=(0,-1667,0) range union id=71 id=12 ; out nds, leftsel node apply force=(0,-3333,0) range union id=46 id=13 ; in nds, leftsel node apply force=(0,-1667,0) range union id=79 id=24 ; out nds, rt.sel node apply force=(0,-3333,0) range union id=54 id=25 ; in nds, rt.将均布载荷转换为等效节点力4.2、简支梁（shell单元）承受两个相等集中载荷32 history add id=1 unbalhistory add id=10 sel node ydisp id=19 ; displ at center; moment, left thirdhistory add id=20 sel recover sres Mx surfX 1,0,0 cid=59; shear, left thirdhistory add id=30 sel recover sres Qx surfX 1,0,0 cid=59solve ratio=1e-7list sel node disp range id=19save shell0return 4.2、简支梁（shell单元）承受两个相等集中载荷33 4.2、简支梁（shell单元）承受两个相等集中载荷34 sel shell id=1 elemtype=cst range y -0.1 0.1 ;no-crossdiagsel shell id=1 crossdiag elemtype=dkt range y -0.1 0.1Crossdiag vs no-crossdiag 4.2、简支梁（shell单元）承受两个相等集中载荷35 局部坐标系！看弯矩到底应该看哪一个？？36 4.2、简支梁（shell单元）承受两个相等集中载荷37 4.2、简支梁（shell单元）承受两个相等集中载荷38 39 4.2、简支梁（shell单元）承受两个相等集中载荷40 surfX Xx Xy XzThe surfx vector (Xx, Xy, Xz) enables a surface coordinate system to be generated for all nodes used by the shell-type SELs in the optional range. The surface coordinate system, x'y'z', has the following properties: (1) z' is normal to the surface; (2) x' is the projection of the given surfx vector onto the surface; and (3) y' is orthogonal to x' and z'. The z'-direction is found at each node by taking the average normal direction of all shell-type SELs in the range. If the surfx vector is aligned at z' at any node, then processing stops and an error message is displayed. To proceed, designate a different surfx vector, or restrict the range of shell-type SELs considered.The surface coordinate system can be queried with the command LIST sel recover surface and the FISH function nd_ssys. It can also be set for an individual node with the FISH function nd_ssysx. It can be visualized with the sel geometry plot item by setting the systemtype switchword.The validity of the surface system at a particular node can be queried with the FISH function nd_svalid. The surface system at a node automatically becomes invalid under the following conditions: (1) large-strain update; or (2) creation or deletion of a shell-type SEL that uses the node. Validity must be reestablished with the SEL recover surface command.Sel Recover keyword 41 5、关于link两种：Node-ZoneNode-Node42 创建一个新的link，link的源节点为sid，而联接目标为node或zone。
Id为新link的ID号Sid是已经存在的节点（作为源节点）的ID号，可选关键字target用于确定目标对象（node或zone）默认的目标对象为zone对于zone目标对象，如果tid没有定义，将会使用与源节点距离在delta范围内的非空zone；否则，如果tid定义了，如果该tid所指示的单元为非空zone，且该zone的边界距离在delta之内，就会建立link对于node目标对象，tid就必须定义了，且两个节点必须彼此很靠近由delta确定如果不能确定源和目标对象，就会报错，且该命令不会对模型产生任何作用新link的attachment条件设置为6个自由度均为“rigid”Side1，side2关键字对于确定嵌入式liner的哪个面上产生linksel可选参数与必选参数！！可选参数与必选参数！！43 5、关于link以预应力锚杆的托盘模拟为例sel cable id=1 beg 0, 0, 0 end 0 ,29, 0 nseg 10 sel cable id=1 beg 0,29,0 end 0,35,0 nseg 6sel cable id=1 prop emod 2e10 ytension 310e3 xcarea 0.0004906 &gr_coh 1 gr_k 1 gr_per 0.0785 range cid 1,10sel cable id=1 prop emod 2e10 ytension 310e3 xcarea 0.0004906 &gr_coh 10e5 gr_k 2e7 range cid 11,17sel delete link range id 1 ；这里删除的是谁的id？Sel link id=100 1 target zonesel link attach xdir=rigid ydir=rigid zdir=rigid xrdir=rigid yrdir=rigid zrdir=rigid range id 100sel cable id=1 pretension 60e3 range cid 1,1044 FLAC3D数值模拟计算实例pBeam单元基坑开挖支护p以前述extrude的模型为例，进行隧道和基坑开挖支护45 bulk(Pa)shear(Pa)fric(ᵒ)coh(Pa)tens(Pa)1e80.3e8351e31e3开挖区域材料力学参数计算模型几何边界pBeam单元基坑开挖支护FLAC3D数值模拟计算实例数值模拟计算实例46 基坑垮塌过程pBeam单元基坑开挖支护FLAC3D数值模拟计算实例数值模拟计算实例47 基坑垮塌过程pBeam单元基坑开挖支护FLAC3D数值模拟计算实例数值模拟计算实例48 基坑垮塌过程pBeam单元基坑开挖支护FLAC3D数值模拟计算实例数值模拟计算实例49 基坑垮塌过程pBeam单元基坑开挖支护FLAC3D数值模拟计算实例数值模拟计算实例50 pBeam单元基坑开挖支护基坑垮塌过程FLAC3D数值模拟计算实例数值模拟计算实例51 pBeam单元基坑开挖支护基坑垮塌过程FLAC3D数值模拟计算实例数值模拟计算实例52 pBeam单元基坑开挖支护基坑垮塌过程FLAC3D数值模拟计算实例数值模拟计算实例53 pBeam单元基坑开挖支护基坑垮塌过程FLAC3D数值模拟计算实例数值模拟计算实例54 pBeam单元基坑开挖支护基坑垮塌过程FLAC3D数值模拟计算实例数值模拟计算实例55 pBeam单元基坑开挖支护FLAC3D数值模拟计算实例数值模拟计算实例56 pBeam单元基坑开挖支护最大不平衡力FLAC3D数值模拟计算实例数值模拟计算实例57 pBeam单元基坑开挖支护最大不平衡力FLAC3D数值模拟计算实例数值模拟计算实例58 Interface59 Interface-概述岩土工程中涉及到很多的接触问题，比如挡土墙与墙后填土之间的接触、桩与土接触、土石坝中混凝土防渗墙与土体之间的接触等。
FLAC/FLAC3D提供了接触面单元，可以分析一定受力条件下两个接触的表面上产生错动滑移、分开与闭合概述FLAC 和 FLAC3D中的接触面单元可以用来模拟岩体中的节理、断层； 地基与土体之间的接触； 矿仓与仓储物的接触面； 相互碰撞物体之间的接触面； 空间中的障碍边界（即固定的不变形的边界）条件FLAC3D中建立接触面单元应遵循以下原则： 小的表面与大的表面相连时，接触面应建立在小的表面上； 如果两相邻的网格有不同的密度，接触面应建立在密度大的区域上； 接触面单元尺寸通常应该等于或小于相连的目标面的尺寸； 使用 Attach 命令连接的两个表面不应再建立接触面60 Interface-概述FLAC 和 FLAC3D中的接触面采用的是无厚度接触面单元，接触面本构模型采用的是库仑剪切模型FLAC3D中接触面的基本理论FLAC3D中接触面单元由一系列三节点的三角形单元构成，接触面单元将三角形面积分配到各个节点中，每个接触面节点都有一个相关的表示面积每个四边形区域面用两个三角形接触面单元来定义，然后在每个接触面单元顶点上自动生成节点，当另外一个网格面与接触面单元相连时，接触面节点就会产生。
FLAC3D中接触面是单面的，认识这一点很重要，这点与二维 FLAC 中所定义的双面接触面不同可以把接触面看作“收缩带”，可以在指定面上拉伸，从而导致接触面和与之可能相连的其它任何面的相互刺入变得敏感接触面单元可以通过接触面结点和实体单元表面（称为目标面）之间来建立联系接触面法向方向所受到的力由目标面方位所决定在每个时间步计算中，首先得到接触面节点和目标面之间的绝对法向刺入量和相对剪切速度，再利用接触面本构模型来计算法向力和切向力的大小61 Interface-概述接触面单元、接触面节点以及节点表示面积的示意图62 为何要分离网格后再“移来移去”63 64 1、手册中给出的接触面建立方法——So-called “移来移去法”gen zone radcyl p0 (0,0,0) p1 (8,0,0) p2 (0,0,-5) p3 (0,8,0) &p4 (8,0,-5) p5 (0,8,-5) p6 (8,8,0) p7 (8,8,-5) &p8 (.3,0,0) p9 (0,.3,0) p10 (.3,0,-5) p11 (0,.3,-5) &size 3 10 6 15 ratio 1 1 1 1.15gen zone radcyl p0 (0,0,-5) p1 (8,0,-5) p2 (0,0,-8) p3 (0,8,-5) &p4 (8,0,-8) p5 (0,8,-8) p6 (8,8,-5) p7 (8,8,-8) &p8 (.3,0,-5) p9 (0,.3,-5) p10 (.3,0,-8) p11 (0,.3,-8) &size 3 6 6 15 ratio 1 1 1 1.15 fillgen zone reflect dd 270 dip 90group zone clay ;group clay (old command）65 1、手册中给出的接触面建立方法——So-called “移来移去法”interface 1 face range cylinder end1 (0,0,0) end2 (0,0,-5.1) radius .31 &cylinder end1 (0,0,0) end2 (0,0,-5.1) radius .29 notinterface 2 face range cylinder end1 (0,0,-4.9) end2 (0,0,-5.1) radius .31gen zone cylinder p0 (0,0,6) p1 (.3,0,6) p2 (0,0,1) p3 (0,.3,6) p4 (.3,0,1) p5 (0,.3,1) &size 3 10 6& ;写简写的时候要注意新老版本的区别gen zone cylinder p0 (0,0,6.1) p1 (.3,0,6.1) p2 (0,0,6) p3 (0,.3,6.1) &p4 (.3,0,6) p5 (0,.3,6) &size 3 1 6gen zone reflect dd 270 dip 90 range z 1 6.1group zone pile range z 1 6.1pauseini zposition add -6.0 range group pile ;ini z add -6.0 range group pile(old command）66 1、手册中给出的接触面建立方法——So-called “移来移去法”67 68 69 70 “移来移去”法; Create Basegen zone brick size 3 3 3 &p0 (0,0,0) p1 (3,0,0) p2 (0,3,0) p3 (0,0,1.5) &p4 (3,3,0) p5 (0,3,1.5) p6 (3,0,4.5) p7 (3,3,4.5)group zone Base; Create Top - 1 unit high for initial spacinggen zone brick size 3 3 3 &p0 (0,0,2.5) p1 (3,0,5.5) p2 (0,3,2.5) p3 (0,0,7) &p4 (3,3,5.5) p5 (0,3,7) p6 (3,0,7) p7 (3,3,7)group zone Top range group Base not;; Create interface elements on the top surface of the baseinterface 1 face range plane norm (-1,0,1) origin (1.5,1.5,3) dist 0.1; Lower top to complete geometryini zpos add -1.0 range group Top71 2、接触面建立方法——So-called “倒来倒去法”我们最终的目的就是在中心小块体与外围网格之间建立接触面。
1.分开建立网格2.建立inner网格及其表面的Interface3.导入外围mesh4.赋予材料属性，测试接触面是否发生了作用72 小练习：1.三种建立接触面的方法计算结果是否相同（只要接触面有响应，肯定是相同的！）2.如果将接触面建立在外部网格的内表面，然后移入小块体，结果是否相同呢？3.不加接触面跟加了接触面，模型的响应（位移、应力）有何区别？NrangeExample-7.1 union nrange73 3、切割模型的方法——实际上是分离连续网格（原来网格连续，通过共用节点（GridPoint）传递力，分离后通过接触面来传递gen zone brick size 3 3 3group zone inner range x 1 2 y 1 2 z 2 3group zone out range group inner notgenerate separate face group aa range group inner group out ;very different from old versioninterface 1 wrap first group inner ;second group out ;interface 1 permeability oninterface 1 maxedge 0.574 Wrap first keyword . . . second keyword . . .Interface elements are created on all zone faces belonging to the range specified after the keyword first. The tokens following first are a range descriptor as though normally used following a range keyword. Optionally, a second range may be given following the second keyword. Interface elements are created along the zone boundary between the first and second ranges. If no second range is given, the default second range is the entire model, indicating that the entire boundary of first range will be used. Note that for an interface element to be created, an exact match must exist between gridpoints in space on either side of the boundary, although the faces themselves do not have to match exactly. To separate one group from another, see the GENERATE separate command.For example, the following command would find the twinned faces between group "rock" and group "soil," and put interface elements on the "rock" faces. Only faces with centroids within the range x 50.0 75.0 would be considered.interface 1 wrap first group rock second group soil range x 50.0 75.0 接触面建立在第1个类组上面，若不指定第2个类组，就默认第2个类组为整个模型，边界面为整个第1个类组的表面。
75 Generate Separate ZonesSeparate Zonesseparate face keyword separates (unmerges) the internal faces specified by the range. The gridpoints of the face are duplicated, and a new surface face is created. New faces and gridpoints get copies of all group and extra variable assignments belonging to the original face and gridpoint.Note that faces may be restricted by giving two group range elements, therefore indicating that the face must be on the surface of the first and the second group. For instance, if "Fred" and "George" are group names assigned to zones, thenrange group Fred group George ;生成的应该是一个共享面will select faces that are connected to both Fred and George. Also note that faces can be selected by group directly.（此外，（此外，face也可以通过也可以通过group直接选中）直接选中）The following keywords can be used to affect the behavior产生由产生由range定义的内部定义的内部 face。
面上的节点复制，并创建新的表面面上的节点复制，并创建新的表面76 The following keywords can be used to affect the behavior.Clearattach By default, an error occurs if any gridpoint that has an attach condition associated with it is found among those to be separated. However, if the clearattach keyword is supplied, then the separation of gridpoints will occur regardless. In addition, FLAC3D will remove any attach conditions connected to gridpoints affected.Group name Newly created faces will be assigned the group name name in the specified slot. The default slot is 1. If add is specified, the name will be added to the first available slot. FLAC3D determines which face will be new and which one will be old by using the origin keyword.originx y zspecifies a location in space used to determine new versus old face assignments. If the vector from the origin to the face centroid is in the same direction as the outward face normal, it is the old face. If it is in the opposite direction, then it is new, and the group name specified is assigned. By default, the origin is (0,0,0).如果从origin指向面中心的向量与面的外法线方向相同，那么就作为old面。
77 检验是否分离开了网格并形成了接触面？？78 新版本中的Facegroup法; select as regional joining Base and Topgroup face Interface internal range group Base group Top; Separate Interface faces, calling the newwly created faces NewFacesgen separate face origin (0,0,0) group NewFaces range group Interface;; Create interface elements on the top surface of the baseinterface 1 face range group Interface79 80 实例——BinFlowSlip; Assign group names to faces that will become interfacesgroup face Int1 internal range group Material group Bin ...plane orig (0,0,0) normal (1,-1,0) abovegroup face Int2 internal range group Material group Bin ...plane orig (0,0,0) normal (1,-1,0) below; Separate those faces, giving newly created surface faces a new group namegen separate face origin (0,0,0) group NewFaces range group Int1 or Int2; Created interfaces on Int1 and Int2 facesinterface 1 face range group Int1interface 2 face range group Int2; Subdivide interface elements a little for better contact detectionint 1 maxedge 0.55int 2 maxedge 0.5581 1. FLAC3D基本操作方法82 1. FLAC3D基本操作方法p软件界面FLAC3D5.00 软件界面83 1. FLAC3D基本操作方法p软件界面FLAC3D5.00 软件界面84 1. FLAC3D基本操作方法p命令流操作方法newgen zone brick size 6 6 6在利用FLAC3D软件进行数值模拟时，主要是通过命令流来实现的。
命令流文件一般以txt或dat格式存储，并在命令输入窗口、菜单栏或快捷图标通过call命令进行调用FLAC3D命令流文件需要遵循一定的格式和语法要求，在满足这些要求的前提下，命令流文件的编写又是十分自由灵活的，可根据用户个人的需求自由编写85 1. FLAC3D基本操作方法p数值计算一般流程u建立模型通过外部导入或在FLAC3D中直接建模的方式建立计算模型u材料参数Ø材料本构Ø材料力学参数u边界条件及初始条件Ø速度边界Ø应力边界Ø水头边界86 1. FLAC3D基本操作方法p建立模型p在FLAC3D中，模型的建立通过关键词generate来实现，其基本格式为： gen keywords1 keywords2 keywords3 …For example: gen zone brick size 6 6 6 p0 0 0 0 p1 6 0 0 p2 0 6 0 … p3 0 0 6 p4 6 6 0 p5 0 6 6 … p6 6 0 6 p7 0 0 6 … (ratio 1 1 1 ) (dim 2 2 2) (fill)87 1. FLAC3D基本操作方法p建立模型FLAC3D建模的基本思路为“堆积木”，即首先建立各种形状的网格单元，最后将建立的网格单元组合在一起，生成可用于数值计算的整体模型。
FLAC3D内置13种不同形状的网格，包括brick(砖形), cshell(圆柱状壳网格), cylinder(圆柱状网格), cylint(圆柱状交叉网格), dbrick(退化砖形网格), pyramid(锥形网格), radbrick(砖形辐射网格), radcylinder(圆柱状辐射网格), radtunnel(平行六边形状辐射网格), retrahedron(辐射网格), tunint(砖形交叉网格), uwedge(均匀楔形网格), wedge(楔形网格).88 1. FLAC3D基本操作方法p网格形状89 1. FLAC3D基本操作方法基本操作方法p brickgen zone brick size 6 8 8plot zone90 1. FLAC3D基本操作方法基本操作方法p dbrickgen zone dbrick size 6 6 6 p0 0 0 0 p1 6 0 0 p2 0 6 0 … p3 0 0 6 p4 6 6 0 p5 0 6 6 p6 6 0 6plot zone91 1. FLAC3D基本操作方法基本操作方法p wedgegen zone wedge size 6 6 8 p0 0 0 0 p1 6 0 0 p2 0 6 0 p3 0 0 6plot zone92 1. FLAC3D基本操作方法基本操作方法p uwedgegen zone uwedge size 6 6 6 p0 0 0 0 p1 6 0 0 p2 0 6 0 p3 0 0 6plot zone93 1. FLAC3D基本操作方法基本操作方法p pyramidgen zone pyramid size 6 6 6 p0 0 0 0 p1 6 0 0 p2 0 6 0 … p3 0 0 6 p4 6 6 0plot zone94 1. FLAC3D基本操作方法基本操作方法p tetrahedrongen zone tetrahedron size 6 6 6 p0 0 0 0 p1 6 0 0 p2 0 6 0 … p3 0 0 6 plot zone95 1. FLAC3D基本操作方法基本操作方法p cylindergen zone cyl size 6 4 10 p0 0 0 0 p1 6 0 0 p2 0 6 0 p3 0 0 6plot zone96 1. FLAC3D基本操作方法基本操作方法p radbrickgen zone radbrick size 3 3 3 6 p0 0 0 0 p1 6 0 0 p2 0 6 0 … p3 0 0 6 dim 2 4 2 ratio 1 1 1 1.2 (fill)plot zone97 1. FLAC3D基本操作方法基本操作方法p radcylindergen zone radcylinder size 6 6 6 12 … p0 0 0 0 p1 6 0 0 p2 0 6 0 p3 0 0 6 dim 2 2 2 2 … ratio 1 1 1 1.2 (fill)plot zone98 1. FLAC3D基本操作方法基本操作方法p通过对称生成网格gen zone radbrick ...p0 (0,0,0) p1 (10,0,0) p2 (0,10,0) p3 (0,0,10) ...size 3,5,5,7 ...ratio 1,1,1,1.5 ...dim 1 4 2 fillgen zone reflect dip 0 dd 90gen zone reflect dip 90 dd 90plot zone99 1. FLAC3D基本操作方法基本操作方法p生成渐变网格gen zone brick size 10 10 10 p0 0 0 0 &p1 9 0 0 p2 0 9 0 p3 0 0 9 rat 0.8 0.8 0.8100 1. FLAC3D基本操作方法基本操作方法p分组(group)gen zone brick size 10 10 10group 1 range z 0 1group 2 range z 1 2 x 2 8save s1model null range group 2101 1. FLAC3D基本操作方法p材料参数 在FLAC3D中，对材料参数是通过两个命令关键词来实现的：MODLE和PROPERTY.FLAC3D中内嵌11种本构模型：Ønull model 空模型Øelastic, isotropic model 各向同性线弹性模型ØMohr-Coulomb model 摩尔库伦模型ØHoek-Brown model 霍克布朗模型ØDruck-Prager model 德鲁克普拉格模型null model常被用于开挖过程的模拟，elastic model赋予材料的线弹性力学行为，mohr model赋予材料塑性力学行为102 1. FLAC3D基本操作方法p材料参数 对于线弹性本构模型(elastic model), 需要指定的力学参数包括体积模量(bulk)和剪切模量(shear)，或者杨氏模量(young)和泊松比(poisson).对于摩尔库伦塑性模型(mohr model), 需要指定的力学参数包括：Ø体积模量和剪切模量，或杨氏模量和泊松比Ø内摩擦角和剪胀角Ø粘聚力Ø抗拉强度如果上述参数没有被指定，那么默认情况下该参数为0.103 1. FLAC3D基本操作方法p边界条件边界条件通常通过命令APPLY和FIX/FREE来施加，初始条件通过命令INITIAL来施加。
Ø在x=10的面上施加10MPa的压应力x分量apply sxx -10e6 range x 10Ø施加线性变化的应力apply syy = -20e6 grad 0,0,20e5 range y -20 z 0,10Y方向正应力按z坐标线性变化：z = 10时σyy = 0 ，z = 0时，σyy = −20 × 106 Pa.当使用关键词gradient 时，应满足如下关系式：S = S(o) + gxx + gyy + gzz104 1. FLAC3D基本操作方法p边界条件Ø法向速度约束 fix x range x 0Ø固定约束 fix x y z range x 0Ø施加初始应力状态 ini sxx -50e6 syy -40e6 szz -10e6Ø设置重力加速度set grav 0 0 -9.81105 2 FISH函数介绍106 2. Fish 函数pFish 函数结构 def function-name ;定义函数名 …… ;函数语句 end ;函数结束的标志 pFish函数举例 def func b=10 a=b+10 end @func print @a107 2. Fish 函数pFish 函数与变量对于Fish函数和变量需要注意一下几点：Ø函数和变量的赋值遵循数据类型的法则；Ø变量和函数名的命名规则不能以数字开头，不能含有中文，不能含有下列字符： . , * / + - ^ = < > # ( ) [ ] @ ; ‘ ’ “ ”Ø变量和函数名不能与FLAC3D, FISH保留字相冲突；Ø对变量赋值时，不能将当前函数的函数名放在“=”右边，这样会形成递归调用；Ø变量和函数名是全局的；Ø对FISH函数和变量的引用和通过如下命令来实现：•PRINT用于查看函数和变量的数值•HISTORY可对函数和变量的数值进行记录•SET可对变量进行赋值108 2. Fish 函数pFish 函数数据类型Ø整型：-2147483648~-2147483648Ø浮点型：10-300~10300Ø字符型：以(’)为分界符，常用于保存时文件名的定义Ø指针型：表示单元和节点的存储地址109 2. Fish 函数pFish 函数语句Ø选择语句 CASEOF 表达式 …默认语句 CASE n1 …表达式值为n1时的语句 CASE n2 …表达式值为n2时的语句 ENDCASEØ条件语句 IF 条件表达式 … ENDIF注：IF语句内部可继续嵌套条件语句Ø循环语句 LOOP var (exp1, exp2) … ENDLOOP或 LOOP WHILE 条件表达式 … ENDLOOPØ命令语句 COMMAND … ENDCOMMAND110 2. Fish 函数p单元遍历和节点遍历Ø基本变量 zone_head和gp_head:分别表示单元和节点的头指针 z_next()和gp_next:分别表示下一个单元（节点）Ø单元遍历程序框架： p_z=zone_head loop while p_z # null ;语句 p_z=z_next(p_z) endloopØ节点遍历程序框架： p_gp=gp_head loop while p_z # null ;语句 p_gp=z_next(p_gp) endloop111 3 FLAC3D数值模拟计算实例112 3.FLAC3D数值模拟计算实例p洞室模型的生成p基坑开挖过程模拟p路基施工过程模拟p荷载引起地基土体的超孔隙水压力模拟p边坡安全系数求解p隧道开挖及支护模拟p矿体开挖回填模拟p具有自由水面的稳定流113 p洞室模型的生成洞室模型的生成3.FLAC3D数值模拟计算实例数值模拟计算实例114 p洞室模型的生成洞室模型的生成3.FLAC3D数值模拟计算实例数值模拟计算实例gen zone radcylinder size 5 20 20 15 p0 0 0 615 p1 add 30 0 0 p2 add 0 50 0 ... p3 add 0 0 30 p4 add 30 50 0 p5 add 0 50 30 ... p6 add 30 0 30 p7 add 30 50 30 ... dim 10 10 10 10 ratio 1 1 1 1.25 group 1gen zone brick size 15 20 15 p0 10 0 600 p1 add 20 0 -25 p2 add 0 50 0 ... p3 add 0 0 15 p4 add 20 50 -25 p5 add 0 50 15 ... p6 add 20 0 15 p7 add 20 50 15 ... ratio 1.25 1 1 group 2gen zone brick size 20 20 15 p0 0 0 575 p1 add 30 0 0 p2 add 0 50 0 ... p3 add 0 0 25 p4 add 30 50 0 p5 add 0 50 25 ... p6 add 10 0 25 p7 add 10 50 25 ... ratio 1 1 0.8 group 3 gen zone reflect orig 0 0 0 norm 1 0 0115 p洞室模型的生成洞室模型的生成3.FLAC3D数值模拟计算实例数值模拟计算实例116 p基坑开挖过程模拟体积模量体积模量剪切模量剪切模量内摩擦角内摩擦角2.67e7Pa0.3e7Pa35ᵒ粘聚力粘聚力抗拉强度抗拉强度初始密度初始密度3e5Pa1e3Pa2700kg/m3力学参数计算模型及监测点布置3.FLAC3D数值模拟计算实例数值模拟计算实例117 p基坑开挖过程模拟各开挖步下最大主应力分布3.FLAC3D数值模拟计算实例数值模拟计算实例118 p基坑开挖过程模拟监测节点位移及总体位移响应3.FLAC3D数值模拟计算实例数值模拟计算实例119 p路基施工过程模拟几何边界计算模型土层名称 ρ/ kg·m3 c/ kPaφ(ᵒ) E/MPaν回填土 150010158.00.33粘土180020204.00.33土层物理力学参数AB3.FLAC3D数值模拟计算实例数值模拟计算实例120 p路基施工过程模拟路基回填过程中铅直向应力云图3.FLAC3D数值模拟计算实例数值模拟计算实例121 p路基施工过程模拟监测点z向位移监测点z向位移填筑结束时沉降云图填筑结束时水平位移云图3.FLAC3D数值模拟计算实例数值模拟计算实例122 p荷载引起地基土体的超孔隙水压力模拟体积模量0.5GPa剪切模量0.3Gpa粘聚力10kPa摩擦角15ᵒ流体模量2GPa力学参数几何边界3.FLAC3D数值模拟计算实例数值模拟计算实例123 p荷载引起地基土体的超孔隙水压力模拟Ø缓慢加载def ramp stress = min(1.0,float(step)/200.0)end@rampapply nstress = -40e3 hist @stress range x -.1 3.1 z 9.9 10.13.FLAC3D数值模拟计算实例数值模拟计算实例124 p荷载引起地基土体的超孔隙水压力模拟荷载作用下的孔压云图荷载作用下的瞬时沉降3.FLAC3D数值模拟计算实例数值模拟计算实例125 p边坡安全系数求解杨氏模量泊松比重度粘聚力摩擦角抗拉强度14.0MPa0.320.0kN/m310kPa30ᵒ0力学参数3.FLAC3D数值模拟计算实例数值模拟计算实例126 p边坡安全系数求解塑性区分布3.FLAC3D数值模拟计算实例数值模拟计算实例127 p边坡安全系数求解最大剪切应变率云图3.FLAC3D数值模拟计算实例数值模拟计算实例128 p隧道开挖及支护模拟FLAC3D数值计算模型3.FLAC3D数值模拟计算实例数值模拟计算实例129 p隧道开挖及支护模拟3.FLAC3D数值模拟计算实例数值模拟计算实例130 p隧道开挖及支护模拟半断面开挖完成3.FLAC3D数值模拟计算实例数值模拟计算实例131 p隧道开挖及支护模拟半断面开挖完成（初期支护）3.FLAC3D数值模拟计算实例数值模拟计算实例132 p隧道开挖及支护模拟全断面开挖完成3.FLAC3D数值模拟计算实例数值模拟计算实例133 p隧道开挖及支护模拟全断面开挖完成（二次衬砌）3.FLAC3D数值模拟计算实例数值模拟计算实例134 p隧道开挖及支护模拟洞顶铅直位移洞壁水平位移洞底铅直位移最大不平衡力3.FLAC3D数值模拟计算实例数值模拟计算实例135 p矿体开挖回填模拟矿体开挖回填模拟开挖方向框中数字表示开挖顺序计算模型概述：边界尺寸400m×400m, 模型顶部为地表，底部埋深400m. 矿体边界尺寸为50m×10m, 位于模型的正中心。
矿体分10步开挖完成，每步开挖5m，采用“隔一采一”的方式进行矿体开挖顺序如左下图所示计算模型共分为14500单元，17766节点3.FLAC3D数值模拟计算实例数值模拟计算实例136 p矿体开挖回填模拟初始状态铅直应力及位移3.FLAC3D数值模拟计算实例数值模拟计算实例137 p矿体开挖回填模拟矿体开挖回填模拟第一步开挖完成3.FLAC3D数值模拟计算实例数值模拟计算实例138 p矿体开挖回填模拟矿体开挖回填模拟第二步开挖完成3.FLAC3D数值模拟计算实例数值模拟计算实例139 p矿体开挖回填模拟矿体开挖回填模拟第三步开挖完成3.FLAC3D数值模拟计算实例数值模拟计算实例140 p矿体开挖回填模拟矿体开挖回填模拟第四步开挖完成3.FLAC3D数值模拟计算实例数值模拟计算实例141 p矿体开挖回填模拟矿体开挖回填模拟第五步开挖完成3.FLAC3D数值模拟计算实例数值模拟计算实例142 p矿体开挖回填模拟矿体开挖回填模拟第六步开挖完成3.FLAC3D数值模拟计算实例数值模拟计算实例143 p矿体开挖回填模拟矿体开挖回填模拟第七步开挖完成3.FLAC3D数值模拟计算实例数值模拟计算实例144 p矿体开挖回填模拟矿体开挖回填模拟第八步开挖完成3.FLAC3D数值模拟计算实例数值模拟计算实例145 p矿体开挖回填模拟矿体开挖回填模拟第九步开挖完成3.FLAC3D数值模拟计算实例数值模拟计算实例146 p矿体开挖回填模拟矿体开挖回填模拟第十步开挖完成3.FLAC3D数值模拟计算实例数值模拟计算实例147 p矿体开挖回填模拟矿体开挖回填模拟第十一步开挖完成3.FLAC3D数值模拟计算实例数值模拟计算实例148 p矿体开挖回填模拟矿体开挖回填模拟矿体顶板沉降位移3.FLAC3D数值模拟计算实例数值模拟计算实例149 p矿体开挖回填模拟矿体开挖回填模拟矿体底板回弹位移3.FLAC3D数值模拟计算实例数值模拟计算实例150 p矿体开挖回填模拟矿体开挖回填模拟地表沉降位移3.FLAC3D数值模拟计算实例数值模拟计算实例151 p具有自由水面的稳定流具有自由水面的稳定流几何边界示意图permeability (k) 10−10 (m/sec)/(Pa/m)porosity (n) 0.3water density (ρw) 1000 kg/m3water bulk modulus (Kf ) 1000 Pasoil dry density (ρ) 2000 kg/m3L = 9 mh1 = 6 mh2 = 1.2 m3.FLAC3D数值模拟计算实例数值模拟计算实例152 p具有自由水面的稳定流具有自由水面的稳定流3.FLAC3D数值模拟计算实例数值模拟计算实例153 p具有自由水面的稳定流具有自由水面的稳定流3.FLAC3D数值模拟计算实例数值模拟计算实例154 p具有自由水面的稳定流Ø单宽流量可根据Duouit’s formula求解Ø理论解与数值解对比Ø理论解 1.920e-6m3/sØ数值解 1.912e-6m3/s3.FLAC3D数值模拟计算实例数值模拟计算实例155 Thank You!156 。