付清凱,劉衛(wèi)其,張 登

(1.南水北調(diào)中線干線工程建設管理局天津分局，天津 300000；2.河北工業(yè)大學土木與交通學院，天津 300401)

0 引 言

天津干線工程是南水北調(diào)中線一期工程的重要組成部分，是解決天津市水資源緊缺、保證國民經(jīng)濟可持續(xù)發(fā)展的重要戰(zhàn)略性工程。天津干線目前整體運行狀態(tài)良好，基本處于安全狀態(tài)下，但在變形縫部位出現(xiàn)了多處滲水的情況，為了解箱涵在滲漏條件下的運行狀態(tài)，針對滲漏處箱涵結構力學特性進行了研究[1-2]。

1 箱涵滲漏模擬簡介

1.1 涵滲漏情況

箱涵在運行過程中，部分地區(qū)出現(xiàn)表面出水現(xiàn)象，滲漏水流出表面，影響了箱涵周圍的土體，給箱涵安全運行帶來隱患。經(jīng)勘察分析發(fā)現(xiàn)，于箱涵變形縫處出現(xiàn)滲漏點，導致箱涵內(nèi)水在承壓條件下發(fā)生滲漏，初步排查，箱涵變形縫處滲漏點分布無明顯規(guī)律。

1.2 計算模擬方案

1.2.1 箱涵結構

天津段的三孔箱涵，單節(jié)箱涵長度15m，孔口尺寸為4.4×4.4m，箱涵材質(zhì)鋼筋混凝土，計算時采取C30混凝土參數(shù)，箱涵斷面計算模型圖如圖1。

圖1 箱涵模型斷面圖

1.2.2 模擬方案

運用FLAC3D靈活地進行計算模擬[3]，其單元材料采用線性或非線性模型，材料發(fā)生屈服后，能夠發(fā)生相應的變形或移動，在采用顯式拉格朗日算法和混合-離散分區(qū)技術[4]的情況下，模擬能夠非常準確地表現(xiàn)材料的塑性破壞和流動[5]。

為充分模擬箱涵滲漏情況，將箱涵滲漏點分布于變形縫上，其滲漏點位(頂部、側(cè)面、底部、角部)為四個，具體位置如圖2。模型土體的左右邊界及下邊界為法向零位移，零速率的不可透水邊界，土體上邊界無約束可透水的邊界，箱涵邊界為不可透水邊界。

圖2 箱涵滲漏點

結合參考地質(zhì)勘察報告中的土層參數(shù)，給出每個地區(qū)進行計算模擬時的參數(shù)如表1。若考慮回填土，參考相關的工程經(jīng)驗，可以考慮取其滲透系數(shù)為正常條件下的2-3倍。

表1 土層計算參數(shù)表

2 模擬結果

2.1 三維模擬模型

計算時在水流方向取15m長度，橫向40m計算寬度，垂直地平方向取20m，其三維模型如圖3。

圖3 三維模型圖

2.2 滲漏模擬分析

進行滲漏模擬時，以單一點位滲漏為前提，主要分析不同點位滲漏對箱涵周圍土體及箱涵結構的影響，給出圖2中A位置在多種工況下的模擬分析結果。

2.2.1 非回填土條件下A位置滲漏模擬結果及分析

該情況下的模擬分析主要考慮，箱涵中間底部開始滲漏，箱涵周圍土當做非回填土，非回填土的滲透系數(shù)相差較小，滲漏水主要在重力作用下向下滲流，直至相對穩(wěn)定狀態(tài)。

如圖4，滲漏時間50.2d，飽和區(qū)域主要分布在箱涵下部。箱涵中間部位底部產(chǎn)生裂縫，產(chǎn)生滲漏，滲漏范圍逐漸向下向周圍擴散，由于周圍土的滲透系數(shù)相差不大，滲漏水沒有明顯的滲流傾向，因此主要在重力作用下向下流動，最終達到滲流范圍擴散很慢時，計算停止。如圖5，土體浸水軟化造成土體抗剪強度及黏聚力下降，在地應力重分部及相應的土體重力作用下，軟化土體產(chǎn)生蠕動，導致箱涵上部產(chǎn)生最大豎向位移為4.27cm，位移分布區(qū)域呈現(xiàn)明顯對稱形狀。

圖4 土體飽和度分布圖

圖5 箱涵及土體豎向位移圖

如圖6，滲漏水的滲流范圍相對穩(wěn)定時，箱涵結構應力云圖，最大主應力為1.65Mpa，位于箱涵中間部位，箱涵處于安全狀態(tài)。

圖6 箱涵結構應力云圖

2.2.2 回填土條件下A位置滲漏模擬結果及分析

該情況下的模擬分析主要考慮，箱涵中間底部開始滲漏，箱涵周圍土體是回填土，回填土相應的滲透系數(shù)相對于原狀土來說更高，因此滲漏水容易在回填土區(qū)域產(chǎn)生滲流。

如圖7，為滲漏時間進行到23.5d時的土體飽和度分布圖，此時飽和區(qū)域主要分布在相應的回填土區(qū)域。箱涵中間部位底部產(chǎn)生裂縫，滲漏水向滲透系數(shù)較大的區(qū)域流動，逐漸向上發(fā)展，直至地表出水。如圖8，當箱涵底部發(fā)生滲漏后，滲漏水向箱涵周圍回填土區(qū)域逐漸發(fā)展，土體浸水軟化造成土體抗剪強度及黏聚力下降，在地應力重分部及土體重力作用下，軟化土體產(chǎn)生蠕動，最大位移位于箱涵角點外部，為4.52cm，位移區(qū)域呈對稱形狀。

圖7 土體飽和度分布圖

圖8 箱涵及土體豎向位移圖

如圖9，為滲漏水流出地表后的箱涵結構應力云圖，最大應力為1.15MPa，出現(xiàn)在箱涵中間部位，此時箱涵處于安全狀態(tài)。

圖9 箱涵結構應力云圖

2.2.3 多位置多工況條件下滲漏模擬結果

在進行箱涵模擬計算時，除考慮了如圖2的多個滲漏位置外，還考慮了土質(zhì)及底部墊層的影響，現(xiàn)將各種工況下的模擬結果匯總?cè)绫?，前文中的算例分別對應編號3、1。

表2 多工況模擬計算匯總表

3 結 語

箱涵滲漏時滲漏水對周圍土體產(chǎn)生作用，使土體強度變低，進而造成土體發(fā)生變形及沉降，對箱涵結構產(chǎn)生影響。在一定時間內(nèi)，土體產(chǎn)生位移對箱涵本身位移的影響不大，箱涵的位移相對很小，并且其應力普遍小于2MPa，處于混凝土應力限值以下，因此可以判斷，在文中模擬的情況下，箱涵處于安全狀態(tài)。在長時間滲漏情況下，對箱涵結構的影響不可忽略，建議進行補漏修補操作。