□蘇 丹 梁春雨 戴菊英

（黃河勘測規(guī)劃設計有限公司）

基于Ansys平臺對水工隧洞復雜斷面襯砌結構計算方法

□蘇 丹 梁春雨 戴菊英

（黃河勘測規(guī)劃設計有限公司）

水工建筑物如隧洞、廠房、水池、水閘及泵房等均要考慮地基與基礎間相互作用，線性彈性地基模型較好地反映了地基特性又能便于分析不同條件下地基與基礎間相互作用，ANSYS為實現(xiàn)這種數值計算提供了一個優(yōu)秀的平臺，通過簡單條件下地基梁數值計算，得到了和理論解一致的結果。文章結合實例舉列說明在實際工程中的應用情況，表明BEAM54單元和BEAM3+COMBIN14單元較好地解決了這一類問題。

地基梁；地基與基礎；ANSYS；BEAM54COMBIN14

0 引言

隧洞支護結構計算模型主要有：載荷—結構模型與圍巖—結構模型。前者是傳統(tǒng)的結構力學模型，它將襯砌結構和圍巖分開來考慮，支護結構是承載主體，圍巖作為載荷的來源和支護的彈性支承；后者是將襯砌結構和圍巖視為一體，作為共同承載的隧洞結構體系，由于在確定圍巖初始應力場，以及表示材料非線性的各種參數確定方面比較復雜，限制了它的應用。目前人們對載荷—結構模型計算方法進行了諸多探討，尚存在如下問題：多數計算針對定制斷面進行的，通用性差；對于多曲線斷面法向、水平或豎直向荷載施加困難。

ANSYS是工程領域中應用最為廣泛的商用軟件，它通用性好，兼容性強，參數化設計方便。文章利用ANSYS平臺，對復雜斷面和荷載型式的支護結構計算方法進行了探討。

1 局部彈性地基理論

1867年前后，溫克爾（E.Winkler）對地基提出如下假設：地基表面任一點的沉降與該點單位面積上所受的壓力成正比。公式如下：

式中：y—地基的沉陷（m）；k—地基系數（kPalm）。其物理意義為：使地基產生單位沉陷所需的壓強；p為單位面積上的壓力強度。這個假設實際上是把地基模擬為剛性支座上一系列獨立的彈簧。由于計算簡便，基本能夠反映地基變形情況，在實際工程中得到了廣泛的應用。在實際運用中對彈性地基梁還做了如下假設：一是梁各點與地基緊密相貼；二是不考慮梁與基礎之間的摩擦力；三是地基梁符合平面假定。

2 基礎梁模擬

ANSYS平臺下，模擬彈性地基梁有兩種方法：一是利用BEAM54單元，該單元是二維梁單元中系列中的一個，引入地基系數后可以用來分析彈性地基梁；在進行支護結構分析時，BEAM54單元采用符合E.Winkler假定的局部彈性地基理論，BEAM54單元在模擬彈性地基梁時，需要定義常態(tài)實常數外，還需要定義彈性地基系數EFS。需要注意的是，彈性地基系數單位是N/m2，不同于地質參數中基床系數單位N/m3。二是利用彈簧單元COMBIN14和BEAM3梁單元。支護結構與地基之間相互作用用COMBIN14來模擬，支護結構用BEAM3來模擬。

一般的，地下支護結構的地基梁可以是平放，也可以是豎放，地基介質可以是巖石、粘土等固體材料。對于地基梁在加載后，會出現(xiàn)脫空現(xiàn)象，這時需要對所有梁段進行篩選，再次設定彈性地基梁，按照預先設定收斂條件，進行多次迭代求解。一般的收斂條件可以各截面前后兩次內力比值來控制，小于一個較小的數，迭代結束。

3 借助APDL參數化語言編程

APDL是ANSYS的參數化設計語言，它提供一般程序語言的功能，利用它可以實現(xiàn)參數化建模、施加參數化荷載與求解以及參數化后處理結果顯示，從而實現(xiàn)參數化有限元分析的全過程。程序流程如圖1所示。

4 理論解與數值解比較

為了驗證方法的準確性，對普通基礎梁的理論解與數值解做一對比驗算。

梁的幾何尺寸：長20 m,寬1 m，高0.50 m；單位巖石抗力系數為40 MPa/m,混凝土彈性模量2.80 GPa,泊松比0.30，跨中受集中力P=10 kN,力偶M=15 kN.m。

驗算結果如表1所示。表1表示，理論解與數值解比值在98%以上。表明方法1與方法2均能模擬彈性地基梁問題。兩種方法所求內力見圖2～圖7，表明兩種方法計算結果一致，由于方法二是用彈性簧來模擬地基與基礎相互作用，所以要想達到方法一的精度，需要更高的網格密度。

圖1 支護結構內力計算程序框圖

表1 普通彈性地基梁理論解與數值解對比表

圖2 方法一地基梁彎矩圖

圖3 方法二地基梁彎矩圖

圖4 方法一地基梁變形圖

圖5 方法二地基梁變形圖

圖6 方法一地基梁剪力圖

圖7 方法二地基梁剪力圖

5 工程應用

某一引水隧洞已知隧洞埋深200 m，地下水深50 m，混凝土標號為C25，單位巖石抗力系數1.20 GPa。該斷面為馬蹄形斷面，頂拱、側墻及底拱均為圓弧形，側墻與底拱間用圓弧連接，對于這種復雜斷面形式，由于目前常用的程序均為定制的，這種斷面無法求解，利用作者編制的程序可以解決這個問題。

分析：由于斷面為多弧段組合而成，需要在前處理時借助局部柱坐標系來設置彈簧，后處理時借助局部坐標系判別是否迭代；對于復雜荷載的施加，可借用單元遍歷自動施加，提高工作效率。采用兩種方法計算運行期工況，荷載組合：襯砌自重、內外水壓力、彈性抗力。一是兩種方法迭代最終變形圖顯示頂拱與左右下側墻三處弧段侵入圍巖，其它位置均與圍巖脫離，這種現(xiàn)象表明，在支護結構計算分析時，要特別注意區(qū)分基礎梁與一般梁，否則結果是不正確的，根據諸多計算分析，這種不正確的結果是不保守的；二是兩種方法計算支護結構的內力結果無論數值或者方向都是一致的，再次驗證了程序的可靠性，在其它工程可推廣使用；三是該工況彎矩和軸力表明，底拱與側墻下部小圓弧三處截面是控制截面，對傳統(tǒng)的馬蹄形斷面在側墻與底拱用弧段聯(lián)接對改善結構受力是有益處的。

6 結語

基于通用的有限元ANSYS平臺，模擬彈性地基梁既方便又準確，結果通用性好，能夠解決工程設計中遇到的復雜結構形狀、復雜圍巖彈抗（不同位置不同彈抗）、復雜荷載類別等問題，并且具有強大的前后處理功能，極大提高工作效率。

［1］徐志英.巖石力學［M］.北京：中國水利水電出版社，2011。

S157.2

A

1673-8853（2017）10-0068-03

2017-8-11

編輯：劉長垠