复杂构型PBX截面密度分布CT测试方法

杨雪海, 张伟斌, 杨仍才, 王依全

(中国工程物理研究院化工材料研究所, 四川 绵阳 621999)

1 引 言

X射线CT(X-ray computed tomography,X-CT)是无损测试材料密度分布或均匀性的重要手段[1],通过测试材料的衰减系数或CT灰度值(以下简称CT值),建立衰减系数或CT值与密度的数学关系,获得被测材料的体密度或局部密度,已用于药粒[2]、粉体[3]、木材[4]、C/C复合材料[5]、聚丙烯酸酯泡沫[6]、石墨[7]等的内部密度均匀性研究测试。

高聚物粘结炸药(Polymer Bonded Explosive,PBX)是军事应用的重要材料,其密度分布或均匀性可反映内部组成和细观结构的差异,对研究PBX的起爆及传爆特性、力学性能、安全性能具有重要意义[8-11]。PBX密度分布或均匀性测试研究主要采用CT技术[12-14],本课题组已采用测试CT值相对差的方法间接表征炸药晶体[15]、PBX[16]密度的均匀性。由于PBX试件构型复杂,受工业X-CT射线硬化的影响,复杂构型材料截面形状和面积的变化会导致CT值变化,而材料内部密度的差异也会引起CT值改变,所以难于判断CT值改变的真正原因。目前,CT测试的标准密度件基本采用圆柱形,与复杂构型待测件的形状、尺寸存在差异,如果直接利用圆柱形标准密度件的CT值,容易给复杂构型待测件的局部密度计算带来较大误差,因此复杂构型材料的局部密度测试成为难题,目前尚无CT测试方法报道。

鉴于上述问题,本实验采用双标准同步CT扫描实验模型,研究了复杂构型PBX与标准密度件的截面CT值分布规律,讨论了X射线能量波动、不同截面形状面积变化对CT值测试的影响,通过标准件的密度计算PBX的截面密度,并利用PBX体密度对计算结果进行了线性修正。

2 实验原理与测试模型

X射线穿透材料前后的光强服从Lambert-Beer′s定理[17]:

(1)

式中,I0为入射X射线强度,I为出射X射线强度,Δs-1;x为X射线穿过的厚度,cm;μ为材料的线衰减系数,cm-1;μ与X射线能量、材料化学组成和密度相关。

当X射线穿过材料某截面,根据该截面在不同方向的投影数据(即I0/I的变化)可以重建获得该截面的线衰减系数μ分布(CT图像),CT图像以灰度值方式显示在屏幕上,CT值分布代表了线衰减系数μ分布,当材料组成相似时CT值与材料密度在一定范围内呈线性关系,通过测试材料截面CT值,利用标准密度材料进行标定,可以获得待测材料的截面密度。

图1为CT同步扫描测试模型,1#、2#为标准密度件,0#为待测密度件,它们的体密度分别为ρ1、ρ2、ρ0,可以采用排水法测得。待测件与标准件具有相同回转构型,三者圆心处于同一圆环上且呈正三角形排列,三角形中心与CT旋转台中心重合,当任一高度位置进行CT扫描时,三者均具有相同的截面形状和面积。

图1CT同步扫描测试模型

Fig.1CT model of synchronous scanning

对于组成相同的0#、1#、2#材料,在任意i截面上三者的截面密度分别为ρ0,i、ρ1,i、ρ2,i,截面CT值分别为CT0,i、CT1,i、CT2,i,CTair,i为空气CT值,0#、1#、2#间的截面密度与截面CT值满足下面公式:

(2)

分别测试0#、1#、2#不同高度的截面CT值,减去不同截面的CTair,i后,根据1#、2#的截面密度(分别等于其体密度),用式(3)计算获得0#不同高度的截面密度。

ρ0,i=ρ1,i+(ρ2,i-ρ1,i)×(CT0,i-CT1,i)/(CT2,i-CT1,i)

(3)

对于组成存在差异的0#、1#、2#材料,CT值与密度ρ的线性关系会受到影响,根据式(3)计算的待测件截面密度与真实密度会存在偏差。由于真实的待测件截面密度加权平均后应该等于其体密度,而体密度可以采用排水法测得,可利用计算出的待测件的体密度ρcal与排水法测得的体密度ρ0间的偏差系数k对ρcal进行线性校正,获得修正后的待测件各截面密度ρcor。

即对于0#待测件,

(4)

设

(5)

有

ρcor,i=k×ρcal,i

(6)

(7)

3 实验方法与CT值提取

设计加工复杂构型的PBX试件(以下简称PBX,体密度为1.890 g·cm-3),标准密度件选择聚甲基丙烯酸甲酯(PMMA,体密度为1.189 g·cm-3)与聚四氟乙烯(PTFE,体密度为2.171 g·cm-3)。采用9 MeV CT(密度分辨率为0.3%)在样件不同高度位置进行水平扇束扫描,本实验选择15个高度,切片厚度为0.2 mm,扫描高度间距为10 mm,扫描后重建获得PBX与PMMA、PTFE的截面CT图像,不同扫描高度样件截面的形状与面积有所变化,图2为不同高度的典型截面CT图像。

a. 140 mmb. 120 mm

c. 80 mmd. 30 mm

图2不同高度的截面CT图像

Fig.2CT images of cross-section at different height

采用MATLAB[18]进行图像格式转换、灰度校正、图像降噪和边缘提取,对不同高度位置的截面CT图像进行数据的阈值分割与提取,圆与圆环的数据区域分割方式见图3,计算每个圆或圆环上8个区域的CT值,将其平均值作为该圆或圆环的CT值,采用相同方法分别获取PMMA、PTFE 、PBX在不同高度位置截面的CT值。

a. circleb. circular ring

图3两种典型截面CT图像MATLAB分割示意图

Fig.3Schematic diagram of CT image segmentation of two typical cross-section by MATLAB

4 实验结果与讨论

4.1 不同扫描时段对CT值的影响

在不同的扫描时段,X射线能量会产生波动,对CT值测试具有一定影响。为考察不同扫描时段,CT能量波动对CT值测试的影响,对PMMA、PTFE、PBX试件选择一个高度截面进行4次截面扫描,每次扫描间隔时间约25 min,采用MATLAB编程,获取第1～4次扫描的PMMA、PTFE、PBX的8个分割区域的CT值,结果见表1。空气灰度值为整个图像空白区域的灰度均值,区域1～8的CT值已减去空气灰度值CTair。计算第1～4次扫描环带CT值的平均值,环带CT值与平均值的相对偏差。

从表1看出，在第1～4次扫描过程中，区域1～8由于距旋转中心距离不一，PMMA、PTFE、PBX的CT值均有变化，且三种材料的变化规律基本一致，每种材料的环带CT值在4次扫描中变化均不大，与平均值的相对偏差为-0.377%～0.214%，空气CT值的变化也不大，这表明在第1～4次扫描过程中，X射线能量输出比较稳定，对PMMA、PTFE、PBX环带CT值的影响幅度小于0.4%，某截面上任一时刻X射线扫描获得的PMMA、PTFE、PBX截面CT值用于密度值计算均具有代表性。

表1不同扫描时段样件的CT值

Table1CT values of samples at different time

samplescanningtimePBX1234PTFE1234PMMA1234airblank346634563450348134663456345034813466345634503481region1186731855518611186112036820263203472038812918128291288112870region2176241753117645176831913019048191281920311611115621162011643region3169671691816967170231883918729187871882511528114411151211502region4182981819918264182482048920350204202038513166130811313713128region5188471874018793187622011519983200872009913015129481299212978region6170291696117074171141872818645187011877411444113931141011428region7172831723017302173711940419279193601935211710116681172211753region8185441844218507184912044820294203762039212936128511291112931Cross-section179081782217895179131969019574196511967712291122221227312279meanvalue178851964812266RD/%0.130-0.3520.0580.1560.214-0.3770.0130.1490.203-0.3620.0580.107

Note: The CT value of cross-section is the average of all region 1～8. The mean value is the average of cross-section in four scanning times. The RD (relative deviation) is calculated between cross-section value and mean value.

4.2 PMMA、PTFE、PBX 截面CT值的获取

利用MATLAB软件编程，逐层获取PMMA、PTFE、PBX在不同高度的截面CT值，分别计算三种材料的截面CT值的平均值，以及截面CT值与平均值的相对偏差，结果见表2。

从表2看出，随着扫描高度变化，PMMA、PBX、PTFE试件的截面CT值也会发生变化，三种材料增大或减小的趋势和幅度基本一致。这表明: (1)采用同步CT扫描模式，X射线能量波动和形状面积变化三种材料CT值测试影响的程度基本相同，可以在测试模型中作为系统误差进行修正，不会对截面密度分布的测试造成影响; (2)受不同时间段X射线能量波动与截面形状面积改变的影响，样件的截面CT值会产生变化，并且不同高度截面的CT值与平均值的相对偏差多数大于1%，而表1中X射线能量波动影响导致截面CT值变化的偏差幅度小于0.4%，说明不同高度截面的CT值变化主要源于截面形状面积的变化。

表2不同扫描高度PMMA、PTFE、PBX试件的截面CT值和相对偏差

Table2CT value and relative deviation of the cross-section of PMMA, PTFE and PBX components on different height

height/mmPMMACTvalueRD/%PTFECTvalueRD/%PBXCTvalueRD/%101370012.29210757.66194488.9820125652.98200742.55184043.133012056-1.19195760.01178820.214011716-3.9819153-2.1617478-2.065011523-5.5618649-4.7316963-4.946011657-4.4618819-3.8617136-3.977011592-4.9918667-4.6416993-4.778011714-3.9918845-3.7317165-3.819011867-2.7419067-2.6017355-2.7510011984-1.7819233-1.7517521-1.82110122690.56196670.47178860.23120124041.66198721.52180841.34130126033.29202153.27183873.04140126333.54202893.65184303.28150127374.39204264.35185433.91Meanvalue12201/19575/17845/

Note: The RD is calculated between CT value at different height and mean value of every material.

4.3 PBX 截面密度分布计算

PMMA、PTFE可视为密度均匀体，各个截面密度等同体密度，利用PMMA和PTFE的截面CT值与截面密度，依据表2数据和式(3)分别计算得到PBX在不同扫描高度的截面密度，结果见表3。根据PBX 15个截面密度ρcal,i和式(4)，计算获得PBX截面平均密度(即体密)ρcal为1.941 g·cm-3，而排水法测试的PBX体密度ρ0为1.890 g·cm-3，二者之间的相对偏差为: (1.941-1.890)/1.890×100%=2.70%。

这种偏差主要来自于标准件与待测件的组成差异。利用式(5)计算出k=0.9737，利用式(6)进行线性修正获得PBX各截面的修正密度ρcor,i，各截面修正密度的均值满足式(7)，计算PBX截面修正密度与排水法测得体密度(ρ0)1.890 g·cm-3间的相对偏差，结果见表3。

从表3看出，利用PMMA、PTFE的CT值-密度关系以及k值线性修正获得的PBX截面修正密度ρcor,i，密度值在PBX的体密度1.890 g·cm-3附近，与体密度的最大偏差为0.74%，截面密度分布具有较好的均匀性，表明采用同步扫描方式并利用排水法密度修正测试复杂构型试件截面密度的方法可行，如果采用更多高度位置的截面扫描，获得的修正k值更加准确，获得的截面密度将更接近真实的截面密度。此外，获得的复杂构型PBX截面密度随着高度增加有减小的趋势，初步分析可能是PBX成型过程导致内外型面与里层部位的密度差异，以后将进一步研究其成因。

表3不同扫描高度PBX的截面密度

Table3PBX cross-section density on different height

height/mmρcal/g·cm-3ρcor/g·cm-3ρ0/g·cm-3RD/%101.9541.9031.8900.740201.9531.9021.8900.619301.9501.8991.8900.474401.9501.8991.8900.476501.9391.8881.890-0.100601.9401.8901.890-0.019701.9391.8881.890-0.100801.9401.8891.890-0.049901.9381.8871.890-0.1601001.9391.8891.890-0.0781101.9351.8841.890-0.3101201.9361.8851.890-0.2431301.9351.8851.890-0.2801401.9331.8821.890-0.4151501.9311.8801.890-0.520meanvalue1.9411.8901.8900.002

5 结 论

(1) 建立了双标准CT同步对称截面扫描模型与实验技术，标准密度件采用与PBX相同的构型，可以获得较为准确的截面CT值，较好解决截面形状、面积变化产生的CT值测试差异问题，提高了PBX截面密度测试的准确性。

(2) X射线能量波动对截面CT值测试的影响较小，不同截面CT值变化主要受截面形状和面积的影响。

(3) 采用双标准CT同步对称截面扫描模型计算获得PBX的截面密度与体密度，采用排水法测得的PBX体密度计算k值进行线性修正的方法，可以获得复杂构型PBX不同高度的截面密度。不同高度截面密度与排水法测得PBX体密度1.890 g·cm-3的最大相对偏差为0.74%。

参考文献:

[1] Perkins D E, Martz H E, Hester L O, et al. Computed tomography experiments of Pantex high explosives[R].UCRL-CR-110256,Lawrence Livermore National Laboratory,Livermore,CA 94551, USA.

[2] 文尚刚, 王胜强, 黄文斌, 等. 密度对压装B炸药燃烧转爆轰性能的影响[J]. 火炸药学报, 2006, 29(5): 5-8.

WEN Shang-gang, WANG Sheng-qiang, HUANG Wen-bin, et al. The effect of density in composition B on deflagration-detonation-transition behavior[J].ChineseJournalofExplosives&Propellants, 2016, 29(5): 5-8.

[3] 王玮, 王建灵, 郭炜, 等. 装药密度及尺寸对RDX基含铝炸药爆压爆速的影响[J]. 含能材料, 2010, 18(5): 563-567.

WANG Wei, WANG Jian-ling, GUO Wei, et al. Effect of charge density and size on detonation pressure and detonation velocity of RDX-based aluminized explosive[J].ChineseJournalofEnergeticMaterals(HannengCailiao), 2010, 18(5): 563-567.

[4] Roger A Morris, Lowell D Harris .High-precision density measurements using tomography[R]. LA-UR-90-2988,Los Alamos National Laboratory,NM 87545,USA.

[5] 王功庆, 杨锦清, 曾寄萍. 炮弹体内炸药密度的测定[J]. 电子学与探测技术, 1999, 19(3): 172-177.

WANG Gong-qing, YANG Jin-qing, ZENG Ji-ping. ICT nondestructive testing for high explosive density[J].NuclearElectronics&DetectionTechnology, 1999, 19(3): 172-177.

[6] 杨文海, 何得昌, 徐军培. γ射线工业CT技术在高能炸药密度检测中的应用[J].火炸药学报, 2001, 24(3): 33-34, 72.

YANG Wen-hai, HE De-chang, XU Jun-pei. Determination of local TNT density in cannon projectile casing by γ ray non destructive inspection[J].ChineseJournalofExplosives&Propellants, 2001, 24(3): 33-34, 72.

[7] 张洪伟, 潘保青.基于ICT技术的弹药无损检测设备[J].飞行器测控学报, 2006, 25(1): 76-79.

ZHANG Hong-wei, PAN Bao-qing. ICT-based ammunition NDT instrument[J].ChineseJournalofSpacecraftTT&CTechnology, 2006, 25(1): 76-79.

[8] ASTM E 1935-1997(Reapproved 2008), Standard test method for calibrating and measuring CT density[S].ASTM International,100 Barr Harbor Drive,PO BOX C700,West Conshohocken,PA 19428-2959,United States.

[9] Sinka I C, Burch S F, Tweed J H, et al.Measurement of density variations in tablets using X-ray computed tomography[J].InternationalJournalofPharmaceutics, 2004 (271): 215-224.

[10] Ohashi J, Arakawa H.Use of X-ray CT scanner for density distribution measurements in powder compacts[J].JJapanSocPowderandPowderMetallurgy, 1995, 42(9): 1032-1036.

[11] 韩书霞, 于雷, 孔超, 等. 基于计算机断层扫描技术的木材密度检测研究[J]. 森林工程, 2007, 23(1): 19-23.

HAN Shu-xia, YU Lei, KONG Chao, et al. Research of wood density testing based on computer tomography technology[J].ChineseJournalofForestEngineering, 2007, 23(1): 19-23.

[12] 刘红林, 金志浩, 郝志彪,等. ICT技术测试炭/炭复合材料内部密度分布[J]. 无损检测, 2007, 29(12): 726-727,735.

LIU Hong-lin, JIN Zhi-hao, HAO Zhi-biao, et al. A new method for testing the density of carbon/carbon composites[J].ChineseJournalofNondestrctiveTesting, 2007, 29(12): 726-727,735.

[13] 单雯雯, 张林, 徐家云. 聚丙烯酸酯泡沫密度均匀性的射线检测技术[J]. 强激光与离子束, 2007, 19(1): 71-74.

SHAN Wen-wen, ZHANG Lin, XU Jia-yun. Detecting technology for density uniformity of trimethylo-lpropane trimethacrylate polymer foams[J].ChineseJournalofHighPowerLaserandParticleBeams, 2007, 19(1): 71-74.

[14] 李俊杰, 韩焱, 王黎明. 基于ICT的石墨密度均匀性检测方法研究[J].无损检测, 2008, 30(3):163-164,192.

LI Jun-jie, HAN yan, WANG Li-ming. Study on the method of graphite density uniformity inspection based on X-ray ICT[J].ChineseJournalofNondestrctiveTesting, 2008, 30(3): 163-164,192.

[15] 张伟斌，田勇，杨仍才, 等. RDX晶体颗粒压制密度分布的μCT实验研究[J]. 含能材料, 2012, 20(5): 565-570.

ZHANG Wei-bin, TIAN Yong, YANG Reng-cai, et al. Density distribution of pressed RDX crystals with a cone-bean micro-focus computed tomography[J].ChineseJournalofEnergeticMaterals(HannengCailiao), 2012, 20(5): 565-570.

[16] 杨雪海，张伟斌，戴斌, 等. 微焦点CT在含能材料密度测试中的应用研究[J]. CT理论与应用研究, 2009, 18(4): 61-67.

YANG Xue-hai, ZHANG Wei-bin, DAI Bin, et al. Research on energetic material density by X-ray micro computed tomography[J].ChineseJournalofComputerizedTomographyTheoryandApplications, 2009, 18(4): 61-67.

[17] 庄天戈.CT原理与算法应用[M].上海交通大学出版社, 1992.

[18] The MathWorks,Inc. 3 Apple Hill Drive[CP],Natick,MA 01760-2098,USA.