CN106777711B - 建立车载变风量空调系统风量预测模型的方法 - Google Patents
建立车载变风量空调系统风量预测模型的方法 Download PDFInfo
- Publication number
- CN106777711B CN106777711B CN201611198508.2A CN201611198508A CN106777711B CN 106777711 B CN106777711 B CN 106777711B CN 201611198508 A CN201611198508 A CN 201611198508A CN 106777711 B CN106777711 B CN 106777711B
- Authority
- CN
- China
- Prior art keywords
- air
- variable
- model
- vehicle
- fresh
- 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.)
- Active
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Air Conditioning Control Device (AREA)
- Air-Conditioning For Vehicles (AREA)
Abstract
建立车载变风量空调系统风量预测模型的方法，属于计算机数值模拟技术领域，尤其是针对车载变风量空调进行模拟，进而得到总风量与新风量的定量关系式的方法。首先确定车载变风量空调系统新风、回风、总送风和风机全压的数学关系，然后建立包括空调机组关键部件和变频风机整车数值计算模型；设计并计算工况，进行整车的气流流动仿真计算；通过变频风机的实际工作点求解变风量空调运行时的新风量、回风量和总风量；利用多元非线性回归分析方法，建立变风量空调系统预测模型。
Description
技术领域
本发明属于计算机数值模拟技术领域，尤其是针对车载变风量空调，特别是轨道车辆使用的变风量空调进行模拟，进而得到总风量与新风量的定量关系式的方法。
背景技术
随着我国交通事业的发展，轨道车辆，包括地铁和铁路客车车内环境舒适性已经成为了一个新兴的研究领域。 GB/T 12817-91《铁道客车通用技术条件》、EN14750-114750-1《Railway Applications Air Conditioning for Urban and Suburban Rolling StockPart I comfort parameters》等标准对车厢内的微风速和温度舒适性等参数进行了规定。
轨道车辆客室内的舒适性与总送风量、额定制冷量以及新风量等因素有关，为达到舒适性要求，需要根据车内负荷变化而调整总风量和新风量；而传统的轨道车辆空调设计一般只考虑总送风量和额定制冷量，无法满足标准规定的舒适性要求。
目前国内外是以车辆内部和送风道的气流组织、以及地铁车站通风系统为研究对象进行仿真和设计。国外对地铁车辆数值仿真的重点是发生火灾时地铁车站的烟气流动和控制、地铁车站通风系统最恶劣工况时的安全通风模式、高速列车车厢内的振动和噪声等问题。如Choquart Clodic等采用CFD软件计算了稳态工况下车辆内部气流流场、压力场的分布及热舒适状况，但是计算模型中未考虑空调系统。国内研究主要集中于地铁车厢和动车组车厢内的气流组织以及送风道的气流组织。周生存等介绍了地铁车辆空调风道设计的特殊性及国内地铁A型车空调的风道设计方法。刘洋等介绍了国内B型车的地铁车辆风道系统，对地铁车辆风道系统出风均匀性设计进行了分析。杨柳通过对风道系统的改进，将风道和车厢合并在一起建立物理模型，进一步进行数值模拟，计算得出风道的优化对车厢内气流组织、人体舒适度的影响。
对空调系统内部气流组织的研究，尤其是对轨道车辆上使用的变风量空调的研究为空白。
发明内容
为了解决现有技术的不足，本发明采用数值模拟的研究方法，以有限体积法为理论基础，以空调机组、风道和车厢作为整体研究对象，建立整车数值计算模型，使用流体计算软件Fluent进行仿真计算并验证计算结果，进而建立风量调节数学模型，获得空调系统内部新风量/回风量与总风量的定量关系，以及新风口与回风口的阀门开度变化对新风量和回风量的影响规律。
本发明采用的技术方案是：建立车载变风量空调系统风量预测模型的方法，包括以下步骤：
1）确定车载变风量空调系统新风、回风、总送风和风机全压的数学关系；
2）确定空调机组关键部件的模型及阻力特性；
3）确定变频风机的模型和性能参数；
4）建立整车数值计算模型；
5）设计并计算工况，进行整车的气流流动仿真计算；
6）通过变频风机的实际工作点求解变风量空调运行时的新风量、回风量和总风量；
7）利用多元非线性回归分析方法，建立变风量空调系统预测模型。
有益效果：本发明对优化轨道车辆上的空调机组设计以及提高车辆舒适性有很大的指导意义：首先，可以利用发明得到的空调机组关键部件阻力特性指导空调设计，对阻力较大的部件进行优化，降低空调机组阻力；其次，可以对不同的空调系统进行简化计算和试验验证，满足车辆舒适性对于回风量和新风量配比的要求。
另外，在轨道车辆运行过程中，变风量空调控制系统可根据变风量空调系统风量预测模型，以及室内温度、载客量、季节和新风温度的不同，调节送风机运转速度以及风阀的开度，改变送风系统的送风量以及新风量，进而提高车厢旅客的舒适性以及空调系统的节能性。
附图说明
图1是车载变风量空调系统的内部气流流动路径物理模型；
图2是风机实际工作点确定示意图；
图3是风机变风量工作点确定示意图；
图4-a为不分单元的蒸发器模型，图4-b为分单元的蒸发器模型；
图5为多单元多孔介质模型；
图6是蒸发器的物理模型；
图7是蒸发器实验值拟合方程；
图8是风阀的物理模型；
图9-1、9-2是车厢结构示意图；
图10是车载变风量空调系统整体物理模型；
图11是总风量相对偏差图；
图12是新风量相对偏差图；
图13是回风量相对偏差图。
具体实施方式
建立车载变风量空调系统风量预测模型的方法，包括以下步骤：
1）确定车载变风量空调系统新风、回风、总送风和风机全压的数学关系；
2）确定空调机组关键部件的模型及阻力特性；
3）确定变频风机的模型和性能参数；
4）建立整车数值计算模型；
5）设计并计算工况，进行整车的气流流动仿真计算；
6）通过变频风机的实际工作点求解变风量空调运行时的新风量、回风量和总风量；
7）利用多元非线性回归分析方法，建立变风量空调系统预测模型。
上述空调机组关键部件包括蒸发器、新风口、新风滤网、混合风滤网、新风阀、回风阀，各部件的模型及阻力特性公式如下：
蒸发器、新风口、滤网、新风滤网、混合风滤网用多孔介质模型，既采用Fluent软件多孔介质仿真的方法，获得速度与压降关系，确定各阻力单元的阻力系数。
阻力特性公式：
蒸发器：y = 7.3099x2 + 9.4747x
新风口：y = 5.43x2 + 3.071x
新风滤网：新风滤网的阻力特性由试验确定：y = 0.699x2 + 1.5733x
混合风滤网：混合风滤网的阻力由试验确定：y = 5.9811x2 + 13.349x
新风阀、回风阀实体建模，
阻力特性公式：
新风阀、回风阀：y=0.092x1 2-11.951x1-4.633x2 2+56.481x2+254.32
式中： x—风速（m/s），x1—风阀角度（°），x2—风速（m/s），y—阻力（Pa）。
离心风机使用Fluent风扇模型，
三个频率下离心风机的风量与压力关系如下所示：
40Hz：y = -0.000007 x2 + 0.015873 x + 389.233
50Hz：y = -0.000009 x2 + 0.016788 x + 586.381
60Hz：y = -0.000011 x2 - 0.003629 x + 804.036
式中：x—风量（m3/h），y—风机全压（Pa）。
整车数值计算模型中包括空调、风道、车厢；空调包括蒸发器、新风口、新风滤网、混合风滤网、新风阀、回风阀。
在步骤5）中，设计工况时，选取三个变量参数：新风阀角度、回风阀角度、风机频率，风阀角度选取从全开到全闭22.5°一个间隔，在风机工作范围内选取5个风机全压。
在步骤7）中，使用SPSS软件的多元非线性回归分析方法分别建立总风量和新风量的回归关系，从而获得相应的回归方程。
通过上述步骤得到的变风量空调系统预测模型为：
总风量模型：y=0.048x1 2-3.767x1-0.205x2 2+27.281x2-2.186x3 2+263.885x3-3991.739
式中：y—总风量（m3/h），x1—新风阀开度（°），x2—回风阀开度（°），x3—风机频率（Hz），40~60Hz。
新风量模型：y=-0.052x1 2+5.643x1+0.14x2 2-18.927x2-0.22x3 2+27.151x3+1147.158
式中：y—新风量（m3/h），x1—新风阀开度（°），x2—回风阀开度（°），x3—风机频率（Hz），40~60Hz。
为了验证模型的准确性，还包括模型验证步骤。
下面结合附图对本发明做进一步详细说明。
理论基础，确定车载变风量空调系统新风、回风、总送风和风机全压的数学关系：
车载变风量空调系统的内部气流流动路径物理模型如图1所示。新风通过新风入口进入空调机组后，与来自车厢顶部的回风充分混合，依次流经蒸发器、离心风机、送风道、风道、车厢后，一部分气流（回风）通过空调机组的回风入口进入机组内部，另一部分通过车顶上部空间的排风口排出。
根据上述物理模型，建立定量描述空调气流在整个轨道车辆空调机组及车厢内部流动过程中流量与压力变化的稳态数学模型。基于该模型，可获得新风量、回风量、总风量与风机全压的数学关系。
根据质量守恒定理，新风量与回风量之和为总风量；根据恒定总流伯努利方程，风机提供的全压用来克服气流在整个空调系统中（如图1所示）的流动阻力；同时，风机的总流量与全压紧密相关。
基于上述理论建立以下方程：流量平衡方程、能量（机械能）平衡方程、风机性能曲线方程。
a. 流量平衡方程 qmf=qmv+qmr。
式中，qmf：总风量，qmv：新风量，qmr：回风量。
b. 能量平衡方程 Pf=ΣΔp。
式中，P f：风机全压，ΣΔp：管路总阻力损失。
c. 风机性能曲线方程 Pf=f(qmf)。
风机实际工作点如图2所示，M点是风机的实际工作点，即风机性能曲线（表示风机全压随总流量的变化关系）与管路性能曲线（表示系统阻力随总流量的变化关系）的交点，通过该点可以获得当前工作状态下的总风量和风机全压。
图3是风机变风量工作点示意图。a是初始条件的风机性能曲线，A是初始条件的管路性能曲线，这时风机性能曲线a与管路性能曲线A的交点是1点。增加风机的频率，风机性能曲线从a改为b，这时风机性能曲线b与管路性能曲线A的交点是2点；增大阀门的开度，管路性能曲线从A改为B，这时管路性能曲线B与风机性能曲线a的交点是3点。
关键部件气流流动仿真：
空调机组蒸发腔由离心风机、蒸发器、新风阀、回风阀、滤网等关键部件以及其他一些部件组成，在这些关键部件中离心风机提供动力，蒸发腔及风道和车厢内的空气流动主要靠离心风机的压力决定；当气流流过蒸发器、风阀、滤网等部件时会产生压力损失，压力损失大小与风速以及这些部件的形状尺寸有关；因此，需要研究不同部件的阻力特性以便于了解这些部件在气流流动中所起的作用。
本发明中，空调机组关键部件包括蒸发器、新风口、新风滤网、混合风滤网、新风阀、回风阀。
针对蒸发器、新风口、新风滤网、混合风滤网，由于对其直接建模的工作量巨大，本实施例采用Fluent软件多孔介质仿真的方法，获得速度与压降关系，确定各阻力单元的阻力系数。
下面以蒸发器为例，描述建模过程：
实际蒸发器由大量翅片和铜管构成，翅片间距很小，直接建立蒸发器的完整模型难度很大，因此，选取单个翅片建立数值模型，根据气流通过蒸发器的流动规律，使用周期性边界条件，计算蒸发器在不同风速下的阻力。
蒸发器翅片密集，对空调内部流动气流的影响主要体现在两方面：（1）“导流”作用。气流通过蒸发器时主要流动方向应是平行于翅片方向；（2）“阻力”作用。气流经过蒸发器的翅片后会有明显的阻力损失。
根据蒸发器对气流流动的主要影响，采用多孔介质模型对其内部气流流动进行合理仿真。蒸发器多孔介质模型具体仿真方案如下：
（1）图4-a为不分单元的蒸发器模型，图4-b为分单元的蒸发器模型，由图可知，仿真模型图4-b中每个单元保持较小间距，可使气流在每个单元内部的流动方向与单元平行（模拟“导流”作用）；
（2） 各单元之间的壁面设置“厚度”（考虑翅片总厚度对流速的影响），使仿真模型迎风面积与蒸发器实际迎风面积相同；
（3） 确定多孔介质阻力特性（模拟“阻力”作用）。需要的关键参数基于变风量风道实验结果获得，具体求解方法见本说明书下面描述。
根据上述方案步骤，建立更为接近真实情况的蒸发器气流流动仿真模型，即“多单元多孔介质模型”如图5所示。
将“多单元多孔介质”仿真模型应用于蒸发器气流流动仿真，建立物理模型，如图6所示。
Fluent仿真计算时，蒸发器物理模型按多孔介质边界条件作简化处理，多孔介质生成的压力梯度与速度（速度平方和）成正比，即：
Δp/L=1/2*ρC2v2+μv *1/k
其中，1/k：粘性阻力系数，1/m2；
C2：惯性阻力系数，1/m；
ρ：密度，kg/m3；
μ：动力粘度，kg/(ms)；
Δp：压降，Pa；
L：长度，m；
v：速度，m/s。
由蒸发器的风道平台实验值求解速度与压降的二次曲线，蒸发器的速度与压降关系及拟合方程如图7所示。
将实验值拟合方程得到的二次多项式系数与上述公式的系数对应，解得：
k=2.5604×10-7 m2
C2=101.13 1/m
将获得的阻力系数代入“多单元多孔介质模型”进行仿真计算，仿真值与实验值的偏差见下表。
表1 仿真值与实验值偏差比较
| 速度(m/s) | 仿真值压降（Pa） | 实验值压降（Pa） | 偏差（%） |
| 1.0055 | 19.88 | 17.5 | 13.62 |
| 1.2569 | 27.9 | 23.5 | 18.73 |
| 1.5083 | 37.14 | 32.5 | 14.26 |
| 1.6759 | 43.96 | 38.0 | 15.7 |
| 1.8435 | 51.33 | 43.9 | 16.93 |
从上表可知，由实验值获得的速度压降关系直接代入“多单元多孔介质模型”会有比较大的偏差，采取将实验值压降乘以校正系数的方法，重新获得速度压降关系再代入蒸发器多孔介质模型，校正后的仿真值与实验值的偏差见下表。
表2校正后仿真值与实验值偏差比较
| 速度(m/s) | 仿真值压降（Pa） | 实验值压降（Pa） | 偏差（%） |
| 1.0055 | 16.86 | 17.5 | -3.66 |
| 1.2569 | 23.66 | 23.5 | 0.68 |
| 1.5083 | 31.49 | 32.5 | -3.1 |
| 1.6759 | 37.29 | 38.0 | -1.87 |
| 1.8435 | 43.54 | 43.9 | -0.82 |
由上表可知仿真结果与实验结果最大相对偏差仅为-3.66%，表明建立的仿真模型是合理的。
最后确定蒸发器的阻力特性公式：y = 7.3099x2 + 9.4747x，x—风速（m/s），y—阻力（Pa）。
新风口、新风滤网、混合风滤网阻力单元采用与蒸发器类似的仿真方法，应用多孔介质边界条件获得每个部件的阻力特性。阻力特性公式：
新风口：y = 5.43x2 + 3.071x
新风滤网：新风滤网的阻力特性由试验确定：y = 0.699x2 + 1.5733x
混合风滤网：混合风滤网的阻力由试验确定：y = 5.9811x2 + 13.349x
式中：x—风速（m/s），y—阻力（Pa）。
新风阀、回风阀阻力单元采用实体建模的方法进行Fluent仿真计算，由于风阀的实际结构比较复杂，物理建模时将其简化为平面的形状，如图8所示，风阀偏转方向与风道实验时相同。
风阀的阻力与风速和阀门角度有关，其阻力特性公式如下：
y=0.092x1 2-11.951x1-4.633x2 2+56.481x2+254.32
式中： x1—风阀角度（°），x2—风速（m/s），y—阻力（Pa）。
确定变频风机的模型和性能参数：
对离心风机进行性能试验检测，可以获得风机的全压值和流量值，当风机的频率改变时，风机全压值和流量值的相应关系也会改变，即不同频率下有不同的风机性能曲线。
三个频率下离心风机的风量与压力关系如下所示：
40Hz：y = -0.000007 x2 + 0.015873 x + 389.233
50Hz：y = -0.000009 x2 + 0.016788 x + 586.381
60Hz：y = -0.000011 x2 - 0.003629 x + 804.036
式中：x—风量（m3/h），y—风机全压（Pa）。
建立整车数值计算模型：
本实施例中，变风量空调机组采用前出下回型机组。车厢内部空间尺寸为：长约32.2m（含司机室），宽约2.4m，高约2.2m；送风主风道总长28252mm（不包括司机室送风软风道长度），主风道总宽1500mm，主风道总高142.6mm。新风通过新风入口进入空调机组后，与来自车厢顶部的回风充分混合，依次流经蒸发器、离心风机、送风道、风道、车厢后，一部分气流（回风）通过空调机组的回风入口进入机组内部，另一部分通过车顶上部空间的排风口排出。车厢结构如图9-1、9-2所示。
物理模型：由于车厢内部结构复杂，车内障碍物较多，完全按照实际结构建立物理模型存在很大的困难，同时生成网格时网格数量巨大，对计算机硬件要求较高，并且仿真计算周期很长，所以在建立车厢物理模型时，对车厢结构做出合理简化。
车厢、风道和变风量空调机组沿长度、宽度方向对称，仿真中取四分之一为研究对象，分别以坐标轴的X、Y、Z轴作为模型的长、宽、高方向，按照实际尺寸建模，在仿真模拟时采用对称模拟方法，在长度方向和宽度方向设置对称面。车载变风量空调系统整体物理模型如图10所示。
网格划分：在划分网格时，根据物理模型结构特点，采用结构和非结构两种类型单元生成网格，且风口网格局部加密，以提高计算精度。车载变风量空调系统整体网格质量（等角倾斜率分布，Equisize Skew）如下表所示，车载变风量空调系统整车模型依次按长、宽、高方向进行网格划分。
表3车载变风量空调系统整体网格质量
| 起始值 | 终值 | 数量 | 百分比（%） |
| 0 | 0.1 | 2157825 | 55.17 |
| 0.1 | 0.2 | 1105874 | 28.28 |
| 0.2 | 0.3 | 425188 | 10.87 |
| 0.3 | 0.4 | 172027 | 4.4 |
| 0.4 | 0.5 | 42079 | 1.08 |
| 0.5 | 0.6 | 7845 | 0.2 |
| 0.6 | 0.7 | 155 | 0.0 |
| 0.7 | 0.8 | 0 | 0.0 |
| 0.8 | 0.9 | 0 | 0.0 |
| 0.9 | 1 | 0 | 0.0 |
| 总计 | — | 3910993 | 100 |
注：等角倾斜率最大值为0.68
由上表可知，绝大部分网格单元的网格质量高，99.8%的网格单元网格质量在0.5以下，最大网格质量为0.68。
边界条件：车载变风量空调系统整车模型均采用以下计算边界条件及其计算方法：
（1）以新风入口为计算入口边界，设为压力入口，与外界大气相连，压力设为0Pa，Fluent参数设置面板中的水力直径计算方法如下：水力直径：Dh=4s/c（m），其中：s为湿截面面积（m2），c为湿周周长（m）。
（2）以废排出口为计算出口边界，设为压力出口，与外界大气相连，压力设为0Pa。
（3）空调机组、风道、排风末端有内部风机，Fluent参数设置面板中设为fan边界条件，需要设置风机压头值。
仿真模型中采用的边界条件如下表所示。
表4仿真模型边界条件
| 编号 | 边界类型 | 水利直径(m) | 湍流强度（%） | 表压力（Pa） |
| Inlet | Pressure-inlet | 0.35 | 5 | 0 |
| Outlet1 | Pressure-outlet | 0.03 | 5 | 0 |
| Outlet2 | Pressure-outlet | 0.03 | 5 | 0 |
| Outlet3 | Pressure-outlet | 0.03 | 5 | 0 |
| Outlet4 | Pressure-outlet | 0.03 | 5 | 0 |
| Outlet5 | Pressure-outlet | 0.03 | 5 | 0 |
| Outlet6 | Pressure-outlet | 0.03 | 5 | 0 |
| Outlet7 | Pressure-outlet | 0.03 | 5 | 0 |
| Outlet8 | Pressure-outlet | 0.03 | 5 | 0 |
| Outlet9 | Pressure-outlet | 0.03 | 5 | 0 |
| Outlet10 | Pressure-outlet | 0.03 | 5 | 0 |
| Outlet11 | Pressure-outlet | 0.03 | 5 | 0 |
| Outlet12 | Pressure-outlet | 0.03 | 5 | 0 |
| Outlet13 | Pressure-outlet | 0.03 | 5 | 0 |
| Outlet14 | Pressure-outlet | 0.03 | 5 | 0 |
| Outlet15 | Pressure-outlet | 0.03 | 5 | 0 |
| Outlet16 | Pressure-outlet | 0.03 | 5 | 0 |
| Outlet17 | Pressure-outlet | 0.03 | 5 | 0 |
| Outlet18 | Pressure-outlet | 0.03 | 5 | 0 |
| Outlet19 | Pressure-outlet | 0.03 | 5 | 0 |
| Outlet20 | Pressure-outlet | 0.03 | 5 | 0 |
| Outlet21 | Pressure-outlet | 0.03 | 5 | 0 |
| Outlet22 | Pressure-outlet | 0.03 | 5 | 0 |
| Outlet23 | Pressure-outlet | 0.03 | 5 | 0 |
| Outlet24 | Pressure-outlet | 0.03 | 5 | 0 |
| Outlet25 | Pressure-outlet | 0.03 | 5 | 0 |
| Outlet26 | Pressure-outlet | 0.03 | 5 | 0 |
| Outlet27 | Pressure-outlet | 0.03 | 5 | 0 |
| Outlet28 | Pressure-outlet | 0.03 | 5 | 0 |
| Outlet29 | Pressure-outlet | 0.03 | 5 | 0 |
| Outlet30 | Pressure-outlet | 0.03 | 5 | 0 |
仿真计算：将物理问题描述为数学问题的过程中，首先进行简化假设，进而建立相关物理变量的微分方程，最后在相应的定解条件下求解微分方程。以下是建立车载变风量空调系统数学模型而作的相关简化：
（1）空调内部气流流动为充分发展的三维稳态湍流流动；
（2）忽略重力影响；
（3）可将各个壁面视为绝热，计算时不考虑空调气流与外界的换热。
算法选择及其设置：本实施例采用SIMPLEC算法，求解压力与速度耦合方程。
仿真结果：车载变风量空调系统整车模型建立后，开始仿真计算。当新风阀与回风阀开启角度改变时，整个城轨车辆空调系统的新风量与回风量会相应改变，当风机全压改变时，整个城轨车辆空调系统的总风量也会相应改变，因此需要确定新风阀、回风阀开启角度及风机全压。
设计并计算工况，进行整车的气流流动仿真计算：
考虑到计算工况数，计算时把风阀角度从全开到全闭平均分为4个间隔，即22.5°一个间隔；同时，在风机性能曲线上在风机工作范围内选取5个风机全压，以上三个参数进行排列组合，共计4*4*5=80组工况，如下表所示。本实施例中，当计算到50组工况时，拟合的多元非线性回归关系已经非常明显，R2=0.997，因此剩余的30组工况未计算，共计算了50组工况。
表5新回风阀角度及风机全压
| 新风阀角度（°） | 22.5 | 45 | 67.5 | 90 | — |
| 回风阀角度（°） | 22.5 | 45 | 67.5 | 90 | — |
| 风机全压（Pa） | 304 | 372 | 497 | 541 | 576 |
通过变频风机的实际工作点求解变风量空调运行时的新风量、回风量和总风量：
1、对车载变风量空调系统整车模型进行迭代计算，将新风阀、回风阀调节至每一种设计角度，风机全压调整到每一种设计值时，就会得到一组仿真结果。
2、确定总风量：
（1）管路性能曲线方程: ∆p=a1qmf 2+b1qmf ，∆p：管路总阻力损失，qmf：总风量。
在车载变风量空调系统中，风机提供的全压用来克服气流在整个空调系统中的流动阻力，如果整个系统没有气流，即总流量为零，管路总阻力损失也为零，因此管路性能曲线应通过坐标原点。
对车载变风量空调系统整车模型进行迭代计算，可以获得新风阀与回风阀在设计角度时的一组全压值和流量值，采用拟合的方法确定管路性能曲线方程，从而得到式中参数a1、b1的数值。
（2）风机性能曲线方程：∆p=a2qmf 2+b2qmf+c1
对离心风机进行性能试验检测，可以获得不同频率下风机的全压值和流量值，同样能采用拟合的方法确定风机性能曲线方程，从而得到式a2、b2、c1的数值。风机性能曲线具体参数由设备方提供。
三个频率下离心风机的风量与压力关系如下所示：
40Hz：y = -0.000007 x2 + 0.015873 x + 389.233
50Hz：y = -0.000009 x2 + 0.016788 x + 586.381
60Hz：y = -0.000011 x2 - 0.003629 x + 804.036
式中：x—风量（m3/h），y—风机全压（Pa）。
基于上述方程（管路性能曲线方程、风机性能曲线方程），风机性能曲线与管路性能曲线的交点是风机实际工作点，从该点可知车载变风量空调系统整车模型的总风量和风机压头。求解得：
qmf ={(b2-b1)±sqrt[(a1-a2)2-4(a1-a2)(-c1)]}/[2(a1-a2)]
∆p=a1qmf 2+b1qmf
3、确定新风量与回风量
新风量与回风量可以通过新回风比（回风量/新风量）确定，而新回风比会根据总风量、新风阀开度和回风阀开度的变化而变化，函数关系表示为：
新回风比=f（总风量、新风阀开度、回风阀开度）。
因此在车载变风量空调系统整车模型中，当新风阀和回风阀开度确定时，新回风比只是总风量的单值函数，通过建立这两者的回归关系，从而得到新风量和回风量。
n=f(qmf，新风阀开度，回风阀开度)
qmf= qmv+ qmr
n= qmr/qmv
qmf：总风量，qmv：新风量，qmr：回风量，n：新回风比。
求解得：
qmv=qmf/(1+n)，qmr=qmf-qmv。
通过上述过程计算的结果如下表：
表6 变风量运行实际工作点
车载变风量空调系统整车模型的全部计算工况的变风量工作点确定后，对所有计算结果进行仿真验证。验证方法为在Fluent整车仿真模型中，将每种计算工况调整到对应的计算风机压头值，比较总风量、新风量和回风量的相对偏差。部分计算结果与仿真结果偏差比较如下表（风机工作在40Hz下）：
表7 40Hz仿真结果与计算结果比较
由上表可知所有结果的相对偏差都在2%以下，说明变风量数学求解的计算值与Fluent迭代计算的仿真值的吻合程度较好。
利用多元非线性回归分析方法，建立变风量空调系统预测模型：
将整车模型变风量工作点的总风量、新风量和回风量进行数学求解后，为了得到更加通用的模型，需要建立各变量间的函数关系，即总风量=f（新风阀开度、回风阀开度、风机频率），新风量=f（新风阀开度、回风阀开度、风机频率），回风量为总风量与新风量之差。这样在已知任何新风阀开度、回风阀开度和风机频率时，通过函数求解能够获得对应的总风量、新风量和回风量。
本发明利用SPSS软件的多元非线性回归分析方法分别建立总风量和新风量的回归关系，从而获得相应的回归方程。
总风量模型如下所示：
y=0.048x1 2-3.767x1-0.205x2 2+27.281x2-2.186x3 2+263.885x3-3991.739。
式中：y—总风量（m3/h），x1—新风阀开度（°），x2—回风阀开度（°），x3—风机频率（Hz），40~60Hz。
下表为总风量回归模型的假设检验，从表中可知R2为0.997，能较好地解释总风量与新风阀开度、回风阀开度及风机频率间的多元非线性回归关系。
表8 总风量回归模型的假设检验
| 源 | 平方和 | df | 均方 |
| 回归 | 5.589E8 | 7 | 79838855.670 |
| 残差 | 14805.174 | 23 | 643.703 |
| 未更正的总计 | 5.589E8 | 30 | |
| 已更正的总计 | 5284729.030 | 29 |
a. R2 = 1 -（残差平方和）/（已更正的平方和）= .997。
新风量模型如下所示：
y=-0.052x1 2+5.643x1+0.14x2 2-18.927x2-0.22x3 2+27.151x3+1147.158。
式中：y—总风量（m3/h），x1—新风阀开度（°），x2—回风阀开度（°），x3—风机频率（Hz），40~60Hz。
下表为新风量回归模型的假设检验，从表中可知R2为0.958，能较好地解释新风量与新风阀开度、回风阀开度及风机频率间的多元非线性回归关系。
表9 新风量回归模型的假设检验
| 源 | 平方和 | df | 均方 |
| 回归 | 71633937.723 | 7 | 10233419.675 |
| 残差 | 21470.664 | 23 | 933.507 |
| 未更正的总计 | 71655408.386 | 30 | |
| 已更正的总计 | 5284729.030 | 29 |
a. R2 = 1 -（残差平方和）/（已更正的平方和）= .998。
模型验证：车载变风量空调系统整车模型的预测模型确定后，可以求解所有计算工况的总风量、新风量和回风量。预测模型结果与变风量计算结果偏差比较见下面三个表。
表10 40Hz预测结果与计算结果比较
表11 50Hz预测结果与计算结果比较
表12 60Hz预测结果与计算结果比较
由上表可知总风量、新风量及回风量的相对偏差都在10%以下，说明利用SPSS软件建立的总风量多元非线性回归方程、新风量多元非线性回归方程是可以接受的，能够利用此预测模型求解变风量风机的实际工作点。
实验验证：首先将新风阀、回风阀调整到90度，风机50Hz，送风口加孔板、回风口加胶带模拟风道和车厢阻力，调整孔板和胶带尺寸使试验的总风量值达到该工况的计算值。然后调整新风阀和回风阀角度及风机频率，分别记录各个工况下的新风量和回风量，并与Fluent计算值比较。
由实验结果（如图11、12、13所示）可知，总风量偏差都在10%以内，回风量偏差都在20%以内。新风量有少量点大于20%偏差，主要位于新风阀22.5°的工况，分析原因如下：
1、试验时没有废排装置，仅依靠在回风口粘贴胶带施加阻力的方法，不能保证新风阀小角度时的进风；
2、建立模型时对风阀结构进行了一些简化，导致与试验结果存在差异。
Claims (6)
1.建立车载变风量空调系统风量预测模型的方法，其特征在于包括以下步骤：
1）确定车载变风量空调系统新风、回风、总送风和风机全压的数学关系；
2）确定空调机组关键部件的模型及阻力特性；
3）确定变频风机的模型和性能参数；
4）建立整车数值计算模型；
5）设计并计算工况，进行整车的气流流动仿真计算；
6）通过变频风机的实际工作点求解变风量空调运行时的新风量、回风量和总风量；
7）利用多元非线性回归分析方法，建立变风量空调系统预测模型；
步骤2）中，空调机组关键部件包括蒸发器、新风口、新风滤网、混合风滤网、新风阀、回风阀，各部件的模型及阻力特性公式如下：
蒸发器、新风口、滤网、新风滤网、混合风滤网采用多孔介质模型，新风阀、回风阀实体建模，
阻力特性公式：
蒸发器：y = 7.3099x2 + 9.4747x
新风口：y = 5.43x2 + 3.071x
新风滤网：y = 0.699x2 + 1.5733x
混合风滤网：y = 5.9811x2 + 13.349x
新风阀、回风阀：y=0.092x1 2-11.951x1-4.633x2 2+56.481x2+254.32
式中： x为风速，单位是m/s，x1为风阀角度，单位是°，x2为风速，单位是m/s，y为阻力，单位是Pa；
步骤7）中得到的变风量空调系统预测模型为：
总风量模型：y=0.048x1 2-3.767x1-0.205x2 2+27.281x2-2.186x3 2+263.885x3-3991.739
式中：y为总风量，单位是m3/h，x1为新风阀开度，单位是°，x2为回风阀开度，单位是°，x3为风机频率，单位是Hz，取值范围：40~60Hz；
新风量模型：y=-0.052x1 2+5.643x1+0.14x2 2-18.927x2-0.22x3 2+27.151x3+1147.158
式中：y为新风量，单位是m3/h，x1为新风阀开度，单位是°，x2为回风阀开度，单位是°，x3为风机频率，单位是Hz，取值范围：40~60Hz。
2.根据权利要求1所述的建立车载变风量空调系统风量预测模型的方法，其特征在于在步骤3）中，离心风机使用Fluent风扇模型，
三个频率下离心风机的风量与压力关系如下所示：
40Hz：y = -0.000007 x2 + 0.015873 x + 389.233
50Hz：y = -0.000009 x2 + 0.016788 x + 586.381
60Hz：y = -0.000011 x2 - 0.003629 x + 804.036
式中：x为风量，单位是m3/h，y为风机全压，单位是Pa。
3.根据权利要求1所述的建立车载变风量空调系统风量预测模型的方法，其特征在于在步骤4）中，整车数值计算模型中包括空调、风道、车厢；空调包括蒸发器、新风口、新风滤网、混合风滤网、新风阀、回风阀。
4.根据权利要求1所述的建立车载变风量空调系统风量预测模型的方法，其特征在于在步骤5）中，设计工况时，选取三个变量参数：新风阀角度、回风阀角度、风机频率，风阀角度选取从全开到全闭22.5°一个间隔，在风机工作范围内选取5个风机全压。
5.根据权利要求1所述的建立车载变风量空调系统风量预测模型的方法，其特征在于在步骤7）中，使用SPSS软件的多元非线性回归分析方法分别建立总风量和新风量的回归关系，从而获得相应的回归方程。
6.根据权利要求1所述的建立车载变风量空调系统风量预测模型的方法，其特征在于还包括模型验证步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611198508.2A CN106777711B (zh) | 2016-12-22 | 2016-12-22 | 建立车载变风量空调系统风量预测模型的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611198508.2A CN106777711B (zh) | 2016-12-22 | 2016-12-22 | 建立车载变风量空调系统风量预测模型的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106777711A CN106777711A (zh) | 2017-05-31 |
| CN106777711B true CN106777711B (zh) | 2019-09-17 |
Family
ID=58899409
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611198508.2A Active CN106777711B (zh) | 2016-12-22 | 2016-12-22 | 建立车载变风量空调系统风量预测模型的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106777711B (zh) |
Families Citing this family (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109059217B (zh) * | 2018-05-04 | 2021-09-28 | 中国中元国际工程有限公司 | 一种基于运行曲线的变风量空调系统总风量控制方法 |
| CN108763721B (zh) * | 2018-05-23 | 2022-09-30 | 特灵空调系统(中国)有限公司 | 空调系统充注量的仿真方法 |
| CN109241634B (zh) * | 2018-09-12 | 2020-08-21 | 珠海格力电器股份有限公司 | 空调吹风广角确定方法、装置以及空调 |
| CN109323427B (zh) * | 2018-11-02 | 2023-08-25 | 青岛海尔空调器有限总公司 | 一种空调器室内机的控制方法、装置、空调器及存储介质 |
| CN110569515A (zh) * | 2018-12-24 | 2019-12-13 | 中国电力科学研究院有限公司 | 用于计算综合管廊风机压损的分析方法 |
| CN109798646B (zh) * | 2019-01-31 | 2021-03-30 | 上海真聂思楼宇科技有限公司 | 一种基于大数据平台的变风量空调控制系统和方法 |
| CN111710034B (zh) * | 2020-06-24 | 2022-06-21 | 四川长虹空调有限公司 | 空调器室外机百叶窗内出风量仿真方法及系统 |
| CN111859662B (zh) * | 2020-07-17 | 2023-09-22 | 泰州市都瑞堡船舶设备有限公司 | 一种变风量空调末端单元内均速管流量计的布设方法 |
| CN112052997B (zh) * | 2020-09-07 | 2021-09-03 | 西安建筑科技大学 | 一种变风量空调系统事件触发预测控制方法 |
| CN112444410A (zh) * | 2020-10-27 | 2021-03-05 | 北京新能源汽车技术创新中心有限公司 | 一种车辆鼓风机的风量测试方法 |
| CN113264075B (zh) * | 2021-06-24 | 2023-04-07 | 中车株洲电力机车有限公司 | 空调机组、轨道车辆空调系统及其风量控制方法 |
| CN113486953A (zh) * | 2021-07-06 | 2021-10-08 | 西门子（上海）电气传动设备有限公司 | 变频器滤网更换时间的预测方法、装置及计算机可读介质 |
| CN114414833A (zh) * | 2021-12-30 | 2022-04-29 | 苏州悠远环境科技有限公司 | 一种基于压差法监测ffu出风面风速的方法 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1224823A (zh) * | 1997-11-07 | 1999-08-04 | 兰迪斯及斯特法有限公司 | 具有前馈和反馈控制的室温控制装置及方法 |
| EP1887673A1 (en) * | 2006-08-11 | 2008-02-13 | Abb Research Ltd. | Parameter estimation of a thermal model of a power line |
| CN102353119A (zh) * | 2011-08-09 | 2012-02-15 | 北京建筑工程学院 | 一种vav变风量空调系统控制方法 |
| CN104359160A (zh) * | 2014-11-07 | 2015-02-18 | 深圳海悟科技有限公司 | 一种节能空调器 |
| CN105216813A (zh) * | 2015-09-08 | 2016-01-06 | 石家庄国祥运输设备有限公司 | 高速列车用空调新风阀精确调控方法 |
| CN106250621A (zh) * | 2016-07-29 | 2016-12-21 | 西南交通大学 | 列车空调能耗计算系统 |
-
2016
- 2016-12-22 CN CN201611198508.2A patent/CN106777711B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1224823A (zh) * | 1997-11-07 | 1999-08-04 | 兰迪斯及斯特法有限公司 | 具有前馈和反馈控制的室温控制装置及方法 |
| EP1887673A1 (en) * | 2006-08-11 | 2008-02-13 | Abb Research Ltd. | Parameter estimation of a thermal model of a power line |
| CN102353119A (zh) * | 2011-08-09 | 2012-02-15 | 北京建筑工程学院 | 一种vav变风量空调系统控制方法 |
| CN104359160A (zh) * | 2014-11-07 | 2015-02-18 | 深圳海悟科技有限公司 | 一种节能空调器 |
| CN105216813A (zh) * | 2015-09-08 | 2016-01-06 | 石家庄国祥运输设备有限公司 | 高速列车用空调新风阀精确调控方法 |
| CN106250621A (zh) * | 2016-07-29 | 2016-12-21 | 西南交通大学 | 列车空调能耗计算系统 |
Non-Patent Citations (1)
| Title |
|---|
| 基于RWI指标的地铁空调变风量分析;张圆圆;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》;20120415;第C033-159/摘要、正文18-44页 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106777711A (zh) | 2017-05-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106777711B (zh) | 建立车载变风量空调系统风量预测模型的方法 | |
| Tao et al. | Numerical and experimental study on ventilation panel models in a subway passenger compartment | |
| Moghaddam et al. | Small-scale single-panel liquid-to-air membrane energy exchanger (LAMEE) test facility development, commissioning and evaluating the steady-state performance | |
| Calautit et al. | CFD and wind tunnel study of the performance of a uni-directional wind catcher with heat transfer devices | |
| CN105115738A (zh) | 顶置式空调系统的性能测试系统及测试方法 | |
| CN102102897B (zh) | 分离式空调设备的温度控制方法 | |
| Guan et al. | Geometric optimization on active chilled beam terminal unit to achieve high entrainment efficiency | |
| Wang et al. | Prediction of vertical thermal stratification of large space buildings based on Block-Gebhart model: Case studies of three typical hybrid ventilation scenarios | |
| CN107023966A (zh) | 一种地铁站空调冷却水出水温度设定值优化方法 | |
| Yang et al. | Research on flow resistance characteristics with different window/door opening angles | |
| CN206974675U (zh) | 新能源汽车热管理系统开发匹配试验台 | |
| CN105938058A (zh) | 一种复合测量车载空调及电池热管理性能的焓差试验装置 | |
| CN109114998A (zh) | 一种机械通风逆流式冷却塔消雾改造的设计计算方法 | |
| CN105973582A (zh) | 一种空调风道高低温载流疲劳试验装置 | |
| Gao et al. | High-efficiency diffuser based on a normalized evaluation index of jet length and resistance | |
| CN103217310B (zh) | 一种食品冷柜性能测试装置 | |
| He et al. | Control strategy analysis of multistage speed compressor for vehicle air conditioning based on particle swarm optimization | |
| CN203083830U (zh) | 整车空调系统模拟试验装置 | |
| CN108756996B (zh) | 一种矿井热湿环境模拟实验平台 | |
| CN204988748U (zh) | 顶置式空调系统的性能测试系统 | |
| CN104008300B (zh) | 一种地下空间通风量算法 | |
| CN203572666U (zh) | 一种冷链冷藏设备试验系统 | |
| CN105956316A (zh) | 一种汽车发动机舱前端模块空气进风量计算方法 | |
| Vali | Modeling a run-around heat and moisture exchanger using two counter/cross flow exchangers | |
| Tao et al. | Comparison on simulation and experiment of supply air through metro vehicle air conditioning duct |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |