SDOF Equation of Motion
$$\ddot{x}+ 2\zeta\omega_n\dot{x}+ \omega_n^2 x = 0$$Logarithmic decrement:
$$\delta = \frac{2\pi\zeta}{\sqrt{1-\zeta^2}}, \quad \omega_n = \sqrt{k/m}$$Adjust mass, stiffness and damping ratio to compare underdamped, critically damped and overdamped responses in real time. Compute logarithmic decrement δ and 2%/5% settling times interactively.
Logarithmic decrement:
$$\delta = \frac{2\pi\zeta}{\sqrt{1-\zeta^2}}, \quad \omega_n = \sqrt{k/m}$$The fundamental equation governing a Single Degree-of-Freedom (SDOF) damped oscillator is Newton's second law, balancing inertia with the spring and damping forces.
$$\ddot{x}+ 2\zeta\omega_n\dot{x}+ \omega_n^2 x = 0$$Here, x(t) is the displacement, ẋ is velocity, and ẍ is acceleration. ωₙ = √(k/m) is the natural frequency (rad/s), and ζ (zeta) is the dimensionless damping ratio. The term 2ζωₙ encapsulates the damping force's magnitude.
For underdamped systems (ζ < 1), the response oscillates and decays. The rate of decay is quantified by the Logarithmic Decrement, which relates the ratio of consecutive oscillation peaks to the damping ratio.
$$\delta = \ln\left(\frac{x_n}{x_{n+1}}\right) = \frac{2\pi\zeta}{\sqrt{1-\zeta^2}}$$δ is the logarithmic decrement. xₙ and xₙ₊₁ are the amplitudes of two successive peaks. This equation is crucial because it allows engineers to determine the system's damping (ζ) simply by measuring the vibration decay, without needing to know m or k directly.
Automotive Suspension Tuning: Engineers carefully select damping ratios (typically between 0.2 and 0.4) for car shock absorbers. The goal is to be underdamped enough to absorb road bumps comfortably but damped enough to prevent excessive bouncing after a pothole, which you can simulate by adjusting ζ around 0.3.
Seismic Design of Buildings: During an earthquake, buildings act as giant vibrating structures. Added damping systems (like tuned mass dampers or viscous wall dampers) increase ζ to rapidly dissipate the seismic energy, reducing sway and preventing structural damage. This is the difference between an overdamped and underdamped collapse scenario.
Aerospace Component Testing: Satellite solar panels and antennae must withstand launch vibrations. Engineers use the log decrement method on test data to verify that the built-in damping is sufficient to prevent resonant vibrations from tearing the components apart in flight.
Consumer Electronics Design: The vibration of a spinning hard drive head or the haptic feedback in a smartphone is carefully damped. Too little damping (low ζ) causes ringing and slow settling, corrupting data or feeling "mushy." Too much damping (high ζ) makes the response sluggish.
When you start using this simulator, there are a few points that are easy to misunderstand. The first one is the tendency to think "the larger the damping ratio ζ, the faster the vibration always settles." While this is true up to ζ=1.0 (critical damping), increasing ζ beyond that (overdamping) actually increases the "settling time" it takes to return to equilibrium. For example, setting ζ=2.0 means there's no oscillation, but the system just lumbers back slowly. The fastest convergence happens at ζ=1.0, which is why it's called "critical."
The second point is overlooking the interaction when tweaking parameters individually. The damping ratio ζ is calculated from three parameters: mass m, damping coefficient c, and stiffness k (ζ = c / (2√mk)). So, you might think, "Doubling stiffness k makes vibration faster, so I should also double damping c." Actually, this leaves the value of ζ unchanged (try plugging it into the formula!). While the vibration speed (ω_n) changes, altering the response, the relative "braking strength" remains the same. Keep this formula in mind when you change parameters.
The third is a practical pitfall. The "settling time" in the simulation and the "time until vibration subsides" measured on an actual machine do not always match. The simulation calculates based on ideal initial conditions and a mathematical definition (e.g., the time for amplitude to fall within 2% of the initial value). But in the field, there's measurement noise, tiny sustained vibrations, and timing offsets in when measurement begins. It's important not to take simulation results at face value but to use them as an indicator, thinking "the theoretical value is around this."
The concept of this SDOF damped vibration is actually foundational to a surprisingly wide range of fields. First up is Control Engineering. Here, the response of an SDOF system is modeled as a "second-order lag system," forming the basic shape for evaluating the stability and responsiveness of feedback control systems. The concepts of this vibration (overshoot) and damping directly come into play even when tuning PID controller parameters.
Next is Acoustical Engineering. The vibration of speaker diaphragms or the strings/membranes of musical instruments is precisely damped vibration. Specifically, the length of a sound's "reverberation" or the "thud" of an impact sound is determined by the damping ratio ζ and the natural frequency ω_n. In high-end audio, the sound quality hinges on how quickly unwanted vibrations of the cabinet (the box) are damped (i.e., using materials with a high damping ratio).
Another, perhaps surprising, field deeply related is Electrical Circuits. The differential equation describing the transient response of an RLC series circuit (Resistor, Inductor, Capacitor) has exactly the same mathematical form as the equation of motion for an SDOF vibration system. Mass m corresponds to inductance L, damping coefficient c to resistance R, and stiffness k to the reciprocal of capacitance 1/C. Therefore, "underdamping," "critical damping," and "overdamping" in vibration appear directly as "oscillatory," "critical," and "non-oscillatory" transient responses in circuits. It's fascinating how completely different physical phenomena can be described by the same equations.
Once you're comfortable with this simulator, your next step could be to explore the world of "Multi-Degree-of-Freedom Vibration Systems (MDOF)." Real-world structures are modeled as continuous bodies with infinite degrees of freedom or as multi-degree-of-freedom systems with multiple coupled masses. Here, each vibration mode (e.g., a building's 1st mode, 2nd mode...) will have its own natural frequency and damping ratio. What you learned with SDOF forms the basis for "modal decomposition," which looks at each of these mode vibrations individually.
If you want to deepen the mathematical background, look into solving differential equations, paying particular attention to the nature of the roots (solutions) of the characteristic equation. The form of the solution to the equation of motion $m\ddot{x}+ c\dot{x}+ kx = 0$ changes depending on whether the roots $s$ of the characteristic equation are real or complex ($s = -\zeta\omega_n \pm \omega_n\sqrt{\zeta^2-1}$). The real part of these "roots" corresponds to the decay rate, and the imaginary part to the frequency. This understanding will definitely help when you learn about more complex systems or active control systems (state-space representation).
Finally, if you want to learn about simulation technology itself, investigate numerical analysis methods, especially time integration methods (Euler method, Runge-Kutta methods, Newmark-β method, etc.). You'll understand how this simulator solves the differential equations and draws the graphs behind the scenes. In practical CAE, these advanced numerical solvers are essential for handling nonlinear damping or large deformations. The best shortcut to a deep understanding is to start by writing your own simple program using the most basic Euler method to calculate SDOF vibration.