Next Article in Journal
Enhanced Soft 3D Reconstruction Method with an Iterative Matching Cost Update Using Object Surface Consensus
Previous Article in Journal
Anomaly Detection of Water Level Using Deep Autoencoder
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

Elastic Properties Measurement Using Guided Acoustic Waves

Institute of Sensor and Actuator Technology (ISAT), Coburg University of Applied Sciences and Arts, Am Hofbräuhaus 1b, 96450 Coburg, Germany
*
Author to whom correspondence should be addressed.
Submission received: 1 September 2021 / Revised: 29 September 2021 / Accepted: 6 October 2021 / Published: 8 October 2021
(This article belongs to the Section Physical Sensors)

Abstract

:
Nondestructive evaluation of elastic properties plays a critical role in condition monitoring of thin structures such as sheets, plates or tubes. Recent research has shown that elastic properties of such structures can be determined with remarkable accuracy by utilizing the dispersive nature of guided acoustic waves propagating in them. However, existing techniques largely require complicated and expensive equipment or involve accurate measurement of an additional quantity, rendering them impractical for industrial use. In this work, we present a new approach that requires only a pair of piezoelectric transducers used to measure the group velocities ratio of fundamental guided wave modes. A numerical model based on the spectral collocation method is used to fit the measured data by solving a bound-constrained nonlinear least squares optimization problem. We verify our approach on both simulated and experimental data and achieve accuracies similar to those reported by other authors. The high accuracy and simple measurement setup of our approach makes it eminently suitable for use in industrial environments.

1. Introduction

The field of guided acoustic waves (GAWs) has attracted increasing interest in the research community over the past decade [1]. In addition to various applications in nondestructive testing [2,3] and structural health monitoring [4,5], guided waves have been successfully used to measure liquid properties [6], monitor biogenic deposits in tubes [7], and even to accelerate electrochemical reactions during the battery charging process [8].
Guided acoustic waves are dispersive, i.e., their phase and group velocities depend on frequency [9]. For particular material and geometry, this dependence is usually represented by dispersion curves, which have proven to be a useful tool for the development of GAW applications [10,11]. Consequently, various methods have been developed to compute dispersion curves for different scenarios, including cylindrical [12,13,14] and multilayer systems [15,16,17], as well as waveguides surrounded by liquids [18,19,20].
Using the opposite approach, the dispersive nature of GAWs can be utilized to determine the elastic properties of thin structures such as sheets, plates, and tubes. The procedure first requires the measurement of group or phase velocities at multiple frequencies. A numerical model is then fitted to these measurements by optimizing the model parameters. Recent research in this area has shown that this strategy can determine the elastic properties of thin structures with remarkable accuracy [21,22,23,24,25]. However, existing studies differ widely in terms of how the velocities are measured.
One of the earliest works on this subject was published by Rogers [21] in the mid-nineties. He used a pair of variable-beam contact transducers to measure phase velocities of different Lamb-type modes, by considering the phase shift over a certain path. Sale et al. [22] used a pair of piezoelectric transducers at a fixed distance to measure the group velocities of zero order guided wave modes excited by a laser. In another publication [23] a speckle interferometer was used to measure the wavelength of multiple guided wave modes, excited using a transducer wedge. A similar approach was reported [24], where a laser-scanning vibrometer was used to measure the wavenumbers of different Lamb-type modes. Gao et al. [25] used a laser ultrasonic setup for both, excitation and detection of different Lamb modes.
Most of the described techniques either require complicated and expensive equipment or involve accurate measurement of an additional quantity, rendering them less suitable for use outside the laboratory environment. In this work, we introduce a new simple and precise approach that requires only a pair of piezoelectric transducers used to measure the group velocities ratio of fundamental Lamb modes. A numerical model based on spectral collocation method is used to fit the measured data, and both longitudinal and transverse wave velocities and the thickness of the inspected structure are estimated by solving a bound-constrained nonlinear least squares optimization problem.

2. Materials and Methods

In the following section, we introduce the measurement setup, the numerical model used to fit the measured data and the optimization approach used to extract the elastic properties of the inspected structure.

2.1. Measurement Setup

The measurement setup involves a pair of identical piezoelectric disc transducers (PIC255, diameter 8 mm, thickness 0.5 mm, PI Ceramic GmbH, Lederhose, Germany), attached to the surface of the test specimen at a fixed distance d = 100 mm, as shown in Figure 1. The transmitting transducer is excited by a signal generator (33500B Series, Keysight Technologies, Inc., Santa Rosa, CA, USA), which provides five-cycle sine bursts modulated by a Hanning window. A digital storage oscilloscope (HDO6034 (Teledyne LeCroy, Chestnut Ridge, NY, USA) is used to record the Lamb-type guided acoustic wave signals arriving at the receiving transducer.
A typical signal recorded by the oscilloscope is demonstrated in Figure 2. In order to extract the propagation times t A 0 and t S 0 of both fundamental Lamb modes, the envelope of the signal is calculated and the time shifts of the maxima relative to the excitation signal are determined. We note that, depending on the application, other approaches such as those based on computation of the wavelet transform [26] or the short-time Fourier transform [27] may also be considered for this purpose.
Once the propagation times have been determined, group velocities c A 0 and c S 0 are usually calculated by considering the distance d between the transmitter and the receiver
c A 0 = d t A 0   and   c S 0 = d t S 0 .
However, accurate measurement of the distance d not only requires additional equipment, but also affects the measurement uncertainty of the group velocities. Our approach overcomes this issue using a simple trick: instead of estimating the velocities directly, we calculate their ratio and thus eliminate the distance d in Equation (1), yielding
y = c S 0 c A 0 = t A 0 t S 0 .

2.2. Numerical Model

A popular approach used to calculate dispersion curves of a plane waveguide, is to solve the analytical frequency equations using root-finding techniques [1]. However, for complicated problems the root-finding process becomes increasingly challenging, particularly if multilayer systems need to be examined [28]. For these scenarios, matrix methods are usually preferred.
Matrix methods compose the individual solutions of each layer into a system matrix and the general solution, which satisfies the boundary and the interface conditions, is then calculated for the entire system. This is done by finding the roots of the characteristic function, which is represented by the determinant of the system matrix (a detailed treatment of matrix methods can be found in [15]). However, matrix methods are prone to numerical errors and are usually associated with high computational costs [29].
In this work, we employ an alternative popular method recently introduced by various authors [28,29,30,31,32,33], which allows the dispersion curves to be calculated by numerically solving the underlying differential equations using the spectral collocation method. This method provides a good trade-off between precision, implementation effort, and computation time, and can be easily applied to cylindrical geometries or multilayer systems [28,30]. However, the spectral collocation method is not always robust. A detailed discussion on the robustness of this method in the calculation of dispersion curves of acoustic waveguides has already been provided by other authors (see, e.g., references [29,30,31]). For applications where robustness is critical, the approach may need to be adapted.

2.2.1. Derivation of Differential Equations

To introduce the numerical model, we consider a planar waveguide as shown in Figure 3. Following [34], we express the displacement field in terms of potentials Φ and Ψ and thus obtain two uncoupled wave equations
2 Φ x 1 2 + 2 Φ x 2 2 = 1 c L 2 2 Φ t 2 , 2 Ψ x 1 2 + 2 Ψ x 2 2 = 1 c T 2 2 Ψ t 2 ,
for longitudinal and transverse waves, which propagate independently in unbounded solid medium with velocities c L and c T , respectively. The coupling of both wave equations is provided by the boundary conditions at the surfaces of the waveguide, requiring both stress components σ 21 and σ 22 to be equal to zero
σ 21 = μ u 2 x 1 + u 1 x 2 = 0 , σ 22 = λ u 1 x 1 + u 2 x 2 + 2 μ u 2 x 2 = 0 ,
where λ and μ are Lamé constants and the displacement components u 1 and u 2 are defined as follows
u 1 = Φ x 1 + Ψ x 2 , u 2 = Φ x 2 Ψ x 1 .
For solutions of the wave equations we consider plane harmonic waves of the type
Φ ( x 1 , x 2 , t ) = ϕ ( x 2 ) e i ( k x 1 ω t ) , Ψ ( x 1 , x 2 , t ) = i ψ ( x 2 ) e i ( k x 1 ω t ) ,
where k is the wavenumber, t the time, ω the angular frequency and i the imaginary number. By substituting (6) into wave Equation (3) we obtain
2 x 2 2 k 2 ϕ = ω 2 c L 2 ϕ , 2 x 2 2 k 2 i ψ = ω 2 c T 2 i ψ .
By analogy, substituting (6) into (4) delivers
i σ 21 μ = 2 k x 2 ϕ + 2 x 2 2 + k 2 i ψ = 0 , σ 22 2 μ = ( α + 1 ) 2 x 2 2 α k 2 ϕ k x 2 i ψ = 0 ,
where
α = λ 2 μ 1 2 c L 2 c T 2 1 .

2.2.2. Spectral Collocation Method

The spectral collocation method (SCM) is a numerical approach to solving differential equations. Its main idea is to expand the solution of the differential equation by a linear combination of some global polynomial basis functions that satisfy the exact solution on a set of N pre-assigned collocation points. This allows the evaluation of derivatives at such collocation points by multiplication with a pre-computed matrix. For bounded non-periodic domains, Chebyshev differentiation matrices are typically used [35].
Following this scheme, we replace the m-th order differential operators in Equation (7) with the corresponding Chebyshev differentiation matrices D ( m ) (see, e.g., [36] for implementation details) and thus obtain a system of equations
L ϕ i ψ = ω 2 M ϕ i ψ ,
where
L = D ( 2 ) k 2 · I 0 0 D ( 2 ) k 2 · I ,
M = 1 / c L 2 · I 0 0 1 / c T 2 · I
and I is an N × N identity matrix.
Similarly, from Equation (8) it follows that
S ϕ i ψ = 0 ,
where
S = 2 k D ( 1 ) D ( 2 ) + k 2 · I ( α + 1 ) D ( 2 ) α k 2 · I k D ( 1 ) .
At this point, the wave equations and boundary conditions are both represented by a set of 2 N equations corresponding to N Chebyshev collocation points defined on the interval [ 1 , 1 ] . Using coordinate transformation
D h 2 , h 2 ( m ) = 2 h m D [ 1 , 1 ] ( m ) ,
the solution is introduced on the problem domain.
Finally, (10) and (13) are combined into a system of equations by replacing the rows 1, N, N + 1 and 2 N of the matrix [ L ] with the corresponding rows of the matrix [ S ] and setting these rows to zero in the matrix [ M ] .
Once the eigenvalues ω 2 are calculated for different k values by solving the matrix eigenvalue problem (10) using the QZ algorithm [37], the group velocities are calculated as follows
c = d ω d k .
The above equations combine to give the numerical model M , which is used to fit the measured data, as described in the next section.

2.3. Optimization

In our approach, we use the numerical model described in the previous section to fit our measured data. The elastic properties and thickness of the test specimen are determined by solving an optimization problem that can be formulated as follows: Consider n measurements y i acquired at frequencies f i such that f i f j i j and i , j { 1 , 2 , , n } . We wish to describe the relationship between the frequency f and the measured quantity y using a numerical model M ( f ; ϑ ) whose parameters form vector ϑ = ( c L , c T , h ) T . As a measure for the discrepancy between our model and the observed data, we define the residuals
r i ( ϑ ) = y i M ( f i ; ϑ )
and finally obtain a nonlinear least squares problem
min ϑ i = 1 n r i 2 ( ϑ ) ,
which can be solved using standard techniques such as the Levenberg–Marquardt algorithm [38].
From a practical point of view, it may be useful to define a search region for the model parameters to prevent the algorithm from getting stuck in a local optimum. In such cases, the task is transformed into a bound-constrained optimization problem, which can be solved efficiently using the STIR (Subspace Trust Region Interior Reflective) method [39].

3. Results

3.1. Simulated Data

To verify our approach, we first applied it to simulated data. Guided acoustic wave signals in the frequency range 110 to 200 k Hz were provided by a finite element model of an h = 1 mm thick aluminium sheet. The elastic properties of aluminium ( c L = 6.35 km / s , c T = 3.10 km / s , ρ = 2.8 g / c m 3 ) were chosen according to the literature [1]. According to Adamou [30], the number of collocation points N required to determine the group velocities of j lowest Lamb modes with an accuracy of 10 digits can be calculated with the formula N 2 j + 10 . For our numerical model, we used a grid of N = 16 collocation points, which proved to be a good trade-off between accuracy and computation time. The parameter bounds were chosen as follows: 5.00 km / s   c L 7.00 km / s , 2.50 km / s   c T 3.50 km / s , and 0.90 mm d 1.10 mm.
The group velocities ratios derived from the simulated acoustic signals at different frequencies as well as the corresponding solution found by optimization are shown in Figure 4. The good agreement between the data and the solution indicates that the algorithm was able to find an optimum. In fact, the relative error between the determined parameters and the reference was about 1%, as can be seen in Table 1.
Since optimization results typically depend on the initial solution, or in our case on the search region specified for the solution, we repeated the optimization process multiple times with some randomly chosen upper and lower parameter bounds. In each case, the relative error between the determined parameters and the reference was about 1%.

3.2. Experimental Data

We also applied our approach to experimental data. For this purpose, an aluminium sheet of h = 1 mm thickness was investigated. A pair of piezoelectric disc transducers were attached to its surface as described in Section 2.1. The excitation signals in the frequency range 110 to 200 k Hz were provided by a signal generator and the received signals were recorded using a digital storage oscilloscope. For the numerical model we used a grid of N = 16 collocation points and the search region for the solution was set in the same way as it was for the simulated data.
To verify the experimental results, additional reference measurements were performed. The thickness of the aluminium sheet was measured using a micrometer screw gauge and the longitudinal wave velocity was obtained by time-of-flight measurements with an ultrasonic transducer (M2017 Delay Line Transducer, Olympus Europa SE & Co. KG, Hamburg, Germany) operating at 30 M Hz .
The group velocities ratios obtained from the measured acoustic signals and the corresponding solution are shown in Figure 4. Compared to the simulated data, some measurement points deviate slightly from the solution found. However, the obtained parameters still show good agreement with the reference, as can be seen in Table 2. Although the relative error is slightly higher in this case, the result can still be termed accurate.

4. Discussion

This paper presents a method for measuring the elastic properties of thin structures such as tubes, plates and sheets, that takes advantage of the dispersive nature of guided acoustic waves. The measurement setup consists of a signal generator, a digital storage oscilloscope and a pair of piezoelectric transducers, which are used to measure the ratio of group velocities at different frequencies. The measured data is then fed into a numerical model based on the spectral collocation method in order to determine the elastic properties and the thickness of the investigated structure via optimization. Validation of the method with simulated and experimental data resulted in a measurement error of 1.1% and 3.4%, respectively, which is similar to the results reported by other authors [21,22,23,24,25].
To determine the wave velocities from the signals measured with a piezoelectric receiver, the distance between the transmitter and the receiver must be known exactly. This requires additional equipment, complicates the measurement setup and raises the possibility of additional measurement uncertainties, which may have a negative impact on the accuracy of the measurement principle. Our approach overcomes this issue by determining the ratio of the group velocities, which can be calculated directly from the propagation times of the two fundamental Lamb modes, significantly simplifying the entire measurement setup.
Laser ultrasonic methods for detection of guided acoustic waves usually rely on the measurement of the wave number rather than the wave velocity. This also avoids the need to determine the distance between the transmitter and receiver. However, such methods require expensive and complicated equipment and are poorly suited for the characterisation of curved or weakly reflective materials. Our approach, on the other hand, uses low-cost components and allows the characterisation of arbitrarily shaped structures and various materials, opening up a wide field of potential applications.
For the determination of the propagation velocity of the two fundamental Lamb modes, a five-cycle sine burst modulated by a Hanning window is used. On the one hand, the propagation time of the wave can be determined relatively accurately using such a compact signal. On the other hand, such a signal has a much broader spectrum than signals with a higher number of cycles, which means that unwanted side effects such as dispersion become more prominent. In the application studied in this article, the Lamb waves are excited at frequencies with a relatively flat course of the frequency-dependent phase velocity, so that dispersion effects only play a minor role. The relatively small distance between transmitter and receiver also has a positive impact. However, depending on the application, such effects can have a negative impact on the measurement method, which is why the use of alternative signal patterns should be investigated in the future.
For future investigations, the method could be adapted for rotationally symmetrical structures such as pipes. In this case, the measurement setup could be reduced to a single transducer that can act as both a transmitter and a receiver of circumferential Lamb waves. In this way, the measurement of wear or deposits in pipes could be provided, enabling continuous online monitoring of industrial plants. However, in contrast to many non-destructive tests performed on pipes and sheet material, the approach here is not aimed at detecting defects but at determining material properties. A defect-free part of the structure must always be selected so that the measurements are not affected by these non-representative inhomogeneities. In this context, weld metal must also be considered a defect, as the homogeneity of the material is no longer given and the evaluation of the measurement results can lead to misleading values.

Author Contributions

Conceptualization, V.F. and F.B.; methodology, V.F. and F.B.; software, V.F. and A.B.; validation, V.F., A.B.; formal analysis, V.F.; investigation, V.F.; resources, A.B.; data curation, V.F.; writing—original draft preparation, V.F.; writing—review and editing, K.S.D.; visualization, V.F.; supervision, K.S.D.; project administration, K.S.D.; funding acquisition, K.S.D. All authors have read and agreed to the published version of the manuscript.

Funding

This research was co-financed by the European Regional Development Fund (EFRE) as part of the “Innovationsterminal für Unternehmen in der ehemaligen nordbayerischen Grenzregion (InnoTerm)” project.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Rose, J.L. Ultrasonic Guided Waves in Solid Media, 2nd ed.; Cambridge University Press: New York, NY, USA, 2014. [Google Scholar]
  2. Su, Z.; Ye, L.; Lu, Y. Guided Lamb Waves for Identification of Damage in Composite Structures: A Review. J. Sound Vib. 2006, 295, 753–780. [Google Scholar] [CrossRef]
  3. Workman, G.L.; Kishoni, D.; Moore, P.O. Ultrasonic Testing, 3rd ed.; American Society for Nondestructive Testing: Columbus, OH, USA, 2007. [Google Scholar]
  4. Raghavan, A.; Cesnik, C.E.S. Review of Guided-Wave Structural Health Monitoring. Shock Vib. Digest 2007, 39, 91–114. [Google Scholar] [CrossRef]
  5. Croxford, A.J.; Wilcox, P.D.; Drinkwater, B.W.; Konstantinidis, G. Strategies for Guided-Wave Structural Health Monitoring. Proc. R. Soc. 2007, 463, 2961–2981. [Google Scholar] [CrossRef]
  6. Lindner, G. Sensors and Actuators Based on Surface Acoustic Waves Propagating Along Solid–Liquid Interfaces. J. Phys. D: Appl. Phys. 2008, 41, 123002. [Google Scholar] [CrossRef]
  7. Tietze, S.; Singer, F.; Lasota, S.; Ebert, S.; Landskron, J.; Schwuchow, K.; Drese, K.S.; Lindner, G. Monitoring of Soft Deposition Layers in Liquid-Filled Tubes with Guided Acoustic Waves Excited by Clamp-on Transducers. Sensors 2018, 18, 526. [Google Scholar] [CrossRef] [Green Version]
  8. Tietze, S.; Schlemmer, J.; Lindner, G. Influence of surface acoustic waves induced acoustic streaming on the kinetics of electrochemical. Micro/Nano Mater. Devices Syst. 2013, 8923, 89231B. [Google Scholar]
  9. Redissi, A.; Miller, S. Experimental Characterization of the Propagation of Guided Acoustic Waves in Pipe Strings. J. Acoust. Soc. Am. 2018, 143, 1385–1391. [Google Scholar] [CrossRef]
  10. Groth, E.B.; Iturrioz, I.; Clarke, T.G.R. The Dispersion Curve Applied in Guided Wave Propagation in Prismatic Rods. Lat. Am. J. Solids Struct. 2018, 15. [Google Scholar] [CrossRef]
  11. Groth, E.B.; Clarke, T.G.R.; Schumacher da Silva, G.; Iturrioz, I.; Lacidogna, G. The Elastic Wave Propagation in Rectangular Waveguide Structure: Determination of Dispersion Curves and Their Application in Nondestructive Techniques. Appl. Sci. 2020, 10, 4401. [Google Scholar] [CrossRef]
  12. Barshinger, J.N.; Rose, J.L. Guided Wave Propagation in an Elastic Hollow Cylinder Coated with a Viscoelastic Material. IEEE Trans. Ultrasonics Ferroelectr. Frequency Control 2004, 51, 1547–1556. [Google Scholar] [CrossRef]
  13. Liu, Y.-C.; Hwang, Y.-F.; Huang, J.-H. Dispersion Relations and Modal Patterns of Wave in a Cylindrical Shell. In Wave Processes in Classical and New Solids; Giovine, P., Ed.; InTech: Rijeka, Croatia, 2012. [Google Scholar]
  14. Sarkar, A.; Sonti, V.R. Simplified Dispersion Curves for Circular Cylindrical Shells Using Shallow Shell Theory. J. Sound Vib. 2009, 322, 1–7. [Google Scholar] [CrossRef]
  15. Lowe, M.J.S. Matrix Techniques for Modeling Ultrasonic Waves in Multilayered Media. IEEE Trans. Ultrasonics Ferroelectr. Frequency Control 1995, 42, 525–542. [Google Scholar] [CrossRef]
  16. Maghsoodi, A.; Ohadi, A.; Sadighi, M. Calculation of Wave Dispersion Curves in Multilayered Composite-Metal Plates. Shock Vib. 2014, 2014, 1–6. [Google Scholar] [CrossRef]
  17. Vaziri Astaneh, A.; Guddati, M.N. Efficient Computation of Dispersion Curves for Multilayered Waveguides and Half-Spaces. Comput. Methods Appl. Mechan. Eng. 2016, 300, 27–46. [Google Scholar] [CrossRef] [Green Version]
  18. Shao, G.-Z.; Li, Q.-C.; Liang, Z.-Q. A Study on Dispersion Curves of Guided Wave in Layered Media with Overlying Liquid Surface. Chin. J. Geophys. 2007, 50, 783–789. [Google Scholar] [CrossRef]
  19. Bao, X.L.; Franklin, H.; Raju, P.K.; Überall, H. The Splitting of Dispersion Curves for Plates Fluid-Loaded on Both Sides. J. Acoust. Soc. Am. 1997, 102, 1246–1248. [Google Scholar] [CrossRef]
  20. Maess, M.; Wagner, N.; Gaul, L. Dispersion Curves of Fluid Filled Elastic Pipes by Standard FE Models and Eigenpath Analysis. J. Sound Vib. 2006, 296, 264–276. [Google Scholar] [CrossRef]
  21. Rogers, W.P. Elastic Property Measurement Using Rayleigh-Lamb Waves. Res. Nondestruct. Eval. 1995, 6, 185–208. [Google Scholar] [CrossRef]
  22. Sale, M.; Rizzo, P.; Marzani, A. Guided Waves Based Approach for the Reconstruction of the Elastic Moduli of Plates. In Proceedings of the 2009 IEEE International Ultrasonics Symposium, Rome, Italy, 20–23 September 2009; pp. 1499–1502. [Google Scholar]
  23. Deán, J.L.; Trillo, C.; Doval, A.F.; Fernández, J.L. Determination of Thickness and Elastic Constants of Aluminum Plates from Full-Field Wavelength Measurements of Single-Mode Narrowband Lamb Waves. J. Acoust. Soc. Am. 2008, 124, 1477–1489. [Google Scholar] [CrossRef] [PubMed]
  24. Ponschab, M.; Kiefer, D.A.; Rupitsch, S.J. Simulation-Based Characterization of Mechanical Parameters and Thickness of Homogeneous Plates Using Guided Waves. IEEE Trans. Ultrasonics Ferroelectr. Frequency Control 2019, 66, 1898–1905. [Google Scholar] [CrossRef] [PubMed]
  25. Gao, W.; Glorieux, C.; Thoen, J. Laser Ultrasonic Study of Lamb Waves: Determination of the Thickness and Velocities of a Thin Plate. Int. J. Eng. Sci. 2003, 41, 219–228. [Google Scholar] [CrossRef]
  26. Lašová, S.; Zemčík, R. Determination of Group Velocity of Propagation of Lamb Waves in Aluminium Plate Using Piezoelectric Transducers. Applied and Computational Mechanics 2017, 11. [Google Scholar] [CrossRef]
  27. Pei, N.; Bond, L.J. Higher Order Acoustoelastic Lamb Wave Propagation in Stressed Plates. J. Acoust. Soc. Am. 2016, 140, 3834. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Karpfinger, F.; Valero, H.-P.; Gurevich, B.; Bakulin, A.; Sinha, B. Spectral-Method Algorithm for Modeling Dispersion of Acoustic Modes in Elastic Cylindrical Structures. Geophysics 2010, 75, 19–27. [Google Scholar] [CrossRef]
  29. Quintanilla, F.H.; Lowe, M.J.S.; Craster, R.V. Modeling Guided Elastic Waves in Generally Anisotropic Media Using a Spectral Collocation Method. J. Acoust. Soc. Am. 2015, 137, 1180–1194. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Adamou, A.T.I.; Craster, R.V. Spectral Methods for Modelling Guided Waves in Elastic Media. J. Acoust. Soc. Am. 2004, 116, 1524–1535. [Google Scholar] [CrossRef]
  31. Lan, B. Non-Iterative, Stable Analysis of Surface Acoustic Waves in Anisotropic Piezoelectric Multilayers Using Spectral Collocation Method. J. Sound Vib. 2018, 433, 16–28. [Google Scholar] [CrossRef]
  32. Quintanilla, F.H.; Fan, Z.; Lowe, M.J.S.; Craster, R.V. Guided Waves’ Dispersion Curves in Anisotropic Viscoelastic Single- and Multi-Layered Media. Proc. R. Soc. A: Math. Phys. Eng. Sci. 2015, 471, 20150268. [Google Scholar] [CrossRef]
  33. Zharnikov, T.; Syresin, D.; Hsu, C.-J. Application of the Spectral Method for Computation of the Spectrum of Anisotropic Waveguides. J. Acoust. Soc. Am. 2013, 133, 3456. [Google Scholar] [CrossRef]
  34. Achenbach, J.D. Wave Propagation in Elastic Solids, 8th ed.; Elsevier: Amsterdam, The Netherlands, 1999. [Google Scholar]
  35. Boyd, J.P. Chebyshev & Fourier Spectral Methods; Springer: Berlin, Germany, 1989. [Google Scholar]
  36. Trefethen, L.N. Spectral Methods in MATLAB; Society for Industrial and Applied Mathematics: New York, NY, USA, 2000. [Google Scholar]
  37. Moler, C.B.; Stewart, G.W. An Algorithm for Generalized Matrix Eigenvalue Problems. SIAM J. Numer. Anal. 1973, 10, 241–256. [Google Scholar] [CrossRef]
  38. Marquardt, D.W. An Algorithm for Least-Squares Estimation of Nonlinear Parameters. J. Soc. Ind. Appl. Math. 1963, 11, 431–441. [Google Scholar] [CrossRef]
  39. Branch, M.A.; Cloeman, T.F.; Li, Y. A Subspace, Interior, and Conjugate Gradient Method for Large-Scale Bound-Constrained Minimization Problems. SIAM J. Sci. Comput. 1999, 21, 1–23. [Google Scholar] [CrossRef]
Figure 1. A schematic diagram of the measurement setup, consisting of a signal generator (GEN), a digital storage oscilloscope (DSO) and two transducers, one for transmitting (T) and one for receiving (R) the guided acoustic wave signals, attached to the surface of the test specimen at a fixed distance (d).
Figure 1. A schematic diagram of the measurement setup, consisting of a signal generator (GEN), a digital storage oscilloscope (DSO) and two transducers, one for transmitting (T) and one for receiving (R) the guided acoustic wave signals, attached to the surface of the test specimen at a fixed distance (d).
Sensors 21 06675 g001
Figure 2. Determination of propagation times of the fundamental Lamb modes.
Figure 2. Determination of propagation times of the fundamental Lamb modes.
Sensors 21 06675 g002
Figure 3. Geometry of the waveguide.
Figure 3. Geometry of the waveguide.
Sensors 21 06675 g003
Figure 4. Group velocities ratios (circles) derived from the acoustic signals at different frequencies and the corresponding fitted model (solid line) for (a) simulated data and (b) experimental data.
Figure 4. Group velocities ratios (circles) derived from the acoustic signals at different frequencies and the corresponding fitted model (solid line) for (a) simulated data and (b) experimental data.
Sensors 21 06675 g004
Table 1. Results for simulated data.
Table 1. Results for simulated data.
ParameterReferenceSolutionError (%)
c L ( m / s )6350.006418.441.08
c T ( m / s )3100.003121.980.71
h ()1.001.010.90
Table 2. Results for experimental data.
Table 2. Results for experimental data.
ParameterReferenceSolutionError (%)
c L ( m / s )6718.756876.132.34
c T ( m / s )-3151.11-
h (mm)1.000.973.40
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Fairuschin, V.; Brand, F.; Backer, A.; Drese, K.S. Elastic Properties Measurement Using Guided Acoustic Waves. Sensors 2021, 21, 6675. https://0-doi-org.brum.beds.ac.uk/10.3390/s21196675

AMA Style

Fairuschin V, Brand F, Backer A, Drese KS. Elastic Properties Measurement Using Guided Acoustic Waves. Sensors. 2021; 21(19):6675. https://0-doi-org.brum.beds.ac.uk/10.3390/s21196675

Chicago/Turabian Style

Fairuschin, Viktor, Felix Brand, Alexander Backer, and Klaus Stefan Drese. 2021. "Elastic Properties Measurement Using Guided Acoustic Waves" Sensors 21, no. 19: 6675. https://0-doi-org.brum.beds.ac.uk/10.3390/s21196675

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