高土石壩地震反應(yīng)特性的動(dòng)力彈塑性分析

高土石壩地震反應(yīng)特性的動(dòng)力彈塑性分析

0國(guó)內(nèi)在高土石壩抗震研究中取得的進(jìn)展中國(guó)擁有豐富的水利水電資源。水利水電工程的規(guī)模和數(shù)量在世界上是第一位的，其中大部分是土石水庫(kù)。在歷次地震中,國(guó)外曾發(fā)生土石壩因地震而潰壩的嚴(yán)重災(zāi)害,在我國(guó)也有土石壩因地震而發(fā)生滑坡、震陷、裂縫等災(zāi)害。近年來,隨著西部大開發(fā)的實(shí)施,高壩大庫(kù)的建設(shè)也越來越多,其中高土石壩占了很大比例,同時(shí)這些地區(qū)的河谷條件、地質(zhì)條件和環(huán)境條件極度復(fù)雜,且一般位于高烈度地震活動(dòng)區(qū)。因而,高土石壩抗震研究工作的迫切性和重要性越來越突出。近十幾年來,國(guó)內(nèi)在高土石壩抗震研究方面有了新的進(jìn)展,取得了一些成果。本文利用FLAC(二維)商用軟件,對(duì)強(qiáng)震區(qū)高土石壩的地震反應(yīng)特性進(jìn)行系統(tǒng)深入的研究,著重討論了壩體的剪應(yīng)變、永久位移以及加速度的變化情況,從而能更好對(duì)地震作用下高土石壩的抗震安全進(jìn)行準(zhǔn)確評(píng)價(jià),為合理的抗震措施提供理論基礎(chǔ)。1拉格朗日有限差的基本原則1.1密度uuit的計(jì)算ρ?u˙i?t=?σij?xj+ρgi(1)ρ?u˙i?t=?σij?xj+ρgi(1)其中,σij為應(yīng)力張量;ρ為連續(xù)體密度;xj為j向坐標(biāo);gi為重力加速度(體力)。1.2節(jié)點(diǎn)運(yùn)動(dòng)速度的估計(jì)連續(xù)體變形時(shí)內(nèi)部各點(diǎn)要滿足以下位移協(xié)調(diào)關(guān)系:e˙ij=12[?u˙i?x˙j+?u˙j?x˙i](2)e˙ij=12[?u˙i?x˙j+?u˙j?x˙i](2)其中,e˙e˙ij為應(yīng)變率張量;u˙u˙i為速度分量。對(duì)定義在任意閉合域或單元上的連續(xù)函數(shù)F,由高斯定理得:∫sFnids=∫A?F?xidA(3)∫sFnids=∫A?F?xidA(3)其中,A為求解域或單元的面積;s為A的邊界;ni為s的單位外法線矢量。定義梯度?F?xi?F?xi的平均值為:<?F?xi>=1A∫A?F?xidA(4)<?F?xi>=1A∫A?F?xidA(4)對(duì)于有N條邊的多邊形,(4)式可以表示成對(duì)N條邊求和的形式,即:<?F?xi>=1A∑NFˉˉˉiniΔsi(5)<?F?xi>=1A∑ΝFˉiniΔsi(5)其中,FˉˉˉFˉi為第i條邊上的平均值,假設(shè)F在第i條邊上線性分布,FˉˉˉFˉi還可以直接用第i條邊兩端點(diǎn)函數(shù)值的平均值代替。對(duì)幾何關(guān)系(2)式應(yīng)用高斯定理,可以由單元節(jié)點(diǎn)的速度求得單元的平均應(yīng)變?cè)隽?<Δeij>=12[<?u˙i?xj>+<?u˙j?xi>]Δt(6)<Δeij>=12[<?u˙i?xj>+<?u˙j?xi>]Δt(6)應(yīng)用高斯定理運(yùn)動(dòng)方程(1)式可以表示成:ρu¨i=1A∑<σij>njΔsi+ρgi(7)ρu¨i=1A∑<σij>njΔsi+ρgi(7)其中,∑<σij>njΔsi為節(jié)點(diǎn)周圍單元作用在該節(jié)點(diǎn)上的集中力。(7)式右邊即為節(jié)點(diǎn)的不平衡力。利用中心差分格式,有:u¨i(t)=u˙(t+Δt/2)?u˙(t?Δt/2)Δt(8)u¨i(t)=u˙(t+Δt/2)-u˙(t-Δt/2)Δt(8)可以從前一時(shí)間步的節(jié)點(diǎn)速度u˙(t?Δt/2)u˙(t-Δt/2),求得當(dāng)前時(shí)間步的節(jié)點(diǎn)速度u˙(t+Δt/2)u˙(t+Δt/2)。這樣,(6)式可以由節(jié)點(diǎn)速度梯度得到單元的平均應(yīng)變?cè)隽?(7)式從單元應(yīng)力可以得到節(jié)點(diǎn)的運(yùn)動(dòng)速度,而單元的應(yīng)力應(yīng)變關(guān)系由本構(gòu)方程定義。FLAC程序分析問題的流程,如圖1所示。2計(jì)算模型的構(gòu)建2.1上下游壩坡計(jì)算采用上、下游完全對(duì)稱的260m高心墻堆石壩,壩底長(zhǎng)952m,上下游壩坡對(duì)稱為1∶1.8,壩頂長(zhǎng)度方向取單位寬度計(jì)算。采用位移邊界條件,模型底面限制豎向位移,側(cè)面均被約束。計(jì)算中沒有考慮上游水荷載和壩體內(nèi)水的滲流作用,按總應(yīng)力計(jì)算。FLAC網(wǎng)格如圖2所示。2.2材料剪切模量應(yīng)力的彈塑性模型分析計(jì)算中堆石體采用常用的Mohr-Coulomb(以下簡(jiǎn)稱M-C)彈塑性模型,對(duì)分析過程中材料剪切模量的變化做了2點(diǎn)考慮:①初始切線模量按照(9)式隨平均應(yīng)力變化,見圖3;②計(jì)算過程中,單元的切線模量還隨單元剪應(yīng)變按照?qǐng)D4模量衰減曲線變化。(1)an9G=CPa(σ0′Pa)n(9)G=CΡa(σ0′Ρa)n(9)其中,G、σ0′、Pa采用同一量綱;Pa為大氣壓;σ0′為平均應(yīng)力;C為系數(shù),n為指數(shù),其值可由試驗(yàn)結(jié)果確定。(2)其他參數(shù)體積彈性模量K=13G/6;粘聚力c、摩擦角φ、拉伸強(qiáng)度σt按表1取值。2.3反映加速度問題模型計(jì)算以人工地震波為主,人工地震波選用擬合糯扎渡壩址區(qū)100a超越概率為2%的基巖期望反應(yīng)譜曲線得到的人造波,峰值加速度為0.288g(2.83m/s2),見圖5所示,本計(jì)算只考慮水平向加速度輸入。3結(jié)論分析3.1高堆石壩局部塑性變形地震作用下,雖然高堆石壩宏觀破壞性態(tài)表現(xiàn)為壩頂1/5范圍內(nèi)沿壩坡淺層的滑落、坍塌,但這一破壞形式往往是由于壩體內(nèi)部累積了較大的塑性變形引起的,因此分析壩體局部塑性變形分布很有意義。采用M-C彈塑性模型和拉格朗日更新坐標(biāo)系,對(duì)高堆石壩進(jìn)行了大變形分析,如圖6所示。由圖6可見,大壩的剪應(yīng)變以壩頂處以及壩坡上部為較大,最大剪應(yīng)變發(fā)生在壩頂處,達(dá)到12%。3.2部分壩體地震長(zhǎng)期位移以往振動(dòng)臺(tái)試驗(yàn)表明,地震作用下堆石壩主要破壞性態(tài)是:壩體遭遇強(qiáng)烈振動(dòng)后,靠近壩頂?shù)膲纹卤砻鎵K石最先出現(xiàn)松動(dòng),隨著地震波作用時(shí)間延長(zhǎng),壩坡淺層塊石會(huì)發(fā)生滑動(dòng)、最后導(dǎo)致壩頂坍塌。從連續(xù)體的角度考慮,這主要是由于壩體發(fā)生了較大的永久變形引起的。因此,高土石壩動(dòng)力計(jì)算時(shí)應(yīng)重視對(duì)永久變形的分析,特別是重視對(duì)靠近壩頂局部區(qū)域的地震永久變形分析。為了進(jìn)行定量比較,圖7示意了上游、中軸線、心墻下游和下游4條特征線,圖8分別給出了它們的地震永久位移,包括水平向和豎向永久位移。從圖8可以看出:(1)各特征線的水平永久位移曲線沿高程規(guī)律分布,壩體上部1/5范圍內(nèi),水平永久位移變化較大,最大值發(fā)生在下游壩坡,大小約為2.86m;(2)各特征線的豎向永久位移曲線沿高程逐漸遞增,最大值發(fā)生在下游壩坡,大小約為1.93m;(3)對(duì)比圖8a和圖8b可知,強(qiáng)震作用下,壩體的最大水平和豎向永久位移均發(fā)生在壩頂及附近,所以應(yīng)加強(qiáng)壩頂區(qū)堆石體的穩(wěn)定,特別是下游壩坡的變形較大。3.3加速度變化情況地震作用下壩體最大加速度分布是工程設(shè)計(jì)十分關(guān)心的問題,圖9給出沿壩中軸線的加速度放大倍數(shù)。從圖9可以看出,在壩高4/5以下,壩中軸線的加速度放大倍數(shù)沿高程方向成遞增分布,最大值以及變化不規(guī)律出現(xiàn)在壩體上部1/5范圍內(nèi),最大放大倍數(shù)達(dá)到3.9,同樣表明壩體頂部在強(qiáng)震作用下的不穩(wěn)定性,應(yīng)采取相應(yīng)的抗震措施緩解加速度變化較大。3.4實(shí)際變形圖的放借助地震作用后壩體的網(wǎng)格變形圖可以更直觀地考察壩體的永久位移變化情況,為了方便觀察,將實(shí)際的變形圖放大了10倍,如圖10所示。從圖10可以看出,大壩在強(qiáng)震作用下,靠近頂部壩坡表面的單元變形較大,坍塌明顯,并且向下游傾斜,應(yīng)采取相應(yīng)的抗震措施抑制頂部壩坡表面單元的坍塌,下游壩坡可適當(dāng)減緩壩坡或加面板。4壩體極限震化分析通過用FLAC軟件對(duì)強(qiáng)震區(qū)高土石壩的地震反應(yīng)特性的計(jì)算分析,得出以下結(jié)論:(1)大壩的剪應(yīng)變以壩頂處以及壩坡上部為較大,最大剪應(yīng)變發(fā)生在壩頂處。(2)壩體的最大水平和豎向永久位移均發(fā)生在壩頂及附近。(3)大壩在強(qiáng)震作用下,靠近頂部壩坡表面的單元變形較大,坍塌明顯,并且向下游傾斜。(4)在壩高4/5以下,壩中軸線的加速度放大倍數(shù)沿高程方向成遞增分