葉柏陽，周松松，孫益松

(1.江蘇淮源工程建設(shè)監(jiān)理有限公司， 江蘇 淮安 223005；2.淮安市水利勘測設(shè)計研究院有限公司，江蘇 淮安 223005)

水利工程建設(shè)為灌區(qū)水資源高效利用、工業(yè)經(jīng)濟發(fā)展以及人類生存均提供了重要支撐作用，但水利工程的安全穩(wěn)定是其發(fā)揮重大作用的前提保障，因而有許多水利工程師持續(xù)致力研究水工建筑的安全運營[1- 3]。除靜力荷載下結(jié)構(gòu)安全穩(wěn)定以外，而地震作為對水利設(shè)施產(chǎn)生較大威脅的自然災(zāi)害，揭示地震荷載下水工建筑結(jié)構(gòu)動力特性很有必要[4- 6]。國內(nèi)外較多學(xué)者基于室內(nèi)振動臺試驗，開展室內(nèi)水工模型試驗，研究了水工結(jié)構(gòu)在模擬地震荷載破壞過程中應(yīng)力變形表現(xiàn)[7- 9]。在現(xiàn)場以微震儀器等監(jiān)測設(shè)備開展水利工程運營過程中監(jiān)測，分析工程結(jié)構(gòu)動力響應(yīng)狀態(tài)，為工程動力穩(wěn)定性評價提供重要參考[10- 11]。數(shù)值仿真為一種較為高效研究手段，利于水利工程等復(fù)雜工況研究[12- 14]，本文將基于模態(tài)仿真分析理論，利用數(shù)值軟件，研究水工結(jié)構(gòu)地震動力特性，為結(jié)構(gòu)抗震設(shè)計提供重要參考。

1 結(jié)構(gòu)動力模態(tài)分析理論

針對進水塔結(jié)構(gòu)開展動力抗震分析主要分為2個方面，一方面為進水塔結(jié)構(gòu)自振特性，另一方面為地震動荷載作用下動力響應(yīng)特性，故而本文針對2個方面分別采用相適應(yīng)的模態(tài)分析計算理論。地震動荷載作用下結(jié)構(gòu)運動方程服從以下方程[15- 16]：

(1)

引入結(jié)構(gòu)地面加速度響應(yīng)值a，則式(1)可變換為：

(2)

在數(shù)值仿真體系中，以模態(tài)法假定地震動荷載作用方式與作用節(jié)點，其等效表達式為：

(3)

為求解地震動荷載下運動方程，仿真體系中引入單元節(jié)點協(xié)調(diào)質(zhì)量矩陣表達式為：

(4)

結(jié)構(gòu)運動過程中會受到阻尼影響，本文以Rayleigh表達式作為阻尼系數(shù)方程，如下所述：

[C]=α[M]+β[K]

(5)

式中，α、β—固有屬性參數(shù)系數(shù)值，其具體表達式為：

(6)

式中，ωi、ωj—自振頻率；ζi、ζj—阻尼系數(shù)參數(shù)。

以HHT快速變換作為數(shù)值仿真迭代求解手段[17]，進而求解獲得泵站進水塔地震動力響應(yīng)特征。另一方面，進水塔結(jié)構(gòu)自振特性以式(1)作為變化參量，運動方程表達式為：

(7)

引入外參數(shù)換算，進而求解式(7)，獲得：

([K]-ω2[M]){φ}=0

(8)

式中，{φ}—振幅。

結(jié)構(gòu)自振狀態(tài)下頻率參數(shù)可表述為：

ω1≤ω2≤ω3≤…≤ωn

(9)

以表達式綜合表述自振特性，進水塔自振頻率解集為：

(10)

從該解集方程式中，可知進水塔自振固有頻率變化特征，即ω1、ω2、ω3……等固有頻率的振動模態(tài)性。本文將利用上述模態(tài)分析計算理論，探討進水塔自振特性與塔群動力抗震響應(yīng)特征。

2 工程概況

華中地區(qū)某水利樞紐工程建設(shè)有抽水泵站，主要面向區(qū)域內(nèi)農(nóng)業(yè)用水、工業(yè)用水調(diào)控，在枯水期為農(nóng)業(yè)灌區(qū)生產(chǎn)提供重要過渡性水資源，上游蓄水庫正常運營期水位約為600m，總庫容超過5000萬m3，建設(shè)有總長約為150km的輸水干渠，渠首聯(lián)通泵站出水池，通閘流量約為0.7m3/s。樞紐工程還涉及混凝土面板堆石壩、溢洪道及下游消力池等水工設(shè)施，抽水泵站與堆石壩軸線平行，另外進水塔主要功能是使引水隧洞水流進入抽水泵站，其支撐排墩為5根混凝土立柱，直徑約為2.6m，高度為13m。另進水塔內(nèi)設(shè)置有弧型鋼閘門，以調(diào)控水流量，直徑約為1.6m，設(shè)置有液壓式啟閉機作為調(diào)控設(shè)施，精確控制閘門開度，提升水資源利用效率。進水塔頂高程為607m，實際塔高為25m，支撐排墩以橫向梁連接，連接構(gòu)件位于塔頂下6m處，厚度約為4.5m，間隔6m布設(shè)梁結(jié)構(gòu)。地區(qū)內(nèi)包括泵站在內(nèi)的所有水工建筑均采用Ⅰ級水利建筑物設(shè)計，但仍然不可忽視泵站工程在地震動荷載作用下其安全運營關(guān)乎地區(qū)水利安全，因而本文考慮對泵站所屬進水塔開展動力抗震分析。

工程場地內(nèi)僅在面板堆石壩左岸坡可見褶皺構(gòu)造，延伸長度為1.6km左右，破碎帶方向與水流方向相反。進水塔支撐排墩所在場地表面覆蓋土層為第四系人工活動填土層，農(nóng)業(yè)灌區(qū)內(nèi)生產(chǎn)種植土亦是該土層，顆粒松散性較大，含水量中等，約為32%，厚度亦較??；另在基巖上覆土層中另有一層砂礫土，粒徑為1.6～6.8mm，級配良好，主要為上游水流沖擊搬運沉積形成砂石土，目前作為渠道襯砌結(jié)構(gòu)所在持力層材料；基巖為弱風(fēng)化灰?guī)r，強度較高，室內(nèi)土工試驗測試表明孔隙度最低僅為0.1%，適合作為大型水工建筑結(jié)構(gòu)的承載層，水介質(zhì)滲透系數(shù)達10-19m2。泵站所用材料均為C25、C35素混凝土，距離工程現(xiàn)場所在地2.1km即為混凝土生產(chǎn)養(yǎng)護場，以泵車運輸及澆筑，保證原材料性能滿足工程要求。

3 進水塔結(jié)構(gòu)動力抗震分析

3.1 模型建立

按照進水塔平面設(shè)計圖資料，利用ANSYS建立數(shù)值分析模型，如圖1(a)所示，該模型包括進水塔與下方地基，此主要為了探討由地基傳輸?shù)卣饎幽芰恐辽喜拷Y(jié)構(gòu)過程中進水塔結(jié)構(gòu)響應(yīng)特征。在仿真計算體系中空間坐標(biāo)系x、y、z正方向分別確定為ANSYS數(shù)值計算中水流順向、結(jié)構(gòu)豎向、水流垂向，地基傳輸?shù)卣饎雍奢d界面采用零質(zhì)量地基模型。在前述參數(shù)設(shè)定基礎(chǔ)上，劃分模型網(wǎng)格單元數(shù)586246個，節(jié)點數(shù)48126個，劃分后進水塔結(jié)構(gòu)網(wǎng)格模型單元如圖1(b)所示。

圖1 數(shù)組模型

3.2 自振特性分析

基于模態(tài)分析理論計算獲得進水塔自振作用下各階次自振頻率，如圖2所示。從圖中可知，自振頻率值與計算階次為正相關(guān)遞增變化，第1階次自振頻率為2.297Hz，而第10階次相比第1階次增大了5.05倍；對比后10階次與前10階次之間自振頻率關(guān)系可知，第20階次相比第10階次自振頻率增幅為56.7%，表明自振頻率與計算階次遞增幅度逐漸放緩。

圖2 自振頻率與計算階次關(guān)系曲線

圖3為各計算階次下進水塔自振振型分布特征云圖，從圖中可知，在低階次下進水塔振型傾向于水流方向，如第1階次中z向參與系數(shù)達7.001，是y向參與系數(shù)的3個量級之多，表明低階次下自振振動以水流垂直方向為主；此現(xiàn)象在第3階次中亦是如此，水流垂直方向上的參與系數(shù)遠超其他2個方向，表明進水塔垂直水流方向上剛度與強度均應(yīng)重點關(guān)注。當(dāng)計算階次為第5階次時，此時仿真模態(tài)計算可認(rèn)為達到迭代中期階段，此時進水塔振動方向傾向于進水塔的豎直方向，此階態(tài)下y向參與系數(shù)達9.055，反而是z方向參與系數(shù)的3個量級之多，即中等計算階次下振型在y、z方向上發(fā)生逆轉(zhuǎn)性變化。當(dāng)達到第10階次時，振型分布中以x向參與系數(shù)最大，達1.933，其中y向為負(fù)參與系數(shù)，為-0.244，此時振型以順?biāo)鞣较驗橹?，但與此同時亦在發(fā)生向下沉降變形，表現(xiàn)為組合振型振動狀態(tài)。在第11階次以后，如第15階次、第20階次，振型傾向均為垂直水流方向，分析造成這種現(xiàn)象并不是由于順?biāo)鞣较蛏显诟哂嬎汶A次下抗震性能較差，而是進水塔排墩剛度相比進水塔整體要低，其設(shè)置的聯(lián)系梁剛度、強度、抗震性能均低于進水塔，故而會造成排墩在高計算階次下表現(xiàn)垂直水流方向上顯著振動，進而引起進水塔高計算階次的振型變化。綜上分析表明，隨計算階次遞增，進水塔振型從z向、y向逆轉(zhuǎn)性主導(dǎo)變化，再以x、y向組合振動主導(dǎo)，最后以低剛度的排墩導(dǎo)致了進水塔為z向主導(dǎo)。

3.3 動力抗震分析

為體現(xiàn)進水塔動力抗震特性，本文引入地震動荷載作用在泵站水工建筑上，從前述地震動模態(tài)分析理論中可知，地震動作用實質(zhì)上可分為2個部分，包括有慣性力荷載與動水荷載，動水荷載以附加質(zhì)量法表征，而慣性力荷載采用擬靜力法，各方向上時程曲線如圖4所示。

基于仿真模態(tài)體系啊求解計算獲得進水塔各方向位移最大值，并給出各方向最大位移所在節(jié)點處位移時程曲線，如圖5所示。從時程曲線可看出，x、y最大正向位移所處時間節(jié)點基本一致，均為8.5s左右，但z向位移最大值為10.5s左右，表明水流垂直方向上位移最大值應(yīng)是受排墩影響，進而滯后于x、y向。圖6為各方向位移最大值所對應(yīng)位移分布云圖，從圖6可看出，x正向最大位移為0.0271m，所處區(qū)域為進水塔頂部右側(cè)區(qū)域，且其負(fù)向最大位移亦位于塔頂，達0.0304m，方向上不論是正位移亦或是負(fù)位移，進水塔結(jié)構(gòu)中自塔頂至塔底部，均為逐漸降低態(tài)勢。y正方向最大位移為0.008m，處于進水塔迎水側(cè)，由迎水側(cè)至背水側(cè)，位移量值逐漸降低，y負(fù)方向最大位移為0.0172m，其塔身位移變化趨勢與正向一致，且正、負(fù)向最大位移時間節(jié)點對應(yīng)的位移具有對稱性。z正向最大位移為相比x、y正向分別增大了2.4、10.5倍，達0.0917m，表明垂直水流方向上柔度變形較大，抗震設(shè)計需要特別關(guān)注，z正、負(fù)向最大位移所處時間節(jié)點的塔身位移分布為對稱，負(fù)向最大位移為0.0712m，且各方向最大位移均出現(xiàn)在塔頂區(qū)域，表明塔頂受地震影響作用較強，剛度、強度均需加固。

圖3 進水塔自振振型分布特征云圖

圖4 地震動荷載時程曲線

圖5 各方向最大位移所在節(jié)點處位移時程曲線

圖6 各方向正負(fù)向最大位移節(jié)點時分布

圖7為進水塔典型部位代表方向應(yīng)力時程曲線，從圖中可看出，時程曲線應(yīng)力值基本以負(fù)向應(yīng)力為主，即壓應(yīng)力占據(jù)主導(dǎo)作用，典型部位中最大拉應(yīng)力為排墩聯(lián)系梁x方向，達13.5MPa，最大壓應(yīng)力亦為排墩聯(lián)系梁x方向，達12.4MPa，從材料允許應(yīng)力值來看，進水塔壓應(yīng)力處于安全狀態(tài)，但拉應(yīng)力顯著過高，會引起排墩聯(lián)系梁破壞。

圖8為各方向最大拉、壓應(yīng)力所處時間節(jié)點下進水塔應(yīng)力分布，從圖中可看出，各方向中最大拉應(yīng)力為x向，達17.3MPa，位于排墩縱向聯(lián)系梁，最大壓應(yīng)力亦是在該處，達15.2MPa。分析聯(lián)系梁應(yīng)力分布可知，縱向聯(lián)系梁應(yīng)力值在高程方向上逐漸降低，頂部拉應(yīng)力相比高程583m處聯(lián)系梁降低了48.6%；但排墩橫向聯(lián)系梁拉應(yīng)力為先增后減，峰值拉應(yīng)力為13.5MPa，位于高程590m處，即圖7(d)所示節(jié)點應(yīng)力時程曲線，同樣高程583m處拉應(yīng)力相比前者降低了11.1%。綜合地震動荷載作用下進水塔應(yīng)力特征可知，不論是橫向聯(lián)系梁亦或是縱向上，排墩聯(lián)系梁框架結(jié)構(gòu)為進水塔中受拉最顯著區(qū)域，極易發(fā)生破壞；另一方面，進水塔底板z方向拉應(yīng)力亦較大，最大拉應(yīng)力可達13.2MPa，筆者分析此主要受動水壓力影響而引起；為保證進水塔結(jié)構(gòu)抗震性能，應(yīng)對排墩聯(lián)系梁框架結(jié)構(gòu)、塔底板結(jié)構(gòu)加固剛度與強度，增強抗震穩(wěn)定性。

圖7 典型部位代表方向應(yīng)力時程曲線

圖8 各方向最大拉、壓應(yīng)力所處時間節(jié)點下應(yīng)力云圖

4 結(jié)論

(1)自振頻率與階次為正相關(guān)，但遞增幅度逐漸減小，10階次相比1階次增大了5.05倍，而20階次與10階次幅度為56.7%；振型在低階次與中等階次下為z、y交替主導(dǎo)，直至x、y組合主導(dǎo)，高階次下為z向主導(dǎo)。

(2)進水塔x、y最大正向位移所處時間節(jié)點基本一致，均位于8.5s，但z向滯后至10.5s，且z正向最大位移均高于x、y正向，最大正負(fù)向位移節(jié)點時塔體位移分布為對稱式。

(3)進水塔最大拉、壓應(yīng)力均處于x向縱向梁，該梁應(yīng)力在高程上遞減，但橫向梁拉應(yīng)力為先增后減，峰值拉應(yīng)力為13.5MPa；抗震設(shè)計應(yīng)著重對排墩聯(lián)系梁框架、塔底板結(jié)構(gòu)加固。