Next Article in Journal
Nonlinear Optical Limiting and Radiation Shielding Characteristics of Sm2O3 Doped Cadmium Sodium Lithium Borate Glasses
Next Article in Special Issue
Damage Characterization of Carbon Fiber Composite Pressure Vessels Based on Modal Acoustic Emission
Previous Article in Journal
Mortar Bond Strength: A Brief Literature Review, Tests for Analysis, New Research Needs and Initial Experiments
Previous Article in Special Issue
Testing the Dispersion of Nanoparticles in a Nanocomposite with an Ultra-Low Fill Content Using a Novel Non-Destructive Evaluation Technique
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Binarization Mechanism Evaluation for Water Ingress Detectability in Honeycomb Sandwich Structure Using Lock-In Thermography

1
Eco-Sustainable Energy Research Institute, Kongju National University, 1223-24 Cheonan-daero, Seobuk-gu, Cheonan-si 31080, Korea
2
Department of Mechanical Engineering, School of Engineering, Kathmandu University, Dhulikhel, Kathmandu 45200, Nepal
3
Department of Future Convergence Engineering, Kongju National University, 1223-24 Cheonan-daero, Seobuk-gu, Cheonan-si 31080, Korea
4
Division of Mechanical & Automotive Engineering, Kongju National University, 1223-24 Cheonan-daero, Seobuk-gu, Cheonan-si 31080, Korea
*
Author to whom correspondence should be addressed.
Materials 2022, 15(6), 2333; https://doi.org/10.3390/ma15062333
Submission received: 20 January 2022 / Revised: 14 March 2022 / Accepted: 19 March 2022 / Published: 21 March 2022
(This article belongs to the Special Issue Non-Destructive Evaluation of Composite Materials)

Abstract

:
The growing use of composite honeycomb structures in several industries including aircraft has demonstrated the need to develop effective and efficient non-destructive evaluation methods. In recent years, active thermography has attracted great interest as a reliable technology for non-destructive testing and evaluation of composite materials due to its advantages of non-contact, non-destructive, full-area coverage, high speed, qualitative, and quantitative testing. However, non-uniform heating, low spatial resolution, and ambient environmental noise make the detection and characterization of defects challenging. Therefore, in this study, lock-in thermography (LIT) was used to detect water ingress into an aircraft composite honeycomb sandwich structure, and the phase signals were binarized through the Otsu algorithm. A square composite honeycomb with dimensions of 210 mm × 210 mm along with 16 different defective areas of various sizes in groups filled with water by 25%, 50%, 75%, and 100% of the cell volume was considered. The sample was excited at multiple modulation frequencies (i.e., 1 Hz to 0.01 Hz). The results were compared in terms of phase contrast and CNR according to the modulation frequency. In addition, the detectability was analyzed by comparing the number of pixels of water ingress in the binarized image and the theoretical calculation.

1. Introduction

A carbon fiber reinforced plastic (CFRP) laminate having a strong anisotropy consists of a sandwich structure. The core of this structure is sandwiched between the outer skins on both sides. It offers excellent strength to weight ratios and is mainly used in the aviation industry, specifically in aircraft flooring, doors, wing flaps and rudders; automobiles; and building panels [1]. It is used as a honeycomb sandwich panel in which a CFRP skin is bonded to a honeycomb core [2,3,4].
Defects in the sandwich structure occur due to various causes. Aircraft are often damaged after the commencement of operations [5]. Honeycomb composites are hollow structural materials, making them sensitive to liquid intrusion. Water and hydraulic oil are common types of liquid intrusion that can be caused by poor sealing or skin damage [6].
Numerical analysis and non-destructive testing are applied as methods for preventing and maintaining critical structural defects. Many numerical methods have been proposed in previous studies to predict the nucleation and propagation of failure in honeycomb structures [7]. The purpose of these studies is to develop a phase field model that can predict both the defect propagation path and the mechanical response of fiber-reinforced composite laminates [8]. Alternatively, to detect the presence of water in the aircraft structure early before fatal damage due to this problem, there has been considerable interest in various non-destructive testing (NDT) methods in the field, and research has been conducted in this area [9]. Over the years, ultrasonic [10], radiological, and neutron radiography, helium mass spectrometer leak detection, and hot water leak testing techniques have proven to be effective in controlling and detecting the water ingress inside honeycomb structures [11]. However, the ultrasonic inspection takes a lot of inspection time due to the inspection for each point, and the productivity is very low. Additionally, radiation and neutron radiographic examination techniques are rarely used in the field due to complex problems such as size, double-sided installation, and strict radiation precautions. Hot water leak testing has the disadvantage of potentially allowing water to penetrate components, and helium mass spectrometry techniques tend to be expensive to process and significantly dependent on user experience. In this context, infrared thermal imaging (IRT) has very great advantages in that it is fast, non-contact, non-destructive, has a fast inspection speed, large inspection area, is safe, and is portable.
IRT collects and analyzes thermal information using an infrared camera, a non-contact thermal imaging device based on the fact that all objects with an absolute temperature of 0K emit infrared energy. The infrared energy detected by the infrared camera is converted and can be processed to create a thermal image. IRT is generally divided into two approaches: passive thermography and active thermography [12]. LIT, also known as modulated thermal image, is an advanced NDT technology in which an object being inspected is heated with modulated thermal energy, causing the surface temperature to change periodically [13,14]. Basic concept and theory for the LIT have been frequently mentioned through previous literature [15,16,17,18]. The application of IRT NDT technology in the field of aircraft was first proposed in the 1960s [19]. In the 1990s, pulsed thermography was applied to detect corrosion and detachment of B737 testbed aircraft [20,21]. Since then, various IRT methods and data processing algorithms have been developed and applied to the detection of water ingress on honeycomb sandwich composites. Previous studies have been conducted to quantify the water ingress and defects such as delamination, impact damage, and fatigue failure on sandwich composites and structures [22,23,24]. In addition, it has been reported that PT and LIT were compared and that the phase angle of LIT produced a higher signal-to-noise ratio than PT [25]. PT is more susceptible to non-uniform heating, surface emissivity changes, surface morphology, and noise than LIT. The eddy current thermography has been introduced as another IRT technique to detect and characterize defects in honeycomb sandwich panels [26]. Furthermore, vibration and IRT combined neural networks have also proposed a technique for identifying the type, location, and extent of damage on sandwich composites [27]. Recently, the LIT algorithm was applied to test the water ingress of honeycomb sandwich panels, and SNR improvement studies were conducted [14].
This study presents an experimental study of the process of detecting water ingress in a typical composite honeycomb sandwich panel using LIT. The main purpose of this study was to perform a quantitative evaluation of LIT in reference samples under different conditions with various sizes and water ingress rates. Typically, a four-point algorithm is used to process LIT experimental data. The second purpose is to process LIT experimental data using binarization and compare its performance with an Otsu algorithm in terms of water ingress content. This study also explored the influence of water ingress content as well as the modulation frequency obtained with the method above-mentioned. Finally, the relationship between phase contrast and contrast to noise ratio (CNR) and its effect on the binarized image is presented. The remainder of the paper is composed as follows. In Section 2, the basic concept of LIT is provided. Subsequently, the theory of the four-point method and Otsu algorithm is provided in Section 3. In Section 4, the details of the test sample and experimental process are described. Section 5 presents the experimental results, and the conclusions of this study are indicated in Section 6.

2. Lock-In Thermography (LIT)

LIT is typically applied to amplitude and phase measurements to evaluate underlying defects. The phase and amplitude information calculated by the processing of the thermal images recorded for each pixel are stored in the form of a 2D matrix, which is then converted to images known as phase and amplitude images. The amplitude image shows the total temperature rise of the system during the power cycle and the phase image provides the power to the device and the time delay between subsequent heating of the surfaces. LIT is a technique to obtain the amplitude and phase images for quantitatively analyzing the size and depth of defects.
The periodical heat transfer through the semi-infinite (planar) body resulting in the thermal wave can be expressed by Equation (1) [28,29].
T z , t = A ( z ) cos [ ω t ( z ) ] = T 0 e z μ cos   [ ω t 2 π z λ ]
where A(z) resembles the thermal amplitude; ω = 2 π f [rad/s] is expressed as an angular frequency; t is the time; ( z ) is the phase shift; T0 [°C] is an initial change in temperature due to a heat source; z [mm] is the defect depth; λ = 2 π μ [m] is expressed as the thermal wavelength; and μ [m] is the thermal diffusion length and can be calculated from the thermal diffusivity and the modulation frequency using Equation (2) [28,30].
μ = 2 α ω = α π f
where f [Hz] is the modulation frequency.
Equation (2) revealed that the thermal diffusion length depends on the modulation frequency and thermal diffusivity (i.e., thermophysical properties: thermal conductivity, density, and specific heat capacity) of the material.
Considering Equation (1), the phase shift ( z ) can be defined by Equation (3), which indicates that the phase is associated with defect depth (z) without the influence of other parameters [15].
= z μ
In contrast, amplitude A(z) can be defined by Equation (4), which indicated that amplitude is associated with the initial temperature T 0 and moreover to some other parameters [28].
A ( z ) = T 0 e z μ

3. Image Processing

3.1. Phase and Amplitude Computation

For the experimental LIT data, the four-point method and the Fourier transform computes the amplitude and phase angle data from the recorded temperature–time history of each pixel, then stored in the form of 2D matrices and converted to amplitude and phase images. Amplitude and phase images are related to the thermal diffusivity and the propagation time, respectively.

3.2. Four-Point Method

If S1, S2, S3, and S4 are four equidistant thermal images, as shown in Figure 1, in a constant period, then the phase ( ) and amplitude (A) images can be expressed by Equations (5) and (6), respectively [31,32].
= tan 1 ( S 1 S 3 S 2 S 4 )
A = ( S 1 S 3 ) 2 + ( S 2 S 4 ) 2
The four-point method is fast and valid for sinusoidal stimulation, but is affected by noise. The increased number of cycles results in the denoising of the signal.

3.3. Binarization Using the Otsu Algorithm

There are many techniques to analyze thermal images, but the simplest and easiest method among them is to use a threshold to binarize the image. Through binarization, automated detection can be characterized for various defects. This is used in the image pre-processing stage to separate the background from the image opening to only extract pixels with data (brightness, etc.) above a certain level, or to simplify the entire image information [33].
In the Otsu algorithm, the threshold is calculated based on the gray scale. An image is divided into two using an image histogram, and an optimal threshold value can be estimated. In the binary image, [0, k] is classified as ‘class 1’, and [k, 1] is classified as ‘class 2’ based on the threshold value k. Through this process, binary images can be acquired based on a threshold value. In general, converting a 2D thermal image to a binary image can characterize the characteristics of the defects [34,35].
In order to classify into two binary images, it is necessary to calculate an optimal threshold value. If it is an M × N image with L intensity levels such as 0, 1, 2, …, L − 1, pixels with intensity values within [0, k] are classified as class 1, and intensity values within [k + 1, L + 1] are classified as class 2. The probability that a pixel is classified into class 1 or 2 is defined by Equations (7) and (8) [36]:
P 1 ( k ) = i = 0 k P i
P 2 ( k ) = 1 P 1 ( k )
The average intensity values of pixels classified into classes 1 and 2 are defined as Equations (9) and (10).
m 1 ( k ) = 1 P 1 ( k ) i = 0 k iP i
m 2 ( k ) = 1 P 2 ( k ) i = k + 1 L 1 iP i
There are mean intensity values up to the k level, which for all images is
m G = P 1 m 1 + P 2 m 2
In order to calculate the optimal threshold value, the Otsu algorithm should allow the concept of between-class variance. The equation of between-class variance is defined as Equation (12).
σ b 2 = ( m G P 1 m ) 2 P 1 ( 1 P 1 )
Calculating the optimal k value is a simple principle, although the optimal k value can be calculated only by substituting all k values in the intensity range [0, L − 1].

4. Methods and Materials

4.1. Test Sample

The study was conducted by applying a honeycomb sandwich sample fabricated by the EZ composite in Korea. The sample was applied with an aramid honeycomb core sandwich panel with CFRP on the front and glass fiber reinforced plastic (GFRP) composite face skin on the back. A heat-resistant phenolic resin was coated on the surface of the specimen to improve the strength and thermal properties. The CFRP skin on the front side of the sample was made of SKYFLEX WSN3KY and Mitsubishi Rayon TR-30, and the GFRP skin on the back side of the sample was made of GEP7628 and E-glass. Figure 2 and Figure 3 show the schematic diagrams and geometrical information of the honeycomb sandwich samples considering various water ingress conditions. The sample is a square plate with a size of 210 mm × 210 mm. In the sample, 16 types of water ingress under various conditions from 10 mm × 10 mm to 40 mm× 40 mm were considered. In addition, the water ingress rate was classified into four groups (filled with water by 25%, 50%, 75%, and 100% of the cell volume.)

4.2. Experimentation

The experiment was conducted with lock-in thermography, which has light as a heat source from the outside. The surface of the sample was heated with a halogen lamp with a sinusoidal heat source. Thermal response from the surface due to the possession of a heat source was investigated with an infrared camera and temperature data were recorded. The heat source was controlled by power amplifiers and function generators and two lamps of 1 kW halogen were applied as external heat sources. The applied frequency was set by decreasing the frequency from 1 Hz to 0.01 Hz. The function generator (Agilent 33210A, Agilent Technologies, Santa Clara, CA, USA) was used to produce sinusoidal waves, and the infrared camera used to detect radiation from the surface of the specimen was applied with the FLIR SC645 model (640 × 480 pixels, 7.5 to 13 μm, Teledyne FLIR, Wilsonville, OR, USA). The sample can be placed on the infrared camera screen using a lens of 41.3 mm. Infrared images and temperature data were obtained using FLIR R&D commercial software. The frame of the camera was set at 50 frames per second for each frequency. The acquired data were processed through MATLAB 2021b software coding. A schematic of the general LIT test configuration is depicted in Figure 4.

5. Results and Discussions

5.1. Water Ingress Detectability in Phase Image

In this study, only the phase angle images obtained from the four-point method were considered in the analysis because phase angle images were relatively free from non-uniform heating, optical absorption, and emissivity variation when compared with the amplitude images. A series of thermal imaging experiments were conducted while reducing the modulation frequency from 1 Hz to 0.01 Hz. In the experiment, a total of two cycles of modulation frequency was applied, and all pixels of the thermal images were processed with phase through Equation (5). Figure 5 shows the resulting images of the phase images at 1 Hz, 0.5 Hz, 0.2 Hz, 0.1 Hz, 0.07 Hz, 0.05 Hz, 0.03 Hz, and 0.01 Hz.
This study shows that the possibility of detecting water ingress in the phase images depends on the applied modulation frequency. As shown in Figure 6, it can be seen that no water ingress was detected at the modulation frequency phase of 1 Hz. It can be seen that when the modulation frequency was reduced to 0.5 Hz, the contrast with the background was improved, and all water ingress began to be detected. In addition, as the modulation frequency was reduced to 0.03 Hz or less, the image resolution was lowered by moisture present therein, and the sharpness of defects deteriorated. This is due to strong external energy due to the different thermal conductivity between water and composites in the surface transient heat transfer mode.

5.2. Phase Contrast and CNR Trend

In LIT, it is important to analyze the effect of the modulation frequency of an external heating source. Therefore, in this study, the optimal frequency was determined through the phase contrast and noise contrast trends, and their relationship was analyzed. First, in order to determine the optimal modulation frequency, the contrast between the water ingress and the adjacent sound area was analyzed, and the contrast analysis is defined as Equation (13) [14].
Δ C = C d C s
where C d is the phase angle of defective area (water ingress), and C s is the phase angle of the sound area.
D1, A1, A3, and D3 with fixed dimensions of 40 mm × 40 mm and water ingress rates of 100%, 75%, 50%, and 25%, respectively, were selected for analysis. As shown in Figure 6, for one, 20 × 20 pixels were applied for the water ingress and sound areas with dimensions of 20 mm × 20 mm or more, and 10 × 10 pixels were applied for the rest.
Figure 7 and Figure 8 show the phase contrast and CNR trends of D1, A1, A3, and D3, and Table 1 shows the total data of phase contrast and CNR. Figure 7 shows the phase contrast according to the modulation frequency. It can be seen that the phase contrast increases from 0.01 Hz to 0.1 Hz as the frequency increases, and then decreases again after 0.1 Hz. Therefore, phase detection can be estimated as the optimal frequency at 0.1 Hz. The maximum phase contrast of modulation frequency in 0.1 Hz, D1 was 0.3502 rad, A1 was 0.2247 rad, A3 was 0.1641 rad, and D4 was 0.1288 rad.
Second, to analyze the relationship between phase contrast and CNR, values according to the modulation frequencies were compared. CNR was the same as a signal to noise ratio (SNR) used in the electrical signal description. The CNR metric is used to differentiate the two regions of interest (ROI). CNR is the ratio of the size of the signal (defects) to the size of the noise (sound). Image processing is performed to increase the SNR, and the higher the CNR value, the better the detection ability of the fault. CNR is defined as Equation (14) [37,38].
CNR = 20 log 10 ( | DROI mean SROI mean | σ )
where DROImean is the arithmetic mean of the defective area; SROImean is the arithmetic mean of the sound area; and σ is the standard deviation of the defective area.
Figure 8 shows the CNR trend according to the modulation frequency of D1, A1, A3, and D3. CNR also showed a similar trend to phase contrast. CNR also showed a similar trend to phase contrast. The best contrast was shown at a modulation wave frequency of 0.1 Hz, and the maximum CNR value was 31.4413 dB for D1, 22.8425 dB for A1, 22.3637 dB for A3, and 21.9970 dB for D3.
Comparing the relationship between phase contrast and CNR, as shown in Table 1, it can be seen that the best CNR was obtained where the maximum phase contrast occurred (i.e., the better the phase contrast, the higher the CNR is obtained). Additionally, the closer the phase contrast is to blind frequency (0 rad), the poorer the CNR quality, and it may have a negative value. For example, the phase contrast values of the modulation frequency 0.07 Hz with negative CNR values of D1, A1, A3, and D3 were −0.0226 rad, 0.0039 rad, 0.0197 rad, and 0.0152 rad, respectively. Furthermore, at 1 Hz, they were 0.0433 rad, 0.0473 rad, −0.0089 rad, and −0.0301 rad, respectively.

5.3. Binary Image Evaluation

For water ingress detection of the phase angle images, binarization processing using the Otsu algorithm was performed on images to which median filtering was applied. Figure 9 shows the results of image binarization through the Otsu algorithm processing of Equations (7)–(12). The binarization was processed for 0.2 Hz and 0.05 Hz, which showed relatively good defect detection potential including 0.1 Hz, which showed the maximum phase contrast. The optimum threshold value for the phase angle for binarization was 4. As shown in Figure 9, all water ingress was detected at 0.05 Hz, 0.1 Hz, and 0.2 Hz. However, unlike 0.1 Hz, at 0.05 Hz and 0.2 Hz, the shape of the defect was non-uniform and very different from that of a square. It depends on the phase contrast and CNR.
The size can be estimated by the number of pixels occupied by the defect in the image. In general, to measure the size of a defect, an actual length corresponding to a unit pixel was first obtained according to the angle of view of an infrared camera using an indicator. Next, the size of the defect was calculated by multiplying the length corresponding to the unit pixel by the number of pixels occupied by the defect. The size of the defect D [ mm ] can be expressed as Equation (15) [39].
D = M × L P
where L is the length of the indicator; P is the number of fractions corresponding to the length of the indicator; and M is the number of pixels occupied by the defect.
In order to quantitatively evaluate the binarized image, the approximate error was analyzed by comparing the number of pixels occupied by each water ingress in the image and the numerical value calculated by Equation (15). As shown in Figure 10, the number of pixels in the image was measured by counting the number of pixels representing level 1 in each area. Table 2 shows the number of pixels in the image for all water ingress and the number of pixels by theoretical calculation. For example, D1 with a water ingress rate of 100% and a size of 40   mm × 40 mm was counted with 7333 pixels in the binarized image. In an actual 210 mm sample, 40 mm becomes 80 mm in a two-dimensional image of 420 × 420 pixels, so theoretically, it becomes 6400 pixels. Figure 11 shows the change in the measurement error of 40   mm × 40 mm and 20   mm × 20 mm of D1 and A3, according to the water ingress rate. In Figure 11, it can be seen that the higher the water ingress rate, the negative error changes to a positive error. This is because the higher the water ingress rate, the more active the excessive heat conduction phenomenon as external energy heating occurs on the surface and the difference in thermal conductivity between water and CFRP. Conversely, the smaller the size and the lower the water ingress rate, less heat energy is supplied and heat exchange does not occur, so the defect detection is not perfect, resulting in a negative error. For example, D3 with a water ingress rate of 25% had a negative error because 5078 pixels were counted in the image, 1322 pixels less than the 6400 pixels in the theoretical calculation. Moreover, in D1 of 100% water ingress rate, 7333 pixels were counted in the image, and there were 933 pixels more than the 6400 pixels in the theoretical calculation, so it became a positive error.

6. Conclusions and Future Works

This study conducted LIT experimental tests to examine the water ingress of composite honeycomb sandwich panels containing various sizes by the water ingress rate levels. Non-destructive testing by lock-in thermography reduced the effect of non-uniform heating. The main observations of this study are as follows. First, all water ingress defects under various conditions were detected by the LIT technique, and the performance of detectability was excellent. Second, the concept of LIT data processing and the Otsu algorithm for binarization were reviewed. The optimum frequency was derived by analyzing the relationship between the phase contrast and CNR according to the modulation frequency; a negative CNR value was obtained as the phase contrast was closer to 0 rad. Third, the binarization by the Otsu algorithm was evaluated for detectability by simplifying the defect shape. As a result, the binarization image by the Otsu algorithm had good detectability of water penetration, but it was greatly affected by the thermal conductivity of the material. Finally, due to the significant difference in thermal conductivity between two materials, caused by the active heat transfer on the surface, and the stronger the heat source, the greater the effect.
Through this study, it is proposed that the water ingress inside the composite material can be detected using lock-in thermography. In future studies, we will consider advanced noise removal techniques based on thermal properties. In addition, data and image processing studies will be considered for predictive simulation and experimental studies of non-linear defects such as damage and cracks. Furthermore, extension studies of LIT such as deep learning-based automated defect recognition, classification, and detection enhancement will be performed.

Author Contributions

Conceptualization, W.K., Y.C. and R.S.; Data curation, R.S.; Formal analysis, R.S. and S.L.; Funding acquisition, W.K.; Investigation, W.K., Y.C., R.S. and S.L.; Methodology, Y.C., R.S. and S.L.; Project administration, W.K.; Resources, W.K.; Software, Y.C., R.S. and S.L.; Supervision, W.K.; Validation, W.K. and R.S.; Visualization, R.S.; Writing—original draft, Y.C. and R.S.; Writing—review & editing, S.L. and W.K. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Research Foundation of Korea (NRF–2022R1A2C1011413) funded by the Korean Government, Ministry of Education, Science and Technology (MEST).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Due to the nature of this research, participants of this study did not agree for their data to be shared publicly and are only available upon reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ciampa, F.; Mahmoodi, P.; Pinto, F.; Meo, M. Recent Advances in Active Infrared Thermography for Non-Destructive Testing of Aerospace Components. Sensors 2018, 18, 609. [Google Scholar] [CrossRef] [Green Version]
  2. Vinson, J. The Behavior of Sandwich Structures of Isotropic and Composite Materials; Routledge: London, UK, 2018. [Google Scholar]
  3. Vinson, J.R. Sandwich structures: Past, present, and future. In Sandwich Structures 7: Advancing with Sandwich Structures and Materials; Springer: Dordrecht, The Netherlands, 2005; pp. 3–12. [Google Scholar]
  4. Wang, Z. Recent Advances in Novel Metallic Honeycomb Structure. Compos. Part B Eng. 2019, 166, 731–741. [Google Scholar] [CrossRef]
  5. Sfarra, S.; Ibarra-Castanedo, C.; Avdelidis, N.P.; Genest, M.; Bouchagier, L.; Kourousis, D.; Tsimogiannis, A.; Anastassopoulous, A.; Bendada, A.; Maldague, X.; et al. A Comparative Investigation for the Nondestructive Testing of Honeycomb Structures by Holographic Interferometry and Infrared Thermography. J. Phys. Conf. Ser. 2010, 214, 012071. [Google Scholar] [CrossRef] [Green Version]
  6. Roy, R.; Nguyen, K.H.; Park, Y.B.; Kweon, J.H.; Choi, J.H. Testing and Modeling of Nomex™ Honeycomb Sandwich Panels with Bolt Insert. Compos. Part B Eng. 2014, 56, 762–769. [Google Scholar] [CrossRef]
  7. Bui, T.Q.; Hu, X. A Review of Phase-Field Models, Fundamentals and their Applications to Composite Laminates. Eng. Fract. Mech. 2021, 248, 107705. [Google Scholar] [CrossRef]
  8. Zhang, P.; Hu, X.; Bui, T.Q.; Yao, W. Phase Field Modeling of Fracture in Fiber Reinforced Composite Laminate. Int. J. Mech. Sci. 2019, 161–162, 105008. [Google Scholar] [CrossRef]
  9. Birman, V.; Kardomateas, G.A. Review of Current Trends in Research and Applications of Sandwich Structures. Compos. Part B Eng. 2018, 142, 221–240. [Google Scholar] [CrossRef]
  10. Hsu, D.K. Non-destructive evaluation (NDE) of aerospace composites: Ultrasonic techniques. In Non-Destructive Evaluation (NDE) of Polymer Matrix Composites; Elsevier: Amsterdam, The Netherlands, 2013; pp. 397–422. [Google Scholar]
  11. Giguere, J.S.R. Damage Mechanisms and Non-Destructive Testing in the Case of Water Ingress in CF-18 Flight Control Surfaces Defence R&D Canada; Defense and Civil Institute of Evironmental Medicine: Toronto, ON, Canada, 2000. [Google Scholar]
  12. Chung, Y.; Shrestha, R.; Lee, S.; Kim, W. Thermographic Inspection of Internal Defects in Steel Structures: Analysis of Signal Processing Techniques in Pulsed Thermography. Sensors 2020, 20, 6015. [Google Scholar] [CrossRef]
  13. Ranjit, S.; Choi, M.; Kim, W. Quantification of Defects Depth in Glass Fiber Reinforced Plastic Plate by Infrared Lock-in Thermography. J. Mech. Sci. Technol. 2016, 30, 1111–1118. [Google Scholar] [CrossRef]
  14. Shrestha, R.; Choi, M.; Kim, W. Thermographic Inspection of Water Ingress in Composite Honeycomb Sandwich Structure: A Quantitative Comparison among Lock-in Thermography Algorithms. Quant. Infrared Thermogr. 2021, 18, 92–107. [Google Scholar] [CrossRef]
  15. Vavilov, V.P.; Taylor, R. Theoretical and Practical Aspects of the Thermal Nondestructive Testing of Bonded Structures. In Research Techniques in Nondestructive Testing; Academic Press: Cambridge, MA, USA, 1982; Volume 5, pp. 239–279. [Google Scholar]
  16. Lahiri, B.B.; Bagavathiappan, S.; Reshmi, P.R.; Philip, J.; Jayakumar, T.; Raj, B. Quantification of Defects in Composites and Rubber Materials using Active Thermography. Infrared Phys. Technol. 2012, 55, 191–199. [Google Scholar] [CrossRef]
  17. Wu, D.; Busse, G. Lock-in Thermography for Nondestructive Evaluation of Materials. Rev. Générale De Therm. 1998, 37, 693–703. [Google Scholar] [CrossRef]
  18. Pickering, S.; Almond, D. Matched Excitation Energy Comparison of the Pulse and Lock-in Thermography NDE Techniques. NDT E Int. 2008, 41, 501–509. [Google Scholar] [CrossRef]
  19. Kutzscher, E.W.; Zimmermann, K.H. A Scanning Infrared Inspection System Applied to Nondestructive Testing of Bonded Aerospace Structures. Appl. Opt. 1968, 7, 1715–1720. [Google Scholar] [CrossRef]
  20. Favro, L.D.; Ahmed, T.; Han, X.; Wang, L.; Wang, X.; Kuo, P.K.; Thomas, R.L. Thermal wave imaging of disbonding and corrosion on aircraft. In Review of Progress in Quantitative Nondestructive Evaluation; Springer: Boston, MA, USA, 1996; pp. 1747–1753. [Google Scholar]
  21. Thomas, R.L.; Han, X.; Favro, L.D. Thermal Wave Imaging of Aircraft for Evaluation of Disbonding and Corrosion. In Proceedings of the Ecndt’ 98, 7th European Conference on Non-Destructive Testing, Copenhagen, Denmark, 26–29 May 1998; pp. 126–130. [Google Scholar]
  22. Meola, C.; Carlomagno, G.M.; Squillace, A.; Vitiello, A. Non-Destructive Evaluation of Aerospace Materials with Lock-in Thermography. Eng. Fail. Anal. 2006, 13, 380–388. [Google Scholar] [CrossRef]
  23. Vavilov, V.P.; Klimov, A.G.; Nesteruk, D.; Shiryaev, V.V. Detecting Water in Aviation Honeycomb Structures by using Transient Infrared Thermographic NDT. In Proceedings of the Thermosense Xxv, AEROSENSE 2003, Orlando, FL, USA, 21–25 April 2003; SPIE: Bellingham, WA, USA, 2003; pp. 345–355. [Google Scholar]
  24. Vavilov, V.; Marinetti, S.; Pan, Y.; Chulkov, A. Detecting Water Ingress in Aviation Honeycomb Panels: Qualitative and Quantitative Aspects. Polym. Test. 2016, 54, 270–280. [Google Scholar] [CrossRef]
  25. Ibarra-Castanedo, C.; Piau, J.; Guilbert, S.; Avdelidis, N.P.; Genest, M.; Bendada, A.; Maldague, X.P. Comparative Study of Active Thermography Techniques for the Nondestructive Evaluation of Honeycomb Structures. Res. Nondestr. Eval. 2009, 20, 1–31. [Google Scholar] [CrossRef]
  26. He, Y.; Tian, G.; Pan, M.; Chen, D. Non-Destructive Testing of Low-Energy Impact in CFRP Laminates and Interior Defects in Honeycomb Sandwich using Scanning Pulsed Eddy Current. Compos. Part B Eng. 2014, 59, 196–203. [Google Scholar] [CrossRef] [Green Version]
  27. Just-Agosto, F.; Serrano, D.; Shafiq, B.; Cecchini, A. Neural Network Based Nondestructive Evaluation of Sandwich Composites. Compos. Part B Eng. 2008, 39, 217–225. [Google Scholar] [CrossRef]
  28. Meola, C.; Carlomagno, G.M. Recent Advances in the use of Infrared Thermography. Meas. Sci. Technol. 2004, 15, R27. [Google Scholar] [CrossRef]
  29. Montanini, R. Quantitative Determination of Subsurface Defects in a Reference Specimen made of Plexiglas by Means of Lock-in and Pulse Phase Infrared Thermography. Infrared Phys. Technol. 2010, 53, 363–371. [Google Scholar] [CrossRef]
  30. Fedala, Y.; Streza, M.; Sepulveda, F.; Roger, J.-P.; Tessier, G.; Boué, C. Infrared Lock-in Thermography Crack Localization on Metallic Surfaces for Industrial Diagnosis. J. Nondestruct. Eval. 2013, 33, 335–341. [Google Scholar] [CrossRef]
  31. Maldague, X. Theory and Practice of Infrared Technology for Nondestructive Testing; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2001. [Google Scholar]
  32. Czichos, H. Handbook of Technical Diagnostics: Fundamentals and Application to Structures and Systems; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2013. [Google Scholar]
  33. Lee, S.J.; Kim, W.T. Binarization Mechanism for Detectability Enhancement using Lock-in Infrared Thermography. J. Korean Soc. Nondestruct. Test. 2021, 41, 79–86. [Google Scholar] [CrossRef]
  34. Manda, M.P.; Kim, H.S. A Fast Image Thresholding Algorithm for Infrared Images Based on Histogram Approximation and Circuit Theory. Algorithms 2020, 13, 207. [Google Scholar] [CrossRef]
  35. Dong, Y.X. Review of Otsu Segmentation Algorithm. In Advanced Materials Research; Trans Tech Publications Ltd.: Bäch, Switzerland, 2014; pp. 1959–1961. [Google Scholar]
  36. Talab, A.M.A.; Huang, Z.; Xi, F.; HaiMing, L. Detection Crack in Image using Otsu Method and Multiple Filtering in Image Processing Techniques. Optik 2016, 127, 1030–1033. [Google Scholar] [CrossRef]
  37. Švantner, M.; Muzika, L.; Chmelík, T.; Skála, J. Quantitative Evaluation of Active Thermography using Contrast-to-Noise Ratio. Appl. Opt. Opt. Technol. Biomed. Opt. 2018, 57, D49–D55. [Google Scholar] [CrossRef]
  38. Netzelmann, U.; Müller, D. Modified Pulse-Phase Thermography Algorithms for Improved Contrast-to-Noise Ratio from Pulse-Excited Thermographic Sequences. NDT Amp E Int. Indep. Nondestruct. Test. Eval. 2020, 116, 102325. [Google Scholar] [CrossRef]
  39. Ranjit, S.; Kang, K.; Kim, W. Investigation of Lock-in Infrared Thermography for Evaluation of Subsurface Defects Size and Depth. Int. J. Precis. Eng. Manuf 2015, 16, 2255–2264. [Google Scholar] [CrossRef]
Figure 1. Principle of the four-point method in lock-in thermography.
Figure 1. Principle of the four-point method in lock-in thermography.
Materials 15 02333 g001
Figure 2. Schematic of the honeycomb-cored sandwich panel considering various water ingress conditions.
Figure 2. Schematic of the honeycomb-cored sandwich panel considering various water ingress conditions.
Materials 15 02333 g002
Figure 3. Aramid honeycomb-cored sandwich panel as a test sample. (a) The front face skin with CFRP and (b) the rear face skin with GFRP.
Figure 3. Aramid honeycomb-cored sandwich panel as a test sample. (a) The front face skin with CFRP and (b) the rear face skin with GFRP.
Materials 15 02333 g003
Figure 4. Experimental configuration of a lock-in thermography inspection system.
Figure 4. Experimental configuration of a lock-in thermography inspection system.
Materials 15 02333 g004
Figure 5. Phase images by the four-point method: (a) 1 Hz, (b) 0.5 Hz, (c) 0.2 Hz, (d) 0.1 Hz, (e) 0.07 Hz, (f) 0.05 Hz, (g) 0.03 Hz, (h) 0.01 Hz.
Figure 5. Phase images by the four-point method: (a) 1 Hz, (b) 0.5 Hz, (c) 0.2 Hz, (d) 0.1 Hz, (e) 0.07 Hz, (f) 0.05 Hz, (g) 0.03 Hz, (h) 0.01 Hz.
Materials 15 02333 g005
Figure 6. Representation of selected areas for contrast and CNR computation. The black dashed line enclosed area represents the water ingress areas and the red solid line enclosed area represents the sound areas.
Figure 6. Representation of selected areas for contrast and CNR computation. The black dashed line enclosed area represents the water ingress areas and the red solid line enclosed area represents the sound areas.
Materials 15 02333 g006
Figure 7. Phase contrast with the four-point method as a function of modulation frequency.
Figure 7. Phase contrast with the four-point method as a function of modulation frequency.
Materials 15 02333 g007
Figure 8. Phase-contrast with four-point method as a function of modulation frequency.
Figure 8. Phase-contrast with four-point method as a function of modulation frequency.
Materials 15 02333 g008
Figure 9. Results of binarized images through Otsu algorithm processing: (a) 0.05 Hz, (b) 0.1 Hz, (c) 0.2 Hz.
Figure 9. Results of binarized images through Otsu algorithm processing: (a) 0.05 Hz, (b) 0.1 Hz, (c) 0.2 Hz.
Materials 15 02333 g009
Figure 10. A binarized image at a modulation frequency of 0.1 Hz to calculate the number of pixels in a region.
Figure 10. A binarized image at a modulation frequency of 0.1 Hz to calculate the number of pixels in a region.
Materials 15 02333 g010
Figure 11. Error changes according to water ingress rate.
Figure 11. Error changes according to water ingress rate.
Materials 15 02333 g011
Table 1. Quantitative data of the phase contrast and CNR according to the modulation frequency.
Table 1. Quantitative data of the phase contrast and CNR according to the modulation frequency.
Modulation Frequency (Hz)
10.50.20.10.070.050.030.01
Phase contrast
(rad)
D10.04330.06260.17600.3502−0.02260.08310.1368−0.1028
A10.04730.04720.13010.22470.00390.05790.03610.0800
A3−0.00890.04080.09690.16410.01970.0452−0.0389−0.0547
D3−0.03010.01840.07050.12880.01520.0333−0.0433−0.0905
CNR
(dB)
D1−17.16406.806921.986831.4413−1.090822.653431.355621.3124
A1−19.81380.734415.506522.84250.890517.446620.38122.638
A3−37.06971.412114.820522.3637−6.416317.507717.440719.6345
D3−10.5926−2.588913.694521.997−20.377315.854315.85817.8589
Table 2. Estimation of the number of pixels and error for water ingress.
Table 2. Estimation of the number of pixels and error for water ingress.
Water Ingress No.PixelsError (%)
ImageTheory
A17059640010.30
A215911600−0.56
A3699164009.23
A415051600−5.93
B14400360022.22
B23019250020.76
B3364036001.11
B424272500−2.92
C1108690020.67
C2362400−9.5
C3753900−16.33
C4324400−19.00
D17333640014.58
D2175616009.75
D350786400−22.22
D46321600−47.94
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chung, Y.; Shrestha, R.; Lee, S.; Kim, W. Binarization Mechanism Evaluation for Water Ingress Detectability in Honeycomb Sandwich Structure Using Lock-In Thermography. Materials 2022, 15, 2333. https://doi.org/10.3390/ma15062333

AMA Style

Chung Y, Shrestha R, Lee S, Kim W. Binarization Mechanism Evaluation for Water Ingress Detectability in Honeycomb Sandwich Structure Using Lock-In Thermography. Materials. 2022; 15(6):2333. https://doi.org/10.3390/ma15062333

Chicago/Turabian Style

Chung, Yoonjae, Ranjit Shrestha, Seungju Lee, and Wontae Kim. 2022. "Binarization Mechanism Evaluation for Water Ingress Detectability in Honeycomb Sandwich Structure Using Lock-In Thermography" Materials 15, no. 6: 2333. https://doi.org/10.3390/ma15062333

APA Style

Chung, Y., Shrestha, R., Lee, S., & Kim, W. (2022). Binarization Mechanism Evaluation for Water Ingress Detectability in Honeycomb Sandwich Structure Using Lock-In Thermography. Materials, 15(6), 2333. https://doi.org/10.3390/ma15062333

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop