Modélisation et Analyse des Systèmes Asservis
Modélisation et Analyse des Systèmes Asservis
Question 2 : Déterminer sous forme canonique, les fonctions de transfert des systèmes modélisés par les équations
différentielles ci-dessous. Puis dessiner ces 3 fonctions dans 3 blocs en séries.
dθ(t)
6) ωs (t) = r ωe (t) 7) ω(t) = 8) τω̇ m (t) + ωm (t) = K um (t)
dt
Des essais sont réalisés sur un système de portes en verre rétractables permettant de filtrer
et sécuriser, à l’aide de badge, l’accès aux bureaux d’une administration ou d’une
entreprise.
Le système, encore en phase de validation, est soumis à différents essais correspondant à
différents réglages.
Pour chacun d’entre eux, l’évolution de la distance parcourue latéralement par la porte par
rapport à la position de départ 𝑥𝑐 (0− ) = 0 𝑚 est mesurée en fonction du temps en
seconde.
Pour chacun de ces essais, la consigne de déplacement latéral imposée est de 𝑥𝑐 (0+ ) =
𝑥𝑐0 = 1 𝑚.
Les réponses 𝑥(𝑡) obtenus sont présentés ci-après.
Question 1 : Évaluer, dans chacun des cas et en utilisant les critères proposés, les performances du système de portes
rétractables.
Question 2 : Donner le réglage qui permet d’avoir le système le plus précis. Donner aussi celui qui permet d’avoir le
système le plus rapide.
Question 3 : Indiquer les risques liés au fait d’avoir une valeur du 1 er dépassement trop élevée sur un tel système.
Temps de réponse à 5% :
Erreur statique :
Réglage N°2
1er dépassement absolu :
Temps de réponse à 5% :
Erreur statique :
Réglage N°3
1er dépassement absolu :
Temps de réponse à 5% :
Erreur statique :
Chariot x(t)
ωv(t)
θv(t)
Capteur
ωm(t)
Interface Homme- Calculateur
Machine
Moteur
électrique Réducteur
poulie/courroie
xc(t) uc(t) um(t)
uv(t)
Variateur
On suppose que toutes les conditions initiales sont nulles, et que toutes les constantes sont positives.
Question 2 : Déterminer la fonction de transfert d’une intégration.
Question 3 : Déterminer la fonction de transfert du système vis-écrou.
Question 4 : À l’aide de la description ci-dessus, déterminer les fonctions de transfert de chaque constituant (sauf l’IHM),
puis compléter le second schéma-bloc page suivante en faisant apparaitre ces dernières à l’intérieur des
blocs ainsi que les flux transmis d’un bloc à l’autre.
Question 5 : Déterminer la fonction de transfert de cet asservissement (Ne pas remplacer 𝐾𝐼𝐻𝑀 ) et la mettre sous forme
𝐾𝑐𝑎𝑝 𝑝𝑎𝑠
canonique (Vous utiliserez 3 couleurs pour 𝐾𝐼𝐻𝑀 , et ).
2𝜋 2𝜋
Question 6 : Conclure sur les performances de précision. Que vaut le gain statique pour un système précis ?
Question 7 : Redéterminer la fonction de transfert de l’interface homme-machine 𝐾𝐼𝐻𝑀 à l’aide d’un raisonnement sur la
proportionnalité des flux que l’on souhaite.
Question 8 : Sur quel constituant intervenir pour améliorer la performance de rapidité ?
𝑑𝜃𝑒𝑐ℎ(𝑡)
𝜃𝑒𝑐ℎ (𝑡) + 𝜏𝑒𝑐ℎ = 𝐾𝑒𝑐ℎ 𝑞(𝑡)
𝑑𝑡
𝑑𝜃(𝑡)
− loi de transfert de chaleur dans l’enceinte avec 𝜃(𝑡) et 𝜃𝑒𝑐ℎ (𝑡) en ° : 𝜃(𝑡) + 𝜏𝑒𝑛𝑐 = 𝐾𝑒𝑛𝑐 𝜃𝑒𝑐ℎ (𝑡)
𝑑𝑡
On suppose que toutes les conditions initiales sont nulles. L’entrée du système est l’angle d’ouverture de la vanne 𝛼(𝑡) en
rad et la sortie, la température de l’enceinte 𝜃(𝑡).
Déterminer les fonctions de transfert de la vanne, de l’échangeur et de l’enceinte.
Question 1 : Dessiner la portion de schéma-bloc représentant le système et faisant intervenir uniquement les
composants précédemment définis.
Afin de réguler la température de l’enceinte, on choisit d’asservir cette température en motorisant la vanne.
Vanne Enceinte
θc(t)
Interface H/M Moteur α(t) Echangeur
Comparateur
Correcteur
Hacheur q(t) θech(t)
θ(t)
Capteur
Pompe
Un capteur de gain 𝐾𝑐𝑎𝑝 installé dans l’enceinte, mesure la température 𝜃(𝑡) en ° et la traduit en une tension 𝑢𝑚𝑒𝑠 (𝑡) en
V.
La tension 𝑢𝑚𝑒𝑠 (𝑡) en V est comparée à la tension image de consigne 𝑢𝑐 (𝑡) en V issue d’une interface homme-machine.
L’image de l’erreur 𝜀(𝑡) = 𝑢𝑐 (𝑡) − 𝑢𝑚𝑒𝑠 (𝑡) est corrigée et traduite en une tension de commande 𝑢𝑣 (𝑡) en V pour le
𝑈 (𝑝)
hacheur. La fonction de transfert du correcteur est 𝑣 = 1 + 𝜏𝑐 𝑝. Remarque 𝜏𝑐 sera choisi égal à 𝜏𝑒𝑛𝑐 .
𝜀(𝑝)
Le hacheur de gain 𝐾𝑣 fournit la tension d’alimentation du moteur 𝑢𝑚 (𝑡) en V. Le moteur agit sur l’angle d’ouverture de la
𝛼(𝑝) 𝐾 1
vanne. La fonction de transfert du moteur est = 𝑚 .
𝑈𝑚 (𝑝) 1+𝜏𝑚 𝑝 𝑝
Question 2 : Déterminer la structure du schéma-bloc modélisant cet asservissement, en identifiant les différents
composants (nom sous les blocs) et en précisant leur fonction de transfert à l’intérieur des blocs, ainsi que
les grandeurs avec leur unité transmises d’un bloc à l’autre.
Question 3 : En déduire sous forme canonique la fonction de transfert globale modélisant cet asservissement.
Question 4 : Conclure sur la performance de précision.
Exercice 8 : SCHEMAS-BLOCS
Question 1 : Déterminer les fonctions de transfert des schémas-blocs suivants. Vous noterez H au lieu de H(p). (On ne
peut pas mettre H(p) sous forme canonique car on n'a pas leur expression)
1)
E(p)
S(p)
+ H1(p) + H2(p) H3(p) - H4(p)
+
- -
H5(p)
2)
H3(p)
E(p)
S(p)
- H1(p) H2(p) +
+ +
- -
Question 2 : Déterminer les fonctions de transfert des schémas-blocs suivants. Les mettre sous forme canonique.
Déterminer le gain, l'ordre et la classe.
1)
E(p) S(p)
10 - 2 2 1
+ +
- p+3 p+1 p 100
K
E(p) Cr(p)
U(p) Ωm(p)
I(p)
- 1 - 1
+ K +
R + Lp f + Jp
+ +
Iref(p)
-
Ac +
3) + Y(p)
d
c
X(p) Z(p)
- K
+ + + G -
p +
- -
Question 3 : Déterminer les fonctions de transfert du schéma-bloc suivant. Vous noterez H au lieu de H(p).
1) Il est parfois moins cher d’acheter 2 petits moteurs, qu’un seul gros. On a alors une motorisation en parallèle :
X(p) Y(p)
H1(p) + H2(p)
-
H3(p)
H6(p)
U(p)
W(p)
-
H4(p) + H5(p)
On s'intéresse à une chaine de puissance composée d'un motoréducteur entraînant une charge.
R L
i(t)
Jm Jr
Cm (t) r
fm
u(t) e(t)
ωm (t) Jc
Cr (t)
fc
MCC
ωr (t)
réducteur
Grandeurs électriques :
Equation électrique :
Loi des mailles et loi d'Ohm :
𝑑𝑖
𝑢(𝑡) = 𝑅 𝑖(𝑡) + 𝐿 (𝑡) + 𝑒(𝑡)
𝑑𝑡
Equations de couplage électromagnétique :
Maxwell-Faraday :
𝑒(𝑡) = 𝐾𝑒 𝜔𝑚 (𝑡)
𝐶𝑚 (𝑡) = 𝐾𝑡 𝑖(𝑡)
Equation mécanique :
Principe fondamental de la dynamique appliqué à l'arbre moteur :
𝑑𝜔𝑚
𝐶𝑚 (𝑡) − 𝐶𝑟𝑒 (𝑡) = 𝐽𝑒𝑞 (𝑡) + 𝑓𝑒𝑞 𝜔𝑚 (𝑡)
𝑑𝑡
A l'instant initial :
𝑖(0− ) = 0 A
𝜔𝑚 (0− ) = 0 rad/s
Question 1 : Expliquer brièvement le nom de chacune des quatre équations et ce qu'elles représentent.
Question 2 : En précisant l'hypothèse utilisée, appliquer la transformée de Laplace aux 4 équations du moteur. On
notera que Ω est la majuscule de ω.
Question 3 : Reproduire et compléter le schéma-bloc suivant :
Cre (p)
Ωm (𝑝)
U(p)
+
-
+
-
théorème utilisé.
La fonction de type 𝐻1 (𝑝) correspond au fonctionnement de type suiveur (consigne variable) du système, alors que la
fonction 𝐻2 (𝑝) correspond au fonctionnement de type régulateur (consigne constante).
Cre (p)
U(p) Ωm (𝑝)
-
F1 (p) + F2 (p)
On vérifie que les deux fonctions ont le même dénominateur, c'est de lui que dépend la stabilité du système.
Question 10 : En suivant une démarche d'identification de 𝐻1 (𝑝) avec un deuxième ordre de classe 0 de la forme
𝐾
2𝑧
1+ 𝑝+
1
𝑝2
montrer que :
𝜔0 𝜔0 2
𝐾𝑡
𝐾=
𝐾𝑒 𝐾𝑡 +𝑅𝑓𝑒𝑞
1 𝑅𝐽𝑒𝑞 +𝐿𝑓𝑒𝑞 𝐾 𝑒𝑠𝑡 𝑙𝑒 𝑔𝑎𝑖𝑛 𝑠𝑡𝑎𝑡𝑖𝑞𝑢𝑒 en (rad/s)/V
𝑧=
2
√𝐽𝑒𝑞 L(𝐾𝑒 𝐾𝑡 +𝑅𝑓𝑒𝑞 ) 𝑧 𝑒𝑠𝑡 𝑙𝑒 𝑓𝑎𝑐𝑡𝑒𝑢𝑟 𝑑′𝑎𝑚𝑜𝑟𝑡𝑖𝑠𝑠𝑒𝑚𝑒𝑛𝑡
𝜔0 𝑒𝑠𝑡 𝑙𝑎 𝑝𝑢𝑙𝑠𝑎𝑡𝑖𝑜𝑛 𝑝𝑟𝑜𝑝𝑟𝑒 en rad/s
𝐾𝑒 𝐾𝑡 +𝑅𝑓𝑒𝑞
𝜔0 = √
{ 𝐽𝑒𝑞 L
𝐿
Soit 𝑇𝑒 = : la constante de temps électrique
𝑅
𝑅𝐽𝑒𝑞
et 𝑇𝑚 = : la constante de temps mécanique
𝐾𝑒 𝐾𝑡 +𝑅𝑓𝑒𝑞
Hypothèse : - On suppose que 𝑇𝑒 ≪ 𝑇𝑚 , et donc on néglige l'inductance L = 0 𝑚𝐻 (la tension aux bornes de
l’inductance est faible vis-à-vis des tensions électriques présentes dans le montage).
- On néglige le frottement 𝑓𝑒𝑞 = 0 Nm/(rad/s) (le couple dû aux frottements visqueux est très faible
devant le couple électromagnétique).
- Le champ magnétique généré par le stator est réalisé par un aimant permanent 𝐾𝑒 = 𝐾𝑡 = 𝐾.
Le dénominateur de la fonction de transfert peut alors se mettre sous la forme (1 + 𝑇𝑒 𝑝)(1 + 𝑇𝑚 𝑝) = 1 + (𝑇𝑚 +
𝑇𝑒 ) 𝑝 + 𝑇𝑚 𝑇𝑒 𝑝2 ≈ 1 + 𝑇𝑚 𝑝 + 𝑇𝑚 𝑇𝑒 𝑝2 ≈ 1 + 𝑇𝑚 𝑝 au vu des valeurs numériques. On approxime alors le moteur par une
fonction de transfert du 1er ordre.
1 𝑅
𝐾 𝐾2
Question 11 : En utilisant les hypothèses simplificatrices, montrer que Ω𝑚 (𝑝) = 𝑈(𝑝) − 𝐶𝑟𝑒 (𝑝)
1+𝑇𝑚 𝑝 1+𝑇𝑚 𝑝
Question 12 : Donner les unités de 𝐾𝑒 et 𝐾𝑡 . Montrer que 𝐾𝑒 et 𝐾𝑡 ont la même unité, pour ce faire, on pourra écrire les
unités d’une puissance électrique, puis les unités d’une puissance mécanique de rotation.
On soumet le moteur à courant continu PARVEX RX 630 E à un échelon de 𝑈0 = 24 𝑉. On néglige la perturbation
𝐶𝑟𝑒 (𝑝) = 0
ainsi
𝐾𝑚 𝐾𝑚
Ω𝑚 (𝑝) = 𝑈(𝑝) ≈ 𝑈(𝑝)
(1 + 𝑇𝑒 𝑝)(1 + 𝑇𝑚 𝑝) 1 + 𝑇𝑚 𝑝
𝐾𝑒 = 52 𝑉/1000 𝑡𝑟/𝑚𝑖𝑛 = 0,5 𝑉/𝑟𝑎𝑑/𝑠 ; 𝐾𝑡 = 0,5 𝑁𝑚/𝐴 ; 𝐿 = 2,6 𝑚𝐻 ; 𝐽𝑚 = 5. 10−3 𝑘𝑔. 𝑚2 ; 𝑅 = 2,5 Ω
Réducteur : 𝑘 = 15 ; 𝐽𝑟 = 2,5. 10−5 𝑘𝑔. 𝑚2
Charge : 𝐽𝑐 = 20 𝑘𝑔. 𝑚2
𝐽𝑒𝑞 = 9,4 10−2 𝑘𝑔. 𝑚2 ; 𝑓𝑒𝑞 = 0,01 𝑁𝑚/(𝑟𝑎𝑑/𝑠)
Question 13 : Calculer 𝑇𝑒 et 𝑇𝑚 . A-t-on 𝑇𝑒 ≪ 𝑇𝑚 ?
Pour une fonction de transfert du 1er ordre, 𝑡𝑟5% est le triple de la constante de temps, et donc ici 𝑡𝑟5% = 3𝑇𝑚 .
Question 14 : Calculer 𝐾𝑚 et 𝑡𝑟5% . (page 25 du cours)
Question 15 : A l’aide du théorème de la valeur finale (page 13 et 14) déterminer 𝜔𝑚∞ .
Question 16 : Tracer 𝜔𝑚 (𝑡). Préciser tous les points particuliers et leurs valeurs numériques. (page 25 du cours)
Question 17 : Evaluer la performance de précision du MCC.
S, f, m
V, B Fm (t)
Fr (t)
q(t)
Δp(t)
Grandeurs mécaniques :
S : section utile du piston [𝑚2 ]
x(t) : position du piston [𝑚]
𝐹𝑚 (𝑡) : force motrice de la tige [𝑚]
𝐹𝑟 (𝑡) : force résistante de la charge [𝑚]
𝑓 : coefficient de frottement visqueux [𝑁/(𝑚/𝑠)]
𝑚 : masse de la tige du vérin [kg]
Grandeurs hydrauliques :
q(t) : débit volumique [𝑚3 /s]
V : volume total du fluide dans le vérin [𝑚3 ]
B : coefficient de compressibilité [𝑃𝑎]
𝛥𝑝(𝑡) : différence de pression entre les orifices d'admission et de refoulement [𝑃𝑎]
Equations hydraulique :
Principe de conservation de la masse dans un fluide compressible :
𝑑𝑥 𝑉 𝑑𝛥𝑝
𝑞(𝑡) = 𝑆 (𝑡) + (𝑡)
𝑑𝑡 2𝐵 𝑑𝑡
Relation entre la force et la pression :
𝐹𝑚 (𝑡) = 𝑆 𝛥𝑝(𝑡)
Equation mécanique :
Principe fondamental de la dynamique appliqué à la tige du vérin :
𝑑𝑥 𝑑2 𝑥
𝐹𝑚 (𝑡) − 𝐹𝑟 (𝑡) − 𝑓 (𝑡) = 𝑚 2 (𝑡)
𝑑𝑡 𝑑𝑡
A l'instant initial :
𝑥(0− ) = 0 𝑚
𝛥𝑝(0− ) = 0 𝑃𝑎
Question 1 : Expliquer brièvement le nom de chacune des trois équations et ce qu'elles représentent.
Question 2 : Quelle est la différence entre un vérin simple effet et double effet ?
Question 3 : En précisant l'hypothèse utilisée, appliquer la transformée de Laplace aux 3 équations du vérin.
Meet two robotic bartenders who know how to shake up your night out. They can mix, muddle, and stir it up too.
With moves as fluid as the Pimm’s in your cup, they can create an almost endless combination of cocktails, whether it’s a
classic Manhattan or a custom order of your own design. Just order by app on the
nearby tablets and watch your bionic mixologist do its thing. Designed and powered by
the minds at Makr Shakr, the Bionic Bar® is making history at sea.
https://sciencesindustrielles.com/TPOpale/Ericc/co/Robot_Ericc_1.html
La CdCF annonce les performances suivantes :
Un client du bar parvient à accéder au plan de travail et tape sur l'extrémité du robot.
Question 7 : Comment peut-on modéliser cette sollicitation extérieure ?
Question 8 : Sans faire de calculs, tracer l’allure de la réponse à cette sollicitation.
Exercice 12 : CAMERA
La camera PTZ étanche IP68 ZOOM 28X existe en deux versions : noir et camouflage Otan.
Cette camera est dotée d'un socle aimanté ce qui permet de la positionner sur un véhicule.
Elle est commandée en position angulaire à l’aide de deux moteurs à courant continu.
Le comportement de la caméra en orientation suivant l’axe vertical, est modélisé par la fonction de transfert :
𝛩(𝑝) 9800
=
𝛩𝑐 (𝑝) 10000 + 600𝑝 + 35𝑝2
𝜃𝑐 (𝑡) représente l’angle de consigne en ° par rapport au plan horizontal ;
𝜃(𝑡) représente l’angle atteint en ° par rapport au plan horizontal.
𝒛
𝒛
z = 0,69
𝐾𝐸0 avec 𝐾𝐸0 z>1
D2
0,63𝐾𝐸0
D2
𝜏 𝑇𝑃 𝑇𝑃 𝜏1 𝜏2
2
CI02 TD02 Systèmes asservis 27/11/2023 14/16
𝑧=1
Exercice 14 : CHARIOT AUTONOME D’AMAZON
Amazon pèse 22% des dépenses on line en France, sur un marché de 40,3 M€. L’entreprise doit son essor
économique à sa robotisation et son automatisation extrême ainsi qu’à son service client. Son principale concurrent
français, Rakuten (anciennement C Discount) pèse lui 10%.
Des colis livrés par Amazon sont déplacé par des chariots autonomes.
https://sciencesindustrielles.com/TPOpale/Chariot%20filoguid%c3%a9/co/Chariot_filog
uide_1.html
Pour limiter les effets dynamiques qui pourraient endommager le contenu des
colis, on désire contrôler l’accélération 𝛾(𝑡) d’un point du colis à partir d’une consigne
d’accélération 𝛾𝑐 (𝑡).
Premier réglage
𝛤(𝑝) 𝑎1
Suite à un premier réglage, le comportement du système est modélisé par = .
𝛤𝑐 (𝑝) 1+𝑏1 𝑝
Avec 𝑎1 = 0,95et 𝑏1 = 0,7 𝑠.
Deuxième réglage
𝛤(𝑝) 𝑎2
Suite au deuxième réglage, le comportement du système est modélisé par = .
𝛤𝑐 (𝑝) 1+𝑏2 𝑝+𝑐2 𝑝2
2
Avec 𝑎2 = 0,98; 𝑏2 = 1 𝑠 et 𝑐2 = 0,2 𝑠 .
Troisième réglage
𝛤(𝑝) 𝑎3
Suite au troisième réglage, le comportement du système est modélisé par = .
𝛤𝑐 (𝑝) 1+𝑏3 𝑝+𝑐3 𝑝2
Avec 𝑎3 = 1; 𝑏3 = 0,62 𝑠 et 𝑐3 = 0,2 𝑠 2 .