FLAC3D軟件的案例
flac3d軟件批量計算
關于abaqus等軟件的批量計算問題，網(wǎng)上資料較多，而針對flac3d軟件的批處理計算網(wǎng)上卻鮮有報道。近日略有空閑，就出一期關于flac3d軟件的批量計算，以應對大量的參數(shù)化計算。話不多說，直接上干貨。
下面給出3種批量計算的方法，itasca系列其他軟件也可以參考實現(xiàn)。
（1）windows系統(tǒng)批處理的方法
（2）flac3d內(nèi)置python方法
（3）flac3d內(nèi)置fish方法
上述三種批量計算的方法，本人在7.0版本上實測均可運行，如有需要文中代碼和模型文件，可私信，需要的人多的話，免費公開。
展開 基于FLAC3D軟件的人工凍結(jié)溫度場分析
FLAC3D軟件作為有限差分軟件，除廣泛應用于巖土體的力學問題分析外，還可以用于凍結(jié)溫度場的分析。本文后續(xù)演示了采用人工地下凍結(jié)法進行煤礦井筒開挖的分析代碼，是我在碩士階段編寫的，雖然代碼還有不少值得完善的地方。但里面的包含了諸如坐標數(shù)據(jù)的讀入、自動截圖和數(shù)據(jù)文件的自動取名保存等等。 希望對大家學習FLAC3D有幫助。
凍結(jié)方式為豎向凍結(jié)，冷媒是鹽水。
圖1 人工凍結(jié)法凍結(jié)壁形成示意圖
代碼：
Exported from Notepad++
; ----------------------------------------
; Description: 豎井溫度場分析
; Author: Qiao Cheng
; ----------------------------------------
new
config thermal
; --------------------
; 定義幾何模型參數(shù)
; --------------------
def Geocons ; geometrical constants
eps = 1.e-3
meps = -1 * eps
;z coordinate of top
z1 = 1.
z1m = z1 - eps
z1p = z1 + eps
;
out_b = 21.
展開 DxfToAnsys下載
在進行巖土方面的數(shù)值模擬時，一般采用FLAC3D軟件進行計算分析。對于復雜的模型建立，FLAC3d很難做到，因此，需要應用其他軟件進行輔助建模(cad-ansys-flac)。具體的步驟如下所示：
1、在cad里進行畫出模型的平面圖：
2、cad導入ansys軟件時有兩種方法：
(1)借用軟件dxftoansys，在cad里面將文件另存為XXX.dxf
用dxftoansys軟件打開XXX.dxf，選擇折線的形式。確定后會出現(xiàn)一些命令流，復制這些命令流，粘貼到ansys命令流框里。之后ansys會顯示圖形出來。ansys顯示的所有的線，將線畫成面即可進行網(wǎng)格的劃分。
(2)不借用軟件dxftoansys，在cad里將圖形建立成面域，將圖形輸出，輸出的后綴名為XXX.sat。打開ansys軟件，點擊file>import>.sat。即可打開XXX.sat,是個面的模型，點擊前視圖，后面進行網(wǎng)格劃分。
3、ansys進行網(wǎng)格劃分
再進行網(wǎng)格劃分時候，首先要定義模型屬性 點擊preferences>點擊structure，展開preprocessor>element type>add/edit/delete>add 輸入plane42點擊apply輸入solid45點擊close。
展開 Itasca軟件(FLAC3D, 3DEC, PFC)推出Linux版本---一個更好的體驗?
在原來的FISH語句中，下面這段代碼是沒有問題的：
fish define in_def(msg,default) local xx = io.in(msg+'('+'default:'+string(default)+'):') if type(xx) = 3 in_def = default else in_def = xx endifend
但不知什么時候軟件作了改變，導致運行這段程序時出現(xiàn)錯誤：“Cannot use operator + between types Float and String.” 意思是浮點數(shù)和字符串不可以用+號相連接， 這意味著不能使用string()函數(shù)顯式地把數(shù)字轉(zhuǎn)換為字符串, 必須使用下面的string.build函數(shù)才不出錯。string.build()函數(shù)類似于Python的format()定位函數(shù)，技術上講，這種改變是一種進步。
xx = io.in(string.build('%1 (default:%2):',msg,default))
測試顯示在3DEC和FLAC3D中都存在這個問題，但在UDEC中還未作出改變，string()和string.build()都可以使用，整合時期感覺弄得挺混亂的，估計以后UDEC也會作出相應的改變。
展開 
要學好FLAC3D，必須學好FISH
要學好FLAC3D，必須學好FISH，F(xiàn)ISH身為FLAC3D的內(nèi)嵌語言，也就是是一種嵌入式編程語言，其編寫短小精干實時，同時兼顧大家操作習慣，不但可以嵌入命令流文件里工作，而且還可以引用FLAC3D本身的任何命令，突破了一般標準程序代碼的限制，實現(xiàn)了讀者對FLAC3D軟件的完美控制。
FISH手冊.pdf
【4月2日項目懸賞】
【單號5268】
預算范圍：2000-5000
使用軟件： Abaqus
需求描述：開發(fā)Lemaitre損傷的umat程序
立即搶單
【單號5292】
預算范圍：800
使用軟件：Flac3D
需求描述：用Flac3D軟件，采用CRD工法做隧道的開挖。比較簡單的開挖
立即搶單
【單號5356】
預算范圍：500-1000
使用軟件：ABAQUS
需求描述：需要一位熟練使用ABAQUS的工程師在線幫助，做傳熱及振動方面的，由于項目涉密，只能我自己做，沒做過熱振分析相關內(nèi)容，故需要在線指導.
立即搶單
【單號5337】
預算范圍：4000+
使用軟件：Hypermesh/Abaqus
需求描述：空氣彈簧強度剛度計算 1.附件的Excel表里幾種工況的靜剛度和動剛度。 2.找到裝配壓力和拉脫力的關系。 3.提議列出來的三種材料，那種防滑線效果最好。
展開 【3月14日項目懸賞】
（查看附件請點擊“立即搶單”）
立即搶單
【單號5242】
預算范圍：1000
使用軟件：Fluent
需求描述：基于未反應核收縮模型的液固兩相反應（A+B+C+D=E+F，其中D、F 是固相，其余為液相）
立即搶單
【單號5257】
預算范圍：3000
使用軟件：Fluent
需求描述：水上系泊浮體在波浪作用下的動力特性分析
立即搶單
【單號5268】
預算范圍：2000-5000
使用軟件：Abaqus
需求描述：開發(fā)Lemaitre損傷的umat程序
立即搶單
【單號5273】
預算范圍：2000
使用軟件：abaqus
需求描述：使用abaqus模擬出該圖上的曲線（查看圖片請點擊“立即搶單”）
立即搶單
【單號5292】
預算范圍：800
使用軟件：Flac3D
需求描述：用Flac3D軟件，采用CRD工法做隧道的開挖。比較簡單的開挖
立即搶單
【單號5303】
預算范圍：100~1000
使用軟件： autodyn ;abaqus
需求描述：如何將autodyn計算后的模型導出來，供其他軟件如abaqus進行后續(xù)的計算？
展開 一個flac3d隧道數(shù)值計算前后處理全過程
15、運行“flac3d隧道后處理助手”程序，打開“flac3dmn.dat”，進行相關繪圖設置后，繪制內(nèi)力圖；或者生成cad腳本文件（*.scr），在cad中輸入命令script，打開相應的腳本文件，生成內(nèi)力圖。（注：在cad中輸入命令script前，必須關閉“對象捕捉，正交，對象追蹤等輔助做圖選項，否則會造成內(nèi)力圖形混亂。
16、若需要對圍巖應力場和位移場等進行后處理，比如做應力場和位移場等值線等，可使用dynamax編寫的“Flac3d2Tecplot.dat”，導出數(shù)據(jù)，在tecplot中進行處理。
基于FLAC3D軟件求解隧道襯砌強度安全系數(shù)的方法.rar
fish及實例.rar
展開 【3月21日項目懸賞】
（查看附件請點擊“立即搶單”）
立即搶單
【單號5242】
預算范圍：1000
使用軟件：Fluent
需求描述：基于未反應核收縮模型的液固兩相反應（A+B+C+D=E+F，其中D、F 是固相，其余為液相）
立即搶單
【單號5257】
預算范圍：3000
使用軟件：Fluent
需求描述：水上系泊浮體在波浪作用下的動力特性分析
立即搶單
【單號5268】
預算范圍：2000-5000
使用軟件：Abaqus
需求描述：開發(fā)Lemaitre損傷的umat程序
立即搶單
【單號5273】
預算范圍：2000
使用軟件：abaqus
需求描述：使用abaqus模擬出該圖上的曲線（查看圖片請點擊“立即搶單”）
立即搶單
【單號5292】
預算范圍：800
使用軟件：Flac3D
需求描述：用Flac3D軟件，采用CRD工法做隧道的開挖。比較簡單的開挖
立即搶單
【單號5303】
預算范圍：100~1000
使用軟件： autodyn ;abaqus
需求描述：如何將autodyn計算后的模型導出來，供其他軟件如abaqus進行后續(xù)的計算？
展開 基于FLAC3D的地下硐室通道錨網(wǎng)噴支護分析
目前絕大多數(shù)數(shù)值計算軟件并沒有直接的相應內(nèi)置單元來模擬錨桿，而FLAC3D內(nèi)置了cable單元用于模擬錨桿或錨索支護結(jié)構(gòu)，同時其強大的Fish語言允許用戶通過編程實現(xiàn)大量錨桿的自動創(chuàng)建、連接的建立、屬性參數(shù)的設置等工作，代碼一旦編寫完成只需少量修改即可反復使用，可以分析不同支護參數(shù)下的支護效果以便進行參數(shù)優(yōu)化。
本文作者之前在國內(nèi)兩家著名仿真論壇公布了錨網(wǎng)噴支護的部分代碼，獲得了大量的點擊和關注，當時公布的代碼已經(jīng)進行了簡化，并未公開其中錨桿與噴層混凝土間連接的相關操作。在本文中，將較詳細的介紹該結(jié)構(gòu)單元節(jié)點上的連接的自動刪除和新建方法，并公開相應的核心計算代碼，該部分內(nèi)容對于在模擬中獲得真實的支護效果非常重要。
本文在接下來的第1節(jié)將介紹FLAC3D中結(jié)構(gòu)單元的部分基礎知識，尤其是單元間連接建立的一般方法，然后介紹在需要建立大量連接時，如何利用內(nèi)置函數(shù)實現(xiàn)連接的建立。在第二節(jié)將介紹一個實例，在實例中具體演示連接的自動建立。第三節(jié)將給出該實例對應的完整代碼。有了前面各節(jié)的鋪墊，讀者將更容易理解和掌握第三節(jié)的代碼。
1 FLAC3D中連接的建立
1.1 基礎介紹
FLAC3D中的結(jié)構(gòu)單元主要包括：錨桿（錨索）單元(cable)、梁單元(beam)、殼體單元(shell)、襯砌單元(liner)和土工格柵(geogrid)。這些結(jié)構(gòu)單元一般在普通實體單元(zone)表面（圖3）或內(nèi)部建立（圖4），結(jié)構(gòu)單元在創(chuàng)建時單元節(jié)點(node)會自動與鄰近的實體單元建立連接（link）以實現(xiàn)與周圍實體單元的相互作用。結(jié)構(gòu)單元的單元與單元之間也通過節(jié)點(node)上的link實現(xiàn)共同作用。因此，link有兩種類型，即：節(jié)點-實體單元(node-zone)（圖4）、節(jié)點-節(jié)點（node-node）。
展開 基于FLAC3D的地下硐室通道錨網(wǎng)噴支護 附FLAC-3D深基坑的開挖與支護的命令流下載
下載地址：FLAC-3D深基坑的開挖與支護的命令流
降雨強度及持時對邊坡穩(wěn)定性影響研究
趙衡等[5]利用FLAC3D軟件對某路塹邊坡進行數(shù)值模擬分析，得出邊坡破壞方式為對稱破壞，并提出斜坡穩(wěn)定性極限平衡計算方法。喬翔等[6]針對某公路邊坡的剖面模型，采用極限平衡法對坡體不同部位進行穩(wěn)定性分析，并根據(jù)受力分析提出合理的邊坡加固方案。劉勇等[7]以改良的極限平衡法為基礎，結(jié)合室內(nèi)測試和數(shù)值模擬技術，揭示降雨對邊坡安全系數(shù)產(chǎn)生顯著影響的影響因素，如降雨強度和降雨時長等，并計算了在降雨入滲作用下，非飽和土質(zhì)邊坡的穩(wěn)定性。
大量邊坡失穩(wěn)工程案例表明，邊坡土體在遭受降雨入滲作用時很容易出現(xiàn)失穩(wěn)破壞。雨水的滲透會使坡面土體飽和度增加，導致其抗剪能力下降，嚴重影響邊坡穩(wěn)定性。本文以某水庫邊坡為研究對象，分析降雨強度和持時對邊坡安全系數(shù)和變形滲流特征的影響，為降雨條件下邊坡的支護設計方案提供參考。
1 工程概況
某水庫邊坡模型示意圖如圖1所示，坡面前緣長10m, 后緣長20m, 橫向投影23m。邊坡主要巖性為坡積土，由泥巖露頭和早期淺成巖等巖性組成含少量碎石，碎石粒徑在5～15mm之間，含量約30%,坡積土下部為花崗巖基巖，厚約15～25m。斜坡土體的物理和力學參數(shù)見表1。根據(jù)斜坡幾何尺寸和物理力學參數(shù)，利用ABAQUS軟件對對邊坡滲流及穩(wěn)定性的影響進行分析。
圖1 邊坡模型示意圖
表1 斜坡面土體物理力學參數(shù)
2 邊坡降雨入滲數(shù)值模擬
2.1 降雨強度對邊坡滲流的影響
在ABAQUS數(shù)值模擬中，降雨強度可以表示為單位流量q(m/h)。通過設置降雨邊界條件實現(xiàn)降雨入滲模擬。斜坡前后及地面采用固定約束邊界，以限制邊坡在該方向上的位移。邊坡除邊界位移條件限制外，還需設置孔隙水壓邊界條件。將邊坡前緣孔壓設置為零，降水滲透基底處孔隙水壓不為零，透水界線設置為離坡底位置10米處，其它區(qū)域設置為不透水界線。
展開 基于FLAC3D的雙孔隧道圍巖穩(wěn)定性數(shù)值模擬分析
摘 要：針對雙孔隧道在施工過程中產(chǎn)生的復雜應力場與位移場，通過有限差分軟件FLAC3D數(shù)值模擬，采用Drucker-Prager準則作為巖土體塑性屈服準則，計算得到雙孔隧道開挖及支護后隧道的應力場及位移場分布規(guī)律。進而分析得到了雙孔隧道需要支護的關鍵位置及錨桿最小設計長度要求。結(jié)果表明：采用噴錨支護能夠有效阻止塑性區(qū)的擴大，對提高圍巖承載能力影響顯著。
關鍵詞：FLAC3D;噴錨支護;雙孔隧道;巖土工程;圍巖;
雙孔隧道作為現(xiàn)在主要的公路隧道形式，由于其獨特的構(gòu)造以及復雜的受力形式，成為了現(xiàn)在隧道建設者研究的重點與難點，尤其當兩個隧道間距較小時，在圍巖上覆荷載與支護反力共同作用下就會形成更加復雜的應力場。
目前隧道穩(wěn)定類的研究方法主要有解析法和數(shù)值模擬兩大類。對于隧道穩(wěn)定類問題，有大量學者對此點進行了相關研究。趙明華等[1]基于有限元極限分析雙孔隧道的穩(wěn)定性，得出影響隧道穩(wěn)定的相關因素。李揚等[2]基于有限元軟件Midas對淺埋雙孔隧道開挖順序進行了相關研究。但是他們都是采用有限元軟件進行了相關研究工作。也有學者采用有限差分軟件[3,4,5,6,7,8]進行了相關研究，但是他們都是研究單孔隧道或者煤礦等，在開挖及支護后圍巖的應力及位移分布。采用有限差分軟件FLAC3D進行雙孔隧道開挖及支護研究的則相對更少。
因此本文采用有限差分軟件FLAC3D對雙孔隧道在開挖和支護兩種工況下，進行了相關的數(shù)值模擬，同時分析了開挖后以及采用錨噴支護后隧道塑性區(qū)分布特點，縱向應力場、位移場分布規(guī)律。
1 圓形隧道理論模型
關于圍巖塑性區(qū)發(fā)展理論，基于Mohr-Coulomb準則計算的Fenner公式比較經(jīng)典，但其未考慮中間主應力對于巖石強度的影響。而中間主應力對巖石強度的影響程度為20%～50%[9],顯然不考慮中間主應力對于計算結(jié)果是有影響的。
展開 基于流固耦合的水底隧道仰拱受力分析與優(yōu)化
本文以廈門東通道海底隧道為工程背景,采用三維有限差分軟件FLAC3D進行流固耦合分析,研究仰拱曲率半徑及仰拱部位的防排水措施對圍巖位移及襯砌結(jié)構(gòu)受力的影響,優(yōu)化水底隧道斷面形式,以指導大型跨江海的水底隧道的設計和施工。
基于流固耦合的水底隧道仰拱受力分析與優(yōu)化.pdf
『原創(chuàng)』基于FLAC3D的連拱隧道開挖與支護
基于FLAC3D的連拱隧道開挖與支護
摘 要：本文在探討連拱隧道兩種常見的施工方法的基礎上，利用有限差分軟件FLAC3D對某處工程問題所引起地表沉降進行了數(shù)值模擬，編寫了自動開挖數(shù)據(jù)記錄程序，實現(xiàn)了連拱隧道全自動開挖和數(shù)據(jù)記錄，得到了較好的結(jié)果，同時也給工程實際提供了有價值的參考。
關鍵詞：連拱隧道 數(shù)值模擬 有限差分 開挖與支護
論文已經(jīng)發(fā)表
模型如下
