Flac3D7.0
關注創建者：Flac3D程小杰 創建時間：2022-10-10
Flac3D7.0的實例教程
Msmin取值區間一般為[0.05，0.1]。
圖2 土體動剪切模量退化行為截斷處理
在V7.0版本中，FLAC3D通過引入如下命令及選項來反映對土體動剪切模量退化行為的截斷處理：
zone dynamic damping hysteretic keyword
keyword：
reduction-minimum：模量退化系數Ms的截斷值。
波動信號處理
動力分析中因分析方法的不同通常對待輸入波動信號的成分組成有規定要求，因此在分析之初有必要依據頻譜分析方法對信號做預處理，處理方法以濾波和基線修正為主：
濾波：目前動力分析方法主要包括頻域法和時域法兩種，基于中心差分算法的時域方法的FLAC3D動力分析要求數值模型中單元最大邊長不超過待輸入波動信號在其中傳播最小波長的1/10。因此，分析中常對波動信號最大頻率作截斷處理，以期求得計算效率與求解精度兩者間的合理平衡；
基線修正：波動信號頻繁特征復雜，其中的長周期成分可能會導致地基模型整體出現平動漂移現象，這部分對工程安全無影響的位移應予以剔除，此即為基線修正的基本原理。
FLAC3D V7.0含向導式插件工具Dynamic Input Wizard來滿足對波動信號開展上述處理的要求，該工具通過程序菜單選項Tools->Dynamic Input Wizard實現調用，起始界面如圖3所示。考慮篇幅原因，此處不對其使用方法作進一步介紹。
圖3 波動信號處理工具Dynamic Input Wizard
下載地址：FLAC 3D實用教程
展開 Msmin取值區間一般為[0.05，0.1]。
圖2 土體動剪切模量退化行為截斷處理
在V7.0版本中，FLAC3D通過引入如下命令及選項來反映對土體動剪切模量退化行為的截斷處理：
zone dynamic damping hysteretic keyword
keyword：
reduction-minimum：模量退化系數Ms的截斷值。
波動信號處理
動力分析中因分析方法的不同通常對待輸入波動信號的成分組成有規定要求，因此在分析之初有必要依據頻譜分析方法對信號做預處理，處理方法以濾波和基線修正為主：
濾波：目前動力分析方法主要包括頻域法和時域法兩種，基于中心差分算法的時域方法的FLAC3D動力分析要求數值模型中單元最大邊長不超過待輸入波動信號在其中傳播最小波長的1/10。因此，分析中常對波動信號最大頻率作截斷處理，以期求得計算效率與求解精度兩者間的合理平衡；
基線修正：波動信號頻繁特征復雜，其中的長周期成分可能會導致地基模型整體出現平動漂移現象，這部分對工程安全無影響的位移應予以剔除，此即為基線修正的基本原理。
FLAC3D V7.0含向導式插件工具Dynamic Input Wizard來滿足對波動信號開展上述處理的要求，該工具通過程序菜單選項Tools->Dynamic Input Wizard實現調用，起始界面如圖3所示。考慮篇幅原因，此處不對其使用方法作進一步介紹。
圖3 波動信號處理工具Dynamic Input Wizard
下載地址：講一下Flac3D的局部坐標
展開 為了與FLAC的本構模型比較，可參考下述鏈接：
IMASS---FLAC3D和3DEC新的本構模型(2)
FLAC2D---過去，現在和將來
FLAC3D 7.0 新特性簡介(P3)---新的本構模型
(2) 用戶自定義的塑性、流動規則和硬化模型。
(3) 混凝土模型: Bresler-Pister; Ottosen; William–Warnke
Bresler and Pister對同時受到扭轉和軸向壓縮的混凝土圓柱體進行了試驗，從這些實驗結果中得出的斷裂準則如下所示：
Ottosen模型如下所示:
(4) 脆性損傷模型(Brittle damage models)
(5) Hoek-Brown模型
(6) 蠕變模型 (Creep models)
(7) 預先定義和用戶定義的非線性彈性材料: Ramberg-Osgood, Hardin-Drnevich, Duncan-Chang等。
這些材料模型也可以耦合創建的任何新方程，并與其它物理場, 如熱傳導、流體流動和多孔介質中的溶質傳輸等進行耦合。
3 Comsol在巖土力學中的應用
Nasir(2014) 使用耦合的thermo-hydro-mechanical-chemical (THMC) 模型研究了氣候變化影響下核廢料場地沉積巖的孔隙度和滲透率變化；Zhou(2014) 模擬了水庫水和降雨對三峽邊坡滑動體的滲透；Zhang (2015)模擬了隧道圍巖塑性區的發展; Cao (2016)模擬了采礦引起的斷裂導致瓦斯遷移；Wang (2018) 模擬了各向異性巖體的邊坡損傷帶；Beya (2019) 模擬了在永久性凍土采礦中水泥漿充填加熱固化條件的過程；Liu (2020) 模擬了斷裂花崗巖中地熱的傳播。
展開 fish set @high_unbal = 5e4fish set @low_unbal = 2e4fish set @high_vel = 2
伺服控制函數在UDEC和FLAC中寫起來比較繁瑣，主要原因是這些2D程序還沒有完全轉換到新的FISH格式(FISH: Loop語句的進化; FLAC3D 7.0 新特性簡介(P2)---FISH的顯著改進); FLAC3D和3DEC的伺服控制函數寫起來很簡單，示例如下:
fish def _servo while_stepping if block.unbal > unbal_limit then vel_ = 0.98*vel_ loop foreach local gp points block.gp.vel.app.z(gp) = vel_ end_loop endif if block.unbal < 0.8*unbal_limit then vel_ = 1.02*vel_ loop foreach gp points block.gp.vel.app.z(gp) = vel_ end_loop endifend
3 應變軟化材料的三軸壓縮試驗
下圖所示的是一個應變軟化材料的三軸壓縮試驗結果。其中左圖顯示的是計算的平均垂直應力sigmav和平均垂直應變ev之間的關系，右圖顯示的是最大不平衡力與時步之間的關系。比較引言中的單軸應力應變曲線，可以發現在三軸條件下應變軟化的行為更加復雜。
展開 兩種計算采用了不同的動態邊界條件：
在自由振動(free vibration)中，使用粘性邊界(Viscous boundary)條件【PLAXIS 機器地基動力分析---Part II】，這種設置適用于動態源的位置位于網格內的問題【風速(Wind Velocity)計算】；對于地震荷載，使用自由場(Free-field)【發布PLAXIS 2D V22.02.00 新的改進要點；FLAC3D 7.0 新特性簡介(P1)---速度提升】和符合基礎邊界(Compliant base boundary)條件，此選項是地震分析的首選，在模型底部施加動態輸入，符合基礎邊界條件用來吸收向下的波，從而使波只向上傳播。
本題的主要內容包括：
(1) 進行動態計算
(2) 定義動態邊界條件(自由場、符合基礎和粘性）
(3) 使用動態放大系數(dynamic multipliers)定義地震
(4) 模擬結構的自由振動
(5) 使用具有小應變剛度的硬化土(Hardening Soil)模型對滯回行為(hysteretic behaviour)進行模擬
(6) 評估傅里葉頻譜(Fourier spectrum)的固有頻率
我們側重地震載荷的動力分析。
2 物理模型
該建筑物由5層樓板和一
個地下室組成。
建筑物寬10米，每層高3m，因此5層高共15m, 地下室高2m，樓板和墻體的重量按
5kN/m^2估算。
建筑物的地基置于15m厚的松散砂層上，其下的持力層是致密的砂層，厚度為25m，如下圖所示。
展開 
Flac3D7.0的相關專題、標簽、搜索
Flac3D7.0的最新內容
兩種計算采用了不同的動態邊界條件：
在自由振動(free vibration)中，使用粘性邊界(Viscous boundary)條件【PLAXIS 機器地基動力分析---Part II】，這種設置適用于動態源的位置位于網格內的問題【風速(Wind Velocity)計算】；對于地震荷載，使用自由場(Free-field)【發布PLAXIS 2D V22.02.00 新的改進要點；FLAC3D 7.0
或者說，以往FF邊界原則上不適用于利用其開展2.5維FLAC3D模型的動力響應分析；FLAC3D V7.0則對該環節予以了改進。
或者說，以往FF邊界原則上不適用于利用其開展2.5維FLAC3D模型的動力響應分析；FLAC3D V7.0則對該環節予以了改進。
fish set @high_unbal = 5e4fish set @low_unbal = 2e4fish set @high_vel = 2
伺服控制函數在UDEC和FLAC中寫起來比較繁瑣，主要原因是這些2D程序還沒有完全轉換到新的FISH格式(FISH: Loop語句的進化; FLAC3D 7.0 新特性簡介(P2)---FISH的顯著改進); FLAC3D和3DEC的伺服控制函數寫起來很簡單
《本構模型(Constitutive Models)選擇》《IMASS---FLAC3D和3DEC新的本構模型(1)》《FLAC3D 7.0 新特性簡介(P3)---新的本構模型》《FLAC2D---過去，現在和將來》。迄今為止，FLAC3D/3DEC已經內置了超過35種本構模型，當模擬一個問題時，我們不可能試驗所有的本構模型。
本構模型
FLAC3D/3DEC包含了線性和非線性的本構模型[FLAC3D 7.0 新特性簡介(P3)---新的本構模型]，通過下面的代碼調用。
define find_dip sum = 0 loop foreach fp fracture.list fdip = fracture.dip(fp) io.out(fdip) sum = sum + fdipend_loop
新的循環算法只使用一句代碼，目前在3DEC，FLAC3D和PFC 7.0 中能夠使用，但還沒有改進到UDEC
define find_dip sum = 0 loop foreach fp fracture.list fdip = fracture.dip(fp) io.out(fdip) sum = sum + fdipend_loop
新的循環算法只使用一句代碼，目前在3DEC，FLAC3D和PFC 7.0 中能夠使用，但還沒有改進到UDEC
現在Itasca的3DEC、FLAC3D和PFC從7.0版本開始為網絡許可證用戶免費推出Linux版本，運行在Ubuntu 20.04 Linux操作系統。Ubuntu是一個非常流行的開源Linux操作系統，可以在Windows系統下安裝，實現雙系統運行；也可以使用可啟動的U盤，不需要把系統安裝在硬盤上。
2 L版比W版運行速度快？
1 引言
這個筆記是對大約一年半前《FLAC3D與Python的集成》系列的更新，當時是在FLAC3D 6.0環境下進行的，盡管基本的操作流程和6.0版本相同，但目前的FLAC3D 7.0發生了很大變化，不僅增強了Python的運行環境(Jupyter QtConsole 4.3.1,Python 3.6.1, IPython 6.2.1,matplotlib backend: Qt5Agg
