1. 引言

改革开放以来，我国的民航事业得到了快速发展，飞机数量直线上升，航空运量持续快速增长，航线网络不断扩大，机队运输能力显著增强，机场、空管等基础设施建设取得了重大进展。航空运输在实施“一带一路”战略、实现小康社会中发挥的作用将越来越大。据预测，在未来20年时间里，中国国内航空客运周转量年平均增长率将达7.1%。到2033年，中国国内航空市场将占同期全球航空客运量的11.9%。由此飞机的维修量急剧上升，并对飞机维修车间的室内环境提出了更高的要求，比如，飞机退漆、喷漆是飞机维修和定期检修的重要项目，飞机在退漆和喷漆时都需要适宜的温度条件。如果达不到喷漆的适宜温度或者飞机表面温度场变换较大，都将导致飞机退漆困难、喷漆不均匀、漆干后飞机表面不同部位存在色差 [1] [2] 。因此，采用先进的技术手段，创造适宜的飞机维修车间室内环境对于提高飞机维修水平，延长飞机寿命均具有重要的意义。

据了解，地板辐射采暖可以使供暖房间内底部温度较高，这是由于底层空气靠近辐射采暖板，并随着高度的上升温度逐渐降低，但下部热空气的上升会使房间顶部温度较高。根据地板辐射采暖这种分层供暖的特性，人员活动在2 m以下区域范围内的使用较好，在欧洲寒冷地区这种特性在广场、足球场等得到了很好的利用。由于飞机维修车间是属于高大空间，辐射供暖适用。

2. 机库辐射供暖模型的建立

2.1. 机库基本情况

本文以山东某飞机维修车间为例来进行分析。该机库由机库维修大厅、北侧附楼、东侧办公区组成，其总建筑面积为25,300 m² 。其中机库维修大厅建筑面积为9600 m²，长度为120 m，宽80 m，网架下弦高25 m。北侧附楼地下一层、地上3层，建筑面积为1500 m²。东侧办公建筑面积14,200 m²。该维修机库主要承担飞机C检及以下各级定检和航线的维护等任务，可同时容纳三架B747飞机入库维修。图1为机库大厅首层平面图 [3] 。

2.2. 供暖空调系统

机库维修大厅供暖系统：低温地板辐射配合暖风机供暖，供回水温度为，室内设计温度为14℃，室内体感温度可以达到16℃，供暖设计热负荷为2100 kw。

附楼供暖系统：使用自备锅炉房提供的热水进行供暖，其热负荷为1300 kw。

机库维修大厅空调系统：设置2套直接蒸发式空调，空调室外机设在附楼屋顶，室内机设在地下室空调机房内，空调管通过通行管沟接至空调地井，用专用接头通过软管接至飞机内部人员工作区附近。

在油箱和地沟等位置设置了机械通风系统用于快速排放出在检测和维修过程中散发出对人体和环境有害的物质和气体，保证检修工作的顺利进行。

附楼空调系统：夏季制冷使用多联机，机库工作人员使用可移动式空调器，其制冷室外机部分设在附楼屋顶上面。

2.3. 机库几何模型的建立

为研究机库辐射供暖条件下温度场、速度场建立几何模型，尺寸为的长方体大空间。机库大门尺寸为。

2.4. 湍流模型

湍流流动是在自然界广泛存在的一种流动方式也是最为常见的一种流动方式，其流动过程中的换热方式是对流换热。飞机维修车间的流态经过计算属于紊流。所以本文选择的模型是三维紊流不可压缩模型 [4] [5] 。

经分析研究，飞机维修车间内质量、能量、动量守恒控制方程的通用形式如下 [6] [7] 。

(1)

式中：为非稳态项；是对流项；为扩散项，为通用变量，可以代表，，，等求解变量；为广义扩散系数；为广义源项。

表1中，为空气的有效粘性系数，且；为空气的层流粘性系数；为空气的湍流粘性系数；、、分别为方向上的速度；为空气的压力；h为空气的焓；分别为方向上的重力加速度；为空气的参考密度；为h的当量普朗特数；为单位体积的发热量。

下面是辐射积分方程的具体形式：

(2)

式中是位置向量，是散射方向，是沿程长度，是吸收系数，是折算系数，是散射系数，是斯蒂芬–玻耳兹曼常数(W/m²∙K⁴)，是辐射强度，是当地温度，是相位函数，是空间立体角。

对通用方程式(1)的源项进行线性化处理， ()，则离散方程为如下形式：

(3)

本文采用Reynolds时均方程的湍流数值模拟方法，湍流模型采用的数学表达式如下：

(4)

式中，表示速度，表示离墙体的最近距离。

将室内空气湍流流动的控制微分方程(1)表示为：

(5)

式中，代表流动的速度、焓(温度)等物理量。

3. 基于辐射供暖的机库热环境分析

地板辐射采暖与室内环境是通过对流和辐射传热进行换热，室内空气之间的导热可以基本忽略。辐射换热量主要取决于维护结构表面、加热板、人体及室内飞机等设备的表面温度，各表面的几何形状，相对位置及其辐射黑度等特性。对流换热量主要取决于地板附近空气对流作用的强弱。值得注意的是，人体与周围环境换热量的一半以上是以辐射方式进行的。

本文以上述方程为基础，结合实测数据，应用CFD软件，对地板辐射、顶板辐射以及两者有机结合方式条件下，飞机维修车间内的温度场进行分析。

3.1. 地板辐射供暖结果分析

从机库流场的温度分布图(图2)可以看出，机库大厅内的温度分布较为均匀，即飞机检修和人体活动的区域温度分布较为均匀，并且温度云图呈对称形状，这是由于机库的热源是温度较为均匀的辐射地板。

由于密度差引起的大空间自然对流影响，机库内整体的垂直温度分布情况是地面附近温度由于靠近辐射板，温度比较高，这个高度正好人体脚部的位置，人体脚部对于环境温度最为敏感，而脚部温度恰好比较高，所以这种供暖方式对身处辐射环境中工作的人员健康较为有利的 [8] 。但随着高度上升温度先是逐渐下降然后慢慢上升，到达机库顶部时温度梯度变大。

由图3可见，机库内人员活动区域和维修部的流速较小，仅在各个壁面附近的边界层内有较高的速度分布，这是由于地板辐射采暖的主要换热方式是辐射传热，经研究得知，机库辐射供暖的辐射换热量占总供热量的65%左右，对流换热量占总换热量的30%左右 [9] 。

上升的气流到达机库顶部时，开始向下运动，与继续上升的气流形成小规模的涡流，使得维修车间内气流形式为两个大循环。由于靠近机库大门的位置，温度较低，流速较快。因此，为减小气流下降速度，提高机库内温度场的均匀性，应该对机库大门采取一定的措施，减少冷风的渗透，这样也有利于节能 [10] 。

3.2. 顶板辐射供暖结果分析

从图3还可看出，机库大厅内的温度有明显的分层现象。整个机库大厅内的温度分布与地板辐射供暖温度场相反，人体活动和飞机检修的高度温度较低，随着高度的增加温度逐渐上升，顶棚的温度最高。地面和墙壁由于他们自身的温度较低，和它们对流换热后，温度也降低。人体活动的区域温度梯度为1℃/m。机库内飞机由于体型庞大，在喷漆、退漆过程中需要分布一致的温度环境，因此这种形式的供暖方式是不适合的。

从截面X = 40和X = 80的云图可以看出，由于受到顶棚辐射板加热的效果，顶棚处的空气温度较高，空气密度比其它地方要小，这些密度较低的气体沿着顶板向四周扩散，由于热量集中于机库较高的部位，属于无效热量，造成了大量的能源浪费。由截面Z = 40的云图可以看出，在机库大厅的外墙温度较低，这个低温区域沿着外墙延伸到顶板。这个区域的形状同截面的外形比较吻合，说明模拟的结果与实际结果较为接近。

3.3. 两种辐射采暖结果分析

在X-Z面的俯视图划分如下的网格，取出五个测点进行分析室内的温度场，这样取的几个测点比较有代表性，如图4所示。

由图5得到不同的两种温度场特性，地板辐射采暖的机库其温度场特点是靠近地面辐射板的边界层温度较高，然后随着高度的增加逐渐降低，机库整体温度是下部温度较高，辐射地板的表面与其周围的换热温差较小，机库上部的温差较大，这是由于大空间对流导致的。而顶板辐射采暖的温度场特性是上部温度较高，从上到下呈倒温度梯度，温度分层现象比较明显 [11] ，机库下部温度较低，人员工作和飞机周围温度整体比地板辐射采暖的情况下要低，这是由于它和地板辐射采暖的热源位置不同引起的，由于顶板的温度较高导致它与机库内的空气换热温差大，温度场分布不均匀。

经过上面的分析，可以看出，在飞机维修车间应该采用地板辐射采暖与顶板辐射采暖相结合的办法，这样就会使整个机库的温度场分布均匀，使喷漆、退漆在较为稳定的温度环境中进行。

4. 结论

1) 使用地板辐射供暖的机库室内温度分布比较均匀，工作区的温度梯度较小，为1℃/m，机库顶端由于大空间对流换热，温度梯度较大。

2) 采用顶板辐射采暖的机库内，温度分层现象明显，顶部温度梯度较大，工作区整体温度较地板辐

射采暖要低，基本满足要求。但是其温度较高的区域集中在机库的顶端，造成大量能源的浪费，在机库等高大空间内使用不合适。

3) 应该采用地板辐射供暖与顶板辐射供暖相互结合的方式，这样就会使整个机库的温度场分布均匀。尽管机库内的飞机体型庞大，但也能使喷漆、退漆的过程都处在较为稳定的温度环境中，可有效消除辐射的不对称性，从而不会导致飞机机身的色彩变化。

参考文献