1. 引言
十九大以来,随着全球信息化工能产业的飞速发展,大量从事在冶金、机械制造工业一线的工作人员长期面临着高温环境的损害。高温带来的热危害会威胁其生命安全和身心健康,因此设计专用的服装已成为各国进行高温防护的重要手段,同时也是我国经济发展过程中保障生产人员生命财产安全,提高产能效力的重要举措 [1]。
热防服能够在高温环境下促进人体热量散发,防止热中暑、烧伤和灼伤等,长期以来,对于热防服的导热过程和热传导下厚度的研究一直是一个热门问题。早期的研究通过Pdepe算法的皮肤表层热传导推演较为准确地表达了空气层和单层的固态热传导模型 [2];Torvi课题组在考虑不同辐射场强的影响下,结合空气层厚度的变换情况,建立了基于冷却时差的外层皮肤灼伤和内部热导模型 [3];部分研究从应用数学角度考虑分数阶偏微分方程的优化求解,推导出了多层非稳态热传递的温度分布精确表达式 [4];同时,优化反演模型的提出进一步验证了热对流的单值性分析的有效性 [5]。纵观现有研究成果,在热防服导热模型的构建上已取得丰硕结论,但不足之处在于未充分考虑不同层之间边界条件对热传导的影响偏差,同时现代算法的模拟预测在特定工作环境中的最优厚度设计和分析等方面仍有提升空间。
因此本文主要考虑了多层热防服层间的导热微分关系和边界条件影响,利用热力学定律将热防服导热场量化为有内热源的多层圆筒壁非稳态导热过程,建立基于Logistic的非稳态导热模型,求解出各层间时间维度的温度分布;进一步利用有限差分法求出空间维度上各层间的温度分布。之后根据实际设计中对特定工作环境和外界温度条件要求,创新地提出暂态电路等效法,结合基尔霍夫方程 [6],建立非线性规划,得到了最优热防服厚度设计数值,同时结合Multisim进行仿真检验,进一步证明了算法模型的有效性。
2. 基于Logistic非稳态导热模型的温度分布
2.1. 热传导定性分析
通过实验研究已知热防服每层材料的参数值以及实验假人表面外层温度5400 s内的变化情况如图1,结合文献 [7],明确了初始温度为37℃的假人人体始终受到来自外界环境温度热量交换,同时热防服有三层不同的织物层,因此本文中将其合理的量化为通过多层圆筒壁导热的有内热源模型分析。
![](//html.hanspub.org/file/19-1542134x10_hanspub.png)
Figure 1. Temperature change outside the skin of dummy
图1. 假人皮肤外侧温度变化
基于微分思想,将织物层的物料微元化处理,根据能量守恒定律单位时间净导入微元体热量
加上微元体内热源生成的热量
应等于微元体焓的增加量
,得到
时间内,沿厚度x轴方向,经
表面导出的热量和微元体净热量表示为:
(1)
(2)
根据傅里叶变换,上式代入能量守恒方程式可得一维导热的微分方程:
(3)
2.2. 边界条件对导热系统温度分布影响
在非稳态导热过程中,由于各层热量传递不均匀,即环境(边界条件)对系统温度分布一直起作用,将导热过程表示为图2。为此合理假设热量在织物层I、II、III、空气层IV层之间传递时没有热对流和热辐射现象发生,且热对流和热辐射只发生在空气流体与I层的接触边界上。
![](//html.hanspub.org/file/19-1542134x19_hanspub.png)
Figure 2. Thermal protective clothing heat transfer simulation diagram
图2. 热防服导热模拟图
由于I层接触面上单位面积与周围流体之间产生对流换热,采用相似原则对传热过程的参数进行归类处理,引入无量纲数,即准数。无量纲温度:
,无量纲坐标:
,将导热方程和边界条件无量纲化处理表达式如下:
(4)
流体力学中,根据动能守恒定律,转化得到相似准数傅里叶准数Fo和毕渥数Bi,其模型表达为:
(5)
其中Fo表征了题中假设的多层圆筒壁导热的有内热源系统的导热性能和储热性能的对比关系,是给定系
统的动态特征量,Bi是Bi织物层内部导热热阻与物体表面对流换热热阻的比值,
是
三个参数的函数。
2.3. 模型构建和温度分布求解
根据图1拟合假人表层温度变化走势,通过多层圆筒壁的非稳态导热过程中,由于热量的转移和各层的物理参数影响,进程随着时间的推进,温度达到蓄热饱和后,趋于稳定值,结合文献 [8] 引入Logistic曲线作为非稳态导热的模型基础,考虑到各层热阻越大则最终温度越低,反之热阻越小最终温度越高,而毕渥数是衡量热阻大小的量,结合公式(4) (5),得到非稳态导热模型表达式为:
(6)
结合非稳态导热变换规律,模型中需要给予时间
一个衡量标准,根据传热学引入与非稳态相关的相似准数——傅里叶准数作为时间(
)的载体,改进模型。考虑到不同介质之间的热传递的边界条件,会产生能量的损耗,为此引入衰减指数
,对热能的损耗进行限定。得到改进后的模型表达式为:
(7)
多层圆筒壁存在内热源的导热模型反映的是各层在非稳态条件下的温度变化率,所以
,得到最终温度表达式为
。
i) 时间维度的温度分布
将假人外侧环境温度代入模型,利用matlab软件绘制得到各层材料中心温度,即随时间变化的各层温度分布情况如图3。
![](//html.hanspub.org/file/19-1542134x33_hanspub.png)
Figure 3. Temperature distribution of each layer
图3. 各层温度分布图
ii) 空间维度的温度分布
当该模型中沿厚度x的方向上有一定的温度差时,在厚度x方向就会有一定的热量传递。从宏观上看,单位时间内通过垂直x方向的单位面积的热量q与温度的空间变化率成正比 [9]。此时考虑时空维度上的非稳态导热过程的温度分布,通过有限差分法并结合matlab绘制时空维度温度分布图如图4。
![](//html.hanspub.org/file/19-1542134x34_hanspub.png)
Figure 4. Space-time temperature distribution
图4. 时空维度温度分布图
结论:1) 通过有限差分法求得的时间维度的各层温度分布走势与基于Logistic的非稳态导热模型求得的结果相一致,说明基于Logistic的非稳态导热模型对温度随时间规律的处理是合理的;2) 通过基于Logistic的非稳态导热模型求得的同一时刻各层中心点的温度与有限差分法求得的空间维度的各层温度分布相符,并且符合热力学中对服装材料参数特性的评估,进一步说明了模型的合理性和可靠性。
3. 基于基尔霍夫的等效暂态电路模型
3.1. 多层圆筒壁下电路类比分析
非稳态导热模型在热动力学层面考虑到实际情况中不同层表面情况的不稳定性和热学变量的复杂性,为更精确的解决问题,将一维导热与电路模型进行类比分析。对于多层圆筒壁模型如图5,设织物层的单位横截面积为s,单位体积的内热源为q,微元厚度在时间
的温度t,织物层与流体的换热系数为h,热流量为Q。
![](//html.hanspub.org/file/19-1542134x36_hanspub.png)
Figure 5. Multilayer cylindrical wall model
图5. 多层圆筒壁模型图
当环境温度为
,在任一时刻
,对织物层任一截面出的微元厚度
列出能量平衡式:
(8)
按傅里叶导热定律得到
,令
,
,
,
结合公式(8)求得
等效织物层微分方程和定解条件为:
(9)
3.2. 非稳态热传递等效电路模型建立
常见的
电路中,当外接电压源为恒定的直流电源时,电路存在暂稳态特性 [10];故将传热学中热量传递过程中热量与温度差的关系转换为电学中欧姆定律形式,类比为具有电流源,无电感的均匀传输线路,设单位长度传输线下电容
,电阻
,漏电导
,电流源
,外加电势源(理想恒压源)
,外电源电阻
,距起点x处传输线的电位u,在x截面的电流i。由电路的基尔霍夫定律(KCL)建立方程:
(10)
(11)
(12)
将(11)、(12)二式代入(10)有:
(13)
对于偏微分方程组(13),利用Wolfram Mathematica进行精确解的计算,得到结果如下:
(14)
且
为毕渥准数,则由
和
得到特征方程如下:
(15)
4. 基于非线性规划的最优厚度研究
4.1. 第三边界条件下单层最优厚度求解
结合热防服层间关系,单II层材料厚度的研究可看做是第三边界条件和类比电路模型约束下的非线性规划问题。考虑实际生产情况下 [11],当环境温度为65℃、III层的厚度为3.6 mm、IV层的厚度为5.5 mm时,确保工作60分钟时,人体皮肤外侧温度不超过47℃,且超过44℃的时间不超5分钟。确定II层的最优厚度。
当织物层厚度增大,其热阻增大,相对导热能力变差,热防服织物系统的防热能力增强,寻找满足约束条件下的最小厚度值即为最优厚度,由此建立导热模型的约束条件S.T [12],非线性规划函数如下:
(16)
此时,相对应的类比电路的等效模型中
,则(14)式变形为:
(17)
这即是多层圆筒壁存在内热源模型在第三边界条件下的受热解 [13],其类比电路为图中电流源
及漏电导开路(即
),利用电学软件Multisim仿真如图6:
![](//html.hanspub.org/file/19-1542134x68_hanspub.png)
Figure 6. Multisim simulation diagram
图6. Multisim仿真图
将导热模型的复杂非线性规划转换至类比电路模型中,结合公式(14)和上述基尔霍夫方程,求得
。即当专用防护服材料参数第IV层材料参数厚度和环境温度等约束条件确定的情况下,为带来长久工作的持续保护,确保工作60分钟时,满足约束条件
,
、
、
下,第II层最优厚度为8.7 mm尚可。
4.2. 第一、三边界条件下多层最优厚度求解
考虑满足当环境温度为80℃时、III层的厚度为3.6 mm,确保工作30分钟时,假人皮肤外侧温度不超过47℃,且超过44℃的时间不超过5分钟,确定II层和IV层的最优厚度。由于空气层和织物层的层间密度与热传导率相差甚远 [14],且第IV层空气层与相邻的第III层织物层以及皮肤表层之间无热对流现象,在织物层未知的基础上需要考虑第IV层与相邻层之间的第一边界条件,为此建立新的非线性规划如下:
(18)
此时,
,由(15)可知
,得到
、
、
,由此可
将(14)式改为:
(19)
这即是多层圆筒壁存在内热源模型在第一类边界条件下非对称受热解,根据假设已知,假人处在环境温度中时,应该为对称受热,此时将模型进行修改,令
。则(14)式变为:
(20)
此时类比电路时令理想电势源
,利用电学软件Multisim仿真如图7。将导热模型的复杂非线性规划转换至类比电路模型中,结合公式(14)和基尔霍夫方程,求得
,
。即当热防护服设计中的决策变量增加时,充分考虑人体温度的物理条件和空气层产生的单值性条件约束下,保证环境温度为80℃时,确保工作30分钟时,满足约束条件
,
、
、
下II层最优厚度为11.6 mm,IV层的最优厚度为4.5 mm。
![](//html.hanspub.org/file/19-1542134x90_hanspub.png)
Figure 7. Multisim simulation diagram
图7. Multisim仿真图
5. 总结
在以人为本的热能源高新产业安全发展背景下,热防服表层温度分布和有效防护厚度的研究成为了保护人员生命安全和推动突破复杂热传导难题的关键。本文依据热传导定律巧妙地将多层热防服结构抽象为多层圆筒壁模型以此研究非稳态导热过程,通过实验数据充分考虑边界条件的影响建立了基于Logistic的非稳态导热模型,结合有限差分法准确求得了时空维度下的不同结构层的温度分布。之后摒弃传统算法求解,类比
电路暂态特性与非稳态导热过程之间的相似性,巧妙地转换为基尔霍夫方程下的等效电路模型,结合外在约束条件建立非线性规划模型,得到了热防服厚度的最优解,模型的良好迁移性与有效性为企业的生产研发提供了新思路。
基金项目
国家自然科学基金资助项目(60362001);云南省自然科学基础基金资助项目(61461053)。
NOTES
*通讯作者。