王倩妮，蘇宇航，郭廣平

（中航工業(yè)北京航空材料研究院，北京 100095）

航空發(fā)動機(jī)空心渦輪葉片在高溫、高壓的惡劣環(huán)境中工作，有很高的強(qiáng)度要求，壁厚是保證葉片強(qiáng)度的一個重要參數(shù)，需要進(jìn)行準(zhǔn)確的測量［1］。目前超聲測厚技術(shù)是航空高溫合金空心葉片主要測量手段。隨著航空葉片越來越多地使用單晶材料，而單晶材料各向異性的特點(diǎn)影響了超聲測厚準(zhǔn)確性。另外，對于曲率較大的測量點(diǎn)由于難以接收超聲波反射信號而無法測量壁厚。工業(yè)CT 技術(shù)成為了航空發(fā)動機(jī)空心渦輪葉片壁厚測量的又一選擇。

工業(yè)CT 技術(shù)是通過若干投影數(shù)據(jù)，利用重建算法獲得試件截面圖像的一種無損檢測技術(shù)，其檢測過程不受零件材料、形狀的限制，得到的CT 圖像不受其他區(qū)域干擾，在零件的CT 圖像上可以直觀地進(jìn)行尺寸測量。高質(zhì)量的圖像是高精度工業(yè)CT尺寸測量的基礎(chǔ)，采用的測量方法對于測量結(jié)果的準(zhǔn)確性同樣有很大影響。

1 尺寸測量方法介紹

受工業(yè)CT 系統(tǒng)的固有特性和空間分辨率的限制，被測零件的CT 圖像邊緣并不存在銳利的邊界，只有灰度值的過渡區(qū)域?；诠I(yè)CT 圖像的尺寸測量方法根本在于確定被測零件與背景的邊界。目前常用的尺寸測量方法主要有基于灰度值變化曲線的半高法和和基于圖像分割的半自動測量法。

1.1 半高法

半高法是根據(jù)CT 圖像沿尺寸測量方向的灰度值變化曲線，將曲線波峰和波谷的中間點(diǎn)定為被測零件與背景的分界點(diǎn)進(jìn)行尺寸測量的方法，具體操作步驟為：①在CT 圖像上沿待測尺寸方向劃一條直線S，如圖1所示。②得到沿此直線的灰度值變化曲線，如圖2所示。③取曲線上最大灰度值和兩側(cè)最小灰度值的中間點(diǎn)（A、B 兩點(diǎn)）作為零件與背景的分界點(diǎn)。④測量A、B 兩點(diǎn)的間距作為所測部位的厚度測量值。

圖2 圖1中直線S處灰度值變化曲線圖

圖2為理想的灰度值變化曲線圖，實際的CT圖像由于受到探測器響應(yīng)不一致、射線散射、射線束硬化等因素影響，曲線的波峰、波谷并不像圖2中平直。而且受到零件結(jié)構(gòu)影響，兩側(cè)波谷處灰度值可能不完全一樣，此時可以分別根據(jù)兩側(cè)的灰度值分別確定被測零件與背景的邊界。

1.2 半自動測量法

半自動測量法是根據(jù)圖像分割算法，將CT 圖像中被測零件與背景區(qū)分開，得到被測零件的邊界，在此基礎(chǔ)上進(jìn)行尺寸測量的方法。受被測零件結(jié)構(gòu)及材料復(fù)雜性及多樣性的限制，并沒有一個對所有零件都適用的圖像分割方法。目前常見的圖像分割算法大致可分為：基于邊緣的分割算法、基于閾值的分割算法、基于區(qū)域的分割算法、基于偏微分方程的分割算法［2］。

筆者使用VGStudio MAX 軟件的表面提取工具進(jìn)行圖像分割。VGStudio MAX 是工業(yè)CT 圖像顯示和分析專用軟件，其表面提取工具分為標(biāo)準(zhǔn)和高級兩種模式。標(biāo)準(zhǔn)表面提取方法根據(jù)CT 圖像直方圖，測定背景峰值和材料峰值，自動計算材料邊界的灰度值，是一種基于閾值的圖像分割方法；高級模式是在標(biāo)準(zhǔn)模式的基礎(chǔ)上，對已測定邊界周圍的像素進(jìn)行解析，補(bǔ)償因射束硬化或其他偽影造成的局部偏差，得到更為準(zhǔn)確的邊界。與標(biāo)準(zhǔn)模式相比，高級模式會耗費(fèi)更多的處理時間，可根據(jù)圖像質(zhì)量分析是否需要使用高級模式。

2 工業(yè)CT掃描

筆者使用標(biāo)準(zhǔn)試塊組成了兩種尺寸試驗件：一個模擬葉片封閉內(nèi)腔結(jié)構(gòu)，由厚度0.5～3.5mm 的13片標(biāo)準(zhǔn)試塊組成，結(jié)構(gòu)及各試塊編號如圖3 所示，編號為1～13號的試塊厚度分別為1.00，3.50，3.50，0.75，2.00，2.00，0.50，1.00，1.00，3.50，0.50，0.50，2.00 mm；另一個為由20 片試塊組成（由薄到厚編號為1～20）的校準(zhǔn)試塊，如圖4所示，試塊1～20 號的具體厚度為0.50，0.60，0.70，0.80，0.90，1.00，1.10，1.20，1.30，1.40，1.50，1.60，1.70，1.80，1.90，2.00，2.25，2.45，2.65，3.00 mm。在工業(yè)CT尺寸測量過程中，為了避免CT 圖像質(zhì)量不佳而導(dǎo)致較大的測量誤差，經(jīng)常放置一個校準(zhǔn)試塊與被測零件一起進(jìn)行CT 掃描，必要時用以校正尺寸測量結(jié)果。

將模擬葉片和校準(zhǔn)試塊放置在視場中間和外側(cè)分別進(jìn)行了CT 掃描。工業(yè)CT 系統(tǒng)和掃描主要參數(shù)為：射線源能量為6 MeV，空間分辨率為2LP／mm，圖像矩陣大小為4096×4096，重建視場半徑為300mm，切層厚度為1mm，三代扇束掃描。

3 測量結(jié)果

3.1 半高法測量結(jié)果

半高法測量使用軟件IPT＿IPS，為了減少人工劃線不穩(wěn)定性帶來的誤差，每個試塊測量3次，測得尺寸D1、D2、D3，取三者平均值作為該試塊尺寸測量值D測量＝，測量誤差絕對值，測量結(jié)果見表1～表2。33個試塊工測量值中，6個測量值的誤差大于0.05mm，結(jié)果誤差較大的測量值像中于厚度為0.50mm 的試塊。

表1 模擬葉片測量結(jié)果 mm

表2 校準(zhǔn)試塊測量結(jié)果 mm

3.2 半自動測量結(jié)果

使用VGStudio MAX 軟件進(jìn)行半自動測量，分別用標(biāo)準(zhǔn)模式和高級模式進(jìn)行了圖像分割，再使用VGStudio MAX 自帶的測量工具“卡尺”完成了尺寸測量，測量結(jié)果見表3～表4。表中D測量－標(biāo)準(zhǔn)和α標(biāo)準(zhǔn)為標(biāo)準(zhǔn)分割模式測量結(jié)果，D測量－高級和α高級為高級分割模式測量結(jié)果。根據(jù)測量結(jié)果，基于標(biāo)準(zhǔn)模式圖像分割的測量值有14 個誤差大于0.05 mm，基于高級模式圖像分割的測量值只有4個誤差大于0.05mm，且均出現(xiàn)在厚度為0.50mm 的試塊。

表3 模擬葉片測量結(jié)果 mm

表4 校準(zhǔn)試塊測量結(jié)果 mm

4 測量結(jié)果分析

4.1 不同測量方法測量結(jié)果對比

圖5，6為使用不同方法測量模擬葉片及標(biāo)準(zhǔn)試塊的測量誤差比較。從圖中可以看出，半高法和使用高級模式分割算法的半自動法得到的測量值比較接近。在測量小于0.6mm 的小尺寸試塊時，兩種方法均出現(xiàn)了較大的誤差。在待測尺寸不小于0.6mm時，半高法測量誤差90%在0.05mm 之內(nèi)，使用高級模式分割算法的半自動法測量誤差100%在0.05mm 之內(nèi)。使用標(biāo)準(zhǔn)模式的半自動法測量偏差值無明顯規(guī)律，且44%測量偏差值在0.05mm以上。

圖5 模擬葉片不同測量方法對比

圖6 模擬葉片不同測量方法對比

4.2 測量值校正

在實際尺寸測量檢測中為避免過大的誤差，經(jīng)常會使用經(jīng)校準(zhǔn)尺寸標(biāo)準(zhǔn)試塊，并對測量值進(jìn)行校正。最簡單的測量值校正算法是線性插值算法，將測量值記為D測量，選取校準(zhǔn)試塊中與D測量最接近的兩片，這兩片試塊實際值記為Ti、Ti＋1、測量值記為Si、Si＋1，則

選用高級模式圖像分割的半自動測量法，以校準(zhǔn)試塊尺寸及測量值（見表4）校正模擬葉片測量值（見表3），測量值與校正值與實際尺寸的誤差對比見圖7。與測量值相比，校正值的誤差基本都在0.05mm以下，避免了小尺寸測量誤差過大的情況；但是對于尺寸大于0.75mm 原本測量精度較高的點(diǎn)，不但沒有減少誤差，反倒由于引入了校準(zhǔn)試塊測量的誤差，導(dǎo)致誤差擴(kuò)大。因此，僅推薦對小于0.75mm的尺寸測量值進(jìn)行修正。

圖7 測量值與校正值誤差對比

表5 模擬葉片尺寸測量校正結(jié)果 mm

5 結(jié)論

通過對工業(yè)CT 尺寸測量半高法和半自動法測量結(jié)果的對比分析，得出：①準(zhǔn)確圖像分割算法是半自動測量法的基礎(chǔ)，只有進(jìn)行了正確的邊界提取，才能得到高精度的測量數(shù)據(jù)；②半高法和半自動測量法都能得到較為準(zhǔn)確的測量數(shù)據(jù)，但是半自動測量法的效率高于半高法；③零件在CT 掃描視場中的位置對于尺寸測量結(jié)果影響較小，因此，在射線束能量足夠且盡量避免零件間散射線影響的前提下，適當(dāng)增加單次CT 掃描的零件數(shù)量是可行的；④由于高能工業(yè)CT 系統(tǒng)對小尺寸（小于0.8 mm）測量誤差較大，可針對小尺寸測量增加數(shù)據(jù)校正。

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

［2］劉玲慧.工業(yè)CT／DR 圖像缺陷檢測算法研究［D］.重慶：重慶大學(xué)，2009.

［3］GB／T 29067－2012 無損檢測 工業(yè)計算機(jī)層析成像（CT）圖像測量方法［S］.