arxiv: 2604.02940 · v1 · submitted 2026-04-03 · 🌊 nlin.CD · math.DS

Recognition: no theorem link

Self-excited oscillations in multi-degree-of-freedom systems subjected to discontinuous forcing

Arunav Choudhury , R. Ganesh

Authors on Pith no claims yet

Pith reviewed 2026-05-13 18:51 UTC · model grok-4.3

classification 🌊 nlin.CD math.DS

keywords self-excited oscillationslimit cyclesstability-axis-flipping bifurcationmethod of averagingmulti-degree-of-freedom systemsdiscontinuous forcingmodal stability

0 comments

The pith

A stability-axis-flipping bifurcation controls which natural mode dominates the long-term response in multi-degree-of-freedom systems driven by discontinuous state-dependent forcing.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper shows that linear multi-degree-of-freedom oscillators subject to discontinuous forcing sustain stable limit cycles at each natural frequency. The method of averaging reduces the equations to slow-flow dynamics whose phase-plane equilibria represent these modal cycles. A stability-axis-flipping bifurcation exchanges stability between neighboring equilibria, so that initial conditions alone decide which frequency persists. Because the reduced equations yield explicit stability boundaries, the work supplies maps that predict the realized mode for given parameters. The same bifurcation structure is shown to persist when the number of degrees of freedom increases.

Core claim

In linear multi-degree-of-freedom systems driven by discontinuous state-dependent forces, averaging produces slow-flow equations whose equilibria correspond to modal limit cycles; a stability-axis-flipping bifurcation occurs when the stability axes in the phase plane invert, transferring the attracting basin from one modal cycle to another and thereby selecting the observed steady-state frequency.

What carries the argument

The stability-axis-flipping (SAF) bifurcation in the slow-flow phase plane, which inverts the stability of two modal equilibria and thereby governs the exchange of attracting basins between natural modes.

Load-bearing premise

The separation between fast oscillatory time scales and slow amplitude evolution remains large enough for the averaging approximation to stay accurate across the regimes of interest.

What would settle it

Direct numerical integration of the original discontinuous equations that records a switch in the dominant oscillation frequency exactly at the parameter value where the averaged slow-flow model predicts the stability-axis-flipping bifurcation.

read the original abstract

This study investigates the existence and stability of limit cycles resulting from self-excited oscillations in linear multi-degree-of-freedom systems subjected to discontinuous, state-dependent forcing. Using the method of averaging and slow-flow phase-plane analysis, analytical expressions are derived for the amplitudes and stability boundaries of limit cycles in a two-degree-of-freedom system. The analysis demonstrates that stable limit cycles may exist in all natural modes, with the steady-state response governed by initial conditions in regimes of multistability. A central contribution of this work is the identification and analytical characterization of the stability-axis-flipping (SAF) bifurcation, which serves as the governing mechanism for the exchange of stability between modes. The framework is then systematically extended to systems with higher degrees of freedom, confirming that the SAF bifurcation remains a universal feature, even under varying feedback configurations. The steady-state dynamics, summarized through stability maps and validated by numerical simulations, delineate the existence and stability regions of modal limit cycles as functions of key system parameters. These results provide efficient criteria for guiding optimization studies to mitigate or generate limit cycles at targeted frequencies in flexible mechanical structures.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

The paper pins down a stability-axis-flipping bifurcation via averaging that governs modal stability exchange in multi-DOF systems with state-dependent discontinuous forcing, with analytical expressions and numerics that hold up in the tested regimes.

read the letter

The main thing here is the identification of the stability-axis-flipping bifurcation in the slow-flow equations. It acts as the switch that moves stability from one mode to another as parameters change, and the authors derive explicit amplitude and boundary expressions for the two-mode case before showing the same feature persists in higher dimensions under different feedback setups. The stability maps and initial-condition dependence for multistability come out cleanly from the phase-plane analysis, and the numerical checks line up with the predictions in the regimes they plot. That gives the work a practical edge for vibration control questions in mechanical structures. The averaging step itself is standard, but the soft spot is the lack of explicit error bounds or scale-separation checks right at the SAF points. If the discontinuity jumps are not small or frequencies get close, the averaged vector field could lose the Lipschitz property the derivation needs, and the claimed universality might narrow. The numerics probably catch this in the examples shown, but tighter bounds would strengthen the boundaries they report. This is for engineers and applied dynamicists who need analytical criteria for limit-cycle amplitudes in systems with discontinuous forcing. A reader already comfortable with averaging methods will get the most out of the explicit SAF characterization and the higher-DOF extension. It deserves a serious referee because the analytical program is concrete, the numerics provide an independent check, and the central claim is falsifiable on its own terms.

Referee Report

2 major / 1 minor

Summary. The manuscript analyzes self-excited limit cycles in linear multi-DOF systems under discontinuous state-dependent forcing. Using the method of averaging and slow-flow phase-plane analysis, it derives analytical expressions for amplitudes and stability boundaries in 2DOF systems, identifies the stability-axis-flipping (SAF) bifurcation as the mechanism governing modal stability exchange, extends the framework to higher DOF while claiming universality of SAF, and validates results via numerical simulations and stability maps.

Significance. If the derivations hold, the work supplies analytical criteria for predicting and controlling modal limit cycles in systems with discontinuous nonlinearities, which is relevant for vibration mitigation in flexible mechanical structures. The explicit identification of SAF as a governing bifurcation and the extension to arbitrary DOF constitute a systematic contribution, strengthened by the combination of closed-form expressions and direct numerical validation.

major comments (2)

[Method of averaging and slow-flow analysis] Averaging derivation and slow-flow equations: the assumption that averaging remains valid near SAF points (where stability exchange occurs) lacks explicit error bounds, Lipschitz verification of the averaged vector field, or scale-separation checks for the discontinuity strength. This is load-bearing because the analytically obtained stability boundaries and the claimed universality of SAF in higher-DOF extensions rest directly on the slow amplitude equations remaining accurate at the flip locus.
[Higher-DOF extension] Extension to higher DOF: the claim that SAF remains universal under varying feedback configurations is stated without detailing how the bifurcation condition generalizes when additional modes and forcing terms are introduced; the 2DOF derivation does not automatically transfer without further assumptions on modal coupling or discontinuity locations.

minor comments (1)

[Stability maps] Stability maps and figures: the curves corresponding to the SAF bifurcation should be explicitly labeled and distinguished from other stability boundaries to improve readability of the existence regions.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive and detailed report. We address each major comment point by point below, providing the strongest honest defense of the manuscript while acknowledging where additional clarification or revision is warranted. Revisions have been made to improve rigor and clarity.

read point-by-point responses

Referee: [Method of averaging and slow-flow analysis] Averaging derivation and slow-flow equations: the assumption that averaging remains valid near SAF points (where stability exchange occurs) lacks explicit error bounds, Lipschitz verification of the averaged vector field, or scale-separation checks for the discontinuity strength. This is load-bearing because the analytically obtained stability boundaries and the claimed universality of SAF in higher-DOF extensions rest directly on the slow amplitude equations remaining accurate at the flip locus.

Authors: We agree that the manuscript would benefit from a more explicit discussion of the averaging assumptions near the SAF points. The derivation follows the standard method of averaging for systems with periodic forcing and state-dependent discontinuities, with the slow-flow equations obtained by integrating over the fast period. Numerical validation across a wide parameter range, including trajectories passing near the flip loci, shows close quantitative agreement between analytical stability boundaries and direct integration, supporting practical accuracy. However, we acknowledge the absence of explicit a priori error bounds or Lipschitz constants for the averaged vector field. We have added a dedicated paragraph in Section 2.2 discussing scale separation (requiring the discontinuity strength parameter to remain O(1) relative to the natural frequencies) and citing relevant averaging theorems for discontinuous systems. Full derivation of uniform error bounds near the bifurcation would necessitate a separate technical analysis and is noted as a limitation for future work. revision: partial
Referee: [Higher-DOF extension] Extension to higher DOF: the claim that SAF remains universal under varying feedback configurations is stated without detailing how the bifurcation condition generalizes when additional modes and forcing terms are introduced; the 2DOF derivation does not automatically transfer without further assumptions on modal coupling or discontinuity locations.

Authors: The extension rests on the modal projection of the discontinuous forcing onto the linear eigenbasis, which yields a set of slow amplitude equations whose structure is preserved for arbitrary N. The SAF bifurcation is identified as the point where the sign of the effective modal damping (or the relevant eigenvalue of the averaged Jacobian) changes for one mode while remaining stable for others. We have revised the higher-DOF section (now Section 4) to include an explicit generalization: the bifurcation condition is expressed as the vanishing of the determinant of the 2x2 modal stability matrix for each mode, with entries depending on the modal participation factors, the locations of the discontinuity surfaces, and the feedback gains. This holds under the maintained assumptions of weak nonlinearity, orthogonal linear modes, and discontinuities acting on a subset of coordinates. Additional text clarifies that strong modal coupling or coincident discontinuity locations would require a modified analysis, but the core SAF mechanism transfers directly. A new figure illustrates the stability map for a 3DOF example. revision: yes

Circularity Check

0 steps flagged

No circularity: standard averaging applied to explicit system equations yields independent stability boundaries

full rationale

The derivation begins from the stated linear multi-DOF equations with discontinuous state-dependent forcing, applies the method of averaging to obtain slow-flow equations, and extracts amplitudes plus stability boundaries directly from the resulting autonomous system. The SAF bifurcation is identified as a feature of those averaged equations (exchange of stability between modes) without any parameter being fitted to the target result or any quantity defined in terms of itself. Extension to higher DOF simply repeats the same averaging procedure on the larger system; no self-citation chain, uniqueness theorem imported from prior work, or ansatz smuggled via citation is invoked as load-bearing. The chain therefore remains self-contained against the original differential equations and does not reduce to its inputs by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 1 invented entities

The central claims rest on the applicability of the method of averaging to linear MDOF systems with discontinuous state-dependent forcing and on the existence of a slow-flow phase plane that captures amplitude and phase evolution.

axioms (2)

domain assumption The underlying system is linear between discontinuities and the forcing is strictly state-dependent and discontinuous.
Stated in the abstract as the class of systems under study.
domain assumption Time-scale separation holds so that averaging produces a valid slow-flow approximation.
Implicit in the choice of method of averaging and slow-flow phase-plane analysis.

invented entities (1)

stability-axis-flipping (SAF) bifurcation no independent evidence
purpose: Mechanism that governs exchange of stability between natural modes under variation of system parameters.
Newly identified and analytically characterized in this work; no independent evidence outside the paper is provided.

pith-pipeline@v0.9.0 · 5493 in / 1435 out tokens · 55255 ms · 2026-05-13T18:51:03.524684+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

22 extracted references · 22 canonical work pages

[1]

Modelling of Earthquake-Induced Pounding Between Adjace nt Structures With a Non-Smooth Contact Dynamics Method,

Langlade, T., Bertrand, D., Grange, S., Candia, G., and d e la Llera, J. C., 2021, “Modelling of Earthquake-Induced Pounding Between Adjace nt Structures With a Non-Smooth Contact Dynamics Method,” Eng. Struct. , 241, p. 112426

work page 2021
[2]

The Dynam ics of Flow- Induced Flutter of a Thin Flexible Sheet,

Kumar, D., Arekar, A. N., and Poddar, K., 2021, “The Dynam ics of Flow- Induced Flutter of a Thin Flexible Sheet,” Phys. Fluids, 33(3), p. 034131

work page 2021
[3]

Aeroelastic Analysis and Nonlinear Character ization of Three-Degree-of-Freedom Systems With Discontinuous No nlinearities,

Vasconcellos, R., Bouma, A., Marques, F., Abdelkeﬁ, A., and Hajj, M. R., 2022, “Aeroelastic Analysis and Nonlinear Character ization of Three-Degree-of-Freedom Systems With Discontinuous No nlinearities,” Int. J. Non-Linear Mech. , 145, p. 104100

work page 2022
[4]

Analysis of Flutter-Induced Limit Cycle Oscillations in Gas- Turbine Structures With Friction, Gap, and Other Nonlinear Contact Interfaces,

Petrov, E. P ., 2012, “Analysis of Flutter-Induced Limit Cycle Oscillations in Gas- Turbine Structures With Friction, Gap, and Other Nonlinear Contact Interfaces,” ASME J. Turbomach., 134(6), p. 061018

work page 2012
[5]

Nayfeh, A. H. and Mook, D. T., 2024, Nonlinear Oscillations , John Wiley & Sons, Hoboken, NJ

work page 2024
[6]

Quenching of Self-Excited Vib rations in Multi Degree-of-Freedom Systems: Application to Stick-Slip Mit igation in Drilling,

Ramakrishnan, G., 2019, “Quenching of Self-Excited Vib rations in Multi Degree-of-Freedom Systems: Application to Stick-Slip Mit igation in Drilling,” Proceedings of the SPE/IADC Drilling Conference and Exhibi tion, The Hague, The Netherlands, doi: 10.2118/194115-MS, SPE Paper No. SPE-194115-MS

work page doi:10.2118/194115-ms 2019
[7]

Diﬀerential Equations With Disc ontinuous Right-Hand Side,

Filippov, A. F., 1960, “Diﬀerential Equations With Disc ontinuous Right-Hand Side,” Matematicheskii Sbornik, 93(1), pp. 99–128, (in Russian)

work page 1960
[8]

A., Vitt, A

Andronov, A. A., Vitt, A. A., and Khaikin, S. E., 1966, Theory of Oscillators , Pergamon Press, Oxford, UK

work page 1966
[9]

R., and Kowalczy k, P ., 2008, Piecewise-Smooth Dynamical Systems: Theory and Applicati ons, Vol

di Bernardo, M., Budd, C., Champneys, A. R., and Kowalczy k, P ., 2008, Piecewise-Smooth Dynamical Systems: Theory and Applicati ons, Vol. 163 of Applied Mathematical Sciences, Springer, London, UK

work page 2008
[10]

Analytical Approx imations for Stick-Slip Vibration Amplitudes,

Thomsen, J. J. and Fidlin, A., 2003, “Analytical Approx imations for Stick-Slip Vibration Amplitudes,” Int. J. Non-Linear Mech. , 38(3), pp. 389–403

work page 2003
[11]

Stick-Slip Vibrations Induced by Alternate Friction Mode ls,

Leine, R. I., van Campen, D. H., de Kraker, A., and van den Steen, L., 1998, “Stick-Slip Vibrations Induced by Alternate Friction Mode ls,” Nonlinear Dyn., 16(1), pp. 41–54

work page 1998
[12]

Construction of a Stati stical Model for the Dynamics of a Base-Driven Stick-Slip Oscillator,

Lima, R. and Sampaio, R., 2017, “Construction of a Stati stical Model for the Dynamics of a Base-Driven Stick-Slip Oscillator,” Mech. Syst. Signal Process. , 91, pp. 157–166

work page 2017
[13]

Limit Cycle Oscil- lation and Dynamical Scenarios in Piecewise-Smooth Nonlin ear Systems With Two-Sided Constraints,

Cao, D.-X., Zhou, X.-X., Guo, X.- Y ., and Song, N., 2024, “Limit Cycle Oscil- lation and Dynamical Scenarios in Piecewise-Smooth Nonlin ear Systems With Two-Sided Constraints,” Nonlinear Dyn., 112(12), pp. 9887–9914

work page 2024
[14]

Multi -Harmonic Analy- sis of Dry Friction Damped Systems Using an Incremental Harm onic Balance Method,

Pierre, C., Ferri, A. A., and Dowell, E. H., 1985, “Multi -Harmonic Analy- sis of Dry Friction Damped Systems Using an Incremental Harm onic Balance Method,” ASME J. Appl. Mech. , 52(4), pp. 958–964

work page 1985
[15]

Eﬃcient Hybrid Symbolic-Numeric Computational Method for Piecewise Line ar Systems With Coulomb Friction,

Shahhosseini, A., Tien, M.-H., and D’Souza, K., 2023, “ Eﬃcient Hybrid Symbolic-Numeric Computational Method for Piecewise Line ar Systems With Coulomb Friction,” ASME J. Comput. Nonlinear Dyn. , 18(7), p. 071004. 12 / PREPRINT FOR REVIEW Transactions of the ASME

work page 2023
[16]

Analysis of Dual-Frequency Self- Excited Oscillations in Relay Feedback Systems,

Rehan, A., Boiko, I., and Zweiri, Y ., 2025, “Analysis of Dual-Frequency Self- Excited Oscillations in Relay Feedback Systems,” J. Franklin Inst. , 362(13), p. 107905

work page 2025
[17]

Self-Oscillations in Dynamic Systems,

Aguilar, L. T., Boiko, I., Fridman, L., and Iriarte, R., 2015, “Self-Oscillations in Dynamic Systems,” Self-Oscillations in Dynamic Systems: A New Methodology via Two-Relay Controllers, Springer, Cham, Switzerland

work page 2015
[18]

SPE/IADC Drilling Conference and Exhibition of SPE/IADC Drilling Conference and Exhibition

2015, Recent Advances in Soft Torque Rotary Systems , Vol. SPE/IADC Drilling Conference and Exhibition of SPE/IADC Drilling Conference and Exhibition

work page 2015
[19]

Generating Self- Excited Oscillations Via Two-Relay Controller,

Aguilar, L. T., Boiko, I., Fridman, L., and Iriarte, R., 2009, “Generating Self- Excited Oscillations Via Two-Relay Controller,” IEEE Tran s. Autom. Control, 54(2), pp. 416–420

work page 2009
[20]

On the Nonlinear Dynamics of Self-Sustained Limit-Cycle Oscillations in a F lapping-Foil Energy Harvester,

Wang, E., Ramesh, K., Killen, S., and Viola, I. M., 2018, “On the Nonlinear Dynamics of Self-Sustained Limit-Cycle Oscillations in a F lapping-Foil Energy Harvester,” J. Fluids Struct. , 83, pp. 339–357

work page 2018
[21]

Thompson, J. M. T. and Stewart, H. B., 2002, Nonlinear Dynamics and Chaos , John Wiley & Sons, Chichester, UK

work page 2002
[22]

Leine, R. I. and Nijmeijer, H., 2013, Dynamics and Bifurcations of Non-Smooth Mechanical Systems, Vol. 18, Springer, Berlin, Heidelberg. Journal of PREPRINT FOR REVIEW / 13

work page 2013