固體火箭發(fā)動(dòng)機(jī)界面脫粘切向CT檢測①

陳慶貴，盧洪義，齊　強(qiáng)，于光輝，張宗偉，黎國寶

(1.海軍航空工程學(xué)院，煙臺(tái)　264001;2.海軍航空工程學(xué)院青島校區(qū)，青島　266041;3.91213部隊(duì)，煙臺(tái)　264003；4.海軍駐貴陽地區(qū)軍事代表辦事處航空發(fā)動(dòng)機(jī)辦公室，平壩　561102)

?

固體火箭發(fā)動(dòng)機(jī)界面脫粘切向CT檢測①

陳慶貴1，2，盧洪義1，齊強(qiáng)1，于光輝2，張宗偉3，黎國寶4

(1.海軍航空工程學(xué)院，煙臺(tái)264001;2.海軍航空工程學(xué)院青島校區(qū)，青島266041;3.91213部隊(duì)，煙臺(tái)264003；4.海軍駐貴陽地區(qū)軍事代表辦事處航空發(fā)動(dòng)機(jī)辦公室，平壩561102)

為了實(shí)現(xiàn)固體火箭發(fā)動(dòng)機(jī)界面脫粘的快速檢測，研究了窄角扇束工業(yè)CT檢測固體火箭發(fā)動(dòng)機(jī)界面脫粘的掃描方案和圖像重建算法，提出了平移/旋轉(zhuǎn)和只旋轉(zhuǎn)掃描方案。采用窄角扇束工業(yè)CT，對某型固體火箭發(fā)動(dòng)機(jī)進(jìn)行掃描，并采用濾波反投影算法、凸集投影法和基于凸集投影的總變差最小化方法，對獲得的切向投影數(shù)據(jù)進(jìn)行重建。重建結(jié)果表明，平移/旋轉(zhuǎn)掃描方式較只旋轉(zhuǎn)掃描方式更適合該窄角扇束工業(yè)CT檢測固體火箭發(fā)動(dòng)機(jī)的界面脫粘?；谕辜队暗目傋儾钭钚』椒ㄖ亟ǖ膱D像質(zhì)量優(yōu)于濾波反投影算法和凸集投影法重建的圖像質(zhì)量，能夠滿足固體火箭發(fā)動(dòng)機(jī)脫粘缺陷的檢測要求。

窄角扇束工業(yè)CT；固體火箭發(fā)動(dòng)機(jī)；界面脫粘；切向重建；凸集投影總變差最小化

0　引言

固體火箭發(fā)動(dòng)機(jī)的粘接界面包括殼體與絕熱層、絕熱層與襯層、襯層與推進(jìn)劑之間的粘接界面，或者是殼體與襯層、襯層與推進(jìn)劑之間的粘接界面。界面脫粘是固體火箭發(fā)動(dòng)機(jī)的常見缺陷，嚴(yán)重影響固體火箭發(fā)動(dòng)機(jī)的結(jié)構(gòu)完整性和工作可靠性[1]。為了確保固體火箭發(fā)動(dòng)機(jī)的安全使用，需要對其進(jìn)行無損探傷檢測。目前，用于固體火箭發(fā)動(dòng)機(jī)脫粘檢測的方法包括超聲、射線、激光、紅外、微波和工業(yè)CT。與其他方法相比，工業(yè)CT檢測固體火箭發(fā)動(dòng)機(jī)的脫粘缺陷具有方便、直觀的特點(diǎn)[2-3]。

國內(nèi)外關(guān)于工業(yè)CT檢測固體火箭發(fā)動(dòng)機(jī)脫粘缺陷的研究相對較少。國內(nèi)方面，史源源研究了三代工業(yè)CT用于大型固體火箭發(fā)動(dòng)機(jī)的脫粘檢測方案和濾波反投影(FBP)圖像重建算法，并對所提出的方案進(jìn)行了計(jì)算機(jī)仿真驗(yàn)證[4-5]。國外方面，美國Omega公司的Gupta N K等提出了基于面檢測的切向CT檢測方案，并研制了相應(yīng)的切向CT系統(tǒng)用于檢測固體火箭發(fā)動(dòng)機(jī)的脫粘。與傳統(tǒng)CT不同，該切向CT系統(tǒng)采用的線陣探測器系統(tǒng)的布置方向與固體火箭發(fā)動(dòng)機(jī)的旋轉(zhuǎn)軸平行，可實(shí)現(xiàn)固體火箭發(fā)動(dòng)機(jī)粘接界面的軸向檢測[6-7]。

由于固體火箭發(fā)動(dòng)機(jī)切向檢測獲得的投影數(shù)據(jù)是不完備的，對其投影數(shù)據(jù)的重建屬于不完全投影重建，切向重建，又可稱作外部重建。關(guān)于切向重建算法，包括奇異值分解(SVD)法、Lambda Tomography算法和迭代法[8]。Quinto E T 研究了切向重建的奇異值分解算法[9-10]，通過一種有界外插方法補(bǔ)充緊支撐函數(shù)的空元素來實(shí)現(xiàn)切向投影數(shù)據(jù)的重建。Quinto E T 還研究了Lambda Tomography 算法[11]在切向重建方面的應(yīng)用，通過補(bǔ)充投影數(shù)據(jù)實(shí)現(xiàn)切向投影數(shù)據(jù)的Lambda 重建。熊亮等則研究了切向投影數(shù)據(jù)的迭代算法IRR法和ML-EM法，首先利用已知的投影數(shù)據(jù)估計(jì)缺失的投影數(shù)據(jù)，然后通過不斷迭代重建物體的圖像[12]。Zeng L等針對管道的切向投影數(shù)據(jù)，提出了基于子區(qū)域平滑的凸集投影總變差最小化(SA-TVM-POCS)迭代算法[13]，并進(jìn)行了計(jì)算機(jī)仿真驗(yàn)證，取得了較好的重建效果。

目前，工業(yè)CT檢測固體火箭發(fā)動(dòng)機(jī)的界面脫粘采用的是全斷面掃描成像，通過掃描并重建固體火箭發(fā)動(dòng)機(jī)的整個(gè)斷面，來判斷是否存在界面脫粘，缺點(diǎn)是檢測時(shí)間長、檢測成本昂貴。因此，本文結(jié)合窄角扇束工業(yè)CT的特點(diǎn)，提出了固體火箭發(fā)動(dòng)機(jī)界面脫粘的快速檢測方案，同時(shí)還研究了基于凸集投影的總變差最小化(POCS-TVM)算法用于窄角扇束工業(yè)CT切向重建的可行性。

1　切向檢測方案

窄角扇束工業(yè)CT檢測固體火箭發(fā)動(dòng)機(jī)時(shí)，發(fā)動(dòng)機(jī)呈“臥式”狀態(tài)。對固體火箭發(fā)動(dòng)機(jī)的粘接界面進(jìn)行CT檢測時(shí)，可供選擇的檢測方案有如下兩種：一種方案是采用平移/旋轉(zhuǎn)(TR)掃描方式，如圖1所示。掃描開始前，使射線源最右側(cè)射線與發(fā)動(dòng)機(jī)的殼體相切；然后，射線源和探測器組以某一固定步長平移掃描至射線源最左側(cè)射線掃描完襯層部分，完成一次切向投影數(shù)據(jù)的采集。接著，固體火箭發(fā)動(dòng)機(jī)繞其旋轉(zhuǎn)中心轉(zhuǎn)過某一固定角度，射線源和探測器組回到起始掃描位置，以相同步長再次對發(fā)動(dòng)機(jī)的粘接界面進(jìn)行掃描，直至完成360°范圍內(nèi)投影數(shù)據(jù)的采集為止。另一種方案是采用只旋轉(zhuǎn)(RO)掃描方式，如圖2所示。掃描時(shí)，射線源和探測器組固定不動(dòng)，使射線源最左測射線與發(fā)動(dòng)機(jī)的外層殼體相切，且令發(fā)動(dòng)機(jī)的粘接界面位于射線源和探測器組的掃描范圍內(nèi)。然后，固體火箭發(fā)動(dòng)機(jī)繞其旋轉(zhuǎn)中心以某一固定角度旋轉(zhuǎn)，發(fā)動(dòng)機(jī)每旋轉(zhuǎn)1次完成1次切向投影數(shù)據(jù)的采集，直至完成360°范圍內(nèi)投影數(shù)據(jù)的采集為止。

圖1　平移/旋轉(zhuǎn)切向檢測示意圖

圖2　只旋轉(zhuǎn)切向檢測示意圖

2　切向投影數(shù)據(jù)的獲取

當(dāng)采用平移/旋轉(zhuǎn)掃描方案時(shí)，需要確定射線源和探測器組的起始掃描位置和終止掃描位置。這里計(jì)算窄角扇束中心射線S0D0(射線源與探測器組中間探測器之間的連線)的起始掃描位置和終止掃描位置。記S0D0與最右側(cè)射線間的夾角為γ0，固體火箭發(fā)動(dòng)機(jī)殼體的外層半徑為R，襯層半徑為r，射線源到發(fā)動(dòng)機(jī)旋轉(zhuǎn)中心的垂直距離為D， 如圖3所示。

圖3　平移/旋轉(zhuǎn)掃描方式下起始和終止掃描位置示意圖

開始掃描時(shí)，最右側(cè)射線與發(fā)動(dòng)機(jī)的外層殼體相切。所以，窄角扇束中心射線S0D0距發(fā)動(dòng)機(jī)旋轉(zhuǎn)中心的距離OC為

OC=OB+BC

=OA/cosγ0+S0Ctanγ0

(1)

結(jié)束掃描時(shí)，最左側(cè)射線掃描完固體火箭發(fā)動(dòng)機(jī)的襯層部分，此時(shí)窄角扇束中心射線S0D0距發(fā)動(dòng)機(jī)旋轉(zhuǎn)中心的距離OE為

OE=OB-BE

(2)

平移/旋轉(zhuǎn)掃描方式下，固體火箭發(fā)動(dòng)機(jī)每旋轉(zhuǎn)1次，射線源和探測器組的起始掃描位置和終止掃描位置都相同。發(fā)動(dòng)機(jī)繞其旋轉(zhuǎn)中心旋轉(zhuǎn)1周，完成360°的掃描，便可獲得所需的切向投影數(shù)據(jù)。

圖4　只旋轉(zhuǎn)掃描方式下起始和終止掃描位置示意圖

對于只旋轉(zhuǎn)掃描方式，調(diào)節(jié)射線源至發(fā)動(dòng)機(jī)旋轉(zhuǎn)中心的垂直距離 ，使得窄角扇束射線張成的扇束角能夠覆蓋發(fā)動(dòng)機(jī)的粘接界面，且使最左側(cè)射線與發(fā)動(dòng)機(jī)的外層殼體相切，如圖4所示。

此時(shí)，窄角扇束中心射線S0D0距發(fā)動(dòng)機(jī)旋轉(zhuǎn)中心的距離OB為

OB=OA-AB=R-Dtanγ1

(3)

式中γ1表示中心射線S0D0與最左側(cè)射線的夾角。

發(fā)動(dòng)機(jī)繞其旋轉(zhuǎn)中心旋轉(zhuǎn)1周，完成360°的掃描，便可獲得所需的切向投影數(shù)據(jù)。

3　切向重建算法

由于獲得的切向投影數(shù)據(jù)是不完備的，這里選用基于凸集投影的總變差最小化算法進(jìn)行圖像重建。設(shè)重建的CT圖像為Ⅰ×Ⅰ維的X0，將其排列為一維向量X。POCS-TVM算法的具體實(shí)現(xiàn)步驟[14-15]如下：

步驟一凸集投影過程。

(2)對圖像進(jìn)行迭代運(yùn)算：

(4)

式中λ為松弛因子；j為射線編號(hào)，j=1,2,…,L(L為投影射線總數(shù))；pj為與射線j對應(yīng)的投影值；Wj為與射線j對應(yīng)的權(quán)因子矩陣。

步驟二總變差最小化過程。

(2)計(jì)算增量因子dm：

(5)

(3)計(jì)算全變分梯度大小及梯度方向：

(6)

(7)

式中n=1,2,…,N；N表示總變差最小化過程的最大迭代次數(shù)。

(4)沿全變分梯度下降的方向，對圖像進(jìn)行迭代修正：

(8)

式中α為調(diào)節(jié)因子，反復(fù)計(jì)算全變分梯度，并迭代修正圖像，直至n=N為止。

(9)

4　實(shí)驗(yàn)驗(yàn)證

由于受機(jī)械掃描系統(tǒng)的限制，實(shí)驗(yàn)采用的窄角扇束工業(yè)CT系統(tǒng)無法實(shí)現(xiàn)固體火箭發(fā)動(dòng)機(jī)粘接界面的切向檢測，采用如下方法驗(yàn)證所提出的固體火箭發(fā)動(dòng)機(jī)界面脫粘切向CT檢測方案和基于凸集投影的總變差最小化算法用于切向重建的有效性。采用450 kV射線源的窄角扇束工業(yè)CT，對某型固體火箭發(fā)動(dòng)機(jī)進(jìn)行全斷面檢測，從測得的全局投影數(shù)據(jù)中，按平移/旋轉(zhuǎn)和只旋轉(zhuǎn)掃描方式選取相對應(yīng)的切向投影數(shù)據(jù)，相關(guān)實(shí)驗(yàn)參數(shù)如表1所示。

表1　CT檢測實(shí)驗(yàn)參數(shù)設(shè)置

CT檢測獲得的全局投影數(shù)據(jù)和圖像重建結(jié)果如圖5所示，相應(yīng)的切向投影數(shù)據(jù)如圖6所示。為了比較圖像重建質(zhì)量，分別采用了濾波反投影(FBP)算法[16-17]、POCS算法和POCS-TVM算法，對切向投影數(shù)據(jù)進(jìn)行重建。其中，λ=0.25，α=0.002。切向重建獲得的圖像如圖7和圖8所示。

(a)CT檢測數(shù)據(jù)　　　　　　(b)重建圖像

(a)平移/旋轉(zhuǎn)掃描　　　　(b)只旋轉(zhuǎn)掃描

(a)FBP　　　　　　(b)POCS　M=20　　　　(c)POCS　M=50　　　　(d)　　　POCS-TVM

(e)　　　POCS-TVM 　　　　　(f)　　　POCS-TVM　　　　　(g)　　　POCS-TVM　　　　　(h)　　　POCS-TVM

(a)FBP　　　　　(b)POCS　M=20　　　　(c)POCS　M=50　　　　(d)　　　POCS-TVM

(e)　　　POCS-TVM 　　　　　(f)　　　POCS-TVM　　　　　(g)　　　POCS-TVM　　　　　(h)　　　POCS-TVM

由圖7和圖8的重建結(jié)果可看出，F(xiàn)BP算法的圖像重建質(zhì)量最差；POCS算法的圖像重建質(zhì)量次之，POCS-TVM算法的圖像重建質(zhì)量最好。TR掃描方式下的圖像重建質(zhì)量較RO掃描方式下的圖像重建質(zhì)量好。RO掃描模式下的重建圖像偽影嚴(yán)重，影響對固體火箭發(fā)動(dòng)機(jī)脫粘缺陷的判斷。TR掃描方式下FBP算法和POCS算法的重建圖像也存在偽影，而POCS-TVM算法在取合適的凸集投影迭代次數(shù)M和總變差最小化迭代次數(shù)N時(shí)重建圖像質(zhì)量較好，在重建圖像中，能清晰觀察到脫粘缺陷。

5　結(jié)論

(1)對于本實(shí)驗(yàn)中的窄角扇束工業(yè)CT，由于發(fā)動(dòng)機(jī)每次旋轉(zhuǎn)的角度為10°，使得采用只旋轉(zhuǎn)(RO)掃描方式獲得的投影數(shù)據(jù)量較少，重建圖像的質(zhì)量較差；而平移/旋轉(zhuǎn)(TR)掃描方式獲得的投影數(shù)據(jù)量較多，采用同一重建算法重建的圖像質(zhì)量較RO方式下的圖像重建質(zhì)量好。因此，TR掃描方式較RO掃描方式更適合該窄角扇束工業(yè)CT系統(tǒng)檢測固體火箭發(fā)動(dòng)機(jī)的界面脫粘。

(2)凸集投影-總變差最小化算法對固體火箭發(fā)動(dòng)機(jī)粘接界面的投影數(shù)據(jù)進(jìn)行切向重建是可行的，在松弛因子λ=0.25、調(diào)節(jié)因子α=0.002的情況下，選擇凸集投影的迭代次數(shù)為M=100、總變差最小化的迭代次數(shù)N=5時(shí)，能獲得最佳的圖像重建質(zhì)量，能夠滿足固體火箭發(fā)動(dòng)機(jī)界面脫粘的檢測要求。

[1]盧洪義,楊興根.程衛(wèi)平.小張角扇形射線束CT檢測固體發(fā)動(dòng)機(jī)局部缺陷[J].推進(jìn)技術(shù),2003,24(3): 251-253.

[2]盧洪義,朱敏,于光輝,等.固體火箭發(fā)動(dòng)機(jī)三維可視化故障診斷技術(shù)[M].北京：國防工業(yè)出版社，2014.

[3]艾春安,桑元,徐志高,等.基于X射線膠片的固體火箭發(fā)動(dòng)機(jī)斷面圖像重建[J].無損檢測，2011,33(3): 1-4.

[4]史源源,王召巴,金永.大型固體火箭發(fā)動(dòng)機(jī)脫粘CT檢測技術(shù)[J].彈箭與制導(dǎo)學(xué)報(bào)，2006,26(1): 928-929.

[5]史源源.切向CT重建技術(shù)研究[D].太原：中北大學(xué)，2005.

[6] Gupta N K,Hughes S H,Isaacson B.Tangential C T,a computed tomography method developed for industrial inspection[C]//16th World Conference on Nondestructive Testing.Canada,2004.

[7]Isaacson B,Gupta N K.Tangential CT as an NDT inspection tool[C]//ASNT Fall Conference and Quality Testing Show.USA:2000.

[8]金鑫,陳志強(qiáng),李亮,等.CT外部重建算法發(fā)展綜述[C]//第十四屆全國核電子學(xué)與核探測技術(shù)學(xué)術(shù)年會(huì)文集,2008:756-760.

[9]Quinto E T.Singular value decompositions and inversion methods for the exterior Radon transform and a spherical transform[J].Journal of Mathematical Analysis and Applications,1983,95: 437-448.

[10]Quinto E T.Tomographic reconstructions from incomplete data-numerical inversion of the exterior Radon transform[J].Inverse Problems,1988,(4): 867-876.

[11]Quinto E T.Local algorithms in exterior tomography[J].Journal of Computational and Applied Mathematics,2007,199:141-148.

[12]熊亮,張麗,陳志強(qiáng),等.切向CT重建算法的研究[J].CT理論與應(yīng)用研究，2007,16(1):10-18.

[13]Zeng L,Liu B D,Liu L H,et al.A new iterative reconstruction algorithm for 2D exterior fan-beam CT[J].Journal of X-Ray Science and Technology,2010,18:267-277.

[14]李超.基于壓縮感知理論的CT圖像重建算法研究[D].西安:西安電子科技大學(xué),2013.

[15]徐志高,艾春安,馬靜.固體火箭發(fā)動(dòng)機(jī)CT圖像代數(shù)重建算法研究[J].固體火箭技術(shù)，2013,36(5)：697-700.

[16]Dennerlein F,Noo F,Hornegger J,et al.Fan-beam filtered backprojection reconstruction without backprojection weight[J].Physics in Medicine and Biology,2007,52: 3227-3240.

[17]陳慶貴,劉趙云,盧洪義,等.基于窄角扇束掃描的固體火箭發(fā)動(dòng)機(jī)CT圖像重建[J].固體火箭技術(shù),2014,37(6): 877-880.

(編輯：薛永利)

Tangential CT inspection of interface debonding of SRM

CHEN Qing-gui2,LU Hong-yi1,QI Qiang1,YU Guang-hui2,ZHANG Zong-wei3,LI Guo-bao4

(1.Naval Aeronautical and Astronautical University,Yantai264001,China;2.Qingdao Branch of Naval Aeronautical and Astronautical University,Qingdao266041,China;3.Troop 91213,Yantai264003,China;4.Aero-engine Office of Military Representatives Office of Navy in Guiyang,Pingba561102,China)

To inspect the interface debonding of solid rocket motor(SRM),rapid scanning schemes and image reconstruction algorithms for the inspection of SRM interface debonding based on narrow fan beam industrial computed tomography (ICT) were investigated.Translation/rotation and rotation only scanning schemes were proposed.One SRM was scanned by narrow fan beam ICT.Filtered backprojection (FBP) algorithm,projection onto convex sets (POCS) algorithm and projection onto convex sets-total variation minimization (POCS-TVM) algorithm were applied to reconstruct images by using obtained tangential projections.Results show that the translation/rotation scanning scheme is more effective than the rotation only scanning scheme for the narrow fan beam ICT to inspect interface debonding of SRM.The reconstructed image quality of POCS-TVM algorithm is better than those of FBP and POCS algorithms,which can meet the need of inspection of SRM's interface debonding.

narrow fan beam ICT;solid rocket motor;interface debonding;tangential reconstruction;projection onto convex sets-total variation minimization

2014-10-20；

2015-05-10。

國家自然科學(xué)基金項(xiàng)目(51005242)。

陳慶貴(1987—)，男，博士生，研究方向?yàn)榘l(fā)動(dòng)機(jī)測試?yán)碚撆c技術(shù)。E-mail：cqgccc12345@163.com

盧洪義，男，教授/博導(dǎo)。E-mail：ytluhongyi@sina.com

V435

A

1006-2793(2016)03-0347-06

10.7673/j.issn.1006-2793.2016.03.010