武連發(fā)，羅亞軍，熊建國(guó)，余　凱

( 珠海格力電器股份有限公司，珠海 519070 )

?

關(guān)于高溫多聯(lián)機(jī)用制冷劑“氣體的比容、粘度”的擬合計(jì)算

武連發(fā)，羅亞軍，熊建國(guó)，余凱

( 珠海格力電器股份有限公司，珠海 519070 )

在開(kāi)發(fā)R410a高溫多聯(lián)機(jī)過(guò)程中，前期會(huì)進(jìn)行連接管(過(guò)熱氣體管、過(guò)冷液體管)的壓降仿真計(jì)算，其中經(jīng)常會(huì)涉及到查找制冷劑物性參數(shù)，或需要在計(jì)算中直接根據(jù)已有的幾個(gè)參數(shù)來(lái)快速計(jì)算制冷劑的其他物性參數(shù)。其中使用比較多但較難查到的是“過(guò)熱氣體的比容、運(yùn)動(dòng)粘度，過(guò)冷液體的比容、運(yùn)動(dòng)粘度”四個(gè)參數(shù)。經(jīng)過(guò)研究，在已知飽和壓力、過(guò)熱度、過(guò)冷度等參數(shù)下，可以擬合比容、運(yùn)動(dòng)粘度。在實(shí)際空調(diào)運(yùn)行范圍內(nèi)，其誤差達(dá)到1%以下。采用VB編程，編制了快速計(jì)算軟件《R410A氣體比容粘度計(jì)算V1.0》，方便其他工程人員使用。[關(guān)鍵詞]高溫多聯(lián)機(jī)；仿真；R410A；比容；粘度；氣體

1　引言

在開(kāi)發(fā)R410a高溫多聯(lián)機(jī)過(guò)程中，前期會(huì)進(jìn)行連接管(過(guò)熱氣體管、過(guò)冷液體管)的壓降理論計(jì)算，其中經(jīng)常會(huì)涉及到查找制冷劑物性參數(shù)，或需要在計(jì)算中直接根據(jù)已有的幾個(gè)參數(shù)來(lái)快速計(jì)算制冷劑的其他物性參數(shù)。另外，制冷劑熱力性質(zhì)的計(jì)算是制冷空調(diào)裝置仿真的重要基礎(chǔ)[1]。這就需要我們能找到各個(gè)參數(shù)的計(jì)算公式。制冷劑熱力性質(zhì)的計(jì)算是相關(guān)領(lǐng)域計(jì)算機(jī)輔助設(shè)計(jì)的主要基礎(chǔ)之一，其計(jì)算精度、速度和穩(wěn)定性將直接影響、甚至決定整個(gè)計(jì)算機(jī)輔助設(shè)計(jì)軟件的計(jì)算精度、速度和穩(wěn)定性[2]。

其中在理論計(jì)算中，使用比較多的是過(guò)熱氣體的比容、運(yùn)動(dòng)粘度。這兩個(gè)參數(shù)“過(guò)熱氣體的比容、運(yùn)動(dòng)粘度”，目前可提供的擬合公式較少，且較難找到可查找這兩個(gè)參數(shù)的軟件。

本文就R410A制冷劑“氣體的比容、粘度”的擬合計(jì)算，進(jìn)行相關(guān)探討。所采用的擬合軟件為Origin軟件，Origin軟件功能強(qiáng)大，是科學(xué)和工程研究人員必備軟件之一[3]。

2　關(guān)于氣體比容的擬合過(guò)程

采用的方法是，先確定此點(diǎn)的“壓力、過(guò)熱度”，然后，先擬合同壓力“TE”下飽和氣體的比容“v0”；之后再根據(jù)“不同的壓力TE、不同的過(guò)熱度△TE”，擬合不同的比容增量“△v”。那么氣體的比容“v1”=“v0+△v”。

2.1同壓力“TE”下飽和氣體的比容“v0”擬合所用數(shù)據(jù)如表1。

表1　同壓力下飽和氣體比容擬合所用數(shù)據(jù)

采用Origin軟件，進(jìn)行擬合得到如下公式：

v0= 32.7345+B1*TE^1+B2*TE^2+B3*TE^3

其中：B1=-1.06623、B2=0.01992、B3= -2.16162E-4。擬合誤差在0.5%以下；擬合曲線如圖1。

圖1　同壓力下飽和氣體的比容擬合曲線

2.2“不同的壓力TE、不同的過(guò)熱度△TE”，擬合不同的比容增量“△v”，所用數(shù)據(jù)如表2。

表2　不同壓力不同過(guò)熱度擬合不同比容增量所用數(shù)據(jù)

采用Origin軟件，進(jìn)行擬合得到如下公式：每隔5℃過(guò)熱度的比容增量平均值=(“△v”/△TE)*5=(1.06146+B4*TE^1+B5*TE^2)

其中：B4=-0.02599、B5=3.82381E-4。擬合誤差在1.1%以下；擬合曲線如圖2。

圖2　不同壓力、不同過(guò)熱度擬合不同比容增量

2.3氣體的比容擬合公式為

“v1”=“v0+△v”=32.7345+B1*TE^1+B2*TE^2+B3*TE^3+(1.06146+B4*TE^1+B5*TE^2)/5*△TE

其中：B1=-1.06623、B2=0.01992、B3=-2.16162 E-4、B4=-0.02599、B5=3.82381E-4。擬合誤差在1%以下；

擬合曲線如圖3。

圖3　氣體的比容擬合曲線

3　關(guān)于氣體運(yùn)動(dòng)粘度的擬合過(guò)程

方法同上述的氣體比容的擬合方法，先確定此點(diǎn)的“壓力、過(guò)熱度”，然后，先擬合同壓力“TE”下飽和氣體的運(yùn)動(dòng)粘度“ν0”；之后再根據(jù)“不同的壓力TE、不同的過(guò)熱度△TE”，擬合不同的粘度增量“△ν”。那么氣體的運(yùn)動(dòng)粘度“ν1”=“ν0+△ν”。

3.1同壓力“TE”下飽和氣體的運(yùn)動(dòng)粘度“ν0”擬合結(jié)果如下：

ν0=0.38483+B6*TE^1+B7*TE^2+B8*TE^3

其中：B6=-0.01044、B7=1.77662E-4、B8=-1.73737E-6。擬合誤差在0.5%以下；擬合曲線如圖4。

圖4　同壓力下飽和氣體的運(yùn)動(dòng)粘度擬合曲線

3.2“不同的壓力TE、不同的過(guò)熱度△TE”，擬合不同的粘度增量“△ν”，結(jié)果如下：

每隔5℃過(guò)熱度的運(yùn)動(dòng)粘度增量平均值=(“△ν”/△TE)*5=(0.02288+B9*TE^1+B10*TE^2)

其中：B9=-5.74048E-4、B10=7.7619E-6。擬合誤差在1.5%以下；擬合曲線如圖5。

圖5　不同壓力、不同過(guò)熱度擬合不同粘度增量

3.3氣體的運(yùn)動(dòng)粘度擬合公式為：

“ν1”=“ν0+△ν”=0.38483+B6*TE^1+B7*TE^2+B8*TE^3+(0.02288+B9*TE^1+B10*TE^2)/5*△TE

其中：B6=-0.01044、B7=1.77662E-4、B8=-1.73737E-6、B9=-5.74048E-4、B10=7.7619E-6。擬合誤差在1%以下；擬合曲線如圖6。

圖6　氣體的運(yùn)動(dòng)粘度擬合曲線

4　計(jì)算軟件的編制

為利于其他工程人員快速計(jì)算這4個(gè)參數(shù)，采用Visual Basic6.0編程。Visual Basic是可視化的“BASIC”語(yǔ)言工具[4]。筆者編制了快速計(jì)算軟件《R410A氣體比容粘度計(jì)算V1.0》，其他工程人員只需要簡(jiǎn)單輸入幾個(gè)參數(shù)，軟件就會(huì)自動(dòng)計(jì)算所需的比容、粘度等參數(shù)。

軟件分為三個(gè)部分，輸入?yún)?shù)、輸出參數(shù)、說(shuō)明。輸入?yún)?shù)部分，輸入氣體對(duì)應(yīng)的飽和壓力(單位℃)、氣體過(guò)熱度(單位℃)；輸入?yún)?shù)部分，在按中間“計(jì)算”按鈕后，可根據(jù)已輸入的參數(shù)自動(dòng)計(jì)算氣體比容、氣體運(yùn)動(dòng)粘度；說(shuō)明部分，列舉了輸入?yún)?shù)的范圍和版本信息。

軟件界面如圖7。

圖7　軟件界面

5　結(jié)論

綜上分析，采用“氣體飽和壓力”加“氣體過(guò)熱度”，來(lái)擬合氣體比容、氣體運(yùn)動(dòng)粘度公式。在實(shí)際空調(diào)運(yùn)行范圍內(nèi)，其誤差達(dá)到1%以下。為利于其他工程人員快速計(jì)算這些參數(shù)，采用VB編程，編制了快速計(jì)算軟件《R410A氣體比容粘度計(jì)算V1.0》，方便其他工程人員查詢。

[1] 吳志剛，丁國(guó)良.制冷劑熱力性質(zhì)的快速計(jì)算Ⅰ.計(jì)算方法[J].上海交通大學(xué)學(xué)報(bào)，2006，40(2)：297-300

[2] 張春路，丁國(guó)良，李灝.制冷劑過(guò)熱氣體熱力性質(zhì)的隱式擬合方法[J].工程熱物理學(xué)報(bào)，2000，21(5)：533-536

[3] 肖信.Origin 8.0實(shí)用教程：科技作圖與數(shù)據(jù)分析[M].北京：中國(guó)電力出版社，2009

[4] 楊培添.Visual Basic程序設(shè)計(jì)教程[M].北京：人民郵電出版社，2006

Coupling Calculation about Gaseous Specific Volume and Viscidity of Refrigerant for High Sensitive Heat VRF

WU Lianfa，LUO Yajun，XIONG Jianguo，YU Kai

(Gree Electric Appliances,Inc of Zhuhai，Zhuhai519070，China))

When highsensitive heat VRFwith refrigerant of R410A was exploited,the simulae calculation about pressure loss in gas pipe or liquid pipe was done in the prophase.In this process,lots of characteristic parameters about refrigerant were often demanded,at the same time,other characteristic parameters about refrigerant were calculated frequently according some existing parameters.Four parameters of those that specific volume or viscosity of super heating refrigerant and super cooling refrigerant were used regularly,yet found difficulty.By researching,specific volume and viscosity could be found out on condition that saturated temperature and super-cooling temperature or super-heating temperature was known.In range of air conditioning running,the precision was low of 1%.The software of specific volume and viscosity about R410A with gas was worked by VB,which could be helpful to others.

Highsensitive heat VRF；Simulate；R410A；Specificvolume；Viscosity；Gas

2016-4-24

國(guó)家科技支撐計(jì)劃課題(2014BAJ02B01)

武連發(fā)(1981-)，男，工程師，研究方向：多聯(lián)機(jī)的系統(tǒng)節(jié)能與產(chǎn)品研發(fā)。Email：lianfawu@sina.com

ISSN1005-9180(2016)03-070-04

TB61+2文獻(xiàn)標(biāo)示碼：Adoi:10.3696/J.ISSN.1005-9180.2016.03.013