2010 - LISCOUETcorr-1
2010 - LISCOUETcorr-1
2010 - LISCOUETcorr-1
En vue de l'obtention du
Titre :
Conception préliminaire des actionneurs électromécaniques- Approche hybride,
directe/inverse
JURY
Jonathan LISCOUËT
A Susan
REMERCIEMENTS
Ces travaux ont été réalisés à l’Institut Clément Ader de Toulouse sous la direction de
Jean-Charles MARÉ et Marc BUDINGER dans le contexte du contrat Européen DRESS
(Distributed Electrical nose gear Steering System) et du projet ANR C6E2 (Simulations pour
l’Ingénierie de la Conception des Systèmes Electriques Embarqués). Je tiens à remercier ici
tous ceux qui ont contribués à ce que ces travaux s’effectuent dans des conditions optimales,
motivantes et enrichissantes durant ces 35 mois de doctorat.
M. Jean-Pierre Garcia, pour m’avoir fait l’honneur de prendre part au jury auquel il a
apporté son expérience et sa vision industrielle éclairée.
Je tiens aussi à remercier les membres de l’Institut Clément Ader pour l’ambiance
amicale et conviviale dans laquelle j’ai passé ces dernières années. Par ordre alphabétique, je
tiens plus particulièrement à remercier : Annie Cazeaux, Ahmad et Battoul Batikh, Christine
Barrot, Christophe Paris, Daniel Ortega, Fabien Hospital, Feriel Samouda, Nicolas Laurien,
Romain Canivenc, Stéphane Orieux, Sandrine Geoffroy, Toufic El-Halabi et Wissam Karam.
Finalement, mes derniers remerciements vont à Susan à qui je dédie ce mémoire, car
c’est elle qui m’a donné l’envie de prendre cette voie dans laquelle je m’épanouis.
TABLE DES MATIERES
CHAPITRE 2 : RECHERCHE
D’ARCHITECTURES
ACRONYME....................................................................................................................................... 21
NOMENCLATURE ............................................................................................................................ 22
INDICES. ............................................................................................................................................. 22
2.1. APPROCHE TRADITIONNELLE (BOTTOM-UP) .............................................................. 23
2.2. APPROCHE HYBRIDE (TOP-DOWN/BOTTOM-UP) ........................................................ 23
2.3. CAS TEST : ACTIONNEUR DU PROJET DRESS ............................................................... 26
2.3.1. Fonctions et contraintes principales ....................................................................................... 26
2.3.2. Architectures fonctionnelles................................................................................................... 30
2.3.3. Architectures conceptuelles ................................................................................................... 33
2.3.4. Architectures organiques........................................................................................................ 36
2.4. ETUDE DE FIABILITE AXEE SECURITE........................................................................... 41
2.4.1. Exigences de fiabilité ............................................................................................................. 41
2.4.2. AMDEC et arbres de défaillance ........................................................................................... 42
2.4.3. Distribution exponentielle...................................................................................................... 43
2.4.4. Taux de défaillance ................................................................................................................ 45
2.4.5. Temps d’exposition à une défaillance.................................................................................... 45
2.4.6. Résultats intermédiaires ......................................................................................................... 46
2.4.7. Spécification des taux de défaillance ..................................................................................... 46
i
TABLE DES MATIERES
CHAPITRE 3 : DIMENSIONNEMENT EN
PUISSANCE
ACRONYME....................................................................................................................................... 58
NOMENCLATURE ............................................................................................................................ 58
INDICES .............................................................................................................................................. 59
EXPOSANT ......................................................................................................................................... 59
3.1. APPROCHE DE MODELISATION ........................................................................................ 61
3.1.1. Dimensionnement en puissance à l’aide de la simulation inverse ......................................... 61
3.1.2. Modélisation paramétrique..................................................................................................... 62
3.1.3. Choix de l’approche de modélisation..................................................................................... 65
3.2. LOIS D’ECHELLE .................................................................................................................... 66
3.2.1. Principe de base...................................................................................................................... 66
3.2.2. Exemple de composant électrique : les moteurs brushless..................................................... 68
3.2.3. Exemple de composant mécanique simple : le réducteur de vitesse...................................... 76
3.2.4. Exemple de composant mécanique complexe : les vis à billes ou à rouleaux ....................... 81
3.2.5. Validation............................................................................................................................... 85
3.2.6. Bilan et tableaux de synthèse ................................................................................................. 91
3.3. GRANDEURS DIMENSIONNANTES .................................................................................... 95
3.3.1. Grandeurs dimensionnantes ................................................................................................... 95
3.3.2. Couple nominal des réducteurs de vitesse.............................................................................. 95
3.3.3. Couple RMS et modèles thermiques des moteurs électriques ............................................... 96
3.4. IMPLEMENTATION DANS DYMOLA/MODELICA.......................................................... 99
3.5. CAS TEST : COMMANDE DE VOL PRIMAIRE (AILERON)......................................... 102
3.5.1. Modèle de charge et profil de mission ................................................................................. 103
3.5.2. Dimensionnement en puissance ........................................................................................... 104
3.5.3. Résultats ............................................................................................................................... 106
3.6. CAS TEST : ORIENTATION DE TRAIN D’ATTERRISSAGE ........................................ 108
3.6.1. Profil de mission .................................................................................................................. 108
3.6.2. Dimensionnement en puissance ........................................................................................... 109
3.6.3. Résultats ............................................................................................................................... 110
3.7. CONCLUSION ......................................................................................................................... 112
REFERENCES .................................................................................................................................. 114
ACRONYME..................................................................................................................................... 118
NOMENCLATURE .......................................................................................................................... 118
INDICES ............................................................................................................................................ 118
EXPOSANT ....................................................................................................................................... 119
ii
TABLE DES MATIERES
iii
TABLE DES MATIERES
CONCLUSION
CONCLUSION.................................................................................................................................. 191
REFERENCES .................................................................................................................................. 197
ANNEXES
ANNEXE A : STATISTIQUES POUR LA FIABILITE .................................................................III
ANNEXE B : MODES ET TAUX DE DEFAILLANCE CONSIDERES POUR L’ETUDE DE
FIABILITE DE L’EXEMPLE D’ACTIONNEUR DE DRESS. .............................XI
ANNEXE C : ARBRES DE DEFAILLANCE DE L’EXEMPLE DE L’ACTIONNEUR DE
DRESS .........................................................................................................................XV
ANNEXE D : MODELE DE LA VIS SANS FIN (EN ANGLAIS)……………….……………..XIX
iv
CHAPITRE 1 – INTRODUCTION GENERALE
CHAPITRE 1
INTRODUCTION GENERALE
ACRONYMES....................................................................................................................................... 1
NOMENCLATURE .............................................................................................................................. 1
1.1. LE PROJET DRESS .................................................................................................................... 3
1.2. LE PROJET C6E2........................................................................................................................ 4
1.3. SPECIFICITES DE CONCEPTION DES ACTIONNEURS A SOURCE DE PUISSANCE
ELECTRIQUE.............................................................................................................................. 5
1.4. METHODOLOGIE DE CONCEPTION CLASSIQUE ........................................................... 7
1.5. ETAT DE L’ART DES OUTILS LOGICIELS D’AIDE A LA CONCEPTION.................... 9
1.5.1. Conception orientée simulation................................................................................................ 9
1.5.2. Outils de modélisation et de simulation physique.................................................................. 11
1.5.3. Simulation directe et inverse.................................................................................................. 13
1.6. METHODOLOGIE DE CONCEPTION PROPOSEE ET PLAN DE THESE.................... 14
REFERENCES .................................................................................................................................... 18
ACRONYMES
NOMENCLATURE
1
CHAPITRE 1 – INTRODUCTION GENERALE
Les pistes envisagées pour répondre aux défis du transport aérien mentionnés
précédemment sont les suivantes :
Conception écologique prenant en compte toutes les phases de vie d’un avion
(production, exploitation et fin de vie ou recyclage).
2
CHAPITRE 1 – INTRODUCTION GENERALE
Les travaux présentés ici s’inscrivent dans le cadre général du développement d’avions
plus électriques. Le but de cette approche est de converger vers des avions « tout
électriques », dans lesquels la source de puissance secondaire serait exclusivement distribuée
par réseau électrique. Dans cette optique, un effort de recherche important a été engagé depuis
les années 90. Dans [Maré 2010] l’auteur rappelle ainsi, qu’au niveau européen, cet effort de
recherche s’est d’abord concentré sur les équipements de puissance (en particulier les
actionneurs de commandes de vol) avec, par exemple, les projets de recherche Epica (1993-
1997) et Elisa (1999-2001). Dans un deuxième temps, cet effort de recherche s’est porté sur
les systèmes de puissance (air conditionné, génération électrique, commande de vol et trains
d’atterrissage) avec le projet Power Optimized Aircraft (POA, 2001-2007), qui a impliqué 46
partenaires différents pour un montant total de l’ordre de 100 million d’euros. Par la suite,
l’effort de recherche a été porté au niveau aéronef (réseau de transmission de puissance et
intégration) avec le projet More Open Electrical Technologies (MOET, 2006-2009), qui a
impliqué 62 partenaires différents pour un financement total de l’ordre de 37 million d’euros.
L’activité de recherche et de développement autour de la thématique de l’avion plus électrique
apparaît aussi clairement à travers le nombre important d’ateliers et de conférences, qui y sont
dédiés (par exemple : [DéciElec 2009 ; MOET 2009 ; SEE 2009]). Au niveau industriel, et
dans le secteur du transport commercial, les derniers modèles A380 et B787 des avionneurs
de premier rang Airbus et Boeing intègrent déjà une technologie plus électrique que leurs
prédécesseurs. Cette tendance se retrouve aussi dans d’autres secteurs que l’aéronautique, et
notamment dans l’automobile [Larminie 2003].
Dans l’objectif de présenter de façon générale le contexte dans lequel les travaux
exposés dans ce rapport ont été conduits, ce chapitre commence par décrire dans les sections
1.1 et 1.2 les deux projets de recherche qui les ont financés. Puis les spécificités de la
conception des actionneurs électriques embarqués de forte puissance sont abordées dans la
section 1.3. La section 1.4 décrit la méthodologie de conception classique et la section 1.5
l’état de l’art des outils de modélisation pour la conception. Finalement, après la présentation
du contexte des travaux et des problématiques rencontrées, la section 1.6 présente la
méthodologie mise en place pour y répondre.
1
Site internet du projet DRESS : www.dress-project.eu
3
CHAPITRE 1 – INTRODUCTION GENERALE
Dispositif à Moteurs
source de électriques
puissance
hydraulique
Vis-sans-fin
Réducteur
Embrayage de vitesse
Actionneur de DRESS
Figure 1.1 : Projet DRESS, remplacement du dispositif à source de puissance hydraulique actuel
d’orientation de la direction de train avant par un dispositif à source de puissance électrique et plus
fiable (photographie du prototype fournie par Equip’Aéro).
Avec une fiabilité très élevée, le nouveau système doit permettre de diriger l’avion lors
d’atterrissages dans des conditions de visibilité réduite (par exemple : atterrissage CAT III C
[AC 120-28D]). De cette façon, le projet DRESS s’inscrit dans l’objectif global
d’automatisation des phases d’approche, d’atterrissage, de manœuvres au sol et de décollage.
A terme, une plus grande automatisation des opérations de l’avion au sol permettrait de
réduire les temps d’attente en l’air et l’engorgement au sol en cas de mauvaise visibilité.
4
CHAPITRE 1 – INTRODUCTION GENERALE
L’objectif du projet C6E2 est de compléter les outils de simulation pour l’ingénierie
du développement des systèmes électriques embarqués. Cet objectif est motivé par
l’importance croissante que prennent les systèmes électriques dans les applications
embarquées : avion tout électrique, véhicules hybrides, commandes électriques, systèmes
motorisés, etc. Pour maîtriser ces systèmes et optimiser les performances, il faut repenser le
processus de conception, pour aboutir à une approche intégrée. Le tout dans des délais de plus
en plus réduits et avec une préoccupation de valeur accrue pour le produit final. La simulation
joue un rôle clef dans cette perspective.
Partant de cette expression de besoin, le projet C6E2 se structure autour de trois axes,
ayant chacun leurs objectifs techniques propres :
D’une manière très générale, l’objectif de conception des actionneurs est de répondre à
un besoin de puissance exprimé au niveau de la charge à piloter (Figure 1.2). Ce besoin peut
prendre différentes formes selon le domaine d’application. Par exemple, en robotique ce
besoin prend souvent la forme d’un trièdre force ( F ), position ( x ), vitesse ( x& ). Cette
représentation est particulièrement utile dans ce cas où certaines caractéristiques de la charge
dépendent de sa position. Par exemple : l’inertie d’un bras de robot peut varier fortement en
fonction de sa position angulaire. D’un autre côté, dans le cas des applications considérées
dans ces travaux (orientation de train avant d’atterrissage, commandes de vol, etc.),
généralement le besoin en puissance est représenté dans le plan puissance force ( F ), vitesse
( x& ), qui est une représentation classique pour les applications industrielles.
5
CHAPITRE 1 – INTRODUCTION GENERALE
5.0E3
0.0E0
Force F (N)
-5.0E3
-1.0E4
-1.5E4
-2.0E4
-2.5E4
-3.0E4
0 1 2 3 4 5 6 7 8 9 10
Temps (s)
0.12
Position x (m) 0.10
0.08
0.06
0.04
0.02
0.00
-0.02
0 1 2 3 4 5 6 7 8 9 10
Temps (s)
Vitesse x (m/s)
0.2
0.1
.
0.0
-0.1
-0.2
0 1 2 3 4 5 6 7 8 9 10
Temps (s)
Figure 1.2 : Exemple de besoin de puissance exprimé sous la forme d’un profil de mission force,
position et vitesse.
Régime
SOA Régime
intermittent ou
permanent
Charge Safe Operating
Force F (N)
Force F (N)
Area (SOA)
0 0
Charge
-0.50 -0.25 0.00 0.25 0.50 -0.50 -0.25 0.00 0.25 0.50
0. 0.
Vitesse x (m/s)
Figure 1.3 : Exemples de besoin et de capacités de puissance exprimés dans le plan puissance force,
vitesse pour des dispositifs à source de puissance hydraulique (a) et électrique (b).
La figure 1.3-(a) montre que le besoin de puissance tracé dans le plan puissance force
( F ), vitesse ( x& ) est bien inclus dans l’enveloppe des capacités de l’actionneur hydraulique.
Cette enveloppe est représentative à la fois des capacités en fonctionnement intermittent que
permanent de l’actionneur. Par conséquent, l’actionneur est bien capable de répondre au
besoin.
D’un autre côté, comme le montre la figure 1.3-(b) la technologie électrique est
caractérisée par des capacités en régime intermittent différentes des capacités en régime
6
CHAPITRE 1 – INTRODUCTION GENERALE
permanent. Par exemple, si l’on considère un moteur électrique sans balais, son couple
maximal admissible est limité en fonctionnement intermittent par le risque de
démagnétisation des aimants permanents. Ce dernier couple est généralement bien plus élevé
que le couple maximal admissible en régime permanent, qui est imposé par la limite
d’endurance thermique du bobinage du moteur. Le couple maximal admissible en régime
permanent est aussi appelé couple RMS, en référence à la méthode de calcul de la moyenne
carrée (en Anglais : « Root Mean Square ») utilisée pour l’évaluer.
Nous avons pris la liberté ici de désigner les capacités de puissance en régime
intermittent par le terme Safe Operating Area (SOA). Typiquement utilisé en électronique, le
terme SOA définit les conditions d’opération limites à l’intérieur desquelles un composant
peut opérer pendant un temps limité sans risque de s’endommager.
La conception système ;
La conception spécifique ;
L’intégration système.
Le point de départ du cycle de conception est constitué des exigences du cahier des
charges, qui traduisent les attentes et les besoins que doit satisfaire le système à développer.
7
CHAPITRE 1 – INTRODUCTION GENERALE
Exigences Produit
me
Conc
ystè
eptio
s
Validation
ion
ns
g ra t
ystèm
Inté
e Conception spécifique
Ingénierie mécanique
Ingénierie électrique
Tech. de l’information
Modélisation et simulation
Entrée
Proto Proto Pdt …
…
Exigences de labo fonct° pilote
8
CHAPITRE 1 – INTRODUCTION GENERALE
C’est pourquoi, il est préférable de réaliser des boucles d’itération le plus tôt
possible dans le cycle de conception pour accélérer ce dernier.
La figure 1.4 met en évidence le support apporté par les outils logiciels de
modélisation et de simulation lors des différentes phases de conception décrites
précédemment. Ces outils logiciels sont généralement spécifiques à chaque discipline
technique.
9
CHAPITRE 1 – INTRODUCTION GENERALE
Exigences Produit
Conception de
Paramètres 3D Paramètres
détail :
localisés 3D localisés
CAO, élts finis
Modèles
géométriques
Pour un même niveau de détail, les différentes formes de modèles sont liées aux
différentes disciplines rencontrées (mécanique, électrique, électromagnétisme, etc.). Cette
particularité soulève le problème de l’intégration des différents types de modèles dans un
même environnement de simulation, pour pouvoir prendre en compte les interactions entre les
différentes disciplines techniques et appliquer une approche de conception multidisciplinaire.
10
CHAPITRE 1 – INTRODUCTION GENERALE
Les approches de modélisation graphiques les plus connues sont les graphes de liaison
(en Anglais : « Bond graph ») et les graphes linéaires. Les graphes de liaison sont basés sur la
modélisation multidisciplinaire structurée [Vergé 2004], qui en fait un outil très puissant et
aussi très répandu. Cependant, cette représentation est réservée à une description purement
comportementale et ne facilite pas une modélisation topologique intuitive (Figure 1.6). De
plus, ces graphes sont caractérisés par un formalisme particulier, difficilement absorbable par
le monde industriel, et supportent mal la modélisation de systèmes non linéaires [Nfonguem
2006]. D’un autre côté, les graphes linéaires sont eux aussi multidisciplinaires et reflètent
fidèlement la topologie du système modélisé. Le langage de modélisation VHDL-AMS a été
développé sur cette deuxième approche.
Les deux grands types de langage de modélisation sont les langages procéduraux
(séquentiels) et les langages déclaratifs. Avec un langage procédural, un modèle de simulation
temporel (dynamique) est défini par une procédure ou séquence prédéfinie de calculs des
dérivées des différentes variables d’état. Avec un langage déclaratif, le modèle est défini par
un jeu d’équations, qui établissent le lien entre les différentes variables d’états et leurs
11
CHAPITRE 1 – INTRODUCTION GENERALE
dérivées. Un compilateur convertit ces fonctions en une expression adaptée aux algorithmes
de simulation. Parmi les principaux langages déclaratifs, on distingue [Paredis 2001 ; Ferretti
2004]: Omola, Easy5, Modelica et VHDL-AMS.
Charge Frottement
inertielle visqueux
Moteur
électrique Légende
i
Ligne associée à la
L R C Ω transmission de
U puissance électrique
J
K Ligne associée à la
fv
transmission de
puissance mécanique
Bond Graph
1:L 1:J
U C
1:i GY:K 1:Ω
i Ω
R:R R:fv
Figure 1.6 : Modèle Bond Graph d’un exemple de dispositif électromécanique simple.
12
CHAPITRE 1 – INTRODUCTION GENERALE
Pour illustrer la différence entre la simulation inverse et directe, nous considérons ici
l’exemple d’une inertie (moment d’inertie J) en rotation pure. L’inertie se déplace sous l’effet
des moments extérieurs C1 et C2. Par conséquent, durant une simulation directe les variables
d’entrée sont les efforts extérieurs C1 et C2 et la variable de sortie est le déplacement angulaire
θ de l’inertie pure. A l’opposé, durant une simulation inverse, le déplacement angulaire est
une variable d’entrée avec l’un des deux efforts extérieurs, l’autre effort extérieur étant la
variable de sortie. Dans le cadre d’un modèle déclaratif ce dispositif peut être défini par la
seule équation suivante basée sur le principe fondamental de la dynamique (PFD) en
convention de signe récepteur :
d 2 θ(t )
J = C1 (t ) + C 2 (t ) ( 1.1 )
dt 2
θ(t ) = ∫∫ (C (t ) + C (t ))dt
1
1 2 ( 1.2 )
J
13
CHAPITRE 1 – INTRODUCTION GENERALE
d 2 θ(t )
C1 (t ) = J − C 2 (t ) ( 1.3 )
dt 2
14
CHAPITRE 1 – INTRODUCTION GENERALE
(a) (b)
Processus classique de Produit Processus proposé de Produit
Exigences développement virtuel virtuel Exigences développement virtuel virtuel
Pour atteindre une juste pré-spécification des composants dès la phase de conception
système, la méthodologie développée, dans le cadre du présent travail de thèse, se décompose
de la façon séquentielle suivante (Figure 1.8) :
Recherche d’architectures ;
Analyse modale.
Les boucles d’itération implémentées dans la phase de conception système sont liées
aux cas test utilisés pour illustrer la méthodologie et les outils mis en place dans le cadre des
présents travaux. Un accent est mis sur les critères de fiabilité et de masse dans le cadre de
l’actionneur du projet DRESS. A ces deux aspects, s’ajoute l’intégration géométrique, dans le
cas de l’actionneur des commandes de vol électriques. Finalement, l’actionneur de poussée
vectorielle du premier étage du lanceur spatial VEGA [Carnevale 2007 ; ESA 2009] est
utilisé, pour illustrer la prise en compte des limitations des performances en boucle fermée.
Le chapitre 1 présente de façon générale le contexte dans lequel les travaux présentés
ont été conduits. Les problématiques abordées et les principaux axes suivis pour y répondre y
sont aussi décrits.
15
CHAPITRE 1 – INTRODUCTION GENERALE
Le chapitre 5 traite de la prise en compte des modes naturels dans l’évaluation des
performances en boucle fermée. Des règles de conception de haut niveau sont proposées ici
pour permettre d’évaluer les performances de rapidité en boucle fermée sans développer de
commande. Cet aspect est illustré avec un cas test caractérisé par des performances
dynamiques importantes : l’actionneur de poussée vectorielle du premier étage du lanceur
spatial VEGA.
16
CHAPITRE 1 – INTRODUCTION GENERALE
Cas test
Chapitre 2 : DRESS
Recherche
d’architectures, Recherche
Recherche
Etude de fiabilité d’architectures
d’architectures Etude de
architecturale. fiabilité
Architectures architecturale
Ph candidates
ase
Chapitre 3 : de Bilan de DRESS
Dimension
Dimension
c
Dimensionnement en masse Commandes
onc
-nement
-nement
puissance. Intégration de vol
ept
en
en puissance
puissance
Chapitre 4 : géométrique
ion
Composants Phase de
pré-spécifiés conception
spécifique
17
CHAPITRE 1 – INTRODUCTION GENERALE
REFERENCES
[AC 120-28D] U.S. Departement of Transportation Federal Aviation Administration. Criteria for
Approval of Category III Weather Minima for Takeoff, Landing, and Rollout. Advisory
Circular. AC 120-28D, 1999.
[Carnevale 2007] C. Carnevale and P. D. Resta. Vega Electromechanical Thrust Vector Control
Development. 43rd AIAA/ASME/SAE/ASEE Joint Propulsion Conference & Exhibit
Cincinnati, OH, USA, 8-11 July 2007.
[ESA 2009] European Space Agency. "European Space Agency Website - Launch Vehicles." from
http://www.esa.int/, 2009.
[Ford 2005] T. Ford. "More-electric aircraft." Aircraft Engineering and Aerospace Technology 77(1),
2005.
[Jasnocha 2001] U. Jasnocha et S. Haasb. "A collaborative environment based on distributed object-
oriented databases " Computers in Industry 29(1-2): 51-61 2001.
[Larminie 2003] J. Larminie and J. Lowry. Electric vehicle technology explained,Wiley, 2003.
[Liscouet - Hanke 2008] S. Liscouet - Hanke. A model-based methodology for integrated preliminary
sizing and analysis of aircraft power system architectures. Laboratoire de Génie Mécanique de
Toulouse (LGMT) INSA-UPS. Toulouse, Université de Toulouse. Ph. D., 2008.
[Maré 2010] J.-C. Maré. Towards more electric drives for embded applications: (re)discovering the
advantages of hydraulics. 7th International Conference on Fluid Power (IFK) Aachen,
Germany, 22 - 24 March 2010.
[MOET 2009] MOET Project Consortium. More Electric Aircraft Forum, Barcelona, Spain, 8-11
September 2009.
[Nangia 2007] R. Nangia. A vision for highly Fuel-efficient Commercial Aviation. EUCASS 2007:
2nd European conference for Aerospace Sciences. Brussels, Belgium, July 1-6, 2007.
[Otter 1999] M. Otter, H. Elmqvist et S. E. Mattsson. Hybrid Modeling in Modelica based on the
Synchronous Data Flow Principle. CACSD’99. Hawaii, USA, August 22-26, 1999.
18
CHAPITRE 1 – INTRODUCTION GENERALE
[Pahl 2004] G. Pahl, W. Beitz, J. Feldhusen and K.H. Grote. Engineering Design: A Systematic
Approach,L. Blessing, London, Springer, 2004.
[Pelegrin 2008] M. Pelegrin. La sécurité des vols. Les entretiens de Toulouse, rencontre aérospatiale.
Toulouse, France, Académie de l'aire et de l'espace, Collège de Polytechnique, 2008.
[SEE 2009] SEE, Aerospace Vallée et 3AF. Vers des aéronefs encore plus électriques, Toulouse, 28 et
29 janvier 2009.
[VDI 2206] VDI Richtlinien. Design methodology for mechatronic systems VDI-Richtlinien.
Düsseldorf, Beuth Verlag, 2004.
[Vergé 2004] M. Vergé et D. Jaume. Modélisation structurée des systèmes avec les Bond
Graphs.Automatique,Technip, Paris, 2004.
19
CHAPITRE 1 – INTRODUCTION GENERALE
20
CHAPITRE 2 – RECHERCHE D’ARCHITECTURES
CHAPITRE 2
RECHERCHE
D’ARCHITECTURES
ACRONYMES..................................................................................................................................... 21
NOMENCLATURE ............................................................................................................................ 22
INDICES .............................................................................................................................................. 22
2.1. APPROCHE TRADITIONNELLE (BOTTOM-UP) .............................................................. 23
2.2. APPROCHE HYBRIDE (TOP-DOWN/BOTTOM-UP) ........................................................ 23
2.3. CAS TEST : ACTIONNEUR DU PROJET DRESS ............................................................... 26
2.3.1. Fonctions et contraintes principales ....................................................................................... 26
2.3.2. Architectures fonctionnelles................................................................................................... 30
2.3.3. Architectures conceptuelles ................................................................................................... 33
2.3.4. Architectures organiques........................................................................................................ 36
2.4. ETUDE DE FIABILITE AXEE SECURITE........................................................................... 41
2.4.1. Exigences de fiabilité ............................................................................................................. 41
2.4.2. AMDEC et arbres de défaillance ........................................................................................... 42
2.4.3. Distribution exponentielle...................................................................................................... 43
2.4.4. Taux de défaillance ................................................................................................................ 45
2.4.5. Temps d’exposition à une défaillance.................................................................................... 45
2.4.6. Résultats intermédiaires ......................................................................................................... 46
2.4.7. Spécification des taux de défaillance ..................................................................................... 46
2.5. CONCLUSION ........................................................................................................................... 53
REFERENCES .................................................................................................................................... 54
ACRONYMES
21
CHAPITRE 2 – RECHERCHE D’ARCHITECTURES
NOMENCLATURE
INDICES
A, B, C Evénement A, B ou C
c Conceptuel
D Mode d’opération dégradé
f Fonctionnel
G, g Grippage
L Dirigeabilité libre et
incontrôlée du train avant
d’atterrissage
o Organique
22
CHAPITRE 2 – RECHERCHE D’ARCHITECTURES
Le but de cette phase de conception est d’identifier, en nombre réduit, les architectures
potentiellement les plus intéressantes vis-à-vis des fonctions que le système à développer doit
assurer et des contraintes qu’il doit respecter.
Les travaux présentés dans cette section ont été conduits dans le cadre du projet
DRESS, qui a pour but de développer un dispositif d’actionnement à source de puissance
électrique pour l’orientation du train avant d’un avion monocouloir. Plus exactement, ces
travaux concernent le développement de la chaîne de transmission de puissance
électromécanique du dispositif. Les nouvelles exigences fonctionnelles et les contraintes de
sécurité et de fiabilité sont très sévères et ne permettent pas d’appliquer une approche de
génération d’architectures traditionnelle. Pour répondre au mieu aux besoins de ce projet, une
approche innovatrice a été développée.
23
CHAPITRE 2 – RECHERCHE D’ARCHITECTURES
Sous-
fonction
Figure 2.1 : Décomposition de la fonction principale d’un système en sous-fonctions, extrait de [Pahl
2004].
Le principal problème soulevé par une approche descendante est le nombre total de
solutions identifiées, qui croît exponentiellement avec le niveau de détail. Afin d’éviter
l’explosion du nombre de solutions, l’approche développée ici combine la génération
descendante de solutions avec un filtrage à chaque niveau de détail, pour obtenir un nombre
réduit d’architectures pertinentes en fin de processus. Le filtrage est réalisé sur la base de
critères qualitatifs issus de l’état de l’art technologique et des contraintes de conception (par
exemple : prohibition de l’utilisation d’éléments pyrotechniques). Les contraintes de
conception et les limites technologiques actuelles sont traduites à chaque niveau de détail de
la conception, pour pouvoir éliminer toutes les solutions irréalisables ou incompatibles avec
l’application, au fur et à mesure de l’avancement du processus. Le filtrage des architectures
est donc un processus montant. Pour cette raison, l’approche globale présentée est dite
hybride, car à la fois descendante pour la génération d’architectures et montante pour la prise
en compte de la réalité technologique. La figure 2.2 illustre cette approche globale de
génération d’architectures hybride décomposée en trois niveaux de détail distincts :
fonctionnel, conceptuel et organique.
La figure 2.3 illustre l’approche proposée ici et les différents niveaux de détail
considérés avec un exemple d’architecture d’actionneur électromécanique issu du cas test
traité dans la section suivante.
24
CHAPITRE 2 – RECHERCHE D’ARCHITECTURES
Nombre de Analyse
descendante
Analyse
Approche
solutions fonctionnelle
fonctionnelle
Architectures
Architectures
7860 fonctionnelles
fonctionnelles
Niveau Filtre
fonctionel
Architectures
Architectures
120 fonctionnelles
fonctionnelles filtrées
filtrées
Architectures
Architectures
1920 conceptuelles
conceptuelles
Niveau Filtre
conceptuel
Architectures
Architectures
16 conceptuelles
conceptuelles filtrées
filtrées
Approche
montante
Architectures
Architectures
80 technologiques
technologiques
Niveau Filtre
organique
(technologique)
3 Architectures
Architectures Technologie
Technologie
technologiques
technologiques filtrées
filtrées
Transformer puissance
transform electrical Moteur Moteur
Source de puissance
mechanical power Moteur Moteur
électrique en puissance
power into mechanical brushless brushless
mécanique
input elec. rotatif elec. rotatif
mécanique cylindrique cylindrique
Transformer puissance
adapt mechanical Etage
mech. Etage
mech. Transform Transform Réducteur Réducteur
RR power RR power
mécanique
power 1 méca.
stage 21 méca.
stage 21 ‘ -ation RR -ation RR cycloïdal cycloïdal
adaptation adaptation
Isolertowing
défaillance Embrayage Embrayage
ensure modede& Etage
mech. Etage
mech. Embrayage
clutch Embrayage
clutch électro- électro-
grippage en rotation
jam tolerance méca.
stage 32 méca.
stage 32 ‘ mécanique
clutch mécanique
clutch
Combiner voies de
sum mechanical Interface
interfacede l’actionneur
with landing Sommation
torque de couple
summation on
puissance Vis-sans-fin Vis-sans-fin
powers avec legear
tubeleg
tournant sur le tube
NLGLtournant
Figure 2.3 : Exemple d’une architecture aux différentes étapes (fonctionnelle, conceptuelle et
organique) de la méthodologie de génération proposée.
25
CHAPITRE 2 – RECHERCHE D’ARCHITECTURES
Comme le montrent la figure 2.2 et la figure 2.3, entre les niveaux fonctionnels et
conceptuels, un cadre architectural est défini. Ce cadre est une première définition physique
des architectures et sert de lien entre la description fonctionnelle et la description conceptuelle
des solutions. Ce lien est nécessaire, car, contrairement à la description fonctionnelle, la
description conceptuelle est basée sur une combinaison de solutions (conceptuelles) données.
Le cadre fonctionnel permet ainsi de procéder à la distribution et à la permutation des
solutions conceptuelles, pour engendrer les différentes combinaisons possibles.
26
CHAPITRE 2 – RECHERCHE D’ARCHITECTURES
analyses peuvent être conduites efficacement avec des approches classiques [Tassinari 1995 ;
ARP 4761]. Le but des travaux présentés ici n’est pas de détailler ces analyses. Cependant, il
est intéressant de situer le contexte de développement de l’actionneur DRESS, pour mieux
appréhender les tenants et les aboutissants des fonctions, des contraintes et des exigences
principales, qui guident le processus de génération d’architectures.
Roulements
Tube
tournant
Fluide
Jambe du d’amortissement
train
Essieu
Roue
27
CHAPITRE 2 – RECHERCHE D’ARCHITECTURES
Le premier mode d’opération est assuré par les deux fonctions principales classiques
mentionnées précédemment. Le deuxième mode d’opération impose une fonctionnalité
supplémentaire à l’actionneur :
Figure 2.5 : Atterrissage d’urgence du vol JetBlue Airways Flight 292. Photographies extraites de
[Scott 2005].
La figure 2.5 illustre l’atterrissage d’urgence du vol JetBlue Airways Flight 292. Le 21
septembre 2005, ce vol de transport commercial en provenance de Burbank (CA, É.-U.) et à
destination de la ville de New-York (NY, É.-U.) a dû être détourné sur Los-Angeles pour un
atterrissage d’urgence avec une orientation de la direction du train d’atterrissage avant
bloquée à 90° (perpendiculaire à la direction de vol). Le blocage de la direction de
l’orientation du train a résulté de la séquence d’événements suivants [NTSB 2008]:
Le pilote procède à une sortie des trains d’atterrissage. Durant cette procédure, le train
avant se déploie 1,5 seconde avant les trains centraux et la fermeture des portes.
28
CHAPITRE 2 – RECHERCHE D’ARCHITECTURES
Après une tentative de recentrage sur une durée de 0,5 s, le système de contrôle
considère qu’il y a défaillance du système d’orientation de la direction du train avant
et le neutralise. Le train avant se retrouve ainsi bloqué avec une direction bloquée à
90°.
Malgré le risque encouru durant cette procédure d’urgence, l’avion a pu être posé sans
qu’aucun des 140 passagers ne soit blessés. Après cet incident, l’avionneur a édité une
procédure permettant de réinitialiser en vol le système de contrôle de la direction de
l’orientation du train avant d’atterrissage. La procédure de test de l’amortisseur du train a été
modifiée pour réduire la fatigue mécaniques induite au niveau des prises de la came de
centrage. Finalement, la prise supérieure et le tube plongeant de la came de centrage ont été
reconçus pour être plus robustes.
Tube
plongeant
Tube
glissant
Prise
supérieure
Prise
inférieure
Came de centrage
Figure 2.6 : Came de centrage de train avant d’avion de transport monocouloir en position désengagée
(gauche) et engagée (droite).
29
CHAPITRE 2 – RECHERCHE D’ARCHITECTURES
est requis de pouvoir basculer vers une orientation du train d’atterrissage en mouvement libre,
afin de permettre au pilote d’atterrir et de diriger l’avion au sol à l’aide d’un moyen alternatif
comme par exemple le freinage différentiel. En effet, bien que ce dernier mode de
fonctionnement soit hasardeux, il est aussi le mode de fonctionnement ultime en cas de
défaillances multiples (en Anglais : « ultimate reversionary mode »). En d’autres termes,
l’actionneur doit inclure la fonctionnalité suivante :
Pour des raisons de sécurité l’actionneur peut être redondant et constitué de plusieurs
voies de puissance.
Active – active. En condition normale d’opération, les deux voies sont actives. En cas
de défaillance d’une voie, celle-ci est désactivée et la voie saine assure seule les
fonctions de l’actionneur.
30
CHAPITRE 2 – RECHERCHE D’ARCHITECTURES
Les deux voies peuvent être dimensionnées pour fournir ensemble la puissance requise
par l’actionneur. En cas de défaillance d’une voie, l’actionneur reste fonctionnel (en
Anglais : « fail functional »), mais ses performances sont dégradées (par exemple :
effort maximum transmis divisé par deux).
Une première voie, dite principale, est dimensionnée pour être capable de fournir seule
toute la puissance de l’actionneur. La deuxième voie, secondaire, est sous-
dimensionnée et ne peut seule fournir qu’une partie de la puissance requise par
l’actionneur. Cette dernière approche, caractéristique de la configuration active-
passive, permet de limiter la masse de l’actionneur, au prix de performances dégradées
en cas de défaillance de la voie de puissance principale.
31
CHAPITRE 2 – RECHERCHE D’ARCHITECTURES
(a) (b)
Avant filtrage Après filtrage Légende
Transformer
puissance élec.
E0 E0 E0 en puissance Ei Etage mécanique
méca. E0 E0
Transmission de
puissance mécanique
E1 E1 E1 Transformer E1 E1
puissance Combiner ou séparer
méca. voie(s) de puissance
E2 E2 E2 E2 E2 mécanique(s)
E3 E3 E3 E3 E3
Isoler défaillance
E4 E4 E4 E4 E4
Permettre la
E5 E5 E5 réversibilité E5 E5
Tube
tournant
Figure 2.7 : Architecture fonctionnelle globale avant (a) et après (b) filtrage par rapport aux exigences
de sécurité.
(a) (b)
(c) (d)
32
CHAPITRE 2 – RECHERCHE D’ARCHITECTURES
Le nombre total nf1 d’architectures différentes est donné par le calcul combinatoire
suivant basé sur l’architecture fonctionnelle globale de la figure 2.7-a :
Permutatio n une ou deux sources
des fonctions de puisance mécanique
} }
nf1 = 5! 2{5 2 = 7860 ( 2.1 )
Présence ou non
d'un chemin unique
entre chaque étage
Par conséquent, seules les solutions avec moins deux voies de puissances
totalement redondantes sont acceptables du point de vue des exigences de sécurité.
33
CHAPITRE 2 – RECHERCHE D’ARCHITECTURES
en puissance a montré que pour la présente application, un moteur linéaire serait trop lourd ou
trop encombrant.
Par conséquent seuls les moteurs électriques rotatifs sont considérés pour la suite
de l’étude.
Sachant que seuls les moteurs rotatifs sont considérés, uniquement deux familles
de transformations mécaniques sont possibles : rotation-rotation-rotation (RRR) et
rotation-translation-rotation (RTR).
Les voies de puissance mécanique redondantes peuvent être combinées par une
sommation de vitesse ou d’effort. En cas de grippage d’une des voies de puissance, la
sommation de vitesse permet d’assurer l’orientation du tube tournant, alors que la sommation
d’efforts entraîne un grippage. D’un autre côté, dans le cas où le chemin de puissance d’une
des voies est coupé (par exemple : rupture d’un arbre de transmission), la sommation d’effort
permet de continuer à diriger la roue avec la voie saine, alors que la sommation de vitesse
cause une orientation libre et incontrôlée de la direction de la roue. C’est pourquoi, les
éléments de sécurité assurant la fonction « Empêcher la propagation de défaillances jusqu’à
l’orientation de la direction de la roue», sont respectivement basés sur le principe de freinage
ou de débrayage pour les sommations de vitesse et d’effort. La sommation de vitesse est
technologiquement réalisable avec un mécanisme à train épicycloïdal, qui serait monté
directement sur le tube tournant. Cependant, cette solution implique un grand nombre
d’éléments rotatifs en contact avec ce tube, ce qui est inacceptable du point de vue du risque
de grippage.
L’état de l’art technologique des principes de sommation de vitesse fait que cette
solution n’est pas acceptable vis-à-vis du risque de grippage. La combinaison des voies
de puissance redondantes est donc réalisée par sommation d’effort et la fonction
« Empêcher la propagation de défaillances jusqu’à l’orientation de la direction de la
roue » est assurée par le débrayage de la voie défaillante.
1
Une défaillance dormante est une défaillance effective, qui pas encore été détectée.
34
CHAPITRE 2 – RECHERCHE D’ARCHITECTURES
E1 E1
Transformer • Transmission méca. :
puissance méca. RRR
E2 E2 TR-RR-RR RR-RR-RR
TR-RT-TR RR-RT-TR
TT-TR-RR RT-TR-RR
E3 E3 Isoler TT-TT-TR RTR
RT-TT-TR
défaillance
• Débrayage
E4 E4
Permettre la • Freinage
réversibilité
E5 E5
• Sommation d’efforts
Combiner voies de
puissance méca. • Sommation de vitesses
35
CHAPITRE 2 – RECHERCHE D’ARCHITECTURES
Le moteur à induction est basé sur un assemblage simple, robuste et bon marché. De
plus, avec un rotor à cage, ce moteur est compact, notamment pour les applications à hautes
vitesses. Le moteur CC est simple à commander et minimise les coûts associés à la
commande. Cependant, ces avantages sont relativisés par la présence de balais, qui requiert
une protection contre un environnement extérieur agressif, comme c’est le cas de DRESS. De
plus, l’usure rapide des balais implique de fréquents remplacements. Les moteurs brushless
ont un rendement élevé et une forte densité de puissance. Ces moteurs nécessitent l’utilisation
d’une commande significativement plus compliquée et coûteuse que pour les moteurs CC.
D’un autre côté, l’absence de balais permet de réduire significativement les besoins en
maintenance et réduit le risque de défaillance (par exemple : arcs électriques), notamment en
opérant dans un environnement agressif. Pour ces dernières raisons, le moteur brushless est le
standard actuel pour les applications aéronautiques. La technologie brushless présente de
nombreuses qualités en adéquation avec le projet DRESS. De plus, le partenaire responsable
de la production du moteur [HDD 2009] est spécialisé dans cette technologie. C’est pourquoi
la technologie brushless a été sélectionnée ici. Parmi les moteurs brushless émergent deux
catégories distinctes : les moteurs cylindriques, favorisant la vitesse, et les moteurs annulaires,
favorisant le couple (Figure 2.10). On note aussi la possibilité d’avoir, en alternative à
l’architecture cylindrique (ou à champ radial) classique, une architecture moteur de type
discoïde (ou à champ axial) qui présentent des avantages de compacité dans certaines
configurations. Cependant, le partenaire du projet DRESS responsable de la production du
moteur électrique est spécialisé dans les moteurs brushless cylindriques et ce sont ces derniers
qui par conséquent ont été sélectionnés ici.
36
CHAPITRE 2 – RECHERCHE D’ARCHITECTURES
Figure 2.10 : Moteurs brushless cylindriques (gauche) [HDD 2009], et annulaires (droite) [ETEL
2007].
37
CHAPITRE 2 – RECHERCHE D’ARCHITECTURES
raison d’un rendement réduit de l’Harmonic Drive pour une température d’opération
extrêmement basse, comme c’est le cas dans DRESS.
Figure 2.11: Réducteur cylcoïdal [Sumitomo 2005] (a), réducteur Harmonic Drive [Harmonic Drive
2005] (b) et réducteur épicycloïdal [APEX Dynamics 2009] (c).
38
CHAPITRE 2 – RECHERCHE D’ARCHITECTURES
Pour cette raison, dans les solutions RRR, une attention particulière doit être
accordée à l’aspect sécurité et fiabilité pendant la conception des composants en contact
avec le tube tournant de la jambe du train.
Ligne primitive de la
Rouleau crémaillère et axe de
Vis à
rouleaux rotation de la vis
Actionneur
installé sur la
jambe du train
Cercle avant
primitif de la EMA
roue du tube
tournant
Actionneur
voie de de
puissance 1
Roue du tube
tournant
Actionneur
voie de de
puissance 2
Vue de dessus Vue de dessus
Actionneur embrayé Actionneur débrayé
Figure 2.12 : Exemple de système de débrayage de l’interface de l’actionneur d’avec le tube tournant,
avec une architecture RTR constituée d’une crémaillère, d’une vis-à-rouleaux et d’un moteur brushless
cylindrique (images fournies par SPAB, SAAB et Messier-Dowty).
39
CHAPITRE 2 – RECHERCHE D’ARCHITECTURES
RRR (E1)
RTR
Réducteur Réducteurs Réducteur Direct-
épicyclique Harmo. Dr. Cycloïdal drive
40
CHAPITRE 2 – RECHERCHE D’ARCHITECTURES
Après filtrage, ce nombre de solutions est réduit à trois, dont deux sont caractérisées
par une transmission RRR et une par une transmission RTR. La figure 2.15 représente
graphiquement les trois grands types d’architectures organiques obtenues à la fin du
processus.
Figure 2.15 : Représentations graphiques des trois grands types d’architectures organiques.
41
CHAPITRE 2 – RECHERCHE D’ARCHITECTURES
L’AMDEC est une méthode qui vise à identifier les modes de défaillance des éléments
du système pour aboutir à leurs effets, et ainsi identifier les différents événements à risques
que peut rencontrer le système. Comme les événements à risques du système sont déduits des
modes de défaillance de ses organes, cette méthode est dite déductive ou montante.
L’AMDEC seule se prête difficilement à l’étude des effets de la combinaison de défaillances.
Pour cette raison, cette méthode est généralement accompagnée d’arbres de défaillance.
2
Failure Mode, Effect and Criticity Analysis.
3
Fault Tree Analysis.
42
CHAPITRE 2 – RECHERCHE D’ARCHITECTURES
Evénement
étudié
A
Porte OU
Evénement
« Evénement A étudié à
intermédiaire
lieu si défaillance C à lieu
OU si défaillances D ET E
B
ont lieu simultanément. »
C
Porte ET
Défaillances
élémentaires D E
Les arbres de défaillance mettent en évidence les différents scénarii pouvant conduire
à l’événement redouté. Ces arbres permettent aussi de procéder efficacement aux calculs de
fiabilité. Une porte ET signifie que l’effet (événement de niveau supérieur) résulte de
l’apparition simultanée des causes (évènements de niveau inférieur). D’après la théorie des
ensembles la probabilité que deux évènements indépendants A et B surviennent
simultanément pour causer un évènement C s’exprime par :
FC = A∩B = FA FB ( 2.6 )
De la même façon, une porte OU signifie que l’effet résulte de l’apparition d’une des
causes. D’après la théorie des ensembles la probabilité que deux évènements indépendants A
ou B surviennent pour causer un évènement C est donnée par :
FC = A∪B = FA + FB − FA FB ( 2.7 )
43
CHAPITRE 2 – RECHERCHE D’ARCHITECTURES
Mauvaise conception ;
Mauvaise installation ;
Mauvaise utilisation ;
Mauvaise maintenance ;
Défaut de fabrication ;
Les processus de défaillance liés aux origines précédentes sont multiples et souvent
complexes, voire imprévisibles (par exemple : présence d’une poussière dans l’engrenage
d’un réducteur de vitesse lors du processus de fabrication). En conséquence, quand on prend
en considération l’ensemble des causes de défaillance, il est très difficile de prédire
l’évolution du mécanisme général de défaillance et celles ci semblent survenir de façon
aléatoire ou catalectique [Schwob 1969 ; Dhillon 1988]. Le taux de défaillance correspondant
est donc constant.
F (t ) = 1 − e − λt (2.8 )
F (t ) = λt ( 2.9 )
La forme simplifiée précédente permet notamment de faciliter les calculs liés aux
arbres de défaillance décrits dans la section précédente. En appliquant la distribution
exponentielle aux calculs de probabilité avec une porte OU (cf. section 2.4.2) on observe le
résultat intéressant suivant :
λC
647 4
8
FC = A∪B (t ) = 1 − e −(λ A + λ B )t ≈ (λ A + λ B )t ( 2.10 )
44
CHAPITRE 2 – RECHERCHE D’ARCHITECTURES
Les temps d’opération de chaque composant ont été obtenus à partir d’un profil de
mission classique pour le type d’avion concerné par le projet DRESS (Figure 2.17),
accompagné d’une analyse des modes opératoires de l’actionneur.
Phases Remor- Roulage Décollage Train Train Train Atterriss- Roulage Remorqu-
quage avant sorti avant rentré avant sorti age age
Activité du
dispositif OFF ON ON ON OFF ON ON ON OFF
d‘orientation
de la roue
Temps
1/3 HO 3/400 HO 1/500 HO 1 HV 1/10 HO 6/500 HO 1/3 HO
(heure)
Cycle de vol
Figure 2.17: Description des phases d’une mission type de l’avion considéré pour le développement de
DRESS.
45
CHAPITRE 2 – RECHERCHE D’ARCHITECTURES
tous ses composants sont considérés comme n’étant exposés à des défaillances que durant la
phase de vol. Finalement, les composants situés en amont de l’embrayage (par exemple :
réducteur de vitesse, moteurs électriques, électroniques de puissance, commande) sont
supposés exposés uniquement pendant les phases de roulage, de décollage et d’atterrissage. A
partir des hypothèses précédentes, on obtient un temps d’exposition équivalent à 25% de
l’heure de vol pour les composants situés en aval de l’embrayage (inclus) et à 22% pour les
composants situés en amont.
Ces pourcentages d’activité permettent de faire le lien entre les exigences de fiabilité
exprimées en heure de vol (HV) et les taux de défaillance des composants, qui sont exprimés
en heure d’opération (HO). De cette façon, on passe d’un taux de défaillance exprimé en /HV
à un taux de défaillance exprimé en /HO à l’aide des coefficients a25 = 0,25 et a22 = 0,22 pour
les composants situés respectivement en aval et en amont de l’embrayage.
Tableau 2.2 : Taux d’apparition des modes de défaillance pour les trois architectures considérées.
Comme le montrent les résultats obtenus précédemment, avec des valeurs numériques
de taux de défaillances issues d’une base de données standard, les objectifs de fiabilité du
projet (10-9 /HV) ne sont atteints par aucune solution architecturale proposée. C’est pourquoi,
dans la sous-section suivante, et dans une approche de spécification des composants, nous
allons chercher à obtenir les valeurs numériques de taux de défaillance qui permettraient
d’atteindre ces objectifs.
46
CHAPITRE 2 – RECHERCHE D’ARCHITECTURES
A A
B B
C C
D E D E
Figure 2.18 : Approche traditionnelle montante (gauche) de calcul la probabilité d’apparition d’un
événement (A) à l’aide d’un arbre de défaillance, et approche montante et descendante (droite) pour la
spécification de taux de défaillance de composants (C).
L’objectif est de converger vers l’expression des taux de défaillance à spécifier pour
les composants interfaçant avec le tube tournant de la jambe du train, ainsi que les
composants situés en amont de la chaîne de transmission électromécanique de l’actionneur
(par exemple : électronique de puissance, commande). Cette méthode, qui consiste à spécifier
d’une manière descendante ou inductive les objectifs de fiabilité est en phase avec une
approche de conception plus systématique.
Pour appliquer cette méthode, il faut préalablement identifier les taux de défaillance
des composants autres que ceux spécifiés. Comme ces derniers composants ont une influence
limitée sur la fiabilité de l’actionneur, il convient de supposer les taux de défaillance
correspondant à une conception et une vie opérationnelle classiques, renseignés par une base
de données électronique [RIAC 2006] ou par les partenaires industriels du projet DRESS.
47
CHAPITRE 2 – RECHERCHE D’ARCHITECTURES
F0
Grippage en rotation
du tube tournant
F1.1 F2.1
Grippage voie Grippage voie
de puissance 1 de puissance 2
F1.21 F2.21
Grippage non F1.22 Grippage non F2.22
isolé en amont Grippage d’un isolé en amont Grippage d’un
de l’embrayage composant en de l’embrayage composant en
contact avec le contact avec le
tube tournant tube tournant
F1.43 F2.43
F1.41 F1.42 Grippage causé F2.41 F2.42 Grippage causé
Grippage du Grippage du par une défaillance Grippage du Grippage du par une défaillance
réducteur moteur réducteur moteur en amont du
en amont du
moteur moteur
La probabilité de grippage FG1.22 d’un composant en contact avec le tube tournant est
exprimée par heure de vol par :
où λ1.22 est le taux de défaillance (constant) par grippage d’un composant en contact avec le
tube tournant.
48
CHAPITRE 2 – RECHERCHE D’ARCHITECTURES
entraîne une révision de l’actionneur à la fin du vol. FG1.21 est la probabilité d’avoir
simultanément un composant défaillant situé en amont de l’embrayage susceptible de causer
un grippage de la voie de puissance et une défaillance de l’embrayage en position embrayée.
Cette probabilité de défaillance peut être calculée en faisant la somme des probabilités
d’apparition des différents scénarii correspondants [Schwob 1969] :
Probabilit é de
défaillance à l'instant t
T 678
FG1.21 (t ) = 2 ∫ f1.31 (t ) F 32 (t à t + 1) dt ≈ 2a 25 λ1.31λ1.32T ( 2.12 )
11.42
4 43 4
0 Probabilit é de défaillanc e
entre les instants t et t +1
où λG1.41 est le taux de défaillance par grippage du réducteur, λG1.42 est le taux de défaillance
par grippage du moteur et λG1.43 est le taux de défaillance des composants situés en amont du
moteur causant un grippage de la voie de puissance.
49
CHAPITRE 2 – RECHERCHE D’ARCHITECTURES
Les calculs font clairement apparaître un taux de défaillance maximal admissible des
composants situés en amont du moteur pouvant aller jusqu’à 0,9 /HO. Ce dernier résultat
confirme que grâce à la présence d’un embrayage, le taux de fiabilité des composants situés
en amont du moteur n’influe quasiment pas l’objectif de défaillance par grippage au niveau de
la jambe. D’un autre côté, pour des composants situés en amont du moteur et supposés
parfaits (probabilité nulle de grippage de la voie de puissance), au mieux le taux de
défaillance maximal admissible des composants d’interface est 2 10-9 /HO, ce qui équivaut à
des composants qui ne sont par nature pas sujets au grippage.
T
FL 0 (t ) = ∫ f L1.1 (t )FL1.1 (t à t + 1)dt ≈ λ2L1.1T ( 2.15 )
0
où f1.1(t) est la densité de probabilité (fonction de masse) d’avoir une perte de la transmission
de puissance dans une première voie de l’actionneur, F1.1(t à t+1) est la probabilité de perte de
la transmission de puissance de la deuxième voie pendant l’heure de vol suivant la défaillance
précédente et λL1. 1 est le taux de défaillance (constant) par perte de la transmission de
puissance d’une voie
avec
( ) (
λ L1.1 = a 25 λ L1.21 + λ L1.22 + a22 λ g1.31 + λ g1.32 + λ g1.33 ) ( 2.16 )
où λL1. 21 est le taux de défaillance par ouverture intempestive de l’embrayage, λL1.22 est le
taux de défaillance par perte de la transmission d’effort du composant interfaçant avec le tube
tournant de la jambe, λg1.31 est le taux de défaillance général du réducteur, λg1.32 est le taux de
défaillance général du moteur et λg1.33 est le taux de défaillance général des composants situés
en amont du moteur.
( ( ) (
λ L 0 = a 25 λ L1.21 + λ L1.22 + a22 λ g1.31 + λ g1.32 + λ g1.33 )) 2
( 2.17 )
50
CHAPITRE 2 – RECHERCHE D’ARCHITECTURES
( ) (
λ D 0 = 2a25 λ L1.21 + λ L1.22 + 2a 22 λ g1.31 + λ g1.32 + λ g1.33 ) ( 2.18 )
51
CHAPITRE 2 – RECHERCHE D’ARCHITECTURES
2,0,E-06
Inadmissible
1,5,E-06
1,0,E-06
Admissible
5,0,E-07
0,0,E+00
0,0,E+00 4,0,E-07 8,0,E-07 1,2,E-06 1,6,E-06 2,0,E-06
Taux de défaillance de perte de capacité de transmission d’effort des composants en
contact avec le tube tournant (/HO)
Figure 2.20 : Taux maximal admissible de défaillance par perte de la transmission d’effort des
éléments d’interface avec le tube tournant en fonction du taux maximal admissible de défaillance
générale des composants situés en amont du moteur, pour un taux du mode d’opération dégradé de 5
10-6 /HV.
Un taux de défaillance par grippage maximal admissible de 2 10-9 /HO pour les
composants en contact avec le tube tournant ;
Un taux de défaillance générale de l’ordre 10-6 /HO pour les composants situés en
amont du moteur.
L’étude de fiabilité précédente a aussi montré que l’objectif de fiabilité de 10-9 /HV
pour le taux d’apparition du mode d’opération dégradé n’est pas réaliste. Les calculs ont
permis d’identifier un objectif plus réaliste de 5 10-6 /HV, à partir duquel les deux dernières
spécifications précédentes ont été obtenues. Il s’agit ici d’un résultat important, que nous
avons pu obtenir très tôt dans le projet grâce à l’application de la méthode proposée.
52
CHAPITRE 2 – RECHERCHE D’ARCHITECTURES
2.5. CONCLUSION
Les objectifs et les exigences du projet DRESS sont très élevés, notamment du point
de vue sécurité, fiabilité et opérationnel (par exemple : réversibilité du mouvement de la
jambe lors du remorquage, déconnection totale de l’actionneur en mode de défaillance
ultime), et n’ont pas permis d’appliquer efficacement une méthodologie traditionnelle de
conception, basée sur une approche déductive (montante) et sur le retour d’expérience. C’est
pourquoi, nous avons développé ici une méthode hybride (descendante et montante) pour
engendrer par énumération et tri un nombre réduit d’architectures solutions compatibles avec
les diverses exigences du projet. En partant, d’un nombre très élevé (7860 au niveau
fonctionnel) d’architectures solutions, cette méthodologie a permis d’identifier trois
architectures types au niveau organique, qui correspondent aux exigences du projet en général
et aux exigences de sécurité et de fiabilité en particulier.
53
CHAPITRE 2 – RECHERCHE D’ARCHITECTURES
REFERENCES
[APEX Dynamics 2009] AF Series Catalogue, 2009. http://www.apexdyna.com/
[ARP 4761] Society of Automotive Engineering International. SAE ARP 4761 "Guidelines and
Methods for Conducting the Safety Assessment Process on Civil Airborne Systems and
Equipment", 1996.
[Currey 1985] N. S. Currey. Aircraft Landing Gear Design: Principles and Practices.AIAA Education
Series,J. S. Przemieniecki, Wright-Patterson Air Force Base, Ohio, AIAA, 1985.
[Dhillon 1988] B. S. Dhillon. Mechanical reliability: theory, models and applications, American
institute of aeronautics and astronautics inc., 1988.
[ETEL 2007] ETEL Motion technology. Torque Motors Data Sheets, TMB Series, 2007.
http://www.etel.ch/
[Hecht 2004] H. Hecht. Systems reliability and failure prevention, Boston, Artech House, 2004.
[JAR 25] Joint Aviation Authorities. JAR-25 Large Aeroplanes. Advisory Material - AMJ, 1994
[MIL-HDBK 217] Department of Defence, United States of America. Military Handbook Reliability
Prediction of Electronic Equipment, 1991.
[Monclair 2008] P. Monclair. Latest developments in electrical landing systems. 26th Congress of
International Council of the Aeronautical Sciences. Anchorage, Alaska, USA, 14 -19
September 2008.
[Multon 2006] B. Multon, H. Ben Ahmed, M. Ruellan et G. Robin. "Comparaison du couple massique
de diverses architectures de machines tournantes synchrones à aimants." Société de
l'Electricité, de l'Electronique et des Technologies de l'Information et de la Communication
(SEE)(8): 85-93, 2006.
[Pahl 2007] G. Pahl, W. Beitz, L. Blessing, J. Feldhusen, K.-H. Grote et K. Wallace. Engineering
Design : A Systematic Approach, London, Springer-Verlag London Limited, 2007.
54
CHAPITRE 2 – RECHERCHE D’ARCHITECTURES
[Pahl 2004] G. Pahl, W. Beitz, J. Feldhusen and K.H. Grote. Engineering Design: A Systematic
Approach,L. Blessing, London, Springer, 2004.
[Redex 2006] Redex Andantex. SRP, Servo-Reducteurs - Robotique & Haute technologie, 2006.
[RIAC 2006] Reliability Information Analysis Center. RIAC Automated Databook, v.2.2. Utica, NY,
USA,2006.
[Schwob 1969] M. Schwob, G. Peyrache. Traité de fiabilité, Masson & Cie, 1969.
[Scott 2005] J. Scott. "JetBlue A320 Emergency Landing." from http://www.aerospaceweb.org/, 2005.
[Tassinari 1995] R. Tassinari. L'analyse fonctionnelle. Collection A savoir, Afnor, Paris, 1995.
55
CHAPITRE 2 – RECHERCHE D’ARCHITECTURES
56
CHAPITRE 3 – DIMENSIONNEMENT EN PUISSANCE
CHAPITRE 3
DIMENSIONNEMENT EN
PUISSANCE
ACRONYMES........................................................................................................................ 58
NOMENCLATURE ............................................................................................................... 58
INDICES ................................................................................................................................. 59
EXPOSANT ............................................................................................................................ 59
3.1. APPROCHE DE MODELISATION ........................................................................... 61
3.1.1. Dimensionnement en puissance à l’aide de la simulation inverse ............................ 61
3.1.2. Modélisation paramétrique........................................................................................ 62
3.1.3. Choix de l’approche de modélisation........................................................................ 65
3.2. LOIS D’ECHELLE ....................................................................................................... 66
3.2.1. Principe de base ........................................................................................................ 66
3.2.2. Exemple de composant électrique : les moteurs brushless........................................ 68
3.2.2.1. Paramètre de définition ...................................................................................... 69
3.2.2.2. Paramètres d’intégration géométrique et de masse ........................................... 71
3.2.2.3. Paramètres d’opération limite............................................................................ 71
3.2.2.4. Paramètres de simulation................................................................................... 74
3.2.3. Exemple de composant mécanique simple : le réducteur de vitesse......................... 76
3.2.3.1. Paramètre de définition ...................................................................................... 76
3.2.3.2. Paramètre d’intégration géométrique et de masse............................................. 76
3.2.3.3. Paramètres d’opération limite............................................................................ 77
3.2.3.4. Paramètres de simulation................................................................................... 77
3.2.3.5. Combinaison d’étages de réduction ................................................................... 78
3.2.4. Exemple de composant mécanique complexe : les vis à billes ou à rouleaux .......... 81
3.2.4.1. Paramètre de définition ...................................................................................... 83
3.2.4.2. Paramètres d’intégration géométrique et de masse ........................................... 83
3.2.4.3. Paramètres d’opération limite............................................................................ 83
3.2.4.4. Paramètres de simulation................................................................................... 84
3.2.5. Validation.................................................................................................................. 85
3.2.6. Bilan et tableaux de synthèse .................................................................................... 91
3.3. GRANDEURS DIMENSIONNANTES ....................................................................... 95
3.3.1. Grandeurs dimensionnantes ...................................................................................... 95
3.3.2. Couple nominal des réducteurs de vitesse................................................................. 95
3.3.3. Couple RMS et modèles thermiques des moteurs électriques .................................. 96
3.4. IMPLEMENTATION DANS DYMOLA/MODELICA............................................. 99
3.5. CAS TEST : COMMANDE DE VOL PRIMAIRE (AILERON)............................ 102
3.5.1. Modèle de charge et profil de mission .................................................................... 103
3.5.2. Dimensionnement en puissance .............................................................................. 104
3.5.3. Résultats .................................................................................................................. 106
57
CHAPITRE 3 – DIMENSIONNEMENT EN PUISSANCE
ACRONYMES
NOMENCLATURE
58
CHAPITRE 3 – DIMENSIONNEMENT EN PUISSANCE
INDICES
a Axial J Joule
abs Absolu (régime intermittent) max Maximal
aim Aimant permanent min Minimal
Amb Ambiant n Ecrou
centr Centrifuge nom Nominal
cond Conducteur nut Ecrou
cont Continu (régime permanent) per Permanent (régime permanent)
crête Crête r Radial
decl Déclassé en température ref Référence
e Entrée RMS « Root Mean Square »
em Electromagnétique s Sortie
eq Equivalent tot Total
f Fin tr Transiant (régime intermittent)
fer Fer vib Vibration
EXPOSANT
‘ Paramètre étudié
* Rapport d’échelle
Moyenne
59
CHAPITRE 3 – DIMENSIONNEMENT EN PUISSANCE
Faut-il utiliser une transmission directe (en Anglais : « direct drive ») ou indirecte ?
Dans le cas d’une transmission indirecte, quel(s) rapport(s) de réduction faut-il choisir
pour la transmission mécanique ?
60
CHAPITRE 3 – DIMENSIONNEMENT EN PUISSANCE
61
CHAPITRE 3 – DIMENSIONNEMENT EN PUISSANCE
Profil de
Propagation fonctionnelle de puissance
mission
Angle Réseau
électrique
Couple
embarqué
Charge Source de
Propagation du dimensionnement en puissance
puissance
62
CHAPITRE 3 – DIMENSIONNEMENT EN PUISSANCE
Légende
Flux de paramètre sortant
Flux de puissance fonctionnelle
direct
Flux de puissance inverse
(dimensionnement)
ds,dn,db : Diamètres de vis, écrou and palier
Modèle de vis-à-billes Fnom : Force nominale de sortie
Nœud de F(t) : Force axiale sur la vis
(5) strk p Fnom Lh Nœud de Frated : Force nominale de sortie
simulation simulation
Fnom p strk dimensionnante
rotative Données
Données de
de translative Fmax : Force maximale admissible de sortie
référence
référence Modèle
Modèle de
de (4) i : Intégrité de la plage de définition
Modèle de (2) Js mn ηd ηi p dim.
dim. (booléen)
… réducteur Modèle Js : Moment d’inertie de la vis
Modèle
Modèle de
de simulation
simulation d’aileron
de vitesse ls,ln,lb : Longueurs de vis, écrou et palier
F (t),ω(t) [Fnom]ref [p]ref Lh : Durée de vie
mn : Inertie en translationde l’écrou
(3) Fmax ωmax strk … Lh … Fnom p mtot : Masse totale de la vis-à-billes
Calculer (10)
Vérifier
Vérifier SOA
SOA Calculer Calculer
Calculer Vérifier
Vérifier la
la p : Pas de la vis
force
force nom.
nom. fiabilité
fiabilité plage de déf.
plage de déf. R : Fiabilité
strk : Course
v : Intégrité de la SOA (booléen)
Dimensionnement(6) v ∆F ∆ω (8) Frated (9) R (7) i ∆F : Marge de force de sortie
en puissance vis ∆ω : Marge de vitesse de vis
ls,ln,lb mtot
Intégration (1) ηd, ηi : Rendements directs et indirects
écrou ds,dn,db
palier ω(t) : Vitesse angulaire de la vis
ωmax : Vitesse angulaire maximale
admissible de la vis
5 paramètres de simulation ;
63
CHAPITRE 3 – DIMENSIONNEMENT EN PUISSANCE
3 paramètres de masse ;
64
CHAPITRE 3 – DIMENSIONNEMENT EN PUISSANCE
Une approche concepteur pure consiste à traduire les processus et les règles de
conception sous la forme de programmes informatiques, qui guident le concepteur tout au
long de la définition du composant. L’avantage de cette approche est qu’elle capitalise le
savoir faire de conception. Cependant, elle nécessite l’intervention des spécialistes de chaque
domaine concerné pour développer et manipuler les modèles. De plus, en phase de conception
préliminaire, les détails requis par cette approche (par exemple : caractéristiques des
matériaux) ne sont généralement pas encore disponibles.
L’approche proposée est basée sur les lois d’échelle pour lier entre eux les différents
paramètres d’un modèle. L’un des principaux avantages de l’utilisation des lois d’échelle est
l’absence de bases de données à développer et à actualiser. A la place, les lois d’échelle ne
requièrent qu’une référence par modèle de composant pour refléter l’état de l’art d’une
technologie donnée. Ainsi, l’outil développé est plus simple à manipuler et à actualiser. De
plus, l’établissement des lois d’échelle basées sur les phénomènes physiques présents dans les
composants, permet de faire ressortir les mécanismes dominant la conception d’un système.
Un autre avantage est la possibilité d’atteindre le juste niveau de détail par des hypothèses de
65
CHAPITRE 3 – DIMENSIONNEMENT EN PUISSANCE
conception adéquates. Par exemple, l’hypothèse d’une gamme de produit obtenue par une
variation homothétique des dimensions d’un composant est généralement mise en avant
[Multon 2006]. En effet, cette hypothèse simplifie les relations entre les différents paramètres
et est adapté à la phase de conception préliminaire, où un bas niveau de détails est requis.
Ainsi l’utilisation des lois d’échelle se montre particulièrement avantageuse. Cependant,
l’établissement de ces lois sur la base des phénomènes physiques, qui dominent le
dimensionnement des composants, est une tâche difficile, qui nécessite une connaissance
approfondie de ces composants et de leurs technologies, ainsi que des contraintes physiques
les dimensionnant et limitant leur domaine de fonctionnement.
l * = l' l ( 3.1 )
66
CHAPITRE 3 – DIMENSIONNEMENT EN PUISSANCE
2r’
2r
l
l'
r* = l* ( 3.2 )
V ' πr 2l
V* = = 2
= r *2l * = l *3 ( 3.3 )
V πr ' l '
M = ρ m dV ⇒ M * = l *3
∫
( 3.4 )
J = ∫ r dM ⇒ J = l
2 * *5
GI 0
K= ⇒ K * = l *3 ( 3.5 )
l
Dans certains cas, il peut être difficile d’identifier et d’écrire d’une manière simple les
équations liant les paramètres représentatifs d’un problème donné. Dans ces cas particuliers, il
est toujours possible de procéder à une analyse dimensionnelle pour établir les lois d’échelle
[Taylor 1974 ; Jufer 1996 ; Szirtes 1997 ; Barenblatt 2003].
67
CHAPITRE 3 – DIMENSIONNEMENT EN PUISSANCE
Aimants permanents
Bobinage
Stator
Rotor
Bloc élémentaire
Epaisseur
Figure 3.4 : Moteurs brushless cylindrique (gauche) et annulaire (droite), basé sur [Multon 2006].
Les lois d’échelle développées pour les moteurs brushless cylindriques sont basées
sur une évolution homothétique de toutes leurs dimensions, ainsi qu’un nombre de
paires de pôles constant.
D’un autre côté, la conception des moteurs brushless annulaires est basée sur un
assemblage de blocs élémentaires (figure 3.4). Ils ont de ce fait un nombre de paires de
pôles variant proportionnellement avec le diamètre du moteur, ainsi qu’une épaisseur
d’anneau demeurant constante lors d’une variation du diamètre.
68
CHAPITRE 3 – DIMENSIONNEMENT EN PUISSANCE
où B est l’induction dans l’entrefer, Jel la densité de courant dans les conducteurs, r la distance
par rapport à l’axe de rotation et V est le volume des conducteurs.
Cette intégrale donne en notation relative par rapport à une géométrie de référence :
Les performances d’un moteur électrique sont essentiellement limitées par la capacité
à évacuer la chaleur engendrée par les phénomènes dissipatifs. On distingue deux grands
types de pertes de puissance dans les moteurs : les pertes liées à l’effet Joule (en Anglais :
« copper losses ») et les pertes fer (en Anglais : « iron losses »). L’effet Joule est la
manifestation thermique du passage d’un courant électrique à travers une résistance
électrique. Les pertes fer rassemblent les pertes par hystérésis et les pertes par courant de
Foucault (en Anglais : « Eddy current losses »). Les pertes par hystérésis sont liées au
phénomène d’hystérésis entre l’induction et le chams magnétique. Le courant de Foucault est
le courant induit dans un conducteur par la présence d’un champ magnétique variable. Au
couple nominal et à vitesse nulle, les pertes du moteur sont essentiellement les pertes PJ liées
à l’effet Joule dans les conducteurs :
PJ = ∫ ρ el J el dV
2
( 3.8 )
Pour une résistivité constante, la variation des pertes Joule relatives lors d’un
changement d’échelle peut être mise sous la forme suivante :
∆T = Rth PJ ( 3.10 )
69
CHAPITRE 3 – DIMENSIONNEMENT EN PUISSANCE
avec
LA
B ≈ Br ( 3.14 )
LA + e f
On suppose qu’il n’y a pas de chute de force magnétomotrice dans la culasse et que la
perméabilité relative de l’aimant est égale à 1. On a donc B* = 1 si LA* = ef* (similitude
géométrique) et Br* = 1 (même technologie d’aimant). On note que pour de tout petits
moteurs, il est difficile de maintenir e* = LA*, d’où une certaine limite pour l’évolution dans
les toutes petites dimensions. Par conséquent, les expressions (3.7) et (3.13) permettent
d’écrire :
70
CHAPITRE 3 – DIMENSIONNEMENT EN PUISSANCE
Cem
*1 3.5
, nom cas du moteur brushless cylindrique
l = *1 3
*
( 3.16 )
Cem,nom cas du moteur brushless annulaire
C em
*3 3.5
, nom cas du moteur brushless cylindrique
M = *2 3
*
( 3.17 )
C em, nom cas du moteur brushless annulaire
71
CHAPITRE 3 – DIMENSIONNEMENT EN PUISSANCE
Légende
Figure 3.5 : Zones d’opération typiques des moteurs brushless cylindriques (haut) et annulaires (bas).
A très haute vitesse, certains éléments du rotor, comme par exemple les aimants
permanents, risquent de se décoller ou de se détériorer sous l’effet centrifuge. Ce
phénomène limite mécaniquement la vitesse maximale ωtr,max que peut atteindre le
moteur en régime transitoire.
72
CHAPITRE 3 – DIMENSIONNEMENT EN PUISSANCE
La contrainte mécanique σcentr imposée par les effets centrifuges est exprimée par
l’équation suivante :
σ centr =
(ρ m ,aim (
Vaim ) ω 2 raim ) ( 3.19 )
S centr
où ρm,aim est la masse volumique des aimants permanents, Vaim leur volume, Scentr la section
contrainte par les effets centrifuges, ω est la vitesse de rotation du moteur et raim est la
distance entre les aimants et l’axe de rotation.
qui permet d’écrire
*−1 3.5
*−1
C em , nom cas du moteur brushless cylindrique
ω *
tr , max =l = *−1 3 ( 3.21 )
C em, nom cas du moteur brushless annulaire
Pour certains moteurs brushless, la limite de vitesse n’est pas imposée par le
risque de décollement des aimants mais par un risque de surchauffe en régime
permanent (Figure 3.5). Quand le moteur atteint la vitesse maximale en régime
permanent ωper,max, l’échauffement n’est dû qu’aux pertes fer.
L’équation générale des pertes fer donnée dans [Grellet 1989] permet d’exprimer leur
évolution sous la forme suivante :
où fel est la fréquence électrique du courant traversant les conducteurs, l’exposant b est une
grandeur expérimentale et Mcond est la masse des conducteurs. D’après [Grellet 1989] b a pour
valeur moyenne 1,5.
73
CHAPITRE 3 – DIMENSIONNEMENT EN PUISSANCE
(3.12) et (3.22) représentatives de l’évolution des pertes liées à l’effet Joule et des pertes fer,
qui caractérisent les deux précédents points de fonctionnement, il vient :
La fréquence électrique dans les conducteurs est liée à la vitesse de rotation du moteur
par le nombre de paires de pôles. Le nombre de paires de pôles est supposé constant durant un
changement d’échelle du moteur cylindrique et croissant proportionnellement avec le
périmètre du moteur annulaire. La vitesse maximale admissible du moteur limitée par les
pertes fer peut donc s’écrire :
* −1
l *−1 b C 3,5b cas du moteur brushless cylindrique
ω *per ,max = *−1 = em,nom ( 3.25 )
l *−1 3 cas du moteur brushless annulaire
C em
,nom
Cem
*5 3, 5
, nom cas du moteur brushless cylindrique
J = *4 3
*
( 3.26 )
Cem,nom cas du moteur brushless annulaire
74
CHAPITRE 3 – DIMENSIONNEMENT EN PUISSANCE
En régime permanent, les pertes totales Ptot sont la somme des pertes Joule et des
pertes fer :
PJ P fer
} }
Ptot = αCem + βωb
2
( 3.27 )
où α et β sont respectivement les coefficient des pertes Joule et des pertes fer.
*−5 3, 5
l *−5 C em,nom cas du moteur brushless cylindrique
α * = *−4 = *− 4 3 ( 3.28 )
l C em,nom cas du moteur brushless annulaire
2
C em ,nom
β=α b
( 3.29 )
ω per ,max
l *3 C em
*3 3, 5
,nom cas du moteur brushless cylindrique
β = *2 + b = * 2 + b
*
( 3.30 )
l C em ,3nom cas du moteur brushless annulaire
*−2 3, 5
C em , nom cas du moteur brushless cylindrique
R = *−2 3
*
th ( 3.31 )
C em,nom cas du moteur brushless annulaire
et
C em
*3 3.5
, nom cas du moteur brushless cylindrique
C = *2 3
*
th ( 3.32 )
C em,nom cas du moteur brushless annulaire
75
CHAPITRE 3 – DIMENSIONNEMENT EN PUISSANCE
Cem
*1 3.5
cas du moteur brushless cylindrique
τ = ,nom
*
th ( 3.33 )
1 cas du moteur brushless annulaire
D’où
F * = l *2 ⇔ C * = l *3 ( 3.34 )
l * = C *1 3 ( 3.35 )
M * = C* ( 3.36 )
76
CHAPITRE 3 – DIMENSIONNEMENT EN PUISSANCE
*
Cmax = C* ( 3.37 )
Les pertes par friction visqueuses augmentent avec la vitesse de rotation du réducteur
causant une augmentation de sa température interne. Ce qui a pour conséquence de réduire
l’épaisseur du film lubrifiant et d’amplifier les phénomènes de friction. Au-delà d’une
température critique, ce phénomène cause des dommages irréversibles. Cette limite de vitesse
ωmax liée une contrainte thermique est commune à de nombreux autres composants
mécaniques comme les roulements, qui sont par ailleurs présents en nombre dans un réducteur
de vitesse. Sur la base d’un rendement de fonctionnement supposé constant, dans une
première approximation, la puissance dissipée est proportionnelle à la puissance mécanique
transmise. En se basant sur une dissipation de la puissance thermique par conduction, la
contrainte d’un échauffement maximal admissible constant lors d’un changement d’échelle
conduit à l’expression suivante :
D’un autre côté, le modèle de simulation peut être réduit aux phénomènes parasites suivants :
Inertie totale des éléments rotatifs ramenée à l’arbre rapide (arbre d’entrée).
ηi = 2 − 1 ηd ( 3.39 )
77
CHAPITRE 3 – DIMENSIONNEMENT EN PUISSANCE
53
J * = C* ( 3.40 )
Figure 3.6 : Exemple de l’architecture d’un réducteur de vitesse composé de trois étages pour un
rapport de réduction de 1000.
1n
ktot ≤ k j ,max ( 3.41 )
C’est pourquoi, les étages supplémentaires sont dimensionnés pour avoir la même
durée de vie que le dernier étage vis-à-vis du phénomène de fatigue. En d’autres termes, on
base le dimensionnement des étages de réduction supplémentaires sur la durée de vie du
dernier étage. La durée de vie Lh d’un étage de réduction vis-à-vis de la fatigue peut être
78
CHAPITRE 3 – DIMENSIONNEMENT EN PUISSANCE
exprimée en fonction du couple nominal Cs,nom sur l’arbre lent et de la vitesse nominale ωe,nom
sur l’arbre rapide correspondant à l’utilisation sur le cycle de mission de la manière suivante
[Budynas 2007] (cf. section 4.2.3.2 du chapitre 4) :
Lh ∝ C s−,nom
a
ω e−,1nom ( 3.42 )
L’égalisation de la durée de vie du jème étage de réduction, pour tout j ≥ 2, avec celle
du premier étage en prenant en compte le rapport de réduction de chaque étage ainsi que le
rendement conduit à l’équation suivante :
a
C 1− nj
C sa,nom, j ωe ,nom, j = s ,jnom,n k ω e ,nom,n ( 3.43 )
k 1− n η n − j
d
où Cs,nom,j est le couple nominal sur l’arbre lent de l’étage de réduction j, ωe,nom,j est la vitesse
sur l’arbre rapide de l’étage de réduction j, Cs,nom,n est le couple nominal sur l’arbre lent du
dernier étage de réduction, k est le rapport de réduction du réducteur, ηd est le rendement
direct d’un étage de réduction, ωe,nom,1 est la vitesse nominale de l’arbre rapide du dernier
étage de réduction.
C s ,nom,3
a
C s ,nom, 2 ω e ,nom, 2 =
a
10ω e ,nom,3
10η d
a
( 3.44 )
a Cs ,nom,3
C s ,nom,1ω e ,nom,1 = 100ω e,nom,3
2
100 η d
Chaque étage de réduction est choisi parmi une même gamme de réducteurs, qui sont
tous construis pour une vitesse nominale donnée identique : la vitesse nominale de référence.
Par conséquent, la vitesse nominale de l’arbre rapide est la même pour chaque étage. En
réalité, en fonction du rapport de réduction, l’arbre rapide du premier étage tourne plus vite
que celui du dernier étage. Ce qui impacte sa durée de vie, comme le reflète l’expression
précédente, qui permet d’écrire :
C s*,nom,n
C *
s , nom , j = j 1
, j≥2 ( 3.45 )
1− 1−
n a n− j
k η d
79
CHAPITRE 3 – DIMENSIONNEMENT EN PUISSANCE
* C s*,nom,3
C s ,nom, 2 ≈
4,6η d
( 3.46 )
* C s*,nom,3
C ≈
s ,nom,1 21,5η 2
d
La contrainte de fatigue dans les étages supplémentaires peut rapidement être dépassée
par les contraintes de fabrication ou même thermiques. Pour cette raison, il existe
généralement une taille minimale lmin pour le dimensionnement des étages de réduction :
l *j = C *j1,nom
3
, l j ≥ lmin ( 3.47 )
La longueur totale du réducteur est la somme des longueurs de ses différents étages de
réduction. L’analyse détaillée des catalogues constructeurs montre une tendance à rationaliser
la fabrication en imposant des diamètres extérieurs identiques à tous les étages de réduction,
c'est-à-dire égale au diamètre du premier étage en général. Les conditions précédentes
fournissent l’expression suivante de la masse Mj d’un étage de réduction :
M *j = C1*,2nom
3 *
lj ( 3.48 )
La masse totale du réducteur est la somme des masses de ses différents étages de
réduction. De la même façon que pour la masse, l’expression du moment d’inertie d’un étage
de réduction peut s’écrire :
J *j = C1*,4nom
3 *
lj ( 3.49 )
n J *j
J =∑
*
eq j −1
( 3.50 )
j =1 2
n
k
Le rendement global direct est le produit des rendements des différents étages :
η d ,tot = η nd , j ( 3.51 )
D’un côté, la limite de couple maximal est donnée par le premier étage de réduction,
où cette contrainte est maximale. D’un autre côté, la limite de vitesse est donnée par le dernier
étage ou celle-ci est maximale. Par conséquent, il est proposé ici que les expressions de
80
CHAPITRE 3 – DIMENSIONNEMENT EN PUISSANCE
La vis,
L’écrou et
Le palier.
Une vis à billes ou à rouleaux est caractérisée par différents modes de fonctionnement,
selon que la vis ou l’écrou est en rotation [Karam 2007]. On se concentre ici sur le cas
particulier où la vis tourne, l’écrou se translate et l’effort axial résultant sur la vis est transmis
au palier, qui maintient l’assemblage en position.
Ecrou
Vis
Palier
81
CHAPITRE 3 – DIMENSIONNEMENT EN PUISSANCE
Du point de vue des phénomènes parasites, nous proposons d’inclure dans le modèle :
C dyn = X 1 Fr + Y1 Fa , si Fa Fr ≤ e
( 3.52 )
C dyn = X 2 Fr + Y2 F , si Fa Fr > e
où X1, X2, Y1, Y2 et e sont des coefficients constants, qui permettent de prendre en compte la
géométrie du roulement dans le calcul de la charge dynamique équivalente. Dans le cas des
roulements à billes, Y2 et e ne sont plus constants et varient en fonction de l’intensité de
l’effort radial [Budynas 2007 ; SKF 2009].
Transformation
Fr
rotation-translation
pv Fa Légende
Vis Js Palier
Transmission de
Ca, θa puissance rotative 1D
Figure 3.8 : Modèle de simulation proposé pour les vis à bille et à rouleaux.
82
CHAPITRE 3 – DIMENSIONNEMENT EN PUISSANCE
Par conséquent, les évolutions des dimensions et des masses des différents éléments
mécaniques constituant les vis à billes et à rouleaux, à l’exception de la vis, sont obtenues de
la même façon que celles des réducteurs de vitesse et tout composant mécanique simple :
l * = Fnom
*1 2
( 3.53 )
et
M * = Fnom
*3 2
( 3.54 )
M l = ∫ ρ m .dS ⇒ M l* = l *2 = Fnom
*
( 3.55 )
La masse de la vis est ensuite calculée en multipliant sa masse linéique par sa longueur
totale. La longueur totale de la vis est la somme de sa course, de la longueur de l’écrou et de
la longueur du palier. Finalement, la masse totale de la vis à billes ou à rouleaux est obtenue
en faisant la somme des masses de ces différents composants (vis, écrou et palier).
83
CHAPITRE 3 – DIMENSIONNEMENT EN PUISSANCE
0 = Fnom
* *
Fmax, ( 3.56 )
Un autre facteur limitant l’effort maximal de sortie est le risque de flambage (en
Anglais : « buckling ») de la vis. L’effort maximal qui peut être appliqué à la vis sous l’angle
du risque de flambage est obtenu à l’aide d’un calcul usuel du premier mode de déformation
(Euler) d’un cylindre monté sur deux supports simples. A cause d’une possible imprécision
dans le montage, le résultat du calcul précédemment mentionné est généralement pondéré par
un coefficient compris entre 3 et 5. A cela un facteur de montage est généralement ajouté, afin
d’adapter le résultat à différentes configurations de montages possibles.
Les lois d’échelle, ci-après, reflètent les limites de vitesse précédentes en incluant les
facteurs de sécurité déjà présents dans la référence. La limite de vitesse de la vis à bille ou à
rouleaux est donnée par la limite la plus basse, c'est-à-dire la plus restrictive.
*1 2 *−1
ω*max,vib = k m Fnom
et ω*max = min(ω*max,vib ; ω*max,nut )
ls
* *−1 2
( 3.57 )
ω max,nut = Fnom
où ωmax est la limite de vitesse de la vis, ωmax,vib est la limite de vitesse de la vis due aux
vibrations, ωmax,nut est la limite de vitesse de la vis due à la vitesse maximale admissible des
éléments roulants de l’écrou, ls est la longueur totale de la vis et km est le facteur de montage.
J n* = Fnom
*5 2
( 3.58 )
84
CHAPITRE 3 – DIMENSIONNEMENT EN PUISSANCE
J l = ∫ ρ m r 2 dS ⇒ J l* = l *4 = Fnom
*2
( 3.59 )
3.2.5. Validation
Les lois d’échelle établies dans les sections 3.2.2, 3.2.3 et 2.4 et résumées dans la
section 3.2.6 sont comparées ici aux données de catalogues constructeur [Danaher 2000 ; HD
2005 ; Sumitomo 2005 ; Parvex 2006 ; ETEL 2007 ; APEX 2009 ; SKF 2009] pour
validation. Les figures 3.9 à 3.13 montrent les résultats dans des graphes log-log, qui ont
l’avantage de pouvoir couvrir de nombreuses décades et de représenter les lois en puissance
par des droites.
Les figures 3.9 à 3.11 valident les lois d’échelle établies pour la masse et l’inertie des
réducteurs de vitesse et des moteurs brushless. La dispersion observée en début de gamme des
moteurs annulaires traduit le fait que le constructeur n’a pas conçu toute sa gamme de produit
à partir d’un changement d’échelle homothétique. En effet l’étude du catalogue correspondant
montre que le constructeur n’a pas allongé la structure des moteurs homothétiquement avec
les autres dimensions, mais a préféré étirer chaque moteur de diamètre donné, pour couvrir
une plus grande gamme de besoins. Comme le mentionne B. Multon dans [Multon 2006], ce
choix répond à la demande du marché pour des fortes capacités de couple tout en cherchant à
minimiser les coût de production. En effet, pour accroître le couple, en conservant les mêmes
dimensions radiales (mêmes tôles notamment), le constructeur peut allonger la structure du
moteur. Il peut également changer les dimensions radiales. De la même façon, on observe que
l’épaisseur de l’anneau n’est pas constante mais croît discrètement avec le diamètre. En outre,
la figure 3.9 montre que les différentes technologies de réducteurs ont des couples massiques
comparables, même si les réducteurs Harmonic Drive se distinguent sur ce point et sont suivis
de près par les réducteurs cycloïdaux. Les données concernant les moteurs cylindriques
incluent le corps du moteur, ce qui n’est pas le cas des moteurs annulaires qui s’en retrouvent
sensiblement avantagés du point de vue de la donnée de masse.
85
CHAPITRE 3 – DIMENSIONNEMENT EN PUISSANCE
1,E+03
Lois d'échelle
PARVEX série NX Moteurs
Kollmorgen série BH *
ETEL série TMB annulaires
APEX série AF
1,E+02 Harmonic Drive série CSG
SUMITOMO FC series
References
scaling danaher
Masse (kg)
Moteurs
1,E+01 **
cylindriques
Réducteurs
1,E+00
de vitesse
Figure 3.9 : Masse en fonction du couple nominal électromagnétique des moteurs brushless et du
couple nominal sur l’arbre lent des réducteurs de vitesse.
1,E+02
Lois d'échelle
Moment d'inertie arbre rapide (kg.m²)
PARVEX - NX annulaires
1,E+00 APEX - AF
Harmonic Drive - CSG
1,E-01
SUMITOMO - FC
References
scaling danaher
1,E-02
Moteurs
1,E-03
cylindriques
1,E-04
1,E-05 Réducteurs
de vitesse
1,E-06
1,E-07
1,E-01 1,E+00 1,E+01 1,E+02 1,E+03 1,E+04
Couple nominal électromagnétique (Nm) - Moteurs brushless
Couple nominal sur arbre lent (Nm) - Réducteurs de vitesse
Figure 3.10 : Moment d’inertie du rotor en fonction du couple nominal électromagnétique des moteurs
brushless et moment d’inertie de l’arbre rapide en fonction du couple nominal sur l’arbre lent des
réducteurs de vitesse.
86
CHAPITRE 3 – DIMENSIONNEMENT EN PUISSANCE
1,E+07
1,E+05
1,E+04
1,E+03
1,E+00 1,E+01 1,E+02 1,E+03 1,E+04
Couple nominal sur arbre lent (Nm) - Réducteurs de vitesse
Figure 3.11 : Élasticité de torsion réfléchie sur l’arbre lent en fonction du couple nominal sur l’arbre
lent des réducteurs de vitesse.
1,E+04
Moteurs
annulaires
Constante de temps thermique (s)
Moteurs
cylindriques
1,E+03
Lois d'échelle
PARVEX - NX
Kollmorgen - BH
ETEL - TMB
References
ref parvex
1,E+02
1,E-01 1,E+00 1,E+01 1,E+02 1,E+03 1,E+04
Couple nominal électromagnétique (Nm)
Figure 3.12 : Constante de temps thermique des moteurs brushless en fonction du couple nominal
électromagnétique.
La figure 3.12 montre que les données constructeur suivent globalement les tendances
tracées par les lois d’échelle pour les paramètres de simulation thermiques des moteurs
87
CHAPITRE 3 – DIMENSIONNEMENT EN PUISSANCE
brushless. On remarque sur la figure 3.12 que la gamme de moteurs annulaires est caractérisée
par environ huit constantes de temps différentes. Ce qui s’explique par le fait que cette
gamme n’a pas été basée sur un changement d’échelle homothétique, mais pour huit
épaisseurs d’anneaux distinctes sur toute la gamme.
De la même façon, la figure 3.13 confirme globalement les lois d’échelles établies
pour la vitesse de rotation maximale. On remarque cependant que les moteurs annulaires
ayant un couple nominal de plus de 40 Nm environ ne sont pas limités en vitesse par leurs
pertes fer. L’observation du catalogue constructeur montre que ces moteurs sont en réalité
limités par la tension maximale à leurs bornes. La limite de vitesse imposée par cette tension
maximale est représentée en pointillé sur la figure 3.13. On peut supposer ici que ce choix
résulte de la rationalisation de la gamme d’électronique de puissance associée aux moteurs.
Les lois d’échelle développées ici, ont pour but de représenter les caractéristiques intrinsèques
des composants. De ce fait, il est préférable de ne considérer que la loi d’échelle donnant la
limite en vitesse des moteurs annulaires vis-à-vis des pertes fer.
1,E+04
Vitesse maximale arbre rapide (rad/s)
Réducteurs
Vitesse maximale rotor (rad/s)
de vitesse
1,E+03
Moteurs
cylindriques
1,E+02
Figure 3.13 : Vitesse maximale du rotor des moteurs brushless en fonction du couple nominal
électromagnétique et vitesse maximale de l’arbre rapide des réducteurs de vitesse en fonction du
couple nominal sur l’arbre lent.
En plus de valider les lois d’échelle présentées, les figures précédentes mettent en
avant les caractéristiques supérieures de la technologie Harmonic Drive pour les réducteurs de
vitesse, notamment du point de vue masse et élasticité de torsion. On remarque cependant,
que les réducteurs Harmonic Drive sont suivis de près par les réducteurs cycloïdaux, qui
apparaissent comme une alternative intéressante pour les actionneurs hautes performances.
De la même façon que pour les moteurs brushless et les réducteurs de vitesse, les lois
d’échelle établies ici pour les vis à billes et à rouleaux sont comparées aux données
88
CHAPITRE 3 – DIMENSIONNEMENT EN PUISSANCE
constructeurs [Rexroth 2008 ; SKF 2009] dans les figures 3.14 à 3.17. Il apparaît
distinctement une limite à la validité des lois d’échelle pour les roulements et les paliers. Cette
limite se situe à de très hautes valeurs de l’effort dimensionnant pour les roulements, et à de
très basses valeurs pour les paliers. A un dimensionnement particulièrement petit ou grand, les
contraintes dimensionnantes peuvent changer. Par exemple, avec une charge dynamique
dimensionnante inférieure à 6 103 N, les paliers sont bien plus petits qu’indiqué par la loi
d’échelle établie ici. Puisqu’il est difficile et coûteux de fabriquer de très petits composants,
qui risquent de ne concerner qu’un marché de niche, cet écart par rapport à la loi d’échelle
établie peut être expliqué par des limites de taille dues aux contraintes de fabrication. Afin, de
prendre en compte ce type de limite de la validité des lois d’échelle établies, une vérification
des paramètres de définition en fonction des largeurs de gammes des constructeurs a été
implémentée dans la structure générale proposée pour la modélisation des composants Figure
3.2-(10).
1,E+03
Loi d'échelle
SKF-Roulements à bille, une rangée
1,E+02 SKF-Roulements à billes à contact oblique
Référence
1,E+01
Référence
Masse (kg)
Roulements
1,E+00
à billes
1,E-01
1,E-02
Largeur de la gamme de
produit autour de la loi
1,E-03 Limite supérieure
d'échelle tracée
de la loi d'échelle
= - 59%, + 62%.
1,E-04
1,E+02 1,E+03 1,E+04 1,E+05 1,E+06 1,E+07
Capacité de charge dynamique (N)
Figure 3.14 : Masse des roulements à une rangée de billes en fonction de la charge dynamique.
89
CHAPITRE 3 – DIMENSIONNEMENT EN PUISSANCE
1,E+02
Loi d'échelle
BoshRexroth - LA
SKF - BLRU Paliers
Reference
1,E+01
Largeur de la gamme de produit
autour de la loi d'échelle tracée
Masse (kg)
= - 46%, + 51%.
1,E+00
Référence
1,E-01
Limite inférieure de
la loi d'échelle
1,E-02
1,E+02 1,E+03 1,E+04 1,E+05 1,E+06
Force nominale de sortie (N)
Figure 3.15 : Masse des paliers en fonction de la force nominale de sortie des vis à billes et à rouleaux.
1,E+03
Lois d'échelle
BoschRexroth vis à billes
SKF vis à rouleaux
1,E+02 Référence
ref ball screw
Référence
vis à billes
Largeur de la gamme de Ecrous de la
Masse (kg)
Figure 3.16 : Masse des écrous en fonction de la force nominale de sortie des vis à billes et à rouleaux.
90
CHAPITRE 3 – DIMENSIONNEMENT EN PUISSANCE
1,E+03
Loi d'échelle
BoschRexroth - Vis à billes
SKF - Vis à rouleaux
Reference
ref ball screw
1,E+02
Masse linéique (kg/m)
Référence
vis à billes Vis de la vis
Largeur de la gamme de à rouleaux
1,E+01 produit autour de la loi
d'échelle tracée =
Largeur de la gamme de
+ 119%, - 62%.
produit autour de la loi
Vis de la vis à d'échelle tracée =
1,E+00 billes + 181%, - 70%.
Référence vis à
rouleaux
1,E-01
1,E+02 1,E+03 1,E+04 1,E+05 1,E+06 1,E+07
Force nominale de sortie (N)
Figure 3.17 : Masse linéique des vis en fonction de la force nominale de sortie pour les vis à billes et à
rouleaux.
Bien que les figures précédentes valident les lois d’échelle préalablement développées,
il apparaît fréquemment une gamme de produits dispersés autour de la droite (dans un plan
log-log) définie par la loi d’échelle et la référence. Cette apparente dispersion est due à la
tendance des constructeurs à étendre leur gamme de produits en faisant varier certaines
dimensions indépendamment des autres. Cependant, sur une grande échelle de variation, les
proportions sont globalement maintenues et les lois d’échelle sont suivies. Dans tous les cas,
toutes les lois d’échelle présentées conduisent à un produit existant ou concevable.
Finalement, il est important de garder à l’esprit que la validité des lois développées couvre
généralement plus de trois décades, alors que l’exploration de l’espace de solution d’un
actionneur en couvre rarement plus d’une.
91
CHAPITRE 3 – DIMENSIONNEMENT EN PUISSANCE
,nom = l ,nom = l
Couple électromagnétique * *3.5 * *3
nominal
Nm Cem Cem (3.15)
Paramètres d’intégration
Longueur et diamètre m l * = Cem
*1 3.5
,nom l * = Cem
*1 3
,nom (3.16)
Masse kg M =C * *3 3.5
em,nom M =C * *2 3
em,nom (3.17)
Paramètres de simulation
Moment d’inertie kg.m2 J * = Cem
*5 3.5
,nom J * = Cem
*4 3
,nom (3.13)
*−5 3.5 *−4 3
Coefficient de pertes Joule W/(Nm)2 α =C
*
em ,nom α =C
*
em ,nom (3.28)
*( 2 +b ) 3
Coefficient de pertes fer W/(rad/s)b β =C
* *3 3.5
em ,nom β =C
*
em ,nom (3.30)
*−2 3.5 *−2 3
Resistance thermique K/W R =C
*
th em ,nom R =C
*
th em ,nom (3.31)
−1
*−1 3
Vitesse max. continue (pertes
fer)
rad/s
ω *
=C
*
3, 5b ω*cont ,max = Cem ,nom (3.25)
cont ,max em ,nom
92
CHAPITRE 3 – DIMENSIONNEMENT EN PUISSANCE
C s*,nom,n
C s*,nom, j = j 1
, j≥2 (3.45)
1− 1−
k n a
η nd− j
Paramètres d’intégration
Longueur et diamètre m d * = Cs*1,nom
3
(3.35)
l j = Cs*,1nom , l j ≥ l j, min
* 3
(3.47)
ns
ltot = ∑ l j (-)
i =1
M *j = Cs*,2nom
3
.l *j , j ≥ 2 (3.48)
n
M tot = ∑ M j (-)
i =1
Paramètres de simulation
Moment d’inertie kg.m2 J n* = Cs*,5nom
3
(3.40)
J =C
*
i
*4 3 *
s ,nom i .l , i ≥ 2 (3.49)
n *
J
J eq* = ∑ j
j −1 (3.50)
j =1 2
n
k
Élasticité de torsion Nm/rad K =C
*
n
*
s , nom (-)
K *j = Cs*,4nom
3
.l *j −1 , j ≥ 2 (-)
Kj
K j , eq = j −1 (-)
2
n
k
n n
K eq = ∏ K j , eq ∑K j , eq , j≥2 (-)
j =1 j =1
93
CHAPITRE 3 – DIMENSIONNEMENT EN PUISSANCE
Tableau 3.3 : Lois d’échelle établies pour les roulements et les vis à billes et à rouleaux.
Paramètres d’intégration
Longueur, diamètre m l * = Fnom
*1 2
(3.53) l * = Fnom
*1 2
(3.53)
Masse kg M * = Fnom
*3 2
(3.54) M * = Fnom
*3 2
(3.54)
M l = Fnom
**
Masse par unité de longueur kg/m - - (3.55)
Paramètres de simulation
Moment d’inertie kg.m2 - - J * = Fnom
*5 2
(3.58)
Moment d’inertie par unité de
J l = Fnom
* *2
kg.m - - (3.59)
longueur
Paramètres d’opération limite
Force maximale N
*
Fmax = Fnom
*
(3.56)
*
Fmax = Fnom
*
(3.56)
*1 2 *−1
*−1 2
ω*max,vib = k m .Fnom .ls
Vitesse angulaire maximale rad/s ω *
=F (3.57)
*−1 2
(3.57)
ω*max,nut = Fnom
max nom
Tableau 3.4 : Paramètres des moteurs brushless, réducteurs de vitesse et vis à billes et à rouleaux.
D’un côté, le développement des lois d’échelle des moteurs brushless et des réducteurs
de vitesse a permis d’illustrer l’application de la méthode proposée avec des composants
simples, dans une première approximation. D’un autre côté, l’exemple des lois d’échelle pour
le modèle des vis à billes et à rouleaux a permis de montrer comment l’approche proposée
peut aussi être efficacement appliquée à des composants plus complexes, dans le sens où ces
composants sont composés de plusieurs éléments qu’il faut dimensionner et simuler en partie
séparément.
94
CHAPITRE 3 – DIMENSIONNEMENT EN PUISSANCE
1/ a
tf
f nom = 1
(t f − t0 )v t∫
f (t )a
v (t ) dt
( 3.60 )
0
où t0 est le temps au début du cycle de mission, tf est le temps à la fin du cycle de mission, et
v est la vitesse moyenne.
95
CHAPITRE 3 – DIMENSIONNEMENT EN PUISSANCE
tf
C2 (t ) ω1 (t ) dt
1
=3 ∫
3
(t f − t0 )ω1 t0
C 2, nom ( 3.61 )
où ω1 est la vitesse de l’axe rapide, C2 est le couple sur l’axe lent et ω1 est la vitesse moyenne
de l’axe rapide sur le cycle de mission.
Le résultat précédent est confirmé par les formules de calcul renseignées dans les
catalogues constructeurs de réducteurs cycloïdaux (par exemple : [Sumitomo 2005]).
Profil (A)
Cem
T°eq
T° ∆T …
…
T°Amb
0 t
Figure 3.18 : Couple électromagnétique et température interne d’un moteur brushless en fonction du
temps pour un profil de mission d’apparence continue.
96
CHAPITRE 3 – DIMENSIONNEMENT EN PUISSANCE
tf
Cem (t ) dt
1
= ∫
2
C RMS
(t f − t0 ) t0 ( 3.62 )
Profil (B)
Cem
T°max T°
T°Amb ∆T
0 t
Figure 3.19 : Couple électromagnétique et température interne d’un moteur brushless en fonction du
temps pour un profil de mission d’apparence discontinue.
Dans ce dernier cas, le calcul de couple RMS traditionnel n’est pas représentatif de la
température de fonctionnement du moteur. Cependant, il n’y a pas aujourd’hui d’équivalent
au couple RMS pour refléter la température du moteur dans le cas d’un profil de mission de
type (B). Pour palier à cette lacune, il est proposé d'étendre au cas du couple moteur requis
l'approche proposée par [Nfonguem 2006] pour le calcul de la puissance requise. Ce qui
conduit à définir le couple filtré CRMS,τth image de la température maximale du moteur durant
un profil de mission donné. Dans une première approximation, le comportement thermique du
moteur est associé à un modèle du premier ordre (cf. section 3.2.2). Sur cette base, le calcul
du couple CRMS,τth est :
−1 C em ( p )
2
= max L
C RMS , τ th
(3.63 )
1 + τ th . p
97
CHAPITRE 3 – DIMENSIONNEMENT EN PUISSANCE
*
CRMS = ∆T *1 2
( 3.64 )
où CRMS,decl est le couple RMS après déclassement en température, Tnom est la température
nominale d’opération du moteur, Tamb est la température ambiante et Tamb,ref est la température
ambiante de référence fournie par le constructeur.
Comme il a été mentionné ci-avant, les calculs du couple RMS moyenné traditionnel
et filtré sont tout deux basés sur une prédominance des pertes Joules dans le comportement
thermique du moteur. Ce cas est le plus général, cependant pour des applications hautes
vitesses par rapport aux caractéristiques du moteur, il se peut que les pertes fer jouent aussi un
rôle significatif dans l’échauffement du moteur. Dans ce dernier cas, le calcul du couple RMS
ne suffit donc pas.
Pour prendre en compte à la fois les pertes Joule et les pertes fers dans l’échauffement
du moteur, nous proposons d’implémenter les deux modèles thermiques distincts suivants
dans le modèle des moteurs électriques :
98
CHAPITRE 3 – DIMENSIONNEMENT EN PUISSANCE
Tableau 3.5 : Approches traditionnelles et proposées de prise en compte des phénomènes thermiques
dans le dimensionnement des moteurs électriques, selon le type de profil de mission.
99
CHAPITRE 3 – DIMENSIONNEMENT EN PUISSANCE
Electronique Réducteur
de puissance Moteur
Vis à
Réseau brushless rouleaux
électrique Aileron
Bras de levier
Electronique
de puissance
Moteur
brushless Vis à
Réducteur rouleaux
Figure 3.20 : Modélisation d’un actionneur électromécanique pour les surfaces de commandes de vol
primaires d’un avion (aileron).
(b)
(a)
(d)
(2)
(1)
(3)
(6)
(4)
(5)
(c)
Figure 3.21 : Modèle du réducteur de vitesse (a), son interface graphique (b) et sa structure interne (c),
(d).
100
CHAPITRE 3 – DIMENSIONNEMENT EN PUISSANCE
La partie (c) de la figure 3.21 montre la structure interne du modèle de réducteur, qui
suit l’architecture générale développée dans la section 3.1.2 et illustrée par la figure 3.2. On
distingue :
(2) Les lois d’échelle, qui calculent tous les paramètres nécessaires à la simulation
à partir des paramètres de définition.
Les dommages cumulés durant le cycle mission peuvent prendre la forme de fatigue,
d’usure, de vieillissement thermique, (etc.) selon la nature et le mode de fonctionnement du
composant. Dans la section 4.2 du chapitre 4, nous proposons de ramener ces différents types
de dommages à deux approches de calcul génériques : une pour les dommages mécaniques
(fatigue, usure) et une pour le vieillissement thermique.
101
CHAPITRE 3 – DIMENSIONNEMENT EN PUISSANCE
C (Nm)
SOA
Cmax
Marge de
couple
ωmax ω (rad/s)
Figure 3.22 : Exemple de profil de mission et de SOA dans le plan puissance quatre quadrants couple-
vitesse.
Aileron
SHA-Aileron
Aileron
Figure 3.23 : Actionneur servo-hydraulique des commandes de vol primaires (aileron) d’un Airbus
A320, images extraites de [Vorst 2007].
102
CHAPITRE 3 – DIMENSIONNEMENT EN PUISSANCE
de l’actionneur actif, celui-ci passe en mode amortisseur, pendant que l’actionneur jusque là
en mode amortisseur passe en mode actif, pour assurer la continuité de l’accomplissement de
la mission. Avec la technologie servo-hydraulique, l’amortissement est réalisé en reliant les
deux chambres du cylindre hydraulique par un orifice calibré. De cette façon, l’actionneur
amortissant n’empêche pas l’actionneur actif d’agir sur la surface de contrôle. En cas de
défaillances multiples impactant les deux actionneurs d’une même surface, l’amortissement
permet d’éviter le phénomène néfaste de battement de l’aileron (en Anglais : « flutter »). Dans
le cas d’une solution électromécanique, le mode amortisseur nécessite de débrayer la voie de
puissance concernée et implique la présence d’un moyen d’amortissement. Cependant, pour
des raisons de simplification, le débrayage et l’amortissement ne sont pas considérés dans
cette étude.
(2) (1)
(3)
Les composants mécaniques sont dimensionnés par rapport aux efforts et aux vitesses
maximaux, qu’ils doivent transmettre, et à la fatigue, qu’ils subissent tout au long de leur vie
opérationnelle. D’un autre côté, le moteur électrique est dimensionné par rapport au couple et
à la vitesse maximaux, qu’il doit fournir et atteindre, ainsi que par les contraintes thermiques
103
CHAPITRE 3 – DIMENSIONNEMENT EN PUISSANCE
(par exemple : température interne et couple RMS). Par conséquent, deux profils de mission
distincts sont utilisés pour le dimensionnement de l’actionneur (Figure 3.25) :
Un profil représentatif des efforts et des vitesses maximales de l’actionneur, ainsi que
de la fatigue des éléments mécaniques et
Profil deprofil
Mission mission
: Maximum“mécanique”
effort and speed ProfilMission
de mission “thermique”
profile : Thermal constraint
100 100
force (%)(%max)
force (%)(%max)
80 80
outputnorm.
outputnorm.
60 60
sortie
sortie
40 40
Force deNorm.
Force deNorm.
20 20
0 0
Figure 3.25 : Profils de mission mécaniques (gauche) et thermiques (droite) normalisés de l’actionneur
dans le plan puissance force-vitesse.
104
CHAPITRE 3 – DIMENSIONNEMENT EN PUISSANCE
vitesse) ainsi que les valeurs nominales des variables dimensionnantes (force
de sortie nominale) sont aussi fournies.
En suivant la procédure précédente, les composants sont dimensionnés les uns après
les autres en remontant la chaîne d’actionnement selon le principe de la vague dimensionnante
(cf. section 3.1.1). La méthodologie de dimensionnement proposée présente de nombreuses
itérations. Cependant, grâce au nombre limité de paramètres à manipuler par le concepteur et
à l’aide de l’outil de développé, ces itérations se font très rapidement.
Sens de progression du
Performances
dimensionnement
… validées?
oui non
(5.a) Composant (5.b) Itération
suivant
(1) Modélisation, (2) paramétrisation
(3)
(3) Simulation
Simulation inverse
inverse
105
CHAPITRE 3 – DIMENSIONNEMENT EN PUISSANCE
3.5.3. Résultats
Les résultats du dimensionnement en puissance sont résumés dans le tableau 3.6. A
partir de ces résultats, il est possible de réaliser des analyses efficaces de masse et
d’intégration géométriques. La masse totale de l’actionneur est donnée par la somme des
masses de chacun de ces composants, que fournissent les différents modèles correspondants.
Le tableau 3.6 inclut aussi la description des composants de référence utilisés pour les lois
d’échelle, ainsi que les composants commerciaux les plus proches des résultats obtenus. Il
apparait clairement, que malgré des références généralement très éloignées des résultats du
dimensionnement (par exemple : le couple nominal du réducteur de vitesse de référence est
plus de vingt fois plus grand que celui dimensionné), les caractéristiques géométriques et de
masse des composants obtenues par lois d’échelle correspondent assez précisément aux
composants du commerce.
A partir des dimensions fournies dans le tableau 3.6, la géométrie de l’actionneur peut
être projetée dans le profil de l’aile et ainsi permettre de vérifier son intégrabilité (Figure
3.27). La longueur maximale de l’actionneur est donnée par la distance entre les pivots
d’accrochage avec la structure porteuse (pivot 1) et d’accrochage à la charge (pivot 2). La
figure 3.27 montre que le montage en ligne de tous les composants ne permet pas d’intégrer
l’actionneur entre ces deux pivots. Une solution d’intégration pourrait être de monter le
moteur et le réducteur en parallèle de la vis à rouleaux à l’aide d’un engrenage à train droit,
qui serait placé entre la sortie du réducteur et l’entrée de la vis à rouleaux.
106
CHAPITRE 3 – DIMENSIONNEMENT EN PUISSANCE
Tableau 3.6 : Dimensionnement de l’actionneur pour la configuration de base (bras de levier de 0,045
m, durée de vie de 48 000 heures et configuration d’opération actif/amortisseur).
Palier
Moteur Réducteur
Ecrou
brushless épicyclique
Vis
Pivot 3
Bras de levier
Pivot 2
Pivot 1 Profil
aileron
Position angulaire
maximale.
Position angulaire
Profil aile
minimale.
0,1 m Position
horizontale
Figure 3.27 : Géométrie de l’actionneur dimensionné (cf. tableau 3.6) dans le profil de l’aile.
107
CHAPITRE 3 – DIMENSIONNEMENT EN PUISSANCE
On suppose que l’actionneur réalise en moyenne 35 000 fois le profil de mission décrit
dans cette section tout au long de sa vie opérationnelle. Ce qui correspondant à deux aller-
retour par jour pendant 23 années d’exploitation. Le profil de mission donne l’angle
d’orientation de la roue du train avant en fonction du temps ainsi que du couple requis. Les
trois principales phases d’utilisation du train avant sont représentées par :
108
CHAPITRE 3 – DIMENSIONNEMENT EN PUISSANCE
Couple (N.m)
Vitesse angulaire (rad/s)
Figure 3.28 : Modèle de charge de la direction du train avant et profil de mission implémentés dans
Dymola/Modelica pour la simulation inverse.
Afin de réduire la masse des moteurs électriques, leur couple d’opération a été
minimisé en maximisant le rapport de transmission global de la partie mécanique de
l’actionneur. De cette façon, le rapport de transmission global est défini pour que les moteurs
atteignent leur vitesse maximale durant le profil de mission. Cette démarche a été appliquée
en prenant en compte les limites d’opération des différents composants (notamment en
vitesse) et les exigences fonctionnelles, par exemple : la vis a rouleaux, la vis sans fin et le
réducteur droit doivent être réversibles afin de permettre le remorquage, ce qui limite leurs
rapports de transmission maximal admissible.
109
CHAPITRE 3 – DIMENSIONNEMENT EN PUISSANCE
Architecture A
Réducteur Electronique de
Embrayage cycloïdal Moteur puissance
brushless
Vis sans fin
Tube
tournant
Architecture B
Tube tournant
Architecture C
3.6.3. Résultats
Les résultats du dimensionnement sont comparés dans le tableau 3.7 :en termes de
masse pour chaque composant et chaque architecture. Afin d’illustrer la pertinence des lois
d’échelle utilisées, quelques composants de l’architecture (B) sont détaillés et comparés à un
produit industriel équivalent dans les tableaux 3.8 à 3.10.
110
CHAPITRE 3 – DIMENSIONNEMENT EN PUISSANCE
Résultats Composant
obtenus à commercialisé le
l’aide des lois plus proche :
d’échelle et PARVEX NX310
Modelica [Parvex 2006]
Couple nominal 1,9 Nm 2 Nm
Masse 2 kg 2,1 kg
Moment d’inertie 10 10-5 kgm² 8 10-5 kgm²
Diamètre 7,4 cm 7,1 cm
Longueur 14 cm 14,7 cm
Remarque : Le couple nominal du moteur ayant servi de référence est de 4 Nm.
Résultats Composant
obtenus à l’aide commercialisé le
des lois plus proche :
d’échelle et FC-A25
Modelica [Sumitomo 2005]
Couple nominal axe lent 315 Nm * 349 Nm
Rapport de réduction 179 179**
Vitesse axe rapide max. 5200 tr/min 5050 tr/min
Masse 6,8 kg 7,5 kg
Moment d’inertie 10 10-5 kgm² 13 10-5 kgm²
Longueur 56 mm 70 mm
Diamètre 134 mm 155 mm
Remarque : Le couple nominal axe lent du réducteur ayant servi de référence est de 1390 Nm.
* Dimensionné sur l’effort maximal transitoire.
** Pour le composant sélectionné, le catalogue constructeur ne fournit qu’un rapport de transmission maximal de
119. Il est donc nécessaire de demander au constructeur d’étendre la capacité de réduction de se produit à 179,
comme c’est le cas pour d’autres produits de ce constructeur.
111
CHAPITRE 3 – DIMENSIONNEMENT EN PUISSANCE
Résultats Composant
obtenus à commercialisé le plus
l’aide des lois proche :
d’échelle et GDF134P
Modelica [TELCOMEC 2009]
Couple maximal 630 Nm* 600 Nm
Masse 6,7 kg 6 kg
Moment d’inertie 18 10-3 kgm² 13 10-3 kgm²
Diamètre 146 mm 134 mm
Longueur 74 mm 83 mm
Remarque : Le couple maximal de l’embrayage ayant servi de référence est de 2000 Nm.
* Dimensionné sur l’effort maximal transitoire.
La simulation inverse montre que les pertes fer sont négligeables devant les pertes
Joule. De plus, les précédents dimensionnements du moteur électrique ont été basés sur
l’hypothèse de la répétition du profil de mission jusqu’à atteindre la stabilisation thermique de
ce dernier. En conséquence, les précédents dimensionnements prennent en compte le couple
RMS moyenné (cf. Tableau 3.5). Si on considère qu’entre le décollage et l’atterrissage, le
moteur a le temps de se refroidir, il convient alors de dimensionner le moteur sur l’évolution
en temps réel de sa température de fonctionnement durant le profil de mission. Dans ce cas, le
couple RMS filtré dimensionnant le moteur est réduit de 25% par rapport au couple RMS
moyenné précédent.
3.7. CONCLUSION
L’utilisation de la modélisation et de la simulation offre un fort potentiel de réduction
des coûts de développement et d’amélioration de la qualité de conception. Notamment, en
permettant d’étudier et de tester un plus grand nombre de concepts d’architecture et de
décisions techniques majeures avant de ne passer à la réalisation d’un prototype physique.
Cependant, malgré les nombreux travaux portant sur le prototypage virtuel, il y a un manque
de modèles de simulation adaptés à la conception préliminaire.
112
CHAPITRE 3 – DIMENSIONNEMENT EN PUISSANCE
113
CHAPITRE 3 – DIMENSIONNEMENT EN PUISSANCE.
REFERENCES
[Ahmed 2004] H. Ben Ahmed, B. Multon et P.E. Cavarec. Actionneurs linéaires directs et indirects :
performances limites Journées du club EEA : « Avion et électricité » Université Cergy-
Pontoise 18-19 mars 2004.
[Bals 2003] J. Bals, G. Hofer, A. Pfeiffer et C. Schallert. Object-Oriented Inverse Modelling of Multi-
Domain Aircraft Equipment Systems and Assessment with Modelica. 3rd International
Modelica Conference, Linköping, Sweden, Peter Fritzson, 2003.
[Cellier 1991] F. E. Cellier. Continuous system modeling, New York, Springer Verlag, 1991.
[ETEL 2007] ETEL Motion technology. Torque Motors Data Sheets, TMB Series, 2007.
http://www.etel.ch/
[Ferretti 1999] G. Ferretti, S. Filippi, C. Maffezzoni, G. Magnani et P. Rocco. "A modular approach to
dynamic virtual reality modeling of robotic systems." IEEE Robotics and Automation
Magazine 6(4): 13-23, 1999.
[Grellet 1989] G. Grellet. "Pertes dans les machines tournantes." Techniques de l'Ingénieur, 1989.
[Groß 2001] H. Groß, J. Hamann et G. Wiegärtner. Motion Diagrams. Electrical Feed Drives in
Automation: Basics, Computation, Dimensioning, MCD corporate Publishing, SIEMENS:
260-290, 2001.
[HD 2005] Harmonic Drive. Units CSG-2UH Digital Catalogue, 2005. http://www.harmonicdrive.net/
[Jasnocha 2001] U. Jasnocha et S. Haasb. "A collaborative environment based on distributed object-
oriented databases " Computers in Industry 29(1-2): 51-61 2001.
[Jufer 1996] M. Jufer. Design and Losses - Scaling Law Approach. Nordic Research Symposium
Energy Efficient Electric Motors and Drives, Skagen, Denmark, 1996.
[Kim 1990] S. W. Kim et J. S. Park. "Computer aided optimum motor selection for d.c. servo drive
systems." International Journal of Machine Tools and Manufacture 30(2): 227-236, 1990.
[Lacroux 1994] G. Lacroux. Les actionneurs électriques pour la robotique et les asservissements,
Paris, Technique et documentation Lavoisier, 1994.
[Multon 2006] B. Multon, H. Ben Ahmed, M. Ruellan et G. Robin. "Comparaison du couple massique
de diverses architectures de machines tournantes synchrones à aimants." Société de
l'Electricité, de l'Electronique et des Technologies de l'Information et de la Communication
(SEE)(8): 85-93, 2006.
[Pahl 2007] G. Pahl, W. Beitz, L. Blessing, J. Feldhusen, K.-H. Grote et K. Wallace. Engineering
Design : A Systematic Approach, London, Springer-Verlag London Limited, 2007.
[Rexroth 2008] Rexroth Bosch Group. Precision ball screw assemblies, 2008.
http://www.boschrexroth-us.com/
[Ricci 1983] G. Ricci. "Weight and rated characteristics of machines: positive displacement pumps,
motors and gear sets." Mechanism and Machine Theory 18(1): 1-6, 1983.
[Ricci 1992] G. Ricci. "Mass and rated characteristics of planetary gear reduction units." Meccanica
27(1): 35-45, 1992.
[Roos 2007] F. Roos. Towards a methodology for integrated design of mechatronic servo systems,
KTH, Machine Design: viii-203, 2007.
[Roos 2006] F. Roos, H. Johansson et J. Wikander. "Optimal selection of motor and gearhead in
mechatronic applications." Mechatronics 16(1): 63-72, 2006.
[Szirtes 1997] T. Szirtes. Applied Dimensional Analysis and Modeling, Toronto, Ontario Canada,
McGraw-Hill Professional, 1997.
[Taylor 1974] E. S. Taylor. Dimensional analysis for engineers, Oxford, Clarendon Press, 1974.
[van de Straete 1999] H. J. van de Straete, J. De Schutter et R. Belmans. "An efficient procedure for
checking performance limits in servo drive selection and optimization." Mechatronics,
IEEE/ASME Transactions on 4(4): 378-386, 1999.
[van de Straete 1998] H.J. van de Straete, P. Degezelle, J. De Schutter et R. J. M. Belmans. "Servo
motor selection criterion for mechatronic applications." Mechatronics, IEEE/ASME
Transactions on 3(1): 43-50, 1998.
[Vorst 2007] R. van der Vorst. A guide to space and aircraft actuators operational vibration testing.
LMS Conference Europe. Stuttgart, Allemagne, 17-18 avril.
CHAPITRE 4
117
CHAPITRE 4 - EFFET DE LA DUREE DE VIE ET DE LA FIABILITE SUR LE DIMENSIONNEMENT
ACRONYMES
NOMENCLATURE
INDICES
118
CHAPITRE 4 - EFFET DE LA DUREE DE VIE ET DE LA FIABILITE SUR LE DIMENSIONNEMENT
EXPOSANT
Moyenne
119
CHAPITRE 4 - EFFET DE LA DUREE DE VIE ET DE LA FIABILITE SUR LE DIMENSIONNEMENT
Les limites d’opération continue incluent les notions de durée de vie et de fiabilité.
Certaines applications, et notamment les applications aéronautiques, peuvent avoir des
exigences particulières en termes de durée de vie et de fiabilité. Il est donc important de
prendre en compte ces exigences dans l’approche de dimensionnement proposée, sous la
forme de paramètres ou de résultats de conception. En d’autres termes, il s’agit ici de
développer une capacité de conception fiabiliste à l’aide d’un dimensionnement probabiliste.
Dans cette optique, la section 4.1 de ce chapitre décrit les principaux types de
dommages rencontrés dans les différents composants d’actionneurs électromécaniques. Puis,
la section 4.2 fait le lien entre ces dommages les grandeurs dimensionnantes et la durée de vie
des composants. Dans la section 4.3, nous faisons le lien entre la durée de vie et la fiabilité.
Ce qui permet au final de lier le dimensionnement et la durée de vie à la fiabilité en fonction
des dommages causés durant le profil de mission. Les liens précédents sont illustrés avec des
exemples de composants mécaniques (réducteurs de vitesse et vis sans fin) et électriques
(moteurs brushless cylindriques). La section 4.4 est consacrée à l’implémentation des
méthodes de calcul développées dans l’outil logiciel présenté dans le chapitre précédent. Dans
la section 4.5, les nouvelles capacités de l’outil sont appliquées aux cas test des commandes
de vol et de la direction de train avant d’atterrissage. L’exemple des commandes de vol
illustre l’influence des paramètres de conception (par exemple : durée de vie) sur le
dimensionnement. D’un autre côté, l’exemple de l’orientation de train avant d’atterrissage
illustre le dimensionnement en fiabilité des composants critiques des architectures A et C (cf.
chapitre 2 section 2.4.7.1). Finalement, les résultats obtenus, la méthode et l’outil logiciel
développés sont discutés dans la section 4.6.
120
CHAPITRE 4 - EFFET DE LA DUREE DE VIE ET DE LA FIABILITE SUR LE DIMENSIONNEMENT
entretenu dans ces conditions et opère dans une gamme de températures raisonnables, alors
seule la fatigue métallique pourra être cause de défaillance.
La fatigue caractérise la présence et la propagation de microfissures dans le matériau
d’un dispositif soumis à des contraintes ou à des déformations répétées. La présence de
microfissures tend à réduire la résistance aux efforts. La fin du processus mécanique de
fatigue est caractérisée par une rupture, qui peut conduire à une défaillance fonctionnelle de
du dispositif. Par exemple : la rupture de la dent d’un engrenage (Figure 4.2) peut conduire à
une augmentation critique du jeu et de la précision de la transmission de mouvement, voire à
un grippage.
Comme l’illustre la figure 4.1, l’endurance d’un dispositif à la fatigue est caractérisée
par les trois domaines distincts suivants [Duprat 1997]:
Nous nous concentrons ici sur les domaines d’endurance limitée et illimitée.
R
défaillances =
90
observée %
Limite de
fatigue
log(L)
Dommages
cumulés
Figure 4.1 : Charge de service σ (force, couple, contrainte mécanique, etc.) en fonction de la durée de
vie L (heures, jours, nombre de tours, etc.) d’un composant mécanique soumis au phénomène de
fatigue.
121
CHAPITRE 4 - EFFET DE LA DUREE DE VIE ET DE LA FIABILITE SUR LE DIMENSIONNEMENT
Figure 4.2 : Rupture classique de la dent d’un engrenage due à un phénomène de fatigue induit par les
efforts de flexion causés par la transmission de couple, extrait de [Neal 2009].
Initialement, les surfaces ont des accélérations et des vitesses relatives nulles et
la majorité de la charge est portée par les aspérités de ces surfaces. Les
aspérités de contact sont semblables à des jonctions, qui doivent être rompues
avant de permettre tout mouvement. Pour cette raison, des forces de friction
initiales sont généralement relativement élevées.
Avec une vitesse relative croissante, le fluide lubrifiant se propage entre les
surfaces, porte une part plus importante de la charge normale appliquée, et les
forces de frictions diminuent en conséquence. Dans l’idéal, le film lubrifiant
permet d’éviter tout contact solide entre les surfaces en mouvement. Dans la
réalité, le film peut être mince à certains endroits et il y a une proportion
significative de contacts métal-métal, qui causent une usure importante.
L’usure conduit à l’accumulation de débris et peut accentuer l’exposition de la
surface dégradée à la corrosion. En elles mêmes, l’accumulation de débris et la
corrosion sont des facteurs d’accélération du phénomène d’usure (par
exemple : par usures abrasives et adhésives).
122
CHAPITRE 4 - EFFET DE LA DUREE DE VIE ET DE LA FIABILITE SUR LE DIMENSIONNEMENT
Finalement, en fin de vie, l’usure devient sévère et peu aboutir à une mise hors
d’usage ou à une défaillance fonctionnelle.
Les composants proposés dans le commerce ont des durées de vie qui varient d’un
type de composant et d’un constructeur à l’autre. De plus, il est rare que la durée de vie
« catalogue » corresponde à l’application considérée. Dans le cas des composants
mécaniques, la durée de vie est caractérisée par un nombre de cycles de charge (par exemple :
nombre de révolutions), mais on préfère généralement l’exprimer en unité de temps (heures)
pour une vitesse moyenne d’opération donnée. La problématique précédente s’applique aussi
à la vitesse moyenne d’opération. C’est pourquoi, dans cette section, nous nous attachons
dans un deuxième temps à adapter la charge nominale dimensionnante en fonction de la durée
123
CHAPITRE 4 - EFFET DE LA DUREE DE VIE ET DE LA FIABILITE SUR LE DIMENSIONNEMENT
de vie et, pour les composants mécaniques, de la vitesse moyenne d’opération, afin d’aider le
concepteur dans la sélection d’un composant catalogue compatible avec l’application.
Comme mentionné précédemment, nous nous concentrons ici sur les domaines
d’endurance limitée et illimitée. Dans une première approximation, il convient de considérer
que toute contrainte en-dessous de la limite de fatigue ne cause pas de dommage. De plus, le
graphe log-log de la figure 4.1 fait ressortir la caractéristique suivante pour le domaine
d’endurance limitée [Budynas 2007] :
Pour une fiabilité donnée (par exemple : 90%), la quantité de dommage cumulée est
fixe [Budynas 2007]. Ce qui permet d’écrire :
1/ p
tf tf
tf
Q = ∫ q (t )dt ⇔ σ p
nom , fat (t f − t 0 )v = ∫ σ(t ) v(t ) dt ⇔ σ nom , fat
p
= 1
(t f − t 0 )v t∫
σ (t ) p
v (t ) dt
( 4.3 )
t0 t0 0
où q(t) est le dommage instantané, σ (t) est l’effort de charge instantané, v(t) est la vitesse
instantanée et v est la vitesse moyenne sur tout le cycle.
124
CHAPITRE 4 - EFFET DE LA DUREE DE VIE ET DE LA FIABILITE SUR LE DIMENSIONNEMENT
(
σ *nom , fat = L*h,10 v )
* −1 p
( 4.4 )
4.2.2. Usure
Comme on a vu dans la section 4.1.2, le processus d’usure fait intervenir de nombreux
phénomènes dont les interactions peuvent être complexes à prendre en compte. Afin de
pouvoir évaluer la dégradation par usure induite par le profil de mission, il convient, dans une
première approximation, d’adapter le modèle classique de calcul du débit de débris formulé
par Achard [Peterson 1980], pour ne prendre en compte que les effets de la variation de la
charge et de la vitesse d’opération. Dans cette optique, nous proposons l’expression de la
durée de vie suivante :
1/ α
tf tf
σ αnom ,us (t f − t 0 )v
(β +1)
= ∫ σ(t ) v(t ) = ( ) ( )
(β +1) 1 (β +1)
dt ⇔ σ nom ,us (β +1) ∫
α α
σ t v t dt ( 4.6 )
(t f − t 0 )v
t0 t0
β +1
*−
σ *nom ,us = L*h−,10
1α
v α ( 4.7 )
où σnom,us est la charge équivalente (par exemple : force, couple, contrainte mécanique) vis-à-
vis du phénomène d’usure.
125
CHAPITRE 4 - EFFET DE LA DUREE DE VIE ET DE LA FIABILITE SUR LE DIMENSIONNEMENT
approche de calcul générique pour les phénomènes de fatigue et d’usure est développée. Dans
un deuxième temps, afin de préparer les calculs de fiabilité décris dans la section 4.3, cette
section traite aussi de l’évaluation de la durée de vie effective des composants.
4.2.3.1. Généralisation
Le modèle de calcul d’usure proposé précédemment a l’avantage de permettre de
reproduire le modèle de fatigue en choisissant une valeur de β égale à zéro, c'est-à-dire en
annulant l’influence de la vitesse sur l’intensité des dommages. Par conséquent, nous
proposons ici d’utiliser le modèle d’usure pour l’approche générale de calcul des dommages,
qui couvre la fatigue et l’usure. A partir du modèle d’usure la quantité de dommage est
calculée à partir du profil de mission de la manière suivante :
tf
∫ σ(t ) v(t )
β +1
Q=
α
dt ( 4.8 )
t0
tf
∫ σ(t ) v(t )
Q 1 β +1
q= =
α
dt ( 4.9 )
t f − t0 t f − t0 t0
1 β +1
−
σ nom = q α ν α ( 4.10 )
Finalement, afin de pouvoir procéder aux calculs de fiabilité décrits dans la section
4.3, on cherche à calculer la durée de vie Lh,10 d’un composant du commerce correspondant
au dimensionnement spécifié. A partir de la quantité de dommages équivalente par unité de
temps, le calcul de ce paramètre s’écrit :
126
CHAPITRE 4 - EFFET DE LA DUREE DE VIE ET DE LA FIABILITE SUR LE DIMENSIONNEMENT
tf
∫ C (t ) ω (t ) dt
1
=3
3
(t f − t0 )ω1
C 2,nom 2 1 ( 4.12 )
t0
( )
13
* −1 3 L ω1
C 2*,nom = L*h ,10 ω 1 ⇔ C 2,nom ,dim = C 2,nom h ,10 ( 4.13 )
Lh ,10,ref ω1,ref
144 42444 3
aL
où C2,nom,dim est le couple nominal dimensionnant adapté à la durée de vie spécifiée Lh,10 et à la
vitesse moyenne d’opération ω1 du dispositif, aL est le facteur de longévité, Lh,10,ref et
ω 1,nom ,ref sont la durée de vie et la vitesse moyenne d’opération correspondant au composant
du commerce.
1
tf
q 3
∫ C (t ) ω (t ) dt
1
q= =
3
et C 2,nom ( 4.14 )
t f − t0
2 1
t0 ω1
Les valeurs des coefficients α et β varient selon les constructeurs et les différents types
de produits qu’ils proposent. Par exemple : dans [CMD 2005], on observe que α vaut 10/3 et
127
CHAPITRE 4 - EFFET DE LA DUREE DE VIE ET DE LA FIABILITE SUR LE DIMENSIONNEMENT
que β est quasi-nul. Ce qui équivaut à un phénomène de fatigue. On note, que cette valeur de
α est semblable à celle des roulements à rouleaux ou coniques [Budynas 2007], qui sont
présents dans les vis-sans-fin et absorbent d’importants efforts pour le guidage de la vis. Dans
[Framo 2005], les roue et vis sont fournies indépendamment des roulements, on observe que α
vaut 2 et que β vaut -1,33. En effet, comme il a été mentionné dans la section 4.1.2, avec une
vitesse croissante, le fluide lubrifiant se propage entre les surfaces et porte une part plus
importante de la charge appliquée, et les forces de friction diminuent en conséquence. C’est
pourquoi, dans une certaine mesure, une vitesse plus élevée réduit les dommages d’usure.
En prenant l’exemple des vis sans fin proposées dans [Framo 2005], où α est pris égal
à 2 et β à -1,33, le calcul des dommages mécaniques par unité de temps prend la forme
suivante :
tf
∫ C (t ) ω (t )
1 − 0 , 33
q=
2
dt ( 4.15 )
t f − t0
2 1
t0
où C2(t) est le couple instantané sur la roue et ω1(t) est la vitesse instantanée de la vis.
0 ,165
C 2,nom = q 0,5 ω1 ( 4.16 )
φ
⇔ log(t ) = log(E ) +
φ
t = Ee kT
( 4.17 )
kT
où t est la durée de la réaction, T est la température à laquelle celle-ci se produit, E est une
constante obtenue expérimentalement, φ est l’énergie d’activation et k = 0,8617 10-4 (eV/K)
est la constante de Boltzmann.
φ 1 1
−
k T Tref
L =e
*
h
( 4.18 )
128
CHAPITRE 4 - EFFET DE LA DUREE DE VIE ET DE LA FIABILITE SUR LE DIMENSIONNEMENT
φ
−
Lh e kT
=V ( 4.19 )
tf φ tf φ
− 1 −
V = ∫e kT
dt ⇒ υ = ∫ e kT dt ( 4.20 )
t0
t f − t0 t 0
φ
Teq = − ( 4.21 )
k log(υ )
129
CHAPITRE 4 - EFFET DE LA DUREE DE VIE ET DE LA FIABILITE SUR LE DIMENSIONNEMENT
Tableau 4.1 : Classes d’isolation thermique des moteurs électriques (Lh = 20 000 h).
A partir des informations fournies par le tableau 4.1, les équations décrites dans la
section précédentes permettent le calcul de la durée de vie effective d’un moteur en fonction
de la température de son bobinage (Figure 4.4). Par exemple, pour un moteur de classe
d’isolation B (standard industriel), on obtient :
1 1
tf 9280 −
1 − 9280 9280
Teq 403,15
T (t )
υ= ∫ = − =
*
e dt , T et L e ( 4.23 )
t f − t 0 t0 log(υ )
eq h
1,E+06
1,E+05
Durée de vie log(heure)
Classe A
1,E+04 Classe B
Classe F
Classe H
1,E+03
1,E+02
90 110 130 150 170 190 210 230
Température de bobinage (°C)
Figure 4.4 : Durée de vie effective des moteurs électriques en fonction de la température du bobinage
pour les classes d’isolations thermiques A, B, F et H.
130
CHAPITRE 4 - EFFET DE LA DUREE DE VIE ET DE LA FIABILITE SUR LE DIMENSIONNEMENT
t = kChU − p Ch ( 4.24 )
où U est la tension appliquée aux bornes du dispositif, kCh et pCh sont des constantes
expérimentales, avec pCh qui est fonction de la température et de l’isolant considéré.
Les trois variables de fiabilité les plus couramment utilisées sont [Hecht 2004] :
Dans cette section, dans un premier temps, le modèle générique choisi pour les calculs
de fiabilité (distribution de Weibull) est présenté. Puis, l’approche de calcul proposée, pour
évaluer les différentes variables de fiabilité d’un dispositif en fonction de sa durée de vie
opérationnelle, est décrite. Pour ce faire, nous nous appuyons sur les résultats des calculs de
durées de vie effective développés dans la section précédente. Finalement, les expressions
obtenues sont illustrées avec les applications aux réducteurs de vitesse et aux moteurs
électriques.
1
Le MTBF (« Mean Time Between Failure ») est généralement confondu avec le MTTF (« Mean Time To
Failure »). En réalité le MTBF est la somme du MTTF et du temps moyen de réparation MTTR (« Mean Time To
Repair »). La confusion est due au fait que le temps moyen de réparation est généralement négligeable devant le
temps moyen de défaillance.
131
CHAPITRE 4 - EFFET DE LA DUREE DE VIE ET DE LA FIABILITE SUR LE DIMENSIONNEMENT
b
b −1 t
bt −
f (t ) = e θ
( 4.25 )
θ θ
b
t b −1
− bt 1
R(t ) = e θ
, λ(t ) = et MTTF = θΓ + 1 ( 4.26 )
θ θ θ
+∞
Γ(n ) = ∫ e − x x n −1dx ( 4.27 )
0
où R est la probabilité de succès pour durée de vie spécifiée, Rref est la probabilité de succès de
référence correspondant à la durée de vie effective.
132
CHAPITRE 4 - EFFET DE LA DUREE DE VIE ET DE LA FIABILITE SUR LE DIMENSIONNEMENT
θ = Lh,10 (− log(Rref ))
−1 b
( 4.29 )
log(Rspec )
1b
ici Lh,spec est la durée de vie spécifiée, qui correspond à la fiabilité Rspec.
1, 5
t
R(t ) = 0,9 50000
MTTF ≈ 2,24.10 h
5
( 4.31 )
λ(t ) ≈ 1,4.10 t
−8
133
CHAPITRE 4 - EFFET DE LA DUREE DE VIE ET DE LA FIABILITE SUR LE DIMENSIONNEMENT
Les valeurs de aR et aRc résumées dans le tableau 4.2, sont respectivement obtenues à
l’aide de l’expression (4.30) et (4.2), où b est prise égale à 1,5, p est prise égale à 3 et la
fiabilité de référence à 0,9 :
log(Rdes )
1 1, 5
a R = et a RC = (a R )
-1 3
( 4.32 )
log(0,9 )
Tableau 4.2 : Correction de la durée de vie effective par aR ou du couple dimensionnant par aRC en
fonction de la fiabilité Rdes désirée pour les réducteurs de vitesse.
Rdes 0,90 0,91 0,92 0,93 0,94 0,95 0,96 0,97 0,98 0,99
aR 1 0,93 0,86 0,78 0,70 0,62 0,53 0,44 0,33 0,21
aRC 1 1,02 1,05 1,09 1,13 1,17 1,24 1,31 1,45 1,68
t
R(t ) = 0,98 20000
MTTF ≈ 9,9.10 h
5
( 4.33 )
λ = 10 −6 h −1
134
CHAPITRE 4 - EFFET DE LA DUREE DE VIE ET DE LA FIABILITE SUR LE DIMENSIONNEMENT
Dans cette section, les calculs de fiabilité sont illustrés avec les exemples de
composants mécaniques et électriques, afin de présenter l’implémentation des deux grands
types de modèles de fiabilité proposés : mécanique et thermique. Dans un premier temps,
l’implémentation des calculs de fiabilité mécanique est illustrée avec le modèle des vis à
billes et vis à rouleaux, dont les différents éléments sont particulièrement sensibles au
phénomène de fatigue. Dans un deuxième temps, l’implémentation des calculs de fiabilité
thermique est illustrée avec le modèle des moteurs brushless cylindriques, dont le bobinage
est sensible aux dégradations thermiques.
135
CHAPITRE 4 - EFFET DE LA DUREE DE VIE ET DE LA FIABILITE SUR LE DIMENSIONNEMENT
(9) Mesure de la force axiale appliquée sur l’écrou, la vis, le palier, etc.
Rvr = Rv R p ( 4.34 )
A partir du résultat précédent, il est possible de procéder aux calculs des autres
variables de fiabilité (MTTF et λ) à l’aide des équations (4.31). Le système vis-écrou et les
roulements partagent les mêmes paramètres de fatigue (α=3 et β=0) que les réducteurs.
136
CHAPITRE 4 - EFFET DE LA DUREE DE VIE ET DE LA FIABILITE SUR LE DIMENSIONNEMENT
(12) (14)
(18)
(16)
(10)
(4)
(7) (8)
Figure 4.5: Modèle des vis à billes et à rouleaux implémenté dans Dymola/Modelica.
137
CHAPITRE 4 - EFFET DE LA DUREE DE VIE ET DE LA FIABILITE SUR LE DIMENSIONNEMENT
138
CHAPITRE 4 - EFFET DE LA DUREE DE VIE ET DE LA FIABILITE SUR LE DIMENSIONNEMENT
(24)
(2)
(23) (18) (20)
(16)
(17) (19)
(6) (8)
(5)
(4) (9)
(1)
(7) (14)
(12) (10)
(15)
(13) (11)
(3)
(21)
(22)
Figure 4.6: Modèle de simulation des moteurs électriques + variateur de vitesse implémenté dans
Dymola/Modelica.
139
CHAPITRE 4 - EFFET DE LA DUREE DE VIE ET DE LA FIABILITE SUR LE DIMENSIONNEMENT
Plus précisément, l’exemple des commandes de vol illustre l’influence des principaux
paramètres de conception (cinématique de charge, durée de vie, configuration de
fonctionnement) sur le dimensionnement de l’actionneur. D’un autre côté, l’exemple de
l’orientation du train d’atterrissage avant, illustre une approche basée sur le
surdimensionnement, pour atteindre de plus hauts niveaux de fiabilité. Ce dernier exemple
illustre aussi les limites d’une telle approche.
Tableau 4.3 : Masses de l’actionneur pour différentes longueurs de bras de levier, durées de vie et
pourcentages d’activité.
140
CHAPITRE 4 - EFFET DE LA DUREE DE VIE ET DE LA FIABILITE SUR LE DIMENSIONNEMENT
141
CHAPITRE 4 - EFFET DE LA DUREE DE VIE ET DE LA FIABILITE SUR LE DIMENSIONNEMENT
L’exigence d’un taux de défaillance de 10-9 /HO, donnée dans le chapitre 2, n’est pas
compatible avec la réalité d’un taux de défaillance des composants mécaniques, qui est
variable (croissant) dans le temps. C’est pourquoi, cette exigence est interprétée ici comme
étant la valeur du taux de défaillance atteinte à la fin de la vie du composant. Les composants
dimensionnés en fiabilité sont constitués de plusieurs éléments mécaniques (par exemple :
vis-écrous et palier pour la vis à rouleaux). Pour les besoins de la démonstration, l’exigence
du taux de défaillance a été répartie équitablement entre ces différents éléments. Les
roulements de la vis sans fin sont dimensionnés indépendamment. Ce qui permet de prendre
en compte leur limite de fatigue (cf. Figure 4.1). La fiabilité des roulements est ainsi garantie
en les dimensionnant pour une durée de vie infinie. Le tableau 4.4 résume les résultats
obtenus pour le dimensionnement en fiabilité de la vis à rouleaux et de la vis sans fin pour une
durée de vie complète (35 000 cycles), puis divisée par deux (17 500 cycles).
Le tableau 4.4 montre que la vis sans fin est moins impactée par le
surdimensionnement en fiabilité que la vis à rouleaux. La vis sans fin est composée d’une vis,
d’une roue et de leurs roulements. La vis et la roue sont dimensionnées par l’effort maximal à
transmettre, qui conduit à un taux de fiabilité satisfaisant (~8 10-10 /HO) par rapport à
l’objectif. Le dimensionnement en fiabilité de la vis sans fin n’impacte donc que ses
roulements, d’où la variation limitée de sa masse totale. Pour la même raison, la variation de
la durée de vie n’influence pas le dimensionnement de la vis sans fin. D’un autre côté, le
dimensionnement en puissance de la vis a rouleaux conduit à un taux de défaillance (~8 10-7
/HO) éloigné de l’objectif de fiabilité, d’où l’impact conséquent du surdimensionnement.
Les fiabilités et les durées de vie L10 des composants atteignent des ordres de
grandeur relativement importants par rapport aux données des catalogues. Il
convient de vérifier auprès du constructeur, s’il peut effectivement garantir des
142
CHAPITRE 4 - EFFET DE LA DUREE DE VIE ET DE LA FIABILITE SUR LE DIMENSIONNEMENT
fiabilités (de 99,97% à 99,99%) aussi élevées pour de longues durées de vie
(de 12,4 106 à 334 106 de tours).
La réduction par moitié de la durée de vie spécifiée n’a pas d’impact sur le
dimensionnement de la vis sans fin, et un impact limité sur le
dimensionnement de la vis à rouleaux. D’un autre côté, diviser la durée de vie
par deux implique une augmentation significative des temps et coûts de
maintenance.
4.6. CONCLUSION
Certaines applications, en particulier aéronautiques, ont des exigences de fiabilité et de
durée de vie très élevées. Il est donc important de pouvoir prendre en compte ces aspects dès
la phase de conception préliminaire.
Pour ce faire, dans un premier temps, nous avons développé dans ce chapitre une
approche de calcul de durée de vie et de fiabilité en fonction du dimensionnement et des
dommages imposés durant le profil de mission. Dans un deuxième temps, ces approches de
calculs ont été implémentées dans l’outil logiciel d’aide à la conception, déjà présenté dans le
chapitre précédent, et développé dans un environnement Dymola/Modelica. Puis, l’outil
réalisé a été utilisé pour le dimensionnement des cas test des actionneurs électriques des
commandes de vol et de la direction du train d’atterrissage. En plus de résultats intéressants,
ces cas test ont démontré l’intérêt de l’outil pour procéder à un dimensionnement en fiabilité
rapide et efficace et pour étudier l’influence de la durée de vie sur le dimensionnement. De
cette façon, pour chacun des cas test considérés, le dimensionnement en fiabilité de
l’actionneur a pris environ 45 minutes.
143
CHAPITRE 4 - EFFET DE LA DUREE DE VIE ET DE LA FIABILITE SUR LE DIMENSIONNEMENT
REFERENCES
[Brancato 1992] E. L. Brancato. "Estimation of Lifetime Expectancies of Motors." IEEE Electrical
Insulation Magazine 8(3): 5-13, 1992.
[CAVEX 2005] Flender Drives and Automation. Worm Gear Units, 2005. http://flender.com/
[Duprat 1997] D. Duprat. "Fatigue et mécanique de la rupture des pièces en alliage léger." Techniques
de l'Ingémieur, traité Génie mécanique, 1997.
[Dynabox 2008] Girard Transmission - Precision gearheads. Right angle servo gearheads, 2008.
http://www.girard-transmissions.com/
[Gieras 1997] J. F. Gieras et M. Wing. Maintenance. Permanent Magnet motor Technology - Design
and Applications. M. Decker: 373-379, 1997.
[Hecht 2004] H. Hecht. System Reliability & Failure Prevention, Artech House, 2004.
[Lacroux 1994] G. Lacroux. Les actionneurs électriques pour la robotique et les asservissements,
Paris, Technique et documentation Lavoisier, 1994.
[Neal 2009] Neal Consulting Engineers. "Gear Failures: Tooth breakage." from www.tribology.co.uk,
2009.
[Peterson 1980] M. B. Peterson et W. O. Winer. Wear Control Handbook, New York, American
Society of Mechanical Engineers (ASME), 1980.
[Schwob 1969] M. Schwob, G. Peyrache. Traité de fiabilité, Masson & Cie, 1969.
[Sharif 2006] K. J. Sharif, H. P. Evans et R. W. Snidle. "Wear modelling in worm gears." IUTAM
Symposium on Elastohydrodynamics and Micro-Elastohydrodynamics 134: 371-383, 2006.
144
CHAPITRE 5 - EVALUATION DES PERFORMANCES EN BOUCLE FERMEE
CHAPITRE 5
EVALUATION DES
PERFORMANCES EN BOUCLE
FERMEE
ACRONYMES................................................................................................................................... 146
NOMENCLATURE .......................................................................................................................... 146
INDICES ............................................................................................................................................ 146
5.1. PROBLEMATIQUE DE CONCEPTION.............................................................................. 149
5.1.1. Démarche adoptée................................................................................................................ 149
5.1.2. Boucle de position................................................................................................................ 151
5.1.3. Boucle de vitesse.................................................................................................................. 153
5.1.4. Boucle de courant................................................................................................................. 155
5.1.5. Boucles de vitesse et de courant........................................................................................... 158
5.1.6. Boucle de position et de vitesse ........................................................................................... 160
5.1.7. Conclusion partielle ............................................................................................................. 161
5.2. MODELISATION ET ANALYSE.......................................................................................... 162
5.2.1. Modèles évolutifs ................................................................................................................. 162
5.2.2. Influence du modèle d’élasticité de la transmission mécanique .......................................... 164
5.2.3. Prise en compte des interfaces mécaniques.......................................................................... 165
5.2.4. Implémentation dans Dymola/Modelica .............................................................................. 168
5.3. CAS TEST : ACTIONNEUR DE POUSSEE VECTORIELLE DU PREMIER ETAGE DE
PROPULSION DU LANCEUR EUROPEEN VEGA........................................................... 169
5.3.1. Cahier des charges................................................................................................................ 170
5.3.2. Profil de mission filtré.......................................................................................................... 171
5.3.3. Implémentation du modèle de charge dans Dylmola/Modelica........................................... 174
5.3.4. Recherche du rapport de réduction global minimisant la masse de l’actionneur à l’aide des
lois d’échelle ........................................................................................................................ 176
5.3.4.1. Recherche du rapport de réduction minimisant le couple crête d’un moteur donné pour
un point de fonctionnement critique donné ................................................................. 178
5.3.4.2. Recherche du rapport de réduction minimisant la masse du moteur, pour un moteur
choisi en fonction du couple crête ............................................................................... 180
5.3.4.3. Influence de la vitesse maximale admissible du moteur sur le choix du rapport de
réduction...................................................................................................................... 181
5.3.4.4. Conclusion partielle .................................................................................................... 182
5.3.5. Dimensionnement en puissance ........................................................................................... 183
5.3.6. Analyse modale.................................................................................................................... 184
5.3.6.1. Spécification des composants en fonction de l’objectif de rapidité de l’asservissement
en position ................................................................................................................... 185
145
CHAPITRE 5 - EVALUATION DES PERFORMANCES EN BOUCLE FERMEE
ACRONYMES
NOMENCLATURE
INDICES
146
CHAPITRE 5 - EVALUATION DES PERFORMANCES EN BOUCLE FERMEE
I Intégral R Résonance
M, m Moteur Red Réducteur
max Maximal ref Référence
N Notch T Rotulage de la tuyère
n Propre u Unitaire
opt Optimal vb Vis à billes
p Proportionnel VR Vis à rouleaux
147
CHAPITRE 5 - EVALUATION DES PERFORMANCES EN BOUCLE FERMEE
Une commande, généralement classique, c’est à dire élaborée à partir des approches
fréquentielles de l’automatique linéaire.
Des capteurs, qui permettent de mesurer les grandeurs à doser (par exemple : position,
vitesse, accélération, effort transmis à la charge) dans le cadre d’un asservissement en
boucle fermée.
(+ dissipation)
Vs
Ordre
Ωm
Figure 5.1 : Architecture matérielle usuelle d’un actionneur électromécanique avec sa commande.
148
CHAPITRE 5 - EVALUATION DES PERFORMANCES EN BOUCLE FERMEE
Dans le cadre d’une conception préliminaire, notre objectif ici est d’évaluer, voire de
spécifier les composants électriques, mécaniques et les interfaces, pour atteindre les objectifs
et les performances en boucle fermée du cahier des charges, sans entrer dans le détail des
capteurs et de la commande. En d’autres termes, on cherche à définir la chaîne de
transmission de puissance de l’actionneur, pour qu’elle soit naturellement bonne vis-à-vis des
performances en boucle fermée, sans connaître les détails des capteurs et de la commande.
Dans cette optique, dans la section 5.1 de ce chapitre, on cherche à propager les
objectifs de performance globaux de l’actionneur (vitesse d’asservissement en position) au
niveau des composants. La section 5.2 présente l’approche de modélisation proposée, pour
procéder à l’analyse modale. L’implémentation de cette approche est illustrée avec l’exemple
du modèle de simulation de la vis à rouleaux développé dans l’environnement
Dymola/Modelica. Puis, dans la section 5.3 l’approche logicielle mise au point est illustrée
avec un exemple industriel : l’actionneur de contrôle vectoriel de poussée du lanceur européen
VEGA (en Italien : « Vettore Europeo di Generazione Avanzata1 »). Finalement, les résultats
obtenus, la méthode et l’outil logiciel développés sont discutés dans la section 5.4.
1
En Français : « lanceur européen nouvelle génération ».
149
CHAPITRE 5 - EVALUATION DES PERFORMANCES EN BOUCLE FERMEE
I ΩM θCh
Kp Kp V
θcom + + + + +
K
- - KI + - KI +
p p
Boucle de courant
Boucle de vitesse
Boucle de position
Une boucle imbriquée est plus rapide qu’une boucle imbriquante. Ainsi, sur la figure
5.2, la boucle de courant est plus rapide que la boucle de vitesse, et la boucle de vitesse est
plus rapide que la boucle de position. En d’autres termes, une boucle externe est limitée en
rapidité par sa boucle interne. Dans la pratique, les actionneurs industriels asservis en
position, semblables à celui de la figure 5.2, sont typiquement caractérisés par [Ellis 2004]:
Une boucle de courant avec une bande passante de l’ordre de 800 Hz.
Une boucle de vitesse avec une bande passante de l’ordre de 75 à 200 Hz.
Une boucle de position avec une bande passante de l’ordre de quelques Hz à 30 Hz.
On cherche ici à étudier l’influence des modes naturels de l’actionneur sur la bande
passante de l’asservissement en position en boucle fermée. L’objectif final est de pouvoir
passer d’un objectif global, la rapidité de l’asservissement en position en boucle fermée, à un
objectif local, la spécification des composants. Pour ce faire, dans un premier temps, nous
négligerons les effets des boucles de correction internes (boucle de vitesse moteur et boucle
de courant) pour étudier l’influence de la transmission mécanique [Moscrop 1999 ; Roos
2006]. Dans un deuxième temps, nous étudierons l’influence de ces boucles internes. De cette
manière, il sera possible d’établir une approche de spécification des propriétés des parties
mécaniques et électriques de l’actionneur vis-à-vis de l’objectif de rapidité d’asservissement
en position, et en prenant en compte leurs influences respectives.
150
CHAPITRE 5 - EVALUATION DES PERFORMANCES EN BOUCLE FERMEE
Comme mentionné précédemment, dans cette première étude, on idéalise les boucles
de courant et de vitesse. Ce qui revient à les assimiler à des gains purs. La figure 5.3
représente le modèle de simulation de l’actionneur électromécanique considéré, dans cette
première étape de l’étude de l’influence des modes naturels sur la rapidité de l’asservissement
en position.
ΩM k
JM JCh Transmission
mécanique : Hθ(p)
+ -
Filtre réjecteur
2ξ p2 fv
1+ p+ 2
θ Ch,com + ωn ω n Ω M ,com ΩM 1 θM + + 1 Ω Ch 1 θ Ch
KN KΩ k
- 2ξ p2 JM p - + J ch p p
1+ N p + 2
ωn ωn
Boucle de vitesse
moteur fermée
151
CHAPITRE 5 - EVALUATION DES PERFORMANCES EN BOUCLE FERMEE
2ξ
1+ p
θ Ch ( p ) 1 ωn
H θ ( p) = = (5.1 )
Ω M ( p) p 2ξ p2
1+ p+ 2
ωn ωn
avec
k f ω
ωn = et ξ = v n ( 5.2 )
J Ch 2k
Un filtre réjecteur est un filtre coupe bande du deuxième ordre avec deux zéros et deux
pôles. Ce filtre atténue le signal autour d’une pulsation ωN donnée, dite pulsation de rejection.
Pour des performances optimales, les paramètres du filtre sont idéalement alignés sur ceux du
système à commander (ωN=ωn et ξ H = ξ ). Il en découle la fonction de transfert en boucle
ouverte du système asservi suivante :
2ξ H p2 2ξ 2ξ
1+ p+ 1+ p 1+ p
ωN ω ωn 2
ωn
FTBOθ ( p ) = N Ω
K K K K
N
= N Ω ( 5.3 )
p 2ξ p 2ξ p22
p 2ξ p2
1+ N p + 2 1+ p+ 2 1+ N p + 2
ωN ωN ωn ωn ωn ωn
où KΩ est le gain pur représentatif de la boucle de vitesse fermée, KN est le gain proportionnel
du filtre réjecteur et ζN est le coefficient d’amortissement du dénominateur du filtre réjecteur.
Typiquement, ζN est pris égal à 0,7, afin de minimiser le temps de réponse à l’échelon.
KN est choisi de manière à obtenir une marge de phase de 45° de la commande en position.
ωu = ω n (ξ 2
N )
+1 − ξ N ≈
ωn
2
( 5.4 )
152
CHAPITRE 5 - EVALUATION DES PERFORMANCES EN BOUCLE FERMEE
Comme l’illustre la figure 5.4, la fréquence de gain unitaire boucle ouverte est proche
de la bande passante du dispositif en boucle fermée. C’est pourquoi, dans une première
approximation, il convient de considérer que ces deux grandeurs sont semblables. En
conséquence, le résultat précédent illustre clairement l’influence directe de la pulsation propre
sur la bande passante.
Réponse en boucle
ouverte
Gain (dB)
Réponse en boucle
fermée Pulsation de gain
unitaire en boucle
ouverte
Fréquence (Hz)
Figure 5.4 : Exemple de diagramme de Bode (gain) de la réponse en boucle ouverte et en boucle
fermée d’un système d’ordre deux.
En conclusion, nous avons montré ici que, dans le cas où la boucle de vitesse
fermée à une bande passante très élevée et n’influe pas sur la boucle de position, la
majoration de la bande passante obtenue pour la boucle de position fermée est deux fois
plus petite que la pulsation propre du mécanisme. Nous pouvons donc spécifier
l’élasticité du mécanisme ou l’inertie de la charge en fonction de l’objectif de rapidité de
l’asservissement en position.
153
CHAPITRE 5 - EVALUATION DES PERFORMANCES EN BOUCLE FERMEE
Kp
ΩM , com + + + ΩM
Icom Cem 1 1 θM
KI KC
- KI + - JM p p
p
∆Ω +
+ fv
Correcteur Boucle de -
courant ∆θ +
PI +
k
-
Transmission
électromécanique: HΩ (p) 1 ΩCh 1 θCh
JCh p p
2ξ p2
1+ p+ 2
ΩM ( p) ω n , AR ω n, AR
H Ω ( p) =
1
=
(
Cem ( p ) p J M + J Ch ) 1+ 2ξ p2
p+ 2
( 5.5 )
ω n,R ω n,R
avec
k J J k
ω n,R = , J eq = M Ch et ω n , AR = ( 5.6 )
J eq J M + J Ch J Ch
où ωn,R est la pulsation de résonnance, ωn,AR est la pulsation d’anti-résonnance et Jeq est le
moment d’inertie équivalent de la transmission électromécanique exprimés au niveau de l’axe
moteur.
ω n, R > ω n, AR ( 5.7 )
La fonction de transfert HΩ(p) est caractérisée par une intégration pure combinée à une
paire de pôles et de zéros complexes. La paire de zéros complexes crée un phénomène d’anti-
résonnance, qui est suivi par le phénomène de résonnance causé par la paire de pôles
complexes. L’intégration pure cause un déphasage initial de -90°. La paire de zéros fait passer
le déphasage à +90° à la pulsation ωn,AR. Finalement, le phénomène de résonnance ramène le
déphasage à -90° à la pulsation ωn,R. La réponse fréquentielle (gain et phase) de la précédente
fonction de transfert HΩ(p) est illustrée par un diagramme de Bode sur la figure 5.6.
154
CHAPITRE 5 - EVALUATION DES PERFORMANCES EN BOUCLE FERMEE
-40 dB/decade
ωn,AR
+90°
Phase (°)
-90°
Fréquence (Hz)
Figure 5.6 : Réponse fréquentielle (gain et phase) de HΩ(p) dans un diagramme de Bode.
155
CHAPITRE 5 - EVALUATION DES PERFORMANCES EN BOUCLE FERMEE
Tension aux
FEM
bornes du moteur
E
I com + + Vcom KH - I
V 1
Kp
- KI + 1+ τ H p + R + Lp
p
U com ( p )
H el ( p ) =
KH 1R
= ( 5.8 )
I ( p) 1 + τ H p 1 + τ el p
τ el = L R ( 5.9 )
où Ucom est la commande de tension produite par le correcteur PI de la boucle de vitesse, I est
le courant moteur, τH est la constante de temps caractérisant le retard pur du hacheur, KH est le
gain pur du hacheur, R est la résistance électrique du moteur et L son inductance.
K
C Ω ( p ) = K p 1 + i ( 5.10 )
p
156
CHAPITRE 5 - EVALUATION DES PERFORMANCES EN BOUCLE FERMEE
[Parvex 2006]). Ce qui nous conduit à prendre une constante de temps τH égale à moitié de
leur somme soit également 125 µs. D’un autre côté, la constante de temps électrique du
moteur est typiquement de l’ordre de 0,01s. Par conséquent, c’est la constante de temps
électrique qui domine et qui est compensée par le correcteur PI de la boucle de vitesse. On en
déduit l’expression du gain intégral suivante :
1
Ki = =R L ( 5.11 )
τ el
I com ( p ) K p
FTBOI ( p ) =
1
= , avec K=KpKH ( 5.12 )
I ( p) L p (1 + τ H p )
I com ( p ) FTBOI ( p )
FTBFI ( p ) =
1
= = ( 5.13 )
I ( p) 1 + FTBOI ( p ) 2ξ p2
1 + el p + 2
ω el ω el
avec
K Lωel
ωel 2 = et K = ( 5.14 )
Lτ H 2ξ el
où ξel est le coefficient d’amortissement de la boucle de courant fermée et ωel est la pulsation
propre de cette boucle.
Généralement le gain du correcteur, soit ici K, est choisi pour obtenir un coefficient
d’amortissement égal à 0,7 afin d’obtenir un asservissement stable et rapide. Ce qui conduit à
l’expression de la fréquence propre de la boucle de courant fermée suivante :
1
ω el ≈ ( 5.15 )
1,4 τ H
157
CHAPITRE 5 - EVALUATION DES PERFORMANCES EN BOUCLE FERMEE
+ 1
Ω M ,com + I Cem 1 ΩM
Kp 2ξ p2 KC
- KI + 1 + el p + 2 f v,M + JM p
p ωel ωel
Figure 5.8 : Schéma bloc de la boucle de vitesse fermée prenant en compte la boucle de courant
fermée.
I com ( p )
FTBFI ( p ) =
1 1
= = ( 5.16 )
I ( p) 2ξ 1 + τ el p
1 + el p
ω el
avec
2ξ el
τ el = ( 5.17 )
ω el
158
CHAPITRE 5 - EVALUATION DES PERFORMANCES EN BOUCLE FERMEE
Ω( p ) K
H Ω ( p) =
1 1
= c ( 5.18 )
I com ( p ) f v , M 1 + τ el p 1 + τ M p
avec
JM
τM = ( 5.19 )
f v,M
1 f v,M
Ki = = ( 5.20 )
τM JM
Ω( p ) KcK p
FTBOΩ ( p ) =
1
= ( 5.21 )
Ω com ( p ) J M p (1 + τ el p )
Ω( p ) FTBOΩ ( p )
FTBFΩ ( p ) =
1
= = ( 5.22 )
Ω com ( p ) 1 + FTBOΩ ( p ) 2ξ Ω p2
1+ p+ 2
ωΩ ωΩ
avec
1 ωel J
ωΩ = = et K p = 2 M (5.23 )
2ξ Ω τ el 2ξ Ω 4ξ Ω τ el K c
159
CHAPITRE 5 - EVALUATION DES PERFORMANCES EN BOUCLE FERMEE
2ξ
1 1+ p
θ Ch,com - ΩM KN ωn θ Ch
2ξ Ω p2
+ 1+ p+ 2 p 2ξ N p2
ωΩ ωΩ 1+ p+ 2
ωn ωn
Boucle fermée de
vitesse moteur Influence de la transmission
mécanique et de la charge
Figure 5.9 : Schéma bloc de la boucle de position fermée incluant la boucle de vitesse.
θ Ch ( p ) K
FTBOθ ( p ) =
1 1
= N (5.24 )
θ Ch,com ( p ) p 2ξ p 2ξ 2
p2
1+ N p + 2 1+ Ω p + 2
ωn ωn ωΩ ωΩ
Pour obtenir un asservissement en position rapide et stable, nous imposons une marge
de phase φm de 45° comme critère de stabilité. Ce qui conduit à l’égalité suivante :
ω ω
2ξ N θ ,u 2ξ Ω θ ,u
ωn ωΩ
− 90° − arctan − arctan = −180° + φ m = −135° ( 5.25 )
ω θ2,u ωθ2,u
1− 2 1− 2
ωn ωΩ
160
CHAPITRE 5 - EVALUATION DES PERFORMANCES EN BOUCLE FERMEE
ωn /ωθ,u
16
14
Mécanisme de la transmission
12
10
0
0 2 4 6 8 10 12 14 16
Boucle de vitesse ωΩ/ωθ,u
Figure 5.10 : Pulsation propre relative du mécanisme ωn/ωθ,u en fonction de la bande passante relative
de la boucle de vitesse fermée ωΩ/ωθ,u.
161
CHAPITRE 5 - EVALUATION DES PERFORMANCES EN BOUCLE FERMEE
ω ω
2ξ N θ ,u 2ξ Ω θ ,u
ωn ωΩ ωΩ =
ω el
ω el ≈
1
arctan
ω2 + arctan ω2 = 45° (5.25) 2ξ Ω
(5.23)
1,4 τ H
(5.15)
1 − θ2,u 1 − θ2,u
ωn ωΩ
Dans cette section, nous étudions dans un premier temps les besoins et les
problématiques de simulation spécifiques au dimensionnement en puissance, à l’analyse
modale et à la synthèse de la commande. Dans un deuxième temps, nous proposons une
approche de modélisation, dont l’objectif et d’éviter d’avoir à réaliser et manipuler un modèle
de simulation différent pour chacun des types d’études considérés ici.
162
CHAPITRE 5 - EVALUATION DES PERFORMANCES EN BOUCLE FERMEE
d’une vis-à-billes, la course de l’écrou peut être limitée par la modélisation de butées au
niveau de son déplacement. Au cours de la simulation inverse, le déplacement imposé de la
charge peut correspondre à un déplacement de l’écrou, qui est plus important que la course de
la vis à bille limitée par les butées. Dans ce cas, nous proposons de ne pas prendre en compte
les butées pour la simulation inverse, et à la place de renseigner le concepteur sur la course
totale parcourue par l’écrou durant la simulation.
Nous avons montré dans la section 5.1 de ce chapitre, qu’il est possible d’évaluer les
performances de rapidité de l’asservissement en position, à partir des modes naturels de
l’actionneur. Nous proposons ici d’estimer les modes naturels d’un dispositif en se basant sur
la linéarisation du modèle de simulation physique. Ce qui implique de procéder par
simulation directe. Dans le cadre de l’analyse modale, le modèle de simulation de la vis à
billes doit être configuré, pour être représentatif de l’influence de ce composant sur les modes
naturels de la transmission mécanique. Pour ce faire, contrairement à la simulation pour le
dimensionnement en puissance, il est nécessaire ici de prendre en compte les élasticités du
composant. D’un autre côté, il n’est pas nécessaire d’inclure les effets des frottements pour
évaluer les modes naturels tant que l’amortissement est faible [Spinnler 2001]. Ces derniers
ne sont donc pas pris en compte ici. Le principe de la linéarisation est d’obtenir une
représentation linéaire du modèle autour d’un point de fonctionnement donné. Des
discontinuités durant simulation peuvent conduire à des aberrations dans les résultats de la
linéarisation (par exemple : gain infini). Il est donc préférable ici de ne pas prendre en compte
de phénomènes discontinus dans le modèle de simulation.
163
CHAPITRE 5 - EVALUATION DES PERFORMANCES EN BOUCLE FERMEE
grande flexibilité de configuration du modèle à l’aide d’une structure globale (Figure 5.14),
qui inclut l’ensemble des phénomènes physiques et des niveaux de détail à considérer durant
le dimensionnement en puissance, l’analyse modale et la synthèse de la commande. Afin de
passer aisément d’une étude à l’autre, nous proposons ici des modèles, qui disposent d’une
interface graphique permettant à l’utilisateur de choisir les différents phénomènes et niveaux
de détail pris en compte durant la simulation.
Configuration A Configuration B
ΩM k ΩM k
Inertie Inertie de
réducteur charge
Figure 5.11 : Modèles de simulation de la transmission électromécanique dans le cas d’une élasticité
placée côté moteur (A) et côté charge (B).
164
CHAPITRE 5 - EVALUATION DES PERFORMANCES EN BOUCLE FERMEE
Nous proposons ici d’évaluer l’écart relatif e entre les pulsations propres de la
transmission mécanique dans les deux cas de placement de l’élasticité identifiés
précédemment et illustrés sur la figure 5.11. Les pulsations propres de la transmission
mécanique dans les configurations A et B sont obtenus à partir de l’équation (5.2) et
s’expriment de la manière suivante :
k k
ωn , A = et ω n , B = ( 5.26 )
J Ch J Ch + J Re d
A partir des expressions précédentes, on obtient l’écart relatif entre les pulsations
propres suivant :
ω n, A − ω n ,B 1
e= =1− ( 5.27 )
ω n, A 1 + J Re d J Ch
Puisqu’aujourd’hui, on ne sait pas quel placement est le plus fidèle à la réalité, nous
proposons ici d’avoir une approche conservatrice dans l’estimation de la bande passante de la
boucle de position fermée, et de placer l’élasticité du côté du moteur dans les modèles de
simulation des composants de transmissions mécaniques (configuration B). Cependant, le
placement proposé n’est pas le cas le plus conservateur pour toutes les architectures
d’actionneur possibles. Il serait donc intéressant, dans de futurs travaux, de permettre à
l’utilisateur de tester les différents placements de l’élasticité dans le modèle.
165
CHAPITRE 5 - EVALUATION DES PERFORMANCES EN BOUCLE FERMEE
sont développés sur le principe du bilan d’effort nul, qui est caractérisé par les deux
hypothèses suivantes :
La somme des efforts appliqués au composant, c’est à dire entrant et sortant des ports
d’interface mécanique, est nulle.
Le mouvement relatif d’un composant par rapport à un support fixe est la différence
entre son mouvement propre et la déflection de son accrochage (ou suspension) par
rapport au support fixe.
Interface mécanique
translationnelle
Transformation de
mouvement Rotation-
Translation
Transmission de k p
Interface
puissance
J mécanique
mécanique rotationnelle
rotationnelle
Moment d’inertie
Élasticité en torsion
équivalent des éléments
équivalente Transmission de
en mouvement de la vis à
puissance mécanique
bille
translationnelle
Légende
Transmission mécanique Transmission mécanique
translationnelle rotationnelle
Figure 5.12 : Représentation cinématique du modèle de simulation des vis à billes et vis à rouleaux
avec ports d’interfaces mécaniques rotationnels et translationnels.
La figure 5.13 illustre, avec l’exemple des modèles développés pour les moteurs sans
balais et les vis à billes et à rouleaux, comment l’implémentation d’interfaces mécaniques
dans les modèles de simulation des composants permet de couvrir l’ensemble des cas
basiques d’asservissement en position d’actionneurs électromécaniques.
166
CHAPITRE 5 - EVALUATION DES PERFORMANCES EN BOUCLE FERMEE
Elasticité en (c)
translation
d’accrochage kanc
Masse de
avec la l’actionneur
MA
structure
porteuse
Cem kvb p
JM Jvb
kacc
Moment d’inertie
du moteur et de la
Élasticité en torsion de
charge
l’accrochage à la charge
Figure 5.13 : Modélisations d’actionneurs caractérisés par une transmission souple en torsion (a), un
accrochage moteur souple en torsion (b) et un accrochage vis à billes souple en translation (c).
167
CHAPITRE 5 - EVALUATION DES PERFORMANCES EN BOUCLE FERMEE
(7) Mesure de la force axiale appliquée sur l’écrou, la vis, le palier, etc.
2
Le passage d’un quadrant à l’autre est rendu continu à l’aide d’une fonction tangente hyperbolique.
168
CHAPITRE 5 - EVALUATION DES PERFORMANCES EN BOUCLE FERMEE
Modèle global
(2a) (7)
(3a) (8a)
(4) (5)
(1)
(2b) (9)
(3b)
(8b)
(2c)
(6)
Synthèse de la commande
Figure 5.14 : Modèle de simulation physique reconfigurable des vis à billes et vis à rouleaux
développé dans Dymola/Modelica.
169
CHAPITRE 5 - EVALUATION DES PERFORMANCES EN BOUCLE FERMEE
de dix mètres, le P80 est chargé de 88 tonnes de propergol solide. Ceci lui permet de disposer
d’une poussée maximale de 3000 kN et d’un temps de combustion de 107 secondes. Alors
que sur Ariane 5 le pilotage vectoriel de la poussée est assuré par des dispositifs à source de
puissance hydraulique, sur le P80 cette tâche est assurée par des dispositifs à source de
puissance électrique (en Anglais : « Power By Wire »).
Accrochage
1° étage de
propulsion
Actionneur
électromécanique
Tuyère
Accrochage
+/- 6°
La tuyère est montée via une liaison cardan sur l’étage de propulsion, ce qui permet un
mouvement à point fixe avec deux degrés de liberté angulaire. Les mouvements de la tuyère
sont contrôlés par deux actionneurs électromécaniques placés dans des plans respectivement
orthogonaux, et agissants chacun sur un degré de liberté (Figure 5.15). L’objectif de cette
étude est de dimensionner ces actionneurs électromécaniques de manière à répondre au mieux
aux spécifications présentées dans la sous-section suivante. Pour ce faire, dans une première
approximation, nous négligeons l’interaction entre les deux degrés de liberté de mouvement
de la tuyère.
170
CHAPITRE 5 - EVALUATION DES PERFORMANCES EN BOUCLE FERMEE
100
Position angulaire norm. (%max)
24
12
-12
-24
-100
0 20 40 60 80 100
Temps (%temps de mission)
Figure 5.16 : Profil typique durant un vol de la consigne de position angulaire de la tuyère en fonction
du temps.
Amortissement
visqueux de rotulage :
Élasticité selon
353 Nm/(rad/s)
l’axe de rotulage :
291 950 Nm/rad
Moment d’inertie de
la tuyère selon l’axe
de rotulage :
818 kgm²
171
CHAPITRE 5 - EVALUATION DES PERFORMANCES EN BOUCLE FERMEE
vitesses discontinues, et donc des pics de puissance illimités au niveau de l’actionneur. Pour
pouvoir procéder au dimensionnement, il est donc nécessaire de filtrer le profil de mission sur
la base des performances dynamiques spécifiées pour l’actionneur.
θ réel ( p )
G( p ) =
1
= ( 5.28 )
θ com ( p ) 2ξ p2
1+ p+ 2
ωf ωf
2,9
t 5% = ( 5.29 )
ωf
Tableau 5.3 : Temps de réponses du cahier des charges avec pulsation propres, accélérations et
vitesses maximales du système minimal équivalent du deuxième ordre.
172
CHAPITRE 5 - EVALUATION DES PERFORMANCES EN BOUCLE FERMEE
A l’aide d’un générateur d’état à partir de la consigne d’ordre deux (Figure 5.18), il est
possible de mettre le système minimal précédent sous la forme d’un schéma bloc faisant
intervenir uniquement des intégrateurs et des gains. Il est alors simple de choisir et de limiter
les variables d’état, qui décrivent le système à partir de cette forme. Le filtre obtenu permet de
transformer le profil de mission carré fourni en profil de mission filtré au deuxième ordre et
représentatif du comportement spécifié dans le cahier des charge (Figure 5.19). Le profil de
mission ainsi rendu réaliste peut être exploité pour le dimensionnement en puissance.
Cependant, pour pouvoir procéder au filtrage, il faut dans un premier temps identifier les
accélérations et les vitesses maximales correspondant aux exigences du cahier des charges et à
renseigner dans générateur d’état.
2ξ ω f
Figure 5.18 : Générateur d’état à partir de la consigne d’ordre deux représenté sous la forme d’un
schéma bloc.
Dans un premier temps, à partir des temps de réponses spécifiés en fonction d’une
consigne échelon dans le cahier des charges, nous avons extrait la pulsation propre du
système minimal d’ordre 2 équivalent. Dans un deuxième temps, nous injectons cette
pulsation propre dans le générateur d’état à partir de la consigne d’ordre deux, pour identifier
l’accélération maximale permettant de satisfaire les différents temps de réponses. Finalement,
de la même manière, l’accélération maximale identifiée est injectée dans le générateur d’état
et nous procédons à l’identification de la limite de vitesse maximale. Les résultats obtenus
sont résumés dans le tableau 5.3.
173
CHAPITRE 5 - EVALUATION DES PERFORMANCES EN BOUCLE FERMEE
6
Consigne
5
Temps (s)
174
CHAPITRE 5 - EVALUATION DES PERFORMANCES EN BOUCLE FERMEE
Mode de simulation
(7) (8) inverse
(4)
(1)
(2) (3)
(5)
(6)
Mode de simulation
directe
175
CHAPITRE 5 - EVALUATION DES PERFORMANCES EN BOUCLE FERMEE
l’actionneur normalisée
150
Force de sortie de
100
(% max.)
50
0
-50
-100
-150
0 10 20 30 40 50 60 70 80 90 100 110
150
Vitesse de sortie de
100
(% max.)
50
0
-50
-100
-150
0 10 20 30 40 50 60 70 80 90 100 110
176
CHAPITRE 5 - EVALUATION DES PERFORMANCES EN BOUCLE FERMEE
Mesure de variables
de flux thermique
Réducteur de
vitesse mécanique Modèle de rotulage
de la tuyère
177
CHAPITRE 5 - EVALUATION DES PERFORMANCES EN BOUCLE FERMEE
Minimiser : M M ( N )
Sous les contraintes :
(P )Ω M ,max, j (N ) ≤ Ω M ,max ( 5.30 )
≡ fonctionnement dans le SOA
C
M ,max, j
(N ) ≤ C M ,max
N ∈ Ν *
Transmission mécanique +
bras de levier par rapport à
l’axe de rotulage Élasticité selon
l’axe de rotulage
CM
Moment d’inertie
du moteur
Moment d’inertie de
la tuyère selon l’axe
de rotulage Amortissement
visqueux de rotulage
178
CHAPITRE 5 - EVALUATION DES PERFORMANCES EN BOUCLE FERMEE
J
(
C M = N J M + Ch2 &θ& + f v θ& + kθ
1
) ( 5.31 )
N N
avec
N Re d
N= l ( 5.32 )
pv
1 f v θ& + kθ
N opt = J CH + &θ&
( 5.33 )
JM
On remarque ici que pour une charge inertielle pure, l’expression précédente
correspond bien au critère classique de choix du rapport de réduction Nopt,classique :
J Ch
N opt ,classique = ( 5.34 )
JM
179
CHAPITRE 5 - EVALUATION DES PERFORMANCES EN BOUCLE FERMEE
C M,crête 1
0
1 6 N opt 11 16
Rapport de réduction
Figure 5.24 : Exemple de couple crête en fonction du rapport de réduction pour un point de
fonctionnement critique et un moteur donnés.
avec
θ& f v + θk
A = J Ch + ( 5.36 )
&θ&
*−1
C M* ,crête = N opt ( 5.37 )
180
CHAPITRE 5 - EVALUATION DES PERFORMANCES EN BOUCLE FERMEE
L’expression précédente montre que plus le rapport de réduction est élevé, plus le
couple crête du moteur et donc sa masse peuvent être diminués. En effet, bien que
l’augmentation du rapport de réduction accroisse l’accélération &θ& du moteur, cela entraîne
aussi une réduction du moment d’inertie du rotor, qui au final tend à réduire le couple crête.
D’un autre côté, comme nous l’avons montré dans la section 3.2.2.3 du chapitre 3 la
vitesse maximale admissible du moteur ΩM,max est liée à son dimensionnement et donc à son
couple crête. En illustrant la présente démonstration avec les moteurs sans balais cylindriques,
nous obtenons l’expression suivante :
Ω *M ,max = C M*−,1crête
3, 5
= N opt
*1 3, 5
( 5.39 )
Comme l’illustre la figure 5.25, au-delà d’un rapport de réduction Nmax/Ω, la vitesse du
moteur au point de fonctionnement donné est supérieure à la vitesse maximale admissible du
moteur. En d’autres termes, nous observons ici que la vitesse maximale du moteur limite le
rapport de réduction de la transmission mécanique.
3, 5
1 3, 5
N Ω M ,max, ref 2, 5
1
Ω M ,max, ref max/ Ω = N max/ Ω θ& ⇔ N max/ Ω = ( 5.40 )
N opt ,ref θ& N 1 2, 5
opt ,ref
181
CHAPITRE 5 - EVALUATION DES PERFORMANCES EN BOUCLE FERMEE
A
N opt ,ref = ( 5.41 )
J M ,ref
7 1
Ω M ,max, ref 5 J M ,ref 5
N max/ Ω = 2 ( 5.42 )
θ& A
4
Vitesse angulaire (rad/s)
ΩM,max
1 Zone SOA
moteur
ΩM
0
N max/ Ω
1 5 9 13 17
Rapport de réduction
(10) Le moteur est choisi de manière à ce que son couple maximal admissible
corresponde au couple crête de fonctionnement au point critique.
Avec cette méthode, il est aussi possible de prendre en compte des marges de vitesse
et de couple. Pour ce faire, il suffit de recalculer la vitesse et le couple maximal admissibles
du moteur de référence en prenant en compte ces marges.
182
CHAPITRE 5 - EVALUATION DES PERFORMANCES EN BOUCLE FERMEE
Entre autres, la figure 5.26 montre que seul le moteur est impacté par la variation de
l’accélération angulaire maximale de la charge. En effet, l’effet de l’inertie de la charge
ramenée à la vis à rouleaux et au réducteur de vitesse est négligeable. De plus, l’inertie de ces
composants de la transmission mécanique est aussi très réduite. C’est pourquoi, seul le
moteur, dont l’inertie est significative, est sensible à la variation de l’accélération de la
charge. De cette façon, le gain (accroissement) de masse de l’actionneur par rapport à la
référence (700 °/s2) est de -25 % pour une accélération angulaire de 600 °/s2 et de +47 % pour
une accélération angulaire de 860 °/s2.
183
CHAPITRE 5 - EVALUATION DES PERFORMANCES EN BOUCLE FERMEE
0,1 m
Moteur sans Réducteur
Accélération balais épicyclique
angulaire Palier Ecrou
maximale : Vis Masse totale :
Figure 5.26 : Masses totales et enveloppes géométriques de l’actionneur dimensionné pour une
accélération angulaire maximale de 600, 700 et 860 °/s2.
184
CHAPITRE 5 - EVALUATION DES PERFORMANCES EN BOUCLE FERMEE
par les modes naturels de l’actionneur. Nous étudions ici ces modes naturels à l’aide du
modèle de simulation de l’actionneur (Figure 5.22), et vérifions qu’une commande classique
de type boucles imbriquées peut permettre d’atteindre la bande passante requise au niveau de
la boucle de position fermée (41,43 rad/s).
La figure 5.28 montre clairement, que la boucle de position seule est caractérisée par
un premier mode naturel égal à ~6138,67 rad/s (~977 Hz). Ce mode de résonnance correspond
à un mouvement d’oscillation synchrone des inerties de la transmission (réducteur et vis à
rouleaux) et de la tuyère avec des déformations au niveau des élasticités du réducteur de
vitesse et du rotulage de la tuyère. Un calcul approximatif de la pulsation propre de ce mode
de résonnance confirme le résultat de la simulation :
k Re d + kT
ω θ ,1 = ≈ 975 Hz ( 5.43 )
J Re d + J VR + J T
où kRed est l’élasticité du réducteur réduite à l’axe moteur, kT est l’élasticité selon l’axe de
rotulage réduite, JRed le moment d’inertie du réducteur, JVR le moment d’inertie de la vis à
rouleaux réduit et JT le moment d’inertie du rotulage de la tuyère réduit.
185
CHAPITRE 5 - EVALUATION DES PERFORMANCES EN BOUCLE FERMEE
977 Hz
Phase (°)
Déphasage de -90°
Déphasage de
-270°
Fréquence (Hz)
Figure 5.28 : Diagramme de Bode (gain et phase) de la réponse en position de la tuyère en fonction de
la vitesse du moteur, dans le cas où on néglige l’influence des boucles de correction internes (vitesse et
courant).
186
CHAPITRE 5 - EVALUATION DES PERFORMANCES EN BOUCLE FERMEE
5.4. CONCLUSION
Dans les chapitres 3 et 4 nous avons traité du dimensionnement en puissance et de
l’aspect fiabilité. Dans ce cinquième chapitre nous avons adressé les performances
d’asservissement en position en boucle fermée sans entrer dans les détails de conception
d’une commande.
Pour ce faire, nous avons établi des liens entre les performances des différentes
boucles de correction d’une commande à boucles imbriquées typique. Dans un deuxième
temps, une approche de modélisation permettant d’étendre la capacité de l’outil logiciel
développé dans les chapitres 3 et 4 a été présentée. L’objectif de cette approche est d’étendre
187
CHAPITRE 5 - EVALUATION DES PERFORMANCES EN BOUCLE FERMEE
la capacité de l’outil proposé à l’analyse modale. Cette nouvelle capacité de l’outil à ensuite
été appliquée à l’étude du dispositif de contrôle vectoriel de poussée du premier étage de
propulsion du lanceur spatial VEGA. Dans cette étude, l’objectif de rapidité de la boucle de
position fermée a pu être propagé jusqu’aux performances de rapidité du hacheur. De cette
façon, nous avons démontré ici, que l’approche développée permet au concepteur de spécifier,
dès la phase de conception préliminaire, les différents composants d’un actionneur
électromécanique non seulement vis-à-vis de l’aspect puissance et fiabilité mais aussi
performance de rapidité en boucle fermée.
188
CHAPITRE 5 - EVALUATION DES PERFORMANCES EN BOUCLE FERMEE
REFERENCES
[Ahmed 2004] H. Ben Ahmed, B. Multon et P.E. Cavarec. Actionneurs linéaires directs et indirects :
performances limites Journées du club EEA : « Avion et électricité » Université Cergy-
Pontoise 18-19 mars 2004.
[Carnevale 2007] C. Carnevale and P. D. Resta. Vega Electromechanical Thrust Vector Control
Development. 43rd AIAA/ASME/SAE/ASEE Joint Propulsion Conference & Exhibit
Cincinnati, OH, USA, 8-11 July 2007.
[Coelingh 2000] H. J. Coelingh. Design support for motion control systems, University of Twente,
Enschede, The Netherlands. Ph.D.: xii-218, 2000.
[Dequidt 2000] A. Dequidt, J.-M. Castelain and E. Valde. "Mechanical pre-design of high
performance motion servomechanisms." Mechanism and Machine Theory 35: 1047-1063,
2000.
[Ellis 2004] G. Ellis. Control system design guide, Elsevier Academic Press, 2004.
[ESA 2009] European Space Agency. "European Space Agency Website - Launch Vehicles." from
http://www.esa.int/, 2009.
[Grellet 1989] G. Grellet. "Pertes dans les machines tournantes." Techniques de l'Ingénieur, 1989.
[Jufer 2001] M. Jufer. Global design of an electric drive. Electromotion 2001. Bologna, Italy. 8: 55-
58, June 19-20, 2001.
[Moscrop 2001] J. W. Moscrop, C. Cook and P. Moll. Control of servo systems in the presence of
motor-load inertia mismatch. EICON'01: The 27th Annual Conference of the IEEE Industrial
Electronics Society. 1: 351-356, 29 November -2 December 2001.
[Moscrop 1999] J. W. Moscrop, C. Cook and F. Naghdy. Comparative Study of Feedback Sensor
Locations dor Servo Systems with Motor/Load Inertia Mismatch. 1999 Australian Conference
on Robotics & Automation, Brisbane, Australia March 20 - April 1, 1999.
[Paredis 2001] C. J. J. Paredis, A. Diaz-Calderon, R. Sinha and P. K. Khosla. "Composable Models for
Simulation-Based Design." Engineering with Computers 17(2): 112-128, 2001.
[Park 1989] J.S. Park and S.W. Kim. "Optimum Speed Reduction Ratio for DC Servo Drive Systems."
Int. J. Mach. Tools Manufact. 29(2): 207, 1989.
189
CHAPITRE 5 - EVALUATION DES PERFORMANCES EN BOUCLE FERMEE
[Richard 1998] W. Richard and Jr. Armstrong. Load to Motor Inertia Mismatch: Unveiling The Truth.
Drives and Controls Conference 1998. Telford England, 1998.
[Roos 2006] F. Roos and J. Wikander. The influence of gear ratio on performance of
electromechanical servo systems Mechatronic Systems. 4th IFAC Symposium on Mechatronic
Systems, Heidelberg, Germany, 2006.
[Saner 2004] C. Saner. A design guideline for robust high performance axes. The IEEE International
Conference on Mechatronics, ICM '04, 2004.
[Townsend 1989] W.T. Townsend and J. K. Salisbury. Mechanical bandwidth as a guideline to high-
performance manipulatordesign. IEEE International Conference on Robotics and Automation,
14-19 May 1989.
[Van Amerongen 2003] J. van Amerongen and P. Breedveld. "Modelling of physical systems for the
design and control of mechatronic systems." Annual Reviews in Control 27(1): 87-117, 2003.
[Van de Straete 1999] H.J. van de Straete, J. De Schutter and R. Belmans. "An efficient procedure for
checking performance limits in servo drive selection and optimization." Mechatronics,
IEEE/ASME Transactions on 4(4): 378-386, 1999.
[Van de Straete 1998] H.J. van de Straete, P. Degezelle, J. De Schutter and R. J. M. Belmans. "Servo
motor selection criterion for mechatronic applications." Mechatronics, IEEE/ASME
Transactions on 3(1): 43-50, 1998.
190
CONCLUSION
CONCLUSION
191
CONCLUSION
Dans un contexte très compétitif, les secteurs aéronautiques et automobiles font face à
de nombreux défis. Des innovations rapides doivent conduire à des performances plus
élevées, un meilleur respect de l’environnement, une plus grande sûreté d’opération, ainsi
qu’une réduction des coûts de développement et d’opération. Parmi les nombreuses pistes
envisagées pour atteindre ces objectifs, ces travaux s’inscrivent dans le cadre de l’effort de
recherche pour le développement de systèmes embarqués plus électriques. Dans ce contexte et
pour les deux projets de recherche finançant ces travaux (projets DRESS et C6E2), nous
développons et étudions des actionneurs électriques embarqués, asservis en position et ayant
une puissance typique allant d’une centaine de watts à quelques kilowatts. Afin de réduire les
temps et coûts de développement de ces dispositifs, nous proposons ici une méthodologie de
conception innovatrice. Sur la base du processus de développement en V, cette méthodologie
a pour objectif de faciliter et d’automatiser les boucles d’itération dès la phase de conception
préliminaire. Classiquement, ces itérations sont réalisées durant la phase montante du
processus de développement. Il nous faut donc transférer une partie de la connaissance et de
la compétence de la phase montante vers la conception préliminaire.
La première étape du processus de conception préliminaire, tel qu’il est abordé ici, est
la recherche de solutions architecturales. Le chapitre 2 traite de cette étape, pour laquelle nous
proposons une approche hybride (montante et descendante) en alternative à l’approche
traditionnelle intuitive et montante. L’approche traditionnelle consiste à recombiner de
manière intuitive des solutions existantes pour obtenir de nouvelles fonctionnalités. Ce qui
permet d’utiliser des solutions connues et réduit les risques de développement. D’un autre
côté, l’aspect intuitif ne garantit pas d’identifier la meilleure solution possible. L’approche
innovatrice proposée combine une génération systématique des solutions architecturales à
partir du besoin fonctionnel (approche descendante) avec un filtrage par rapport aux
exigences du projet et de l’état de l’art technologique (approche montante) et sur des critères
qualitatifs. De cette façon, l’approche proposée conduit à un nombre réduit de solutions
architecturales compatibles avec les principales exigences du projet et l’état de l’art
technologique. Cette approche est appliquée avec succès au projet DRESS, qui est caractérisé
par des exigences très élevées (sûreté, fiabilité et opérationnel). Dans un deuxième temps, en
raison d’exigences de sûreté et de fiabilité très élevées (taux de défaillance de l’ordre de 10-9
par heure de vol), une étude de la fiabilité des solutions identifiées est menée au niveau
architectural. Cette étude permet de conclure sur la pertinence des solutions étudiées vis à vis
de l’aspect fiabilité et de mettre en évidence l’extrême criticité de certaines exigences du
cahier des charges.
192
CONCLUSION
représentatives des phénomènes dimensionnants est une tâche difficile. Elle nécessite une
connaissance approfondie des composants et de leurs technologies. L’approche de
modélisation proposée est appliquée à une méthode de conception développée ici et
accompagnée de la réalisation d’un outil logiciel associé dans l’environnement
Modelica/Dymola. Finalement, l’efficacité de l’approche de conception et de l’outil
développé est illustrée et confirmée avec le dimensionnement d’actionneurs pour le projet
DRESS et pour des commandes de vol primaires (aileron). En plus de résultats intéressant,
cette étude montre la rapidité et l’efficacité de la méthode de conception et de l’outil logiciel
développés pour le dimensionnement de différentes architectures. Ce qui permet d’étudier
plus rapidement un plus grand nombre d’alternatives pour un système donné.
193
CONCLUSION
Une composante essentielle et innovatrice des méthodes d’évaluation proposées ici est
l’élaboration et l’utilisation de lois d’échelle représentatives de phénomènes physiques
dimensionnants, pour lier entre eux les principaux paramètres des modèles (définition,
simulation, etc.). Dans ces travaux, nous avons développé des lois pour des composants
d’actionnements électromécaniques. Par exemple : réducteur de vitesse, moteurs brushless.
Cette méthode est cependant limitée à la représentation de composants simples. Ce qui nous
conduit dans certains cas à combiner des modèles pour aboutir à des composants plus
complexes. De cette façon, le modèle de la vis à rouleaux ou à billes proposé est l’assemblage
des modèles de la vis, de l’écrou et du palier. Cet aspect limite l’application de cette méthode
à un niveau organique de représentation d’architecture et n’est pas applicable à des modèles
de plus haut niveau. Par exemple, cette méthode ne peut pas être appliquée directement pour
réaliser un modèle global d’actionneur électromécanique, dans le cadre d’une évaluation au
niveau aéronef ou véhicule.
194
CONCLUSION
Une dernière limite caractérisant les méthodes proposées se situe au niveau de la prise
en compte des phénomènes parasites dans les modèles de simulation. Par exemple : Les
phénomènes dissipatifs et les élasticités des composants mécaniques. Dans ces travaux, ces
phénomènes sont abordés un haut niveau de détail. De cette façon, la dissipation est
généralement ramenée à un effort ou couple de frottement calculé à partir d’un rendement fixe
représentatif d’un fonctionnement en conditions nominales. Le choix du positionnement des
élasticités a été réalisé sur une approche conservatrice appliquée à un cas simple (inertie-
raideur-inertie) mais pas représentatif de l’ensemble des cas rencontrés en actionnement
électromécanique. Or l’expérience montre que ces phénomènes dissipatifs peuvent avoir une
influence importante sur l’évaluation d’architectures (ici en termes de puissance et de
performances de vitesse en boucle fermée). Il est donc important de pouvoir investiguer ce
point dans de futurs travaux.
Finalement, avant de présenter les pistes pour de futurs travaux. Il est important de
noter que malgré les limites mentionnées précédemment, les méthodes et l’implémentation
proposés ont permis de procéder efficacement à l’indentification et à l’évaluation rapide et
efficace de nouveaux systèmes d’actionnement électromécaniques (exemples des actionneurs
du projet DRESS, de commandes de vol primaires et de VEGA). Un avantage majeur de ces
méthodes et de cette implémentation apparaît aussi dans le cas de la re-conception d’un
système, dans le but de l’adapter à une application similaire, mais avec un facteur d’échelle
différent. Dans ce dernier cas, les efforts mis à contribution pour réaliser les modèles pour la
première conception du système sont intégralement capitalisés et la re-conception fortement
accélérée.
Un axe ayant pour objectif de réaliser des interfaces entre l’approche de conception
préliminaire proposée ici avec les phases de développement antérieures et
postérieures.
195
CONCLUSION
Les incertitudes des évaluations conduites sont aussi un aspect important [Malak
2007], qui pourrait relativiser les résultats des évaluations et qui serait par conséquent
intéressant de prendre en compte dans de futurs travaux.
Toujours au niveau de l’axe d’avancement vertical, il est aussi prévu de créer un pont
entre le niveau de modélisation par paramètres localisés des modèles que nous avons réalisés
ici avec le niveau de modélisation 3D. De cette façon, le passage de la phase de conception
préliminaire à la phase avale de conception spécifique sera significativement facilité et
accéléré.
Finalement, une dernière piste de recherche, qui ne se situerait ni sur l’axe vertical ni
sur l’axe horizontal mentionnés précédemment, mais sur un troisième axe indépendant, serait
de mettre à contribution l’outil logiciel proposé ici pour « nourrir » des méthodologies et des
moyens de simulation se plaçant au niveau avion en aéronautique [Liscouet - Hanke 2008], et
au niveau véhicule dans l’industrie automobile. Ce qui permettrait d’évaluer rapidement et
efficacement des solutions non plus à un niveau local (système) mais à un niveau global
(aéronef, véhicule).
196
CONCLUSION
REFERENCES
[Budinger 2008] M. Budinger, J. Liscouët, S. Orieux et J.-C. Maré. Automated Preliminary Sizing of
Electromechanical Actuator Architectures. ELECTRIMACS 2008 Conference, Québec,
Canada, 2008.
[Liscouet - Hanke 2008] S. Liscouet - Hanke. A model-based methodology for integrated preliminary
sizing and analysis of aircraft power system architectures. Laboratoire de Génie Mécanique de
Toulouse (LGMT) INSA-UPS. Toulouse, Université de Toulouse. Ph. D., 2008.
[Malak 2007] R. J. Malak et C. J. J. Paredis. Using Parameterized Pareto Sets to Model Design
Concepts. ASME International Mechanical Engineering Congress and Exposition
(IMECE2007), Seattle, WA, USA,
[Maré 2010] J.-C. Maré. Towards more electric drives for embded applications: (re)discovering the
advantages of hydraulics. 7th International Conference on Fluid Power (IFK) Aachen,
Germany, 22 - 24 March 2010.
197
CONCLUSION
198
ANNEXES
ANNEXES
i
ANNEXES
ii
ANNEXES
ANNEXE A :
STATISTIQUES POUR LA
FIABILITE
Notation
Nomenclature
iii
ANNEXES
L’objectif de cette annexe est de fournir les notions de base nécessaires à une bonne
compréhension des calculs de fiabilité réalisés dans les chapitres 2 et 4. Pour ce faire, nous
présentons dans un premier temps les principaux modèles mathématiques utilisés en fiabilité.
Dans un deuxième temps, nous illustrons ces modèles avec les deux lois de fiabilité
considérées dans ces travaux : la loi exponentielle et la loi de Weibull.
t
F (t ) = ∫ f (t )dt ( A.1 )
0
t
F (t ) + R(t ) = 1 ⇔ R (t ) = 1 − ∫ f (t )dt ( A.2 )
0
iv
ANNEXES
F (t → t + ∆t ) F (t + ∆t ) − F (t ) R(t ) − R(t + ∆t )
λ(t ) = = = ( A.3 )
R (t )∆t R (t )∆t R (t )∆t
Pour un intervalle de temps tendant vers 0, le taux de défaillance est dit instantané. A
partir de l’équation précédente, on obtient ainsi l’expression du taux de défaillance instantané
suivant :
R(t ) − R(t + ∆t )
z (t ) = lim
dR(t )
=− ( A.4 )
∆t →0 R (t ) R (t )dt
Le début de vie ou rodage. Durant cette phase de vie, les composants marqués par un
défaut de fabrication sont progressivement écartés. Ce qui se traduit par un taux de
défaillance décroissant.
La vie utile. Une fois tous les dispositifs marqués par un défaut de fabrication écartés,
les dispositifs « sains » opèrent avec un taux de défaillance relativement constant
(légèrement croissant pour les dispositifs mécaniques).
Temps (h)
Figure A.1 : Variation du taux de défaillance en fonction du temps pour les trois principales phases de
vie d’un dispositif.
La moyenne des temps de bon fonctionnement, aussi appelé temps moyen entre panne
MTTF (en Anglais : « Mean Time To Failure »), est la moyenne arithmétique du temps entre
pannes d'un système. Par conséquent, dans le cas d’un système non réparable, cette grandeur
s’exprime en fonction de la fonction densité de probabilité de défaillance de la manière
suivante :
v
ANNEXES
+∞
MTTF = ∫ tf (t )dt ( A.5 )
0
dF (t ) dR(t )
f (t ) = =− ( A.6 )
dt dt
Ce qui permet de réécrire (A.5), pour obtenir le temps moyen entre panne en fonction
de la fonction fiabilité de la manière suivante :
=0
+∞ 647 48 + ∞ +∞
MTTF = − ∫ tdR(t ) = − [tR(t )]0 + ∫ R(t )dt = ∫ R(t )dt
+∞
( A.7 )
0 0 0
Dans le cas d’un dispositif réparable, on parle généralement de Mean Time Between
Failure (MTBF) pour désigner le temps moyen entre panne [Hecht 2004]. Cette grandeur est
la somme du temps moyen d’opération avant l’apparition d’une panne et du temps moyen de
réparation Mean Time To Repair (MTTR) comme l’illustre la Figure A.2.
MTBF
MTTR MTTR
MTTF
Réparation
Opération
Panne Panne
Temps (h)
Figure A.2 : Représentation séquentielle des temps moyens de réparation (MTTR) et des temps
moyens entre panne (MTTF et MTBF), pour un dispositif réparable.
En général, le temps moyen de réparation est très petit devant le temps moyen entre
panne MTTF. Ce qui conduit à l’égalité suivante :
vi
ANNEXES
t
dR (t )
= − z (t )dt ⇒ ln R(t ) = − ∫ z (t )dt + constante ( A.9 )
R(t ) 0
t
− ∫ z ( t ).dt
R(t ) = e 0
( A.11 )
dF (t ) d (1 − R (t )) dR(t )
t
− ∫ z ( t ).dt
f (t ) = = =− ⇒ f (t ) = z (t )e 0
( A.12 )
dt dt dt
z (t ) = λ = constante ( A.13 )
vii
ANNEXES
t t
− ∫ z ( t ).dt − ∫ z ( t ).dt
f (t ) = z (t )e 0
= λe − λt
, R(t ) = e 0
= e −λt et F (t ) = 1 − e − λt ( A.14 )
+∞ +∞ +∞
e − λt 1 − λt
MTTF = ∫ R(t )dt = ∫ e dt = = ( A.15 )
0 0 − λ 0 λ
0,1 0,1
0,04 0,04
λ 5λ
λ
0,02 0,02
0 0
0 0,1 0,2 0,3 0,4 0,5 0 0,1 0,2 0,3 0,4 0,5
λt λt
Probabilité de défaillance
(1)
1,2
1
10λ
0,8
5λ
0,6
0,4
t =1/λ=MTTF
0,2 λ
0
0 0,2 0,4 0,6 0,8 1 1,2
λt
viii
ANNEXES
b
b −1 t
bt −
f (t ) = e θ
( A.16 )
θ θ
où θ est le temps caractéristique avant défaillance (F(θ) = 0,63) et b est le paramètre de forme
de la distribution.
b b
b −1 t t
bt − −
1
z (t ) = , R(t ) = e θ
, F (t ) = 1 − e θ
et MTTF = θΓ + 1 ( A.17 )
θ θ θ
+∞
Γ(n ) = ∫ e − x x n −1dx ( A.18 )
0
ix
ANNEXES
Probabilité de défaillance
(1)
1,2
0,8 b = 0,5 b =1
0,6
b = 1,5
b = 3,4
0,4
0,2 t =θ
0
0 0,5 1 1,5 2
t /θ
Références
[Hecht 2004] H. Hecht. Systems reliability and failure prevention, Boston, Artech House, 2004.
[Schwob 1969] M. Schwob, G. Peyrache. Traité de fiabilité, Masson & Cie, 1969.
x
ANNEXES
ANNEXE B :
MODES ET TAUX DE
DEFAILLANCE CONSIDERES
POUR L’ETUDE DE FIABILITE
DE L’EXEMPLE
D’ACTIONNEUR DE DRESS.
xi
ANNEXES
Composant Mode de défaillance Taux de Taux de Source des données Environnement des
fonctionnel défaillance défaillance données source
(/HO)* (/HV)**
-5 -6
Vis sans fin général 1,17 10 2,97 10 RIAC, Engrenage, vis Général
-6 -6
grippage 5,83 10 1,48 10 serré/bloqué (50%)
-6 -6
perte de trans. de couple 5,83 10 1,48 10 jeu excessif (50%)
-6 -7
dégradé 2,91 10 7,42 10 serré/bloqué (50%)
-7 -7
Engrenage droit général 5,70 10 1,45 10 RIAC, Roue Militaire (AIF)
-8 -9
grippage 3,71 10 9,44 10 grippé/bloqué (6,5%)
-8 -9
perte de trans. de couple 3,71 10 9,44 10 rupture (6,5%)
-7 -7
dégradé 4,98 10 1,27 10 usé,collé,griffé,déplacé,bruyant (87,3%)
-9 -10
Crémaillère général 1,07 10 2,73 10 Données DRESS
-10 -10
grippage 5,35 10 1,36 10 supposé (50%)
-10 -10
perte de trans. de couple 5,35 10 1,36 10 supposé (50%)
dégradé 0 0 supposé (50%)
-6 -7
Vis à rouleaux général 2,27 10 5,78 10 RIAC, vis à billes, assemblage Sol (GM)
-8 -8
grippage 9,74 10 2,48 10 serré/bloqué (1/3.13%)
-7 -7
perte de trans. de couple 3,95 10 1,01 10 fracture (17,4%)
-6 -7
dégradé 1,87 10 4,77 10 collé, serré, courbé, opération dégradée (82,4%)
-6 -6
Embrayage général 5,06 10 1,29 10 RIAC, Embrayage Aéronautique (AUF)
électromagnétique grippage 0 0 supposé
-6 -7
perte de trans. de couple 1,81 10 4,60 10 pas de mouvement (35,7%)
-6 -7
dégradé 1,82 10 4,64 10 opération dégradée (7,1%)
-6 -7
bloqué ouvert 1,81 10 4,60 10 perte de trans. de couple
-8 -9
bloqué fermé 1,00 10 2,55 10 (car naturellement ouvert)
-7 -7
Réducteur général 5,70 10 1,45 10 RIAC, Roue Militaire (AIF)
-8 -9
de vitesse grippage 3,71 10 9,44 10 grippé/bloqué (6,5%)
-8 -9
perte de trans. de couple 3,71 10 9,44 10 rupture (6,5%)
-7 -7
dégradé 4,98 10 1,27 10 usé,collé,griffé,déplacé,bruyant (87,3%)
xii
ANNEXES
Composant Mode de défaillance Taux de Taux de Source des données Environnement des
fonctionnel défaillance défaillance données source
(/HO)* (/HV)**
-6 -6
Moteur brushless général 8,40 10 1,83 10 Moteur, synchrone Sol, commercial (GF)
-7 -7
grippage 4,80 10 1,04 10 10% de la défaillance générale attribué
au rouelement (57,1%)
-6 -7
perte de trans. de couple 3,61 10 7,86 10 rupture de câble électrique, surcharge
électrique, défaillance transistor
(42,9%)
-6 -6
dégradé 7,92 10 1,72 10 défaillance générale de roulement,
rupture de câble, surcharge électrique,
défaillance transistor (42,9%)
-6 -6
Electronique de général 5,08 10 1,10 10 Convertisseur Général
-5 -6
puissance contaminé 2,00 10 4,35 10 contamination + interrupteurs de puiss.
-7 -7
pas de puiss. en sortie 8,43 10 1,83 10 pas de puiss. en sortie
-6 -7
dégradé 3,96 10 8,61 10 dépassement de la spécification,
surchauffe
-6 -6
Systéme de général 6,50 10 1,66 10 Données DRESS
-6 -6
débrayage grippage 6,50 10 1,66 10 supposé
-6 -7
bloqué ouvert 1,00 10 2,55 10 typique pour équipement électronique
-6 -7
bloqué fermé 1,00 10 2,55 10 typique pour équipement électronique
* HO = Heure d’opération
** HV = Heure de vol
xiii
ANNEXES
xiv
ANNEXES
ANNEXE C :
ARBRES DE DEFAILLANCE DE
L’EXEMPLE DE
L’ACTIONNEUR DE DRESS
xv
ANNEXES
FL1.1 FL2.1
Mouvement libre et incontrôlé Mouvement libre et incontrôlé
de la voie de puissance 1 de la voie de puissance 2
xvi
ANNEXES
FD0
Mode d’opération dégradé
FD1.1 FD2.1
Débrayage de Débrayage de
la voie de puissance 1 la voie de puissance 1
xvii
ANNEXES
xviii
ANNEXES
ANNEXE D :
MODÈLE DE LA VIS SANS FIN
(EN ANGLAIS)
NOMENCLATURE
xix
ANNEXES
INDICES
xx
ANNEXES
The present annexe describes the proposed simulation model and scaling laws
established for the worm gear. Then, the scaling laws are resumed in a table before being
validated by comparison with manufacturer data.
Figure D.1 shows a classical worm-gear configuration with its different geometrical
parameters: dsc the nominal screw diameter, lsc the screw length, p the pitch, λ the lead angle,
φn the pressure angle, lwh the wheel width, dwh,in the internal wheel diameter, dwh,ex the external
nominal wheel diameter and ∆dwh the difference between these two last diameters.
xxi
ANNEXES
mechanical transformation, which returns the axial and radial forces in the rotational axes and
takes into account the losses due to the friction between the screw meshes and gear teeth,
along with the associated thermal dissipation. The axial and radial forces are absorbed by
bearings, whose model has already been described in the chapter 3. A particular effort is put
into modelling the friction as it can affect the operation of the worm-gear significantly. The
proposed model is illustrated in Figure D.2.
Legend
Fr,s
Fa,s Thermal connector
1D rotational power transmission
screw Jsc wheel-bearings 2D thermal power transmission
wheel 3D effort transmission
Ta, θa Jwh
Tb, θb
screw inertia
p,ds,dw,ηd,ηi
screw-bearings wheel inertia
Fa,w
90° rotation-rotation Fr,w
transformation with losses
In a receptor convention, the overall simulation model for the worm-gear is described
by the following equations:
( )
k Ta − J sc &θ& a + Tb − J wh &x&b + T f = 0
( D.1 )
θ a = kθ b
where Ta is the (input) torque at the screw, Jsc is the screw moment of inertia, θa is the screw
angular position, k is the speed reduction ratio, Tb is the (output) torque at the wheel, Jwh is the
wheel moment of inertia, Tf is the frictional torque representative of the power transmission
losses and θb is the wheel angular position.
The wheel and screw torques are due to tangential forces, which are also radial and
axial forces on the rotational axes:
where Fsc,r is the radial force on the screw axis, Fwh,ax is the axial force on the wheel axis, Fwh,r
is the radial force on the wheel axis, Fsc,ax is the axial force on the screw axis and dwh is the
nominal wheel diameter. The axial and radial forces must be compensated by bearings in
order to avoid any undesired relative motion between the wheel and the screw. Consequently,
a tri-dimensional force connector implemented into the worm-gear simulation model, in order
xxii
ANNEXES
to interface with the bearing ones. The losses due to the friction at the meshing are assimilated
to a four-quadrant rotational friction.
( ) ( )
k (η d − 1) Ta − J sc &θ& a , if the screw is driving : Ta − J sc &θ& a θ& a > 0
( ) ( )
T f = ∑ − (η i − 1) Tb − J wh &θ&b , if the wheel is driving : Tb − J wh &θ&b θ& b < 0 ( D.3 )
0, if braking or no power transmited
The friction model requires the direct (power flows from screw to wheel) and indirect
(power flows from wheel to screw) efficiencies, in order to calculate the friction torque. The
ηd,wg direct and ηi,wg indirect efficiencies are given by the frictional losses occurring at the
meshing between the screw mesh and wheel teeth. Thus, these frictions depend on the worm-
gear geometry (tree dimensional) and the µ friction coefficient [Budynas 2007]:
cos φ n − µ tan λ
ηd , wg = cos φ + µ cot λ
n
( D.4 )
η = cos φ n − µ cot λ
cos φ n + µ tan λ
i , wg
Typically, the friction coefficient is assumed constant and equal to 0.05. Actually, it
depends on the vs translational meshing velocity and ranges between 0.145 and 0.016 for a
classical iron-bronze design. Near the zero speed this friction coefficient is significantly
higher. This particularity causes the lifetime of the worm-gear to be potentially sensitive to
the number of “starts” occurring during the mission cycle. It is therefore relevant to take into
account the variation of the friction coefficient as a function of the meshing velocity, e.g. by
applying standard calculation methods [Hohn 1998]. In the absence of compliance and
backlash there is no relative parasitic motion between the screw and the wheel. Therefore, the
pitch of the screw is a function of the transmission ratio and wheel diameter:
p = πd wh k (D.5 )
The losses of the worm-gear are dominated by the friction between the meshes. These losses
can be important and must be dissipated to avoid an overheating of the mechanism. Typically,
xxiii
ANNEXES
l * = ∆d * = T *1 3
wh wh wh
*
d sc = l wh
*
( D.7 )
*
l sc = d wh ,ext
*
Prior to developing more specific scaling laws for a variable internal diameter of the
wheel, it is interesting to consider a homothetic scaling of all the dimensions. Indeed, the
latter results can be compared with manufacturer data and thus give confidence in developing
more specific scaling laws.
l wg = Twh
* *1 3
( D.8 )
xxiv
ANNEXES
M wg * = Twh *
* ( D.9 )
J wg = Twh
*5 3
where Mwg is the mass of the worm-gear box and Jwg the screw moment of inertia.
S wh = l wh ∆d wh
* * *
( D.10 )
In the same way the Ftan,wh tangential force applied to the wheel is given by the
transmitted torque as follows:
*
Twh
Ftan, wh =
*
*
( D.11 )
d wh ,ext
For a constant mechanical stress in the wheel, (D.10) and (D.11) yields
*2 / 3
d wh , ext .∆d wh = Twh
* *
( D.12 )
Substituting the result of (D.12) into (D.7) gives the remaining geometrical parameters
of the worm-gear (lwh, ds, ls). The mass and inertia of the wheel are then approximated as
π
(
m wh = ρ m, wh l wh 4 d wh ,ext − d wh ,in
2 2
)
( D.13 )
J = ρ
wh m , wh l wh
π
12
(
d wh ,ext − d wh ,in
3 3
)
where mwheel is the wheel mass and ρm,wheel is the wheel mass density, which can be calculated
from the reference data geometrically. Typically, the wheel is made of bronze, with a density
ranges between 8 400 and 9 200 kg/m3. In the same way, the mass and inertia of the screw are
obtained from its geometry and its mass density. Typically, the screw is made of iron, with a
density of 7 860 kg/m3.
xxv
ANNEXES
The speed limitation of the worm-gear is given by its bearings, as described in the
chapter 3.
β +1
*−
Twh , nom, de − rated = Twh , nom ω sc α ( D.15 )
where Twh,nom,de-rated is the de-rated nominal wheel torque, ωsc* is the evolution of the nominal
screw speed and α and β are respectively experimental torque and speed factors (e.g., α = 2
and β = -1,33).
xxvi
ANNEXES
Table D.1: Proposed scaling laws established for the worm gears.
Mass kg
M wg = Twh
* * (D.9)
mwh = ρ m, wh l wh
π
4
(
d wh ,ext − d wh ,in
2 2
) (D.13)
π 2
msc = ρ m , sc l sc d sc -
4
Simulation parameters
Moment of
Inertia
kg.m2
J wg = Twh
* *5 3 (D.9)
J wh = ρ m,wh .lwh .
π
12
(
d wh ,ext − d wh ,in
3 3
) (D.13)
π -
J sc = ρ m ,sc l sc
3
d sc
12
Direct - cos φ n − µ (vs ) tan λ
efficiency ηd , wg (t ) =
cos φ n + µ (vs ) cot λ
(D.4)
Operational limit
parameters
Maximum N or Nm
Twh ,max = Twh ,nom (D.14)
* *
torque
xxvii
ANNEXES
1,E+00
Scaling laws Prod. range =
References + 26%; - 11%
DYNASET Gear diameter
DYNASET Worm length
DYNASET Gear width
DYNASET Worm diameter
scal worm length
Dimensions (m)
Prod. range =
+ 2%; - 0%
1,E-01
Figure 3: Worm-gear wheel dimensions as a function of the nominal output (wheel) torque.
1,E+04
Scaling laws Prod. range=
CAVEX CUW01 K88 series + 68%; - 32%
DYNABOX
References
1,E+03
Mass (kg)
1,E+02
Figure 4: Worm-gear box mass (incl. housing) as a function of the nominal output (wheel) torque.
xxviii
ANNEXES
1,E+01
Scaling laws Prod. range =
References
1,E+00 DYNABOX + 0%; - 59%
CAVEX K88 CUW01 series
ref cavex inertia
1,E-01
Inertia (kg.m²)
1,E-02
1,E-03
Input speed = 2000 rpm
1,E-04 DYNABOX
reduction ratio = 10.25:1
Références
[ANSI/AGMA 1983] ANSI/AGMA. Specification for worm gearing. Metric units: 28, 30/09/83.
[Girard 2003] Dynabox XL - Right angle servo gearheads for heavy duty applications, G. T.-P.
gearheads, 2003. http://www.girard-transmissions.com/
[Girard 2008] Dynabox - Right angle servo gearheads, G. T.-P. gearheads, 2008. http://www.girard-
transmissions.com/
[Hohn 1998] B.-R. Hohn and K. Streingrover. DIN 3996: Calculation of load capacity of cylindrical
worm gear pairs with shaft angle Σ = 90°. D. I. f. N. e.V., Beuth Verlag, 1998.
xxix
ANNEXES
xxx
Titre : Conception préliminaire des actionneurs électromécaniques - Approche hybride,
directe/inverse
Résumé : L’objectif de cette thèse est de proposer des méthodes innovantes de conception
préliminaire d’actionneurs électromécaniques (EMA) et de les implémenter sous la forme
d’outils logiciels rapidement disponibles pour les industriels. Cet objectif répond à une
demande forte de l’industrie, en particulier en aéronautique dans le cadre du développement
d’avions plus électriques. Dans un premier temps, cette thèse propose une méthode hybride
(montante et descendante) de recherche systématique d’architectures solutions et de sélection
vis-à-vis des exigences du cahier des charges et de l’état de l’art technologique. Dans un
deuxième temps, des méthodes d’évaluation d’architectures en termes de puissance,
d’intégration (enveloppe géométrique et masse), de fiabilité et de performances en boucle
fermée sont proposées. L’implémentation de ces méthodes se base sur la modélisation
acausale et la combinaison de simulations numériques inverses et directes. Des lois d’échelle,
représentatives des phénomènes physiques dimensionnants, sont établies pour réduire la
complexité d’utilisation des modèles et l'intervention d’experts de domaine dans les phases
préliminaires. Les méthodes proposées et leur implémentation dans l'environnement de
simulation Modelica/Dymola ont été appliquées avec succès aux exemples d’actionneurs
électromécaniques d’orientation d’un train avant d’atterrissage, de commandes de vol
primaires et de contrôle de la poussée vectorielle du premier étage de propulsion d’un lanceur
spatial. De cette façon, la durée de la phase d’évaluation d’architectures a pu passer d’un
ordre de grandeur en jour à un ordre de grandeur en heure.
Discipline : Génie Mécanique
Mots-clés : Actionneur électromécanique (EMA), avion plus électrique, conception
préliminaire, conception orientée simulation, fiabilité, lois d’échelle, méthode
montante/descendante, Modelica, modélisation acausale, simulation inverse.
Jonathan LISCOUËT