（蘭州理工大學(xué) 石油化工學(xué)院，蘭州 730050）

0 引言

渦旋壓縮機(jī)相比于其他類型的壓縮機(jī)具有高效、低振和低噪等諸多優(yōu)點(diǎn)而廣泛應(yīng)用于動(dòng)力工程、交通運(yùn)輸以及空調(diào)制冷等領(lǐng)域。作為渦旋壓縮機(jī)的主要部件，轉(zhuǎn)子系統(tǒng)的工作狀態(tài)直接影響著渦旋壓縮機(jī)的穩(wěn)定性。制造裝配誤差、材料變形及磨損等問(wèn)題造成的運(yùn)動(dòng)副間隙，必然導(dǎo)致渦旋壓縮機(jī)的動(dòng)不平衡。在渦旋壓縮機(jī)運(yùn)行期間，含間隙運(yùn)動(dòng)副元素會(huì)發(fā)生接觸失真的現(xiàn)象，接觸時(shí)產(chǎn)生碰撞，引起非線性沖擊載荷，造成運(yùn)動(dòng)副的破壞以及影響壓縮機(jī)系統(tǒng)載荷的正常傳遞，降低了渦旋壓縮機(jī)的工作效率和可靠性。

國(guó)內(nèi)外已有諸多針對(duì)渦旋壓縮機(jī)轉(zhuǎn)子系統(tǒng)的相關(guān)研究，李超等［1］針對(duì)變?nèi)萘繙u旋壓縮機(jī)曲軸系統(tǒng)，運(yùn)用多體動(dòng)力學(xué)和有限元方法，計(jì)算獲得了曲軸各軸承載荷及曲軸應(yīng)力應(yīng)變隨壓縮機(jī)轉(zhuǎn)速的變化規(guī)律，并對(duì)曲軸最大應(yīng)力應(yīng)變位置進(jìn)行了安全性校核；邱海飛［2］通過(guò)建立的渦旋壓縮機(jī)主軸傳動(dòng)系統(tǒng)等效力學(xué)模型，應(yīng)用參數(shù)化建模與動(dòng)態(tài)優(yōu)化方法，獲得主軸系統(tǒng)自振頻率和振型，為改進(jìn)主軸的動(dòng)平衡性能，對(duì)主軸平衡重進(jìn)行了優(yōu)化設(shè)計(jì)；余洋等［3］考慮曲軸軸徑處存在運(yùn)動(dòng)副間隙條件下的動(dòng)渦旋盤動(dòng)力特性，分析了動(dòng)渦旋盤的運(yùn)動(dòng)軌跡以及位移變化規(guī)律，為運(yùn)動(dòng)副間隙影響的渦旋壓縮機(jī)動(dòng)平衡設(shè)計(jì)提供了參考；黃華軍等［4］針對(duì)小曲拐防自轉(zhuǎn)機(jī)構(gòu)與動(dòng)渦旋盤連接處存在運(yùn)動(dòng)副間隙的情況，分析了間隙數(shù)目與間隙大小對(duì)動(dòng)渦旋盤及小曲拐接觸力的影響規(guī)律，總結(jié)了間隙位置的不同對(duì)渦旋壓縮機(jī)動(dòng)平衡的影響。

上述研究分析了渦旋壓縮機(jī)轉(zhuǎn)子系統(tǒng)的振動(dòng)、軸承受力以及運(yùn)動(dòng)副間隙對(duì)動(dòng)渦旋盤和防自轉(zhuǎn)機(jī)構(gòu)小曲拐運(yùn)動(dòng)特性的影響。但有關(guān)運(yùn)動(dòng)副間隙對(duì)曲軸動(dòng)態(tài)特性的影響、渦旋壓縮機(jī)含間隙運(yùn)動(dòng)副對(duì)其瞬態(tài)動(dòng)力性能的影響等方面的研究較少。本文以十字滑環(huán)為防自轉(zhuǎn)機(jī)構(gòu)的渦旋壓縮機(jī)轉(zhuǎn)子系統(tǒng)為研究對(duì)象，在考慮運(yùn)動(dòng)副間隙的作用時(shí)采用非線性彈簧阻尼模型和庫(kù)侖摩擦力模型，基于多體動(dòng)力學(xué)分析軟件，并結(jié)合有限元方法，詳細(xì)分析了運(yùn)動(dòng)副間隙下轉(zhuǎn)子系統(tǒng)的瞬態(tài)動(dòng)力特性，可為渦旋壓縮機(jī)加工制造和裝配以及提高轉(zhuǎn)子系統(tǒng)動(dòng)力特性提供理論依據(jù)。

1 渦旋壓縮機(jī)轉(zhuǎn)子系統(tǒng)模型

轉(zhuǎn)子系統(tǒng)的組成如圖1所示，曲軸由主、副軸承支承，曲柄銷處由滾針軸承將曲軸與動(dòng)渦旋盤連接，在運(yùn)行過(guò)程中，電機(jī)通過(guò)帶輪帶動(dòng)曲軸驅(qū)動(dòng)動(dòng)渦旋盤作公轉(zhuǎn)平動(dòng)，實(shí)現(xiàn)與靜渦旋盤的嚙合，完成對(duì)氣體的吸氣、壓縮及排氣過(guò)程［5］。

圖1 渦旋壓縮機(jī)轉(zhuǎn)子系統(tǒng)三維模型Fig.1 Three-dimensional model of rotor system for scroll compressor

2 運(yùn)動(dòng)副間隙力學(xué)建模

在渦旋壓縮機(jī)運(yùn)行期間，運(yùn)動(dòng)副間隙的存在使得軸與軸承之間的理想接觸產(chǎn)生了松動(dòng)，造成接觸失真現(xiàn)象，軸與軸承之間從接觸狀態(tài)到自由運(yùn)動(dòng)狀態(tài)，再到接觸狀態(tài)，整個(gè)過(guò)程包含著接觸與碰撞。旋轉(zhuǎn)副間隙中含有潤(rùn)滑油，可將這種接觸考慮為彈簧阻尼模型，根據(jù)文獻(xiàn)［4］與文獻(xiàn)［9］建立接觸碰撞模型。

2.1 運(yùn)動(dòng)副間隙的矢量模型

如圖2所示，在平面旋轉(zhuǎn)鉸內(nèi)引入一個(gè)間隙矢量來(lái)表達(dá)間隙矢量模型［6-8］。

圖2 間隙矢量模型示意Fig.2 Schematic diagram of joint clearance vector model

間隙矢量可以反映運(yùn)動(dòng)副間隙元素之間的相對(duì)位置，可有效處理間隙運(yùn)動(dòng)副連接點(diǎn)相對(duì)位置的變化，間隙大小用軸承半徑與曲柄銷半徑之差來(lái)描述，則：

式中 c——間隙；

rj——軸承半徑；

ri——曲柄銷軸半徑。

定義XOY為廣義坐標(biāo)系，間隙矢量定義在一個(gè)局部浮動(dòng)笛卡爾坐標(biāo)系xjojyj中，以軸承的回轉(zhuǎn)中心為間隙矢量的基準(zhǔn)起始點(diǎn)，矢量方向指向曲柄銷和軸承相對(duì)運(yùn)動(dòng)的接觸點(diǎn)，終點(diǎn)為曲柄銷的回轉(zhuǎn)中心。由圖可得，曲柄銷中心與驅(qū)動(dòng)軸承中心的間隙矢量eij可表示為：

間隙矢量的大小可表示為：

式中 ex，ey——軸承中心到曲柄銷中心距離在X和Y方向的分量。

曲柄銷與驅(qū)動(dòng)軸承之間的穿透深度δ為：

2.2 運(yùn)動(dòng)副間隙的接觸碰撞力模型

接觸碰撞模型將碰撞過(guò)程歸結(jié)為“自由運(yùn)動(dòng)-接觸變形”兩種狀態(tài)的隨機(jī)變換，通過(guò)考慮接觸表面的阻尼和彈性因素，比較真實(shí)地描述了碰撞過(guò)程。本文采用的是一種間隙非線性連續(xù)接觸碰撞力模型［9］，是基于 Lankarani-Nikravesh 接觸模型與彈性模型的融合。其表達(dá)式為：

式中 Fn——接觸碰撞力；

K——非線性剛度系數(shù)；

δ——接觸點(diǎn)的法向穿透距離；

n——冪指數(shù)，取決于碰撞材料的種類，一般金屬材料為1.5；

C——阻尼系數(shù)；

E*——接觸表面的彈性模量；

μi，μj——2 種材料的泊松比；

Ei，Ej——2種材料的楊氏模量；

Cr——恢復(fù)系數(shù)，代表能量耗散，其值小于1，為碰撞前后的速度之比。

2.3 運(yùn)動(dòng)副間隙的摩擦力模型

為了進(jìn)一步真實(shí)地描述運(yùn)動(dòng)副間隙，必須考慮接觸時(shí)的切向摩擦力，采用修正的庫(kù)侖摩擦模型［9］，切向摩擦力在可以表示為：

式中 vt——接觸點(diǎn)相對(duì)滑移速度的切向分量。

此模型中的摩擦因數(shù)并不是一個(gè)常數(shù)，而是一個(gè)關(guān)于切向滑移速度的動(dòng)態(tài)摩擦函數(shù)，表達(dá)式為：

式中 vs——靜摩擦臨界速度；

vd——最大動(dòng)摩擦臨界速度；

μd——?jiǎng)幽Σ料禂?shù)；

μs——靜摩擦系數(shù)。

2.4 ADAMS接觸碰撞力函數(shù)

ADAMS軟件內(nèi)置的接觸碰撞力函數(shù)是IMPACT函數(shù)，其采用非線性彈簧阻尼器來(lái)計(jì)算接觸碰撞力，IMPACT接觸碰撞力的計(jì)算公式為：

Cmax——最大阻尼系數(shù)；

dmax——在最大阻尼系數(shù)前提下的穿透邊界深度。

3 渦旋壓縮機(jī)轉(zhuǎn)子系統(tǒng)動(dòng)力學(xué)分析

3.1 建立轉(zhuǎn)子系統(tǒng)動(dòng)力學(xué)模型

根據(jù)表1，2所示的轉(zhuǎn)子系統(tǒng)參數(shù)，通過(guò)SOLIDWORKS建立渦旋壓縮機(jī)轉(zhuǎn)子系統(tǒng)三維模型，驅(qū)動(dòng)軸承材料為GCr15，密度為7.83 g/mm3，彈性模量219 GPa，泊松比為0.3。其型號(hào)為IS4215，外徑62 mm，內(nèi)徑50 mm，滾針個(gè)數(shù)為55，滾針直徑3.5 mm，長(zhǎng)度為26.25 mm。然后導(dǎo)入ADAMS軟件中建立動(dòng)力學(xué)模型，模型如圖3所示。

表1 型線及渦旋盤參數(shù)Tab.1 Profile and parameters of scroll member

表2 轉(zhuǎn)子系統(tǒng)各部件參數(shù)Tab.2 Parameters of each component of the rotor system

圖3 轉(zhuǎn)子系統(tǒng)動(dòng)力學(xué)模型Fig.3 Dynamic model of rotor system

3.2 施加載荷及約束

根據(jù)運(yùn)動(dòng)副間隙對(duì)渦旋壓縮機(jī)動(dòng)平衡的影響，所建立的模型中曲柄銷和驅(qū)動(dòng)軸承連接處采用接觸副連接，同時(shí)設(shè)置庫(kù)侖摩擦力參數(shù)。根據(jù)轉(zhuǎn)子的真實(shí)運(yùn)動(dòng)情況在VIEW中對(duì)各個(gè)部件施加約束以及旋轉(zhuǎn)驅(qū)動(dòng)并設(shè)置相應(yīng)的轉(zhuǎn)速函數(shù)。

3.3 動(dòng)力學(xué)仿真結(jié)果與分析

由于曲軸與驅(qū)動(dòng)軸承之間存在間隙，兩者之間存在碰撞，引起轉(zhuǎn)子系統(tǒng)的非理想狀態(tài)運(yùn)行，因此針對(duì)驅(qū)動(dòng)軸承和曲軸進(jìn)行測(cè)量。設(shè)置4個(gè)間隙梯度，分別為0（理想間隙），0.01，0.03，0.05 mm。

圖4（a）（b）分別示出了曲軸的質(zhì)心軌跡及局部放大，由圖可知理想情況下曲軸的軌跡為一標(biāo)準(zhǔn)圓形，當(dāng)存在間隙時(shí)質(zhì)心軌跡出現(xiàn)了折點(diǎn)，而且間隙越大同一個(gè)軌跡中折點(diǎn)之間的距離會(huì)相應(yīng)增加，曲軸的質(zhì)心軌跡相對(duì)于理想情況的偏移越大。

圖4 曲軸位移Fig.4 Crankshaft displacement

圖5（a）（b）分別示出了驅(qū)動(dòng)軸承分別在水平與垂直方向上的速度變化規(guī)律，理想情況時(shí)速度曲線為光滑的正余弦曲線，存在間隙時(shí)速度曲線沿著理想情況下的曲線上下波動(dòng)，隨著間隙的增加速度曲線波動(dòng)的峰值增大，但波動(dòng)頻率減小。

圖6所示（a）（b）分別示出了驅(qū)動(dòng)軸承在水平與垂直方向上的加速度變化，間隙為0時(shí)加速度曲線為正余弦曲線，間隙情況下的曲線沿著理想情況曲線上下劇烈的波動(dòng)，并且曲線峰值隨著間隙的增加而增大。

圖5 驅(qū)動(dòng)軸承速度Fig.5 Driving bearing velocity

圖6 驅(qū)動(dòng)軸承加速度Fig.6 Driving bearing acceleration

圖7示出了運(yùn)動(dòng)副間隙處的接觸力隨時(shí)間的變化規(guī)律，理想情況下接觸力約為2 200 N，間隙條件下的接觸力出現(xiàn)了高頻率的波動(dòng)，隨著間隙的增大，接觸力的瞬時(shí)峰值也逐漸增大，間隙值為0.05 mm時(shí)，曲軸與驅(qū)動(dòng)軸承之間的接觸力最大可達(dá)8 000 N左右，隨著轉(zhuǎn)動(dòng)周期的增加，接觸力逐漸減小并在4 000 N的范圍內(nèi)上下波動(dòng)。

圖7 曲柄銷與驅(qū)動(dòng)軸承之間的接觸力Fig.7 Contact force between crank pin and driving bearing

從圖8可知，不含間隙時(shí)驅(qū)動(dòng)軸承角速度約為303 rad/s，第一個(gè)轉(zhuǎn)動(dòng)周期內(nèi)，角速度呈現(xiàn)大幅度的波動(dòng)，在間隙值為0.05 mm時(shí)，最大值約為304.4 rad/s，隨后波動(dòng)峰值逐漸減小，并沿著303 rad/s上下波動(dòng)。

圖8 驅(qū)動(dòng)軸承角速度Fig.8 Driving bearing angular velocity

圖9示出了驅(qū)動(dòng)軸承角加速度隨時(shí)間的變化，理想情況下驅(qū)動(dòng)軸承角加速度為零，間隙的存在對(duì)角加速度的影響十分敏感，角加速度的變化十分劇烈且瞬時(shí)峰值較大。在間隙值為0.05 mm時(shí)，穩(wěn)定后的角加速度在3500 rad/s2內(nèi)波動(dòng)。

圖9 驅(qū)動(dòng)軸承角加速度Fig.9 Angular acceleration of driving bearing

間隙接觸力、位移、速度、加速度、角速度以及角加速度的值均會(huì)隨著間隙的增加而增大，但其波動(dòng)的頻率均會(huì)隨著間隙的增加而減小。這是由于間隙增加時(shí)，驅(qū)動(dòng)軸承與曲軸之間從接觸狀態(tài)到自由運(yùn)動(dòng)狀態(tài)的時(shí)間也相應(yīng)地增加，接觸力變化的頻率也就越小，但接觸力的大小將變大，而力的增大導(dǎo)致了速度、加速度、角速度和角加速度等也隨之增大。速度增大將導(dǎo)致曲軸與驅(qū)動(dòng)軸承再次接觸碰撞前構(gòu)件動(dòng)能增加，碰撞力增大；加速度越大，相同時(shí)間內(nèi)構(gòu)件速度的增速就越大，曲軸與驅(qū)動(dòng)軸承再次碰撞時(shí)的碰撞力就越大；該過(guò)程是連續(xù)交替循環(huán)的，會(huì)形成惡性循環(huán)，使構(gòu)件之間的磨損加劇，間隙再次增大，因此，需要對(duì)渦旋壓縮機(jī)進(jìn)行動(dòng)平衡以控制間隙處于合理的范圍，以防降低渦旋壓縮機(jī)的效率甚至失效。

4 渦旋壓縮機(jī)轉(zhuǎn)子系統(tǒng)瞬態(tài)特性分析

瞬態(tài)動(dòng)力學(xué)分析是對(duì)隨時(shí)間任意變化的載荷作用下，結(jié)構(gòu)動(dòng)力響應(yīng)過(guò)程的分析，動(dòng)力學(xué)通用方程如下：

式中 ［M］——質(zhì)量矩陣；

［C］——阻尼矩陣；

［K］——?jiǎng)偠染仃嚕?/p>

4.1 轉(zhuǎn)子系統(tǒng)有限元模型建立

將轉(zhuǎn)子系統(tǒng)三維模型導(dǎo)入到軟件中，進(jìn)行材料屬性的定義，然后選擇合適的網(wǎng)格進(jìn)行劃分，得到371 470個(gè)節(jié)點(diǎn)以及83 767個(gè)單元，如圖10所示。

圖10 轉(zhuǎn)子系統(tǒng)網(wǎng)格劃分Fig.10 Mesh division of the rotor system

4.2 施加載荷及約束

根據(jù)轉(zhuǎn)子系統(tǒng)的運(yùn)動(dòng)特點(diǎn)，選擇合理的載荷及約束是模擬的必要條件。首先，由于主、副軸承起支撐作用，本文采用軸承單元進(jìn)行代替，如圖11所示。

圖11 軸承單元簡(jiǎn)化Fig.11 Simplified diagram of bearing unit

曲柄銷處施加遠(yuǎn)端約束，曲軸尾端帶輪處施加柱面約束，均限制軸向位移。然后在曲柄銷處設(shè)置角加速度約束，并將獲得的驅(qū)動(dòng)軸承角加速度數(shù)據(jù)輸入，以此作為瞬態(tài)分析的邊界條件對(duì)曲軸進(jìn)行分析；設(shè)置轉(zhuǎn)速范圍，輸入不同的間隙下的剛度系數(shù)，以此對(duì)曲軸進(jìn)行臨界轉(zhuǎn)速計(jì)算分析。

4.3 有限元仿真結(jié)果與分析

通過(guò)WORKBENCH軟件計(jì)算分析可獲得不同間隙值下應(yīng)力分布情況，如圖12所示。

在曲軸與動(dòng)渦旋盤連接的曲拐處、曲軸與帶輪連接的軸肩位置均發(fā)生了應(yīng)力集中現(xiàn)象，并且曲軸與動(dòng)渦旋盤連接的曲拐處的應(yīng)力最大，這是因?yàn)榍N與主軸連接處軸徑突變和曲柄銷偏心導(dǎo)致曲柄銷旋轉(zhuǎn)離心力增加而引起應(yīng)力集中。圖12（a）所示，間隙為0.01 mm時(shí)曲軸的等效應(yīng)力為25.3 MPa；圖12（b）所示，間隙為0.03 mm時(shí)曲軸的等效應(yīng)力為30.6 MPa；圖12（c）所示，間隙為0.05 mm時(shí)曲軸的等效應(yīng)力為36.4 MPa。隨著間隙的增大，發(fā)生應(yīng)力集中部位的應(yīng)力也相應(yīng)增大，這是因?yàn)椴煌g隙下曲軸與驅(qū)動(dòng)軸承之間的接觸力載荷發(fā)生了變化，接觸力的大小隨間隙增大而增大，使應(yīng)力集中現(xiàn)象更加嚴(yán)重，而且曲軸的受力過(guò)程是連續(xù)交替循環(huán)的，這會(huì)給曲軸的安全運(yùn)行帶來(lái)更大的隱患。

圖12 不同間隙下曲軸應(yīng)力分布云圖Fig.12 Crankshaft stress distribution nephogram with different clearances

間隙為 0.01，0.03，0.05 mm 時(shí)，曲柄銷與驅(qū)動(dòng)軸承的接觸剛度分別為：2.90×107，2.08×107，1.35×107N/m。經(jīng)過(guò)WORKBENCH軟件計(jì)算分析可畫(huà)出不同間隙下曲軸系統(tǒng)的坎貝爾圖，如圖13所示分別為0.01，0.03，0.05 mm間隙時(shí)曲軸的坎貝爾圖，圖中的三角形分別表示不同間隙下曲軸系統(tǒng)的臨界轉(zhuǎn)速，由圖可知同一間隙下該曲軸有2個(gè)臨界轉(zhuǎn)速，而不同間隙情況下的第一階臨界轉(zhuǎn)速均小于9.13 rad/s，對(duì)應(yīng)的第二階臨界轉(zhuǎn)速分別為 1 447.4，1 444.3，1 438.6 rad/s。

圖13 不同間隙下曲軸坎貝爾圖Fig.13 Campbell diagram of crankshaft with different clearances

為了保證轉(zhuǎn)子正常運(yùn)行，需要使轉(zhuǎn)子工作轉(zhuǎn)速避開(kāi)臨界轉(zhuǎn)速，其安全裕度一般為10%～20%以上。結(jié)合圖可知，由于第1個(gè)臨界轉(zhuǎn)速很小，轉(zhuǎn)子可以很快跨躍，不會(huì)產(chǎn)生共振，對(duì)于第2個(gè)臨界轉(zhuǎn)速而言，間隙小于0.05 mm時(shí)，第2個(gè)臨界轉(zhuǎn)速為1 450 rad/s（約為13 800 r/min），遠(yuǎn)大于渦旋壓縮機(jī)正常運(yùn)行轉(zhuǎn)速范圍2 500～4 000 r/min，因此此種情形下也不會(huì)引起共振現(xiàn)象。然而當(dāng)渦旋壓縮機(jī)向高轉(zhuǎn)速發(fā)展時(shí)，就必須考慮間隙對(duì)該臨界轉(zhuǎn)速的影響，否則轉(zhuǎn)速越大，間隙越大時(shí)，越接近該臨界轉(zhuǎn)速，轉(zhuǎn)子發(fā)生共振的可能性就越大，圖14所示為曲軸的臨界轉(zhuǎn)速隨間隙大小的變化情況。

圖14 不同間隙下曲軸臨界轉(zhuǎn)速Fig.14 Crankshaft critical speed with different clearances

由圖可知，當(dāng)間隙小于0.05 mm時(shí)，隨著間隙的增大，曲軸系統(tǒng)的臨界轉(zhuǎn)速隨之緩慢降低，是因?yàn)殚g隙的存在破壞了曲柄銷與驅(qū)動(dòng)軸承的接觸剛度并使其發(fā)生變化，接觸剛度隨間隙的增大而減小，從而使曲軸系統(tǒng)的臨界轉(zhuǎn)速隨之變小。當(dāng)間隙大于0.05 mm時(shí)，曲軸臨界轉(zhuǎn)速隨著間隙增大降低越快，極大地增加了渦旋壓縮機(jī)發(fā)生共振的可能性，需要通過(guò)嚴(yán)格控制間隙大小和施加預(yù)緊力等措施以增加接觸剛度，以避免轉(zhuǎn)子共振現(xiàn)象。

5 結(jié)論

（1）同一間隙下，曲軸的位移波動(dòng)最小，驅(qū)動(dòng)軸承角速度的波動(dòng)較大，角加速度的波動(dòng)最大。

（2）不同間隙時(shí)，曲軸及驅(qū)動(dòng)軸承的位移、角速度、角加速度以及間隙接觸力均呈現(xiàn)不同程度的波動(dòng)，隨著間隙增大，波動(dòng)幅值越大，但波動(dòng)頻率降低。

（3）曲軸與驅(qū)動(dòng)軸承以及帶輪連接處均發(fā)生了應(yīng)力集中現(xiàn)象，而間隙的存在使曲軸的應(yīng)力集中現(xiàn)象加劇，隨著間隙的增大，曲軸所受的應(yīng)力也逐漸增大。

（4）間隙的存在改變了曲柄銷與驅(qū)動(dòng)軸承的接觸剛度，間隙增大，剛度減小，導(dǎo)致曲軸系統(tǒng)的臨界轉(zhuǎn)速降低，增加了轉(zhuǎn)子共振的危險(xiǎn)。