余瑩瑩，徐曉楠（中國人民武裝警察部隊學院，河北廊坊，065000）

?

木質(zhì)類人造板材炭化痕跡實驗及數(shù)值重構研究

余瑩瑩＊，徐曉楠
（中國人民武裝警察部隊學院，河北廊坊，065000）

摘要：利用錐形量熱儀測定不同輻射熱通量和不同燃燒時間下木質(zhì)類人造板材的炭化深度，建立木質(zhì)類人造板材炭化痕跡的數(shù)學模型，并嵌入到FDS（Fire Dynamics Simulator）源程序中，與自身功能相結合，實現(xiàn)重構木質(zhì)類人造板材炭化痕跡的功能。同時實驗研究了不同油盤尺寸、不同燃燒時間、不同火源位置和不同材料種類對炭化痕跡的影響，并利用新編譯的FDS程序?qū)嶒炦M行重構，對比重構結果和實驗結果，證明該程序可以用于重現(xiàn)木質(zhì)類人造板材的炭化痕跡。

關鍵詞：人造板材；炭化痕跡；數(shù)值重構；FDS；火災調(diào)查

0　引言

木質(zhì)類人造板材廣泛應用于室內(nèi)裝修裝潢和家具的制造過程中，受到火場作用后，形成的炭化痕跡是火災現(xiàn)場中常見的痕跡物證之一，可直接或間接證明火災發(fā)生的時間、起火點位置、起火原因和蔓延路線等［1，2］。但由于火災燃燒過程的復雜性、火災本身的破壞性、不同火場的特殊性等原因，給火災調(diào)查工作帶來了極大的困難［3］，而火災數(shù)值重構技術作為火災調(diào)查工作中的一種新技術，能夠直觀準確地再現(xiàn)火災過程以及得到火災現(xiàn)場勘驗無法得到的數(shù)據(jù)，發(fā)展前景十分廣闊。國外，關于火災數(shù)值重構的研究主要是針對火災蔓延和煙氣流動過程，Madrzykowski等人［4］對1999年5月30日華盛頓某街道別墅火災進行重構，結果表明，火災中地下室的玻璃推拉門打開，使氧氣進入了高溫受限防火分區(qū)，導致60s后發(fā)生轟然，產(chǎn)生的煙氣高速沖向地下室樓梯。國內(nèi)，不僅研究了火災蔓延和煙氣流動過程的數(shù)值重構，還研究了火災痕跡的數(shù)值重構，同課題組的韓［5］研究了PVC板熱蝕痕跡的數(shù)值重構、吳［6］研究了防火板煙熏痕跡的數(shù)值重構、施［7］研究了石膏板脫水痕跡的數(shù)值重構，但木質(zhì)類人造板材炭化痕跡的數(shù)值重構還未研究，且炭化痕跡的形成機理與前幾種痕跡均不同。因此，本文以木質(zhì)類人造板材炭化痕跡為研究對象，通過實驗和理論相結合的方法建立數(shù)學模型，并將數(shù)學模型嵌入到FDS源程序中，實現(xiàn)對木質(zhì)類人造板材炭化痕跡的重構，同時，設計了實驗驗證新編譯的FDS程序的有效性，將該技術用于預測火場中木質(zhì)類人造板材炭化痕跡，與火場中的其它證據(jù)相印證，為火調(diào)人員反推起火點提供幫助。

1　木質(zhì)類人造板材炭化痕跡數(shù)學模型的建立與實現(xiàn)

1．1 數(shù)學模型的建立

1．1．1 木質(zhì)類人造板材受熱炭化實驗

按照ISO5660將刨花板、細木工板或密度板裁成尺寸為100mm×100mm的試樣進行錐形量熱實驗，為了減小模型建立時產(chǎn)生的誤差，選擇炭化深度差別較大的時間參數(shù)進行實驗，具體安排見表1，實驗過程中，記錄試樣表觀狀態(tài)隨時間的變化，實驗結束后，測量試樣的炭化深度（以此值表征試樣的炭化程度）。

表1　實驗安排表Table 1 The schedule of experiment

1．1．2 結果分析

火災烈度是衡量火災破壞性的指標，對火災中各構件受損程度具有重要的影響，可近似地用接收的凈輻射熱通量（q·″net）和暴露時間（t）的乘積來計算［8］。當火災烈度為0MJ／m2，炭化深度為0mm，對實驗數(shù)據(jù)進行曲線擬合時應過原點，擬合曲線見圖1，相應的數(shù)學關系式見式（1）。

式中：D－炭化深度，mm；S－火災烈度，MJ／m2。

由式（1）可知，炭化深度與火災烈度均符合D＝a·Sb＝a·（q·″net·t）b的關系，其中a，b是與材料性質(zhì)有關的系數(shù)，q·″net是隨時間t的變化量，因此，木質(zhì)類人造板材炭化痕跡的數(shù)學模型用式（2）表示：

式中：k－與實驗尺度有關的系數(shù)。根據(jù)實驗過程記錄的材料表面開始炭化的時間，求得刨花板、密度板、細木工板形成炭化痕跡的邊界閾值分別為：火災烈度為1．71MJ／m2、1．67MJ／m2、1．43MJ／m2時，炭化深度為1．71mm、1．62mm、1．68mm。

圖1　炭化深度與火災烈度關系的擬合圖Fig．1 Fitting chart of relationship between carbonization depth and fire severity

1．2 數(shù)學模型的實現(xiàn)

在Visual Studio 2010操作環(huán)境下，利用Fortran語言［9］對FDS 6．0．1源程序（共有32個文件）［10］中的read．f90、dump．f90、wall．f90等文件進行修改并嵌入式（2），經(jīng)過不斷調(diào)試，式（2）與FDS源程序成功耦合，編譯成了新的FDS程序，新編譯的FDS程序能夠計算并顯示木質(zhì)類人造板材炭化痕跡。

2　木質(zhì)類人造板材炭化痕跡實驗

2．1 實驗材料

實驗用刨花板、密度板或細木工板作為壁面材料，購自于廊坊市建材市場，裁成實驗所需尺寸1200mm×800 mm，厚度分別為16 mm、8 mm、15mm；正方形油盤由廊坊市捷強五金加工廠制成，高度為50mm；燃料為93汽油。

2．2 實驗裝置

先將石膏板固定在搭建好的實驗臺上，再將人造板材固定在石膏板上，在實驗臺底部放置電子天平記錄燃料質(zhì)量的變化（用于設定模擬中火源的熱釋放速率），油盤放置在電子天平上用于盛放汽油，用計時器記錄燃燒時間，用攝像機記錄炭化痕跡的形成過程，用照相機記錄最終形成的炭化痕跡，用卷尺測量炭化痕跡的最大寬度、高度，并沿著炭化痕跡中心位置從下至上每隔10cm處，用炭化深度測量儀測量一次炭化深度，實驗裝置見圖2。

2．3 實驗步驟

圖2　實驗裝置圖Fig．2 Experimental apparatus

由式（2）可知，木質(zhì)類人造板材上形成的炭化痕跡受壁面接收的凈輻射熱通量、燃燒時間、材料性質(zhì)等因素的影響，而壁面接收的凈輻射熱通量與油盤尺寸和位置有關，通過控制變量法研究各影響因素對炭化痕跡的影響。

（1）不同油盤尺寸。油盤邊長分別為13cm、15cm、20cm、25cm、30cm，位于底部中央距壁面2．0cm，為保持燃燒時間的一致性，待燃燒時間到120s時，人為將火熄滅，壁面材料為刨花板，實驗編號為1～5。

（2）不同燃燒時間。燃燒時間分別為100s、120s、140s、160s、180s，為了便于觀察形成的炭化痕跡，選取邊長為20cm的油盤，位于底部中央距壁面2．0cm，壁面材料為刨花板，實驗編號為6～10。

（3）不同油盤位置。油盤分別位于底部中央距壁面1．5cm、2．0cm、2．5cm和距離墻角2．0cm，同樣選取邊長為20cm的油盤，燃燒時間為120s，壁面材料為刨花板，實驗編號為11～14。

（4）不同材料種類。壁面材料分別為刨花板、密度板和細木工板，油盤邊長為20cm，位于底部中央距壁面2．5cm，燃燒時間為120s，實驗編號為15 ～17．

3　木質(zhì)類人造板材炭化痕跡數(shù)值重構

利用新編譯的FDS程序重構人造板材炭化痕跡實驗，計算空間有兩個面全部敞開，與實驗的設置保持一致。根據(jù)網(wǎng)格獨立性測試及時間步長的調(diào)整［11，12］，在1．2m×0．8m×1．2m的模擬計算空間內(nèi)共設了60×48×60個網(wǎng)格。材料的熱物性參數(shù)設置參照NIST發(fā)布的典型材料熱物性參數(shù)和《SFPE消防工程手冊》，模擬時間設置與實驗時間相同，實驗室門是關閉的，風速設為0m／s，實驗裝置模擬圖見圖3。

圖3　實驗裝置模擬圖Fig．3 Simulation diagram of experimental apparatus

4　結果與討論

4．1 油盤尺寸的影響

matchShapes值由三個與Hu矩有關的函數(shù)計算得到，用于比較兩個圖像形狀的相似性，一般認為，當matchShapes值小于0．4，兩個圖像在形狀上較相似［13，14］。表2是1～5實驗和重構的炭化痕跡圖，matchShapes值分別為0、0．361291、4．44089×10－16、0．321281和0．130638；實驗得到的痕跡的最大寬度分別為0．0cm、12．5cm、15．0cm、15．3cm、18．1cm，最大高度分別為0．0cm、27．4cm、46．1cm、68．9cm、85．4cm，重構得到的痕跡的最大寬度分別為0．0cm、11．6cm、14．9cm、15．7cm、17．4cm，最大高度分別為0．0cm、27．8cm、46．7cm、70．1cm、86．0cm；實驗和重構得到的痕跡的最大寬度、高度均與油盤邊長呈線性關系，見式（3）和（4），最大高度的增長速率快于寬度的，錐形形狀越來越明顯；圖4為1～5實驗和重構的炭化深度值，隨著油盤尺寸的增大，炭化深度值增大，且炭化深度值在第2個或第3個探測開始出現(xiàn)，主要原因油盤與壁面存在一定距離，且底部有總高度為16cm的電子天平和隔熱材料；證明了新編譯的FDS程序可以較好的重現(xiàn)炭化痕跡的形狀、尺寸、隨油盤尺寸的變化規(guī)律、炭化深度等特征。

式中：W－痕跡的最大寬度，cm；H－痕跡的最大高度，cm；L－油盤邊長，cm。

表2　1～5實驗和重構的炭化痕跡圖Table 2 Experimental and reconstructed carbonization trace of 1～5

圖4　1～5實驗和重構的炭化深度值Fig．4 Experimental and reconstructed ofcarbonization depth of 1～5

4．2 燃燒時間的影響

表3　6～10實驗和重構的炭化痕跡圖Table 3 Experimental and reconstructed carbonization trace of 6～10

圖5　6～10實驗和重構的炭化深度值Fig．5 Experimental and reconstructed of carbonization depth of 6～10

式中：t－燃燒時間，s。

4．3 火源位置的影響

11～14實驗和重構的炭化痕跡圖見表5，matchShapes值分別為0．31095、4．44089×10－16、1．07429×10－2和6．81211×10－2；11～14實驗得到的痕跡的最大寬度分別為13．4cm、15．0cm、3．1cm、15．3cm，最大高度分別為63．4cm、46．1cm、9．2cm、90．4cm；重構得到的痕跡的最大寬度分別為13．7cm、14．9cm、2．9cm、15．4cm，最大高度分別為64．1cm、46．7cm、8．9cm、90．6cm；對比11～13的結果，隨著火源與壁面距離的增大，痕跡尺寸逐漸減小，顏色也逐漸變淺，對比12和14的結果，壁面受墻角的交叉輻射效應，接收的累積熱量增多，痕跡尺寸明顯增大，且顏色明顯變深；圖6為11～14實驗和重構的炭化深度值，火源距壁面越近，炭化深度值越大，火源位于墻角，炭化深度明顯增大，炭化深度值仍開始出現(xiàn)于第2個或第3個探測點．

表4　11～14實驗和重構的炭化痕跡圖Table 4 Experimental and reconstructed carbonization trace of 11～14

4．4 材料種類的影響

圖6　11～14實驗和重構的炭化深度值Fig．6 Experimental and reconstructed ofcarbonization depth of 11～14

15～17實驗和重構的炭化痕跡圖見表5。 matchShapes值分別為1．07429×10－2、0．327619和0．234081；15～17實驗得到的痕跡的最大寬度分別為3．1cm、12．1cm、17．9cm，最大高度分別為9．2cm、45．1cm、71．9cm；重構得到的痕跡的最大寬度分別為2．9cm、12．3cm、16．5cm，最大高度分別為8．9cm、45．4cm、72．3cm；根據(jù)1．1．2的結果，相同的火災烈度，炭化深度由小變大的順序為刨花板、密度板、細木工板，穩(wěn)定性的順序相反，15～17重構和實驗的結果表明，炭化痕跡的尺寸隨材料穩(wěn)定性減弱而逐漸增大、顏色逐漸加深；圖7為15～17實驗和重構的炭化深度值，材料穩(wěn)定性越好，炭化深度值越小，炭化深度值仍開始出現(xiàn)于第2個或第3個探測點。

表5　15～17實驗和重構的炭化痕跡圖Table 5 Experimental and reconstructed carbonization trace of 15～17

圖7　15～17實驗和重構的炭化深度值Fig．7 Experimental and reconstructed of carbonization depth of 15～17

5　結論

（1）炭化痕跡的最大寬度、高度均與油盤邊長呈線性關系，且高度的增長速率快于寬度的；炭化痕跡的最大寬度與燃燒時間呈Logistic的關系，最大高度與燃燒時間呈Allometric1的關系，隨燃燒時間的增加，高度明顯增加，而寬度幾乎沒有變化。

（2）火源位于中間位置，隨著與壁面的距離增大，炭化痕跡尺寸、炭化深度都逐漸減小，顏色也逐漸變淺；火源位于墻角位置，炭化痕跡尺寸、炭化深度有明顯的增大，且顏色明顯變深。

（3）三種材料的穩(wěn)定性順序為：刨花板＞密度板＞細木工板，炭化痕跡的尺寸、炭化深度隨材料穩(wěn)定性減弱而逐漸增大，顏色逐漸加深；

（4）利用新編譯的FDS程序?qū)嶒炦M行重構，分別從痕跡的形狀、尺寸、隨影響因素的變化規(guī)律、炭化深度四個方面對比重構結果和實驗結果，發(fā)現(xiàn)兩者具有較好的一致性，驗證了新編譯的FDS軟件的有效性，將該軟件應用于消防工作，不僅可以作為一種定量手段來描述炭化痕跡的形成過程，還可以用來反推起火點位置，較好地輔助火災調(diào)查工作。

參考文獻

［1］Georgali B，Tsakiridis PE．Microstructure of fire－damaged concrete：a case study［J］．Cement and Concrete Composites，2005，27（2）：255－259．

［2］Setién J，et al．Evolution of mechanical behavior in a structural steel subjected to high temperatures［J］．Engineering Failure Analysis，2002，9（2）：191－200．

［3］張偉．淺談數(shù)值重構技術在火災調(diào)查中的應用［J］．科技視界，2012，（25）：373－374．

［4］Madrzykowski D，Vettori RL．Simulation of the dynamics of the fire at 3146Cherry Road NE，Washington D C，May 30，1999［DB／OL］．http：／／fire．nist．gov／bfrlpubs／fire00／PDF／f00043．pdf，2000－04．

［5］韓娜．數(shù)值重構技術在火災調(diào)查中的應用研究［D］．廊坊：中國人民武裝警察部隊學院，2012．

［6］吳迪．利用數(shù)值重構技術重現(xiàn)煙熏痕跡的應用研究［D］．廊坊：中國人民武裝警察部隊學院，2013．

［7］施照成．紙面石膏板受熱痕跡的實驗與數(shù)值重構研究［D］．廊坊：中國人民武裝警察部隊學院，2014．

［8］Chu Nguong．Calcination of gypsum plasterboard under fire exposure［EB／OL］．http：／／www．civil．Canterbury．a(chǎn)c．nz／fire／Pdfreports／Chu04．pdf，2014－02－11．

［9］陳科．FORTRAN完全自學手冊［M］．北京：機械工業(yè)出版社，2009．

［10］McGrattan K，et al．Fire dynamics simulator user’s guide （sixth edition）［M］．Washington DC：National Institute of Standards and Technology，2013．

［11］Hartzell GE．Engineering analysis of hazards to life safety in fires：the fire effluent toxicity component［J］．Safety Science，2001，38（3）：147－155．

［12］Fire Code Reform Center Limited．Fire engineering guidelines［M］．1st ed．Sydney：Fire Code Reform Center Limited，1996：17－21．

［13］丁明躍．不變矩算法研究［J］．數(shù)據(jù)采集與處理，1992，7 （1）：1－9．

Keyword：Man－made wood board；Carbonization trace；Numerical reconstruction；FDS；Fire investigation

［14］Hu MK．Visual pattern recognition by moment invariants ［J］．IRE Transactions on Information Theory，1962：179－187．

Experiment and numerical reconstruction of carbonization trace of man－made wood board

YU Yingying，XU Xiaonan
（Chinese People’s Armed Police Force Academy，Langfang 065000，China）

Abstract：Carbonization depth of man－made wood board is determined under different radiation heat fluxes and different burning time using the cone calorimeter．A mathematical model of man－made wood board carbonization trace is established，and is embedded into FDS（Fire Dynamics Simulator）source program．By this means，the man－made wood board carbonization can be reconstructed．The influence of pool size，burning time，fire source position and wood species on the carbonization trace is also studied．By comparing the results of experiment and reconstruction，the newly compiled FDS is justified to be able to reconstruct the carbonization depth of man－made wood board．

通訊作者：余瑩瑩，E－mail：994673524＠qq．com．

作者簡介：余瑩瑩（1991．7－），女，漢，在讀研究生，主要從事火災調(diào)查、火災數(shù)值模擬方向的學習和研究。

基金項目：公安部應用創(chuàng)新計劃項目（2011YYCXWJXY134）；中國人民武裝警察部隊學院優(yōu)秀碩士學位論文培育工程資助（201403）。

收稿日期：2015－01－06；修改日期：2015－02－04

DOI：10．3969／j．issn．1004－5309．2015．01．06

文章編號：1004－5309（2015）(－)0040－07

文獻標識碼：A

中圖分類號：X928．7；X915．5