automatique cours et exercices, Exercises of Automatic Controls

Download automatique cours et exercices and more Exercises Automatic Controls in PDF only on Docsity! Automatique Linéaire 1 1A ISMIN Automatique linéaire 1 – J.M. Dutertre – 2016 Automatique Linéaire 1 1A ISMIN www.emse.fr/~dutertre/enseignement.html - 2015 2 Sommaire. I. Introduction, définitions, position du problème. p. 3 I.1. Introduction. p. 3 I.2. Définitions. p. 5 I.3. Position du problème. p. 6 II. Modélisation des systèmes linéaires. p. 11 II.1. Système du premier ordre. p. 12 II.2. Système du second ordre. p. 21 II.3. Systèmes d’ordre supérieur à 2. p. 31 III. Stabilité des systèmes asservis. p. 33 III.1. Schéma général d’un asservissement. p. 33 III.2. Interprétation géométrique du passage de la boucle ouverte à la boucle fermée. p. 36 III.3. Réponse impulsionnelle d’un système bouclé en régime linéaire. p. 40 III.4. Le critère de Routh-Hurwitz (critère algébrique). p. 41 III.5. Les critères géométriques de stabilité. p. 44 IV. Performances des systèmes asservis. p. 50 IV.1. Précision. p. 50 IV.2. Rapidité des systèmes. p. 56 V. Correction des systèmes asservis. p. 59 V.1. Introduction. p. 59 V.2. Correction proportionnelle et dérivée (P.D.) – Correction à avance de phase. p. 63 V.3. Correction proportionnelle et intégrale (P.I.) – Correction à retard de phase. p. 66 V.4. Correction proportionnelle intégrale et dérivée (P.I.D.). p. 68 V.5. Modèle du second ordre. p. 72 Bibliographie. p. 74 Annexe 1 – Signaux type. p. 75 Annexe 2 – Transformée de Laplace. p. 77 Annexe 3 – Systèmes linéaires du second ordre. p. 80 Annexe 4 – Abaque de Black-Nichols. p. 84 Automatique Linéaire 1 1A ISMIN www.emse.fr/~dutertre/enseignement.html - 2015 5 I.2. Définitions. Cette partie rappelle, ou donne, un certain nombre de définitions permettant d’aborder rigoureusement la suite du cours. Définition 1 : On appelle modèle d’un système (ou processus) la loi qui relie l’entrée (cause) à la sortie (effet). Définition 2 : On distingue deux régimes dans le comportement des systèmes : • le régime permanent ou établi, caractérisant la réponse stabilisée du système à une entrée quelconque, • le régime transitoire, caractérisant l’évolution de la réponse avant que le régime permanent ne soit atteint. Le régime statique est le régime permanent dans le cas ou l’entrée est constante. Définition 3 : Un système est causal si sa sortie y(t) à un instant t0 ne dépend que des valeurs de son entrée u(t) pour t ≤ t0 (cf. figure I.3). Fig. I.3 – Système. Un système causal ne répond pas avant d’être excité (système non anticipatif). Les systèmes physiques temporels réalisables sont causals. Un signal x(t) est causal si ∀t < 0 x(t) = 0. En pratique un signal temporel est toujours causal, à condition de bien choisir l’origine des temps. Définition 4 : Un système à temps invariant a un modèle identique à tout instant (un retard τ ne change pas la loi du modèle) : € u(t) système⎯ → ⎯ ⎯ y(t) u(t −τ) système⎯ → ⎯ ⎯ y(t −τ) Automatique Linéaire 1 1A ISMIN www.emse.fr/~dutertre/enseignement.html - 2015 6 Définition 5 : Un système est dit instantané si à un instant donné sa sortie ne dépend que de l’excitation à cet instant : y(t) = a.u(t) Dans tous les autres cas, il est dit, à mémoire ou dynamique, par exemple pour : y(t) = a.u(t-τ) ou : y(t) = a.u(t) + b.y’(t) Définition 6 : Un système est stable si et seulement si toute entrée bornée génère une sortie bornée. Un système physique est stable s’il retourne spontanément vers son état d’équilibre lorsqu’il en est écarté. Il est instable si sa sortie n’a pas de valeur fixe (asymptotiquement) lorsque son entrée est nulle. Définition 7 : Un système est linéaire s’il satisfait au principe de superposition : € a.u1(t) + b.u2(t) syst. linéaire⎯ → ⎯ ⎯ ⎯ a.y1(t) + b.y2(t) Ce cours traite des systèmes causals, linéaires et à temps invariant : les S.L.T.I. Les systèmes étudiés sont analogiques, leurs signaux d’entrée et de sortie sont continus à la fois en temps et en amplitude. La relation qui lie leur entrée et leur sortie2 est dès lors une équation différentielle linéaire à coefficients constants. I.3. Position du problème. a – La commande automatique ou comment remplacer l’homme. La finalité de l’automatique, telle que nous venons de la définir, est de remplacer l’homme ou de suppléer à ses limites dans la conduite d’un système (cf. figure I.4 si l’on revient à l’exemple de la partie I.1 concernant une automobile). La problématique se réduit dès lors à l’étude et à la modélisation du système considéré dans le but d’élaborer une commande automatique. 2 Au singulier, on se limitera en effet à l’étude des systèmes monovariables, c’est-à-dire ayant une entrée et une sortie. Automatique Linéaire 1 1A ISMIN www.emse.fr/~dutertre/enseignement.html - 2015 7 Fig. I.4 – Remplacer l’homme. Pour se faire, l’utilisation d’une rétroaction est nécessaire. Le système est placé en boucle fermée, ce qui introduit les notions d’asservissement (lors de la poursuite d’une consigne variable) et de régulation (concernant la compensation de perturbations3 externes). Le paragraphe suivant donne l’exemple de l’asservissement d’un chauffage central individuel. b – La boucle d’asservissement. Considérons le système de chauffage central d’un logement représenté figure I.5. Fig. I.5 – Système de chauffage central. Avec : θ température intérieure, T température de l’eau chaude envoyée dans les radiateurs, θe température extérieure (considérée comme une perturbation). T est réglée par le chauffagiste pour obtenir une température de consigne donnée θc = 19°C. Cependant, le réglage est à refaire à chaque variation de θe. 3 On définira comme étant une perturbation une entrée du système imprévisible et/ou sur laquelle on ne peut agir. Automatique Linéaire 1 1A ISMIN www.emse.fr/~dutertre/enseignement.html - 2015 10 potentielles d’instabilité sont le retard à l’exécution des ordres reçus, ou pire l’existence d’une contre réaction positive (ε = yc + y). La stabilité des systèmes asservis est étudiée au chapitre III. La précision s’exprime par l’écart ε entre la consigne yc et la sortie y du système. Dans le cas d’une loi de commande proportionnelle du type u = K.ε, l’obtention d’une bonne précision nécessite d’avoir un gain élevé (en effet pour obtenir une valeur de commande u donnée, K devra être d’autant plus importante que ε sera faible). De même, une perturbation n sera d’autant plus efficacement corrigée (erreur résiduelle faible) que K sera grand. Or, on a vu qu’un grand K peut être source d’instabilité. D’où le fait (à mémoriser) que la stabilité et la précision soient des paramètres potentiellement contradictoires. La troisième qualité essentielle d’un asservissement est sa rapidité. La rapidité d’un processus peut se mesurer par son temps de réponse à un échelon de commande comme défini au IV.2. Les notions de précision et de rapidité des systèmes font l’objet du chapitre IV. D’une façon générale, la synthèse d’un asservissement résulte d’un compromis stabilité – précision – rapidité. L’automatisation des processus requiert l’utilisation d’une boucle fermée (rétroaction), celle- ci est nécessaire afin de pouvoir : - stabiliser un système instable en boucle ouverte, - compenser des perturbations externes, - compenser des incertitudes liées au système lui-même (vieillissement, imprécision du modèle, etc.). Automatique Linéaire 1 1A ISMIN www.emse.fr/~dutertre/enseignement.html - 2015 11 II. Modélisation des systèmes linéaires. La caractéristique statique (c'est-à-dire la relation entre l’entrée et la sortie en régime permanent) d’un système linéaire est une droite (cf. figure II.1). Fig. II.1 – Caractéristique statique d’un système linéaire. Cela ne doit pas amener de confusion avec le comportement dynamique du système en régime transitoire. La figure II.2 représente la réponse en sortie d’un système linéaire à un échelon u(t) sur son entrée (on constate bien que le tracer de y(t)en fonction de u(t) ne serait pas une droite). Fig. II.2 – Réponse de la sortie d’un système linéaire à un échelon en entrée. Nous avons déjà énoncé précédemment que l’équation liant la sortie et l’entrée d’un système linéaire, continu et invariant dans le temps est une équation différentielle linéaire à coefficients constants. Cette équation traduit aussi bien le comportement dynamique du système que son comportement statique (il suffit pour cela d’annuler les dérivées). Les parties suivantes traitent de la modélisation et du comportement des systèmes du premier ordre, du second ordre et d’ordre supérieur. Cependant, les systèmes physiques (réels) ne sont pas nécessairement linéaires. Il est néanmoins souvent possible de les étudier avec les outils classiques de l’automatique linéaire après avoir linéarisé leur comportement autour d’un point de repos. Cette façon de procéder est familière aux électroniciens, ils l’utilisent par exemple, pour étudier les transistors en amplification ; la figure II.3 donne l’exemple de la linéarisation autour d’un point de repos, M0(V0, I0), d’une diode. Le modèle linéaire obtenu permet d’étudier les variations de id et vd autour de M0, dès lors que ces grandeurs restent dans le domaine de validité du modèle. Automatique Linéaire 1 1A ISMIN www.emse.fr/~dutertre/enseignement.html - 2015 12 Fig. II.3 – Linéarisation d’un système autour d’un point de repos. Il y a deux façons d’obtenir le modèle linéaire d’un système : - par la mise en équation du système à partir de ces lois physiques (équations électriques, mécaniques, etc.), - ou par identification, le modèle étant déterminé expérimentalement en étudiant la réponse du système à des stimuli classiques. II.1. Système du premier ordre. Définition 8 : Un système est dit du 1er ordre si la relation entre son entrée et sa sortie est une équation différentielle du 1er ordre. Exemple : établir l’équation différentielle du circuit RC de la figure II.4. Fig. II.4 – Circuit RC. Les équations électriques du système sont : Automatique Linéaire 1 1A ISMIN www.emse.fr/~dutertre/enseignement.html - 2015 15 c – Réponse indicielle. La réponse indicielle, c'est-à-dire à un échelon (d’amplitude A0), d’un système du 1er ordre est (on note Γ(t) l’échelon unitaire5) : € u(t) = A0 .Γ(t) € y(t) = K .A0 1− e − t /τ( ) Eq. II.4 ⇓ 1 ⇑ 3 € Y (p) = K .A0 p. 1+τ.p( ) La réponse indicielle est représentée figure II.5. Fig. II.5 – Réponse indicielle d’un système du 1er ordre. La valeur finale atteinte en régime permanent par y(t) est K fois la valeur de l’entrée (K est le gain statique). Pour t = τ, y(τ) atteint 63% de la valeur finale. Le temps de réponse à 5% (le temps au bout duquel y(t) approche la valeur finale à 5% près, et y reste) est t = 3τ. La tangente à l’origine (cf. Fig. II.5) a une pente de KA0/τ, on observe effectivement une cassure assez nette de y(t) qui est caractéristique de la réponse indicielle d’un système du 1er ordre. d – Réponse à une rampe. La réponse à une rampe de pente a d’un système du premier ordre est : € u(t) = at.Γ(t) € y(t) = K .a t −τ +τ.e− t /τ( ) Eq. II.5 ⇓ 1 ⇑ 3 € U(p) = a p2 € Y (p) = K .a p2 . 1+τ.p( ) 5 Le lecteur trouvera en annexe 1 quelques rappels sur les signaux type. Automatique Linéaire 1 1A ISMIN www.emse.fr/~dutertre/enseignement.html - 2015 16 Le tracé de la réponse à une rampe est donné figure II.6. Fig. II.6 – Réponse à une rampe d’un système du 1er ordre. En régime permanent (t >> τ) on a y(t) = Ka(t - τ) : la sortie tend vers une rampe de pente : Ka. Pour K = 1, y(t) suit la rampe d’entrée avec un retard τ. La différence entre l’entrée et la sortie est appelée erreur de trainage, εt, elle vaut εt = aτ. Pour K ≠ 1, les pentes étant différentes, εt tend vers l’infini (divergence). e – Réponse harmonique. La réponse harmonique d’un système est sa réponse à une sinusoïde permanente, u(t) = Um.sin(ω.t) , le régime transitoire étant éteint. Rappel : La réponse harmonique d’un système linéaire (quel que soit son ordre) est une sinusoïde de même pulsation, d’amplitude Ym, déphasée d’un angle ϕ par rapport à l’entrée : y(t) = Ym.cos(ω.t +ϕ ) Les signaux étant périodiques, l’analyse de la réponse harmonique se fait en complexe (p = jω), et plus précisément en étudiant H(jω). Soit : € Y ( jω ) = H( jω ).U( jω ) = K 1+ jωτ .U( jω ) H( jω) = K 1+ jωτ ⎛ ⎝ ⎜ ⎞ ⎠ ⎟ qui donne : Automatique Linéaire 1 1A ISMIN www.emse.fr/~dutertre/enseignement.html - 2015 17 € Ym = Y ( jω) = K 1+ω 2τ 2 .Um (Module) € ϕ = Arg Y ( jω)( ) = −arctan ωτ( ) (Argument) L’étude des propriétés fréquentielles des systèmes linéaires (c'est-à-dire leur réponse à une action sinusoïdale permanente dont on fait varier la fréquence) permet d’en déduire les propriétés dynamiques temporelles (c’est-à-dire leur évolution dans le temps en fonction des actions subies) comme nous le verrons par la suite. C’est la raison pour laquelle on attache tant d’importance à cette étude. En général les paramètres étudiés sont le gain et le déphasage : € Gain = Ym Um = H( jω ) € Déphasage =ϕ = Arg H( jω)( ) que l’on représente sous forme de diagramme de Bode, de Black, ou de Nyquist (cf. ci-après). Diagramme de Bode. Le diagramme de Bode d’une fonction de transfert comporte deux courbes : • son module exprimé en décibels (dB), € HdB = 20.logH( jω ) • et sa phase (ou argument), € Arg H( jω)( ) tracées en fonction de la pulsation ω (axe gradué suivant une échelle logarithmique). Pour un système linéaire du 1er ordre : € H( jω ) = K 1+ jωτ On en déduit : € HdB = 20.logH( jω ) = 20.logK − 20.log 1+ω 2τ 2 € ArgH = −arctan(ωτ ) La pulsation de coupure à -3 dB, ωc, est la pulsation pour laquelle le gain exprimé en dB est inférieur de 3 dB au gain statique (gain pour ω = 0). Les tracés du gain et de la phase sont donnés figure II.7 en pointillés rouge, pour un axe des abscisses gradué en pulsation réduite ω / ωc. On trouve ωc = 1/τ. Automatique Linéaire 1 1A ISMIN www.emse.fr/~dutertre/enseignement.html - 2015 20 Fig. II.9 - Représentation de Nyquist de la réponse harmonique d’un système du 1er ordre. A la pulsation de coupure pour ω = ωc = 1/τ : € H( jω) = K / 2 € Arg H( jω)( ) = −π /4 f – Relation temps – fréquence. La rapidité de la réponse d’un système linéaire du 1er ordre est liée à sa fréquence de coupure (c’est-à-dire à sa bande passante, telle que, fc = 1/(2πτ) ; d’après ωc = 1/τ ). Le temps de montée, tm , d’un système soumis à un échelon étant le temps mis par la sortie pour passer de 10% à 90% de sa valeur finale est une façon d’exprimer cette rapidité. Or on démontre que tm = 2,2.τ On en déduit que : € tm . fc = 0,35 Eq. II.6 Ainsi plus la bande passante d’un système sera large (fc élevée) plus il sera rapide (tm faible), et inversement. Automatique Linéaire 1 1A ISMIN www.emse.fr/~dutertre/enseignement.html - 2015 21 II.2. Système du second ordre. Définition 10 : Un système est dit du second ordre si la relation entre son entrée et sa sortie est une équation différentielle du 2ème ordre. La forme générale de l’équation différentielle d’un système du deuxième ordre d’entrée u et de sortie y est (on prendra toujours un second membre indépendant de u’(t)) : € d2y(t) dt 2 + 2mω 0. dy(t) dt +ω 0 2.y(t) = Kω 0 2 .u(t) Eq. II.7 Avec : K gain statique, m coefficient d’amortissement (parfois noté ξ), ω0 pulsation propre non amortie. a – Fonction de transfert. Par application de la TL à l’équation II.7 (en prenant des conditions initiales nulles) il vient : Eq. II.8 b – Réponse indicielle. La réponse indicielle (à un échelon unitaire), d’un système du 2nd ordre est : € u(t) = Γ(t) ⇓ 1 € U(p) = 1 p € Y (p) = K p. 1+ 2m ω 0 p + p 2 ω 0 2 ⎛ ⎝ ⎜ ⎞ ⎠ ⎟ = Kω 0 2 p. p2 + 2mω 0p +ω 0 2( ) La dernière étape de détermination de y(t) nécessite l’étude de trois cas en fonction de m : • m > 1, régime apériodique. Le discriminant réduit de l’équation caractéristique du dénominateur est alors € Δ = (m2 −1)ω 0 2 > 0 Automatique Linéaire 1 1A ISMIN www.emse.fr/~dutertre/enseignement.html - 2015 22 Le dénominateur possède donc deux pôles réels distincts : € p1,2 = −mω 0 ±ω 0 m 2 −1 négatifs (car € p1 + p2 = −2mω 0 < 0 et p1.p2 =ω 0 2 > 0) tels que : € Y (p) = Kω 0 2 p. p − p1( ). p − p2( ) On cherche alors à exprimer Y(p) sous la forme : € Y (p) = α p + β p − p1 + γ p − p2 afin de calculer aisément la TL inverse de Y(p). On obtient : € Y (p) = K . 1 p + p2 p1 − p2 . 1 p − p1 + p1 p2 − p1 . 1 p − p2 ⎡ ⎣ ⎢ ⎤ ⎦ ⎥ D’où par TL-1 : € y(t) = K . 1+ p2 p1 − p2 .ep1t + p1 p2 − p1 .ep2t ⎡ ⎣ ⎢ ⎤ ⎦ ⎥ Eq. II.9 On a : € limt→+∞ y(t) = K (p1,2 pôles réels négatifs) y'(0) = 0 la tangente à l’origine est nulle Dans l’hypothèse ou m >> 1, on a € p1 − p2 = 2ω 0 m 2 −1 grand en valeur absolue. C’est-à-dire € p1 << p2 et donc € p1 p2 − p1 << p2 p1 − p2 , ainsi le terme devient très rapidement négligeable devant d’où € y(t) ≅K . 1+ p2 p1 − p2 .ep1t ⎡ ⎣ ⎢ ⎤ ⎦ ⎥ ce qui correspond à la réponse d’un système du 1er ordre. On parle alors de pôle dominant, ici p1, c’est le pôle qui a le plus d’influence sur le comportement du système (l’influence de p2 étant comparativement négligeable). Les pôles sont représentés figure II.10, le pôle dominant est le pôle le plus proche de l’axe imaginaire. Automatique Linéaire 1 1A ISMIN www.emse.fr/~dutertre/enseignement.html - 2015 25 On note Tp la pseudo période des oscillations amorties : (pour m = 0, Tp = T0 = 2π/ω0, d’où le nom de pulsation propre non amortie de ω0). Le premier dépassement correspond à : et d’une façon générale Le temps de réponse à 5%, tr5%, est donné par l’abaque de la figure II.13. Il est minimal pour € m ≅ 0,7et vaut € tr5% 2π ω 0 = 0,44 . Fig. II.13 – Temps de réponse à 5% d’un système linéaire du 2nd ordre. Pour m petit on a l’approximation (les oscillations durent longtemps) : € tr5% 2π ω 0 = 32πm Pour m < 0, un calcul similaire au cas du régime apériodique nous amène à : € y(t) = K . 1+ p2 p1 − p2 .ep1t + p1 p2 − p1 .ep2t ⎡ ⎣ ⎢ ⎤ ⎦ ⎥ Automatique Linéaire 1 1A ISMIN www.emse.fr/~dutertre/enseignement.html - 2015 26 Or Re(p1,2) = -mω0 > 0 donc € limt→+∞ e Re( p1,2 ).t = +∞ Les termes exponentiels divergent, le système est instable. La figure II.14 donne une synthèse de la réponse indicielle d’un système du second ordre en fonction de son coefficient d’amortissement m. On y retrouve les régimes apériodique, critique et pseudopériodique (stable et instable). Figure II.14 – Synthèse de la réponse indicielle d’un système du 2nd ordre en fonction de m. Automatique Linéaire 1 1A ISMIN www.emse.fr/~dutertre/enseignement.html - 2015 27 c – Réponse harmonique. La réponse harmonique d’un système du 2nd ordre est sa réponse à une sinusoïde permanente, u(t) = Um.cos(ω.t), le régime transitoire étant éteint. Partant de l’équation II.8 (p = jω) on trouve : € H( jω ) = K 1+ 2 jm ω ω 0 + j ω ω 0 ⎛ ⎝ ⎜ ⎞ ⎠ ⎟ 2 Eq. II.13 Etude asymptotique : Pour € ω →0 € H( jω ) = K d’où € HdB = 20logK € Arg H = 0 Pour € ω →+∞ d’où € Arg H = −π La figure II.15 donne le tracé asymptotique de HdB pour K = 1. La valeur de m n’a d’influence sur le tracé de HdB qu’au voisinage de ω0. Fig. II.15 – Tracé asymptotique du gain d’un système linéaire du 2nd ordre. Pour m ≥ 1, le dénominateur de H(jω) a deux racines réelles positives telles que : € H( jω ) = K 1+ jω ω1( ). 1+ jω ω 2( ) avec € ω1,2 = mω 0 ±ω 0 m 2 −1 € ω1.ω 2 =ω 0 2 Automatique Linéaire 1 1A ISMIN www.emse.fr/~dutertre/enseignement.html - 2015 30 Le tracé permet de retrouver et d’identifier un certain nombre de points caractéristiques : le régime statique (ω = 0), la résonance (pour ω = ωr on obtient MdB), la pulsation propre ω0 (elle donne QdB), la pulsation de coupure ωc à -3 dB, et l’asymptote vers -π pour € ω →∞ . Représentation de Nyquist. La figure II.19 donne la représentation de Nyquist d’un système linéaire du 2nd ordre pour différentes valeurs du coefficient d’amortissement : 0,3 < m < 5 (le gain statique est unitaire). Le tracé est obtenu à partir de l’expression : € H( jω ) = K . 1− ω ω 0( ) 2( ) 1− ω ω 0( ) 2( ) 2 + 4m2 ω ω 0( ) 2 − 2 j K .m.(ω ω 0) 1− ω ω 0( ) 2( ) 2 + 4m2 ω ω 0( ) 2 Elle-même calculée à partir de l’équation II.13. Pour € ω →+∞ on a € H( jω )→0 par les valeurs négatives (le gain tend vers 0 et la phase vers -π). Fig. II. 19 – Lieu de Nyquist d’un système linéaire du 2nd ordre (0,3 < m < 5). Automatique Linéaire 1 1A ISMIN www.emse.fr/~dutertre/enseignement.html - 2015 31 II.3. Systèmes d’ordre supérieur à 2. D’une façon générale, l’équation différentielle représentative d’un système linéaire d’ordre supérieur à 2 (S.L.T.I.) peut s’écrire : € a0.y(t) + a1. dy(t) dt + ...+ an . dny(t) dtn = b0.u(t) + b1. du(t) dt + ...+ bm . dmu(t) dtm Eq. II.14 avec ai, bi coefficients constants réels, n ≥ m , pour les systèmes physiques réalisables (c’est-à-dire respectant le principe de causalité), n est l’ordre du système. On en déduit l’expression de la fonction de transfert correspondante (conditions initiales nulles) : € H(p) = N(p) D(p) = bm p m + bm−1p m−1 + ...+ b1p + b0 an p n + an−1p n−1 + ...+ a1p + a0 Eq. II.15 Les racines du numérateur, N(p), sont les zéros de la fonction de transfert H(p), et, les racines du dénominateur, D(p), ses pôles. Les coefficients étant réels, les n pôles (p1 à pn) de H(p) sont soit réels, soit complexes conjugués deux à deux. Ainsi, le dénominateur peut s’écrire sous la forme : € D(p) = an .(p − p1).(p − p2)....(p − pn ) D’où la possibilité d’exprimer H(p) comme une somme d’éléments simples : € H(p) = Ai p − pii=1 n ∑ avec Ai complexe Soit une décomposition additive en sous-systèmes du 1er ordre (pour les pôles réels) et du 2ème ordre (pour les pôles complexes conjugués). Dès lors, la réponse du système complet est la superposition des réponses de chacun des sous- systèmes qui le composent (par application du principe de superposition), comme illustré figure II.20 par la réponse à un échelon. Automatique Linéaire 1 1A ISMIN www.emse.fr/~dutertre/enseignement.html - 2015 32 Fig. II.20 – Réponse indicielle d’un système linéaire d’ordre supérieur à 2. En termes de stabilité, il suffit d’un seul pôle à partie réelle positive pour entraîner l’instabilité de l’ensemble. Les pôles dominants sont situés à proximité de l’axe imaginaire. Pour un pôle réel cela correspond à une constante de temps élevée ; pour un couple de pôles complexes conjugués à un coefficient d’amortissement faible. La réponse globale du système dépend principalement des pôles dominants. Automatique Linéaire 1 1A ISMIN www.emse.fr/~dutertre/enseignement.html - 2015 35 D’où € FTBF = A(p) 1+ A(p).B(p) Eq. III.2 Soit : € FTBF = A(p) 1+ FTBO La figure III.3 présente le cas particulier du retour unitaire. Fig. III.3 – Cas du retour unitaire. On a alors : € FTBF = H(p) = T(p) 1+T(p) FTBFretour unitaire = FTBO 1+ FTBO ⎛ ⎝ ⎜ ⎞ ⎠ ⎟ Structure équivalente à un retour unitaire. Il est possible de ramener tout système asservi au cas du retour unitaire (cf. figure III.4) à partir de la FTBF du système réduit (le système réduit correspond au système pour lequel la sortie est Xr(p) au lieu de Y(p). Fig. III.4 – Structure équivalente à un retour unitaire. Avec : € FTBFr = Hr(p) = A(p).B(p) 1+ A(p).B(p) = Tr(p) 1+Tr(p) Automatique Linéaire 1 1A ISMIN www.emse.fr/~dutertre/enseignement.html - 2015 36 III.2. Interprétation géométrique du passage de la boucle ouverte à la boucle fermée. Dans le cas d’un asservissement à retour unitaire tel que € H(p) = T(p) 1+T(p) Avec € T(p) =GBO (p).e jφ BO (p ) tel que GBO soit le gain et φBO soit la phase en boucle ouverte. Soit € H(p) = GBO (p).e jφ BO (p ) 1+GBO (p).e jφ BO (p ) = GBO (p) GBO (p) + e − jφ BO (p ) € H( jω ) = GBO (ω) GBO (ω) + cosφBO (ω ) − j.sinφBO (ω) En notant € H( jω ) =GBF (ω).e jφ BF (ω ) On obtient € GBF (ω) = GBO (ω) GBO 2(ω ) +1+ 2.GBO (ω ).cosφBO (ω) € φBF (ω) = arctan sinφBO (ω ) GBO (ω) + cosφBO (ω) Ainsi, connaissant le FTBO (GBO, φBO) pour un ω donné on en déduit la FTBF (GBF, φBF). a – Abaque de Black-Nichols. L’abaque de Black-Nichols permet de repérer par un système de doubles coordonnées les valeurs de la FTBO et de la FTBF correspondante (pour un retour unitaire uniquement) dans le plan de Black. Dans le système de coordonnées rectangulaires (φBO, GBO dB), on trace les courbes isomodules GBF dB = cte (en traits continus) et isophases φBF = cte (en pointillés) de la FTBF. Courbe isomodule GBF dB = cte , tracée d’après : € T( jω ) = H( jω ) 1−H( jω ) = GBF .e jφ BF 1−GBF .e jφ BF € T( jω ) = GBF −GBF + cosφBF − j sinφBF Automatique Linéaire 1 1A ISMIN www.emse.fr/~dutertre/enseignement.html - 2015 37 Soit € GBO = GBF GBF 2 − 2.GBF .cosφBF +1 € φBO = arctan sinφBF cosφBF −GBF Le tracé des courbes isomodules (traits continus verts sur la figure III.5) se fait pour GBF dB = cte , dans le repère (φBO, GBO dB) en faisant varier φBF. Le tracé des courbes isophases (traits pointillés rouges, fig. III.5) se fait pour φBF = cte , dans le repère (φBO, GBO dB) en faisant varier GBF dB. Voir également l’abaque donné en annexe 4. Fig. III.5 – Courbes isomodules et isophases de l’abaque de Black-Nichols. Automatique Linéaire 1 1A ISMIN www.emse.fr/~dutertre/enseignement.html - 2015 40 La méthode à suivre est illustrée figure III.9 : partant de l’isomodule correspondant au gain statique (HBF dB(0) en rouge), on cherche le point d’intersection de la FTBO avec l’isomodule HBF dB(0) - α dB (en bleu) qui donne la pulsation de coupure correspondante ωc. Fig. III.9 – Détermination de la bande passante à -α dB en BF. III.3. Réponse impulsionnelle d’un système bouclé en régime linéaire. L’étude de la réponse impulsionnelle (réponse à un Dirac) d’un système bouclé permet d’aborder la question de la stabilité. D’après la définition 6, un système stable écarté de son point de repos doit y retourner. En exprimant la FTBF, H(p), sous forme d’une fraction rationnelle (cf. équation II.15, partie II.3) : € H(p) = N(p) D(p) = bm p m + bm−1p m−1 + ...+ b1p + b0 an p n + an−1p n−1 + ...+ a1p + a0 (tq n ≥ m) Et d’après : € xc (t) = δ(t) TL⎯ → ⎯ XC (p) =1 On a : € Y (p) = H(p).1 = bm p m + bm−1p m−1 + ...+ b1p + b0 an p n + an−1p n−1 + ...+ a1p + a0 Automatique Linéaire 1 1A ISMIN www.emse.fr/~dutertre/enseignement.html - 2015 41 Ce qui nous conduit à (cf. partie II.3) : € Y (p) = A1 p − p1 + ...+ An p − pn Ai complexe Ainsi, lors du retour dans le domaine temporel, on a pour chaque pôle pi de H(p) : € Ai p − pi TL−1⎯ → ⎯ Aie pi t Pour que ces exponentielles ne divergent pas vers +∞ il faut que la partie réelle de chaque pi soit strictement négative. Conclusion : un système de transmittance H(p) est stable si et seulement si tous ses pôles sont à partie réelle strictement négative (c’est-à-dire que les zéros de D(p) = 1 + T(p) sont à partie réelle négative). On remarque donc qu’il suffit de connaitre les zéros de 1 + T(p) pour conclure sur la stabilité d’un système. D’où la règle fondamentale de l’automatique : La connaissance de la FTBO permet de conclure sur la stabilité du système en boucle fermée. En dehors du calcul direct des pôles (par ordinateur), il existe deux grandes familles de critère pour étudier la stabilité d’un système : - les critères algébriques, - les critères géométriques. III.4. Le critère de Routh-Hurwitz (critère algébrique). Considérons un système de FTBF6 : € H(p) = N(p) D(p) = N(p) an p n + an−1p n−1 + ...+ a1p + a0 L’étude du polynôme caractéristique D(p) = 0 (ou polynôme d’Hurwitz) permet de conclure sur la stabilité du système (le critère est énoncé ci-après sans être démontré). 6 L’attention du lecteur est attirée sur le fait qu’il s’agit du seul critère basé sur l’étude de la FTBF (les critères géométriques présentés ci-après reposent sur l’étude de la FTBO). Automatique Linéaire 1 1A ISMIN www.emse.fr/~dutertre/enseignement.html - 2015 42 Le critère de stabilité de Routh se décompose en deux conditions : • Une condition nécessaire : la stabilité exige que tous les coefficients ai soient de même signe et non nuls. • Une condition nécessaire et suffisante : le système est stable (i.e. les zéros de D(p), c'est-à-dire les pôles de H(p), sont tous à partie réelle strictement négative) si et seulement si tous les termes de la 1ère colonne du tableau de Routh sont de même signe. Construction du tableau de Routh : € pn an an−2 an−4 .... pn−1 an−1 an−3 an−5 .... pn−2 an−1an−2 − anan−3 an−1 = b1 an−1an−4 − anan−5 an−1 = b2 b3 = ... pn−3 p1 p0 b1an−3 − an−1b2 b1 = c1 c2 = ... Les deux premières lignes du tableau de Routh sont obtenues en reportant les coefficients du polynôme caractéristique (les emplacements vides correspondent à la valeur zéro). Les coefficients des lignes suivantes sont calculés à partir des coefficients des deux lignes immédiatement supérieure et correspondant plus précisément à la 1ère colonne et à la colonne suivante (le cadre et le gamma inversé gris clairs ajoutés au tableau illustrent le calcul de b1). Exemple 1 : conclure quant à la stabilité du système de FTBF : € H(p) = 1 p5 + p4 − 2p3 + 2p2 + p + 5 Solution : les coefficients associés au polynôme d’Hurwitz ne sont pas tous de même signe. La condition nécessaire de stabilité n’est pas vérifiée. Le système associé est donc instable. Exemple 2 : conclure quant à la stabilité du système de FTBF : € H(p) = K 1+τ.p( ) p4 + 5p3 + 3p2 + 6p +1 …. Automatique Linéaire 1 1A ISMIN www.emse.fr/~dutertre/enseignement.html - 2015 45 Fig. III.11 – Contour de Bromwich. D’après le théorème de Cauchy, si p décrit (C), le contour de Bromwich, alors 1 + T(p) décrit une courbe fermée (Γ) faisant N tours dans le sens horaire autour de l’origine avec : N = Z+ - P+ tel que P+ = pôles à partie réelle > 0 de 1 + T(p) P+ = pôles à partie réelle > 0 de T(p) P+ = PBO+ pôles instables de la boucle ouverte et Z+ = zéros à partie réelle > 0 de 1 + T(p) Z+ = pôles à partie réelle > 0 de Z+ = PBF+ pôles instables de la boucle fermée Soit N = PBF+ - PBO+ Ainsi, connaissant PBO+ et N, on en déduit le nombre de pôles instables de la FTBF (i.e. la stabilité). Point critique (-1,0) : plutôt que de regarder le nombre de tours de (Γ) = {1 + T(p)} autour de l’origine, on trace (Γ ‘ ) = {T(p)} (la FTBO) et on regarde le nombre de tours autour du point critique (-1,0). Ce résultat est obtenu par une simple translation de -1 selon l’axe réel. On est alors à même d’énoncer le : Automatique Linéaire 1 1A ISMIN www.emse.fr/~dutertre/enseignement.html - 2015 46 Critère de stabilité de Nyquist : lorsque p décrit le contour de Bromwich (C) dans le sens horaire, T(p) (la FTBO) décrit une courbe (Γ) dans le plan complexe. Le système est stable en boucle fermée si et seulement si le nombre de tours de (Γ) autour du point critique -1 comptés dans le sens horaire est égal à moins le nombre de pôles instables de T(p) (pôles Re>0 de la FTBO). Pour un système stable en boucle ouverte, le système est stable en boucle fermée si (Γ) n’entoure pas le point critique. Rmq : pour un système instable en BF, le nombre de pôles instables est PBF+ = N + PBO+. La connaissance de la boucle ouverte permet de conclure sur la stabilité en boucle fermée. Comportement à l’origine : le contour d’exclusion à l’origine (pour p = 0) permet d’éviter les problématiques d’existence de la FTBO dans le cas où elle peut s’exprimer sous la forme : € T(p) = a pα .1+ b1p + ...+ bm p m 1+ a1p + ...+ an p n avec α > 0 Pour p ≈ 0 (au voisinage de 0) : Ainsi, en parcourant le demi-cercle de rayon r, l’argument de p passe de -π/2 à +π/2 en faisant un demi-tour autour de l’origine dans le sens antihoraire, et donc l’argument de T(p) passe de απ/2 à -απ/2 en tournant dans le sens horaire. Ainsi, on passe de T(0-) à T(0+) en faisant une rotation de απ dans le sens horaire. L’exemple suivant, qui sera corrigé en TD, illustre ce comportement. Exemple : Etudier la stabilité du système de FTBO, € T(p) = K p.(p +1).(p + 2) , par application du critère de Nyquist. Automatique Linéaire 1 1A ISMIN www.emse.fr/~dutertre/enseignement.html - 2015 47 b – Le critère du revers. Le critère du revers est une simplification du critère de Nyquist pour les systèmes simples. Conditions suffisantes devant être vérifiées par la FTBO pour pouvoir appliquer le critère du revers : - système stable en boucle ouverte, - ordre de la FTBO > 1, - T(p) à phase minimale (i.e. pas de zéro à Re>0), - le rapport des coefficients de plus bas degré du numérateur et du dénominateur est positif. Critère du revers : si un système vérifie les conditions suffisantes exposées précédemment et si le lieu de Nyquist de la FTBO, décrit dans le sens des pulsations croissantes (0+ à +∞), laisse le point critique à sa gauche, alors, le système sera stable en boucle fermée. La figure III .12 illustre l’application du critère du revers pour des systèmes instable et stable. Fig. III.12 – Illustration du critère du revers. Critère du revers dans le plan de Black : si un système en boucle ouverte vérifie les conditions suffisantes énoncées précédemment, et, si son lieu dans le plan de Black parcouru dans le sens des pulsations croissantes (0+ à +∞) laisse le point critique à sa droite alors le système est stable en boucle fermée. Automatique Linéaire 1 1A ISMIN www.emse.fr/~dutertre/enseignement.html - 2015 50 IV. Performances des systèmes asservis. Dans ce chapitre on considérera que les systèmes étudiés sont stables (et à retour unitaire). Les deux critères de performance étudiés sont la précision et la rapidité (cf. I.3.c). IV.1. Précision. Définition 11 : Estimer la précision d’un système asservi c’est mesurer ou prédire l’évolution temporelle de l’écart entre la consigne d’entrée et la sortie du système : ε(t) = yc(t) – y(t). Le but étant de minimiser ε(t). Le système est susceptible d’évoluer sous l’effet d’une modification de la consigne yc(t) ou de l’apparition de perturbations extérieures n(t). Fig. IV.1 – Schéma bloc d’un système à retour unitaire. D’après la figure IV.1 on peut écrire : € Y (p) = T(p). Yc (p) −Y (p)[ ] + B(p).N(p) € Y (p). 1+T(p)[ ] = T(p).Yc (p) + B(p).N(p) € Y (p) = T(p) 1+T(p) FTBF ! " # $ # .Yc (p) + B(p) 1+T(p) .N(p) Ainsi, d’une façon générale, on peut décomposer l’étude en deux : - d'une part l’étude de la poursuite : évolution de l’erreur pour les variations de la consigne en l’absence de perturbations, - et d’autre part l’étude en régulation : évolution de l’erreur en présence de perturbations pour une consigne fixe. Automatique Linéaire 1 1A ISMIN www.emse.fr/~dutertre/enseignement.html - 2015 51 a – Précision statique en poursuite – Erreur en régime permanent. L’erreur en régime permanent est : € limt→+∞ε(t) = εs Donc, d’après le théorème de la valeur finale : € εs = limt→+∞ε(t) = limp→ 0 pε(p) Avec € ε(p) =Yc (p) −Y (p) € ε(p) =Yc (p) − T(p) 1+T(p) .Yc (p) € ε(p) = Yc (p) 1+T(p) D’où € εs = limp→ 0 p.Yc (p) 1+T(p) Eq. IV.1 L’erreur statique, εs, dépend du signal de consigne et de la FTBO. Dans le cas où l’on peut écrire la FTBO sous la forme : € T(p) = a pα .1+ b1p + ...+ bm p m 1+ a1p + ...+ an p n (système de classe α, i.e. nb d’intégrations pures) On a alors : € limp→ 0T(p) = limp→ 0 K pα La FTBO correspond alors aux exemples de diagrammes de Bode asymptotique (gain uniquement) donnés figure IV.2. Fig. IV.2 – Allures du gain de systèmes de classe 0, 1 et 2. Automatique Linéaire 1 1A ISMIN www.emse.fr/~dutertre/enseignement.html - 2015 52 Réponse à un échelon de consigne (erreur de position) : € yc (t) = Γ(t) TL⎯ → ⎯ Yc (p) =1/ p D’où € εs = limp→ 0 1 1+T(p) ≅ limp→ 0 1 1+ K pα et donc pour α = 0 : et pour α ≥ 1 : Un système qui possède au moins un intégrateur (α ≥ 1) en boucle ouverte a une erreur de position nulle. Réponse à une rampe (erreur de trainage) : € yc (t) = t TL⎯ → ⎯ Yc (p) =1/ p 2 D’où pour α = 0 € εs →+∞ pour α = 1 € εs =1 K pour α ≥ 2 Un système qui possède au moins deux intégrateurs (α ≥ 2) en boucle ouverte a une erreur de traînage nulle. Réponse à une parabolique (erreur en accélération) : € yc (t) = t 2 2 TL⎯ → ⎯ Yc (p) =1/ p 3 D’où € εs = limp→ 0 1 p2 . 1 1+ K pα pour α = 0 ou 1 € εs →+∞ pour α = 2 € εs =1 K pour α ≥ 3 Un système qui possède au moins trois intégrateurs (α ≥ 3) en boucle ouverte a une erreur d’accélération nulle. Automatique Linéaire 1 1A ISMIN www.emse.fr/~dutertre/enseignement.html - 2015 55 Le tableau IV.2 synthétise les valeurs de εs pour les différentes classes des blocs A et B. Tab. IV.2 – Précision statique en régulation (échelon unitaire). On déduit de la dernière colonne du tableau précédent que l’erreur statique engendrée par une perturbation assimilable à un échelon est nulle s’il existe au moins une intégration en amont de la perturbation. d – Critères de performance. Plusieurs critères de performance permettent d’estimer l’erreur (et la rapidité de sa disparition), pour les présenter, on considère un système stable dont l’erreur statique en réponse à un échelon est nulle, cf. figure IV.4. Fig. IV.4 – Visualisation de l’erreur et de l’erreur quadratique. Un premier critère consiste à intégrer le carré de l’erreur : € I = ε 2 0 +∞ ∫ (t)dt La meilleure précision est obtenue pour I faible. Automatique Linéaire 1 1A ISMIN www.emse.fr/~dutertre/enseignement.html - 2015 56 On rencontre d’autres façons de calculer le critère : € I = ε(t) .dt 0 +∞ ∫ afin de minimiser l’influence excessive du 1er lobe de ε2 Ou encore : € I = t.ε(t) .dt 0 +∞ ∫ (ITAE, Integral of Time Multiplied Absolute value of Error) pour introduire une pondération temporelle du poids des lobes. IV.2. Rapidité des systèmes. On cherche à obtenir des systèmes asservis une réponse rapide aux variations de la consigne et une aptitude à effacer rapidement les perturbations. Le temps de réponse à 5% donne une bonne évaluation de la rapidité d’un système, il exprime le temps mis par le processus soumis à un échelon pour atteindre sa valeur de régime permanent à ±5% près (et y rester). La figure IV.5, illustre la mesure du temps de réponse à 5% pour un 1er et un 2nd ordre. Fig. IV.5 – Illustration de la mesure du temps de réponse à 5%. Le temps de réaction est lié à l’inertie du système considéré, c'est-à-dire à ses propriétés physiques (réaction rapide d’un micro moteur, dynamiques lentes d’un réacteur nucléaire). Elles ne sont pas modifiables. Dans le cas d’un 1er ordre, on a établi la relation liant le temps de montée et la fréquence de coupure, € tm . fc = 0,35 (cf. §II.1.f) ; le système étant d’autant plus rapide (tm faible) que sa bande passante est large. D’une façon générale : l’amélioration de la rapidité (propriété temporelle) d’un système passe par l’élargissement de sa bande passante (propriété fréquentielle). Pour préciser cette notion, étudions l’effet du bouclage à retour unitaire d’un système du 1er ordre T(p) (cf. figure IV.6). Automatique Linéaire 1 1A ISMIN www.emse.fr/~dutertre/enseignement.html - 2015 57 Fig. IV.6 – Asservissement d’un système du 1er ordre7. Calcul de la FTBF : € H(p) = K .T(p) 1+K .T(p) = ... = K1 1+τ1p H(p) 1er ordre avec € K1 = K1+K et τ1 = τ 1+K On en déduit les pulsations de coupure en boucle ouverte : € ω c =1/τ et en boucle fermée : € ω c1 =1/τ1 = (1+K).ω c Ainsi, le passage en boucle fermée entraîne une multiplication par 1 + K de la pulsation de coupure et, par conséquence, une division par 1 + K du temps de montée (cf. € tm . fc = 0,35). On notera, que les propriétés intrinsèques du processus physique T(p) ne sont pas modifiées par le passage en boucle fermée, T(p) a toujours une pulsation de coupure égale à 1/τ. Cependant, le système bouclé se comporte globalement comme s’il avait une bande passante de 1/τ1. C’est la commande générée en entrée de T(p) par la boucle d’asservissement qui, du fait de son amplitude, permet d’accélérer la réponse comme illustré figure IV.7 pour la boucle ouverte (à gauche) et pour la boucle fermée (à droite). La forme de u(t), la loi de commande du système en boucle fermée, appelle deux remarques : - en pratique, la pointe de commande de u(t) est limitée, si on atteint la saturation, le système a un comportement non linéaire et l’étude perd sa validité mathématique, - un à-coup de démarrage trop important sur u(t) n’est pas toujours souhaitable, par exemple, pour éviter un vieillissement trop rapide d’un système mécanique, ou pour des raisons de surconsommation d’énergie. 7 Le lecteur attentif aura relevé que l’utilisation de la transmittance de Laplace avec des signaux temporels est abusive. Automatique Linéaire 1 1A ISMIN www.emse.fr/~dutertre/enseignement.html - 2015 60 On a bien vu au IV.2 que la correction proportionnelle permet d’améliorer la rapidité d’un système d’autant plus efficacement que K est élevé. Mais cela finit par poser un problème de stabilité comme illustré dans le plan de Nyquist à la figure V.2 pour des valeurs de K croissantes. Fig. V.2 – Correction proportionnelle croissante dans le plan de Nyquist. L’augmentation de K s’accompagne d’une proximité croissante du point critique et donc d’un risque d’instabilité croissant. La figure V.3 illustre avec un diagramme de Bode l’apparition d’une instabilité liée à une constante de proportionnalité trop élevée. Fig. V.3 – Illustration d’une correction proportionnelle dans le diagramme de Bode et du risque d’instabilité associé. On notera que la correction proportionnelle a un effet sur l’ensemble de la FTBO (quelle que soit la pulsation). Automatique Linéaire 1 1A ISMIN www.emse.fr/~dutertre/enseignement.html - 2015 61 c – La commande intégrale. Une loi de commande intégrale permet d’obtenir un démarrage progressif et un effet persévérant. Loi de commande intégrale : € u(t) = 1 Ti . ε(x)dx 0 t ∫ avec Ti constante de temps d’intégration. La figure V.4 illustre l’effet de la correction intégrale pour un échelon de consigne. La commande du système, u(t), est alors une rampe de pente initiale E0/Ti, la pente va diminuer au fur et à mesure que le système va rejoindre la consigne (et donc que l’erreur va décroître), pour finalement se stabiliser à la valeur atteinte lorsque l’erreur s’annule (dans l’hypothèse où c’est le cas). On parle de commande persévérante car la loi de commande u(t) garde une valeur constante non nulle après annulation de l’erreur. Fig. V.4 – Réponse d’un correcteur intégral à un échelon de consigne. d – La commande dérivée. La dérivée de l’erreur ε apporte une information supplémentaire sur son sens d’évolution (augmentation ou diminution) et sur sa dynamique. D’où l’intérêt d’ajouter un terme dérivé au terme proportionnel d’une loi de commande pour exploiter cette information, soit : € u(t) = K ε(t) +Td . ʹ ε (t)[ ] Ainsi, si l’erreur est croissante (ε’ > 0), on obtient un terme correctif qui s’ajoute à la correction proportionnelle afin d’accélérer sa correction, et à l’inverse, si l’erreur est décroissante (ε’ < 0), le terme correctif est retranché ce qui permet de rejoindre en douceur la consigne. Cette propriété trouve une application lorsque tout dépassement de la consigne doit Automatique Linéaire 1 1A ISMIN www.emse.fr/~dutertre/enseignement.html - 2015 62 absolument être évité ; l’exemple de l’accostage à quai d’un navire présenté figure V.5 en est une parfaite illustration. Fig. V.5. – Illustration de la commande dérivée – Accostage d’un navire. L’erreur, ε, représente la distance entre la navire et le quai. Lors de l’accostage une stratégie optimale consiste à inverser les moteurs peu avant de toucher le quai de façon à l’aborder avec une vitesse nulle. Sur la partie droite de la figure sont reportées, K.ε, la composante proportionnelle, et Td.ε’, la composante dérivée de la loi de commande u. Le terme dérivé est toujours négatif (et décroissant en valeur absolue), il permet d’obtenir la loi de commande négative désirée (inversion des moteurs) à partir du point d’inversion I, et au final un accostage sans dépassement. Automatique Linéaire 1 1A ISMIN www.emse.fr/~dutertre/enseignement.html - 2015 65 Fig. V.8 – Effet dans le plan de Black d’une correction à avance de phase. Le choix des paramètres du correcteur (τ et a) est essentiel, si la contrainte € 1 τ <ω r < 1 aτ n’est pas respectée, la correction apportée peut être au mieux sans effet sur les marges de stabilité et au pire provoquer une instabilité du système comme représenté figure V.9. Fig. V.9 – Instabilité due à un mauvais réglage du correcteur à avance de phase. Automatique Linéaire 1 1A ISMIN www.emse.fr/~dutertre/enseignement.html - 2015 66 V.3. Correction proportionnelle et intégrale (P.I.) – Correction à retard de phase. La transmittance de Laplace d’un correcteur proportionnel et intégral est : € C(p) = K . 1+ 1τ.p ⎛ ⎝ ⎜ ⎞ ⎠ ⎟ Eq. V.2 Le terme intégral apporte une amélioration de la précision (cf. partie IV.1). La transmittance de Laplace d’un correcteur à retard de phase est : € C(p) = K . 1+τ.p 1+ bτ.p avec b >1 Eq. V.3 Ce correcteur apporte un gain et un retard pour les basses fréquences. Les diagrammes de Bode de ces deux correcteurs sont donnés figure V.10 (tel que K=1 pour le correcteur P.I.). Fig. V.10 – Corrections proportionnelle intégrale et à retard de phase (diagrammes de Bode). Le correcteur PI permet une augmentation importante du gain en boucle ouverte aux basses fréquences ce qui permet d’améliorer la précision ; sans modification de la FTBO à proximité du point critique et donc sans dégradation des marges de stabilité (si τ est bien choisi). La figure V.11 illustre dans le plan de Black cet effet pour . On constate bien que la FTBO n’est modifiée que pour les basses fréquences (inférieures à 1/τ) et que de ce fait les marges de stabilité ne sont pas modifiées (la FTBO ne change pas à proximité du point critique). Automatique Linéaire 1 1A ISMIN www.emse.fr/~dutertre/enseignement.html - 2015 67 Fig. V.11 – Effet d’un correcteur P.I. dans le plan de Black (réglage correct). Pour un mauvais choix de τ ( ), les marges de stabilité sont dégradées (cf. figure V.12). Fig. V.12 – Dégradation des marges de stabilité pour un mauvais réglage du P.I. La figure V.13 illustre l’effet d’un correcteur à retard de phase dans le plan de Black pour un réglage correcte des paramètres ( ). Automatique Linéaire 1 1A ISMIN www.emse.fr/~dutertre/enseignement.html - 2015 70 Fig. V.15 – Effet d’un correcteur P.I.D. dans le plan de Black. Exemple : la figure V.16 présente la correction d’un système par un P.I.D. dans le plan de Black. Avant correction (en bleu), on relève Mφ = 25° et MG = 4,75 dB, ce qui est trop faible et € FTBO dB ω =0 =12 dB soit un gain statique Ks = 4 synonyme d’une précision limitée. Le but, en choisissant une correction P.I.D., est d’obtenir une erreur de position nulle (terme intégral), d’améliorer les marges de stabilité et d’avoir un dépassement indiciel inférieur à 10% (ce qui correspond à un facteur de résonance MdB = 0,5 dB en assimilant le système à un 2nd ordre, cf. annexe 3). Le pivot du P.I.D. est choisi en ω = 0,15 rad/s en prenant τd = 4 s et τi = 12 s. Pour K = 1, courbe verte, on obtient MdB = 1,1 dB (isomodule tangente à la FTBO) soit un dépassement maximal de 15%. En diminuant K on abaisse la courbe, en choisissant K = 0,4 (courbe rouge) on atteint un facteur de résonance quasi nul (dépassement proche de zéro). Automatique Linéaire 1 1A ISMIN www.emse.fr/~dutertre/enseignement.html - 2015 71 Fig. V.16 – Exemple P.I.D. (Plan de Black). La figure V.17 présente les résultats précédents dans le domaine temporel. Fig. V.17 – Correction P.I.D. dans le domaine temporel. Automatique Linéaire 1 1A ISMIN www.emse.fr/~dutertre/enseignement.html - 2015 72 V.5. Modèle du second ordre. Un critère général de bonne qualité d’un asservissement est la possibilité d’assimiler sa FTBF à un système du second ordre de coefficient d’amortissement 0,5 < m < 0,7, de pulsation propre non amortie élevée et de gain statique proche de un (cela correspond à un temps de réponse à 5% faible, à l’absence de déplacement ou à un dépassement faible, à une pseudo- période Tp élevée ayant pour conséquence un nombre d’oscillations faible). Une méthode usuelle de synthèse d’un correcteur consiste à le choisir tel qu’il modifie la FTBO du système dans le plan de Black de façon à approcher une FTBO corrigée donnant en boucle fermée un comportement similaire à un 2nd ordre répondant aux critères précédents. La lecture de la FTBO dans le plan de Black permet grâce aux tables données en annexe 3 d’en déduire un modèle équivalent du 2nd ordre en boucle fermée et de dessiner la réponse temporelle à un échelon. En effet, l’abaque de Black-Nichols permet de retrouver le facteur de résonance en boucle fermée puis le coefficient d’amortissement à partir des tables ; l’abaque donne également la pulsation de résonance en boucle fermée, on en déduit la pulsation propre non amortie (connaissant m). Exemple. Etablir le modèle du 2nd équivalent de la FTBO représentée figure V.18 (en bleu), et les principaux paramètres permettant de tracer sa réponse temporelle à un échelon. En ω = 0, la FTBO coupe l’isomodule à -1,4 dB on en déduit € FTBF dB ω =0 = −1,4 dB . La FTBO est tangente à l’isomodule 3,6 dB, d’où € ω r BF = 0,29 rad /s et € FTBF dB ω =ω r BF = 3,6 dB . On en déduit le facteur de résonance : € MdB = FTBF dB ω =ω r BF − FTBF dB ω =0 = 3,6 − −1,4( ) = 5dB Puis grâce au tableau de l’annexe 3 : m = 0,3 et € ω r BF /ω 0 = 0,91 d’où ω0 = 0,32 rad/s. Calcul du gain statique : € 20logK = −1,4 dB ⇒ K = 0,85. Automatique Linéaire 1 1A ISMIN www.emse.fr/~dutertre/enseignement.html - 2015 75 Annexe 1 – Signaux type. Rampe unitaire causale. Echelon unitaire Γ(t) - Fonction de Heaviside. Impulsion unitaire δ(t) – Impulsion de Dirac. Automatique Linéaire 1 1A ISMIN www.emse.fr/~dutertre/enseignement.html - 2015 76 Signal sinusoïdal (périodique). Signal causal retardé. Produit de convolution. pour des signaux causals : Pour un S.L.T.I, en notant h(t) sa réponse impulsionnelle on a: Automatique Linéaire 1 1A ISMIN www.emse.fr/~dutertre/enseignement.html - 2015 77 Annexe 2 – Transformée de Laplace. Transformée de Laplace monolatérale. € TL f (t)[ ] = F(p) = f (t).e− pt dt 0 +∞ ∫ € p =σ + jω Linéarité. Convolution. € a. f (t) + b.g(t) TL⎯ → ⎯ a.F(p) + b.G(p) Fonction de transfert – H(p). € y(t) = (u∗h)(t) TL⎯ → ⎯ Y (p) =U(p).H(p) H(p) = Y (p) U(p) Dérivation en temps. Intégration. € TL f (θ)dθ 0 t ∫[ ] = 1p .F(p) Dérivation en p. Translation en p. € dF(p) dp = TL −t. f (t)[ ] € TL eαt . f (t)[ ] = F(p −α) Théorème du retard temporel. € g(t) = f (t −τ) → TL g(t)[ ] = e−pτ .F(p) Théorème de la valeur initiale. Théorème de la valeur finale. € limp→+∞p.TL f (t)[ ] = f (0+) € limp→ 0p.TL f (t)[ ] = limt→+∞ f (t) = f (+∞) Transformées de Laplace usuelles. Dirac : € δ(t) TL⎯ → ⎯ 1 € sin(ωt) TL⎯ → ⎯ ω p2 +ω 2 Échelon : € Γ(t) TL⎯ → ⎯ 1 p € cos(ωt) TL⎯ → ⎯ p p2 +ω 2 Rampe : € e−a.t .cos(ωt) TL⎯ → ⎯ p + a p + a( )2 +ω 2 € e−a.t TL⎯ → ⎯ 1 p + a € e−at .sin ωt( ) TL⎯ → ⎯ ω p + a( )2 +ω 2 € t.e−a.t TL⎯ → ⎯ 1 p + a( )2 Automatique Linéaire 1 1A ISMIN www.emse.fr/~dutertre/enseignement.html - 2015 80 Annexe 3 – Systèmes linéaires du second ordre. Pulsation de résonance € ω r =ω 0 1− 2.m 2 Pulsation de coupure € ω c =ω 0 1− 2.m 2 + 1+ 1− 2.m2( )2 Facteur de résonance € M = 1 2m 1−m2 Facteur de qualité Temps de montée € tm = 1 ω 0 1−m 2 . π − arccosm( ) Temps de réponse à n% (m<0,7) € tr ≅ 1 ω 0 .m .ln 100 n ⎛ ⎝ ⎜ ⎞ ⎠ ⎟ Temps de pic € tpic = π ω 0 1−m 2 Pseudo-période Dépassement € D% =100.e −π .m 1−m 2 Nombre d’oscillations complètes Automatique Linéaire 1 1A ISMIN www.emse.fr/~dutertre/enseignement.html - 2015 81 m D % MdB m 0,1 0,15 0,2 0,25 0,3 0,35 0,4 0,45 0,5 0,55 0,6 0,65 1,68 1,74 1,81 1,88 1,97 2,06 2,16 2,28 2,42 2,58 2,77 3 30 20 14 11 10,1 7,9 7,7 5,4 5,3 5,3 5,2 5 3,16 3,18 3,21 3,24 3,29 3,35 3,43 3,52 3,63 3,76 3,93 4,13 6,31 6,36 6,41 6,49 6,59 6,71 6,86 7,04 7,26 7,52 7,85 8,27 73 62 53 44 37 31 25 21 16 12,6 9,5 6,8 0,99 0,98 0,96 0,94 0,91 0,87 0,82 0,77 0,71 0,63 0,53 0,39 1,54 1,53 1,51 1,48 1,45 1,42 1,37 1,33 1,27 1,21 1,15 1,08 1,56 1,56 1,57 1,59 1,61 1,63 1,67 1,72 1,8 1,93 2,17 2,74 14 10,5 8,1 6,3 4,8 3,6 2,7 1,9 1,2 0,7 0,3 0,1 0,1 0,15 0,2 0,25 0,3 0,35 0,4 0,45 0,5 0,55 0,6 0,65 0,7 3,29 3 4,4 8,8 4,6 0,14 1,01 7,14 0 0,7 0,75 0,8 0,85 0,9 0,95 3,66 4,16 4,91 6,17 9,09 3,1 3,4 3,7 4 4,1 4,75 5,24 5,96 7,21 10,06 9,5 10,5 11,93 14,41 20,12 2,84 1,52 0,63 0,15 0,01 - - - - - 0,94 0,87 0,81 0,75 0,69 - - - - - - - - - - 0,75 0,8 0,85 0,9 0,95 Automatique Linéaire 1 1 Ts réponse fréquentielle 2" ORDRE wee) srr ws " -1 0 1 2s 0, 0dB fou 0,5 =o point d'inflexion 4 6 8 10 réponse indicielle systéme du 2" ordre 1 144 24(2 a \oy me www.emse.fr/~dutertre/enseignement.html - 2015