A versatile algorithm for estimating natural frequencies with high accuracy
GilbertRainer Gillich^{1} , Dorian Nedelcu^{2} , ConstantinIoan Barbinta^{3} , Nicoleta Gillich^{4}
^{1, 2, 3, 4}Universitatea “Eftimie Murgu” din Resita, Resita, Romania
^{1, 2}Corresponding authors
Vibroengineering PROCEDIA, Vol. 27, 2019, p. 3742.
https://doi.org/10.21595/vp.2019.20946
Received 11 August 2019; accepted 28 August 2019; published 30 September 2019
41st International Conference on Vibroengineering in Leipzig, Germany, September 30  October 1, 2019
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 interline 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.
 An algorithm for postprocessing 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 superresolution.
 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, superresolution.
1. Introduction
Features extracted from the vibration response are used since decades to characterize structures [13]. One of the main beneficiaries of structural identification is the field of nondestructive detection techniques. Currently, there are numerous methods for detecting defects based on vibration measurement [48]. 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 [1113] or three spectral lines [1416], 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 [1820]. 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 interline 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 nondestructive testing by standard algorithms based on the Fourier transforms is the short length of the acquired signals, which results in a raw frequency resolution $\mathrm{\Delta}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 interline 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 simpletouse 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
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 ${f}_{Ny}$ found in the DFT is the half of the total number of samples $N$, and usually much bigger than the biggest frequency ${f}_{max}$ which is targeted in the analysis.
For the selected ${f}_{max}$ the DFT is calculated and plotted as shown in Fig. 2. In the plot area boxes indicating the frequencies that are local maxima (${f}_{1}$, ${f}_{2}$ and ${f}_{3}$ 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
By selecting a frequency, let’s say ${f}_{1}$, 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 $2{T}_{1}=2/{f}_{1}$ 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 upperright 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 ${f}_{1}$, 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 ${f}_{1}$ 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
a)
b)
c)
d)
e)
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/s^{2}, 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}_{real}$ [Hz]  ${f}_{DFT}$ [Hz]  ${f}_{DFT1}$ [Hz]  ${f}_{DFT+1}$ [Hz]  $n$  ${f}_{est}$ [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 ${f}_{ireal}$ and ${f}_{iDFT}$  Estimated frequency  Error between ${f}_{ireal}$ and ${f}_{iest}$ 
${f}_{ireal}$ [Hz]  ${f}_{iDFT}$ [Hz]  [%]  ${f}_{iest}$ [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 vibrationbased 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 socalled noise) and may slightly affect the precision of the estimation. By increasing the signal length this effect is diminished.
References
 Zenzen R., Belaidi I., Khatir S., Wahab M. A. A damage identification technique for beamlike and truss structures based on FRF and Bat algorithm. Comptes Rendus Mécanique, Vol. 346, Issue 12, 2018, p. 12531266. [Publisher]
 Gillich G. R., Praisach Z. I., OnchisMoaca 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. 7681. [Search CrossRef]
 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. 386401. [Publisher]
 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. 335352. [Publisher]
 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. 151154. [Search CrossRef]
 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. 89114. [Publisher]
 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. 101108. [Search CrossRef]
 Zhou Y. L., Abdel Wahab M. Cosine based and extended transmissibility damage indicators for structural damage detection. Engineering Structures, Vol. 141, 2017, p. 175183. [Publisher]
 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. 8287. [Search CrossRef]
 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. 23112329. [Publisher]
 Grandke T. Interpolation algorithms for discrete Fourier transforms of weighted signals. IEEE Transactions on Instrumentation and Measurement, Vol. 32, 1983, p. 350355. [Publisher]
 Quinn B.G. Estimating frequency by interpolation using Fourier coefficients. IEEE Transactions on Signal Processing, Vol. 42, 1994, p. 12641268. [Publisher]
 Jain V. K., Collins W. L., Davis D. C. Highaccuracy analog measurements via interpolated FFT. IEEE Transactions on Instrumentation and Measurement, Vol. 28, 1979, p. 113122. [Publisher]
 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. 4348. [Publisher]
 Voglewede P. Parabola approximation for peak determination. Global DSP Magazine, Vol. 3, Issue 5, 2004, p. 1317. [Search CrossRef]
 Jacobsen E., Kootsookos P. Fast, accurate frequency estimators. IEEE Signal Proccessing Magazine, Vol. 24, Issue 3, 2007, p. 123125. [Publisher]
 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. 1518. [Search CrossRef]
 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. 233242. [Publisher]
 Onchis Moaca D., Gillich G. R., Frunza R. Gradually improving the readability of the timefrequency spectra for natural frequency identification in cantilever beams. European Signal Processing Conference, 2012, p. 809813. [Search CrossRef]
 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. 637652. [Search CrossRef]
 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. 123128. [Publisher]
Cited By
Lecture Notes in Civil Engineering
Nicoleta Gillich, David Lupu, Codruta Hamat, GilbertRainer Gillich, Dorian Nedelcu

2021
