TP Automatique Licence
TP Automatique Licence
TP Automatique Licence
T.P. Automatique
Licence EEA
Année Universitaire 2003/2004
3
O. Horn, E. Losson, M. Kreutner et Y. Morère TP Automatique Licence EEA
4
TP 1 : Etude de la régulation en
vitesse d’un moteur
1.1 Objectif du TP
Identifier un système réel du deuxième ordre par sa réponse indicielle et sa réponse harmo-
nique.
J
f
v ω
5
O. Horn, E. Losson, M. Kreutner et Y. Morère TP Automatique Licence EEA
d’où :
kV (p) − (R + Lp)Cr (p)
Ω(p) = (1.1)
(f + Jp)(R + Lp) + k 2
1.3 Manipulation
L’étude est faite sur un banc de régulation de vitesse Tergane 25. Ce banc est constitué d’un
moteur à courant continu et d’une électronique associée permettant une commande en boucle
ouverte ou une régulation en boucle fermée.
Les appareils de mesures nécessaires sont :
– un oscilloscope double trace à mémoire,
– un générateur de signaux basse fréquence.
kv0 − Rc0
ω0 = (1.2)
f R + k2
6
O. Horn, E. Losson, M. Kreutner et Y. Morère TP Automatique Licence EEA
1
T = f
gΩ(f )
a A
V (f )
x
Fig. 1.2 – Mesure du gain et du déphasage avec un oscilloscope
|F (f )| = Aa
Tracer sur papier semi-logarithmique le diagramme de Bode de F en gain (20log|F (f )|) et
de phase (Arg(F (f ))). A partir des asymptotes du diagramme en gain donner les valeursde K
et τ (voir figure 1.3).
Indications théoriques :
K L−1 Kt −t/τ
−→ e u(t)
(1 + τ p)2 τ2
· µ ¶ ¸
K L−1 t
−→ K 1− 1+ 2 e −t/τ
u(t)
p(1 + τ p)2 τ
la seconde transformée de Laplace correspond à la réponse indicielle du système de la première
transformée. Donc, si V (t) = u(t) soit V (p) = p1 , alors :
lim gΩ(t) = K
t→∞
7
O. Horn, E. Losson, M. Kreutner et Y. Morère TP Automatique Licence EEA
20Log10 (K)
−40 dB/dec
1
1 f= 2πτ 10 100 f en Hz
0
−90
Phase (en deg)
−180
– à t = 3τ , gΩ(t) = 0.8K
– à t = 4τ , gΩ(t) = 0.9K
Comparer les valeurs de ξBF , ωnBF et KBF mesurées avec les valeurs théoriques obtenues en
F (p)
calculant la FTBF : F (p) = 1+F (p)
.
Quel est l’apport sur les performances du systèmes, lors du passage de la boucle ouverte à la
boucle fermée ?
8
O. Horn, E. Losson, M. Kreutner et Y. Morère TP Automatique Licence EEA
K
0.91K
0.8K
Amplitude
3τ 4τ temps en s
Amplitude
KBF
tm temps en s
9
Cette page est laissée blanche intentionnellement
TP 2 : Etude de l’asservissement de
position d’une charge
2.1 Objectif du TP
Mettre au point un asservissement de position par identification des paramètres du système
et introduction d’un correcteur à avance de phase.
Moteur
R, L, J, f, k
Réducteur de
vitesse 1/N
Poulie de renvoi
charge
Potentiomètre
Poulie Motrie (d) x(t) de recopie
11
O. Horn, E. Losson, M. Kreutner et Y. Morère TP Automatique Licence EEA
gAkkx d L J
En posant : K = N (f R+k2 )
, Te = R
, Tm = 2 , on obtient :
f + kR
K
F (p) =
p(1 + Te p)(1 + Tm p)
En constatant de plus que Te (constante de temps électrique) est très faible devant Tm
(constante de temps mécanique), on peut approcher la FTBO par :
K
F (p) =
p(1 + Tm p)
Cette approximation est admissible en automatique si Te < Tm /5.
En admettant ce modèle de l’asservissement en position, on se propose dans ce TP d’identifier
cette fonction de transfert, puis de la corriger afin d’en améliorer les performances.
2.3 Manipulation
2.3.1 Etude qualitative
Sur le banc de travail Tergane connecter tous les éléments de la chaı̂ne directe comme
l’indique le synoptique dessiné sur la maquette. Mettre le préamplificateur sur l’entrée 10.
Observer la réponse du système en regardant le déplacement de la charge :
– en envoyant un échelon d’entrée dont on fera varier l’amplitude et pour différentes valeurs
de l’atténuateur.
– en envoyant des sinusoı̈des à l’entrée (à l’aide d’un GBF) de différentes fréquences et de
différentes amplitudes.
Décrire qualitativement les résultats obtenus et les justifier.
Pour identifier la FTBF, il suffit donc de determiner ξ et ωn qui sont en relation avec le
temps de réponse et le dépassement.
12
O. Horn, E. Losson, M. Kreutner et Y. Morère TP Automatique Licence EEA
d
Amplitude
tm Temps (sec)
– Dépassement : D% = d/M
– Temps de montée : tm
Quel est l’effet de la variation de g sur le dépassement D%, et le temps de monté tm (donner
la valeur limite de g au delà de laquelle le système est instable2 .
Identification
Placer l’atténuateur sur 0.5 (milieu de l’echelle).
– à partir des mesure de D et tm identifier le système (valeur de ξ et ωn ) au moyens des
abaques.
– en déduire la fonction de transfert en boucle ouverte du système (K et Tm ).
avec T = 2.5ms et a = 4
Il peut être réaliser par le circuit passif de la figure 2.4.
Avec pour valeur :
– condensateur C = 1µF
– resistance R1 = 10kΩ
– resistance R2 = 3.3kΩ
2
il faudra certainement passer le préamplificateur à 100.
13
O. Horn, E. Losson, M. Kreutner et Y. Morère TP Automatique Licence EEA
R1
Entrée Sortie
Correcteur C R2 Correcteur
14
TP 3 : Etude et commande d’un
processus thermique
3.1 Objectif du TP
Étudier le fonctionnement d’une régulation en température et déterminer empiriquement
les paramètres d’un correcteur PID.
3.2 Présentation
Il s’agit d’étudier un système comprenant un processus thermique et son matériel de com-
mande (correcteur PID). Principe de l’appareil : l’air est aspirée dans l’atmosphère par une
pompe centrifuge qui le souffle sur une source de chaleur, traverse l’enceinte cylindrique puis
est rejeté dans l’atmosphère.
Le processus (figure 3.1) consiste à élever la température soufflée dans l’enceinte à la valeur
désirée (set value) comprise entre la température ambiante et 60˚C, et le rôle de la régulation
est de mesurer la température de l’air dans le cylindre, de la comparer avec la température
désirée et d’engendrer un signal de commande qui module la puissance fournie à la résistance
chauffante.
L’étude du système sera faite autour d’une position d’équilibre fixée à 40˚C. Les échelons de
température que l’on appliquera seront de ±5˚C qui correspond à ±2V injecté sur l’entrée D
du schéma. Le système est en boucle ouverte si les points x et y ne sont pas reliés. Si les points
x et y sont connectés, on peut établir une commande proportionnelle (bande proportionnelle)
ou introduire un circuit correcteur entre les point A et B.
3.3 Manipulation
Tout d’abord brancher le système et le laisser chauffer pour atteindre 40˚C (la températures
mesurée doit être égale à la température affichée).
Par la suite, toutes les mesures de transitoires se feront dans le sens “froid → chaud ”.
Attention : Un échelon de tension négatif provoque une variation de tempé-
rature positive
15
O. Horn, E. Losson, M. Kreutner et Y. Morère
Sonde de
Température
Ventilateur
de chauffe
filament
Fig. 3.1 – Schéma Bloc du Principe de l’asservissement de température. Processus
W
Bande
Proportionnelle Signal de Sortie
Alimentation en %
x y Pont
Mesure
16
C Valeur affichée
A B D INT
Commande EXT TS
Commande Dépassement Externe Niveau de Valeur
Chauffage Perturbation affichée
O. Horn, E. Losson, M. Kreutner et Y. Morère TP Automatique Licence EEA
e−τ p
F (p) = K
1 + Tp
17
O. Horn, E. Losson, M. Kreutner et Y. Morère TP Automatique Licence EEA
18
TP 4 : Stabilité des systèmes asservis
4.1 Objectif du TP
Simuler la réponse d’un système asservi afin de déterminer ses performances, puis calculer
un correcteur en fonction de critères de stabilité et vérification de la correction.
4.2 Introduction
Soit le système asservi suivant :
10
G(p) = ¡ p 3
¢ (4.1)
1+ 5
19
O. Horn, E. Losson, M. Kreutner et Y. Morère TP Automatique Licence EEA
1 1 + aT p
C(p) = √ , avec : a>1 (4.2)
a 1 + Tp
20
O. Horn, E. Losson, M. Kreutner et Y. Morère TP Automatique Licence EEA
µ ¶
a−1
Φm = arcsin (4.4)
a+1
4. Relever également la réponse indicielle du système asservi afin d’en déduire le temps de
montée et le dépassement.
Diagrammes de Bode
Gain (dB)
20log(a)
c
de
B/
0d
+2
20log( a )
0
Phase (deg)
φm
0
1 1 1 pulsation (rad/sec)
aT aT T
1+aT p
Fig. 4.2 – Diagramme de Bode de 1+T p
.
21
O. Horn, E. Losson, M. Kreutner et Y. Morère TP Automatique Licence EEA
1 1 + aT p
C(p) = √
a 1 + Tp
µ ¶
a−1
φm = arcsin
a+1
1 + sin φm
a =
1 − sin φm
1
ωcp = √
aT
1
T = √
ωcp a
20log|G(j ω)|
1
20log
a
20log|CG(j ω)|
ωcp
ω
– M φG : marge de phase de G,
– M φCG : marge de phase de CG (système corrigé),
– M φCG = M φG + φm .
22
O. Horn, E. Losson, M. Kreutner et Y. Morère TP Automatique Licence EEA
ωcp
ΜφCG
Arg(G(j ω))
Arg(CG(j ω))
ΜφG
−π
23
Cette page est laissée blanche intentionnellement
TP 5 : Performances des systèmes
asservis et calcul de correcteurs P, PI
et PID.
5.1 Objectif du TP
Visualiser l’action des correcteurs à action proportionnelle et à action intégrale en se basant
sur des critères de précision, de rapidité ou encore de stabilité, puis faire le calcul de correcteurs
PI et PID à l’aide du tracé du lieu d’Evans du système régulé.
5.2 Introduction
Soit le système asservi de la figure 5.1 :
Entrée Sortie
C(p) G(p)
1
G1 (p) =
1+p
1
G2 (p) =
(1 + p)(1 + p2 )(1 + p3 )
1
G1 (p) =
1+p
25
O. Horn, E. Losson, M. Kreutner et Y. Morère TP Automatique Licence EEA
1 ki
C(p) = =
Ti p p
K
F T BF (p) = 2ξ p2
1+ ωn
p + 2
ωn
√
avec ωn = ki et ξ = 2√1ki .
– Pour quelles valeurs de ξ et donc de ki , est-il possible d’observer un dépassement sur la
réponse indicielle du système corrigé en boucle fermée.
Le premier dépassement observé sur la réponse indicielle est donné par la formule :
à !
−ξπ
d% = exp p
1 − ξ2
On peut vérifier grâce à cette équation qu’une valeur de ξ = 0.45 permet d’obtenir un
dépassement d’environ 20%.
– Calculer la valeur de ki pour laquelle le dépassement est de 20%. Observer la réponse
indicielle obtenue pour cette valeur de ki vérifier la valeur du dépassement et relever
l’erreur de position, le temps de montée et le temps de maintien à 2%.
– Observer l’évolution du dépassement, du temps de montée et du temps de maintien en
fonction de la valeur de ki .
– Quels sont les avantages et les inconvénients du correcteur à action intégrale ?
1
G2 (p) =
(1 + p)(1 + p2 )(1 + p3 )
26
O. Horn, E. Losson, M. Kreutner et Y. Morère TP Automatique Licence EEA
27
Cette page est laissée blanche intentionnellement
Notice d’utilisation du programme
auto1
29
O. Horn, E. Losson, M. Kreutner et Y. Morère TP Automatique Licence EEA
NYQUIST
En appuyant sur ce bouton, on obtient la visualisation de la représentation de Nyquist du
système.
BLACK
En appuyant sur ce bouton, on obtient la visualisation de la représentation de Black du
système.
NB : pour ces deux derniers types de représentations en fréquence, les marges
peuvent également être visualisées tout comme dans le diagramme de Bode.
Représentations temporelles
REPONSE INDICIELLE
En appuyant sur ce bouton, on peut visualiser la réponse du système lorsqu’on lui soumet un
échelon en entrée. En cliquant sur le bouton droit de la souris lorsque le curseur est situé dans
la fenêtre graphique (cadre blanc contenant un graphe) on voit apparaı̂tre le menu suivant :
– Plot type
– Systems
– Characteristics
– Zoom
– Grid
Si vous cliquez sur Characteristics vous voyez apparaı̂tre quatre options :
– Peak Response
– Settling Time
– Rise Time
– Steady State
En cliquant avec le bouton gauche de la souris sur Peak Response, un cercle plein permet de
localiser la valeur maximale du signal sur la réponse. En cliquant sur ce cercle, vous pouvez
visualiser les coordonnées du point en question. En cliquant sur Settling Time, il est possible
d’obtenir le temps de maintien défini ici par le temps nécessaire pour que la réponse se stabilise
à plus ou moins 2% de la valeur finale.
Rise Time donnera le temps de montée défini ici par le temps écoulé lorsque le signal passe de
10% à 90% de sa valeur finale.
Les échelles de temps et de fréquence peuvent être modifiée en utilisant l’option Tools du menu,
puis Reponse preference.
REPONSE IMPULSIONNELLE
En appuyant sur ce bouton, on peut visualiser la réponse du système lorsqu’on lui sou-
met une impulsion de Dirac en entrée. Les mêmes paramètres que ceux définis précédemment
peuvent être obtenus sur cette représentation.
NB : il sera possible de visualiser les coordonnées de points sur les courbes en
cliquant à l’endroit souhaité avec le bouton gauche de la souris.
30
O. Horn, E. Losson, M. Kreutner et Y. Morère TP Automatique Licence EEA
stable ou non.
NB : si le système est instable les pôles ont nécessairement des valeurs réelles
positives.
CALCUL DU CORRECTEUR
En appuyant sur ce bouton, on peut déterminer la valeur du gain du correcteur en fonction
de critères de stabilité du système (cf. Remarque particulière concernant cette option à la fin
de la notice)
31