Open access peer-reviewed chapter

Graphical Method for Robust Stability Analysis for Time Delay Systems: A Case of Study

By Gerardo Romero, David Lara, Irma Pérez and Esmeralda Lopez

Submitted: November 21st 2015Reviewed: March 17th 2016Published: July 6th 2016

DOI: 10.5772/63158

Downloaded: 777

Abstract

This chapter presents a tool for analysis of robust stability, consisting of a graphical method based on the construction of the value set of the characteristic equation of an interval plant that is obtained when the transfer function of the mathematical model is connected with a feedback controller. The main contribution presented here is the inclusion of the time delay in the mathematical model. The robust stability margin of the closed-loop system is calculated using the zero exclusion principle. This methodology converts the original analytic robust stability problem into a simplified problem consisting on a graphic examination; it is only necessary to observe if the value-set graph on the complex plane does not include the zero. A case of study of an internal combustion engine is treated, considering interval uncertainty and the time delay, which has been neglected in previous publications due to the increase in complexity of the analysis when this late is considered.

Keywords

  • robust stability
  • robustness margin
  • polynomials of Kharitonov
  • value set
  • interval uncertainty

1. Introduction

The automatic control systems that are applied on the industrial process or on any production plant are previously analyzed, before they can be physically implemented. This analysis consists of obtaining some properties of interest such as qualitative and quantitative that are very helpful to achieve good performance in the closed-loop control system. The analysis is realized with the intention to predict the behavior that the control system will have when this is implemented. However, most of the time, this prediction lacks precision such that analysis is strongly based on the mathematical model of the physical process, and this does not represent in an exact way its dynamic behavior. This problem has been of great interest for scientific researchers in the last years, because it is desirable to get the best results when the control system is implemented but depends on these properties. A way of solving this problem is through the consideration of uncertainty on the mathematical model of the physical process. This uncertainty can be dynamic (see references [13]) or parametric (see references [46]). This chapter describes the qualitative property analysis of robust stability of a system. A case of study of an internal combustion engine is analyzed considering parametric uncertainty in the mathematical model. This case of study process is taken from references [5, 7], where conditions are obtained to verify the robust stability of the control system.

It is important to mention that in these research works several simplifications of the mathematical model were made, which may affect the real behavior of the system; one simplification performed by the authors is the cancellation of time delay in a part of the mathematical model, which has an influence that affects the stability property. The main contribution presented here is the consideration of the time delay on the mathematical model of an internal combustion engine, which is taken into account to obtain the property of robust stability. The methodology used is based on the application of the value-set concept to the particular case of the internal combustion engine (see reference [5]); to be more precise, it consists of the characterization of the value set of the resulting characteristic equation of the closed-loop system when the controller is connected. This controller assigns the poles in a position previously defined. Using this characteristic and applying the zero exclusion principle [8], it is possible to obtain robust stability conditions through a visual inspection of a graphic in the complex plane.

This chapter is organized as follows: Section 2 presents the elements used to verify the robust stability, Section 3 provides an abstract regarding obtaining the mathematical model of the internal combustion motor, in Section 4 the problem statement is set, in Section 5 the robust stability tools are used in simulation to verify robustness margin, and finally in Section 6 the conclusions of this work are presented.

2. Preliminary

In this section, we give the value-set concept, which is widely used to verify robust stability (see references [5, 9, 10]). The methodology that will be used consists of obtaining the value set for the characteristic equation, which results from the interval plants including time delay that are defined using the following definition:

Definition 1. A transfer function type interval plant is composed in the following manner:

g(s,q,r)=n(s,q)d(s,r)=i=0m[qi,qi+]sisn+i=0n1[ri,ri+]siqQ,rRE1

Now, if the interval plant includes the time delay, then it can be denoted by the following expression:

g(s,q,r,eτs)=g(s,q,r)eτsτ[0,τmax]E2

Notice that m < n, then g(s,q, r) is a set of strictly proper rational functions, Q and R are sets that represent the parametric uncertainty and are defined as follows:

R{r=[r1rn1]T:ririri+}Q{q=[q1qn1]T:qiqiqi+}E300

These kinds of sets are called as boxes. Equation (2) represents a type of systems known as interval plants. Both the numerator and the denominator of the transfer function have coefficients of uncertain values, which reside in a closed interval. The main interest in the study of this chapter is the robust stability analysis of feedback control systems, as the one depicted in Figure 1, whose main process is represented by an interval plant, with time delay and negative unitary gain feedback. The characteristic equation of this feedback system can be expressed in terms of the numerator, denominator, and time delay as follows:

p(s,q,r,eτs)=d(s,r)+n(s,q)eτsE3

Figure 1.

Uncertain interval plant with time delay in negative feedback.

The term quasi-polynomial is used to describe these types of functions. Note that the above characteristic Eq. (3) is not a single equation but rather a family of an infinite number of characteristic equations. Then, this whole family must be considered if robust stability verification is carried out. This family is defined as follows:

Pτ{p(s,q,r,eτs)=qQ;rR;τ[0,τmax]}E4

The robust stability property is guaranteed if and only if the following equation is satisfied:

p(s,q,r,eτs)0sC+E5

where C+is used to denote the set of complex numbers with real part positive or equal to zero. From Eq. (3), it can be comprehended that the robust stability property of the dynamic system is very hard to verify using analytical methods, because this equation contains an endless number of equations. The goal of this work is to show a simple method to validate the robust stability property of time delay systems. The contribution of this work is based on the value-set characterization of the family of characteristic equations Pτ. The value set is defined here as follows:

Definition 2. The value set, denoted as Vτ (ω), of the characteristic equation of an interval polynomial with time delay is the set of complex numbers obtained by substituting s = in the polynomial:

Vτ(ω){p(s,q,reτs):s=jω,ωR,qQ;rR;τ[0,τmax]}E6

where Qand Rrepresent the set containing all possible values of the uncertainty of the parameters qi and ri expressed in a vector form as elements of the vector q and r. It is clear that the value set of Pτ is a set of complex numbers plotted on the complex plane for values of qi, ri, ω, and τ inside the defined boundaries. An important result that is applied in this chapter is the characterization of the value set for a characteristic equation as the one considered in the previous definition; this result is represented in references [11] and [12] with the following lemma:

Lemma 1. For each frequency ω, the value set Vτ(ω) is formed by octagons, where each one changes their geometry in function of the time delay τ. The coordinates in the complex plane of the vertices or corners of each octagon are given by the following formulas:

υi+1=di+1(jω)+nk(jω)ejωτE7
υi+5=di+1(jω)+nh(jω)ejωτE8

where

i=0,1,2,3,E1000
k=(γ+i)mod4+1E1100
h=(γ+i+1)mod4+1E1200

The term mod4 represents the whole module base four operation, for example, mod4 (3). γ can take integer values 0, 1, 2, and 3 depending on ωτ, see reference [9]:

r{02nπωτπ2+2nπ1π2+2nπωτπ+2nπ2π+2nπωτ3π2+2nπ33π2+2nπωτ2π+2nπE1300

The corresponding Kharitonov polynomials for the numerator and denominator of the interval plant are denoted by ni (s) y di (s), respectively.

An example of a particular value set for fixed frequency and time delay is presented in Figure 2.

Figure 2.

Value set for a fixed frequency and time delay.

From the definition given for the value set, one can conclude that it includes all the values that the infinite family Pτ can take when it is evaluated in s = . Then, if the complex plane origin is contained in the value set Vτ(ω), this means that Pτ has roots located over the imaginary axis for some values of ωR. This, in fact, causes instability in the time delay feedback system. Consequently, the value-set technique can be used as an instrument that serves to validate the robust stability property. A question that arises in this point is how to show that Pτ does not have roots on the right half plane when a sweep over is done. The answer to this question is sustained in the following result known as the zero exclusion principle, see references [5, 11, 12]. This result will be applied to verify robust stability property.

Theorem 1: Suppose that p(s, q, r, e‒τs) has at least one member stable for τ = 0, then p(s, q, r, e‒τs) is robustly stable only if the value set satisfies:

0VT(ω)ω0E9

Proof: See reference [8]

From the previous theorem, the robust stability problem is transformed into a problem where it only needs to verify that the value-set plot just avoids the zero of the complex plane.

3. Case of study mathematical model

In this section, the internal combustion motor is modeled. The idle mode operation condition is considered, which means the vehicle engine is running without accelerating; this is because the original objective is to reduce the fuel consumption when the vehicle is on and stopped, such as waiting for the green of the traffic light when circulating in the city. This fuel economy can be achieved increasing the air ratio of the fuel mixture, but this action causes instability in engine operation, resulting in a variation in the angular speed of the motor shaft. The mathematical model was compiled from article [7] and it is divided into three parts for better comprehension: (a) manifold chamber, (b) internal combustion chamber, and (c) rotational motion system.

3.1. Manifold chamber

The rate of chance of the pressure in the manifold chamber is affected by the current chamber pressure, the opening position of the throttle valve d (t), which controls the incoming air mass flow in a proportional way, and the outgoing flow that is proportional to the motor angular velocity n (t). The output of the chamber is the relative air pressure p (t). Then, the equation that gives the relationship between these two inputs and the output is given as the first-order differential equation as follows:

p˙(t)+k2p(t)=k1d(t)k3n(t)E10

where k1, k2, y k3 are the respective proportionality constants. Applying the Laplace transform to previous equation, the following transfer function corresponding to the manifold chamber is obtained:

p(s)=1s+k2[k1d(s)k3n(s)]E11

3.2. Internal combustion chamber

The combustion chamber produces the necessary torque to move the motor shaft. The torque generation subsystem can be modeled in a simple way having as inputs the spark advance (forward position of the rotor) a(t), the relative air pressure on the chamber p(t), the motor angular velocity n(t), and the fuel flow f (t). These variables contribute to the torque generation in a linear form. There is a time delay τd called induction-to-power-stroke delay, which affects the fuel control and the chamber pressure variables; this delay τd depends on the motor speed and the number of cylinders activated independently (denoted by nc) as given by the following formula:

τd=120ncn(t)E12

Thus, the Laplace transform of the engine torque delivery Te (s) generated for the combustion block is represented by the following equation:

Te(s)=eτds[k4p(s)+k5n(s)+kff(s)]+k6a(s)E13

Most of the time the delay is neglected, but when parameter uncertainty is considered, this has an influence on the control system stability, which is the reason this has been taken into account for analysis in this work.

3.3. Shaft rotational dynamics

Finally, the equations that represent the rotational can be obtained using Newton’s second law for rotational movement as follows:

Jn˙(t)=Te(t)TL(t)k7n(t)E14

where J represents the rotational inertia, TL is the torque of external load including its disturbances, and k7 is an attenuation constant of the viscous friction that depends on the temperature, type of lubricant, and the gear that the motor uses. Therefore, the Laplace transform of the engine speed n (s) can be described as follows:

n(s)=1Js+k7[Te(t)TL(s)]E15

Then, considering the transfer functions for each subsystem, the complete block diagram representing the mathematical model of the internal combustion motor can be obtained and implemented in Simulink MATLAB, for numerical simulation as shown in Figure 3.

Figure 3.

System block diagram in Simulink MATLAB.

The interest in this case of study is the relation between the variables d (s) and n (s) taking into account the time delay; then using block diagram algebra, the corresponding transfer function is given as follows:

g(s)=k1k4eτdsJs2+(k2J+k7)s+k2k7+[k3k4k2k5k5s]eτdsE16

The previous transfer function will be considered as the mathematical model of the internal combustion motor.

4. Problem statement

It is important to mention that the engine mathematical model is expressed in terms of the parameters ki and J. These parameters depend on the motor operating point; therefore, they cannot be considered as constant in the transfer function and they will be considered as uncertain terms that change depending on the operation point, but around the nominal parameters. The nominal values considered here are the same values as [7]: k1 = 3.4329, k2 = 0.1627, k3 = 0.1139, k4 =0.2539, k5 = 1.7993, k7 =1.8201, and the inertia J = 1. Making a change on the variables, the transfer function on the mathematical model can be simplified as follows:

g(s)=a1eτdss2+a2s+a3+(a4s+a5)eτdsE17

where the uncertain parameters are the ones included in ai, the inertia J will be considered as a fixed number and equal to 1. Comparing the transfer functions (16) and (17), it can be observed that the new nominal parameters take the following values: a1 = 0.8716, a2 = 1.9828, a3 = 0.2961, a4 = 1.7993, a5 = ‒0.2638.

A controller can be connected in the closed-loop system, as shown in Figure 4, with two objectives: the first one to regulate the angular speed and the second to improve the performance of the process. This controller is given by the following function:

c(s)=50.0194s+26.3065s2+9.8165s+33.1664E18

which assigns the system poles of the feedback control, considering a unitary feedback [13], at the following position: {−1, −2, −3, −4} on the complex plane. The controller was designed with the nominal parameters and time delay equal to zero. The characteristic equation of the closed-loop control system considering the uncertain parameters and the time delay has the following structure:

p(s,q,r,eTds)=d(s,r)+n(s,q)eτdsE2400

where

d(s,r)=s4+(a2+9.82)s3+a3+9.82a2+33.17)s2+(9.82a3+33.17a2)s+33.17a3E19
n(s,q)=a4s3+(a59.82a4)s2+(9.82a533.17a4+50.02a1)s+33.17a5+26.31a1E20

The problem considered in this work is to determine the property of robust stability of the internal combustion motor when uncertainty is included in the new nominal parameters ai. This property is directly related to the characteristic equation.

Figure 4.

Control for angular velocity of the internal combustion motor.

5. Robust stability verification

First, the characteristic equation is considered without uncertainty, taking into account the nominal parameters and the time delay, which means:

p(s,q,eTds)=s4+r1s3+r2s2+r3s+r4+(q1s3+q2s2+q3s+q4)eτdsE2700

where r1 = 11.80, r2 = 52.93, r3 = 68.67, r4 = 9.82, q1 = −1.80, q2 = −17.93, q3 = −18.67, q4 = 14.18.

Using Lemma 1, the corresponding value set for ω[0,.5]yτd[0,9.9]is obtained, which is represented on the graphic shown in Figure 5.

Figure 5.

Value set without uncertainty in the parameters.

From the previous graphic, it can be appreciated that the value set does not reach the zero of the complex plain, but is very close; thus, it can be determined that the maximum delay the control system supports to preserve the stability is 9.9 s. Several numerical simulations can be run for different time delay limits to appreciate more clearly if the value set reaches the origin, for instance, considering the time delays τ1=2.5,τ2=9.9,τ3=13, which correspond to the stable, oscillatory, and unstable system responses, respectively, as shown in Figures 68, respectively.

Figure 6.

Transient response for stable behavior.

Figure 7.

Transient response for oscillatory behavior.

Figure 8.

Transient response for unstable behavior.

Finally, considering an uncertainty of 10% in each of the parameters, the following characteristic equation is obtained, where an overestimation was made to obtain interval polynomials with time delay:

p(s,q,r,eTds)=s4+r1s3+r2s2+r3s+r4+(q1s3+q2s2+q3s+q4)eτdsE2800

where r1 = [11.60,11.99],, r2 = [50.95, 54.90], r3 = [61.80,75.53], r4 = [8.84,10.80], q1 = −1.98, −1.62], q2 = [−19.72, −16.13], q3 = −29.25, −8.08], q4 = [11.01,17.35].

Once again, using the lemma, the value set of the characteristic equation is obtained now with uncertainty in the coefficients; This can be appreciated in the graphic of the figure 9.

Figure 9.

Value set with uncertainty in the parameters.

The previous value set was obtained for a value range of ω[0,1.5], τd[0,2.8]and like the previous case it can be observed that it barely reaches the zero of the complex plain; therefore, it can be assumed that the maximum time delay the control system supports is 2.8 s. Therefore, it can be clearly appreciated that parameter uncertainty decreases the delay margin that the control system supports and it is important to take this into account when the robust stability property is being verified. It is necessary to note that due to the overestimation about the parameters to represent the characteristic equation as a polynomial delay interval, the results obtained only warranty enough robust stability conditions; nevertheless, the next simulation shows how instability is presented for each of the values contained in the uncertainty, on the control system, when having a delay of τd = 4.1 s (see Figure 10).

Figure 10.

Step response with uncertainty in both parameters and time delay.

6. Discussion and conclusion of results

A robust stability verification methodology based on a graphic method was mentioned in this chapter. First, the mathematical model with nominal approximated parameters must be obtained; then, a control algorithm connected in feedback with the plant is analyzed when in the plant of the closed-loop system, uncertainty is introduced in both parameters and time delay. The main tools used are the value-set graph and the zero exclusion principle, then the analytical problem to verify the robust stability is transformed into a graphical problem, which consist of checking if the value-set plot for frequency interval does not include the origin. The case of study of an internal combustion motor was presented. The time delay appears when modeling the corresponding section of the engine combustion chamber of the motor. It is important to note that by considering the uncertainty in the parameters of the mathematical model, the time delay influences more in the control stability system; so, it can become significant even when it is very small and therefore the inclusion in the analysis is important. Then, with the tools presented here, a robustness margin for parameters and time delay can be obtained.

© 2016 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Gerardo Romero, David Lara, Irma Pérez and Esmeralda Lopez (July 6th 2016). Graphical Method for Robust Stability Analysis for Time Delay Systems: A Case of Study, Robust Control - Theoretical Models and Case Studies, Moises Rivas López and Wendy Flores-Fuentes, IntechOpen, DOI: 10.5772/63158. Available from:

chapter statistics

777total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Robust Observer-Based Output Feedback Control of a Nonlinear Time-Varying System

By Chieh-Chuan Feng

Related Book

First chapter

Substructuring Method in Structural Health Monitoring

By Shun Weng, Hong-Ping Zhu, Yong Xia and Fei Gao

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us