Material properties of PZT and CFRP lamina
1. Introduction
To improve the safety and reliability of various engineering structure, it is essential to develop efficient techniques for non-destructive damage detection or structural health monitoring. Lamb wave can travel a long distance in plate-like and shell-like structures made of materials even with high attenuation ratio (e.g. Carbon Fibre/Epoxy Polymer composites). To take this advantage, many researchers have recently explored the possibility of using Lame waves for damage identification [1]. To date, many developed Lamb wave-based techniques are generally based on so called two-stage prediction models by which the difference in the signals between a defective structure and a benchmark (intact structure) can be evaluated. Then, the residual error is easy to be defined no matter what information extracted from the signals is used, such as the information in time domain [2, 3] or frequency domain [4, 5]. Therefore, a benchmark or baseline signal is essential for the detection, which is very reliable and suitable for monitoring the propagation of damage. Also, tremendous efforts have been put to the delamination identification, which could be treated as a problem of inverse pattern recognition using calibrated numerical methods such as artificial neural network [6]. The interaction between Lamb waves and delaminationhas also been investigated numerically and theoretically [7-9]. However, the complex wave scattering phenomenon in a delamination area has not been clearly understood in these studies.
In this chapter,a technique for delamination identification in laminated composites using the zero-order mode of Lamb waves, i.e., S0 mode and A0 mode, without referring to the baseline data, is described. Through measuring the propagation speed of a wave and the traveling time of a reflected wave from thedelamination, the delamination position can be accurately identified. Moreover, to understand the complex interaction of Lamb waves with a long delamination damage, the numerical simulations have been carried out.
2. Detective technique without the baseline data
2.1. Experiments and numerical analyses
In general,comparedwith S0mode, the attenuation of A0 mode in structures is more severe.Here, S0mode corresponds to an axialdeformation mode while A0 mode corresponds to a flexuraldeformation mode.Therefore, it should be paid more attention to the detective capability of both modes for a delamination damage in the different interfaces along the thickness direction. Moreover, the experimental setup of S0mode isslightly different with that of A0 mode.
2.1.1. Materials and experimental procedure
For S0mode,as shown in Fig. 1, a CFRP laminated composite beam of stack sequence of [010/9012/010] was used. An artificial delaminationdamage with different lengths for several cases, i.e., 30 mm, 20 mm and 10 mm, respectively, was intentionally created at the interface between the 10th and 11th plies by inserting a Teflon filmwitha thickness of 25μm. A PZT actuatorwas attached on the top surface of the left end of the beam. The PZT actuator had a diameter of 10 mm and a thickness of 0.5 mm. In this case, the generated waves from the actuator merged into the reflected waves from the left end of beam, which leads to a simpler signal. And the same PZT unit was used as a sensor to pick up the reflected wave from the delamination damage. Because the propagation speed of S0mode is much higher than that of A0 mode, confirmed as 4 times higher in the used composites from our testing, it is expected that in a specified time domain, no A0 mode will be collected due to its slow speed, making the analysis of S0 modes much easier.
For the A0 mode, a similar composite beam was used, as shown in Fig. 2. Two kinds of stack sequence of the laminated beam were used, i.e., [010/9012/010] and [012/04/04/012]. Note that [012/04/04/012] is a unidirectional laminate, and it can be simply expressed as [032]. To simplify the description of delamination later, this expression is used. Like the case of S0 mode, a delamination damage with different lengths, i.e., 30 mm, 20 mm and 10 mm, was intentionally created at the interface between two plies by inserting a Teflon film with a thickness of 25 μm. The distance between the center of the delamination and the left end of beam is 790 mm. To generate the A0 mode, two actuators were attached on the top and bottom surfaces of the beam with applied out-of-phase voltages because the difference of signals in the two actuators can produce the pure
where
Without the baseline data, it is still easy to estimate the arrival of the reflections from the beam boundaries, provided that the wave propagation speed and the dimensionsof the beam are known. Then, if the reflected signal from the delamination is not overlapped with the reflections from the boundaries, it can be simply collected by the sensor. Therefore, it is still a technical challenge when the delamination is located close to the beam boundaries and the reflected signalofthe delaminationis completely overlapped with the reflections from the boundaries.
Material properties | |
PZT | E11=62 GPa, E33=49 GPa, d33=472 pC/N, d31=-210 pC/N, |
CFRP lamina | E11=115 GPa, E22=E33=9 GPa, G12=G23=5.5 GPa, G13=3 GPa, |
2.1.2. Finite element analysis procedure
To explore the wave propagation in laminates with a long delamination case, a three-dimensional8-noded brick hybrid element proposed by the authors [10] and the explicit time integration algorithm were used in the finite element simulations without considering the dynamic contact effects like those in some previous studies [11-13].For the delamination area, one node used on the intact interface along the through-thickness direction, and double nodes were used on the delamination interface, i.e., one belongs to the elements of upper delaminated portion, and the other belongs to the elements of lower delaminated portion. The contact effects in the delamination area were neglected. The same signal stated in Eq. (1) was also used here. Furthermore, the actuator and the sensor were discretized using 3D brick elements. If an electrical field is applied on the actuator, the tensional or compressive strains will be generated in it from a proper relation, which connects the applied voltage and the generated internal strains. The induced stresses from the generated strains can be used tocalculate the elemental nodal forces, which form axial force and bending moment in the beam simultaneously. Then, both S0and A0 modes can be automatically generated. The material properties of CFRP are shown in Table 1.
2.2. Results of S0 mode
2.2.1.Experimental. analysis
Firstly, to evaluate the wave propagation speed, an intact beam with two sensors attached,was employed. The distance between the two sensors was 420 mm. By using the wavelet transformation technique [14], the arrival times of the incident wavesto the two sensors were determined and the wave propagation speed of S0 mode at
A comparison of the signals from an intact beam and a delaminated one (a 30 mm delamination damage) for the case in Fig. 1 is shown in Fig. 3. It is clear that there is a reflected wave in the signals from the delaminated beam, which is located between the incident and the reflected waves. However, no clear reflected wave could be detected when the delamination length was reduced to 10 mm. No A0 mode can be observed in Fig. 3, since its incident wave is not fast enough to reach the sensor within 300μs. For a signal with
2.2.2. Numerical simulation
A comparison between the experimental result and numerical simulationfor the case of Fig. 1 is shown in Fig.5. A good agreement between two results is observed. For convenience, the amplitude of the first arrival S0 mode in the simulated waves was calibrated using the experimental data. The attenuation coefficients of the CFRP in the numerical model were determined by matching the amplitudes of the reflected waves from the right end of the beam to the experimental data. When the actual sensor thickness, i.e., 0.5 mm, was used in numerical simulations, a small reflected wave from the sensor could be observed immediately next to the incident wave. This was also confirmed by the experimental result. Corresponding to a 0.05 mm sensor thickness, however, no obvious reflection from the sensor could be observed in the numerical simulation. Therefore, the sensor thickness of 0.05 mm was eventually chosen in the simulations. Note that the selection of sensor thickness does not affect the simulation results as the calculated amplitudes will be calibrated using the experimental data. Similar to the experimental results, the numerical simulation showed that the reflected wave from the 10 mm delamination was very weak.
To understand the effect of delamination position along the through-thickness direction on the propagation of a Lamb wave, a 30 mm delamination damage was created in the mid-plane of the laminated beam, i.e., between the 16th and the 17th plies. With the signal of
To examine the capability of the numerical simulation to identify the location of delamination, the length of delamination Lwas increased from 30to 90mm. The length of the beam was also increased up to 1500 mm. Under a condition o
2.3. Results of A0 mode
To confirm the effectiveness of present excitation and sensing techniques for A0 mode (Fig. 2), numerical simulation was firstly conducted for the laminate with [010//906/906/010] and 30 mm delamination under a condition o
The results of the laminate [010//906/906/010] with 30 mm and 10 mm delamination damages are shown in Fig 9. Sampling time for both numerical and experimental resultswas 1.0×10-4 ms. For the delamination of the lengths of 30 mm and 10mm, the reflections from the delamination can be clearly identified. In Figs. 9(a) and 9(b), it is interesting to note that the reflection signal from the 10 mm delamination is stronger than that from the 30 mm delamination. The reason may be that the reflections from the two ends of delamination are more easily overlapped in a shorter delamination case, e.g., the 10 mm one, resulting in a higher total reflection. Therefore, the intensity of the reflection from the delamination does not certainly depend on the length of the delamination. When the delamination is located on the mid-plane of the laminates, i.e., [010/906//906/010], the results are shown in Fig. 10. Similar to Fig. 9, the reflection can be identified clearly. Also, the reflection from 10 mm delamination is still stronger than that of 30 mm. For the short delamination, e.g., smaller than 30 mm, reflections from the two ends of the delamination should be overlapped if we consider the duration time of a 5 cycle signal (
From Fig. 9 to Fig. 11, a good agreement between experimental and numericalresults is observed. For convenience, the amplitude of the first arrival A0 mode in the simulated waves was calibrated using the experimental data. The attenuation coefficients of the CFRP in the numerical model were determined by matching the amplitude of the reflected waves from the right end of the beam to that of experimental data. Naturally, the material properties of CFRP in Table 1 were also adjusted slightly within a reasonable range to match the numerical results to the experimental ones.
To determine the wave propagation speed of A0 mode, two sets of distant sensor pairs were attached on an intact laminated beam. The wave speed was determined from the distance and the difference of arrival times between the two sets of sensor pairs. For the laminate with [010/906/906/010], the experimental and numerical wave speeds were 1555 m/s and 1470 m/s, respectively. For the case of [012/04/04/012], the experimental and numerical wave speeds were 1723 m/s and 1616 m/s, respectively. From these results, it can be found that the wave speed of [012/04/04/012] is only slightly higher than that of [010/906/906/010]. The reason is that the wave speed of A0 mode is determined by the bending stiffness of laminates, which is dominated by the fibre orientation of outer layers near top and bottom surfaces of laminates. In both cases, 0° degree layers are located near the two surfaces of laminates, which results in the comparatively small difference of wave speeds in both cases. For the case of [010/906/906/010], if 90° degree layers were located near the two surfaces of laminates, it would be a completely different story.
With the knowledge of the wave speeds, the delaminationpositions for various situations were evaluated, as shown in Table 2. The distances between identified positions and the actual centers of delamination are listed. In this table, the negative values denote that the identified positions are located on the left side of the delamination centers, and positive values represent that the identified positions are on the right side of the delamination centers. In this table, it is clear that the delamination positions in almost all laminates have been successfully detected compared with the beam length. For the laminate with [010//906/906/010] and 20 mm delamination, the reflection is not very clear due to strong noises in experiments. As a result, only numerical estimation is presented in this table. The similar case can be found in the laminate with [012/04//04/012] and 30 mm delamination, as shown in this table.
In summary, unlike S0 mode, A0 mode can be used for different applications. However, the amplitude of reflection from delaminationmay not depend on the length of the delamination. The amplitude of reflection from short delaminationis higher than that from long delaminationdue to overlapping of reflections from the ends of short delamination. The disadvantage of A0 mode is the short travelling distance due to its high attenuation in CFRP materials. Therefore, if possible, both the actuators and sensors should be placed close to the potential damage sites.
Unit (mm) | [0/90//90/0] | [0//90/90/0] | [0/0//0/0] | [0//0/0/0] | |
30mm | Num. | -11.0 | 4.0 | 0.0 | 18.0 |
Exp. | -19.0 | 39.0 | -11.0 | ||
20mm | Num. | -12.0 | -5.0 | -14.0 | 18.0 |
Exp. | 9.0 | -11.0 | 4.0 | ||
10mm | Num. | 2.0 | 2.0 | 1.0 | 1.0 |
Exp. | 2.0 | -2.0 | 10.0 | -17.0 |
3. Interaction between Lamb waves and delamination
3.1. S0 mode
To gain a better understanding of this complex interaction, a beam with a length of 1500 mm and a 200 mm long delamination damage was examined using the numerical simulation under a condition of
At the time domain of 139μs, the S0 mode arrives at the left end of the delamination and then reflects from or transmits from the left end, as shown in Fig. 12(a). At 182μs, the S0 mode passes through the right end of the delamination. Fig. 12(b) shows two separate reflected modes, i.e., A0 and S0 from the left end of the delamination. In principle, it is easy to distinguish a S0 from an A0 mode by the propagation speed and wave length. However, the amplitude of the reflected S0 mode is smaller than that of the A0 mode. This can be explained by the fact that thebending deformation in a laminated beam is much higher as compared to the axial deformation. Also, when the S0 mode travels through the left end of delamination, a new transmitted A0 mode is generated as shown in Fig. 12(b). The reflected waves from the right end of the delamination, consisted both A0 and S0 modes, can be also observed in Fig. 12(b). Corresponding to 236μs, the incident A0 mode from the actuator does not reach the sensor, andthe reflected S0 mode from the left side of the delamination has already passed through the sensor. However, at the same time point in Fig. 13, this reflected S0 mode cannot be detected by the sensor due to its small amplitude. Also, the reflected S0 mode from the right end of the delamination, which is surpassing the reflected A0 mode from the left end of the delamination, can be observed, as shown in Fig. 12(c). A new transmitted A0 mode and the transmitted S0 mode from the right end of the delamination can be clearly observed, Fig. 12(c). At 303μs, the incident A0 mode does not arrive at the sensor. As indicated in Fig. 12(d), the reflected S0 mode from the right end of the delamination completely surpasses the reflected A0 mode from its left end. This reflected S0 mode has already passed through the sensor position. From Fig. 13, the arrival time of the reflected wave from the delamination is identified as around 275μs. The reflected wave actually arrives at the sensor betweenthe time points shown in Figs. 12(c) and12(d), i.e., 236μs and 303μs. Therefore, the reflected wave signaldetected clearly by the sensor is considered to be the reflected S0 mode from the right end rather than the left end of the beam, due to its higher intensity. Based onthe numerical simulations above, the very complex interaction between the different signal modes and the boundaries of delamination can be identified. When only a single S0 mode passes through the delamination, four modes will be generated at one end of the delamination, including two transmitted A0 and S0 modes, and two reflected A0 and S0 modes. Further study is required to understand why the stronger reflected A0 and S0 modes. However, the reflected S0 mode received by a sensor is from the reflections are from the right end of delamination. One possible explanation is that the delaminated region is of lower bending stiffness, which can be considered to be a softer region. When a wave propagates from a harder or intact region into a softer region, the reflections become weak. In contrast, when the wave propagates from a softer region into a harder one, a stronger reflection is expected.
3.2. A0 mode
Like S0 mode, abeam with a length of 1200 mm and a 200 mm long delamination case was examined using the numerical simulation under a condition of
Firstly, for the laminate with[010//906/906/010],at the time domain of 318μs, two A0 mode waves were excited by the actuator pair, which propagates to the left and right directions, independently. The A0 mode wave propagating to the right direction arrives at the left end of the delamination and then reflects from or transmits from the left end, as shown in Fig. 14(a), respectively. At 433μs, the transmitted A0 mode from the left end of delamination passes through the right end. Fig. 14(b) shows thereflected wave from the right end of delamination. By comparing Fig. 14(b) and Fig. 15 at the same time domains, we can identify that the reflected A0 mode from the left end of delamination passes through the sensor. The reflected S0 mode due to mode change can also be identified. In principle, it is easy to distinguish a S0 from an A0 mode by using phase information of signals from the two sensors. At 488μs, the reflected and transmitted waves at the right end of the delamination separate to each other, as shown in Figs. 14(c). Also, the amplitude of the reflected wave is much higher than that of the reflected wave at the left end of the delamination, confirmed in Fig. 14(b). This phenomenon is similar to that obtained for S0 mode. Fig. 14(d) shows the transmitted and reflected waves at the left end of the delamination when the reflected wave from the right end of the delamination arrives at the left end at 621μs. It is worthwhile mentioning that although the amplitude of reflected wave from the right end of delamination is very strong, the reflection at another end (left end) of the delamination can significantly reduce the amplitude of the transmitted wave, which is expected to be monitored by the sensors. Therefore, wave mode change and multiple reflections may create complicated interactions between theA0 mode and the delamination.
In Fig. 15, the amplitude of the reflected wave from the right end of the delamination is still higher than that of the reflected wave from the left end. Similar results can be observed in the laminate with[012//04/04/012]. In consequence, as shown in Table 2, more identified delamination positions are close to the right end of delamination for the laminates with [010//906/906/010] and [012//04/04/012]. This phenomenon is similar to the result obtained for S0 mode in the case of [010//906/906/010].
For the case of [010/906//906/010], the signals of two sensors are plotted in Fig. 16 with the enlarged picture for the reflection from the delamination. In Fig. 16, there is no S0 mode and the mode change does not occur as the delamination is located on the mid-plane of the laminates. In contrast to the cases of [010//906/906/010] and [012//04/04/012],it is interesting to note that the amplitude of the reflected wave from the left end of delamination is higher than that of the right end. Similar results are obtained in the case of [012/04//04/012]. Corresponding to Table 2, more identified delamination positions are close to the left end of the delamination. The physical meaning of why the stronger reflection occurs at the different ends of delamination for different delamination cases in Figs 15 and 16 is not clear. The position of the stronger reflection seems to mainly depend on the following two factors: bending stiffness of upper and lower delaminated layers and their difference, and Lamb wave deformation mode.
4. Conclusion
A technique for identifying a delamination damage in cross-ply laminated composite beams has been developed by using the zero-order mode of Lame waves. The delamination position can be identified if the arrival time of the arrival time of a reflected wave from a delamination and the wave propagation speed of Lame waves are known. One of the advantages of this approach is the baseline data acquired from intact beams are not required.
And, extensive finite element simulations have been conducted to investigate the interaction of Lame waves with various delamination cases. For S0 mode, when a single S0 mode propagates into a delamination, mode change happens and reflected S0 mode with high amplitude is generated from the end point of the delamination where the waves propagate out; for A0 mode, it has been found that the situation is more complex than S0 mode. For the laminates with stack sequence of [010//906/906/010] and [012//04/04/012,], mode change also happens and the reflected wave from the right end of delamination is slight higher than that from the left end. However, for the case of [010/906//06/010] and [012/04//04/012], there is no mode change when a single A0 mode passes through delamination, and the stronger reflection from the delamination occurs at the left end of the delamination.
References
- 1.
Guided Lamb waves for identification of damage in composite structures: A review. J Sound & VibrSu Z. Q. Ye L. Lu Y. 2006 295 753 80 - 2.
Proceedings of DETC, Pittsburgh, PA, USA;Valdes S. H. D. Soutis C. A. structural health. monitoring system. for laminated. composite 2001 - 3.
Damage diagnosis using time series analysis of vibration signals. Smart Mater and StructSohn H. Farrar C. R. 2001 10 1 6 - 4.
Ihn JB, Chang FK. Built-in diagnostics for monitoring crack growth in aircraft structures. In: Chang Fu-Kuo (ed.), Proceedings of the international workshop on structural health monitoring. Stanford University, USA;2001 - 5.
Automated methodology to locate notches with Lamb waves. Acoustics Research Letters OnlineHurlesbaus S. Niethammer M. Jacobs L. J. Valle C. 2001 2 97 102 - 6.
Assessment of delamination in composite beams using shear horizontal (SH) wave mode. Compos Sci TechnolSu Z. Yang C. Pan N. Ye L. Zhou L. M. 2007 67 244 51 - 7.
94 2240 6 Guo N, Cawley P. The interaction of Lamb waves with delaminations in composite laminates. J Acoust Soc Am 1993; 94: 2240-6 - 8.
Wang CH, Rose LRF. Wave reflection and transmission in beams containing delamination and inhomogeneity. J Sound & Vibration2003 264 851 72 - 9.
Compos StructMahapatra D. R. Gopalakrishnan S. A. spectral finite. element model. for analysis. of axial-flexural-shear. coupled wave. propagation in. laminated composite. beams 2003 59 67 88 - 10.
Int J Numer Meth EngrgCao Y. P. Hu N. Lu J. Fukunaga H. Yao Z. H. A. . D. brick element. based on. Hu-Washizu variational. principle for. mesh distorsion. 2002 53 2529 - 11.
Impact analysis of composite laminates with multiple delaminations. Int J Impact EngHu N. Sekine H. Fukunaga H. Yao Z. H. 1999 22 633 48 - 12.
Comput & StructHu N. A. solution method. for dynamics. contact problems. 1997 63 1053 63 - 13.
Low-velocity impact-induced damage of continuous fiber-reinforced composite laminates. Part II. verification and numerical investigation, Compos part ALi C. F. Hu N. Cheng J. G. Sekine H. Fukunaga H. 2002 33 1063 72 - 14.
Wavelet analysis of plate wave propagation in composite laminates. Compos StructJeong H. Jang Y. S. 2000 49 443 50 - 15.
Nayfeh AH, Chimenti DE. Elastic wave propagation in fluid-loaded multi-axial anisotropic media. J Acoust Soc Am1991 89 542 9 - 16.
Quantitative damage detection in cross-ply laminates using Lamb wave method.Compos Sci and TechnolToyama N. Noda J. Okabe T. 2003 63 1473 79 - 17.
Development of smart composite structures with small-diameter fiber Bragg grating sensors for damage detection: Quantitative evaluation of delamination length in CFRP laminates using Lamb wave sensing. Compos Sci TechnolTakeda N. Okabe Y. Kuwahara J. Kojima S. Ogisu T. 2005 65 2575 87