Skip to main content
Sandbox Physics

M068 · Vibration design

Parametric Resonance Playground

A stiffness actuator periodically changes the spring. Adjust modulation depth, frequency and phase, release the cart and audit the signed pump work. Integrate two independent basis releases over one period to measure Floquet multipliers; build and select a stability map.

Interactive modelParametric Resonance Playground
Acquired time0 s0\,\mathrm{s}
Primary displacement0 m0\,\mathrm{m}
Primary velocity0 m s−10\,\mathrm{m\,s^{-1}}
Mechanical energy0 J0\,\mathrm{J}
Integrated signed actuator work0 J0\,\mathrm{J}
Independent viscous dissipation0 J0\,\mathrm{J}
Independent energy ledger defect0 J0\,\mathrm{J}
Position change with half step · first four cycles0 m0\,\mathrm{m}
Actual integration step0 ms0\,\mathrm{ms}
Model domain / record statusReady\text{Ready}
Current stiffness00
Largest Floquet multiplier magnitude00
Asymptotic linear growth rate00
Numerical linear stability00
Monodromy determinant defect00
Multiplier change with half step00
Computed map points00
Map acquisition status00

Physics tutorial

A pump changes the rule of the spring

BackgroundA time-dependent stiffness can transfer energy even when no independent force drives the mass. The phase of that change matters to a finite release.

Why it mattersNIST DLMF 28.29 defines the monodromy matrix and its Floquet eigenvalues for Hill equations. The damped version here adds a determinant identity and independent work integration.

Start with the essentials

Focus question
How can one-period measurements predict long-term linear growth?
One-sentence intuition
Evolve two basis states to reconstruct the period map; compare its eigenvalue magnitudes with unity.

Core mathematical model

Modulated spring law

mx¨+2mζω0x˙+k0[1+hcos⁡(Ωt+ϕ)]x=0,ω0=k0/mm\ddot x+2m\zeta\omega_0\dot x+k_0[1+h\cos(\Omega t+\phi)]x=0,\quad \omega_0=\sqrt{k_0/m}

Depth stays below unity, so the instantaneous stiffness remains positive.

Actuator power

Ppump=12k˙(t)x2,E=12mx˙2+12k(t)x2P_{\mathrm{pump}}=\tfrac12\dot k(t)x^2,\quad E=\tfrac12m\dot x^2+\tfrac12k(t)x^2

Pump work can have either sign. Integrate it independently from viscous heat.

One-period map

[x(T)x˙(T)/ω0]=Φ(T)[x(0)x˙(0)/ω0],T=2π/Ω\begin{bmatrix}x(T)\\\dot x(T)/\omega_0\end{bmatrix}=\Phi(T)\begin{bmatrix}x(0)\\\dot x(0)/\omega_0\end{bmatrix},\quad T=2\pi/\Omega

Two independent normalized basis releases construct the matrix columns.

Growth and determinant

γ=T−1ln⁡ρ(Φ),det⁡Φ=e−2ζω0T\gamma=T^{-1}\ln\rho(\Phi),\quad \det\Phi=e^{-2\zeta\omega_0T}

A positive growth rate is a linear instability; the determinant supplies an independent numerical check.

Small-depth principal threshold

Ω≃2ω0,hcrit≃4ζ\Omega\simeq2\omega_0,\quad h_{\mathrm{crit}}\simeq4\zeta

This is a weak-modulation, weak-damping asymptotic relation, not the map calculation.

Common difficulties

Zero is still a solution

Typical misconceptionAn unstable multiplier means every release grows.

Better mental modelThe exact zero state remains zero; specially aligned states can initially shrink.

Finite phase sensitivity

Typical misconceptionChanging the starting phase changes the stability tongue.

Better mental modelIt changes startup projection; monodromy spectra are invariant under a shift of phase.

Linear growth does not saturate

Typical misconceptionThe domain stop is a physical saturation.

Better mental modelThe model has no nonlinear saturation; its growth continues beyond the displayed record.

Run the experiment

  1. 01

    Find a tongue

    Acquire the principal-tongue and detuned releases.

    What to observe: A positive Floquet rate aligns with long-term growth.
  2. 02

    Balance pumping and damping

    Increase damping and inspect work and heat.

    What to observe: The net energy change follows the independent ledger.
  3. 03

    Test phase and zero

    Change modulation phase, then select exact zero release.

    What to observe: Finite amplification changes while multiplier magnitude persists; the zero state remains motionless.
  4. 04

    Build the map

    Generate the map, select a cell, then halve the time step.

    What to observe: The selected period calculation has determinant and refinement diagnostics.