典型胶合板的着火时间研究
典型胶合板的着火时间研究
叶茂昌，王金广
（厦门市消防支队，福建厦门361012）
摘 要：为更加简单方便地求得材料的着火时间，考虑材料对热的吸收率和热损，在经典着火理论的基础上对一维导热方程的边界条件进行修正，对方程进行离散处理并且用编程求解材料的表面温度。代入杨木胶合板和沙比利胶合板的热物性参数，计算不同的辐射强度下胶合板表面温度随着时间的变化，得出着火时间，与实验结果相差不大。一维导热方程分析材料的着火时间是可行性的。3mm的胶合板吸收率约为0.25。
关键词：胶合板；一维导热方程；着火时间；吸收率；热损
中图分类号：X913.4 ，TK121 文献标志码：B
文章编号：1009-0029（2017）01-0031-03
胶合板材料因美观舒适而被广泛运用于家具和装饰品中。但胶合板的火灾危险性较高，会提供大量的火灾荷载，一旦发生危险会造成严重的后果。而着火时间是火灾发生和发展的一个关键因素。国内外很多学者对胶合板的热力学性质、物理性质以及化学组分进行了分析和模拟，也有很多学者提出了各种材料的着火时间的模型。这些学者的工作都很深入而且细致，但是通用性不强而且比较复杂。笔者以锥型量热仪的实验数据为基础，用一维导热的控制方程研究常用木材的着火时间来得到一个通用性强而且简单的计算公式，为火灾科学基础研究提供数据和理论支持。
1 着火模型
1．1 经典着火理论和修正
在恒定的辐射热流下热厚型材料内部温度的变化是由一维热传导方程控制的，如式（1）～式（3）所示。
假设材料初始的温度和环境温度一致，且不考虑材料各个位置热物性的差异。辐射热流的照射下上边界有热流的流入，大小假定为锥型量热仪标定的辐射功率。下边界由于热厚型材料在着火前热穿透厚度没有达到材料的厚度，可认为在材料着火前没有热流的流出。
由上述方程，可以推导出着火时间和辐射热流的关系，如式（4）所示。
实际上锥型量热仪的辐射功率是用热流计测定的，把这个标定的数值作为经典理论中流入材料热流，与实际流入材料上边界的热流是有些差异的，因为材料不可能对热流完全吸收，同时当加热时间过长且材料比较薄的情况下是有热流从下边界流出的，因此需要考虑材料对热的吸收率β和对外的辐射热损
。即上边界条件如式（5）所示。
修正后同样可以推导出着火时间与辐射热流的关系，如式（6）所示。
1．2 方程离散处理
为了使规律更加简单而便于运用，可以认为材料表面温度达到点燃温度就立即着火，因此需要由上述导热方程求解出材料表面温度随时间的变化。直接求数学解析式比较麻烦，可以通过离散方法来求解数值解。
对时间项进行向后差分，如式（7）所示。
对空间项进行中心差分，如式（8）所示。
左右两边整理可得式（9）。
合并相同项可得式（10）。
上述方程只适用于第二个到倒数第二个点，还需补充离散方程。引入边界条件，如式（11）、式（12）所示。
上述方程式全隐式格式，需要联立所有的方程组依据第n层时间层的温度数据来求解第（n＋1）层的温度数据。矩阵如式（13）所示。
1．4 数值求解
用三对角阵求解上述数值结果是一种很好的解法。离散后可以通过编程进行求解计算，然后输出每个时间步长的材料表面温度。如式（14）的矩阵。
可以通过编程迭代计算每个节点每个时间步长的温度数据。
2 实验数据拟合
有关胶合板的锥型量热燃烧特性已经有许多学者做了研究工作，实验数据规律基本一致，笔者着重研究用简单一维导热方程预测材料着火时间的可行性，故在此引用其他研究工作者的数据进行计算。对杨木和沙比利材质的胶合板进行了锥型量热燃烧实验，所做的实验研究和笔者所需要的相符，得到着火时间如表1所示。
调研文献资料可知，胶合板的导热率平均值为k＝0.125 W／（ｍ·Ｋ）、比热容c＝2.5 kJ／（kg·K），杨木胶合板的密度为388 kg／㎥，沙比利胶合板的密度为568 kg／㎥。着火温度约为300℃。为方便作图和发现胶合板点燃的相关规律，需将方程简化为可用直线拟合的表达式，如式（15）所示。
拟合作图结果，拟合关系如式（16）、式（17）所示。
通过以上直线拟合公式，可求得实验中两种胶合板对辐射热流的吸收率和热损，如式（18）、式（19）所示。
从以上结果看，胶合板对于热的吸收率比较低，这是因为王磊实验中材料厚度只有3 mm，在长时间的加热过程中热流已经从底部有所流出，再加上自身的热辐射，因此热流的吸收较少。从材料的本身考虑，沙比利的颜色比杨木要深，对热流的吸收和耗散比杨木胶合板高，数值大小是较为合理的。
３ 数值计算结果分析
在辐射热流的作用下，假设胶合板表面温度达到着火温度即被点燃，不考虑内部热解和碳化等复杂的反应机理。为了进行数值计算，需要对时间步长和空间步长进行划分。对于时间步长，可以取Δt＝1ｓ，而空间步长把3 mm厚的材料划分为101个节点，即每个节点间隔为＝0.03 mm，代入式（15）可以计算得到：
在不同的辐射强度下，可以根据吸收率和热损计算出由边界条件流入胶合板的热流，如表2所示。
把以上数据分别输入程序，计算可得胶合板表面温度随着时间的变化情况，结果如图1、图2所示。
依照前面所描述的简化考虑，可认为胶合板表面温度达到着火温度即被点燃。因此，以温度纵坐标为300℃作一条直线与各个曲线相交，可得杨木胶合板和沙比利胶合板的着火时间，如表3所示。
对比分析实验测量值和数值计算结果，发现计算结果和实验值吻合较好，只有几秒钟的误差，这充分说明了考虑热吸收和热损后用一维导热方程预测材料的着火时间是非常可靠的。造成数值计算误差的主要因素是材料的热物性参数没有精确测量而是采用平均的数据带入程序计算的，如果这些热物性参数更加精确，所得的结果也会更加接近实验测量的着火时间。
4 结 论
（1）以一维导热方程分析材料的着火时间能够很简单简洁而又方便地得到与实验值相差不大的结果。
（2
上一篇 : 锂电池热失控火灾与变动环境热失控实验
下一篇 : 非连续介质热解挥发分对其着火的影响
