Models for the Control of Technical Systems Motion Taking into Account Optimality Conditions Olga N. Masina Olga V. Druzhinina Alexey A. Petrov Bunin Yelets Federal Research Center Bunin Yelets State University “Computer Science State University Communards str. 28, and Control” Communards str. 28, 399770 Yelets, Russia of RAS 399770 Yelets, Russia olga121@inbox.ru Vavilov str. 44, building 2, xeal91@yandex.ru 119333 Moscow, Russia ovdruzh@mail.ru Abstract The problems of search of the optimal motion parameters for dynamical systems modeled by differential inclusions are considered. The models are studied for the cases of two and three degrees of freedom. The nonlinear generalized model is constructed for the case of six degrees of freedom taking into account influence of aerodynamic forces. The algorithm for the search of optimal motion parameters are developed with using of logical controllers and artificial neural networks. The results can find applications in the problems of synthesis and research of nondeterministic mathematical models of technical systems, in particular, in the problems of aircrafts motion control. 1 Introduction The design of mathematical models and the development of methods for their research are necessary for an adequate description of the functioning of the new classes of technical systems. The problem of optimal parameters search for the controlled technical systems motion is an actual problem arising in the modeling of many classes of the controllable technical systems. The questions of mathematical modeling for some classes of such systems are considered in [Masina, 2006], [Masina & Druzhinina, 2016], [Shestakov et al., 2014], [Druzhinina et al., 2017], [Druzhinina et al, 2017]. The examples of such systems are aircraft control systems, space vehicles, autocontrolled cars, compact unmanned flying vehicles including multicopters. The problem of optimal control becomes especially relevant for the model of the control of several unmanned vehicles so called swarm. In present paper we use the theory of differential inclusions and optimal control theory for description of the processes of dynamics and control of the technical systems. The practical application of the theory of differential inclusions is connected with the possibility of combining the results of this theory with the Pontryagin maximum principle, with stability theory and with methods of dynamic programming [Pontryagin, 1985], [Blagodatskikh & Filippov, 1985], [Aleksandrov et al., 2005], [Tolstonogov, 1986], Copyright ⃝ c by the paper’s authors. Copying permitted for private and academic purposes. In: Yu. G. Evtushenko, M. Yu. Khachay, O. V. Khamisov, Yu. A. Kochetov, V.U. Malkova, M.A. Posypkin (eds.): Proceedings of the OPTIMA-2017 Conference, Petrovac, Montenegro, 02-Oct-2017, published at http://ceur-ws.org 386 [Borisovich et al., 2011], [Merenkov, 2000], [Shestakov, 1990], [Shestakov & Druzhinina, 2002] [Masina, 2008], [Druzhinina & Masina, 2009]. In the case when we describe the model using differential inclusions, it is convenient to write assumptions on the controlled system using the multivalued function f (x, u), where x is the trajectory of the system, u is the control function. A wide class of objects, including control systems with variable control region, as well as control systems with state constraints [Pontryagin, 1985], [Blagodatskikh & Filippov, 1985], [Aleksandrov et al., 2005] can be described using the differential inclusions. The existence and stability of solutions of differential inclusions were considered in [Merenkov, 2000], [Tolstonogov, 1986], [Borisovich et al., 2011], [Masina, 2008] and in other works. The synthesis and the sequential analysis of different models of controlled systems described by differential inclusions are performed in the present paper. The mathematical models of optimal control are constructed and investigated taking into account the tracking of the point moving on the plane. First, we consider a simple two-dimensional model without taking into account the variability. Secondly, a three-dimensional model is constructed and studied taking into account the variability. Thirdly, a three-dimensional model taking into account the aerodynamic forces is constructed and studied. The algorithms of the optimal parameters search for three-dimensional models are proposed. The software package that allows selecting optimal parameters of the system motion in case of changing input data is developed. 2 The Optimal Parameters Search of the Two-Dimensional Model Initially, we consider the two-dimensional model of the controlled dynamical system with Cartesian coordinates x, y and starting point (0, 0). Let the object motion is carried out in two stages: in the time intervals (0, t1 ) and (t1 , t2 ) for the first and second stages respectively. At the first stage the object with mass m moves in the plane xoy under the influence of the constant vector thrust (p, q) until the maximum height h. At the second step the object moves under the influence of vector thrust (b, s), the boundary conditions have the form x(t2 ) = l, ẋ(t2 ) = 0, y(t2 ) = 0, ẏ(t2 ) = 0. Selecting of vector thrust values is made from strictly positive intervals p ∈ (p1 , p2 ), q ∈ (q1 , q2 ), b ∈ (b1 , b2 ), s ∈ (s1 , s2 ). The considered model takes into account the impact of the gravitational force with the magnitude of the acceleration of gravity g. The differential inclusions describing the system have the form: { mẍ ∈ p, 0 < t ≤ t1 , (1) mÿ ∈ q − mg, { mẍ ∈ −b, t < t ≤ t2 . (2) mÿ ∈ s − mg, 1 The single-valued implementations of the model (1), (2) can be represented in the form of ordinary differential equations systems { mẍ = p, 0 < t ≤ t1 , (3) mÿ = q − mg, { mẍ = −b, t < t ≤ t2 . (4) mÿ = s − mg, 1 The optimality criterion is ∫ t1 ∫ t2 (p + q)dt + (b + s)dt → min . (5) 0 t1 From (3), (4) taking into account the optimality criterion (5) we obtained that the parameters p, q, b, s, t1 , t2 in case of single-valued realization have the form 387 2lτ (τ − 1)2 gm 2hτ 2 (τ − 1)2 gm p= , q= + mg, 4l(τ − 1) + 2hτ 2 2 4l(τ − 1)2 + 2hτ 2 2lτ (τ − 1)gm 2hτ 2 (2τ − 1)gm b= , s = + mg, (6) 4l(τ − 1) + 2hτ 2 4l(τ − 1)2 2hτ 2 √ 1 4l(τ − 1)2 + 2hτ 2 τ = 1 + √ , t1 = , t 2 = τ t1 . 3 τ 2 (τ − 1)2 g The multivaluence that can be caused by the rarefied medium resistance factor in the described model takes into account when constructing a set of the single-valued realizations of the form (3), (4) in case when these implementations describe the motion of the object within a circle. The graph images of the possible trajectories of the object motion are presented in [Masina, 2006] when the endpoint is localized in the circle of positions of the plane xoy. The stability analysis of the nominal motion determined by the criterion of optimality is carried out [Masina, 2006]. The computational package Mathematica is used to find the formulas (6). Model (1), (2) is a strongly simplified model for describing the motion of the controlled object since it does not take into account the possibility of considering a mobile end point (variability condition). Besides, model (1), (2) does not take into account air resistance and requires generalization and clarification. 3 The Synthesis of the Three-Dimensional Dynamic Model The two-dimensional model modification is considered in case when the endpoint moves in the additional plane xoz. The algorithm for constructing the motion trajectory taking into account the optimal conditions, which consists of three steps is developed. Step 1. Object moves in the plane xoz under the influence of a constant thrust vector before reaching the point (0, h, 0). During the motion time t1 end point moves to the point of the positions circle. Step 2. Out of the all possible trajectories of the object motion we select the one that leads to an intermediate destination point. For the remaining time the object will reach the plane xoz taking into account the direction to the center of the circle of positions. Next we consider the motion for the time t/2 and establish a new circle of positions and select the trajectory leading to its center. Step 3. If the radius of the positions circle is larger than the radius of the conformity, we repeat the step 2. Wherein the time of the motion, when the object passes half of the route in the center of the circle of positions is halved. As a result of several repetitions of step 2, we get a circle of positions with a radius of less than σ that allows to continue to motion in the center of the circle of positions with a guaranteed hit at the endpoint. The two-dimensional model is modified for the implementation of the algorithm of the dynamical system motion. Let us choose the radius of the projection endpoint r and the radius conformity σ. In step 2 a sequential solution of equations system under condition of endpoint t2 motion along axis Ox is carried out. The system describing the object motion is converted to the form ẍ = −R1 , x(t1 ) = a1 , ẋ(t1 ) = c1 , ÿ = −S1 , y(t1 ) = a2 , ẏ(t1 ) = c2 , (7) z̈ = Z1 , z(t1 ) = a3 , ż(t1 ) = c3 , where R1 , S1 , Z1 are the acceleration components acting on the x, y, z axis, a1 , a2 , a3 are the final values of previous iterations. The derivatives c1 , c2 , c3 can be approximately calculated by the formula xk − xk−1 ci = , tk − tk−1 where k is the number of iteration of the algorithm. Given the conditions ẋ(t2 ) = 0, ẏ(t2 ) = 0, ż(t2 ) = δ, a1 = h2 , t2 = τ t1 from (7) we obtain: c1 c2 c3 R1 = − , S1 = − , Z1 = − + δ. t1 (1 − τ ) t1 (1 − τ ) t1 (τ − 1) The trajectory of motion considering optimal values R1 , S1 , Z1 corresponding to the single-valued implementation under conditions l = 10, h = 5, g = 9.8, m = 1, σ = 1, a1 = 0, a2 = 0, a3 = 0, c1 = 0, c2 = 0, c3 = 0 has the view shown in Figure 1. 388 Figure 1: The trajectory of motion taking into account variability The set of single-valued realizations of the model describes the motion of the system taking into account the optimality and variability. Thus, the model of the dynamical system described by differential inclusions is constructed taking into account optimality and variability. The software package in the built-in language of the Octave 4.0 system is elaborated on the basis of the algorithm [Masina, 2006]. 4 Generalization of the Three-Dimensional Model and the Optimal Parameters Search Next, we will consider the modification of the three-dimensional model taking into account the influence of aerodynamic forces. Let the object of control has the initial coordinates (0, y0 , 0); the initial velocity is (x¯0 , y¯0 , 0). The motion is being considered at the finite time interval T ∈ [t0 , t∆ ]. The endpoint moves according to a certain law, so that its position can only be determined in each moments ti , (i = 1, 2, ..., n) of motion. The center of the object mass is affected by influencing the thrust vector A as well as the force of gravity. The aerodynamic force acting on the object consists of three components – the drag, the lifting force, directed orthogonally to the motion vector, and the aerodynamic moment. The value of the axial position (α, β, γ) is also considered. Let us assume that the motion is accessible if the coordinate y is decreasing continuously and equals zero at the end point L(lx , lz , 0), as well as stable if the object axis coincides with the direction of motion. Object motion equations have the form 1 ẍ = (Ax + (ϕx (α, β, γ, ẋ) + ψx (α, β, γ, ẏ, ż))ẋ2 ), m 1 ÿ = (Ay − g + (ϕx (α, β, γ, ẏ) + ψy (α, β, γ, ẋ, ż))ẏ 2 ), m 1 z̈ = (Az + (ϕx (α, β, γ, ż) + ψz (α, β, γ, ẋ, ẏ))ż 2 ), m 1 (8) α̈ = (Bα + µα (α, β, γ, ẋ, ẏ, ż)), κα 1 β̈ = (Bα + µβ (α, β, γ, ẋ, ẏ, ż)), κβ 1 γ̈ = (Bγ + µγ (α, β, γ, ẋ, ẏ, ż)), κγ where ϕ, ψ, µ are the drag parameters, the lifting force and aerodynamic moment, respectively defined for each object form experimentally, A, B are the controlled actions in the respective measurements, κ is the axial moment 389 of inertia. Optimality criterion is ∫ t∆ (|Ax | + |Ay | + |Az | + |Bα | + |Bβ | + |Bγ |)dt → min . t0 From the algorithm in the Section 3 let us represent the time interval T as a sum of subintervals species ∑ n t∆ − ti T = . i=0 2 The final trajectory E of system can be represented as a sum of intermediate trajectories in each interval [ti , ti+1 ] ∈ T as ∑ n E= Mi (t), ti ≤ t < ti+1 , i=0 where Mi is vector (xi , yi , zi ). The considered restrictions have the form xi (t∆ ) = li , yi (t∆ ) = 0, zi (t∆ ) = ki , Ṁi (t∆ ) = 0, where li and ki are the results of measurement of endpoint displacement. The search for optimal values of Ai and Bi is the search minimization functional min(Ai , Bi ) = = f (Pi , Vi , ti , t∆ , ϕ, ψ, µ) where Pi , Vi are initial coordinates and velocity vectors respectively. The algorithm of optimal parameters finding using artificial neural networks [Haykin, 2006] and logical controllers [Gostev, 2008] for model (8) is designed. The graph of the algorithm of the optimal parameters research includes components such as a fuzzyfication, a rule base, a defuzzyfication, a neural network, sets of input and output parameters of the model (8). The analysis of the nonlinear model to assess training error criterion is performed. The training process means weighting coefficients tuning. The training criterion is the following   2   2 li 0 g(ti )  0  − Mi (t∆ ) + z(ti )  0  − Ṁi (t∆ ) → min, ki 0 where Mi , Ṁi calculated along the optimal values Ai , Bi . The functions g(t) and z(t) in each interval [ti , ti+1 ] are analogues to the membership functions, which allows us to carry out braking when approaching the endpoint. The advantage of the proposed algorithm is the opportunity to study the model in case when analytical research is difficult. Besides, the algorithm is characterized by the insignificant computational complexity. In the considered case it is possible to use Takagi-Sugeno controllers [Tanaka & Wang, 2001], [Druzhinina & Masina, 2009], which allow to carry out the “downward” method of modeling. The adjustment of the parameters with such modeling occurs in the course of a more detailed study of the dynamic system. 5 Conclusions The considered models show different approaches to the search of optimal parameters of control for technical systems. The first approach is connected with construction of the dynamical model in two-dimensional space taking into account the multivalence and criteria of an optimality. The second approach is the development of the first approach and second approach is based on the construction of model in three-dimensional space taking into account the multivalence and variability. The third approach is the generalization of the first and second approaches in case when air resistance is taken into account. The algorithms for finding of the optimal parameters of control including on the basis of artificial neural networks are proposed. These algorithms can be implemented as a computer programs. The obtained results can be used for a wide range of applications, in particular, the problems of autopiloting car for the destinations taking into account topography and external conditions, the problems of autonomous control for reusable spacecraft in long-term research missions with regard to optimality of motion in gravitational complex configuration fields. 390 References [Masina, 2006] Masina, O.N. (2006). The problems of the motion control of transport systems. Transport: science, technology, control, 12, 10–12. [Masina & Druzhinina, 2016] Masina, O.N. & Druzhinina, O.V. (2016). On optimal control of dynamical systems described by differential inclusions. Proceedings of the VII International conference on optimization methods and applications ”Optimization and application” (OPTIMA-2016) held in Petrovac, Montentgro, September 25 – October 2. (pp. 104–105). Moscow: Dorodnicyn Computing Centre of FRC CSC RAS. [Shestakov et al., 2014] Shestakov, A.A., Druzhinina, O.V. & Masina, O.N. (2014). Safety assessment of railway vehicles motion on the base of generalized technical stability and stability by Zhukovsky. Transport: science, technology, control, 2, 3–8. [Druzhinina et al., 2017] Druzhinina, O.V., Masina, O.N. & Petrov, A.A. (2017). Model of motion control of transport system taking into account conditions of optimality, multivaluence and variability. Transport: science, technology, control, 4, 3–9. [Druzhinina et al, 2017] Druzhinina, O.V., Masina, O.N. & Petrov, A.A. (2017). Approach elaboration to solution of the problems of motion control of technical systems modeled by differential inclusions. Information-measuring and controlling systems, 15, 4, 64–72. [Pontryagin, 1985] Pontryagin L.S. (1985). The mathematical theory of optimal processes and differential games. Proceedings of Steklov Mathematical Institute, 169, 119–158. [Blagodatskikh & Filippov, 1985] Blagodatskikh, V.I. & Filippov, A.F. (1985). Differential inclusions and optimal control. Proceedings of Steklov Mathematical Institute, 169, 194–252. [Aleksandrov et al., 2005] Aleksandrov, V.V., Boltyanskii, V.G., Lemak, S.S., Parusnikov, N.A. & Tikhomirov, V.M. (2005). Optimal motion control. Moscow: Fizmatlit. [Tolstonogov, 1986] Tolstonogov, A.A. (1986). Differential inclusions in Banach space. Novosibirsk: Nauka. [Borisovich et al., 2011] Borisovich, Yu.G., Gelman, B.D., Myshkis, A.D. & Obukhovskii, V.V. (2011). Introduction to the theory of multivalued mappings and differential inclusions. Moscow: URSS. [Merenkov, 2000] Merenkov, Yu. N. (2000). Stability-like properties of differential inclusions, fuzzy and stochastic differential equations. Moscow: Peoples Friendship University of Russia. [Shestakov, 1990] Shestakov, A.A. (1990). Genralized direct Lyapunov method for systems with distributed parameters. Moscow: Nauka. [Shestakov & Druzhinina, 2002] Shestakov, A.A. & Druzhinina, O.V. (2002). Generalized direct Lyapunov method for the analysis of stability and attraction in general time systems. Sb. Math. 193, 10, 1411– 1441. [Masina, 2008] Masina, O.N. (2008). On the existence of solutions of differential inclusions. Differential Equations, 44, 6, 845–847. [Haykin, 2006] Haykin, S. (2006). Neural networks: complete course. Moscow: Williams. [Gostev, 2008] Gostev, V.I. (2008). Fuzzy controllers in automatic control systems. Kiev: Radioamateur. [Tanaka & Wang, 2001] Tanaka, K. & Wang, H.O. (2001). Fuzzy control systems design and analysis: a linear matrix inequality approach. N.Y.:Wiley. [Druzhinina & Masina, 2009] Druzhinina, O.V., & Masina, O.N. (2009). Methods of stability research and controllability of fuzzy and stochastic dynamical systems. Moscow: Dorodnicin Computing Centre of RAS. 391