樂 英,王 穎

(華北電力大學 能源動力與機械工程學院，河北 保定 071000)

基于NURBS插補算法的汽輪機葉片數(shù)控加工*

樂 英,王 穎

(華北電力大學 能源動力與機械工程學院，河北 保定 071000)

為了提高汽輪機葉片的加工精度，文章通過分析汽輪機葉片的結構特點，采用一種NURBS曲線插補算法對其進行插補仿真。文章所采用算法由兩部分組成：速度規(guī)劃和參數(shù)計算。首先速度規(guī)劃采用簡化的五段S曲線加減速控制方法，保證了葉片在高速加工過程中加速度的連續(xù)，使機床運行平穩(wěn)，然后利用牛頓迭代法來計算插補參數(shù)，得到更精確的插補參數(shù)，進一步提高了葉片的加工精度和加工速度。汽輪機葉片的插補仿真表明，該算法有高的穩(wěn)定性和運算精度，并且使機床振動減小，速度波動小，保證了葉片的高質量加工。

汽輪機葉片；NURBS插補算法；五段S曲線加減速控制；牛頓迭代法

0 前言

汽輪機葉片的結構設計復雜，是一種典型的復雜自由曲面，是汽輪機中的關鍵部件，其加工質量直接影響到機組的運行效率和可靠性[1]。傳統(tǒng)的數(shù)控加工系統(tǒng)，一般采用很多小的離散化的線段或者圓弧逼近曲線，這種方法會造成曲線輪廓不平滑，造成進給速度的不連續(xù)和波動，從而破壞工件表面的光滑性，而且大量的程序增加了CAD/CAM和計算機數(shù)控系統(tǒng)之間的通訊負擔，影響插補的實時性。

目前，研究NURBS插補算法的國內外學者較多。Bedi[2]提出均勻參數(shù)插補算法，將參數(shù)增量視為常數(shù)，首次實現(xiàn)了NURBS曲線的實時直接插補。Shipital和Yang[3-4]分別采用一階和二階泰勒展開式，應用泰勒展開法計算插補參數(shù)使參數(shù)曲線直接插補應用在實踐加工中成為可能，實現(xiàn)了恒定的進給速度插補，提高了計算精度，減小了速度波動。任杰青[5]等提出基于預估誤差補償?shù)腘URBS曲線預估-校正插補算法，降低了運算的復雜度并且利用基于割線法的校正公式減少了迭代計算的次數(shù)，進一步提高了計算精度。孫海洋[6]等在研究了參數(shù)曲線的數(shù)控插補原理后，發(fā)現(xiàn)一些直接插補算法的局限性，提出了利用牛頓迭代法計算插補參數(shù)，減小了進給速度波動，但是參數(shù)初值的計算相當復雜。上述這些算法都沒有考慮加減速過程中的速度、加速度的變化，難以保證插補精度。Nam[7]等在考慮了加速度的約束之后，提出了一種限制加加速度的加減速控制方法，有效控制了輪廓誤差，但是這種方法的實現(xiàn)比較復雜。周德勝[8]等全面考慮了弓高誤差和法向加速度等參數(shù)的影響，提出一種基于NURBS曲線插補的五段S曲線加減速方法，解決了減速點預測的困難，實時性較好。潘海鴻[9]等提出一種七段式非對稱S曲線加減速控制方法，考慮了工程應用的實際問題，提高了加工速度，但是計算復雜。

上述算法大多是應用在自由曲線曲面上，本文以汽輪機葉片為研究對象，根據(jù)汽輪機葉片的結構特點，對其采用NURBS曲線插補算法，該算法首先進行速度規(guī)劃，速度規(guī)劃采用簡化的五段S曲線加減速控制算法，然后進行參數(shù)計算，參數(shù)計算利用牛頓迭代法，迭代初值采用簡單線性遞推方程計算。對汽輪機葉片截面的仿真結果表明，汽輪機葉片采用該算法進行插補保證了速度和加速度的連續(xù)，加加速度實現(xiàn)了有界變化，從而使得機床運行平穩(wěn)，葉片的插補精度得到了提高，進給速度波動限制在理想范圍內，進一步提高了葉片的加工質量。

1 NURBS曲線插補

一條k次NURBS曲線的表達式為[10-11]：

(1)

其中,wi(i=0,1,2,…,n)稱為權或權因子(weights)。di(i=0,1,2,…,n)為控制頂點，與權因子wi(i=0,1,2,…,n)一一對應，順序相連形成控制多邊形。Ni,k(u)是由節(jié)點矢量U=[u0,u1,…,un+k+2]按德布爾-考克斯遞推公式(2)推導的k次規(guī)范B樣條基函數(shù)。

(2)

NURBS插補主要用弦線逼近弧線，這種逼近會造成實際進給速度與理想進給速度之間存在差異，從而導致速度波動，速度波動的大小采用δi表示，如公式(3)所示[12]：

(3)

2 速度規(guī)劃算法

2.1NURBS曲線自適應速度調整

在實時插補過程中，弓高誤差、進給速度以及曲率半徑之間的關系，采用如圖1所示的方法進行計算：

圖1 利用圓弧計算弦高誤差示意圖

ρi為u=ui處的曲率半徑。弓高誤差的計算如下式所示：

(4)

由最大弓高誤差ER限定下的進給速度V(ui)為：

(5)

式(5)表明進給速度V(ui)與弓高誤差ER和曲率半徑ρi有關，其變化關系如下：

(6)

其中，F(xiàn)為編程進給速度。

2.2 五段S曲線加減速規(guī)劃

速度的自適應調整可以得到滿足精度要求的速度，但在曲線曲率半徑變化較大的區(qū)域，比如曲率半徑較小的區(qū)域，速度變化較大，可能會超出機床的加減速能力，對機床產(chǎn)生沖擊[13]。

為調整小曲率半徑區(qū)域曲線進給速度，使得加速度和加加速度實現(xiàn)連續(xù)有界變化，在曲線曲率尖角處采用五段S曲線加減速控制方法。

五段S曲線加減速控制算法將速度規(guī)劃分為加加速段、減加速段、勻速段、加減速段和減減速段五個階段。以減速段為例：

假設t=0時為減速過程的開始時刻，Vs為起始速度，Ve為終止速度，A為最大加速度，J為最大加加速度。

圖2 五段S曲線加速度控制方法減速過程的速度、加速度和加加速度曲線

從圖2可以看出，加減速階段的速度方程曲線為:

(7)

當t=t1時，A1(t1)=-Jt1=-A，有：

(8)

減速階段的速度方程為：

(9)

由于整條速度曲線是連續(xù)的，有：

V1(t)=V2(t)

(10)

將式(7)～(9)代入式(10)可得：

(11)

設減速階段的速度變化量為ΔV，即ΔV=Ve-Vs,當:

(12)

則機床加工的加速度和加加速度能夠達到最大值。若

(13)

則機床加工的加速度不能達到最大值A，設此時達到的加速度值為Ai，由式(11)可得：

(14)

如果

(15)

則機床加工的加加速度不能達到最大值J，設此時達到的加加速度值為Ji，由式(11)可得：

(16)

因此，在計算減速段的速度方程時，先根據(jù)式(12)～式(16)計算出減速階段能達到的加速度和加加速度，然后根據(jù)公式(7)和式(9)得到減速段的方程。加速段過程的計算方法與減速段相同。

2.3 插補參數(shù)計算

為了減小截斷誤差，降低速度波動，本文采用牛頓迭代法計算插補參數(shù)：

(17)

式中，ui—當前插補參數(shù);ui+1—下一個插補參數(shù)。

則問題轉化為求ξ=ξ*，使得：F(ξ*) =0。

牛頓迭代法的計算原理如圖3所示。

圖3 牛頓迭代法示意圖

設F(ξ)=0的近似根是ξk，過點Pk引切線，點Pk是曲線F(ξ)上橫坐標為ξk的點，該切線與ξ軸的交點的橫坐標為ξk+1，將ξk+1作為ξ*的新的近似值，因此可以得到迭代公式：

(18)

通常迭代法的迭代終止條件為：

(19)

ε為迭代精度。在滿足式(19)時，為了同時滿足速度精度的要求，設置速度波動率δ的上限值作為迭代終止條件，并且規(guī)定了最大迭代次數(shù)的值K。

這里采用公式(20)所表示的線性遞推法來計算迭代參數(shù)初值：

ui+1=2.5ui-2ui-1+0.5ui-2

(20)

為了避免出現(xiàn)插補反轉，利用公式：

ui+1=2ui-ui-1=ui+(ui-ui-1)

(21)

式(21)可以保證ui+1>ui,不會出現(xiàn)插補反轉的問題。

3 試驗仿真

以某國產(chǎn)汽輪機組通流部分低壓缸動葉片為例，葉身長為432mm，沿葉高方向有11個關鍵截面，如圖4所示，各截面形狀如圖5所示。

圖4 葉片橫截面分布圖

圖5 葉片11個截面示意圖

本文以葉片中間的截面為研究對象，葉片中間的截面采用NURBS擬合如圖6所示。

圖6 葉片高度z=216mm的截面擬合圖

從圖中可以看出，葉片的曲率半徑變化較大，進汽邊和出汽邊存在曲率半徑很小的尖角。根據(jù)上圖所示的葉片橫截面，利用MATLAB軟件采用上述方法對其進行NURBS曲線插補仿真。仿真時的編程速度F=100mm/s，插補周期T=0.001s，最大加速度A=0.008m/s2，最大加加速度J=0.004m/s3，允許的最大弓高誤差ER=0.001μm，迭代精度E=10-6，最大迭代次數(shù)K=5，最大速度波動率Δ=2×10-4%。

仿真試驗得到的弓高誤差、速度曲線、加速度曲線及加加速度曲線分別如圖7～圖10所示。

圖7 弓高誤差圖

圖8 速度圖

圖9 加速度圖

圖10 加加速度圖

仿真結果表明，本文所采用的插補算法使葉片所有插補點的弓高誤差、進給速度、加速度和加加速度都控制在允許范圍內，從而保證了葉片的加工精度滿足要求，同時使得速度曲線實現(xiàn)平滑過渡，加速度曲線在加速段和減速段實現(xiàn)連續(xù)變化，并且加加速度實現(xiàn)有界變化，保證了葉片加工過程中機床的平穩(wěn)運行，進一步提高了葉片的加工精度和加工效率。

圖11 速度波動率圖

仿真結果得到的速度波動率如圖11所示，從圖中可以看出，本文采用牛頓迭代法進行插補參數(shù)的計算使得所有插補點的速度波動率均控制在10-4%量級范圍內，說明相對于泰勒展開法計算所得的參數(shù)值更準確，速度更平穩(wěn)，有效降低了速度波動。

綜上所述，表明該算法穩(wěn)定好，且運算精度高，減小了機床的振動，保證了加工過程的穩(wěn)定性，大幅度提高了葉片的加工精度和加工效率。

4 結論

汽輪機葉片是汽輪機的重要部分，它的加工精度直接影響到汽輪機的工作效率，為了提高汽輪機葉片的加工質量，本文采用NURBS曲線插補算法對葉片截面進行插補仿真。算法分為兩個部分：速度規(guī)劃和參數(shù)計算。速度規(guī)劃部分采用簡化的5段S曲線加減速控制方法，而參數(shù)計算部分應用牛頓迭代法來計算插補參數(shù)，克服了傳統(tǒng)方法利用泰勒公式進行計算計算量大且不準確的缺點。仿真試驗結果表明，該算法能夠提高葉片的加工效率，保證葉片插補過程中速度和加速度連續(xù)，從而提高了葉片的加工精度，并且使加加速度連續(xù)有界，避免了刀具在運行過程中對機床產(chǎn)生較大的沖擊，保證機床可以平穩(wěn)運行，有效減小了速度波動，從而滿足葉片加工中高質量、高效率的加工要求。

[1] 樂英．提高汽輪機葉片加工效率的方法[J]．汽輪機技術，2012，54(3)：234-236，186．

[2] BEDI S，Ali S，QUAN N． Advanced interpolation techniques for CNC machines[J]．ASME Journal of Engineering for Industry，1993，S115(8)：329-336．

[3] SHIPIYAL NI M，KOREN Y，Lo C C．Real-time curve interpolation[J]．Computer-Aided Design， 1994，26(11)：832-838．

[4] Yang D C H，Kong T．Parametric interpolator for precision CNC machining parametric curves[J]．Computer-Aided Design．1994，26(3)：225-234．

[5] 任杰青，劉凱，趙東標．NURBS基于預估誤差補償?shù)腘URBS曲線插補算法[J]．機械科學與技術，2015，34(8)：1037-1041．

[6] 孫海洋，范大鵬，李玲．一種參數(shù)曲線實時數(shù)控插補計算新方法[J]．國防大學學報，2008，30(3)：122-127．

[7] NAM S H，YANG M Y．A study on a generalized parametric interpolator with real-time jerk-limited acceleration[J]．Computer-Aided Design，2004，36(1)：27-36．

[8] 周勝德，梁宏斌，喬宇．基于NURBS曲線插補的五段S曲線加減速控制方法研究[J]．組合機床與自動化加工技術，2011(4)：37-41．

[9] 潘海鴻，楊增啟，陳琳，等．一種適用于前瞻的高精度7段式非對稱S曲線加減速控制 [J]．機械科學與技術,2015，34(7)，1024-1030．

[10] 董海濤，潘海鴻，黃麗宇，等．一種優(yōu)化的NURBS曲線插補算法[J]．計算機集成制造系統(tǒng)，2014，20(9)：2172-2177．

[11] 施法中．計算機輔助幾何設計與非均勻有理B樣條[M]．北京：北京航空航天大學出版社，1994．

[12] 王允森，蓋榮麗，孫一蘭，等．基于牛頓迭代法的NURBS曲線插補算法[J]．組合機床與自動化加工技術，2013(4)：13-17．

[13] 潘海鴻，楊微，陳琳，等．全程S曲線加減速控制的自適應分段NURBS曲線插補算法[J].中國機械工程，2010，21(2)：190-195．

(編輯 李秀敏)

NURBS Interpolation Algorithm for Turbine Blade in NC Machining

YUE Ying, WANG Ying

(School of Energy, Power and Mechanics Engineering, North China Electric Power University, Baoding Hebei 071000, China)

In order to improve the machining accuracy of turbine blade, a NURBS curve interpolation algorithm is used according to the structure features of the turbine blade，which included two parts of velocity planning and parameter calculating. First five phased S-curve acceleration and deceleration control method was adopted in velocity planning module．The accelerate speed is continuous when blade in high speed machining process using the method and the CNC machine runs smoothly．Then the interpolation parameters are calculated using Newton iterative method. It gets more accurate interpolation parameters，the machining precision and speed of the blade are more improved. The simulation test for interpolation of turbine blade shows that the algorithm has good stability, high computing precision. It can satisfy the requirements of the real-time interpolation，reduce machine tool vibration and limit the speed fluctuation in the real-time interpolation to the ideal level.

turbine blade；NURBS interpolation algorithm；five phased S-curve acceleration and deceleration control；Newton iterative method

1001-2265(2017)01-0078-04

10.13462/j.cnki.mmtamt.2017.01.021

2016-03-24；

2016-05-06

河北省自然科學基金項目(E2014502042)；中央高?；究蒲袠I(yè)務費專項資金項目(11QJ61)

樂英(1971—)，女，河南信陽人，華北電力大學(保定)副教授，博士，研究方向為復雜曲面造型及數(shù)控加工的研究，(E-mail)yueying71@163.com。

TH165;TG659

A