廖袖鋒，魏奇科，胡李俐，肖正直，丁鑫

（1重慶市合川區(qū)建設(shè)工程質(zhì)量監(jiān)督站，重慶401520；2中冶建工集團(tuán)有限公司，重慶400084；3西華師范大學(xué)基建處，四川南充637009；4重慶大學(xué)土木工程學(xué)院，重慶400045）

基于V2F湍流模型的坡屋面體型系數(shù)研究

廖袖鋒1，魏奇科2，胡李俐3，4，肖正直4，丁鑫4

（1重慶市合川區(qū)建設(shè)工程質(zhì)量監(jiān)督站，重慶401520；2中冶建工集團(tuán)有限公司，重慶400084；3西華師范大學(xué)基建處，四川南充637009；4重慶大學(xué)土木工程學(xué)院，重慶400045）

在土木工程中，兩方程k-ε系列湍流模型在結(jié)構(gòu)風(fēng)荷載模擬中有著廣泛的應(yīng)用。但是這些湍流模型都需結(jié)合壁面函數(shù)使用，從而使得在壁面附近以及鈍體轉(zhuǎn)角處計(jì)算精度不高。該文在°penF°AM平臺(tái)上，基于V2F湍流模型對(duì)低矮坡屋面風(fēng)場(chǎng)進(jìn)行數(shù)值模擬，研究坡屋面體型系數(shù)的變化規(guī)律。經(jīng)過(guò)非線(xiàn)性修正的V2F湍流模型在壁面附近及鈍體轉(zhuǎn)角處都有很高的計(jì)算精度，并可以不再使用壁面函數(shù)。通過(guò)對(duì)坡屋面體型系數(shù)的模擬，分析了地貌、湍流強(qiáng)度對(duì)體型系數(shù)的影響，并與規(guī)范進(jìn)行了對(duì)比。該文最后還研究了坡屋面建筑之間的干擾效應(yīng)及其對(duì)體型系數(shù)的影響，為工程設(shè)計(jì)提供了參考。

低矮坡屋面；體型系數(shù)；風(fēng)荷載；°penF°AM；V2F湍流模型

引言

風(fēng)荷載是結(jié)構(gòu)設(shè)計(jì)的重要荷載，處理不當(dāng)會(huì)造成危及生命財(cái)產(chǎn)的安全事故[1]。在我國(guó)西部以高山大嶺和山地為主，到處分布著大量的低矮建筑。由于復(fù)雜山地環(huán)境對(duì)低矮建筑的風(fēng)荷載的作用復(fù)雜，屋面受風(fēng)荷載影響而破壞較為常見(jiàn)，建筑地形地貌對(duì)其結(jié)構(gòu)的影響越來(lái)越受到重視[2]。

數(shù)值模擬是近年來(lái)發(fā)展較快的一種結(jié)構(gòu)風(fēng)工程研究辦法，而要獲得精度較高的建筑表面風(fēng)壓就對(duì)數(shù)值求解方法提出了較高要求[3]。湍流在壁面附近區(qū)域受壁面的約束，湍流發(fā)展不充分，屬于低雷諾數(shù)流體。所以如標(biāo)準(zhǔn)k-ε模型、RNG k-ε模型等常見(jiàn)高雷諾數(shù)模型就無(wú)法對(duì)近壁區(qū)進(jìn)行直接求解，必須結(jié)合壁面函數(shù)使用，但是壁面函數(shù)畢竟是一種經(jīng)驗(yàn)簡(jiǎn)化后的函數(shù)，還是有一定的局限性。V2F模型是經(jīng)過(guò)非線(xiàn)性修正的湍流模型，它在近壁區(qū)為各向異性的，并且需要在近壁區(qū)加密網(wǎng)格，不使用壁面函數(shù)，在近壁區(qū)域可以視其為低雷諾數(shù)模型。本文在OpenFOAM平臺(tái)中使用能夠直接處理低雷諾數(shù)近壁區(qū)的V2F模型進(jìn)行數(shù)值模擬以求得相對(duì)精度較高的結(jié)果。

1 V2F湍流模型

V2F湍流模型[4]它是建立在二階矩基礎(chǔ)上的湍流模型。它基于雷諾應(yīng)力模型和直接數(shù)值模擬，在標(biāo)準(zhǔn)k-ε湍流模型基礎(chǔ)上，增加了速度尺度變量υ2和橢圓松弛方程f。其中速度尺度變量是垂直于壁面方向的湍流強(qiáng)度的表示量。橢圓松弛方程用來(lái)反映近壁面處流動(dòng)和非局域效應(yīng)的。該湍流模型為四方程模型，其控制方程為：

式中，L為長(zhǎng)度尺度，T為時(shí)間尺度，其表達(dá)式為：

另外，該湍流模型也對(duì)湍流粘性系數(shù)進(jìn)行了修正，其表達(dá)式為：

式（1）-（4）中其余各常數(shù)分別為：C1ε=1.3+0.25/d為計(jì)算點(diǎn)到壁面的距離；Cμ=0.19，C1=1.4，C2= 0.3，CT=6，σε=1.3，CL=0.3，Cη=70。有別于原始模型，修正后的模型參數(shù)均與到壁面的距離以及局部雷諾數(shù)無(wú)關(guān)。雖然該模型多了兩個(gè)控制方程，但是在計(jì)算過(guò)程中不需要計(jì)算壁面距離或者壁面函數(shù)，因此計(jì)算量增加有限，計(jì)算效率可以接受。不少研究[5]已經(jīng)證明了該模型具有出色的湍流預(yù)測(cè)能力和較強(qiáng)的普適性。它對(duì)壁面附近的分離、撞擊給予準(zhǔn)確的預(yù)測(cè)，近年來(lái)得到了廣泛的使用。

2 坡屋面風(fēng)荷載模擬

2.1 模型建立及計(jì)算邊界

坡屋面建筑分布廣泛，并且其結(jié)構(gòu)形式多樣，綜合考慮實(shí)際尺寸以及建筑長(zhǎng)寬比等因素[7]，本文確定研究對(duì)象的平面尺寸為12m×6m，屋檐高度為5m，屋脊高度根據(jù)坡度變化，屋面坡角分別為15°、20°、25°、30°、40°、45°、50°、60°。在數(shù)值模擬中，為提高計(jì)算效率，網(wǎng)格在距離建筑較遠(yuǎn)處稀疏，較近處密集。為了滿(mǎn)足V2F湍流模型近壁面計(jì)算處理的要求，壁面處使用三菱柱邊界層網(wǎng)格進(jìn)行加密，最小網(wǎng)格約為0.005，如圖1所示。

圖1網(wǎng)格示意圖

根據(jù)《建筑結(jié)構(gòu)荷載規(guī)范》[8]的規(guī)定，入口風(fēng)速剖面和湍流度剖面分別為：

其中α為地貌指數(shù)，v10和I10分別為10m高度處的風(fēng)速和湍流強(qiáng)度。

在數(shù)值模擬中，入口為速度入口，壓力為zero Gradient；出口條件為壓力出口，壓力定為0，其余量設(shè)置為zero Gradient；左右兩側(cè)及頂部為對(duì)稱(chēng)邊界，設(shè)置為symmetry邊界；底部為固定邊界。

2.2 屋面壓力分布

在B類(lèi)地貌下，分別對(duì)8種不同屋面坡度的模型進(jìn)行計(jì)算，得到各屋面的表面平均壓力系數(shù)云圖，如圖2。

圖2坡屋面平均風(fēng)壓系數(shù)分布圖

由圖2可以看出，在迎風(fēng)面坡度較小時(shí)兩端分布有低負(fù)壓區(qū)，隨著坡度增加低負(fù)壓區(qū)向中部轉(zhuǎn)移；在風(fēng)壓變化為正壓后，兩區(qū)域?yàn)楦邏簠^(qū)并且開(kāi)始合并分布于迎風(fēng)面左側(cè)區(qū)域。在背風(fēng)面風(fēng)壓始終為負(fù)，并且隨著坡度增加在屋面中心區(qū)域形成一低氣壓區(qū)，容易對(duì)屋面覆面結(jié)構(gòu)造成破壞。

2.3 地面粗糙度和屋面坡度影響分析

《建筑結(jié)構(gòu)荷載規(guī)范》中對(duì)結(jié)構(gòu)風(fēng)荷載的計(jì)算是通過(guò)體型系數(shù)來(lái)體現(xiàn)的。體型系數(shù)是將得到的建筑物表面風(fēng)壓值除以該點(diǎn)前方未受擾動(dòng)的氣流動(dòng)壓經(jīng)無(wú)量化得到的風(fēng)壓系數(shù)的加權(quán)平均。低矮坡屋面一般用在A、B、C類(lèi)地貌下，現(xiàn)對(duì)8種坡度下屋面在三類(lèi)地貌下進(jìn)行數(shù)值模擬，得到體型系數(shù)分布規(guī)律如圖3所示。

圖3體型系數(shù)對(duì)比圖

從迎風(fēng)面體型系數(shù)圖3(a)對(duì)比可以看出：建筑的體型系數(shù)比較規(guī)律并且是與地面粗糙度有關(guān)。迎風(fēng)面壓力最?。w型系數(shù)為0）的屋面坡度為25°～30°；分別對(duì)3種不同地貌下的迎風(fēng)面體型系數(shù)進(jìn)行擬合，其結(jié)果如下式，式中x為坡度：

從背風(fēng)面體型系數(shù)圖3(b)對(duì)比可以看出：背風(fēng)面體型系數(shù)同樣隨坡角增大發(fā)生變化，但是變化幅度不大，與《建筑結(jié)構(gòu)荷載規(guī)范》的建議值-0.5相比，差值在0.1以?xún)?nèi)。

2.4 湍流強(qiáng)度影響分析

在《建筑結(jié)構(gòu)荷載規(guī)范》中，湍流強(qiáng)度對(duì)體型系數(shù)的影響沒(méi)有得到體現(xiàn)。本節(jié)以B類(lèi)地貌粗糙度下，I10分別為5%、10%、14%及20%時(shí)，對(duì)不同坡度坡屋面的體型系數(shù)進(jìn)行數(shù)值模擬。其中，I10為14%是規(guī)范中規(guī)定的B類(lèi)地貌湍流強(qiáng)度值。

圖4不同迎風(fēng)面坡度、湍流強(qiáng)度下的體型系數(shù)

從圖4可知體型系數(shù)隨著湍流強(qiáng)度的增大而增大，坡度相同時(shí)，體型系數(shù)與湍流強(qiáng)度呈正相關(guān)關(guān)系。某些工況下，體型系數(shù)的取值略有超過(guò)現(xiàn)行荷載規(guī)范，如坡度為15°、湍流強(qiáng)度為5%時(shí)迎風(fēng)面體型系數(shù)為-0.616，超過(guò)了規(guī)范的-0.6。在湍流強(qiáng)度為20%時(shí)，也有部分工況超過(guò)了規(guī)范取值。在復(fù)雜山地環(huán)境下，湍流發(fā)展較為充分，湍流強(qiáng)度比較大，局部可達(dá)20%～30%。所以在高湍流度地區(qū)，需適當(dāng)提高體型系數(shù)的取值。對(duì)于特別復(fù)雜的山區(qū)及重要建筑，則應(yīng)結(jié)合風(fēng)洞試驗(yàn)或數(shù)值模擬的結(jié)果確定體型系數(shù)。

圖5不同背風(fēng)面坡度、湍流強(qiáng)度下的體型系數(shù)

從圖5可以看出背風(fēng)面體型系數(shù)隨著湍流強(qiáng)度的增加，體型系數(shù)不斷減小。雖然表現(xiàn)出不斷減小的趨勢(shì)，但是數(shù)值上波動(dòng)不大，且在規(guī)范規(guī)定的值-0.5以?xún)?nèi)，故而可按照規(guī)范取值。

3 坡屋面的干擾效應(yīng)

前面的數(shù)值模擬都是針對(duì)單個(gè)建筑進(jìn)行的，但在實(shí)際工程中，一般都是以建筑群形式出現(xiàn)的。建筑之間有相互干擾效應(yīng)，這對(duì)于坡屋面的體型系數(shù)的影響也是值得研究的。如1965年英國(guó)Ferrybridge熱電廠的8座冷卻塔群中的3個(gè)在暴風(fēng)中被吹毀，調(diào)查研究表明干擾相應(yīng)使上游兩塔之間產(chǎn)生加速流，加大了后方的風(fēng)荷載。在《建筑結(jié)構(gòu)荷載規(guī)范》中對(duì)高層建筑群的干擾作用作了相應(yīng)的說(shuō)明，通過(guò)干擾系數(shù)來(lái)考慮干擾效應(yīng)的影響，但是對(duì)低矮建筑規(guī)范就沒(méi)有給出相關(guān)要求?，F(xiàn)在通過(guò)兩棟雙坡屋面的數(shù)值模擬來(lái)討論其影響。房屋模型的幾何尺寸仍然與前面一致，屋面坡度取為15°。定義干擾系數(shù)Ic為：

式中，μsi為受擾建筑的體型系數(shù)，μs為單體建筑的體型系數(shù)。

干擾建筑與受擾建筑之間的相互關(guān)系如圖6示例。

圖6建筑布置位置示意圖

3.1 串列排布的影響

對(duì)建筑物前后串列排布進(jìn)行研究，Sy=0，Sx分別取0.5～8b，共14個(gè)位置。其中b為建筑的寬度。由于前方建筑的遮擋效應(yīng)，后方建筑的風(fēng)荷載作用復(fù)雜，從圖7所示的流線(xiàn)圖可以看出。

圖7建筑間的速度云圖流線(xiàn)圖

體型系數(shù)及其干擾系數(shù)見(jiàn)圖8。由于受擾建筑位于干擾建筑的后方，所以受擾建筑的表面風(fēng)壓一般要小于無(wú)干擾作用，干擾系數(shù)均小于1。

從圖8可以看出迎風(fēng)面體型系數(shù)從0.5～1.5b干擾較小，干擾系數(shù)約0.1；間距為8b時(shí)，前方建筑對(duì)后方建筑的干擾作用依然存在，干擾系數(shù)為0.7。背風(fēng)面的體型系數(shù)與迎風(fēng)面有所不同，呈雙折線(xiàn)變化。間距為0.5b逐漸增大到2b時(shí)，體型系數(shù)逐漸變?。婚g距從2b增加到8b，體型系數(shù)逐漸增大。雖然從變化曲線(xiàn)上看規(guī)律明顯，但其數(shù)值差值并不大，最大差值僅為0.07；隨著間距的增大，干擾系數(shù)先減小再增大，干擾系數(shù)值在0.67～0.82之間變化。

圖8前后串列時(shí)體型系數(shù)和干擾系數(shù)

3.2 并列排布的影響

并列排布時(shí)，Sx=0，Sy分別為0.5～4b，共8個(gè)位置。與串列不同，并列的兩棟建筑物相鄰較近時(shí)，流場(chǎng)存在相互干擾疊加的情況，干擾系數(shù)一般大于1。迎風(fēng)面和背風(fēng)面的體型系數(shù)及干擾系數(shù)見(jiàn)圖9。

圖9并列排布時(shí)體型系數(shù)和干擾系數(shù)

從圖9可知：在干擾建筑的影響下，受擾建筑迎風(fēng)面和背風(fēng)面的干擾系數(shù)均大于1，并且隨著距離的增大而逐漸減小，當(dāng)距離大于4b時(shí)干擾系數(shù)接近1。當(dāng)Sy=0.5b時(shí)，迎風(fēng)面體型系數(shù)達(dá)到最大值，并且隨著距離增大而迅速減小；背風(fēng)面體型系數(shù)在Sy=1b時(shí)達(dá)到最大值。

3.3 干擾系數(shù)分布圖

在10b×6b區(qū)域選取多個(gè)點(diǎn)進(jìn)行計(jì)算，分別得到不同位置處迎風(fēng)面和背風(fēng)面的干擾系數(shù)。其中b的意義同上，為建筑的寬度。從圖10可以看出，干擾建筑物給后方的受擾建筑的影響主要體現(xiàn)在體型系數(shù)有所減小，并且越近越明顯；在兩側(cè)時(shí)體型系數(shù)會(huì)增大，但是增大的比例最大為22%。

圖10干擾系數(shù)分布圖

4 結(jié)論

本文基于V2F湍流模型，在OpenFOAM平臺(tái)中對(duì)低矮坡屋面建筑的風(fēng)荷載進(jìn)行了數(shù)值模擬，得出以下結(jié)論。

（1）經(jīng)過(guò)非線(xiàn)性修正的V2F湍流模型可以不再使用壁面函數(shù)，并在壁面附近及鈍體轉(zhuǎn)角處都有很好的模擬效果。

（2）分別對(duì)A、B、C三類(lèi)地貌粗糙度下不同坡度的坡屋面進(jìn)行的分析表明，不同地貌粗糙度下體型系數(shù)有明顯不同。論文給出了不同坡度下，迎風(fēng)面體型系數(shù)計(jì)算公式。不同工況時(shí)，背風(fēng)面的體型系數(shù)變化不大，可按照規(guī)范建議取為-0.5。

（3）湍流強(qiáng)度的增加會(huì)增加坡屋面的吸力，特別是在迎風(fēng)面，因此在復(fù)雜山地環(huán)境等高湍流度區(qū)域應(yīng)適當(dāng)增加迎風(fēng)面的體型系數(shù)。

（4）坡屋面建筑的互相干擾作用時(shí)，干擾建筑會(huì)對(duì)兩側(cè)的風(fēng)壓起增強(qiáng)作用（體型系數(shù)增大），后方由于建筑物遮擋，體型系數(shù)會(huì)變小。

[1]黃本才，汪從軍.結(jié)構(gòu)抗風(fēng)分析原理及應(yīng)用[M].上海：同濟(jì)大學(xué)出版社，2001.

[2]陳平.地形對(duì)山地丘陵風(fēng)場(chǎng)的數(shù)值研究[D].浙江：浙江大學(xué)，2007.

[3]Murakami S.Current status and future trends in computational wind engineering[J].Journal of Wind Engineering and Industrial Aerodynamics，1997:67，68:3-34.

[4]Gianluca Iaccarino，Paul Durbin.Unsteady 3D ANS simulations using the v2f model[J].Center for Turbulent Research.(Stanford)-Annual Research Briefs，2000:263-269.

[5]Lien F.S，Kalitzin G.Computations of transonic flow with the v2f turbulence model[J].International Journal of Heat and Fluid Flow 2001(22):53-61.

[6]Castro.I.P，Robins A G.The flow around a surfacemounted cube in uniform and turbulent streams[J].Journal of Fluid Mechanics，1977，79:307-335.

[7]付國(guó)宏.低層房屋風(fēng)荷載特性及抗臺(tái)風(fēng)設(shè)計(jì)研究[D].浙江：浙江大學(xué)，2002.

[8]住房和城鄉(xiāng)建設(shè)部.GB50009-2012建筑結(jié)構(gòu)荷載規(guī)范[S].北京：中國(guó)建筑工業(yè)出版社，2012.

責(zé)任編輯：孫蘇，李紅

Study on Slope Roof Shape Factors Based on V2F Turbulence Model

The 2-equation k-ε serial turbulence model is extensively applied in the structural wind load simulation in civil engineering.However,these turbulence models all need to be used in combination with wall functions,causing low computational accuracy near wall surfaces and at the corners of blunt bodies.This paper studies the change rules of slope-roof shape factor based on the numerical simulation of low slope roof's wind field by V2F turbulence model and on the platform of OpenFOAM.V2F turbulence model of non-linear revision presents high computational accuracy near wall surfaces and at the corners of blunt bodies,and the wall function is no longer needed.Through the simulation of slope roof shape coefficient,it analyzes the influence of terrain,turbulence intensity to the shape coefficient,and the analysis result is compared with the standard values.This paper also provides the interfering effect among slope-roof buildings,and the influence to shape coefficient,which is for the reference of engineering design.

low slope roof;shape factor;wind load;OpenFOAM;V2F turbulence model

TU311.3

：A

：1671-9107（2016）11-0049-05

10.3969／j.issn.1671-9107.2016.11.049

2016-08-29

廖袖鋒(1983-)，男，重慶人，研究生，高級(jí)工程師，主要從事建筑節(jié)能與建設(shè)工程管理相關(guān)工作。