2. 吉林大学 汽车工程学院,吉林 长春 130025;
3. 哈尔滨固泰电子有限责任公司,黑龙江 哈尔滨 150060;
4. 吉林建筑大学 市政与环境工程学院,吉林 长春 130118
2. College of Automotive Engineering, Jilin University, Changchun 130025, China;
3. Harbin Gutai Electronic Co. Ltd., Harbin 150060, China;
4. College of Municipal and Environmental Engineering, Jilin University of Construction, Changchun 130118, China
随着社会的进步,我国汽车工业的飞速发展,汽车已经成为了日常生活中必不可少的一部分,人们对汽车的性能的要求也越来越高,在关注动力性、经济性的同时,汽车乘员热舒适性,也越来越引起人们的重视,它直接影响驾乘人员的身心健康,关系着行驶安全[1].热舒适感觉是人体对热环境的生理与心理因素的综合反映[2].根据美国供暖制冷空调工程师学会的定义[3],人体的热舒适性是由周围热环境等客观因素和人体自身机能等主观因素[4]共同决定的.
重型车辆主要用于执行特大任务工作,我国的重型车辆在设计初期受前苏联的影响较大,主要强调强劲的动力性,高度的可靠性,却忽视了乘员舱内乘员的舒适性.本文研究的重型车辆其外部车体由钢板构成,材料的导热系数较大,乘员舱内空间狭小,外界环境的变化对乘员舱影响较大.寒冷的冬季,在我国东北地区乘员舱的温度将低至-30 ℃,发动机工作产生的热量、乘员工作产生的热量远远不足以保证乘员的热舒适性,将影响到驾驶员的反应操纵能力,因此重型车辆中安装采暖系统刻不容缓.
目前国内外在重型车乘员舱内热环境研究少,特别是寒区冬季重型车乘员舱内乘坐热环境改进优化方面几乎处于空白,而在冬季乘员舱内温度很低,应在乘员舱内布置合理的采暖设备,从而改善乘员的热舒适性.本文从以下几方面开展工作:
(1) 设计一种采暖系统,利用燃油预热器对发动机的冷却水进行预热,将高温冷却水送入到采暖系统的换热器内,为了保证采暖系统的热风充分的被利用,设计了一种两通道的风道结构形式.
(2) 对某重型车乘员舱内部进行结构简化,进行三维建模,对乘员舱的采暖系统进行结构、参数的优化设计,利用正交设计试验,分析出风口倒角半径、风道轴线与X轴角度、出风口高度、回风口形式这四个试验因素对乘员舱内监测点温度值的影响,选择最优水平.对结构优化后的乘员舱内的热环境进行数值模拟.
1 重型车乘员舱采暖系统优化设计 1.1 采暖系统设计原采暖系统使用旧燃油预热器为乘员舱内供暖,原采暖系统图如图 1所示,在外界气温低于5 ℃时,在启动发动机之前加温发动机本体、机油、传动油,保证发动机能平稳启动.在发动机平稳启动后,将发动机的高温冷却液引入到加热器前面的乘员取暖器上,并用风扇将采暖系统换热器的热量吹入到乘员舱内.重型车处于停止状态时,低温的冷却液在水泵的作用下,进入到燃油预热器中,进行加热,启动风扇,为乘员舱供暖.然而目前的燃油预热器存在问题,在实际操作过程中,启动燃油预热器后温度上升迅速,温度不可控,容易发生开锅的现象,同时燃油预热器的工作寿命为3 000 h,频繁的操作燃油预热器容易发生问题.因此在实际操作过程中,往往不启动燃油预热器.
本文设计了一种复合式采暖系统改善乘员舱内的乘坐热环境,采暖系统图如图 2所示,燃油预热器工作过程中约有70%的热量通过燃烧区金属壁面传导到冷却水区的冷却液,约25%的热量以废气形式排出,还有5%的能量以辐射形式散发到周围空气中.
循环冷液带走冷却区域内散热量,高温的冷却液从燃油预热器的出水口流出,进入到散热器内,散热器翅片起到强化传热作用,冷却液流经换热器芯体与外界冷空气进行热交换,空气温度升高,低温冷却液在水泵的作用下吸入到燃油加热器的冷却区域,冷却水吸收燃烧产生的热量,反复循环工作,冷却液的温度控制在85 ℃左右,当冷却液温度超过85 ℃后,燃油加热器的油泵停止工作.高温的冷却液送入到散热器芯体后,在风扇的作用下,被加热的空气吹入到乘员舱内,将乘坐热环境进行优化改进.
1.2 风道结构设计设计了两种类型的风道结构,具体结构如图 3所示.类型1风道结构设置为两个出风口的风道,左右两侧风道是对称布置的,入口截面尺寸为150 mm×110 mm,倒角半径为30 mm,出口截面为正六边形,边长为108 mm,倒圆角为10 mm.类型2风道结构设置为四个出风口的风道,左右两侧的风道呈对称布置,小的入口截面尺寸为48 mm×48 mm,大入口截面尺寸为150 mm×55 mm,出口截面分为115mm×179 mm的六边形和边长为56 mm的正六边形,倒圆角为10 mm.
通过调节采暖系统的档位控制出风口的速度,在平均环境温度为14.2 ℃时的试验数据如图 4所示.每隔60 s取一次数据进行数据采集,该实验数据作为对风道结构数值模拟的边界条件.
以1档的边界条件对2种类型的风道结构进行数值模拟,选择出合理的风道结构安装与采暖系统中.
根据图 5可知,类型1的风道结构在拐弯处产生负压,产生涡流现象,将出风口设置为渐扩形,减小拐弯处的压降.类型2的风道结构最高压力升高,主要是由于风道形状及扭曲程度造成风道管壁阻滞造成的,而设置成4通道的类型2风道结构时,低压区降低,因此在风道处理过程中,应尽量将风道设置为平滑曲面,减少较高或较低压区的存在,从而减小由于风道结构对出风口风速的影响.
如图 6所示,在拐弯处会产生气体涡流现象,类型1比类型2的风道结构的气体的涡流现象有明显的改善.根据上述分析发现在换热器的风扇处于一档工作时,类型1风道的气体流动效果更好.
综合对比后,再分别在2档、3档条件下对风道结构进行数值模拟,经过比较后,类型1的风道结构更适合安装在采暖系统中.
2 重型车乘员舱热环境优化设计仿真分析 2.1 乘员舱几何模型验证 2.1.1 模型建立对某重型车乘员舱按1:1的尺寸比例利用CATIA进行三维建模,由于主要分析乘员舱内部的流场及温度场,研究对象为乘员舱内流场,而发动机舱、车轮等结构对乘员舱内流场的数值模拟没有影响,因此可以忽略.为了减少网格数量,提高网格质量,减少模拟时间,对车头、车尾进行简化.因此忽略车外零部件的影响,主要对乘员舱内部进行建模,并且对内部结构进行简化,如图 7所示.
坐标系定位如下:X表示车体宽度方向,从乘员舱内的圆柱的圆心指向右侧;Y表示车体长度方向,从乘员舱内的圆柱的圆心指向前方;Z表示车体高度方向,从乘员舱内地板表面上的圆柱的圆心指向上方.
2.1.2 乘员舱热环境实验分析对乘员舱内的热环境进行寒区冬季热环境实验,获得乘员舱内温度变化的情况,利用实验数据作为数值模拟的边界条件和验证数值模拟结果合理性监测点.
实验测试用车一辆,其长宽高分别为:7.66 m、3.4 m和3.2 m,车体表面材料为钢板,车辆定员3人(驾驶员、乘员1和乘员2).成年男子,衣着厚棉衣和棉裤.
根据实验标准要求,15 min内温度的相对变化率在5%以内,可结束实验.本次实验多点记录仪数据采集的时间间隔为30s.实验前将重型车在室外停放了一夜.实验时将车辆停放在阳光充足的场地,且重型车满足以下条件:
(1) 重型车内表面干净整洁,设备良好完整,停放在阳光充足的场地;
(2) 室外大气温度约-9.7 ℃左右,天气晴朗,阳光充足,微风;
(3) 各个舱门处于打开状态.
为体现出的工作过程中乘员舱温度分布实际情况,在实验过程中各个乘员舱门处于打开状态.乘员舱内的布点情况如图 8所示.
其中测点:1~8测点用于测量贴近动力舱的壁面温度;9测点用于测量乘员2后方电台表面温度;10~12测点用于测量乘员1周围的仪表表面温度;13~14测点分别测量乘员1、乘员2头部周围温度;15~16测点用于测量炮架上测点温度;17~18测点分别测量乘员1、乘员2脚部周围温度;19测点用于测量炮台上温度;20~21测点分别测量乘员2周围仪表表面温度;22~23测点用于测量驾驶员后方空气温度;24测点用于测量驾驶员头部周围空气温度;25~26测点用于测量驾驶员脚部温度.
2.1.3 边界条件及相关参数设置乘员舱内流场的流动是三维非稳态流动,为了保证模拟结果与实验结果吻合,准确的流动模型和边界条件是十分必要的.
边界条件[5]是指在求解域的边界上所求解的变量或一阶导数随地点时间变化的规律,边界条件是求解计算区域的必要条件,合理的边界条件是保证模拟准确性的重要因素.
靠近动力舱壁面为定温边界条件,温度-6.4 ℃;进出口门为压力出口,压力为-0.05 Pa,温度值为-9.7 ℃;乘员舱壁面为对流换热边界条件,温度为-9.7 ℃,对流换热系数为10.8 W·(m-2·K);出风口为压力出口,温度11.2 ℃,压力为-0.05 Pa;回风口为速度入口,温度为3 ℃,v=0.2 m·s-1.
2.1.4 模拟结果与实验结果验证在数值模拟过程中监测驾驶员、乘员1和乘员2的头部及脚部周围测点、乘员舱壁面测点、仪表表面测点进行监测温度值变化情况,将数值模拟结果与实验结果进行对比,具体数值如图 9所示.
根据实验值和模拟值的对比可以看出,除个别点以外,实验值和模拟值误差在10%之内,误差值基本符合实验标准要求,测点温度实验值与模拟值均在合理误差范围内.说明建立的三维乘员舱模型是基本符合要求,因此可以利用该模型对乘员舱内的模型进行优化设计数值模拟.
2.2 乘员舱热环境的正交试验优化设计影响乘员舱内流场、温度场分布的因素包括出风口位置、出风口形状(角度)以及回风口形状(面积)等,为分析各个因素的乘员舱内流场、温度场分布影响规律,本文采取正交试验设计方法[6]进行研究.通过正交试验设计进行优化设计时,需要考虑三个参数[7]:① 试验指标:研究的内容是改善乘员舱内的热环境,那么乘员舱舱内测点的温度就是试验指标; ② 试验因素:A出风口倒角半径、B风道轴线与X轴角度、C出风口高度、D回风口类型; ③ 水平:在正交试验中,试验因素可取不同的值,这种不同的取值状态称为因素水平,简称为水平.通常用下标1、2、3……表示.
2.2.1 正交试验设计针对所设计的采暖系统,选择最优化的出风口位置、出风口角度及回风口面积,对整车内的乘员、驾驶员头部和脚部温度进行试验,温度即为试验指标.具体参数如表 1所示.
本次试验选取的水平数为4,因素为4,若是进行全面试验,则需要进行256次试验.常用的正交试验表中的4水平的有L16(45),L64(421),没有4水平4因素的正交试验表,又不存在交互作用,综合考虑选择L16(45),即仅需进行16次试验即可完成优化设计.根据表 2的试验方案分别进行数值模拟,试验方案是按照L16(45)组织试验.
如表 3所示获得在不同的正交试验组合时获得的乘驾人员头部、脚部测点温度值.确定试验因素的优水平和最优水平组合.
本文设计的正交试验中,涉及试验因素分别为A出风口倒角半径、B风道轴线与X轴之间角度、C出风口高度和D回风口类型.每个试验因素有4个水平,以试验因素A为例,分别为A1、A2、A3和A4.计算出各试验某一因素同一水平所对应的试验指标之和ki值、试验指标之和的平均值ki值、极差试验指标R值.计算的结果如表 4~表 9所示.
考察A因素对4组试验结果的影响,根据表 7可知:k2 < k3 < k1 < k4(出风口倒角半径),k3 < k4 < k1 < k2(风道轴线与X轴夹角),k2 < k3 < k4 < k1(出风口高度),k3 < k2 < k4 < k1(回风口类型),由于驾驶员脚部测点的温度应当尽量的高一些,因此可以判断出因素A、B、C和D的优水平为A4、B2、C1、D1.所以优水平组合为A4B2C1D1.
2.2.2 试验结果分析本文优化设计涉及4个试验指标,单独分析出的优化条件不一致,需要考虑因素的影响的主次顺序,综合考虑,确定出各个因素最佳的参数条件.采取综合平衡法确定最优水平.
由于驾驶员、乘员头部的舒适温度为12~15 ℃,因此监测点温度应尽量低一些.因此驾驶员头部、乘员1头部和乘员2头部温度越低越好,而驾驶员脚部、乘员1脚部和乘员2脚部温度越高越好.
其中乘员1头部、驾驶员脚部测点温度A为主要因素,驾驶员头部、乘员2头部测点温度以B为主要因素,乘员1脚部测点温度以D为主要因素、乘员2脚部测点温度以C为主要因素.
对于因素A:是乘员1头部、驾驶员脚部测点温度影响的第一主要因素.A1对乘员1头部、乘员2头部测点温度较好,A2对驾驶员头部测点温度较合理,A4对驾驶员、乘员1、乘员2脚部测点温度较合理.因此不能确定出A的水平值,应进行计算对比.(A1与A2对比)当取A1时,驾驶员头部测点温度降低3.28%(有利),乘员1头部测点温度上升1.38%(不利),乘员2头部测点温度上升1.22%(不利),驾驶员脚部测点温度升高7.71%(有利),乘员1脚部测点温度上升7.14%(有利),乘员2脚部测点温度下降0.34%(不利).(A1与A2对比)当取A2时,驾驶员头部测点温度上升3.40%(不利),乘员1头部测点温度下降1.36%(有利),乘员2头部测点温度下降1.23%(有利),驾驶员脚部测点温度下降8.34%(不利),乘员1脚部测点温度下降7.70%(不利),乘员2脚部测点温度上升0.34%(有利).A1与A2效果更好,选择A1.
(A1与A4对比)当取A1时,驾驶员头部测点温度不变,乘员1头部测点温度上升4.89%(不利),乘员2头部测点温度上升3.2%(不利),驾驶员脚部测点温度下降7.71%(不利),乘员1脚部测点温度下降5%(不利),乘员2脚部测点温度下降8.5%(不利).很明显A1与A4相比,A4效果更加理想.
对于因素B:对于因素B是驾驶员头部测点、乘员2头部测点温度的影响的第一主要因素.(B3与B4对比)当取B3时,驾驶员头部测点温度上升5.19%(不利),乘员1头部测点温度下降0.75%(有利),乘员2头部测点温度下降5.54%(有利),驾驶员脚部测点温度下降0.3%(不利),乘员1脚部测点温度下降0.37%(不利),乘员2脚部测点温度下降2.05%(不利).当取B4时,驾驶员头部测点温度下降4.94%(有利),乘员1头部测点温度上升0.75%(不利),乘员2头部测点温度上升5.86%(不利),驾驶员脚部测点温度上升0.3%(有利),乘员1脚部测点温度上升0.37%(有利),乘员2脚部测点温度上升2.01%(有利).经对比B4比B3效果更理想.
对于因素C:对于乘员2头部测点温度、驾驶员脚部测点温度而言,因素C排在最后,因此是次要因素.对于驾驶员头部测点温度、乘员1头部测点温度,因素C是第二重要因素.而且取C3分别为两个温度值的优水平,因此,应取C3值.
对于因素D,是乘员1脚部测点温度的主要因素,D3为优水平,同时,对于乘员2头部测点温度,因素D是第二重要因素,而且取D3为优水平.因此,应选择D3值.
因此,本次正交试验优化后的参数条件为A4B4C3D3.由因素水平可知,此时出风口倒角半径为60mm,风道轴线与X轴之间的夹角为110°,出风口高度为570mm,回风口形式为Z-6-280×5,最后以该参数作为优化设计的结构参数.
2.3 最终优化方案数值模拟 2.3.1 边界条件车体表面壁面边界条件为对流换热,温度为-15 ℃,对流换热系数为10.8 W·m-2·k-1;车门壁面边界条件为对流换热,温度为-15 ℃,对流换热系数为14.5 W·m-2·k-1风道表面、暖风机壳体表面的壁面边界条件为绝热;回风口为压力出口,出口压力为-0.05 Pa,为了壁面产生回流现象;在1档时,出风口为速度入口,温度为26.6 ℃,速度为0.25 m·s-1;在2档时,出风口的温度为34.6 ℃,速度为0.56 m·s-1;在3档时,出风口的温度为43.6℃,速度为0.76 m·s-1.
2.3.2 监测截面位置如图 10所示,乘员舱内建立9个监测截面位置,X轴方向建立三个监测面,分别为X=-0.5 m,X=0 m以及X=0.5 m,这3个截面位置分为为乘员1、驾驶员和乘员2座椅纵向中线的截面.Y轴方向建立4个监测面,分别为Y=-1.52 m、Y=0 m和Y=0.38 m以及Y=0.7 m.
采用不均匀系数作为乘员舱热环境内气流分布的评价指标,分析在最优化试验设计结果的温度的均匀性情况.温度和速度对应的不均匀系数计算公式[8]如下:
${k_T} = \frac{{{\sigma _T}T}}{{\bar T}}$ | (1) |
${k_v} = \frac{{{\sigma _v}}}{{\bar v}}$ | (2) |
式中:T为温度平均值;σT为温度均方根偏差;v为速度平均值;σv为速度均方根偏差.
为了分析乘员舱内温度场和流场分布均匀性情况,计算不同截面的温度不均匀系数、速度不均匀系数.一般情况下希望不均匀系数小于1,这样可以获得较为合理的乘员舱流场分布.利用CFD数值模拟的软件的报告功能获得不同截面的表面平均值和表面标准偏差数值,计算出不同表面的温度、速度不均匀系数.
2.3.4 数值模拟结果分析通过数值模拟得出乘员舱在建立的9个监测截面位置不同档位的温度分布云图和温度、速度不均匀系数图,通过分析可知环境温度为-15 ℃时,在二档工况下乘员舱内的热环境更适宜于乘员的乘坐,而且能够基本满足驾乘人员的热舒适性.
同时对驾驶员、乘员1、乘员2头部脚部的温度进行监测,监测测点温度随时间的变化情况,其中驾驶员头部测点坐标(0,-1.72.0.52),驾驶员脚部坐标(-0.45,-3.02, 0.24),乘员1头部坐标(0.52, 0.2, 0.95),乘员1脚部坐标(0.63,-0.63, 0.515),乘员2头部坐标(-0.52, 0.2, 0.95),乘员2脚部坐标(-0.63,-0.63, 0.515).
记录数值模拟过程中0~600 s乘员舱内温度变化情况.在采暖系统温度工作后,乘驾人员头部、脚部温度迅速上升,记录在不同时刻温度值的变化情况.数值模拟过程中的时间步长为1 s,最大内部迭代次数为10步,最大物理时间为600 s.具体结果如图 11所示.采暖系统稳定工作10分钟左右后,开始对驾乘人员头部、脚部测点温度监测.
由图 11可知,驾驶员头部测点温度在300 s后温度趋于稳定为15.4℃,驾驶员脚部测点温度在420 s后温度趋于稳定,其值为10.1 ℃;乘员1头部测点温度在360 s后趋于稳定,其数值为17.4 ℃,480 s后乘员1脚部温度稳定在9.8 ℃;乘员2头部测点温度在420 s后温度值稳定在17.2 ℃,脚部测点温度在480 s后温度在8.8 ℃.
使用旧采暖系统时乘员舱内热环境情况极差,使用旧采暖系统时测点周围均处于0 ℃以下,这对驾驶员、乘员的热感觉影响很大,而采用新采暖系统时乘员舱的测点温度能够基本满足驾乘人员的热舒适性.
3 结论设计一种采暖系统,利用燃油预热器对发动机的冷却水进行预热,将高温冷却水送入到采暖系统的换热器内,达到改善乘员舱内热环境的目的,同时设计一种结构合理的风道结构,使采暖系统中换热器的散热量充分吹到人体周围,达到改善乘员舱乘坐热环境的目的.
为改善乘员舱内热环境,对采暖系统的结构参数进行优化设计,针对回风口形状、出风口倒角半径、风道轴向与X轴角度、出风口高度等四个因素,根据正交试验理论,确定出最优化方案设计,即出风口倒角半径为60 mm,风道轴线与X轴之间的夹角为110°,出风口高度为570 mm,回风口形式为Z-6-280×5.对最优化方案进行CFD数值模拟,说明最优化设计的采暖系统对改善驾驶员、乘员1和乘员2乘坐热感觉起到了明显的作用.
[1] |
邓峰. 汽车内流场与风窗玻璃除霜除雾研究[D]. 长沙: 湖南大学, 2009. DENG Feng.Research on airflow field in car cabin and defrost and defog of Windshield[D].Changsha:Hunan University, 2009. |
[2] |
熊可嘉, 杨志刚. 轿车空调车室内流动换热计算[J]. 同济大学学报(自然科学版), 2011, 39(5): 716 XIONG Kejia, YANG Zhigang. Simulation of Flow Field and Heat transfer in Passenger Compartment[J]. Journal of Tongji University (Natural Science), 2011, 39(5): 716 |
[3] |
ANSI/ASHRAE. Thermal environmental conditions for human occupancy[M]. Atlanta: American Society of Heating, Refrigerating and Air Conditioning Engineering Inc, 1992
|
[4] |
Alahmer A, Mayyas A, Mayyas A A, et al. Vehicular thermal comfort models:a comprehensive review[J]. Applied Thermal Engineering, 2011, 31(6-7): 995 DOI:10.1016/j.applthermaleng.2010.12.004 |
[5] |
温正, 石良辰, 任毅. FLUENT流体计算应用教程[M]. 北京: 清华大学出版社, 2009 WEN Zheng, SHI Liangchen, REN Yi. FLUENT fluid computing applications[M]. Beijing: Tsinghua University Press, 2009 |
[6] |
百度. 正交试验设计及结果分析[EB/OL]. [2014-10-30]. http://wenku.baidu.com/view/48df8302a6c30c2259019ec0.html?re=view. Baidu.Design and result analysis of orthogonal experiment.[DB/OL].[2014-10-30].http://wenku.baidu.com/view/48df8302a6c30c2259019ec0.html?re=view. |
[7] |
任露泉. 试验优化技术[M]. 北京: 机械工业出版社, 1987 REN Luquan. Test optimization techniques[M]. Beijing: China Machine Press, 1987 |
[8] |
李先庭, 赵彬. 室内空气流动数值模拟[M]. 北京: 机械工业出版社, 2009 LI Xianting, ZHAO Bin. Numerical simulation of indoor air flow[M]. Beijing: China Machine Press, 2009 |