Torus and Subharmonic Motions of a Forced Vibration System in 1:5 Weak Resonance

The Neimark-Sacker bifurcation of a forced vibration system is considered in this paper. The series solution to the motion equation is obtained, and the Poincaré map is established. The fixed point of the Poincaré map is guaranteed by the implicit function theorem. The map is transformed into its normal form at the fifth-order resonance case. For some parameter values, there exists the torus T1. Furthermore, the phenomenon of phase locking on the torus T1 is investigated and the parameter condition under which there exists subharmonic motion on the torus T1 is determined.


Introduction
In this article, we investigate the torus and subharmonic motions of the following system: _ y 1 = −y 2 , _ y 2 = y 1 + a 1 μy 2 + by 2 2 + sy 3 1 + cy 3 2 where μ, ε, and δ are small parameters; f is a 2π periodic function; and a 1 , s, b, and c are constants. ðq + δÞ is the frequency of the external force. If q is a positive integer with q ≥ 5, we say that the system (1) is in 1 : q weak resonance. For q being a positive integer satisfying q ≤ 4, the system (1) is referred to as 1 : q strong resonance. h:o:t: represents the "higher-order terms" than those that have been written out, and the same is true below. There are some mechanical models whose dynamical behaviors can be described by Equation (1), for example, the system shown in Figure 1, see Ref. [1] for more information. In this paper, we investigate abstractly system (1) and give a method for analyzing its Neimark-Sacker bifurcation. The value of relevant param-eters and constants depends on some specific mechanical models whose dynamics can be described by Equation (1). Therefore, we do not introduce the given parameters or constants here, and only choose several sets of values for numerical simulations in Section 4. For ε = 0, Equation (1) undergoes the Hopf bifurcation under certain conditions, and then, for ε ≠ 0 and jεj, sufficiently small torus or the qth order subharmonic motions can occur to Equation (1). The problem of 1 : q resonance of a closed orbit in R 3 (or in C × S 1 ) leads to the study of the Z q -equivariant planar vector (see Refs. [2,3]) whose versal unfolding has been studied for q ≠ 4 and been conjectured for q = 4 by Arnold [2]. Bifurcation sequence inventory at 1 : 4 resonance has been presented by Krauskopf [3]. Gambaudo [4] considered the general study of the periodic perturbation of a one family of autonomous differential equations in the plane satisfying conditions for a generic Hopf bifurcation. Iooss [5] investigated the subharmonic motion in the 1 : 3 and 1 : 4 resonance case. Wan [6] analyzed the Neimark-Sacker bifurcation in the 1 : 4 strong resonance case for the planar map. The width of the resonance tongue at a distance σ from the unit circle given by Arnold [2,7] is of the order of σ ðq−2Þ/2 . Shilnikov et al. [8] and Iooss [9] computed the Arnold tongue in weak resonance case for the planar map. The Neimark-Sacker bifurcation of an oscillator with dry friction was observed in Ref. [1]. Periodic-impact motions and bifurcations of vibroimpact systems near the 1 : 4 strong resonance point are considered in Ref. [10]. Results on other types of forced vibrations can be found in some literatures, see, for example, Refs. [11][12][13][14][15][16].
Judging from the above statement, we know that there is a lack of efficient criteria depending on the coefficients of the original differential equations, based on which we can talk about the asymptotic behaviors of trajectory. In this paper, we will restrict our attention to the case of 1 : 5 resonance, namely, q = 5 in Equation (1), and obtain criteria. For higher-order resonance, the procedure is essentially the same as this case but needs tedious computation. This paper is organized as follows. In Section 2, the Poincaré map is established according to power series solution to Equation (1). The map is further transformed into its normal form. In Section 3, the Neimark-Sacker bifurcation is investigated. When there is a circle bifurcating from the fixed point, the phenomenon of phase locking on the invariant circle is studied and the parameter region in which subharmonic motion can occur is determined. In Section 4, choosing a set of parameters, the theoretical results stated above are verified by numerical simulations.

The Numerical Simulations
To illustrate the results stated above, numerical simulations will be presented in this section. As mentioned in the introduction, f is a 2π periodic function. For simplicity, let f ðð5 + δÞtÞ = sin ðð5 + δÞtÞ in Equation (1), which is a simple form of f ðð5 + δÞtÞ. It follows from (6)- (14) that Furthermore, we take the set of parameters a 1 = 2, b = 1, s = 0, c = −4: We can calculate that Re ðα 0 λ 0 Þ = −1:8850 < 0 and the derivation of Re ðλ 1 Þ with respect to parameters μ and ε is nonzero, which means that the supercritical Neimark-Sacker bifurcation takes place for map (25) [17]. 4 Advances in Mathematical Physics

Advances in Mathematical Physics
Choosing ε = −0:3, μ = −0:2, and δ = 0, we get Re ðλ 1 Þ = −0:1919 and then assert that map (25) possesses a stable fixed point (see Figures 2(a) and 2(b)), namely, a stable periodic solution of Equation (1) (shown in Figures 2(c) and 2(d)). If the dynamic behaviors of the model in Figure 1 can be described by Equation (1) with the present parameters, the period of vibration of the mass body is the same as the external force.
Choosing ε = −0:3 and μ = 0:2, we have Re ðλ 1 Þ = 0:3108 and then assert that map (25) possesses a stable invariant circle [17], namely, a stable torus motion of Equation (1). Because the limitation of map (25) on the stable invariant circle is a circle diffeomorphism, the trajectory on the torus is quasiperiodic or subharmonic motion, which depends on expression (24). As will be investigated below.
For δ = 0, By straight computation, we obtain  Figures 4(c) and 4(d)). If the dynamic behaviors of the model in Figure 1 can be described by Equation (1) with the present parameters, the vibration of the mass body is subharmonic, whose period is five times than that of the external force.
Because the analysis method of this paper is for system (1) with abstract coefficients, it can be applied in other mechanical models whose dynamics can be described by Equation (1), for example, the forced Van der Pol equations [6], the forced dry friction system [18], the vibration of railway bow net, and the forced vibration of cantilevered flowconveying pipe.

Conclusions
In this paper, we study the Neimark-Sacker bifurcation of a forced vibration system by theoretical analysis and numerical simulations in the 1 : 5 resonance case. The Poincaré map is established by the analytical method. By means of analyzing the map, it is shown that there exist quasiperiodic and subharmonic solutions on the torus. Numerical simulations agreed with the theoretical results. It is certain that the method applied in this paper can be applied to some other analogous systems.

Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
The author declares that they have no conflicts of interest.