吳永飛
(沈陽市蘇家屯區(qū)農(nóng)業(yè)技術(shù)推廣與行政執(zhí)法中心,沈陽 110101)
重力壩因其自身特點在壩工設(shè)計中廣泛應(yīng)用,高重力壩壩體斷面較大,對基巖適用性要求低于同高度拱壩。當(dāng)前壩工界重力壩設(shè)計要點為盡可能減少揚壓力,通過不斷改進(jìn)重力壩結(jié)構(gòu)形式,改良混凝土工藝和施工工序,以便能夠加快施工進(jìn)度,降低工程造價等[1],進(jìn)而針對性降低荷載帶來的不利影響,減少重力壩應(yīng)用的技術(shù)壁壘。為重力壩提供更為廣泛的應(yīng)用空間。
有限元分析中通常采用8結(jié)點六面體單元,示意圖見圖1。通過坐標(biāo)變換可得曲面空間的8結(jié)點等參元,坐標(biāo)變換關(guān)系式為:
(1)
圖1 空間8結(jié)點等參元
則單元的位移函數(shù)為:
(2)
式中:ui、vi、wi為結(jié)點位移。
由虛功原理得單元剛度矩陣為:
(3)
(4)
作用于單元結(jié)點上的等效節(jié)點力為:
{F}e=[K]e{δ}e
(5)
外力作用于單元結(jié)點上的等效荷載為:
1)若外力為體積力:
(6)
2)若外力為分布面力:
(7)
整個結(jié)構(gòu)的平衡方程為:
[K]{δ}={P}
(8)
某水電站為西北省份的梯級電站,樞紐建筑物主要包括:混凝土重力壩、泄洪底孔、電站廠房等,主要建筑物設(shè)計為3級,臨時建筑物5級,安全級別II級。最大壩高78.0m,壩頂全長174.0m。水庫正常蓄水位1895.00m,總庫容0.63億m3[3],具有日調(diào)節(jié)性能。
大壩建基面為弱風(fēng)化花崗巖體,花崗巖為中期侵入巖漿巖,巖質(zhì)堅硬,抗壓強(qiáng)度高,抗?jié)B性強(qiáng),彈性模量大?;炷两Y(jié)構(gòu)體采用SOLID45單元模擬[4],共模擬單元總數(shù)94722,節(jié)點36943個。有限元模型模擬結(jié)果見圖2。
圖2 混凝土重力壩基巖
基巖花崗巖彈性模量E=1.84×104MPa;混凝土抗剪強(qiáng)度f'=0.85,C'=0.8MPa;泊松比=0.21[5]。
混凝土容重γc=25kN/m3;泊松比c=1/6;
泥沙干容重γs=11.2 kN/m3;孔隙率0.53;內(nèi)摩擦角20°。
荷載組合包括:施工在建:自重荷載。
正常運行:自重+外水壓力+內(nèi)水壓力+泥沙壓力+揚壓力+側(cè)向外水壓力。
圖3 工況一基巖X方向應(yīng)力等值線圖
圖4 工況二基巖Y方向應(yīng)力等值線圖
圖5 工況二基巖Z方向應(yīng)力等值線圖
由圖示結(jié)果可得,基巖橫向和豎向應(yīng)力變化值較大,軸向應(yīng)力值較小。工況一順河向最大正應(yīng)力1.026MPa,出現(xiàn)在大壩基巖面接觸部位,豎直方向最大正應(yīng)力出現(xiàn)在壩基偏上10m位置。
基巖應(yīng)力變化結(jié)果見表1,位移變化結(jié)果見表2。
表1 基巖應(yīng)力變化結(jié)果
表2 基巖位移變化結(jié)果
從表1和表2中可得,基巖各向變形均不大,最大值為0.798mm,位于豎向向下變形。工況一的基巖最大拉應(yīng)力值發(fā)生在順河向,各向變形均較小,工況1、2的豎向和軸向最大拉應(yīng)力相差不多。
文章在有限元分析原理基礎(chǔ)上,結(jié)合實際重力壩工程建立壩體基巖三維有限元模型,計算不同荷載組合的壩體基巖應(yīng)力變形,得出基巖豎向應(yīng)變最大,分析其為受到大壩自重和水壓力荷載作用下的正常形變,成果符合混凝土壩設(shè)計規(guī)范規(guī)定,同時驗證了該重力壩建基面是安全合理的。