李遇春，邸慶霜

（同濟(jì)大學(xué) 土木工程學(xué)院 水利工程系，上海 200092）

渡槽在地震或脈動(dòng)風(fēng)作用下，渡槽結(jié)構(gòu)會(huì)發(fā)生振動(dòng)，渡槽內(nèi)的流體也會(huì)發(fā)生晃動(dòng)，晃動(dòng)的流體對(duì)支撐結(jié)構(gòu)會(huì)起減振作用還是會(huì)增加結(jié)構(gòu)的振動(dòng)是一個(gè)令人關(guān)注的問(wèn)題?，F(xiàn)有的一些研究認(rèn)為［1-3］，渡槽中的流體具有TLD（調(diào)頻液體阻尼器）效應(yīng)，在分析渡槽結(jié)構(gòu)地震反應(yīng)時(shí)，流體不能作為固定質(zhì)量考慮，當(dāng)考慮流體晃動(dòng)效應(yīng)時(shí)，結(jié)構(gòu)的地震反應(yīng)一般會(huì)比流體作為固定質(zhì)量時(shí)要低，因而認(rèn)為流體的晃動(dòng)對(duì)結(jié)構(gòu)有利，即晃動(dòng)的流體具有TLD效應(yīng)。文獻(xiàn)［1-3］所論述的流體TLD效應(yīng)實(shí)際上指的是流體晃動(dòng)的減震效應(yīng)，現(xiàn)有的流體TLD效應(yīng)研究都是基于特定的算例得到的計(jì)算結(jié)論，所得到的某些結(jié)論也不盡相同，這些結(jié)論是否在一般更為廣泛的情況下都適用?有必要做進(jìn)一步的研究。

渡槽結(jié)構(gòu)工程師需要對(duì)槽內(nèi)流體的晃動(dòng)效應(yīng)有充分的理解，才能合理估計(jì)結(jié)構(gòu)的地震與或脈動(dòng)風(fēng)響應(yīng)，使得所設(shè)計(jì)的渡槽結(jié)構(gòu)具有足夠的安全儲(chǔ)備，如果誤認(rèn)為流體的晃動(dòng)任何時(shí)候都具有減振效應(yīng)，那將會(huì)對(duì)結(jié)構(gòu)安全構(gòu)成不利影響。研究渡槽中流體晃動(dòng)效應(yīng)對(duì)結(jié)構(gòu)的振動(dòng)影響，應(yīng)明確以下二個(gè)問(wèn)題：① 什么是TLD效應(yīng)?② 什么時(shí)候流體對(duì)結(jié)構(gòu)振動(dòng)具有抑制作用?什么時(shí)候會(huì)具有放大作用?本文將從理論上對(duì)以上兩個(gè)問(wèn)題展開(kāi)分析與討論。

1 什么是TLD效應(yīng)?

TLD為英文Tuned Liquid Damper的縮寫，中文譯文為調(diào)頻液體阻尼器，它實(shí)際上是一個(gè)盛有液體的容器。TLD本質(zhì)上是一個(gè)動(dòng)力吸振器，人們可以調(diào)節(jié)容器內(nèi)液體的晃動(dòng)頻率使之與結(jié)構(gòu)的基頻接近，當(dāng)結(jié)構(gòu)振動(dòng)時(shí)，TLD能從結(jié)構(gòu)中吸收（并消耗）振動(dòng)能量，這時(shí)TLD中的流體會(huì)晃動(dòng)起來(lái)，由于結(jié)構(gòu)的部分振動(dòng)動(dòng)能傳遞給了液體，因而結(jié)構(gòu)的振動(dòng)幅度會(huì)減小，這樣就達(dá)到結(jié)構(gòu)減振的目的。

由于TLD構(gòu)造簡(jiǎn)單、價(jià)格低廉，被廣泛應(yīng)用于高層（高聳）建筑的風(fēng)振控制［4-5］，TLD的設(shè)計(jì)必須滿足二個(gè)基本條件：① TLD中的流體晃動(dòng)頻率應(yīng)調(diào)到與結(jié)構(gòu)振動(dòng)主頻率（基頻）一致，這可以通過(guò)調(diào)節(jié)流體的斷面尺寸來(lái)實(shí)現(xiàn)；② TLD（子系統(tǒng)）的設(shè)計(jì)質(zhì)量應(yīng)足夠小，一般為建筑結(jié)構(gòu)總質(zhì)量的1% ～2%［4-5］，因?yàn)槿绱诵〉南鄬?duì)質(zhì)量對(duì)整個(gè)建筑結(jié)構(gòu)的主頻率（基頻）基本沒(méi)有影響，當(dāng)TLD被安裝在建筑結(jié)構(gòu)上并調(diào)頻以后，結(jié)構(gòu)的主頻率基本保持不變（或者只有細(xì)微的變化），這樣才能保證TLD的晃動(dòng)頻率與結(jié)構(gòu)主頻率一致。實(shí)際的TLD通常為淺水型，這時(shí)TLD中流體的大多數(shù)會(huì)晃動(dòng)起來(lái)，流體得到了充分的利用，這樣設(shè)計(jì)的TLD具有較高的效率。

對(duì)于渡槽而言，裝滿水的槽身支撐在結(jié)構(gòu)（混凝土排架、槽墩等）上，看起來(lái)很像一個(gè)“TLD”，然而這個(gè)“TLD”與建筑上所采用的TLD有很大的不同，首先由于渡槽輸送水的質(zhì)量與結(jié)構(gòu)本身的質(zhì)量相當(dāng)，有的甚至大于結(jié)構(gòu)本身的質(zhì)量，當(dāng)水體的質(zhì)量加入后，這時(shí)結(jié)構(gòu)（包含流體）體系的主頻會(huì)大幅下降，如果利用渡槽中的水體晃動(dòng)來(lái)達(dá)到結(jié)構(gòu)減振的目的，那么必須調(diào)節(jié)水體的晃動(dòng)頻率與結(jié)構(gòu)主頻一致，而這一點(diǎn)實(shí)際上很難做到，因?yàn)楫?dāng)水體注入渡槽后，結(jié)構(gòu)主頻已發(fā)生了很大改變；其次渡槽內(nèi)的水體通常為深水型，這時(shí)槽內(nèi)流體的只有較少的一部分會(huì)晃動(dòng)起來(lái)，大多數(shù)的流體會(huì)像固定質(zhì)量一樣不參與晃動(dòng)。

2 渡槽流體－結(jié)構(gòu)系統(tǒng)橫向動(dòng)力特性

首先需要對(duì)渡槽結(jié)構(gòu)-流體系統(tǒng)進(jìn)行模態(tài)分析才能了解系統(tǒng)的動(dòng)力特性，從而進(jìn)一步揭示流體對(duì)結(jié)構(gòu)的振動(dòng)效應(yīng)。

以下列舉四個(gè)工程實(shí)例，其渡槽結(jié)構(gòu)型式分別為雙聯(lián)排架式、雙懸臂式、拱桁架式與斜拉式渡槽（圖1～圖4）。這些渡槽的結(jié)構(gòu)尺寸與計(jì)算取值見(jiàn)文獻(xiàn)［6-7］，本文采用ANSYS程序?qū)@四種渡槽結(jié)構(gòu)系統(tǒng)進(jìn)行建模，建模的方法見(jiàn)文獻(xiàn)［8］，流體單元采用Fluid80，槽身結(jié)構(gòu)采用殼單元Shell63，框架結(jié)構(gòu)采用梁?jiǎn)卧狟eam4，斜拉索采用Link10單元（斜拉渡槽建模與動(dòng)力特性分析將另文說(shuō)明）。流體與結(jié)構(gòu)的交界面采用所謂的滑動(dòng)邊界條件，即在交界面上，強(qiáng)制流體與結(jié)構(gòu)的法向位移相等，而兩者的切向位移不加以約束。

圖1為雙聯(lián)排架式渡槽，為了說(shuō)明問(wèn)題簡(jiǎn)單計(jì)，僅取了一榀雙聯(lián)排架加以說(shuō)明；圖2為雙懸臂式渡槽的一個(gè)獨(dú)立節(jié)段；圖3為拱桁架式渡槽；圖4為斜拉渡槽。其中圖1～圖3的渡槽在過(guò)去均出現(xiàn)過(guò)風(fēng)毀事故［6］。采用Ansys程序?qū)@四種渡槽結(jié)構(gòu)進(jìn)行空間流-固耦合模態(tài)分析，可得到流-固耦合系統(tǒng)一階與二階橫向振動(dòng)振型（見(jiàn)圖1～圖4），它們分別是：① 同相位振型，這時(shí)體系的振動(dòng)頻率較低，結(jié)構(gòu)振動(dòng)的方向與流體的晃動(dòng)方向一致；② 異相位振型，這時(shí)體系的振動(dòng)頻率較高，結(jié)構(gòu)振動(dòng)的方向與流體的晃動(dòng)方向相反。

圖4 斜拉渡槽一階與二階橫向振動(dòng)振型Fig.4 The 1stand 2ndtransverse fluid-structureinteraction modes of a cable-stayed aqueduct

3 流－固耦合作用對(duì)橫向動(dòng)力特性的影響

結(jié)構(gòu)-流體系統(tǒng)橫向同相位與異相位振型受流體-結(jié)構(gòu)相互耦合的影響，流體-結(jié)構(gòu)相互耦合的強(qiáng)弱取決于支撐結(jié)構(gòu)的剛度與流體質(zhì)量的多寡，為了展示這種影響，現(xiàn)以文獻(xiàn)［7］中的雙懸臂渡槽（其橫向流-固耦合振型見(jiàn)圖2）為例加以說(shuō)明，排架柱截面為250 mm×400 mm，排架梁截面為250 mm×250 mm，靜水截面尺寸為2.05 m×2.0 m，結(jié)構(gòu)跨度為16 m。為了說(shuō)明問(wèn)題簡(jiǎn)單起見(jiàn)，對(duì)于圖2的雙懸臂渡槽，僅通過(guò)改變支撐結(jié)構(gòu)剛度來(lái)考察流體-結(jié)構(gòu)相互耦合的強(qiáng)弱對(duì)體系固有頻率的影響，結(jié)構(gòu)的剛度主要取決于結(jié)構(gòu)的布置與截面尺寸，為簡(jiǎn)單計(jì)，這里僅改變支撐結(jié)構(gòu)（雙懸臂排架）的彈性模量E來(lái)修改結(jié)構(gòu)剛度。設(shè)支撐結(jié)構(gòu)的參考彈性模量為E0=3.0×1010（N/m2），取彈性模量E的變化范圍為：E=（0.013 ～5.0）E0，采用ANSYS程序，通過(guò)空間流-固耦合模態(tài)分析，可分別得到系統(tǒng)的一階（同相位）振動(dòng)頻率ω1、二階（異相位）振動(dòng)頻率ω2隨結(jié)構(gòu)剛度變化的曲線（見(jiàn)圖5）。

圖5 雙懸臂渡槽一階（同相位）頻率ω1、二階（異相位）ω2隨結(jié)構(gòu)剛度變化的影響ig.5 Influences of the 1st（in-phase）and 2nd（out-phase）frequencies ω1，ω2due to variation of structural stiffness for a double-cantilever aqueduct

由圖5（a）可以看出，當(dāng)支撐結(jié)構(gòu)剛度較小時(shí)，流體-結(jié)構(gòu)的相互作用較為強(qiáng)烈，結(jié)構(gòu)剛度較小的改變會(huì)引起流-固耦合系統(tǒng)一階（同相位）振動(dòng)頻率ω1發(fā)生較大的改變；當(dāng)支撐結(jié)構(gòu)剛度較大時(shí)，流體-結(jié)構(gòu)的相互作用較弱，頻率ω1隨結(jié)構(gòu)剛度改變的影響很小，此時(shí)結(jié)構(gòu)的相對(duì)振幅很?。◣缀醪粍?dòng)），體系的振動(dòng)主要表現(xiàn)為流體的晃動(dòng)，體系的振動(dòng)頻率ω1趨向于流體的一階固有晃動(dòng)頻率ωl=3.78（rad/s）。

由圖5（b）可以看出，流-固耦合體系的二階（異相位）振動(dòng)頻率ω2隨結(jié)構(gòu)的剛度E近似地呈線性變化。

對(duì)于其它結(jié)構(gòu)型式的渡槽，更多的計(jì)算分析表明，渡槽流-固耦合作用影響對(duì)橫向動(dòng)力特性具有上述（圖5）相同的結(jié)論，一般的情況是：ω1≤ωl≤ω2。

4 流體的晃動(dòng)對(duì)結(jié)構(gòu)振動(dòng)的影響分析

4.1 槽內(nèi)的固定質(zhì)量、晃動(dòng)質(zhì)量效應(yīng)

為了更好地理解流體晃動(dòng)的作用，可將渡槽中的流體按等效機(jī)械模型［9］考慮。以圖6（a）的矩形渡槽為例，則槽內(nèi)流體可以簡(jiǎn)化為如圖6（b）的固定質(zhì)量M0與一個(gè)質(zhì)量-彈簧振子（M1，K1）（其中M1為晃動(dòng)質(zhì)量，K1為等效的彈簧），質(zhì)量-彈簧振子對(duì)應(yīng)流體的一階晃動(dòng)模態(tài)，由于一階晃動(dòng)模態(tài)在一般動(dòng)力反應(yīng)中占有絕對(duì)優(yōu)勢(shì)地位，因此作為工程應(yīng)用，一般只需取一階流體晃動(dòng)模態(tài)即可滿足計(jì)算精度要求，這里高階的流體晃動(dòng)模態(tài)忽略不計(jì)。圖6（b）的模型也通常稱之為Housner模型［10］。根據(jù)Housner的計(jì)算公式，對(duì)于不同高寬比的矩形渡槽，表1給出了固定質(zhì)量M0、晃動(dòng)質(zhì)量M1占總質(zhì)量的份額，表中M為總質(zhì)量，由表1可以看出，對(duì)于常見(jiàn)的高寬比0.6～0.7，槽內(nèi)固定質(zhì)量M0占多數(shù)，而晃動(dòng)質(zhì)量M1占少數(shù)，即對(duì)于大多數(shù)的（矩形）渡槽而言，槽內(nèi)的水體為深水型，槽內(nèi)流體的只有較少的一部分會(huì)晃動(dòng)起來(lái)，大多數(shù)的流體會(huì)像固定質(zhì)量一樣不參與晃動(dòng)。

固定質(zhì)量與槽體結(jié)構(gòu)同步運(yùn)動(dòng)，會(huì)加大結(jié)構(gòu)的動(dòng)力慣性反應(yīng)。槽內(nèi)晃動(dòng)質(zhì)量M1（子系統(tǒng)）使得結(jié)構(gòu)系統(tǒng)的橫向振動(dòng)模態(tài)擴(kuò)展為如前所述的兩個(gè)振型：一個(gè)為（一階）同相位振型，另一個(gè)為（二階）異相位振型。流體晃動(dòng)質(zhì)量M1在某些情況下對(duì)結(jié)構(gòu)能起到減振的作用，對(duì)結(jié)構(gòu)安全有利，但在某些情況下則會(huì)加劇結(jié)構(gòu)的振動(dòng)，對(duì)結(jié)構(gòu)安全不利。

圖6 矩形渡槽內(nèi)流體的等效機(jī)械模型Fig.6 Equivalent mechanical model of fluid in rectangular aqueduct

表1 矩形渡槽內(nèi)不同高寬比下的固定質(zhì)量M0、晃動(dòng)質(zhì)量M1的份額（M為總質(zhì)量）Tab.1 Percentages of fixed mass M0and sloshing mass M1among total mass with varying aspect ratios for a rectangular aqueduct

4.2 流體對(duì)結(jié)構(gòu)的振動(dòng)效應(yīng)

渡槽中的流體對(duì)支撐結(jié)構(gòu)是否具有減振效應(yīng)?還是會(huì)加劇結(jié)構(gòu)的振動(dòng)?主要取決于結(jié)構(gòu)系統(tǒng)的幅頻響應(yīng)特性，仍以文獻(xiàn)［7］中的雙懸臂渡槽為例加以說(shuō)明，結(jié)構(gòu)阻尼比系數(shù)取為5%，取地面加速度脈沖激勵(lì)為：

結(jié)構(gòu)剛度按E=3.0×1010（N/m2）（較?。?、E=1.5×1011（N/m2）（較大）兩種情形考慮，在式（1）的脈沖激勵(lì)下，首先求解結(jié)構(gòu)頂部位移的脈沖響應(yīng)，然后對(duì)脈沖響應(yīng)進(jìn)行Fourier變換可分別得到結(jié)構(gòu)頂部位移的幅頻響應(yīng)曲線見(jiàn)圖7。

圖7（a）為結(jié)構(gòu)剛度較小的情形，由圖可以看出，流體-結(jié)構(gòu)體系出現(xiàn)了兩個(gè)共振頻率，共振峰值對(duì)應(yīng)的頻率值為3.07（rad/s），4.91（rad/s），它們正好分別位于結(jié)構(gòu)體系的一階（同相位）頻率ω1=3.23（rad/s）與二階（異相位）頻率ω2=4.71（rad/s）附近。在一階（同相位）頻率ω1處的結(jié)構(gòu)響應(yīng)幅值較大，而在二階（異相位）頻率ω2處的結(jié)構(gòu)響應(yīng)幅值較小，在ω3=4.30（rad/s）處，結(jié)構(gòu)響應(yīng)幅值為極小值。圖7（b）為結(jié)構(gòu)剛度較大的情形，由圖可以看出，流體-結(jié)構(gòu)體系也出現(xiàn)了兩個(gè)共振頻率，共振峰值對(duì)應(yīng)的頻率值為3.68（rad/s），9.20（rad/s），它們正好在結(jié)構(gòu)體系的兩個(gè)固有振動(dòng)頻率（同相位頻率ω1=3.67（rad/s），異相位頻率ω2=9.31（rad/s）附近，在第一個(gè)固有振動(dòng)頻率ω1處的結(jié)構(gòu)響應(yīng)幅值較小，而第二個(gè)固有振動(dòng)頻率ω2處的結(jié)構(gòu)響應(yīng)幅值則較大，在ω3=4.30（rad/s）處，結(jié)構(gòu)響應(yīng)幅值為極小值。

圖7（a）、（b）有一個(gè)共同的特點(diǎn)，它們均在ω3=4.30（rad/s）處的結(jié)構(gòu)響應(yīng)取極小值，當(dāng)再改變結(jié)構(gòu)剛度時(shí)，結(jié)構(gòu)響應(yīng)幅值均在ω3處取極小值，這表明ω3與結(jié)構(gòu)本身無(wú)關(guān)，ω3僅與槽內(nèi)的流體截面尺寸有關(guān)。

圖7 結(jié)構(gòu)頂部位移幅頻響應(yīng)曲線Fig.7 Amplitude-frequency response curves of structural-top displacements

通過(guò)以上分析知，只要槽內(nèi)水體橫截面尺寸一定，不論結(jié)構(gòu)如何變化，水體均在頻率ω3=4.30（rad/s）處具有減振效應(yīng)。進(jìn)一步分析計(jì)算表明，其他結(jié)構(gòu)形式的渡槽均具有上述相同的現(xiàn)象，只是ω3的值有所不同而已。

對(duì)于一個(gè)確定的渡槽流體-結(jié)構(gòu)體系，流體什么時(shí)候具有減振效應(yīng)?什么時(shí)候會(huì)加劇結(jié)構(gòu)的振動(dòng)?取決于外動(dòng)力荷載的頻譜特性。對(duì)于脈動(dòng)風(fēng)作用而言，其卓越頻率遠(yuǎn)低于流體晃動(dòng)頻率，渡槽結(jié)構(gòu)剛度愈小，結(jié)構(gòu)的脈動(dòng)風(fēng)響應(yīng)愈大，這時(shí)流體-結(jié)構(gòu)相互作用愈強(qiáng)烈，結(jié)構(gòu)體系的一階橫向同相位振動(dòng)頻率處于較低的位置，同相位模態(tài)在脈動(dòng)風(fēng)響應(yīng)中起到主要作用，這時(shí)流體晃動(dòng)與結(jié)構(gòu)的振動(dòng)方向一致，流體的晃動(dòng)會(huì)加劇結(jié)構(gòu)的振動(dòng)，對(duì)結(jié)構(gòu)會(huì)起不利的作用；當(dāng)渡槽結(jié)構(gòu)剛度較大時(shí)，結(jié)構(gòu)的脈動(dòng)風(fēng)響應(yīng)很小，通?？梢院雎圆挥?jì)。對(duì)于地震作用而言，當(dāng)場(chǎng)地的卓越頻率與ω1（或ω2）接近時(shí)，結(jié)構(gòu)會(huì)發(fā)生共振反應(yīng)，流體的晃動(dòng)會(huì)加劇結(jié)構(gòu)的振動(dòng)，但當(dāng)場(chǎng)地的卓越頻率與ω3接近時(shí)，流體才具有減振作用（所謂的TLD效應(yīng)），結(jié)構(gòu)反應(yīng)較小。根據(jù)圖7的位移幅頻響應(yīng)曲線可以看出，當(dāng)動(dòng)力荷載的頻率很小或很大時(shí)，結(jié)構(gòu)體系的動(dòng)力反應(yīng)很小。

另外需要指出的是，本文所討論的流體的晃動(dòng)都，程，系統(tǒng)的固有頻率ω1，ω2與流體、結(jié)構(gòu)的振幅無(wú)關(guān)，現(xiàn)有的研究表明［11］：當(dāng)流體晃動(dòng)幅度較大時(shí)，流體的運(yùn)動(dòng)呈現(xiàn)非線性，系統(tǒng)的固有頻率隨流體振幅發(fā)生變化。對(duì)于TLD而言，當(dāng)流體晃動(dòng)幅度較大時(shí)，流體調(diào)頻效應(yīng)被打破，其減振效應(yīng)會(huì)大打折扣，因此水體晃動(dòng)的幅度越大，TLD減振效應(yīng)未必越大。

5 結(jié)論

（1）渡槽體系的橫向模態(tài)都呈現(xiàn)兩個(gè)基本形態(tài)：一階（同相位）振型，這時(shí)體系的振動(dòng)頻率ω1較低，支撐結(jié)構(gòu)振動(dòng)的方向與流體的晃動(dòng)方向一致；二階（異相位）振型，這時(shí)體系的振動(dòng)頻率ω2較高，支撐結(jié)構(gòu)振動(dòng)的方向與流體的晃動(dòng)方向相反。

（2）渡槽內(nèi)的流體質(zhì)量可以分為固定質(zhì)量與晃動(dòng)質(zhì)量?jī)刹糠?，其中流體中的固定質(zhì)量與槽體結(jié)構(gòu)同步運(yùn)動(dòng)，會(huì)加大支撐結(jié)構(gòu)的動(dòng)力慣性反應(yīng)，可能會(huì)對(duì)結(jié)構(gòu)安全不利。

（3）當(dāng)外動(dòng)力荷載的卓越頻率與ω1（或ω2）接近時(shí)，結(jié)構(gòu)體系會(huì)發(fā)生共振反應(yīng)，晃動(dòng)的流體會(huì)加劇支撐結(jié)構(gòu)的振動(dòng)。

（4）當(dāng)外動(dòng)力荷載的卓越頻率與ω3接近時(shí)，晃動(dòng)的流體對(duì)支撐結(jié)構(gòu)起到減震的作用。

為什么渡槽內(nèi)的流體在固定的頻率ω3處具有減振效應(yīng)，并且ω3不隨結(jié)構(gòu)的改變而改變，其內(nèi)在物理機(jī)制需要作進(jìn)一步的研究。目前渡槽的設(shè)計(jì)，由于考慮溫度等的影響，槽體與墩體間通常采用施加盆式支座或其他的柔性連接，關(guān)于盆式支座或其他的柔性連接對(duì)渡槽結(jié)構(gòu)流—固耦合振動(dòng)的影響需要作進(jìn)一步的研究。

［1］張俊發(fā)，劉云賀，王 穎，等.渡槽-水體系統(tǒng)的地震反應(yīng)分析［J］.西安理工大學(xué)學(xué)報(bào)，1999，15（4）：46-51.

［2］吳 軼，莫海鴻，楊 春.三維排架渡槽中水的調(diào)頻液體阻尼效應(yīng)［J］.水利學(xué)報(bào)，2005（9）：1115-1120.

［3］季日臣，夏修身，陳堯隆，等.考慮流一固藕合梁式矩形渡槽橫向地震響應(yīng)研究［J］.地震學(xué)報(bào)，2007，29（3）：328-333.

［4］Wakahara T，Ohyama T，F(xiàn)ujii K.Suppression of wind-induced vibration of a tall building using tuned liquid damper［J］.Journal of Wind Engineering and Industrial Aerodynamics，1992，43（1-3）：1895-1906.

［5］Sun L M，F(xiàn)ujino Y，Pacheco B M，et al.Modeling of turned liquid damper［J］. JournalofWind Engineeringand Industrial Aerodynamics，1992，41（44）：1883-1894.

［6］ Li Y C，Luo L R.Wind-destroyed mechanism analyses of three aqueductbridges［C］.Proceedingsofthe 12th International Conference on Wind Engineering， Cairns，Australia，2007：1535-1542.

［7］李遇春.某雙懸臂渡槽風(fēng)致破壞原因分析［J］.同濟(jì)大學(xué)學(xué)報(bào)，2008，36（11）：1485-1489.

［8］李遇春，樓夢(mèng)麟，尚 偉，等.大型渡槽抗震分析中流體的位移有限元模式［J］.水利學(xué)報(bào)，2003（2）：93-97.

［9］Ibrahim R A. Liquid sloshingdynamics（theoryand applications）［M］.Cambridge University Press，2005.

［10］ Housner G W.Dynamic pressure on accelerated containers［J］.Bulletin of the Seismological Society of America，1957，47（1）：15-35.

［11］ Ikeda T，Nakayama N.Nonlinear vibration of a structures caused by water sloshing in a rectangular tank［J］.Journal of Sound and Vibration，1997，201（1）：23-41.