In the present paper we consider and study numerically two systems based on model FitzHugh–Nagumo neurons, where in the presence of periodic modulation of parameters it is possible to implement chaotic dynamics on the attractor in the form of a Smale–Williams solenoid in the stroboscopic Poincaré map. In particular, hyperbolic chaos characterized by structural stability occurs in a single neuron supplemented by a time-delay feedback loop with a quadratic nonlinear element.

The main goal of the article is to suggest a two-dimensional map that could play the role of a generalized model similar to the standard Chirikov–Taylor mapping, but appropriate for energy-conserving nonholonomic dynamics. In this connection, we consider a Chaplygin sleigh on a plane, supposing that the nonholonomic constraint switches periodically in such a way that it is located alternately at each of three legs supporting the sleigh. We assume that at the initiation of the constraint the respective element (“knife edge”) is directed along the local velocity vector and becomes instantly fixed relative to the sleigh till the next switch. Differential equations of the mathematical model are formulated and an analytical derivation of mapping for the state evolution on the switching period is provided. The dynamics take place with conservation of the mechanical energy, which plays the role of one of the parameters responsible for the type of the dynamic behavior. At the same time, the Liouville theorem does not hold, and the phase volume can undergo compression or expansion in certain state space domains. Numerical simulations reveal phenomena characteristic of nonholonomic systems with complex dynamics (like the rattleback or the Chaplygin top). In particular, on the energy surface attractors associated with regular sustained motions can occur, settling in domains of prevalent phase volume compression together with repellers in domains of the phase volume expansion. In addition, chaotic and quasi-periodic regimes take place similar to those observed in conservative nonlinear dynamics.

A generalized model with bifurcations associated with blue sky catastrophes is introduced. Depending on an integer index $m$, different kinds of attractors arise, including those associated with quasi-periodic oscillations and with hyperbolic chaos. Verification of the hyperbolicity is provided based on statistical analysis of intersection angles of stable and unstable manifolds.

We investigate strange nonchaotic self-oscillations in a dissipative system consisting of three mechanical rotators driven by a constant torque applied to one of them. The external driving is nonoscillatory; the incommensurable frequency ratio in vibrational-rotational dynamics arises due to an irrational ratio of diameters of the rotating elements involved. It is shown that, when losing stable equilibrium, the system can demonstrate two- or three-frequency quasi-periodic, chaotic and strange nonchaotic self-oscillations. The conclusions of the work are confirmed by numerical calculations of Lyapunov exponents, fractal dimensions, spectral analysis, and by special methods of detection of a strange nonchaotic attractor (SNA): phase sensitivity and analysis using rational approximation for the frequency ratio. In particular, SNA possesses a zero value of the largest Lyapunov exponent (and negative values of the other exponents), a capacitive dimension close to 2 and a singular continuous power spectrum. In general, the results of this work shed a new light on the occurrence of strange nonchaotic dynamics.

A nonautonomous system with a uniformly hyperbolic attractor of Smale – Williams type in a Poincaré cross-section is proposed with generation implemented on the basis of the effect of oscillation death. The results of a numerical study of the system are presented: iteration diagrams for phases and portraits of the attractor in the stroboscopic Poincaré cross-section, power density spectra, Lyapunov exponents and their dependence on parameters, and the atlas of regimes. The hyperbolicity of the attractor is verified using the criterion of angles.

We investigate strange nonchaotic self-oscillations in a dissipative system consisting of three mechanical rotators driven by a constant torque applied to one of them. The external driving is nonoscillatory; the incommensurable frequency ratio in vibrational-rotational dynamics arises due to an irrational ratio of diameters of the rotating elements involved. It is shown that, when losing stable equilibrium, the system can demonstrate two- or three-frequency quasi-periodic, chaotic and strange nonchaotic self-oscillations. The conclusions of the work are confirmed by numerical calculations of Lyapunov exponents, fractal dimensions, spectral analysis, and by special methods of detection of a strange nonchaotic attractor (SNA): phase sensitivity and analysis using rational approximation for the frequency ratio. In particular, SNA possesses a zero value of the largest Lyapunov exponent (and negative values of the other exponents), a capacitive dimension close to “2” and a singular continuous power spectrum. In general, the results of this work shed a new light on the occurrence of strange nonchaotic dynamics.

We consider the movement of Chaplygin sleigh on a plane that is a solid body with imposed nonholonomic constraint, which excludes the possibility of motions transversal to the constraint element (“knife-edge”), and complement the model with an attached mass, periodically oscillating relatively to the main platform of the sleigh. Numerical simulations indicate the occurrence of either unrestricted acceleration of the sleigh, or motions with bounded velocities and momenta, depending on parameters. We note the presence of phenomena characteristic to nonholonomic systems with complex dynamics; in particular, attractors occur responsible for chaotic motions. In addition, quasiperiodic regimes take place similar to those observed in conservative nonlinear dynamics.

From analysis of time series obtained on the numerical solution of a plane problem on the motion of a body with an elliptic cross section under the action of gravity force in an incompressible viscous fluid, a system of ordinary differential equations approximately describing the dynamics of the body is reconstructed. To this end, coefficients responsible for the added mass, the force caused by the circulation of the velocity field, and the resisting force are found by the least square adjustment. The agreement between the finitedimensional description and the simulation on the basis of the Navier–Stokes equations is illustrated by images of attractors in regular and chaotic modes. The coefficients found make it possible to estimate the actual contribution of different effects to the dynamics of the body.

For a Chaplygin sleigh on a plane, which is a paradigmatic system of nonholonomic mechanics, we consider dynamics driven by periodic pulses of supplied torque depending on the instant spatial orientation of the sleigh. Additionally, we assume that a weak viscous force and moment affect the sleigh in time intervals between the pulses to provide sustained modes of the motion associated with attractors in the reduced three-dimensional phase space (velocity, angular velocity, rotation angle). The developed discrete version of the problem of the Chaplygin sleigh is an analog of the classical Chirikov map appropriate for the nonholonomic situation. We demonstrate numerically, discuss and classify dynamical regimes depending on the parameters, including regular motions and diffusive-like random walks associated, respectively, with regular and chaotic attractors in the reduced momentum dynamical equations.

Computer verification of hyperbolicity is provided based on statistical analysis of the angles of intersection of stable and unstable manifolds for mechanical systems with hyperbolic attractors of Smale – Williams type: (i) a particle sliding on a plane under periodic kicks, (ii) interacting particles moving on two alternately rotating disks, and (iii) a string with parametric excitation of standing-wave patterns by a modulated pump. The examples are of interest as contributing to filling the hyperbolic theory of dynamical systems with physical content.

Examples of mechanical systems are discussed, where quasi-periodic motions may occur, caused by an irrational ratio of the radii of rotating elements that constitute the system. For the pendulum system with frictional transmission of rotation between the elements, in the conservative and dissipative cases we note the coexistence of an infinite number of stable fixed points, and in the case of the self-oscillating system the presence of many attractors in the form of limit cycles and of quasi-periodic rotational modes is observed. In the case of quasi-periodic dynamics the frequencies of spectral components depend on the parameters, but the ratio of basic incommensurate frequencies remains constant and is determined by the irrational number characterizing the relative size of the elements.

Dynamical equations are formulated and a numerical study is provided for selfoscillatory model systems based on the triple linkage hinge mechanism of Thurston–Weeks–Hunt–MacKay. We consider systems with a holonomic mechanical constraint of three rotators as well as systems, where three rotators interact by potential forces. We present and discuss some quantitative characteristics of the chaotic regimes (Lyapunov exponents, power spectrum). Chaotic dynamics of the models we consider are associated with hyperbolic attractors, at least, at relatively small supercriticality of the self-oscillating modes; that follows from numerical analysis of the distribution for angles of intersection of stable and unstable manifolds of phase trajectories on the attractors. In systems based on rotators with interacting potential the hyperbolicity is violated starting from a certain level of excitation.

Dynamical equations are formulated and a numerical study is provided for selfoscillatory model systems based on the triple linkage hinge mechanism of Thurston – Weeks – Hunt – MacKay. We consider systems with a holonomic mechanical constraint of three rotators as well as systems, where three rotators interact by potential forces. We present and discuss some quantitative characteristics of the chaotic regimes (Lyapunov exponents, power spectrum). Chaotic dynamics of the models we consider are associated with hyperbolic attractors, at least, at relatively small supercriticality of the self-oscillating modes; that follows from numerical analysis of the distribution for angles of intersection of stable and unstable manifolds of phase trajectories on the attractors. In systems based on rotators with interacting potential the hyperbolicity is violated starting from a certain level of excitation.

Results are reviewed concerning the planar problem of a plate falling in a resisting medium studied with models based on ordinary differential equations for a small number of dynamical variables. A unified model is introduced to conduct a comparative analysis of the dynamical behaviors of models of Kozlov, Tanabe–Kaneko, Belmonte–Eisenberg–Moses and Andersen–Pesavento–Wang using common dimensionless variables and parameters. It is shown that the overall structure of the parameter spaces for the different models manifests certain similarities caused by the same inherent symmetry and by the universal nature of the phenomena involved in nonlinear dynamics (fixed points, limit cycles, attractors, and bifurcations).

Results are reviewed relating to the planar problem for the falling card in a resisting medium based on models represented by ordinary differential equations for a small number of variables. We introduce a unified model, which gives an opportunity to conduct a comparative analysis of dynamic behaviors of models of Kozlov, Tanabe – Kaneko, Belmonte – Eisenberg – Moses and Andersen – Pesavento – Wang using common dimensionless variables and parameters. It is shown that the overall structure of the parameter spaces for the different models shows certain similarities caused obviously by the same inherent symmetry and by universal nature of the involved phenomena of nonlinear dynamics (fixed points, limit cycles, attractors, bifurcations). In concern of motion of a body of elliptical profile in a viscous medium with imposed circulation of the velocity vector and with the applied constant torque, a presence of the Lorenz-type strange attractor is discovered in the three-dimensional space of generalized velocities.

For a solid body of convex form moving on a rough horizontal plane that is known as a rattleback, numerical simulations are used to discuss and illustrate dynamical phenomena that are characteristic of the motion due to a nonholonomic nature of the mechanical system; the relevant feature is the nonconservation of the phase volume in the course of the dynamics. In such a system, a local compression of the phase volume can produce behavior features similar to those exhibited by dissipative systems, such as stable equilibrium points corresponding to stationary rotations; limit cycles (rotations with oscillations); and strange attractors. A chart of dynamical regimes is plotted in a plane whose axes are the total mechanical energy and the relative angle between the geometric and dynamic principal axes of the body. The transition to chaos through a sequence of Feigenbaum period doubling bifurcations is demonstrated. A number of strange attractors are considered, for which phase portraits, Lyapunov exponents, and Fourier spectra are presented.

We consider an autonomous system of partial differential equations for a onedimensional distributed medium with periodic boundary conditions. Dynamics in time consists of alternating birth and death of patterns with spatial phases transformed from one stage of activity to another by the doubly expanding circle map. So, the attractor in the Poincaré section is uniformly hyperbolic, a kind of Smale–Williams solenoid. Finite-dimensional models are derived as ordinary differential equations for amplitudes of spatial Fourier modes (the 5D and 7D models). Correspondence of the reduced models to the original system is demonstrated numerically. Computational verification of the hyperbolicity criterion is performed for the reduced models: the distribution of angles of intersection for stable and unstable manifolds on the attractor is separated from zero, i.e., the touches are excluded. The example considered gives a partial justification for the old hopes that the chaotic behavior of autonomous distributed systems may be associated with uniformly hyperbolic attractors.

We outline a possibility of implementation of Smale–Williams type attractors with different stretching factors for the angular coordinate, namely, $n=3,\,5,\,7,\,9,\,11$, for the maps describing the evolution of parametrically excited standing wave patterns on a nonlinear string over a period of modulation of pump accompanying by alternate excitation of modes with the wavelength ratios of $1:n$.

We outline a possibility of chaotic dynamics associated with a hyperbolic attractor of the Smale–Williams type in mechanical vibrations of a nonhomogeneous string with nonlinear dissipation arising due to parametric excitation of modes at the frequencies $\omega$ and $3\omega$, when the pump is supplied by means of the string tension variations alternately at frequencies of $2\omega$ and $6\omega$.

We study numerically the dynamics of the rattleback, a rigid body with a convex surface on a rough horizontal plane, in dependence on the parameters, applying methods used earlier for treatment of dissipative dynamical systems, and adapted here for the nonholonomic model. Charts of dynamical regimes on the parameter plane of the total mechanical energy and the angle between the geometric and dynamic principal axes of the rigid body are presented. Characteristic structures in the parameter space, previously observed only for dissipative systems, are revealed. A method for calculating the full spectrum of Lyapunov exponents is developed and implemented. Analysis of the Lyapunov exponents of the nonholonomic model reveals two classes of chaotic regimes. For the model reduced to a 3D map, the first one corresponds to a strange attractor with one positive and two negative Lyapunov exponents, and the second to the chaotic dynamics of quasi-conservative type, when positive and negative Lyapunov exponents are close in magnitude, and the remaining exponent is close to zero. The transition to chaos through a sequence of period-doubling bifurcations relating to the Feigenbaum universality class is illustrated. Several examples of strange attractors are considered in detail. In particular, phase portraits as well as the Lyapunov exponents, the Fourier spectra, and fractal dimensions are presented.

We suggest a simple two-dimensional map, parameters of which are the trace and Jacobian of the perturbation matrix of the fixed point. On the parameters plane it demonstrates the main universal bifurcation scenarios: the threshold to chaos via period-doublings, the situation of quasiperiodic oscillations and Arnold tongues. We demonstrate the possibility of implementation of such map in radiophysical device.

A model system is proposed, which manifests a blue sky catastrophe giving rise to a hyperbolic attractor of Smale–Williams type in accordance with theory of Shilnikov and Turaev. Some essential features of the transition are demonstrated in computations, including Bernoulli-type discrete-step evolution of the angular variable, inverse square root dependence of the first return time on the bifurcation parameter, certain type of dependence of Lyapunov exponents on control parameter for the differential equations and for the Poincaré map.

A non-autonomous flow system is introduced, which may serve as a base for elaboration of real systems and devices demonstrating the structurally stable chaotic dynamics. The starting point is the map of the sphere composed of four stages of sequential continuous geometrically evident transformations. The computations indicate that in a certain parameter range the map posseses an attractor of Plykin type. Accounting the structural stability intrinsic to this attractor, a modification of the model is undertaken, which includes a variable change with passage to representation of instantaneous states on the plane. As a result, a set of two non-autonomous differential equations of the first order with smooth coefficients is obtained explicitly, which has the Plykin type attractor in the plane in the Poincaré cross-section. Results of computations are presented for the sphere map and for the flow system including portraits of attractors, Lyapunov exponents, dimension estimates. Substantiation of the hyperbolic nature of the attractors for the sphere map and for the flow system is based on a computer procedure of verification of the so-called cone criterion; in this context, some hints are applied, which may be useful in similar computations for some other systems.

We consider the dynamics of small perturbations of stable two-frequency quasiperiodic orbits on an attracting torus in the quasiperiodically forced Hénon map. Such dynamics consists in an exponential decay of the radial component and in a complex behaviour of the angle component. This behaviour may be two- or three-frequency quasiperiodicity, or it may be irregular. In the latter case a graphic image of the dynamics of the perturbation angle is a fractal object, namely a strange nonchaotic attractor, which appears in auxiliary map for the angle component. Therefore, we claim that stable trajectories may approach the attracting torus either in a regular or in an irregular way. We show that the transition from quasiperiodic dynamics to chaos in the model system is preceded by the appearance of an irregular behaviour in the approach of the perturbed quasiperiodic trajectories to the smooth attracting torus. We also demonstrate a link between the evolution operator of the perturbation angle and a quasiperiodically forced circle mapping of a special form and with a Harper equation with quasiperiodic potential.

Amplitude equations are obtained for a system of two coupled van der Pol oscillators that has been recently suggested as a simple system with hyperbolic chaotic attractor allowing physical realization. We demonstrate that an approximate model based on the amplitude equations preserves basic features of a hyperbolic dynamics of the initial system. For two coupled amplitude equations models having the hyperbolic attractors a transition to synchronous chaos is studied. Phenomena typically accompanying this transition, as riddling and bubbling, are shown to manifest themselves in a specific way and can be observed only in a small vicinity of a critical point. Also, a structure of many-dimensional attractor of the system is described in a region below the synchronization point.

A method is suggested for computation of the generalized dimensions for a fractal attractor associated with the quasiperiodic transition to chaos at the golden-mean rotation number. The approach is based on an eigenvalue problem formulated in terms of functional equations with coeficients expressed via the universal fixed-point function of Feigenbaum–Kadanoff–Shenker. The accuracy of the results is determined only by precision of representation of the universal function.

A method is suggested for the computation of the generalized dimensions of fractal attractors at the period-doubling transition to chaos. The approach is based on an eigenvalue problem formulated in terms of functional equations with coefficients expressed in terms of Feigenbaum's universal fixed-point function. The accuracy of the results depends only on the accuracy of the representation of the universal function.

We analyse dynamics generated by quadratic complex map at the accumulation point of the period-tripling cascade (see Golberg, Sinai, and Khanin, Usp. Mat. Nauk. V. 38, № 1, 1983, 159; Cvitanovic; and Myrheim, Phys. Lett. A94, № 8, 1983, 329). It is shown that in general this kind of the universal behavior does not survive the translation two-dimensional real maps violating the Cauchy–Riemann equations. In the extended parameter space of the two-dimensional maps the scaling properties are determined by two complex universal constants. One of them corresponds to perturbations retaining the map in the complex-analytic class and equals $\delta_1 \cong 4.6002-8.9812i$ in accordance with the mentioned works. The second constant $\delta_2 \cong 2.5872+1.8067i$ is responsible for violation of the analyticity. Graphical illustrations of scaling properties associated with both these constants are presented. We conclude that in the extended parameter space of the two-dimensional maps the period tripling universal behavior appears as a phenomenon of codimension $4$.