雷 克，張其林，龐紹華

（1.同濟(jì)大學(xué) 土木工程學(xué)院，上海 200092；2.上海同磊土木工程技術(shù)有限公司，上海 200092）

建筑玻璃幕墻是指由玻璃和支承結(jié)構(gòu)體系共同組成的建筑外圍護(hù)結(jié)構(gòu)，按照支承體系的不同可分為點(diǎn)支、全玻式和框式玻璃幕墻等［1］.玻璃幕墻已成為現(xiàn)代建筑的標(biāo)志.伴隨著能源危機(jī)的加重和中國(guó)城市化進(jìn)程的推進(jìn)，建筑能耗占我國(guó)社會(huì)總能耗的比例呈增長(zhǎng)態(tài)勢(shì)，且增速越來(lái)越快，從1996年的19%增加到了2008年的23%［2］.而建筑能耗中的很大一部分用于改善室內(nèi)的熱環(huán)境［3］.建筑幕墻作為現(xiàn)代建筑中使用較多的圍護(hù)結(jié)構(gòu)，其熱工性能的優(yōu)劣將直接影響室內(nèi)熱環(huán)境，從而決定建筑能耗.

國(guó)外關(guān)于建筑玻璃幕墻熱工性能的研究起步較早［4-7］，國(guó)內(nèi)學(xué)者也對(duì)玻璃幕墻熱工性能計(jì)算進(jìn)行了相關(guān)研究［8-10］.我國(guó)在參考國(guó)內(nèi)外幕墻熱工計(jì)算研究成果和標(biāo)準(zhǔn)體系的基礎(chǔ)上，頒布了《建筑門(mén)窗玻璃幕墻熱工計(jì)算規(guī)程》（JGJ／T151—2008）［11］，包括了玻璃幕墻熱工性能的相關(guān)計(jì)算方法，但關(guān)于計(jì)算理論未詳細(xì)說(shuō)明.

本文對(duì)框式玻璃幕墻傳熱機(jī)理進(jìn)行了研究，基于傳熱學(xué)理論，運(yùn)用差分法和有限單元法分別推導(dǎo)并建立了玻璃體系和框體系的熱工分析數(shù)值計(jì)算模型，并采用Visual C＋＋和ObjectARX.對(duì)AutoCAD進(jìn)行了二次開(kāi)發(fā)，研發(fā)了適用于我國(guó)規(guī)程的框式玻璃幕墻熱工性能分析計(jì)算軟件.研究結(jié)果表明，本文的計(jì)算理論和所編制的軟件能夠準(zhǔn)確地計(jì)算玻璃幕墻熱工性能，并且能夠推廣應(yīng)用到其他形式的幕墻體系中.

1 玻璃體系傳熱數(shù)值計(jì)算模型

傳熱系數(shù)是反應(yīng)玻璃熱工性能的重要指標(biāo).在幕墻設(shè)計(jì)過(guò)程中，采用不同型號(hào)的玻璃組成的玻璃體系傳熱系數(shù)往往未知.由于玻璃的傳熱過(guò)程可近似地認(rèn)為是一維熱傳導(dǎo)，因此本文基于一維熱傳導(dǎo)理論［12］，以雙層玻璃體系且封閉空腔填充空氣為例，運(yùn)用差分法建立玻璃體系傳熱系數(shù)的數(shù)值計(jì)算模型.

雙層玻璃由玻璃與封閉空腔組成，如圖1所示.當(dāng)傳熱過(guò)程處于穩(wěn)態(tài)時(shí)，玻璃體系的熱過(guò)程包括玻璃內(nèi)部的導(dǎo)熱、玻璃與環(huán)境、玻璃與玻璃之間的輻射換熱以及空腔內(nèi)的對(duì)流換熱.假設(shè)熱能由室內(nèi)傳向室外，根據(jù)能量守恒原理針對(duì)上述熱過(guò)程建立雙層玻璃傳熱數(shù)值計(jì)算模型如式（1）所示.

式中：［q］為熱流密度矩陣；［C］為熱系數(shù)矩陣，包含導(dǎo)熱系數(shù)、對(duì)流和輻射換熱系數(shù)；Δ［T］為溫差矩陣.將式（1）詳細(xì)展開(kāi)如式（2）所示.

式中：hc，in和hc，out分別為室內(nèi)和室外的對(duì)流換熱系數(shù)，W·（m2·K）-1；hr，in和hr，out分別為室內(nèi)和室外的輻射換熱系數(shù)，W·（m2·K）-1；hc，v和hr，v分別為封閉空腔的對(duì)流和兩玻璃界面之間的輻射換熱系數(shù)，W·（m2·K）-1；Tin，Tout分別為室內(nèi)外環(huán)境溫度，K；T11，T12，T21，T22分別表示雙層玻璃體系各界面的室內(nèi)、外表面的溫度，K；λg1和λg2為玻璃體系各玻璃板的導(dǎo)熱系數(shù)，W·（m2·K）-1；dg1和dg2為封閉空腔的厚度，m；q1～q5表示各熱過(guò)程中的熱流密度，W·m-2.

圖1 雙層玻璃熱過(guò)程計(jì)算模型Fig.1 Thermal process of double glazing

通常，雙層玻璃與室內(nèi)、外之間的對(duì)流換熱系數(shù)可以測(cè)得.封閉空腔內(nèi)的氣體由于相鄰兩玻璃界面的溫差驅(qū)使出現(xiàn)對(duì)流現(xiàn)象.封閉空腔內(nèi)的對(duì)流換熱系數(shù)hc，v與氣體努塞爾數(shù)Nu、氣體層的導(dǎo)熱系數(shù)λa和氣體層的厚度da有關(guān)，計(jì)算如式（3）所示.相關(guān)參數(shù)具體計(jì)算過(guò)程可參見(jiàn)文獻(xiàn)［13］.

輻射換熱系數(shù)是關(guān)于輻射面的輻射溫度值和發(fā)射率的函數(shù).由于玻璃面板表面溫度未知，所以式（2）中的輻射換熱系數(shù)為未知量.玻璃與室內(nèi)環(huán)境之間的輻射換熱系數(shù)計(jì)算公式如（4）所示.

式中：εin為室內(nèi)環(huán)境發(fā)射率；ε11為玻璃室內(nèi)一側(cè)界面的的表面發(fā)射率；σ為Boltzmann常數(shù).

根據(jù)能量守恒定律，當(dāng)雙層玻璃熱過(guò)程處于穩(wěn)態(tài)時(shí)，上述各熱過(guò)程的熱流密度q必相等.計(jì)算前，假定固體材料的熱屬性不隨溫度變化而變化，室內(nèi)、外的發(fā)射率、環(huán)境溫度和輻射溫度、對(duì)流換熱系數(shù)和玻璃的導(dǎo)熱系數(shù)、發(fā)射率均為已知值，其余量未知.計(jì)算時(shí)，由已知的室內(nèi)、外溫度值假設(shè)出玻璃體系各界面的溫度初值.此時(shí)可計(jì)算出玻璃體系與室內(nèi)、外的輻射換熱系數(shù)，也可計(jì)算出氣體相應(yīng)的熱特性參數(shù)，進(jìn)而計(jì)算出封閉空腔的對(duì)流換熱系數(shù)和兩片玻璃之間的輻射換熱系數(shù).求解式（2），可求得各熱過(guò)程中的熱流密度q.將所得的熱流密度的平均值帶入式（2）反算出玻璃體系各界面的溫度值.更新未知的參數(shù)，重復(fù)上述計(jì)算步驟，通過(guò)多次迭代計(jì)算直到收斂，便可求得各個(gè)表面的溫度值和熱流密度.玻璃體系的傳熱系數(shù)U由式（5）求解.至此，雙層玻璃的傳熱系數(shù)數(shù)值計(jì)算模型已經(jīng)建立.3層及3層以上的玻璃體系數(shù)值模型可以依此類推.

2 玻璃體系光學(xué)參數(shù)計(jì)算模型

2.1 單層玻璃光學(xué)性能計(jì)算模型

玻璃的光學(xué)性能同樣是反應(yīng)玻璃幕墻節(jié)能的一項(xiàng)重要指標(biāo).玻璃光學(xué)性能的指標(biāo)通常包括可見(jiàn)光透射比τv和反射比ρv；太陽(yáng)光直接透射比τsol與反射比ρsol；太陽(yáng)光直接吸收比Asol和總透射比g以及玻璃的遮陽(yáng)系數(shù)Cg等.玻璃體系一般由兩片及多片玻璃組合而成，其光學(xué)特性由單片玻璃的光學(xué)特性所決定，首先應(yīng)對(duì)單片玻璃的光學(xué)參數(shù)計(jì)算進(jìn)行研究.

單片玻璃在一定波長(zhǎng)范圍內(nèi)的透射比光譜數(shù)據(jù)τ（λ）與反射比光譜數(shù)據(jù)ρ（λ）可以通過(guò)試驗(yàn)測(cè)定.可見(jiàn)光波長(zhǎng)范圍通常在380～780nm之間，而太陽(yáng)光波長(zhǎng)在300～2500nm之間.對(duì)單片玻璃的光譜數(shù)據(jù)在可見(jiàn)光波長(zhǎng)范圍內(nèi)或太陽(yáng)光波長(zhǎng)范圍內(nèi)進(jìn)行積分，可求得其可見(jiàn)光透射比τv與反射比ρv、太陽(yáng)光直接透射比τsol與反射比ρsol和太陽(yáng)光直接吸收比

單片玻璃的太陽(yáng)光總透射比g是指太陽(yáng)輻射能量透過(guò)玻璃進(jìn)入室內(nèi)的熱量，既包括直接透過(guò)的部分，也包括玻璃吸收太陽(yáng)輻射后向室內(nèi)的二次傳熱.因此，單片玻璃的太陽(yáng)光透射比g與室內(nèi)外兩側(cè)的對(duì)流換熱系數(shù)有關(guān)，計(jì)算式如式（6）所示.遮陽(yáng)系數(shù)是指單片玻璃的太陽(yáng)光透射比g與通過(guò)3mm標(biāo)準(zhǔn)單層平板白玻璃的太陽(yáng)光透射比的比值.遮陽(yáng)系數(shù)Cg的計(jì)算如公式（7）所示.

式中：hin，hout分別為室內(nèi)室外的對(duì)流換熱系數(shù).

2.2 雙層玻璃光學(xué)性能計(jì)算模型

玻璃體系的可見(jiàn)光透射比光譜數(shù)據(jù)τ（λ）與反射比光譜數(shù)據(jù)ρ（λ）需由單片玻璃光譜數(shù)據(jù)來(lái)合成.本文再次以雙層玻璃為例（如圖1所示），假定室外只有太陽(yáng)輻射，室外和室內(nèi)環(huán)境反射比為0，根據(jù)能量守恒原理，采用差分法建立玻璃體系光譜數(shù)據(jù)合成計(jì)算數(shù)值模型，如式（8）所示.

式中：［S］是由雙層玻璃中各玻璃的光譜數(shù)據(jù)組成的光譜矩陣；［In］是各間層的光能矩陣；［Is］輻射能矩陣.

對(duì)［In］在一定波長(zhǎng)范圍內(nèi)進(jìn)行多次求解，得到各間層的光譜數(shù)據(jù).再按式（10）擬合多層玻璃的光譜數(shù)據(jù).將式（8）展開(kāi)如式（9）所示.

式中：“＋”和“-”分別表示輻射流向室外和流向室內(nèi)；τ1（λ）和τ2（λ）分別為玻璃1和2的透過(guò)比光譜數(shù)據(jù)；ρ11（λ），ρ12（λ），ρ21（λ）和ρ22（λ）分別為玻璃體系各界面的反射比光譜數(shù)據(jù)，分別表示各間層流向室外和室內(nèi)的光能數(shù)據(jù)；Is（λ）表示輻射光能值.

當(dāng)?shù)玫诫p層玻璃的光譜數(shù)據(jù)后，可按照單層玻璃的計(jì)算方法計(jì)算雙層玻璃的光學(xué)性能參數(shù).3層及3層以上的玻璃體系光學(xué)參數(shù)數(shù)值計(jì)算模型可以依此類推.

3 框體系熱工性能數(shù)值計(jì)算模型

3.1 穩(wěn)態(tài)熱傳導(dǎo)微分方程

框式玻璃幕墻是由玻璃體系與支承結(jié)構(gòu)體系構(gòu)成，包括明框與隱框兩種類型.框是框式玻璃幕墻的重要組成部分，因此，在其整體熱工性能中框的熱工性能具有重要的作用.框通常是由金屬型材和斷熱膠條等材料組成，為了能夠起到很好的保溫隔熱作用，其斷面設(shè)計(jì)構(gòu)造較為復(fù)雜，包括較多的空腔等，這直接導(dǎo)致框架體系內(nèi)部的傳熱過(guò)程和溫度場(chǎng)分布比玻璃體系的要復(fù)雜.框的傳熱過(guò)程可近似為沿著框斷面的一個(gè)二維熱傳導(dǎo)過(guò)程.為了能夠?qū)蝮w系的熱工性能和內(nèi)部的傳熱機(jī)理以及溫度場(chǎng)分布進(jìn)行分析，本文以二維穩(wěn)態(tài)熱傳導(dǎo)理論為基礎(chǔ)，利用有限單元法建立框體系的傳熱數(shù)值計(jì)算模型.

根據(jù)二維穩(wěn)態(tài)熱傳導(dǎo)理論建立的無(wú)內(nèi)熱源的熱傳導(dǎo)微分方程如式（11）所示［15］.

式中：kx和ky是材料兩個(gè)主方向（x，y方向）的導(dǎo)熱系數(shù)；T為溫度.

本文利用平面三節(jié)點(diǎn)三角形單元對(duì)體系進(jìn)行離散，設(shè)定每個(gè)節(jié)點(diǎn)僅有一個(gè)溫度自由度.二維穩(wěn)態(tài)熱傳導(dǎo)的有限元求解格式為

式中：［K］為三角形單元熱傳導(dǎo)矩陣；［T］為三角形單元溫度矩陣；［P］為溫度載荷矩陣.

推導(dǎo)出三角形單元二維熱傳導(dǎo)矩陣［K］計(jì)算式為

式中：A為三角形單元面積；b和c分別為三角形單元兩節(jié)點(diǎn)的縱、橫坐標(biāo)差值；i，j，k分別為節(jié)點(diǎn)編號(hào)［16］.

熱力學(xué)邊界條件包括對(duì)流、輻射等.單元熱力學(xué)邊界條件不同，熱傳導(dǎo)矩陣［K］和溫度載荷矩陣［P］都需要進(jìn)行相應(yīng)的修正.下文將僅以對(duì)流邊界為例進(jìn)行詳細(xì)推導(dǎo).

3.2 對(duì)流換熱邊界條件

當(dāng)單元邊界為對(duì)流換熱時(shí)，二維穩(wěn)態(tài)熱傳導(dǎo)微分方程為

式中：hf為對(duì)流換熱系數(shù)；Tf為外界流體溫度；nx，ny分別為法向向量.

本文以單元ij邊與外界進(jìn)行對(duì)流換熱為例，推導(dǎo)單元熱傳導(dǎo)矩陣和單元溫度載荷矩陣的修正式如下：

式中：lij表示三角形單元的各邊長(zhǎng).同理，當(dāng)三角形單元ik邊和jk邊與外界進(jìn)行對(duì)流換熱時(shí)，單元熱傳導(dǎo)矩陣和溫度載荷矩陣的修正式也可以推出.

當(dāng)對(duì)流換熱系數(shù)hf無(wú)窮大時(shí)，可認(rèn)為流體與物體邊界的溫度相等，即常溫邊界條件.在計(jì)算時(shí)，hf越大，計(jì)算結(jié)果越準(zhǔn)確.

4 玻璃幕墻熱工計(jì)算軟件研發(fā)

4.1 軟件研發(fā)

以Visual C＋＋為開(kāi)發(fā)工具，在前文所建立的理論推導(dǎo)和數(shù)值計(jì)算模型的基礎(chǔ)上，基于ObjectARX.技術(shù)對(duì)AutoCAD進(jìn)行了二次開(kāi)發(fā)，實(shí)現(xiàn)了玻璃幕墻熱工性能分析的軟件化.相比LBNL系列熱工軟件，本文所開(kāi)發(fā)軟件各模塊統(tǒng)一度更高，包括玻璃體系的光學(xué)、傳熱和結(jié)露性能分析、框體系的二維有限元計(jì)算以及幕墻整體熱工性能計(jì)算等，并且建立了中國(guó)玻璃數(shù)據(jù)庫(kù)，包括了中國(guó)所使用的大多數(shù)玻璃種類，更能適合我國(guó)實(shí)際情況.由于基于AutoCAD軟件平臺(tái)，因此，相對(duì)于其他國(guó)外軟件具有更好的建模功能，并實(shí)現(xiàn)了對(duì)材料進(jìn)行直接填充等功能，在操作方面更加便捷.利用CGAL（Computational Geometry Algorithm Library）提供的三角化程序，對(duì)框體系進(jìn)行了有限元網(wǎng)格的自適應(yīng)劃分，可對(duì)框體系在對(duì)流、輻射以及多種邊界耦合作用下的的熱工性能和溫度場(chǎng)進(jìn)行有限元數(shù)值模擬分析.后處理可實(shí)現(xiàn)繪制模型等溫線、溫度云圖和熱流矢量圖等功能.所編制的軟件不僅可以應(yīng)用于框式玻璃幕墻熱工性能計(jì)算中，同時(shí)也可應(yīng)用于全玻式幕墻和非透明幕墻等多種幕墻體系中.軟件流程圖如圖2所示.

圖2 玻璃幕墻軟件計(jì)算程序流程Fig.2 The flow chart of software calculation program

4.2 軟件對(duì)比驗(yàn)證

4.2.1 玻璃體系熱工性能算例驗(yàn)證及分析

以間層填充為空氣、氬氣、氪氣和氙氣的3層玻璃體系為例，將采用本文方法所計(jì)算的熱工性能結(jié)果與美國(guó)伯克利大學(xué)的LBNL系列熱工軟件計(jì)算結(jié)果進(jìn)行對(duì)比.計(jì)算時(shí)，采用相同的玻璃型號(hào)，確保材料的一致性.玻璃厚度為3mm，氣體間層厚度為12 mm.設(shè)室內(nèi)、外環(huán)境溫度分別為20℃和-20℃，對(duì)流換熱系數(shù)分別為1.8和16.0W·（m2·K）-1，假設(shè)室內(nèi)、室外平均輻射溫度等于室內(nèi)、室外環(huán)境溫度.計(jì)算結(jié)果如表1，2所示.

表1 3層玻璃傳熱系數(shù)計(jì)算結(jié)果對(duì)比Tab.1 The contrast of 3-layer glass system thermal transmittance

表2 3層玻璃光學(xué)參數(shù)計(jì)算結(jié)果對(duì)比（內(nèi)填空氣）Tab.2 The contrast of 3-layer glass system optic spara meters

由表1和表2的3層玻璃傳熱系數(shù)和光學(xué)參數(shù)計(jì)算結(jié)果對(duì)比表明，本文與LBNL軟件計(jì)算結(jié)果十分接近，從而驗(yàn)證了本文的正確性.同時(shí)，玻璃體系的傳熱系數(shù)隨著惰性氣體惰性的增強(qiáng)而減小，也就意味著保溫隔熱性能越來(lái)越好.

4.2.2 框體系熱工性能算例驗(yàn)證

對(duì)一塊導(dǎo)熱系數(shù)為0.18W·（m·K）-1的正方形板分別采用通用有限元軟件ANSYS，LBLN系列軟件熱工模塊THERM和本文計(jì)算方法進(jìn)行了二維熱傳導(dǎo)的計(jì)算.正方形板尺寸為400mm×400 mm；上邊界對(duì)流換熱系數(shù)為6W·（m2·K）-1，流體溫度為20℃；下邊界對(duì)流換熱系數(shù)為8W·（m2·K），流體溫度為10℃；左邊界溫度為15℃；右邊界溫度為30℃.正方形板尺寸及節(jié)點(diǎn)編號(hào)如圖3所示.

對(duì)正方形板在設(shè)置好的邊界條件下的溫度場(chǎng)進(jìn)行計(jì)算，3種軟件計(jì)算結(jié)果如表3所示.由對(duì)比結(jié)果可知，本文各節(jié)點(diǎn)溫度值計(jì)算結(jié)果與THERM計(jì)算結(jié)果十分接近，與ANSYS計(jì)算結(jié)果誤差也在1%左右，從而驗(yàn)證了本文框體系數(shù)值計(jì)算方法的正確性.

圖3 正方形板尺寸及節(jié)點(diǎn)編號(hào)Fig.3 The dimensions and points number of square plate

表3 節(jié)點(diǎn)溫度值計(jì)算結(jié)果Tab.3 The point temperature calculation results

由于框式幕墻的框體系內(nèi)部構(gòu)造較復(fù)雜，為了進(jìn)一步驗(yàn)證本文框體系熱工性能計(jì)算理論的正確性，對(duì)在夏熱冬冷地區(qū)實(shí)際工程中所采用的幕墻框在夏、冬季以及過(guò)度季節(jié)環(huán)境邊界條件下的熱工性能和溫度場(chǎng)分布進(jìn)行模擬計(jì)算，并與THERM軟件計(jì)算結(jié)果進(jìn)行對(duì)比.圖4為某玻璃幕墻框截面尺寸圖.模型中材料包括鋁合金、聚氨酯密封膠、硅酮結(jié)構(gòu)膠和浮法玻璃.玻璃體系為雙層玻璃，玻璃的厚度為6mm，封閉空腔中的氣體為空氣，厚度為12mm.材料的導(dǎo)熱系數(shù)和表面發(fā)射率如表4所示.

表4 材料熱工性能參數(shù)Tab.4 Thermal property of materials

將各種材料填充到框截面模型中，如圖所5示.對(duì)截面進(jìn)行網(wǎng)格劃分如圖6所示.邊界條件類型包括室內(nèi)、室外和絕緣邊界，其中上、下表面分別定義為室內(nèi)和室外邊界，兩側(cè)為絕緣邊界，不同季節(jié)的邊界條件如表5所示.分別計(jì)算各季節(jié)下截面最大、最小溫度值以及截面的傳熱系數(shù)，并與THERM軟件結(jié)果進(jìn)行對(duì)比，對(duì)比結(jié)果如表6所示.從表6中的數(shù)據(jù)可知，本文計(jì)算結(jié)果與THERM計(jì)算結(jié)果相差在2%以內(nèi)，進(jìn)一步表明本文計(jì)算理論的合理性.在不同季節(jié)下的框截面溫度場(chǎng)云圖如圖7所示，可知不同季節(jié)下框截面的溫度場(chǎng)分布大不相同.在夏季和過(guò)渡季節(jié)，框表面的溫度均高于室內(nèi)，在冬季則低于室內(nèi).

圖4 框截面尺寸（單位：mm）Fig.4 The dimensions of frame section（unit：mm）

圖5 模型材料填充Fig.5 The materials filling model

圖6 模型網(wǎng)格劃分Fig.6 The model element mesh

表5 不同季節(jié)的邊界條件Tab.5 Boundary conditions of different seasons

表6 對(duì)比驗(yàn)證結(jié)果Tab.6 Comparison of the calculated results with THERM results

圖7 不同季節(jié)邊界條件溫度場(chǎng)云圖Fig.7 The temperature nephogram in different seasons

5 結(jié)論

分別針對(duì)框式玻璃幕墻的玻璃體系的傳熱、光學(xué)性能和框體系的熱工性能和溫度場(chǎng)計(jì)算理論進(jìn)行了研究，建立了各種參數(shù)的數(shù)值計(jì)算模型，編制了計(jì)算程序，并與多種軟件計(jì)算結(jié)果進(jìn)行對(duì)比.通過(guò)上述研究，得到如下結(jié)論：

（1）建立了玻璃體系傳熱和光學(xué)參數(shù)數(shù)值計(jì)算模型，編制了玻璃體系熱工性能計(jì)算程序.

（2）對(duì)框體系在不同環(huán)境邊界條件下的傳熱計(jì)算進(jìn)行了理論推導(dǎo)，驗(yàn)證了計(jì)算方法的有效性，并對(duì)某實(shí)際工程中的框體系在不同季節(jié)下的熱工性能和溫度場(chǎng)分布進(jìn)行了數(shù)值模擬.

（3）以 Visual C＋＋為開(kāi)發(fā)工具，基于ObjectARX.技術(shù)對(duì)AutoCAD進(jìn)行了二次開(kāi)發(fā)，實(shí)現(xiàn)了框式玻璃幕墻熱工性能分析的軟件化，相比國(guó)外其他幕墻熱工軟件具有建模功能強(qiáng)、各模塊統(tǒng)一度高、更適應(yīng)我國(guó)規(guī)范等優(yōu)點(diǎn)，同時(shí)，軟件不僅可以應(yīng)用于框式玻璃幕墻熱工性能計(jì)算中，也可應(yīng)用于全玻式幕墻和非透明幕墻等多種幕墻體系中.

［1］張其林.玻璃幕墻結(jié)構(gòu)［M］.濟(jì)南：山東科技出版社，2006.ZHANG Qilin.Glass curtain wall structures［M］.Ji’nan：Shandong Science and Technology Press，2006.

［2］清華大學(xué)建筑節(jié)能研究中心.中國(guó)建筑節(jié)能年度研究發(fā)展報(bào)告2011［M］.北京：中國(guó)建筑工業(yè)出版社，2011.Tsinghua University Building Energy Research Center.2011 annual report on China building energy efficiency ［M］.Beijing：China Architecture ＆Building Press，2011.

［3］清華大學(xué)建筑節(jié)能研究中心.中國(guó)建筑節(jié)能年度研究發(fā)展報(bào)告2010［M］.北京：中國(guó)建筑工業(yè)出版社，2010.Tsinghua University Building Energy Research Center.2010 annual report on China building energy efficiency ［M］.Beijing：China Architecture and Building Press，2010.

［4］Adelqui Fissore，Nestor Fonseca.Experimental study of the thermal balance of a window，design description［J］.Building and Environment，2007，42（9）：3309.

［5］Shih N J，Huang Y S.An analysis and simulation of curtain wall reflection glare［J］.Building and Environment，2001，36（1）：619.

［6］International Organization for Standardization.ISO 15099：2003 Thermal performance of windows，doors and shadingdevices-Detailed calculations ［S］. London： International Organization for Standardization Press，2003.

［7］National Fenestration Rating Council. NFRC 100：2004 Procedure for determining fenestration product U-factors［S］.Silver Spring：National Fenestration Rating Council Press，2004.

［8］YU Kezhi，DING Guoliang，CHEN Tianji.Simulation of air curtains for vertical display cases with a two-fluid model［J］.Applied Thermal Engineering，2007，27（9）：2583.

［9］Chen Y G.Parametric evaluation of refrigerated air curtains for thermal insulation ［J］.International Journal of Thermal Sciences，2009，48（7）：2988.

［10］劉濤.基于有限單元法的玻璃幕墻熱工性能的分析和設(shè)計(jì)［D］.上海：同濟(jì)大學(xué)土木工程學(xué)院，2010.LIU Tao.Analysis and design of thermal performance of curtain walls based on FE ［D］.Shanghai：College of Civil Engineering of Tongji University，2010.

［11］中華人民共和國(guó)住房和城鄉(xiāng)建設(shè)部.JGJ／T151—2008建筑門(mén)窗玻璃幕墻熱工計(jì)算規(guī)程［S］.北京：中國(guó)建筑工業(yè)出版社，2008.Ministry of Housing and Urban-Rural Development of the People’s Republic of China.JGJ／T 151-2008 Calculation specification for thermal performance of windows，doors，and glass curtain-walls［S］.Beijing：China Architecture ＆Building Press，2008.

［12］楊世銘，陶文銓.傳熱學(xué)［M］.北京：高等教育出版社，2010.YANG Shiming，TAO Wenquan.Heat transfer［M］.Beijing：Higher Education Press，2010.

［13］Wdght L J.A correlation to quantify convective heat transfer between vertical window glazings［J］.ASHRAE Transactions，1996，102（1）：940.

［14］中華人民共和國(guó)國(guó)家質(zhì)量技術(shù)監(jiān)督局.GB／T2680—94建筑玻璃可見(jiàn)光透射比、太陽(yáng)光直接透射比、太陽(yáng)能總透射比、紫外線透射比及有關(guān)窗玻璃參數(shù)的測(cè)定［S］.北京：中國(guó)標(biāo)準(zhǔn)出版社，1994.The State Bureau of Quality and Technical Supervision of the People’s Republic of China.GB／T2680-94 Determination of light transmittance，solar direct transmittance，total solar energy transmittance and ultraviolet transmittance for glass in building and related glazing factors［S］.Beijing：Standards Press of China，1994.

［15］王勖成.有限單元法［M］.北京：清華大學(xué)出版社，2010.WANG Xucheng.Finite element method ［M］.Beijing：Tsinghua University Press，2010.

［16］孔祥謙.有限單元法在傳熱學(xué)中的應(yīng)用［M］.北京：科學(xué)出版社，1998.KONG Xiangqian.The application of the finite element method in heat transfer［M］.Beijing：Science Press，1998.