工程设计学报, 2024, 31(3): 357-367 doi: 10.3785/j.issn.1006-754X.2024.03.205

可靠性与保质设计

送风参数对老化箱内部温湿度场的影响研究

孙敏,,, 卢丰源, 赵宇轩, 王青春, 陈忠加,

北京林业大学 工学院,北京 100083

Study on influence of air supply parameters on temperature and humidity field inside aging chamber

SUN Min,,, LU Fengyuan, ZHAO Yuxuan, WANG Qingchun, CHEN Zhongjia,

School of Technology, Beijing Forestry University, Beijing 100083, China

通讯作者: 陈忠加(1981—),男,内蒙古乌海人,副教授,硕士生导师,博士,从事农业装备工程与技术研究,E-mail: chenzhongjia@bjfu.edu.cn

收稿日期: 2023-10-16   修回日期: 2023-12-14  

基金资助: 国家自然科学基金资助项目.  51475255

Received: 2023-10-16   Revised: 2023-12-14  

作者简介 About authors

孙 敏(1998—),女,湖南娄底人,硕士生,从事热力学仿真研究,E-mail:1061758591@qq.com,https://orcid.org/0009-0008-3094-6995 , E-mail:1061758591@qq.com

摘要

利用老化箱对机械电子元件进行老化筛选,能有效提高产品的可靠性。老化箱内部温湿度场的均匀性决定了其整体性能,这对老化箱功能有至关重要的影响。为了提高老化箱的性能,基于流体力学原理,利用CFD(computational fluid dynamics,计算流体力学)仿真软件对老化箱简化模型进行仿真模拟,并根据仿真结果确定老化箱内部速度场与温度场的分布情况,以优选不同送风温度、送风速度和送风角度组合下的送风方案。仿真试验采用正交试验法,以送风温度、送风速度和送风角度为试验因素,以能量利用系数、温度不均匀系数和相对湿度不均匀系数为评价指标。通过对仿真试验结果进行极差和方差分析可知,送风角度对能量利用系数的影响最显著,送风温度对温度不均匀系数和相对湿度不均匀系数的影响最显著。针对3个评价指标,得到3种最优方案:送风温度为90 ℃、送风速度为10 m/s、送风角度为0°,送风温度为90 ℃、送风速度为12 m/s、送风角度为-10°,以及送风温度为90 ℃、送风速度为8 m/s、送风角度为0°。最后,开展了老化箱内部温湿度测量实验,并与仿真结果进行了对比。结果表明,仿真结果与实验结果的相对误差较小,均在合理范围内,验证了仿真试验的可靠性和有效性。研究以提高老化箱性能和能量利用率为导向,可为同类老化箱的设计及其送风参数的设置提供参考。

关键词: 计算流体力学 ; 老化箱 ; 送风参数 ; 正交试验 ; 不均匀系数

Abstract

Using aging chamber to screen mechanical and electronic components can effectively improve the reliability of products. The uniformity of the temperature and humidity field inside the aging chamber determines its overall performance, which has a crucial impact on the functionality of the aging chamber. With the goal of enhancing the performance of the aging chamber, the simulation was conducted on the simplified model of aging chamber by using CFD (computational fluid dynamics) software based on the fluid dynamics principle. According to the simulation results, the distribution of velocity field and temperature field inside the aging chamber was determined to optimize the air supply scheme under the combination of different air supply temperature, air supply speed and air supply angle. The orthogonal test method was adopted in the simulation test, with air supply temperature, air supply speed and air supply angle as test factors, and energy utilization coefficient, temperature non-uniformity coefficient and relative humidity non-uniformity coefficient as evaluation indexes. Through range and variance analysis of simulation results, it could be seen that the influence of air supply angle on energy utilization coefficient was the most significant, and the influence of air supply temperature on temperature non-uniformity coefficient and relative humidity non-uniformity coefficient was the most significant. For the three evaluation indexes, three optimal schemes were obtained: air supply temperature of 90 ℃, air supply speed of 10 m/s, air supply angle of 0°, air supply temperature of 90 ℃, air supply speed of 12 m/s, air supply angle of -10°, and air supply temperature of 90 ℃, air supply speed of 8 m/s, air supply angle of 0°. Finally, the temperature and humidity measurement experiment inside the aging chamber was carried out and compared with the simulation results. The results showed that the relative errors between the simulation results and the experimental results were small and within a reasonable range, which verified the reliability and effectiveness of the simulation test. The research is oriented to improve the performance and energy utilization of aging chamber, which can provide reference for the design of similar aging chambers and the setting of air supply parameters.

Keywords: computational fluid dynamics ; aging chamber ; air supply parameter ; orthogonal test ; non-uniformity coefficient

PDF (3199KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

孙敏, 卢丰源, 赵宇轩, 王青春, 陈忠加. 送风参数对老化箱内部温湿度场的影响研究[J]. 工程设计学报, 2024, 31(3): 357-367 doi:10.3785/j.issn.1006-754X.2024.03.205

SUN Min, LU Fengyuan, ZHAO Yuxuan, WANG Qingchun, CHEN Zhongjia. Study on influence of air supply parameters on temperature and humidity field inside aging chamber[J]. Chinese Journal of Engineering Design, 2024, 31(3): 357-367 doi:10.3785/j.issn.1006-754X.2024.03.205

随着机械设备的应用场景不断拓展,机械设备电子元件的老化问题逐渐凸显。通过老化箱筛选出性能稳定的电子元件是目前最常见的一种老化筛选方法。开展老化筛选的重要前提是保证老化箱内部流场的均匀性。CFD(computational fluid dynamics,计算流体力学)仿真模拟是一种高效的数值模拟方法[1],可有效探究老化箱内部温湿度场的分布情况。

近年来,国内外学者针对烘房、发电机房等设施的送风参数优化开展了大量研究。王振文等[2]对热泵烘房的送风工艺参数进行了优化,结果表明,排风速度为6 m/s时热泵烘房内的气流分布较均匀。曹如玥[3]对蜂箱调温过程进行了数值模拟,结果显示:送风速度对蜂箱流体域内湍流强度的影响不显著;通风孔径对湍流强度的影响最为显著,其次是通风位置。林海等[4]综合考虑了多项工艺参数,通过仿真优化得到:当初始温度为25 ℃时,送风温度为207 ℃、送风速度为5 m/s的参数组合下汽车车身烘房的能耗最低。魏亚兴[5]选取了10个对矿井通风巷热交换产生影响的因素,以温度作为评价指标开展了正交试验,结果显示,进风温度对巷道终点温度的影响较大,且两者呈正相关。杜林昕等[6]为探究进风量、进风角度和进风口位置对烤房内部气流分布情况的影响,设计了相应的正交试验,结果表明,上述3个送风参数改变均会对烤房内部的气体速度产生显著影响,影响的主次顺序为进风口位置、进风角度和进风量。Li等[7]选择送风口高度、送风口宽度、送风速度和热源强度作为影响因素,通过开展正交试验来优化发电机房的通风性能,经极差分析和方差分析可知,送风口高度对发电机房内部的流场和通风性能有显著影响。

综上所述,送风温度、送风口位置、送风口高度对烘箱等热循环设备内部流场的影响显著。基于此,本文将探究送风温度、送风速度和送风角度这3个参数对老化箱内部温湿度场的影响。通过设计9组正交试验来系统性地探讨不同送风温度、送风速度和送风角度下老化箱内部温度场、湿度场的分布情况,旨在为老化箱工艺参数的优化设计提供理论依据。

1 老化箱内部流场计算理论

1.1 湍流模型

一般情况下,对于老化箱内部的流体流动,通常认为雷诺数Re<2 300时为层流,Re>2 300时为湍流。雷诺数Re的计算式如下:

Re=VρLμ=VDrυ

式中:V为流道内流体的平均速度,m/s;ρ为流体密度,kg/m3L为特征长度,m;μ为动力黏度,Pa·s;Dr为流道特征尺寸,m;υ为运动黏度,m2/s。

老化箱的送风口为矩形,对于非圆形管道内的流动,L=Dr。根据本文所设置的最大风速12 m/s,计算得到老化箱内部的流体流动为湍流。在Fluent软件中调用湍流模型,将已知的最大风速代入马赫数Ma的计算式,得到Ma=0.035。Ma0.3,则老化箱内部的流体可视为不可压缩流体[8]

1.2 基本控制方程

在流体力学领域,流体流动和传热传质均遵循三大基本守恒定律,即质量守恒、动量守恒和能量守恒定律。

1)质量守恒方程。

流体流动在时间和空间尺度上的求解方法基于质量守恒方程。质量守恒是指流入和流出某个微元体的流体质量相等,可表示为:

ρt+ρux+ρvy+ρwz=SM

式中:t为时间,s;u、v、wxyz方向上流体的速度分量,m/s;SM 为质量源项,kg/(m3·s)。由于老化箱内会发生水的相变,因此引入水蒸气质量源项,表示水由气相转化为液相所产生的质量变化。

通过前文分析可知,老化箱内部的流体为不可压缩流体,即流体密度为常数,因此质量守恒方程可简化为:

ρux+ρvy+ρwz=SM

2)动量守恒方程。

动量守恒是指一个微元体的始末动量变化量等于该运动过程中微元体所受力之和。Naiver-Stokes方程(N-S方程)是描述黏性不可压缩流体动量守恒的运动方程。本文老化箱内部的不可压缩流体在xyz方向上的动量方程[9]可简化为:

ρut+uux+vuy+wuz=ρFx-px+μ2ux2+2uy2+2uz2ρvt+uvx+vvy+wvz=ρFy-py+μ2vx2+2vy2+2vz2ρwt+uwx+vwy+wwz=ρFz-pz+μ2wx2+2wy2+2wz2

式中:p为微元体上的压力,N;FxFyFz为作用在微元体上的xyz方向分力,N。

3)能量守恒方程。

能量守恒是指流动过程中流体的总能量保持不变。能量守恒方程可表示为:

ρTt+ρuTx+ρvTy+ρwTz=xkcpTx+ykcpTy+zkcpTz+ST

式中:cp 为比热容,kJ/(kg·℃);T为温度,℃;k为流体传热系数,W/(m2·℃);ST 为黏性耗散项[10],W。

1.3 多相流模型

在计算流体力学中,多相流模型通常将流场中的各相介质看作连续介质,即相与相之间的流动是不间断的。本文选用流体体积(volume of fluid, VOF)模型来描述老化箱内部水蒸气的冷凝现象。VOF模型主要用于流体射流、流体中气泡运动、气液界面的稳态和瞬态处理等,其通过求解一组动量方程并跟踪整个计算域内每种流体的容积比来模拟多种无法互相融合的流体[11-12]

2 老化箱送风参数的正交试验方案

为了研究送风参数对老化箱内部温湿度场的影响,选择送风温度、送风速度、送风角度作为试验因素(分别记为因素A、因素B、因素C),每个因素选取3个水平(见表1),开展三因素三水平正交试验。在送风相对湿度为80%的条件下,采用SPSS软件中的L9(33)正交表设计正交试验方案,如表2所示。

表1   老化箱送风参数的正交试验因素水平表

Table 1  Factor and level table for orthogonal test of air supply parameters of aging chamber

水平因素

A(送风温度)/

B(送风速度)/(m/s)C(送风角度)/(º)
1508-10
270100
3901210

新窗口打开| 下载CSV


表2   老化箱送风参数的正交试验方案

Table 2  Orthogonal test scheme of air supply parameters of aging chamber

序号因素方案
A/℃B/(m/s)C/(º)
1508-10A1B1C1
250100A1B2C2
3501210A1B3C3
47080A2B1C2
5701010A2B2C3
67012-10A2B3C1
790810A3B1C3
89010-10A3B2C1
990120A3B3C2

新窗口打开| 下载CSV


3 老化箱内部流场仿真模型构建

3.1 老化箱三维模型

老化箱内部置有循环风扇、加湿器和加热管,其工作原理如图1所示。加湿器释放的水蒸气与经加热管加热的空气混合后形成湿热气体;在循环风扇的作用下,湿热气体由送风口进入工作室,经热量传递后,气体由回风口离开工作室,从而实现热风循环。在气体循环的过程中,高温和高湿的环境使得电子元件加速老化[13-14]

图1

图1   老化箱工作原理

Fig.1   Working principle of aging chamber


本文老化箱采用侧送侧回热风循环方式。与其他送风方式相比,侧送侧回送风方式更能保证老化箱内气体的速度、温度和湿度不均匀系数满足要求[15-16]。以老化箱内部空间为研究对象,采用SolidWork软件构建其1∶1三维模型,忽略工作室内对气体流动影响很小的紧固件和结构件。所构建的老化箱三维简化模型如图2所示。其中:箱体内部长1 600 mm,宽1 600 mm,高1 600 mm;箱体上方设有4个送风口,箱体下方设有4个回风口,送风口尺寸为170 mm×115 mm,回风口尺寸为680 mm×165 mm;工作面距离箱底400~1 200 mm。

图2

图2   老化箱三维简化模型

Fig.2   Simplified 3D model of aging chamber


3.2 网格划分

本文的老化箱三维模型较为规则,以及采用结构化网格的计算方式耗时短且计算效率高,故采用结构化网格对老化箱三维模型进行划分。考虑到网格数量过多会增加计算时长,在保证仿真结果精度不受影响的前提下,选择合适尺寸的网格单元对老化箱三维模型进行划分,最终的网格数量为160万个,如图3所示。

图3

图3   老化箱三维模型网格划分

Fig.3   Grid division of aging chamber 3D model


3.3 边界条件设置

鉴于本文所研究的老化箱内部的气流为不可压缩流体[17-18]且已知送风口尺寸和气流速度,将进口边界条件设为速度进口(velocity-inlet),出口边界条件设为压力出口(pressure-outlet)。因回风口与大气连通,将出口表面的压力设为0 Pa。

此外,虽然在设计老化箱箱体壁面时加入了保温材料,但是箱内的热量仍会通过壁面传递到外界,故将壁面的热边界条件设为固定温度壁面。

4 老化箱内部流场仿真结果分析

4.1 正交试验结果分析

根据表2所示的正交试验方案,利用CFD仿真软件对老化箱内部的速度场和温度场进行仿真模拟。不同送风方案下老化箱内部流场区域距底面600 mm处水平面的流线图和温度云图分别如图4图5所示。其中:图5中温度量程的上下限取自老化箱内温度的最高值和最低值,以便更好地观察箱内温度分布的均匀性。

图4

图4   老化箱内距底面600 mm处水平面的流线图

Fig.4   Streamline diagram of horizontal plane at a distance of 600 mm from bottom inside aging chamber


图5

图5   老化箱内距底面600 mm处水平面的温度云图

Fig.5   Temperature cloud map of horizontal plane at a distance of 600 mm from bottom inside aging chamber


图4可以看出,在9组送风方案下,老化箱内部的气流运动趋势基本一致,即气流从送风口进入老化箱后向下移动,在接触到壁面后往上散开,总体上气流运动大致呈上下来回移动。通过对比发现,在A3B1C3、A3B2C1、A3B3C2这3组方案下,老化箱内距底面600 mm处水平面上的气流分布相对较均匀。

图5可以看出,在9组送风方案下,老化箱内部的温度分布大致相同,箱内中间温度高于壁面温度。通过对比方案A1B1C1/A1B2C2/A1B3C3以及方案A1B1C1/A2B1C2/A3B1C3可知,当送风角度逐渐增大时,老化箱内低温区的面积逐渐减小,高温区集中在中心区域。此外,从图5中还可以看出,不同送风速度下老化箱内部的温度场基本一致,这与文献[19]中分析得到的送风速度对新型单侧供热烘干房的内部温度场的影响结果吻合。由此说明,单侧送风方式下送风速度对老化箱内部温度分布的影响较小。

对不同送风方案下老化箱内部流场的仿真结果进行处理,得到对应的老化箱性能评价指标,结果如表3所示。

表3   不同送风方案下老化箱的性能评价指标

Table 3  Performance evaluation indexes of aging chamber under different air supply schemes

序号方案能量利用系数温度不均匀系数相对湿度不均匀系数
1A1B1C10.909 60.077 10.125 2
2A1B2C21.058 60.078 80.119 7
3A1B3C31.045 00.080 20.127 5
4A2B1C21.028 10.054 30.104 8
5A2B2C31.029 10.054 10.106 4
6A2B3C10.938 80.050 90.106 8
7A3B1C31.061 90.040 20.092 0
8A3B2C10.983 40.038 80.098 7
9A3B3C21.051 30.039 50.099 0

新窗口打开| 下载CSV


4.2 正交试验结果极差分析

利用极差分析法来分析送风温度、送风速度和送风角度这3个因素对老化箱的能量利用系数、温度不均匀系数和相对湿度不均匀系数的影响程度。令Ki 表示各因素各水平下各评价指标值之和,ki 表示各因素各水平下各评价指标值的平均值,R表示ki 值的极差。R值越大,表示该因素对相应评价指标的影响越大。各因素对老化箱能量利用系数影响的极差分析结果如表4所示。为更加直观地观察老化箱能量利用系数随各因素各水平的变化趋势,绘制相应的位级图,如图6所示。

表4   各因素对老化箱能量利用系数影响的极差分析结果

Table 4  Range analysis results of influence of each factor on energy utilization coefficient of aging chamber

分析指标因素A因素B因素C
K13.013 22.999 62.831 8
K22.996 03.071 13.138 0
K33.096 63.035 13.136 0
k11.004 40.999 90.943 9
k20.998 71.023 71.046 0
k31.032 21.011 71.045 3
R0.033 50.023 80.102 1
优水平A3B2C2

新窗口打开| 下载CSV


图6

图6   老化箱能量利用系数的影响因素位级图

Fig.6   Level diagram of influencing factors of energy utilization coefficient of aging chamber


表4可知,因素C对应的极差最大,即送风角度是影响老化箱能量利用系数的主要因素;不同因素对能量利用系数的影响顺序为C>A>B。从图6中可直观看出,因素C(送风角度)的水平变化对老化箱能量利用系数的影响最大,而因素A(送风温度)和因素B(送风速度)对能量利用系数的影响相对不明显。通过观察发现,当送风角度从负角度变为正角度后,老化箱的能量利用系数显著提高。根据图6,在侧重考虑老化箱的能量利用效率时,送风参数的最优组合为A3B2C2,即送风温度为90 ℃、送风速度为10 m/s、送风角度为0º,在该组合下老化箱运行过程中的能量利用效率较高。

各因素对老化箱温度不均匀系数影响的极差分析结果如表5所示。同样,为更加直观地观察老化箱温度不均匀系数随各因素各水平的变化趋势,绘制相应的位级图,如图7所示。

表5   各因素对老化箱温度不均匀系数影响的极差分析结果

Table 5  Range analysis results of influence of each factor on temperature non-uniformity coefficient of aging chamber

分析指标因素A因素B因素C
K10.236 10.171 60.166 8
K20.159 30.171 70.172 6
K30.118 50.170 60.174 5
k10.078 70.057 20.055 6
k20.053 10.057 20.057 5
k30.039 50.056 90.058 2
R0.039 20.000 30.002 6
优水平A3B3C1

新窗口打开| 下载CSV


图7

图7   老化箱温度不均匀系数的影响因素位级图

Fig.7   Level diagram of influencing factors of temperature non-uniformity coefficient of aging chamber


表5可知,因素A对应的极差最大,即送风温度是影响老化箱温度不均匀系数的主要因素;不同因素对温度不均匀系数的影响顺序为A>C>B。观察图7可知,老化箱的温度不均匀系数随因素A呈线性变化,即送风温度越高,箱内温度分布的均匀性越好,且送风温度对温度不均匀系数的影响最为显著;因素B和因素C对温度不均匀系数的影响较小,即温度不均匀系数基本不随送风速度和送风角度的水平变化而变化。综上,在侧重考虑老化箱的温度分布均匀性时,送风参数的最优组合为A3B3C1,即送风温度为90 ℃、送风速度为12 m/s、送风角度为-10°,该组合下老化箱内部温度分布均匀,有利于电子元件的老化筛选。

各因素对老化箱相对湿度不均匀系数影响的极差分析结果如表6所示。同样,为更加直观地观察老化箱相对湿度不均匀系数随各因素各水平的变化趋势,绘制相应的位级图,如图8所示。

表6   各因素对老化箱相对湿度不均匀系数影响的极差分析结果

Table 6  Range analysis results of influence of each factor on relative humidity non-uniformity coefficient of aging chamber

分析指标因素A因素B因素C
K10.372 40.322 00.330 7
K20.318 00.324 80.323 5
K30.289 70.333 30.325 9
k10.124 10.107 30.110 2
k20.106 00.108 30.107 8
k30.096 60.111 10.108 6
R0.027 50.003 80.002 4
优水平A3B1C2

新窗口打开| 下载CSV


图8

图8   老化箱相对湿度不均匀系数的影响因素位级图

Fig.8   Level diagram of influencing factors of relative humidity non-uniformity coefficient of aging chamber


表6可知,因素A对应的极差最大,即送风温度是影响老化箱相对湿度不均匀系数的主要因素;不同因素对相对湿度不均匀系数的影响顺序为A>B>C。观察图8可知,相对湿度不均匀系数随因素A呈线性变化,即送风温度越高,相对湿度分布的均匀性越好。在侧重考虑老化箱的湿度分布均匀性时,送风参数的最优组合为A3B1C2,即送风温度为90 ℃、送风速度为8 m/s、送风角度为0°,该组合下老化箱内的相对湿度不均匀系数较小,即相对湿度分布均匀,有利于电子产品的老化筛选。

综合以上分析可得:送风角度为影响老化箱能量利用系数的主要因素;送风温度对温度不均匀系数的影响最为显著,且对相对湿度不均匀系数的影响也十分显著,均呈线性相关性,即送风温度的提高能有效提高箱内温度和相对湿度分布的均匀性;送风速度对3个评价指标的影响均不明显。

4.3 正交试验结果方差分析

为进一步分析送风温度、送风速度和送风角度这3个因素对老化箱能量利用系数、温度不均匀系数和相对湿度不均匀系数的影响程度,利用方差分析法对正交试验结果进行处理。定义显著度P,当P<0.01时,表明该因素对评价指标的影响极为显著,用“**”表示;当P<0.05时,表明该因素对评价指标的影响显著,用“*”表示;当P>0.05时,表明该因素对评价指标无显著影响,用“—”表示。同时,引入因素贡献率来量化各因素对不同评价指标的影响程度,其定义为某一因素的离差平方和与总离差平方和之比。F值为正态分布表中的边界值。

各因素对老化箱能量利用系数影响的方差分析结果如表7所示。表中:空白表示不存在此项分析结果,下文同。为更加直观地分析各因素对老化箱能量利用系数的贡献率,绘制柱状图,如图9所示。

表7   各因素对老化箱能量利用系数影响的方差分析结果

Table 7  Variance analysis results of influence of each factor on energy utilization coefficient of aging chamber

方差来源离差平方和自由度均方P显著性
总和0.024 508 10380.012 254 051
因素A0.001 930 19620.000 965 0981.881 5720.347 033
因素B0.000 852 05620.000 426 0280.830 5910.546 272
因素C0.020 700 00920.010 350 00420.178 5500.047 218*
误差0.001 025 84220.000 512 921

新窗口打开| 下载CSV


图9

图9   各因素对老化箱能量利用系数的贡献率

Fig.9   Contribution rate of each factor to energy utilization coefficient of aging chamber


表7可知,因素C(送风角度)的变化对老化箱能量利用系数有显著影响,因素A(送风温度)和因素B(送风速度)对能量利用系数的影响不显著。由图9可知,因素C的贡献率为84.46%,对能量利用系数的影响最大,因素A和因素B的贡献率均低于8%,即这2个因素对能量利用系数的影响较小。综上,方差分析结果与极差分析结果相吻合,均表明送风角度对老化箱能量利用系数的影响最大。

各因素对老化箱温度不均匀系数影响的方差分析结果如表8所示。同样,为更直观地分析各因素对老化箱温度不均匀系数的贡献率,绘制柱状图,如图10所示。

表8   各因素对老化箱温度不均匀系数影响的方差分析结果

Table 8  Variance analysis results of influence of each factor on temperature non-uniformity coefficient of aging chamber

方差来源离差平方和自由度均方P显著性
总和0.002 390 04180.001 195 019
因素A0.002 376 96020.001 188 4801 128.303 7970.000 886**
因素B0.000 000 24720.000 000 1230.117 0890.895 184
因素C0.000 010 72720.000 005 3635.091 7720.164 156
误差0.000 002 10720.000 001 053

新窗口打开| 下载CSV


图10

图10   各因素对老化箱温度不均匀系数的贡献率

Fig.10   Contribution rate of each factor to temperature non-uniformity coefficient of aging chamber


表8可知,因素A(送风温度)对老化箱温度不均匀系数有显著影响;因素B(送风速度)和因素C(送风角度)对温度不均匀系数的影响很小,与极差分析结果相吻合。由图10可知,因素B和因素C对温度不均匀系数的贡献率接近于0,因此可忽略这2个因素对温度不均匀系数的影响;因素A对温度不均匀系数的贡献率为99.45%,表明该因素对温度不均匀系数有着决定性作用。

各因素对老化箱相对湿度不均匀系数影响的方差分析结果如表9所示。同样,为更加直观地分析各因素对老化箱相对湿度不均匀系数的贡献率,绘制柱状图,如图11所示。

表9   各因素对老化箱相对湿度不均匀系数影响的方差分析结果

Table 9  Variance analysis results of influence of each factor on relative humidity non-uniformity coefficient of aging chamber

方差来源离差平方和自由度均方P显著性
总和0.001 243 42180.000 621 709
因素A0.001 177 72720.000 588 86335.002 7740.027 776*
因素B0.000 023 08720.000 011 5430.686 1500.593 067
因素C0.000 008 96020.000 004 4800.266 2970.789 704
误差0.000 033 64720.000 016 823

新窗口打开| 下载CSV


图11

图11   各因素对老化箱相对湿度不均匀系数的贡献率

Fig.11   Contribution rate of each factor to relative humidity non-uniformity coefficient of aging chamber


表9可知,因素A对应的P<0.05,说明送风温度对老化箱相对湿度不均匀系数有显著影响;从总体上看,3个因素对相对湿度不均匀系数的影响程度的次序为A>B>C,即送风温度与老化箱内部相对湿度的均匀性有很大关联性。由图11可知,因素A对相对湿度不均匀系数的贡献率高达94.72%,即送风温度对相对湿度不均匀系数有决定性影响;因素B、因素C和误差的贡献率均低于5%,说明这3个因素对相对湿度均匀性的影响很小。

5 实验验证

为验证上述仿真结果的可靠性与有效性,对老化箱内部的温度和湿度进行测量,并对比实验结果与仿真结果。

5.1 实验方案

所采用的温湿度检测系统主要由PT100温度传感器、HM100湿度传感器、Keysight34972A数据采集器和监控电脑组成。根据GB/T 10586—2006《湿热试验箱技术条件》[20],在老化箱内部布点并测量相应点位的温度和湿度,如图12所示。

图12

图12   老化箱实物图及箱内测点布置示意

Fig.12   Physical drawing of aging chamber and layout of measuring points in chamber


1)为保证测量分布全面,分别在老化箱工作室内进行3个水平面的布点测量。第1个水平面设置在距最底层样品架上方10 mm处,第2个水平面设置在工作室内的几何中心处,第3个水平面设置在距离顶面1/10工作室高度的位置处。

2)根据GB/T 10586—2006中的要求,共布置13个温度测点(I,II,…,XIII)和4个相对湿度测点(甲、乙、丙、丁),测点分布在3个平面上,中心测点位于工作室几何中心处,其余测点位于距离工作室壁面1/10边长处。

温度和湿度的测量步骤如下:实验开始时送风口持续输送温度为80 ℃和90 ℃、相对湿度为80%和90%的气流,使老化箱内部的温度和相对湿度达到规定值并稳定30 min;安装在相应测点处的温湿度传感器每间隔1 min测量一组温湿度数据并存储在数据采集器中,共测量30组数据。

5.2 实验结果

当送风温度分别为80 ℃和90 ℃时,老化箱内各测点温度的实测值与仿真值的对比如表10所示。由表10可知,老化箱内温度实测值与仿真值一致,两者之间的最大误差为2.3 ℃,误差在合理范围内,由此验证了仿真结果的可靠性与有效性。当送风温度为80 ℃时,上层4个测点的最大温度差为1.7 ℃,中层5个测点的最大温度差为1.5 ℃,下层4个测点的最大温度差为0.9 ℃;当送风温度为90 ℃时,上层4个测点的最大温度差为1.4 ℃,中层5个测点的最大温度差为1.1 ℃,下层4个测点的最大温度差为0.4 ℃。通过计算可得,当送风温度分别为80 ℃和90 ℃时,老化箱内的平均温度分别为78.0 ℃和87.7 ℃,箱内温度与实际送风温度有一定差距,这是因为老化箱存在热耗散现象,热量通过壁面或未密封完全的空隙向外界传递热量。

表10   老化箱内部温度的实测值与仿真值对比 (℃)

Table 10  Comparison of measured and simulated temperature inside aging chamber

测点送风温度为80 ℃送风温度为90 ℃
实测值仿真值实测值仿真值
I78.676.786.986.5
II77.377.887.087.0
III76.976.188.387.1
IV78.578.687.986.7
V77.278.188.185.8
VI78.777.887.186.7
VII77.678.288.286.8
VIII77.878.087.787.2
IX78.177.587.285.8
X78.876.788.186.5
XI78.478.987.887.0
XII77.978.788.288.1
XIII78.577.887.986.7

新窗口打开| 下载CSV


当相对湿度分别为80%和90%时,老化箱内各测点相对湿度的实测值与仿真值的对比如表11所示。由表11可知,实际测量得到的相对湿度与仿真得到的相对湿度基本一致,两者之间的误差不超过1.3%,验证了仿真结果的可靠性与有效性。

表11   老化箱内部相对湿度的实测值与仿真值对比 (%)

Table 11  Comparison of measured and simulated relative humidity inside aging chamber

测点相对湿度为80%相对湿度为90%
实测值仿真值实测值仿真值
80.381.689.989.4
80.580.590.090.0
80.081.090.090.0
80.281.290.190.4

新窗口打开| 下载CSV


6 结 论

本文选择送风温度、送风速度、送风角度这3个因素设计了三因素三水平正交试验,研究了各个因素对老化箱内部温度均匀性、相对湿度均匀性和速度均匀性的影响,相较于单因素研究更具综合性和实际意义。

通过对正交试验结果进行极差分析和方差分析可知,对于能量利用系数,各因素的影响次序为送风角度、送风温度、送风速度,即送风角度的影响最为显著。对于温度均匀性,各因素的影响次序为送风温度、送风角度、送风速度,即送风温度的影响最为显著,其对温度不均匀系数的贡献率高达99.45%。对于相对湿度不均匀系数,各因素的影响次序为送风温度、送风速度、送风角度,即送风温度的影响最为显著。综上,针对3个不同的评价指标,最优送风方案分别为:送风温度为90 ℃、送风速度为10 m/s、送风角度为0°,送风温度为90 ℃、送风速度为12 m/s、送风角度为-10°,以及送风温度为90 ℃、送风速度为8 m/s、送风角度为0°。

实验验证结果表明,老化箱内温度、湿度的仿真值与实测值之间的差异较小,其中温度的最大误差为2.3 ℃,相对湿度的最大误差为1.3%,误差均在合理范围内,由此验证了仿真结果的可靠性与有效性。

参考文献

马超群张凯柴麟.

机械式自动垂直钻具执行机构内部流场规律研究

[J].工程设计学报,2022293):384-393. doi:10.3785/j.issn.1006-754X.2022.00.044

DOI:10.3785/j.issn.1006-754X.2022.00.044      [本文引用: 1]

MA C QZHANG KCHAI Let al.

Research on internal flow field law of mechanical automatic vertical drilling tool actuator

[J]. Chinese Journal of Engineering Design, 2022293): 384-393.

DOI:10.3785/j.issn.1006-754X.2022.00.044      [本文引用: 1]

王振文吴敏徐新民.

热泵烘房结构及参数优化仿真设计

[J].农业机械学报,202051():464-475. doi:10.6041/j.issn.1000-1298.2020.S1.055

DOI:10.6041/j.issn.1000-1298.2020.S1.055      [本文引用: 1]

WANG Z WWU MXU X Met al.

Optimal simulation design of structure and parameter in heat pump drying room

[J]. Transactions of the Chinese Society for Agricultural Machinery, 202051(): 464-475.

DOI:10.6041/j.issn.1000-1298.2020.S1.055      [本文引用: 1]

曹如玥.

基于CFD的蜂箱调温过程数值模拟及实验研究

[D].昆明昆明理工大学2021.

[本文引用: 1]

CAO R Y.

Numerical simulation and experimental research of beehive temperature adjustment process based on CFD

[D]. KunmingKunming University of Science and Technology2021.

[本文引用: 1]

林海叶永伟陈建.

基于CFD的汽车车身烘房节能研究

[J].机械制造,20104810):71-73. doi:10.3969/j.issn.1000-4998.2010.10.025

DOI:10.3969/j.issn.1000-4998.2010.10.025      [本文引用: 1]

LIN HYE Y WCHEN Jet al.

Research on energy saving of oven for automobile body based on CFD

[J]. Machinery, 20104810): 71-73.

DOI:10.3969/j.issn.1000-4998.2010.10.025      [本文引用: 1]

魏亚兴.

矿井通风井巷热交换规律数值试验研究

[D].长沙中南大学2012.

[本文引用: 1]

WEI Y X.

Numerical test research on the law of mine ventilation roadway heat exchange

[D]. ChangshaCentral South University2012.

[本文引用: 1]

杜林昕黄亚宇查蕾蕾.

基于CFD的密集式烤房气流均匀性研究

[J].农业装备与车辆工程,20215912):126-129. doi:10.3969/j.issn.1673-3142.2021.12.027

DOI:10.3969/j.issn.1673-3142.2021.12.027      [本文引用: 1]

DU L XHUANG Y YZHA L Let al.

Research on air flow uniformity of intensive barn based on CFD

[J]. Agricultural Equipment & Vehicle Engineering, 20215912): 126-129.

DOI:10.3969/j.issn.1673-3142.2021.12.027      [本文引用: 1]

LI J XLI A GZHANG Cet al.

Analysis and optimization of air distribution and ventilation performance in a generator hall using an innovative attached air supply mode

[J]. Building and Environment, 2022216108993.

[本文引用: 1]

白志鹏.

可压缩流体的瓦斯抽采管路流阻计算方法研究

[J].煤,2016254):8-11. doi:10.3969/j.issn.1005-2798.2016.04.003

DOI:10.3969/j.issn.1005-2798.2016.04.003      [本文引用: 1]

BAI Z P.

Study on calculation method of gas drainage pipeline flow resistance of compressible fluid

[J]. Coal, 2016254): 8-11.

DOI:10.3969/j.issn.1005-2798.2016.04.003      [本文引用: 1]

陈忠加卢丰源雷雯雯.

基于强迫对流的热风干燥烘房送风速度及温度优选

[J].农业工程学报,202238():37-46. doi:10.11975/j.issn.1002-6819.2022.z.005

DOI:10.11975/j.issn.1002-6819.2022.z.005      [本文引用: 1]

CHEN Z JLU F YLEI W Wet al.

Air

velocity and temperature optimization of hot air drying room based on forced convection[J]. Transactions of the Chinese Society for Agricultural Machinery, 202238(Supplement 1): 37-46.

DOI:10.11975/j.issn.1002-6819.2022.z.005      [本文引用: 1]

任塞峰唐汝宁李军军.

化学分析实验室换气次数数值模拟分析

[J].建筑热能通风空调,2021407):70-74. doi:10.3969/j.issn.1003-0344.2021.07.018

DOI:10.3969/j.issn.1003-0344.2021.07.018      [本文引用: 1]

REN S FTANG R NLI J J.

Research on numerical simulation of air exchange rate in chemical analysis laboratory

[J]. Building Energy & Environment, 2021407): 70-74.

DOI:10.3969/j.issn.1003-0344.2021.07.018      [本文引用: 1]

周鹏.

高炉风口回旋区燃烧特性的数值模拟研究

[D]. 赣州江西理工大学2022. doi:10.1016/j.fuel.2021.122490

DOI:10.1016/j.fuel.2021.122490      [本文引用: 1]

ZHOU P.

Numerical study of combustion characteristics in the raceway of an ironmaking blast furnace

[D]. GanzhouJiangxi University of Science and Technology2022.

DOI:10.1016/j.fuel.2021.122490      [本文引用: 1]

李强孟利民周琪.

基于机械虹膜机构的新型摆动泵及其流场分析

[J].机电工程,20234010):1607-1615. doi:10.3969/j.issn.1001-4551.2023.10.017

DOI:10.3969/j.issn.1001-4551.2023.10.017      [本文引用: 1]

LI QMENG L MZHOU Qet al.

Flow characteristics of new swing pump based on mechanical iris mechanism

[J]. Journal of Mechanical & Electrical Engineering, 20234010): 1607-1615.

DOI:10.3969/j.issn.1001-4551.2023.10.017      [本文引用: 1]

吴昊龙铁明蓝文清.

高低温环境试验箱循环风机优化设计

[J].环境技术,2022406):184-190. doi:10.3969/j.issn.1004-7204.2022.06.041

DOI:10.3969/j.issn.1004-7204.2022.06.041      [本文引用: 1]

WU HLONG T MLAN W Q.

Optimization of blower in enviornmental test chamber

[J]. Environmental Technology, 2022406): 184-190.

DOI:10.3969/j.issn.1004-7204.2022.06.041      [本文引用: 1]

高祥马玫.

双轴型热空气交换试验箱的研制

[J].合成材料老化与应用,2022514):109-111. doi:10.3969/j.issn.1671-5381.2022.4.hccllhyyy202204035

DOI:10.3969/j.issn.1671-5381.2022.4.hccllhyyy202204035      [本文引用: 1]

GAO XMA M.

Development of biaxial hot air exchange test chamber

[J]. Synthetic Materials Aging and Application, 2022514): 109-111.

DOI:10.3969/j.issn.1671-5381.2022.4.hccllhyyy202204035      [本文引用: 1]

陈忠加雷雯雯王青春.

基于温度和速度均匀性的侧送风烘房设计及仿真

[J].农业工程学报,20213719):18-26. doi:10.11975/j.issn.1002-6819.2021.19.003

DOI:10.11975/j.issn.1002-6819.2021.19.003      [本文引用: 1]

CHEN Z JLEI W WWANG Q C.

Design and simulation of side air supply drying room based on temperature and velocity homogeneity

[J]. Transactions of the Chinese Society for Agricultural Machinery, 20213719): 18-26.

DOI:10.11975/j.issn.1002-6819.2021.19.003      [本文引用: 1]

方一鸣.

侧送风气流组织设计计算

[J].城市建筑,20138):268. doi:10.3969/j.issn.1673-0232.2013.08.238

DOI:10.3969/j.issn.1673-0232.2013.08.238      [本文引用: 1]

FANG Y M.

The design and calculation of side air supply airflow organization

[J]. Urbanism and Architecture, 20138): 268.

DOI:10.3969/j.issn.1673-0232.2013.08.238      [本文引用: 1]

GUI X LLUO X BWANG X Pet al.

Computational fluid dynamic (CFD) investigation of thermal uniformity in a thermal cycling based calibration chamber for MEMS

[J]. Heat and Mass Transfer, 20155112): 1705-1715.

[本文引用: 1]

张伟严平.

热风干燥烘房气流速度场均匀性优化

[J].轻工机械,2019371):83-87. doi:10.3969/j.issn.1005-2895.2019.01.016

DOI:10.3969/j.issn.1005-2895.2019.01.016      [本文引用: 1]

ZHANG WYAN P.

Optimization of airflow velocity field uniformity in hot air drying kiln

[J]. Light Industry Machinery, 2019371): 83-87.

DOI:10.3969/j.issn.1005-2895.2019.01.016      [本文引用: 1]

李会知侯孟言邢金超.

单侧供热汽车烘干房气流组织数值模拟研究

[J].山东农业大学学报(自然科学版),2021521):84-90. doi:10.3969/j.issn.1000-2324.2021.01.015

DOI:10.3969/j.issn.1000-2324.2021.01.015      [本文引用: 1]

LI H ZHOU M YXING J C.

Study on numerical simulation for airflow distribution in automotive baking room with a single side heating

[J]. Journal of Shandong Agricultural University (Natural Science Edition), 2021521): 84-90.

DOI:10.3969/j.issn.1000-2324.2021.01.015      [本文引用: 1]

中国工业机械联合会. 湿热试验箱技术条件: [S].北京中国标准出版社20061-12.

[本文引用: 1]

China Machinery Industry Federation. Specifications for damp heat testing chambers: [S]. BeijingStandards Press of China20061-12.

[本文引用: 1]

/