Flac3D本構模型
關注創建者：Flac3D程小杰 創建時間：2022-10-10
Flac3D本構模型的實例教程
<p>基于fish語言的flac3d的interface界面單元非線性本構模型開發實例</p>
</p><p>本文將介紹一種使用FLAC3D模擬頂板巖層預裂卸壓效果，主要包含：（1）數值模型建立、預制斷裂接觸帶和定義爆破軟化區；（2）原巖應力場平衡；（3）巷道循環開挖及錨桿支護；（4）采用雙屈服本構模型模擬采空區；（6）對目標巖層實施斷裂</p><p><br></p><figure style="text-align: center;" class="ql-align-center">
<figure class="figure-image" contenteditable="false" data-img="https://img.jishulink.com/202507/attachment/102ed06e8bae4b72a59703055ce3d9db.png" style="display: inline-block;">
<img src="https://img.jishulink.com/202507/attachment/102ed06e8bae4b72a59703055ce3d9db.png" data-mobile-src="https://img.jishulink.com/202507/attachment/102ed06e8bae4b72a59703055ce3d9db.png?image_process=/format,webp" data-pc-src="https://img.jishulink.com/202507/attachment/102ed06e8bae4b72a59703055ce3d9db.png?
展開 <p>本案例提供一個VUMAT子程序，用于近似復刻Abaqus自帶3維摩爾庫倫本構退化Tresca準則。通過一個椎體勻速貫入案例予以驗證VUMAT子程序有效性和效率。方便用戶在此子程序基礎上進行更為復雜的土體強度參數調整，構建更為復雜的本構模型。附三個時刻Mises應力場對比結果。</p><p><br></p><figure style="text-align: center;" class="ql-align-center">
<figure class="figure-image" contenteditable="false" data-img="https://img.jishulink.com/202509/attachment/d3e5056748bc4317b9f2d39379fe90a3.png" style="display: inline-block;" data-regular="true">
<img src="https://img.jishulink.com/202509/attachment/d3e5056748bc4317b9f2d39379fe90a3.png" style="" width="379" data-mobile-src="https://img.jishulink.com/202509/attachment/d3e5056748bc4317b9f2d39379fe90a3.png?image_process=/format,webp" data-pc-src="https://img.jishulink.com/202509/attachment/d3e5056748bc4317b9f2d39379fe90a3.png?
展開 對于纖維增強復合材料的模擬，在<a href="/major/ABAQUS中，集成了二維Hashin失效準則與多種損傷演化準則，但缺少三維的復合材料本構模型。
參考一篇已發表的SCI文章，使用Fortran語言建立三維平紋織物復合材料彈塑性、漸進損傷本構模型-Vumat子程序。平紋織物復合材料在1方向和2方向絲束性能近似相同。
該程序是博士期間學習復材子程序的小部分總結，編程結構并不是非常漂亮及完美，但確保能順利運行，且單元驗證結果與理論公式一致，介意請勿拍。
程序中塑性迭代部分并非主流的牛頓-拉夫遜和梯度下降方法，但經過驗證能夠適用于該模型，介意請勿拍。
附件內容：1. inp算例模型(低速沖擊工況，1/4模型，層間使用cohesive element) 2. 子程序 3 .使用方法 4.參考論文名稱
首先介紹該子程序的使用方法與效果
1. 在ABAQUS中建立三維復合材料模型，這里建立一個簡單的方塊。賦給材料方向，1,2方向分別表示絲束的方向，3方向表示垂直于1，2的方向，也就是面外方向。
2. 建立材料屬性
3. 建立顯示Explicit計算時間步，時間0.005，在場輸出中勾選輸出 SDV和 STATUS.
4. 劃分網格，賦給Explicit 3D stress單元類型，邊界條件根據需要設定即可，此處不再贅述。此處劃分為一個單元，使用12方向往復加載卸載。建立Job，提交模型前在Job中選擇該子程序，設置雙精度計算。
5. 查看結果，等效塑性應變在卸載時沒有變化，再次加載時剪切應力按照原來的路徑返回，剪切損傷在卸載時也保持不變。
6. 將該子程序應用在低速沖擊模型中，可以順利運行。
接下來簡要介紹該子程序的相關理論，子程序、參考的論文名稱以及輸入材料參數的對應含義打包在附件中。
展開 工程實踐中最流行最可靠的分析程序是FLAC，它可以進行1D,2D和3D場地反應分析，其它一些基于有限元的數值程序，例如Plaxis, GTS NX, DIANA, Quake/W等也能進行簡單的二維分析。參考：
Dynamic Modeling Considerations
Comparison of FLAC3D to SHAKE for a Layered, Linear-Elastic Soil Deposit
Comparison of FLAC3D to SHAKE for a Layered, Nonlinear-Elastic Soil Deposit
非線性模型可以通過使用先進的FLAC3D本構模型，CYSoil, CHSoil, Soft-Soil, NorSand, P2PSand, Finn等更準確地捕捉經歷高應變場地，例如軟土場地和液化場地的反應。而當土中的剪切應變較低時，例如地面運動較弱或場地土較硬，等效的線性模型也能得到相似的結果，從實踐的角度來講，在這種情形下不必使用非線性模型。根據問題的復雜性和重要性，線性模型和非線性模型都可以使用更高的維度2D或3D。
3 工程場地地震安全性評價
對于處在地震活動區域的工程場地，必須進行地震安全性評價。國內《工程場地地震安全性評價》(GB 17744-2005)推薦使用的是一維等效線性方法，場地典型土剪切模量和阻尼比的取值如下表所示。
根據鉆孔測得的剪切波速和密度，可以得到不同樣本和不同超越概率水平下的峰值加速度以及地震動相關反應譜。有一點需要說明的是進行場地反應分析時，國內建議的鉆孔深度好像是20m(還沒找到正式的出處)，在北美建議的鉆孔深度通常是30m。
展開 
Flac3D本構模型的相關專題、標簽、搜索
Flac3D本構模型的最新內容
<p>本案例提供一個VUMAT子程序，用于近似復刻Abaqus自帶3維摩爾庫倫本構退化Tresca準則。通過一個椎體勻速貫入案例予以驗證VUMAT子程序有效性和效率。方便用戶在此子程序基礎上進行更為復雜的土體強度參數調整，構建更為復雜的本構模型。附三個時刻Mises應力場對比結果。</p><p><br></p><figure style="text-align: center;" class="ql-align-center
<p> 堅硬頂板由于強度高、整體性好、難以垮落，在回采過程中形成大跨度懸露結構，導致應力積聚。當頂板發生突然垮斷時，會瞬間釋放大量能量，對巷道和工作面造成劇烈沖擊，誘發沖擊地壓等災害。因此，預裂卸壓技術（如高位定向長孔預裂爆破）被廣泛應用于煤礦采場應力調控實踐中。</p><p>本文將介紹一種使用FLAC3D模擬頂板巖層預裂卸壓效果，主要包含：（1）數值模型建立
<p>基于fish語言的flac3d的interface界面單元非線性本構模型開發實例</p>
對于纖維增強復合材料的模擬，在<a href="/major/ABAQUS中，集成了二維Hashin失效準則與多種損傷演化準則，但缺少三維的復合材料本構模型。
參考一篇已發表的SCI文章，使用Fortran語言建立三維平紋織物復合材料彈塑性、漸進損傷本構模型-Vumat子程序。平紋織物復合材料在1方向和2方向絲束性能近似相同。
該程序是博士期間學習復材子程序的小部分總結，編程結構并不是非常漂亮及完美
本構模型，CYSoil, CHSoil, Soft-Soil, NorSand, P2PSand, Finn等更準確地捕捉經歷高應變場地，例如軟土場地和液化場地的反應。
