沈　希，熊　奇，顧江萍，金華強，徐　鳴

(浙江工業(yè)大學 機械工程學院，浙江 杭州 310014)

?

基于排氣壓力的壓縮機轉速測量方法研究

沈希，熊奇，顧江萍，金華強，徐鳴

(浙江工業(yè)大學 機械工程學院，浙江 杭州 310014)

摘要：壓縮機的轉速是測量壓縮機性能的一個重要指標之一，它影響著制冷系統(tǒng)的制冷量、效能比以及功耗、吸排氣壓力等，但是對于全封閉式的壓縮機傳統(tǒng)的測量方法(如編碼器)已不能適用.為了快速精確測量全封閉壓縮機轉速，提出了一種基于排氣壓力的周期性波動來測量壓縮機轉速的方法，通過LABVIEW軟件平臺對排氣壓力信號進行數(shù)據采集、信號分析、結果顯示.該方法不僅適用于小型制冷往復式壓縮機，也適用于周期性排氣壓力變化的壓縮機，具有一定的推廣意義.

關鍵詞：往復式壓縮機;轉速測量;LABVIEW; Chirp-Z

壓縮機作為制冷系統(tǒng)的核心部件，其技術性指標直接影響著整個制冷系統(tǒng)的性能，而壓縮機的運行本質上是電動機的轉動，也就是壓縮機轉速.轉速影響著制冷系統(tǒng)的制冷量、效能比以及功耗、吸排氣壓力等，因此壓縮機轉速是一項重要的技術指標.目前，小型制冷壓縮機都是全封閉型的，傳統(tǒng)的測量方法存在一定的缺陷.其中，直接法[1]采用在壓縮機內部裝傳感器來測轉速，但這無疑會對壓縮機造成一定的損傷，不符合測量的原則；另一種是間接法[2]，目前主要有測壓縮機的振動頻率和壓縮機電流.前者由于壓縮機運行的同時有其他接近壓縮機振動的頻率干擾,不易分辨;后者在測量裝置和精度上有一定的局限性[3].針對上述方法的這些局限，筆者根據壓縮機和制冷系統(tǒng)排氣壓力波動，提出了一種運用數(shù)字信號處理方法進行精確轉速測量的方法.

1壓縮機轉速測量原理

制冷壓縮機的工作過程是通過曲柄連桿機構和活塞在氣缸內實現(xiàn)，活塞的往復運動會使氣缸內的壓力呈現(xiàn)周期性變化，根據壓強-容積圖的關系和現(xiàn)有的科學實驗[4]可得在曲柄轉過的一個周期內，氣缸內的壓力有周期性變化.同理可得，在壓縮機排氣口會有一個微弱的周期性壓力信號[5]，根據這個壓力信號即可求得壓縮機的轉速.

通過采集系統(tǒng)得到這個微弱的周期性壓力信號后必須進行信號處理，常用的信號頻域處理方法有很多，比如快速傅里葉變換(FFT)、線性調頻Z變換(CZT).

1.1　快速傅里葉變換

快速傅里葉變換[6]是頻譜分析最常用的處理方法，但它存在一定的缺點[7].頻譜泄漏，由于對信號進行截斷，有限區(qū)間以外的值都為零，使得信號未能在頻譜中完全顯示；柵欄效應，進行快速傅里葉變換時是用一個個離散的點去逼近連續(xù)頻譜值，因此只能觀察到有限個頻譜值.若采樣頻率為10 kHz，采樣點數(shù)為10 000，那么頻率分辨率就只有1 Hz，也就是說轉速誤差會有60 r/min，這對轉速的測量誤差實在太大.若增大采樣點數(shù)，則會增加采樣時間，無法進行快速測量.

1.2　線性調頻Z變換

作為數(shù)字信號處理的一種算法，它是對頻譜的局部“細化”，能夠看到快速傅里葉變換看不到的頻率信號[8].CZT是沿Z平面上的一段螺線作等分角抽樣[9],z的取樣點為

zk=AW-k

(1)

式中：A=A0ejθ0;W=W0e-jφ0；k為0到M的實數(shù)；M為欲分析的復頻譜的點數(shù).

M不一定等于N，這其實是能夠細化頻譜的關鍵.zk各點在Z平面的值如圖1所示，它正是沿著這段螺線來獲得特定區(qū)間的頻譜.對zk進行z變換可得

(2)

(3)

與FFT相比，CZT的頻譜分辨率不受采樣頻率和采樣個數(shù)的影響，精度更高，但會增加處理器資源.筆者通過采集排氣壓力信號，運用數(shù)字信號處理方法快速得到壓縮機的精確轉速.

圖1　CZT的變換路徑Fig.1　Transformation path of CZT

2具體實施方案

2.1　硬件平臺的搭建

測試系統(tǒng)采用研華工控機，研華PCI-1715U高速數(shù)據采集卡，雁樺GY403型壓力傳感器及LABVIEW開發(fā)平臺，對壓縮機系統(tǒng)中的排氣口壓力進行采集.

由于排氣壓力信號的頻率在50 Hz左右，為了能夠采集到它的變化就必須是高頻響的快變式傳感器，因此選用了GY403型壓力傳感器，它的精度較高，頻率響應能達到10 K，同時傳感、變送一體化的封裝，省去了后續(xù)的電橋檢測及放大等信號調理環(huán)節(jié)，大大減少了電磁干擾的影響.

為了準確的采集到信號，數(shù)據采集卡也有一定的要求，要選用分辨率較高的數(shù)據采集卡，因此采用研華PCI-1715U高速數(shù)據采集卡.PCI-1715U高速數(shù)據采集卡提供了32個模擬量輸入通道，采樣頻率能達到500 kS/s、12位分辨率及2 500 VDC的直流隔離保護.根據奈奎斯特采樣定理，而且往復式壓縮機的轉速大概在2~3 kr/min，因此綜合考慮采樣頻率設置為10 k/s.測量裝置示意圖如圖2所示，壓縮機在運行的時候，壓力傳感器不斷地高速采集排氣壓力信號并通過高速采集卡送入計算機[10].

1—壓縮機；2—排氣壓力表；3—傳感器；4—冷凝器；5—膨脹閥；6—蒸發(fā)器；7—平衡閥；8—吸氣閥；9—吸氣壓力表圖2　硬件系統(tǒng)框圖Fig.2　Hardware system diagram

2.2　軟件平臺的實現(xiàn)

圖3是測試中軟件的主要流程，軟件平臺采用LABVIEW[11]，它是一種基于圖形化的編程語言和開發(fā)環(huán)境，與以往的那些基于文本的高級編程語言不同，采用一種全新的數(shù)據流圖編程方式，即用形象的圖標和連線來代替文本，讓人們從復雜的程序中解脫出來，研華公司提供了多種開發(fā)例程，方便用戶快速開發(fā)應用.

由于MATLAB具有強大的數(shù)學運算能力，因此本實驗是將LABVIEW中采集到的數(shù)據導入到MATLAB進行頻譜分析，這樣就大大減小處理時間[12].目前，實現(xiàn)在LABVIEW中調用MATLAB的方法主要有MATLAB Script、基于Active X技術、基于COM技術、動態(tài)鏈接庫技術(DLL).使用MATLAB Script可以很方便的實現(xiàn)多輸入多輸出的功能，但是對MATLAB控制能力較弱.MATLAB Script對輸入輸出數(shù)據有明確的要求，只支持Real，Complex，1-D Array of real，1-D Array of complex，2-D Array of real，2-D Array of complex，String和Path這8種數(shù)據類型，使用MATLAB Script不能脫離MATLAB單獨運行.基于COM技術和動態(tài)鏈接庫技術(DLL)實現(xiàn)的混合編程的方法都具有可以脫離MATLAB環(huán)境、代碼運算效率高、占用資源少等優(yōu)點，但它們都適合開發(fā)大型應用程序，在實際使用中開發(fā)難度大，需要較多的專業(yè)知識.由于MATLAB強大的計算能力，綜合考慮，筆者采用MATLAB Script方式，在MATLAB的環(huán)境下調試好腳本程序，然后將代碼復制到LABVIEW的MATLAB節(jié)點中去，并且使用MATLAB Script能夠直接在MATLAB節(jié)點中修改程序.

由于壓縮機啟動時速度很快，因此啟動轉速的測量本方法并不適用.當硬件設備運行穩(wěn)定后開始數(shù)據采集，但采集到的數(shù)據不能直接用來分析，首先需要預處理，尤其是消除趨勢項，這樣便于后期的信號變換.在測試系統(tǒng)中，采集到的壓力信號數(shù)據由于存在放大器隨溫度變化產生的零點漂移、傳感器頻率范圍外的低頻性能的不穩(wěn)定以及傳感器周圍的環(huán)境干擾往往會偏離基線,甚至偏離基線的大小還會隨時間變化.趨勢項的存在,會使時域中的相關分析和頻率中的功率譜分析產生大的誤差,甚至使低頻譜完全失去真實性，而這無疑會對轉速的測量產生較大的誤差，故采用最小二乘法來消除多項式趨勢.

在濾波方面，高頻信號對排氣壓力信號在頻譜分析中并不影響，使用低頻濾波即可.濾波是信號分析領域的重要組成部分，根據濾波器所處理信號不同分為模擬濾波器和數(shù)字濾波器.以往的濾波器大多采用模擬電路技術進行設計，但是模擬電路技術存在電壓漂移、溫度漂移和噪聲等很多難以解決的問題，而采用數(shù)字技術則避免很多類似的難題，并且數(shù)字濾波器在濾波效果方面也有很多突出的優(yōu)點，都是模擬技術所不能及的.巴特沃斯低通數(shù)字濾波器具有最大平坦幅頻響應的特性，且具有良好的線性相位特性等優(yōu)點，目前應用廣泛.巴特沃斯低通濾波器的階數(shù)越高，幅頻特性衰減的速度越快，就越接近理想幅頻特性，但實現(xiàn)起來電路越復雜.LABVIEW基于良好的圖形化編程方式，可以直接調用濾波模塊，對于該系統(tǒng)可采用3階巴特沃斯低通濾波器，設置相應的低截止頻率、采樣頻率和濾波階數(shù)，這樣就省去了搭建硬件的麻煩，大大降低了開發(fā)成本.之后就是信號變換，得到信號的頻率后乘以60，即為壓縮機的轉速，最后結果在LABVIEW平面上顯示.軟件流程如圖3所示.

圖3　軟件系統(tǒng)測量流程Fig.3　Software measurement process

3實驗結果分析

壓縮機運行時腔內的壓力是在不斷變化的，根據曲柄滑塊機構的原理可得在排氣口處會有一個減弱的壓力波動信號，圖4(a)為排氣口處采集到的壓力波動信號，可知排氣口存在微弱的周期性波動，振幅在0.008 V左右.

圖4　時域及頻域結果Fig.4　The time domain and frequency domain result

圖4(b，c)是壓縮機排氣口壓力信號經過FFT和CZT后的結果圖，實驗在標準工況下進行，排氣壓力為7 bar，電源頻率50 Hz.在FFT頻譜圖中頻率落在49 Hz上，可得轉速為2 940 r/min，它的分辨率為1 Hz，誤差達到60 r/min.在線性調頻Z變換后的頻譜圖中，采樣時間為1 s，分辨率為0.01，誤差為0.6 r/min，在1 s內相比于FFT精確度大大提高，得出的頻率值為48.87 Hz，即轉速為2 932.2 r/min，這與廠家提供的壓縮機額度轉速十分接近.

為了進一步驗證測量到的壓縮機轉速的精度,通過使用變頻壓縮機在不同的轉速情況下進行測量，系統(tǒng)采集到的數(shù)據會經過預處理和平均值處理，因此可以保證系統(tǒng)的魯棒性.從表1中分析采集1 s后的結果可知：由于頻率分辨率的限制，F(xiàn)FT不能準確的反應信號的精確頻率，雖然加大快速傅里葉變換的個數(shù)能提高分辨率，但是會增加測量時間(不補零).在相同的時間內，線性調頻Z變換得到了更加精確的轉速.為了進一步驗證實驗結果，進行了2 s的實驗，由表中的數(shù)據可以看出，隨著時間的增大，兩種變換的精確率都有所提高，F(xiàn)FT的頻率分辨率為0.5 Hz，因此頻率只可能落在0.5 Hz的整數(shù)倍上，相比與1 Hz有了一定的提高.對于CZT隨著時間的增加，轉速精度提高，但沒有FFT明顯.在工程應用中，由于工況的復雜性，比如存在一些噪聲，變頻干擾，要想獲得更加精確的頻率可以適當?shù)脑黾訙y量時間.

表1　轉速比較

4結論

通過對排氣壓力信號的時域波形分析及變頻驅動壓縮機實驗，驗證了此方法測量轉速的可行性.對CZT算法和FFT算法比較，CZT算法能夠快速得到精確轉速.由于實際工況的復雜性(如受到干擾等影響)，要想獲得更加準確的轉速，可以適當增加采樣時間以確保測量結果的準確性.筆者介紹的轉速測量方法不僅適用于小型往復式制冷壓縮機，也適用于排氣口有周期性壓力波動的壓縮機，因此具有廣泛的適用性，有一定的工程應用價值.

參考文獻：

[1]DEMAY M B, FLESCH C A, ROSA A P. Indirect measurement of hermetic compressor speed through externally-measurable quantities[J]. International journal of refrigeration,2011,34(5):1268-1275.

[2]席靜珠,趙國軍,計時鳴,等.全封閉式壓縮機轉速測量方法的研究[J].強度與環(huán)境,1993,15(3):38-42.

[3]沈希,王烈鑫,席靜珠,等.全封閉壓縮機轉速的測量[J].機電工程,1995(4):16-18.

[4]陳楠.活塞式制冷壓縮機PV圖測試系統(tǒng)研究[D].杭州：浙江工業(yè)大學,2012.

[5]黨錫淇,夏永源.往復式壓縮機管路內壓力脈動的測量[J].壓縮機技術,1979(4):28-33.

[6]童莉葛,李紅英,王立,等.頻譜分析法在測量封閉式制冷壓縮機轉速中的應用[J].測試技術學報,2007(3):202-206.

[7]楊麗娟,張白樺,葉旭楨.快速傅里葉變換FFT及其應用[J].光電工程,2004,31(1):1-3.

[8]徐建軍,胡光東,李錦明.線性調頻Z變換在信號頻譜分析中的應用[J].應用基礎與工程科學學報,2009,17(6):966-972.

[9]單傳賓.基于Chirp-Z變換算法在雷達物位計上的應用與研究[D].天津：天津理工大學,2014.

[10]徐丹陽,童建平,隋成華,等.光柵光譜儀中的高速數(shù)據采集系統(tǒng)[J].浙江工業(yè)大學學報,2014,42(6):690-693.

[11]肖瑾.基于LABVIEW的封閉式制冷壓縮機轉速測量[J].企業(yè)導報,2010,10(6):294-296.

[12]徐稼軒,邱陽,朱因遠,等.壓縮機管道固有頻率的敲擊響應FFT分析法[J].壓縮機技術,1983,36(3):31-35.

(責任編輯：劉巖)

Research of compressor speed measurement method based on discharge pressure

SHEN Xi, XIONG Qi, GU Jiangping, JIN Huaqiang, XU Ming

(College of Mechanical Engineering, Zhejiang University of Technology, Hangzhou 310014, China)

Abstract:The speed of the compressor is one of the most important index for measuring the compressor performance. It effects the refrigerating capacity of the refrigeration system, efficiency, power consumption, discharge pressure and so on. But for hermetic compressor, traditional measurement methods (such as encoder) can not be applied. In order to measure hermetic compressor speed quickly and accurately, proposed a method based on discharge pressure that had cyclical fluctuations to measure the rotational speed of compressor based on the LABVIEW software platform of discharge pressure signal data acquisition, signal analysis, display results. This method is not only suitable for small reciprocating refrigeration compressor, also applies to compressor that has periodic discharge pressure change, which has a certain popularization significance.

Keywords:reciprocating compressor; speed measurement; LABVIEW; Chirp-Z

中圖分類號：TP39

文獻標志碼：A

文章編號：1006-4303(2016)01-0011-04

作者簡介：沈希(1966—),男,浙江杭州人,教授,主要從事制冷系統(tǒng)檢測和控制的研究,E-mail:sx@zjut.edu.cn.

基金項目：國家自然科學基金資助項目(51076143)

收稿日期：2015-09-16