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

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

2、單元(liner),FLAC3D中的結(jié)構(gòu)單元是巖土工程中實際結(jié)構(gòu)的一種“抽象”，即采用簡單的單元形式來模擬復(fù)雜的結(jié)構(gòu)體。 結(jié)構(gòu)單元由結(jié)構(gòu)節(jié)點(node)和結(jié)構(gòu)構(gòu)件(SELs)構(gòu)成。 結(jié)構(gòu)單元中的節(jié)點(node)可以與周圍的實體網(wǎng)格(zone)或其它結(jié)構(gòu)節(jié)點建立連接(link)，通過連接實現(xiàn)巖土體或結(jié)構(gòu)與其它結(jié)構(gòu)發(fā)生相互作用。 注意：結(jié)構(gòu)節(jié)點并不是簡單地與實體網(wǎng)格的節(jié)點(gridpoint)建立聯(lián)系，也不能建立node與gridpoint之間的link,1、結(jié)構(gòu)單元的類型,4,實操應(yīng)用,梁單元 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)點，第二種方式適合建立復(fù)雜曲線結(jié)構(gòu)單元（但是要注意它不會自動建立link！若不手動link就無任何作用）,5,實操應(yīng)用

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

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

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

7、結(jié)構(gòu)單元建模的基本原則。,2、結(jié)構(gòu)單元的建模方法注意事項,13,實操應(yīng)用,梁單元 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單位長度上水泥漿粘結(jié)力cg gr_fric水泥漿的摩擦角g gr_k單位長度上水泥漿剛度kg gr_per水泥漿外圈周長Pg slide大變形滑動標志 slide_tol大變形滑動容差 yco

8、mp抗壓強度(力) density密度 thexp熱膨脹系數(shù),3、結(jié)構(gòu)單元的參數(shù)取值,14,實操應(yīng)用,3、結(jié)構(gòu)單元的參數(shù)取值,某些結(jié)構(gòu)單元參數(shù)的取值要視具體情況而定，根據(jù)經(jīng)驗且必要時調(diào)整參數(shù)通過試算來確定。,15,實操應(yīng)用,4、結(jié)構(gòu)單元實例分析,4.1、簡支梁（beam單元）承受兩個相等集中載荷 4.2、簡支梁（shell單元）承受兩個相等集中載荷,16,實操應(yīng)用,Simple Beam Two Equal Concentrated Loads,4.1、簡支梁（beam單元）承受兩個相等集中載荷,17,實操應(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、簡支梁（beam單元）承受兩個相等集中載荷,18,實操應(yīng)用,根據(jù)理論公式計算得到：,4.1、簡支梁（beam單元）承受兩個相等集中載荷,19,實操應(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、簡支梁（beam單元）承受兩個相等集中載荷,20,實操應(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、簡支梁（beam單元）承受兩個相等集中載荷,如何設(shè)置結(jié)構(gòu)單元的跟蹤變量！,21,實操應(yīng)用,Beam_concent_loads_Example1.3,4.1、簡支梁（beam單元）承受兩個相等集中載荷,22,實操應(yīng)用,撓度計算,4.1、簡支梁（beam單元）承受兩個相等集中載荷,23,實操應(yīng)用,剪力、彎矩計算,24,實操應(yīng)用,剪力、彎矩計算,這是節(jié)點力！,25,實

14、操應(yīng)用,26,實操應(yīng)用,梁單元局部坐標系：x軸從節(jié)點1到節(jié)點2，y軸在橫截面中,4.1、簡支梁（beam單元）承受兩個相等集中載荷,27,實操應(yīng)用,彎矩矢量的指向，右手法則！,4.1、簡支梁（beam單元）承受兩個相等集中載荷,28,實操應(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é)：梁單元的常用命令,4.1、簡支梁（

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

17、個相等集中載荷,33,實操應(yīng)用,4.2、簡支梁（shell單元）承受兩個相等集中載荷,34,實操應(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、簡支梁（shell單元）承受兩個相等集中載荷,35,實操應(yīng)用,局部坐標系！看彎矩到底應(yīng)該看哪一個？,36,實操應(yīng)用,4.2、簡支梁（shell單元）承受兩個相等集中載荷,37,實操應(yīng)用,4.2、簡支梁（shel

18、l單元）承受兩個相等集中載荷,38,實操應(yīng)用,39,實操應(yīng)用,4.2、簡支梁（shell單元）承受兩個相等集中載荷,40,實操應(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,實操應(yīng)用,5、關(guān)于link,兩種： Node-Zone Node-Node,42,實操應(yīng)用,創(chuàng)建一個新的link，link的源節(jié)點為sid，而聯(lián)接目標為node或zo

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

24、rigid”。 Side1，side2關(guān)鍵字對于確定嵌入式liner的哪個面上產(chǎn)生link。,sel,可選參數(shù)與必選參數(shù)！,43,實操應(yīng)用,5、關(guān)于link,以預(yù)應(yīng)力錨桿的托盤模擬為例 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,實操應(yīng)用,2、接觸面建立方法So-called “倒來倒去法”,我們最終的目的就是在中心小塊體與外圍網(wǎng)格之間建立接觸面。 分開建立網(wǎng)格 建立inner網(wǎng)格及其表面的Interface 導(dǎo)入外圍mesh 賦予材料屬性，測試接觸面是否發(fā)生了作用。,72,實操應(yīng)用,小練習(xí)： 三種建立接觸面的方法計算結(jié)果是否相同（只要接觸面有響應(yīng)，肯定是相同的?。?如果將接觸面建立在外部網(wǎng)格的內(nèi)表面，然后移入小塊體，結(jié)果是否相同呢？ 不加接觸面跟加了接觸面，模型的響應(yīng)（位移、應(yīng)力）有何區(qū)別？,Nrange Exa

26、mple-7.1 union nrange,73,實操應(yīng)用,3、切割模型的方法實際上是分離連續(xù)網(wǎng)格（原來網(wǎng)格連續(xù)，通過共用節(jié)點（GridPoint）傳遞力，分離后通過接觸面來傳遞。,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,實操應(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個類組上面，若不指定第2個類組，就默認第2個類組為整個模型，邊界面為整個第1個類組的表面。,75,實操應(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)該是一個共享面 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 產(chǎn)生由range定義的內(nèi)部 face。面上的節(jié)點復(fù)制，并創(chuàng)建

34、新的表面。,76,實操應(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,實操應(yīng)用,檢驗是否分離開了網(wǎng)格并形成了接觸面？,78,實操應(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,實操應(yīng)用,80,實操應(yīng)用,實例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,實操應(yīng)用,1. FLAC3D基本操作方法,82,實操應(yīng)用,1. FLAC3D基本操作方法,軟件界面,FLAC3D5.00 軟件界面,83,實操應(yīng)用,1. FLAC3D基本操作方法,軟件界面,FLAC3D5.00 軟件界面,84,實操應(yīng)用,1. FLAC3D基本操作方法,命令流操作方法,new gen zone brick size 6 6 6,在利用FLAC3D軟件進行數(shù)值模擬時，主要是通過命令流來實現(xiàn)的。命令流文件一般以txt或dat格式存儲，并在命令輸入窗口、菜單欄或快捷圖標通過call命令進行調(diào)用。,FLAC3D命令流文件需要遵循一定的格

42、式和語法要求，在滿足這些要求的前提下，命令流文件的編寫又是十分自由靈活的，可根據(jù)用戶個人的需求自由編寫。,85,實操應(yīng)用,1. FLAC3D基本操作方法,數(shù)值計算一般流程,建立模型 通過外部導(dǎo)入或在FLAC3D中直接建模的方式建立計算模型。 材料參數(shù) 材料本構(gòu) 材料力學(xué)參數(shù) 邊界條件及初始條件 速度邊界 應(yīng)力邊界 水頭邊界,86,實操應(yīng)用,1. FLAC3D基本操作方法,建立模型 在FLAC3D中，模型的建立通過關(guān)鍵詞generate來實現(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,實操應(yīng)用,1. FLAC3D基本操作方法,建立模型 FLAC3D建模的基本思路為“堆積木”，即首先建立各種形狀的網(wǎng)格單元，最后將建立的網(wǎng)格單元組合在一起，生成可用于數(shù)值計算的整體模型。 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,實操應(yīng)用,1. FLAC3D基本操作方法,網(wǎng)格形狀,89,實操應(yīng)用,1. FLAC3D基本操作方法,brick,gen zone brick size 6 8 8 plot zone,90,實操應(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,實操應(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,實操應(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,實操應(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,實操應(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,實操應(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,實操應(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,實操應(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,實操應(yīng)用,1. FLAC3D基本操作方法,通過對稱生成網(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,實操應(yīng)用,1. FLAC3D基本操作方法,生成漸變網(wǎng)格,gen zone brick size 10 10 10 p0 0 0 0 定義函數(shù)名

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

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

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

52、_z # null ;語句 p_gp=z_next(p_gp) endloop,111,實操應(yīng)用,3 FLAC3D數(shù)值模擬計算實例,112,實操應(yīng)用,3.FLAC3D數(shù)值模擬計算實例,洞室模型的生成 基坑開挖過程模擬 路基施工過程模擬 荷載引起地基土體的超孔隙水壓力模擬 邊坡安全系數(shù)求解 隧道開挖及支護模擬 礦體開挖回填模擬 具有自由水面的穩(wěn)定流,113,實操應(yīng)用,洞室模型的生成,3.FLAC3D數(shù)值模擬計算實例,114,實操應(yīng)用,洞室模型的生成,3.FLAC3D數(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,實操應(yīng)用,洞室模型的生成,3.FLAC3D數(shù)值模擬計算實例,116,實操應(yīng)用,基坑開挖過程模擬,力學(xué)參數(shù),計算模型及監(jiān)測點布置,3.FL

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