宋明毅,吳偉烽,李直

(西安交通大學能源與動力工程學院, 710049, 西安)

?

汽車空調(diào)壓縮機氣閥運動規(guī)律模擬

宋明毅,吳偉烽,李直

(西安交通大學能源與動力工程學院, 710049, 西安)

為了解決斜盤式汽車空調(diào)壓縮機因進、排氣閥片斷裂導致失效的問題,建立了汽車空調(diào)壓縮機工作過程的三維流固耦合模型,分析了壓縮機內(nèi)部工作過程和氣閥的運動規(guī)律,并采用汽車空調(diào)壓縮機性能實驗進行了驗證。研究結果表明:在進氣壓力為0.18 MPa、排氣壓力為1.5 MPa、轉(zhuǎn)速為2 000 r/min時,壓縮機的平均進氣壓力損失為0.046 MPa、平均排氣壓力損失為0.11 MPa、壓縮機的容積效率為64.98%;以進氣閥片是否延遲為選擇原則,閥片厚度存在一個最佳值0.281 mm,若大于該厚度,則進氣閥會關閉延遲,該結果與增加舌簧閥的閥片厚度可以縮短延遲關閉的結果相反。該結果對小型制冷壓縮機舌簧閥的設計與優(yōu)化具有指導意義。

汽車空調(diào)壓縮機;三維流固耦合;氣閥運動

汽車空調(diào)壓縮機是一種變轉(zhuǎn)速的壓縮機,其主軸與發(fā)動機之間一般通過皮帶輪連接,因此其轉(zhuǎn)速變化與發(fā)動機轉(zhuǎn)速變化一致。由于汽車空調(diào)壓縮機具有往復式壓縮機的特點,因此變轉(zhuǎn)速對壓縮機影響最大的是氣閥的可靠性和經(jīng)濟性。對氣閥的研究與試驗表明,設計良好的氣閥,其流動阻力損失約為壓縮機的軸功率的3%～7%,而設計較差的,其流動阻力損失約為壓縮機的軸功率的15%～20%[1]。按照單機功率為2 kW計算,目前中國的汽車保有量約有2 260萬輛,若每天有1/3的汽車開啟空調(diào)一個小時,其能量損失達226萬 kW。更為重要的是,實際工程中常出現(xiàn)壓縮機氣閥的斷裂現(xiàn)象,這與變轉(zhuǎn)速的工況是極為相關的。因此,設計經(jīng)濟性好、可靠性高的氣閥對汽車空調(diào)壓縮機至關重要。

已有很多學者對氣閥及其運動規(guī)律進行了研究,其中,國外學者對氣閥的研究相對較早,文獻[2]中指出,Costagliola基于質(zhì)量彈簧系統(tǒng)建立了閥片的準靜態(tài)方程,在此基礎之上,Costagliola、Linke等人通過理論和實驗研究,得到了經(jīng)典氣閥理論。Touber基于經(jīng)典氣閥理論建立了壓縮過程的數(shù)學模型,對系統(tǒng)的參數(shù)可以進行快速的估算,Glen利用動態(tài)應變片建立了閥片運動的測量系統(tǒng)[3];Buliga利用激光儀建立了閥片運動的測量系統(tǒng)[4];Ludu利用電子內(nèi)窺儀直接測量了閥片的運動規(guī)律[5];2006年Habing考慮流體滯后和閥片速度的影響,以消除理論預測與實驗的偏差[6]。在國內(nèi)學者中,林梅等人建立的平板系統(tǒng)的振幅最接近實際測量值[7]。目前的研究中主要基于多維瞬態(tài)控制體方程,忽略了氣閥與被壓縮氣體的相互作用[8]。李蓉通過流固耦合技術與實驗對比,得到小型制冷壓縮機更準確的氣閥運動規(guī)律[9]。2006年Kim等人通過實驗驗證ADINA流固耦合技術的可靠性[10],2010年Choi等人也通過模擬與實驗的對比驗證ADINA流固耦合技術可以很好的應用到壓縮機氣閥運動規(guī)律的模擬中[11]。綜上所述,考慮氣閥與被壓縮氣體之間的相互作用是必要的,但同時由于汽車空調(diào)壓縮機體積小、工況復雜,要通過實驗直接測量其p-V圖及其氣閥的運動規(guī)律極其困難。

另外一方面,目前沒有標準和成熟的理論可以指導實際工程中如何選取最優(yōu)的氣閥厚度。

鑒于流固耦合技術已經(jīng)成熟地應用于其他往復式壓縮機如冰箱壓縮機的研究中,本文擬利用ADINA軟件建立三維流固耦合模型對閥片厚度不同的汽車空調(diào)壓縮機氣閥運動規(guī)律進行了模擬。由于缺少汽車空調(diào)壓縮機氣閥運動規(guī)律的直接測試結果,本文采用壓縮機性能測試的結果對計算模型進行了對比,總體證明了計算模型的準確性。

1 計算模型與方法

1.1 計算工況

本文計算的模型為一臺變排量壓縮機,7缸單作用,缸徑31.61 mm,計算與實驗工況下的行程為28.9 mm,排量為163 mL/r。

汽車空調(diào)壓縮機的閥片為舌簧閥。在計算模型中,閥片的變形即閥片的運動規(guī)律是流固耦合計算的關鍵內(nèi)容。實際機型采用的閥片材料為SANDVIK特種鋼帶,其特性如表1所示。

表1 計算模型中的閥片材料特性[12]

制冷工質(zhì)為R134a,實驗工況下,壓縮機轉(zhuǎn)速為2 000 r/min,蒸發(fā)溫度為8.72 ℃,排氣溫度為81.96 ℃,過熱度為9.94 ℃,過冷度為5 ℃。對應工況下,斜盤傾角角度全開壓縮機行程為最大行程,進氣壓力為0.18 MPa,排氣壓力為1.5 MPa。

圖1為活塞位于上止點的幾何模型,主要有進氣腔、排氣腔和活塞端面,其中活塞端面位于上止點。

圖1 活塞位于上止點的幾何模型

圖2為活塞位于下止點的幾何模型,其中活塞端面位于下止點。

圖2 活塞位于下止點的幾何模型

1.2 計算方法

壓縮機氣閥運動會影響氣缸內(nèi)流體的運動,對流場的影響不能忽略,流體的運動反過來又會影響閥片的運動,這屬于雙向流固耦合問題。ADINA軟件具備強大的流固耦合分析功能,因此本文基于ADINA建立汽車空調(diào)壓縮機工作過程的三維流固耦合模型,對其工作過程進行理論分析,對閥片運動規(guī)律、特性進行深入研究,活塞端面設置為移動壁面Leader-Follower,并涉及參數(shù)化動網(wǎng)格(ALE技術)的使用[13]。

流場網(wǎng)格數(shù)為15 991,結構場網(wǎng)格數(shù)為8 721,全部為規(guī)則的六面體網(wǎng)格,時間步長為0.01 ms。

1.3 迭代耦合雙向計算方式

迭代耦合雙向計算的流體方程及結構方程分別為

(1)

(2)

2 計算結果與分析

單個時間步長內(nèi)最大的迭代次數(shù)為500,如果收斂就不繼續(xù)迭代,進而計算下一步,當?shù)?00次后還是不能收斂,計算終止,并給出不收斂的可能原因。采用FCB-I單元計算,可大大縮短時間,精度不會受到影響[14]。計算時間為8 h,模型運用標準k-ε方程求解。

2.1 結果驗證

本文模型的計算結果與實驗結果比較如表2所示。

表2 特征值的計算結果與實驗結果

由表2可知,容積效率相對誤差為3.55%,制冷系數(shù)(COP)的相對誤差為3.52%,都在可以接受的范圍內(nèi),說明計算結果是準確的。

2.2 壓縮機的內(nèi)部工作過程

壓縮機分為4個工作過程,吸氣、壓縮、排氣、膨脹,每個工作過程都對應于壓縮機的工作過程p-V曲線圖。圖3為數(shù)值模擬壓縮機工作過程p-V曲線圖。

圖3 壓縮機工作過程P-V圖

壓縮機的具體工作過程對應圖3如下:da為吸氣過程;ab為壓縮過程;bc為排氣過程;cd為膨脹過程。壓縮機的各個系數(shù)如下:容積系數(shù)λV=88.26%;壓力系數(shù)λp=88.7%;溫度系數(shù)λt=83%;進氣系數(shù)λs=88.26%×88.7%×83%=64.98%,實驗值是62.75%。

瞬時進氣壓力損失和瞬時排氣壓力損失也是壓縮機的重要參數(shù),為得到這些參數(shù),現(xiàn)取壓縮機進氣結束以及排氣開始時的壓力分布圖作為參考如圖4所示。

(a)進氣結束壓力分布圖

(b)排氣開始壓力分布圖圖4 進排氣壓力分布圖

圖4中得出瞬時進氣壓力損失為0.069 MPa;瞬時排氣壓力損失為0.58 MPa,排氣壓力損失遠大于進氣壓力損失。這是因為排氣閥片的厚度,彈性力遠大于進氣閥片。

2.3 進氣閥的運動規(guī)律

進氣閥片的設計優(yōu)先考慮運動特性以及應力分布,圖5為閥片網(wǎng)格劃分圖。為獲得進氣閥片的相關數(shù)據(jù),在圖5中進氣閥片的邊緣處取一點A,然后進行升程,對速度、應力進行數(shù)據(jù)統(tǒng)計,仿真結果如圖6所示。

圖5 進氣閥片網(wǎng)格劃分圖中A點位置

(a)進氣閥片升程速度曲線圖

(b)進氣閥片升程應力曲線圖圖6 進氣閥模型的升程速度曲線及接觸應力圖

如圖6a所示,撞擊速度最大值約為5 m/s。較大的撞擊速度,將會帶來更大的噪聲。根據(jù)動量守恒方程,在閥片末端自由處的微元上進行計算可得到該處的接觸應力在10-1N數(shù)量級。從圖6b還可以看出,閥片末端在撞擊氣缸頂部的時刻,應力值幾乎為0,所以并不會反彈,從閥片運動的視頻中也可以證明這一點,閥片回落的過程相對緩慢。

延遲關閉將會降低壓縮機排氣量,進而降低壓縮機的效率,根據(jù)之前文獻中描述的可以增加閥片的彈性力來減小延遲關閉角,對此增加閥片厚度,并進行了模擬分析,并選取如圖7中閥片頭部中心處B點作為參考點,仿真結果如圖8所示。

圖7 進氣閥片網(wǎng)格劃分圖中B點位置

由圖8可知,閥片厚度h為0.281 mm時,進氣閥關閉角為227.48°;閥片厚度為0.305 mm時,關閉角為228.20°;閥片厚度為0.381 mm時,關閉角為231.12°;閥片厚度為0.481 mm時,關閉角為248.24°。當厚度增加到0.481 mm時,出現(xiàn)嚴重的顫振現(xiàn)象。通過以上數(shù)據(jù)可以發(fā)現(xiàn)增加彈性力后,反而延遲現(xiàn)象更加嚴重,這與通常的理論預測結果相反。分析原因如下:閥片厚度增加時,氣閥將要關閉時的通流面積相對減小,導致同一時刻氣缸內(nèi)的壓力減小,閥片上的作用力相應增大,導致其延時關閉。為了進一步揭示其原因,首先給出不同閥片厚度進氣閥即將關閉的過程中氣缸內(nèi)壓力的變化,如圖9所示。

圖8 不同厚度閥片的升程曲線

圖9 氣閥關閉時刻閥片厚度對氣缸內(nèi)壓力的影響

從以上數(shù)據(jù)可以看出,在同一轉(zhuǎn)角因不同閥片厚度而導致氣缸內(nèi)的壓力不同,厚度越大,氣缸內(nèi)的壓力越小。這是因為閥片越厚,其回落時對應的轉(zhuǎn)角越小,導致進入氣缸內(nèi)的氣體減少,從而使得最終將要關閉時氣缸內(nèi)壓力更小,入口與氣缸內(nèi)壓力差就會更大,從而導致其延遲關閉。閥片厚度為0.305 mm時所對應的壓力與厚度為0.381 mm時所對應的壓力差約為4 kPa,閥片厚度差為0.076 mm,而閥片的進氣側的壓力都為入口壓力0.18 MPa,根據(jù)上述情況,基于ADINA的結構模塊在以下3種工況下進行模擬分析。工況1:在厚度為0.305 mm的進氣閥片的進氣口端施加壓力1 kPa;工況2:在厚度為0.305 mm的進氣閥片的進氣口端施加壓力5 kPa;工況3:在厚度為0.381 mm的進氣閥片的進氣口端施加壓力1 kPa。并在閥片的頂端取一點C作為位移參考點,如圖10所示。

圖10 進氣閥片位移參考點

對參考點C在3種工況下的變形量進行統(tǒng)計,結果見表3。

表3 不同工況下的參考點C位移

注:Δx為與工況1中同一參考點位移的差值。

從表3中的數(shù)據(jù)可以明顯看出,因壓力增加引起的位移要大于因厚度減小而增加的位移。所以,厚度增加(即彈簧力增加)將導致進氣閥片延遲關閉現(xiàn)象更加嚴重。同時可以看出,閥片厚度增加后,時間截面變小,這將導致壓縮機的效率下降。

通過數(shù)值模擬技術,還可以得到閥片整體運動的彈性變化模式,在同一時刻對不同厚度的舌簧閥沿長度方向相同位置取點,統(tǒng)計閥片升程方向的位移,如圖11所示。

圖11 同時刻不同厚度閥片的長度方向位移

從圖11可以明顯看出,隨著厚度的增加,閥片的升程面積減小。這意味著進氣截面積減小,導致氣流絕對速度增加,流體壓降增加,功耗增大。

文獻[9]還指出:減小進氣閥片厚度對減小能耗具有積極的意義;減小厚度可以增大排氣量和延長壽命;但是過小的厚度會使閥片根部應力增加,閥片還會撞擊活塞,對結構造成損傷。

進氣閥片的受力狀況也是實際工程中非常關注的問題,尤其閥片撞擊升程限制器的那一刻,現(xiàn)截取閥片撞擊升程限制器那一刻的應力分布,如圖12所示。

圖12 進氣閥片初次撞擊升程限制器時的應力分布圖

從圖12中可以得到最大應力值為500 MPa,沒有超出失效率在5%時的許用應力值為680 MPa,閥片撞擊升程限制器部位的應力值并不大。

2.4 排氣閥運動規(guī)律

圖13 排氣閥片網(wǎng)格劃分圖中D點位置

圖14 排氣閥片升程-速度曲線圖

排氣閥片的運動規(guī)律對壓縮機的性能影響較大,因此有必要分析其運動特性以及應力特性。在排氣閥網(wǎng)格劃分圖上選取參考點D作為參考對象,如圖13所示。然后,對D點進行升程,并統(tǒng)計速度數(shù)據(jù)。圖14為排氣閥片升程速度曲線。由圖14可以看出,在排氣過程中,閥片波動較大,非常不穩(wěn)定,導致流動阻力損失變大,容積效率下降[15],并且會加劇氣流脈動,噪聲。由于排氣閥片的升程只有0.3 mm,因此會產(chǎn)生類似小孔噴注現(xiàn)象,激發(fā)高頻噪聲[16]。

圖15為排氣閥片最大升程處的應力云圖,從圖中可以得到最大應力值為815 MPa,超出了失效率在5%時的許用應力值680 MPa,因此閥片在長期的開閉工作狀態(tài)下容易斷裂。

圖15 排氣閥片最大升程處的應力云圖

3 結 論

本文采用ADINA軟件對變排量汽車空調(diào)壓縮機在最大排氣量的工況下的內(nèi)部工作過程和氣閥運動規(guī)律進行了三維流固耦合數(shù)值模擬,并將計算結果與性能試驗結果進行了比較驗證,計算結果的誤差在5%以內(nèi)。

根據(jù)本文對氣閥運動規(guī)律的分析得出如下結論。

(1)本文算例中,當進氣壓力為0.18 MPa、排氣壓力為1.5 MPa、轉(zhuǎn)速為2 000 r/min時,壓縮機的最大瞬時進氣阻力損失為0.069 MPa,平均進氣壓力損失為0.046 MPa,最大瞬時排氣損失為0.58 MPa,平均排氣壓力損失為0.11 MPa,壓縮機的容積效率為64.98%。通過數(shù)值模擬,可以看出目前汽車空調(diào)壓縮機的效率并不高。

(2)本文的研究發(fā)現(xiàn),閥片厚度不僅使阻力損失增大,還可能導致氣閥關閉延遲。本文算例中,當閥片厚度從0.281 mm增加時,閥片關閉角開始延遲,當厚度到達0.481 mm時,閥片關閉角從227.48°延遲到248.24°。這與通常的認識,即增加彈簧力(舌簧閥的閥片厚度)能減小延遲關閉相反。這是由于當進氣閥片彈性力增大后,閥片回落提前,導致進氣截面積減小,進氣量減少,從而進氣口和氣缸內(nèi)的壓力差會更大,活塞移動到上止點時,氣缸內(nèi)的壓力沒有達到進氣口的壓力值,因而在壓縮過程的初始階段,使用較厚閥片的氣缸內(nèi)壓力升高緩慢,在較大壓差的作用下導致進氣閥片延遲關閉現(xiàn)象更加嚴重。采用更厚的閥片,則導致出現(xiàn)顫振現(xiàn)象。

[1] 林梅, 吳業(yè)正. 壓縮機自動閥 [M]. 西安: 西安交通大學出版社, 1991: 2.

[2] 王楓. 活塞式制冷壓縮機的舌簧閥研究綜述 [J]. 流體機械, 2010, 39(1): 42-51. WANG Feng. Study on the reed valve of piston refrigeration compressor [J]. Fluid Machinery, 2010, 39(1): 42-51.

[3] GLEN P J R, WILLIAM P H. Optical measurement of discharge valve modal parameters for a rolling piston refrigeration compressor [J]. Measurement, 2003, 33(1): 75-84.

[4] BULIGAN G, PAONE N, REVEL G M, et al. Valve lift measurement by optical techniques in compressors [C]∥Sixteenth International Compressor Engineering Conference. Purdue, IN, USA: Purdue University, 2002: 5-13.

[5] LUDU A, BETTO A, REGNER G. Endoscope video of compressor valve motion and pressure measurement assist simulations for design improvement [C]∥Fifteenth International Compressor Engineering Conference. Purdue, IN, USA: Purdue University, 2000: 443-450.

[6] HABING R A. An experimental method for validating compressor valve vibration theory [J]. Journal of Fluids and Structures, 2006, 2(5): 683-679.

[7] 林梅, 姚峻, 潘樹林. 簧片閥三種力學模型的建立及比較 [J]. 應用力學學報, 1997, 14(1): 118-123. LIN Mei, YAO Jun, PAN Shulin. Establishment of three mathematical models and comparison for the reed valves [J]. Chinese Journal of Applied Mechanics, 1997, 14(1): 118-123.

[8] 鄭傳祥, 李蓉, 王亮, 等. 小型制冷壓縮機吸氣閥的流固耦合分析和優(yōu)化研究 [J]. 工程設計學報2014, 21(1): 69-74. ZHENG Chuanxiang, LI Rong, WANG Liang, et al. Fluid-structure interaction analysis and optimization study of suction valves in small refrigeration compressors [J]. Chinese Journal of Engineering Design, 2014, 21(1): 69-74.

[9] 李蓉. 制冷壓縮機過程模擬及節(jié)能技術開發(fā)研究 [D]. 杭州: 浙江大學, 2014.

[10]KIM J, WANG S, PARK S, et al. Valve dynamic analysis of a hermetic reciprocating compressor [C]∥Eighteenth International Compressor Engineering Conference. Purdue, IN, USA: Purdue University, 2006: 1-9.

[11]CHOI Y S, LEE J H, JEONG W B, et al. Dynamic behavior of valve system in linear compressor based on fluid-structure interaction [J]. Journal of Mechanical Science and Technology, 2010, 24(7): 1371-1377.

[12]Sandvik Materials Technology Company. Sandvik 20C flapper valve steel. [EB/OL]. [2012-10-31]. http:∥www.doc88.com/p-1813461662795.html.

[13]岳戈, 陳晨, 馬野, 等. ADINA流體與流固耦合功能的高級應用 [M]. 北京: 人民交通出版社, 2010: 61-63.

[14]ADINA Research and Development, Inc. Theory and modeling guide [M]. Watertown, USA: ADINA Research and Development, Inc. 2005: 8.

[15]NAGATA S, NOZAKI T, AKIZAWA T. Analysis of dynamic behavior of suction valve using strain gauge in reciprocating compressor [C]∥Twentieth International Compressor Engineering Conference. Purdue, IN, USA: Purdue University, 2010: 1273.

[16]仇穎, 李紅旗, 呂亞東. 應用小孔噴注噪聲原理分析全封閉活塞式制冷壓縮機吸排氣閥噪聲 [J]. 流體機械, 2007, 35(1): 25-28. QIU Ying, LI Hongqi, Lü Yadong. Analysis of suction and discharge valve noise in a hermetic reciprocating compressor by using jet noise theory [J]. Fluid Machinery, 2007, 35(1): 25-28.

(編輯 劉楊 苗凌)

Numerical Simulation of the Valve Motion in Automotive Air Conditioning Compressors

SONG Mingyi,WU Weifeng,LI Zhi

(School of Energy and Power Engineering, Xi’an Jiaotong University, Xi’an 710049, China)

A three-dimensional fluid-structure interaction model for the working process of the automotive air conditioning compressors is proposed and its internal work process and the motion of the valve plate are analyzed to solve the problems of valve fracture in automotive air conditioning compressors. Experiments on an automotive air conditioning compressor are performed to validate the model, and the results show a reasonable agreement and that the average loss of the suction pressure is 0.046 MPa, the average loss of the discharge pressure is 0.11 MPa, and the volumetric efficiency is 64.98% when the suction pressure is 0.18 MPa, the discharge pressure is 1.5 MPa and the rotation speed is 2 000 r/min. Take the delay shut-off of valve or not as the selection principle, the results also show that an optimal thickness of the suction valve plate is 0.281 mm, and when the thickness of the valve plate exceeds the optimal value, the shut-off time of the valve is delayed with the increasing of the thickness. These results are opposite to the traditional theory that increasing the thickness of the valve leads to forward shut-off and provide a guidance to the design and optimization of reed valves used in small refrigeration compressors.

automotive air conditioning compressor; three-dimensional fluid-structure interaction; valve motion

2015-02-02。

宋明毅(1989—),男,碩士生;吳偉烽(通信作者),男,副教授,博士生導師。

時間:2015-10-26

10.7652/xjtuxb201512023

TB652;U270.38

A

0253-987X(2015)12-0144-07

網(wǎng)絡出版地址:http:∥www.cnki.net/kcms/detail/61.1069.T.20151026.2044.002.html