Poincaré Map-Based Continuation of Periodic Orbits in Dynamic Discontinuous and Hysteretic Systems
Abstract A numerical algorithm is proposed to compute variation of periodic solutions and their codimension-one bifurcations in discontinuous and hysteretic systems in the relevant control parameter space. For dynamic systems with discontinuities and hysteresis, some components of the associated vector fields are nondifferentiable. Therefore, one cannot resort on classical numerical tools based on the evaluation of the Jacobian of the vector field for path-following analyses. Using the pertinent state space, periodic orbits are sought as the fixed points of a Poincaré map based on an appropriate return time. The Jacobian of the map is computed numerically via either a forward or a central finite-difference scheme and a Newton-Raphson procedure is used to determine the fixed points. The continuation scheme is a pseudo-arclength algorithm based on arclength parameterization. The eigenvalues of the Jacobian of the map — Floquet multipliers — are computed to ascertain the stability of the periodic orbits and the associated bifurcations. The procedure is used to construct frequency-response curves of a bilinear, a Masing-type, and a Bouc-Wen oscillator in the primary and superharmonic frequency ranges for various excitation levels. The proposed numerical strategy proves to be very effective in capturing a rich class of solutions and bifurcations — including jump phenomena, pitchfork (symmetry-breaking), and period-doubling.