CN103868624A - Room air conditioner refrigerating capacity collecting error correcting method - Google Patents

Abstract

The invention discloses a room air conditioner refrigerating capacity collecting error correcting method capable of being used for correcting refrigerating capacity test results of room air conditioners at different refrigerating capacity grades under the refrigerating working condition. According to the method, 14 groups of room air conditioner refrigerating capacity data under the test working condition are collected, and each constant coefficient value in a refrigerating capacity regression equation is confirmed through a Bessel's formula, so a refrigerating capacity deviation correcting equation is determined; the refrigerating capacity of a model machine at any wet and dry bulb temperature in the national standard is collected, and a refrigerating capacity deviation correcting value is determined through the refrigerating capacity deviation correcting equation, so a regression correcting value from the actual tested working condition to the nominal working condition is obtained. The invention provides the refrigerating capacity collecting error correcting method applied to room air conditioner performance evaluation, the method has the advantages that the evaluation accuracy can be improved, the room air conditioner performance can be really reflected, and the predicted error division caused by artificial control on energy efficiency levels of the room air conditioners can be avoided; meanwhile, the operation parameter setting of variable-frequency air conditioner equipment can be guided, and the regulation on the room air state can be more precise.

Description

A kind of Cooling Capacity For The Room Air Conditioner Acquisition Error modification method

Technical field

The present invention relates to a kind of refrigerating capacity Acquisition Error modification method, is a kind of method that is applied to the correction of Cooling Capacity For The Room Air Conditioner Acquisition Error specifically.

Background technology

Room air conditioner Energy Efficiency Standard is the Energy Efficiency Standard of the consumer worked out the earliest of China, and front and back have experienced three revisions, and the related specifications of the U.S., European Union and Japan has been followed in the formulation of Energy Efficiency Standard and revision.The revision of 2000 and twice Energy Efficiency Standard in 2004 is mainly reflected in the rise of efficiency limit value, and existing room air conditioner Energy Efficiency Standard, its efficiency grade is changed to three grades of efficiencies of existing GB GB12021.3-2010 by the Pyatyi efficiency of former GB GB12021.3-2004: existing Energy Efficiency Standard is directly deleted three, four and Pyatyi efficiency grade in former Energy Efficiency Standard, and former Energy Efficiency Standard I and II efficiency grade is redefined to two, three grades of efficiency grades into existing Energy Efficiency Standard.Waiting under the condition of step-length based on this, using the existing Energy Efficiency Standard secondary efficiency first order of one-level as act.std of boosting.

For the correct efficiency grade of judging room air conditioner, in " the thermodynamics sophistication of room air conditioner is analyzed " literary composition, to introduce thermodynamics sophistication room air conditioner is carried out to efficiency evaluation, this is a kind of room air conditioner energy efficiency analysis method for air meriting attention.Meanwhile, room air conditioner Energy Efficiency Standard upgrades revision, corresponding, and its performance standard also needs revision, and matches, and on to greatest extent, avoids because measuring former thereby causing the mistake of room air conditioner efficiency grade to divide.

Specified refrigerating capacity and Energy Efficiency Ratio need reflect the performance of room air conditioner under nominal condition point, and this is also the foundation of room air conditioner efficiency ranking.When Cooling Capacity For The Room Air Conditioner and Energy Efficiency Ratio are surveyed, the impact of tested person condition and human factor, the actual condition point of performance of room air conditioners test tends to be offset in tolerance.Correlation standard the franchise of wet and dry bulb temperature be respectively ± 1 ℃ and ± 0.5 ℃.This test result that often causes refrigerating capacity and Energy Efficiency Ratio is not at nominal condition point, but under a certain operating point, records in wet and dry bulb tolerance.Based on this test result, room air conditioner is carried out to efficiency ranking, may cause the mistake of room air conditioner efficiency grade to be divided.On the one hand, under nominal condition, can reach the room air conditioner of (not reaching) three grades of efficiency limit values, may test operating mode and be positioned at wet and dry bulb temperature tolerance lower limit (upper limit) due to reality, can not reach (reaching) three grades of efficiency limit values and be divided into by mistake, and three grades of threshold values that efficiency limit value is the room air conditioner market access.On the other hand, along with instrument and meter precision improves constantly, the actual test of manual control operating point is positioned at wet and dry bulb temperature high tolerance or lower limit, also may cause room air conditioner that predictable wrong division can occur.

Summary of the invention

Technical matters: for the problems referred to above, the invention provides one and be applied to room air conditioner, can obtain wet and dry bulb temperature in national standard tolerance arbitrarily under operating point refrigerating capacity to the recurrence modified value of nominal condition, improve the Cooling Capacity For The Room Air Conditioner Acquisition Error modification method of test and appraisal accuracy.

Technical scheme: Cooling Capacity For The Room Air Conditioner Acquisition Error modification method of the present invention, comprises the following steps:

1) be captured in the Cooling Capacity For The Room Air Conditioner data under operating condition of test, described operating condition of test is that dry-bulb temperature is 35 ℃ outside holding chamber, wet-bulb temperature is 24 ℃, inside difference pulpit, wet and dry bulb temperature is in 3 large classes, under 14 kinds of operating modes, carry out the measurement of Cooling Capacity For The Room Air Conditioner, the deviation that records refrigerating capacity under refrigerating capacity and nominal condition is Δ W _i(i=1,2,, 14), i represents different test operating modes:

Wherein, first kind operating mode is that indoor wet-bulb temperature is constant, changes indoor dry-bulb temperature; Operating mode 1～5 indoor wet-bulb temperature is 19 ℃, corresponding 26 ℃, 26.5 ℃, 27 ℃, 27.5 ℃, 28 ℃ respectively of indoor dry-bulb temperatures;

Equations of The Second Kind operating mode is that indoor dry-bulb temperature is constant, changes indoor wet-bulb temperature; Operating mode 6～10 indoor dry-bulb temperatures are 27 ℃, corresponding 18.5 ℃, 18.8 ℃, 19 ℃, 19.3 ℃, 19.5 ℃ respectively of indoor wet-bulb temperature;

The 3rd class is indoor wet and dry bulb temperature limit couple variations, and operating mode 11,12 indoor dry-bulb temperatures are 26 ℃, corresponding 18.5 ℃, 19.5 ℃ respectively of indoor wet-bulb temperature; Operating mode 13,14 indoor dry-bulb temperatures are 28 ℃, corresponding 18.5 ℃, 19.5 ℃ respectively of indoor wet-bulb temperature;

2) by the following regression equation of each floor data difference substitution gathering in described step 1), obtain the overdetermined equation group of 14 hexa-atomic linear function compositions, then take Bessel Formula as criterion solves described system of equations, obtain every constant coefficient W of regression equation ₀, a, b, c, d, e value:

$\mathrm{\ΔW}(t_1,t_2)=W_0+{\mathrm{at}}_1+{\mathrm{bt}}_2+{\mathrm{ct}}_1^2+{\mathrm{dt}}_2^2+{\mathrm{et}}_1{t}_2,$

Wherein, t ₁represent indoor dry-bulb temperature and the nominal condition point temperature difference of indoor dry-bulb temperature down under tested operating mode, t ₂represent indoor wet-bulb temperature and the nominal condition point temperature difference of indoor wet-bulb temperature down under tested operating mode;

W ₀represent that in regression equation, constant term, a represent temperature difference item t of indoor dry-bulb temperature ₁coefficient, b represent temperature difference item t of indoor wet-bulb temperature ₂coefficient, c be indoor dry-bulb temperature secondary temperature difference item t ² ₁coefficient, d represent indoor wet-bulb temperature secondary temperature difference item t ² ₂coefficient, e represent indoor wet and dry bulb temperature coupling temperature difference item t ₁t ₂coefficient;

3) in wet and dry bulb temperature tolerance, the refrigerating capacity drift correction amount under any practical running operating point, can calculate according to following formula:

ΔW(t ₁,t ₂)=W ₀+ζΛQt+(Qt) ^TΛQt；

Wherein, t is temperature difference vector, t=(t ₁, t ₂) ^t, T represents vectorial transposition;

Λ is weight matrix, Λ=diag (θ ₁, θ ₂), diag represents diagonal matrix, θ ₁represent the weighing factor of the indoor dry bulb temperature difference to refrigerating capacity deviation, θ ₂represent the weighing factor of indoor wet-bulb depression to refrigerating capacity deviation;

Q is temperature transition matrix, $Q=\left(\begin{array}{cc}{\mathrm{\α}}_{11}& {\mathrm{\α}}_{12}\\ {\mathrm{\α}}_{21}& {\mathrm{\α}}_{22}\end{array}\right),$ α ₁₁represent the influence coefficient of the indoor dry bulb temperature difference to indoor dry-bulb temperature tolerance under test operating mode, α ₂₁represent the influence coefficient of the indoor dry bulb temperature difference to indoor wet bulb bulb temperature tolerance under test operating mode, α ₁₂represent the influence coefficient of indoor wet-bulb depression to indoor dry-bulb temperature tolerance under test operating mode, α ₂₂represent the influence coefficient of indoor wet-bulb depression to indoor wet-bulb temperature tolerance under test operating mode;

ζ be temperature difference item with respect to the high order temperature difference item relative effect coefficient vector to refrigerating capacity, ζ=(ζ ₁, ζ ₂); The temperature term that ζ 1 represents the indoor dry bulb temperature difference is the relative effect coefficient to refrigerating capacity deviation with respect to high order temperature term, ζ ₂a temperature term of expression indoor wet-bulb depression is the relative effect coefficient to refrigerating capacity deviation with respect to high order temperature term;

4) the refrigerating capacity drift correction value Δ W (t obtaining according to described step 3) ₁, t ₂) the refrigerating capacity Acquisition Error of tested room air conditioner is revised:

W _t=W _m-ΔW(t ₁,t ₂)；

Wherein, W _mfor the refrigerating capacity measured value under any operating mode of tested room air conditioner in wet and dry bulb temperature tolerance; Δ W (t ₁, t ₂) be the drift correction value under this operating mode calculating in described step 3), W _tfor refrigerating capacity under this operating mode returns modified value.

The present invention revises the refrigerating capacity measured value under any operating mode in wet and dry bulb tolerance, obtain refrigerating capacity under this operating mode recurrence modified value to nominal condition, can be used for the actual operation parameters adjustment of convertible frequency air-conditioner: in the time that refrigerating capacity recurrence modified value is greater than refrigerating capacity measured value, suitably turn down compressor operating frequency; In the time that refrigerating capacity recurrence modified value is less than refrigerating capacity measured value, suitably heighten compressor operating frequency; In the time that refrigerating capacity recurrence modified value equals refrigerating capacity measured value, do not need to adjust.

In the inventive method, described step 2) in Bessel Formula be:

ε=k σ, wherein ε represents to monitor residual error, and k represents to monitor coefficient, and σ represents the standard deviation of refrigerating capacity actual measurement deviation, $\mathrm{\σ}=\sqrt{\frac{\underset{i=1}{\overset{14}{\mathrm{\Σ}}}{[{\mathrm{\ΔW}}_i-\frac{1}{14}\underset{i=1}{\overset{14}{\mathrm{\Σ}}}{\mathrm{\ΔW}}_i]}^2}{13}},(i=\mathrm{1,2},...,14);$

Bessel Formula decision condition is: monitoring coefficient k value 0.65,2,3 successively, in the time of k=0.65, has at least 50% refrigerating capacity actual measurement deviation to meet relational expression ε _i≤ ε; In the time of k=2, have at least 95% refrigerating capacity actual measurement deviation to meet relational expression ε _i≤ ε; In the time of k=3, have at least 100% refrigerating capacity actual measurement deviation to meet relational expression ε _i≤ ε;

ε _irepresent the measurement error of refrigerating capacity under different operating points and the residual error of equation deviation,

ε _i=| Δ W _i-Δ W _i(t ₁, t ₂); Δ W _irepresent the difference of refrigerating capacity measured value and nominal value; Δ W _i(t ₁, t ₂) represent the difference of refrigerating capacity functional value and nominal value; I represents different test operating modes, i=1, and 2 ..., 14.

In the inventive method, in described step 3), ${\mathrm{\θ}}_1=\frac{{\mathrm{\λ}}_1}{\left|{\mathrm{\λ}}_1\right|+\left|{\mathrm{\λ}}_2\right|},{\mathrm{\θ}}_2=\frac{{\mathrm{\λ}}_2}{\left|{\mathrm{\λ}}_1\right|+\left|{\mathrm{\λ}}_2\right|};$

${\mathrm{\α}}_{11}=e\sqrt{\left|{\mathrm{\λ}}_1\right|+\left|{\mathrm{\λ}}_2\right|}/\left|\right|p_1\left|\right|,{\mathrm{\α}}_{12}=(c-d-\sqrt{{(c-d)}^2+e^2})\sqrt{\left|{\mathrm{\λ}}_1\right|+\left|{\mathrm{\λ}}_2\right|}/\left|\right|p_1\left|\right|;$

${\mathrm{\α}}_{21}=e\sqrt{\left|{\mathrm{\λ}}_1\right|+\left|{\mathrm{\λ}}_2\right|}/\left|\right|p_2\left|\right|,{\mathrm{\α}}_{22}=(c-d-\sqrt{{(c-d)}^2+e^2})\sqrt{\left|{\mathrm{\λ}}_1\right|+\left|{\mathrm{\λ}}_2\right|}/\left|\right|p_2\left|\right|;$

${\mathrm{\ζ}}_1=\frac{({\mathrm{a\α}}_{22}-{\mathrm{b\α}}_{12})}{{\mathrm{\θ}}_1({\mathrm{\α}}_{11}{\mathrm{\α}}_{22}-{\mathrm{\α}}_{12}{\mathrm{\α}}_{21})},{\mathrm{\ζ}}_2=\frac{({\mathrm{b\α}}_{11}-{\mathrm{a\α}}_{21})}{{\mathrm{\θ}}_2({\mathrm{\α}}_{11}{\mathrm{\α}}_{22}-{\mathrm{\α}}_{12}{\mathrm{\α}}_{21})};$

Wherein, λ ₁, λ ₂for the eigenwert of secondary temperature difference item matrix of coefficients B, p ₁, p ₂for the proper vector of secondary temperature difference item matrix of coefficients B, $B=\left(\begin{array}{cc}c& e/2\\ e/2& d\end{array}\right),{\mathrm{\&lambda;}}_1=\frac{c+d-\sqrt{{(c-d)}^2+e^2}}{2},{\mathrm{\&lambda;}}_2=\frac{c+d-\sqrt{{(c-d)}^2+e^2}}{2}$ $\left|{\mathrm{\&lambda;}}_1\right|<\left|{\mathrm{\&lambda;}}_2\right|;\left|\right|p_1\left|\right|={(e^2+{[c-d-\sqrt{{(c-d)}^2+e^2}]}^2)}^{\frac{1}{2}},\left|\right|p_2\left|\right|={(e^2+{[c-d+\sqrt{{(c-d)}^2+e^2}]}^2)}^{\frac{1}{2}}).$

Beneficial effect: the present invention compared with prior art, has the following advantages:

1. the refrigerating capacity measured value in the limited discrete wet and dry bulb temperature tolerance of the inventive method application, describes the continuous variation characteristic of refrigerating capacity about wet and dry bulb temperature.

2. the standard deviation of surveying deviation function deviation take refrigerating capacity is as object, take the numerical relation of the residual sum standard deviation of test value and equation value as target, each test operating point is carried out to dynamic monitoring, with specific features and the accuracy thereof of determining that refrigerating capacity deviation changes with wet-bulb depression continuity.

The inventive method can obtain wet and dry bulb temperature in national standard tolerance arbitrarily under operating point refrigerating capacity to the recurrence modified value of nominal condition, thereby instruct the setting of the concrete operational factors of air-conditioning equipment such as convertible frequency air-conditioner, equipment operation is more accurate;

The inventive method can obtain wet and dry bulb temperature in national standard tolerance arbitrarily under operating point refrigerating capacity to the modified value of nominal condition, improve test and appraisal accuracy, reflect more really performance of room air conditioners, and reference be provided to the revision of national standard;

5. refrigerating capacity drift correction formula is once determined, the Cooling Capacity For The Room Air Conditioner measured value of follow-up identical refrigerating capacity grade is directly determined by this correction formula, avoids repeatedly the huge consumption of the human and material resources that revision test brings;

6. along with instrument and meter precision improves constantly, avoid the actual test of manual control operating point to be positioned at wet and dry bulb temperature high tolerance or lower limit, thereby cause room air conditioner efficiency grade that predictable wrong division occurs.

Accompanying drawing explanation

Fig. 1 is the process step figure of the inventive method.

Embodiment

Below in conjunction with drawings and Examples, the present invention will be further described:

Cooling Capacity For The Room Air Conditioner Acquisition Error modification method of the present invention can the refrigeration capacity test result under cooling condition be revised the room air conditioner of different refrigerating capacity grades.Be modified to embodiment with the refrigerating capacity measured result under three class room air conditioners (refrigerating capacity is between between 7.1kW and 14.0kW) cooling condition below and describe, concrete steps are as follows:

Step 1) is captured in the Cooling Capacity For The Room Air Conditioner data under operating condition of test, described operating condition of test is that dry-bulb temperature is 35 ℃ outside holding chamber, wet-bulb temperature is 24 ℃, inside difference pulpit, wet and dry bulb temperature is in 3 large classes, under 14 kinds of operating modes, carry out the measurement of Cooling Capacity For The Room Air Conditioner, the deviation that records refrigerating capacity under refrigerating capacity and nominal condition is Δ W ₁=-0.0185kW, Δ W ₂=-0.0058kW, Δ W ₃=0kW, Δ W ₄=0.0049kW, Δ W ₅=0.0075kW, Δ W ₆=-0.2157kW, Δ W ₇=-0.0843kW, Δ W ₈=0kW, Δ W ₉=0.1133kW, Δ W ₁₀=0.1548kW, Δ W ₁₁=-0.22kW, Δ W ₁₂=0.1392kW, Δ W ₁₃=-0.2098kW, Δ W ₁₄=0.1951kW:

Wherein, first kind operating mode is that indoor wet-bulb temperature is constant, changes indoor dry-bulb temperature; Operating mode 1～5 indoor wet-bulb temperature is 19 ℃, corresponding 26 ℃, 26.5 ℃, 27 ℃, 27.5 ℃, 28 ℃ respectively of indoor dry-bulb temperatures;

Equations of The Second Kind operating mode is that indoor dry-bulb temperature is constant, changes indoor wet-bulb temperature; Operating mode 6～10 indoor dry-bulb temperatures are 27 ℃, corresponding 18.5 ℃, 18.8 ℃, 19 ℃, 19.3 ℃, 19.5 ℃ respectively of indoor wet-bulb temperature;

The 3rd class is indoor wet and dry bulb temperature limit couple variations, and operating mode 11,12 indoor dry-bulb temperatures are 26 ℃, corresponding 18.5 ℃, 19.5 ℃ respectively of indoor wet-bulb temperature; Operating mode 13,14 indoor dry-bulb temperatures are 28 ℃, corresponding 18.5 ℃, 19.5 ℃ respectively of indoor wet-bulb temperature;

Step 2) the each floor data gathering in described step 1) is distinguished to the following regression equation of substitution:

$\mathrm{\&Delta;W}(t_1,t_2)=W_0+{\mathrm{at}}_1+{\mathrm{bt}}_2+{\mathrm{ct}}_1^2+{\mathrm{dt}}_2^2+{\mathrm{et}}_1{t}_2,$

Wherein, t ₁represent indoor dry-bulb temperature and the nominal condition point temperature difference of indoor dry-bulb temperature down under tested operating mode, t ₂represent indoor wet-bulb temperature and the nominal condition point temperature difference of indoor wet-bulb temperature down under tested operating mode;

W ₀represent that in regression equation, constant term, a represent temperature difference item t of indoor dry-bulb temperature ₁coefficient, b represent temperature difference item t of indoor wet-bulb temperature ₂coefficient, c be indoor dry-bulb temperature secondary temperature difference item t ₁ ²coefficient, d represent indoor wet-bulb temperature secondary temperature difference item

coefficient, e represent indoor wet and dry bulb temperature coupling temperature difference item t ₁t ₂coefficient;

Under the 3rd, 6 operating modes, t ₁, t ₂,

,

, t ₁t ₂, Δ W _i(t ₁, t ₂) value be 0, therefore obtain the overdetermined equation group of 12 hexa-atomic linear functions compositions:

W ₀-a+c=-0.0185；

W ₀-0.5a+0.25c=-0.0058；

W ₀+0.5a+0.25c=0.0049；

W ₀+a+c=0.0075；

W ₀-0.5b+0.25d=-0.2157；

W ₀-0.2b+0.04=-0.0843；

W ₀+0.3b+0.09d=0.1133；

W ₀+0.5b+0.25d=0.1548；

W ₀-a-0.5b+c+0.25d+0.5e=-0.22；

W ₀-a+0.5b+c+0.25d-0.5e=0.1392；

W ₀+a-0.5b+c+0.25d-0.5e=-0.2098；

W ₀+a+0.5b+c+0.25d+0.5e=0.1951

Solve above-mentioned system of equations, obtain every constant coefficient W of regression equation ₀, a, b, c, d, e value be respectively 0.0019,0.19,0.37,0.0027,0.011,0.0098;

Meanwhile, adopt the convergence of Bessel Formula checking refrigerating capacity regression equation:

The residual error of test value and equation value is ε _i=| Δ W _i-Δ W _i(t ₁, t ₂) |, Δ W _irepresent the difference of refrigerating capacity measured value and nominal value, Δ W _i(t ₁, t ₂) represent the difference of refrigerating capacity functional value and nominal value; Trying to achieve according to above-mentioned formula the residual error that each operating point is corresponding is respectively: ε ₁=0.01582, ε ₂=0.0025, ε ₃=0.00681, ε ₄=0.00669, ε ₅=0.01079, ε ₆=0.01058, ε ₇=0.00075, ε ₈=0.01078, ε ₉=0.0106, ε ₁₀=0.00709, ε ₁₁=0.00146, ε ₁₂=0.00511, ε ₁₃=0.00863, ε ₁₄=0.01483

σ represents the standard deviation of refrigerating capacity actual measurement deviation,

$\mathrm{\&sigma;}=\sqrt{\frac{\underset{i=1}{\overset{14}{\mathrm{\&Sigma;}}}{[{\mathrm{\&Delta;W}}_i-\frac{1}{14}\underset{i=1}{\overset{14}{\mathrm{\&Sigma;}}}{\mathrm{\&Delta;W}}_i]}^2}{13}}=0.012,(i=\mathrm{1,2},...,14);$

Bessel Formula is ε=k σ, respectively when monitoring coefficient k value 0.65,2,3 successively, in the time of k=0.65, has 57% refrigerating capacity actual measurement deviation to meet relational expression ε _i≤ ε, (i represents different test operating modes, i=1, and 2 ..., 14); In the time of k=2, there is 100% refrigerating capacity actual measurement deviation to meet relational expression ε _i≤ ε; In the time of k=3, there is 100% refrigerating capacity actual measurement deviation to meet relational expression ε _i≤ ε.Therefore the convergence of refrigerating capacity regression equation, W ₀, a, b, c, d, e value be respectively 0.0019,0.19,0.37,0.0027,0.011,0.0098;

$\begin{array}{c}\mathrm{\&Delta;W}(t_1,t_2)=W_0+{\mathrm{at}}_1+{\mathrm{bt}}_2+{\mathrm{ct}}_1^2+{\mathrm{dt}}_2^2+{\mathrm{et}}_1{t}_2\\ =0.0019+{0.19t}_1+{0.37t}_2+0.0027{t_1}^2+0.011{t_2}^2+0.0098t_1{t}_2\end{array}$

Step 3) is in wet and dry bulb temperature tolerance, and (specified refrigerating capacity and Energy Efficiency Ratio are respectively 12kW, 3.32kWkW to a room air conditioner model machine of the present invention ^-1) actual measurement refrigerating capacity W _m=12.2783kW (operating condition of test point: indoor dry-bulb temperature is that 26.7 ℃, wet-bulb temperature are 18.9 ℃), the refrigerating capacity drift correction amount under this operating mode, can calculate according to following formula:

ΔW(t ₁,t ₂)=W ₀+ζΛQt+(Qt) ^TΛQt；

Wherein, t is temperature difference vector, t=(t ₁, t ₂) ^t, T represents vectorial transposition;

Λ is weight matrix, Λ=diag (θ ₁, θ ₂), diag represents diagonal matrix, θ ₁represent the weighing factor of the indoor dry bulb temperature difference to refrigerating capacity deviation, θ ₂represent the weighing factor of indoor wet-bulb depression to refrigerating capacity deviation;

Q is temperature transition matrix, $Q=\left(\begin{array}{cc}{\mathrm{\&alpha;}}_{11}& {\mathrm{\&alpha;}}_{12}\\ {\mathrm{\&alpha;}}_{21}& {\mathrm{\&alpha;}}_{22}\end{array}\right),$ α ₁₁represent the influence coefficient of the indoor dry bulb temperature difference to indoor dry-bulb temperature tolerance under test operating mode, α ₂₁represent the influence coefficient of the indoor dry bulb temperature difference to indoor wet bulb bulb temperature tolerance under test operating mode, α ₁₂represent the influence coefficient of indoor wet-bulb depression to indoor dry-bulb temperature tolerance under test operating mode, α ₂₂represent the influence coefficient of indoor wet-bulb depression to indoor wet-bulb temperature tolerance under test operating mode;

ζ be temperature difference item with respect to the high order temperature difference item relative effect coefficient vector to refrigerating capacity, ζ=(ζ ₁, ζ ₂); ζ ₁a temperature term of the expression indoor dry bulb temperature difference is the relative effect coefficient to refrigerating capacity deviation with respect to high order temperature term, ζ ₂a temperature term of expression indoor wet-bulb depression is the relative effect coefficient to refrigerating capacity deviation with respect to high order temperature term;

θ ₁, θ ₂, α ₁₁, α ₁₂, α ₂₁, α ₂₂, ζ ₁, ζ ₂determine by the following method:

Coefficient c, the d of secondary temperature term in refrigerating capacity regression equation, e are carried out to quadratic form and change to obtain Second Order with Constant Coefficients matrix $B=\left(\begin{array}{cc}c& e/2\\ e/2& d\end{array}\right),$ Determine its eigenvalue λ ₁, λ ₂and characteristic of correspondence vector p ₁, p ₂:

${\mathrm{\&lambda;}}_1=\frac{c+d+\sqrt{{(c-d)}^2+e^2}}{2},{\mathrm{\&lambda;}}_2=\frac{c+d+\sqrt{{(c-d)}^2+e^2}}{2},$

$p_1={(e,c-d-\sqrt{{(c-d)}^2+e^2})}^T,p_2={(e,c-d+\sqrt{{(c-d)}^2+e^2})}^T;$

, ${\mathrm{\&theta;}}_1=\frac{{\mathrm{\&lambda;}}_1}{\left|{\mathrm{\&lambda;}}_1\right|+\left|{\mathrm{\&lambda;}}_2\right|},{\mathrm{\&theta;}}_2=\frac{{\mathrm{\&lambda;}}_2}{\left|{\mathrm{\&lambda;}}_1\right|+\left|{\mathrm{\&lambda;}}_2\right|},$ Weight matrix Λ=diag (θ ₁, θ ₂);

${\mathrm{\&alpha;}}_{11}=e\sqrt{\left|{\mathrm{\&lambda;}}_1\right|+\left|{\mathrm{\&lambda;}}_2\right|}/\left|\right|p_1\left|\right|,{\mathrm{\&alpha;}}_{12}=(c-d-\sqrt{{(c-d)}^2+e^2})\sqrt{\left|{\mathrm{\&lambda;}}_1\right|+\left|{\mathrm{\&lambda;}}_2\right|}/\left|\right|p_1\left|\right|,$

${\mathrm{\&alpha;}}_{21}=e\sqrt{\left|{\mathrm{\&lambda;}}_1\right|+\left|{\mathrm{\&lambda;}}_2\right|}/\left|\right|p_2\left|\right|,{\mathrm{\&alpha;}}_{22}=(c-d-\sqrt{{(c-d)}^2+e^2})\sqrt{\left|{\mathrm{\&lambda;}}_1\right|+\left|{\mathrm{\&lambda;}}_2\right|}/\left|\right|p_2\left|\right|,$ Temperature transition matrix $Q=\left(\begin{array}{cc}{\mathrm{\&alpha;}}_{11}& {\mathrm{\&alpha;}}_{12}\\ {\mathrm{\&alpha;}}_{21}& {\mathrm{\&alpha;}}_{22}\end{array}\right);$ ${\mathrm{\&zeta;}}_1=\frac{({\mathrm{a\&alpha;}}_{22}-{\mathrm{b\&alpha;}}_{12})}{{\mathrm{\&theta;}}_1({\mathrm{\&alpha;}}_{11}{\mathrm{\&alpha;}}_{22}-{\mathrm{\&alpha;}}_{12}{\mathrm{\&alpha;}}_{21})},{\mathrm{\&zeta;}}_2=\frac{({\mathrm{b\&alpha;}}_{11}-{\mathrm{a\&alpha;}}_{21})}{{\mathrm{\&theta;}}_2({\mathrm{\&alpha;}}_{11}{\mathrm{\&alpha;}}_{22}-{\mathrm{\&alpha;}}_{12}{\mathrm{\&alpha;}}_{21})},$ Temperature difference item is the relative effect coefficient vector ζ=(ζ to refrigerating capacity with respect to high order temperature difference item ₁, ζ ₂);

Be that refrigerating capacity drift correction formula is:

$\begin{array}{c}\mathrm{\&Delta;W}(t_1,t_2)=0.0019+\left(\mathrm{5.57,3.54}\right)\left(\begin{array}{cc}0.03& \\ & 0.97\end{array}\right)\left(\begin{array}{cc}0.42& 0.91\\ 0.91& -0.42\end{array}\right)\left(\begin{array}{c}{t}_1\\ t_2\end{array}\right)\\ +{\left(\left(\begin{array}{cc}0.42& 0.91\\ 0.91& -0.42\end{array}\right)\left(\begin{array}{c}{t}_1\\ t_2\end{array}\right)\right)}^T\left(\begin{array}{cc}0.03& \\ & 0.97\end{array}\right)\left(\begin{array}{cc}0.42& 0.91\\ 0.91& -0.42\end{array}\right)\left(\begin{array}{c}{t}_1\\ t_2\end{array}\right)\end{array}$

Calculate the refrigerating capacity drift correction amount under this operating mode according to above-mentioned solving result:

$\begin{array}{c}\mathrm{\&Delta;W}\left(\mathrm{26.7,18.9}\right)=0.0019+\left(\mathrm{5.57,3.54}\right)\left(\begin{array}{cc}0.03& \\ & 0.97\end{array}\right)\left(\begin{array}{cc}0.42& 0.91\\ 0.91& -0.42\end{array}\right)\left(\begin{array}{c}26.7\\ 18.9\end{array}\right)\\ +{\left(\left(\begin{array}{cc}0.42& 0.91\\ 0.91& -0.42\end{array}\right)\left(\begin{array}{c}26.7\\ 18.9\end{array}\right)\right)}^T\left(\begin{array}{cc}0.03& \\ & 0.97\end{array}\right)\left(\begin{array}{cc}0.42& 0.91\\ 0.91& -0.42\end{array}\right)\left(\begin{array}{c}26.7\\ 18.9\end{array}\right)=-0.09146\mathrm{kW};\end{array}$

The refrigerating capacity drift correction value Δ W (t that step 4) obtains according to described step 3) ₁, t ₂) the refrigerating capacity Acquisition Error of tested room air conditioner is revised:

W _t=W _m-ΔW(t ₁,t ₂)=12.2090kW+0.09146kW=12.3005kW

Be dry-bulb temperature be 26.7 ℃, wet-bulb temperature be at 18.9 ℃ refrigerating capacity to return modified value be 12.3005kW.

Wherein, W _mfor the refrigerating capacity measured value under any operating mode of tested room air conditioner, be 12.2090kW; Δ W (t ₁, t ₂) be the drift correction value under this operating mode calculating in described step 3), for-0.09146kW; W _tfor refrigerating capacity under this operating mode returns modified value.

Claims (3)

1. a Cooling Capacity For The Room Air Conditioner Acquisition Error modification method, is characterized in that, the method comprises the following steps:

1) be captured in the Cooling Capacity For The Room Air Conditioner data under operating condition of test, described operating condition of test is that dry-bulb temperature is 35 ℃ outside holding chamber, wet-bulb temperature is 24 ℃, inside difference pulpit, wet and dry bulb temperature is in 3 large classes, under 14 kinds of operating modes, carry out the measurement of Cooling Capacity For The Room Air Conditioner, the deviation that records refrigerating capacity under refrigerating capacity and nominal condition point is Δ W _i(i=1,2,, 14), i represents different test operating modes:

Wherein, first kind operating mode is that indoor wet-bulb temperature is constant, changes indoor dry-bulb temperature; Operating mode 1～5 indoor wet-bulb temperature is 19 ℃, corresponding 26 ℃, 26.5 ℃, 27 ℃, 27.5 ℃, 28 ℃ respectively of indoor dry-bulb temperatures;

Equations of The Second Kind operating mode is that indoor dry-bulb temperature is constant, changes indoor wet-bulb temperature; Operating mode 6～10 indoor dry-bulb temperatures are 27 ℃, corresponding 18.5 ℃, 18.8 ℃, 19 ℃, 19.3 ℃, 19.5 ℃ respectively of indoor wet-bulb temperature;

The 3rd class is indoor wet and dry bulb temperature limit couple variations, and operating mode 11,12 indoor dry-bulb temperatures are 26 ℃, corresponding 18.5 ℃, 19.5 ℃ respectively of indoor wet-bulb temperature; Operating mode 13,14 indoor dry-bulb temperatures are 28 ℃, corresponding 18.5 ℃, 19.5 ℃ respectively of indoor wet-bulb temperature;

2) by the following regression equation of each floor data difference substitution gathering in described step 1), obtain the overdetermined equation group of 14 hexa-atomic linear function compositions, then take Bessel Formula as criterion solves described system of equations, obtain every constant coefficient W of regression equation ₀, a, b, c, d, e value:

$\mathrm{\&Delta;W}(t_1,t_2)=W_0+{\mathrm{at}}_1+{\mathrm{bt}}_2+{\mathrm{ct}}_1^2+{\mathrm{dt}}_2^2+{\mathrm{et}}_1{t}_2,$

Wherein, t ₁represent indoor dry-bulb temperature and the nominal condition point temperature difference of indoor wet and dry bulb temperature down under tested operating mode, t ₂represent indoor wet-bulb temperature and the nominal condition point temperature difference of indoor wet-bulb temperature down under tested operating mode;

W ₀represent that in regression equation, constant term, a represent temperature difference item t of indoor dry-bulb temperature ₁coefficient, b represent temperature difference item t of indoor wet-bulb temperature ₂coefficient, c be indoor dry-bulb temperature secondary temperature difference item t ² ₁coefficient, d represent indoor wet-bulb temperature secondary temperature difference item t ² ₂coefficient, e represent indoor wet and dry bulb temperature coupling temperature difference item t ₁t ₂coefficient;

3) in wet and dry bulb temperature tolerance, the refrigerating capacity drift correction amount under any practical running operating point, can calculate according to following formula:

ΔW(t ₁,t ₂)=W ₀+ζΛQt+(Qt) ^TΛQt；

Wherein, t is temperature difference vector, t=(t ₁, t ₂) ^t, T represents vectorial transposition;

Λ is weight matrix, Λ=diag (θ ₁, θ ₂), diag represents diagonal matrix, θ ₁represent the weighing factor of the indoor dry bulb temperature difference to refrigerating capacity deviation, θ ₂represent the weighing factor of indoor wet-bulb depression to refrigerating capacity deviation;

Q is temperature transition matrix, $Q=\left(\begin{array}{cc}{\mathrm{\&alpha;}}_{11}& {\mathrm{\&alpha;}}_{12}\\ {\mathrm{\&alpha;}}_{21}& {\mathrm{\&alpha;}}_{22}\end{array}\right),$ α ₁₁represent the influence coefficient of the indoor dry bulb temperature difference to indoor dry-bulb temperature tolerance under test operating mode, α ₂₁represent the influence coefficient of the indoor dry bulb temperature difference to indoor wet bulb bulb temperature tolerance under test operating mode, α ₁₂represent the influence coefficient of indoor wet-bulb depression to indoor dry-bulb temperature tolerance under test operating mode, α ₂₂represent the influence coefficient of indoor wet-bulb depression to indoor wet-bulb temperature tolerance under test operating mode;

ζ be temperature difference item with respect to the high order temperature difference item relative effect coefficient vector to refrigerating capacity, ζ=(ζ ₁, ζ ₂); ζ ₁a temperature term of the expression indoor dry bulb temperature difference is the relative effect coefficient to refrigerating capacity deviation with respect to high order temperature term, ζ ₂a temperature term of expression indoor wet-bulb depression is the relative effect coefficient to refrigerating capacity deviation with respect to high order temperature term;

4) the refrigerating capacity drift correction value Δ W (t obtaining according to described step 3) ₁, t ₂) the refrigerating capacity Acquisition Error of tested room air conditioner is revised:

W _t=W _m-ΔW(t ₁,t ₂)；

Wherein, W _mfor the refrigerating capacity measured value under any operating mode of tested room air conditioner in wet and dry bulb temperature tolerance; Δ W (t ₁, t ₂) be the drift correction value under this operating mode calculating in described step 3), W _tfor refrigerating capacity under this operating mode returns modified value.

2. Cooling Capacity For The Room Air Conditioner Acquisition Error modification method according to claim 1, is characterized in that described step 2) in Bessel Formula be: ε=k σ, wherein ε represents to monitor residual error, k represents to monitor coefficient, and σ represents the standard deviation of refrigerating capacity actual measurement deviation $\mathrm{\&sigma;}=\sqrt{\frac{\underset{i=1}{\overset{14}{\mathrm{\&Sigma;}}}{[{\mathrm{\&Delta;W}}_i-\frac{1}{14}\underset{i=1}{\overset{14}{\mathrm{\&Sigma;}}}{\mathrm{\&Delta;W}}_i]}^2}{13}},(i=\mathrm{1,2},...,14);$

Described Bessel Formula decision condition is: monitoring coefficient k value 0.65,2,3 successively, in the time of k=0.65, has at least 50% refrigerating capacity actual measurement deviation to meet relational expression ε _i≤ ε; In the time of k=2, have at least 95% refrigerating capacity actual measurement deviation to meet relational expression ε _i≤ ε; In the time of k=3, have at least 100% refrigerating capacity actual measurement deviation to meet relational expression ε _i≤ ε;

ε _irepresent the measurement error of refrigerating capacity under different operating modes and the residual error of equation deviation, ε _i=| Δ W _i-Δ W _i(t ₁, t ₂); Δ W _irepresent the difference of refrigerating capacity measured value and nominal value; Δ W _i(t ₁, t ₂) represent the difference of refrigerating capacity functional value and nominal value; I represents different test operating modes, i=1, and 2 ... 14.

3. Cooling Capacity For The Room Air Conditioner Acquisition Error modification method according to claim 1, is characterized in that, in described step 3), ${\mathrm{\&theta;}}_1=\frac{{\mathrm{\&lambda;}}_1}{\left|{\mathrm{\&lambda;}}_1\right|+\left|{\mathrm{\&lambda;}}_2\right|},{\mathrm{\&theta;}}_2=\frac{{\mathrm{\&lambda;}}_2}{\left|{\mathrm{\&lambda;}}_1\right|+\left|{\mathrm{\&lambda;}}_2\right|};$

${\mathrm{\&alpha;}}_{11}=e\sqrt{\left|{\mathrm{\&lambda;}}_1\right|+\left|{\mathrm{\&lambda;}}_2\right|}/\left|\right|p_1\left|\right|,{\mathrm{\&alpha;}}_{12}=(c-d-\sqrt{{(c-d)}^2+e^2})\sqrt{\left|{\mathrm{\&lambda;}}_1\right|+\left|{\mathrm{\&lambda;}}_2\right|}/\left|\right|p_1\left|\right|;$

${\mathrm{\&alpha;}}_{21}=e\sqrt{\left|{\mathrm{\&lambda;}}_1\right|+\left|{\mathrm{\&lambda;}}_2\right|}/\left|\right|p_2\left|\right|,{\mathrm{\&alpha;}}_{22}=(c-d-\sqrt{{(c-d)}^2+e^2})\sqrt{\left|{\mathrm{\&lambda;}}_1\right|+\left|{\mathrm{\&lambda;}}_2\right|}/\left|\right|p_2\left|\right|;$

${\mathrm{\&zeta;}}_1=\frac{({\mathrm{a\&alpha;}}_{22}-{\mathrm{b\&alpha;}}_{12})}{{\mathrm{\&theta;}}_1({\mathrm{\&alpha;}}_{11}{\mathrm{\&alpha;}}_{22}-{\mathrm{\&alpha;}}_{12}{\mathrm{\&alpha;}}_{21})},{\mathrm{\&zeta;}}_2=\frac{({\mathrm{b\&alpha;}}_{11}-{\mathrm{a\&alpha;}}_{21})}{{\mathrm{\&theta;}}_2({\mathrm{\&alpha;}}_{11}{\mathrm{\&alpha;}}_{22}-{\mathrm{\&alpha;}}_{12}{\mathrm{\&alpha;}}_{21})};$

Wherein, λ ₁, λ ₂for the eigenwert of secondary temperature difference item matrix of coefficients B, p ₁, p ₂for the proper vector of secondary temperature difference item matrix of coefficients B, $B=\left(\begin{array}{cc}c& e/2\\ e/2& d\end{array}\right),{\mathrm{\&lambda;}}_1=\frac{c+d-\sqrt{{(c-d)}^2+e^2}}{2},{\mathrm{\&lambda;}}_2=\frac{c+d-\sqrt{{(c-d)}^2+e^2}}{2}$ $\left|{\mathrm{\&lambda;}}_1\right|<\left|{\mathrm{\&lambda;}}_2\right|;\left|\right|p_1\left|\right|={(e^2+{[c-d-\sqrt{{(c-d)}^2+e^2}]}^2)}^{\frac{1}{2}},\left|\right|p_2\left|\right|={(e^2+{[c-d+\sqrt{{(c-d)}^2+e^2}]}^2)}^{\frac{1}{2}}).$

Images