flac3d5.0结构单元教程（课件分享）

1、作者：作者：Dr.FengDr.Feng flac3d5.0结构单元教程（课件分享） flac3d5.0结构单元教程（课件分享） StructuralElement 2 flac3d5.0结构单元教程（课件分享） FLAC3D结构单元 1.结构单元的类型 2.结构单元的建模方法 3.结构单元的参数取值 4.结构单元实例分析 5.关于link 3 flac3d5.0结构单元教程（课件分享） FLAC3D中包含六种形式的结构单元，可以分成两类：中包含六种形式的结构单元，可以分成两类： 线型结构单元线型结构单元： 梁单元梁单元(beam) 锚索单元锚索单元(cable) 桩单元桩单元(pile) 壳

2、型结构单元：壳型结构单元： 壳单元壳单元(shell) 土工格栅土工格栅(geogrid) 衬砌单元衬砌单元(liner) FLAC3D中的结构单元是岩土工 程中实际结构的一种“抽象”， 即采用简单的单元形式来模拟复 杂的结构体。 结构单元由结构节点(node)和结 构构件(SELs)构成。 结构单元中的节点(node)可以与 周围的实体网格(zone)或其它结 构节点建立连接(link)，通过连接 实现岩土体或结构与其它结构发 生相互作用。 注意：结构节点并不是简单地与 实体网格的节点(gridpoint)建立联 系，也不能建立node与gridpoint 之间的link 1、结构单元的类型

3、 4 flac3d5.0结构单元教程（课件分享） 梁单元梁单元 sel beam id 1 beg 4 0 -1 end 5 0 -2 nseg 4 sel node id=1 0 0 0 sel node id=2 2 0 0 sel node id=3 4 0 -1 sel node id=4 5 0 -2 sel beamsel cid=1 id=1 node 1 2 ; sel beamsel cid=2 id=1 node 2 3 sel beamsel cid=3 id=1 node 3 4 桩单元桩单元 sel pile id 1 beg 0 0 0 end 0 0 10 nse

4、g 4 2、结构单元的建模方法 两种建模方式各 有各的优点，第 二种方式适合建 立复杂曲线结构 单元（但是要注但是要注 意它不会自动建意它不会自动建 立立linklink！若不！若不 手动手动linklink就无任就无任 何作用何作用） 5 flac3d5.0结构单元教程（课件分享） 锚索单元锚索单元 sel cable id 1 beg 4 0 -1 end 5 0 -2 nseg 4 2、结构单元的建模方法 6 flac3d5.0结构单元教程（课件分享） 建立梁单元，并显 示单元坐标系！ 2、结构单元的建模方法线型结构单元 起始点坐标并给定分段数目的方法； 7 flac3d5.0结构单元教

5、程（课件分享） ID号相同，共用Node，ID不同，各个ID对应的结构单元有各自独立的 node。除非设置联系，否则即使节点位于同一位置也不会传递力。 结构单元的显示！GUI操作和命令操作（manual）！ 调整好显示效果后可以将显示的命令文件另存出来，以备下次使用。 （最适用于几何模型相同，参数不同的，不同工况分析的比较） 2、结构单元的建模方法线型结构单元 8 flac3d5.0结构单元教程（课件分享） 先建立节点再联接成单元的方法； 2、结构单元的建模方法线型结构单元 9 flac3d5.0结构单元教程（课件分享） 壳单元壳单元 2、结构单元的建模方法壳型结构单元 10 flac3d5.

6、0结构单元教程（课件分享） def set_vals global ptA = 25.0 * sin( 40.0*degrad ) ; global ptB = 25.0 * cos( 40.0*degrad ) end set_vals generate zone cylinder p0=( 0.0, 0.0, 0.0 ) delete all zones sel node init zpos add -25.0 2、结构单元的建模方法壳型结构单元 11 flac3d5.0结构单元教程（课件分享） 通过附着在实体网格表面来生成shell单 元。 The shells can then be

7、repositioned if ecessary by using the SEL node init command 2、结构单元的建模方法壳型结构单元 12 flac3d5.0结构单元教程（课件分享） FLAC3D是岩土工程的专业软件，因此一般很少用来做专门的是岩土工程的专业软件，因此一般很少用来做专门的 结构分析。在涉及到结构单元的问题中，往往都要考虑结构与结构分析。在涉及到结构单元的问题中，往往都要考虑结构与 周围的实体单元的相互作用。在结构单元的建模时要特别注意周围的实体单元的相互作用。在结构单元的建模时要特别注意 一个基本原则：一个一个基本原则：一个zone至多包含一个至多包含一个

8、structure node! 因此在建立线型结构单元时，要特别注意因此在建立线型结构单元时，要特别注意nseg变量的大小。变量的大小。 nseg太小则会导致计算不精确，而太大就会违反结构单元建模太小则会导致计算不精确，而太大就会违反结构单元建模 的基本原则。的基本原则。 2、结构单元的建模方法注意事项 13 flac3d5.0结构单元教程（课件分享） 梁单元梁单元 emod弹性模量，弹性模量，E nu泊松比，泊松比， xcarea横截面积，横截面积，A xciy梁结构梁结构y轴惯性矩轴惯性矩, Iy xciz梁结构梁结构z轴惯性矩，轴惯性矩， Ix xcij极惯性矩，极惯性矩，J densi

9、ty密度，密度， 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热膨胀系数热膨胀系

10、数 3、结构单元的参数取值 14 flac3d5.0结构单元教程（课件分享） 3、结构单元的参数取值 某些结构单元参数的取值要视具体情况而定，根据经验且必要 时调整参数通过试算来确定。 15 flac3d5.0结构单元教程（课件分享） 4、结构单元实例分析 4.1、简支梁（beam单元）承受两个相等集中载荷 4.2、简支梁（shell单元）承受两个相等集中载荷 16 flac3d5.0结构单元教程（课件分享） Simple Beam Two Equal Concentrated Loads 4.1、简支梁（beam单元）承受两个相等集中载荷 17 flac3d5.0结构单元教程（课件分享） A

11、 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 b

12、etween 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) as 4.1、简支梁（beam单元）承受两个相等集中载荷 18 flac3d5.0结构单元教程（课件分享） 载荷(N) 载荷距支 座的距离 (m) 铰支座之间 的距离(m) 弹性模 量(Pa) 惯性矩 (m4) 惯性矩 (m5) 最大挠度(m) PaLEIyIzmax 10000392.00E+11 2.00E-04

13、2.00E-04 0.006468750 根据理论公式计算得到： 4.1、简支梁（beam单元）承受两个相等集中载荷 19 flac3d5.0结构单元教程（课件分享） new title 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=3 sel beam id=1 begin=

14、( 3, 0, 0) end=( 6, 0, 0) nseg=4 sel beam id=1 begin=( 6, 0, 0) end=( 9, 0, 0) nseg=3 ; = ; Assign beam properties sel beam id=1 prop emod=2e11 nu=0.30 = = ; Specify model boundary conditions (including applied loads) sel node fix z xr yr ; restrict all non-beam modes sel node fix y range id=1 ; sel

15、 node fix y range id=9 ; ; rollers at beam ends sel node apply force=(0.0,-1e4,0.0) range id=2 ; apply point loads sel node apply force=(0.0,-1e4,0.0) range id=5 ; 4.1、简支梁（beam单元）承受两个相等集中载荷 20 flac3d5.0结构单元教程（课件分享） ; = = ; Setup histories for monitoring behavior. history add id=10 sel node ydisp id=

16、7 history add id=30 sel beamsel moment mz end2 cid=1 ; moment, right of SEL-1 history add id=31 sel beamsel moment mz end1 cid=2 ; moment, left of SEL-2 ; = = ; Bring the problem to equilibrium solve ratio=1e-7 save equal-concent-loads ; = = ; Print out beam responses. list sel beam force list sel b

17、eam moment list sel node disp range id=7 return 4.1、简支梁（beam单元）承受两个相等集中载荷 如何设置结构单元的跟踪 变量！ 21 flac3d5.0结构单元教程（课件分享） Beam_concent_loads_Examp le1.3 4.1、简支梁（beam单元）承受两个相等集中载荷 22 flac3d5.0结构单元教程（课件分享） 挠度计算 4.1、简支梁（beam单元）承受两个相等集中载荷 23 flac3d5.0结构单元教程（课件分享） 剪力、弯矩计算 24 flac3d5.0结构单元教程（课件分享） 剪 力、 弯 矩 计 算 这

18、是节点力！ 25 flac3d5.0结构单元教程（课件分享） 26 flac3d5.0结构单元教程（课件分享） 梁单元局部坐标系：x轴从节点1到节点2，y轴在横截面中 4.1、简支梁（beam单元）承受两个相等集中载荷 27 flac3d5.0结构单元教程（课件分享） 弯矩矢量的指向，右手法则！弯矩矢量的指向，右手法则！ 4.1、简支梁（beam单元）承受两个相等集中载荷 28 flac3d5.0结构单元教程（课件分享） List sel beam nodal forces: components are displayed in terms of the beam local coordin

19、ate systems. These are the forces exerted by the nodes on the beamSEL. 小结：梁单元的常用命令 4.1、简支梁（beam单元）承受两个相等集中载荷 history sel beamsel cid怎么找？（坐标or鼠标 information？） 29 flac3d5.0结构单元教程（课件分享） sel Sel node 命令 针对所有的结构单元 30 flac3d5.0结构单元教程（课件分享） Sel node fix keyword . . . 31 flac3d5.0结构单元教程（课件分享） new title Simp

20、le Beam (modeled using shellSELs) gen zone brick size 12,3,1 shell 宽度为单位1（z方向上） sel shell id=1 crossdiag elemtype=dkt range y -0.1 0.1 sel shell id=1 prop iso=(2e11, 0.0) thick=0.133887 delete zone sel node fix x y xr yr range x=(-0.1, 0.1) ; support at left end hinge 铰支座 sel node fix y xr yr range

21、x=( 8.9, 9.1) ; support at right end roller 辊轴支 座 sel node fix z xr yr ; restrict non-beam deformation modes sel node apply force=(0,-1667,0) range union id=71 id=12 ; out nds, left sel node apply force=(0,-3333,0) range union id=46 id=13 ; in nds, left sel node apply force=(0,-1667,0) range union i

22、d=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 flac3d5.0结构单元教程（课件分享） history add id=1 unbal history add id=10 sel node ydisp id=19 ; displ at center ; moment, left third history add id=20 sel recover sres

23、Mx surfX 1,0,0 cid=59 ; shear, left third history add id=30 sel recover sres Qx surfX 1,0,0 cid=59 solve ratio=1e-7 list sel node disp range id=19 save shell0 return 4.2、简支梁（shell单元）承受两个相等集中载荷 33 flac3d5.0结构单元教程（课件分享） 4.2、简支梁（shell单元）承受两个相等集中载荷 34 flac3d5.0结构单元教程（课件分享） sel shell id=1 elemtype=cst ra

24、nge y -0.1 0.1 ;no-crossdiag sel shell id=1 crossdiag elemtype=dkt range y -0.1 0.1 Crossdiag vs no- crossdiag 4.2、简支梁（shell单元）承受两个相等集中载荷 35 flac3d5.0结构单元教程（课件分享） 局部坐标系！看弯矩到底 应该看哪一个？ 36 flac3d5.0结构单元教程（课件分享） 4.2、简支梁（shell单元）承受两个相等集中载荷 37 flac3d5.0结构单元教程（课件分享） 4.2、简支梁（shell单元）承受两个相等集中载荷 38 flac3d5.0结

25、构单元教程（课件分享） 39 flac3d5.0结构单元教程（课件分享） 4.2、简支梁（shell单元）承受两个相等集中载荷 40 flac3d5.0结构单元教程（课件分享） surfX Xx Xy Xz The 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, xyz, has the follo

26、wing 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

27、 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 a

28、lso 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 automati

29、cally 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 flac3d5.0结构单元教程（课件分享） 5、关于link 两种： Node-Zone Node-Node 42 fla

30、c3d5.0结构单元教程（课件分享） 创建一个新的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确定

31、。如果不能确 定源和目标对象，就会报错，且该命令不会对模型产 生任何作用。新link的attachment条件设置为6个自 由度均为“rigid”。 Side1，side2关键字对于确定嵌入式liner的哪个面 上产生link。 sel 可选参数与必选参数！可选参数与必选参数！ 43 flac3d5.0结构单元教程（课件分享） 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 6 sel cable id=1 prop

32、emod 2e10 ytension 310e3 xcarea 0.0004906 group clay (old command） 65 flac3d5.0结构单元教程（课件分享） 1、手册中给出的接触面建立方法So-called “移来移去法” interface 1 face range cylinder end1 (0,0,0) end2 (0,0,-5.1) radius .31 写简写的时候要注意新老版本的区别 gen zone cylinder p0 (0,0,6.1) p1 (.3,0,6.1) p2 (0,0,6) p3 (0,.3,6.1) ini z add -6.0 r

33、ange group pile(old command） 66 flac3d5.0结构单元教程（课件分享） 1、手册中给出的接触面建立方法So-called “移来移去法” 67 flac3d5.0结构单元教程（课件分享） 68 flac3d5.0结构单元教程（课件分享） 69 flac3d5.0结构单元教程（课件分享） 70 flac3d5.0结构单元教程（课件分享） “移来移去”法 ; Create Base gen zone brick size 3 3 3 Create Top - 1 unit high for initial spacing gen zone brick size

34、3 3 3 ; Create interface elements on the top surface of the base interface 1 face range plane norm (-1,0,1) origin (1.5,1.5,3) dist 0.1 ; Lower top to complete geometry ini zpos add -1.0 range group Top 71 flac3d5.0结构单元教程（课件分享） 2、接触面建立方法So-called “倒来倒去法” 我们最终的目的就是在中心小块体与外围网格之间建立接触面。 1. 分开建立网格 2. 建立i

35、nner网格及其表面的Interface 3. 导入外围mesh 4. 赋予材料属性，测试接触面是否发生了作用。 72 flac3d5.0结构单元教程（课件分享） 小练习： 1. 三种建立接触面的方法计算结果是否相同（只要接触面有响应，肯定是 相同的！） 2. 如果将接触面建立在外部网格的内表面，然后移入小块体，结果是否相 同呢？ 3. 不加接触面跟加了接触面，模型的响应（位移、应力）有何区别？Nrange Example-7.1 union nrange 73 flac3d5.0结构单元教程（课件分享） 3、切割模型的方法实际上是分离连续网格（原来网格连续，通过 共用节点（GridPoint

36、）传递力，分离后通过接触面来传递。 gen zone brick size 3 3 3 group zone inner range x 1 2 y 1 2 z 2 3 group zone out range group inner not generate separate face group aa range group inner group out ;very different from old version interface 1 wrap first group inner ;second group out ;interface 1 permeability on int

37、erface 1 maxedge 0.5 74 flac3d5.0结构单元教程（课件分享） 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.

38、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

39、. 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 comm

40、and 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个类组 为整个模型，边

41、界面为整个第1个类组的表面。 75 flac3d5.0结构单元教程（课件分享） Generate separate 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 belongi

42、ng 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, then range group Fred group Ge

43、orge ;生成的应该是一个共享面 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 flac3

44、d5.0结构单元教程（课件分享） 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 gridpoint

45、s 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.

46、FLAC3D determines which face will be new and which one will be old by using the origin keyword. originx y z specifies 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 t

47、he 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 flac3d5.0结构单元教程（课件分享） 检验是否分离开了网格并形成了接触面？ 78 flac3d5.0结构单元教程（课件分享） 新版本中的Facegroup法 ; select as regional joining Base

48、 and Top group face Interface internal range group Base group Top ; Separate Interface faces, calling the newwly created faces NewFaces gen separate face origin (0,0,0) group NewFaces range group Interface ; ; Create interface elements on the top surface of the base interface 1 face range group Inte

49、rface 79 flac3d5.0结构单元教程（课件分享） 80 flac3d5.0结构单元教程（课件分享） 实例BinFlowSlip ; Assign group names to faces that will become interfaces group face Int1 internal range group Material group Bin . plane orig (0,0,0) normal (1,-1,0) above group face Int2 internal range group Material group Bin . plane orig (0,0

50、,0) normal (1,-1,0) below ; Separate those faces, giving newly created surface faces a new group name gen separate face origin (0,0,0) group NewFaces range group Int1 or Int2 ; Created interfaces on Int1 and Int2 faces interface 1 face range group Int1 interface 2 face range group Int2 ; Subdivide i

51、nterface elements a little for better contact detection int 1 maxedge 0.55 int 2 maxedge 0.55 81 flac3d5.0结构单元教程（课件分享） 1. FLAC3D基本操作方法 82 flac3d5.0结构单元教程（课件分享） 1. FLAC3D基本操作方法 p软件界面 FLAC3D5.00 软件界 面 83 flac3d5.0结构单元教程（课件分享） 1. FLAC3D基本操作方法 p软件界面 FLAC3D5.00 软件界 面 84 flac3d5.0结构单元教程（课件分享） 1. FLAC3D基本操

52、作方法 p命令流操作方法 new gen zone brick size 6 6 6 在利用FLAC3D软件进行数值模拟时，主要是通过命令流来实现 的。命令流文件一般以txt或dat格式存储，并在命令输入窗口、菜单 栏或快捷图标通过call命令进行调用。 FLAC3D命令流文件需要 遵循一定的格式和语法要 求，在满足这些要求的前 提下，命令流文件的编写 又是十分自由灵活的，可 根据用户个人的需求自由 编写。 85 flac3d5.0结构单元教程（课件分享） 1. FLAC3D基本操作方法 p数值计算一般流程 u建立模型 通过外部导入或在FLAC3D 中直接建模的方式建立计算模 型。 u材料参数

53、 材料本构 材料力学参数 u边界条件及初始条件 速度边界 应力边界 水头边界 求解计算求解计算 边界 条件 及初 始条 件 建立 模型 材料 参数 86 flac3d5.0结构单元教程（课件分享） 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

54、0 0 6 (ratio 1 1 1 ) (dim 2 2 2) (fill) 87 flac3d5.0结构单元教程（课件分享） 1. FLAC3D基本操作方法 p建立模型 FLAC3D建模的基本思路为“堆积木”，即首先建立各种形状的网格单元，最后将建立的网格单元组合在一起，生成可用于数值计算的整体模型。 FLAC3D内置13种不同形状的网格，包括brick(砖形), cshell(圆柱状壳网格), cylinder(圆柱状网格), cylint(圆柱状交叉网格), dbrick(退化砖形网格), pyramid(锥形网格), radbrick(砖形辐射网格), radcylinder(圆柱状

55、辐射网格), radtunnel(平行六边形状辐射网格), retrahedron(辐射网格), tunint(砖形交叉网格), uwedge(均匀 楔形网格), wedge(楔形网格). 88 flac3d5.0结构单元教程（课件分享） 1. FLAC3D基本操作方法 p网格形状 89 flac3d5.0结构单元教程（课件分享） 1. FLAC3D基本操作方法基本操作方法 p brick gen zone brick size 6 8 8 plot zone 90 flac3d5.0结构单元教程（课件分享） 1. FLAC3D基本操作方法基本操作方法 p dbrick gen zone db

56、rick 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 plot zone 91 flac3d5.0结构单元教程（课件分享） 1. FLAC3D基本操作方法基本操作方法 p wedge gen zone wedge size 6 6 8 p0 0 0 0 p1 6 0 0 p2 0 6 0 p3 0 0 6 plot zone 92 flac3d5.0结构单元教程（课件分享） 1. FLAC3D基本操作方法基本操作方法 p uwedge gen zone uwedge size 6 6 6 p

57、0 0 0 0 p1 6 0 0 p2 0 6 0 p3 0 0 6 plot zone 93 flac3d5.0结构单元教程（课件分享） 1. FLAC3D基本操作方法基本操作方法 p pyramid gen 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 0 plot zone 94 flac3d5.0结构单元教程（课件分享） 1. FLAC3D基本操作方法基本操作方法 p tetrahedron gen zone tetrahedron size 6 6 6 p0 0 0 0 p1 6 0 0 p2 0

58、 6 0 p3 0 0 6 plot zone 95 flac3d5.0结构单元教程（课件分享） 1. FLAC3D基本操作方法基本操作方法 p cylinder gen zone cyl size 6 4 10 p0 0 0 0 p1 6 0 0 p2 0 6 0 p3 0 0 6 plot zone 96 flac3d5.0结构单元教程（课件分享） 1. FLAC3D基本操作方法基本操作方法 p radbrick gen 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

59、1.2 (fill) plot zone 97 flac3d5.0结构单元教程（课件分享） 1. FLAC3D基本操作方法基本操作方法 p radcylinder gen 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 zone 98 flac3d5.0结构单元教程（课件分享） 1. FLAC3D基本操作方法基本操作方法 p通过对称生成网格 gen zone radbrick . p0 (0,0,0) p1 (10,0,0) p2

60、 (0,10,0) p3 (0,0,10) . size 3,5,5,7 . ratio 1,1,1,1.5 . dim 1 4 2 fill gen zone reflect dip 0 dd 90 gen zone reflect dip 90 dd 90 plot zone 99 flac3d5.0结构单元教程（课件分享） 1. FLAC3D基本操作方法基本操作方法 p生成渐变网格 gen zone brick size 10 10 10 p0 0 0 0 定义函数名 ;函数语句 end ;函数结束的标志 pFish函数举例 def func b=10 a=b+10 end func p