J-0
00m
00j
00h
00min
00s

Version interactive avec LaTeX compilé

Télécharger PDF

X ENS Modélisation PSI 2014

Notez ce sujet en cliquant sur l'étoile
0.0(0 votes)
Logo x-ens
🔗 Explorer
Banque
X-ENS
Année
2014
Filière
PSI
Matière
Physique Modelisation
X-ENS PSIX-ENS 2014PSI Physique ModelisationPhysique Modelisation 2014
2025_09_04_11d2a0f8dadc6748a781g

Banque commune École Polytechnique - InterENS

PSI

Session 2014

Épreuve de Modélisation

Durée: 5 heures
Aucun document n'est autorisé
Aucune calculatrice n'est autorisée
Si, au cours de l'épreuve, un candidat repère ce qui lui semble être une erreur d'énoncé, il le signale sur sa copie et poursuit sa composition en expliquant les raisons des initiatives qu'il est amené à prendre.

Régulation de niveau en coulée continue

Le secteur de l'acier est un des secteurs stratégiques industriels, qui en fait une industrie moderne à forte valeur ajoutée. Les nouveaux aciers de plus en plus hauts de gamme nécessitent des techniques de fabrication innovantes pour limiter au maximum les défauts à l'intérieur du matériau.
L'acier liquide en fusion peut être coulé suivant deux procédés :
  • la coulée en lingots,
  • la coulée continue.
La coulée en lingots tend à disparaître au profit de la coulée continue. Ce dernier type de coulée permet des gains de productivité importants car elle mobilise moins de personnel et la fabrication est continue. De plus, elle permet de mieux maîtriser les caractéristiques mécaniques du produit final.
Dans le procédé de coulée continue, l'acier en fusion contenu dans la poche est coulé dans une lingotière sans fond (voir Figure 1). Au contact de ses parois maintenues à basse température, l'acier commence à se solidifier en formant une peau. Des rouleaux extracteurs tirent ensuite l'acier vers le bas, où il passera dans un système de refroidissement à eau pour achever la solidification. À la sortie de la machine, la barre d'acier est découpée en tronçons appelés brames, qui continueront leur chemin dans la suite de la chaîne de fabrication.
Figure 1 - Description de la coulée continue
Les principaux éléments constitutifs de la machine de coulée (voir Figure 1 et Figure 2) sont énumérés ci-dessous :
  • le répartiteur est un grand récipient d'une capacité de 20 tonnes, situé entre la poche et la lingotière. Il permet d'alimenter plusieurs lingotières pour réaliser plusieurs coulées en continu,
  • la lingotière sans fond est un moule à base rectangulaire. Elle permet la solidification en surface de l'acier au contact des parois. En effet, un flux d'eau assure le refroidissement des bords de la lingotière et donc la formation d'une peau solide de l'acier. La lingotière est alimentée en acier par une busette à 2 orifices immergée dans le métal liquide.
    Le métal en fusion est brassé électromagnétiquement afin que la phase liquide reste homogène, la surface en phase liquide est protégée de l'oxydation à l'aide du laitier.
Figure 2 - Tête de la machine de coulée continue
Le débit du métal fourni par la busette est contrôlé par une quenouille qui vient ouvrir plus ou moins l'entrée de la busette, en fonction de la vitesse de déplacement de l'acier en sortie fixée par les rouleaux.
Le contrôle de la coulée est basé sur l'asservissement du niveau de métal dans la lingotière. Celui-ci doit rester le plus stable possible pour obtenir de bonnes propriétés de l'acier. En effet, il est très facile de constater des fissures, des trous sur des barres d'acier qui auraient subi ses fluctuations de niveau (voir Figure 3).
Pour ce faire, il est essentiel de modéliser le système étudié selon la démarche décrite ci-après et illustrée sur la Figure 4.
La démarche de modélisation suivie dans cette étude se compose symboliquement de deux domaines : le domaine physique où l'on fait des mesures sur un système réel et le domaine virtuel où l'on réalise des simulations sur un système virtuel (analytique, numérique,...) censées représenter le comportement du système réel. Le domaine physique est nécessaire car il est indispensable de s'appuyer sur l'observation des phénomènes réels pour proposer et valider une modélisation du système, de son milieu extérieur et de son comportement. Le solveur est un outil qui permet de résoudre les équations et qui donne le résultat de la simulation. La modélisation est pertinente si l'écart entre les réponses expérimentales et les résultats issus de la simulation est faible. Une fois le modèle validé, le domaine virtuel permet de prédire le bon fonctionnement avant la phase d'intégration sur le système réel. L'étude présentée ici s'intéresse principalement au domaine virtuel de simulation.
Figure 3 - Fissure dans un acier coulé avec une fluctuation de niveau
Figure 4 - Présentation de la démarche de modélisation
Le sujet est composé de 4 parties :
  • modélisation de la coulée, des différents éléments constitutifs du système et des perturbations;
  • identification des paramètres;
  • estimation des perturbations;
  • compensation des perturbations;
Une synthèse conclura l'étude. Même si certaines parties peuvent être traitées séparément, il est préférable de suivre le questionnement proposé.

I Modélisation du système de contrôle de la coulée

La structure de contrôle du niveau d'acier dans la lingotière est donnée sur la Figure 5.
On retrouve :
  • la lingotière qui permet de fabriquer la barre d'acier ; on notera la section de la lingotière, le débit entrant et le débit sortant,
  • la busette d'injection de l'acier dans la lingotière,
  • la quenouille qui permet de réguler le débit d'acier entrant dans la busette; on note sa position,
  • l'actionneur qui permet de déplacer la quenouille dont la consigne est notée ,
  • le capteur de niveau qui donne une image notée du niveau réel ,
  • la structure de contrôle qui régule le niveau réel par rapport à la consigne .
Figure 5 - Représentation de la structure de contrôle du niveau d'acier
Il est maintenant nécessaire de modéliser chacun des constituants afin de calculer les paramètres du correcteur.

I. 1 Analyse de la coulée dans la lingotière

On suppose que l'acier est un fluide incompressible de masse volumique . Dans la lingotière, l'acier s'écoule comme de l'eau. On le caractérise par une viscosité dont on donne l'ordre de grandeur à la température de la coulée : .
L'écoulement de l'acier n'étant pas parfait, les effets visqueux sont modélisés par une force
volumique : où est la vitesse de la particule fluide, la viscosité cinématique définie par et l'opérateur laplacien. De plus, le frottement de la paroi de la lingotière sur les particules fluides est modélisé par une force volumique définie par :
où est le nombre de Reynolds caractéristique de cet écoulement, la distance de la particule fluide à la paroi. On considère un point de l'espace repéré par ses coordonnées cartésiennes . On définit la dérivée particulaire d'une grandeur en et à l'instant par
Si la grandeur est une grandeur vectorielle, la dérivée particulaire sera définie pour chacune de ses composantes scalaires cartésiennes.
Q1. Rappeler dans le cas général l'équation de conservation de la matière. Que devient cette équation dans le cas de l'acier?
Q2. Interpréter chacun des termes composant la dérivée particulaire.
Q3. Justifier par des considérations qualitatives la forme de la force de frottement .
Q4. Établir, en appliquant le principe fondamental de la dynamique à une particule fluide, l'équation de Navier-Stokes :
où est la dérivée particulaire du vecteur vitesse .
La Figure 6 représente la carte des vecteurs vitesse dans la lingotière.
Figure 6 - Champ des vitesses moyennes.
Q5. Rappeler la signification physique du nombre de Reynolds. Grâce à la Figure 6, évaluer le nombre de Reynolds Re en précisant la démarche suivie. Sous certaines conditions à préciser, identifier les termes qui peuvent être négligés dans l'équation de Navier-Stokes. Commenter.
L'équation de Navier-Stokes n'a pas de solution analytique simple. Elle peut être résolue par des méthodes numériques après discrétisation du problème. Les algorithmes utilisés sont basés sur la méthode par éléments finis et requièrent des itérations successives.
La définition de la dérivée particulaire permet d'écrire l'accélération particulaire sous la forme :
où est appelé vecteur tourbillon. On note et les composantes du vecteur tourbillon (dans la base cartésienne).
La résolution numérique permet de calculer le vecteur tourbillon et en particulier la composante suivant l'axe qui s'exprime par :
et sont les composantes de la vitesse suivant les axes et (voir figure Figure 6).
L'acier fondu est recouvert en surface de laitier pour le protéger de l'oxydation.
Q6. Expliquer en quoi il est important d'éviter une composante trop importante du vecteur tourbillon à l'interface acier/laitier.
La Figure 7 représente sur la première carte le niveau moyen dans la lingotière et sur la deuxième carte l'écart-type de cette grandeur.
Figure 7 - Niveau moyen d'acier dans la lingotière et écart-type (vue de dessus).
Q7. Commenter la Figure 6 et la Figure 7. Conclure sur la nécessité de contrôler le niveau d'acier dans la lingotière. Déterminer une zone dans laquelle un capteur sensible au niveau moyen de l'acier pourrait être placé.
Cette étude nous a permis de comprendre certaines spécificités de l'écoulement dans la lingotière et la nécessité d'un système de contrôle de la coulée. Cependant les modèles utilisés pour cette analyse restent lourds à mettre en oeuvre et la précision obtenue n'est pas en adéquation avec celle nécessaire pour contrôler correctement le niveau d'acier dans la lingotière.

I. 2 Modélisation du comportement du système de contrôle de coulée

Le modèle précédent est bien trop complexe pour l'intégrer dans une démarche de conception de la partie commande du contrôleur.
On cherche donc à déterminer les fonctions de transfert de chacun des constituants afin de compléter le schéma bloc de la Figure 8.
Figure 8 - Modèle par schéma-bloc de la chaîne d'action

I.2.1 Modélisation macroscopique de l'écoulement dans la lingotière

On suppose que la lingotière est un parallélépipède dans lequel l'écoulement se fait avec un débit entrant et un débit sortant . La section du parallélépipède traversée par l'acier est notée . La hauteur d'acier dans le parallélépipède est notée .
Q8. Déterminer l'équation différentielle en temps reliant les différentes quantités : , et .
Q9. En déduire la fonction de transfert de la lingotière è où est la variable de Laplace. Vous préciserez les hypothèses utilisées.

I.2.2 Modélisation de l'acquisition de la hauteur de coulée

Le capteur utilisé est un capteur à courants de Foucault. Un schéma simplifié du dispositif est présenté à la figure Figure 9. Ce capteur, placé au-dessus de la surface libre de l'acier en fusion et du laitier, est composé d'une bobine parcourue par un courant à haute fréquence qui donne naissance à un champ magnétique variable . L'acier est alors soumis à ce champ magnétique dont l'intensité dépend de la distance capteur/surface libre de l'acier. L'acier est alors le siège de courants de Foucault.
Le laitier est un milieu diélectrique non magnétique sans charges ni courants libres. L'acier, métal conducteur de conductivité , est supposé linéaire, homogène et isotrope de permittivité et de perméabilité . La conductivité dépend du type d'alliage utilisé. Cette conductivité est de l'ordre de .
Figure 9 - Capteur à courants de Foucault
Q10. Énoncer les équations de Maxwell vérifiées par le champ électrique et par le champ magnétique en un point quelconque de l'acier et à l'instant . Établir une relation de dispersion simplifiée grâce à une hypothèse sur la fréquence de travail .
Q11. Après avoir introduit la distance caractéristique , exprimer le vecteur densité de courant en un point de l'acier à un instant quelconque . Donner la signification physique de .
Q12. Exprimer une condition sur la lingotière ou sur l'acier pour que la réponse du capteur soit indépendante de l'épaisseur de l'acier. Préciser si l'application de cette condition remet en cause l'utilisation d'un tel capteur en coulée continue.
Ce type de capteur peut être modélisé par un circuit couplé par mutuelle induction (voir Figure 10). La bobine excitatrice réelle est modélisée par une bobine d'inductance et de résistance . Cette bobine est alimentée par une source de tension sinusoïdale de fréquence . L'acier en fusion, siège des courants de Foucault, est modélisé par une inductance et une résistance . On note le coefficient de mutuelle inductance entre les deux circuits. vérifie la relation où est le coefficient de couplage.
Ce capteur est alimenté par une source de tension sinusoïdale, d'amplitude et de pulsation . On note le nombre complexe tel que . On associe à chaque grandeur électrique sa représentation complexe où est l'amplitude complexe et telle que .
Figure 10 - Modélisation électrique du capteur de niveau
Q13. Indiquer de quel paramètre géométrique dépend principalement le coefficient . On note le temps caractéristique des variations du niveau moyen de la surface libre de l'acier. Indiquer dans quel cadre, les variations de avec le temps pourront être négligées devant les variations des différentes grandeurs électriques.
Q14. Ecrire les équations électriques vérifiées par les amplitudes complexes des courants et et par la tension . En déduire une relation entre et .
Q15. En supposant que l'acier est un bon conducteur et en particulier que , modéliser le circuit côté excitation et montrer que
Q16. Expliquer les variations de avec le coefficient de couplage . Donner un encadrement de .
Ce capteur est placé dans un circuit oscillant. La fréquence propre de ce circuit est liée à la valeur de l'inductance . Le signal est ensuite mis en forme par un conditionneur qui est ici indispensable. En effet, la relation entre le signal électrique en sortie du capteur et le niveau de l'acier est complexe et non linéaire. De plus, l'amplitude des variations de l'inductance ne permet pas de linéariser le signal électrique image du niveau dans la lingotière. Un circuit électronique est donc placé en sortie du conditionneur permettant d'obtenir, après traitement, une relation simple entre le niveau mesuré proportionnel au signal électrique obtenu en sortie et le niveau dans la lingotière :
Q17. Conclure sur les avantages d'un tel capteur. On pourra en particulier comparer l'utilisation de ce type de capteur à un capteur inductif (à entrefer variable par exemple).

I.2.3 Modélisation de l'actionneur de la quenouille

La fonction de transfert de l'actionneur est assimilée à un modèle du premier ordre de gain unitaire et de constante de temps .

I.2.4 Modélisation de la busette

I.2.4.1 Influence de la hauteur de la busette Nous souhaitons maintenant exprimer le débit d'acier sortant de la busette et alimentant la lingotière. On suppose que l'acier en fusion est incompressible. On note le débit en sortie de la busette, é la différence d'altitude entre la sortie de la busette et la surface libre du répartiteur. On note de plus l'accélération de la pesanteur. Un schéma simplifié de l'écoulement est proposé à la Figure 11.
Nous supposons dans un premier temps que l'écoulement est parfait. On considère un tube de courant unidimensionnel (représenté sur la Figure 11). et sont des sections droites de ce tube de courant à l'interface air/acier dans le répartiteur et à l'entrée de la busette. On notera et les surfaces élémentaires associées. L'écoulement est supposé stationnaire. L'acier en fusion traverse et avec des vitesses respectives et . On note aussi et les pressions respectives toujours au niveau de sections et la masse volumique de l'acier.
Figure 11 - Modèle simplifié de l'écoulement de l'acier du répartiteur vers la busette (le schéma n'est pas à l'échelle).
Q18. En considérant que l'énergie totale du fluide se conserve dans un tube de courant, démontrer la forme suivante de l'équation de Bernoulli
où sont respectivement la pression, la vitesse et la surface élémentaire en un point d'ordonnée du tube de courant.
Q19. On prend maintenant en compte la viscosité de l'acier et les modifications du profil de vitesse associées. Représenter sur un schéma où figure un tube de courant, les vecteurs vitesse en différents points d'une même section de ce tube de courant. En considérant un tube de courant qui s'appuie sur les surfaces latérales de la busette, intégrer la relation de Bernoulli sur une section . Interpréter le terme , et le relier au débit . Dans la suite, on notera .
Q20. On s'intéresse maintenant à l'intégrale sur la surface du troisième terme . On définit alors une vitesse moyenne d'énergie cinétique telle que :
Montrer que . En déduire une nouvelle expression de la relation de Bernoulli qui prend en compte la viscosité de l'acier :
Q21. Entre et , on suppose que la variation de pression est négligeable devant les effets de la pesanteur. En appliquant la relation de la question précédente entre et retrouver la formule suivante
é
où représente la section effective de passage de l'acier au niveau de la busette, et le coefficient de perte de charge due à la viscosité de l'acier. Exprimer en fonction des vitesses précédemment définies. Donner la valeur du coefficient dans le cas d'un écoulement parfait.
En pratique, de l'argon est injecté dans la busette pour éviter qu'elle ne se bouche. Cela induit un retard au niveau du débit entrant dans la lingotière. Ce retard est bénéfique pour la régulation du niveau. En première approximation, il vaut :
éé
avec la hauteur entre le haut de la busette et la surface libre d'acier dans la lingotière (voir Figure 5).
En pratique, ce retard est identifié sur le système réel pour plus de précision.
Seule la section apparente est à déterminer précisément.
Figure 12 - Géométrie de la zone de contrôle du débit entrant dans la busette
I.2.4.2 Influence de la position de la quenouille par rapport à la busette La section apparente est la section à travers laquelle passe le fluide avant d'entrer dans la busette. La quenouille obstrue plus ou moins l'entrée de la busette, on note sa position . quand la quenouille est en contact avec la busette (voir Figure 12).
On supposera pour simplifier que la section apparente est un tronc de cône dont les différents paramètres sont donnés sur la Figure 13. On note :
  • le rayon de la sphère du bout de quenouille, le centre de cette sphère,
  • le rayon de la busette et le rayon de l'orifice d'entrée de la busette, est le centre du cercle dans le plan de représentation,
  • est la position de la quenouille,
  • est l'ouverture, cela représente la longueur du côté du cône,
  • quand alors ,
  • est une longueur constante,
  • le demi angle au sommet du cône.
Q22. Par la méthode de votre choix, montrer que la surface apparente (surface latérale du tronc de cône) peut s'écrire sous la forme où vous exprimerez la constante .
Figure 13 - Modélisation de la surface apparente
Q23. Écrire une relation géométrique lorsque et et en déduire la valeur de en fonction de .
Q24. Montrer maintenant que où vous exprimerez la constante . Vous pourrez supposer que et sont petits devant les autres longueurs en effectuant un développement limité à l'ordre 1 .
La relation précédente est valable pour où est la valeur limite.
Q25. Justifier l'existence de cette valeur limite et préciser ce que vaut la section apparente au delà de cette valeur.
I.2.4.3 Fonction de transfert de la busette D'après ce qui précède, le comportement de la busette peut être modélisé par l'équation suivante
On rappelle que la définition de la transformée de Laplace est :
Q26. Déterminer la fonction de transfert de la busette .

II Validation et Identification du modèle de la busette

La modélisation du système, détaillée dans la partie précédente, a permis de mettre en place un modèle pour chaque constituant et de proposer un modèle pour les perturbations. En pratique, certains paramètres sont inconnus et il faut les déterminer pour pouvoir dimensionner le correcteur correctement.
Les constantes inconnues vont être identifiées sur des réponses expérimentales.
On va s'intéresser uniquement à l'identification des paramètres de la busette connaissant l'excitation et la vitesse de déplacement de l'acier à la sortie de la lingotière et en mesurant la niveau d'acier dans la lingotière (voir Figure 14).
Figure 14 - Modèle adopté pour l'identification
L'excitation réelle ne peut pas être assimilée dans la réalité à un échelon, une rampe... c'est pourquoi les mesures vont être réalisées de manière discrète à l'aide d'une entrée de type séquence binaire pseudo-aléatoire.
Les mesures obtenues seront utilisées pour réaliser l'identification en minimisant une fonction erreur obtenue par une représentation d'état du système.

II. 1 Génération de la séquence pseudo-aléatoire

Le principe est d'utiliser un registre à décalage qui stocke un mot binaire . Å chaque pas de temps , le bit va être transmis pour générer une entrée:
  • si alors
  • si alors


  • quant à lui, sera remplacé par le résultat d'une opération de «OU EXCLUSIF» entre certains bits.
Figure 15 - Registre à décalage
Une des caractéristiques de cette technique est que la séquence se répète au bout d'un certain nombre de pas de temps, c'est le caractère «pseudo» de la méthode. La taille de la séquence dépend du nombre de bits utilisés. En général, la séquence est initialisée avec des 1 partout.
Q27. Pour un registre de 4 bits et un retour tel qu'au pas de temps , reproduire et compléter le tableau ci-dessous. Déterminer la taille de la séquence. correspond à la valeur envoyée en entrée, pour le pas de temps .
1 1 1 1 -
0 1 1 1 1
0 0 1 1 1

II. 2 Linéarisation et domaine de validité - Modèle d'état

Afin de réaliser la minimisation, il est nécessaire de mettre en place un modèle d'état linéaire du système.
Le problème est que le système n'est pas linéaire à cause de la fonction exponentielle de la Figure 14. Il faut proposer une approximation sous forme de fraction rationnelle de la fonction exponentielle.
Plusieurs approximations sont possibles :
  • développement de Taylor
  • approximation de Padé
Le développement de Taylor n'est pas satisfaisant sur des grandes plages de fréquences, c'est pourquoi l'approximation de Padé est utilisée classiquement.
L'approximation de Padé à l'ordre de la fonction est la fraction rationnelle calculée à l'aide de la formule suivante :
é
Q28. Déterminer l'approximation de Padé à l'ordre 1,2 et 3 .
Dans le cadre de la fonction retard, il faut approximer la fonction retard .
Q29. Déterminer le gain de la fonction de Padé d'ordre de la fonction retardée et le comparer à la valeur théorique attendue.
Sur la Figure 16 est tracée la phase en fonction de la pulsation pour différents ordres de l'approximation de Padé.
Q30. Sachant que le domaine d'identification de la phase est entre 0 et , choisir un ordre d'approximation en justifiant votre choix afin d'obtenir une précision inférieure à .
Figure 16 - Comparaison de la phase en fonction de la pulsation entre la phase exacte de la fonction exponentielle et différents ordres d'approximation de Padé.
Indépendamment de la question précédente, on retient une approximation de Padé d'ordre 1 pour approximer l'exponentielle.
Le modèle d'état consiste à relier des variables d'état, qui caractérisent l'état du système (ces variables peuvent être des grandeurs physiques ou des grandeurs virtuelles), aux variables
de commande (entrées réelles du système). Le modèle est mis sous une forme matricielle :
avec :
  • : la matrice d'état
  • : la matrice de commande
    : une variable virtuelle telle que
    Q31. Déterminer les matrices et , en commençant par exprimer la relation de transfert à partir de la Figure 14, puis en prenant en compte l'approximation de Padé à l'ordre 1, puis en repassant dans le domaine temporel.

II. 3 Identification des paramètres caractéristiques de la busette

Les mesures expérimentales obtenues avec la technique d'excitation pseudo-aléatoire permettent de connaître en fonction de et au cours du temps. Les mesures étant discrètes, on a pour chaque instant , les et .
L'objectif est de déterminer les paramètres inconnus et en minimisant une erreur définie par
é
é la valeur obtenue grâce au modèle d'état et le nombre d'échantillons temporels.
On suppose que la fonction est convexe.
Q32. Expliquer comment déterminer é. Proposer un algorithme permettant de déterminer le couple de valeurs ( ) permettant de minimiser l'erreur .
Q33. Expliquer précisément la démarche d'identification des paramètres en synthétisant les différentes étapes vues depuis le début de cette partie et en insistant sur les limitations.
Cette partie a permis de mettre en place une technique d'identification des paramètres du modèle basé sur l'écriture d'un modèle d'état et la minimisation d'une erreur quadratique.
Pour avoir un modèle correct, il convient encore de connaître plus précisément les perturbations. C'est l'objet de la partie suivante qui développe une méthode d'estimation des perturbations.

III Modélisation des perturbations

Il y a plusieurs perturbations qui peuvent altérer le fonctionnement du système. Les deux principales sont :
  • le bouchage
  • le gonflement
Le bouchage se situe au niveau de la busette : un dépôt de métal se colle sur la paroi de la busette ce qui l'obstrue et réduit le débit, puis le dépôt s'arrache réaugmentant le débit.
Q34. En notant le débit induit par cette perturbation, proposer une modification du schéma bloc de la Figure 8.
Le gonflement correspond à un mouvement de la peau du métal passant entre les rouleaux extracteurs. Ce phénomène cyclique va provoquer une modification du débit sortant de la lingotière, le métal étant encore liquide entre les peaux du métal (voir Figure 17).
Figure 17 - Illustration du gonflement
Pour estimer la variation de débit, on va se placer dans la phase statique ( vue de la Figure 17). On supposera que :
  • la déformation de la peau dans la zone notée entre les rouleaux et correspond à un triangle
  • de base, la distance entre les rouleaux et , notée
  • de hauteur
  • la hauteur oscille sinusoïdalement à la fréquence d'amplitude
  • l'épaisseur du bloc d'acier (suivant l'axe de rotation des rouleaux) est noté
  • la section du bloc d'acier est toujours notée
  • la vitesse d'avance de l'acier est notée .
Q35. Déterminer le débit sortant en l'absence de perturbations en fonction de et . Déterminer la variation de débit sortant due au gonflement de la zone noté .
Q36. En notant la variation totale de débit dû au gonflement des différentes zones, exprimer le débit sortant . Proposer une modification du schéma-bloc de la Figure 8 permettant de prendre en compte la perturbation due au gonflement.
L'objectif est maintenant de proposer une méthode pour estimer les perturbations afin de les prendre en compte correctement dans le modèle. La méthode d'estimation est une nouvelle fois basée sur les modèles d'état. Seule la perturbation due au bouchage est prise en compte pour le moment.
Figure 18 - Modèle avec la perturbation due au bouchage
Le modèle d'état adopté est avec :
: le vecteur d'état défini par avec
  • : une variable virtuelle telle que
  • hypothèse :
  • : la matrice d'état défini par
  • : la matrice de commande
    : le vecteur de commande défini par
    Q37. Déterminer la matrice de commande .
Pour déterminer la sortie, on définit également la matrice d'observation telle que
Q38. Déterminer la matrice d'observation .
Un système est dit observable si, à partir de l'observation des entrées et des sorties du système, le vecteur d'état peut être entièrement connu, notamment les conditions initiales.
Pour cela, la condition nécessaire et suffisante d'observabilité est que la matrice d'observa-
bilité soit de rang , avec la dimension du vecteur d'état .
Q39. Calculer la matrice d'observabilité et conclure quant à l'observabilité du système.
Si le système est observable, alors on peut définir un observateur dynamique qui permet
d'estimer l'état du système :
En notant :
  • le vecteur d'état estimé,
  • le vecteur des sorties du système, ici ,
  • le vecteur des sorties estimées,
  • la matrice de gain, choisie pour assurer la convergence de vers .
En combinant les deux équations, on obtient : .
Cette technique permet ainsi d'obtenir en fonction du niveau d'acier dans la lingotière, de la position de la quenouille et de la vitesse de déplacement de l'acier à la sortie de lingotière (voir Figure 19).
Figure 19 - Entrées-sortie de l'observateur de bouchage
Cette technique a été utilisée pour réaliser les simulations numériques suivantes sur le système en boucle fermée. La niveau de consigne d'acier est de -100 mm et l'asservissement est réalisé par un correcteur de type PID.
Les résultats de la simulation sont donnés sur la Figure 20. Sur la Figure 20(a), les entrées de l'observateur : le niveau de l'acier et la position de la quenouille . Sur la Figure 20(b), les sorties :
  • le débit de bouchage réel
  • le débit de bouchage estimé
  • le débit maximum théorique en sortie de busette
  • le débit entrant dans la lingotière réelle .
Figure 20 - Résultats de simulation illustrant les variations des entrées et des sorties de l'observateur de bouchage
Q40. Commenter les courbes de la Figure 20 et conclure quant à la qualité de l'observateur.
Dans cette partie, une démarche d'estimation des perturbations a été mise en place. En pratique, l'effet des perturbations n'étant pas découplable, il faudrait réaliser la même démarche en prenant en compte toutes les perturbations en même temps. Cela permettrait d'obtenir un observateur global qui estimerait chacune des perturbations en même temps (voir Figure 21).
Dans la partie suivante, cet estimateur global va être utilisé afin de compenser les perturbations dans le modèle, c'est-à-dire prédire la perturbation et modifier la commande en conséquence pour en atténuer son effet.
Figure 21 - Entrées-sortie de l'observateur de bouchage

IV Compensation des perturbations

La plupart des asservissements implantés sur les machines actuelles ont été calculés sans prendre en compte les sources de perturbations de manière précise. Dans cette partie, nous allons voir comment modifier localement le processus afin de prendre en compte l'effet des perturbations afin de les compenser et d'en atténuer leurs effets sur le niveau d'acier dans la lingotière.
Le principe de la compensation est illustré uniquement sur le gonflement de la barre d'acier au contact des rouleaux extracteurs en sortie de lingotière, mais le principe est identique avec toutes les perturbations.
Le principe de la compensation de la perturbation de gonflement est illustré sur la Figure 22. À partir des données du processus, l'observateur peut reconstruire à chaque instant le débit de gonflement estimé et sa dérivée . Ces valeurs vont permettre d'estimer à l'aide d'un prédicateur la valeur suivante du débit de gonflement et sa dérivée . Et enfin en inversant la relation de transfert de l'actionneur, on détermine la valeur à imposer à l'actionneur pour corriger le débit entrant à l'instant .
Figure 22 - Principe de la compensation de perturbations

IV. 1 Prédiction du gonflement

Comme nous l'avons vu dans la partie précédente pour l'observateur de bouchage, nous pouvons écrire l'équation d'observabilité pour l'observateur global (les variables sont indicées
avec un g sans italique) :
et
Q41. Montrer que la solution de l'équation différentielle précédente à l'instant peut se mettre sous la forme :
avec . Exprimer en fonction des différentes données.
Indication : la résolution de l'équation différentielle matricielle peut être réalisée de manière similaire à la résolution classique des équations différentielles ordinaires, vous ferez notamment apparaître des exponentielles de matrices qui seront explicitées dans la question suivante à l'aide d'un développement en série entière.
peut être calculée une fois pour toute avec un logiciel de simulation, mais doit être calculé en temps réel pendant l'asservissement.
Q42. En supposant que est constant sur l'intervalle d'intégration à , déterminer la solution sous la forme d'un développement en série entière.
Notons la matrice d'observation du débit , c'est-à-dire que . De même, notons la matrice d'observation de la variation du débit .
Q43. Déterminer les relations permettant d'obtenir et .

IV. 2 Inversion de l'équation de l'actionneur et de la lingotière

À partir du débit de gonflement estimé en , on souhaite déterminer quelle doit être la commande à imposer à l'actionneur pour modifier en conséquence le débit entrant dans la lingotière.
On rappelle que la fonction de transfert reliant le débit entrant dans la lingotière à la consigne de l'actionneur (en l'absence de ) est :
Q44. Déterminer la relation dans le domaine temporel permettant de déterminer la modification de consigne en fonction de et en travaillant par analogie avec la fonction de transfert précédente.
Figure 23 - Résultats expérimentaux d'un essai de coulée avec compensation activée ou non.

IV. 3 Validation expérimentale

L'essai expérimental se déroule dans les conditions suivantes :
  • la consigne de niveau est de 200 mm ,
  • de 0 à 200 s , l'écoulement se met en place et atteint son régime permanent à ,
  • la procédure de compensation n'est activée que pendant les intervalles 400 à 600 s et 800 à 1000 s .
    La Figure 23(a) donne l'évolution du niveau dans la lingotière , la Figure 23(b) donne l'évolution de la consigne de compensation , la Figure 23(c) donne l'évolution de la position de la quenouille et la Figure 23(d) donne l'évolution du débit estimé . Q45. Commenter ces différentes courbes à partir de .

V Synthèse

Q46. Réexpliquer en 10 à 15 lignes la démarche de modélisation qui a été adoptée dans le sujet ainsi que les résultats clés qui ont été mis en évidence.
X ENS Modélisation PSI 2014 - Version Web LaTeX | WikiPrépa | WikiPrépa