1. Introduction
Structural deterioration occurs during the course of the structure’s lifespan. Detection of this deterioration at an early stage is invaluable for maintenance planning. Furthermore, this may allow an extension of the lifetime, which in turn reduces life-cycle costs. Thus detection of the deterioration has attracted interest in the field of structural health monitoring (SHM). A range of techniques making use of different damage-sensitive parameters such as vibration-based SHM [
1], strain-based SHM [
2], guided wave-based SHM [
3] etc. have been developed.
GW-based techniques have been popularly employed for plate-like structures, as they offer the possibility to reliably detect and localize the damage. Several different sensors such as the lead zirconate titanate (PZT), macro fiber composites, optical fiber sensors, acoustic non-contact transducers, as well as a laser Doppler vibrometer, have been used for the sensing. The optical fiber sensors such as the fiber Bragg grating (FBG) sensors and Fabry Perot interferometers offer low-weight alternatives for sensing [
4,
5]. The FBG sensors also offer several other advantages such as immunity to magnetic and electric fields, small size, ability to be multiplexed etc.
FBG sensors were traditionally used for strain and temperature measurement in the wavelength division multiplexing (WDM) approach wherein the change in the reflected wavelength is used to determine the measure. The sensitivity of the FBG sensors in the WDM configuration is not sufficient for GW measurements. So they are used in the edge filtering configuration. The edge filtering approach has been discussed in detail by Soman et al. [
6].
In addition to the edge filtering approach, the sensitivity of the FBG sensors can be further enhanced by employing it in the remote bonding configuration. The remote bonding configuration is shown in
Figure 1a. In the remote bonding configuration the Symmetric (S) and Antisymmetric (A) modes encounter the bond with the optical fiber and get converted into optical fiber longitudinal (L) and flexural (F) modes. The conversion of plate modes in the optical fiber is shown in
Figure 1b. The phenomenon is discussed in detail in [
6,
7]. The optical fiber acts as an excellent wave guide for the mechanical waves and these fiber modes travel along the fiber with minimal attenuation. Along the fiber, as the FBG is in free condition (without constraints), it undergoes greater deformation leading to a higher sensitivity. In the directly bonded configuration, a shear lag effect is observed which results in a lower response in the directly bonded FBG than in the remotely bonded FBG. The phenomenon was studied in detail by Huang et al. [
8]. Wee et al. quantified the improved sensitivity and reported that the remote bonding configuration yields an 81% increase in sensitivity for some conditions [
9].
Some researchers have leveraged the added sensitivity of the FBG sensors in the remote configuration for damage detection [
9,
10]. Yu et al. [
11,
12] proposed the use of FBG in remote configuration for sensing in extremely high temperatures. The remote FBG may be paired up with non-contact excitation techniques and actuation and sensing may be achieved in spite of the high temperature of the structure to be monitored. The same non-contact actuation was paired up with a remote FBG in an improved adhesion configuration [
13]. The remote FBG is located equidistant from two perpendicular bonds on the same fiber. The two bonds allow overcoming of the directional sensitivity of the FBG sensors. The full wavefield information is then used for detecting damage in the structure. In both the described studies, the full wave field measurements utilized make the approach time-consuming and not suitable for in-service SHM. In the studies using sparse arrays of sensors or actuators, work by Wee et al. [
14] which uses the remote bonded FBG for self-referencing stands out. They use two bonds equidistant from the remote FBG on either side. The interference of the waves coupled at the two bonds is then used for detecting damage. This method is limited only to damage detection in a limited region of the structure. For damage localization, Yu et al. [
11]. used the difference in the arrivals of the A0 and S0 wave from an acoustic event. This allowed determination of the distance between the acoustic event and the bond location, but is not useful for knowing the damage location precisely. Hence based on the literature review, there is no work where the remotely bonded FBGs are used for damage localization. This paper addresses this lack of research in the area and presents a proof of concept for the use of FBG sensors in the remote configuration for the damage localization.
The work builds on the work by Soman et al. [
15,
16]. who used the ellipse mapping algorithm for damage detection and localization using directly bonded FBG sensors. Soman et al. [
15,
16] discretized the plate into smaller pixels and calculated the damage index (DI) for all the pixels using a brute force method. This technique leads to additional computational costs and delays in the measurements and the determination of the condition of the structure. In order to reduce the computational load. Balasubramaniam et al. [
17]. proposed a multi-level approach where the regions of the damage were identified in the first step and the damage localization was carried out only in the identified region. This approach indeed reduces the computational load but is not always possible as it needs a large number of actuator-sensor pairs. The use of additional actuators or sensors increases the cost significantly and hence is not desirable and may not be possible on large structures. Hence this paper proposes a damage location algorithm based on particle swarm optimization (PSO). The damage localization problem is posed as an optimization problem with the objective to maximize the damage index. The paper was inspired by similar work by Chen et al. [
18,
19] where the evolutionary algorithm was used for damage localization. The key difference with the above work is the development of a different objective function, which for the current problem using FBG sensors in remote configuration performs better than the objective function developed by Chen et al. [
19].
The paper provides a proof of concept for the use of FBG sensors in the remote bonding configuration for damage localization. A PSO-based algorithm is developed to determine the damage location while reducing the computational load as compared to the damage mapping technique. The paper also presents a new objective function for damage localization and compares it with a metric previously developed in the literature. The validation is performed on a simple aluminum plate with simulated damage. The presented results indicate that indeed the PSO-based algorithm works efficiently and successfully localizes damage with reduced computation time. The reduction in computation time makes the method more suitable for real-time SHM.
The rest of the article is organized as follows: the next section highlights the methodology for damage localization including the main features of the PSO, the objective functions, and their implementation. The
Section 3 outlines the setup used for the validation of the methodology. The
Section 4 presents the results for the damage localization along with some sensitivity and comparative studies. The last section draws some conclusions based on the presented results and identifies areas for further research.
2. Methodology
Damage localization through an ellipse imaging method can be done with a pixelated approach, where a damage index is calculated on every discrete pixel on a grid using the Equation (
1) [
20]. This approach is termed the brute-force approach.
where
is the damage index at a particular pixel,
is the number of actuator-sensor paths,
is the timing difference between the measured and expected time-of-flight (TOF), and
is a sensitivity factor. This approach has been shown to be effective but is computationally expensive and time consuming.
In order to develop a computationally efficient approach, damage localization is treated as a function optimization problem. As a result, the damage may be localized in fewer damage index (DI) calculations. This method involves creating an objective function that returns the likelihood of a particular location being the damaged location, and using existing mathematical optimization tools to find the location that results in the highest probability.
Particle swarm optimization (PSO) is a common optimization technique that is by design simple and computationally efficient. PSO involves a group, or swarm, of candidate solutions, known as particles, that have a position and velocity [
21]. Each particle will move a small amount every iteration, attracted to both its individual best position it has found so far, and the overall best position in the swarm. The attraction to each particle’s individual best is intended to cause particles to find and fully explore the search space, and the attraction to the global best is intended to achieve convergence. This tendency for each particle to interact with its own previous locations and the best performing location gives PSO algorithms an element of swarm intelligence not found in techniques like genetic algorithms or gradient descent. For a 2D localization problem, the PSO is ideally suited and has fewer artificial constraints. This allows a more thorough search of the possible solution space and efficiency in the implementation. Hence the PSO was chosen as the algorithm for the optimization in this paper.
2.1. Particle Swarm Optimization and Its Implementation
The basic procedure for a PSO algorithm is as follows. The flowchart is provided in
Figure 2.
1. Initialize Swarm: A population of particles called a swarm is created, with positions randomly distributed uniformly throughout the search space. The particle’s fitness is found by passing its position into the objective function up for optimization. In our research, standard elliptical imaging TOF methods are used to find ellipses focused on each actuator-sensor (AS) pair, and the fitness is based on the particle’s proximity to these ellipses which corresponds to the likelihood that the particle’s location is the damaged location. The particle with the highest global fitness is found and its location recorded.
2. Update Particle Positions: The change in a particle’s position from one iteration to the next is known as the particle’s velocity. Each particle’s velocity is determined at each iteration, and then added to the particle’s position to get the particle’s next position. Each particle’s velocity is determined as a weighted vector sum of 3 vectors: the particle’s previous velocity, the displacement vector to the location with the highest fitness that particle has found so far, and the displacement to the location with the highest global fitness found so far. The weight given to its previous velocity is known as its momentum coefficient, and the weights given to the individual maxima and global maxima are known as acceleration coefficients. Adjusting these acceleration coefficients changes each particle’s attraction to its individual and global maxima. Note that on the first iteration, as each particle’s individual maxima is its current position, the particle with the highest fitness does not move.
3. Evaluate Fitness: The fitness is evaluated for the new position of each particle. The fitness value is then used to determine the global maximum and the individual maximum which are then used to check for convergence as well as determination of the weighting coefficients. The fitness function used in this research is discussed in more detail later.
4. Check Terminating Condition: If the maximum global fitness value exceeds a certain criterion, or if the average particle velocity is below a certain threshold, or if the maximum number of iterations has been reached, the algorithm returns the position with the best fitness found as the optimal solution.
The number of particles, also known as swarm size, along with the momentum and acceleration coefficients, are input factors that need to be tuned for best performance.
2.2. Objective Function
The objective function is based on the candidate location’s proximity to ellipses focused on each actuator-sensor pair. The measured time of flight (TOF) is compared with the expected time, and the difference is passed through a negative exponential function as shown in Equation (
1). To combat the less than perfect overlap of ellipses from separate actuator-sensor (AS) pairs, this value is calculated for the particle’s position as well as a few positions nearby. These values are averaged to provide the final particle fitness. This helps to smooth the steep slope of the negative exponential scaling function when the timing difference is very small, however, it comes at the cost of a more computationally expensive objective function, which usually must be counteracted by using fewer particles.
In addition to the above-mentioned objective function, one other objective function was studied. The objective function is based on the work by Chen et al. [
19] who used the number of ellipses passing through a point on the structure as an objective function. They applied this metric on a sensor network consisting of 6 actuator sensors, leading to 15 AS pairs. The suitability of the metric for the damage localization in our case is studied and compared with the exponential metric.
2.3. Experimental Setup
The experimental validation was conducted on a 61 cm × 60.5 cm × 0.08 cm 6061 aluminum plate, with weights added to simulate damage. The additional weights offer a repeatable and reversible way of introducing wave scatterers in the structure and are commonly used. The added mass is a scatterer for the waves and is considered equivalent to drilled holes [
23,
24]. The weights used were rare-earth magnets on either side of the plate to ensure firm attachment. The magnets were 2 cm in diameter, and weighed a total of 25 g. Three piezoelectric transducers (PZT) for wave excitation were bonded to the plate using cyanoacrylate adhesive. The PZTs, driven by an arbitrary waveform generator (Agilent 33500B Series) and voltage amplifier (Krohn-Hite Model 7602M Wideband) set to an output voltage of 60 V, generated bursts of a 300 kHz excitation signal to generate
and
Lamb waves. The excitation signal was a 5.5 wavelength Hanning windowed sine wave.
3 polyimide-coated optical fibers were bonded radially to the first PZT at a distance of 16 cm in accordance with the procedure by Wee et al. [
9], with fiber Bragg grating (FBG) sensors in the remotely bonded configuration in each fiber. The distance from the fiber bonding point to the FBG is different for each fiber, and are 31 cm, 71 cm, and 80.25 cm respectively. The ends of each fiber were placed in index matching gel to reduce reflections from the end of the fiber. A point to note is that the aim of the experiments was to provide a proof-of-concept of the method. The deployment of the number of actuators, sensors, and their location was based on prior experience of the authors. In reality, the deployment of the actuators and sensors needs to be optimized to ensure reliable performance. Several techniques tackling this problem may be found in the literature [
15,
24,
25]. The optimization of the network is identified as an area of future work.
A tunable laser source (Photonetics/NetTest TUNICS-Plus) is used for generating the laser wavelength at the midpoint of the rising edge of the FBGs, and the signals reflected from the FBG were captured using a circulator, photo-detector (Newport 1544-A Near-IR Fiber-Optic Receiver), and oscilloscope (Keysight DSO5032A 5000 Series). The photo of the setup is shown in
Figure 3 and the schematic of the experimental setup is shown in
Figure 4.
For each of the 7 damage cases as well as the healthy condition, the FBG sensor response is recorded for each actuator-sensor pair, with only one actuator active at a time. In order to remove any biases or high-frequency noise, a bandpass filter in the range of 150 kHz to 450 kHz was applied to the raw data. The filtered signals were then used for further processing. For each case the residual signal, which is the difference between the healthy case and each damaged case is found. The Hilbert envelope is fitted on the residual and smoothed. The envelope and smoothing allow us to overcome minor misalignment in the healthy and damaged signal. This leads to a more repeatable and robust identification of the peaks. The peaks in the smoothed and enveloped residual are identified using the ’findpeaks’ command in MATLAB. These peaks are used in the elliptical fitness function described previously.
A point to note is that the optical fibers show directional sensitivity, when the excitation is in an axial direction to the fiber, coupling of the A and S modes is more dominant. The domination of the coupling of S and A modes lasts until around 75°. When the actuation is more than 75° or near perpendicular to the adhesive bond the coupling of the SH waves into the optical fiber is more dominant. When the SH waves are dominant, the identification of the direct arrival of the A0 wave is difficult. This makes signal processing more challenging. So, only the AS pairs that lead to A and S mode coupling are chosen for the signal processing namely A1 with all three Bonds (Bond 1, Bond 2, and Bond 3), A2 with Bond 3, and A3 with Bond 1. The A0 mode is chosen as it has been shown to be better for damage localization. The A0 wave has a smaller wavelength and hence may be used for detecting small damage.
4. Conclusions
GW-based SHM using FBG sensors is a field that is receiving renewed interest in the last few years. The paper provides a proof of concept of the use of FBG sensors in remote configuration for the localization of damage in an aluminum plate. To the best knowledge of the authors, this work is the first time remotely bonded FBG sensors have been used for damage localization. In order to reduce the time needed for the localization of the damage, a PSO-based algorithm has been implemented making use of an innovative objective function. The parameters for the PSO have been chosen based on engineering judgment to ensure a thorough global as well as local search. The PSO indeed allows the damage localization in significantly lower time as opposed to the brute force algorithm.
The performance of the proposed objective function is qualitatively and quantitatively better than the objective function proposed in the literature. The objective function based on the number of ellipses works well where a large number of AS pairs are present. As the FBG sensors are passive sensors, the number of AS pairs are severely reduced. Hence, the new objective function developed takes into consideration the passive nature of the FBG sensors and hence is more suitable for use with the FBG sensors. It should be conceded that the study is applied for the proof of concept on a simple aluminum plate, additional complexities such as anisotropy of the structure and the presence of more complex structural features need to be taken into account before the above technique may be used for damage localization in real structures and is identified as an area of future work. Although the proposed cost function yields the correct damage location, there is a chance that the PSO gets stuck in local maxima, in order to avoid such scenarios, a multi-objective optimization needs to be developed and is planned for future work.