The hyperchaotic MathieuVan der Pol system is an important autonomous system model, with four state variables and four Lyapunov exponents, three of which are positive. Also called the TriChaos, the MathieuVan der Pol system was mathematically modeled via linear coupling of Mathieu and Van der Pol nonlinear and nonautonomous systems. In this work, a nonideal system based on the MathieuVan der Pol System is modeled considering its parameters as uncertain, which makes it closer to reality. Numerical simulations are presented demonstrating that the system has a chaotic behavior with three positive Lyapunov exponents. Given such unstable and unpredictable behavior, the linear optimal control design is proposed to reduce the chaotic movement of this system to a fixed point. The simulation results show that the identification by Linear Optimal Control is very effective.
Keywords: Nonideal system, uncertain parameters, optimal linear control
Humanity is going through a very fast technological revolution, which makes well established companies and startups strive daily to invest in the search for introducing new, customizable, autonomous and more efficient products in the work market. For such purposes, mechanical engineering contributes with projects of advanced technological products, such as the manufacturing of electric and autonomous vehicles in the automotive sector. Therefore, with a simple software update, costs related to maintenance for changing oil, filters and other parts of cars powered by internal combustion engines, will be replaced.
The challenge of creating a new product requires developing projects that are straightforward concerning the objectives to be achieved and all stages to be completed, from their beginning, development, until their end. During the execution phase of a mechanical project, the probability of mistakes happening is high, leading to undesirable effects, that is, negative effects, or even chaotic ones. Considering this, managing the risks effectively may decrease the occurrence of a chaotic event that hinders the development of a project. Formally speaking, a chaotic event, or simply chaos, represents the irregularity of a dynamical system. Its essential characteristics are: sensitivity to initial conditions, which means that its evolution in time may be altered by small perturbations, and unpredictability, which implies that it is impossible to predict the posterior evolution of the system, even knowing its state during a longtime interval.
A priori, it is possible to state that a chaotic response can be understood based on a series of contractionexpansionbending transformations, and among the tools applied to characterize chaos, the Lyapunov exponents stand out, which are one of the invariants of the system, and may be used even for evaluating other invariants [1].
Due to the possibility of a chaotic event happening during the development of projects, the alternative is to seek to mitigate this disorder state, that is, to manage the chaos that is inherent to the system, checking the weaknesses that unexpectedly arise, so as not to harm the production stage. Nowadays, Chaos Engineering contributes to carrying out experiments in the infrastructure and pointing out the weaknesses related to the system before they arise. This empirical verification process leads to invulnerable systems, and provides reliability for a complex system to operate under turbulent situations during the operational process of production. In 1990, Grebogi, Edward Ott and James Yorke demonstrated the possibility of controlling chaos [2]. Ever since then, controlling chaotic systems has been of great interest in mechanical engineering, and implies in designing state feedback control laws which stabilize the chaotic system around the instable equilibrium points. Many control techniques may be used to control a system, and each of them presents peculiarities and different results, especially in practice. To stabilize chaotic systems, in many works feedback linear control was used [10,11,18]. In the present paper, we apply the Optimal Linear Control (OLC) methodology, proposed by Rafikov and Balthazar [18], which ensures the application of linear control to nonlinear systems. Among the main characteristics of Optimal Linear Control, we highlight the use of the statespace representation model, the obtaining of optimal control signals through the Algebraic Ricatti Equation (ARE), and the objective function, parameterized by the matrices Q and R to weigh the state and control vectors, respectively [18].
This work aims at proposing a new numerical model for the MathieuVan der Pol system coupled to a source with a rotative unbalance [4,5,6,7] and considering uncertain parameters; afterwards, a control technique is applied to stabilize the whole proposed system.
Nowadays, the applications in mechanical engineering require that the systems operate with optimal parameters under certain operation conditions, which demonstrates the need for developing realistic numerical models which consider uncertainties in the parameters and inputs of the system properly. In mechanical engineering, the uncertainties have been considered in several areas, such as structure dynamics [13], rotor dynamics [14], robot manipulators [15], among others. In the process of modeling dynamical systems, uncertainties may be considered based on a probabilistic approach, through predefined intervals, or even the possibility theory [16].
The MathieuVan der Pol autonomous system is a model with four state variables and three positive Lyaponuv exponents, proposed by ShihYu et al. [3]. It is a hyperchaotic system that was modeled via linear coupling of the nonlinear and nonautonomous systems of Mathieu and Van der Pol [3], whose equations of motion are as follows:

(1) 
where , , and are the four states of the system and and are the parameters of the MathieuVan der Pol System [3,24].
The next section deals with nonideal excitations, as well as presents a proposal for a new model called NonIdeal MathieuVan der Pol System.
A system is said to be ideal when the excitation does not influence the response of the system, ie, the ideal excitation has an ideal power source [8,9,10,11,12]. On the other hand, when the excitation is influenced by the response of the system, the excitation is said to be nonideal [19]. A nonideal excitation can be modeled as shown in Figure 1 [4,11].
Figure 1. Nonideal power source [11] 
This representation on the based in Kononenko model and can be mathematically modeled according to equations (2), where is a DC motor mass, is unbalanced mass, is the dumping, is the rigidity, is the eccentricity distance from the unbalanced mass to the torque source, is the moment of inertia of the unbalanced mass, and the state variables are representation for the beam tip position and for the unbalanced mass angular position. The net torque is a function of the angular velocity and described as for [6,11,17], where a is the net torque applied by the DC motor and b is the resistive net torque constant

(2) 
The first characteristic of a nonideal system is that it presents an equation that describes the interaction of the energy source with the ideal dynamic system, and the second characteristic is that the dynamics of a nonideal system approaches the ideal case as the supplied power increases. For the present work the adimensional nonideal mathematical model configuration [6,9] of the nonideal power source is given by equations (3)

(3) 
The angular position is the nonideal excitation response, , are motor torque constants, is the distance from the unbalanced mass to the rotation center of the DC motor and is related to the moment of inertia of the system. All parameters are dimensionless positive constants. The term is due to the interaction between the dynamical system and an energy source, for example, a DC motor with limited power supply. is responsible for the Nonidealization of the system. The parameter is the applied torque constant and depends on the initial conditions, and is the resistive net torque and has no influence of initial conditions and is considered as internal damping of the DC Motor.
Making , ,,, and , and, has

(4) 
Figures 2 and 3 illustrate the behavior of the adopted dynamical model, using numerical values , and [24] for the MathieuVan der Pol System and , , and for the nonideal excitation [6].
Figure 3. Phase Portrait. (a) Mathieu system. (b) Nonideal excitation. (c) Van der Pol system 
Figure 4 illustrates the six Lyapunov exponents, ; ; ; and , , demonstrating the presence of chaos with one positive Lyapunov exponent. The total time for computing the Lyapunov exponents was with a timestep of [25].
Figure 4. Dynamics of Lyapunov expoents 
Aiming to minimize vibrations and reduce the oscillatory motion caused in the system , in the following section we propose the application of Linear Optimal Control to reduce this chaotic behavior to a fixed point.
To consider the effect of the parameter uncertainties on the performance of the controller they are added to the state and will be associated with the difference between real values and the parameters of the mathematical model with a random error of 20%, where will be the normally distributed random functions [20]. The nominal values of the parameters are given below.
The real unknown parameters of system are supposed to be as follows: , , , , , , , , , and , where are normally distributed random functions [20]. Numerical simulation results are shown in the next section, where the robustness of the methods is assessed considering that the system parameters present uncertainties.
The Optimal Linear Control design for the Nonideal MathieuVan der Pol system reduces the oscillatory movement to a small stable orbit. Due to its simplicity in terms of configuration and implementation, the linear state feedback control is especially attractive [18].
It should be stressed that this approach is analytical and it can be used without dropping any nonlinear term. Rewriting the equations of motion (4), in a state form:

(5) 
Considering a vector function , witch characterizes the desired trajectory, and the control vector , consisting of two parts: for the feedforward and for a linear feedback, in such way that,

(6) 
where is a constant matrix, and driving the deviation of the trajectory of system (6) to the desired one .
The control is optimal, in order to transfer the nonlinear system (4) from any initial to final state , minimizing the functional, where the symmetric matrix is evaluated through the solution of the matrix Ricatti differential equation [19,20,21,22]

(7) 
satisfying the final condition .
In addition, with the feedback control (4), there exists a neighborhood ,, of the origin such that if , the solution , of the controlled system (7) is locally asymptotically stable, and . Finally, if then the solution , of the controlled system (7) is globally asymptotically stable.
Using the theorem [20] the dynamic error y can be minimized ().
From equation (4):

(8) 
where the function of control is defined in equation (5).
We will obtain
, , ,,
where the controllability matrix of the system for the pair [] is obtained by .
Thus, . Then the Matrix
is given by
and (an optimal control)
).
The trajectories of the system with control can be seen in Figures 5 and 6. According to the optimal control verification [23], the function is numerically calculated by , where is positive definite.
In this work, the dynamics of a new model called the NonIdeal MathieuVan Der Pol System with chaotic behavior is analyzed through the Lyapunov exponents. Three positives exponentes was found, which demonstrates the existence of a chaotic behavior and characterizes an hyper chaotic behavior. The chaotic behavior detected in the system indicates that the system is unstable and its construction, development and efficiency are not possible, and this behavior is not desired To stabilize such chaotic behavior, a control method based on the Linear Optimal Control was implemented as a strategy for reducing the chaotic motion of the system to a desired fixed point applied to any initial condition imposed on the system.
Figures 5 and 6 illustrated the effectiveness of the control strategy applied to solve these problems. The potential results of this work may be used, in future research in projects of advanced technological products, such as the manufacturing of electric and autonomous vehicles in the automotive sector.
The authors would like to thank Universidade Estadual Paulista ”Júlio de Mesquita Filho”, especially the Faculty of Engineering of Ilha Solteira, for the space provided to carry out the work, as well as the Laboratory of Complex Systems. We would also like to thank the Universidade Federal do Estado de Mato Grosso for the support given to the work.
[1] Nayfeh, A.H., Mook, D.T. Nonlinear oscillations. John Wiley & Sons, New York, 1979.
[2] Ott, E., Grebogi, C., Yorque, J.A. Controlling chaos. Phys. Rev. Lett., 66:11961199, 1990.
[3] ShihYu, L., ShengChieh, H., ChengHsiung, Y., ZhengMing, G. Generating trichaos attractors with three positive Lyapunov exponents in new four order system via linear coupling. Nonlinear Dyn., 69:805816, 2012.
[4] Balthazar, J.M., Mook, D.T., Weber, H.I., Brasil, R.M.L.R.F., Fenili, A., Belato, D., Felix, J.L.P. An overview on nonideal vibrations. Meccanica, 38(6):613621, 2003.
[5] Dantas, M.J.H., Balthazar, J.M. On local analysis of oscillations of a nonideal and nonlinear mechanical model. Meccanica, 39:313330, 2004.
[6] Chavarette, F.R. On an optimal linear control of a chaotic nonideal duffing system. Applied Mechanics and Materials, 138139:5055, 2012.
[7] Souza, S.L.T., Caldas, I.L., Viana, R.L., Balthazar, J.M. Control and chaos for vibroimpact and nonideal oscillators. Journal of Theoretical and Applied Mechanics, 463:641664, 2008.
[8] Belato, D., Weber, H.I., Balthazar, J.M., Mook, D.T. Chaotic vibrations of a nonideal electromechanical system. International Journal of Solids and Structures, 38:16991706, 2001.
[9] Piccirillo, V., Balthazar, J. M., Pontes JR, B.R., Felix, J.L.P. On a nonlinear and chaotic nonideal vibrating system with shape memory alloy (SMA). Journal of Theoretical and Applied Mechanics, 463:597620, 2008.
[10] Chavarette, F.R., Balthazar, J.M., Felix, J.L.P. Remarks on an optimal linear control design applied to a nonideal and an ideal structure coupled to an essentially nonlinear oscillator. Journal of Computational and Nonlinear Dynamics, 5:19, 2010.
[11] Chavarette, F.R. Control design applied to a nonideal structural system with behavior chaotic. International Journal of Pure and Applied Mathematics, 863:487500, 2013.
[12] Tsuchida, M., Guilherme, K.L., Balthazar, J.M. On chaotic vibrations of a nonideal system with two degrees of freedom: 1:2 resonance and sommerfeld effect. Journal of Sound and Vibration, 282:12011207, 2005.
[13] Ghanem, R.G., Spanos, P.D. Stochastic finite elements: A spectral approach. Courier Dover Publications, 1991.
[14] Koroishi, E.H., Cavalini Jr, A.A., de Lima, A.M.G. , Stefen Jr, V. Stochastics modeling of exible rotors. Journal of the Brazilian Society of Mechanical Sciences and Engineering, 34:597603, 2012.
[15] LaraMolina, F.A., Rosario, J.M., Dumur, D., Wenger, P. Generalized predictive control of parallel robots. Proceeding of Robot Motion and Control, 422, Lecture Notes in Control and Information Sciences, 159169, 2012.
[16] Moller, B., Beer, M. Fuzzy randomness, uncertainty in civil engineering and computational mechanics. SpringerVerlag, 2004.
[17] Tusset, A.M., Balthazar, J.M., Chavarette, F.R., Felix, J.L.P. On energy transfer phenomena, in a nonlinear ideal and nonideal essential vibrating systems, coupled to a (MR) magnetorheological damper. Nonlinear Dyn., 694:122, 2012.
[18] Rafikov, M., Balthazar, J.M. On control and synchronization in chaotic and hyperchaotic system via linear control feedback. Communications in Nonlinear Science and Numerical Simulation, 13:12461255, 2008.
[19] Chavarette, F.R., Peruzzi, N.J., Balthazar, J.M., Barbanti, L., Damasceno, B.C. On an optimal linear control applied to a nonideal load transportation system, modeled whit periodic coefficients. Applied Mechanics and Materials, 5254:1318, 2011.
[20] Chavarette, F.R. Optimal linear control to parametric uncertainties in a micro electro mechanical system. International Journal of Pure and Applied Mathematics, 834:539548, 2013.
[21] Frahm, H. Adevice for damping vibrations of bodies. US Patent 989958, 1911.
[22] Viquié, R., Kerschen, G. Nonlinear vibration absorber coupled to a nonlinear primary system: A tuning methodology. Journal of Sound and Vibration, 326:780793, 2009.
[23] Cvetianin, L. Dynamics of the nonideal mechanical systems: a review. Journal of the Serbian Society for Computational Mechanics, 42:7586, 2010.
[24] Chavarette, F.R., Peruzzi, N. J., Lopes, M.L.M., Cossi, A.M. A note of optimal linear control technical applied in the chaos stabilization in the Mathieuvan der Pol autonomous oscillator. International Journal of Applied Mathematics, 25:861869, 2012.
[25] Wolf, A., Swift, J.B., Swinney, H.L., Vastano, J.A. Determining Lyapunov exponents from a time series. Physica D: Nonlinear Phenomena, 16(3):285317, 1985.
Published on 02/09/19
Accepted on 19/08/19
Submitted on 17/04/19
Volume 35, Issue 3, 2019
DOI: 10.23967/j.rimni.2019.08.001
Licence: CC BYNCSA license
Are you one of the authors of this document?