滑塊式球形轉(zhuǎn)子泵的振動(dòng)抑制研究

李世林，白國(guó)長(zhǎng)

(鄭州大學(xué)機(jī)械與動(dòng)力工程學(xué)院，河南鄭州 450000)

0 前言

球形泵屬于回轉(zhuǎn)式容積泵，其轉(zhuǎn)子的形狀為球形，可實(shí)現(xiàn)超微型化、高可靠性的容積式動(dòng)力機(jī)械，在壓縮機(jī)和液壓泵領(lǐng)域得到了相當(dāng)廣泛的應(yīng)用。滑塊式球形動(dòng)力機(jī)械理論體系及知識(shí)產(chǎn)權(quán)體系，是由王陸一、李正平[1]原創(chuàng)發(fā)明的、具有顛覆性意義的動(dòng)力機(jī)械原理的開(kāi)拓性發(fā)明成果；利用此原理研制出來(lái)的球形微型液壓泵結(jié)構(gòu)簡(jiǎn)單，運(yùn)行穩(wěn)定、變工況性能好，能實(shí)現(xiàn)超小體積、高壓力，可用于機(jī)器人、工程機(jī)械、機(jī)床工業(yè)、航空航天、VR裝備等產(chǎn)業(yè)。

但是由于滑塊式球形泵的旋轉(zhuǎn)盤(pán)做往復(fù)式運(yùn)動(dòng)，在工作時(shí)會(huì)產(chǎn)生陀螺力矩，是球形泵機(jī)械振動(dòng)的主要來(lái)源。而且，滑塊式球形泵的容積效率高，一周雙排量的特性，產(chǎn)生較大的流量脈動(dòng)，必然會(huì)造成振動(dòng)加劇，降低驅(qū)動(dòng)的平順性，影響泵的性能和使用壽命，對(duì)滑塊式球形泵的振動(dòng)進(jìn)行研究尤為必要。目前國(guó)內(nèi)外對(duì)球形泵的研究相對(duì)較少，西安交通大學(xué)的GUAN等[2]對(duì)一種偏置轉(zhuǎn)子式球形泵的運(yùn)動(dòng)學(xué)和動(dòng)力學(xué)進(jìn)行了深入研究，結(jié)果顯示由于沒(méi)有往復(fù)運(yùn)動(dòng)的部件，球形泵的噪聲只有傳統(tǒng)電磁泵的1/4。南洋理工大學(xué)的CHEN等[3]研究了柱塞泵的動(dòng)態(tài)振動(dòng)特性，發(fā)現(xiàn)機(jī)械振動(dòng)幅值隨著轉(zhuǎn)速和扭矩的增大而增大。康考迪亞大學(xué)的BAHR等[4]建立了一個(gè)數(shù)學(xué)模型來(lái)研究恒功率調(diào)節(jié)斜盤(pán)軸向柱塞泵的振動(dòng)特性，發(fā)現(xiàn)在橫向力矩的激勵(lì)作用下，比例閥芯以與力矩相同的頻率周期性振動(dòng)。本文作者研究的滑塊式球形泵存在往復(fù)運(yùn)動(dòng)件，采用雙泵抵消力矩的方法減小機(jī)械振動(dòng)，并研究雙泵對(duì)流量脈動(dòng)的影響。

1 球形轉(zhuǎn)子泵的振動(dòng)機(jī)制

球形轉(zhuǎn)子泵的核心是轉(zhuǎn)子機(jī)構(gòu)，主軸上沿徑向有一條滑槽，旋轉(zhuǎn)盤(pán)下端的滑塊插入滑槽形成平面副，旋轉(zhuǎn)盤(pán)通過(guò)銷(xiāo)軸與活塞連接，活塞頂部的軸與缸蓋孔連接，活塞的軸相對(duì)于垂直方向有一個(gè)角度α，稱(chēng)之為活塞角。主軸驅(qū)動(dòng)旋轉(zhuǎn)盤(pán)，旋轉(zhuǎn)盤(pán)帶動(dòng)活塞繞其自身軸線旋轉(zhuǎn)。而旋轉(zhuǎn)盤(pán)的運(yùn)動(dòng)是兩種運(yùn)動(dòng)組合的結(jié)果：一個(gè)是圍繞著主軸軸線的旋轉(zhuǎn)，另一個(gè)是相對(duì)滑槽面的往復(fù)擺動(dòng)，造成機(jī)械振動(dòng)。

缸蓋、缸體和球形轉(zhuǎn)子形成密封空間，當(dāng)主軸轉(zhuǎn)動(dòng)時(shí)，容積腔體積變化，起到吸排液體的作用[2]。由于球形轉(zhuǎn)子運(yùn)動(dòng)的非線性，流量并不是均勻的，產(chǎn)生流量脈動(dòng)。

圖1 球形轉(zhuǎn)子泵

2 球形泵動(dòng)力學(xué)分析

2.1 動(dòng)力學(xué)建模

根據(jù)滑塊式球形泵的結(jié)構(gòu)和工作原理，其結(jié)構(gòu)簡(jiǎn)化模型見(jiàn)圖2。活塞繞軸線2旋轉(zhuǎn)；旋轉(zhuǎn)盤(pán)的運(yùn)動(dòng)可分解為兩種：(1)繞旋轉(zhuǎn)軸1做角速度為ω的勻速轉(zhuǎn)動(dòng)，其為絕對(duì)運(yùn)動(dòng)；(2)滑塊在滑槽平面以點(diǎn)O為圓心做往復(fù)鐘擺運(yùn)動(dòng)，軌跡為A″-A1-A′，其為相對(duì)運(yùn)動(dòng)。建立兩個(gè)原點(diǎn)都位于球心點(diǎn)O的空間直角坐標(biāo)系xyz和x′y′z′，x′y′z′坐標(biāo)系是xyz坐標(biāo)系繞x軸旋轉(zhuǎn)角度α得到[2，5]。

A是旋轉(zhuǎn)盤(pán)滑塊中心，B是銷(xiāo)軸軸線的外端點(diǎn)，銷(xiāo)軸半徑r，且OA=OB=R，R是球形轉(zhuǎn)子的半徑。假設(shè)滑塊中心在z軸上時(shí)為球形泵的參考起始狀態(tài)C-A1-B1-D，驅(qū)動(dòng)主軸從C-O′-A1轉(zhuǎn)到C-O″-A，在平面xoy上的角位移記為θ，銷(xiāo)軸從O-B1旋轉(zhuǎn)到O-B，在平面x′oy′上角位移記為φ。OB與z軸夾角為σ1，OA與z軸夾角為旋轉(zhuǎn)盤(pán)擺角σ2。

圖2 球形轉(zhuǎn)子泵簡(jiǎn)化模型

在xyz坐標(biāo)系中，點(diǎn)A的位置可以表示為

(1)

點(diǎn)B在x′y′z′坐標(biāo)系的位置可以表示為

(2)

將式(2)通過(guò)坐標(biāo)變換方法用xyz坐標(biāo)系表示為

(3)

同時(shí)zB可以用σ1表示為

zB=-Rcosσ1

(4)

聯(lián)合求解式(3)和式(4)可以得到：

σ1=arccos(sinφsinα)

(5)

又因?yàn)镺B與OA呈90°，所以σ1+σ2=90°，解得擺角σ2為

(6)

當(dāng)旋轉(zhuǎn)主軸以等角速度ω旋轉(zhuǎn)時(shí)，活塞也將以角速度ω1運(yùn)轉(zhuǎn)，且ω1的大小會(huì)隨著旋轉(zhuǎn)主軸的轉(zhuǎn)角θ的改變而改變[6]，則有：

(7)

將上式(7)對(duì)時(shí)間t進(jìn)行求導(dǎo)，且dθ/dt=ω，可以得到活塞的角速度ω1：

(8)

對(duì)式(8)的角速度ω1再次求導(dǎo)就是活塞的角加速度ε1：

(9)

將式(7)代入式(6)中，并對(duì)σ2關(guān)于時(shí)間t求導(dǎo)可以獲取滑塊擺角的角速度ω2和角加速度ε2：

ω2=

(10)

(11)

2.2 轉(zhuǎn)子力矩分析

2.2.1 活塞的慣性力偶矩

周期性變化的角速度ε1的存在，使得活塞需要承受相應(yīng)的慣性力偶矩M2的作用，力矩方向沿活塞軸2，則根據(jù)活塞的運(yùn)動(dòng)特性，活塞的慣性力偶矩M2[7]可以表示為

(12)

式中：J2是活塞對(duì)轉(zhuǎn)軸2的轉(zhuǎn)動(dòng)慣量。

2.2.2 旋轉(zhuǎn)盤(pán)陀螺效應(yīng)分析

當(dāng)高速轉(zhuǎn)動(dòng)部件的轉(zhuǎn)動(dòng)軸在空間中改變方位時(shí)，轉(zhuǎn)動(dòng)部件必對(duì)約束部件作用一個(gè)附加力偶，影響系統(tǒng)的穩(wěn)定性，這種現(xiàn)象稱(chēng)為陀螺效應(yīng)，附加力偶的矩稱(chēng)為陀螺力矩[8-9]。由上一節(jié)的分析可知，旋轉(zhuǎn)盤(pán)的運(yùn)動(dòng)存在陀螺力矩，建立研究其陀螺力矩的動(dòng)力學(xué)模型，如圖3所示。

圖3 旋轉(zhuǎn)盤(pán)陀螺效應(yīng)動(dòng)力學(xué)模型

旋轉(zhuǎn)盤(pán)滑塊的點(diǎn)A繞著其旋轉(zhuǎn)軸3，以角速度ω2做自轉(zhuǎn)運(yùn)動(dòng)，點(diǎn)A的軌跡為圖3圓弧4上的a0-a1段；旋轉(zhuǎn)軸3繞旋轉(zhuǎn)軸1以角速度ω做勻轉(zhuǎn)速運(yùn)動(dòng)，此為旋轉(zhuǎn)軸3的公轉(zhuǎn)，旋轉(zhuǎn)軸3在過(guò)原點(diǎn)O且與軸1垂直的平面5上。旋轉(zhuǎn)盤(pán)在運(yùn)動(dòng)過(guò)程中有圖3所示的陀螺力矩M1產(chǎn)生，與軸線1和軸線3垂直。

根據(jù)陀螺力矩的計(jì)算公式推導(dǎo)出該球形泵的陀螺力矩M1：

M1=J1×ω×ω2

(13)

式中：J1是旋轉(zhuǎn)盤(pán)對(duì)旋轉(zhuǎn)軸3的轉(zhuǎn)動(dòng)慣量。

研究式(13)可知，當(dāng)J1為定值時(shí)，旋轉(zhuǎn)盤(pán)的陀螺力矩M1由旋轉(zhuǎn)盤(pán)滑塊的擺角速度ω2和ω決定，而根據(jù)式(10)，ω2的大小跟ω有關(guān)，所以陀螺力矩M1的大小取決于主軸轉(zhuǎn)速ω的值。將ω和ω2代入，得到陀螺力矩：

M1=J1·

(14)

2.2.3 空間力矩的合成

利用解析法將陀螺力矩M1和慣性轉(zhuǎn)矩M2在xyz坐標(biāo)系上分解[7]：

(15)

(16)

根據(jù)合力矢投影定理：

(17)

則可以得到合力偶矩矢的大小：

(18)

a、b、c分別是合力矩矢M與x、y、z軸的正向夾角，則合力偶矩矢的方向：

(19)

力矩分布曲線如圖4所示，活塞慣性轉(zhuǎn)矩為正弦周期分布。在轉(zhuǎn)角θ為90°和270°時(shí)的陀螺力矩M1是最小值為0；在0°和180°時(shí)達(dá)到最大值，且M1(0°)=M1(180°)。合成力矩M基本上與陀螺力矩M1曲線重合，說(shuō)明球形泵整體所受力矩的大小主要取決于旋轉(zhuǎn)盤(pán)陀螺力矩M1的值。

圖4 力矩曲線

3 球形泵容積

根據(jù)主軸轉(zhuǎn)角、轉(zhuǎn)盤(pán)偏角、活塞轉(zhuǎn)角3個(gè)參數(shù)，可推導(dǎo)出工作腔(容積腔)張角，進(jìn)而求得工作腔容積的變化規(guī)律[10-12]。該模型工作腔結(jié)構(gòu)比較規(guī)范，工作腔容積與夾角成正比，可以直接通過(guò)夾角計(jì)算出體積，如式(20)所示：

V=Kβ

(20)

式中：V為工作腔容積，mm3；K為轉(zhuǎn)換系數(shù)，mm3/rad ；β為工作腔張角，rad。

K代表工作腔內(nèi)每單位角度內(nèi)的容積，表征了工作腔的容納能力，大小與球形泵的結(jié)構(gòu)有關(guān)，K可以確定為

K=π/270(R2-r2)3/2

(21)

由空間機(jī)構(gòu)運(yùn)動(dòng)、動(dòng)力分析得工作腔張角：

(22)

4 減振措施

在球形泵工作過(guò)程中，慣性力矩和流量脈動(dòng)是振動(dòng)產(chǎn)生的主要源頭。可以通過(guò)減小主軸轉(zhuǎn)速ω和活塞角α，縮小球形轉(zhuǎn)子半徑R以降低轉(zhuǎn)動(dòng)慣量，但會(huì)影響泵的容積率，增大加工難度。作者采用成對(duì)球形轉(zhuǎn)子泵對(duì)稱(chēng)安裝的措施，來(lái)抵消力矩，并研究對(duì)流量脈動(dòng)的影響。轉(zhuǎn)子陀螺力矩M1比活塞慣性轉(zhuǎn)矩M2大得多，優(yōu)先減小陀螺力矩。

如圖5所示，兩泵對(duì)稱(chēng)安裝，分別建立直角坐標(biāo)系x1y1z1和x2y2z2。兩陀螺力矩矢量方向相反，且大小相等，全部抵消。兩慣性轉(zhuǎn)矩大小相同，分別沿O1D1和O2D2方向。經(jīng)空間矢量分解，合成整體力矩。

(23)

圖5 雙泵安裝示意

4.1 ADAMS振動(dòng)仿真分析

球形轉(zhuǎn)子泵的力矩和運(yùn)動(dòng)都是周期性的，存在周期變化的激振力，導(dǎo)致系統(tǒng)發(fā)生受迫振動(dòng)。ADAMS的Vibration插件可以分析系統(tǒng)在頻域的強(qiáng)迫振動(dòng)響應(yīng)特性，通過(guò)采用虛擬激振器來(lái)替代物理實(shí)驗(yàn)裝置。同時(shí)，也可以利用ADAMS/View進(jìn)行時(shí)域的振動(dòng)分析[13]。設(shè)置激振力是正弦掃頻信號(hào)(文中的轉(zhuǎn)速n=6 000 r/min，頻率f=n/60=100 Hz，周期T=1/f=0.01 s)，然后運(yùn)行分析[13-14]。經(jīng)過(guò)后處理，系統(tǒng)主軸質(zhì)心點(diǎn)處的位移、加速度的時(shí)域響應(yīng)和頻域響應(yīng)如圖6所示。

圖6為振動(dòng)的幅頻和相頻曲線，3個(gè)波動(dòng)區(qū)域即3階固有頻率[12]，工作頻率100 Hz位于2階頻率和3階頻率之間，不會(huì)因?yàn)楣舱駧?lái)重大損害，且振動(dòng)幅值響應(yīng)很小，保證了系統(tǒng)運(yùn)行的穩(wěn)定性。

圖6 頻域響應(yīng)

圖7(a)、(b)分別為單泵和雙泵振動(dòng)位移的時(shí)域響應(yīng)，周期為0.01 s的正弦曲線。單泵時(shí)的振動(dòng)位移幅值為0.91 mm，雙泵時(shí)的振動(dòng)位移幅值為0.000 7 mm，機(jī)械振動(dòng)大幅度降低。

圖7 振動(dòng)位移響應(yīng)曲線

4.2 流量曲線分析

對(duì)容積V用插值法取得工作室容積的導(dǎo)數(shù)，即球形泵的流量變化，可以用來(lái)分析流量波動(dòng)的大小[12]。單泵和雙泵合成時(shí)的流量曲線如圖8所示。

流量脈動(dòng)的大小可以用波動(dòng)率表示，波動(dòng)率δ=(Q最大流量-Q平均流量)/Q平均流量×100%。由圖8可以得到單泵的流量波動(dòng)率約45%，相比之下，雙泵組合的波動(dòng)率只有約 18%，相比降低了60%。

圖8 流量波動(dòng)曲線

通過(guò)仿真分析，對(duì)比單泵和雙泵的機(jī)械振動(dòng)和流量脈動(dòng)，證明了雙泵對(duì)機(jī)械振動(dòng)和流量脈動(dòng)的抑制成效顯著。實(shí)際工程項(xiàng)目中的多泵采用這種方式，可以很大程度上減小系統(tǒng)振動(dòng)，提高工作的穩(wěn)定性。

5 結(jié)論

振動(dòng)會(huì)對(duì)球形泵造成損傷，并且?guī)?lái)噪聲，因此研究泵的機(jī)械振動(dòng)和流量脈動(dòng)是必須的。對(duì)球形泵進(jìn)行動(dòng)力學(xué)建模，求解旋轉(zhuǎn)盤(pán)和活塞的角速度和角加速度，得到活塞的慣性力矩和旋轉(zhuǎn)盤(pán)的陀螺力矩。將球形轉(zhuǎn)子的力矩進(jìn)行合成分析可知，球形泵所受慣性力矩的大小主要取決于旋轉(zhuǎn)盤(pán)的陀螺力矩。本文作者在不改變結(jié)構(gòu)參數(shù)的前提下，采用多泵減振的方案，抵消力矩和流量脈動(dòng)。通過(guò)ADAMS對(duì)球形轉(zhuǎn)子泵進(jìn)行振動(dòng)仿真，頻域曲線證明了球形泵在激振條件下受迫振動(dòng)的穩(wěn)定性。建立球形泵的容積腔變化公式，插值法取得容積變化曲線的導(dǎo)數(shù)即球形泵流量，并使用MATLAB生成單泵和雙泵的流量曲線。對(duì)比單泵和雙泵的機(jī)械振動(dòng)位移和流量波動(dòng)率，雙泵狀態(tài)的振動(dòng)位移幅值和流量波動(dòng)率皆大幅度減小。綜上結(jié)果分析，采用多泵安裝既抵消了大部分力矩，又減小了流量脈動(dòng)，增加球形泵的工作壽命和平順性。