王衛(wèi)華，黃漢杰

(中國空氣動(dòng)力研究與發(fā)展中心，四川綿陽621000)

屋面雪荷載分布風(fēng)洞試驗(yàn)研究

王衛(wèi)華*，黃漢杰

(中國空氣動(dòng)力研究與發(fā)展中心，四川綿陽621000)

為了預(yù)測屋面雪荷載分布，采用細(xì)石英砂粒子模擬風(fēng)吹雪進(jìn)行風(fēng)洞試驗(yàn)研究，通過降雪模擬和均勻重分布試驗(yàn)，獲得了幾種典型屋面的積雪系數(shù)分布，并與中國荷載規(guī)范進(jìn)行了比較。結(jié)果表明:當(dāng)試驗(yàn)風(fēng)速達(dá)到閾值風(fēng)速后，測量的風(fēng)洞流場有效氣動(dòng)粗糙長度隨著風(fēng)速增加而增大，顯示出風(fēng)吹雪邊界層流場的典型特性;試驗(yàn)獲得的階梯形屋面積雪系數(shù)均小于中國荷載規(guī)范值;單跨雙坡屋面的迎風(fēng)屋面積雪系數(shù)可能會(huì)超過荷載規(guī)范值，特別是20°屋面，而10°屋面積雪系數(shù)隨時(shí)間發(fā)展逐漸減小，并最終在規(guī)范值之內(nèi)。在模擬降雪試驗(yàn)中，雙跨雙坡屋面的第一個(gè)迎風(fēng)屋面積雪系數(shù)相對較大，但小于荷載規(guī)范值;在均勻重分布試驗(yàn)中，雙跨雙坡屋面的最大積雪系數(shù)均在第一個(gè)背風(fēng)屋面的屋脊后方，其值接近1．5，超過了荷載規(guī)范值。

建筑物;屋面;雪荷載;風(fēng)洞試驗(yàn)

0 引言

由屋面積雪導(dǎo)致結(jié)構(gòu)損毀的事故時(shí)有發(fā)生，如1978年美國哈特福德市政體育館雪災(zāi)事故［1］，2006年波蘭卡托維茨國際展覽會(huì)廳以及莫斯科巴斯曼室內(nèi)市場的雪災(zāi)事故［2-3］等。在中國，雪災(zāi)事故也十分頻繁，如2007年東北雪災(zāi)，僅遼寧省全省就有92家企業(yè)的300多處鋼結(jié)構(gòu)廠房遭到不同程度的損壞，許多鋼架廠房屋面板結(jié)構(gòu)坍塌或整體垮塌［4］;2008年中國南方特大雪災(zāi)，造成直接經(jīng)濟(jì)損失高達(dá)1111億元［5-6］。許多雪災(zāi)事故調(diào)查表明，屋面局部雪荷載過大是導(dǎo)致結(jié)構(gòu)損毀的重要原因［7-8］。

風(fēng)致積雪分布的研究方法包括實(shí)地觀測［9-10］、理論分析［11］、數(shù)值模擬［12-13］和風(fēng)洞試驗(yàn)研究［14-15］等，一般主要采用數(shù)值模擬和風(fēng)洞試驗(yàn)方法。風(fēng)吹雪是一種氣固兩相流，其內(nèi)在相互作用機(jī)理十分復(fù)雜。在目前的研究中，無論是數(shù)值模擬還是風(fēng)洞試驗(yàn)，還有許多地方需要改進(jìn)和完善。如數(shù)值模擬中的簡化模型與實(shí)際情況有較大差別，風(fēng)洞試驗(yàn)中模型相似準(zhǔn)則難以完全滿足等等。盡管如此，已有不少研究者通過風(fēng)洞試驗(yàn)獲得了有意義的結(jié)果。O′Rourke［14］等人利用水槽及核桃殼粉粒材料進(jìn)行屋面積雪分布特性研究。Thiis［15］等人對一類通風(fēng)屋面防滲雪進(jìn)行風(fēng)洞試驗(yàn)研究，試驗(yàn)在冷風(fēng)洞中進(jìn)行，采用人工生成雪粒，研究了不同屋檐結(jié)構(gòu)對防滲雪效果。Kimbar［16］等人根據(jù)雪粒運(yùn)動(dòng)方程，分析了風(fēng)洞試驗(yàn)?zāi)Ｐ拖嗨茀?shù)，并在此基礎(chǔ)上對體育場屋蓋雪荷載分布進(jìn)行風(fēng)洞試驗(yàn)研究。在國內(nèi)，王衛(wèi)華［17］等人通過風(fēng)洞試驗(yàn)研究了不同風(fēng)向下幾種典型屋面的積雪重分布特性，Zhou［18］等人利用風(fēng)洞試驗(yàn)研究了不同跨度的平屋面積雪重分布特性等。

隨著全球氣候環(huán)境持續(xù)惡化，極端天氣帶來的雪災(zāi)事故將越來越多。特別是近幾十年來，隨著中國城市建設(shè)規(guī)模的發(fā)展壯大，各類用途的大跨度屋面結(jié)構(gòu)不斷涌現(xiàn)，大量工業(yè)化廠房、庫房更是鱗次櫛比。因此，尋求在設(shè)計(jì)階段充分預(yù)測結(jié)構(gòu)的雪荷載分布已變得越來越重要。本文針對幾種典型屋面結(jié)構(gòu)，采用風(fēng)洞試驗(yàn)方法，以細(xì)石英砂作為模擬雪粒材料，通過模擬實(shí)際降雪和均勻降雪重分布等方法，獲得屋面的積雪分布，測量了流場參數(shù)，分析了屋面積雪系數(shù)分布特性，并與中國荷載規(guī)范［19］進(jìn)行了比較。

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

1．1 模型及相似參數(shù)

試驗(yàn)?zāi)Ｐ桶?種典型屋面結(jié)構(gòu)：1種階梯形屋面、2種單跨雙坡屋面及2種雙跨雙坡屋面，模型縮尺比為1∶40，如圖1所示。

圖1 試驗(yàn)?zāi)Ｐ褪疽鈭DFig．1 Scheme of test models

風(fēng)吹雪風(fēng)洞試驗(yàn)的模型相似準(zhǔn)則較為復(fù)雜。無量綱分析表明，相關(guān)模型相似參數(shù)多達(dá)十幾種，不可能全部都能滿足。目前大多風(fēng)洞試驗(yàn)中，主要考慮的重要相似參數(shù)包括以下幾種［20］：

式中：g為重力加速度，L為模型特征尺度，U為參考速度，u＊為流場摩阻速度，u＊t為粒子運(yùn)動(dòng)閾值摩阻速度，wf為粒子均勻沉降速度，ν為流體運(yùn)動(dòng)粘性系數(shù)。除此之外還有一些其他參數(shù)，如粒子休止角等。

公式（1）中，第1個(gè)為弗勞德數(shù)，第2個(gè)為速度比，這2個(gè)參數(shù)由風(fēng)吹粒子運(yùn)動(dòng)方程經(jīng)無量綱化分析得到；第3個(gè)參數(shù)也為速度比，是由粒子輸運(yùn)方程經(jīng)無量綱化分析得到；最后1個(gè)參數(shù)為流場粗糙長度雷諾數(shù)，為滿足流動(dòng)一致性條件，Kind［20］建議其值只要大于下限值30，一般風(fēng)洞試驗(yàn)中可以滿足。

雷諾數(shù)參數(shù)與弗勞德數(shù)是相互矛盾的，使用密度較大的粒子可以滿足雷諾數(shù)相似而不能滿足弗勞德數(shù)，反之亦然。許多試驗(yàn)結(jié)果表明，滿足雷諾數(shù)而忽略弗勞德數(shù)，試驗(yàn)結(jié)果與實(shí)測符合最好［20］，因此風(fēng)吹雪風(fēng)洞試驗(yàn)中一般使用密度較大的粒子。

本試驗(yàn)中雪粒模擬采用100目的石英砂粒子，粒子外形不甚規(guī)則，可近似為球形，直徑約為0．15mm（±0．05mm），休止角約31°，均勻沉降速度約2．0m／s。模型粒子屬性及主要參數(shù)如表1所示。

表1 粒子屬性及主要參數(shù)Table 1 Particle properties and typical similitude parameters

1．2 設(shè)備及試驗(yàn)方法

試驗(yàn)在中國空氣動(dòng)力研究與發(fā)展中心低速所一座小型下吹式直流風(fēng)洞中進(jìn)行，該風(fēng)洞試驗(yàn)段長3m，截面為0．8m×0．8m的切角正方形，常用試驗(yàn)風(fēng)速為0～35m／s，圖2為風(fēng)洞試驗(yàn)段輪廓。

圖2 風(fēng)洞試驗(yàn)段輪廓Fig．2 Picture of wind tunnel test section

試驗(yàn)前，在整個(gè)風(fēng)洞試驗(yàn)段地板覆蓋一層石英砂粒子，厚度約30mm，以模擬來流地面積雪邊界條件。石英砂粒子層表面鋪蓋平整，以避免局部堆積對邊界層流場的影響。

屋面積雪系數(shù)測量分為降雪模擬試驗(yàn)和均勻重分布試驗(yàn)。降雪模擬試驗(yàn)中，通過模擬降雪獲得吹雪氣象條件下的屋面初始積雪分布狀態(tài)；均勻重分布試驗(yàn)則在屋面初始均勻覆蓋約20mm厚的粒子層。然后進(jìn)行吹風(fēng)試驗(yàn)，獲得若干時(shí)間內(nèi)風(fēng)致屋面積雪重分布特性。模擬降雪試驗(yàn)時(shí)風(fēng)速為2．5m／s，重分布試驗(yàn)時(shí)風(fēng)速為4．6m／s，試驗(yàn)風(fēng)速以屋面高度為參考。

為實(shí)現(xiàn)降雪模擬，在模型前方風(fēng)洞上壁安裝一粒子降落裝置，該裝置為漏斗形狀，漏斗下口處有一閥門，可以控制粒子的下漏速率。

流場測量采用KANOMAX風(fēng)速計(jì)，該儀器為可伸縮熱敏探頭，不易受氣流中細(xì)小粒子影響，探頭尺寸較細(xì)，對流場影響小。

屋面粒子分布深度測量采用數(shù)顯深度卡尺及輔助工具，測量位置在屋面中心對稱面（x軸向），階梯形屋面只測量低屋面粒子分布。試驗(yàn)前，在模型屋面表面可粘貼一層砂紙，以保證屋面初始粗糙狀態(tài)。圖3為風(fēng)洞試驗(yàn)照片。

圖3 模型風(fēng)洞試驗(yàn)照片F(xiàn)ig．3 Model in the test

1．3 流場參數(shù)

積雪邊界層流場是完全氣動(dòng)粗糙的，水平平均風(fēng)速用對數(shù)律描述為：

式中：κ為馮卡門常數(shù)，取0．41；z0為氣動(dòng)粗糙長度。上式可寫為：

利用流場測量得到的風(fēng)速剖面數(shù)據(jù)以及最小二乘法，可計(jì)算出u＊及z0值。

在風(fēng)洞試驗(yàn)段地板覆蓋粒子層后，當(dāng)風(fēng)速達(dá)到一定值后（即粒子起動(dòng)閾值風(fēng)速），粒子開始運(yùn)動(dòng)。在風(fēng)吹雪（或沙）中，粒子發(fā)生躍移運(yùn)動(dòng)時(shí)，躍起的粒子從氣流中吸收能量，使周圍流速降低，相當(dāng)于流場有效氣動(dòng)粗糙長度增大了。流場有效氣動(dòng)粗糙長度與摩擦速度平方成正比［21］：

式中：g為重力加速度；C為比例系數(shù)，可利用流場測量數(shù)據(jù)擬合得到。

粒子起動(dòng)閾值風(fēng)速通過觀測得到。當(dāng)風(fēng)洞開車速控表讀數(shù)為6．4時(shí)，觀測有少量粒子運(yùn)動(dòng)，考慮到視覺誤差，確定速控表讀數(shù)6．3時(shí)為粒子起動(dòng)閾值風(fēng)速，此時(shí)對應(yīng)參考風(fēng)速為4．4m／s。

1．4 屋面積雪系數(shù)

中國建筑荷載規(guī)范［19］中，屋面雪荷載設(shè)計(jì)值計(jì)算為：

式中：sk為雪荷載標(biāo)準(zhǔn)值（k N／m2）；μr為屋面積雪分布系數(shù)；s0為基本雪壓（k N／m2）。中國建筑荷載規(guī)范中，基本雪壓定義為50年或100年重現(xiàn)期最大雪壓，μr為屋面水平投影面積上的雪荷載與s0的比值，實(shí)際就是地面基本雪壓換算到屋面雪荷載的系數(shù)。

本文屋面積雪系數(shù)與上述略有差別。在降雪模擬風(fēng)洞試驗(yàn)中，屋面積雪系數(shù)cs定義為屋面粒子沉積厚度h與水平地面平均降落粒子厚度h0的比值，其中h0計(jì)算為降落粒子總質(zhì)量與粒子密度及平均降落面積的比值。在均勻重分布試驗(yàn)中，cs定義與上述類似，只是此時(shí)h0取為屋面初始覆蓋粒子平均厚度。表2列出了部分降雪模擬試驗(yàn)中h0值。

表2 降雪模擬試驗(yàn)h0值（單位：mm）Table 2 The h0values in snowfall simulation tests（Unit：mm）

由cs及μr定義可看出，二者物理意義是相同的。如果把風(fēng)洞試驗(yàn)降雪量或初始屋面覆蓋深度由原型當(dāng)?shù)鼗狙簱Q算確定，則二者就完全一致了。

2 結(jié)果與分析

2．1 流場測量

流場測量取4種風(fēng)速，分別為風(fēng)洞速控表讀數(shù)6．0、6．3、6．5及6．8，對應(yīng)參考風(fēng)速分別為4．2、4．4、4．6和4．8m／s。圖4顯示了風(fēng)速剖面測量結(jié)果，圖中同時(shí)給出了擬合曲線。由式（3）計(jì)算出流場摩阻速度分別為0．201、0．216、0．225和0．238m／s，即u＊t＝0．216m／s。由圖可看出，在來流u＊＝0．201m／s時(shí)，氣動(dòng)粗糙長度z0最小，隨來流風(fēng)速增大，當(dāng)u＊＞u＊t時(shí)，粒子開始運(yùn)動(dòng)，此時(shí)有效氣動(dòng)粗糙長度隨之增大，表現(xiàn)出風(fēng)吹雪邊界層流場的典型特性。

圖5給出了有效氣動(dòng)粗糙長度擬合結(jié)果，擬合公式為：

由擬合公式可得式（4）中C＝0．0078。

圖4 測量的風(fēng)速剖面Fig．4 Wind speed profiles

圖5 有效氣動(dòng)粗糙長度擬合Fig．5 Effective aerodynamic－rough length

2．2 積雪系數(shù)

試驗(yàn)結(jié)果圖中橫坐標(biāo)以階梯形屋面高差或雙坡屋面的斜坡跨距進(jìn)行無量綱化。圖6顯示了階梯形屋面3次降雪模擬試驗(yàn)結(jié)果，圖中同時(shí)給出了Tsuchiya［9］的3次觀測結(jié)果平均值及中國荷載規(guī)范［19］設(shè)計(jì)值。由圖可見，屋面積雪系數(shù)均在中國荷載規(guī)范值以內(nèi)。3次試驗(yàn)結(jié)果的變化規(guī)律基本相似，主要在高低屋面交界附近（x＜1）差別略大。在x＝0處，試驗(yàn)1和3的最大值約1．25，試驗(yàn)2的最大值超過1．5。試驗(yàn)2的2次重分布時(shí)間結(jié)果變化不大，表明屋面粒子分布接近穩(wěn)定狀態(tài)。試驗(yàn)結(jié)果與Tsuchiya的觀測結(jié)果有一定差別，在x＝0．3附近，實(shí)測積雪系數(shù)突然減小至0．2，在x＝0．6附近達(dá)到最大值，約1．25，而試驗(yàn)結(jié)果則平緩下降。這可能是由于一方面存在試驗(yàn)誤差，如粒子休止角誤差以及模擬參數(shù)誤差等；另一方面由于觀測時(shí)間持續(xù)較長，受環(huán)境變化影響大，如風(fēng)速、風(fēng)向、溫度以及濕度等的改變。此外，為了防止模型被積雪掩埋，觀測模型與地面有一定距離（0．5m），因此屋面積雪不可避免受下洗流動(dòng)影響，在屋面交界處形成較強(qiáng)的回流渦，造成積雪突然減小。

圖7為10°單跨雙坡屋面3次降雪模擬試驗(yàn)結(jié)果。由圖可見，迎風(fēng)屋面系數(shù)較小，在x＝－0．2附近最大值約0．5，并隨時(shí)間發(fā)展進(jìn)一步減??；屋面最大積雪系數(shù)在背風(fēng)屋面靠屋脊后方，最大系數(shù)約為1．2，并隨時(shí)間發(fā)展有后移趨勢；在背風(fēng)屋面屋檐附件，試驗(yàn)1和2的結(jié)果接近，最大積雪系數(shù)值約為0．9，但試驗(yàn)3的結(jié)果較小，值約為0．5，這可能是由于模擬降雪速度差別帶來的影響。與規(guī)范值相比，10°單跨雙坡屋面的積雪系數(shù)分布形態(tài)與規(guī)范值基本一致，并在規(guī)范值之內(nèi)。

圖6 階梯形屋面積雪系數(shù)Fig．6 Snow coefficients distribution on stepped roof

圖7 10°單跨雙坡屋面積雪系數(shù)Fig．7 Snow coefficients distribution on 10°gable roof

圖8 為10°雙跨雙坡屋面在2個(gè)時(shí)間段內(nèi)的降雪模擬試驗(yàn)結(jié)果。由圖可見，隨時(shí)間發(fā)展，屋面大部分積雪分布變化很小，表明已接近穩(wěn)定狀態(tài)；在x＜0時(shí)屋面積雪系數(shù)相對較大，第1個(gè)迎風(fēng)斜坡屋面（x＜－1）最大積雪在屋面中部，積雪系數(shù)接近規(guī)范值1；在x＝0附近由于屋面凹陷，運(yùn)動(dòng)粒子大多沉積于此，積雪量最多，最大系數(shù)接近規(guī)范的情況1。規(guī)范的情況2在此處取最大值2，安全系數(shù)相對較高，分布形態(tài)與試驗(yàn)結(jié)果也較吻合。

圖9顯示了2種單跨雙坡屋面均勻重分布試驗(yàn)結(jié)果及中國荷載規(guī)范值。由圖可看出，10°單跨雙坡屋面，在迎風(fēng)屋面前部（x＜－0．5）積雪系數(shù)較小，而迎風(fēng)屋面中后部局部系數(shù)較大，最大值接近1，超過了中國荷載規(guī)范值的0．75，但隨著時(shí)間發(fā)展有所減小。這是由于在開始時(shí)，迎風(fēng)屋面前部粒子受風(fēng)力侵蝕不斷后移，使屋面后部積雪系數(shù)初始增加，但隨著遷移粒子數(shù)減少，屋面后部受風(fēng)力侵蝕逐漸增大，積雪系數(shù)有所減小。背風(fēng)屋面靠屋檐附近積雪系數(shù)值最大，接近規(guī)范值的1．25。20°單跨雙坡屋面與10°屋面類似，在迎風(fēng)屋面中后部，最大積雪系數(shù)接近1，大于荷載規(guī)范值；最大積雪系數(shù)在背風(fēng)屋面的屋脊后方，最大值接近規(guī)范的1．25。可以看出，試驗(yàn)的2個(gè)時(shí)段內(nèi)背風(fēng)屋面積雪系數(shù)變化較小，表明此處粒子分布已趨于穩(wěn)定狀態(tài)。

圖8 10°雙跨雙坡屋面積雪系數(shù)Fig．8 Snow coefficients distribution on 10°dual－span gable roof

圖9 2種單跨雙坡屋面積雪系數(shù)Fig．9 Snow coefficients distribution on gable roofs

圖10 顯示了2種雙跨雙坡屋面均勻重分布試驗(yàn)結(jié)果及中國荷載規(guī)范值。由圖可見，10°雙跨雙坡屋面的第1個(gè)迎風(fēng)屋面（x＜－1）靠近屋脊處，積雪系數(shù)相對較大，接近規(guī)范值1，但隨著時(shí)間發(fā)展明顯減??；屋面最大積雪系數(shù)在第1個(gè)背風(fēng)屋面的屋脊后方（－1＜x＜－0．5），最大系數(shù)超過規(guī)范值的1．4，而在屋面中間凹陷處積雪系數(shù)變化平緩，為1．1左右，這與降雪模擬試驗(yàn)結(jié)果有差別；在第2個(gè)背風(fēng)屋面（x＞1），最大積雪系數(shù)接近規(guī)范值1，但隨著時(shí)間發(fā)展明顯減小。20°雙跨雙坡屋面的第1個(gè)迎風(fēng)屋面（x＜－1）初始積雪系數(shù)相對較大，最大系數(shù)接近規(guī)范值1，并隨著時(shí)間發(fā)展略有減小；同10°屋面類似，屋面最大積雪系數(shù)在第1個(gè)背風(fēng)屋面的屋脊后方（－1＜x＜－0．5），初始時(shí)積雪系數(shù)在規(guī)范值以內(nèi)，約為1．3，但隨著時(shí)間發(fā)展，屋面粒子重分布后，最大積雪系數(shù)超過規(guī)范值，達(dá)到了1．5，并沿流向迅速減小，到第2個(gè)迎風(fēng)屋面（0＜x＜1）處，粒子幾乎完全被侵蝕，積雪系數(shù)為0；在第2個(gè)背風(fēng)屋面的屋脊后方（1＜x＜1．5），初始屋面積雪系數(shù)接近0．9，但隨后也很快減小。

圖10 兩種雙跨雙坡屋面積雪系數(shù)Fig．10 Snow coefficients distribution on dual－span gable roofs

3 結(jié) 論

（1）在邊界層流場測量中，當(dāng)達(dá)到粒子起動(dòng)的閾值風(fēng)速后，有效氣動(dòng)粗糙長度隨著風(fēng)速增加而增大，顯示出風(fēng)吹雪邊界層流場典型特性。

（2）階梯形屋面模擬降雪試驗(yàn)，最大積雪系數(shù)在屋面交界處，積雪系數(shù)分布形態(tài)接近規(guī)范的情況1，但與觀測值存在誤差；單跨雙坡屋面的最大積雪系數(shù)在背風(fēng)屋面，但在中國荷載規(guī)范值之內(nèi)，在均勻重分布試驗(yàn)中，2種單跨雙坡屋面的迎風(fēng)屋面積雪系數(shù)均超過了中國荷載規(guī)范值。

（3）雙跨雙坡屋面第1個(gè)迎風(fēng)屋面積雪系數(shù)相對較大，但均不大于荷載規(guī)范值，并隨時(shí)間發(fā)展略有減?。?0°雙跨雙坡屋面模擬降雪試驗(yàn)，最大系數(shù)在屋面中部凹陷位置，與荷載規(guī)范的情況2較符合。重分布試驗(yàn)中，2種屋面的最大系數(shù)均在第1個(gè)背風(fēng)屋面的屋脊后方，值接近1．5，超過了荷載規(guī)范值；20°雙跨雙坡屋面后半部分屋面積雪系數(shù)較小。

（4）采用本文的風(fēng)洞試驗(yàn)方法對屋面雪荷載分布預(yù)測，可得到屋面最大雪荷載分布位置，從而為結(jié)構(gòu)荷載設(shè)計(jì)提供參考。

［1］Feld J，Carper K．Construction Failure［M］．2nd ed．New York：John Wiley ＆Sons，1997．

［2］Perski Z，Leijen F van，Hanssen R．Applicability of psinsar for building hazard identification：Study of the 29 January 2006 Katowice Exhibition Hall collapse and the 24 February 2006 Moscow Basmanny Market collapse［C］／／Proceedings of the Envisat Symposium 2007，Monteux，Switzerland，2007．

［3］Lomovtsev E，Grishin A．Market collapse－Moscow′s largest－ever man－made disaster［J］．Current Digest of the Russian Press，2006，58（9）：4－5．

［4］王元清，胡宗文，石永久，等．門式剛架輕型房屋鋼結(jié)構(gòu)雪災(zāi)事故分析與反思［J］．土木工程學(xué)報(bào)，2009，42（3）：65－70．

Wang Y Q，Hu Z W，Shi Y J，et al．Analysis and reflection on snow disaster accidents of steel structures of light－weight buildings with portal frames［J］．China Civil Engineering Journal，2009，42（3）：65－70．

［5］藍(lán)聲寧，鐘新谷．湘潭輕型鋼結(jié)構(gòu)廠房雪災(zāi)受損分析與思考［J］．土木工程學(xué)報(bào)，2009，42（3）：71－75．

Lan S N，Zhong X G．Damage diagnoses and lessons learnt from the failure of light－steel structure by heavy snow in Xiangtan［J］．China Civil Engineering Journal，2009，42（3）：71－75．

［6］周昌農(nóng)．南方雪災(zāi)中輕鋼結(jié)構(gòu)房屋災(zāi)情調(diào)查及原因分析［J］．城市建設(shè)，2009，32：286－287．

Zhou C N．Damage diagnoses and lessons learnt from the failure of light－steel structure by heavy snow in South China［J］．Chengshi Jianshe Yu Shangye Wangdian，2009，32（46）：286－287．

［7］張德海，南波，舒錚．雪災(zāi)后某網(wǎng)架破壞分析及加固［J］．沈陽建筑大學(xué)學(xué)報(bào)（自然科學(xué)版），2010，26（1）：62－67．

Zhang D H，Nan B，Shu Z．Destroying analysis and reinforcement research on space truss structure after snowstorms［J］．Journal of shenyang Jianzhu University（Natural Science），2010，26（1）：62－67．

［8］顧明，黃友欽，趙明偉．風(fēng)雪共同作用下門式剛架廠房的動(dòng)力穩(wěn)定性［J］．同濟(jì)大學(xué)學(xué)報(bào)，2011，39（9）：1266－1272．

Gu M，Huang Y Q，Zhao M W．Dynamic instability of lightweight steel structures with portal frames under wind and snow loads［J］．Journal of Tongji University，2011，39（9）：1266－1272．

［9］Tsuchiya M，Tomabechi T，Hongoa T，et al．Wind effects on snowdrift on stepped flat roofs［J］．Journal of Wind Engineering and Industrial Aerodynamics，2002，90：1881－1892．

［10］Gordon M，Taylor P A．Measurements of blowing snow，Part I：Particle shape，size distribution，velocity，and number flux at Churchill，Manitoba，Canada［J］．Cold Regions Science and Technology，2009，55（1）：63－74．

［11］O′Rourke M，Degaetano A，Tokarczyk J D．Analytical simulation of snow drift loading［J］．Journal of Structural Engineering，2005，131：660－667．

［12］Tominaga Y，Okaze T，Mochida A．CFD modeling of snowdrift around a building：An overview of models and evaluation of a new approach［J］．Building and Environment，2011，46（4）：899－910．

［13］王衛(wèi)華，廖海黎，李明水．基于時(shí)變邊界屋面積雪分布數(shù)值模擬［J］．西南交通大學(xué)學(xué)報(bào)，2013，48（5）：851－967．

Wang W H，Liao H L，Li M S．Numerical simulation of windinduced roof snow distributions based on time variable boundary［J］．Journal of Southwest Jiao－tong University，2013，48（5）：851－967．

［14］O′Rourke M，Degaetano A，Tokarczyk J D．Snow drifting transport rates from water flume simulation［J］．Journal of Wind Engineering and Industrial Aerodynamics，2004，92：1245－1264．

［15］Thiis T K，Barfoed P，Delpech P，et al．Penetration of snow into roof constructions—Wind tunnel testing of different eave cover designs［J］．Journal of Wind Engineering and Industrial Aerodynamics，2007，95（9－11）：1476－1485．

［16］Kimbar G，F(xiàn)laga A．Wind tunnel tests of snow load redistribution on large span flat roofs［C］．The 13th International Conference on Wind Engineering，Amsterdam，2011．

［17］王衛(wèi)華，廖海黎，李明水．風(fēng)致屋面積雪分布風(fēng)洞試驗(yàn)研究［J］．建筑結(jié)構(gòu)學(xué)報(bào)，2014，35（5）：143－149．

Wang W H，Liao H L，Li M S．Wind tunnel test on wind－induced roof snow distribution［J］．Journal of Building Structures，2014，35（5）：143－149．

［18］Zhou X Y，Kang L Y，Yuan X M，et al．Wind tunnel test of snow redistribution on flat roofs［J］．Cold Regions Science and Technology，2016，127：49－56．

［19］中華人民共和國建設(shè)部．GB 50009－2012建筑結(jié)構(gòu)荷載規(guī)范［S］．北京：中國建筑工業(yè)出版社，2012．

［20］Kind R J．Snowdrifting：a review of modeling methods［J］．Cold Regions Science and Technology，1986，12（3）：217－228．

［21］Beyers J H M，Harms T M．Outdoors modelling of snowdrift at SANAE IV Research Station，Antarctica［J］．Journal of Wind Engineering and Industrial Aerodynamics，2003，91：551－569．

Study of snow-load distribution on roof by w ind tunnel test

Wang Weihua*，Huang Hanjie
(China Aerodynamics Research and Development Center，Mianyang Sichuan 621000，China)

In order to predict the snow-load distribution on a roof，wind tunnel testswere conducted with fine quartz-sand to simulate snow particles．The snow coefficients of several typical roofs were obtained through snowfall simulation and redistribution tests，and were compared with the load design of building structures．The results show that，when the wind speed exceeds the threshold value，the efficient aerodynamic-rough-length increases as the wind speed increases;the snow coefficients of the stepped roof obtained in the tests are within the load specifications;for the single-span gable roof，the snow coefficients of the windward roofmay exceed the load specifications，particularly for the 20°roof，while the 10°roof snow coefficients decrease with time and finally fall into the range of load specifications;for the double-span gable roof，in snowfall simulation tests，the largest coefficient lies on the firstwindward roof，the value ofwhich is less than the load specifications，while in the simulation redistributions test，the largest coefficient lies on the first leeward roof，the value is approximately 1．5，larger than the load specifications．

building;roof;snow-load;wind tunnel test

V211．73

A

（編輯：李金勇）

1672-9897(2016)05-0023-06

10．11729/syltlx20160039

2016-03-07;

2016-06-03

*通信作者E-mail:wwhblue@163．com

Wang W H，Huang H J．Study of snow-load d istribution on roof by wind tunnel test．Journalof Experiments in Fluid Mechanics，2016，30(5):23-28．王衛(wèi)華，黃漢杰．屋面雪荷載分布風(fēng)洞試驗(yàn)研究．實(shí)驗(yàn)流體力學(xué)，2016，30(5):23-28．

王衛(wèi)華（1972－），男，安徽合肥人，博士，高級工程師。研究方向：結(jié)構(gòu)風(fēng)工程。通信地址：四川省綿陽市二環(huán)路南段6號（621000）。E－mail：wwhblue＠163．com