節(jié)理巖體模擬的案例
Ubiquitous節(jié)理巖體模擬(Ubiquitous Joint Rock Mass Modelling)
1 引言
Ubiquitous節(jié)理巖體模擬最初在FLAC中引入，Sainsbury等人(2008)在SRM中使用了UJRM，從而在FLAC3D模型中考慮了巖體強(qiáng)度和各向異性, 能夠研究完整巖石內(nèi)聚力的逐漸減弱和Ubiquitous節(jié)理破壞, [Sainsbury, B. et al. (2008) Analysis of cave behaviour using a Synthetic Rock Mass (SRM) – Ubiquitous Joint Rock Mass (UJRM) modelling technique.]，他們的方法與Pierce等人(2007)和Mas Ivars等人(2008)的方法不同，沒有使用PFC3D，而是用FLAC3D評(píng)估來自較小模型模擬的巖體屬性。Lavoie(2011) 提出了一個(gè)UJRM的解析解[Thierry Lavoie (2011) An Analytical Geomechanical Upscaling Approach for Modeling Jointed Rock Mass Behaviour Using Ubiquitous Joints. UBC MS. Thesis. 186p.]。這個(gè)筆記簡(jiǎn)要討論了ADONIS的Ubiquitous本構(gòu)模型，并與FLAC/Slope中的Ubiquitous模型作了比較。
2 ADONIS模擬
一個(gè)簡(jiǎn)單的Ubiquitous節(jié)理巖體邊坡如下圖所示。
(1) 項(xiàng)目設(shè)置
在項(xiàng)目設(shè)置中，主要設(shè)置合適的單位和重力加速度。
(2) 網(wǎng)格劃分
邊坡的幾何形狀可以在ADONIS內(nèi)直接輸入坐標(biāo)，選擇6節(jié)點(diǎn)的三角形單元，最大邊長取0.75。
(3) 邊界條件
邊界使用xyfix命令固定三邊，邊坡面不能固定。
展開 擴(kuò)展有限元(XFEM)模擬巖體節(jié)理網(wǎng)絡(luò)(DFN | Joint Network）
1 引言
節(jié)理巖體邊坡穩(wěn)定性分析是巖石工程數(shù)值模擬最具挑戰(zhàn)性的問題之一，因?yàn)?em>節(jié)理的分布模式及其強(qiáng)度對(duì)邊坡穩(wěn)定性有重大影響，而這兩個(gè)因素在很大程度上都是不確定的。當(dāng)進(jìn)行數(shù)值模擬分析時(shí)，有兩種不同途徑處理節(jié)理化巖體的應(yīng)力和應(yīng)變行為，最常用的方法是離散元法，離散元法把節(jié)理作為邊界來處理，由節(jié)理切割而成的巖塊可以是剛性體，也可以是變形體，這種方法最廣泛應(yīng)用的程序是UDEC和3DEC，本公眾號(hào)有近200篇文章討論了這種方法；另一種方法是把巖體作為連續(xù)性介質(zhì)，節(jié)理當(dāng)作另一種材料模型，但節(jié)理置于有限元網(wǎng)格的節(jié)點(diǎn)上，當(dāng)巖體發(fā)生變形時(shí)，節(jié)理隨著單元網(wǎng)格移動(dòng)，但不能象離散元那樣節(jié)理產(chǎn)生分離【Voronoi多邊形在有限元中的應(yīng)用】，如下左圖所示。不過，這種處理方法在復(fù)雜的幾何結(jié)構(gòu)中會(huì)表現(xiàn)出顯著的缺點(diǎn)，因?yàn)樵谶@種情形下節(jié)理不能完全適應(yīng)網(wǎng)格的劃分。為了克服這一限制，引入了隱式的節(jié)理網(wǎng)絡(luò)巖體模擬方法---擴(kuò)展有限元方法(XFEM), XFEM能夠單獨(dú)識(shí)別節(jié)理，與網(wǎng)格劃分無關(guān)。節(jié)理不必置于有限元網(wǎng)格的節(jié)點(diǎn)上，可以自由分布在所研究的域內(nèi)，如下右圖所示。
2 XFEM工作機(jī)理
擴(kuò)展有限元法(XFEM)【Abaqus 2021 擴(kuò)展有限元 XFEM新功能；[最新文獻(xiàn)]錨桿加固全飽和斷裂多孔介質(zhì)的數(shù)值過程】是21世紀(jì)初開發(fā)的一種新的數(shù)值方法[Belytschko T., Black T., 1999, Elastic crack growth in finite elements with minimal remeshing, International Journal for Numerical methods in Engineering, 45, 601-620.]
展開 利用3DEC仿真三維規(guī)則節(jié)理巖體的文章
個(gè)別元素法于三維巖體力學(xué)行為之應(yīng)用
楊長義 陳志民 陳錦清 淡江大學(xué)土木工程學(xué)研究所 中興顧問社大地力學(xué)研究中心
摘 要本文利用三維個(gè)別元素分析法程序(3DEC)，仿真三維規(guī)則節(jié)理巖體在單軸與真三軸應(yīng)力下之變形與強(qiáng)度特性。主要研究結(jié)果如下：(1) 3DEC可用以定性分析三維節(jié)理巖體力學(xué)機(jī)制，利用該程序可簡(jiǎn)易探討任何節(jié)理分布位態(tài)下之力學(xué)行為，免除物理模模型仿真試驗(yàn)之困難；(2) 多軸應(yīng)力下巖體之異向性行為亦可透過3DEC仿真分析，定性上均與物理現(xiàn)象相互一致；(3)在節(jié)理間距、勁度比較大的巖盤較需要比對(duì)二維與三維分析之差異。
一、前 言
自然界巖體多處于真三軸應(yīng)力狀態(tài)下，以往受限于分析工具與實(shí)驗(yàn)設(shè)備，巖石工程之分析大多局限于二向度分析，對(duì)三維巖體行為之仿真則較少[1]。例如目前可用于分析具大變形特性之離散巖體的程序如DDA[2]、UDEC[3]均局限于對(duì)二維問題的解析；而3DEC[4]程序系以個(gè)別元素法(distinct element method)在UDEC基礎(chǔ)下發(fā)展而成之?dāng)?shù)值分析程序，正可用以仿真三維節(jié)理巖體之力學(xué)行為：3DEC可將巖體視為由許多離散的完整巖塊所組成，各完整巖塊可以仿真成剛體或可變形體；而各完整巖塊間為節(jié)理所分隔。(1)在對(duì)節(jié)理的仿真方面，主要依據(jù)位移－作用力法則，計(jì)算在節(jié)理面上之剪應(yīng)力及正向應(yīng)力，以作為個(gè)別巖塊之邊界應(yīng)力條件，因此可仿真巖塊大位移與轉(zhuǎn)動(dòng)之情況。(2) 3DEC在仿真可變形巖塊時(shí)，系根據(jù)「edge」指令程序自行將三維巖塊再細(xì)分成許多四面體狀次級(jí)塊體(sub-block)，次級(jí)塊體可以為任意形狀。
展開 舊文新讀: 節(jié)理化巖石的合成巖體模型[A SRM Model for Jointed Rock]
在過去, 預(yù)測(cè)節(jié)理化巖體強(qiáng)度和變形最廣為接受的方法是基于半經(jīng)驗(yàn)估計(jì)的Hoek-Brown破壞準(zhǔn)則, 然而Hoek-Brown準(zhǔn)則沒有強(qiáng)調(diào)巖體的脆性行為, 這也部分地造成Hoek-Brown準(zhǔn)則的估算有時(shí)會(huì)引起人們的懷疑, 例如<Hoek的巖體變形模量經(jīng)驗(yàn)估計(jì)---Is it reliable ?>. 而在塊體崩落法和分段崩落法中, 考慮巖石從峰值強(qiáng)度到巖石的完全破碎過程是必須的.
SRM使用PFC3D建立的BPM模型表示原巖, 用離散斷裂網(wǎng)絡(luò)DFN表示節(jié)理, 如下圖所示, 來估算巖體的強(qiáng)度屬性和變形屬性, 從而預(yù)測(cè)巖體的脆性斷裂行為.
建模方法是把斷裂信息(即斷裂幾何形狀和特性)疊加到BPM模型上。即BPM對(duì)完整巖石建模，通過修改斷裂接觸點(diǎn)處的接觸模型引入斷裂的力學(xué)行為。由于PFC模型本質(zhì)上是離散的，因此破壞可能在完整的BPM區(qū)域和沿?cái)嗔衙姘l(fā)生。
SRM克服了早期工作中存在的模型尺寸和表示節(jié)理的限制，允許快速構(gòu)造和測(cè)試直徑為10-100米的中等到嚴(yán)重節(jié)理化的巖石樣本，這些巖石含有數(shù)千個(gè)非貫通性節(jié)理。SRM模擬用于估算巖體的峰前屬性(模量、損傷閾值、峰值強(qiáng)度)和峰后屬性(脆性、殘余強(qiáng)度、破碎)，并用于分析大規(guī)模邊值問題(例如巖體邊坡穩(wěn)定性)。SRM模擬允許考慮三維大型復(fù)雜非貫通性節(jié)理網(wǎng)絡(luò)以及塊體斷裂，包括不完整節(jié)理對(duì)塊體強(qiáng)度的影響。
一個(gè)通用的SRM數(shù)值模擬步驟如下:(1) 根據(jù)現(xiàn)場(chǎng)數(shù)據(jù)建立離散斷裂網(wǎng)絡(luò)DFN; (2) 對(duì)構(gòu)造的DFN模型進(jìn)行了隨機(jī)抽樣，按恒定的高寬比分離出N個(gè)立方樣本進(jìn)行模擬; (3) 對(duì)每個(gè)立方樣本進(jìn)行強(qiáng)度試驗(yàn), 并記錄了每個(gè)樣本的全部應(yīng)力應(yīng)變行為。這種方法為確定大規(guī)模巖體樣本的復(fù)雜構(gòu)成行為提供了一種方法。這在實(shí)驗(yàn)室中通常很難實(shí)現(xiàn)或不可能實(shí)現(xiàn)。
展開 
利用3DEC仿真三維規(guī)則節(jié)理巖體的文章
摘 要本文利用三維個(gè)別元素分析法程序(3DEC)，仿真三維規(guī)則節(jié)理巖體在單軸與真三軸應(yīng)力下之變形與強(qiáng)度特性。主要研究結(jié)果如下：(1) 3DEC可用以定性分析三維節(jié)理巖體力學(xué)機(jī)制，利用該程序可簡(jiǎn)易探討任何節(jié)理分布位態(tài)下之力學(xué)行為，免除物理模模型仿真試驗(yàn)之困難；(2) 多軸應(yīng)力下巖體之異向性行為亦可透過3DEC仿真分析，定性上均與物理現(xiàn)象相互一致；(3)在節(jié)理間距、勁度比較大的巖盤較需要比對(duì)二維與三維分析之差異。
一、前 言
自然界巖體多處于真三軸應(yīng)力狀態(tài)下，以往受限于分析工具與實(shí)驗(yàn)設(shè)備，巖石工程之分析大多局限于二向度分析，對(duì)三維巖體行為之仿真則較少[1]。例如目前可用于分析具大變形特性之離散巖體的程序如DDA[2]、UDEC[3]均局限于對(duì)二維問題的解析；而3DEC[4]程序系以個(gè)別元素法(distinct element method)在UDEC基礎(chǔ)下發(fā)展而成之?dāng)?shù)值分析程序，正可用以仿真三維節(jié)理巖體之力學(xué)行為：3DEC可將巖體視為由許多離散的完整巖塊所組成，各完整巖塊可以仿真成剛體或可變形體；而各完整巖塊間為節(jié)理所分隔。(1)在對(duì)節(jié)理的仿真方面，主要依據(jù)位移－作用力法則，計(jì)算在節(jié)理面上之剪應(yīng)力及正向應(yīng)力，以作為個(gè)別巖塊之邊界應(yīng)力條件，因此可仿真巖塊大位移與轉(zhuǎn)動(dòng)之情況。(2) 3DEC在仿真可變形巖塊時(shí)，系根據(jù)「edge」指令程序自行將三維巖塊再細(xì)分成許多四面體狀次級(jí)塊體(sub-block)，次級(jí)塊體可以為任意形狀。每個(gè)次級(jí)塊體可配合所指定之材料組成律及外力情況，計(jì)算巖塊之受力及應(yīng)力分布情況；每個(gè)次級(jí)塊體的節(jié)點(diǎn)有三個(gè)自由度，以計(jì)算這些次級(jí)塊體上節(jié)點(diǎn)之運(yùn)動(dòng)情形，然后配合材料組成律計(jì)算次級(jí)塊體上之應(yīng)力應(yīng)變關(guān)系，則可得塊體間之作用力，接著配合邊界所產(chǎn)生的接觸力計(jì)算得新合力與加速度，以作為下一時(shí)階計(jì)算可變形巖塊之邊界條件。
個(gè)別元素法于三維巖體力學(xué)行為之應(yīng)用.doc
展開 應(yīng)變軟化模型(Strain-Softening and IMASS)
block zone cmodel assign strain-softeningblock zone property dens 2500 bulk 1.19e10 shear 1.1e10block zone property cohesion 2.72e5 friction 44 tension 2e5block zone property table-cohesion 'coh' table-friction 'fri'table 'coh' add (0,2.72e5) (1e-4,2e5) (2e-4,1.5e5) (3e-4,1.03e5) (1,1.03e5)table 'fri' add (0,44) (1e-4,42) (2e-4,40) (3e-4,38) (1,38)
table的使用方法參考下面的鏈接：
輸入和輸出table的完整說明
輸入和輸出table: table export/import
增加和改進(jìn)table: table add/insert | table.value
與外部數(shù)據(jù)交互: table命令與table FISH
與應(yīng)變軟化Mohr-Coulomb模型類似，雙線性的應(yīng)變軟化(Bilinear Strain-Softening) Ubiquitous-Joint模型[Ubiquitous節(jié)理巖體模擬(Ubiquitous Joint Rock Mass Modelling)]能夠模擬巖體中分布的軟弱面。
展開 二維有限元分析軟件ADONIS新功能評(píng)述
關(guān)于這個(gè)軟件的部分功能可參看以下過去的文章：
免費(fèi)二維有限元分析軟件ADONIS計(jì)算步驟
邊坡剪切強(qiáng)度折減分析(Shear Strength Reduction Analysis)
板樁墻(Sheet Pile Wall)模擬---FEM中的界面元
Ubiquitous節(jié)理巖體模擬(Ubiquitous Joint Rock Mass Modelling)
2 Swell Model
本軟件的Swell模型模仿了FLAC3D的Swell Model，因此下面的討論基于FLAC3D使用的約定。Swell Model的基礎(chǔ)是Mohr-Coulomb模型的非關(guān)聯(lián)剪切和關(guān)聯(lián)拉伸流動(dòng)規(guī)則，與之不同的是，通過耦合潤濕應(yīng)變(wetting strains)與潤濕前的模型狀態(tài)，考慮了潤濕引起的變形(wetting-induced deformations)，但屈服和勢(shì)能函數(shù)、塑性流動(dòng)規(guī)則以及應(yīng)力修正與Mohr-Coulomb模型的假設(shè)相同。FLAC3D的Swell Model是專為由濕潤引起的土的變形發(fā)展的，與Plaxis和RS中的Swelling Rock Model有差別，這個(gè)模型的應(yīng)用參看如下論文：
[1] Noorany, I. et al. (1999) Prediction of soil slope deformation due to wetting.
[2] Pathak, Y. et al. (2003) Wetting-induced deformation of geosynthetic reinforced slopes with expansive soils.
展開 COMSOL Multiphysics?模擬裂隙對(duì)巖體滲流的影響 ￥30
研究了理想裂隙引起的流動(dòng)擾動(dòng)及其幾何排列對(duì)裂隙巖體穩(wěn)態(tài)滲流（有效滲透率、壓力分布、流線等）的影響；
參考論文：Fluid flow partitioning between fractures and a permeable rock matrix,
https://doi.org/10.1029/2003GL019027
論文原圖：
COMSOL Multiphysics?計(jì)算云圖（由于參數(shù)不同，云圖存在些許差別）：
為了方便對(duì)比，不同裂隙幾何、參數(shù)影響結(jié)果都放在同一個(gè)模型不同研究中了。
技能包含：模型計(jì)算設(shè)置方式；計(jì)算云圖后處理方式（參考論文）；等效滲透率計(jì)算方式
展開 基于DYNA的含裂隙巖體爆破裂紋及擴(kuò)展及損傷模擬 ￥48.9
巖體中的裂隙直接影響爆破應(yīng)力波在巖體中的傳播，進(jìn)而影響爆破效果。本案例研究了含裂隙巖體爆破中裂紋的擴(kuò)展及損傷過程，Ls-dyna模擬了爆炸應(yīng)力波在裂隙巖體中的傳播特性。結(jié)果表明:爆破過程中應(yīng)力波分布不均勻，主要向自由表面?zhèn)鞑ラ]合型宏觀裂隙阻礙爆炸應(yīng)力波的傳播，且在裂隙處會(huì)止裂，裂紋及損傷會(huì)繞過裂隙處，模擬結(jié)果如下 ：
圖1 含裂隙巖體爆破裂紋及擴(kuò)展有限元模型
圖2 含裂隙巖體爆破裂紋及擴(kuò)展過程
圖3 含裂隙巖體爆破裂紋及等效應(yīng)力波傳播過程
本案例適用于研究爆炸、沖擊、侵徹動(dòng)力學(xué)的朋友，下面附上該模擬的K文件，大家有疑問可以在私信我，歡迎交流！
展開 基于LS-DYNA的準(zhǔn)二維巖體爆破裂紋的模擬(附K文件） ￥19.8
巖體中的炸藥在爆炸瞬間釋放出大量的爆炸能量，產(chǎn)生爆炸沖擊波和應(yīng)力波，以動(dòng)載的形式作用周圍巖體并使得巖石損傷和破碎。一般認(rèn)為，巖體的爆破裂紋是由兩種不同類型的動(dòng)力荷載作用下產(chǎn)生的:爆炸應(yīng)力波和爆炸氣體。首先是炸藥起爆后在孔洞周圍產(chǎn)生的應(yīng)力沖擊波的作用，其次是持續(xù)時(shí)間較長的爆炸性氣體的作用。當(dāng)炸藥起爆時(shí)，壓力急劇增大，爆破孔周圍巖體被壓碎形成破碎區(qū)，而周圍巖體受到拉應(yīng)力產(chǎn)生裂縫，形成裂隙擴(kuò)展區(qū)。本案例巖石模型采用003號(hào)材料*MAT_PLASTIC_KINEMATIC，裝藥方式采用空氣不耦合裝藥，通過定義失效準(zhǔn)則，使得巖體產(chǎn)生破碎及擴(kuò)展裂隙，模擬結(jié)果如下
本案例適用于研究爆炸、沖擊、侵徹動(dòng)力學(xué)的朋友，下面附上該模擬的K文件，大家有疑問可以在私信我，歡迎交流！
展開 ABAQUS中非均質(zhì)煤巖體的水力壓裂模擬 ￥80
今天給大家?guī)硪粋€(gè)有關(guān)水力壓裂模擬的案例，廢話少說直接上論文：
非均質(zhì)彈塑性煤體水壓致裂裂紋形態(tài)研究 - 中國知網(wǎng) (cnki.net)
模擬裂紋擴(kuò)展的方法有很多，但我覺得FDEM方法是模擬縫網(wǎng)形態(tài)的最好方法，這篇論文使用ABAQUS通過二次開發(fā)實(shí)現(xiàn)FDEM方法，相關(guān)方法星辰北極星大佬早就使用過了，這里我們主要解決兩個(gè)問題：
1、材料非均質(zhì)性的實(shí)現(xiàn)方法
2、多射孔多注入點(diǎn)的集中流量注入實(shí)現(xiàn)
模擬的基本模型很簡(jiǎn)單，辦公室的工作站只夠帶動(dòng)的二維的；因?yàn)榭紤]了地層的非均質(zhì)性，所以計(jì)算不易收斂，這個(gè)二維模型為100w自由度，40核工作站算一步10分鐘左右，單步增量時(shí)長1e-5到-7的樣子，所以還請(qǐng)大家斟酌后建模。具體模型如下圖：15m*15m的地層中有一簇射孔，詳細(xì)尺寸圖中已經(jīng)給出啦。
這一簇射孔在際水力壓裂中是一個(gè)分隔段，我們需要將中間四個(gè)注液點(diǎn)耦合到一個(gè)節(jié)點(diǎn)上，給這個(gè)節(jié)點(diǎn)設(shè)置集中注液流量，被耦合的四個(gè)點(diǎn)根據(jù)孔隙壓力自己平衡每個(gè)點(diǎn)的注液流量。
展開 
基于LS-DYNA的巖體微差爆破動(dòng)態(tài)響應(yīng)模擬 ￥38.79
基于流固耦合算法在動(dòng)態(tài)分析軟件LS-DYNA中研究爆炸沖擊波的傳播規(guī)律在巖石介質(zhì)和周圍巖體的爆破振動(dòng)影響下的爆破孔同時(shí)起爆和微差爆破兩個(gè)工況時(shí)。結(jié)果表明:雙孔同時(shí)起爆初期，損傷破碎區(qū)擴(kuò)展與單孔爆破相似，爆炸沖擊波,彼此是重疊的,兩炮孔中間縱向單元和藥柱內(nèi)外兩側(cè)橫向近區(qū)單元的壓力和等效應(yīng)力隨爆心距的增大而減小，而自由面上單元呈現(xiàn)出先增后減的變化趨勢(shì)， 微差起爆可緩解爆破振動(dòng)和改善爆破效果，模擬結(jié)果對(duì)比如下：
圖1 雙孔同時(shí)起爆時(shí)等效應(yīng)力變化過程
圖2 雙孔延期起爆時(shí)等效應(yīng)力變化過程
本案例適用于研究爆炸、沖擊、侵徹動(dòng)力學(xué)的朋友，下面附上該模擬的K文件，大家有疑問可以在私信我，歡迎交流！
展開 COMSOL基于漿液黏度時(shí)空變化的水平裂隙巖體注漿擴(kuò)散數(shù)值模擬 ￥210
速凝類漿液的雙液混合注漿方式及其黏度時(shí)變特性導(dǎo)致漿液擴(kuò)散區(qū)內(nèi)黏度空間分布不均勻。基于此，認(rèn)為速凝類漿液流型為具有黏度時(shí)變性的賓漢流體，研究其在靜水條件下水平裂隙中的注漿擴(kuò)散過程，建立恒定注漿速率條件下考慮漿液黏度時(shí)空變化的水平裂隙注漿擴(kuò)散理論模型，推導(dǎo)漿液擴(kuò)散區(qū)內(nèi)的黏度及壓力時(shí)空分布方程，進(jìn)而得到注漿壓力與注漿時(shí)間及漿液擴(kuò)散半徑的關(guān)系。
基于顆粒流PFC的巖體爆破破壞效應(yīng)模擬 ￥55
影響巖體爆破破壞工程的因素眾多，節(jié)理特征、鉆孔孔徑、應(yīng)力波峰值等。目前針對(duì)巖石爆破的數(shù)值模擬采用的計(jì)算手段有LS-dyna和離散元方法等。本算例采用顆粒流PFC對(duì)巖石爆破過程進(jìn)行模擬。
分別就單點(diǎn)爆破、單點(diǎn)增大爆破壓時(shí)、三點(diǎn)同時(shí)爆破、三點(diǎn)微差爆破這四種工況進(jìn)行了仿真計(jì)算。
首先建立模型，在邊界墻體的伺服功能下平衡模型：
刪除邊界墻體，對(duì)左右兩側(cè)邊界附近及底部邊界附近處的顆粒運(yùn)動(dòng)進(jìn)行約束模擬邊界，對(duì)邊界顆粒施加荷載，吸收掉入射的波動(dòng)能量，以模擬無限介質(zhì)：
單點(diǎn)爆破的結(jié)果和爆破壓力的波形如下：
單點(diǎn)爆破增大炮孔壓力后的結(jié)果
三點(diǎn)同時(shí)爆破結(jié)果如下：
三點(diǎn)微差爆破結(jié)果及爆破壓力的波形如下：
具體建模思路及完整代碼（含基本注釋）如下：
展開 基于LS-DYNA的準(zhǔn)二維巖體爆破裂紋模擬2 ￥39.9
炸藥在巖體中爆炸瞬間釋放出大量的爆炸能量，產(chǎn)生爆炸沖擊波和應(yīng)力波，以動(dòng)載的形式作用周圍巖體并使得巖石損傷和破碎。一般認(rèn)為，巖體的爆破裂紋是由兩種不同類型的動(dòng)力荷載作用下產(chǎn)生的:爆炸應(yīng)力波和爆炸氣體。首先是炸藥起爆后在孔洞周圍產(chǎn)生的應(yīng)力沖擊波的作用，其次是持續(xù)時(shí)間較長的爆炸性氣體的作用。當(dāng)炸藥起爆時(shí)，壓力急劇增大，爆破孔周圍巖體被壓碎形成破碎區(qū)，而周圍巖體受到拉應(yīng)力產(chǎn)生裂縫，形成裂隙擴(kuò)展區(qū)。本案例巖石模型采用272號(hào)材料*MAT_RHT，模型算法采用ALE(任意拉格朗日-歐拉算法），由材料自帶的損傷參數(shù)獲得巖體產(chǎn)生破碎去及擴(kuò)展裂隙過程，模擬結(jié)果如下
圖1 巖體斷面裂紋擴(kuò)展模擬結(jié)果與試驗(yàn)結(jié)果比較
圖2 巖體斷面裂紋擴(kuò)展及損傷變化模擬結(jié)果
本案例適用于研究爆炸、沖擊、侵徹動(dòng)力學(xué)的朋友，下面附上該模擬的K文件，大家有疑問可以在私信我，歡迎交流！
展開 