M1 CSy – module P8

PROJET DE SIMULATION AVEC MATLAB

∗
Commande de la position angulaire d’une antenne
–
Christophe Calmettes & Jean-José Orteu

On considère le système représenté sur la figure 1 et constitué :

– d’une antenne parabolique dont on souhaite asservir la position angulaire ;

– d’un moteur permettant, grâce à un jeu d’engrenages, de faire pivoter l’antenne ;
– d’un module comprenant la partie différentielle (comparateur) et la partie puissance du système ;
– d’un potentiomètre de sortie (capteur) et d’un potentiomètre d’entrée (consigne).

Fig. 1 – Le système étudié – éléments constitutifs

Lors de ce travail, les questions devant être traitées à l’aide du logiciel Matlab sont repérées par le
sigle M .

I Du système au modèle
Le schéma de principe de l’asservissement est représenté sur la figure 2. La notation des différents
signaux mis en jeu est la suivante :
∗
sujet inspiré de l’ouvrage de Norman S. Nise “Control Systems Engineering”, Sixth edition, WILEY.

1
– sortie de l’asservissement (variable à commander) exprimée en radians : θo (t) ;
– signal de commande (entrée de l’amplificateur de puissance) exprimé en volts : vp (t) ;
– consigne (signal de référence) exprimé en radians : θi (t) ;
– on notera respectivement θm (t) et ωm (t) la position et la vitesse angulaire du moteur.

Fig. 2 – Schéma de principe

Les potentiomètres d’entrée et de sortie de l’asservissement sont supposés linéaires et de gain Kpot .
L’amplificateur couplé au comparateur est également linéaire et de gain K variable ; il pourra être
ajusté pour régler la composante proportionnelle de la commande.
L’amplificateur de puissance se comporte comme un système du premier ordre de fonction de
transfert :
K1
p+a

L’ensemble moteur-charge est modélisé par un premier ordre avec intégrateur ; on note
Km
p(p + am )

sa fonction de transfert.
Enfin, Kg représente le gain de réduction du jeu d’engrenages.

Les valeurs numériques des paramètres sont reportées dans le tableau 1 (sans unité).

1. A partir du schéma de principe (figure 2), dessiner un schéma-blocs détaillé de l’asservisse-

ment : y faire apparaı̂tre toutes les fonctions de transfert du système ainsi que les différents
signaux mis en jeu.
Le module differential preamplifier sera représenté par un comparateur couplé avec un cor-
recteur proportionnel de gain K variable.

2
 Kpot 0.318
K variable
K1 125
a 100
Km 2.083
am 1.71
Kg 0.1

Tab. 1 – Paramètres de l’asservissement

II Analyse en boucle ouverte (de l’ensemble amplification-moteur-

charge)
On s’intéresse au comportement dynamique de l’ensemble amplification-moteur-charge lorsqu’on
applique une tension sur le signal de commande vp (t) ; le signal d’intérêt est ici la vitesse angulaire
du moteur ωm (t)∗ .

2. Ecrire la fonction de transfert M (p) = ΩVm (p)

p (p)
.
Quel est son ordre ?
Donner la valeur des zéros, des pôles et du gain statique de M (p).

3. M Afficher les pôles et les zéros dans le plan complexe.

4. L’ordre de ce modèle est-il raisonnablement simplifiable ? Justifier.

Dans l’ensemble amplification-moteur-charge, par quelle fonction de transfert l’amplificateur
de puissance peut-il être remplacé ?

5. Calculer l’expression de ωm (t) du modèle simplifié lorsque vp (t) est un échelon unitaire.

6. M Tracer l’allure de la réponse indicielle du modèle M (p). La comparer à celle du modèle

simplifié.

7. Justifier et commenter d’un point de vue physique l’allure obtenue pour ωm (t).

8. M Tracer l’allure de la position angulaire θm (t) du procédé (modèle non simplifié) soumis à
un échelon. Commenter le résultat obtenu.

III Asservissement
Cette partie concerne l’analyse de l’asservissement de la figure 2 ; on note θi (t) le signal de consigne.

III.1 Schémas-blocs et fonctions de transfert

9. Dessiner le schéma-blocs de l’asservissement en considérant la correction proportionnelle de
gain K (gain variable du pré-amplificateur). Spécifier le nom et la nature des différents signaux
∗
on notera Ωm (p) la transformée de Laplace de ωm (t)

3
 mis en jeu.

10. Sans calcul, quelle projection peut-on faire sur les performances en précision (erreurs de
position, de vitesse et d’accélération) de cet asservissement ? Justifier.

11. Ecrire la fonction de transfert T (p) en boucle ouverte et la fonction de transfert F (p) en
boucle fermée.

On introduit Ts (p) et Fs (p), fonctions de transfert obtenues lorsqu’on applique la simplification vue
à la question 4 du paragraphe II .

12. M Définir T (p), Ts (p), F (p) et Fs (p) sous Matlab pour K = 1.

III.2 Analyse temporelle du modèle simplifié

On considère, dans ce paragraphe, le modèle simplifié du système.

13. Mettre Fs (p) sous forme canonique et en extraire les paramètres représentatifs. Les exprimer
en fonction des paramètres du système.

14. Donner la valeur de ces paramètres pour K = 800.

En déduire une estimation des caractéristiques dynamiques de l’asservissement (premier
dépassement, temps de montée et temps de réponse à 5%) attendues dans le cas d’une réponse
indicielle.

15. M Retrouver ces résultats à l’aide de Matlab.

16. Quelle tension est délivrée par le potentiomètre d’entrée pour une consigne de 1 tour ?

17. M Pour K = 800, tracer sur un même graphe les réponses indicielles de l’asservissement
simplifié et de l’asservissement non simplifié. Commenter.

18. Pour quelle valeur de K a-t-on un premier dépassement de 15%.

Justifier votre réponse en basant votre démarche sur le modèle simplifié ?

19. M Simuler et commenter.

III.3 Analyse de l’asservissement dans le lieu des racines pour une correction
proportionnelle
On considère l’asservissement non simplifié de fonctions de transfert T (p) (boucle ouverte) et F (p)
(boucle fermée).

20. En utilisant le critère de Routh, déterminer la condition sur K pour que l’asservissement soit
stable.

21. M Dessiner le lieu des racines.

22. Quels est (sont) le(s) pôle(s) dominant(s) ? Justifier.

4
 23. M Quelles valeurs Kc1 et Kc2 de K (Kc1 < Kc2 ) correspondent à des changements notables
dans la dynamique de cet asservissement ?

24. Compléter (pour les caractéristiques dynamiques, donner le sens de variation en le justifiant)
et commenter le tableau suivant :

K 0 Kc1 Kc2 +∞

comportement dynamique

temps de réponse

temps de montée

premier dépassement

25. M A partir du lieu des racines en donnant la démarche :

– déterminer la valeur du gain Kc3 pour lequel l’amortissement égale 0.7. Estimer le temps
de montée et le temps de réponse attendus ;
– déterminer la valeur du gain Kc4 pour lequel le système a un temps de réponse de 4 secondes.
Estimer le temps de montée et le premier dépassement attendus ;
– déterminer la valeur du gain Kc5 pour lequel le système a un temps de montée de 0.2
secondes. Estimer le temps de réponse et le premier dépassement attendus.

26. M Tracer la réponse indicielle pour ces trois cas et commenter chacun d’eux.

III.4 Analyse de l’asservissement dans le plan de Bode

On considère l’asservissement non simplifié de fonctions de transfert T (p) et F (p).

27. M Quel lieu de transfert permet de faire l’analyse de la stabilité d’un asservissement par le
critère géométrique des marges ?
Le tracer dans le plan de Bode et dans le plan de Nyquist pour K = 1.

28. A partir du lieu dans le plan de Bode, mesurer la marge de gain. Que peut-on déduire de
cette valeur ?
Retrouver ce résultat dans le plan de Nyquist.

29. M On considère les valeurs de gain Kci (i = 1, . . . , 5) de la partie III.3 et on note :

– Mφ la marge de phase et MG la marge de gain ;

– ξ le coefficient d’amortissement et ωn la pulsation naturelle du ou des pôles dominants ;
– D1 le premier dépassement et tm le temps de montée, mesurés sur la réponse indicielle de
l’asservissement ;
– BP et Q la bande passante et le coefficient de surtension (résonance) mesurés en boucle
fermée (le tracé des lieux de transfert en boucle fermée sera donné dans le plan de Bode).

5
 K Kc1 = Kc3 = Kc4 = Kc5 = Kc2 = ր

Mφ

D1

ωn

BP

tm

Remplir le tableau avec les valeurs mesurées demandées et donner, dans la dernière colonne,
le sens de variation de chaque paramètre lorsque K augmente.
Discuter des paramètres liés au degré de stabilité de l’asservissement d’une part, et de ceux
liés à la montée en régime (raideur) d’autre part.

30. M Quelle marge de gain a-t-on sur le modèle simplifié ? Ce résultat était-il prévisible ?

III.5 Performances en précision

31. Calculer l’erreur de position et l’erreur de vitesse (traı̂nage).

32. Pour quelle valeur de gain Kc6 a-t-on une erreur de vitesse de 4% ?

33. M En vous aidant du lieu des racines du paragraphe III.3, déduire, pour K = Kc6 , les
caractéristiques dynamiques de la réponse indicielle.

34. M Tracer, pour cette valeur, la réponse indicielle et la réponse à une rampe de l’asservisse-
ment.
Y porter les grandeurs remarquables.

IV Synthèse de lois de commande

Le gain variable K du pré-amplificateur constituera, pour la synthèse, la composante proportionnelle
des différents correcteurs.

IV.1 Synthèse par approche fréquentielle dans le plan de Bode

Dans cette partie, la synthèse des correcteurs D1 (p), D2 (p) et D3 (p) sera effectuée avec le modèle
non simplifié dans le plan de Bode. Le choix du correcteur et la démarche d’obtention de ses
paramètres doivent être clairement explicités.

6
Correcteur D1 (p)
Le cahier des charges est le suivant :

3 temps de montée compris entre 0.2 et 0.4 s ;

3 premier dépassement de 5%.

35. Quel type de correcteur permet de répondre au cahier des charges ? Justifier.

36. Faire la synthèse du correcteur D1 (p).

37. M Simuler le comportement de l’asservissement corrigé. Commenter.

Correcteur D2 (p)
Le cahier des charges est le suivant :

3 erreur de vitesse nulle ;

3 premier dépassement de 5%.

38. Quel type de correcteur permet de répondre au cahier des charges ? Justifier.

39. Faire la synthèse du correcteur D2 (p).

40. M Simuler le comportement de l’asservissement corrigé. Commenter.

Correcteur D3 (p)
Le cahier des charges est le suivant :

3 erreur de vitesse nulle ;

3 premier dépassement de 5% ;
3 temps de montée compris entre 0.2 et 0.4 s.

41. Quel type de correcteur permet de répondre au cahier des charges ? Justifier.

42. Faire la synthèse du correcteur D3 (p).

43. M Simuler le comportement de l’asservissement corrigé. Commenter.

IV.2 Méthode du modèle – Placement de pôles

Le modèle de départ est ici le modèle simplifié.
On considère le correcteur D4 (p) :
p+α
D4 (p) = K
p+β

7
On souhaite que l’asservissement corrigé se comporte comme un système du second ordre de pôles
−2 ± 1.6j.

44. Quelles caractéristiques dynamiques sont associées à ces pôles ?

45. Calculer les paramètres du correcteur qui répond au cahier des charges.

46. M Simuler le comportement de l’asservissement corrigé. Commenter.

47. M Vérifier que le correcteur D4 (p) remplit également sa fonction dans le cas du modèle non
simplifié.

IV.3 Commande Numérique

On souhaite, pour des raisons de flexibilité, commander le système par un calculateur. On considère
le modèle non simplifié et la période d’échantillonnage notée Te est, dans un premier temps, imposée
et fixée à 25ms.

48. Dessiner le schéma-blocs de l’asservissement dans ce cas.

49. M Calculer la fonction de transfert échantillonnée de l’ensemble bloqueur + amplificateur de

puissance + moteur.

Correction Proportionnelle
On considère ici les trois cas Kc3 , Kc4 et Kc5 du paragraphe III.3.

50. M Définir, pour chaque cas, la fonction de transfert en boucle fermée.

51. M Tracer pour chaque cas, la réponse indicielle de l’asservissement. Comparer qualitati-
vement les allures obtenues dans le cas numérique avec celles des réponses indicielles de
l’asservissement analogique de la partie III.3.
Conclure sur le choix de la période d’échantillonnage.

52. M Définir la fonction de transfert en boucle ouverte dans le cas de l’asservissement échantil-
lonné.
A l’aide du lieu des racines, retrouver les nouvelles valeurs de Kc1 et de Kc2 correspondant à
des changements notables dans la dynamique de l’asservissement (l’outil grid disponible dans
la figure Matlab est adapté à ce travail).

53. M Que deviennent ces deux valeurs lorsque la période d’échantillonnage est doublée ? Lors-
qu’elle est multipliée par 4 ? Conclure.

8
Discrétisation de correcteur continu
On souhaite mettre en œuvre sous forme numérique le correcteur analogique calculé à la question
IV.2.

54. M Calculer, par la méthode de Tustin, la fonction de transfert équivalente au correcteur

D4 (p).

55. M Tracer la réponse indicielle de l’asservissement.

Comparer qualitativement les caractéristiques dynamiques obtenues avec celles de l’asservis-
sement analogique de la question IV.2.

56. M Analyser l’effet d’une augmentation de la période d’échantillonnage sur les caractéristiques
de l’asservissement.

Synthèse directe d’un correcteur numérique par placement de pôles

La synthèse est réalisée à partir du modèle simplifié. On se propose de régler un régulateur avance
de phase numérique de structure :
z−A
D5 (z) = K
z−B
On souhaite que le système en boucle fermé se comporte comme un système du second ordre avec
un facteur d’amortissement de 0.7 et un temps de pic (instant du premier dépassement) de 10 Te .

57. M Calculer la fonction de transfert échantillonnée de l’ensemble bloqueur + amplificateur de

puissance + moteur simplifié.

58. M Calculer ses pôles.

59. Choisir le paramètre A de façon à compenser un des deux pôles.

60. Calculer les pôles de la fonction de transfert en boucle fermée qui satisfont au cahier des
charges.

61. En déduire le réglage des paramètres K et B.

62. M Tracer la réponse indicielle de l’asservissement. Commenter.