田 野，寧 濤，陳志同

（1.北京航空航天大學機械工程及自動化學院復雜曲面先進加工技術實驗室，北京 100191；2.北京航空航天大學機械工程及自動化學院飛行器制造工程實驗室，北京 100191）

復雜曲面類零件被廣泛應用于航空發(fā)動機、直升機螺旋槳、汽輪機、汽車精密模具等關鍵部件的生產(chǎn)制造中。航空發(fā)動機具有高復雜度和高精密性的特點，其制造技術是國防實力的重要衡量標準。航空發(fā)動機葉片的質(zhì)量很大程度上決定了發(fā)動機的核心性能。葉片曲面的光順造型對后續(xù)的加工流程極為重要，光順的曲面造型更有利于生成平滑的刀軌，提高加工精度與加工效率。因此葉片的曲面造型不是僅停留在對數(shù)據(jù)點的簡單擬合來獲得其大致外形，更重要的是對數(shù)據(jù)點進行一系列的優(yōu)化，自動快速地獲得光順曲面。國內(nèi)外學者關于曲面的光順造型開展了一系列研究。

在葉片光順曲面造型方面，根據(jù)葉片幾何特征間約束的表達可分為基于曲線變形的曲面造型方法與基于曲線擬合的曲面造型方法。曲線變形方法包含七段圓弧法、非剛性配準法、分段高斯映射法等。七段圓弧法最早由Mohaghegh 等[1]提出，通過單層二維截面上數(shù)據(jù)點的劃分來拼接構造七段圓弧以獲得完整截面輪廓線。從特征約束的角度出發(fā)，Mohaghegh 等[2]提出了葉片曲面重構的3 個基本步驟來構造葉片曲面。Piya 等[3]提出分段高斯映射法并在機翼上應用且取得了成功。Li等[4]提出的非剛性變形配準方法可在一定的精度條件下重構多層截面輪廓線。Zhao 等[5]提出基于自由變形的葉片截面線重構方法重構截面輪廓線。曲線擬合方法中較有代表性的方法為Werghi 等[6–7]提出的直線和圓弧特征在葉片截面輪廓線中的約束擬合法。Ke 等[8]提出了基于B 樣條等自由曲線的約束擬合方法。Khameneifar 等[9]給出了基于數(shù)據(jù)點不確定度的最小波動線擬合方法。

在曲面造型這一領域，國內(nèi)也進行了相當多的研究工作。莫堃等[10]提出了結合流場約束的隱式曲面重構方法，保證了重構曲面的幾何精度要求。呂學庚[11]提出一種航空發(fā)動機葉片流曲面重構方法，降低了葉片曲面重構的不確定性。曲學軍等[12]利用統(tǒng)計學知識給出了平面輪廓線數(shù)據(jù)的B 樣條曲面擬合算法并對平面輪廓線掃描數(shù)據(jù)進行B 樣條曲面擬合。張云[13]給出了一種改進的中弧線–厚度分布方法，可用于精鍛葉片自適應重構。施恒濤等[14]提出了一種基于多項式的曲率連續(xù)前緣造型方法，通過給定葉片前后緣和葉身交點過渡部分的中弧線和厚度分布的各階導數(shù)來保證截面輪廓線曲率連續(xù)。周承號等[15]對點云數(shù)據(jù)進行數(shù)據(jù)處理，并對數(shù)據(jù)進行了流場的數(shù)據(jù)模擬，完成了渦輪葉片的曲面造型。黃攀[16]基于葉片二維翼型及三維葉片泛函集成表達理論,利用NURBS 節(jié)點控制向量構建了風力機葉片曲面模型。

上述研究中依舊存在一些問題，具體表現(xiàn)在：（1）通常將每層截面輪廓線作為整體來進行優(yōu)化，生成的整張封閉曲面往往需要在調(diào)整數(shù)據(jù)點位置時，同時控制其逼近精度，導致算法異常復雜，影響在實際工程中的應用；（2）忽視了參數(shù)均勻性對葉片曲面造型的影響，除了沿截面輪廓線的周向參數(shù)外，葉片曲面沿徑向的參數(shù)均勻性也會很大程度上影響曲面質(zhì)量。

為了進一步解決由于參數(shù)不均勻?qū)е碌那娌还忭槅栴}，本文提出了一種基于高參數(shù)均勻性的葉片曲面光順造型方法：（1）通過生成4 張曲面以代替整張封閉曲面，實現(xiàn)了對葉片造型過程的解耦，避免葉片邊緣過渡部分與葉身的過度約束；（2）以弧長參數(shù)化方式構造關鍵曲線，在保證造型精度的前提下使多層截面輪廓線光順，達到葉片周向和徑向同時參數(shù)均勻?qū)R的效果，提升了曲面的光順程度；（3）用戶輸入數(shù)據(jù)點，應用優(yōu)化算法即可高效自動地生成光順葉片曲面，極大程度上簡化了葉片的建模流程。通過對某型號發(fā)動機葉片數(shù)據(jù)點進行曲面造型與實際加工，對提出方法進行了驗證。

1 基于高參數(shù)均勻性的葉片曲面光順造型方法

在航空發(fā)動機葉片模型構造過程中，可能會面臨兩方面問題。（1）初始數(shù)據(jù)點通常沿徑向分層分布，每一層數(shù)據(jù)點均位于同一平面上，層與層之間相互平行。在構建葉片截面輪廓線時，如直接用封閉曲線擬合，在首尾相接處極易出現(xiàn)不光順現(xiàn)象，邊緣部分曲率變化越劇烈，得到封閉曲線的精度越低，進而導致構造的曲面發(fā)生扭曲或自交叉。（2）目前曲線的參數(shù)化方式多為積累弦長參數(shù)化，在數(shù)據(jù)點曲率變化較大的部分，實際弦長與實際弧長的差值變大，故需要對積累弦長參數(shù)化方法進行改進以提高參數(shù)均勻性，具有高均勻性的參數(shù)可以保證葉片曲面各層截面輪廓線的周向和徑向參數(shù)準確對齊，避免曲面起伏或局部厚度變化而導致整體造型結果不理想，本文提出一種葉片曲面光順造型方法。算法輸入為葉片初始數(shù)據(jù)點信息，通過數(shù)據(jù)預處理進行初步區(qū)域劃分，進而利用樣條擬合以及非線性方程組的迭代求解獲取關鍵曲線參數(shù)，進而結合中弧線–變半徑厚度圓的曲線包絡理論生成光順葉片曲面，整體構造流程如圖1所示，需要說明的是，每個前置條件衍生出的所有后置條件均需全部滿足，才能進入下一流程。

圖1 葉片光順曲面造型流程Fig.1 Process of blade smooth surface modeling

1.1 光順葉片截面輪廓線的快速生成

1.1.1 誘導曲線的構造

本文采用中弧線–厚度分布函數(shù)的形式來生成葉片截面輪廓線，在輸入僅有數(shù)據(jù)點的前提下，首先將數(shù)據(jù)點大致分為葉盆、葉背、前緣、后緣4 個部分，葉盆曲線與葉背曲線統(tǒng)稱為葉身曲線，對前緣和后緣部分的數(shù)據(jù)點分別進行最小二乘圓弧擬合，得到其擬合圓圓心；接著對葉身部分的數(shù)據(jù)點進行樣條擬合，以獲得葉身曲線的大致輪廓；最后由兩擬合圓圓心以及其連線的中垂線與葉身曲線的兩個交點的中點生成中弧線的誘導曲線。誘導曲線生成算法具體步驟如下。

（1）將數(shù)據(jù)點分為左右兩個半?yún)^(qū)，需要先獲取這組數(shù)據(jù)點集的最小二乘直線。對于同一平面上的數(shù)據(jù)點集，最小二乘直線是指已知的數(shù)據(jù)點集到該直線距離的平方和最小，即具有最小二乘解。在獲得最小二乘直線LALB后，可以獲取數(shù)據(jù)點集中距離最小二乘直線距離最大的測量點，記為PA，過PA做Y軸的平行線LP，找到數(shù)據(jù)點集中距離LP最近的一點，記為PB；以PA和PB將數(shù)據(jù)點分為左右半?yún)^(qū)，申請兩組容器pntsL和pntsR 來儲存左右半?yún)^(qū)的數(shù)據(jù)。

（2）為將數(shù)據(jù)點分為葉盆、葉背、前緣、后緣4 個區(qū)域，需要選取4 個特征點。首先找到數(shù)據(jù)點的精確公切線，由于僅尋找公切線，會有很多可行解，因此需要加上一個限制條件：兩個切點分別位于左右半?yún)^(qū)的下側。得到的精確公切線CT 和切點PCL、PCR，如圖2所示，選取PCL、PCR作為兩個特征點。利用已獲取的公切線，遍歷數(shù)據(jù)點集，尋找左右半?yún)^(qū)兩端在公切線方向向量上投影的極值點作為另外兩個特征點，文中取x坐標的極大值點PM與極小值點PN，圖3給出了后緣部分的特征點PN。獲取了全部4 個特征點后，將數(shù)據(jù)點集分為4 部分并存儲在4組容器中，后緣、葉盆、前緣和葉背部分數(shù)據(jù)點分別對應存儲容器pnts1、pnts2、pnts3 和pnts4，點區(qū)間的選取規(guī)則為左閉右開。

圖2 截面輪廓線數(shù)據(jù)點的精確公切線Fig.2 Exact common tangent of contour data points

圖3 后緣部分的特征點Fig.3 Feature points of trailing edge

（3）利用步驟（2）存儲容器pnts1和pnts3 中的點進行最小二乘圓弧擬合，得到葉片前緣與后緣部分的最小二乘擬合圓弧。值得注意的是，進行一次擬合很難滿足精度要求，因此需要給定一定的容差δ來進行迭代擬合，容差由葉片造型允許誤差決定。得到滿足精度要求的最小二乘擬合圓弧后，可以獲得兩側的圓弧圓心，記為O1和O2。

（4）利用步驟（2）存儲容器pnts2和pnts4 中的點，以樣條插值方式得到近似的葉身曲線，記為CA和CB。連接O1和O2，得到線段并計算出其中垂線T，中垂線T與CA、CB分別交于PU和PD兩點，取PU和PD的中點，記為PO；以O1、PD、O23 點插值出誘導曲線LG，如圖4所示。

圖4 精確中弧線的獲取Fig.4 Generation of accurate camber curve

1.1.2 弧長參數(shù)化方式光順中弧線的獲取

根據(jù)1.1.1 節(jié)中生成的誘導曲線，利用改進的牛頓迭代法生成初始中弧線，進而對中弧線進行弧長參數(shù)化，最后進行給定容差的光順處理，可以得到容差范圍內(nèi)的光順中弧線以及厚度分布函數(shù)。

誘導曲線的參數(shù)區(qū)間Interval 為u∈[0，1]，在誘導曲線上按照一定的步長取點，過這些點分別作誘導曲線的垂線，內(nèi)切圓圓心一定在垂線上。葉背曲線CA的含參數(shù)表達式為CA（s），葉盆曲線CB的含參數(shù)表達式為CB（t），O（x0，y0）為誘導曲線上的點，將其作為內(nèi)切圓圓心的初始點，內(nèi)切圓半徑為r0，其初始值設為左側最小二乘圓弧的半徑，直線VC為過O點的誘導曲線的垂線且滿足方程Ax0+By0+C=0，Q（s0）與R（t0）分別為VC與CA和CB的交點，將其作為切點初始值；CA'（S0）和CB'（t0）為曲線在該處的切矢，則可列出非線性方程組，即

通過具有預估–校正格式的改進牛頓迭代法可求解上述非線性方程組，初始值上文已經(jīng)給出。通過試驗得知，上文的初始值迭代運算次數(shù)一般在3 次以內(nèi)，可以快速收斂至近似解?？梢缘玫揭幌盗袃?nèi)切圓圓心{Oi}以及內(nèi)切圓半徑{ri}，i=0，1，…，n，其中，n為按給定步長選取的參數(shù)點數(shù)量，進而利用3 次樣條插值的方法獲得精確中弧線，其參數(shù)表達式為MCac（u）。對生成的精確中弧線進行光順處理和弧長參數(shù)化，首先利用改進的能量法，對中弧線進行給定容差的光順處理，得到更加平滑的中弧線，按示例葉片的加工要求，建模允許誤差應比最大允許加工誤差降低一個數(shù)量級，因此其偏差控制在0.001 mm 以內(nèi)，針對不同葉片的加工要求，光順偏差可靈活調(diào)整，圖5給出了光順前后的中弧線曲率梳分析。

圖5 中弧線曲率梳分析Fig.5 Analysis of curvature combs of camber curve

目前曲線的參數(shù)化多為積累弦長參數(shù)化方式，即用弦長近似代替弧長，這種參數(shù)化方式在參數(shù)的末端邊界部分會產(chǎn)生較大的積累誤差。為了消除由于算法本身所帶來的誤差，對光順中弧線進行重新參數(shù)化，即由積累弦長參數(shù)化轉(zhuǎn)變?yōu)榛￠L參數(shù)化。具體分為以下4 個步驟。

（1）計算出光順中弧線的總弧長，然后按照上文給定的參數(shù)步長在中弧線上取點，記錄下這些數(shù)據(jù)點{Di}，以及對應中弧線上的參數(shù){ui}，其中，i=0，1，…，n。

（2）依次計算從中弧線起始點到數(shù)據(jù)點{Di}的弧長與曲線總弧長的比值，以這個比值來代替原先的參數(shù)，記為{ui'}，i=0，1，…，n。

（3）計算{ui'}與{ui}的差值，當所有差值均小于給定容差 （針對文中算例，容差取10–6mm）時弧長參數(shù)化結束，令ui=ui'，否則回到步驟（1），直到滿足差值小于給定容差，其中，i=0，1，…，n。

（4）用給定點串和參數(shù)的方式插值步驟 （3）中最終得到的數(shù)據(jù)點，其參數(shù)和節(jié)點矢量均固定，最終得到弧長參數(shù)化的光順中弧線MC（u）。

1.1.3 光順葉片截面輪廓線生成

在生成弧長參數(shù)化的光順中弧線后，首先利用中弧線–厚度分布函數(shù)的方式構造葉身曲線，然后在葉片截面輪廓線的前緣和后緣部分生成過渡曲線，并保證過渡曲線與葉身曲線G1連續(xù)，最終得到光順的葉片截面輪廓線。具體實現(xiàn)步驟如下。

（1）求解厚度分布函數(shù)。對內(nèi)切圓半徑集合{ri}，其中，i=0，1，…，n以及1.1.2 節(jié)中最終得到的中弧線上的參數(shù)集合{ui}進行3 次樣條插值可以得到r–u曲線，即厚度分布函數(shù)曲線。

（2）厚度分布函數(shù)曲線的光順與弧長參數(shù)化。利用改進的能量法，對步驟 （1）中生成的r–u曲線進行給定容差的光順處理，得到更加平滑的中弧線，其偏差控制在0.001 mm以內(nèi)。利用與1.1.2 中相同的弧長參數(shù)化方法對厚度分布函數(shù)曲線進行弧長參數(shù)化。圖6給出了光順前后的厚度分布函數(shù)曲線以及以曲率為參數(shù)的曲率梳分析。

圖6 厚度分布函數(shù)曲線曲率梳分析Fig.6 Analysis of curvature combs of thickness distribution function curve

（3）生成葉身曲線。利用1.1.2中得到的中弧線MC（u）以及步驟（2）中得到的光順厚度分布函數(shù)曲線生成葉身曲線。首先計算一系列內(nèi)切圓圓心，中弧線起點記為O0MC，其參數(shù)u0MC=0，終點記為O1MC，其參數(shù)u1MC=0，取參數(shù)中點=（u0MC+u1MC），然后得到對應內(nèi)切圓圓心，計算到O0MC與O1MC連線之間的距離，當距離大于給定容差時 （針對文中算例，容差取10–4mm），記分別為O3MC、u3MC，O3MC將中弧線分為兩部分，對這兩部分中弧線重復上述過程，直至所有相鄰3 個內(nèi)切圓圓心均小于給定容差，計算結束。

最終得到內(nèi)切圓圓心集合{OiMC}，以及對應中弧線上參數(shù){uiMC}，i=0，1，…，p。其中，p為最終得到的內(nèi)切圓圓心數(shù)量。

通過上述計算得到的厚度分布函數(shù)曲線方程以及一系列內(nèi)切圓方程可以看作單參數(shù)曲線族F（x，y，c）=0，其中，c為參數(shù)，對應中弧線上的參數(shù)u，F(xiàn)（x，y，c）對x，y，c連續(xù)可微。由包絡曲線的定義可以得到

在式（2）中消去參數(shù)而得到的方程F（x，y）=0，所表示的曲線即為單參數(shù)曲線族F（x，y，c）=0 的包絡曲線。依此求解可得到葉盆和葉背上的點，然后利用樣條插值得到葉身曲線，如圖7所示。

圖7 葉身曲線Fig.7 Profile of blade

（4）生成過渡曲線。在得到葉身曲線后，需要生成前緣和后緣部分的過渡曲線，可得到完整的葉片截面輪廓線；下文以生成葉片的前緣過渡曲線為例，后緣部分的過渡曲線生成方法完全相同。分別取葉身曲線的左端點M和N作為生成過渡曲線的兩個控制頂點，并生成端點在曲線處的切線TM和TN，過1.1.1 節(jié)中步驟 （3）得到的擬合最小二乘圓弧的所有數(shù)據(jù)點做兩條切線的投影，在TM和TN上找到距離葉身左端點最遠的投影點VA與VB，以此作為擬合過渡曲線的另兩個控制頂點。

以4 個控制頂點M、N、VA、VB生成3 次Bezier 曲線作為初始過渡曲線，為保證生成的過渡曲線與葉身曲線的連續(xù)性，將3 次Bezier 曲線進行兩次升階操作，即5 次Bezier 曲線；對5 次Bezier 曲線以1.1.1 節(jié)中步驟（3）得到的擬合最小二乘圓弧的數(shù)據(jù)點為參考點進行迭代的逼近與光順處理，直至過渡曲線與葉身曲線的搭接部分誤差小于10–8mm，以得到與葉身曲線G1連續(xù)的過渡曲線，如圖8所示。

圖8 過渡曲線的生成Fig.8 Generation of transition curve

值得指出的是，若直接用擬合最小二乘圓弧的所有數(shù)據(jù)點進行迭代逼近與光順處理，會面臨點串越界的情況，即部分數(shù)據(jù)點位于M和N的右側，因此在進行逼近和光順之前，對擬合最小二乘圓弧的數(shù)據(jù)點集進行了處理，剔除掉越界的點。

1.2 光順葉片曲面的生成

1.2.1 葉身曲面造型

在得到光順的葉片截面輪廓線后，通過曲面放樣可獲得葉盆和葉背兩張曲面。曲面放樣中沿v方向參數(shù)化的具體方法：對1.1.3 節(jié)中生成的葉身曲線進行反算，得到其B 樣條的控制頂點，如圖9所示，以葉盆曲線的一部分為例，RB 和RB'為相鄰兩層葉盆曲線的一部分，{CPi}和{CPi'}（i=0，1，…，6）為對應的控制頂點，求解相鄰兩層對應的控制頂點之間的實際距離占所有不同層對應控制頂點的實際距離之和的比例作為v方向的控制參數(shù)。由于同一層截面線上的控制頂點均在同一平面上，且經(jīng)過1.1.3 節(jié)中步驟（3）的算法處理，所有的葉身曲線非常接近于弧長參數(shù)化，極大地減小了積累弦長參數(shù)化引起的不同層截面線之間的u方向誤差，因此可以確保v方向的參數(shù)對齊且不會出現(xiàn)參數(shù)曲線扭曲的情況，生成的放樣曲面如圖10所示。

圖9 反算B 樣條的控制頂點Fig.9 Inverse computation of B-splines control points

圖10 放樣曲面Fig.10 Lofted surface

1.2.2 前后緣過渡曲面造型

拾取兩張曲面，兩條脊線和多條過渡曲線進行融合操作。以前緣過渡曲面生成為例，其具體步驟如下。

（1）如圖11所示，拾取葉盆曲面SE 和葉背曲面SX，以兩曲面的左邊界SC 與SC'作為前緣過渡曲面的脊線，相鄰兩層過渡曲線為TC0與TC1。TC0與SC、SC'交于IP0、IP0'并取連線的中點O0作為原點，過O0作平面PV垂直于線段IP0IP0'，過IP0作SC的切線T0，過IP0'作SC'的切線T0'，將T0與T0'求平均并投影到平面PV，投影所得向量為x軸，以平面PV中x軸的垂線方向為z軸，利用右手定則得到y(tǒng)軸建立局部坐標O0xyz。同樣的方法用于TC1建立局部坐標系O1xyz。

圖11 融合曲面Fig.11 Blend surface

（2）取曲線SC 在IP0與IP1之間的部分取參數(shù)中點，得到IP2，對SC'采用相同的策略，獲得IP2'，用步驟（1）中的方法建立局部坐標系O2xyz，將TC0與TC1在O2yz面內(nèi)的投影曲線進行加權平均，并對其添加限制條件，令其兩個端點分別與IP2、IP2'兩點重合，且端點處與SE、SX 與O2yz平面的交線保持G1連續(xù)，利用迭代逼近的方法得到插值曲線TC2。

（3）對全部相鄰的過渡曲線重復步驟 （2），直至滿足密度要求，得到曲線組{TC0，TC1，…，TCn}，其中，n為給定曲線密度，以得到的曲線組生成平滑過渡曲面。

文中算法均在CAM5 曲面造型系統(tǒng)中進行開發(fā)，在葉片的光順造型過程中，用戶選擇航空發(fā)動機葉片的數(shù)據(jù)點文件，系統(tǒng)讀取數(shù)據(jù)點后，在用戶操作下，自動進行葉片曲面光順造型相關操作，包含誘導曲線的構造、弧長參數(shù)化方式光順中弧線的獲取、光順葉片截面輪廓線生成、葉身曲面造型以及前后緣過渡曲面造型，直到完成葉片曲面光順造型。

2 結果分析

2.1 理論分析

假設給出某條曲線的參數(shù)表達式C（u），a≤u≤b，其中，a∈R，b∈R，曲線上的數(shù)據(jù)點集為{pi}，其中，i=0，1，…，n，則按照積累弦長參數(shù)化所得到的參數(shù)區(qū)間為

對于曲線C（u），弧長參數(shù)表示式為

式中，s為弧長參數(shù)；a為曲線初始參數(shù)左邊界；t為積分參數(shù)；s（t）為一個確定的極限，與曲線整體長度、數(shù)據(jù)點集分布、曲線局部的曲率變化率無關，因此其精度高于積累弦長參數(shù)化。本文的優(yōu)化方法構造了弧長參數(shù)化的光順中弧線以及葉身曲線，選擇合適的曲面放樣方法得到葉身曲面。在邊緣過渡部分，由于曲線整體長度較短，且脊線參數(shù)分布合理 （為葉身曲面的邊界），故采取積累弦長參數(shù)化方式，為盡量消除曲率變化率對精度的影響，在參數(shù)化時采取了加密的策略。數(shù)據(jù)點的參數(shù)經(jīng)上述步驟完全確定，在4 張曲面上的數(shù)據(jù)點沿雙方向等參數(shù)線曲率變化平滑，避免了曲面的扭曲現(xiàn)象，達到光順的效果。

2.2 造型實例分析

為驗證本文提出的葉片曲面光順造型方法的正確性，選取具有20層截面輪廓線的某型號航空發(fā)動機葉片數(shù)據(jù)點為例，進行了葉片曲面光順造型試驗驗證，葉片的加工最大允許誤差為–0.03～+0.05 mm。通過對造型結果以及實際加工效果的分析，驗證了所提出方法的正確性。結果分析從以下4 方面進行。

（1）精度分析。計算所有數(shù)據(jù)點到生成曲面的距離，為簡潔起見，圖12中僅展示葉背曲面的100 個數(shù)據(jù)點及后緣過渡曲面的100 個數(shù)據(jù)點。其中數(shù)據(jù)點到葉身曲面的最大距離為0.008915 mm，到前緣后緣過渡曲面最大距離為0.006173 mm，均滿足葉片建模的精度要求。

圖12 數(shù)據(jù)點到曲面的距離Fig.12 Distance of data point to surface

（2）曲面曲率分析。分別對初始葉片曲面和經(jīng)過手工調(diào)整的模型曲面以及利用文中算法生成的光順葉片曲面進行曲率分析。圖13中曲面由左至右分別是UG 初始建模曲面、人工調(diào)整的曲面和光順造型曲面。從分析結果可以看出，UG 直接通過插值曲線和曲線組方式建模，將葉片模型作為一個整體來進行曲面建模，其在前緣和后緣部分有曲率突變并延伸到了葉身部分 （圖13（b）），人工調(diào)整的模型曲面有同樣的問題，而通過文中方法進行的葉片光順曲面建模是由葉盆、葉背、前緣、后緣4個部分構成，通過精確的分區(qū)避免了曲率突變問題，較劇烈的曲率變化被嚴格地限制在邊緣過渡曲面部分，在邊緣和葉身過渡部分曲率變化更加平滑，曲面的整體也更加光順。

圖13 曲面的曲率分析Fig.13 Curvature analysis of surfaces

（3）截面輪廓線曲率梳分析。截面輪廓線的曲率梳分析最能說明曲線的光順性，以曲率為參數(shù)對3 個模型均進行了相同參數(shù)位置的截面輪廓線曲率梳生成。初始模型的截面輪廓線在前后緣過渡部分保持了G1連續(xù)，但葉身曲線曲率變化不平滑；手工調(diào)整模型是根據(jù)加工經(jīng)驗，通過人工修改數(shù)據(jù)點使其符合加工要求，可以看到，其葉身曲率變化不平滑，且在前后緣過渡部分無法保證G1連續(xù)。雖然可以生成較為平滑的刀軌，但是對于曲面建模而言具有很大的風險與不確定性。通過分析對比 （圖14）可以得到結論：應用文中算法進行快速建模的光順葉片曲面模型，其截面輪廓線曲率變化平滑，且在前緣后緣部分保持G1連續(xù)，具有良好的幾何性質(zhì)。

圖14 截面輪廓線曲率梳分析Fig.14 Curvature comb of profile curves

（4）刀軌平滑性分析。應用控制變量法，分別在3 個模型中生成的葉片曲面上相同區(qū)域選取相同的加工參數(shù)生成刀軌，在初始模型中生成的刀軌有回折，從而引起刀軸突變，因此試件加工后表面會有振紋，加工效果很差。而應用文中算法生成的葉片曲面上的刀軌平滑無回折，加工效果較好且滿足加工要求，如圖15和16 所示。由此可以得到結論：利用文中優(yōu)化方法快速自動生成的光順葉片曲面具有良好的參數(shù)均勻性及光順度，減少了刀軌的波動性，從而大幅提高加工質(zhì)量、減少機床磨損。

圖15 刀具軌跡的生成Fig.15 Generation of toolpaths

3 結論

圖16 經(jīng)過加工的葉片F(xiàn)ig.16 Machined blade

（1）本文以航空發(fā)動機葉片為研究對象，對葉片曲面造型中存在的問題進行分析，將葉片特征之間的關系解耦與特征曲面均勻參數(shù)化相結合，提出了基于高參數(shù)均勻性的葉片曲面光順造型方法?？紤]到葉片目前的整體造型過程較為煩瑣，將優(yōu)化算法集合于曲面造型系統(tǒng)中，用戶可以通過簡易操作自動高效地獲取光順的葉片曲面。

（2）本文提出優(yōu)化算法通過生成4 張曲面以代替整張封閉曲面，實現(xiàn)了對葉片造型過程的解耦，避免了葉片邊緣過渡部分和葉身部分的過度約束，在使算法流程簡化的同時保證了造型精度，利于在實際工程中的應用。

（3）在葉片造型過程中實現(xiàn)特征曲面雙方向均勻參數(shù)化，有效提升了曲面的光順程度，便于后續(xù)的加工。在葉片曲面上具有相同v參數(shù)的點可以構成等v參數(shù)線族，在實際加工中可以進一步研究此曲線族與加工方法的關系。文中提出的方法還可以推廣至其他具有復雜曲面特征的典型航空類零件的曲面造型。

（4）在輸入數(shù)據(jù)點后，用戶只需要進行簡單的操作即可應用優(yōu)化算法快速自動地生成光順葉片曲面，解決了葉片曲面造型整體流程煩瑣的問題。將優(yōu)化方法應用于某型號葉片的實際加工，數(shù)據(jù)分析和加工效果驗證了方法的正確性。