flac3d常見問題總結(jié)

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

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

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

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

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

6、ange annulus center xyz（中心坐標(biāo)）radius r1 r2（兩半徑）想要表達(dá)的是球面時(shí)，r1r2取值相差小點(diǎn)即可，如radius 2.9 3.1 這時(shí)范圍就是半徑為3的面。2。圓柱的環(huán)面施加力或約束時(shí)也用到該命令。命令range cylinder end1 xyz end2 xyz radius r(end是表示柱體斷面的中心點(diǎn)坐標(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、時(shí)的線的垂線與y軸正向的夾角，順為正，說白話，就是面的指向方向面是朝那個(gè)方向的意思。dip面與xoy面的夾角，與z軸負(fù)向?yàn)檎?，origin表示面經(jīng)過xyz點(diǎn)range normal xyz 法線方向?yàn)閤yz的面range not 沒試用過（范圍外的東西，想想用到圓上取范圍非常好）現(xiàn)在回到建模上面說過建模一般只建一半然后用reflect命令（鏡像）鏡像首先得有個(gè)面吧，就目前我用的面都是xoy面yoz面。估計(jì)大家也很少用到斜面鏡像吧面的命令上面range里說過了這里就不講了命令：reflect normal xyz origin xyz （根據(jù)法線和過一點(diǎn)建的面）最常用至于實(shí)際使用過程中有人問道：

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

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

10、一半然后鏡像對吧尺寸，自己在草紙上畫好，點(diǎn)坐標(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。整體的建模加顯示，的命令就完了合一塊 建立一個(gè)新的txt 開頭加一個(gè)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 顯示最小值，圖上小于這個(gè)

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

15、rav 0 0，-9.81ini dens 1000設(shè)置重力就一定要設(shè)置密度不考慮重力就不要設(shè)計(jì)密度。兩者共存的，設(shè)計(jì)程序的時(shí)候一定要記的set mech force 50solve用SET命令設(shè)定了最大不平衡力的上限值為50N,即：計(jì)算的最大不平衡力一旦小于此值，我們便人為模型達(dá)到了平衡狀態(tài)。求解過程終止。設(shè)置這個(gè)目的是在一定精度要求下，加快計(jì)算（值越大，算的量越少，花的時(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 。如這個(gè)命令，顯示圖像的時(shí)候命令有很多，那么我們想顯示別的圖像，又不想吧這個(gè)弄好的圖像丟掉則用plot create name建立一個(gè)圖像，在建立新的圖像時(shí)plot create name1（新的名）顯示以前的圖像那么輸入如下命令：plot current name。另一種方法，輸入：plot show name這里有一點(diǎn)最好開始的時(shí)候?qū)⒕W(wǎng)格剖分的相對稀疏一些（10001500個(gè)

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

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

19、令顯示不同的色了）如何查看剖面上任一點(diǎn)的位移值？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初始化命令，如初始化計(jì)算體的應(yīng)力狀態(tài)等；apply邊界條件限制命令，如施加邊界的力、位移等約束等。initial的應(yīng)力狀態(tài)會隨計(jì)算過程的發(fā)生而發(fā)生改變，一般體力需要初始化，而apply施加的邊界條件不會發(fā)生變化如果想看模型在計(jì)算過程中的變化情況

20、Set movie avi step 1 file 1.avi；定義動畫，動畫記錄的頻率為1時(shí)步一頻率（計(jì)算量大的時(shí)候就沒必要1時(shí)步的記錄，可以10,100時(shí)步的記錄，減少電腦的使用）Movie start；動畫開始hist n= 5 ：N(Nstep)=5指定了每迭代計(jì)算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é)點(diǎn)Z方向的位移（這個(gè)用的多，查看收斂量時(shí)就是這個(gè)命令，上面兩個(gè)個(gè)人人認(rèn)為用處不大，也很少

21、用過）如何在程序中進(jìn)行暫停，并可恢復(fù)計(jì)算？答：在命令中加入pause命令，用continue進(jìn)行繼續(xù)個(gè)人認(rèn)為這個(gè)命令非常好用在我們分步求解中想得到某一個(gè)過程中的結(jié)果，不用等到全求完，還可以在分布求解錯(cuò)誤的時(shí)候就進(jìn)行改正。而不是等到結(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í)慣需要,例如樁載荷試驗(yàn)曲線,沉降曲線等都不是畫在第一象限內(nèi).最近看有很多對后處理中的收斂有很多不同的想法，自己也總結(jié)一下plot sk dis scale 0.07顯示位移矢量箭頭的后面的scale默認(rèn)的是0.053越大箭頭越大magf也是放大變形的意思，但是不能放在sk dis的后面，計(jì)算機(jī)會出錯(cuò)回到主題，收斂數(shù)據(jù)，我們想得到一條收斂量隨時(shí)間變化，或者不同應(yīng)力直間關(guān)系的圖線方法1用set log on然后把數(shù)據(jù)導(dǎo)入excel或者其他的處理軟件進(jìn)行處理2.用hist i

23、d=1（不寫id的話，電腦會自己按的順序進(jìn)行排序）gp（點(diǎn)，如果是體的話用zone）zdis xyz（xyz點(diǎn)處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ù)多計(jì)算文件又太大，浪費(fèi)空間。設(shè)置步長的命令set hist_rep n(n就是步長)上面說的p

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

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

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