Les conceptions récentes de turbomachines ont tendance à évoluer vers des liaisons entre étages de plus en plus souples et des niveaux d'amortissement faibles, donnant lieu à des configurations où les modes sont susceptibles de présenter des niveaux de couplages inter-étages forts. En général, les ensembles disques aubes multi-étagés n'ont aucune propriété de symétrie cyclique d'ensemble et l'analyse doit porter sur un modèle de la structure complète donnant lieu à des calculs très coûteux. Pour palier ce problème, une méthode récente appelée symétrie cyclique multi-étages peut être utilisée pour réduire le coût des calculs des rotors composés de plusieurs étages, même lorsque les étages ont un nombre différent de secteurs. Cette approche profite de la symétrie cyclique inhérente à chaque étage et utilise une hypothèse spécifique qui aboutit à des sous-problèmes découplés pour chaque ordre de Fourier spatial. La méthodologie proposée vise à étudier l'effet des incertitudes sur le comportement dynamique des rotors en utilisant l'approche de symétrie cyclique multi-étages et l'expansion en Chaos Polynomial. Les incertitudes peuvent découler de l'usure des aubes, des changements de température ou des tolérances de fabrication. En première approche, seules les incertitudes provenant de l'usure uniforme de l'ensemble des aubes sont étudiées. Celles-ci peuvent être modélisées en considérant une variation globale des propriétés du matériau de l'ensemble des aubes d'un étage particulier. L'approche de symétrie cyclique multi-étages peut alors être utilisée car l'hypothèse de secteurs identiques est respectée. La positivité des matrices aléatoires concernées est assurée par l'utilisation d'une loi gamma très adaptée à la physique du problème impliquant le choix des polynômes de Laguerre comme base pour le chaos polynomial. Dans un premier temps des exemples numériques représentatifs de différents types de turbomachines sont introduits dans le but d'évaluer la robustesse de la méthode de symétrie cyclique multi-étages. Ensuite, les résultats de l'analyse modale aléatoire et de la réponse aléatoire obtenus par le chaos polynomial sont validés par comparaison avec des simulations de Monte-Carlo. En plus des résultats classiquement rencontrés pour les fréquences et réponses forcées, les incertitudes considérées mettent en évidence des variations sur les déformées modales qui évoluent entre différentes familles de modes dans les zones de forte densité modale. Ces variations entraînent des modifications sensibles sur la dynamique globale de la structure analysée et doivent être considérées dans le cadre des conceptions robustes.

Cette thèse a pour objectif de mener une optimisation des corrections de forme des engrenages cylindriques, droits et hélicoïdaux. Le travail se décompose en quatre parties principales. Dans la première partie, on présente un état de l'art sur les différents types de corrections de forme proposées dans la littérature. Une analyse des travaux d'optimisation, menés jusqu'à présent, est conduite. La deuxième partie est focalisée sur une approche déterministe visant à cerner l'influence des corrections de dentures sur les principaux critères de performance. Dans ce contexte, on propose un développement analytique qui caractérise les fluctuations d'erreur de transmission quasi-statique permettant d'obtenir des relations approchées originales. En présence de plusieurs paramètres de corrections, un algorithme génétique est utilisé afin d'identifier, en un temps réduit, les solutions optimales. Nous proposons, en troisième partie, une étude probabiliste pour caractériser les corrections robustes. Ainsi, on définit une fonction objectif de robustesse faisant intervenir des paramètres statistiques. Après une étape de validation, l'estimation de ces paramètres est effectuée en utilisant les formules de quadrature de Gauss. Plusieurs études paramétriques sont ensuite menées et qui reflètent entre autre l'influence des classes de qualité, la forme de la correction, etc. Enfin, on a conduit une optimisation multicritère en utilisant un algorithme d'optimisation spécifique : " NSGA-II ".

Ce travail vise à prédire le comportement dynamique de transmissions par engrenages supportées par des paliers hydrodynamiques, semblables à celles utilisées dans le domaine de la propulsion navale. Un modèle couplé de transmission mécanique est présenté ; il traite la plupart des interactions possibles entre les engrenages, les arbres et les paliers. Un élément d'engrenages à dentures larges, avec raideur non linéaire et dépendante du temps, est combiné à des éléments finis de poutre alors que la contribution des paliers est introduite par la résolution directe de l'équation de REYNOLDS et d'un model thermique simplifié. A cause des jeux de fonctionnement des paliers, une attention particulière a été portée à la définition de l'état de référence et des degrés de liberté. La réponse du système est calculée par la combinaison d'un schéma d'intégration temporelle, d'une méthode de NEWTON-RAPHSON et d'un algorithme de contact normal unilatéral de manière à ce que les conditions de contact aux paliers et entre les dents soient simultanément traitées. Les résultats de simulation sont comparés à des mesures réalisées sur un banc d'essai de haute précision possédant un train simple étage en dentures droites et hélicoïdales, supporté par des paliers lisses de type hydrodynamique. Les comparaisons fournissent un bon accord et aboutissent à la validation du modèle de couplage tant sur les aspects globaux que locaux. Bon nombre de résultats sont également présentés, qui montrent que certains paramètres, bien souvent ignorés dans les modèles de la littérature, tels que la position de la zone d'alimentation des paliers, la température de lubrifiant dans les paliers, les accouplements avec l'environnement extérieur peuvent fortement modifier l'équilibre statique et le comportement dynamique du système, jusqu'aux contacts entre les dents.

De nombreuses hypothèses sont classiquement utilisées pour décrire le comportement du fluide dans un contact lubrifié : film continu, viscosité constante dans l'épaisseur, film mince, fluide newtonien... Or, certaines s'avèrent erronées dès lors que l'on s'intéresse aux contacts Elasto- HydroDynamiques fortement glissants ou à la répartition du lubrifiant en sortie de contact. Une approche numérique originale, basée sur un retour aux équations de la mécanique des fluides générale et prenant en compte le couplage fluide/solide et les effets thermiques sont proposés ici dans le but d'apporter des éléments physiques supplémentaires aux modélisations usuelles. Dans un premier temps, l'influence des effets thermiques sur l'évolution du frottement dans les contacts Thermo-EHD est mise en évidence. La présence d'un minimum de frottement pour le cas du glissement pur est expliquée par l'analyse des transferts thermiques entre le lubrifiant et les solides. L'origine des modifications locales d'épaisseur de film observées et l'existence même d'une épaisseur de film lubrifiant pour les cas de vitesse d'entraînement nulle sont alors reliées à la présence d'un fort gradient de viscosité dans l'épaisseur de film. Une comparaison qualitative avec des éléments expérimentaux de la littérature est réalisée, validant les tendances obtenues. Dans un second temps, l'écoulement à surface libre du lubrifiant en périphérie du contact est étudié expérimentalement puis numériquement par une méthode à interface diffuse. Le rôle des effets capillaires est analysé et les résultats numériques confrontés à des résultats issus de la littérature. Un bon accord est obtenu tant qualitativement que quantitativement. Validé par l'étude numérique diphasique (air/lubrifiant) réalisée, un modèle analytique simplifié est alors développé, prédisant une loi de répartition du lubrifiant en sortie de contact. La zone de sortie des contacts EHD est ensuite traitée par un modèle de cavitation vaporeuse et la prise en compte nécessaire de l'air environnant est discutée. Enfin, une première modélisation tridimensionnelle de l'écoulement à surface libre du lubrifiant autour d'un contact ponctuel est réalisée mettant en avant l'influence des effets capillaires et la faisabilité d'une telle approche.

Les ensembles turbo-alternateurs des centrales électriques sont de grandes machines tournantes de plus de 50 mètres de long et de plusieurs centaines de tonnes. Lors du fonctionnement normal d'une telle machine, une probabilité non nulle existe d'un détachement accidentel d'une aube. Dans une telle situation, un balourd important est généré et du contact apparaît entre les parties tournantes et non tournantes de la machine. Il est alors capital de pouvoir simuler efficacement la dynamique de ce type d'évènement faisant intervenir de fortes non linéarités dans le système. Cette thèse a été réalisée dans le cadre du projet ANR (Agence Nationale de la Recherche) IRINA (SImulation et maîtRise des rIsques en coNception des mAchines tournantes) et en particulier entre le LaMCoS (LAboratoire de Mécanique des Contacts et des Structures) de l'INSA de Lyon et le département AMA (Analyses Mécaniques et Acoustiques) d'EDF R et D à Clamart. Elle a pour objectif de mettre au point une technique rapide de simulation du comportement des lignes d'arbres de machines tournantes en cas de présence de non linéarité de type contact entre rotor et stator. Pour atteindre cet objectif, une double démarche a été mise en place. La première consiste à mettre au point des modèles simplifiés afin de réduire le nombre de degrés de liberté du problème. De surcroît, une technique de réduction de modèle adaptée au cas de non linéarité localisée est utilisée afin de réduire encore plus la taille du système à résoudre. La seconde démarche consiste à mettre au point une technique de résolution rapide du système réduit afin d'obtenir la solution encore plus rapidement. Pour cela, au lieu d'utiliser les traditionnelles techniques d'intégration temporelle directe, c'est la méthode de la balance harmonique qui est mise à profit. Cette technique permet d'obtenir directement la réponse stabilisée du système grâce à une résolution des équations dans le domaine fréquentiel. Dans ce cadre, une maquette numérique a été mise au point mettant en oeuvre les fonctionnalités citées. Cette dernière permet de reproduire les phénomènes physiques périodiques ainsi que quasi-périodiques et de déterminer leur stabilité. Des études paramétriques sur des exemples de problèmes de contact rotor-stator viennent illustrer cette démarche. Enfin, une application sur un cas industriel de groupe turbo alternateur EDF est présentée.

Cette thèse est consacrée è la modélisation et a la simulation de l'ensemble PISTON /SEGMENT/ CHEMISE, d'un moteur thermique. Dans un moteur, le piston muni de segments est en mouvement relatif dans la chemise. Ce système est lubrifié par projection d'huile, les segments et la surface de la chemise "contrôlent" la quantité de lubrifiant restant sur la paroi. Pour modéliser ce phénomène, le modèle conservatif parabolique-hyperbolique d'Elrod-Adams (P − θ), prenant en compte la cavitation (présence de gaz dans le fluide, aspect multiphasique de l'écoulement), est couplé a un modèle de contact, celui de Greenwood-Tripp (basé sur une approche statistique de la surface). En se restreignant a l'écoulement et au contact, nous négligeons beaucoup de phénomènes physiques pouvant intervenir comme le grippage ou les effets thermiques. Dans une première partie, les phénomènes mécaniques sont étudiés, de la réalisation de la surface des chemises è la modélisation de la physique de la lubrification du contact segment/chemise. Cette étude expose les simplifications effectuées et le choix des conditions aux limites pour la prise en compte de la cavitation en tant que problème a frontière libre, notamment au niveau du débit d'entrée pour simuler un fonctionnement "normal". Dans une deuxième partie, la mise en oeuvre et la modification de l'algorithme "classique" pour des surfaces mesurées sont décrites en détail. Nous présentons les performances simulées de différentes surfaces, et nous les comparons è des essais effectués par Renault. Dans une troisième partie, des modifications sont introduites au modèle (P − θ) utilisé précédemment afin de rendre l'écoulement cohérent avec le modèle de Navier-Stokes. Bien que ce nouveau problème puisse être posé en dimension 2, l'étude mathématique ne concerne que la dimension 1. On se ramène a un système dynamique dont les inconnues sont les frontières libres. L'unicité de la solution impose l'introduction d'un paramètre artificiel. Le comportement local et/ou global en temps est déterminé par les paramètres géométriques et les conditions aux limites. L'introduction de deux paramètres supplémentaires permet de démontrer l'existence et l'unicité de la solution. Dans une dernière partie, nous comparons le modèle de cavitation proposé avec le modèle P − θ.

Un modèle de trains planétaires est proposé afin de tenir compte de l'influence d'erreurs de position et d'excentricités en lien avec d'éventuels montages 'flottants' sur le comportement dynamique d'une transmission. La formulation théorique repose sur le formalisme des torseurs de déplacements infinitésimaux pour simuler à la fois les erreurs géométriques et les degrés de liberté du modèle. Une des propriétés principales de cette approche est que la géométrie des engrènements et les excitations correspondantes sont couplées aux degrés de liberté, conduisant ainsi à des excitations complexes présentant des modulations d'amplitude et de phase. Les résultats de simulation sont comparés avec des mesures sur banc d'essai et un très bon accord est obtenu en terme de partage de charge entre les satellites, validant ainsi le modèle de contact développé. Enfin, des résultats d'études paramétriques portant sur le rôle de certaines erreurs ainsi que sur l'apport éventuel de solaire et/ou satellites flottants dans des applications grandes vitesses concluent ce travail de thèse.

L'allègement des structures imposé par les réductions de coût se traduit par des structures de plus en plus souples qui les rendent de plus en plus sensibles aux vibrations. Le contrôle des vibrations devient donc un enjeu majeur dans de nombreuses applications industrielles et les limites des matériaux imposent maintenant un recours au contrôle actif de plus en plus fréquent. L'évolution des structures au cours du temps (viellisement, conditions aux limites, architecture, ...) pose le problème de la robustesse du contrôle. Par ailleurs, l'actionnement de plus en plus présent dans le domaine mécanique constitue à la fois une source supplémentaire de vibrations, mais aussi de contrôle et d'évolution d'architecture des structures. La thèse s'intéresse au contrôle actif autoadaptatif des vibrations permettant de maintenir automatiquement la performance et la stabilité des structures évolutives. Il s'agit donc de s'affranchir de la connaissance des causes et des informations sur les évolutions. La méthode proposée s'appuie sur un développement modal permettant de limiter le nombre de composants de contrôle et de cibler les modes à contrôler en limitant l'énergie de contrôle. Ainsi, il est nécessaire de reconstruire les caractéristiques du modèle modal indispensables pour réactualiser le contrôle en figeant seulement une structure de modèle. S'affranchissant à la fois des causes d'évolution de la structure et utilisant seulement une structure de modèle, la méthode est généralisable à toute application en mécanique des structures. La méthode proposée, basée sur l'utilisation d'un identificateur exploitant à la fois excitation et réponse de la structure, prend en compte les limites imposées par le contrôleur. Le modèle constitue le lien qui doit être établi entre identificateur et contrôle pour permettre la réactualisation. Par ailleurs, un compromis entre l'objectif d'atténuation des vibrations et les performances de l'identification est alors nécessaire du fait du couplage identification/contrôle apparaissant dans la boucle fermée. Ce compromis est également conditionné par le matériel utilisé. La méthode proposée est exploitée sur une structure discrète mettant en évidence une inversion de formes modales au cours de son évolution qui déstabilise un contrôle figé. Le choix opéré pour répondre aux différents compromis cités ci dessus a conduit à l'utilisation d'un contrôleur classique (LQG) et un identificateur basé sur la méthode des sous-espaces (N4SID). Cette application sur une structure simple a permis de caractériser un certain nombre de limites physiques : la bande passante, densité modale, vitesse d'évolution, Le contrôle modal autoadaptatif proposé s'avère robuste en performance et efficace lorsque la réactualisation est systématique. Une variante conditionnelle, toujours basée sur l'analyse de la réponse de la structure, est enfin proposée pour optimiser le processus de réactualisation afin de suivre plus efficacement les évolutions.

La corrosion sous contrainte (CSC) est un phénomène synergique d'endommagement qui résulte d'un processus de corrosion (dissolution, adsorption) et d'une rupture mécanique (fissuration). Les mécanismes de couplage mécano-électrochimique en pointe de fissure nécessaires à la compréhension du phénomène sont encore mal connus puisqu'ils dépendent du système d'étude (métal/milieu agressif) et font intervenir de nombreux facteurs mécaniques et électrochimiques. Dans cette thèse, nous nous proposons d'étudier les interactions réciproques entre la dissolution et l'état de contrainte mécanique en pointe de fissure (facteurs d'intensité des contraintes) pour le cas du Zircaloy-4 en milieu aqueux halogéné. Les éprouvettes sont d'abord pré-fissurées par fatigue à l'air avec la technique du Load-Shedding, ce qui permet alors de maîtriser le facteur d'intensité des contraintes résiduel de pré-fissuration par fatigue. Ensuite, une pré-oxydation thermique est réalisée pour produire une couche de film passif sur les surfaces de l'éprouvette. Les réactions électrochimiques sont alors concentrées en pointe de fissure qui induit également une concentration des effets mécaniques. Des techniques de corrélation d'images sont développées dans le but d'identifier les facteurs d'intensité des contraintes et d'estimer la longueur de fissure en temps réel. Ainsi des essais originaux de CSC, permettant d'imposer les facteurs d'intensité en pointe de fissure, sont conduits et les résultats montrent alors clairement les effets synergiques entre l'évolution des mesures mécaniques et des courants de dissolution. D'autre part, l'existence d'un facteur d'intensité des contraintes seuil de propagation en CSC, présentant une forte dépendance à l'histoire de chargement et à la variation de l'état de contrainte locale en pointe de fissure (variation des facteurs d'intensité des contraintes), est mise en évidence. Ceci montre alors que la plasticité doit évoluer pour que la dissolution se produise. Ainsi, le temps caractéristique d'arriver de nouvelles dislocations en pointe de fissure ne doit pas dépasser la durée caractéristique des réactions de passivation. Enfin sur la base des résultats expérimentaux, un modèle de loi de propagation phénoménologique de fissure en CSC est proposé et ses paramètres sont identifiés et validés à partir de différents essais de CSC.

L'intégrité des portées d'étanchéité des robinets des centrales REP nécessite l'emploi de revêtements durs comme par exemple le Norem02, base Fe, déposé par soudage PTA (Plasma Transferred Arc). La prédiction de leur tenue en service nécessite la connaissance de leur état mécanique en fin de fabrication. Pour cela, la simulation numérique du procédé de dépôt doit être effectuée. La microstructure du Norem02 est essentiellement constituée d'une matrice austénitique parsemée d'ilôts de ferrite et de carbures. Son évolution et ses éventuelles transformations allotropiques sont étudiées durant des cycles thermiques représentatifs de ceux subis par le matériau au cours du dépôt, amenant à la conclusion que l'évolution métallurgique de l'alliage n'a que très peu d'influence sur ses propriétés mécaniques. En s'appuyant sur des essais jusqu'à 1000°C de dilatométrie libre, de traction isotherme, ainsi que de traction-compression, l'évolution des caractéristiques thermomécaniques est déterminée en fonction de la température et une loi de comportement élasto-plastique à écrouissage cinématique non linéaire est identifiée. Des mesures expérimentales de températures, déplacements, déformations et contraintes résiduelles lors d'essais de dépôt PTA sur des maquettes représentatives de la géométrie réelle d'un robinet permettent l'identification d'un chargement thermique équivalent. Enfin, des simulations numériques sont conduites en utilisant le code éléments finis développé au sein d'EDF R et D Code_Aster. Elles permettent la validation du caractère prédictif des modèles développés par comparaison entre les résultats mécaniques simulés et expérimentaux. Une excellente corrélation est ainsi observée entre les températures mesurées et simulées tout au long du procédé. De la même façon, à l'issue du refroidissement, les distorsions, déformations et contraintes résiduelles expérimentales et simulées sont en adéquation, validant ainsi la stratégie de modélisation proposée dans ces travaux.