黃鵬

(廣州大學(xué)-淡江大學(xué)工程結(jié)構(gòu)災(zāi)害與控制聯(lián)合研究中心）

1 CFD 分析TLD 模型的理論原理

TLD 水箱內(nèi)液體三維晃動(dòng)問題可以由Navier-Stok es 方程和連續(xù)性方程來描述，并采用雷諾時(shí)均的湍流數(shù)值模擬（RANS）方法，其中Navier-Stokes 方程可用笛卡爾張量形式表示為：

上述控制方程可以采用基于有限體積法的數(shù)值模擬算法來加以求解，本文利用Fluent 軟件中基于多相流的基本模型，采用流體體積函數(shù)（Volume Of Fluid,VOF）方法捕獲晃動(dòng)液體的自由表面。VOF 方法被普遍采用來跟蹤和捕獲液體晃動(dòng)問題中的自由表面。在VOF 方法中，流體（相）之間沒有相互穿插，所有相的體積分?jǐn)?shù)之和為1。在單元中, 如果第q 相流體的體積分?jǐn)?shù)記為αq,則可能出現(xiàn)以下三個(gè)情況：

⑴αq=0：第q 相流體在單元中是空的。

⑵αq=1：第q 相流體在單元中是充滿的。

⑶0＜αq＜1：單元中包含了第q 相流體和一相或者其他多相流體的界面。

2 模型算例

為了驗(yàn)證本文采用的基于多相流的基本模型CFD數(shù)值模擬方法，模擬內(nèi)置立柱矩形TLD 系統(tǒng)中液體晃動(dòng)特征的可行性，本文采用Love 等的內(nèi)置立柱TLD 水箱振動(dòng)臺試驗(yàn)進(jìn)行對比分析。試驗(yàn)水箱尺寸如圖1 所示。

圖1 水箱平面圖及立面圖

水箱的長寬高為L×B×H=889mm×785mm×785mm，水深h=133.35mm，立柱尺寸為ax×ay=117.75mm×117.75mm，且立柱中心點(diǎn)分別位于長和寬的1/4、1/2、3/4 倍距離位置處。試驗(yàn)在水箱底部施加單向水平正弦激勵(lì)為：

式⑶中：X0為位移激勵(lì)幅值；ω 為外加激勵(lì)對應(yīng)的頻率，TLD 系統(tǒng)對應(yīng)的各階固有頻率為：

如式⑷中，n=1 則對應(yīng)為TLD 系統(tǒng)的一階晃動(dòng)頻率；g 為重力加速度，取9.8(m/s2)；h 為水深。

3 計(jì)算結(jié)果

從圖2 (a)- (b) 中可看到，在底部激勵(lì)較?。╔0/L=0.3%）時(shí)，側(cè)壁波高和水平控制力隨著激頻比的增大而增大，在激頻比達(dá)到0.92 左右時(shí)，側(cè)壁波高和水平控制力達(dá)到最大值，隨后隨著激頻比的增大，側(cè)壁波高和水平控制力逐漸減小，可看到本文數(shù)值模擬的結(jié)果與試驗(yàn)結(jié)果較為一致，吻合精度較高。圖2(c)(d)所示為側(cè)壁波高幅值隨激頻比變化的數(shù)值模擬結(jié)果與試驗(yàn)結(jié)果在底部激勵(lì)較大時(shí)（X0/L=0.6%）的對比。上述兩圖中可看到兩種不同外激勵(lì)幅值得出的數(shù)值模擬結(jié)果與試驗(yàn)結(jié)果基本一致，均在激頻比達(dá)到0.92 左右時(shí)，側(cè)壁波高和水平控制力達(dá)到最大值。由圖2 中的各圖可看到，本文所建立的CFD 數(shù)值模擬分析模型能較好地模擬內(nèi)置立柱TLD 水箱液體晃動(dòng)時(shí)所產(chǎn)生的水平控制力和側(cè)壁波高，可以為高層及超高層建筑頂部安裝有內(nèi)置立柱TLD 水箱的振動(dòng)控制效益分析提供參考。

4 結(jié)論

本文采用計(jì)算流體動(dòng)力學(xué)（CFD）方法，對內(nèi)置立柱的矩形TLD 水箱系統(tǒng)，在底部正弦位移激勵(lì)時(shí)其中液體晃動(dòng)的振動(dòng)特征，建立了數(shù)值計(jì)算仿真分析模型，基于已有試驗(yàn)數(shù)據(jù)，驗(yàn)證了該模型及數(shù)值模擬技術(shù)的正確性，可以得到以下結(jié)論：

⑴通過與已有試驗(yàn)結(jié)果的對比分析，結(jié)果表明本文建立的基于有限體積法的CFD 數(shù)值仿真分析模型可較好地模擬內(nèi)置立柱的矩形TLD 水箱在不同激頻比作用下其中液體的晃動(dòng)特性，能較好地預(yù)測箱內(nèi)液體運(yùn)動(dòng)時(shí)的側(cè)壁波高和水平控制力。

⑵通過對內(nèi)置立柱TLD 的數(shù)值模擬研究發(fā)現(xiàn)，TLD系統(tǒng)在外激勵(lì)與固有頻率相近時(shí)，將出現(xiàn)共振現(xiàn)象，側(cè)壁波高和水平控制力達(dá)到最大值，出于工程經(jīng)濟(jì)的考慮，應(yīng)使TLD 系統(tǒng)盡可能的出現(xiàn)共振，實(shí)現(xiàn)效益最大化。

圖2 數(shù)值模擬結(jié)果與試驗(yàn)結(jié)果的對比

⑶內(nèi)置立柱會(huì)導(dǎo)致TLD 固有頻率發(fā)生變化。本文算例中內(nèi)置立柱導(dǎo)致TLD 的固有頻率減小約8%。為了讓TLD 系統(tǒng)發(fā)揮出其最大的減振效果，在調(diào)諧TLD 系統(tǒng)頻率應(yīng)考慮內(nèi)置立柱對其固有頻率的影響。