文档介绍：高温作业专用服装设计探析
 
 
欧伟铭 皇甫晓相
摘 要：借2018年全国大学生数学竞赛A题的背景与数据，探讨了高温作业专用服装设计问题，对于人体与外界对流换热的情况，导出对应热传导方程及其初边值条件并利用差分法进行数值求解，分别得出热防护服内部温度随时间和空间的分布。从减轻防护服重量角度出发，解得在相关约束下热防护服各层的最优厚度。
关键词：热传导方程;差分法;广度优先;穷举法;初边值问题
：TB ：Adoi：.1672-
1 问题重述
（1）根据已知条件建立数学模型，计算温度分布。
（2）在环境温度、IV层厚度、工作时间确定，皮肤外侧的温度低于47℃，且高于44℃的时间不超过 5 分钟时，求II层最优厚度。
2 模型的建立与求解
问题一模型的建立与求解
问题一已知各层织物的厚度、环境温度、工作时间及皮肤外侧的测量温 度，要求计算温度分布。求解步骤如下：
（1）以热防护服 I 层与外界环境的边界为原点，热量传递方向为 x 轴建立坐标系。
（2）根据热传导方程建立各层织物的非稳态热传导方程并确定各层织物的边界条件建立温度分布模型。
（3）通过显式欧拉法逐步迭代，求出皮肤外侧温度关于空间位置及时间关系，从而得到温度分布。
（4）假人中心抽象为恒温热源，将皮肤内层环境看作流体，增加对流换热的情况考虑，来精确皮肤表层温度的分布。
在环境-热防护服-皮肤外侧系统中，原点为I层与外界环境的边界，x的正向为热量传递的方向。由于本题已假定热传递仅沿垂直于皮肤方向进行，可视为一维导热问题，故可忽略y、z方向分量。又假人体内的温度控制在37℃恒定，故为0。
因此，在本题中热传导方程可简化为：
a为热扩散率。
热传导方程组的构建：热量从环境中以热传导的方式由I层逐步向内层传递，每层织物材料均有两个导热界面。根结合初始、边界条件建立各层织物的热传导方程组。
I层热传导方程：因本题仅考虑热传导一种传递作用，由 0 时刻的初始温度，得到初始条件及边界条件，对于 I 层：
II 、III层热传导方程组：
II层既与I层相邻，又与III层相邻，且两种导热界面均只有热传导一种作用。而IV层是空气层，在狭小的空间内，我们可以忽略空气的流动，即不考虑热对流作用，则III层的导热情况与II层类似。对于II 、III层有：
其余层对于边界条件的处理基本相同，不再赘述。
在假人体内中心，存在37℃的恒温热源，以热对流形式向周围“吹着 37℃的热风”，使皮肤内层温度稳定在 37℃。
附件数据的分析：
用Excel 对附件5400组数据进行分析。1645s前，皮肤外侧温度随时间变化而上升;1645s后，皮肤外侧温度不随时间的变化而改变。故前者为非稳态导热过 程，后者为稳态导热过程。
皮肤边界节点离散方程的建立：
考虑皮肤表层边界及外侧部分。皮肤表层受到内部流体的热对流作用，对流换热系数为h。此时边界点 N 代表宽度为x/ 2 的元体。
对该元体应用能量守恒定律：
对流换热系数h的求解：
本题存在稳态和非稳态导热两个阶段，针对不同的阶段，对流换热系数具有不同形式：
以上两式分别为稳定导热和非稳态导热情况下对流换热系数的表达式。