基于多维曲线拟合算法的空调装置制冷量测试方法

摘要

本发明公开了基于多维曲线拟合算法的空调装置制冷量测试方法，属于空调制冷量测试方法技术领域。测试方法依据的原理为：建立空调装置的基准样机物理模型，以此为研究对象，通过建立空调装置基准样机的制冷量与各个影响参数之间的多维曲线拟合数学模型，采用包括名义工况在内的多组样本数据进行训练仿真，获得该空调装置在任意工况下的制冷量数学模型。通过该数学模型，可以获得任意工况下空调装置基准样机的制冷量，作为同类型被测空调装置制冷量的理论值。在相同测试工况下，获得被测空调装置实际制冷量。利用该测试方法制作的制冷量测试装置或仪器可以用于空调装置的制冷量测试，广泛用于交通车辆、船舶、民用与工业应用的空调装置制冷量测试。

说明书

技术领域：

本发明涉及基于多维曲线拟合算法的空调装置制冷量测试方法，属于空调制冷量测试方法技术领域。

背景技术：

空调装置制冷量测试主要依照相关标准，按照焓差法原理进行。这种依照相关标准的测试方法离不开复杂的工况控制和风量控制等硬件设备，几乎被局限于实验室应用。但是，对于空调装置制冷量测试的需求，是在空调装置的应用过程中普遍存在的。例如，车载空调装置的应用极为广泛，但因故障率远高于普通的家用空调器，所以，需要经常离线或在线测试其制冷量，而这种对制冷量的测试往往需要在非标准工况下进行的。

由于空调装置制冷量的测试属于多参数综合测试，限于目前测试原理没有突破，所以，至今尚未脱离相关标准硬件设备条件的空调装置制冷量测试理论及其相应的测试仪器。

以下是传统空调试验方法与相关设施配置要求，《TB/T2432-93铁路客车单元式空调机组试验方法》对客车机组制冷量测试的通用方法及标准装置做了相关规定，图1为标准中规定的测试工况。机组制冷量测试常用的方法有风管热平衡法、空气焓差法和房间热平衡法。其中，空气焓差法是应用最为广泛并被空调行业普遍认可的一种测试方法。

标准中规定的空调装置性能试验台有风洞式、环路式及房间式三种，其中的环路式空调装置试验台结构原理如图2所示，整个装置分为相互独立的室内侧和室外侧实验室，并通过房间空调器维持稳定的内部环境以满足检测时的工况要求。检测过程中，通过安装在蒸发器一侧的送、回风风管上的温湿度及压差测量装置与传感器，测得规定工况下蒸发器送、回风的焓值及流量，再通过下式即可确定机组的制冷量。

式中：Q为蒸发器制冷量，W；

Lo为室内侧空气流量，m³/s；

h_n为室内送风空气焓值，J/kg干空气；

h_c为室内回风空气焓值，J/kg干空气；

V为室内侧空气比容，m³/kg；

为室内侧空气含湿量，kg/kg干空气。

传统的空调装置性能检测方法必须在特定的实验室内按照规定的工况进行，其它应用条件下几乎无法完成测试。除了测试环境的严苛外，在实际操作中，外界环境的扰动会导致试验工况的不稳定，甚至造成检测室内的空气温湿度值偏离了规范的允许值，故检测中需对试验工况进行时时监控和调整，这种工况调整中消耗大量的时间及能源。因此，传统的测试方法是不能脱离实验室专门测试系统使用的。

发明内容：

针对上述问题，本发明要解决的技术问题是：1、建立空调装置的基准样机物理模型，并以此为研究对象，建立空调装置基准样机的制冷量与各个影响参数之间的多维曲线的数学模型。2、通过该数学模型，提供一种新型空调装置制冷量测试方法。

本发明基于多维曲线拟合算法的空调装置制冷量测试方法依据的基本原理为：建立空调装置的基准样机物理模型，以此为研究对象，通过建立空调装置基准样机的制冷量与各个影响参数之间的多维曲线拟合数学模型，并采用包括名义工况在内的多组样本数据进行训练仿真，获得该空调装置在任意工况下的制冷量数学模型，通过该数学模型，可以获得任意工况下空调装置基准样机的制冷量，作为同类型被测空调装置制冷量的理论参考值，在相同测试工况下，利用该模型结果与实际测试数据的分析比较，获得被测空调装置实际制冷量。

它的具体测试方法表述为：a：把空调装置视为热力系统，该系统的输入参数包括：蒸发器回风温度t_h、蒸发器回风相对湿度蒸发器通风量V_l、冷凝器进风温度t_c、冷凝器通风量V_c及大气压力P_a,输出参数为制冷量和功耗。以下将建立空调装置制冷量与各个输入参数的多维曲线拟合模型；

由于决定空调装置制冷量的主要参数包括：蒸发器回风温度t_h、蒸发器回风相对湿度蒸发器通风量V_l、冷凝器进风温度t_c、冷凝器通风量V_c及大气压力P_a,因此对于指定的一台空调装置来说，只要确定上述六个参数，机组的能耗及制冷量也是唯一确定的，其中蒸发器回风温度t_h、蒸发器回风湿度及冷凝器送风温度t_c对机组制冷量的影响最为明显，且由于大气压力变化较小，其变化对空调装置制冷量的影响可以忽略，蒸发器、冷凝器通风量在测量精度要求的范围内可以作为常量处理，即认为在实际测试过程中空调装置冷凝器、蒸发器的通风量及大气压力基本维持不变,因此，制冷量在测量过程中的主要变量是：蒸发器回风温度t_h、蒸发器回风湿度及冷凝器进风温度t_c；

上述变量中的蒸发器回风相对湿度主要通过蒸发器凝结水影响翅片管蒸发器的管外壁对流换热系数，从而影响蒸发器的传热量。根据传热学理论和湿空气热质交换理论，当蒸发器回风相对湿度由低至高达到一定值蒸发器表面将出现干工况至湿工况的转变；因此，蒸发器回风湿度的影响表现在两个区域：干工况区和湿工况区。在这两个区内，蒸发器回风湿度的影响相对不显著；故可以对相对湿度的影响进行分段研究；在相对湿度各分段中，空调装置制冷量的主要影响因素可减少为蒸发器回风温度t_h及冷凝器进风温度t_c两个，实现了空调装置制冷量模型降维的目的。将蒸发器回风相对湿度按照临界值分区后，即可将模型降为不同回风相对湿度区段的三维模型；

临界相对湿度的数学模型的建立及干湿分区的判别：

在蒸发器表面温度10～11℃、大气压力101325Pa时，借助Matlab软件分析拟合，得到蒸发器回风温度t_h与相应干、湿分区的临界相对湿度数学模型(公式1)，模型的适用范围：10℃≤t_h≤50℃；

式中：

p₁＝6.2029e-009,p₂＝-1.5991e-006,p₃＝0.00018189

p₄＝-0.012068,p₅＝0.51418,p₆＝-14,p₇＝199.02

在实际工况下，当测得的相对湿度大于时，认定该工况点位于湿工况区，当测得的相对湿度小于时，认定该工况点位于干工况区；

空调装置制冷量多维曲线模型的建立：

利用Matlab软件，对两个相对湿度分区中制冷量与蒸发器回风温度t_h、冷凝进风温度t_c的三维曲线进行拟合，下面以某种单元式空调机组为例，给出具体拟合结果:

该机组在大气压力B(101325Pa)、送风量6000m³/h、冷凝风量11800m³/h的条件下，在蒸发器表面温度10～11℃时，得到制冷量与蒸发器回风温度、冷凝进风温度的关系如下：

Q_o＝p₀₀+p₁₀t_c+p₀₁t_h+p₂₀t_c²+p₁₁t_ct_h+p₀₂t_h²+p₃₀t_c³+p₂₁t_c²t_h+p₁₂t_ct_h²+p₀₃t_h³(公式2)

湿工况区:

p₀₀＝-26.9,p₁₀＝9.023,p₀₁＝-6.165,p₂₀＝-0.2055,p₁₁＝-0.1919p₀₂＝0.4461,p₃₀＝0.001953,p₂₁＝0.0003906,p₁₂＝0.002734p₀₃＝-0.007161

干工况区:

p₀₀＝-16.91,p₁₀＝6.02,p₀₁＝-4.223,p₂₀＝-0.1047，p₁₁＝-0.2256，p₀₂＝0.4188，p₃₀＝0.001302，p₂₁＝-0.0007813，p₁₂＝0.004687，p₀₃＝-0.007813

b：采用包括名义工况在内的多组样本数据进行训练仿真，获得该空调装置在任意工况下的制冷量数学模型，通过该数学模型，可以获得任意工况下空调装置基准样机的制冷量，并将其作为同类型被测空调装置制冷量的理论参照值，在相同测试工况下，利用该模型结果与实际测试数据的分析比较，获得被测空调装置实际制冷量；

测试工况下(任意运行工况)被测空调装置的实际制冷量Q_o,s：

Q_o,s＝V_l·ρ(h_h,j-h_h,c)>

测试工况下(任意运行工况)被测空调装置的理论制冷量Q_o,t:

则两式相除，得：

即

实际测试工况下被测空调装置的理论制冷量Q_o,t由上述制冷量数学模型(见公式2)可求出，又由公式4知：

$h_{h,j}^0-h_{h,c}^0=\frac{Q_{o,t}}{V_l·\rho }$

上式右侧均为已知参数，则理论焓降可以求得，此外，公式5右侧实际焓降h_h,j-h_h,c可在测试过程中计算得到，因此，利用公式5可以直接求得被测空调装置的实际制冷量；

在实际测量过程中测得被测空调装置的运行参数，通过上述测量空调装置制冷量的新方法，获得实际测试工况下空调装置的实际制冷量Q_o,s。将测得的运行参数代入公式2，可求得被测空调装置的理论制冷量Q_o,t，则被测空调装置的相对制冷量由下式可求：

由公式6可评判被测空调装置在任意测试工况下空调装置制冷能力是否正常。

本发明的有益效果：利用该测试方法制作的制冷量测试装置或仪器可以用于空调装置的制冷量测试，广泛用于交通车辆、船舶、民用与工业应用的空调装置制冷量测试。其技术较为先进，在国际领域中遥遥领先。

附图说明：

为了易于说明，本发明由下述的具体实施及附图作以详细描述。

图1为背景技术中现有领域中铁路客车空调装置试验工况图；

图2为背景技术中现有领域中基于空气焓差法的环路式空调装置性能测试装置图；

图3为本发明算法中空调装置系统的输入、输出图；

图4为本发明算法中蒸发器表面湿空气处理过程焓湿图；

图5为干、湿分区示意图。

具体实施方式：

本具体实施方式采用以下技术方案：本具体实施方式基于多维曲线拟合算法的空调装置制冷量测试方法依据的基本原理是：建立空调装置的基准样机物理模型，以此为研究对象，通过建立空调装置基准样机的制冷量与各个影响参数之间的多维曲线拟合数学模型，并采用包括名义工况在内的多组样本数据进行训练仿真，获得该空调装置在任意工况下的制冷量数学模型。通过该数学模型，可以获得任意工况下空调装置基准样机的制冷量，作为同类型被测空调装置制冷量的理论值。在相同测试工况下，利用该模型结果与实际测试数据的分析比较，获得被测空调装置实际制冷量。

参照以下实施方式来进行具体说明：空调装置制冷量的关联因素众多，除内部参数及其它次要因素外，主要包括蒸发器回风温度t_h、蒸发器回风相对湿度蒸发器通风量V_l、冷凝器进风温度t_c、冷凝器通风量V_c及大气压力P_a，对于指定的一台空调装置来说，只要确定上述六个参数，机组的能耗及制冷量也是唯一确定的；

将空调装置视为热力系统，整个过程的输入输出如图3所示，试验表明，蒸发器回风温度t_h、蒸发器回风湿度及冷凝器送风温度t_c对机组制冷量的影响最为明显，由于大气压力变化较小，其变化对空调装置制冷量的影响可以忽略；蒸发器、冷凝器通风量参数在测量精度要求的范围内可以作为常量处理，即在实际测试过程中认为空调装置冷凝器和蒸发器的送风量及大气压力基本维持不变。因此，制冷量(或功率)在测量过程中的变量主要是：蒸发器回风温度t_h、蒸发器回风湿度及冷凝器进风温度t_c；

上述变量中的蒸发器回风相对湿度主要通过蒸发器凝结水影响翅片管蒸发器的管外壁对流换热系数，从而影响蒸发器的传热量，根据传热学理论和湿空气热质交换理论，当蒸发器回风相对湿度达到一定值蒸发器回风湿度的影响表现在两个区域，干工况区和湿工况区，在这两个区内，蒸发器回风湿度的影响相对不显著，蒸发器表面温度t_e,b对应的饱和水蒸气分压与蒸发器回风温度t_h下的相对湿度处理过程焓湿图如图4所示，对相对湿度进行分段研究，达到在相对湿度各分段中空调装置制冷量的影响因素降维的目的，将蒸发器回风相对湿度按照分区后，即可将模型降为不同回风相对湿度区段的三维模型，利用Matlab软件，在各相对湿度分区中，对制冷量与蒸发器回风温度t_h、冷凝进风温度t_c的三维曲线进行拟合，以某种待测试的KLD40型单元式空调机组为例，给出具体拟合结果:

临界相对湿度的数学模型的建立及干湿分区的判别：

在蒸发器表面温度10～11℃、大气压力101325Pa时，借助Matlab软件分析拟合，得到蒸发器回风温度t_h与相应干、湿分区的临界相对湿度数学模型(公式1)，模型的适用范围：10℃≤t_h≤50℃

式中：

p₁＝6.2029e-009,p₂＝-1.5991e-006,p₃＝0.00018189

p₄＝-0.012068,p₅＝0.51418,p₆＝-14,p₇＝199.02

在实际工况下，当测得的相对湿度大于时，认定该工况点位于湿工况区，当测得的相对湿度小于时，认定该工况点位于干工况区，分区示意图如图5。

该机组在大气压力B(101325Pa)、送风量6000m³/h、冷凝风量11800m³/h的条件下。结合设计分析和实验研究成果，在蒸发器表面温度10～11℃时，得到制冷量与蒸发器回风温度、冷凝进风温度的关系如下：

Q_o＝p₀₀+p₁₀t_c+p₀₁t_h+p₂₀t_c²+p₁₁t_ct_h+p₀₂t_h²+p₃₀t_c³+p₂₁t_c²t_h+p₁₂t_ct_h²+p₀₃t_h³(公式2)

湿工况区:

p₀₀＝-26.9,p₁₀＝9.023,p₀₁＝-6.165,p₂₀＝-0.2055,p₁₁＝-0.1919

p₀₂＝0.4461,p₃₀＝0.001953,p₂₁＝0.0003906,p₁₂＝0.002734

p₀₃＝-0.007161

干工况区:

p₀₀＝-16.91,p₁₀＝6.02,p₀₁＝-4.223,p₂₀＝-0.1047，p₁₁＝-0.2256，

p₀₂＝0.4188，p₃₀＝0.001302，p₂₁＝-0.0007813，p₁₂＝0.004687，

p₀₃＝-0.007813

通过对基准样机以模型中的三个输入参数作变工况性能数据训练，不断修正各参数的权值与阈值，可获得被测空调机组的精度更高的制冷量数学模型，在此基础上，可进一步得到各种工况下的制冷量、耗电量的测量输出；例如，对该型号空调机组测量了一组数据：相对湿度蒸发器回风温度t_h＝27℃、冷凝进风温度t_c＝27℃，此时，可得到该机组理论制冷量Q_o,s＝43.76kW。

KLD40型单元式空调机组在大气压力B(101325Pa)、送风量6000m³/h、冷凝风量11800m³/h的条件下,名义工况(蒸发器回风温度29℃，湿球温度23℃，相对湿度60％，冷凝进风温度35℃)下制冷量为40.7kW，用本模型计算得到蒸发器回风温度29℃，相对湿度60％，冷凝进风温度35℃时的理论制冷量为41.21kW，相对误差为1.3％。所以，该空调机组制冷量模型准确度高。

测试工况下(任意运行工况)被测空调装置的实际制冷量Q_o,s：

Q_o,s＝V_l·ρ(h_h,j-h_h,c)>

测试工况下(任意运行工况)被测空调装置的理论制冷量Q_o,t:

则两式相除，得：

即

实际测试工况下被测空调装置的理论制冷量Q_o,t由上述制冷量数学模型(见公式2)可求出，又由公式4知：

$h_{h,j}^0-h_{h,c}^0=\frac{Q_{o,t}}{V_l·\rho }$

上式右侧均为已知参数，则理论焓降可以求得，此外，公式5右侧实际焓降h_h,j-h_h,c可在测试过程中计算得到，因此，利用公式5可以直接求得被测空调装置的实际制冷量。

在实际测量过程中测得被测空调装置的运行参数，通过上述测量空调装置制冷量的新方法，获得实际测试工况下空调装置的实际制冷量Q_o,s。将测得的运行参数代入公式2，可求得被测空调装置的理论制冷量Q_o,t，则被测空调装置的相对制冷量由下式可求：

由公式6可评判被测空调装置在任意测试工况下空调装置制冷能力是否正常。

以上显示和描述了本发明的基本原理和主要特征和本发明的优点。本行业的技术人员应该了解，本发明不受上述实施例的限制，上述实施例和说明书中描述的只是说明本发明的原理，在不脱离本发明精神和范围的前提下，本发明还会有各种变化和改进，这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。