CN108168996A - 一种基于声发射信号的pbx损伤演化ct原位分析方法 - Google Patents
一种基于声发射信号的pbx损伤演化ct原位分析方法 Download PDFInfo
- Publication number
- CN108168996A CN108168996A CN201711332249.2A CN201711332249A CN108168996A CN 108168996 A CN108168996 A CN 108168996A CN 201711332249 A CN201711332249 A CN 201711332249A CN 108168996 A CN108168996 A CN 108168996A
- Authority
- CN
- China
- Prior art keywords
- crackle
- pbx
- acoustic emission
- scan
- image
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 73
- 230000006378 damage Effects 0.000 title claims abstract description 53
- 238000011065 in-situ storage Methods 0.000 title claims abstract description 39
- 238000011161 development Methods 0.000 title claims abstract description 25
- 238000002591 computed tomography Methods 0.000 claims abstract description 105
- 239000000463 material Substances 0.000 claims abstract description 40
- 230000008569 process Effects 0.000 claims abstract description 25
- 238000012544 monitoring process Methods 0.000 claims abstract description 20
- 238000004458 analytical method Methods 0.000 claims abstract description 14
- 238000004445 quantitative analysis Methods 0.000 claims abstract description 14
- 238000012800 visualization Methods 0.000 claims abstract description 11
- 208000037656 Respiratory Sounds Diseases 0.000 claims description 90
- 238000012360 testing method Methods 0.000 claims description 58
- 230000018109 developmental process Effects 0.000 claims description 21
- 230000009466 transformation Effects 0.000 claims description 19
- 239000011159 matrix material Substances 0.000 claims description 15
- 239000002360 explosive Substances 0.000 claims description 8
- 238000012512 characterization method Methods 0.000 claims description 6
- 230000001360 synchronised effect Effects 0.000 claims description 6
- 230000008859 change Effects 0.000 claims description 5
- 230000001186 cumulative effect Effects 0.000 claims description 5
- 238000001514 detection method Methods 0.000 claims description 5
- 239000000284 extract Substances 0.000 claims description 5
- 239000000203 mixture Substances 0.000 claims description 4
- 230000011218 segmentation Effects 0.000 claims description 4
- 230000003321 amplification Effects 0.000 claims description 3
- 238000005422 blasting Methods 0.000 claims description 3
- 230000002708 enhancing effect Effects 0.000 claims description 3
- 239000008187 granular material Substances 0.000 claims description 3
- 239000004519 grease Substances 0.000 claims description 3
- 238000003199 nucleic acid amplification method Methods 0.000 claims description 3
- 229920001296 polysiloxane Polymers 0.000 claims description 3
- 230000002123 temporal effect Effects 0.000 claims description 3
- 229940099259 vaseline Drugs 0.000 claims description 3
- 239000011230 binding agent Substances 0.000 claims 1
- 230000000977 initiatory effect Effects 0.000 abstract description 3
- 229920000642 polymer Polymers 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 3
- 206010011376 Crepitations Diseases 0.000 description 2
- 230000001066 destructive effect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000004880 explosion Methods 0.000 description 2
- 230000000149 penetrating effect Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 201000004569 Blindness Diseases 0.000 description 1
- 240000002853 Nelumbo nucifera Species 0.000 description 1
- 235000006508 Nelumbo nucifera Nutrition 0.000 description 1
- 235000006510 Nelumbo pentapetala Nutrition 0.000 description 1
- 244000131316 Panax pseudoginseng Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000002485 combustion reaction Methods 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 239000003380 propellant Substances 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 238000002604 ultrasonography Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
- G01N23/02—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
- G01N23/04—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material
- G01N23/046—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material using tomography, e.g. computed tomography [CT]
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N3/02—Details
- G01N3/06—Special adaptations of indicating or recording means
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/0058—Kind of property studied
- G01N2203/006—Crack, flaws, fracture or rupture
- G01N2203/0062—Crack or flaws
- G01N2203/0064—Initiation of crack
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/0058—Kind of property studied
- G01N2203/006—Crack, flaws, fracture or rupture
- G01N2203/0062—Crack or flaws
- G01N2203/0066—Propagation of crack
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/02—Details not specific for a particular testing method
- G01N2203/025—Geometry of the test
- G01N2203/0252—Monoaxial, i.e. the forces being applied along a single axis of the specimen
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/02—Details not specific for a particular testing method
- G01N2203/06—Indicating or recording means; Sensing means
- G01N2203/0641—Indicating or recording means; Sensing means using optical, X-ray, ultraviolet, infrared or similar detectors
- G01N2203/0647—Image analysis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/02—Details not specific for a particular testing method
- G01N2203/06—Indicating or recording means; Sensing means
- G01N2203/0658—Indicating or recording means; Sensing means using acoustic or ultrasonic detectors
Landscapes
- Health & Medical Sciences (AREA)
- Immunology (AREA)
- General Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Pathology (AREA)
- General Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Biochemistry (AREA)
- Engineering & Computer Science (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Pulmonology (AREA)
- Radiology & Medical Imaging (AREA)
- Theoretical Computer Science (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
一种基于声发射信号的PBX损伤演化CT原位分析方法，首先搭建CT原位加载及损伤检测实验系统，该系统由单轴原位加载系统、CT扫描系统及声发射监测系统组成；同时实施力学加载和声发射信号的采集，根据声发射信号特征确定CT扫描载荷节点，进行不同加载阶段CT图像采集；提取CT图像中裂纹形态，并进行三维可视化和定量分析；结合材料加载曲线、声发射信号特征、不同加载阶段CT图像裂纹三维特征，定量描述PBX材料受载过程中的损伤演化；本发明方法结合声发射技术可在线监测的特点，利用CT原位观测的方法可以实现对加载过程中PBX材料损伤演化的定量描述，具有准确捕捉裂纹萌生起始点、损伤演化过程可视化的优点，易实现、易操作、效率高，可广泛用于PBX材料损伤演化的分析中。
Description
技术领域
本发明涉及PBX材料损伤演化研究领域，具体涉及一种基于声发射信号的PBX损伤演化CT原位分析方法。
背景技术
PBX(Polymer Bonded Explosive)是聚合物粘结炸药的简称，聚合物粘结炸药是一种以高能炸药颗粒和高聚物粘结剂为主要成分粘结而成的颗粒填充复合材料，在常规武器战斗部及火箭推进剂中具有广泛应用。PBX在使用中不仅承受着复杂的结构载荷，同时作为一种含能材料，还需保证良好的燃烧爆炸功能，因此PBX同时具有特殊功能材料和结构部件的重要作用。
PBX的服役环境非常复杂，在生产、加工、运输、储存和发射过程中，均处于不同加载速率、应力状态之中。PBX在服役过程中会产生裂纹等多种形式的损伤，这些损伤的存在一方面使材料的力学性能劣化，并可能最终导致材料的破坏；另一方面，损伤对材料“热点”的形成也有重要的影响，甚至会影响炸药的感度、燃烧和爆炸性能。因此对于PBX的损伤及力学性能的研究对于PBX材料的结构设计及安全性能评估等都具有重要意义。
一般认为，PBX发生破坏的主要原因是材料内部初始微小裂纹和缺陷在不同的应力作用下发生了一系列的累积性变化而导致最终的断裂破坏，要弄清初始微小裂纹发展成破坏性贯穿裂纹的过程以及扩展机理，就必须从PBX材料在不同应力状态下微裂纹的萌生、发展全过程入手，系统完整地研究PBX材料整个破坏过程的损伤演化规律。考虑PBX材料组分的复杂性及功能的特殊性，目前对于PBX材料损伤的观测手段有显微观察、超声、声发射、CT(计算机断层成像)、中子小角散射等，其中声发射和CT技术广泛应用于工程中PBX的损伤监测和检测中。CT技术可以观测到材料内部的结构形貌及变化过程，为研究不同载荷阶段材料内部的损伤演化过程，传统的方法是将试件加载到某个特定的载荷阶段，停止加载，卸载后进行离线CT扫描。一方面，CT扫描的应力选择依靠经验，扫描载荷节点的选择与材料的损伤状况没有直观的对应关系，扫描载荷节点的盲目选择不仅导致实验成本的增加，而且无法捕捉到微裂纹萌生及扩展过程；另一方面，试件卸载后的应力状态与加载过程中有很大区别，卸载后的裂纹形态与受载状态下的裂纹形态也不同。考虑声发射具有实时观测的特点，依据CT原位加载过程中声发射信号特征确定CT扫描载荷点，通过对CT图像裂纹的三维可视化和定量分析，建立声发射信号特征和裂纹三维特征参数之间的定量关系，有望实现PBX损伤演化过程的三维定量描述，解释PBX材料破坏的本质特征。
针对以上问题，本发明提出了一种基于声发射信号的PBX损伤演化CT原位分析方法，依据PBX材料CT原位加载过程中声发射信号特征确定CT扫描载荷节点，通过对CT扫描图像裂纹的三维可视化和定量分析，建立声发射信号特征和裂纹三维特征参数之间的定量关系，为PBX材料损伤演化提供定量分析方法。
发明内容
为了解决上述现有技术存在的问题，本发明的目的在于提供一种基于声发射信号的PBX损伤演化CT原位分析方法，首先搭建CT原位加载及损伤检测实验系统，该系统由单轴原位加载系统、CT扫描系统及声发射监测系统组成；同时实施力学加载和声发射信号的采集，根据声发射信号特征确定CT扫描载荷节点，进行不同加载阶段CT图像采集；提取CT图像中裂纹形态，并进行三维可视化和定量分析；结合材料加载曲线、声发射信号特征、CT图像裂纹三维特征，定量描述PBX材料受载过程中的损伤演化；本发明方法结合声发射技术可在线监测的特点，利用CT原位观测的方法可以实现对加载过程中PBX材料损伤演化的定量描述，具有准确捕捉裂纹萌生起始点、损伤演化过程可视化的优点，易实现、易操作、效率高，可广泛用于PBX材料损伤演化的分析中。
为达到以上目的，本发明采用如下技术方案：
一种基于声发射信号的PBX损伤演化CT原位分析方法，包括如下步骤：
步骤1：搭建CT原位加载及损伤检测实验系统，该系统主要包括三部分：包括由射线源5、载物台4和探测器6组成的CT扫描系统，由声发射传感器7、信号放大器8、声发射仪9和计算机10依次连接组成的声发射监测系统以及单轴原位加载系统；搭建方法为：在PBX试件2表面涂抹凡士林作为耦合剂与声发射传感器7连接，声发射传感器7输出端通过信号放大器8与声发射仪9连接；在PBX试件2端部及单轴原位加载系统的上压板1和下压板3上涂抹硅脂进行润滑，减小上压板1和下压板3对PBX试件2的端部约束，将PBX试件2放置在下压板3中心位置并使上压板1与PBX试件2完全接触，将单轴原位加载系统固定于CT扫描系统的载物台4上；
步骤2：基于声发射信号的PBX试件单轴加载CT图像采集，具体步骤如下：
1)通过断铅试验确定合适的声发射监测参数值，包括采样频率、放大器增益、检测门槛、峰值鉴别时间PDT、波击鉴别时间HDT和波击锁闭时间HLT；
2)设置CT扫描系统扫描参数，包括扫描电压、扫描电流、放大倍数、曝光时间和扫描帧数；
3)开启射线源5对PBX试件2进行第一次扫描，记录PBX试件2的初始状态；
4)同时运行声发射监测系统和单轴原位加载系统，同步实施声发射信号的采集和对PBX试件2的加载；以声发射信号持续时间作为声发射基本参数监测PBX试件2的加载过程，当首次监测到声发射信号持续时间大于0时，停止加载并进行载荷保持，进行第二次CT扫描，扫描过程中声发射监测系统保持采集状态；
5)第二次CT扫描结束后继续进行加载，当再次监测到声发射信号持续时间大于0时，停止加载并进行载荷保持，进行第三次CT扫描，如此重复加载及扫描过程，直到试件发生破坏，进行最后一次CT扫描；
步骤3：基于CT图像处理的损伤可视化及定量分析，具体步骤如下：
1)PBX试件2的CT扫描图像中灰度值最高的区域为炸药颗粒，灰度值较低的区域为高聚物粘结剂，灰度值最低的区域为裂纹，根据CT扫描图像特点，在MATLAB环境下进行CT扫描图像的预处理；首先利用“imadjust”函数对CT扫描图像进行灰度变换；然后利用“imtophat”函数和“imbothat”函数对灰度变换后的图像分别进行高帽变换和低帽变换；其次利用高帽变换后的图像与低帽变换后的图像相减，进行图像的增强；最后利用最大类间方差法对增强后的图像进行局部阈值分割，提取出裂纹的二维图像；
2)将步骤1)中提取出来的裂纹的二维图像进行三维堆叠，得到裂纹三维图像，即裂纹在空间上的三维形态；
3)根据步骤2)获得的裂纹三维图像获取裂纹三维特征参数，包括裂纹最大宽度、裂纹面面积、裂纹体积、裂纹各向异性度、裂纹三维计盒维数；
4)利用体视学的方法对步骤2)中获得的裂纹三维图像进行定量的分析：首先定义一个球形的分析区域，用一系列不同方向上的射线穿过该分析区域，得到每个方向上的平均截断长度；使各个方向上的直线穿过当前分析区域坐标系的原点，每个方向上线段的长度在数值上等于该方向上的平均截断长度，得到平均截断长度在三维空间上的分布；用一个三维椭球拟合平均截断长度在空间上的分布，椭球的方程写作
其中：x1,x2,x3为椭球面上一点在当前分析区域坐标系中的三个坐标分量，A,B,C,D,E,F为椭球方程的系数；
将椭球方程的系数写作一个表征材料各向异性的矩阵M：
计算矩阵M的特征值和特征向量，矩阵M的特征向量表示椭球三个主轴的方向，即与当前分析区域定义坐标主轴的夹角；矩阵M的特征值表示每个主轴方向平均截断长度的相对大小；用矩阵M的最大特征值和最小特征值定义表征材料各项异性程度的标量参数DA：
DA从0到1变化来表示裂纹分布的各向异性程度，DA值越大，裂纹空间分布的各向异性越大；
5)计算步骤3)中裂纹三维计盒维数：首先读入三维裂纹图像，存为三维矩阵Q，构成样本空间；用一系列边长为ε的正方体盒子覆盖裂纹所在空间，计算覆盖裂纹所在空间需要的盒子数N；对log(N)和log(1/ε)进行线性拟合，所得直线的斜率即为该裂纹三维空间的计盒维数；计盒维数值越大，说明裂纹的三维空间分布越复杂；
步骤4：建立声发射信号特征和裂纹三维特征参数之间的关系，定量描述材料受载过程中的损伤演化，具体步骤如下：
1)以声发射事件累计数作为声发射特征参数分析PBX试件2加载过程中的声发射信号特征，根据时间上的同步关系，以PBX试件2所受载荷为横坐标，声发射事件累计数为纵坐标，绘制声发射事件累计数与PBX试件2所受载荷的关系图；
2)对应于步骤1)中得到的关系图，将CT扫描图像与CT扫描载荷节点对应，得到不同加载阶段对应的声发射事件累计数及CT扫描图像中裂纹的三维特征参数；对应声发射事件累计数与CT扫描图像中裂纹的裂纹最大宽度、裂纹面面积、裂纹体积、裂纹各向异性度及裂纹三维计盒维数，定量描述PBX试件2在加载过程中声发射事件累积数及裂纹三维特征参数的变化过程。
所述的CT扫描为CT原位载荷保持扫描。
所述的CT为纳米焦点或者微米焦点工业CT。
所述的CT扫描电压为100kV，扫描电流为100μA。
监测过程所用的声发射基本参数采用声发射信号持续时间，分析过程所用的声发射特征参数采用声发射事件累计数。
和现有技术相比，本发明的优点如下：
1)本发明针对PBX损伤演化过程的分析，结合CT原位技术和声发射技术的优点，依据PBX试件CT原位单轴加载过程中声发射信号特征准确确定CT扫描载荷节点，精确捕捉裂纹萌生及扩展过程，建立声发射信号特征和裂纹三维特征参数的定量关系，可以实现PBX材料损伤演化过程的三维定量描述，解释PBX材料破坏的本质特征。
2)本发明提出的CT原位加载并同时进行声发射信号采集和CT图像采集的方法能够避免传统离线检测方法中卸载后裂纹形态发生改变的缺点，能够反映材料真实受载状态下的损伤状况；
3)本发明建立了适合PBX材料CT扫描图像的裂纹提取及定量分析方法，通过对CT扫描图像的处理，提取出PBX试件2中的裂纹并进行可视化和定量分析，裂纹参数的定量化为CT扫描图像与声发射信号特征的结合提供了基础。
附图说明
图1为本发明方法所用CT原位加载及损伤检测实验系统示意图。
图2为本发明方法中CT图像裂纹提取过程，其中：图2a为CT扫描图像原图，图2b为高帽变换后的图像与低帽变换后的图像相减得到的图像，图2c为局部阈值分割后的图像。
图3为本发明方法中裂纹三维可视化图像。
图4为本发明方法建立的载荷、声发射事件累计数与裂纹三维特征参数定量关系示意图。
具体实施方式
如图1所示，本发明方法的检测步骤为：首先搭建CT原位加载及损伤检测实验系统，该系统由单轴原位加载系统、CT扫描系统及声发射监测系统组成；同时实施力学加载和声发射信号的采集，根据声发射信号特征确定CT扫描载荷节点，进行不同加载阶段CT图像采集；提取CT扫描图像中裂纹形态，并进行三维可视化和定量分析；结合材料加载曲线、声发射信号特征、CT图像裂纹三维特征，定量描述PBX试件2受载过程中的损伤演化。
下面结合图1、图2、图3、图4和具体实施例对本发明作进一步的详细描述。
本发明基于声发射信号的PBX损伤演化CT原位分析方法，包括如下步骤：
步骤1：搭建如图1所示的CT原位加载及损伤检测实验系统，主要包括三部分：由射线源5、载物台4和探测器6组成的CT扫描系统，由声发射传感器7、信号放大器8、声发射仪9和计算机10依次连接组成的声发射监测系统以及单轴原位加载系统；在PBX试件2表面涂抹凡士林作为耦合剂与声发射传感器7连接，声发射传感器7输出端通过信号放大器8与声发射仪9连接；在PBX试件2端部及单轴原位加载系统的上压板1和下压板3上涂抹硅脂进行润滑，减小上压板1和下压板3对PBX试件2的端部约束，将PBX试件2放置在下压板3中心位置并使上压板1与PBX试件2完全接触，将单轴原位加载系统固定于CT扫描系统的载物台4上；
步骤2：基于声发射信号的PBX试件单轴加载CT图像采集，具体步骤如下：
1)通过断铅试验确定合适的声发射监测参数值，包括采样频率、放大器增益、检测门槛、峰值鉴别时间PDT、波击鉴别时间HDT和波击锁闭时间HLT；
2)设置CT扫描系统扫描参数，包括扫描电压、扫描电流、放大倍数、曝光时间和扫描帧数；
3)开启射线源5对PBX试件2进行第一次扫描，记录PBX试件2的初始状态；
4)同时运行声发射监测系统和单轴原位加载系统，同步实施声发射信号的采集和对PBX试件2的加载；以声发射信号持续时间作为声发射基本参数监测PBX试件2的加载过程，当首次监测到声发射信号持续时间大于0时，停止加载并进行载荷保持，进行第二次CT扫描，扫描过程中声发射监测系统保持采集状态；
5)第二次CT扫描结束后继续进行加载，当再次监测到声发射信号持续时间大于0时，停止加载并进行载荷保持，进行第三次CT扫描，如此重复加载及扫描过程，直到试件发生破坏，进行最后一次CT扫描；
步骤3：基于CT图像处理的损伤可视化及定量分析，具体步骤如下：
1)PBX试件2的CT扫描图像如图2a所示，图中灰度值最高的区域为炸药颗粒，灰度值较低的区域为高聚物粘结剂，灰度值最低的区域为裂纹，根据CT扫描图像特点，在MATLAB环境下进行CT扫描图像的预处理；首先利用“imadjust”函数对CT扫描图像进行灰度变换；然后利用“imtophat”函数和“imbothat”函数对灰度变换后的图像分别进行高帽变换和低帽变换；其次利用高帽变换后的图像与低帽变换后的图像相减，进行图像的增强，得到如图2b所示的的图像；最后利用最大类间方差法对图像进行局部阈值分割，提取出裂纹的二维图像如图2c所示；
2)将步骤1)中提取出来的裂纹的二维图像进行三维堆叠，得到裂纹三维图像如图3所示，即裂纹在空间上的三维形态；
3)根据步骤2)获得的裂纹三维图像获取裂纹三维特征参数，包括裂纹最大宽度、裂纹面面积、裂纹体积、裂纹各向异性度、裂纹三维计盒维数；
4)利用体视学的方法对步骤2)中获得的裂纹三维图像进行定量的分析：首先定义一个球形的分析区域，用一系列不同方向上的射线穿过该分析区域，得到每个方向上的平均截断长度；使各个方向上的直线穿过当前分析区域坐标系的原点，每个方向上线段的长度在数值上等于该方向上的平均截断长度，得到平均截断长度在三维空间上的分布；用一个三维椭球拟合平均截断长度在空间上的分布，椭球的方程写作
其中：x1,x2,x3为椭球面上一点在当前分析区域坐标系中的三个坐标分量，A,B,C,D,E,F为椭球方程的系数；
将椭球方程的系数写作一个表征材料各向异性的矩阵M：
计算矩阵M的特征值和特征向量，矩阵M的特征向量表示椭球三个主轴的方向，即与当前分析区域定义坐标主轴的夹角；矩阵M的特征值表示每个主轴方向平均截断长度的相对大小；用矩阵M的最大特征值和最小特征值定义表征材料各项异性程度的标量参数DA：
DA从0到1变化来表示裂纹分布的各向异性程度，DA值越大，裂纹空间分布的各向异性越大；
5)计算步骤3)中裂纹三维计盒维数：首先读入三维裂纹图像，存为三维矩阵Q，构成样本空间；用一系列边长为ε的正方体盒子覆盖裂纹所在空间，计算覆盖裂纹所在空间需要的盒子数N；对log(N)和log(1/ε)进行线性拟合，所得直线的斜率即为该裂纹三维空间的计盒维数；计盒维数值越大，说明裂纹的三维空间分布越复杂；
步骤4：建立声发射信号特征和裂纹三维特征参数之间的关系，定量描述材料受载过程中的损伤演化，具体步骤如下：
1)以声发射事件累计数作为声发射特征参数分析PBX试件2加载过程中的声发射信号特征，根据时间上的同步关系，以PBX试件2所受载荷为横坐标，声发射事件累计数为纵坐标，绘制声发射事件累计数与PBX试件2所受载荷的关系图；
2)对应于步骤1)中得到的关系图，将CT扫描图像与CT扫描载荷节点对应，得到不同加载阶段对应的声发射事件累计数及CT扫描图像中裂纹三维特征参数；对应声发射事件累计数与CT扫描图像中裂纹的裂纹最大宽度、裂纹面面积、裂纹体积、裂纹各向异性度及裂纹三维计盒维数，定量描述PBX试件2在加载过程中声发射事件累积数及裂纹三维特征参数的变化过程；如图4所示，随着载荷的增加，PBX试件2加载过程中声发射事件累积数随载荷的增加分为四个阶段；PBX试件2在载荷a点开始进行加载，在载荷b点裂纹开始萌生，随着载荷的增加，在载荷c点裂纹开始扩展，并在载荷d点达到能量临界释放点，裂纹开始迅速扩展并贯通整个PBX试件2，在载荷e点发生宏观破坏，a～e点对应着不同的CT扫描图像及裂纹三维特征参数。
作为本发明的优选实施方式，所述的CT扫描为CT原位载荷保持扫描。
作为本发明的优选实施方式，所述的CT为纳米焦点或者微米焦点工业CT。
作为本发明的优选实施方式，所述的CT扫描电压为100kV，扫描电流为100μA。
作为本发明的优选实施方式，监测过程所用的声发射基本参数采用声发射信号持续时间，分析过程所用的声发射特征参数采用声发射事件累计数。
Claims (5)
1.一种基于声发射信号的PBX损伤演化CT原位分析方法，其特征在于：包括如下步骤：
步骤1：搭建CT原位加载及损伤检测实验系统，该系统主要包括三部分：包括由射线源5、载物台4和探测器6组成的CT扫描系统，由声发射传感器7、信号放大器8、声发射仪9和计算机10依次连接组成的声发射监测系统以及单轴原位加载系统；搭建方法为：在PBX试件2表面涂抹凡士林作为耦合剂与声发射传感器7连接，声发射传感器7输出端通过信号放大器8与声发射仪9连接；在PBX试件2端部及单轴原位加载系统的上压板1和下压板3上涂抹硅脂进行润滑，减小上压板1和下压板3对PBX试件2的端部约束，将PBX试件2放置在下压板3中心位置并使上压板1与PBX试件2完全接触，将单轴原位加载系统固定于CT扫描系统的载物台4上；
步骤2：基于声发射信号的PBX试件单轴加载CT图像采集，具体步骤如下：
1)通过断铅试验确定声发射监测参数值，包括采样频率、放大器增益、检测门槛、峰值鉴别时间PDT、波击鉴别时间HDT和波击锁闭时间HLT；
2)设置CT扫描系统扫描参数，包括扫描电压、扫描电流、放大倍数、曝光时间和扫描帧数；
3)开启射线源5对PBX试件2进行第一次扫描，记录PBX试件2的初始状态；
4)同时运行声发射监测系统和单轴原位加载系统，同步实施声发射信号的采集和对PBX试件2的加载；以声发射信号持续时间作为声发射基本参数监测PBX试件2的加载过程，当首次监测到声发射信号持续时间大于0时，停止加载并进行载荷保持，进行第二次CT扫描，扫描过程中声发射监测系统保持采集状态；
5)第二次CT扫描结束后继续进行加载，当再次监测到声发射信号持续时间大于0时，停止加载并进行载荷保持，进行第三次CT扫描，如此重复加载及扫描过程，直到试件发生破坏，进行最后一次CT扫描；
步骤3：基于CT图像处理的损伤可视化及定量分析，具体步骤如下：
1)PBX试件2的CT扫描图像中灰度值最高的区域为炸药颗粒，灰度值较低的区域为高聚物粘结剂，灰度值最低的区域为裂纹，根据CT扫描图像特点，在MATLAB环境下进行CT扫描图像的预处理；首先利用“imadjust”函数对CT扫描图像进行灰度变换；然后利用“imtophat”函数和“imbothat”函数对灰度变换后的图像分别进行高帽变换和低帽变换；其次利用高帽变换后的图像与低帽变换后的图像相减，进行图像的增强；最后利用最大类间方差法对增强后的图像进行局部阈值分割，提取出裂纹的二维图像；
2)将步骤1)中提取出来的裂纹的二维图像进行三维堆叠，得到裂纹三维图像，即裂纹在空间上的三维形态；
3)根据步骤2)获得的裂纹三维图像获取裂纹三维特征参数，包括裂纹最大宽度、裂纹面面积、裂纹体积、裂纹各向异性度、裂纹三维计盒维数；
4)利用体视学的方法对步骤2)中获得的裂纹三维图像进行定量的分析：首先定义一个球形的分析区域，用一系列不同方向上的射线穿过该分析区域，得到每个方向上的平均截断长度；使各个方向上的直线穿过当前分析区域坐标系的原点，每个方向上线段的长度在数值上等于该方向上的平均截断长度，得到平均截断长度在三维空间上的分布；用一个三维椭球拟合平均截断长度在空间上的分布，椭球的方程写作
其中：x1,x2,x3为椭球面上一点在当前分析区域坐标系中的三个坐标分量，A,B,C,D,E,F为椭球方程的系数；
将椭球方程的系数写作一个表征材料各向异性的矩阵M：
计算矩阵M的特征值和特征向量，矩阵M的特征向量表示椭球三个主轴的方向，即与当前分析区域定义坐标主轴的夹角；矩阵M的特征值表示每个主轴方向平均截断长度的相对大小；用矩阵M的最大特征值和最小特征值定义表征材料各项异性程度的标量参数DA：
DA从0到1变化来表示裂纹分布的各向异性程度，DA值越大，裂纹空间分布的各向异性越大；
5)计算步骤3)中裂纹三维计盒维数：首先读入三维裂纹图像，存为三维矩阵Q，构成样本空间；用一系列边长为ε的正方体盒子覆盖裂纹所在空间，计算覆盖裂纹所在空间需要的盒子数N；对log(N)和log(1/ε)进行线性拟合，所得直线的斜率即为该裂纹三维空间的计盒维数；计盒维数值越大，说明裂纹的三维空间分布越复杂；
步骤4：建立声发射信号特征和裂纹三维特征参数之间的关系，定量描述材料受载过程中的损伤演化，具体步骤如下：
1)以声发射事件累计数作为声发射特征参数分析PBX试件2加载过程中的声发射信号特征，根据时间上的同步关系，以PBX试件2所受载荷为横坐标，声发射事件累计数为纵坐标，绘制声发射事件累计数与PBX试件2所受载荷的关系图；
2)对应于步骤1)中得到的关系图，将CT扫描图像与CT扫描载荷节点对应，得到不同加载阶段对应的声发射事件累计数及CT扫描图像中裂纹三维特征参数；对应声发射事件累计数与CT扫描图像中裂纹的裂纹最大宽度、裂纹面面积、裂纹体积、裂纹各向异性度及裂纹三维计盒维数，定量描述PBX试件2在加载过程中声发射事件累积数及裂纹三维特征参数的变化过程。
2.根据权利要求1所述的一种基于声发射信号的PBX损伤演化CT原位分析方法，其特征在于：所述的CT扫描为CT原位载荷保持扫描。
3.根据权利要求1所述的一种基于声发射信号的PBX损伤演化CT原位分析方法，其特征在于：所述的CT为纳米焦点或者微米焦点工业CT。
4.根据权利要求1所述的一种基于声发射信号的PBX损伤演化CT原位分析方法，其特征在于：所述的CT扫描电压为100kV，扫描电流为100μA。
5.根据权利要求1所述的一种基于声发射信号的PBX损伤演化CT原位分析方法，其特征在于：监测过程所用的声发射基本参数采用声发射信号持续时间，分析过程所用的声发射特征参数采用声发射事件累计数。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711332249.2A CN108168996B (zh) | 2017-12-13 | 2017-12-13 | 一种基于声发射信号的pbx损伤演化ct原位分析方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711332249.2A CN108168996B (zh) | 2017-12-13 | 2017-12-13 | 一种基于声发射信号的pbx损伤演化ct原位分析方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108168996A true CN108168996A (zh) | 2018-06-15 |
| CN108168996B CN108168996B (zh) | 2019-03-12 |
Family
ID=62525157
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711332249.2A Active CN108168996B (zh) | 2017-12-13 | 2017-12-13 | 一种基于声发射信号的pbx损伤演化ct原位分析方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108168996B (zh) |
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111351810A (zh) * | 2020-03-30 | 2020-06-30 | 山东省分析测试中心 | 一种含多尺度缺陷金属断裂失效行为的分析方法 |
| CN111678993A (zh) * | 2020-05-26 | 2020-09-18 | 中国电子产品可靠性与环境试验研究所（（工业和信息化部电子第五研究所）（中国赛宝实验室）） | 无损检测系统以及方法 |
| CN111707548A (zh) * | 2020-06-30 | 2020-09-25 | 中国矿业大学 | 一种基于裂隙和能量演化的锚杆支护效果评价方法 |
| CN111999391A (zh) * | 2020-06-24 | 2020-11-27 | 天津科技大学 | 声发射监测装置及其对材料结构的疲劳损伤分析方法 |
| CN112800623A (zh) * | 2020-07-22 | 2021-05-14 | 北京理工大学 | 基于反应演化模型的带泄压结构弹药的反应烈度评估方法 |
| CN113063789A (zh) * | 2021-03-18 | 2021-07-02 | 山东建筑大学 | 一种盾构隧道管片内部微破坏检测装置及方法 |
| CN113917562A (zh) * | 2021-09-27 | 2022-01-11 | 中国科学院武汉岩土力学研究所 | 深埋软弱夹层宏细观结构表征及三维空间构造方法及装置 |
| CN113984500A (zh) * | 2021-10-27 | 2022-01-28 | 中国航发北京航空材料研究院 | 一种材料失效监测装置及其分类方法 |
| CN114813798A (zh) * | 2022-05-18 | 2022-07-29 | 中国工程物理研究院化工材料研究所 | 用于表征材料内部结构及成分的ct检测装置和成像方法 |
| CN115043118A (zh) * | 2021-09-09 | 2022-09-13 | 海沃机械（中国）有限公司 | 用于垃圾转运车的拉臂部件的损伤识别方法和装置 |
| CN116642410A (zh) * | 2023-07-26 | 2023-08-25 | 西安交通大学 | 一种非接触式cfrp结构损伤监测系统及方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4468965A (en) * | 1982-11-05 | 1984-09-04 | Union Carbide Corporation | Test method for acetylene cylinders |
| CN101285796A (zh) * | 2008-04-30 | 2008-10-15 | 湘潭大学 | 热障涂层损伤及其失效过程的声发射实时检测方法 |
| CN102109498A (zh) * | 2009-12-28 | 2011-06-29 | 天津工业大学 | 一种对三维编制复合材料的无损检测系统及检测分析方法 |
| CN103852523A (zh) * | 2014-03-27 | 2014-06-11 | 河南省第二公路工程有限公司 | 基于声发射技术判别预应力混凝土梁结构损伤类型的系统及方法 |
| CN103995054A (zh) * | 2014-05-19 | 2014-08-20 | 中国矿业大学 | 一种基于声发射检测含水煤岩内部结构损伤分析的方法 |
-
2017
- 2017-12-13 CN CN201711332249.2A patent/CN108168996B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4468965A (en) * | 1982-11-05 | 1984-09-04 | Union Carbide Corporation | Test method for acetylene cylinders |
| CN101285796A (zh) * | 2008-04-30 | 2008-10-15 | 湘潭大学 | 热障涂层损伤及其失效过程的声发射实时检测方法 |
| CN102109498A (zh) * | 2009-12-28 | 2011-06-29 | 天津工业大学 | 一种对三维编制复合材料的无损检测系统及检测分析方法 |
| CN103852523A (zh) * | 2014-03-27 | 2014-06-11 | 河南省第二公路工程有限公司 | 基于声发射技术判别预应力混凝土梁结构损伤类型的系统及方法 |
| CN103995054A (zh) * | 2014-05-19 | 2014-08-20 | 中国矿业大学 | 一种基于声发射检测含水煤岩内部结构损伤分析的方法 |
Non-Patent Citations (1)
| Title |
|---|
| 史攀飞等: "《基于CT图像的混凝土材料声发射损伤检测优化》", 《江西建材》 * |
Cited By (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111351810A (zh) * | 2020-03-30 | 2020-06-30 | 山东省分析测试中心 | 一种含多尺度缺陷金属断裂失效行为的分析方法 |
| CN111678993A (zh) * | 2020-05-26 | 2020-09-18 | 中国电子产品可靠性与环境试验研究所（（工业和信息化部电子第五研究所）（中国赛宝实验室）） | 无损检测系统以及方法 |
| CN111678993B (zh) * | 2020-05-26 | 2024-01-05 | 中国电子产品可靠性与环境试验研究所（（工业和信息化部电子第五研究所）（中国赛宝实验室）） | 无损检测系统以及方法 |
| CN111999391A (zh) * | 2020-06-24 | 2020-11-27 | 天津科技大学 | 声发射监测装置及其对材料结构的疲劳损伤分析方法 |
| CN111707548A (zh) * | 2020-06-30 | 2020-09-25 | 中国矿业大学 | 一种基于裂隙和能量演化的锚杆支护效果评价方法 |
| CN111707548B (zh) * | 2020-06-30 | 2021-08-27 | 中国矿业大学 | 一种基于裂隙和能量演化的锚杆支护效果评价方法 |
| CN112800623A (zh) * | 2020-07-22 | 2021-05-14 | 北京理工大学 | 基于反应演化模型的带泄压结构弹药的反应烈度评估方法 |
| CN112800623B (zh) * | 2020-07-22 | 2022-07-05 | 北京理工大学 | 基于反应演化模型的带泄压结构弹药的反应烈度评估方法 |
| CN113063789A (zh) * | 2021-03-18 | 2021-07-02 | 山东建筑大学 | 一种盾构隧道管片内部微破坏检测装置及方法 |
| CN115043118A (zh) * | 2021-09-09 | 2022-09-13 | 海沃机械（中国）有限公司 | 用于垃圾转运车的拉臂部件的损伤识别方法和装置 |
| CN115043118B (zh) * | 2021-09-09 | 2023-10-13 | 海沃机械（中国）有限公司 | 用于垃圾转运车的拉臂部件的损伤识别方法和装置 |
| CN113917562B (zh) * | 2021-09-27 | 2023-02-28 | 中国科学院武汉岩土力学研究所 | 深埋软弱夹层宏细观结构表征及三维空间构造方法及装置 |
| CN113917562A (zh) * | 2021-09-27 | 2022-01-11 | 中国科学院武汉岩土力学研究所 | 深埋软弱夹层宏细观结构表征及三维空间构造方法及装置 |
| CN113984500A (zh) * | 2021-10-27 | 2022-01-28 | 中国航发北京航空材料研究院 | 一种材料失效监测装置及其分类方法 |
| CN114813798A (zh) * | 2022-05-18 | 2022-07-29 | 中国工程物理研究院化工材料研究所 | 用于表征材料内部结构及成分的ct检测装置和成像方法 |
| CN116642410A (zh) * | 2023-07-26 | 2023-08-25 | 西安交通大学 | 一种非接触式cfrp结构损伤监测系统及方法 |
| CN116642410B (zh) * | 2023-07-26 | 2023-10-10 | 西安交通大学 | 一种非接触式cfrp结构损伤监测系统及方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108168996B (zh) | 2019-03-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108168996B (zh) | 一种基于声发射信号的pbx损伤演化ct原位分析方法 | |
| Xing et al. | Full-field measurement and fracture characterisations of rocks under dynamic loads using high-speed three-dimensional digital image correlation | |
| Xie et al. | Fractal property of spatial distribution of acoustic emissions during the failure process of bedded rock salt | |
| Yang et al. | An experimental investigation on the damage of granite under uniaxial tension by using a digital image correlation method | |
| Loeffler et al. | Detection and segmentation of mechanical damage in concrete with X-Ray microtomography | |
| Tashman et al. | Microstructure characterization for modeling HMA behaviour using imaging technology | |
| Lo et al. | Microstructural and mechanical characterization of variability in porous advanced ceramics using X-ray computed tomography and digital image correlation | |
| Cala et al. | Identification of microstructural properties of shale by combined use of X-ray micro-CT and nanoindentation tests | |
| D’Accardi et al. | Quantitative analysis of thermographic data through different algorithms | |
| Wisner et al. | The use of nondestructive evaluation methods in fatigue: A review | |
| Daryabor et al. | Image fusion of ultrasonic and thermographic inspection of carbon/epoxy patches bonded to an aluminum plate | |
| Zhang et al. | Plastic deformation of submicron-sized crystals studied by in-situ Kikuchi diffraction and dislocation imaging | |
| Li et al. | Detection and characterization of mechanical impact damage within multi-layer carbon fiber reinforced polymer (CFRP) laminate using passive thermography | |
| Ai et al. | Experimental and numerical study on the fracture characteristics of concrete under uniaxial compression | |
| Nagy et al. | Volume and surface area distributions of cracks in concrete | |
| Hall et al. | Strain localisation in sand under triaxial loading: characterisation by x-ray micro tomography and 3D digital image correlation | |
| Ai et al. | Experimental study of fracture characterizations of rocks under dynamic tension test with image processing | |
| Kaur et al. | Physical insights into principal component thermography | |
| Wisner et al. | Fatigue damage precursor identification using nondestructive evaluation coupled with electron microscopy | |
| Luong et al. | Characterization of mechanical damage in granite | |
| Muc et al. | Prediction of subsurface defects through a pulse thermography: experiments vs. numerical modeling | |
| Viot et al. | Impact test deformations of polypropylene foam samples followed by microtomography | |
| Zuo et al. | Measuring micro/meso deformation field of geo-materials with SEM and digital image correlation method | |
| Blasone et al. | Investigation of Spalling Damage in Ultra-High Performance Concrete Through X-ray Computed Tomography | |
| Xiaodong et al. | Tensile strength of sea ice using splitting tests based on the digital image correlation method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |