曾祥華

（上?？睖y設(shè)計研究院，上海 200434）

1 工程概況

大溪水庫建于1958年11月，位于溧陽市區(qū)西南的觀山、黃家山、愚公山一帶低山丘陵區(qū)，庫區(qū)距溧陽市13 km。水庫集水面積為90 km2，水庫總庫容1.03億m3。水庫以防洪、灌溉為主，兼顧供水、養(yǎng)殖等，與溧陽沙河水庫、宜興市橫山水庫一起成為太湖湖西宜溧山區(qū)的3座大型水庫，屬太湖西南河水系上游的骨干防洪工程。

大溪水庫為大（2）型水庫，樞紐工程等別為Ⅱ等，主要建筑物為2級建筑物。新建溢洪閘布置于烏龜山與前山之間的埡口位置，為3孔寬頂堰結(jié)構(gòu)，工作閘門為3扇弧形鋼閘門。

2 有限元分析模式與分析目的

巖基上的溢洪閘底板與基巖間存在結(jié)構(gòu)面，閘室邊墩與回填土體亦非連續(xù)介質(zhì)，在傳統(tǒng)的結(jié)構(gòu)有限元計算分析[1-3]中，底板和基巖以及邊墩與回填土體的結(jié)構(gòu)單元往往共用節(jié)點，導(dǎo)致當(dāng)邊墩受側(cè)向土壓力作用時，閘室底板與基巖之間以及閘墩與側(cè)向土體間可能出現(xiàn)拉應(yīng)力，應(yīng)力分析結(jié)果不能反應(yīng)結(jié)構(gòu)的真實受力狀態(tài)，甚至?xí)霈F(xiàn)錯誤的結(jié)果。本文在對大溪水庫溢洪閘的有限元分析中采用了接觸模型[4]，接觸單元可以較好地模擬混凝土與基巖、土體等介質(zhì)接觸面的壓緊和張開效應(yīng)。此種模型更接近結(jié)構(gòu)的真實變形狀態(tài)，計算結(jié)果也更為合理。

此外，溢洪閘本身的動力響應(yīng)也是結(jié)構(gòu)分析的重點及難點，本文所提的動力響應(yīng)既包含對地震動的響應(yīng)，也包括對溢洪時高速水流的脈沖振動響應(yīng)。動力響應(yīng)分析的前提是計算出結(jié)構(gòu)的各階自振頻率，進而分析閘室的某一階自振頻率是否會與脈沖水流發(fā)生共振。傳統(tǒng)的地震擬靜力分析方法不能計算閘室結(jié)構(gòu)的自振特性以及在地震作用下的動態(tài)響應(yīng)等問題，本文通過基于反應(yīng)譜的振型分解法[4]來計算溢洪閘對地震的動力響應(yīng)，分析結(jié)果更為合理。

3 有限元模型的建立

3.1 基本計算假定

1）假定庫水為不可壓縮流體，庫水對閘室的動力作用簡化為附加質(zhì)量；

2）假定溢洪閘混凝土為各向同性材料，基巖為均勻的無質(zhì)量各向同性材料；

3）邊墩兩側(cè)填土采用Drucker-Prager模型。

3.2 有限元模型范圍及單元剖分

取整個溢洪閘及邊墩兩側(cè)填土和部分地基建立三維有限元模型，實體結(jié)構(gòu)及土體全部采用計算精度較高的八節(jié)點六面體實體單元，主要結(jié)構(gòu)部位沿厚度方向剖分不少于3層，溢洪閘網(wǎng)格尺寸不大于0.5m，土體部位的網(wǎng)格尺寸適當(dāng)放大。

如何模擬土體與混凝土結(jié)構(gòu)之間的相互作用目前仍然是一個非常棘手的問題，一般采用無厚度的Goodman單元來模擬土體與結(jié)構(gòu)之間的相互作用[5]，本次計算采用面-面接觸模型，在土體與邊墩及底板之間的接觸面設(shè)置接觸單元，以模擬土體與混凝土結(jié)構(gòu)之間的接觸效應(yīng)。有限元模型見圖1。

圖1 閘室及地基有限元模型圖

圖中有限元模型應(yīng)用的坐標(biāo)系為：X向為順?biāo)鞣较颍琘向為鉛直向上，Z向為向右岸。有限元模型約束情況為：地基底面為三向約束，上、下游面和側(cè)面均為法向約束。

3.3 材料物理力學(xué)參數(shù)

計算模型中閘室混凝土和基巖采用均質(zhì)、各向同性的連續(xù)線彈性材料，邊墩兩側(cè)填土采用Drucker-Prager模型[5]材料。Drucker-Prager模型是在考慮靜水壓力的Drucker-Prager屈服準(zhǔn)則的基礎(chǔ)上建立起來的理想彈塑性模型，對于巖土，砂礫石等材料應(yīng)考慮靜水壓力、材料粘性及內(nèi)摩擦角的影響。Drucker-Prager本構(gòu)模型的材料特性有4個參數(shù)：K、G、a、k。K、G為彈性常數(shù)，可由彈性模量E和泊松比μ換算得出，塑性參數(shù)a、k由Mohr-Coulomb準(zhǔn)則的材料參數(shù)粘聚力c和內(nèi)摩擦角Φ換算得出。

有限元模型中采用的材料物理力學(xué)參數(shù)指標(biāo)見表1。

3.4 計算荷載及組合

靜力分析計算了閘室在正常擋水、設(shè)計洪水和校核洪水3個工況下的應(yīng)力應(yīng)變。荷載包括結(jié)構(gòu)自重、設(shè)備重量、上游水壓力、浪壓力、邊墩外側(cè)水壓力、揚壓力和弧門推力等。

貝多芬生活在古典主義與浪漫主義的過渡階段。他的作品不僅展現(xiàn)了古典義時期音樂的精髓，把古典主義音樂發(fā)展到了極致，同時他還為浪漫主義音樂開啟了先河。尤其在他晚期的作品中更是明顯的體現(xiàn)出了浪漫主義的傾向。也正因如此，后人稱是“集古典主義之大成，開浪漫主義之先河”的音樂大師。本文主要從貝貝多芬生活在古典主義與浪漫主義的過渡階段。他的作品不僅展現(xiàn)了古典主多芬的音樂思想、創(chuàng)作風(fēng)格、審美原則三方面對貝多芬的浪漫主義音樂進行探究。

地震動力分析計算采用的水位組合為上游正常擋水，下游無水。動力計算中庫水對閘室的作用[6]采用抗震規(guī)范[7]推薦的維斯特伽特(Westergaard)公式進行計算。

4 有限元靜力分析結(jié)果

閘室結(jié)構(gòu)在正常擋水工況下的第一主應(yīng)力云圖見圖2，其主要部位應(yīng)力分析如下。

在側(cè)向土壓力及邊墩外側(cè)水壓力的作用下，兩邊墩下游側(cè)底部外側(cè)面豎直向（Y向）出現(xiàn)較大的拉應(yīng)力，約為0.9MPa左右。胸墻處的邊墩外側(cè)面產(chǎn)生了順?biāo)飨颍╔向）1.0 MPa左右的拉應(yīng)力，邊墩下游側(cè)中部內(nèi)側(cè)面產(chǎn)生了豎直向（Y向）1.9 MPa左右的拉應(yīng)力。底板下游底面產(chǎn)生了橫河向（Z向）1.2 MPa左右的拉應(yīng)力。在設(shè)計洪水和校核洪水工況下，由于閘室內(nèi)水壓力的作用，上述拉應(yīng)力分別減小了0.2、0.1、0.4和0.3 MPa左右。胸墻基本都為壓應(yīng)力，且數(shù)值較小。其他部位的拉應(yīng)力亦較小。

表1 有限元計算材料參數(shù)

圖2 正常蓄水工況下X向第一主應(yīng)力云圖(MPa)

邊墩的牛腿根部由于受弧門水推力的作用及外側(cè)土、水壓力的作用，出現(xiàn)較大的拉應(yīng)力，圖3為牛腿根部順?biāo)鞣较蚶瓚?yīng)力沿閘墩厚度的分布圖，從圖中可以看出，混凝土表面拉應(yīng)力達到了3.0MPa，沿閘墩厚度方向遞減，邊墩內(nèi)表面0.23 m范圍內(nèi)的拉應(yīng)力均較大，超過1.2MPa。

圖3 牛腿根部第一主應(yīng)力應(yīng)力沿閘墩厚度的分布

5 動力反應(yīng)譜分析

5.1 模態(tài)分析

模態(tài)分析也被稱為振型疊加法動力分析，是線性結(jié)構(gòu)系統(tǒng)地震分析中最常用而且最有效的分析方法，是其他動力分析的基礎(chǔ)，包括反應(yīng)譜分析和時程分析。本文通過模態(tài)分析計算溢洪閘的各階振型及其特征周期、振型參與系數(shù)等動力參數(shù)，采用特征向量法求解。

在以往的水工結(jié)構(gòu)振動分析[8-9]中，多采用結(jié)構(gòu)自振頻率是否避開機組的強振頻率或外部激勵頻率的共振校核法。但是水工結(jié)構(gòu)均是多構(gòu)件組成的復(fù)雜結(jié)構(gòu)，各結(jié)構(gòu)構(gòu)件的頻率相差很大、頻率分布范圍廣,共振校核經(jīng)常難以達到預(yù)期效果。為全面分析溢洪閘自振特性，確定閘墩、胸墻上部結(jié)構(gòu)等部位的振型及其特征周期，本文采用Block lanczos法[10]計算了閘室結(jié)構(gòu)的前200階自振頻率，表2列出了前10階自振頻率、周期和各向振型參與系數(shù)。閘室結(jié)構(gòu)的前4階振型形狀見圖4～7。根據(jù)各階振型確定的結(jié)構(gòu)各主要構(gòu)件的自振頻率見表3。

表2 溢洪閘前10階自振頻率及振型參與系數(shù)

表3 閘室主要部位自振頻率

由模態(tài)分析結(jié)果可以看出，閘室結(jié)構(gòu)的第1階頻率（基頻）為8.467Hz，其前八階頻率在8.5～30.0Hz之間。

5.2 反應(yīng)譜計算

圖4 閘室結(jié)構(gòu)第一階振型（上游視圖）

圖5 閘室結(jié)構(gòu)第二階振型（上游視圖）

圖6 閘室結(jié)構(gòu)第三階振型（頂視圖）

圖7 閘室結(jié)構(gòu)第四階振型（右岸視圖）

圖8 設(shè)計反應(yīng)譜

6 反應(yīng)譜分析計算結(jié)果

從模態(tài)分析的結(jié)果可知，下游側(cè)中墩及上游側(cè)邊墩的自振周期在0.1 s和Tg(0.3 s)之間，因此地震影響系數(shù)最大，所以地震作用效應(yīng)很大。

邊墩下游側(cè)底部外側(cè)面產(chǎn)生了豎直向（Y向）1.8 MPa左右的拉應(yīng)力，胸墻處的邊墩外側(cè)面產(chǎn)生了順?biāo)飨颍╔向）1.1 MPa左右的拉應(yīng)力，邊墩下游側(cè)中部內(nèi)側(cè)面產(chǎn)生了豎直向（Y向）1.9 MPa左右的拉應(yīng)力。底板下游外側(cè)產(chǎn)生了橫河向（Z向）1.3 MPa左右的拉應(yīng)力。下游側(cè)交通橋頂面與邊墩接觸部位產(chǎn)生了橫河向（Z向）1.5 MPa左右的拉應(yīng)力。

由于地震作用，閘墩與底板的連接部位，交通橋與閘墩的連接部位等區(qū)域都產(chǎn)生了應(yīng)力集中，但應(yīng)力值都較小。

7 結(jié)語

1）對該溢洪閘的三維有限元分析表明，接觸分析模型能更合理的模擬閘室結(jié)構(gòu)的受力狀況，邊墩及底板的拉應(yīng)力極值較傳統(tǒng)算法大20%左右；

2）三維有限元靜、動力計算結(jié)果表明，邊墩，底板以及牛腿附近拉應(yīng)力較大，個別部位已超過混凝土設(shè)計抗拉強度，其余各構(gòu)件拉應(yīng)力較小，整個閘室結(jié)構(gòu)應(yīng)力分布均勻，溢洪閘閘室結(jié)構(gòu)設(shè)計基本合理，雖然有些部位拉應(yīng)力超過了混凝土的抗拉強度，但是通過結(jié)構(gòu)配筋都能滿足設(shè)計要求；

3）閘室結(jié)構(gòu)交通橋結(jié)構(gòu)初步選定為預(yù)制形式，通過有限元動力分析，邊墩的地震反應(yīng)較大，為此將下游交通橋改為現(xiàn)澆形式，使其與閘墩的連接更為可靠；

4）應(yīng)用大型通用有限元軟件可以拓展設(shè)計人員對復(fù)雜水工建筑物結(jié)構(gòu)計算的能力，計算模型比傳統(tǒng)結(jié)構(gòu)力學(xué)的方法更全面，計算結(jié)果也更為合理。

[1] 王平，何順賓，李永紅．福堂水電站攔河閘壩三維有限元靜力分析和動力響應(yīng)分析[J]．水電站設(shè)計，2002（6）：37-41．

[2] 彭成山，袁淮中，王甜．南水北調(diào)工程七里河水閘的三維有限元分析[J]．華北水利水電學(xué)院學(xué)報，2010（6）：12-13．

[3] 孟旭央，陳軍芳．西霞院工程泄洪閘閘室結(jié)構(gòu)應(yīng)力分析[J]．人民黃河，2006（8）：63-64．

[4] 朱伯芳．有限單元法原理與應(yīng)用[M]．2版．北京：中國水利水電出版社，1998：517-519．

[5] 錢家歡，殷宗澤.土工原理與計算[M]．2版．北京：中國水利水電出版社，1996：54-105．

[6] 董毓新，李彥碩．水電站建筑物結(jié)構(gòu)分析[M]．大連：大連理工大學(xué)出版社，1995：157．

[7] 中華人民共和國國家經(jīng)濟貿(mào)易委員會．DL 5073—2000水工建筑物抗震設(shè)計規(guī)范[S]．北京：中國電力出版社，2000:18．

[8] 李瓚，陳興華，鄭建波，等．混凝土拱壩設(shè)計[M]．北京：中國電力出版社，2000：563-569．

[9] 李彥碩．水電站建筑物的振動[M]．大連：大連理工大學(xué)出版社，1989：13-22．

[10]CRAIG.R R J.結(jié)構(gòu)動力學(xué)[M]．常嶺，李振邦,譯．北京：人民交通出版社，1981：259-319．