柳樹搖 李紅 李楨

摘 要：進(jìn)水塔結(jié)構(gòu)復(fù)雜，在高地震烈度區(qū)表現(xiàn)為復(fù)雜的動(dòng)力特性，塔體結(jié)構(gòu)的穩(wěn)定性直接關(guān)系到結(jié)構(gòu)的安全性及設(shè)計(jì)的合理性。文章采用三維動(dòng)力有限元的方法，通過在人工邊界上輸入實(shí)測地震波來模擬地震作用，得到了進(jìn)水塔結(jié)構(gòu)動(dòng)位移和動(dòng)應(yīng)力響應(yīng)，并基于有限元計(jì)算結(jié)果進(jìn)行了進(jìn)水塔穩(wěn)定性分析。計(jì)算結(jié)果表明，拉拉山水電站進(jìn)水塔在地震作用下是穩(wěn)定的。

關(guān)鍵詞：進(jìn)水塔；有限元；動(dòng)力分析

1 工程概況

拉拉山水電站位于甘孜州巴塘縣境內(nèi)，電站為引水式開發(fā)，該電站位于巴楚河中游河段上，工程河段上起松多，下止措松龍，河段長23.5km，天然落差225m。電站閘址位于松多鄉(xiāng)松多村，正常蓄水位3004.00m，相應(yīng)庫容106萬m3，電站裝機(jī)容量96MW，具有日調(diào)節(jié)性能。電站進(jìn)水塔緊靠排污閘布置于左岸，進(jìn)口分水角為64°51'59"。進(jìn)水口由喇叭口、攔污柵、進(jìn)水閘組成，總長36.2m。

2 巖體及混凝土靜、動(dòng)力學(xué)參數(shù)

3 進(jìn)水塔抗震分析及評價(jià)

3.1 計(jì)算方法及假定

利用ANSYS有限元分析程序建立三維有限元模型，地震作用采用時(shí)程分析法。本次分析計(jì)算中所利用的時(shí)間歷程分析理論的基本假定有：（1）庫水假定為不可壓縮流體，因此庫水對進(jìn)水塔的動(dòng)力作用即相當(dāng)于附加質(zhì)量。（2）地基為均勻體，并人為地截取一條邊界作為地震輸入邊界。在計(jì)算中假定地震運(yùn)動(dòng)沿人工邊界均勻輸入。（3）假定人工邊界范圍以內(nèi)的地基是“無質(zhì)量地基”，在形成整個(gè)系統(tǒng)的特征矩陣時(shí)，地基單元只考慮彈性，不考慮質(zhì)量，以消除波的傳播效應(yīng)，避免人為的放大作用。

3.2 地震波的選取

地震波是一個(gè)頻帶較寬的非平穩(wěn)隨機(jī)振動(dòng)，受各種因素影響而變化，采用時(shí)程法對結(jié)構(gòu)進(jìn)行地震分析時(shí)，輸入地震波不同，所得的地震反應(yīng)相差較大，因此，合理選擇地震波進(jìn)行直接動(dòng)力分析是保證計(jì)算可靠性的重要保證。

拉拉山水電站的場地類別屬于III類，工程場地的地震基本烈度為Ⅷ度，進(jìn)水塔未來50年超越概率為10%基巖地震動(dòng)峰值加速度1.85m/s2，特征周期為0.4s，進(jìn)水塔抗震分析計(jì)算中選擇場地實(shí)測波進(jìn)行計(jì)算，本次計(jì)算考慮水平向地震波。地震波正向峰值加速度1.60m/s2，負(fù)向峰值加速度為-1.85m/s2。加速度時(shí)程曲線見圖1。

3.3 計(jì)算模型

計(jì)算模型模擬范圍包括塔體部分和地基部分，其中塔體部分的下部邊界取至底板下表面，高程2988.00m；平臺高程3005.50m，排架頂部高程3016.00m，計(jì)算模擬的總高度為28m。地基部分的上下游長約22m，左右兩側(cè)長約19.5m，與塔體兩側(cè)一致；地基深度取25m，即延伸約1.5倍最大塔高。坐標(biāo)系的X軸取順?biāo)鞣较颍灾赶蛏嫌螢檎?；Y軸取垂直水流方向，指向右岸為正；Z軸取豎直方向，向上為正。符合右手螺旋定則。塔體及地基的三維有限元模型見圖2。

3.4 計(jì)算結(jié)果及分析

3.4.1 位移場分析

經(jīng)比較，在t=13.22s時(shí)，進(jìn)水塔結(jié)構(gòu)綜合位移值較大。從整個(gè)時(shí)間歷程上看，順?biāo)飨蛭灰芔x：對于塔體頂部，X向動(dòng)位移最大值為-0.003mm，位于塔體頂部，出現(xiàn)時(shí)間為t=13.22s；X向動(dòng)位移最小值為-0.01mm，出現(xiàn)時(shí)間為t=5.46s；對于排架，X向動(dòng)位移最大值為2.84mm，位于啟閉機(jī)房頂部，出現(xiàn)時(shí)間為t=13.22s；X向動(dòng)位移最小值-2.29mm，位于啟閉機(jī)房頂部處，出現(xiàn)時(shí)間為t=11.46s。豎向位移Uz：對于塔體頂部，Z向動(dòng)位移最大響應(yīng)為-0.162mm，出現(xiàn)時(shí)間為t=13.22s；Z向動(dòng)位移最小響應(yīng)為-0.163mm，出現(xiàn)時(shí)間為t=11.38s；對于排架，Z向動(dòng)位移最大響應(yīng)為-0.33mm，出現(xiàn)時(shí)間為t=11.46s；Z向動(dòng)位移最小響應(yīng)-0.55mm，出現(xiàn)時(shí)間為t=13.22s。

從時(shí)間點(diǎn)（塔體t=13.22s）的位移等值線圖上看，順?biāo)飨蛭灰芔x：對于塔體和排架整體，X向動(dòng)位移最大響應(yīng)為2.92mm，位于排架頂部；上部排架X向動(dòng)位移明顯大于下部塔體。垂直水流向位移Uy：對于塔體和排架整體，Y向動(dòng)位移最大響應(yīng)為0.23mm，位于明渠平臺邊緣；除去排架頂部，Y向動(dòng)位移在-0.14mm～0.21mm范圍內(nèi)。豎向位移Uz：對于塔體和排架整體，Z向動(dòng)位移最大響應(yīng)為0.65mm，位于排架頂?shù)闹胁?；Z向動(dòng)位移分布趨勢為從上而下逐漸減小。綜合位移Usum：對于塔體和排架整體，綜合位移最大響應(yīng)為2.92mm，位于塔體頂?shù)闹胁浚痪C合位移以豎向位移為主，分布趨勢同豎向位移相似。進(jìn)水塔排架啟閉機(jī)房頂部動(dòng)位移最大，分布趨勢從上而下逐漸減小，整體動(dòng)位移較小。

3.4.2 應(yīng)力場分析

經(jīng)比較，在t=11.87s時(shí)，進(jìn)水塔結(jié)構(gòu)動(dòng)應(yīng)力響應(yīng)較大。從代表點(diǎn)的整個(gè)時(shí)間歷程上看，在t=11.87s時(shí)，進(jìn)水塔動(dòng)應(yīng)力響應(yīng)最大，此時(shí)整個(gè)結(jié)構(gòu)以受壓為主；各應(yīng)力關(guān)鍵點(diǎn)中σ1最大值為1.15MPa，位于進(jìn)水口胸墻處d點(diǎn)，出現(xiàn)時(shí)間為t=11.87s；σ3最小值為-1.54MPa，位于排架底部e點(diǎn)，出現(xiàn)時(shí)間為t=7.54s。地震作用效應(yīng)的計(jì)算結(jié)果表明，在塔體與排架結(jié)構(gòu)結(jié)合部位、中間各層板梁結(jié)構(gòu)、部分載面突變處會(huì)產(chǎn)生較高的應(yīng)力集中，需要配筋處理，增強(qiáng)混凝土的抗拉性能。

從時(shí)間點(diǎn)（塔體t=11.87s）的應(yīng)力等值線圖上看，第一主應(yīng)力σ1：在地震作用下，第一主應(yīng)力主要為拉應(yīng)力，最大值為1.4MPa，為應(yīng)力集中，未超過C25混凝土動(dòng)態(tài)抗拉強(qiáng)度2.31MPa；除去應(yīng)力集中區(qū)域，塔架大部分區(qū)域應(yīng)力在1.0MPa以內(nèi)；塔體第一主應(yīng)力值較小。第三主應(yīng)力σ3：在地震作用下，第三主應(yīng)力主要為壓應(yīng)力，最大壓應(yīng)力為-5.52MPa，發(fā)生在進(jìn)水塔胸墻部位；較大壓應(yīng)力主要出現(xiàn)在上部排架柱的底部和各梁柱的連接處，除去應(yīng)力集中區(qū)域，進(jìn)水塔整體第三主應(yīng)力值較小，大部分區(qū)域壓應(yīng)力值在-3.5MPa以內(nèi)，未超過C25混凝土動(dòng)態(tài)抗壓強(qiáng)度21.71MPa。Z向正應(yīng)力σz：塔體沿Z向以受壓為主，壓應(yīng)力最小值為-3.1MPa；較大壓應(yīng)力主要發(fā)生在排架柱底部；塔體Z向應(yīng)力值較小，大部分區(qū)域應(yīng)力值在-2.3MPa～-0.1MPa范圍內(nèi)。巖基豎直向正應(yīng)力σz：塔體底部基礎(chǔ)豎直向全部承受壓應(yīng)力，最大豎直向壓應(yīng)力-0.42MPa，發(fā)生塔體兩側(cè)明渠底部巖基，進(jìn)水塔巖基動(dòng)態(tài)承載力標(biāo)準(zhǔn)值為靜態(tài)承載力標(biāo)準(zhǔn)值的1.5倍，則動(dòng)態(tài)承載力標(biāo)準(zhǔn)值為1.5MPa～1.8MPa，巖基最大豎向壓應(yīng)力值小于巖基動(dòng)態(tài)承載力標(biāo)準(zhǔn)值，滿足規(guī)范要求。

3.4.3 基于有限元計(jì)算結(jié)果的進(jìn)水塔穩(wěn)定計(jì)算

基于前面有限元?jiǎng)恿Φ挠?jì)算結(jié)果，對進(jìn)水塔沿建基面的抗滑穩(wěn)定進(jìn)行計(jì)算分析。在地震波的作用下，進(jìn)水塔沿建基面的抗滑穩(wěn)定安全系數(shù)為6.20～8.60，均大于規(guī)范規(guī)定的抗滑穩(wěn)定安全系數(shù)1.05，滿足規(guī)范要求，說明進(jìn)水塔在地震工況作用下是安全穩(wěn)定的。

3.5 進(jìn)水塔抗震時(shí)程分析及評價(jià)

進(jìn)水塔頂部動(dòng)位移最大，分布趨勢從上而下逐漸減小，整體動(dòng)位移不大，規(guī)律性合理；根據(jù)前面的有限元?jiǎng)恿τ?jì)算結(jié)果，進(jìn)水塔動(dòng)應(yīng)力均小于混凝土強(qiáng)度設(shè)計(jì)值，基礎(chǔ)豎向應(yīng)力均小于基礎(chǔ)承載力；采用基于有限元靜力計(jì)算成果的剛體極限平衡法，將進(jìn)水塔在實(shí)測波作用下，沿建基面的抗滑安全穩(wěn)定系數(shù)大于規(guī)范允許安全系數(shù)，進(jìn)水塔沿建基面安全穩(wěn)定。

4 結(jié)束語

文章采用三維有限元方法，對拉拉山水電站進(jìn)水塔進(jìn)行了三維動(dòng)力分析，分析了其位移和應(yīng)力響應(yīng)，對分析了其抗滑穩(wěn)定安全性。分析結(jié)果表明，進(jìn)水塔在地震工況下，其結(jié)構(gòu)應(yīng)力、基礎(chǔ)承載力及抗滑穩(wěn)定均滿足規(guī)范要求。為確保工程安全，建議加強(qiáng)地震動(dòng)監(jiān)測，密切關(guān)注震前、震后進(jìn)水塔結(jié)構(gòu)的位移和應(yīng)力變化。同時(shí)，進(jìn)水塔塔身結(jié)構(gòu)體型力求簡單對稱，質(zhì)量和剛度變化平緩，減少應(yīng)力集中，并確保足夠剛度。塔體上部和后部與巖體緊貼，并在基礎(chǔ)部位布設(shè)固結(jié)灌漿，增加整體性。

參考文獻(xiàn)

[1]SL203-97.水工建筑物抗震設(shè)計(jì)規(guī)范[S].北京：水利電力出版社，1998.

[2]樂成軍，任旭華，邵勇，等.水電站進(jìn)水塔結(jié)構(gòu)抗震設(shè)計(jì)研究[J].水力發(fā)電，2009，35（5）：86-89.

[3]范書立，陳健云，陳明陽，等.水電站進(jìn)水塔不同抗震驗(yàn)算方法對比分析[J].人民長江，2011，42（17）：51-54.

[4]趙曉西，王宗敏，朱太山，等.ANSYS軟件在高進(jìn)水塔靜動(dòng)力分析中的應(yīng)用[J].計(jì)算機(jī)應(yīng)用，2005，31（3）：72-73.