姬棟宇

（湖南城建職業(yè)技術(shù)學(xué)院，湖南 湘潭 411101）

大型涵洞式薄殼渡槽具有受力性能好［1-2］、施工方便、溫度應(yīng)力低［3-4］等優(yōu)點，成為近年來渡槽理論和應(yīng)用研究的熱點問題。 由于大型涵洞式薄殼屬于典型的薄壁結(jié)構(gòu)，在結(jié)構(gòu)的設(shè)計和施工過程中，對其應(yīng)力應(yīng)變特性的研究有利于指導(dǎo)涵洞式薄殼渡槽設(shè)計和施工運營， 筆者針對沙河涵洞式渡槽結(jié)構(gòu)體型大、受力復(fù)雜等問題，采用有限單元法對沙河涵洞式渡槽結(jié)構(gòu)進(jìn)行了仿真分析。

1 工程概況

沙河渡槽是南水北調(diào)中線一期工程總干渠上的大型河渠交叉建筑物，項目位于河南魯山縣城東，沙河渡槽是一個組合建筑物， 跨越河道及低洼地帶采用梁式渡槽、跨越灘地采用涵洞式渡槽。沙河渡槽全長為7590m，渡槽由4部分組成，第一部分為沙河板梁式渡槽、第二部分為涵洞式渡槽（第一段）、第三部分為大郎河板梁式渡槽、第四部分為涵洞式渡槽（第二段）。 第一段渡槽的典型斷面為涵洞，有兩個涵洞并排設(shè)置雙渡槽輸水。 其凈高5700mm，凈寬5800mm，底板厚900mm， 頂板厚850mm。 涵洞式渡槽凈寬11800mm，凈高8000mm，底板厚850mm，側(cè)墻從底部厚1100mm，過渡到頂部厚400mm。

2 計算模型

2.1 模型參數(shù)

槽身采用C40混凝土，E1=32.5GPa； 涵底墊層混凝土為C15，彈性模量E2=17.5GPa；墊層下基礎(chǔ)為4m深的壤土， 彈性模量E3=11MPa； 壤土下為2.5m天然砂，彈性模量E4=25MPa；最下部為天然卵石層，彈性模量E5=40MPa。

2.2 模型單元

沙河涵洞式渡槽結(jié)構(gòu)模型采用8節(jié)點等參塊體單元計算，參塊體單元具有8個節(jié)點，每個節(jié)點具有3個平動自由度。 沙河涵洞式渡槽結(jié)構(gòu)的三維有限元計算模型如圖1。 該參塊體單元具有塑性、 蠕變、膨脹、應(yīng)力硬化、大變形、大應(yīng)變等功能。

圖1 渡槽結(jié)構(gòu)有限元計算模型

2.3 溫度荷載

在溫度荷載分析過程中主要有如下步驟： ①通過實測得出構(gòu)件表面溫度分布情況， 通過溫度升降的工況分析得出相應(yīng)的溫度應(yīng)力值； ②溫度分析塊體單元呈3個方向進(jìn)行熱傳導(dǎo)， 對8節(jié)點單元進(jìn)行三維熱分析，實現(xiàn)勻速熱流的傳遞。

2.4 混凝土徐變和應(yīng)力松弛

涵洞式渡槽在運營期， 構(gòu)件的變形隨著荷載的增大，當(dāng)構(gòu)件應(yīng)力未超過材料的彈性極限時，首先構(gòu)件將會產(chǎn)生彈性變形， 當(dāng)構(gòu)件應(yīng)力超過材料的彈性極限后，構(gòu)件發(fā)生塑性變形。同時由于荷載的持續(xù)作用，構(gòu)件在應(yīng)力不變狀態(tài)下會發(fā)生徐變。

在計算溫度應(yīng)力過程中，一般會進(jìn)行假定，即假定徐變的應(yīng)力和應(yīng)變成正比，可以得到時間t的總應(yīng)變?yōu)椋?/p>

在施工工況溫度應(yīng)力計算時，松弛系數(shù)Kp=0.3。

2.5 計算工況

考慮到渡槽結(jié)構(gòu)在運行過程中的受力特點，主要考慮了以下10種計算工況［5］：其中水槽設(shè)計深度用S表示，涵洞設(shè)計水深用HD表示，滿槽水深用M表示，渡槽自重用Z表示，河道設(shè)計水深揚壓力用Y表示，風(fēng)荷載用F表示，溫升應(yīng)力值用WS表示，溫降應(yīng)力值用WJ表示，人活動荷載用R表示，檢修荷載用JX表示，涵洞設(shè)計深度用H表示，渡槽內(nèi)水荷載用NH表示：

分析工況1考慮了：Z+S （槽內(nèi)水深6.79m）+HD（洞內(nèi)水深2m）+Y+F+WS+R；

分析工況2考慮了：Z+M（槽內(nèi)水深8m）+HD+Y+F+WS+R；

分析工況3考慮了：Z+S+F+WS+R；

分析工況4考慮了：Z+S+F+WJ+R；

分析工況5考慮了：Z+M+F+WS+R；

分析工況6考慮了：Z+M+F+WJ+R；

分析工況7考慮了：Z+F+S（一槽設(shè)計水深，一槽無水）+WS+JX；

分析工況8考慮了：Z+F+S（一槽設(shè)計水深，一槽無水）+WJ+JX；

分析工況9考慮了：Z+WS；

分析工況10考慮了：Z+WJ。

3 渡槽結(jié)構(gòu)分析

3.1 分析路徑

對沙河涵洞式渡槽結(jié)構(gòu)進(jìn)行有限元分析時，定義如圖2的分析路徑，下面分別給出路徑應(yīng)力分布情況。

圖2 渡槽結(jié)構(gòu)分析路徑示意圖

3.2 應(yīng)力分析

在如圖2的分析路徑上，各工況下渡槽和涵洞的縱向應(yīng)力變化曲線如圖3～圖7。

從圖3～圖4可知， 涵洞式渡槽底板在縱向呈現(xiàn)“連續(xù)梁” 的性質(zhì)主要表現(xiàn)： 從路徑2的分析結(jié)果可見，在底板上側(cè)跨中位置處受壓，“支座”位置處受拉；從路徑4的分析結(jié)果可見，在底板下側(cè)側(cè)跨中位置處受拉，“支座”位置處受壓；通過分析可知，分析工況5下出現(xiàn)拉滓拉max＝1.227MPa， 分析工況2下出現(xiàn)滓壓max＝0.334MPa?？梢姾词蕉刹鄣装逦恢蒙媳砻婵缰谐霈F(xiàn)壓區(qū)，上表面支座出現(xiàn)出現(xiàn)拉區(qū)，下表面跨中出現(xiàn)拉區(qū)，下表面支座出現(xiàn)壓區(qū)，此分析結(jié)果可為設(shè)計參考。

圖3 路徑2上渡槽縱向（Y向）應(yīng)力變化曲線（MPa）

圖4 路徑4上渡槽縱向（Y向）應(yīng)力變化曲線（MPa）

圖5 路徑7上涵洞縱向（X向）應(yīng)力變化曲線（MPa）

圖6 路徑10上涵洞縱向（X向）應(yīng)力變化曲線（MPa）

圖7 路徑13上涵洞縱向（X向）應(yīng)力變化曲線（MPa）

從圖5～圖7可知，受渡槽側(cè)墻的影響，涵洞頂部的縱向在每個渡槽的中部受拉，端部受壓，最大拉應(yīng)力為1.206MPa（工況6），最大壓應(yīng)力為0.839MPa（工況2）。 涵洞底部上表面沿縱向在所有工況下基本受壓，但在涵洞的內(nèi)端部形成拉應(yīng)力集中區(qū)。在涵底部下表面沿縱向在所有工況下受較小的拉應(yīng)力， 但在涵洞的內(nèi)端部形成壓力應(yīng)力集中區(qū)。

3.3 變形分析

考慮了S，Z，Y，F(xiàn)，WS，WJ，HD，NH等單荷載因素對涵洞式渡槽的沉降量、 沉降量差及槽頂橫向相對位移的影響情況。 通過計算結(jié)果分析表明，Z和NH是影響沉降量的主要因素，WS不是運營期的結(jié)構(gòu)沉降的影響因素。 Z和NH是影響涵洞式渡槽頂部橫向相對位移的主要因素， 且并排相鄰的兩個槽頂有趨于靠近的位移變化；在WS和WJ作用下，溫升過程中渡槽兩側(cè)墻位移使其有趨于分開趨勢， 溫降過程中兩槽頂呈靠近有趨于靠近的位移變化。 在分析工況6下， 涵洞式渡槽槽頂橫向相對位移達(dá)為S*max=3.445cm；在分析工況2下，沉降量Smax=-7.417cm，沉降量差為ΔSmax=1.003cm； 經(jīng)過計算分析得出工況2、工況5、工況6為應(yīng)變最大的3種工況，其最具代表性工況2的應(yīng)用云圖如圖8～圖9。

圖8 工況2下渡槽整體豎向位移云圖（m）

圖9 工況2下渡槽槽身橫向位移云圖（m）

4 結(jié)語

（1）工況5在運營期間環(huán)向拉應(yīng)力最大，因此在對涵洞式渡槽結(jié)構(gòu)設(shè)計時，可把工況5作為控制工況。

（2）根據(jù)涵洞式渡槽結(jié)構(gòu)應(yīng)力分布特點，在涵洞式渡槽底板縱向配筋時， 應(yīng)使受力鋼筋符合連續(xù)梁的配筋原則，即在跨中鋼筋靠下側(cè)布置，在支座處鋼筋靠上側(cè)布置。

（3）從應(yīng)力分析結(jié)果可以得出，涵洞式渡槽上的拉應(yīng)力和壓應(yīng)力較小，通過配筋加強，能夠滿足抗裂設(shè)計的要求。

（4）在各種單荷載計算過程因素中，對涵洞式渡槽的沉降量和槽頂橫向位移量影響較大的因素主要是涵洞式結(jié)構(gòu)自重和渡槽輸送水時的荷載。