鄧育林 郭慶康 何雄君

(武漢理工大學(xué)交通學(xué)院 武漢 430063)

考慮流固耦合作用的橋梁深水矩形空心高墩振動特性分析*

鄧育林 郭慶康 何雄君

(武漢理工大學(xué)交通學(xué)院 武漢 430063)

以一典型深水矩形空心高墩為例，采用完全數(shù)值法，基于勢流理論建立了橋墩-水流固耦合數(shù)值分析模型，分析了墩身僅外有水、僅內(nèi)有水，以及內(nèi)、外均有水等3種情況下橋墩振動特性的變化規(guī)律.結(jié)果表明，墩-水耦合作用使得橋墩各階振型周期有所增大，且墩身內(nèi)、外均有水時增幅最大；當(dāng)水深小于墩高一半時，橋墩低階振型周期變化很小，但高階振型周期變化較大；上部結(jié)構(gòu)質(zhì)量雖然會減弱流固耦合作用對橋墩低階振型的影響，但會增大對橋墩高階振型的影響.由于高階振型對高墩底部剪力響應(yīng)貢獻較大,因此，水深不大時應(yīng)當(dāng)重視橋墩高階振型周期的變化.

橋梁；深水橋墩；矩形空心高墩；流固耦合；振動特性

0 引 言

大量研究表明，水的存在會改變水下結(jié)構(gòu)的動力特性，如周期、振型、阻尼比等參數(shù)，進而影響結(jié)構(gòu)在地震等作用下的動態(tài)響應(yīng).對于跨江、河、水庫等深水橋梁，由于這些橋梁的墩身和基礎(chǔ)大多位于深水中，不可避免的要受周圍水體的影響，因此，深入了解墩-水耦合作用下對結(jié)構(gòu)振動特性的影響規(guī)律，是研究其在如地震荷載等動態(tài)荷載作用的響應(yīng)規(guī)律的前提與基礎(chǔ).

雖然水-結(jié)構(gòu)動力相互作用的研究已有80多年的歷史，但研究主要針對儲液罐和壩體等結(jié)構(gòu)，而有關(guān)橋墩-水動力相互作用的研究相對較少.主要的研究方法有解析法、半解析半數(shù)值法和完全數(shù)值法等3類.其中，解析法、半解析半數(shù)值法均需先得到動水附加質(zhì)量和附加阻尼的解析解，然后集中到結(jié)構(gòu)上，雖然可以降低計算的復(fù)雜程度，但是解的準確性受到近似及簡化的限制[1-2]，而完全數(shù)值法是對結(jié)構(gòu)和水體進行精確模擬和分析有效數(shù)值手段，可用來分析各類復(fù)雜的流固耦合問題，如楊吉新等[3-7]采該方法取得較好的研究成果.

但目前分析大多針對實體橋墩或群樁基礎(chǔ)，而針對空心橋墩，特別是空心墩內(nèi)水的影響研究較少，為此，文中以典型深水橋梁矩形空心高墩為例，利用有限元分析軟件ADINA，基于勢流體單元的流固耦合完全數(shù)值法，建立了深水矩形空心高墩-內(nèi)、外水流固耦合數(shù)值分析模型，分析了高墩結(jié)構(gòu)在墩身僅外有水、僅內(nèi)有水，以及內(nèi)、外均有水3種情況下橋墩振動特性的變化規(guī)律，特別是高階振型的變化，并對墩頂附加質(zhì)量的影響進行了討論.

1 水-結(jié)構(gòu)耦合系統(tǒng)模態(tài)分析原理

水下結(jié)構(gòu)模態(tài)數(shù)值分析可采用基于勢的流體單元來模擬結(jié)構(gòu)-水的流固耦合作用.根據(jù)流體力學(xué)基本理論，勢流體為無旋、無粘、無熱傳遞、不可壓或輕微可壓縮的介質(zhì).當(dāng)不考慮結(jié)構(gòu)阻尼且流體速度很小時，基于勢流體的結(jié)構(gòu)動力特征值問題基本方程為

(1)

式中：F(j)=-iφ(j),其中，φ(j)為流體節(jié)點速度勢向量；K為結(jié)構(gòu)剛度矩陣；M為結(jié)構(gòu)質(zhì)量矩陣；U，F(xiàn)分別為結(jié)構(gòu)和流體節(jié)點位移列向量；MFF，CFU，KFF分別為勢流體的質(zhì)量、剛度和阻尼矩陣；(KUU)S為由于流固耦合而產(chǎn)生的剛度矩陣；ωj為耦合系統(tǒng)第j階圓頻率.

式(1)將動水效應(yīng)產(chǎn)生的附加質(zhì)量、附加阻尼，以及其對系統(tǒng)剛度的影響都計入，通過求解式(1)即可得到結(jié)構(gòu)-水耦合系統(tǒng)的各階頻率，但式(1)為一非標準特征值問題，在ADINA程序中可采用行列式搜索法或者蘭所思迭代法進行求解.

2 墩-水耦合數(shù)值分析模型

以某深水矩形橋墩為背景，墩高為60 m，墩身截面為長7.5 m、寬4.5 m、壁厚0.8 m的矩形空心截面，沿墩高方向截面無變化.利用ADINA軟件建立橋墩-水耦合系統(tǒng)完全數(shù)值分析模型，見圖1.為消除外水域邊界的影響，水域范圍取半徑為60 m的圓柱體.采用20節(jié)點3維實體單元對墩身進行離散，采用20節(jié)點基于勢的流體單元模擬墩身內(nèi)、外水.假定水體無波動，具有理想的無窮遠邊界，水體表面設(shè)為自由液面，四周設(shè)為無限遠邊界，底部為固壁邊界；模型材料參數(shù)：墩身為C40混凝土，密度2 500 kg/m3，彈性模量32.5 GPa，泊松比0.2；水體密度1 000 kg/m3，體積模量2.2 GPa.邊界條件：假設(shè)橋墩墩底固結(jié)，墩頂自由；水域表面為自由液面，水域四周為無限遠邊界，水域底面采用剛壁界面；水-結(jié)構(gòu)接觸面為流固耦合界面.

2.2 網(wǎng)格劃分與計算工況

分析表明，流固耦合計算結(jié)果是否收斂與水體豎向網(wǎng)格尺寸，以及流固耦合邊界兩側(cè)水體及結(jié)構(gòu)網(wǎng)格大小直接相關(guān).魏凱等[8]對流固耦合分析中網(wǎng)格劃分問題有較詳細的探討.本文在前人研究的基礎(chǔ)上，通過不斷調(diào)整網(wǎng)格密度，以及試算的方法，確定模型網(wǎng)格尺寸大小，以保證計算結(jié)果收斂性，限于篇幅，文中不再闡述具體過程.

對于處于深水中的空心橋墩，為平衡墩身受力往往會在墩身開孔，使墩身內(nèi)部也充滿水體.為此文中分析針對以下3種情況進行探討：工況1，只考慮空心墩墩身內(nèi)水的影響；工況2，只考慮空心墩墩身外水的影響；工況3，同時考慮空心墩內(nèi)、外水的影響.

考慮水域深度的變化，取水深變化范圍為0～55 m，對于工況3而言，假定內(nèi)、外水域頂面齊平.

3 深水空心高墩振動特性結(jié)果分析

3.1 不考慮上部結(jié)構(gòu)的影響

為分析水深對矩形高墩結(jié)構(gòu)振動特性的影響，引入量綱一的量參數(shù)：相對水深Rw=(水深Hw/墩高Hp)；周期增幅RT=(結(jié)構(gòu)有水自振周期Tw-結(jié)構(gòu)無水自振周期T0)/結(jié)構(gòu)無水自振周期T0.由于橋墩縱橋向與橫橋向振動特性以及動力響應(yīng)并不耦聯(lián)，故可對其分別進行探討.研究表明，高階振型對高墩結(jié)構(gòu)的地震響應(yīng)貢獻較大，因此，文中除了探討墩-水耦合作用對結(jié)構(gòu)低階振動的影響，還將分析其對結(jié)構(gòu)高階振動的影響.

圖2為工況3、無水條件下橋墩結(jié)構(gòu)前3階振型圖，圖3為工況3、水深45 m時深水橋墩結(jié)構(gòu)縱橋向前3階振型圖.由圖2～3可知，有水條件下橋墩結(jié)構(gòu)的振型圖與無水情況下基本一致，對于橫橋向或其他水深條件下也是類似，限于篇幅，不再列舉.

圖2 無水條件下橋墩各階振型矢量圖(工況3)

圖4為不同水深條件下，工況1～3計算得到的橋墩結(jié)構(gòu)縱橋向前3階振型周期增幅變化曲線；圖5為不同水深條件下，工況1～3計算得到的橋墩結(jié)構(gòu)橫橋向前3階振型周期增幅變化曲線.

圖3 深水橋墩振型矢量圖(工況3)

圖4 深水橋墩前3階振動周期變化(縱橋向)

圖5 深水橋墩前3階振動周期變化(橫橋向)

由圖4～5可知，有水條件下，橋墩各階振型周期均有所增大，且隨著水深的增大而增大；不同工況下橋墩同一階振型周期變化曲線基本一致，但影響程度，內(nèi)、外均有水時振動周期增幅最大，僅外有水次之，而僅內(nèi)有水時最??；對于同一工況下不同階振型，振動周期變化規(guī)律并不完全相同，第1階振型周期增幅速度隨著水深增大而加快，即曲線中間沒有拐點，第2階振型周期增幅速度隨著水深增大為先增大而后減小，即曲線中間有一個拐點，第3階振型周期增幅則更為復(fù)雜，曲線中間有2個拐點，對照相應(yīng)的振型圖可以看出，振型周期增幅曲線圖與振型圖十分類似；對于同一工況下，墩-水耦合作用對結(jié)構(gòu)同一方向第2階、第3階振型周期的影響大于對第1階振動周期的影響.

對比橋墩3個方向各階振型周期的變化幅值可以看出，墩-水耦合作用對縱橋向的影響大于對橫橋向的影響，但工況1相差不大，如水深55 m時，縱橋向前3階振動周期增幅依次為12.8%，14.2%，15.4%，而橫橋向?qū)?yīng)工況周期最大增幅降低至12.1%，13.5%，14.2%；工況2和工況3相差較大，如工況3在水深55 m時，縱橋向前3階振動周期增幅依次為45.6%，50.0%，51.2%，而橫橋向?qū)?yīng)工況周期最大增幅降低至27.5%，30.5%，31.6%；此因橋墩為矩形截面，垂直于順橋向截面邊長(7.5 m) 大于橫橋向(4.5 m) ，地震波沿順橋向作用時，順橋向動水壓力大于橫橋向，故動水作用對橋墩縱橋向周期影響率大于橫向橋.

以上還顯示，當(dāng)相對水深小于0.5 m時，無論縱橋向還是橫橋向，橋墩第1階振型周期的增幅很小，而第2階、第3階振型周期的增幅則較大，如當(dāng)相對水深為0.5時，工況3計算的增幅分別達到了30%和20%左右.可見，當(dāng)水深不大時，雖然墩-水耦合作用對結(jié)構(gòu)第1階振型周期的影響很小，但對結(jié)構(gòu)第2,3階振型周期影響較大.由于高階振型對高墩結(jié)構(gòu)地震響應(yīng)貢獻較大，因此，在低水位時也不可忽略墩-水耦合作用的影響.

4.2 上部結(jié)構(gòu)附加質(zhì)量的影響

實際橋梁結(jié)構(gòu)中墩頂承載著梁體及橋面系質(zhì)量，其將影響橋墩的振動特性.為考慮上部結(jié)構(gòu)的影響，文中將上部結(jié)構(gòu)簡化為墩頂附加集中質(zhì)量，假定為1 500 t.與上節(jié)類似，探討墩頂附加質(zhì)量的影響，經(jīng)計算得出的不同水深條件下縱橋向和橫橋向，工況1～3橋墩結(jié)構(gòu)前3階振型周期變化曲線見圖6～7.

圖6 深水橋墩前3階振動周期變化

圖7 深水橋墩前三階振動周期變化(橫橋向)

由圖6～7可知，考慮墩頂附加質(zhì)量時，有水條件下橋墩縱橋向與橫橋向前3階振型周期變化規(guī)律與不考慮上部結(jié)構(gòu)質(zhì)量時大體一致，但對于同一工況且相同水深條件下，上部結(jié)構(gòu)質(zhì)量的存在明顯降低了結(jié)構(gòu)第1階振型周期的增幅，如工況1～3，不考慮上部結(jié)構(gòu)質(zhì)量時第1階振動周期最大增幅依次為12.8%，35.9%，45.6%，而考慮上部結(jié)構(gòu)質(zhì)量后，對應(yīng)工況周期最大增幅降低至3.8%，11.2%，14.7%；而第2、3階振動周期的增幅反而有所增大，如工況3，不考慮上部結(jié)構(gòu)質(zhì)量時其第2、3階振動周期最大增幅依次為42.8%，43.4%，而考慮上部結(jié)構(gòu)質(zhì)量后，對應(yīng)周期最大增幅增長至46.9%，47.3%.可見，考慮上部結(jié)構(gòu)質(zhì)量的影響，雖然會減弱流固耦合作用對低階振型的影響，但會增其大對結(jié)構(gòu)高階振型的影響，因此，分析深水高墩結(jié)構(gòu)地震響應(yīng)時不可忽略墩-水耦合作用對結(jié)構(gòu)高階振型的影響.

5 結(jié) 論

1) 墩-水耦合作用使得橋墩結(jié)構(gòu)各階振型周期有所增大，但各階振型周期的變化規(guī)律并不完全相同，第1階振型周期增幅速度隨著水深增大單調(diào)加快，而第2階振型周期增幅隨著水深增大為先增大后減小，第3階振型周期增幅隨著水深增大變化更為復(fù)雜，其變化圖與對應(yīng)結(jié)構(gòu)的振型圖類似.

2) 相對僅考慮外域水體作用而言，同時考慮內(nèi)、外水域時動水壓力對橋墩振動特性的影響增大，所以當(dāng)橋墩存在內(nèi)域水體時其內(nèi)域水體的動水壓力作用不容忽視.

3) 當(dāng)水深不大時，橋墩結(jié)構(gòu)第1階振型周期變化很小，但其第2，3階振型周期變化較大，特別是考慮上部結(jié)構(gòu)質(zhì)量時更為明顯.由于高階振動對高墩結(jié)構(gòu)動力響應(yīng)的貢獻較大，因此，在低水位時也不可忽略墩-水耦合作用的影響.

4) 考慮上部結(jié)構(gòu)質(zhì)量的影響，雖然會減弱流固耦合作用對橋墩結(jié)構(gòu)低階振型的影響，但會增大對結(jié)構(gòu)高階振型的影響，因此，分析深水高墩結(jié)構(gòu)動力響應(yīng)時不可忽略墩-水耦合作用對結(jié)構(gòu)高階振型的影響.

[1]GOYAL A., CHOPRA A K. Earthquake analysis of intake-outlet towers including tower-water-foundation-soil intcraction[J]. Earthquake Engineering & Structural Dynamics,1989,8(3):325-344.

[2]杜修力,趙娟,趙密.大直徑深水圓柱結(jié)構(gòu)動水壓力的時域算法[J].土木工程學(xué)報,2012(2):316-320.

[3]楊吉新,張可,黨慧慧.基于ANSYS的流固耦合動力分析方法[J].船海工程,2008,37(6):85-89.

[4]江輝,楚芹,王寶喜,等.近、遠場地震下深水橋墩動力響應(yīng)特性對比研究[J]．振動與沖擊,2011,33(2):58-66.

[5]楊萬理,李喬.墩-水耦合計算模式及深水橋墩動力響應(yīng)研究[J].地震工程與工程振動,2012,32(3):130-137.

[6]魏凱，伍勇吉，徐燦,等.橋梁群樁基礎(chǔ)-水耦合系統(tǒng)動力特性數(shù)值模擬[J]．工程力學(xué)，2011，28(1):195-200.

[7]WEI K, YUAN W C. Experimental and numerical assessment of the three-dimensional modal dynamic response of bridge pile foundations submerged in water[J]. Journal of Bridge Engineering，2013，18(3):1032-1041.

Study on Vibration Characteristics of Deep-water Rectangular Hollow High Bridge Pier Considering Fluid-structure Interaction

DENG Yulin GUO Qingkang HE Xiongjun

(SchoolofTransportation,WuhanUniversityofTechnology,Wuhan430063,China)

In this paper, a typical deep-water rectangular hollow high pier is taken as an example and the full numerical method is used to build the numerical models of pier-water interaction system based on the potential flow theory. The vibration characteristics of the deep-water pier are investigated for three scenarios including only with water outside the hollow pier, only with water inside the hollow pier, and with water both inside and outside the hollow pier. The research results indicate that the vibration periods of the pier increase due to the effect of pier-water interaction, with the maximum increase corresponding to the case of water both inside and outside the hollow pier. As the water lever is smaller than half of the pier height, the periods of low-order mode vibration of the pier change slightly, while the periods of high-order mode vibration change significantly. Although the mass of superstructure could mitigate the effect of fluid-structure interaction on the low-order modes of vibration of the pier, it would cause large effects on the high-order modes of vibration of the pier since the high-order modes of vibration contribute more significantly to the bottom shear of the pier, Therefore, great attention should be paid to the effect of pier-water coupling on the high-order modes of vibration of the pier when the water depth is not large enough.

bridge; deepwater bridge pier; rectangular hollow high pier; fluid-structure interaction; vibration characteristics

2016-09-28

*國家自然科學(xué)基金項目資助(51378406)

U442.5+5

10.3963/j.issn.2095-3844.2016.06.006

鄧育林(1977—)：男，博士，副教授，主要研究領(lǐng)域為橋梁抗震