In the work, a PID design method with desired settling time is proposed for plants with
transfer functions of the following types: first-order system, second-order system and integral
first-order system. The method provides better performances than both the magnitude optimum
method and the symmetric optimum method in the sense of overshoot and settling time. In
addition, the settling time can be achieved as expected for a larger frequency band. Some
examples for DC motor control are used to illustrate the method.
Future works will focus on robustness, effects of model error and other types of transfer
functions, may be higher-order system with time delay
11 trang |
Chia sẻ: honghp95 | Lượt xem: 534 | Lượt tải: 0
Bạn đang xem nội dung tài liệu Settling time assignment with pid - Nguyen Doan Phuoc, để tải tài liệu về máy bạn click vào nút DOWNLOAD ở trên
Vietnam Journal of Science and Technology 55 (3) (2017) 357-367
DOI: 10.15625/2525-2518/55/3/8808
SETTLING TIME ASSIGNMENT WITH PID
Nguyen Doan Phuoc, Nguyen Hoai Nam
Hanoi University of Science and Technology, 1 Dai Co Viet Street, Hanoi
*Email: nam.nguyenhoai@hust.edu.vn
Received: 24 October 2016; Accepted for publication: 21 February 2017
ABSTRACT
In this paper, the settling time of closed-loop systems controlled with PID, which is
established by using magnitude and symmetric optimum methods, is analyzed. These methods
have some limitations such as limited frequency bands, existence of overshoot and larger settling
time. Thus, our aim is to improve these performances of the closedloop system for a certain class
of transfer functions. The proposed method will be tested on a DC motor by simulation. It can be
concluded that the method provides a larger range of frequency, smaller settling time and no
overshoot.
Keywords: PID controller, magnitude optimum criterion, symmetric optimum criterion, first-
order system, first-order plus integral system.
1. INTRODUCTION
Many PID tuning methods have been developed since 1940s. These methods can be
classified into groups: practical approach [1 - 4]; optimal approach [5 - 8]; gain and phase
margin specification based approach [9 - 13] and others [14 - 18]. However, there have not been
yet any method that concurrently dealing with both the settling time and the overshoot of the
closed-loop system.
It is known that when plants, having transfer functions of first-order, second-order and
third-order systems the magnitude optimum method based parameter determination , p Ik T and
DT [5 - 8] for the PID controller seems to be one of the most effective methods, validated in
reality. However, there has been no fundament until now to conclude that the magnitude
optimum method is the only one to provide PID controller parameters to make the closed-loop
system stable and the settling time short in the sense that the closed-loop system with transfer
function
( ) ( )( )
1 ( ) ( )= +c
R s G s
G s
R s G s
has
Nguyen Doan Phuoc, Nguyen Hoai Nam
358
{ } ( ) 1R+Ω = ∈ ≈cG jω ω (1)
in small range of frequency, which is large enough.
A method, proposed in this work, has an approach that the closed-loop system, obtained as
follows,
/
1( )
1
=
+
cG s
T s
(2)
has as small /T as possible (when /T is smaller, Ω is larger), which will demonstrate the
conclusion.
2. PID CONTROLLER DESIGN
2.1. Control of First-Order Systems
Firstly, let’s assume that the plant has a form of first-order system and a PI controller is
used. Then, the transfer function of the plant is represented as follows
( )
1
=
+
k
G s
Ts
and the PI controller is
( )11( ) 1 + = + =
p I
p
I I
k T s
R s k
T s T s
Thus, the open-loop transfer function is
( )
( )
1( ) ( ) ( )
1
+
= =
+
p I
o
I
k k T s
G s R s G s
T s Ts
(3)
and the closed-loop transfer function is
2
( ) 1( )
1 ( ) 11 1
+
= =
+
+ + +
o I
c
o I
I
p p
G s T s
G s
G s TT
s T s
k k k k
(4)
Since the first-order system (2) with unit gain is always stable, the steady state error is zero and
the settling time is smaller if the time constant is shorter, in this case the two parameters pk and
IT are determined such that the closed-loop transfer function (4) has the expected form as (2).
In order to do that, pk and IT are chosen to cancel out a zero by a pole. The zero is 1= − Is T .
This is equivalent to
2 11 1 0
+ + + =
I
I
p p
TT
s T s
k k k k
as
1
= −
I
s
T
or
( )0 1= − + +I p I pT T k k T k k
Settling time assignment with PID
359
Thus .=IT T
Next, the closed-loop system becomes
/
2
1 1( )
111 1
+
= =
+
+ + +
I
c
I
I
p p
T s
G s
T sTT
s T s
k k k k
with the selected parameters as above, if the time constant /T satisfies the following equation:
( )( ) ( )( )
( )
/ / 2
2
2
11 1 1 1 1 1
1
1
+ + = + + = + + +
+
= + +
I
I I
p p
p
p p
TT
T s T s T s Ts s T s
k k k k
T k kT
s s
k k k k
⇔ ( ) ( )2/ / 2 2 11 1++ + + = + +p
p p
T k kT
T T s T Ts s s
k k k k
⇔
( )/ 1+
= − =
p
p p
T k k T
T T
k k k k
This implies
( )2
1 1( )
1 1 1
+
= =
+ + + +
I
c
p I I p
p
T s
G s
Tk k TTs T k k s s
k k
Hence, in the case of =IT T the bigger pk is chosen, the smaller settling time of closed-
loop system will be. Furthermore, the closed-loop system is stable and has no overshoot.
Obviously, in comparison to the performance of the magnitude optimum criterion method,
the proposed PI controller provides a more flexible band of frequency (Ω). When the parameter
pk is bigger, the frequency range Ω becomes wider:
lim ( ) 1,
→∞
= ∀
p
c
k
G s s
Therefore, the desired settling time 5%T of closed-loop system can be assigned by choosing
5%
ln 20
=p
T
k
kT
which is clearly obtained from its unit step response
/( ) 1 −= − t Th t e or /5%0.95 1 −= − T Te
2.2. Control of Second-Order Systems
When the plant has the form of a second-order system, we use a PID controller. The plant model
Nguyen Doan Phuoc, Nguyen Hoai Nam
360
and the controller are represented, respectively, as follows
− Plant: ( )( )1 2( ) 1 1= + +
k
G s
Ts T s
− PID controller:
( )( )1 11( ) 1 + + = + + =
p A B
p D
I I
k T s T s
R s k T s
T s T s
where + =A B IT T T and =A B I DT T TT , it means that whenever AT and BT are known we can
obtain IT and DT . The open-loop system is
( )( )
( )( )
( ) ( )
( ) ( )
/
1 2 1 2
1 1 1 1( ) ( ) ( )
1 1 1 1
+ + + +
= = =
+ + + +
p A B p A B
o
I A
k k T s T s k k T s T s
G s R s G s
T s Ts T s T s Ts T s
(5)
where
/
=
p A
p
I
k T
k
T
Thus, if it is chosen that 2=BT T then the transfer function (5) has the same form as (3) shown
in the first case
( )
( )
/
1
1( )
1
+
=
+
p A
o
A
k k T s
G s
T s Ts
where the parameters IT and pk are replaced by the parameters AT and
/
pk , respectively. So,
with the set of parameters in the case 1: 1=AT T and
/
pk arbitrary, the closed-loop transfer
function becomes a first-order system with unity gain and the time constant
/ 1 1
/= = =
I I
p A pp
T TT T
T
kk T kkk k
and therefore, the bigger pk is chosen, the smaller its time constant
/T will be. Thus, we have
1 2
1 2
1 2
and = + = + = =
+
A B
I A B D
I
T T TT
T T T T T T
T T T
In summary, the larger pk is optionally selected the smaller the settling time associated with
/T
will be. The assignment of desired settling time 5%T for the closed-loop system can be easily
realized with
( )1 2
5%
ln 20+
=p
T T
k
kT
In addition, there is no overshoot in the closed-loop system.
2.3. Extended Method to Integral First-Order Systems
Settling time assignment with PID
361
For the basic control theory of linear systems, it is well known that the symmetric optimum
method is applied to find the parameters , p Ik T and DT of the PID controller
1( ) 1 = + +
p D
I
R s k T s
T s
(6)
for integral first-order systems and integral second-order systems, in order to make the closed-
loop system stable and keep the phase margin larger. However, the disadvantage of the method is
big overshoot. Although the overshoot can be decreased with the usage of an input filter, this
causes an additional cost.
The proposed method for determination of PID parameters later will guarantee that the
closed-loop system has no overshoot and arbitrary small settling time. It is similar to first-order
or second-order systems, the idea of this approach is also to find parameters , p Ik T and DT of
the PID controller for the integral first-order plant
( )( ) 1= +
k
G s
s Ts
(7)
such that the closed-loop system has the following form:
( )( )1 2
1( )
1 1
+
=
+ +
t
c
m m
Ts
G s
T s T s
(8)
Consider a closed-loop system, which includes an integral first-order system (7) and a PID
controller (6)
( )( )1 11( ) 1 + + = + + =
p A B
p D
I I
k T s T s
R s k T s
T s T s
with + =A B IT T T and =A B I DT T TT . The open-loop system is as follows
( )( )
( )2
1 1( ) ( ) ( )
1
+ +
= =
+
p A B
o
I
kk T s T s
G s R s G s
T s Ts
Thus, if =BT T then = +I AT T T and the closed-loop system is
( )
( ) ( )2
1( )( )
1 ( ) 1
+
= =
+ + + +
p Ah
c
h A p B
kk T sG s
G s
G s T T s kk T s
So, the closed-loop system will be the same as the expected form (8) if there the following
equation must be satisfied
( )
( ) ( ) ( )( )2 1 2
1 1
,
1 11
+ +
= ∀
+ ++ + +
p A t
m mA p A
kk T s Ts
s
T s T sT T s kk T s
⇔ ( ) ( )( ) ( ) ( ) ( )21 21 1 1 1 1 , + + + = + + + + ∀ p A m m A p A tkk T s T s T s T T s kk T s Ts s
Nguyen Doan Phuoc, Nguyen Hoai Nam
362
Controlled
Object PID
Input
filter
Figure 1. Control with PID and input filter.
⇔ ( )
( )
1 2 1 2
1 2
1 2
+
+ + = +
+ = + +
+
=
A
A m m m m A t
p
A t A m m
A t
m m A
p
T T
T T T T T T T
kk
T T T T T
T T T
T T T
kk
⇔
( )
1 2
1 2
1 2
+
=
= +
+
=
A
m m
p
t m m
A t
m m A
p
T T
T T
kk
T T T
T T T
T T T
kk
⇔
1 2
1 2
+
=
= = +
A
m m
p
A t m m
T T
T T
kk
T T T T
This implies that
( )
1 2
1 2 1 2
1 2
1 2
1 2
+ + +
= =
= + = + +
+
= =
+ +
A m m
p
m m m m
I A m m
m mA B
D
I m m
T T T T T
k
kT T kT T
T T T T T T
T T TT T
T
T T T T
(9)
It can be seen that with the above equations, 1 mT and 2mT are freely chosen. The expected
transfer function (8) becomes
( )
( )( )
1 2
1 2
1( )
1 1
+ +
=
+ +
m m
c
m m
T T s
G s
T s T s
Thus, with an addition of an input filter as illustrated in Fig.1:
( )1 2
1( )
1
=
+ +m m
M s
T T s
to obtain the closed-loop system:
( )( )1 2
1( ) ( )
1 1
=
+ +c m m
G s M s
T s T s
we always can choose 1mT and 2mT such that this closed-loop system has a settling time 5%T as
desired. Obviously, there are two possible combinations of time constants 1mT and 2.mT The
first case is 1 2= =m m mT T T and the second
one is 1 2≠m mT T .
Now some steps, about how to choose these
time constants such that the desired settling
time can be achieved, will be shown. Firstly,
Settling time assignment with PID
363
let’s consider the case one. With the usage of input filter as shown in Fig. 1, the closed-loop
transfer function is as follows
( )2
1( )
1
=
+
c
m
G s
T s
So, the unit step response is
/ /( ) 1 − −= − −m mt T t T
m
t
h t e e
T
Thus, the desired settling time 5%T will be the solution to the following equation
/ / 0.05 0− −+ − =m mt T t T
m
t
e e
T
Let 5%= mx T T . Then
( ) (1 ) 0.05 0−= + − =xf x e x (10)
As 0, (0) 0.95= =x f and when = ∞x it is obtained ( ) 0.05.∞ = −f In addition,
/ ( ) 0, 0.−= − xf x e x x That means there must be a unique solution to the equation (10).
The solution is approximately equal to 4.7439. Hence,
5%
1 2= 4.7439
= =m m m
T
T T T (11)
Then, the PID parameters are determined from Eq. (9).
For the second case, the closed-loop system with input filter is
( )( )1 2
1( )
1 1
=
+ +c m m
G s
T s T s
Without loss of generality, let 1 2m s mT k T= where 0 1.sk< < Doing the same procedure as the
first case, the following equation is obtained
1 2/ /
1 2
1 2
0.05 0
− −
−
− =
−
m mt T t T
m m
m m
T e T e
T T
(12)
Obviously, the expected settling time will be the solution to the equation (12). Let
5% 2= mx T T , then substitute it and 1mT into the equation (12) to have
/
( , ) 0.05 0
1
− −
−
= − =
−
sx k x
s
s
s
k e e
f x k
k
(13)
It can be easily verified that
(0, ) 0.95, ( , ) 0.05= ∞ = −s sf k f k and / ( , ) 0<sf x k .
This implies the equation (13) must have only one solution with respect to .sk The Tab. 1 shows
solutions for some values of sk .
Table 1. Solution with respect to sk .
sk 0.01 0.3 0.5 0.7 0.99
*x 3.005 3.352 3.6761 4.067 4.7202
Nguyen Doan Phuoc, Nguyen Hoai Nam
364
Thus,
5%
2 *=m
T
T
x
and 5%1 * .m s
T
T k
x
= (14)
Hence, given a desired settling time, first we have to choose sk , then time constants 1mT and
2mT are determined by using Eq. 14, and finally the PID parameters are calculated from Eq. 9. It
can be concluded that as 1→sk then
* 4.7439→x . Thus, the second case becomes the first case.
3. SIMULATION RESULTS
3.1. DC Motor Speed Control
Given a speed transfer function of a DC motor as follows [19]
( )
1
=
+
k
G s
Ts
where 100=k , 1.9568=T , the input is the current applied to the DC motor and the output is
the speed of the rotor. The I controller parameter is determined by using the magnitude optimum
method is 2 391.4= =IT kT , whereas the proposed method provides values for the PI
controller as 1.9568= =IT T and 0.02=pk ( 0.1=pk for the second case) which can be
arbitrarily chosen such that the settling time is smaller than any expected value. Simulation
results for two cases are shown in Fig. 2 and Fig. 3. For the first case, the reference is a step
function 0.5=r rad/s, but for the second case, the reference is a sine function
0.5sin(0.3 )=r t . It can be seen that the proposed method gives much better performances in
terms of settling time and overshoot. Especially, when the reference signal is a sine function, the
magnitude optimum method failed to keep the stability and other performances of the close-loop
system.
To compare performances obtained by using different methods, the IAEs (integral absolute
error) at a set point change and IAEd at process disturbance change are used [4]. IAEs is
calculated as follow
0 5 10 15 20 250
0.1
0.2
0.3
0.4
0.5
seconds
ra
d/
s
The speed of rotor
Magnitude optimum method
Proposed method
Setpoint signal
Figure 2. Responses to a step function.
0 5 10 15 20 25 30-0.5
0
0.5
seconds
ra
d/
s
The speed of rotor
Magnitude optimum
Proposed method
Reference signal
Figure 3. Responses to a sine function.
Settling time assignment with PID
365
IAE ( )
f
i
t
s
t
e t dt∫= (15)
where e is the error between the plant output and the set point, it is the initial time, which is the
time of the step change, ft is an appropriate final time. For the IAEd, the formula is the same as
Eq. (15) but the initial time will be the time of a disturbance change.
The IAEs and IAEd for the two cases are shown in Tab. 2. A step disturbance with amplitude
of 0.01 is applied to the input of the plant at the time of 30 seconds to measure the IAEd. In
conclusion, the proposed method provides better performances in term of disturbance
compensation and set point change for both cases.
Table 2. IAE and IAEd for the two cases.
Methods Case 1 Case 2 IAE IAEd IAE IAEd
Magnitude Optimum 2.23 4.265 9.254 15.046
Proposed Method 0.4892 0.9788 0.5405 0.7585
3.2. DC Rotor Position Control
In this section, the PID controllers are designed to control the angular position of the rotor
[19] which has transfer function as below
( )( ) 1= +
k
G s
s Ts
where the values of parameters are the same as in the previous section. The PI controller
parameters are chosen by using the symmetric optimum method as follows
9, 9 17.61= = =Ia T a and
1 0.0017= =pk
kT a
but our method gives
( )
1 2
1 2
1 2
1 2
1 2
0.1575, 2.8,
0.5893
+ +
= = = + + =
+
= =
+ +
m m
p I m m
m m
m m
D
m m
T T T
k T T T T
kT T
T T T
T
T T T
where the desired settling time 5% 2=T seconds and 1 2 0.4216.= =m mT T For the second case,
if 0.5=sk and 5% 2=T , then 2 0.5441=mT and 1 0.2720.=mT Thus, 37.4724=pk ,
=2.7729IT and 0.5759.=DT
For zero overshoot, an input filter ( ) M s is used to cancel the zero of the close-loop
system out. The two proposed PID controllers provide similar performances as shown in Fig. 4.
The results show that the proposed method gives better performances, such as smaller settling
time as desired. When parameters 1mT and 2mT increase the settling time is larger but the
Nguyen Doan Phuoc, Nguyen Hoai Nam
366
overshoot is smaller and vice versa. The IAEs and IAEd are shown in Tab. 3, where a step
disturbance with amplitude of 0.01 is applied to the input of the plant at the time of 70 seconds
to measure the IAEd. Again, the IAEs and IAEd proved that the proposed method gives much
better performance than the symmetric optimum method.
0 10 20 30 40 50 60 700
1
2
3
4
5
6
seconds
ra
d
The angular position of rotor
Symmetric optimum
Proposed method
Setpoint signal
Figure 4. Responses to a step function.
Table 3. IAEs and IAEd.
Methods Performance
IAEs IAEd
Symmetric Optimum 24.66 104.94
Proposed Method 0.7858 0.2080
4. CONCLUSIONS AND FUTURE WORKS
In the work, a PID design method with desired settling time is proposed for plants with
transfer functions of the following types: first-order system, second-order system and integral
first-order system. The method provides better performances than both the magnitude optimum
method and the symmetric optimum method in the sense of overshoot and settling time. In
addition, the settling time can be achieved as expected for a larger frequency band. Some
examples for DC motor control are used to illustrate the method.
Future works will focus on robustness, effects of model error and other types of transfer
functions, may be higher-order system with time delay.
REFERENCES
1. Ziegler J. G. and Nichols N. B. - Optimum settings for automatic controllers, Trans.
ASME 64 (1942) 759-765.
2. Oldenbourg R. C. and Hans S. - A uniform approach to the optimum adjustment of control
loops, Trans. ASME 18 (1954) 1-9.
3. Hang C. C., Astrom K. I. and Ho W. K. - Refinements of the Ziegler-Nichols tuning
formula, IEE Proc. D. 138 (1991) 111-118.
Settling time assignment with PID
367
4. Finn H. and Bernt L. - Relaxed Ziegler-Nichols Closed-Loop Tuning of PI Controllers,
Modeling, Identification and Control 34 (2013) 83-97.
5. Kessler C. - Das symmetrische optimum, Regelungstechnik 11 (1958) 395-400.
6. Jeffrey W. U. and Mohammed S. - Magnitude and symmetric optimum criterion for the
design of linear control systems: What is it and how does it compare with the others?
IEEE Trans. Indus. Appl. 26 (1990) 489-497.
7. Konstantinos G. P. and Nikolaos I. M. - Extending the symmetrical optimum criterion to
the design of PID type-p control loops, Journal of Process Control 22 (2012) 11-25.
8. Konstantinos G. P., Nikolaos D. T. and Nikolaos I. M. - Revisiting the magnitude
optimum criterion for robust tuning of PID type-I control loops, Journal of Process
Control 22 (2012) 1063-1078.
9. Astrom K. J. and Hagglund T. - Automatic tuning of simple regulators with specifications
on phase and amplitude margins, Automatica 20 (1984) 645-651.
10. Weng K. H., Chang C. H. and Lisheng S. C. - Tuning of PID controllers based on gain
and phase margin specifications, Automatica 31 (1995) 497-502.
11. Sheng Y. C. and Ching C. T. - Tuning of PID controllers based on gain and phase margin
specification using fuzzy neural network, Fuzzy Sets and Systems 101 (1999) 21-30.
12. Cominos P. and Munro N. - PID controllers: recent tuning methods and design to
specification, IEE Proc. Control Theory Appl. 149 (2002) 46-53.
13. Keyu L. - PID tuning for optimal closed-loop performance with specified gain and phase
margins, IEEE Trans. Control Syst. Tech. 21 (2013) 1024-1030.
14. Tyreus B. D. and Luyben W. L. - Tuning PI controllers for integrator/dead time processes,
Ind. Eng. Chem. 31 (1992) 2626-2628.
15. Astrom K. J. and Hagglund T. - PID controllers: Theory, design and tuning, ISA 1995.
16. Astrom, K. J. and Hagglund T. - Advanced PID control, The Instrumentalion, Systems,
and Automation Society, 2005.
17. Phuoc, N. D. - Linear Control Theory, Science and Technology Publisher, 2007.
18. Ntogramatzidis L. and Ferrante A. - Exact tuning of PID controllers in control feedback
design, IET Cont. Theo. Appl. 5 (2011) 565-578.
19. Siri W. and El S. - Identification and control of a DC motor using back-propagation neural
networks. IEEE Transactions on Energy Conversion 6 (1991) 663-669.
Các file đính kèm theo tài liệu này:
- 8808_37438_1_pb_8116_2061747.pdf