CN113916389A - 汽车车内温度计算方法 - Google Patents

汽车车内温度计算方法 Download PDF

Info

Publication number
CN113916389A
CN113916389A CN202111409688.5A CN202111409688A CN113916389A CN 113916389 A CN113916389 A CN 113916389A CN 202111409688 A CN202111409688 A CN 202111409688A CN 113916389 A CN113916389 A CN 113916389A
Authority
CN
China
Prior art keywords
temperature
vehicle
air
heat exchange
energy
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202111409688.5A
Other languages
English (en)
Inventor
李玉钦
张良
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai Pff Electronic Technology Co ltd
Original Assignee
Shanghai Pff Electronic Technology Co ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shanghai Pff Electronic Technology Co ltd filed Critical Shanghai Pff Electronic Technology Co ltd
Priority to CN202111409688.5A priority Critical patent/CN113916389A/zh
Publication of CN113916389A publication Critical patent/CN113916389A/zh
Priority to CN202210670990.4A priority patent/CN115031861A/zh
Priority to CN202210670989.1A priority patent/CN115031860A/zh
Priority to CN202210670999.5A priority patent/CN115165127A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01KMEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
    • G01K1/00Details of thermometers not specially adapted for particular types of thermometer
    • G01K1/14Supports; Fastening devices; Arrangements for mounting thermometers in particular locations
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60HARRANGEMENTS OF HEATING, COOLING, VENTILATING OR OTHER AIR-TREATING DEVICES SPECIALLY ADAPTED FOR PASSENGER OR GOODS SPACES OF VEHICLES
    • B60H1/00Heating, cooling or ventilating [HVAC] devices
    • B60H1/00642Control systems or circuits; Control members or indication devices for heating, cooling or ventilating devices
    • B60H1/00735Control systems or circuits characterised by their input, i.e. by the detection, measurement or calculation of particular conditions, e.g. signal treatment, dynamic models
    • B60H1/00807Control systems or circuits characterised by their input, i.e. by the detection, measurement or calculation of particular conditions, e.g. signal treatment, dynamic models the input being a specific way of measuring or calculating an air or coolant temperature
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T10/00Road transport of goods or passengers
    • Y02T10/80Technologies aiming to reduce greenhouse gasses emissions common to all road transportation technologies
    • Y02T10/88Optimized components or subsystems, e.g. lighting, actively controlled glasses

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Thermal Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Mechanical Engineering (AREA)
  • Air-Conditioning For Vehicles (AREA)
  • Radiation Pyrometers (AREA)
  • Glass Compositions (AREA)

Abstract

本发明涉及汽车智能化技术领域,公开了一种汽车车内温度计算方法,包括如下步骤:(1)分析初始状态车内温度数据,得到车内各部分实时温度;(2)输出当前温度;(3)计算环境与车内换热量Q1;(4)将Q1换算成车内温度变化值;(5)根据变化值得到车内当前温度;(6)在预定的时间间隔内循环步骤(2)‑(5)。本发明充分分析汽车内外环境的因素,精准计算出汽车车厢各部位的实时温度,为精确调控提供数据支撑。

Description

汽车车内温度计算方法
技术领域
本发明涉及汽车智能化技术领域,尤其涉及一种汽车车内温度计算方法。
背景技术
随着生活水平的提高,汽车的乘坐舒适性要求越来越高。汽车车厢内的温度控制直接关系到乘客的乘坐体验,现有的汽车温度控制一般通过车内空调进行调节,而反馈信息是通过车内和空调出风口的温度传感器得到的,通过简单的温度控制,实现车内温度的实时调节。
现有这种温控方式,无法精准地对车厢内各部分的温度进行精确计算,经常出现头部温度太高腿部温度太低或者类似的情况出现,影响乘坐体验。另外,由于车厢环境受车外天气、车速、车内仪器仪表的影响,也会对车厢温度产生不容忽视的影响,这些因素不加考虑的话,并不能很好地对车厢温度实现真正控制。
发明内容
本发明的目的是为了解决上述问题,提供一种汽车车内温度计算方法,充分分析汽车内外环境的因素,精准计算出汽车车厢各部位的实时温度,为精确调控提供数据支撑。
本发明采取的技术方案是:
一种汽车车内温度计算方法,其特征是,包括如下步骤:
(1)分析初始状态车内温度数据,得到车内各部分实时温度;
(2)输出当前温度;
(3)计算环境与车内换热量Q1;
(4)将Q1换算成车内温度变化值;
(5)根据变化值得到车内当前温度;
(6)在预定的时间间隔内循环步骤(2)-(5),
所述步骤(1)中的车内温度数据通过车内温度传感器获得,
所述步骤(3)中的环境包括车辆环境、车辆空调和车辆内饰,
所述步骤(4)中,通过热力学能量温度公式Q=C×M×ΔT,得到车内温度变化值,
所述步骤(5)中,将车内温度变化值与当前温度的运算得到新的当前温度。
进一步,所述步骤(2)中的输出温度,包括车辆前排左右侧头部呼吸点位置温度和后排头部呼吸点温度。
进一步,所述步骤(3)中的车辆环境的换热量包括车体换热量、车体漏风引入的能量、人体发热量、电器零件发热量。
进一步,所述车体换热量包括非玻璃材质车体和玻璃材质车体的换热量,所述玻璃材质车体的换热包括玻璃的热交换和玻璃透光的换热量,所述非玻璃材质车体的换热量和玻璃材质车体的玻璃热交换,是通过热力学能量温度公式计算的,所述和玻璃透光的换热量,通过确定阳光照射信息进行计算,包括直射换热量和散射换热量。
进一步,所述步骤(1)中,所述车内温度数据通过车内温度传感器和空调出风温度传感器的计算得到,当环境温度大于20度时,选择两者之中数值较大数据作为车内温度数据;当环境温度小于0度时,选择两者之中数值较小数据作为车内温度数据;当环境温度位于0至20度之间时,将车内温度传感器的温度作为车内温度数据。
进一步,阳光直射换热量和散射换热量由实时阳光入射强度与阳光高度角通过三角函数关系计算获得。
进一步,车体的车顶温度根据阳光强度、车体颜色吸收系数和最高温度比例系数计算得到,所述
Figure BDA0003373843850000031
Figure BDA0003373843850000032
进一步,所述玻璃材质车体通过玻璃表面温度和室内温度的温差进行热交换,所述玻璃表面的温度通过阳光强度、投影面积、玻璃透光率和留存率计算得到,所述留存率=3×e-0.09×车速-0.5。
进一步,车辆空调的等效出风温度DAT_e=(DAT_r-DAT_(e_prv))×((Vair*ρ_air)/(Vair_Dat)×Rate_space)+
DAT_(e_prv),其中,DAT_r为实时出风温度,DAT_(e_prv)为上一时刻等效出风温度,Vair为每秒空调的出风风量,Vair_Dat为接受空调空间的风量,ρ_air为空气密度,Rate_space为接受空调的空间和整个车载空间的比例。
进一步,根据空调出风的模式,从空调箱直接送入车量后排的能量分别为吹面模式25%、吹面吹脚模式30%、吹脚模式25%、吹脚除模式30%、除霜模式32%、吹面除霜模式35%、全吹模式45%。
本发明的有益效果是:
(1)从各个方位考虑汽车车厢温度影响因素,通过合理计算,得到车厢各部位的精确温度;
(2)综合考虑车内车外各热量源的特性,逐项进行分析计算;
(3)通过定时实时计算,使汽车车厢内实现动态温度控制的最佳状态。
附图说明
附图1是双驱阳光传感器的特性曲线图。
具体实施方式
下面对本发明汽车车内温度计算方法的具体实施方式作详细说明。
从总体上来讲,本发明的汽车车内温度计算方法,是通过热传递原理进行计算的,具体包括如下步骤:
(1)分析初始状态车内温度数据,得到车内各部分实时温度,车内数据通过车内温度传感器获得,温度传感器可以是车内各部位的温度传感器以及空调出风位置的温度传感器。
(2)所当前温度作为此时刻的实时温度,并输出当前温度。
(3)进入下一轮温度计算后,先计算环境与车内换热量Q1,环境因素包括车辆环境、车辆空调和车辆内饰等等各种因素。
(4)将Q1换算成车内温度变化值,可通过热力学能量温度公式Q=C×M×ΔT,得到车内温度变化值,各变化值根据具体情况逐一进行分析计算。
(5)根据变化值得到车内当前温度,将车内温度变化值与当前温度的运算得到新的当前温度,这个过程针对车厢内各部位单独进行计算。
(6)在预定的时间间隔内循环步骤(2)-(5),实时温度实时控制,达到精准控温的效果。
下面具体介绍温度计算过程。考虑各种影响车内温度的因素,分别进行计算。
计算过程中的输入量:车内温度传感器测量值及有效信号、环境温度传感器测量值及有效信号、左侧实际出风温度测量值及有效信号、右侧实际出风温度测量值及有效信号、后排实际出风温度测量值及有效信号、车内空气湿度、车外空气湿度、空调前排出风风量、空调后排出风风量、各个位置的座椅占用位置、内外循环风门开度、车速、左侧阳光强度值、右侧阳光强度值、系统时间、GPS信号、各个车门开启信号、前排左侧电子元器件发热量、前排右侧电子元器件发热量、后排电子元器件发热量、前排空调出风模式、鼓风机工作状态反馈、阳光高度角、阳光与车辆行驶方法的水平面夹角,等等。
计算结果的输出量:前排左侧头部呼吸点拟合温度、前排右侧头部呼吸点拟合温度、后排头部呼吸点拟合温度。
计算过程包括数据预处理、左侧区域温度拟合、右侧区域温度拟合、后侧温度拟合和校准等几部分。
1、数据预处理
1.1、将对应汽车的所有参数进行初始化处理,并将数据限制在合理的数值范围内。
1.2、依据汽车前排出风风量的分配比例,计算出前排左侧出风风量和前排右侧出风风量。
1.3、出风温度等效化处理:
由于每一秒内,经过空调加热或冷却的风并不是立即和周边空气进行热交互,所以,算法获取的实时出风温度不能代表当前和车舱内大部分空气换热时的温度。所以需要获取一个等效的换热温度作为当前的出风温度。所以假设,被空调送出的空气在一定的空间体积内被推送,这个空间标记为第一空间,当推送至这个空间边界外才与其他空气或车内物体进行热交互,第一空间边界外的空间,标记为第二空间。在第一空间内,空调提供的能量仅在空调送出的空气内部进行热交换。
第一空间的体积和空调的送风模式有关,在不同的送风模式下,第一空间在温度拟合的第二空间的占比(Rate_space)是不同的,为经验数据。大致比例如下:
吹面模式35%;
吹面吹脚50%;
吹脚模式67%;
吹脚除霜50%;
除霜40%;
吹面除霜40%;
全吹模式60%;
关闭模式100%。
通过每秒空调的出风风量Vair和空间①Vair_Dat的比例关系,以及上一时刻等效出风温度DAT_(e_prv)和实时出风温度DAT_r的温度计算出,当前时刻的等效出风温度DAT_e。计算公式如下:
DAT_e=(DAT_r-DAT_(e_prv))×((Vair*ρ_air)/(Vair_Dat)×Rate_space)+DAT_(e_prv)
其中ρ_air为空气密度
1.4、阳光强度计算
由于早期控制器计算能力的限制,导致阳光传感器需要通过光学特性进行一定的物理模拟,模拟阳光在不同的角度下的阳光入射强度输入。但是由于现代半导体技术的进步,处理器的运算能力已经十分充沛。软件可以在获知阳光入射角度的前提下直接拟合出车内的阳光入射能量。所以对于阳光强度的输入,已经从阳光传感器的输入转化为对于直接的阳光强度输入。在这个过程中,原有阳光传感器的光学特性在本算法中成为一种额外修正情况,需要针对传感器输入的阳光长度进行逆向修正,从而获取真正的阳光强度。
本算法中,对于阳光传感器的选择不存在双区单区的选择。但是阳光传感器一样需要进行阳光输入数据的修正。
本示例中,参见附图1,针对双驱阳光传感器进行修正。典型的双区阳光传感器的如图1所示。从图中可以看出双驱阳光传感器的特性曲线呈M型。传感器输入峰值出现在单侧。因此可以通过输入的阳光传感器数值,配合阳光角度反向求解出实际的阳光输入强度。由于左区阳光传感器和右区阳光传感器的曲线是对称的。所以在实际使用中,对两组逆向后的阳光强度需要选择数值较大的数据作为真实的直接阳光输入。
1.5、漏风量计算
当车辆车门关闭的时候,车辆内部几乎是密闭空间。外界空气无法直接进入车内,因此只需要考虑车辆缝隙造成的漏风量。但是如果车门打开时,车外空气可以直接进入车内,所以需要针对车门打开时进行独立的漏风量计算。为了简化这个漏风量的计算,一般使用一个固定标准量作为车门打开时的漏风量。
车门关闭时的漏风量可以通过整车密封性实验数据获得。
车门开启时的漏风量可以在经验数据的基础上,通过实验数据矫正后得到。
1.6、风量透过率以及透过风量计算。
由于前排空调送出的风经过本区换热后,会流经后排空间进行2次换热,从而影响后排温度的计算。由于前后排空间并非封闭,所以从前排流经后排的空气可以被分为两种状态。第一种,充分和前排物体换热温度基本与后排温度一致的空气。第二种,经由空调鼓风机吹出,在本区内部分换热,流动到后排区间和后排空气充分换热的空气。针对第一种类型的空气,由于前后排的温差较小,不做计算与考虑。针对第二种类型的空气,则首先需要计算出空气的风量,简称透过风量。
针对透过风量的计算,需要考虑当前区域的空调出风模式以及空调出风风量。整个计算过程可以分为三步,第一步,依据空调的出风模式查表获得空气基本透过率,第二步,通过本区空调的出风风量查找获得风量的修正率。第三步将空气基本透过率乘以风量修正率从而获得风量透过率。
将实时的空调出风量乘以风量透过率,就可以计算出经由前排空调送出的且在本区空间内部分换热后流至后排空间换热的风量。
空气基本透过率经验值:
吹面模式 35%;
吹面吹脚 30%;
吹脚模式 25%;
吹脚除霜 30%;
除霜 32%;
吹面除霜 35%;
全吹模式 45%;
关闭模式 100%。
风量修正率表格:
风量 比例
0 0.5
50 0.6
100 0.75
150 1
200 1.10
250 1.15
300 1.18
350 1.20
400 1.22
450 1.25
500 1.15
2、左侧区域温度拟合。
2.1、车内温度初始化。
在模型调用的初期,需要针对模型内部的参数进行初始化。其中,软件拟合的车内温度其自身也需要初始化。初始化数据的选择一般采用环境温度或者车内其它的温度传感器,其中出空温度传感器就是一个比较好的选择。帮我们空调系统自身包含车内温度传感器时,模型调用初期这个温度值也是一个非常好的选择。
首先,需要针对车内传感器是否有效进行判断,如果车内温度传感器不存在,则直接将环境温度的数值赋值给车内温度变量。其后需要对车内出风温度传感器进行判断,如果传感器有效,出风温度传感器数值赋值给车内温度变量,如果传感器无效,则保持车内温度变量不变。
需要对环境温度进行判断,如果环境温度大于20度,选择车内温度变量和出风温度传感器读数两者之中数值较大数据作为初始温度。如果环境温度小于0度则选择车内温度变量和出风温度传感器读数两者之中数值较小的数据作为初始化温温度。如果环境温度在0度和20度之间,则直接选用车内温度变量的数值作为初始化温度。
2.2、车体环境热交换能量计算。
车体环境热交换的能量可以被分为以下几个部分:
车身框架热交换能量;
车身玻璃侧导入能量;
车身漏风引起的热交换能量;
乘客的发热量;
电器零件的发热量。
综合以上能量,可以计算出整个车体环境热交换的能量数值。
车外的空气密度计算需要2个输入参数,环境温度和车外空气湿度。通过环境温度可以查表获得干空气密度和饱和湿度下空气密度。
然后将车外空气湿度比例乘以干空气密度和饱和湿度下空气密度的差值,再加上干空气密度,即可获得当前湿度下的空气密度。
车内空气密度计算需要2个输入参数,车内温度和车内空气湿度。其中车内温度应选用拟合车内温度。通过车内温度可以查表获得干空气密度和饱和湿度下空气密度。
然后将车外空气湿度比例乘以干空气密度和饱和湿度下空气密度的差值,再加上干空气密度,即可获得当前湿度下的空气密度。
车外空气焓值计算,空气焓用符号i表示。湿空气焓值等于1kg干空气的焓值与dkg水蒸气焓值之和。
湿空气焓值计算公式化:i=(1.01+1.84d)t+2500d(kj/kg干空气)。
式中:t—空气温度℃;
d—空气的含湿量kg/kg干空气;
1.01—干空气的平均定压比热kj/(kg.K);
1.84—水蒸气的平均定压比热kj/(kg.K);
2500—0℃时水的汽化潜热kj/kg。
车内空气焓值计算,空气焓用符号i表示。湿空气焓值等于1kg干空气的焓值与dkg水蒸气焓值之和。
湿空气焓值计算公式化:i=(1.01+1.84d)t+2500d(kj/kg干空气)。
式中:t—空气温度℃;
d—空气的含湿量kg/kg干空气;
1.01—干空气的平均定压比热kj/(kg.K);
1.84—水蒸气的平均定压比热kj/(kg.K);
2500—0℃时水的汽化潜热kj/kg;
车内外空气焓差计算,将车外空气焓值减去车内空气焓值,即可获得车内外空气焓差。
车身框架数据准备。算法中需要将整车部分参数作为算法的重要参数,用来计算车身框架热交换能量。这些数据包含了整车中,整个温度拟合过程中所有需要参数,在左侧的温度拟合时,需要针对左侧的请教进行相关配置。
在计算前排左侧的车内温度时,右侧为开放空间,相关的车门面积,玻璃面积应被赋值为零。左右共有的部分,以实际使用计算空间做划分。
这些数据包括:
车顶面积;
车顶热交换率;
前左侧车门面积;
前左侧车门热交换率;
后左侧车门面积;
后左侧车门热交换率;
前右侧车门面积;
前右侧车门热交换率;
后右侧车门面积;
后右侧车门热交换率;
车身防火墙面积;
车身防火墙热交换率;
车身后侧面积;
车身后侧热交换率;
车身正前方框架面积;
车身正前方框架热交换率;
车身底盘面积;
车身底盘热交换率;
车身颜色修正系数;
车体玻璃数据准备。
算法中需要将玻璃部分参数作为算法的重要参数,用来计算车身玻璃测导入能量。
这些数据包括:
挡风玻璃水平面投影面积;
玻璃热交换率;
前左侧玻璃水平面投影面积;
前左侧玻璃热交换率;
前右侧玻璃水平面投影面积;
前右侧玻璃热交换率;
后左侧玻璃水平面投影面积;
后左侧玻璃热交换率;
后右侧玻璃水平面投影面积;
后右侧玻璃热交换率;
后挡风玻璃水平面投影面积;
后挡风玻璃热交换率;
天窗玻璃面积;
天窗玻璃热交换率;
挡风玻璃面积;
前左侧玻璃面积;
前右侧玻璃面积;
后左侧玻璃面积;
后右侧玻璃面积;
后挡风玻璃面积;
挡风玻璃垂直面投影面积;
前左侧玻璃垂直面投影面积;
前右侧玻璃垂直面投影面积;
后左侧玻璃垂直面投影面积;
后右侧玻璃垂直面投影面积;
玻璃折射率;
仪表台面积;
左侧人体模型面积;
右侧人体模型面积。
其中各个玻璃的水平投影面积或垂直投影面积,应被理解为阳光通过玻璃后,在车身的内部水平方向或垂直方向形成的光斑面积。这些光斑面积与玻璃本身的面积以及玻璃倾斜角度相关。同时,有和太阳高度角和阳光入射角有关。因此,算法中也因将玻璃倾斜角度,太阳高度角和阳光入射角作为一个重要的输入。
阳光长度数据准备:
算法中需要将阳光本体特性作为算法的重要参数用来计算车身玻璃侧导入能量。
这些数据包括:
水平方向上的阳光入射强度;
垂直方向上的阳光强度;
水平方向上的阳光散射强度;
垂直方向上的阳光散射强度;
透明玻璃的阳光透射率;
有色玻璃的阳光透射率;
玻璃矫正系数;
前挡风玻璃水平方向阳光入射强度;
后挡风玻璃水平方向阳光入射强度;
左侧玻璃水平方向阳光入射强度;
右侧玻璃水平方向阳光入射强度。
阳光在水平方向和垂直方向的入射强度与太阳高度角和阳光入射夹角有关,因此需要引入这2个参数作为重要的输入。
此时用于计算的各个方向上强度分量的阳光总强度,就是数据预处理过程中经过逆向处理后的直接阳光强度输入。
其中垂直方向上的阳光强度和水平方向上的阳光强度可以通过直接阳光入射强度与阳光高度角通过三角函数关系计算获得。
水平方向上的阳光散射强度公式如下:
Figure BDA0003373843850000161
其中,地球半径=6371KM;
大气层厚度=1000KM;
大气透射率=0.7;
垂直方向上的阳光散射强度公式如下:
Figure BDA0003373843850000162
其中,地球半径=6371KM;
大气层厚度=1000KM;
大气透射率=0.7。
车体热交换能量计算。
车体热交换能量包括以下几个部分:
车顶热交换能量;
前侧车体热交换能量;
两侧车体热交换能量;
后侧车体热交换能量;
防火墙热交换能量;
底盘热交换能量。
车顶热交换能量:
车顶温度和车速有关同时也和日照强度有关。车顶的换热量需要通过车顶温度和车内的温差以及车顶的换热率来计算。
车顶温度的计算方式如下:
需要通过车速计算出车顶的最高温度比例系数,公式如下:
Figure BDA0003373843850000171
利用阳光强度,车体颜色吸收系数和最高温度比例系数,相乘得到当前车速下的最高玻璃温度。
利用车顶最高温度和前一次的车顶温度的温差,乘以温升比例系数,计算出当前的车顶温度。
前侧车体热交换能量:
利用环境温度和车内温度计算出温差,将温差乘以前侧车体面积,在乘以前侧车体换热系数。得到前侧车体热交换能量。
两侧车体热交换能量:
两侧车体热交换能量包括以下部分:
前排左侧热交换能量;
前排右侧热交换能量;
后排左侧热交换能量;
后排右侧热交换能量;
热交换能量的主要思路:首先计算车门表面的温度,再通过车门表面温度和室内温度的温差进行热交换的量计算。
在计算车门表面温度之前首先先要计算车门表面的阳光功率输入,车门表面的质量,然后利用换热公式计算出每1秒钟的车门表面的温升,加上上一次出门表面温度获得当前时刻的车门表面温度。
热交换公式如下:
Q=C×M×ΔT
其中Q为阳光输入能量;
C是车门金属的比热容;
M是车门金属的总质量;
ΔT是车门表面的前后时刻的温差。
利用车门表面温度和车内温度计算出温差,将温差乘以单侧车体面积,在乘以单侧车体换热系数。得到单侧车体热交换能量。
刚分别求出各测的热交换能量之后,将这些能量求和,就可以获得两侧车体热交换能量。
后侧车体热交换能量:
后侧车体热交换能量计算方式与两侧车体热交换能量计算方式一致。
防火墙热交换能量:
防火墙热交换能量计算时需要用到发动机舱内温度。如果车辆为燃油发动机配置,发动机舱对温度可以使用发动机水温代替。如果车辆为新能源车辆,发动机舱温度可以使用电机温度或电机排水水温。如不精确考虑,可以使用环境温度代替。
利用防火墙表面温度和车内温度计算出温差,将温差乘以防火墙面积,在防火墙结构换热系数。得到防火墙热交换能量。
底盘热交换能量:
底盘热交换能量计算方式和前侧车体热交换能量的计算方式一致。
当计算获得各个部分的热交换能量后,将这些能量求和就可以获得车体与车内空气换热的能量。
由于车内空气受空调驱使在不断的流动,因此需要考虑不同的空调出风量情况下实际的热交换能量。在这个算法中,使用一个等效系数,控制或计算热交换能量的数值。
等效系数的计算:
等效系数的计算方式:通过获取空调的出风总量查表获得对应的等效系数。
风量kg/m3 系数
0 1
50 1
100 1
150 1.1
200 1.7
250 1.9
300 1.9
将求和后的能量乘以等效系数就可以获得实时的车体热交换能量。
车身玻璃侧导入能量:
车身玻璃侧导入能量分3个部分:
玻璃温度热传导能量;
阳光直射输入能量;
阳光散射输入能量。
玻璃温度热传导能量:热交换能量的主要思路:首先计算玻璃表面的温度,再通过玻璃表面温度和室内温度的温差进行热交换的量计算。
玻璃表面的温度和车速有关,所以先计算阳光能量在不同车速下的留存率,计算方式如下:
留存率=3×e-0.09×车速-0.5
在计算玻璃表面温度之前,首先要计算玻璃表面的阳光功率(水平方向上的阳光强度)输入,阳光能量等于水平方向上的阳光强度乘以玻璃水平投影面积加上垂直方向上的阳光强度乘以玻璃垂直投影面积的和,再乘以玻璃的整体的光透过率((100%-折射率)*透过率)与留存率的乘积。
计算玻璃的质量(玻璃面积*玻璃厚度*玻璃密度),
利用换热公式计算出每1秒钟的玻璃表面的温升,加上一次出玻璃面温度获得当前时刻的玻璃表面温度。
热交换公式如下:
Q=C×M×ΔT
其中Q为阳光输入能量;
C是车门玻璃的比热容;
M是车门玻璃的总质量;
ΔT是玻璃表面的前后时刻的温差。
利用玻璃表面温度和车内温度计算出温差,将温差乘以单侧玻璃面积,在乘以单侧玻璃换热系数。得到单侧玻璃热交换能量。
刚分别求出各测的热交换能量之后,将这些能量求和,就可以获得整个车体中玻璃的热交换能量。
阳光直射输入能量:
阳光直射输入的能量通过内饰发热,热传递给车内空气。而这种热传递的途径依据车辆结构不同而有所不同。其中可以归类为2种:仪表板上层和挡风玻璃之间的窄空间中低风速的对流换热,以及开阔空间较高风速的强制换热。
针对对流换热的能量计算,需要通过仪表板和前挡风玻璃之间的空气为计算载体,计算狭窄空间和开阔空间之间的空气热交换。
对于狭窄空间的空气,主要热交换能量有:
仪表板发热量;
通过挡风玻璃和外交的热交换能量;
与开阔空间的空气热交换能量。
仪表板发热量的计算方式:
仪表板发热量=垂直方向的阳光强度×挡风玻璃透过率×玻璃校准系数×仪表水平面积×仪表台的阳光能量留存系数
通过仪表板发热量减去前一次挡风玻璃和外交的热交换能量,可以反向求解出仪表表面的温度:
公式如下:
Q=C×M×ΔT
其中Q为热量差
C是仪表板发热部分的比热容;
M是仪表板发热部分的质量;
ΔT是仪表表面温度前后时刻的温差。
将温差加上前一时刻的仪表表面的温度就可以获得当前时刻的仪表表面的温度。
通过挡风玻璃和外交的热交换能量计算方式:
挡风玻璃和外交的热交换能量
=(环境温度-前一时刻窄空间空气温度)×仪表面积÷cos(挡风玻璃倾斜角度)×挡风玻璃热交换率
由于与开阔空间的空气热交换能量等于仪表板发热量减去挡风玻璃和外交的热交换能量,且利用公式:Q=C×M×ΔT,可以反向求解出当前时刻窄空间空气温度。其中C是空气的比热容,M是空气质量。Q是能量差。ΔT是窄空间空气温度前后时刻的温差。将前一时刻窄空间空气温度加上温差就是当前时刻窄空间空气温度。
当当前时刻窄空间空气温度计算完毕,可以实现对于狭窄空间的空气种主要热交换能量的迭代计算。
从计算过程中,获取传递给开阔空间的能量作为真正传递给车内空气的能量,用于阳光直射输入能量的计算。
开阔空间较高风速的强制换热:
需要考虑前左侧,前后侧,后左侧,后右侧,车顶玻璃侧,后侧挡风玻璃的换热。
换热方式可以采用阳光直射后发热的简单模式完成计算。采用公式:
阳光入射的能量=水平面上光斑面积×水平面上阳光强度×玻璃的透过率×玻璃校准系数
将开阔空间换热能量和狭窄空间换热能量相加后就是:阳光直射输入能量。
阳光散射输入能量:
需要考虑挡风玻璃侧,前左侧,前后侧,后左侧,后右侧,车顶玻璃侧,后侧挡风玻璃的散射输入。
换热方式可以采用阳光直射后发热的简单模式完成计算。采用公式:
阳光入射的能量
=玻璃垂直平面投影面积×阳光散射强度
×玻璃的透过率×玻璃校准系数
这些能量相加就可以获得阳光散射输入能量。
将玻璃温度热传导能量,阳光直射输入能量,阳光散射输入能量相加,就可以获得车身玻璃侧导入能量。
车身漏风引起的热交换能量:
车身漏风引起的热量交换计算方式如下:
车身漏风引起的热交换能量
=车身漏风风量×空气密度×车内外空气焓差
乘客的发热量:
在计算车内乘客的发热量时,需要考虑乘客的数量,以及乘客的位置。如果是驾驶员则选用驾驶员的发热量,如果是普通乘客则使用普通乘客的发热量。
车内乘客在进入或退出车辆时,对于车辆内部空气的热交换并非是直接影响。所以需要针对乘客进出车辆时的热力交换来进行热传递过程的拟合。一般情况下,车辆空调使用过程中空调的出风口是避免对着乘客吹的,虽然出现了强制对流换热,但是整个热交换过程从实验结果来看,基本属于线性过程。因此在算法中使用了一个线性的过渡过程模拟热量逐步传递至车内的过程。
首先,将乘客在单位时间内的发热量处以一定的步长(上升步长),计算出单步递增发热量。上升步长可以通过实验获得。
接着,将乘客在单位时间内的发热量除以一定的步长(下降步长),计算出但不下降发热量。下降步长可以通过实验获得。
设置一个变量,用来计算乘客发热量的输出。将这个变量在前一时刻的数值和当前时刻的乘客发热的输入数值进行比较。如果判定输入的乘客能量在增加,则在变量上累加单步递增发热量。如果判定输入的乘客能量在下降,则在变量上累加单步下降发热量。将累加后的发热量赋值给变量本身,并输出作为当前的乘客发热量。
电器零件的发热量:
电器零件的发热量是由上级程序直接给出,所以可以直接被利用,不需要做相关的处理。
将车身框架热交换能量,车身玻璃侧导入能量,车身漏风引起的热交换能量,乘客的发热量,电器零件的发热量求和后,就获得了车体环境热交换能量。
2.3、空调提供的能量。
在计算空调提供的能量时,并非计算空调本身消耗的能量。可以通过计算空调出风口送出的空气与车内空气进行热交换的能量即可。
数据预处理:
空调箱送出空气的密度;
通过环境温度查出干空气的密度和饱和湿空气的密度。
用饱和湿空气的密度减去干空气的密度乘以湿度,再加上干空气密度即可获得当前湿度下的空气密度。
透风温度计算:
在计算后排的温度时,要考虑前排送入的能量。所以需要知道前排温度拟合区域像后排透过的空气总量(风量)以及透过空气的温度。
首先通过空调出风的模式计算出从空调箱直接送入后排的能量比例。这个比例是一个经验值。
Figure BDA0003373843850000251
Figure BDA0003373843850000261
利用上一次的空调出风口的能量,乘以系数,除以漏风量除以空气比热容,在加上上一次的透风温度。既可以获得当前情况下的透风温度。
空调出风口送出的能量:
空调送出总能量:
将车内拟合温度减去出风口温度,差值乘以出风口风量,乘以空气比热容,就是空调送出的从能量。
通过漏风送往后排的能量:
将透风温度减去出风口温度,差值乘以漏风风量,乘以空气比热容,就是通过漏风送往后排的能量。
将空调送出总能量减去漏风送往后排的能量,就是空调出风口送出的能量。
2.4、内饰换热的能量
内饰换热分2个部分,塑料件换热和金属件换热。
塑料件换热热量:
由于车内空气流动将直接影响车内的空气换热速度。所以需要通过风量信息确定温度拟合区域的塑料零件的换热率比例。这种比率是以稳定换热情况下的换热率为基准的比率。比率的数据可以通过试验获取或标定。
用车内的拟合温度减去前一次的内饰温度,并将差值乘以内饰面积,再乘以基础换热率,再乘以换热比率。及获得内饰和车内空气进行热交换的热量。
用内饰和空气的热交换热量除以内饰塑料件的比热容,除以材料密度,除以材料评价厚度,除以评价面积,并加上上一次的内饰温度即可获得当前时刻的内饰温度。
由于内饰温度计算属于迭代计算。当前时刻的内饰温度将在下一时刻,成为前一时刻的内饰温度。
选择内饰和车内空气进行热交换的热量输出。
金属件换热热量:
金属零件的能量计算方式和内饰换热的方式一致,只是内部的材料特性参数有变化。
输出塑料件换热热量与金属件换热热量。
2.5、车内温度的拟合计算。
首先计算系统空气测的总换热量。
将车体环境热交换的能量,空调提供的能量,内饰换热的能量求和。
将总能量除以空气比热容,除以拟合空间内的空气体积,除以空气密度,加上前一次拟合的车内温度。即得到当前时刻的车内拟合温度。
由于温度计算属于迭代计算。当前时刻的车内拟合温度将在下一时刻,成为前一时刻的车内拟合温度。
输出车内拟合温度。
2.6、由于实际使用过程中左侧和右侧的温度偏差不大。所有不针对左侧右侧的空气热量交换进行计算。如果一定需要增加计算,可以在计算换热总能量的过程中,增加一个左右侧热交换能量。
这个能量计算的方式为:左侧车内拟合温度减去右侧车内拟合温度,差值乘以左右空间相接的垂直面积,乘以空气换热系数。
2、右侧车内温度拟合。
拟合方式同左侧拟合方式。
需要将车辆的车体参数,玻璃参数等物理特性参数从左侧的设定调整为右侧设定。
由于实际使用过程中左侧和右侧的温度偏差不大。所有不针对左侧右侧的空气热量交换进行计算。如果一定需要增加计算,可以在计算换热总能量的过程中,增加一个左右侧热交换能量。
这个能量计算的方式为:右侧车内拟合温度减去左侧车内拟合温度,差值乘以左右空间相接的垂直面积,乘以空气换热系数。
3、后侧温度拟合。
后侧温度拟合方式基本与左侧的拟合方式一致,仅在计算空调输出能量的方式有不同。
后区空调输出能量:
后区空调输出的能量需要考虑一下几个部分:
前排左侧透风带来的能量;
前排右侧透风带来的能量;
空调后排出风口送出的能量;
前排左侧透风带来的能量。
前左侧漏风温度减去后排车内拟合温度,差值乘以空气比热容,乘以左侧漏风风量,即为前排左侧透风带来的能量。
前排右侧透风带来的能量:
前右侧漏风温度减去后排车内拟合温度,差值乘以空气比热容,乘以右侧漏风风量,即为前排右侧透风带来的能量。
空调后排出风口送出的能量:
后排空调出风口的出风温度减去后排车内拟合温度,差值乘以空气比热容,乘以后排出风风量,即为空调后排出风口送出的能量。
将前排左侧透风带来的能量,前排右侧透风带来的能量,后排出风口送出的能量求和,获得后区空调输出能量。
当完成后侧温度拟合后,整车的主要位置的温度拟合均完成。
以上仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。

Claims (10)

1.一种汽车车内温度计算方法,其特征在于:包括如下步骤:
(1)分析初始状态车内温度数据,得到车内各部分实时温度;
(2)输出当前温度;
(3)计算环境与车内换热量Q1;
(4)将Q1换算成车内温度变化值;
(5)根据变化值得到车内当前温度;
(6)在预定的时间间隔内循环步骤(2)-(5),
所述步骤(1)中的车内温度数据通过车内温度传感器获得,
所述步骤(3)中的环境包括车辆环境、车辆空调和车辆内饰,
所述步骤(4)中,通过热力学能量温度公式Q=C×M×ΔT,得到车内温度变化值,
所述步骤(5)中,将车内温度变化值与当前温度的运算得到新的当前温度。
2.根据权利要求1所述的汽车车内温度计算方法,其特征在于:所述步骤(2)中的输出温度,包括车辆前排左右侧头部呼吸点位置温度和后排头部呼吸点温度。
3.根据权利要求1所述的汽车车内温度计算方法,其特征在于:所述步骤(3)中的车辆环境的换热量包括车体换热量、车体漏风引入的能量、人体发热量、电器零件发热量。
4.根据权利要求3所述的汽车车内温度计算方法,其特征在于:所述车体换热量包括非玻璃材质车体和玻璃材质车体的换热量,所述玻璃材质车体的换热包括玻璃的热交换和玻璃透光的换热量,所述非玻璃材质车体的换热量和玻璃材质车体的玻璃热交换,是通过热力学能量温度公式计算的,所述和玻璃透光的换热量,通过确定阳光照射信息进行计算,包括直射换热量和散射换热量。
5.根据权利要求4所述的汽车车内温度计算方法,其特征在于:所述步骤(1)中,所述车内温度数据通过车内温度传感器和空调出风温度传感器的计算得到,当环境温度大于20度时,选择两者之中数值较大数据作为车内温度数据;当环境温度小于0度时,选择两者之中数值较小数据作为车内温度数据;当环境温度位于0至20度之间时,将车内温度传感器的温度作为车内温度数据。
6.根据权利要求4所述的汽车车内温度计算方法,其特征在于:阳光直射换热量和散射换热量由实时阳光入射强度与阳光高度角通过三角函数关系计算获得。
7.根据权利要求4所述的汽车车内温度计算方法,其特征在于:车体的车顶温度根据阳光强度、车体颜色吸收系数和最高温度比例系数计算得到,所述
Figure FDA0003373843840000021
Figure FDA0003373843840000022
8.根据权利要求4所述的汽车车内温度计算方法,其特征在于:所述玻璃材质车体通过玻璃表面温度和室内温度的温差进行热交换,所述玻璃表面的温度通过阳光强度、投影面积、玻璃透光率和留存率计算得到,所述留存率=3×e-0.09×车速-0.5。
9.根据权利要求4所述的汽车车内温度计算方法,其特征在于:车辆空调的等效出风温度DAT_e=(DAT_r-DAT_(e_prv))×((Vair*ρ_air)/(Vair_Dat)×Rate_space)+DAT_(e_prv),其中,DAT_r为实时出风温度,DAT_(e_prv)为上一时刻等效出风温度,Vair为每秒空调的出风风量,Vair_Dat为接受空调空间的风量,ρ_air为空气密度,Rate_space为接受空调的空间和整个车载空间的比例。
10.根据权利要求4所述的汽车车内温度计算方法,其特征在于:根据空调出风的模式,从空调箱直接送入车量后排的能量分别为吹面模式25%、吹面吹脚模式30%、吹脚模式25%、吹脚除模式30%、除霜模式32%、吹面除霜模式35%、全吹模式45%。
CN202111409688.5A 2021-11-25 2021-11-25 汽车车内温度计算方法 Pending CN113916389A (zh)

Priority Applications (4)

Application Number Priority Date Filing Date Title
CN202111409688.5A CN113916389A (zh) 2021-11-25 2021-11-25 汽车车内温度计算方法
CN202210670990.4A CN115031861A (zh) 2021-11-25 2022-06-15 仪表台和挡风玻璃构成的空间的空气温度计算系统及方法
CN202210670989.1A CN115031860A (zh) 2021-11-25 2022-06-15 基于阳光输入能量计算的变色玻璃控制方法及系统
CN202210670999.5A CN115165127A (zh) 2021-11-25 2022-06-15 车辆玻璃温度计算系统及方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111409688.5A CN113916389A (zh) 2021-11-25 2021-11-25 汽车车内温度计算方法

Publications (1)

Publication Number Publication Date
CN113916389A true CN113916389A (zh) 2022-01-11

Family

ID=79247989

Family Applications (4)

Application Number Title Priority Date Filing Date
CN202111409688.5A Pending CN113916389A (zh) 2021-11-25 2021-11-25 汽车车内温度计算方法
CN202210670989.1A Pending CN115031860A (zh) 2021-11-25 2022-06-15 基于阳光输入能量计算的变色玻璃控制方法及系统
CN202210670999.5A Pending CN115165127A (zh) 2021-11-25 2022-06-15 车辆玻璃温度计算系统及方法
CN202210670990.4A Pending CN115031861A (zh) 2021-11-25 2022-06-15 仪表台和挡风玻璃构成的空间的空气温度计算系统及方法

Family Applications After (3)

Application Number Title Priority Date Filing Date
CN202210670989.1A Pending CN115031860A (zh) 2021-11-25 2022-06-15 基于阳光输入能量计算的变色玻璃控制方法及系统
CN202210670999.5A Pending CN115165127A (zh) 2021-11-25 2022-06-15 车辆玻璃温度计算系统及方法
CN202210670990.4A Pending CN115031861A (zh) 2021-11-25 2022-06-15 仪表台和挡风玻璃构成的空间的空气温度计算系统及方法

Country Status (1)

Country Link
CN (4) CN113916389A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115626021A (zh) * 2022-06-15 2023-01-20 上海普法芬电子科技有限公司 汽车虚拟传感器系统

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115626021A (zh) * 2022-06-15 2023-01-20 上海普法芬电子科技有限公司 汽车虚拟传感器系统
CN115626021B (zh) * 2022-06-15 2024-10-11 上海普法芬电子科技有限公司 汽车虚拟传感器系统

Also Published As

Publication number Publication date
CN115165127A (zh) 2022-10-11
CN115031860A (zh) 2022-09-09
CN115031861A (zh) 2022-09-09

Similar Documents

Publication Publication Date Title
US9428035B2 (en) Vehicle HVAC noise control system
US20100019050A1 (en) Automatic Climate Control for a Vehicle
US5988517A (en) Method and system for controlling an automotive HVAC system based on the principle of HVAC work
CN106705348A (zh) 用于改进的热舒适性的基于模型的自动气候控制系统
JPS588419A (ja) 自動車用空調装置
WO2018176598A1 (zh) 一种车辆用太阳能空调的温度补偿方法及系统
KR102017534B1 (ko) 자동차의 내부 온도 조절 방법 및 관련 에어­컨디셔닝 시스템
US5209079A (en) Control apparatus for air conditioner used for vehicles
Ružić et al. Thermal interaction between a human body and a vehicle cabin
JPH08238927A (ja) 曇り防止のために自動車用hvacシステムを制御する方法および制御システム
CN103963598B (zh) 汽车自动空调控制方法
EP0060463B1 (en) Solar radiation heat influx sensor for an automotive vehicle
CN104589951B (zh) 汽车空调循环风门的控制方法及其控制系统
JP2004256101A (ja) 熱流束解析を利用した2ゾーン自動環境制御アルゴリズム
CN1667383A (zh) 内部温度检测方法和装置
CN113916389A (zh) 汽车车内温度计算方法
CN104986011A (zh) 一种汽车双驱自动空调的温度风门控制系统
Shimizu et al. Analysis on air–conditioning heat load of a passenger vehicle
Lahimer et al. Potential of solar reflective cover on regulating the car cabin conditions and fuel consumption
KR101940729B1 (ko) 자동차의 내부 온도의 멀티존 조절 방법 및 관련 에어­컨디셔닝 시스템
Orzechowski et al. Evaluation of thermal conditions inside a vehicle cabin
CN104875577B (zh) 一种全工况双区温度校准方法及其汽车空调控制系统
CN115626021B (zh) 汽车虚拟传感器系统
KR20100017678A (ko) 개선된 일광 보정을 구비한 차량용 공기 조화 시스템
JP3182983B2 (ja) 空調装置

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
WD01 Invention patent application deemed withdrawn after publication
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20220111