鄭文華

(中國建筑科學研究院地基基礎研究所，北京 100013)

1 前言

隨著計算機技術的廣泛應用，巖土工程數值分析軟件在大型水利、土建、礦山、交通、市政、石油等工程中的設計分析中起著越來越重要的作用。在巖土工程數值分析過程中，為了更好的模擬巖土體之間的接觸問題不可避免的要用到接觸單元［1］。

Interface 作為FLAC3D 中的接觸面單元，是FLAC3D 中自帶的一種無厚度單元，其本構模型采用庫侖剪切模型，可用于分析一定受力條件下兩個接觸的表面上的錯動滑移、分開與閉合，能很好地模擬巖錨梁體與周圍巖體的接觸變形［2］。

可用于模擬節(jié)理、斷層以及不同塊體間的接觸面等界面的滑移或分離問題。FLAC3D 中已有的interface 建模方法［3］“分離法”和“移動組合法”難以滿足復雜模型的建模需要，因此，國內學者在實際模擬分析中也探討了部分新的interface 建模方法，如文獻［4］中提到的“依次建模法”。但是該方法仍只局限于接觸問題的簡單分析，即在模型建立之初就必須進行interface 單元的建模，且僅適用于可以依次建模的分析問題。

對于類似隧道開挖支護過程，接觸面由無到有的全過程復雜模型interface 建模方法，國內文獻中至今仍無詳細論述。因此，本文從FLAC3D 建模方法的角度出發(fā)，探討了適用于巖土工程全過程數值模擬分析的復雜模型interface 建模方法，以期為今后的FLAC3D 中interface 的建模方法提供一種參考。

2 interface 單元已有的建模方法

2.1 FLAC3D 軟件簡介［5］

2.2 interface 界面單元

如圖1 所示，在FLAC3D 中interface 接觸單元是由一系列三角形單元組成的，每個單元由3 個節(jié)點來定義。一般而言，interface 接觸單元是依附于固定的面，兩個三角形單元可用來定義任意四邊形的接觸面。在每個接觸單元的頂點處會自動生成接觸節(jié)點。接觸節(jié)點會自動探測接觸的發(fā)生過程，并通過法向剛度、剪切剛度或滑動性質等來表達這種接觸關系。具體表達方式見圖2。

圖1 接觸單元模型

圖2 接觸單元本構模型［3］

其力學方程為:

2.3 FLAC3D 中已有的interface 單元建模方法

(1)分離法

先建立兩個完整的模型，將需要建立接觸面的模型網格分離開，最后設置接觸面單元。

(2)移來移去法

將有接觸面的兩個模型分開建立在不同的位置，然后在其中一個模型上建立接觸單元，最后將另一塊模型移動到該模型的對應位置，實現接觸關系的建立。

圖4 移來移去法模型示例

(3)依次建模法［4］

加息對股市到底屬于利好還是利空？巴菲特曾經說過：“即使我左邊坐著美國財政部長，右邊坐著美國聯(lián)邦儲備局局長，他倆偷偷告訴我今年美國將減多少稅和加多少息，我也不會根據他們的信息去買賣股票?！蓖顿Y大師在美國本土尚且如此，何況是遠隔重洋呢？有效市場假說理論認為，所有消息都會被提前消化，其影響難以確定。事實上，僅憑有限的分析數據，真很難界定單個因素在所有驅動因子中的權重。

按照先后順序，先建立模型一，然后進行接觸面單元的建立與設置，最后建立另一塊模型。

以上三種方法都可以建立簡單模型的interface單元，但是對于大型曲面的interface 單元的建立則難以實現。

圖5 依次建模法模型示例

3 復雜模型中interface 建模方法

針對復雜模型的接觸面單元建立過程困難的問題，需借助于其他的FLAC3D 前處理軟件在接觸部位先建立曲面，然后進入FLAC3D 的操作界面，進行interface 單元的建立，即“先實體面后接觸面法”。

3.1 Ansys 中復雜模型實體面的建立

利用Ansys 的強大功能，在需要進行接觸面模擬的界面兩側建立體模型，劃分網格。給接觸面兩側的體模型賦不同的材料屬性，隨后將模型信息導入FLAC3D 中。在此過程中最好保證接觸面部位存在實體面。

圖6 ansys 中的實體面

3.2 FLAC3D 中接觸面單元(interface)的建立

將Ansys 中已建好的模型導入FLAC3D 中，可先進行接觸問題之前的其他力學問題的模擬分析。當需要進行接觸分析時，針對需建立interface 單元的模型:

(1)進行接觸面兩側模型的分組

語句:gro contact1 range gro 5

gro contact2 range gro 6

(2)將接觸面所在的模型分離

語句:gen separate contact1

(3)進行接觸面單元的設置

語句:interface 1 wrap contact1 contact2

在完成接觸面模型的建立后，進行相關參數設置及下一步工作。

4 應用實例簡介

4.1 工程簡介

某隧道工程，采用盾構施工法。上部為軟土層，下部為全風化硬層。在模擬隧道開挖施工的過程中，需考慮盾構管片之間以及管片與圍巖之間的接觸問題，需建立interface 接觸單元［6-10］。

4.2 地鐵隧道模型

首先利用Ansys 軟件建立隧道模型，建模過程中保證了在接觸面處存在實體面。如圖7 所示。

圖7 隧道模型

4.3 interface 界面單元建模

在模型導入FLAC3D 軟件以后，對于接觸面兩側的模型進行分組(在圖中以不同顏色表示)，然后進行模型分離以及接觸單元的設置。建模結果如下:

(1)管片縱向接觸面

圖8 縱向接觸面模型

(2)環(huán)向接觸面

圖9 環(huán)向接觸面模型

5 結束語

本文從FLAC3D 軟件的實際運用角度出發(fā)，在總結現有interface 建模方法的基礎上，提出了基于Ansys 的FLAC3D 復雜模型interface 建模方法。同時，將該方法運用到實際的工程數值計算過程中，驗證了該建模方法的有效性、實用性，對今后的FLAC3D 的接觸方法建模具有重要的參考價值。

［1］劉波，韓彥輝.FLAC 原理、實例與應用指南［M］.北京:人民交通出版社，2005.

［2］馬震岳，凌影，金長宇，江權.某水電站地下廠房巖錨吊車梁三期澆筑可行性分析［J］.水電能源科學，2010，28(3):66-69.

［3］FLAC3D User's Manual.

［4］姜慶紅，范慧鶴，李紅英，李利彬.FLAC3D 中INTERFACE 建模新方法在采空區(qū)模擬中的應用［J］.計算機與現代化，2008，(10):113-115.

［5］彭文斌.FLAC3D 實用教程［M］.北京:機械工業(yè)出版社，2007.

［6］黃建明.盾構管片計算模型的選擇［J］.鐵道建筑，2004，(6):29-31.

［7］鄧飛皇.地鐵營運期盾構隧道管片位移響應分析［J］.科學技術與工程，2009，9(13):3684-3688.

［8］劉遠亮.建筑物施工隊地鐵隧道影響的三位數值模擬分析［J］.中國水運，2011，(6):255-259.

［9］陳立，趙敏，方俊林.軸向荷載作用下灌注樁數值模擬［J］.河南工程學院學報(自然科學版)，2010，22(1):25-27.

［10］曾亞武，劉繼國，蔡元奇.深基坑開挖引起的土體變形模擬［J］.武漢大學學報(工學版)，2006，39(5):53-56.