parametric-vibration

name	parametric-vibration
description	Parametric vibration — Mathieu equation stability (Strutt/Ince-Strutt diagram), Hill's equation (Floquet theory), principal parametric resonance at Ω ≈ 2ωₙ, combination resonances, parametric instability in rotating shafts (Duffing), columns under pulsating loads, cable vibration (Faraday waves), beam with time-varying stiffness, multiple-scales perturbation solution, Bolotin's regions, damping threshold for instability, and engineering applications (compressors, tension members, rotating machinery).
metadata	{"priority":7,"promptSignals":{"phrases":["parametric vibration","parametric excitation","Mathieu equation","parametric resonance","Strutt diagram","Hill equation"],"minScore":3}}

Parametric Vibration — Complete Skill

Governing Equation — Mathieu Equation

Standard Form

Mathieu equation: ẍ + (δ + ε cos τ) × x = 0 [τ = Ωt = normalized time; δ = (ωₙ/Ω/2)² = static frequency ratio; ε = parametric excitation amplitude]

Physical meaning: Equation of motion when the system's stiffness (or another parameter) varies periodically with time Occurs when: axial load on beam varies, shaft rotational speed fluctuates, cable tension oscillates

Physical form (beam with pulsating axial load): mẍ + c ẋ + [k₀ - P_dynamic cos(Ωt)] × x = 0 Convert: δ = k₀ / (mΩ²/4); ε = P_dynamic / (mΩ²/4) ωₙ = √(k₀/m); Ω_excitation = forcing frequency of parameter variation

Key difference from forced vibration: Forced: ẍ + 2ζωₙẋ + ωₙ²x = F₀ cos Ωt [resonance at Ω = ωₙ] Parametric: resonance at Ω ≈ 2ωₙ (principal) and Ω ≈ ωₙ, 2ωₙ/3, ... (higher order) Parametric resonance: amplitude grows exponentially; cannot be bounded by damping alone in unstable zones

Stability Analysis — Strutt (Ince-Strutt) Diagram

Stability Regions

Strutt diagram (δ-ε plane): Stability boundary: curves in (δ, ε) space separating stable and unstable regions Stable: bounded motion; Unstable: exponentially growing amplitude

Principal instability zone (first unstable region): Centered at δ = 1/4 (i.e., Ω = 2ωₙ) Width of instability zone at small ε: ε_boundary ≈ ±2√(δ) × ε/2 → approximately |ε| > 4ζ√δ (with damping)

Instability zones at: δ = n²/4, n = 1, 2, 3, ... (i.e., Ω = 2ωₙ/n for integer n) n = 1 (Ω = 2ωₙ): principal parametric resonance — widest zone; most dangerous n = 2 (Ω = ωₙ): secondary resonance — narrower zone n = 3 (Ω = 2ωₙ/3): tertiary resonance — very narrow; often neglected

Approximate boundaries (small ε, Bolotin): Principal zone boundary: δ = 1/4 ± ε/2 + O(ε²) In physical parameters (Ω ≈ 2ωₙ): instability when |P_dynamic/P_cr| > 4ζ (damping threshold)

Damping effect: Damped Mathieu equation: ẍ + 2ζωₙẋ + (δ + ε cosτ)x = 0 With damping, instability zones shrink; minimum ε for instability: ε_min = 4ζ × δ^(1/2) [at principal resonance δ = 1/4; ε_min ≈ 2ζ] For small damping: system can still be unstable even with ζ > 0 if ε is large enough

Floquet Theory (General Periodic Coefficients — Hill's Equation)

Hill's equation: ẍ + p(τ)x = 0 [p(τ) periodic with period T]

Floquet theorem: x(τ + T) = M × x(τ) [M = monodromy matrix; 2×2 for second-order equation] Stability: |trace(M)| ≤ 2 → stable; |trace(M)| > 2 → unstable

Characteristic multipliers (Floquet multipliers): ρ₁ × ρ₂ = 1 (for Mathieu; conservative); if |ρ| > 1 → unstable Characteristic exponents: ρ = e^(μT); μ = Floquet exponent (complex in general)

Computation: integrate two linearly independent solutions over one period numerically → construct M → compute trace

Multiple-Scales Perturbation Solution

Asymptotic Analysis (Small ε)

Near principal resonance (Ω ≈ 2ωₙ): Let Ω = 2ωₙ + εσ (detuning parameter σ) Two-timing: x(τ; ε) = x₀(T₀,T₁) + ε×x₁(T₀,T₁) + ... T₀ = τ, T₁ = ετ

Solvability condition at O(ε): a' = -ζωₙ a + (ε/4ωₙ) × a × sin(2γ) γ' = -(σ/2) + (ε/4ωₙ) × cos(2γ) [a = amplitude; γ = phase angle in slow time]

Steady-state amplitude: Steady state (a' = 0, γ' = 0): non-trivial solution when ε²/16ωₙ² > ζ² + σ²/ωₙ² → instability boundary At exact resonance (σ = 0): instability requires ε > 4ζωₙ [parametric excitation exceeds damping]

Frequency response: In stable zone near resonance: amplitude bounded; response at ωₙ (not Ω) In unstable zone: response grows at rate Im(μ) per period

Physical Applications

Column Under Pulsating Axial Load (Bolotin)

Equation (Euler beam, pulsating load): EI × ∂⁴w/∂z⁴ + [P₀ + P_t cos(Ωt)] × ∂²w/∂z² + ρA × ∂²w/∂t² = 0

Galerkin reduction (first mode, simply supported): Leads to Mathieu: q̈ + ω₁²[1 - (P₀/P_cr) - (P_t/P_cr)cos(Ωt)] × q = 0 ω₁ = π²/L² × √(EI/(ρA)); P_cr = π²EI/L²

Instability conditions (Bolotin): Principal zone: Ω ≈ 2ω₁√(1 - P₀/P_cr) Upper boundary of instability: P_t/P_cr = (Ω/ω₁)² - 4 + ... (function of detuning)

Design check: Compute P₀/P_cr (static ratio) and P_t/P_cr (dynamic ratio) Plot on Strutt diagram → if point falls in unstable region → increase damping or change excitation frequency

Rotating Shaft with Asymmetric Stiffness

Jeffcott rotor with asymmetric shaft stiffness: (k₁ ≠ k₂ in two principal directions) In rotating frame: parametric excitation at 2× rotational frequency (Ω_exc = 2ω_rot) Instability zone: ω_rot ≈ (k₁+k₂)/(2m×4) → not simply at ωₙ/2

Practical example: Two-pole generator rotor (stiffness differs due to pole geometry) Parametric resonance: instability at ω_rot ≈ ωₙ/2 (half critical speed region) Mitigation: symmetric solid rotor design; add damping; avoid operation near ωₙ/2

Cable Vibration (Faraday Waves)

Taut cable with oscillating tension: T(t) = T₀ + T₁ cos(Ωt); cable wave speed c = √(T/ρ_L) First natural frequency: f₁ = c/(2L) = (1/2L)×√(T₀/ρ_L)

Parametric resonance condition: Ω ≈ 2f₁ = (1/L)×√(T₀/ρ_L) [principal resonance] Application: stay cables, mooring lines under wave-induced tension variations

Mitigation: increase cable sag ratio (increases damping); add viscous dampers at cable anchors; install cross-ties

Compressor Connecting Rod

Axial force in connecting rod: F(t) = m_piston × r × ω² × (cosθ + λcos2θ + ...) [sinusoidal + harmonics due to kinematics] Connecting rod acts as column under pulsating compressive load Risk: parametric resonance if rod natural frequency ≈ 2× crankshaft speed

Crank speed check: f_exc = crankshaft RPM / 60; principal resonance if f_rod,bending ≈ 2 × f_exc

Damping Requirements

Minimum Damping to Prevent Instability

Threshold condition (principal resonance, Ω = 2ωₙ): ε_cr = 2ζ × (1 + σ²/ζ²)^(1/2) ≈ 2ζ [at exact resonance σ = 0]

Required damping ratio: ζ_min = ε/2 = (P_t/P_cr)/2 [for Bolotin column; must have ζ ≥ ζ_min to remain stable]

Engineering damping sources: Material damping: ζ ≈ 0.001–0.01 (metals); insufficient for high ε Structural joints: ζ ≈ 0.005–0.02 Fluid damping: ζ ≈ 0.002–0.05 (depends on viscosity and gap geometry) Dedicated dampers: ζ = 0.05–0.30 (tuned mass damper, viscous dampers)

Rule of thumb: If ε = P_t/P_cr > 0.1 → significant instability risk unless ζ > 0.05

Standards and References

Reference	Scope
Bolotin (1964) Dynamic Stability of Elastic Systems	Definitive treatment of parametric excitation in structural mechanics
Nayfeh & Mook (1979) Nonlinear Oscillations	Multiple scales for Mathieu/Hill equations
Mathieu (1868)	Original equation for elliptic membrane vibration
Floquet (1883)	General theory of linear periodic ODEs
Genta (2005) Dynamics of Rotating Systems	Shaft parametric instability
ASME B31.3	Piping design (includes cable vibration appendix)

Output

Provide: physical system description (what parameter varies, variation amplitude ε and frequency Ω [rad/s]), static frequency ratio δ = (2ωₙ/Ω)² and location on Strutt diagram (stable or unstable zone), principal instability zone bounds (ε_boundary vs. detuning σ/ωₙ), damping ratio ζ [%] (measured or assumed) and minimum damping threshold ζ_min = ε/2, Floquet multipliers ρ₁,ρ₂ (if numerical integration performed; |ρ|>1 → unstable), multiple-scales amplitude response (if within stable zone: steady-state amplitude a/X_static vs. Ω/ωₙ), Bolotin instability boundaries (P₀/P_cr vs. P_t/P_cr stability chart for pulsating column loads), response growth rate μ [1/s] in unstable zone (time to 10× amplitude), mitigation recommendations (detuning Ω away from 2ωₙ; damping addition; stiffness modification), and applicable reference (Bolotin 1964, Nayfeh & Mook 1979).