Regulation Automatique PDF
Regulation Automatique PDF
Regulation Automatique PDF
lments de rgulation
automatique
La rgulation automatique est une technique de l'ingnierie orant les mthodes et
les outils ncessaires la prise de contrle d'une ou plusieurs grandeurs physiques
d'un systme en vue d'en imposer le comportement. Avec le qualicatif
automa-
tique, on admet qu'aucune intervention humaine n'est ncessaire pour atteindre cet
objectif.
Ce chapitre a pour but de sensibiliser l'tudiant ce qu'est la rgulation automatique
et lui donner une information susante pour aborder des problmes simples.
G(s)
H(s)
aboutissant au sommateur
Celui-ci
peut faire la somme ou la dirence des deux signaux ; dans le premier cas, on parlera
de raction positive et, dans le deuxime, de raction ngative ou plus simplement
de contre-raction.
L'criture des quations correspondantes et leur rsolution nous permettent de remplacer l'ensemble du schma fonctionnel par un seul bloc : la fonction de transfert
en boucle ferme
Gf (s).
91
3.
contre-raction
Y (s) =
G(s)
W (s)
1 + G(s) H(s)
(3.1)
Gf (s)
G(s)
Y (s)
=
W (s)
1 + G(s) H(s)
(3.2)
Cette fonction de transfert est susamment importante pour qu'on lui donne le
nom de
formule de l'automaticien .
G(s) H(s)
forme la
Go (s)
F (s)
= G(s) H(s)
E(s)
(3.3)
H(s)
est un
Go (s) =
Y (s)
= G(s),
E(s)
Gf (s) =
Go (s)
1 + Go (s)
(3.4)
92
2008 freddy.mudry@gmail.com
Ecart
Consigne
Commande
Grandeur rgle
Perturbation
v(t)
e(t)
w(t)
u(t)
x(t)
Rgulateur
Comparateur
Actionneur
Processus
y(t)
Grandeur rgle
mesure
Capteur
Systme rgler
lment
Fonction
Exemple
Comparateur
Amplicateur
direntiel
Rgulateur
le signal de commande
Organe de
commande
Processus
L'installation asservir
Il forme une image
Capteur
y(t)
Amplicateur
de puissance
Moteur
aussi dle
Signal
Rgulateur PID
x(t)
Remarques
Dynamo
tachymtrique
Exemple
w(t)
application donne
Grandeur
rgle
x(t)
Grandeur
rgle
y(t)
mesure
Commande
u(t)
v(t)
e(t)
[rad/sec]
L'image de la
vitesse [V]
Tension de
l'actuateur
commande [V]
cart
rotation
Perturbation
La vitesse de
y(t)
Couple extrieur
[Nm]
Tension d'cart
[V]
2008 freddy.mudry@gmail.com
93
3.
exemple de composants et signaux est donn dans le cas particulier d'un asservissement de la vitesse de rotation d'un moteur lectrique. On notera qu'avec le schma
adopt, le
systme rgler
u(t)
y(t).
E(s)
U(s)
Gc(s)
Ga1(s)
Ga2(s)
Y(s)
Ga1 (s)
Ga2 (s)
Gc (s) ;
v(t) ;
v(t).
systme rgler
Y (s)
= Ga1 (s) Ga2 (s)
Ga (s)
U (s) v(t)=0
2. Le
(3.5)
(3.6)
Dans ce cas,
3.
Y (s)
Go (s)
Gw (s)
=
W (s) v(t)=0 1 + Go (s)
dont le but est de suivre la consigne
94
w(t).
2008 freddy.mudry@gmail.com
(3.7)
rgulation de maintien
Ga2 (s)
Y (s)
Gv (s)
=
V (s) w(t)=0 1 + Go (s)
y(t)
(3.8)
w(t)
gale la consigne
v(t).
1.5
w(t)
y(t)
1
0.5
0
0
10
12
14
16
18
20
1.5
w(t)
v(t)
y(t)
1
0.5
0
0
10
temps
12
14
16
18
20
La gure 3.4 illustre les rponses temporelles obtenues dans les deux modes de rgulation automatique. Dans la ralit, les deux modes coexistent obligatoirement car le
rgulateur ragit toute forme d'erreur quelle qu'en soit la cause (consigne variable
(3.9)
(3.10)
Ka
2008 freddy.mudry@gmail.com
Ga (s)
(gure 3.5).
95
3.
W(s)
E(s)
U(s)
Ka
Y(s)
Ga(s)
Ga (s)
1
Y (s)
=
U (s)
1 + s
Y (s)
1
= Ka
E(s)
1 + s
Go (s)
Elle est caractrise par son gain
Ka
et sa constante de temps
Gf (s)
En crivant
Gf (s)
1
Ka 1+s
Ka
Y (s)
=
=
1
W (s)
1 + Ka + s
1 + Ka 1+s
Gf (s) =
Ka
1
1 + Ka 1 + s 1+K
a
Kf =
Ka
< 1,
1 + Ka
Comportement transitoire
Ka
et une
f =
<
1 + Ka
tmes est dcrit par les ples de la fonction de transfert. Dans ce cas,
Gf (s)
tant
p1 =
1 + Ka
Le ple tant partie relle ngative, le systme asservi sera toujours stable suivant
la valeur du gain, on obtient les rponses indicielles illustres par les graphes de la
gure 3.6.
96
2008 freddy.mudry@gmail.com
0.8
0.8
0.6
0.6
0.4
0.4
0.2
0.2
y(t)
b.o.
0
0
Ka = 4
0
3
0.8
0.8
0.6
0.6
0.4
0.4
y(t)
0.2
0.2
Ka = 8
0
0
Ka = 12
0
3
temps
temps
Im
Re
0
-1/
Fig. 3.6.: Rponses indicielles d'un systme asservi d'ordre 1 et dplacement du ple
2008 freddy.mudry@gmail.com
97
3.
p1
Ka
aug-
Ka
Ga (s)
1
Y (s)
=
U (s)
s (1 + s )
Y (s)
1
= Ka
E(s)
s (1 + s )
Go (s)
ferme donne
Gf (s)
En crivant
Gf (s)
Y (s)
W (s)
1
s (1+s )
1
Ka s (1+s
)
Ka
1+
s2
Ka
+ s + Ka
Gf (s) =
1
1+s
1
Ka
s2 Ka
=
1 + 2
s
n
s
n
2
Gf (s) caractrise
par :
1. Un gain DC gal 1 grce au terme intgrateur
Go (s)
Kf = 1
2. Une pulsation naturelle variable qui augmente avec le gain
r
n =
Ka
98
n
1
=
2 Ka
2 Ka
2008 freddy.mudry@gmail.com
Ka
Ka
1.5
10
0.5
y(t)
15
b.o.
0
0
10
Ka = 0.2
0
15
1.5
0.5
0.5
10
15
y(t)
1.5
Ka = 1
0
0
10
Ka = 4
0
15
temps
10
15
temps
Im
p1
p2
p1
Re
0
-1/
p
2
Fig. 3.7.: Rponses indicielles d'un systme asservi d'ordre 2 et dplacement des
2008 freddy.mudry@gmail.com
99
3.
Comportement transitoire
< 1)
Ka .
Puis,
est constant
1
n =
2 Ka
Ka
1
=
On peut dcrire plus prcisment le comportement transitoire du systme en calculant ses ples qui sont les racines du dnominateur
D(s) = s2 + 2n s + n2
Ces ples sont au nombre de 2 et valent
p1,2
q
p
= n (n )2 n2 = n n 2 1
faibles gains , on a
1
Ka
1
4
p1 = n + n
une lente :
p2 = n n
une rapide :
Dans le cas de
p
2 1 < 0
2 1 < p1 < 0
forts gains , on a
<1
Ka >
1
4
p1,2 = n j n
p
1 2
bf =
1
1
=
= 2
|Re (p1,2 )|
n
p = |Im (p1,2 )| = n
1 2
Les deux ples tant partie relle ngative, on en conclut que le systme asservi
sera toujours stable quelle que soit la valeur de
100
Ka
suprieure 0.
2008 freddy.mudry@gmail.com
Ka
1,
p1,2
en
c'est--dire
Ka 1/(4 ),
1/
1
< p1,2 < 0
Leur dplacement respectif conduit une rponse indicielle qui tend de plus
en plus rapidement vers la valeur asymptotique.
2. Pour
< 1,
c'est--dire
Ka > 1/(4 ),
Ka
Re (p1,2 ) = n =
Im (p1,2 ) = n
p
1
1 2 =
1
2
r
Ka
1
4
Ka
Ga (s) =
1
(1 + s 1 ) (1 + s 2 ) (1 + s 3 )
Go (s)
1
Y (s)
= Ka
E(s)
(1 + s 1 ) (1 + s 2 ) (1 + s 3 )
Ka
1,2,3 .
ferme donne
2008 freddy.mudry@gmail.com
101
3.
On voit ainsi que le systme d'ordre 3 boucl sur lui-mme conduit une fonction
de transfert
Gf (s)
1. Un gain DC
Kf =
Ka
1 + Ka
Ka
Ka .
p1,2,3
en
permet de voir que les ples partent des ples en boucle ouverte
/3.
Une illustration
en est prsente dans la gure 3.8. On constate ainsi, que dans le cas de systmes
d'ordre suprieur 2, la contre-raction peut conduire des ples partie relle
positive. Le systme devient alors instable. Une illustration des rponses indicielles
est donne la gure 3.8.
3.2.4. Conclusions
Des trois exemples ci-dessus, on peut tirer les conclusions gnrales suivantes. Dans
un systme asservi :
1. L'ordre du systme asservi reste le mme que celui du systme en boucle
ouverte.
2. Le gain en boucle ferme
Kf
Ka ;
il tend vers 1 si
Ka
102
Ka
Kf
Ka .
2008 freddy.mudry@gmail.com
vaut 1.
Im
p1
p3
Re
-1/3
-1/2
-1/1
p2
1.5
2
b.o.
Ka = 1
1.5
y(t)
1
1
0.5
0.5
0
0
10
10
4
6
temps
10
4
Ka = 4
Ka = 12
1.5
y(t)
2
1
0
1
0.5
2
0
0
4
6
temps
10
Fig. 3.8.: Dplacement des trois ples b.f. dans le plan complexe et rponses indi-
2008 freddy.mudry@gmail.com
103
3.
N = 14 et un rendement = 75%
15 [V], on demande :
et que le
1. Dessinez les schmas technologique et fonctionnel du systme asservi comprenant l'amplicateur, le moteur (modle d'ordre 1), le rducteur et le capteur
de position.
2. Calculez la valeur numriques des paramtres du moteur
le gain
Kmot , mot
ainsi que
du capteur de position.
Ka
optimum.
5. Calculez le temps de rglage de la position du disque.
Solution
1) Schmas technologique et fonctionnel du systme asservi
104
2008 freddy.mudry@gmail.com
2008 freddy.mudry@gmail.com
105
3.
106
2008 freddy.mudry@gmail.com
y(t)
1.2
1.2
0.8
0.8
0.6
0.6
0.4
0.4
0.2
0.2
Ka = 15
y(t)
0.05
0.1
0.15
0.2
0.8
0.8
0.6
0.6
0.4
0.4
0.2
0.2
Ka = 45
0
0
0.05
0.1
temps [sec]
0.15
Ka = 30
0
0.05
0.1
0.2
Ka = 60
0
0.2
0.15
0.05
0.1
temps [sec]
0.15
0.2
2008 freddy.mudry@gmail.com
107
3.
Gc (s) =
Gc (s)
2+s
5s
et du processus
Ga (s) =
Ka
Ga (s)
1
1 + 2s + 3s2
indicielle ait un dpassement d'environ 10% (ceci est un bon compromis entre le
temps de monte et la dure du rgime transitoire).
Initialisation du chier
% fichier : xple_asserv.m
% fmy - janvier 2006
% Analyse d'un asservissement
clear all ; close all ; clc ;
format compact ; format short g ;
Donne des fonctions de transfert
% fonctions de transfert
num = [1,2] ; den = [5,0] ;
Gcs = tf(num,den) ;
num = 1 ; den = [3,2,1] ;
Gas = tf(num,den) ;
% fonction de transfert en b.o.
Ka = 1 ;
Gos = Ka *Gcs * Gas
Dans la fentre de commandes, Matlab ache alors le rsultat suivant :
Transfer function :
s + 2
--------------------15 s^3 + 10 s^2 + 5 s
Dynamique du systme en b.o.
Go (s)
systme :
108
2008 freddy.mudry@gmail.com
pk = roots(den_Go)
tau_k = 1 ./ abs(real(pk))
T_k = 2*pi ./ abs(imag(pk))
t_trans = 5*tau_k
N_osc = t_trans ./ Tp_k
On obtient alors les rsultats suivants :
zk =
-2
tau_k =
Inf
3
3
t_trans =
Inf
15
15
La constante de temps
pk =
-0.333 j0.471
T_k =
Inf
13.329
13.329
N_osc =
NaN
1.1254
1.1254
k
et la priode
Tk
Avec
Ka = 1,
le calcul de la fonction
Transfer function :
15 s^4 + 40 s^3 + 25 s^2 + 10 s
-----------------------------------------------------225 s^6 + 300 s^5 + 265 s^4 + 140 s^3 + 50 s^2 + 10 s
Ce rsultat est surprenant, car on sait que la boucle d'asservissement ne change
pas l'ordre du systme qui dans notre cas vaut 3. Le rsultat ci-dessus vient du fait
qu'une simplication est possible entre le numrateur et le dnominateur. On peut
le voir en demandant l'achage de
Gf (s)
avec la fonction
zpk
zpk(Gfs)
qui ache le rsultat suivant
Zero/pole/gain :
0.0667 s (s+2) (s^2 + 0.6667s + 0.3333)
--------------------------------------------------------------s (s+0.4431) (s^2 + 0.2236s + 0.3009) (s^2 + 0.6667s + 0.3333)
2008 freddy.mudry@gmail.com
109
3.
Gfs = minreal(Gfs)
On obtient alors le rsultat attendu :
Transfer function :
0.06667 s + 0.1333
----------------------------------s^3 + 0.6667 s^2 + 0.4 s + 0.1333
z1 = 2 [1/sec]
p1 = 0.443 [1/sec]
p2,3 = 0.1118 j 0.537 [1/sec]
On en dduit les temps caractristiques suivants
1 =
2,3 =
ttrans =
Tp,2,3 =
Nosc =
Du thorme des valeurs limites, on tire les valeurs initiale et nale de sa rponse
indicielle
y(t 0) = s Y (s)|s
0.1333
1
=
=0
= s Gf (s)
s
s3 s
s
y(t ) = s Y (s)|s0
1
0.1333
= s Gf (s)
=
=1
s
0.1333
s0
Gc (s), la valeur
y() est gale la consigne. Par contre, cette rponse indicielle est loin d'tre
optimale puisque l'on peut compter presque quatre priodes d'oscillation pendant
la dure transitoire ; ce qui est conrm par le deuxime graphe de la gure 3.9. On
est donc amen devoir rduire le gain
Ka
oscillation et, ventuellement, rduire la dure de rglage ou, ce qui est quivalent,
son temps d'tablissement.
110
2008 freddy.mudry@gmail.com
w(t)
% reponse temporelle
npts = 1000 ; tmax = 30 ;
dt = tmax/npts ;
tt = 0 :dt :tmax-dt ;
% saut de consigne d'amplitude A
A = 2;
wt = A*ones(size(tt)) ;
wt(1) = 0 ;
On calcule ensuite la rponse temporelle un signal quelconque avec la fonction
lsim
indicielles
figure ;
Ka = [0.2,0.3,0.5,0.65] ;
for k1 = 1 :length(Ka)
Gfs = Ka(k1)*Gos / (1+Ka(k1)*Gos) ;
Gfs = minreal(Gfs) ;
yt = lsim(Gfs,wt,tt) ;
subplot(2,2,k1) ;
plot(tt,yt,'LineWidth',2) ;
axis([-0.5,tmax,-0.05,2.5]) ; grid on ;
texte = ['K_a = ', num2str(Ka(k1),2)] ;
if (k1 == 1) | (k1 == 3), ylabel('y(t)'), end ;
if (k1 == 3) | (k1 == 4), xlabel('temps'), end ;
legend(texte,4) ;
end ;
2008 freddy.mudry@gmail.com
111
3.
1.5
0.5
0.4
y(t)
dy(t) / dt
1
0.3
0.2
0.5
0.1
0
b.o.
y(t)
10
20
30
temps
40
50
1.2
1.2
0.8
0.8
0.6
0.6
0.4
0.4
0.2
10
20
30
temps
40
50
0.2
Ka = 0.2
0
0
y(t)
b.f.: Ka = 1
10
20
30
40
Ka = 0.3
0
50
1.2
1.2
0.8
0.8
0.6
0.6
0.4
0.4
0.2
10
20
30
40
50
0.2
Ka = 0.5
0
0
10
20
30
temps
40
Ka = 0.65
0
50
10
20
30
temps
40
50
Fig. 3.9.: Rponses indicielles d'un systme asservi avec son optimisation
112
2008 freddy.mudry@gmail.com
u2(t)
u3(t)
Dtecteur
de phase
u2(t)
u4(t)
Filtre
VCO
passe-bas
u1(t)
10
u1(t)
5
u (t)
2
0
u3(t)
5
0
3
2.5
u (t)
4
2
1.5
1
0.5
0
3
temps [ms]
Les lments ncessaires la ralisation d'un tel systme sont prsents avec les
signaux dans la gure 3.10. On y trouve un dtecteur de phase, un ltre passe-bas
et un oscillateur command en tension (VCO). Avec le circuit PLL, on asservit la
frquence de l'oscillateur interne de manire le synchroniser au signal de rfrence
appliqu en entre. Ce circuit est donc un systme contre-raction dont le signal
verrouillage de phase montre trs bien que, dans un premier temps (3msec environ),
le PLL tente de s'accrocher au signal d'entre. Ds lors qu'il y est parvenu, il se
comporte comme un systme linaire.
2008 freddy.mudry@gmail.com
113
3.
3.5.1. Dmodulation FM
An de rendre les choses plus concrtes, on peut encore dcrire et observer les
signaux rencontrs dans une opration de modulationdmodulation (gure 3.11 ).
Ces signaux correspondent l'mission et rception de signaux FM. On y trouve :
1. Le message transmettre reprsent ici par une sinusode de frquence
fm =
1kHz.
2. Le signal envoy par l'metteur, de frquence centrale
message
signal mis
Uin, pll
Uout, pll
Uxor
message retard
signal dmodul
2
temps [msec]
Il est bien clair que les frquences d'un systme FM rel sont direntes puisque
les metteurs de radio FM transmettent les signaux audio (50Hz 20kHz) avec une
frquence centrale d'environ 100MHz.
114
2008 freddy.mudry@gmail.com
S OUR C E
F OLLOWE R
VC O in
P C A in
@ F R E QUE NC Y f
14
3
P HAS E
C OMPAR ATOR
P C B in
2 OR 13
P C 1out
OR
P C 2out
E XTE R NAL
LOWPAS S
F ILTE R
S F out
10
RSF
9
11
7
C IA
R1
VC O out
@ F R E QUE NC Y Nf = f
VC O
6
12
C IB
R2
CI
E XTE R NAL
N
C OUNTE R
R3
(a)
INP UT
OUTP UT
C2
2fC
OUTP UT
2 fL
R3 C2
Typically:
R3
R4 C2
R4
(R 3
C2
6N
N
2
f
fmax
100N
f R C
3, 000 ) C 2
4 2
fmax 2
f = fmax fmin
NOT E : S ometimes R 3 is s plit into two s eries res is tors each R 3 2. A capacitor C C is then placed from the midpoint to ground. T he value for
C C s hould be s uch that the corner frequency of this network does not s ignificantly affect n. In F igure B , the ratio of R 3 to R 4 s ets the
damping, R 4
(0.1)(R 3) for optimum res ults .
L OWPA S S F ILT E R
F ilter B
F ilter A
Definitions :
0.707
K KVCO
NR 3 C 2
N n
2K K V C O
1
R 3C 2S
F (s )
K KVC O
NC 2 (R 3 R 4 )
0.5
F (s )
n (R 3 C 2
N
)
K KVC O
R 3C 2S
1
S (R 3 C 2 R 4 C 2 )
Waveforms
P has e C omparator 1
P has e C omparator 2
V DD
P C A in
VS S
V DD
P C A in
VS S
V OH
V OH
P C B in
V OL
V OH
P C 1out
V OL
V OH
VC O in
P C B in
V OL
V OH
LD
V OL
V OH
P C 2out
V OL
V OH
V OL
VC O in
V OL
1
Note: for further information,
s ee:
=
V SSSon,
) New York, 1966.
O input
(1) fmin
F. G ardner,
P has eLock Techniques ,(VJ C
ohn
Wiley =
and
R 2 (C 1 + 32 pF )
(2) G . S . Mos
chytz, Miniature R C F ilters Us ing P has eLocked Loop, B S T J , May, 1965.
(3) G arth Nas h, P has eLock Loop Des ign F undamentals , AN535, Motorola Inc.
1
(4)fmax
A. B .=P rzedpels ki, P has eLocked
Articles
, AR
(V CDes
input
= V DD
) 254, reprinted by Motorola Inc.
+ fmin Loop
O ign
R 1 (C 1 + 32 pF )
F igure 3. G eneral P has eL oc ked L oop C onnec tions and Waveforms
Where: 10K
R1
1M
10K
R2
1M
100pF
C1
.01 F
2008 freddy.mudry@gmail.com
115
3.
Dtecteur de phase
Le dtecteur de phase I, dont on s'occupera ici, est une simple porte XOR recevant
deux signaux carrs (gure 3.13). Son signal de sortie est un signal rectangulaire
dont le rapport cyclique dpend de la dirence de phase entre les signaux d'entre.
La tension moyenne, obtenue par ltrage passe-bas, est donc une mesure directe de
ce dphasage
U3, moy =
VDD
(3.11)
R4 )
et d'une capacit
C2
(gure 3.14).
Caractristique du VCO
Le domaine de frquence dans lequel travaillera le VCO est choisi l'aide des composants externes
R1 , R2
et
C1 .
gure 3.15a . On y voit que pour une tension d'entre comprise entre
la frquence varie entre
fmin
et
fmax
et que
f0
Vmin
et
VDD ,
VDD .
A la lecture des donnes fournies par les fabricants, on constate que les relations
permettant d'valuer la caractristique du VCO sont trs variables. On peut cependant retenir les relations suivantes valables pour une alimentation unipolaire
VDD = 5 [V ]
fmin '
1
R2 (C1 + Cp )
Vmin ' 2 V
1
R1 (C1 + Cp )
(3.12)
(3.13)
10 k R1,2 1 M, 100 pF C1 10 nF et Cp ' 35 pF. On admet gnralement que l'oscillateur est au repos lorsque la tension d'entre vaut VDD /2. Il oscille
alors sa frquence de repos f0 .
avec
116
2008 freddy.mudry@gmail.com
u1(t)
A
u2(t)
XOR
14
u3(t)
C
t
0
T/2
R3
R3
R4
U3
C2
U4
U3
U4
C2
fvco
VDD
fmax
8
9
11
u4(t)
16
4
7
VCO
12
u2(t)
f0
R1
fmin
Vmin VDD / 2
C1
U4
(a)
0
R2
(b)
VDD
Fig. 3.15.: Caractristique d'un VCO (a) et son schma de ralisation (b)
1
u1
G1(s)
u3
G2(s)
u4
Kvco
2
u2
2008 freddy.mudry@gmail.com
117
3.
(t)
Z
(1 (t) 2 (t)) dt =
(t) =
(t)dt
G1 (s)
(s)
1
= [sec]
(s)
s
(3.14)
gain :
u3
K
(3.15)
rad
U4 (s)
G2 (s)
U3 (s)
u3 (t)
V
(3.16)
Kvco
u4
rad/sec
(3.17)
u1 (t),
u2 (t),
u3 (t),
u4 (t).
118
2008 freddy.mudry@gmail.com
phases 1 (t) et 2 (t) ou, leur dirence, le dphasage (t) = 1 (t) 2 (t).
Contrairement l'habitude prise en analyse frquentielle, il faut bien voir que ces
pulsations et phases sont des
les grandeurs indiques au-dessus des lignes du schma fonctionnel sont celles sur
lesquelles nous portons notre attention alors que les tensions mentionnes au-dessous
ne sont que leurs correspondants.
De plus, il est important de comprendre que le comparateur
et le gain
l'intgrateur
G1 (s)
G2 (s) =
1 + s/4
1 + sC2 R4
=
1 + sC2 (R3 + R4 )
1 + s/3
(3.18)
2 (s)
Go (s)
= G1 (s) K Kvco G2 (s)
(s)
rad/sec
Kbo K Kvco
(3.19)
rad/sec
Kbo
(3.20)
sec
2 (s)
1 1 + s/4
Go (s)
= Kbo
(s)
s 1 + s/3
rad/sec
rad/sec
(3.21)
Go (s)
2 (s)
Gf (s)
=
1 (s)
1 + Go (s)
rad/sec
rad/sec
2008 freddy.mudry@gmail.com
(3.22)
119
3.
Gf (s) =
Kbo 1+s/4
1+s/3
s
4
+ Ksbo 1+s/
1+s/3
Gf (s) =
=
s (1 + s/3 ),
on obtient
Kbo (1 + s/4 )
s (1 + s/3 ) + Kbo (1 + s/4 )
Kbo (1 + s/4 )
+
Kbo + s 1 + Kbo
4
Kbo ,
s2
3
on obtient la forme canonique
(1 + s/4 )
1
1
1 + s Kbo + 4 +
Gf (s) =
Gf (s)
1+
Gf (s) =
(3.23)
s2
Kbo 3
est de la forme :
s
4
1 + 2 sn +
s
n
(3.24)
2
4 ,
un d-
et son coecient
Kbo 3
1
1
1
2 =
= n
+
Q0
Kbo 4
n =
(3.25)
(3.26)
La rponse indicielle d'un systme dcrit par cette fonction de transfert dpend
fortement du coecient d'amortissement
pour direntes valeurs de
3.6.5. Calcul de 3 et 4
Dans l'tude gnrale des rgimes transitoires, on a vu que l'on choisit gnralement
compris entre 0.5 et 1 pour avoir un bon compromis entre un temps de monte
rapide et un temps d'tablissement court (gure 3.17). Le temps de rponse ou
d'tablissement du rgime permanent 5% prs vaut alors
t5% ' 3 bf =
120
3
n
2008 freddy.mudry@gmail.com
(3.27)
1.6
= 0.25
1.4
1.2
= 0.5
105%
1
95%
= 0.707
0.8
=1
0.6
0.4
0.2
5
6
temps * n []
10
Kbo
et du VCO,
Kbo K Kvco
(3.28)
sec
t5%
dpend
et
4 .
t5%
et
n '
Portant
3
t5%
(3.29)
dans les quations (3.25) et (3.26), on obtient alors les deux pulsations
recherches
3 =
4 =
2
1
n Kbo
n2
Kbo
1
=
(3.30)
n Kbo
2 Kbo n
2008 freddy.mudry@gmail.com
(3.31)
121
3.
VDD = 5 [V]
son domaine de frquences : 0 < f < 25 [kHz]
un temps de rponse de t5% ' 1 [msec]
un coecient d'amortissement ' 0.25 ou Q0 = 1/(2) ' 2.
sa tension d'alimentation :
1
u1
1
s
u3
u4
1+s4
1+s3
Kvco
2
u2
Circuit XOR
VDD
u1(t)
3 XOR
2
14
u3(t)
R3
8
9
11
u4(t)
16
4
7
VCO
12
R4
R1
R2
C2
122
2008 freddy.mudry@gmail.com
C1
u2(t)
2008 freddy.mudry@gmail.com
123
3.
3.7.4. Oscillateur
124
2008 freddy.mudry@gmail.com
2008 freddy.mudry@gmail.com
125
3.
126
t5% ' 1 ms
2008 freddy.mudry@gmail.com
PARAMETERS:
f min = 0
VDD = 5V
Kv co = 52.4k
U1A
0v
R1
VDD
fm
xor
3
in
7486
f iltre
55k
sin
IN+ OUT+
IN- OUT-
Kv co
VDD/3.5
VDD
E1
integr
R2
3k
1Meg
1Meg
V3
LP
VOFF = 0
VAMPL = 1
FC = 15K
MOD = 0
FM = 1K
out
0
R3
1k
EVALUE
sin(6.28*f min*time + V(integr))
C1
10n
XOR avec
adaptation de
niveau
Filtre
passe-bas
VCO
20
15
Uin, pll
10
Uout, pll
5
0
Uxor
0
0.5
1.5
2.5
3.5
2.5
2
ULP
1.5
1
0.5
0
0.5
1.5
2.5
3.5
1.5
ULP
1
Partie linaire normalise de la rponse indicielle
0.5
0
0.5
1.5
2
temps [msec]
2.5
3.5
2008 freddy.mudry@gmail.com
127
3.
uin (t).
(t)),
(t)
(t)
on a
Z
uvco (t) = A sin((t))
avec
(t)
(t) dt
0
0
uin (t))
(t) = 2 f (t) = 2 (fmin + Kvco
il vient
Z
(t)
Z
(t) dt = 2
t
0
(fmin + Kvco
uin (t)) dt
Z
(t) = 2fmin t + Kvco
uin (t) dt
(3.32)
0
C'est cette dernire quation qui est introduite dans la partie VCO de la simulation
du PLL o l'on a
E1
Evalue dont
VDD .
La rponse indicielle illustre par la gure 3.18 montre trs bien les tentatives d'accrochage successives jusqu'en
fa = fa2 fa1
128
f1 . Pour que
2008 freddy.mudry@gmail.com
f2
pourra varier dans un domaine plus large tout en restant synchronise avec
f1 .
fs = fs1 fs2
fs1 est la frquence de sortie du domaine
que fs2 est celle de sortie par le bas.
o
f2
domaine d
accrochage
f0
f1
0
fs2
fa1
f0
fa2
fs1
domaine de
synchronisation
Fig. 3.19.: Fonctionnement d'un PLL en mode non-synchronis
2008 freddy.mudry@gmail.com
129
3.
3.10. Exercices
Reg 1
temps
et
Ka
et de deux constantes de
dcrit par
1
1
(1 + s 1 ) (1 + s 2 )
Go (s) = Ka
1. Que valent
n ,
et le gain
Kbo
Gw (s) ;
b) que valent
c)
d)
Rp. :
Ka = 3;
Reg 2
Ka ' 1 /2
Go (s) = Ka
1
s (s + 10)
Ka
Kbf ?
Rp. :
= 0.59;
Reg 3
Ka = 71.5;
Kbf = 1;
Go (s) = Ka
1
1
(s + a) (s + 25)
que l'on boucle sur lui-mme avec une raction ngative, trouvez la valeur du gain
Ka
et la position du ple
> 0.59;
Reg 4
t5%
a => 35 [1/sec];
(s + z)
1
(s + a) s(s + 3)
que l'on boucle sur lui-mme avec une raction ngative, trouvez les valeurs de
Ka , z
et
d'tablissement
130
t5%
d'une seconde.
2008 freddy.mudry@gmail.com
3.10. Exercices
Reg 5
Kg = 2 V/1000
Ka = 5 [V/V],
1
1 + s mot
avec
Kmot = 5 [(rad/sec)/V]
mot
10 [msec]
Kbo .
w(t)
vaut
10 [V].
w(t)
e(t)
u(t)
Ka
Gmot(s)
u(t)
(t)
Kg
Reg 6
Ka
lorsque
K1 = 3
rad/sec
V
Kg = 2
mV
rad/sec
1 = 10 msec
Kg
Ka
K1
1+s1
1
s
Reg 7
2008 freddy.mudry@gmail.com
131
3.
I0 }. En dcrivant
Y (t) = Y0 + y(t)
I(t) = I0 + i(t)
F (Y (t), I(t)) = F (Y0 , I0 ) + F (y(t), i(t))
' F0 + k1 y(t) + k2 i(t)
En portant ce rsultat dans l'quation de Newton, on obtient
mY (t) = m
y (t) = mg + F0 + k1 y(t) + k2 i(t)
Comme autour du point de fonctionnement la force
F0
quilibre le poids
mg ,
on
m
y (t) = +k1 y(t) + k2 i(t)
I0 + i(t)
F(Y(t),I(t))
Y0 + y(t)
mg
Pour maintenir la sphre une hauteur donne, il faut la placer dans une boucle de
rgulation contenant un correcteur qui gnre le courant circulant dans la bobine.
Ce correcteur, de type proportionnel-driv, est dcrit par l'quation
i(t) = Kp
de(t)
e(t) + Td
dt
Ga (s) = Y (s)/I(s)
132
2008 freddy.mudry@gmail.com
et montrez
3.10. Exercices
2. Calculez la fonction de transfert du correcteur
Gc (s) = I(s)/E(s).
Gf (s) =
Kp est susamment grand. Remarque : un systme d'ordre 2 est stable si tous les coecients du dnominateur
G(s)
Kp
et
Td
de manire ce que le systme soit susamment rapide (n ) avec un amortissement satisfaisant ( ). Pour le voir, crivez
et montrez que
KP
et
Td
Gf (s)
manire suivante
Kp =
m n2 + k1
,
k2
Td =
et
de la
2n m
Kp k 2
m = 0.012 kg,
et que l'on souhaite avoir
Kp
et
Td
k1 = 12
N
m
n = 20 rad/sec
et
k2 = 0.66
= 0.5,
N
A
t5% ?
PLL 0
avec
Z2 /Z1
Z2 = R 2
et
Z1 = R1 //C .
Pour les exercices qui suivent, on considre un circuit PLL de type 4046B
prend
PLL 1
On applique un circuit XOR deux signaux carrs de mme priode 200 [s]
dcals de 20 [s].
1. Dessinez le signal de sortie du XOR ; que valent le rapport cyclique et la priode
de ce signal ? Que vaut le gain du comparateur XOR ?
2. Le circuit XOR est reli un ltre retard de phase ralis avec
R3 = 45 [k],
R4 = 5 [k] et C2 = 5 [nF]. Que vaut la tension moyenne que l'on peut mesurer
C2 ? Que vaut la
PLL 2
C2 .
calculez son
2008 freddy.mudry@gmail.com
133
3.
PLL 3
Kbo
vaut
105 [1/sec].
Dessinez son schma fonctionnel puis calculez les trois composants pour que ce
circuit ait un temps de rponse de
100 [sec]
0.5.
PLL 4
3.3 [k]
Nosc ?
PLL 5
f2 (t)
un saut de frquence
f1 (t),
que valent
ttrans
et
f2 (t).
et 150 kHz.
1. Dessinez la caractristique du VCO puis calculez les trois composants ncessaires son fonctionnement.
2. Le PLL reoit en entre un signal carr dont la frquence varie entre 80 et
120 kHz ; que vaut la tension d'entre du VCO ?
PLL 6
et un
' 0.5.
Simulez votre circuit avec une frquence d'entre de 30 kHz. Combien de temps
faut-il au PLL pour qu'il entre dans sa zone de fonctionnement linaire ? Son temps
de rponse correspond-il celui souhait ?
PLL 7
On veut raliser un circuit PLL dont la frquence de sortie est deux fois
0.5.
134
2008 freddy.mudry@gmail.com
3.10. Exercices
Quelques rponses :
PLL 1
1
2
PLL 2
PLL 3
PLL 4
2
3
4
PLL 5
1
2
PLL 6
VCO
ltre
Spice
PLL 7
VCO
ltre
K = VDD /
= 0.25 msec
0
Kvco = 4.59 [kHz/V]
R3 = 21 [k], R4 = 6.7 [k]
Kvco = 20.9 103 [(rad/sec)/V]
t5% = 0.18 [msec]
Nosc ' 1.4
R1 = 10 [k], R2 = 20 [k]
2.9 [V ] Uvco 4.1 [V ]
R1 = 25 [k], R2 = 100 [k]
R3 = 85 [k], R4 = 7.6 [k]
ts ' 0.46 [ms]
R1 = 50 [k], R2 = 25 [k]
R3 = 170 [k], R4 = 15 [k]
0
0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
0.1
0.2
0.3
0.4
0.5
0.6
temps [msec]
0.7
0.8
0.9
2008 freddy.mudry@gmail.com
135
3.
136
2008 freddy.mudry@gmail.com
Bibliographie
Cours de rgulation automatique, eivd-iAi, 2004
Michel Girard, Boucles verrouillage de phase, McGraw-Hill, 1988
[1] M. Etique :
[2]
[3] MicroSim Application Notes, version 8.0, june 1997, pp. 303-305
137
Bibliographie
138
2008 freddy.mudry@gmail.com