陳 就 夏林海 沙 豐 吳長波 陳劍平
中國民用航空局第二研究所航材技術研究中心,成都,610041
CP型單螺桿壓縮機具有體積小、容量大、單級壓縮比大和效率高等其他類型壓縮機無可比擬的優點,能夠較好地滿足軍事、醫藥等領域對壓縮機高性能和環保的要求[1]。國內外關于單螺桿壓縮機的研究主要集中在壓縮機幾何結構優化、壓縮性能、流場分析等結構和特性方面[2-5]。單螺桿壓縮機中的核心零件單螺桿齒槽結構為球面包絡,相比其他類型的齒類零件,其齒槽結構更加復雜,加工難度也更大。目前國內常見的CP型單螺桿(以下簡稱“單螺桿”)的加工方法主要有仿形車削法和圓柱二次包絡銑削法[1]。仿形車削主要通過成形車刀模擬星輪齒形進行螺旋切削,車刀刃越靠近前端其切削周期越長,局部磨損也越嚴重,另外,車刀刀臺轉速慢,單次切削空行程也很大,導致仿形車削存在加工精度差和效率低等問題。圓柱二次包絡銑削法結合了先進的伺服控制技術,采用成形銑刀對左右齒面分別銑削,該方法相比仿形車削法具有加工效率高、刀具磨損均勻等優點,但仍存在左右齒面分別銑削時出現接刀痕和螺桿齒槽的分度誤差等缺陷。
單次包絡成形加工方法是具有螺旋齒面特征零件高精度和高效率的齒形加工方法之一,采用回轉刀具(如銑刀、砂輪等)加工,既可以實現成形銑削,也可以實現成形磨削,如滾刀、蝸桿、齒輪、螺桿等零件的單次包絡成形,對此國內學者有深入的研究[6-10],研究結果表明這些零件都能通過單次包絡加工而成。而對于單螺桿齒槽可否采用單次包絡成形加工,還有待研究。
本文根據單螺桿的螺旋齒面結構,結合空間嚙合原理,建立了單次包絡成形加工數學模型,分析了理論加工的誤差,并根據計算得出的誤差值,探究了單次包絡成形應用于單螺桿加工的可行性。
為了便于計算,建立固定坐標系σ(OXYZ)、工件(單螺桿)坐標系σ1(O1X1Y1Z1)及刀具坐標系σ2(O2X2Y2Z2),如圖1所示。在起始狀態,刀具坐標系σ2(O2X2Y2Z2)與固定坐標系σ(OXYZ)重合,工件坐標系σ1(O1X1Y1Z1)沿固定坐標系σ(OXYZ)偏置位移(0,A,0)T。

圖1 計算坐標系Fig.1 Calculation coordinates system
刀具在單螺桿齒槽的成形運動有兩個:單螺桿繞X1軸的轉動,角速度為ω1;刀具隨刀臺繞Z2軸的轉動,角速度為ω2。為便于描述,以右旋齒形為例進行說明,令ω1=(1,0,0)T,根據單螺桿與星輪的聯動關系,有ω2=(0,0,k/z)T,其中,k為單螺桿的頭數,z為與單螺桿相配星輪的齒數。成形刀具繞自身軸線Y2做回轉運動與齒槽成形無關,令成形刀具的自轉角速度為0。
令成形刀具回轉軸線Y2與固定坐標系軸線Y重合,且單螺桿在X1O1Y1軸截面上的某個齒形關于Y2軸對稱的時刻稱為初始狀態,該截面齒形為基準齒形,如圖2所示,其中B為單螺桿齒槽寬,D為單螺桿外徑,h為齒槽深。對基準齒形作離散處理,所有離散點將分別在成形加工的某一時刻進入嚙合,通過空間嚙合原理確定刀具與螺桿接觸時的空間嚙合點。由嚙合原理可知,這些空間嚙合點既是單螺桿齒槽回轉面上的點,也是成形刀具回轉面上的點,將這些點分別繞Y2軸回轉一定角度返回到X2OY2軸截面上,即可得到成形刀具的截面刃形。

圖2 加工初始狀態示意圖Fig2 The initial state diagram of machining
將基準齒形離散成N個離散點,記任一離散點P在工件坐標系σ1中的徑矢為RW,截面切矢為τW,由于P點及其截面切矢τW均在軸截面X1O1Y1上,因此對基準齒形方程上的P點求導,有
(1)
式中,f(x,y)為X1O1Y1軸截面上的基準齒形方程。
設P點進入嚙合時,螺桿轉動角度為β,刀臺轉動角度為θ,由前述聯動關系可知:
(2)
通過坐標變換可知,在固定坐標系中的徑矢和切矢分別為
(3)
(4)
j=(0,1,0)T

在刀具坐標系中的徑矢和切矢分別為
(5)
(6)

在固定坐標系σ(OXYZ)中,P點處的成形運動有兩個,分別為隨螺桿繞X1軸轉動的線速度v1和隨刀臺繞Z2軸轉動的線速度v2:
v1=ω1×(R-Aj)
(7)
v2=ωT×R
(8)
在坐標系σ(OXYZ)中,當P點成為嚙合點時,該點處的相對運動速度為
v12=v1-v2
(9)
單螺桿齒槽面由成形刀具單次包絡而成,且齒槽基準齒形是包絡面上的曲線,由解析幾何可知,在固定坐標系σ中,P點處的法矢n垂直于P點的相對運動速度v12和截面切矢τ:
n=v12×τ
(10)
轉至刀具坐標系σ2(O2X2Y2Z2)中即為
(11)
P點為嚙合點時,由嚙合原理可知[11],在刀具坐標系σ2中,回轉刀具上的徑矢RT、法矢nT及回轉刀具的回轉軸線方向向量j三者的混合積為0,即
(RT,nT,j)=0
(12)
解式(12)可以得到P點成為嚙合點時螺桿的轉動角度β,代入式(3)中可以得到P點在固定坐標系σ中的嚙合點徑矢。
由于P點是基準齒形中離散點中的任意一點,故對所有離散點按以上方法進行計算就可以得到基準齒形中每一個離散點在進入嚙合時螺桿的轉動角度βi(i=1~N,N為計算的離散點個數)。將每一個βi代入徑矢R計算式中,可以得到一系列離散點在固定坐標系σ中的嚙合點徑矢矩陣,進而得到單螺桿基準齒形在單次包絡成形時的空間嚙合線。
將上述每一個計算得到的βi值依次代入式(5)中,就可以得到空間嚙合線在刀具坐標系σ2中的徑矢矩陣。由于空間嚙合線也在回轉刀具回轉面上,故依次將空間嚙合線上所有離散點繞Y2回轉至回轉刀具軸截面(X2OY2)上,設第i點的回轉角度為αi,將該點回轉到刀具軸截面(XTOYT)上的徑矢為
(13)

同理可以得到該點的截面法矢:
(14)
回轉后RS(i)是截面(XTOYT)上的點,其Z向分量為0,令
Z(RS(i))=0
(15)
求解式(15)可得出第i點的回轉角度αi,代入RS(i)和nS(i)計算式中,即可求得第i點的截面徑矢和截面法矢,也即對應基準刃形上的任意點P都可以在回轉刀具截面刃形上得到一點F與之對應。對所有離散點進行求解即可得到成形刀具截面刃形上所有離散點的截面徑矢和截面法矢,對徑矢點進行曲線擬合,即可得到回轉刀具的截面刃形。
成形誤差指單螺桿在過X1軸的任意軸截面上齒形離散點的理論徑矢和單次包絡成形徑矢之間的齒形偏差。根據單螺桿齒槽面的成形原理,理論徑矢可以通過基準齒形上的離散點沿特定的螺旋線掃掠得到。齒槽截面成形包絡徑矢方程則需要根據上述計算得到的成形刀具截面徑矢和法矢,結合理論徑矢所在的齒形軸截面,運用嚙合原理反推得到。
進行誤差計算時,設過螺桿軸線的任一齒槽軸截面與基準齒形所在軸截面之間的夾角為γ,記該齒槽軸截面為γ軸截面。γ軸截面齒形上Q點由基準齒形上任意點P掃掠得到,將該點回轉至螺桿所在的齒槽截面上,可得到螺桿在工件坐標系σ1中的螺旋面徑矢方程:
(16)
由于P點為基準齒形上的任意一點,因此將基準齒形上的所有離散點代入式(16)中,就可以得到單螺桿γ軸截面上齒形的所有離散點的理論徑矢方程。
在包絡成形中,單螺桿γ軸截面上的齒形離散點由前述計算得到的成形刀具截面刃形離散點通過嚙合計算得到。由前文可知,成形刀具截面刃形離散點由基準刃形離散點通過嚙合計算得到,也即基準刃形的離散點之間存在一一映射關系。也就是說基準刃形上的離散點P通過其對應的空間嚙合點回轉到刀具截面刃形上,對應刀具截面刃形上的F點,同理,通過嚙合計算,也可以得到與F點相對應的單螺桿γ軸截面上的E點,如圖3所示。

圖3 E點嚙合計算說明圖Fig.3 E point meshing calculation diagram
當E點進入嚙合時,E點與F點重合,記此時單螺桿轉角為β1,對應刀臺轉角為θ1,F點繞Y2軸旋轉至E點位置,轉動角度為α1,根據聯動關系有
(17)
此時,E點在刀具坐標系中的徑矢和法矢分別為
(18)
(19)

在固定坐標系中的徑矢和法矢分別為
(20)
(21)

隨螺桿及刀臺的轉動速度分別為
v3=ω1×(RE-Aj)
(22)
v4=ω2×RE
(23)
嚙合時其相對運動速度為
v34=v3-v4
(24)
依據嚙合原理,有
nv34=0
(25)
E點在工件坐標系中的徑矢為
(26)
又由于E點在γ軸截面齒形上,根據幾何關系有

(27)
聯合式(25)、式(27)可求出E點成為嚙合點時相應的轉角α1和φ1,代入RE表達式中即可得到對應γ軸截面上E點的徑矢。將成形刀具截面上所有點代入運算,即可得到單螺桿上對應γ軸截面齒形的包絡線。同理對每一個螺桿轉角γ進行求解即可得到螺桿齒形的包絡成形曲面。
由前述計算分別得到了單螺桿基準齒形上的P點對應γ軸截面齒形上的理論點Q和包絡成形點E的徑矢。Q點和E點的距離就是對應P點在γ軸截面上的成形誤差:
(28)
通過對基準刃形上所有離散點進行計算就可以得到任意夾角γ下的齒形成形誤差分布。
應用上述單次包絡成形計算方法,對一種常見規格的單螺桿齒槽分別計算其包絡成形加工刀具的截面刃形和單螺桿單次包絡成形誤差。在計算中,本文只對包絡成形加工的可行性作定性分析,為簡化計算過程,不考慮成形包絡點落在單螺桿幾何體以外的情況。
某型單螺桿部分設計參數如下:螺桿外徑D=300 mm,齒深h=40 mm,齒寬B=30 mm,齒根圓弧半徑r=3 mm,頭數k=6;螺桿星輪中心距A=250 mm,相配星輪齒數z=11。
取基準齒形離散點個數N=49,利用MATLAB軟件對齒形進行離散計算,離散后所有離散點的基準齒形見圖4。

圖4 離散點的基準齒形Fig.4 Base tooth profile of discrete dots
基準齒形所有離散點進入嚙合時,在固定坐標系σ中的空間嚙合位置如圖5所示。

圖5 離散點的空間嚙合位置Fig.5 Space meshing position of discrete dots
由圖5可以看出,齒根線兩端點在嚙合時,嚙合點位置產生了突變,對應到基準齒形上為圓弧段與齒根線的交點位置。此時,需檢查求解式(12),判斷變量β1在其有效的成形轉角范圍內是否為唯一解。由包絡成形加工方法可知,成形刀具的直徑應不大于所在齒槽位置的寬度,否則在實際應用中會產生干涉。根據解析幾何,以左齒面的交點為例進行說明,所有離散點的轉角β1的解空間應滿足:
(29)
根據式(12),令
f(β1)=(RT,nT,j)
(30)
可得函數f(β1)和變量β1之間的關系曲線,如圖6所示。f(β1)在解空間內單調遞增,僅有一個解,說明齒根線兩端點在嚙合時,不存在除計算出的嚙合點位置之外的其他可能的空間嚙合點位置,也即基準齒形在包絡成形時,該點處確實存在嚙合位置突變。

圖6 函數f(β1)在解空間內的曲線圖Fig.6 Curve graph of function f(β1)in solution space
根據離散點的空間嚙合位置,進一步求得基準齒形對應的包絡成形的成形刀具截面刃形,如圖7所示。由圖7可以看出,由于齒根線兩端點的空間嚙合點位置發生了突變,導致成形刀具截面刃形在對應點處產生了突出的尖角,該刃形是無法應用于實際的。如果要應用單次包絡成形進行單螺桿齒槽的成形加工,則需要優化設計參數。不斷增大齒根圓弧半徑r值,然后分別進行計算,結果表明,增大齒根圓弧半徑能有效緩解刃形的尖角突出的問題,此處給出了具有代表性的r=7,10,15 mm(此時基準齒形無直線齒根線)對應的成形刀具截面刃形,分別如圖8~圖10所示。

圖7 成形刀具截面刃形(基準齒形)Fig.7 Curve graph of forming tool profile(base tooth profile)

圖8 成形刀具截面刃形(r=7 mm)Fig.8 Curve graph of forming tool profile(r=7 mm)

圖9 成形刀具截面刃形(r=10 mm)Fig.9 Curve graph of forming tool profile(r=10 mm)
由圖8和圖9可以看出,如果基準齒形上的齒根線為直線,那么該齒根線對應的成形刀具截面刃形內凹,在該段刃形的兩端就會產生曲率突變,在實際應用中,該兩點位置就會產生急劇磨損,影響加工精度,不是包絡成形理想的基準齒形。而由圖10可以看出,將齒根線直接改成了圓弧后,離散點的嚙合空間位置(圖11)和成形刀具截面刃形均平滑過渡,是包絡成形理想的基準齒形。為此,下面將對r=15 mm(齒根為圓弧)的截面齒形進行包絡成形誤差分析。

圖10 成形刀具截面刃形(r=15 mm)Fig.10 Curve graph of forming tool profile(r=15 mm)

圖11 離散點的空間嚙合位置(r=15 mm)Fig.11 Space meshing position of discrete dots(r=15 mm)
在成形誤差分析中,基于單螺桿齒槽關于Y1軸的對稱性,計算Y1軸齒槽長的一側(圖2中Y1軸的左側)的誤差即可。根據聯動和幾何關系可得到該側齒槽轉角γ的取值范圍:
(31)
將設計參數代入計算可得:0≤γ≤96.47°。從0°開始每隔5°取一個誤差計算截面,最后一個截面取γ=96.47°,得到的理論和單次包絡成形曲面對比圖見圖12,局部放大圖見圖13。

圖12 理論和包絡成形曲面對比圖Fig.12 Curve graph of comparison between theory and envelope forming surface

圖13 理論和包絡成形齒形對比局部放大圖Fig.13 Partial enlarged detail curve graph of comparison between theory and envelope forming surface
對每一個γ軸截面上理論和包絡成形齒形之間的誤差值進行計算,計算結果和齒形誤差分布如圖14所示,最大誤差為1.061 mm。

圖14 不同γ軸截面上齒形誤差分布情況Fig.14 The tooth profile error distribution on different axis sections with correspond γ angles
從以上計算結果可以看出:
(1)當單螺桿基準齒形的齒根線為直線時,其包絡成形刀具的截面刃形在對應齒根線左右兩端點位置會出現突出尖角,該兩點處左右兩側的曲率變化量和齒根圓角半徑值有直接關系,當齒根圓角增大時,曲率變化量會減小,刃形曲線會更加平滑,但突出的尖角始終存在。
(2)當單螺桿基準齒形的齒根線為圓弧時,其包絡成形刀具的截面刃形平滑,無尖角突出,但單次包絡成形導致的齒形誤差值較大,最大位置的誤差超過了1 mm。
(3)優化后的基準齒形成形誤差分布特點為:齒形誤差值在齒形的中心點處為零,向兩側先增大后減小,且兩側的誤差值不對稱,對于右旋螺桿,其左齒面誤差值較大;隨著轉動角度γ的不斷增大,γ軸截面上的齒形點進入嚙合時,刀臺的轉動角度也不斷增大,相對于基準齒形上的點,其相對速度的變化量也逐漸變大,導致誤差也逐漸增大。
(1)對CP單螺桿齒槽的單次包絡成形方法進行了理論研究,依據嚙合原理建立了采用單次包絡成形加工刀具截面刃形、理論齒槽曲面徑矢及包絡成形齒槽曲面徑矢的計算方法與公式。
(2)利用該計算方法對一種常規的單螺桿齒槽進行了計算,結果表明該方法加工單螺桿齒槽曲面成形誤差較大,且離齒槽中心點越遠誤差值越大,無法通過單次包絡成形加工方法將CP型單螺桿齒槽一次加工成形。