Noise Cancellation of Helicopter Blade Deformations Measurement by Fiber Bragg Gratings

The work presents data treatment methods aimed at eliminating the noise in the strain sensor data induced by vibrations of the helicopter blade in flight conditions. The methods can be applied in order to enhance the metrological performance of the helicopter weight estimation system based on the deformation measurement of the main rotor blades. The experimental setup included a composite plate fixed to the vibrating stand on the one end, with six fiber-optic strain sensors attached to its surface. In this work, the procedure of the optimal linear smoothing (POLS) and 3D-invariant methods were used to obtain monotone calibration curves for each detector, thereby making it possible to distinguish the increase of load applied to the free end of the plate with an increment of 10 g. The second method associated with 3D invariants took into account 13 quantitative parameters defined as the combination of different moments and their intercorrelations up to the fourth-order inclusive. These 13 parameters allowed the calculation of the 3D surface that can serve as a specific fingerprint, differentiating one set of initial data from another one. The combination of the two data treatment methods used in this work can be applied successfully in a wide variety of applications.


Introduction
Helicopter main rotor blades are essential components that define the flight performance and safety of the helicopter. Therefore, monitoring of their state is of the utmost importance during helicopter exploitation. Moreover, a sophisticated aeroelastic analysis of the rotor blades is necessary at the early stages of helicopter development, since the blades are subjected to significant displacements during flight and the aerodynamic conditions around the rotor are highly unsteady and complex [1]. In order to perform the abovementioned tasks, real-time blade deformation measurements are required, which enable the blade wear or damage detection during helicopter operation as well as the evaluation of blade displacements during development testing procedures.
The existing approaches to helicopter blade deformation measurement can be divided into four main categories: the systems based on the detection of reflected light [2][3][4], the systems based on laser units [5], the systems based on electrical resistance strain gauges [6][7][8], and the fiber-optic strain measurement systems [9][10][11].
The methods of the first group capture the light reflected from the rotor blades in order to determine their deformation. Thus, the system presented in [2] uses light from the external source that is modulated by the rotating blades to define the amplitude of blade flapping, its velocity, bending, etc. The system described in [3] comprises CCD detectors installed in the helicopter fuselage and is able to determine the blade pitch angle and flapping. The work [4] presents an approach to blade deformation monitoring based on the digital image correlation (DIC) technique, in which the surface of the blade is painted with a high-contrast random dot pattern. Two high-resolution digital cameras are separated by a certain distance and view the same area on the blade so that the parallax is used to measure three-dimensional displacements of the blade surface. In general, the main advantage of the systems based on the reflected light detection is that they do not require the structural modification of the blade. On the other hand, the performance of such systems depends on the conditions of the ambient environment, such as weather conditions, altitude, time of day, disturbances from external light sources.
In order to increase the measurement accuracy of blade deformation, in comparison with the aforesaid methods, systems based on laser units have been proposed. Thus, in [5], Projection Moiré Interferometry (PMI) with infrared laser source is applied to measure dynamic rotor blade deflections. Similarly, to the previous group of approaches, the laser units-based methods are realized without blade structural modifications, while not being affected by disturbances from external visible light sources. However, the accuracy of such approaches also depends on the weather conditions and optical density of the environment.
Electrical resistance strain gauges have been applied for the development of rotor blade strain measurement systems for more than four decades [6]. The NASA/Army UH-60A Airloads Program was completed in 1994 [7,8] and performed over 200 flights, collecting data from the instrumentation system installed on the main rotor of the helicopter. The system included 25 strain gauges, accelerometers, pressure transducers, temperature sensors, and acquired a vast amount of data, including deformation to estimate the airloads. Unlike the methods that detect optical radiation reflected from the blades, the systems based on electrical resistance strain gauges are immune to the influence of weather conditions and generally provide higher accuracy and resolution of measurements. However, the implementation of electrical resistance strain gauges requires the corresponding modifications of the blade structure. In addition, despite the low weight and compact dimensions of the gauges, the installation of a high number of such sensors may lead to a significant increase in the overall blade weight due to the metal wires connecting the gauges.
At present, a considerable amount of research activity is aimed at the development of fiber-optic sensor systems for blade deformation measurement. The main advantages of FBG-based techniques are low weight, immunity to electromagnetic interference, small dimensions (the optical fiber diameter is 125 µm), and the absence of an electrical power supply to the sensors [9]. Moreover, optical fibers can be embedded into composite structures during fabrication [10] or attached to the surface of the blade [11]. The experimental setup presented in [12] includes 60 fiber-optic sensors based on fiber Bragg gratings (FBG) nonuniformly allocated along four optical fibers so that one-third of the sensors are positioned close to the blade root where the strain gradients are higher. The system is used to measure blade deformation in vertical and lateral directions. The system discussed in [13] uses 16 FBGs interrogated by optical frequency domain reflectometry (OFDR) and is capable of retrieving the first and second bending modes of the blade.
The application of blade deformation measurements to aeroelastic analysis requires the detection of blade oscillatory deformations during aerodynamic testing procedures. However, an alternative application for such systems can be found in the vertical displacement definition of the blades in order to estimate the gross vehicle weight of a helicopter. The development of such systems significantly contributes to the safety of helicopters, since, at present, the weight is generally determined by manually tracking the weight of cargo, passengers, fuel, etc., which does not guarantee adequate accuracy. The work [14] proposes a system that calculates the weight of a helicopter by measuring and averaging the distances between the fuselage and the underside of each blade of the main rotor while in hover or level unaccelerated forward flight. The system comprises two laser distance measuring devices at the front and rear of the fuselage. Predetermined calibration data are used to compensate for the effects of vertical and horizontal drag forces on the fuselage due to the spinning rotor and forward velocity. The task of blade vertical displacement measurement can be also accomplished using strain sensors installed on the blades [14], similarly to the previously mentioned developments. In this application, the oscillatory behavior of the blade vertical displacements caused by aeroelastic effects may deteriorate the accuracy and resolution of the weight measurement. The problem is illustrated in Figure 1, in which Figure 1a represents a qualitative graph of a signal from the strain sensor located on a steady blade under the load that increased over time in five fixed steps, while Figure 1b shows the qualitative graph of the signal from the same sensor on the oscillating blade under the same sequence of loads. As it can be seen from this example, the oscillations significantly complicate the task of distinguishing the loads acting on the blade.
Sensors 2021, 21, x FOR PEER REVIEW cargo, passengers, fuel, etc., which does not guarantee adequate accuracy. The wo proposes a system that calculates the weight of a helicopter by measuring and ave the distances between the fuselage and the underside of each blade of the main rotor in hover or level unaccelerated forward flight. The system comprises two laser di measuring devices at the front and rear of the fuselage. Predetermined calibratio are used to compensate for the effects of vertical and horizontal drag forces on the fu due to the spinning rotor and forward velocity. The task of blade vertical displac measurement can be also accomplished using strain sensors installed on the blade similarly to the previously mentioned developments. In this application, the osci behavior of the blade vertical displacements caused by aeroelastic effects may dete the accuracy and resolution of the weight measurement. The problem is illustrated ure 1, in which Figure 1a represents a qualitative graph of a signal from the strain located on a steady blade under the load that increased over time in five fixed steps Figure 1b shows the qualitative graph of the signal from the same sensor on the osci blade under the same sequence of loads. As it can be seen from this example, the o tions significantly complicate the task of distinguishing the loads acting on the blad The aim of the current work is to propose the method for the elimination of n the strain measurement data caused by the blade oscillations, thereby enhancing th rological performance of the helicopter weight measurement system.

Deformation Measurement Using Fiber Bragg Gratings
A fiber Bragg grating (FBG) is a segment of an optical fiber that reflects par wavelengths of optical radiation and transmits the others. FBGs are manufactured scribing a systematic variation of the refractive index of the fiber core. The central length of the reflected light is called the Bragg wavelength and is defined according following generally accepted formulation: where neff is the effective refractive index of the grating, and Λ is the grating perio the period of the refractive index variation. The working principle of FBG is based on the change of either the effective refr index (neff) or the grating period (Λ) and the resulting shift of its central waveleng when the sensor is subjected to physical fields, such as strain, temperature, pressu midity, etc. [15].  The aim of the current work is to propose the method for the elimination of noise in the strain measurement data caused by the blade oscillations, thereby enhancing the metrological performance of the helicopter weight measurement system.

Deformation Measurement Using Fiber Bragg Gratings
A fiber Bragg grating (FBG) is a segment of an optical fiber that reflects particular wavelengths of optical radiation and transmits the others. FBGs are manufactured by inscribing a systematic variation of the refractive index of the fiber core. The central wavelength of the reflected light is called the Bragg wavelength and is defined according to the following generally accepted formulation: where n eff is the effective refractive index of the grating, and Λ is the grating period, i.e., the period of the refractive index variation. The working principle of FBG is based on the change of either the effective refractive index (n eff ) or the grating period (Λ) and the resulting shift of its central wavelength (λ B ) when the sensor is subjected to physical fields, such as strain, temperature, pressure, humidity, etc. [15]. Figure 2a shows the schematic representation of an FBG, and Figure 2b presents the typical spectral response of an FBG for transmitted and reflected light. The most common techniques for FBG multiplexing and interrogation, such as wavelength [16], time [17], frequency [18], polarization [19], and spatial [20] division are implemented using complex electro-optical devices, including spectrum analyzers, tunable Fabry-Perot interferometers, diffraction gratings, etc. The common disadvantage of the mentioned interrogation techniques is the complexity and high cost of devices required for their implementation. In order to simplify the FBG interrogation, the usage of addressed fiber Bragg structures (AFBSs) has been proposed [21,22]. An AFBS is a type of FBG, the spectral response of which has two ultra-narrowband components. The frequency spacing between the components is called the address frequency of the AFBS, it is in the microwave range and unique for each sensor in the system, and does not change when the AFBS is subjected to strain or temperature variations. The microwave-photonic approach for AFBS interrogation utilizes a broadband light source, an optical filter with predefined linear slope spectral response, and a photodetector. The central wavelength shift of AFBS is defined using the amplitude of the beating signal generated at the photodetector at the address frequencies of the sensors. This approach is further developed with the introduction of multiaddressed fiber Bragg structures (MAFBSs) [23], which have three or more ultra-narrowband components in their spectral response. The combination of address frequencies expands the sensor capacity of the measurement system and increases the accuracy of the MAFBS central wavelength definition.
The FBG interrogation device used in the current work implemented conventional wavelength-division multiplexing of FBG sensors, where each grating had a different Bragg wavelength so that their spectral responses did not superimpose on each other during measurements. The device was based on the Ibsen ® I-MON spectrum analyzer and was operated by custom software.

Description of the POLS Method
In order to decrease the noise evoked by vibro-oscillations, we use the procedure of the optimal linear smoothing (POLS) that was successfully applied earlier in papers [24][25][26][27][28][29] covering different regions of physics.
The basic formula is given by the following expression:  The most common techniques for FBG multiplexing and interrogation, such as wavelength [16], time [17], frequency [18], polarization [19], and spatial [20] division are implemented using complex electro-optical devices, including spectrum analyzers, tunable Fabry-Perot interferometers, diffraction gratings, etc. The common disadvantage of the mentioned interrogation techniques is the complexity and high cost of devices required for their implementation. In order to simplify the FBG interrogation, the usage of addressed fiber Bragg structures (AFBSs) has been proposed [21,22]. An AFBS is a type of FBG, the spectral response of which has two ultra-narrowband components. The frequency spacing between the components is called the address frequency of the AFBS, it is in the microwave range and unique for each sensor in the system, and does not change when the AFBS is subjected to strain or temperature variations. The microwave-photonic approach for AFBS interrogation utilizes a broadband light source, an optical filter with predefined linear slope spectral response, and a photodetector. The central wavelength shift of AFBS is defined using the amplitude of the beating signal generated at the photodetector at the address frequencies of the sensors. This approach is further developed with the introduction of multiaddressed fiber Bragg structures (MAFBSs) [23], which have three or more ultra-narrowband components in their spectral response. The combination of address frequencies expands the sensor capacity of the measurement system and increases the accuracy of the MAFBS central wavelength definition.
The FBG interrogation device used in the current work implemented conventional wavelength-division multiplexing of FBG sensors, where each grating had a different Bragg wavelength so that their spectral responses did not superimpose on each other during measurements. The device was based on the Ibsen ® I-MON spectrum analyzer and was operated by custom software.

Description of the POLS Method
In order to decrease the noise evoked by vibro-oscillations, we use the procedure of the optimal linear smoothing (POLS) that was successfully applied earlier in papers [24][25][26][27][28][29] covering different regions of physics.
The basic formula is given by the following expression: Here, K(t) is a smoothing kernel given in the form of the Gaussian function. The value w determines the width of the smoothing window. We want to remark the following attractive features of this expression: (a) This expression is linear with respect to the function y i subjected to the smoothing procedure. Therefore, it does not contain any additional distortions evoked by the possible treatment of a function y i ; (b) If w >> 1, then (as it is easily seen from Expression (2)) the smoothed function Ysm j coincides with its mean value; (c) If w becomes close to the zero value, then the smoothing kernel K(t) coincides practically with delta function δ(x j − x i ) and therefore, in this case, Ysm j (0) ∼ = y j . These attractive features allow the selection of high-frequency oscillations and leave only low-frequency "trend" that is proved to be useful for practical applications.
In this paper, we apply the simple Expression (2) for the elimination of high-frequency oscillations, and the remaining trend allows us to select the desired parameters that differentiate the oscillations without load from oscillations with a load.

Description of 3D-DGI Method
In this section, we describe the mathematical details associated with the derivation of the complete discrete geometrical invariant (DGI) in 3D space. We remind here that preliminary results based on the application of the incomplete DGI form of the fourth order in 3D space are outlined recently in [29]. Let us consider the complete form of the fourth-order power law.
In Expression (3), the upper indices define the combination of the variables y α (α = 1, 2, 3) fixing the location of an arbitrary point M(y 1 , y 2 , y 3 ) in 3D space, the low indices determine the values of the power law exponents that correspond to the algebraic form of the fourth order. The choice of the sign's combination (±) before the constants in (3) will be explained below. Three random sequences are determined by the values r αk (α = 1, 2, 3; k = 1, 2, . . . , N). Expression (1) represents itself the complete form of the fourth order that contains the combination of three variables associated with an arbitrary point M(y 1 , y 2 , y 3 ) and three arbitrary sequences r αk . The desired DGI is obtained from the following requirement: In order to remove in Expression (4) the cubic terms, we introduce the variables and nullify the linear terms. This requirement helps us to separate the desired variables Y α from each other and keep only the terms of the second and fourth orders, correspondingly.
In order to decrease the number of constants in (3) and derive the DGI not depending on some additional constants, one defines three key ratio constants R (α,β) , with (α,β) = (1,2), (1,3), (2,3), It is convenient also to introduce the following notations for the integer moments and their intercorrelations and present them as In the result of the introduced notations (6) and (7), the system of linear equations for the finding of unknown ratios R (α,β) from the nullification requirement of the entering linear terms accepts the following form: The linear system of equations helps to reduce three moments (Q 333 , Q 222 , Q 111 ) and seven intercorrelations of the third order (Q 332 , Q 322 , Q 221 , Q 211 , Q 331 , Q 311 , Q 321 ) to the calculation of three unknown ratios R (α,β) only. We should notice also that the combination of the algebraic signs in (3) is chosen in that way for the keeping of the partial solution R = 1 of system (8) in the case when all three random sequences r αk are identical to each other, i.e., r 1k = r 2k = r 3k . It is natural to define it as the case of spherical symmetry. If only two sequences coincide with each other (for example, r 1k = r 2k = r 3k ), then we deal with the case of the cylindrical symmetry. In this case, the linear system (8) is reduced to a couple of linear equations relatively the variables R (1,2) = R (1,3) = R (2,3) . After averaging procedure applied to Expression (4), the structure of the fourth-order form can be rewritten as The fourth-and the second-order forms entering to the left-hand side can be presented as The constants A αβ figuring in Expression (10b) are defined as The constant I 4 (defined by 3 moments and 12 intercorrelations of the fourth order) figuring in the right-hand side of (9) is defined as It is interesting to notice that in the case of the spherical symmetry (r 1k = r 2k = r 3k ), all correlations coincide with each other and the value of I 4 equals zero. The form of the fourth order (9) admits the separation of the variables in the spherical system of coordinates. If one accepts the following conventional notations: then substitution of these variables into (9) leads to the following biquadratic equation relatively the unknown radius R(θ,ϕ): The desired solution (R(θ, ϕ) > 0) is written as The polynomials P 2,4 (θ,ϕ) entering in (14) are defined by the following expressions: P 4 (θ, ϕ) = sin 4 θ · cos 4 ϕ + sin 4 θ · sin 4 ϕ + cos 4 θ+ +R (1,2) sin 2 θ sin ϕ cos ϕ cos 2 θ − sin 2 θ 2 (sin ϕ + cos ϕ) 2 + +R (1,3) sin θ cos θ cos ϕ sin 2 θ sin 2 ϕ − 1 2 (sin θ cos ϕ + cos θ) 2 + +R (2,3) sin θ cos θ sin ϕ sin 2 θ cos 2 ϕ − 1 2 (sin θ sin ϕ + cos θ) 2 (15a) P 2 (θ, ϕ) = A 11 sin 2 (θ) cos 2 (ϕ) + A 22 sin 2 (θ) sin 2 (ϕ)+ +A 33 cos 2 (θ) + A 12 sin 2 (θ) sin(ϕ) cos(ϕ)+ +A 13 sin(θ) cos(θ) cos(ϕ) + A 23 sin(θ) cos(θ) sin(ϕ).
The last Expressions (13)-(15) determine the final form of the DGI in 3D space. It includes three surfaces determined by Expression (13). The further analysis shows that Expression (14b) equals zero (because I 4 = 0) in the case of the coincidence of three compared random sequences (r 1k = r 2k = r 3k ). The radius R(θ,ϕ) can contain the complex expression when the integrand in (14b) becomes negative. It accepts the negative values when the constant I 4 (which in most cases is defined by Expression (12)) becomes negative. In this case, it is convenient to rewrite expression (13) in the form y 1 = y 1 + |R(θ, ϕ)| sin θ cos ϕ, y 2 = y 2 + |R(θ, ϕ)| sin θ sin ϕ, y 3 = y 3 + |R(θ, ϕ)| cos θ, |R(θ, ϕ)| = [Re(R(θ, ϕ))] 2 + [Im(R(θ, ϕ))] 2 , 0 ≤ θ < π, 0 ≤ ϕ < 2π. The curves defined by Equation (16) facilitate considerably the numerical analysis of initial data. As it follows from this preliminary analysis, this 3D surface is determined by the combination of 13 parameters: three moments of the first order <y α >, α = 1, 2, 3 from (3), six correlators of the second-order A αβ from (11), three reduced correlators R (α,β) from (8), and invariant of the fourth-order I 4 from (12). We want to stress here again that the final Expressions (14) and (15) do not use any model and are determined completely by the measured data together with their measurement errors. Finishing this section, one can say that this method can be applied for the reduction of initial data. It is necessary to note that the dimension of the radius [R(θ,ϕ)] coincides with the dimension of initial data y α .
This reduction procedure can be divided into the following stages: This expression, in spite of its simplicity, is really effective for the comparison of the statistical closeness of a pair of random functions belonging to the given/another sampling participating in the comparison operation. Really, if the function Q 1,2 (s) is located in the interval [1,2], then the pair random functions are statistically close to each other. In the case when Q 1,2 (s) ∈ [0,1), one can conclude that the pair random functions compared are statistically different. Besides this important parameter (16), one can take into account the symmetry of the random function y 1 (m). Any random function located in the rectangle M × Range[y(m)] crosses the line <Pr(m)>, coinciding with its mean value. Therefore, for evaluation of the symmetry of a random function, one can introduce the value If the value Sm(y) is located near zero (Sm(y) ≈ 0), then the line <y> divides the rectangle M × Range[y(m)] into two almost equal parts. In other cases, the value Sm(y) ∈ [−0.5, 0.5] determines the measure of asymmetry. After the application of Expression (17) for comparing similar columns (belonging to the same parameter Pr s ), one can receive finally the vector of the length S = 13 that contains information about the statistical closeness of two matrices compared. It is interesting to notice that simple Expression (17) can be used also for comparison each successive measurement with another one in the given rectangle matrix [N × M]. If one compares the vectors y m forming the columns of the initial matrix with each other, then in the result of application (17), one can obtain the symmetrical matrix U(m 1 ,m 2 ) (m 1,2 = 1, 2, . . . , M) with elements located in the interval 0 ≤ U(m 1 ,m 2 ) ≤ 2. Only elements located in interval 1 ≤ U(m 1 ,m 2 ) ≤ 2 will correspond to a "good" experiment, while the elements from interval 0 ≤ U(m 1 ,m 2 ) < 1 should be considered as possible "outliers" and correspond to "bad" experiment.
This normalization procedure makes all parameters Pr s statistically close to each other with mean value equaled zero mean(Prn s ) = 0 and with the Range(Prn s ) = 1. If one integrates Expression (19) for each parameter one can receive the statistically different curves JP s = Integral(Prn s ) for each initial parameter (s = 1, 2, . . . , S). The distributions of the ranges of these integral curves P 1 = Range(JP s ), together with the distribution of asymmetries P 2 = Range(Sm(JP s )) calculated with the help of Expression (18), give finally the matrix containing [S × 2] containing (S = 13) rows and two columns only. If we calculate the ranges of these two columns, we obtain finally two values only that can characterize the initial matrix [N × M]. If we have a set of matrices [N × M] q (q = 1, 2, . . . , Q), then this simple and general procedure allows us to select the "best" one having minimal values of these two key parameters P 1,2 . It will characterize the stability of the initial sequence and their minimal values will serve as a criterion for the selection of the "best" TLS among other TLS(s). We want to emphasize here that this final stage of treatment of "big" matrices differs from the procedure used in the paper [29]. Earlier, one of us (RRN) had the set of rectangle matrices that can be characterized as "normal/reference" ones and compared them with "strange/tested" matrices associated with defects. Data that will be analyzed below do not contain this information. Therefore, we propose another procedure described above for the selection of the "best" data expressed in the form of rectangle matrices.
Concluding this section, it is necessary to mention the following. The POLS is the universal tool that can be applicable for smoothing all available data. It possesses some remarkable features that make it really universal for many applications, which include the following:

1.
The POLS is a linear tool, and it does not distort initial data; 2.
When the value of the smoothing window tends to zero (w → 0), then the smoothed replica coincides with the initial data. In another limiting case, when w >> 1, the smoothed replica coincides with its arithmetic mean.
Another new instrument as the 3D geometrical invariants can be defined as the "universal" tool as well, by virtue of the following features: 1.
Thanks to 13 universal parameters defining the feature space, it allows the comparison of the different random sequences having different natures; 2.
It can be applied to analyses of the TLS and, therefore, this tool forms a universal platform that cannot contain treatment errors; 3.
Section 3, given above, gives an example for its application to real data.
We should also note that both methods keep the units of the initial data.

Experimental Setup
In order to imitate the process of helicopter blade deformation measurement, an experimental setup was designed that included a composite plate with six strain detectors based on fiber Bragg gratings attached to the upper surface of the plate using epoxy adhesive. The FBGs used in the experiment were manufactured at the Department of Radiophotonics and Microwave Technologies of KNRTU-KAI (Kazan) using a common fabrication technique based on the Lloyd interferometer. Figure 3a  During the experimental procedure, the plate was fixed on one end to the vibrating stand, while the opposite end remained free, mimicking the vibrations of a helicopter blade. To the free end of the plate, a load with predefined mass was attached at point A (see Figure 3b). The values of the Bragg wavelength shifts of all FBG detectors were acquired for seven different load cases (L = 0, 50, 60, . . . , 100 g). It must be noted that in this preliminary research, the authors did not pursue the goal to compare the vibrations with the real-life reference, which will be the aim of the subsequent studies. The measurements were taken during a certain time interval (~8 s) so that 800 points of data were registered with the rate of~100 Hz for each load case. The data were collected using an FBG interrogation device developed at KNRTU-KAI based on the Ibsen I-MON 512 interrogation monitor.

Treatment Procedure
Based on the treatment procedure explained in Section 3, we received seven matrices. Each matrix included 6 columns and 800 data points corresponding to different loads (L = 0, 50, 60, …, 100 g), which were subjected to the same vibrations measured in nm as the wavelength shifts of the fiber-optic detectors. Each column inside the matrix 800 × 6 included the recorded vibrations associated with the fixed detector. Each detector shown in Figure 3b was tightly associated with a number of columns (their correspondence to each other shown below), and it was located on the plate imitating the helicopter blade.
In order to select the desired value of the smoothing window w figuring in Expression (2), we chose a simple criterion that the low-frequency trends obtained after POLS should be almost uncorrelated with each other. If we fix the interval of correlations [−0.5 × 10 −3 , 0.5 × 10 −3 ] for the smoothed trends and require that where Ysm(L,w) is determined by Expression (2)

Treatment Procedure
Based on the treatment procedure explained in Section 3, we received seven matrices. Each matrix included 6 columns and 800 data points corresponding to different loads (L = 0, 50, 60, . . . , 100 g), which were subjected to the same vibrations measured in nm as the wavelength shifts of the fiber-optic detectors. Each column inside the matrix 800 × 6 included the recorded vibrations associated with the fixed detector. Each detector shown in Figure 3b was tightly associated with a number of columns (their correspondence to each other shown below), and it was located on the plate imitating the helicopter blade.
In order to select the desired value of the smoothing window w figuring in Expression (2), we chose a simple criterion that the low-frequency trends obtained after POLS should be almost uncorrelated with each other. If we fix the interval of correlations [−0.5 × 10 −3 , 0.5 × 10 −3 ] for the smoothed trends and require that where Ysm(L,w) is determined by Expression (2)   In order to stress their difference, the same picture Figure 4 (b) is placed on the right-hand side for detector D-2.2 (with load 10 g) and without load (red lines), correspondingly. For other detectors, the obtained pictures look the same and therefore they are omitted. Figure 5a demonstrates the results of the application POLS for w = 150 and detector D-1.1. Figure 5b shows the same result for D-2.2 at w = 150, as well. Analysis of these curves shows that for construction of the calibration curves showing the dependence of each detector with respect to the applied load one can select the differences between mean values corresponding to these smoothed curves.   In order to stress their difference, the same picture Figure 4 (b) is placed on the right-hand side for detector D-2.2 (with load 10 g) and without load (red lines), correspondingly. For other detectors, the obtained pictures look the same and therefore they are omitted. Figure 5a demonstrates the results of the application POLS for w = 150 and detector D-1.1. Figure 5b shows the same result for D-2.2 at w = 150, as well. Analysis of these curves shows that for construction of the calibration curves showing the dependence of each detector with respect to the applied load one can select the differences between mean values corresponding to these smoothed curves.  In order to stress their difference, the same picture Figure 4 (b) is placed on the right-hand side for detector D-2.2 (with load 10 g) and without load (red lines), correspondingly. For other detectors, the obtained pictures look the same and therefore they are omitted. Figure 5a demonstrates the results of the application POLS for w = 150 and detector D-1.1. Figure 5b shows the same result for D-2.2 at w = 150, as well. Analysis of these curves shows that for construction of the calibration curves showing the dependence of each detector with respect to the applied load one can select the differences between mean values corresponding to these smoothed curves.    It is interesting to analyze also 13 parameters (given by Expressions (8), (11), and (12)) that form the desired 3D surface. In order to compare these parameters with each other, we form the corresponding SRAs after applying the POLS procedure for different loads. Figure 7a demonstrates 3 curves calculated for the load L = 0 g. In order to notice the possible monotone behavior of these curves, we use the corresponding SRAs for different loads. On the right-hand side (Figure 7b), we place similar triple curves corresponding to L = 50 g. Other curves for the loads L = 60-100 g look similar and therefore are not shown. In order to compare these parameters with each other, we take the mean values for the correlators of the first order (corresponding to the gravity center), then correlators of the second, third, and fourth orders, correspondingly. For us, it is important to detect a monotone behavior among these correlators. Figure 8a demonstrates six calibration It is interesting to analyze also 13 parameters (given by Expressions (8), (11), and (12)) that form the desired 3D surface. In order to compare these parameters with each other, we form the corresponding SRAs after applying the POLS procedure for different loads. Figure 7a demonstrates 3 curves calculated for the load L = 0 g. In order to notice the possible monotone behavior of these curves, we use the corresponding SRAs for different loads. On the right-hand side (Figure 7b), we place similar triple curves corresponding to L = 50 g. Other curves for the loads L = 60-100 g look similar and therefore are not shown. It is interesting to analyze also 13 parameters (given by Expressions (8), (11), and (12)) that form the desired 3D surface. In order to compare these parameters with each other, we form the corresponding SRAs after applying the POLS procedure for different loads. Figure 7a demonstrates 3 curves calculated for the load L = 0 g. In order to notice the possible monotone behavior of these curves, we use the corresponding SRAs for different loads. On the right-hand side (Figure 7b), we place similar triple curves corresponding to L = 50 g. Other curves for the loads L = 60-100 g look similar and therefore are not shown. In order to compare these parameters with each other, we take the mean values for the correlators of the first order (corresponding to the gravity center), then correlators of the second, third, and fourth orders, correspondingly. For us, it is important to detect a monotone behavior among these correlators. Figure 8a demonstrates six calibration In order to compare these parameters with each other, we take the mean values for the correlators of the first order (corresponding to the gravity center), then correlators of the second, third, and fourth orders, correspondingly. For us, it is important to detect a monotone behavior among these correlators. Figure 8a demonstrates six calibration curves calculated for the mean values of the first-order correlators. All curves are almost monotone, and we see an increase in these curves with respect to the applied load.  Analysis of these figures shows that they have different sensitivity to the applied load. Only the mean correlators of the first order having monotonic property can be used for calibration purposes. It presents an interesting possibility to give also the plots of 3D surfaces at different loads. It will reflect a unique combination of the 13 correlation parameters and this surface will serve as a specific fingerprint facilitating an initial analysis of the available data. Figure 10a shows the 3D surface for the L = 0 g, while Figure 10b presents the surface for L = 50 g. The 3D surface for the L = 90 g is demonstrated in Figure 11a. Figure 10b shows the surface for L = 100 g. In all these 3D surfaces, the unit of the OZ axis is conserved, and it is given in the units coinciding with initial data (values of vibrations in our case). As one can notice from the comparison of these figures, they are different in their values and have different forms. Analysis of these figures shows that they have different sensitivity to the applied load. Only the mean correlators of the first order having monotonic property can be used for calibration purposes. It presents an interesting possibility to give also the plots of 3D surfaces at different loads. It will reflect a unique combination of the 13 correlation parameters and this surface will serve as a specific fingerprint facilitating an initial analysis of the available data. Figure 10a shows the 3D surface for the L = 0 g, while Figure 10b presents the surface for L = 50 g. The 3D surface for the L = 90 g is demonstrated in Figure  11a. Figure 10b shows the surface for L = 100 g. In all these 3D surfaces, the unit of the OZ axis is conserved, and it is given in the units coinciding with initial data (values of vibrations in our case). As one can notice from the comparison of these figures, they are different in their values and have different forms.  Analysis of these figures shows that they have different sensitivity to the applied load. Only the mean correlators of the first order having monotonic property can be used for calibration purposes. It presents an interesting possibility to give also the plots of 3D surfaces at different loads. It will reflect a unique combination of the 13 correlation parameters and this surface will serve as a specific fingerprint facilitating an initial analysis of the available data. Figure 10a shows the 3D surface for the L = 0 g, while Figure 10b presents the surface for L = 50 g. The 3D surface for the L = 90 g is demonstrated in Figure  11a. Figure 10b shows the surface for L = 100 g. In all these 3D surfaces, the unit of the OZ axis is conserved, and it is given in the units coinciding with initial data (values of vibrations in our case). As one can notice from the comparison of these figures, they are different in their values and have different forms.
(a) (b) Figure 10. (a) The 3D surface for the L = 0 g; (b) the 3D surface for L = 50 g. Figure 10. (a) The 3D surface for the L = 0 g; (b) the 3D surface for L = 50 g.

Discussion and Conclusions
The proposed data treatment method based on the procedure of the optimal linear smoothing (POLS) allowed us to eliminate the noise induced by vibrations of the experimental plate, imitating the helicopter blade in flight conditions. As the result of the treat- Figure 11. (a) The 3D surface for the L = 90 g; (b) the 3D surface for L = 100 g.

Discussion and Conclusions
The proposed data treatment method based on the procedure of the optimal linear smoothing (POLS) allowed us to eliminate the noise induced by vibrations of the experimental plate, imitating the helicopter blade in flight conditions. As the result of the treatment procedure, monotonic calibration curves were obtained for each detector, making it possible to distinguish the increase of load applied to the plate with an increment of 10 g, which otherwise would be impracticable. It should be noted that the calibration curves comprising mean values of POLS-treated wavelength shift data of the detectors (Figure 6a) demonstrated higher mean deformation of the detector 2.2 (and even higher in case of the detector 2.3) than the one of the detector 2.1, which does not correspond to the a priori assumption that the detectors 2.1 and 1.1 (which are the nearest to the fixed end of the plate) should be subjected to the highest deformations, while 2.3 and 1.3 should experience the lowest strain. More detailed analysis based on a complete set of data allows concluding that detector 2.1 had lower sensitivity, in comparison with others, which can be caused by relatively weak adhesive bonding to the surface of the plate.
Highly likely, this peculiarity is related to the chosen material and formation in it the deformation waves.
However, the calibration curves calculated for mean values of the first-order correlators ( Figure 8a) do agree with the mentioned assumption, showing the higher values at the detectors 1.1 and 2.1, and the lowest ones at the detectors 1.3 and 2.3.
As one can infer from the content given above, the second method associated with 3D invariants takes into account 13 quantitative parameters associated with a combination of different moments and their intercorrelations up to the fourth-order inclusive. They reflect more sophisticated behavior of the moments and their intercorrelations evoked by the applied load. Figures 8b and 9a,b show their peculiarities. These 13 parameters allow the calculation of the desired 3D surface that can serve as a specific fingerprint, differentiating one set of initial data from another one. Figure 10a,b and Figure 11a,b vividly demonstrate these peculiarities.
Therefore, the combination of these two methods that were used to analysis of real load data can be applied successfully to other data, because they are general, do not use a specific model, and are free from treatment errors.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author. The data are not publicly available due to the rules of our contract conditions with our customer.