Design of an Automatic Defect Identification Method Based ECPT for Pneumatic Pressure Equipment

,


Introduction
As one of the most significant energy storage and transmission kind of equipment, the pneumatic pressure equipment has been widely used in modern aerospace industries, such as manned spacecraft and ground wind tunnel testing in Figure 1.Moreover, given its wide range of applications, the equipment possesses a variety of specifications that make it flexible.The flexibility of their specifications adds complexity to problem of inspecting and testing this important equipment to identify defects.Researchers have studied methods to characterize pressure vessel properties in recent years.J. Lee et al. proposed a method of estimating the design pressure in a prismatic storage vessel for liquefied natural gas [1].H. Al-Gahtaniet et al. studied the local pressure testing of spherical vessels with nozzles [2].J. Proczka et al. reported guidelines for the efficient design and sizing of Small-Scale Compressed Air Energy Storage (SS-CAES) pressure vessels [3].P. Blanc studied the residual burst pressure of composite cylinders after mechanical impacts on their cylindrical part [4].
These studies, however, do not explore the impact of defects related to the on-going use of pressure vessels made of diverse materials in a wide range of operating environments.
The complicated usage environments and diversity of materials cause the pneumatic pressure equipment to generate fatigue cracks, corrosion pits, and other defects.Therefore it is important to detect the potential defects for the pressure equipment.Nondestructive testing and evaluation (NDT&E) [5][6][7][8][9][10] are important to ensure the safety of manufacturing environments and equipment operation.K. Schabowicz et al. presented some nondestructive methods for the testing of concrete members and fibre cement boards [11].Furthermore, the criteria for special specimen inspection are distinct, so various detection methods are needed including surface nondestructive testing, appearance inspection, and physicochemical test.The surface nondestructive testing mainly includes magnetic particle testing and penetrant testing.Due to the advantages of surface nondestructive testing, like high sensitivity and easy identification, this method been broadly adopted.However, there are many problems that need to be resolved in its practical application, like high labor intensity, low efficiency, pollution of the environment, poor safety, and so on.In recent years, ECPT has been developing rapidly.In addition to not damaging the body, the method is fast, is efficient, and can solve the issues of regular methods.So in this paper, the eddy current pulsed thermography is implemented in the detection of defects of pneumatic pressure equipment.Researchers also have made valuable contributions in processing data [12][13][14][15], which can be applied in (NDT&E).P. Zhu et al. presented an improved feature extraction algorithm for automatic defect identification based on eddy current pulsed thermography [16].H. Yu utilized the normalized fuzzy weighting functions to dispose the records of experiment [17].
Nevertheless, the above-mentioned researches possess low efficiency and limitations in processing data from lowalloy steel material.Besides, there are many methods, such as Independent Component Analysis (ICA); the processed feature information will be distorted after restructuring, because of data normalization, lending to no visual effect in detection.In recent years, adaptive technique has gained a lot of interest to enhance the efficiency of processing data [18][19][20][21].To make the detection system be more useful for detecting the defect, it may own the ability of adapting property.A novel characteristic identification algorithm that uses the similarity of the typical transient thermal responses () with the mixing vectors (the vectors of the pseudo-inverse matrix of the demixing matrix) in ICA is proposed.It solves the poor accuracy and low efficiency of the ICA efficiently by using the choice of a known message.The proposed algorithm has four steps.First the thermal imager collects the thermal image sequence of the pneumatic pressure equipment.Second, the  are separated into parts using threshold values with a variable interval search.The variable interval search in an infrared image sequence can reduce the repetitive computation and hold the typical feature of the pneumatic pressure equipment.Third, the membership matrix and the clustering center are calculated to classify the acquired .Moreover, the rule of the largest distance of two class is used to identify the typical .Fourth, the typical  are transformed into a two-dimension matrix in linear time.The main features of the infrared image sequence in the pneumatic pressure equipment can be extracted by the typical .Experimental results indicate the algorithm can select the typical feature more precise than ones of the ICA.Furthermore, the processing time of the presented algorithm is shorter than ones of the ICA.Therefore, the new algorithm can extract the main image information of the pneumatic pressure equipment accurately and efficiently in comparison to the ICA.

Background: Pneumatic Pressure
Equipment and Nondestructive Testing . .Pneumatic Pressure Equipment.The manufacture of high pressure vessels is complex, involving the coordination of subject knowledge and professional technology in many industries, including metallurgy, corrosion and protection, mechanical processing, safety protection, chemical engineering, and testing.The recent advances in technologies across various industries, including nondestructive testing, have enabled considerable progress in pressure vessel manufacturing.
(1) Varieties and specifications: the types of material can be mainly divided into carbon steel, low-alloy steel, and a small amount of stainless steel.Moreover, the main specification of external diameter has dozens of kinds.
(2) Complex usage environment: pressure vessels and pipes often suffer from high pressure, some special aerodynamic erosion, and environmental corrosion in the open air, which usually lead to cracks, perforation, fatigue damage, and corrosion of the equipment.
(3) Devastating consequences of defects: accidents not only endanger the safety of personnel, equipment, and the plant, but also lead to delay of test task, failures, and so on.
. .Nondestructive Testing.Surface nondestructive testing is the key test method of pressure equipment, such as pressure piping and pressure vessels.Surface nondestructive testing basically includes magnetic particle testing and penetrant testing.The method possesses numerous advantages such as high defect detection rate and high sensitivity.But there are also many problems, including high labor intensity, long overhaul cycle, low efficiency, pollution of the environment, and poor safety.So there exist a great many obstructors for development of detection skill.In flammable explosive environment, because of the safety risks, the traditional types of surface nondestructive testing methods (magnetic particle, penetrant) cannot be implemented in the field.
In recent years, surface detection by eddy current has been (developing) rapid.In addition to not damaging the body, the method is fast and efficient and can achieve large area rapid detection reducing manpower and material resources.The surface detection by eddy current is a method that utilizes thermal infrared image analysis.
The temperature field distribution is based on a vortex phenomenon that can be measured using infrared detection.The device is a thermal infrared imager with high speed and high resolution.Moreover, through the analysis of the infrared thermal image sequence, we could detect the changes of electromagnetic and thermal characteristics of the structure and material defects.It can detect the surface and near surface by using the eddy current effect; it also can detect deeper defects.In the light of Joule's law, the eddy current is turned into Joule heat in the testing piece generating high temperature and low temperature zones.The infrared thermal imager collects these temperature changes in the specimen.The collected data are analyzed to compute anticorrelation results and detect the presence of defects.
Eddy density can be indicated as   =  × (/), in which  expresses the magnetic vector bit and  indicates the electrical conductivity of materials.In the depth direction, the eddy current density is   () =   (0) ⋅  −√ , in which  is the depth,  is the permeability of testing material, and  is the frequency of alternating current flowing through an exciting coil.The depth of eddy density has an attenuation to the 1/ of the strength surface, which is called as skin depth  = 1/√.According to Joule's law, the heat power is 2 and E is electric field intensity.The Joule heat Q is transmitted inside the material, and the heat propagation equation is   (/) − Δ(  ∇) = , in which  is material density,   is specific heat capacity of materials,   is thermal conductivity of materials, and  is temperature of materials.The crack depth is larger than eddy current skin depth in this specimen, so in which  is height, H is magnetic field intensity, and The thermal power in the unit volume of the two sides of the crack is   = 0.5(1/)∫  2   , when the eddy skin depth is smaller than the plate thickness ( ≪ ), the thermal power   =  2 /(2).Recognizing the defects change the heat distribution in the parts, or zones of a specimen, the principal component analysis enhances the zones where the thermal response and the eigenvector has similar trends.So, in the enhancement of defects in a principal component, the amplitude of other regions will be inhibited and distributed in a very small range.In this paper, a new algorithm in view of principal component selection is presented to automatically extract the principal components containing defects.

Introduction of Proposed Algorithm in ECPT
First of all, some definitions are introduced in the Appendix to make the algorithm clear.Then, the proposed algorithm is explicit as follows.
Step .3D matrix  saves the initial image sequence.Every pixels value of thermal images is saved in the matrix .The third-dimension of  is the time axis t.
Step . .Set  ( = 1, 2, ⋅ ⋅ ⋅  ⋅ ⋅ ⋅ ) temperature thresholds (), ( = 1, 2, ⋅ ⋅ ⋅ ) from large to small, at the same time, the column which includes the highest temperature point (  ,   ,   ) is divided into  + 1 parts.After the row is divided, the temperature value of the column also should be cut apart rationally in order to extract the characteristics of the defect more accurately.
Step . .Set temperature thresholds    , for each part, and the highest temperature point   (  ,   , :), ( = 1, 2, ⋅ ⋅ ⋅ , ), will be found in every part.The number of each block is regarded as   , ( = 1, 2, ⋅ ⋅ ⋅ ).Moreover, the same method as that of column is employed to find the row variable interval   .The specific procedure of searching for interval is in the following.(d) Calculate the PCC of (  (  ,   , :) with   (  , , :)) (e) until their PCC is less than    .(f) The number of (  (  ,   , :) whose PCCs with   (  , , :)) are more than    is regarded as   .
(h) The same process as the searching for column interval to find the row interval.

Complexity (i)
The row interval is regarded as   .
Step . .What is more, set the threshold DD, initialize  = 1,  = 1,  = 1, and carry out the following procedures: Pneumatic pressure equipment includes high, medium, and low pressure storage containers and inlet/outlet gas pipelines, as shown in Figure 2. The pneumatic pressure equipment stores and conveys compressed air.The main features of these pressure vessels and the significance of their defect inspection are reflected in the following aspects.
The specific calculation process is shown as Figure 3.
Remark .For the interval   , the suitable method to seek out is searching for the length of region with largest temperature variation.What is more, the  with largest peak value is always contained in the region with the largest temperature variation.As analysis of Step 1, the length of region with largest temperature variation can be found by the coordinate value of   .For   , the homologous setting rule is analogous to   .The proposed algorithm can include the total typical temperature variations.And the representative temperature variations express the feature in the corresponding image sequence.characteristic pick-up algorithm is put forward to handle data.In the proposed algorithm, the appropriate interval values of column and row   , ( = 1, 2, ⋅ ⋅ ⋅ ),   , ( = 1, 2, ⋅ ⋅ ⋅ ) are set to reduce the repeated calculation.Moreover, the variable interval also can save the significant features.At the same time the proposed algorithm possesses less redundant computation than ICA.
Step .Divide (:, ) into  parts.The specific calculation process is shown as Figure 4.
(a) Set the cluster number L. And, at the same time, initialize the cluster center   ; moreover set the number of iterations , and set the weighting coefficient ; finally set terminating iterative threshold .(e) The membership maximization rule is used to deblur all characteristic information weight and get the category of each information weight.  = arg max(  (  )).
Remark .Making use of the method of computing membership function and updating the clustering center in algorithm can extract feature information accurately and effectively.Taking advantage of the COV to classify data is called hard division.Hard division divides every object into a certain category strictly, and each class is unrelated to each other.However, the actual defect information objective exists intermediation in form and category.Moreover there is no definite boundary conditions to distinguish the class.Therefore, classifying the categories according to the membership and cluster centers of each class can be more accurate.The feature information weights in each category represent the optimal weight of the feature information, making the extracted features more accurate and more reliable.
Step .The mean value of the  ℎ classification is Remark .‖ * ‖  is P-norm, ‖ * ‖ ∞ is infinity-norm, and ‖ * ‖ :  →  satisfies the following: (1) positivity: ‖‖ ≥ 0, and ‖‖ = 0 ⇐⇒  = 0(2) positive homogeneity: ‖‖ = ||‖‖(3) trigonometric inequality: Remark .The purpose of the proposed algorithm in Step 1 to Step 3 is to select the typical thermal responses.These typical responses have less intimately connection with each other in ECPT.The noncorrelation is larger and the thermal responses are more representative.Moreover, the selected accurate data is classified by Steps 2 and 3; the cluster center and membership function are carried out to classify the defect feature.The classification simplifies the data of image processing and makes the postprocessing more precise.Step .The 3 initial image sequence matrix  will be transformed into 2 matrix .Remark .The final purpose of the proposed algorithm is to extract the defect feature and ignore redundant information.The mean value of selected region is obtained by the similarity criterion of the method.Moreover the calculation of the uniform measurement divides the region into defect feature and other areas.
Remark .In order to highlight the efficiency and reasonability of the proposed method, the comparison of the ICA and the proposed algorithm is shown in Figure 6.In the first part "Data whitening", the data in the ICA is dealt with by normalization and then is computed to obtain the whitening vectors according to the corresponding feature vectors and feature matrix.These processes of calculation take a lot of time to dispose data and reduce efficiency of processing.However, the new algorithm does not have these computational procedures, to avoid the redundancy calculation.Hence, the amount of data in the proposed algorithm after processing of pretreatment and whitening is much less than that in the ICA, so, in the postprocessing, the proposed algorithm avoids redundancy calculation and possesses higher efficiency than the ICA.2) In the second part "Data analysis and comparison", compared with the extracted defect feature information graph, the amplitude of the graph extracted by the new algorithm contains the physical characteristics of the defect, and the content is closer to the defect feature information.Thus the proposed algorithm is more accurate and reasonable than the ICA.

Remark .
Step 1 realizes the segmentation of data block and the calculation of column interval and row interval.Step 2 and Figure 3 show the process of the variable interval search.
Step 3 is the correlation value classification.Step 4 selects the  typical  based on the distance of two classes and extracts the typical features.

Experimental Design
The setup of the experiment is shown in Figure 7.It includes five functional parts: induction heater; coil; PC; timing trigger; and IR camera.The induction heater produces high frequency alternating current for coil excitation.A rectangular coil applies directional excitation to heat the sample.The timing trigger controls the time for heating the sample.The IR camera records the thermal image sequence of the sample.The analysis of experimental data is represented below.In view of the devices volume and degree of damage, specimen 1 and specimen 2, shown in Figure 9, are extracted from the pneumatic pressure equipment to analyze.
The defects have been detected using the method of ECPT and the algorithm presented in Section 3. The two specimens are analyzed as follows.For sample 1, the parameter index of specimen 1 is presented in Table 1.
The process of Steps 3 and 4 is represented in  11.It is obvious that the features of the selected  are similar to the ICs.In Figure 12, the difference of  is shown.Comparing the proposed algorithm result with the result of ICA in Figure 13, the trends of the red curves are similar to the blue curves, respectively.Moreover, Figure 14 illustrates the precise and accurate defect feature processed by the algorithm.When comparing the proposed algorithm result with the result of the ICA in Figure 15, the final result of the algorithm ignores more redundant information than the ICA result.At the same time, the defect feature extracted by the proposed algorithm is more precise and accurate than the ICA result.Consequently, the algorithm has selected the typical thermal responses and extracted main features successfully.The proposed method can not only extract the main features like the ICA, but also reduce the processed time substantially.Its biggest advantage is the efficiency.Figure 16 shows the processing time of the ICA and the proposed algorithm.It is obvious that the proposed algorithm needs less time to complete the feature extraction process.For sample 2, the parameter index of specimen 2 is presented in Table 2.

CLAS could be chosen
The process of Steps 3 and 4 is represented in .9435,0.9678, and 0.9878, respectively.It is obvious that the features of the selected  are similar to the ICs.In Figure 19, the difference of  is shown.Comparing the proposed algorithm result with the result of ICA in Figure 20, the trends of the red curves are similar to the blue curves, respectively.Moreover, Figure 21 illustrates the precise and accurate defect feature processed by the algorithm.When comparing the proposed algorithm result with the result of the ICA in Figure 22, the final result of the algorithm ignores more redundant information than the ICA result.At the same time, the defect feature extracted by the proposed algorithm is more precise and accurate than the ICA result.Consequently, the algorithm has selected the typical thermal responses and extracted main features successfully.The proposed method can not only extract the main features like the ICA, but also reduce the processed time substantially.Its biggest advantage is the efficiency.Figure 23 shows the processing time of the ICA and the proposed algorithm.It is obvious that the proposed algorithm needs less time to complete the feature extraction process.

Conclusion and Future Work
In this paper, an accurate and more efficient algorithm in ECPT is proposed.The validity and efficiency of the proposed method are demonstrated with experimental results.The physical meaning of ECPT and the mathematical foundation of ICA are integrated in this proposed method.The results contribute to advancing the use of ECPT to detect defects as follows.
(1) The primary features of the thermal image sequences can be extracted by the proposed approach.Meanwhile the main features can be utilized to detect defects.
(2) Experimental results show that whitening preprocedure is time-consuming.The proposed algorithm has no data     whitening procedure.So the proposed algorithm is more effective.
(3) After processing of pretreatment and whitening, the amount of data in the proposed algorithm is much less than that of in ICA.Hence, in the postprocessing, the proposed algorithm avoids redundancy calculation and possesses higher efficiency than the ICA.So the proposed algorithm is more accurate and reasonable than the ICA.
Future work will pay more attention to how to enhance the efficiency and accuracy of the novel algorithm.In the experiment, some external interference will bring plenty of noise to influence the precision.Moreover, the change of threshold values in the algorithm is also anticipated to impact the efficiency of the detection.Therefore, it is worth researching how to reduce the external interference and to adjust the threshold value.

Appendix Notations of the Proposed Algorithm
To make the algorithm explicit and clear, the mathematics definitions in the algorithm are shown in Table 3.

CLAS could be chosen
So in the  th ,the  th CLAS could be chosen

Figure 2 :
Figure 2: Pneumatic pressure equipment in China Aerodynamics Research & Development Center.

Figure 4 :
Figure 4: The specific process of Step 2.

Figure 8
Figure 8 shows the damaged pneumatic pressure equipment.All the test materials in this paper are provided by the China Aerodynamics Research & Development Center, which broadly supports research on aerodynamic equipment.In view of the devices volume and degree of damage, specimen 1 and specimen 2, shown in Figure9, are extracted from the pneumatic pressure equipment to analyze.The defects have been detected using the method of ECPT and the algorithm presented in Section 3. The two specimens are analyzed as follows.For sample 1, the parameter index of specimen 1 is presented in Table1.The process of Steps 3 and 4 is represented in Figure 10.In Step 3, set p=1, b=1, c=100,  = 10 −5 .There are 46, 13, 35, and 2  in the corresponding parts.Through the algorithm of Step 3, in the 1  part, compared with other , the 2  possesses the maximum noncorrelation,   2 = 1.3124.The same as in the 2  part, the 11 ( ℎ) possesses the maximum noncorrelation,   11 2 = 1.3876.In the 3  part, the 35 ℎ possesses the maximum noncorrelation,   35 3 = 1.3097.In the 4 ℎ , the 1  possesses the maximum noncorrelation,   1 4 = 1.3951. 1 ,  2 ,  3 , and  4 have been extracted by Steps 3 and 4. The extraction results of proposed algorithm and the ICA result are shown in Figure11.It is obvious that the features of the selected  are similar to the ICs.In Figure12, the difference of  is shown.Comparing the proposed algorithm result with the result of ICA in Figure13, the trends of the red curves are similar to the blue curves, respectively.Moreover, Figure14illustrates the precise and accurate defect feature processed by the algorithm.When comparing the proposed algorithm result with the result of the ICA in Figure15, the final result of the algorithm ignores more redundant information than the ICA result.At the same time, the defect feature extracted by the proposed algorithm is more precise and accurate than the ICA result.Consequently, the algorithm has selected the typical thermal responses and extracted main features successfully.The proposed method can not only extract the main features like the ICA, but also reduce the processed time substantially.Its biggest advantage is the efficiency.Figure16shows the processing time of the ICA and the proposed algorithm.It is obvious that the proposed algorithm needs less time to complete the feature extraction process.

Figure 10 .
Figure 8 shows the damaged pneumatic pressure equipment.All the test materials in this paper are provided by the China Aerodynamics Research & Development Center, which broadly supports research on aerodynamic equipment.In view of the devices volume and degree of damage, specimen 1 and specimen 2, shown in Figure9, are extracted from the pneumatic pressure equipment to analyze.The defects have been detected using the method of ECPT and the algorithm presented in Section 3. The two specimens are analyzed as follows.For sample 1, the parameter index of specimen 1 is presented in Table1.The process of Steps 3 and 4 is represented in Figure 10.In Step 3, set p=1, b=1, c=100,  = 10 −5 .There are 46, 13, 35, and 2  in the corresponding parts.Through the algorithm of Step 3, in the 1  part, compared with other , the 2  possesses the maximum noncorrelation,   2 = 1.3124.The same as in the 2  part, the 11 ( ℎ) possesses the maximum noncorrelation,   11 2 = 1.3876.In the 3  part, the 35 ℎ possesses the maximum noncorrelation,   35 3 = 1.3097.In the 4 ℎ , the 1  possesses the maximum noncorrelation,   1 4 = 1.3951. 1 ,  2 ,  3 , and  4 have been extracted by Steps 3 and 4. The extraction results of proposed algorithm and the ICA result are shown in Figure11.It is obvious that the features of the selected  are similar to the ICs.In Figure12, the difference of  is shown.Comparing the proposed algorithm result with the result of ICA in Figure13, the trends of the red curves are similar to the blue curves, respectively.Moreover, Figure14illustrates the precise and accurate defect feature processed by the algorithm.When comparing the proposed algorithm result with the result of the ICA in Figure15, the final result of the algorithm ignores more redundant information than the ICA result.At the same time, the defect feature extracted by the proposed algorithm is more precise and accurate than the ICA result.Consequently, the algorithm has selected the typical thermal responses and extracted main features successfully.The proposed method can not only extract the main features like the ICA, but also reduce the processed time substantially.Its biggest advantage is the efficiency.Figure16shows the processing time of the ICA and the proposed algorithm.It is obvious that the proposed algorithm needs less time to complete the feature extraction process.
 st ,the  th CLAS could be chosen So in the  rd ,the  th CLAS could be chosen So in the  th ,the  th CLAS could be chosen So in the  nd ,the  th

Figure 15 :
Figure 15: The specimen defect, the ICA result, and the final result.

Figure 16 :
Figure 16: The time comparison for specimen 1.
st part For the  nd part So in the  st ,the  th CLAS could be chosen So in the  nd ,the  th CLAS could be chosen For the  rd part For the  th part So in the  rd ,the  th

Table 3
the all number of columns in D  the total number of the images at the t axis (, , :) the transient thermal response (  ,   ,   ) the peak temperature point (  ,   ,   ) = [(, , )]    the temperature thresholds in searching for the variable interval      the temperature thresholds in searching for the variable interval     the number of each block in searching for the variable interval     the number of each block in searching for the variable interval