FLAC3D中隨機裂隙模型實現技術研究

翟會超,任鳳玉,張國建

(1.東北大學，遼寧 沈陽 110006；2.遼寧科技大學，遼寧 鞍山 114051 )

1 前言

在工程巖體中常存在節理裂隙。而裂隙弱面的發育程度、分布密度、力學性質等決定著結構體的形狀、方位和大小，同時也是控制巖體穩定性、巖體強度的重要因素。在FLAC3D數值分析中，常用Fish語言對模型實體內不同空間位置的單元賦予隨機材料參數，或通過雙線性應變硬化模型、多節理模型、分界面處理節理等方法來實現對巖體非線性研究。但以上方法較難描述巖體工程中賦存的宏觀隨機裂隙分布狀況。然而，巖體的變形破壞特征與宏觀裂隙密切相關。針對這一問題，文中采取Autocad、Ansys聯合建模技術，實現宏觀隨機裂隙模型在FLAC3D中的數值分析。

2 隨機裂隙模型實現技術

若生成線型宏觀隨機裂隙模型，首先對工程巖體進行結構面調查分析，獲取優勢結構面組的統計資料，進而確定宏觀隨機裂隙產狀、分布密度等空間特征，運用VB程序完成隨機裂隙面組在Autocad中的實現；其次，利用Ansys強大的建模功能，對宏觀隨機裂隙模型按照研究的工程范圍進行整體優化設計；最后，通過轉換程序將Ansys模型導入FLAC3D中進行數值模擬。主要建模技術流程見圖1。

每個裂隙的空間分布主要受傾角、傾向、裂隙長度、空間坐標四個參數控制，而隨機裂隙群組在此基礎上還要考慮裂隙的分布規律、密度等。因此，為確立空間隨機裂隙群，就要以裂隙的產狀、分布規律等相關參數為自變量，以控制裂隙方位的空間坐標、裂隙面總數為因變量。運用VB語言按照變量之間的函數關系編寫隨機裂隙面群生成程序，調用Autocad，并在其中實現隨機裂隙空間分布。

圖1 隨機裂隙模型建模過程

文中裂隙的形狀為簡單的正四邊形，利用四邊形角點的坐標便可控制裂隙的空間位置。隨機裂隙的空間分布特征以現場實際結構面調查、分析優勢結構面組為依據。生成的裂隙全部為優勢結構面組類，在空間位置上均勻分布，在產狀特征上以優勢結構面組的研究結果為根據，即文中生成的隨機裂隙是調查分析后的優勢裂隙面。

根據裂隙空間位置與裂隙面產狀、尺寸的關系，即：

Ψ(X,Y,Z)&Φ(dd,dip,l)

式中：dd為傾向；dip為傾角；l為裂隙邊長。

由點生線、由線生面的方法畫出隨機裂隙面，再循環n次得到n條隨機裂隙。

依據工程巖體中實際宏觀裂隙的分布條件，使用已編程序在Autocad中生成三維空間隨機裂隙面組，見圖2。

其次，把cad圖形文件保存為dxf格式，運用

轉換程序形成Ansys識別的命令流lgw格式文件，并調用lgw格式文件在Ansys里形成隨機裂隙組模型。然后利用Ansys獨特的強大建模技術進行整體模型設計，在通過布爾運算進行宏觀隨機裂隙的挖除，進而劃分網格。

圖2 Autocad中生成的隨機裂隙面

最后，把Ansys網格單元模型中的節點及單元數據以lis類型文件導出，運用已編程序將lis文件轉換成FLAC3D可以識別導入的flac3d類型文件。通過命令：

impgrid

將網格模型調入FLAC3D進行數值分析。

3 簡單隨機裂隙模型舉例

以簡單巷道開挖為例，試驗宏觀隨機裂隙模型在FLAC中的模擬應用。本模型尺寸為：40m×30m×20m的立方體，當中有一近4m×3m尺寸的拱形巷道。隨機裂隙模型總體節點數22860個，單元數量129107個；非裂隙模型總節點數1467個，單元數量1135個。模型采用摩爾庫侖塑性破壞準則進行分析。材料參數見表1。

表1 巖石材料參數

模型依靠重力平衡，然后再挖去巷道范圍應力場二次重分布引起巖體塑性變形，巷道開挖后兩種模型在y=15m剖面的塑性區范圍見圖3。

圖3 塑性區分布圖

圖5中巷道開挖引起圍巖受拉破壞，兩種模型巷道頂板及底板受開挖影響的破壞范圍基本一致。但隨機裂隙模型在裂隙周圍還存在拉剪塑性區，這加大了破壞面積，說明現實工程巖體的穩定性受裂隙分布的影響很大。根據斷裂理論，巖體內存在的節理裂隙、空隙等弱面，在一定程度上降低了工程結構的承載能力，易使工程失穩破壞。

4 結論

根據工程巖體內存在裂隙的特點，研究應用VB程序對autocad進行二次開發，在其cad繪圖窗口生成隨機裂隙面組，然后將其轉入Ansys中進行模型整體設計，最后導入FLAC3D中進行模擬分析。從裂隙模型應用結果來看隨機裂隙的存在使得巖體內塑性范圍加大。根據斷裂理論，巖體內存在的節理裂隙、空隙等弱面，在一定程度上降低了工程結構的穩定性。同時，隨機裂隙模型比非裂隙模型的單元網格數量大數十倍，這是因為隨機裂隙模型不易按規整方式劃分網格的緣故，這將耗費掉一定的計算機內存。

[1] 丁萬濤.隨機裂隙對節理巖體穩定性影響研究及其在海底隧道中的應用[D].濟南：山東大學,2008.

[2] 彭文斌.FLAC 3D實用教程[M].北京：機械工業出版社,2008.

[3] 張晉西. Visual Basic與AutoCAD二次開發[M].北京：清華大學出版社, 2002.

[4] 邢靜忠,李軍,等.ANSYS的建模方法和網格劃分[J].中國水運，2006.6(9):116-118.

[5] 李新平,朱瑞庚,夏元友.裂隙分布對地下硐室穩定性的影響研究[J].金屬礦山，1997(3):15-18.

[6] 程靳,趙樹山.斷裂力學[M].北京：科學出版社，2008.