Introduction à la commande par ordinateur

1. Introduction à la commande par ordinateur Question étudiée dans ce cours : comment un ordinateur embarqué peut-il prendre les commandes d’un processus physique ? Les notions abordées : • processus physique : entrées, sorties, équations, temps continu • Système asservi : capteurs, actionneurs, loi de contre réaction • Discrétisation : échantillonnage, blocage, temps discret • Loi de commande : boucle ouverte, retour d’état, commande optimale • Performance : rapidité, précision, stabilité, consommation La page http://cours.polytech.unice.fr/intcom contient tous les documents utilisés en cours, et liens vers les outils utilisés

2. Contenu du cours, projet à réaliser, et synthèse demandée • Le coursdes six premières semaines aborde les notions citées dans un con-texte déterministe où entrées, sorties et équations sont connues parfaitement • Les six semaines suivantes, on ajoute (J. Le Roux) la notion de bruit d’ob-servation, les informations des capteurs sont entachées d’incertitudes. • L’animation Lunar Lander (Flash, ActionScript) illustre le cours en simu-lant l’alunissage en pilotage manuel d’un module lunaire sous forme de jeu. • Un projet est réalisé en groupes de deux élèves durant les séances de travaux dirigés, il s’agit de modéliser un processus physique commandé par ordinateur, par exemple l’alunissage automatisé d’un Lunar Lander. Technologie utilisée : Html5 et javaScript, ou Processing, ou python, ou ? • Une synthèse personnelle est demandée à chacun, sur les notions du cours • Pour la séance de TD de cette semaine • Télécharger Scilab pour les calculs nécessaires, discrétisation, loi de commande • Télécharger Lunar Lander, exécutable et source, pour une étude détaillée

3. On dispose d’une illustration des notions du cours avec l’animation Lunar Lander qui simule l’alunissage d’un module lunaire en 2D Animation pédagogique tirée du site http://phet.colorado.edu/ et enrichie de plusieurs lois de commandes, Lunar Lander implémente : • des équations tirées de la relation fondamentale de la dynamique • Un modèle discrétisé à la période d’échantillonnage Te= 40 ms • un tableau de bord avec l’état des capteurs et actionneurs, ... • Plusieurs lois de commandes, ‘m’ manuel, ‘e’ pour retour d’état, • des paramètres de commande transmis par Scilab dans des fichiers texte • … on demande une analyse détaillée des fonctionnalités dans le premier td

4. Planning prévisionnel du cours des six premières semaines • Définir les notions utilisées, sur l’exemple simple du processus de remplissage d’une cuve • TD : expérimenter Scilab, Analyser Lunar Lander, apprendre html5 ou … • Écrire les équations de Lunar Lander en utilisant la seconde loi de Newton, les mettre sous la forme de la représentation d’état, et discrétiser ces équations pour le projet à l’aide de Scilab • Commande en boucle ouverte de Lunar Lander • Commande par retour d’état, commande de vitesse, commande de position • Commande optimale linéaire quadratique à horizon fini, résolution de l’équation récurrente de Ricatti • Rattrapages, finition, et bilan Débutons avec l’exemple du processus de remplissage d’une cuve

5. Qu’est ce que l’Automatique ? « l’essence des bonnes machines est de se gouverner par elles-mêmes, autant qu’il est possible et sans le secours de l’intelligence humaine »J.V. Poncelet (1826), extrait du premier ( ?) cours d’Automatique (cité dans ‘Eléments d’Automatique’, de Faurre et Robin, chez Dunod).

6. Exemple de la clepsydre, une contre-réaction antique L’horloge à eau du mécanicien grec Ktésibios, revue et améliorée par les physiciens arabes au 9ème siècle, était connue des égyptiens, et également des amérindiens. Il y a diverses versions et principes. Dans la version suivante, elle met en œuvre un régulateur de niveau dans le réservoir intermédiaire afin de réguler le débit dans le réservoir horloge : Tiré de ‘Eléments d’Automatique’, de Faurre et Robin, chez Dunod

7. Autre exemple: le régulateur à boules de James Watt Le régulateur à boules de Watt schématisé ci-contre est l’un des multiples mécanismes ingénieux développés au 18ème siècle durant la révolution industrielle. Il s’agit de stabiliser la vitesse de rotation du moteur à vapeur, en régulant la pression de la vapeur dans la chaudière. Sinon, ce processus est instable. C’est donc une application du principe de contre réaction, et Watt est présenté par les anglo-saxons comme le père des automatismes (d’après le livre déjà cité de Faurre et Robin). On règle la « tringlerie » pour maintenir la vitesse autour de la consigne wc tringlerie Les boules s’écartent quand la vitesse croît pointeau

8. Processus, entrées, sorties, relation d’entrée sortie Schéma de remplissage (manuel): Processus physique de remplissage d’une cuve avec fuite : d(t) entrée, f(t) perturbation h(t) sortie Règle: capteur Robinet: actionneur Relation d’entrée sortie, la variation du volume V=S*h dans la cuve est la différence d-f robinet d(t), débit h(t), hauteur S f(t), fuite

9. d(t) 1 0 t h(t) 1m 1s t 0 Équation du processus de remplissage, test de réponse indicielle, instabilité Entrée Bornée Sortie Bornée (EBSB) • t, le temps, varie continûment • Prenons S=1 m2, et f=0 litre par seconde cond. init. Le processus étudié est instable au sens E.B.S.B. car on a trouvé une entrée bornée pour laquelle la sortie est non bornée … la cuve déborde

10. Soit hc, la hauteur de liquide souhaitée, ou consigne, ou encore référence : Loi de commande par tout ou rien (non linéaire) : Loi de commande linéaire : a(t) actionneur r(t) capteur d(t) h cuve Capteur, actionneur, loi de commande, contre-réaction Actionneur : électrovanne Commande a(t) : entrée électrique qui commande l’ouverture de l’électrovanne et donc le débit d(t) Capteur : convertit hauteur h(t) en tension r(t) On peut maintenant automatiser le processus On fait pour simplifier : a(t)=d(t) et r(t)=h(t) d(t) est la commande k est le gain de contre réaction hc-h est l’erreur d’asservissement h(t) est le retour

11. Équation du système bouclé avec la loi de commande linéaire • l’équation du système asservi avec la loi de commande linéaire : • Équation caractéristique du système bouclé : • Une racine de l’équation caractéristique, ou pôle, ou encore valeur propre du système bouclé • La règle : un pôle réel négatif est stable au sens EBSB

12. Test indiciel du système bouclé, précision, temps de réponse • Si h(0)=0 et hc(t)=hc constante Solution homogène (sans second membre) Solution particulière avec second membre hc(t) hc D’où la solution : t 0 La cuve ne déborde plus sauf si la consigne donnée dépasse la hauteur de la cuve : h(t) hc • hc=est la valeur finale • le gain statique du système bouclé vaut 1 • La précision est parfaite, pas d’erreur h(t)hc • Temps de réponse à 5% : tr5%=3/k, car 0.95hc 0 t tr5%=3/k

13. Test du retour à l’équilibre, stabilité asymptotique • Placé à l’écart de ses conditions d’équilibre, à consigne nulle, comment un processus ou un système revient-il (ou pas) à l’équilibre L’équilibre est défini comme suit : la sortie et toutes les dérivées de la sortie sont nulles. Ici, c’est : • Ici, le système bouclé revient à l’équilibre : • le retour à l’équilibre dure tr5%=3/k • il est apériodique (pas de dépassements) h(t) h0 Stabilité asymptotique : c’est quand le processus revient à l’équilibre C’est plus contraignant que la stabilité au sens EBSB (cf. processus seul ici) 0.05h0 0 t tr5%=3/k

14. Horloge d’échantillonnage, échantillonnage de l’entrée, blocage de la sortie Commande par ordinateur implique : • Rythmée par une période d’échantillonnage Te, • blocage de l’entrée d, BoZ • échantillonnage de la sortie h(t) On passe en temps discret, les signaux changent ou sont mesurés aux instants nTe d(t) entre t= nTe et t= (n+1)Te, d(t)=d(nTe) constante, Blocage d’ordre zéro d0 d1 d2 d3 … 0 Te t h(t) … hn=h(nTe) échantillonnage h2 h1 t 0 2Te Te

15. Discrétisation du processus de remplissage de cuve • Discrétiser, c’est calculer la relation récurrente entre la commande dn et la condition initiale hn à l’instant nTe et la sortie hn+1à l’instant (n+1)Te • Pour cela, on tient compte du bloqueur d’ordre zéro, en intégrant l’équation différentielle différencielle du processus à entrée constante • Par exemple entre t=0 et t=Te, • On généralise entre nTe et (n+1)Te :

16. Commande par ordinateur, loi de commande programmée : • Équation caractéristique et racines (pôles, ou encore valeurs propres) : • Pour un système en temps discret, la stabilité EBSB s’obtient si tous les pôles sont de module inférieur à un • Test du retour à l’équilibre : • Loi de commande : • Équation du système bouclé

17. Calculer la réponse indicielle du système bouclé • La consigne est égale à hc constante

18. Comportement indiciel selon le gain de contre réaction • Racine a= 0.25 • Racine caractéristique a= 0.5

19. Comportement indiciel selon le gain de contre réaction • Racine a= 0; (réponse pile ?) • Racine caractéristique a= - 0.5

20. Valeur finale et temps de réponse du système bouclé : • Sur la réponse indicielle, on lit que h(t) tend vers hc quand t tend vers l’infini • hc=est la valeur finale • Le gain statique du système bouclé vaut 1 • La précision est parfaite, pas d’erreur d’asservissement • Le régime transitoire dure à peu près tr5%=3/k, car h(3/k) égale environ 0.95hc La cuve ne déborde plus sauf si la consigne dépasse la hauteur de la cuve