These
These
These
Spécialité : Mécanique
Par
Wissam KARAM
Ingénieur en génie mécanique
Université libanaise, faculté de génie
- Au professeur Jean-Charles MARE pour m’avoir proposé un doctorat financé sur fonds
propres, pour sa rigueur scientifique, pour sa disponibilité même durant les weekends et
pour ses magnifiques qualités humaines et son hospitalité.
- Au professeur Alain BERLIOZ pour avoir présidé la commission d’examen et géré les
débats.
- Tous les thésards et membres de l’équipe ASCOM qui furent à la fois collègues et
amis : Gérardo, Gustave, Jonathan, Susan, Batoul, Nicolas, Mario …
- Mes parents qui m’ont soutenu lorsque j’ai décidé de faire mes premiers pas vers la
recherche.
- Mes divers colocataires qui ont rendu mes soirées et week-ends agréables.
- Madame Dominique MARE pour son accueil et son hospitalité sans égal.
Remerciements ........................................................................................................................ 3
Introduction............................................................................................................................. 9
Chapitre 1 .............................................................................................................................. 12
Etat de l’art......................................................................................................................... 12
I. Généralités ........................................................................................................................... 12
II. Vérins électromécaniques .................................................................................................. 12
1. Domaine d’utilisation des vérins à vis ............................................................................. 13
2. Différents modes de motorisation .................................................................................... 14
3. Chaîne de transmission moteur-vis .................................................................................. 14
4. Blocage de la tige contre la rotation................................................................................. 16
III. Caractéristiques des vis à glissement................................................................................ 16
IV. Caractéristiques des vis à roulement................................................................................. 18
1. Vis à billes........................................................................................................................ 19
2. Vis à rouleaux .................................................................................................................. 21
3. Précharge.......................................................................................................................... 24
4. Rigidité............................................................................................................................. 24
V. Modélisation des moteurs synchrones ............................................................................... 25
1. Equations électriques ....................................................................................................... 25
2. Couple électromagnétique................................................................................................ 26
3. Equations mécaniques...................................................................................................... 26
VI. Conclusion ........................................................................................................................ 27
Chapitre 2 .............................................................................................................................. 28
Architecture et modélisation .............................................................................................. 28
I. Introduction ......................................................................................................................... 28
II. Architectures d’un EMA .................................................................................................... 28
1. Réducteur : Rotation/Rotation ......................................................................................... 29
2. Transformateur : Rotation/Translation............................................................................. 30
3. Moteurs électriques .......................................................................................................... 31
4. Générations des architectures candidates......................................................................... 32
5. Organigramme récapitulatif ............................................................................................. 33
III. Dimensionnement énergétique.......................................................................................... 35
1. Usinage............................................................................................................................. 36
2. Durée de vie ..................................................................................................................... 36
3. Rendement ....................................................................................................................... 37
4. Couple et inertie équivalente............................................................................................ 39
IV. Modélisation dynamique d’un EMA ................................................................................ 41
1. Bond-graph....................................................................................................................... 41
2. Modélisation de la vis à rouleaux .................................................................................... 42
3. Modélisation du variateur ................................................................................................ 44
V. Dynamiques d’un EMA ..................................................................................................... 47
VI. Conclusion ........................................................................................................................ 49
Chapitre 3 .............................................................................................................................. 50
Banc d’essais...................................................................................................................... 50
I. Généralités ........................................................................................................................... 50
1. Application ou objectif..................................................................................................... 50
2. Hypothèses et choix ......................................................................................................... 50
II. Géométrie du banc d’essai ................................................................................................. 52
1. Structure porteuse ............................................................................................................ 55
2. Raideur d’ancrage ............................................................................................................ 56
3. Arbre et inertie ................................................................................................................. 60
III. Accessoires de l’actionneur hydraulique .......................................................................... 62
1. Accumulateurs ................................................................................................................. 63
2. Capteurs de pression ........................................................................................................ 63
3. Volume............................................................................................................................. 63
4. Fuite réglable.................................................................................................................... 63
5. Sécurité............................................................................................................................. 63
IV. Capteurs et acquisition...................................................................................................... 65
1. Grandeurs à mesurer ........................................................................................................ 66
2. Traitement des signaux .................................................................................................... 66
V. Conclusion ......................................................................................................................... 70
Chapitre 4 .............................................................................................................................. 71
Caractérisation expérimentale............................................................................................ 71
I. Introduction ......................................................................................................................... 71
II. Mesures des différentes raideurs ........................................................................................ 71
1. Structure du banc et ancrage ............................................................................................ 71
2. EMA................................................................................................................................. 73
3. Raideur globale du banc vue par l’EMA ......................................................................... 75
III. Modélisation sous charge nulle de l’EMA........................................................................ 77
1. Approche de modélisation ............................................................................................... 77
Références............................................................................................................................ 161
Abstract................................................................................................................................ 166
Introduction
La consultation et la lecture d’une thèse sont en général plus aisées lorsqu’on connaît les
motivations et le goût des personnes impliquées, en particulier le doctorant et son directeur de
thèse.
L’objet du présent travail de recherche est d’évaluer les facteurs limitatifs (jeu, frottement,
raideur…) des actionneurs électromécaniques utilisés en commande d’effort vis à vis des
performances de puissance, de commande et de durée de vie. Il s’agit ainsi de doter les
fournisseurs de bancs d’essais de connaissances nouvelles dont ils ont besoin pour sélectionner
au mieux leur technologie de chargement.
Le point de départ de ce travail de thèse est caractérisé par l’absence d’expérience sur les
vérins électromécaniques ainsi que l’absence de tout dispositif expérimental au sein du
laboratoire. Comme le montre ce mémoire, le travail a été conduit avec la volonté
d’apporter méthode et connaissance aux intégrateurs de banc d’essais (vision externe de
l’actionneur) plutôt qu’aux concepteurs d’actionneurs.
Tout d’abord une étude bibliographique, nécessaire à tout travail de recherche, porte sur les
transmissions par vis à billes et vis à rouleaux ainsi que su les différents types de moteurs
synchrones. Ce travail permet de caractériser les différentes architectures d’un actionneur
électrique et de définir les différentes combinaisons vis-écrou & moteur.
Une fois l’étude externe accomplie (architecture, géométrie …), nous passons à la modélisation
des différents constituants, essentiellement le système vis-écrou et le moteur. Les modèles ainsi
obtenus sont intégrés dans les environnements de simulation Matlab-Simulink [ 20] et
AMESim [ 19] pour simuler le fonctionnement de l’actionneur.
Une application type, visée tout au long de la thèse, est le test au sol d’un actionneur
aéronautique ayant les points de fonctionnement suivants : 50 kN d’effort bloqué, 30 mm/s de
vitesse à vide et 25 kN à 26 mm/s nominal. Ceci correspond à un actionneur d’aileron de la
famille des Airbus A320 qui doit, pour être qualifié, satisfaire les exigences d’essais suivantes :
- Exigence d’essai # 4 : sur charge à 0.6 mm/s, 9 kN statique + 0.5 kN/mm (3 106
cycles).
La deuxième phase de l’étude est consacrée à la conception et la mise en œuvre d’un banc
d’essai. Un massif châssis métallique supporte face à face un actionneur hydraulique et un
prototype d’actionneur électromécanique, sujet principal de cette thèse. Le vérin
hydraulique représente le vérin aéronautique à tester, il est asservi en position, tandis que
l’actionneur électrique doit être asservi en effort pour reproduire les efforts aérodynamiques
agissant sur l’aileron. Des raideurs d’ancrages ont été prévues coté vérin hydraulique pour
reproduire la déformation de la structure avion. Un ensemble de capteurs (pression, effort,
déplacement …), équipés de conditionneurs et d’afficheurs, transmettent les informations
Afin de déterminer la meilleure méthode pour commander l’actionneur électrique, des essais
préliminaires servent à déterminer le transfert du système. Ceci compense le manque
d’informations entourant le variateur et permet de prédire son fonctionnement.
Ensuite, d’autres essais sont effectués pour déceler d’éventuels problèmes et les dynamiques
non prises en compte dans les modélisations. Ceci permet de perfectionner l’asservissement et
de valider les modèles numériques utilisés par la suite pour la création d’un prototype d’essai
virtuel.
Enfin, une étude analytique et une modélisation dans l’environnement d’AMESim d’un vérin
hydraulique, similaire en terme de puissance, permet de comparer les technologies hydraulique
et électromécanique en asservissement d’effort.
Quant au mémoire, il est divisé en cinq chapitres. Le premier chapitre est consacré à l’état de
l’art de l’actionnement électromécanique pour notre application. Le second chapitre analyse la
conception des actionneurs industriels et propose une méthodologie de dimensionnement et
choix. La conception et le fonctionnement du banc d’essai, mis en œuvre au cours de la thèse,
font l’objet du troisième chapitre. Le quatrième chapitre enchaine avec une caractérisation
expérimentale du vérin électromécanique et des différentes parties du banc ainsi que
l’élaboration d’un modèle virtuel de simulation et de conduite. L’asservissement en effort et
l’étude dynamique sont exposés dans le chapitre cinq qui se termine par la comparaison entre
les technologies électrohydraulique et électromécanique. Enfin, une conclusion générale met
en avant les points clefs relatifs à l’utilisation d’actionneurs électromécaniques comme
générateurs d’effort. Elle suggère aussi des compléments de travaux qui relèvent d’avantage du
domaine de l’automatique.
Chapitre 1
ETAT DE L’ART
I. Généralités
L’actionneur électrique développe des actions mécaniques (effort et/ou déplacement), sur un
corps quelconque, à partir d’une source de puissance électrique. Les mouvements de
translation sont l’objet d’étude de cette présente thèse. Les moteurs linéaires à caractère
industriel sont encore limités à quelques kilo Newtons ([ 15]…) et requièrent d’habitude un
guidage externe de type prismatique. Pour ces raisons, les générateurs d’efforts utilisent un
moteur électrique combiné avec un transformateur rotation-translation, TF, qui assure le
passage d’un mouvement rotatif à un mouvement de translation. Il se trouve que le dispositif le
plus utilisé est le système vis-écrou.
Moteur
~ Variateur
électrique rotation
TF
translation
Réseau
électrique EMA
Fig 1- 1: Principaux composants d’un EMA
La figure ci-dessus montre les trois parties distinctes (électrique et mécanique) qui composent
l’actionneur électrique qui sera nommé dans la suite EMA (electro-mechanical actuator).
D’après les fournisseurs, les nouveaux vérins électriques allient la puissance de l’hydraulique à
la vitesse du pneumatique. La plupart de ces vérins, qui sont construits autour des vis à
rouleaux (il existe aussi des versions avec vis à billes et vis à filet trapézoïdal), repoussent les
limites des performances des vérins linéaires. Ils sont conçus pour des applications à hautes
fréquences d’utilisation demandant une longue durée de vie, une accélération rapide ou une
charge élevée.
Les avantages multiples des vérins à vis leur permettent des utilisations très variées selon les
configurations choisies (type de moteur, de vis et d’accouplement). Ils sont utilisés dans
différents domaines dont on cite : médical, télécommunication, machines agricoles, industrie,
machines outils, aéronautique et aérospatiale, manutention, levage. L’effort maximal
développé varie entre 10 N et 400 KN. Quant à la vitesse maximale, elle est comprise entre 100
et 1500 mm/s. La Fig 1- 2présente les références catalogue (F max , V max ), effort / vitesse
maximaux, des principaux fournisseurs et les classent suivant l’architecture du vérin (les
architectures seront détaillées ultérieurement, page 32). Il est important de noter que l’effort
maximal ne peut pas être développé à la vitesse maximale.
300
MOOG SKF
50
SKF Exlar
Force (KN)
Force (KN)
240 40
180 30
120 20
60 10
0 0
0 250 500 750 1000 1250 1500 0 200 400 600 800 1000 1200
Vmax (mm/s) Vmax (mm/s)
Les moteurs à courant continu (12-36 V DC) sont utilisés pour les vérins de faible puissance
principalement dans le domaine médical (positionnement des lits d’hôpitaux…). Les moteurs
asynchrones fournissent des puissances élevées mais sont peu adaptés au positionnement
précis. Les moteurs « Brushless » synchrones sont utilisés là où un positionnement très précis
est requis en combinaison avec une puissance importante et une grande accélération. Les
nouveaux moteurs synchrones, à rotors creux, permettent quant à eux d’obtenir un
encombrement réduit parce que la vis, l’écrou et la tige du vérin peuvent être insérés à
l’intérieur du rotor.
Le couple produit par le moteur est transmis à la vis via une chaîne de transmission. Elle est
conçue de manière à fournir un rapport de réduction déterminé avec une position bien
particulière du moteur dû à l’encombrement ou juste à un accouplement qui autorise les défauts
d’alignement et amortit les vibrations des deux axes. Plusieurs configurations géométriques
sont proposées comme ci-après.
a) Axes perpendiculaires
Dans cette configuration l’axe du moteur est perpendiculaire à l’axe de la vis. La connexion
entre les deux se fait au moyen d’engrenages coniques ou d’un système à vis sans fin (vis à
l’extrémité de l’axe du moteur et roue à l’extrémité d’axe de la vis).
b) Axes parallèles
Dans cette configuration l’axe du moteur est parallèle à l’axe de la vis. Dans la plupart des cas,
la transmission est confiée à une courroie mais on peut aussi trouver des transmissions par train
d’engrenages.
c) Axes en lignes
Dans cette configuration l’axe du moteur et l’axe de la vis sont en ligne. Moteur et vis sont
alors reliés soit par un accouplement élastique, soit par un train d’engrenages planétaires où
l’un des axes est la couronne et l’autre le soleil (on obtient alors un rapport de réduction élevé
dans un encombrement réduit).
L’axe de la vis (ou l’écrou si la vis se translate) peut aussi servir comme un axe pour le moteur
(pas de réduction) ce qui permet de diminuer la longueur du système et d’obtenir une rigidité
de liaison maximale. On parle alors de vérins compacts et d’entrainement direct.
Par définition, un vérin est équivalent à une liaison prismatique motorisée. Lors de la rotation
de la vis dans l’écrou, les efforts de frottement, bien que minimes dans le cas des vis à
roulements, créent sur l’écrou un couple selon son axe. Ceci est une caractéristique spécifique
des vérins électriques pour lesquels la fonction d’antirotation est souvent rejetée sur le bâti et
sur la charge. La liaison mécanique entre le vérin et l’extérieur ne peut plus être réalisée par
des liaisons rotules qui doivent être remplacées par des liaisons pivot qui rendent le système
hyperstatique. Par contre, lorsque la charge ne peut supporter que les efforts axiaux, il faut
incorporer dans le vérin un système qui bloque la rotation de la tige de vérin et reprend le
couple créé par frottement et effet inertiel (Annexe 1).
Tige du vérin
Guidage
rectiligne
Galets
Vis
Ecrou
Structure du vérin
La Fig 1- 9 montre l’une des méthodes utilisées pour bloquer l’écrou contre la rotation. Elle
consiste à utiliser 2 galets dans des guidages rectilignes. Ils doivent reprendre le couple produit
par le frottement. Suivant l’importance du projet, ces galets peuvent être équipés de roulement
à aiguilles (au lieu de bagues en bronze ou en polyamide) pour diminuer autant que possible les
pertes par frottement et les jeux.
Dans la suite de l’étude, plusieurs développements feront référence à la vis à glissement. Il est
donc important de s’y attarder un peu [ 21]. La Fig 1- 10 représente une vis à glissement à filet
carré, ayant un pas p [m], un diamètre moyen d [m] et un angle d’hélice φ [rd].
écrou
p
d
Fig 1- 10 : Vis à glissement à filet carré
La rotation d’un élément (vis ou écrou) entraîne la translation de l’autre d’une quantité x
proportionnelle à l’angle de rotation θ [rd] et à l’avance par tour l [m] (lead). Il est intéressant
de noter que l’avance par tour n’est autre que le pas du filet (distance crête à crête) multiplié
par le nombre d’entrées de la vis n v :
θ
x= l
2π (1- 1)
l = p nv (1- 2)
l
tan ϕ = (1- 3)
πd
La figure ci-dessous présente les réactions exercées par l’écrou sur la vis au niveau des
surfaces de glissement lorsqu’une charge axiale F [N] est appliquée à la vis.
Résistant Aidant
Filet de Filet de
F la vis
F la vis
Fres µN
y µN Faid
φ N l φ N l
Ecrou Ecrou
x
Cr Ca
πdm + πdm
Fd ⎛ l + πμd ⎞
Charge résistante ⇒ Cres = ⎜ πd − μl ⎟
2 ⎝ ⎠ (1- 4)
Fd ⎛ −l + πμd ⎞
Charge aidante ⇒ Caid = −
2 ⎜⎝ πd + μl ⎟⎠ (1- 5)
On remarque que le couple aidant C aid peut changer de signe et être soit négatif soit positif.
Ceci dépend en effet de l’angle d’hélice et du coefficient de frottement µ. Pour 0 <µ< tg φ,
l’effort de frottement n’est pas suffisamment important pour retenir la charge si elle est aidante,
elle se déplace alors d’elle même (d’où la nomination de charge aidante) et le couple externe
C aid est alors un couple de freinage. Par contre lorsque µ> tg φ, on dit que la vis « s’auto
verrouille », le système vis-écrou est alors irréversible.
Pour information, le Tableau 1 donne les valeurs pratiques du coefficient de frottement entre
les divers matériaux composant les vis à glissement.
Les vis à roulement (billes ou rouleaux satellites) remplacent le frottement de glissement par
un frottement de roulement. Ceci entraîne une diminution considérable de la perte d’énergie à
travers la transmission et réduit énormément l’échauffement et l’usure. Ces vis conservent
donc la même précision quelle que soit la vitesse de fonctionnement. Du point de vue
dissipation, la vis à roulement est donc équivalente à une vis de glissement ayant un coefficient
de frottement µ plus faible.
Les principaux avantages des vis à roulement par rapport aux vis à filet trapézoïdal sont :
1. Vis à billes
Une vis à billes comprend essentiellement une vis, un écrou fileté intérieurement et un
ensemble de billes circulant dans le chemin formé par les filetages de l’écrou et de la vis.
Après avoir effectué un certain nombre de rotation autour de la vis, les billes sont ramenées à
leur position initiale dans l’écrou par un canal de retour. Ce dernier peut présenter différents
aspects suivant les constructeurs (tuyau externe, interne ou multilignes). Suivant le pas de la
vis, on peut trouver des vis à filets multiples pour augmenter le nombre de points de contacts et
par la suite la charge admissible pour une même taille d’écrou (sans changer le diamètre des
billes). Les billes qui circulent continuellement diminuent le frottement entre la vis et l’écrou et
transmettent la charge. Par suite du faible frottement de roulement, le rendement est élevé,
typiquement 90 % et les vis à billes conservent leur précision beaucoup plus longtemps que les
vis conventionnelles.
Les billes en circulant, diminuent le frottement entre la vis et l’écrou mais créent un frottement
entre elles-mêmes et une quantité de bruit perceptible surtout dans le canal de retour. Certains
des constructeurs ont modifié la structure conventionnelle pour réduire ces défauts. Ci-dessous
sont présentées deux des méthodes utilisées pour améliorer le rendement et diminuer l’usure
des billes.
La première consiste à intercaler entre les billes des cages en polyamide qui remplacent le
frottement acier-acier (bille contre bille) par un frottement acier-polyamide (bille contre cage)
beaucoup plus faible. Ces cages permettent aussi d’éliminer le bruit de collision entre les billes
(intensité du bruit réduite de 5 dB) et d’y intercaler un film d’huile ou de graisse. Cette
structure n’altère pas la charge maximale supportée car le nombre de billes travaillant reste
invariant.
Fig 1- 14 : Structure conventionnelle [ 11] Fig 1- 15 : Structure avec cage à billes [ 11]
La deuxième solution consiste à intercaler entre les billes une bille morte. Son diamètre étant
inférieur à celui des autres billes, elle roule sans toucher la vis ni l’écrou. Il n’y a alors plus de
frottement entre billes portantes. L’inconvénient de cette méthode est que le nombre de billes
portantes est réduit de moitié. Il faut alors un écrou plus long ou un diamètre de bille plus
grand. Les billes portantes peuvent être réalisées en céramique ce qui diminue le coefficient de
frottement entre les billes et les rails et autorise une augmentation de la pression de contact.
Fig 1- 16 : Structure conventionnelle [ 41] Fig 1- 17 : Structure avec billes mortes [ 41]
2. Vis à rouleaux
Dans ce genre de vis, le corps roulant intermédiaire est constitué par un rouleau fileté. La
cinématique de ce système est beaucoup plus complexe. L’ensemble filetage de l’écrou,
filetage des rouleaux satellites et filetage de la vis forment un train d’engrenages planétaires
hélicoïdaux à grand angle d’hélice.
Vis à Vis à
rouleaux billes
Points de contact
supplémentaires
Le nombre de points de contact entre les composants est supérieur à celui présent dans les
systèmes vis à billes, ce qui permet un chargement plus important pour un même diamètre et
une même taille d’écrou.
Du fait de leur aptitude à déplacer des fortes charges pendant des milliers d’heures dans des
conditions très difficiles, les vis à rouleaux satellites sont particulièrement bien adaptées aux
applications les plus exigeantes. La robustesse de l’écrou leur permet de résister aux chocs et le
système de synchronisation des rouleaux assure un fonctionnement fiable même dans les
environnements difficiles. Enfin, leurs grandes avances par tour et la symétrie de l’écrou
permettent d’obtenir des vitesses linéaires élevées.
La vis filetée (1) présente un filetage à entrées multiples, de pas relativement grand. L’angle
sur les flancs est de 90° et le profil est triangulaire. L’écrou (2) possède un filetage intérieur
identique à celui de la vis. Les rouleaux (3) possèdent un filetage à une entrée, dont l’angle
d’hélice correspond à celui de l’écrou. Il ne se produit ainsi aucun déplacement axial entre
l’écrou et les rouleaux. Une recirculation des rouleaux n’est donc pas nécessaire. Les flancs des
filets des rouleaux sont bombés. Les filets des rouleaux ont un pas égal au pas apparent de la
vis et de l’écrou (pas de la vis / nombre des entrées).
Les rouleaux présentent à chacune des extrémités un pivot cylindrique et une denture. Les
pivots sont montés dans les alésages des porte-rouleaux (5). Ainsi, les rouleaux sont maintenus
à des distances régulières. Les porte-rouleaux sont disposés flottants dans les écrous et sont
maintenus axialement par des joncs (6). Les dentures des rouleaux s’engrènent dans celles des
couronnes (4) fixées dans l’écrou. Les rouleaux sont ainsi guidés parallèlement à l’axe et un
fonctionnement parfait est assuré.
Ce type de vis filetée (7) présente des pas très fins et est utilisé lorsque l’on a besoin d’une très
grande précision de positionnement associée à une grande rigidité et à une capacité de charge
élevée. Les principaux éléments sont les mêmes, la vis (7), l’écrou (8) et les rouleaux (9) qui
sont guidés et maintenus à distance dans une cage (10). La vis possède un filet à une ou deux
entrées seulement (sinon le filet serait très fin et ne supporterait pas la charge) avec un profil
triangulaire. L’angle sur les flancs est de 90°.
L’écrou possède un filetage intérieur identique à celui de la vis. Les rouleaux ne possèdent pas
de filet, mais des gorges disposées perpendiculairement à l’axe de la vis (on peut dire qu’ils ont
des filets avec un pas nul). La distance entre les gorges correspond au pas apparent de la vis et
de l’écrou. Lors d’une rotation de la vis ou de l’écrou les rouleaux se déplacent axialement
dans l’écrou. Après un tour complet chaque rouleau est ramené dans la position initiale par
deux cames (11) fixées aux extrémités de l’écrou. Cette recirculation des rouleaux est rendue
possible par une rainure longitudinale dans l’écrou.
3. Précharge
Les écrous présentent la possibilité d’être préchargés, ce qui permet d’éliminer le jeu vis-écrou
et d’augmenter la rigidité tout en améliorant l’exactitude de positionnement. La figure suivante
présente la déformation axiale, au niveau de la liaison entre la vis et l’écrou, en fonction de la
charge appliquée pour les écrous simples et les écrous préchargés.
Ecrou préchargé
Charge axiale
La précharge est appliquée en plaçant deux écrous l’un contre l’autre et en les tournant l’un par
rapport à l’autre en direction opposée. Pour éviter une réduction significative de la durée de
vie, la précharge ne doit pas dépasser 1/3 de la charge axiale réelle maximale. En fonction de
l’application, les systèmes d’écrou sont préchargés de 7 % à 10 % de la capacité de charge
dynamique. Dans le cas des vis à rouleaux inversées (filetage sur l’écrou), la précharge n’est
pas mise en œuvre.
4. Rigidité
La rigidité du système vis-écrou influe la dynamique de l’actionneur. Elle aussi est affectée par
les autres pièces telles que paliers, supports, boîtiers d’écrou, etc. La rigidité axiale équivalente
k eq [N/m] est alors la combinaison de plusieurs rigidités en série :
1 1 1 1 1
= + + +
keq k pa kec kv kve (1- 6)
Avec :
a) Rigidité vis-écrou
La rigidité de la liaison entre la vis et l’écrou k ve est difficile à calculer. Cela est dû au fait
qu’on a une raideur de contact entre la vis et les rouleaux (ou billes) k vr [N/m] et une raideur de
contact entre ces derniers et l’écrou k er [N/m].
F Rouleau
Vis
Ecrou
Rer
F
Rvr
1 1 1
= +
kve ker kvr (1- 7)
La raideur équivalente k ve est d’habitude fournie par les constructeurs sous forme de formule
empirique ou d’un tableau en fonction du diamètre de la vis et du type de l’écrou (taille,
précharge, matériau…). Elle est considérée comme constante pour une vis donnée.
Le moteur asynchrone associé à des convertisseurs statiques simples permet d’obtenir des
performances très voisines de celles du moteur à courant continu [ 23]. Pour faire varier la
vitesse du moteur synchrone [ 31], il faut varier la fréquence de ses tensions et courants
d’entrée. Le plus souvent on utilise un variateur composé de deux convertisseurs en cascade :
un redresseur alimenté par le réseau fournit une tension continue, qui par la suite est convertie
en tension alternative de fréquence variable à l’aide d’un onduleur.
1. Equations électriques
Les équations électriques d’un moteur synchrone sont données par le triplet ci-dessous qui
présente les relations relatives à chacune des phases. L’inductance L s tient compte de
l’inductance propre d’une phase et de l’inductance mutuelle entre phases comme proposé dans
la référence [ 30].
dI1
Ls = V1 − Rs I1 + K ω sin( θe )
dt
dI
Ls 2 = V2 − Rs I 2 + K ω sin( θe − 2 π / 3)
dt
dI
Ls 3 = V3 − Rs I 3 + K ω sin( θe − 4 π / 3)
dt (1- 8)
Où :
2. Couple électromagnétique
Le couple électromagnétique C e [Nm] exercé sur le rotor par les bobinages du stator est la
somme des couples produits par chacune des phases. Son expression est alors donnée comme
suit :
3. Equations mécaniques
dω
Jr = Ce − C r
dt (1- 10)
avec J r [kgm2] l’inertie du rotor et C r [Nm] l’ensemble des couples résistants (frottement,
charge …).
VI. Conclusion
Etant donné le très faible nombre d’articles traitant de la commande en effort des EMA
linéaires, à l’exception de ceux qui présentent la commande des actionneurs des minirobots
dont la gamme de puissance est totalement différente de la notre, nous avons consacré ce
premier chapitre à l’état de l’art des EMA.
Tout d’abord, nous avons montré des différents avantages mettant en valeur les actionneurs
électriques. Ensuite, nous avons exposé les diverses fonctions et composants rentrant dans la
conception des EMA. Les différents types de motorisation, de transmission et de
transformation du mouvement (ou de l’effort) ont alors été présentés avec leurs
caractéristiques, principes de fonctionnement et modèles usuels.
Ce chapitre est important à nos yeux car il permet d’apporter une vision globale sur les EMA
linéaires, avant le développement de considérations plus approfondies.
Chapitre 2
ARCHITECTURE ET MODELISATION
I. Introduction
Ce chapitre est dédié à la présentation des actionneurs disponibles chez les différents
fournisseurs et à la méthodologie de choix. Trois principaux aspects sont abordés : la sélection
de l’architecture d’EMA la plus adaptée, le dimensionnement en puissance et enfin la
modélisation dynamique préliminaire pour le prototypage numérique.
Les actionneurs d’efforts dynamiques doivent fonctionner dans des conditions très sévères.
Tout d’abord, ils doivent développer des efforts alternatifs pouvant changer de signe. Toute
présence de jeu peut engendrer des chocs qui à leur tour vont décroître la durée de vie et la
performance du système. Ensuite, tout déphasage dans la commande entre l’effort désiré et la
consigne réalisé change les caractéristiques de l’effort produit, produisant par exemple un effet
inertiel parasite au lieu d’un effet élastique pur. En conséquence, l’inertie équivalente des
pièces en mouvement, les efforts de frottement, les raideurs et les jeux doivent être
soigneusement pris en compte pendant la phase de conception afin d’assurer une cohérence
entre les choix de conception et les performances visées en boucle fermée. Quand ceci est bien
étudié, l’actionneur devient facile à contrôler.
Le chapitre précédent a exposé d’une manière générale quelques composants des actionneurs
électriques. La présente partie a comme but de cibler les différentes architectures possibles,
surtout celles disponibles en standard sur le marché et de relever les avantages et les
inconvénients de chacune d’elles. Pour cela, une étude architecturale des principaux éléments
d’un EMA est conduite afin d’analyser les choix de conception.
1. Réducteur : Rotation/Rotation
Les architectures et les géométries des EMA ont une grande influence sur la performance
finale. Dans les applications peu exigeantes, les réducteurs sont utilisés pour adapter les
caractéristiques nominales du moteur (couple nominal C n et vitesse angulaire nominale ω n ) sur
la charge (force nominale F n et vitesse linéaire nominale V n ). Alors une application nécessitant
un grand effort et une petite vitesse peut être motorisée par un moteur ayant un faible couple
mais une grande vitesse angulaire. Les principaux types de transmission mis en jeu
sont présentés avec leurs avantages et leurs inconvénients dans le tableau ci dessous:
Avantages Inconvénients
- Conception en tandem - Présence de jeu
Engrenages droits - Axes du moteur et de la vis parallèles - Largeur du système est
- Longueur du système réduite plus importante
- Axes du moteur et de la vis - Présence de jeu
Engrenages coniques
perpendiculaire - Taille importante
- Axes du moteur et de la vis
- Grande inertie équivalente
Roue et vis sans fin perpendiculaire
- Irréversibilité
- Grand rapport de réduction
- Grand rapport de réduction
- Grande inertie équivalente
Engrenages épicycloïdaux - Réducteur très compact
- Présence de jeu
- Axes en ligne
- Pas de jeu
- Transmission raide
- Besoin d’un grand couple
Liaison directe - Petite inertie
moteur
- Système très compact
- Axes identiques
Il sera montré ultérieurement que les réducteurs induisent d’énormes inerties et masses
équivalentes. Pour ces raisons, un entraînement direct est souhaitable car il permet d’obtenir
une transmission réversible, rigide et sans jeu. Il reste à noter que cette solution, bien
qu’avantageuse au niveau des performances, va entraîner l’installation d’une puissance motrice
bien supérieure (jusqu’à 10 fois) à la puissance requise. Ceci vient du fait que pour les
solutions industrielles, les mêmes moteurs, ayant des vitesses de rotation standard (2800 ou
4200 tr/min), sont employés pour toutes les applications, peu importe la grandeur de la vitesse
désirée. Par contre, pour les applications aéronautiques, la conception sur mesure des moteurs
peut être envisagée afin limiter autant que possible le poids et la puissance inutile.
2. Transformateur : Rotation/Translation
Avec leurs rendements médiocres, les vis à friction sont en général irréversibles et présentent
des tailles excessives lorsqu’une longue durée de vie est combinée à un effort dynamique
élevé. Remplaçant alors le glissement par du roulement, les vis à billes ou à rouleaux, devenues
depuis quelques temps des composants mécaniques standard chez différents fournisseurs
d’EMA (SKF-Transroll [ 7], Thomson-Danaher [ 1], Rollvis [ 8], Nook [ 13], THK [ 11],
Spiracon [ 6]…), peuvent offrir de bien meilleures prestations.
Ecrou en
rotation
a b
Ecrou Vis
c d
Vis en
rotation
Vis Ecrou Tube
Les conceptions « a » et « c » ne sont pas applicables pour les vis à billes (recirculation des
billes impossible au niveau de la vis). Quant aux vis à rouleaux, leur course est limitée à cause
des contraintes d’usinage de l’écrou. Par contre, ces conceptions permettent un bon guidage de
la tige du vérin vu que la partie externe est circulaire et que les parties filetées sont internes et
bien protégées de l’environnement externe.
La conception « b » permet d’importantes courses mais requiert une enveloppe pour protéger
la vis contre la pollution et contenir le lubrifiant. En plus, les guidages linéaires sont difficiles à
intégrer et nécessitent un alésage tout au long de la vis pour y placer un guidage interne.
Dans la conception « d », l’écrou est lié à la charge via un tube rigide. Ce tube présente la
solution à la plupart des problèmes rencontrés avec la conception « b » (pollution, lubrification
et guidage). La facilité et la simplicité du guidage de la tige dans cette conception permettent
de grandes courses (limitées que par le flambage et la vibration de la tige). Lorsqu’une grande
course est envisagée, pour mieux tenir et guider la vis, un roulement à bille est placé à
l’extrémité de celle-ci afin de former une liaison linéaire annulaire entre la vis et la tige du
vérin et assurer ainsi l’isostatisme de la vis (Fig 2- 5, g).
Il faut noter que dans la plupart des conceptions, la charge assure la fonction d’antirotation.
Ceci a deux inconvénients importants qui ne concernent pas les actionneurs hydrauliques. Tout
d’abord, les embouts à rotule ne peuvent pas être utilisés pour connecter l’EMA à la charge et
au banc d’essai. Ensuite, la charge doit être conçue de telle sorte à pouvoir supporter le couple
de friction ainsi que le couple de démarrage entre la vis et l’écrou transmis par la tige de
l’actionneur. On peut citer comme exemple : vérins électriques des inverseurs de poussée de
l’A380 (Fig 2- 2).
3. Moteurs électriques
En écartant les moteurs à balais qui ne peuvent pas supporter les contraintes d’une durée de vie
élevée, il reste les moteurs synchrones qui offrent la solution la plus efficace pour un
fonctionnement correct et identique dans les quatre quadrants de puissance comparés à leurs
confrères asynchrones. Ils ont une taille réduite et fournissent des couples élevés même à
vitesse nulle. Ils peuvent être fournis avec des rotors, à technologie innovatrice, creux
permettant des conceptions compactes ou des rotors standards pleins. La plupart des moteurs
sont équipés de frein de sécurité à manque de courant et intègrent le résolveur qui réalise la
mesure de position nécessaire à la commutation électronique.
Roulement Connecteur
Frein
Stator
Résolveur
Rotor
Roulement Stator
Rotor à
aimants
permanents
Fig 2- 3 : Motor à rotor creux [ 15] Fig 2- 4 : Motor à rotor plein [ 14]
Des exemples génériques d’actionneurs [ 37] à entraînement direct conçus par différents
fournisseurs sont montrés dans la Fig 2- 5.
c d
Vis filetée & écrou en rotation
e f
Ecrou fileté et en rotation
Comme première approche, on remarque que ces différentes conceptions sont le fruit de
l’union de l’un des types de systèmes vis-écrou avec l’une des technologies de moteurs
synchrones. Quant aux efforts axiaux, ils sont repris par une paire juxtaposée de roulements à
billes à contact oblique. Un petit roulement supplémentaire assure l’isostatisme, de la partie en
rotation (rotor et écrou), en bloquant les déplacements transversaux.
Les conceptions en ligne f et g, bien que leur système vis-écrou soit de type différent,
augmentent la taille de l’actionneur. Ceci résulte d’un positionnement en ligne du moteur et du
système vis-écrou. Les conceptions concentriques c, d et e profitent de l’axe creux du rotor
pour rendre l’actionneur plus compact. Cette fois-ci, la tige du vérin, lorsqu’elle se rétracte
vient se placer à l’intérieur du rotor, ci qui diminue la longueur totale de l’actionneur.
Quant aux systèmes vis-écrous utilisés, les architectures c et d y sont basées sur la
conception « b ». Par ailleurs ces deux architectures ne présentent pas de solution au problème
de pollution mais chacune a résolu le problème de guidage d’une manière différente.
L’architecture c utilise un guidage interne à la vis tandis que l’autre possède un épaulement, à
son bout non articulé, et joue le rôle d’un guidage externe. Ces deux architectures profitent du
vide à l’intérieur de leur vis pour y loger un capteur de déplacement linéaire à technologie
LVDT (potentiomètre inductif linéaire variable).
Le reste des architectures (e,f et g) est du point de vue technologique plus simple et plus
facile à réaliser. Les tiges des vérins étant lisses peuvent se contenter d’un guidage primitif en
bronze tandis que leur mécanisme est bien protégé à l’intérieur des carters. Il reste à noter que
l’architecture g, bien qu’ayant la taille la plus grande, est la seule à pouvoir offrir des courses
importantes (plusieurs mètres).
5. Organigramme récapitulatif
Transmission
Vis
Ecrou
Moteur
Tige
Filetage carré
Rouleaux
Planétaire
Reprendre
l’effort axial : Par exemple : n paires de
roulement à billes à
Roulements contact oblique
Dimensionner « en puissance » un EMA consiste à choisir, parmi des éléments standards, les
composants à intégrer, principalement le moteur et la vis (le réducteur est éliminé à cause de
son influence sur l’inertie ramenée). Le domaine opérationnel du moteur est d’habitude défini
par les données du fournisseur (C max , ω max ), C max (Nm) étant le couple qui correspond à
l’intensité maximale à vitesse nulle et ω max (rd/s) est la vitesse du moteur lorsque ce dernier est
alimenté avec la tension maximale à couple nul. En conséquence, la capacité du moteur en
termes de puissance se trouve limitée par un domaine quasi rectangulaire, sauf dans les cas où
il y a une limitation de puissance (les isopuissances sont des hyperboles dans le plan
puissance). Il est intéressant de noter que, hormis la différence de rendements, l’EMA a un
fonctionnement identique dans les quatre quadrants de puissance.
C
Cint
P = Pint
Service intermittent
Cmax
Service permanent P = Pmax
ωmax ω
Au premier ordre, la transformation d’une puissance électrique en une puissance mécanique est
fixée par trois paramètres : ceux du moteur (C max , ω max ) et l’avance par tour de la vis l. Notons
que, l’actionneur peut fonctionner dans des domaines intermittents et fournir des efforts plus
importants. Ces domaines n’étant pas clairement définis par les fournisseurs du moteur et de
l’électronique de puissance, ne sont pas pris en compte dans la phase de dimensionnement en
régime établi. Au besoin, un dimensionnement thermique équivalent ([ 32] et [ 33]) permet de
mettre à profit le domaine de fonctionnement intermittent.
1. Usinage
La paire (l, d) doit être choisie parmi un ensemble de valeurs standards relatives aux unités SI
ou anglo-saxonnes pour écarter toute possibilité d’interférence entre l’outil d’usinage et la vis
(respectivement l’écrou).
2. Durée de vie
L’actionneur doit être capable d’accomplir le nombre de cycles demandés avec une fiabilité
donnée. La durée de vie dépend directement de la capacité dynamique de la vis à rouleaux qui
à son tour dépend de la paire (l, d). Les vis ayant des diamètres plus grands possèdent des
surfaces de contact et de travail plus importantes ce qui se traduit par un accroissement de la
capacité de charge dynamique.
La durée de vie, d’un écrou simple, en mètre parcourus, correspondant à une fiabilité de 90%,
notée L 10m [m], est calculée grâce à la formule ci-dessous. Celle ci est commune aux systèmes
à éléments roulants, elle donne la durée de vie en fonction de la capacité dynamique théorique
C dyn [N] ainsi que de l’effort moyen F m [N] :
3
⎛C ⎞
L10 m = l ⎜ dyn ⎟ 106
⎝ Fm ⎠ (2- 1)
Il est important de noter que dans cette formule, l’utilisation de l’avance par tour l permet
d’obtenir une durée de vie exprimée en distance parcourue et non pas en nombre de tours.
Lorsqu’une fiabilité différente est requise, des facteurs correcteurs sont ajoutés à la fonction ci-
dessus.
Pour les cycles complexes, l’effort moyen équivalent est calculé en utilisant les efforts
minimaux et maximaux, F mi et F Mi [N], ainsi que l’amplitude A i [m] et le nombre N i , de
chacun des cycles :
∑N A
3
i i
i (2- 2)
Quant à la distance correspondante parcourue pendant ces cycles, elle est donnée par :
dist = ∑ 4 N i Ai
i (2- 3)
Contrainte 2 : Le dimensionnement du système vis-écrou doit alors satisfaire les
On déduit de ce qui précède que deux vis à rouleaux (ou a billes) ayant des capacités
dynamiques identiques peuvent effectuer les cycles prévus avec des fiabilités différentes. En
effet, ceci dépend de l’avance par tour l de chacune des vis.
Quant aux écrous préchargés, ils permettent une transmission sans jeu. Cependant la précharge
diminue la durée de vie qui doit alors être calculée pour chaque partie de l’écrou [ 7].
3. Rendement
Les vis à rouleaux à haute performance réduisent les coûts initiaux (taille des moteurs) ainsi
que les coûts de fonctionnement (consommation d’énergie). Outre ces raisons économiques, la
discontinuité des rendements autour de la vitesse nulle ainsi que l’irréversibilité de la
transmission nuisent à la conception d’une commande performante. Une attention spéciale doit
être alors attribuée à l’évaluation des rendements d’un système vis-écrou.
- Rendement direct (η d ) pour les charges résistantes, quand la puissance est transmise du
moteur à la charge,
- Rendement inverse (η i ) pour les charges aidantes avec une transmission réversible,
quand la puissance est transmise de la charge au moteur qui joue le rôle de frein,
- Pseudo rendement (η p ) pour les charges aidantes avec une transmission irréversible,
quand le moteur et la charge unissent leurs puissances pour vaincre les frottements (on
montre facilement que la réversibilité est assurée lorsque ηd > 0.5 [ 21]).
Le calcul détaillé du rendement des vis à rouleaux étant complexe, on les assimile en pratique à
des vis de glissement ayant les mêmes diamètres d et avance par tour l mais présentant un
coefficient de frottement µ modifié par la présence du roulement. Ainsi, leurs rendements sont
exprimés comme suit [ 37] :
1 1 1
ηd = ηp = =
1+ μ /β ηi = 1 − μ / β 1 − μ / β ηi (2- 4)
où β = l / πd est la tangente de l’angle d’hélice φ.
Le coefficient de frottement, dans les conditions dynamiques, est d’habitude indiqué par les
fournisseurs comme étant égal à 0.01 pour les vis à rouleaux (0.003 à 0.01 pour les vis à billes
contre 0.04 à 0.11 pour les vis à glissement). Il n’y a pas d’informations précises concernant les
coefficients de frottement de démarrage qui dans certains cas peuvent atteindre le double de la
valeur dynamique. A partir des équations ci-dessus, nous avons tracé les rendements direct et
indirect ainsi que le pseudo rendement d’une vis à rouleaux en fonctions de l’angle d’hélice φ
pour µ = 0.01.
1.2
Pseudo rendement
1
Rendement direct
0.8
Rendement
Rendement inverse
0.6
0.4
ηi >0. 9 ηd
0.2
Reversible
0
0 0.5 1 1.5 2 2.5 3 3.5 4
Nous proposons de représenter sur la Fig 2- 10 les trois contraintes précédentes dans le plan
(l, d) sur lequel nous superposons les combinaisons (l, d) disponibles dans les catalogues et
symbolisées par des marques triangulaires. Pour choisir la paire la plus convenable parmi les
valeurs satisfaisantes, des critères additionnels doivent être pris en considération, spécialement
ceux qui introduisent l’adaptation moteur/vis. C’est l’objet du paragraphe suivant.
»
d
η
9
0.
Domaine satisfaisant
>
η
i
les 3 contraintes
«
p
su
ite
m
Li
Limite inf
« duré
durée de vie »
inf
Limite »
a e
g
« usin
l (mm)
Fig 2- 10 : Contraintes agissant sur la sélection de (l, d)
Pour les charges résistantes, le couple moteur nécessaire C n [Nm] est proportionnel à l’avance
par tour l :
Fl
Cn = = g (l )
2 π ηd (l ) (2- 5)
De ce point de vue, choisir le plus petit l possible diminue le couple moteur nécessaire à la
réalisation de l’actionneur. D’un autre coté, la masse équivalente en translation m e [kg] de
l’actionneur doit être minimisée autant que possible. En effet, d’importantes valeurs de m e
peuvent causer des chocs nocifs lorsque l’actionneur rencontre un obstacle comme les butées
mécaniques ou consommer du couple moteur lors des phases d’accélération.
Par exemple, une référence d’actionneur pouvant fournir 50 kN d’effort bloqué avec 3 mm
d’avance par tour produit une masse équivalente de 37000 kg ! En ce qui concerne la
commande, une grande inertie accroît le temps de réponse et diminue la fréquence naturelle de
l’actionneur. Finalement, la « masse équivalente » en translation (nommée aussi « masse
réfléchie en sortie ») de l’actionneur est donnée par :
2
⎛ 2π ⎞ 1
meq = J r (l ) ⎜ ⎟ + mt = f ( 2 )
⎝ l ⎠ l (2- 6)
Avec :
Comme le montre Fig 2- 11, à chaque diamètre de vis correspond une avance par tour l
optimale. La première paire (l, d) ne faisant pas partie du domaine satisfaisant (Fig 2- 10) est
éliminée au profit de la seconde paire (30, 3) qui optimise globalement le dimensionnement du
système vis-écrou.
2,2E+06
d =36 mm
2,0E+06
1,8E+06
C0 me (Nm.kg)
d =30 mm
1,6E+06
Avance par
1,4E+06 tour optimale
d =25 mm
1,2E+06
1,0E+06
8,0E+05
0,0 1,0 2,0 3,0 4,0 5,0 6,0 7,0
l (mm)
Pour conclure, le choix des composants d’un EMA commence par la détermination des
composants acceptables comme il vient d’être vu et se termine par une sélection basée sur un
compromis entre taille-prix et performance dynamique. Les aspects raideurs ne sont pas
considérés dans cette approche préliminaire.
La modélisation de l’EMA, comme toute modélisation, passe par plusieurs étapes. L’approche
choisie ici passe tout d’abord par l’établissement d’un Bond-graph qui sert à situer chaque
composant par rapport au système global et définit les interfaces ainsi que les échanges de
puissance entre les différents composants. Ensuite, la modélisation de chacun des composants
s’effectue d’une manière distincte en respectant les exigences préalablement définies par le
Bond-graph. A la fin, le modèle obtenu est implémenté selon les besoins dans les
environnements de simulation, comme Matlab-Simulink ou AMESim. L’utilisation simultanée
des deux logiciels relève en fait de la capacité du premier à fournir des boites à outils pour la
CAO de l’Automatique dans des approches « signal », tandis que le second permet la
modélisation réaliste et robuste des composants ou phénomènes selon une approche
« physique » multiports.
1. Bond-graph
La structure du modèle de l’EMA est présentée par le Bond-graph de la Fig 2- 13. Elle est
organisée de manière à reproduire la combinaison fonctionnelle entre l’électronique de
puissance, le moteur électrique et la transmission par vis-écrou. Les liaisons de puissance
fonctionnelles sont mises en valeurs (en gras). Conformément aux architectures e et f de la
Fig 2- 5, le moteur entraîne l’écrou tandis que la tige se translate. Les barres de causalités sont
distribuées sans aucun conflit et fixent les interfaces des sous modèles. Le capteur d’effort est
modélisé comme un ressort combiné à un amortissement structural. De même, d’autres
raideurs (roulements, ancrage sur la structure ou la charge, tige …) peuvent être considérées si
elles influencent significativement le comportement dynamique. Par ailleurs, le rendement
mécanique est traduit par une perte de couple C f , influencée par la vitesse de la charge.
Variateur Moteur
Vis-
Vis-écrou Jv
(L, R, KT) (Jr, fm)
kant
(r, Cfrot, kv)
me mt Charge
kr kF
Une modélisation réaliste du système vis-écrou est d’une importance particulière pour obtenir
des simulations précises de l’EMA. On peut remarquer sur le Bond-Graph (Fig 2- 13) qu’un
modèle à quatre entrées-sorties a été développé englobant à la fois la rotation et la translation
de la vis et l’écrou. Ce choix permet la modélisation de la raideur axiale ainsi que de la raideur
du système antirotation qui pourront significativement affecter la réponse dynamique de
l’actionneur (traitée ultérieurement). Cependant, et comme il a déjà été dit, la plupart des EMA
dotés d’une haute performance n’assurent pas la fonction d’antirotation. La charge doit alors
assurer cette fonction et être capable de reprendre le couple de frottement résultant entre
l’écrou et la vis. C’est une autre raison pour laquelle il est alors très important d’incorporer cet
effet au sous modèle de la charge.
Le modèle de la vis à rouleaux que nous proposons est divisé en deux parties : tout d’abord un
écrou rigide sans jeu avec pertes par frottement et ensuite la raideur de la transmission k t ,
rejetée du coté de la translation. La dernière peut inclure les effets de jeu pour les écrous qui ne
sont pas préalablement chargés (typiquement 20 µm). Comme mentionné précédemment, la vis
à rouleaux est assimilée à une vis à filet rectangulaire avec un coefficient de frottement
équivalent. Par conséquent, un modèle de perte simple, équation (2- 8), est proposé en se
basant sur les rendements direct, inverse et pseudo de l’équation (2- 4). La variable
ψ = sgn( FV ) définit le quadrant opérationnel (+1 pour les charges résistantes et –1 pour les
charges aidantes) tandis que δ introduit la réversibilité de la vis (+1 si réversible et –1 dans le
cas contraire). La transformation couple / force est alors modélisée au moyen de l’équation
suivante :
l ⎡1 + ψ 1 − ψ δ ⎤
C=F + δηi ⎥
2 π ⎢⎣ 2ηd 2 ⎦ (2- 8)
Pour être réaliste, la valeur de δ est liée au rendement direct dont la valeur permet d’identifier
la réversibilité de la vis (η d < 0.5 ==> vis irréversible). En plus, la fonction signe, introduite
par Ψ et discontinue aux alentours de la puissance nulle, produit une discontinuité qui est
incompatible avec une bonne simulation. Nous la remplaçons alors par la fonction tangente
hyperbolique qui lisse la courbe et assure la continuité autour de la puissance nulle,
ψ = tanh( FV / P0) avec P 0 [W] comme puissance de seuil.
Une description plus interne, basée sur le principe fondamental de la dynamique (Fig 1- 11), du
fonctionnement des vis-écrous introduit explicitement le coefficient de frottement µ dans le
modèle dont l’équation devient :
d ⎡ tg ϕ + μψ ⎤
C=F
2 ⎢⎣ 1 − μψtg ϕ ⎥⎦ (2- 9)
Cette fois ci, nous sommes capables de reproduire l’influence de la vitesse de glissement V g
[m/s] sur le frottement. Ceci permet alors d’inclure les effets visqueux, de Coulomb et de
Stribeck ([ 26] et [ 27]) sur le coefficient de frottement µ par la relation:
μ = μC + (μ S − μC ) exp ⎡⎣ −(Vg / Vs ) 2 ⎤⎦ + f v VS
(2- 10)
d cos ϕ
Vg = ωve
2 (2- 11)
Avec :
Le premier graphe présente la vitesse angulaire de la vis et le chargement appliqué à son écrou.
Le rendement correspondant est donné par le second graphe. Il est situé entre le rendement
direct et le rendement indirect à l’exception des zones où il est altéré par l’effet de Stribeck au
voisinage de la vitesse nulle. Comme le coefficient de frottement µ est combiné avec le signe
de la puissance Ψ dans l’équation (2-10), le modèle n’est pas capable de reproduire les pertes
par frottement à vitesse nulle. Cependant, l’effet de séparation au voisinage d’une vitesse nulle
peut être reproduit par un choix approprié de P 0 et V s . Si nécessaire, des modèles de
frottement plus complexes comme ceux de Dahl [ 28] ou Lu-Gre [ 29] peuvent être utilisés afin
de produire un frottement statique pur. Il faut juste se rappeler que rendre le modèle plus
complexe entraîne l’utilisation de paramètres de frottement additionnels, souvent inconnus.
3. Modélisation du variateur
De nos jours, les variateurs industriels offrent une multitude de stratégies de commande :
couple, accélération, vitesse ou position, tout en incluant des filtres et des contrôleurs PID. Ils
peuvent également englober des interfaces graphiques pour faciliter l’interface homme-
machine.
« Le variateur DPM 50/80 peut être modélisé comme un système E/S dans lequel on entre une
référence correspondant à une demande en courant 10 V = 50 A et duquel on sort un courant
proportionnel à la consigne. La bande passante du variateur est de 600 à 800 Hz, ce qui est
équivalent à un filtre passe bas du second ordre. Le coefficient d'amortissement est de 0,65. »
Une autre caractéristique intéressante mais également mal documentée et peu définie est le
pouvoir de bénéficier d’un courant impulsionnel élevé. Sa valeur est quasiment le double du
courant maximal et ceci pendant une durée d’une seconde suivie obligatoirement par une
période de refroidissement des composants électroniques. Mais une incertitude persiste sur la
manière avec laquelle se fait la saturation et par la suite le passage du courant impulsionnel
vers le courant maximal permanent.
Tenant compte de ce qui vient d’être dit, le variateur sera alors assimilé à un amplificateur de
puissance qui reçoit une consigne de la part des calculateurs et délivre la puissance désirée
avec la dynamique indiquée ci-dessus. Cette consigne peut représenter la tension ou le courant
désiré ce qui nous pousse à envisager deux modèles distincts comme le montre la Fig 2- 15.
Dynamique
Tension U
désirée électrique du
moteur
Limitation de tension
-540 < U < 540
Courant
I’ I
Courant
Transformation
électrique Î mécanique
Limitation de courant
I < Imax I < Iimp
Couple
Courant ωn ξ
désiré
Dynamique équivalente
du variateur
Nous avons alors dans le premier cas une limitation de tension suivie par la dynamique du
moteur tandis que dans le second cas le variateur est modélisé par le filtre passe-bas proposé
par le service technique. A ces modèles s’ajoute une limitation du courant qui alimente le
moteur afin d’éviter un échauffement excessif de celui-ci et des composants électroniques.
La limitation de tension n’a rien de compliqué, elle est ici modélisée par une saturation à
l’extérieur de la bande –540 à 540 V. Par contre la limitation de courant doit tenir compte de
deux limites : le courant maximal autorisé et le courant impulsionnel tout en interdisant des
nouvelles impulsions avant un certain temps de refroidissement.
Pour la modélisation de ceci, plusieurs saturations, intégrateur et « switch » ont été utilisées. Le
graphe ci-dessous (Fig 2- 16) montre le fonctionnement du limiteur de courant lorsqu’il est
sujet à un courant sinusoïdal dont l’amplitude est supérieure au courant impulsionnel limite.
On peut remarquer que lorsque le courant dépasse l’intensité maximale, il continue à croître
pendant un certain temps (limité dans cette simulation à 0.2 s), il est alors considéré comme un
courant impulsionnel. Une fois la durée maximale de l’impulsion atteinte, le limiteur sature le
courant fourni au courant maximal pendant le refroidissement des composants électroniques (la
constante thermique du moteur est elle beaucoup plus élevée, 30 mn environ). Lorsque cette
durée est écoulée (limitée dans cette simulation à 1s), le variateur permet au courant de sortie
de dépasser le courant maximal tout en restant inférieur ou égal au courant impulsionnel et le
cycle recommence de nouveau.
0.2 s 1s 0.2 s
I’
On cherche ici à évaluer les dynamiques naturelles selon que l’EMA est commandé en tension
ou en courant. La capacité d’un actionneur, à opérer comme un système de chargement
dynamique, est facilement évaluée, dans un premier temps, quand il est supposé asservir en
effort une charge bloquée (impédance mécanique infinie). Ce résultat est un points-clé
important pour traiter, dans un deuxième temps, le cas le plus difficile qui est la conception de
la commande en effort sur une charge d’impédance mécanique complexe.
L, R, KT r = l / 2π
kel
Stator Rotor mt
kr Jr kve kF
Roulement Ecrou Vis-Ecrou Tige Capteur d’effort
Les dynamiques de l’EMA contrôlé en tension commencent par un premier ordre dont les
origines sont similaires à celles du vérin hydraulique. Comme pour l’EMA contrôlé en courant,
le second ordre est fixé par l’inertie du rotor et les différentes raideurs. Le contrôle en tension
introduit une raideur électrique k el dont la valeur (1.06 1010 N/m) est 70 fois plus grande que la
raideur globale de l’EMA. Cette fois, la raideur équivalente résulte de l’association en parallèle
de k el , k r et k F série k ve . Comme la raideur k F du capteur d’effort est beaucoup plus grande que
k ve , il est important de noter que cette fois ci k r et k ve sont combinés en parallèle et non pas en
série. Une fois de plus, les modes de la tige et de l’écrou sont bien plus grands que 1 kHz.
Finalement, la fréquence critique de la performance en boucle ouverte est fixée par le mode
électromécanique.
1 1
f = r k1 / J r (7.6 Hz) f = r k2 / J r (68 Hz)
2π 2π
k 1 : raideur équivalente de k ve , k r et k F en k 2 : raideur équivalente de k el , k r et (k F série k s ) en parallèle
série k el : raideur électrique K E KT / Lr ²
VI. Conclusion
Au cours du chapitre architecture et modélisation nous avons en premier proposé une méthode
multicritère de dimensionnement énergétique pour l’adaptation moteur vis. L’approche
proposée fournit la paire optimale (d, l) sur la base de la tenue mécanique, de l’usinabilité, de
la minimisation des inerties réfléchies et du couple moteur. Dans une seconde partie, nous
avons structuré le modèle de l’EMA puis proposé une modélisation générique du système vis-
écrou (transformateur avec pertes mécaniques) intégrant en un seul modèle les effets de
quadrant de fonctionnement et les diverses natures de frottement. Pour finir, nous avons mis en
évidence les principales différences de comportement dynamique entre les actionneurs
électromécaniques commandés en tension et ceux qui sont commandés en courant. Les
résultats de ce chapitre nous fournissent des points de repères importants qui permettent
d’aborder plus en détail l’analyse du propotype d’EMA et la conception du banc d’essai
associé.
Chapitre 3
BANC D’ESSAIS
I. Généralités
Avant de procéder à la description du fonctionnement du banc d’essais ainsi que ses différents
composants, il est intéressant de rappeler l’application visée, les différentes hypothèses et les
grandeurs intéressantes à mesurer. Cet ensemble d’informations constitue en effet le cahier des
charges nécessaire à la conception du banc d’essai.
1. Application ou objectif
« Une application type, visée tout au long de la thèse, est le test au sol d’un actionneur
aéronautique ayant les limites suivantes : 50 kN d’effort bloqué, 30 mm/s de vitesse à vide et
25 kN à 26 mm/s comme point de fonctionnement nominal ».
Par ailleurs, l’environnement aéronautique et l’application visée imposent eux aussi leurs
propres contraintes au cahier des charges résumées comme suit :
2. Hypothèses et choix
Comme dans toute conception mécanique, nous sommes parfois contraints à faire des choix,
suggérer des solutions et des hypothèses qui nous semblent bien adaptées, physiquement et
théoriquement, aux problèmes rencontrés. Les différents choix que nous avons pris dans le
cadre de cette étude sont présentés en détails ci après.
a) Actionneur à tester
Le but principal de cette thèse est d’évaluer la capacité à tester un actionneur aéronautique en
utilisant des moyens électromécaniques. Ceci est équivalent à dire que l’important dans ce
travail est d’évaluer le fonctionnement du dispositif de test et non pas le composant à tester lui-
même. L’actionneur aéronautique à tester, vu son prix relativement cher, est alors remplacé par
un vérin hydraulique de banc d’essai (plus lourd) ayant la même course, un effort nominal
identique et un comportement dynamique semblable. Le vérin hydraulique est censé effectuer
les mêmes déplacements (vitesse et amplitude) et réagir d’une manière identique aux différents
chargements qui lui seront imposés.
b) Chargement
Les efforts aérodynamiques qui constituent le chargement réel de l’actionneur à tester sont
reproduits par un actionneur électromécanique. Il doit être capable de suivre de près
l’actionneur à tester car tout déphasage ou défaut d’asservissement se traduisent par d’énormes
efforts à cause des raideurs élevées de la chaîne de transmission. Ceci implique alors la
présence d’une connexion rigide entre l’actionneur à tester, Servo-controlled Hydraulic
Actuator SHA, et l’actionneur de chargement, EMA.
c) Inertie
La masse équivalente de l’aileron s’élève à 600 kg. C’est une masse très grande pour être
manipulée en sécurité dans un laboratoire de recherche. Il nous faut alors reproduire cet effet
inertiel par un dispositif peu encombrant et sûr. La solution adoptée pour reproduire cet effet
est explicitée par la suite.
d) Raideur d’ancrage
d’ancrage et de la structure du banc, nous avons décidé de concevoir le banc de telle manière
que sa raideur structurelle soit dix fois plus grande que celle de la raideur d’ancrage
(égale à, 1.4 107 N/m).
Les deux vérins, bien que l’un soit asservi en position (vérin à tester) et l’autre en effort
(EMA), exercent en fin de compte tous les deux des efforts importants (50 kN chacun).
L’effort créé par l’un doit être redirigé, par un dispositif mécanique, vers la tige de l’autre tout
en passant par la charge à asservir.
Le fait que nous ayons deux vérins suggère la présence de deux dispositions possibles pour la
mise en place des vérins :
La sélection a été basée sur une étude minutieuse des avantages et des inconvénients, de
chacune des architectures, présentés brièvement dans le tableau ci-dessous.
3 - Efforts quasi-linéaires et opposés pour les petites courses - Course limitée par le bras de
(< ± 25 mm) levier de la charge
- Charge légère - Variation de l’effet inertiel en
- Possibilité d’avoir des vérins de longueurs différentes fonction de l’angle de la charge
- Possibilité de remplacer un vérin par un autre de taille - Châssis soumis à un moment de
différente par modification de l’emplacement de l’équerre flexion
4 - Compacte - Course limitée par le bras de
- Charge légère levier de la charge
- Châssis chargé en traction seulement - Variation de l’effet inertiel en
fonction de l’angle de la charge
- Déformation de la charge sous
l’effet des efforts appliqués
- La reprise d’un énorme effort par
le guidage en rotation
Il est clair que dans les architectures compactes d et f, les efforts, n’ayant pas la même ligne
d’action, produisent des moments entraînant des déformations et des contraintes
supplémentaires dans la structure. L’utilisation de l’architecture d est alors quasi-impossible
vu les risques d’arc-boutement qui dans le meilleur des cas sont remplacés par d’énormes
efforts de frottement. Quant aux architectures en ligne c et e, elles éliminent le bras de levier
entre les deux tiges des vérins mais accroissent la longueur du banc qui peut désormais
recevoir des tailles différentes de vérins par simple modification de l’emplacement des tenons
arrière. A part ça, les architectures en ligne permettent aussi une meilleure accessibilité, vision
et manipulation des pièces en mouvement par opposition à l’architecture compacte f qui
cache les deux vérins.
Le choix entre ces deux architectures repose en fait sur deux paramètres relatifs aux essais à
réaliser :
- la masse de la charge.
L’architecture e convient bien alors aux essais de faibles amplitudes, afin de limiter
l’inclinaison des vérins et la rotation de la charge, dans lesquels une grosse masse inertielle est
reproduite par une inertie rotative et un bras de levier. Par opposition, l’architecture c
convient plutôt aux essais de grandes amplitudes et faibles charges.
Nos essais faisant partie de la première catégorie, nous avons choisi alors de développer et
concevoir intégralement l’architecture e représentée en 3D dans la figure ci-dessous.
Barre d’
d’inertie
EM
A Raideur d’
d’ancrage
Stru
ctur
e po SVHHA
rteu
Équerre de se
fixation
Capteur
d’effort Axe de
rotation
Nous notons la présence des équerres qui fixent les vérins parallèlement à la structure porteuse
ainsi qu’une pièce, nommée arbre de rotation, qui sert à lier les tiges des vérins et forme avec
la barre d’inertie la charge du système.
1. Structure porteuse
Comme spécifié précédemment, la structure porteuse doit tenir les deux vérins tout en étant dix
fois plus rigide que la raideur d’ancrage. La structure est en effet soumise aux efforts des vérins
qui se traduisent essentiellement par un moment de flexion et un effort de traction le long de la
poutre comme le montre la figure suivante dans le cas d’un chargement statique.
z Equerre de fixation
x y F F F F
a
F Structure porteuse F
F
M=aF M=aF F
Fig 3- 3 : Efforts agissant sur la structure
Les points de fixations des vérins étant situés à une distance a de l’axe neutre, ils sont
influencés par la moindre flexion de la structure. Afin de respecter les contraintes de raideurs
évoquées précédemment, il faut augmenter autant que possible le moment d’inertie de la
structure par rapport à l’axe des x. Pour des raisons à la fois esthétiques et fonctionnelles, nous
avons préféré utiliser deux profilés rectangulaires, soudés côte à côte, au lieu des poutrelles
IPE (poutre dont la section est en forme de I). En effet, pour un même moment d’inertie I xx , les
poutrelles IPE pèsent 30% de moins que les profilés rectangulaires sélectionnés mais en
contrepartie elles sont plus hautes, ce qui augmente la valeur de a et par la suite diminue la
rigidité, et moins larges ce qui rend difficile la fixation des équerres et de l’axe de rotation. Il
reste à noter que comme la raideur est le principal facteur dimensionnant, les contraintes
normales dans la structure ne dépassent pas les 20 Mpa, ce qui est largement inférieur à la
contrainte maximale admissible (> 200 Mpa).
Afin de ne pas créer des concentrations de contraintes sur les profilés (épaisseur 8 mm), ces
derniers sont renforcés par d’épaisses plaques en acier aux emplacements des paliers et des
équerres.
Pion de
positionnement
Positions possibles
Profilé
Profilés
rectangulaires Vis de fixation
Plaque de
renforcement
Comme le montre la figure ci-dessous, au niveau des équerres, ces plaques sont suffisamment
longues pour permettre différentes positions et par la suite plusieurs tailles et courses de vérins.
La taille du vérin plus sa demi-course peut être comprise alors entre 365 et 855 mm. Les
équerres sont positionnées au moyen d’un pion qui sert aussi comme une liaison pivot (rotation
limitée au jeu présent dans les trous de fixations des vis) ayant pour but d’assurer le
parallélisme de la raideur d’ancrage et de l’axe de rotation. Une fois que les équerres sont en
position, elles sont fixées par 8 vis sur la structure porteuse. Le calcul des vis de fixation des
équerres ainsi que les vis de fixation des paliers est présenté dans l’annexe 3.
2. Raideur d’ancrage
Pour notre présente application, seule la raideur correspondante à l’actionneur à tester est
active tandis que l’autre, présente pour des applications ultérieures, est bloquée en traction et
compression par une entretoise. L’EMA est alors ancré rigidement au banc.
En effet, l’ancrage ayant une rigidité de 1.4 107 N/m se déforme de 3.6 mm sous le chargement
maximal de 50 kN. En plus, une exigence géométrique impose que sa largeur ne doit pas
dépasser celle du banc et le point d’application de la force doit être situé au plus près de l’axe
neutre de la structure (pour limiter le moment de flexion).
< 300 mm
± 50 kN
Vérin
Le plus petit
possible
± 3.6 mm
Vu la symétrie du banc, le tenon du vérin sera fixé au milieu de la pièce formant l’élasticité qui
à son tour sera fixée à l’équerre par les deux autres extrémités. Les chargements étant bien
dynamiques et accompagnés par des efforts parasites, les fixations de type encastrement sont
les plus appropriées. Nous obtenons alors, l’équivalent de deux poutrelles, positionnées de part
et d’autre du tenon, qui travaillent en parallèle. La relation suivante permet de calculer leur
déformée (flèche) :
M ( x)
y '' =
EI ( x ) (3- 1)
Avec :
- y´´(x ) [m-1] est la dérivée seconde de la flèche suivant la longueur de la poutrelle,
- M(x) [Nm], le moment fléchissant à la section d'abscisse x,
- E [MPa], le module d'élasticité (module d'Young),
- I(x) [m4], le moment d'inertie de la section.
Pour faciliter l’usinage, la section de ces poutrelles est choisie rectangulaire dans un premier
temps. Alors le moment d’inertie est donné par la formule suivante :
b( x ) h 3 ( x )
I ( x) =
12 (3- 2)
Avec :
Devant l’impossibilité de satisfaire les trois contraintes avec des poutres à section constante,
nous avons suggéré de concevoir des poutres à section (hauteur, épaisseur) variable. Ceci
modifie la répartition des contraintes au sein de la poutre, Fig 3-6, et permet d’obtenir une
déformation plus importante sans toutefois dépasser la contrainte maximale admissible.
350
Contrainte maximale admissible, coef de sécurité = 2
300
Contrainte de Von Mises (MPa)
250
200
150
100
Contraintesadditionnelles
Contraintes additionnelle
Î déformation
==> plusgrande
déformation plus grande
50
Section constante, flèche: 1.09 mm
Section variable, flèche: 1.8 mm
0
0 10 20 30 40 50 60 70 80 90
Abscisse sur la demi poutrelle (mm)
Fig 3- 6 : Comparaison entre les contraintes des poutres à sections variable et constante
A cet effet, nous avons utilisé le solveur d’Excel pour optimiser l’épaisseur et la hauteur des
plaques en aluminium tout en respectant les contraintes d’élasticité, d’espace et de rupture.
Malheureusement, ceci n’a permis d’obtenir que la moitié de la déformation désirée. Alors,
comme le montre la figure ci-dessous, nous avons composé l’élasticité d’ancrage par deux
lames en aluminium, travaillant en série et fixées d’une manière rigide (encastrées) entre des
plaques en acier.
Equerre de fixation
Pour conclure en ce qui concerne la géométrie des poutrelles, le moment d’inertie a suivi
l’évolution du moment fléchissant. Nous trouvons alors les poutrelles épaisses au niveau des
fixations (moment fléchissant maximal) tandis qu’elles sont minces et allongées à mi-chemin
entre les bords du tenon et les extrémités de la poutrelle.
h(x)
F
b(x)
Bien que l’élasticité d’ancrage, présentée ci haut, soit rigide (3.6 mm sous 50 kN) dans le sens
principal de fonctionnement (axe des y), elle ne l’est pas dans les autres sens (5 mm sous le
poids du vérin). Pour empêcher tous déplacements transversaux (parasites) suivant les axes x et
z, un axe est fixé sur le tenon du vérin et guidé au niveau de l’équerre par une douille à billes
pour former un pivot glissant qui reprend les efforts radiaux parasites. Ce dispositif présente
l’inconvénient de créer des efforts de frottements parasites dont la valeur maximale calculée est
inférieure à 50 N, ce qui est négligeable devant les efforts mis en jeu. Ce dispositif ainsi que le
sens des axes sont présentés dans la vue 3D ci-dessous.
Plaque en aluminium
Fixation pour reproduire la
vérin raideur avion Capteur de Douille
déplacement
Axe
y
Axe + douille pour empêcher les Pion
déplacements transversaux
3. Arbre et inertie
Arbre de fixation
des vé
vérins
Arbre de
rotation
b
Buté
Butée
mécanique
L’arbre, présenté dans la figure ci-dessus a comme but principal de créer une liaison pivot
entre les tiges des 2 actionneurs. Le rayon b (82 mm) au niveau duquel sont accrochés les 2
vérins est quatre fois plus grand que l’amplitude maximale (20 mm) des déplacements de ce
vérin qui permet de limiter la variation du bras de levier (4 % maxi) et l’inclinaison des vérins
(1.5° maxi) comme le montre la figure ci dessous.
< 0.04 b
d1 d2
Le second objectif de l’arbre de rotation est de créer avec la structure du banc une butée
mécanique qui, lors d’une défaillance des butées logicielle et électrique, limite la course par un
crash afin de ne pas détériorer les butées internes des vérins.
Tige du
Tige du
vérin 2
vérin 1
α
centre de
rotation
Butée mécanique
z Fente pour amortir le choc
En effet, quatre pièces métalliques en fer doux, évidées afin d’amortir le choc en transformant
l’énergie cinétique en une déformation plastique, sont intercalées de part et d’autre de l’axe de
rotation. La hauteur z de ces pièces est fonction de l’angle de rotation maximal α qui à son tour
dépend de la course maximale envisagée.
b
m m
Centre de rotation
1 mètre
Bien que la barre avec les masselottes ne pèse que 250 N environ, la disposition des masses lui
donne une grande inertie (4 kgm2). Pendant les phases d’accélération, de décélération et à la
limite pendant les chocs, l’arbre doit transmettre à la barre d’inertie des couples énormes
pouvant atteindre les 4000 Nm. Les accouplements ordinaires se trouvent incapables de
transmettre ce couple avec un diamètre d’arbre de 55 mm. Par ailleurs, si nous augmentons le
diamètre de l’axe, ceci va nous obliger à choisir des paliers plus grands, alors b augmente et
par la suite I augmente. Nous avons résolu ce problème en utilisant un moyeu de serrage
conique spécial (Koniclamp) qui peut transmettre 6500 Nm avec ce même diamètre d’axe. Ce
type d’accouplement est formé d’une partie conique, qui une fois serrée transmet à l’arbre le
couple par friction, et d’une couronne qui se visse à la barre d’inertie. La figure ci-dessous
présente une vue 3D des différents éléments évoqués dans ce paragraphe.
Barre
d’inertie
Capteur de
Position des rotation
masselottes
Accouplement
Palier
Pour un bon fonctionnement, une facilité d’asservissement et une bonne sécurité, l’actionneur
hydraulique, formé par un vérin symétrique et une servovalve alimenté à pression constante,
doit être équipé de quelques accessoires présentés ci-dessous.
Nous avons regroupé ces fonctions sur une plaque manifold prototype, conçue par nos soins.
1. Accumulateurs
2. Capteurs de pression
L’asservissement nécessite deux capteurs de pression qui mesurent la pression dans chacune
des chambres du vérin afin d’éliminer l’erreur statique sous charge et de compenser la non
linéarité de la servovalve.
3. Volume
La fréquence du premier ordre est fixée par la capacité hydraulique C h [m3/Pa]. Celle-ci peut
être artificiellement augmentée, en utilisant des volumes morts, afin d’améliorer la réjection
des perturbations liés à la vitesse de la charge [ 39].
4. Fuite réglable
5. Sécurité
La sécurité du matériel utilisé ainsi que la sécurité des personnes présentes aux alentours du
banc d’essai est primordiale. Ceci implique la présence d’une électrovalve qui permet de court-
circuiter les deux chambres du vérin (effet d’une roue libre). De plus, deux limiteurs de
pression sont utilisés afin d’éviter les surpressions dans les chambres et de limiter l’effort
hydrostatique maximal.
Pour ne pas encombrer notre système, nous avons choisi un montage sur bloc foré. Les
électrovalves, fuites réglables et limiteurs de pressions sont ainsi de type cartouche munies de
filetages standards T-162A et T-13A. Le bloc foré s’intercale en sandwich entre la servovalve
et le vérin selon une implantation normalisée ISO 10372-04.
Comme la plaque embase n’est fixée que par les quatre vis de la servovalve, il faut alors
qu’elle soit de taille restreinte afin d’éviter d’avoir de grosse partie en porte à faux. Les divers
perçages sont alors disposés suivants deux plans, suffisamment distincts pour éviter que la
différence de pression (200 bars) crée des fissures ou des grandes contraintes. Les perçages
sont répartis symétriquement pour limiter autant que possible l’utilisation de bouchon. Les
accumulateurs de pression, étant de petite taille et de faible masse, sont fixés directement sur la
plaque embase.
Fuite
réglable
Electrovalve
de fuite
Capteur de
pression Pb
Electrovalve Vers volume Limiteur de
pression ds A Bouchon
de Va Va
Emplacement des vis
ISO 10372-04
P
Bouchon A B Bouchon
Accumulateur Accumulateur
100 bars 2 bars
1. Grandeurs à mesurer
Dans nos essais, nous effectuons principalement des asservissements en position (SHA) et en
effort (EMA). Nous en déduisons alors qu’en premier lieu nous avons besoin de :
Des capteurs supplémentaires sont utilisés pour faire des contrôles avancés ou pour des simples
visualisations. Nous pouvons citer à ce sujet :
Les capteurs utilisés mesurent des grandeurs très diverses (force, rotation, déplacement,
pression …). Alors chaque capteur a sa propre technique de mesure qui nécessite une
alimentation bien définie et fournit des signaux qui doivent être conditionnés d’une manière
standard.
Se basant sur le mode de fonctionnement des capteurs, nous pouvons les répartir en quatre sous
ensembles comme le montre le tableau ci-dessous :
Une fois conditionnés, les signaux des capteurs sont transmis via une carte d’acquisition à un
contrôleur embarqué. Ce dernier exécute les calculs des lois de commande qui y sont
préalablement programmées et contrôle les deux actionneurs. D’autre part, le contrôleur
transmet les signaux provenant des capteurs à l’ordinateur pour une visualisation et une
sauvegarde qui permet un traitement ultérieur des mesures.
La figure ci-après présente une schématisation du banc dans laquelle nous montrons les deux
actionneurs, la partie visualisation préliminaire (numérique en sortie des conditionneurs), la
partie acquisition et commande ainsi que les différents emplacements des capteurs.
Fig 3- 19 : Diverses photos du banc d’essais, conçu et mis en œuvre par nos soins
V. Conclusion
Le travail ainsi réalisé a permis de concevoir, d’intégrer et finalement de disposer d’un banc
« bien né » qui s’est avéré fort utile pour l’identification et la commande de l’EMA. A noter
que l’intégralité des équipements de mesure et de commande ont été câblés, assemblés et testés
par nos soins.
Chapitre 4
CARACTERISATION EXPERIMENTALE
I. Introduction
Dans le chapitre précédent, nous avons présenté les différentes parties du banc d’essais, leur
conception ainsi que le fonctionnement désiré. Mais, avant d’entamer la phase
d’asservissement, il est primordial de vérifier l’exactitude des différentes valeurs numériques.
Ces dernières sont soit visées pendant la conception (raideur…) soit données par les
fournisseurs des composants (inertie…). A défaut des bonnes valeurs, la détermination et
l’exploitation des valeurs réelles permettent, dans le pire des cas, de prédire les limitations
matérielles du banc d’essais (dues aux jeux, à la souplesse …).
Différentes raideurs rentrent en jeu dans le fonctionnement de notre banc d’essais. Certaines
sont structurelles, comme celles du banc et de la transmission, tandis que d’autres sont liées
aux composants comme les raideurs de l’EMA et du capteur d’effort. A ceci s’ajoute la raideur
d’ancrage dont le but est de reproduire l’élasticité de la structure aéronautique, support de
l’actionneur en essai. Les sous paragraphes suivants ont pour but alors de présenter les
méthodes de quantification de chacune de ces raideurs.
Le prototype d’EMA a été utilisé pour charger afin de mesurer la raideur k b [N/m] effective du
banc ainsi que celle de la raideur d’ancrage k a [N/m]. La Fig 4- 1 montre les moyens utilisés
pour la mesure de chacune des raideurs. La déformation de la raideur d’ancrage a été mesurée
de deux manières distinctes : comparateurs et capteur de position résistif (Fig 4- 1, a). De
l’autre coté, la raideur de la structure du banc a été mesurée par un comparateur fixé sur des
barres à profil carré (Fig 4- 1, c), coulissant l’une dans l’autre et fixées sur les équerres à la
hauteur des articulations des vérins.
La Fig 4- 2 montre que les déformations sont bien linéaires. La raideur du banc k b est bien
conforme aux exigences ci-dessus tandis que la raideur d’ancrage est inférieure de 23 %. Il est
important de rappeler que la déformation de la structure du banc présentée ci-dessous est
mesurée entre les deux équerres. Alors la raideur ressentie par chacun des vérins, au niveau de
leur chape arrière, est le double de la raideur globale du banc (1.5 108 N/m).
2. EMA
L’élasticité de l’EMA peut dégrader d’une manière significative les performances en boucle
fermée. En pré-chargeant le système vis-écrou, le fabriquant peut éliminer les jeux, mais réduit
malheureusement la durée de vie ainsi que le rendement mécanique.
La raideur de l’actionneur peut être estimée en se basant sur les catalogues des fournisseurs ([
7] et [ 18]) (raideurs de la vis, raideurs des roulements …). Afin d’obtenir des résultats plus
réalistes tout en ayant une vision externe, nous sommes cependant obligés de conduire
quelques essais [ 38]. Deux des procédures sont suggérées ci dessous :
a) Passive
En l’absence d’un moyen spécifique pour retenir la rotation du rotor, le frein de maintien en
position est utilisé pour bloquer le rotor. Ceci permet d’éteindre l’électronique de puissance et
d’éviter toute nuisance due aux interférences électromagnétiques. Dans cette approche, l’EMA
est chargé par une force externe tout en mesurant son extension par des capteurs très sensibles
aux petits déplacements. Cette méthode est limitée par la capacité du frein (9 Nm dans notre
cas). Par ailleurs, si la raideur en torsion du frein n’est pas suffisamment importante, elle peut
amener à sous estimer la raideur d’EMA.
b) Active
Nous pouvons noter que la raideur globale est trois fois moindre aux alentours de l’effort nul.
La raideur mesurée présente alors une zone de mouvement perdu qui est représentative d’une
certaine précharge partielle du système vis-écrou. Ceci peut énormément influencer la réponse
sous faible amplitude.
Pour cette raison, il est important d’inclure cette non-linéarité de raideur dans le modèle de
l’EMA, utilisé par la suite pour réaliser un modèle virtuel du banc sous AMESim. Nous
proposons la formulation mathématique suivante pour représenter cet effet :
x0
F = ( k0 − k1 tanh )x
x (4- 1)
Où :
L’élasticité de l’environnement de l’EMA peut elle aussi dégrader d’une manière significative
ses performances en boucle fermée. Cet environnement élastique est composé de la structure
du banc, dont la mesure de la raideur a déjà été présentée, et des différents composants de
fixation (chapes des vérins…) et de transmission (embout à rotule…). Les raideurs de ces
composants étant difficiles à mesurer d’une manière individuelle, nous avons bloqué l’interface
d’ancrage du vérin à tester et utilisé l’EMA pour charger l’ensemble. Le montage, à
l’exception du comparateur, est identique à celui présenté dans la Fig 4- 1, c. Le déplacement
relatif de la tige de l’EMA, mesuré avec le capteur potentiométrique, correspond alors à la
déformation globale du banc.
La Fig 4- 4 présente la déformation globale du banc sous différents chargements générés par
des couples EMA à profils sinusoïdaux. L’utilisation d’amplitudes différentes sert à s’assurer
de la linéarité de la déformation sous les grands chargements, tandis que la faible amplitude
permet de focaliser sur le comportement aux alentours de l’effort nul. D’un autre coté, la
différence de fréquences permet de déceler, s’il a lieu, l’amortissement structural.
Comme conclusion, nous pouvons signaler que la raideur locale à effort élevé k gb [Nm] est
bien supérieure (500 fois) à celle qui règne aux alentours de la position nulle. La déformation
mesurée décèle alors une zone de jeu de 0.36 mm autour de la force nulle qui est due au cumul
des jeux (3 liaisons pivots et une rotule). Ceci peut énormément influencer le fonctionnement
sous faible chargement, créer des chocs et du bruit lors du changement du signe de l’effort et
augmenter les dépassements maximaux pour les échelons partant d’un effort nul. Par ailleurs,
on observe que l’amortissement structural est faible (peu d’effet de phase en fonction de la
fréquence).
Pour ces essais, nous proposons un modèle de représentation similaire à celui utilisé pour la
raideur de l’EMA :
x0
F = ( k0 − k1 tanh )x
x (4- 2)
L’identification paramétrique fournit les paramètres suivants :
Notre prototype d’EMA, Fig 4- 5, est composé d’un variateur 58/80 ampères et d’un moteur
synchrone d’une puissance de 6 kW qui actionne directement une vis a rouleaux inverse (pas
de 3 mm, 30 mm de diamètre). Avec ses 170 mm de course, il est capable de développer 50 kN
comme effort permanent bloqué et une vitesse de 210 mm/s à vide. La puissance maximale
qu’il peut fournir est de 5.1 kW (35 kN @ 146 mm/s). Un frein de maintien en position est
inclus pour des raisons de sécurité, afin de tenir la charge en place lorsque le moteur n’est pas
alimenté.
350 mm
1. Approche de modélisation
Dans les applications nécessitant une commande en position, la boucle de contrôle de haut
niveau est appliquée à une variable d’état de faible dynamique (position, intégrale double de
l’accélération). Ainsi, pour les asservissements en position ordinaires, les positions simulées et
expérimentale peuvent être très proches l’une de l’autre, même si les variables d’états liées à
l’accélération sont mal prédites par le modèle (ceci est dû à l’effet de filtre exercé par la double
intégration). Contrairement à ceci, l’asservissement en effort est appliqué à des variables d’état
de grande dynamique. Alors, la précision et la performance de prédiction est beaucoup plus
sensible aux erreurs de modélisation.
Quand la technologie hydraulique est utilisée pour effectuer un chargement dynamique, l’effort
de modélisation se focalise sur le frottement des joints et le comportement de la servovalve au
voisinage de l’ouverture nulle. Par ailleurs, grâce au grand gain en pression des servovalves, la
prédiction des performances est peu dépendante de la modélisation des frottements. A
La modélisation des frottements est bien documentée dans [ 26] mais l’identification des
paramètres et des irréversibilités est généralement omise et peu traitée car la plupart des
applications s’intéressent à l’asservissement en position. Ceci s’applique spécialement aux vis
à rouleaux dont les seules informations disponibles sont fournies par le fabriquant sous forme
de rendement mécanique, sans aucune distinction entre les modes opérationnels (direct ou
inverse).
Ce modèle ne prend pas en considération la précharge axiale du système vis-écrou, qui sert à
réduire/annuler le jeu présent dans la transmission [ 7]. La précharge axiale recommandée
s’élève à F p = 35% de l’effort axial maximal. De cette manière, sous un chargement extérieur
F p , le contact sur la partie chargée aurait atteint sa rigidité maximale tandis que l’autre coté
serait complètement libre. Le couple C p [Nm] requis pour compenser cette précharge est
calculé en fonction du rapport de réduction et du rendement direct pratique comme le montre la
formule suivante :
l⎛ 1 ⎞
CP = ⎜ − 1⎟ Fp
π ⎝ ηP ⎠ (4- 3)
1
ηdP = 0.9ηd =
1 + μπd / l (4- 4)
Où :
Pendant le mouvement, la force de précharge est équilibrée par le couple moteur et agit comme
une charge F p opposée au mouvement. Dans le but d’améliorer et de mettre à jour le modèle
des pertes de l’actionneur, il est préféré d’inclure la force de précharge dans le modèle au
niveau du couple moteur. De cette manière, il est alors possible de le combiner avec les pertes
d’efforts dues aux roulements et au guidage linéaire, car de toute façon il n’y a aucun intérêt et
aucune possibilité pour l’utilisateur à distinguer chaque source de perte séparément (Fig 4- 6).
D’après ce qui précède, le couple moteur global C 0 [Nm] sous charge nulle, incluant la force
de précharge, est modélisé par l’équation (4- 5) comme une fonction de la vitesse du rotor et de
cinq paramètres relatifs aux différents types de frottement : Coulomb, visqueux, quadratique et
Stribeck.
Le modèle proposé dans l’équation (4- 5) comporte cinq paramètres. Des paramètres dissipatifs
comme ceci ne peuvent pas être déterminés analytiquement mais doivent être identifiés à partir
de relevés expérimentaux. Pour cela, il est très important de définir un processus adéquat
cohérent avec les contraintes des essais.
Dans notre cas, le but du travail est de développer un processus d’identification des
paramètres avec une vue de consommateur au lieu d’une approche de concepteur
d’EMA. Cette décision élimine toute possibilité d’insertion de capteurs au sein de l’EMA
et assure des mesures non intrusives. De la même manière, aucune modification de
l’électronique de puissance n’est permise. Enfin, la méthodologie proposée est basée sur
les points de vues présentés ci-dessus afin d’être applicable par les concepteurs de
système qui intègrent des composants industriels.
Le but des essais est l’acquisition de suffisamment de données pour l’identification des
différents coefficients de perte de couple. Afin de distinguer le modèle de frottement
additionnel du modèle du système vis-écrou, toutes les expériences sont conduites avec une
charge externe nulle et avec le minimum d’inertie externe. Les pertes sous charge nulle sont
alors identifiées en utilisant deux approches différentes, sachant que toutes les deux, font
fonctionner le vérin dans le mode de commande en vitesse. La consigne de vitesse est fixée en
fonction des limitations physiques imposées par la course limitée de la tige et par la tension ou
le courant d’alimentation. Différentes amplitudes et fréquences sont combinées afin d’obtenir
un vaste domaine d’identification.
a) Vitesse constante
Dans ce cas, grâce au générateur de STIMULI présent dans le logiciel de contrôle (fourni avec
l’EMA), l’EMA est forcé à suivre un profil de vitesse trapézoïdale. Le couple moteur est alors
mesuré dans la partie à vitesse constante.
b) Vitesse en sinus
Dans la deuxième approche, l’EMA est forcé à suivre un profil de vitesse sinusoïdal qui
produit une position et une accélération en cosinus. Cette expérience permet aussi d’identifier
l’inertie rotative équivalente J r dont le couple moteur est proportionnel à son accélération
angulaire.
Le flux de données sur le banc d’essai est exposé sur la Fig 4- 7. Deux problèmes majeurs
doivent être surmontés pour conduire les expériences.
La première difficulté concerne le variateur de puissance qui génère les signaux nécessaires à
la commande du rotor en position et en vitesse et qui estime également le couple
électromagnétique. Même si ces données peuvent être visionnées au sein du logiciel, seulement
un seul signal est accessible de l’extérieur par une sortie analogique. Pour cette raison, nous
avons décidé de ne pas séparer le variateur de puissance de la partie électrique du moteur car
l’utilisateur n’est pas capable d’identifier séparément chacun des deux modèles. D’après, les
données des fournisseurs, le couple électromagnétique répond à la consigne de courant selon
une dynamique du second ordre dont la fréquence propre est voisine de 600 Hz. A l’échelle des
dynamiques de notre application, nous pouvons donc considérer que le couple
électromagnétique est proportionnel à la consigne de courant.
La deuxième difficulté est liée aux interférences électromagnétiques générées par les fortes
intensités mises en jeu ainsi que par les hautes fréquences de commutation. Une solution très
efficace peut être trouvée en synchronisant les acquisitions avec le signal d’horloge qui gère
l’électronique de commutation. Malheureusement, ce signal n’est pas une sortie standard sur
les générateurs de puissances industriels. Un câblage totalement symétrique de tous les
capteurs offre aussi une solution intéressante. Mais de nouveau, plusieurs conditionneurs de
capteurs n’offrent pas cette possibilité. Pour nos mesures, un filtrage passe-bas a été donc
appliqué (un filtre analogique RC à 330 Hz plus un filtre digital à 10 Hz après l’acquisition).
Lorsque les mesures sont prises à basse fréquence, l’influence du filtrage sur l’amplitude reste
négligeable. Par contre ce filtrage peut modifier les mesures d’une manière significative aux
alentours du changement de signe de la vitesse où le couple de frottement de Coulomb est
généralement discontinu. Les mesures effectuées à vitesse constante ne sont donc pas affectées
par ce type de problème ni par les effets inertiels. Le cas de la vitesse sinusoïdale nécessite la
procédure suivante :
• 1 è r e étape. Calcul de l’accélération
Même filtrée, la vitesse mesurée est très bruitée et ne permet aucune dérivation pour obtenir
l’accélération d’une manière utilisable (Fig 4- 8, haut à droite). Pour cette raison, nous avons
identifié le signal de vitesse pur à partir de celui mesuré en utilisant la méthode des moindres
carrés. La position et l’accélération sont alors obtenues par intégration ou dérivation analytique
de la vitesse sinusoïdale idéale.
• 2 è m e étape. Compensation du déphasage
Le couple moteur, équation (4- 5), doit être calculé en utilisant la vitesse effective. Ceci est
accompli en ajoutant une avance de phase à la vitesse mesurée qui compense le déphasage créé
lors de la mesure. Puis, le couple calculé est filtré afin d’être comparé au couple mesuré.
Les résultats expérimentaux les plus intéressants pour les pertes de couple sous charge nulle
sont résumés sur la Fig 4- 8.
La figure du haut à gauche présente les essais à vitesse constante tandis que les trois autres sont
relatives à des essais effectués avec des vitesses sinusoïdales. Les résultats expérimentaux sont
tracés en vert et les mesures réelles, non soumises au déphasage sont tracés à titre indicatif en
bleu (nous pouvons noter la discontinuité du couple aux alentours de la vitesse nulle).
Les cinq paramètres de frottement ainsi que l’inertie des parties tournantes de l’EMA ont été
identifiés en utilisant la méthode des moindres carrés, appliquée aux pertes de couple mesurées
d’un coté et simulé/filtré de l’autre.
Les tracés, présentant les mesures effectuées à vitesse sinusoïdale, montrent que cet unique
quintuplé de paramètres donne une prédiction précise pour les différentes amplitudes et
fréquences.
Par contre, nous pouvons remarquer sur le graphe du haut à gauche, que cet ensemble de
paramètres sous estime les pertes de couple à vitesse constante. En augmentant de 70 % le
coefficient de frottement visqueux pour ces mesures, nous obtenons des résultats satisfaisants,
tout en laissant inchangés les autres paramètres. Il est aussi intéressant de noter sur ce tracé que
le coefficient de Stribeck est négatif. Ceci est confirmé par la valeur de C 0S bien que ce soit
une situation peu connue. On observe également cette particularité (au pic de démarrage près)
sur les réducteurs Harmonic-Drive [ 45]. Lié aux équations (2- 9), (2- 10) et (2- 11), le couple
de démarrage C 0C + C 0S pourrait servir à l’estimation de la précharge axiale. En considérant un
coefficient de frottement équivalent pour les systèmes vis-écrou égal à 0.01 (valeur usuelle qui
mène à un rendement direct de 0.76 et un rendement indirect de 0.68), nous obtenons F p = 3
kN qui représente seulement 6 % de l’effort externe maximal.
Afin de prédire d’une manière exacte la réponse en effort, le modèle de l’EMA doit être le plus
représentatif du prototype d’essais. Divers essais sous chargement ont été réalisés sur des
charges bloquées ou mobiles. Le but est de déterminer le modèle de vis (fonction de l’effort
Les rendements de la vis ont été recherchés de deux manières différentes et distinctes qui
constituerons les sujets des paragraphes suivants.
a) Charge élastique
Une expérience typique qui pointe cet effet est représentée sur la Fig 4- 9. Une commande en
couple en boucle ouverte (20 Nm statique + 10 Nm sinusoïdal @ 50 rd/s) est appliquée à
l’EMA qui est connecté directement à la raideur d’ancrage du vérin de commande de vol. La
réponse idéale de l’EMA (sans pertes dues aux frottements et avec l’effet inertiel) est tracée
comme référence. Nous notons deux zones sur le graphe. La première correspond à la réponse
à l’échelon de 20 Nm, filtrée par un premier ordre afin d’éviter le choc dû aux jeux lors du
démarrage qui pourra résulter en une valeur finale différente. La deuxième correspond à
l’application du couple en sinus.
• Modèle a) Coefficient de frottement unique
Le modèle de pertes dans la transmission est pris des équations (2- 10) et (4- 5) tout en utilisant
les paramètres identifiés dans le paragraphe précédent (III-3). Le seul paramètre qui reste à
identifier est le coefficient de frottement global µ du vis-écrou.
• Modèle b) Coefficients liés au quadrant de
fonctionnement
Dans ce cas, le coefficient µ prend deux valeurs différentes qui dépendent du type de
chargement (aidant, résistant).
La Fig 4- 10 présente le modèle développé dans l’environnement AMESim pour effectuer les
simulations mentionnées ci dessus.
b) Charge mobile
Cette fois-ci le chargement est assuré par le SHA asservi en effort, ce qui permettra de
combiner la vitesse du rotor et l’effort exercé sur la tige du vérin. Une campagne de 14 essais
consécutifs a été réalisée. Durant chacun des essais, l’EMA a accompli des va-et-vient à vitesse
50
Partie 1
40
20
10
0
1 2 3 4 5 6 7 8 9 10 11 12 13 14
D’un autre coté, le vérin hydraulique imposait un effort prédéfini qui changeait à la fin de
chaque cycle de stimuli. Cette campagne d’essais a permis alors de balayer le fonctionnement
de l’EMA sous différents aspects :
30 15
0 0
-10 -5
-20 -10
-30 -15
0 50 100 150 0 50 100 150
Tem ps (s ) Tem ps (s)
20000 15
Mesuré
Effort (N)
Désiré
10000 5
5000
0
0
-5
0 50 100 150 0 50 100 150
Temps (s) Temps (s)
La Fig 4- 12 présente les résultats d’un essai type réalisé avec une vitesse rotor de 20 rd/s
(essai 4, Fig 4- 11). Les graphes du haut permettent de mieux comprendre le fonctionnement
du générateur de stimuli. Celui du bas à gauche présente le profil d’effort désiré, en rouge, et
l’effort réel, en bleu, influencé par le sens du déplacement. Quant au dernier graphe il présente
la valeur estimée du couple électromagnétique fournie par le variateur de puissance. Il
permettra, une fois comparé à l’effort mesuré de déterminer les pertes de couple ainsi que le
rendement de l’EMA.
L’approche de modélisation présentée dans la Fig 4- 6 permet d’écrire les équations suivantes
lorsque l’EMA fonctionne sous charge :
Cm = C0 + Cv (4- 6)
l 1 l
Cv = F ou Cv = F ηi
2π ηd 2π (4- 7)
Où :
Quant aux résultats, un traçage brut sans traitement supplémentaire des rendements montre des
traits discontinus dont la valeur dépasse l’unité dans la plupart des cas. En effet, le point faible
dans cette approche est le modèle de frottement. Il est influencé par plusieurs phénomènes
différents comme la température ainsi que le rodage du système vis écrou qui s’améliore au fur
et à mesure. Face à ce problème, nous avons identifié les coefficients de correction (Tableau 7)
nécessaires aux paramètres identifiés dans III-3 de telle façon que les courbes de rendements se
regroupent pour former un seul faisceau.
C 0C fv fq C 0S ωr
La Fig 4- 14 présente les deux faisceaux de rendements direct et inverse du système vis-écrou.
Chacune des couleurs différentes correspond à la vitesse de chacun des essais, conformément à
la légende présentée dans la Fig 4- 13. Ceci nous permet alors d’identifier l’influence de la
vitesse (comparaison des couleurs d’une même partie d’essais) et du réchauffement de l’EMA
(essais effectués à la même vitesse mais à des moments différents, partie 1 & 2).
Quant aux rendements tracés, nous notons le fait que les rendements augmentent avec l’effort
externe et deviennent constants lorsque ce dernier dépasse les 15 kN. Les valeurs auxquelles
culminent les rendements coïncident bien avec ceux retrouvés sous un chargement élastique
dans le paragraphe précédent (IV-1, a). Comme dernière remarque, nous pouvons signaler que
les courbes de rendement inverses démarrent avec des efforts de l’ordre de 2 ou 3 kN qui, en
réalité, délimitent une zone de pseudo-rendement que nous n’avons pas dessiné.
Zone de
pseudo-rendement
Il s’est avéré que cette approche additive, produisant des bons résultats à vide, n’était pas
satisfaisante sous charge à cause de la présence de pertes, autres que celles propres au vis-
écrou puisqu’elle nécessitait l’utilisation de coefficients de correction. Pour cette raison, nous
avons décidé d’élaborer un modèle global de l’EMA qui regroupe en son sein toutes les pertes
(moteur, liaisons, guidage et transformation de mouvement).
2. Modèle global
Dans cette partie, nous nous basons sur la même campagne d’essais, utilisées dans le
paragraphe précédent (IV-1, b), pour établir un modèle global de l’EMA incluant
conjointement la vitesse et l’effort externe de charge. Une nouvelle approche est alors suggérée
dans laquelle le couple moteur se divise en une partie dédiée pour vaincre les frottements, C f
[Nm], et une autre pour déplacer la charge, C ch [Nm], tout en considérant que la transformation
du mouvement s’effectue sans perte. Les équations du couple moteur et du rendement global
η g deviennent alors comme suit :
l
Cm = C f + Cch avec Cch = Fext
2π (4- 8)
sgn( ωmCm )
C
ηg = ch
Cm
(4- 9)
Afin de déterminer le meilleur modèle de frottement qui décrit la perte de couple dans l’EMA,
il est judicieux de regarder l’influence de l’effort externe ainsi que la vitesse de tige (ou du
rotor) sur le frottement. Pour une meilleure lisibilité et compatibilité des valeurs, l’effort de
frottement est tracé en fonction de l’effort externe pour plusieurs vitesses angulaires du rotor
EMA (Fig 4- 15).
Vu que les différentes courbes sont quasi parallèles et évoluent avec la même allure, nous
pouvons stipuler que les effets de la vitesse et du chargement sont découplés. Pour le modèle
de la vitesse, nous pouvons alors nous baser sur le modèle élaboré précédemment pour le
fonctionnement sous charge nulle. Une seule différence existe ; c’est l’absence du frottement
visqueux car la vitesse des essais réalisés (actionneur de commande de vol) est limitée au quart
de la vitesse maximale de l’EMA. Quant à l’influence de l’effort externe, il est clair que le
frottement y augmente d’une façon quasi linéaire mais la pente varie en fonction du quadrant
de fonctionnement.
Compte tenu de tout ce qui précède, nous proposons alors le modèle de frottement suivant :
Vitesse Partie 1 6786 -4069 21.45 0.201 -0.131 292 21.4 3.6
Partie 2 6462 -3979 55.6 0.21 -0.129 251 22.8 5.1
&
Partie 1&2 7590 -4702 70.55 0.218 -0.13 669 30.1 9.5
Il apparaît une grande différence pour les vitesses de rotation de référence ω r des parties 1 et 2
des essais. La seule explication à ce phénomène est l’effet de la température de l’EMA sur le
frottement de Stribeck. Ce dernier requiert alors plus de vitesse pour se saturer ce qui diminue
à basse vitesse l’effort de frottement non dépendant de l’effort externe. Par contre, cette
différence de paramètre entre les deux parties augmente considérablement l’erreur lorsque l’on
les regroupe. Néanmoins, une erreur moyenne voisine de 10 % sur la prédiction du frottement
doit être considérée comme acceptable vus les différents phénomènes qui y agissent.
Pseudo-rendement
Rendement inverse
Rendement direct
Pour une meilleure visibilité des effets de la vitesse et de la présence du pseudo-rendement, les
rendements calculés d’après le modèle de frottement identifié sont tracés sous forme lissés
pour trois vitesses différentes du moteur (Fig 4- 17). On notera sur la figure de droite l’effort
de réversibilité qui est largement dépendant de la vitesse et qui peut atteindre 20 % de l’effort
maximum de l’EMA (50 kN).
Ce paragraphe est dédié à la validation du modèle de frottement de l’EMA ainsi que les
différents modèles de raideurs (globale du banc & EMA). Il complète la construction du
prototype virtuel de l’EMA pour le besoin de la synthèse de la commande dans le chapitre
suivant.
Une approche progressive est mise en place. Dans un premier temps l’interface d’ancrage et le
vérin de commande de vol à tester sont remplacés par une barre rigide afin d’éliminer toutes
les souplesses additionnelles extérieures à l’EMA qui doit donc fonctionner sur charge
bloquée.
La partie motrice dans les deux modèles est identique. Elle est composée du moteur CC
équivalent au moteur synchrone et d’une boucle de commande en courant PI qui produit la
même dynamique de second ordre (fréquence 600 Hz & 0.7 d’amortissement) que celle fournie
par le service technique de Parvex (Chapitre 2, III-3). A ceci s’ajoute une saturation en tension
qui limite la tension fournie au moteur à la tension de bus maximale (± 540 V).
TF : l / 2π
Modèle de
frottement
Deux approches différentes sont utilisées pour les parties purement mécaniques. Dans la Fig 4-
18, les raideurs du banc (sous charge élevée) et de l’EMA sont représentées par leur valeur
linéaire équivalente (k gb et k EMA , Fig 4- 4 et Fig 4- 3) auxquelles est associé le montant du jeu
présent dans la transmission (± 0.18 mm). Le système vis écrou est remplacé par un
Moteur CC
équivalent
TF : l / 2π
Masse de la
tige de l’EMA
Dans la deuxième approche (Fig 4- 19), le banc d’essai est divisé en deux parties distinctes
présentes de part et d’autre du rotor. Ceci permet d’enrichir les dynamiques du système en
rajoutant les masses de la tige et du corps de l’EMA et permet de rajouter les raideurs en
torsion de ces dernières. Contrairement à ce qui a été fait dans la première approche, les
raideurs utilisées ici sont basées sur les modèles présentés dans les équations (4- 1) et (4- 2)
afin de lisser le déplacement aux alentours de l’effort nul. Enfin, vu la division du banc en deux
parties, sa raideur a été partagée en deux quantités (k’ = 2k) placée de part et d’autre du rotor.
Il reste encore à noter la présence, dans les 2 modèles, de trois filtres successifs, représentatifs
des divers filtrages effectués lors de l’acquisition des mesures.
2. Essais réalisés
Ce paragraphe est dédié à la comparaison des valeurs mesurées sur le banc réel avec celles
simulées sur le banc virtuel. Afin de couvrir le domaine de fonctionnement du banc, trois types
de consignes différentes ont été appliqués et sont présentées dans les paragraphes suivants :
a) Consigne en échelons
Cette consigne consiste à demander au moteur une succession d’échelons de couple. L’EMA
fonctionne donc en boucle ouverte d’effort. La boucle de couple est assurée par la commande
rapprochée de l’EMA au sein du variateur. Le signal excite les dynamiques du système qui
peuvent alors être identifiées grâce, entre autre, à la mesure de l’effort appliqué à la charge.
D’un autre coté, ce type de consigne permet de déterminer le temps de réponse T r [ms] entre le
couple appliqué et l’effort développé par la tige de l’actionneur.
Deux types différents d’échelons (Fig 4- 20) sont utilisés. Le premier consiste à appliquer une
consigne de couple donné, se stabiliser puis rebrousser chemin vers une consigne nulle, tandis
que le second consiste à passer d’une valeur de consigne à son opposé afin de traverser le jeu.
A part la présence du jeu qui apparaît nettement, deux phénomènes très importants sont à
remarquer. Le premier est l’influence du jeu sur les dépassements maximaux, noté dans la suite
DM. En fait lors du passage dans le jeu, tout le couple se transforme en énergie cinétique (plus
de frottement et d’effort externe) qui se restitue en comprimant d’avantage les raideurs. Le
deuxième est la très faible réversibilité de l’EMA sous les grandes charges qui se manifeste
dans le premier graphe où l’effort se bloque à une valeur proche de son DM, très distant de
l’effort idéal prévu. La comparaison des allures des différentes réponses indique que les deux
modèles AMESim ont des comportements très proches hormis la cassure au niveau de l’effort
nul. Une bonne concordance est à noter avec la mesure dans la partie à efforts positifs (allure,
valeur du DM et de l’effort lors du blocage) tandis que dans la partie à efforts négatifs les
résultats diffèrent sensiblement.
1 2
b) Consigne en sinus
Cette consigne demande au moteur de réaliser un profil de couple sinusoïdal centré à zéro.
Ceci permet de déterminer l’effort de démarrage nécessaire pour vaincre les frottements
statiques et de déceler s’il y a lieu une avance saccadée (Stick-Slip), présente en général
lorsque la vitesse est faible.
1 2
La comparaison des modèles AMESim (Fig 4- 21) montre toujours des résultats concordants
entre eux. L’effort mesuré présente un léger retard de phase par rapport aux résultats simulés.
Les allures dans la partie à efforts résistants sont similaires en temps de réponse et forme tandis
qu’elles diffèrent légèrement dans la partie à efforts aidants. Effectivement, en réalité, le
déblocage se produit de façon plus précoce (peut être à cause des vibrations haute fréquence
générées par le découpage de la tension de bus). Néanmoins, même avec la présence des
différentes non-linéarités, les modèles sont bien représentatifs du fonctionnement réel.
c) Balayage de fréquence
Nous avons réalisé le balayage en imposant des couples sinusoïdaux à une amplitude
constante, mais dont la fréquence décroit linéairement avec le temps, comme le montre
l’équation ci-dessous :
C* = Cmoy + A sin γt γ = γ 0 − αt
(4- 11)
Où :
Il est intéressant de noter que d’autres essais ont été réalisés avec des fréquences initiales
supérieures à 4 Hz afin de s’assurer de l’absence d’autres fréquences propres.
Quant aux essais présentés dans la Fig 4- 22 ci-dessous, le premier est centré à zéro. Les jeux
aident à l’amplification des oscillations, de sorte que l’effort dépasse, si l’on ne prend pas
garde, la plage mesurable par le capteur (± 50 kN). L’information recueillie du capteur d’effort
étant tronquée pour cet essai, nous avons utilisé l’oscillation de la position du rotor pour
déterminer la fréquence à laquelle l’amplitude atteint son maximum. Pour éviter ce problème,
le second essai est centré à 5 Nm et son amplitude est réduite afin que le couple oscille entre 0
et 10 Nm. Ce qui est le plus remarquable dans cet essai est l’allure quasi plate de la partie basse
de la courbe d’effort (non symétrique à la partie du haut qui est courbée). Ceci est dû en fait au
choc qui se produit lors du passage dans le jeu. Comme l’énergie emmagasinée est faible, le
système ne fait que rebondir avant que le couple ne le rappelle de nouveau dans l’autre sens.
Les fréquences auxquelles se produisent les grandes oscillations sont identiques dans les deux
essais (3.3 Hz environ). Il reste à noter que ces fréquences de résonnance sont inférieures à
celle obtenue en assimilant le système à un ensemble masse ressort (à partir de k gb ) du fait de
l’augmentation de la période due à la perte de raideur dans le jeu.
A première vue, lorsque le second essai (Fig 4- 22, 2) est comparé à l’essai simulé nous notons
que les résultats concordent peu. Une brève réflexion permet de déceler la source du problème.
En effet l’essai de balayage en fréquence et les essais en échelons sont réalisés sur une charge
bloquée dont la position n’a pas changée entre deux essais. Sous l’effet des chocs la graisse,
présente entre les différents composants du vis-écrou, est alors chassée ce qui augmente les
efforts de frottement. En augmentant le frottement de Coulomb de 7590 N à 8700 N, soit une
augmentation de 15%, dans le modèle de frottement préalablement élaboré nous notons que les
réponses coïncident de nouveau (Fig 4- 23, 1).
La Fig 4- 23, 2 prouve l’effet de la température sur l’amplitude de l’effort lors du balayage. En
effet, au fur et à mesure que les essais se déroulent, la température du moteur et du système de
transmission augmente, ce qui diminue entre autre, la viscosité de la graisse et entraine une
diminution du frottement visqueux. Comme le montrent les agrandissements a et b effectués
aux niveaux des parties extrêmes, la succession des essais entraine un accroissement des
amplitudes des oscillations. Cet effet prend fin lorsque le système atteint sa température
normale de fonctionnement (voir Fig 5- 23 dans le chapitre commande).
Pour conclure cette partie, il ressort que les essais réalisés et présentés permettent de valider le
modèle de frottement élaboré, pour lequel seuls quelques réglages et modifications de
paramètres ont été nécessaires en fonction des essais à réaliser. Ce modèle constitue donc une
base solide pour la synthèse de la commande qui sera traitée dans le chapitre suivant.
Le travail présenté dans ce chapitre constitue un passage indispensable entre les phases de
conception, conduites par les fournisseurs (structure de l’EMA…) et l’utilisateur (banc…), et
la phase d’asservissement constituant l’objectif de cette étude. Il a permis de s’assurer que les
pièces ou fonctions réalisées correspondent bien aux spécifications (raideurs …) et sont
représentatives de l’étude menée. D’un autre coté, il a permis d’identifier des modèles
paramétriques, macroscopiques simples, lorsque la modélisation physique est complexe à
mettre en oeuvre ou peu prédictive. En particulier, nous avons proposé et validé un modèle
paramétrique global des pertes mécaniques de l’EMA qui combine les effets de la vitesse, de
l’effort transmis et du quadrant de fonctionnement. A l’issue de ce chapitre, nous disposons
désormais :
C* Dynamique à 600Hz Cm 2π - F
+ +
l
saturations courant et tension
Filtrage :
100 Hz, 330 Hz, 500 Hz
Phase 1
Modélisation Hypothèses et Usinage, assemblage
Besoins et objectifs préliminaire conception et implémentation
Approches de modélisation
s
s ai
es
es
Rendements du vis-écrou
êm
m
Chapitre 5
ASSERVISSEMENT EN EFFORT DE
L’EMA
I. Introduction
1. Fréquence naturelle
Variateur
C* I
kEMA k ½ banc
Jr r Charge
Il est clair que sur charge bloquée, le premier mode est dû au système masse ressort formé par
l’inertie du rotor (inertie dominante) et les raideurs du système qui l’entoure. La fréquence
naturelle théorique f n-th [Hz] et l’amortissement ξ correspondant à un tel système du second
ordre sont alors facilement calculés par l’équation suivante :
1 k3 fv / r2
f n -th = ξ=
2π J r / r2 2 k3 J r / r 2 (5- 1)
Où :
- k3 [N/m] raideur de l’EMA en série avec la raideur globale du banc (2.52 107 N/m)
L’application numérique basée sur les valeurs mesurées dans le chapitre précédent donne une
fréquence de 4.14 Hz.
Etant donné que le système est fortement non linéaire (raideurs dépendantes du point de
fonctionnement), nous sommes contraints de mesurer la fréquence naturelle réelle afin de
valider ou corriger la valeur calculée. Deux des méthodes utilisées dans la dernière partie du
chapitre précédent peuvent contribuer à ce sujet, mais malheureusement chacune d’elle
présente un défaut qui la rend inefficace. En effet le balayage en fréquence est perturbé par la
faible raideur (ou le jeu) aux alentours de l’effort nul ce qui diminue la fréquence naturelle.
Quant aux réponses aux échelons de couples, les oscillations qui d’habitude permettent de
calculer la fréquence et l’amortissement sont quasiment inexistantes à cause de la faible
réversibilité du système et de l’important niveau de frottement non linéaire (Coulomb).
Ceci révèle qu’en boucle ouverte il est très difficile de retrouver expérimentalement la
fréquence naturelle du système. La solution alternative consiste donc à déterminer la fréquence
naturelle en boucle fermée et en déduire analytiquement la fréquence naturelle en boucle
ouverte. En effet, l’utilisation d’une commande en effort proportionnelle P (Nm/N) diminue
l’amortissement ce qui permet l’observation des oscillations lors des échelons d’effort.
P/r ξ
FTBF ′ = où ξ′ =
f / r2 J / r2 2 1 + P / r , ω′n = ωn 1 + P / r (5- 2)
1+ P / r + v s+ r s
k3 k3
2. Essais réalisés
Les essais réalisés en commande proportionnelle d’effort sont utilisés pour identifier la
fréquence naturelle du système.
La Fig 5- 2 montre que pour des échelons d’efforts réalisés à partir d’un effort nul, les
oscillations apparaissent mais les dépassements maximums sont supérieurs à 100% et les
pulsations propres dépendent de l’amplitude de l’échelon. Ceci traduit une forte non linéarité
due à la présence du jeu qui augmente l’énergie emportée par l’inertie et qui réduit la raideur
au voisinage de l’effort nul.
ω1=54.1 rd/s
ω1=48.3 rd/s
ω2=53.2 rd/s
ω2=46.9 rd/s
ω3=51.5 rd/s
ω3=42.7 rd/s
1 2
Il est donc important que l’essai n’impose pas la traversée du jeu si l’on veut identifier le mode
naturel dans la région à raideur constante, Fig 5- 3.
ω1=47.6 rd/s
ω1=56.6 rd/s
1 2
Dans ces conditions, les pulsations propres sont quasi indépendantes de l’amplitude de
l’échelon et les dépassements maximums restent représentatifs d’un système du second ordre.
Il persiste une évolution de l’amortissement due aux pertes mécaniques dans la vis (effet
fonction de la charge et de la vitesse). Ces résultats permettent d’identifier rapidement les
pulsations propres et amortissements apparents, par exemple, à partir de la pseudo-période et
du premier dépassement en appliquant les formules suivantes relatives aux systèmes du second
ordre. Bien entendu, on pourrait ici utiliser des méthodes de modèles, basées sur les moindres
carrés, en ligne ou hors ligne, pour identifier sur l’ensemble de la réponse et pas seulement sur
le dépassement maximum DM et la période T [ 24].
2π
T=
ωd , DM = e
−π ω2n / ω2d −1 ξ = 1 − ω2d / ω2n
, (5- 3)
Où :
Les équations (5- 2) et (5-3) permettent alors de déterminer la fréquence naturelle du système
(EMA sur charge bloquée) : 4.16 Hz pour les essais avec un gain P = 0.001 Nm/N et 4.33 Hz
pour les essais avec un gain P = 0.0015 Nm/N. Ces fréquences sont proches l’une de l’autre et
coïncident bien avec la fréquence naturelle théorique estimée au paragraphe précédent. D’un
autre coté, l’amortissement varie entre 8 et 32 % en fonction de la valeur de l’échelon. Ceci est
expliqué par l’augmentation de l’effort de frottement et du niveau d’effort avec l’effort externe.
Il reste à noter que ces essais sont fortement limités par le système physique dont la stabilité
impose une limitation sur le gain proportionnel.
Ces résultats autorisent donc à assimiler l’EMA à un système de second ordre pour la
synthèse de la commande. Il faut cependant garder à l’esprit que le fonctionnement
autour de l’effort nul réduit sensiblement la raideur moyenne et donc la pulsation propre.
Afin de réaliser la commande en effort la plus adéquate, pour un système donné, il faut tout
d’abord analyser ce dont on dispose naturellement (cotés positif et négatif). Ensuite, vient la
deuxième étape, consacrée à l’application et l’analyse de différentes approches de commande
en fonction du cas du chargement (typiquement, charge bloquée ou charge mobile).
1. Analyse du système
Cette partie est dédiée à l’évaluation des différents effets perturbateurs, l’analyse de la fonction
de transfert en boucle ouverte et la caractérisation du gain de boucle.
a) Effets perturbateurs
EMA
Comme le montre la Fig 5- 4, notre système inclut les effets de la chaine de mesure
(conditionnement, filtrage et acquisition). Ainsi, il y a trois dynamiques parasites d’ordre 2,
une dynamique parasite du premier ordre, un échantillonnage et une saturation :
En résumé, nous pouvons signaler que ces effets, bien qu’ils n’aient pas un caractère nocif,
créent des légers déphasages qui altèrent la commande. Leur prise en compte permet alors
d’améliorer le réalisme du prototype virtuel du banc.
Remarque importante : Les correcteurs que nous mettrons en œuvre ont une existence dans le
domaine des temps continu. Pour cette raison, nous préférons représenter dans le domaine
continu l’effet de la commande numérique (retard pur ½ période d’échantillonnage) plutôt que
de représenter l’actionneur tel qu’il est perçu par le calculateur dans le domaine du temps
discret.
Il est donc particulièrement important de bien maitriser ces effets secondaires pour
borner la réduction des performances qu’ils imposent à l’EMA commandé en effort.
Problème numérique dû à
l’approximation de Padé
Fig 5- 5 : Influence des effets secondaires sur la réponse harmonique (approche linéaire)
Il est important de noter qu’en addition avec ces effets présentés ci-dessus il existe d’autres qui
influent plus ou moins la commande. On donne comme exemple la quantification des signaux
de capteurs ainsi que la précision de ces derniers.
Des conclusions importantes pour la synthèse de la commande peuvent être tirées, compte tenu
de la forme de la fonction de transfert en boucle ouverte de l’EMA qui est de classe 0. En
commande proportionnelle pure :
- l’écart permanent, (ou erreur statique) e s , en boucle fermée à une consigne d’effort est
constant et directement lié au gain de boucle G b :
1
es = F*
1 + Gb (5- 4)
- l’écart de trainage à la rampe de consigne d’effort est infini,
1
esf = F frot
1 + Gb (5- 5)
Le modèle seul de l’EMA étant du type second ordre, il ne permet pas de montrer l’instabilité
en boucle fermée. Les phénomènes secondaires (filtrage et commande numérique) ont alors un
rôle important sur la condition de stabilité absolue en boucle fermée.
Démarrage, Cm > CC
es= f (Gbi)
EMA frein
Une même consigne de forme sinusoïdale est imposée avec des gains proportionnels différents
pour les figures de droite et de gauche. Cette approche met clairement en évidence :
L’asservissement en effort sur charge bloquée, bien qu’il constitue un cas de chargement non
exigé dans le cahier des charges, permet de s’affranchir des perturbations de vitesse engendrées
par la charge (Fig 5- 7) et ainsi d’élaborer une référence de performances.
x&
Charge
Commande C* Variateur EMA F
F*
Les essais effectués à très faible vitesse (essais sinus ou en rampe sur charge bloquée) ont
montré que la tige du vérin avance par à-coups. Ceci résulte de l’interaction entre la raideur
globale k g de l’EMA, l’inertie de son rotor et le frottement de Coulomb de la vis. Cet effet est
pénalisant car il modifie l’allure des efforts requis en leur donnant une forme d’escalier (Fig 5-
9). Nous avons essayé de traiter ce problème en créant un effet de frottement visqueux par
retour de la vitesse du rotor (Fig 5- 8). Il est intéressant de noter que, vu la raideur élevée de
l’EMA par rapport aux raideurs externes, l’angle rotor θ et le déplacement de la tige x peuvent
être considérés comme proportionnels.
Vis écrou réel
Frottement x&
F
C* θ θ1 x
Jr k TF
C F
x&
Vis-écrou parfait
G
La Fig 5- 9 présente trois essais de commande en couple réalisés avec différents gains G de
retour en vitesse rotor, dont l’origine des temps a été artificiellement décalée pour obtenir une
meilleure visualisation. Les gains de retour en vitesse sont exprimés en Nm/(mm/s) afin
d’obtenir un lien direct avec la vitesse de tige de l’EMA.
0 Nm/(mm/s)
0.25 Nm/(mm/s)
0.5 Nm/(mm/s)
Efforts mesurés
D’après ce qui vient d’être dit, l’introduction d’un frottement visqueux artificiel par retour de
vitesse du rotor n’améliore pas le comportement et rétrécit le domaine de fonctionnement de
l’actionneur.
b) Commande PID
• Réglage par essai‐correction
Compte tenu des fortes non linéarités présentes dans le système à commander, une approche
numérique par essai-correction a été mise en œuvre dans une première phase. Elle constitue
une base pour la synthèse de la commande.
Dans le cas présent, nous avons commencé par deux réglages distincts correspondant aux
zones à raideur faible (aux alentours de l’effort nul, avec présence du jeu) et à raideur élevée
(effort transmis supérieur à 10 kN).
La Fig 5- 10 illustre l’évolution des réglages effectués dans la zone à faible raideur. Tout
d’abord nous pouvons remarquer que l’effet du jeu est prépondérant vu les grands DM, même
en présence d’un bon amortissement (une seule oscillation).
1 2
Il faut noter que l’erreur permanente dépend de la position d’arrêt de la tige à cause du fort
frottement de Coulomb. Il est donc très important de procéder à des essais pour des valeurs
finales et des excursions d’effort variables pour explorer toute la plage d’erreur due à ce
frottement. Ainsi, la Fig 5- 10, 1 semble indiquer que la commande proportionnelle est
capable, à elle seule, de fournir un bon résultat (e s < 10%) alors qu’un second essai (Fig 5- 10,
2) réalisé avec une consigne d’effort inférieure de 10%, montre une erreur plus grande (29 %).
Enfin nous vérifions l’effet déstabilisateur de l’utilisation de l’action intégrale qui augmente le
DM et le temps de réponse à 5 %. Ceci est accentué par la présence du jeu dont la traversée à
effort nul se traduit par une forte accélération du rotor qui emmagasine de l’énergie cinétique.
La Fig 5- 11 présente les réglages effectués dans la zone à grande raideur. Contrairement aux
essais précédents, nous notons ici la présence de plusieurs oscillations (limitées à deux par nos
critères) ainsi que la valeur de l’erreur qui est conforme à celle calculée analytiquement en
fonction du gain de boucle (24 % pour un gain de commande de 0.0015 Nm/N correspondant à
un gain G b de 3.139). Les actions D et I sont efficaces pour éliminer les oscillations et l’erreur
statique. L’EMA présente alors un temps de réponse à 5 % voisin de 90 ms, ce qui correspond
à un système du premier ordre de fréquence de coupure 5.3 Hz.
• Méthode empirique de Ziegler‐Nichols
Pour les systèmes dont la réponse est assimilable à une dynamique du second ordre, la méthode
de Ziegler-Nichols [ 25] propose une manière empirique pour déterminer en boucle fermée les
gains du correcteur PID. L’action consiste à augmenter progressivement le gain proportionnel
jusqu’à obtenir des oscillations auto-entretenues. Le relevé du gain proportionnel P 0
correspondant et de la période T 0 des oscillations permettent de définir les gains recherchés :
Fig 5- 10 & Fig 5- 11). Nous notons que le système est très énergique ce qui se traduit par de
grands DM et un temps de réponse accru. Ce type de réponse est inacceptable pour notre
application.
• Régla ge numérique du correcteur PID
Nous proposons ici de régler le correcteur PID par une approche numérique minimisant le
temps de réponse à 5% à partir de la simulation temporelle du modèle non-linéaire de l’EMA.
Un algorithme efficace consiste à balayer systématiquement les plages de réglage qui satisfont
le critère de stabilité. Pour limiter les dépassements, le critère est pénalisé lorsque le
dépassement est supérieur à la limite acceptée.
Sans limitation
F* + 5% de courant
F* - 5% Réelle, mesurée
Tr à 5% = 45 ms Avec limitation
de courant
Influence du jeu
sur les DM
L’examen de la réponse dans la zone à grande raideur (figure de droite) valide l’approche
proposée. La saturation en courant, prise en compte pour l’optimisation n’altère pas la stabilité
de la commande bien qu’elle se manifeste pendant près de 50 % de la réponse à l’échelon. La
réponse avec traversée du jeu (figure de gauche) montre toujours un fort dépassement pour les
mêmes raisons que précédemment. Bien que le modèle utilisé pour la recherche du meilleur
réglage ne prenne pas en compte le jeu, il n’y aurait aucune difficulté à l’y intégrer. Cette
approche a donc permis de diviser par deux le temps de réponse à 5 % du système, désormais
égal à 45 ms contre 90 ms préalablement.
L’action dérivée de la commande PID s’applique à l’erreur, donc de façon égale à l’effort
mesuré et à la consigne. Ceci devient pénalisant quand la consigne est de type échelon dont la
dérivée est infinie. Une alternative consiste alors à appliquer séparément l’action dérivée à la
consigne (pour régler compensation de vitesse de consigne qui diminue le trainage) et au retour
(pour régler l’action dérivée sur le retour qui améliore la stabilité). Cette solution a un temps de
réponse à 5 % supérieur de 25 % (60 ms contre 45 ms) à celui retrouvé avec la commande PID.
Malgré cette perte en temps de réponse, elle est choisie car elle évite de saturer la commande et
rend le système plus stable (pas de dépassement lors des échelons sous grande raideur).
DM différents causés
Force (N) par des positions
initiales différentes
dans la zone de jeu Tr à 5%
= 60 ms
1 2 3 4 5
Ces essais couvrent les différentes zones de fonctionnement (raideur faible et raideur élevée,
amplitude faible ou élevée). L’asymétrie de la réponse du premier graphe est imputée aux
poids de l’EMA qui engendre un effort perturbateur à cause de la géométrie bielle-manivelle
adoptée pour la structure du banc (le zéro d’effort ne correspond pas à un effort perturbateur
nul). Les phénomènes d’échelle rendent ce défaut moins perceptible sur les autres graphes.
A noter que ces relevés sont effectués avec compensation d’erreur statique sans effet de
quadrant, conformément au paragraphe suivant. A l’exception du premier essai, l’erreur qui en
résulte est constante et de l’ordre de 200 N. Cette valeur représente 0.4 % de l’effort maximal
et correspond à un signal de 40 mV à la sortie du conditionneur d’effort (la valeur du bruit est
proche de 20 mV).
Compte tenu du caractère non intégrateur en boucle ouverte, l’erreur statique vis-à-vis de la
consigne est importante. Pour les gains adoptés dans les paragraphes ci-dessus, elle représente
typiquement 10 % de la valeur de la consigne d’effort. L’amélioration de la précision basse
fréquence par une action intégrale est limitée en pratique par les cycles limites non linéaires
engendrés par les fortes non linéarités statiques comme le frottement de Coulomb ou le jeu.
C’est pourquoi toute action d’anticipation, hors boucle donc sans effet sur sa stabilité, lui est
préférable. En ce sens, il est facile de compenser l’erreur permanente vis-à-vis de la consigne
sur la base de l’équation (5- 4). L’effort de consigne est alors « pré-accentué » d’un facteur
1+1/G b .
Le déplacement de la charge crée des perturbations sur la commande en effort. Tout d’abord,
l’effort de frottement, comme montré précédemment est fonction de la vitesse de charge. Donc,
dès que la vitesse de charge dépasse une certaine valeur, l’effort de frottement ne correspond
plus à un frottement de Coulomb ce qui rend le gain de boucle variable en fonction de la
vitesse. D’un autre coté, un mouvement de charge avec vitesse variable nécessite l’accélération
ou la décélération du rotor du moteur.
Processus à
commander x&
? Frottement
F frot
Charge
ε Commande, C* Système de 2nd ordre, + +
+- F + +
F* H(s) gain statique 2π/l idéal F F
La Fig 5- 14 présente une schématisation du système commandé, sur laquelle on trouve les
notations suivantes :
- ε [N] : erreur statique entre la consigne d’effort F* [N] et l’effort mesurée F [N]
- F idéal [N] : effort idéal produit par le vérin en l’absence des perturbations externes
- H(s) : transfert de la commande appliquée
⎧ε = F * − F
⎪
⎨ Fidéal = εH ( s )2π / l
⎪F = F (5- 8)
idéale + F frot + ?( x )
⎩ &
En négligeant l’influence de la vitesse de charge nous pouvons écrire l’équation ci-dessous qui
donne une estimation ε̂ de l’erreur statique, fonction de la consigne d’effort et de l’effort de
frottement Fˆ frot estimé grâce au modèle élaboré et validé dans le chapitre 4.
F * − Fˆ frot
εˆ =
1 + Gb (5- 9)
Afin de confirmer l’hypothèse qui nie l’influence éventuelle d’une perturbation liée à la vitesse
de charge autre que celle dépendante et engendrée par l’effort de frottement, des essais sont
réalisés à vitesse constante (Fig 5- 15) pour s’affranchir de l’accélération des différentes
inerties. La commande utilisée est proportionnelle de gain 0.002 Nm/N. La consigne d’effort
non nulle permet de se tenir à l’écart de la zone de jeu.
Fonctionnement
moteur
Fonctionnement Fonctionnement
frein frein
Fonctionnement
moteur
Les graphes de la Fig 5- 15 montrent que l’erreur statique estimée concorde bien avec l’erreur
mesurée dans les quatre quadrants de fonctionnement. Ceci confirme que la commande du
moteur n’est pas sensible à la vitesse de charge. Ces essais sont limités à 50 % de la vitesse
nominale exigée (30mm/s) à cause de la présence de bruit (pas d’amortissement dans la
commande utilisée pour cet essai).
b) Structure du correcteur
T
UT
x
F* + + x*
+- S + EMA Charge SHA
UR F
R
• Correction par anticipation
2π
Dynamique du rotor : (Cm − J r ω
&) − F frot = F
l (5- 10)
⎧ F = F * −ε
⎪
⎨ C m = C * + Ca
⎪C * = H ( s ) ε (5- 11)
⎩
Où :
2π 2π
ε(1 + H ( s )) = F * + F frot + ω
& Jr − Ca
l l (5- 12)
Alors, afin d’obtenir une erreur statique nulle, on en déduit qu’il faut appliquer un couple
d’anticipation de la forme suivante (Fig 5- 17) :
l
Ca = ( F * + Fˆ frot ) + J rω
&
2π (5- 13)
Rappel : F frot = ⎡ Fc + Fs e( − ω / ωs ) + Fext ( c + d sgn( ωF )) ⎤ sgn( ω)
⎣ ⎦ (5- 14)
Il est important de noter que le frottement utilisé pour réaliser ce couple d’anticipation doit être
calculé en fonction de la consigne d’effort F* afin de rester en harmonie avec l’hypothèse qui
stipule que l’erreur statique est nulle. Par ailleurs, l’utilisation de F* à la place de F évite un
nouveau bouclage du système qui pourrait altérer sa stabilité.
Dynamique
Anticipation (correcteur [T]) équivalente du SHA
Modèle de frottement x̂& d 1 x*
dt 1 + τs
Fˆ frot
F* + l ω̂
& d2 &x&
+ ++ Jr SHA
2π dt 2
2 F frot
⎛ l ⎞ Frottement x& +
Jr ⎜ ⎟
⎝ 2π ⎠ Charge
Ca
C* + Cm Système de 2 ordre
nd
ε H(s) -
+- + gain statique 2π/l ++
F* F
Processus
Il faudra néanmoins noter que si le SHA ne répond pas comme prévu (pression ou débit
insuffisants) à la consigne de commande en position, l’anticipation devient inadaptée et pourra
endommager le système. Ceci peut être évité en comparant en temps réel la position estimée et
la position réelle du SHA afin de déceler tout disfonctionnement et reporter le problème à
l’utilisateur sous forme de signal sonore ou autre…
• Analyse vibratoire
L’analyse vibratoire du système englobant l’EMA n’est pas de la même simplicité que dans le
cas d’un banc ne comportant que des vérins hydrauliques. En fait l’EMA présente deux
comportements totalement différents suivant que son rotor est en mouvement ou pas (Fig 5-
18) :
- Premier cas : le rotor est immobile. L’EMA se comporte alors comme une masse inerte
(sa propre masse) suspendue entre deux ressorts, qui dans ce cas agissent en parallèle,
- Second cas : le rotor est en mouvement (charge mobile). Dans ce cas l’inertie rotor
produit à travers le système vis-écrou une masse équivalente en translation de plusieurs
tonnes. Dans ce cas, les raideurs présentes aux deux extrémités de l’EMA agissent cette
fois ci en série. En effet, la rotation du rotor engendre des efforts et des déplacements
antagonistes aux extrémités de l’EMA. Les raideurs sont alors comprimées ou tendues
simultanément. Le déplacement de la tige est alors accru.
L’analyse vibratoire du banc d’essai est conduite à partir du modèle simplifié de la figure ci-
dessus.
Tout d’abord il est très important de rappeler que le comportement du système est fortement
non linéaire à cause de divers efforts de frottement :
Ces efforts étant naturellement élevés, ils amortissent énormément les oscillations produites
par les raideurs de l’EMA, d’ancrage et la raideur hydraulique du SHA. Les oscillations
décelables lors des essais, et qui peuvent dégrader nos asservissement, sont ceux liés à la
structure du banc qui, quant à elle, a un faible coefficient d’amortissement (de l’ordre de 1%).
Le premier mode mal amorti étant compris dans la bande passante (quelques Hz) de l’EMA est
atténué en utilisant un retour en dérivé de l’effort mesuré comme présenté au paragraphe
suivant. Par contre, l’EMA est incapable de rejeter les perturbations d’effort à 33 Hz (hors
bande passante). Il réagit donc pratiquement comme s’il fonctionnait en boucle ouverte. Une
solution consiste alors à éviter d’exciter ce mode, au travers de la commande de l’EMA, en
introduisant dans le correcteur série [S] un filtre réjecteur.
• Correction R‐S
Comme l’a montré le paragraphe précédent, traitant de la commande en effort sur une charge
bloquée, le correcteur R-S combinant une action série proportionnelle et un retour en dérivé
d’effort offre une solution pour une commande performante. Cette approche sera alors utilisée
aussi pour la commande en effort sur les charges mobiles. Néanmoins, les raideurs mises en
jeux et les inerties ont changé (présence du SHA et de la raideur d’ancrage). Il faudrait alors il
faudrait évaluer de nouveau les gains P et D F qui offrent une bonne réponse indicielle du
système dans toutes les plages de fonctionnement (faibles et fortes raideurs). Les nouveaux
gains sont alors : P = 0.002 Nm/N et D F = 0.00007 Nm/Ns-1 (contre 0.005 Nm/N et 0.0001
Nm/Ns-1 précédemment). La présence de la charge impose une réduction de 60 % du gain
proportionnel qui se traduit par un rallongement du temps de réponse T r de 60 ms à 90 ms.
• Conclusion sur la commande
Compte tenu de ce qui vient d’être dit et argumenté dans les paragraphes précédents le
correcteur de type RST est mis en place pour la commande en effort de l’EMA avec le SHA.
c) Essais réalisés
• Validation du choix des correcteurs
Cette première partie d’essai présente l’intérêt de chacun des correcteurs utilisés ainsi que le
comportement global du système dans toute la plage de fonctionnement ce qui permet de
valider les choix réalisés.
Ce premier essai (Fig 5- 19) présente la réponse indicielle à un échelon d’effort dans la zone à
forte raideur (5 kN à 15 kN au temps t = 6 s). La boucle PD F permet d’obtenir un temps de
réponse à 5 % de 90 ms. Cette valeur dépasse de 33 % celle obtenue avec le même type de
commande sur charge bloqué (60 ms) et correspond à une bande passante de 5.3 Hz à -3 dB
(modèle premier ordre équivalent). De plus, on superpose une pulsation de vitesse engendrée
par la commande du SHA (sinus 10 mm, 0.5 Hz). Sur la partie droite de la figure nous pouvons
remarquer l’effet du filtre trou utilisé pour éliminer l’excitation des oscillations à 33 Hz.
Tr à 5% = 90ms Oscillations à 33 Hz
Le second essai (Fig 5- 20) présente l’effet de l’anticipation basée sur les consignes d’effort et
de position (F* et x*). Cette anticipation permet de réduire l’erreur d’asservissement d’environ
5 kN à environ 500N (soit 2 % de l’effort nominal exigé) dans les cadrans moteurs et freins de
l’EMA.
Position sinus
(5 mm, f Hz) Erreur : indépendante de la fréquence
de la consigne de position f
Le dernier relevé (Fig 5- 23), présente un essai de démarrage à froid qui montre la nécessité de
prévoir un fonctionnement de quelques minutes afin que le frottement atteigne son régime
permanent. Ce problème a déjà été évoqué durant l’identification des paramètres du modèle de
frottement au chapitre précédent (la vitesse de référence de Stribeck varie en fonction du
montant des cycles effectués).
Mauvaise anticipation
de Ffrot Variation de l’effort peu significative Î l’effort de
frottement a atteint son régime permanent
Position sinus
(10 mm, 0.5 Hz)
• Respect des exigences
Cette deuxième partie d’essai vise à vérifier que l’EMA, muni de la commande préalablement
réalisée, respecte les exigences d’essais présentées dans l’introduction de la thèse. Il est à
remarquer que les exigences comportant une plage d’amplitude (±A 0 mm à ± A 1 mm) ont
seulement été testées aux deux amplitudes extrêmes (minimum et maximum).
Le premier essai (Fig 5- 24) présente la réalisation de l’exigence # 1 par l’EMA. On peut noter
un écart de 300 N (1.2 % de l’effort nominal exigé) entre la consigne d’effort (courbe bleue) et
l’effort réalisé (courbe rouge) qui dépend de la vitesse de charge. Au besoin, un facteur de
correction spécifique à cet essai (égal à 0.82) permet de réduire cet écart moyen à 50 N. Il
persiste une oscillation non amortie de 250 N d’amplitude crête à crête à la fréquence du mode
naturel de l’ensemble banc-EMA.
(0.82)
Mode naturel bouclé de l’EMA
La Fig 5- 25 est liée à l’exigence # 2 dans laquelle le chargement est extrême. Dans le graphe
de gauche l’amplitude de déplacement utilisée est l’amplitude maximale. La pression
d’alimentation du verin SHA en essai étant volontairement limitée, la position de consigne ne
peut plus être réalisée par le SHA lors des efforts importants. L’effort antagoniste produit par
l’EMA embarque le SHA qui n’est plus alors capable d’asservir la charge en position.
(0.7) (0.7)
Problème dû au grand
décalage entre la
position asservie et sa
consigne
Pour les mêmes raisons, l’asservissement de la charge en position à haute fréquence (1.3 Hz,
graphe de droite) présente un déphasage entre la consigne et la position réelle ce qui se
répercute sur l’asservissement en effort dont l’anticipation est calculée en fonction de la
consigne de position et de la réponse, hors saturation du SHA. Il n’y aurait de problème
particulier à incorporer ces limitations de puissance dans le modèle du SHA utilisé pour la
compensation (au niveau du gain de servovalve). Par ailleurs, on peut remarquer sur le graphe
de gauche que l’effort produit (courbe rouge) suit bien la consigne mais présente un léger
déphasage. En fait ce déphasage apparent est dû à un écart dépendant de la vitesse comme celui
retrouvé dans l’essai de l’exigence # 1. Il est éliminé en appliquant un coefficient de correction
de 0.7 sur l’anticipation de frottement.
La Fig 5- 26 présente l’exigence # 3 dans laquelle le chargement est modéré. Dans le graphe de
gauche on peut remarquer une bonne corrélation entre l’effort et sa consigne. Comme dans les
essais précédents, l’erreur d’asservissement reste acceptable. Eventuellement, un coefficient de
correction spécifique à cet essai permet d’améliorer encore l’efficacité de l’anticipation de
frottement. D’un autre coté, sur le graphe de droite, on peut remarquer nettement les
perturbations liées au déplacement de charge surtout lors du changement du signe de la vitesse.
Les phénomènes d’échelle rendent ce défaut moins perceptible sur les autres graphes (dans cet
essai l’erreur, voisine de 350 N, représente 17.5 % de l’amplitude de l’effort demandé mais
seulement 1.4 % de l’effort nominal).
Ce dernier tracé présente le comportement de l’EMA face à l’exigence # 4 qui est un essai à
très basse vitesse. Le plus important à noter en ce qui concerne cet essai est l’avance saccadée
(due à la basse vitesse) présentée préalablement dans l’asservissement sur charge bloquée.
L’erreur statique reste faible devant le niveau d’effort à réaliser.
• Conclusion
Les essais réalisés et présentés dans les deux paragraphes précédents permettent de valider
globalement la méthode de commande et les modèles de comportement de l’EMA. D’un autre
coté, nous pouvons noter que l’EMA arrive à bien respecter les exigences d’essais au prix de
l’introduction éventuelle de coefficients de correction. Ceci n’introduit pas une limitation
sévère dans le cas d’essais avec consignes connues à l’avance et répétitives (plusieurs centaines
de cycles). Ce réglage en fonction du type d’essai est courant et quasi incontournable, même
lorsque le chargement est assuré par un vérin électrohydraulique. De plus, il est facile de mettre
en ouvre une boucle d’adaptation lente du facteur de correction.
On cherche ici à mettre en évidence les différences entre les solutions électromécanique et
électrohydraulique utilisées comme générateurs d’effort. A cet effet, le SHA considéré est un
vérin linéaire symétrique à commande par servovalve débit.
1. Caractéristiques de puissance
Afin de réaliser une comparaison entre EMA et SHA, nous avons dimensionné un SHA
produisant les même vitesse à vide et effort bloqué que l’EMA étudié. Sur la base d’une
pression d’alimentation de 210 bars et d’un rendement mécanique de 95 %, la section utile du
vérin hydraulique serait de 26.4 cm2. La vitesse à vide imposerait un débit nominal de
servovalve de 19 l/mn. Les caractéristiques de puissance en régime permanent sont tracées sur
la Fig 5- 28.
La figure Fig 5- 29 présente les phénomènes majeurs qui déterminent le comportement statique
et dynamique des EMA et SHA. Le flux fonctionnel est tracé en gras.
I C
Jr kt
C* Ordre 2, Ce Cm F
600 Hz SE 1 ω TF 1 0 x&
l
2π
EMA R
Ffrot
Isv Ordre 2,
200 Hz R
ΔPd
ΔPv ΔPv F
SE 1 0 TF 1 x&
Ps Qv Qfonc
Qcomp
C RF
frot
SHA ∫ dt
a) Effet inertiel
Dans l’EMA, l’inertie rotor agit à travers le transformateur vis-écrou pour produire en sortie
une masse équivalente considérable de 37000 kg (à comparer à la masse de charge de 600 kg).
Dans le SHA, il n’y a pas d’effet inertiel significatif.
b) Effet capacitif
En comparaison, l’EMA est bien plus rigide avec une élasticité de transmission du système vis-
écrou 7 fois plus grande (3.2 108 N/m).
c) Effet dissipatif
Les pertes par frottement du SHA sont modérées ou négligeables selon la technologie utilisée
pour l’étanchéité (joint ou paliers hydraulique). Le rendement mécanique est alors typiquement
compris entre 90 et 99 %. Les pertes par frottement sont affectées par la vitesse de tige et les
pressions dans les chambres [ 40].
Les pertes mécaniques dans un EMA à vis à rouleaux sont importantes et entraînent une
mauvaise réversibilité. Comme le chapitre 4 l’a montré, les rendements direct et inverse restent
inférieurs à 80 et 50 % respectivement et les pertes mécaniques dépendent de la vitesse de tige
et de l’effort transmis.
D’un point de vue énergétique, la commande du SHA par valve présente un rendement
médiocre puisque la pression d’alimentation non nécessaire à vaincre la charge est gaspillée en
perte de charge au passage des orifices du distributeur, créant la dépendance vitesse-effort.
3. Dynamiques naturelles
Pour l’EMA, l’importante inertie rotor ramenée à la charge se combine aux raideurs internes ou
externes à l’EMA pour produire un mode naturel à basse fréquence : 15 Hz sur charge
parfaitement bloquée, 4 Hz sur l’ensemble banc et vérin de test. Pour la présente application,
l’élasticité interne de l’EMA a une influence très faible et peut être négligée. Il faut tenir
compte que l’EMA a deux modes naturels différents suivant que le rotor est immobile ou pas
(Fig 5- 18). Le couple électromagnétique répond à la consigne de couple avec une dynamique
du 2ème ordre à 600 Hz est insensible à la vitesse du rotor. Le transfert couple de consigne Æ
effort produit n’est pas intégrateur.
La commande en effort des SHA est bien documentée [ 39]. Les solutions proposées, Fig 5- 30,
b, ont été mises en œuvre avec profit sur de nombreux bancs d’essais d’actionneurs [ 40]. Elle
est basée sur un modèle linéaire en boucle ouverte de type 2ème ordre (servovalve) plus
intégrateur (capacité hydraulique). La performance en boucle fermée est donc directement liée
à la dynamique de la servovalve. La stabilité en boucle fermée est améliorée par un retour en
dérivée d’effort transmis et éventuellement un filtre réjecteur en série, calé sur le mode mal
amorti du banc. Eventuellement, une fuite réglable entre les chambres du vérin permet
d’amortir les modes hautes fréquences que ne peut contrôler la servovalve. A cause du
Dynamique
a équivalente du SHA
1
Modèle de frottement x̂& s x*
1 + τs
Fˆ frot
F* + l ω̂
& &x&
+ ++ Jr s2 SHA
2π 2 F frot en test
⎛ l ⎞
Jr ⎜ ⎟
Frottement x&
⎝ 2π ⎠ +
Ca Charge
ε C* + Filtre Cm Système de 2 ordre
nd
-
+- P + réjecteur gain statique 2π/l ++
F* F
Processus (EMA)
1 +DFs
b Dynamique
équivalente du SHA
Compensation vitesse
x̂& 1
x*
1 + τs
S SHA
x&
Qfonc en test
nd +
Isv Système de 2 Non Qd F Charge
ε + Non ordre : linéarité ++ s ++
+- P +
F* linéarité Servovalve tiroir F frot
inverse Capteur ΔPv 1
Filtre tiroir s
ΔP Ch
réjecteur
Processus (SHA)
1 +DFs
Le gain en pression borné de la servovalve est par contre difficile à compenser. Dans le cas
présent il introduirait une erreur d’environ 150 N pour une consigne de 10 kN. La non linéarité
dominante, produite par la caractéristique hydraulique du distributeur est efficacement
compensée par une non linéarité inverse dans l’action série. La rapidité typique obtenue en
simulation sur la présente application est caractérisée par un temps de réponse de 40 ms.
Dans ce chapitre, nous avons montré comment une approche conventionnelle pouvait être mise
en œuvre pour la commande en effort des EMA (Fig 5- 30, a). Pour la commande en effort,
l’EMA est considéré comme un système du second ordre. La dynamique du variateur et la
raideur propre de l’EMA n’introduisent pas de limitation de performance. Pour la précision en
boucle fermée, l’absence de caractère intégrateur impose l’introduction d’une action intégrale
dans la commande. Cette solution dégrade rapidement la stabilité et introduit des cycles limites
non linéaires en présence des fortes non-linéarités des pertes mécaniques. Il est donc préférable
d’anticiper l’effet des phénomènes perturbateurs (frottement et inertie rotor). Bien que le
frottement soit complexe, nous avons pu observer qu’il est reproductible même après plusieurs
mois de service. La compensation de l’inertie est incontournable. En pratique, son efficacité est
limitée par la capacité à élaborer un signal d’accélération sans bruit et sans déphasage. Hors
effet de phase, l’erreur d’asservissement typique est de 500 N pour une compensation à
paramètres invariants. Comme pour le SHA, un retour en dérivée de l’effort améliore la
stabilité tandis qu’un filtre trou permet de rejeter le premier mode mal amorti. Dans la présente
application, l’instabilité en boucle fermée apparaît à 4 Hz. Ce mode limite la rapidité en boucle
fermée et produit finalement un temps de réponse à l’échelon d’effort qui avoisine 90 ms. La
réponse à très basse vitesse est encore affectée par le mode à 4 Hz et par l’avance saccadée.
Pour améliorer d’avantage les performances des commandes plus évoluées pourraient être
mises en œuvre :
Conclusion générale
Le travail qui nous a été confié il y a trois ans traduisait une évolution des thématiques
abordées dans l’équipe de recherche « ASCOM », approche système en conception mécanique.
Il s’agissait au travers de cette thèse de développer les compétences et de se doter en moyens
d’essais concernant les actionneurs électromécaniques pour les systèmes embarqués.
Le sujet de l’étude est remonté des besoins des industriels qui s’interrogeaient sur la faisabilité
de solutions électromécaniques pour la réalisation de générateurs de force statique et
dynamique à haute puissance (quelles technologies mises en œuvre ? quels types
d’architectures disponibles ? quelles méthodologies de dimensionnement ? quelles dynamiques
naturelles et quelles conséquences sur la conception des bancs d’essai ? quelles performances
avec des commandes classiques ?).
L’application type, support de l’étude, est la génération d’effort antagoniste par vérin linéaire
pour les essais de qualification d’un actionneur d’aileron de gabarit Airbus A320. L’actionneur
d’effort sélectionné est un prototype fabriqué par SKF qui combine un moteur brushless à rotor
creux et son variateur à une vis à rouleaux inversés.
Après un premier chapitre consacré à la revue critique des technologies existantes, nous avons
étudié en détail dans le chapitre 2 les architectures fonctionnelles et géométriques ainsi que la
transmission vis-écrou. Nous avons proposé une approche originale multicritère pour le
dimensionnement préliminaire, mécanique et énergétique, du système vis-écrou et pour son
adaptation au moteur électrique. Nous avons également proposé une architecture de modèle de
l’EMA et en particulier du système vis-écrou en se référant à un modèle de vis à glissement
équivalent. Le modèle ainsi obtenu a permis d’évaluer a priori les dynamiques en présence.
La conception globale puis détaillée du banc d’essai a été présentée dans le troisième chapitre.
Le choix de l’architecture géométrique, la reproduction de l’environnement avion, la sécurité et
la sureté de fonctionnement ont piloté cette phase d’étude que nous avons intégralement prise
en charge.
système du second ordre fait clairement apparaître les effets perturbateurs linéaires (inertie
rotor) et fortement non-linéaire (frottement du système vis-écrou) ainsi que les dynamiques
parasites introduites par le filtrage des mesures et la commande numérique. Le modèle original
proposé pour le frottement prend en compte les effets de Coulomb, de Stribeck, de charge
extérieure et de cadran de fonctionnement. Tout au long de cette partie, nous nous sommes
efforcés de proposer des méthodologies d’essais non intrusives et qui peuvent être facilement
mises en œuvre par les intégrateurs de banc dans un environnement industriel.
Dans le dernier chapitre, la commande en effort de l’EMA a été abordée en mettant à profit la
progression pas à pas et les connaissances acquises dans les chapitres précédents. La validation
du modèle de conduite a été réalisée en boucle fermée proportionnelle. Elle a confirmé
l’influence du cadran de fonctionnement sur le gain statique de l’EMA et la faible fréquence
naturelle du premier mode oscillatoire due à l’effet de l’inertie rotor. Dans un premier temps, le
réglage de la commande sur charge bloquée a permis d’élaborer une référence de performance.
La synthèse de la commande en effort en présence du vérin à tester a été effectuée sur la base
d’un correcteur à structure RST. L’action série proportionnelle est dotée d’un filtre réjecteur.
L’action parallèle est de type proportionnel et dérivée. L’action T a fait l’objet d’une étude
détaillée. Elle compense en premier l’erreur statique de poursuite. La réjection de la
perturbation d’inertie est élaborée à partir de l’estimation de l’accélération de charge, calculée
par dérivation filtrée de la vitesse fournie par un modèle de l’actionneur en essai en réponse à
sa consigne de position. La réjection des perturbations de frottement est mise en œuvre grâce à
une estimation des pertes mécaniques, sur la base du modèle de frottement qui a été validé à
l’aide du prototype virtuel dans l’environnement de simulation AMESim. Pour clôturer ce
chapitre, une comparaison des solutions SHA et EMA est proposée du point de vue du
dimensionnement énergétique des phénomènes et des dynamiques naturelles et de la spécificité
de leur commande en effort.
A l’issue de ce travail, il apparaît que la prochaine étape doit être confiée à un spécialiste
d’automatique afin de développer et de mettre en œuvre des stratégies de commande plus
évoluées, tout en conservant le caractère non intrusif vis-à-vis de l’EMA et du SHA. En
particulier, il est nécessaire de s’intéresser à l’élaboration des grandeurs dérivées de l’effort et
de la position, par exemple à l’aide d’observateurs d’état et / ou de filtre de Khalman. Il serait
également intéressant de procéder à l’auto-adaptation basse fréquence des paramètres de
frottement utilisés pour l’anticipation. Le découplage de la commande en effort vis-à-vis de
l’actionneur à tester pourrait également apporter une amélioration sensible des performances.
Annexe 1
• Etude de l’isostatisme du système
Pour mieux comprendre la nécessité d’avoir un système qui reprend les efforts de rotation de la
tige de l’EMA nous allons faire une comparaison entre l’isostatisme d’un SHA typique et celui
d’un EMA.
1°_Actionneur hydraulique
2
Corps Tige Charge
1 3 ou
0
Dans le cas ci dessus, nous avons deux liaisons rotules et une liaison pivot glissant tandis que
la charge est d’habitude fixée au bâti par une liaison glissière ou une liaison pivot Î ΣL ij = 15.
Quant aux mobilités, il y a la rotation de la tige du vérin et la rotation de son corps (mobilités
internes) ainsi que le déplacement tige-charge (mobilité utile). Ceci nous amène à conclure que
le système est isostatique, h = 0.
2°_Actionneur électromécanique
3
0 2 Charge
4 ou
1
Dans ce schéma cinématique, nous avons repris le schéma précédent et remplacé l’actionneur
hydraulique par un EMA. Le rotor, solidaire de la vis forme avec son stator une liaison pivot.
Pour assurer la durée de vie du vis-écrou, son intégration dans l’EMA est telle que les paliers et
guidages externes doivent reprendre tous les efforts sur la vis, excepté l’effort axial. Dans ces
conditions, la transmission vis-écrou peut être assimilée à une liaison ponctuelle. D’un autre
coté, la tige du vérin, solidaire de l’écrou forme avec le stator une liaison pivot-glissant.
On calcul alors : ΣL ij = 21
Quant aux mobilités, le déplacement tige-charge est toujours la mobilité utile du système tandis
la rotation de la tige du vérin et la rotation de son corps ne peuvent plus être considérées
comme des mobilités internes car elles influencent le fonctionnement.
En appliquant la formule présentée ci haut, nous pouvons conclure que le système est
hypostatique d’ordre deux, h = - 2.
En effet, la rotation de la tige influence son déplacement axial. Dans le pire des cas, si la
vitesse de rotation de la tige est égale à la vitesse du couple écrou-rotor, il n’y aura plus de
déplacement axial. Il faut alors bloquer la tige contre la rotation (ceci ajoute une contrainte de
liaison). Ceci peut se faire de deux manières différentes :
- Remplacer la liaison pivot glissant par une liaison glissière entre la tige et le corps du
vérin.
D’un autre coté, le couple électromagnétique agissant entre le stator et le rotor oblige chacun
des deux à tourner dans des sens opposés. Afin que le rotor tourne avec la bonne vitesse, il faut
bloquer la rotation du stator en reprenant le couple qui y agit dessus. Ceci ajoute alors une
deuxième contrainte de liaison et rend notre système isostatique.
Ce blocage peut se faire aussi de deux manières qui y sont identiques à celles présentées ci haut
(cardan entre stator et bâti ou liaison glissière à la place du pivot glissant). Il faudra faire
attention sur le fait que l’utilisation de la liaison glissière ne règle pas les deux problèmes
simultanément car elle bloque la rotation entre la tige et le stator tandis que l’EMA reste
libre de tourner par rapport au bâti. Ceci se produira lors des phases transitoires
(démarrage, arrêt …). Lorsque le rotor change de vitesse, le stator a tendance à tourner dans le
sens contraire afin de conserver la quantité de mouvement du système.
• Etude de l’isostatisme de notre banc d’essai
4
3 6
2 Charge 7
8’
1 8
5 0
0
Fig A1- 3 : Schéma cinématique de l’EMA utilisé dans notre banc d’essai
Pour une meilleure compréhension du schéma ci-dessus nous allons commencer par énumérer
le nom et le fonctionnement des différentes pièces relatives à notre banc et présente dans ce
schéma.
- 5 : Paire de paliers, elle sert à tenir la charge et lui permettre un mouvement de rotation
(liaison pivot)
- 1, 2 et 3 : EMA
- 6 et 7 : SHA
L’équerre étant fixée au bâti par un système de pion et vis, présenté en détail dans l’annexe 3,
nous pouvons considérer que sa liaison avec le bâti est une liaison pivot (avant le serrage des
vis). Quant à la paire de paliers, sa liaison avec le bâti est considérée comme étant une liaison
plan (avant le serrage des vis).
En ce qui concerne les mobilités et liaisons des différents éléments on se réfère aux
paragraphes précédents.
Î ΣL ij = 48
Avant le serrage des vis, le système est hypostatique d’ordre deux, h = - 3. Ceci permet alors
un bon positionnement des pièces sans aucun risque de les déformer. Une fois que l’alignement
des tiges est assuré, nous serrons les vis qui fixent les paliers, la liaison plan se transforme en
un encastrement et augmente de trois le nombre des inconnues de liaison. Notre système paraît
alors isostatique. Mais en découplant le calcul des cotés gauche (EMA) et droit (SHA), nous
retrouvons que la fixation de l’EMA est hyperstatique d’ordre un tandis que la fixation de
l’SHA est hypostatique d’ordre un.
Enfin, lorsque nous serrons les vis des équerres la fixation de l’EMA devient hyperstatique
d’ordre deux et celle de l’SHA devient isostatique. Il faudra noter que l’hyperstaticité dont
nous parlons provient du fait que nous avons utilisé des liaisons pivots pour fixer l’EMA au
lieu d’utiliser des joints de cardan. C’est un choix de conception assurant une plus grande
simplicité, l’hyperstatisme étant absorbé par les différents jeux et élasticités des liaisons.
Annexe 2
RIGIDITE DE LA VIS
• Roulement fixe sur une extrémité de la vis
Dans ce premier cas, un seul ensemble de roulements assure la fixation de la vis sur une
extrémité (l’autre restant libre), ex : vérin ayant une faible course. La liaison est considérée
comme un encastrement. Les paragraphes ci-dessous donnent la rigidité du système en fonction
du type de la vis : à billes ou à rouleaux.
1°_Vis à billes
E π ( d 0 − 0.71 Dw )
2
kv =
4 n (A2- 2)
Fig A2- 1
Avec :
2°_Vis à rouleaux
E π d 22
Rv =
4 n (A2- 3)
Avec :
Dans ce second cas, deux ensembles de roulements assure la fixation de la vis aux deux
extrémités, ex : vis de positionnement du chariot d’une machine outils. La liaison, bien
qu’hyperstatique (2 encastrements), augmente la rigidité du système mais rend la rigidité
fonction de la position de l’écrou. Les paragraphes ci-dessous donnent la rigidité du système
pour les deux types de vis.
k v = k v1 + k v 2
F (A2- 4)
Fig A2- 6
1°_Vis à billes
E π ( d − 0.71Dw ) 2 ns
kv =
4 n ns − n (A2- 5)
Fig A2- 7
Avec :
( d − 0.71 Dw )
2
kv min = Eπ
ns (A2- 6)
Fig A2- 8
Comme le montre la Fig A2- 8, la rigidité est minimale au milieu de la vis : n = ns / 2 , et elle
tend vers l’infini aux extrémités.
2°_Vis à rouleaux
Eπ d 22 n
kv =
4 ns ( n − ns ) (A2- 7)
De même, la rigidité est minimale au milieu de la vis :
d 22
kv min = E π
ns (A2- 8)
Annexe 3
PRINCIPE DE FONCTIONNEMENT
D’UNE VIS A ROULEAUX SATELLITES
1. Structure
La Fig A3- 1 présente une vue détaillée d’une vis à rouleaux satellites ([ 34] et [ 35]). La vis et
l’écrou ont des filetages à entrées multiples et un profil triangulaire (angle de 90°) tandis que
les rouleaux ont un filetage unique formé par des arcs circulaires. Les engrenages, situé à
chaque extrémité du rouleau, sont engagés avec les engrenages internes situés sur les deux
cotés de l’écrou. Deux attaches assurent le bon positionnement des rouleaux (espacement
identique, axe parallèle à celui de la vis).
1- Vis 5- Attache
2- Ecrou 6- Anneau
3- Rouleau 7- Goupille
4- Engrenage interne 8- Engrenage externe Détails du rouleau
9- Fausset
2. Mécanisme
Le fonctionnement d’une vis à rouleaux satellites, est semblable à celui d’un train
d’engrenages planétaires. Cependant contrairement au train d’engrenages, la rotation des
composants d’une vis à rouleaux s’accompagne d’une translation plus ou moins importante.
a) Rotation
Considérons le cas où la vis est en train de tourner et l’écrou est bloqué contre la rotation.
Ecrou
dv [m] diamètre de la vis
de [m] diamètre de l’écrou
Vis
dr [m] diamètre du rouleau
dm [m] diamètre moyen, égal à dr + dv
En utilisant les symboles présentés dans la Fig A3- 2, la relation entre la vitesse de rotation du
rouleau ω’ et la vitesse de rotation de la vis ω est exprimée comme suit :
ω ' d e ωr d r
Ecrou fixe ⇒ − =0
2 2 (A3- 1)
ω ' d v ωr d r ω d v
Vis en rotation ⇒ + =
2 2 2 (A3- 2)
dv ω dv ω k
D'où: ω′ = ω = =
de + dv 2 dv + dr 2 k + 1 (A3- 3)
k étant le rapport entre le diamètre de la vis et celui du rouleau (k = d v /d r )
Il semble intéressant de noter que dans une utilisation réelle la vitesse angulaire ω est la vitesse
relative entre la vis et l’écrou. Ceci est équivalent à dire que si la vis est fixe et l’écrou tourne,
la relation entre ω et ω’ reste inchangée mais cette fois ci ω est la vitesse de rotation de l’écrou.
b) Translation
• Rouleaux / Ecrou
Tout calcul fait, le déplacement axial a 1 du rouleau, relativement à l’écrou, pour un tour de la
vis est donné par :
1 kp n k +2
a1 = (ω′pne ± ωr p ) = ( e ± )
ω 2 k +1 k +1 (A3- 4)
Avec :
Lorsque les rouleaux sont solidaires de l’écrou (pas de recyclage, vis longue), le déplacement
axial entre les rouleaux et l’écrou est nul. Ceci implique des filetages de sens inverses ainsi
qu’une relation entre les diamètres des rouleaux et de l’écrou.
de
ne = k + 2 =
dr (A3- 5)
• Rouleaux / Vis
Tout calcul fait, le déplacement axial a 2 du rouleau, relativement à la vis, pour un tour de la
vis vaut :
1 k ( k + 2) k
a2 = ( ωr p ± ω′pnv ± ωnv p ) = p( ± nv + nv )
ω 2( k + 1) 2( k + 1) (A3- 6)
On constate que le déplacement axial du rouleau relativement à la vis dépend de l’avancement
par tour de la vis pn v ainsi que du rapport des diamètres k. Afin de ressembler le plus possible
aux vis à glissement traditionnelles, il vaut mieux choisir n v = k + 2 de telle façon à éliminer la
composante en k. D’après l’équation (A3- 6), on remarque qu’il faut choisir le même nombre
de filets pour l’écrou et la vis. Il reste à noter que le nombre de rouleaux n’est limité que par les
contraintes d’espace, ce qui présente un bon intérêt pour la capacité de charge.
Annexe 4
MODELISATION DU MOTEUR
SYNCHRONE
dI
U=L + RI + K e ωm
dt (A4- 1)
Ce = K T I (A4- 2)
Avec :
Le choix de ce modèle n’est pas trop éloigné de la réalité vu que le moteur synchrone a un
fonctionnement semblable au moteur CC commandé par tension d’induit à flux inducteur
constant. Par ailleurs, la plupart des fournisseurs de moteurs synchrones ne présentent sur leurs
documentations que les paramètres des moteurs CC équivalents. Le tableau ci dessous, fourni
par la société Parvex [ 16] en est un bon exemple. Il présente les paramètres du moteur qui
équipe notre actionneur EMA. Ces mêmes paramètres seront utilisés dans nos simulations. Il
reste à noter que la résistance et l’inductance équivalentes sont égales au tiers des valeurs
mesurées entre les phases du moteur synchrone.
Le moteur qu’on vient de décrire peut être commandé par tension ou par courant (par vitesse
ou par couple dans l’interface du logiciel Parvex Motion Explorer). A ce stade, il suffit de
noter que la commande en courant élimine la dynamique du moteur et le transforme en un
simple gain entre le courant et le couple conformément à l’équation (A4- 2) mais d’un autre
coté, elle ajoute une dynamique au sein du variateur (passage de U à I).
Fig A4- 1 : Paramètres du moteur CC équivalent du prototype d’EMA fourni par SKF
Il est intéressant de noter que la phase transitoire du démarrage du moteur n’a pris que 0.15 s
ce qui annonce une rapidité de réponse pendant les asservissements ultérieurs.
Annexe 5
• Equerres
Chaque équerre du banc d’essais est positionnée au moyen d’un pion de diamètre 30 mm et est
fixée par deux rangées de quatre vis chacune. L’effort appliqué par le tenon arrière de
l’actionneur, n’appartenant pas au plan de contact entre l’équerre et la structure porteuse, crée
dans les vis en plus des contraintes de cisaillement des contraintes normales dont les valeurs
varient en fonction de l’emplacement des vis.
y
F F
x
c
O
F1
F2
F3
F4
F6 F6 F6 F6
F F5 F
F6 F6 F6 F6
a b b b a
Fig A5- 1 : Schématisation de l’équerre avec les efforts qui y sont appliqués
1°_Contraintes de cisaillement
Comme le montre la figure ci dessus, la reprise des efforts selon l’axe des x est assurée par le
pion et les huit vis. En écrivant le bilan des efforts statiques appliqués à l’équerre selon l’axe
des x, nous obtenons :
8 F6 + F5 = F (A5- 1)
En consultant de près la conception de l’équerre, on remarque que le pion est placé avec un
ajustement juste (dans l’équerre et la structure) tandis que les vis sont entourées d’un petit jeu
permettant de pivoter l’équerre afin d’assurer un bon parallélisme avec l’arbre de rotation.
Alors, dans le cas le plus défavorable où nous négligeons les efforts de frottement, c’est le pion
qui bloque le déplacement plan de l’équerre.
La valeur des contraintes de cisaillement agissant sur le pion est donné par :
F F
τ= = = 70 MPa
S πd 2 / 4 (A5- 2)
Cette contrainte est 3 fois plus petite que la limite de la rupture en cisaillement d’un acier de
qualité moyenne. A ceci s’ajoute le fait que le serrage des vis, fixant l’équerre à la structure
porteuse, crée un grand effort normal entre les 2 surfaces qui à son tour se traduit par un effort
de frottement qui empêche tout déplacement horizontal de l’équerre. On en conclut que le pion
et l’effet de frottement sont largement capables de tenir l’équerre en place.
Quant à la contrainte de matage sur la surface du pion, elle est voisine de 55 MPa, valeur
acceptable.
2°_Contrainte normale
La fonction principale des vis est la reprise et / ou la production d’efforts normaux. La force
agissant sur chacune des vis fixant notre équerre dépend à la fois de la force F qui lui est
appliquée et de l’emplacement de la vis par rapport à la ligne de basculement passant par O.
Compte tenu du nombre important des vis, notre système est hyperstatique. Ceci nécessite une
résolution par élément finis où des hypothèses réalistes pour réduire le nombre des inconnues.
Supposons que l’équerre et le banc sont bien plus rigides que les vis. Alors, comme le montre
la figure ci-dessous, sous l’action de la force F, les vis sont les seuls éléments qui se
déforment.
y
F
x
O
ε1 ε2 ε3 ε4
Tout calcul fait : F 4 = 33.4 kN. Comme il y a deux rangés de vis, la vis la plus sollicitée reprend
16.7 kN. Elle est alors soumise à une contrainte normale de :
F 16700
σ= = = 53MPa
S π 202 / 4 (A5- 3)
La résistance à la rupture des vis de classe moyenne (classe 8.8) est de 800 MPa. Sachant que
les vis sont, d’habitude lors du serrage, précontraintes à 70 % de leur limite de rupture (560
MPa), la contrainte maximale, produite par la reprise des efforts du vérin, ne présente que 10 %
de la précontrainte. Ceci permet alors une très bonne résistance à la fatigue et une très bonne
fixation des équerres.
• Palie rs et butée de crash
Les paliers, comme spécifié dans la description du banc d’essai, servent à la fixation de l’arbre
de rotation. Chacun des deux paliers est fixé par deux vis M16. Des épaulements empêchent
tout déplacement horizontal et évitent que les vis reprennent des contraintes de cisaillement.
Pendant le fonctionnement normal des vérins, hors arrivée en butée, les paliers sont faiblement
sollicités. En effet les efforts selon l’axe des x produit par les deux vérins se compensent (hors
effet inertiel de l’arbre de rotation) et les efforts selon l’axe des y, repris par les paliers, sont
très faibles vu l’inclinaison très réduite des vérins (maxi 1.5°).
L’arbre de rotation est conçu de manière à former avec la structure du banc une butée
mécanique qui lors d’une défaillance des butées électriques évite que les vérins arrivent à leur
propre fin de course. Le calcul des contraintes est effectué pour le cas le plus contraignant :
lorsque les deux actionneurs forcent dans le même sens.
La figure ci-dessous présente les différents efforts externes appliqués au système lorsque celui-
ci touche la buté mécanique ainsi que les interactions entre ses différents composants :
- F act1 : effort produit par le premier actionneur, appliqué sur l’arbre de rotation
- F act2 : effort produit par le second actionneur, appliqué sur l’arbre de rotation
y
Fact1 Fact2
x
Fpy
d
Fb Fpx
A
Fv1 Fv2
c
a b
Fig A5- 3 : Principaux efforts résultant d’une arrivée en buté
Sachant que:
F act1 = F act2 = 50 kN
La vis la plus sollicitée doit pouvoir alors supporter 42.8 kN (F v2 /2, car 2 paliers). Les
contraintes normales agissant dans cette vis lors de l’arrivé en buté s’élèvent à :
F 42800
σ= = = 213 MPa
A π 162 / 4 (A5- 4)
En faisant le même raisonnement que précédemment (classe …), nous pouvons considérer que
les vis M16 sont suffisantes pour fixer les paliers. Il semble important de noter que le
dimensionnement est réalisé dans les conditions extrêmes : crash et défaillance totale de toute
commande.
Annexe 6
CAPTEURS
Cette annexe présente d’une manière plus complète et détaillée les capteurs utilisés pour
effectuer les différentes mesures et acquisitions (sensibilité, étendue de mesure …).
• Résolveur
1°_ Points forts :
- Précision comprise entre 12’ et 20’ d’angle (valeurs basées sur la gamme des résolveurs
Parvex, référence précise non communiquée)
2°_Point faible :
• Capteurs potentio métriques
Ces capteurs sont utilisés pour la mesure de l’élongation de la tige de l’EMA (Mégatron,
RC20-100), de la compression de la raideur d’ancrage du vérin en essai (Mégatron, RC13-25)
et de la position angulaire de la charge (Mégatron, MCP40). Ce type a été choisi grâce à son
immunité vis-à-vis de la pollution électromagnétique.
1°_Points forts :
- Bonne linéarité, 0.1% (0.2 pour le RC13-25) de la pleine échelle (respectivement 100
mm, 25 mm et 360°)
- Résolution infinie
2°_Points faibles :
- Signal bruité malgré le blindage des câbles de mesures Î impossible de dériver les
signaux sans un filtrage poussé (filtre numérique de second ordre à 10 Hz)
• Capteur LVDT
Ce capteur de type LVDT (Sensorex, SX12W100) intégré dans le corps du SHA, mesure la
position relative tige/corps. Il y est fourni avec d’une manière standard.
1°_Points forts :
- Résolution infinie
2°_Points faibles :
- Dérivée bruitée
• Capteur d’effort
Ce capteur (HBM, U10M-50kN) fixé sur la tige de l’EMA, permet de mesurer l’effort appliqué
à la charge et d’effectuer la commande par la suite. Il est réalisé avec des jauges résistives en
pont Wheatstone complet.
1°_Points forts :
- Taux de bruit très faible < 8 mV à la sortie du conditionneur (sur une plage de ± 10V)
- Très haute rigidité axiale (30 µm/50 kN), pas d’influence sur la rigidité globale du banc
- Permet de transmettre des couples très élevés (1270 Nm), cette qualité est très
importante vu le fait qu’il fait partie de la chaîne qui assure le blocage en rotation la
tige du vérin
2°_Points faible :
• Cap teurs de pression
Ces capteurs (P901) fixés sur la plaque embase du SHA, mesurent la pression dans chacune
des chambres du vérin afin de réaliser une commande poussée. Ils sont réalisés avec des jauges
résistives en pont Wheatstone complet.
1°_Points forts :
- Taille réduite
2°_Point faible :
• Con ditionneurs
Les conditionneurs sont nécessaires pour les capteurs de types inductif (LVDT) et pont
Wheatstone (effort et pression). Ils sont choisis dans la gamme des conditionneurs Sensorex
dont les séries sont respectivement 3120 et 3310. Ils requièrent une alimentation stabilisée (24
V) et fournissent une tension de sortie standard (± 10 V). Le signal de sortie est filtré à 500 Hz.
Nomenclature
es erreur statique N
f fréquence naturelle du vérin Hz
F force axiale appliquée à la tige de l’EMA N
F' effort conditionné et filtré N
F'' effort conditionné et filtré, filtré par le filtre RC N
F* consigne d'effort N
F1 à F4 effort normaux supportés par les vis 1 à 4 N
F2 force exercée par l'actionneur à tester N
F act1 effort produit par le premier actionneur N
F act2 effort produit par le second actionneur N
F aid force aidante N
Fb réaction de la butée sur l'arbre de rotation N
Fc frottement de Coulomb N
Fe force externe agissante sur l'écrou N
F frot effort de frottement N
Fˆ frot effort de frottement estimé N
F idéal effort idéal produit par le vérin en l'absence des perturbations externes N
Fm effort moyen appliqué à l'écrou N
F max effort maximal fourni par un EMA kN
F mi effort minimal durant le cycle "i" N
F Mi effort maximal durant le cycle "i" N
f n-th fréquence naturelle théorique Hz
Fp effort de précharge estimé N
F px réaction de l'arbre de rotation sur le palier suivant l'axe des x N
F py réaction de l'arbre de rotation sur le palier suivant l'axe des y N
fq coefficient de frottement quadratique Nm/(rd/s)²
F res force résistante N
Fs frottement de Stribeck N
FTBF fonction de transfert en boucle fermée
FTBO fonction de transfert en boucle ouverte
fv coefficient de frottement visqueux Nm/(rd/s)
Fv force externe agissante sur la vis N
F v1 réaction du palier sur la première paire de vis N
F v2 réaction du palier sur la deuxième paire de vis N
Gb gain de boucle du système -
G bd gain de boucle en prenant en compte le rendement direct -
G bi gain de boucle en prenant en compte le rendement inverse -
h degré d'hyperstatisme du système -
H(s) transfert de la commande appliqué
h(x) hauteur de la poutrelle à la section d'abscisse x m
I courant A
I' courant limité A
I(x) moment d’inertie de la section m4
I eq courant équivalent d'alimentation A
Ii courant dans l’enroulement « i » A
I imp courant impulsionnel maximal A
I max courant maximal A
J inertie de la charge
Jr inertie des parties tournantes de l'actionneur kg.m²
K paramètre liée au champ magnétique du moteur -
k rapport entre les diamètres de la vis et du rouleau m/m
ηd rendement direct -
η dP rendement direct pratique -
ηi Rendement inverse -
ηp Pseudo rendement (charge aidante et transmission irréversible) -
γ pulsation de la sinusoïde rd/s
γ0 pulsation à l'origine des temps rd/s
Références
[ 26 ] : Olson H., Astrom K.J., Canudas de Wit C., Friction Models and Friction
Compensation, November 28, 1997, www.control.lth.se/~kja/friction.pdf
[ 30 ] : Chami M., The use of the CHDN to model a permanent magnet synchronous
motor powered by ultracapacitors, Vehicle power and propulsion (IEEE, Oct
2004)
[ 31 ] : Aghili F., Buehler M., Hollerbach J. M., Optimal Commutation Laws for
Torque Control of Synchronous Motors, Proceedings of the American Control
Conference, pp 1968-1973, Albuquerque, New Mexico, June 1997
[ 33 ] : Pasch K. A., On the drive systems for high performance machines, ASME J.
Mechanisms, Transmissions, Automation, Design, Vol 106, No 1 Mars 1984
[ 34 ] : Otsuka J, Osawa T., Study planetary roller screw comparison static stiffness
vibration characteristics those ball screw, Bull. Japan Soc. Of Prec. ENGG.,
Vol. 23, No. 3 (Sept. 1989), pp 217-223
[ 38 ] : Karam W, Maré J-C., Modelling and measure of the compliance and friction
losses of screwjack electromechanical actuators, Power Transmission and
Motion Control (PTMC 2007), pp 145-157, Hadleys Ltd, 2007
[ 40 ] : Maré J-C., Dynamic loading systems for ground testing of high speed
aerospace actuators, Aircraft engineering and aerospace technology
(international journal), Vol. 78 No 4, 2006, Emerald ltd
[ 41 ] : Perni F, Pizzoni L., Ballscrews for the more electric aircraft, Recent advances
in aerospace actuation systems and components, Nov 2004, Toulouse France
[ 44 ] : Plummer A., Robust electrohydraulic force control for test rigs, Recent
advances in aerospace actuation systems and components, June 2007, Toulouse
France
[ 45 ] : Han J.D., Event-based discontinuous control for harmonic drive system with
the damping enhancement of acceleration feedback, Proceedings of the 3rd
world congress on intelligence control and automation, June 2000, Heifei China
Résumé
• Titre :
• Résumé :
• Mots clés :
Abstract
18B
• T itle :
178B
Investigation into the electromechanical actuator when used for high performance force control
• A bstract :
179B
The recent progress in power electronics and power drives enable the development of
electromechanical actuators that become attractive when price and control easiness are
considered. Subsequently, they can replace hydraulic actuators in applications requiring
medium efforts (few kN) and a friendly environment. On the other hand, certification and
qualification, of aeronautical or automobile actuators, request performance tests covering their
domain of operation. This requires high performance actuators to generate the loading efforts.
The main objective of this research work is evaluating the limitation of electromechanical
actuators used for force control, with respect to power performance, control and life
expectancy. It is organized in order to provide advanced knowledge to test bench designers and
users when high performance load simulation is required.
• K ey words :
180B