田 苗, 李 俊

(1. 東華大學(xué) 服裝·藝術(shù)設(shè)計學(xué)院， 上海 200051； 2. 東華大學(xué) 功能防護服裝研究中心， 上海 200051；3. 東華大學(xué) 現(xiàn)代服裝設(shè)計與技術(shù)教育部重點實驗室， 上海 200051)

?

數(shù)值模擬在熱防護服裝性能測評中的應(yīng)用

田 苗1,2,3, 李 俊1,2,3

(1. 東華大學(xué) 服裝·藝術(shù)設(shè)計學(xué)院， 上海 200051； 2. 東華大學(xué) 功能防護服裝研究中心， 上海 200051；3. 東華大學(xué) 現(xiàn)代服裝設(shè)計與技術(shù)教育部重點實驗室， 上海 200051)

對熱防護服裝性能合理的評價在公共安全領(lǐng)域具有重要意義。數(shù)值模擬成為繼試驗研究之后相關(guān)領(lǐng)域的重要研究手段。從模型的建立及驗證2方面出發(fā)，主要對傳熱模型、皮膚燒傷預(yù)測模型和火場環(huán)境仿真等相關(guān)研究進行回顧，并對典型模型的發(fā)展過程、特點和不足進行歸納。根據(jù)相關(guān)技術(shù)的發(fā)展，預(yù)測了熱防護服裝數(shù)值模擬的發(fā)展趨勢。未來將在完善熱濕傳遞模型的基礎(chǔ)上，發(fā)展全尺度的數(shù)值模擬，并結(jié)合人體熱調(diào)節(jié)等模型對防護服進行多方位的評價。通過在優(yōu)化熱防護服裝方面的應(yīng)用，充分體現(xiàn)數(shù)值模擬研究的實用價值。

數(shù)值模擬； 熱防護服； 性能測評； 皮膚燒傷； 傳熱模型

隨著計算機技術(shù)的發(fā)展，數(shù)值計算方法的研究得到了極大地推動，從而使數(shù)值模擬的規(guī)模和精度不斷優(yōu)化。數(shù)值方法逐漸成為與理論分析、實驗研究相并列的三大科學(xué)研究手段之一，可以對許多原本無法或很難用理論分析求解的復(fù)雜問題進行模擬求解，比物理實驗具有更大的自由度和靈活度；同時通過實現(xiàn)虛擬，現(xiàn)實環(huán)境、可視分析環(huán)境和協(xié)同仿真環(huán)境，數(shù)值模擬可以在一定程度上替代物理實驗，節(jié)約試驗成本，加快研究進度，擴大研究范圍[1]。

在目前服裝熱防護性能的測評中，普遍使用物理試驗的方法，即通過TPP(Thermal Protective Performance，熱防護性能)試驗或燃燒假人試驗評估服裝的熱防護性能。這種試驗具有破壞性及不可重復(fù)性，試驗結(jié)果的準確度會受到環(huán)境和人為因素的影響，并且成本較高。由于服裝本身屬性特殊，服裝領(lǐng)域相關(guān)的數(shù)值模擬研究起步較晚，且多為國外學(xué)者所關(guān)注，國內(nèi)相關(guān)研究鮮見。隨著服裝領(lǐng)域?qū)嶒炑芯矿w系的逐漸完善，研究方法亟待突破，而數(shù)值模擬正是關(guān)鍵所在。服裝領(lǐng)域?qū)⒊蔀槔^暖通空調(diào)、機械制造、材料科學(xué)和航空學(xué)等領(lǐng)域[2-3]之后，數(shù)值模擬應(yīng)用的重要舞臺。

熱防護服裝數(shù)值模擬研究的目的是預(yù)測隨時間變化的織物、皮膚、空氣層中溫度、熱流的分布，燒傷等級的量化，及達到皮膚燒傷所需的時間等[4]。通過探索防護服中熱傳遞的物理現(xiàn)象，提高服裝的性能，并應(yīng)用在消防安全領(lǐng)域。數(shù)值模擬研究的關(guān)鍵問題是模型的建立和驗證。

1 數(shù)值模型的建立

進行數(shù)值模擬，首先要建立反映問題各變量之間的微分方程及相應(yīng)的定解條件。處在火災(zāi)環(huán)境時，人體傷害的判定基于皮膚燒傷，造成燒傷的直接原因是火場環(huán)境，而在火場和人體皮膚之間的唯一屏障是消防服系統(tǒng)，因此，將數(shù)值模型的建立劃分為織物-空氣層-皮膚熱傳遞模型，皮膚模型和皮膚燒傷預(yù)測模型以及火場環(huán)境仿真3個部分進行回顧和評價。

1.1 織物-空氣層-皮膚熱傳遞模型

數(shù)值研究的環(huán)境條件主要分為2種：強熱流閃火條件和中等熱流熱輻射條件。閃火主要發(fā)生在石油化工行業(yè)，時間少于5 s但伴隨強熱流[5]。由于許多消防員的燒傷是長時間暴露于中等輻射熱源引起的[6]，因此，有學(xué)者開展了熱輻射條件下的相關(guān)研究。另外，數(shù)學(xué)模型只能在一定程度上滿足實際條件，為了簡化模型，一般需要進行假設(shè)。

1.1.1 模型假設(shè)

在模型維數(shù)和儀器形態(tài)的假設(shè)方面，Torvi[5]，Mell[6]，Ghazy[7]，Song[8]，Chitrphiromsri[9]等學(xué)者均采用了基于平壁假設(shè)的一維傳熱模型。Crown等[10]實驗研究證明臺式測試不能獲得織物熱收縮對其熱防護性能的影響，因此，朱方龍等[11]假設(shè)模擬人體的圓柱體儀器內(nèi)是一維徑向傳熱。Sawcyn等[12]則運用近似二維的方法改進防護面料臺式測試中空氣層的熱傳遞模型。

對傳熱模型進行的一般假設(shè)為：1) 織物為輻射灰體，質(zhì)量輸運傳遞可忽略不計；2) 輻射熱流量能夠滲透到面料內(nèi)部一定深度；3) 織物面料的熱屬性與溫度相關(guān)；4) 人體體核溫度保持37 ℃恒溫等。另外，不同學(xué)者根據(jù)各自的研究領(lǐng)域也進行了其他方面的假設(shè)，如：Torvi[5]，Song[8]，Chitrphiromstri等[9]通常不考慮空氣層中熱傳遞形式之間的相互作用；Ghazy等[7]則假設(shè)空氣層中傳導(dǎo)和輻射產(chǎn)生耦合作用等。

1.1.2 微分方程的建立

在建立數(shù)值模型時，根據(jù)特定的邊界條件，一般將模型分為織物、空氣層和皮膚3個部分。在織物暴露于熱源的表面發(fā)生輻射和對流熱傳遞?？諝鈱又械膫鳠岱绞奖容^復(fù)雜，當其厚度小于6.4 mm時，傳熱方式為傳導(dǎo)和輻射，厚度為6.4 mm或更大時，空氣發(fā)生自然對流[5]，因此，建立空氣層中的微分方程時要考慮厚度的影響。

Torvi[5]利用有限元法建立了織物-空氣層-皮膚系統(tǒng)的微分方程

(1)

式(1)的初始和邊界條件為

t=0時，

(2)

x=0，t>0時，

(3)

x=Lb，t>0時，

(4)

式中：CA為顯熱容，J/K；qrad為輻射熱流，W/m2；k為導(dǎo)熱率，W/(m·℃)；γ為消光系數(shù)，1/m；T為溫度，K；t為時間，s；Ti(x)為織物-空氣層-皮膚系統(tǒng)的初始溫度分布；Ta為環(huán)境溫度，K；h為對流傳熱系數(shù)W/(m2·K)。

Sawcyn等[12]開發(fā)二維改進模型，模擬熱防護織物和傳感器之間水平空氣層的熱傳遞。將底部邊界作為一系列等溫矩形計算輻射傳熱，每個單元(j)到傳感器的輻射傳熱通過下式計算(忽略封閉空間與輻射發(fā)射率相關(guān)的散射和吸收)：

(5)

式中：A和ε分別為表面積(m2)和發(fā)射率。

考慮到空氣層中傳導(dǎo)和輻射之間的相互作用，Ghazy等[4]得出當空氣層厚度小于6.4 mm時一維瞬態(tài)傳導(dǎo)-輻射傳熱方程。

(6)

Song[8]、Chitrphiromstri等[9]使用了Torvi[5]的方法建立了微分方程和邊界條件。Mell等[6]認為各層織物之間輻射是雙向的，但未將織物的光學(xué)性能(透射率、反射率)用作模型的輸入?yún)?shù)。朱主龍等[11]建立模型時則未考慮源項，因為溫度未達到發(fā)生反應(yīng)的要求。

1.1.3 微分方程的求解

Torvi[5]使用有限元法解決微分方程，有限元計算機程序在Microsoft?QuickBASICTM4.5[13]中編寫。Song等[8]使用有限差分法解決假人表面每個傳感器處的微分方程。對于非線性的輻射方程，采用Gauss-Seidel逐點迭代的方法解決。Ghazy等[14]則運用有限體積法對織物、空氣層和皮膚的能量方程、邊界條件進行離散，使用完全隱式法進行時間離散。Sawcyn[15]通過內(nèi)部代碼FASTEST3D (Flow Analysis by Solving Transport Equations Simulating Turbulence)解決控制方程。代碼基于有限體積法并寫入曲線坐標，SIMPLE算法解決動量方程，中心差分法離散發(fā)散熱流，代數(shù)方程則通過Stone的半隱式程序求解。

可以發(fā)現(xiàn)，數(shù)值模型的建立和計算方法逐漸從有限差分法和有限元法，發(fā)展到有限體積法[16]。

1.2 皮膚模型和皮膚燒傷預(yù)測模型

1.2.1 皮膚模型

皮膚指身體表面包在肌肉外面的組織，主要承擔(dān)保護身體、排汗、感覺冷熱和壓力的功能。人體皮膚由表皮、真皮、皮下組織組成。常見的皮膚模型主要有Pennes皮膚模型[17]和熱波皮膚模型[18]。

Pennes模型是目前最常用的皮膚傳熱模型，其傳熱方程形式為

(7)

式中：ρs為皮膚的密度，kg/m3；cps為皮膚的比熱容，J/(kg·℃)；λs為皮膚的導(dǎo)熱系數(shù)，W/(m·℃)；ρb為血液的密度，kg/m3；cpb為血液的比熱容，J/(kg·℃)；λb為血液的導(dǎo)熱系數(shù)，W/(m·℃)；ωb為血流灌注率。方程的內(nèi)邊界條件設(shè)定為體核37 ℃恒溫，初始條件為皮膚外、內(nèi)表面初始溫度值之間呈(34～37) ℃線性分布。

Pennes模型的建立基于經(jīng)典的傅里葉定律，主要針對皮膚瞬時熱流暴露，但器官系統(tǒng)需要熱弛豫時間τ來積聚傳遞到皮膚內(nèi)部所需的熱量。Liu等[18]提出了與皮膚模擬傳感器表面熱流相結(jié)合的熱波皮膚模型。

(8)

熱波皮膚模型代表的是有限速度熱傳遞，與傳統(tǒng)的熱傳導(dǎo)方程相比預(yù)測更加真實，但在描述快速傳熱的合理性方面存在一定爭議。

1.2.2 皮膚燒傷預(yù)測模型

一般通過估算穿過面料造成皮膚燒傷的總傳熱量來評價服裝的熱防護性能[14]。臺式測試利用Stoll等[19]提出的燒傷準則評價面料樣本的TPP，而假人測試則是利用Henriques燒傷積分[20]評價整件服裝不同位置的TPP。

Stoll二級燒傷準則是出現(xiàn)較早的皮膚燒傷預(yù)測模型。在恒定熱流的暴露下，銅片熱流計溫度歷史曲線和Stoll曲線相交，相交點的橫坐標即為達到二級燒傷所需的時間。這種預(yù)測二級燒傷的方法無需復(fù)雜的數(shù)學(xué)運算，即可迅速預(yù)測燒傷。

Henriques燒傷積分模型是目前應(yīng)用最為廣泛的皮膚燒傷模型。Henriques發(fā)現(xiàn)皮膚的損傷可以通過化學(xué)反應(yīng)過程來表示，提出一階阿倫尼烏斯(Arrhenius)方程，用于表示組織損傷速率。

(9)

式中：Ω為皮膚燒傷程度，無量綱；P為頻率因子，s-1；E為皮膚的活化能，J/mol；R為摩爾常數(shù)，8.31 J/(mol·℃)；t為皮膚溫度高于44 ℃的時間，s；T為皮膚下深度為x的絕對溫度，K。對式(9)進行積分，得到Ω。皮膚溫度T>44 ℃且Ω=0.53 時，皮膚一級燒傷；皮膚溫度T>44 ℃且Ω≥1 時，皮膚二級燒傷。該模型適用于長時間低熱流條件下皮膚表層燒傷評價，而對于短時間高熱流條件下是否適用尚無定論。

另外，Autrique等[21]通過對Henriques、Mehta 和 Wong、Weaver和Stoll等皮膚燒傷預(yù)測模型的靈敏度分析，發(fā)現(xiàn)多數(shù)模型參數(shù)都未進行很好的定義，并認為影響模型準確度的關(guān)鍵參數(shù)為表皮和真皮層的體積熱容、消光系數(shù)和皮膚厚度。對于燒傷預(yù)測，活化能是保證預(yù)測有效性的關(guān)鍵參數(shù)。

1.3 火場環(huán)境仿真

消防員作業(yè)時需頻繁出入火災(zāi)現(xiàn)場，模擬火場環(huán)境可以更加真實可靠地評價消防服的熱防護性能。消防科學(xué)與火災(zāi)科學(xué)在該方向的研究成果較為顯著。

火災(zāi)試驗是一種破壞性試驗，目前能夠進行火場環(huán)境仿真的為TPP測試及燃燒假人試驗，它們分別從織物和服裝的層次完成火場環(huán)境中的測試。

作為模型實驗的一種補充研究方法，計算機數(shù)值模擬得到了越來越多的應(yīng)用。得益于計算機硬件技術(shù)的快速發(fā)展，計算流體動力學(xué)(CFD)成為火災(zāi)數(shù)值模擬的主要手段[22]。耿鵬[23]針對室內(nèi)火災(zāi)的過程利用FLUENT、CFAST軟件，以典型的建筑室內(nèi)結(jié)構(gòu)及火災(zāi)荷載為例，研究室內(nèi)火災(zāi)的煙氣熱量分布、建筑火災(zāi)中煙氣的溫度和煙氣流動速度等。根據(jù)計算得到火災(zāi)和煙氣在建筑物內(nèi)的發(fā)展情況，確定建筑物的火災(zāi)安全狀況，并分析有關(guān)消防設(shè)施的功能、人員的安全疏散。火災(zāi)發(fā)生時會形成密集湍流，并伴有燃燒和輻射熱。在服裝領(lǐng)域，Jiang等[24]利用CFD軟件FrontFlow/Red(FFR)模擬燃燒室中火焰燃燒時的氣流和熱傳遞。FFR綜合運用了Large Eddy Simulation(LES)和Reynolds Averaged Navier-Stokes(RANS)模擬湍流，可以通過超級計算機進行高效的并行計算。

1.4 模型的發(fā)展過程、特點與不足

選取具有代表性的歷年模型，對其發(fā)展過程以及模型的特點與不足進行了歸納，如表1所示。

表1 模型發(fā)展過程

熱防護服裝數(shù)值模擬研究主要分為3個階段：第1階段單純考慮織物-空氣層-皮膚系統(tǒng)的熱傳遞建立傳熱模型。第2階段出現(xiàn)了熱濕傳遞模型(在普通織物常溫狀態(tài)下的研究比較常見)。由于濕傳遞機制復(fù)雜，并且熱傳遞模型仍然存在問題，因此，有學(xué)者仍致力于完善傳熱模型。對于傳熱模型的優(yōu)化集中在對空氣層中傳熱的討論：一種認為空氣層中的傳熱方式會產(chǎn)生相互作用，進行傳導(dǎo)-輻射耦合模擬；另一種假設(shè)傳熱方式相互獨立，但利用二維的方法進行模擬。第3階段，開始利用CFD進行服裝整體的熱防護性能研究。

2 數(shù)值模型的驗證

由于數(shù)值模型求解的問題比較復(fù)雜，如非線性方程的數(shù)值求解方法在理論上不夠完善，需要通過試驗來加以驗證。在服裝領(lǐng)域的數(shù)值模擬主要采用數(shù)學(xué)方法和試驗方法來進行模型驗證。

2.1 數(shù)學(xué)方法驗證

Torvi[5]利用幾個簡化的案例，如：半無限大固體[25]，兩層半無限大固體，集總熱容分析和擁有變導(dǎo)熱系數(shù)的半無限大固體等，分別驗證計算機程序的各個部分。

Mell等[6]比較了數(shù)值模擬結(jié)果與簡化問題的精確解。通過比較檢測數(shù)值模擬的效果，確定一定精確程度的網(wǎng)格解析度。確認數(shù)值模型是否可以在材料界面很好地解決傳熱系數(shù)不連續(xù)的問題。與傳導(dǎo)熱流不同，輻射熱流直接模擬獲得，它的精確度取決于輻射熱傳遞物理模型的準確性以及光學(xué)性能參數(shù)的應(yīng)用，而不是微分計算的準確性。

早期研究數(shù)值模擬的學(xué)者，如Torvi[5]和Mell等[6]均使用數(shù)學(xué)方法計算精確解或比較研究的方式對模型進行了驗證。Song[8]、Chitrphiromsri[9]和Sawcyn[15]等的數(shù)值模擬方法多是從Torvi[5]的成果演化而來，因此，可以認為這些學(xué)者只是在文獻中省略了數(shù)學(xué)方法驗證的過程。

2.2 實驗方法驗證

實驗驗證是目前數(shù)值模擬當中常用的驗證方法，根據(jù)標準慣例設(shè)計實驗對結(jié)果進行驗證，也有少數(shù)學(xué)者，如Torvi[5]和Sawcyn[15]等將流場可視化實驗應(yīng)用到數(shù)值模擬的驗證當中。

2.2.1 臺式測試或假人測試

在TPP測試中，Torvi[5]利用紅外溫度計和熱電偶測量織物和襯墊的溫度。根據(jù)實驗結(jié)果，定性地討論了織物的熱響應(yīng)。朱方龍等[11]在同等受熱條件下分別模擬和測試了芳綸Metamax?織物和阻燃棉織物的表面溫度。Ghazy等[7]使用皮膚模擬材料及固定的空氣層厚度6.4 mm，模擬TPP測試[9]，比較測試、預(yù)測的皮膚表面溫度隨時間的變化趨勢，判斷模擬結(jié)果的準確性。

Song等[8]開發(fā)的數(shù)值模型可以預(yù)測強熱環(huán)境中通過服裝的熱傳遞。利用PyroMan燃燒假人系統(tǒng)，對比燒傷等級的比例和分布來驗證模型。Jiang等[24]開發(fā)了AFFECTION程序，計算通過服裝和人體皮膚的輻射、傳導(dǎo)綜合熱傳遞。通過模擬ISO6942對AFFECTION進行驗證，通過模擬ISO/FDIS 13506，進行假人層次燒傷實驗，驗證模擬結(jié)果。

2.2.2 流場可視化實驗

最早使用流場可視化方法的是Torvi[5]，他證實當空氣層厚度為6.4 mm或更大時，空氣發(fā)生自然對流，并利用流場可視化的方法描述了空氣層中的氣流形態(tài)。Sawcyn[15]進行流場可視化試驗時對儀器進行了適當?shù)母脑?，輔助確定了模型中產(chǎn)生對流傳熱的部分。

數(shù)學(xué)驗證的方法比較準確，但會受到模型假設(shè)的影響，且輻射傳熱部分不易驗證；試驗驗證結(jié)果比較直觀，但會受到試驗環(huán)境、人為因素等的影響，并且在試驗結(jié)果的評價上沒有統(tǒng)一的標準。因此，利用數(shù)學(xué)驗證和試驗驗證相結(jié)合的方法，可以相對完整地對模型進行驗證。但是，無論是TPP臺式測試還是燃燒假人測試，測試結(jié)果與模擬結(jié)果之間吻合度的判定都需要建立相對規(guī)范的標準。

3 熱防護服裝數(shù)值模擬的發(fā)展趨勢

在過去的50年中，計算機技術(shù)的發(fā)展使其性能得到了空前的提高，推動了數(shù)值計算方法的研究，使數(shù)值模擬的規(guī)模和精度不斷優(yōu)化。這為開展服裝熱防護領(lǐng)域的數(shù)值模擬研究創(chuàng)造了良好的條件。軟件和平臺的標準化，令通用設(shè)備上所實現(xiàn)的計算模型更容易擴展到大規(guī)模并行環(huán)境[26]，云計算[27]的出現(xiàn)提供了更為靈活的數(shù)值計算平臺。另外，人體熱生理模型[28-29]和室內(nèi)火災(zāi)[30-31]等研究為熱防護服性能的數(shù)值模擬提供了理論基礎(chǔ)，而3D掃描儀及CFD的應(yīng)用[32]則為相關(guān)研究提供了技術(shù)支持。熱防護服裝數(shù)值模擬的發(fā)展趨勢主要有：

1) 服裝-空氣層-皮膚系統(tǒng)熱濕傳遞模型準確度的提高。數(shù)值模擬中參數(shù)的敏感度研究說明，已有研究使用的某些輸入?yún)?shù)并不精確，因此，在建立新的模型時，要對經(jīng)驗參數(shù)進行論證，判斷其是否有效。采用準確的邊界條件，選取合適燃燒模型，考慮服裝熱收縮對熱傳遞的影響等都是提高模型準確度的途徑。

2) 服裝層次全尺度模擬規(guī)模的提高。目前服裝領(lǐng)域的數(shù)值模擬多局限于織物層次，僅有少數(shù)服裝層次數(shù)值模擬的研究。隨著計算機技術(shù)的發(fā)展，CFD技術(shù)的應(yīng)用及學(xué)科交叉的深入，服裝層次全尺度的數(shù)值模擬的規(guī)模和精度將得到提高。

3) 人體熱調(diào)節(jié)模型的應(yīng)用。目前對于人體熱調(diào)節(jié)模型以及熱舒適性的評價已經(jīng)發(fā)展到一定階段，隨著在高溫環(huán)境預(yù)測準確度的提高，有必要將其引入到熱防護服的舒適性評價中。從人體熱感覺的角度，完善熱防護服在火災(zāi)環(huán)境中的數(shù)值模擬，使其不僅可以進行熱防護性能的評價，而且可以進行熱舒適性能的評價。

4) 室內(nèi)火災(zāi)與煙氣流動模擬。火災(zāi)環(huán)境中，消防服或其他建筑材料在受熱或燃燒時會產(chǎn)生一定的煙氣，煙氣的危害性主要是其毒性和遮光性。對于燃燒假人的室內(nèi)火災(zāi)煙氣流動的模擬，是從安全性的角度考慮對消防員的影響，更加全面地評估服裝的防護性能，減少安全隱患。

4 結(jié) 語

熱防護服裝系統(tǒng)是火場或高溫環(huán)境中對消防員等作業(yè)人員產(chǎn)生保護的唯一屏障，是服裝領(lǐng)域的研究熱點之一。熱防護服裝性能數(shù)值模擬的目的并不是完全取代所有真實試驗，而是為熱傳遞和燒傷過程提供有效的信息，從而減少實驗的人力耗費，并指導(dǎo)進一步探索的方向。本文主要對20世紀90年代以來，國內(nèi)外熱防護服裝領(lǐng)域的數(shù)值模擬進行了回顧。以數(shù)值模擬的應(yīng)用方法為主線，總結(jié)了數(shù)值模擬在服裝熱防護性能測評中的作用；結(jié)合相關(guān)研究和技術(shù)的發(fā)展，概述了該方向未來的發(fā)展趨勢。目前的研究主要集中在小規(guī)模的織物層次，解決織物-空氣層-皮膚系統(tǒng)的熱濕傳遞問題，并有少量對于服裝層次或全尺度燃燒假人層次的研究。未來的研究趨勢將是在完善熱濕傳遞模型的基礎(chǔ)上，開展全尺度的數(shù)值模擬；結(jié)合人體熱調(diào)節(jié)等模型對防護服進行多方位的評價，從而更好地得到反饋并應(yīng)用在熱防護服裝的優(yōu)化方面，充分體現(xiàn)數(shù)值模擬研究的實用價值。

[1] 何雅玲, 王勇, 李慶. 格子Boltzmann方法的理論及應(yīng)用[M]. 北京: 科學(xué)出版社, 2009: 1-2. HE Yaling, WANG Yong, LI Qing. Lattice Boltzmann Method: Theory and Application[M]. Beijing: Science Press, 2009: 1-2.

[2] KAYNAKLI O, KILIC M. Investigation of indoor thermal comfort under transient conditions[J]. Building and Environment, 2005, 40 (2): 165-174.

[3] KILIC M, KAYNAKLI O, YAMANKARADENIZ R. Determination of the required core temperature for thermal comfort with steady-state energy balance method[J]. Heat Mass Transfer, 2006, 32 (2): 199-210.

[4] GHAZY A, BERGSTROM D J. Numerical simulation of transient heat transfer in a protective clothing system during a flash fire exposure[J]. Numerical Heat Transfer, 2010, 58(9):702-724.

[5] TORVI D A. Heat transfer in thin fibrous materials under high heat flux conditions[D]. Edmonton: University of Alberta, 1997: 1-277.

[6] MELL W E, LAWSON J R. A heat transfer model for firefighters′ protective clothing[J]. Fire Technology, 2000, 36:39-68.

[7] GHAZY A, BERGSTROM D J. Influence of the air gap between protective clothing and skin on clothing performance during flash fire exposure[J]. Heat Mass Transfer, 2011, 47(10):1275-1288.

[8] SONG G W, BARKER R L, HAMOUDA H, et al. Modeling the thermal protective performance of heat resistant garments in flash fire exposures[J]. Textile Research Journal, 2004, 74(12):1033-1040.

[9] CHITRPHIROMSRI P. Modeling of thermal performance of firefighter protective clothing during the intense heat exposure[D]. USA: North Carolina State University, 2004: 1-150.

[10] CROWN E M, DALE J D, BITNER E. A comparative analysis of protocols for measuring heat transmission through flame resistant materials: capturing the effects of thermal shrinkage[J]. Fire and Materials, 2002, 26: 207-213.

[11] 朱方龍，王秀娟，張啟澤，等. 火災(zāi)環(huán)境下應(yīng)急救援防護服傳熱數(shù)值模擬[J]. 紡織學(xué)報, 2009, 30(4): 106-110. ZHU Fanglong, WANG Xiujuan, ZHANG Qize, et al. Numerical model of heat transfer in protective clothing for emergency rescue in firefighting[J]. Journal of Textile Research, 2009, 30(4): 106-110.

[12] SAWCYN C M J, TORVI D A. Improving heat transfer models of air gaps in bench top tests of thermal protective fabrics[J]. Textile Research Journal, 2009, 79: 632-644.

[13] Wikipedia the free encyclopedia. QuickBASIC [DB/OL]. [2013-09-28]. http://zh.wikipedia.org/wiki/Microsoft_QuickBASIC.

[14] GHAZY A, BERGSTROM D J. Numerical simulation of the influence of fabric′s motion on protective clothing performance during flash fire exposure[J]. Heat Mass Transfer, 2013, 49(6):775-788.

[15] SAWCYN C M J. Heat transfer model of horizontal air gaps in bench top testing of thermal protective fabrics[D]. Saskatoon: University of Saskatchewan, 2003:1-134.

[16] 王福軍. 計算機流體動力學(xué)分析-CFD軟件原理與應(yīng)用[M]. 北京：清華大學(xué)出版社，2012:4. WANG Fujun. Analysis of Computer Fluid Dynamics: The Theory and Application of CFD Software[M]. Beijing: Tsinghua University Press, 2012:4.

[17] PENNES H H. Analysis of tissue and arterial blood temperature in resting human forearm[J]. Journal of Applied Physiology, 1948, 1(2): 93-122.

[18] LIU J, XU C, LISA X X. New Thermal wave aspects on burn evaluation of skin subjected to instantaneous heating[J]. Transactions on Biomedical Engineering, 1999, 46(4): 420-428.

[19] STOLL A M, CHIANTA M A. Method and rating system for evaluation of thermal protection[J]. Aerospace Medical, 1969, 40: 1232-1238.

[20] HENRIQUES F C. Studies in thermal injuries: V: the predictability and the significance of thermally induced rate processes leading to irreversible epidermal injury[J]. Archives of Pathology, 1947, 43: 489-502.

[21] AUTRIQUE L, LORMEL C. Numerical design of experiment for sensitivity analysis-application to skin burn injury prediction[J]. Biomedical Engineering, 2008, 55(4): 1279-1290.

[22] 王春，方惠惠，黃顯青，等. 地鐵火災(zāi)數(shù)值模擬研究最新進展[J]. 消防技術(shù)與產(chǎn)品信息, 2013 (3):33-35. WANG Chun, FANG Huihui, HUANG Xianqing, et al. Proceedings of numerical study on subway fire[J]. Fire Technique and Products Information, 2013 (3):33-35.

[23] 耿鵬. 室內(nèi)火災(zāi)數(shù)值模擬[D]. 太原：中北大學(xué), 2008: 1-49. GENG Peng. Numerical simulation of indoor fire[D]. Taiyuan: North University of China, 2008: 1-49.

[24] JIANG Y Y, YANAI E, NISHIMURA K, et al. An integrated numerical simulator for thermal performance assessments of firefighters′ protective clothing[J]. Fire Safety Journal, 2010,10(45):314-326.

[25] INCROPERA F P, DEWITT D P, BERGMAN T L, et al. Fundamentals of Heat and Mass Transfer[M]. 6 th ed. Beijing: Chemical Industry Press, 2012: 176-180.

[26] DONGARRA J, FOSTER I, FOX G, et al. Review of Parallel Computation[M]. Beijing: Electronic Industry Press, 2005: 62-65.

[27] 徐煒民, 嚴允中. 計算機系統(tǒng)結(jié)構(gòu)[M]. 3版. 北京: 電子工業(yè)出版社, 2010: 336- 339. XU Weimin, YAN Yunzhong. Computer System Architecture[M]. 3rd ed. Beijing: Electronic Industry Press, 2010: 336-339.

[28] TANABE S, KOBAYASHI K, NAKANO J, et al. Evaluation of thermal comfort using combined multi-node thermoregulation (65MN) and radiation models and computational fluid dynamics (CFD)[J]. Energy and Buildings, 2002, 34(6): 637-646.

[29] LI F Z, LI Y. Effect of clothing material on thermal responses of human body[J]. Modeling and Simulation in Materials Science and Engineering, 2005, 13(6): 809-827.

[30] CHEN C J, HSIEH W D, HU W C, et al. Experimental investigation and numerical simulation of a furnished office fire[J]. Building and Environment, 2010, 45 (12): 2735-2742.

[31] WU J X, MIAO Q, WANG E Y, et al. Numerical study on indoor fire of leaked gaseous fuel[J]. Fire Safety Science, 2005(8): 681-691.

[32] TALUKDAR P, TORVI D A, SIMONSON C J, et al. Coupled CFD and radiation simulation of air gaps in bench top protective fabric tests[J]. International Journal of Heat and Mass Transfer, 2010, 53(1): 526-539.

Application of numerical simulation on performance evaluation of thermal protective clothing

TIAN Miao1,2,3, LI Jun1,2,3

(1.FashionandArtInstitute,DonghuaUniveristy,Shanghai200051,China; 2.ProtectiveClothingResearchCenter,DonghuaUniveristy,Shanghai200051,China; 3.KeyLaboratoryofClothingDesign&Technology,MinistryofEducation,DonghuaUniversity,Shanghai200051,China)

Evaluating the performance of thermal protective clothing properly is significant for public safety. Following experiment research, numerical simulation has become an important means in related fields. From the aspects of development and validation of models, this paper reviews the current research status, such as heat transfer models, skin burn predictive models and fire simulation. The development process, characteristic and insufficient points of representative models are concluded. According to the development of related technologies, trends on numerical simulation on performance evaluation of thermal protective clothing are predicted. Full-scale simulation will be carried out based on the improvement of heat and moisture model. All-round evaluation of protective clothing will be realized combined with human thermal regulation model. Through the application on the optimization of thermal protective clothing, the practical value of numerical simulation can be expressed adequately.

numerical simulation; thermal protective clothing; performance evaluation; skin burn; heat transfer model

10.13475/j.fzxb.201501015807

2013-11-12

2014-10-06

國家自然科學(xué)基金資助項目(51106022)；中央高?；究蒲袠I(yè)務(wù)費專項基金資助項目(14D110715117118)

田苗(1989—)，女，博士生。研究方向為服裝功能與舒適性。李俊，通信作者，E-mail: lijun@dhu.edu.cn。

TS 941.73

A