駱秀萍，易 平

(新疆兵團(tuán)勘測(cè)設(shè)計(jì)院(集團(tuán))有限責(zé)任公司，烏魯木齊 830002)

1 工程概況

烏拉泊水庫(kù)是烏魯木齊河中游的一座攔河水庫(kù)，被稱為“烏市頭上一盆水”。水庫(kù)始建于1959年，1961年建成投入運(yùn)行。1975年被列為全國(guó)37座重點(diǎn)病險(xiǎn)水庫(kù)之一，1978年又被列為全國(guó)13座重點(diǎn)抗震加固水庫(kù)之一。由于壩體防滲型式的復(fù)雜性(壩體壤土斜墻與混凝土防滲墻之間的飽和體存在超靜孔隙水壓力，有液化可能，水庫(kù)抗震設(shè)防存在隱患)，水庫(kù)的抗震安全直接影響烏魯木齊市各族人民生命財(cái)產(chǎn)安全和城市發(fā)展，因此復(fù)核水庫(kù)大壩的抗震安全性是十分必要的[1-6]。鑒于烏拉泊水庫(kù)壩體結(jié)構(gòu)的復(fù)雜性，壩體的抗震穩(wěn)定性動(dòng)力分析委托中國(guó)水利水電科學(xué)研究院巖土工程研究所完成。

烏拉泊水庫(kù)由大壩、放水涵洞、溢洪道等組成，總庫(kù)容7 000萬(wàn)m3，最大壩長(zhǎng)1 050 m，最大壩高26 m，工程設(shè)防烈度為8度，大壩按1級(jí)建筑物設(shè)計(jì)。0+028.00 m～0+092.00 m壩段為有水平鋪蓋的壤土心墻砂殼壩，0+092.00 m～0+806.80 m段為有水平鋪蓋的壤土斜墻砂殼壩，0+806.80 m～1+050.00 m段為人工填筑砂礫石均質(zhì)壩。0+028.00 m～0+806.80 m壩段，在壩體上游馬道設(shè)置有防滲墻，基礎(chǔ)深入基巖以下1 m。

2 抗震穩(wěn)定分析

抗震穩(wěn)定分析[7-13]包括2項(xiàng)內(nèi)容：① 按照《水工建筑物抗震設(shè)計(jì)規(guī)范》[14]的規(guī)定，對(duì)壩坡進(jìn)行常規(guī)的瑞典圓弧擬靜力法抗震穩(wěn)定分析；② 根據(jù)有限元靜、動(dòng)應(yīng)力計(jì)算的結(jié)果，對(duì)壩坡進(jìn)行應(yīng)力滑弧穩(wěn)定分析。本文針對(duì)烏拉泊水庫(kù)進(jìn)行了設(shè)計(jì)洪水位工況下大壩上、下游坡的抗震穩(wěn)定計(jì)算[15]。

在應(yīng)力滑弧穩(wěn)定分析中，首先，將擬定的滑弧體分成若干個(gè)垂直條塊，條塊寬度一般不大于有限元單元寬度的1/2。確定條塊底部滑弧中點(diǎn)的坐標(biāo)值，找出形心與某滑弧中點(diǎn)坐標(biāo)最近的單元，取此單元的應(yīng)力作為該滑弧中點(diǎn)的應(yīng)力值，當(dāng)滑弧通過(guò)多個(gè)單元時(shí)，以滑弧中點(diǎn)所在單元的應(yīng)力值代表滑弧的平均應(yīng)力值，以避免繁瑣計(jì)算。然后，對(duì)滑弧中點(diǎn)的應(yīng)力值進(jìn)行坐標(biāo)變換，即把由直角坐標(biāo)系表示的有效徑向正應(yīng)力、環(huán)向正應(yīng)力和切向剪應(yīng)力σx′、σy′與τxy變換為以滑弧中心為極點(diǎn)的極坐標(biāo)系表示的σθ′、σr′與τrθ。其中，σθ′、σr′分別為極坐標(biāo)系下的有效徑向正應(yīng)力和環(huán)向正應(yīng)力，τrθ為切向剪應(yīng)力。類似地，由有限元?jiǎng)討?yīng)力分析計(jì)算得出的σxd、σyd與τxyd變換出τxθd，將切向動(dòng)剪應(yīng)力τxθd乘以0.65的系數(shù)與切向剪應(yīng)力τrθ相加，即為該滑弧條塊上的滑動(dòng)力。

根據(jù)各條塊滑弧中點(diǎn)的初始徑向正應(yīng)力σr′和初始剪應(yīng)力比α=τrθ/σr′，可確定該滑弧的靜強(qiáng)度、動(dòng)力總剪強(qiáng)度為：

τf=c+σr′tanφ

(1)

τf′=τf(σr′,α=τrθ/σr′)

(2)

式中：c為土的黏聚力；φ為土的內(nèi)摩擦角。

在應(yīng)力穩(wěn)定分析中，對(duì)各條塊采用動(dòng)力總剪強(qiáng)度，當(dāng)動(dòng)力總剪強(qiáng)度高于靜強(qiáng)度值時(shí)，采用靜強(qiáng)度值?；∶嫔系恼w抗震穩(wěn)定安全系數(shù)為：

(3)

式中:ΔL為各條塊的滑弧長(zhǎng)度。

在分析中，首先在大范圍搜索安全系數(shù)最小的滑弧圓心位置和滑弧半徑，逐步縮小搜索步長(zhǎng)，找出安全系數(shù)最小的可能滑動(dòng)面。

抗震穩(wěn)定計(jì)算中使用的壩基、壩體材料的密度與靜強(qiáng)度指標(biāo)，見(jiàn)表1。在應(yīng)力穩(wěn)定分析中，采用Dr=0.72的地震總應(yīng)力抗剪強(qiáng)度。

表1 抗震穩(wěn)定靜力計(jì)算參數(shù)表

注：對(duì)水下土石料，密度為浮密度；對(duì)水上土石料，密度為天然密度。

二維穩(wěn)定計(jì)算的結(jié)果以0+150.00 m為典型斷面，計(jì)算成果見(jiàn)表2。

表2 0+150.00 m斷面壩坡最小安全系數(shù)表

根據(jù)三維有限元靜、動(dòng)力計(jì)算得出的靜、動(dòng)應(yīng)力，得出高壩段0+150.00 m、0+200.00 m斷面、S形壩段0+450.00 m、0+505.00 m以及0+600.00 m斷面的靜、動(dòng)應(yīng)力。以下再采用應(yīng)力穩(wěn)定分析法對(duì)壩體的上、下游壩坡進(jìn)行計(jì)算，得出了5個(gè)斷面應(yīng)力穩(wěn)定法的壩坡最小安全系數(shù)，詳見(jiàn)表3。

表3 應(yīng)力穩(wěn)定法壩坡最小安全系數(shù)表

3 方法評(píng)述

從表2、3與圖1可以看出：

(1) 當(dāng)基巖輸入加速度峰值為2.00 m/s2時(shí)，采用瑞典圓弧法和應(yīng)力穩(wěn)定分析，得出上、下游壩坡的安全系數(shù)均大于1.30；當(dāng)基巖輸入加速度峰值為3.18 m/s2時(shí)，采用應(yīng)力穩(wěn)定分析，得出上、下游壩坡的安全系數(shù)均大于1.20，滿足規(guī)范要求。

(2) 與傳統(tǒng)的瑞典圓弧法計(jì)算結(jié)果相比，由應(yīng)力穩(wěn)定分析得到的安全系數(shù)略高，危險(xiǎn)滑弧略淺，這是由于按照瑞典圓弧法進(jìn)行擬靜力法的穩(wěn)定計(jì)算時(shí)，在壩的上部區(qū)需施加較大的水平慣性力，導(dǎo)致壩體上部的抗震安全性較低。而應(yīng)力穩(wěn)定分析直接采用了靜、動(dòng)應(yīng)力，能更客觀地反映滑弧面上的應(yīng)力狀態(tài)。特別是應(yīng)力穩(wěn)定分析中動(dòng)應(yīng)力直接應(yīng)用了動(dòng)力反應(yīng)分析的結(jié)果，可以考慮不同地震輸入的影響；此外，在應(yīng)力穩(wěn)定分析中，還可以利用動(dòng)力試驗(yàn)得到的地震總應(yīng)力抗剪強(qiáng)度，這都是傳統(tǒng)的瑞典圓弧法所無(wú)法采用的。

圖1 抗震穩(wěn)定分析的危險(xiǎn)滑弧位置圖

綜上所述，在高烈度地區(qū)建壩采用應(yīng)力穩(wěn)定分析進(jìn)行抗震計(jì)算，對(duì)降低工程造價(jià)、保證大壩抗震安全是十分必要的。