Explore parametric resonance with the Mathieu equation. Adjust modulation parameters to see the stability diagram, time-domain response, and phase portrait. Discover why pumping a swing at twice the natural frequency causes runaway oscillation.
Parameters
δ (natural-freq² proxy)
ε (modulation depth)
Damping γ
Initial displacement x₀
Simulation time
T
Stable
Results
Floquet |μ₁|
—
Floquet |μ₂|
—
Growth rate
—
Resonance order
—
Stability
Red: unstable region · Gray: stable region · Yellow cross: current (δ, ε)
Primary resonance occurs near $\delta \approx n^2$ for $n = 1, 2, 3, \ldots$ Damping shrinks the unstable wedges.
What is a Parametric Oscillator?
🙋
What exactly is a "parametric" oscillator? I thought oscillators just had a spring constant and mass.
🎓
Great question! In a normal oscillator, the restoring force is constant. A parametric oscillator is different because one of its key parameters—like the stiffness of the spring—is varied in time. In this simulator, that's the "Modulation amplitude" (q) and "Modulation frequency" (ωp). Try setting q to zero in the controls above; you'll see it's just a normal damped oscillator. Now increase q—you're literally pumping energy in by wiggling the spring's stiffness.
🙋
Wait, really? So you can make it oscillate wildly just by changing the stiffness, not by pushing it directly? That seems weird.
🎓
Exactly! It's counter-intuitive but true. The most famous example is a playground swing. You don't push it directly back and forth; you pump your legs, which effectively changes the pendulum length (a parameter) at the right moment. In our simulator, the "principal resonance" happens when the modulation frequency ωp is about twice the natural frequency ω₀. Try setting ωp to be exactly 2 * ω₀ and watch the amplitude grow even with a tiny q.
🙋
Okay, so what's the "stability diagram" for? It looks like a map with shaded and white regions.
🎓
That's the key to engineering with these systems! The diagram plots the stability parameter 'a' against the modulation strength 'q'. If your (a, q) point is in a white "unstable" region, the system experiences parametric resonance and the oscillations grow without bound. If you're in a shaded "stable" region, the oscillations die out. Now, slide the damping ratio ζ up and down. See how the unstable regions shrink? Damping is your friend for suppressing this instability.
Physical Model & Key Equations
The core physics is described by the Mathieu equation, a special case of a differential equation with periodic coefficients. It governs the motion of a damped oscillator whose natural frequency is being modulated sinusoidally in time.
Variables:
• $x(t)$: Displacement of the oscillator.
• $\omega_0$: Natural frequency of the unmodulated system (set by the "Natural frequency" slider).
• $\zeta$: Damping ratio (controls how quickly free oscillations decay).
• $q$: Modulation amplitude (strength of the stiffness variation).
• $\omega_p$: Modulation frequency (how fast the stiffness is varied).
To analyze stability, we often transform the equation. A common form uses dimensionless parameters $a$ and $q$, which define the famous Mathieu stability chart you see in the simulator.
$$a = 4\left(\frac{\omega_0}{\omega_p}\right)^2$$
Physical Meaning: The parameter $a$ is essentially a scaled ratio of the natural frequency to the pumping frequency. The stability of the system (whether $x(t)$ grows or decays) depends entirely on where the point $(a, q)$ lies on this chart. The principal parametric resonance occurs near $a \approx 1$, which corresponds to $\omega_p \approx 2\omega_0$.
Frequently Asked Questions
Adjust δ and ε so that they fall within the red region (unstable region) on the stability chart. For example, setting ε to around 2.0 and δ to around 1.0 will cause typical parametric resonance, and you can observe the amplitude rapidly increasing in the time history response.
Increasing the damping γ reduces the unstable region and suppresses the occurrence of resonance. When γ is 0, the unstable region appears wide, but if you increase γ to about 0.5, the system may become stable even with the same δ and ε. Observe the changes in the chart in real time.
It models phenomena where the stiffness or capacitance of a system changes periodically, such as the motion of pumping a swing, parametric amplification in electromagnetic circuits, and the dynamic stability of bridges. The Mathieu equation is the fundamental equation for analyzing these parametric resonances.
When parametric resonance occurs, the system absorbs energy from the outside, causing the amplitude to increase exponentially. This is more pronounced when damping is small. In reality, nonlinear effects and damping suppress the divergence, but this simulator assumes a linear model.
Real-World Applications
Rotating Machinery & Driveshafts: A long, spinning shaft can have a slight asymmetry or varying cross-section. As it rotates, its effective bending stiffness changes periodically, creating a parametric excitation. If the rotation speed matches twice a natural bending frequency, catastrophic vibrations (parametric resonance) can occur, leading to failure.
Offshore Platform Legs: The legs of an offshore oil platform are submerged in water. As waves pass, the changing water level periodically submerges and exposes parts of the leg, altering its effective stiffness and mass. Engineers must check that wave frequencies don't trigger parametric instability in the structure.
Microelectromechanical Systems (MEMS): Tiny mechanical resonators in sensors often use parametric pumping to achieve ultra-high gain and sensitivity. By modulating the electrical voltage on an electrode (which changes electrostatic stiffness), designers can create very precise parametric amplifiers that outperform direct forcing methods.
Parametric Amplifiers in Electronics: At microwave and radio frequencies, direct amplification adds noise. A parametric amplifier uses a nonlinear reactance (like a varactor diode) whose capacitance is "pumped" at a high frequency. A weak input signal can be amplified with remarkably low added noise, crucial for deep-space communication and radio astronomy.
Common Misconceptions and Points to Note
When you start using this simulator, there are a few common pitfalls. First, you might think "larger ε always leads to instability," but it's not that simple. Near the tip of the δ=1 wedge, the system is indeed unstable for small ε, but as you increase ε further, "islands" of stability can appear. For example, try slowly varying ε from 0 to 2 with δ=1 and γ=0. You should see the state alternate between unstable → stable → unstable. This happens because excessive modulation can actually desynchronize the timing of energy injection.
Next, remember that simulation "divergence" is not equivalent to real-world "failure". Since this tool uses a linear model, divergence indicates theoretically unbounded growth. However, real structures always have nonlinearities (e.g., a spring cannot extend beyond its limit), so the amplitude will settle at a finite value or lead to a different failure mode. Don't immediately judge a red zone in the simulation as unacceptable; instead, view it as a danger zone where you need to carefully examine the nonlinear response.
Finally, a tip on parameter settings: keep in mind that "simulating with γ (damping) near zero is often unrealistic". Most mechanical structures have some inherent damping. While γ=0 is fine for understanding principles in research, for practical applications, get into the habit of setting a realistic damping coefficient (e.g., around 0.01 to 0.05 for steel structures), estimated from materials and joints, before looking at the chart. Otherwise, you might end up with an overly conservative (heavy, costly) design.
Set the detuning parameter (delta) using vDeltaNum slider: positive values move away from the primary resonance frequency ω₀.
Adjust the modulation depth (epsilon) via vEpsNum: higher values increase parametric forcing amplitude; typical range 0.1–2.0 for industrial applications.
Configure damping ratio (vDampNum) between 0–0.3; structural systems rarely exceed 0.2 critical damping.
Enter initial displacement (vX0Num) in normalized units; observe whether the response enters the instability region of the Mathieu stability diagram.
Run simulation and examine the phase portrait and time-domain trace to confirm parametric resonance or stable oscillation.
Worked Example
A rotating shaft with 50 mm eccentricity modulates its stiffness at 2× operating speed (ω_p = 2ω₀). Set delta = 0, epsilon = 0.35, damping = 0.05, x₀ = 0.1 mm. The simulator shows exponential growth in the phase portrait until reaching amplitude ≈ 8.5 mm around t = 12 seconds; this matches the Mathieu instability region (a,q) ≈ (0, 0.175). Adding damping = 0.12 stabilizes motion to ±2.3 mm, confirming that critical damping prevents parametric failure in flexible rotordynamics.
Practical Notes
Parametric resonance occurs near δ = 0 with sufficient ε; avoid operating near ω = ω₀/2 in machinery with modulated stiffness (belts, gears, cam-followers).
The stability diagram (Ince–Strutt chart) reveals forbidden and allowed zones; simulator highlights these in real time as you adjust delta and epsilon.
High-frequency PWM-driven systems (power electronics, solenoids) exhibit parametric behavior at modulation sidebands; use the phase portrait to detect nonlinear energy transfer.