黑 燦，伍鶴皋，石長征，傅 丹(武漢大學(xué)水資源與水利水電工程科學(xué)國家重點(diǎn)實(shí)驗(yàn)室，武漢 430072)

0 引 言

地震由于其突發(fā)性和難預(yù)測性，在眾多自然災(zāi)害中堪稱群災(zāi)之首。我國是地震多發(fā)國，并且我國的水電開發(fā)多集中在西南地震高發(fā)區(qū)，汶川地震水電工程震害調(diào)查發(fā)現(xiàn)，水電站地面廠房下部大體積混凝土的震損較輕微，上部結(jié)構(gòu)的震損較重[1,2]，而上部結(jié)構(gòu)的地震響應(yīng)與屋面結(jié)構(gòu)形式關(guān)系密切。近些年隨著水電站水輪發(fā)電機(jī)組單機(jī)容量不斷加大，廠房跨度也隨之增大，網(wǎng)架結(jié)構(gòu)的廠房屋面優(yōu)勢愈來愈顯著，已成為未來的發(fā)展趨勢[3,4]。因此，對(duì)于網(wǎng)架屋面水電站廠房的抗震問題應(yīng)有足夠的重視。而對(duì)這種上部為鋼網(wǎng)架下部為混凝土的混合結(jié)構(gòu)體系進(jìn)行設(shè)計(jì)時(shí)，往往將網(wǎng)架和下部支承結(jié)構(gòu)分開獨(dú)立設(shè)計(jì)，沒有考慮兩者之間的相互作用，可能導(dǎo)致設(shè)計(jì)偏于不安全[5]。

目前，對(duì)于網(wǎng)架和下部支承結(jié)構(gòu)相互作用的研究多針對(duì)體育館、火車站等空曠建筑[6,7]，對(duì)水電站廠房的研究較少。在研究體育館等大跨度混合結(jié)構(gòu)的抗震性能時(shí)發(fā)現(xiàn)，網(wǎng)架屋面結(jié)構(gòu)在地震作用下的破壞多為網(wǎng)架桿件破壞，或下部支承連接處破壞導(dǎo)致網(wǎng)架塌落[8]。并且有研究表明包含網(wǎng)架的結(jié)構(gòu)在地震作用下的響應(yīng)與傳統(tǒng)結(jié)構(gòu)有很大不同[9]，所以有必要將網(wǎng)架放置于下部混凝土支承上，對(duì)兩者的相互作用進(jìn)行研究。在水電領(lǐng)域，研究網(wǎng)架屋面廠房的抗震性能常采用簡化的模型進(jìn)行有限元模擬計(jì)算：最初一般將網(wǎng)架結(jié)構(gòu)簡化為板單元，而板單元不但不能真實(shí)模擬網(wǎng)架本身的動(dòng)力特性，還將影響墻體的動(dòng)力特性；隨著研究的深入，將網(wǎng)架用連接在上下游墻之間的剛性二力桿進(jìn)行模擬[10]，但剛性二力桿的剛度顯然比網(wǎng)架實(shí)際剛度大很多，這種剛度差異對(duì)下部混凝土結(jié)構(gòu)地震響應(yīng)的影響還有待研究；另外有些研究中由于資料不足，只能采用不考慮網(wǎng)架的上部開洞模型進(jìn)行簡化計(jì)算。不同的簡化方式對(duì)廠房的動(dòng)力特性有何種影響、其計(jì)算精度能否滿足要求，這些都是值得研究的問題。此外，網(wǎng)架與下部支承連接方式的不同對(duì)廠房地震響應(yīng)及網(wǎng)架受力的影響如何、采用何種連接方式比較合理，也需要進(jìn)一步探討。本文將以某水電站地面廠房為工程背景，借助ANSYS軟件，采用時(shí)程分析法，對(duì)網(wǎng)架屋面的模擬方法和支承方式對(duì)地面廠房抗震性能的影響展開研究。

1 計(jì)算方案和條件

某水電站地面廠房順河向長52.00 m，單個(gè)機(jī)組段寬26.50 m，高64.35 m，水平向設(shè)計(jì)地震加速度0.176 g。以一個(gè)中間標(biāo)準(zhǔn)機(jī)組段為研究對(duì)象建立整體有限元模型，模型包括廠房和地基，基巖范圍由廠房向上、下游側(cè)分別取100 m，基巖深度約為130 m。

在計(jì)算范圍內(nèi)，對(duì)主廠房上下游墻(排架柱)、風(fēng)罩、機(jī)墩、蝸殼以及尾水管等均按實(shí)際尺寸進(jìn)行模擬，網(wǎng)架采用桿單元LINK180，蝸殼、尾水管鋼襯、座環(huán)等鋼結(jié)構(gòu)采用殼單元SHELL63，混凝土和墊層采用實(shí)體單元SOLID45，機(jī)組重量采用質(zhì)量單元MASS21模擬在相應(yīng)位置。

根據(jù)以往研究，水電站地面廠房結(jié)構(gòu)在順河向抗震性較差，本文主要關(guān)注廠房在順河向的地震響應(yīng)。計(jì)算采用時(shí)程分析法，選用Koyna實(shí)測波，最大幅值調(diào)整為0.176 g，時(shí)間長度10 s，步長0.02 s。結(jié)構(gòu)的阻尼比根據(jù)《水工建筑物抗震設(shè)計(jì)規(guī)范》DL 5073-2000[11]取0.05。

對(duì)于屋面網(wǎng)架，參考工程中常用的網(wǎng)架模擬方式以及連接方式，本文共對(duì)比分析了以下5種方案，以期為網(wǎng)架式屋面結(jié)構(gòu)的計(jì)算提供參考，具體計(jì)算方案見表1。

表1 計(jì)算方案Tab.1 Calculation schemes

注：①方案A、B、C中的網(wǎng)架自重按照《空間網(wǎng)格結(jié)構(gòu)技術(shù)規(guī)程》JGJ7-2010[12]進(jìn)行估算；②方案C彈性連桿剛度根據(jù)實(shí)際工程情況確定。

整體模型網(wǎng)格和廠房模型網(wǎng)格如圖1和圖2所示，圖1、圖2均為方案E的模型，不同方案模型的差別僅在于網(wǎng)架模擬及支承方式的不同。網(wǎng)架模型如圖3所示。

圖1 整體模型網(wǎng)格Fig.1 The overall model meshes

圖2 廠房模型網(wǎng)格Fig.2 Power house model meshes

圖3 網(wǎng)架模型Fig.3 Grid model

計(jì)算邊界條件為模型底部基巖的各個(gè)面施加法向約束，蝸殼進(jìn)水管處鋼管橫截面施加法向(管軸向)約束，其他邊界均為自由邊界。

模型涉及的材料共4種，分別為混凝土、鋼材、墊層和基巖，機(jī)組段各部位的混凝土設(shè)計(jì)強(qiáng)度等級(jí)均為C25，廠房基巖采用無質(zhì)量地基，其變形模量采用Ⅲ1類巖體對(duì)應(yīng)的變形模量值，所有材料參數(shù)見表2。

表2 模型材料計(jì)算參數(shù)Tab.2 Model material calculation parameters

2 地震響應(yīng)對(duì)比分析

考慮到網(wǎng)架模擬方式的不同對(duì)廠房順河向的剛度影響較大，本節(jié)將采用時(shí)程法對(duì)比分析各方案在順河向地震作用下的響應(yīng)。

2.1 混凝土結(jié)構(gòu)地震響應(yīng)

為分析廠房混凝土結(jié)構(gòu)在不同網(wǎng)架模擬方式和支承方式下的地震響應(yīng)，本節(jié)整理了典型部位的順河向動(dòng)位移和豎向動(dòng)應(yīng)力，各特征點(diǎn)位置見圖4。相應(yīng)特征點(diǎn)的順河向動(dòng)位移峰值和豎向動(dòng)應(yīng)力峰值見圖5、圖6。

1-上游立柱頂部；2-下游立柱頂部；3-上游立柱與牛腿交接處；4-下游立柱與牛腿交接處；5-上游墻與發(fā)電機(jī)層樓板交接處；6-下游墻與發(fā)電機(jī)層樓板交接處；7-風(fēng)罩內(nèi)側(cè)；8-定子基礎(chǔ)內(nèi)側(cè)；9-上游墻水輪機(jī)層高層外側(cè)；10-尾水平臺(tái)樓板跨中圖4 特征點(diǎn)示意圖Fig.4 Typical locations of dynamic displacement

圖5 特征點(diǎn)順河向動(dòng)位移峰值Fig.5 The peak downstream directional dynamic displacement of typical locations

圖6 特征點(diǎn)豎向動(dòng)應(yīng)力峰值Fig.6 The peak vertical dynamic stress of typical locations

5種方案各特征點(diǎn)的地震響應(yīng)隨著高程降低基本上呈下降趨勢，上下游立柱頂部(1點(diǎn)、2點(diǎn))以及上游立柱與牛腿交接處(3點(diǎn))的順河向動(dòng)位移較大，上下游立柱與牛腿交接處(3點(diǎn)、4點(diǎn))以及上游墻與發(fā)電機(jī)層樓板交接處(5點(diǎn))的豎向動(dòng)應(yīng)力較大。下游立柱與牛腿交接處(4點(diǎn))的動(dòng)位移以及下游墻與發(fā)電機(jī)層樓板交接處(6點(diǎn))的動(dòng)應(yīng)力，由于下游副廠房的支撐作用，其地震響應(yīng)比上游側(cè)相應(yīng)響應(yīng)量要小。比較5種方案的結(jié)果可以發(fā)現(xiàn)，網(wǎng)架的模擬方式對(duì)廠房下部大體積混凝土影響很小，對(duì)上游立柱的影響比較明顯。

為了更準(zhǔn)確、更全面地反映地震作用下的全過程情況，提取了1、2、8點(diǎn)的順河向動(dòng)位移時(shí)程曲線見圖7。

圖7 特征點(diǎn)順河向動(dòng)位移時(shí)程Fig.7 The downstream directional dynamic displacement time-history curve of typical locations

由圖5-7可知：

(1)位于廠房上部結(jié)構(gòu)的1點(diǎn)和2點(diǎn)，各方案的動(dòng)位移時(shí)程曲線相差較大，位于下部大體積混凝土結(jié)構(gòu)的8點(diǎn)，各方案的動(dòng)位移時(shí)程曲線相差較小，說明不同的網(wǎng)架模擬方式及支承方式對(duì)廠房上部結(jié)構(gòu)的地震響應(yīng)影響較大，對(duì)下部大體積混凝土結(jié)構(gòu)影響較小。

(2)采用簡化模型計(jì)算的方案B、C與按照實(shí)際尺寸模擬計(jì)算的方案D、E的時(shí)程曲線相差很大，特別是下游側(cè)2點(diǎn)，方案B、C的地震響應(yīng)明顯大于方案D、E的地震響應(yīng)，峰值相差達(dá)50%～90%。除方案B外，各方案1點(diǎn)的動(dòng)位移均明顯大于2點(diǎn)的動(dòng)位移，這是由于上游立柱沒有水平支撐，鞭梢效應(yīng)更為明顯；而方案B將網(wǎng)架簡化為剛性二力桿，1點(diǎn)和2點(diǎn)的順河向動(dòng)位移時(shí)程基本一樣，這與實(shí)際情況不符。說明將網(wǎng)架簡化為連桿支撐，特別是簡化為剛性二力桿時(shí)，計(jì)算結(jié)果與實(shí)際情況相差較大，故不建議將網(wǎng)架簡化為連桿支撐進(jìn)行計(jì)算。

(3)方案A與方案D、E的時(shí)程曲線較為接近，特別是與方案E的時(shí)程曲線基本一樣，各特征點(diǎn)的響應(yīng)峰值與方案E響應(yīng)峰值的差別除3點(diǎn)動(dòng)位移(15%)和5點(diǎn)動(dòng)應(yīng)力(17%)較大外，其余均在10%左右或者更小，并且方案A的響應(yīng)峰值除4點(diǎn)的動(dòng)位移和6點(diǎn)的動(dòng)應(yīng)力稍小于(差別小于2%)方案E的相應(yīng)響應(yīng)峰值外，其余均大于方案E的響應(yīng)峰值。說明采用不考慮網(wǎng)架的模型簡化計(jì)算時(shí)與實(shí)際的差別相對(duì)較小，特別是當(dāng)網(wǎng)架結(jié)構(gòu)與下部支承為一端固定鉸支座一端滾軸支座連接時(shí)，采用不考慮網(wǎng)架的簡化模型能基本滿足精度要求。

2.2 網(wǎng)架結(jié)構(gòu)地震響應(yīng)

按照實(shí)際尺寸模擬網(wǎng)架時(shí)，兩端固定鉸支座方案和一端固定鉸支座一端滾軸支座方案的網(wǎng)架軸向應(yīng)力最大值均出現(xiàn)在跨中下弦桿處。本節(jié)整理了方案D和方案E在順河向地震作用下網(wǎng)架跨中下弦桿的軸向應(yīng)力時(shí)程曲線，見圖8。并考慮自重及屋面荷載所引起的軸向應(yīng)力，將地震和靜力共同作用的組合結(jié)果列于表4。

圖8 跨中下弦桿地震作用下軸向應(yīng)力時(shí)程Fig.8 The axial stress time-history curve of mid-span bottom chord bar under seismic excitation

由圖8和表3可知：

(1)在單獨(dú)地震作用下，方案D的軸向應(yīng)力明顯大于方案E的軸向應(yīng)力，方案D網(wǎng)架的動(dòng)應(yīng)力較大主要是由于上下游立柱順河向相對(duì)變形較大，說明相比較于兩端均為固定鉸支座的連接方式，一端為固定鉸支座一端為滾軸支座時(shí)網(wǎng)架的順河向抗震性能較好。但應(yīng)注意滾軸支座的水平約束較弱，在自重及屋面荷載等靜荷載作用下產(chǎn)生的軸向應(yīng)力較大，考慮靜力和地震組合作用后，一端固定鉸支座一端滾軸支座連接的網(wǎng)架桿件軸向應(yīng)力大于兩端均為固定鉸支座的網(wǎng)架桿件軸向應(yīng)力，在具體設(shè)計(jì)時(shí)應(yīng)注意靜力及地震的共同作用。

(2)E方案中網(wǎng)架桿件軸向應(yīng)力的組合結(jié)果達(dá)407 MPa，即使是應(yīng)力較小的D方案也達(dá)308 MPa，若網(wǎng)架采用屈服強(qiáng)度較低的Q235鋼材，在靜力和地震的共同作用下網(wǎng)架很容易進(jìn)入塑性階段，甚至有可能拉斷導(dǎo)致屋頂垮塌，即使采用Q345鋼材，網(wǎng)架也有發(fā)生塑性屈服的可能。所以，對(duì)于網(wǎng)架屋面的水電站廠房，網(wǎng)架結(jié)構(gòu)是抗震中的薄弱部位，并且不同支承方式的網(wǎng)架受下部支承的影響也不同，在抗震設(shè)計(jì)中應(yīng)著重進(jìn)行分析。

3 結(jié) 語

(1)不同的網(wǎng)架模擬方式及支承方式對(duì)廠房下部大體積混凝土結(jié)構(gòu)的地震響應(yīng)影響較小，但對(duì)上部結(jié)構(gòu)的影響較大；將網(wǎng)架簡化為連桿支撐，特別是簡化為剛性二力桿時(shí)，計(jì)算結(jié)果與實(shí)際情況相差較大，不建議將網(wǎng)架簡化為連桿支撐進(jìn)行計(jì)算。

(2)相比較于連桿支撐(方案B和C)，不考慮網(wǎng)架的簡化模型(方案A)與實(shí)際網(wǎng)架一端固定鉸支座一端滾軸支座(方案E)的計(jì)算結(jié)果的差別相對(duì)較小，因此當(dāng)缺少網(wǎng)架資料，特別是當(dāng)網(wǎng)架結(jié)構(gòu)與下部支承為一端固定鉸支座一端滾軸支座連接時(shí)，采用不考慮網(wǎng)架的簡化模型能基本滿足精度要求，可采用不考慮網(wǎng)架的模型簡化計(jì)算。

(3)網(wǎng)架結(jié)構(gòu)是水電站廠房抗震中的薄弱部位，并且其地震響應(yīng)受下部支承影響較大，在研究網(wǎng)架的抗震特性時(shí)，建議建立包含下部支承的整體模型。一端固定鉸支座一端滾軸支座的網(wǎng)架連接方式相對(duì)于兩端均為固定鉸支座的連接方式，前者網(wǎng)架的地震響應(yīng)較小，但考慮靜力的組合作用后，前者的網(wǎng)架桿件軸力反而較大，故選取網(wǎng)架與下部支承的連接方式時(shí)應(yīng)注意地震和靜力的組合作用。

□

[1] 晏志勇，王 斌，周建平，等. 汶川地震災(zāi)區(qū)大中型水電工程震損調(diào)查與分析[M]. 北京：中國水利水電出版社，2009.

[2] 李錦成，于麗偉. 河床式水電站地震荷載分布規(guī)律研究 [J]. 中國農(nóng)村水利水電，2007，(3)：75-78.

[3] 董石麟，邢 棟，趙 陽. 現(xiàn)代大跨空間結(jié)構(gòu)在中國的應(yīng)用與發(fā)展 [J]．空間結(jié)構(gòu)，2012,18(1)：3-16．

[4] 馬 濤，曾 旭，蔡 波. 水電站廠房屋頂結(jié)構(gòu)型式研究 [J]．水資源與水工程學(xué)報(bào)，2014，25(6):197-202．

[5] 曹 資，張 超，張毅剛，等. 網(wǎng)殼屋蓋與下部支承結(jié)構(gòu)動(dòng)力相互作用研究 [J]. 空間結(jié)構(gòu)，2001，7(2)：19-26.

[6] 張明文. 平板網(wǎng)架與下部支承結(jié)構(gòu)整體工作分析研究[D]. 西安：西安建筑科技大學(xué)，2010.

[7] 田承昊，米宏廣，王彥芳，等. 哈大鐵路客運(yùn)專線鞍山站站房屋蓋與下部支承結(jié)構(gòu)的動(dòng)力相互作用研究 [J]. 鐵道標(biāo)準(zhǔn)設(shè)計(jì)，2012，(5)：154-157.

[8] Toshitsugu Saka, Yoshiya Taniguchi. Damage to Spatial Structures by the 1995 Hyogoken-Nanbu Earthquake in Japan [J]. International Journal of Space Structures， 1997，3(12)：125-133.

[9] Hassan A.Moghaddam. Seismic Behaviour of Space Structures [J]. International Journal of Space Structures，2000，15(2)：119-135.

[10] 馬震岳，張運(yùn)良，陳 婧，等. 水電站廠房和機(jī)組耦合動(dòng)力學(xué)理論及應(yīng)用[M]. 北京：中國水利水電出版社，2013.

[11] DL 5073-2000，水工建筑物抗震設(shè)計(jì)規(guī)范[S].

[12] JGJ7-2010，空間網(wǎng)格結(jié)構(gòu)技術(shù)規(guī)程[S].