Commande d'une suspension magnétique : Différence entre versions

Semaine 2(02.01-02.07)

Un objet magnétique est soumis à l'attraction d'un électroaimant dans lequel on contrôle l'intensité, I. Un capteur de position optique permet de mesurer la distance séparant l'objet de électroaimant. Ainsi on peut agir sur la distance mesurée, x, entre cet objet et électroaimant via le courant I, la force d'attraction étant une fonction croissante de l'intensité. La commande du système consiste à contrôler l'intensité, I, en fonction de la distance x.

La distance séparant l'objet de électroaimant est notée x et l'axe x est dirigé positivement " vers le bas "
Le capteur de position associé à son électronique fourni le signal Vx=-a.x ( Vx diminue si x augmente ) donc Vx augmente si x diminue c'est à dire si l'objet monte. Vx représente donc l'altitude de l'objet.
La bobine de électroaimant possède la résistance R et l'inductance L. On négligera l'influence de la position de l'objet ( donc de x ) sur la valeur L.
La masse de l'objet suspendu est M.
L'accélération de la pesanteur est g.
Le courant dans la bobine est I.
La force d'attraction créée par le courant I sur l'objet est F=F(I,x).
La force d'attraction dans la direction x est donnée par la relation : F = − dW où W dx est l'énergie emmagasinée par le système soit:W=1／2＊L* I^2 si le matériau magnétique n'est pas saturé.
Soit:F=1/2*I^2*(−dL/dx)

Le contrôle de l'altitude de l'objet s'effectue en contrôlant la force d'attraction c'est à dire le courant, pour ce faire nous modifions la tension (Vk) sur la bobine de électroaimant. Il apparaît donc plus judicieux de régler directement le courant dans la bobine, ainsi la commande sera plus aisée en asservissant le courant plutôt qu'en agissant sur la tension d'alimentation de électroaimant. Nous allons donc étudier le cas où le courant est asservi.
Nous utiliserons donc toujours une commande en courant de électroaimant plutôt qu'une commande en tension car ainsi le modèle est seulement du second ordre au lieu du troisième ce qui conduira à une commande plus aisée et à des performances dynamiques meilleures.

Semaine 3(02.08-02.12)

A L'aide du document"Mise en Equation de la suspension magnétique.pdf", nous avons bien compris ce système, Ce système est divisé en 3 parties. Pour la partie électro-aimant, on injecte d'abord une tension qui va surement créer un courant électrique, le champ magnétique ainsi créé provoque une force qui soutient le ballon dans la position désirée.

selon cette image,La masse est soumise à deux forces opposées :

Semaine 4(02.22-02.28)

Semaine 5(02.29-03.06)

Cette semaine, nous avons essayé de utiliser la correcteur PID pour stabiliser notre système. Il faut déterminer les paramètres, nous avons cherché les informations sur Internet, nous avons plusieurs méthodes pour déterminer ces paramètres.

    
%opts = pidtuneOptions('CrossoverFrequency',32,'PhaseMargin',90);  
    
%[C, info] = pidtune(T, 'pid', opts)

Grâce aux formules dessus, nous trouvons les valeur kp et wi. Selon la théorie, nous calculons les valeurs de kp, ki, kd qui sont presque parfait:

Semaine 6(03.07-03.13)

Cette semaine, nous avons essayé d'utiliser le simulink avec les paramètres que nous avons trouvé, mais nous n’avons pas arrivé au bon résultat, notre système est encore instable. Car nous nous sommes trompé le pidtune et le pid. Donc, le professeur nous conseille d'utiliser l’autre méthode (retour d’état) pour stabiliser notre système.

Avec k1 = 0.24, te = L/R = 0.001s
On néglige la gravité mg:
Avec tm = m/b, k2 = -a/b=-25.8, tm = 0.05
On pose:

Semaine 7(03.14-03.20)

On peut trouver les poles du sys p1 = -1/te, p2 = -1/tm, p3 = 1/tm;
Comme p3 est positive, donc le système est instable, on doit déplacer les pôles p2 et p3.
On pose:
Z= 0.707, wn = 35
On utilise le matlab pour trouver les deux pole, p2 = -24.7450+24.7525i P3 = -24.7450-24.7525i.

Et utilise ce schéma pour tester


On peut trouver que notre système est bien stable et correspondant, mais il y a aussi un peu de dépassement.

Semaine 8(03.21-03.27)

On attaque l'obeservateur du systeme cette semaine.
Un conservateur est pour le but deobserver les états du système que nous pouvons pas les mesurer directement dans le système.
Par exemple, la variation du courant de la machine et de la vitesse du ballon ne peuvent pas être mesuré. C’est la raison pour laquelle nous avons besoin d’un observateur.

L’observateur dans notre schéma bloc est effectivement une copine du système réel. Cependant, les variable sont pas les valeurs en temps réel. Il faut donc bien choisir les paramètres de la matrice L pour que les erreurs entre les valeur estimés et les valeurs réels soient plus petites.
Et voici ce que nous avons observé sur la variation du courant, de la vitesse et de la position

Semaine 9(03.28-04.03)

Nous essayons une méthode avancement du phase pour ce système.
En considérant un système dont la fonction de transfert en boucle ouverte est H(jw). En boucle fermée, il aura une fonction de transfert : H(jw)/(1+H(jw))
Comme notre système est 3ème ordre, et instable.
Nous pouvons obtenir le Bode diagram suivant:

Quand magnitude(db) = 0db, et le phase > -180, le système est stable. Cependant, pour notre système, quand le magnitude = 0db, le phase < 180, donc , notre système est instable.
Donc, nous utilisons la méthode avance de phase afin de stabiliser notre système. Nous désirons apporter un déphasage de 30°, alors -181°+30° = -151° > -180°. Notre système deviendras stable.

Nous pouvons utiliser ces formules pour calculer a, et en plus Wc = 1/(ta). Nous pouvons alors avoir le t.
Quand les valeurs de t et a sont terminés, le correcteur avance phase est fixé.
Nous les appliquent sur matlab.
Nous avons alors un nouveau Bode diagram:

notre système est maintenant stable grâce à .

Semaine 10(04.18-04.24)

Le prof nous a posé la question sur les équations qui définissent le système de la suspension magnifique.
Donc cette semaine nous avons utilisé deux méthodes différents a fin de vérifier le système.
le but est de vérifier si nous pouvons définir le même système avec deux méthodes différents.



on remarque bien que les deux méthodes définissent le même système.

Semaine 11(04.24-04.30)

Nous attaquons alors la partie manipulation de notre machine.
C'est la partie plus important mais aussi plus difficile.
Dans la partie de la simulation, nous avons toujours un bon résultat, puis que nous simulons dans un système bien linéaire.
Nous avons rencontré des gros soucis au début puisque nous avons essayé plusieurs façons à simuler le système réel (les équations d'état) sous matlab.
Il sont tous une approximation du système réel mais certain méthode n'est pas assez précis à cause de certains négligences ou approximation de notre calculation.


Nous finalement fixons à la méthode de troisième ordre avec les équations calculé par la fonction transfert après avoir exécuté plusieurs essais.

En jouent sur le temps d'échantillonnage, l'amortissement, les pôles de l'observateurs, nous avons à la fin réussi à stabiliser notre système.

rapport