劉金鋼 卜 昆 楊小寧 程云勇 董一巍

西北工業(yè)大學(xué)現(xiàn)代設(shè)計(jì)與集成制造技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室，西安，710072

0 引言

渦輪葉片是航空發(fā)動機(jī)的關(guān)鍵部件之一，由于其內(nèi)部冷卻結(jié)構(gòu)復(fù)雜，且長期服役在復(fù)雜循環(huán)熱應(yīng)力環(huán)境條件下，故對其氣動外形、性能強(qiáng)度及尺寸精度要求高，而葉片的壁厚是關(guān)系到其性能和精度的一個(gè)重要參數(shù)，因此對葉片壁厚的準(zhǔn)確測量有著非常重要的意義。在葉片的檢測方法中，傳統(tǒng)的方法主要有超聲波法［1］和電渦流法，其中超聲波法受超聲波測厚儀影響，存在測量盲區(qū)和測量下限的限制，而電渦流法受渦流測厚儀法向找準(zhǔn)、測頭半徑補(bǔ)償，以及工件尺寸大小等因素的影響，測量精度低，無法達(dá)到公差要求。在實(shí)際壁厚測量中，一般需要測量葉身的3～5個(gè)關(guān)鍵截面，大多使用超聲波測厚儀沿葉片表面測量點(diǎn)的法線方向進(jìn)行測量，該方法不易確定法線方向，誤差較大，容易放大壁厚。文獻(xiàn)［2－5］提出的基于ICT（industry computerized tomography）圖像的壁厚檢測方法，通過像素標(biāo)定技術(shù)對截面內(nèi)外輪廓進(jìn)行測量，其測量精度較高，但花費(fèi)時(shí)間長，成本高。

為了克服傳統(tǒng)壁厚測量的不足，本文以CT（computerized tomography）測量數(shù)據(jù)為基礎(chǔ)，提出一種基于中弧線的空心渦輪葉片壁厚計(jì)算方法，采用等距法計(jì)算外輪廓中弧線，根據(jù)壁厚定義計(jì)算空心葉片壁厚。與傳統(tǒng)方法相比，該方法測量精度高，適用范圍廣。

1 CT數(shù)據(jù)的采集和處理

空心渦輪葉片的壁厚是指內(nèi)外輪廓之間的厚度。針對渦輪葉片內(nèi)輪廓無法測量的特點(diǎn)，采用實(shí)驗(yàn)室的錐束CT測量設(shè)備進(jìn)行切片圖像采集，通過檢測分析軟件對切片圖像進(jìn)行圖像處理，可得到渦輪葉片的點(diǎn)云模型。

1.1 葉片CT數(shù)據(jù)與葉片CAD配準(zhǔn)

葉片CT測量數(shù)據(jù)配準(zhǔn)定位精度的高低直接影響壁厚測量精度，要準(zhǔn)確測量空心葉片的壁厚必須選擇精度高的配準(zhǔn)方法。目前針對配準(zhǔn)問題主要有三類方法：基于標(biāo)記的配準(zhǔn)方法、基于特征的配準(zhǔn)方法［6］、基于表面點(diǎn)集的配準(zhǔn)方法。根據(jù)渦輪葉片點(diǎn)云數(shù)據(jù)的特點(diǎn)，本文采用第三類方法，選擇SVD＿ICP算法。該算法是一種經(jīng)典的模型匹配算法［7］。首先在模型之間建立一些最近點(diǎn)的點(diǎn)對關(guān)系，然后以距離平方和為目標(biāo)函數(shù)，不斷地迭代，求解兩模型之間的旋轉(zhuǎn)矩陣R和平移矩陣T，使得目標(biāo)函數(shù)取得最小值，最后使兩個(gè)模型緊密貼合。

1.2 截面線獲取

經(jīng)過配準(zhǔn)后，渦輪葉片點(diǎn)云模型的測量坐標(biāo)系統(tǒng)與渦輪葉片CAD模型的設(shè)計(jì)坐標(biāo)系統(tǒng)對應(yīng)起來，要計(jì)算渦輪葉片對應(yīng)點(diǎn)處厚度，必須獲得對應(yīng)點(diǎn)處的截面。采用Geomagic Qualify 11軟件的平面截取功能，得到渦輪葉片的截面點(diǎn)云數(shù)據(jù)，運(yùn)用二次提取法［8－10］對截面線數(shù)據(jù)點(diǎn)進(jìn)行排序、粗提取和精提取等處理，最后得到CT截面線輪廓，如圖1所示。

圖1 CT截面線輪廓提取

2 葉片壁厚計(jì)算

2.1 空心渦輪葉片壁厚相關(guān)定義

為了減小誤差，提高精度，對內(nèi)外輪廓的截面線進(jìn)行擬合，計(jì)算外輪廓的中弧線，求取內(nèi)輪廓（不包括加強(qiáng)肋的外輪廓）上的點(diǎn)P2到中弧線的最短距離的點(diǎn)O，以點(diǎn)O為圓心，作點(diǎn)O的內(nèi)切圓，內(nèi)切圓與外輪廓的切點(diǎn)為P1，連接P1和P2，P1和P2間的距離即為切點(diǎn)P1處的壁厚，如圖2所示。

圖2 空心葉片壁厚定義示意圖

2.2 基于中弧線計(jì)算壁厚

白瑀等［11］、王軍偉［12］提出了兩種計(jì)算中弧線的方法：等半徑法和等角度法。這兩種方法都要計(jì)算切線或法線，要是沒有曲線的連續(xù)模型，就必須用差商來代替，而這樣處理誤差會較大。采用內(nèi)切圓法計(jì)算時(shí)，法線方向可能與葉盆（或葉背）曲線存在多個(gè)交點(diǎn)，采用等角度法計(jì)算時(shí)，必須將前一個(gè)內(nèi)切圓的半徑作為下一個(gè)內(nèi)切圓半徑的初始值，存在誤差傳播。中弧線上的點(diǎn)一定是截面線上某兩個(gè)不同點(diǎn)的法線的交點(diǎn)，而且它到這兩個(gè)點(diǎn)的距離相等，因此中弧線上的點(diǎn)就是截面線等距線的自交點(diǎn)?；诖?，本文采用等距線的方法求解中弧線。

利用等距線法求中弧線，首先運(yùn)用二次提取法得到截面線輪廓，然后通過調(diào)整葉盆和葉背偏置距離得到等距線交點(diǎn)，最后對這些葉盆和葉背交點(diǎn)進(jìn)行排序、插值得到光順的中弧線。具體求解步驟如下。

（1）由于有理三次Bezier樣條曲線可以精確地表示二次曲線和計(jì)算等距線，因此將葉身外截面線的葉背Cb（θ）離散為n＋1個(gè)插值點(diǎn)，構(gòu)造的n段G2連續(xù)分段有理三次Bezier樣條曲線記為Bi（u）（i＝1，2，…，n），其沿內(nèi)法線方向的單位法向量為Ni（u），則外截面線的葉背等距線為

式中，ρ為偏置距離。

（2）同理，將葉身外截面線的葉盆Cp（θ）離散為n＋1個(gè)插值點(diǎn)，構(gòu)造的n段G2連續(xù)分段有理三次Bezier樣條曲線記為Pi（u），其沿外法線方向的單位法向量記為Wi（u），則外截面線葉盆等距線為

（3）求解下式，以判斷外輪廓線葉身和葉背等距線是否有交點(diǎn)：

分三種情況分析：①當(dāng)偏置距離ρ小于最小壁厚時(shí)，沒有交點(diǎn)。② 當(dāng)取ρ為最小壁厚時(shí)，令ρ←ρ＋Δρ，其中，Δρ為偏置距離的增大步長，可取為葉片數(shù)據(jù)測量精度的1／20。增大偏置距離ρ，求解式（3）得到交點(diǎn)，該點(diǎn)即為內(nèi)切圓圓心，此時(shí)ρ即為內(nèi)切圓半徑。③當(dāng)偏置距離ρ增大到一定程度時(shí)，葉盆和葉背的等距線就不存在交點(diǎn)，此時(shí)ρ應(yīng)該是葉片最大厚度ρmax，可用葉片設(shè)計(jì)厚度ρmax作為偏置距離的上限，當(dāng)ρ≥ρmax時(shí)，求解結(jié)束。

（4）將求解方程得到的所有內(nèi)切圓圓心記為Oj（j＝1，2，…，m），經(jīng)過處理后得到中弧線，如圖3所示。將內(nèi)切圓半徑記為Rj，求解下式得到內(nèi)切圓與截面線葉背Cb（θ）切點(diǎn)：

切點(diǎn)記為Pbj（xbj，ybj）。求解下式得到內(nèi)切圓與截面線葉盆Cp（θ）切點(diǎn)：

此切點(diǎn)記為點(diǎn)Ppj（xpj，ypj）。

圖3 處理后的中弧線示意圖

（5）將內(nèi)輪廓截面線離散成大量的數(shù)據(jù)點(diǎn)，記為Nk（xk，yk）（k＝0，1，…，h），求解點(diǎn)Nk到中弧線的最短距離Dk，得到中弧線上的點(diǎn)，即內(nèi)切圓圓心Oj，根據(jù)步驟（4）的切點(diǎn)Pbj（xbj，ybj）和Ppj（xpj，ypj），求解下式：

如果Lj1≥Lj2，則Lj2為Pbj（xbj，ybj）點(diǎn)的壁厚；如果Lj1＜Lj2，則Lj1為Ppj（xpj，ypj）點(diǎn)的壁厚。

3 空心渦輪葉片壁厚算法驗(yàn)證與分析

本文討論的基于中弧線的壁厚計(jì)算方法，是在UG平臺下基于UG open API函數(shù)庫編程實(shí)現(xiàn)的。以空心渦輪葉片的CT點(diǎn)云數(shù)據(jù)為對象，分別選擇五個(gè)截面進(jìn)行壁厚計(jì)算，得到相應(yīng)的壁厚，其三維分布如圖4所示。

圖4 某復(fù)雜空心渦輪葉片各截面壁厚三維分布圖

其五個(gè)截面的壁厚分布如圖5～圖9所示。由圖5～圖9可知，空心渦輪葉片的壁厚在葉背凸峰和葉盆凹峰處最大，向前緣和后緣方向依次減小，在沒有內(nèi)輪廓的地方突然增大，這是由于存在加強(qiáng)肋的緣故。

圖5 某復(fù)雜空心渦輪葉片截面1壁厚分布圖

圖6 某復(fù)雜空心渦輪葉片截面2壁厚分布圖

圖7 某復(fù)雜空心渦輪葉片截面3壁厚分布圖

圖8 某復(fù)雜空心渦輪葉片截面4壁厚分布圖

圖9 某復(fù)雜空心渦輪葉片截面5壁厚分布圖

從表1的對比結(jié)果可知，在后緣位置法線方法誤差較大，中弧線方法比法線方法平均誤差減小了1.077%，綜合比較得出中弧線方法與實(shí)際測量結(jié)果相對誤差較小，與實(shí)際結(jié)果一致性較好。

表1 渦輪葉片壁厚計(jì)算結(jié)果對比表

4 結(jié)語

本文介紹了基于CT測量數(shù)據(jù)的空心渦輪葉片壁厚計(jì)算實(shí)施過程，提出了一種基于中弧線的空心葉片壁厚計(jì)算方法，通過對空心葉片CT切片數(shù)據(jù)的采集與處理，建立了相應(yīng)的點(diǎn)云數(shù)據(jù)測量模型，并與葉片CAD模型進(jìn)行配準(zhǔn)，實(shí)現(xiàn)了空心葉片各截面線的提取和處理，進(jìn)而對壁厚計(jì)算進(jìn)行了詳細(xì)介紹。使用本文介紹的壁厚計(jì)算方法可以在葉片（主要是葉身）任意橫截面上計(jì)算渦輪葉片的壁厚。與傳統(tǒng)方法相比，本文方法在計(jì)算范圍、精度和自動化程度方面對提高空心葉片壁厚精度及可靠性具有一定的實(shí)用性。

［1］李澤，李由佳.某機(jī)高壓渦輪空心葉片超聲波測壁厚［J］.紅旗技術(shù)，1997（4）：33－36.

［2］程云勇，張定華，卜昆，等.基于工業(yè)CT測量數(shù)據(jù)的空心渦輪葉片三維壁厚分析［J］.機(jī)械強(qiáng)度，2009，31（5）：791－796.

［3］王慶勝，張定華，程云勇.基于ICT切片圖像的航空發(fā)動機(jī)渦輪葉片壁厚檢測［J］.計(jì)算機(jī)工程與應(yīng)用，2005（2）：203－204.

［4］傅健，路宏年，侯濤.基于ICT圖像的航空發(fā)動機(jī)渦輪葉片壁厚尺寸精密測量方法［J］.航空學(xué)報(bào)，2006，27（2）：325－330.

［5］程云勇，張定華，毛海鵬，等，一種基于工業(yè)CT的航空發(fā)動機(jī)渦輪葉片生產(chǎn)檢測系統(tǒng)關(guān)鍵技術(shù)研究［J］.制造技術(shù)與機(jī)床，2004（1）：27－30.

［6］胡明昊，任明武，楊靜宇.一種快速實(shí)用的特征點(diǎn)匹配算法［J］.計(jì)算機(jī)工程，2004，30（9）：23－27.

［7］劉晶.葉片數(shù)字化檢測中的模型配準(zhǔn)技術(shù)及應(yīng)用研究［D］.西安：西北工業(yè)大學(xué)，2006.

［8］徐岳琳.空心渦輪葉片蠟?zāi)１诤窆鈱W(xué)檢測技術(shù)研究［D］.西安：西北工業(yè)大學(xué)，2011.

［9］Huang M C，Tai C C.The Pre－processing of Data Points for Curve Fitting in Reverse Engineering［J］.Int.J.Adv.Manuf.Technol.，2000，16（9）：635－642.

［10］Li D Z，Campbell J，Li Y Y.Filling System for Investment Cast Ni－base Turbine Blades［J］.Journal of Materials Processing Technology，2004，148（3）：310－316.

［11］白瑀，張定華，任軍學(xué)，等.葉片高質(zhì)量造型方法研究［J］.機(jī)械科學(xué)與技術(shù)，2003（3）：447－449.

［12］王軍偉.葉片類曲面造型中的參數(shù)網(wǎng)格優(yōu)化技術(shù)研究［D］.西安：西北工業(yè)大學(xué)，2003.