# Calculate amplitude(displacement) from accelerometer

1. May 6, 2013

### Pivskid

Hi there
Given a result from an accelerometer, mounted on a vibrating machinery.
I would like to be able to calculate the physical amplitude. That is the actual movement of the vibration.
I'm sure it should be possible, but my attempt to manipulate the formula of acceleration has not resulted in valid results.

Given is a sinusodial measurement of 20hz, 4g "peak to peak".

The acceleration fomula is fairly simple: a=v/t
But how do I convert it into using hz and g as input values and output the displacement (amplitude) ?

2. May 6, 2013

### VantagePoint72

The formula $v = at$ is only true for constant acceleration. In vibrational motion, you have a changing force (e.g. changing linearly in simple harmonic motion) and so acceleration is not constant. In the general case, you need to integrate: $v(t) = \int _{t_0} ^t a(t) dt + v(t_0)$. Since you will have a set of numerical data points, you will need to do use some technique of numerical integration. With $t_0$ chosen so that $v(t_0)=0$ (the turn-around point), you can ignore the integration constant. Since integrating over a single cycle of periodic motion is equivalent wherever in the cycle you start, you can just thrown this term out. Now, having $v$ as a function of time, you can numerically integrate it again to get displacement as a function of time.

3. May 6, 2013

### Pivskid

Can you outline for me more specifically how it would look like ?
I assume it is the "length" (time) from turn-around-point(=t0) to peak acceleration, that is necessary.
That time should be =(1/hz)/4
How to get from there ? ...integrate from t0 to (1/hz)/4 ?

4. May 6, 2013

### AlephZero

The usual practical way to do this is to assume the vibration is simple harmonic motion, not to integrate the acceleration numerically.

If the displacement is $A\sin \omega t$, the acceleration is $-\omega^2 A \sin \omega t$.

You need to convert "g" into units of meters/sec^2 or inches/sec^2, and $\omega = 2 \pi f$ .

5. May 6, 2013

### Pivskid

If ω=2π*frequency
Then what is t ? 1/f ?

6. May 6, 2013

### AlephZero

t is time.

If the maximum displacement is $A$, and the maximum acceleration is $\omega^2 A$.

7. May 6, 2013

### VantagePoint72

Given this a physical experiment that is actually measuring the acceleration, that would rather seem to defeat the purpose. Why bother even using an accelerometer if you're just going to assume SHM?

8. May 6, 2013

### AlephZero

First, most "vibrating machinery" does move in SHM, to a good approximation. It is responding to a steady-state sinusoidal applied force.

Second, mesuring relatively large accelerations is cheap and simple, compared with measuring small displacements directly.

For example at 50 Hz, or 3000 RPM for a rotating machine, a peak accleration of 1g corresponds to a peak displacement of about 0.01mm.

9. May 6, 2013

### VantagePoint72

Yes, I am perfectly aware of that. However, the OP asked how to calculate displacement given a series of acceleration measurements. Making assumptions—no matter how well founded—about the form the displacement will take defeats the entire purpose of doing an empirical measurement. This would be like if one asked, "How can I accurately measure the variation in the earth's gravitational field near its surface?" and you responded, "Well, we usually assume the gravitational field is constant near the earth's surface, which it is to good approximation. Just use $F=mg$." It's completely circular.

I don't understand what you're saying. I am not proposing measuring the peak displacement directly, I am explaining how to calculate it from the acceleration measurements, exactly as the OP requested. As you have said, measuring accelerations of the order OP will be looking at is easy and cheap, and hence provides a nice way of directly determining the (otherwise hard to measure) displacements.

Numerical integration is easy and, if the accelerometer has sufficient time resolution, very accurate. It will give OP his answer directly without having to make any assumptions about the linearity of material. In fact, it will provide a way of checking to what degree the motion deviates from SHM, if that's something OP is interested to know.

I can't for the life of me understand why someone would advocate relying on a theoretical (and imperfect) model for determining something you could actually directly determine using the available tools. It makes no sense at all.

Last edited: May 6, 2013
10. May 6, 2013

### VantagePoint72

Your acceleration measurements should be fairly cyclic (I say "fairly" because friction means they will decay over time). So, the easiest thing to do will be to pick out the peaks in the acceleration profile. These will also correspond to the peak displacement since the restorative force in most oscillatory motion is generally at a peak at the maximum displacement. Hence, if you numerically integrate (twice) from one crest to the next you will get the displacement between the two points of greatest displacement from equilibrium. That is, twice the amplitude of oscillation.

This looks to be a nice overview of the techniques of numerical integration for definite integrals. Software packages like MATLAB and Python have built in modules for doing this as well.

Last edited: May 6, 2013
11. May 6, 2013

### AlephZero

The OP also said
Which part of "sinusoidal" don't you understand?

The OP never said anything about "a series of acceleration measurements" either. You don't need digitallly sampled data to do this sort of vibration monitoring, in the real world.

12. May 6, 2013

### cjl

Because it is entirely possible (and I would guess even fairly likely) that the error introduced in the numerical integration (depending on the sample rate and resolution of the equipment) may in fact be much larger than the error introduced by assuming SHM and calculating based on the peak (or RMS) acceleration value.

13. May 7, 2013

### Pivskid

Hi again
I found the solution.
http://www.spaceagecontrol.com/calcsinm.htm

Very simple, and it conforms to my actual "real-life" measurements.

Now I just wonder how the sin() and integration appeared on the scene.