Harmonic Motion Simulator (SHM, Damped, Driven)

The fundamental equation governing all these motions is Newton's second law for a mass on a spring, including damping and driving forces. This is the Driven, Damped Harmonic Oscillator equation:

Here, $m$ is mass, $c$ is the damping constant (related to γ by $\gamma = c/(2m)$), $k$ is the spring constant (related to ω₀ by $\omega_0 = \sqrt{k/m}$), and $F_0$ is the amplitude of the driving force.

For the steady-state solution (after initial transients die out), the amplitude of oscillation $X$ is given by the formula central to this simulator:

This is the Amplitude Response Function. $X$ is the maximum displacement, $\omega_0$ is the natural frequency, $\omega$ is the driving frequency, and $\gamma$ is the damping coefficient. The denominator gets smallest when $\omega \approx \omega_0$, causing the resonant amplitude peak. The width of this peak is directly proportional to $\gamma$ and inversely proportional to the Q-factor.

Tuning Circuits in Radios: An LC circuit is an electrical harmonic oscillator. By adjusting the circuit's natural frequency (tuning) to match the frequency of a specific radio station (the driving signal), the circuit resonates, picking up that signal strongly while rejecting others. The Q-factor of the circuit determines how selective the radio is between closely spaced stations.

Seismic Engineering & Building Design: Buildings have natural sway frequencies. Engineers must design structures so their natural frequencies do not match the dominant frequencies of earthquakes or strong winds to avoid destructive resonance. Damping systems (like tuned mass dampers in skyscrapers) are added to increase effective γ and lower the Q, dissipating vibrational energy.

Atomic Force Microscopy (AFM): A sharp tip on a cantilever beam is oscillated near its resonance. When the tip interacts with a sample surface, the effective damping (γ) and resonant frequency (ω₀) shift. Monitoring these changes allows the microscope to map surface topography with atomic-scale resolution, a direct application of driven harmonic oscillator physics.

Automotive Suspension Design: A car's suspension is a damped oscillator (spring and shock absorber). The damping coefficient γ is carefully chosen for a low Q-factor. This ensures the car doesn't resonate uncomfortably on bumpy roads—the suspension absorbs and quickly dissipates the energy from bumps, providing a smooth ride.

First, a common mistake is reversing the relationship between the "damping coefficient γ" and the "Q factor". γ is the "coefficient that reduces vibration," so the larger it is, the faster the vibration settles. On the other hand, the Q factor is essentially its inverse; the smaller γ is (less damping), the larger the Q factor becomes. For example, decreasing γ from 0.1 to 0.01 makes the Q factor about 5 times larger, and the resonance peak becomes sharper and higher. In practical work, when you hear "a system with a high Q factor," you can interpret it as "a system with low damping that reacts sensitively to a specific frequency."

Next, a point where people often mistakenly think the "natural angular frequency ω₀" does not change even if parameters are altered. In this simulator, ω₀ is a value determined by "mass m" and "spring constant k." Therefore, changing the forced vibration frequency ω_d does not change ω₀ itself. However, in the real world, large amplitudes can cause material deformation changing k (nonlinearity), or changes in mass distribution altering the effective m. In other words, the resonance point itself can move. Keep in mind that the simulator uses a linear model.

Finally, avoid trying to observe the "steady state" of forced vibration immediately. Right after the external force is applied, it's in a "transient state" with unstable vibration. Especially in systems with low damping (high Q factor), it can take a very long time to settle into a steady state. For instance, starting a simulation with γ=0.05 and ω_d=ω₀, it might take dozens of cycles for the amplitude to reach its maximum value. When comparing with experimental data, always look at the waveform from the "steady state" after sufficient time has passed.