Open Educational Resources

Review of single and multi-degree of freedom (mdof) systems: Center of Percussion and Vibration of a Compound Pendulum

The concept of the center of percussion is usually discussed using the analysis of the changes in angular momentum and linear momentum associated with a body impacting a rotating object and essentially determining the point at which there is no change in acceleration at the pivot point due to this impact. (The usual example is a ball hitting a baseball or cricket bat.) For the compound pendulum shown the natural frequency is p = (mgd/J_0)^{1/2} where J_0 is the moment of inertia of the body about the axis of rotation at O. Instead of having a distributed mass, let the mass of the pendulum be concentrated at a single point located at a distance L from the pivot point O. The natural frequency of this simple pendulum is (g/L)^{1/2}. If the natural frequencies of these two pendulums are the same, then the point at which m is concentrated is called the center of oscillation and is identical to the center of percussion. If instead, we consider rotation about the center of oscillation (percussion) then the original pivot at O becomes the center of percussion (oscillation).

Early automobile manufacturers did not realize this and as a result, cars had much more uncomfortable rides because impacts at one set of wheels (front or rear) were felt at the other. One of the first automobile designs to incorporate this idea was the very futuristic Chysler Airflow (circa 1930) which was years ahead in its design in a number of aspects. (So far ahead that it was a commercial failure.) For the illustration of the automobile shown the relationship between the position of the center of gravity and the positions of the front and rear wheels is given by r^* = (ab)^{1/2} where r^* is the radius of gyration of the vehicle about its center of gravity and I = a + b is the wheelbase (a and b are the horizontal distances from the front and rear wheels to the center of gravity). From a vibration perspective if a car is designed with a mass distribution to meet these criteria (almost all recent automobiles are) then we can examine what happens at one set of wheels individually for certain situations (e.g. hitting a bump with the front wheels only).

Using the diagram below show that the condition for the center of percussion to be at the front wheels with the pivot at the rear wheels that r^* = (ab)^1/2 and that this is also true for the pivot point moving to the front wheels with the impact at the rear ones.

    \[J = m(r^*)^2\]

Where J is the moment of inertia about G and r^* is the radius of gyration about G:

    \[\begin{split}+ \circlearrowleft  \sum M_0 &= -mgd\theta \\&= J \ddot{\theta}+md\ddot{\theta}(d)\end{split}\]

Therefore:

    \[(J+md^2)\ddot{\theta} + mgd\theta = 0\]

    \[J_0 = J + md^2\]

Therefore:

    \[\begin{split}p &= \sqrt{\frac{mgd}{m(r^*)^2+md^2}} \\&= \sqrt{\frac{gd}{(r^*)^2+d^2}}\end{split}\]

Equivalent Pendulum

    \[p = \sqrt{\frac{g}{L}}\]

Therefore choose:

    \[L = \frac{(r^*)^2+d^2}{d}\]

    \[ \ell = a+b \]

For application to the auto above \ell = L (wheelbase). For rotation about B:

    \[d = b\]

Therefore:

    \[\begin{split} L &= a+b \\&= \frac{(r^*)^2+b^2}{b} \\&= \frac{(r^*)^2}{b} +b \end{split}\]

Therefore:

    \[r^* = \sqrt{ab}\]

For rotation about A:

    \[d = a\]

Therefore:

    \[a+b = \frac{(r^*)^2}{a} +a \]

    \[r^* = \sqrt{ab}\]

Determination of Damping Ratio

We know that for any SDOF viscously damped system the equation for free vibration is:

    \[m_e\ddot{x} + C_e\dot{x} + k_ex = 0\]

and the system parameters are the natural frequency:

    \[p = \sqrt{\frac{k_e}{m_e}}\]

and the damping ratio:

    \[\zeta = \frac{C_e}{C_c} = \frac{Ce}{Z_{mep}}\]

We can determine the damping ratio from how rapidly the vibration decays. To do this we use the so-called logarithmic decrement. Consider the damped response:

    \[ x(t) = e^{-\zeta pt}\big[X\sin{(\sqrt{1-\zeta^2}pt+\phi)}\big]\]

which for \zeta < 1 looks like:

These peaks occur approximately when:

    \[\sin{\big(\sqrt{1 - \zeta^2}pt + \phi\big)} = 1\]

Therefore if the first occurs at t_0 then the next time it is at:

    \[t_1 = t_0 + \frac{2\pi}{p\sqrt{1-\zeta^2}}\]

Therefore:

    \[\begin{split} \frac{X_0}{X_1} &= \frac{e\strut^{-\zeta pt_0}X}{e\strut^{-\zeta pt_1}X}  \\&= \frac{e\strut^{-\zeta pt_0}} {e\strut^{-\zeta p\big(t_0 + \frac{2\pi}{p \sqrt{1-\zeta^2}}\big)}} \\&= e\strut^{\frac{2\pi\zeta}{\sqrt{1-\zeta^2}}} \end{split}\]

and this is the same ratio for \frac{X_1}{X_2}, \frac{X_2}{X_3}, etc.

Therefore this creates the logarithmic decrement equation:

    \[ \begin{split} \ln{\frac{X_0}{X_1}} &= \frac{2\pi\zeta}{\sqrt{1-\zeta^2}} \\&\equiv \delta \end{split}\]

If we take n cycles:

    \[\frac{X_0}{X_n} = \frac{X_0}{X_1}\frac{X_1}{X_2}\frac{X_2}{X_3}..... \frac{X_{n-1}}{X_n}\]

    \[\begin{split}\ln {\frac{X_0}{X_n}} &= \ln {\frac{X_0}{X_1}} + \ln {\frac{X_1}{X_2}}..... \\&= n\frac{2\pi\zeta}{\sqrt{1-\zeta^2}} \\&= n\delta \end{split}\]

Therefore:

    \[\begin{split} \delta &= \frac{1}{n}\ln{\frac{X_0}{X_n}} \\&= \frac{2\pi\zeta}{\sqrt{1-\zeta^2}} \end{split}\]

Therefore if we measure \frac{X_0}{X_n} we can easily determine the damping ratio \zeta.

If \zeta<<1, \delta \approx2\pi\zeta

\zeta0.050.10.20.30.40.5
\frac{X_0}{X_1}1.371.873.617.2115.537.6

It does not take a large \zeta to quickly damp out vibrations.

Leave a Reply

Your email address will not be published.