Preprint
Article

This version is not peer-reviewed.

Research on CFRP Defects Recognition and Localization Based on Metamaterial Sensors

A peer-reviewed article of this preprint also exists.

Submitted:

01 December 2024

Posted:

02 December 2024

You are already at the latest version

Abstract
The wide application of carbon-fiber-reinforced polymer(CFRP) has brought with it an urgent demand for their reliability and stability. However, considering factors such as the manufacturing process, material characteristics, and use environment, defects may exist in composite materials, posing a threat to engineering safety. Metamaterial sensors are widely used in non-destructive testing of microwave sensors due to their high resolution and high sensitivity. This paper proposes a defect identification and location method based on principal component analysis and support vector machine. The trained model is used to classify the dimensionally reduced data, and the reconstructed defect binary image is obtained. Simulation and physical experiment results show that the method used in this article can effectively identify and locate defects in carbon fiber composite materials.
Keywords: 
;  ;  ;  

I. Introduction

With the rapid progress of science and technology, CFRP has become an indispensable material in the industrial field, especially the aerospace field, due to its properties such as lightness, high stiffness, fatigue resistance and pressure resistance. The popularity of CFRP has brought about an urgent need for its safety. However, defects can occur in carbon fiber composites during the manufacturing process and use. Therefore, defect detection of composite materials is of great significance. In the industry, commonly used non-destructive testing technologies include ultrasonic testing, eddy current testing, thermal imaging testing and microwave testing. However, ultrasonic testing cannot penetrate highly porous materials and is limited by coupling agents, while also requiring extensive data interaction. Eddy current testing is widely used for surface crack and corrosion detection, but due to the penetration limitations of low-frequency electromagnetic waves, it is not effective when detecting low-loss dielectric materials.
A variety of microwave non-destructive testing technologies for composite material testing have been reported [1]. Although some promising results have been achieved, traditional microwave technologies still face some challenges, including high data complexity, space Poor image quality and blurry defect shapes, etc. Distance changes and optimal frequency selection can cause these problems, which in turn reduce the geometric measurement accuracy of defects [2]. Therefore, artificial intelligence technology is used in microwave non-destructive testing to solve the challenges based on microwave technology. Since machine learning technology has high capabilities in training and solving complex and nonlinear data, the ability to accurately classify defective and defect-free areas using machine learning classifiers will be reflected in measurement accuracy. Therefore, the use of intelligent classifiers can process and detect raw data while ensuring reliability and speed [3], and at the same time, it can improve the resolution of defect detection while automating the processing process [4]. In addition, artificial intelligence can solve the complexity problem when collecting data through non-linear classification [5].
An earlier attempt was made in [6] to introduce intelligent microwave non-destructive testing, using a combination of SVM and ANN classifiers to find small changes in reflection coefficients. A corrosion defect imaging system using microwave sensors and AI models was implemented in [7]. Feature extraction is performed using the supervised PCA algorithm. Implement artificial intelligence model through SVM algorithm. The trained model shows good generalization ability in defect representation. Frequency sweep-based features will lead to an increase in the amount of operations and calculations, and are very time-consuming in the initial stage of inspection, so machine learning feature selection is used in [8] to reduce the sweep frequency. Furthermore [9] developed and implemented an image gradient-based feature extraction method that includes a filter-based method and parameters from Weibull fitting. Liu et al. proposed a new loss function design that incorporates multiple scattering using near-field quantities to reconstruct the dielectric constant [10].
This paper is based on the non-destructive testing of CFRP using metamaterial sensors, studies defect identification and positioning technology, and proposes defect identification and positioning technology based on principal component analysis and support vector machines. Use a metamaterial sensor to conduct two-dimensional plane scanning of composite materials, and use principal component analysis to reduce the dimensionality of the original scan data obtained; then use a support vector machine to train the dimensionally reduced data with labels to obtain a classification model . Finally, the trained model is used to perform binary classification of defect-free and defective experimental data, and the reconstructed defect binary image is obtained, finally realizing the identification and positioning of defects.

II. Equivalent Model

A. Conductivity

Typically, research generally agrees that unidirectional carbon fiber composite laminates have anisotropic electrical properties. That is, the conductivity along the direction perpendicular to the fibers is relatively small, while the conductivity of the composite along the fiber direction can be calculated using the fiber volume fraction and the total single fiber conductivity. For an ideal unidirectional composite, the composite resistivity σ c p parallel to the fibers can be calculated as:
σ c p = σ f / v f
where v f is the volume fraction of the fiber; σ f is the resistivity of the carbon fiber. If the conductivity in all three directions is considered, the conductivity tensor of the carbon fiber composite plate can be expressed as:
σ = σ l cos 2 θ + σ t sin 2 θ σ l σ t / 2 sin 2 θ 0 σ l σ t / 2 sin 2 θ σ l cos 2 θ + σ t sin 2 θ 0 0 0 σ c p
where σ l represents the longitudinal conductivity of the layer; σ t represents the transverse conductivity of the layer; σ c p is the cross-layer conductivity; θ represents the direction of each ply.

B. Three-Layer Uniform Model

The characteristic of the three-layer uniform model is that the dielectric coefficient ε e f f of the central layer is constant throughout the layer. When microwaves propagate into fiber composite panels, they pass through three separate zones. The microwaves enter the panel first and see a uniform layer. This layer has constant material properties corresponding to the material properties of the matrix or resin surrounding the fibers. The thickness t of this layer corresponds to the distance from the edge of the panel to the edge of the fiber, and its dielectric coefficient is ε m . The central layer is the fiber layer. The wave finally enters the uniform layer again, which is the same as the first layer (i.e., constant material properties ε m and thickness t ). At any given cross-section in the uniform layer, the material properties approximate the effective material properties of the periodic laminate media shown in Figure 1.
This type of medium behaves like an anisotropic but homogeneous material with the following tensor permittivity and magnetic permeability:
ε = ε x 0 0 0 ε y 0 0 0 ε z
μ = μ x 0 0 0 μ y 0 0 0 μ z
where:
ε x 1 = 1 g ε a 1 + g ε b 1
μ x 1 = 1 g μ a 1 + g μ b 1
ε y = ε z = 1 g ε a + g ε b
μ y = μ z = 1 g μ a + g μ b
Fiber composite panels are structures composed of fibers aligned in only one direction, embedded in the center of a matrix or resin. Some of the composite structures analyzed here consist of two or more multilayer panels. Typically, the fibers in each panel of a composite structure can be oriented in different directions. Considering a multi-panel composite structure, when a three-layer model is used for this multi-panel structure, then each panel will be represented by three layers. Therefore, the multi-panel composite will have a periodic arrangement with carbon fiber and epoxy layers interlaced as shown in the model shown in Figure 1. At the same time, the fiber directions of different layers of carbon fiber are different, and the material parameters are also different, showing anisotropy. The light-colored part is the epoxy resin layer, and the dark-colored part is the carbon fiber layer. Its structure shows periodic characteristics.

C. Carbon Fiber Composite Structural Analysis

This paper established a carbon fiber composite structure simulation model in HFSS. The simulation model is shown in Figure 2. The composite material equivalent simulation model adopts a three-layer uniform model. Figure 2 shows the overall simulation model of carbon fiber composite material. Figure 1 shows the side view of carbon fiber composite material, that is, the layering situation of the composite material. It can be seen from the figure that each layer of material is mainly composed of carbon fiber layer and FR4 epoxy resin plate. Carbon fiber and FR4 epoxy resin plates are alternately arranged to form the overall simulation structure of carbon fiber composite materials. The materials of each layer of the model and their laying angles are shown in Table 1. Considering that in carbon fiber composite materials, different layers of carbon fibers are anisotropic, in the simulation model, different layers of carbon fibers are laid at different angles, namely along the x-axis and along the y-axis.

III. Methods of Defect Identification and Location

A. Analysis of Non-Destructive Testing of Resonant Sensors

Microwave resonance sensors utilize the propagation characteristics of microwaves to achieve detection by detecting the interaction between the target to be measured and the sensor. The sensor works in the microwave frequency range. When there is a defect in the target to be measured, it will affect the distribution of the field near the resonant sensing structure, causing changes in the relevant parameters of the sensing structure. Generally, it is reflected in the scattering parameter amplitude and resonance frequency. influence, thereby realizing the detection of the target to be tested. The detection schematic diagram is shown in Figure 3. Typically, the behavior of a resonant structure can be equivalent to a resonant circuit. In resonant sensor detection, if there are defects in the material being tested, the defects will affect the resonant sensing structure, causing the equivalent circuit of the original resonant sensing structure to change. After excitation of the resonant sensing structure using a two-port network, the changes caused by defects can be converted into detectable observations, such as the resonant frequency of the scattering parameters and their resonance amplitude.

B. Analysis of the Impact of Defects on Resonant Sensors

Assume that before defects are introduced in the resonant cavity, let ε 0 and μ 0 represent the dielectric constant and magnetic permeability of the resonant cavity in the initial state, E 0 and H 0 represent the electric field and magnetic field respectively before the defects are introduced, and ω 0 represents the initial state of the resonant cavity. Resonant frequency, V 0 represents the volume of the resonant cavity and S 0 represents the inner surface area of the resonant cavity; when introducing the target to be measured, let ε m and μ m represent the relative dielectric constant and relative permeability of the target to be measured respectively, and its volume is V m ; After introducing the perturbation of defects in the target to be measured, let ε m + ε and μ m + μ represent the relative permittivity and relative permeability of the defects in the target to be measured respectively, and the defect volume is V d . At this time, the electric field and magnetic field in the resonant cavity become E 1 and H 1 respectively, and the resonant frequency becomes ω 1 . Then the resonant frequencies ω 0 and ω 1 have the following relationship before and after introducing defect disturbance:
ω 1 ω 0 ω 1 = V d Δ ε E 1 E 0 * + Δ μ H 1 H 0 * d V V 0 V m ε 0 E 1 E 0 * + μ 0 H 1 H 0 * d V + V m ε m E 1 E 0 * + μ m H 1 H 0 * d V
Generally speaking, assuming that the changes in the dielectric constant and magnetic permeability of the resonant cavity after being affected by defects are relatively small, that is, ε and μ are relatively small, the electric and magnetic fields E 0 and H 0 of the original resonant cavity can be used to approximately replace E 1 and H 1 , so the above formula can be rewritten as:
ω 1 ω 0 ω 1 = V d Δ ε E 0 2 + Δ μ H 0 2 d V V 0 V m ε 0 E 0 2 + μ 0 H 0 2 d V + V m ε m E 0 2 + μ m H 0 2 d V

C. Principal Component Analysis

The algorithm flow of principal component analysis mainly consists of the following steps:
1)
The original data set is formed into a matrix X with n rows and m columns by column, indicating that the data set contains n samples, and each sample contains m variables;
2)
Demean each row (i.e., each sample) of the data set X, that is, subtract the mean of the current row, X ~ = X μ = X i = 1 m x i / n ;
3)
Solve the covariance matrix V of the demeaned matrix X ~ , V = X ~ X ~ T / n ;
4)
Solve the eigenvalue λ i of the covariance matrix V and its corresponding eigenvector w i ;
5)
Arrange the eigenvectors into a new matrix according to the size of the eigenvalues, and take the eigenvectors corresponding to the first k eigenvalues Λ = d i a g λ 1 , λ 2 , , λ k as the basis of the submatrix, that is, the matrix of the required k principal component is P = W d T X ~ ;
6)
The data set X = W P + μ is reconstructed based on the extracted k principal components, which achieves data dimensionality reduction.
Figure 4 shows the contribution rate of each principal component obtained by using PCA to process the original data obtained by scanning the measured material with the sensor, and the cumulative contribution rate of the principal components obtained by ranking according to the contribution rate. The original data is frequency sweep data with a frequency sweep range of 0.5GHz to 1GHz, that is, the S 21 amplitude of each frequency point. As can be seen from the figure, after PCA processing of the original data, the contribution rate of the first principal component is 74.9878%, the contribution rate of the second principal component is 16.096%, and the cumulative contribution rate of the first two principal components reached 91.0838%. The contribution rate up to the 23rd principal component can be basically ignored. Therefore, using PCA can realize feature extraction of the original data and characterize the original data without losing too much information.

D. Support Vector Machine Classifier

In order to obtain the model, the first step is to collect data. This article uses sensors to obtain frequency sweep data of the tested material, and then uses PCA to reduce the dimension. Then the data set is divided into a training set and a test set, and the data is classified. Unification; the next step is to choose an appropriate kernel function. The quadratic polynomial kernel function used in this article is. At the same time, cross-validation is used to find the best parameters →C and →ξ; then the best parameters →C and →ξ are used to train the entire training set to derive the model.
The SVM model, which uses a quadratic polynomial as the kernel function, classifies two types of problems, and the results are shown in Figure 5. Figure 5 is obtained from an SVM classifier applied to data reduction, providing a visual image of the learned two-dimensional model. Defective scans (red circles) and defective scans (blue triangles) are classified.

IV. Experiments and Results

A. Sensor Settings

The resonant sensing structure of the metamaterial sensor is shown in Figure 6a, which consists of nested CSRR etched on the ground plate of the microstrip transmission line and wave ports at both ends, combined to form a dual port sensor model. The sensing structure uses nested CSRR. Load CSRR on a grounding copper plane with a thickness of 35μm, using a 0.1mm thick FR4 epoxy resin board as the substrate. The top plane includes a 1.86mm wide microwave transmission line as the electromagnetic field source between ports 1 and 2. Figure 6b illustrates the setup of the resonant sensing structure model embedded with CSRR. In the CSRR unit structure, l1 = w1 = 12mm, l2 = w2 = 4.4mm, and g = 0.2mm. The length of the base lCSRR is 60mm, and the width lCSRRis 20mm. The simulation excitation method adopts lumped port excitation.

B. Internal Defects in the MUT

This paper constructed a carbon fiber composite material model of 80mm × 80mm with a thickness of 1.9mm. The thickness of the carbon fiber layer of the composite material is 0.5mm, and the thickness of the FR4 epoxy resin layer is 0.1 mm. The internal bubble defect is located in the FR4 epoxy resin layer. The bubble defect area is 5mm × 5mm × 0.05mm, and the depth is 0.625mm. The lifting distance from the resonant sensor to the composite material is 0.7mm.
When the sensor scans the tested composite material, it covers an area of 20 mm × 20 mm, with a step size of 1 mm, and constructs a 441-pixel image in a two-dimensional scanning mode. Figure 7a,b show the top and side views of the tested composite materials. For each scan position, the frequency is scanned in the frequency range 0.5GHz to 1.0GHz in 2MHz increments. The scanned results are shown in Figure 8. Figure 8 shows an example of the amplitude-frequency plot of the S 21 parameter recorded at each scanning position, indicating that when the CSRR is close to the defect, an obvious frequency shift occurs.
The results of the scattering matrix are imported into the MATLAB environment for post-processing, and the simulation results are shown in Figure 9. The yellow rectangular area in Figure 9a shows the defective area, and the blue and green backgrounds correspond to the rest of the scanned area. The S 21 resonant frequency shift at the center of the internal defect area is increased to 10MHz. However, when the sensor moves away from the defect area, the S 21 resonant frequency offset starts to decrease.
The defect reconstruction image of the tested composite material is shown in Figure 10. Black represents the defective area and white represents the defect-free area. The reconstructed image is based on the data obtained by the metamaterial sensor. PCA is used for data dimensionality reduction. Considering the binary classification situation, SVM is used to implement the model. The size of the reconstructed defect is basically consistent with the size of the simulation setting. At the same time, due to the use of binary images, the imaging of the defect can be more prominent. The trained model can better characterize defects in carbon fiber composite materials.

C. Experimental Results and Analysis of Aircraft Skin

This paper uses aircraft skin composite material as the MUT, which is also a kind of carbon fiber composite material, and its structure is shown in Figure 11. Figure 11 shows the top view and cross-sectional view of the material, with white dots indicating the location of defects. The length of this material is 6.7 cm, and the width of the top is 3cm. The defect size is 1cm.
The experimental platform includes a network analyzer, a resonance sensor, and the tested material. The experimental platform is shown in Figure 12. In the experiment, the metamaterial sensor was placed on the top of the material being tested, and a lifting distance was maintained between the sensor and the top plane of the material being tested. A sensor is used to perform a two-dimensional planar scan of the area near the defect on top of the material being tested.
For each scanning position, the frequency is scanned in the frequency range from 0.364GHz to 1.220GHz, and the scanning results are shown in Figure 13. The experiment uses a vector network analyzer to obtain frequency sweep data, and also imports the original data obtained from the VNA into the Matlab environment for processing. The simulation results are shown in Figure 14. The S 21 resonant frequency shift in the defective area is increased to 6MHz. The three-dimensional stereoscopic image constructed in Figure 14b can more intuitively show that the S 21 resonant frequency shift at the defective position increases, while the S 21 resonant frequency shift at the defect-free position is small or almost no shift.
parameter for all scanning positions.
The reconstruction result of the defect was measured as shown in Figure 14c. The size of the reconstructed defect is basically consistent with the size of the actual object, and a binary image is used, so the imaging of the defect can be more prominent. There are also some redundant points on the black image boundary. According to analysis, it is considered that the model is not accurate enough or the sensor resolution is insufficient due to insufficient data sets for training the model. There are also problems trained based on the carbon fiber composite simulation model. The problem of mismatch between classification model and real objects.

V. Conclusion

This paper constructs simulation models of carbon fiber composite materials and metamaterial sensors, uses metamaterial sensors to scan defects near the tested material defects, and uses defect recognition and localization technology based on PCA and SVM to process the original scanning data to obtain reconstructed binary images of defects, achieving defect recognition and localization. The experimental results indicate that the defect recognition and localization technology based on principal component analysis and support vector machine used in this article can effectively characterize defects.

References

  1. T. Wang, D. Wu, W. Chen, and J. Yang. Detection of delamination defects inside carbon fiber reinforced plastic laminates by measuring eddy-current loss. Composite Structures 2021, 268, 114012. [Google Scholar] [CrossRef]
  2. S. Kharkovsky and R. Zoughi, “Microwave and millimeter wave nondestructive testing and evaluation - Overview and recent advances. IEEE Instrumentation & Measurement Magazine 2017, 10, 26–38. [Google Scholar]
  3. M. S. M. Naqiuddin, M. S. Leong, and L. M. Hee, “Ultrasonic signal processing techniques for Pipeline: A review. MATEC Web of Conferences 2019, 255. [Google Scholar]
  4. S. Sambath, P. Nagaraj, and N. Selvakumar, “Automatic defect classification in ultrasonic NDT using artificial intelligence. Journal of Nondestructive Evaluation 2011, 30, 20–28. [Google Scholar] [CrossRef]
  5. R. D. Tipones and J. C. D. Cruz, “Design and development of a material impact tester using neural network for concrete ratio classification,” in Proc. IEEE 13th International Colloquium on Signal Processing & its Applications (CSPA), 2017, pp. 106–111.
  6. Ali, B. Hu, and O. Ramahi, “Intelligent Detection of Cracks in Metallic Surfaces Using a Waveguide Sensor Loaded with Metamaterial Elements. Sensors 2015, 15, 11402–11416. [Google Scholar] [CrossRef]
  7. A. Ali, A. A. Ali, A. Albasir, and O. M. Ramahi, “Microwave sensor for imaging corrosion under coatings utilizing pattern recognition. Proc. IEEE International Symposium on Antennas & Propagation 2016, 951–952. [Google Scholar]
  8. Moomen, A. Ali, and O. Ramahi, “Reducing Sweeping Frequencies in Microwave NDT Employing Machine Learning Feature Selection. Sensors 2016, 16, 559. [Google Scholar] [CrossRef]
  9. Y. D. Kotriwar, O. Y. D. Kotriwar, O. Elshafiey, L. Peng, et al. Gradient feature-based method for defect detection of carbon fiber reinforced polymer materials. Proc. IEEE International Conference on Prognostics and Health Management (ICPHM) 2023, 246–252. [Google Scholar]
  10. Z. Liu, M. Roy, D. K. Prasad, et al. Physics-guided loss functions improve deep learning performance in inverse scattering. IEEE Transactions on Computational Imaging 2022, 8, 236–245. [Google Scholar] [CrossRef]
  11. Zhaoxuan zhu,. Ph.D. degree,His main research interests include aviation energy testing and testing technology, aviation energy automation, non-destructive testing, and high-speed data acquisition and processing technology.
  12. Rui Han. Graduate. Main research directions: aviation energy testing and detection technology, non-destructive testing.
Figure 1. Simulation model of carbon fiber composite materials.
Figure 1. Simulation model of carbon fiber composite materials.
Preprints 141430 g001
Figure 2. Carbon fiber composite overall simulation model.
Figure 2. Carbon fiber composite overall simulation model.
Preprints 141430 g002
Figure 3. Schematic diagram of resonance sensor detection.
Figure 3. Schematic diagram of resonance sensor detection.
Preprints 141430 g003
Figure 4. (a) The contribution rate of each principal component; (b) Principal component cumulative contribution rate based on contribution rate arrangement.
Figure 4. (a) The contribution rate of each principal component; (b) Principal component cumulative contribution rate based on contribution rate arrangement.
Preprints 141430 g004
Figure 5. Decision result.
Figure 5. Decision result.
Preprints 141430 g005
Figure 6. (a) Resonant sensing structure; (b) Structural parameter settings.
Figure 6. (a) Resonant sensing structure; (b) Structural parameter settings.
Preprints 141430 g006
Figure 7. (a) Defect top view; (b) Defect side view.
Figure 7. (a) Defect top view; (b) Defect side view.
Preprints 141430 g007
Figure 8. S 21 parameter frequency-amplitude plot of 441 scan points near internal defects.
Figure 8. S 21 parameter frequency-amplitude plot of 441 scan points near internal defects.
Preprints 141430 g008
Figure 9. (a) S 21 resonant frequency offset size at each scan position; (b) 3D plot of S 21 resonant frequency shift at each scan position.
Figure 9. (a) S 21 resonant frequency offset size at each scan position; (b) 3D plot of S 21 resonant frequency shift at each scan position.
Preprints 141430 g009
Figure 10. Reconstructed image of internal defects of the MUT.
Figure 10. Reconstructed image of internal defects of the MUT.
Preprints 141430 g010
Figure 11. Aircraft skin. (a) Top view; (b) Cross section diagram.
Figure 11. Aircraft skin. (a) Top view; (b) Cross section diagram.
Preprints 141430 g011
Figure 12. Testing experimental platform.
Figure 12. Testing experimental platform.
Preprints 141430 g012
Figure 13. Frequency amplitude plot of S 21
Figure 13. Frequency amplitude plot of S 21
Preprints 141430 g013
Figure 14. (a) S 21 resonant frequency offset size at each scan position; (b) 3D plot of S 21 resonant frequency shift at each scan position; (c) Reconstructed images of aircraft skin composite defects.
Figure 14. (a) S 21 resonant frequency offset size at each scan position; (b) 3D plot of S 21 resonant frequency shift at each scan position; (c) Reconstructed images of aircraft skin composite defects.
Preprints 141430 g014
Table 1. Model Ply Material and Laying Angle.
Table 1. Model Ply Material and Laying Angle.
Number of layers Material name Laying angle
1 FR4-epoxy non
2 carbon fiber 0
3 FR4-epoxy non
4 carbon fiber 90
5 FR4-epoxy non
6 carbon fiber 0
7 FR4-epoxy non
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

Subscribe

Disclaimer

Terms of Use

Privacy Policy

Privacy Settings

© 2025 MDPI (Basel, Switzerland) unless otherwise stated