Shaping Energy Exchange with Gyroscopic Interconnections: a geometric approach

Gyroscopic interconnections enable redistribution of energy among degrees of freedom while preserving passivity and total energy, and they play a central role in controlled Lagrangian methods and IDA-PBC. Yet their quantitative effect on transient energy exchange and subsystem performance is not well characterised. We study a conservative mechanical system with constant skew-symmetric velocity coupling. Its dynamics are integrable and evolve on invariant two-tori, whose projections onto subsystem phase planes provide geometric description of energy exchange. When the ratio of normal-mode frequencies is rational, these projections become closed resonant Lissajous curves, enabling structured analysis of subsystem trajectories. To quantify subsystem behaviour, we introduce the inscribed-radius metric: the radius of the largest origin-centred circle contained in a projected trajectory. This gives a lower bound on attainable subsystem energy and acts as an internal performance measure. We derive resonance conditions and develop an efficient method to compute or certify the inscribed radius without time-domain simulation. Our results show that low-order resonances can strongly restrict energy depletion through phase-locking, whereas high-order resonances recover conservative bounds. These insights lead to an explicit interconnection-shaping design framework for both energy absorption and containment control strategies, while taking responsiveness into account.

💡 Research Summary

The paper investigates how gyroscopic (skew‑symmetric velocity) interconnections shape energy exchange between subsystems in a simple two‑degree‑of‑freedom conservative mechanical system. Starting from the linear equations ¨q + n ż + q = 0, ¨z − n q̇ + z = 0, the authors rewrite the dynamics in Hamiltonian form with the standard quadratic Hamiltonian H = ½(‖x‖²+‖p‖²) and a constant skew‑symmetric matrix J. Because Jᵀ = −J, the total energy is conserved, while the gyroscopic coupling redistributes energy between the q‑ and z‑subsystems without injection or dissipation.

By introducing the complex variable u = q + i z, the system decouples into two modal frequencies Ω₁ and Ω₂ that satisfy Ω₁Ω₂ = 1. The frequencies depend on the coupling parameter n, and the ratio Ω₁/Ω₂ is rational if and only if there exist coprime integers τ, σ such that Ω₁Ω₂ = τ/σ. This condition is equivalent to the algebraic relation n² = (τ − σ)²/(τσ). When the ratio is rational the flow on the invariant two‑torus is periodic, and the projections onto the phase planes (q, q̇) and (z, ż) become closed Lissajous curves.

The authors derive the exact convex hull of the projected trajectory as the Minkowski sum of two ellipses E_{Ω₁} and E_{Ω₂}. Using support functions they obtain a closed‑form boundary expression (equation 11). For n = 0 the hull reduces to a single ellipse; for n ≠ 0 it is the sum of two unequal ellipses and generally non‑elliptical.

A central contribution is the definition of the “resonant inscribed radius” r_res, the radius of the largest origin‑centred circle that fits inside the projected Lissajous curve. Because the subsystem energy is H_q = ½(q²+q̇²), r_res²/2 provides a guaranteed lower bound on the energy that the q‑subsystem can never drop below, regardless of initial impulse magnitude.

The paper distinguishes low‑order resonances (τ + σ ≤ M for some modest M) from high‑order resonances (τ + σ → ∞, τσ → 1). Theorem 1 shows that r_res = 0 (complete energy depletion) occurs precisely when the difference δ = τ − σ satisfies δ ≡ 2 (mod 4). The proof uses product‑to‑sum identities and 2‑adic valuations, revealing that at such resonances the two modal ellipses become exact negatives of each other at a common phase, cancelling both displacement and velocity.

When δ ≠ 2 (mod 4), Theorem 2 provides a systematic way to locate the global minimum of the radius. The Lissajous curve is partitioned into “lobes” bounded by successive zeros of q(θ) and q̇(θ). Lobes in which |cos β(θ)| (with β = (τ − σ)θ/2) attains values arbitrarily close to zero are identified as “envelope‑minimising”. Within each such lobe the squared radius R(θ) = q²+q̇² is strictly decreasing up to the unique interior point where q̇ = 0, then strictly increasing, guaranteeing a unique local minimum. By evaluating R at these critical points the global r_res can be computed. The authors also present a certified numerical algorithm that evaluates the minimum with a uniform error bound, avoiding time‑domain simulation.

From a control‑design perspective, the coupling parameter n acts as a tunable knob that selects a resonance class, thereby shaping the attainable depth of energy exchange and the associated beat‑time (responsiveness). Low‑order resonances enable strong, fast energy transfer (useful for absorption or damping), but they may also allow complete depletion (r_res = 0) if the arithmetic condition holds. High‑order resonances recover the conservative bound (r_res close to the initial amplitude) and provide robustness against excessive energy loss, at the cost of slower exchange. The paper formulates two complementary design problems: (i) maximize energy absorption from a host subsystem (minimise r_res) and (ii) guarantee a non‑zero retained energy (ensure r_res ≥ r_min). Both problems are solved by selecting n to achieve the desired resonance pair while respecting a prescribed maximum beat‑time.

The authors illustrate the theory with numerical examples, showing how different choices of n produce dense versus sparse phase‑plane trajectories, how the inscribed radius varies with the resonance order, and how the proposed design framework can be applied to vibration mitigation, dual‑mass gyroscope mode‑matching, and energy funneling to a single actuator.

In summary, the paper provides a rigorous geometric and number‑theoretic analysis of gyroscopic interconnections, introduces a practical performance metric (the inscribed radius), supplies exact resonance conditions and certified computation methods, and translates these insights into actionable design guidelines for both energy‑absorption and containment strategies in passive, port‑Hamiltonian systems.