基于工業(yè)計算機(jī)斷層成像的裝藥底隙無損檢測方法研究

呂寧，徐更光

（北京理工大學(xué)爆炸科學(xué)與技術(shù)國家重點實驗室，北京100081）

基于工業(yè)計算機(jī)斷層成像的裝藥底隙無損檢測方法研究

呂寧，徐更光

（北京理工大學(xué)爆炸科學(xué)與技術(shù)國家重點實驗室，北京100081）

在彈藥裝藥底隙檢測中，運用密度對比法結(jié)合體積效應(yīng)，將底隙測量轉(zhuǎn)化為彈底密度變化的測量，建立彈底掃描斷層中心部位計算機(jī)斷層成像（CT）值變化與底隙值之間的數(shù)學(xué)模型，解決底隙測量理論與實際操作中的難題。預(yù)置底隙測量實驗結(jié)果表明：該模型可完成對0.20～0.50 mm底隙的測量，最大相對誤差小于3.5%，可實現(xiàn)對裝藥底隙的定性判斷與定量檢測。

兵器科學(xué)與技術(shù)；工業(yè)計算機(jī)斷層成像；密度對比；底隙；定量檢測

0　引言

在彈藥生產(chǎn)過程中，由于生產(chǎn)因素影響、凝固裝藥與鋼制殼體膨脹系數(shù)不同，使得藥柱與彈底之間存在間隙，也稱為底隙。底隙尺寸一般在0.10～2.00 mm之間，雖然尺寸很小，但嚴(yán)重影響發(fā)射安全性：發(fā)射過程中，藥柱底部與殼體之間的空氣受到高沖擊、高過載壓縮，溫度升高，加熱周圍裝藥形成熱點［1］，同時伴隨藥柱的軸向移動，產(chǎn)生很大的摩擦力和沖擊力，使得裝藥的分解速率提高，增加了膛炸的可能性。大量研究工作表明［2-4］，當(dāng)?shù)紫冻叽缧∮?.40 mm時，發(fā)射過程中彈底主應(yīng)力與無底隙時基本一致。當(dāng)?shù)紫冻叽绱笥?.40 mm時彈底主應(yīng)力開始變化，超過0.50 mm后，最大主應(yīng)力顯著升高。因此，美軍標(biāo)規(guī)定底隙的最大允許尺寸為0.381 mm［5］.

現(xiàn)有的底隙檢測方法為抽樣檢測法，根據(jù)生產(chǎn)量的大小，人為確定抽樣比例，剖解樣品，用樣品的質(zhì)量表征同批次產(chǎn)品的質(zhì)量。在生產(chǎn)的初期工藝缺陷出現(xiàn)概率較高，采用抽樣檢測很可能接受有缺陷的產(chǎn)品，給后續(xù)加工、使用帶來風(fēng)險［6］。因此，在不損傷產(chǎn)品的前提下，高精度、高效率地檢測出裝藥中是否存在底隙、確定底隙尺寸，并判定產(chǎn)品質(zhì)量的檢測方法研究就顯得意義重大。工業(yè)計算機(jī)斷層成像（CT）采用非傳統(tǒng)的成像方法，通過非接觸、非破壞的方式得到炮彈內(nèi)部裝藥檢測斷層各點的射線衰減系數(shù)，并建立相應(yīng)的斷層圖像，以此獲得裝藥疵病的檢測數(shù)據(jù)。但因CT系統(tǒng)本身的噪聲和偽影的存在，底隙尺寸小且與斷層面平行，無法直接準(zhǔn)確從斷層圖像觀測到軸向尺寸數(shù)據(jù)［3］。針對這些問題，本文擬通過采用密度檢測中常用的對比法，將底隙的尺寸檢測問題轉(zhuǎn)化為密度的對比，從側(cè)面反映底隙的存在并測量尺寸。

1　檢測原理

在實際生產(chǎn)中，無法在實彈中預(yù)置底隙，故采用如圖1所示底隙模擬彈作為檢測對象。

圖1　底隙模擬彈示意圖Fig.1 Model of projectile with base separation

該模擬彈為某型號實彈殼體靠近彈底部位徑向剖開，加工完成匹配的上部殼體后重新注藥，上部殼體與藥柱一體，方便取出。采用預(yù)置尺寸環(huán)狀墊片構(gòu)造底隙，這樣既保證了底隙檢測的真實性又方便預(yù)置不同的底隙尺寸。如在上、下部之間加裝厚度為k的環(huán)狀墊片后，藥柱底部升高同樣距離，底隙尺寸即為k.

在工業(yè)CT掃描中，準(zhǔn)直后的射線并非理想的無限薄，其厚度由準(zhǔn)直器和探測器決定，所以在裝藥質(zhì)量檢測中重建的斷層圖像數(shù)據(jù)為一定厚度裝藥密度的平均結(jié)果，CT斷層掃描的這種特性被稱為“體積效應(yīng)”［7-8］。

合理的運用體積效應(yīng)，選擇斷層厚度，對彈底可能存在底隙的位置由下至上進(jìn)行連續(xù)多層掃描，測量斷層中心區(qū)域平均密度。因鋼質(zhì)殼體、底隙中的空氣、藥柱的密度各不相同，斷層每移動一個步長，3種物質(zhì)在斷層中的比例也會隨之改變。從而影響斷層表征的平均密度，通過測量密度變化，建立相應(yīng)的CT值變化曲線，可實現(xiàn)底隙有無的定性判斷與尺寸大小的定量檢測。

2　數(shù)學(xué)模型

2.1CT值校正

底隙因位置在彈底鋼殼和裝藥高低密度交界處，檢測有其特殊性。對中大口徑彈藥，射線在該部位透射距離在100 mm以上，且鋼質(zhì)殼體厚度大、密度高，與裝藥密度差較大，二者對射線作用的規(guī)律不相同，射線在透射殼體過程中產(chǎn)生硬化效應(yīng)和散射，導(dǎo)致該區(qū)域藥柱對應(yīng)CT值高于真實值，與密度對應(yīng)關(guān)系產(chǎn)生變化［9］，將影響基于密度對比的裝藥底隙測量精度，故應(yīng)分析殼體對裝藥CT值影響規(guī)律并予以消除。

取圖1所示模擬彈的空殼體進(jìn)行多次掃描，結(jié)果取平均值，選取某一位置斷層數(shù)據(jù)，發(fā)現(xiàn)殼體內(nèi)部CT值并不為0，呈現(xiàn)規(guī)律如圖2所示。

圖2　空彈殼內(nèi)CT分布曲線圖Fig.2 Distribution curve of CT in empty shell

對圖2所示數(shù)據(jù)擬合，得到空殼體該位置斷層CT值分布CT′與殼體內(nèi)徑r關(guān)系式為

式中：乘積因子α=33.64；指數(shù)因子β=10.53.

因彈體形狀特殊，在彈底部位藥室內(nèi)徑緩慢變化，對CT值的影響規(guī)律也不盡相同。對壁厚均為20 mm，內(nèi)徑20～60 mm的鋼制殼體掃描，測得空殼體內(nèi)徑向CT值分布如圖3所示。

圖3　殼體內(nèi)徑對CT值影響曲線圖Fig.3 Curves of CT values affected by internal diameter of shell

由圖3可見，壁厚不變的情況下，內(nèi)徑的改變并未對殼體內(nèi)中心區(qū)域CT值造成較大影響，內(nèi)徑40～60 mm時，中心區(qū)域半徑5 mm內(nèi)CT值基本一致。因此，在中大口徑彈藥底隙檢測時，可忽略彈底殼體內(nèi)徑變化對CT值的影響。

對內(nèi)徑均為60 mm，壁厚15～30 mm的鋼質(zhì)殼體進(jìn)行測量，測得空殼體內(nèi)CT值分布如圖4所示。

圖4　殼體厚度對CT值影響曲線圖Fig.4 Curves of CT values affected by thickness of shell

為簡化模型，提高運算速度，采用與（1）式相同的形式對圖4數(shù)據(jù)進(jìn)行分析，結(jié)果如表1所示。

綜上所述，中大口徑裝藥底隙檢測中，對CT值造成影響的主要因素是殼體厚度變化，CT值校正公式如下：

式中：CT1、CT0為校正前、后CT值；系數(shù)α、β隨所在斷層殼體壁厚變化而改變。

為驗證上述校正法的可行性，將該方法應(yīng)用于殼體底部位置多個斷層中心對應(yīng)密度值測量，因CT值與密度呈正比，故可通過密度測量的精度驗證校正（2）式，結(jié)果如表2所示。

表1　殼體厚度對CT值影響曲線擬合結(jié)果Tab.1 Curve fitting coefficients of CT values affected by thickness of shell

表2　斷層中心密度測量結(jié)果Tab.2 Measured results of central area density in cross-section

由表2結(jié)果可知，運用上述CT值校正方法后，彈底部位密度測量的最大絕對誤差為0.008 g/cm3，最大相對誤差為0.477%，說明CT值校正（2）式準(zhǔn)確可靠，能基本消除殼體厚度變化對CT值的影響，為基于密度對比的高精度底隙測量打下基礎(chǔ)。

2.2數(shù)學(xué)模型建立

2.2.1無底隙標(biāo)準(zhǔn)曲線

在試驗中，將無底隙裝藥由下至上，斷層中心區(qū)域的CT值掃描曲線定義為標(biāo)準(zhǔn)曲線，如圖5所示。并以其為基準(zhǔn)，分析對比不同預(yù)置底隙裝藥的掃描曲線。

圖中x軸為掃描斷層位置沿藥柱軸向位移，y軸為不同斷層中心區(qū)域?qū)?yīng)的CT值。從中可以看出，隨著x的增大，CT值逐漸變小，并最終趨于平穩(wěn)。為更好分析說明，將曲線分為3個部分：A為殼體區(qū)，B為底隙區(qū)，C為裝藥區(qū)。在彈底殼體區(qū)，射線透射鋼殼路徑較長，衰減較多，信噪比低，CT值存在一定的波動，但因該區(qū)域斷層主要為高密度的鋼制殼體，底隙及裝藥對CT值的影響較小，故該區(qū)域曲線呈相對平緩的下降趨勢。在底隙區(qū)，殼體的影響逐漸減小，斷層中底隙與裝藥所占的比例逐漸擴(kuò)大，平均密度也隨之下降，隨著x的增大，CT值呈明顯的減小趨勢。底隙尺寸不同，相同位置斷層內(nèi)的平均密度也不盡相同，則不同預(yù)置底隙裝藥掃描曲線在該區(qū)域應(yīng)呈現(xiàn)相同趨勢、不同程度的變化，由此可實現(xiàn)底隙有無的定性判斷與尺寸大小的定量檢測。在裝藥區(qū)，掃描斷層中殼體影響消失，低密度裝藥成為主要成分，隨著x的增大，裝藥比例繼續(xù)增大，并最終充滿整個斷層，斷層中心CT值也趨于穩(wěn)定。

圖5　無底隙標(biāo)準(zhǔn)曲線Fig.5 Standard curve without base separation

2.2.2預(yù)置底隙曲線

通過加裝不同厚度墊片，預(yù)置底隙尺寸為k1和k2的模擬彈（k1＜k2），掃描結(jié)果與標(biāo)準(zhǔn)曲線對比如圖6所示。

圖6　預(yù)置k底隙曲線圖Fig.6 Curves of base separation of size k

從圖6可以看出：不同底隙的模擬彈掃描曲線，隨著掃描斷層位移的變化，曲線變化趨勢相同。在殼體區(qū)與裝藥區(qū)，所有曲線基本一致，但在底隙區(qū)，因底隙尺寸不同，斷層中心區(qū)域平均密度改變，曲線逐漸分開，呈變化趨勢一致，但變化程度不同的發(fā)展，且底隙越大，曲線下降趨勢越明顯，與標(biāo)準(zhǔn)曲線之間的差別就越大。隨著斷層位移的改變，曲線發(fā)展至近裝藥區(qū)，又呈現(xiàn)彼此相交的趨勢。則在曲線分開發(fā)展、變化差異最大的區(qū)域［x1，x2］，所圍成的積分面積就越大，這與之前的分析預(yù)測相吻合?？山⒌紫冻叽鏺與積分面積Sk的函數(shù)關(guān)系，通過測量反映密度變化的CT值曲線，計算Sk值，從而求得底隙尺寸k的大小。

2.2.3函數(shù)關(guān)系建立

預(yù)置無底隙模擬彈，選取曲線在底隙區(qū)內(nèi)任一點x，可知該點對應(yīng)斷層應(yīng)由鋼殼和裝藥組成，設(shè)密度分別為ρ1、ρ2，在斷層所占比例分別為αx、βx（0≤αx≤1，αx+βx=1）。則該斷層平均密度為

預(yù)置底隙尺寸為k的模擬彈，選取與前述無底隙模擬彈相同的位置x掃描，則所得斷層除了鋼和裝藥外，應(yīng)有部分空氣。假設(shè)鋼所占比例不變，仍為αx，空氣所占比例設(shè)為γx，若忽略空氣密度，則該斷層平均密度為

對比可知，因底隙存在造成的該點密度差為

則對應(yīng)底隙區(qū)間［x1，x2］上各個點的密度差累加為

式中：A、B為待定系數(shù)，受所測彈種影響，隨殼體材料，裝藥性質(zhì)不同而改變。

運用圖1所示模擬彈，預(yù)置0.10～0.50 mm底隙，多次測量平均結(jié)果如圖7所示。

觀察可知，除0.10 mm處數(shù)據(jù)誤差較大外，其他各點線性關(guān)系良好，因底隙檢測標(biāo)準(zhǔn)為0.381 mm，遠(yuǎn)大于0.10 mm，故可對0.10 mm底隙定量測量不做要求，擬合其他各點得線性關(guān)系式如下：

當(dāng)Sk=0時，由（8）式可得k=-0.002 44 mm，出現(xiàn)底隙為負(fù)數(shù)情況。結(jié)合0.10 mm底隙測量結(jié)果，說明（8）式不能滿足小于等于0.10 mm底隙的定量測量，只能結(jié)合曲線對比圖，實現(xiàn)底隙有無的定性判斷。

綜上所述，針對該型號炮彈，運用密度對比法檢測底隙的模型已經(jīng)建立成功。在實際檢測過程中，分別對無底隙彈體和預(yù)置底隙模擬彈彈底部位由下至上、連續(xù)多層掃描，得到如前所述標(biāo)準(zhǔn)曲線與底隙曲線，求得積分面積Sk值，后利用（8）式即可求得底隙尺寸k.整個過程可通過編程實現(xiàn)自動化操作。

圖7　模型彈預(yù)置底隙測量結(jié)果Fig.7 Measured results of base separations of projectile

3　檢測結(jié)果與分析

為驗證密度對比法測量底隙的可靠性，預(yù)置多組0.20～0.50 mm的底隙模擬彈，運用前述模型進(jìn)行測量，得到結(jié)果如圖8所示。

圖8　底隙測量結(jié)果Fig.8 Measured results of base separation

從圖8可以看出，0.20～0.50 mm間的底隙測量值與真實值符合較好，最大絕對誤差為0.01 mm，最大相對誤差為3.49%.這說明檢測模型可較好完成0.20～0.50 mm底隙的測量。

為進(jìn)一步驗證模型的準(zhǔn)確性和可靠性，對0.35～0.40 mm范圍內(nèi)的底隙進(jìn)行了多次測量，得到結(jié)果如表3所示。

表3　0.35～0.40 mm底隙測量結(jié)果Tab.3 Measured results of 0.35～0.40 mm base separations

從表3可以看出，模型在檢測0.35～0.40 mm底隙時，保持較高準(zhǔn)確性，最大絕對誤差為0.01 mm，最大相對誤差為2.86%，能夠穩(wěn)定可靠地實現(xiàn)底隙有無的定性判斷與尺寸大小的定量檢測。

4　結(jié)論

本文將密度檢測中常用的對比法應(yīng)用于底隙檢測，提出底隙尺寸測量的對比模型，將底隙檢測的問題轉(zhuǎn)化為斷層平均密度的測量，解決了理論與實際操作中的難題?？梢钥闯觯?/p>

1）模型可實現(xiàn)0.20～0.50 mm的小尺寸底隙測量，最大絕對誤差0.01 mm，最大相對誤差3.49%，能夠滿足實際生產(chǎn)中的要求，實現(xiàn)裝藥底隙的定性判斷與定量檢測。

2）底隙尺寸小于等于0.10 mm時，該模型測量誤差相對較大，僅能實現(xiàn)底隙有無的定性判斷。但因0.10 mm小于國軍標(biāo)規(guī)定的0.381 mm的最大底隙尺寸，故不影響該檢測模型推廣使用。

3）檢測方法簡單、重復(fù)性好，易于操作。測量結(jié)果直觀明了，避免了人工判斷的主觀誤差。

4）CT值校正（2）式中系數(shù)α、β隨所在斷層殼體壁厚變化而改變。線性關(guān)系（7）式中系數(shù)A、B隨殼體材料、裝藥性質(zhì)不同而改變。因此在不同型號彈種推廣應(yīng)用中，需重新標(biāo)定。

5）該檢測方法，目前在文中所述模型彈檢測中得到驗證并保持較高精度。但在實彈應(yīng)用中，仍需解決不同批次炮彈，殼體材料和結(jié)構(gòu)加工誤差對CT值的影響問題。這也是本文后續(xù)研究、完善基于密度對比的底隙檢測方法的主要方向。

（

）

［1］Starkenberg J.Ignition of solid high explosive by the rapid compression of an adjacent gas layer［C］∥Proceedings of 7th International Symposium on Detonation.White Oak，MD：Naval SurfaceWeapons Center，1981：3-16.

［2］李文彬，王曉鳴，趙國志，等.裝藥底隙對彈底應(yīng)力及發(fā)射安全性影響研究［J］.彈道學(xué)報，2001，13（3）：64-67. LI Wen-bin，WANG Xiao-ming，ZHAO Guo-zhi，et al.The research of the effect of base gap on the stress of explosives and the lunching safety［J］.Journal of Ballistics，2001，13（3）：64-67.（in Chinese）

［3］盧艷平，王玨，喻洪麟.基于工業(yè)CT成像的彈藥小底隙測量方法研究［J］.兵工學(xué)報，2009，30（12）：1714-1717. LU Yan-ping，WANG Jue，YU Hong-lin.Study on measurement method for small bottom gap of ammunition by industrial CT imaging［J］.Acta Armamentarii，2009，30（12）：1714-1717.（in Chinese）

［4］周培毅，徐更光，王廷增.炸藥裝藥在后座沖擊條件下的點火模型研究［J］.火炸藥學(xué)報，2000，1（1）：1-5. ZHOU Pei-yi，XU Geng-guang，WANG Ting-zeng.Ignition models of explosive charge subjected to setback impact［J］.Chinese Journal of Explosives&Propellants，2000，1（1）：1-5.（in Chinese）

［5］United States Department of Defense.MIL-P-63252A Projectile， 155 mm，HE，M795 loading，assembling and packing；military specification［S］.US：United States Department of Defense，1995.

［6］盧國勝，施軍，王凱.中大口徑榴彈采用ICT檢驗注裝藥底隙方案研究［J］.彈箭與制導(dǎo)學(xué)報，2010，30（1）：121-123. LU Guo-sheng，SHI Jun，WANG Kai.Research of checking base separation to be injected of loading explosive charge in the large and middle calibers HE shells through an ICT plan［J］.Journal of Projectiles，Rockets，Missiles and Guidance，2010，30（1）：121-123.（in Chinese）

［7］Wells K，Chiverton J，Partridge M，et al.Quantifying the partial volume effect in PET using Benford's law［J］.IEEE Transactions on Nuclear Science，2007，54（5）：1616-1625.

［8］Dehmeshki J，Ye X J，Amin H，et al.Volumetric quantification of atherosclerotic plaque in CT considering partial volume effect［J］. IEEE Transactions on Medical Imaging，2007，26（3）：273-282.

［9］Kropas C V，Moran T J，Yancey R N.Effect of composition on density measurement by X-ray computed tomography［J］.Materials Evaluation，1991，49（4）：487-489.

Inspection of Base separation of Charge by Industrial CT Imaging

LYU Ning，XU Geng-guang
（State Key Laboratory of Explosion Science and Technology，Beijing Institute of Technology，Beijing 100081，China）

A volume effect of industrial computed tomography（CT）imaging and a contrast method for density inspection are used to build an algebraic model for base separation inspection，which improves the base separation inspection method in theoretical and practical ways.The size information is qualitatively and quantitatively analyzed with high precision by the fluctuation of inspected CT values.The experimental results of base separation inspection show that the density contrast model is suitable for the inspection of 0.20～0.50 mm base separation，and the relative error is less than 3.5%.

ordnance science and technology；industrial computed tomography；density contrast；base separation；quantitative inspection

TJ410.6

A

1000-1093（2015）01-0157-06

10.3969/j.issn.1000-1093.2015.01.023

2014-02-10

呂寧（1988—），男，博士研究生。E-mail：rex17@163.com；徐更光（1932—2015），男，中國工程院院士，博士生導(dǎo)師。