Experimental Modal Analysis
Experimental Modal Analysis
(OP)
Hi everybody,
Can anyone comment my work, since I'm new in this are? I'm trying to perform SIMO analysis in MAtlab.
I have done an analysis of Euler-Bernouli beam subjected to square force in the middle of the beam in Matlab. Using Wilson's Theta numerical integration method I have calculated the beam response, e.g. accelerations, velocities and displacements for the applied excitation. Time step for calculation is 0.002sec. Load magnitude is 25kN and durations is 0.5 sec. The time sampling interval is 5 sec. The force is applied in the middle of the beam, what coincides with the 6th node.
Now I want to obtain point FRF. I know that FRF is the ratio of the beam response to applied loading, i.e. displacements to square force in frequency domain, but here are some unclarities how to obtain this in Matlab.
I'm not sure what should be the sampling frequency. Is it 1 over time step 0.002 for which I've calculated the beam response, or is it 1 over the time sampling interval of 5sec? I've used the time sampling interval, but I don't know is it correct.
Then, I have performed the fft of the displacements for the node where the load is acting and only for translational DOF. Fft of the force I've also done. The absolute values of both are calculated.
Here is the code which I have used for obtain the FRF's:
%% Calculation of beam response in the middle and applied force in frequency domain
Fs = 1/T; %Sampling frequency in Hertz
timeinterval = time; %time = (0:dt:T)
% FFT of displacement at 6th node for translational DOF and force vector
DISP = fft(displace(11,:));
P = fft(Pn);
% Calculating absolute values of displacements and force in frequency
% domain
DISPabs = abs(DISP);
Pabs = abs(P);
displength = length(DISPabs)/2;
Displacement = DISPabs(1:displength);
Force = Pabs(1:displength);
% Calculating FRFs and plotting
FRF66 = Displacement./Force;
Nyquist = Fs/2;
frequency = linspace(0,Nyquist,displength);
figure
semilogy(frequency,FRF66);
grid on
title ('FRF 66')
xlabel('Frequency (Hz)')
ylabel('Amplitude (dB)')
I have also attached the picture with the results of the beam response,
and picture of the FRF at node 6
Would it be reasonable to calculate the power spectral density with periodogram or welch method instead of doing just fft?
Can anybody comment my work, because I don't know is this exact, am I going in the right way?
Please help!
Can anyone comment my work, since I'm new in this are? I'm trying to perform SIMO analysis in MAtlab.
I have done an analysis of Euler-Bernouli beam subjected to square force in the middle of the beam in Matlab. Using Wilson's Theta numerical integration method I have calculated the beam response, e.g. accelerations, velocities and displacements for the applied excitation. Time step for calculation is 0.002sec. Load magnitude is 25kN and durations is 0.5 sec. The time sampling interval is 5 sec. The force is applied in the middle of the beam, what coincides with the 6th node.
Now I want to obtain point FRF. I know that FRF is the ratio of the beam response to applied loading, i.e. displacements to square force in frequency domain, but here are some unclarities how to obtain this in Matlab.
I'm not sure what should be the sampling frequency. Is it 1 over time step 0.002 for which I've calculated the beam response, or is it 1 over the time sampling interval of 5sec? I've used the time sampling interval, but I don't know is it correct.
Then, I have performed the fft of the displacements for the node where the load is acting and only for translational DOF. Fft of the force I've also done. The absolute values of both are calculated.
Here is the code which I have used for obtain the FRF's:
%% Calculation of beam response in the middle and applied force in frequency domain
Fs = 1/T; %Sampling frequency in Hertz
timeinterval = time; %time = (0:dt:T)
% FFT of displacement at 6th node for translational DOF and force vector
DISP = fft(displace(11,:));
P = fft(Pn);
% Calculating absolute values of displacements and force in frequency
% domain
DISPabs = abs(DISP);
Pabs = abs(P);
displength = length(DISPabs)/2;
Displacement = DISPabs(1:displength);
Force = Pabs(1:displength);
% Calculating FRFs and plotting
FRF66 = Displacement./Force;
Nyquist = Fs/2;
frequency = linspace(0,Nyquist,displength);
figure
semilogy(frequency,FRF66);
grid on
title ('FRF 66')
xlabel('Frequency (Hz)')
ylabel('Amplitude (dB)')
I have also attached the picture with the results of the beam response,

and picture of the FRF at node 6

Would it be reasonable to calculate the power spectral density with periodogram or welch method instead of doing just fft?
Can anybody comment my work, because I don't know is this exact, am I going in the right way?
Please help!





RE: Experimental Modal Analysis
"square force " - what is that?
"The time sampling interval is 5 sec" Bet it isn't, sample interval =1 /sampling frequeny, i think yoiu are talking about your frame length in seconds, hence you should have a resuolution of 0.2 seconds and a nyquist frequency of 1/.002/2=250 Hz
Do you understand Welch? If you do, great, use it.
If not then you are trusting somebody else's black box. Why not take the FT of the input signal and the output signal, and look at the complex ratio of them? That surely is the essence of a transfer function (even if computationally it is inefficient compared with welch), it can be checked for reasonableness, done by hand, and so on.
Cheers
Greg Locock
New here? Try reading these, they might help FAQ731-376: Eng-Tips.com Forum Policies http://eng-tips.com/market.cfm?
RE: Experimental Modal Analysis
Thank you for your answer.
I'm at a learning stage of EMA, that's why I'm struggling that much.
"square force" in my example is force of 25N magnitude and total duration of 0,5 sec. Since I don't have a real recordings from an impact, I have used force of a square shape just for learning purposes.
I have misunderstood the terms sampling interval and my frame length, but now everything is clear.
I'll try with Welch and see what I get. Would it be reasonable to use the exponential window here? I have seen from the literature that this type of window is used for data from transient signals. What about other types for PSD estimation in Matlab, as pburg, pyulear or peig? Which of them are most used when dealing with data of transient signals of short duration?
I have now the complex ratio of the input and output
Maybe I ask stupid questions, sorry, I just want to understand the point.
Cheers!
RE: Experimental Modal Analysis
You desperately need to read a coherent introduction to modal, Ewins is the obvious candidate.
Cheers
Greg Locock
New here? Try reading these, they might help FAQ731-376: Eng-Tips.com Forum Policies http://eng-tips.com/market.cfm?
RE: Experimental Modal Analysis
The FRF is very simple FRF = 1/(Mass(omega0²-omega²)+j*Damper)
Use a Bode diagram for vizualization.
Then, you can add more oscillators.
Beam is a continuous system, i.e. there are an infinite number of oscillators!
From the plot, you can try to retrieve the damper value using the -3db bandpass method, for example...
RE: Experimental Modal Analysis
Cheers
Greg Locock
New here? Try reading these, they might help FAQ731-376: Eng-Tips.com Forum Policies http://eng-tips.com/market.cfm?
RE: Experimental Modal Analysis
The -3dB method is very efficient and helful when the system can be analysed using SIMO technique.
SIMO implies that modes are well separated.
If you can isolate one mode from all the others, then -3db method is a very good tool for estmating damping.
The problem with E2015 is that I don't understand the relationship between his goal (to do SIMO analysis) and the choice of his model.
I think that you should ask to yourself :
What is the objective ? To do SIMO analysis for ...
What is the best structure for this ? Discrete systems or continuous systems ?
If continuous system is choosen, what is the best model ?
What should be the parameters for this model ? Structure (Beam, plate, others), Boundary conditions, Excitation...
What is the best method for solving? Analytical ou numerical method?
RE: Experimental Modal Analysis
Absolute nonesense.
Cheers
Greg Locock
New here? Try reading these, they might help FAQ731-376: Eng-Tips.com Forum Policies http://eng-tips.com/market.cfm?
RE: Experimental Modal Analysis
1. One single excitation can't inject enough energy in the whole structure.
2. Modes are coupled (they are not well separated, several modes exist at the same frequency)
In the second case, you can interpret this from the point of view of the information theory : You don't have enough information about the dynamical behavior of the structure. But by exciting the structure in differents points, you get more information in view to determine the modal model of the structure.
Please don't try to be smarter than me ...
RE: Experimental Modal Analysis
You might be wondering how we validated the extracted mode shapes, frequencies and damping, quite simply we resynthesised the original FRFs, a process which is perhaps the forerunner of some of the more modern fitting techniques. Other validation techniques included performing a modal from a second driving point and checking reciprocity, and typically we'd start a project by doing several mini surveys from different driving points and deciding which driving point gave the best quality data and modal excitation.
Cheers
Greg Locock
New here? Try reading these, they might help FAQ731-376: Eng-Tips.com Forum Policies http://eng-tips.com/market.cfm?
RE: Experimental Modal Analysis
I have done the roving hammer method with eight impacts. Now my results show some sense what means that peaks of the FRF are visible and they correspond to values of natural frequencies. The values of first two peaks are almost identical to the first natural frequencies and they are tapered the most, the third peak value is also almost the same as the third natural frequency but the peak isn't that much tapered.
How can I make other peaks visible on my FRF curve? Should I apply greater force? Anny suggestions on this topic.
I have attached the final FRF curve which I obtained with calculating the ratio of auto power spectrum of output and cross power spectrum of the input and output using cpsd function in MATLAB with hanning window length equal to number of samples.
Also on my plot there are some anomalies. I have zoomed this on the following picture. Is this happening because of the type of window I used? Any suggestions how could I smooth the curve.
Comments and suggestions are welcome.
Many thanks!
Emina B.
RE: Experimental Modal Analysis
Cheers
Greg Locock
New here? Try reading these, they might help FAQ731-376: Eng-Tips.com Forum Policies http://eng-tips.com/market.cfm?
RE: Experimental Modal Analysis
As you can see, tapering decreases with peak number.
It means that the damping is certainly viscous or hysteretic.
You can try to verify for each peak if (Delta f)/f = Constant ( where f is the frequency of the peak and Delta f is the passband located at Ymax/sqrt(2) ).