# Acoustic emission for on-line monitoring of damage in various application fields.

AbstractStructural damage in materials is a "top priority" to be investigated by material scientists all over the world. The future or limit of a material to be used in structural applications depends upon it, as safety and environmental protection are important issues in national policies. The acoustic emission (AE) technique offers the possibility to listen to the damage process when and where it occurs. This paper focuses on a number of cases investigated at the department MTM using AE to analyze damage in various application fields. The tools handled to interpret the AE signals, the relation between AE and the residual properties and the additional understanding of a particular material behavior will be outlined. The fatigue damage of CFRP, the curing of concrete and the corrosion of carbon steel and stainless steel will be discussed as case studies.

Introduction

The acoustic emission technique is a "passive" method, which detects the elastic energy released when microstructural changes in a material or a process take place. The AE sensors, mostly PZT transducers, resonant or broadband, are used to pick up the elastic stress waves from the material and convert them into an electrical signal. This electrical signal is further processed by filters, amplifiers and a signal processing unit on a PC. Depending on the application the correct frequency bandwidth has to be chosen to ensure the best sensitivity of the detection system for the phenomena to be measured. The sampling rate also depends upon the latter. The further tuning of the signal-to-noise ratio (filters/amplifiers) can be done to optimize the monitoring conditions. For the analysis of the AE data it can be of help that with the AE signal also a recording is done of the time of loading or the loading itself (mechanical or thermal). From the AE signal itself the following parameters or AE characteristics are recorded: amplitude, duration, AE (ringdown) counts, rise time, energy, slope, frequency, modal content of the AE signal, AR (auto regressive) and wavelet parameters.

Case study 1: Fatigue damage monitoring of carbon fiber-reinforced epoxy laminates with AE

Tension-tension fatigue tests have been performed on quasi-isotropic (0[degrees], +45[degrees], 90[degrees], -45[degrees])s carbon-epoxy laminates at different R-ratios. The typical damage development curve of this type of laminate lay-up is depicted in Fig. 1. From this curve it is said that the first increase in damage is caused by matrix cracking (stage 1). Matrix cracking is then followed by a steady state growth of delamination in stage II and in the final stage III fibers are broken leading to the final failure of the composite.

[FIGURE 1 OMITTED]

In Fig. 2, where the total counts of the "filtered" AE are depicted as a function of fatigue cycles of a fatigue test run at [S.sub.max] = 70% of the UTS and R = 0.3 at a frequency of 3 Hz, the shape of Fig. 1 can be identified as well. Elaborate studies of the damage development, using replica's, X-ray radiography with penetrant and ultrasonic C-scan identified matrix cracking in stage I, more matrix cracking and delamination in stage II and fiber failures in stage III. The filter consisted of only counting the AE counts in the upper 20% of the sinusoidal loading cycle, which hereby excludes the AE signals generated by crack closure.

[FIGURE 2 OMITTED]

Taking this evidence into account an attempt was made to model the fatigue behavior of the composite using the AE counts as an input.

For stage I it was obtained that: [n.sub.I] - [n.sub.ini] + [n.sub.f],

where [n.sub.ini] is the cumulative counts during quasi-static loading and of the cumulative counts during fatigue loading. Here, dn/dN = [A.sub.1] [[sigma].sub.max.sup.2c] / [cn.sup.c-1], so

[n.sub.I] - ([A.sub.1] [[sigma].sub.max.sup.2c] N + [n.sub.ini.sup.c])[sup.1/c]

and [A.sub.1] and c constants to be defined as a function of the applied stresses.

For stage II the variation of the cumulative counts with fatigue cycles [dn.sub.II]/dN = [A.sub.2] ([DELTA] [sigma]/[sigma]u).sup. [lambda]]

with [A.sub.2] and [lambda] the constants to be defined by plotting [dn.sub.II]/dN as a function of [DELTA][sigma]/ [[sigma].sub.u], being the fatigue stress amplitude divided by the tensile strength of the composite. It was found that for this composite [A.sub.2] = 32.3 and [lambda] = 10.2.

For stage III a problem arises when the final damage development is not located in between the monitored region as was the case for a 'certain' number of fatigue tests. [1]

The results presented indicate clear that the damage development in a composite can be monitored with the acoustic emission technique whenever the whole critical area is being monitored and the modeling of non-linear behavior (D) can be done based on AE data from the critical region when the modeling parameter is related to physical damage. Once the material behavior is characterized by such curves, monitoring the slope of the AE cumulative counts with the number of fatigue cycles can identify the damage stage the quasi-isotropic composite material is in.

AE system used: Vallen AMS-3 AE System, DECI 375M AE-sensors

Case study 2: Cure monitoring of concrete with AE

Another application of the use of AE is the on-line monitoring of curing processes in freshly poured concrete. The ultimate goal for monitoring concrete during its setting is to predict the long-term behavior of concrete and to interpret the dependence of this behavior on the initial composition and initial drying conditions of the concrete. [2]

Preliminary monitoring tests were performed in the laboratory, using a curing cell (Fig. 3) of 200 x 150 x 100 mm, instrumented with 2 AE sensors, a temperature sensor (T) and 2 ultrasonic transducers (U). Concrete of the following composition was used: 1.88 kg sand, 3.57 kg aggregates, 1.35 kg CEM I 42.5R, 0.45 1 water and 13.5 ml Rheolmild. The curing cell was placed in an environmental chamber (Heraeus-Votsch Climate Cabinet VLK 04/150) to isolate it from variations in outside temperature and humidity. However, the climate conditions were not controlled, allowing variations in temperature inside the concrete sample.

[FIGURE 3 OMITTED]

The curing process during the first 72 hours was monitored, measuring the cumulative counts of the AE events due to cracking, which is caused by the setting of the material during hardening. Apart from real curing related AE events, the received AE signals (Fig. 4a: original data) also consist of mechanical noise and EMI signals. The high frequency signals (Fig. 4b) can easily be distinguished from the low frequency AE events (Fig. 4c) and were filtered out by applying a low pass filter on both channels, generating the "filtered" AE data (Fig. 4a).

[FIGURE 4 OMITTED]

Simultaneously with the AE measurements, the temperature inside the concrete cell was registered and the evolution in the compressional wave velocity during hardening was derived from the time of flight of an ultrasonic compressional pulse, which was sent through the sample every 10 minutes. Figure 5(a) illustrates the correlation between the cumulative counts of the filtered AE events and the temperature variations inside the concrete sample. A first increase in the AE events after 9 hours corresponds with an increased inside temperature. Figure 5(b) shows that this increase in AE events occurs when the compressional wave velocity is almost stabilized. The first, and most important, increase in temperature is caused by an acceleration of the exothermal hydration reactions, which reach a peak value after 11 hrs, when the hydration reaction is at its maximum efficiency. This acceleration is also noticed in the evolution of the compressional wave velocity in time, indicating a stage of increased hardening. When the hydration reaction is at its maximum efficiency, the concrete sample is almost completely cured, and the velocity remains at a more stable level after 11 hrs. The first increase in the cumulative counts of the filtered AE events can thus be explained by intensification in microcrack generation at the final phase of the main concrete setting and shrinkage process (<20 hrs). At later stages, microcracking is caused by further curing and by temperature variations due to day/night cycles, which are also detected in the cumulative counts of the filtered AE events.

[FIGURE 5 OMITTED]

AE system used: Fracture Wave Detector, Digital Wave Co., VS375-M 340 AE-sensors

Case study 3: Corrosion monitoring of carbon steel and stainless steel by means of advanced data analysis and AI techniques on AE signals

Many failures in chemical process installations are due to corrosion processes. Literature [3] mentions that annually 3 to 5% of the G.D.P. in industrialized countries is lost directly or indirectly due to corrosion processes. Conventional non-destructive testing techniques are most frequently used to evaluate corrosion in chemical installations [4]. All these techniques are typically applied when the installation is taken out of service during periodic inspections. These techniques are insufficient when a considerable amount of corrosion has taken place between two periodic inspections.

A continuous corrosion monitoring system obtrudes. Acoustic emission offers the possibility to listen continuously to the materials while damage is developing. Typical AE signals could offer operators the alarm to inspect installations while damage is occurring. By immediate response to the occurring damage one can avoid that aggressive chemicals are liberated such that risk for people and environment is minimized.

This case study focuses on identifying 3 most important corrosion processes for chemical industry by means of the acoustic emission technique: uniform corrosion, pitting and stress corrosion cracking (SCC). Corrosion processes were induced on carbon steel (DIN 1.0038) and stainless steel (DIN 1.4541). The number of selected signals for each corrosion process is:

* 102 pitting events (brackish water + Fe[Cl.sub.3] 1% at 60[degrees]C)

* 43 SCC signals in stainless steel (Ca[Cl.sub.2] 40% at 80[degrees]C)

* 42 SCC signals in carbon steel (Ca([NO.sub.3])[sub.2] 60% at boiling point)

* 81 uniform corrosion signals in carbon steel ([H.sub.3][PO.sub.4] 10% at room temperature)

The generated AE signals in the different corrosion processes are highly dependent on the geometry of the material. In order to guarantee reproducible results and to facilitate up scaling from laboratory to chemical plant installations the research focused on a fixed geometry of a reference probe to monitor the signals. The corrosion damage on the probe needs to be representative for the damage on the chemical installation. This implies that the material of the probe underwent the same treatment as the material of the installation. Critical stress locations in the installations need to be identified, the same stresses need to be applied to the reference probe. The experimental set-up has been published in [5, 6].

Two reference probes are partially submerged in the corrosive environment. When transient signals are detected simultaneously these are likely due to system influences (flow of the chemical or pressurized air) or due to environmental influences (electromagnetic interference, mechanical vibrations, etc.). The third probe is placed in a non-corrosive environment. Simultaneous transients on all 3 probes are due to environmental influences. Transients that could be related to system or environmental influences are not retained within the data analysis.

Data analysis

Different steps were undertaken to achieve an automatic classification of the different corrosion processes. First, different filtering (denoising) techniques were tried out to increase the signal-to-noise ratio. The filtering techniques applied are: spectral subtraction [7], NLMS (normalized least-mean-squares adaptive filtering) [8], signal subspace filtering and an enhanced form of signal subspace filtering [9].

Second, the denoised time signals are transformed for feature extraction. The transformations applied are: FFT (fast Fourier transform), AR model (auto regressive model) [10] and a wavelet transform [11]. The FFT allows characterizing the time signals in frequency domain, the AR model in time domain and the wavelet transform in frequency-time domain. The best results in classification were obtained in frequency domain (FFT). 512 samples were selected in each time series and the FFT was applied. This results in 512 samples in the frequency domain. Only half of the spectrum is retained because real time series result in a symmetric spectrum, the second half contains no additional information.

Third, the information in the frequency domain, 256 frequencies with corresponding amplitudes, is compressed in a few parameters. The principal component analysis (PCA) [12] was applied to extract important parameters from the frequency domain. PCA finds projections of the data such that variance in data is maximal after each projection, with the constraint that projection axes are orthogonal.

In a fourth step the visualization of the data is done by estimating the probability density function (PDF) of the data in 2D space, after projecting the frequency data on the 2 most important principal axes (vectors) of the PCA. If different corrosion processes fall apart in different clusters, a classification machine should be successful in separating different corrosion processes. Different clusters can be detected by means of the maxima in the PDF. This allows determining qualitatively the quality of the transform on the time series, after PCA, independent of the chosen classification procedure. The density was estimated by means of an unsupervised neural network: Kohonen self-organizing map [12], see Figs. 6 and 7.

[FIGURE 6 OMITTED]

[FIGURE 7 OMITTED]

The kernel widths of the Gaussian kernels, which are centered on the neurons, are tuned automatically by using the least-squares cross-validation method [13]. Only for the FFT different corrosion classes resulted in separate clusters, see Fig. 7. The wavelet transform merged all data in I cluster. The AR model merged pitting and SCC in 1 cluster and resulted in a separate cluster for uniform corrosion. From this it can be concluded qualitatively that the FFT is the most appropriate transform to distinguish between different corrosion classes.

For quantitative assessment of separability, the least squares support vector machine (LS-SVM) [14] classification procedure is used. The function to be minimized in LS-SVM's is defined by:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (1)

such that [y.sub.k].([W.sup.T].[phi]([c.sub.k])+b)=1-[e.sub.k], k=1, ... ,N (2)

For a binary classification problem the classifier y(c) takes the form:

y(e) = sign([w.sup.T].[phi](c) + b) (3)

The w in (1) is related to the margin between different classes in the feature space. Minimizing 1/2[w.sup.T] .w results in maximizing the margin, see Fig. 8.

In SVM's the input space is defined by the data points: the [c.sub.k]. The feature space is defined by the [phi]([c.sub.k]), see Fig. 8. The 9-mapping is a transformation to a higher dimensional feature space. In this higher dimensional feature space a linear separation is made between the different classes.

[FIGURE 8 OMITTED]

The second term in (1) tries to minimize the classification error. When [lambda] is large less training examples will be misclassified, but this might result in over-fitting. The margin will allow for generalization. [y.sub.k] in (2) is the class label, for a first class this is chosen +1 for a second class -1. The 1 on the right hand side in (2) is a target value, and in case of overlapping distributions between classes some error [e.sub.k] is allowed on the target value.

The constrained optimization problem (1)-(2) can be solved easily by defining the Lagrangian of the problem:

L(w,b,e;a) = [J.sub.p](w,e) - [N.summation over (k=1)] [[alpha].sub.k] {[y.sub.k][[W.sup.T].[phi]([c.sub.k]) + b] -1 + [e.sub.k] (4)

In the general case the y-mapping can become infinite dimensional and a so-called dual problem is solved. The dual problem is defined by the search for the optimal Lagrangian multipliers: [[alpha].sub.1], [[alpha].sub.2] ... [[alpha].sub.N] of the constrained optimization problem, where a = {[[alpha].sub.1], [[alpha].sub.2] ... [[alpha].sub.N]}. Solving the primal problem corresponds to searching for the optimal w. It can be shown that optimal a is found by solving a set of linear equations. The advantage of the dual problem is to circumvent the problem of the possible infinite [phi]-mapping. The classifier in the dual space takes the form:

y(C)=Sign [[N.summation over (k=1)] [[alpha].sub.k].[y.sub.k].K(c,[c.sub.k]) + b] (5)

with K(c,[c.sub.k]) = [phi][(c).sup.T].[phi]([c.sub.k]) (6)

By solving the dual problem all cp disappear and are replaced by kernel functions K(c,[c.sub.k]) of the form in (6). This is known as the kernel trick. In practice different kernels can be chosen: linear, polynomial, radial basis function (RBF), etc. The RBF kernel takes the form:

K(c,[c.sub.k]) = exp (-[parallel]c Z -[c.sub.k] [parallel].sub.2.sup.2]/[[sigma].sup.2] (7)

For the classification problem both the linear kernel and the RBF-kernel were tried out, results were slightly better with the RBF-kernel. Classification results were always better for the data after the FFT compared to the wavelet transform and the AR-model. This could be expected while the FFT was the only transform that resulted in separate clusters for each corrosion class. The separation between different corrosion classes: SCC, pitting, uniform corrosion and time series in absence of corrosion could be achieved with about 2% classification error on average. The classification performances are obtained using the 10-fold cross-validation technique. Best results are obtained when using the spectral subtraction filtering technique. Figure 9 shows an example of the decision boundaries computed by the LS-SVM. Figure 10 shows an example of the decision boundaries computed from a Bayesian classifier. In this case density estimations are obtained from the four different classes and the c's are assigned to the class with maximum a posterior class probability:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (8)

[C.sub.i] represents the class and [[theta].sub.i] is the set of parameters that determines the class conditional density, P(c|[[theta].sub.i], [C.sub.i]). [theta] is the set of [[theta].sub.i], parameters: [theta]= {[[theta].sub.1], [[theta].sub.2], ... [[theta].sub.N]}, where N is the number of classes.

[FIGURE 9 OMITTED]

[FIGURE 10 OMITTED]

Density estimations are obtained by the minimum message length (MML) criterion for Gaussian mixture models (GMM) [15]. Observe that the decision boundaries obtained in Figs. 9 and 10 are very similar. There is no statistically significant difference in classification performance obtained between both classification paradigms for this particular data set.

All programs were partially developed in the C-programming language (Kohonen self-organizing map) and MATLAB. The LS-SVM was available as a software tool, developed at the Department of Electrical Engineering (ESAT), at the Katholieke Universiteit, Leuven.

AE system used: Fracture Wave Detector, Digital Wave Co., B1025 Digital Wave AE-sensors

Conclusions

Modeling the material behavior based on NDT data is gaining more interest as the confidence in the techniques increases, as the relation with physical damage is proven. The on-line monitoring of fatigue in CFRP laminates and of the curing of concrete with AE results in cumulative AE counts which after filtering offer the input data to model the damage behavior of the material. The AE technique offers parameters linked with physical damage to do the modeling with.

Sophisticated signal processing of the AE signals is also presented in order to extract features from the time signals; subsequently these features are used as input for a classification machine in order to identify different corrosion processes.

Case study 1: Research funded by the EC human capital and mobility fellowship,

Case study 2: Research funded by the FWO, Flemish Fund for Scientific Research,

Case study 3: Research funded by a Ph.D. grant of the Institute for the Promotion of

Innovation through Science and Technology in Flanders (IWT-Vlaanderen).

References

[1.] Tsamtsakis, D. and Wevers, M., "Acoustic emission to model the fatigue behavior of quasi-isotropic carbon-epoxy laminate composites", INSIGHT - Non-Destructive Testing and Condition Monitoring, Journal of British Institute of Non-Destructive Testing, 41 (8), August 1999, 513-516.

[2.] Lacouture, J.-C., "Modelisation de l'evolution des coefficients mecaniques du beton pendant la prise. Liaison avec les mesures ultrasonores lineaires et non lineaires" These, UFR De Physique, 2002, Universite Paris 7-Denis Diderot, 201p.

[3.] National Bureau of Standard/US Department of Commerce, "Economic effects of metallic corrosion in the United States", NBS, Special Publication 511-1, Washington DC, 1975

[4.] Rudlin, JR., "Review of existing and possible techniques for corrosion under insulation and wall thickness measurement in steel pressure containments", Insight, 39 (6), 1997, 413-416.

[5.] Winkelmans, M. and Wevers, M., "Non-destructive testing for corrosion monitoring in chemical plants", Journal of Acoustic Emission, 20, 2002, 206-217 (CD-ROM).

[6.] Winkelmans, M., Ph.D., "Fusion of non-destructive inspection techniques for corrosion monitoring in chemical process installations", Katholieke Universiteit Leuven, ISBN 905682-484-8, UDC 620.179:620.193.4, University of Leuven, 2004.

[7.] Boll, S.F., "Suppression of acoustic noise in speech using spectral subtraction", IEEE Transactions on Acoustics, Speech and Signal Processing, ASSP-27 (2), April 1979.

[8.] Moonen, M., Proudler, L, Introduction to adaptive signal processing.

[9.] Hu, Y., Loizou, C., "A subspace approach for enhancing speech corrupted by colored noise", IEEE Signal Processing Letters, 9 (7), July 2002.

[10.] Chatfield, C., The analysis of time series, Chapman & Hall, 1987.

[11.] Daubechies, L, Ten lectures on wavelets, Society for Industrial and Applied Mathematics, Philadelphia, Pennsylvania 1999.

[12.] Ham, F.M., Kostanic, L, Principles of neurocomputing for science & engineering, McGrawHill Higher Education, 2001.

[13.] Silverman, B.W., Density estimation for statistics and data analysis, Chapman & Hall, pp. 106-107, 1992.

[14.] Suykens, J.A.K., Van Gestel, T., De Brabanter, J., De Moor, J., Vandewalle J., Least squares support vector machines, World Scientific Pub. Co., Singapore, 2002.

[15.] Figueirido, M.A.T. and Jain, A.K., "Unsupervised learning of finite mixture models," IEEE Trans. On Pattern Analysis and Machine Intelligence, 24, (3) 381-396, March 2002.

MARTINE WEVERS (1), GERT VAN DIJCK (1, 2), WENDY DESADELEER (1), MARK WINKELMANS(3) and KOEN VAN DEN ABEELE (4)

(1) Katholieke Universiteit Leuven, Department of Materials Engineering, Kasteelpark Arenberg 44, B-3001 Leuven, Belgium; (2) Katholieke Universiteit Leuven, Computational Neuroscience Research Group, Campus Gasthuisberg, Herestraat 49, B-3000 Leuven, Belgium; (3) BASF AG, Carl-Bosch-Stra[beta]e 38, D-67056 Ludwigshafen, Germany; (4) Katholieke Universiteit Leuven, Afdeling Kortrijk, Subfaculty of Sciences, Sabbelaan 53, B-8500 Kortrijk, Belgium

Printer friendly Cite/link Email Feedback | |

Author: | Wevers, Martine; Van Dijck, Gert; Desadeleer, Wendy; Winkelmans, Mark; Van Den Abeele, Koen |
---|---|

Publication: | Journal of Acoustic Emission |

Date: | Jan 1, 2004 |

Words: | 3749 |

Previous Article: | Rods and tubes as AE waveguides. |

Next Article: | The effect of waveguide material and shape on acoustic emission transmission characteristics part 1: traditional features. |