高 霖 王明振,2

(1.重慶文理學院建筑工程學院，重慶 402160； 2.中國地震局工程力學研究所 中國地震局地震工程與工程振動重點實驗室，黑龍江 哈爾濱 150080)

地面式鋼筋混凝土圓形水池液固耦合模態(tài)分析★

高 霖1王明振1,2

(1.重慶文理學院建筑工程學院，重慶 402160； 2.中國地震局工程力學研究所 中國地震局地震工程與工程振動重點實驗室，黑龍江 哈爾濱 150080)

對地面式鋼筋混凝土圓形水池進行了考慮液固耦合效應的模態(tài)分析，并使用ADINA有限元軟件，分析了儲液狀態(tài)、儲液高度、水池半徑以及液體表面重力波對模態(tài)的影響，對比了儲液晃動頻率的理論和有限元計算結果，推導得到了水池儲液晃動基本周期計算公式。

水池，液固耦合，ADINA，模態(tài)分析，晃動周期

0 引言

給水系統(tǒng)是典型的城市生命線系統(tǒng)，是城市健康運營的物質基礎。生命線系統(tǒng)的功能以整個系統(tǒng)的運營狀態(tài)為體現，其中每一個環(huán)節(jié)對整個系統(tǒng)的功能均有不可忽略的作用。給水系統(tǒng)功能有取水、儲水、凈水和輸配水。給水系統(tǒng)中，完成儲水和凈水功能角色的則是水池，水池是給水系統(tǒng)中的重要構筑物，對整個系統(tǒng)的功能有重要影響。歷次地震震害表明：水池結構在強烈地震作用下易發(fā)生不同程度的破壞。結構的模態(tài)特性對其在地震作用下的響應有重要影響。通過模態(tài)分析可以得到結構的頻率、振型及模態(tài)參與質量等動力特性，由動力特性根據反應譜可確定結構的地震影響系數。1957年，Housner教授對剛性儲液結構的模態(tài)和液動壓力進行了分析[1]，后續(xù)較多專家學者也對儲液結構模態(tài)進行了分析[2-5]，但大多數結果是針對剛性儲液結構給出的。我國GB 50032—2003室外給水排水和燃氣熱力工程抗震設計規(guī)范中尚未給出儲液水池的模態(tài)計算方法，不便于儲液水池的地震響應計算。同時，由于水池這類結構在正常使用時包括固體結構和液體兩相，內部液體的存在會影響水池結構的模態(tài)特性，因此在研究儲液結構模態(tài)特征時必須考慮兩相間液固耦合效應的影響。

綜上所述，為了得到合理的儲液水池模態(tài)分布特征，本文使用ADINA有限元軟件和理論公式法對地面式鋼筋混凝土圓形水池進行考慮液固耦合效應的模態(tài)分析，研究模態(tài)分布的規(guī)律，為計算儲液水池的地震響應奠定基礎。

1 液固耦合水池有限元模型

水池結構尺寸的確定主要參照相關規(guī)范、圖集和實際工程案例，所分析水池的結構形式均為地面式、無頂蓋的圓形水池。將水池的結構形式、名稱、容量、底厚、壁厚、儲液半徑、最大儲液高度及配筋等基本參數列于表1中。為分析不同儲液量對水池振動模態(tài)的影響，同時考慮水池的日常工作狀態(tài)，選取從無水至70%儲液量之間共8個儲液工況(如表2所示)對水池A，B，C分別進行模態(tài)分析。為便于分析結果的表述，使用編號“C-70%”表示容量2 000 m3水池、儲水高度為2.45 m的工況，其他工況以此類推。ADINA中液體材料的密度為1 000 kg/m3，體積模量2.3×109Pa，阻尼比0.16%[6]。

表1 圓形水池基本參數

表2 儲水量與儲水高度對應表

使用ADINA Parasolid幾何建模方式建立水池模型。液體采用3D-Fuild勢流體單元，可模擬與結構單元的耦合問題。池體結構采用3D-Solid單元，為混凝土材料，使用Truss單元中的Rebar選項設置鋼筋[7]。模型網格劃分時兼顧計算精度和速度。根據研究對象的特點，選用行列式搜索法(Determinant-Search Method)進行模態(tài)分析，同時得到模態(tài)參與因子，用以表征某階模態(tài)對結構振動的貢獻程度。

2 模態(tài)計算結果分析

進行液固耦合系統(tǒng)模態(tài)分析時，首先對模型施加重力荷載進行靜力分析，在模態(tài)分析時選擇靜力計算結果進行重啟動，進行模態(tài)和參與因子計算。3種不同尺寸水池、8種儲液狀態(tài)、是否考慮表面重力波影響2種情況，共計算48個工況。對比分析各工況的計算結果，分析不同因素對模態(tài)的影響程度和趨勢。

2.1 儲液狀態(tài)對模態(tài)的影響

計算無液和儲液狀態(tài)時水池的模態(tài)，分析儲液狀態(tài)對模態(tài)的影響。無液時，水池的振型主要為池壁的振動，且振動頻率較高。當水池內儲有液體時，無論儲液量的多少，液固耦合系統(tǒng)包括剛體模態(tài)、儲液低階晃動模態(tài)和結構高階振動模態(tài)三種。模態(tài)計算時第一階模態(tài)為剛體模態(tài)，其頻率非常小，無限接近于零，主要表現為流體所有節(jié)點上的壓力值相同，固體結構在均布內壓下的變形，剛體模態(tài)的參與質量百分比幾乎為零。以A-50%工況為例，圖1為其剛體模態(tài)。從第二階開始直至結構高階振動模態(tài)之間全部為儲液晃動模態(tài)，儲液晃動模態(tài)的頻率相對較小，主要表現為池內液體的晃動，水池結構雖在晃動壓力下變形，但相對液體晃動而言，池壁的振動響應是微小的。圖2為A-50%工況的儲液晃動一階模態(tài)。結構振動模態(tài)頻率較高，不易被激發(fā)，主要表現為池壁的振動，池內液體為脈沖模態(tài)，結構一階振動模態(tài)如圖3所示。當儲液量不同時，出現結構振動模態(tài)的階數略有差異。

對比分析各工況模態(tài)計算結果可知：液體的存在改變了結構的振動模態(tài)，由單一的結構振動轉變?yōu)榘▌傮w模態(tài)、儲液低階晃動模態(tài)和結構高階振動模態(tài)三種。不同儲液量時，相同階數三種模態(tài)的振型幾乎完全相同；儲液時的結構高階振型與無液時的結構振型幾乎完全相同，且出現的順序基本相同。儲液水池的儲液晃動模態(tài)很容易被激發(fā)，且其振型參與質量較大，對結構的振動響應影響較大。由于結構的旋轉對稱性，振動模態(tài)大多成對出現，僅振型方向不同。

2.2 儲液高度對模態(tài)的影響

將A水池不同儲液高度下的剛體模態(tài)頻率、前三階儲液晃動頻率和結構高階振動頻率列于表3中。

表3 A水池不同儲液高度時若干階模態(tài)頻率

Hz

從表3中可以看出：儲液時的第一階模態(tài)均為剛體模態(tài)，其頻率無限小接近于零；隨著儲液高度的增大，同一階次的儲液晃動頻率逐漸增大，且儲液高度越大，數值增大越緩慢；對比結構振動頻率的變化趨勢可以分析液體對結構振動的影響，同一階次結構振動頻率隨著儲液高度的增大而不斷減小，說明儲液高度越大，液體對結構振動的影響程度越大，使得液固耦合效應越顯著。

2.3 液體表面重力波對模態(tài)的影響

表面重力波存在于液體和氣體兩種流體的分界面上，以表面波的形式存在，并以重力為恢復力[8]。前述計算水池動力特性時均考慮了液體表面重力波，即液面為自由液面，其約束為Free Surface。當約束住液面的Z向位移，即約束了液面的自由晃動，也即不考慮液體表面重力波，通過ADINA計算可直接得到結構振動模態(tài)。將水池A儲水高度為10%，30%，50%和70%的考慮和不考慮液體表面重力波時的結構振動頻率列于表4中。

表4 考慮和不考慮液體表面重力波時水池A的結構振動頻率 Hz

從表4中可以看出：儲液高度相同時，對于同一階結構振動頻率，考慮液體表面重力波比不考慮時的略小，但相差甚微，說明表面重力波對結構振動模態(tài)有一定的影響，但影響程度較小。同時發(fā)現考慮與不考慮表面重力波時的模態(tài)振型基本相同，這與文獻[9]的研究結果相同。因此，求解結構模態(tài)時可以不考慮液體表面重力波的影響。由于約束了液面的Z向位移，不考慮液體表面重力波時得不到儲液晃動模態(tài)。

2.4 水池半徑對模態(tài)的影響

將儲液高度為70%的A，B和C水池前三階儲液晃動頻率和結構振動頻率列于表5中。計算時均考慮液體表面重力波。

對于儲液晃動模態(tài)，從表5中可以看出，儲液晃動頻率數值隨水池半徑的增大而減??；對比圖2，圖4和圖5，水池半徑不同時同一階儲液晃動振型基本相同。對于結構振動模態(tài)，頻率數值也隨水池半徑的增大而減小，水池半徑越大，液固耦合效應越顯著；水池半徑不同時，結構振動振型被激發(fā)的順序不同。水池半徑對儲液晃動和結構振動頻率均有影響，水池半徑對儲液晃動振型影響較小，對結構振動振型影響較大。

表5 儲液高度70%的水池A，B和C的模態(tài)頻率計算結果 Hz

3 儲液晃動頻率有限元法和理論法計算結果對比

GB 50032—2003室外給水排水和燃氣熱力工程抗震設計規(guī)范中尚未給出儲液水池晃動基本周期計算公式。文獻[3]中給出了基于如下假定的儲液結構內液體振動頻率計算公式，見式(1)，認為儲液結構為絕對剛體，液體為無粘、無旋、不可壓縮的理想液體，且做小振幅運動。

(1)

由式(1)計算所分析水池的儲液晃動頻率，列于表6中，并與ADINA計算結果進行對比。

表6 儲液晃動頻率理論與有限元計算結果對比 Hz

從表6中可以看出：不同儲液高度、不同水池半徑下，儲液晃動前三階頻率的理論和ADINA計算結果相差很小，也說明了所建立的有限元模型合理可靠。由于有限元計算時儲液結構具有一

定的彈性模量，而并非剛性，液面為自由液面且考慮表面重力波的影響，使得兩種方法計算的頻率值有一定的差異?；谏鲜龇治?，推導得到圓形水池儲液晃動基本周期計算公式，見式(2)。儲液晃動基本周期是計算儲液結構地震響應的重要參數。

(2)

其中，T1為儲液晃動基本周期，s；R為儲液半徑，m;HW為儲液高度，m。

4 結語

本文對地面式鋼筋混凝土圓形水池進行了考慮液固耦合效應的模態(tài)分析，主要得到以下5條結論：1)儲液結構由于內部液體的存在，改變了結構的動力特性。儲液結構模態(tài)包括剛體模態(tài)、以液體晃動為主的儲液低階晃動模態(tài)和以結構振動為主的結構高階振動模態(tài)三種。2)水池半徑相同時，儲液高度越大，同一階次儲液晃動頻率越大，同一階次結構振動頻率越小。儲液高度越大，液固耦合效應越顯著。3)液體表面重力波對儲液結構模態(tài)影響較小。不考慮液體表面重力波時，ADINA軟件可直接求得結構振動模態(tài)。4)儲液高度相同時，水池半徑越大，同一階次儲液晃動頻率和結構振動頻率越小，儲液晃動振型基本相同，但結構振動振型被激發(fā)的順序不同。5)有限元法和理論法計算儲液晃動頻率相差很小，驗證了所建立有限元模型的可靠性。推導得到了水池儲液晃動基本周期計算公式，可用于水池結構地震響應計算。

[1] Housner GW.Dynamic pressures on accelerated fluid containers[J].Bulletin of the Seismological Society of America,1957,47(1):15-35.

[2] Haroun MA.Vibration studies and tests of liquid storage tanks[J].Earthquake Engineering and Structural Dynamics,1983,11(2):179-206.

[3] 居榮初,曾心傳.彈性結構與液體的耦聯振動理論[M].北京：地震出版社,1983.

[4] 項忠權,李清林.立式儲罐抗震[M].北京：地震出版社,1990.

[5] 陳志濤.卵形消化池流—固耦合系統(tǒng)的ADINA有限元分析[D].天津:天津大學,2007.

[6] 沈世杰.貯水池動水壓力計算探討[J].特種結構,1995,12(4):4-7.

[7] ADINA北京代表處.ADINA土木工程分析實例手冊[M].北京:ADINA代表處,2005.

[8] 高 霖.地面式鋼筋混凝土水池自愈、滲漏試驗及地震響應分析[D].哈爾濱:中國地震局工程力學研究所,2015.

[9] 舒亞俐.考慮液固耦合振動的貯液結構動力特性和地震響應研究[J].特種結構,2013,30(1):69-75.

The modal analyses considering liquid-structure coupling effect of ground reinforced concrete circular tanks★

Gao Lin1Wang Mingzhen1,2

(1.CollegeofArchitectureEngineering,ChongqingUniversityofArtsandSciences,Chongqing402160,China; 2.InstituteofEngineeringMechanics,ChinaEarthquakeAdministration,KeyLaboratoryofEarthquakeEngineeringandEngineeringVibration,ChinaEarthquakeAdministration,Harbin150080,China)

The modal of ground reinforced concrete tanks are analyzed considering liquid-structure coupling effects. The influence of liquid state, liquid height, tank radius and surface gravity waves on modal are analyzed by using ADINA software. And the theory and finite element calculation results of sloshing frequency are compared. And the sloshing foundation period formula of circular reinforced concrete tank is derived.

tank, liquid-structure coupling, ADINA, modal analysis, sloshing period

1009-6825(2017)02-0047-03

2016-11-08★：重慶文理學院校級科技項目(項目編號：R2015JJ06)；國家自然科學基金資助項目(項目編號：51678544)；國家科技支撐計劃課題(課題編號：2015BAK17B05)

高 霖(1988- )，女，博士，講師

P315.8

A