有限元法在吉林省中部城市引松供水工程中的應用

王 倩，侯 薇，鞠 曄

（1.吉林省水利水電勘測設計研究院，吉林 長春 130000；2.水利部松遼水利委員會，吉林 長春 130021）

有限元法在吉林省中部城市引松供水工程中的應用

王 倩1，侯 薇1，鞠 曄2

（1.吉林省水利水電勘測設計研究院，吉林 長春 130000；2.水利部松遼水利委員會，吉林 長春 130021）

有限元法具有考慮巖土介質的非均勻性、各項異性、非連續性和材料與幾何非線性，且能適應于各種實際邊界條件的優點，在隧洞分析中比常規計算方法具有一定的優越性。本文針對吉林省中部城市引松供水工程中的有壓隧洞，采用有限元法計算在不同圍巖類別條件下開挖隧洞所引起的位移變形、應力和應變值，分析圍巖穩定性，為設計、開挖及維護隧洞提供可行性方法。

有限元；有壓隧洞；圍巖穩定性；引松供水工程

1 概 述

吉林省中部城市引松供水工程初步擬定從第二松花江流域豐滿水庫壩上左岸取水，途經溫德河、岔路河、飲馬河、雙陽河到達叢家嶺分水樞紐，分別由長春干線、遼源干線、四平干線及各附屬支線供水給長春市、四平市、遼源市及所屬的九臺市、德惠市、農安縣、公主嶺市、梨樹縣、伊通縣、東遼縣、長春雙陽區等11個市、縣、區的城區，以及供水線路附近可直接供水的25個鎮。

工程總體布局由1條總干線、1處分水樞紐、3條干線、干線3個調節水庫、13條支線、支線7個調節水庫和各線路上相應交叉及附屬建筑物組成。輸水線路總長550.6 km，其中輸水總干線起點為豐滿水庫左岸，終點為馮家嶺分水樞紐，全長110 km，有壓隧洞長91 km，占總干線長度的83%。

因此，如何合理、科學、精確地進行隧洞結構計算成為中部城市引松供水工程項目的關鍵問題之一。

2 有限元技術

2.1 方 法

隧洞屬于地下修筑的結構物。由于對地下結構的特性認識不充分，在設計方法上多數是沿用地面結構的設計方法，但是這種方法與實際情況相差很大。隨著科學技術的提高，人們對地下結構的認識有了進一步的深入，已經認識到地下結構是由周邊圍巖和支護結構兩者共同組成，并相互作用的結構體系，即地下結構＝支護結構+周邊圍巖。

隧洞的動靜力學計算是一項比較困難的課題，地層巖土介質和隧洞結構相互作用相當復雜。只有那些具有規則幾何形狀和理想的材料特性，且荷載形式與邊界條件是簡單的線彈性體系，才能得到較為精確的解答。但是，對于非線性巖土體內的連續或不連續介質和任意幾何外形的隧洞結構，其力學計算必須借助于近似的數值方法：有限元法、邊界元法、有限元—邊界元耦合法。

對于隧洞，主要關心的區域是隧洞附近，可用有限元法。有限元法的優點在于可以考慮巖土介質的非均勻性、各項異性、非連續性和材料與幾何非線性，且能適用于各種實際的邊界條件；缺點在于要求分析區域的幾何、物理連續性[1]。

2.2 采用的強度理論

非線性有限元方法分析隧洞在開挖過程中的變形、應力及應變，考慮了巖石本身的非線性本構關系，能夠模擬施工過程，可適用于任意復雜的邊界條件，但由于巖石本身的復雜性，導致采用不同的屈服準則對有限元計算結果有很大影響。

經典強度理論是針對無摩擦的金屬材料建立起來的，而且假設抗拉強度等于抗壓強度，故該理論對于巖石這種具有摩擦且拉壓強度相差懸殊的材料并不適用[2]。

目前巖土材料常采用莫爾-庫侖屈服準則（M-C）和廣義米賽斯屈服準則。莫爾-庫侖屈服準則較好地反映了巖土材料拉壓不等的特性，應用最為廣泛，但也存在諸多缺點，例如，它在三維應力空間中的屈服面存在棱角奇異點而導致數值計算不收斂。為此前人對其做了大量的修正，總體上看，這些修正準則將在π平面上的六角形屈服曲線抹圓，雖較好地解決了莫爾-庫侖屈服準則棱角不收斂的問題，但表達式往往過于復雜，不便于應用。與莫爾-庫侖屈服準則及其眾多的修正準則不同，廣義米賽斯屈服準則在π平面上是一個圓，且表達式簡單，便于數值計算。廣義米賽斯準則的一種德魯克-普拉格準則可由莫爾-庫侖準則基于關聯流動法則推出，在π平面上是莫爾-庫侖準則的內切圓。

以平面應變理想彈塑性為條件，并且基于非相關聯流動法則，也可以將莫爾-庫侖屈服準則轉化為具有相同的數學表達形式的廣義米賽斯屈服準則，推導過程省略。說明平面應變條件下莫爾-庫侖匹配的廣義米賽斯屈服準則在關聯流動法則和非關聯流動法則條件下都與莫爾-庫侖屈服準則形式一致。

采用德魯克-普拉格準則（D-P準則）既可以考慮中間主應力和靜水壓力對剪切屈服或強度的影響，也克服了莫爾-庫侖屈服準則（M-C）存在棱角奇異點而導致數值計算不收斂的缺點。可見，D-P準則作為有限元分析的強度理論是科學、合理的。

3 計算實例

吉林省中部引松供水工程主體工程91 km隧洞，底坡為1/6 000，設計引水流量38 m3/s，開挖直徑7 m。隧洞施工采用以TBM施工為主，鉆爆法施工為輔的聯合施工方法。前段61.8 km采用3臺掘進機施工，每臺控制在18～20 km。

3.1 計算模型

該工程隧洞大部分屬于長距離深埋隧洞，在對隧洞進行開挖過程的有限元模擬分析中，為了使模型更加接近實際，計算模型采用邊界尺寸是洞徑的7～8倍，巖體單元采用PLANE42來模擬，襯砌單元采用BEAM3來模擬，共9 527個單元。

考慮到隧洞施工過程的最不利情況，計算中采用全斷面開挖施工方法，對不加任何支護手段的裸洞進行圍巖穩定性分析，通過對洞室周邊圍巖的應力場、位移場以及塑性區的對比分析，從而揭示開挖隧洞周邊圍巖的應力應變規律，進而為隧洞進行合理支護設計做參考。

3.2 模型參數

根據隧洞圍巖的物理力學性質，在此次有限元計算當中，采用了彈塑性的非線性有限元法。圍巖材料的本構模型采用 Drucker-Prager（D-P）模型，以計算隧洞結構與地層在開挖過程中發生的非線性變形特性。采用的各類圍巖參數依據鉆爆施工方法確定，如表1所示。

3.3 計算成果及分析

以Ⅳ類圍巖的隧洞有限元分析為例，洞室開挖后的等效應力、水平位移及豎向位移見圖1～3，計算結果列于表2。

計算結果表明，隧洞開挖后Ⅱ，Ⅲ類圍巖拱頂變形很小，僅有1 mm左右，由于此次計算選取的Ⅳ類圍巖參數比較保守，隧洞開挖后圍巖拱頂變形不到5 mm，Ⅴ類圍巖則幾乎達到25 mm，隨著圍巖巖性力學性能的下降，等效應力值逐漸減小，而等效塑性應變值逐漸增大，到Ⅴ類圍巖達到11.56%，可見Ⅱ，Ⅲ類圍巖基本可以達到自穩，Ⅳ類圍巖頂拱兩側部位等效塑性應變達4.45%，說明Ⅳ，Ⅴ類圍巖均開始出現塑性變形，應重點對這些部位加強支護，以防發生崩塌事故。

表1 各類圍巖物理力學參數

圖1 Ⅳ類圍巖洞室開挖后等效應力圖

圖2 Ⅳ類圍巖洞室開挖后水平位移圖

圖3 Ⅳ類圍巖洞室開挖后豎向位移圖

4 結語

由于目前吉林省中部引松供水工程仍然處于可研階段，關于隧洞有限元分析也處于初探狀態，隧洞有限元分析還存在幾個問題：1）目前有限元模型參數依據鉆爆法施工的經驗值，缺乏試驗性數據及掘進機施工監測數據；2）對于過河及淺埋地段的Ⅴ類圍巖采用無粘結預應力襯砌有待進一步進行有限元計算；3）處于高埋深且節理發育段Ⅴ類圍巖的處理方式要進一步采用有限元方法模擬。

表2 各類圍巖計算結果

［1］郝文化.ANSYS土木工程應用實例［M］.北京：中國水利水電出版社，2005.

［2］薛守義.高等土力學［M］.北京：中國建材工業出版社，2007.

Application of FEM in water diversion project from Songhua river to middle cities of Jilin province

WANG Qian,HOU Wei,JU Ye

The FEM has advantages of considering heterogeneity,anisotropy,discontinuity,material nonlinearity and geometric nonlinearity of rock and soil media,adapts to various actual boundary conditions,and is superior to routine method in tunnel analysis.According to the pressure tunnel of water diversion project from Songhua river to middle cities of Jilin province,the paper calculates displacement deformation,stress and strain value caused by tunnel excavation for different class of surrounding rock with FEM,analyzes the surrounding rock stability,and provides creditable data for design,excavation and maintenance of tunnel.

finite element;pressure tunnel;surrounding rock stability;water diversion project from Songhua river

TV67

B

1002－0624（2012）06－0005－03

2011-09-22