Fermer X

Contents

On considère dans un premier temps un moteur linéaire de masse M, se déplaçant le long d’un rail rectiligne, relié à un corps de masse m par une tige flexible de masse négligeable que l’on peut assimiler à un ensemble ressort-amortisseur, le ressort ayant une raideur k et l’amortisseur un coefficient d’amortissement r. L’abscisse du moteur (resp. de la masse) est notée x (resp. z). Le moteur produit une force F qui est la variable de commande. Le système s’écrit, en première approximation,
M  ¨x = F −  k(x − z) − r(˙x − ˙z)
m z¨=  k(x − z) + r(x˙−  ˙z)
(1)


PIC


Figure 1: Simple masse reliée élastiquement à un moteur linéaire.

On désire générer des déplacements arrêt-arrêt rapides de l’ensemble masse-moteur pour positionner la masse de façon précise et par conséquent sans oscillation. Notons que dès que la force F n’est pas nulle, le moteur se déplace et, à cause de la tige flexible, engendre une dynamique oscillante de la masse.

Dans ce qui suit, nous prenons M = 4 kg pour la masse du moteur, m = 0.15 kg pour la masse transportée, k = m(4 2π)2 N/m pour la raideur de la tige (qui correspond à une fréquence propre de 4 Hz), r = 2 0.01√km--- Ns/m pour le coefficient d’amortissement de la tige (qui correspond à une constante de temps d’amortissement de -1-
0.01 s).

1ère partie : Étude du déplacement sans tenir compte de la masse transportée

Dans un premier temps, on veut effectuer le déplacement sans tenir compte de la masse transportée. Pour cela, on fait comme si le système était donné par

M  ¨x = F
les contributions de la masse transportée étant considérées comme des perturbations.

Ainsi, pour un déplacement arrêt-arrêt partant de x0 au repos à l’instant 0 et arrivant à x1 au repos à l’instant T, avec la force et sa dérivée nulles aux instants 0 et T, on peut prendre comme trajectoire de référence

                              (                                      )
                        ( t )4          ( t )      ( t )2      ( t)3
xref(t) = x0 + (x1 − x0 ) --    35 − 84   --  + 70   --   − 20  --
                          T               T          T          T
(2)

On vérifie alors que la référence de la force Fref est donnée par

                       (         ) (  ) (       (  )     (   )2     (   )3)
                         x1-−-x0-    t-           t-       -t         t-
Fref = M  ¨xref = 420 M      T 2       T    1 − 4   T   + 5  T    − 2   T      .

Pour rendre le déplacement du moteur insensible aux perturbations, on ajoute à la force de référence une boucle PID :

                                             ∫ t
F  = Fref − kP (x − xref) − kD (x˙− x˙ref) − kI   (x (τ) − xref(τ))d τ
                                              0
(3)

Question 1 Exécuter le fichier Smasse0.sci. Pour un déplacement de x0 = 0 à x1 = 0.1 m, régler la durée T telle que le maximum du module de la force de référence ne dépasse pas 1.5g.

Question 2 Régler les gains kP , kD et kI du PID pour que l’influence de la masse transportée soit presque insensible sur la position du moteur. Quelle amplitude d’oscillations de la masse transportée obtient-on ?

Question 3 Peut-on modifier la fréquence d’oscillation et la constante de temps d’amortissement de la masse transportée ?

2ème partie : Étude du déplacement en tenant compte de la masse transportée

On veut essayer d’empêcher les oscillations de la masse sans toutefois observer sa position. Pour cela, on va utiliser le modèle (1) pour chercher des trajectoires de consigne de déplacement du moteur et de la force telles que le moteur et la masse transportée arrivent au point final au repos à l’instant T.

Question 4 Exprimer x, z et F en fonction d’une sortie plate du système (1) et de ses dérivées.

Question 5 En déduire une trajectoire arrêt-arrêt pour le moteur et la masse transportée partant de x0 = z0 = 0 à l’instant t = 0 à x1 = z1 = 0.1 m à l’instant t = T calculé précédemment.

Question 6 Exécuter le fichier Smasse1.sci. Régler les gains kP , kD et kI du PID, pour assurer la stabilité et la robustesse d’erreur de suivi dans le cas d’une erreur sur les conditions intiales de l’ordre de 5% sur z, ou sur les masses, raideurs et amortissements de l’ordre de 0.01%.

Question 7 En reprenant la méthode de la première partie, régler la durée T1 pour que la loi de commande (3) avec (2), où T est remplacée par T1, produise des oscillations de la masse transportée de même amplitude qu’aux questions 2.2 et 2.3 avec T.

3ème partie : Même étude pour deux masses transportées

Question 8 Reprendre l’étude précédente en ajoutant une seconde masse m= 0.04 kg reliée au moteur par une seconde tige flexible dont les caractéristiques sont k= m(7 2π)2 N/m (fréquence propre de 7 Hz) et r= 2 0.012√ -′--′
  k m Ns/m (constante de temps d’amortissement de  1
0.012 s). L’étude du déplacement du moteur sans tenir compte des masses transportées se fait à l’aide du fichier Dmasses0.sci, et celle tenant compte des masses transportées à l’aide du fichier Dmasses1.sci. On pourra modifier T et les gains du PID en faisant exécuter les fichier Scilab.


PIC


Figure 2: Deux masses reliées élastiquement à un moteur linéaire.

Dans ce cas, le système est donné par

M ¨x = F −  k(x − z) − r(˙x − ˙z) − k′(x − z′) − r′(˙x − z˙′)

mz′¨=′  k(x′ − z) +′ r(˙x −′ ˙z) ′
m ¨z =  k(x − z ) + r(x˙− z˙).
(4)

On pourra se servir des formules

        (       ′)     (        ′     ′)     (    ′     ′)           ′
x = y +   r-+  r-  ˙y +   m-+  m--+ -rr-  ¨y +   mr--+-rm--  y(3) + mm---y(4)
        ( k    k′)     ( k    k′ ) kk ′           kk ′            kk′
          r-   r′       m-′   rr′-      rm-′ (3)
z = y +   k +  k′  ˙y +   k′ + kk ′ y¨+  kk′ y
         ( r   r′)     ( m    rr′)      mr ′
z′ = y +   --+ -′  y˙+   -- + ---′ y¨+  --′ y (3)
           k   k         k    kk        kk
et, pour générer des trajectoires arrêt-arrêt :
                            (
                     ( t )7              (  t)          ( t)2          ( t )3
y(t) = x0 + (x1 − x0 )  --    1716 −  9009  --  +  20020   --   − 24024   --
                       T                   T              T              T)
                                         ( t)4         ( t)5        ( t)6
                                +  16380   --   − 6006   --   + 924  --     .
                                           T             T           T
L'École des Ponts ParisTech est membre fondateur de

L'École des Ponts ParisTech est certifiée 

ParisTech ParisEst ParisTech

 

Copyright 2014
École des Ponts ParisTech
All rights reserved