flac3d5.0結(jié)構(gòu)單元教程【精制材料】

1、FLAC3D 5.0培訓(xùn)日程安排,FLAC3D V5.0界面操作 FLAC3D基本操作方法vs應(yīng)用流程； FLAC3D內(nèi)置Fish語(yǔ)言的應(yīng)用； FLAC3D結(jié)構(gòu)單元vs接觸單元； FLAC3D滲流模塊 其他,1,實(shí)操應(yīng)用,StructuralElement,2,實(shí)操應(yīng)用,FLAC3D結(jié)構(gòu)單元,結(jié)構(gòu)單元的類(lèi)型 結(jié)構(gòu)單元的建模方法 結(jié)構(gòu)單元的參數(shù)取值 結(jié)構(gòu)單元實(shí)例分析 關(guān)于link,3,實(shí)操應(yīng)用,FLAC3D中包含六種形式的結(jié)構(gòu)單元，可以分成兩類(lèi)： 線型結(jié)構(gòu)單元： 梁?jiǎn)卧?beam) 錨索單元(cable) 樁單元(pile) 殼型結(jié)構(gòu)單元： 殼單元(shell) 土工格柵(geogrid) 襯砌

2、單元(liner),FLAC3D中的結(jié)構(gòu)單元是巖土工程中實(shí)際結(jié)構(gòu)的一種“抽象”，即采用簡(jiǎn)單的單元形式來(lái)模擬復(fù)雜的結(jié)構(gòu)體。 結(jié)構(gòu)單元由結(jié)構(gòu)節(jié)點(diǎn)(node)和結(jié)構(gòu)構(gòu)件(SELs)構(gòu)成。 結(jié)構(gòu)單元中的節(jié)點(diǎn)(node)可以與周?chē)膶?shí)體網(wǎng)格(zone)或其它結(jié)構(gòu)節(jié)點(diǎn)建立連接(link)，通過(guò)連接實(shí)現(xiàn)巖土體或結(jié)構(gòu)與其它結(jié)構(gòu)發(fā)生相互作用。 注意：結(jié)構(gòu)節(jié)點(diǎn)并不是簡(jiǎn)單地與實(shí)體網(wǎng)格的節(jié)點(diǎn)(gridpoint)建立聯(lián)系，也不能建立node與gridpoint之間的link,1、結(jié)構(gòu)單元的類(lèi)型,4,實(shí)操應(yīng)用,梁?jiǎn)卧?sel beam id 1 beg 4 0 -1 end 5 0 -2 nseg 4 sel node

3、 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 nseg 4,2、結(jié)構(gòu)單元的建模方法,兩種建模方式各有各的優(yōu)點(diǎn)，第二種方式適合建立復(fù)雜曲線結(jié)構(gòu)單元（但是要注意它不會(huì)自動(dòng)建立link！若不手動(dòng)link就無(wú)任何作用）,5,實(shí)操應(yīng)用

4、,錨索單元 sel cable id 1 beg 4 0 -1 end 5 0 -2 nseg 4,2、結(jié)構(gòu)單元的建模方法,6,實(shí)操應(yīng)用,建立梁?jiǎn)卧?，并顯示單元坐標(biāo)系！,2、結(jié)構(gòu)單元的建模方法線型結(jié)構(gòu)單元,起始點(diǎn)坐標(biāo)并給定分段數(shù)目的方法；,7,實(shí)操應(yīng)用,ID號(hào)相同，共用Node，ID不同，各個(gè)ID對(duì)應(yīng)的結(jié)構(gòu)單元有各自獨(dú)立的node。除非設(shè)置聯(lián)系，否則即使節(jié)點(diǎn)位于同一位置也不會(huì)傳遞力。,結(jié)構(gòu)單元的顯示！GUI操作和命令操作（manual）！ 調(diào)整好顯示效果后可以將顯示的命令文件另存出來(lái)，以備下次使用。（最適用于幾何模型相同，參數(shù)不同的，不同工況分析的比較）,2、結(jié)構(gòu)單元的建模方法線型結(jié)構(gòu)單元,8

5、,實(shí)操應(yīng)用,先建立節(jié)點(diǎn)再聯(lián)接成單元的方法；,2、結(jié)構(gòu)單元的建模方法線型結(jié)構(gòu)單元,9,實(shí)操應(yīng)用,殼單元,2、結(jié)構(gòu)單元的建模方法殼型結(jié)構(gòu)單元,10,實(shí)操應(yīng)用,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、結(jié)構(gòu)單元的建模方法殼型結(jié)構(gòu)單元,11,實(shí)操應(yīng)用,通

6、過(guò)附著在實(shí)體網(wǎng)格表面來(lái)生成shell單元。,The shells can then be repositioned if ecessary by using the SEL node init command,2、結(jié)構(gòu)單元的建模方法殼型結(jié)構(gòu)單元,12,實(shí)操應(yīng)用,FLAC3D是巖土工程的專(zhuān)業(yè)軟件，因此一般很少用來(lái)做專(zhuān)門(mén)的結(jié)構(gòu)分析。在涉及到結(jié)構(gòu)單元的問(wèn)題中，往往都要考慮結(jié)構(gòu)與周?chē)膶?shí)體單元的相互作用。在結(jié)構(gòu)單元的建模時(shí)要特別注意一個(gè)基本原則：一個(gè)zone至多包含一個(gè)structure node! 因此在建立線型結(jié)構(gòu)單元時(shí)，要特別注意nseg變量的大小。nseg太小則會(huì)導(dǎo)致計(jì)算不精確，而太大就會(huì)違反

7、結(jié)構(gòu)單元建模的基本原則。,2、結(jié)構(gòu)單元的建模方法注意事項(xiàng),13,實(shí)操應(yīng)用,梁?jiǎn)卧?emod彈性模量，E nu泊松比， xcarea橫截面積，A xciy梁結(jié)構(gòu)y軸慣性矩, Iy xciz梁結(jié)構(gòu)z軸慣性矩，Ix xcij極慣性矩，J density密度， pmoment塑性矩，Mp thexp熱膨脹系數(shù)，t ydirection矢量Y,錨索單元 emod彈性模量， E xcarea橫截面積，A gr_coh單位長(zhǎng)度上水泥漿粘結(jié)力cg gr_fric水泥漿的摩擦角g gr_k單位長(zhǎng)度上水泥漿剛度kg gr_per水泥漿外圈周長(zhǎng)Pg slide大變形滑動(dòng)標(biāo)志 slide_tol大變形滑動(dòng)容差 yco

8、mp抗壓強(qiáng)度(力) density密度 thexp熱膨脹系數(shù),3、結(jié)構(gòu)單元的參數(shù)取值,14,實(shí)操應(yīng)用,3、結(jié)構(gòu)單元的參數(shù)取值,某些結(jié)構(gòu)單元參數(shù)的取值要視具體情況而定，根據(jù)經(jīng)驗(yàn)且必要時(shí)調(diào)整參數(shù)通過(guò)試算來(lái)確定。,15,實(shí)操應(yīng)用,4、結(jié)構(gòu)單元實(shí)例分析,4.1、簡(jiǎn)支梁（beam單元）承受兩個(gè)相等集中載荷 4.2、簡(jiǎn)支梁（shell單元）承受兩個(gè)相等集中載荷,16,實(shí)操應(yīng)用,Simple Beam Two Equal Concentrated Loads,4.1、簡(jiǎn)支梁（beam單元）承受兩個(gè)相等集中載荷,17,實(shí)操應(yīng)用,A simply supported beam is loaded by two

9、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.

10、 The maximum deflection of the beam,max, occurs at the center and is given by AISC (1980, p. 2-116) as,4.1、簡(jiǎn)支梁（beam單元）承受兩個(gè)相等集中載荷,18,實(shí)操應(yīng)用,根據(jù)理論公式計(jì)算得到：,4.1、簡(jiǎn)支梁（beam單元）承受兩個(gè)相等集中載荷,19,實(shí)操應(yīng)用,new title Simple Beam - Two Equal Concentrated Loads Symmetrically Placed ; = ; Create the grid, insure that nodes wi

11、ll exist at third points. sel beam id=1 begin=( 0, 0, 0) end=( 3, 0, 0) nseg=3 sel beam id=1 begin=( 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 ,4.1、簡(jiǎn)支梁（beam單元）承受兩個(gè)相等集中載荷,20,實(shí)操應(yīng)用,; = ; Setup hi

12、stories for monitoring behavior. history add id=10 sel node ydisp id=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-loa

13、ds ; = ; Print out beam responses. list sel beam force list sel beam moment list sel node disp range id=7 return,4.1、簡(jiǎn)支梁（beam單元）承受兩個(gè)相等集中載荷,如何設(shè)置結(jié)構(gòu)單元的跟蹤變量！,21,實(shí)操應(yīng)用,Beam_concent_loads_Example1.3,4.1、簡(jiǎn)支梁（beam單元）承受兩個(gè)相等集中載荷,22,實(shí)操應(yīng)用,撓度計(jì)算,4.1、簡(jiǎn)支梁（beam單元）承受兩個(gè)相等集中載荷,23,實(shí)操應(yīng)用,剪力、彎矩計(jì)算,24,實(shí)操應(yīng)用,剪力、彎矩計(jì)算,這是節(jié)點(diǎn)力！,25,實(shí)

14、操應(yīng)用,26,實(shí)操應(yīng)用,梁?jiǎn)卧植孔鴺?biāo)系：x軸從節(jié)點(diǎn)1到節(jié)點(diǎn)2，y軸在橫截面中,4.1、簡(jiǎn)支梁（beam單元）承受兩個(gè)相等集中載荷,27,實(shí)操應(yīng)用,彎矩矢量的指向，右手法則！,4.1、簡(jiǎn)支梁（beam單元）承受兩個(gè)相等集中載荷,28,實(shí)操應(yīng)用,List sel beam nodal forces: components are displayed in terms of the beam local coordinate systems. These are the forces exerted by the nodes on the beamSEL.,小結(jié)：梁?jiǎn)卧某Ｓ妹?4.1、簡(jiǎn)支梁（

15、beam單元）承受兩個(gè)相等集中載荷,history sel beamsel,cid怎么找？（坐標(biāo)or鼠標(biāo)information？）,29,實(shí)操應(yīng)用,sel,Sel node 命令 針對(duì)所有的結(jié)構(gòu)單元,30,實(shí)操應(yīng)用,Sel node fix keyword . . .,31,實(shí)操應(yīng)用,new title Simple Beam (modeled using shellSELs) gen zone brick size 12,3,1 in nds, rt.,將均布載荷轉(zhuǎn)換為等效節(jié)點(diǎn)力,4.2、簡(jiǎn)支梁（shell單元）承受兩個(gè)相等集中載荷,32,實(shí)操應(yīng)用,history add id=1 unba

16、l history add id=10 sel node ydisp id=19 ; displ at center ; moment, left third history add id=20 sel recover sres 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、簡(jiǎn)支梁（shell單元）承受兩

17、個(gè)相等集中載荷,33,實(shí)操應(yīng)用,4.2、簡(jiǎn)支梁（shell單元）承受兩個(gè)相等集中載荷,34,實(shí)操應(yīng)用,sel shell id=1 elemtype=cst range 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、簡(jiǎn)支梁（shell單元）承受兩個(gè)相等集中載荷,35,實(shí)操應(yīng)用,局部坐標(biāo)系！看彎矩到底應(yīng)該看哪一個(gè)？,36,實(shí)操應(yīng)用,4.2、簡(jiǎn)支梁（shell單元）承受兩個(gè)相等集中載荷,37,實(shí)操應(yīng)用,4.2、簡(jiǎn)支梁（shel

18、l單元）承受兩個(gè)相等集中載荷,38,實(shí)操應(yīng)用,39,實(shí)操應(yīng)用,4.2、簡(jiǎn)支梁（shell單元）承受兩個(gè)相等集中載荷,40,實(shí)操應(yīng)用,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 following properties: (1) z

19、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 pr

20、ocessing 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 indiv

21、idual 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 u

22、nder 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,實(shí)操應(yīng)用,5、關(guān)于link,兩種： Node-Zone Node-Node,42,實(shí)操應(yīng)用,創(chuàng)建一個(gè)新的link，link的源節(jié)點(diǎn)為sid，而聯(lián)接目標(biāo)為node或zo

23、ne。Id為新link的ID號(hào)。Sid是已經(jīng)存在的節(jié)點(diǎn)（作為源節(jié)點(diǎn)）的ID號(hào)，可選關(guān)鍵字target用于確定目標(biāo)對(duì)象（node或zone）。默認(rèn)的目標(biāo)對(duì)象為zone。 對(duì)于zone目標(biāo)對(duì)象，如果tid沒(méi)有定義，將會(huì)使用與源節(jié)點(diǎn)距離在delta范圍內(nèi)的非空zone；否則，如果tid定義了，如果該tid所指示的單元為非空zone，且該zone的邊界距離在delta之內(nèi)，就會(huì)建立link。 對(duì)于node目標(biāo)對(duì)象，tid就必須定義了，且兩個(gè)節(jié)點(diǎn)必須彼此很靠近。由delta確定。如果不能確定源和目標(biāo)對(duì)象，就會(huì)報(bào)錯(cuò)，且該命令不會(huì)對(duì)模型產(chǎn)生任何作用。新link的attachment條件設(shè)置為6個(gè)自由度均為“

24、rigid”。 Side1，side2關(guān)鍵字對(duì)于確定嵌入式liner的哪個(gè)面上產(chǎn)生link。,sel,可選參數(shù)與必選參數(shù)！,43,實(shí)操應(yīng)用,5、關(guān)于link,以預(yù)應(yīng)力錨桿的托盤(pán)模擬為例 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 emod 2e10 ytension 310e3 xcarea 0.0004906 Create Base gen zone brick size 3 3 3 Lower top to co

25、mplete geometry ini zpos add -1.0 range group Top,71,實(shí)操應(yīng)用,2、接觸面建立方法So-called “倒來(lái)倒去法”,我們最終的目的就是在中心小塊體與外圍網(wǎng)格之間建立接觸面。 分開(kāi)建立網(wǎng)格 建立inner網(wǎng)格及其表面的Interface 導(dǎo)入外圍mesh 賦予材料屬性，測(cè)試接觸面是否發(fā)生了作用。,72,實(shí)操應(yīng)用,小練習(xí)： 三種建立接觸面的方法計(jì)算結(jié)果是否相同（只要接觸面有響應(yīng)，肯定是相同的?。?如果將接觸面建立在外部網(wǎng)格的內(nèi)表面，然后移入小塊體，結(jié)果是否相同呢？ 不加接觸面跟加了接觸面，模型的響應(yīng)（位移、應(yīng)力）有何區(qū)別？,Nrange Exa

26、mple-7.1 union nrange,73,實(shí)操應(yīng)用,3、切割模型的方法實(shí)際上是分離連續(xù)網(wǎng)格（原來(lái)網(wǎng)格連續(xù)，通過(guò)共用節(jié)點(diǎn)（GridPoint）傳遞力，分離后通過(guò)接觸面來(lái)傳遞。,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

27、wrap first group inner ;second group out ;interface 1 permeability on interface 1 maxedge 0.5,74,實(shí)操應(yīng)用,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 descri

28、ptor 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, indicat

29、ing 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 GE

30、NERATE 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 ra

31、nge x 50.0 75.0,接觸面建立在第1個(gè)類(lèi)組上面，若不指定第2個(gè)類(lèi)組，就默認(rèn)第2個(gè)類(lèi)組為整個(gè)模型，邊界面為整個(gè)第1個(gè)類(lèi)組的表面。,75,實(shí)操應(yīng)用,Generate Separate Zones 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 copie

32、s 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 na

33、mes assigned to zones, then range group Fred group George ;生成的應(yīng)該是一個(gè)共享面 will select faces that are connected to both Fred and George. Also note that faces can be selected by group directly.（此外，face也可以通過(guò)group直接選中） The following keywords can be used to affect the behavior 產(chǎn)生由range定義的內(nèi)部 face。面上的節(jié)點(diǎn)復(fù)制，并創(chuàng)建

34、新的表面。,76,實(shí)操應(yīng)用,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 w

35、ill 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. FLA

36、C3D 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 the

37、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,實(shí)操應(yīng)用,檢驗(yàn)是否分離開(kāi)了網(wǎng)格并形成了接觸面？,78,實(shí)操應(yīng)用,新版本中的Facegroup法,; select as regional joining Base and Top group face Interface internal

38、 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 Interface,79,實(shí)操應(yīng)用,80,實(shí)操應(yīng)用,實(shí)例BinFlowSlip,;

39、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,0) normal (1,-1,0) below ; Separate those faces, giving newly created s

40、urface 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 interface elements a little for better contact detection int 1 maxedge 0.

41、55 int 2 maxedge 0.55,81,實(shí)操應(yīng)用,1. FLAC3D基本操作方法,82,實(shí)操應(yīng)用,1. FLAC3D基本操作方法,軟件界面,FLAC3D5.00 軟件界面,83,實(shí)操應(yīng)用,1. FLAC3D基本操作方法,軟件界面,FLAC3D5.00 軟件界面,84,實(shí)操應(yīng)用,1. FLAC3D基本操作方法,命令流操作方法,new gen zone brick size 6 6 6,在利用FLAC3D軟件進(jìn)行數(shù)值模擬時(shí)，主要是通過(guò)命令流來(lái)實(shí)現(xiàn)的。命令流文件一般以txt或dat格式存儲(chǔ)，并在命令輸入窗口、菜單欄或快捷圖標(biāo)通過(guò)call命令進(jìn)行調(diào)用。,FLAC3D命令流文件需要遵循一定的格

42、式和語(yǔ)法要求，在滿足這些要求的前提下，命令流文件的編寫(xiě)又是十分自由靈活的，可根據(jù)用戶(hù)個(gè)人的需求自由編寫(xiě)。,85,實(shí)操應(yīng)用,1. FLAC3D基本操作方法,數(shù)值計(jì)算一般流程,建立模型 通過(guò)外部導(dǎo)入或在FLAC3D中直接建模的方式建立計(jì)算模型。 材料參數(shù) 材料本構(gòu) 材料力學(xué)參數(shù) 邊界條件及初始條件 速度邊界 應(yīng)力邊界 水頭邊界,86,實(shí)操應(yīng)用,1. FLAC3D基本操作方法,建立模型 在FLAC3D中，模型的建立通過(guò)關(guān)鍵詞generate來(lái)實(shí)現(xiàn)，其基本格式為： gen keywords1 keywords2 keywords3 For example: gen zone brick size 6

43、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,實(shí)操應(yīng)用,1. FLAC3D基本操作方法,建立模型 FLAC3D建模的基本思路為“堆積木”，即首先建立各種形狀的網(wǎng)格單元，最后將建立的網(wǎng)格單元組合在一起，生成可用于數(shù)值計(jì)算的整體模型。 FLAC3D內(nèi)置13種不同形狀的網(wǎng)格，包括brick(磚形), cshell(圓柱狀殼網(wǎng)格), cylinder(圓柱狀網(wǎng)格), cylint(圓柱狀交叉網(wǎng)格), dbrick(退化磚形網(wǎng)

44、格), pyramid(錐形網(wǎng)格), radbrick(磚形輻射網(wǎng)格), radcylinder(圓柱狀輻射網(wǎng)格), radtunnel(平行六邊形狀輻射網(wǎng)格), retrahedron(輻射網(wǎng)格), tunint(磚形交叉網(wǎng)格), uwedge(均勻楔形網(wǎng)格), wedge(楔形網(wǎng)格).,88,實(shí)操應(yīng)用,1. FLAC3D基本操作方法,網(wǎng)格形狀,89,實(shí)操應(yīng)用,1. FLAC3D基本操作方法,brick,gen zone brick size 6 8 8 plot zone,90,實(shí)操應(yīng)用,1. FLAC3D基本操作方法,dbrick,gen zone dbrick size 6 6 6 p

45、0 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,實(shí)操應(yīng)用,1. FLAC3D基本操作方法,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,實(shí)操應(yīng)用,1. FLAC3D基本操作方法,uwedge,gen zone uwedge size 6 6 6 p0 0 0 0 p1 6 0 0 p2 0 6 0 p3 0 0 6 plot zone,93,實(shí)操應(yīng)用,1. FLAC3D基本操作方

46、法,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,實(shí)操應(yīng)用,1. FLAC3D基本操作方法,tetrahedron,gen zone tetrahedron size 6 6 6 p0 0 0 0 p1 6 0 0 p2 0 6 0 p3 0 0 6 plot zone,95,實(shí)操應(yīng)用,1. FLAC3D基本操作方法,cylinder,gen zone cyl size 6 4 10 p0 0 0 0 p1 6 0 0 p2 0 6 0 p3 0 0 6

47、 plot zone,96,實(shí)操應(yīng)用,1. FLAC3D基本操作方法,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 1.2 (fill) plot zone,97,實(shí)操應(yīng)用,1. FLAC3D基本操作方法,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

48、 zone,98,實(shí)操應(yīng)用,1. FLAC3D基本操作方法,通過(guò)對(duì)稱(chēng)生成網(wǎng)格,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 fill gen zone reflect dip 0 dd 90 gen zone reflect dip 90 dd 90 plot zone,99,實(shí)操應(yīng)用,1. FLAC3D基本操作方法,生成漸變網(wǎng)格,gen zone brick size 10 10 10 p0 0 0 0 定義函數(shù)名

49、;函數(shù)語(yǔ)句 end ;函數(shù)結(jié)束的標(biāo)志 Fish函數(shù)舉例 def func b=10 a=b+10 end func print a,107,實(shí)操應(yīng)用,2. Fish 函數(shù),Fish 函數(shù)與變量 對(duì)于Fish函數(shù)和變量需要注意一下幾點(diǎn)： 函數(shù)和變量的賦值遵循數(shù)據(jù)類(lèi)型的法則； 變量和函數(shù)名的命名規(guī)則不能以數(shù)字開(kāi)頭，不能含有中文，不能含有下列字符： . , * / + - = # ( ) ; “ ” 變量和函數(shù)名不能與FLAC3D, FISH保留字相沖突； 對(duì)變量賦值時(shí)，不能將當(dāng)前函數(shù)的函數(shù)名放在“=”右邊，這樣會(huì)形成遞歸調(diào)用； 變量和函數(shù)名是全局的； 對(duì)FISH函數(shù)和變量的引用和通過(guò)如下命令來(lái)實(shí)現(xiàn)

50、： PRINT用于查看函數(shù)和變量的數(shù)值 HISTORY可對(duì)函數(shù)和變量的數(shù)值進(jìn)行記錄 SET可對(duì)變量進(jìn)行賦值,108,實(shí)操應(yīng)用,2. Fish 函數(shù),Fish 函數(shù)數(shù)據(jù)類(lèi)型 整型：-2147483648-2147483648 浮點(diǎn)型：10-30010300 字符型：以()為分界符，常用于保存時(shí)文件名的定義 指針型：表示單元和節(jié)點(diǎn)的存儲(chǔ)地址,109,實(shí)操應(yīng)用,2. Fish 函數(shù),Fish 函數(shù)語(yǔ)句,選擇語(yǔ)句 CASEOF 表達(dá)式 默認(rèn)語(yǔ)句 CASE n1 表達(dá)式值為n1時(shí)的語(yǔ)句 CASE n2 表達(dá)式值為n2時(shí)的語(yǔ)句 ENDCASE,條件語(yǔ)句 IF 條件表達(dá)式 ENDIF 注：IF語(yǔ)句內(nèi)部可繼續(xù)

51、嵌套條件語(yǔ)句,循環(huán)語(yǔ)句 LOOP var (exp1, exp2) ENDLOOP 或 LOOP WHILE 條件表達(dá)式 ENDLOOP,命令語(yǔ)句 COMMAND ENDCOMMAND,110,實(shí)操應(yīng)用,2. Fish 函數(shù),單元遍歷和節(jié)點(diǎn)遍歷 基本變量 zone_head和gp_head:分別表示單元和節(jié)點(diǎn)的頭指針 z_next()和gp_next:分別表示下一個(gè)單元（節(jié)點(diǎn)）,單元遍歷程序框架： p_z=zone_head loop while p_z # null ;語(yǔ)句 p_z=z_next(p_z) endloop,節(jié)點(diǎn)遍歷程序框架： p_gp=gp_head loop while p

52、_z # null ;語(yǔ)句 p_gp=z_next(p_gp) endloop,111,實(shí)操應(yīng)用,3 FLAC3D數(shù)值模擬計(jì)算實(shí)例,112,實(shí)操應(yīng)用,3.FLAC3D數(shù)值模擬計(jì)算實(shí)例,洞室模型的生成 基坑開(kāi)挖過(guò)程模擬 路基施工過(guò)程模擬 荷載引起地基土體的超孔隙水壓力模擬 邊坡安全系數(shù)求解 隧道開(kāi)挖及支護(hù)模擬 礦體開(kāi)挖回填模擬 具有自由水面的穩(wěn)定流,113,實(shí)操應(yīng)用,洞室模型的生成,3.FLAC3D數(shù)值模擬計(jì)算實(shí)例,114,實(shí)操應(yīng)用,洞室模型的生成,3.FLAC3D數(shù)值模擬計(jì)算實(shí)例,gen zone radcylinder size 5 20 20 15 p0 0 0 615 p1 add 30

53、 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 1 gen 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

54、 1.25 1 1 group 2 gen 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 0,115,實(shí)操應(yīng)用,洞室模型的生成,3.FLAC3D數(shù)值模擬計(jì)算實(shí)例,116,實(shí)操應(yīng)用,基坑開(kāi)挖過(guò)程模擬,力學(xué)參數(shù),計(jì)算模型及監(jiān)測(cè)點(diǎn)布置,3.FL

55、AC3D數(shù)值模擬計(jì)算實(shí)例,117,實(shí)操應(yīng)用,基坑開(kāi)挖過(guò)程模擬,各開(kāi)挖步下最大主應(yīng)力分布,3.FLAC3D數(shù)值模擬計(jì)算實(shí)例,118,實(shí)操應(yīng)用,基坑開(kāi)挖過(guò)程模擬,監(jiān)測(cè)節(jié)點(diǎn)位移及總體位移響應(yīng),3.FLAC3D數(shù)值模擬計(jì)算實(shí)例,119,實(shí)操應(yīng)用,路基施工過(guò)程模擬,幾何邊界,計(jì)算模型,土層物理力學(xué)參數(shù),A,B,3.FLAC3D數(shù)值模擬計(jì)算實(shí)例,120,實(shí)操應(yīng)用,路基施工過(guò)程模擬,路基回填過(guò)程中鉛直向應(yīng)力云圖,3.FLAC3D數(shù)值模擬計(jì)算實(shí)例,121,實(shí)操應(yīng)用,路基施工過(guò)程模擬,監(jiān)測(cè)點(diǎn)z向位移,監(jiān)測(cè)點(diǎn)z向位移,填筑結(jié)束時(shí)沉降云圖,填筑結(jié)束時(shí)水平位移云圖,3.FLAC3D數(shù)值模擬計(jì)算實(shí)例,122,實(shí)操應(yīng)用,荷載引起地基土體的超孔隙水壓力模擬,力學(xué)參數(shù),幾何邊界,3.FLAC3D數(shù)值模擬計(jì)算實(shí)例,123,實(shí)操應(yīng)用,荷載引起地基土體的超孔隙水壓力模擬 緩慢加載 def ramp stress = min(1.0,float(step)/200.0) end ramp apply nstress = -40e3