基于ANSYS的尾調(diào)室襯砌結(jié)構(gòu)有限元分析

陸周祺

（上?？睖y(cè)設(shè)計(jì)研究院有限公司，上海200434）

0 引言

尾調(diào)室是水電站引水系統(tǒng)的重要組成部分，它可以減小尾水管道中的水擊壓強(qiáng)，從而改善機(jī)組的運(yùn)行條件。尾調(diào)室襯砌的變形和破壞是影響其安全的重要因素，傳統(tǒng)的結(jié)構(gòu)分析方法存在較多的簡(jiǎn)化和假定，不能考慮結(jié)構(gòu)間的相互作用，計(jì)算成果的精度不能得到保證。而采用有限元分析方法，可以對(duì)尾調(diào)室襯砌結(jié)構(gòu)進(jìn)行三維建模計(jì)算，采用適當(dāng)?shù)倪吔鐥l件，可以較為精確地得到各個(gè)部位的受力情況[1-3]。因此，本文采用有限元分析軟件ANSYS，對(duì)尾調(diào)室襯砌結(jié)構(gòu)進(jìn)行有限元分析計(jì)算，而計(jì)算得到的內(nèi)力也為襯砌配筋設(shè)計(jì)提供了依據(jù)。

厄瓜多爾某水電站總裝機(jī)253 MW，主要建筑物包括取水口、廠房、大壩等。水電站尾調(diào)室與廠房平行布置，尾調(diào)室斷面為城門洞形，從頂部往底部開挖，其頂部的通風(fēng)洞兼施工支洞，邊開挖邊支護(hù)，支護(hù)方案采用掛網(wǎng)噴混凝土+全長黏結(jié)型砂漿錨桿+預(yù)應(yīng)力錨索+預(yù)應(yīng)力錨桿。尾調(diào)室采用鋼筋混凝土襯砌，襯砌厚0.4 m，襯砌后凈寬10.0 m，凈長41.3 m，拱頂至底面凈高36.5 m，錨桿直徑32.0 mm，長8.0 m，間距1.5 m。

1 計(jì)算模型

利用通用有限元程序ANSYS 對(duì)尾調(diào)室襯砌結(jié)構(gòu)進(jìn)行線彈性分析，襯砌斷面尺寸為41.3 m×10.0 m×36.5 m（長×寬×高），厚度為0.4 m，襯砌采用殼單元（Shell281）進(jìn)行模擬[4，5]。

由于襯砌能夠抵御的外水壓力十分有限，在采用排水措施降低外水壓力的同時(shí)，將支護(hù)所用的錨桿預(yù)留適當(dāng)長度，澆筑在襯砌內(nèi)，與襯砌共同承受外部靜水壓力，加強(qiáng)其抵御外壓的能力。由于襯砌沿錨桿方向的位移被錨桿限制，所以用相應(yīng)方向的位移支座來模擬錨桿。

圖1 為根據(jù)襯砌結(jié)構(gòu)尺寸建立的三維有限元模型，圖2 為根據(jù)上述邊界條件設(shè)定的位移支座模擬錨桿。

2 計(jì)算參數(shù)及工況

2.1 計(jì)算參數(shù)

尾調(diào)室地表高程約935 m，根據(jù)鉆孔資料判斷覆蓋層厚度在14.5～37.0 m，其下為基巖，基巖主要為安山巖和巖屑凝灰?guī)r。根據(jù)鉆孔資料，地下水位埋深一般10～20 m，為覆蓋層中的孔隙性潛水?；炷烈r砌及巖石的材料參數(shù)如表1 所示。

圖1 尾調(diào)室有限元計(jì)算模型

圖2 位移支座模擬錨桿圖

表1 土石料的滲透系數(shù)

2.2 計(jì)算工況及荷載

對(duì)于襯砌的設(shè)計(jì)，參考EM 1110-2-2901 規(guī)范中規(guī)定的荷載組合和相應(yīng)的荷載系數(shù)[6]。襯砌所受的荷載組合如表2 所示。

表2 設(shè)計(jì)荷載組合及荷載系數(shù)表

對(duì)表2 中除自重外的各項(xiàng)荷載說明如下：

1）巖石荷載。由于開挖支護(hù)時(shí)間較長，受開挖擾動(dòng)的圍巖已趨于穩(wěn)定，應(yīng)力基本釋放完成。另外，襯砌厚度相對(duì)調(diào)壓室尺寸很小，在外水壓力作用下，襯砌呈現(xiàn)脫離圍巖趨勢(shì)，所以計(jì)算襯砌時(shí)，不考慮巖石荷載的影響。

2）運(yùn)行靜水壓力。指的是正常運(yùn)行時(shí)的內(nèi)水壓力，減去在正常運(yùn)行條件下的最小外水壓力。

3）瞬態(tài)靜水壓力。指的是發(fā)生水擊的最大瞬間內(nèi)水壓力減去最小外水壓力。

4）外部靜水壓力。指的是作用在放空的尾調(diào)室上的最大外部靜水壓力。

對(duì)于荷載組合U1 和U2，襯砌在內(nèi)水壓力作用下向外擠壓圍巖，尾調(diào)室圍巖類別為Ⅲ類，圍巖與襯砌共同承受內(nèi)水壓力，同時(shí)考慮與外水壓力的組合，因此不屬于襯砌配筋計(jì)算的控制工況。對(duì)比荷載組合U3 和U4，其中U3 組合的外水壓力數(shù)值更大，荷載系數(shù)也很大。

故尾調(diào)室襯砌的結(jié)構(gòu)計(jì)算以荷載組合U3 作為控制工況，因?yàn)椴豢紤]巖石壓力的影響，所以承受的荷載為襯砌的自重、作用在放空的尾調(diào)室上的最大外部靜水壓力。

用位移支座模擬錨桿時(shí)，在外水壓力作用下，每個(gè)支座所受到的力不能超過錨桿的最大抗拉強(qiáng)度，否則錨桿屈服，支座失效。尾調(diào)室錨入襯砌的錨桿直徑為32 mm，屈服強(qiáng)度設(shè)計(jì)值為420 MPa。據(jù)此計(jì)算，錨桿即每個(gè)支座能承受的最大拉力為：420×（0.016）2×3.14×103=337.6 kN，當(dāng)每個(gè)支座上所承受的力超過337.6 kN時(shí)，錨桿將出現(xiàn)破壞屈服，支座也就同時(shí)失效，最大外部靜水壓力將直接作用于襯砌。由于錨桿間距為1.5 m，相當(dāng)于每個(gè)支座承受面積1.5×1.5=2.25 m2的水壓力。根據(jù)錨桿能承受的最大拉力及所承受水壓力的面積，推算單根錨桿可承受的最大水壓為337.6/2.25/1.4=107.2 kPa，大約相當(dāng)于10 m的最大外部靜水壓力。

根據(jù)上述計(jì)算分析結(jié)果，采取一定的排水措施后，最大外部靜水壓力需要控制在不大于10 m。此次計(jì)算中，最大外部靜水壓力按10 m 考慮。

3 計(jì)算成果

根據(jù)上述計(jì)算工況及參數(shù)，利用ANSYS 對(duì)尾調(diào)室襯砌進(jìn)行有限元分析計(jì)算，得到各個(gè)部位的彎矩值，從而為襯砌配筋設(shè)計(jì)提供了依據(jù)。尾調(diào)室各個(gè)部位各個(gè)方向的彎矩極值如表3 所示。

表3 尾調(diào)室襯砌有限元內(nèi)力計(jì)算結(jié)果匯總表kN·m

底板水平方向和豎直向彎矩均為由邊墻4 向邊墻2 從最大正彎矩變?yōu)樽畲筘?fù)彎矩，與邊墻4 連接處為最大正彎矩40.2 kN·m。由于邊墻2 與底板連接處開孔，在開孔處彎矩達(dá)到最大負(fù)值76.9 kN·m，未開孔處的彎矩遠(yuǎn)小于開孔處彎矩，說明在外部靜水壓力及自重作用下，底板與邊墻2 和邊墻4 連接處彎矩?cái)?shù)值最大。而由于存在開孔，應(yīng)力集中，在邊墻2 開孔處產(chǎn)生了更大的彎矩，因而在襯砌結(jié)構(gòu)設(shè)計(jì)時(shí)，需要對(duì)開孔部位更加關(guān)注。

邊墻1彎矩由頂向下從最大正彎矩19.6 kN·m變?yōu)樽畲筘?fù)彎矩49.1 kN·m，由于邊墻1 為結(jié)構(gòu)完整，不存在開孔，整體彎矩分布較為均勻，不存在突變現(xiàn)象。

邊墻2 由于頂部存在通風(fēng)洞開孔，底部存在泄水洞開孔，墻身上下都有開孔，彎矩變化規(guī)律與邊墻1完全不同，正彎矩極值出現(xiàn)在邊墻中心處，即兩個(gè)開孔中間，最大正彎矩為16.6 kN·m，彎矩由中間向上下兩側(cè)減小變?yōu)樨?fù)彎矩，最大負(fù)彎矩出現(xiàn)在邊墻2與底板連接處，最大負(fù)彎矩為53.2 kN·m。

邊墻3下側(cè)有3個(gè)進(jìn)水洞開孔，彎矩分布規(guī)律與邊墻1類似，由頂向下從正彎矩變?yōu)樨?fù)彎矩，在兩開孔之間與底板連接處達(dá)到最大負(fù)彎矩49.6 kN·m。而最大正彎矩為24.1 kN·m，出現(xiàn)在兩開孔頂部圓弧之間，并不是邊墻3 頂部與頂板連接處，這同樣是由于開孔，兩開孔之間結(jié)構(gòu)縮小，應(yīng)力集中，產(chǎn)生了更大的彎矩，最大正彎矩和最大負(fù)彎矩均出現(xiàn)在開孔周圍。

邊墻4 彎矩分布規(guī)律與邊墻1 完全一致，與頂板連接處為最大正彎矩17.5 kN·m，與底板連接處為最大負(fù)彎矩51.8 kN·m。

頂板與邊墻1 和邊墻2 連接處存在最大負(fù)彎矩24.4 kN·m，彎矩由邊墻1 和邊墻2 兩側(cè)往中間變小。頂板與邊墻3 和邊墻4 連接處存在最大正彎矩19.3 kN·m，同樣由邊墻3 和邊墻4 兩側(cè)向中間變小，說明頂板與邊墻連接處彎矩最大，中間彎矩相對(duì)較小。相比較于底板與邊墻，頂板彎矩極值為最小，由于不同方向結(jié)構(gòu)尺寸不同，最大負(fù)彎矩及最大正彎矩值均出現(xiàn)在水平方向。

另外由于錨桿的作用，邊墻、底板和頂板與錨桿連接處的彎矩均較小，這是因?yàn)榻r(shí)考慮襯砌沿錨桿方向的位移被錨桿限制，用相應(yīng)方向的位移支座來模擬錨桿，限制相應(yīng)方向的位移。所以，襯砌與錨桿連接處彎矩均較小，由于錨桿位移支座的作用，錨桿周圍的應(yīng)力重新分布，越靠近錨桿越小，遠(yuǎn)離錨桿應(yīng)力相應(yīng)變大，根據(jù)有限元分析結(jié)果，錨桿位移支座的影響范圍約為1/4 錨桿間距，錨桿形成支點(diǎn)約束了襯砌的位移，從而減小了襯砌的跨度，提高了襯砌本身的承載力，對(duì)襯砌抵抗外水壓力有相當(dāng)大的作用。

根據(jù)頂板、邊墻、底板的彎矩分布，可以發(fā)現(xiàn)彎矩值大體上遵循以下規(guī)律：從上至下，由正彎矩極值往負(fù)彎矩極值變化，并且最大值出現(xiàn)在底板與邊墻連接處，分布規(guī)律與傳統(tǒng)結(jié)構(gòu)分析方法一致，說明此次模型有限元分析成果較為合理。根據(jù)上述計(jì)算得到的最大彎矩值76.9 kN·m，對(duì)400 m 厚襯砌結(jié)構(gòu)進(jìn)行配筋計(jì)算，經(jīng)過計(jì)算，按最小配筋率配筋就能滿足計(jì)算結(jié)果，最終采用配筋結(jié)果為φ 18@200 mm[7，8]。

4 結(jié)論

1）調(diào)壓室放空檢修時(shí)外水壓力作用于襯砌時(shí)是最不利的工況，應(yīng)采取有效的排水措施來降低外水壓力。

2）尾調(diào)室襯砌結(jié)構(gòu)跨度大而壁厚只有0.4 m，承受外水壓力的能力有限，因此將錨桿錨入襯砌中，形成支點(diǎn)約束襯砌的位移，減小了襯砌的跨度，提高了襯砌抵抗外水壓力的能力。

3）為減少錨桿受力，在調(diào)壓室放空時(shí)需要控制放空速度，增加必要的排水廊道，將外壓水頭限制在一定數(shù)值以內(nèi)。建議采取以下措施：沿尾調(diào)室布置一層排水廊道，并施打長排水孔收集巖體中的裂隙水；對(duì)尾調(diào)室周邊的巖體進(jìn)行固結(jié)灌漿，形成一個(gè)封閉的防滲圈，阻止防滲圈以外的裂隙水作用到襯砌外壁上；在襯砌上布置間距1.5 m×1.5 m 入巖1.0 m 的短排水孔，將襯砌與巖體之間的水排入尾調(diào)室內(nèi)。

4）通過排水減壓和錨桿錨入襯砌后，襯砌所受的彎矩較小，單米最大彎矩為77 kN，采用φ 18@200 的配筋即可滿足規(guī)范要求。