柳樹搖 李紅 李楨

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

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

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 巖體及混凝土靜、動力學(xué)參數(shù)

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

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

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

3.2 地震波的選取

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

拉拉山水電站的場地類別屬于III類，工程場地的地震基本烈度為Ⅷ度，進(jìn)水塔未來50年超越概率為10%基巖地震動峰值加速度1.85m/s2，特征周期為0.4s，進(jìn)水塔抗震分析計(jì)算中選擇場地實(shí)測波進(jìn)行計(jì)算，本次計(jì)算考慮水平向地震波。地震波正向峰值加速度1.60m/s2，負(fù)向峰值加速度為-1.85m/s2。加速度時程曲線見圖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時，進(jìn)水塔結(jié)構(gòu)綜合位移值較大。從整個時間歷程上看，順?biāo)飨蛭灰芔x：對于塔體頂部，X向動位移最大值為-0.003mm，位于塔體頂部，出現(xiàn)時間為t=13.22s；X向動位移最小值為-0.01mm，出現(xiàn)時間為t=5.46s；對于排架，X向動位移最大值為2.84mm，位于啟閉機(jī)房頂部，出現(xiàn)時間為t=13.22s；X向動位移最小值-2.29mm，位于啟閉機(jī)房頂部處，出現(xiàn)時間為t=11.46s。豎向位移Uz：對于塔體頂部，Z向動位移最大響應(yīng)為-0.162mm，出現(xiàn)時間為t=13.22s；Z向動位移最小響應(yīng)為-0.163mm，出現(xiàn)時間為t=11.38s；對于排架，Z向動位移最大響應(yīng)為-0.33mm，出現(xiàn)時間為t=11.46s；Z向動位移最小響應(yīng)-0.55mm，出現(xiàn)時間為t=13.22s。

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

3.4.2 應(yīng)力場分析

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

從時間點(diǎn)（塔體t=11.87s）的應(yīng)力等值線圖上看，第一主應(yīng)力σ1：在地震作用下，第一主應(yīng)力主要為拉應(yīng)力，最大值為1.4MPa，為應(yīng)力集中，未超過C25混凝土動態(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混凝土動態(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)水塔巖基動態(tài)承載力標(biāo)準(zhǔn)值為靜態(tài)承載力標(biāo)準(zhǔn)值的1.5倍，則動態(tài)承載力標(biāo)準(zhǔn)值為1.5MPa～1.8MPa，巖基最大豎向壓應(yīng)力值小于巖基動態(tài)承載力標(biāo)準(zhǔn)值，滿足規(guī)范要求。

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

基于前面有限元動力的計(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)水塔抗震時程分析及評價

進(jìn)水塔頂部動位移最大，分布趨勢從上而下逐漸減小，整體動位移不大，規(guī)律性合理；根據(jù)前面的有限元動力計(jì)算結(jié)果，進(jìn)水塔動應(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)行了三維動力分析，分析了其位移和應(yīng)力響應(yīng)，對分析了其抗滑穩(wěn)定安全性。分析結(jié)果表明，進(jìn)水塔在地震工況下，其結(jié)構(gòu)應(yīng)力、基礎(chǔ)承載力及抗滑穩(wěn)定均滿足規(guī)范要求。為確保工程安全，建議加強(qiáng)地震動監(jiān)測，密切關(guān)注震前、震后進(jìn)水塔結(jié)構(gòu)的位移和應(yīng)力變化。同時，進(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)水塔靜動力分析中的應(yīng)用[J].計(jì)算機(jī)應(yīng)用，2005，31（3）：72-73.