丁未

（貴州省水利水電工程咨詢有限公司 貴州貴陽 556600）

基于三維地質(zhì)模型的巖土工程有限元自動建模方法

丁未

（貴州省水利水電工程咨詢有限公司 貴州貴陽 556600）

巖土工程改造中，要建立三維地質(zhì)模型進行綜合評估，掌握區(qū)域地質(zhì)結(jié)構(gòu)分布規(guī)律，提前對地質(zhì)災害及工程建設提供科學的指導。隨著計算機技術普及化發(fā)展，三維地質(zhì)模型在巖土結(jié)構(gòu)勘察與分析中得到應用，通過構(gòu)建模型平臺實現(xiàn)了工程改造科學化。結(jié)合“三角形有限元法”特點，本文分析了基于三維地質(zhì)模型的巖土建模分析方式，對巖土施工改造提供科學的指導依據(jù)。

三維模型；巖土工程；有限元；建模

巖土環(huán)境對區(qū)域工程建設有直接性應想，不同巖土構(gòu)造會引起多方面的地質(zhì)問題，影響了工程改造作業(yè)的一體化。鑒于傳統(tǒng)地質(zhì)規(guī)劃與發(fā)展趨勢，地質(zhì)模型建設關系著巖土工程改造質(zhì)量水平，對現(xiàn)場施工作業(yè)流程控制起到了引導作用。早期模型分析僅限于平面層次，對多角度巖土結(jié)構(gòu)識別缺乏準確性，利用三維模型分析巖土結(jié)構(gòu)才是發(fā)展的必然趨勢。

1 三維地質(zhì)模型

“地質(zhì)災害”是指在自然或者人為因素的作用下形成的，具有破壞性性質(zhì)的地質(zhì)作用或現(xiàn)象，包括：崩塌、滑坡、泥石流、地震、火山、水土流失等，這些都對地區(qū)穩(wěn)定發(fā)展造成了嚴重威脅。為了更好地防治地質(zhì)災害，設定三維地質(zhì)模型可實現(xiàn)災害的綜合分析。從空間布局來說，三維模式是利用X、Y、Z等三軸坐標進行建模的平臺，從多個角度呈現(xiàn)研究對象的動態(tài)形象。相比于傳統(tǒng)二維模型分析，三維模型研究更具精確性特點，可在短區(qū)域內(nèi)對工程模型精細化管理，提高項目建設與分析的針對性。新時期“三角形有限元法”在巖土結(jié)構(gòu)改造中發(fā)揮重要作用，利用三維模型實現(xiàn)巖土建模及在線分析，對區(qū)域、規(guī)模、承載、厚度等指標進行全面研究，掌握了區(qū)域地質(zhì)結(jié)構(gòu)的實際情況，結(jié)合有限元建模分析處理對現(xiàn)場施工給予指導。

2 三角形有限元法的基本求解過程

2.1 假設單元的位移函數(shù)

巖土工程測量結(jié)果需要以圖像形式呈現(xiàn)出來，才能為巖土企業(yè)提供形象的數(shù)據(jù)分析，掌握符合巖土工程實況的開發(fā)平臺。面對傳統(tǒng)地質(zhì)分析模式存在的不足，新時期巖土工程測量要堅持科技創(chuàng)新原則，引入三維模型作為技術支撐平臺，設定更加多元化的巖土圖像處理中心，客觀地反映出巖土工程動態(tài)。如圖1所示為一端部承受豎向集中荷載作用的懸臂梁，將此問題簡化為二維平面問題，梁取單位厚度，固定端簡化為兩個不動鉸支座，已知懸臂梁的長度為8m，高度為1m，材料為鋼材，它的彈性模量取2.06e11Pa，泊松比取0.25，右端處作用的豎向荷載為100kN。二維域Ω被離散成如圖1所示的九個三角形單元，后續(xù)編程中的結(jié)點號和單元號都以此為基礎，任取一三角形單元，設其結(jié)點編碼為i，j，m，以逆時針編碼為正向（如圖 2）。

圖1 任意區(qū)域的三角形單元離散

圖2 3結(jié)點三角形單元

該三角形單元在x、y方向的位移函數(shù)表達式為：

其中Ni（i=i，j，m)稱為單元的插值函數(shù)或形函數(shù)

上式（i，j，m）表示下標輪換，即 i→j，j→m，m→i。

（1）和（2）式簡記為：

（5）式的矩陣形式是：

2.2 單元應變矩陣和彈性矩陣

為了避免地質(zhì)分析不當造成數(shù)據(jù)結(jié)果失誤，可引入虛擬模型作為前期操控處理平臺?？偨Y(jié)地質(zhì)分析技術缺陷，以數(shù)字模型為虛擬化控制平臺，提出更為先進的地質(zhì)分析技術方案。虛擬模型為用戶提供動脈地質(zhì)分析平臺，本質(zhì)上并非正式的地質(zhì)分析操作，而是在虛擬環(huán)境下對圖像進行多種處理。因此，地質(zhì)分析采用虛擬模型具有精確化特性，避免正式操作中產(chǎn)生錯誤，提高了巖土圖像的處理質(zhì)量。本次將位移函數(shù)（6)式代入平面問題的幾何方程得到單元的應變表達式

其中B稱為應變矩陣：

將（3）式代入上式得：

對于平面應力問題，單元的彈性矩陣如下：

2.3 整體平衡方程

根據(jù)最小勢能原理或變分原理，對總勢能取駐值，即δπ=0得：

其中，單元剛度矩陣：

單元載荷列陣：

本文分析的問題只是在懸臂梁右端某一結(jié)點處作用一豎向荷載，荷載列陣相對比較簡單，如果涉及到分布荷載，需要局部坐標到整體坐標的轉(zhuǎn)換，形成整體荷載列陣。

2.4 組集單元矩陣到整體方程

按對號入座的原則對單元平衡方程進行組集得到整體求解的結(jié)構(gòu)控制方程：

2.5 引入邊界條件

三維模型在地質(zhì)分析系統(tǒng)中，與安全管理系統(tǒng)交互并協(xié)助控制管理地質(zhì)分析。例如，數(shù)字勘察功能隱患將對巖土結(jié)果產(chǎn)生誤導作用，影響了數(shù)字勘察結(jié)果的準確性，操作層可根據(jù)用戶指令狀態(tài)，及時調(diào)整巖土測量三維模型的處理模式。例如，邊界條件包括外力邊界條件和位移邊界條件，對于本文分析的問題，外力邊界條件主要是在右端點處有一個y方向的集中荷載。位移邊界條件主要是左端的兩個結(jié)點處在x、y兩個方向位移為零。

2.6 求解整體離散方程

地質(zhì)分析是我國巖土行業(yè)不可缺少的工程要素，巖土工程測量必須依賴于各項科技術應用，才能創(chuàng)造更為豐厚的產(chǎn)業(yè)收益。地質(zhì)分析是一項復雜的工藝流程，以圖像模型為中心執(zhí)行處理方案，從根本上解決了傳統(tǒng)地質(zhì)分析的不足。例如，在求出單元剛度矩陣組集而成的整體剛度矩陣和單元等效荷載列陣以后，求解整體離散方程就變成求解線性方程組的問題了。經(jīng)過組集單元矩陣到整體矩陣形成14式中的整體剛度矩陣K以及整體荷載列陣F。

3 巖土工程監(jiān)理工作內(nèi)容

隨著巖土工程建設不斷發(fā)展，三維模型技術用于地質(zhì)分析成為行業(yè)發(fā)展趨勢，將其作為主導技術可擴大巖土工程測量范圍，構(gòu)建更為穩(wěn)定的巖土測量平臺。

3.1 數(shù)字勘察

在數(shù)字勘察技術應用階段，要考慮巖土區(qū)域?qū)嶋H巖土情況，同時構(gòu)建數(shù)字化勘察模型平臺，及時掌握巖土工程構(gòu)造特點，借助三維模型技術獲取更多的巖土測量數(shù)據(jù)，指導巖土數(shù)字勘察工作有序進行。未來更好地完成項目建設，結(jié)合有限元自主建模分析法，總監(jiān)理工程師必須掌握本工程施工階段的安全控制要點，并讓監(jiān)理人員加以了解，檢查施工單位安全管理人員是否合格。

3.2 災害治理

結(jié)合三維地質(zhì)模型分析結(jié)果，對崩塌、滑坡、泥石流、地震、火山、水土流失等常見地質(zhì)災害進行綜合分析，及時采取針對性的治理方案，提高災害防控與治理效果。一般來說，巖土監(jiān)理與施工是并行操作的，在施工過程中指導災害防治施工，在防治施工中配合監(jiān)理工作。施工單位要結(jié)合具體災害類型，采取地基加固、邊坡支護、樁基礎改造等方式，做好地質(zhì)災害治理工作。

4 結(jié)論

總之，巖土工程涵蓋了地上、地下和水中的各類工程，對地區(qū)現(xiàn)代化建設與改造具有重要意義。為了擺脫傳統(tǒng)施工流程的不足，利用三維地質(zhì)模型進行有限元分析，借助建模平臺掌握巖土構(gòu)造特點，科學指導了區(qū)域工程改造方案，降低了各種因素帶來的風險隱患?；谌S地質(zhì)建模平臺下，巖土工程有限元建模分析精度更加準確，比傳統(tǒng)二維模型分析更具可視化、精確化等特點，值得在未來巖土工程分析中普及使用。

[1]張芳，朱合華，寧民霞.適于海量數(shù)據(jù)的三維地層建模方法[J].巖石力學與工程學報，2016（S1）.

[2]徐能雄，武雄，汪小剛，賈志欣，段慶偉.基于三維地質(zhì)建模的復雜構(gòu)造巖體六面體網(wǎng)格剖分方法[J].巖土工程學報，2016（08）.

[3]李金國，王 慧，劉 紅.基于Hypermesh的有限元前置優(yōu)化設計[J].機械工程師，2005（07）.

[4]夏艷華，白世偉，倪才勝.某水利樞紐廠房開挖三維可視化與數(shù)值模擬耦合研究[J].巖土力學，2015（06）.

[5]王純祥，白世偉.三維地層信息系統(tǒng)與有限元方法集成研究[J].巖石力學與工程學報，2014（21）.

TU452

A

1004-7344（2016）17-0156-02

2016-5-29

丁 未（1985-），男，侗族，助理工程師，碩士研究生，主要從事巖土工程方面的工作。