(山東科技大學(xué) 機械電子工程學(xué)院,山東 青島 266590)
近年來,冰箱制冷壓縮機市場發(fā)展迅速[1],在全球范圍內(nèi),約有99%家用冰箱使用活塞式制冷壓縮機[2]?;钊綁嚎s機具有結(jié)構(gòu)簡單、易批量生產(chǎn)和成本低等優(yōu)點,然而也存在一些尚未解決的問題,比如運行過程中會產(chǎn)生較大的噪聲和機械振動[3]。壓縮機的核心部件之一是曲柄連桿機構(gòu),其不平衡力是產(chǎn)生振動噪聲的根源。很多學(xué)者對曲柄連桿機構(gòu)產(chǎn)生的不平衡慣性力進行了研究,如康獻民等[4]利用主軸傾角和不平衡率的研究來調(diào)整最大慣性力。王秋曉等[5]通過分析偏轉(zhuǎn)角度來研究不平衡率和一階、二階慣性力。Pont等[6]探討了曲軸質(zhì)量偏心角對曲軸穩(wěn)定性的影響和偏移對壓縮機行為的影響,閔劍青等[7]對曲軸質(zhì)心的位置對慣性力平衡的影響進行研究。
很多學(xué)者運用多體動力學(xué)軟件Adams來研究曲柄連桿機構(gòu)。吳迪等[8]和王金友等[9]通過該軟件建立了曲柄連桿機構(gòu)的動力學(xué)仿真模型,分析曲柄連桿機構(gòu)的運行情況。張小明等[10]和張冬香等[11]運用Adams軟件得到慣性力與平衡重參數(shù)的變化曲線,通過優(yōu)化平衡重的參數(shù),可以減小干擾力,實現(xiàn)慣性力的平衡。蔣立軍[12]運用該軟件研究了偏心距和曲柄連桿機構(gòu)產(chǎn)生的干擾力和力矩的關(guān)系,表明改變偏心距可以減小往復(fù)慣性力。夏江敏等[13]對曲柄連桿機構(gòu)進行振動分析,得到加速度頻率響應(yīng)曲線。這些研究運用多體動力學(xué)軟件,通過對影響參數(shù)的優(yōu)化,減小了慣性力。
上述研究以減少曲柄連桿機構(gòu)受到的慣性力為目的來調(diào)整機構(gòu)和優(yōu)化參數(shù),但這些參數(shù)對壓縮機振動的影響尚未有進一步研究。本研究以某型號壓縮機為研究對象,運用理論分析和多體動力學(xué)仿真方法,研究曲柄連桿結(jié)構(gòu)參數(shù)對壓縮機振動的影響。實驗驗證了優(yōu)化后的參數(shù)能夠降低壓縮機的振動。
某型號壓縮機為單缸臥式壓縮機,其曲柄連桿機構(gòu)的主軸為直立布置,采用偏心軸形式,氣體只在活塞一側(cè)進行壓縮,簡化后的物理模型[14]如圖1所示。
圖1 曲柄連桿機構(gòu)動力學(xué)模型Fig.1 Dynamic model of crank-connecting rod mechanism
設(shè)曲柄OA繞曲軸中心線O以角速度ω作逆時針轉(zhuǎn)動,曲柄半徑OB=r,連桿長L,曲軸偏置e,平衡塊的質(zhì)心到曲軸中心線O的距離為ρ,連桿與氣缸中心線的夾角為連桿擺動角β,曲柄轉(zhuǎn)角α。
設(shè)曲柄連桿機構(gòu)旋轉(zhuǎn)質(zhì)量為mr,其包括連桿大頭等效質(zhì)量和曲軸銷質(zhì)量,旋轉(zhuǎn)慣性力Fr,質(zhì)徑積為Pr=mrr;往復(fù)質(zhì)量為mj,其包括連桿小頭等效質(zhì)量、活塞銷質(zhì)量和活塞質(zhì)量,一階往復(fù)慣性力Fj,質(zhì)徑積為Pj=mjr;平衡塊的質(zhì)量為mw,質(zhì)徑積為Pw=mwρ。
根據(jù)牛頓第二定律,在x,y軸上建立平衡方程[15-16],有
(1)
帶入質(zhì)徑積得
(2)
簡化后得
(3)
合成后,剩余慣性力為
(4)
根據(jù)對曲柄連桿機構(gòu)的受力分析,曲柄連桿機構(gòu)所受的往復(fù)慣性力Fj通過連桿和曲柄傳遞到壓縮機的主軸承座上,使壓縮機產(chǎn)生垂向振動,而且在傳遞過程中,會產(chǎn)生傾覆力矩Mj,使壓縮機產(chǎn)生扭轉(zhuǎn)振動。
Mj=Fjtanβ(lcosβ+rcosα)=Pjω2tanβ(lcosβ+rcosα)。
(5)
壓縮機的旋轉(zhuǎn)慣性力Fr由曲柄傳至主軸承座,使壓縮機產(chǎn)生垂向和橫向振動。因此,當(dāng)壓縮機運轉(zhuǎn)時,這些慣性力是引起壓縮機振動的根源[17-18]。減少慣性力是降低壓縮機振動的重要途徑,從式(4)可以看出,剩余慣性力會隨著曲軸的旋轉(zhuǎn)而周期性變化,所以通過改變影響參數(shù),可以減少剩余慣性力。
該壓縮機曲柄半徑r=10.3 mm,連桿長L=41 mm,曲軸偏置e=2.95 mm,旋轉(zhuǎn)質(zhì)徑積Pr=623.665 g·mm,往復(fù)質(zhì)徑積Pj=474.83 g·mm,平衡塊質(zhì)徑積Pw=906.93 g·mm??紤]到壓縮機的空間尺寸,從式(4)可以看出,在影響剩余慣性力的參數(shù)中,曲軸和連桿的尺寸難以改變,因此對平衡塊質(zhì)徑積Pw和曲軸偏置e進行分析。
根據(jù)式(1)分析可得,平衡塊能夠完全平衡旋轉(zhuǎn)慣性力和部分往復(fù)慣性力,所以平衡塊的質(zhì)徑積在旋轉(zhuǎn)慣性質(zhì)徑積與旋轉(zhuǎn)、往復(fù)質(zhì)徑積的和之間,即Pw在Pr到Pr+Pj之間。對Pw進行賦值,在Matlab軟件上進行分析[19],得到剩余慣性力F和Pw的關(guān)系,如圖2所示。
由圖2可以看出,給定平衡塊的參數(shù)后,剩余慣性力隨著曲軸轉(zhuǎn)角的變化呈正弦曲線變化。但是,當(dāng)Pw取值范圍在750~950 g·mm時,剩余慣性力變化較小。因此在這一區(qū)間,可以通過仿真分析,取得Pw最優(yōu)值。
該型號壓縮機的曲軸回轉(zhuǎn)中心線與氣缸中心線不相交,因此為偏心曲柄連桿機構(gòu)。根據(jù)偏心方向不同,分為正偏心機構(gòu)和負偏心機構(gòu)。該壓縮機采用正偏心機構(gòu),在活塞下行時連桿擺角較小,使得做功行程中活塞側(cè)推力有所減小。
根據(jù)式(1),對曲軸偏置e進行賦值,運用Matlab軟件進行數(shù)值計算,得到剩余慣性力和曲軸偏置e的關(guān)系,如圖3所示。
圖2 剩余慣性力F與Pw的關(guān)系Fig.2 Relationship between residual inertial forceFandPw
圖3 剩余慣性力F與e的關(guān)系Fig.3 Relationship between residual inertia forceFande
從圖3中可以看出,給定曲軸偏置的數(shù)值后,剩余慣性力隨著曲軸轉(zhuǎn)角的變化呈正弦曲線,在曲軸偏置數(shù)值緩慢增大過程中,正弦曲線的峰值也在逐漸增大,即剩余慣性力也隨之增大。
圖4 壓縮機三維實體模型Fig.4 3D solid model of compressor
利用Pro/Engineer三維造型軟件建立壓縮機實體三維模型[20],將模型導(dǎo)入到ADAMS/VIEW中。底座的彈簧與大地連接,曲軸與電機同步轉(zhuǎn)動,添加固定副,平衡塊用固定副安裝在曲軸上,曲軸和連桿大頭用點-線副連接,連桿大頭和活塞銷以及活塞銷和活塞的連接處添加旋轉(zhuǎn)副,活塞和氣缸之間用滑移副連接[21],仿真模型如圖4所示。
對壓縮機多體動力學(xué)模型進行仿真,在活塞運動一側(cè)添加一個最大振幅為10,相位角為0,水平方向的正弦激振力作為輸入。在壓縮機底部電機與外界連接的4個彈簧處添加4個垂直方向的輸出通道,通過4個輸出通道的頻率響應(yīng),分析剩余慣性力參數(shù)對壓縮機振動的影響。
對壓縮機的平衡塊質(zhì)徑積參數(shù)的理論分析可知,當(dāng)Pw在750~950 g·mm之間取值時,剩余慣性力變化較小。因此,對Pw間隔10 g·mm進行取值,分別進行振動仿真分析,得到壓縮機內(nèi)部的振動傳遞到殼體的頻率響應(yīng)。圖5為Pw影響的加速度幅頻響應(yīng)曲線。
圖5 Pw影響的加速度幅頻響應(yīng)曲線Fig.5 Acceleration amplitude-frequency response curve affected byPw
圖5每個通道加速度幅頻曲線的4條曲線分別是:當(dāng)平衡塊質(zhì)徑積Pw的取值為當(dāng)前參數(shù)906 g·mm時,頻率響應(yīng)的曲線為l1;當(dāng)Pw的取值為740、770、790、870和940 g·mm時,頻率響應(yīng)的曲線均為l2;當(dāng)Pw的取值為760、780、820、840、860、920和930 g·mm,頻率響應(yīng)的曲線均為l3;當(dāng)取Pw其余值時,頻率響應(yīng)的曲線為l4。
從圖5可以看出,平衡塊質(zhì)徑積Pw對壓縮機的振動產(chǎn)生較大影響,改變平衡塊質(zhì)徑積參數(shù)可以改變壓縮機的加速度頻率響應(yīng)。當(dāng)平衡塊質(zhì)徑積為曲線l3代表的參數(shù)時,在第二和第三通道比原有參數(shù)產(chǎn)生的加速度幅值有明顯降低。但在Pw取值范圍750~950 g·mm中,有些參數(shù)取值可以減小剩余慣性力,但壓縮機的振動并未降低。因此,通過振動仿真分析,優(yōu)化平衡塊質(zhì)徑積參數(shù),可以更有效地降低壓縮機振動。
通過對壓縮機實際測量,得到曲軸偏置為2.95 mm。為分析曲軸偏置對振動的影響,在0~3.5 mm選取7個參數(shù),e1=0.45 mm,e2=0.95 mm,e3=1.45 mm,e4=1.95 mm,e5=2.45 mm,e6=2.95 mm,e7=3.45 mm,分別對每個參數(shù)進行振動仿真分析。壓縮機底部的加速度幅頻響應(yīng)曲線如圖6所示。
圖6 曲軸偏置e影響的加速度幅頻響應(yīng)曲線Fig.6 Acceleration amplitude-frequency response curve affected by crankshaft offsete
從圖6中可以看出,曲軸偏置對壓縮機的振動影響較大。當(dāng)曲軸偏置為e2時,加速度幅值最小。當(dāng)曲軸偏置為e4時,壓縮機的加速度幅值與原偏置e6接近相同。當(dāng)曲軸偏置為e3時,除在第3個通道產(chǎn)生的加速度比原來的偏置稍大外,在其他輸出通道,加速度幅值均有較大幅度降低。當(dāng)曲軸偏置1.95和1.45 mm時,加速度幅值可以降至50 dB,在第一通道降到20 dB,有效減小了壓縮機的振動。
為驗證仿真結(jié)果,在半消聲室內(nèi)對壓縮機振動進行測試,如圖7所示。進氣閥口、排氣閥口分別與進氣管道、出氣管道連接,3個加速度傳感器分別放置在壓縮機的X、Y和Z軸,壓縮機正常運轉(zhuǎn),測量壓縮機殼體的振動。
圖8(a)為平衡塊質(zhì)徑積對壓縮機振動影響的加速度幅頻響應(yīng)曲線,L1為3.2中仿真優(yōu)化后的平衡塊質(zhì)徑積參數(shù),L2為壓縮機原有平衡塊質(zhì)徑積參數(shù)。圖8(b)為曲軸偏置對壓縮機振動影響的幅頻響應(yīng)曲線,e1為3.3中仿真優(yōu)化后的曲軸偏置參數(shù),e2為壓縮機原有的曲軸偏置參數(shù)。
圖7 壓縮機殼體振動實驗Fig.7 Compressor shell vibration experiment
圖8 實驗測得加速度幅頻響應(yīng)曲線Fig.8 Experimentally measured acceleration amplitude-frequency response curve
從圖8(a)可以看出,優(yōu)化后的幅頻響應(yīng)曲線更趨平緩,加速度幅值降低。在頻率為50 Hz時,加速度幅值由112 dB減小為105 dB。從圖8(b)可以看出,曲軸偏置經(jīng)優(yōu)化后,加速度幅值明顯降低,平均降幅約為5 dB。
本研究采用多體動力學(xué)仿真和實驗相結(jié)合的方法,以減小壓縮機振動加速度幅頻響應(yīng)為目標,對壓縮機曲柄連桿機構(gòu)的平衡塊質(zhì)徑積和曲軸偏置參數(shù)進行了分析。與單純依靠減小不平衡慣性力來降低壓縮機振動噪聲的傳統(tǒng)方法相比,該方法通過優(yōu)化曲柄連桿機構(gòu)的結(jié)構(gòu)參數(shù),能夠更加有效地降低壓縮機的振動與噪聲,實驗結(jié)果驗證了該方法的可行性。該方法為冰箱壓縮機曲柄連桿結(jié)構(gòu)的參數(shù)設(shè)計,降低壓縮機的振動噪聲,提供了可行性的途徑。