Polycopié
Polycopié
Polycopié
G11(p)
E(p) +
G12(p)
+
+
G2(p)
Y(p)
H(p)
Chapitre
s(t)
La tension e(t) reprsente le signal dentre et la tension s(t) aux bornes du condensateur C
reprsente le signal de sorite. La relation entre le signal e(t) et le signal de sortie s(t) peut tre
dcrite par lquation diffrentielle suivante :
( = )
Cours LAA311
)(
+ )(
Enseignant : O. GUENOUNOU
(I. 1)
1
Impulsion de Dirac
Cours LAA311
Enseignant : O. GUENOUNOU
Cours LAA311
Enseignant : O. GUENOUNOU
Rflexion
Effet de
laction
Action
Observations
Cours LAA311
Enseignant : O. GUENOUNOU
Sortie
Chaine directe
(Ou daction)
Chaine de retour
(Ou dobservation)
Entre
Erreur
Correcteur
Actionneur
Ecart
Systme
dynamique
Sortie
Capteur
Cours LAA311
Enseignant : O. GUENOUNOU
(t )
e(t )
s(t )
(t )
e(t )
s(t )
Cours LAA311
Enseignant : O. GUENOUNOU
15
10
1.5
5
1
0
0.5
0
-5
0
-10
0
1.5
Systme instable
0.5
10
Cours LAA311
Enseignant : O. GUENOUNOU
Chapitre
Transforme de Laplace
II
On note :
f(t)
().
TL
(II. 1)
= +
Original
Image
si |(|). existe do :
On dmontre quil existe un nombre tel que lintgrale soit absolument convergente pour
> et ne le soit pas pour < .
Ce nombre sappelle labscisse de convergence absolue ou abscisse de sommabilit.
II.1.2 Exemple :
-
0
() =
1
<0
>0
(II. 2)
Cours LAA311
Enseignant : O. GUENOUNOU
= )(
1
1
1
= = 0 =
{ = )(u(t)} =
( = )sin ( )
= )(
sin( ).
> )( 0
( )
( ) ( )
1
( )
= )(
+
2 ( ( ) + )
1
1
1
= )(
=
2 ( ) ( + )
+
{ = )(sin ( = })
> )( 0
1. Linarit
})({ = )(
{. ( )+ . ( = }). )(+ . )(
})({ = )(
2. Drivation temporelle
{ })({ = )(( = }). )( (0)
Plus gnralement
})({ = )( () =
)(
(II. 3)
(II. 4)
(II. 5)
(II. 6)
lim ( = )lim . )(
5. Translation (a> 0)
Cours LAA311
Enseignant : O. GUENOUNOU
(II. 7)
(II. 8)
(II. 9)
9
6. Thorme du retard
Soit f(t) =0 pour t< 0 (fonction causale)
({ })({ = )( )} = ()
(II. 10)
})({ = )(( )= . (. )
Et plus gnralement
({) ()} =
()
()
Soit :
Alors
= )( )( = )(
(). ( )
})({ = )( })({ = )(
(II. 11)
(II. 12)
(II. 13)
(II. 14)
(II. 15)
Pour retrouver loriginale dune fonction F(p) donne, on dcompose cette fonction
(en gnral, une fraction rationnelle en p) en lments simples dont en prendra loriginal dans
la table de transformation.
Exemple
() =
Cours LAA311
( + )( + )
( = ){ = })( {})(
Enseignant : O. GUENOUNOU
10
= )(
1
1
+ 2 + 3
1
1
( = )
+ 2
+ 3
( = )exp(2 )+ exp(3)
>0
()
)(
Chacune des valeurs de p qui annule soit le numrateur, soit le dnominateur, sappelle une
racine. Chacune des valeurs de p qui annule le numrateur est appele un zro. Chacune des
valeurs de p qui annule le dnominateur est appele un ple.
Lide consiste alors mettre lexpression du dnominateur sous forme dun produit de
facteurs o apparat chacune des racines :
( = )( )( ) . ( ) ( )
Ceci permettra de dcomposer F(p) sous la forme dune somme de fractions simples. Dans
tous les cas qui vont suivre nous considrerons que N(p) est dordre polynomial infrieur
D(p). Si ce nest pas le cas, il faudra dabord extraire la partie entire par division
polynomiale. Puis traiter uniquement la partie fractionnaire.
Exemple :
() =
+ +
+
() = + + +
partie entire
partie fractionnaire
()
=
+
+ +
(II. 16)
( ). . . ( ) ( ) ( )
( )
( )
Cours LAA311
Enseignant : O. GUENOUNOU
11
= ( ) )(
(II. 17)
Exemple :
>0
+ 1
6
+ 1
+ 1
=
6 ( 2)( + 3)
+ 1
=
+
+
( 2)( + 3)
2 + 3
+ 1
1
=
( 2)( + 3)
6
+ 1
3
= lim ( 2) ( = )lim
=+
(+ 3)
10
+ 1
2
= lim ( + 3) ( = )lim
=
( 2)
15
= lim = )(lim
4) Exprimer la solution en p:
+ 1
11 3 1
2 1
(= )
=
+
( 2)( + 3)
6 10 2 15 + 3
5) Exprimer la solution dans le domaine temporel :
1 3 2
= )( +
6 10
15
pour t>0 ( ne pas oublier)
()
=. . +
+
+. . +
+.. (II. 18)
( ) ( )
( )
( )
. . ( ) . .
Cours LAA311
Enseignant : O. GUENOUNOU
12
( ) )(
( )!
(II. 19)
Exemple :
4 + 4
=
+ +
+
( 2) ( 1)
( 2)( 1)
4 + 4
=2
( 2)( 1)
1
4 + 4
[
= lim
( = ])lim
= 3
( 2)( 1)
1!
= lim ( = )lim
4 + 4
= 1
( 1)
4 + 4
= lim ( 1) ( = )lim
= 1
( 2)
= lim ( 2) ( = )lim
4 + 4
3 2
1
1
( = )
= +
( 2)( 1)
2 1
( = )3 + 2 exp(2 ) )(,
>0
Cours LAA311
Enseignant : O. GUENOUNOU
13
c) Systme invariant
Un systme est dit invariant sil obit la loi suivante :
)(
)(
()
()
)
+
(
=
+ + ()
Soit encore :
( ) + + +
=
= )(
)(
+ + +
Cours LAA311
Enseignant : O. GUENOUNOU
(II. 20)
14
H(p) est la fonction de transfert du systme (transmittance), elle sexprime par le rapport de
deux polynmes en p construit partir des coefficients de lquation diffrentielle.
E(p)
S(p)
H(p)
h(t)
(t)
0
h(t)
t
-Chaque systme est caractris par une rponse impulsionnelle h(t) qui est la rponse du
systme limpulsion.
-Lexpression de la sortie s(t) correspondante une entre e(t) scrit :
= )( ()( ) = ( )()
= )( ()( ) = ()( )
Nous avons :
= )( ()( ) = ( )()
()
= {()( = })
)(
Cours LAA311
()( = ). {(})
Enseignant : O. GUENOUNOU
(II. 21)
15
Note :
La fonction de transfert H(p) dun systme est la transforme de Laplace de sa rponse
impulsionnelle h(t).
Exercices sur le chapitre 2
Exercice N01
Quelles sont les transformes de Laplace des signaux suivants :
Impulsion de Dirac
fonction Echelon
fonction Rampe
-1
fonction Signe
0
T
Porte causale
T
2T
Toit causal
2T
Exercice N02
Trouver les transformes de Laplace des fonctions suivantes ( laide de la table).
f3 t sin(2t) ,
f1 t t 2 3t 4 u t ,
f 2 t t 3 exp(t) ,
f5 t sin(2t) exp(t) ,
f6 t cos(t) exp(3t) ,
f4 t cos(t)
Exercice N03
Donner la transforme de Laplace dune fonction borne 0, et priodique 0, T
pour t> 0 (on considre T ).
Cours LAA311
Enseignant : O. GUENOUNOU
16
fT t
.
0
2T
T
Fonction Priodique
priodiques suivantes :
5
A
T/2
3T/2
2T
4 5
10
Onde Triangulaire
Onde Carre
Exercice N04
Retrouver les originales des fonctions rationnelles en p suivantes
F p
p 1
p2 p 1
, F p
p 1 p 1
2
, F p
6 p3
p 1 p 1 p 3
, F p
p 1 p 1
2
e 3 p
Exercice N05
Soit un systme de 2me ordre rgi par lquation diffrentielle suivante :
..
y (t ) 5 y (t ) 6 y (t ) sin( 2t )
.
y (0 ) 2
..
et y (0) 1
y (t ) y (t ) 6 y (t ) 12t 20
.
y (0) 0
Cours LAA311
et y (0) 0
Enseignant : O. GUENOUNOU
17
Exercice N06
-Calculer pour chaque montage ci-dessous, la fonction de transfert F(p).
R
e(t)
Cours LAA311
R1
s(t)
e(t)
Enseignant : O. GUENOUNOU
R2
s(t)
18
Chapitre
III
INTRODUCTION
Il existe un nombre important de systmes rels (Circuits RC, RLC, systme Masse
Ressort, etc) qui peuvent tre dcrits par des simples quations diffrentielles du premier et
du second ordre. La complexit dun systme est en ralit est d la multitude de soussystmes dordre infrieur ou gal 2 qui le composent. Par exemple un systme de troisime
ordre peut tre dcompos, en trois sous-systmes du premier ordre ou en un systme du
premier ordre et un systme du second ordre. Il est donc trs important de comprendre et de
maitriser les comportements et les caractristiques des systmes du premier et du second
ordre.
III.1 SYSTEME DU PREMIER ORDRE
Un systme est dit du premier ordre sil est rgi par une quation diffrentielle de
premier ordre de la forme :
)(
+ )( = )(
(III. 1)
En supposant que les conditions initiales soient nulles (CI=0), il est possible de calculer la
fonction de transfert G(p) du systme en appliquant la transforme de Laplace aux deux
membres de l'quation (III.1):
[ )( (0)] + ()( = )
( [)+ 1] = )(
()
=
)(1 +
= )(
1 +
(III. 2)
K : Gain statique
: Constante de temps
Cours LAA311
Enseignant : O. GUENOUNOU
19
(= )
()( = )()( = )
1 +
(III. 3)
1 +
= )(
= )(
1
+
(III. 4)
Point darriv :
Cours LAA311
[ ]
[ ]
Enseignant : O. GUENOUNOU
20
0,37
0,05
ds(t)/dt
0,37 0,05
s(t)
+
0
0
)(
{ = })(
(= )
=
1 + + 1
+
=
1
1
+
= )( 1 0
(III. 5)
Cours LAA311
[ ]
Enseignant : O. GUENOUNOU
21
Point darriv :
lim = )(lim = )(
[ ]
s(t)
ds(t)/dt
)(
=
0
(III. 6)
-Il est important de remarquer que plus est faible, plus le systme est rapide.
Influence de la constante de temps sur la rponse indicielle
1
0.9
=1s
=3s
0.8
=4s
0.7
=2s
s(t)
0.6
0.5
0.4
0.3
0.2
0.1
0
10
15
Temps[s]
Figure III.3: Influence de la constante de temps () sur la rponse dun systme du premier
ordre
Cours LAA311
Enseignant : O. GUENOUNOU
22
( = )
=
1 + + 1
(= )
+
+
=
+ +
1
+ 1
( = ) + >0
(III. 7)
Enseignant : O. GUENOUNOU
(III. 8)
23
En supposant que les conditions initiales soient nulles, il est possible de calculer la fonction
de transfert G(p) du systme en appliquant la transforme de Laplace aux deux membres de
l'quation (III.8):
1
2
[ )( (0)] + ()( = )
)( (0) (0)+
= )(
avec :
1
2
+
+ 1()( = )
1 2
+ + 1
(III. 9)
Pour > 1
4
( 1)
=
1
Pour =
= )(
2
20
+ 1 = (1 + 1 ()1 + 2 )
1
1
=
= )(
2
20
+ 1 = (1 + )2
Cours LAA311
Enseignant : O. GUENOUNOU
24
La figure ci-dessous donne la position des ples dans le plan complexe en fonction de
la valeur du coefficient d'amortissement.
Comme dans le cas du systme du premier ordre, lentre applique est un chelon
damplitude E0
= )(
(1 + 1 ()1 + 2 )
(= )
(1 + 1 ()1 + 2 )
La rponse temporelle s(t) correspondante secrit :
1
= )( 1
>0
(III. 10)
avec :
2
1
+ =
Etude de s(t)
lim = )(
() =
() = 0
Puisque > ; alors :
1
1
<
>
() > 0 pour >0, s(t) est monotone
d'o quelques que soient les valeurs de et ,
croissante sur , donc pas de dpassement. La rponse est alors qualifie dapriodique
puisquelle ne prsente aucun dpassement relativement la valeur finale. Plus le facteur
damortissement est grand, plus le temps de rponse est consquent.
() = 0
Cours LAA311
Enseignant : O. GUENOUNOU
25
)( 1
>0
(III. 11)
(= )
avec
(1 + )2
0
La rponse temporelle s(t) correspondante secrit :
= )([1 (1 + ) ]
>0
(III. 12)
Etude de s(t)
() = 0
lim = )(
() =
() > 0 >0 ,
() = 0
Cours LAA311
Enseignant : O. GUENOUNOU
26
(= )
avec:
1 2 2
(2 + 2 0 + 20 )
2 + + 1
0
0
+ 2 + = ( )( )
= + 1
= 1
1
avec :
cos () =
sin() = 1
(III. 13)
Enseignant : O. GUENOUNOU
(III. 14)
27
La pseudo-pulsation
= 1
(III. 15)
=
qui se produisent aux instants :
(III. 16)
(III. 17)
1
Le premier extremum, valeur maximale de s(t) est appel le dpassement (exprim en %)
linstant tpic
D% = 100
=
(III.18)
100
90
80
70
D%
60
50
40
30
20
10
0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
Coefficient d'amotissement
0.8
0.9
Le temps de monte tm est dfini comme l'instant pour lequel la valeur finale est atteinte pour
la premire fois.
( ) =
(III. 19)
La rsolution de l'quation (III.19), nous donne :
=
arctan
(III. 20)
Enseignant : O. GUENOUNOU
28
Soit :
( )( = ) ( ),
()( = ). (),
()( = ). ()
alors :
(= )
)(
()
)(
)(
=
1 + )(. )(
)(1 + )(. )(
= )(
)(
1 + )(. )(
(III. 21)
Ou encore :
= )(
)(
)(
=
1 + )(. )(1 + )(. )(
Enseignant : O. GUENOUNOU
29
= )(
()( = ). )(. )(
()
= )(. )(
)(
)( = )( = )(. )( = )( . )(
alors :
= )(
=
)(
)(
=
1 + )(. )(1 + ()
)(
1 + )(
(III. 22)
)(
1 + )(
Remarque :
Tous les systmes boucls retour non unitaire peuvent se ramener sous forme de systmes
retour unitaire.
Cours LAA311
Enseignant : O. GUENOUNOU
30
Systmes en srie
Systmes en parallle
E(p)
G(p)
S(p)
Cours LAA311
Enseignant : O. GUENOUNOU
31
G(p)
1
G(p)
S(p)
Z(p)
Ces quivalences sont gnralement utilises pour rorganiser des schmas-blocs qui
prsentent des difficults pour calculer la fonction de transfert.
III.3.2 Systmes plusieurs entres
Principe de superposition
Un systme dynamique est linaire si le principe de superposition peut tre appliqu.
Ainsi la rponse s(t) dun systme linaire d plieurs entres appliques simultanment est
gale la somme des rponses de chaque entre applique sparment.
Exemple :
Soit calculer la sortie S(p) du systme, deux entres, reprsent ci-dessous :
Enseignant : O. GUENOUNOU
32
Schma quivalent
Schma quivalent
B(p)
S2(p)
G(p).H(p)
(= )
)(
)(
1 + )(. )(
( = )( )+ (= )
(= )
1
)(
1 + )(. )(
)()(+ )(
1 + )(. )(
1
1
1
, F2 ( p )
et F3 ( p ) 2
2
p 1
10 p 11 p 1
p p 1
Figure 1
Figure 2
Figure 3
Exercice N02
- Dterminer la fonction de transfert du circuit ci-dessous et montrer quil sagit dun systme
du premier ordre.
- Pour C=2F et R1=R2 =1M, dterminer le gain statique K et la constant de temps .
- Donner lexpression de la rponse indicielle.
Cours LAA311
Enseignant : O. GUENOUNOU
33
Exercice N03
- Dterminer la fonction de transfert du systme mcanique (figure a) et montrer quil sagit
dun systme du second ordre. Quand une force F=2N est applique linstant t=0, la masse
se met osciller comme illustr par la figure (b).
2.5
1.5
0.5
0.5
1.5
2.5
3.5
Exercice N04
La figure (a) donne la rponse indicielle (chelon unitaire) du systme boucl de la figure (b).
-
Cours LAA311
Enseignant : O. GUENOUNOU
34
(a)
(b)
Exercice N05
Simplifier les schmas bloc suivants puis calculer les fonctions de transfert Y(p)/X(p) :
Cours LAA311
Enseignant : O. GUENOUNOU
35
Exercice N06
Trouver lexpression de la sortie en fonction des diffrentes entres pour les systmes
suivants:
D(p)
F1(p)
E(p)
H1(p)
H2(p)
S(p)
E2(p)
+
E1(p) +
-
F1(p)
F2(p)
H1(p)
S(p)
H2(p)
E3(p)
Cours LAA311
Enseignant : O. GUENOUNOU
36
Chapitre
IV
Analyse frquentielle
Introduction
Lanalyse frquentielle
(= )
()
)(
o E(p) et S(p) sont les transformes dans le domaine de Laplace des fonction temporelles :
e(t) et s(t).
On appelle fonction de transfert complexe ou transmittance, la fonction obtenue en
remplaant la variable de Laplace p par le terme jw (imaginaire pur).
()
(= )
(IV. 1)
) (
Cours LAA311
Enseignant : O. GUENOUNOU
37
Diagramme damplitude
log( )
()
00
log( )
Diagramme de phase
= +
= +
Il suffit donc dajouter les diagrammes des fonctions T1 et T2 aussi bien sur le diagramme
damplitude que sur le diagramme de phase.
Dfinitions
Gain en dcibels (dB)
() = 20 log()
(IV. 2)
Exemples :
() = 2 () = 20 log(2) = 6
() = 10 () = 20 log(10) = 20
Octave
Cest un intervalle dont les deux bornes sont distantes dun rapport de 2.
Exemples : [1 2], [10 20] et [5 10].
Dcade
Cest un intervalle tel que ses bornes sont distantes dun rapport de 10.
Exemples : [1 10], [10 100] et [5 50].
Gain statique
Gain lorsque 0 , celui-ci doit tre fini.
Cours LAA311
Enseignant : O. GUENOUNOU
38
soit :
() =
1 +
1 +
= |)(| =
1 + ()
Cours LAA311
Enseignant : O. GUENOUNOU
(IV. 3)
(IV. 4)
(IV. 5)
39
(IV. 6)
La courbe du module est approxime par une droite horizontale (asymptote horizontale)
dquation :
= 20 log
lim ( ) = lim (20 log 10 log(1 + ()))
Nous avons :
lim (1 + ()) = () (on nglige le nombre 1 devant )
= 20 log 20 log
On dit que la pente est de -20dB par dcade (le module diminue de 20dB pour une
augmentation d'un facteur 10 de ). Cest aussi quivalent une pente de -6dB par octave (le
module diminue de 6dB pour une augmentation d'un facteur 2 de ).
Les deux asymptotes se croisent pour =c, cette pulsation est appele pulsation de cassure :
1
=
(IV. 7)
= 20 log
= 20 log 20 log
1
= 1 =
Pour =c
|( )| =
=
(IV. 8)
1 + 2
= 20 log 20 log 2 = 20 log 3
Pour la pulsation de cassure l'cart par rapport l'asymptote est de -3dB.
Pour les pulsations double et moiti de c , l'cart est de -1dB.
Diagramme de phase
La phase de la fonction complexe T(j) est :
() = arctan()
(IV. 9)
Le diagramme asymptotique a la forme dune marche descalier, il nest pas suffisamment
prcis pour reprsenter lvolution de la phase.
Quelques points de la courbe
( ) = arctan(1) = 45
1
= arctan = 26,56
2
2
(2 ) = arctan(2) = 63,43
Cours LAA311
Enseignant : O. GUENOUNOU
40
20 log( K )
c
de
B/
0d
-2
C 1/
log( )
C 1/
log( )
(1 )
=
=
1 +
1+
1+
1 +
> 0 [0, +[
1
+
= ()=
< 0 [0, +[
1 +
Point de dpart (w=0)
= lim =
= ()=
(IV. 10)
(IV. 11)
= lim = 0
= lim = 0
Cours LAA311
Enseignant : O. GUENOUNOU
41
c
Figure IV.6 : Diagramme de Nyquist dun systme du premier ordre
IV.3.3 Diagramme de Black
On retrouve sur le trac les points caractristiques dfinis prcdemment. Le lieu a une
asymptote horizontale aux basse frquences (=0) et une asymptote verticale aux hautes
frquences ( ).
( = 0) = 20 log
( = 0) = 0
( ) =
20 log( K )
( ) = 90
0
( ) = 20 log 3
c
)
( = 45
1 2
+ + 1
(= )
(IV. 12)
2
2 +
+ 1 = (1 + 1 ()1 + 2 )
1 + 1 (1 + 2 )
La reprsentation prend la forme ci-dessous :
() =
Cours LAA311
Enseignant : O. GUENOUNOU
(IV. 13)
42
1 2
20 log( K )
d
B/
0d
-2
ec
2
- 40
log( )
/de
dB
c
log( )
Points du diagramme
( = 0) = 20 log
( = 0) = 0
( ) =
( ) = 180
( ) = 90
Pour = 1
= )(
2
2 +
+ 1 = (1 + )2
20 log( K )
log( )
log( )
Enseignant : O. GUENOUNOU
43
Points du diagramme
( = 0) = 20 log
( = 0) = 0
( ) =
( ) = 180
( ) = 90
Pour < 1
Les racines sont complexes conjugues = + 1 et =
Posons =
= |)(| =
(IV. 14)
1 + 2
() =
1 + 2
(1 ) + 4
= ([1 ) + 4 ] = ()
(0) = 1 , lim () = 0
1
() = [4(1 ) + 8][(1 ) + 4 ]
2
() = 2 [ (1 2)][(1 ) + 4 ]
() = 0 = 0 (1 2) = 0
car (1 2)
Cours LAA311
Enseignant : O. GUENOUNOU
44
0
f(x)
| |)(K
<
+
+
-
||)( 20log(K)
X
f(x)
||)(
0
0
-
+
+
-
1 2
= 1 2
+
0
K
||)(
0
20 log
20log (K)
Rsonance
( ) = 20 log
(IV. 16)
21
Facteur de surtension Q
Le facteur de surtension exprim en dB est la diffrence entre la valeur du gain la rsonance
( ( ) ) et le gain statique (0):
1
= 20 log
20 log = 20 log
21
21
Le facteur de rsonance est alors :
1
=
21
Cours LAA311
Enseignant : O. GUENOUNOU
(IV. 17)
(1.18)
45
Asymptotes
( = 0) = 20 log . La droite = 20 log est une asymptote horizontale.
Nous avons :
= |)(| =
= ([1 ) + 4 ]
(1 ) + 4
Le terme 1 +
()
2(2 1) 1
= 1 +
+
= 20 log = 20 log( ) 40 log
Diagramme de phase
La phase de la fonction de transfert est :
2
arctan
si 1 > 0
() =
(IV. 19)
2
+ arctan
si
1
<
0
1
Nous avons :
( = 0) = 0 , ( ) = et ( = ) =
On peut vrifier que la fonction () est ngative sur R+. La fonction () est monotone
dcroissante, do le tableau ci-dessous :
0
0
Cours LAA311
Enseignant : O. GUENOUNOU
46
AdB max
1
2
20 log( K )
1
2
log( )
log( )
[1 2]
() =
=
1 + 2 (1 ) + 4
(IV. 20)
[1 ]
= )(=
(1 ) + 4
2
(IV. 21)
= lim = 0
= lim = 0
= ( 1) =
=
(1 ) + 4 2
Les cordonns du point dintersection du lieu avec laxe imaginaire est 0, et ce point
nous avons = .
Cours LAA311
Enseignant : O. GUENOUNOU
47
1,5
0
0.5
0.8
0
0
(0 )
Cours LAA311
Enseignant : O. GUENOUNOU
48
Soit :
= Arg[T(j)]
posons
= Arg[F(j)]
()
1 + ()
et
et
(IV. 22)
A = |T(j)|
B = |F(j)|
1 + A + 2A cos
sin
si A + cos > 0
arctan
A
+
cos
sin
= arctan
si A + cos < 0
A + cos
si A + cos = 0
2
() et B() sont des fonctions des deux variable A et , paramtres en .
(IV. 23)
(IV. 24)
Tracer le lieu de T(j) dans le plan de BLACK (utiliser les coordonnes rectangulaires
de labaque pour reporter les points)
Lire sur les courbes isom odules et isophases de labaque les valeurs des points
appartenant F(j).
Cours LAA311
Enseignant : O. GUENOUNOU
49
40
0 dB
30
0.25 dB
0.5 dB
20
1 dB
-1 dB
3 dB
10
-3 dB
6 dB
-6 dB
-10
-12 dB
-20 dB
-20
-30
-40
-40 dB
-360
-315
-270
-225
-180
-135
-90
-45
Cours LAA311
Enseignant : O. GUENOUNOU
50
20
20
(1 + )
15
, (= )
, (= )
, (= )
(1 + )
(1 + 10)
(1 + 10()1 + 100)
+
30
50
15
, = )(
, = )(
1.6 + 1
+ + 1
(1 10()1 + 100)
Exercice N02
Soit un systme asservis S caractris une fonction de transfert en boucle ouverte T(p).
La figure ci-dessous reprsente le lieu asymptotique de Bode de la fonction T(p).
Sachant que le systme S est un systme dphase minimal, dterminer lexpression de T(p)
T ( jw) dB
Exercice N0 03
Soit le systme dcrit par la fonction de transfert F(p) :
1
(= )
(1 + ()1 + 10)
1. Dterminer la valeur du coefficient damortissement , de la pulsation naturelle n , et
du gain statique Ks.
2. Tracer les trois diagrammes (Bode, Nyquist et Black) de T(p).
Exercice N0 04
Reprsenter le diagramme de Nyquist de la fonction T(p)
(= )
10
(1 + 0.5()1 + 5)
Exercice N0 05
Soit la fonction de transfert T(p) avec a R+
(= )
1 +
1 +
Cours LAA311
Enseignant : O. GUENOUNOU
51
Chapitre
V.1 Stabilit
V.1.1 Dfinition :
Un systme est dit stable si sa sortie tend vers zro lorque lentre sannule.
+ = )(()
)(
+ +
(V. 1)
(V. 2)
+ +
= )(
+
+ + +
(V. 3)
1.
( = ) e
Cours LAA311
Enseignant : O. GUENOUNOU
(V. 4)
52
p 0
p 0
2. ( = ), . ()! e
p 0
4. (( = ))! ,sin + e
Enseignant : O. GUENOUNOU
53
Conclusion :
Un systme linaire est strictement stable si tous les ples de sa fonction de transfert sont
partie relle ngative
V.1.2 Stabilit des systmes en boucle ferme
Soit le systme en boucle ferme reprsent ci-dessous avec G(p) et H(p) sont
respectivement les fonctions de transfert de la chane directe et de la chane de retour.
()
)(
)(
=
=
)(
1 + )( )(
1 + ()
1 + ( = )0
(V. 5)
Si T(p) est la FTBO associe F(p), le systme est stable si seulement si lquation
1+( = )0 na que des solutions (racines) parties relles ngatives.
La stabilit de la BO nimplique pas ncessairement celle de la BF et inversement.
La connaissance des ples de la fonction de transfert dun systme (BO, BF) permet de
vrifier aisment la stabilit de ce dernier. Cependant il est souvent difficile de calculer les
ples notamment ceux des systmes boucls.
Une alternative consiste donc tudier la stabilit sans passer par le calcul des ples. On
distingue deux approches :
Approche algbrique ;
Approche gomtrique.
Cours LAA311
Enseignant : O. GUENOUNOU
54
(V. 6)
Condition 1.
Une condition ncessaire de stabilit est que les coefficients ai du polynme soient
tous prsents et tous de mme signe. Cette condition est gnralement vrifie pour la
plupart des systmes physiques
Condition 2.
Dans le cas o la condition 1 est vrifie, on construit le tableau de Routh partir des
coefficients du polynme.
Une condition ncessaire et suffisante pour que le systme soit stable est que tous les
coefficients de la premire colonne soient de mme signe.
Remarque :
Le nombre de changement de signe est gal au nombre de racines partie relle positive
V.3.1. 1 Cas particuliers
1ier Cas :
Il sagit dexaminer le cas o tous les termes dune ligne sont nuls. Cela signifie que p = jw
est racine de D( p). Le systme est donc juste oscillant et lon peut crire :
Cours LAA311
( = )( + )()
Enseignant : O. GUENOUNOU
(V. 7)
55
Puisque tous les coefficients de D(p) sont prsents et de mme signe, on construit alors la
table de ROUTH.
11
18
2
2
18
18
0
0
0/2
0/0
0/0
18
On constate, aprs avoir construit la table de Routh, que les lments de premire colonne
sont de mme signe.
= )( + 2 + 11 + 18 + 18 = ( + 9)( + 2 + 2) = 0 = 3
Il sagit du cas o le pivot dune ligne non nulle est gal zro.
Exemple : Etudier la stabilit de la fonction de transfert F(p) admettant comme
dnominateur :
= )( + 5 + 4 + 20 + 1
Puisque tous les coefficients de D(p) sont prsents et de mme signe, on construit alors la
table de ROUTH.
Cours LAA311
Enseignant : O. GUENOUNOU
56
5
0/
20
1
0
0
<0
Il y a deux changements de signe, alors D(p)=0 admet deux racines partie relle positive. Le
systme est alors instable.
Sous MATLAB il possible de calculer les racines de D(p) et ce en utilisant les deux
instructions suivantes :
%%%%%%%%%%%%%
D=[1 5 4 20 1]
roots(D)
%%%%%%%%%%%%
Les ples sont : -4.9931, 0.0218 + 1.9918i, 0.0218 - 1.9918i et -0.0505. Il y a effectivement
deux ples partie relle positive.
V.1.4. Critres gomtriques (Graphiques)
Lobjectif de ces critres est dobtenir des renseignements sur la stabilit dun systme en
boucle ferme partir du lieu de sa fonction de transfert en boucle ouverte T(p).
V.1.4.1. Critres du revers
Les critres du revers sont valables pour les systmes dphase minimal, cest dire les
systmes pour lesquels la fonction de transfert en boucle ouverte T(p) na ni de ples ni de
zros partie relle positive.
V.1.4.1.1. Critre du revers dans le plan de NYQUIST
Un systme dphase minimal sera stable en boucle ferme si le lieu de transfert de
NYQUIST (de T(p)) parcouru dans le sens des pulsations () croissantes laisse le point -1
(point critique) sur sa gauche.
Enseignant : O. GUENOUNOU
57
AdB
AdB
A >0dB
0dB
c0
Log(w)
0dB
c0
Log(w)
c0
0dB
Log(w)
A <0dB
Stable
180
Oscillant
Log(w)
180
Instable
180
Log(w)
Log(w)
180
-180
-180
-180
180
Cours LAA311
Enseignant : O. GUENOUNOU
(V. 8)
(V. 9)
58
=P-Z
Cours LAA311
Enseignant : O. GUENOUNOU
59
( + )( + )( + )
Cours LAA311
Enseignant : O. GUENOUNOU
60
Limage du segment [0, +[, est le trac de la boucle ouverte T(p) dans le plan de
Nyquist dcrit dans le sens des croissants.
Limage du segment ], 0] est symtrique de la prcdente par apport laxe rel.
Lmage du demi-cercle de rayon infini est rduite un point (en gnral lorigine) si le
systme est physiquement ralisable.
0
R
0
0
j
Figure V.13: Contour de Nyquist en prsence de ple lorigine ou de ples imaginaires purs
Des exemples de ples lorigine seront tudis durant la sance de TD.
V.1.5. Conditions de stabilit pratique
Un systme la limite de la stabilit est mal amorti. Son bon fonctionnement n'est pas
assur car une faible modification de ses caractristiques peut le rendre instable. Les lieux de
transfert peuvent tre obtenus par modlisation ou exprimentalement. Quelle que soit la
mthode utilise, ces lieux ne sont pas connus de faon exacte.
Cours LAA311
Enseignant : O. GUENOUNOU
61
dB
0dB
dB
180 o
log( )
log( )
= + ( )
= ( )
(V. 10)
(V. 11)
Enseignant : O. GUENOUNOU
62
dB
0dB
180 o
).
(,
= 20 log() , =
Cours LAA311
Enseignant : O. GUENOUNOU
63
V.2. PRECISION
Considrons le systme asservi ci-dessous:
B(p)
(p)
E(p) +
G1(p)
+
G2(p)
S(p)
H(p)
Figure V.17: Systme asservi avec perturbation
Ce systme asservi est conu afin de faire suivre la sortie s(t) une loi dtermine par l'entre
e(t). La qualit d'un tel systme se jugera :
par sa stabilit
par la prcision avec laquelle la loi est suivie
Afin d'valuer la prcision de ce systme, on dfinit tout instant l'cart (ou l'erreur) () .
Le cas idal de point de vue prcision est d'avoir chaque instant () = . Cependant
l'entre e(t) peut varier au cours du temps et le systme est soumis divers perturbations.
l'cart n'est pratiquement jamais nul.
Remarques:
o Minimiser ou annuler l'cart quant l'entre e(t) varie au cours du temps, c'est
rsoudre un problme de poursuite.
o Minimiser ou annuler l'cart quant l'entre e(t) est constante, c'est rsoudre un
problme de rgulation.
V.2.1 Etude avec perturbation nulle "Prcision vis vis l'entre"
Dans ce cas-l, le schma de la figure ci-dessus se simplifie et on obtient le schma-bloc cidessous:
Cours LAA311
Enseignant : O. GUENOUNOU
64
Expression de l'cart ()
E(p)
1 + T(p)
E(p)
1 + T(p)
(V. 12)
E(p)
1 + T(p)
= lim = )(lim
E(p)
1 + T(p)
(V. 13)
(V. 14)
l'expression de est :
= lim
= lim
1 + T(p)
1 + T(p)
Cours LAA311
1 + T(p)
= lim
Enseignant : O. GUENOUNOU
(V. 15)
65
( = )
1 + + + + +
(1 + + + + + )
(V. 16)
()
/ (0) = (0) = 1
)(
= lim
= lim
1 + T(p)
= lim
)(
= lim
)( + )(
()
1+
)(
(V. 17)
1+
= 0
Enseignant : O. GUENOUNOU
66
( = ). . = )( )(
= lim
= lim
1 + T(p)
p[1 + T(p)]
)(
= lim
= lim
)( + )(
()
p 1 +
)(
= lim
+
(V. 18)
= +
=
= 0
Pour annuler l'erreur de position, il faut au moins deux (2) intgrations dans la
fonction de transfert boucle ouverte T(p) (2)
Si T(p) est de classe 1, l'erreur de vitesse diminue quand le gain K augmente.
Cet cart est dfini pour une entre sous forme parabolique.
( = ). . = )( )(
2
Enseignant : O. GUENOUNOU
67
= lim
(V. 19)
= +
=
= 0
Pour annuler l'erreur de position, il faut au moins deux (3) intgrations dans la
fonction de transfert boucle ouverte T(p) (2)
Si T(p) est de classe 2, l'erreur en acclration diminue quand le gain K augmente.
Tableau rcapitulatif
( = )G(p). H(p) =
classe 0
=0
Ecart de vitesse
1+
Ecart en acclration
Ecart de position
()
avec (0) = (0) = 1
)(
classe 1
classe 2
classe 3
=1
=2
=3
Cours LAA311
Enseignant : O. GUENOUNOU
68
V2.2. Etude avec entre nulle "Prcision vis vis des perturbations"
Dans ce cas-l, le schma de la figure devient :
G(p)
B(p)
1 + G(p)H(p)
(V. 20)
(V. 21)
Il suffit donc d'avoir (p) = 0 pour que la sortie du systme ne sera pas affecte par la
perturbation.
Perturbation en chelon:
La perturbation b(t) est un chelon ( = )= )( )(
L'cart en rgime permanent d la perturbation est gale :
= lim p
Notons:
( = )
()
()
= lim p
H(p)G(p)
B(p)
1 + G(p)G(p)H(p)
et ( = )()
()
()
H(p)G(p)
. 2 . 2 (). 1
= +
1 2 1 (). 2 ( )+ 1 . 1 (). 2 . 2 ()
1 + G(p)G(p)H(p)
=
. 2 . 1
1 +2
Cours LAA311
+ 1. 2
Enseignant : O. GUENOUNOU
(V. 22)
69
= 0
. 2
1 + 1. 2
> 0
= 0
= 1
> 0
Conclusion
En rponse une perturbation en chelon, un systme comportant au moins une intgration en
amont du point d'entre de cette perturbation, aura un cart nul.
Ainsi il est clair que linsertion dun correcteur doit tre faite juste derrire le
comparateur.
Perturbation en rampe
Par un raisonnement similaire aux prcdents, on aurait:
. 2 . 1 1
1 +2 + 1 . 2
(V. 23)
Rsultat :
= 0
= 1
(p) =
(V. 24)
E(p)
1 + T(p)
Enseignant : O. GUENOUNOU
70
(t) =
( [1 + ( )])
1 + ( )
(V. 25)
existe une limitation des signaux d'entre en vitesse et en acclration, il possible d'obtenir
une condition caractrisant le prcision exige.
Le signal d'entre est alors dfini par :
(V. 26)
Parmi les signaux satisfaisant ces conditions, le signal sinusodal est de loin le plus simple.
Soit ( = ) ( )un tel signal.
d'o :
<
1 + ()
1 + ()
On obtient:
|1 + (|)
<
Notons que cette condition est ncessaire et peut dans certains cas n'tre pas suffisante.
Cours LAA311
Enseignant : O. GUENOUNOU
71
V.3 Rapidit
V.3.1 Rponse transitoire
Considrons un systme peu amorti dun ordre suprieur deux (approxim par deux
ples complexes conjugus dominants), soumis une entre chelon, dont la rponse est
donne ci-dessous:
s(t)
D1
+5%
sfinale=KE0
-5%
tm
tr5%
Cours LAA311
Enseignant : O. GUENOUNOU
72
FdB max
FdB 0
-3dB
c 3dB
c0
c 3dB
c 0
log( )
()
=
1 + )(
K
1 p
()
()
+1
= )(
=
=
=
1 + )( + )(+ 1 + 1 +
+1
Enseignant : O. GUENOUNOU
73
20 log(
K
)
K 1
c
de
B/
0d
-2
C
log( )
1 K
Remarque: : La bande passante est dfinie sur le diagramme de gain de la fonction F(p) et
non sur le diagramme de gain de la T(p) utilise pour ltude de la stabilit.
Il est alors ncessaire de raliser un lien entre les bandes passantes de la FTBO et de la FTBF.
Dune manire gnrale, on montre que la bande passante -3 dB dun systme en BF peut
tre approche par la pulsation de coupure 0 dB de sa FTBO (si K>>1).
Dans le cas particulier du systme du 1ier ordre, plus est faible, plus tr5% =3, le temps de
rponse 5% est faible et donc plus le systme est rapide. La rponse harmonique du
systme du 1ier ordre, montre que plus est faible, plus la bande passante 0,
est
leve. On peut dduire que le systme est dautant plus rapide que la bande passante est
leve. Cest le cas pour tous les autres systmes physiques.
Cours LAA311
Enseignant : O. GUENOUNOU
74
( = )()()
+ 5 + 8 + 5
(1 + )
(+ 1)(0.5 + 1)
Un systme boucl est caractris par une fonction de transfert en boucle ouverte T(p) :
(= )
(+ 3)
Exercice N03
Soit le systme en boucle ferme de la figure avec :
T ( p)
K
( P 1) 3
u(t)
T ( jw)
e(t)
T(p)
y(t)
K
(1 jw) 3
Cours LAA311
Enseignant : O. GUENOUNOU
75
Exercice N04
Soit le schma-bloc dun systme asservi de la figure ci-dessous o e(t) est lentre et b(t) une
perturbation mesurable.
b(t)
(p)
e(t) +
10
(p+3)
10
(p+1)
s(t)
-En appliquant le critre de Routh, montrer que le systme est toujours stable pour K positif.
- Pour K=1, tracer le diagramme de Bode de la fonction de transfert en boucle ouvert T(p).
- Dterminer la valeur du Gain K correspondant une marge de phase de 45 degrs.
- Calculer en fonction de K les erreurs de position et de vitesse.
II. Pour amliorer la prcision du systme boucl, on ajoute comme indiqu par la figure cidessous un compensateur de fonction de transfert = )(
(.)
Cours LAA311
Enseignant : O. GUENOUNOU
76
VI. Correction
Chapitre
VI
Correction
INTRODUCTION
Au cours du chapitre prcdent nous avons tudi les performances dun systme asservi
savoir la stabilit et la prcision. Nous avons montr que pour annuler lerreur statique, la
FTBO doit comporter un ou plusieurs intgrateurs et doit prsenter un gain statique K le plus
lev possible. Par contre pour obtenir un bon degr de stabilit il faut une marge de phase et
une marge de gain suffisantes, do un gain statique faible.
Il est clair que le respect des spcifications du cahier de charges ne peut tre en gnral
obtenu par un simple rglage dun gain ou dinsertion dintgrateurs. Dans la plus part des
cas, il est indispensable dinsrer un correcteur (compensateur, rgulateur).
VI.1 Structures de correction
VI.1.1 Correction srie
Ce type de correcteur est insr dans la chane directe juste aprs le comparateur, il est en
srie avec le systme et dlivre un signal de commande
()( = )
(VI. 1)
Enseignant : O. GUENOUNOU
(VI. 2)
77
VI. Correction
()
Gnralement la fonction
approximatives de
des perturbations.
Cours LAA311
()
= )(
()
1
()
(VI. 3)
Enseignant : O. GUENOUNOU
78
VI. Correction
La sortie S(p) peut tre considre comme la superposition de deux schmas blocs cidessous :
(= )
Cours LAA311
)(
1 + ()()()
Enseignant : O. GUENOUNOU
(VI. 4)
79
VI. Correction
( = ))((= )
( = )( )+ (= )
)()(()()
1 + ()()()
(VI. 5)
)(
)()(()()
(VI. 6)
1
()()
(VI. 7)
On notera que, dans ce cas, le systme asservi suit parfaitement la loi dentre. Comme le cas
de la compensation des perturbations, souvent la fonction
() ()
ralisable. La compensation ne sera donc pas parfaite dans les rgimes transitoires.
VI.2 Classes des correcteurs
On peut dfinir deux classes de correcteurs:
1. Les correcteurs spcifiques qui sont dtermins partir de la fonction de transfert en
boucle ferme recherche. Il s'agit de dterminer la fois la structure et les paramtres
du correcteurs ( structure+paramtres)
2. Les correcteurs classiques, dont la structure est bien dfinie et dont seuls les
paramtres sont dterminer.
VI.3 Correcteurs classiques (sries)
On distingue trois actions de base : Proportionnelle (P), Intgrale (I) et Drive (D),
qui permettent, individuellement, damliorer telle ou telle performance. Elles sont simples
raliser mais, en gnral, dgradent dautres performances.
Type daction
Signal de
Commande m(t)
Proportionnelle
(P)
Intgrale (I)
Drive (D)
( = )()
(= )
1
()
( = )
()
Fonction de
Effet
transfert
= )(
= )(
= )(
Cours LAA311
Enseignant : O. GUENOUNOU
80
VI. Correction
(VI. 8)
1
Td
Inconvnients
-
un
correcteur
P.D,
laction
drive
ntant
pas
physiquement
ralisable
(Deg(N)>Deg(D)) on adopte :
= )( 1 +
Cours LAA311
1 + (+ )
=
1 +
1 +
Enseignant : O. GUENOUNOU
81
VI. Correction
1 +
avec a > 1
(VI. 9)
1 +
avec a > 1
1
1
et =
1 1 /a
m
(VI. 10)
(VI. 11)
2 1 /
1 1 /a
2 1 /
1 + ()
1 + ()
()
1
= 0 =
Cours LAA311
Enseignant : O. GUENOUNOU
(VI. 12)
82
VI. Correction
log( ) = log
1
1
1
1 1
1
1
1
= log = log = log + log
2
2
2
1
= [log( ) + log( )]
2
1 +
1 + 1/
( 1)
1 + 1 1/
4 + ( 1)
1 + 1/
( 1)
( + 1)
1
+ 1
2 + ( 1)
+ 1
1
= arcsin
+ 1
Notons que :
|( )| =
(VI. 13)
1
4 + ( 1) =
+ 1
1
1 + sin
( + 1) sin = 1 =
+ 1
1 sin
=
19,5
30
1 + sin
1 sin
(VI. 14)
9
10
11
12
13
59
L'intrt de ce correcteur est visible sur son diagramme de phase : il permet de remonter la
phase la pulsation , d'o il est possible d'amliorer la marge de phase du systme
(VI. 15)
Cours LAA311
Enseignant : O. GUENOUNOU
(VI. 16)
(VI. 17)
83
VI. Correction
1 + ()
(VI. 18)
(VI. 19)
() = 90 + arctan()
log
1
Ti
log
Ce correcteur admet une forme approche ds que l'intgration n'est pas pure:
avec a>>1
= )(1 +
(1 + )
(1 + )
Enseignant : O. GUENOUNOU
84
VI. Correction
= )(
(1 + )
(1 + )
(VI. 20)
avec a > 1
Les mmes calculs que pour le correcteur avance de phase sont utiliss pour tudier
le comportement frquentielle de ce correcteur.
Pour la pulsation =
Notons que :
1 sin
1 + sin
(VI. 21)
19,5
30
10
11
12
13
59
1
a
(t)
+
1
Ti p
m(t)
Td p
Figure VI.9: Structure d'un correcteur PID parallle
Un correcteur P.I.D thorique peut tre dfini par une fonction de transfert de la forme:
Cours LAA311
Enseignant : O. GUENOUNOU
85
VI. Correction
= )(+
avec T > T
(VI. 22)
(1 + )
avec T > T
(VI. 23)
avec T > T
(VI. 24)
Comme l'action drive n'est pas ralisable physiquement, un PID rel peut tre obtenu en
1 + +
+ =
1 + +
= )(
C(p) est une fraction rationnelle dont les racines du numrateur sont relles si:
T 4T
En considrant l'quation (IV.26) , C(p) peut tre mis sous une forme :
Cours LAA311
Enseignant : O. GUENOUNOU
(VI. 25)
(VI. 26)
86
VI. Correction
= )(
avec:
=
1 + 1
(1 + )(1 + )
On note que T + T = Ti et T T = Ti T
1 1
(VI. 27)
4
1
T
K
Ti
1 / T1
et =
1
T
(VI. 28)
1
KTd
1 / T2
1
TiTd
Cours LAA311
Enseignant : O. GUENOUNOU
87
VI. Correction
E(p) +
S(p)
Systme
Tosc
-Noter Kosc la valeur du plus petit gain permettant de maintenir les oscillations.
-Noter Tosc la priode des oscillations.
- Rgler les paramtres du correcteur comme indiqu par le tableau ci-dessous
PI
PID
Srie
Ti
Td
Parallle
Mixte
,
, .
.
,
Cours LAA311
Enseignant : O. GUENOUNOU
88
VI. Correction
PI
PID
Ti
Td 0
3,3 T1
2 T1
0,5 T1
(1 + ()1 + ( )1 + )
avec
Cours LAA311
Enseignant : O. GUENOUNOU
(VI. 29)
89
VI. Correction
(1 + ()1 + ( )1 + )
(VI. 30)
T(p)
1 + T(p)
Avec le critre mplat, on cherche obtenir pour le systme rgl la bande passante la plus
large possible en vitant que le facteur de rsonance dpasse 1. La rponse harmonique en
boucle ouverte (T(p)) correspondante cet objectif est donne par la figure ci-dessous.
-1
1
1
2T p
1
TP
log(w)
-2
1
2T p(1 + Tp)
(VI. 31)
Pour obtenir donc T(p) T(p), on aligne la (ou les) pulsation(s) de coupure du correcteur
(1/T1 , le cas chant 1/ T2) sur la (ou les) pulsations de coupures du systme rgler (1/1 , le
cas chant 1/ 2).
La fonction de transfert du rgulateur est alors :
= )(
(1 + ()1 + )
(1 + ()1 + )
(1 + ()1 + ( )1 + )
T(p) =
Cours LAA311
(1
Enseignant : O. GUENOUNOU
+ )
(VI. 32)
(VI. 33)
90
VI. Correction
1
=
2
= 2
(1 + ( = )1 + Tp) T = < 1
T =
(VI. 34)
(VI. 35)
(= )
1 +
. )(
1 +
1 +
()
1 +
Si la FTBO ne possde pas dintgration, lcart de position ne sera pas nul, il faut
prvoir soit un intgrateur dans le correcteur, soit un gain suffisant.
Enseignant : O. GUENOUNOU
91
VI. Correction
fonction FTBO corrige, note 0c . Le problme est que lon ne connat que la pulsation de
coupure 0dB de la fonction FTBOnc non corrige, note 0. Choisissons une valeur de 0c
lgrement suprieure 1,50.
correcteur, suppos bien plac est . Le cahier des charges impose gnralement la marge
1 + sin
1 sin
, il faut que
Etape 5 : Vrification
Cours LAA311
Enseignant : O. GUENOUNOU
92
VI. Correction
( p)
Figure 1
1. A partir du schma, exprimer la fonction de transfert C(p) du correcteur sous la
forme :
= )(
()
= +
+
()
1 + 1 +
Avec = 20 10 . et = 80 10 .
avec :
= )(
Cours LAA311
10
,
(1 + 0,1)
Enseignant : O. GUENOUNOU
>0
93
VI. Correction
1. En appliquant le critre de Routh, montrer que le systme reste toujours stable pour
K>0.
2. Dterminer la valeur de K pour avoir une marge de phase de 45.
3. Quelle est la condition sur K pour avoir une erreur de vitesse infrieure 0,01
4. Pour K=10, tracer le diagramme de Bode de la FTBO.
5. Peut-on avoir au mme temps une marge de phase de 45 et une erreur de vitesse
infrieure 0.01
6. Pour avoir une marge de phase de 45 et une erreur de vitesse infrieure 0.01, nous
introduisant un correcteur avance de phase dans la chane directe. Dterminer les
paramtres du correcteur.
7. On dcide dutiliser un correcteur retard de phase, dterminer les paramtres de ce
correcteur.
Exercice 03 (En relation avec le critre mplat)
Soit le systme asservi retour unitaire de la figure ci-dessous:
= )(
1
1 +
Avec <
1. Donner lexpression de F(p) FTBF, puis dduire les valeurs du gain, de la pulsation
propre et du coefficient damortissement.
2. Donner la condition pour avoir un systme boucl sans rsonnance.
3. Donner la condition pour avoir un systme oscillatoire amorti, dans ce cas donner
lexpression du premier dpassement et la sortie s(t) (lentre e(t)= 5u(t)).
4. Tracer le diagramme de Bode de la FTBO (Il faut tenir compte de la condition <
)
Cours LAA311
Enseignant : O. GUENOUNOU
94
VI. Correction
avec :
10
>0
(1 + 2()1 + 0,1)
En utilisant le critre mplat, donner les correcteur correspondant.
= )(
En utilisant la mthode de Ziegler-Nichols, donner les paramtres des trois correcteurs (P, PI
et PID) correspondants au systme dont la rponse indicielle est donne ci-dessous :
10
9
8
reponse du systme
7
6
5
4
3
2
1
0
Cours LAA311
Enseignant : O. GUENOUNOU
95