Next Article in Journal
Two-Field Excitation for Contactless Inductive Flow Tomography
Previous Article in Journal
Advancing DIEP Flap Monitoring with Optical Imaging Techniques: A Narrative Review
Previous Article in Special Issue
Development of a Predictive Model for Evaluation of the Influence of Various Parameters on the Performance of an Oscillating Water Column Device
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Natural Frequency Transmissibility for Detection of Cracks in Horizontal Axis Wind Turbine Blades

Department of Mechanical and Industrial Engineering, University of Toronto, 5 King’s College Road, Toronto, ON M5S 3G8, Canada
*
Author to whom correspondence should be addressed.
Sensors 2024, 24(14), 4456; https://doi.org/10.3390/s24144456
Submission received: 18 March 2024 / Revised: 2 July 2024 / Accepted: 4 July 2024 / Published: 10 July 2024
(This article belongs to the Special Issue Real-Time Monitoring Technology for Built Infrastructure Systems)

Abstract

:
Defects on horizontal axis wind turbine blades are difficult to identify and monitor with conventional forms of non-destructive examination due to the blade’s large size and limited accessibility during continuous operation. This article examines both strain and acceleration transmissibility as methods of continuous damage detection on wind turbine blades. A scaled 117 cm offshore wind turbine blade was first designed, 3D printed, and modelled numerically in ANSYS. Transverse cracks were deliberately introduced to the blade at 10 cm intervals along its leading edge. Subsequent changes in the transmissibility, relative to an undamaged baseline model, were measured using different variable combinations at the blade’s first three natural frequencies. Experimental results indicated that strain transmissibility was able to locate a 1.0 cm defect at a range of 70–110 cm from the blade hub using the amplitudes of the first natural frequency of vibration. The numerical model was able to simulate the strain experimental results and was determined to be valid for future defect characterization. Acceleration transmissibility was unable to experimentally identify defects sized at 1.0 cm and below but was able to identify 1.0 cm sized defects numerically. It was concluded that transmissibility is viable for continuous damage detection on blades but that further research into other defect types and locations is required prior to conducting full-scale testing.

1. Introduction

1.1. Objectives

This research focuses on the use of natural frequency transmissibility as an online method of damage detection on horizontal axis wind turbine blades. Given a wind turbine blade’s large size and limited accessibility while in service, challenges exist with online monitoring and defect characterization while damage is still in its infancy. It is important for the wind turbine blades to be continuously monitored as many damage mechanisms for wind turbines are random in nature and are not time-based. Extensive repairs, replacements, and wind turbine failure can be mitigated with an appropriate structural health monitoring program combined with timely maintenance.
Our investigation centred on changes in the vibration transmissibility between pairs of sensors due to a transverse defect added to the leading edge of the blade. Single defects were induced at different locations along the blade’s length and then assessed using both strain and acceleration transmissibility. The transmissibility changes between pairs of sensors were assessed as functions of several different variables, including modal frequency, signal type, and sensor orientation. The intent was to determine parameter settings that maximized the damage detection capabilities of transmissibility. This was completed by comparing the results of different variable combinations using a confusion matrix analysis. The objective of this research was to determine an approximate range of defect size and location on the blade for which the defect could be recognized with the best performing variable combination. An additional objective of this research was the development and validation of a numerical model to allow for additional types of damage to be characterized by transmissibility methods in future work.
This project expands on past transmissibility research by focusing on the transmissibility changes at each natural frequency individually as opposed to utilizing the entire frequency spectrum or a summation of the natural frequencies. It develops a damage identification algorithm to decide if a defect has occurred and determine the location. It provides an analysis of the effect of the defect location on the ability to identify the damage with transmissibility. It also examines the effects of the sensor orientation as well as the signal type used for analysis. Furthermore, this project utilizes smaller, more realistic-sized defects that would not be immediately detrimental to blade integrity and would be difficult to identify with conventional structural health monitoring methods.

1.2. Background and Literature Review

1.2.1. Wind Turbine Blades

Horizontal axis wind turbines (HAWT) have become a symbol of renewable energy generation in Canada and across the world. In 2022, approximately 6.6% of the electricity demand in Canada was supplied by wind power [1]. Many projects are currently in development to increase the amount of energy produced by wind in Canada. This includes onshore projects that are currently under construction such as the 200 MW Buffalo Atlee project in Alberta [2] and the 200 MW Bekevar Wind Facility in Saskatchewan [3]. The provinces of Nova Scotia and Newfoundland have recently made steps toward the development of offshore wind in Canada by amending existing legislation to regulate offshore renewable energy projects [4]. Nova Scotia has also announced plans to provide leases for 5 GW of offshore wind energy by 2030 [5]. Wind turbines are, therefore, an important part of Canada’s current and future energy outlook. A typical large-scale horizontal axis wind turbine is shown in Figure 1a.
Many wind turbines used in large-scale energy production, such as Vestas [6,7], Siemens Gamesa [8], GE [9], and Enercon [10], are designed for a 20–30 year lifespan before decommissioning, replacement in kind or repowering is required. In Canada, approximately one-third of the wind turbines listed in the Canadian Wind Turbine Database are set to reach an age of 20 or older by 2030 [11]. As turbines approach their end of life, the components are more likely to experience damage and degradation. Inspection and timely maintenance of wind turbine components are therefore important to identify and mitigate damage to keep the turbine online and produce electricity.
Wind turbine blades are prone to damage, which can affect the overall availability and efficiency of the turbine over time. The author in [12] referenced blade data from the insurer GCube, stating that blades had an annual failure rate of 0.54%, with over 40% of wind turbine insurance claims occurring as a result of blade failure. Replacement blade costs are estimated at USD 80,000–USD 300,000 as of 2023 [13]. Blade repair costs are significantly cheaper at estimates of USD 20,000 for a typical repair. Repair cost is, however, dependent on the severity of the damage.
Several reviews of wind turbine blade damage types and causes were conducted in [14,15,16,17]. Common damage seen in blades includes leading edge erosion, leading and trailing edge bond splits, spar cracking/debonding, fatigue cracking, and delaminations. Damage can occur resulting from regular operation, with leading edge erosion and fatigue cracking probability becoming more prevalent the longer the blade is in service. Operating incidents, including high wind events, hail, icing, lightning, debris and bird strikes, may also contribute to blade damage. Some defects, such as delaminations, voids, fibre misalignment, lack of fusion, kissing bonds, and wrinkles, can occur in the initial manufacturing process [18]. These imperfections can propagate and become problematic once in service. Defects can be both external or internal to the blade, as shown in Figure 1b [19].
Figure 1. (a) Horizontal axis wind turbine at a wind farm near Pincher Creek, Alberta, Canada. (b) Typical defects shown on a wind turbine blade’s cross section. Reproduced with permission from [19]. Katnam et al., Journal of Adhesion; published by Taylor & Francis Ltd., 2014.
Figure 1. (a) Horizontal axis wind turbine at a wind farm near Pincher Creek, Alberta, Canada. (b) Typical defects shown on a wind turbine blade’s cross section. Reproduced with permission from [19]. Katnam et al., Journal of Adhesion; published by Taylor & Francis Ltd., 2014.
Sensors 24 04456 g001

1.2.2. Inspection Methods

Vestas reports that an inspection costs 0.001 percent of a blade replacement and can reduce repair costs by upwards of 75% if the damage is caught early [20]. The majority of inspections that occur on wind turbine blades in-service are visual. Visual inspections can be completed with technicians via crane or rope access. Inspections can also be remotely completed with high powered cameras and drones. Industrial companies have been able to recognize defects as small as a few millimetres with remote visual methods [21,22]. Artificial intelligence has also been used in research to identify defects from remote visual inspection using image recognition [23]. Although effective for the detection of surface damage both externally and internally, visual inspections can be time consuming and cannot detect subsurface damage within the blade’s structure.
Passive thermography can examine the external and internal components of the blade via camera or drone and can be effective at identifying many different types of defects, especially when combined with visual inspection. The research group in [24,25] examined a blade segment exposed to different ambient conditions. Their work was successful at recognizing different types of defects using thermograms, but they noted the dependence on time of day, temperature and weather conditions. Different types of defects were individually more or less apparent on the thermogram, depending on these parameters at the time of the inspection. Multiple inspections would, therefore, be required throughout the day in ideal weather conditions to properly examine the blade for all damage types.
Ultrasonic testing (UT) is another method that can be used for defect detection on wind turbine blades. This includes conventional pulse-echo and phased array [26], guided wave [27] and automated methods [28]. Although UT is well-established and can identify many types of defects, ultrasound is impractical in service, given the amount of time required to perform a complete blade scan. Tests of the conventional pulse-echo in [26] found an average scan rate of 11:48 m i n / f t 2 on composite panels representing different parts of a blade. Phased array and automated methods would be much faster to use but would still take significant time, given the size of most wind turbine blades. Guided wave UT was effective at detecting disbonds on wind turbine blades, as shown in [27], and is ideal for longer distances. Guided wave UT is still primarily in the research phase for wind turbine blades. In general, UT is more prevalent in industry for the identification of factory defects with ultrasound before the blade is put into service.
The above non-destructive examination methods, although effective, only represent the blade’s condition at a moment in time. As a wind turbine is in constant operation and many of its damage mechanisms are not time dependent, a continuous inspection method is preferable. Continuous methods of blade monitoring include acoustic emission, auditory monitoring with microphones and vibration. SCADA data analysis can also be used to scan abnormalities in the blade using operating data, for example, blades exposed to lightning strikes [29].
Acoustic emission (AE) has been the subject of numerous research studies, such as [30,31]. There are also commercial products available that use AE on blades, such as Sensoria, offered by MISTRAS [32]. Although effective as a continuous monitoring method at identifying propagating defects as well as defect location, AE struggles with assessing defect characterization and severity. It is also prone to noise from the environment and requires complicated signal processing methods to assess the signal [33].
Research in [34] examined passive acoustic monitoring in wind turbine blades. This method utilized microphones measuring changes in the acoustic transmissions within the structure when exposed to different through-wall defects. This method was able to identify different severities but noted different results for a defect based on its location on the airfoil’s cross section. This research is promising for a structural health monitoring method as it is low-cost and easy to implement. This method would be reliant on defects being through-wall, so damage severity may be more severe at the time of discovery versus other methods.
Vibration condition monitoring systems are available in the industry for blades, such as Weidmüller’s BLADEcontrol system [35]. Most of the research on vibration monitoring of wind turbine blades is passive-based and relies on excitation occurring due to the wind. In [36,37], operational modal analysis techniques were used for damage detection in blade simulations. Active vibration monitoring using an actuator was also conducted on an operating Vestas turbine in [38] and was able to recognize a trailing edge crack.
Vibration condition monitoring programs often focus on monitoring the amplitudes of a raw signal or changes to the natural frequencies of the component. Identification of large changes in the amplitude of a raw vibration signal can provide protection from operating events as the turbine can be shutdown to minimize significant damage. Large structures, such as wind turbine blades, tend to have minimal natural frequency changes at the lower modes when exposed to damage. Damage detection is, therefore, difficult unless the damage is severe. Furthermore, locating the defect can be a challenge with more common vibration methods.
The inspection methods discussed above have their limitations. Visual inspection, thermography and ultrasound are primarily discontinuous and periodic inspections may fail to capture blade damage until it has propagated to a major repair or failure. Continuous methods such as SCADA analysis, acoustic emission, and vibration are effective but require significant damage to alert the operator of an issue. Characterization of the defect is also difficult with continuous methods. AE is the only method of the three that can determine location but is prone to environmental noise. A subset of vibration, transmissibility, aims to provide better sensitivity to damage than traditional vibration analysis using amplitudes and natural frequencies. Transmissibility is a continuous method of NDE, making it highly advantageous for equipment in operation. It aims to identify and localize the defect and provide that information to the wind turbine operator to assess the risk to the blade’s integrity.

1.2.3. Transmissibility Theory of Single and Multiple Force Functions

The transmissibility theory of multi-degree-of-freedom (MDOF) systems, as well as the application to damage detection, is discussed extensively in [39,40,41,42,43]. Vibration transmissibility describes the relationship between an output response and an input function. In an MDOF system, transmissibility is calculated by determining the ratio of the output responses between two locations on the same structure. Research has shown that the transmissibility between a pair of sensors can be independent of the applied force input; this suggests its potential as a tool for non-destructive evaluation (NDE) of in-service engineering structures.
Consider the wind turbine blade in Figure 2. In this figure, X i ω and X j ω are the measured output responses at positions i and j resulting from an applied forced F k ω at position k.
X i ω = H i k ω F k ω
X j ω = H j k ω F k ω
T i j k ( ω ) = X i ( ω ) X j ( ω ) = H i k ( ω ) F k ( ω ) H j k ( ω ) F k ( ω ) = H i k ( ω ) H j k ( ω )
The transmissibility, T i j k , is defined as the ratio of the responses between i and j in Equation (3) for an MDOF system. By dividing the two responses, the force F k ω cancels out. The transmissibility, in this case, is, therefore, a ratio of the frequency response functions (FRF) or transfer functions, H, of the material at i and j with respect only to the location of the force F k ω .
This method is useful in cases where the force magnitude is constantly changing, but the location is a known constant. However, the above method only applies to singular forces, so it is limited in its potential applications. Problems arise with this method when the location is unknown, changing, or multiple forces are used. Consider the wind turbine blade in Figure 3, which has two forces applied at two different locations.
X i ω = H i k ω F k ω + H i l ω F l ω
X j ω = H j k ω F k ω + H j l ω F l ω
T i j k l ω = X i ω X j ω = H i k ω F k ω + H i l ω F l ω H j k ω F k ω + H j l ω F l ω
The transmissibility between i and j in this scenario is much more complex as force cancellation does not occur in Equation (6). Information is required on the applied force magnitude and location for transmissibility to be used for damage detection in this scenario.
In [44], the transmissibility was investigated under different force application locations. This research found that when the transmissibility between two sensor outputs was plotted in the frequency domain under the same undamaged condition but with different forces, the transmissibilities converged to the same value at the natural frequencies of the system. This was found to be true using both singular and multiple forces. By restricting the transmissibility analysis to the natural frequencies as in Equation (7) below, the transmissibility becomes completely independent of the magnitude of the applied forces. In addition, the FRFs also become independent of the applied force locations k and l. A system with multiple unknown forces can now be assessed using only the output responses, making this ideal for online damage detection.
T i j ω n = X i ω n X j ω n = H i ω n F k ω n + F l ω n H j ω n F k ω n + F l ω n = H i ω n H j ω n

1.2.4. Damage Detection Using Transmissibility

Transmissibility is dependent on the local mass, stiffness, and damping of the material, making it ideal for damage detection. Changes to the structure due to the presence of flaws would result in a localized change in the transmissibility. Damage detection using transmissibility methods is of interest due to its potential to assess both linear and non-linear types of defects [45]. There are many methods that can be used to assess the changes in transmissibility resulting from a defect. Some examples include the Damage Indicator method [46], Transmissibility Damage Indicator (TDI) [39], and the Damage Index method [45].
Many damage detection methods utilize controlled forces and a large range of frequencies in their research. As such, research often follows the theory described in Equations (1)–(3). As a follow-up to their work in [44], [47] examined logarithmic transmissibility changes using small frequency bands around the natural frequencies. A defect on a cantilever beam was identified using this method at different vibration modes. However, it was noted that some modes were more amenable to identification of the defect than others. The overall importance of frequency band selection was also noted in this work. During numerical testing, different forces applied to the healthy beam revealed that false positive indications of a defect appeared in wider frequency ranges that were not narrowly focused on the natural frequencies. These false positives masked the defect location once the defect was implanted in the structure.
Different types of transmissibility can also be used to assess damage. Much of the literature discussed in this paper focuses on acceleration transmissibility, although other types, such as velocity [48] and strain [49] transmissibility, have also been used. The researchers in [49] argued that strain transmissibility is more sensitive than the other types due to the strain’s relationship with curvature.
Other researchers have focused on the relative merits of various types of signal analysis with transmissibility. Signal amplitudes are predominantly used in the literature. The phase angle was used successfully, along with the signal amplitude in [46], to identify a 5% stiffness reduction added to different elements of a cantilever beam simulation. Research has also been conducted to remove the requirement of transmissibility analysis being conducted in the frequency domain. In [50], strain transmissibility was used to determine damage in the time domain on a simply supported beam set-up. Wavelet transforms have also been used in transmissibility analysis for damage detection in [51].

1.2.5. Transmissibility in Wind Turbine Blades

The transmissibility concept has been applied to damage detection or the characterization of wind turbine blades in a number of different ways. Early research into the topic in [48] assessed velocity transmissibility on a blade using a scanning laser Doppler vibrometer over a large frequency bandwidth. This research was able to recognize an added mass simulating damage but noted a lot of random noise present in the frequency response functions.
A full-size blade undergoing static loading and fatigue testing at the first natural frequency was examined in [52]. An adhesive joint debond between the spar and the shell was successfully identified with accelerometers and FBGs using transmissibility changes of adjacent sensors. This method was also able to monitor the defect’s propagation as it increased in severity. In this experiment, only one defect was present, located in the middle of the blade.
Wavelet energy transforms were used for damage detection on wind turbine blades in [53]. This method was able to determine that two of the three blades tested were damaged, but it did not specify the damage location. These blades were also used in [54] to compare four different system identification models using transmissibility with a similar outcome.

2. Materials and Methods

2.1. Experimental Set-Up

2.1.1. Blade Design and Fabrication

A wind turbine blade was designed based on the NREL 15MW reference turbine [55]. This blade was scaled from its original length of 117 m to 117 cm (1:100 scale) to maintain practicality for experimental laboratory testing. Additional changes to the scaled blade design included increasing the chord and wall thicknesses to assist with fabrication. A baseplate was also designed to attach the blade to an optical table for vibration testing.
The blade was 3D printed using the fused filament fabrication (FFF) method. This method allowed for convenient and repeatable fabrication and testing. FFF thermoplastic parts are also easily repairable so that one blade can be repaired and reused many times during experimental testing. The blade was fabricated with a Raise3D Pro2 printer [56] and Raise3D Premium PLA filament [57] (Raise3D, Lake Forest, CA, USA). It was printed in the longitudinal direction to best simulate uniaxial glass fibre. A 0.1 mm layer height and 100% infill were used to best represent the solid material used in the numerical model. The fully fabricated blade assembly is shown in Figure 4.

2.1.2. Blade Joint Adhesion

Due to limitations with print build size, the blade was fabricated in five separate pieces, which were then attached together. Extensive coupon testing was completed to ensure that a blade composed of multiple attached components performed similarly to a blade fabricated in one piece. As the blade was to be subjected only to elastic vibration testing, Young’s modulus (and density) were the key material properties required for the development of an ANSYS numerical model. All other parameters, such as temperature, were assumed to be constant. Full-size 3D printed coupons were compared with half-sized coupons attached together with adhesive; the objective was to confirm the rule of mixtures and that a blade printed as a single unit would show the same vibration characteristics as a blade composed of several segments attached together.
Preliminary testing found that several factors affected the characteristic Young’s modulus of the lab specimens consisting of two half-coupons joined together. The alignment of the half coupons was important in achieving a consistent cross-sectional area and symmetrical stress distribution. The surface roughness, cleanliness, and contact force at the specimen joint were crucial in achieving a sufficient bond. Print location and platform support type were important in maintaining print quality and minimizing variations among coupons. Therefore, all coupons were printed in the same location on the printer with the same ”raft” platform support to minimize differences between prints. A consistent optimal adhesion procedure was also developed using cyanoacrylate glue [58]. The contact surface was first roughened with sandpaper and cleaned with alcohol to remove contaminants. The contact surfaces were assembled and aligned with elastic bands on a backing plate. The adhesive was applied and cured for 24-h under constant force from the elastics. After the elastics and backing plate were removed, a second layer of adhesive was applied to ensure the contact surfaces were fully covered. Once cured, the adhesive cap was sanded until flush with the base material.
Tensile testing was performed on 5 single-print coupons and 6 glued coupon sets. Figure 5 presents the resultant specimen-averaged Young’s modulus for each coupon. Young’s modulus had averages of 2984 MPa for glued coupons and 3013 MPa for full-sized coupons. The glued coupons had a larger variation in Young’s modulus results, with a standard deviation of 197 MPa compared to a standard deviation of 100 MPa for the full-sized coupons. This was due to imperfections implanted in the blade during the gluing process. Overall, these data show that the adhesion procedure results in no significant change to Young’s modulus compared to a single-print specimen, with a slightly larger variation with the glued samples.

2.1.3. Sensors

Two different types of sensors were used in the experimental set-up: eight MPU-6050 micro electro-mechanical system (MEMS) accelerometers [59], positioned 10 cm to 80 cm away from the hub at 10 cm intervals; and eleven fibre Bragg grating (FBG) strain sensors, positioned 10 cm to 110 cm away from the hub at 10 cm intervals. An additional FBG was attached to a 3D-printed coupon separate from the blade for temperature compensation. The accelerometers were attached to the midpoint of the blade’s chord with elastic bands and 3 M tape. The FBGs were spliced together into a single fibre, which was placed in a groove at the midpoint of the blade. The splices were protected with a 3D printed splice cover and attached to the blade with 3 M tape next to the FBGs. Each FBG was adhered to the blade using paraffin wax and tape. Figure 6 and Figure 7a,b show the sensor layout and typical sensor assemblies, respectively. This set-up allowed an assessment of acceleration and strain transmissibility as well as a comparison of the effectiveness of MEMS accelerometers with FBGs.

2.1.4. Defect Inclusion and Repair Procedure

Defects were induced, one at a time, at 10 cm intervals along the blade length from 10 cm to 110 cm away from the blade hub. The defects were in the form of 0.5 cm or 1.0 cm long simulated cracks created using a FEIN Multimaster saw [60]. All “cracks” were approximately 0.1 cm wide and were cut transverse to the blade on the leading edge of the blade’s inside curve. Cracks were chosen for this experiment given their prevalence in the wind industry as either an initial damage to the blade or the result of a different defect type that is allowed to propagate over time. Cracks are also easier to define and quantify than other defect types and are able to be simulated both numerically and experimentally. Figure 7c shows a typical defect.
After the blade’s vibration characteristics corresponding to each defect were measured, the defect was repaired using a MYNT3D printing pen [61] (MYNT3D, Salt Lake City, UT, USA) before the next defect was induced. The repair procedure involved roughening the defect area with sand paper. The region was then cleaned with alcohol to remove contaminants. Raise3D premium PLA filament was applied with the 3D printing pen at 230 °C to promote bonding with the base material. The resulting repair was then ground flush with the blade.

2.1.5. Experimental Testing Procedure

The blade was attached to an optical table using the 3D-printed baseplate and 4 screws. For each blade condition, the blade was struck with an impact hammer at a random location to provide excitation. The vibration response was measured by both sensor arrays using an Arduino Uno R3 [62] (Arduino, Monza, Italy) for the accelerometers and a combination of a laser and interrogator from Ibsen Photonics [63,64] (Ibsen Photonics, Farum, DK) for the FBGs. Separate impact tests were conducted for the accelerometers and the FBGs for each blade condition. The data for each test were filtered and converted to the frequency domain using a fast Fourier transform (FFT) in a Python 3.10.4 script. Five tests were completed for a baseline condition (undamaged blade) as well as for each defect condition. The first three natural frequencies were determined for each test. The natural frequencies selected for transmissibility analysis were based on the median value of the five trials. The transmissibilities of the five trials for each possible pair of sensors were averaged at each median natural frequency for the real, imaginary, amplitude and phase. The transmissibility change between the baseline and each defect condition was then determined for each sensor pair.

2.2. ANSYS Numerical Model

2.2.1. Model Set-Up

The CAD model used for 3D printing the blade and baseplate was added to ANSYS Workbench 2022 R1 [65] as a single blade structure. The sensors were modelled as distributed mass components on both sides of the blade. An overview of the simulated blade is shown in Figure 8.

2.2.2. Material

A custom isotropic linear elastic material was created in ANSYS for Raise3D Premium PLA. The Young’s modulus and density were determined using averages obtained experimentally from the 3D-printed adhered coupons discussed in Section 2.1.2. A Poisson’s ratio of 0.33 was used based on research in [66]. A summary of these material values is shown in Table 1.

2.2.3. Defect Inclusion

Transverse through-wall slots of 0.5 cm, which were 1.0 cm in length and 0.1 cm in width, were added to the model. Each defect was added in 10 cm intervals (in line with the sensors) at the leading edge on the upper side of the blade. Figure 8 shows a typical 1.0 cm defect.

2.2.4. Mesh

Tetrahedron elements that were 5 mm, 1 mm, and 10 mm were used for the body, sensor locations, and baseplate, respectively. Defects were primarily modelled using 1 mm mesh elements, but some were subsequently adjusted to achieve convergence. A typical blade mesh is shown in Figure 8.
To test whether these mesh sizes were appropriate, the modal analysis results for the undamaged ANSYS model were compared with the average experimental results of the pristine baseline blade. This comparison is shown in Table 2. All three simulated natural frequencies were within 1% of the experimental results, and the mesh geometry was, therefore, deemed adequate.

2.2.5. Numerical Simulations Using Modal and Harmonic Analyses

The numerical transmissibility simulations were performed using two features of the ANSYS workbench. First, modal analysis was used to determine the natural frequencies. The harmonic response function was then used to apply an input force of random magnitude and location to the simulated blade. The output response was calculated at frequencies from 0 to 50 Hz in increments of 0.05 Hz to capture the first three natural frequencies. A damping ratio of 1.2% was used as suggested by [67]. Transmissibility values were then calculated for each possible sensor pair using the real, imaginary, amplitude and phase responses at each natural frequency.

2.3. Total Method

Transmissibility research such as [47,49] primarily assessed transmissibility changes using pairs of sensors that were adjacent to each other. Comparing adjacent pairs was found to be effective with varying damage identification methods.
When a flaw is implemented, transmissibility changes also occur with other pairs of sensors within the sensor matrix. For example, if a defect occurs near sensor i, a change in transmissibility will be found in all sensor pairs associated with sensor i, not just between sensor i and its adjacent sensors. The Total Method in Equation (9) uses a variation of the Damage Indicator from [46] in Equation (8) and sums all the transmissibility changes from all the available pairs for each sensor at each natural frequency. This formula yields a value, Y i , for the total transmissibility change for a sensor in the sensor matrix, where m is the total number of sensors that can be paired with sensor i.
D i j ( ω n ) = T i j D a m a g e d ω n T i j U n d a m a g e d ω n T i j U n d a m a g e d ω n
Y i ω n = j = 1 m D i j ω n
With the Total Method, the transmissibility changes are amplified as data from multiple sensors are used. This results in a higher total transmissibility change than when using just the adjacent sensors, which makes it easier to identify a defect. The Total method was chosen for this experiment as it provides a higher sensitivity to small defect identification while suppressing any noise effects emanating from a single sensor. This was crucial for the defect sizes chosen for the evaluation of the transmissibility concept on the wind turbine blade.

3. Results

3.1. Strain Transmissibility

3.1.1. Variables

Non-destructive defect detection on the model turbine blade was assessed by means of strain transmissibility using combinations of variables from four different categories: defect size, vibration mode, strain orientation, and signal type. A total of 24 different variable combinations were used to analyze the transmissibility changes experimentally along a single axis with the FBGs. The ANSYS model was able to examine the strain in all three axes, so 72 variables were tested numerically. These variable combinations are listed in Table 3 and were all assessed to determine which best determined the existence and location of a flaw. Further analysis of the viability of the strain transmissibility concept was then pursued for this optimal choice of the four variables.

3.1.2. Defect Identification Threshold and Confusion Matrix Development

Qualitative analysis of the strain transmissibilities found that damage could be identified and located in select models by identifying the sensor with the largest total method value. An example of the experimental results is shown in Figure 9. In this figure, a 1.0 cm defect is added to a location 100 cm from the blade hub, and the total method values are calculated at each of the 11 sensors using Equations (8) and (9). The maximum value for this dataset is noted to occur at 100 cm—the location of the defect.
A damage threshold value is required to properly analyze the blade condition with different defect sizes in different locations. This value can be used to simplify the damage detection algorithm to a binary response function. If the total method value at a sensor is above the threshold, damage is assumed to be present at this location and a value of one is returned. If the total method value at a sensor is below the threshold, damage is assumed to be absent at that location and a value of zero is returned. An ideal threshold value should be small enough to identify defects in as many locations as possible. It must, however, also be a large enough value to minimize false positive identifications of defects at other sensor locations.
A confusion matrix analysis was used to determine an ideal threshold value for each set of variables listed in Table 3 and identify the combination that could most reliably detect damage with minimal occurrence of false positives. The confusion matrix determines the true positives (TP), true negatives (TN), false positives (FP) and false negatives (FN) when comparing a set of data to an ideal case [68]. Eleven defects were assessed for each possible combination of variables. The results from each variable combination were compared to an 11 × 11 matrix of zeros with ones on the diagonal (ones at each damage location in the ideal case and zeros in all other locations).
The “accuracy” and “precision” score metrics from Scikit-Learn [69] were used to identify the best overall model. The accuracy equation looks at the overall predictability. The precision gives an indication of the potential presence of false positives when compared to true positives. These equations are shown below.
A c c u r a c y = T P + T N T P + T N + F P + F N
P r e c i s i o n = T P T P + F P

3.1.3. Experimental Results

Experimental results for detecting blade damage will be examined first to characterize the effect of each of the four identified variables on the output of the confusion matrix.

0.5 cm Results

Figure 10 summarizes the results of the confusion matrix outputs for each set of parameters in Table 3 for the 0.5 cm defect. The best performing set of variables in this plot is the amplitude mode 1 set, which achieves accuracy scores above 0.9. The challenge with identifying the 0.5 cm defect is shown in the precision plot. In Figure 10b, Amplitude mode 1 also has the highest precision score with a maximum score of 0.67 occurring between threshold values of 10 and 11. This score indicates a high percentage of false positives occurring relative to the number of true positives the model can identify. Furthermore, the precision score decreases as the threshold value increases, which is an indication that the maximum total method values for the 0.5 cm defect are very small. True positives become false negatives as the damage threshold is increased above the maximum total method value. The precision score decreases as it is dependent on the number of true positives, as shown in Equation (11). Although the amplitude mode 1 set of parameters can identify a small range of defects (90–100 cm) at a threshold of 10–11, it is not recommended to use due to the low precision score. Overall, the 0.5 cm defect was determined to be too small to be accurately and precisely identified with any of the strain transmissibility methods.

1.0 cm Results

The larger-sized defect proved easier to identify, with higher accuracy and precision scores compared to the 0.5 cm defect. The 1.0 cm results are shown in Figure 11.
The maximum accuracy score in Figure 11 occurs at a threshold between 12 and 24 for the mode 1 amplitude model. This model performs the best out of all the available options. The precision score for this model increases to a value of 1 after a threshold value of 12, indicating that no false positives are occurring with this model.
The Total Method transmissibilities for all 11 defect locations for the mode 1 amplitude variable combination is shown in Figure 12 for the 1.0 cm defect. This plot depicts all the positive Total Method results. All negative sensor results were set to zero to aide in plot visualization. In this figure, maximum total method values are noted at a number of sensors corresponding to the damage location. This maximum value increases as the defect moves away from the hub and toward the tip of the blade. Closer to the blade hub, the peak values are less prominent, and the damage location is difficult to identify compared to nearby sensor locations. The implementation of a threshold removes these locations from the resulting defect identification range as the maximum total method value is too small to reliably identify the damage location at this defect size.
Figure 12 shows that defects can be successfully located at a range of 70 cm to 110 cm for the 1.0 cm defect utilizing a threshold value of 12. At any threshold value between 12 and 24, this location range can be identified. This model is both accurate and precise per the confusion matrix analysis and is, therefore, suitable for detecting defects in this range. This is a promising result, as many wind turbine defects occur close to the tip of the blade.

3.1.4. Numerical Results

The numerical model results developed with ANSYS are examined next. First, a comparison is made between the numerical model and the highest performing experimental model to assess the feasibility of the ANSYS model for defect detection. The numerical results of the other variable combinations are also reviewed.

Comparison with Experimental Results

Given that experimentally the best result was based on the mode 1 amplitude variable combination with sensors oriented in the Z direction, the ANSYS numerical model was first assessed under these conditions to determine if it achieved similar results. The numerical results showed that like the experimental model, a peak value occurred at the sensor closest to the damage. Significantly larger total method values occurred in the numerical model compared to the experimental results. This is shown in Figure 13 for the 1.0 cm defect size.
A confusion matrix analysis found that the numerical model required a higher threshold to eliminate the false positives versus the experimental model. This was likely due to the larger total method values produced overall by the numerical model at all possible sensor locations. A threshold value of 18 obtained the highest accuracy score of 0.975 and the highest precision score of 1.0. This is shown in Figure 14. The larger total method values in the numerical model also yielded an increased defect range. A 1.0 cm defect is identifiable in ANSYS at locations ranging from 40 cm to 110 cm from the blade hub; this range encompasses the successful experimental result range of 70 cm to 110 cm. The ANSYS numerical model is therefore a useful tool for planning experimental NDT studies to identify defects with a minimum size of 1.0 cm. Caution should be taken when assessing defects located in the 40 cm to 60 cm range numerically, as the expected experimental values will be much smaller and may not exceed the damage threshold for identification.

Other Variables

The amplitude mode 1 model performed the best overall experimentally and was also a high performing model numerically. However, the highest performing model numerically was the imaginary mode 1 variable combination as shown in Figure 14 for the 1.0 cm defect. The corresponding experimental results were poor for this set of parameters. Further improvements to the experimental methods and signal analysis may provide better results when using the imaginary portions of the signal in future research.
The numerical model was able to obtain results for the 0.5 cm defect at some locations, in part due to the overall higher total method values outputted by the numerical model. Like the experimental model, the threshold value with the highest accuracy and precision was only useful at identifying a few damage locations. Although more accurate and precise than the experimental model, the numerical model is also not recommended for predicting defects at a size of 0.5 cm or less. The reliable damage location range is not large enough for this model to be effective.
The X and Y strain orientations were also examined with the numerical model. In the confusion matrix analysis of the 1.0 cm defect, both directions performed similarly to the Z direction, with both the amplitude and imaginary responses in the first mode. In a direct comparison of the outputs of these models, the X-orientation model produced slightly higher total method values at the damaged sensors, indicating that this orientation is potentially more sensitive to this type of defect. Further experimental testing is recommended to assess the strain transmissibility results in the X direction.

3.2. Acceleration Transmissibility

3.2.1. Variables

Flaw detection via acceleration transmissibility was assessed using the same variable categories as those for strain transmissibility: defect size, vibration mode, acceleration orientation, and signal type. Table 4 presents the variables examined. A total of 72 different variable combinations were used in ANSYS numerical analysis. However, initial experimental results with MPU-6050 accelerometers showed that the acceleration in the Z direction had the least excitation; it was therefore eliminated from further data collection and analysis to reduce the sampling time. As a consequence of that decision, only 48 variable combinations were assessed experimentally for acceleration in the X and Y orientations. In addition, defect locations for the experimental tests ranged only from 10 cm to 80 cm away from the blade hub, as the maximum number of MEMS accelerometers was limited to eight with the available experimental testing equipment.

3.2.2. Experimental Results

Experimentally, neither of the 0.5 cm and 1.0 cm defect sizes were detectable at any location on the model blade using the variable combinations described in Table 4. The data acquired from the MPU-6050 sensors were random, and no patterns were noted using qualitative methods. The assessment of the ideal variable combination was carried out using ANSYS simulations, and the differences between the experimental and numerical results for acceleration are discussed in Section 4.3.

3.2.3. Numerical Results

The numerical results were examined qualitatively to determine if a pattern existed for locating defect damage using acceleration transmissibility and the total method as described in Section 2.3. Several different patterns emerged from the ANSYS numerical transmissibility results that could be used to indicate the presence and location of a defect. An example is shown in Figure 15, where the total method value for each sensor was calculated using Equations (8) and (9).
For the results shown in Figure 15, the transmissibility readings are positive for all the sensors located at a distance greater than 50 cm—the location of the damage. Negative transmissibility values are present for all sensors from 10 cm to 50 cm. For the acceleration data, the defect’s location can be identified by determining where the data intercept the X-axis. Most variable combinations with tangible results followed the pattern in Figure 15. However, some combinations showed an inverted pattern, where the sensors before the damaged sensor were positive, and negative values occurred after the sensor closest to the damage location. Both of these pattern types were assessed with a confusion matrix for all the variable combinations shown in Table 4. The results presented below are for the pattern shown in Figure 15, as this pattern was more prevalent and had the highest accuracy scores overall.
Based on the qualitative analysis, an initial damage threshold of zero was used to determine the damage location. All total method values that were below zero were set to zero, and all total method values above zero were set to one and compared to an ideal results matrix. In an ideal 8 × 8 matrix assuming perfect results, all sensors before and including the damaged sensor yield a value of zero (below the threshold). All sensor locations after the defect yield a value of one (above the threshold). The accuracy scores of the resulting confusion matrices were calculated using Scikit-Learn and Equation (10). The accuracy score was used as the deciding factor for the best overall model as the acceleration results rely on the outputs of multiple sensors to identify each damage location instead of a singular sensor output. The more sensors that match the ideal matrix outputs overall, the easier it will be to identify the damage location. The accuracy scores of all 72 variable combinations are shown in Figure 16.
Figure 16 shows that only a few select variable combinations have accuracy scores in excess of 0.7. This indicates that very few combinations can clearly identify the defect location at any given defect size. The majority of the variable combinations assessed produced random results. The best performing set of parameters was X amplitude mode 2, with a defect size of 1.0 cm. This variable combination had the highest accuracy score of 0.875, indicating that it would be the most reliable at identifying the location of a flaw numerically. The equivalent variable combination had an accuracy score of 0.75 with a 0.5 cm defect. This indicates that a relationship exists between defect size and accuracy for this set of parameters, where the larger the defect, the higher the accuracy score. The detailed results for the X amplitude mode 2 parameter combination are shown in Figure 17 for the 1.0 cm flaw size. Further analysis of the results for the detection of 0.5 cm flaws was not pursued due to the relatively poor success rate for locating defects of this size.
In Figure 17, the total method values change from negative to positive at the sensor closest to the damage when a 1.0 cm defect is added. Each intercept location is marked with a red star. The location of the defect can, therefore, be determined. A defect location range of 10–70 cm can be correctly identified with the X amplitude mode 2 results.
Although the location of the 1.0 cm defect can be determined through numerical modelling, an additional threshold was required to minimize the probabilities of “false positives” so that this model could be of practical use. The maximum Total Method values were examined for the X amplitude mode 2 parameter set, and a value of 1.4 was selected as the threshold to indicate if a defect of 1.0 cm or greater had occurred. This was the highest possible threshold that would capture the presence of a defect over the full location range of 10 cm to 70 cm. Given that the overall transmissibility changes due to a defect are very small and, therefore, more susceptible to noise, the threshold needs to be as high as possible to minimize the likelihood of false alarms.
Based on numerical modelling, overall cumulative transmissibility values are much smaller using acceleration versus strain. This makes the presence of a defect more difficult to assess with the acceleration model. The acceleration algorithm is also much more complex, with two thresholds required to identify that the blade has a defect, as well as determine its location. The primary advantage of utilizing acceleration transmissibility versus strain transmissibility is in locating the defect.
The strain transmissibility results were more sensitive to defect locations closer to the tip of the blade. As shown in Figure 17, the numerical acceleration transmissibility can be used to locate a defect situated in the range of 10–70 cm, with more prominent transmissibility changes occurring when the defect location is closer to the blade hub. The acceleration transmissibility results would be useful for defects in this range as they are not covered by the range of strain transmissibility. Acceleration transmissibility, in tandem with strain transmissibility, can increase the overall defect identification range on the blade. However, further research is required experimentally to properly validate the numerical acceleration results.

4. Discussion

4.1. Transmissibility Analysis vs. Natural Frequency Analysis

NDT, through vibration analysis, often relies on monitoring changes to the natural frequencies of the structure caused by damage. This method is not sensitive enough to detect the 1.0 cm blade damage considered in this study. For all defect sizes and locations tested, a maximum change of only 0.58% was noted at the first three natural frequencies from the experiments. The ANSYS simulations also indicated a very small change in the natural frequencies of 0.32% across the first three natural frequencies. These frequency changes are too minute to detect small defects. Much larger defects could significantly shift the natural frequencies and signal potential issues, albeit at a stage that would necessitate more extensive repairs.
The signal amplitudes can also be used in vibration monitoring to ensure that blade damage does not exceed a level that could compromise the blade’s structural integrity. However, for a system where the applied force magnitude and location are entirely random, it is difficult to appropriately compare the amplitudes of two different signal outputs to identify damage.
This study demonstrated that the transmissibility method may be able to locate a 1.0 cm defect at certain locations along the length of the blade and can be considered a worthwhile alternative to traditional vibration monitoring methods. The downside to transmissibility is the number of sensors required to perform this analysis. Regular vibration monitoring can be completed by one or two sensors placed in each direction to assess amplitude or natural frequency changes. Transmissibility analysis requires an array of sensors across the structure to properly identify and locate the defect.

4.2. Transmissibility Analysis at the Natural Frequencies

As discussed in Section 1.2.3, the transmissibilities at the natural frequencies can be used for defect detection regardless of the applied blade excitation force characteristics. The theory suggests that all vibration modes offer a similar potential to characterize a flaw and that it is possible to detect any size of a defect. In practice, certain vibration modes proved more effective than others in detecting defects, and a vibration sensor’s sensitivity and orientation also influenced its flaw detection ability. The results from this project found that a minimum defect size was required in order for the changes in transmissibility at the natural frequencies to be detected by the vibration sensors.
The ANSYS numerical simulation found that the first and third vibrational modes of the blade are primarily associated with displacements oriented in the flap-wise Y direction. The second mode primarily featured displacements along the edgewise X direction. When performing NDT through transmissibility analysis, it must be considered that the application of various random forces will not equally excite the blade in all directions.
Repeatable and accurate measurements of natural frequencies associated with the second and third vibrational modes were very difficult in a laboratory setting, likely due to a lack of strong excitation of these modes by the applied impact force. In practice, field testing with the given wind turbine geometry and weather patterns can determine which vibrational modes are adequately excited for acceleration and strain transmissibility analysis. In our set-up, mode 1 was clearly the best candidate as it was the only mode to yield useful damage information experimentally while using strain transmissibility.
Figure 18 shows the ANSYS model strain amplitude transmissibilities in the X, Y, and Z directions between sensors at 10 cm and 20 cm from the blade hub of an undamaged blade subjected to two different random forces.
In the first mode, the strain transmissibilities remain relatively constant around the natural frequency in all three directions, regardless of the magnitude and location of the applied forces. Transmissibility variations are slightly larger in the third vibration mode but still converge to the natural frequency in all orientations. This suggests that when using modes one and three, correct identification of the natural frequency is not critical. The X and Y strain transmissibility pairs do not converge at the second natural frequency, which demonstrates the importance of the sensor’s measurement direction. The X and Y directions fail to identify defects in the second mode as they do not follow the theory in Equation (7), and the transmissibility values depend on the applied force. The Z direction transmissibilities converge close to the second natural frequency but diverge at frequencies that are immediately adjacent. Therefore, to use the mode two transmissibility in the Z direction, one must accurately identify the natural frequency, as inherent transmissibility variations around the natural frequency can obscure those caused by damage. Experimentally, accurate identification of the second mode natural frequency was difficult; this explains the inability to detect defects using strain transmissibility with the second vibration mode.

4.3. Strain Transmissibility vs. Acceleration Transmissibility

The strain transmissibility results were easier to interpret than the acceleration results for the transverse defects tested experimentally and modelled numerically. A one-step algorithm could consistently locate the defect under all orientations, modes and signal types by identifying the sensor with the highest total method value. In contrast, the acceleration-based models required a two-step process: first to identify the presence of a defect and then to locate it. Different damage identification patterns emerged in the qualitative analysis of the acceleration numerical model depending on the selected variables for analysis. This added an additional challenge in creating an algorithm for the acceleration transmissibility analysis that would be suitable for all variables.
Experimentally, only the strain transmissibility was effective at locating defects. The acceleration transmissibility was unable to identify damage at any of the locations tested. The numerical models showed that acceleration transmissibility was viable for damage detection, but the overall total method values achieved were miniscule for both the 0.5 cm and 1.0 cm defects.
As shown in comparison between the numerical and experimental methods for strain transmissibility analysis, the total method values in the numerical model were significantly larger than those obtained through experiments. The same numerical model was used for the acceleration analysis. As strain and acceleration are related through a time derivative and Hooke’s Law, this suggests that the true acceleration results are perhaps much smaller than the numerical prediction. Experimental noise and artifacts may have obscured the defects. This indicates that the minimum detectable defect size is potentially larger than 1.0 cm for acceleration transmissibility analysis or that a more sensitive accelerometer is required. Further research is required, particularly with larger defects and alternative sensor types, to experimentally validate the acceleration transmissibility method.
The fibre optic technology in the case of strain sensing was also advantageous over the MEMS-based acceleration sensors as multiple FBGs were attached in series on a single fibre, and the final set-up was much lighter than its accelerometer counterpart. The sampling rate of the fibre optic interrogator could also support a much higher frequency range. The MPU-6050 accelerometers, while having the advantages of lower cost and being multi-directional, were limited by the maximum achievable sampling rate of the Python and Arduino operating codes. As a result, the number of accelerometers was limited to eight, and only two of three acceleration orientations could be sampled. These changes allowed the equipment to sample at a rate high enough to analyze the first three natural frequencies without being subject to aliasing.

4.4. Total Method Effectiveness

The Total method for processing vibration signals was found to be effective for locating transverse defects at sensor locations. It is feasible as a less complex transmissibility calculation than other methods found in research. This method utilized signals from all sensors to amplify the transmissibility changes at a defect location, making it easier to identify defects than using only adjacent sensor pairs. As each output is dependent on all the sensors within the network, this method can be susceptible to erroneous readings from a defective sensor, which can mask the defect.
A direct comparison between the Total Method and other transmissibility calculations was outside of the scope of this project. Further research is recommended comparing the Total method to other signal analysis methods found in the literature to determine the most effective transmissibility calculation technique for damage detection.

4.5. Comparison with Other Transmissibility Research

Natural frequency-based transmissibility originally discussed in [44,47] was validated in this research for damage detection. This experiment assessed the transmissibility changes at each exact natural frequency instead of examining changes at a frequency range around each natural frequency. The effects of sensor orientation, signal type, defect size, location and defect size were also examined in this research. The damage imparted to the test structure was also not large enough to shift any of the lower vibrational modes in this case.
It was argued in [49] that strain transmissibility is more sensitive to damage than acceleration transmissibility. The experiments in this paper aligned with this argument for both the numerical and experimental results. Strain transmissibility was more sensitive to the smaller defect size overall, but this could potentially be explained by the difference in sensors used experimentally. Fewer steps were also required to analyze the strain experimental data versus the acceleration data, making the strain transmissibility method more practical for use.
In [52], the focus of the transmissibility experiments was on the propagation of a defect in the middle of the blade resulting from a fatigue test. The research in this paper expanded on this by examining defects in different locations and by looking at multiple natural frequencies. This research also developed a threshold that could be used to assess that damage had occurred, which would enable this transmissibility method to be used for online monitoring.

5. Conclusions

This project was able to identify a transverse defect on a 3D-printed wind turbine blade using transmissibility analysis. This research found that the amplitude results were the best at identifying the defect using strain transmissibility at the first natural frequency. This variable combination was able to detect experimentally a 1.0 cm defect at a range of 70 cm to 110 cm away from the blade hub. This range was determined using a damage threshold value, which identifies whether a defect is present. A numerical model was completed and verified with the strain experimental model in the Z direction.
Acceleration transmissibility was unable to identify a defect experimentally. Numerically, the 1.0 cm defect was best identified with the amplitude results in the X direction at the second natural frequency. It was found that the numerical model required a two-step algorithm to both identify that a defect had occurred and determine the location using acceleration transmissibility.
Overall, this research shows that transmissibility is promising as an online method of damage detection for wind turbine blades as a small-sized defect was able to be identified at a natural frequency and thus be independent of the unknown applied forces to the wind turbine blade during operation. The transmissibility analysis used in this project found limitations in the size of the defect it could identify with the scaled blade. The effectiveness of transmissibility when examining other defect parameters is currently unknown. This research also found that the natural frequencies of the blade were not equally successful at identifying the defect. Careful natural frequency selection is required to use the transmissibility method. Future work is currently in progress testing additional crack sizes, orientations and locations on the blade to assess the effectiveness of transmissibility and further identify its advantages and limitations. Additional research is also required to test other defect types at different locations and orientations on the blade to assess their relationship with the natural frequencies prior to conducting tests on in-service blades.

Author Contributions

Conceptualization, R.H., A.S. and F.A.; methodology, R.H., A.S. and F.A.; software, R.H.; validation, R.H.; formal analysis, R.H.; investigation, R.H.; resources, F.A.; data curation, R.H. and F.A.; writing—original draft preparation, R.H.; writing—review and editing, A.S. and F.A; visualization, R.H. and F.A.; supervision, A.S. and F.A.; project administration, F.A.; funding acquisition, A.S. and F.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Canadian Natural Sciences and Engineering Research Council (NSERC), grant number RGPIN-2019-04096, and The University of Toronto.

Data Availability Statement

The data are available on request.

Acknowledgments

The authors would like to acknowledge the undergraduate and graduate researchers in the Decisionics Lab at the University of Toronto for their contributions and assistance with the execution of this experiment.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Kilpatrick, R.; Lafreniere, K.; Caceres, A. IEA Wind TCP Canada 2022. 2022. Available online: https://iea-wind.org/wp-content/uploads/2023/10/Canada_2022.pdf (accessed on 22 November 2023).
  2. Natural Resources Canada. Federal Government Invests in 12 New Wind, Solar and Smart-Grid Projects with Alberta Indigenous and Industry Partners. Calgary, 18 September 2023. Available online: https://www.canada.ca/en/natural-resources-canada/news/2023/09/federal-government-invests-in-12-new-wind-solar-and-smart-grid-projects-with-alberta-indigenous-and-industry-partners.html (accessed on 29 November 2023).
  3. SaskPower. Bekevar Wind and Power Line Project. Available online: https://www.saskpower.com/Our-Power-Future/Infrastructure-Projects/Construction-Projects/Planning-and-Construction-Projects/Bekevar-Wind-and-Power-Line-Project (accessed on 7 July 2024).
  4. Natural Resources Canada. Building Offshore Renewables in Newfoundland and Labrador and Nova Scotia. 30 May 2023. Available online: https://www.canada.ca/en/natural-resources-canada/news/2023/05/building-offshore-renewables-in-newfoundland-and-labrador-and-nova-scotia.html (accessed on 29 November 2023).
  5. Premier’s Office Natural Resources and Renewables. Province Sets Offshore Wind Target. 22 September 2022. Available online: https://novascotia.ca/news/release/?id=20220920003 (accessed on 29 November 2023).
  6. Vestas. V236-15.0 MW. Available online: https://www.vestas.com/en/energy-solutions/offshore-wind-turbines/V236-15MW (accessed on 7 July 2024).
  7. Vestas. Vestas Refurbished Turbines. Available online: https://www.vestas.com/en/products/refurbished-turbines (accessed on 29 November 2023).
  8. Siemens Gamesa Renewable Energy. A Clean Energy Solution—From Cradle to Grave Environmental Product Declaration SG 8.0-167 DD. Zamudio. Available online: https://www.siemensgamesa.com/en-int/-/media/siemensgamesa/downloads/en/sustainability/environment/siemens-gamesa-environmental-product-declaration-epd-sg-8-0-167.pdf (accessed on 12 September 2023).
  9. General Electric. Sierra Wind Platform. Available online: https://www.ge.com/renewableenergy/wind-onshore/sierra-platform/ (accessed on 29 November 2023).
  10. Enercon. Sustainability Report 2022. 2022. Available online: https://cdn.prod.website-files.com/64c38ca9b1a2e59bd5b7d64a/65c214deae43f54857116e27_ENERCON_Sustainability_Report_2022_low1.pdf (accessed on 7 July 2024).
  11. Natural Resources Canada. Canadian Wind Turbine Database. Available online: https://open.canada.ca/data/en/dataset/79fdad93-9025-49ad-ba16-c26d718cc070 (accessed on 10 September 2023).
  12. Bayar, T. Wind Turbine Blade Design Shift Up a Gear. Power Engineering International, 15 February 2015. Available online: https://www.powerengineeringint.com/renewables/strategic-development/sharpening-the-blade/ (accessed on 3 December 2023).
  13. Coultate, J. What’s the Next Big Thing in Wind Turbine Monitoring? Power Engineering, 3 June 2023. Available online: https://www.power-eng.com/om/whats-the-next-big-thing-in-wind-turbine-monitoring/#gref (accessed on 3 December 2023).
  14. Katsaprakakis, D.A.; Papadakis, N.; Ntintakis, I. A Comprehensive Analysis of Wind Turbine Blade Damage. Energies 2021, 14, 5974. [Google Scholar] [CrossRef]
  15. Mishnaevsky, L. Root Causes and Mechanisms of Failure of Wind Turbine Blades: Overview. Materials 2022, 15, 2959. [Google Scholar] [CrossRef]
  16. Wang, W.; Xue, Y.; He, C.; Zhao, Y. Review of the Typical Damage and Damage-Detection Methods of Large Wind Turbine Blades. Energies 2022, 15, 5672. [Google Scholar] [CrossRef]
  17. Ashwill, T.; Ogilvie, A.; Paquette, J. Blade Reliability Collaborative: Collection of Defect, Damage and Repair Data; Sandia National Lab.: Albuquerque, NM, USA; Livermore, CA, USA, 2013. [Google Scholar] [CrossRef]
  18. Hayman, B.; Wedel-Heinen, J.; Brøndsted, P. Materials Challenges in Present and Future Wind Energy. MRS Bull. 2008, 33, 343–353. [Google Scholar] [CrossRef]
  19. Katnam, K.B.; Comer, A.J.; Roy, D.; da Silva, L.F.M.; Young, T.M. Composite Repair in Wind Turbine Blades: An Overview. J. Adhes. 2015, 91, 113–139. [Google Scholar] [CrossRef]
  20. Vestas. Blade Services Expert Blade Service Solutions for All Your Turbines. 2019. Available online: https://nozebra.ipapercms.dk/Vestas/Communication/Productbrochure/PartsandRepair/blade-services1/ (accessed on 3 December 2023).
  21. Phase One. Drone-Captures High-Resolution Imagery for Wind Turbine Inspection with Phase One Cameras. Available online: https://geospatial.phaseone.com/case-studies/ati_wind-turbine-inspection-2/ (accessed on 9 December 2023).
  22. ABJ Renewables. Vestas Wind Turbine Inspection in Canada. Available online: https://abjdrones.com/vestas-wind-turbine-inspection-in-canada-with-drones/ (accessed on 9 December 2023).
  23. Reddy, A.; Indragandhi, V.; Ravi, L.; Subramaniyaswamy, V. Detection of Cracks and damage in wind turbine blades using artificial intelligence-based image analytics. Measurement 2019, 147, 106823. [Google Scholar] [CrossRef]
  24. Worzewski, T.; Krankenhagen, R.; Doroshtnasir, M.; Röllig, M.; Maierhofer, C.; Steinfurth, H. Thermographic inspection of a wind turbine rotor blade segment utilizing natural conditions as excitation source, Part I: Solar excitation for detecting deep structures in GFRP. Infrared Phys. Technol. 2016, 76, 756–766. [Google Scholar] [CrossRef]
  25. Worzewski, T.; Krankenhagen, R.; Doroshtnasir, M. Thermographic inspection of wind turbine rotor blade segment utilizing natural conditions as excitation source, Part II: The effect of climatic conditions on thermographic inspections—A long term outdoor experiment. Infrared Phys. Technol. 2016, 76, 767–776. [Google Scholar] [CrossRef]
  26. Roach, D.; Rice, T.; Paquette, J. Probability of Detection Study to Assess the Performance of Nondestructive Inspection Methods for Wind Turbine Blades; Sandia National Lab.: Albuquerque, NM, USA; Livermore, CA, USA, 2017. [Google Scholar] [CrossRef]
  27. Muñoz, C.Q.G.; Marquez, F.P.G.; Crespo, B.H.; Makaya, K. Structural health monitoring for delamination detection and location in wind turbine blades employing guided waves. Wind Energy 2019, 22, 698–711. [Google Scholar] [CrossRef]
  28. Ye, G.; Neal, B.; Boot, A.; Kappatos, V.; Selcuk, C.; Gan, T.-H. Development of an Ultrasonic NDT System for Automated In-Situ Inspection of Wind Turbine Blades. In Proceedings of the EWSHM—7th European Workshop on Structural Health Monitoring, Nantes, France, 8–11 July 2014; IFFSTTAR, Inria, Universite de Nantes: Nantes, France, 2014. Available online: https://inria.hal.science/hal-01020459/file/0156.pdf (accessed on 7 July 2024).
  29. Matsui, T.; Yamamoto, K.; Sumi, S.; Triruttanapiruk, N. Detection of Lightning Damage on Wind Turbine Blades Using the SCADA System. IEEE Trans. Power Deliv. 2021, 36, 777–784. [Google Scholar] [CrossRef]
  30. Van Dam, J.; Bond, L.J. Acoustic emission monitoring of wind turbine blades. In Proceedings of the SPIE Smart Structures and Materials + Nondestructive Evaluation and Health Monitoring, San Diego, CA, USA, 8–12 March 2015; Meyendorf, N.G., Ed.; [Google Scholar] [CrossRef]
  31. Tang, J.; Soua, S.; Mares, C.; Gan, T.-H. An experimental study of acoustic emission methodology for in service condition monitoring of wind turbine blades. Renew. Energy 2016, 99, 170–179. [Google Scholar] [CrossRef]
  32. Sensoria Edge-to-Edge Intelligence. Technology & Features. Available online: https://sensoria.com/technology/ (accessed on 10 December 2023).
  33. TWI Ltd. Structural Health Monitoring for Wind Turbine Blades. Available online: https://www.twi-global.com/media-and-events/insights/acoustic-emission-solution-for-structural-health-monitoring-of-wind-turbine-blades (accessed on 10 December 2023).
  34. Beale, C.; Willis, D.J.; Niezrecki, C.; Inalpolat, M. Passive acoustic damage detection of structural cavities using flow-induced acoustic excitations. Struct. Health Monit. 2020, 19, 751–764. [Google Scholar] [CrossRef]
  35. Weidmüller. Elevate Your Blade Performance. Available online: https://www.weidmueller.com/int/company/markets_industries/wind/blade_condition_monitoring.jsp#wm-1086239 (accessed on 11 December 2023).
  36. Di Lorenzo, E.; Petrone, G.; Manzato, S.; Peeters, B.; Desmet, W.; Marulo, F. Damage detection in wind turbine blades by using operational modal analysis. Struct. Health Monit. 2016, 15, 289–301. [Google Scholar] [CrossRef]
  37. Pacheco-Chérrez, J.; Probst, O. Vibration-based damage detection in a wind turbine blade through operational modal analysis under wind excitation. Mater. Today Proc. 2022, 56, 291–297. [Google Scholar] [CrossRef]
  38. Tcherniak, D.; Mølgaard, L.L. Active vibration-based structural health monitoring system for wind turbine blade: Demonstration on an operating Vestas V27 wind turbine. Struct. Health Monit. 2017, 16, 536–550. [Google Scholar] [CrossRef]
  39. Maia, N.M.M.; Almeida, R.A.B.; Urgueira, A.P.V.; Sampaio, R.P.C. Damage detection and quantification using transmissibility. Mech. Syst. Signal Process. 2011, 25, 2475–2483. [Google Scholar] [CrossRef]
  40. Maia, N.M.M.; Silva, J.M.M.; Ribeiro, A.M.R. The transmissibility concept in multi-degree-of-freedom systems. Mech. Syst. Signal Process. 2001, 15, 129–137. [Google Scholar] [CrossRef]
  41. Maia, N.M.M.; Urgueira, A.P.V.; Almei, R.A.B. Whys and Wherefores of Transmissibility. In Vibration Analysis and Control—New Trends and Developments; IntechOpen Ltd.: London, UK, 2011. [Google Scholar] [CrossRef]
  42. Yan, W.J.; Zhao, M.Y.; Sun, Q.; Ren, W.X. Transmissibility-based system identification for structural health Monitoring: Fundamentals, approaches, and applications. Mech. Syst. Signal Process. 2019, 117, 453–482. [Google Scholar] [CrossRef]
  43. Chesné, S.; Deraemaeker, A. Damage localization using transmissibility functions: A critical review. Mech. Syst. Signal Process. 2013, 38, 569–584. [Google Scholar] [CrossRef]
  44. Devriendt, C.; Guillaume, P. The use of transmissibility measurements in output-only modal analysis. Mech. Syst. Signal Process. 2007, 21, 2689–2696. [Google Scholar] [CrossRef]
  45. Johnson, T.J.; Adams, D.E. Transmissibility as a Differential Indicator of Structural Damage. J. Vib. Acoust. 2002, 124, 634–641. [Google Scholar] [CrossRef]
  46. Zhang, H.; Schulz, M.J.; Naser, A.; Ferguson, F.; Pai, P.F. Structural health monitoring using transmittance functions. Mech. Syst. Signal Process. 1999, 13, 765–787. [Google Scholar] [CrossRef]
  47. Devriendt, C.; Presezniak, F.; De Sitter, G.; Vanbrabant, K.; De Troyer, T.; Vanlanduit, S.; Guillaume, P. Structural Health Monitoring in Changing Operational Conditions Using Tranmissibility Measurements. Shock. Vib. 2010, 17, 153273. [Google Scholar] [CrossRef]
  48. Ghoshal, A.; Sundaresan, M.J.; Schulz, M.J.; Pai, P.F. Structural health monitoring techniques for wind turbine blades. J. Wind. Eng. Ind. Aerodyn. 2000, 85, 309–324. [Google Scholar] [CrossRef]
  49. Cheng, L.; Busca, G.; Roberto, P.; Vanali, M.; Cigada, A. Damage Detection Based on Strain Transmissibility for Beam Structure by Using Distributed Fiber Optics. In Structural Health Monitoring & Damage Detection; Niezrecki, C., Ed.; Springer International Publishing: Cham, Switzerland, 2017; Volume 7, pp. 27–40. [Google Scholar]
  50. Zou, Y.; Lu, X.; Yang, J.; Wang, T.; He, X. Structural Damage Identification Based on Transmissibility in Time Domain. Sensors 2022, 22, 393. [Google Scholar] [CrossRef] [PubMed]
  51. Fan, Z.; Feng, X.; Zhou, J. A novel transmissibility concept based on wavelet transform for structural damage detection. Smart Struct. Syst. 2013, 12, 291–308. [Google Scholar] [CrossRef]
  52. Yang, W.; Lang, Z.; Tian, W. Condition Monitoring and Damage Location of Wind Turbine Blades by Frequency Response Transmissibility Analysis. IEEE Trans. Ind. Electron. 2015, 62, 6558–6564. [Google Scholar] [CrossRef]
  53. Wang, X.; Zhang, L.; Heath, W.P. Wavelet Energy Transmissibility Analysis for Wind Turbine Blades Fault Detection. In Proceedings of the 2020 IEEE Applied Power Electronics Conference and Exposition (APEC), New Orleans, LA, USA, 15–19 March 2020; pp. 152–156. [Google Scholar] [CrossRef]
  54. Wang, X.; Zhang, L.; Heath, W.P. Wind turbine blades fault detection using system identification-based transmissibility analysis. Insight Non-Destr. Test. Cond. Monit. 2022, 64, 164–169. [Google Scholar] [CrossRef]
  55. Gaertner, E.; Rinker, J.; Sethuraman, L.; Zahle, F.; Anderson, B.; Barter, G.; Abbas, N.; Meng, F.; Bortolotti, P.; Skrzypinski, W.; et al. Definition of the IEA 15-Megawatt Offshore Reference Wind. Mar. 2020. Available online: https://www.nrel.gov/docs/fy20osti/75698.pdf (accessed on 14 May 2023).
  56. Raise3D. Raise3D Pro2 3D Printer. Available online: https://www.raise3d.com/products/pro2-3d-printer/ (accessed on 24 July 2023).
  57. Raise3D. Raise3D Premium PLA Technical Data Sheet. Version 5.1. Available online: https://s1.raise3d.com/2023/02/Raise3D-Premium_PLA_TDS_V5.1_EN.pdf (accessed on 10 July 2023).
  58. Glue, K. Products. Available online: https://www.krazyglue.com/products.html (accessed on 9 March 2024).
  59. InvenSense Inc. MPU-6000 and MPU-6050 Product Specification. PS-MPU-6000A-00. Available online: https://invensense.tdk.com/wp-content/uploads/2015/02/MPU-6000-Datasheet1.pdf (accessed on 12 August 2023).
  60. FEIN. FEIN Multimaster. Available online: https://fein.com/en_ca/machines/oscillating-multi-tools/multimaster/ (accessed on 25 August 2023).
  61. MYNT3D. MYNT3D 3d Pen PRO. Available online: https://www.mynt3d.com/products/3d-printing-pen (accessed on 25 August 2023).
  62. Arduino. Arduino Uno R3. Available online: https://docs.arduino.cc/hardware/uno-rev3 (accessed on 12 August 2023).
  63. Ibsen Photonics. I-MON 512 High Speed. Available online: https://ibsen.com/product/i-mon-512-high-speed/ (accessed on 30 July 2023).
  64. Ibsen Photonics. Ibsen Photonics-DL-BP1-1501A SLED. Available online: https://ibsen.com/product/light-source-for-fbg-sensing/ (accessed on 10 July 2023).
  65. ANSYS. Elevate Engineering Design 2022 Product Releases & Updates. Available online: https://www.ansys.com/en-gb/products/release-highlights-backup (accessed on 10 March 2024).
  66. Ferreira, R.T.L.; Amatte, I.C.; Dutra, T.A.; Bürger, D. Experimental characterization and micrography of 3D printed PLA and PLA reinforced with short carbon fibers. Compos. B Eng. 2017, 124, 88–100. [Google Scholar] [CrossRef]
  67. Medel, F.; Abad, J.; Esteban, V. Stiffness and damping behavior of 3D printed specimens. Polym. Test. 2022, 109, 107529. [Google Scholar] [CrossRef]
  68. Larner, A.J. The 2 × 2 Matrix; Springer International Publishing: Cham, Switzerland, 2021. [Google Scholar] [CrossRef]
  69. Scikit-Learn. Available online: https://scikit-learn.org/stable/index.html (accessed on 15 February 2024).
Figure 2. Wind turbine blade with a force applied at location k.
Figure 2. Wind turbine blade with a force applied at location k.
Sensors 24 04456 g002
Figure 3. Wind turbine blade with separate forces applied at locations k and l.
Figure 3. Wind turbine blade with separate forces applied at locations k and l.
Sensors 24 04456 g003
Figure 4. Front and side view of fully assembled wind turbine blade and baseplate.
Figure 4. Front and side view of fully assembled wind turbine blade and baseplate.
Sensors 24 04456 g004
Figure 5. Young’s moduli of (a) uniform full-size 3D printed coupons and (b) half-sized coupons glued together.
Figure 5. Young’s moduli of (a) uniform full-size 3D printed coupons and (b) half-sized coupons glued together.
Sensors 24 04456 g005
Figure 6. The layout of MPU-6050 accelerometers (A) and fiber Bragg gratings (FBGs) placed along the length of a wind turbine blade. The sensor number corresponds to its location on the blade in 10 cm intervals.
Figure 6. The layout of MPU-6050 accelerometers (A) and fiber Bragg gratings (FBGs) placed along the length of a wind turbine blade. The sensor number corresponds to its location on the blade in 10 cm intervals.
Sensors 24 04456 g006
Figure 7. (a) Typical MPU-6050 assembly on the blade. (b) Typical FBG assembly on the blade. (c) Typical 0.5 cm “crack” defect added to the blade. This defect was added 100 cm away from the blade hub.
Figure 7. (a) Typical MPU-6050 assembly on the blade. (b) Typical FBG assembly on the blade. (c) Typical 0.5 cm “crack” defect added to the blade. This defect was added 100 cm away from the blade hub.
Sensors 24 04456 g007
Figure 8. Wind turbine blade in ANSYS 2022 R1 with acceleration and strain sensor locations. A typical meshed 1.0 cm defect is shown in the top left corner. This defect is in line with a simulated sensor. The orientation axis is also shown on the bottom right of the figure.
Figure 8. Wind turbine blade in ANSYS 2022 R1 with acceleration and strain sensor locations. A typical meshed 1.0 cm defect is shown in the top left corner. This defect is in line with a simulated sensor. The orientation axis is also shown on the bottom right of the figure.
Sensors 24 04456 g008
Figure 9. Experimental (Z direction) mode 1 amplitude results for a 1.0 cm transverse defect. The defect is located at 100 cm.
Figure 9. Experimental (Z direction) mode 1 amplitude results for a 1.0 cm transverse defect. The defect is located at 100 cm.
Sensors 24 04456 g009
Figure 10. (a) Accuracy results and (b) precision results using the 0.5 cm defect size and the vibration mode and signal type variables from Table 3. All results are from the FBGs orientated in the Z direction.
Figure 10. (a) Accuracy results and (b) precision results using the 0.5 cm defect size and the vibration mode and signal type variables from Table 3. All results are from the FBGs orientated in the Z direction.
Sensors 24 04456 g010
Figure 11. (a) Accuracy results and (b) precision results using the 1.0 cm defect size and the vibration mode and signal type variables from Table 3. All results are from the FBGs orientated in the Z direction.
Figure 11. (a) Accuracy results and (b) precision results using the 1.0 cm defect size and the vibration mode and signal type variables from Table 3. All results are from the FBGs orientated in the Z direction.
Sensors 24 04456 g011
Figure 12. Mode 1 amplitude experimental results for the 1.0 cm defect at 11 locations 10 cm to 110 cm away from the hub. A damage threshold of 12 is added, and the identifiable defect range of 70 cm to 110 cm is shown on the wind turbine blade. Note that the data points are connected in these plots only to aide with visualization, and all negative sensor values are set to zero.
Figure 12. Mode 1 amplitude experimental results for the 1.0 cm defect at 11 locations 10 cm to 110 cm away from the hub. A damage threshold of 12 is added, and the identifiable defect range of 70 cm to 110 cm is shown on the wind turbine blade. Note that the data points are connected in these plots only to aide with visualization, and all negative sensor values are set to zero.
Sensors 24 04456 g012
Figure 13. Comparison of the Z direction ANSYS numerical model results with the experimental results for the 1.0 cm defect and the mode 1 amplitude parameters. Defects are located at 11 locations 10 cm to 110 cm away from the hub. The defect ranges of 40–110 cm numerically and 70–110 cm experimentally are shown on the wind turbine blade. Note that the data points are connected in these plots only to aide with visualization and all negative sensor values are set to zero.
Figure 13. Comparison of the Z direction ANSYS numerical model results with the experimental results for the 1.0 cm defect and the mode 1 amplitude parameters. Defects are located at 11 locations 10 cm to 110 cm away from the hub. The defect ranges of 40–110 cm numerically and 70–110 cm experimentally are shown on the wind turbine blade. Note that the data points are connected in these plots only to aide with visualization and all negative sensor values are set to zero.
Sensors 24 04456 g013
Figure 14. (a) Accuracy results and (b) precision results using the 1.0 cm defect size and the vibration mode and signal type variables from Table 3. All results are from the ANSYS numerical model with simulated sensors orientated in the Z direction.
Figure 14. (a) Accuracy results and (b) precision results using the 1.0 cm defect size and the vibration mode and signal type variables from Table 3. All results are from the ANSYS numerical model with simulated sensors orientated in the Z direction.
Sensors 24 04456 g014
Figure 15. Numerical mode 1 results (imaginary signal component) for a 1.0 cm transverse defect with sensor orientation in the Y direction. The defect is located at 50 cm from the blade hub.
Figure 15. Numerical mode 1 results (imaginary signal component) for a 1.0 cm transverse defect with sensor orientation in the Y direction. The defect is located at 50 cm from the blade hub.
Sensors 24 04456 g015
Figure 16. Accuracy score comparison of all 72 variable combinations for the ANSYS numerical model.
Figure 16. Accuracy score comparison of all 72 variable combinations for the ANSYS numerical model.
Sensors 24 04456 g016
Figure 17. X amplitude mode 2 results for the 1.0 cm defect for acceleration transmissibility. The horizontal line is the identification threshold where damage has occurred and is set to 1.4. The intercept locations are marked with red stars to identify damage locations. The identifiable defect range of 10 cm to 70 cm is shown on the wind turbine blade.
Figure 17. X amplitude mode 2 results for the 1.0 cm defect for acceleration transmissibility. The horizontal line is the identification threshold where damage has occurred and is set to 1.4. The intercept locations are marked with red stars to identify damage locations. The identifiable defect range of 10 cm to 70 cm is shown on the wind turbine blade.
Sensors 24 04456 g017
Figure 18. ANSYS strain amplitude transmissibilities calculated between the 10 cm and 20 cm sensors with two different random forces applied and plotted with respect to frequency. Each plot contains the X, Y and Z directional transmissibilities at the (a) first, (b) second and (c) third natural frequencies. Each natural frequency is shown as a dotted vertical line.
Figure 18. ANSYS strain amplitude transmissibilities calculated between the 10 cm and 20 cm sensors with two different random forces applied and plotted with respect to frequency. Each plot contains the X, Y and Z directional transmissibilities at the (a) first, (b) second and (c) third natural frequencies. Each natural frequency is shown as a dotted vertical line.
Sensors 24 04456 g018
Table 1. PLA parameters used in ANSYS.
Table 1. PLA parameters used in ANSYS.
Young’s Modulus (MPa)Density (kg/m3)Poisson’s Ratio
298411280.33 [66]
Table 2. Natural Frequencies for the baseline undamaged ANSYS simulation and experiments.
Table 2. Natural Frequencies for the baseline undamaged ANSYS simulation and experiments.
Vibration ModeANSYS Natural Frequency (Hz)Experimental Natural Frequency (Hz)
18.9988.990
222.23422.133
327.12227.150
Table 3. Strain Transmissibility Variables.
Table 3. Strain Transmissibility Variables.
Variable TypeFBG Experimental
Variables Examined
ANSYS Numerical
Variables Examined
Defect Size0.5 cm, 1.0 cm0.5 cm, 1.0 cm
Vibration Mode Number1, 2, 31, 2, 3
Strain OrientationZX, Y, Z
Signal TypeReal, Imaginary, Amplitude, PhaseReal, Imaginary, Amplitude, Phase
Table 4. Acceleration Transmissibility Variables.
Table 4. Acceleration Transmissibility Variables.
Variable TypeMPU-6050 Experimental
Variables Examined
ANSYS Numerical
Variables Examined
Defect Size0.5 cm, 1.0 cm0.5 cm, 1.0 cm
Vibration Mode Number1, 2, 31, 2, 3
Strain OrientationX, YX, Y, Z
Signal TypeReal, Imaginary, Amplitude, PhaseReal, Imaginary, Amplitude, Phase
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Henderson, R.; Azhari, F.; Sinclair, A. Natural Frequency Transmissibility for Detection of Cracks in Horizontal Axis Wind Turbine Blades. Sensors 2024, 24, 4456. https://doi.org/10.3390/s24144456

AMA Style

Henderson R, Azhari F, Sinclair A. Natural Frequency Transmissibility for Detection of Cracks in Horizontal Axis Wind Turbine Blades. Sensors. 2024; 24(14):4456. https://doi.org/10.3390/s24144456

Chicago/Turabian Style

Henderson, Rachel, Fae Azhari, and Anthony Sinclair. 2024. "Natural Frequency Transmissibility for Detection of Cracks in Horizontal Axis Wind Turbine Blades" Sensors 24, no. 14: 4456. https://doi.org/10.3390/s24144456

APA Style

Henderson, R., Azhari, F., & Sinclair, A. (2024). Natural Frequency Transmissibility for Detection of Cracks in Horizontal Axis Wind Turbine Blades. Sensors, 24(14), 4456. https://doi.org/10.3390/s24144456

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop