基于NURBS曲線的螺桿轉(zhuǎn)子型線反向設(shè)計(jì)研究

黃海楠HUANG Hai-nan 何雪明 - 張 榮 姜振鋼 -

(1. 江南大學(xué)江蘇省食品先進(jìn)制造裝備技術(shù)重點(diǎn)實(shí)驗(yàn)室，江蘇 無錫 214122；2. 江南大學(xué)理學(xué)院，江蘇 無錫 214122)

雙螺桿壓縮機(jī)因?yàn)樾阅軆?yōu)良、運(yùn)行平穩(wěn)、噪聲低、體積小、維護(hù)方便等，被廣泛應(yīng)用于制冷[1-2]、食品[3-4]、石油、化工等行業(yè)，其核心零件為內(nèi)部的一對(duì)陰陽轉(zhuǎn)子。轉(zhuǎn)子型線作為螺桿轉(zhuǎn)子形狀的重要特征，很大程度上決定了壓縮機(jī)的工作性能及轉(zhuǎn)子的加工制造性能。

目前轉(zhuǎn)子型線的設(shè)計(jì)方法使用比較廣泛的是正向設(shè)計(jì)：首先設(shè)計(jì)出陰轉(zhuǎn)子(或者陽轉(zhuǎn)子)的型線，再通過解析包絡(luò)法推導(dǎo)出另一個(gè)轉(zhuǎn)子的型線。Stosic等[5]根據(jù)齒條的嚙合原理提出了由齒廓線到轉(zhuǎn)子型線的設(shè)計(jì)方法，并開發(fā)了N型線；Su等[6]設(shè)立轉(zhuǎn)子接觸線長度、泄露三角形面積等性能參數(shù)作為目標(biāo)，實(shí)現(xiàn)對(duì)轉(zhuǎn)子型線的優(yōu)化；Wu等[7]提出了基于標(biāo)準(zhǔn)齒條生成系統(tǒng)的螺桿轉(zhuǎn)子輪廓設(shè)計(jì)方法，并通過對(duì)比論證闡釋了該方法的優(yōu)越性；何雪明等[8]提出了以自由曲線作為組成曲線的轉(zhuǎn)子設(shè)計(jì)方法，豐富了轉(zhuǎn)子型線線型的種類。

正向設(shè)計(jì)雖然比較成熟，但要判斷所設(shè)計(jì)的轉(zhuǎn)子型線及壓縮機(jī)性能的優(yōu)劣，必須在得到陰陽轉(zhuǎn)子型線后，繼續(xù)推導(dǎo)得到嚙合線再做出判斷，設(shè)計(jì)過程十分繁瑣。針對(duì)正向設(shè)計(jì)的多個(gè)弊端，Zaytsev等[9-10]提出了反向設(shè)計(jì)轉(zhuǎn)子型線的“嚙合線法”，該方法是從嚙合線出發(fā)推導(dǎo)出轉(zhuǎn)子輪廓線，并通過實(shí)例驗(yàn)證了該方法的可行性；卜嘯華等[11]推導(dǎo)了當(dāng)嚙合線曲線類型為圓弧、內(nèi)擺線和外擺線時(shí)的螺桿轉(zhuǎn)子型線計(jì)算公式；何雪明等[12]分別采用直線、圓弧、三次B樣條曲線為嚙合線段進(jìn)行陰陽轉(zhuǎn)子的反向設(shè)計(jì)。由此，嚙合線推導(dǎo)轉(zhuǎn)子型線的反向設(shè)計(jì)方法逐漸發(fā)展起來。

在反向設(shè)計(jì)的線型選擇方面，現(xiàn)在已經(jīng)發(fā)展到比較高級(jí)的B樣條曲線，但B樣條曲線的根基是多項(xiàng)式參數(shù)曲線，并不能表達(dá)基本的二次曲線，如圓弧、橢圓弧等。而NURBS曲線擁有更高的自由度和局部修改性，可以精確地表達(dá)圓弧、橢圓弧等二次曲線。本研究擬將NURBS曲線運(yùn)用到反向設(shè)計(jì)中，旨在實(shí)現(xiàn)設(shè)計(jì)曲線線型的統(tǒng)一，并利用靈活的曲線修改手段進(jìn)行轉(zhuǎn)子型線優(yōu)化設(shè)計(jì)。

1 嚙合線法設(shè)計(jì)原理

已知雙螺桿轉(zhuǎn)子的嚙合線的方程式是x0、y0關(guān)于參數(shù)t的方程，通過嚙合線法求解轉(zhuǎn)子輪廓線，需要建立如圖1所示的坐標(biāo)系，包括陽轉(zhuǎn)子靜坐標(biāo)系O1X1Y1和陰轉(zhuǎn)子靜坐標(biāo)系O2X2Y2、隨轉(zhuǎn)子轉(zhuǎn)動(dòng)的陽轉(zhuǎn)子動(dòng)坐標(biāo)系O1x1y1和陰轉(zhuǎn)子動(dòng)坐標(biāo)系O2x2y2，以及一個(gè)固結(jié)于機(jī)殼上的嚙合線坐標(biāo)系O0x0y0。其中陽轉(zhuǎn)子逆時(shí)針轉(zhuǎn)動(dòng)，陰轉(zhuǎn)子順時(shí)針轉(zhuǎn)動(dòng)。圖1中Ra1是陽轉(zhuǎn)子齒頂圓半徑，Rf2是陰轉(zhuǎn)子齒根圓半徑，滿足Ra1+Rf2=H。以下下標(biāo)為1都是指代陽轉(zhuǎn)子坐標(biāo)系中的參數(shù)，下標(biāo)為2的指代陰轉(zhuǎn)子坐標(biāo)系中的參數(shù)。

圖1 轉(zhuǎn)子法中的坐標(biāo)關(guān)系Figure 1 Relationships of coordinates in the rotor method

根據(jù)轉(zhuǎn)子間的嚙合傳動(dòng)原理可得：

(1)

φ1+φ2=(1+i)φ1=kφ1，

(2)

R1+R2=H，

(3)

式中：

φ2、φ1——陰、陽轉(zhuǎn)子旋轉(zhuǎn)角度，(°)；

n2、n1——陰、陽轉(zhuǎn)子的轉(zhuǎn)速，r/min；

ω2、ω1——陰、陽轉(zhuǎn)子轉(zhuǎn)動(dòng)角速度，rad/s；

R2、R1——陰、陽轉(zhuǎn)子節(jié)圓半徑，mm；

z2、z1——陰、陽轉(zhuǎn)子齒數(shù)；

i——陰轉(zhuǎn)子和陽轉(zhuǎn)子的傳動(dòng)比；

H——陰、陽轉(zhuǎn)子中心距，mm。

根據(jù)上述的等式以及坐標(biāo)轉(zhuǎn)動(dòng)中的幾何關(guān)系，可以推導(dǎo)出下列幾項(xiàng)坐標(biāo)系之間的變換關(guān)系：

(1) 陽轉(zhuǎn)子動(dòng)坐標(biāo)系O1x1y1與嚙合線坐標(biāo)系O0x0y0關(guān)系：

(4)

(2) 陰轉(zhuǎn)子動(dòng)坐標(biāo)系O2x2y2與陽轉(zhuǎn)子動(dòng)坐標(biāo)系O1x1y1之間的變換關(guān)系：

(5)

根據(jù)包絡(luò)條件可以得到φ1關(guān)于參數(shù)t的等式[13]：

(6)

2 嚙合線法的NURBS曲線包絡(luò)條件

2.1 包絡(luò)條件的求解

一段k次NURBS曲線的函數(shù)定義為[14]：

(7)

式中：

k——曲線次數(shù)；

Pi——控制點(diǎn)；

wi——控制點(diǎn)Pi的權(quán)因子，決定了控制點(diǎn)偏離曲線的程度，所有wi>0；

Ni,k(u)——定義在非周期且非均勻節(jié)點(diǎn)矢量U={a,…,a,uk+1,…,um-p-1,b,…,b}上的k次B樣條基函數(shù)。

一段嚙合線的NURBS曲線段的參數(shù)方程為：

0≤u≤1。

(8)

將式(8)代入式(6)中可以得到：

(9)

(10)

φ0是嚙合的起始角度，對(duì)于嚙合線的第一段曲線，取φ0=0。

f(u)=

(11)

由于被積函數(shù)f(u)的原函數(shù)無法用初等函數(shù)表達(dá)，因此需要采用數(shù)值積分法，本設(shè)計(jì)方法采用的是龍貝格(Romberg)求積公式[15-16]，也稱為逐次分半加速法。龍貝格求積公式用的是一種外推算法，相比于其他計(jì)算方式，其計(jì)算量相對(duì)較小、精度較高，且易于編程計(jì)算。龍貝格求積公式為：

(12)

其中兩個(gè)初始值的計(jì)算公式為：

(13)

(14)

由此得到φ1值后，再將值代入式(4)、(5)就得到了陰轉(zhuǎn)子和陽轉(zhuǎn)子對(duì)應(yīng)的齒曲線段。

2.2 NURBS曲線的求導(dǎo)

(15)

(16)

而A(u)與K(u)的導(dǎo)數(shù)都可以轉(zhuǎn)換成關(guān)于基函數(shù)的導(dǎo)數(shù)，由式(17)～(20)求得：

(17)

(18)

(19)

(20)

式(17)～(20)中的基函數(shù)的一階導(dǎo)數(shù)和二階導(dǎo)數(shù)可以由式(21)求得。

(21)

(22)

可以看到式(21)～(22)是類似于求解基函數(shù)的一個(gè)遞推公式，可以從低到高一階一階地計(jì)算出基函數(shù)導(dǎo)數(shù)的值。由此便可計(jì)算出NURBS曲線的一階導(dǎo)數(shù)和二階導(dǎo)數(shù)。

3 NURBS曲線嚙合線的設(shè)計(jì)與優(yōu)化

本研究針對(duì)用NURBS曲線對(duì)復(fù)盛型線[17]的嚙合線進(jìn)行反向設(shè)計(jì)，并通過誤差分析，驗(yàn)證基于NURBS曲線的反向設(shè)計(jì)方法的可行性，并在此基礎(chǔ)上，利用NUBRS的局部可調(diào)性對(duì)嚙合線進(jìn)行了優(yōu)化設(shè)計(jì)。

3.1 NURBS曲線嚙合線設(shè)計(jì)流程

NURBS曲線設(shè)計(jì)嚙合線設(shè)計(jì)流程如圖2所示，首先必須確定陰陽轉(zhuǎn)子的中心距，然后利用NURBS直接設(shè)計(jì)出嚙合線，根據(jù)NURBS曲線的數(shù)據(jù)推導(dǎo)出陰、陽轉(zhuǎn)子型線數(shù)據(jù)，如果符合設(shè)計(jì)要求，則輸出型線數(shù)據(jù)，不符合則需要修改設(shè)計(jì)參數(shù)，重新計(jì)算。

圖2 NURBS曲線反向設(shè)計(jì)流程Figure 2 Reverse design process with NURBS curve

3.2 復(fù)盛型線的反向設(shè)計(jì)

復(fù)盛型線轉(zhuǎn)子型線由圓弧、橢圓弧及其對(duì)應(yīng)的包絡(luò)線設(shè)計(jì)而成，見圖3。圓弧以參數(shù)方程的形式表達(dá)為：

(23)

式中：

r——圓弧的半徑，mm；

c1、c2——圓心橫坐標(biāo)和縱坐標(biāo)，mm。

橢圓弧以參數(shù)方程的形式表達(dá)為：

(24)

式中：

a、b——橢圓的長短軸長，mm；

c1、c2——橢圓中心橫坐標(biāo)和縱坐標(biāo)，mm。

復(fù)盛型線的各段齒曲線的詳細(xì)設(shè)計(jì)參數(shù)見表1。其嚙合線由陰陽轉(zhuǎn)子齒曲線通過嚙合條件式計(jì)算所得，表2給出了嚙合線曲線段和轉(zhuǎn)子齒曲線段的對(duì)應(yīng)關(guān)系，生成的嚙合線見圖4。

根據(jù)復(fù)盛的嚙合線數(shù)據(jù)，利用高精度的NURBS曲線來代替原有的嚙合線[見圖5(a)～(d)]進(jìn)而得到了該嚙合線每一段的控制點(diǎn)與權(quán)因子；然后將每一段嚙合線的端點(diǎn)作為NURBS曲線的節(jié)點(diǎn)，利用累積弦長法計(jì)算出整段型線的節(jié)點(diǎn)矢量，使其成為一段式的嚙合線[見圖5(e)]，最終設(shè)計(jì)的這段嚙合線只由控制點(diǎn)、權(quán)因子和節(jié)點(diǎn)矢量決定。進(jìn)一步利用上文推出的NURBS曲線包絡(luò)條件，得到了對(duì)應(yīng)的陰陽轉(zhuǎn)子型線，生成的型線誤差數(shù)據(jù)見表3。

圖3 復(fù)盛型線轉(zhuǎn)子型線Figure 3 Rotor lines of Fusheng Profile

分析表3中數(shù)據(jù)可以發(fā)現(xiàn)：嚙合線的最大精度誤差出現(xiàn)在A0B0段，為0.092 μm；陽轉(zhuǎn)子的最大精度誤差在D1E1段，為0.442 μm；陰轉(zhuǎn)子的最大精度誤差出現(xiàn)在D2E2段，為0.443 μm?？梢钥闯鲭m然各部分都存在精度誤差，但都控制在0.5 μm以下，是能接受的范圍。所以該實(shí)例驗(yàn)證了結(jié)合NURBS曲線和嚙合線法設(shè)計(jì)轉(zhuǎn)子型線方法的可行性。

圖4 復(fù)盛型線嚙合線Figure 4 Meshing lines of Fusheng Profile

表1 復(fù)盛型線各段參數(shù)取值Table 1 Value of the parameters of the Fusheng Profile

表2 復(fù)盛型線設(shè)計(jì)相關(guān)數(shù)據(jù)Table 2 Design relevant data of Fusheng Profile

圖5 獲取復(fù)盛嚙合線的NURBS曲線參數(shù)過程

Figure 5 Process of obtaining the NURBS curve parameters of the complex meshing line is obtained

表3 NURBS曲線嚙合線誤差數(shù)據(jù)Table 3 Error data of NURBS curve meshing line

3.3 嚙合線的優(yōu)化

根據(jù)NURBS曲線的局部修改性，通過移動(dòng)控制點(diǎn)可以調(diào)整嚙合線形狀，對(duì)嚙合線進(jìn)行優(yōu)化。對(duì)嚙合線中b、c、d 3點(diǎn)進(jìn)行了微調(diào)，陰陽轉(zhuǎn)子的齒曲線發(fā)生相應(yīng)的變化，見圖6；轉(zhuǎn)子型線性能參數(shù)發(fā)生相應(yīng)的變化，其規(guī)律總結(jié)見表4[13]。所以可以根據(jù)這些變化規(guī)律對(duì)嚙合線的線段進(jìn)行修改，以達(dá)到優(yōu)化型線性能的目的。

在不改變?cè)修D(zhuǎn)子結(jié)構(gòu)的情況下，如陰陽轉(zhuǎn)子的齒頂圓弧大小、齒數(shù)比等參數(shù)，主要以減小泄漏三角形的面積和增大面積利用系數(shù)為目標(biāo)對(duì)型線進(jìn)行優(yōu)化。即不移動(dòng)圖6中的c點(diǎn)位置，變化b和d 2點(diǎn)的位置。根據(jù)表4的規(guī)律，將嚙合線的控制點(diǎn)a5、a6和b1～b4向外移動(dòng)，使得嚙合線交點(diǎn)下移，但不能超過陰陽轉(zhuǎn)子齒頂圓的范圍；然后移動(dòng)控制點(diǎn)c2～c5，使得嚙合線最高點(diǎn)下移，控制點(diǎn)移動(dòng)對(duì)比見圖7。優(yōu)化嚙合線后，陰陽轉(zhuǎn)子齒曲線也發(fā)生了相應(yīng)的變化，優(yōu)化前后的對(duì)比見圖8。陣列優(yōu)化后的陰陽轉(zhuǎn)子齒曲線得到完整的轉(zhuǎn)子型線，見圖9。

圖6 嚙合線頂點(diǎn)與陰陽轉(zhuǎn)子型線變化關(guān)系

Figure 6 Relationship between the vertices of the meshing line and the yin-yang rotor line

表4型線性能參數(shù)與各頂點(diǎn)移動(dòng)方向關(guān)系?

Table 4 Relationship between the linear performance parameters and the direction of each vertex movement

性能參數(shù)移動(dòng)方向bcd齒心角增加 陰/陽+x/+y-x+x/-y齒厚增加 陰陽-x/-y+x/+y+x-x-x/+y+x/-y面積利用系數(shù)增加/+x/+y-x+x/-y

? “陰”表示陰轉(zhuǎn)子，“陽”表示陽轉(zhuǎn)子，“+x”表示向x軸正方向移動(dòng)，“-x”表示向x軸負(fù)方向移動(dòng)，“+y”、“-y”擁有類似的含義。

將優(yōu)化前后的轉(zhuǎn)子性能參數(shù)進(jìn)行整理得到表5。由表5可知，優(yōu)化后的型線增加了陰轉(zhuǎn)子的厚度，陰轉(zhuǎn)子的齒間面積減少了16.212 mm2，但陽轉(zhuǎn)子齒間面積增加了32.678 mm2，所以面積利用系數(shù)提高了1%，而泄漏三角形的面積減少了5%。綜上，優(yōu)化后的轉(zhuǎn)子性能得到了提升。

4 結(jié)論

(1) 利用龍貝格求積公式來解決嚙合線法中原函數(shù)不能用初等函數(shù)表達(dá)的問題，實(shí)現(xiàn)了NURBS曲線的反向型線設(shè)計(jì)。

圖7 型線改進(jìn)優(yōu)化前后控制點(diǎn)對(duì)比

Figure 7 Comparison of control points before and after optimization of meshing lines

圖8 型線優(yōu)化前后對(duì)比Figure 8 Comparison of line optimization before and after optimization

圖9 優(yōu)化后的陰陽轉(zhuǎn)子型線Figure 9 Yin-yang rotor profile after optimization

表5 優(yōu)化前后性能參數(shù)對(duì)比Table 5 Comparison of performance parameters before and after optimization

(2) 針對(duì)典型型線復(fù)盛型線的嚙合線進(jìn)行反向設(shè)計(jì)，嚙合線用一整段NURBS曲線表達(dá)，只由控制點(diǎn)、權(quán)因子和節(jié)點(diǎn)矢量決定，得到了誤差在0.5 μm以內(nèi)的型線，驗(yàn)證了用NURBS進(jìn)行型線設(shè)計(jì)的可行性。增添了雙螺桿壓縮機(jī)轉(zhuǎn)子型線設(shè)計(jì)的曲線類型。

(3) 根據(jù)嚙合線頂點(diǎn)與型線性能的變化關(guān)系，利用NURBS曲線的局部修改性對(duì)型線進(jìn)行優(yōu)化。最終得到優(yōu)化后的型線，增加了陰轉(zhuǎn)子齒厚，增大了面積利用系數(shù)，減小了泄漏三角形的面積，整體性能得到了提升。

(4) 結(jié)合NURBS曲線和“嚙合線法”反向設(shè)計(jì)轉(zhuǎn)子型線，優(yōu)勢(shì)在于NURBS曲線的便捷的局部調(diào)整性，目前該方法只通過調(diào)整控制點(diǎn)的方式調(diào)控曲線，為了最大化利用NURBS曲線的形狀調(diào)整工具，今后應(yīng)將調(diào)整權(quán)因子的手段添加進(jìn)來。

(5) 螺桿壓縮機(jī)轉(zhuǎn)子對(duì)加工精度要求非常高，所以在制造過程中應(yīng)選擇高精度的銑削或磨削加工方式，在加工過程中應(yīng)特別關(guān)注可避免的加工誤差的產(chǎn)生，如刀具的安裝、磨損等產(chǎn)生的加工誤差。

[1] 郭耀君, 謝晶, 朱世新, 等. 超低溫制冷裝置的研究現(xiàn)狀和進(jìn)展[J]. 食品與機(jī)械, 2015, 31(1): 238-243.

[2] STOSIC N, SMITH I K, KOVACEVIC A. A twin screw combined compressor and expander for CO2refrigeration systems[C] // STOSIC N. International Compressor Engineering Conference, Proceedings of the 2002. Purdue: Purdue University, 2002: 1 591-1 594.

[3] 傅廷振, 林雨晨. 無油空壓機(jī)技術(shù)為食品行業(yè)帶來純凈氣源[J]. 食品安全導(dǎo)刊, 2016(7): 72-73.

[4] 阿特拉斯-科普柯公司. 阿特拉斯.科普柯公司在中國推出全新系列高效能壓縮機(jī)[J]. 包裝與食品機(jī)械, 2012, 30(6): 70.

[5] STOSICN, SMITH I K, KOVACEVIC A. Screw compressors: Mathematical modeling and performance calculation[M]. Berlin: Springer, 2005: 32-34.

[6] SU Shyh-haur, TSENG Ching-huan. Synthesis and optimization for rotor profiles in twin rotor screw compressors[J]. Journal of Mechanical Design, 2000, 122(4): 543-552.

[7] WU Yu-ren, FONG Zhang-hua. Rotor profile design for the twin-screw compressor based on the normal-rack generation method[J]. Journal of Mechanical Design, 2008, 130(4): 42 601-42 608.

[8] 何雪明, 戴進(jìn), 劉洪園. 基于自由曲線的螺桿轉(zhuǎn)子型線的正反向設(shè)計(jì)[J]. 中國機(jī)械工程, 2012, 23(22): 2 752-2 757.

[9] ZAYTSEV D, FERRIRA C A I. Profile generation method for twin screw compressor rotors based on the meshing line[J]. International Journal of Refrigeration, 2005(28): 744-755.

[10] ZAYTSEV D, FERRIRA C A I. Screw compressor rotor profile generation method based on pre-defined meshing line[C]// 3rd International Compressor Technique Conference(ICTC 2001). Wuxi: [s.n.], 2001: 117-127.

[11] 卜嘯華, 鄒友金, 孟爾鏞. 螺桿壓縮機(jī)型線設(shè)計(jì)的逆命題法[J]. 上海理工大學(xué)學(xué)報(bào), 1986(3): 43-51.

[12] 何雪明, 童潔, 戴進(jìn), 等. 雙螺桿轉(zhuǎn)子型線的反向優(yōu)化設(shè)計(jì)及應(yīng)用[J]. 機(jī)械設(shè)計(jì)與研究, 2015, 31(4): 151-157.

[13] 何雪明, 姜振鋼, 武美萍, 等. 雙螺桿壓縮機(jī)轉(zhuǎn)子嚙合線與型線關(guān)系研究[J]. 中國機(jī)械工程, 2017, 28(7): 763-770.

[14] PIEGL L A, TILLER W. 非均勻有理B樣條[M]. 2版. 趙罡, 穆國旺, 王拉柱, 等, 譯. 北京: 清華大學(xué)出版社, 2010: 86-93.

[15] 邢誠, 王建強(qiáng), 賈志強(qiáng), 等. 多種數(shù)值積分方法比較分析[J]. 城市勘測(cè), 2010(1): 104-106.

[16] 劉玉娟, 陳應(yīng)祖, 盧克功. 龍貝格積分法及其應(yīng)用編程[J]. 重慶科技學(xué)院學(xué)報(bào): 自然科學(xué)版, 2007, 9(1): 97-99.

[17] LEE H T. Screw-rotor machine with an ellipse as a part of its male rotor: US, A4890992[P]. 1990-01-02.