The Detection of Impact Damage to the Edges of CFRP Plates Using Extensional Ultrasonic Edge Waves

Extensional edge waves propagate along the edges of plates, with low attenuation in the propagation direction and amplitude decreasing rapidly (within one or two wavelengths) in the direction perpendicular to the plate edge. This makes them an ideal candidate for inspecting the edges of plate-like structures. Here, finite-element models and experiments are used to investigate the propagation and scattering of extensional edge waves in composite plates and application to damage detection is demonstrated. Piezoceramic transducers attached to the edge of a 4-mm-thick carbon-fibre-reinforced polymer (CFRP) plate were used to excite 140-kHz edge waves and damage detection demonstrated using two experimental configurations: pitch-catch and two-transducer pulse-echo. Damage due to edge-on impacts of 5 J and 10 J were detected in both configurations. A mass-addition damage model was used to investigate the effect of damage location. Detection was specific to damage within 33 mm (1.5 wavelengths) of the plate edge with propagation unaffected by features beyond that distance. The time of arrival of reflected signals in pulse-echo mode was accurately predicted using the edge-wave group velocity indicating that this configuration can be used for locating damage on edges. The localisation of edge waves means that they can specifically detect damage at the edges of structures, and attenuation measurements indicate that their useful propagation distance is large (5.5 m).


Introduction
Carbon-fibre-reinforced polymer (CFRP) laminates are increasingly being used in the construction of high performance vehicle panels, particularly fixed-wing aircraft primary structures. However, their potential advantages in terms of weight savings and high yield strength have yet to be fully exploited due to the need for conservative design due to the design risks associated with their susceptibility to impact damage from unanticipated discrete events, which can lead to severe subsurface delamination. In-plane impacts on the edge of the laminate are particularly dangerous as they generate more damage than out-of-plane impact for comparable energies [1,2]. One consequence of this is that the potential for lowvelocity impacts to the edges of wing-panel stiffeners due to dropped tools during manufacture or maintenance and the B Charles R. P. Courtney c.r.p.courtney@bath.ac.uk 1 Department of Mechanical Engineering, University of Bath, Bath BA2 7AY, UK difficulty in inspecting these structures leads to conservative design [2]. Improved non-destructive testing (NDT) and structural health monitoring (SHM) methods for detecting damage at the edge of panels, or other thin structures, is desirable in order to allow weight savings and reduced inspection costs. Ultrasonic edge waves are a guided wave mode that propagates along the edges of thin structures. This behaviour makes them an excellent candidate for inspecting the edges of plates, but they have not been applied to anisotropic layered structures, such as composites. In this work the detection of damage on the edges of CFRP panels using ultrasonic edge waves is modelled numerically and demonstrated experimentally.
The propagation of ultrasonic guided waves (particularly Lamb waves) in plates over large distances and subsequent scattering by damage has been widely investigated [3,4], but has seen limited application near plate edges, where reflection from the edge complicates the interaction with damage. Edge waves are guided waves that propagate parallel to the edge of a structure, with amplitude that decays perpendicular to the edge and so have the potential to address this limita-tion in guided wave inspection. At low frequencies (such that the longitudinal wavelength exceeds the thickness) two edge wave modes are present, a symmetrical extensional mode and an antisymmetric flexural mode [5]. Their low-dispersion and low-attenuation behaviour makes both modes attractive for NDT applications. The localisation of the waves to the plate edge allows transmission over substantial distances and simplifies inspection by focusing transmission along the edge, without attempting to inspect the more geometrically complex parts of the inspected structure. For example, provided that the height of a bonded stringer exceeds the penetration depth of the applied edge waves, the geometry, bonding and any filler at the base of the stringer should not affect the transmission, simplifying the inspection of the edges of such components.
Flexural edge waves have been applied to the inspection of structures that taper to a zero-thickness edge. These waves, often referred to as 'wedge waves' were first discovered in the finite element modelling of wedge shaped wave guides in the 1970s [6]. In 2010 an experimental investigation of the interaction of wedge waves with a crack demonstrated that they had potential in NDT and applications for tool edge inspection were proposed [7]. Recent work has experimentally and numerically demonstrated their suitability for inspecting long features with various tapered profiles, with particular application to turbine sealing fins [8].
For flat ended profiles, such as the edges of plates, symmetric extensional edge waves have attracted more attention than flexural waves. These waves are also referred to as 'Rayleigh edge waves', due to their similarity to Rayleigh waves on surfaces. Originally observed in models of rectangular waveguides [9], damage detection using extensional edge waves was demonstrated using laser generated ultrasound in a pitch-catch configuration to detect pitting in steel rail bases in 2012 [10]. Vien et al. numerically modelled the interaction of extensional edge waves in isotropic plates with cracks [11] and developed and experimentally demonstrated the scattering of edge waves by cracks on race-course shaped holes in aluminium plates [12]. Small cracks (less than one fifth of the wavelength) were shown to produce appreciable scattering (into both edge waves and plate modes) with amplitudes varying linearly with crack length. Extensional modes are investigated in this work due to their ease of excitation by piezoelectric plates and high velocity and low dispersion compared to flexural edge waves.
Progress towards application of edge waves in composites has been limited to numerical calculations of their behaviour [13,14]. Here generation and detection of edge waves in CFRP composite material is demonstrated experimentally and their use for the detection of on-edge delaminations investigated. The experimental work is supported by finite element modelling directed to: ensuring that the edge waves are correctly identified, elucidating their char-acteristics (principally penetration depth) and clarifying the interaction with the delamination inferred from the experimental results.
Section 2 introduces the general behaviour and analytical solutions for edge waves in isotropic materials, Sect. 3 describes the experimental setup and Sect. 4 describes the finite element model of the system. Experimental and numerical results are presented in Sect. 5 and discussed in Sect. 6.

Background
Propagating edge waves can have symmetric and antisymmetric (referring to the symmetry of the material motion through the thickness of the plate) modes of various orders in a manner similar to Lamb waves in plates. At low frequencies only the zeroth-order modes are present and, as the ratio of thickness to wavelength tends to zero, they become the thin-plate extensional edge waves and flexural edge waves, respectively. In this work the frequency regime is such that only these lowest-order modes are considered and the excitation is tuned to excite only the extensional/symmetrical mode. This section outlines the theoretical behaviour of extensional edge waves in homogeneous isotropic materials, followed by the expected deviations from this behaviour for composite materials as background to the new numerical modelling and experimental results of this work. Figure 1 illustrates the extensional edge-wave mode shape in isotropic materials. Figure 1a defines the system of axes. Propagation is in the positive x-direction. Figure 1b shows the amplitude of displacement oscillation in each direction as the distance from the edge increases; examples of the displacement are shown inset. In each case the wave propagates in the positive x-direction (parallel to the plate edge) with elliptical motion in the x y-plane that decreases with increasing distance, y, from the edge.
In an isotropic material thin enough that the plane stress approximation applies the complex material displacement in the direction of propagation can be shown to be [5]: and the displacement perpendicular to the surface is: v(x, y) = e iξ e k pl x k pl e −γ k pl y + Ae −δk pl y .
In each case the time dependence e −iωt is omitted.
is the axial-plate-wave velocity for a material of Young's modulus E, density ρ and Poisson's ratio ν and k pl is the corresponding wavenumber. The parameters A, γ and (1) and (2). Solid line is displacement parallel to x-axis, dashed line is parallel to y-axis. Displacement paths inset. All distances scaled to the wavelength of the edge wave, λ e δ are given by: where k e is the edge-wave wavenumber and the relative wavenumber of the edge wave, ξ e = k e k pl , is the real solution to The amplitudes of the displacements in Eqs. (1) and (2) are shown in Fig. 1b. The extensional edge wave is nondispersive with a velocity of It is notable that ξ e depends only on ν. A straightforward numerical solution of Eq. (4) for Poisson's ratio in the range 0 ≤ ν ≤ 0.5 indicates that ξ e varies monotonically from 1.61 to 2.15 (ξ e = 1.86 for ν = 0.31). So for most practical materials this means that the velocity of the edge wave is considerably slower than the axial plate wave (and hence the low frequency S 0 Lamb waves that tend to it for thin plates) and faster than the A 0 Lamb waves that tend to zero velocity at low frequencies.
Equations (1) and (2) indicate that the material displacement forms an ellipse with axes aligned in the x and y-directions. The amplitude and aspect ratio vary with depth from the surface according to the parenthetical terms in Eqs. (1) and (2). For motion perpendicular to the plate edge, v(x, y), the amplitude of motion always has the same sign. The behaviour of the motion parallel to the plate edge has a zero crossing near the surface which leads to a change in the direction of material rotation at . This ranges between 0.17 λ e and 0.26 λ e in the region 0 ≤ ν ≤ 0.5. In Fig. 1b this occurs at y = 0.20λ e , where the zero in the amplitude of the displacement in the x-direction can be seen.
Of particular interest for inspection methods is the depth of penetration of the edge wave, i.e. the distance from the edge at which the wave amplitude becomes negligible. From Fig. 1b v(x, y) is always greater than u(x, y) so we can take that as the amplitude of interest. For a material with ν = 0.31, v(x, y) drops to 10% of its peak value at y = 1.26λ e , which is marked on Fig. 1b.
Analytical investigations of edge waves in anisotropic materials are more limited, but their behaviour has been investigated numerically. Ohyoshi considered extensional edge waves in unidirectional fibre-reinforced composites using numerical methods [13]. The general behaviour (even in highly anisotropic materials) was broadly similar to that in isotropic materials, with waves propagating along the edge of the plate with elliptical displacement that decayed with distance from the edge. The main differences observed were that varying the principle axis of the material relative to the edge caused a variation in the edge wave velocity (but it remained non-dispersive) and caused the major axis of the elliptical displacement to become angled to the perpendicular. Cerv et al. used numerical modelling validated with experimental edge wave velocity measurements to investigate the behaviour of edge waves in unidirectional and cross-ply composites [14]. It was confirmed that the waves were non-dispersive, but that velocity depended on the orientation of the fibres relative to the edge.
In summary, the properties that are of particular interest for the development of a structural health monitoring technique are that extensional edge waves are non-dispersive, propagate along the edge of the plate with the largest displacements in the direction perpendicular to the edge and with amplitudes that decay over a distance of the order of a wavelength in the direction perpendicular to the edge. The non-dispersive behaviour is important as it allows long distance transmission of pulses without distortion, the displacement direction informs how best to generate and detect edge waves and influences how they will interact with damage and the localisation at the plate edge defines the region that can be inspected.

Experimental Method
Experiments were performed to determine: the group velocity and attenuation of the edge waves in undamaged CFRP plates, the effect of edge-on impact damage and side-on impact damage on the edge waves, and the effect of damage location on the detected edge wave signals. Methods of using piezoceramic transducers are well known for Lamb waves and have previously been applied successfully to aircraftrelated composite plates [15][16][17][18][19], here edge-mounted piezoceramic transducers were designed to preferentially generate and receive edge waves by tuning their dimensions and frequencies.
Two approaches to damage detection were investigated: pitch-catch and pulse echo. Figure 2 shows the two configurations. In each case the piezoceramic plates were permanently attached to the edge of a CFRP plate to generate and detect edge waves, the damage detection methods differed only in the positioning of the transducers.
In the pitch catch configuration two widely-separated transducers are used, with one transmitting a pulse that propagates along the edge until it is detected by the second transducer. Any damage occurring between the transducers causes scattering of the edge wave and reduces the amplitude of the detected wave.
In the two-transducer pulse-echo mode two transducers are positioned close together with the first sending a pulse. The second transducer detects this pulse almost immediately, followed by later pulses due to reflected waves. In this configuration scattering due to edge damage leads to detectable reflected edge waves. Baseline subtraction is used to identify new features in the received signal. By measuring reflected signals at transducers and calculating the time of flight (TOF), the damage location can be determined. This approach assumes that each echo in the signal was distinct from all others and isolated from interference and coherent noise [20].
In order to demonstrate damage detection, tests were performed using samples impacted on the edge with a drop impactor using energies of 5 J, 10 J and 24 J. To allow investigation of the effect of damage location, two physical models of damage were used: a mass added to the edge of the plate (to simulate on-edge damage) and a clamp across the plate (to simulate side-on damage). Both models were calibrated to the response of real damage created by impact.
Impact tests were conducted on 300 × 94 × 4mm samples at three different incident energy levels: 5 J, 10 J and 24 J . The 5 J and 10 J on-edge impacts were applied in the plane of the plate at 70 mm from the transducers, the damage location, sample dimensions and transducer locations are shown in Fig. 3a, b. On the third sample a 24 J near-edge impact was applied perpendicular to the surface of a sample at a distance of approximately 40mm away from the centre point of the plate edge, as shown in Fig. 3c. The edge-on impact energies were selected to give energy/thickness levels (1.25 J/mm and 2.5 J/mm) that are associated of BVID in out-of-plane impacts and comparable with energies for which the effect  [21]) and found to be significant. The out-of-plane energy was selected to induce significant damage to demonstrate the technique's insensitivity to defects away from the edge.
After the impact tests, the damaged area was inspected using a CT scan and an ultrasonic phased-array scan to determine its extent. The phased array (Diagnostic Sonar: 128 0.5-mm-width elements, 5 MHz-centre-frequency) was mechanically scanned along the surface of the plate (with coupling gel), while the elements were scanned in the direction perpendicular to the motion, to provide a C scan of the damaged region.
Each test piece was instrumented using soft-doped PZT plates (NCE51 from Noliac, Kvistgaard, Denmark) of thickness 2 mm and length 7 mm and width 4 mm as shown in  A similar configuration, with a 280 mm transducer separation, was used to detect the three different impact damages (on 300 × 94 mm samples) and is shown in Fig. 3.

Pulse-Echo Configuration
For experimental simplicity, a two-transducer method was used, where two piezoceramic plates placed close together were used to make the pulse-echo measurements. This allowed the same instrumentation to be used for both pulseecho and pitch-catch measurements with only the locations of the transducers altered. Figure 3b shows the pulse-echo setup used to detect the response due to the 5 J and 10 J edge-on impacts. A pair of transducers acting as the transmitter and receiver are located at the centre of the plate edge (on the 300 × 94 mm samples) with a 5-mm gap between them and pulse-echo measurements taken before and after the impact damage was introduced.
A similar arrangement (see Fig. 2b) was applied to detecting the location of simulated damage (a 206 g weight) at each of three equally spaced distances (62.5 mm, 125 mm and 187.5 mm) on the pate edge (of the larger, 600 × 300 × 4 mm plate).

Attenuation and Signal-to-Noise Ratio (SNR)
A scanning laser vibrometer (Polytec, PSV-400) was used to measure the amplitude of the propagating wave at increasing distances from the excitation transducer in order to determine the attenuation. A thin layer of a metallic silver water-based paint was applied to the plate edge in order to improve the consistency of the signal response. Assuming linear attenuation, the amplitude of the decaying wave is: where A (z) is the amplitude of the wave at distance z, and α is the attenuation coefficient [22]. An exponential curve was fitted, using a non-linear least squares method, to the amplitude as a function of distance to give the attenuation coefficient, α.

Signal Processing
The time domain signal was captured at a sample rate of 25 MHz, averaged over 100 repetitions and a digital lowpass tenth-order Butterworth filter with a cut-off frequency of 250 kHz applied to reduce unwanted noise in the pulseecho and pitch-catch measurements. The signal envelope was calculated using the Hilbert Transform and used to extract the peak amplitude and pulse arrival time.
For pulse-echo measurements baseline subtraction was used, whereby the residual signal was obtained by subtracting the baseline signal from the received signal. The baseline was recorded when the structure was in a healthy condition, so that the residual signal only includes information about damage.

Finite Element Modelling Method
The propagation and scattering of edge waves by delaminations was investigated using finite element (FE) models developed in Abaqus/CAE [23]. The numerical analysis was carried out by the Abaqus/Explicit solver.
C3D8R elements (8-node linear brick element) were used to discretise each layer of the composite layup. 16   laminae thickness (0.25 mm), and in the other two dimensions a mesh size of 1 mm by 1 mm was selected to be sufficiently small to resolve the smallest wavelength (11 mm for the highest frequency in the dispersion curves, 16 mm for the damaged-plate simulations) and capture the desired highest frequency response. A time increment of 1¯s was selected using the Courant-Friedrichs-Levy (CFL) condition [23]. The total number of elements was 450,000 and run times were of the order of 1 hour on a standard desktop computer.
The sensors were represented by single nodes across the plate thickness since they have relatively small size with respect to the distance from the transmitter. The excitation signals generated by the PZT transducer were simulated as displacements perpendicular to the plate edge. The imposed displacement is a 6-cycle sinusoidal Hann-windowed tone burst at 140 kHz, which is the same waveform as the electrical signal considered in the experiments (as further discussed in Sect. 5.1). The schematic illustrations for the setup of FE simulations are shown in Fig. 5.
To model multiple delaminations in composite laminates, the plate model was partitioned into two pairs of sub-plates, where a delaminated region (area: 30 × 30 mm 2 ) was defined on the inner contact surfaces of sub-plates (the location is shown in Fig. 5). The inner surfaces at the delamination region were excluded from the connection using a surfaceto-surface tie constraint, in conjunction with the master-slave relationship in Abaqus [23]. This produced discountinuties between the layers of the composite laminated plate. The delaminations were arranged parallel to each other and sym-

Results
The following results start with comparison of the predicted group-velocity dispersion curves with those obtained experimentally to ensure that the correct mode is detected (Sect. 5.1). This is followed by numerical prediction of the shape of the edge wave and its penetration depth in the CFRP sample (Sect. 5.2). Section 5.3 experimentally evaluates the attenuation of the edge wave. Having determined how the edge wave propagates in the sample the ability to detect damage is demonstrated: in Sect. 5.4 established nondestructive methods (ultrasonic C scan and x-ray CT imaging) are used to characterise the induced damage before detection with edge waves is demonstrated in pitch-catch (Sect. 5.5) and pulseecho configurations (Sect. 5.6). Figure 6 shows the group velocity dispersion curves for edge waves and Lamb waves. The experimental results were obtained from the time-of-flight of the first arrival, using the pitch-catch configuration described in Sect. 3 for edge waves with 15 mm diameter PZT disk transducers bonded to the plate face to detect Lamb waves. The central frequency of a 6cycle Hann-windowed tone burst signal was increased from 80 kHz to 300 kHz in steps of 10 kHz. A smaller number of similar experiments were performed numerically using the finite-element model implemented in Abaqus and described in Sect. 4. During the group velocity measurements the edge wave at 140 kHz had the highest amplitude and so was selected for the damage detection application; a 6-cycle Hann-windowed 140-kHz was used in the following sections. This frequency corresponds to a group velocity of 3050 ms −1 and a wavelength of 22 mm giving ratio of transducer length in the propagation to wavelength of 0.18. Figure 7 shows the displacement amplitude, obtained from the numerical models described in Sect. 4, 90 μs after application of the source displacement at the centre of the plate edge. The localised edge waves of interest are apparent on the upper edge. The use of displacements applied perpendicular to the edge as the source also leads to extensional plate waves travelling into the plate and shear horizontal waves at approximately 45 • to the plate edge. The wave field is not symmetrical about the centreline of the plate, in particular the edge waves are both angled in the same direction to the edge. This is due to the absence of reflection symmetry in the lay-up. Maio et al. [24] found similar dependencies of velocity on laminate stacking sequence for Lamb waves in composite laminated plates.

Edge-Wave Mode Shape
The presence of other wave modes as the edge wave decays makes exact measurement of the penetration depth of the edge wave difficult, but a dashed line at y = 55 mm = 2.5λ in Fig. 7 indicates the approximate penetration depth of the edge wave into the plate. Experimental Data Exponential fit Fig. 8 Laser vibrometer measurements on attenuation over distance along the plate edge Figure 8 shows the pulse amplitude as a function of distance with Eq. 6 fitted using a non-linear least squares method. The unattenuated amplitude A 0 with 95% confidence bounds is 0.040±0.001 and the attenuation coefficient

Attenuation
The measured noise of the processed signal received at 0.5 m from the source was 0.0004% of the signal amplitude (signal to noise ratio 108 dB). Assuming that the noise remains constant with distance and the signal attenuation continues at the above rate, useful information can be obtained using transducers with large separations. The noise will be 1% of the signal (signal to noise ratio: 40 dB) at 5.5 m transducer separation for pitch-catch measurements. Figure 9 shows the damage introduced to the 300 × 94 × 4 mm plates, after baseline measurements had been taken at 140 kHz. The on-edge impacts (at 5-J and 10-J) introduced visible dents on the impacted surface as seen in Fig. 9a, b. The near-edge CT image in Fig. 9c shows a 2D slice near (at a depth of 1.7 mm) the edge of the plate. Delaminations can be observed for both 5-J and 10-J impacts. For the 5-J impact, delaminations extend approximately 17 mm along the edge and for the 10-J impact, they extend 38 mm. These lengths correspond to 0.75 λ and 1.7 λ respectively at 140 kHz. The delaminations occurred at symmetrical positions about the mid-thickness of the specimen at depths of 1 mm and 3 mm from the top surface.

Characterisation of Damage
The 24 J side-on impact perpendicular to the plate surface produced damage that showed little change on the surface, but considerable internal disruption. Figure 9c shows an ultra-

Pitch-Catch Method
The effect of damage on the amplitude of transmitted edge waves in the pitch-catch configuration described in Sect. 3.2 was assessed by comparing the envelope of the received signal of the baseline (before damage) with that after the impact damage was introduced. In each case the signal was a 6-cycle

(b)
Baseline 5J Impact Baseline peak 5J Impact peak Fig. 10 Effect of on-edge impact damage on pitch-catch responses; Comparison of the normalised envelope for the composite specimen (size: 300×94×4 mm) obtained from a FE model with multiple delaminations; and from the experiments for the impacted plates at b 5 J, and c 10 J incident energy levels, located at approximately 70 mm from the transducers 140-kHz toneburst. Figure 10 compares the envelope of the first arrival of the pitch-catch signal responses in three scenarios; the FE simulation and experimental values for the two edge-impacted plates. In each case the signals are scaled to give a peak value of the undamaged signal equal to one. The FE model results shown in Fig. 10a indicate a 48% decrease in the out-of-plane displacement of the local maximum of the first arrival at 125 μs, and then a second local maximum at 158 μs reduced by 33%. Figure 10b, c display the amplitude of the envelopes from the signal responses for 5 J and 10 J on-edge impacts and their respective baselines. The changes in first-arrival peak amplitude are 4% and 41% respectively.
There is a small peak at 90 μs corresponding to the S 0 mode reflected at the back face and on to the receiver. Though small in this configuration, the S 0 is sufficiently large in the pulse-echo configuration to necessitate the use of baseline subtraction for those measurements. Another peak appears 30 μs after the direct edge wave. In the model result (Fig.  10a) this peak is reduced by the defect in a similar manner to the direct pulse, but it is largely unaffected by the impact damage in the experimental results. This time of flight corresponds to the edge wave initially travelling away from the receiver before reflecting from the plate edge and equivalently the pulse that passes the receiver before returning from the other edge. Beyond this second arrival the wave field becomes increasingly complex due to multiple reflections from the ends of the plates and mode conversions at reflection.
To investigate the effect of defect location, a damage proxy was required. After preliminary tests, a 206 g weight was selected as giving results most similar to those seen for the 10 Simulated damage 24J Impact damage Fig. 12 Effect of applying load and 24-J side-on near-edge impact damage on plate surface on pitch-catch signals; The percentage difference of the local maxima extracted from the envelopes with a G-clamp, applied to plate surface with varying distances from the edge J impact damage. Figure 11a presents two different envelopes from the signal responses, with and without the 206g weight. A decrease of 33.5% in the peak amplitude was observed relative to the baseline. Repeating this measurement (replacing the same mass at the same point each time) resulted in a distribution of peak amplitudes with a standard deviation of 0.6%. As in the modelled results the second edge-wave peak is reduced by the presence of the mass. The effect of varying the weight location is shown in Fig. 11b: the percentage change in peak amplitude due to the weight at different locations varied between 31 and 37%. Similarly a proxy was needed to investigate sensitivity to damage with increasing distance from the plate edge. Figure  12 shows the effect on the pitch-catch signal of a clamp across the plate at increasing distances from the edge. Where the clamp is located at the edge it reduces the transmitted signal by 12%. This effect decreases to 2% as the clamp approaches 35 mm from the edge, which is approximately one-and-half wavelengths (33 mm) from the edge. At distances greater than 45 mm (2λ) the percentage change settles at 0.6%. At 10 mm from the edge the resulting change in amplitude is similar to the 5% change seen in the the 24 J near-edge impact damage, which extends from 10 mm from the plate edge.

Pulse-Echo Experimental Results
In the pulse echo method of detecting damage, baseline subtraction was used to produce a residual signal that measures the waves scattered by damage in the sample back to the transducer location (source and receiver were located close to each other). Figure 13 demonstrates the detection of the impact damage using the pulse-echo response. After the damage has been introduced, a peak occurs near the expected arrival time for a defect 70 mm from the transducers (45 μs, marked with a vertical black line) in both Fig. 13a, b. For the 5 J impact the peak amplitude is at 41.5 μs, corresponding to a reflector at a distance 63 mm to the transducer. For the 10 J impact the maximum signal occurs at 45 μs with another peak at 31 μs, which corresponds to a distance from of 47 mm. Due to the broadly one-dimensional propagation, there are also changes in the signal at times after the arrival of the reflected signal, as all reflections from the end of the plate also pass through the damaged region. The simulation results in Fig. 13c also confirms that the local maximum aligned with the arrival time of the reflected pulse from the discontinuities. Figure 14 presents the residual signals left after baseline subtraction using the 206-g mass damage model. The vertical black solid lines represent the expected TOF from the mass location, calculated as the time taken to cover the known distance to the simulated damage and back at the group velocity of the edge wave (3050 m/s at 140 kHz). Figure 14a-c show that the arrival time of the wave packet varies with the distance of the weight as predicted by the group velocity. In each case a large residual signal occurs just after 150 μs (independent of defect location). This corresponds to the expected arrival time of the reflection of the plate edge behind the defect (163 μs). This reflected pulse passes through the damaged region twice leading to a substantial decrease in amplitude (of the order of that observed in pitch-catch measurements) which appears as a large pulse in the residual signal.

Discussion
Using transducers in a pitch-catch configuration, it was seen that edge-on impacts cause an appreciable reduction in amplitude in a signal received after passing through the damage region, indicating that plate edge damage can be detected by this method. For a 10 J impact this decrease amounted to 41% of the amplitude of the first arrival pulse in the undamaged case. This reflects the significant delaminations and matrix cracking observed in Fig. 9. For the 5 J case, with less extensive damage, the reduction in amplitude was 4%, which is less obvious in the plotted data, but exceeds the variation seen in repeated measurements of 0.6%. Using an edge mounted mass, selected to give a similar response to the 10 J damage it was shown that the change in the signal was approximately the same regardless of the position of the simulated damage between the transducers. A second-arrival edge wave signal is received shortly after the first due to reflections from the plate edge. In the model and simulated damage results this behaves as might be expected, with a reduction in amplitude when the simulated damage is added. For the plate with 10-J-impact damage the second arrival seems unaffected by the damage. This behaviour is unexpected and may relate to the symmetrical location of the transducers leading to a complex interaction of two or more paths.
Pulse-echo measurements demonstrated that edge waves successfully detected impact damage at both 5 J and 10 J. The arrival time of the first reflected pulse in each damage case (41.5 μs and 31 μs, respectively) can be converted to a distance using the group velocity and gives distances of 63 mm and 47 mm. These values agree well with the separation of the transducers and the nearest edge of the damage (61 mm and 51 mm respectively) obtained from the images in Fig. 9a-c. The prediction of the damage location was further explored using simulated damage (added mass) and the time of arrival of the reflected waves accurately can be predicted using the time of flight, at a range of positions.
In order to demonstrate the localisation of the waves at the plate edge and the specificity of the technique to damage at the edge, simulated damage consisting of a clamp tightened across the plate was used. 24 J impact damage centred on 40 mm from the edge resulted in a delamination extending to 10 mm from the edge was used to indicate that the effect of the clamp is comparable to real damage. As the clamp head was moved away from the edge, the change in the signal reduced to reach a level of approximately 0.5% at approximately 45 mm (≈ 2λ) and did not decrease further. This is of the order that might be expected from solutions to Eqs. (2) and (4), which analytically determine the oscillation of the edge wave perpendicular to the plate in isotropic materials and indicate that the amplitude drops to 1% of its peak (on-edge) value at 2.2λ. The experimental results are also compatible with the penetration depth obtained from the FE model of the quasiisotropic plate (see Fig. 7) which is approximately 2.5λ.
The finite-element model uses a line source resulting in significant generation of Lamb waves that travel into the plate, as seen in Fig. 7. The experimental method uses transducers of finite width (4 mm) in the edge-wave propagation direction and the frequency of operation (140 kHz) is selected to maximise the edge-wave amplitude. This provides sufficient specificity in generating edge waves (rather than Lamb waves) to demonstrate edge damage detection in the pitchcatch configuration. However, the generation of some Lamb Waves which are reflected at the opposing edge of the panel to the transducers, requires that baseline subtraction is used for the pulse-echo measurements. There is potential for baselinefree inspection if greater mode-specificity can be achieved. Wedges and inter-digital transducers are methods that merit investigation to achieve this.
The one-dimensional propagation of the waves along the edge allows large transmission distances, with only the attenuation of the waves reducing their amplitude, in contrast to the reduction of amplitude due to two dimensional spreading in guided waves. The attenuation coefficient of 1.3±0.1 suggests that transducer separations of up to 5.5 m are possible.

Conclusions
The experimental and numerical results presented provide demonstrate the viability of damage detection with extensional edge waves in composite components and provide an investigation of the propagation behaviours relevant to edgedamage detection.
The attenuation of the edge wave was measured experimentally and showed that edge waves are suitable for NDT and SHM applications over long distances (up to 5.5 m). Both pitch-catch and two-transducer pulse-echo configurations were shown to be able to detect 5 J and 10 J on-edge impact damage, with the pulse-echo response demonstrating the ability to locate the damage. Pitch-catch and pulse-echo detection methods were demonstrated over a range of damage locations using added mass as a damage model. The edgewave detection method was specific to damage near (within 33 mm) to the edge and detectability did not vary with location along the edge.
The localised inspection of edges using permanantly attached sensors allows regions to be monitored that are difficult to reach with Lamb waves and has the potential to allow targeted monitoring of important parts of structures while avoiding wave propagation in areas that are of less interest, but complicate the signals significantly. The use of extensional edge waves to inspect anisotropic materials significantly broadens the range of materials and geometries that can be inspected with guided edge waves.
is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/.