盧華喜，李 軍，周葉威，梁平英

（華東交通大學(xué)土木建筑學(xué)院，江西 南昌 330013）

高速列車(chē)引起的深水橋墩流固耦合的振動(dòng)分析

盧華喜，李 軍，周葉威，梁平英

（華東交通大學(xué)土木建筑學(xué)院，江西 南昌 330013）

為研究水體對(duì)橋墩結(jié)構(gòu)振動(dòng)特性的影響，以ANSYS為計(jì)算平臺(tái)，建立了橋墩-水流固耦合有限元模型，計(jì)算了不同幾何尺寸和淹沒(méi)比情況下橋墩的自振頻率，分析了橋墩在高速列車(chē)作用下，不同水深對(duì)其振動(dòng)特性的影響。分析后得出結(jié)論：當(dāng)水體深度小于墩高50%時(shí)，橋墩的自振頻率降低不明顯；當(dāng)水體高度大于墩高50%時(shí)，橋墩自振頻率出現(xiàn)明顯降低。高速列車(chē)作用下，橋墩墩頂縱向位移出現(xiàn)極值的時(shí)間，隨淹沒(méi)比的增加向后推遲；橋墩墩頂縱向加速度隨淹沒(méi)比的增加而增加。因此，高速列車(chē)作用下，水體對(duì)于涉水鐵路橋梁橋墩的自振與動(dòng)力振動(dòng)特性有著明顯的影響，且這種影響不可忽視。

自振頻率；激振力；流固耦合；縱向位移；縱向加速度；

當(dāng)前涉水鐵路橋梁的建設(shè)，基本都是在向深水和大跨度發(fā)展，在已建和在建的涉水鐵路橋梁中，橋墩入水深度最深已達(dá)168 m。當(dāng)前的涉水鐵路橋梁設(shè)計(jì)中，大部分將橋墩所處的環(huán)境視為靜態(tài)環(huán)境來(lái)考慮，即將橋墩周?chē)牧黧w考慮為空氣而不是水體［1］。但實(shí)際上，涉水鐵路橋梁的橋墩特別是深水橋墩都是浸沒(méi)在不同深度的水中，在高速列車(chē)的作用下，這些處于深水中的橋墩會(huì)發(fā)生一定振動(dòng)和變形，并引起周?chē)w的晃動(dòng)，水體又以動(dòng)水壓力的形式反作用于橋墩，改變橋梁墩身的振動(dòng)和變形狀態(tài)，這種作用與反作用伴隨著高速列車(chē)的通過(guò)。橋墩作為橋梁的下部承載結(jié)構(gòu)，對(duì)整個(gè)橋梁的安全和穩(wěn)定起著至關(guān)重要的作用。因此，深入研究橋墩在考慮水體的影響下自振與動(dòng)力振動(dòng)特性，對(duì)橋墩乃至整座橋梁結(jié)構(gòu)安全和穩(wěn)定是至關(guān)重要的。本文主要分析在高速列車(chē)作用下，橋墩與外部水體流固耦合的振動(dòng)問(wèn)題，使用ANSYS有限元軟件建立橋墩-水流固耦合有限元模型，分析得出不同水深情況下橋墩的自振和動(dòng)力振動(dòng)特性。

1 有限元模型的建立

為了研究水體對(duì)橋墩振動(dòng)特性的影響程度，本文選用幾種幾何尺寸不同的橋墩模型進(jìn)行模擬分析，橋墩的材料屬性取彈性模量E=3.0×1010Pa，泊松比ν=0.2，密度ρ1=1 460 kg·m-3，且橋墩底部為固定端，梁體對(duì)墩頂面有一定的約束作用，故墩頂施加彈性約束。由于橋墩周?chē)虻膶挾认鄬?duì)橋墩可以認(rèn)為是無(wú)限大的，所以在建立水體模型時(shí)，取水域半徑為橋墩界面半徑的5倍較為合適［2］，水體的材料屬性取密度ρ2=1 000 kg·m-3，水中的聲速取V=1 460 m·s-1，吸收系數(shù)β=0.8。

利用ANSYS有限元軟件中的Solid45單元來(lái)模擬實(shí)體橋墩部分，F(xiàn)luid30流體單元來(lái)模擬水體部分。當(dāng)KEYOPT（2）＝0時(shí)，流體單元與結(jié)構(gòu)單元接觸，可建立與橋墩結(jié)構(gòu)相接觸的水體模型；當(dāng)KEYOPT（2）＝1時(shí)，流體單元與結(jié)構(gòu)單元沒(méi)有接觸，可建立沒(méi)有與橋墩結(jié)構(gòu)直接接觸的水體模型。Fluid30單元還必須和Fluid130單元聯(lián)合使用，用Fluid130單元構(gòu)造包圍Fluid30單元的無(wú)限外殼，此外殼被認(rèn)為是第二階吸收邊界。當(dāng)壓力波（水波）到達(dá)這個(gè)邊界時(shí)，將會(huì)被吸收，然后以最小的反射吸收到流體域內(nèi)［3-9］。選擇Fluid30流體單元作為水體模型，模擬與橋墩結(jié)構(gòu)流固耦合的模態(tài)頻率分析時(shí)，應(yīng)該滿足以下假定：假設(shè)水體是可壓縮的、非粘性的、無(wú)平均流量且不考慮水體的熱傳導(dǎo)性。模型建立完成后，在橋墩與水體的接觸面上施加FSI流固耦合標(biāo)簽，即可實(shí)現(xiàn)流固耦合，在水體外圍施加壓力為零的邊界條件。流固耦合有限元模型和模型截面類(lèi)型分別見(jiàn)圖1、圖2。

圖1 橋墩-水流固耦合有限元模型Fig.1 The finite element model of bridge pier-water fluid-solid interaction

2 水深對(duì)橋墩自振特性的影響

圖2 模型截面類(lèi)型Fig.2 The section type of model

為了分析只考慮不同水深對(duì)橋墩自振特性的影響，故選用不同半徑和墩高不同的橋墩模型進(jìn)行模擬計(jì)算，這樣可以排除橋墩尺寸對(duì)結(jié)果的影響，以及可以得到較為普遍的結(jié)論。所以，?、贅蚨瞻霃絉1=1.0 m，橋墩高度分別取H=10，20，30，40，50 m；② 橋墩半徑R1=2.0 m，橋墩高度分別取H=60，70，80，90，100 m。根據(jù)各種材料的屬性，以及橋墩的幾何尺寸、流體域的大小建立ANSYS有限元模型。為突出研究水深對(duì)橋墩自振頻率的影響，選橋墩淹沒(méi)比h/H（水深/墩高）分別為0，0.1，0.2，0.3，0.4，0.5，0.6，0.7，0.8，0.9，1.0。對(duì)模型進(jìn)行模態(tài)化分析，可以得到橋墩在不同淹沒(méi)比下的自振頻率，見(jiàn)表1和表2。

表1 R1=1.0 m時(shí)橋墩自振頻率Tab.1 The natural frequency of bridge piers when R1=1.0 m Hz

表2 R1=2.0 m時(shí)橋墩自振頻率Tab.2 The natural frequency of bridge piers when R1=2.0 m

由表1和表2可以看出，不同橋墩尺寸下，橋墩自振頻率隨著淹沒(méi)比的增加呈減小趨勢(shì)，與橋墩半徑和高度無(wú)關(guān)。當(dāng)橋墩淹沒(méi)比h/H＜0.5時(shí)，橋墩自振頻率無(wú)明顯減小，最大減小幅度為0.5%，可以忽略不計(jì)；當(dāng)橋墩淹沒(méi)比h/H＞0.5時(shí)，橋墩自振頻率出現(xiàn)明顯減小，其減小幅度隨著淹沒(méi)比的增加而增加；當(dāng)橋墩淹沒(méi)比h/H=1.0即橋墩被完全淹沒(méi)時(shí)，其自振頻率最小，相對(duì)無(wú)水情況下的自振頻率，最大減小幅度達(dá)到16.5%，故不可忽略。橋墩自振頻率變化曲線，見(jiàn)圖3和圖4。

圖3 R1=1.0 m時(shí)橋墩自振頻率變化曲線Fig.3 The curve of the natural frequency when R1=1.0 m

圖4 R1=2.0 m時(shí)橋墩自振頻率變化曲線Fig.4 The curve of the natural frequency when R1=2.0 m

3 橋墩的動(dòng)力振動(dòng)特性分析

3.1 縱向激振力及模型加載［10］

本文主要分析車(chē)速為120 km·h-1時(shí)，車(chē)橋結(jié)構(gòu)體系對(duì)橋墩的縱向激振力引起橋墩-水流固耦合的振動(dòng)問(wèn)題，縱向激振力時(shí)程圖，見(jiàn)圖5。這種激振力分布作用在墩頂上，根據(jù)圣維南原理，可將分布力等效為集中荷載作用在墩頂?shù)闹行?，橋?水流固耦合模型加載圖，見(jiàn)圖6。

圖5 激振力時(shí)程Fig.5 Time history curve of exciting force

圖6 模型加載圖Fig.6 The picture of model loading

3.2 水深對(duì)墩頂縱向位移的影響

為了研究水體對(duì)橋墩墩頂縱向位移的影響，選用橋墩的幾何尺寸分別為①橋墩半徑R1=2.0 m，墩高H=90 m，100 m，②橋墩半徑R1=2.5 m，墩高H=90 m，100 m。由第3結(jié)可以看出，當(dāng)橋墩淹沒(méi)比h/H=1.0時(shí)，水體對(duì)橋墩自振頻率的影響最為突出，故橋墩的淹沒(méi)比h/H分別取0，1.0即無(wú)水和滿水，這樣有利于分析和對(duì)比，橋墩墩頂縱向位移的時(shí)程曲線，見(jiàn)圖7、圖8。

圖7 位移時(shí)程Fig.7 Time history curves of displacement

圖8 位移時(shí)程Fig.8 Time history curves of displacement

由圖7和圖8可以看出，在高速列車(chē)作用下，同一墩高時(shí)，墩頂縱向位移隨著半徑的增加而減??；同一半徑時(shí)，墩頂縱向位移隨著墩高的增加而增加；同一半徑和墩高時(shí)，墩頂縱向位移最大值和最小值，隨著淹沒(méi)比的增加向后延遲。

3.3 水深對(duì)墩頂縱向加速度的影響

為了研究水深對(duì)墩頂縱向加速度的影響，橋墩模型尺寸同4.2節(jié)。在橋墩模型頂部施加縱向激振力，可以分析得到墩頂縱向加速度的大小，其時(shí)程曲線，見(jiàn)圖9、圖10。

圖9 加速度時(shí)程Fig.9 Time history curves of acceleration

圖10 加速度時(shí)程Fig.10 Time history curves of acceleration

由以上各圖可知，在高速列車(chē)作用下，同一墩高時(shí)，墩頂縱向加速度的大小隨半徑的增大而減??；同一半徑時(shí)，墩頂縱向加速度的大小隨墩高的增加呈減小趨勢(shì)，減小幅度并不明顯；同一半徑和墩高時(shí)，墩頂縱向加速度的大小隨淹沒(méi)比的增加而增加。這是由于高速列車(chē)作用產(chǎn)生的激振力引起橋墩的振動(dòng)，橋墩將這種振動(dòng)部分傳遞給周?chē)w，并引起水體的晃動(dòng)，水體又以動(dòng)水壓力的形式反作用于橋墩，改變橋梁墩身的振動(dòng)和變形狀態(tài)。

4 結(jié)語(yǔ)

本文建立了不同幾何尺寸和橋墩淹沒(méi)比情況下橋墩-水流固耦合有限元模型，分析了水深對(duì)橋墩自振特性的影響；通過(guò)在墩頂施加縱向激振力，研究了不同淹沒(méi)比情況下，橋墩的動(dòng)力振動(dòng)特性。并得出以下結(jié)論：①當(dāng)橋墩淹沒(méi)比h/H＜0.5時(shí)，水體對(duì)橋墩自振頻率的影響十分微小，可以忽略不計(jì)。②當(dāng)橋墩淹沒(méi)比h/H≥0.5時(shí)，橋墩自振頻率隨著橋墩淹沒(méi)比的增加出現(xiàn)明顯降低，最大降低幅度達(dá)到16.5%。③當(dāng)橋墩淹沒(méi)比h/H=1.0即橋墩被完全淹沒(méi)時(shí)，橋墩自振頻率達(dá)到最小值，且橋墩自振頻率隨淹沒(méi)比的變化規(guī)律與橋墩幾何尺寸無(wú)關(guān)。④墩頂縱向位移出現(xiàn)最大值和最小值的時(shí)間，隨著橋墩淹沒(méi)比的增加向后延遲，且與橋墩幾何尺寸無(wú)關(guān)。⑤墩頂縱向加速度隨淹沒(méi)比的增加呈增大的趨勢(shì)，且與橋墩幾何尺寸無(wú)關(guān)。

[1]屈愛(ài)平,高淑英.梁-橋墩-樁基礎(chǔ)的動(dòng)力特性研究[J].西南交通大學(xué)學(xué)報(bào),2001,36(6):641-644.

[2]柳春光,齊念.考慮流固耦合作用的深水橋墩地震響應(yīng)分析[j].防突減災(zāi)工程學(xué)報(bào),2009,29(4):433-437.

[3]楊萬(wàn)理,李喬.深水橋墩流固耦合下的自振特性研究[J].四川建筑科學(xué)研究,2012,38(3):164-170.

[4]錢(qián)若軍,董石麟,袁行飛,等.流固耦合理論研究進(jìn)展[J].空間結(jié)構(gòu),2008,14(1):3-15.

[5]楊吉新,雷凡,李昆,等.水下橋墩結(jié)構(gòu)的振動(dòng)分析[J].世界橋梁,2009,(3):40-42.

[6]崔濤.高速列車(chē)流固耦合振動(dòng)及運(yùn)行安全性研究[D].西南交通大學(xué),2011.

[7]賴(lài)偉,王君杰,韋曉,等.橋墩地震動(dòng)水效應(yīng)的水下振動(dòng)臺(tái)試驗(yàn)研究[J].地震工程與工程振動(dòng),2006,26(6):164-171.

[8]傅作新.結(jié)構(gòu)與水體的動(dòng)力相互作用問(wèn)題.[J].水利水運(yùn)科學(xué)研究,1982,(2):104-119.

[9]高芒芒,李永強(qiáng),許北軍,等.高速列車(chē)作用下的蕪湖長(zhǎng)江大橋車(chē)橋耦合振動(dòng)分析[J].中國(guó)鐵道科學(xué),2001,22(5):34-40

[10]羅文俊,雷小燕,練松良,等.車(chē)輛-高架橋耦合系統(tǒng)豎向振動(dòng)分析車(chē)輛軌道新模型[J].華東交通大學(xué)學(xué)報(bào),2013,30(2):1-8

Vibration Analysis of Deepwater Bridge Pier Based on Fluid-Solid Interaction Caused by High-speed Train

Lu Huaxi,Li Jun,Zhou Yewei,Liang Pingying
(School of Civil Engineering and Architecture,East China Jiaotong University,Nanchang 330013,China)

To research the effect of water on the vibration characteristics of bridge pier structure,this study built a bridge pier-water fluid-solid interaction finite element model based on the computing platform of ANSYS,calcu?lated the natural frequency of the bridge pier at different water depths,and then analyzed effects of water depths on the bridge pier vibration characteristics under the role of high-speed train.The conclusions drawn from the analysis were:when the water depth was less than 50%height of the pier,the natural frequency of bridge pier did not decrease significantly;when the water depth was higher than 50%height of the pier,the natural frequency of bridge pier reduced obviously;under high-speed train,with the increase of the water height the time of the longitu?dinal displacement extreme appeared lingeringly and the longitudinal acceleration of the pier increased.There?fore,it concludes that water structure has obvious effects on the static and dynamic vibration characteristics of bridge pier,which can not be neglected.

natural frequency;exciting force;fluid-solid interaction;longitudinal displacement;longitudinal acceler?ation

U443.22

A

1005-0523（2014）02-0026-06

2013-09-20

國(guó)家自然科學(xué)基金項(xiàng)目（50968007）；江西省青年科學(xué)家培養(yǎng)對(duì)象計(jì)劃（2010DQ01800）

盧華喜（1976—），男，副教授，博士，研究方向?yàn)榈鼗A(chǔ)—部結(jié)構(gòu)動(dòng)力機(jī)組作用研究。