Dynamics: Damped Free Vibration

In an ideal world, a swinging pendulum would oscillate forever, and a car’s suspension would bounce indefinitely after hitting a pothole. Reality, of course, is far less tidy. Every real-world oscillatory system loses energy, causing its motion to gradually subside—a phenomenon known as damping. Understanding damped free vibration is crucial for engineers designing everything from earthquake-resistant skyscrapers and comfortable vehicle suspensions to precise micro-electromechanical sensors. It is the bridge between the perfect, perpetual motion of ideal theory and the controlled, decaying response of functional engineering systems.

The Viscous Damping Model and Equation of Motion

To analyze damping mathematically, we need a tractable model. The most common and analytically convenient model is viscous damping. This model assumes the damping force is directly proportional to the velocity of the mass and acts in the opposite direction. You can visualize it as the force exerted by a piston moving through a thick oil: the faster you try to push it, the greater the resistive force.

Consider a simple mass-spring-damper system with mass $m$, spring stiffness $k$, and a viscous damper with damping coefficient $c$. Applying Newton’s second law, the sum of the forces equals mass times acceleration. The forces are the spring force ($-kx$) and the damping force ($-c\dot{x}$). This yields the equation of motion for a damped free vibration:

$$m\ddot{x}+c\dot{x}+kx=0$$

This is a homogeneous, second-order linear differential equation with constant coefficients. Its solution describes how the system responds after being displaced and released, with no external force applied—hence "free" vibration. The character of this response depends entirely on the relationship between the damping coefficient $c$ and the system's inherent stiffness and mass.

Damping Ratio and Critical Damping

Instead of working directly with the damping coefficient $c$, it’s more insightful to use a dimensionless measure called the damping ratio, denoted by $\zeta$ (zeta). The damping ratio compares the actual damping $c$ to a special benchmark value known as the critical damping coefficient, $c_c$.

Critical damping $c_c$ is defined as the minimum amount of viscous damping that prevents oscillatory motion. It is calculated from the system's mass and stiffness:

$$c_c=2\sqrt{km}=2m{\omega }_n$$

Here, ${\omega }_n=\sqrt{k/m}$ is the system's undamped natural frequency in radians per second. The damping ratio is then:

$$\zeta =\frac{c}{c_c}=\frac{c}{2\sqrt{km}}=\frac{c}{2m{\omega }_n}$$

This single parameter $\zeta$ becomes the key that unlocks the three distinct types of system response.

Classifying the Response: Underdamped, Critically Damped, and Overdamped

The value of the damping ratio $\zeta$ categorizes the system's behavior into three regimes. Let's solve the characteristic equation associated with our differential equation: $m{s}^2+cs+k=0$. Dividing by $m$ and substituting ${\omega }_n$ and $\zeta$, we get $s^2+2\zeta {\omega }_{n}s+{\omega }_n^2=0$. The roots of this equation are:

$$s=-\zeta {\omega }_n\pm {\omega }_n\sqrt{{\zeta }^2-1}$$

The nature of these roots determines the system's motion.

1. Underdamped ($0<\zeta <1$): The most common case for structural and mechanical systems designed to oscillate. Here, the term under the square root is negative, leading to complex roots. The solution is an oscillatory motion that decays exponentially over time:

$$x(t)=e^{-\zeta {\omega }_{n}t}[A\cos ({\omega }_{d}t)+B\sin ({\omega }_{d}t)]$$ The frequency of this decaying oscillation is the damped natural frequency, ${\omega }_d={\omega }_n\sqrt{1-{\zeta }^2}$, which is always less than ${\omega }_n$.

1. Critically Damped ($\zeta =1$): The roots are real, equal, and negative. This represents the fastest possible return to equilibrium without oscillation. The solution takes the form:

$$x(t)=(A+Bt)e^{-{\omega }_{n}t}$$ Critical damping is often a design target for systems like door closers or electrical meters, where overshoot is undesirable.

1. Overdamped ($\zeta >1$): The roots are real, distinct, and negative. The system returns to equilibrium without oscillating, but does so more slowly than a critically damped system. The response is a sum of two decaying exponentials:

$$x(t)=A{e}^{s_{1}t}+B{e}^{s_{2}t}$$ where $s_1$ and $s_2$ are the two real, negative roots.

The Logarithmic Decrement: Measuring Damping Experimentally

For an underdamped system, we often need to determine the damping ratio $\zeta$ from experimental data, such as a recorded displacement-time plot. The logarithmic decrement $\delta$ provides a powerful method to do this. It is defined as the natural logarithm of the ratio of any two successive peak amplitudes in the free vibration decay.

$$\delta =\ln \left(\frac{x(t)}{x(t+T_d)}\right)$$

where $T_d=2\pi /{\omega }_d$ is the damped period. Substituting the underdamped solution, the relationship simplifies beautifully to:

$$\delta =\frac{2\pi \zeta }{\sqrt{1-{\zeta }^2}}$$

For lightly damped systems ($\zeta <0.2$), the approximation $\delta \approx 2\pi \zeta$ is often used. To improve accuracy, you can measure the amplitude ratio over $n$ cycles: $\delta =(1/n)\ln (x_0/x_n)$. This experimental technique is fundamental in modal testing of structures, machinery diagnostics, and validating analytical models against real-world behavior.

Common Pitfalls

• Confusing Damping Coefficient with Damping Ratio: A common error is to state that a system is "underdamped" because $c$ is small. This is not necessarily true. You must compute the dimensionless ratio $\zeta =c/c_c$. A small $c$ value on a very light, stiff system (resulting in a tiny $c_c$) could still produce a high $\zeta$ and an overdamped response.
• Misapplying the Damped Frequency Formula: Remember that the damped natural frequency ${\omega }_d={\omega }_n\sqrt{1-{\zeta }^2}$ is only valid for underdamped systems ($\zeta <1$). For critically damped and overdamped systems, there is no oscillatory frequency.
• Incorrect Logarithmic Decrement Calculation: When using $\delta =(1/n)\ln (x_0/x_n)$, ensure $x_0$ and $x_n$ are peak amplitudes measured n complete cycles apart. Using amplitudes from non-consecutive peaks without dividing by the number of cycles $n$ will yield an incorrect, inflated estimate of $\zeta$.
• Overlooking the Exponential Envelope: When sketching an underdamped response, students often draw a sinusoid with linearly decreasing amplitude. The correct amplitude decay is exponential, governed by the term $e^{-\zeta {\omega }_{n}t}$. The peaks should lie on the curve defined by $\pm X{e}^{-\zeta {\omega }_{n}t}$, not on a straight line.

Summary

• Viscous damping models resistance proportional to velocity, leading to the fundamental equation of motion: $m\ddot{x}+c\dot{x}+kx=0$.
• The damping ratio $\zeta =c/c_c$ is the key parameter, where critical damping $c_c=2\sqrt{km}$ is the threshold for oscillatory motion.
• System response is classified as underdamped (decaying oscillation, $\zeta <1$), critically damped (fastest non-oscillatory return, $\zeta =1$), or overdamped (slow non-oscillatory return, $\zeta >1$).
• The logarithmic decrement $\delta$ provides a practical method to measure the damping ratio $\zeta$ experimentally from the decay of free vibration peaks in an underdamped system.
• Mastery of these concepts allows engineers to predict, measure, and design the dynamic response of systems to meet specific performance criteria, such as quick settling time, vibration isolation, or energy dissipation.