On delayed loss of stability in one mechanical problem E.V. Shchetinina1 1 Samara National Research University, 34 Moskovskoe Shosse, 443086, Samara, Russia Abstract Consider a double pendulum under the external force. This model can be described by the multi-speed systems of ordinary differential equations. We study the existence of stable motions depending on the parameters of the system and on the value of the external force. Critical values of the parameters are defined. Ключевые слова: slow-fast systems; delayed loss of stability; double pendulum 1. Introduction Numerical investigation of dynamical systems can be very complicated due to the high dimensions of systems and presence of big and small parameters. In such cases it is reasonable to use combination of numerical and analytical methods. By analytical methods we can decrease the dimension of the model, restrict the domain of the investigation, define the different types of behavior depending on the parameters values. And then by using numerical methods we are able to investigate specific solutions and manifolds of them. One such analytical method is the method of the integral manifolds. With the help of it we are able to decrease the number of variables, to restrict the domains of the investigation, to predict different types of behaviour of the solutions to the system. In this paper we consider a double pendulum with elastic hinges under the external force. We assume that the force is slowly changing in time. Our goal is to investigate the behavior of the solutions for different values of the parameters of the system and to find the influence of the force growth to the pendulum motions. 2. Integral manifolds Integral manifolds method is an efficient tool for studying complicated dynamical systems. It was developed by many authors. This method is quiet successful for studying multi-scaled systems (see, e.g. [3]). The main ideas if this method are as follows. Consider the singularly perturbed system: dx dy =f(x,y,a,ε),ε dt =g(x,y,a,ε), (1) dt with 𝑥 ∈ 𝑅 ,y ∈ 𝑅 , t ∈ 𝑅, a is a parameter, 0<ε ≪ 1. Here x is a slow variable, y is a fast variable. Integral manifold of such a system is an invariant set of the system. We are interested in the integral manifolds of the form y=h(x, ε), with h smoothly depending on ε. This type of manifolds is called slow integral manifolds. The motion on this manifolds is provided by the equation 𝑥˙ =f(x,h(x,ε),ε). Suppose that the degenerated equation g(x,y,a,0)=0 has an isolated root y=h0(x,a). The surface defined by the relation y=h0(x,a) is called slow manifold. Consider the Jacobi matrix 𝜕 𝑔⁄𝜕 𝑦(𝑥, y, 𝑎) with y=h0(x,a). If all eigenvalues of the Jacobi matrix are in the left open complex half-plane, then the slow manifold is attractive. If there exists at least one eigenvalue in the right complex half-plane, then the slow manifold is repelling. In the ε-neighbourhoods of the attracting and repelling slow manifolds there exist attracting and repelling slow integral manifolds. The surface with the condition that the eigenvalues of the Jacobi matrix are on the imaginary axis is called the surface of change of attractivity. Due to the presence of the additional parameter a we are able to glue together attracting and repelling slow integral manifolds at one point. Thus we obtain the solution that follows first an attractive slow manifold and then repelling slow manifold. Such a solution is called a canard solution. Also there exists another type of change of attractivity. Suppose that in the spectrum of the Jacobi matrix there exists a pair of complex conjugated eigenvalues which cross an imaginary axes from the left to the right with nonvanishing speed. In this case the trajectories of the system starting in the small neighbourhood of the attracting part of the slow manifold follow it until the point of change of attractivity. But after crossing this point the trajectories do not leave this small neighbourhood immediately. The loss of attractivity is delayed: for some time the trajectories stay near repelling part of the slow manifold and then jump away. This phenomenon was described in [2]. Let us mention that the presence of the additional parameter allow us to change the time of the attractivity loss delay [4]. Ziegler pendulum We consider the Ziegler system which describes a double pendulum with elastic hinges under the external force P (fig. 1) [5]: 3rd International conference “Information Technology and Nanotechnology 2017” 184 Mathematical Modeling / E.V. Shchetinina Fig. 1. Ziegler pendulum. The behaviour of such a system can be described by the Lagrange equations. In the case that the spring rigidity is high, then the system has a big parameter. Suppose that the strength of the force P is slowly changing according to the law 𝑝 = 𝑝 + 𝜀𝜏, where 0<ε ≪ 1. It means that the system becomes multispeed. By using the integral manifold theory we obtain the following system of ordinary differential equations 𝑝˙=ε, 𝜑 = −𝜀 𝜑 +κ (1+µ)cos(𝜑 ) + cos(𝜑 ) +p sin(χφ − 𝜑 ) − sin (𝜒 − 1)𝜑 , 𝜑 = −𝜀 𝜑 +κ 2cos(𝜑 ) + (1+µ)cos(𝜑 ) +p 2sin (𝜒 − 1)𝜑 − sin(χφ − 𝜑 ) . Here κ is a small parameter inversely related to the spring rigidity, µ=m ⁄𝑚 . The system possesses the stationary state 𝜑 =φ = 0. Linearization of the fast subsystem near this manifold leads to the form 𝑝˙=ε, 𝜑 = −𝜀 𝜑 (𝑝 − 1)+φ 𝑝(1 − 2𝜒) , 𝜑 = −𝜀 pφ − 𝜑 1+p(3𝜒 − 2) . The characteristic equations is λ+1 − 3⁄2 𝑝(1 − 𝜒) +p ⁄4 (1 − 𝜒)(9𝜒 − 5) = 0. From this, it follows that for 0<χ< 5⁄9 the last equation has two real eigenvalues, and for 5⁄9 <χ<1 it has a pair of complex conjugated eigenvalues. In the case of real eigenvalues there exists the point pcr defined by 𝑝cr = 3⁄2 − 1⁄2 (5 − 9𝜒)⁄(1 − 𝜒) Such that for ppcr one eigenvalue becomes positive. Therefore the slow manifold 𝜑 = 𝜑 = 0 is attracting for ppcr . Therefore the solution 𝜑 = 𝜑 = 0 is a canard solution. All other solutions, starting for ppcr the real part of them becomes positive. It means that the solutions of the system possesses delay of loss of attractivity. The trajectories starting at the point with p5/9. References [1] Gorelov GN, Sobolev VA, Shchepakina EA. Singularly perturbed combustion models. Samara, 1999; 198 p. (in Russian) [2] Neishtadt AI. Persistance of stability loss for dynamical bifurcations, I. Differents. Uravn, 1987; 23(12): 2060–2067. II. Differents. Uravn, 1988; 24(2): 226–233. (in Russian) [3] Strygin VV, Sobolev VA. Separation of motions by the integral manifolds. M.: Nauka, 1988; 256 p. (in Russian) [4] Shchetinina EV. Integral manifolds for slow-fast systems and delayed loss of stability. Vestnik SGU 2010; 6(80): 93–105. (in Russian) [5] Ziegler H. Die Stabilitaetskriterien der Elastomechanik. Ingenieur-Archiv, 1952; Band XX: 49–56. 3rd International conference “Information Technology and Nanotechnology 2017” 186