陳肖笑

(西南交通大學(xué)土木工程學(xué)院，四川成都 610031)

索網(wǎng)結(jié)構(gòu)在建筑形式的表達(dá)上有其獨(dú)特的美感且兼具有自重輕、易折疊、收縮比大的優(yōu)點(diǎn)，這也使其被設(shè)計(jì)師們所青睞，并廣泛應(yīng)用于建筑結(jié)構(gòu)中[1]，但這類結(jié)構(gòu)質(zhì)量小、阻尼小、柔性大的特點(diǎn)致其對(duì)風(fēng)荷載十分敏感[2]。馬鞍形索網(wǎng)屋蓋結(jié)構(gòu)多應(yīng)用于大型戶外公共建筑(如體育館)等露天設(shè)施，是一種常見的大跨度索網(wǎng)結(jié)構(gòu)體系，且包括看臺(tái)、室內(nèi)場(chǎng)館等附屬建筑。對(duì)于此類位于大氣層邊界的大跨度索網(wǎng)屋蓋結(jié)構(gòu)，其鈍體繞流非常復(fù)雜，而現(xiàn)有風(fēng)荷載規(guī)范多適用于簡(jiǎn)單屋面形式，因此，研究該類結(jié)構(gòu)的風(fēng)荷載特性并且總結(jié)相應(yīng)的風(fēng)荷載規(guī)律，對(duì)實(shí)際工程中類似結(jié)構(gòu)的設(shè)計(jì)和分析是有意義的。

目前，針對(duì)于大跨度建筑結(jié)構(gòu)的風(fēng)荷載研究已趨于成熟，研究方法主要包括理論分析、戶外實(shí)測(cè)、風(fēng)洞試驗(yàn)和數(shù)值風(fēng)洞四種。理論分析受制于Navier-Stokes方程解的不確定性，只能適用于簡(jiǎn)單模型的處理，針對(duì)大跨度結(jié)構(gòu)的復(fù)雜繞流目前沒有進(jìn)展；戶外實(shí)測(cè)主要針對(duì)建成的建筑結(jié)構(gòu)，且現(xiàn)場(chǎng)環(huán)境因素復(fù)雜，較難控制；風(fēng)洞試驗(yàn)是目前研究結(jié)構(gòu)風(fēng)荷載特性的主要方法之一，其實(shí)驗(yàn)條件可以得到很好的人為控制，準(zhǔn)確度較高，不足之處是模型存在縮尺，需要通過相似準(zhǔn)則對(duì)原尺模型進(jìn)行考察；數(shù)值風(fēng)洞得益于計(jì)算流體力學(xué)(CFD)的快速發(fā)展，現(xiàn)已廣泛應(yīng)用于建筑風(fēng)工程，其能通過構(gòu)建足尺的計(jì)算模型模擬真實(shí)的建筑風(fēng)場(chǎng)，彌補(bǔ)了實(shí)驗(yàn)風(fēng)洞的不足，同時(shí)兼具周期短，費(fèi)用低的特點(diǎn)，現(xiàn)已成為相關(guān)行業(yè)的重要研究手段。

既有文獻(xiàn)對(duì)大跨度結(jié)構(gòu)的風(fēng)場(chǎng)研究多采用風(fēng)洞試驗(yàn)與數(shù)值模擬相結(jié)合的手段。在風(fēng)洞試驗(yàn)方面，1998年，K.Suresh.Kumar等人研究了不同場(chǎng)地條件下低矮建筑表面的譜分布特征，并對(duì)不同的幾何特征模型提出適宜的分析模型[3]。2004年，顧明等人對(duì)一座將建于強(qiáng)風(fēng)區(qū)的懸挑屋蓋(雨篷)結(jié)構(gòu)模型進(jìn)行了風(fēng)洞試驗(yàn), 研究了大跨度懸挑屋蓋上的風(fēng)荷載特性。不同工況下屋蓋表面的風(fēng)壓分布特性表明, 周邊建筑對(duì)所測(cè)建筑的風(fēng)荷載有較大的干擾影響[4]。在數(shù)值模擬方面，2006年，王叢軍等基于商業(yè)計(jì)算流體力學(xué)軟件，利用Baseline雷諾應(yīng)力模型對(duì)越南國家體育場(chǎng)大懸挑屋蓋的風(fēng)壓進(jìn)行了預(yù)測(cè)[5]；2010年，周晅毅等基于風(fēng)洞試驗(yàn)的結(jié)果得到某航站樓屋蓋的表面風(fēng)壓系數(shù)特性，再結(jié)合數(shù)值方法模擬了該航站樓的各個(gè)風(fēng)向角下的風(fēng)速流場(chǎng)特性，進(jìn)一步闡述了大跨屋面表面風(fēng)壓規(guī)律[6]。陳亞楠?jiǎng)t采用標(biāo)準(zhǔn)k-ε湍流模型對(duì)環(huán)狀大懸臂屋蓋風(fēng)荷載和風(fēng)場(chǎng)進(jìn)行模擬分析，分別分析了風(fēng)向角、屋蓋傾角、看臺(tái)、開洞情況等參數(shù)對(duì)屋蓋風(fēng)壓的影響作用[7]。

本文以成都露天公園大跨屋蓋結(jié)構(gòu)為例，基于商業(yè)計(jì)算流體力學(xué)軟件ANSYS CFX 19.0，研究建筑模型在不同風(fēng)向角系數(shù)下的表面風(fēng)壓系數(shù)，闡述了各風(fēng)向角下風(fēng)荷載特性及典型測(cè)點(diǎn)的風(fēng)壓系數(shù)變化和空間分布規(guī)律，該研究對(duì)類似結(jié)構(gòu)的設(shè)計(jì)具備參考意義。

1 計(jì)算模型

1.1 幾何模型

擬建的成都露天音樂公園大跨屋蓋結(jié)構(gòu)(圖1)，長跨約200 m，短跨約70 m，高度約40 m。舞臺(tái)上方結(jié)構(gòu)采用雙曲箱型鋼拱斜拉鋼索形成馬鞍雙曲面的膜結(jié)構(gòu)屋蓋。結(jié)構(gòu)跨度約為180 m，建筑跨度約為200 m，斜拱結(jié)構(gòu)中心線計(jì)算高度47.5 m。膜結(jié)構(gòu)為鋼索上下雙層膜形式，下層膜結(jié)構(gòu)最低點(diǎn)高度約為35 m。

圖1 幾何模型

1.2 湍流模型

目前, 基于雷諾平均Navier-Stokes方程 (RANS) 的湍流模型可分為兩大類：一類是求解湍流黏性系數(shù)或渦黏性系數(shù)的模型，稱為渦旋黏性模型；另一類是求解雷諾應(yīng)力輸運(yùn)方程的模型，稱為雷諾應(yīng)力模型。

數(shù)值計(jì)算采用雙方程SSTk-ω湍流模型，這是一類較為典型的渦旋黏性模型。k-ω模型的提出旨在克服k-ε模型的諸多弱點(diǎn)，如不能準(zhǔn)確模擬逆壓流動(dòng)、在黏性底層中的計(jì)算存在數(shù)值困難等[8]。但標(biāo)準(zhǔn)k-ω模型對(duì)自由流的ω值非常敏感，故Menter等提出了剪切應(yīng)力輸運(yùn)(Shear Stress Transport，SST)模型進(jìn)行了改進(jìn)[9]，該模型基于對(duì)k-ε模型和k-ω模型加權(quán)平均的方法，在壁面附近采用k-ω模型，在遠(yuǎn)離壁面的區(qū)域采用k-ε模型，降低了原先k-ω模型中邊界層外的自由流對(duì)ω的敏感度。

雙方程SSTk-ω的湍動(dòng)能方程(1)及比耗散率方程(2)可寫作[9]。

(1)

(2)

1.3 近壁區(qū)邊界層網(wǎng)格優(yōu)化

由于黏性作用，近壁區(qū)的流動(dòng)具有很大的速度梯度和壓力梯度，近壁區(qū)邊界層網(wǎng)格質(zhì)量直接決定了計(jì)算精度和計(jì)算效率，故準(zhǔn)確還原近壁區(qū)流動(dòng)是保證數(shù)值模擬精度的關(guān)鍵[10]。為描述近壁區(qū)的流動(dòng)，既有文獻(xiàn)多引入無量綱參數(shù)y+，表征離壁面最近的網(wǎng)格節(jié)點(diǎn)到壁面的距離，進(jìn)一步得到邊界層網(wǎng)格尺寸。

(3)

式中：ut為壁面剪切速度，ρ為流體密度，μ為動(dòng)力黏度。

壁面剪切速度ut可根據(jù)特征速度與模型特征長度確定：

(4)

根據(jù)計(jì)算模型幾何特性取得的邊界層參數(shù)見表1。

表1 邊界層參數(shù)

1.4 邊界條件

大氣邊界層是指地球表面值數(shù)百到一千米的大氣層，這個(gè)范圍內(nèi)的風(fēng)特性對(duì)建筑物風(fēng)荷載的影響較為顯著。A. G. Davenport[11]提出平均風(fēng)速沿高度變化的規(guī)律可用指數(shù)函數(shù)予以描述，即：

(5)

按照我國規(guī)范[12]，B類地區(qū)α=0.16，參考風(fēng)速取為10 m處22 m/s，得到的平均風(fēng)速剖面如圖2所示。

圖2 平均風(fēng)速剖面

2 計(jì)算結(jié)果與分析

2.1 風(fēng)壓系數(shù)分布

實(shí)際計(jì)算中，在0～180 °的區(qū)間內(nèi)，每隔15 °劃分風(fēng)向角，并提取模型表面的風(fēng)壓和速度。分析時(shí)，常采用無量綱的風(fēng)壓系數(shù)來描述風(fēng)對(duì)建筑物的作用：

(6)

式中：VH為屋蓋頂部高度處的平均風(fēng)速；pH為該高度處參考靜壓；pi為給定點(diǎn)處的風(fēng)壓；ρ為空氣密度。

此處以0 °和180 °風(fēng)向角為例(圖3)?？梢钥吹斤L(fēng)向角為0 °時(shí)，看臺(tái)蓋板完全處于負(fù)壓，后部由于表面上凸，呈現(xiàn)較大的負(fù)壓，前端下沉位置由于風(fēng)速較小，風(fēng)壓接近0，兩翼由于構(gòu)造較為復(fù)雜且存在空洞，出現(xiàn)小渦。而頂棚底板由于表面平整。整體呈現(xiàn)正壓，且沿流動(dòng)方向呈下降趨勢(shì)，中間處風(fēng)壓分布較為均勻。

風(fēng)向角為180 °時(shí)，看臺(tái)蓋板迎風(fēng)邊出現(xiàn)較大的負(fù)壓，沿流動(dòng)方向風(fēng)壓系數(shù)呈現(xiàn)增長趨勢(shì)，這同樣是由于前部(背風(fēng)側(cè))風(fēng)速較小造成的。頂棚底板則完全出現(xiàn)負(fù)壓，背風(fēng)側(cè)形成兩個(gè)對(duì)稱的大渦，出現(xiàn)明顯的風(fēng)吸現(xiàn)象。

頂棚蓋板的風(fēng)壓分布也是較為規(guī)律的，0 °時(shí)迎風(fēng)側(cè)表面風(fēng)速較高，負(fù)壓較大，180 °時(shí)則在中心出現(xiàn)正壓，周圍產(chǎn)生負(fù)壓。由于膜結(jié)構(gòu)表面不平整，故產(chǎn)生了眾多小渦，這通過云圖的不連續(xù)性也有所體現(xiàn)。

(a)看臺(tái)蓋板(0°)

(b)看臺(tái)蓋板(180°)

(c)頂棚底板(0°)

(d)頂棚底板(180°)圖3 表面風(fēng)壓系數(shù)云圖(來流方向?yàn)橄路?

2.2 典型部位風(fēng)壓系數(shù)隨風(fēng)向角的變化

圖4給出了頂棚底板和看臺(tái)蓋板典型部位的風(fēng)壓系數(shù)隨風(fēng)向角改變的變化規(guī)律，其中前部為0 °風(fēng)向角時(shí)的迎風(fēng)側(cè)，尾部為180 °風(fēng)向角時(shí)的迎風(fēng)側(cè)。

可以看到，頂棚底板在風(fēng)向角接近0 °或180 °時(shí)，風(fēng)壓系數(shù)變化趨于平穩(wěn)(圖4(a))。而在45～105 °的區(qū)間內(nèi)，從正壓陡跌至負(fù)壓，這是由于頂棚前部開敞，來流在前檐發(fā)生分離，形成了較大的向上壓力，但這一現(xiàn)象隨風(fēng)向角的增大而減弱。當(dāng)風(fēng)向角大于90 °時(shí)，底板處的氣流已形成渦旋，出現(xiàn)負(fù)壓。由于大風(fēng)向角時(shí)來流被看臺(tái)結(jié)構(gòu)遮擋，故此時(shí)風(fēng)壓變化不大。

頂棚蓋板由于膜結(jié)構(gòu)表面不平整，風(fēng)向角變化時(shí)曲線存在震蕩，規(guī)律性不夠明顯，需要通過修改幾何模型進(jìn)一步研究其變化規(guī)律。

看臺(tái)底板由于表面流速較高，始終處于負(fù)壓(圖4(b))。值得注意的是在風(fēng)向角為105 °時(shí)，前部(下沉位置中心)風(fēng)壓陡降，應(yīng)當(dāng)是該部位出現(xiàn)了局部的渦旋。尾部風(fēng)壓系數(shù)的變化規(guī)律表現(xiàn)出的規(guī)律性較強(qiáng)，總體呈下降趨勢(shì)，由于結(jié)構(gòu)有所傾斜，在大風(fēng)向角時(shí)產(chǎn)生的鈍體繞流導(dǎo)致風(fēng)壓有所升高。

(a)頂棚底板

(b)看臺(tái)蓋板圖4 結(jié)構(gòu)典型位置處風(fēng)壓系數(shù)隨風(fēng)向角的變化曲線

3 結(jié)論

本文主要利用計(jì)算流體力學(xué)方法，基于ANSYS CFX軟件與雙方程SSTk-ω湍流模型，引入y+值評(píng)估網(wǎng)格無關(guān)性，分析了不同風(fēng)向角下馬鞍形索網(wǎng)屋蓋及其附屬建筑的風(fēng)荷載特性。

(1)風(fēng)壓系數(shù)隨風(fēng)向角變化呈現(xiàn)較強(qiáng)的規(guī)律性，變化曲線未出現(xiàn)明顯的震蕩，可看出數(shù)值模型預(yù)測(cè)類似建筑的風(fēng)荷載是較為可靠的。

(2)仿真結(jié)果顯示風(fēng)向角在90 °附近時(shí)結(jié)構(gòu)呈現(xiàn)較大的負(fù)壓，且變化較為明顯，設(shè)計(jì)時(shí)應(yīng)予以注意。

(3)看臺(tái)前部處在高于90 °的風(fēng)向角下時(shí)，風(fēng)壓產(chǎn)生突變，需要通過風(fēng)洞實(shí)驗(yàn)等方法進(jìn)一步證實(shí)。