Understanding Modal Test

This chapter describes the calculation princeples and recorded data channels.

Signal and Data flow

Fig. 24 Signal and Data flow

Calculation modes

In this chapter, the calculation modes for all relevant results will be explained.

FRF - Frequency Response function

The FRF is calculated between each excitation/response pair in mono or triaxial mode, depending on the sensor. These calculations are available as single-event transfer functions and averaged transfer functions for each excitation/response pair of the same direction. The data can be displayed as amplitude spectrum, phase spectrum, real part, and imaginary part in the instrument properties tab of the modal test screen. The information is stored as a complex frequency spectrum and can be exported in various formats; see chapter Recording and exporting data .

FRF Calculation - Transfer function & bode plot

In the ”Trigger & FRF” Fig. 16 tab, one calculation type has to be chosen to define a modal test setup. The variables A and B stand for the complex spectra of the input (excitation) and output (response). Additionally, A’ and B’ present the complex conjugate spectra of the input and output channels.

H(f) = \frac{B(\omega)}{A(\omega)}

A(\omega) Complex spectra of input channel.

B(\omega) Complex spectra of output channel.

In case the output signal is noisy, the H_1 algorithm can be used:

H_1(\omega) =
\frac{B(\omega)}{A(\omega)}
*
\frac{A'(\omega)}{A'(\omega)}
=
\frac{S_{AB}(\omega)}{S_{AA}(\omega)}
=
\frac{G_{AB}(\omega)}{G_{AA}(\omega)}

A'(\omega) Conjugate complex spectra of input channel

B'(\omega) Conjugate complex spectra of output channel

If the input signal is noisy, the H_2 algorithm can be used:

H_2(\omega) =
\frac{B(\omega)}{A(\omega)}
*
\frac{B'(f)}{B'(f)}
=
\frac{S_{BB}(\omega)}{S_{AB}(\omega)}
=
\frac{G_{BB}(\omega)}{G_{AB}(\omega)}

In case both input and output signal are noisy, H_3 or H_v can be used:

H_3(\omega) =
\frac{H_1(\omega)+H_2(\omega)}{2}

H_v(\omega) =
\sqrt{H_1(\omega)}
*
\sqrt{H_2(\omega)}

Coherence

The Coherence function compares the response function of the active group (excitation and response) with previously recorded response functions and returns values between 0 and 1. A value of 1 indicates an identical response function. Typically, measurements with a coherence below 0.8 are rejected and must be repeated.

Coherence calculation for two signals:

\gamma_{XY}^{2}(f) =
\frac{\left|G_{XY}(f)\right|^{2}}
{\left\langle G_{XX}(f)\right\rangle * \left\langle G_{YY}(f)\right\rangle}

G_{XY}(f) Cross-spectral density

G_{XX}(f) and G_{YY}(f) Auto spectral density

Coherence calculation of a signal group and another signal:

\gamma^{2}_{[A]B}(f) =
\frac
{G^{H}_{[A]B}(f) * G^{-1}_{[A][A]}(f) * G_{[A]B}(f)}
{G_{BB}(f)}

The following applies to this formula:

\gamma^{2}_{[A]B}(f) \in [0,1]

A value of 1 indicates a linear relationship between all signals in the selected group [a_{1}(t);\ldots;a_{n}(t)] and the signal b(t).

Whereas a value of 0 indicates no correlation between the signals a_{1}(t);\ldots;a_{n}(t) and b(t).

MIF - Mode Indicator functions

The Mode Indicator Function (MIF) evaluates the Frequency Response Function (FRF) of all excitation pathways. Values approaching 0 indicate the presence of a mode, while values approaching 1 indicate the absence of a mode. Additionally, the MIF for each excitation–FRF pair is calculated in the Excitation tab.

\mathrm{MIF}(f)
=
\frac{\sum_{i=1}^{N}\mathrm{Re}\{H(F)\}^{2}}
     {\sum_{i=1}^{N}|H(F)|^{2}}

In the current implementation of the modal test, no model validation is implemented.

SDOF circle fit

The SDOF (Single Degree of Freedom) circle fit is a numerical method to estimate the natural frequency and loss factor. To calculate the loss factor, a spectral region, see Fig. 25 , is chosen where one natural frequency is assumed and transformed into a Nyquist plot, see Fig. 26 . In the Nyquist diagram, the point at which the circle crosses the Y-axis represents the exact natural frequency.

Frequencies with the amplitude fit method using the amplitude response of an ideal single-mass oscillator

Fig. 25 Frequencies with the amplitude fit method using the amplitude response of an ideal single-mass oscillator

Nyquist plot for SDOF with the locations of the natural frequency :math:`\omega_0` and other frequencies :math:`\omega_a` and :math:`\omega_b`

Fig. 26 Nyquist plot for SDOF with the locations of the natural frequency \omega_0 and other frequencies \omega_a and \omega_b

Subsequently, two known frequencies and their angle can be used to calculate the loss factor. The spectral region of interest must only contain one natural frequency.

\eta
=
\frac{\omega_{0}^{2}-\omega_{b}^{2}}
     {\omega_{0}^{2}\left[\tan\left(\frac{\theta}{2}\right)+\tan\left(\frac{\phi}{2}\right)\right]}
\approx
\frac{2(\omega_{a}-\omega_{b})}
     {\omega_{0}\left[\tan\left(\frac{\theta}{2}\right)+\tan\left(\frac{\phi}{2}\right)\right]}

\eta Loss factor

\omega_0 Natural frequency (Eigenfrequenz)

\omega_a Frequency above and

\omega_b Frequency below the natural frequency \omega_0

H(\omega) Frequency response

An example of the OXYGEN implementation and application is shown in chapter SDOF circle fit .