A versatile algorithm for estimating natural frequencies with high accuracy

Gilbert-Rainer Gillich1 , Dorian Nedelcu2 , Constantin-Ioan Barbinta3 , Nicoleta Gillich4

1, 2, 3, 4Universitatea “Eftimie Murgu” din Resita, Resita, Romania

1, 2Corresponding authors

Vibroengineering PROCEDIA, Vol. 27, 2019, p. 37-42. https://doi.org/10.21595/vp.2019.20946
Received 11 August 2019; accepted 28 August 2019; published 30 September 2019

Copyright © 2019 Gilbert-Rainer Gillich, et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Creative Commons License
Abstract.

Detection of damage requires an accurate estimation of the natural frequencies of the monitored structure. This paper introduces an algorithm implemented in Python which improves the frequency readability by increasing the number of spectral lines without requiring a signal extension in the time domain. We achieve this by overlapping several spectra calculated from the acquired signal repeatedly shortened. In this way, the overlapped spectrum gets an increased number of spectral lines. The dense mesh of spectral lines permits us to obtain a fine frequency resolution without being necessary an extension of the signal in the time domain. The high density of the spectral lines ensures a sufficient number of points on the main lobes that permits performing an efficient quadratic polynomial interpolation to find the maximizer. It represents the amplitude of the real frequency and is typically located on an inter-line position, thus cannot be found by standard frequency estimation. We implemented the algorithm in Python and tested it successfully for generated signals, containing one or more harmonics, with known frequencies.

A versatile algorithm for estimating natural frequencies with high accuracy

Highlights
  • An algorithm for post-processing signals in order to increase the spectral line density is proposed.
  • It is based on the repetitive truncation of the signal and performance of a DFT.
  • All obtained spectra are superposed, finally being obtained an overlapped spectrum with super-resolution.
  • The envelope around the targeted frequency's amplitudes is found and its maximum is calculated.
  • The frequency that corresponds to the maximum amplitude is the true frequency.

Keywords: vibration, natural frequency, structure, frequency estimation, Python, algorithm, quadratic interpolation, Discrete Fourier Transform, overlapped spectrum, super-resolution.

1. Introduction

Features extracted from the vibration response are used since decades to characterize structures [1-3]. One of the main beneficiaries of structural identification is the field of non-destructive detection techniques. Currently, there are numerous methods for detecting defects based on vibration measurement [4-8]. The most common modal parameter used to detect defects is the natural frequency, because it can be easily measured and involves the use of simple and robust equipment. The mathematical relationship between the damage dimensions and the frequency drop is proposed in [9]. Because the damage has a low influence on the frequency changes, in order to observe the occurrence of damage as soon as possible [10], accurate frequency estimation is requested. The simplest way to achieve this target is using an interpolation method based on two spectral lines [11-13] or three spectral lines [14-16], but the results have been proved to be unsatisfactory for detecting damage in the early stage [17].

In prior research, we have developed several techniques to improve the frequency readability [18-20]. In our latest research, we have developed a structure excitation procedure that provides the highest amplitude for the targeted vibration mode [21]. In this paper, we complete the procedure with an original algorithm based on quadratic polynomial interpolation that helps finding the real frequency at an inter-line position. The novelty of the algorithm consists in applying the interpolation to spectral lines belonging to three different spectra obtained by stepwise shortening the acquired signal and finding the abscissa of the maximum obtained from the contrived polynomial function. The algorithm is implemented in Python as an application (PyFEST) and we tested its functionality for a signal containing three harmonics with known frequencies and amplitudes. We managed to estimate both frequencies and amplitudes with high accuracy and found that the application is easy to use and versatile, allowing the user to choose the optimal input parameters during the frequency evaluation process.

2. The proposed algorithm and its implementation in Python

The main disadvantage of estimating the natural frequencies for non-destructive testing by standard algorithms based on the Fourier transforms is the short length of the acquired signals, which results in a raw frequency resolution Δf. While this resolution is somehow controllable for the Discrete Fourier Transform (DFT) but for the Fast Fourier Transform (FFT) it isn’t because the length of the analyzed signal depends on the specific number of samples that must be a power of 2. In consequence, the real frequency is typically located at a spectral inter-line position, and therefore an accurate estimation is not possible.

To overcome this inconvenient, we developed an algorithm which performs a quadratic interpolation based on three points that belong to three spectra, each spectrum being obtained for the original signal with a different time length. To be convinced that we have chosen the right points for interpolation, we execute a signal truncation until it is shortened to the length corresponding to two periods T of the target frequency f. At each step, two samples are extracted from the signal, resulting in a large number of truncations for each T. In this way we obtain numerous amplitudes associated with the corresponding frequencies, which usually form three distinct curves: the first is a portion of a cycle, the second represents a whole cycle and the last again a portion of a cycle. Finally, we extract three points from the curve with the biggest amplitude which has a maximum and two neighbors and perform the interpolation. The achieved maximum of the polynomial represent the amplitude and its associated abscissa the real frequency. We implemented the algorithm in Python and a simple-to-use application resulted. The actions performed to run the program are presented in Table 1, along with the displayed results.

In the first step, the user select a CSV (preferably due to the increased data import speed) or Excel file containing the analyzed signal. The signal is displayed as shown in in Fig. 1.

Table 1. The main steps of the algorithm

Step
Action
Result
1
Select input data file
The chart representing the vibration signal is displayed together with a box to specify the upper frequency limit
2
Select the upper frequency limit
The selected portion of the DFT is displayed and frequencies of the maximums and their neighboring frequencies are highlighted
3
Select the targeted frequency
The overlapped spectrum, consisting of the maxima resulted from the DFT calculations for the repeatedly truncated signal grouped by vibration cycles, is displayed
4
Select proper number of cycles/periods
The frequency corresponding to the maximum calculated by interpolation and the amplitudes resulting from DFT and recalculated for the real signal are displayed

Fig. 1. The signal to be analyzed and the box for choosing the upper frequency limit to be analyzed

 The signal to be analyzed and the box for choosing the upper frequency limit to be analyzed

Along with the original signal representing the acquired or generated data, a box that permits the user to select the upper frequency limit is displayed in the center of the diagram. The frequency selected here will define the number of spectral lines for which the DFT is calculated. This option is available to save computing time, because the highest frequency fNy found in the DFT is the half of the total number of samples N, and usually much bigger than the biggest frequency fmax which is targeted in the analysis.

For the selected fmax the DFT is calculated and plotted as shown in Fig. 2. In the plot area boxes indicating the frequencies that are local maxima (f1, f2 and f3 in this case) and the frequencies for the neighbor spectral lines are displayed. Here, the user can select a frequency to be analyzed.

Fig. 2. DFT obtained for the imported signal and boxes indicating the frequencies that are local maxima and their neighbors

 DFT obtained for the imported signal and boxes indicating the frequencies  that are local maxima and their neighbors

By selecting a frequency, let’s say f1, the application perform DFT for the signal containing the original number of samples N. Afterward, two samples are cut consecutively from the signal until the original signal is shortened with 2T1=2/f1 seconds. For all achieved signals the DFT is performed again and the biggest amplitude is plotted in a new graph, see Fig. 3(a), 3(c) and 3(e). The points are grouped in function of the cycle/period to which they belong. This information is displayed on the upper-right corner of the plot area. Note that, the DFT values are calculated just for the spectral lines that belong to the frequency domain framed by the neighbors of the targeted frequency f1, again to save computing effort.

The user can select the proper number of cycles in the box located in the middle of the plot area. The implicit value is the median value found for the number of cycles that is 7 for the presented case in which f1 is targeted. By selecting one of the suggested values, the application sorts all points belonging to the given cycle (Fig. 3(b)) and selects the maximizer and its two neighbors, which are used for interpolation. The polynomial function is found and its maximum is calculated. This maximum represents the DFT amplitude and its abscissa the real frequency. In order to find the real amplitude, the application divides the DFT amplitude by N/2. All three features, the real frequency, DFT amplitude and real amplitude are indicated at the top of the plot area in Fig. 3(b), 3(d) and 3(f).

For the frequency estimation process is crucial to select the right number of cycles. In order to obtain reliable points for interpolation, it is extremely important to achieve at least eight truncated signals, and consequently eight maximizer, that belong to one cycle/period. Also, the number of cycles contained in the original signal should be five, but preferably bigger than ten.

The estimation process must be repeated for all targeted frequencies; to this aim, steps 1 to 4 should be resumed.

Fig. 3. The overlapped DFT in the selected frequency rage and the sorted values for the chosen cycle

 The overlapped DFT in the selected frequency rage and the sorted values for the chosen cycle

a)

 The overlapped DFT in the selected frequency rage and the sorted values for the chosen cycle

b)

 The overlapped DFT in the selected frequency rage and the sorted values for the chosen cycle

c)

 The overlapped DFT in the selected frequency rage and the sorted values for the chosen cycle

d)

 The overlapped DFT in the selected frequency rage and the sorted values for the chosen cycle

e)

 The overlapped DFT in the selected frequency rage and the sorted values for the chosen cycle

f)

3. Results and discussion

The first test concerns harmonic signals with the frequencies indicated in Table 2. All signals have the amplitude A= 1 mm/s2, and are generated with N= 501 samples and have a time length t= 1.25 seconds. The selected parameters, the estimated frequencies and the estimation errors are also presented in Table 2. It can be observed that the errors obtained by involving the Python application based on the original frequency estimation algorithm lead to excellent results, the errors being less than ±0,2 %, and in absolute values the frequencies estimated with this application are much closer than those obtained from a standard DFT.

A second approach considers a test signal composed of three harmonics with the frequencies and amplitudes indicated in Table 3. This signal simulates the first three modes of vibration of a beam. It is generated with N= 501 samples and again has a length of time t= 1.25 seconds.

Table 2. The test results if one harmonic is considered in the experiment

Test number
Real frequency
DFT estimation
Lower neighbor
Upper
neighbor
Number of cycles
Estimated frequency
Error
f r e a l [Hz]
f D F T [Hz]
f D F T - 1 [Hz]
f D F T + 1 [Hz]
n
f e s t [Hz]
[%]
1
6.12
6.40
5.6
7.2
7
6.1142
0.094
2
15.78
16.00
15.2
16.8
19
15.8060
–0.164
3
25.34
25.60
24.8
26.4
31
25.3874
–0.187

Table 3. The test results for the signal composed by three harmonics with different amplitudes

Mode number
Real
frequency
Frequency from DFT estimation
Error between fi-real and fi-DFT
Estimated frequency
Error between fi-real and fi-est
f i - r e a l [Hz]
f i - D F T [Hz]
[%]
f i - e s t [Hz]
[%]
1
6.12
6.40
–4.575
6.1088
0.183
2
17.93
17.60
1.840
17.9596
–0.165
3
29.00
28.80
0.689
29.0613
–0.211

The first harmonic component intentionally coincides with the first sinusoidal signal tested because we intended to see if the estimated frequencies in the composed signals are affected by the other harmonics. Analyzing the results in Tables 2 and 3 we observe that the frequency of 6.12 Hz is precisely identified, even if there is a very small difference between the results obtained for the simple signal and the one composed of three harmonics.

The frequencies of the other two components are also accurately estimated, which proves that the algorithm can be used with confidence in identifying the natural frequencies of beams and plates in mechanical or civil structures.

4. Conclusions

The proposed algorithm for estimating the natural frequencies of structures, for which we have designed the PyFEST application in Python, allows an accurate estimation of the natural frequencies of beams or plates. This application makes possible early assessment of cracks when using the vibration-based detection methods. We tested the application for signals composed by single or multiple harmonics with known frequencies and found that in all cases the frequencies and amplitudes are estimated with high accuracy, the errors being less than ±0.22 %. In fact, the frequencies can be read with an error occurred by the second digit after the dot even for high frequencies, which makes PyFEST extremely reliable.

We also found that, if the signal contains more harmonic components the energy of these components is easily spread along the spectrum (the so-called noise) and may slightly affect the precision of the estimation. By increasing the signal length this effect is diminished.

References

  1. Zenzen R., Belaidi I., Khatir S., Wahab M. A. A damage identification technique for beam-like and truss structures based on FRF and Bat algorithm. Comptes Rendus Mécanique, Vol. 346, Issue 12, 2018, p. 1253-1266. [Publisher]
  2. Gillich G. R., Praisach Z. I., Onchis-Moaca D., Gillich N. How to correlate vibration measurements with FEM results to locate damages in beams. Proceedings of the 4th WSEAS International Conference on Finite Differences – Finite Elements – Finite Volumes – Boundary Elements, 2011, p. 76-81. [CrossRef]
  3. Janeliukstis R., Ručevskis S., Kaewunruen S. Mode shape curvature squares method for crack detection in railway prestressed concrete sleepers. Engineering Failure Analysis, Vol. 105, 2019, p. 386-401. [Publisher]
  4. Sha G., Radzieński M., Cao M., Ostachowicz W. A novel method for single and multiple damage detection in beams using relative natural frequency changes. Mechanical Systems and Signal Processing, Vol. 132, 2019, p. 335-352. [Publisher]
  5. Chioncel C. P., Gillich G. R., Mituletu I. C., Gillich N., Pelea I. Visual method to recognize breathing cracks from frequency change. Romanian Journal of Acoustics and Vibration, Vol. 12, Issue 2, 2015, p. 151-154. [CrossRef]
  6. Dahak M., Toua N., Kharoubi M. Damage detection in beam through change in measured frequency and undamaged curvature mode shape. Inverse Problems in Science and Engineering, Vol. 27, Issue 1, 2019, p. 89-114. [Publisher]
  7. Gillich G. R., Minda P. F., Praisach Z. I., Minda A. A. Natural frequencies of damaged beams – a new approach. Romanian Journal of Acoustics and Vibration, Vol. 9, Issue 2, 2012, p. 101-108. [CrossRef]
  8. Zhou Y. L., Abdel Wahab M. Cosine based and extended transmissibility damage indicators for structural damage detection. Engineering Structures, Vol. 141, 2017, p. 175-183. [Publisher]
  9. Praisach Z. I., Minda P. F., Gillich G. R., Minda A. A. Relative frequency shift curves fitting using FEM modal analyses. Proceedings of the 4th WSEAS International Conference on Finite Differences – Finite Elements – Finite Volumes – Boundary Elements, 2011, p. 82-87. [CrossRef]
  10. Gillich G. R., Maia N., Mituletu I. C., Praisach Z. I., Tufoi M., Negru I. Early structural damage assessment by using an improved frequency evaluation algorithm. Latin American Journal of Solids and Structures, Vol. 12, Issue 12, 2015, p. 2311-2329. [Publisher]
  11. Grandke T. Interpolation algorithms for discrete Fourier transforms of weighted signals. IEEE Transactions on Instrumentation and Measurement, Vol. 32, 1983, p. 350-355. [Publisher]
  12. Quinn B.G. Estimating frequency by interpolation using Fourier coefficients. IEEE Transactions on Signal Processing, Vol. 42, 1994, p. 1264-1268. [Publisher]
  13. Jain V. K., Collins W. L., Davis D. C. High-accuracy analog measurements via interpolated FFT. IEEE Transactions on Instrumentation and Measurement, Vol. 28, 1979, p. 113-122. [Publisher]
  14. Ding K., Zheng C., Yang Z. Frequency estimation accuracy analysis and improvement of energy barycenter correction method for discrete spectrum. Journal of Mechanical Engineering, Vol. 46, Issue 5, 2010, p. 43-48. [Publisher]
  15. Voglewede P. Parabola approximation for peak determination. Global DSP Magazine, Vol. 3, Issue 5, 2004, p. 13-17. [CrossRef]
  16. Jacobsen E., Kootsookos P. Fast, accurate frequency estimators. IEEE Signal Proccessing Magazine, Vol. 24, Issue 3, 2007, p. 123-125. [Publisher]
  17. Minda P. F., Praisach Z. I., Gillich N., Minda A. A., Gillich G. R. On the efficiency of different dissimilarity estimators used in damage detection. Romanian Journal of Acoustics and Vibration, Vol. 10, Issue 1, 2013, p. 15-18. [CrossRef]
  18. Gillich G. R., Mituletu I. C., Praisach Z. I., Negru, Tufoi I. M. Method to enhance the frequency readability for detecting incipient structural damage. Iranian Journal of Science and Technology. Transactions of Mechanical Engineering, Vol. 41, Issue 3, 2017, p. 233-242. [Publisher]
  19. Onchis Moaca D., Gillich G. R., Frunza R. Gradually improving the readability of the time-frequency spectra for natural frequency identification in cantilever beams. European Signal Processing Conference, 2012, p. 809-813. [CrossRef]
  20. Gillich G. R., Mituletu I. C., Negru I., Tufoi M., Iancu V., Muntean F. A method to enhance frequency readability for early damage detection. Journal of Vibration Engineering and Technologies, Vol. 3, Issue 5, 2015, p. 637-652. [CrossRef]
  21. Gillich G. R., Mituletu I. C., Nedelcu D., Hamat C. O. A procedure for an accurate estimation of the natural frequencies of structures. Vibroengineering Procedia, Vol. 19, 2018, p. 123-128. [Publisher]