Numerical simulation of the transition to chaos in a dissipative Duffing oscillator with two-frequency excitation

A mathematical modeling technique is proposed for oscillation chaotization in an essentially nonlinear dissipative Duffing oscillator with two-frequency excitation on an invariant torus in $\mathbb R^2$. The technique is based on the joint application of the parameter continuation method, Floquet stability criteria, bifurcation theory, and the Everhart high-accuracy numerical integration method. This approach is used for the numerical construction of subharmonic solutions in the case when the oscillator passes to chaos through a sequence of period-multiplying bifurcations. The value of a universal constant obtained earlier by the author while investigating oscillation chaotization in dissipative oscillators with single-frequency periodic excitation is confirmed.