0% ont trouvé ce document utile (0 vote)

36 vues459 pages

Analyse des pôles de transfert

Transféré par

Nous prenons très au sérieux les droits relatifs au contenu. Si vous pensez qu’il s’agit de votre contenu, signalez une atteinte au droit d’auteur ici.

Formats disponibles

Téléchargez aux formats PDF, TXT ou lisez en ligne sur Scribd

0% ont trouvé ce document utile (0 vote)

36 vues459 pages

Analyse des pôles de transfert

Transféré par

Zakari Yae

Nous prenons très au sérieux les droits relatifs au contenu. Si vous pensez qu’il s’agit de votre contenu, signalez une atteinte au droit d’auteur ici.

Formats disponibles

Téléchargez aux formats PDF, TXT ou lisez en ligne sur Scribd

Analyse et Correction des Systèmes Continus

Florent Nageotte

Laboratoire des Sciences de l’Image de l’Informatique et de la Télédétection

UMR CNRS - UDS 7005
Equipe Automatique Vision et Robotique
Contact : Nageotte@[Link]
http : //eavr .u − [Link] /wiki/[Link]/Page personnelle de Florent Nageotte

Licence ESA 3ème année

Année Universitaire 2010/2011

(UDS) 1 / 409
Automatique

système
signaux d’entrée signaux de sortie
phénomènes
-physiques
-électroniques
-mécaniques
-chimiques
-etc.

Automatique [6] [2] [3] [1] [4] [5]

I Modélisation, représentation, analyse et correction/régulation des systèmes
I Objectif : Donner comportement désiré au système

(UDS) 2 / 409
Automatique
Synonymes
I Contrôle / commande de systèmes
I Asservissement / régulation de systèmes

Dans l’industrie
I Régulation des procédés industriels
I Contrôle / commande des machines / robots
I Systèmes avancés de régulation dans le développement des produits
technologiques de pointe (aéronautique, spatial, etc.)
I Modélisation des systèmes biologiques, économiques

Dans la vie courante

I Automobile : régulateurs de vitesse, de trajectoire
I Avions : régulateurs d’assiette
I Satellites : régulation de position, d’orientation, de vitesse
I Transports : tous les systèmes autonomes (métros, portes, etc.)
I Informatique : disques durs
(UDS) 3 / 409
Objectifs de l’asservissement

I Suivre un signal appelé ”consigne” au plus près

I La consigne correspond à l’évolution temporelle désirée du signal de sortie
I =⇒ comportement idéal : sortie = consigne (impossible en général)
I En pratique on cherche à obtenir sortie = consigne en régime permanent pour des
entrées particulières + réponse transitoire rapide du système

(UDS) 4 / 409
Automatique
Exemple
Température d’une douche
I système non réglable : un seul robinet
I système réglable manuellement : deux robinets (chaud/froid)
I système auto-régulé : mélangeur avec capteur

Objectifs
I Stabilisation de la température → stabilité
I température demandée → erreur nulle
I en peu de temps → rapidité du système
I sans passer par 50˚ → Dépassement raisonnable
Mais aussi : comportement dynamique (marge de phase, amortissement,
dépassement), rejet de perturbations, suivi de trajectoires, etc.
Commande du processus en fonction
I de la différence entre ce qui est demandé et mesuré
I des caractéristiques du système (modèle)
(UDS) 5 / 409
Nécessité des systèmes asservis

Pourquoi asservir un système ?

1. Parce que la boucle ouverte est instable
2. Parce que la BO est oscillante (peu amortie)
3. Pour améliorer le comportement du système (précision, rapidité, stabilité, etc.)

(UDS) 6 / 409
Exemples d’asservissement (position du MCC)
Moteur à courant continu
I Signal d’entrée (commande) : tension d’induit
I (Signal de) Sortie : position du moteur
I Système instable en boucle ouverte

Réponse du système

Comment stabiliser la sortie à θ = 2 rad ? Une seule solution : asservir le signal de

(UDS) 7 / 409
Exemples d’asservissement

Modèle de disque dur

I Signal d’entrée (commande) : tension à l’entrée du disque
I (Signal de) Sortie : position de la tête de lecture
I Système oscillant en boucle ouverte

(UDS) 8 / 409
Exemples d’asservissement (positionnement de la tête d’un HD)

Réponse indicielle en BO

Comment obtenir une erreur de positionnement nulle sans oscillation ? Une seule
solution : asservir le signal de position

(UDS) 9 / 409
Exemples d’asservissement (positionnement de la tête d’un HD)

Réponse indicielle en BF avec correcteur bien réglé (échelon d’amplitude 0.5)

(UDS) 10 / 409
Comportement en fonction de la sortie commandée

Le comportement d’un système bouclé dépend du signal qui est asservi

Moteur asservi en courant asservi en vitesse asservi en position

Application : porte électrique de tramway
I Comment faut-il asservir un tel dispositif ?
I Avantages / Inconvénients ?

(UDS) 11 / 409
Cadre de l’étude

Systèmes LTI mono-variables

I Systèmes linéaires invariants mono-variables (1 entrée et 1 sortie)
I Etude fréquentielle
I Etude pôles / zéros

(UDS) 12 / 409
Exemple simple de système asservi

+
A
v1 ∈ −15 ; 15

v1 - v2
 Av1 si A A
−15
v2 = −15 si v1 < A
15 si v2 > 15

A

A
+ v2 = A(v1 − v2 ) =⇒ v2 = v1 ∼ v1
A
1+A
v1 - v2 Modèle approché valable
I si A 1
I si on ne considère pas les limites
dynamiques de l’AOP (slew rate)

(UDS) 13 / 409
Automatique et automatismes

Automatique Automate programmable Automatisme

Système (parmi d’autres)

permettant de mettre en
œuvre une stratégie de
Analyse du système et commande Méthodes pour l’analyse
des besoins (cahier des + gestion des com- du fonctionnement et la
charges), choix d’une munications (entrées - programmation des auto-
méthode de commande sorties), protocole mates
+ gestion de la lo-
gique (butées, prise de
référence, etc.)

(UDS) 14 / 409
Plan du cours

Systèmes linéaires et transformée de Laplace

Réponse temporelle des systèmes linéaires

Représentation des systèmes : schémas fonctionnels et graphes de fluence

Réponse fréquentielle des systèmes linéaires

Stabilité des systèmes

Stabilité des systèmes asservis

Précision et sensibilité des systèmes asservis

(UDS) 15 / 409
Plan du cours

Introduction

Correction série

Correction parallèle

Correction par boucles imbriquées

Conclusion

(UDS) 16 / 409
Systèmes linéaires et transformée de Laplace

Plan du cours
Systèmes linéaires et transformée de Laplace
Quelques exemples de systèmes
Modélisation des systèmes
Propriétés des systèmes
Transformée de Laplace
Fonction de transfert
Pôles et zéros des systèmes linéaires
Réponse temporelle des systèmes linéaires
Représentation des systèmes : schémas fonctionnels et graphes de fluence
Schémas fonctionnels
Graphes de fluence
Réponse fréquentielle des systèmes linéaires
Fonction de transfert harmonique
Diagrammes de Bode
Diagrammes de Nyquist
Diagrammes de Black
Stabilité des systèmes
Stabilité des systèmes asservis
Critère de Nyquist
Lieu d’Evans
Précision et sensibilité des systèmes asservis
Précision des systèmes asservis
´
(UDS)
` 17 / 409
Systèmes linéaires et transformée de Laplace Quelques exemples de systèmes

Plan du cours
Systèmes linéaires et transformée de Laplace
Quelques exemples de systèmes
Modélisation des systèmes
Propriétés des systèmes
Transformée de Laplace
Fonction de transfert
Pôles et zéros des systèmes linéaires
Réponse temporelle des systèmes linéaires
Représentation des systèmes : schémas fonctionnels et graphes de fluence
Schémas fonctionnels
Graphes de fluence
Réponse fréquentielle des systèmes linéaires
Fonction de transfert harmonique
Diagrammes de Bode
Diagrammes de Nyquist
Diagrammes de Black
Stabilité des systèmes
Stabilité des systèmes asservis
Critère de Nyquist
Lieu d’Evans
Précision et sensibilité des systèmes asservis
Précision des systèmes asservis
´
(UDS)
` 18 / 409
Systèmes linéaires et transformée de Laplace Quelques exemples de systèmes

Système électrique : un transformateur

C M
R1 i2(t ) L3
i1(t )

v1(t ) L1 v2(t )
L2
R2

I Entrée : tension v1
I Sortie : tension v2

(UDS) 19 / 409
Systèmes linéaires et transformée de Laplace Quelques exemples de systèmes

Transformateur

Relation entre entrée et sortie :

1. Boucle amont Z t
1 di1 di2
v1 = R1 i1 + i1 (t)dt + L1 +M
C 0 dt dt
2. Boucle aval
di2 di1
v2 = (L2 + L3 ) +M
dt dt
3.
v2 = −R2 i2
v1 imposé et connu : 3 équations à 3 inconnues (i1 , i2 , v2 )

d 2 v1 R2 L2 + L3 dv2 d 2 v2 dv 3
MR2 2
= v2 +(R1 R2 + ) +(R1 (L2 +L3 )+L1 R2 ) 2 +(L1 (L2 +L3 )−M 2 ) 32
dt C C dt dt dt
Equation différentielle à coefficients constants d’ordre 3

(UDS) 20 / 409
Systèmes linéaires et transformée de Laplace Quelques exemples de systèmes

Système électro-mécanique : Moteur à courant continu

Convertisseur d’énergie électrique en énergie mécanique

Excitation indépendante, inducteur à aimant permanent

I Entrée : tension u
I Sortie : vitesse de rotation ω ou position angulaire θ

(UDS) 21 / 409
Systèmes linéaires et transformée de Laplace Quelques exemples de systèmes

Moteur à courant continu

I Equation électrique

di
u = Ri + L +e
dt

R L
i

u
e

(UDS) 22 / 409
Systèmes linéaires et transformée de Laplace Quelques exemples de systèmes

Moteur à courant continu

I Equation mécanique

dω
J = Γ − fω
dt

ω
u

I Equations de couplage électro-mécanique

e = Ke ω et Γ = Km i

(UDS) 23 / 409
Systèmes linéaires et transformée de Laplace Quelques exemples de systèmes

Moteur à courant continu

Relation entre entrée et sortie :

dω d 2ω
Km u = (Ke Km + Rf )ω + (Lf + RJ) + LJ 2
dt dt
Equation différentielle à coefficients constants d’ordre 2

(UDS) 24 / 409
Systèmes linéaires et transformée de Laplace Quelques exemples de systèmes

Régulation de niveau

C
h0 + ∆ h qs

I entrée : débit d’entrée qe défini par la vanne d’entrée

I sortie : hauteur de liquide dans la cuve h
I C section de la cuve
I R résistance hydraulique de la vanne de sortie

(UDS) 25 / 409
Systèmes linéaires et transformée de Laplace Quelques exemples de systèmes

Régulation de niveau
I Pour une position fixée de la vanne de sortie, le débit de sortie est relié à la
hauteur de liquide dans la cuve par une relation non-linéaire
hauteur de liquide

débit de sortie

(UDS) 26 / 409
Systèmes linéaires et transformée de Laplace Quelques exemples de systèmes

Régulation de niveau
I Pour une position donnée de la vanne de sortie, le débit de sortie est relié à la
hauteur de liquide dans la cuve par une relation non-linéaire
hauteur de liquide

h0 + ∆ h α
h0

qs0 qs0 + ∆qs débit de sortie

I Relation linéaire pour de petites variations autour d’un point de fonctionnement

∆h = R∆qs avec R = tan α

Linéarisation autour d’un point de fonctionnement
A l’état d’équilibre : h = h0 et dh = 0 ⇐⇒ qs0 = qe0
(UDS) 26 / 409
Systèmes linéaires et transformée de Laplace Quelques exemples de systèmes

Régulation de niveau

I Variation de la hauteur
Cd∆h = (qe − qs )dt
Relation entre entrée et sortie :

d∆h
C = qe0 R + ∆qe R − qs0 R − δqs R = ∆qe R − ∆qs R
dt
d∆h
=⇒ RC + ∆h = ∆qe R
dt
Equation différentielle à coefficients constants d’ordre 1 valable uniquement autour du
point de fonctionnement h0 , c’est-à-dire pour de petites variations de qe

(UDS) 27 / 409
Systèmes linéaires et transformée de Laplace Quelques exemples de systèmes

Systèmes à retard

Radiateur électrique
I Entrée : tension d’alimentation U
I Sortie : température du radiateur θ

d 2θ dθ
=a + bθ(t) + cU(t)
dt 2 dt

I Sortie : température à l’autre bout de la pièce Θ

I Temps de propagation : τ
Θ(t) = θ(t − τ )
2
d Θ(t + τ ) dΘ(t + τ )
=a + bΘ(t + τ ) + cU(t)
dt 2 dt

(UDS) 28 / 409
Systèmes linéaires et transformée de Laplace Modélisation des systèmes

Plan du cours
Systèmes linéaires et transformée de Laplace
Quelques exemples de systèmes
Modélisation des systèmes
Propriétés des systèmes
Transformée de Laplace
Fonction de transfert
Pôles et zéros des systèmes linéaires
Réponse temporelle des systèmes linéaires
Représentation des systèmes : schémas fonctionnels et graphes de fluence
Schémas fonctionnels
Graphes de fluence
Réponse fréquentielle des systèmes linéaires
Fonction de transfert harmonique
Diagrammes de Bode
Diagrammes de Nyquist
Diagrammes de Black
Stabilité des systèmes
Stabilité des systèmes asservis
Critère de Nyquist
Lieu d’Evans
Précision et sensibilité des systèmes asservis
Précision des systèmes asservis
´
(UDS)
` 29 / 409
Systèmes linéaires et transformée de Laplace Modélisation des systèmes

Obtention d’un modèle mathématique du comportement

Recherche de relations mathématiques entre signaux d’entrée(s) et de sortie(s)

Type de modèle recherché
I Modèle physique complet : limité à des systèmes simples
I Modèle physique simplifié. Ex : applications des forces en un point, négligence de
certains courants ou forces
I Modèle entrée - sortie de type boı̂te noire : pas de connaissance a priori sur le
fonctionnement du système

Obtention des paramètres du modèle

I Par calcul : connaissance de toutes les caractéristiques du système
I Par identification fréquentielle (ex : diagrammes de Bode)
I Par identification entrée - sortie : recherche des paramètres qui minimisent l’erreur
entre sortie mesurée et sortie estimée

(UDS) 30 / 409
Systèmes linéaires et transformée de Laplace Propriétés des systèmes

Plan du cours
Systèmes linéaires et transformée de Laplace
Quelques exemples de systèmes
Modélisation des systèmes
Propriétés des systèmes
Transformée de Laplace
Fonction de transfert
Pôles et zéros des systèmes linéaires
Réponse temporelle des systèmes linéaires
Représentation des systèmes : schémas fonctionnels et graphes de fluence
Schémas fonctionnels
Graphes de fluence
Réponse fréquentielle des systèmes linéaires
Fonction de transfert harmonique
Diagrammes de Bode
Diagrammes de Nyquist
Diagrammes de Black
Stabilité des systèmes
Stabilité des systèmes asservis
Critère de Nyquist
Lieu d’Evans
Précision et sensibilité des systèmes asservis
Précision des systèmes asservis
´
(UDS)
` 31 / 409
Systèmes linéaires et transformée de Laplace Propriétés des systèmes

Systèmes linéaires

u1
i k
u2
j

Définition
Système Σ à n entrées, m sorties
I Signal u1 (t) sur l’entrée i crée un signal y1 (t) sur la sortie k
I Signal u2 (t) sur l’entrée j crée un signal y2 (t) sur la sortie k
Σ est linéaire si la combinaison des signaux αu1 et βu2 produit le signal αy1 + βy2 .

Cas d’un système mono-variable

I Signal u1 (t) sur l’entrée crée un signal y1 (t) sur la sortie
I Signal u2 (t) sur l’entrée crée un signal y2 (t) sur la sortie
Σ est linéaire si le signal αu1 + βu2 produit la réponse αy1 + βy2 .
Linéarité =⇒ on peut appliquer le théorème de superposition

(UDS) 32 / 409
Systèmes linéaires et transformée de Laplace Propriétés des systèmes

Systèmes linéaires

u1 systeme
u2

Cas d’un système mono-variable

(UDS) 32 / 409
Systèmes linéaires et transformée de Laplace Propriétés des systèmes

Linéarité

I Tous les systèmes décrits par des équations différentielles sont linéaires
I Les systèmes linéaires auxquels on ajoute un retard sont linéaires

(UDS) 33 / 409
Systèmes linéaires et transformée de Laplace Propriétés des systèmes

Exemple de non-linéarités

Eléments dont les courbes caractéristiques statiques d’entrée - sortie ne sont pas des
droites passant par l’origine
I Offsets
I Saturations
I Hystérésis
I Seuils, zones mortes
I Caractéristiques non-linéaires du processus physique

(UDS) 34 / 409
Systèmes linéaires et transformée de Laplace Propriétés des systèmes

Traitement des non-linéarités

I Fonctionnement en régime linéaire

I Approximation des petites variations
I Linéarisation fréquentielle (Approximation de la 1ère harmonique)
I Linéarisation par boucle de retour

(UDS) 35 / 409
Systèmes linéaires et transformée de Laplace Propriétés des systèmes

Systèmes invariants

Définition
Soit le système Σ
Le signal u(t) sur l’entrée i crée la réponse y (t) sur la sortie j. Σ est invariant (dans le
temps) si le signal u(t − τ ) crée la réponse y (t − τ ), ∀i, j.
Les équations qui régissent le comportement physique du système sont invariantes.
Exemples de sources de non-invariance
I usure
I dépendance de la température ambiante

La non-invariance sera considérée si les relations entrée - sortie sont modifiées

pendant le fonctionnement du système.

(UDS) 36 / 409
Systèmes linéaires et transformée de Laplace Propriétés des systèmes

Invariance

I Les systèmes décrits par des équations différentielles à coefficients constants

sont des systèmes linéaires invariants
I Les systèmes à retard constant sont des systèmes invariants

(UDS) 37 / 409
Systèmes linéaires et transformée de Laplace Propriétés des systèmes

Causalité

Définition
Un système Σ est causal si les sorties ne varient pas tant que tous les signaux
d’entrée sont constants. En d’autres termes, la réponse d’un système causal ne
précède pas son excitation.

Cas d’un système mono-variable linéaire

Système causal
u(t) = 0∀t < T =⇒ y (t) = 0∀t < T

Tous les systèmes physiquement réalisables sont causaux

(UDS) 38 / 409
Systèmes linéaires et transformée de Laplace Propriétés des systèmes

Causalité

Exemple de filtre non-causal : filtre moyenneur sans déphasage

Z t+T
1
y (t) = u(τ )dτ
2T t−T

Pour T = 1s

0 si t < 0
u(t) = Γ(t) =
1 si t ≥ 0
alors y (−0.5) = 0.25

(UDS) 39 / 409
Systèmes linéaires et transformée de Laplace Propriétés des systèmes

Causalité

I Un signal f (t) est causal ssi f (t) = 0∀t < 0.

Ne pas confondre système causal et signal causal
I Tous les signaux physiques réels peuvent être rendus causaux par décalage de
l’origine de la mesure du temps

Propriétés
I La réponse d’un système causal à un signal d’entrée causal est un signal causal

(UDS) 40 / 409
Systèmes linéaires et transformée de Laplace Propriétés des systèmes

Réponse impulsionnelle d’un système

Impulsion unité
I Signal impulsion unité : u(t) = δ(t) impulsion de Dirac

δ(t) = 0 ∀t 6= 0
Z +∞ Z +
δ(t)dt = δ(t)dt = 1
−∞ −

I Signal à énergie finie dissipée en un temps infiniment court (non réalisable

physiquement)
I Propriétés : Z ∞
δ(t)f (t)dt = f (0)
−∞
Z ∞ Z ∞
δ(t) ∗ f (t) = δ(t − τ )f (τ )dτ = δ(τ )f (t − τ )dτ = f (t)
−∞ −∞

Réponse impulsionnelle
I g(t) réponse impulsionnelle du système Σ
I Réponse du système lorsqu’un signal impulsion u(t) = δ(t) est appliqué à son
entrée (UDS) 41 / 409
Systèmes linéaires et transformée de Laplace Propriétés des systèmes

Fonction de poids des systèmes linéaires invariants

La réponse d’un système linéaire mono-variable invariant à un signal d’entrée u(t)
causal est donnée par :

y (t) = g(t) ∗ u(t)

Démonstration
Z +∞ Z +∞
u(t) = δ(t − τ )u(τ )dτ = δ(t − τ )u(τ )dτ
−∞ 0

I g(t) est la réponse de Σ à δ(t)

I g(t − τ ) est la réponse de Σ à δ(t − τ )
Par linéarité Z +∞ Z t
y (t) = g(t − τ )u(τ )dτ = g(t − τ )u(τ )dτ
−∞ 0

La sortie à l’instant t est la pondération de l’entrée jusqu’à l’instant t par la réponse

impulsionnelle g(t).
I g(t) est aussi appelée fonction de poids
I Pour les systèmes invariants, la connaissance de g(t) permet de déterminer la

réponse du système à tout signal d’entrée

(UDS) 42 / 409
Systèmes linéaires et transformée de Laplace Propriétés des systèmes

Fonction de poids des systèmes linéaires invariants

g (t )

g (t − τ )

u (τ )

(UDS) 43 / 409
Systèmes linéaires et transformée de Laplace Propriétés des systèmes

Réponse indicielle

I Réponse d’un système à une entrée de type échelon unité

0 si t < 0
u(t) = Γ(t) =
1 si t ≥ 0
La réponse indicielle est l’intégrale de la réponse impulsionnelle
démonstration

Z t Z t
y(t) = g(t − τ )u(τ )dτ = g(t − τ )dτ
0 0
Z 0 Z 0
= g(τ 0 )(−dτ 0 ) = − g(τ 0 )dτ 0
t t
Z t
= g(τ 0 )dτ 0
0

(UDS) 44 / 409
Systèmes linéaires et transformée de Laplace Transformée de Laplace

Plan du cours
Systèmes linéaires et transformée de Laplace
Quelques exemples de systèmes
Modélisation des systèmes
Propriétés des systèmes
Transformée de Laplace
Fonction de transfert
Pôles et zéros des systèmes linéaires
Réponse temporelle des systèmes linéaires
Représentation des systèmes : schémas fonctionnels et graphes de fluence
Schémas fonctionnels
Graphes de fluence
Réponse fréquentielle des systèmes linéaires
Fonction de transfert harmonique
Diagrammes de Bode
Diagrammes de Nyquist
Diagrammes de Black
Stabilité des systèmes
Stabilité des systèmes asservis
Critère de Nyquist
Lieu d’Evans
Précision et sensibilité des systèmes asservis
Précision des systèmes asservis
´
(UDS)
` 45 / 409
Systèmes linéaires et transformée de Laplace Transformée de Laplace

Transformée de Laplace

Transformée de Laplace d’une fonction f (t)

Z +∞
L {f (t)} = F (s) = f (t)e−ts dt
0

I Fonction de la variable complexe s = σ + jω

I Transformée de Laplace monolatérale
I Valide sur un domaine de définition dans le plan complexe

(UDS) 46 / 409
Systèmes linéaires et transformée de Laplace Transformée de Laplace

Exemples
I f (t) = e−at Γ(t)

Z +∞ Z +∞
F (s) = e−at e−ts dt = e−(a+s)t
0 0
∞
1 −(a+s)t 1
= − e =
s+a 0 s + a

Domaine de définition : σ > −a

I f (t) = Γ(t)
1
F (s) =
s
Valable ssi σ > 0
I f (t) = tΓ(t)
1
F (s) =
s2
Valable ssi σ > 0

(UDS) 47 / 409
Systèmes linéaires et transformée de Laplace Transformée de Laplace

Propriétés de la TL

I Linéarité : L {αf (t) + βg(t)} = αL {f (t)} + βL {g(t)}

n o
I Dérivation : L dfdt(t) = sF (s) − f (0)
n o
Si f n’est pas définie en 0 L dfdt(t) = sF (s) − f (0+ ) avec f (0+ ) = limt→> 0 f (t)
nR o
t
I Intégration : L 0 f (τ )dτ = F (s)s

I Retard pour une fonction causale L {f (t − τ )} = e−τ s F (s)

Attention ! : valable que si f(t) est causale
I Produit de convolution : Si f (t) et g(t) sont causales
Z +∞ Z +∞
f (t) ∗ g(t) = f (t − τ )g(τ )dτ = f (t − τ )g(τ )dτ
−∞ 0

Z +∞
L {f (t) ∗ g(t)} = L f (t − τ )g(τ )dτ
0
= G(s)F (s)

(UDS) 48 / 409
Systèmes linéaires et transformée de Laplace Transformée de Laplace

Propriétés de la TL

I Théorème de la valeur finale :

lim f (t) = lim sF (s)

t→∞ s→0

Valable uniquement si 0 est dans le domaine de convergence de la transformée

de Laplace
I Théorème de la valeur initiale :

lim f (t) = lim sF (s)

t→0+ s→∞

Translation de la variable de Laplace : L e−at f (t) = F (s + a)

(UDS) 49 / 409
Systèmes linéaires et transformée de Laplace Transformée de Laplace

Théorème de la valeur finale

lim f (t) = lim sF (s)

t→∞ s→0

I Valable uniquement si 0 est dans le domaine de convergence de la transformée

de Laplace
I On verra qu’il faut pour cela que tous les pôles de F (s) soient à partie réelle
strictement négative, sauf un pôle possible en 0.
I Exemple : f (t) = e2t

(UDS) 50 / 409
Systèmes linéaires et transformée de Laplace Transformée de Laplace

Transformée de Laplace des fonctions périodiques

I f (t) périodique de période T , nulle pour t < 0

+∞
X
f (t) = f0 (t − nT )
n=0

avec f0 (t) non nulle sur [0, T ], nulle partout ailleurs

Z T
L {f0 (t)} = f0 (t)e−ts dt = F0 (s)
0

( +∞ ) ∞
X X
L {f (t)} = L f0 (t − nT ) = L {f0 (t − nT )}
n=0 n=0
∞
X 1
= F0 (s)e−nTs = F0 (s)
1 − e−Ts
n=0

I Exemple : TL de f (t) = sin ωt

(UDS) 51 / 409
Systèmes linéaires et transformée de Laplace Transformée de Laplace

Relation avec la transformée de Fourier

I f (t) fonction causale

R +∞ R +∞
I Transformée de Fourier de f : F {f (t)} = −∞
f (t)e−jωt dt = 0
f (t)e−jωt dt
I Pour des fonctions causales : TL =⇒s=jω TF

(UDS) 52 / 409
Systèmes linéaires et transformée de Laplace Transformée de Laplace

Transformée de Laplace inverse

3 façons de calculer l’inverse d’une transformée de Laplace

I Formule d’inversion
I Tables de transformées
I Décomposition de F (s) en éléments simples

(UDS) 53 / 409
Systèmes linéaires et transformée de Laplace Transformée de Laplace

Formule d’inversion

I
Z c+j∞
1
L−1 {F (s)} = f (t) = F (s)est ds
2πj c−j∞

avec c tel que tous les pôles de F (s) sont à gauche de la droite σ = c
I On peut en déduire que :

1
L {f (t)g(t)} = F (s) ∗ G(s)
2πj

(UDS) 54 / 409
Systèmes linéaires et transformée de Laplace Transformée de Laplace

Formule d’inversion

jω

Calcul P
par la méthode des
résidus :
f (t) = pi Res F (s)est s=p

i

Γ I Résidu en un pôle simple en s = a

n o
Res F (s)est = lim (s + a)F (s)est
0 c a s→a

σ I Résidu en un pôle de multiplicité m en s = a

n o 1 d m−1
Res F (s)est = lim (s + a)m F (s)
a s→a (m − 1)! dsm−1

Méthode peu utilisée en pratique

(UDS) 55 / 409
Systèmes linéaires et transformée de Laplace Transformée de Laplace

Utilisation des tables de transformées

(UDS) 56 / 409
Systèmes linéaires et transformée de Laplace Transformée de Laplace

Décomposition en éléments simples

Pm
j=0 bj sj N(s) N(s)
F (s) = Pn = =
i=0 ai si an Πni=1 (s − pi ) an Πki=1 (s − pi )qi

I pi pôles complexes de F (s) d’ordre de multiplicité qi

I On peut mettre F (s) sous la forme :
qi
k X
X Aij
F (s) =
(s − pi )j
i=1 j=1

A B
I somme de termes s−pi
, (s−pi )2
, etc.
I Calcul de la TL inverse de chaque terme
I Par linéarité de la TL (et donc de son inverse) on peut obtenir l’inverse de la
fonction complète

(UDS) 57 / 409
Systèmes linéaires et transformée de Laplace Transformée de Laplace

Transformée de Laplace inverse

Exemple :
4(s + 2)
F (s) =
(s + 4)(s + 1)

(UDS) 58 / 409
Systèmes linéaires et transformée de Laplace Transformée de Laplace

Utilisation de la TL pour résoudre des équations différentielles

i=n j=m
X d i f (t) X d j g(t)
ai i
= bj
dt dt j
i=0 j=0

Fonction g(t) donnée, fonction f (t) recherchée.

Résolution directe dans l’espace du temps
I Résolution de l’équation homogène pour une entrée nulle (sans terme de droite)
solutions connues pour les équations d’ordre 1 et 2
I Résolution de l’équation complète sous conditions initiales nulles
Recherche d’une solution particulière ou méthode de variation des constantes
I Calcul des coefficients inconnus à partir des conditions initiales

Pas de méthode simple et systématique pour des équations d’ordre supérieur à 2

(UDS) 59 / 409
Systèmes linéaires et transformée de Laplace Transformée de Laplace

Utilisation de la TL pour résoudre des équations différentielles

i=n j=m
X d i f (t) X d j g(t)
ai i
= bj
dt dt j
i=0 j=0

j=m j−1
i=n i−1
! !
X X d i−1−k f (t) X X d j−1−k g(t)
ai si F (s) − sk (0) = bj sj G(s) − (0)
dt i−1−k dt j−1−k
i=0 k=0 j=0 k=0

i=n i=n i−1 j=m m j−1

X X X d i−1−k f (t) X X X d j−1−k g(t)
ai si F (s) − ai sk i−1−k
(0)) = bj sj G(s) − bj (0)
dt dt j−1−k
i=0 i=1 k =0 j=0 j=1 k=0

j−1−k
k d i−1−k f
Pj=m Pj=m Pj−1
sk d dt j−1−k
g(t) Pi=n Pi−1
G(s)( j=0 bj sj ) − j=1 bj k =0 (0) i=1 ai k=0 s dt i−1−k (0)
F (s) = Pi=n + Pi=n
i=0 ai si i=0 ai s
i

dk f
F (s) est complètement définie si on dispose des n conditions initiales dt k
(0) pour
k ∈ [0, n − 1]
=⇒ calcul de f (t) par transformée de Laplace inverse

(UDS) 60 / 409
Systèmes linéaires et transformée de Laplace Transformée de Laplace

Résolution systématique

I Application de la transformée de Laplace à chaque terme de l’équation

diférentielle
I Résolution de l’équation algébrique obtenue
I Application de la transformée de Laplace inverse à la solution

(UDS) 61 / 409
Systèmes linéaires et transformée de Laplace Transformée de Laplace

Exemple de résolution

d 2y dy
+2 + y = Γ(t)
dt 2 dt
avec y (0) = −1 et y 0 (0) = 2

(UDS) 62 / 409
Systèmes linéaires et transformée de Laplace Transformée de Laplace

Exemple de résolution

d 2y dy
+2 + y = Γ(t)
dt 2 dt
avec y (0) = −1 et y 0 (0) = 2

y (t) = (1 − 2e−t )Γ(t)

(UDS) 62 / 409
Systèmes linéaires et transformée de Laplace Fonction de transfert

Plan du cours
Systèmes linéaires et transformée de Laplace
Quelques exemples de systèmes
Modélisation des systèmes
Propriétés des systèmes
Transformée de Laplace
Fonction de transfert
Pôles et zéros des systèmes linéaires
Réponse temporelle des systèmes linéaires
Représentation des systèmes : schémas fonctionnels et graphes de fluence
Schémas fonctionnels
Graphes de fluence
Réponse fréquentielle des systèmes linéaires
Fonction de transfert harmonique
Diagrammes de Bode
Diagrammes de Nyquist
Diagrammes de Black
Stabilité des systèmes
Stabilité des systèmes asservis
Critère de Nyquist
Lieu d’Evans
Précision et sensibilité des systèmes asservis
Précision des systèmes asservis
´
(UDS)
` 63 / 409
Systèmes linéaires et transformée de Laplace Fonction de transfert

Fonction de transfert

I Système Σ modélisé par une équation différentielle à coefficients constants

i=n j=m
X d i y (t) X d j u(t)
ai i
= bj
dt dt j
i=0 j=0

I Signal d’entrée causal : u(t) = 0∀t < 0

I Système Σ au repos : y (t) = 0∀t < 0

U(s)( j=m j
P
j=0 bj s )
Y (s) = Pi=n i
i=0 ai s

Fonction de transfert : Pj=m j

Y (s) j=0 bj s
F (s) = = Pi=n
U(s) i=0 ai s
i

Fraction rationnelle en s

(UDS) 64 / 409
Systèmes linéaires et transformée de Laplace Fonction de transfert

Fonction de transfert

Fonction de transfert ou transmittance d’un système mono-variable

Définition
Rapport entre la transformée de Laplace du signal de sortie et la transformée de
Laplace du signal d’entrée pour des conditions initiales (CI) toutes nulles.
I Définition valable même si T (s) n’est pas une fraction rationnelle
I Les conditions initiales sont nulles si le système est au repos (sortie nulle) et si le
signal d’entrée appliqué est causal

(UDS) 65 / 409
Systèmes linéaires et transformée de Laplace Fonction de transfert

Fonction de transfert

I La fonction de transfert d’un système linéaire invariant est la transformée de

Laplace de sa réponse impulsionnelle
Valable même si le système n’est pas représenté par une équation différentielle

Démonstration

Y (s) = L {y (t)} = L {g(t) ∗ u(t)}

= L {g(t)} U(s)

Donc
Y (s)
F (s) = = L {g(t)}
U(s)

(UDS) 66 / 409
Systèmes linéaires et transformée de Laplace Fonction de transfert

Fonction de transfert d’un système avec retard

F (s)

u (t ) v (t ) y (t )
G(s)
τ

L {y (t)} = L {v (t − τ )} = e−τ s L {v (t)}

V (s) = G(s)U(s)
=⇒ F (s) = e−τ s G(s)

(UDS) 67 / 409
Systèmes linéaires et transformée de Laplace Pôles et zéros des systèmes linéaires

Plan du cours
Systèmes linéaires et transformée de Laplace
Quelques exemples de systèmes
Modélisation des systèmes
Propriétés des systèmes
Transformée de Laplace
Fonction de transfert
Pôles et zéros des systèmes linéaires
Réponse temporelle des systèmes linéaires
Représentation des systèmes : schémas fonctionnels et graphes de fluence
Schémas fonctionnels
Graphes de fluence
Réponse fréquentielle des systèmes linéaires
Fonction de transfert harmonique
Diagrammes de Bode
Diagrammes de Nyquist
Diagrammes de Black
Stabilité des systèmes
Stabilité des systèmes asservis
Critère de Nyquist
Lieu d’Evans
Précision et sensibilité des systèmes asservis
Précision des systèmes asservis
´
(UDS)
` 68 / 409
Systèmes linéaires et transformée de Laplace Pôles et zéros des systèmes linéaires

Pôles et zéros

I Fonction de transfert du système :

Pm
N(s) bj sj
T (s) = = Pi=0
n
D(s) i=0 ai s
i

I zéros du systèmes : zi i ∈ {1, · · · , m} ∈ C/N(zi ) = 0

m zéros complexes si m est le degré de N(s)
I pôles du systèmes : pi i ∈ {1, · · · , n} ∈ C/D(pi ) = 0
n pôles complexes si n est le degré de D(s)
I Equation D(s) = 0 appelée équation caractéristique du système
I Les coefficients de la fonction de transfert sont réels =⇒ les pôles et zéros sont
réels, ou complexes conjugués. On ne peut pas trouver de pôle (ou de zéro)
complexe sans son conjugué

(UDS) 69 / 409
Systèmes linéaires et transformée de Laplace Pôles et zéros des systèmes linéaires

Pôles et zéros

Fonction de transfert sous forme pôles - zéros

bm Πmj=1 (s − zj ) Πm
j=1 (s − zj )
T (s) = n
= K n
an Πi=1 (s − pi ) Πi=1 (s − pi )
K = bamn est le gain de la forme pôles - zéros
A ne pas confondre avec le gain statique

(UDS) 70 / 409
Systèmes linéaires et transformée de Laplace Pôles et zéros des systèmes linéaires

Ordre des systèmes linéaires

Définition
L’ordre d’un système linéaire représenté par une équation différentielle est l’ordre de
l’équation différentielle donné par l’ordre de la plus haute dérivée du signal de sortie,
lorsque l’ordre minimal de la dérivation du signal d’entrée ou de sortie est égal à 0.

Exemples
I

d 2 v1 R2 L2 + L3 dv2 d 2 v2 dv 3
MR2 2
= v2 +(R1 R2 + ) +(R1 (L2 +L3 )+L1 R2 ) 2 +(L1 (L2 +L3 )−M 2 ) 32
dt C C dt dt dt

I
d∆h
RC + ∆h = ∆qe R
dt
I
d 2u du di
RC + =R
dt 2 dt dt

(UDS) 71 / 409
Systèmes linéaires et transformée de Laplace Pôles et zéros des systèmes linéaires

Pôles et zéros

I L’ordre d’un système donné par sa fonction de transfert F (s) est le nombre de
pôles de F (s) comptés avec leur ordre de multiplicité, une fois toutes les
simplifications pôles / zéros faites
Exemple :
10(s + 4)
F (s) =
(s + 10)2 s
10(s + 4)(s + 2)
F (s) =
(s + 10)2 s(s + 2)
I Les systèmes physiques réels ont au moins autant de pôles que de zéros. On dit
qu’ils sont propres. Le plus souvent ils seront strictement propres : nbre de pôles
n > nbre de zéros m. Cette propriété est liée à la vitesse finie des systèmes réels.

(UDS) 72 / 409
Systèmes linéaires et transformée de Laplace Pôles et zéros des systèmes linéaires

Représentation par la carte des pôles et des zéros

10(s + 4)
F (s) =
(s + 10 + 20j)(s + 10 − 20j)s

Im(s)

Re(s)

I La carte des pôles et des zéros contient toute l’information du système à

l’exception de son gain
(UDS) 73 / 409
Systèmes linéaires et transformée de Laplace Pôles et zéros des systèmes linéaires

Conclusion

Représentation des systèmes linéaires invariants

I Equation différentielle à coefficient constant
I Fonction de transfert
I Carte des pôles et zéros

(UDS) 74 / 409
Réponse temporelle des systèmes linéaires

Plan du cours
Systèmes linéaires et transformée de Laplace
Quelques exemples de systèmes
Modélisation des systèmes
Propriétés des systèmes
Transformée de Laplace
Fonction de transfert
Pôles et zéros des systèmes linéaires
Réponse temporelle des systèmes linéaires
Représentation des systèmes : schémas fonctionnels et graphes de fluence
Schémas fonctionnels
Graphes de fluence
Réponse fréquentielle des systèmes linéaires
Fonction de transfert harmonique
Diagrammes de Bode
Diagrammes de Nyquist
Diagrammes de Black
Stabilité des systèmes
Stabilité des systèmes asservis
Critère de Nyquist
Lieu d’Evans
Précision et sensibilité des systèmes asservis
Précision des systèmes asservis
´
(UDS)
` 75 / 409
Réponse temporelle des systèmes linéaires

Calcul de la réponse d’un système

I Résolution temporelle : résolution des équations différentielles

I Résolution temporelle : à partir de la réponse impulsionnelle y(t) = g(t) ∗ u(t)
I Résolution par la transformée de Laplace
espace temporel espace de Laplace (fréquentiel)

équations L
équations
différentielles algébriques

résolution

L−1
y (t ) Y (s)

I Résolution en utilisant la fonction de transfert (si C.I. sont nulles)

(UDS) 76 / 409
Réponse temporelle des systèmes linéaires

Fonction de transfert

Calcul de la réponse d’un système au repos donné par sa fonction de transfert F (s), à
un signal d’entrée u(t)
I Calcul de U(s) = L {u(t)}
I Calcul de Y (s) = F (s)U(s)
I Calcul de y (t) = L−1 {Y (s)}
Valable uniquement si les conditions initiales sont nulles

(UDS) 77 / 409
Réponse temporelle des systèmes linéaires

Exemples de calcul de réponse

2
I Système F (s) = s+2
signal d’entrée : u(t) = Γ(t) avec y (0) = 0
2
I Système F (s) = s+2
signal d’entrée : u(t) = Γ(t) avec y (0) = 2

(UDS) 78 / 409
Réponse temporelle des systèmes linéaires

Réponse forcée et réponse propre

I système modélisé par une équation différentielle :

i=n j=m
X d i y (t) X d j u(t)
ai i
= bj
dt dt j
i=0 j=0

Pj=m Pj=m Pj−1 k d j−1−k u(t) Pi−1 k d i−1−k y

bj sj )− j=1 bj
Pi=n
U(s)( s (0) i=1 ai s (0)
j=0 k =0 dt j−1−k k =0
Pi=n dti
i−1−k
Y (s) = Pi=n
a s i + a s
i=0 i i=0 i
= réponse forcée + réponse propre

I La réponse propre est la solution de l’équation différentielle homogène (entrée

nulle). Elle ne dépend que des conditions initiales de l’état du système (valeurs de
la sortie y et de ses dérivées à t = 0).
I La réponse forcée du système est la solution de l’équation différentielle avec
second membre lorsque le système est initialement au repos.
I La réponse totale du système est la somme de la réponse propre et de la réponse
forcée

(UDS) 79 / 409
Réponse temporelle des systèmes linéaires

Régime transitoire et régime permanent

I Régime permanent : régime atteint lorsque l’effet des conditions initiales (sur y et
sur u) n’apparaı̂t plus sur le signal de sortie

lim y (t)
t→∞

I La réponse transitoire est la partie du signal de sortie qui tend vers 0 quand
t → ∞.
Attention : il y a une réponse transitoire même lorsque les conditions initiales sont
toutes nulles
I La réponse totale du système est composée de la réponse transitoire et du régime
permanent

(UDS) 80 / 409
Réponse temporelle des systèmes linéaires

Décomposition de la réponse

2
Exemple : G(s) = s+2
, u(t) = Γ(t) et y (0) = 2

(UDS) 81 / 409
Réponse temporelle des systèmes linéaires

Caractéristiques de la réponse indicielle

1.5

105% D1 = 0.25
100% 1.2
95%
90%
1

0.5

10%
0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7

tm t1
t5%
(UDS) 82 / 409
Réponse temporelle des systèmes linéaires

Caractéristiques temporelles d’un système

Caractéristiques temporelles données sur la réponse indicielle

I Valeur finale : y∞
I Temps d’établissement ou temps de réponse à n% (souvent n=2 ou n=5)
I Temps de montée : tm
I Dépassement : D1 , en pourcentage : 100 yD1
∞

I Temps du 1er maximum : t1

I Notion de rapidité du système : tn% , tm

I Notion de stabilité du système : D1
I Notion de gain en régime permanent : y∞

(UDS) 83 / 409
Réponse temporelle des systèmes linéaires

Gain statique des systèmes

I Rapport entre amplitude du signal de sortie obtenu en régime permanent et

amplitude d’un échelon d’entrée.
I
y (t)
Gstatique = lim lorsque u(t) = UΓ(t)
t→∞ u(t)

(UDS) 84 / 409
Réponse temporelle des systèmes linéaires

Calcul du gain statique

y (t)
Gstatique = lim lorsque u(t) = UΓ(t)
t→∞ u(t)
limt→∞ y (t)
Gstatique = lorsque u(t) = UΓ(t)
U

Gstatique = lim y (t) lorsque u(t) = Γ(t)

t→∞

En utilisant le théorème de la valeur finale, limt→∞ y(t) = lims→0 sY (s). Or

Y (s) = F (s)U(s) = 1s F (s) où F (s) est la fonction de transfert du système.
Donc : Gstatique = lim F (s)
s→0

à condition que les pôles de F (s) soient tous à partie réelle strictement négative.

(UDS) 85 / 409
Réponse temporelle des systèmes linéaires

Systèmes d’ordre 1 sans zéro

I Relation entrée - sortie : y (t) + τ dy
dt
= Ku(t)
K
K
I Fonction de transfert : T (s) = 1+τ s
= τ
s+ τ1

I τ est la constante de temps du système

(UDS) 86 / 409
Réponse temporelle des systèmes linéaires

Systèmes d’ordre 1 sans zéro

I Relation entrée - sortie : y (t) + τ dy
dt
= Ku(t)
K
K
I Fonction de transfert : T (s) = 1+τ s
= τ
s+ τ1

I τ est la constante de temps du système

Exemples
i (t ) R L
i (t )

e(t ) v (t ) C
e (t ) v (t ) R

Entrée : e(t), sortie : v (t)

(UDS) 86 / 409
Réponse temporelle des systèmes linéaires

Systèmes d’ordre 1 sans zéro

I Relation entrée - sortie : y (t) + τ dy
dt
= Ku(t)
K
K
I Fonction de transfert : T (s) = 1+τ s
= τ
s+ τ1

I τ est la constante de temps du système

Exemples
i (t ) R L
i (t )

e(t ) v (t ) C
e (t ) v (t ) R

Entrée : e(t), sortie : v (t)

I 1er cas : v (t) + RC dv
dt
= e(t) alors τ = RC
I 2ème cas : v (t) + L dv
R dt
= e(t) alors τ = L
R

(UDS) 86 / 409
Réponse temporelle des systèmes linéaires

Réponse indicielle

u(t) = Γ(t) C.I. nulles de la sortie

Résolution
t
I y (t) = αe− τ + K
I C.I. : y (0) = 0
t
I y (t) = K (1 − e− τ )

Si τ > 0
I Valeur finale : K, Gain statique : K
I Dépassement : nul
I tm = 2.2τ
I t5% = 3τ
K
I Pente à l’origine τ

(UDS) 87 / 409
Réponse temporelle des systèmes linéaires

Système d’ordre 1 sans zéro

I Exemple :
2
T (s) = 1+0.1s

(UDS) 88 / 409
Si τ < 0
I Réponse divergente

I Exemple :
2
T (s) = 1−0.1s
Réponse temporelle des systèmes linéaires

Systèmes d’ordre 2 sans zéro

d2y
I Relation entrée - sortie : ωn2 y(t) + 2ζωn dy
dt
+ dt 2
= K ωn2 u(t)
On impose ζ > 0 et ωn > 0
I ζ est le facteur d’amortissement
I ωn est la pulsation naturelle du système
I Fonction de transfert :
K K ωn2
T (s) = =
1+ 2 ωζn s + 1 2
ωn2
s (s − p1 )(s − p2 )

I Pôles
p
p1,2 = (−ζ ± j 1 − ζ 2 )ωn = −ζωn ± jωd si 0 < ζ ≤ 1
p
p1,2 = (−ζ ± ζ 2 − 1)ωn si ζ ≥ 1

|p1 p2 | = ωn2
p
Si ζ < 1, ωd = ωn 1 − ζ 2 est la pseudo-pulsation √
1−ζ 2
Si ζ < 1 on a |p1,2 | = ωn et arg (p1,2 ) = π ± arctan ζ = π ± ψ avec ζ = cos ψ

(UDS) 90 / 409
Réponse temporelle des systèmes linéaires

Systèmes d’ordre 2 sans zéro

Exemple
i (t ) R L

e (t ) C
v (t )

Entrée : e(t), sortie : v (t)

2
I v (t) + RC dv
dt
+ LC ddt 2v = e(t)

(UDS) 91 / 409
Réponse temporelle des systèmes linéaires

Réponse indicielle

K K ωn2
T (s) = =
1+ 2 ωζn s + 1 2
ωn2
s (s − p1 )(s − p2 )

u(t) = Γ(t) C.I. nulles de la sortie

Résolution temporelle
I Equation caractéristique associée : ωn2 + 2ζωn x + x 2 = 0
I Solutions de l’eq. carac. :
p
x1,2 = (−ζ ± j 1 − ζ 2 )ωn si 0 < ζ ≤ 1
p
x1,2 = (−ζ ± ζ 2 − 1)ωn si ζ ≥ 1
I Solution complète :
 p
−ζωn t
 K + βe
 sin ( 1 − ζ 2 ωn t + ϕ) si 0 < ζ < 1
−ζωn t
K + (β + γt)e
√ √ si ζ = 1
2 2
K + βe(−ζ+ ζ −1)ωn t + γe(−ζ− ζ −1)ωn t

si ζ > 1


(UDS) 92 / 409
Réponse temporelle des systèmes linéaires

Système sur-amorti : ζ > 1

I Système à pôles réels négatifs :
p
p1,2 = (−ζ ± ζ 2 − 1)ωn
I Réponse indicielle :
K
y (t) = (p2 (1 − ep1 t )Γ(t) − p1 (1 − ep2 t )Γ(t))
p2 − p1
100
Exemple : T (s) = s2 +15s+50

(UDS) 93 / 409
Réponse temporelle des systèmes linéaires

Système sur-amorti : ζ > 1

Caractéristiques de la réponse indicielle
I Réponse apériodique
I Valeur finale : K, Gain statique : K
I Pas de dépassement
I Pente à l’origine nulle

Approximation du pôle dominant

T (s) = 2
0.0016s2 +0.012s+1
, ζ = 0.15, ωn = 25rad.s−1

(UDS) 99 / 409
Réponse temporelle des systèmes linéaires

Régime critique : ζ = 1

I Un pôle double réel négatif :

p1,2 = −ωn
I Réponse indicielle :
y (t) = K (1 − ep1 t + p1 tep1 t )Γ(t)
2
Exemple : T (s) = 0.0016s2 +0.08s+1

(UDS) 100 / 409

Réponse temporelle des systèmes linéaires

Régime critique : ζ = 1

Caractéristiques de la réponse indicielle

I Réponse apériodique
I Valeur finale : K, Gain statique : K
I Pas de dépassement
I Pente à l’origine nulle
3.36
I Temps de montée tm ∼ ωn
4.5
I t5% ∼ ωn

(UDS) 101 / 409

Réponse temporelle des systèmes linéaires

Effet du facteur d’amortissement

T (s) = 2
avec ωn = 25rad.s−1
1+2s ωζ +s2 12
n ωn

(UDS) 102 / 409

Réponse temporelle des systèmes linéaires

Choix de ζ
√
2
I Meilleur compromis temps de montée / dépassement : ζ = 2
= 0.707
I D < 5%
I Une seule oscillation : pas de minimum

(UDS) 103 / 409

Réponse temporelle des systèmes linéaires

Effet de la pulsation naturelle

1
T (s) = avec ζ = 0.1
1+2s ωζ +s2 12
n ωn

(UDS) 104 / 409

Réponse temporelle des systèmes linéaires

Abaques de la carte des pôles et des zéros

ζ = 0.1
100 ζ = 0.ζ3 = 0.2
ζ = 0.4 Im(s)
ζ = 0.5
80 ζ = 0.6
ζ = 0.7

60 ζ = 0.8

ζ = 0.9
40

Re(s)
0 −1
ωn = 90radω.sn −=1 70rad .s ωn = 10rad
−1 .s−1
ωn = 50radω.sn −=1 30rad .s
−20

−40

−60

−80

−100
−100 −80 −60 −40 −20 0 20 40 60 80 100

(UDS) 105 / 409

Réponse temporelle des systèmes linéaires

Systèmes d’ordre 2 avec zéro

K (1+τ s)
T (s) = 2
1+2 ωζ s+ ωs 2
n n

I Comme les pôles sont inchangés, la décomposition en éléments simples conduit

aux mêmes éléments de la réponse : K (Aep1 t + Bep2 t + 1)γ(t). Seul le poids
relatif entre les modes propres, donné par les coefficients A et B est modifié.
I =⇒ la tangente à l’origine n’est plus nulle
I La tangente à l’origine prend des valeurs négatives lorsque le zéro est positif
I Si les pôles sont réels, le zéro peut créer des dépassements mais en aucun cas
des oscillations de la sortie
I Les zéros négatifs ont tendance à diminuer le temps d’établissement du système
I Les zéros positifs ont tendance à augmenter le temps d’établissement du
système. Le retard pris au départ n’est jamais rattrappé.
I L’effet des zéros est important lorsqu’ils sont plus près de l’axe imaginaire que les
pôles

(UDS) 106 / 409

Réponse temporelle des systèmes linéaires

Effet des zéros

N(s)
T (s) = 2 avec ζ = 0.7 et ωn = 25rad.s−1
1+2 ωζ s+ ωs 2
n n

(UDS) 107 / 409

Réponse temporelle des systèmes linéaires

Exercice

Identification d’un système à partir de sa réponse indicielle

On donne la réponse indicielle suivante. Identifiez les paramètres du système en
supposant qu’il est du second ordre.

(UDS) 108 / 409

Réponse temporelle des systèmes linéaires

Exercice
Identification d’un système à partir de sa réponse indicielle
On donne la réponse indicielle suivante. Identifiez les paramètres du système en
supposant qu’il est du second ordre.

T (s) = 2
1+0.032s+0.0016s2
, ζ = 0.4, ωn = 25rad.s−1

(UDS) 108 / 409

Réponse temporelle des systèmes linéaires

Conclusion

I La forme des réponses dépend de la position des pôles du système

I Plus les pôles sont à partie réelle (négative) de valeur absolue importante, plus
les systèmes sont rapides
I La constante de temps équivalente à un pôle réel négatif est égale à l’inverse de
l’opposée du pôle : τ = − p1
I La constante de temps équivalente à une paire de pôles complexes conjugués est
égale à l’inverse de l’opposée de la partie réelle des pôles : p = −ζωn + jωd ,
τ = ζω1 n

(UDS) 109 / 409

Représentation des systèmes

Plan du cours
Systèmes linéaires et transformée de Laplace
Quelques exemples de systèmes
Modélisation des systèmes
Propriétés des systèmes
Transformée de Laplace
Fonction de transfert
Pôles et zéros des systèmes linéaires
Réponse temporelle des systèmes linéaires
Représentation des systèmes : schémas fonctionnels et graphes de fluence
Schémas fonctionnels
Graphes de fluence
Réponse fréquentielle des systèmes linéaires
Fonction de transfert harmonique
Diagrammes de Bode
Diagrammes de Nyquist
Diagrammes de Black
Stabilité des systèmes
Stabilité des systèmes asservis
Critère de Nyquist
Lieu d’Evans
Précision et sensibilité des systèmes asservis
Précision des systèmes asservis
´
(UDS)
` 110 / 409
Représentation des systèmes Schémas fonctionnels

Plan du cours
Systèmes linéaires et transformée de Laplace
Quelques exemples de systèmes
Modélisation des systèmes
Propriétés des systèmes
Transformée de Laplace
Fonction de transfert
Pôles et zéros des systèmes linéaires
Réponse temporelle des systèmes linéaires
Représentation des systèmes : schémas fonctionnels et graphes de fluence
Schémas fonctionnels
Graphes de fluence
Réponse fréquentielle des systèmes linéaires
Fonction de transfert harmonique
Diagrammes de Bode
Diagrammes de Nyquist
Diagrammes de Black
Stabilité des systèmes
Stabilité des systèmes asservis
Critère de Nyquist
Lieu d’Evans
Précision et sensibilité des systèmes asservis
Précision des systèmes asservis
´
(UDS)
` 111 / 409
Représentation des systèmes Schémas fonctionnels

Schémas fonctionnels

u (t ) y (t )
g (t )

I Système mono-entrée et mono-sortie, de

U (s) Y (s) fonction de transfert G(s) et de réponse
G (s ) impulsionnelle g(t).
I Signal d’entrée u(t) et signal de sortie y (t).
u (t ) y (t )
I On a : y (t) = g(t) ∗ u(t) et Y (s) = G(s)U(s)
G(s)

(UDS) 112 / 409

Représentation des systèmes Schémas fonctionnels

Schémas fonctionnels

u1(t ) y1(t )
G(s) y (t )
+ u (t ) +
- +
u2(t )
y2(t )

I Systèmes représentés par des blocs et leur fonction de transfert

I Signaux représentés par des traits
I Sens des signaux donné par des flèches
I Somme de signaux représentée par un sommateur
I Comparaison de signaux représentée par un soustracteur ou comparateur

(UDS) 113 / 409

Représentation des systèmes Schémas fonctionnels

Schémas fonctionnels

u1(t ) y1(t )
G(s) y (t )
+ u (t ) +
- +
u2(t )
y2(t )

Ecriture temporelle Ecriture dans l’espace de Laplace

y(t) = y1 (t) + y2 (t) Y (s) = Y1 (s) + Y2 (s)

u(t) = u1 (t) − u2 (t) U(s) = U1 (s) − U2 (s)
y1 (t) = g(t) ∗ u(t) Y1 (s) = G(s)U(s)
y(t) = g(t) ∗ (u1 (t) − u2 (t)) + y2 (t) Y (s) = G(s)(U1 (s) − U2 (s)) + Y2 (s)

(UDS) 114 / 409

Représentation des systèmes Schémas fonctionnels

Mise en cascade de blocs

v (t) = u(t) ∗ g(t)

u (t ) v (t ) y (t )
g (t ) h(t ) y (t) = v (t) ∗ h(t)
y (t) = u(t) ∗ g(t) ∗ h(t) = g(t) ∗ h(t) ∗ u(t)
U (s ) V (s) Y (s)
G(s) H (s)

V (s) = G(s)U(s)
Y (s) = H(s)V (s)
Y (s) = H(s)G(s)U(s) = G(s)H(s)U(s)
Attention

G(s)H(s) = L {g(t)} L {h(t)} = L {g(t) ∗ h(t)} 6= L {g(t)h(t)}

(UDS) 115 / 409

Représentation des systèmes Schémas fonctionnels

Forme usuelle des systèmes asservis

R (s ) (s) U (s ) Y (s)
C (s) G (s )
+
-

H (s )
Ym (s)

Signaux
I R : consigne Blocs
I U : commande I C correcteur
I Y : sortie I G procédé ou processus
I : écart I H capteur
I Ym : mesure

(UDS) 116 / 409

Représentation des systèmes Schémas fonctionnels

Forme usuelle des systèmes asservis

R (s ) (s) U (s ) Y (s)
C (s) G (s )
+
-

R (s) (s) Y (s)

G(s)
+
-

H (s )
Ym (s)

T (s )

Y (s) = G(s)(s)
Ym (s) = H(s)Y (s)
(s) = R(s) − Ym (s)
Y (s) = G(s)(R(s) − H(s)Y (s))
G(s)
Y (s) = R(s)
1 + G(s)H(s)
Y (s) G(s)
T (s) = =
R(s) 1 + G(s)H(s)
(UDS) 118 / 409
Représentation des systèmes Schémas fonctionnels

Calcul des boucles fermées

R (s) (s) Y (s)

G(s)
+
-

H (s )
Ym (s)

T (s )

Formule fondamentale des boucles fermées

I
FT (chaı̂ne directe) G(s)
FTBF (s) = =
1 + FTBO(s) 1 + G(s)H(s)
Si le signal de mesure est soustrait du signal de consigne (comparateur)
I
FT (chaı̂ne directe) G(s)
FTBF (s) = =
1 − FTBO(s) 1 − G(s)H(s)
Si le signal de mesure est ajouté au signal de consigne (additionneur)
(UDS) 119 / 409
Représentation des systèmes Schémas fonctionnels

Calcul des boucles fermées

R (s) (s) Y (s)

G(s)
+
+

H (s )
Ym (s)

T (s )

Formule fondamentale des boucles fermées

Simplification des schémas blocs

I Réduction d’une chaı̂ne
U (s ) Y (s)
G(s) H (s ) ⇐⇒ G(s)H (s)
U (s) Y (s )

I Réduction de branches en parallèle

U (s) + Y (s)
G(s)
-
⇐⇒
G(s) − H (s)
U (s) Y (s)

H (s)

I Réduction d’une boucle

U (s) Y (s)
G(s)
+
-
⇐⇒ G(s)
U (s) 1+G(s)H (s)
Y (s)

H (s)

(UDS) 120 / 409

Représentation des systèmes Schémas fonctionnels

Simplification des schémas blocs

I Déplacement d’un point de prélèvement

En amont d’un bloc
V (s)
V (s) G(s)

U (s) Y (s) U (s) Y (s )

G(s) H (s) ⇐⇒ G (s ) H (s)

En aval d’un bloc

V (s) 1
V (s) H (s)

U (s) Y (s) U (s) Y (s )

G(s) H (s) ⇐⇒ G (s ) H (s)

(UDS) 121 / 409

Représentation des systèmes Schémas fonctionnels

Simplification des schémas blocs

I Déplacement d’un bloc
En amont d’un sommateur (ou comparateur)
U (s) + Y (s) U (s) + Y (s)
G(s) H (s) ⇐⇒ G (s ) H (s)
+ +

H (s )
V (s)

V (s)

En aval d’un sommateur (ou comparateur)

U (s) + Y (s ) U (s ) Y (s)
G (s ) H (s) ⇐⇒ G(s) H (s )
+
+ +

1
V (s) G(s)

V (s )

(UDS) 122 / 409

Représentation des systèmes Graphes de fluence

Plan du cours
Systèmes linéaires et transformée de Laplace
Quelques exemples de systèmes
Modélisation des systèmes
Propriétés des systèmes
Transformée de Laplace
Fonction de transfert
Pôles et zéros des systèmes linéaires
Réponse temporelle des systèmes linéaires
Représentation des systèmes : schémas fonctionnels et graphes de fluence
Schémas fonctionnels
Graphes de fluence
Réponse fréquentielle des systèmes linéaires
Fonction de transfert harmonique
Diagrammes de Bode
Diagrammes de Nyquist
Diagrammes de Black
Stabilité des systèmes
Stabilité des systèmes asservis
Critère de Nyquist
Lieu d’Evans
Précision et sensibilité des systèmes asservis
Précision des systèmes asservis
´
(UDS)
` 123 / 409
Représentation des systèmes Graphes de fluence

Graphes de fluence

Représentation simplifiée des systèmes, utile pour les systèmes complexes.

U1 G Y1
U Y
1 1
−1 1
Y2
U2

I Signaux (variables du systèmes) représentés par des nœuds

I Systèmes représentés par des branches orientées associées à la fonction de
transfert
I Source ou nœud d’entrée : nœud qui n’a que des branches divergentes
I Puits ou nœud de sortie : nœud qui n’a que des branches convergentes

(UDS) 124 / 409

Représentation des systèmes Graphes de fluence

Graphes de fluence

U1 G Y1
U Y
1 1
−1 1
Y2
U2

Expression des puits en fonction des sources

Y (s) = Y1 (s) + Y2 (s)

Y1 (s) = G(s)U(s)
U(s) = U1 (s) − U2 (s)
Y (s) = G(s)(U1 (s) − U2 (s)) + Y2 (s)

(UDS) 125 / 409

Représentation des systèmes Graphes de fluence

Algèbre des graphes de fluence

I Addition en un nœud : la variable en un nœud est égale à la somme des signaux

convergents
X
U A C
W

D Y
V B

W = AU + BV soit W (s) = A(s)U(s) + B(s)V (s)

I Distribution depuis un nœud : la variable en un nœud est transmise sur toutes les
branches issues de ce nœud

X = CW et Y = DW

(UDS) 126 / 409

Représentation des systèmes Graphes de fluence

I Réduction de branches en série

X Y Z X AB Z
⇐⇒
A B

I Réduction de branches en parallèle

A
X Y X A+B
Y
⇐⇒
B

I Absorption d’un nœud

X
U A C U AC X
W
⇐⇒ AD
D Y BC
V B
V BD Y

(UDS) 127 / 409

Représentation des systèmes Graphes de fluence

Les graphes de fluence ne permettent pas de représenter des comparateurs. Le signe

moins du comparateur doit donc être reporté sur la fonction de transfert de la branche
convergente

R (s) (s) U (s) Y (s ) R C Y1 Y

G(s) U
C (s ) 1
+ ⇐⇒ 1 G
-
−1 H

Ym
H (s)
Ym (s)

(UDS) 128 / 409

Représentation des systèmes Graphes de fluence

Notion de parcours

E
D
W X C
A B Y Z
F

I Parcours : succession de branches dans le sens des signaux

Ex : WXYZ, XYYXZ, YX, mais pas XYZX
I Parcours ouvert : Parcours dont chacun des nœuds n’est franchi qu’une seule fois
Ex : WXZ, YXZ mais pas XYXZ
I Parcours fermé ou boucle : parcours qui aboutit au nœud de départ, mais dont
chaque autre nœud n’est franchi qu’une seule fois
Ex : YY, XYX mais pas XYYX

(UDS) 129 / 409

Représentation des systèmes Graphes de fluence

Règle de Mason

Règle de Mason (1956)

On considère un graphe de fluence ayant une entrée (ou source) U et une sortie (ou
puits) Y . La fonction de transfert T = YU vaut :
P
i Pi ∆i
T =
∆
avec :
I Pi transmittance du i eme parcours ouvert entre U et Y
I ∆ déterminant du graphe définit par :
X X X
∆=1− Bi + Bi Bj − Bi Bj Bk + · · ·
où les Bi sont les transmittances des boucles du graphe. Dans les produits de
transmittances de boucles, on ne retient que les boucles disjointes (aucun nœud
en commun
I ∆i est le déterminant du sous-graphe obtenu lorsqu’on a supprimé le parcours qui
donne Pi

(UDS) 130 / 409

Représentation des systèmes Graphes de fluence

Règle de Mason

Exemples : Déterminez la fonction de transfert du système suivant en utilisant la règle

de Mason et la simplification de schéma fonctionnel

-
R + + Y
G1 G2 G3
- +
-

(UDS) 131 / 409

Représentation des systèmes Graphes de fluence

Règle de Mason

Exemples : Déterminez la fonction de transfert du système suivant en utilisant la règle

de Mason et la simplification de schéma fonctionnel

−H2
R G2 G3
G1
1 1 A B C Y

−H1

−1

Y G1 G2 G3
Solution : R
= 1+G1 G2 H2 +G2 G3 H1 +G1 G2 G3

(UDS) 131 / 409

Réponse fréquentielle des systèmes linéaires

Plan du cours
Systèmes linéaires et transformée de Laplace
Quelques exemples de systèmes
Modélisation des systèmes
Propriétés des systèmes
Transformée de Laplace
Fonction de transfert
Pôles et zéros des systèmes linéaires
Réponse temporelle des systèmes linéaires
Représentation des systèmes : schémas fonctionnels et graphes de fluence
Schémas fonctionnels
Graphes de fluence
Réponse fréquentielle des systèmes linéaires
Fonction de transfert harmonique
Diagrammes de Bode
Diagrammes de Nyquist
Diagrammes de Black
Stabilité des systèmes
Stabilité des systèmes asservis
Critère de Nyquist
Lieu d’Evans
Précision et sensibilité des systèmes asservis
Précision des systèmes asservis
´
(UDS)
` 132 / 409
Réponse fréquentielle des systèmes linéaires Fonction de transfert harmonique

Plan du cours
Systèmes linéaires et transformée de Laplace
Quelques exemples de systèmes
Modélisation des systèmes
Propriétés des systèmes
Transformée de Laplace
Fonction de transfert
Pôles et zéros des systèmes linéaires
Réponse temporelle des systèmes linéaires
Représentation des systèmes : schémas fonctionnels et graphes de fluence
Schémas fonctionnels
Graphes de fluence
Réponse fréquentielle des systèmes linéaires
Fonction de transfert harmonique
Diagrammes de Bode
Diagrammes de Nyquist
Diagrammes de Black
Stabilité des systèmes
Stabilité des systèmes asservis
Critère de Nyquist
Lieu d’Evans
Précision et sensibilité des systèmes asservis
Précision des systèmes asservis
´
(UDS)
` 133 / 409
Réponse fréquentielle des systèmes linéaires Fonction de transfert harmonique

Réponse à une entrée sinusoı̈dale

Théorème
En régime permanent, la réponse d’un système linéaire donné par sa fonction de
transfert F (s) à un signal d’entrée sinusoı̈dal u(t) = Asin(ωt) est un signal sinusoı̈dal
d’amplitude A |F (jω)| et déphasé de ϕ(F (jω)) = arg F (jω) :

y (t) = A |F (jω)| sin(ωt + arg F (jω))

à condition que les pôles de F (s) soient tous à partie réelle négative

Remarques
I La sinusoı̈de de sortie a la même pulsation ω que la sinusoı̈de d’entrée
I |F (jω)| représente l’amplification (ou l’atténuation) de l’amplitude de la sinusoı̈de
par le système pour la pulsation ω :
Ampl(y(t))
|F (jω)| =
Ampl(u(t))
I ϕ(F (jω)) = arg F (jω) représente le déphasage entre la sinusoı̈de d’entrée et la
sinusoı̈de de sortie

(UDS) 134 / 409

Réponse fréquentielle des systèmes linéaires Fonction de transfert harmonique

Réponse à une entrée sinusoı̈dale

Démonstration
N(s) Aω
F (s) = et U(s) = 2
Πi (s − pi ) s + ω2
!
N(s) X Bi C C∗
Y (s) = Aω 2 =A + +
(s + ω 2 )Πi (s − pi ) s − pi s − jω s + jω
i
F (jω) ∗ −F (−jω) −F (jω)∗
avec C = 2j
et C = = 2j 2j
Par transformée de Laplace inverse, on obtient :

X
y (t) = Bi epi t Γ(t) + A(Cejωt + C ∗ e−jωt )Γ(t)
i
X |F (jω)| (jωt+arg F (jω))
= Bi epi t Γ(t) + A (e − e−(jωt+arg F (jω)) )Γ(t)
2j
i
X
= Bi epi t Γ(t) + A |F (jω)| sin (ωt + arg F (jω))Γ(t)
i

Si Re(pi ) < 0 ∀i, y (t) ∼t→∞ A |F (jω)| sin (ωt + arg F (jω))Γ(t)

(UDS) 135 / 409

Réponse fréquentielle des systèmes linéaires Fonction de transfert harmonique

Réponse à une entrée sinusoı̈dale

Conséquences
I La connaissance de F (jω), appelée transmittance ou fonction de transfert
harmonique est suffisante pour déterminer la réponse du système en régime
permanent à des entrées sinusoı̈dales
I Puisque tous les signaux peuvent être décomposés en une somme (finie ou
infinie) de sinusoı̈des (série de Fourier ou intégrale de Fourier), par linéarité, la
transmittance harmonique est suffisante pour déterminer le régime permanent de
la sortie du système pour n’importe quelle entrée
I On note g(t) la réponse impulsionnelle du système
Z ∞
F (jω) = g(t)e−tjω dt
0

est la transformée de Fourier de la réponse impulsionnelle du système.

(UDS) 136 / 409

Réponse fréquentielle des systèmes linéaires Fonction de transfert harmonique

Réponse fréquentielle des systèmes linéaires

I Le comportement dynamique des systèmes peut être étudié par la réponse des
systèmes à des excitations sinusoı̈dales, c’est-à-dire par leur transmittance
harmonique.
I On parle de réponse fréquentielle car on étudie le comportement fréquence par
fréquence
I Il existe de nombreuses façons de représenter la réponse fréquentielle d’un
système
I Diagrammes de Bode
I Diagrammes de Nyquist
I Diagrammes de Black
I On sépare souvent module et argument de la fonction de transfert harmonique

(UDS) 137 / 409

Réponse fréquentielle des systèmes linéaires Fonction de transfert harmonique

Fonction de transfert harmonique

K Πm
i=1 (s − zi )
F (s) = avec K > 0
Πnk=1 (s − pk )
K Πm
i=1 (jω − zi )
F (jω) =
Πnk=1 (jω − pk )

Phase
m
X n
X
arg F (jω) = arg (jω − zi ) − arg (jω − pk )
i=1 k =1

(UDS) 138 / 409

Réponse fréquentielle des systèmes linéaires Fonction de transfert harmonique

Fonction de transfert harmonique

jω

|j ω − zi |
arg (j ω − pk )

(UDS) 139 / 409

Réponse fréquentielle des systèmes linéaires Fonction de transfert harmonique

Fonction de transfert harmonique

I Effet d’un pôle réel : module maximum pour ω = 0, décroissance monotone

argument nul pour ω = 0 qui diminue jusqu’à − π2 quand ω → ∞
I Effet d’un zéro réel : module minimum pour ω = 0, croissance monotone
Si le zéro est négatif : argument nul pour ω = 0 qui augmente jusqu’à π2 quand
ω→∞
Si le zéro est positif : argument = π pour ω = 0 qui diminue jusqu’à π2 quand
ω→∞
I Effet d’une paire de pôles conjugués : augmentation du module jusqu’au voisinage
du pôle à partie imaginaire positive (ω ∼ ωd ), puis décroissance forte
Argument nul pour ω = 0 et diminuant jusqu’à −π quand ω → ∞.

(UDS) 140 / 409

Réponse fréquentielle des systèmes linéaires Diagrammes de Bode

Plan du cours
Systèmes linéaires et transformée de Laplace
Quelques exemples de systèmes
Modélisation des systèmes
Propriétés des systèmes
Transformée de Laplace
Fonction de transfert
Pôles et zéros des systèmes linéaires
Réponse temporelle des systèmes linéaires
Représentation des systèmes : schémas fonctionnels et graphes de fluence
Schémas fonctionnels
Graphes de fluence
Réponse fréquentielle des systèmes linéaires
Fonction de transfert harmonique
Diagrammes de Bode
Diagrammes de Nyquist
Diagrammes de Black
Stabilité des systèmes
Stabilité des systèmes asservis
Critère de Nyquist
Lieu d’Evans
Précision et sensibilité des systèmes asservis
Précision des systèmes asservis
´
(UDS)
` 141 / 409
Réponse fréquentielle des systèmes linéaires Diagrammes de Bode

Diagrammes de Bode

I 2 diagrammes
I Diagramme d’amplitude (ou de gain) qui représente |F (jω)| ou le plus souvent
20 log |F (jω)| en fonction de ω
I Diagramme de phase qui représente la phase ϕ(F (jω)) = arg F (jω) en fonction
de ω
I Généralement représentés avec une échelle semi-logarithmique : l’axe des
pulsations est gradué en log ω

(UDS) 142 / 409

Réponse fréquentielle des systèmes linéaires Diagrammes de Bode

Remarques
I Le diagramme en amplitude est une fonction paire :

|F (−jω)| = |F (jω)|

Le diagramme en phase est une fonction impaire :

arg F (−jω) = − arg F (jω)

=⇒ On ne trace en général que pour les valeurs de ω > 0

I Les diagrammes de Bode peuvent être tracés même pour des systèmes dont les
pôles sont à partie réelle négative. On les trace aussi lorsque le système a des
pôles en 0.
I la valeur du diagramme d’amplitude à la pulsation ω fournit l’amplification (ou
l’attenuation) de l’amplitude de la sinusoı̈de d’entrée lorsque sa pulsation vaut ω
I la valeur du diagramme de phase à la pulsation ω fournit le déphasage d’une
sinusoı̈de à la pulsation ω

(UDS) 143 / 409

Réponse fréquentielle des systèmes linéaires Diagrammes de Bode

Diagrammes de Bode : exemples

1000(s + 10)
F (s) =
(s + 100)2

(UDS) 144 / 409

Réponse fréquentielle des systèmes linéaires Diagrammes de Bode

Notes sur les diagrammes de Bode

I Echelle logarithmique ⇐⇒ échelle linéaire pour log (ω)

I une octave ⇐⇒ doublement de la pulsation (ou fréquence)
une décade ⇐⇒ multiplication de la pulsation (ou fréquence) par 10
I La courbe de gain est au-dessus de l’axe 0dB du diagramme en amplitude
lorsque le gain est supérieur à 1 et en dessous lorsque le gain est inférieur à 1

(UDS) 145 / 409

Réponse fréquentielle des systèmes linéaires Diagrammes de Bode

Forme de Bode des fonctions de transfert

K Πmj=1 (s − zj )
T (s) =
Πni=1 (s − pi )
K Πqj=1 (s − zj )Πrk=1 (s2 + 2ζk ωnk s + ωnk
2
)
=
Πti=1 (s − pi )Πvl=1 (s2 + 2ζl ωnl s + ωnl2 )
K2 Πqj=1 (1 + s
ωj
)Πrk =1 (1 + 2 ωζk s + 1
2
ωnk
s2 )
nk
=
sl Πt−l
i=1 (1 +
s
ωi
)Πvl=1 (1 + 2 ωζl s + 1
2
ωnl
s2 )
nl

I
KB Πqj=1 (1 + s
ωj
)Πrk =1 (1 + 2 ωζk s + 1
2
ωnk
s2 )
nk
T (s) =
sl Πt−l
i=1 (1 +
s
ωi
)Πvl=1 (1 + 2 ωζl s + 1
2
ωnl
s2 )
nl

est appelée forme de Bode de la fonction de transfert.

I l, le nombre de pôles en zéro du système est appelé classe du système
I KB est le gain de Bode du système. Si l = 0, alors KB est aussi le gain statique du
système. Si l 6= 0 le gain statique du système est infini.

(UDS) 146 / 409

Réponse fréquentielle des systèmes linéaires Diagrammes de Bode

Forme de Bode des fonctions de transfert

Exercice
Mettre la fonction de transfert suivante sous la forme de Bode

2(s + 4)
T (s) =
s(s + 1 + j)(s + 1 − j)

(UDS) 147 / 409

Réponse fréquentielle des systèmes linéaires Diagrammes de Bode

Diagrammes de Bode des systèmes élémentaires

|F (j ω)|dB
Constante 20 log |K |
I ω
F (s) = K =⇒ F (jω) = K
I

|F (jω)| = |K | =⇒ |F (jω)|dB = 20 log |K | arg F (j ω)

si K > 0
I ω

arg F (jω) = 0 si K ≥ 0
arg F (jω) = −π si K < 0 si K < 0
−180o

(UDS) 148 / 409

Réponse fréquentielle des systèmes linéaires Diagrammes de Bode

|F (j ω)|dB
Pôles à l’origine
1
I l = 1, pente −20dB /dec ω
1 1 l = 2, pente −40dB /dec
F (s) = =⇒ F (jω) =
sl (jω)l l = 3, pente −60dB /dec

I
arg F (j ω)
1
|F (jω)| = l =⇒ |F (jω)|dB = −20l log ω
ω ω
−90o l=1
I pente −20ldB/dec
I
−180o l=2
π
arg F (jω) = −l
2 −270o l=3

(UDS) 149 / 409

Réponse fréquentielle des systèmes linéaires Diagrammes de Bode

|F (j ω)|dB
l = 3, pente 60dB /dec

1 l = 2, pente 40dB /dec

zéros à l’origine ω
l = 1, pente 20dB /dec
I

F (s) = sl =⇒ F (jω) = (jω)l

I
arg F (j ω)
|F (jω)| = ω l =⇒ |F (jω)|dB = 20l log ω l=3
270o
pente 20ldB/dec l=2
180o
I
π
arg F (jω) = l 90o l=1
2
ω

(UDS) 150 / 409

Réponse fréquentielle des systèmes linéaires Diagrammes de Bode

Pôle réel simple négatif (passe-bas 1er ordre)

1 1
F (s) = =⇒ F (jω) =
1 + ωs 1 + jω
ω
0 0

Diagramme d’amplitude
2 !
1 ω
|F (jω)| = r 2 =⇒ |F (jω)|dB = −10 log 1 + ω0
1 + ωω
0

Approximation asymptotique :
I quand ω → 0, |F (jω)|dB ∼ 0
I quand ω → ∞, |F (jω)|dB ∼ −20 log ω + 20 log ω0 =⇒ pente −20dB/dec
I Intersection des asymptotes en (ω = ω0 , GdB = 0)
I Erreur commise en ω = ω0 : 3dB

(UDS) 151 / 409

Réponse fréquentielle des systèmes linéaires Diagrammes de Bode

Pôle réel simple négatif (passe-bas 1er ordre)

1 1
F (s) = =⇒ F (jω) =
1 + ωs 1 + jω
ω
0 0

Diagramme de phase
ω
arg F (jω) = − arctan
ω0
Approximation asymptotique grossière
I quand ω → 0, arg F (jω) ∼ 0
I quand ω → ∞, arg F (jω) ∼ − π2
I arg F (jω) = − π4 quand ω = ω0
Approximation asymptotique plus fine
ω0
I quand ω < 5
, arg F (jω) ∼ 0
I quand ω > 5ω0 , arg F (jω) ∼ − π2
log( ωω )
I quand ω ∈ [ ω50 ; 5ω0 ], arg F (jω) ∼ − π4 (1 + 0
log 5
)

(UDS) 152 / 409

Réponse fréquentielle des systèmes linéaires Diagrammes de Bode

Pôle réel simple négatif

ω0
0

−10
Ampl. (dB)

−20
pente −20dB /dec
−30

−40
−1 0 1 2 3
10 10 10 10 10
w(rad/s)

−20
phase (deg)

−40

−60

−80

−100
−1 0 1 2 3
10 10 10 10 10
w(rad/s)

(UDS) 153 / 409

Réponse fréquentielle des systèmes linéaires Diagrammes de Bode

zéro réel simple négatif (passe-haut non propre du 1er ordre)

s jω
F (s) = 1 + =⇒ F (jω) = 1 +
ω0 ω0

Diagramme d’amplitude
s 2 2 !
ω ω
|F (jω)| = 1+ =⇒ |F (jω)|dB = 10 log 1 +
ω0 ω0
Approximation asymptotique :
I quand ω → 0, |F (jω)|dB ∼ 0
I quand ω → ∞, |F (jω)|dB ∼ 20 log ω − 20 log ω0
I Intersection des asymptotes en (ω = ω0 , GdB = 0)
I Erreur commise en ω = ω0 : 3dB

(UDS) 154 / 409

Réponse fréquentielle des systèmes linéaires Diagrammes de Bode

zéro réel simple négatif (passe-haut non propre du 1er ordre)

s jω
F (s) = 1 + =⇒ F (jω) = 1 +
ω0 ω0

Diagramme de phase
ω
arg F (jω) = arctan
ω0
Approximation asymptotique grossière
I quand ω → 0, arg F (jω) ∼ 0
π
I quand ω → ∞, arg F (jω) ∼ 2
,
π
I arg F (jω) = 4
quand ω = ω0
Approximation asymptotique plus fine
ω0
I quand ω < 5
, arg F (jω) ∼ 0
π
I quand ω > 5ω0 , arg F (jω) ∼ 2
log( ωω )
I quand ω ∈ [ ω50 ; 5ω0 ], arg F (jω) ∼ π
4
(1 + 0
log 5
)

(UDS) 155 / 409

Réponse fréquentielle des systèmes linéaires Diagrammes de Bode

Zéro réel négatif

ω0
50

40
Ampl. (dB)

20
pente 20dB /dec
10

0
−1 0 1 2 3
10 10 10 10 10
w(rad/s)

100

80
phase (deg)

0
−1 0 1 2 3
10 10 10 10 10
w(rad/s)

(UDS) 156 / 409

Réponse fréquentielle des systèmes linéaires Diagrammes de Bode

zéro réel simple positif

s jω
F (s) = 1 − =⇒ F (jω) = 1 −
ω0 ω0
avec ω0 > 0
Diagramme d’amplitude
s 2 2 !
ω ω
|F (jω)| = 1+ =⇒ |F (jω)|dB = 10 log 1 +
ω0 ω0
Identique au zéro positif
Approximation asymptotique :
I quand ω → 0, |F (jω)|dB ∼ 0
I quand ω → ∞, |F (jω)|dB ∼ 20 log ω − 20 log ω0
I Intersection des asymptotes en (ω = ω0 , GdB = 0)
I Erreur commise en ω = ω0 : −3dB

(UDS) 157 / 409

Réponse fréquentielle des systèmes linéaires Diagrammes de Bode

zéro réel simple positif

s jω
F (s) = 1 − =⇒ F (jω) = 1 −
ω0 ω0

(UDS) 158 / 409

Réponse fréquentielle des systèmes linéaires Diagrammes de Bode

Zéro réel positif

40
Ampl. (dB)

20
+1
10

0
−1 0 1 2 3
10 10 10 10 10
w(rad/s)

0
ω0
−20
−1
phase (deg)

−40

−60

−80

−100
−1 0 1 2 3
10 10 10 10 10
w(rad/s)

(UDS) 159 / 409

Réponse fréquentielle des systèmes linéaires Diagrammes de Bode

Paire de pôles complexes conjugués (passe-bas 2eme ordre)

1 1
F (s) = s2
=⇒ F (jω) =
1 + 2 ωζsn + ωn2
1 − ( ωωn )2 + j 2ζω
ωn

Diagramme d’amplitude

1 2ζ 2 − 1 2 ω 4
|F (jω)| = q =⇒ |F (jω)| dB = −10 log (1 + 2 ω + 4)
ω2
(1 − ( ωωn )2 )2 + 4ζ 2 ω ωn2 ωn
2
n

Approximation asymptotique :
I quand ω → 0, |F (jω)|dB ∼ 0
I quand ω → ∞, |F (jω)|dB ∼ −40 log ω + 40 log ωn
I Intersection des asymptotes en (ω = ωn , GdB = 0)
Variation de |F (jω)|
1
I Si ζ > √
2
, |F (jω)| décroı̂t de façon monotone
p
I Si 0 ≤ ζ < √1 , |F (jω)|, présente un maximum pour ω = ωr = ωn 1 − 2ζ 2 qui
2
vaut |F (jω)|max = √1
2ζ 1−ζ 2
|F (jω)|max
On appelle facteur de résonance (ou de qualité) Q = |F (0)|
(UDS) 160 / 409
Réponse fréquentielle des systèmes linéaires Diagrammes de Bode

Paire de pôles complexes conjugués (passe-bas 2ème ordre)

Diagramme de phase

2ζ ωωn
arg F (jω) = − arctan ω2
quand ω < ωn
1− ωn2

2ζ ωωn
= −π − arctan ω2
quand ω > ωn
1− ωn2

Approximation asymptotique grossière

I quand ω → 0, arg F (jω) ∼ 0
I quand ω → ∞, arg F (jω) ∼ −π
I arg F (jω) = − π2 en ω = ωn
Approximation asymptotique plus fine
ωn
I quand ω < 5
, arg F (jω) ∼ 0
I quand ω > 5ωn , arg F (jω) ∼ −π
log( ωω )
I quand ω ∈ [ ω5n ; 5ωn ], arg F (jω) ∼ − π2 (1 + n
log 5
)

(UDS) 161 / 409

Réponse fréquentielle des systèmes linéaires Diagrammes de Bode

Paire de pôles complexes conjugués (passe-bas 2ème ordre)

ωn
20

10
ζ = 0.1
20 log Q
0
Ampl. (dB)

−10 ζ = 0.707
−20

−30

−40

−50
0 1 2
10 10 10
w(rad/s)

0
ζ = 0.1
−50
ζ = 0.707
phase (deg)

−100

−150

−200
0 1 2
10 10 10
w(rad/s)
(UDS) 162 / 409
Réponse fréquentielle des systèmes linéaires Diagrammes de Bode

Tracé des diagrammes asymptotiques

KB Πqj=1 (1 + s
ωj
)Πrk =1 (1 + 2 ωζk s + 1
2
ωnk
s2 )
nk
T (s) =
sl Πt−l
i=1 (1 +
s
ωi
)Πvl=1 (1 + 2 ωζl s + 1
2
ωnl
s2 )
nl

KB Πqj=1 (1 + jω
ωj
)Πrk=1 (1 + 2 ωζk jω − 1
2
ωnk
ω2 )
nk
T (jω) =
(jω)l Πt−l
i=1 (1 +
jω
ωi
)Πvl=1 (1 + 2 ωζl jω − 1
2
ωnl
ω2 )
nl

q r
1 + jω + 1 + 2 ζk jω − 1 ω 2 +
X X
|T (jω)|dB = |KB |dB + 2
ωj dB
ωnk ωnk dB
j=1 k=1

t−l v
1 X 1
X 1
+
+
(jω)l jω ζ 1 2
dB i=1
1 +
ωi

l=1
1 + 2 ω jω − ω2 ω
l
dB nl nl dB

q r
X jω ζkX 1
arg T (jω) = arg KB + arg 1 + + arg 1 + 2 jω − 2 ω 2 +
ωj ωnk ωnk
j=1 k =1
!  
X t−l v
1 1 X 1
arg + arg + arg  
(jω)l
i=1
1 + jω
ωi l=1
1 + 2 ωζl jω − ω12 ω 2
nl nl

(UDS) 163 / 409

Réponse fréquentielle des systèmes linéaires Diagrammes de Bode

Tracé des diagrammes asymptotiques

I Le diagramme de Bode en amplitude est la somme algébrique des diagrammes

de bode en amplitude de termes simples : gain, pôles en zéro, systèmes du 1er et
2ème ordre, zéros réels et complexes conjugués
I Le diagramme de Bode en phase est la somme algébrique des diagrammes de
bode en phase de termes simples : gain, pôles en zéro, systèmes du 1er et 2ème
ordre, zéros réels et complexes conjugués
On obtient les diagrammes asymptotiques complets en sommant les diagrammes
asymptotiques élémentaires

(UDS) 164 / 409

Réponse fréquentielle des systèmes linéaires Diagrammes de Bode

Tracé des diagrammes asymptotiques

2(s+100)
Exemple : Tracer les diagrammes de Bode asymptotiques de F (s) = s(s+10)

(UDS) 165 / 409

Réponse fréquentielle des systèmes linéaires Diagrammes de Bode

Tracé des diagrammes asymptotiques

2(s+100)
Exemple : Tracer les diagrammes de Bode asymptotiques de F (s) = s(s+10)

20
Ampl. (dB)

−20

−40

−60

−80
−1 0 1 2 3 4
10 10 10 10 10 10
w(rad/s)

−90

−100

−110
phase (deg)

−120

−130

−140

−150
−1 0 1 2 3 4
10 10 10 10 10 10
w(rad/s)

(UDS) 165 / 409

Réponse fréquentielle des systèmes linéaires Diagrammes de Bode

Systèmes à déphasage minimal / non-minimal

Définition
Un système est dit à déphasage minimal si tous ses zéros sont à partie réelle négative
ou nulle
Explication :
I Les zéros de partie réelle de valeur absolue a ont le même comportement en
amplitude quel que soit leur signe, mais les zéros négatifs créent un déphasage
de + π2 en haute fréquence alors que les zéros positifs créent un déphasage de
− π2
I =⇒ pour un diagramme d’amplitude donné, on peut avoir des diagrammes de
phase différents selon le signe des zéros
I Parmi tous les diagrammes de phase possibles, celui pour lequel le déphasage
est le plus petit en valeur absolu (donc minimal) est celui pour lequel tous les
zéros sont à partie réelle négative ou nulle
I Un système à déphasage non-minimal est plus lent que le système à déphasage
minimal équivalent de même structure. Notamment, la réponse indicielle d’un
système ayant un zéro positif a une pente à l’origine négative

(UDS) 166 / 409

Réponse fréquentielle des systèmes linéaires Diagrammes de Bode

Systèmes à déphasage minimal / non-minimal

10(s+5)(s+2) −10(s−5)(s+2) 10(s−5)(s−2)

G1 = (s+10)(s+20)(s+100)
, G2 = (s+10)(s+20)(s+100)
, G3 = (s+10)(s+20)(s+100)

(UDS) 167 / 409

Réponse fréquentielle des systèmes linéaires Diagrammes de Bode

Systèmes à déphasage minimal / non-minimal

10(s+5)(s+2) −10(s−5)(s+2) 10(s−5)(s−2)
G1 = (s+10)(s+20)(s+100)
, G2 = (s+10)(s+20)(s+100)
, G3 = (s+10)(s+20)(s+100)

(UDS) 168 / 409

Réponse fréquentielle des systèmes linéaires Diagrammes de Bode

Identification à partir des diagrammes de Bode

(UDS) 169 / 409

Réponse fréquentielle des systèmes linéaires Diagrammes de Bode

Identification à partir des diagrammes de Bode

1000 800
Exacte : G(s) = (s+10)(s+50)
et estimée : F (s) = (s+20)(s+20)

(UDS) 169 / 409

Réponse fréquentielle des systèmes linéaires Diagrammes de Bode

Diagramme de Bode en amplitude des systèmes asservis

G chaı̂ne directe et H chaı̂ne de retour

G(s)
FTBF = F (s) =
1 + G(s)H(s)

G(jω)
F (jω) =
1 + G(jω)H(jω)
|G(jω)|
|F (jω)| =
|1 + G(jω)H(jω)|
1
I Si |G(jω)| |H(jω)| ⇐⇒ |GH| 1, (en basse fréquence si G et H sont de type
1
passe-bas) |F (jω)| ∼ |H(jω)|
1
I Si |G(jω)| |H(jω)| ⇐⇒ |GH| 1, (en haute fréquence si G et H sont de type
passe-bas) |F (jω)| ∼ |G(jω)|
Dans les 2 cas le gain du système bouclé est équivalent au plus petit des gains |G(jω)|
1
et |H(jω)| . En dB, le diagramme de |F (jω)|dB est équivalent pour toute pulsation ω au
plus petit des diagrammes |G(jω)|dB et − |H(jω)|dB .

(UDS) 170 / 409

Réponse fréquentielle des systèmes linéaires Diagrammes de Bode

Diagramme de Bode en amplitude des systèmes asservis

500
G(s) = s(s+100)
et H(s) = 2

(UDS) 171 / 409

Réponse fréquentielle des systèmes linéaires Diagrammes de Bode

Notion de bande passante

I Bande passante d’un système F définie à −NdB : zone de pulsations (ou de

fréquences) dans laquelle |F (jω)|dB − |F |dBmax > −NdB
I Pour un système passe-bas elle est définie par une fréquence limite supérieure
appelée fréquence de coupure.
I Pour un système passe-bas, la largeur de la bande passante est une image de la
rapidité du système : les systèmes rapides peuvent ”transmettre” des signaux de
fréquence élevée, tandis que les systèmes lents ”coupent” ces fréquences en les
atténuant fortement (et en introduisant un déphasage important).

(UDS) 172 / 409

Réponse fréquentielle des systèmes linéaires Diagrammes de Bode

Bande passante en automatique

I Pour un système non asservi, la pulsation de coupure est la pulsation pour

laquelle |FBO (jω)| = 1, ou encore |FBO (jω)|dB = 0.
I Pour un système asservi la pulsation de coupure est la pulsation de coupure de la
boucle ouverte correspondante

Explications
I En général les boucles fermées sont de type passe-bas
I D’après l’approximation du gain harmonique des boucles fermées, la bande
passante correspond à la zone de fréquences pour laquelle le gain du système
bouclé est très proche de son gain statique

(UDS) 173 / 409

Réponse fréquentielle des systèmes linéaires Diagrammes de Bode

Bande passante en automatique

10
I chaı̂ne directe G(s) = s(s+10)
I chaı̂ne de retour H(s) = 0.1

(UDS) 174 / 409

Réponse fréquentielle des systèmes linéaires Diagrammes de Bode

Approximation du pôle dominant

K1 K2 K
Soit G(s) = G1 (s)G2 (s) = 1+ ωs 1+ ωs
= (1+ ωs )(1+ ωs )
avec ω1 ω2
1 2 1 2

Ampl. (dB) 0

−20

−40

−60

−80
−1 0 1 2 3
10 10 10 10 10
ω1 w(rad/s) ω2
0

−50
phase (deg)

−100

−150

−200
−1 0 1 2 3
10 10 10 10 10
w(rad/s)

(UDS) 175 / 409

Réponse fréquentielle des systèmes linéaires Diagrammes de Bode

Approximation du pôle dominant

I On retrouve que le comportement est donné principalement par le pôle le plus lent
p1 = −ω1 (le plus proche de l’axe imaginaire), puisque l’effet de p2 = −ω2
intervient pour des fréquences ω > ω2 pour lesquelles le pôle p1 atténue déjà
largement les signaux d’entrée. D’un point de vue harmonique, le système est
donc globalement équivalent à un système du 1er ordre ayant un pôle en p1 .
I Si on travaille à des fréquences ω < ω2 , le sous-système G2 peut être approché
par son gain statique K2 .

(UDS) 176 / 409

Réponse fréquentielle des systèmes linéaires Diagrammes de Nyquist

Plan du cours
Systèmes linéaires et transformée de Laplace
Quelques exemples de systèmes
Modélisation des systèmes
Propriétés des systèmes
Transformée de Laplace
Fonction de transfert
Pôles et zéros des systèmes linéaires
Réponse temporelle des systèmes linéaires
Représentation des systèmes : schémas fonctionnels et graphes de fluence
Schémas fonctionnels
Graphes de fluence
Réponse fréquentielle des systèmes linéaires
Fonction de transfert harmonique
Diagrammes de Bode
Diagrammes de Nyquist
Diagrammes de Black
Stabilité des systèmes
Stabilité des systèmes asservis
Critère de Nyquist
Lieu d’Evans
Précision et sensibilité des systèmes asservis
Précision des systèmes asservis
´
(UDS)
` 177 / 409
Réponse fréquentielle des systèmes linéaires Diagrammes de Nyquist

Diagrammes de Nyquist
I Représentation de Im(F (jω)) en fonction de Re(F (jω)) pour ω ∈] − ∞, +∞[.
Courbe graduée en ω
I ⇐⇒ tracé de F (jω) dans le plan complexe
I Egalement appelé diagramme polaire : la distance d’un point à l’origine est
donnée par |F (jω)| et l’angle par rapport à l’axe réel par arg F (jω)
1.5
Im(F (j ω))

0.5
Im(F(jω))

0
Re(F (j ω))
arg F (j ω)
−0.5
|F (j ω)|

−1

ω
−1.5
−0.5 0 0.5 1 1.5 2
Re(F(jω))

(UDS) 178 / 409

Réponse fréquentielle des systèmes linéaires Diagrammes de Nyquist

I Symétrique par rapport à l’axe vertical (généralement tracé pour ω > 0)
I Combinaison des 2 diagrammes de Bode sur un seul diagramme
I Courbe graduée en valeurs de ω
I Dépendance de ω non explicite (comme pour les diagrammes de Nyquist)

(UDS) 184 / 409

Réponse fréquentielle des systèmes linéaires Diagrammes de Black

Diagrammes de Black

2 2
T (s) = T (s) = avec ζ = 0.1
1 + 0.1s 1 + 2 ωζn s + 1 2
ωn2
s

et ωn = 25rad.s−1

(UDS) 185 / 409

Réponse fréquentielle des systèmes linéaires Diagrammes de Black

Diagrammes de Black

2(s + 50)
T (s) =
s(s + 10)

(UDS) 186 / 409

Réponse fréquentielle des systèmes linéaires Diagrammes de Black

Conclusion

I Réponse en régime permanent d’un système linéaire à une excitation sinusoidale

est sinusoı̈dale
I Comportement fréquentiel (ou harmonique) du système donné par |F (jω)| et
arg(F (jω))
I Plusieurs types de diagrammes : Bode, Nyquist, Black
I Notion de bande passante, de gain statique et de rapidité du système

(UDS) 187 / 409

Stabilité des systèmes

Réponse impulsionnelle et modes propres d’un système

Système donné par une fonction de transfert :
N(s) N(s)
G(s) = = p q
Πni=1 (s − pi ) Πi=1 (s − ai )Πk=1 (s − bk )2 Πrl=1 (s − cl )(s − cl∗ )
p q r
X Ai X Bk Bk0 X Cl Cl∗
= + + + +
s − ai (s − bk )2 s − bk s − cl s − cl∗
i=1 k=1 l=1

p pôles simples réels, q pôles doubles réels, r paires de pôles complexes conjugués
On pose Cl = ρl ejϕl et cl = σl + jωl

p q r
X X X ∗
g(t) = Ai eai t Γ(t) + (Bk tebk t + Bk0 ebk t )Γ(t) + (Cl ecl t + Cl∗ ecl t )Γ(t)
i=1 k =1 l=1
p q r
X X X
= Ai eai t Γ(t) + (Bk tebk t + Bk0 ebk t )Γ(t) + (2ρl eσl t cos (ωl t + ϕl ))Γ(t)
i=1 k =1 l=1

La réponse impulsionnelle est la somme pondérée de différents modes définis par la

position des pôles du système

(UDS) 189 / 409

Stabilité des systèmes

Modes propres d’un système

pôle réel simple

g(t) = Aept Γ(t)

I si p < 0, g(t) −→ 0 quand t → ∞

I si p > 0, g(t) −→ ∞ quand t → ∞
I si p = 0, g(t) = cste

pôle réel double

g(t) = (Btept + B 0 ept )Γ(t)

I si p < 0, g(t) −→ 0 quand t → ∞

I si p > 0, g(t) −→ ∞ quand t → ∞
I si p = 0, g(t) −→ ∞ quand t → ∞

(UDS) 190 / 409

Stabilité des systèmes

Modes propres d’un système

paire de pôles complexes conjugués

g(t) = CeRe(p)t cos (Im(p)t + ϕ)Γ(t)

I si Re(p) < 0, g(t) −→ 0 quand t → ∞

I si Re(p) > 0, g(t) −→ ∞ quand t → ∞
I si Re(p) = 0, g(t) ∼ C cos (Im(p)t + ϕ) quand t → ∞

Conclusion
I les modes propres associés à des pôles à partie réelle négative s’atténuent
I les modes propres associés à des pôles à partie réelle positive divergent, de
même que les modes propres associés à des pôles multiples à partie réelle nulle
I les modes propres associés à des pôles simples à partie réelle nulle ne divergent
pas mais ne s’atténuent pas

(UDS) 191 / 409

Stabilité des systèmes

Modes propres
Im(s)

pôle simple
pôle simple

pôle double Re(s)

pôle double

(UDS) 192 / 409

Stabilité des systèmes

Théorème de la valeur finale

lim y (t) = lim sY (s) = lim sG(s)U(s)

t→∞ s→0 s→0

valable quand y (t) converge à l’infini.

X
lim y (t) = lim yi (t)
t→∞ t→∞
i

Les modes yi de y convergent si les pôles de Y (s) sont à l’intérieur du cercle unité
sauf un pôle possible en 0.
Conclusion
On peut utiliser le théorème de la valeur finale si les pôles de Y (s) sont à gauche de
l’axe imaginaire à l’exception possible d’un pôle en 0.

(UDS) 193 / 409

Stabilité des systèmes

Notion de stabilité

Stabilité BIBO (Bounded Input Bounded Output)

Un système est BIBO stable si sa réponse y(t) à tout signal d’entrée u(t) borné est
également bornée :

si ∃Usup / |u(t)| < Usup ∀t alors ∃Ysup / |y (t)| < Ysup ∀t

Propriété fondamentale
Un système donné par une fonction de transfert sous forme de fraction rationnelle est
BIBO stable ssi tous les pôles de sa fonction de transfert sont à partie réelle
strictement négative.

(UDS) 194 / 409

Stabilité des systèmes

BIBO stabilité
Démonstration
Z +∞ Z t
y (t) = g(t) ∗ u(t) = g(τ )u(t − τ )dτ = g(τ )u(t − τ )dτ
−∞ 0

Soit u(t) tel que |u(t)| < Usup ∀t

I Si tous les pôles sont à partie réelle strictement négative :
p q
Z r

∞X X X
ai t bk t 0 bk t σl t
|y (t)| < Ai e Γ(t) + Bk te + Bk e Γ(t) + 2ρl e cos (ϕl + ωl )tΓ(t) U

0
i=1 k=1 l=1
< Ysup
I Si un des pôles est à partie réelle positive, l’intégrale correspondante ne converge
pas.
I Si il y a un pôle p nul, le mode correspondant est gi (t) = K et en prenant
u(t) = K1 Γ(t)
Z t Z t
y (t) = gi (τ )u(t − τ )dτ = dτ
0 0

ne converge pas quand t → ∞

(UDS) 195 / 409

Stabilité des systèmes

BIBO stabilité

Démonstration (suite)
I Si il y a une paire de pôles complexes conjugués à partie réelle nulle p et p∗ , le
mode correspondant est gi (t) = C cos (Im(p)t + ϕ). En prenant
1
C
Γ(τ ) si cos (Im(p)(t − τ ) + ϕ) ≥ 0
u(τ ) =
− C1 Γ(τ ) si cos (Im(p)(t − τ ) + ϕ) < 0

qui est bornée, Z t

y (t) = gi (t − τ )u(τ )dτ
0

ne converge pas quand t → ∞

(UDS) 196 / 409

Stabilité des systèmes

I La stabilité d’un système ne dépend que de ses pôles, pas de ses zéros
I Un système constitué de plusieurs sous-systèmes stables en cascade est stable
I Un système constitué de plusieurs sous-systèmes en cascade dont l’un est
instable est instable
I La stabilité d’un système à retard se détermine sur la fonction de transfert dont on
a enlevé l’exponentielle. En effet, le retard pur est un système stable
I Un système qui n’est pas BIBO stable n’a que des points d’équilibre instables
I On ne peut pas amener un système linéaire instable dans un état d’équilibre
différent de son état initial
I =⇒ un système instable est inutilisable en l’état, il faut le stabiliser

(UDS) 197 / 409

Stabilité des systèmes

N(s) N(s)
T (s) = = Pn
D(s) i=0 ai s
i

On peut déterminer la stabilité d’un système à partir de la position de ses pôles.

=⇒ il faut pouvoir déterminer numériquement la valeur des pôles, c’est-à-dire résoudre
l’équation caractéristique D(s) = 0
I Facile pour les systèmes d’ordre 1 et 2
I Les systèmes d’ordre 3 et plus nécessitent un outil de calcul numérique
I Critères pour déterminer la stabilité sans calculer les pôles explicitement : utile
lorsqu’on n’a pas d’outil de calcul numérique ou lorsque l’équation caractéristique
dépend d’un ou plusieurs paramètres

(UDS) 198 / 409

Stabilité des systèmes

Critère de Routh - Hurwitz

Théorème
Les zéros du polynôme caractéristique D(s) = a0 + a1 s + ... + an sn sont à partie réelle
négative ssi :
1. tous les coefficients ai sont de même signe
2. tous les éléments de la première colonne du tableau de Routh ont le même signe

Tableau de Routh
an an−2 an−4 ...
an−1 an−3 an−5 ...
a a −an an−3 an−1 an−4 −an an−5
bn−1 = n−1 n−2
a
bn−3 = an−1
bn−5 ...
n−1
b a −a b bn−1 an−5 −an−1 bn−5
cn−1 = n−1 n−3b n−1 n−3 cn−3 = bn−1
... ...
n−1
c b −b c
dn−1 = n−1 n−3c n−1 n−3 ... ... ...
n−1
... ... ... ...

(UDS) 199 / 409

Stabilité des systèmes

Critère de Routh

I Si le système est d’ordre n, le tableau de Routh a n + 1 lignes

I Le nombre de changements de signe dans la première colonne du tableau
correspond au nombre de pôles à partie réelle positive du système
I Si la première colonne du tableau de Routh a tous ses éléments de même signe,
alors tous les coefficients de D(s) sont de même signe
I Si le premier élément d’une ligne k est nul, cela indique que le système a un pôle
à partie réelle nulle.
I Si le premier élément de la ligne k est nul, on le remplace par pour continuer les
calculs.
I Une ligne peut être multipliée par un coefficient constant avant le calcul de la ligne
suivante dans le but de simplifier les calculs

(UDS) 200 / 409

Stabilité des systèmes

Critère de Routh : exemples

2(s − 2)
T (s) =
s4 + 40s3 + 625s2 + 5750s + 25000

(UDS) 201 / 409

Stabilité des systèmes

Critère de Routh : exemples

2(s − 2)
T (s) =
s4 + 40s3 + 625s2 + 5750s + 25000
Tous les coefficients sont positifs
1 625 25000
4 575 0
bn−1 = 481.25 bn−3 = 25000 0
cn−1 = 367.2078 cn−3 = 0 0
dn−1 = 25000 0 0
Le système est donc stable (pôles en −10, −20, −5 + 10j et −5 − 10j)

(UDS) 201 / 409

Stabilité des systèmes

Critère de Routh

Systèmes du 1er ordre

N(s)
T (s) =
s−p
1
−p
Stable si p < 0

Systèmes du 2ème ordre

N(s)
T (s) =
as2 + bs + c
a c
b 0
c 0
Stable si a, b et c sont de même signe ⇐⇒ ζ > 0 et ωn > 0

(UDS) 202 / 409

Stabilité des systèmes

Application aux systèmes asservis

R (s) (s) Y (s)

G(s)
+
-

H (s )
Ym (s)

T (s )

K (s + 100) 1000
G(s) = et H(s) =
s(s + 10) s + 1000

(UDS) 203 / 409

Stabilité des systèmes

Application aux systèmes asservis

R (s) (s) Y (s)

G(s)
+
-

H (s )
Ym (s)

(UDS) 207 / 409

Stabilité des systèmes asservis

R (s ) (s) Y (s )
G(s)
+
-

H (s)
Ym (s)

T (s)

G
FTBF =
1 + FBO (s)
Un système asservi est stable ssi les zéros de 1 + FBO (s) sont tous à partie réelle
strictement négative.

(UDS) 209 / 409

Stabilité des systèmes asservis

Objectifs
I Déterminer la stabilité de la boucle fermée à partir du comportement de la boucle
ouverte
I Déterminer si un système est stabilisable
I Déterminer des marges de stabilité et de robustesse

Outils
I Critère de Nyquist
I Diagrammes de Bode
I Lieu des racines

(UDS) 210 / 409

Stabilité des systèmes asservis Critère de Nyquist

Plan du cours
Systèmes linéaires et transformée de Laplace
Quelques exemples de systèmes
Modélisation des systèmes
Propriétés des systèmes
Transformée de Laplace
Fonction de transfert
Pôles et zéros des systèmes linéaires
Réponse temporelle des systèmes linéaires
Représentation des systèmes : schémas fonctionnels et graphes de fluence
Schémas fonctionnels
Graphes de fluence
Réponse fréquentielle des systèmes linéaires
Fonction de transfert harmonique
Diagrammes de Bode
Diagrammes de Nyquist
Diagrammes de Black
Stabilité des systèmes
Stabilité des systèmes asservis
Critère de Nyquist
Lieu d’Evans
Précision et sensibilité des systèmes asservis
Précision des systèmes asservis
´
(UDS)
` 211 / 409
Stabilité des systèmes asservis Critère de Nyquist

Critère de Nyquist
Système bouclé avec contre-réaction négative et sans préfiltre

R (s ) (s) Y (s )
G(s)
+
-

H (s)
Ym (s)

T (s)

Objectif
Déterminer la stabilité de la boucle fermée en analysant la boucle ouverte en régime
harmonique

(UDS) 212 / 409

Stabilité des systèmes asservis Critère de Nyquist

Critère de Nyquist

Théorème de Cauchy
Soit un contour fermé Γ du plan complexe orienté dans le sens trigonométrique.
L’image de Γ par la fonction complexe F (s) est un contour fermé du plan complexe qui
entoure l’origine dans le sens trigonométrique n fois, avec n = Z − P, où Z et P sont le
nombre de zéros et de pôles de F à l’intérieur du contour Γ.

Im
Im
Γ F (s )

(UDS) 213 / 409

Stabilité des systèmes asservis Critère de Nyquist

Contour de Nyquist

Im
Contour de Nyquist d’un système continu
Soit un système ayant pour fonction de transfert
en boucle ouverte FBO .
Le contour de Nyquist pour le système continu
FBO est l’axe imaginaire s = jω parcouru dans le Γ
sens croissant de ω → −∞ à ω → +∞, et
complété par un demi-cercle de rayon infini
entourant le demi-espace réel négatif. Dans le Re
cas où FBO a des pôles pi = jωi sur l’axe
imaginaire, le contour évite ces points par la
droite par des demi-cercles dans le sens
trigonométrique de rayon infiniment petit ρ :
z = pi + ρejθ avec θ ∈ [− π2 , + π2 ].

(UDS) 214 / 409

Stabilité des systèmes asservis Critère de Nyquist

Lieu de Nyquist

Lieu de Nyquist
Le lieu de Nyquist d’un système continu de fonction de transfert en boucle ouverte FBO
est l’image du contour de Nyquist de FBO par la fonction complexe FBO .

Im
Im
Γ
FBO (s)

Re
Re -1

(UDS) 215 / 409

Stabilité des systèmes asservis Critère de Nyquist

Propriétés du lieu de Nyquist

Propriété 1
L’image par FBO du demi-cercle de rayon infini entourant le demi-espace réel négatif
est un point M. M est l’origine si la boucle ouverte a plus de pôles que de zéros,
M = bamn (gain en haute fréquence) si la boucle ouverte a autant de zéros que de pôles.
Il suffit donc de tracer l’image de l’axe imaginaire.

Propriété 2
Le lieu de Nyquist est symétrique par rapport à l’axe réel. On se contente en général de
tracer la partie pour ω ∈ [0, +∞] et on obtient la partie pour ω ∈ [−∞, 0] par symétrie.

Propriété 3
Chaque demi-cercle du contour de Nyquist évitant un pôle de multiplicité n sur l’axe
imaginaire parcourus dans le sens trigonométrique sont transformés par FBO en n
demi-cercles de rayon infini parcourus dans le sens anti-trigonométrique.

Propriété 4
Le lieu de Nyquist est identique au diagramme de Nyquist s’il n’y a pas de pôle de
FBO (s) sur l’axe imaginaire

(UDS) 216 / 409

Stabilité des systèmes asservis Critère de Nyquist

Propriétés du lieu de Nyquist

Propriété 5
Lorsqu’on multiplie la boucle ouverte par un gain K, le lieu de Nyquist subi une
homothétie de facteur K et de centre 0.

(UDS) 217 / 409

Stabilité des systèmes asservis Critère de Nyquist

Image des demi-cercles d’évitement

Dém : Qm
j=1 (s − zj )
FBO (s) = K Qn
i=1 (s − pi )
On suppose que p1 = jω est sur l’axe imaginaire. Le demi-cercle évitant p1 s’écrit :
s1 = p1 + ρejθ avec θ ∈ [− π2 , π2 ].
Qm jθ
j=1 (p1 + ρe − zj )
lim FBO (s1 ) = lim K Qn jθ
i=1 (p1 + ρe − pi )
ρ→0 ρ→0
Qm
j=1 (p1 − zj )
= K jθ Qn
ρe i=2 (p1 − pi )

Donc
lim |FBO (s)| = ∞
ρ→0

et
ϕ(FBO (s)) ∼ ϕ0 − θ
où ϕ0 = arg K + j=1 arg (p1 − zj ) − ni=2 arg (p1 − pi ) est indépendant de θ. Rque : si
Pm P

p1 = 0, ϕ0 = 0 ou π.

(UDS) 218 / 409

Stabilité des systèmes asservis Critère de Nyquist

Critère de stabilité de Nyquist

Théorème de Nyquist
Un système continu de fonction de transfert FBO (s) en boucle ouverte est stable en
boucle fermée ssi le lieu de Nyquist de FBO entoure le point −1, appelé point critique,
dans le sens trigonométrique un nombre de fois n égal au nombre de pôles de FBO à
l’extérieur du contour de Nyquist.
Dém :
G(s)
FBF (s) =
1 + FBO (s)
Soit P et Z le nombre total de pôles et zéros de 1 + FBO . On a Z = P. L’image du
contour de Nyquist par FBO est l’image du contour par 1 + FBO translatée de −1 selon
l’axe x. Donc, d’après le théorème de Cauchy, FBO entoure −1 n = Z+ − P+ fois où Z+
et P+ sont le nombre de zéros et de pôles de 1 + FBO à l’intérieur du contour de
Nyquist. Si le système bouclé est stable, alors Z+ = Z − Z− = Z − 0 = Z et finalement
n = Z − P+ = P − P+ = P− . Les pôles de FBO et 1 + FBO sont identiques, et donc P−
est aussi le nombre de pôles de FBO à l’extérieur du contour de Nyquist. CQFD.

(UDS) 219 / 409

Stabilité des systèmes asservis Critère de Nyquist

Exemple
R (s ) (s) U (s ) Y (s)
C (s) G (s )
+
-

Stabilité des systèmes asservis Critère de Nyquist

(UDS) 222 / 409
Stabilité des systèmes asservis Critère de Nyquist

Critère du revers
Attention ! !

−1

Ne s’applique pas !
Dans le doute, appliquer le critère de Nyquist complet ! !
(UDS) 222 / 409
Stabilité des systèmes asservis Critère de Nyquist

Marges de stabilité

Limitation de la notion simple de stabilité

Pôles stables mais proches de l’axe imaginaire ne sont pas souhaitables
I Modes convergeant très lentement =⇒ réponses lentes
I En cas d’erreur sur des paramètres du système, le système peut devenir instable

Dans les conditions du critère du revers

Robustesse et marges de stabilité
I Marge de gain : pour un système bouclé stable, gain conduisant à l’instabilité
1
GM = quand ϕ(FBO (jω)) = −180˚
|FBO (jω)|
I Marge de phase : c’est l’opposé du déphasage qui ajouté à FBO rend le système
en boucle fermé instable.

ϕM = 180 + ϕ(FBO (jω)) quand |FBO (jω)| = 1

Plus la marge de gain GM et la marge de phase ϕM sont importantes plus le système

bouclé est loin de l’instabilité.

(UDS) 223 / 409

Stabilité des systèmes asservis Critère de Nyquist

Marges de stabilité
2000
Exemple : FBO (s) = (s+10)2 (s+5)

1
Im(F(jω))

0 ω −> ∞ +
ω −> 0

−1

−2

−3
−2 −1 0 1 2 3 4
Re(F(jω))

(UDS) 224 / 409

Stabilité des systèmes asservis Critère de Nyquist

Marges de stabilité
2000
Exemple : FBO (s) = (s+10)2 (s+5)

1.5

0.5
Im(F(jω))

0 ω −> ∞
ϕM 1
GM

−0.5

−1

−1.5
−1.5 −1 −0.5 0 0.5 1 1.5
Re(F(jω))

(UDS) 224 / 409

Stabilité des systèmes asservis Critère de Nyquist

Stabilité et marges sur les diagrammes de Bode

Dans les conditions du critère du revers

Marges de stabilité
I la marge de phase est la distance du diagramme de phase à −180˚ lorsque le
diagramme de gain (en dB) coupe l’axe horizontal (pour la dernière fois) (pulsation
ωc )
I la marge de gain (exprimée en dB) est la valeur du diagramme de gain (en dB)
lorsque le diagramme de phase vaut −180˚ (pour la première fois après ωc )

(UDS) 225 / 409

Stabilité des systèmes asservis Critère de Nyquist

Stabilité et marges sur les diagrammes de Bode

2000
Exemple : FBO (s) = (s+10)2 (s+5)

0 GM (dB )
−20
Ampl. (dB)

−40

−60

−80

−100

−120
−1 0 1 2 3
10 10 10 10 10
w(rad/s)

−50

−100
phase (deg)

−150
ϕM
−200

−250

−300
−1 0 1 2 3
10 10 10 10 10
w(rad/s)

(UDS) 226 / 409

Stabilité des systèmes asservis Critère de Nyquist

Stabilité et marges sur les diagrammes de Bode

Stabilité de la BF
La boucle fermée est stable si le déphasage de la BO est supérieur à −π lorsque le
gain en dB est nul (pour la plus haute des fréquences d’annulation). La marge de gain
et la marge de phase sont alors positives
=⇒ Un système d’ordre 1 ou 2 à phase minimale (pôles et zéros à partie réelle
négative avec éventuellement un pôle en 0) et à gain statique positif est stable en BF

(UDS) 227 / 409

Stabilité des systèmes asservis Critère de Nyquist

Stabilité et marges sur les diagrammes de Bode

Attention aux cas particuliers
2
FBO (s) = 500(s+500)
(s+10)3

120

100

80
Ampl. (dB)

−20
−1 0 1 2 3
10 10 10 10 10
w(rad/s)

−50
phase (deg)

−100

−150 ϕM
−200

−250
−1 0 1 2 3
10 10 10 10 10
w(rad/s)

(UDS) 228 / 409

Stabilité des systèmes asservis Critère de Nyquist

Stabilité et marges sur les diagrammes de Bode

Attention aux cas particuliers
2
FBO (s) = 500(s+500)
(s+10)3

120

100
Ampl. (dB) 80

−20
−1 0 1 2 3
10 10 10 10 10
w(rad/s)

−50
phase (deg)

−100

−150 ϕM
−200

−250
−1 0 1 2 3
10 10 10 10 10
w(rad/s)

GM = +∞ mais le système peut devenir instable lorsqu’on diminue le gain !

(UDS) 228 / 409
Stabilité des systèmes asservis Critère de Nyquist

(UDS) 229 / 409

Stabilité des systèmes asservis Critère de Nyquist

Compromis stabilité / rapidité

K
G(s) = et H(s) = 1
s(s + 10)

(UDS) 230 / 409

Stabilité des systèmes asservis Critère de Nyquist

Compromis stabilité / rapidité

L’augmentation de la bande passante de la BO (et donc de la BF) par augmentation du

gain du système s’accompagne généralement d’une diminution de la marge de phase
et de la marge de gain
=⇒ Lorsqu’on boucle un système, il est souvent nécessaire de trouver un compromis
entre stabilité et rapidité, à moins d’utiliser des correcteurs ”stabilisants” (PD, PID, etc.)

(UDS) 231 / 409

Stabilité des systèmes asservis Critère de Nyquist

Compromis stabilité / rapidité

Temporellement : compromis temps de réponse / dépassement

(UDS) 232 / 409

Stabilité des systèmes asservis Lieu d’Evans

Plan du cours
Systèmes linéaires et transformée de Laplace
Quelques exemples de systèmes
Modélisation des systèmes
Propriétés des systèmes
Transformée de Laplace
Fonction de transfert
Pôles et zéros des systèmes linéaires
Réponse temporelle des systèmes linéaires
Représentation des systèmes : schémas fonctionnels et graphes de fluence
Schémas fonctionnels
Graphes de fluence
Réponse fréquentielle des systèmes linéaires
Fonction de transfert harmonique
Diagrammes de Bode
Diagrammes de Nyquist
Diagrammes de Black
Stabilité des systèmes
Stabilité des systèmes asservis
Critère de Nyquist
Lieu d’Evans
Précision et sensibilité des systèmes asservis
Précision des systèmes asservis
´
(UDS)
` 233 / 409
Stabilité des systèmes asservis Lieu d’Evans

Lieu des racines

Système bouclé avec contre-réaction négative et sans préfiltre

Lieu d’Evans ou lieu des racines
Déterminer la position des pôles de la boucle fermée (et donc la stabilité de la boucle
fermée) en fonction du gain de la boucle ouverte K ,

R (s) (s) Y (s )
K G (s )
+
-

H (s)
Ym (s)

T (s)

(UDS) 234 / 409

Stabilité des systèmes asservis Lieu d’Evans

Lieu des racines

Système bouclé avec contre-réaction négative et sans préfiltre

Lieu d’Evans ou lieu des racines
Déterminer la position des pôles de la boucle fermée (et donc la stabilité de la boucle
fermée) en fonction du gain de la boucle ouverte K ,

R (s) (s) Y (s )
G(s)
+
-

KH (s)
Ym (s)

T (s)

(UDS) 234 / 409

Stabilité des systèmes asservis Lieu d’Evans

Lieu d’Evans

NG (s) NH (s) NGH (s)

FBO (s) = KG(s)H(s) = K =K
DG (s) DH (s) DGH (s)
KG(s) N 0 (s)
FBF (s) = =
1 + KG(s)H(s) DGH (s) + KNGH (s)
DGH (s) + KNGH (s) = 0 est l’équation caractéristique de la boucle fermée
Définition
Le lieu d’Evans (ou lieu des racines) représente la localisation dans le plan complexe
I des pôles du système en boucle fermée FBF en fonction du gain K (positif),
c’est-à-dire
I la position des zéros du polynôme caractéristique du système : DGH + KNGH ou
encore
I les racines de l’équation 1 + FBO (s) = 0 en fonction de K.

(UDS) 235 / 409

Stabilité des systèmes asservis Lieu d’Evans

Lieu d’Evans (suite)

Le lieu d’Evans se trace à partir de la position des pôles et zéros de la boucle ouverte
et de règles simples
Qm
(s − zi )
FBO (s) = K Qni=1
j=1 (s − pj )
Qm
(s − zi ) 1
1 + FBO (s) = 0 ⇐⇒ Qni=1 =−
j=1 (s − p j ) K
Un point M d’affixe sM appartient au lieu d’Evans ssi :
1. Q
m (s − z )
i=1 M i 1
Qn =
j=1 (sM − pj ) K
(condition des modules)
2.
n
X n
X
arg(sM − zi ) − arg(sM − pj ) = (2λ + 1)π avec λ ∈ N
i=1 j=1

(condition des angles)

(UDS) 236 / 409

Stabilité des systèmes asservis Lieu d’Evans

Règles de construction
I R1 nombre de branches = n.
Points de départ pour K = 0 : les pôles de FBO (pj )
Points d’arrivée pour K → ∞ : les zéros de FBO (zi )
n − m branches partent à l’infini
I R2 Lieu symétrique par rapport à l’axe réel
I R3 Les points de l’axe réel qui appartiennent au lieu des racines ont un nombre de
points singuliers (pôles et zéros comptés avec leur ordre de multiplicité) à leur
droite et sur l’axe réel impair.
I R4 Asymptotes : les n-m asymptotes des branches partant vers l’infini font avec
l’axe réel des angles :
(2λ + 1)π
αλ = avec λ ∈ [0, n − m − 1]
n−m
Les asymptotes s’intersectent sur l’axe réel au point d’abscisse :
Pn Pm
j=1 pj − i=1 zi
σ=
n−m

(UDS) 237 / 409

Stabilité des systèmes asservis Lieu d’Evans

Règles de construction (suite)

I R5 Points de séparation de branches :
∂(1+FBO (s)) d(FBO (s))
Condition Nécessaire ∂s
= ds
=0
La condition n’est pas suffisante et il faut aussi vérifier que le point appartient au
lieu d’Evans : (1 + FBO (s) = 0)
Angle entre les demi-branches au point d’intersection : Nπ où N est le nombre de
branches qui s’intersectent (2 dans la très grande majorité des cas)
I R6 Angle de départ avec l’axe réel d’une branche en un pôle complexe pk de
multiplicité nk :
m n
1 X X
βk = ( arg(pk − zi ) − arg(pk − pj ) − (2λ + 1)π) avec λ ∈ [0, nk − 1]
nk
i=1 j=1,j6=k

Angle d’arrivée avec l’axe réel d’une branche en un zéro complexe zk de

multiplicité mk :
n m
1 X X
γk = ( arg(zk − pj ) − arg(zk − zi ) − (2λ + 1)π) avec λ ∈ [0, mk − 1]
mk
j=1 i=1,i6=k

(UDS) 238 / 409

Stabilité des systèmes asservis Lieu d’Evans

Règles de construction (fin)

I R7 Graduation du lieu en K : au point sM , K vaut
1
K =−
G(sM )H(sM )
I R8 Intersection du lieu avec l’axe imaginaire : utilisation du critère de Routh pour
déterminer les valeurs de K limite ou résolution de 1 + FBO (jω) = 0

Après la règle 7 on peut tracer le lieu d’Evans à main levée en l’orientant des pôles de
la BO vers les zéros de la BO ou vers les asymptotes en suivant l’augmentation du
gain K aux points particuliers où il a été calculé par la règle 7.
Détermination de la stabilité
Le système en boucle fermée est stable pour les valeurs de K pour lesquelles tous les
pôles de la boucle fermée sont à gauche de l’axe imaginaire.

(UDS) 239 / 409

Stabilité des systèmes asservis Lieu d’Evans

Exemple
R (s ) (s) U (s ) Y (s)
C (s) G (s )
+
-

H (s )
Ym (s)

10K (s + 5)
C(s) = avec K > 0
s(s + 100)
100
G(s) =
s + 10
H(s) = 10

(UDS) 240 / 409

Stabilité des systèmes asservis Lieu d’Evans

Exemple
I R1 3 branches, 2 branches vers l’infini (2 asymptotes)
I R3 Lieu sur l’axe réel : [−100, −10] et [−5, 0]
I R4 Angles des asymptotes avec l’axe réel : ± π2
Les asymptotes s’intersectent sur l’axe réel au point d’abscisse :
0 − 100 − 10 − (−5)
σ= = −52.5
2
I R5
dFBO −2s3 − 125s2 − 1100s − 5000
=
ds D2
dFBO
= 0 ⇐⇒ s = −53.01 ou s = −4.7425 ± 4.9664j
ds
Seul s = −53.01 est valide
I R6 Ne s’applique pas ici
I R7 En s = −53.01, K = 0.2232
I R8 le lieu d’Evans est entièrement dans le demi-espace réel négatif. Le système
en BF est donc stable ∀K > 0 (à vérifier avec le critère de Routh)

(UDS) 240 / 409

Stabilité des systèmes asservis Lieu d’Evans

Exemple

(UDS) 240 / 409

Stabilité des systèmes asservis Lieu d’Evans

Lieu d’Evans
Quelques propriétés
I Sur une portion du lieu d’Evans entre 2 pôles réels sans zéro intermédiaire le
nombre de points de séparation est nécéssairement impair (1, 3, etc.), puisque les
branches vont des pôles vers les zéros
I Sur une portion du lieu d’Evans entre un pôle réel et un zéro réel le nombre de
points de séparation est nécéssairement pair (0, 2, etc.), puisque les branches
vont des pôles vers les zéros

Quelques cas simples

I si le système a seulement 2 pôles réels, le lieu d’Evans est exactement une croix,
le point de séparation des branches se trouve exactement au milieu des 2 pôles
I si le système a 2 pôles réels et un zéro réel plus à gauche, le lieu d’Evans est
constitué de la partie entre les pôles et de l’axe réel négatif à partir du zéro, ces 2
parties étant reliées entre elles par un cercle
I si le système a un pôle, puis un zéro, puis un pôle en s’éloignant de l’axe
imaginaire vers les réels négatifs, le lieu d’Evans est constitué de la portion entre
le premier pôle et le zéro et de l’axe réel négatif à partir du deuxième pôle.
I On retrouve dans ces cas le fait qu’un système du second ordre stable en BO
conduit à un système stable en BF
(UDS) 241 / 409
Stabilité des systèmes asservis Lieu d’Evans

Conclusion

On peut déterminer la stabilité d’une boucle fermée à partir du comportement de la

boucle ouverte
I A partir du comportment en régime harmonique : lieu de Nyquist
I A partir d’une analyse algébrique de la position des pôles : lieu d’Evans
I Attention : pour pouvoir tracer le lieu d’Evans ou appliquer le critère de Routh, il
faut que la FTBO puisse se mettre sous forme de fraction rationnelle =⇒ On ne
peut pas tracer le lieu d’Evans et appliquer le critère de Routh à un système à
retard
I Le critère de Nyquist (ou sa transposition dans le diagramme de Bode) est la
seule méthode permettant de déterminer la stabilité d’un système à retard asservi

(UDS) 242 / 409

Précision et sensibilité des systèmes asservis

Plan du cours
Systèmes linéaires et transformée de Laplace
Quelques exemples de systèmes
Modélisation des systèmes
Propriétés des systèmes
Transformée de Laplace
Fonction de transfert
Pôles et zéros des systèmes linéaires
Réponse temporelle des systèmes linéaires
Représentation des systèmes : schémas fonctionnels et graphes de fluence
Schémas fonctionnels
Graphes de fluence
Réponse fréquentielle des systèmes linéaires
Fonction de transfert harmonique
Diagrammes de Bode
Diagrammes de Nyquist
Diagrammes de Black
Stabilité des systèmes
Stabilité des systèmes asservis
Critère de Nyquist
Lieu d’Evans
Précision et sensibilité des systèmes asservis
Précision des systèmes asservis
´
(UDS)
` 243 / 409
Précision et sensibilité des systèmes asservis Précision des systèmes asservis

Plan du cours
Systèmes linéaires et transformée de Laplace
Quelques exemples de systèmes
Modélisation des systèmes
Propriétés des systèmes
Transformée de Laplace
Fonction de transfert
Pôles et zéros des systèmes linéaires
Réponse temporelle des systèmes linéaires
Représentation des systèmes : schémas fonctionnels et graphes de fluence
Schémas fonctionnels
Graphes de fluence
Réponse fréquentielle des systèmes linéaires
Fonction de transfert harmonique
Diagrammes de Bode
Diagrammes de Nyquist
Diagrammes de Black
Stabilité des systèmes
Stabilité des systèmes asservis
Critère de Nyquist
Lieu d’Evans
Précision et sensibilité des systèmes asservis
Précision des systèmes asservis
´
(UDS)
` 244 / 409
Précision et sensibilité des systèmes asservis Précision des systèmes asservis

Dénomination (rappel)

R (s ) (s) Y (s )
G(s)
+
-

H (s)
Ym (s)

T (s)

I R : consigne
I Y : signal de sortie
I Ym : signal de mesure
I : écart

(UDS) 245 / 409

Précision et sensibilité des systèmes asservis Précision des systèmes asservis

Précision des systèmes asservis

Capacité d’un système à faire suivre au signal de sortie la valeur du signal d’entrée
I Précision dynamique : pour un signal d’entrée quelconque
I Précision en régime permanent : une fois que le régime permanent est établi
I Précision statique : en régime permanent en réponse à une entrée de type
échelon
Dans ce cours, on ne considère que les erreurs en régime permanent
Intérêt
I Si l’erreur en régime permanent pour un type de signal est nulle, alors il suffit
d’imposer ce signal à l’entrée du système pour que la sortie du système prenne
les mêmes valeurs en régime permanent

(UDS) 246 / 409

Précision et sensibilité des systèmes asservis Précision des systèmes asservis

Ordre des erreurs

Erreur d’ordre n
Ecart en régime permanent entre Pconsigne et sortie mesur ée en réponse à une entrée
polynômiale de degré n : r (t) = ni=0 ai t i Γ(t), R(s) = s(n+1)
N

Exemples
1
I Erreur de position : erreur d’ordre 0 (erreur statique) : r (t) = Γ(t), R(s) = s
1
I Erreur de traı̂nage (ou de vitesse) : erreur d’ordre 1 : r (t) = tΓ(t), R(s) = s2
I Erreur d’accélération : erreur d’ordre 2 : r (t) = 12 t 2 Γ(t), R(s) = 1
s3

Erreur relative
limt→∞ (t)
I r (t) = E0 Γ(t) : erreur de position relative : E0
limt→∞ (t)
I r (t) = V0 tΓ(t) : erreur de vitesse relative : V0
limt→∞ (t)
I r (t) = W0 t 2 Γ(t) : erreur d’accélération relative : W0

(UDS) 247 / 409

Précision et sensibilité des systèmes asservis Précision des systèmes asservis

Classe d’un système

Un pôle en 0 dans un système réalise l’intégration de l’entrée du système.

système G(s) d’entrée U(s) et de

sortie Y (s)
1
G(s) =
s
U(s)
Y (s) =
s
Z t
y(t) = u(t)dt
0

Conséquence
La sortie d’un système contenant un intégrateur devient constante lorsque l’entrée du
système devient nulle

Classe d’un système continu asservi

Nombre d’intégrateurs (pôles en 0) dans la boucle ouverte

(UDS) 248 / 409

Précision et sensibilité des systèmes asservis Précision des systèmes asservis

Erreur d’ordre n
On considère un système asservi stable en boucle fermée (FBF )

B(s)
FBO (s) =
A(s)sl
avec l classe du système

N(s)
R(s) =
sn+1
avec n l’ordre de la consigne
R(s) N(s)A(s)sl
(s) = = n+1
1 + FBO (s) s (B(s) + sl A(s))
les zéros de B(s) + sl A(s) sont à partie réelle négative (puisqu’on a supposé que le
système stable en BF)
Th. de la valeur finale : (t) converge à condition que l ≥ n. Dans ce cas,

lim (t) = lim s(s)

t→∞ s→0

N(0)A(0)sl−n
= lim
s→0 B(0) + sl A(0)

(UDS) 249 / 409

Précision et sensibilité des systèmes asservis Précision des systèmes asservis

Bilan

I (t) → 0 ssi l > n

I Si l = n l’erreur tend vers une valeur finie
I Si l < n l’erreur est infinie (erreur non bornée)
B(0)
Cas courants (on pose K = A(0)
, gain de Bode de la BO)
E0 V0 W0
Entrée échelon R(s) = s rampe R(s) = s2
parabole R(s) = s3
E0
classe 0 1+K
∞ ∞
V0
classe 1 0 K
∞
W0
classe 2 0 0 K

I Lorsque l = n, l’erreur diminue quand on augmente le gain du système. Mais on

tend aussi en général à déstabiliser le système.
I =⇒ Compromis à trouver entre précision et stabilité

(UDS) 250 / 409

Précision et sensibilité des systèmes asservis Précision des systèmes asservis

Rejet des perturbations

Erreur d’ordre n par rapport à une perturbation

Valeur en régime permanent de l’écart entre consigne et sortie mesurée lorsque la
consigne est nulle et une entrée de perturbation d’ordre n activée

P (s )

r (t ) +
(t ) + y (t )
G1(s) G2(s)
+
−
ym (t )

H (s)

n ordre de la perturbation, l classe de la BO, l1 nombre d’intégrateurs dans G1 (s) et

l2 = l − l1 nombre d’intégrateurs dans G2 H.

(UDS) 251 / 409

Précision et sensibilité des systèmes asservis Précision des systèmes asservis

Rejet des perturbations (suite)

Num
P(s) = et R(s) = 0
sn+1
B(s)
FBO (s) =
A(s)sl
N(s)
G2 (s)H(s) =
D(s)sl2

−G2 (s)H(s)P(s) −N(s)A(s)sl1 Num

(s) = = n+1
1 + FBO (s) s D(s)(A(s)sl + B(s))
(t) converge à condition que l1 ≥ n. Dans ce cas,

lim (t) = lim s(s)

t→∞ s→0

−N(0)A(0)sl1 Num(0)
= lim
s→0 sn D(0)(A(0)sl + B(0))

(UDS) 252 / 409

Précision et sensibilité des systèmes asservis Précision des systèmes asservis

H (s ) H (s)
Ym (s) Ym (s)

T (s)

(UDS) 257 / 409

Précision et sensibilité des systèmes asservis Précision des systèmes asservis

Effet de la chaı̂ne de retour

Si le système asservi est de classe l, le préfiltre permet d’annuler l’erreur d’ordre l − 1
entre r et y en régime permanent à condition que la dynamique du capteur ne soit pas
excitée par la consigne

20(s + 5)2 4000

G(s) = et H(s) =
s2 (s + 10) s + 1000

(UDS) 258 / 409

Précision et sensibilité des systèmes asservis Précision des systèmes asservis

Importance de la chaı̂ne de retour

Perturbations
L’effet des perturbations sur la chaı̂ne de retour (perturbations de mesure) ne peut pas
être rejeté de l’erreur y (t) − r (t), à moins de pouvoir mesurer ces perturbations

20 1
G(s) = avec H(s) = 1 et V (s) = (mais inconnue)
s(s + 10) s

R (s) (s) Y (s)

G(s)
+
-
Ym (s)
+
H (s)
+

V (s)

(UDS) 259 / 409

Précision et sensibilité des systèmes asservis Précision des systèmes asservis

Conclusion

I Pour annuler une erreur permanente d’ordre n par rapport à la consigne, il faut
n + 1 intégrateurs dans la chaı̂ne directe
I Pour rejeter en régime statique l’effet d’une perturbation d’entrée ou de sortie
d’ordre n, il faut n + 1 intégrateurs entre le comparateur et le point d’application de
la perturbation (on dit n + 1 intégrateurs avant le point d’application de la
perturbation)
I Pour compenser le comportement statique non unitaire de la chaı̂ne de retour
H(s), il faut utiliser un préfiltre de gain H(0)
I Les perturbations de mesure ne peuvent pas être rejetées à moins de pouvoir les
mesurer

(UDS) 260 / 409

Précision et sensibilité des systèmes asservis Sensibilité des systèmes asservis

Plan du cours
Systèmes linéaires et transformée de Laplace
Quelques exemples de systèmes
Modélisation des systèmes
Propriétés des systèmes
Transformée de Laplace
Fonction de transfert
Pôles et zéros des systèmes linéaires
Réponse temporelle des systèmes linéaires
Représentation des systèmes : schémas fonctionnels et graphes de fluence
Schémas fonctionnels
Graphes de fluence
Réponse fréquentielle des systèmes linéaires
Fonction de transfert harmonique
Diagrammes de Bode
Diagrammes de Nyquist
Diagrammes de Black
Stabilité des systèmes
Stabilité des systèmes asservis
Critère de Nyquist
Lieu d’Evans
Précision et sensibilité des systèmes asservis
Précision des systèmes asservis
´
(UDS)
` 261 / 409
Précision et sensibilité des systèmes asservis Sensibilité des systèmes asservis

Sensibilité d’un modèle

Modèle T (α) complexe dépendant d’un paramètre α (complexe) T (α) = |T (α)| ejϕ(α)
définition
Sensibilité de T par rapport à α
dT
d ln T dT α
SαT = = T
dα
=
d ln α α
dα T

définition
Sensibilité du module de T par rapport à α
d|T |
d ln |T | |T | d |T | α
Sα|T | = = dα
=
d ln α α
dα |T |

définition
Sensibilité de l’argument de T par rapport à α
dϕ
d ln ϕ ϕ dϕ α
Sαϕ = = dα
=
d ln α α
dα ϕ

(UDS) 262 / 409

Précision et sensibilité des systèmes asservis Sensibilité des systèmes asservis

Sensibilité d’un système

I Sensibilité d’une fonction de transfert

F (s) fonction de transfert nominale
Dépend d’un certain nombre de paramètres (ex : éléments électronique,
température, etc.)
dF
d ln F dF α
SαF = = F
dα
=
d ln α α
dα F
I Sensibilité d’une réponse fréquentielle
G(jω) réponse harmonique nominale
Sensibilité en fonction de la pulsation
dG(jω)
d ln G(jω) G(jω) dG(jω) α
SαG(jω) = = dα
=
d ln α α
dα G(jω)

(UDS) 263 / 409

Précision et sensibilité des systèmes asservis Sensibilité des systèmes asservis

Sensibilité des boucles ouvertes et fermées

I Sensibilité par rapport à une variation de G(s)
En BO :
Y (s) = G(s)R(s)
Si R(s) ne varie pas
dY dG
=
Y G
En BF :
Y (s) G
T (s) = =
R(s) 1 + GH
dG(1 + GH) − dGGH dG
dT = =
(1 + GH)2 (1 + GH)2
dT
dG 1 + GH G 1
SGT = T
dG
= =
G
(1 + GH)2 G dG 1 + GH
Si R(s) ne varie pas
dY 1 dG
=
Y (1 + GH) G
La sensibilité est diminuée d’un facteur 1 + GH

(UDS) 264 / 409

Précision et sensibilité des systèmes asservis Sensibilité des systèmes asservis

Sensibilité d’une boucle fermée

I Sensibilité par rapport à une variation de H(s)

Y (s) G
T (s) = =
R(s) 1 + GH

G2 dH
dT = −
(1 + GH)2
dT
G2 dH 1 + GH H GH
SHT = T
dG
=− =−
G
(1 + GH)2 G dH 1 + GH
Si R(s) ne varie pas et si |GH| 1
GH
SHY = − ∼ −1
(1 + GH)

dY dH
∼−
Y H
La sensibilité par rapport à la chaı̂ne de retour est très importante

(UDS) 265 / 409

Précision et sensibilité des systèmes asservis Sensibilité des systèmes asservis

Sensibilité d’une boucle fermée

Z1
−
A
+
Vs
Ve

(UDS) 266 / 409

Précision et sensibilité des systèmes asservis Sensibilité des systèmes asservis

A
Vs = Z1
Ve
1+ Z 1+Z 2
A
A(Z 1+Z 2) Z2
Si A 1 Vs = V
AZ 1+Z 1+Z 2 e
∼ (1 + Z1
)Ve

Ve Vs
A
+
-

Z1
β = Z 1+Z 2
Variation de A
dVs 1 dA dA
= ∼
Vs 1 + βA A βA2
Variation de β
dVs βA dβ dβ
=− ∼
Vs 1 + βA β β

(UDS) 267 / 409

Précision et sensibilité des systèmes asservis Sensibilité des systèmes asservis

Conclusion

I Les BF sont peu sensibles aux variations de la chaı̂ne directe

I Il faut prendre garde aux dérives sur la chaı̂ne de retour et choisir les composants
du capteur avec précaution

(UDS) 268 / 409

Introduction

Plan du cours

Introduction

Correction série
Principe de la compensation pôles - zéros
Stabilité interne
Principe de la correction parfaite
Correcteurs série classiques
Réglage à l’aide des diagrammes harmoniques : synthèse fréquentielle
Réglage à partir du lieu d’Evans : synthèse par placement de pôles
Réglage empirique des correcteurs
Réalisation des correcteurs

Correction parallèle

Correction par boucles imbriquées

Conclusion

(UDS) 269 / 409

Introduction

Correction dynamique des systèmes

I Correction : modification du procédé physique par ajout d’un (ou de plusieurs)

bloc(s) appelé(s) correcteur ou régulateur et rétroaction, afin de donner au
système complet un comportement aussi proche que possible d’un comportement
désiré
I Dynamique : le correcteur est un système dynamique

(UDS) 270 / 409

Introduction

Objectifs de la synthèse des correcteurs numériques

Cahier des charges

I Stabilité : système en BF stable
I Performance :
I suivi de consigne (précision)
I rapidité
I comportement selon un modèle (type 2ème ordre)
I rejet des perturbations
I Robustesse : Conserver stabilité et performances malgré :
I incertitudes sur les paramètres
I dynamiques non modélisées
I non-linéarités

Rappel : le choix du correcteur dépend du système à asservir

(UDS) 271 / 409

Correction série

Plan du cours

Introduction

Correction parallèle

Correction par boucles imbriquées

Conclusion

(UDS) 272 / 409

Correction série

Principe de la correction série

w (t ) δ y (t )

+ + + y (t )
r (t ) (t ) u (t )
C (s ) G(s )
+ +
-

ym (t ) +
H (s )
+

v(t)

I ym signal de mesure
I G(s) procédé ou processus
I écart d’asservissement
I H(s) capteur
I u commande
I C(s) correcteur
I w(t) perturbation d’entrée
I r signal de consigne
I v (t) bruit de mesure
I y signal de sortie (à réguler)
I δy (t) perturbation de sortie
(UDS) 273 / 409
Correction série Principe de la compensation pôles - zéros

Plan du cours

Introduction

Correction parallèle

Correction par boucles imbriquées

Conclusion

(UDS) 274 / 409

Correction série Principe de la compensation pôles - zéros

Principe de la compensation pôle - zéro

Positionnement d’un zéro sur un pôle (ou d’un pôle sur un zéro)
Le pôle (ou le zéro) compensé n’apparaı̂t plus dans la FTBO et l’ordre du système est
artificiellement diminué de 1
Cs

Ce
Rs oscilloscope
V Re
sonde Cc
Vmes

(UDS) 275 / 409

Correction série Principe de la compensation pôles - zéros

Compensation pôle - zéro

1
Cs (s + Rs Cs
)
Vmes = Rs +Re
V
(Cs + Cc + Ce )(s + Re Rs (Cs +Cc +Ce )
)

I Si Rs Cs = Re (Cc + Ce ) :
compensation : le pôle et le zéro sont
égaux
Re
Vmes = V
Re + Rs
I Si Rs Cs > Re (Cc + Ce ) :
sur-compensation, le zéro est plus
lent que le pôle
I Si Rs Cs < Re (Cc + Ce ) :
sous-compensation : le zéro est plus
rapide que le pôle

(UDS) 276 / 409

Correction série Stabilité interne

Plan du cours

Introduction

Correction parallèle

Correction par boucles imbriquées

Conclusion

(UDS) 277 / 409

Correction série Stabilité interne

Stabilité interne des systèmes

I Stabilité BIBO insuffisante pour garantir le fonctionnement d’un système en boucle

fermée

P (s )

r (t ) +
(t ) + y (t )
C (s) G(s)
+
−

−10(s + 10)(s − 10) −10

C(s) = et G(s) =
s(s + 100) (s + 10)(s − 10)
Y (s) 100 100
= 2 =
R(s) s + 100s + 100 (s + 98.99)(s + 1.01)
Système BIBO stable

(UDS) 278 / 409

Correction série Stabilité interne

Stabilité interne
Entrée de perturbation P
Y (s) −10s(s + 100)
=
P(s) (s + 10)(s − 10)(s + 98.99)(s + 1.01)

Système non BIBO stable

I Le moindre bruit (même d’approximation numérique) rend le système instable ! !

avec le couple pôle/zéro en 10 sans le couple

(UDS) 279 / 409

Correction série Stabilité interne

Stabilité interne

On utilise une définition plus complète de la stabilité pour les systèmes en boucle
fermée : on parle de stabilité interne
Définition
Un système est stable de manière interne si toutes les fonctions de transfert entre les
entrées externes et tous les points du système ont leurs pôles à partie réelle
strictement négative
Comme toutes les fonctions de transfert ont le même dénominateur, il suffit d’étudier la
BIBO stabilité du transfert consigne / sortie avant toute simplification pôle - zéro.
Conséquences
I On ne peut pas compenser des pôles à partie réelle positive ou nulle
I On ne peut pas compenser des zéros à partie réelle positive ou nulle

(UDS) 280 / 409

Correction série Stabilité interne

Compensation des pôles et zéros instables

La position des pôles et zéros d’un système n’est jamais parfaitement connue
I La fonction de transfert est un modèle mathématique toujours incomplet
I Dynamiques non modélisées
I Non-linéarités non-modélisées
I Variations des paramètres du système (température, vieillissement, etc)
=⇒ On ne peut jamais parfaitement compenser un pôle ou un zéro
Conséquences
I Il ne faut pas compenser les pôles et zéros à partie réelle positive ou nulle
I Il faut en général éviter de compenser des pôles ou zéros complexes trop proches
de l’axe imaginaire

(UDS) 281 / 409

Correction série Stabilité interne

Compensation des pôles et zéros instables

(UDS) 282 / 409

Correction série Principe de la correction parfaite

Plan du cours

Introduction

Correction parallèle

Correction par boucles imbriquées

Conclusion

(UDS) 283 / 409

Correction série Principe de la correction parfaite

Principe de la correction
m 0
parfaite
K Πj=1 (1+τj s)
Procédé : G(s) = sl Πni=1 (1+τi s)
de classe l (0 ou 1)
I Compensation de tous les pôles pi et de tous les zéros zj
I Ajout d’un intégrateur si G est de classe 0
I Choix d’un gain

Kc Πni=1 (1 + τi s)
C(s) = 0
s1−l Πmj=1 (1 + τj s)

KKc
FTBO = C(s)G(s) =
s
1
FTBF =
1 + KKs c

Propriétés de la BF
I Système du 1er ordre que l’on peut rendre aussi rapide que désiré en augmentant
Kc
I Annulation de l’erreur statique
I Rejet des perturbations de sortie constantes
I Rejet des perturbations d’entrée constantes si l = 0

(UDS) 284 / 409

Correction série Principe de la correction parfaite

Exemple de correction parfaite

r (t ) (t ) u (t ) y (t )
C (s) G(s)
+
−

10(s + 20)
G(s) =
(s + 10)(s + 100)
100(s + 10)(s + 100)
C(s) =
s(s + 20)

I procédé initial G : stable, erreur

statique : 0.8, temps d’établissement
0.3s
I FTBF : stable, erreur statique nulle,
constante de temps 0.001s, temps
d’établissement 3ms

(UDS) 285 / 409

Correction série Principe de la correction parfaite

Limitations de la correction parfaite

I Pour être réalisable physiquement, le correcteur doit avoir au moins autant de

pôles que de zéros
I On ne peut pas compenser les pôles et zéros à partie réelle positive ou nulle
I La vitesse réelle de la boucle fermée est déterminée par les limitations de
l’actionnement du système
I On cherche en général à limiter l’ordre du correcteur afin de rendre sa réalisation
physique plus simple

(UDS) 286 / 409

Correction série Principe de la correction parfaite

Système propre

Système propre
I Système ayant au moins autant de pôles que de zéros
I Les systèmes physiques réels sont propres
I Systèmes impropres : bande passante infinie

Ex : correcteur
n m
Kc Π m
j=1 (s − zj )
X Ai X
C(s) = n
= + An + Ai si−n
Πi=1 (s − pi ) (s − pi )
i=0 i=n+1

Les éléments à partir du terme n + 1 sont des différentiateurs purs

Problèmes liés aux différentiateurs purs
I Non réalisables physiquement
I Amplifient les bruits haute fréquence
I Créent des commandes à énergie très élevée non transmissibles par le procédé
(écrêtage (saturation des étages de puissance), butées, etc.) et souvent néfastes

(UDS) 287 / 409

Correction série Principe de la correction parfaite

Effet néfaste des différentiations pures

U(s) = C(s)(s)

I Réponse à un échelon unité : u(t) = δ(t)

I Réponse à du bruit très haute fréquence (t) = η sin ωt
avec η = 1e−6 V , ω = 2πf , f = 1MHz u(t) = ηω cos ωt d’amplitude 2π ! !

(UDS) 288 / 409

Correction série Principe de la correction parfaite

Effet des saturations de l’actionnement

En pratique
saturation
r (t ) (t ) u (t ) y (t )
C (s) G(s)
+
−

(UDS) 289 / 409

Correction série Principe de la correction parfaite

Effet des saturations de l’actionnement

Saturation de l’actionnement entre -10 et +10

FTBF : constante de temps 0.001s mais temps d’établissement 0.25s !

En saturant la commande, on allonge le temps d’établissement ! ! (effet de la saturation
du terme intégral)
On ne peut pas rendre un système physique aussi rapide que l’on souhaite

(UDS) 290 / 409

Correction série Principe de la correction parfaite

Forme d’un correcteur série

Kc Πmj=1 (s − zj )
C(s) =
Πni=1 (s − pi )

I Ordre du correcteur : degré du dénominateur, nombre de pôles

I Il faut que n ≥ m : le correcteur doit être un système propre pour être réalisable
physiquement et on veut filtrer les bruits haute fréquence
I Les pôles et zéros d’un correcteur peuvent être des paires de complexes
conjugués
I Exemple : un correcteur d’ordre 2 possède 2 pôles, aucun, un ou deux zéros et un
gain Kc . Un correcteur d’ordre n a donc 2n + 1 paramètres réglables.

(UDS) 291 / 409

Correction série Correcteurs série classiques

Plan du cours

Introduction

Correction parallèle

Correction par boucles imbriquées

Conclusion

(UDS) 292 / 409

Correction série Correcteurs série classiques

Correcteurs série classiques

Correcteur proportionnel
I Le plus simple C(s) = Kc , ordre 0
I Commande u(t) = Kc (t) proportionnelle à l’écart
I Utilisé lorsque le procédé contient déja un intégrateur et que le compromis
stabilité / rapidité de la boucle fermée est satisfaisant
I Permet une augmentation de la bande passante mais couplée à une diminution de
la marge de phase

(UDS) 293 / 409

Correction série Correcteurs série classiques

Correcteur proportionnel
5
Exemple : G(s) = s(s+10)
et H(s) = 1

Cahier des charges

I Marge de phase ≥ 45o
I Bande passante ≥ 5rad/s
I Erreur statique nulle

(UDS) 294 / 409

Correction série Correcteurs série classiques

Correcteur proportionnel
5
Exemple : G(s) = s(s+10)
et H(s) = 1
Kc = 28

Cahier des charges

I Marge de phase ≥ 45o
I Bande passante ≥ 5rad/s BP = 10 rad/s et ϕM = 45o
I Erreur statique nulle

(UDS) 294 / 409

Correction série Correcteurs série classiques

Correcteurs série classiques

Correcteur proportionnel - dérivé idéal (PD idéal)

S’il y a beaucoup de bruit, il faut prendre τd proche de τc

I Correcteur d’ordre 2
I On remplace la dérivation pure par une dérivation approchée
(1+τc1 s)(1+τc2 s)
I C(s) = Kc (1+τ s)(1+τ s)
avec τd1 < τd2 < τc1 < τc2
d1 d2

(UDS) 300 / 409

Correction série Correcteurs série classiques

Correcteur PDD 2

Représentation

R + U
+
D
+
- +

I PDD 2 idéal : P = Kc , D = Kc (τc1 + τc2 )s et D 2 = Kc τc1 τc2 s2

(UDS) 301 / 409

Correction série Correcteurs série classiques

Correcteur PDD 2
(1+0.1s)2
Comportement fréquentiel C(s) = (1+0.01s)2
(Kc = 1)

(UDS) 302 / 409

Correction série Correcteurs série classiques

Correcteurs série classiques

Correcteur proportionnel - intégral idéal (PI idéal)

I C(s) = Kc 1+τs c s = Kc τc + Kc
s
correcteur d’ordre 1
Kc
Rt
I U(s) = Kc τc (s) + (s) =⇒ u(t) = Kc τc (t) + Kc
s 0
(t)dt Commande
proportionelle à l’écart et à l’intégrale de l’écart
I Utilisé lorsque le procédé ne contient pas d’intégrateur, afin de garantir une erreur
statique nulle, et que le compromis rapidité / stabilité est satisfaisant

(UDS) 303 / 409

Correction série Correcteurs série classiques

Correcteur PI

Représentation

P
R + U
+
- +
I

Kc
I PI idéal : P = Kc τc et I = s

(UDS) 304 / 409

Correction série Correcteurs série classiques

Correcteur PI

1+0.1s
Comportement fréquentiel C(s) = s
(Kc = 1)

(UDS) 305 / 409

Correction série Correcteurs série classiques

Correcteur PI

K
Procédé + capteur : GH(s) = Πni=1 (1+τi s)

Réglage standard
I Compensation du pôle dominant : τc = τ1 (τ1 constante de temps associée au
pôle stable dominant) ou
I τc = τ1 + τ2 + · · · somme des constantes de temps des pôles dominants
I Réglage de Kc pour obtenir la marge de phase requise

(UDS) 306 / 409

Correction série Correcteurs série classiques

Correcteur PI
500
Exemple : G(s) = (s+10)(s+100)
et H(s) = 1

Cahier des charges

I Marge de phase ≥ 45o
I Bande passante ≥ 100rad/s
I Erreur statique nulle

On compense le pôle en −10 et on ajoute un

Correcteur proportionnel - intégral - dérivé réel

I On remplace la dérivation pure par une dérivation réelle
c1 s)(1+τc2 s)
I C(s) = Kc (1+τs(1+τ s)
avec τd τc1 < τc2
d
I Utilisé lorsque le procédé ne contient pas d’intégrateur, afin de garantir une erreur
statique nulle, et pour augmenter la BP à marge de phase constante ou
augmenter la marge de phase à BP constante
I Inutile pour les systèmes d’ordre 1 (un PI est alors le correcteur parfait)

(UDS) 308 / 409

Correction série Correcteurs série classiques

Correcteurs série classiques

Correcteur proportionnel - intégral - dérivé généralisé (PID généralisé)

(1+2ζc τc s+τc2 s2 )
I C(s) = Kc s(1+τd s)
I Permet de compenser des pôles complexes conjugés du procédé
I Utilisé lorsque le procédé n’a pas d’intégrateur et possède une paire de pôles
complexes conjugués dominants

(UDS) 309 / 409

Correction série Correcteurs série classiques

Correcteur PID

Représentation

R + U
+
I
+
- +

Kc
I PID idéal : P = Kc (τc1 + τc2 ), I = s
et D = Kc τc1 τc2
I PID réel : P = Kc (τc1 + τc2 − τd ), I = Kc
s
et D = Kc (τc1 −τ1+τ
d )(τc2 −τd )s
s d

(UDS) 310 / 409

Correction série Correcteurs série classiques

Correcteur PID

(1+0.1s)(1+0.2s)
Comportement fréquentiel C(s) = s(1+0.001s)
(Kc = 1)

(UDS) 311 / 409

Correction série Correcteurs série classiques

I Marge de phase ≥ 45o
I Bande passante ≥ 20rad/s
I Erreur statique nulle

370(s+10)(s+20)
C(s) = s(s+200)
BP = 24rad/s et ϕM = 45o

(UDS) 313 / 409

Correction série Correcteurs série classiques

Limitations des correcteurs classiques

I Saturation du terme intégral

Procédé ayant un intégrateur + Correcteur ayant un intégrateur + Consigne élevée
Régime permanent : y (t) = cste, r (t) = cste −→ (t) = 0 et u(t) = 0 =⇒ Egalité
des surfaces entre y et r et entre r et y
5
Ex : G(s) = s(s+10) et C(s) = 8(s+2)
s
avec commande saturée à ±5V

saturation
r (t ) (t ) u (t ) y (t )
C (s) G(s)
+
−

(UDS) 314 / 409

Correction série Correcteurs série classiques

Limitations des correcteurs classiques

I Dérivation de la consigne
Même avec une dérivation réelle, la dérivation des échelons conduit à des
commandes à énergie très importante
5
Exemple : G(s) = s(s+10) et C(s) = 1000(s+10)
(s+100)

=⇒ on préfère parfois dériver uniquement la mesure

(UDS) 315 / 409

Correction série Correcteurs série classiques

Limitations des correcteurs classiques

R + U
+
I
+
- -

900s
P(s) = 100 et D(s) = s+100

(UDS) 316 / 409

Correction série Correcteurs série classiques

Limitations des correcteurs classiques

Dérivation de la mesure seulement

I Comportement identique en régime permanent
I Comportement identique vis-à-vis des perturbations
I Comportement dynamique différent
I =⇒ A utiliser uniquement pour des consignes de type échelon ou créneau

(UDS) 317 / 409

Correction série Correcteurs série classiques

Choix de la structure du correcteur

Procédé de classe 0
I PI pour annuler l’erreur statique. Permet également de rejeter les perturbations
d’entrée constantes
I PID si la stabilité n’est pas suffisante pour la bande passante désirée
I Correcteurs plus complexes (ordre supérieur) dans les cas plus compliqués

Procédé de classe 1 ou supérieure

I P
I PD pour améliorer la stabilité à la bande passante désirée
I PDD 2 si la stabilité n’est toujours pas suffisante pour la bande passante désirée
I PID si on veut rejeter une perturbation d’entrée constante (rarement suffisant pour
un procédé de classe 1)
I Correcteurs plus complexes (ordre supérieur) dans les cas plus compliqués

(UDS) 318 / 409

Correction série Correcteurs série classiques

Réglage standard des paramètres

I τc du PI ou PD : constante de temps dominante du (procédé + capteur) ou somme

des constantes en l’absence de pôle dominant
I τc1 et τc2 du PID : les 2 plus grandes constantes de temps du (procédé + capteur),
ou la paire de pôles conjugués dominants (PID généralisé, attention si les pôles
sont très oscillants)
τc1 τc1
I τd du PD réel et PID réel : de 10
si le système est bruité à 100
si le système est
peu bruité
I Kc : A déterminer en fonction de la marge de phase souhaitée
Réglage approximatif : K10ci où Kci est le gain du correcteur amenant l’instabilité

(UDS) 319 / 409

Correction série Correcteurs série classiques

Réglages plus fins

I Réglages standards non optimaux

I Réglages standards non utilisables pour des systèmes instables en BO
I Réglages standards ne permettent pas toujours de respecter un cahier des
charges spécifiques
I Nécessitent un modèle complet du système
I =⇒ Nécessité de méthodes de réglage plus fines

(UDS) 320 / 409

Correction série Réglage à l’aide des diagrammes harmoniques : synthèse fréquentielle

Plan du cours

Introduction

Correction parallèle

Correction par boucles imbriquées

Conclusion

(UDS) 321 / 409

Correction série Réglage à l’aide des diagrammes harmoniques : synthèse fréquentielle

Réglage à partir des diagrammes harmoniques

I Diagrammes de Nyquist
I Diagrammes de Black
I Diagrammes de Bode

Réglage par les diagrammes de Bode

I Cahier des charges donné sous forme de marge de phase (ou de gain), de bande
passante et de précision
I Attention ! : Valable uniquement pour les systèmes dont la BO ne contient pas de
pôle à partie réelle strictement positive (vérifie les conditions du critère du revers) :
dans le cas contraire les définitions des marges de stabilité ne sont plus valables

(UDS) 322 / 409

Correction série Réglage à l’aide des diagrammes harmoniques : synthèse fréquentielle

Utilisation des diagrammes de Bode pour régler un correcteur

I Réglage du gain des PI, PD et PID : voir précédemment

I Réglage optimal des correcteurs (PI, avance de phase)
I Réglage des correcteurs sans fonction de transfert modèle

(UDS) 323 / 409

Correction série Réglage à l’aide des diagrammes harmoniques : synthèse fréquentielle

Réglage optimal d’un correcteur PI

Données : diagramme de Bode de la boucle ouverte non corrigée (GH)

Cahier des charges
I Bande passante > ω0
I Marge de phase > ϕM
I Erreur statique nulle
I + rejet éventuel de perturbations constantes

Motivations
I la compensation du pôle dominant n’est pas toujours adéquate (cas de plusieurs
pôles à constantes de temps équivalentes)
I Il est souvent nécessaire de placer le zéro du correcteur aussi loin que possible
de zéro pour éviter une compensation partielle de l’effet de l’intégrateur

(UDS) 324 / 409

Correction série Réglage à l’aide des diagrammes harmoniques : synthèse fréquentielle

Réglages
I Si pour ω = ω0 , ϕ(GH) > − π2 + ϕM il n’est pas nécessaire d’ajouter un zéro dans
le correcteur
Alors C(s) = Ksc et on obtient Kc par
ω0
Kc =
|GH(jω0 )|
Le correcteur a un simple effet intégral.
I Si pour ω = ω0 , ϕ(GH) < − π2 + ϕM la présence du zéro est indispensable
Le zéro doit amener δϕ = − π2 + ϕM − ϕ(GH(jω0 )) en ω = ω0 .
π
I Si δϕ > 2
, un correcteur PI est insuffisant, il faut utiliser un PID.
π
I Si δϕ < 2
tan (δϕ)
On veut que arctan(τc ω0 ) = δϕ =⇒ τc = ω0
.
On obtient enfin Kc par
ω0
Kc = q
|GH(jω0 )| 1 + τc2 ω02
I Si δϕ est proche de π2 le zéro du correcteur doit être placé très proche de zéro. Il a
alors un effet néfaste de compensation de l’intégrateur. Dans ce cas, on préfère
utiliser un correcteur PID.

(UDS) 325 / 409

Correction série Réglage à l’aide des diagrammes harmoniques : synthèse fréquentielle

Réglage optimal d’un PI

Cahier des charges
I bande passante > 10rad/s
I erreur statique nulle
I marge de phase > 45o

(UDS) 326 / 409

Correction série Réglage à l’aide des diagrammes harmoniques : synthèse fréquentielle

Réglage optimal d’un PI (exemple)

pour ω = 10rad/s, ϕ(GH) = −62.5o et GdB = −5.35dB

Le zéro doit apporter δϕ = 17.5o en ω = 10rad/s

tan (δϕ)
τc = = 0.0315s
ω0

ω0
Kc = q = 17.66
|GH(jω0 )| 1 + τc2 ω02

17.66(1 + 0.0315s)
C(s) =
s

(UDS) 327 / 409

Correction série Réglage à l’aide des diagrammes harmoniques : synthèse fréquentielle

Réglage optimal d’un PI

(UDS) 328 / 409

Correction série Réglage à l’aide des diagrammes harmoniques : synthèse fréquentielle

Réglage optimal d’un PI

Cahier des charges
I bande passante > 10rad/s
I erreur statique nulle
I marge de phase > 45o

(UDS) 329 / 409

Correction série Réglage à l’aide des diagrammes harmoniques : synthèse fréquentielle

Réglage optimal d’un PI

pour ω = 10rad/s, ϕ(BO) = −132o et GdB = −16dB
Le zéro doit apporter δϕ = 87o en ω = 10rad/s
Il faut utiliser un PID ! Voilà ce que donne le réglage d’un PI :

tan (δϕ)
τc = = 1.9068s
ω0

ω0
Kc = q = 3.3044
|GH(jω0 )| 1 + τc2 ω02

3.3044(1 + 1.9068s)
C(s) =
s

(UDS) 330 / 409

Correction série Réglage à l’aide des diagrammes harmoniques : synthèse fréquentielle

Réglage optimal d’un PI

Le zéro du PI est trop proche de zéro et l’effet de l’intégrateur est fortement perturbé : il
faut utiliser un PID

(UDS) 331 / 409

Correction série Réglage à l’aide des diagrammes harmoniques : synthèse fréquentielle

Réglage pratique d’un PI

Données : diagramme de Bode de la BO non corrigée

Cahier des charges
I Marge de phase 45o
I erreur statique nulle

Réglage
I On détermine sur le diagramme de Bode la pulsation ω0 pour laquelle ϕ = − π2
I On choisit τc = 1
ω0
ce qui amène la phase de CG à −135o pour ω = ω0
I On détermine Kc pour amener le gain à 1 en ω = ω0
ω0
Kc = √
2 |GH(ω0 )|

(UDS) 332 / 409

Correction série Réglage à l’aide des diagrammes harmoniques : synthèse fréquentielle

Réglage optimal d’un correcteur PD réel

Motivations
I Placer le zéro aussi loin que possible de zéro pour ne pas affecter l’effet intégral
I Placer le pôle à la fréquence minimale pour filtrer au mieux le bruit haute
fréquence

Caractéristiques de l’avance de phase

I Pulsation pour laquelle l’avance de phase amenée par un PD réel est maximale
1+τ cs
C(s) = Kc 1+τ s
=⇒ ϕ = arctan τc ω − arctan τd ω
d

dϕ τc τd
= −
dω 1 + τc2 ω 2 1 + τd2 ω 2
r
dϕ 1
= 0 ⇐⇒ ω =
dω τc τd
q
1
ωmax = τc τd
moyenne géométrique des pulsations associées au pôle et au
zéro : ωc = τ1c et ωd = 1
τd
I Valeur de l’avance maximale ϕmax = ϕ(ωmax )
On montre que sin(ϕmax ) = ττcc −τ
+τ
d
d
(UDS) 333 / 409
Correction série Réglage à l’aide des diagrammes harmoniques : synthèse fréquentielle

Réglage optimal d’un correcteur PD réel

1 + 0.1s
C(s) =
1 + 0.01s
20

15
Ampl. (dB)

0
−1 0 1 2 3
10 10 10 10 10
w(rad/s)

40
phase (deg)

30
ϕmax
20

0
−1 0 1 2 3
10 10 10 10 10
w(rad/s) ωmax

(UDS) 334 / 409

Correction série Réglage à l’aide des diagrammes harmoniques : synthèse fréquentielle

Réglage optimal d’un correcteur PD réel

Données : diagramme de Bode de la boucle ouverte non corrigée comportant un

intégrateur (GH)
Cahier des charges
I Bande passante > ω0
I Marge de phase > ϕM

Principe du réglage optimal

I Régler le correcteur pour que le maximum de phase soit apporté juste à la
fréquence de coupure
I Amener juste l’avance de phase nécessaire

(UDS) 335 / 409

Correction série Réglage à l’aide des diagrammes harmoniques : synthèse fréquentielle

Réglage optimal d’un correcteur PD réel

Réglage
Pour ω = ω0 , il faut amener un déphasage positif valant δϕ = −π + ϕM − ϕ(GH(ω0 ))
I si δϕ > π2 , un correcteur PD est insuffisant : il faut utiliser un correcteur PDD 2 . En
général on préfèrera utiliser un PDD 2 dès que δϕ > π3
I si δϕ < 0, un correcteur proportionnel est suffisant
I si δϕ ∈ [0, π3 ], on peut utiliser un correcteur PD.
Réglage optimal : amener le maximum du déphasage positif en ω0 afin de bien
filtrer le bruit haute fréquence
r
1 τc − τd
=⇒ ωmax = = ω0 et sin(δϕ) = = sin δϕ
τc τd τc + τd
1 τd 1 − sin(δϕ)
=⇒ τc τd = et =
ω02 τc 1 + sin(δϕ)
On obtient enfin Kc par
q
1 + τd2 ω02
Kc = q
|GH(jω0 )| 1 + τc2 ω02
I Si τc est trop grand on a tendance à compenser l’effet de l’intégrateur du procéd
(UDS) 336é.
/ 409
Correction série Réglage à l’aide des diagrammes harmoniques : synthèse fréquentielle

Réglage optimal d’un correcteur PD réel

Exemple :

Cahier des charges

I Bande passante > 10rad/s
I Marge de phase > 60o
I Erreur statique nulle

(UDS) 337 / 409

Correction série Réglage à l’aide des diagrammes harmoniques : synthèse fréquentielle

Réglage optimal d’un correcteur PD réel

ϕ(10rad/s) = −159o et |GH(10rad/s)|dB = −35dB

δϕ = 39o

1
τc τd = = 0.01
ω02
τd 1 − sin(δϕ)
= = 0.2275
τc 1 + sin(δϕ)
=⇒ τc = 0.2096 et τd = 0.0477
q
1 + τd2 ω02
Kc = q = 26.83
|GH(jω0 )| 1 + τc2 ω02

26.83(1 + 0.2096s)
C(s) =
1 + 0.0477s

(UDS) 338 / 409

Correction série Réglage à l’aide des diagrammes harmoniques : synthèse fréquentielle

(UDS) 339 / 409

Correction série Réglage à l’aide des diagrammes harmoniques : synthèse fréquentielle

(UDS) 340 / 409

Correction série Réglage à l’aide des diagrammes harmoniques : synthèse fréquentielle

Démarche générale d’une synthèse fréquentielle

1. Il est nécessaire d’exprimer le cahier des charges en des termes fréquentiels :

I les notions de rapidité doivent être exprimées sous forme de bande passante,
I les notions de stabilité, de dépassement doivent être exprimées sous forme de marge
de phase ou de marge de gain.
I On peut ajouter des contraintes de rejet de perturbations
I et de précision (gain statique, annulation d’erreurs d’ordre n par rapport à la consigne,
etc.)
2. Choix partiel de la structure du correcteur : détermination des intégrateurs à
ajouter dans la boucle ouverte en fonction des contraintes de précision
C(s) = s1n C 0 (s)
3. On trace les diagrammes de Bode en amplitude et en phase de l’ensemble
G(s)H(s) s1n , soit à partir d’une fonction de transfert modèle de GH, soit à partir du
diagramme de Bode mesuré du système
4. Choix de la structure complète du correcteur
5. Réglage des pôles et zéros (choix arbitraires ou calculs optimaux) en fonction des
contraintes de rapidité et de marges de stabilité
6. Détermination du gain du correcteur

(UDS) 341 / 409

Correction série Réglage à l’aide des diagrammes harmoniques : synthèse fréquentielle

Démarche générale d’une synthèse fréquentielle

Remarques
I Méthode de synthèse utilisable pour les systèmes à retard
I A ne pas utiliser si la boucle ouverte contient des pôles à partie réelle positive
I le cahier des charges peut être mal posé et ne pas avoir de solution (cas des
systèmes à retard)

(UDS) 342 / 409

Correction série Réglage à l’aide des diagrammes harmoniques : synthèse fréquentielle

Synthèse fréquentielle : exemple

Diagramme de Bode mesuré de GH (boucle ouverte non corrigée)

cahier des charges

I marge de phase> 60o
I bande passante > 100 rad/s
I Erreur statique nulle
I rejet des pertubations de
sortie constantes

(UDS) 343 / 409

Correction série Réglage à l’aide des diagrammes harmoniques : synthèse fréquentielle

I Contraintes de précision =⇒ il faut un intégrateur dans C(s)

(UDS) 344 / 409

Correction série Réglage à l’aide des diagrammes harmoniques : synthèse fréquentielle

I La phase de GHs
en ωo vaut −256o . =⇒ le zéro d’un PI est insuffisant pour
remonter la phase en −120o =⇒ il faut un PID (PI + PD = PI + avance de phase)
I On place le zéro du PI pour amener suffisament de phase en ωo , mais pas trop
près de 0 pour ne pas compenser l’effet de l’intégrateur : par exemple en s = −10

(UDS) 345 / 409

Correction série Réglage à l’aide des diagrammes harmoniques : synthèse fréquentielle

GH(s+10)
I La phase de s
en ωo vaut −172o . =⇒ il faut amener 52o avec l’avance de
phase
1+0.029s
I Réglage optimal du PD =⇒ 1+0.0034s

(UDS) 346 / 409

Correction série Réglage à l’aide des diagrammes harmoniques : synthèse fréquentielle

GH(s+10)(1+0.029s)
I Le gain de s(1+0.0034s)
vaut −2.25dB en ωo
2.25 11.1(s+34.48)(s+10)
I Réglage de Kc = 10 20 = 1.30 d’où C(s) = s(s+294.1)

(UDS) 347 / 409

Correction série Réglage à partir du lieu d’Evans : synthèse par placement de pôles

Plan du cours

Introduction

Correction parallèle

Correction par boucles imbriquées

Conclusion

(UDS) 348 / 409

Correction série Réglage à partir du lieu d’Evans : synthèse par placement de pôles

Réglage à partir du lieu d’Evans

I Cahier des charges donné sous forme de temps de réponse ou de bande

passante, de dépassement, et de précision
I Attention ! : Valable uniquement pour les systèmes sans retard dont la fonction de
transfert peut être mise sous forme de fraction rationnelle
I Nécessite une fonction de transfert modèle du système contrairement aux
réglages par méthode fréquentielle
I Outils : lieu d’Evans accompagné des abaques ζ constant et ωn constant
I Méthode : placer les pôles de la boucle fermée de sorte à obtenir le
comportement demandé
I Point fort : applicable aux systèmes ayant des pôles ou des zéros à partie réelle
positive en boucle ouverte

(UDS) 349 / 409

Correction série Réglage à partir du lieu d’Evans : synthèse par placement de pôles

Réglage standard d’un PI

I Compensation du pôle dominant

I Ajout d’un intégrateur
I Réglage du gain pour amener les pôles de la BF à la position désirée
20
exemple : G(s) = (s+5)(s+10)

Cahier des charges

I Dépassement < 5%
I temps d’établissement < 1s
I Annulation de l’erreur statique
I Rejet des perturbations de sortie constantes

(UDS) 350 / 409

Correction série Réglage à partir du lieu d’Evans : synthèse par placement de pôles

(UDS) 351 / 409

Correction série Réglage à partir du lieu d’Evans : synthèse par placement de pôles

2.5(s+5)
Solution : C(s) = s

(UDS) 351 / 409

Correction série Réglage à partir du lieu d’Evans : synthèse par placement de pôles

Réglage standard d’un PD

I Compensation du pôle dominant

I Ajout d’un pôle 10 fois plus rapide
I Réglage du gain pour amener les pôles de la BF à la position désirée
20
exemple : G(s) = (s+5)(s+10)s

Cahier des charges

I Dépassement < 5%
I temps d’établissement < 1s
I Annulation de l’erreur statique
I Rejet des perturbations de sortie constantes

(UDS) 352 / 409

Correction série Réglage à partir du lieu d’Evans : synthèse par placement de pôles

(UDS) 353 / 409

Correction série Réglage à partir du lieu d’Evans : synthèse par placement de pôles

103(s+5)
Solution : C(s) = s+50

(UDS) 353 / 409

Correction série Réglage à partir du lieu d’Evans : synthèse par placement de pôles

Réglage standard d’un PID

I Compensation des pôles dominants

I Ajout d’un pôle 10 fois plus rapide que le plus rapide des pôles compensés
I Ajout d’un intégrateur
I Réglage du gain pour amener les pôles de la BF à la position désirée
50
exemple : G(s) = (s+5)(s+10)(s+15)

Cahier des charges

I Dépassement < 5%
I temps d’établissement < 1s
I Annulation de l’erreur statique
I Rejet des perturbations de sortie constantes

(UDS) 354 / 409

Correction série Réglage à partir du lieu d’Evans : synthèse par placement de pôles

(UDS) 355 / 409

Correction série Réglage à partir du lieu d’Evans : synthèse par placement de pôles

194(s+5)(s+10)
Solution : C(s) = (s+100)s

(UDS) 355 / 409

I Il existe un lien exact entre marge de phase et dépassement (voir abaque)

I rappel : pour un système de cette forme, la BF est toujours stable et ϕM ∈]0, 90[

(UDS) 360 / 409

Correction série Réglage à partir du lieu d’Evans : synthèse par placement de pôles

Effet des zéros

Les zéros de la chaı̂ne directe et donc du correcteur sont aussi des zéros de la boucle
fermée
N N
C(s) : DC , G(s) : DG , H = NDH
C G H

NC NG
BF =
DC DG DH + NC NG NH
Les zéros de la BF (donc les zéros de CG) modifient le comportement de la BF par
rapport à une BF ayant les mêmes pôles mais pas de zéros (voir chap 2)
I les zéros à partie réelle positive ralentissent le système et le rendent à déphasage
non-minimal =⇒ à éviter dans le correcteur
I les zéros à partie réelle plus petite que les pôles dominants tendent à créer des
dépassements (phénomène de sur-compensation) =⇒ à éviter
I les zéros à partie réelle négative plus grande que les pôles dominants tendent à
accélérer le système =⇒ peu gênants
Lorsqu’un zéro de la BF (donc de CG) est proche d’un pôle de la boucle ouverte, il a
tendance à compenser son effet (compensation approchée) =⇒ Eviter les zéros trop
proches de 0 qui compensent l’effet des intégrateurs

(UDS) 361 / 409

Correction série Réglage à partir du lieu d’Evans : synthèse par placement de pôles

A propos des contraintes sur le lieu des racines

Attention ! Les abaques d’iso-amortissement et d’iso-pulsation ainsi que les contraintes
de dépassement et de temps de réponse ne sont valables que pour des systèmes
d’ordre 2 sans zéro ou équivalents (autres pôles dominés et zéros à effet négligeable)
Exemples
100K (s+50)
I BO(s) = s(s+10)

Les abaques et les contraintes (t5% < 0.8s et D < 5%) sont valables

(UDS) 362 / 409

Correction série Réglage à partir du lieu d’Evans : synthèse par placement de pôles

10K (s+10)
I BO(s) = s(s2 +10s+41)

Les abaques et les contraintes (t5% < 1s et D < 5%) ne sont pas valables : BF de
type 3ème ordre

(UDS) 363 / 409

Correction série Réglage à partir du lieu d’Evans : synthèse par placement de pôles

10K (s2 +6.94s+34.3)

I BO(s) = s(s2 +8.95s+68)

Les abaques et les contraintes (t5% < 1s et D < 5%) ne sont pas valables : effet
des zéros non négligeable

(UDS) 364 / 409

Correction série Réglage à partir du lieu d’Evans : synthèse par placement de pôles

Exemple de réglage par le lieu d’Evans

100
G(s) = s(s+10)

Cahier des charges

I Temps d’établissement à 5% < 0.2s
I Erreur statique nulle
I Erreur permanente par rapport à une perturbation d’entrée constante nulle
I Dépassement < 25%

Contraintes de précision =⇒ il faut un intégrateur dans C(s)

I PI : pour obtenir le temps d’établissement, il faut que le pôle dominant de la BF
soit < −15 : impossible à obtenir avec un PI
I PID : pour simplifier le système on compense le pôle du procédé et on le remplace
par un pôle 10 fois plus rapide : C(s) = K (s+10)(s−z)
s(s+100)
Réglage de la position du
zéro : en le rapprochant de l’axe imaginaire on stabilise le système mais en le
ralentissant
compromis possible : C(s) = 36.4(s+13)(s+10)
s(s+100)

(UDS) 365 / 409

Correction série Réglage à partir du lieu d’Evans : synthèse par placement de pôles

Remarque : l’approximation du 2nd ordre n’est pas valable. Le dépassement est

principalement du au zéro lent en −13.

(UDS) 366 / 409

Correction série Réglage à partir du lieu d’Evans : synthèse par placement de pôles

Bilan

I Il faut toujours commencer par les correcteurs les plus simples possibles
I Le réglage du gain du correcteur se fait dans la très grande majorité des cas en
dernier, car :
I le gain le modifie pas le diagramme de phase
I le gain ne modifie pas le lieu des racines
I il faut toujours utiliser le nombre minimal d’intégrateur dans le correcteur car les
intégrateurs déstabilisent les BF

(UDS) 367 / 409

Correction série Réglage empirique des correcteurs

Plan du cours

Introduction

Correction parallèle

Correction par boucles imbriquées

Conclusion

(UDS) 368 / 409

Correction série Réglage empirique des correcteurs

Réglage empirique des correcteurs

I Méthodes de réglage pratique

I Ne nécessitent pas de modèle du système
I Réglage directement sur le système physique
I Permet un compromis acceptable entre stabilité et bande passante

Méthodes de Ziegler-Nichols
I Recherche des paramètres d’un PI ou PID

1
I PI : C(s) = K 1+
Ti s
Td s
I PID : C(s) = K 1 + T1s + 1+0.1Td s
i

I Ne permet pas de répondre à un cahier des charges spécifique

(UDS) 369 / 409

Correction série Réglage empirique des correcteurs

1ère Méthode de Ziegler - Nichols

Réglage à partir de la réponse indicielle en BO
I Approximation de la réponse temporelle par un système d’ordre 1 + retard
I Mesure de la réponse indicielle
I Mesure de la pente λ de la tangente au point d’inflexion
I Mesure du temps mort Lr

E0
0.9E0 Lr
I PI : K = λLr
, Ti = 0.3
1.2E0
I PID : K = λLr
, Ti = 2Lr , Td = 0.5Lr
λ

t
Lr

I Nécessite de pouvoir faire fonctionner le système en boucle ouverte (stable en

BO, ou instabilité non destructrice)
I Pas utilisable pour des réponses de type ordre 1 : Lr = 0 ! !
(UDS) 370 / 409
Correction série Réglage empirique des correcteurs

2ème Méthode de Ziegler - Nichols

Réglage par génération d’oscillation en BF

I Système en boucle fermée. On annule les effets intégral et dérivateur (correcteur
proportionnel de gain K )
I On augmente K jusqu’à obtenir une oscillation entretenue en réponse à un
échelon de faible amplitude
I On mesure le gain Ko correspondant et la fréquence d’oscillation fo
I On règle alors (valeurs empiriques)
0.85
I PI : K = 0.45Ko et Ti = fo
1 1
I PID : K = 0.6Ko , Ti = 2fo
, Td = 8fo

Nécessite de pouvoir faire osciller la boucle fermée

I impossible par exemple pour les systèmes d’ordre 1 et 2 à phase minimale en BO
I limitations pratiques : délicat de faire osciller par exemple un processus chimique

(UDS) 371 / 409

Correction série Réglage empirique des correcteurs

2ème Méthode de Ziegler - Nichols

70
G(s) = (s+10)(s+30)(s+50)

Ki ∼ 2400 et f0 ∼ 7.6

(UDS) 372 / 409

Correction série Réglage empirique des correcteurs

2ème Méthode de Ziegler - Nichols

1080(s+8.99)
I PI C(s) = s

(UDS) 373 / 409

Correction série Réglage empirique des correcteurs

15840(s2 +57s+849)
I PID C(s) = s(s+611)

(UDS) 374 / 409

Correction série Réglage empirique des correcteurs

Auto-réglage

D’autres méthodes plus complexes permettent de régler des caractéristiques de la

boucle fermée ou de la boucle ouverte, par exemple la marge de phase souhaitée
Exemple : méthode du Relais (non linéaire)

(UDS) 375 / 409

Correction série Réalisation des correcteurs

Plan du cours

Introduction

Correction parallèle

Correction par boucles imbriquées

Conclusion

(UDS) 376 / 409

Correction série Réalisation des correcteurs

Implémentation électronique des correcteurs

Correcteur Proportionnel

Ve −
R1 Vs

R2
Vs = (1 + )Ve
R1

(UDS) 377 / 409

Correction série Réalisation des correcteurs

Implémentation électronique des correcteurs

Correcteur Proportionnel - Intégral

Ve −
R1 Vs

R2
C2

1
R1C2
(1 + (R1 + R2)C2s)
Vs = Ve
s

(UDS) 378 / 409

Correction série Réalisation des correcteurs

Implémentation électronique des correcteurs

Correcteur Proportionnel - Dérivé réel

+
R1
Ve −
Vs
C1
R2

(1 + R2
R1
)(1 + R1R2(C1+C2)
R1+R2
)
Vs = Ve
1 + R2C2s

(UDS) 379 / 409

Correction parallèle

Plan du cours

Introduction

Correction parallèle

Correction par boucles imbriquées

Conclusion

(UDS) 380 / 409

Correction parallèle

R (s) (s) U (s ) Y (s)

+ G2(s) G3(s)
G1 ( s )
+
- -

1
G1 CH 2
G3H

ω
FTBO

(UDS) 385 / 409

Correction parallèle

Exemple : la correction tachymétrique

R (s) (s) U (s) Y (s)

+ Kv 1
K 1+τ s s
+
- -

vitesse
C (s)
Ω

position
Ym (s)

Cahier des charges

I Erreur statique nulle
I bande passante > ωco
I marge de phase > 45o

(UDS) 386 / 409

Correction parallèle

Exemple

Kv 1
FTBO = K
1 + τ s + Kv C(s) s

I La boucle ouverte corrigée contient un intégrateur =⇒ l’erreur statique sera nulle

I C(s) permet de modifier la constante de temps de la boucle ouverte, donc la
phase à la pulsation ωco .
I le gain K permet de modifier le gain de la boucle ouverte donc la bande passante

(UDS) 387 / 409

Correction parallèle

Réglage graphique asymptotique

Kv
G2 (s) = 1+τ s
et H2 (s) = 1.
I G2 CH2 1
KKv
FTBO =
s(1 + τ s)
I G2 CH2 1
K
FTBO =
sC(s)

(UDS) 388 / 409

Correction parallèle

Réglage graphique asymptotique

||dB

1
τ
ωco ω
Kv
s(1+τ s)

FTBO désirée

−135

√ 2
ωco 1+τ 2 ωco
Il faut augmenter le gain K pour obtenir la bande passante désirée : K = Kv

(UDS) 388 / 409

Correction parallèle

Réglage graphique asymptotique

R (s) (s) U (s) Y (s)

+ Kv 1
K s(1+τ s) s
+
- -

C (s)

position
Ym (s)

Cahier des charges

I Erreur statique nulle
I bande passante > ωco
I marge de phase > 45o

(UDS) 390 / 409

Correction parallèle

Exemple

KKv
FTBO =
s(1 + τ s) + Kv C(s)

I La boucle ouverte corrigée doit contenir un intégrateur =⇒ C(s) doit contenir s en

facteur à son numérateur
I C(s) permet de modifier la constante de temps de la boucle ouverte, donc la
phase à la pulsation ωco .
I le gain K permet de modifier le gain de la boucle ouverte donc la bande passante

(UDS) 391 / 409

Correction parallèle

Réglage graphique asymptotique

I G2 CH2 1
KKv
FTBO =
s(1 + τ s)
I G2 CH2 1
K
FTBO =
C(s)

(UDS) 392 / 409

Correction parallèle

Réglage graphique asymptotique

||dB

1
τ
ωco ω
Kv
s(1+τ s)

FTBO désirée

−135

√ 2
ωco 1+τ 2 ωco
Il faut augmenter le gain K pour obtenir la bande passante désirée : K = Kv

(UDS) 392 / 409

Correction parallèle

Réglage graphique asymptotique

||dB

1
τ
ωco ω
KKv On déduit une forme désirée
s(1+τ s) K
simple de sC(s) qui permet
FTBO désirée d’obtenir la forme désirée de
FTBO
ϕ K
C(s)
= ωsco
ω

−135

(UDS) 392 / 409

Correction parallèle

Réglage graphique asymptotique

||dB

1
τ
ωco ω
KKv
s(1+τ s)
K
sC (s)
D’où on tire : C(s) = ωKsco
FTBO désirée Il faut dériver la sortie,
c’est-à-dire estimer la
ϕ vitesse de rotation

−135

(UDS) 392 / 409

Correction parallèle

Réglage exact

A partir de l’analyse graphique on choisit C(s) = Kc s.

On a alors :
Kv
(1+τ s)s KKv
FTBO = K Kv Kcs
=
1 + s(1+τ s)
s(1 + Kv Kc + τ s)

On veut que ϕ(ωco ) = −135o =⇒

π ωco τ 3π
− − arctan ( )=−
2 1 + Kv Kc 4

=⇒ Kc = ωcoKτv−1
Kc permet de régler la phase à ωco
On veut que |FTBO(ωco )| = 1 =⇒
KKv
p =1
ωco (1 + Kv Kc )2 + (τ ωco )2 )
√ √
ω (1+Kv Kc )2 +(τ ωco )2 ) ω2 τ 2
=⇒ K = co Kv
= coKv
K permet de régler la bande passante

(UDS) 393 / 409

Correction parallèle

Correction tachymétrique

Kc s
En pratique on utilise C(s) = 1+Ts
pour que la dérivation soit réalisable et pour filtrer le
bruit haute fréquence

(UDS) 394 / 409

Correction par boucles imbriquées

Plan du cours

Introduction

Correction parallèle

Correction par boucles imbriquées

Conclusion

(UDS) 395 / 409

Correction par boucles imbriquées

régulateur secondaire
2(s)
R (s ) 1(s) Y (s )
+ G2(s) G1(s)
C1(s) C2(s)
+
- -

régulateur H2(s)
primaire boucle secondaire (interne)

H1(s)
Ym (s) boucle principale (externe)

Asservissement de signaux intermédiaires (ils doivent pouvoir être mesurés ou

estimés)
Exemple : moteur asservi en position (boucle principale), avec asservissement de
vitesse (boucle secondaire) et asservissement de courant (boucle interne)

(UDS) 396 / 409

Correction par boucles imbriquées

Intérêts
I Rejet plus efficace et plus rapide des perturbations sur la boucle interne
I Utilisation de correcteurs plus simples
I Permet souvent d’éviter l’utilisation de termes dérivés dans le régulateur principal
I Amélioration du régime transitoire lorsque la boucle interne est beaucoup plus
rapide que la boucle principale
I Linéarisation d’une partie non linéaire de la chaı̂ne

(UDS) 397 / 409

Correction par boucles imbriquées

Principe de réglage

I Réglage du régulateur de la boucle interne. On cherche à :

I rejeter les perturbations sur la partie G2 du système
I à rendre la boucle interne plus rapide que G1
Usuellement un PI ou PID qui compense les pôles de G2 (on cherche à se
rapprocher du correcteur parfait pour G2 de sorte à rendre la boucle interne de
type 1er ou 2ème ordre)
C2 G2
I On calcule la fonction de transfert de la boucle interne : 1+C2 G2
C2 G2
I Réglage du correcteur C1 sur la boucle ouverte principale : G .
1+C2 G2 1
On cherche
à
I Assurer une erreur statique nulle
I Rejeter les perturbations constantes sur G1
I accélerer le système : compenser les constantes de temps lentes (celles de G1 )
Usuellement un PI ou PID

(UDS) 398 / 409

Correction par boucles imbriquées

Exemple : Moteur à courant continu avec frottements secs
Commande : tension d’induit U, Sortie : vitesse de rotation Ω
On mesure également le courant d’induit I
Cahier des charges
I Rejeter l’effet de la force contre-électromotrice
I Rejeter l’effet des frottements secs
I temps de réponse < 10 ms
I Dépassement < 5%

-
U I Γ Ω
Gelec Ki Gmeca
+ + +

5 50
Gmeca = s+10
Gelec = s+100
Ki = 1

(UDS) 399 / 409

Correction par boucles imbriquées

correction du courant
I correcteur PI (correcteur parfait)
I on veut rendre la boucle de courant très rapide par rapport à la partie mécanique
du système
20(s+100)
I CI (s) = s

-
I∗ U I Γ Ω
CI Gelec Ki Gmeca
+ + +
- +

Calcul de la boucle de courant : I = 1

1+0.001s
I∗

(UDS) 400 / 409

Correction par boucles imbriquées

5000
FTBO =
(s + 10)(s + 1000)

correction de la vitesse
I Rejet des frottements secs =⇒ nécessite un intégrateur
I Réglage d’un PI par le lieu d’Evans
100(s+10)
I CΩ (c) = s

-
Ω∗ I∗ U I Γ Ω
CΩ CI Gelec Ki Gmeca
+ + + +
- +

(UDS) 401 / 409

Correction par boucles imbriquées

Comparaison avec un correcteur série simple

Cahier des charges

I Rejeter l’effet de la force contre-électromotrice
I Rejeter l’effet des frottements secs
I temps de réponse < 10 ms
I Dépassement < 5%
250
En considérant la fcem comme une perturbation G(s) = Gelec Gmeca = (s+10)(s+100)
1000(s+10)(s+100)
Correction avec PID (correcteur parfait) C(s) = s(s+1000)

-
Ω∗ U I Γ Ω
C Gelec Ki Gmeca
+ + +
- +

(UDS) 402 / 409

Correction par boucles imbriquées

Comparatif
Cahier des charges identique
I Rejet de la fcem (Sans frottements secs)

Boucles imbriquées Boucle simple

La boucle simple rejette mal la fcem, car la partie électrique n’est pas
suffisamment rapide pour considérer la perturbation comme constante. Au
contraire, la boucle de courant a linéarisé la partie électrique en rejetant la fcem

(UDS) 403 / 409

Correction par boucles imbriquées

En général la tension d’induit est saturée

-
Ω∗ U Γ Ω
I∗ I
CΩ CI Gelec Ki Gmeca
+ + + +
- +

-
Ω∗ U Γ Ω
I
C Gelec Ki Gmeca
+ + +
- +

(UDS) 404 / 409

Correction par boucles imbriquées

Comparatif

I Comportement du couple avec saturation de la tension

Boucles imbriquées Boucle simple

La boucle de courant permet une montée de la vitesse à courant (couple) constant

(UDS) 405 / 409

Conclusion

Plan du cours

Introduction

Correction parallèle

Correction par boucles imbriquées

Conclusion

(UDS) 406 / 409

[3] Philippe de Larminat.

Automatique appliquée.
Hermès, 2007.

[4] G. F. Franklin, J.D. Powell, and A. Emami-Naeini.

Feedback control of dynamic systems.
Prentice Hall, 2002.

[5] K. Ogata.
Modern control Engineering.
Prentice Hall, 2002.

[6] Eric Ostertag.

Systèmes et asservissements continus.
Technosup - Ellipses, 2004.

(UDS) 409 / 409

Vous aimerez peut-être aussi

Communication en Automatique 2
Pas encore d'évaluation
Communication en Automatique 2
95 pages
Cours sur les Systèmes Asservis
100% (1)
Cours sur les Systèmes Asservis
43 pages
Introduction à l'automatique et asservissement
100% (1)
Introduction à l'automatique et asservissement
111 pages
Asservissement des Systèmes Linéaires
100% (2)
Asservissement des Systèmes Linéaires
190 pages
Asservissement et régulation en automatique
Pas encore d'évaluation
Asservissement et régulation en automatique
28 pages
Introduction aux Systèmes Asservis Linéaires
Pas encore d'évaluation
Introduction aux Systèmes Asservis Linéaires
7 pages
Systèmes Linéaires en Automatique
100% (2)
Systèmes Linéaires en Automatique
99 pages
Cours Auto Complet
Pas encore d'évaluation
Cours Auto Complet
212 pages
Cours d'Automatique et Régulation
100% (3)
Cours d'Automatique et Régulation
82 pages
Automatique Des Systèmes Continus
Pas encore d'évaluation
Automatique Des Systèmes Continus
83 pages
Asservissement et régulation des systèmes
Pas encore d'évaluation
Asservissement et régulation des systèmes
28 pages
Système de Commande Automatique
Pas encore d'évaluation
Système de Commande Automatique
55 pages
Chap 1
Pas encore d'évaluation
Chap 1
28 pages
Cours sur les systèmes asservis en électrotechnique
Pas encore d'évaluation
Cours sur les systèmes asservis en électrotechnique
53 pages
Systèmes Linéaires en Automatique
100% (2)
Systèmes Linéaires en Automatique
99 pages
Asservissement numérique-ISGE - 074124
Pas encore d'évaluation
Asservissement numérique-ISGE - 074124
175 pages
Concepts des asservissements linéaires
Pas encore d'évaluation
Concepts des asservissements linéaires
81 pages
Notions sur les Systèmes Asservis
Pas encore d'évaluation
Notions sur les Systèmes Asservis
41 pages
Théorie des Systèmes Asservis Linéaires
Pas encore d'évaluation
Théorie des Systèmes Asservis Linéaires
115 pages
Asservissement des Systèmes Linéaires
100% (1)
Asservissement des Systèmes Linéaires
53 pages
Chapitre1-2025-2026 (Généralités Sur Les Sytèmes Asservis)
Pas encore d'évaluation
Chapitre1-2025-2026 (Généralités Sur Les Sytèmes Asservis)
7 pages
Introduction à l'Automatique des Systèmes
Pas encore d'évaluation
Introduction à l'Automatique des Systèmes
121 pages
Modélisation et Commande des Systèmes
Pas encore d'évaluation
Modélisation et Commande des Systèmes
18 pages
Génie Automatique : Systèmes Linéaires et Asservis
100% (2)
Génie Automatique : Systèmes Linéaires et Asservis
62 pages
Introduction à l'automatique et systèmes dynamiques
Pas encore d'évaluation
Introduction à l'automatique et systèmes dynamiques
37 pages
Systèmes Asservis en Électrotechnique
Pas encore d'évaluation
Systèmes Asservis en Électrotechnique
71 pages
Systèmes Asservis : Cours et Exercices
100% (2)
Systèmes Asservis : Cours et Exercices
97 pages
Systèmes Asservis Linéaires Continus
Pas encore d'évaluation
Systèmes Asservis Linéaires Continus
69 pages
Automatique des Systèmes Linéaires
Pas encore d'évaluation
Automatique des Systèmes Linéaires
94 pages
Types de systèmes asservis et critères
Pas encore d'évaluation
Types de systèmes asservis et critères
20 pages
Cours sur l'automatique continue
Pas encore d'évaluation
Cours sur l'automatique continue
74 pages
Cours sur l'asservissement des systèmes
Pas encore d'évaluation
Cours sur l'asservissement des systèmes
53 pages
Cours sur l'asservissement des systèmes
Pas encore d'évaluation
Cours sur l'asservissement des systèmes
28 pages
Performances des Systèmes Asservis
Pas encore d'évaluation
Performances des Systèmes Asservis
67 pages
Systèmes Asservis Linéaires Continus
Pas encore d'évaluation
Systèmes Asservis Linéaires Continus
38 pages
Introduction aux systèmes asservis
Pas encore d'évaluation
Introduction aux systèmes asservis
9 pages
Systèmes Linéaires Continus et Invariants
Pas encore d'évaluation
Systèmes Linéaires Continus et Invariants
71 pages
Automatique
Pas encore d'évaluation
Automatique
41 pages
Automatique
Pas encore d'évaluation
Automatique
73 pages
Cours sur l'asservissement des systèmes
88% (8)
Cours sur l'asservissement des systèmes
56 pages
Introduction aux Systèmes Asservis
Pas encore d'évaluation
Introduction aux Systèmes Asservis
24 pages
Introduction aux Systèmes Asservis
Pas encore d'évaluation
Introduction aux Systèmes Asservis
13 pages
Notions de base en automatique et asservissements
Pas encore d'évaluation
Notions de base en automatique et asservissements
71 pages
Asservissement et Boucles de Contrôle
86% (7)
Asservissement et Boucles de Contrôle
82 pages
Introduction aux systèmes asservis
Pas encore d'évaluation
Introduction aux systèmes asservis
8 pages
Cours sur les Systèmes Linéaires Asservis
Pas encore d'évaluation
Cours sur les Systèmes Linéaires Asservis
109 pages
f2138560 Automatique
Pas encore d'évaluation
f2138560 Automatique
125 pages
Introduction à l'automatique et systèmes
Pas encore d'évaluation
Introduction à l'automatique et systèmes
21 pages
Asservissements des Systèmes Linéaires
Pas encore d'évaluation
Asservissements des Systèmes Linéaires
53 pages
Introduction aux systèmes asservis en automatique
Pas encore d'évaluation
Introduction aux systèmes asservis en automatique
108 pages
Systèmes Asservis en Automatique
Pas encore d'évaluation
Systèmes Asservis en Automatique
36 pages
Transformée de Laplace et Systèmes Linéaires
Pas encore d'évaluation
Transformée de Laplace et Systèmes Linéaires
59 pages
Introduction aux systèmes asservis
Pas encore d'évaluation
Introduction aux systèmes asservis
54 pages
Introduction aux systèmes asservis
Pas encore d'évaluation
Introduction aux systèmes asservis
24 pages
Cours sur l'Automatique Linéaire
Pas encore d'évaluation
Cours sur l'Automatique Linéaire
68 pages
Généralités Automatique
Pas encore d'évaluation
Généralités Automatique
8 pages
Cours d'Automatique - Support 2012
100% (1)
Cours d'Automatique - Support 2012
59 pages
Cours d'Automatique à l'ENIT
Pas encore d'évaluation
Cours d'Automatique à l'ENIT
59 pages
Cours de Regulation Et Asservissement
Pas encore d'évaluation
Cours de Regulation Et Asservissement
31 pages
Exercices de révision en chimie
Pas encore d'évaluation
Exercices de révision en chimie
2 pages
Évaluation Diagnostique Physique-Chimie SIBM
Pas encore d'évaluation
Évaluation Diagnostique Physique-Chimie SIBM
4 pages
Évaluation de l'état mécanique en physique
Pas encore d'évaluation
Évaluation de l'état mécanique en physique
2 pages
Modélisation des matériaux ferromagnétiques
Pas encore d'évaluation
Modélisation des matériaux ferromagnétiques
111 pages
Estérification et Hydrolyse : Mécanismes et Contrôle
Pas encore d'évaluation
Estérification et Hydrolyse : Mécanismes et Contrôle
47 pages
Grandeurs physiques et unités SI
Pas encore d'évaluation
Grandeurs physiques et unités SI
1 page
Questions de base en physique-chimie
Pas encore d'évaluation
Questions de base en physique-chimie
4 pages
Comprendre le courant de Foucault
100% (2)
Comprendre le courant de Foucault
17 pages
Cours de Physique Quantique
100% (2)
Cours de Physique Quantique
149 pages
Exercices corrigés sur le cyclotron
Pas encore d'évaluation
Exercices corrigés sur le cyclotron
22 pages
L'expérience d'Oersted et les courants de Foucault
Pas encore d'évaluation
L'expérience d'Oersted et les courants de Foucault
2 pages
Cours de Mécanique Quantique
100% (1)
Cours de Mécanique Quantique
149 pages
Freinage par Courants de Foucault
Pas encore d'évaluation
Freinage par Courants de Foucault
6 pages
Cours de Physique Quantique
100% (1)
Cours de Physique Quantique
149 pages
Cours de Physique Quantique - Chapitre 2
100% (3)
Cours de Physique Quantique - Chapitre 2
33 pages
Cours de Physique Quantique
100% (1)
Cours de Physique Quantique
149 pages
Cours de Physique Quantique - Yamina Achkar
100% (1)
Cours de Physique Quantique - Yamina Achkar
149 pages
TP Physique Important
100% (1)
TP Physique Important
286 pages
Estimations et Échantillonnage en Statistiques
Pas encore d'évaluation
Estimations et Échantillonnage en Statistiques
2 pages
Introduction à l'Apprentissage Automatique
Pas encore d'évaluation
Introduction à l'Apprentissage Automatique
166 pages
Cours d'Algorithmique avec Exercices
Pas encore d'évaluation
Cours d'Algorithmique avec Exercices
120 pages
Commande par retour d'état expliquée
100% (1)
Commande par retour d'état expliquée
7 pages
TD - Chapitre 4 - Correcteur Numerique
Pas encore d'évaluation
TD - Chapitre 4 - Correcteur Numerique
3 pages
Introduction à la Mécanique des Milieux Continus
Pas encore d'évaluation
Introduction à la Mécanique des Milieux Continus
46 pages
Méthodes de Classification des Données
Pas encore d'évaluation
Méthodes de Classification des Données
40 pages
Correction TD Analyse Numérique S4 MIP/PC
Pas encore d'évaluation
Correction TD Analyse Numérique S4 MIP/PC
3 pages
Dérivation Numérique et Intégration
Pas encore d'évaluation
Dérivation Numérique et Intégration
7 pages
Entropie et probabilités en théorie de l'information
Pas encore d'évaluation
Entropie et probabilités en théorie de l'information
8 pages
Machine Learning en Génie Mécanique 2024
Pas encore d'évaluation
Machine Learning en Génie Mécanique 2024
18 pages
Distributions de Gibbs en mécanique statistique
Pas encore d'évaluation
Distributions de Gibbs en mécanique statistique
8 pages
Arbre de Décision pour Classification
Pas encore d'évaluation
Arbre de Décision pour Classification
47 pages
Équilibre de Nash en stratégie mixte
Pas encore d'évaluation
Équilibre de Nash en stratégie mixte
2 pages
Méthode du Simplexe en Recherche Opérationnelle
Pas encore d'évaluation
Méthode du Simplexe en Recherche Opérationnelle
2 pages
Exercices d'algorithmes en MIAGE
Pas encore d'évaluation
Exercices d'algorithmes en MIAGE
2 pages
Représentation d'état des systèmes dynamiques
0% (1)
Représentation d'état des systèmes dynamiques
4 pages
Exposés sur l'Apprentissage Automatique
Pas encore d'évaluation
Exposés sur l'Apprentissage Automatique
1 page
R01 Exercice - Solution
Pas encore d'évaluation
R01 Exercice - Solution
4 pages
Méthodes de Gradient en Analyse Numérique
75% (4)
Méthodes de Gradient en Analyse Numérique
4 pages
Canal de Transmission et Modèles Discrets
Pas encore d'évaluation
Canal de Transmission et Modèles Discrets
8 pages
Système Chaotique et Circuit de Chua
Pas encore d'évaluation
Système Chaotique et Circuit de Chua
84 pages
Annales d'examens corrigés d'Algèbre 2
Pas encore d'évaluation
Annales d'examens corrigés d'Algèbre 2
101 pages
Devoir Surveillé Maths : Suites Terminale ES/L
Pas encore d'évaluation
Devoir Surveillé Maths : Suites Terminale ES/L
2 pages
Comprendre la récursivité en C
Pas encore d'évaluation
Comprendre la récursivité en C
17 pages
Analyse temporelle et stabilité des systèmes
Pas encore d'évaluation
Analyse temporelle et stabilité des systèmes
3 pages
Introduction au Machine Learning et ses usages
Pas encore d'évaluation
Introduction au Machine Learning et ses usages
4 pages
Événements et Probabilités Élémentaires
Pas encore d'évaluation
Événements et Probabilités Élémentaires
80 pages
Probabilités et Échantillonnage en Statistique
Pas encore d'évaluation
Probabilités et Échantillonnage en Statistique
5 pages
Analyse de systèmes de contrôle et régulation
Pas encore d'évaluation
Analyse de systèmes de contrôle et régulation
3 pages