flac3d常見問題總結(jié)

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

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

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

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

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

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

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

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

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

10、一半然后鏡像對吧尺寸，自己在草紙上畫好，點坐標(biāo)都弄清楚了把以前的句子直接套下來gen zone radcylinder size 4 25 6 14（網(wǎng)格數(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這兩行命令，用來顯示指定面上的位移變化云圖。（吧第二行改改，還可以顯示應(yīng)力云圖，等其他）plot con syy y向應(yīng)力云圖plot con szz outline on z向應(yīng)力云圖加顯示網(wǎng)格plot con szz effective on z向有效應(yīng)力plot con szz inter 1e4 云圖顯示增量（后面的值越小，圖上顏色條帶就越細(xì)，）plot con szz min -10e4 顯示最小值，圖上小于這個

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

15、rav 0 0，-9.81ini dens 1000設(shè)置重力就一定要設(shè)置密度不考慮重力就不要設(shè)計密度。兩者共存的，設(shè)計程序的時候一定要記的set mech force 50solve用SET命令設(shè)定了最大不平衡力的上限值為50N,即：計算的最大不平衡力一旦小于此值，我們便人為模型達(dá)到了平衡狀態(tài)。求解過程終止。設(shè)置這個目的是在一定精度要求下，加快計算（值越大，算的量越少，花的時間也越少）說到這里，add與set的區(qū)別add原來有的加新的用addset原來沒有的，設(shè)置用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這里有一點最好開始的時候?qū)⒕W(wǎng)格剖分的相對稀疏一些（10001500個

17、單元左右），先求解一下，看看結(jié)果，然后再根據(jù)大致的結(jié)果有針對性的增加模型的網(wǎng)格數(shù)，以得到更精確的解答。更快的檢驗一下自己 的模型是否建錯了不同的區(qū)域定義不同的本構(gòu)模型model mohr range z=0,5prop ClayeyGravel range z=0,5在定義參數(shù)后面加個range范圍為了避免在輸入命令時總要通過range全局坐標(biāo)域的方法來限定該命令執(zhí)行的區(qū)域，我們可以先通過range命令來給不同的區(qū)域命名，然后就可以用range區(qū)域名的方法來限定范圍了。 這樣就避免了每次輸入命令要限定在該區(qū)域執(zhí)行時，都要輸入一遍坐標(biāo)范圍的麻煩。具體如何給區(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的方向很多人說網(wǎng)格后的云圖有帶齒什么的，改進的方法就是重復(fù)用下面的命令Gen ADJust 用此命令來微調(diào)網(wǎng)格, 使之離散平滑。該命令可以連續(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很多人都是這么做的，我也是學(xué)習(xí)，initial 與 apply 有何區(qū)別？initial初始化命令，如初始化計算體的應(yīng)力狀態(tài)等；apply邊界條件限制命令，如施加邊界的力、位移等約束等。initial的應(yīng)力狀態(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次記錄一次相關(guān)的值（如果不指定，則系統(tǒng)默認(rèn)值為10）hist unbal 命令的第二行指定記錄最大不平衡力（the maximum unbalanced force）hist gp zdisp 4,4,8 記錄坐標(biāo)值為（4，4，8）的節(jié)點Z方向的位移（這個用的多，查看收斂量時就是這個命令，上面兩個個人人認(rèn)為用處不大，也很少

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

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

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

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

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

26、力的問題，我也有個疑問，平衡地應(yīng)力后，水平方向的地應(yīng)力如何確定的？是由本構(gòu)模型來確定的嗎？這個問題現(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來賦予初始地應(yīng)力,但是為何又有set gravity 0,0,-10這樣不是重復(fù)了嗎? ini是的確是模型的內(nèi)力，這個內(nèi)力通過set grav ,然后就像一些帖子上介紹的那樣,讓程序計算不就行了嗎?干嘛還有手動在設(shè)定一次呢?我感覺它這樣的做法有點重復(fù)了. 回答：（1）它這樣做其實只是減少程序的計算步驟.5 _這樣設(shè)置以后 整個模型就已經(jīng)是平衡了（2）設(shè)置重力的原因不是為了平衡減少程序計算步驟，而是因為如果未初始化應(yīng)力的話，在set grav的時候，模型會在自重應(yīng)力下發(fā)生應(yīng)力變形。所以，初始化應(yīng)力的原因是防止模型在自重應(yīng)力下發(fā)生變形（我們認(rèn)為地層自重變形已經(jīng)發(fā)生，當(dāng)然欠固結(jié)等現(xiàn)象是另外一回事）。3）其實就讓模型發(fā)生變形，然后把位移清零也可以的 6。ini命令來對地層進行初始應(yīng)力的施加，而APPLY是7 u(施加應(yīng)力邊界條件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,