基于FLAC3D復雜三維地質建模

潘曦 朱斌

摘要：真實地質體三維數值模型構建是進行巖體工程數值分析面臨的難題，開展地質體三維數值模型構建方法比較研究具有重要意義。以FLAC3D數字化模型為基礎，提出了FLAC3D-Rhino耦合建模方法，詳細闡述了各建模方法具體步驟。構建了桂林市五美路停車場商住樓深基坑土方開挖及工程三維數值模型，分析了深基坑上開挖地層位移和水壓的變化，研究成果對準確構建地質體三維數值模型具有重要指導作用。

Abstract： The construction of three-dimensional numerical model of real geological body is a difficult problem for the numerical analysis of rock mass engineering. It is of great significance to carry out a comparative study of the construction method of three-dimensional numerical model of geological body. Based on the FLAC3D digital model， a FLAC3D-Rhino coupling modeling method is proposed and the detailed steps of each modeling method are elaborated. Three-dimensional numerical model of excavation and engineering for deep excavation of commercial parking building of Wumei Road parking lot in Guilin City was constructed. The change of excavation stratum displacement and water pressure in deep excavation was analyzed. The research results have important guiding role in accurately constructing three-dimensional numerical models of geological bodies.

關鍵詞：三維數值模型；地質體；多軟件耦合；Extrusion

Key words： three-dimensional numerical model；geological body；multi-software coupling；Extrusion

中圖分類號：P628 文獻標識碼：A 文章編號：1006-4311（2018）25-0188-05

0 引言

在工程領域中，有限元、有限差分等數值模擬技術的應用越來越廣泛。當前，在FLAC3D的應用實踐中逐漸顯現出了一些弊端，突出表現在建模工作量較大、軟件前處理功能薄弱、難以有效構建復雜的曲面、三維模型建成后不便修改及構建的數值模型與真實地質體嚴重不符等方面，嚴重影響了計算結果的可靠性[1]。

基于計算機構建的巖體可視化模型，構建三維數值計算模型是大勢所趨。當前巖土、地下水工程常用數值模擬軟件主要有ADINA、ABAQUS、ANSYS、MIDAS、FLAC等，但上述軟件均存在較難構建復雜地質體三維數值模型、靈活性差、前期建模能力較弱等問題。基于此，諸多學者針對上述問題展開了研究，鄭文棠等[3]結合AutoCAD與FLAC3D及3DEC計算功能，構建了復雜邊坡三維地質可視化和數值模型。田樹昆等[4]提出了使用AutoCAD二維轉三維程序生成三維模型線框的方式構建ANSYS數值計算模型。李翔等[5]利用SURPAC、AutoCAD和ARCGIS等軟件生成地質體網格點坐標，并經ANSYS軟件網格化后構建了FLAC3D數值模型。廖秋林等[6]通過分析FLAC3D與ANSYS單元數據關系，運用VB語言編寫了FLAC3D-ANSYS接口程序。馬新根等[7]通過KUBRIX插件實現了Rhino3D模型與FLAC3D的對接。崔芳鵬等[8]通過FLAC3D內嵌Fish語言，編制了可用FLAC3D直接讀取Surfer軟件數據的程序。孫濤等[9]、許國等[10]利用GOCAD軟件建立的地質體三維可視化模型，通過接口程序在ABAQUS、FLAC3D軟件中構建了三維數值模型。郭穎星等[11]通過構建Petrel與ADINA的接口程序，實現復雜網格數據的快速移植；王樹仁等[12]利用Matlab編寫MIDAS/GTSFLAC3D的接口程序，實現了MIDAS/GTS到FLAC3D的轉化。類似上述多軟件耦合構建三維數值模型的方法較多[13，14]，但構建的數值模型網格質量不高，與真實地質體賦存產狀對應性差，且設計工程過程與數值模擬過程差異大。由此可知，FLAC3D的建模優化還有較大空間，理想建模環境應具備以下幾項特點：

①建模過程直觀，使用空間交互式3D建模方法建模，禁止少用或者不用編程命令。

②建模操作簡單，不僅建模難度和工作量降低了，也更便于使用者快速掌握。

③可二次編輯已構建的三維模型，且網格單元劃分可控。

④所建三維模型兼容性強，可直接利用或修改同類三維建模軟件。

針對這些特點進行優化，本文將地質體三維數字化模型與Rhino、FLAC3D等軟件耦合，開展復雜地質體三維數值模型構建方法比較研究，并結合實例來進行驗證。

1 復雜地質體三維數字化模型構建

1.1 工程背景

本工程系五美路停車場商住樓，本工程位于桂林市中心象山區榕湖附近，北臨五美路，南臨臨桂路，東靠交通路，西至桂林市規劃局。本項目屬于一類高層建筑。設計使用年限為50年，抗震設防烈度為六級。地上部分為商業，辦公，住宅綜合體，地下負一至三層均為停車庫、設備用房。建筑高度33.8m。建筑總面積為48103.07m2。商業總建筑面積為10963.18m2，住宅總建筑面積為19216.49m2，辦公總建筑面積為11170.03m2，配套公建總建筑面積為540.23m2。五美路停車場商住樓，嵌固端為地下室頂板，地下3層，地上分縫形成3個塔樓，采用框支-剪力墻結構，建筑高度為33.8米，本工程±0.000絕對標高152.780米。

工程概況如表1。

地下室概況：

原地面標高為-1.3m，地下室平面為不規則圖形，長106.2米，寬80.4米，平面面積約5618.23平方米。3層地下室層高（相對標高）自上而下依次為4.3m（-5.6m）、3.2m（-8.8m）、3.25m（-12.05m），筏板厚度分650mm和2000mm，電梯井位置筏板厚度最大為2000mm，底板開挖面相對標高為-12.8～15.8米。

①本工程采用逆作法施工，要求逆作鋼管立柱樁垂直度偏差控制在1/500以內，且鋼管頂標高低于地面1.3m，鋼管立柱樁施工難度較大。

②本工程的出土工程量大、開挖深度深，因為采用逆作法，出土效率較低。

③土方開挖形式：

若基坑的土方和混凝土結構工程量較大，其所需工期會相應增加，客觀上增大了基坑風險。為了有效控制基坑變形，可利用“時空效應”，將基坑上方開挖和主體結構劃分施工段并采取分塊開挖的方法。

周邊環境：

北臨五美路，北面無建筑物，且與人行道相鄰。支護樁距道路中線距離大約20米左右，市政管線包括路燈電線、給水管、排水管、電纜等。

南為臨桂路，臨桂路路面狹小，市政管線包括路燈電線、給水管、排水管、電纜等。

東靠交通路，并有一棟7層住宅樓一棟九層住宅樓，兩層住宅樓都為框架剪力墻結構，支護樁中心到建筑物地面磚間距最寬為9.41米，最窄距離為2米，市政管線包括路燈電線、給水管、排水管、電纜等。

西臨桂林市城鄉建設局。支護樁中心線距離建設局建筑物為9.28米，且距建筑物6.76米的位置為70年代老式建筑。建筑結構為老式磚混結構，基礎為淺式基礎。

1.2 基于Rhino與Flac3D Extrusion的三維數字化模型構建

現有的有限元、有限差分等大型數值模擬軟件存在前處理功能薄弱、較難構建復雜地質體的三維精細化數值模型等不足，要想建立復雜地質體的三維精細化數字模型，三維建模軟件不可或缺，將三維數字化模型構建方法與數值模型構建方法耦合，建立三維精細化數值模型是必然趨勢。本文利用FLAC3D Extrusion軟件構建地質體三維數字化模型，并通過Rhino軟件與該模型的耦合，開展地質體三維數值模型構建方法研究。

利用FLAC3D軟件分別構建支護樁、止水帷幕、地層等的三維可視化實體模型，如圖1所示。

2 復雜地質體三維數值模型多軟件構建方法

2.1 FLAC3D-Rhino 耦合建模方法

2.1.1 基于FLAC3D Extrusion 塊體模型的構建方法

FLAC3D5.0版本新增添了Extrusion功能，此功能使用擠壓能力生成網格應用于FLAC3D中。是一套通過2d形狀（畫）線性擴展（擠壓）形成第三維度的建模方法。

具體操作為：

①建立模型輪廓圖。在cad中畫出模型輪廓平面圖，這個平面圖相當于要建立模型的橫剖面圖。另存為dxf格式。然后讀入模型輪廓圖。打開FLAC3D5.0版本軟件，新建工程，這個工程是你整個工程的名字，之后打開菜單欄中Panes，點擊New Extrusion，起工程名后，點擊下拉框右邊第一個按鈕及Import，選中之前建立的*.dxf。

②描圖。在flac3d5.0版本中，通過point-Edge Tool來描圖，兩點確定一直線。若是曲線的話，在曲線兩端連接一條直線后，用Contrl point tool控件，將直線中點和要重合的曲線中點重合即可，這時可以看到直線有三點和曲線重合了，這時單擊這個折線，會出現Object properties，在里面的Eage type中選擇curv就可以使這個折線與曲線完全重合了。也可以直接在Extrusion進行1：1構圖。

③網格劃分。本文選擇按線上單元體個數劃分。做法是單擊Autozone控件，選擇下拉列表中model extent，進行調整，點擊確定。再點擊每條線，在object properties中可以改變zones的個數，劃分成合適的網格大小。

④分組。在construction view控件中進行分組，點擊后的視圖中用show groups控件及相應的object properties調整組。用ctrl+左鍵點選在object properties填寫組名之后回車分組，如圖2所示。

⑤拉伸模型。這個概念和cad中概念相似。點擊extrusion控件，可以用point edge tool控件在底部新建幾個直線，若拉伸方向不用分組那么就做這個直線了。之后相應屬性（比如說長度，每份個數都可以在相應的object properties里調整）。拉伸完成后，生成單元體。點擊Extrude按鈕，生成的模型，如圖3所示。

2.1.2 基于FLAC3D 面模型的構建方法

FLAC3D 5.01軟件具有網格點（Grid）通過（Export）功能進行導出。運用文本格式（.txt）方式打開。得知FLAC3D 5.0網格單元的數據形式由4個部分組成。①表示模型生成的時間；②* GRIDPOINTS（點）：表示節點的生成，格式為：節點標志、節點序號、節點（x，y，z）坐標；③* ZONES（單元）：表示單元的生成，格式為：單元標志、單元類型、單元中的節點號；④* GROUPS（分組）：表明單元所屬的組。FLAC3D 遵從的是點（GRIDPOINT）、單元（ZONE）、組（GROUP）自下而上的網格建立模式，即在建立實體模型的同時，軟件自動完成該實體的部分，并以點、單元和組的形式保存下來。Surfer、Rhino 等軟件能夠較好地構建復雜地質體模型，且各軟件在模型構建方面存在各自優勢。上述各軟件間存在一定聯系，可以通過特定格式文件進行相互轉換。本文基于FLAC3D Extrusion三維可視化模型包含的各層面點信息，通過Surfer軟件將點信息轉入Rhino 軟件，生成地表面和地層面，利用FLAC3D程序講構建出的地質體三維數值模型進行計算。也就是說，構建模型真實地貌的實質就是構建每個地層的面。對于本文模型而言，構建真實地貌的地層面相當于仿真其真實高程值，究其根本為模型文件FLCD3格式中GRIDPOINT節點z軸的坐標。為此，本文將勘察報告中剖面圖、柱狀圖和平面圖，運用EXCEL表對勘探點地層高程的層面點信息分類輸入。利用Surfer軟件參考。Rhino軟件形成曲面，導出坐標點更改.flac3d格式文件的z軸坐標進行仿真地層模型的構建。多軟件建模耦合流程圖，如圖4所示。

具體步驟依次如下：

①基于FLAC3D Extrusion通過工程實際要求構建三維模型，將層面體按一定X，Y間距生成點集，見圖5。將所有點坐標信息導出.flac3d格式，并整理為Excel可識別的格式。

②利用Surfer軟件對真實地層高程點信息進行插值，生成“.Grd”格式的等高線文件，然后導入等高線文件輸出地層圖。

③利用Excel軟件處理“.flac3d”格式的模型文件，提出模型一個層面的x，y，z坐標并導入Rhino軟件中形成平面a。再將一層真實地層高程的點信息（如：地表層）通過文件讀取Excel表的層面信息選擇選取起始面（平面a），參考之前Surfer形成的地層圖，對曲面的點間距、U和V方向跨距數、硬度、起始曲面拉力進行調整，形成地層曲面。并對地層曲面提取點坐標，將此曲面的x，y，z軸坐標進行導出并保存至EXCEL表。如圖6、圖7。

④通過網格x，y一一對應原則，進行z坐標的修改。將修改好的曲面坐標替換掉之前模型.flac3d格式文件中，從GRID進行導入，形成一層仿真地層地貌。如圖8。

⑤依據此法修改兩層地層形成模型體，如圖9。

⑥通過層層疊加原理建立仿真地層模型，如圖10。

3 結語

本文在研究基于多軟件與Flac3D耦合的三維數字化模型構建，結合桂林市某深基坑工程的三維數值模型。通過模擬模型可知，此方法能夠保證實驗數據分層信息與結果數據的分層信息相匹配，基本實現了三維數值建模，對實際的生產應用有一定的參考價值。但是，采用本文所提出的方法，精確度仍有進步的空間。

參考文獻：

[1]鄒力，彭雄志.淺談FLAC3D的應用原理、優缺點及改進措施[J].計算機應用，2007（2）：152-156.

[2]劉曉明，羅周全，楊彪，等.復雜礦區三維地質可視化及數值模型構建[J].巖土力學，2010，31（12）：4006-4015.

[3]ZHENG Wen-tang， XU Wei-ya， TONG Fu-guo， et al. 3Dgeological visualization and numerical modeling of complicated slope[J]. Chinese Journal of RockMechanics and Engineering， 2007， 26（8）： 1634-1644.

[4]TIAN Shu-kun， CAO Lan-zhu. The simulation analysis of slope by FLAC3D based on CAD and ANSYS[J].Microcomputer Information， 2010， 16（11）： 259-260.

[5]李翔，王金安，張少杰.復雜地質體三維數值建模方法研究[J].西安科技大學學報，2012，32（6）：676-681.

[6]LIAO Qiu-lin， ZENG Qian-bang， LIU Tong， et al.Automatic model generation of complex geologic body with FLAC3D based on ANSYS platform[J]. Chinese Journal of Rock Mechanics and Engineering， 24（6）：1010-1013.

[7]MA Xin-gen， SUN Jian-dong， ZHAO Hong-ze， et al.Application of FLAC3D rapid modeling in slope problem of open-pit mine based on Rhino3D NURBS[J]. Safety in Coal Mines， 2015， 46（8）： 138-141.

[8]CUI Fang-peng， HU Rui-lin， LIU Zhao-lian， et al. Surfer software platform based complex three-dimen sional geological digital models for preprocessing of FLAC3D[J]. Journal of Engineering Geology， 2008， 16（5）： 699-702.

[9]孫濤，劉亮明，趙義來，等.基于GOCAD 平臺的復雜地質體系的動力學建模研究[J].礦產與地質，2011，25（2）：163-167.

[10]許國，王長海.萬家口水電站復雜地質體三維模型及其數值模型構建[J].武漢大學學報（工學版），2013，46（4）：469-474.

[11]郭穎星，祝愛玉，孫豐成，等.復雜斷層數值建模及程序接口實現[J].地震地磁觀測與研究，2016，37（2）：90-97.

[12]WANG Shu-ren， ZHANG Hai-qing. A coupling modeling method with MIDAS/ GTS-FLAC3D and its application[J]. Journal of Civil， Architectural & Environmental Engineering， 2010， 32（1）： 12-17.

[13]BOCK S. New open-source ANSYS-SolidWorks-FLAC3D geometry conversion programs[J]. Journal of Sustainable Mining， 2015， 14： 124-132.

[14]YANG K M， WANG L B， ZHANG T， et al. Novel modeling on numerical computing the geo-deformation information in coalmine based on the GIS-Excel[J]. High Technology Letters， 2013， 19（2）： 156-161.