DAMAGE DETECTION IN PLATE STRUCTURES USING SPARSE ULTRASONIC TRANSDUCER ARRAYS AND ACOUSTIC WAVEFIELD IMAGING

Similar documents
Multi-Mode and Multi-Frequency Differential Lamb Wave Imaging with in situ Sparse Transducer Arrays

Detection, localization and characterization of damage in plates with an in situ array of spatially distributed ultrasonic sensors

Ultrasonic structural health monitoring: Strategies, issues and progress

ULTRASONIC SIGNAL CHARACTERIZATIONS OF FLAT-BOTTOM HOLES IN

Piezoelectric Wafer Active Sensor Guided Wave Imaging

Co-Located Triangulation for Damage Position

EMBEDDED NON-DESTRUCTIVE EVALUATION FOR DAMAGE DETECTION USING PIEZOELECTRIC WAFER ACTIVE SENSORS

A NOVEL HIGH SPEED, HIGH RESOLUTION, ULTRASOUND IMAGING SYSTEM

NONDESTRUCTIVE EVALUATION OF CLOSED CRACKS USING AN ULTRASONIC TRANSIT TIMING METHOD J. Takatsubo 1, H. Tsuda 1, B. Wang 1

Probability of Detection Assessment of a Guided Wave Structural Health Monitoring System

FATIGUE CRACK DETECTION IN METALLIC MEMBERS USING SPECTRAL

Instantaneous Baseline Damage Detection using a Low Power Guided Waves System

Ultrasonic Guided Wave Testing of Cylindrical Bars

The Development of Laser Ultrasonic Visualization Equipment and its Application in Nondestructive Inspection

ULTRASONIC GUIDED WAVE ANNULAR ARRAY TRANSDUCERS FOR STRUCTURAL HEALTH MONITORING

Properties of Interdigital Transducers for Lamb-Wave Based SHM Systems

CIRCULAR LAMB AND LINEAR SHEAR HORIZONTAL GUIDED WAVE ARRAYS FOR STRUCTURAL HEALTH MONITORING

Title: Reference-free Structural Health Monitoring for Detecting Delamination in Composite Plates

Time Reversal FEM Modelling in Thin Aluminium Plates for Defects Detection

MEASUREMENT OF RAYLEIGH WAVE ATTENUATION IN GRANITE USING

LASER GENERATION AND DETECTION OF SURFACE ACOUSTIC WAVES

Investigation on Sensor Fault Effects of Piezoelectric Transducers on Wave Propagation and Impedance Measurements

AFRL-RX-WP-TP

Keywords: Ultrasonic Testing (UT), Air-coupled, Contact-free, Bond, Weld, Composites

NUMERICAL MODELING OF AIR-COUPLED ULTRASOUND WITH EFIT. D. E. Chimenti Center of Nondestructive Evaluation Iowa State University Ames, Iowa, USA

ULTRASONIC METHODS FOR DETECTION OF MICRO POROSITY IN COMPOSITE MATERIALS

Measurement of phase velocity dispersion curves and group velocities in a plate using leaky Lamb waves

Ultrasonic Guided Waves for NDT and SHM

DETECTION OF CORROSION IN BOTTOM PLATES OF GAS AND OIL TANKS USING GUIDED ULTRASONIC WAVES AND ELECTROMAGNETIC ULTRASONIC (EMAT) TRANSDUCERS

ENHANCEMENT OF SYNTHETIC APERTURE FOCUSING TECHNIQUE (SAFT) BY ADVANCED SIGNAL PROCESSING

ACOUSTO-ULTRASONIC EVALUATION OF HYBRID COMPOSITES USING

Generation Laser Scanning Method for Visualizing Ultrasonic Waves Propagating on a 3-D Object

RapidScan II Application Note General Composite Scanning

Guided Wave Travel Time Tomography for Bends

Hybrid Passive/Active Impact Detection & Localization for Aerospace Structures Seth S. Kessler and Eric B. Flynn. Metis Design Corporation

RELIABILITY OF GUIDED WAVE ULTRASONIC TESTING. Dr. Mark EVANS and Dr. Thomas VOGT Guided Ultrasonics Ltd. Nottingham, UK

A Numerical study on proper mode and frequency selection for riveted lap joints inspection using Lamb waves.

Basic functions of the universal flaw detector GEKKO

Development of Automatic Neural Network Classifier of Defects Detected by Ultrasonic Means

THE ANALYSIS OF ADHESIVE BONDS USING ELECfROMAGNETIC

Research Article An Investigation of Structural Damage Location Based on Ultrasonic Excitation-Fiber Bragg Grating Detection

Keywords: Guided wave, structural health monitoring, HCSS, disbond, damage index. More Info at Open Access Database

Validation of a Lamb Wave-Based Structural Health Monitoring System for Aircraft Applications

A COMPARISON BETWEEN ASTM E588 AND SEP 1927 RELATING RESOLUTION LIMITS AT DETERMINATION OF THE PURITY GRADE

Lamb Wave Dispersion Compensation in Piezoelectric Wafer Active Sensor Phased-Array Applications

CRACK SIZING USING A NEURAL NETWORK CLASSIFIER TRAINED WITH DATA OBTAINED FROM FINI1E ELEMENT MODELS

DETECTION AND SIZING OF SHORT FATIGUE CRACKS EMANATING FROM RIVET HOLES O. Kwon 1 and J.C. Kim 1 1 Inha University, Inchon, Korea

Structural Integrity Monitoring using Guided Ultrasonic Waves

RECENT ADVANCEMENTS IN THE APPLICATION OF EMATS TO NDE

Long Range Ultrasonic Testing - Case Studies

A New Lamb-Wave Based NDT System for Detection and Identification of Defects in Composites

SPARSE ARRAY TOMOGRAPHY SYSTEM FOR CORROSION EXTENT MONITORING H. Bian, H. Gao, J. Rose Pennsylvania State University, University Park, PA, USA

Frequency Considerations in Air-Coupled Ultrasonic Inspection.

DISBOND DETECTION AND CHARACTERIZATION USING HORIZONT ALL Y

REFLECTION AND TRANSMISSION OF LAMB WAVES AT DISCONTINUITY IN PLATE Z. Liu NDT Systems & Services AG, Stutensee, Germany

EFFECT OF SURFACE COATINGS ON GENERATION OF LASER BASED ULTRASOUND

Detection of Cracks at Rivet Holes in Thin Plates Using Lamb-Wave Scanning

Identifying Scatter Targets in 2D Space using In Situ Phased Arrays for Guided Wave Structural Health Monitoring

Lamb Wave Ultrasonic Stylus

ACOUSTIC MICROSCOPY INSPECTION OF GLASS REPAIR TECHNIQUES

Quasi-Rayleigh Waves in Butt-Welded Thick Steel Plate

MEASUREMENT OF SURFACE ACOUSTIC WAVE USING AIR COUPLED TRANSDUCER AND LASER DOPPLER VIBROMETER

vibro-acoustic modulation

New Multi-Technology In-Line Inspection Tool For The Quantitative Wall Thickness Measurement Of Gas Pipelines

INVESTIGATION OF IMPACT DAMAGE OF CARBON FIBER- RAINFORCED PLASTIC (CFRP) BY EDDY CURRENT NON- DESTRUCTIVE TESTING

FATIGUE CRACK GROWTH MONITORING OF AN ALUMINUM JOINT STRUCTURE

Development of Magnetostriction Based Ultrasonic Transducer. For In-situ High Temperature Inspection

NDT-PRO Services expands service offering

Capabilities of Flip Chip Defects Inspection Method by Using Laser Techniques

Identification of Delamination Damages in Concrete Structures Using Impact Response of Delaminated Concrete Section

Laser-Generation Based Imaging of Ultrasonic Wave Propagation on Welded Steel Plates and Its Application to Defect Detection

LAMB WA VB TOMOGRAPHY USING LASER-BASED ULTRASONICS

GUIDED WAVES FOR DAMAGE MONITORING IN PLATES FOR NOTCH DEFECTS

Transducer degradation and high amplitude behavior of broadband piezoelectric stack transducer for vibrothermography

Delamination Detection in Composite Structures using Laser Vibrometer Measurement of Lamb Waves

ON THE DEVELOPMENT OF METHODS AND TECHNIQUES FOR AIRCRAFT STRUCTURAL HEALTH MONITORING

BINDT Telford. Guided Wave Testing and Monitoring Over Long and Short Ranges

REAL-TIME B-SCAN ULTRASONIC IMAGING USING A DIGITAL PHASED. Robert Dunki-Jacobs and Lewis Thomas General Electric Company Schenectady, New York, 12301

EVOLUTION OF THE CRYOGENIC EDDY CURRENT MICROPROBE

Automation of data collection for PWAS-based structural health monitoring

Detection of Internal OR External Pits from Inside OR Outside a tube with New Technology (EMIT)

ON FIBER DIRECTION AND POROSITY CONTENT USING ULTRASONIC PITCH-CATCH TECHNIQUE IN CFRP COMPOSITE SOLID LAMINATES

Spectral Distance Amplitude Control for Ultrasonic Inspection of Composite Components

Ultrasonic Detection of Inclusion Type Defect in a Composite Panel Using Shannon Entropy

FATIGUE CRACK CHARACTERIZATION IN CONDUCTING SHEETS BY NON

AUTOMATED EDDY CURRENT DETECTION OF FLAWS IN SHOT-PEENED

Guided wave based material characterisation of thin plates using a very high frequency focused PVDF transducer

ACOUSTIC MICRO IMAGING ANALYSIS METHODS FOR 3D PACKAGES

FIRST MEASUREMENTS FROM A NEW BROADBAND VIBROTHERMOGRAPHY MEASUREMENT SYSTEM

Pipeline Technology Conference 2010

Detecting 1 st and 2 nd Layer Simulated Cracks in Aircraft Wing Spanwise Splice Standards Using Remote-Field Eddy Current Technique

ULTRASONIC TECHNIQUES TO QUANTIFY MATERIAL DEGRADATION IN

Testing of Buried Pipelines Using Guided Waves

Ultrasonic pulse propagation in a bonded three-layered structure

Laser-Based Guided Wave Propagation and Mode Decomposition in Detecting the Integrity of Structural I-Beams

Non-Destructive Method Based on Rayleigh-Like Waves to Detect Corrosion Thinning on Non- Accessible Areas

High Frequency Ultrasonic Systems with Frequency Ranges of 35 to 200 MHz

Equipment for Attenuation and velocity of ultrasound in solid state materials (transmission), experimental set-up

KAERI Feeder Tube Inspection Using EMAT Generated Circumferential Guided Waves

In-Situ Damage Detection of Composites Structures using Lamb Wave Methods

Transcription:

DAMAGE DETECTION IN PLATE STRUCTURES USING SPARSE ULTRASONIC TRANSDUCER ARRAYS AND ACOUSTIC WAVEFIELD IMAGING T. E. Michaels 1,,J.E.Michaels 1,B.Mi 1 and M. Ruzzene 1 School of Electrical and Computer Engineering, Georgia Institute of Technology, Atlanta, GA - G.W.Woodruff School of Mechanical Engineering, Georgia Institute of Technology, Atlanta, GA - School of Aerospace Engineering, Georgia Institute of Technology, Atlanta, Atlanta, GA -1 ABSTRACT. A methodology is presented for health monitoring and subsequent inspection of critical structures. Algorithms have been developed to detect and approximately locate damaged regions by analyzing signals recorded from a permanently mounted, sparse array of transducers. Followup inspections of suspected flaw locations are performed using a dual transducer ultrasonic approach where a permanently mounted transducer is the source and an externally scanned transducer is the receiver. Scan results are presented as snapshots of the propagating ultrasonic wavefield radiating out from the attached transducers. This method, referred to here as Acoustic Wavefield Imaging (AWI), provides an excellent visual representation of the interaction of propagating ultrasonic waves with the structure. Preflaw and post-flaw ultrasonic waveforms are analyzed from an aluminum plate specimen with artificially induced damage, and the AWI results show the location and spatial extent of all of the defects. INTRODUCTION Structural health monitoring (SHM) is of great interest for critical aerospace and civil infrastructure applications[1][]. The general goal is to identify, locate and characterize defects as they initiate and grow to a critical size during the service lifetime of a structure. Conventional nondestructive evaluation (NDE) methods such as ultrasonic testing, eddy current testing and X-rays can be used with great effectiveness for identifying and locating defects at the component level, but their use on structural assemblies is limited due to access restrictions and the high cost of field deployable instrumentation systems. The challenge is to develop a viable SHM technique that has the sensitivity of high resolution NDE methods but is cost effective for use on complicated structures. A methodology for such an approach is presented in the paper, and results are shown for an intentionally flawed aluminum plate specimen. In previous work, a sparse array of ultrasonic sensors was used to detect the onset of damage in structures [], and pre-flaw and post-flaw ultrasonic waveforms were compared using change detection algorithms to detect and classify damage [] []. Similar change detection algorithms are used in this paper to detect damage, and a time-of-flight algorithm is used for locating suspected damage regions. Awavefield visualization tool, referred to here as acoustic wavefield imaging (AWI), is used to snapshot ultrasonic waves as they propagate through a structure. The AWI method produces images

that are similar to those produced from numerical calculations using the Local Interaction Simulation Approach (LISA) [], and also to those produced using a Scanning Laser Doppler Vibrometer (SLDV) system []. METHODOLOGY A methodology for structural health monitoring and subsequent followup inspection is illustrated in Figure 1. A sparse array of permanently mounted piezoelectric transducers is used for in-situ ultrasonic measurements to determine if damage has occurred. Followup inspections are performed using an external ultrasonic transducer to scan suspected defect areas. A dual transducer method is used whereby the in-situ transducers are pulsed and the external transducer is the receiver. This methodology is an alternative to conventional NDE, and can most likely be implemented without disassembly of the structure for many applications. Ordinarily, followup inspections would be made using conventional pulse-echo and through-transmission ultrasonic inspection methods. Pulseecho methods are difficult to use on complicated structures, and the more frequently used throughtransmission approach usually requires disassembly of the structure to gain access to both sides. Furthermore, through-transmission ultrasonic inspections require precise position synchronization of transducers on each side of the part. However, the two transducer method proposed here using the in-situ transducers as fixed sources and a scanning external receiver simplifies the scanning hardware requirements, and results obtained are similar in sensitivity to those obtainable via throughtransmission ultrasonic inspections. A key aspect of this methodology is the use of change detection algorithms to compare both insitu discrete transducer measurements and followup wavefield scans to the respective data sets from previous inspections. This change detection approach is discussed in papers previously presented by Monitor In-Situ Transducers Damage Change Detection Analysis No Change Approved for Continued Service Usage Change Detected Data From Previous Inspections and Analyses Additional NDE Required To Characterize And Localize Damage Acoustic Wavefield Imaging (AWI) Pulse in-situ transducers Measure the acoustic wavefield using: - Scanned ultrasonic/acoustic transducer - Laser displacement vibrometer Conventional NDE Methods Ultrasonic testing (UT) Eddy current testing (ET) X-ray Flaw Characterization, Data Analysis, and Assessment of Structural Integrity Pulse on in-situ transducers, and record wavefields with an externally scanned receiver Cost effective and easily implemented Detailed/Tedious NDE Methods Larger area inspections required Expensive to implement, and may require disassembly of the structure Repair/Retire FIGURE 1. A methodology for health monitoring and followup inspection of critical structures.

the authors for analyzing data from the in-situ transducers [] [] [], and is demonstrated in this paper for the AWI results. EXPERIMENTAL DETAILS The aluminum plate specimen illustrated in Figure was used for this study, and damage conditions A E refer to various stages of damage progression in the plate caused by introducing artificial flaws. Both hole and notch type flaws were machined in the plate, and locations are shown in Figure along with the locations of four. MHz piezoelectric disc transducers that were permanently bonded to the back surface of the specimen. The following data were recorded for the undamaged baseline condition of the plate and damage conditions A E: (a) ultrasonic waveforms generated and received between the various permanently mounted transducers (in-situ measurements), and (b) ultrasonic waveforms received by scanning a conventional, flat. MHz transducer above the plate while pulsing the permanently mounted transducers (AWI measurements). For the in-situ measurements, waveform sets were recorded from the following transducer pairs (pulser to receiver): 1 to, 1 to, 1 to, to, to, and to. Because of reciprocity, pulsing on 1 and receiving on yields the same ultrasonic waveform as pulsing on and receiving on 1; thus, a total of waveforms were recorded for each of the aforementioned plate conditions. The top two waveforms in Figure are typical of those recorded from each transducer pair. All waveforms were captured using an bit, MHz digitizer and a data window of µseconds. (,1.) (.,1.) Zoomed Area For Waveform Energy Plots Notch. in x. in Transducer Transducer 1 (AWI Source) Hole 1 Diameter.1 in. in 1. in Transducer Notch 1.in x. in 1. in x. in Transducer. in Sparse Array Transducer Locations (. MHz) : Transducer Coordinates * 1 ( in,. in) (. in, in) (1. in,. in) (1. in, in) Artificial Defect Locations and Insertion Order: Defect Designation Size (in) Coordinates Hole 1 A.1 Dia. (1 in, in) Hole 1 B. Dia. (1 in, in) Notch 1 C. x. ( in, in) Notch 1 D 1. x. ( in, in) Notch E. x. (. in, in) FIGURE. defects. Aluminum plate specimen (.1 inch thickness) with permanently mounted transducers and artificial

The AWI measurements were made using an immersion scanning system whereby an ultrasonic transducer was scanned in a.1 inch square grid above the plate. The ultrasonic source was transducer 1 and the external transducer was the receiver. However, because of reciprocity, interchanging the source and receiver transducers does not significantly affect the recorded waveforms, and this was verified experimentally in our study. Thus, for an eventual field application of the AWI method, it may be more appropriate to use the external transducer as the ultrasonic source and record ultrasonic waveform responses on the permanently mounted transducers. Ultrasonic waveforms were recorded at each pixel location. Again, waveforms were captured using an bit, MHz digitizer and a data window of µseconds. Although not shown here, the individual waveforms are very similar in frequency content and detail to those from the in-situ transducer measurements shown in Figure. IN-SITU MEASUREMENT RESULTS Data from the in-situ transducers were used to determine the presence of damage and to estimate the location of the affected region. Specific details are beyond the scope of what can be described in this paper, but the steps are summarized here. Waveforms shown in Figure are typical of those recorded from the in-situ transducers. The basic premise is that these waveforms change when damage occurs in the structure []; this is demonstrated in Figure by subtracting the waveforms that were recorded before and after notch 1 was enlarged in length from. inches to 1. inches. Next, the differenced signal was squared and integrated to produce an accumulated energy waveform, and this waveform was analyzed to determine residual waveform arrival times. Two distinct arrival times may be noted on the waveforms shown in Figure : (a) arrival time of the leading wave packet is approximately µseconds, and (b) residual signal arrival time is approximately µseconds. The first arrival time is used to compute the acoustic wave velocity for the leading wave packet, and this velocity and time of the residual arrival are used to compute the total path length from the source to the suspected flaw and on to the receiver. The loci of possible defect positions on the plate which could have produced this path trace out an ellipse on the surface of the part with the positions of the two transducers as the focal points. 1 Waveform 1: Pulse on Transducer 1, Receiver on Transducer for. in long Notch 1 Amplitude 1 1 Waveform : Pulse on Transducer 1, Receiver on Transducer for 1. in long Notch 1 Amplitude 1 1 Residual Signal (Waveform Waveform 1) Amplitude 1 Accumulated Energy (Computed from Residual Signal) Energy 1 1 1 Time (microseconds) FIGURE. Signal processing steps used for in-situ measurements. 1

FIGURE. Arrival time ellipses projected on the specimen surface. The aforementioned algorithm was applied to each of the six data sets recorded from the various transducers pairs to obtain the total path length, and the resulting ellipses are shown in Figure. Note that the actual flaw location is very near the spatial position where most of the ellipses intersect each other. These computed ellipses were used to estimate the location of the flaw by summing the minimum distance from each point on the plate to each of the ellipses. The resultant area image plot of the summed distances is shown in Figure. The dark region is where the sum of distances is the lowest, and this region agrees very well with the actual flaw location. ACOUSTIC WAVEFIELD IMAGING RESULTS Snapshots of experimentally recorded acoustic wavefields are shown in Figure for material condition E; i.e., after all the artificial flaws were placed in the specimen. Results from four transit times are displayed to show the propagation of the wavefield out from the source transducer. Boundary Y (inches) 1 X (inches) FIGURE. Results from sum-of-minimum-distances algorithm.

reflections are clearly evident, and interactions of the leading wave packet with mounted transducers and damage regions are also visible. The following observations can be made regarding the wavefield images of Figure : (a) at µseconds the wavefront is slightly past notch (1: position), and just arriving at notch 1 (: position), (b) at µseconds the leading wave packet just passed transducer (: position) and is just arriving at hole 1 (: position), (c) at µseconds the wave is interacting with transducer (: position), and (d) at µseconds the wave is interacting with transducer (: position). An accumulated energy algorithm was used to enhance regions on the wavefield image plots where there were strong interactions between the propagating waves and discontinuities in the specimen. Specifically, voltage waveforms recorded at each scan pixel location were processed as follows: (a) energy was computed by squaring the waveform voltage, and (b) these computed energy values were accumulated with propagation time. Results are shown in Figure for the specimen after all the artificial flaws were introduced, and all discontinuities located within the scanned area are clearly visible. After µseconds the leading edge of the energy wavefield is past both notches 1 and and is just arriving at hole 1. After µseconds the leading edge is past transducers and. Next, residual accumulated energy images were formed by subtracting the accumulated energy waveforms before and after insertion of the various defects, and results are shown in Figure. This analysis method removes interactions from all non-changing discontinuities; i.e., only the interactions from the changed defect state are visible. Specifically, Figure shows only the interaction with hole 1 as interactions from transducers and were removed by this differencing algorithm. Figure (b) very clearly shows the response from notch 1 with interactions from hole 1 removed. Figure (c) shows the residual image for enlarging notch 1 from. inches to 1. inches. Note that the center portion of notch 1 image is removed and only interactions from the ends of the lengthened notch are visible. Figure (d) shows the final specimen state containing all the flaws differenced from the specimen baseline with no flaws. All three defects are clearly visible and notches are imaged at their proper angle of geometric orientation on the specimen. 1 (a) µseconds propagation time. (b) 1 µseconds propagation time. 1 (c) µseconds propagation time. (d) 1 µseconds propagation time. FIGURE. Snapshots of the ultrasonic wavefield radiating from transducer 1 after introduction of all artificial flaws.

1 1 1 (a) µseconds propagation time. (b) 1 1 1 µseconds propagation time FIGURE. Accumulated energy plots for the region shown in the zoomed area of Figure. 1 1 1 (a) Condition A - No Flaw Condition (b) 1 1 1 Condition C - Condition B 1 1 1 (c) Condition D - Condition C (d) 1 1 1 Condition E - No Flaw Condition FIGURE. Differenced accumulated energy wavefields for a propagation time of µseconds.

SUMMARY 1. A methodology has been developed for health monitoring and followup inspection of critical structures. It is a two step approach which utilizes a sparse array of permanently mounted ultrasonic transducers for in-situ monitoring, and a rapidly scanned external transducer for detailed analysis of regions with suspected defects.. A sum-of-minimum-distances algorithm has been used to analyze data from the in-situ transducers, and computational results accurately locate the area on the specimen with suspected damage.. Acoustic wavefield imaging (AWI) is demonstrated to be an effective tool for visualizing interactions of propagating ultrasonic waves with structural discontinuities.. All artificial defects introduced in an aluminum plate were successfully imaged by differencing pre-flaw and post-flaw accumulated energy wavefields. REFERENCES 1. S. R. Hall, The effective management and use of structural health data, in Proceedings of of the nd International Workshop on Structural Health Monitoring, pp., 1.. V. Giurgiutiu, A. Zagral, and J. J. Bao, Piezoelectric wafer embedded active sensors for aging aircraft structural health monitoring, Structural Health Monitoring, vol. 1, no. 1, pp. 1 1,.. T. E. Michaels and J. E. Michaels, Sparse ultrasonic transducer array for structural health monitoring, in Review of Progress in Quantitative Nondestructive Evaluation (D. O. Thompson and D. E. Chimenti, eds.), vol. B, (New York), pp. 1 1, American Institute of Physics,.. J. E. Michaels and T. E. Michaels, Ultrasonic signal processing for structural health monitoring, in Review of Progress in Quantitative Nondestructive Evaluation (D. O. Thompson and D. E. Chimenti, eds.), vol. B, (New York), pp. 1 1, American Institute of Physics,.. J. E. Michaels, A. C. Cobb, and T. E. Michaels, A comparison of feature-based classifiers for ultrasonic structural health monitoring, in Proceedings of SPIE Conference on Health Monitoring and Smart Nondestructive Evaluation of Structural and Biological Systems III (T. Kundu, ed.), vol., pp.,.. P. P. Delsanto and M. Scalerandi, A spring model for the simulation of the propagation of ultrasonic pulses through imperfect contact surfaces, Journal of the Acoustical Society of America, vol., no., p... M. Ruzzene, S. M. Jeong, T. E. Michaels, J. E. Michaels, and B. Mi, Simulation of measurement of ultrasonic waves in elastic plates using laser vibrometry, in Review of Progress in Quantitative Nondestructive Evaluation (D. O. Thompson and D. E. Chimenti, eds.), vol., (New York), American Institute of Physics,.