燈泡貫流式機組廠房動力特性及地震響應分析

黃海龍 王顯治

(中國電建集團昆明勘測設計研究院有限公司，云南 昆明 650051)

燈泡貫流式機組廠房動力特性及地震響應分析

黃海龍 王顯治

(中國電建集團昆明勘測設計研究院有限公司，云南 昆明 650051)

采用三維有限單元法，結合在建的老撾某水電站，對燈泡貫流式機組廠房的動力特性及地震響應進行了分析研究，并對地震加速度進行頻譜分析，從頻域角度對結構的地震響應做了初步的探索，為電站設計及類似工程提供參考。

河床式廠房，動力特性，有限元法，燈泡貫流式機組，地震響應

0 引言

隨著中、高水頭電站開發(fā)到一定程度，低水頭水力資源逐漸得到重視，燈泡貫流式水電站以其在投資、電站布置和運行等多方面的優(yōu)越性，已成為開發(fā)低水頭水力資源的優(yōu)良機型[1，2]。燈泡貫流式機組廠房由于其結構較特殊，內(nèi)力情況復雜，動力特性以及在地震作用下的動態(tài)響應等問題，需要通過理論與數(shù)值研究確定。為了較準確的分析燈泡貫流式機組廠房結構動力特性和動力響應，本文結合在建的大型燈泡貫流式水電站廠房，將梁板柱系統(tǒng)、副廠房、尾水管及廠房主體結構作為整體研究，采用結構動力分析方法，借助ANSYS有限元軟件，對該型式廠房的動力特性和地震響應的特點開展了一些初步的研究，以期為設計及類似工程提供參考。

1 工程概況

該水電站位于老撾湄公河一支流上，樞紐主要由河床式廠房和兩岸堤壩及壩上非常溢洪道組成。河床式廠房共布置4臺65 MW燈泡貫流式機組，單機單縫，廠房空間孔洞跨度較大，轉輪直徑達7.2 m。

2 計算模型及材料參數(shù)

2.1 模型范圍

計算模型以廠房標準機組段為對象，該機組段長度21.34 m，上下游寬度為81.13 m，主廠房最大高度59.40 m。基巖向上、下游方向各延伸120 m，底部向下延伸120 m。

模型對廠房流道、周邊混凝土及上部結構等均按實際尺寸建模，兩側混凝土的邊界因機組分縫按自由面考慮，基巖的各邊界施加法向約束。

整體模型和廠房網(wǎng)格分別見圖1，圖2。

2.2 材料參數(shù)

計算模型采用的主要材料參數(shù)如表1所示。

3 動力特性分析

因廠房上、下部結構的質量和剛度差異比較大，采用2種模型進行結構的模態(tài)分析：

1)整體模型，見圖1；2)廠房下部結構模型，見圖2，上部結構用Mass21質量單元模擬。

流道內(nèi)水體用Mass21質量單元模擬，計算采用分塊Lanczos迭代法。

根據(jù)計算成果，將2種模型的前20階自振頻率提取列于表2。

根據(jù)結構自振頻率和相應模態(tài)結果可知:

1)對整體模型，低階振型主要體現(xiàn)進水口及上部排架柱、墻的振型，說明排架柱剛度較低，這主要是排架柱跨度較大而構件截面較小所致。同時，有較多階也都是發(fā)生在這些部位，因此這是廠房結構中抗振的重點處理部位，建議設計時在進水口上部排架柱間布置混凝土墻，增強各向剛度。進水口邊墻從高階(14階)才開始橫河向振動，說明本電站進水口攔污柵下游混凝土橫梁對進口邊墻和中墩支承作用較大，提高了進口高邊墻、高中墩的側向剛度。

2)對廠房下部結構模型，前4階的振型主要表現(xiàn)為下部結構沿不同方向的振動，分別為橫河向、順河向、水平扭轉和鉛直向，說明廠房下部結構水平兩向剛度較小，同時還應特別注意由于結構質量、剛度分布不均勻等引起的扭轉作用對結構抗振性能的削弱。該廠房設計時在各機組段副廠房兩側分縫處設置了混凝土隔墻，增加了副廠房的順河向剛度，使副廠房樓板、柱和四周邊墻形成了整體箱形結構，高階振型也多表現(xiàn)為副廠房的整體振動。

4 地震響應分析

為全面了解廠房結構在地震作用下的動力響應規(guī)律，使結構設計更加經(jīng)濟合理，本節(jié)采用動力時程法，對本廠房結構進行了瞬態(tài)動力時程分析，對其總體抗震性能進行綜合分析與評價。

4.1 計算模型與假定

1)計算假定：

a.動水壓力采用Westergaard公式，以壩面附加質量的形式考慮；

b.為消除基礎對地震的放大作用，以無質量地基底部均勻輸入的方法考慮廠房和地基間的動力相互作用和地震波的輸入；

c.假定結構受水平雙向地震作用。

4.2 地震加速度時程曲線

本文選取典型地震波Taft波前20 s加速度數(shù)據(jù)，并將其加速度最大幅值調(diào)整為0.1g(本電站MDE=0.1g)，調(diào)整后的水平向Taft波見圖3,圖4。

4.3 動力時程分析成果

為分析廠房結構在單獨地震作用下的時程響應情況，選取結構部分特征點(1點～10點)進行分析，位置見圖5，同時將廠房結構各特征點響應參數(shù)最值統(tǒng)計于表3。

分析表3中計算結果可知：

1)在結構動位移響應方面，順河向動位移響應最大，發(fā)生在排架頂點(1點)，峰值為28.17 mm，發(fā)生時刻滯后其順河向激勵峰值發(fā)生時刻(3.72 s)，滯后現(xiàn)象較為明顯；橫河向動位移在進水口頂部(8點)達到峰值21.05 mm，上部排架及胸墻部位也出現(xiàn)較大橫河向動位移，亦是墻體橫河向剛度偏小的緣故。由此可知，地震作用下，動位移響應較大的部位主要分布在廠房上部結構、進水口邊墻及胸墻位置。

2)在結構動加速度響應方面，響應規(guī)律與動位移響應規(guī)律基本一致，結構順河向、橫河向動加速度在上部排架頂點(1點)最大，較輸入波加速度峰值放大約6倍和10倍，主要是廠房上部結構高程較高，各向質量與剛度相對弱小，受到了經(jīng)過下部主體結構放大后的地震加速度，因而動力響應遠超過下部結構的動力響應，發(fā)生了明顯的“鞭梢效應”；另外，由于廠房進水口邊墻橫河向剛度較弱，橫河向加速度在邊墻最頂部(8點)處也較大。

3)在結構動應力響應方面，第一主應力響應峰值發(fā)生在6點處，即發(fā)電機吊物孔角點，數(shù)值為4.03 MPa；另外各向動應力較大值還出現(xiàn)在管型座及進口邊墻底部，超過混凝土動態(tài)設計抗拉強度，因此建議設計時對上述局部部位加強配筋，提高其抗震能力。

4.4 地震加速度頻譜分析

結構在地震作用下的響應，一方面取決于結構自身的動力特性，另一方面還與地震作用的大小及地震波的頻譜特性密切相關。

1)輸入地震波的頻譜特性。

本節(jié)采用快速傅立葉變換法(FFT)，對地震加速度進行頻譜分析。針對調(diào)整后的Taft波進行傅立葉變換，得到了兩個方向地震波變換后的幅度譜，見圖6，圖7。各向地震波相應的主頻和主頻域分布見表4。

由頻譜分析可看出，調(diào)整后的Taft波的頻率分布范圍較寬，優(yōu)勢頻率主要集中在0 Hz～10 Hz，故此地震波能量集中在低頻部分，是一個以低頻為主的地震波。

由廠房動力特性分析知，整體結構的前20階自振頻率分布較密集，且主要集中在1 Hz～10 Hz之間，與地震優(yōu)勢頻率比較接近，尤其是整體結構的前3階自振頻率與地震波的主頻非常接近，振型主要表現(xiàn)為廠房上部排架柱的振動，故這些薄弱部位易在地震作用下產(chǎn)生較大的響應，甚至可能發(fā)生明顯的“鞭梢效應”。

2)地震響應加速度的頻譜分析。

為進一步分析輸入地震波的頻譜特性與結構地震反應的關系，選取廠房結構薄弱部位的特征點加速度響應頻譜進行分析，即廠房進口邊墩頂點(8點)和排架中部頂點(1點)，相應的加速度響應頻譜見圖8，圖9，各特征點加速度響應峰值及對應響應主頻見表5。

分析以上結果可知：

特征點的地震加速度響應主頻均在輸入波的優(yōu)勢頻率范圍內(nèi)(0 Hz～10 Hz)。廠房上部排架的基頻均較低，與輸入波的主頻較接近，且其響應主頻也主要在輸入波能量集中的低頻部分(0 Hz～5 Hz)，因而這些部位的地震響應比較大，與前述結論一致。

副廠房樓板在地震作用下主要關注其豎向振動響應，經(jīng)分析，樓板振動主要出現(xiàn)在第11階后，自振頻率相對較高，其第1階自振頻率約12.09 Hz，錯開了輸入Taft波能量較為集中的低頻部分(0 Hz～5 Hz)，加之本電站副廠房樓板四周均有墻體支撐，形成箱體結構，提高了整體剛度，故地震作用下的振動反應也相對較小，但副廠房是運行人員工作及各類儀器設備布置的場所，其振動幅值也應引起足夠重視。

5 結論與建議

1)燈泡貫流式機組廠房上部排架柱結構部位剛度最薄弱，是廠房結構中抗振的重點關注部位，必須采取有效措施防止發(fā)生共振，可考慮改善結構增加剛度，如柱間砌體墻改為混凝土墻結構，適當縮小柱間距及加大柱斷面；加強廠房上部結構之間的聯(lián)系，如將屋頂桁架與排架柱鉸接，將上下游排架柱相關聯(lián)，提高上部結構順河向剛度。另外，進口中墩及邊墻一般較高，側向剛度較小，其自振頻率亦不高，也應引起抗振重視，建議可在攔污柵下游適當高程處增加橫梁支撐，提高進水口橫河向剛度，進而提高廠房整體抗振性能。

2)燈泡貫流式機組廠房結構在地震作用下，也表現(xiàn)出流道中段大孔洞位置各向剛度較小，地震響應應力相對較大，建議對應力響應較大的部位，尤其是管型座上下支柱區(qū)域，還需承受機組振動等影響，需加強局部配筋，提高其抗振能力，以利燈泡體的安全。

3)通過對輸入地震波和地震響應加速度的頻譜分析，所得結論與動力時程法分析所得的結論是一致的，說明通過對結構自振特性分析和對輸入波頻譜分析，可實現(xiàn)對廠房結構進行抗震薄弱部位定性分析，亦可作為結構抗震分析結論的補充和論證依據(jù)。

[1] 潘家錚.大力開發(fā)水電:中國持續(xù)發(fā)展的必由之路[J].中國三峽，2009(1)：5-12.

[2] 劉志超.開發(fā)低水頭電站的優(yōu)良機型[J].機械研究與應用，1995(3)：39-41.

[3] 練繼建，王海軍，王日宣.河床式水電站廠房結構動力特性研究[J].水利水電技術，2004，35(8)：37-40.

[4] 蘇禮邦，劉云賀，李守義.燈泡貫流式水電站廠房壩段的抗震分析[J].西北農(nóng)林科技大學學報(自然科學版)，2007，35(7)：202-206.

[5] 單智杰.燈泡貫流式水電站廠房三維靜動力分析[D].西安：西安理工大學，2005.

[6] 歐陽金惠，陳厚群，張超然,等.三峽電站15#機組廠房結構動力分析[J].中國水利水電科學研究院學報，2007，5(2)：137-142.

[7] 魏 來，劉馬群，李 奎.信號分析方法及其在地震波處理中的應用[J].洛陽大學學報，2007，22(4)：69-73.

Abstract： Based on a hydropower station which is under-construction in Laos, a dynamic characteristics research and a preliminary discussion of seismic response of the bulb tubular hydropower plant has been conducted by means of 3D finite element analysis method. Besides, frequency spectrum analysis of earthquake acceleration is performed based on the result of dynamic time-history analytical method, then, a preliminary study on the structure seismic response from the perspective of frequency domain is put forward, which will be helpful for the plant design.

Keywords： riverbed-hydropower house, dynamic characteristics, finite element method, bulb tubular unit, seismic response

Studyondynamiccharacteristicsandseismicresponseofthebulbtubularpowerhouse

HuangHailongWangXianzhi

(ChinaPowerGroupKunmingEngineeringCorporationLimited,Kunming650051,China)

1009-6825(2017)25-0043-03

TU318

A

2017-06-23

黃海龍(1987- ),男，碩士，工程師； 王顯治(1983- ),男,碩士,高級工程師