林擁軍,沈艷忱,李明水,羅　楠

(1.西南交通大學(xué)土木工程學(xué)院,四川 成都 610031; 2.南京水利科學(xué)研究院,江蘇 南京 210029)

隨著社會(huì)不斷進(jìn)步,大跨屋蓋結(jié)構(gòu)已廣泛應(yīng)用于試驗(yàn)大廳、展覽館、航空港、體育館、車站等公共建筑中,且多采用整體受力性能好的曲面結(jié)構(gòu),曲面形式對(duì)風(fēng)壓分布影響很大,不同的曲面形式可能會(huì)有不同的風(fēng)壓分布規(guī)律[1-3].除此之外,這些大跨度屋蓋結(jié)構(gòu)還具有質(zhì)量輕、柔性大、自振頻率低的特點(diǎn),對(duì)風(fēng)荷載十分敏感,屋蓋的繞流和空氣動(dòng)力作用較為復(fù)雜[4-6].其中,大跨度翹曲屋蓋結(jié)構(gòu)造型獨(dú)特優(yōu)美,但目前還沒有可供參考的風(fēng)荷載體型系數(shù),主要通過風(fēng)洞試驗(yàn)來確定.

目前國(guó)內(nèi)外研究者對(duì)大跨結(jié)構(gòu)風(fēng)荷載的研究主要有基于現(xiàn)場(chǎng)實(shí)測(cè)、風(fēng)洞試驗(yàn)和數(shù)值模擬3種方法,現(xiàn)場(chǎng)實(shí)測(cè)由于需要耗費(fèi)大量的人力、物力,因此很少采用[7],因此對(duì)結(jié)構(gòu)風(fēng)荷載的研究主要還是采用風(fēng)洞試驗(yàn)和數(shù)值模擬方法.

隨著計(jì)算機(jī)軟硬件水平的飛速進(jìn)步以及流體力學(xué)理論的發(fā)展,計(jì)算流體動(dòng)力學(xué)(computational fluid dynamics,CFD)方法在結(jié)構(gòu)風(fēng)工程中已有所應(yīng)用[8-10].Gloria等[11]對(duì)復(fù)雜建筑進(jìn)行了平均風(fēng)壓的數(shù)值模擬,王振華[8]和KIM[12]等還分別采用基于雷諾時(shí)均方程的標(biāo)準(zhǔn)k-ε、RNGk-ε、可實(shí)現(xiàn)的(realizable)k-ε和雷諾應(yīng)力模型(Reynolds stress model,RSM) 4種湍流模型對(duì)大跨屋蓋表面平均風(fēng)壓分布進(jìn)行了數(shù)值模擬,結(jié)果表明4種模型的模擬結(jié)果差異不大.雖然CFD方法和網(wǎng)格技術(shù)等還有許多需要改進(jìn)之處,但是與傳統(tǒng)風(fēng)洞試驗(yàn)相比,CFD方法不僅可以建立建筑原型尺度,而且還可以根據(jù)研究需要,方便地改變流場(chǎng)和結(jié)構(gòu)的相關(guān)參數(shù),對(duì)研究對(duì)象進(jìn)行全方位多層次分析,從而避免風(fēng)洞試驗(yàn)的不足,且周期短、成本低,已逐漸成為結(jié)構(gòu)風(fēng)工程領(lǐng)域的研究熱點(diǎn)和預(yù)測(cè)建筑物風(fēng)荷載較為有效的方法[13-16].

本文以某大跨翹曲屋蓋試驗(yàn)廳作為研究對(duì)象,首先介紹了風(fēng)洞試驗(yàn),并分析了屋蓋風(fēng)壓分布及門窗開啟狀態(tài)的影響,然后利用流體力學(xué)分析軟件CFX,采用RNGk-ε模型,通過選用較為合理的參數(shù)設(shè)定,對(duì)其進(jìn)行CFD數(shù)值模擬,并將數(shù)值模擬結(jié)果與風(fēng)洞試驗(yàn)結(jié)果進(jìn)行對(duì)比分析,驗(yàn)證了數(shù)值方法模擬計(jì)算大跨翹曲屋蓋結(jié)構(gòu)平均風(fēng)壓分布的適用性.

1　風(fēng)洞試驗(yàn)

1.1　工程概況

某試驗(yàn)廳為大跨翹曲屋蓋結(jié)構(gòu),造型獨(dú)特,結(jié)構(gòu)平面為方形,主跨為150.0 m,高為33.6 m,該建筑中部凹曲,向兩側(cè)具有不同高度的翹曲,結(jié)構(gòu)體系復(fù)雜,圖1為其結(jié)構(gòu)示意圖,50年重現(xiàn)期基本風(fēng)壓值為0.45 kN/m2,100年重現(xiàn)期基本風(fēng)壓值為0.50 kN/m2,《建筑結(jié)構(gòu)荷載規(guī)范》(GB50009—2012)對(duì)于這種造型獨(dú)特的大跨度屋蓋結(jié)構(gòu)的風(fēng)荷載計(jì)算缺乏準(zhǔn)確的體型系數(shù)規(guī)定,也無(wú)參考資料可借鑒,為確保結(jié)構(gòu)的抗風(fēng)安全,應(yīng)采用風(fēng)洞試驗(yàn),研究作用于建筑物上的風(fēng)荷載,為結(jié)構(gòu)設(shè)計(jì)提供依據(jù).

1.2　試驗(yàn)設(shè)備

試驗(yàn)在西南交通大學(xué)風(fēng)工程試驗(yàn)研究中心XNJD-3回流式低速風(fēng)洞進(jìn)行[17],該風(fēng)洞試驗(yàn)段長(zhǎng)為36 m,寬為22.5 m,高為4.5 m,試驗(yàn)風(fēng)速范圍為1.0～16.5 m/s.風(fēng)洞配備了模擬大氣邊界層的裝置,風(fēng)洞底壁設(shè)有可轉(zhuǎn)動(dòng)360°的轉(zhuǎn)盤,以變換試驗(yàn)的風(fēng)向角.采用美國(guó)Scanvalve電子掃描閥進(jìn)行測(cè)壓,Dantec熱線風(fēng)速儀進(jìn)行風(fēng)速測(cè)量.

1.3　試驗(yàn)?zāi)Ｐ图皽y(cè)點(diǎn)

采用剛性模型,綜合考慮結(jié)構(gòu)幾何尺寸和風(fēng)洞試驗(yàn)段尺寸,模型的幾何縮尺比為 1∶75,由金屬管材、復(fù)合材料、有機(jī)玻璃等制成.在屋蓋自身表面、雨棚的上下表面布置測(cè)壓點(diǎn),測(cè)試風(fēng)壓分布.根據(jù)屋蓋外形特征及風(fēng)向需要進(jìn)行測(cè)點(diǎn)布置,在屋蓋邊緣、拐角位置及大挑檐等部位測(cè)點(diǎn)要密些,其他部位要疏一些.823個(gè)測(cè)點(diǎn)布置在建筑表面區(qū)域,50個(gè)測(cè)點(diǎn)布置在門面及雨棚上下表面,圖2為測(cè)點(diǎn)布置的基本情況.

圖1　某試驗(yàn)廳屋蓋結(jié)構(gòu)(單位:m)Fig.1　Roof structure of an experimental hall (unit:m)

圖2　測(cè)點(diǎn)布置情況Fig.2　Arrangement of measuring points

1.4　試驗(yàn)方案

根據(jù)《建筑結(jié)構(gòu)荷載規(guī)范》(GB50009—2012),當(dāng)重現(xiàn)期為100年時(shí),該建筑物所在地的基本風(fēng)壓為0.5 kN/m2,基本風(fēng)速V0=28.6 m/s.所在地地貌為A類,邊界層粗糙度指數(shù)為0.12,風(fēng)剖面及湍流度分布由檔板、尖塔、粗糙元模擬(圖3).

圖3　風(fēng)洞試驗(yàn)?zāi)Ｐ虵ig.3　Wind tunnel test model

試驗(yàn)參考點(diǎn)屋蓋頂面風(fēng)速為8 m/s,每個(gè)測(cè)點(diǎn)采樣時(shí)間為60 s,采樣頻率為200 Hz.所有壓力測(cè)點(diǎn)的脈動(dòng)壓力時(shí)程將同步獲得.為了解門窗開啟狀態(tài)對(duì)結(jié)構(gòu)抗風(fēng)的影響,試驗(yàn)分為5種情況:WC-Ⅰ,門窗全部關(guān)閉;WC-Ⅱ,開一側(cè)門和窗;WC-Ⅲ,開兩側(cè)門關(guān)窗;WC-Ⅳ,關(guān)門開一邊窗;WC-Ⅴ,關(guān)門開兩邊窗.窗開啟時(shí)的情況如圖4所示.

門窗全部關(guān)閉時(shí),試驗(yàn)風(fēng)向按24個(gè)羅盤方向設(shè)置,隔15°設(shè)置一個(gè)試驗(yàn)風(fēng)向,如圖5所示.其余情況試驗(yàn)風(fēng)向按8個(gè)羅盤方向設(shè)置,隔45°設(shè)置一個(gè)試驗(yàn)風(fēng)向.每風(fēng)向重復(fù)測(cè)量2次,所有壓力測(cè)點(diǎn)的脈動(dòng)壓力時(shí)程同步獲得.

圖4　窗開啟時(shí)的試驗(yàn)?zāi)Ｐ虵ig.4　Test model of window opening

圖5　試驗(yàn)風(fēng)向示意Fig.5　Test wind direction

1.5　試驗(yàn)結(jié)果及分析

1.5.1屋蓋風(fēng)壓分布

圖6給出了門窗全部關(guān)閉時(shí),屋脊線風(fēng)壓系數(shù)(Cp)隨風(fēng)向角的變化曲線,各曲線代表的是屋脊線上不同測(cè)點(diǎn)位置.

由圖6可知:在各方向角下,屋面風(fēng)壓分布以吸力為主,和文獻(xiàn)[7-8]的結(jié)論基本一致,說明大跨屋蓋結(jié)構(gòu)具有較為相似的風(fēng)壓分布特性;在凹曲方向,迎風(fēng)面風(fēng)壓大于背風(fēng)面;最大平均負(fù)風(fēng)壓系數(shù)為-1.350,發(fā)生在風(fēng)向角為270°時(shí)來流側(cè)屋面檐口的A1-1點(diǎn);最大平均正風(fēng)壓系數(shù)為0.089,發(fā)生在風(fēng)向角為150°時(shí)來流側(cè)屋面頂部的A1-10點(diǎn);在翹曲方向,較大平均負(fù)風(fēng)壓系數(shù)分別為-0.949、-0.949、-1.028和-1.029,發(fā)生在風(fēng)向角分別為45°、135°、225°和315°時(shí)屋面頂部的A3-11點(diǎn);最大平均正風(fēng)壓系數(shù)為0.298,發(fā)生在風(fēng)向角為180°時(shí)來流側(cè)屋面檐口的A14-9點(diǎn).

根據(jù)體型系數(shù)的定義,風(fēng)洞試驗(yàn)所得到的大跨翹曲屋蓋各分區(qū)體型系數(shù)在-1.3～0.6之間,這一結(jié)果比《建筑結(jié)構(gòu)荷載規(guī)范》(GB50009—2012)對(duì)拱形屋面規(guī)定的-0.8～-0.5要大很多,說明翹曲屋蓋所受風(fēng)力作用比拱形屋面更復(fù)雜,主要原因可能是屋面檐口區(qū)域脫落的旋渦使得該區(qū)域產(chǎn)生了較大的吸力,并提高了該處的分離強(qiáng)度,導(dǎo)致負(fù)壓絕對(duì)值增大,為安全起見,大跨翹曲屋蓋屋面部分的風(fēng)荷載體型系數(shù)建議取不大于-1.3[8-9].

1.5.2門窗開啟對(duì)風(fēng)壓分布的影響

為了解門窗開啟狀態(tài)對(duì)屋蓋風(fēng)壓分布特性的影響,分別選取屋蓋檐口中部測(cè)點(diǎn)A1-1、檐口角部測(cè)點(diǎn)A25-3以及屋蓋頂部中央測(cè)點(diǎn)A1-7進(jìn)行分析,在不同風(fēng)向角下風(fēng)壓系數(shù)的對(duì)比結(jié)果如圖7所示.由圖7可知,門窗的開啟對(duì)試驗(yàn)廳屋蓋外風(fēng)壓系數(shù)影響較小,僅對(duì)開啟位置附近的測(cè)點(diǎn)風(fēng)壓系數(shù)有一定影響,除一側(cè)門窗開啟時(shí),270°方向角上風(fēng)壓系數(shù)變化較大之外,整個(gè)屋蓋上表面風(fēng)壓系數(shù)變化并不大.

不同風(fēng)向角下內(nèi)壓隨門窗開啟狀態(tài)的變化情況如圖8所示.

圖7　門窗開啟狀態(tài)對(duì)風(fēng)壓特性的影響Fig.7　Influence of open status of windows and doors on wind pressure characteristics

圖8　不同風(fēng)向角下內(nèi)壓隨門窗開啟狀態(tài)的變化情況Fig.8　Internal pressure varying according to open status of doors and windows and wind direction

由圖8可知:門窗的開啟對(duì)試驗(yàn)廳的內(nèi)壓有一定影響,不開門窗時(shí),內(nèi)壓接近于0;開一邊門或一邊窗,當(dāng)正吹時(shí),出現(xiàn)較小正壓,最大為0.087 MPa,屋蓋受到向上的升力,應(yīng)予以注意;兩邊同時(shí)開啟時(shí),內(nèi)壓為負(fù)值,屋蓋受到向下的吸力,偏于安全[17].

2　數(shù)值模擬

2.1　控制方程

在建筑結(jié)構(gòu)領(lǐng)域中所涉及到的低速流動(dòng)空氣是具有剪切應(yīng)力的牛頓流體,近地面風(fēng)的馬赫數(shù)一般比較小,可看作不可壓縮流體,其基本控制方程為時(shí)均形式的連續(xù)方程和動(dòng)量方程[7,13]如式(1)、(2).

(1)

(2)

式中:ρ為空氣密度;

μi、μj為時(shí)均速度;

xi、xj為時(shí)均位移;

p為時(shí)均壓強(qiáng);

μ為湍動(dòng)黏度.

2.2　計(jì)算域、網(wǎng)格及邊界條件

根據(jù)該建筑物的實(shí)際尺寸:150 m(長(zhǎng))×140 m(寬)×33.69 m(高),結(jié)合計(jì)算條件,將計(jì)算域取為1 600 m×1 600 m×800 m,阻塞比足以滿足模擬建筑物所處大氣環(huán)境的要求.由于該建筑物外形復(fù)雜,采用Solidwork軟件構(gòu)建其實(shí)體模型,使用流體力學(xué)網(wǎng)格劃分軟件Ansys Icem CFD進(jìn)行部分非結(jié)構(gòu)四面體網(wǎng)格劃分,然后輸入到流體力學(xué)分析軟件CFX中進(jìn)行求解.

為了準(zhǔn)確模擬屋面風(fēng)壓分布,在網(wǎng)格劃分時(shí),建筑物表面及其附近采用加密網(wǎng)格,往外網(wǎng)格尺寸逐漸增大,共劃分體網(wǎng)格2 884 882個(gè),計(jì)算域網(wǎng)格劃分情況如圖9所示.

圖9　計(jì)算域網(wǎng)格劃分Fig.9　Mesh grid of computational domain

入口邊界條件:采用速度入口模擬A類大氣邊界層風(fēng)剖面,數(shù)值模擬采用的平均風(fēng)速剖面與風(fēng)洞試驗(yàn)平均風(fēng)速剖面接近.平均風(fēng)速V1=V10(z/10)α,其中:離地面10 m高度處的平均風(fēng)速V10=25 m/s;粗糙度指數(shù)α=0.12;z為距離地面的高度.

來流湍流特性通過直接給定的湍流動(dòng)能k和湍流耗散率ε得出[18-19].

k=0.5(VzIz)2,

(3)

ε=0.093/4k3/2/l,

(4)

式中:Iz為湍流強(qiáng)度,參考日本規(guī)范[8],Iz=0.1×(z/300)-α-0.05;

l為湍流積分尺度,l=100(z/30)0.5.

出口邊界條件:采用完全發(fā)展出流邊界條件,即流場(chǎng)任意物理量沿出口法向梯度為0.在計(jì)算域頂部和兩側(cè)采用對(duì)稱邊界條件來等價(jià)黏性流動(dòng)中的無(wú)滑移壁面.在建筑物表面和地面采用無(wú)滑移的壁面條件.

2.3　計(jì)算參數(shù)

風(fēng)工程的數(shù)值模擬中涉及到湍流模式的假定,較為廣泛的湍流模型是二方程模型,諸如標(biāo)準(zhǔn)k-ε模型、RNGk-ε模型、Wilcoxk-ω模型、SSTk-ω模型等[10].RNGk-ε模型由于引入了主流時(shí)均應(yīng)變率,可以反映渦流的非各向同性性質(zhì)[20-21],較其他模型有更好的計(jì)算精度,特別是在鈍體繞流的模擬中,可實(shí)現(xiàn)比標(biāo)準(zhǔn)k-ε模型更高的精度[22],因此本文采用RNGk-ε模型.針對(duì)該模型在雷諾數(shù)Re較低時(shí)適用性降低的情況,在計(jì)算時(shí)與非平衡壁面函數(shù)結(jié)合使用[1].

為保證計(jì)算的穩(wěn)定性,并獲得較高的計(jì)算精度,對(duì)流項(xiàng)的離散選用了具有三階精度的二次迎風(fēng)插值格式,速度壓力耦合采用SIMPLEC算法,迭代計(jì)算的收斂標(biāo)準(zhǔn)設(shè)為無(wú)量綱均方根殘差降至10-5以下,且表面風(fēng)壓基本保持不變時(shí),即認(rèn)為流場(chǎng)進(jìn)入穩(wěn)態(tài)[9].

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

圖10為門窗全部關(guān)閉(WC-Ⅰ)時(shí),風(fēng)洞試驗(yàn)和數(shù)值模擬計(jì)算得到的屋面風(fēng)壓系數(shù)等值線情況.

從圖10中可以看出:0°風(fēng)向角下,屋蓋最大風(fēng)壓系數(shù)發(fā)生在第1屋蓋翹曲兩側(cè)邊緣處,計(jì)算風(fēng)壓系數(shù)為-1.27,與風(fēng)洞試驗(yàn)最大負(fù)風(fēng)壓系數(shù)測(cè)點(diǎn)A2-4位置較為接近,試驗(yàn)最大負(fù)風(fēng)壓系數(shù)為-1.06;90°風(fēng)向角下,屋蓋最大風(fēng)壓系數(shù)發(fā)生在來流方向第2屋蓋翹曲邊緣處,計(jì)算風(fēng)壓系數(shù)為-1.24,與風(fēng)洞試驗(yàn)最大負(fù)風(fēng)壓系數(shù)測(cè)點(diǎn)A3-17的位置較為接近,試驗(yàn)最大負(fù)風(fēng)壓系數(shù)為-1.54;0°和90°風(fēng)向角下,迎風(fēng)面均為壓力,而背風(fēng)面、兩側(cè)和屋蓋大部分為吸力,90°風(fēng)向角下風(fēng)壓分布的數(shù)值模擬結(jié)果稍大于風(fēng)洞試驗(yàn)結(jié)果,而0°風(fēng)向角下風(fēng)壓分布的數(shù)值模擬結(jié)果與風(fēng)洞試驗(yàn)結(jié)果吻合較好.

(b) 90°風(fēng)向角圖10　風(fēng)壓系數(shù)等值線(單位:MPa)Fig.10　Wind pressure coefficient contour (unit:MPa)

圖11為0°和90°風(fēng)向時(shí),10 m高度處水平風(fēng)速矢量圖.

(b) 90°風(fēng)向角圖11　不同風(fēng)向時(shí)10 m高度處水平風(fēng)速矢量圖Fig.11　Vectorgraph of horizontal wind speed at 10 m height for different wind directions

從圖11中可以看出:0°和90°風(fēng)向來流在迎風(fēng)屋面處均未發(fā)生流動(dòng)分離,也沒有強(qiáng)的旋渦;在0°風(fēng)向時(shí),風(fēng)流繞過建筑后,在側(cè)面伴隨著分離和漩渦產(chǎn)生,在背風(fēng)面形成了兩個(gè)較大且對(duì)稱的尾渦;90°風(fēng)向時(shí),側(cè)面和背風(fēng)面分離和漩渦均不如0°風(fēng)向明顯.

圖12為0°風(fēng)向時(shí),翹曲向屋脊線豎剖面風(fēng)速矢量圖,從圖12中可以看出:風(fēng)流在翹曲的邊緣有較大的分離,且翹曲表面有較強(qiáng)的漩渦產(chǎn)生,說明這種翹曲屋蓋所受風(fēng)力比普通大跨度屋蓋復(fù)雜得多[2].

圖13為0°風(fēng)向角下,凹曲和翹曲方向屋脊線的風(fēng)壓系數(shù).

圖12　0°風(fēng)向時(shí)翹曲向屋脊線豎剖面風(fēng)速矢量圖Fig.12　Wind speed vectorgraph of vertical cut plane of warpage roof ridge for 0° wind

圖13　0°風(fēng)向角屋脊線上風(fēng)壓系數(shù)Fig.13　Wind pressure coefficient of roof ridge for 0° wind

從圖13中可以看出,由于風(fēng)洞試驗(yàn)采用的是縮尺模型,而數(shù)值方法采用的是全尺模擬,數(shù)值模擬網(wǎng)格劃分質(zhì)量、計(jì)算參數(shù)的設(shè)定以及試驗(yàn)?zāi)Ｐ涂s尺等原因,數(shù)值模擬值與試驗(yàn)值之間存在一定的差別[8],但差別不大,且變化趨勢(shì)一致,表明采用RNGk-ε模型模擬計(jì)算大跨度翹曲屋蓋平均風(fēng)壓分布具有較好的計(jì)算精度.

綜上所述,風(fēng)荷載作用下,由于風(fēng)流在屋蓋翹曲的邊緣有較大的分離,且在翹曲表面有較強(qiáng)的漩渦產(chǎn)生,其受風(fēng)力比普通大跨度屋蓋復(fù)雜得多,對(duì)風(fēng)荷載作用也更為敏感.風(fēng)荷載是結(jié)構(gòu)設(shè)計(jì)的主要控制荷載,無(wú)論是主體結(jié)構(gòu)風(fēng)荷載計(jì)算,還是附屬構(gòu)件的設(shè)計(jì),務(wù)必予以足夠的重視.同時(shí),由于翹曲邊緣的流動(dòng)分離作用,大跨翹曲屋蓋結(jié)構(gòu)中,風(fēng)荷載最不利的位置在翹曲邊緣和屋面頂部區(qū)域,局部風(fēng)壓系數(shù)較大,無(wú)論是結(jié)構(gòu)整體設(shè)計(jì)還是局部構(gòu)件設(shè)計(jì),翹曲邊緣和屋面頂部區(qū)域應(yīng)予以重點(diǎn)注意.

3　結(jié)　論

通過對(duì)大跨翹曲屋蓋結(jié)構(gòu)風(fēng)荷載試驗(yàn)數(shù)據(jù)進(jìn)行分析,以及采用RNGk-ε模型進(jìn)行CFD數(shù)值模擬,并將數(shù)值模擬結(jié)果與風(fēng)洞試驗(yàn)結(jié)果進(jìn)行對(duì)比分析,主要結(jié)論如下:

(1) 門窗的開啟對(duì)大跨翹曲屋蓋外風(fēng)壓系數(shù)影響較小,對(duì)內(nèi)壓有一定影響,開一邊門窗時(shí),屋蓋會(huì)受到向上的升力,兩邊同時(shí)開啟時(shí),內(nèi)壓對(duì)屋蓋有向下的吸力作用.

(2)采用RNGk-ε模型模擬計(jì)算大跨度翹曲屋蓋平均風(fēng)壓分布具有較好的計(jì)算精度,來流方向?yàn)槁N曲向(0°風(fēng)向)時(shí),風(fēng)流繞過建筑后,在來流方向建筑兩側(cè)會(huì)伴隨著分離和漩渦產(chǎn)生,且在背風(fēng)面會(huì)形成大的對(duì)稱尾渦,而來流方向?yàn)榘记?90°風(fēng)向)時(shí),側(cè)面和背風(fēng)面分離和漩渦不明顯.

(3) 風(fēng)流在翹曲的邊緣有較大的分離,在翹曲表面有較強(qiáng)的漩渦產(chǎn)生,其受風(fēng)力比普通大跨度屋蓋復(fù)雜,對(duì)風(fēng)荷載作用也更敏感.

(4) 大跨翹曲屋蓋在各方向角下,屋面風(fēng)壓分布以吸力為主,風(fēng)荷載最不利位置在翹曲邊緣和屋面頂部區(qū)域,無(wú)論是結(jié)構(gòu)整體設(shè)計(jì)還是局部構(gòu)件設(shè)計(jì),翹曲邊緣和屋面頂部區(qū)域都應(yīng)予以重點(diǎn)注意.

參考文獻(xiàn):

[1]樓文娟,孫斌,盧旦,等.復(fù)雜型體懸挑屋蓋風(fēng)荷載風(fēng)洞試驗(yàn)與數(shù)值模擬[J].建筑結(jié)構(gòu)學(xué)報(bào),2007,28(1):107-112.

LOU Wenjuan,SUN Bin,LU Dan,et al.Wind tunnel tests and numerical simulation of wind loads on complex cantilevered roofs[J].Journal of Building Structures,2007,28(1):107-112.

[2]林擁軍,宋長(zhǎng)江,羅楠,等.大跨度單層網(wǎng)殼結(jié)構(gòu)風(fēng)洞試驗(yàn)研究[J].工業(yè)建筑,2013,43(7):130-134.

LIN Yongjun,SONG Changjiang,LUO Nan,et al.Study on wind tunnel tests of large-span single-layer reticulated shell structure[J].Industrial Construction,2013,43(7):130-134.

[3]馬文勇,劉慶寬,尉耀元.具有凹面外形的大跨屋蓋結(jié)構(gòu)風(fēng)荷載分布及風(fēng)洞試驗(yàn)研究[J].振動(dòng)與沖擊,2013,31(22):34-38.

MA Wenyong,LIU Qingkuan,WEI Yaoyuan.Wind load distribution and wind tunnel tests for a curved concave long-span roof[J].Journal of Vibration and Shock,2013,31(22):34-38.

[4]羅楠,廖海黎,李明水.大跨屋蓋結(jié)構(gòu)共振響應(yīng)的簡(jiǎn)化CQC 法[J].西南交通大學(xué)學(xué)報(bào),2012,47(6):915-921.

LUO Nan,LIAO Haili,LI Mingshui.Simplified CQC method for resonant response of long-Span roof structure[J].Journal of Southwest Jiaotong University,2012,47(6):915-921.

[5]黃鵬,蔡玢,全涌,等.基于實(shí)測(cè)的低矮房屋屋面風(fēng)壓極值計(jì)算方法[J].西南交通大學(xué)學(xué)報(bào),2014,49(2):247-253.

HUANG Peng,CAI Bin,QUAN Yong,et al.Peak pressure estimation method of wind loads on low-rise building based on field measurement[J].Journal of Southwest Jiaotong University,2014,49(2):247-253.

[6]UEMATSU Y,WATANABE K,SASAKI A,et al.Wind induced dynamic response and result ant load estimation of a circular flat roof[J].Journal of Wind Engineering & Industrial Aerodynamics,1999,83:251-261.

[7]傅繼陽(yáng),趙若紅,徐安,等.大跨屋蓋結(jié)構(gòu)風(fēng)效應(yīng)的風(fēng)洞試驗(yàn)與原型實(shí)測(cè)研究[J].湖南大學(xué)學(xué)報(bào):自然科學(xué)版,2010,37(9):12-18.

FU Jiyang,ZHAO Ruohong,XU An,et al.Wind tunnel and full scale study of wind effects on a large roof structure[J].Journal of Hunan University:Natural Sciences,2010,37(9):12-18.

[8]王振華,袁行飛,董石麟.大跨度橢球屋蓋結(jié)構(gòu)風(fēng)壓分布的風(fēng)洞試驗(yàn)和數(shù)值模擬[J].浙江大學(xué)學(xué)報(bào):工學(xué)版 ,2007,41(9):1462-1466.

WANG Zhenhua,YUAN Xingfei,DONG Shilin.Wind tunnel experiment and numerical simulation of wind pressure distribution of long span ellipsoidal shell structure[J].Journal of Zhejiang University:Engineering Science,2007,41(9):1462-1466.

[9]顧明,趙雅麗,黃強(qiáng),等.低層房屋屋面平均風(fēng)壓的風(fēng)洞試驗(yàn)和數(shù)值模擬[J].空氣動(dòng)力學(xué)學(xué)報(bào),2010,28(1):82-87.

GU Ming,ZHAO Yali,HUANG Qiang,et al.Wind tunnel tests and numerical simulation of mean wind pressures on roof of low-rise buildings[J].Acta Aerodynamica Sinica,2010,28(1):82-87.

[10]MURAKAMI S.Current status and future trends in computational wind engineering[J].Wind Eng.Indus.Aerodyn,1997,67/68:3-34

[11]GLORIA M,MORET R A,PEDRO M.Experimental and numerical study of wind pressures on irregular plan shapes[J].Journal of Wind Engineering and Industrial Aerodynamics,2005,93:741-756.

[12]KIM S E,FERIT B.Application of CFD to environment flows[J].Journal of Wind Engineering and Industrial Aerodynamics,1999,81:145-158.

[13]GOSMAN A.Developents in CFD for industrial and environmental applications in wind engineering[J].Journal of Wind Engineering and Industrial Aerodynamics,1999(81):21-39

[14]日本建筑學(xué)會(huì).建筑風(fēng)荷載流體計(jì)算指南[M].北京:中國(guó)建筑工業(yè)出版社,2010:19-33

[15]黃本才,汪叢軍.結(jié)構(gòu)抗風(fēng)分析原理及應(yīng)用[M].上海:同濟(jì)大學(xué)出版社,2008:125-159.

[16]王福軍.計(jì)算流體動(dòng)力學(xué)分析——CFD 軟件原理與應(yīng)用[M ].北京:清華大學(xué)出版社,2004:18-23.

[17]林擁軍,羅楠,丁巍,等.某試驗(yàn)廳風(fēng)洞試驗(yàn)研究報(bào)告[R].成都:西南交通大學(xué)風(fēng)工程試驗(yàn)研究中心,2013.

[18]馬駿,周岱,李華鋒等.大跨度空間結(jié)構(gòu)抗風(fēng)分析的 數(shù)值風(fēng)洞方法[J].工程力學(xué),2007,24(7):77-85.

MA Jun,ZHOU Dai,LI Huafeng,et al.Numerical wind tunnel technique for the wind resistance analysis of long span spatial structures[J].Engineering Mechanics,2007,24(7):77-85.

[19]林郁,卓新.開敞式叉筒網(wǎng)殼風(fēng)場(chǎng)數(shù)值模擬與受力分析[J].浙江大學(xué)學(xué)報(bào):工學(xué)版,2004,38(9):1170-1174.

LIN Yu,ZHUO Xin.Numerical simulation and mechanics analysis of wind loading on open latticed groined barrel vault[J].Journal of Zhejiang University:Engineering Science,2004,38(9):1170-1174.

[20]YAKHOT V,ORSZAG S A,THANGAM S.Development of turbulence models for shear flows by a double expansion technique[J].Physics of Fluids A,1991,4(7):1510-1520.

[21]SPECIAL C G,THANGAM S.Analysis of an RNG based turbulence model for separated flows[J].International Journal of Engineering Science,1992,10:1379-1388.

[22]KASPERSKI M,NIEMANN H J.The LRC (load response correlation) method:a general method of estimating unfavorable wind load distributions for linear and nonlinear structural behavior[J].Journal of Wind Engineering and Industrial Aerodynamics,1992,43(3):1753-1763.