余香林， 董　銳， 王守強(qiáng)

(1. 同濟(jì)大學(xué) 土木工程學(xué)院， 上海 200092； 2. 福州大學(xué) 土木工程學(xué)院， 福建 福州 350108；3. 同濟(jì)大學(xué) 土木工程防災(zāi)國(guó)家重點(diǎn)實(shí)驗(yàn)室， 上海 200092)

?

典型屋頂太陽(yáng)能光伏板風(fēng)壓實(shí)測(cè)研究

余香林1， 董銳2,3， 王守強(qiáng)3

(1. 同濟(jì)大學(xué) 土木工程學(xué)院， 上海 200092； 2. 福州大學(xué) 土木工程學(xué)院， 福建 福州 350108；3. 同濟(jì)大學(xué) 土木工程防災(zāi)國(guó)家重點(diǎn)實(shí)驗(yàn)室， 上海 200092)

摘要：對(duì)一棟多功能低層建筑平屋頂上的太陽(yáng)能光伏板進(jìn)行風(fēng)壓實(shí)測(cè)，根據(jù)風(fēng)壓時(shí)程特性選擇適當(dāng)?shù)母怕拭芏群瘮?shù)對(duì)各測(cè)點(diǎn)的風(fēng)壓系數(shù)進(jìn)行非線性回歸分析.運(yùn)用6種概率密度函數(shù)對(duì)各測(cè)點(diǎn)的均值和極值(正極大值和負(fù)極小值)風(fēng)壓系數(shù)時(shí)程以及面積平均后的極值風(fēng)壓系數(shù)時(shí)程進(jìn)行數(shù)值擬合，比較不同概率分布的擬合效果.回歸分析結(jié)果表明：均值風(fēng)壓系數(shù)的概率分布近似無(wú)偏分布，t分布的擬合效果最好，其次是Logistic分布、正態(tài)分布、極值I型分布；負(fù)壓極小值和正壓極大值風(fēng)壓系數(shù)的概率分布分別為左偏分布和右偏分布，極值I型分布的擬合效果最好，其次是Lognormal分布和Gamma分布，而正態(tài)分布的擬合效果最差；最不利風(fēng)壓系數(shù)由極值負(fù)風(fēng)壓控制，即光伏板的風(fēng)力由吸力主導(dǎo)；經(jīng)面積平均后的最不利風(fēng)壓系數(shù)建議取值-2.3.

關(guān)鍵詞：平屋頂； 太陽(yáng)能光伏板； 原型試驗(yàn)； 風(fēng)壓系數(shù)； 極值分布

屋頂太陽(yáng)能光伏板在住宅和辦公建筑中日益興起.作用于光伏板上的風(fēng)壓具有很強(qiáng)的脈動(dòng)性和持續(xù)性，其風(fēng)壓系數(shù)極值和均值的概率分布是確定工程設(shè)計(jì)用風(fēng)壓系數(shù)大小的基礎(chǔ).普通房屋結(jié)構(gòu)風(fēng)壓時(shí)程的概率分布已有不少學(xué)者進(jìn)行過(guò)研究，且將風(fēng)壓系數(shù)編入了相關(guān)荷載規(guī)范[1-3].Holmes等[4]運(yùn)用極值I型分布和廣義極值分布對(duì)風(fēng)壓系數(shù)極值進(jìn)行了擬合，發(fā)現(xiàn)具有較小正形狀參數(shù)的廣義極值分布擬合效果最好.Tieleman等[5]用三參數(shù)Gamma分布很好地?cái)M合了偏度較小的風(fēng)壓系數(shù)時(shí)程.陶玲等[6]得出低矮房屋風(fēng)壓時(shí)程的概率分布與偏度有關(guān)，偏度增大時(shí)分別接近對(duì)數(shù)正態(tài)分布、廣義極值分布，而面積平均后的風(fēng)壓時(shí)程與對(duì)數(shù)正態(tài)分布、Gamma分布和廣義極值分布都吻合較好.由于現(xiàn)場(chǎng)實(shí)測(cè)難度較大，屋頂太陽(yáng)能光伏板的風(fēng)壓分布規(guī)律研究一般通過(guò)風(fēng)洞試驗(yàn)進(jìn)行.操金鑫等[7]對(duì)平屋頂光伏板的凈風(fēng)壓系數(shù)進(jìn)行了風(fēng)洞試驗(yàn)，結(jié)果表明單排電池板的最不利負(fù)風(fēng)壓系數(shù)遠(yuǎn)大于多排布置的情況，且隨電池板傾角和排距的增大而增大.Stathopoulos等[8]通過(guò)風(fēng)洞試驗(yàn)得出光伏板最不利風(fēng)壓系數(shù)由負(fù)壓主導(dǎo)，并建議面積平均后的最大凈風(fēng)壓系數(shù)取-1.7.

盡管日本規(guī)范[9]對(duì)作用于光伏板上的風(fēng)荷載作了明確規(guī)定，但國(guó)內(nèi)荷載規(guī)范[2]并無(wú)此類(lèi)規(guī)定.由于不同地區(qū)的光伏板支撐系統(tǒng)和布局形式存在不同，導(dǎo)致其風(fēng)荷載分布也存在區(qū)別.鑒于此，有必要對(duì)其他地區(qū)光伏板的風(fēng)壓分布進(jìn)行研究，以便保證結(jié)構(gòu)的設(shè)計(jì)安全.風(fēng)洞試驗(yàn)、現(xiàn)場(chǎng)實(shí)測(cè)和CFD數(shù)值模擬是結(jié)構(gòu)抗風(fēng)研究的三種重要手段，其中現(xiàn)場(chǎng)實(shí)測(cè)是目前公認(rèn)的相對(duì)最準(zhǔn)確的研究方法.本文采用現(xiàn)場(chǎng)原型試驗(yàn)對(duì)某典型太陽(yáng)能光伏板風(fēng)壓分布進(jìn)行了研究，并從統(tǒng)計(jì)學(xué)和結(jié)構(gòu)抗風(fēng)設(shè)計(jì)可靠性的角度分析了其風(fēng)壓分布規(guī)律和風(fēng)壓系數(shù)的取值標(biāo)準(zhǔn)，得到了許多有意義的結(jié)論，可為我國(guó)工程設(shè)計(jì)和相關(guān)規(guī)范的制定提供借鑒和參考.

1試驗(yàn)概況

本試驗(yàn)的研究對(duì)象是某低矮建筑平屋頂上的太陽(yáng)能光伏板，距離室外地面的高度為12 m，沿屋頂縱向排布.為充分吸收太陽(yáng)能，光伏板與平屋頂?shù)膴A角為25°.為使現(xiàn)場(chǎng)實(shí)測(cè)數(shù)據(jù)具有代表性，取靠近山墻一側(cè)的6排4列共計(jì)24小塊光伏板(上下兩大塊)作為觀測(cè)對(duì)象，太陽(yáng)能光伏板試驗(yàn)原型如圖1所示.該低矮建筑位于房屋較稀疏的郊區(qū)，場(chǎng)地類(lèi)別為B類(lèi)[2].

圖1　屋頂太陽(yáng)能光伏板實(shí)物圖

在24塊光伏板上布置12個(gè)風(fēng)壓傳感器，每塊光伏板的面積均為1.3 m2(1.58 m×0.83 m).在光伏板頂部靠左位置安裝型號(hào)85 000超聲風(fēng)速儀.光伏板尺寸及風(fēng)壓傳感器測(cè)點(diǎn)布置如圖2所示.

風(fēng)壓傳感器測(cè)得的風(fēng)壓為凈風(fēng)壓，即正負(fù)壓差.光伏板迎風(fēng)面和背風(fēng)面均設(shè)有測(cè)壓孔，通過(guò)測(cè)量正反面的壓差，將測(cè)得的凈風(fēng)壓數(shù)據(jù)傳輸至基于LabVIEW2012的數(shù)據(jù)采集系統(tǒng)中.傳感器在光伏板上的布置如圖3所示.

圖2　太陽(yáng)能電板傳感器布置圖(單位：m)

圖3　風(fēng)壓傳感器布置圖

超聲風(fēng)速儀可測(cè)得實(shí)時(shí)風(fēng)速大小和風(fēng)向.因太陽(yáng)能光伏板具有25°傾角，且風(fēng)速儀布置在光伏板頂部檐口處；此外，風(fēng)速儀的總高度為0.347 m，安裝深度為0.114 m，故風(fēng)速儀探頭距離屋頂表面的豎直高度為5.59×sin25°+0.347-0.114=2.6 m.因此，風(fēng)速儀所測(cè)風(fēng)速和風(fēng)向受建筑物干擾影響很小，可以代表來(lái)流風(fēng)速和風(fēng)向.規(guī)定北風(fēng)(自北朝南吹的風(fēng)) 為0°風(fēng)，東風(fēng) (自東向西吹的風(fēng)) 為90°風(fēng)，沿順時(shí)針?lè)较蜃兓瘯r(shí)風(fēng)向角逐漸增大，變化范圍在0°～360°之間.風(fēng)速儀實(shí)物圖和俯視圖見(jiàn)如圖4和圖5.

圖4　型號(hào)85 000超聲風(fēng)速儀實(shí)物圖

圖5型號(hào)85 000超聲風(fēng)速儀俯視圖

Fig.5Top view of Model 85 000 Ultrasonic Anemometer

風(fēng)速、風(fēng)向及風(fēng)壓等試驗(yàn)數(shù)據(jù)的采樣頻率為100 Hz，每個(gè)樣本的采樣持時(shí)為15 min，樣本容量為90 000，樣本數(shù)量共計(jì)501個(gè).為防止灰塵、雨水等滲入風(fēng)速儀和風(fēng)壓傳感器的測(cè)壓孔中，每次采集數(shù)據(jù)前，均留有5 min的試驗(yàn)準(zhǔn)備時(shí)間，用來(lái)檢查各測(cè)量?jī)x器及線路板是否正常工作.

2光伏板風(fēng)壓系數(shù)分析

因單個(gè)樣本數(shù)據(jù)量較大(90 000×14=1.26×106)，為排除外界環(huán)境等對(duì)實(shí)測(cè)風(fēng)壓和風(fēng)速造成的干擾，運(yùn)用滑動(dòng)平均法對(duì)數(shù)據(jù)進(jìn)行平滑處理(滑動(dòng)平均窗口為50)，再運(yùn)用低通濾波器對(duì)高頻成分進(jìn)行濾波處理，截止頻率取15Hz.風(fēng)壓系數(shù)的計(jì)算公式為

(3)

2.1平均風(fēng)壓系數(shù)

平均風(fēng)壓系數(shù)的概率分布近似對(duì)稱(chēng)分布，故采用正態(tài)分布，t分布和Logistic分布擬合.

2.1.1正態(tài)分布

正態(tài)分布屬無(wú)偏分布，其概率密度函數(shù)(PDF)和相應(yīng)的累積分布函數(shù)(CDF)分別為[10]

-∞

(4)

(5)

式中：σX和μX表示隨機(jī)變量X的標(biāo)準(zhǔn)差和均值.

正態(tài)分布的均值凈風(fēng)壓系數(shù)CPmean參數(shù)回歸分析見(jiàn)表1.取雙側(cè)2.5%分位點(diǎn)對(duì)應(yīng)的風(fēng)壓系數(shù)作為該測(cè)點(diǎn)的平均凈風(fēng)壓系數(shù).

表1　正態(tài)分布的平均風(fēng)壓系數(shù)CPmean參數(shù)分析

2.1.2t分布

t分布屬于抽樣分布，具有厚尾特性.當(dāng)樣本容量足夠大時(shí)，t分布近似于正態(tài)分布.因CPmean的均值不為0(接近于0)，故用三參數(shù)t分布進(jìn)行回歸分析.其PDF和CDF分別為[11]

(6)

(7)

式中:Γ為Gamma函數(shù)；a為位置參數(shù);b為尺度參數(shù);c為自由度數(shù).t分布的均值凈風(fēng)壓系數(shù)CPmean參數(shù)回歸分析見(jiàn)表2.表中，R2為拉入合度.

表2　t分布的平均風(fēng)壓系數(shù)CPmean參數(shù)分析

由表2知，每塊光伏板上排測(cè)點(diǎn)(S10, S11, S12和S4, S5, S6)的CPmean的均值在左側(cè)2.5%分位點(diǎn)最小(絕對(duì)值最大)，分別為-0.98和-1.53.

2.1.3Logistic分布

Logistic分布屬于增長(zhǎng)型分布，亦具有厚尾特性.其PDF和CDF分別為[11]

-∞

(8)

(9)

式中:α是位置參數(shù)，β是尺度參數(shù).Logistic分布的CDF稱(chēng)為L(zhǎng)ogistic函數(shù)，且通常用于Logistic回歸分析中.Logistic分布的均值凈風(fēng)壓系數(shù)CPmean參數(shù)回歸分析見(jiàn)表3.

由表3可知，雙側(cè)2.5%分位點(diǎn)對(duì)應(yīng)的CPmean的絕對(duì)值比正態(tài)分布情況下的絕對(duì)值稍大，而比t分布情況下的絕對(duì)值小很多，表明Logistic分布的厚尾程度介于正態(tài)分布和t分布之間.均值風(fēng)壓系數(shù)CPmean擬合度對(duì)比如圖6所示.

由圖6知，t分布的擬合度最高，其次是Logistic分布和正態(tài)分布，而極值I型分布的擬合度最低.

2.2極值風(fēng)壓系數(shù)

極值風(fēng)壓系數(shù)分為正極大值CPmax和負(fù)極小值CPmin兩種.根據(jù)極值風(fēng)壓的分布特點(diǎn)，選取常見(jiàn)的概率密度函數(shù)(極值I型分布, Gamma分布, 對(duì)數(shù)正態(tài)分布, 正態(tài)分布)對(duì)各測(cè)點(diǎn)的極值風(fēng)壓系數(shù)進(jìn)行參數(shù)回歸分析.

圖6　均值風(fēng)壓系數(shù)CPmean的擬合度對(duì)比

2.2.1極值I型分布

通過(guò)分析各測(cè)點(diǎn)極值風(fēng)壓系數(shù)的概率質(zhì)量函數(shù)(PMF, 用直方圖表示)及統(tǒng)計(jì)量(見(jiàn)表4)可知，CPmax為右偏分布，CPmin為左偏分布(S7測(cè)點(diǎn)例外).

表4　極值風(fēng)壓系數(shù)的統(tǒng)計(jì)量

極值I型(Gumbel)分布適用于支撐集為R的右偏分布(偏度為1.14).其PDF和CDF分別 為[10]

(10)

(11)

式中：α是尺度參數(shù)，且α>0；β是對(duì)應(yīng)于概率密度最大值的變量X的眾值.α，β用X的均值和標(biāo)準(zhǔn)差表示為

(12)

圖7　左偏分布的模型變換

通過(guò)對(duì)極值凈風(fēng)壓系數(shù)CPmax和CPmin進(jìn)行參數(shù)回歸分析可知，頂部整塊光伏板上排測(cè)點(diǎn)(S10, S11, S12)的CPmax的均值最大，為0.4.每大塊光伏板頂排測(cè)點(diǎn)(S10, S11, S12和S4, S5, S6) 的CPmin均值的絕對(duì)值最大，分別為1.75和2.4.

2.2.2Gamma分布

(13)

(14)

(15)

對(duì)于左偏分布且x<0時(shí)，令y=-x>0，可將模型對(duì)稱(chēng)變換為右偏分布.

通過(guò)對(duì)極值凈風(fēng)壓系數(shù)CPmax和CPmin進(jìn)行參數(shù)回歸分析可知，頂部整塊光伏板上排測(cè)點(diǎn)(S10, S11, S12)的CPmax的均值最大，為0.45.每大塊光伏板頂排測(cè)點(diǎn)(S10, S11, S12和S4, S5, S6) 的CPmin的均值的絕對(duì)值最大，分別為2.07和2.75.

2.2.3Lognormal分布

兩參數(shù)Lognormal分布(即對(duì)數(shù)正態(tài)分布)的隨機(jī)變量X的自然對(duì)數(shù)InX服從正態(tài)分布，適用于右偏分布.其隨機(jī)變量X的PDF和CDF分別為[10]

0

(16)

(17)

偏度為

(18)

式中：λX是位置參數(shù)；ζX是尺度參數(shù).λX和ζX可用X的均值和標(biāo)準(zhǔn)差表示為

(19)

通過(guò)對(duì)極值凈風(fēng)壓系數(shù)CPmax和CPmin進(jìn)行參數(shù)回歸分析可知，頂排整塊光伏板上排測(cè)點(diǎn)(S10, S11, S12)的CPmax的均值最大，為0.62.每大塊光伏板頂排測(cè)點(diǎn)(S10, S11, S12和S4, S5, S6)的CPmin的絕對(duì)值最大，分別為-3.54和-2.97，均大于極值I型分布和Gamma分布對(duì)應(yīng)的數(shù)值，這是由Lognormal分布的厚尾特性引起的.極值風(fēng)壓系數(shù)的擬合度對(duì)比如圖8和圖9所示.

由圖8和圖9知，Lognormal,Gamma和極值I型分布的擬合度均較高，正態(tài)分布的擬合度最低.

對(duì)上下兩塊光伏板的各測(cè)點(diǎn)極值風(fēng)壓系數(shù)時(shí)程進(jìn)行實(shí)時(shí)平均，得到面積平均后的正極大值和負(fù)極小值風(fēng)壓系數(shù).各種分布的面積平均后極值風(fēng)壓系數(shù)擬合度對(duì)比如表5所示.概率密度和面積平均后的概率質(zhì)量(直方圖) 比較如圖10～圖13所示.

圖8　正極大值風(fēng)壓系數(shù)CPmax的擬合度對(duì)比

圖9　負(fù)極小值風(fēng)壓系數(shù)CPmin的擬合度對(duì)比

表5　各種分布的面積平均后極值風(fēng)壓系數(shù)擬合度對(duì)比

圖10　面積平均后頂部光伏板的CPmax的概率密度對(duì)比

Fig.10Comparison of probability density for area-averaged maximum pressure coefficientCPmaxof top PV panels

圖11　面積平均后底部光伏板的CPmax的概率密度對(duì)比

Fig.11Comparison of probability density for area-averaged maximum pressure coefficientCPmaxof bottom PV panels

圖12　面積平均后頂部光伏板的CPmin的概率密度對(duì)比

Fig.12Comparison of probability density for area-averaged minimum pressure coefficientCPminof top PV panels

Fig.13Comparison of probability density for area-averaged minimum pressure coefficientCPminof bottom PV panels

由表5可知，除了頂部整塊光伏板的負(fù)極小值風(fēng)壓系數(shù)的擬合度較低外，其他極值風(fēng)壓系數(shù)的擬合度均較高.頂部光伏板負(fù)極小值風(fēng)壓系數(shù)擬合度偏低是因S7測(cè)點(diǎn)的分布為極值I型分布 (擬合度為93%，運(yùn)用正態(tài)分布的擬合度為86%，概率密度和概率質(zhì)量如圖14所示)，且CPmin有正有負(fù)，而Gamma分布和Lognormal分布僅適用于自變量X均為正數(shù)或均為負(fù)數(shù)(經(jīng)Y=-X的對(duì)稱(chēng)變換后，將Y作為新的自變量即為正數(shù))的情況，故不適于用Gamma分布和Lognormal分布擬合.

圖14　測(cè)點(diǎn)S7負(fù)極小值風(fēng)壓系數(shù)CPmin的概率密度對(duì)比

Fig.14Comparison of probability density for negative minimum peak pressure coefficientCPminof S7

3結(jié)果分析

S7的負(fù)極小值風(fēng)壓系數(shù)不宜使用Lognormal分布和Gamma分布擬合，且面積平均后的風(fēng)壓系數(shù)擬合度較極值I型分布更低，故從概率密度函數(shù)的整體擬合效果來(lái)看，極值風(fēng)壓系數(shù)時(shí)程最適合用極值I型分布擬合，而均值凈風(fēng)壓系數(shù)最適合用t分布.基于這兩種分布的各測(cè)點(diǎn)風(fēng)壓系數(shù)分析結(jié)果如表6所示.

表6　各測(cè)點(diǎn)凈風(fēng)壓系數(shù)分析結(jié)果

由表6可知，各測(cè)點(diǎn)的最大凈風(fēng)壓系數(shù)(絕對(duì)值最大)由負(fù)極小值凈風(fēng)壓系數(shù)控制，且均小于0，表明作用于光伏板上的風(fēng)力以吸力為主.因12個(gè)風(fēng)壓測(cè)點(diǎn)均勻分布在24小塊光伏板的典型位置處，為使每塊從屬面積均包含1個(gè)測(cè)點(diǎn)，故等面積劃分成12塊.各測(cè)點(diǎn)的負(fù)極小值凈風(fēng)壓系數(shù)及從屬面積劃分如圖15所示.

圖15　各測(cè)點(diǎn)負(fù)極小值凈風(fēng)壓系數(shù)及從屬面積劃分

取各測(cè)點(diǎn)風(fēng)壓系數(shù)的實(shí)測(cè)值乘以對(duì)應(yīng)的從屬面積進(jìn)行加權(quán)平均，得到整塊光伏板上的風(fēng)壓系數(shù).面積A平均后的凈風(fēng)壓系數(shù)CP用下列公式計(jì)算[13]

(20)

則上式簡(jiǎn)化為

(21)

故頂部和底部整塊光伏板的風(fēng)壓系數(shù)分別為

面積平均后的頂部和底部半塊光伏板的風(fēng)壓系數(shù)分布如圖16所示.

由圖16和計(jì)算結(jié)果可知，每塊光伏板上半部分的風(fēng)壓系數(shù)均比下半部分大，且頂部光伏板的風(fēng)壓系數(shù)大于底部光伏板的風(fēng)壓系數(shù).

將各個(gè)測(cè)點(diǎn)風(fēng)壓系數(shù)對(duì)應(yīng)的從屬面積進(jìn)行加權(quán)平均后求得的凈風(fēng)壓系數(shù)進(jìn)行線性回歸分析，且超越概率取5%，可得出如圖17所示的設(shè)計(jì)用三折線凈風(fēng)壓系數(shù)圖.當(dāng)有效受風(fēng)面積為0.1 m2～0.9 m2時(shí)，最不利凈風(fēng)壓系數(shù)為-2.3.故平屋頂太陽(yáng)能光伏板傾角為25°時(shí)，建議荷載規(guī)范對(duì)光伏板凈風(fēng)壓系數(shù)取值為-2.3.

圖16　凈風(fēng)壓系數(shù)分布圖

圖17　太陽(yáng)能光伏板凈風(fēng)壓系數(shù)建議值

4結(jié)論

對(duì)作用于低層建筑平屋頂太陽(yáng)能光伏板上的風(fēng)壓和風(fēng)速時(shí)程進(jìn)行現(xiàn)場(chǎng)實(shí)測(cè)，合理運(yùn)用常見(jiàn)的概率密度函數(shù)對(duì)風(fēng)壓時(shí)程分布進(jìn)行非線性回歸分析，比較選取最適合風(fēng)壓分布規(guī)律的分布函數(shù)，提出規(guī)范化的凈風(fēng)壓系數(shù)分布模型，得出以下結(jié)論：

(1)極值風(fēng)壓系數(shù)符合極值I型分布，且Lognormal分布和Gamma分布的擬合效果均較好.

(2)均值風(fēng)壓系數(shù)時(shí)程近似無(wú)偏分布，且具有厚尾性質(zhì)，故t分布的擬合效果非常好.

(3)最不利風(fēng)壓由負(fù)壓極小值控制，表明作用于光伏板上的風(fēng)壓由吸力主導(dǎo)，且最不利風(fēng)壓系數(shù)均出現(xiàn)在每塊光伏板的上部，這是由于上下兩塊光伏板之間的空隙較大，且頂部光伏板上半部分位于自由風(fēng)檐口處導(dǎo)致產(chǎn)生更大的吸力.底部光伏板因距離平屋頂檐口較遠(yuǎn)，故所受平均吸力明顯小于頂部光伏板.

(4)頂部整塊光伏板的凈風(fēng)壓系數(shù)明顯大于底部整塊光伏板的凈風(fēng)壓系數(shù).建議取最不利風(fēng)壓系數(shù)(即體型系數(shù)) 作為設(shè)計(jì)參考，經(jīng)面積平均后的最不利風(fēng)壓系數(shù)取值為-2.3.

參考文獻(xiàn):

[1]American Society of Civil Engineers. ASCE/SEI 7-10 Minimum design loads for buildings and other structures[S]. Reston: American Society of Civil Engineers, 2010.

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

Ministry of Housing and Urban-Rural Development of the People’s Republic of China. GB50009-2012 Load code for the design of building structures[S]. Beijing: China Architecture & Building Press, 2012.

[3]International Code Council (ICC). 2012 International building code[S]. Country Club Hills: ICC, 2011.

[4]Holmes J D, Cochran L S. Probability distributions of extreme pressure coefficients [J]. Journal of Wind Engineering and Industrial Aerodynamics, 2003, 91: 893.

[5]Tieleman H W, Ge Z F, Haji M R. Theoretically estimated peak wind loads [J]. Journal of Wind Engineering and Industrial Aerodynamics, 2007, 95: 113.

[6]陶玲, 黃鵬, 顧明,等. 低矮房屋風(fēng)壓時(shí)程的概率分布 [J]. 同濟(jì)大學(xué)學(xué)報(bào): 自然科學(xué)版, 2013, 41 (1): 27.

TAO Ling, HUANG Peng, GU Ming,etal. Probability density distribution of wind pressure time series of low-rise buildings [J]. Journal of Tongji University: Natural Science, 2013, 41 (1): 27.

[7]Cao J X, Yoshida A, Saha P K,etal. Wind loading characteristics of solar arrays mounted on flat roofs [J]. Journal of Wind Engineering and Industrial Aerodynamics, 2013, 123: 214.

[8]Stathopoulos T, Zisis I, Xypnitou E. Local and overall wind pressure and force coefficients for solar panels [J]. Journal of Wind Engineering and Industrial Aerodynamics, 2014, 125: 195.

[9]Japanese Industrial Standards Committee. JIS C 8955:2011 Design guide on structures for photovoltaic array[S]. Tokyo: Japanese Standards Association, 2011.

[10]Halder A, Mahadevan S. Probability, reliability and statistical methods in engineering design [M]. New York: John Wiley & Sons, 2000.

[11]McLaughlin M P. Regress+, appendix A: a compendium of common probability distributions, Version 2.3 [M]. Fairfax: McLean, 1999.

[12]Johnson N L, Kotz S, Balakrishnan N. Continuous univariate distributions (Vol.1)[M]. 2nd ed. New York: John Wiley & Sons, 1994.

[13]張相庭. 結(jié)構(gòu)風(fēng)工程: 理論·規(guī)范·實(shí)踐 [M]. 北京: 中國(guó)建筑工業(yè)出版社, 2006.

ZHANG Xiangting. Structural wind engineering: theory, code and practice [M]. Beijing: China Architecture & Building Press, 2006.

Full-scale Measurement and Investigation of Wind Pressure Loadings on Typical Rooftop Photovoltaic Solar Panels

YU Xianglin1, DONG Rui2, 3, WANG Shouqiang3

(1. College of Civil Engineering, Tongji University, Shanghai 200092, China; 2. College of Civil Engineering, Fuzhou University, Fuzhou 350108, China; 3. State Key Laboratory for Disaster Reduction in Civil Engineering, Tongji University, Shanghai 200092, China)

Abstract：A full-scale measurement of wind pressure loadings was conducted on photovoltaic (PV) solar panels mounted on the flat roof of a multipurpose low-rise building. According to the characteristics of wind pressure time series, several probability density functions were properly selected for nonlinear regression analysis of wind pressure coefficients of each gauging point. Numerical fittings of mean and peak (positive maxima and negative minima) and area-averaged peak pressure coefficients were performed and compared by using six probability distributions. The analysis results show that the probability distributions of negative minima and positive maxima of pressure coefficients are left-skewed and right-skewed, respectively. The mean pressure coefficient time series comply with unbiased distribution and the t distribution fits the measured data best, followed by Logistic, Normal, Type I EDV distribution, respectively; Type I EVD has the best fitting effect for peak values, and then the Lognormal distribution, Gamma distribution and Normal distribution in order; The most unfavorable pressure coefficient is dominated by negative peak pressure, that is, suction forces are predominant over wind forces on PV panels; The tributary area-averaged unfavorable pressure coefficient is recommended to take -2.3.

Key words：flat roof; Photovoltaic (PV) solar panel; full-scale experiment; pressure coefficient; extreme value distribution

文獻(xiàn)標(biāo)志碼：A

中圖分類(lèi)號(hào)：TU312+.1; O211.3

基金項(xiàng)目：橋梁結(jié)構(gòu)抗風(fēng)技術(shù)交通行業(yè)重點(diǎn)實(shí)驗(yàn)室開(kāi)放課題基金(KLWRTBMC14-03) ; 福建省自然科學(xué)基金(2015J05098)

收稿日期：2015—05—21

第一作者： 余香林(1984—)，男，博士生,主要研究方向?yàn)榻Y(jié)構(gòu)抗風(fēng)和穩(wěn)定. E-mail: 6yux@#edu.cn