Home Technology On a program manifold’s stability of one contour automatic control systems
Article Open Access

On a program manifold’s stability of one contour automatic control systems

  • S. S. Zumatov EMAIL logo
Published/Copyright: December 29, 2017
Become an author with De Gruyter Brill

Abstract

Methodology of analysis of stability is expounded to the one contour systems automatic control feedback in the presence of non-linearities. The methodology is based on the use of the simplest mathematical models of the nonlinear controllable systems. Stability of program manifolds of one contour automatic control systems is investigated. The sufficient conditions of program manifold’s absolute stability of one contour automatic control systems are obtained. The Hurwitz’s angle of absolute stability was determined. The sufficient conditions of program manifold’s absolute stability of control systems by the course of plane in the mode of autopilot are obtained by means Lyapunov’s second method.

1 Introduction

The problem of construction of all set of the systems of differential equations represented by an integral curve, as well as problems of construction of systems of differential equations in different aspects were considered by N.P. Erugin, A.S. Galiullin, I.A. Mukhametzyanov, R.G. Mukharlyamov and others. It should be noted that in the process of solution of these problems, being one of basic problems of theory of stability, the problems of construction of the steady systems on a given program manifold grew into an independent theory. A great number of works is devoted to the construction of systems of equations on the given program manifold, possessing properties of stability, optimality and establishment of estimations of indexes’ quality of transient in the neighborhood of a program manifold and to solving of various inverse problems of dynamics (see [1, 11]). The detailed reviews of these works are shown in [2, 4, 9, 11]. The large number of the nonlinear systems of automatic control can be presented as one contour system consisting of two parts: linear part and nonlinear element [12]. The analysis of such systems shows that at the nonlinear elements of certain kind, transients in them can be similar, that allows to generalize the results of private experiment.

The given program Ω (t) is exactly realized only if the initial values of the state vector satisfy the condition ω(t0, x0) = 0. However, this condition cannot always be exactly satisfied. Therefore, in the construction of systems of program motion, the requirement of the stability of the program manifold Ω (t) with respect to the vector function ω should also be taken into account.

The study of these problems is due to the existence of a number of the inverse dynamics problems.

In this paper we study the problem of stability of one contour automatic control system in the neighborhood of a program manifold.

Let the differential equation

η˙=Qη,(1)

where η is n -dimension vector, Q (n × n) is some matrix, possesses by a (ns)-dimension smooth manifold Ω (t), with given the vector equation

ω(t,x)=0,(2)

where ω is sn dimension vector.

On the basis of the criterion that Ω (t) is integral for the system (1), we have

ω˙=ωt+ωηη˙=F(t,η,ω).(3)

The purpose is the construction of one contour systems on the given manifold.

We will notice that at the construction of the systems, except a condition (3), one of the basic requirements is the conditions of stability of the manifold Ω (t) in relation to some given function. The given program Ω (t) is exactly realized only if the initial values of the state vector satisfy the condition ω(t, x) = 0. However, this condition cannot always be exactly satisfied. Therefore, in the construction of systems of program motion, the requirement of the stability of the program manifold Ω (t) with respect to the vector function ω should also be taken into account.

Together with equation (1), we will consider the system of kind

η˙=Qηkf(σ),σ˙=lTPTωρn+1σ,(4)

where η, k and l are n-dimensions vector -columns:

ηT=η1,,ηn,σ=ηn+1;kT=α1,0,,0;

and P (s × n), Q (n × n) are constant matrices:

Q=ρ10000α2ρ20000α3ρ300000αnρn(5)

Nonlinear characteristic f (σ) is continuous on σ and satisfies the conditions

f(0)=0;k1σ2f(σ)σk2σ2σ0.(6)

where k1 and k2 are positive constants.

We will notice that a function f (σ) is essentially a control function on deviation from the program and at ω = σ = 0, and by virtue of conditions (6), the system (4) takes the form of (1). Consequently Ω (t) is integral for the system (4) also.

If

xi=ηi,σ=xn+1,ρi=Ti1>0αi=γiTi1=γiρi>0i1n+1,(7)

then we will get one contour system consisting of n + 1 aperiodic links, closed by a nonlinear feedback connection. Aperiodic links can serve engines of different types (electric, hydraulic, pneumatic etc.), electric generator of direct-current, reservoir with gas etc. [13]. Usually processes in such system are described by equations [14]:

(T1p+1)x1=γ1f(xn),(Tip+1)xi=γixi1,i2n+1,(8)

where ρ=ddt;

Ti, γi are positive constants:

Ti>0γi>0i1n+1.(9)

We will suppose that F = -Aω, (-A (s×s)) is Hurwitz matrix, then differentiating a manifold (2), by virtue of the system (4), we will get

ω˙=Aωωηkσ,(10)
σ˙=lTPTωρn+1σ.

Definition 1

A program manifold Ω (t) is called absolutely stable with respect to a vector-function ω, if it is asymptotically stable on the whole on the solutions of the system (10) at any ω (t0, η0) and function f (σ) satisfying conditions (6).

We will consider the linearized system

ω˙=Aωbhσ,lT=0,,0,αn+1,σ˙=cTωρn+1σ,(11)

which is derived from (10)f (σ) = hσ, b=ωηk,

c = Pl. We will suppose that it is asymptotically stable at∀ h:

k1εhk2+ε,(12)

where ɛ is how pleasingly small positive constant.

Asked, whether a program manifold Ω (t) of the system (10) will be asymptotically stable on the whole with respect to a vector-function ω for any function ∀f (σ) ϵc[k1,k2].

2 Necessary conditions of stability and determination of Hurwitz’s angle

Theorem 1

For absolute stability of program manifold Ω (t) of the system (10) with respect to a vector-function ω in Hurwitz’s angle [k1, k2], it is necessary that the system (1) was asymptotically stable at the condition (12)

Characteristic equation of the system (11) looks like

Δ(h,λ)=Q+λEhbcTρn+1+λ=m=0s+1as(h)λs+1m=0.(13)

We compose the Hurwitz matrix by means of coefficients of characteristic equation

Γm(h)==a1(h)1000a3(h)a2(h)a1(h)00am1(h)am1(h)0am(h),m1s+1(14)

We will introduce the following denotations:

ξm=infξm(h)m1s+1,(15)
ξm(h)=detΓm(h)m1s+1,(16)

then, by virtue of the Hurwitz theorem for asymptotic stability of program manifold Ω (t) with respect to a vector-function ω on the whole it is necessary the inequality

ξm>0m1s+1.(17)

Thus a Hurwitz’s angle is determined as crossing on a line

[k1ε,k2+ε]=m=1s+1ξs(h)>0h.(18)

Thus, for determination of a Hurwitz’s angle it is necessary calculations of polynomials ξm (h).

3 Algorithm of calculation of coefficients am (h)

We will define

detQ+λE=τ=0n1qτλn1τ(q0=1);(19)
detQ+λE(ρn+λ)=s=0namλnm(a0=1);(20)
Q(λ)=γ=0n2Qmλn2γ,(21)

where Q (λ) is the added matrix for ‖ Q + λE‖, coefficients qm and am are determined as:

q1=i=1n1ρi;q2=ρ1ρ2++ρn2ρn1;;qn1=i=1n1ρi;a1=i=1n1ρi;a2=ρ1ρ2++ρn1ρn;;an=i=1nρi=T1,T=i=1nTi.

and a multiplier Δ (h, λ) from (13) we will transform to the kind

Δ(h,λ)=|Q+λE|[(ρn+λ)+hlT(Q+λE)m1]==|Q+λE|(ρn+λ)+hlTQ(λ)b.(22)

By definition of the revers matrix, we have

Q+λEQ(λ)=Δ(h,λ)E,(23)

where E is unit matrix of dimension (n - 1) x (n - 1).

On the basis of formulas (20) - (22) from (23) after equating the coefficients at identical degrees λ, we will get a recurrent formula for determination of matrix- coefficients:

Qm=QQm2+qmEm1n2(Q1=0;Q0=E),(24)

what is more

QQn2=qn2E.(25)

On the basis of (13), (22) and (23) we will determine all the coefficients of the polynomial Δ (h, λ):

am(h)=am+hdmm1n,(26)
ds=lTQs2(Q2=Q1=0;Q0=E).(27)

Thus, a Hurwitz’s angle (18) is determined by virtue of correlations (14) - (16) and (26). From (4) and (5) we will get the simplified expression:

lTQ(λ)b=0,,0,αn×Qij×α1,0,,0T=i=1nαi.(28)

Consequently, we have

d1=d2==dn1=0;dn=i=1nαi=γ/T>0;(29)
γ=i=1nγi>0T>0.(30)

From (26) taking into account inequalities (9)(29), and (30) we obtained

as(h)=as>0s1n1;an(h)=an+hγ/T=T1(1+hγ)>0.(31)

As it is seen from inequalities (31), the necessary coefficient conditions of absolute stability are performed, if constants to time Ti and γi are positive. It remains only to require the inequality (17) by virtue of correlations (31).

4 Algebraic criterion of stability

For the system (10), we build the Lyapunov function of the form

V=ωTLω+β0σf(σ)dσ,(32)

where L = LT > 0; β is parameter, and we find

V˙=ωTCω+2ωTgf+βρn+1σf,(33)

where C is symmetrical three diagonal matrix

C=ATL+LA,g=Lb(β/2)cT(34)

Supposing h (σ) = f (σ) /σ from (33), we get

V˙>0,V˙=ωTCω+2ωTgh(σ)+βρn+1h(σ)σ2.(35)

For that - > 0, it is necessary and sufficiently performing of a generalized Silvester’s inequality:

G[h(σ)]=Cgh(σ)gTh(σ)βρn+1(σ)ε>0,(36)

where ɛ is a small enough positive number. As is generally known, if -A is a Hurwitz’s matrix, then C > 0, consequently, for performing of condition (36), it is enough to require that the inequality

βρn+1hgTC1gε(37)

subject to condition h (σ) > ɛ. We will set a matrix C = cij1m>0 so that hold next correlations:

detcij1mm1s.(38)

Then sufficient conditions of absolute stability C > 0 is equivalent to the necessary conditions (17).

We will define now, that

inf(βρn+1hgTC1g)=gs.(39)

Thus, if correlation L = LT > 0 ∧ β > 0, is hold subject to condition (38), (39), then on the basis of Zubov's method of the seeming linearizing is valid

Theorem 2

If –Q is a Hurwitz’s matrix and there exists L = LT > 0 ∧ β > 0, so that the equality (39) hold, then absolute stability of program manifold (2) of system (10) with respect to vector-function ω in angle (18) follows from asymptotical stability of the system (11) in the same angle.

A parameter β is determined from a condition (39). Putting the value of g in (39), we get a quadratic equation in relation to this parameter

τ1β22τ2β+τ30,(40)

where

τ1=h4(cTC1c)>0;τ2=12(LcTC1b+ρn+1).(41)

For existing of β > 0, it is necessary and sufficiently performing of inequality

τ2>0τ22τ1τ3>0.(42)

We suppose, that L = diagL1,…, Ls‖ > 0,

A=ρ10000α2ρ20000α3ρ300000αsρs,(43)

then we have

C=ATL+LA=c11c21000c21c22c32000c32c3300000cs1,s1cs,s1000cs,s1cs,s.

Here

cii=2ρiLi;cij=αiLj.(44)

We put

c11=2ρ1L1;|cij|12=g1α2L1α2L12ρL2=g2;
|cij|13=g1α2L10α2L12ρL2α3L20α3L22ρ3L3=g3;;|cij|1s=gs.(45)

From here, a recurrent formula is easily set for the coefficients of the matrix L:

Li=(gi+gi2αi2Li12)/2ρigi1i1s(g1=0;g0=1),(46)

where gi determine by the formula (45).

We calculate a vector g from (34) taking into account (45). We have

gT=α1H1,0,,0,αnβHn121.(47)

Theorem 3

Let a function (32), L is diagonal, and a matrix A has a structure (43) and condition (42) is valid. Then a matrix C is determined from (44), (45)and the program manifold (2) of system (10) is absolute stable with respect to vector-function ω in angle (18).

5 Program manifold’s absolute stability of control systems by the course of plane

We will consider now the equations of motion of control system the course of plane in the mode of autopilot [15]:

T2ψ¨+Uψ˙+kψ=0;μ˙=fT(σ)σ=aω+Eω˙+G2ω¨l1μ,(48)

that possessing by smooth manifold Ω (t), given by equation

ω(τ,ψ)=0.(49)

Here constant T2 characterizes the inertness of regulating object, U is the natural damping, k characterizes the action of restorative force, a, E, G2, l are constants of regulator.

Statement of the problem: to determine the condition of stability for the program manifold (49).

Differentiating a program manifold (49) on τ, by virtue of the system (48) we get

ω¨=2ωτ2+22ωτψψ˙+2ωψ2ψ˙2+ωψ˙UT2ψ˙μT2.(50)

We suppose, that [7]

2ωτ2+22ωτψψ˙+2ωψ2ψ˙2+ωψUT2ψ˙kT2ψF(τ,ψ,ω,ω˙),F(τ,ψ,ω,ω˙)=α1ω˙α2ω.(51)

Then we have

ω¨+α1ω˙+α2ω+ωψμT2=0μ˙=fT(σ)σ=aω+Eω˙+G2ω¨l1μ(52)

Let ωψ=h is const. Introducing denotations

ω=η1;ω˙=rη2;μ=iξ,t=τr,r=ii=lT2T2+lG22h;φ(σ)=α2r;a2=α2rp1=aα2G2;p2=(Eα1G2)r.(53)

we reduce the system (52) to the normalized form in dimensionless variables

η˙=Aηδξξ˙=φ(σ)σ=pTηξ,(54)

where

η=η1η2;A=01a2a1;δ=01;p=p1p2;φ(σ)C(o,k).

We transform the system (52) by means of substitute z1 = ɳ2, z2 = -a2ɳ1 - ɳ2 -ξ, to the form

z˙=Azδφ(σ)ξ˙=φ(σ)σ=CTzγζ.(55)

Here

CT=P1a1a2+p2,p1a2;γ=1+p1a2.

A matrix L determine as the following [4]:

l1=c2+c1a12a2+c2a1+c12a1;l2=c12a2;l3=c22a1+c12a1a2

Theorem 4

Let (-A) be a Hurwitz’s matrix.

Then for absolute stability of program manifold (49) of the system (48) it is sufficiently performing of the following conditions

γ>0c>0,c2a1+c1a1a2+p1a2=0c1a2p1a1a2+p2=0.

Theorem 5

For absolute stability of program manifold (49) of the system (48) on the angle (0, k] it is sufficiently performing of the following inequality

π(μ,q)=k1+γq+cT(A2+μE)1(A+qμE)δ>0μ0.

6 Conclusion

The conditions of stability of the one contour automatic control systems are obtained in the neighborhood of program manifold by means of construction of Lyapunov function. These results may be use for construction of stable one contour automatic control systems with respect to zero positions of equilibrium and will be developed for one contour systems with nonlinearity on type backlash, by the zone of insensitivity, satiation.

References

[1] Galiullin A., Some questions of stability of programmatic motion. Kazan.1960Search in Google Scholar

[2] Galiullin A., Mukhametzyanov I., Mukharlyamov R., Review of researches on the analytical construction of the systems programmatic motions, Vestnik RUDN, 1994, 1, 343-353Search in Google Scholar

[3] Maigarin B.G. Stability and qusality of processes of the nonlinear systems of automatic control. Alma-Ata: Nauka, 1980Search in Google Scholar

[4] Zhumatov S.S., Krementulo V.V., Maigarin B.G. Lyapunov’s Second Method in the Problems of Stability of Control of Motion Almaty: Gylym, 1999Search in Google Scholar

[5] Layton R.A. Principle of Analytical System Dynamics. – N.-Y. Springer, 199810.1007/978-1-4612-0597-5Search in Google Scholar

[6] Meiser Peter, Enge Olaf, Freudenberg Heiko and Kielau Gerald. Electromechanical Interactions in Multibody Systems Containing Electromechanical Drives, Multibody System Dynamics, 1997, 1, 281-30210.1023/A:1009729120459Search in Google Scholar

[7] Mukharlyamov R., Differential-Algebraic Equations of Programmed Motions of Lagrangian Dynamical Systems, Mechanics of Solids, 2011, 46, 4, 534-54310.3103/S0025654411040042Search in Google Scholar

[8] Mukharlyamov R., Simulation of Control Processes, Stability and Stabilization of Systems with Program Contraints, Journal of Computer and Systems Sciences International, 2015, 54, 1, 13-2610.1134/S1064230715010116Search in Google Scholar

[9] Llibre J, Ramirez R., Inverse Problems in Ordinary Differential Equations and Applications, Springer International Publishing Switzerland. 201610.1007/978-3-319-26339-7Search in Google Scholar

[10] Andreev A.S., Peregudova O.A., The stabilization of program motions of holonomic mechanical system, Automatic. Telemekh., 2016, 3, 66-8010.1134/S0005117916030048Search in Google Scholar

[11] Zhumatov S., Frequently conditions of convergence of control systems in the neighborhoods of program manifold, Nelineinye kolebania, 2016, 28, 3, 367-375Search in Google Scholar

[12] Dedesh V., Stability and auto oscillations one contour nonlinear system of automatic control, Uchenye zapiski TcAGI, 2010, XLI, 3. 82-9210.1615/TsAGISciJ.v41.i3.70Search in Google Scholar

[13] Becekerskiy V., Popov E., Theory of the systems of automatic control,. Moscow: Nauka, 1975Search in Google Scholar

[14] Trukhan N., On an one contour systems absolutely stabi8lity in Hurwitz’s angle, Avtomatic. and Telemekh., 1968, 11, 5-8Search in Google Scholar

[15] Letov A., Stability of the nonlinear control systems, Moscow: Fizmatgiz, 196210.1063/1.3058196Search in Google Scholar

Received: 2017-8-10
Accepted: 2017-11-14
Published Online: 2017-12-29

© 2017 S. S. Zumatov

This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 License.

Articles in the same Issue

  1. Regular Articles
  2. The Differential Pressure Signal De-noised by Domain Transform Combined with Wavelet Threshold
  3. Regular Articles
  4. Robot-operated quality control station based on the UTT method
  5. Regular Articles
  6. Regression Models and Fuzzy Logic Prediction of TBM Penetration Rate
  7. Regular Articles
  8. Numerical study of chemically reacting unsteady Casson fluid flow past a stretching surface with cross diffusion and thermal radiation
  9. Regular Articles
  10. Experimental comparison between R409A and R437A performance in a heat pump unit
  11. Regular Articles
  12. Rapid prediction of damage on a struck ship accounting for side impact scenario models
  13. Regular Articles
  14. Implementation of Non-Destructive Evaluation and Process Monitoring in DLP-based Additive Manufacturing
  15. Regular Articles
  16. Air purification in industrial plants producing automotive rubber components in terms of energy efficiency
  17. Regular Articles
  18. On cyclic yield strength in definition of limits for characterisation of fatigue and creep behaviour
  19. Regular Articles
  20. Development of an operation strategy for hydrogen production using solar PV energy based on fluid dynamic aspects
  21. Regular Articles
  22. An exponential-related function for decision-making in engineering and management
  23. Regular Articles
  24. Usability Prediction & Ranking of SDLC Models Using Fuzzy Hierarchical Usability Model
  25. Regular Articles
  26. Exact Soliton and Kink Solutions for New (3+1)-Dimensional Nonlinear Modified Equations of Wave Propagation
  27. Regular Articles
  28. Entropy generation analysis and effects of slip conditions on micropolar fluid flow due to a rotating disk
  29. Regular Articles
  30. Application of the mode-shape expansion based on model order reduction methods to a composite structure
  31. Regular Articles
  32. A Combinatory Index based Optimal Reallocation of Generators in the presence of SVC using Krill Herd Algorithm
  33. Regular Articles
  34. Quality assessment of compost prepared with municipal solid waste
  35. Regular Articles
  36. Influence of polymer fibers on rheological properties of cement mortars
  37. Regular Articles
  38. Degradation of flood embankments – Results of observation of the destruction mechanism and comparison with a numerical model
  39. Regular Articles
  40. Mechanical Design of Innovative Electromagnetic Linear Actuators for Marine Applications
  41. Regular Articles
  42. Influence of addition of calcium sulfate dihydrate on drying of autoclaved aerated concrete
  43. Regular Articles
  44. Analysis of Microstrip Line Fed Patch Antenna for Wireless Communications
  45. Regular Articles
  46. PEMFC for aeronautic applications: A review on the durability aspects
  47. Regular Articles
  48. Laser marking as environment technology
  49. Regular Articles
  50. Influence of grain size distribution on dynamic shear modulus of sands
  51. Regular Articles
  52. Field evaluation of reflective insulation in south east Asia
  53. Regular Articles
  54. Effects of different production technologies on mechanical and metallurgical properties of precious metal denture alloys
  55. Regular Articles
  56. Mathematical description of tooth flank surface of globoidal worm gear with straight axial tooth profile
  57. Regular Articles
  58. Earth-based construction material field tests characterization in the Alto Douro Wine Region
  59. Regular Articles
  60. Experimental and Mathematical Modeling for Prediction of Tool Wear on the Machining of Aluminium 6061 Alloy by High Speed Steel Tools
  61. Special Issue on Current Topics, Trends and Applications in Logistics
  62. 10.1515/eng-2017-0001
  63. Special Issue on Current Topics, Trends and Applications in Logistics
  64. The Methodology of Selecting the Transport Mode for Companies on the Slovak Transport Market
  65. Special Issue on Current Topics, Trends and Applications in Logistics
  66. Determinants of Distribution Logistics in the Construction Industry
  67. Special Issue on Current Topics, Trends and Applications in Logistics
  68. Management of Customer Service in Terms of Logistics Information Systems
  69. Special Issue on Current Topics, Trends and Applications in Logistics
  70. The Use of Simulation Models in Solving the Problems of Merging two Plants of the Company
  71. Special Issue on Current Topics, Trends and Applications in Logistics
  72. Applying the Heuristic to the Risk Assessment within the Automotive Industry Supply Chain
  73. Special Issue on Current Topics, Trends and Applications in Logistics
  74. Modeling the Supply Process Using the Application of Selected Methods of Operational Analysis
  75. Special Issue on Current Topics, Trends and Applications in Logistics
  76. Possibilities of Using Transport Terminals in South Bohemian Region
  77. Special Issue on Current Topics, Trends and Applications in Logistics
  78. Comparison of the Temperature Conditions in the Transport of Perishable Foodstuff
  79. Special Issue on Current Topics, Trends and Applications in Logistics
  80. E-commerce and its Impact on Logistics Requirements
  81. Topical Issue Modern Manufacturing Technologies
  82. Wear-dependent specific coefficients in a mechanistic model for turning of nickel-based superalloy with ceramic tools
  83. Topical Issue Modern Manufacturing Technologies
  84. Effects of cutting parameters on machinability characteristics of Ni-based superalloys: a review
  85. Topical Issue Desktop Grids for High Performance Computing
  86. Task Scheduling in Desktop Grids: Open Problems
  87. Topical Issue Desktop Grids for High Performance Computing
  88. A Volunteer Computing Project for Solving Geoacoustic Inversion Problems
  89. Topical Issue Desktop Grids for High Performance Computing
  90. Improving “tail” computations in a BOINC-based Desktop Grid
  91. Topical Issue Desktop Grids for High Performance Computing
  92. LHC@Home: a BOINC-based volunteer computing infrastructure for physics studies at CERN
  93. Topical Issue Desktop Grids for High Performance Computing
  94. Comparison of Decisions Quality of Heuristic Methods with Limited Depth-First Search Techniques in the Graph Shortest Path Problem
  95. Topical Issue Desktop Grids for High Performance Computing
  96. Using Volunteer Computing to Study Some Features of Diagonal Latin Squares
  97. Topical Issue on Mathematical Modelling in Applied Sciences, II
  98. A polynomial algorithm for packing unit squares in a hypograph of a piecewise linear function
  99. Topical Issue on Mathematical Modelling in Applied Sciences, II
  100. Numerical Validation of Chemical Compositional Model for Wettability Alteration Processes
  101. Topical Issue on Mathematical Modelling in Applied Sciences, II
  102. Innovative intelligent technology of distance learning for visually impaired people
  103. Topical Issue on Mathematical Modelling in Applied Sciences, II
  104. Implementation and verification of global optimization benchmark problems
  105. Topical Issue on Mathematical Modelling in Applied Sciences, II
  106. On a program manifold’s stability of one contour automatic control systems
  107. Topical Issue on Mathematical Modelling in Applied Sciences, II
  108. Multi-agent grid system Agent-GRID with dynamic load balancing of cluster nodes
Downloaded on 16.1.2026 from https://www.degruyterbrill.com/document/doi/10.1515/eng-2017-0051/html
Scroll to top button