毛有智

(貴州水投水務(wù)有限責(zé)任公司，貴州貴陽 550003)

0 引言

碾壓混凝土筑壩通過改變常態(tài)混凝土壩需要振搗密實(shí)，代之層面碾壓，可以把大型施工機(jī)械用于壩體層面［1-3］。重力壩壩體結(jié)構(gòu)較為簡(jiǎn)單，不設(shè)置縱縫、不用模板避免了橫縫的形成，相對(duì)別的壩體，碾壓混凝土重力壩更可以加快施工［4］。碾壓混凝土相比常態(tài)混凝土水泥用量少，節(jié)約造價(jià)。隨之而來，我國(guó)近幾十年碾壓混凝土重力壩越來越多，尤其是百米級(jí)高壩［5］。高壩往往更需要考慮地震的作用，一旦遭受地震災(zāi)害，那么將產(chǎn)生不可估計(jì)的后果。因此，碾壓混凝土重力壩的地震動(dòng)力特性的研究意義重大［6-8］。下文采用通用有限元軟件ANSYS分析一120 m的碾壓混凝土重力壩，說明其進(jìn)行地震動(dòng)力反應(yīng)譜分析的過程。

1 反應(yīng)譜方法簡(jiǎn)介

地震的抗震性能分析主要有三種方法:擬靜力法、時(shí)程分析法和反應(yīng)譜方法。

擬靜力法是把地震作用以折算靜荷載的方式來表示。表示后的地震荷載按常規(guī)靜力法對(duì)壩體進(jìn)行應(yīng)力、位移的抗震分析。假定地震作用時(shí)，與地面加速度相同的加速度作用于壩體，求出慣性力后對(duì)大壩安全做出評(píng)價(jià)，是一種比較傳統(tǒng)的做法。

時(shí)程分析法需將地震記錄作用在結(jié)構(gòu)上，對(duì)運(yùn)動(dòng)方程進(jìn)行時(shí)間積分，以求得任意時(shí)刻結(jié)構(gòu)的地震反應(yīng)。結(jié)構(gòu)規(guī)范［9］條文說明的5.1.2條，認(rèn)為時(shí)程分析法應(yīng)作為抗震規(guī)范的補(bǔ)充計(jì)算方法。而水利規(guī)范［10］的4.5.8條對(duì)時(shí)程分析法也進(jìn)行了相應(yīng)規(guī)定。

反應(yīng)譜方法的全稱為振型分解反應(yīng)譜法。先求結(jié)構(gòu)對(duì)應(yīng)各階振型的地震效應(yīng)后，再組合得到總的地震作用。反應(yīng)譜方法又分為單點(diǎn)一致輸入、多點(diǎn)輸入、多向單點(diǎn)輸入和多向多點(diǎn)輸入。本文介紹使用的單質(zhì)點(diǎn)一致輸入，是以單點(diǎn)在實(shí)際地震過程中的反應(yīng)作為基礎(chǔ)，進(jìn)行結(jié)構(gòu)反應(yīng)分析。用反應(yīng)譜將動(dòng)力問題靜力化，簡(jiǎn)化復(fù)雜的地震反應(yīng)。按照實(shí)際地面運(yùn)動(dòng)計(jì)算壩體的反應(yīng)，得到地震譜曲線。這條曲線是實(shí)際壩體的最大反應(yīng)和自身體系自振周期的函數(shù)關(guān)系。

2 碾壓混凝土重力壩模型建立

一碾壓混凝土重力壩如圖1所示?？拷嫌渭皦误w使用常態(tài)混凝土，材料參數(shù):楊氏彈性模量 E=28.5 GPa，泊松比 v=0.167，容重 γ =24.5 kN/m3;碾壓混凝土 1 分區(qū) rcc1:E1=26 GPa，v1=0.163，γ1=24 kN/m3;碾壓混凝土 2 分區(qū) rcc2:E2=25 GPa，v2=0.163，γ2=24 kN/m3;碾壓混凝土 3 分區(qū) rcc3:E3=24 GPa，v1=0.163，γ1=24 kN/m3。壩體上游水位為 100.0 mm，下游水位為30.0 mm。荷載只考慮上、下游靜水壓力以及揚(yáng)壓力作用，壩體底部假設(shè)全部固結(jié)于基巖上，浪壓力、泥沙壓力、動(dòng)水壓力等均不考慮。

ANSYS建模時(shí)選用Plane82單元。此單元為平面8節(jié)點(diǎn)四邊形單元，可處理不規(guī)則邊界，退化為6節(jié)點(diǎn)三角形單元。單元的關(guān)鍵項(xiàng)3需設(shè)置為strain平面應(yīng)變形式。有限元網(wǎng)格模型見圖2，圖中共有857個(gè)單元，1 824個(gè)節(jié)點(diǎn)。加載靜力求解得到X，Y方向最大/最小位移值分別為4.28 mm(受拉)/0;0/－3.45 mm(受壓)。X，Y方向最大/最小應(yīng)力值分別為 0.01 MPa(受拉)/－0.54 MPa(受壓);0/－1.77 MPa(受壓)。從應(yīng)力值看出符合規(guī)范［11］6.3.4 條對(duì)重力壩的應(yīng)力要求。

圖1 碾壓混凝土重力壩模型圖

圖2 碾壓混凝土重力壩ANSYS網(wǎng)格劃分圖

3 模態(tài)分析振型獲取

由1節(jié)可知，復(fù)雜的結(jié)構(gòu)可以轉(zhuǎn)化成若干振型的疊加，而每個(gè)振型又可看作單個(gè)質(zhì)點(diǎn)。本節(jié)就是對(duì)建立好的ANSYS靜力模型進(jìn)行模態(tài)分析，獲得重力壩的振動(dòng)特性，包括其固有的自振頻率和振型，是結(jié)構(gòu)物承受地震荷載的重要?jiǎng)討B(tài)參數(shù)，同時(shí)也是反應(yīng)譜分析的基礎(chǔ)。模態(tài)分析采用子空間迭代法提取9階模態(tài)，可得碾壓混凝土重力壩9個(gè)固有的自振頻率Freq，由小到大，分別為 3.068，7.041，9.656，12.452，18.910，21.594，26.610，29.881 和31.387。限于篇幅只給出前3階振型圖，見圖3～圖5。

4 反應(yīng)譜分析

假設(shè)該碾壓混凝土重力壩所處位置地震設(shè)防烈度為8度，則可根據(jù)《水工建筑物抗震設(shè)計(jì)規(guī)范》［10］查得地震水平加速度值為0.2g。規(guī)范中4.3.4條規(guī)定設(shè)計(jì)反應(yīng)譜的最大代表值 βmax=2.0。Ⅰ類場(chǎng)地類別，對(duì)應(yīng)的特征周期為 Tg=0.2 s。根據(jù)規(guī)范［10］4.3.3條中的圖示公式，得大壩反應(yīng)譜曲線方程為:

其中，T為振動(dòng)周期，由3節(jié)計(jì)算結(jié)果得到，T=1/Freq。

圖3 碾壓混凝土重力壩第1階振型圖

圖4 碾壓混凝土重力壩第2階振型圖

圖5 碾壓混凝土重力壩第3階振型圖

式(1)計(jì)算得 9 個(gè)對(duì)應(yīng)的反應(yīng)譜值分別為 1.289，2.0，2.0，1.803，1.529，1.463，1.376，1.335，1.319。把反應(yīng)頻率和反應(yīng)譜值一起輸入ANSYS，進(jìn)行反應(yīng)譜spectrum分析，結(jié)果再與模態(tài)擴(kuò)展modal expass分析合并，得到地震荷載下壩體反應(yīng)譜分析的主要的1階～5階、7階振型的應(yīng)力、位移云圖，以及對(duì)應(yīng)階的X，Y方向的最大值，見表1。由于6階X，Y方向的應(yīng)力、位移云圖總共有24張，僅給出第3階及其第1主應(yīng)力云圖作為示例(見圖6～圖10)。從圖10看出，最大拉應(yīng)力出現(xiàn)在壩踵，符合實(shí)際情況。

表1 碾壓混凝土重力壩1階～5階、7階振型X，Y位移、應(yīng)力最大值

圖6 碾壓混凝土重力壩第3階振型X方向位移圖

圖7 碾壓混凝土重力壩第3階振型Y方向位移圖

從表1可以看出，位移值滿足《混凝土重力壩設(shè)計(jì)規(guī)范》［11］要求。應(yīng)力值出現(xiàn)最大拉應(yīng)力1.73 MPa，這是規(guī)范［11］所不允許的，可以通過改進(jìn)壩體設(shè)計(jì)，或者施工時(shí)加強(qiáng)這些部位?？傊雺夯炷林亓蔚卣鸱磻?yīng)譜分析獲得的各階振型的最大、最小的應(yīng)力和位移值的具體位置，可以通過ANSYS查詢到，這無論對(duì)設(shè)計(jì)、施工還是后期的運(yùn)行，都有指導(dǎo)意義。

圖8 碾壓混凝土重力壩第3階振型X方向應(yīng)力圖

圖9 碾壓混凝土重力壩第3階振型Y方向應(yīng)力圖

圖10 碾壓混凝土重力壩第3階振型第1主應(yīng)力圖

5 結(jié)語

1)選用Plane82單元，建立了碾壓混凝土重力壩ANSYS有限元法模型。2)以每個(gè)振型作為單個(gè)質(zhì)點(diǎn)，采用子空間迭代法進(jìn)行模態(tài)分析，獲得了碾壓混凝土重力壩9階自振頻率，并給出對(duì)應(yīng)的振型圖。3)按照規(guī)范要求計(jì)算對(duì)應(yīng)的地震荷載反應(yīng)譜。反應(yīng)譜spectrum分析與模態(tài)擴(kuò)展modal expass分析相結(jié)合，得到壩體地震荷載下反應(yīng)譜分析的1階～5階、7階振型下的X，Y方向最大、最小的應(yīng)力和位移值。所得結(jié)果對(duì)碾壓混凝土重力壩的設(shè)計(jì)、施工以及后期運(yùn)行有著重要的工程實(shí)際意義。

［1］ 霍福山，楊建軍.大峽水電站碾壓混凝土重力壩施工［J］.東北水利水電，2011(1):21-24.

［2］ 石義剛，程志華.高碾壓混凝土壩關(guān)鍵施工技術(shù)［J］.中國(guó)水能及電氣化，2011(10):43-46.

［3］ 劉東海，李丙揚(yáng)，崔 博.高碾壓混凝土壩智能碾壓理論研究［J］.中國(guó)工程科學(xué)，2011(12):74-79.

［4］ 焦陽太.龍華口碾壓混凝土筑壩施工技術(shù)［J］.山西水利科技，2011(2):27-29.

［5］ 孫恭堯.高碾壓混凝土重力壩［M］.北京:中國(guó)電力出版社，2004.

［6］ 陳明陽，范書立，陳健云，等.官地碾壓混凝土重力壩抗震措施研究［J］.水力發(fā)電，2010(4):40-42.

［7］ 鄧良軍，張 杰，李雙寶.金安橋水電站碾壓混凝土重力壩抗震研究［J］.工程抗震與加固改造，2005(3):187-191.

［8］ 黃宜勝，常曉林，李建林.切割式橫縫碾壓混凝土重力壩抗震安全性研究［J］.巖土力學(xué)，2008(5):47-51.

［9］ GB 50011-2010，建筑抗震設(shè)計(jì)規(guī)范［S］.

［10］ SL 203-1997，水工建筑物抗震設(shè)計(jì)規(guī)范［S］.

［11］ SL 319-2005，混凝土重力壩設(shè)計(jì)規(guī)范［S］.