黨康寧+王媛+宋志強+劉云賀+張希

摘要：介紹了基于速度勢解析算法求解水電站進水塔在水中彎曲自由振動的方法。分析進水塔的主要受力特點，為引入流體速度勢解析算法對結(jié)構(gòu)進行簡化和假設(shè)。以分布參數(shù)梁體系推導(dǎo)得到無水時進水塔彎曲自由振動頻率方程，根據(jù)速度勢拉普拉斯方程確定了所求解問題的邊界條件，通過內(nèi)外水壓力表達式建立進水塔的彎曲自由振動方程。推導(dǎo)振動方程得到頻率方程，經(jīng)過數(shù)學(xué)軟件編程求解，最終可得到結(jié)構(gòu)振型和頻率。推導(dǎo)過程表明流體速度勢法是將水體轉(zhuǎn)化成與塔體振型相關(guān)的質(zhì)量附加于結(jié)構(gòu)，從而影響結(jié)構(gòu)的振動。與數(shù)值方法對比和探討，分析了兩者異同，提出本文方法適用范圍。結(jié)果表明，本文解法可以引入內(nèi)外含水矩形橫截面進水塔的彎曲自由振動問題的求解，為下一步進水塔地震響應(yīng)分析提供一種新的思路。

關(guān)鍵詞：進水塔；彎曲自由振動；速度勢法；解析方法；流固耦合

中圖分類號：TV312 文獻標(biāo)識碼：A 文章編號：1672-1683（2017）04-0142-06

Abstract：This paper introduces the velocity potential analysis method for solving the free bending vibration of an intake tower in water in a hydroelectric power station.It analyzes the basic mechanical characteristics of an intake tower and makes some simplification and hypothesis for the structure in order to introduce the velocity potential analysis method.The frequency equation of the free bending vibration of a standalone intake tower is derived from the distribution parameter beam system.The boundary conditions of the problem are determined according to the Laplace equation of the velocity potential.The bending vibration equation of the intake tower is established via the expression of internal and external water pressure.The frequency equation is derived from the vibration equation.The structural vibration mode and frequency can be obtained through mathematical software programming.The derivation shows that the fluid velocity potential method converts the water into added mass that is associated with the vibration mode of the tower and thus affects the vibration of the structure.This paper then compares this method with the numerical method，analyzes their similarities and differences，and proposes the scope of application of this method.Results show that this method can be used to solve the free bending vibration problem of intake towers with a rectangular cross-section and water on both the inside and the outside，and provides a new way of thinking for further analysis of the seismic response of intake towers.

Key words：iintake tower；free bending vibration；velocity potential method；analytic method；fluid-solid coupling

水電站進水塔為高聳的懸臂結(jié)構(gòu)，一般情況下，進水塔處于水體包圍中，內(nèi)外均有不超過結(jié)構(gòu)高度的流動水體。因高寬比較大結(jié)構(gòu)抗彎剛度相對較柔，在水庫中受到的風(fēng)、地震等橫向荷載作用[1]起主導(dǎo)作用。地震時結(jié)構(gòu)-流體耦聯(lián)情況下結(jié)構(gòu)的響應(yīng)一直為學(xué)者積極研究的問題，水體對結(jié)構(gòu)的影響機理及處理方法備受關(guān)注[2-3]。

對于水中結(jié)構(gòu)響應(yīng)的研究，從20世紀30年代開始學(xué)者便對水中柱體的彎曲自由振動進行研究，并取得一些成果[4-6]。速度勢解析方法以其理論基礎(chǔ)完善、求解結(jié)果精確、更能反映問題內(nèi)在機理等特點，很早便在固體和流體耦合領(lǐng)域采用。國內(nèi)學(xué)者在20世紀80年代、90年代對水中結(jié)構(gòu)的自由振動問題進行了大量的研究：張悉德對部分埋入水中的懸臂圓柱體和橢圓柱體的彎曲自由振動進行了研究[7，8]，首次給出柱水耦聯(lián)體系振型函數(shù)的精確解和以有限階行列式表示的頻率方程。周叮、朱永誼、陳怡然等學(xué)者[9-14]研究了貯液圓筒、橢圓筒、圓環(huán)形貯液容器等多種截面形狀的水中結(jié)構(gòu)自由振動問題，亦考慮表面波，不同深度等因素的影響，導(dǎo)出了貯液圓筒的振型函數(shù)及固有頻率的精確計算公式。以上研究著眼于解決工程問題，以水中實心或空心圓柱、橢圓柱結(jié)構(gòu)為主，采用速度勢解析方法，結(jié)構(gòu)沒于水中或高于液面，內(nèi)部有水或僅外部有水等情況求解其水中振型函數(shù)及固有頻率。對于水中進水塔的振動問題，采用解析算法的研究較少，董育堅[15]對拉普拉斯表示的動水壓力用數(shù)值法求解該方程，并采用有限元理論詳細推導(dǎo)了矩形截面進水塔的地震動水壓力。

進水塔為內(nèi)外充水，部分淹沒于水中的高聳懸臂式結(jié)構(gòu)，其邊界條件較為簡單，與上述研究結(jié)構(gòu)既有很多相似之處也有不同：進水塔橫截面為矩形，中空薄壁，振動時內(nèi)外動水壓力應(yīng)共同對結(jié)構(gòu)產(chǎn)生影響。因此，如嘗試將流固耦合的解析算法推廣到此類結(jié)構(gòu)上具有重要意義。本文根據(jù)進水塔的主要受力特點，將其等效簡化為受彎曲振動的梁結(jié)構(gòu)，將速度勢解析方法引入求解塔體-水體耦聯(lián)彎曲振動，為其后續(xù)地震作用下結(jié)構(gòu)響應(yīng)提供基礎(chǔ)，也為研究類似水中中空矩形截面結(jié)構(gòu)提供一種新的思路和方法。

1 基本模型

1.1 塔水耦合模型及簡化假設(shè)

文獻[16]通過對進水塔動水壓力的研究表明地震荷載作用下，順?biāo)鞣较騽铀畨毫h大于水體重力對進水塔的影響，且慣性耦合作用主要存在與水流水方向。本文研究對象為考慮平面內(nèi)（沿X軸向）彎曲的進水塔，圖1示意了塔-水耦合體系。

3 本文方法結(jié)果分析和討論

帶入具體數(shù)據(jù)后便能夠得到計算結(jié)果，為分析本文解析方法適用范圍，同時使用有限元軟件進行數(shù)值計算比對分析。以某實際工程的洪洞進水塔為依據(jù)，提取其主要特征尺寸，得到各具體參數(shù)如下（符號代表的含義見圖1）：進水塔塔高度H=81 m，塔體橫截面尺寸：a=2.34 m，b=6.78 m，c=7 m，d=14 m；塔體材料彈性模量E=2.8×1010 Pa，塔體密度ρ0=2 500 kg/m3。塔內(nèi)外水位高度相同均為h=75 m，水的密度為ρ=1 000 kg/m3，水的體積彈性模量取一個大氣壓下20攝氏度時值，為2.18×109 Pa。分別建立無水進水塔模型（見圖2（a））和內(nèi)外有水模型（見圖2（b）），進水塔塔體采用三維實體單元，水體采用三維流體單元，程序會在結(jié)構(gòu)和流體接觸面自動設(shè)置耦合邊界。進水塔底部完全固結(jié)，塔體周圍的流體范圍為四倍塔體寬度或長度的尺寸，內(nèi)外流體表面均設(shè)置自由表面。

圖3為塔水耦合有限元模型進水塔在XZ平面彎曲前三階振型和頻率值。圖4為速度勢法求得的進水塔前三階彎曲振型。從有限元分析的結(jié)果來看，本文公式推導(dǎo)得到的振型均能在有限元孤塔和塔水耦合模型中塔的找到對應(yīng)振型，兩者振型相似。速度勢法和有限元法得到的孤塔模型振型頻率都大于塔水耦合模型，1階降低分別為9.78%，13.31%，可見水的存在對進水塔振動頻率有顯著影響，計算中不可忽略水的作用。

由表1可以看出，有限元分析得到的塔體各階對應(yīng)頻率均小于速度勢法求得值，并且隨著階數(shù)的增大，兩者差別越來越大。其中1階頻率兩種方法差異在5%左右，可認為速度勢法基本能得到滿意的基本振型。分析原因，兩種方法具有以下不同：速度勢法將結(jié)構(gòu)簡化成平面內(nèi)梁單元模型，其實質(zhì)是采用歐拉梁理論，假設(shè)梁只有彎曲變形，而梁對剪切變形完全剛性，即忽略了剪切變形和轉(zhuǎn)動慣量，這對于梁截面尺寸遠小于長度的梁（高跨比<1/10）來說，不會引起顯著的誤差。以上例子中，梁的長度與高度之比為H/2b=5.97，此比值的梁屬于深受彎梁，此時梁內(nèi)的橫向剪切力所產(chǎn)生的剪切變形將引起梁的附加撓度，并使原來垂直于中性面的截面變形后不再與中性面垂直，與假定有所不同。同時，速度勢法將水體處理成與振型函數(shù)有關(guān)的附加質(zhì)量，改變梁沿高度的質(zhì)量分布，水體和塔體的相互作用未能完全考慮，有限元塔水耦合模型中，結(jié)構(gòu)和流體的相互作用可以體現(xiàn)。

為驗證高寬比對結(jié)果的影響，在前述實例的基礎(chǔ)上，將塔高度提高到135.6 m（10倍截面寬度）和271.2 m（20倍截面寬度），水位與前述例子的水位比例相當(dāng)，分別為125 m和250 m。其它參數(shù)不變。表2和表3列出分別列出了有限元法和速度勢法求的前三階振動頻率。

由表2、表3可以看出，當(dāng)塔高增加到截面寬度10倍時，速度勢法和有限元法計算的頻率更加接近，1階頻率相差0.8%。二階和三階相差在20%以內(nèi)。當(dāng)增加到截面寬度20倍時，二階和三階頻率相差不到10%。

由以上分析可以看出，結(jié)構(gòu)的高寬比大于10時，計算結(jié)果較準(zhǔn)確。有限元法中模型節(jié)點存在平動、扭轉(zhuǎn)等自由度，截面可變形，在高階振型時其它方向質(zhì)量參與復(fù)雜，這是兩種方法高階振型差異較多的原因。

4 結(jié)論

本文以進水塔結(jié)構(gòu)為研究對象，抓住其主要受力特點和邊界條件，將其簡化成內(nèi)外充水的懸臂梁結(jié)構(gòu)，在前人研究基礎(chǔ)上將圓柱型式水中結(jié)構(gòu)的彎曲自由振動算法引用到內(nèi)空矩形截面的進水塔結(jié)構(gòu)上，經(jīng)過嚴密的推導(dǎo)求得了進水塔的振動頻率和振型，可以得到如下結(jié)論。

（1）進水塔可簡化為懸臂梁結(jié)構(gòu)，采用流體速度勢法求解其水中彎曲自由振動。速度勢法雖然簡化假設(shè)較多，但不失為研究類似結(jié)構(gòu)提供一種新的思路。經(jīng)過實例分析和探討，可知與有限元法相比，速度勢法求解得到的1階頻率接近，高寬比大于10時，二三階頻率開始接近。但實際進水塔高寬比多在10以內(nèi)，速度勢法求得的1階頻率可以作為參考。

（2）推導(dǎo)得到了水-進水塔耦聯(lián)體系振型函數(shù)精確解和頻率方程，其表達式較為復(fù)雜，編制數(shù)學(xué)求解程序能夠快速求解得到結(jié)果，相比于有限元法省去復(fù)雜建模、計算和處理過程，計算過程明確，結(jié)果有一定參考價值。

（3）從推導(dǎo)過程可以看出，速度勢法求解結(jié)構(gòu)振動時將水的動水壓力作用轉(zhuǎn)化成與振型有關(guān)的質(zhì)量附加于結(jié)構(gòu)之上，增大了結(jié)構(gòu)的分布質(zhì)量，引起了自振頻率的降低，此方法為水體對結(jié)構(gòu)的單向流固耦合。

參考文獻（References）：

[1] 林繼庸.水工建筑物（第5版）[M].北京：中國水利電力出版社，2013.（LIN Ji-yong.Hydraulic structures （Fifth Edition）[M].Beijing：China Water & Power Press，2013 （in Chinese））

[2] Housner，G.W.1957.Dynamic Pressures on Accelerated Fluid Containers[J].Bulletin of the Seismological Society of America，1957，Vol 47，pp 15-35.

[3] 鄭哲敏，馬宗魁.懸臂梁在一側(cè)受有液體作用時的自由振動[J].力學(xué)學(xué)報，1959，3（2）：111-119.（ZHENG Zhe-min，MA Zong-kui.Vibration of cantilever beam placed against water with free surface[J].Acta Mechanica Sinica[J].1959，3（2）：111-119.（in Chinese））

[4] Housner G M.The dynamic behavior of water tanks[J].Bull.Seism.Soc.Am.1963.53（2）：381-387.

[5] Moiseev N N.Introduction to the theory of oscillations of liquid-containing bodies[J].Advances in Applied Mechanics，1964，8：233-289.DOI：10.1007/BF02128917

[6] Khalak A，Williamson C H K.Dynamics of a hydroelastic cylinder with very low mass and damping[J].Journal of Fluids & Structures，1996，10（5）：455–472.DOI：10.1006/jfls.1996.0031

[7] 張悉德.部分埋入水中懸臂圓柱體的彎曲自由振動[J].應(yīng)用數(shù)學(xué)和力學(xué)，1982.7，3（4）：537-546.（ZHANG Xi-de.Free Bending vibration of circular column partially submerged in water[J].Applied Mathematics and Mechanics；1982.7，3（4）：537-546.（in Chinese））

[8] 張悉德.部分埋入水中橢圓柱體彎曲自由振動[J].應(yīng)用數(shù)學(xué)和力學(xué)，1986.9，7（9）：813-820.（ZHANG Xi-de.Free bending vibration of elliptic cylinder partially submerged in water[J].Applied Mathematics and Mechanics；1986.9，7（9）：813-820.（in Chinese））

[9] 周叮.貯液橢圓筒在水中的橫向彎曲自由振動[J].應(yīng)用力學(xué)學(xué)報，1989，6（4）：83-90.（ZHOU Ding.Free transverse vibrations of elliptic reservoir submerged in water[J].Chinese Journal of Applied Mechanics，1989，6（4）：83-90.（in Chinese））

[10] 周叮.貯液圓筒在水中的彎曲自由振動[J].應(yīng)用力學(xué)學(xué)報，1990，11（5）：439-446.（HOU Ding.The free bending vibration of cylindrical tank partially filled with liquid and submerged in water[J].Chinese Journal of Applied Mechanics，1990，11 （5）：439-446.（in Chinese））

[11] 朱永誼，翁智遠.考慮液面晃動和液體可壓縮影[JP+2]響時的柱水耦聯(lián)振動[J].力學(xué)學(xué)報，1989，6 （6）：657-667.（ZHU Yong-yi，WEN Zhi-yuan.The coupled vibration between column and waterconsidering the effects of surface waveand compressibilty of water[J].Acta Mechanica Sinica，1989，6 （6）：657-667.（in Chinese））DOI：10.6052/0459-1879-1989-6-1989-355

[12] 周叮.考慮表面波的貯液圓筒的彎曲自由振動[J].振動與沖擊，1989（3）：32-36.（ZHOU Ding.Free bending vibration of cylindrical tank partially filled with liquid in the consideration of surface wave[J].Journal of Vibration and Shock，1989（3）：32-36.（in Chinese））

[13] 周叮.計及表面波的圓環(huán)形貯液容器的彎曲自由振動[J].應(yīng)用力學(xué)學(xué)報，1994，15（9）：789-796.（ZHOU Ding.Free bending vibration of annular cylindrical tank partially filled with liquid in the consideration of surface wave[J].Chinese Journal of Applied Mechanics，1994，15 （9）：789-796.（in Chinese））

[14] 陳怡然，白燁.考慮表面波部分潛入水中橢圓柱體的扭振分析[J].河北工業(yè)大學(xué)學(xué)報，1992，4（4）：78-83.（CHEN Yiran，BAI Ye.The Analysis of Torsional Vibration of Elliptic Column partially Submerged in Water with Surface Waves[J].Journal of Hebei University of Technology，1992，4（4）：78-83.（in Chinese））

[15] 董育堅.矩形截面進水塔的地震動水壓力計算[J].水利學(xué)報，1982（1）：47-52（DONG Yujian.Calculation of seismic hydrodynamic pressure in the inlet of the rectangular cross section[J].Journal of hydraulic engineering，1982（1）：47-52 （in Chinese））

[16] 黃虎，李異，張建偉 基于流固耦合的高聳進水塔動水壓力分布研究[J].水力發(fā)電，2012，38（6），30-33（HUANG Hu，LI Yi，ZHANG Jian-wei.Study on hydrodynamic pressure distribution of high intake tower based on fluid-structure interaction[J].Water Power，2012，38（6），30-33 （in Chinese）） DOI：10.3969/j.issn.0559-9342.2012.06.009

[17] （美）R.克拉夫，J.彭津.結(jié)構(gòu)動力學(xué)[M].第二版（修訂版）.北京：高等教育出版社，2006.（R.W.Clough and J.Penzien，Dynamics of Structures[M].Second Edition （Revised Edition）.Beijing：Higher Education Press，2006.（in Chinese））

[18] 劉晶波，杜修力.結(jié)構(gòu)動力學(xué)[M].北京：機械工業(yè)出版社，2005.1.（LIU Jing-bo，DU Xiu-li.Dynamics of structures[M].Beijing：China Machine Press，2005.1 （in Chinese））

[19] 譚曉晶.大型儲液罐流固耦合地震反應(yīng)分析[D].中國地震局工程力學(xué)研究所，2007.（TAN Xiao-jing.Seismic response analysis of large liquid storage tank based fluid solid coupling[D].Institute of Engineering Mechanics，China Earthquake Administration，2007 （in Chinese））

[20] 孔瓏.工程流體力學(xué)（第四版）[M].北京：中國電力出 版社，2014.2（KONG Long.The engineering fluid mechanics （Fourth Edition）[M].Beijing：Chinese Power Press，2014.2 （in Chinese）