1. Introduction
Ultrasonic waves are already used to measure the elastic properties and thickness of solid materials. These properties, such as the Young’s modulus and Poisson’s ratio, can easily be estimated by measuring the propagation velocities of longitudinal or shear ultrasonic waves in isotropic materials [
1]. For such measurements, bulk samples with much larger dimensions than the ultrasonic wavelength are usually used. In the case of thin plates and films two different approaches may be applied. The first one is based on the application of high frequency ultrasonic waves incident to the front surface of a film and the analysis of the reflected signals. The second method uses guided Lamb or surface acoustic waves that are excited in the plate [
2]. The frequency of the used ultrasonic waves as a rule is lower than in the first approach.
Measurement of elastic properties of thin films by high frequency ultrasonic waves was analyzed by a few authors [
3,
4]. For example, measurement of elastic properties and density of thermoplastic thin polymer films such as polyethylene terephthalate, polycaprolactone triol, and polyvinyl butyral was discussed in [
3]. The samples of thin films samples were placed on a solid substrate and immersed in water. Measurements were performed by a pulse-echo method using high frequency 50 MHz plane ultrasonic waves reflected from the thin films. The elastic moduli and density were found solving the inverse problem using the reflected waves. The obtained values of mechanical parameters correspond to the values reported by other researchers. However, measurement of the thickness of such films was not discussed.
Application for measurements of a focused ultrasonic beam was discussed in [
4]. In this paper, the authors analyzed a double focus method based on irradiation of thin metallic plates by a high frequency 57 MHz ultrasonic focused beam. The depth scanning of samples changing the distance between a focused ultrasonic transducer was performed. The focused signals reflected from front and back surfaces of the plate were used for measurements. For a better understanding of the obtained results, the authors proposed the ray model of multi-mode focusing. Experiments were performed in thin stainless steel and aluminium alloy plates immersed in water. The authors claim that the proposed multi-model method can be used for simultaneous sound-velocity and thickness measurement if the density of the material is known, however this topic is not elaborated.
Such methods have several limitations. The wavelength must be much smaller than the thickness of the sample what in the case of thin films requires very high frequencies of ultrasonic waves. Attenuation of ultrasonic waves increases with frequency and it can reduce the accuracy of measurements. Samples must be immersed in water, which in some cases is not allowed. The parameters of the samples are measured in the thickness direction at the selected position; therefore, in order to get a spatial distribution of thickness, the C-scan of the sample is required, etc. [
2].
Some of those limitations can be avoided using the guided Lamb waves instead of bulk waves. The advantage of those waves is that they propagate in the plane of a plate and their propagation characteristics are sensitive to properties of the material through which they propagate [
5,
6,
7,
8]. On the other hand, the Lamb waves have a multimode and dispersive nature.
Nowadays, ultrasonic guided waves are widely used not only for detection of various defects in planar structures [
5,
6,
7], but for material characterization purposes also. Such characterization is directed to determination of some elastic properties and a thickness of the sample under testing [
7,
8,
9,
10,
11,
12].
Determination of elastic properties and thickness of thin plates and films usually is based on measurement of propagation velocities of at least two different guided wave modes, for example A
0 and S
0 modes, and reconstruction of those parameters from the measured velocities [
7,
8,
9,
10,
11,
12,
13,
14]. For reconstruction, the theoretical velocity dispersion curves are adjusted to experimental results, and in such way the elastic constants and/or thickness of the thin planar specimen are obtained.
Propagation of Lamb waves in a thin gold film on a fused silica semi-infinite substrate was analyzed in [
8]. In the presented research, a numerical inversion procedure was used to determine elastic properties of the gold layer. For measurements, a scanning acoustic microscopy method was used. It was found that the determined Poisson’s coefficient and the density of the gold film were slightly lower and the Young’s modulus was higher than that of the bulk gold material.
A method based on fitting of the theoretical Lamb waves dispersion curves to experimental data is presented in [
9]. This paper describes determining of the thickness and two independent elastic constants of aluminium plates of a few millimetres thickness. Lamb wave modes at different frequencies are generated by means of the contact prismatic coupling block method. The visualization of the displacements of the Lamb wave was performed using electronic speckle pattern interferometry. From the out-of-plane displacement maps, the wavelength of a single-mode Lamb wave of the known frequency is measured. Numerous Lamb modes at different frequencies are generated, from which the thickness of the plate and two independent elastic constants are determined.
A similar method was used in [
10]. The presented paper described the combination of the simplex method and the semi-analytic finite element (SAFE) algorithm to determine elastic constant of a given structure. The propagation of Lamb waves is monitored in order to extract the group velocity dispersion curves in an aluminium plate. A pulsed laser excited the guided wave, and a pair of broadband piezoelectric transducers was used to detect the Lamb wave. The reconstruction of the elastic properties of the plate is performed until the discrepancy between the SAFE and the experimental dispersion curves is minimized. The reconstruction used the A
0 and S
0 modes in the frequency range between 100 and 650 kHz.
The hybrid computational system for aluminium plate parameter identification is presented in [
11]. The method is based on guided wave measurement and application of artificial neural networks. The Lamb waves were generated by a contact method or laser. The proposed method is non-iterative and using pseudo-experimental data.
Measurement of the thickness of metallic films with an antisymmetric Lamb A
0 wave mode was discussed in [
12]. The Lamb wave was excited by a piezoelectric element and picked up by an electromagnetic acoustic transducer. Therefore, this method is suitable for measurements of metallic films only. It is also necessary to point out that the velocity of the A
0 mode depends not only on a thickness of the film, but on the Young’s modulus also; however, it was not analyzed in this paper.
In most of the presented papers, the attention is focused on measuring a thickness and elastic properties not of plastic but of aluminium plates. The thicknesses of the plates ranged from 1 to 5 mm, and a laser was used to excite or receive the Lamb waves. Aluminum plates are favorable objects for the study of Lamb waves due to the simple excitation, propagation, and reception of these waves.
The guided Lamb waves were also applied for measurement of elastic properties of thin paper and mineral products [
13,
14,
15]. However, in [
13,
14] the two-side access air-coupled method is used—the ultrasonic wave is sent through an air gap to the paper sample and is picked up on the other side of the sample. The elastic constant such as the Young’s modulus is found from the comparison of theoretical and experimentally obtained dispersion curves. For investigation of a mineral paper, rather high-frequency ultrasonic waves in the range of 0.15–2.3 MHz were used. Such a method is not convenient for industrial applications due to two-side access approach.
The one-side approach for on-line measurement of the tensile stiffness of the moving anisotropic paper web in a paper machine is presented in [
15]. The guided S
0 wave mode is excited by a dry friction and picked up by air-coupled ultrasonic transducers. The tensile stiffness index is found from the measured velocity of the S
0 mode. This method was successfully verified in industrial conditions.
Attenuation of a slow A
0 Lamb wave mode in a polyvinyl chloride (PVC) film was measured by us [
16] using cylindrical guided wave propagating inside the film. Application of such wave enabled reduction of diffraction errors and increasing accuracy of attenuation measurements. However, other parameters such as thickness or elastic moduli were not measured.
From the presented review, it follows that almost all publications are devoted to measurement of elastic properties and/or thickness of relatively thick 0.8–3 mm aluminium plates by using contact or semi-contact methods. There are some publications in which are described measurements of composite materials [
17,
18,
19], but there are no publications devoted to measurements of very thin plastic films such as polyvinylchloride (PVC) films using guided Lamb waves. Production of such films only in Europe is 5 million tons annually [
20]. We have found that such parameters of PVC films as thickness claimed by a manufacturer do not correspond to the real observed values and significant variations of a thickness across even a small size A4 format samples reach up to 10%. Therefore, development of ultrasonic measurement methods suitable for monitoring quality of the production becomes relevant. For industrial applications, the methods with one-side approach and without direct contact with a film would be preferable because they could be applied for on-line monitoring.
It is necessary to point out that those properties of plastic materials, and especially of PVC films, are very different from metals. The thickness of such films is in the range from 25 µm up to 200 µm instead of thicknesses in the range of a few millimeters characteristic for metallic plates. Propagation velocities of ultrasonic waves including A
0 and S
0 guided wave modes in plastic films are significantly lower. For example, the velocity of S
0 mode in aluminium plates is 5400 m/s and in PVC films it is only 1600 m/s. The attenuation of bulk longitudinal ultrasonic waves in plastics increases with the frequency and is higher than in metallic plates. For example, in plastics such as polypropylene and PVDF, they are between 2 dB/cm at 300 kHz and 5 dB/cm at 500 kHz [
21,
22]. In the case of guided waves, the frequency dependences of the attenuation are more complicated. Our measurements have shown that, even at the low frequency of 44 kHz, the attenuation coefficient of A
0 mode in PVC 150 µm thickness film is 2 dB/cm [
16]. All of those differences indicate that the methods developed for characterization of metal plates cannot be applied directly for plastic films.
Therefore, the objective of this research was development of the single side access semi-contactless ultrasonic method suitable for simultaneous measurement of the thickness and Young’s modulus of thin plastic films.
The paper is organized as follows. In
Section 2 theoretical analysis of propagation of A
0 and S
0 guided wave modes in PVC film are analyzed, and a measurement algorithm is proposed. In
Section 3, the experimental set-up with which measurements were performed is described and experimental results are presented. In
Section 4, conclusions and discussion of the obtained results are given.
2. Theoretical Analysis
2.1. Properties of the Measured Parameters
Propagation of guided Lamb waves in thin plates and films depends on elastic properties of the material, density, thickness, and frequency of the guided wave. It allows determination of a thickness and some elastic parameters such as Young’s modulus of the structure under a test exploiting the measured phase velocities of different Lamb wave modes.
In most cases, the density of such structures is known or estimated experimentally, the Poisson’s ratio ν is often known and in a practice the influence of its variations on propagation velocities is negligible. Therefore, the main parameters characterizing such structures are their thickness d, and Young’s modulus E.
For investigation of feasibility to measure, simultaneously the thickness and the Young’s modulus we have selected thin PVC film with the thickness in the range up to 250 µm. The properties of the PVC film provided by the manufacturer (United States Plastics Corporation VYCOM, Scranton, PA, USA) are presented in
Table 1 [
20].
First, we have calculated dispersion curves of guided wave modes by the semi-analytic finite element method (SAFE) [
5]. Calculation by this method in the case of a thin film or plate with infinite lateral dimensions of the plate is based on splitting it into a finite number of
M thin layers, each of which is described in one axis direction. In the second axis direction, it is assumed that the plate is infinite. Therefore, the particle displacement of any point in an element is given by
where
N(
ζ) is the matrix of the shape function,
ζ is the variable in the local coordinate system,
U(e) are the nodal displacements of the element,
x-axis is the wave propagation direction,
k is the wavenumber,
ω is the angular frequency, and
t is the time.
After a standard finite element, assembling procedure the following linear system of algebraic equations in the global coordinate system is obtained
where the
U represents the global vector of nodal displacements,
K1,
K2,
K3, and
M are the matrices in the global coordinate system. The solutions of this equation describe the propagation characteristics of the Lamb wave in a film and the results are presented as dispersion curves.
Analysis of the calculated dispersion curves showed that in the 150 µm thickness film at the frequencies lower than 50 MHz only two Lamb wave modes A
0 and S
0 may propagate. The dispersion curves of those modes in the frequency range 0–300 kHz are shown in
Figure 1.
From the results presented in
Figure 1a, it follows that the phase velocity of the A
0 mode strongly depends both on a thickness
d of the film and the frequency of the ultrasonic signal
f. Contrarily, the phase velocity of the S
0 mode in the low frequency range (<300 kHz) almost does not depend either on the frequency or on the thickness of the film. This conclusion is valid for films of different thickness. Taking that into account, it is possible to assume that the velocity of the S
0 mode does not depend on the thickness of the film, but it should depend on elastic properties of the film, particularly on the Young’s modulus. Correspondingly, the velocity of the A
0 mode should depend on both parameters—the thickness and the Young’s modulus.
It is necessary to point out that the A
0 mode phase velocity’s sensitivity to variations of the
f·d product is biggest in the range of low
f·d values—e.g., less than 20. In the case of the 150 µm thickness film, it corresponds to the frequencies lower than 100 kHz (
Figure 1b). Hence, for precise thickness measurements, the frequency of the A
0 mode should be also selected lower than 100 kHz. It means that the thickness and the Young’s modulus of the film can be found exploiting the measured phase velocities of the S
0 and A
0 modes.
For excitation of A
0 and S
0 modes, different frequencies were selected. For the S0 mode excitation the fixed frequency
f = 180 kHz was selected because the dispersion curve of this mode in the low frequency range (<300 kHz) is independent of the frequency and the same measurement results would be obtained at other frequencies. For excitation of the A
0 mode, the lower frequency
f = 50 kHz was chosen because according to the modeling results at this frequency the best sensitivity to variations of the film thickness is obtained [
23].
To check those conclusions, we have performed calculations of the phase velocities of both modes versus the Young’s modulus and the thickness of the film. The results are presented in
Figure 2a,b.
The calculations were performed at the fixed frequency 50 kHz. In
Figure 2a, dependencies of the phase velocities of the A
0 and S
0 modes for the fixed film thickness
d = 150 µm and for the fixed Young modulus
E = 2.937 GPa (
Figure 2b) are presented. The obtained results show that as expected the phase velocity of the S
0 mode strongly depends on the Young’s modulus, but almost is not affected by the thickness
d of the film. The phase velocity of A
0 mode depends on both parameters, but the thickness of the film influences the phase velocity much more than the phase velocity of the S
0 mode.
Based on those conclusions, it is possible to assume that, by knowing the dispersion curves of the A0 and S0 modes in a certain frequency range, it is possible to estimate the thickness and the Young’s modulus of the film under study.
2.2. Measurement Algorithm
For estimation of a film thickness and the Young’s modulus using velocity measurements, the relations of those parameters to the phase velocities of A
0 and S
0 modes are necessary. It is not possible to get them in an analytic way, but they can be found numerically from the set of dispersion curves calculated for films of different thickness and different Young’s modulus. The relations between the film thickness
d, the Young’s modulus
E and the phase velocities of both modes obtained numerically were approximated by third order polynomials:
where
cphA is the phase velocity of the A
0 mode and
cphS is the phase velocity of the S
0 mode. Both velocities in those approximations are in m/s. The subscripts A and S indicate that the parameter is related to A
0 or S
0 mode.
From the performed calculations, it follows that the phase velocity of the S0 mode does not noticeably depend on the thickness d of the film; therefore, it should be exploited only for determination of the Young’s modulus.
The obtained results allow formulating the following algorithm for estimation of the Young’s modulus E and the film thickness d:
Equation (8) is derived from the well-known formula proposed by Cremer using the resonance method [
24]
This formula is valid for homogeneous thin plates in a low frequency range, which in our case is up to 100 kHz.
For thickness measurements, Equations (3) and (8) were proposed. Equation (8) presented in the described measurement algorithm is preferable because it allows to calculate the film thickness faster and easier. Using Equation (3), new approximation curves should be calculated from the obtained experimental results and only then the film thickness could be estimated.
The proposed algorithm consists of several previously investigated and verified algorithms, which are described in detail and investigated in or publications [
24,
25,
26]. The new algorithm is designed as a combination of those algorithms for solving newly raised problems.
2.3. Guided Wave Excitation and Reception Method
For performing measurements of velocities of A0 and S0 modes the excitation and reception methods of those modes should be selected. For practical applications, for example for on-line measurements during manufacturing process, non-contact, e.g., air-coupled methods should be more preferable.
The excitation method very much depends on selection of the frequency
f of guided waves to be excited (
Figure 3).
In the case of A0 mode, there are two different frequency ranges. The first one is the low frequency range in which velocity of the A0 mode is lower than the ultrasound velocity in air cair = 342 m/s. In this frequency region, the A0 mode does not excite a leaky wave in surrounding air. The second one is the range of higher frequencies in which the velocity of this mode is higher than the ultrasound velocity in air and its propagation is accompanied by a leaky wave in air. Excitation of guided waves in those two different frequency ranges by air-coupled ultrasonic transducers is different.
For excitation and reception of ultra-slow A
0 mode, very efficient lead magnesium niobate-lead titanate (PbMg1/3Nb2/3O3-PbTiO3) PMN-32%PT air-coupled transducers and arrays can be used [
18]. For the film of 150 μm thickness, the phase velocity of the A
0 mode is slower than the ultrasound velocity in air in the frequency range up to 300 kHz (
Figure 3). The highest frequency at which
cphA becomes slower than
cair is increasing when the thickness of the film is decreasing. The lower frequency range in the case of thin films is more attractive because, in it, the sensitivity of the phase velocity to thickness variations is higher than at higher frequencies. Consequently, in this frequency range, it is possible to obtain a higher accuracy of thickness measurements.
At higher frequencies air-coupled transducers deflected with respect to the normal to the surface of the film at the angle
αopt given by the Snell’s law can be exploited
where
λair and
λA0(
f) are wavelengths correspondingly in air and in the film.
The phase velocity of the S0 mode cS0(f) always is higher than the ultrasound velocity in air cair; therefore, it can be excited by an air-coupled transducer deflected according to Snell’s law.
The optimal deflection angles of air-coupled ultrasonic transducers versus the frequency for films of two different thicknesses are shown in
Figure 4.
From the presented results, is follows that—in the low frequency range up to 230 kHz in the film of 200 μm thickness—to excite the A
0 mode by a deflected air-coupled transducer is impossible because the required deflection angle with respect to the normal to the film surfaces becomes 90°, which is impossible to realize. The same conclusion is valid for the film of 150 μm in the frequency range up to 305 kHz. Therefore, for excitation of the A
0 mode in this frequency range, PMN-32%PT arrays can be exploited [
18].
Efficiency of air-coupled excitation of S0 mode is lower than of the A0 mode, especially in the lower frequency range, due to much higher velocity of this mode, the wavelength is much longer and consequently the ratio of the excitation area diameter to the wavelength is much smaller. In such a case, the efficiency of excitation reduces. Therefore, for excitation of the S0 mode, we have selected the frequency range higher than 100 kHz. For this, we propose to use air-coupled flat ultrasonic transducer deflected from the normal to the surface of the film at the angle given by Snell’s formula (Equation (10)).
4. Discussion and Conclusions
The aim of this work was to develop and to investigate a one-side access ultrasonic method suitable for simultaneous measurement of the thickness and the Young’s modulus of thin plastic films. To achieve this goal, the excitation, propagation, and reception of the Lamb waves in thin PVC films were studied. It was found that the methods used for characterizing thin metal plates cannot be directly applied to thin plastic films. This is due to the much lower velocities of the guided Lamb waves propagating in plastic films and a higher attenuation of those waves. In the case of thin plastic films, the phase velocity of the A0 mode exploited for thickness measurements becomes slower than the ultrasound velocity in air and there is no leaky wave radiated by the film into the air. In the case of metallic films, in most cases, phase velocity of A0 mode is higher than in surrounding air and its propagation is accompanied by a leaky wave. That requires different excitation and reception methods of ultrasonic guided Lamb waves.
It was found that the phase velocity of the symmetrical S
0 mode does not depend on the thickness of the film but depends on the Young’s modulus. Contrarily, the phase velocity of the A
0 mode depends both on the thickness and on the Young’s modulus of the film. Exploiting measurements of phase velocities of those two modes, the measurement algorithm for evaluation of the thickness and the Young’s modulus was proposed. It was found that the highest sensitivity of the A
0 mode to variations of the film thickness is observed in the low frequency range up to 100 kHz, in which the phase velocity of this mode is lower than the ultrasound velocity in air. It means that, in this frequency range, it is impossible to excite the A
0 mode by conventional means using deflected air-coupled transducers. In order to overcome this problem, we have proposed using high efficiency PMN-PT piezoelectric crystals, which allows exciting the A
0 mode without contact. To excite S
0 mode in a contactless way is problematic, therefore for this purpose we have used contact point type ultrasonic transducers with a hemispherical plastic tip, which provided a reliable dry acoustic contact with a film and left no traces on the film surface when it was moving. In principle, it is possible to use a multi-element phased array with PMN-PT piezoelectric elements for excitation of the S
0 mode in a low frequency range [
29].
The performed measurements demonstrated feasibility to simultaneously measure the thickness and the Young’s modulus of thin PVC films by means of guided ultrasonic waves. It is necessary to point out that the measured thickness of PVC films significantly differed from the values declared by the manufacturer. The method developed looks promising for a contactless on-line monitoring of the mentioned parameters. In this case, for contactless reception of guided waves, the air-coupled array and corresponding signal processing procedure developed by us and described in [
30] could be used.