Научная статья на тему 'SYNTHESIS OF ACTIVE DISTURBANCE REJECTION CONTROL'

SYNTHESIS OF ACTIVE DISTURBANCE REJECTION CONTROL Текст научной статьи по специальности «Электротехника, электронная техника, информационные технологии»

CC BY
196
76
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
GENERAL DISTURBANCES / ESTIMATOR / TRACKING DIFFERENTIAL / NONLINEAR STATE ERROR FEEDBACK

Аннотация научной статьи по электротехнике, электронной технике, информационным технологиям, автор научной работы — Alhelou M., Gavrilov A.I.

Synthesis of Active Disturbance Rejection Control (ADRC) has been discussed in this work. Two main approaches have been presented for the second-order ADRC: Linear ADRC and nonlinear ADRC. A design procedure for the three main components of ADRC: controller, estimator, and disturbance rejection scheme has been presented in each approach. Parameters of the linear approach are tuned using the desired closed loop system bandwidth. Parameters of the nonlinear approach are selected by categorizing them first into usual parameters and key parameters. After that, the key parameters are optimized using Genetic Algorithm (GA). The two approaches have been tested on a quarter-car system that deals with the passenger comfort problem. Simulation results show a good performance and a good compensation for external disturbances and dynamic uncertainties

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Текст научной работы на тему «SYNTHESIS OF ACTIVE DISTURBANCE REJECTION CONTROL»

UDC 303.732.4:629.1

DOI: 10.18698/0236-3933-2020-4-22-41

SYNTHESIS OF ACTIVE DISTURBANCE REJECTION CONTROL

M. Alhelou A.I. Gavrilov

[email protected] [email protected]

Bauman Moscow State Technical University, Moscow, Russian Federation

Abstract

Synthesis of Active Disturbance Rejection Control (ADRC) has been discussed in this work. Two main approaches have been presented for the second-order ADRC: Linear ADRC and nonlinear ADRC. A design procedure for the three main components of ADRC: controller, estimator, and disturbance rejection scheme has been presented in each approach. Parameters of the linear approach are tuned using the desired closed loop system bandwidth. Parameters of the nonlinear approach are selected by categorizing them first into usual parameters and key parameters. After that, the key parameters are optimized using Genetic Algorithm (GA). The two approaches have been tested on a quarter-car system that deals with the passenger comfort problem. Simulation results show a good performance and a good compensation for external disturbances and dynamic uncertainties

Keywords

General disturbances, estimator, tracking differential, nonlinear state error feedback

Received 18.12.2019 Accepted 12.05.2020 © Author(s), 2020

Introduction. Active Disturbance Rejection Control (ADRC) is a control method or a strategy intended to transcend complexity of the traditional control methods such as proportional-integral-derivative (PID). ADRC is inherited from PID and can be considered as a robust control method. It represents all unknown dynamics which are not included in the mathematical model of the controlled system. Moreover; it compensates for the modeling uncertainties and external perturbations in real time.

ADRC has dynamic methodology and unnecessary knowledge of a plant precise model [1-3]. Therefore, various applications use ADRC as a replacement for the classic controllers.

ADRC is a process that aims to reject the unknown internal dynamics of a plant and the external disturbances in real time. It depends on forcing the controlled plant to behave as an n-order integrator system which is easily

controlled by a PD controller, even if the plant is nonlinear and time-variant. For doing so, an estimation for the nonlinear and time-variant internal dynamics has been performed. Moreover; a prediction of the external disturbances is conducted to reject these disturbances from the whole system. The unknown dynamics and the external disturbances could be combined together in "General Disturbances" then estimated and rejected.

ADRC has been widely used since its beginning in the nineties of the last century. It inherits from PID little overshoot, high convergence speed, high accuracy, strong anti-interference ability and other characteristics [4]. It depends on expanding the system model to include everything the user does not include in the mathematical description. Zhu E., et al. in [5] proposed an ADRC control based method for an airship to track a trajectory. Internal dynamic uncertainties and external disturbances are estimated by an Extended State Observer (ESO). ADRC is used to realize the decoupling control for multivariable systems. Simulation results show steady and rapid responses. Results have been compared with a sliding mode control (SMC) results and showed more steady performance.

Li X., et al. in [6] used ADRC in order to improve the control of permanent magnet brushless motor. The authors presented the second order form of ADRC. Discrete equations of the tracking differentiator (TD), the extended state observer (ESO), and the nonlinear-state error feedback (NLSEF) were shown. Simulation results of the motor controlled system, based on ADRC, were displayed. The anti-disturbance capability of the system was presented by adding a square wave at a specific time point and system reaction was observed. Simulation results showed that TD and ESO can ensure the ability of tracking the input and output.

Gao Z., et al. in the United States patent [7] presented multiple designs and methods to control a system using ADRC. The authors presented a precise design for ADRC to control a second order system. ESO and TD parameters were chosen as functions of the system bandwidth. A discrete implementation of ADRC was presented and gain matrices of its components were shown. The authors claimed that a generalized form of ADRC to control an n-th order system can be presented in the same manner. A new design called Predictive Active Disturbance Rejection Controller was proposed. The basic design principles for this controller were explained and a compatible Reduced order ESO (RESO) was used.

Peng C., et al. in [8] proposed an ADRC trajectory tracking control based on Particle Swarm Optimization (PSO) for a quad-rotor (QR) system. PSO is used to adjust and optimize parameters of ADRC so as to enhance the system design efficiency. A dynamic model for QR was presented and the three impediments

of the ADRC controller: TD, ESO and NLSEF were proposed. The TD unit in the altitude z channel was designed based on fast discrete tracking differentiator method as in [9]. New equation based forms for ESO and NLSEF were presented. The quad-rotor closed loop system was divided into four independent channels by ADRC. Controller optimization using PSO was applied to the four channels respectively.

Wang Ch., et al. in [10] addressed the problem of steering Underactuated Surface Vessels (USVs) along given paths. ADRC technique was proposed to eliminate dynamic uncertainty and ocean disturbances induced by unknown wind, waves, and ocean currents. Second-order ADRC was used to control the equivalent rudder angle by the propeller rate. Blocks based MATLAB Simulation was used to simulate the proposed USV system, ADRC parameter optimization using Chaos PSO (CPSO) for clamping force control system in order to achieve higher control precision, designed equations for TD, ESO, and NLSEF, CPSO algorithm to minimize a fitness function J which is a function of time, error, and control signal.

Simulation results showed better performance when using ADRC with PSO optimization algorithm in comparison with using ADRC alone.

ADRC Approaches. In this section, two approaches of second-order ADRC controller are presented.

Linear ADRC. The linear ADRC is described based on the American patent [7]. Fig. 1 shows the contents of a linear ADRC structure.

Vo

LTD

V1 /о el

V2

«2

PD

control

щ л

-Hgb Ma plant

z2

z3

LESO

Fig. 1. Linear ADRC structure

This structure consists of controller, linear estimator (Linear Extended State Observer (LESO)), linear tracking differentiator (LTD) and a disturbance rejection scheme.

Linear Tracking Differentiator (LTD). Using LTD, the reference signal is smoothed and output signals are generated to track the reference signal and its difference. The designed algorithm is as follows:

(1)

vi=v2;

1>2 =-h(v1-vo)-k2v2,

where Vo is the reference signal, Vi is the tracking signal of the system and v2 is the differential signal of the system. Here k\, k2 are adjustable parameters

of LTD. When k\ = r2; k2 = 2r; r > 0, the transition process is not overshoot, and transition time is about To =71 r, r is speed factor.

Linear Extended State Observer. The idea behind ESO is capturing information about the generalized disturbances (uncertainties and external

A A

disturbances /) and the internal dynamics of the system (y, y). System dynamics can be expressed in the following general form:

y = g(t> y> y) + au + w, (2)

where y is the output signal, u is the control input, g(.) is the dynamic of the plant (including the unknown dynamics), w is the external disturbance, and a is the system coefficient. All parts of this dynamics (g(.), a, w) are usually not precisely known. By combining the external and the internal disturbances in one function /(.) the dynamic can be written as:

y = f(t, y, y, w) + au. Equation (3) can be written in state space in the form:

xi=x2; x2= f + au; y = x 1.

The total disturbance is augmented with the states in the following way:

xi=x2; x2=xi+ au; x3=f(t, Xi, x2, w); y = x 1.

This can be expressed in the vector-matrix form as:

x = Axx + Bxu + Exf; y = Cxx,

where

(3)

(4)

(5)

(6)

0 1 0 0 0

Ах = 0 0 1 . вх = а , Сх=[ 1 0 0], Ех = 0

0 0 0 0 1

(7)

A linear ESO (LESO) can be used here to estimate the states x\, X2, x$. This way, LESO can be designed to be:

¿1 = z2 + aie;

¿2 = Z3 + au + a2e; (8)

¿3 = + ase,

where Z\, Z2 and Z3 are the approximated values of states x\, X2 and X3 respectively, ai, a2 and 0C3 are the observer gains, and e = y - Z\ is the estimated error of x\, a is an estimated value for a in equation (2) and it can be chosen empirically in this structure.

__A A

The estimated variables (Z\ = y, z2 = y, Z3 = /) besides the approximated

value a are then used to eliminate the disturbance and control the system as shown in Fig. 1.

Disturbance rejection scheme. The disturbance rejection scheme can be defined as:

A

»0 -Z3 up -f

U =-:-= —:—, (9)

a a

where u is the rejection scheme output and uq is the control part output. By replacing u by its estimated value in equation (2):

У = /(•) + a

f "f\ uq- J

A

a

(10)

If the estimator is good enough to consider that ä » a and f ~ f, then the plant can be written as a second order integrator:

y *u0. (11)

Feedback controller. If a proportional-derivative (PD) controller is used as a feedback controller, the control signal uo can be written in the form

u0(t) = Kp(v0 -y) + K/y, (12)

where vo is the reference input. PD gains can be chosen as

Kp=w2CL, Kd=-l%wCL, (13)

where wcl is the desired closed loop pole and is the desired damping coefficient of the closed loop.

For observing poles, Weso should be placed n times to the left of the closed loop pole to ensure that the dynamics of the observer are faster enough than the dynamics of the controlled plant, where n e [3,10], i.e.,

■WESO = nWcL- (14)

For simplicity, all the observer poles are placed in one location. That means, the characteristic equation of the observer is:

DQ0 = (X + WESO)3 + 3WES0^2 + + WggQ, (15)

where X represents Eigen values and ai, a2, 0L3 are calculated by solving the equation

D(X) = \sI-Ax+LCX\, (16)

where s is Laplace frequency parameter and

(17)

"l 0 0 ai

0 1 0 , L = «2

0 0 1 аз

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

1 =

As a result, the estimator gains can be chosen as

ai =3wESo, a2 = 3w|so, a3 = wlso-

(18)

Nonlinear ADRC. Fig. 2 shows a typical second order nonlinear ADRC scheme.

d{t)

Fig. 2. Nonlinear ADRC structure

Similar to linear structure, ADRC consists of two main loops: the feedback loop and the estimation loop. It has three main components: the Tracking Differential (TD), ESO and a NLSEF.

Tracking Differential. TD is used for arranging the transient process, which can be expressed as:

vi=v2;

V2 = /td(vi - v0, v2, r) = -r sign

Vi - Vo +

V2\V2\ 2 r

(19)

where vo is the input signal, v\ is the tracking signal of vo, v2 is the differential signal of Vi, r is the parameter which determines the tracking speed.

The nonlinear function /td is proposed in [11] to provide the fastest tracking of the reference signal and its derivative subject to the acceleration limit of r.

The parameter r is an application dependent and it is set accordingly to speed up or slow down the transient profile.

Extended State Observer. As for the linear structure, ESO is employed to monitor the output and predict the real-time state of the plant. It can be expressed as follows [12]:

e = z\-y, ¿1 = z2 + die;

. | ai1/2 . . . (20)

z2 = z$ + a2 \e\ signe + au;

■ I *|l/4 . A

Z3 = 0C3 \e\ sign e,

where y is the system output, Z\ is the tracking signal of y, e is the estimated error, z2 is the differential signal of Z\, z$ is the tracking signal of general disturbances, 04, a2 and (X3 are estimator gains, u is the control input, and a is the system coefficient.

Nonlinear State Error Feedback and disturbance rejection process. NLSEF is a nonlinear control strategy which can enhance the precision for control system. It can be expressed as follows:

¿1 = vi - Zi;

e2=v2- z2;

"o = Pi fnl (h, Yi, ti) + p2 fnl {e2, J2, r|); (21)

Z3

u = uo--,

a

where e\ is the estimated system error, e2 is the estimated differential of system error, Pi and P2 are gain coefficients, and z$la is the total disturbance compensation of the system; /„/ should achieve good control effect and high frequency tremor, and it can be chosen as [12]:

fnlie, у, Л) =

\e\ < r|;

r|Y_1' ' ' " (22)

|e|Ysigne, \e\ > rj.

Discrete implementation ofADRC. ADRC is implemented in discrete time. Discretization can be done using Euler or zero order hold (ZOH).

Discretization of linear ADRC. In discrete implementation of the linear ADRC, the ESO pole is determined in Z domain using the following equation [4]:

Dz(z) = I zl - Ф + LpCx I = (z- a)3, a =

(23)

where h is sampling time, Lp is the observer gain matrix, Cx is the plant observation matrix, and O is the plant discrete state matrix. In [7], authors proposed two main discrete ESO.

The first one is Predictive Discrete Extended State Observer (PDESO), for which the following equations are applied:

z(k +1) = Фz(k) + Гu(k) + Lp[y(k) - Zi(k)]; z(k) = z(k).

(24)

The second one is Current Discrete Extended State Observer (CDESO), for which the following equations are applied:

z(k + 1) = Фz(k) + Tu(k) + Lp[y(k) - Zi(fc)];

z(k) = z(k) + Lc[y(k)-zm.

The different matrices for Euler Discretization are

(25)

Ф =

1 h 0 0 1 h О 0 1

Lp =

Г =

О âh О

Lr =

3(1 - а) h

(1 - а)3 h2

1 - а3 (1 - а)2 (2 + а)

h

(1 - ос)3

h2

(26)

These matrices for ZOH Discretization are

h2

Ф =

1 h

О 1 О О

2 h

Lp =

3(1-а) (1-а)2(5 + а) 2 h (1-а)3 h2

(27)

г =

a-

h2

2 ah

0

Lc =

1-a3 , (l-a)2(l + a) 2 h

(1 - a)3/T2

(27)

(28)

Discretization of nonlinear ADRC. Authors in [11] proposed a discrete form of TD in order to avoid introducing considerable numerical errors in a discrete-time implementation. This discrete form is obtained as follows:

vi =vi+hv2; v2=v2+hv, |v|<r,

where v is written in the form

v = fhan(vi - v0,v2,r,h). (29)

Here Vo is the input signal, V\ is the tracking signal of Vo, V2 is the differential signal of vi, r is the parameter which determines the tracking speed, is the sampling time and fhan is:

d = hr2; ao=hv2; yo=vi+ao; a\ = yjd(d + 8|jy|); a2 = flo + signy0

{a\ — d)

Sy —

sign(7o + d) - sign(yo - d)

af = (öo + Уо ~ a-2 )sy + a2;

(30)

5a =

sign(af + d) - sign(af - d)

fhan — *

aj_ d

sign a f

sa - rsignaf.

Authors in [11] proposed a discrete implementation of the extended state observer ESO as follows:

e = Zi - y;

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

fe = /esoCs, 0.5, r|); fei = /Eso(e, 0.25, r|);

Z\ = Z\ + hz2 + aie;

z2=z2+ h{z3 + au) + a2fe; (31)

z3 = z3 + a3/ei;

|ef signe, \e\>r\,

feso(e, С. Л) =

Л >0,

where h is the sampling time and C, is the nonlinearity factor.

Finally, control and rejection scheme parts are the same as for the continuous implementation, i.e.,

h = vi - z\\

e2 =v2 - z2;

"o = Pi/eso (¿i> Yi, r|) + P2/eso(£2, Y2. ^

Z3

U = Uq--.

a

Discrete nonlinear ADRC parameters selection. There are two parameters in TD, [r, h], six in ESO, [ai, a2, 013, ¿¡1, C,2, a] and five in NLSEF, [Pi, P2.Y1.Y2.Tl]- In general there are 13 parameters, these parameters are determined empirically. ADRC suffers from very limited application due to its complicated parameters adjustment. To overcome this problem, key parameters of ADRC are usually adjusted in two stages. In the first stage, parameters are roughly adjusted to estimate the ballpark ranges. In the second stage, further optimization of the parameters with specific algorithm is done to determine the precise adjustment. Genetic Algorithm (GA), and Particle Swarm Optimization (PSO) have been widely used in the optimization process to determine optimal ADRC parameters. In order to facilitate estimation of the parameters, they can be divided into key parameters and common parameters [8]. Key parameters are [ai, a2, 0C3, Pi, P2], on which optimization process is applied.

Case Study: ADRC for quarter-car active suspension system. Comfort is a very important issue for passengers when they travel from one place to another. During traveling, a car may encounter various vibrations and shocks from ground. These vibrations may cause problems to passengers. To overcome these problems, loads from ground should be damped out. Automotive suspension systems are intended to absorb and decrease the shocks and vibrations transferred from the ground to the passengers as well as to the vehicle body. Traditionally, passive suspension system is used to absorb vibrations in many different types of vehicles. Active suspension systems separate actuators to apply controlled forces and provide better ride.

The suspension system shown in Fig. 3 represents vehicle system at each wheel. It consists of a spring ks, a damper bs and a hydraulic actuator Fa. In addition, kt represents tire stiffness and bt represents tire damping. The effective vehicle body mass is expressed by Ms (sprung mass), and Mu (unsprung mass) represents the effective mass for the wheel and the axle. States xu, xs represent the vertical displacements from the static equilibrium for Mu

Fig. 3. Quarter-car active suspension system

and Ms respectively. The road profile is represented by xr. The suspension travel xs —xu is measured and compared to the set point (vo = 0).

The required actuator force is determined by the controller to eliminate errors, and thus, to reduce the vehicle oscillations. The actuator can provide a maximum force of 3 000 N.

Table 1 shows the approximated values of the system parameters.

Table 1

Approximated values of the quarter-car system parameters

Parameter Unit Value Parameter Unit Value

Ms kg 300 k N/m 180 000

М„ kg 50 bs N-s/m 1200

к$ N/m 18 000 bt N-s/m 0

Vibration model of the sprung and unsprung masses can be expressed by the following dynamic equations:

Muxu = ks(xs — xu) + bs(x s — xu) — kt(xu — xr) — bt{xu — xr) — Fa;

- , / x , -x (33)

Msxs = -ks(xs - xu) - bs(xs - xu) + Fa.

However, Mohite, et al. in [13] proposed a nonlinear model for the passive quarter-car system. In their work, tire stiffness kt for the linear system can vary as kti, kt2, kt3 for considering of nonlinearity. In the same manner, ks can vary as ks\, ks2, and bs as bsi, bs2. In this way, equations of motion can be expressed in the following nonlinear form:

Muxu = ksi(xs-xu) + ks2(xs-xuf + bsi(xs-xu) + bs2(xs-xu)2 -- kn(xu-xr)-kt2(xu-xr)2 + kt3(xu-xr)3 - bt{xu-xr)-Fa; (34) Msxs = -ksi(xs-xu) - ks2(xs-xu)3 - bsl(xs-xu) - bs2(xs-xu)2 + Fa.

Table 2 shows the approximated values of the nonlinear system parameters.

Table 2

Approximated values of the quarter-car nonlinear system parameters

Parameter Unit Value Parameter Unit Value

м5 kg 295 kti N/m 42 509

Ми kg 39 kt3 N/m 22 875

ksi N/m 15 302 bsi N-s/m 3 482

kS2 N/m 2 728 bs2 N-s/m 580

kl N/m 60 063 bt N-s/m 0

Control of quarter-car model using linear ADRC. Simulink blocks have been built to represent linear and nonlinear models for simulating control of quarter-car active suspension system. LTD, PD and LESO blocks have been built as in the first section. Closed loop pole wql is usually chosen empirically as a trade-off between states approximation convergence speed and the influence of both sampling time and noise [14]. However; wql =1-7 is found to be a good choice. The estimated value of a from equation (2) is chosen empirically to be a = 0.001.

Linear ADRC is simulated for studying the control of the linear and nonlinear models of quarter-car system. Moreover; using this simulation, a comparison is made between this system performance and the performance of the passive system. Passive suspension system is simulated by applying a constant input Fa = 0 when the controller is not activated. Two types of road disturbance are tested: a double bump disturbance and random disturbance. The double bump is described by equation, where w is the sudden bump, hb is the bump

height, to is the bump initiating time, xt is a time delay, V is the car speed, fi is an arbitrary constant and k is the disturbance wavelength

0.5 fa

1 - cos

2nV(t -t0)

\

/

, to<t<tQ+—;

w(t) = -

0.5 ц/zfe 0 otherwise.

' (2nV(t -t0 1 - cos -

l A.

, to + xt < tf <to + xt + —; (35)

Fig. 4 shows the bump signal that corresponds to the following parameters V = 45 km/h; fa = 0.1 m; to = 4 s; (J, = 1 m and xt = 1 s.

0.10

4.0 4.5 5.0 5.5

Fig. 4. Double bump road disturbance

6.0 Time, s

Random road disturbance is randomly filtered as in [15] which has clear physical meaning and easy computing character:

q(t) = -2nfoq(t) + 2nnoy]Gq(no)v w(t),

(36)

where fo is the filter lower-cut-off frequency, q(t) is the random road input signal, Gq(no) is the road roughness coefficient and w(t) is a Gaussian white noise. The following parameters are considered: an integration filter (/o = 0), the vehicle speed is v = 20 m/s driving is on a class D road for which Gq(no) =

= 1 024 • 10-6 (see [15]).

Fig. 5 shows both responses of passive and linear ADRC control systems to a double bump road disturbance which was shown in Fig. 4. It is obvious that linear ADRC compensates for disturbances and leads the system to enter the stability region faster.

0.10 0.05

a

1 o

I

a

& -0.05 Q

-0.10 -0.15

3.0 3.5 4.0 4.5 5.0 5.5 6.0 6.5 Time, s Fig. 5. ADRC and Passive control systems responses for double bump disturbance

Fig. 6 shows the estimated values for the output y during the process. Fig. 7 shows the estimated values for the output rate y. It is clear that there is a good estimation for the displacement but a poor one for its rate. To improve rate estimation, the output acceleration could be taken as a second output. However, it's not the subject of this research.

0.10 0.08 0.06 a o.o4 | 0.02

& -0.02 p

-0.04 -0.06 -0.08 -0.10

3.0 3.5 4.0 4.5 5.0 5.5 6.0 6.5 Time, s Fig. 6. Linear ADRC ESO estimation for the output signal

Fig. 8 shows both responses of passive and linear ADRC control systems to the random road disturbance described in equation (36). The ability of linear ADRC to compensate for disturbances is clear. To verify the control efficiency for random disturbance, root mean square (RMS) values are calculated for the passive and the controlled system. Using MATLAB, the following values are found for RMS:

RMS for Passive system output is 0.0118;

RMS for linear ADRC controlled system output is 0.009.

This means that the disturbance density is reduced by 25 %.

5 4

3

2 -

$ 1

* 0 -1

-2 -3

-4 3

---Linear ADRC

-Passive

—I

0 3.5 4.0 4.5 5.0 5.5 6.0 6.5 Time, s Fig. 7. Linear ADRC ESO estimation for the output signal rate

0.03

---Linear ADRC

-Passive

-0.03

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Time, s

1 2 3 4 5 6 7 8 Fig. 8. Linear ADRC for random road disturbance

Now the ability of linear ADRC to compensate for dynamic uncertainties could be tested by simulating it for the nonlinear quarter-car system.

Here, the nonlinear equations (34) are simulated in the presence of an actuator model in the form of equation (37), where t = 70 s is the actuator's pole and T = 0.015 s is the actuator's delay.

-T

z -u, -1 500 < u <1 500;

У act ='

1 + T

-1 500, 1 500,

и < -1 500; и > 1 500.

(37)

It can be seen from Fig. 9 that linear ADRC controller slightly compensates for the different nonlinearities of the system. This compensation could be made better by retuning the ADRC parameters or moving to a nonlinear ADRC.

Fig. 9. Linear ADRC vs passive system performance for nonlinear quarter-car system

Control of quarter-car model using nonlinear ADRC. In this section the discrete version of the nonlinear ADRC is used. Parameters are chosen as in Table 3.

Table 3

Tuning parameters for nonlinear ADRC

Key parameters General parameters

ai 3weso г 20

a2 3w|so h 0.01

аз WESO Yi 0.4

ßi 1 Ъ 0.6

ß2 50 Ci 0.5

- - С2 0.25

- - ц 1.5

In Table 3 ai, a2, 0C3 are chosen similarly as in the linear ADRC.

Fig. 10 shows responses of both the passive and the nonlinear ADRC control systems to the double bump road disturbance shown in Fig. 4. In this figure, the ability of nonlinear ADRC to compensate for disturbances and damp them faster is clear.

g a

и ä Ъ

■ i—I

Q

0. 0. 0. 0.

-0. -0. -0. -0. -0. -0.

08

06

04

02

0

02

04

06

08

10

12 3.0

Fig.

i. ---Nonlinear ADRC -Passive

1 1 1 1 1 1 1

3.5

4.0

4.5

5.0

5.5

6.0

6.5 Time, s

10. Double bump response of quarter-car system controlled by nonlinear ADRC

Fig. 11 shows responses of both the passive and the nonlinear ADRC control systems to the random road disturbance described in equation (36). RMS values for passive and nonlinear ADRC systems outputs are: RMS (for Passive system output) is 0.0118; RMS (for nonlinear ADRC controlled system output) is 0.008. This means that disturbance density is reduced by 33 %.

0.03

9 Time, s

Fig. 11. Random road disturbance response of quarter-car system controlled

by nonlinear ADRC

It can be noticed that the output RMS value here is a little better than the RMS value for linear ADRC. This is due to the nonlinear functions used inESO.

Again, the ability of the nonlinear ADRC to compensate for dynamic uncertainties could be tested by simulation in the nonlinear quarter-car system.

Fig. 12 shows that nonlinear ADRC controller compensates for the different nonlinearities of the system. RMS values for passive and nonlinear ADRC systems outputs are:

RMS (for Passive system output) is 0.0081;

RMS (for nonlinear ADRC controlled system output) is 0.006.

-0.020 -1-1-1-1-1-1-1-1-1-

0123 456789 Time, s

Fig. 12. Nonlinear ADRC vs passive system performance for nonlinear quarter-car system

This means that the disturbance density is reduced by 30 % and the nonlinear ADRC shows better performance than the linear one.

Conclusion. In this work, two main approaches of the second order ADRC are presented: linear and nonlinear ADRC. Both approaches consist of two loops: feedback loop and estimation loop. The three main components: controller, estimator and disturbance rejection scheme for each approach are presented. Timing and selection parameters for both approaches are discussed. Discrete implementation forms for both approaches are presented. A simple linear and nonlinear quarter-car models are used as a study case. Active suspension system controlled by ADRC is tested against passive system. Simulation results show a good ADRC performance and efficiency against internal and external disturbances.

REFERENCES

[1] Qing Z., Gao Z. On practical applications of active disturbance rejection control.

Proc. 29th Chinese Control Conf., 2010.

Available at: https://ieeexplore.ieee.org/document/5572922

0.020

-0.015

и

Translated by Authors

[2] Qian L., Li D., Tan W. Design of multi-loop ADRC controllers based on the effective open-loop transfer function method. Proc. 33rd Chinese Control Conf., 2014.

DOI: https://doi.Org/10.l 109/ChiCC.2014.6895546

[3] Xue W., Bai W., Yang S., et al. ADRC with adaptive extended state observer and its application to air-fuel ratio control in gasoline engines. IEEE Trans. Ind. Electron., 2015, vol. 62, no. 9, pp. 5847-5857. DOI: https://doi.org/10.1109/TIE.2015.2435004

[4] Li F., Zhang Z., Armaou A., et al. Study on ADRC parameter optimization using CPSO for clamping force control system. Math. Probl. Eng., 2018, vol. 2018, art. 2159274. DOI: https://doi.Org/10.l 155/2018/2159274

[5] Zhu E., Pang J., Sun N., et al. Airship horizontal trajectory tracking control based on active disturbance rejection control (ADRC). Nonlinear Dynamics, 2014, vol. 75, pp. 725-734.

[6] Li X., Wang S., Wang X., et al. Permanent magnet brushless motor control based on ADRC. MATEC Web Conf., 2016, vol. 40, art. 08003.

DOI: https://doi.org/10.1051/matecconf/20164008003

[7] Gao Z., Tian G. Extended active disturbance rejection controller. Patent US 8180464. Appl. 20.08.2008, publ. 01.01.2009.

[8] Peng C., Tian Y., Bai Y., et al. ADRC trajectory tracking control based on PSO algorithm for a quad-rotor. ICIEA, 2013, pp. 800-805.

DOI: https://doi.Org/10.l 109/ICIEA.2013.6566476

[9] Han J.Q. From PID technique to active disturbance rejection control technique. Control Engineering of China, 2002, vol. 9, no. 3, pp. 13-18.

[10] Wang Ch., Zhang H., Yu Y. USV trajectory tracking control system based on ADRC. CAC, 2017, pp. 7534-7538. DOI: https://doi.org/10.1109/CAC.2017.8244141

[11] Han J. From PID to active disturbance rejection control. IEEE Trans. Ind. Electron., 2009, vol. 56, no. 3, pp. 900-906. DOI: https://doi.org/10.1109/TIE.2008.2011621

[12] Zhou X., Yang C., Zhao B., et al. A high-precision control scheme based on active disturbance rejection control for a three-axis inertially stabilized platform for aerial remote sensing applications. /. Sens., 2018, vol. 2018, art. 7295852.

DOI: https://doi.Org/10.l 155/2018/7295852

[13] Mohite A.G., Mitra A.C. Development of linear and non-linear vehicle suspension model. Mater. Today, 2018, vol. 5, no. 2-1, pp. 4317-4326.

DOI: https://doi.Org/10.1016/j.matpr.2017.l 1.697

[14] Madoriski R., Herman P. An experimental verification of ADRC robustness on a cross-coupled aerodynamical system. IEEE Int. Symp. Industrial Electronics, 2018, pp. 859-863. DOI: https://doi.org/10.1109/ISIE.2011.5984271

[15] Zhou Q. Research and simulation on new active suspension control system. Master Thesis. Lehigh University, 2013.

Alhelou M. — Post-Graduate Student, Department of Automatic Control Systems, Bauman Moscow State Technical University (2-ya Baumanskaya ul. 5, str. 1, Moscow, 105005 Russian Federation).

В Издательстве МГТУ им. Н.Э. Баумана вышло в свет учебное пособие авторов Е.А. Микрина, М.В. Михайлова

«Навигация космических аппаратов по измерениям от глобальных спутниковых навигационных систем»

Рассмотрены вопросы проектирования и разработки сложных многофункциональных систем космической навигации на базе глобальных спутниковых навигационных систем для широкого класса низкоорбитальных, высокоорбитальных и высокоэллиптических космических аппаратов, а также круг вопросов, связанных с созданием бортовых средств навигации для автономного определения орбиты космического аппарата.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

По вопросам приобретения обращайтесь:

105005, Москва, 2-я Бауманская ул., д. 5, стр. 1

+7 (499) 263-60-45

[email protected]

Ьйрв^/ЬтвШ.ргевв

Gavrilov A.I. — Cand. Sc. (Eng.), Assoc. Professor, Department of Automatic Control Systems, Bauman Moscow State Technical University (2-ya Baumanskaya ul. 5, str. 1, Moscow, 105005 Russian Federation).

Please cite this article as:

Alhelou M., Gavrilov A.I. Synthesis of active disturbance rejection control. Herald of the Bauman Moscow State Technical University, Series Instrument Engineering, 2020, no. 4, pp. 22-41. DOI: https://doi.org/10.18698/0236-3933-2020-4-22-41

E.A. Микрин, M.B. Михайлов

Навигация космических аппаратов по измерениям от глобальны* спутниковых навигационных систем

i Надоели баннеры? Вы всегда можете отключить рекламу.