The following article will attempt to explain the basic theory of the frequency response function. This basic theory will then be used to calculate the frequency response function between two points on a structure using an accelerometer to measure the response and a force gauge hammer to measure the excitation.
Fundamentally a frequency response function is a mathematical representation of the relationship between the input and the output of a system.
So for example the frequency response function between two points on a structure. It would be possible to attach an accelerometer at a particular point and excite the structure at another point with a force gauge instrumented hammer. Then by measuring the excitation force and the response acceleration the resulting frequency response function would describe as a function of frequency the relationship between those two points on the structure.
The basic formula for a frequency response function is
Where H(f) is the frequency response function.
And Y(f) is the output of the system in the frequency domain.
And where X(f) is the input to the system in the frequency domain.
Frequency response functions are most commonly used for single input and single output analysis, normally for the calculation of the H1(f) or H2(f) frequency response functions. These are used extensively for hammer impact analysis or resonance analysis.
The H1(f) frequency response function is used in situations where the output to the system is expected to be noisy when compared to the input.
The H2(f) frequency response function is used in situations where the input to the system is expected to be noisy when compared to the output.
Additionally there are other possibilities, but they are outside of the scope of this article.
H1(f) or H2(f) can be used for resonance analysis or hammer impact analysis. H2(f) is most commonly used with random excitation.
The breakdown of H1(f) is as follows,
Where H1(f) is the frequency response function.
And Sxy(f) is the Cross Spectral Density in the frequency domain of X(t) and Y(t)
And where Sxx(f) is the Auto Spectral Density in the frequency domain of X(t).
In very basic terms the frequency response function can be described as
The breakdown of H2(f) therefore is as follows,
Where H2(f) is the frequency response function.
And Syx(f) is the Cross Spectral Density in the frequency domain of Y(t) and X(t)
And where Syy(f) is the Auto Spectral Density in the frequency domain of Y(t)
In very basic terms the frequency response function can be described as
In the following example we will discuss and show the calculation of the H1(f) frequency response function.
The excitation or input would be the force gauge instrumented hammer, as shown in Figure 1 as a time history.
Figure 1: X(t)
In this case the response or output would be the accelerometer, as shown in Figure 2.
Figure 2: Y(t)
However as discussed earlier the frequency response function is a frequency domain analysis, therefore the input and the output to the system must also be frequency spectra. So the force and acceleration must be first converted into spectra.
The first part of the analysis requires the Cross Spectral Density of the input and output, this is Sxy(f). This is calculated using the response as the first input and the excitation as the second input to the Cross Spectral Density Analysis in DATS the result is shown in Figure 3. Were Sxy(f) being calculated for use with H2(f) for example, then the excitation would be the first input and the response the second input to the Cross Spectral Density Analysis in DATS.
Figure 3: Sxy(f)
Next the Auto Spectral Density of the input, or excitation signal is required. This is calculated using the Auto Spectral Density Analysis in DATS, this analysis is sometimes known as Auto Power, the result of which is shown in Figure 4, this is Sxx(f).
Figure 4: Sxx(f)
The Cross Spectrum is then divided by the Auto Spectrum and the resulting frequency response function is shown in Figure 5.
Figure 5: H1(f)
The entire analysis as used in DATS.toolbox is shown in Figure 6, the data flow from the original input and output, force and response, can be seen through to the frequency response function. The DATS software does, of course, provide a single step transfer function analysis. We have deliberately used the long-hand form below to illustrate the steps in this article.
Figure 6: Complete DATS worksheet (Click to expand)
It is necessary to understand that for the purposes of understanding and clarity in this article some important steps have been glossed over, windowing of the input for example, to allow the basic understanding of what makes up the frequency response function.
Hello,
Normally people don’t expect to see an Autospectrum, or Power Spectral Density on a linear scale, and it looks so noisy, that I suspect the impact was recorded using a Hanning weighting – a common mistake.
The end result, – the transfer function is also odd, to my eyes.
I would have expected a bode plot.
I am aware that linear display is useful for fatigue and stress purposes, but the frequency range is too high for this to be the application.
Hi Ken,
Thanks for posting on our blog.
I’d like to respond to your points if I may.
I understand your point of view, the Power Spectral Density is normally shown on a logarithmic scale (Y axis) as you suggest.
However for the purposes of this article it is not, we have shown it as linear. Here in this article we are trying to explain and discuss the basics in a straight forward simple fashion.
With regards to your points about noise and windowing, we apply a window to the force signal. We also apply a window to the response if it is very noisy, but we prefer not to do this.
The window is an exponential window, like a transient window for example. It is not a hanning window or similar.
The structure in question is a very simple structure, if a little noisy, and you’re correct it’s not being used for fatigue analysis.
With regards to the final output, our DATS software can quickly switch between Modulus & Phase and Real & Imaginary. So it is possible show in the Bode format you mention.
what can be infered from figure 5,please explain..
Hello Eapen,
I’m not sure anything can be inferred from Figure 5, it simply shows the Frequency Response Function for this particular test.
how can v obtain auto spectral density of velocity from auto spectral density of the displacement
Hi Kiran,
Thanks for asking a question on our blog.
This conversion is one of the fundamental laws of Newton’s physics.
It is possible to convert from Acceleration to Velocity to Displacement using calculus, specifically integration.
Our Prosig software performs this conversion in an advanced fashion to take account of the constant and remove this error from the results.
You should keep in mind that the original time series is required for this conversion.
If you would like to discuss this feature further please feel free to contact us directly.
Hi
Just say you have several accelerometers on a complex vibrating structure. Each accelerometer has a slightly different frequency spectrum. Lets also pretend that you have a mic at some distance from this vibrating structure, and that you are trying to locate the particular component or part of the structure that is responsible for a radiating a particular tonal frequency. Would the cross spectrum be valuable in identifying which accelerometer is the culprit of this offending sound?
Many thanks
Hello Paul,
Thanks for asking a question on our blog.
We can see that you have a good understanding of what you are trying to do.
In our opinion the Cross Spectrum is a step in the correct direction, but you should go one step further and calculate the Coherence between each vibration source and the microphone response, you can do this very easily with a software package like DATS.
To actually rank each vibration with respect to the microphone you should consider something like the Source Contribution Package, again as part of DATS.
The Source Contribution Analysis uses a method called Singular Value Decomposition. The Singular Value Decomposition computation produces an eigenvector matrix, this matrix is used to derive the cross spectra between the vibration references and the measured sound response. These cross spectra are then used to calculate the Reference Related Auto spectra at the response position. Each Reference Related Auto spectrum is related to the coherent contributions from the particular references and source input.
If you would like to discuss this further please feel free to contact us directly.
Going back to Paul’s hypothetical situation. Suppose that the accellerometers are measuring the start of an event that will be producing the noise that is being measured, but that there is no real reason to expect there to be any frequency correlation between them, would any conventional analysis work then?
For example, if you have 5 drums and a lightbeam sensor for each drum that triggered just before the drumstick hit each one, and that the drums are being hit in a regular sequence. How might one analyse the data (one mic channel and 5 drumstick sensors) to determine which drum was loudest/highest-pitch etc?
Hi Andy,
Thanks for asking a question on our blog.
I think there is some difference between your question and Pauls.
In Paul’s case the vibration responses are all related.
In your case they are not.
Therefore the same analysis would not apply.
I think you are trying to find which drum gives a certain frequency response when hit. The easiest thing to do would be to use a microphone, mounted about the drums and hit each drum in sequence. Starting and stopping a data capture for each drum. Then simply frequency analyse each of these data captures. You might need to do it several times to build up an average for each drum as you may hit it differently each time.
This will give you the full frequency spectrum for each drum.
You wouldn’t need any accelerometers at all.
Sorry, perhaps I took my analagous situation too far. This is actually exactly the same question as Paul’s, but we have a slightly different conception of the nature of the data. I should perhaps have chosen a better analogy, as the trigger pulses occur at 20-150Hz which would be a pretty inhuman rate of drumstick operation.
Hi Andy,
Thanks for the additional information, but I still don’t understand your application and can’t comment on it unless you can explain what your trying to do. Perhaps you could contact us directly to discuss?