Cours: AUTO 3 Asservissement Et Régulation

Université des Frères Mentouri Constantine 1
INSTITUT DES SCIENCES ET TECHNIQUES APPLIQUEES

Département Productique Mécanique et Industrialisation
Cours : AUTO 3
Asservissement et régulation
PMI S4
Année universitaire 2019/2020

Cours : Asservissement et Régulation 2019-2020 2ème Année PMI_ISTA
L’automatique : conduite sans intervention humaine des systèmes.

Plusieurs branches de l’automatique :
a) Automatique de système linéaire : ⤇ Déterministe.
⤇ Stochastique (aléatoire).
⤇ Monovariable.
⤇ Multivariable
⤇ Continus
⤇ Echantillonnés.
b) Automatique de système non linéaire.
 Continus : automatisation en analogique ; avec des composants électroniques,
résistances, capacités, amplificateurs.
Automatisation par calculateur : microprocesseur, micro ordinateurs
Système non linéaire (la même diversité que dans le système linéaire)
Notre intérêt dans ce module est l’étude des systèmes linéaires monovariables déterministe
continus.
Problème de l’automatique :
Moteur
- à vide 3000 tr/min

- charge 50Kg : 1500 tr/min
- charge 100Kg : 1000 tr/min
On remarque que la vitesse varie avec la charge.

Souvent on veut garder la vitesse constante quelle que soit la charge : il faut alors construire un
système automatique.
Système
automatique Moteur
de réglage
Dr. ASSABAA .M Page 1

Chapitre I : Position du problème

I. Introduction
Ce cours est une introduction à l’automatique des systèmes linéaires. Ainsi le premier chapitre nous
permettra de définir l’automatique (position du problème).
Dans le deuxième chapitre nous définirons les systèmes linéaires et leurs modèles.
Ce premier chapitre nous permettra :
- De définir un système automatique.
- De définir les structures possibles d’un système automatique.
- De donner les moyens techniques de mise en œuvre.
- De préciser les étapes de la conception d’un système automatique.
II. Les systèmes automatiques

Aujourd’hui on a très souvent besoin que des grandeurs : physiques, chimiques, économiques,
subissent des variations précises (par exemple : on veut que la grandeur reste constante)
Exemples:
1) Niveau de liquide dans un réservoir.
Objectif : niveau de liquide constant en toute situation.
2) Vitesse d’un moteur.

Objectif : vitesse constante quelle que soit la charge.
3) Température d’un four.

Objectif : Température constante.

Automatiser ces systèmes revient à réaliser les différents objectifs sans intervention humaine.
Comment ?
Quelles sont les solutions ?
Elles reposent essentiellement sur deux structures :
- Boucle ouverte.
- Boucle fermée.
III. Les structures de commande

Les deux structures (ou schéma) possibles peuvent être appliquées l’une ou l’autre suivant les
conditions de fonctionnement du système.
1) Structure en boucle ouverte

Pour maintenir le niveau constant, il faut que à tout moment le liquide qui entre par 1 soit égale à
celui qui sort par 2 ; 1=2.
a) Définitions importantes
Dans le vocabulaire de l’automatique, le niveau de liquide est appelé grandeur à commander ou
bien SORTIE.
L’ouverture de la vanne 1 ou bien le débit en 1 s’appelle grandeur de commande ou ENTRÉE.
Entrée Sortie
Réservoir
(Débit) (Niveau)
b) Problème de la boucle ouverte

Cette structure peut très bien fonctionner lorsque rien ne dérange le système. Mais lorsque une
perturbation se présente (pluie, fuite), cette structure échoue.
2) Structure en boucle fermée
La structure en boucle fermée est envisagée pour régler en particulier le problème des
perturbations.

Son principe est basé sur le fonctionnement dirigé par un opérateur humain. Pour réaliser son
objectif par exemple de maintenir constant le niveau de liquide dans un réservoir, l’opérateur suit
les étapes suivantes :
- Mesurer le niveau (la sortie).
- Décider après comparaison (entre consigne et mesure).
- Agir sur (entrée).
Alors, il faut concevoir un système automatique qui réalise ces trois actions (étapes) à la place de
l’opérateur.
Donc, pour réaliser un système automatique, il faut trouver un moyen technique pour remplacer ces
trois opérations humaines.
Information
Action
Consigne Système
Décision à
commander
Information Mesure
Ainsi, on doit construire trois organes principaux ;

- Un organe de mesure qu’on appelle CAPTEUR.
- Un organe de décision qu’on appelle CONTROLEUR.
- Un organe d’action qu’on appelle ACTIONNEUR.
Il faut ensuite les réaliser. Ces organes et la liaison entre eux sont les moyens techniques de
réalisation.
IV. Les moyens techniques de réalisation

La première décision à prendre concerne le support d’information qui relie les différents organes.
On a trois supports :
a) Courant électrique faible : Dans ce cas les organes du système automatique sont des fils
électrique qui transportent un courant électrique faible, en générale compris entre 4mA et
20mA. Alors dans ce cas le capteur doit délivrer un courant électrique.
Exemple :
Thermocouple, Mesure de niveau.

b) L’organe de décision : le contrôleur est formé de composants électroniques qui peuvent être
de type continu (R, L, C, Amplificateur) ou numérique (calculateur).
Consigne Circuit électrique Courant
Courant
Capteur
Consigne
Calculateur
Capteur
c) La transmission de la décision à l’actionneur se fait par courant électrique (4mA à 20mA),

ce courant est alors amplifie pour pouvoir actionner les différents éléments (domaine de
l’électronique de puissance).
Carte de
Courant faible puissance
(4 à 20mA)
Actionneur
Autre possibilité :
Carte de
Actionneur
Courant faible puissance
hydraulique ou
(4 à 20mA) transformation pneumatique

Boucle fermée :
K
erreur entrée sortie
Consigne + Contrôleur Système
r(t) -
Mesure
V. Méthodologie de conception des systèmes de commande automatique

Quelles sont les méthodes, quelles sont les étapes qui permettent de concevoir (développer) un
système de commande automatique.
1) Les étapes de la conception
On a la grandeur à commander.
Exemples :
Vitesse d’un moteur, Température d’un four, Equilibrage des ailes d’un avion, Attitude d’un
satellite, Concentration d’un produit chimique.
Il faut alors :
 Décider quelle est la grandeur de commande.
 Boucle ouverte ou boucle fermée.
 Si boucle ouverte stop (rare)
 Si boucle fermée alors :
Le schéma de commande est par exemple :
e(t) u(t) y(t)

Consigne + Contrôleur Système
Analyse
r(t) -
Mesure
 Déterminer le modèle mathématique qui donne la relation entre l’entrée (grandeur de

commande) et la sortie (grandeur à commander).
 Les lois ou l’expérience.

 Calculer le contrôleur.
 Trouver la relation entre l’erreur et la grandeur de commande.
Synthèse  Réaliser le contrôleur.
 Faire les essais en simulation ou modèle réduit.
 Implantation sur site.
2) Les méthodes
a) Modélisation
Les méthodes de modélisation soient classées en deux classes :
- Celles qui utilisent les lois connues de la physique, de la chimie, de l’économie,
elles donnent des modèles avec des équations différentielles par rapport au
temps.
- Celles basées sur les essais directs sur le système. On obtient des tableaux ou des
courbes qui donnent le modèle.
b) Calcul des contrôleurs
Il existe plusieurs méthodes de calcul des contrôleurs qu’on peut les grouper en trois catégories :
 Méthodes classiques : travail en générale dans le domaine fréquentiel et temporel basées
sur la transformée de Laplace.
 Méthodes modernes : travail dans le domaine temporel.
 Méthodes post-modernes : combine les deux méthodes.

Chapitre II : Système et Modèle

I. Introduction
Un modèle mathématique dynamique décrit le comportement dans le temps d’un système. En
automatique un modèle dynamique est utilisé pour :
- La simulation.
- Le calcul des contrôleurs.
Cependant, on n’utilise pas le même modèle pour la simulation et le calcul du contrôleur. Pour la
simulation, on utilise un modèle non linéaire aussi précis que possible. Pour le calcul du contrôleur,
dans le domaine de l’automatique linéaire, on utilise un modèle linéaire généralement peu précis,
simple obtenu par simplification ou linéarisation du modèle non linéaire.
Le modèle linéaire simple nous permet d’utiliser des représentations spéciales. Ainsi au lieu
d’utiliser les équations différentielles, on utilise « la fonction de transfert » basée sur la transformée
de Laplace.
Essai 3
Système
Rare
1
2
Modèle de 4 Simplification Modèle linéaire

simulation simple
Linéarisation
Représentation par
Fonction de transfert
II. La représentation par fonction de transfert
Nous avons vu que les modèles mathématiques dynamiques sont représentés par des équations
différentielles ordinaires ou partielles. En particulier pour les modèles linéaires à coefficients
constants, la relation entrée-sortie peut s’écrire :
d n y (t ) d n1 y (t ) u
Système linéaire
y
an  an 1  ...  a0 y (t )  bu(t )
dt n dt n 1

Ces modèles linéaires sont utilisés pour le calcul des contrôleurs. Cependant, les équations
différentielles ne sont pas faciles, pour cela on les transforme en équations algébriques à l’aide de la
transformée de Laplace.
u y
Système linéaire G(P)=y/u
Le modèle algébrique obtenu s’appelle Fonction de transfert.
1) La transformée de Laplace : Définition et propriétés

a) Définition
Soit une fonction f(t), continue de R →R et soit une variable complexe p=σ+jω, la transformée de
Laplace (T.L) de f(t) est une fonction de P définie par :
( )= ( )
Fb(t) s’appelle T.L bilatérale.
Dans l’étude des signaux et systèmes, la variable "t" est le temps et on peut écrire f(t)=0, t<0.
( )= ( )
T.L unilatérale.
F(p) est fonction rationnelle.
+ +⋯+ +
( )=
+ +⋯+ +
b) Propriétés
i) Superposition
Si F1(p) es la T.L de f1(t) ; L(f1(t))= F1(p). et si F2(p) es la T.L de f2(t) ; L(f2(t))= F2(p).
Alors L(α1 f1(t)+α2 f2(t))= α1F1(p)+α2 F2(p).
ii) Décalage dans le temps
Soit f(t) et son décalage f(t-λ)
Calculons f(t-λ), t’=t-λ
( ( − )= ( ) ′

( ( − )=∫ ( ) ( )
′ f(t)
( )
( ( − )= ( ) ′
t
( ( − )= ( )
t
iii) Multiplication par une exponentielle λ
Soit f(t) et soit e-atf(t)
( ( )= ( )
( )= ( )
( ( )
( ( )= ( + )
iv) Dérivation
Soit f(t) et df/dt, L(f(t))=F(p)
= ( ) − (0) (0) = 0
= ( )
C’est pour cela que p est appelé Opérateur dérivation
= ( )− (0) − − −⋯−
Si
(0) = = =0∀
Alors :
= ( )

v) Intégration
( )
( ) = + (0)
( )= ( ) (0) = 0
Alors :
( )
( ) =
1/p : intégration.
vi) Convolution
On définit le produit de convolution par :
( )∗ ( )= ( − ) ( )
( )∗ ( )= ( ) ( )
vii) Terme de la valeur finale
On a :
( )
( ) = ( )=
( )
N(p), D(p) : polynômes.
Si les racines de D(p) (D+(p)=0) ont la partie réelle négatives
Alors
lim ( ) = lim ( )
→ →
viii) Terme de la valeur initiale
De la même façon :
lim ( ) = lim ( )
→ →

2) La transformée inverse de Laplace
La transformée de Laplace transforme une fonction f(t) de variable réelle "t" en une fonction de la
variable complexe "p" notée F(p).
F(p) est rationnelle pour f(t) continue.
( )
( ) = ( )=
( )
La transformée inverse nous donne f(t) lorsque on a F(p).
F(p)=L(f(t)) ; L : transformée de Laplace.
f(t) =L-1(F(p)) ; L : transformée inverse de Laplace.
La transformée inverse est donnée par :
1
( )= ( )
2
En générale, on n’utilise pas cette formule, on détermine les transformées inverses à partir des
tableaux, lorsque F(p) n’est pas donnée directement sur les tables, on transforme F(p) en une
combinaison d’éléments des tables, pour cela, soit :
( )
( )=
( )
( )= + + ⋯+ +
( )= + + ⋯+ +
i) Racine simples
Si D(p)=0 a des racines simples c.à.d.
( ) = ( − )( − ) … ( − )
≠ ≠ ≠⋯≠
Alors :
( )= + +⋯+
( − ) ( − ) ( − )
= [( − ) ( )]|
Exemple :
+3
( )=
( + 2)( + 5)
Calculer f(t) ?
( )= + +
( + 2) ( + 5)

= [( − ) ( )]|
+3 3
= ( ) =
( + 2)( + 5) 10
= [( − ) ( )]|
+3 −1
= ( + 2) =
( + 2)( + 5) 6
= [( − ) ( )]|
+3 −2
= ( + 5) =
( + 2)( + 5) 15
Alors :
3 1 1 1 2 1
( )= − −
10 6 ( + 2) 15 ( + 5)
D’où :
3 1 2
( )= − −
10 6 15
ii) Racines multiples
Si D(p)=0 a des racines multiples.
( )=( − )( − )…( − ) ….
( )= + +⋯+ + +⋯+
( − ) ( − ) ( − ) ( − ) ( − )
1
= (( − ) ( ))
!
= 0,1,2, … , − 1
1
= 0, = [( − ) ( )]
0!
Exemple :
+3
( )=
( + 1)( + 2)
Calculer f(t)=L-1(F(p))
( )= + + +
( + 1) ( + 2)

= [( + 1) ( )]|
+3 +3
= ( + 1) = ⇒ = −2
( + 1)( + 2) ( + 2)
= [( + 2) ( )]|
+3 +3 1
= ( + 2) = ⇒ =−
( + 1)( + 2) ( + 1) 4
1 +3 3
= (( ) ( )) = =
0! ( + 1)( + 2) 2
1 +3 3
= (( ) ( )) = =
1! ( + 1)( + 2) 2
−( + 6 + 7) 7
= ⇒ =−
( + 1)( + 2) 4
( )= + + +
( + 1) ( + 2)
−2 1 1 71 3 1
( )= − − +
( + 1) 4 ( + 2) 4 2
1 −8 1 7 6
( )= − − +
4 ( + 1) ( + 2)
D’où :
1
( ) = [−8 − −7+6 ]
4

ANNEXE :
 Solution des équations différentielles ordinaires linéaires :
Nous allons considérer la solution de l’équation suivante :
d n y (t ) d n1 y (t )
an  an 1  ...  a0 y (t )  bu(t )
dt n dt n 1
Les coefficients ai sont constants
La solution y(t) est la somme de deux termes :

y(t)=yl(t)+yf(t)
i) yl(t) est la solution générale de :
d n y(t ) d n 1 y(t )
an  a n 1  ...  a0 y(t )  0
dt n dt n 1
Equation sans second membre. Elle caractérise le comportement interne du système.
Elle dépende de la structure du système. Elle ne change que si on change la structure du système.
yl   A e i
it
τi: sont les solution de l’équation caractéristique :

an n  an 1 n 1  ...  a0  0
 i  a  jb
 Analyse de la solution :
yl (t )   A e i
it  A1e 1t  A2 e 2 t  ...  An e n t
yl(t) dépend des τi
yl i (t )  Ai e i t  Ai e ( a  jb)t  Ai e at e jbt
Nous avons plusieurs cas :
1- Si b=0
yi(t)
yl i (t )  Ai e at
at
Si ai<0 : yl i (t )  Ai e
Ai
t 0 yi (0)  Ai
t
t   yi (t )  0

yl i (t )  Ai e at
Si ai>0 :
yi(t)
t 0 yi (0)  Ai
t   yi (t )  
Ai
2- Si b≠0
yli (t )  Ai e (ai  jbi )t  Ai e ai t e jbi t
e jbi t : c’est un cosinus de bi t
( ai  jbi )t
Si ai<0 : yli (t )  Ai e  Ai e ai t e jbi t
t 0 yi (0)  Ai
t   yi (t )  0
( ai  jbi )t
Si ai>0 : yli (t )  Ai e  Ai e ai t e jbi t
t 0 yi (0)  Ai
t   yi (t )  
La solution de yli(t) dépend surtout de ai.
Si ai<0 alors yli(t) converge vers 0.
Si ai>0 alors yli(t) diverge vers ∞.
jbi t
Si ai=0 alors yli (t )  Ai e .
La réponse libre totale est la somme des yli(t). il suffit qu’un seul ai>0 pour que yl(t) diverge.

ii) yf(t) est la réponse forcée :
yf(t) dépend de u(t).
yf(t) est une solution particulière de l’équation avec second membre.
d n y (t ) d n 1 y (t )
an  an 1  ...  a0 y (t )  u (t )
dt n dt n 1
y(t)=yl(t)+yf(t)
yl(t) : dépend du système
yf(t) : dépend de l’entrée
Exemple :
u(t)=1,
Solution particulière : yf(t)=K.
Alors : y(t)=yl(t)+K.

III. La fonction de transfert :

Nous avons vu qu’un modèle mathématique est décrit par des équations différentielles, pour
faciliter le travail il est préférable d’avoir des équations algébriques.
La transformée de Laplace permet de transformer le modèle différentiel en modèle algébrique.
Ce modèle est particulièrement intéressant lorsque les conditions initiales sont nulles, on l’appelle
alors la FONCTION de TRANSFERT. Il s’obtient en utilisant la propriété de dérivation de la
transformée de Laplace.
 df (t ) 
L  pF ( p)  f (0)
 dt 
Si f(0)=0
 df (t ) 
L  pF ( p)
 dt 
Dans ce cas P peut être considéré comme l’opérateur (symbolique) dérivé.
Exemple :
d n y (t ) d n1 y (t )
an  an 1
 ...  a0 y (t )  bu(t )
dt n dt n 1
Conditions initiales nulles.
La fonction de transfert est le rapport de Y ( p) soit L(u (t ))  U ( p) , L( y (t ))  Y ( p)
U ( p)
an p nY ( p)  an1 p n1Y ( p)  ...  a1 pY ( p)  a0Y ( p)  bU ( p)

Y ( p) b
  G ( p)
U ( p ) an p n  an1 p n 1  ...  a1 p  a0
Définition :
La fonction de transfert donnée par Y ( p )  G ( p ) est définie pour les systèmes linéaires.
U ( p)
Pour les systèmes non linéaires :
dy(t )
 y 2 (t )  u(t )
dt
 dy(t ) 
L  y 2 (t )   Lu(t )
 dt 
 dy(t ) 
  L y (t )   Lu(t )
2
L
 dt 

Impossible de trouver la fonction de transfert, il faut simplifier le modèle non linéaire pour pouvoir
trouver la fonction de transfert.
Définition :
N ( p)
La fonction de transfert G(p) est généralement rationnelle G ( p) 
D( p )
D(p)=0 : équation caractéristique.
Les racines de D(p) s’appellent les pôles de la fonction de transfert.
Les racines de N(p) s’appellent les zéros de la fonction de transfert.
Les pôles de G(p) correspondent aux racines ( = 1, … , ) de l’équation caractéristique.
Alors si un système a une fonction de transfert G( p )  N ( p )

D( p )
Stable si (convergent) pi=ai+jbi , ai<0, réel(pi)<0 ∀
Instable si (divergent) pi=ai+jbi , ai>0, réel(pi)>0
Calculer si possible les fonctions de transfert des systèmes :
2
d y (t ) dy(t )
1) 4
2
  y (t )  u (t ) ( y0'  0, y0''  0, y0  0)
dt dt
dy(t ) du(t )
2)  14 y(t )  u (t )  3 ( y0'  0, y0  0, u0  0, u0'  0)
dt dt
dy(t )
3)  8 y(t )  u (t  4) ( y0'  0, y0  0, u0  0)
dt
2 2
 d y(t )  '
4) 
2
  y (t )  u (t ) ( y0  0, y0  0)
 dt 
d 2 y(t )
5) 2
 2 y(t )  u (t ) ( y0'  0, y0  1)
dt
d 2 y (t )
6) 2
 y (t )  u(t ) ( y0'  0, y0  0)
dt
dy(t ) du(t )
7)  3 y (t )  u (t )  2 ( y0'  0, y0  0, u0  0)
dt dt

Définition :
Une fonction de transfert G(p) est propre si deg(Num(G(p))< deg(Dén(G(p))
2
Exemple : G ( p )  p  2 p  1 (propre)
3
p 1
Une fonction de transfert G(p) est strictement propre si deg(Num(G(p))<= deg(Dén(G(p))
Les systèmes réels sont propres ou strictement propres.
Les Schémas Blocs :

Généralement un système est composé de plusieurs « sous-systèmes » (ou petits systèmes) reliés
chaque sous systèmes à sa fonction de transfert (entrée sortie).
Un schéma bloc est une représentation du système qui décrit chaque sous-systèmes par sa fonction
de transfert sous forme d’un rectangle G(p) et décrit les connections entre les sous-
systèmes.
Exemple : G1(p) G2(p)
G3(p)
Définition cascade
U1 Y1 U2 Y2 Un Yn
G1(p) G2(p) Gn(p)
U1 Y1 U2 Y2
G1(p) G2(p)
Y2 ( p) Y2 ( p) U 2 ( p) Y2 ( p ) Y1 ( p )
G( p)   .  .  G1 ( p ).G2 ( p)
U 1 ( p ) U 2 ( p) U 1 ( p ) U 2 ( p ) U 1 ( p )
Généralement :
G ( p)  G1 ( p).G2 ( p )...Gn ( p)
où G1(p) G2(p) … Gn(p) sont en cascade.

Chapitre III
Les Systèmes en boucle fermée
Régulateurs et performances
I / Introduction :
Un système de commande en boucle fermée est composé de plusieurs éléments reliés les uns aux
autres comme suit :
Actionneur
Consigne Comparateur Système
+ à
Contrôleur commander
Capteur
Sous forme de schéma bloc :
Système
u’(t) y(t)
Consigne + e(t) u(t)
Contrôleur Actionneur à
r(t) - commander
Capteur
Dans ce cours (et d’une façon générale pour les systèmes linéaires) on considère que :
Capteur
Actionneur

Le schéma standard d’un système en boucle fermée devient alors :
Système
u’(t) y(t)
Consigne + e(t) u(t)
Contrôleur Actionneur à
r(t) - commander
Système à retour unitaire

Définition :
Régulation et Asservissement : lorsque la consigne r(t) est une constante r(t)=R on dit qu’on a un
problème de régulation. La régulation est très courante dans l’industrie des procédés (PROCESS
CONTROL). Dans cette industrie on désire maintenir constantes des grandeurs telle que :
Température, Pression, niveau d’un réservoir, le débit, Concentration.
Le but principal : est la lutte contre les perturbations.
- Lorsque r (t) est variable on a un problème de poursuite ou d’asservissement, on rencontre
ces problèmes en Robotique, Aéronautique, Aérospatiale. Exemple : trajectoire d’un robot.
La régulation est un cas particulier de l’asservissement quand r(t)=constante. A ce niveau nous
pouvons donc dire que le contrôleur devient un régulateur lorsque la consigne est une constante.
- Deux questions peuvent alors être posées et dont les réponses sont l’objet de ce chapitre :
 Qu’est ce qu’un régulateur (ou contrôleur).
 Quel est l’objectif du régulateur dans la boucle.
Autrement dit quelles sont les performances désirées.
Il y a plusieurs types de régulateurs (contrôleurs). Dans ce chapitre nous allons s’intéresser à ce
qu’on appelle les régulateurs classiques parce qu’ils sont les plus anciens et les plus utilisés dans
l’industrie.
II/ les régulateurs classiques :

Dans le schéma bloc ci-dessus l’entrée du régulateur est l’erreur e(t) et sa sortie est le signal de
commande u(t). Ainsi le contrôleur réalise une fonction de la forme :
u (t )  F (e(t ))
e(t) Contrôleur u(t)

Ou bien sous forme de fonction de transfert, si F est linéaire U ( p )  G ( p ) , R : Rationnelle.

R
E ( p)
Il y a plusieurs fonctions F selon le type du régulateur. Les régulateurs classiques sont formés de
combinaison de fonction simple au nombre de trois appelés actions élémentaires. Ces actions
sont :
i. Action proportionnelle P : cette action réalise la fonction :
u(t )  ke(t ) G P
( p)  k 
Cette action améliore la sensibilité.
ii. Action intégrale I : cette action réalise la fonction :
1t  1 
u (t )   e(t )dt  GI ( p)  
i 0   i p 
Cette action élimine complètement l’erreur entre la consigne et la sortie. : constante de temps
d’intégration.
iii. Action dérivée D : cette action réalise la fonction :
de(t )
u (t )   d G (p)  p
D d
dt
C’est action anticipatrice, elle améliore la stabilité.
Les régulateurs classiques sont des combinaisons de ces actions. On distingue :
i. Le régulateur P :
U ( p)
GP ( p )   k , k : gain
E ( p)
Exemple :
Commande de niveau d’un réservoir par régulateur P.
Le modèle du réservoir :
dh qi Contrôleur h(t)
 qi  r h ,
dt
h : le niveau, qi : le débit, hc : la consigne.

Système h(t)
e(t) qi(t)
Consigne + k à
hc(t) - commander
e=hc-h
qi=k.e
Simulation :
Initialisation : h0=1m, hc=2m
Loop :
( + ) = ( ) + (∆ )( ( ) − ∗ ( ( ))
( )= − ( + )
( )= . ( )
end (Loop)
ii. Le régulateur PI :
U ( p) 1 ou 1t
GPI ( p)  k u (t )  ke(t )   e(t )dt
E ( p) i p i 0
Simulation :
Initialisation : h0=1m, hc=2m, somme=0
Loop :
( + ) = ( ) + (∆ )( ( ) − ∗ ( ( ))
( )= − ( + )
= + ( ). ∆
( )= . ( )+
end (Loop)
iii. Le régulateur PID :

U ( p) 1 1t de(t )
GPID ( p)  k   d p ou u (t )  ke(t )   e(t )dt   d
E ( p) i p i 0 dt

Nous remarquons que chaque régulateur est caractérisé par des paramètres :
P : k Gain
PI : k et
PID : k , et
Comment obtenir les valeurs de ces paramètres ? c’est le problème de la conception. Pour cela il
faut :
- Spécifier les performances.
- Développer une méthodologie.
Dans la deuxième partie de ce chapitre nous allons développer et analyser les performances des
systèmes automatiques.
III/ Les performances des systèmes en boucle fermée :

Le calcul des contrôleurs est basé sur des objectifs que doit réaliser le système en boucle fermée.
La conception des systèmes automatique passe alors par les étapes suivantes :
1) A partir des cahiers de charges établir les spécifications des performances.
2) Choisir un ou des régulateurs.
3) Choisir une méthode de calcul des régulateurs.
4) Calculer le régulateur.
5) Si le système en boucle fermée est satisfaisant stop, sinon aller à 3 ou 4.
Nous allons s’intéresser à l’étape 1 ; quelles sont les performances à satisfaire ?
Par sa nature même, l’objectif évident d’un système automatique en boucle fermée est la précision
mesurée par la valeur de l’erreur entre la consigne et la valeur réelle.
Par sa construction la boucle fermée, introduit un problème important. Le problème de la stabilité.
Définissons d’abord la stabilité.
Définition de la stabilité :
Un système est dit stable si a une entrée bornée il répond par une sortie bornée.
Stabilité BIBO (Borned Input Borned Output).
Un système est stable si : écarté de sa position d’équilibre il tend à y revenir.
Entrée bornée : u(t)
∀ ( )<∞
Exemple : u(t)=1 ∀ borné
u(t)=t ∀ non borné
le deuxième objectif des systèmes en boucle fermée est donc d’assurer une bonne stabilité.

1°/ analyse de la stabilité :

Considérons un système linéaire décrit par une équation différentielle :
d n y (t ) d n1 y (t )
an  an 1
 ...  a0 y(t )  bu(t )
dt n dt n1
Nous allons analyser la solution de cette équation
y(t)=yl(t)+yf(t) : yl(t) dépend du système (des racines de l’équation caractéristique).
yf(t) dépend de u(t).
la conclusion de notre étude étant :
yl(t) 0 si les racines de l’équation caractéristique ont des parties réelles négatives.
yl(t) ∞ si au moins une racine a la partie réelle positive.
yl(t) : sinus si une racine est imaginaire et les autres parties réelles négatives.
Conclusion :
Puisque les racines de l’équation caractéristique sont les pôles de la fonction de transfert, on peut
alors énoncer propriété.
Un système est stable si et seulement si les pôles de sa fonction de transfert sont à partie réelle
négative.
Exemple :
d 2 y(t ) dy(t )
2 2
5  y(t )  2u(t )
dt dt
Quelles sont les racines de l’équation caractéristique ? le système est il stable ?
Ecrire la fonction de transfert et donner les pôles.
L’équation caractéristique : 2r2+5r+1=0
 5  25  8
r12 
4
r12 <0 : système stable.
La fonction de transfert :
2 Les pôles : p12   5  25  8
G( p)  2
2 p  5p 1 4
Conséquence de la propriété :
- La stabilité est une caractéristique interne du système.
- Les zéros des fonctions de transfert n’influent pas sur la stabilité. G ( p)  N ( p)

D( p )
Les zéros sont les racines de N(p).

On peut cependant remarquer que si on a des zéros a partie réelle positive on a un comportement
inverse (système a non minimum de phase).
p 1
G( p) 
( p  1)( p  2)
a/ Notion de stabilité relative :

Lorsque un système est stable, on veut connaitre sont degré de stabilité, c.à.d. estimer sa distance
par rapport à l’instabilité.
On peut donc mesurer le degré de stabilité à l’aide des paramètres de la réponse y(t).
Ces paramètres sont en générale :
- Le premier dépassement Md (premier pic (first over shoot).
- Le temps de réponse à : x%.
- Le temps de montée à : s%.

Ces paramètres définissent le régime transitoire qui commence à t=0 et se termine à T=Tr. Le
régime transitoire caractérise le degré de stabilité.
b/ analyse des système du 2ème ordre et du 1er ordre :

Nous allons maintenant analyser le régime transitoire pour des systèmes simple du 1er et du 2ème
ordre.
1) système du 2ème ordre
n2
G ( p) 
p 2  2n p   n
Forme standard du 2ème ordre.
ωn : pulsation naturelle non amortie.
ξ : facteur d’amortissement.
Si u(t)=1 :
G ( p) 1 n2
Y ( p)  
p p p 2  2 n p   n
1
1 n2 
y(t )  L  2 
 p p  2n p  n 
On a trois cas :
- les pôles de G(p) sont réels.
- les pôles de G(p) sont doubles.
- les pôles de G(p) sont complexes.
Les pôles de G(p) sont p12  n  n  2  1
 ξ>1 : p12  n  n  2  1
p12 : réels. Solution amortie :

n  1 t 1 t 
y (t )  1  e  e 
2  2  1   

 ξ<1 : p12  n  jn 1   2
p12 : complexes. Solution pseudopériodique (sous-amortie) :

y (t )  1  e t sin d t   
  n , d  n 1   2 ,   cos1 
 ξ=1 : p12  n
p12 : double réels. Solution critique:

y(t )  1  e   t 1  nt 
n
la reponse indicielle
1
0.9
0.8
0.7
0.6
0.5
0.4
0.3
0.2
0.1
0 500 1000 1500 2000 2500 3000 3500 4000 4500 5000
time

Caractéristique du système sous amortie :
4
- le temps de réponse : Tr2% 

1.8
- le temps de montée : Tm 
n


er 1 2
- le 1 dépassement : Md  e
2) système du 1er ordre

1 forme standard, τ : constante de temps.
G( p ) 
1  p
(un système sous une forme non standard G( p)  1 pour arriver au système on attire a)
pa
u(t) G( p) 
1 y(t)
1  p
Réponse à l’échelon u(t)=1

1 1 
y (t )  L1   
 p 1  p 

A B 
y (t )  L1   
 p 1  p 
1  
y (t )  L1   
 p 1  p 
t

y (t )  1  e 
Plus τ augmente plus l’amortissement augmente.

3°/ carte des pôles et stabilité relative :

La carte des pôles définit la place des pôles d’un système dans le plan complexe.
p3
G( p ) 
( p  1)( p  15)( p  (2  j ))( p  (2  j ))
- Plus on approche de l’axe Imp plus la stabilité « diminue » (les oscillations augmentent).
- Plus on s’éloigne de l’axe Imp plus la stabilité « augmente » (les oscillations diminuent).
Les pôles qui sont près de l’axe Imp s’appellent pôles dominants.
4°/ Critère de Routh :

Pour savoir si un système est stable il faut calculer ses pôles et vérifier s’ils ont tous des parties
réelles négatives. Le critère de Routh est une méthode simple qui permet de déterminer le nombre
de racines à partie réelles positive d’un polynôme.
Soit un polynôme D( p)  an p n  a n1 p n1  ...  a1 p  a0
On construit le tableau de Routh :

an a n 2
a n 4 a n 6
a n 1 a n 3
a  a  a n  an 3 a  a  a n  a n 5 a n 5 a n 7
b1  n 1 n 2 b2  n 1 n 4 an 1  a n 6  an  a n 7
an 1 a n 1 b3  ...
b1  a n 3  a n 1  b2 b1  a n 5  a n 1  b3 an 1
... ...
b1 b1
Le nombre de racines a partie réelle positive de D(p) est égale au nombre de changements de signe
dans la 1ère colonne du tableau.
Exemple :
p3
G( p) 
2 p  3p  p3  2 p  1
5 4
Ce système est stable ?

Appliquant le critère de Routh au dénominateur de G(p).
2 1 2
3 0 1
3 62
0
3 3
4 1
 0
1 1
19
 0
12
1
Deux pôles a partie réelles positive donc système instable.
Cas particulier :
Si on trouve un zéro sur la première colonne (marginalement stable) : une racine imaginaire j.a
Il faut lever l’indétermination (car le zéro deviendra pivot ; division par zéro). On forme le
polynôme auxiliaire avec la ligne précédente.

pn . . . . . . . .
p n 1 . . . . . . . .
p n 2 . . . . . . . .
. . . . . . . . .
. . . . . . . . .
p n i c1 c2 . . . . . .
0 . . . . . . .
. . . . . . . .
. . . . . . . .
. . . . . . . .
Le polynôme auxiliaire s’obtient comme suit :

1) On prend les coefficients de la ligne [c1 , c2 , …]
2) On forme le polynôme de degré pn-i.
3) on dérive ce polynôme.
4) On remplace la ligne du zéro par les éléments du polynôme dérivé.
II / Etude de la précision :
Après avoir étudié la stabilité des systèmes, nous considérons maintenant la précision.
La précision caractérise le régime permanent (ou établi) donc pour un temps assez élevé.
La précision est alors déterminée par l’erreur en régime permanent donnée par :
  lim e(t ) d’où e(t )  r (t )  y(t )
t 
avec : r(t) la consigne du système.

y(t) la sortie du système.
Pour l’étude de la précision on considère un système en boucle fermée à retour unitaire.
r(t) + e(t) y(t)

kG(p)
-
Pour calculer   lim e(t ) nous allons appliquer le théorème de la valeur finale :
t 
lim pE ( p )  lim e(t ) avec E ( p)  Le(t ) 

p 0 t 

Calculons E(p) :
Soit Y ( p)  L y (t )  et R( p)  Lr (t ) 
1
E ( p )  R( p )
1  kG( p)
Remarque :
Contrairement à la stabilité qui ne dépend pas de l’entrée la précision dépend de l’entrée. Pour cela,
nous allons choisir des consignes (entrée) standards :
- Echelon r(t)=1, t>=0.
- Rampe r(t)=a.t, t>=0.
- Parabole r(t)= a.t² ,t>=0.
Nous avons alors :
Lorsque r(t) échelon : εp : erreur de position.
r(t) rampe : εp : erreur de vitesse.
r(t) parabole : εp : erreur en régime d’accélération.
1
r (t )  1  R( p) 
p
1 1 1 1
 p  lim pE( p)  lim pR( p)  lim p  lim
p 0 p 0 1  kG( p) p0 p 1  kG( p) p0 1  kG( p)
1
r (t )  t  R( p) 
p2
1 1 1 1 1
 v  lim pE( p)  lim pR( p)  lim p 2  lim  lim
p 0 p 0 1  kG( p) p0 p 1  kG( p) p0 p  pkG( p) p0 pkG( p)
2
r (t )  t 2  R( p) 
p3
1 2 1 2 1 2
 a  lim pE( p)  lim pR( p)  lim p 3  lim 2  lim 2
p 0 p 0 1  kG( p) p0 p 1  kG( p) p0 p 1  kG( p) p0 p kG( p)

r(t)=1 1 1 k p  lim kG ( p )
 p  lim p  p0
p 0 1  kG( p) 1 k p
r(t)=t 1 1 k v  lim pkG ( p)
 v  lim v  p0
p 0 pkG( p ) kv
r(t)=t² 2 1 p 2 kG( p)
 a  lim 2 a  k a  lim
p  0 p kG ( p ) ka p 0 2
Définition :
Le type d’une fonction de transfert G(p) est égale au nombre de pôles de G(p) au point p=0.
Une fonction de transfert a :
 p  4 
Un zéros au point  p  2 Un zéros au point  p  1 Un zéros au point 
p0 
des pôles aux points  p  3  des pôles aux points  p  0  des pôles aux points  p  0 

     3  2 j 
 p  2   p3   
 type 0   type1   3  2 j 
 type 2 
p2
Type 0 : G( p)  type 0 : position
( p  3)( p  2)
p 1
Type 1 : G ( p)  type 1 : vitesse
p ( p  3)
p4
Type 2 : G ( p)  2
type 0 : accélération.
p ( p  3  2 j )( p  3  2 j )
Nous allons calculer εp, εv, εa pour les systèmes de type 0,1,2.
 Système de type 0 :
k
kG ( p )  , pi  0 i .
( p  p1 )( p  p 2 )...
k
kG ( p)  , pi  0 i .
p ( p  p1 )( p  p 2 )...
k
kG ( p)  2
, pi  0 i .
p ( p  p1 )( p  p 2 )...
Pour un système de type 0 on a :

1 k k
p  , k p  lim 
1 k p p 0 ( p  p )( p  p )....
1 2 p1 . p 2 . p 3 ...
Pour un système de type 1 on a :
k k 1 1
k p  lim   , p   0
p 0 p( p  p1 )( p  p 2 ).... 0 1 k p 
k
k v  lim pkG( p)  lim p 0
p 0 p 0 ( p  p1 )( p  p 2 )....
1 1
v   
kv 0
Type 0 1 2
Erreur
εp kp : constante kp   kp  
1 p 0 p 0
p 
1 k p
εv kv =0 kv : constante kv  
v   1 v  0
v 
kv
εa ka=0 ka=0 ka=constante
a   a   1
a 
ka
Rejet des perturbations :
Considérons le système suivant avec d(t) une perturbation à l’entrée :
d(t)
Gd(p)
Système
u(t) + e(t) y(t)
G(p)
La sortie y(t) est alors donnée par :
Y ( p )  G( p )U ( p )  Gd ( p ).D( p ) 
L(d (t ))  D( p), L(u (t ))  U ( p), L( y(t ))  Y ( p)
Y ( p )  G( p )U ( p)  Gd ( p).D( p )G( p) , avec Gd ( p ).D( p)G ( p) c’est l’effet de perturbation.

Perturbation en boucle fermée :

d(t)
Gd(p)
y(t)
r(t) + e(t) +
Gc(p) G(p)
-
Gc(p) est un contrôleur.
Cherchons la fonction de transfert en boucle fermée, entrée r(t) et d(t) et la sortie y(t).
Y ( p )  G( p )U ( p)  Gd ( p).D( p )G( p)
Y ( p )  G( p )Gc ( p )R( p )  Y ( p)  Gd ( p ).D( p)G( p ) , avec U ( p )  Gc ( p )R( p)  Y ( p)
R( p )G ( p )Gc ( p ) D ( p )G ( p )G d ( p)
Y ( p)  
1  Gc ( p).G ( p ) 1  Gc ( p ).G ( p)
Pour diminuer l’effet de la perturbation, on doit augmenter Gc(p).G(p) : gain de boucle (Loop
Gain).
R( p )G ( p )Gc ( p ) D ( p )G ( p )G d ( p)
Y ( p)  
1  Gc ( p).G ( p ) 1  Gc ( p ).G ( p)
Effet sur la consigne :
R( p )G ( p)Gc ( p)
Y0 ( p ) 
1  Gc ( p ).G ( p)
Si Gc(p).G(p)>>1 alors :
R( p )G ( p )Gc ( p )
Y0 ( p)   Y0 ( p )  R( p)
Gc ( p ).G ( p)

Perturbation sur la sortie : d(t)
Gd(p)
u(t) y1(t) + y(t)

G(p)
d(t)
y(t) : sortie perturbée.
Gd(p)
r(t) + e(t) + y(t)

Gc(p) G(p)
-
R ( p )G ( p )Gc ( p)
Y ( p)   D ( p ).Gd ( p ) , perturbation sans intérêt.
1  Gc ( p ).G ( p)
La perturbation n’étant pas mesurée (capturée) : elle est à l’extérieur de la boucle de commande.
d(t)
Gd(p)
e(t) + y(t)
r(t) +
Gc(p) G(p)
-
R( p)G( p )Gc ( p) D( p).Gd ( p)

Y ( p)  
1  Gc ( p).G ( p) 1  Gc ( p ).G( p)

Chapitre IV
La méthode du lieu des pôles
I. Introduction
Dans ce chapitre nous allons présenter une méthode pour l’analyse et la synthèse des systèmes en
boucles fermée. Cette méthode travaille dans le domaine temporel, cela veut dire que les
performances désirées s’expriment dans le temps, c.à.d. temps de réponse (Tr), temps de monté
(Tm) et le dépassement (Md).
Soit une fonction de transfert :
( − )( − )…( − )
( )=
( − )( − )…( − )
Donc m zéros et n pôles.
Lorsque tous les pôles sont dans la partie gauche du plan complexe, le système est stable.
Im p
Stable
Réel p
Considérons un système en boucle fermée :
r(t) + e(t) y(t)

kG(p)
-
Y ( p) KG ( p )
  H ( p)
R( p ) 1  KG ( p )
Les zéros de H(p) sont les zéros de G(p)

Les pôles de H(p) sont les racines de 1+KG(p)=0

Ces pôles déterminent la stabilité du système en boucle fermée.

La méthode du lieu des pôles développé par EVANS est basé sur le tracé du déplacement des pôles
de H(p) lorsque k varie de 0 à ∞
Exemple :
r(t) + e(t) y(t)

k/(p+1)(p+4)
-
Le système en boucle ouverte est il stable?

Oui, car les pôles en BO sont (-1, -4), ils sont à partie réelle négatifs.
Le système en boucle fermée est il stable ?
Il faut calculer la fonction de transfert en boucle fermée.
K
Y ( p) ( p  1)( p  4) K K
H ( p)    
R( p) 1  K K  ( p  1)( p  4) p 2  5 p  ( K  4)
( p  1)( p  4)
La méthode d’EVANS consiste à tracer le lieu des racines de 1+KG(p)=0.
Ici : p2+5p=(K+4)=0
Les racines sont :
 5  9  4K  5  9  4K
p1  , p2 
2 2
Ces pôles peuvent être réels ou complexe dépendant de K.

9
Ils sont réels si : K ,
4
9
Ils sont complexe si K  ,
4
1) Pour K=0 ; p1=-1, p2=-4, ce sont les pôles en boucle ouverte.
9
2) Pour 0 K  ,alors 2.5  p1  pôle
1, réel.
4
9  5  j 4 K pôle
9 complexe à partie réelle constante.
3) Pour K  , alors p1  ,
4 2
9
4) Pour 0  K  , alors 4  p2  2pôle
.5, réel.
4
9  5  j 4K  9
5) Pour K  , alors p2  pôle, complexe à partie réelle constante
4 2
Le système en boucle fermée reste stable ∀ >0
Im p
K
Branches
p1
P2 p1
-4 -1 Réel p
K=0 K=0
P2
K
Cette méthode directe n’est pas applicable pour des degrés supérieurs à 2. Pour cela EVANS a
développé une méthode pour tracer le lieu des pôles.
Nous allons présenter cette méthode dans la section suivante, ensuite nous verrons comment
l’utiliser pour calculer les contrôleurs.
II. La méthode d’EVANS

Pour tracer le lieu des pôles considérons une équation d’une variable complexe : p=a+jb.
D(p)=0=A+jB
On peut écrire cette équation comme suit :
A=0 et B=0 ; | | = 0, (∠ = 0)
D(p)=-1⤇| | = 1, (∠ = (1 + 2 ))
La méthode d’EVANS est basée sur la solution géométrique : phase et module de 1+KG(p) qui peut
aussi s’écrire : KG(p)=-1
Un point ∈ ℂ est solution de l’équation caractéristique si :
| ( )| = 1 et ∠ ( ) = (1 + 2 )
Exemple :
+1
( )=
+2
+1
1+ ( )= 1+
+2
+1
= −1
+2
+1 +1
= −1 ∠ =−
+2 +2
+1
∠ = − ⤇ ∠( + 1 ) − ∠( + 2 ) = − ⤇ − =−
+2
On cherche tous les points p qui satisfont cette relation angulaire, avec p=a+jb ; p quelconque.
Ce point p n’est
pas solution P P+1 P+2 Im p
ϕ1-ϕ2
ϕ2
-2 -1 ϕ1 +1 +2 Réel p
∠( − (−1)) − ∠( − (−2)) = −
La figure montre que p=a+jb ; pour a≠0 et b≠0
P n’est pas solution 1+KG(p)=0, cependant si b=0 et a≠0, on a la figure suivante :
Im p
P=-0.5
-2 -1 (P-(-1)) (P-(-2)) Réel p
ϕ1-ϕ2=0
∠( − (−1)) − ∠( − (−2)) = −
 Si p>-1
La solution n’est pas satisfaite.

 Si -2<p<-1 Im p
P=-1.5, p+1=-0.5, p+2=-0.5
K=∞
K=0 (P-(-1))(0) π
Réel p
-2 P -1
(P-(-2))(0)
Le lieu des pôles
 Si p<-2 alors :
∠( − (−1)) − ∠( − (−2)) = 0
p n’est pas solution de 1+KG(p)=0
Pour K=0 le pôle de H(p) est p=-2 (c’est le pôle en boucle ouverte)
+1
= −1 ⤇ ( + 1) = − − 2
+2
=0⤇ = −2
→∞⤇ = −1 c’est le zéro en boucle ouverte
Nous avons tracé le lieu des pôles pour un système simple par tâtonnement pour un système
complexe.
EVANS a développé des règles simples pour construire le lieu des pôles de H(p).
II.1. Les règles d’EVANS

Soit le système en boucle fermée (BF) avec :
( − )( − )…( − )
( )=
( − )( − )…( − )
zi, i=1,…,m ; les zéros de G(p)
pi, i=1,…,n ; les pôles de G(p)
pour tracer le lieu des pôles de H(p), on utilise les règles suivantes :
KG ( p)
H ( p) 
1  KG ( p )

Règle N°01 :
Tracer le plan complexe et placer : les pôles de G(p) :x
Les zéros de G(p) :o
Im p
K=0
K=∞ K=0 Réel p
K=0
Puisque le nombre de pôles de H(p)= nombre de pôles de G(p)=n ; le lieu de pôles de H(p) a "n"
branches, une branche commence à K=0 à un pôle de G(p) et elle se termine à un zéros de G(p) à
K→∞.
Lorsque le nombre de zéros est insuffisant, on place les zéros qui manquent à l’infini.
r(t) + e(t) y(t)

K/(p(p+2)(p+5))
-
Les pôles de G(p) :
0, -2, -5
Les zéros de G(p) :
Im p ∞, ∞, ∞
Il existe : 03 branches
OA=-7/3
60°
K=0 K=0 K=0 Réel p

-5 A -2 0
-60°
P2=-0.87

Règle N°02 :
Sur l’axe des réels, le lieu des pôles existe à gauche d’un nombre impair de pôles et de zéros.
Règle N°03 :
Point de sortie des branches de l’axe des réels est la solution de :
d  1 
 0
dp  G ( p ) 
d
 p( p  2)( p  5)
dp
 7  19 p1  3.8 p2  0.87
p1,2 
3
Règle N°04 :
Angle de sortie des branches de l’axe des réels :
90°
Pour deux branches : 90°
Pour 04 branches : 60°
60°
Règle N°05 :
Asymptotes :
Point d’appui sur l’axe des réels :
OA 
 pôles de G( p)  zéros de G( p)
nombre de pôles de G( p)  nombre de zéros de G ( p) 
OA 
2  5  0  0  7
3  0 3
L’angle de l’asymptote avec l’axe des réels :
k 
1  2k 
k  0,  1,  2
nombre de pôles de G( p)  nombre de zéros de G( p)
 
3 zéros à l’infini ⤇3 asymptotes. 0  , 1   ,  1 
3 3
Règle N°06 :
Point de coupure de l’axe imaginaire des racines de 1+KG(p)=0 sont : j :
p   alors
1+KG(jω)=0
On calcul K à l’aide du critère de Routh marginalement stable.
Règle N°07 :
Angle de départ d’un pôle complexe de G(p)
 ( p0 )  180    p
i0
0  pi    p
j
0 zj 
Angle d’arrivée à zéro de G(p)
 ( z0 )   z
i

0  pi  z
j 0

0  z j  180
Exemple 01:
r(t) + e(t) y(t)

K/((p+1)(p+2))
-
Y ( p)
H ( p) 
R( p)
Tracer le lieu des pôles de H(p)
Lieu de pôle Im p
Asymptote
K=0 K=0 Réel p

-2 -1
3
2

Règle N°01 :
On trace le plan complexe et on place : les pôles de G(p) :x
Puisque le nombre de pôles de H(p)= nombre de pôles de G(p)=2 ; le lieu de pôles de H(p) a "2"
branches, une branche commence à K=0 à un pôle de G(p) et elle se termine à un zéros de G(p) à
K→∞.
Lorsque le nombre de zéros est insuffisant, on place les zéros qui manquent à l’infini (2 zéros à ∞).
r(t) + e(t) y(t)

K/((p+1)(p+2))
-
-1, -2
∞, ∞
Règle N°02 :
Règle N°03 :
d  1 
 0
dp  G ( p ) 
d
( p  1)( p  2)  0 (2 p  3)  0  p  3
dp 2
Règle N°04 :
On a deux branches : 90°
Règle N°05 :
Asymptotes :
OA 
 pôles de G( p)  zéros de G( p) OA 
1  2  0  3
nombre de pôles de G( p)  nombre de zéros de G ( p)  2  0 2

k 
1  2k 
k  0,  1,  2
 
2 zéros à l’infini ⤇2 asymptotes. 0  ,  1 
2 2
Règle N°06 :
p   alors
1+KG(jω)=0
On calcul K à l’aide du critère de Routh marginalement stable.
Pas de coupure de l’axe imaginaire des racines de 1+KG(p)=0
Règle N°07 :
 ( p0 )  180    p
i0
0  pi    p
j
0 zj 
 ( z0 )   z
i

0  pi  z
j 0

0  z j  180
Pas de pôles complexes, donc pas d’angle de départ d’un pôle complexe ni angle d’arrivée à un
zéros de G(p).
Exemple 02:
r(t) + e(t) y(t)

K/(p(p+4)2+16))
-
Y ( p)
H ( p) 
R( p)
Tracer le lieu des pôles de H(p)
Les pôles de G(p) : 0, -4+4j, -4-4j
Les zéros de G(p) : ∞, ∞, ∞

Im p
jω0
K=0
-45°
K=0 Réel p
8 0
3
45°
K=0 -jω0
Règle N°01 :
On trace le plan complexe et on place : les pôles de G(p) :x
Puisque le nombre de pôles de H(p)= nombre de pôles de G(p)=3 ; le lieu de pôles de H(p) a "3"
Lorsque le nombre de zéros est insuffisant, on place les zéros qui manquent à l’infini (3 zéros à ∞).
r(t) + e(t) y(t)

K/( (p(p+4)2+16))
-
0, -4+4j, -4-4j
∞, ∞, ∞
Règle N°02 :
Règle N°03 :
d  1 
 0
dp  G ( p ) 
Un pôles à 0 et deux pôles complexe conjugués, donc pas d’intersection avec l’axe des réels.

Règle N°04 :
Pas de branches sortants de l’axe des réels.
Règle N°05 :
Asymptotes :
OA 
 pôles de G( p)  zéros de G( p)
nombre de pôles de G( p)  nombre de zéros de G ( p) 
OA 
4  4 j  4  4 j  0  0  8
3  0 2
k 
1  2k 
k  0,  1,
 
3 zéros à l’infini ⤇3 asymptotes. 0  ,  1  , 1  
3 3
Règle N°06 :
p   alors
1+KG(jω)=0
p[(p+4)2+16]+K= jω[(jω+4)2+16]+K=0
jω[(-ω2+8jω+16)+16]+K=- jω3-8ω2+jω32+K=0
  j 3  j 32  0    32

 8 2  K  0  K  256
Règle N°07 :
 ( p0 )  180    p
i0
0  pi   p
j
0 zj 
 ( z0 )   z
i

0  pi  z
j 0

0  z j  180

G ( p)   
( − )( − )…( − )
( )=
( − )( − )…( − )
( p  z1 )   ( p  z 2 )   ...  ( p  z m )   ( p  p1 )   ( p  p2 )   ...  ( p  p n )    (1)
L’angle de départ d’un pôle p1 est :
( p  p1 )     ( p  z1 )   ( p  z 2 )   ...  ( p  z m )   ( p  p2 )   ...  ( p  pn ) 

Puisque p1 appartient au lieu des pôles on a :
 p1  lim ( p  p1 ) 
p  p1
 n m 
 
 p1  lim ( p  p1 )  lim   
p  p1 p  p1 
i 1

( p  pi )  
j 1

 (p z j) 

 
p1=-4+4j
n m
 p1     ( p  p ) ( p  z )
i 1
i
j 1
j
pas de zéros alors

n
 p1     ( p  p )
i 1
i
 p1    (4  j 4  4  j 4)   (4  j 4)   p1    ( p1  (4  j 4))   ( p1 _ 0) 
 p1    ( j8)   (4  j 4)   p1  180  90  135   p1  45
III. Calcul des contrôleurs par la méthode du lieu des pôles

Nous allons voir dans cette section comment utiliser le lieu des pôles pour calculer les contrôleurs,
ce calcul est basé sur deux points :.
 Les spécifications sont dans le domaine temporel : temps de réponse, temps de monté et
dépassement.

 On approxime les systèmes d’ordre n>2 par un système du 2ème ordre, en considérant les
deux pôles dominants.
Exemple :
Un système d’ordre n a n pôles :
N ( p) N ( p)
H ( p)  
D ( p) an p n  a n 1 p n 1  ...  a0
Les pôles dominants sont les pôles les plus proches de l’axe des imaginaires.
Im p
Réel p
Pôles dominants
pi=-σi+jωd
Dans la réponse libre y (t )  Ae i
pi t
ce sont ces pôles qui influent le plus sur la réponse du
système. les pôles éloignés n’influent pas beaucoup, ils ont une réponse rapide.
On va faire le calcul de contrôleurs sur la base des pôles dominants.
IV. Calcul des régulateurs par la méthode du lieu des pôles

1. Position du problème
Soit un système de fonction de transfert :
Moteur
Réacteur
G(p)
Réservoir
Chaudière
Notre objectif est de construire un système de commande en boucle fermée avec un contrôleur de
fonction de transfert Gc(p)

e(t) y(t)
+
r(t) Gc(p) G(p)
-
Tel que la réponse indicielle de y(t) à r(t) soit définie par les paramètres suivants : temps de
réponse, temps de monté et dépassement par rapport à un système du second ordre.
Le système en boucle fermée H(p) ;
Gc ( p )G ( p )
H ( p)  , n’est pas du 2ème ordre, on prend uniquement les deux pôles dominants.
1  Gc ( p)G ( p)
De plus, on désir que les erreurs en régime permanent soient égales à des valeurs fixées (bien
définies), ce problème est résolut avec la structure de Gc(p) connue.
2. Méthodologie de solution par lieu des pôles
a. Structure de Gc(p)
D’une façon générale s’il s’agit de diminuer ou d’éliminer les erreurs, on utilise un régulateur PI
1  1  T1 p 
GPI ( p )  K   K  

i p   i p 

-1/T1 a
pb
Ou un compensateur à retard de phase : Gc retard ( p)  , ab
pa
a<0, b<0
S’il s’agit d’améliorer la stabilité et le régime transitoire, on utilise un régulateur PD
G PD ( p )  K   d p
pb
Ou un compensateur à avance de phase :
Gc avance ( p )  , ab
pa
a<0, b<0
GPD(p) en pratique impossible à réaliser.
S’il s’agit d’améliorer les erreurs et la stabilité on utilise un régulateur PID.
G PI ( p)  K 
1
 d p 
1  T1 p 1  T2 p 
i p i p
Gc RAV ( p) 
 p  b1   p  b2 
Ou un compensateur à avance et retard de phase
 p  a1   p  a2 
Cependant parfois, on peut utiliser un régulateur proportionnel : Gp(p)=K

Exemple 01:
On considère le système suivant:
1
G( p) 

p ( p  4) 2  16 
Trouver un contrôleur Gc(p) tel que le système en boucle fermée
e(t) y(t)
+
r(t) Gc(p) G(p)
-
Ait un coefficient d’amortissement ξ=0.5

Solution :
1) Tracer le lieu des pôles :
Les pôles de G(p) : 0, -4+4j, -4-4j
Les zéros de G(p) : ∞, ∞, ∞
Im p
ξ=0.5
jω0
K=0
-45°
60°
K=0 Réel p
8 0
60°
3
45°
K=0 -jω0

2) Déterminer les pôles du système du 2ème ordre équivalent aux performances : ξ=0.5
 n2 p1,2   n  j n 1   2
G ( p ) 2% 
p 2  2 n p   n2
Cos(θ)= ξ =0.5
θ
Les pôles du système de 2ème ordre équivalent : p1,2   n  j n 1   2

Ici ξ=0.5, donc p1,2 ne sont pas fixé mais varient suivant : cos(θ)=0.5 ⤇θ=60°
 n n
cos( )   
 2 n2  (1   2 ) n2 2
 n   (1   ) 2
3) Détermination des paramètres

On procède par étapes :
Etape 01 :
Calculer à partir des performances, les pôles du système du 2ème ordre équivalent. Pour cela les
formules suivantes sont utiles :
 
  
 
 1 2  4
Md  e  Tr  (2%) cos( )  0.5

Les pôles du 2ème ordre :
p1,2   n  j n 1   2    n  d  j n 1   2
Etape 02 :
On place les pôles p1 et p2 sur le même repère que le lieu des pôles.
Etape 03 :
Question : le lieu des pôles passe-t-il par p1 et p2
Si oui :
Cela veut dire que p1 et p2 sont des racines de 1+KG(p)=0
Alors il existe un gain K>0 qui satisfait les spécifications. Donc un régulateur proportionnel suffit.

D’une façon générale :

1) Tracer le lieu des pôles pour
KG ( p)
H ( p) 
1  KG ( p )
2) Etant donner les spécifications, on considère un système de 2ème ordre, on calcul ps1, ps2
3) Placer ps1 et ps2 dans le plan complexe.
4) Si ps1, ps2 appartiens au lieu des pôles alors Ǝ K/ 1+KG(ps1)=0 alors un régulateur
proportionnel suffit pour satisfaire les spécifications :
1
K
G ( p s1 )
Si ps1, ps2 n’appartiens pas au lieu des pôles alors il faut un contrôleur dynamique :
pb
Gc ( p ) 
pa
Pour déformer le lieu des pôles et le faire passer par le nouveau lieu des pôles de 1+KG(p)Gc(p)=0
Comment choisir a et b ?
Pour cela :
On sait que si le lieu des pôles ne passe pas par ps1 (ps2)
G ( p s1 )   180
Alors nous allons calculer l’angle θ qui manque pour faire passer le lieu des pôles par ps1 (ps2).
 Si θ>0 le contrôleur doit nous apporter un angle θ°>0 et donc :
G ( p)   0,  p  b    p  a   0
 p  b    p  a    
C’est un régulateur PD ou un contrôleur à avance de phase.
pb
Gc avance ( p )  , ab
pa
On fixe b et on calcul a.
 Si θ<0 le contrôleur doit nous apporter un angle θ°<0 et donc :
G ( p)   0,  p  b    p  a   0
 p  b    p  a    
C’est un régulateur PI ou un compensateur à retard de phase

pb
Gc retard ( p)  , ab
pa
On fixe a et calcul b
Pour l’exemple : 1
G( p) 

p ( p  4) 2  16 
Calcul d’un compensateur qui donne ξ=0.5 pour H(p)
r(t) + e(t) y(t)

KG(p)
-
Im p
ξ=0.5
jω0
K=0
3.4j
-45°
60°
K=0 Réel p
8 2 0
60°
3
45°
K=0 -jω0
Solution :
1) Lieu des pôles : ok
2) Ps1 et ps2 tel que cos(θ)=ξ : c’est un ensemble de points.
Est-ce que le lieu des pôles passe par un de ces points : oui ; alors un régulateur proportionnel
1
satisfait les spécifications, il faut calculer K, K  , p0  2  j 3.4 à partir du graphe ou
G ( p0 )
par ordinateur.

Exemple 02 :
e(t) y(t)
+
r(t) Gc(p) K/(p(p+2))
Gc(p)=1+τd p
Gc(p) est un régulateur PD.
Montrer en utilisant le lieu des pôles que ce régulateur PD améliore la stabilité de ce système.
Lieu de pôle Im p
Asymptote
K=0 K=0 Réel p

-2 0
-1
1) Gc(p)=1
r(t) + e(t) y(t)
K
- p( p  2)
2) Gc(p)=1+τd p= τd((1/ τd)+p)
r(t) + e(t)  1  y(t)

  p 

K d  
d
- p( p  2)
On considère deux cas :

a) -2<-1/ τd
b) -2>-1/ τd

Exemple 03 :
Soit le système suivant :
1
p ( p  1)( p  5)
1) Soit une commande en boucle fermée avec un régulateur P de gain K.

 Tracer le lieu des pôles du système en boucle fermée pour K Є [0, ∞[.
 Trouver le gain KL limite de stabilité.
 Quelles sont les pôles complexes dominants pour K0=0.95KL
 Donner ξ et ωn, déduire le temps de réponse et le premier dépassement approximatifs
pour le système en boucle fermée avec K0.
2) Maintenant on veut avoir ωn=3.3rd/s et tr=3s.
 Un régulateur proportionnel est-il suffisant ?
Si oui calculer K.
pb
Si non, calculer un régulateur (compensateur) dynamique Gc ( p )  K
pa
3) Calculer l’erreur de trainage.
On veut réduire cette erreur de 50%
 En augmentant le gain statique (quand p→0), calculer le gain K et montrer que le
régime transitoire se détériore.
pb
 En utilisant un compensateur dynamique Gc ( p )  , calculer a et b pour ne
pa
pas détériorer le régime transitoire de 2.

Cours: AUTO 3 Asservissement Et Régulation

Transféré par

Droits d'auteur :

Formats disponibles

Cours: AUTO 3 Asservissement Et Régulation

Transféré par

Informations du document

Titre original

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

Cours: AUTO 3 Asservissement Et Régulation

Transféré par

Droits d'auteur :

Formats disponibles

Université des Frères Mentouri Constantine 1

INSTITUT DES SCIENCES ET TECHNIQUES APPLIQUEES

Année universitaire 2019/2020

L’automatique : conduite sans intervention humaine des systèmes.

- à vide 3000 tr/min

On remarque que la vitesse varie avec la charge.

Dr. ASSABAA .M Page 1

Chapitre I : Position du problème

II. Les systèmes automatiques

Objectif : niveau de liquide constant en toute situation.

2) Vitesse d’un moteur.

3) Température d’un four.

Dr. ASSABAA .M Page 2

III. Les structures de commande

1) Structure en boucle ouverte

b) Problème de la boucle ouverte

Dr. ASSABAA .M Page 3

Ainsi, on doit construire trois organes principaux ;

IV. Les moyens techniques de réalisation

Dr. ASSABAA .M Page 4

Consigne Circuit électrique Courant

c) La transmission de la décision à l’actionneur se fait par courant électrique (4mA à 20mA),

Dr. ASSABAA .M Page 5

V. Méthodologie de conception des systèmes de commande automatique

e(t) u(t) y(t)

 Déterminer le modèle mathématique qui donne la relation entre l’entrée (grandeur de

Dr. ASSABAA .M Page 6

Dr. ASSABAA .M Page 7

Chapitre II : Système et Modèle

Modèle de 4 Simplification Modèle linéaire

II. La représentation par fonction de transfert

Dr. ASSABAA .M Page 8

Le modèle algébrique obtenu s’appelle Fonction de transfert.

1) La transformée de Laplace : Définition et propriétés

Fb(t) s’appelle T.L bilatérale.

F(p) est fonction rationnelle.

Alors L(α1 f1(t)+α2 f2(t))= α1F1(p)+α2 F2(p).

ii) Décalage dans le temps

Soit f(t) et son décalage f(t-λ)

Calculons f(t-λ), t’=t-λ

Dr. ASSABAA .M Page 9

Soit f(t) et soit e-atf(t)

Soit f(t) et df/dt, L(f(t))=F(p)

C’est pour cela que p est appelé Opérateur dérivation

Dr. ASSABAA .M Page 10

On définit le produit de convolution par :

vii) Terme de la valeur finale

N(p), D(p) : polynômes.

Si les racines de D(p) (D+(p)=0) ont la partie réelle négatives

viii) Terme de la valeur initiale

Dr. ASSABAA .M Page 11

2) La transformée inverse de Laplace

F(p) est rationnelle pour f(t) continue.

La transformée inverse nous donne f(t) lorsque on a F(p).

F(p)=L(f(t)) ; L : transformée de Laplace.

f(t) =L-1(F(p)) ; L : transformée inverse de Laplace.

La transformée inverse est donnée par :

Dr. ASSABAA .M Page 12

Dr. ASSABAA .M Page 13