flac3d常見問題總結

1、flac3d學習有一段時間，過一段時間沒學究忘了很多，所以學習最好不要放下來。對自己總結，也是學習的一種。flac3d命令很多，用到的不多，就那么幾個，基礎的就不多說了進入正題：個人認為如果學習flac3d最好把flac的命令進行總結一下。做模型就好比拼盤：要地應力，就把地應力的那段命令摘入，要出圖，加入plot的相關命令，要出動畫，加個動畫命令，所以總結好命令后，以后做模型，只要把命令組一條條的對應放入，一個程序就出來了，后期調試一下，改相關參數(shù)就可以了。接下來，我把自己弄好的整理一下，按照我上面的思路進行建模什么人建模都會有自己的喜好，下面說一下我自己建模第一步就是建立一個全新的txt文件

2、（不論命令大小都是如此），這樣的好處是，模型可以修改方便快捷。而且還可以自我查錯建立新的txt后，第一行命令就是new，加new的目的是為了防止其他命令的影響（相當于刷新）再執(zhí)行其他命令前刷新一下。不用想，接下來就是建模了，建立模型，不同人有不同的方法，有的人喜歡建完后在對模型開挖最終得到自己想要的模型。model null 等；有的人喜歡用各種帶缺角的組合成最終的模型。具體建模看各自的喜好，我覺得，不同點是建模的網格劃分有區(qū)別，用帶缺角的組合模型網格劃分更好最好把常用的集中模型記下來組成各種句子，以后用的時候直接插入例如：gen zone brick/radbrick size *,*,*,

3、* ratio *,*,*,*建模的過程中，或同一個文件中有重復應用的命令，就可以用宏命令了，宏命令macro 名稱 命令的內容例：macro shiping plot con zdip outline on用shiping 代替后面的引號里的東西。在建模中只要有相同的都可以用shiping代替，縮減了書寫命令的時間，整體上也能讓人看的更明白注意點是，你命名的符號最好跟falc自帶的命令不要重復，我的處理方法是，名稱后帶數(shù)字如shiping1如果你想修改宏對象，改變命令的內容，在名稱上加單引號就是替換的意思例macro shiping1 新的命令這些很多課本里都有講，不在敘述很多時候我們要對單

4、一的對象進行操作，例如進行模擬運算后，只想單獨的顯示一部分的內容，單一的土層，錨桿，等就要用到groupgroup給對象命名，例：在gen zone brick - size - group name（name隨便取，不要跟flac的已有的命令重合就行）這里的意思就是對這個體進行了命名接下來想把這一塊模型挖掉加model null range group namegroup的用法：1。group直接加到模型的后面相當于給整個模型命名2.group name range x xyy zz給一個范圍的體命名繼續(xù)學習，昨天說到了建模，我是搞巷道的，重點講巷道。巷道建模，用到brick，radbric

5、k，radtunnel，radcylinder這4個模型比較多，一般模型都是建一半，讓后用鏡像命令復制另一半，建模也是大家有各自不同的習慣，總的有兩種，一種用點法，就是吧模型的所有關鍵點都輸入進去，這種方法可以用來建立不規(guī)則的模型，弊端就是工作量大，不當在設點的時候考慮多，輸入的時候也很麻煩。第二種。按照模型只要取三個點確定xyz三周的方法，有點輸入簡單，弊端，只能建立規(guī)則的模型，像要建立，路堤模型，邊坡模型，就只能是按照數(shù)點的方法建模range的命令：書上有，但是我用自己的話說一下他的用途1：取圓球面的范圍，就是同一個中心，不同半徑之間夾的那一部分，用在給球面施加力，約束時用該命令，命令：r

6、ange annulus center xyz（中心坐標）radius r1 r2（兩半徑）想要表達的是球面時，r1r2取值相差小點即可，如radius 2.9 3.1 這時范圍就是半徑為3的面。2。圓柱的環(huán)面施加力或約束時也用到該命令。命令range cylinder end1 xyz end2 xyz radius r(end是表示柱體斷面的中心點坐標)range plane above 表示的是z大于0的上部結構range plane blow 表示的是z小于0的下部結構range dd a dip b origin xyz 表示的是走向角度a，傾向角度b的面dd表示，當面與xoy面相交

7、時的線的垂線與y軸正向的夾角，順為正，說白話，就是面的指向方向面是朝那個方向的意思。dip面與xoy面的夾角，與z軸負向為正，origin表示面經過xyz點range normal xyz 法線方向為xyz的面range not 沒試用過（范圍外的東西，想想用到圓上取范圍非常好）現(xiàn)在回到建模上面說過建模一般只建一半然后用reflect命令（鏡像）鏡像首先得有個面吧，就目前我用的面都是xoy面yoz面。估計大家也很少用到斜面鏡像吧面的命令上面range里說過了這里就不講了命令：reflect normal xyz origin xyz （根據法線和過一點建的面）最常用至于實際使用過程中有人問道：

8、normal -1 0 0 與normal 1 0 0的區(qū)別，我試過兩者的效果是一樣的，沒區(qū)別，雖然方向不同，但表達的是同一個面reflect dd a dip b origin xyz 這個就是用里建立斜面鏡像的，各個符號的意思上面講過，不重復了有一點注意了，鏡像的時候如果是模型鏡像的話，相當于在建立一個模型，故，在reflect的前面要加建模的關鍵字，gen zone reflect。先講點后期出圖的，然后總結一下（plot set rotation 20 0 30 設置觀測角度的plot set center auto 設置觀測中心點plot set magnification 1.5

9、設置放大倍數(shù)的plot show）顯示這一整部分就是我開頭說的要自己整理出來的句段了，等以后用到的時候直接插入，作用就是顯示模型的。任何模型的顯示都可以用他來表示。很好很強大。如果還要顯示其他功能就把相應的命令在插入上面那部分就行了plot set rotation 20 0 30 plot set center autoplot set magnification 1.5 plot add surface yellowplot add axesplot show這個就是我常用來查看模型的命令，大家也可以自己整理一分，總結，做個實例，給大家看看隧道，半圓直墻隧道分兩步一，建模，二，顯示建模 建

10、一半然后鏡像對吧尺寸，自己在草紙上畫好，點坐標都弄清楚了把以前的句子直接套下來gen zone radcylinder size 4 25 6 14（網格數(shù)量，自己看著辦）ratio1 1 1 1.2&p0 0 0 0 p1 50 0 0 p2 0 100 0 p3 0 0 50 dimension 2 2 2這是建上半部分gen zone rad tunnel size 4 25 6 14 ratio 1 1 1 1.2 &p0 0 0 0 p1 0 0 -50 p2 0 100 0 p3 50 0 0 dimension 1.8 2 2建立下半部部分建一半后，用鏡像命令也講過了直接插入，g

11、en zone reflect normal -1 0 0 origin 0 0 0最后一步就是顯示給大家看那一打包顯示直接插入plot set rotation 20 0 30 plot set center autoplot set magnification 1.5 plot add surface yellowplot add axesplot show。整體的建模加顯示，的命令就完了合一塊 建立一個新的txt 開頭加一個newnewgen zone radcylinder size 4 25 6 14 ratio1 1 1 1.2&p0 0 0 0 p1 50 0 0 p2 0 10

12、0 0 p3 0 0 50 dimension 2 2 2gen zone rad tunnel size 4 25 6 14 ratio 1 1 1 1.2 &p0 0 0 0 p1 0 0 -50 p2 0 100 0 p3 50 0 0 dimension 1.8 2 2gen zone reflect normal -1 0 0 origin 0 0 0plot set rotation 20 0 30 plot set center autoplot set magnification 1.5 plot add surface yellowplot add axesplot show

13、整體就這樣了，大家可以粘貼一下看看。plot命令大集合：顯示切片的命令plot set plane ori 0 0 0 norm 0 1 0plot con zd plane這兩行命令，用來顯示指定面上的位移變化云圖。（吧第二行改改，還可以顯示應力云圖，等其他）plot con syy y向應力云圖plot con szz outline on z向應力云圖加顯示網格plot con szz effective on z向有效應力plot con szz inter 1e4 云圖顯示增量（后面的值越小，圖上顏色條帶就越細，）plot con szz min -10e4 顯示最小值，圖上小于這個

14、數(shù)的值都用同一色表示plot con szz max 10e4 顯示最大值，圖上大于這個數(shù)的值都用同一色表示plot fip outside on 顯示力的箭頭建模完事后就要給模型賦值，賦參數(shù)了，關鍵詞，modelmodel mor（賦摩爾庫倫模型，想要什么模型改相應的詞）pro den 2500（密度）bulk2e8體積模量，shear 2e8 （剪切模量）co 2e6（內聚力）fric45（摩擦角）ten1e6 抗拉強度。range.(我們要賦值的范圍)給模型賦值就是hi這么簡單，難在參數(shù)的選取?，F(xiàn)在看邊界，如果是不動面，用fix x range。如果邊界帶速率，或位移用ini。set g

15、rav 0 0，-9.81ini dens 1000設置重力就一定要設置密度不考慮重力就不要設計密度。兩者共存的，設計程序的時候一定要記的set mech force 50solve用SET命令設定了最大不平衡力的上限值為50N,即：計算的最大不平衡力一旦小于此值，我們便人為模型達到了平衡狀態(tài)。求解過程終止。設置這個目的是在一定精度要求下，加快計算（值越大，算的量越少，花的時間也越少）說到這里，add與set的區(qū)別add原來有的加新的用addset原來沒有的，設置用setplot create GravVplot set plane dip=90 dd=0 origin=3,4,0plot s

16、et rot 15 0 20plot set center 2.5 4.2 4.0plot add bound behindplot add bcont szz planeplot add axesplot show 。如這個命令，顯示圖像的時候命令有很多，那么我們想顯示別的圖像，又不想吧這個弄好的圖像丟掉則用plot create name建立一個圖像，在建立新的圖像時plot create name1（新的名）顯示以前的圖像那么輸入如下命令：plot current name。另一種方法，輸入：plot show name這里有一點最好開始的時候將網格剖分的相對稀疏一些（10001500個

17、單元左右），先求解一下，看看結果，然后再根據大致的結果有針對性的增加模型的網格數(shù)，以得到更精確的解答。更快的檢驗一下自己 的模型是否建錯了不同的區(qū)域定義不同的本構模型model mohr range z=0,5prop ClayeyGravel range z=0,5在定義參數(shù)后面加個range范圍為了避免在輸入命令時總要通過range全局坐標域的方法來限定該命令執(zhí)行的區(qū)域，我們可以先通過range命令來給不同的區(qū)域命名，然后就可以用range區(qū)域名的方法來限定范圍了。 這樣就避免了每次輸入命令要限定在該區(qū)域執(zhí)行時，都要輸入一遍坐標范圍的麻煩。具體如何給區(qū)域命名，以及日后如何用區(qū)域名來限定該區(qū)

18、域，可以看下面這個例子：range name trench x=-1,1 y=0,4 z=-2,2model null range trench今天做方案，沒時間，就講幾點得了，1：右手法則，建立模型時，的點的排放時有規(guī)律的，方法如下，拿出右手伸直了4指合并指向就是p0到p1的方向掌心指向就是p0到p2的方向拇指與4指垂直的方向就是p0到p3的方向很多人說網格后的云圖有帶齒什么的，改進的方法就是重復用下面的命令Gen ADJust 用此命令來微調網格, 使之離散平滑。該命令可以連續(xù)使用以增進平滑效果的。plo blo gro可以看到不同的group的顏色分布（例不同的土層定義組后，就可以用此命

19、令顯示不同的色了）如何查看剖面上任一點的位移值？print gp disp range x() y( ) z( )或print gp disp range idplot gp disp range x() y( ) z( )或plot gp disp range id很多人都是這么做的，我也是學習，initial 與 apply 有何區(qū)別？initial初始化命令，如初始化計算體的應力狀態(tài)等；apply邊界條件限制命令，如施加邊界的力、位移等約束等。initial的應力狀態(tài)會隨計算過程的發(fā)生而發(fā)生改變，一般體力需要初始化，而apply施加的邊界條件不會發(fā)生變化如果想看模型在計算過程中的變化情況

20、Set movie avi step 1 file 1.avi；定義動畫，動畫記錄的頻率為1時步一頻率（計算量大的時候就沒必要1時步的記錄，可以10,100時步的記錄，減少電腦的使用）Movie start；動畫開始hist n= 5 ：N(Nstep)=5指定了每迭代計算5次記錄一次相關的值（如果不指定，則系統(tǒng)默認值為10）hist unbal 命令的第二行指定記錄最大不平衡力（the maximum unbalanced force）hist gp zdisp 4,4,8 記錄坐標值為（4，4，8）的節(jié)點Z方向的位移（這個用的多，查看收斂量時就是這個命令，上面兩個個人人認為用處不大，也很少

21、用過）如何在程序中進行暫停，并可恢復計算？答：在命令中加入pause命令，用continue進行繼續(xù)個人認為這個命令非常好用在我們分步求解中想得到某一個過程中的結果，不用等到全求完，還可以在分布求解錯誤的時候就進行改正。而不是等到結果出來關于dd、dip的理解：不同人有不同的理解，我的理解就是dip面與水平面的夾角dd是面與水平面的交線的垂線與y軸的夾角1) 在plot hist m vs n的形式里,m代表y軸,n代表x軸(不管m,n的正負);5 L0 D# _1 1 C- L2) -表示對其值作mirror,如果大家使用AUTOCAD的話就明白mirror的含義了,比如上面的 hist 1

22、 gp ydisp的值是0100,那么vs -1就變成-1000.以此類推.3) 為什么要這么做呢?主要是為了符合巖土工程的習慣需要,例如樁載荷試驗曲線,沉降曲線等都不是畫在第一象限內.最近看有很多對后處理中的收斂有很多不同的想法，自己也總結一下plot sk dis scale 0.07顯示位移矢量箭頭的后面的scale默認的是0.053越大箭頭越大magf也是放大變形的意思，但是不能放在sk dis的后面，計算機會出錯回到主題，收斂數(shù)據，我們想得到一條收斂量隨時間變化，或者不同應力直間關系的圖線方法1用set log on然后把數(shù)據導入excel或者其他的處理軟件進行處理2.用hist i

23、d=1（不寫id的話，電腦會自己按的順序進行排序）gp（點，如果是體的話用zone）zdis xyz（xyz點處z向位移）這里應力應變之類的都是用zone，位移之類的都是用gp最后相要出相應的結果用命令plot hist 1；如果相要的不同的id之間的關系（比如說豎向應力szz與剪應力sxz的關系）命令plot hist 2 v 3 （2,3是對應應力開始設的id號）flac3d程序里面自動設置的步長是10，也就是說每隔10步記錄一次，這樣如果模型小得出的結果就數(shù)據少圖線不光滑，如果模型大，得到的數(shù)據多計算文件又太大，浪費空間。設置步長的命令set hist_rep n(n就是步長)上面說的p

24、lot多是flac里顯示的如果想把結果存到別處去用命令hist write 2 v 3 file namehist.txt(其實2相當于x軸，3相當于y軸)設置影片格式上面見過，溫故而知新set movie avi step 1 （存儲步長，可以任意設，就是每隔幾步存儲的意思）file name.avimovie start分步開挖就是分布計算每多一步（多添一塊土，或多開挖，多掘進等）就計算一下，加一個solve即可如果想保存每一步的結果和圖，那就在每一步的solve后面添加保存命令設立相應的文件，以后調用即可1.設定初始地應力場是十分必要的 沒有地應力場進行計算會使模型在重力作用下發(fā)生過大的

25、變形 從而不能與實際情況相符2。apply是施加的外力 ini是初始的內力。3. 有人問如果改參數(shù)這樣做的話地應力不是被人為改變了嗎？為什么要人為加大它的參數(shù)呢，很不合理？回答：不能改變bulk和shear模量，因為這兩個參數(shù)是用來計算應力的。而con和fric是用來判斷破壞狀態(tài)的。如果將bulk和shear模量改變了，那你在重力條件下計算的應力就不真實了。咱們改大參數(shù)的目的是讓模型在不破壞的條件下得到正確的應力，所以只有將con和fric變大才能達到目的。 當然如果在彈性模型下計算初始應力就沒有這種問題了。因為彈性模型是沒有破壞準側的當然也沒有描述達到破壞狀態(tài)的參數(shù)了。4. 大家正好討論地應

26、力的問題，我也有個疑問，平衡地應力后，水平方向的地應力如何確定的？是由本構模型來確定的嗎？這個問題現(xiàn)在我也沒答案也是我的問題5. gen zone brick size 10,10,10 p1 20,0,0 p2 0,20,0 p3 0,0,20model mohrprop bulk 5e9 shear 3e9 fric 35ini density 2500 set gravity 0,0,-10fix x y z range z -0.1 0.1ini szz = -5.0e6 grad 0,0,2.5e4ini sxx = -2.5e6 grad 0,0,1.25e4ini syy = -

27、2.5e6 grad 0,0,1.25eapply szz = -4.5e6 range z 19.9 20.1apply szz = -5.0e6 range z -0.1 0.17apply 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.1 U

28、# P! y- P$ _+ D它的意思是通過ini szz來賦予初始地應力,但是為何又有set gravity 0,0,-10這樣不是重復了嗎? ini是的確是模型的內力，這個內力通過set grav ,然后就像一些帖子上介紹的那樣,讓程序計算不就行了嗎?干嘛還有手動在設定一次呢?我感覺它這樣的做法有點重復了. 回答：（1）它這樣做其實只是減少程序的計算步驟.5 _這樣設置以后 整個模型就已經是平衡了（2）設置重力的原因不是為了平衡減少程序計算步驟，而是因為如果未初始化應力的話，在set grav的時候，模型會在自重應力下發(fā)生應力變形。所以，初始化應力的原因是防止模型在自重應力下發(fā)生變形（我們認為地層自重變形已經發(fā)生，當然欠固結等現(xiàn)象是另外一回事）。3）其實就讓模型發(fā)生變形，然后把位移清零也可以的 6。ini命令來對地層進行初始應力的施加，而APPLY是7 u(施加應力邊界條件7. 用到ini初始平衡的時候要和重力計算的相同才有用8。gen zone brick size 10 10 10 &p0 0,-25,0 p1 20,-25,0 p2 0,0,0 p3 0,-25,