Next Article in Journal
Wearable Noninvasive Glucose Sensor Based on CuxO NFs/Cu NPs Nanocomposites
Next Article in Special Issue
Optimal Transducer Placement for Deep Learning-Based Non-Destructive Evaluation
Previous Article in Journal
Acoustic Emission and Artificial Intelligence Procedure for Crack Source Localization
Previous Article in Special Issue
Crack Growth Monitoring with Structure-Bonded Thin and Flexible Coils
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Simulation of Layer Thickness Measurement in Thin Multi-Layered Material by Variable-Focus Laser Ultrasonic Testing

1
School of Information Engineering, Wuhan University of Technology, Wuhan 430070, China
2
Shaanxi ERC of NDT and Structural Integrity Evaluation, State Key Laboratory for Strength and Vibration of Mechanical Structures, Xi’an Jiaotong University, Xi’an 710049, China
3
State Key Laboratory of Advanced Technology for Materials Synthesis and Processing, Wuhan University of Technology, Wuhan 430070, China
*
Author to whom correspondence should be addressed.
Sensors 2023, 23(2), 694; https://doi.org/10.3390/s23020694
Submission received: 8 December 2022 / Revised: 30 December 2022 / Accepted: 5 January 2023 / Published: 7 January 2023

Abstract

:
Thin multi-layered materials are widely used in key structures of many high technology industries. To ensure the quality and safety of structures, layer thickness measurement by non-destructive testing (NDT) techniques is essential. In this paper, a novel approach for the measurement of each layer’s thickness in thin multi-layered material is proposed by using ring-shaped laser generated focused ultrasonic bulk waves. The proposed method uses a ring-shaped laser with a variable radius to generate shear waves with variable focus inside the structure. By analyzing the signal characteristics at the ring center when the laser radius varies from zero to maximum, the direct measurement of layer thickness can be realized, considering that only when the focal depth and the layer thickness satisfy the specific relationship, the reflected shear waves converge and form a peak at the ring center. This straightforward approach can increase the pulse-echo SNR and prevent the processing of aliasing signals, and therefore provides higher efficiency and accuracy for the layer thickness measurement. In order to investigate the feasibility of this method, finite element simulations were conducted to simulate the ring-shaped laser generated ultrasonic waves in multi-layered structure in detail. Following the principle of the proposed method, the layer thickness of a bi-layer and 3-layer structure were respectively measured using simulation data. The results confirm that the proposed method can accurately and efficiently measure the layer thickness of thin multi-layered material.

1. Introduction

Thin multi-layered materials such as functional gradient materials, polymer laminates, and thermal barrier coatings are frequently employed in major aerospace engineering, energy, and other high-technology industrial structures [1,2,3]. The thickness of each layer must be measured using non-destructive testing (NDT) methods to assure the quality and safety of such structures because variations in layer thicknesses might reduce the mechanical qualities of the structure, or possibly result in its failure.
For measuring material thickness, several NDT techniques have been reported until now [4,5]. One of the most widely used NDT techniques among them is ultrasonic testing (UT), which has been used extensively for thickness measuring, material performance testing, and defect inspection [6]. However, it is challenging to use standard UT, such as the pulse-echo based approach, to measure the layer thickness of thin multi-layered material. The reasons can be listed as follows: (1) The overlap of ultrasonic echoes from interlayers is severe because layer thickness is narrow (often in the micrometer range). (2) Echoes from deep levels are greatly diminished as a result of the dispersion and reflection of numerous contacts. The severe signal aliasing and poor SNR of interface pulse echoes make it difficult to extract and analyze the time of flight information of ultrasonic signals. Current research focuses on the separation of overlapped pulse echoes using signal processing [7,8,9] and parameter optimization using the inversion algorithm [10,11,12,13] to address the aforementioned issues. However, the method of indirect thickness evaluation is vulnerable to a difficult inversion process or subpar evaluation accuracy.
Laser ultrasonic testing (LUT), a recent addition to the UT family, has demonstrated its capacity to produce ultra-high frequency and broadband ultrasonic waves with a non-contact characteristic. In several earlier studies, LUT has been researched and used for inspecting surface or near-surface defects. In addition, several reports [14,15] have appeared where thickness measurements were made utilizing longitudinal waves produced by lasers. Related investigations that use the pulse-echo approach, nevertheless, have found that signal aliasing issues are still present. Additionally, the thermoelastic regime generates bulk waves with weak signals. This effect significantly reduces the pulse echo’s SNR in deep layers and restricts its use in the assessment of thickness. Research has been conducted on the modulation of pulse laser patterns, such as the ring-shaped laser ultrasonic testing (RLUT) approach, to increase the signal amplitude without damaging the material surface by ablation. In contrast to surface waves produced by a line source, P. Cielo et al. [16] have reported a 20-fold increase in signal level using a ring-shaped pulsed laser to generate convergent surface waves. To find cracks on the bottom surface, X. Wang et al. originally employed ring-shaped laser-generated ultrasound [17]. The theoretical model and finite element simulation of an acoustic wave produced by a ring laser were developed by J. Guan et al. [18]. The benefits of ring focused laser ultrasound in fracture identification over non-focused laser ultrasound were further enhanced by S. Dixon [19]. The superposition of generated bulk waves is made possible at a particular focus point by the use of a ring-shaped laser. SNR and spatial resolution can be significantly increased as a result. The research reports mentioned above, however, only focused on the inspection of surface or interior flaws and cannot be used to assess the thickness of thin material, particularly in multi-layered structures.
In this study, a novel variable-focus ring-shaped laser ultrasonic testing (VRLUT) technique is proposed for measuring the layer thickness of thin layered material. The ultrasonic waves with varied focus are generated inside the structure by repeatedly exciting ring-shaped lasers of various diameters on the material surface. Our research found that only when the focal depth and the layer thickness satisfy the specific relationship do the reflection waves converge at the ring center. This phenomenon can be witnessed from A-scan signals, and subsequently allows the analysis of the quantitative relationship between focal depth and layer thickness. By matching the focal depth and layer interface locations, direct thickness measurement is then possible, avoiding issues with traditional approaches, such as complex parameter inversion and the separation of echo signals with severe aliasing. Additionally, the employment of deep-focused waves can successfully raise the signal-to-noise ratio (SNR) of echoes from deep layers, which have a high degree of attenuation and dispersion.
This paper is organized as follows. The principle of the proposed method is introduced in Section 2. Numerical simulations using the finite element method (FEM) to examine the viability of this approach is introduced in Section 3. Measurement results for various thin stacked metal plates by using simulation data are given in Section 4. Finally, the conclusion is presented in Section 5.

2. Principle of VRLUT-Based Layer Thickness Measurement Method

2.1. Focus Pattern of Laser-Generated Bulk Waves

There is a directivity pattern in the ultrasonic bulk waves produced by lasers. The amplitude of longitudinal and shear waves generated by thermoelastic expansion is a function of angle, and is defined in Equations (1) and (2) [20].
u L sin θ sin 2 θ κ 2 sin 2 θ ( κ 2 2 sin 2 θ ) 2 + 4 sin 2 θ cos θ κ 2 sin 2 θ
u s sin 2 θ ( 1 2 sin 2 θ ) ( 1 2 sin 2 θ ) 2 + 4 sin 2 θ cos θ κ 2 sin 2 θ
where
u L : Amplitude of longitudinal wave,
u L : Amplitude of shear wave,
θ : Angle to the surface normal,
κ = C L / C S : Ratio of the longitudinal wave velocity to the shear wave velocity.
Figure 1 shows the directivity patterns of the laser-generated ultrasonic waves in an iron plate with CL/CS = 1.86. While the amplitude of the shear wave rapidly decreases in other directions, it exhibits a dramatic peak at around 34° normal to the surface. In contrast, the longitudinal wave’s amplitude exhibits a peak at around 69° normal to the surface, but does not rapidly decelerate in other directions.
If the ring laser source is used as a point source assessment, each point source will produce a shear wave and a directive longitudinal wave that overlap at the angle of maximum amplitude. Then, at a particular focal depth below the center of the laser ring, a focal area is generated. The focal area, particularly for shear waves, is so narrow that it can be regarded as a point. As a result, this study focused on the application of shear waves and their primary characteristics.

2.2. Layer Thickness Measurement Method

The designed system and principle for thickness measurement using RLUT are shown in Figure 2. The initial point-like laser beam is generated by a Q-switch solid-state pulsed laser with duration ≤ 8 ns, wavelength of 1064 ns, and maximum pulse energy of 80 mJ. A beam expander initially enlarges the initial point while measuring. After that, the enlarged beam is sequentially passed via a diffractive axicon (DA) and a convex lens (CL). The focused laser source is then created on the specimen surface in the shape of a ring. To change the distance to the specimen, the CL and DA are mounted on two micro-displacement scanning stages with stepping accuracy of 6.5 μm. The ring radius can be modified by changing the DA while keeping the ring width constant because the ring width only depends on the distance between the convex lens and the specimen. Therefore, by gradually increasing the distance between DA and the specimen, it is possible to irradiate the specimen surface with a ring-shaped laser whose radius continually varies from zero to maximum. Through the use of a laser interferometer and a dichroic mirror (DM), the surface displacement brought about by the generated ultrasonic waves is detected. The laser interferometer can be a single-point heterodyne type with bandwidth of at least 500 MHz. A high-speed oscilloscope with bandwidth of 30 GHz and sampling rate of 80 GS/s is used to record and later process the ultrasonic signals that the interferometer observed. The testing should be conducted at room temperature and indoor environment to ensure the accuracy of optical measurement system.
As depicted in Figure 3a, while performing the measurement, the ring-shaped laser, with radius varying from zero to maximum, is generated on the specimen surface to conduct the R-scan procedure. Supposing the reflection angle of reflected bulk waves is absolutely the same as incident angle of initial waves, the superposition of reflected shear waves will be observed at the ring center if the focal depth d, angle of maximum shear amplitude θ, and the ring radius R follow:
2d = Rtanθ,
As shown in Figure 3b, the reflection waves will converge in the center of the laser ring as illustrated by the blue lines, and a peak will be visible in the time-domain waveform of the center while the ring radius R fits the layer thickness h, which means h equaling d. In other cases, peaks will occur at the ring around the center, but not exactly at the center, as illustrated by the red lines.
Figure 4 displays the typical RLUT signal at the ring center when R matching h and R unmatching h. Due to its higher velocity compare to other waves, the surface longitudinal wave (sL) in both situations arrives at the ring center first. Due to the very modest shear wave component in this direction (as shown in Figure 1), the surface shear wave (sS) with low amplitude then emerges. The Rayleigh wave has a substantially bigger amplitude than other waves because of its comparatively high initial amplitude and center overlap. Different phenomena may then be seen: when R matches h, the superposition of reflection waves (2S) occurs at the center, and when R unmatches h, the amplitude is significantly lower than in the former situation. This allows one to determine which R best matches the h after evaluating the signals produced by various radius ring-shaped lasers and contrasting the signals at the center. Then the thickness of the first layer can be calculated by Equation (3). It should be noted that only shear waves are employed in this method since they are straightforward to identify, have a high signal-to-noise ratio, have a smaller maximum amplitude angle, and have a more evident superposition effect than longitudinal waves. However, for a particular ring radius, not only the shear echo but also the longitudinal echo and reflected mode-converted waves may overlap at the center and give a deceptive indication of the matching R. Therefore, the relationship between the flight time and propagation distance is utilized to further verify the calculated h using:
c s 1 t s 1 = 2 ( R / 2 ) 2 + h 1 2 ,
where Csi and tsi are the shear wave velocity and flight time in ith layer, respectively.
Only when the calculated h simultaneously satisfies Equations (3) and (4), the h is regarded as a true value. For the second layer, as shown in Figure 3c, Equation (3) can be written as:
R = 2 ( h 1 tan θ 1 + h 2 tan θ 2 ) ,
According to Snell’s law:
sin θ 1 sin θ 2 = c s 1 c s 2 ,
The thickness of the second layer can be estimated once the first layer’s thickness and the shear wave velocities in all layers are known. Following the flight time of superposition shear waves at the center point, the predicted thickness should also be:
2 h 1 cos θ 1 + c s 2 ( t s 2 2 h 1 c s 1 cos θ 1 ) = 2 ( h 1 cos θ 1 + ( R 2 h 1 tan θ 1 ) 2 + h 2 2 ) ,
The thickness measurement of a bi-layer material is then achieved. In the same way, it is possible to determine each layer’s thickness layer by layer in a multi-layered structure with n layers by:
R = 2 i = 1 n h i tan θ i ,
where hi and θi are respectively the thickness and angle of maximum shear amplitude in ith layer.

3. Simulation of Thickness Measurement by VRLUT

3.1. Simulation Model and Main Parameters

The proposed thickness measurement technique was validated by the analysis of ultrasonic waves using the finite element method (FEM). A commercial software COMSOL was used to conduct the FEM simulation. Two numerical models were developed, as shown in Figure 5, to depict the thickness measurement process and examine the precision of the proposed method. A cylindrical coordinate system is also introduced in our model to reduce computation costs. As a result, the issue can be thought of in the axial symmetric plane of two dimensions. The simulation models contain a W-Cu bi-layer plate and an Al-Fe-Cu three-layer plate. On the top surface of the first layer in both models, a ring-shaped laser is irradiated with energy below the sample damage threshold. The ring radius and half of the irradiation width are respectively denoted by R and w. The thickness of each layer of the n-layer material is denoted by hi (i = 1, 2, …… n) from top to bottom. The spatial distribution f(r) and temporal distribution g(t) of the laser source in the simulation model can be expressed as:
f ( r ) = exp ( ( r R ) 2 / w 2 ) ,
g ( t ) = 8 t 3 / t 0 3 exp ( 2 t 2 / t 0 2 ) ,
where t0 is the rise time of the laser pulse.
According to some reports, the laser energy, ring radius, and width can all have an impact on the laser-generated ultrasonic waves [21,22]. To obtain a significant signal amplitude, the laser energy should be set as high as feasible while yet staying below the thermoelastic mode threshold. The accuracy of the measurement can be increased by employing a thin width ring laser since the frequency of ultrasound has an inverse relationship with ring width. However, the diffraction limit of the optical components places a limit on the ring width, and a smaller laser width would require finer meshing, increasing the cost of the calculation. To balance the calculation cost and the accuracy of the numerical results, the laser pulse width and laser ring width are assumed to be 1 ns and 2 μm respectively. The laser’s peak power density was set at 0.012 W/m−2. During the R-scan procedure, the radius of the laser ranges from 0 to 500 μm with a step of 2 μm. The element size of the grid was set as 0.1 μm in the vicinity of the laser-affected area, while the element size was set at 0.5 μm outside the heat-affected zone. Details of material parameters used in the simulation are listed in Table 1. All parameters are from COMSOL’s built-in materials.

3.2. Calibration of Wave Velocity and Focal Angle

Before the application for thickness assessment, simulation is used to confirm the focal properties of the shear wave produced by the ring-shaped laser. Moreover, the maximum shear amplitude angle might not match the theoretical value because of the attenuation of ultrasonic waves and the fact that the focus region is not a precise point. The velocity of the shear wave is also unknown, according to Table 1. Therefore, it is crucial to calibrate the maximum shear amplitude angle and wave velocity before the measurement procedure. In order to calibrate and determine the wave velocity and focal angle, simulations were firstly conducted by using a one-layer model. In this model, a ring-shaped laser with radius of 50 μm and width of 2 μm is used to irradiate the top surface of a copper plate with thickness of 100 μm. The 3D view of the simulated ultrasonic field inside the plate at different times is shown in Figure 6. The bulk waves have not yet overlapped at time t = 10 ns. Then at t = 15 ns, the longitudinal waves began to overlap; however, because the focus effect was not obvious, the amplitude of the waves in the overlapped area was relatively low. At t = 30 ns, the shear waves overlap, and an obvious focus point is observed. Finally, at t = 50 ns, the focus point vanished after the shear waves split. The above results were helpful in understanding the focus pattern of bulk waves, and it could be clearly observed that shear waves perform better in focus.
The calibration of shear wave velocity and focal angle can be then accomplished by using above simulation results. From the geometrical relationship as shown in Figure 6c, the focal angle θ can be obtained by:
tan θ = f d / R ,
where R is the radius of laser ring and fd is the distance between laser ring center and focal point. The shear wave velocity vs is obtained according to the flight time t:
v s = f d 2 + R 2 / t ,
As the focal point is thought to be located at the coordinate of the maximum peak value, the wave amplitudes along the Z-axis were then extracted to find out the coordinate of focal point. The extracted amplitudes of the shear wave along the Z-axis from 38 ns to 42 ns are shown in Figure 7. The shear wave’s amplitudes exhibit a unimodal distribution, as expected, and the maximum peak-to-peak amplitude is attained at a position close to Z = 36.8 μm with flight time t = 40 ns. The velocity of shear wave vs = 2040 m/s and focal angle θ = 38.0° can then be calculated according to Equations (11) and (12), which is larger than the angle theoretically analyzed by Equation (2). The reason for the error is that the propagation distance of ultrasonic wave in this direction is smaller than that in the theoretical one, which means the attenuation is simultaneously smaller. These calibrated values served as the ideal values in this study. Similarly, the velocity of the shear wave and the focal angle in the aluminum plate is calculated as vs = 3014 m/s and the focal angle θ = 36.4°.

4. Measurement Results

By measuring the thickness of multilayered materials, the performance of the suggested method is examined to validate it. In Section 3.1, the simulation model and computation setup have been shown. The bi-layer model’s typical simulation results for the volumetric ultrasonic displacement field at various times when exposed to a 400 μm radius ring laser are shown in Figure 8. While there is no focus to be seen in Figure 8a,c, Figure 8b shows a large amplitude of displacement in the center of the ring that indicates a good convergence of shear echoes. The time domain signals of the ring center are extracted as illustrated in Figure 9a,b. Only when R varies near 190 μm and 410 μm, which correspond to the thickness of the first layer and second layer, respectively, can the peak value behind the Rayleigh wave, indicating the focal shear echoes, be seen. Then the amplitude of peak value Δd is furtherly compared, as shown in Figure 10. Finally, the optimum matching radius for the two-layer is selected as 190 μm and 412 μm. The thickness is determined to be 121.59 μm and 80.54 μm, respectively, via Equations (3) and (5). Additionally, the computed thickness values are verified by flight duration. The predicted flight times of overlapped shear echoes are estimated by Equations (4) and (6) to be around 151 and 255 ns, which is quite similar to the values seen in Figure 10. The simulation results for the 3-layer model presented in Figure 11 demonstrate a similar tendency. The optimum matching radii for the three layers were selected as 74 μm, 152 μm, and 178 μm, and the thicknesses were calculated as 50.19 μm, 50.30 μm, and 30.36 μm, respectively. The thickness measurement results of each layer are listed in Table 2. The aforementioned results demonstrate that the suggested method can accurately determine the thickness of every layer with an error of less than 1.5%. Additionally, each measurement step takes 200 ns for signal acquisition and 0.2 s for automatically adjusting R, resulting in a total measurement time of less than 1 min for the 250 steps when R ranges from 0 to 500 μm.
It should be noted that the proposed method and all results demonstrated above are valid only when following conditions are met: (1) The material in each layer is isotropic and homogeneous. In this case, the propagation of shear wave in a layer can be regarded as along one direction. (2) The total thickness of tested specimen is much smaller than other dimensions, so that the wave reflection of the specimen edge does not produce a large influence on the signal analysis. Due care should also be taken to ensure the reproducibility of the obtained values: (1) The distribution of laser power density along the ring direction should be uniform. (2) The surface of the specimen should be smooth and clean. In this way, we can ensure that, in each measurement, the amplitude of focal shear echoes at the ring center is able to keep constant and reach maximum when R matches the thickness. In practical applications, an optical homogenizer may be introduced to improve the uniformity of the laser beam when the laser quality is poor. For the tested specimen of which the surface condition is undesirable, proper surface treatment such as polish operation or painting may be desired.

5. Conclusions

In this study, a variable-focus ring-shaped laser ultrasonic testing (VRLUT) technique is proposed for the measurement of each layer’s thickness in thin multi-layered material. The matching feature between layer thickness and focus position is retrieved by examining the signal characteristics at various relative positions between the focus point and the multilayer interface. As a result, it is possible to directly quantify layer thickness using the matching relationship between the radius of laser ring and the thickness of the layer combining with the flight time of focused echoes. The proposed method may deliver a high SNR of measuring echoes while also avoiding the complicated data processing process, which is able to effectively measure the layer thickness of thin multi-layered structures in practice with higher accuracy and speed compared to current approaches. Furthermore, as a fully non-contact testing method, the proposed technique has the potential to be applied for the on-line measurement of material thickness during manufacturing.
The proposed technique is validated with use of simulated signals to assess the layer thickness of a bi-layer and 3-layer structure. The results show that the developed method is well adapted to measure multiple layer thicknesses that are even less than 50 μm. For both the bi-layer and 3-layer structure, their layer thickness was obtained with relative error of less than 1.5%, which proved the validity of proposed method and corresponding algorithm. The above findings suggest that the VRLUT-based approach can offer a very efficient and precise NDT technique for assessing multiple layer thicknesses in thin multi-layered material.
Furthermore, the measurement system will be established, and experimental research needs to be conducted to examine the realism and applicability of these novel methodologies for measuring layer thickness.

Author Contributions

Conceptualization, J.Q. and Z.L.; methodology, J.Q. and C.P.; software, J.Q.; validation, Z.L. and J.Q.; formal analysis, J.Q.; investigation, J.Q., G.L. and Z.L.; resources, G.L.; data curation, J.Q.; writing—original draft preparation, J.Q.; writing—review and editing, Z.L. and C.P.; visualization, J.Q.; supervision, Z.L.; project administration, Z.L.; funding acquisition, J.Q. and Z.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (grant number: No. 52205574, No. 62275205) and the Fundamental Research Funds for the Central Universities (grant number: WUT:2021IVA007B).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhang, L.; Liu, Y.; Zhang, C.; Luo, G.; Liu, H.; Shen, Q. Fabrication and Properties of W-SiCP/Cu Composites by Hot Pressing Sintering. Key Eng. Mater. 2016, 697, 207–210. [Google Scholar] [CrossRef]
  2. Sarkeeva, A.; Kruglov, A.; Lutfullin, R.; Gladkovskiy, S.V.; Zhilyaev, A.P.; Mulyukov, R.R. Characteristics of the mechanical behavior of Ti-6Al-4V multilayer laminate under impact loading. Compos. Part B Eng. 2020, 187, 107838. [Google Scholar] [CrossRef]
  3. Padture, N.; Gell, M.; Jordan, E. Thermal Barrier Coatings for Gas-Turbine Engine Applications. Science 2002, 296, 280–284. [Google Scholar] [CrossRef] [PubMed]
  4. Yang, X.; Verboven, E.; Ju, B.; Kersemans, M. Comparative study of ultrasonic techniques for reconstructing the multilayer structure of composites. NDT E Int. 2021, 121, 102460. [Google Scholar] [CrossRef]
  5. Kolpatzeck, K.; Liu, X.; Haring, L.; Balzer, J.C.; Czylwik, A. Ultra-High Repetition Rate Terahertz Time-Domain Spectroscopy for Micrometer Layer Thickness Measurement. Sensors 2021, 21, 5389. [Google Scholar] [CrossRef]
  6. Islam, M.T.; Rahman, M.N.; Samsuzzaman, M.; Mansor, M.F.; Misran, N. Resonator-Inspired Metamaterial Sensor: Design and Experimental Validation for Measuring Thickness of Multi-Layered Structures. Sensors 2018, 18, 4213. [Google Scholar] [CrossRef] [Green Version]
  7. Chen, J.; Bai, X.; Yang, K.; Ju, B.F. An ultrasonic methodology for determining the mechanical and geometrical properties of a thin layer using a deconvolution technique. Ultrasonics 2013, 53, 1377–1383. [Google Scholar] [CrossRef]
  8. Lu, Z.; Yang, C.; Qin, D.; Luo, Y.; Momayez, M. Estimating ultrasonic time-of-flight through echo signal envelope and modified Gauss Newton method. Measurement 2016, 94, 355–363. [Google Scholar] [CrossRef]
  9. Chapon, A.; Pereira, D.; Toews, M.; Belanger, P. Deconvolution of ultrasonic signals using a convolutional neural network. Ultrasonics 2021, 111, 106312. [Google Scholar] [CrossRef]
  10. Yang, X.; Zhang, C.; Sun, A.; Bai, X.; Ju, B.; Shen, Q. Numerical and experimental analysis of a focused reflected wave in a multi-layered material based on a ray model. Ultrasonics 2018, 86, 41–48. [Google Scholar] [CrossRef]
  11. Ma, Z.; Luo, Z.; Lin, L.; Krishnaswamy, S.; Lei, M. Quantitative characterization of the interfacial roughness and thickness of inhomogeneous coatings based on ultrasonic reflection coefficient phase spectrum. NDT E Int. 2019, 102, 16–25. [Google Scholar] [CrossRef]
  12. Zhang, J.; Xu, D.; Dong, F.; Wang, R.; Yang, Q.; Wang, X.; Xue, R. Solid–liquid interface reconstruction for sandwich structure metal plate via laser-ultrasonic techniques. Rev. Sci. Instrum. 2021, 92, 123003. [Google Scholar] [CrossRef] [PubMed]
  13. Geng, X.; Zhang, C.; Zhang, J.; Luo, G.; Shen, Q. Simultaneous measurement of layer thicknesses in thin layered materials using the phase of ultrasonic reflection coefficient spectrum. Appl. Acoust. 2022, 195, 108835. [Google Scholar] [CrossRef]
  14. Zhang, K.; Li, S.; Zhou, Z. Detection of disbonds in multi-layer bonded structures using the laser ultrasonic pulse-echo mode. Ultrasonics 2019, 94, 411–418. [Google Scholar] [CrossRef] [PubMed]
  15. Yi, K.; Liu, P.; Park, S.; Sohn, H. Femtosecond laser ultrasonic inspection of a moving object and its application to estimation of silicon wafer coating thickness. Opt. Lasers Eng. 2022, 148, 106778. [Google Scholar] [CrossRef]
  16. Cielo, P.; Nadeau, F.; Lamontagne, M. Laser generation of convergent acoustic waves for materials inspection. J. Ultrason. 1985, 23, 55–62. [Google Scholar] [CrossRef]
  17. Wang, X.; Littman, M.G.; McManus, J.B.; Tadi, M.; Kim, Y.S.; Askar, A.; Rabitz, H. Focused bulk ultrasonic waves generated by ring-shaped laser illumination and application to flaw detection. J. Appl. Phys. 1996, 80, 4274–4281. [Google Scholar] [CrossRef] [Green Version]
  18. Guan, J.; Shen, Z.; Ni, X.; Lu, J.; Wang, J.; Xu, B. Numerical simulation of the ultrasonic waves generated by ring-shaped laser illumination patterns. Opt. Laser Technol. 2007, 39, 1281–1287. [Google Scholar] [CrossRef]
  19. Dixon, S.; Harrison, T.; Fan, Y.; Petcher, P.A. Thermoelastic laser generated ultrasound using a ring source. J. Phys. D Appl. Phys. 2012, 45, 175103. [Google Scholar] [CrossRef]
  20. Pei, C.; Yi, D.; Liu, T.; Kou, X.; Chen, Z. Fully noncontact measurement of inner cracks in thick specimen with fiber-phased-array laser ultrasonic technique. NDT E Int. 2020, 113, 102273. [Google Scholar] [CrossRef]
  21. Feng, W.; Yang, D.; Guo, Y.; Chang, Y. Finite element modeling of bulk ultrasonic waves generated by ring-shaped laser illumination in a diamond anvil cell. Opt. Express 2012, 20, 6429–6438. [Google Scholar] [CrossRef] [PubMed]
  22. Zeng, W.; Cai, M.; Lu, T.; Yu, F.; Yao, F.; Qi, S.; Liu, L. Finite element simulation of dual-laser source generated ultrasonic waves for detection of branched-breaking defects. Optik 2021, 225, 165263. [Google Scholar] [CrossRef]
Figure 1. Directivity patterns of laser-generated ultrasonic waves in iron: (a) longitudinal wave and (b) shear wave in thermoelastic regime.
Figure 1. Directivity patterns of laser-generated ultrasonic waves in iron: (a) longitudinal wave and (b) shear wave in thermoelastic regime.
Sensors 23 00694 g001
Figure 2. Schematic diagram of RLUT system for thickness measurement.
Figure 2. Schematic diagram of RLUT system for thickness measurement.
Sensors 23 00694 g002
Figure 3. Schematic diagram of (a) ring-shaped laser irradiation and cross-sectional view of wave propagation in (b) a single-layer material and (c) a bi-layer material.
Figure 3. Schematic diagram of (a) ring-shaped laser irradiation and cross-sectional view of wave propagation in (b) a single-layer material and (c) a bi-layer material.
Sensors 23 00694 g003
Figure 4. Typical RLUT signal: the ring radius R matches the layer thickness h vs. R unmatches h.
Figure 4. Typical RLUT signal: the ring radius R matches the layer thickness h vs. R unmatches h.
Sensors 23 00694 g004
Figure 5. Schematic diagram of the simulation model: (a) 3D view; (b) half cross-section view of the Cu-W bi-layer plate; (c) half cross-section view of the Al-Fe-Cu three-layer plate. The red zone represents the laser irradiated area.
Figure 5. Schematic diagram of the simulation model: (a) 3D view; (b) half cross-section view of the Cu-W bi-layer plate; (c) half cross-section view of the Al-Fe-Cu three-layer plate. The red zone represents the laser irradiated area.
Sensors 23 00694 g005
Figure 6. Cut plane view of the volumetric ultrasonic displacement field at different times: (a) before bulk waves overlap; (b) longitudinal waves overlap; (c) shear waves overlap; (d) after bulk waves overlap.
Figure 6. Cut plane view of the volumetric ultrasonic displacement field at different times: (a) before bulk waves overlap; (b) longitudinal waves overlap; (c) shear waves overlap; (d) after bulk waves overlap.
Sensors 23 00694 g006
Figure 7. Amplitudes of the shear wave along the central axis near the focus time.
Figure 7. Amplitudes of the shear wave along the central axis near the focus time.
Sensors 23 00694 g007
Figure 8. Volumetric ultrasonic displacement field at different times for R = 410 μm.
Figure 8. Volumetric ultrasonic displacement field at different times for R = 410 μm.
Sensors 23 00694 g008
Figure 9. Typical simulation results of displacement at the ring center for the thickness measurement of (a) 1st layer and (b) 2nd layer.
Figure 9. Typical simulation results of displacement at the ring center for the thickness measurement of (a) 1st layer and (b) 2nd layer.
Sensors 23 00694 g009
Figure 10. Simulation results of displacement at the ring center and comparison of peak-to-peak amplitudes of the focal shear echoes when R varies near (a) 190 μm and (b) 410 μm.
Figure 10. Simulation results of displacement at the ring center and comparison of peak-to-peak amplitudes of the focal shear echoes when R varies near (a) 190 μm and (b) 410 μm.
Sensors 23 00694 g010aSensors 23 00694 g010b
Figure 11. Simulation results of displacement at the ring center and comparison of peak-to-peak amplitudes of the focal shear echoes when R varies near (a) 80 μm, (b) 150 μm, and (c) 180 μm.
Figure 11. Simulation results of displacement at the ring center and comparison of peak-to-peak amplitudes of the focal shear echoes when R varies near (a) 80 μm, (b) 150 μm, and (c) 180 μm.
Sensors 23 00694 g011aSensors 23 00694 g011b
Table 1. Material parameters used in simulation cases.
Table 1. Material parameters used in simulation cases.
LayerThermal Conductive Coefficient
λ (W·m−1·K)
Thermal Capacity
c (J·kg−1·K−1)
Density
ρ (kg·m3)
Young’s Modulus
E (GPa)
Poisson’s Ratio
σ
Thermal Expansion Coefficient
α (K−1)
Copper400.038589601100.351.70 × 10−5
Tungsten175.0132178003600.280.45 × 10−5
Aluminum238.09002700700.331.70 × 10−5
Iron76.244078702000.291.22 × 10−5
Table 2. Layer thickness measurement results.
Table 2. Layer thickness measurement results.
LayerMeasured
Thickness
Real
Thickness
Relative
Error
Copper (120 μm)121.59 μm120.00 μm1.33%
Tungsten80.54 μm80.00 μm0.68%
Aluminum50.19 μm50.00 μm0.38%
Iron50.30 μm50.00 μm0.60%
Copper (30 μm)30.36 μm30.00 μm1.20%
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.

Share and Cite

MDPI and ACS Style

Qiu, J.; Li, Z.; Pei, C.; Luo, G. Simulation of Layer Thickness Measurement in Thin Multi-Layered Material by Variable-Focus Laser Ultrasonic Testing. Sensors 2023, 23, 694. https://doi.org/10.3390/s23020694

AMA Style

Qiu J, Li Z, Pei C, Luo G. Simulation of Layer Thickness Measurement in Thin Multi-Layered Material by Variable-Focus Laser Ultrasonic Testing. Sensors. 2023; 23(2):694. https://doi.org/10.3390/s23020694

Chicago/Turabian Style

Qiu, Jinxing, Zhengying Li, Cuixiang Pei, and Guoqiang Luo. 2023. "Simulation of Layer Thickness Measurement in Thin Multi-Layered Material by Variable-Focus Laser Ultrasonic Testing" Sensors 23, no. 2: 694. https://doi.org/10.3390/s23020694

APA Style

Qiu, J., Li, Z., Pei, C., & Luo, G. (2023). Simulation of Layer Thickness Measurement in Thin Multi-Layered Material by Variable-Focus Laser Ultrasonic Testing. Sensors, 23(2), 694. https://doi.org/10.3390/s23020694

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