Chapter II
Chapter II
Where: K '=K ϕ ϕ .
Where:
V a is the applied armature voltage.
Φ is constant [14].
Where:
T is the torque developed by the motor.
J is the equivalent moment of inertia of motor shaft and load referred to the motor.
B is the equivalent coefficient of motor and load referred to the motor shaft.
T L The resistant torque.
Energy storage is associated with the magnetic fields produced by i f and i a and with the
kinetic energy of the rotating parts. The above equations are a set of nonlinear (because of
products i f (t)Ω(t) and i f (t)i a (t )) state equations with state variables i f , i a and Ω . The
solution must be obtained numerically.
II-2-3) Transfer Functions and Block Diagrams:
In the simple linear case of motor response to changes in armature voltage, it is assumed that
the field voltage is constant and steady- state is existing on the field circuit, i.e. I f = constant.
Equations (II−1), (II−2) and (II−3) now become linear as given below:
d
V a (t )=K ' Ω(t)+ R a i a (t )+ L a i (t). (II 4)
dt a
d
T (t)=K ' i a (t )=J Ω( t)+BΩ (t)+T L (t) . (II 5)
dt
Laplace transforming eqs (II 4) and (II 5):
V a (s)=K ' Ω(s)+(Ra + s La )I a (s). (II 6)
K
’
TL
Ω
Ω2 Ω1 Ωn
{
T =K ϕ ϕ I a .
V −R a I a (II
Ω= a .
Kϕ ϕ
12)
Va
Pour T =0 ; I a=0 Ω= : When the value of ϕ is reduced, the speed increases.
Kϕϕ
Va Va
Pour Ω=0; I a= T =K ϕ ϕ : When the value of ϕ is reduced, the torque decreases.
Ra Ra
This setting is bad from a technical point of view; the characteristics being concurrent. It is
good from an economic point of view because the power dissipated in the inductor is very
low compared to the power absorbed; The engine efficiency will not be affected. However,
the following remarks should be made :
With this process, the motor speed can only be increased compared to its nominal speed. If
the load torque is constant T =K ϕ ϕ I a=cst , the current will increase when the flow is reduced
and the motor may heat up.
Under these conditions, the engine must be sized accordingly. It should be noted that this
disadvantage does not occur if the drive is at constant power because by virtue of the
relationship P=V a I a=cst , as the voltage U a is constant, the current I a will remain constant.
T
ϕn
ϕ1
ϕ2
ϕn > ϕ1 > ϕ 2
TL
Ω
Ωn Ω1 Ω2
{
T =K ' I a
V a−R a I a (II
Ω= .
K'
13)
Va ∆Va
Pour T =0 ; I a=0 Ω= ∆ Ω= .
K' K'
∆T
When the voltage value is reduced, the torque decreases. Let's calculate the slope:
∆Ω
2
∆T K'
= =cst .
∆ Ω Ra
U2 U1 Un U1 > U2
TL
Ω
Ω2 Ω1 Ωn
Circuit de
puissance
Ia
Source Va M
Rsh Tachymeter
Imes Ωmes
Command Regulator
(I)
-+ Regulator
(Ω)
-+ Ωref
exterior
command
Duty cycle circuit
ID
IT K Ia Ie
VT
DC
source
V DL Va M Re Ve
V a =V → I T =I a →V T =0 → V D =−V → I D=0.
V a =0 → I T =0 →V T =V →V D=0 → I D =I a .
Figure (II 7): Motor currents and voltages, IGBT and diode.
According to the curves shown in Figure (II 7), the average value of the voltage across the
motor is:
t on
1 t
V a = ∫ V dt= on V =αV . (II 14)
T 0 T
where: α ∈ [0,1] α: is called a duty cycle.
The rotational velocity expression is related proportional to α [1]:
V a−Ra I a α V −Ra I a
Ω= = . (II 15)
K' K'
Speed
sensing
K
’
Figure (II 10): Open-loop block diagram of separately excited DC motor drive.
The response due to a step change in the reference voltage is obtained by setting T L to zero.
From Figure (II 10), we can obtain the speed response due to reference voltage as:
Ω(s) K 2 K ' /(R a B)
=GΩ−V = . (II
V r (s) 2 2
(τ a τ m)s +(τ a +τ m)s +1+ K ' /(R a B)
17)
The response due to a change in load torque T L can be obtained by setting V r to zero. The
block diagram for a step change in load torque disturbance is shown in Figure (II 11):
T L ( s) 1 Ω(s)
B(s
T (s)
τ m+ 1¿
I a( s) 1 E g( s)
K -1 K
R a (s
’ ’
τ a+ 1¿
Figure (II 11): Open-loop block diagram for torque disturbance input.
Vc
K r =α = . (II 20)
V cm
Therefore, a converter can be modelled with a transfer function Gc (s) of a certain gain and
phase delay as described by:
τr
Gc ( s )=K r e .
Kr
G r ( s )= . (II 21)
1+ s τ r
Once a switching device is turned on, its gating signal cannot be changed. There is a delay
until the next device is gated and a corrective action is implemented. The delay time is
generally half of the interval between two switching devices.
T L ( s)
V r( s) Ve
K2
Va 1 Ia
K’
T 1 Ω(
R a (s B (s
Vb E g( s) τ a+ 1¿ τ m+ 1¿
K’
K1
Figure (II 12): Block diagram for closed-loop control of separately excited dc motor.
The closed-loop step response due to a change in reference voltage can be found from Figure
(II 12) with T L =0. The transfer function becomes:
K2K
Ω(s) Ra B
= 2
. (II
V r (s) K +K 1 K2 K
(τ a τ m)s 2 +(τ a +τ m)s +1+
Ra B
22)
The response due to a change in the load torque T L can also be obtained from Figure (II 12)
by setting V r to zero. The transfer function becomes:
Ω(s) 1 s τ a+1
= .
T (s) B K +K1 K 2 K
2
(II
(τ a τ m )s 2+(τ a+ τ m )s+1+
Ra B
23)
Ω(s) Ω(s ) I a ( s)
= × .
V a (s) I a (s) V a ( s)
T1
ia* Vc va 1 ia T 1 Ω
R a + sL a
Kb
B + sJ
Current controlled ve
iam controller converter
Kb
Hc
V a(s)
1 Ia(s) Kb Ω(s
R a + sL a B + B l + sJ
E (s)
K b2
B + B l + sJ
The blocks in Figure (II 14) can be simplified to obtain the following transfer function
relationships:
Ω(s) Kb
= .
I a (s) B (1+s τ m )
I a (s ) 1+s τ m
=K m .
V a (s) (1+ s τ 1)(1+ s τ 2)
Where:
B
Km= 2 .
K + Ra B
b
[( ) √( ) ( )]
2 2
−1 1 1 Ra B Ra B Ra B+ K b
s1 ,2 = ,− = − + ± + −4 . (II 24)
τ1 τ2 2 La J La J J La
Two block representation of the open-loop motor transfer function between the output speed
and the input voltage is shown in Figure (II 15), where total viscous friction Bt =B+ B L.
1+ s m I s( s) K b /B t
K1
(1 + s 1)( 1+ s 2) 1+ s m
The overall closed-loop system with a current feedback loop is shown in Figure (II 16),
where Bt is the total viscous friction of the motor and the load. The converter can be
represented by Eq (II 21). Using proportional-integral type of control, the transfer functions
of the current controller Gc (s) and the speed controller Gs (s) can be represented as:
K c (1+s τ c )
Gc (s)= . (II 25)
s τc
K s (1+s τ s )
Gs (s)= . (II 26)
s τc
Where:
K c and K s are the gains of the current and speed controllers, respectively.
τ c and τ s are the time constants of the current and speed controllers, respectively.
Ω
¿
Ωr G s( s)
ia*
G c( s)
vc
G r( s )
va
K1
1+ s m ia K b /B t
(1 + s 1)( 1+ s 2) 1+ s m
Speed Current Converter
Limiter controller
controller
Ωmr Hc
GΩ (s)
A dc tachogenerator along with a low-pass fitter of a time constant less than 10 ms is often
required. The transfer function of the speed feedback filter can be expressed as:
KΩ
GΩ (s)= .
1+ s τ Ω
Where:
K Ω and τ Ω are the gain and time constants of the speed feedback loop.
Figure (II 17) shows the block diagrams with the current-control loop. i am is the feedback
armature- motor current. The simplified diagram is shown in Figure (II 18).
K c (1 + s c) Vc Kr Va 1+ s m
i a* K1 ia
s c 1+ s r (1 + s 1) (1 + s 2)
iam
Hc
Kc K 1 m
i a* Kr ia
Tc 1+ s 3
Hc
G(s) H (s )=
( Km Kc K r H c
τc ) (1+s τ c )(1+ s τ m)
s (1+ s τ 1 )(1+s τ 2 )(1+ s τ r )
.
To reduce the system to a second order, the following assumptions can be made for practical
motor drives:
K
G(s) H (s )= . (II
(1+ s τ 1 )(1+s τ r )
27)
Km Kc K r H cτ m
Where: K= .
τc
The characteristic equation of the loop gain in Eq. (14.105) is given by:
This gives the natural frequency ω n and damping factor ξ of the current-control loop as:
ω n=
√ 1+ K
τ1 τr
.
τ 1 +τ r
τ1τ r .
ξ=
2 ωn
Setting the damping ratio ξ=0.707 for critically damped, assuming K>>1 and τ 1 > τ r , the gain
of the current controller can be expressed as:
1 τ1 τc 1
K c= .
2 τr Km K r Hc τ m
I a (s) Ki
¿ = . (II
I a (s) (1+ s τ i)
28)
Where:
τ3
τi= .
1+ K 1
τ 3 =τ 1 +τ r.
K i=
1
(K1
H c 1+ K 1
. )
Km Kc Kr H c τm
K 1= .
τc
Replacing the speed controller with its transfer function in Eq (II 26) and approximating
current-control loop in Eq (II 28), the block diagram of the outer speed- control loop is
shown in Figure (II 19). The loop gain function is expressed as:
G(s) H (s )=
( K s Ki Kb HΩ
B τs ) (1+s τ s )
s (1+ s τ i)(1+ s τ m)(1+ s τ Ω )
.
Ωr
¿ K s (1 + s s) ia* Ki ia K b /B t Ω
s s 1+ s i 1+ s m
Ωmr HΩ
1+ sΩ
Assuming 1+τ m ≈ τ m and combining the delay time τ s of the speed controller with the delay
time τ Ω of the speed feedback filter to an equivalent delay time τ e=τ s + τ Ω:
G(s) H (s )=K Ω ( )
K s (1+ s τ s )
τ s s2 (1+ s τ e )
. (II 29)
Where:
τ e=τ s + τ Ω.
K i Kb HΩ
KΩ= .
B τm
The closed-loop transfer function of the speed in response to speed reference signal is given
by:
[ ]
KΩKs
(1+ s τ s)
Ω(s) 1 τs
¿ = . (II
Ω (s) K Ω 3 2 KΩ Ks
s τe + s + s KΩ Ks+
τs
30)
Equation (II 30) can be optimized by making the magnitude of its denominator such that the
coefficients of ω2 and ω4 in the frequency domain are zero. This should widen the bandwidth
operating over a wide frequency range. It can be shown that the gain K s and the time constant
τ s of the speed controller under the optimum conditions are given by:
1
Ks= .
2 KΩ τe
τ s=4 τ e .
Substituting K s and τ s into Eq. (II 30) gives the optimized closed-loop transfer function of
the speed in response to speed input signal as:
Ω(s)
¿ =
1
[
(1+ 4 s τ e )
Ω (s) K Ω 1+ 4 s τ e +8 s 2 τ 2e + 8 s 3 τ 3e
.
]
It can be shown that the corner frequencies for the open-loop gain function H(s)G(s) are
1/(4 τ e ¿ )¿ and 1/τ e. The gain crossover occurs at 1/(2 τ e ¿ )¿ at a slope of the magnitude
response—20 dB/decade. As a result, the transient response should exhibit the most desirable
characteristic for good dynamic behaviour. The transient response is given by:
[ ( )]
−t −t
1 + e 2 τ −2 e 4 τ cos √
1 3t
Ω(t)= e e
. (II 31)
KΩ 4 τe
This gives a rise time of 3.1τ e, a maximum overshoot of 43.4%, and a settling time of 16.5 τ e.
The high overshoot can be reduced by adding a compensating network with a pole-zero in the
speed feedback path, as shown in Figure (II 19). The compensated closed-loop transfer
function of the speed in response to speed input signal is:
Ω(s)
¿ =
1
[ 1
Ω (s) K Ω 1+ 4 s τ e +8 s 2 τ 2e + 8 s 3 τ 3e
.
]
¿
Ωr 1 1 1+ 4 e s
1+ 4 e s HΩ 1+ 4 e s +8 e 2 s2 +8 e 3s3
[ ( )]
−t −t
1 + e 4 τ − e 4 τ sin √
1 2 3t
Ω(t)= e
.e
(II 32)
KΩ √3 4τe
This gives a rise time of 7.6 τ e, a maximum overshoot of 8.1%, and a settling time of 13.3 τ e.
It should be noted that the overshoot has been reduced to approximately 20% of the value,
and the settling time has come down by 19%. But the rise time has increased, and the
designer must make a compromise on the rise time and the overshoot.
Where:
e (t): is the error between the desired and actual speed.
K p : is the proportional gain.
Integral control reduces the time to converge to setpoint and eliminates some of the steady
state error. The higher the integrated charge gain, the shorter the time to converge to setpoint,
but since the integral control cumulates whole previous errors, it may result in an overshoot
of the return value.
Derivative control:
The role of derivative control is to forecast future values based on the control error. The
expression is as follows:
de (t)
U (t )=K d .
dt
Differential control improves tuning time and system stability [16].
(
Gc ( s )=K p 1+
1
Tis )
+s T d . (II 34)
Where:
Kp
T i= : is the integral time.
Ki
Kd
T d= : is the derivative time.
Kp
There are two methods called Ziegler–Nichols tuning rules: the first method and the second
method. We shall give a brief presentation of these two methods.
First Method:
In the first method, we obtain experimentally the response of the plant to a unit-step input, as
shown in Figure (II 21). If the plant involves neither integrator(s) nor dominant complex-
conjugate poles, then such a unit-step response curve may look S-shaped, as shown in Figure
(II 22). This method applies if the response to a step input exhibits an S-shaped curve. Such
step-response curves may be generated experimentally or from a dynamic simulation of the
plant. The S-shaped curve may be characterized by two constants, delay time L and time
constant T. The delay time and time constant are determined by drawing a tangent line at the
inflection point of the S-shaped curve and determining the intersections of the tangent line
with the time axis and line c(t)=K, as shown in Figure (II 22).
C ( s ) K e− Ls
= . (II
U ( s ) T s +1
35)
Ziegler and Nichols suggested to set the values of K p , T i, and T d according to the formula
shown in Table (II1):
Type of Controller Kp Ti Td
T
P ∞ 0
L
T L
PI 0.9 0
L 0.3
T
PID 1.2 2L 0.5 L
L
Table (II1): Ziegler–Nichols Tuning Rule Based on Step Response (First Method).
Notice that the PID controller tuned by the first method of Ziegler–Nichols rules give:
2
1
(
Gc ( s )=K p 1+
1
Tis ) T
(
+T d s =1.2 1+
L
1
2 Ls
+ 0.5 Ls . )
Gc ( s )=0.6 T
(s+ )
s
L
.
(II 36)
Thus, the PID controller has a pole at the origin and double zeros at −1/ L.
Second method:
In the second method, we first set T i=¿ and T d=0. Using the proportional control action only
(see Figure (II 23)), increase K p from 0 to a critical value K cr at which the output first
exhibits sustained oscillations. (If the output does not exhibit sustained oscillations for
whatever value K p may take, then this method does not apply.) Thus, the critical gain K cr and
the corresponding period Pcr are experimentally determined (see Figure (II 24)). Ziegler and
Nichols suggested that we set the values of the parameters K p , T i and T d according to the
formula shown in Table (II2):
Type of Controller Kp Ti Td
P 0.5 K cr ∞ 0
1
PI 0.45 K cr P 0
1.2 cr
PID 0.6 K cr 0.5 Pcr 0.125 Pcr
Table (II2): Ziegler–Nichols Tuning Rule Based on Critical Gain K cr and Critical Period Pcr
(second Method)
Figure (II 24): Sustained oscillation with period Pcr ( Pcr is measured in sec.)
Notice that the PID controller tuned by the second method of Ziegler–Nichols rules give:
2
4
(
Gc ( s )=K p 1+
1
Tis ) (
+T d s ¿ 0.6 K cr 1+
1
0.5 Pcr s
+0.125 Pcr s )
Gc ( s )=0.075 K cr P cr
(s +
P cr
s
)
.
(II 37)
Thus, the PID controller has a pole at the origin and double zeros at s=−4 /P cr .