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.

L'identification et la caractérisation des sources acoustiques restent encore aujourd'hui deux sujets d'importance pour les industriels qui ont besoin de techniques permettant d'identifier des sources acoustiques ou vibratoires sur des surfaces complexes dans un environnement acoustique non contrôlé. La thèse que nous présentons ici s'inscrit dans ce cadre. Nous y étudions les qualités d'une méthode inverse, appelée iPTF (pour inverse Patch Transfer Functions, pour la résolution de ce problème particulier. Nous consacrerons le premier chapitre de notre étude à la synthèse bibliographique des méthodes les plus pertinentes permettant de résoudre des problèmes similaires. Nous soulignerons également les difficultés de ces méthodes pouvant être liées à leurs applications pratiques ou à leurs fondements théoriques. Dans un second chapitre, nous présenterons la méthode iPTF à partir de sa formulation directe, c'est-à-dire de la source vers le bruit rayonné. Celle-ci est une approche par sous-structuration de domaines permettant l'étude des problèmes vibro-acoustiques en basses et moyennes fréquences. Nous montrerons particulièrement la façon dont l'association des deux formulations directe et indirecte permet de définir une méthode pouvant conduire jusqu'à l'identification des trois champs de vitesses, de pressions et d'intensités sur la surface de l'objet source. Notre troisième chapitre sera consacré à la présentation des premiers résultats d'identification dans un cas d'application numérique simple. Nous effectuerons, dans les chapitres quatre et cinq, une étude des principaux paramètres conditionnant les résultats donnés par la méthode. Le premier de ces deux chapitres présentera la mise en place d'un code de calcul permettant de résoudre rapidement le problème de rayonnement d'une structure simple. La méthode ainsi définie sera utilisée dans le chapitre cinq afin de générer de nombreux champs rayonnés présentant des caractéristiques différentes de façon à étudier la stabilité de la méthode iPTF face à la variation de différents paramètres. Un sixième chapitre présentera une étude approfondie faite sur les ondes évanescentes qui prennent une part non négligeable dans le champ rayonné par les structures. Ces ondes évanescentes, par définition, ne peuvent pas être mesurées au delà du champ proche, ce qui peut être la cause d'une part des défauts d'identification rencontrés lors de l'application de notre méthode. L'étude faite ici aura pour but de déterminer l'importance de ces ondes évanescentes dans le champ rayonné. Nous présenterons enfin, dans un dernier chapitre, les résultats de mesures expérimentales réalisées.

Les systèmes de transmission par engrenages sont largement utilisés pour transmettre de la puissance et adapter les vitesses de rotation entre organes moteurs et récepteurs. Dans ce contexte, les engrenages sont fréquemment les organes parmi les plus sensibles de la chaîne cinématique et peuvent être soumis à un grand nombre d'avaries (fatigue de contact, fatigue de flexion, usure...etc.) apparaissant lors du fonctionnement et dont les causes sont multiples. L'objectif de ce travail est, d'une part, de simuler l'usure abrasive et la fatigue de contact conduisant à de l'écaillage (pitting) et, d'autre part, d'analyser les interactions entre ces avaries et le comportement statique et dynamique de transmission par engrenages. A cette fin, un modèle dynamique tridimensionnel d'engrenages de fortes largeurs est couplé à des modèles d'usure et d'avaries de contact. L'usure est simulée en s'appuyant sur le modèle d'Archard modifié afin de tenir compte de l'influence du régime de lubrification. Les usures obtenues après un certain nombre de cycles de chargement sont considérées comme des écarts initiaux additionnels par rapport à la géométrie idéale du flanc de denture. Les phénomènes de fatigue de contact par pitting sont analysés en deux étapes; a) une période d'initiation de fissure simulée en s'appuyant sur plusieurs critères de fatigue multiaxiaux et b) une phase de propagation de fissure traitée par la mécanique linéaire élastique de la rupture. Les sollicitations dynamiques fournies par le modèle dynamique d'engrenages sont utilisées comme données d'entrée pour la simulation des périodes d'initiation puis de propagation. Un grand nombre d'exemples d'application sont présentés et les interactions entre comportement dynamique, usure et fatigue sur des engrenages cylindriques sont analysées.

La conception des engrenages spiro-coniques reste encore très complexe de nos jours car la géométrie des dentures, et donc les performances cinématiques, découle du mode de fabrication de ce type d'engrenage. Le taillage est lié à deux constructeurs principaux : Gleason et Klingelnberg. De nombreux paramètres de réglage des machines influencent directement les surfaces de denture, leur optimisation n'est donc pas intuitive. Avec les progrès réalisés cette dernière décennie par les machines d'usinage à commande numérique et la FAO (Fabrication Assistée par Ordinateur), il devient possible de fabriquer des engrenages spiro-coniques de bonne qualité sur une machine 5 axes. Un modèle numérique a été développé pour générer une géométrie simplifiée de type Gleason, usinée par la suite avec une machine 5 axes. Une étude de métrologie, permettant de comparer les dents usinées avec les modèles CAO, a ensuite été réalisée pour prouver que l'usinage par une machine 5 axes peut être une alternative aux méthodes de taillage classiques. De nouveaux types de géométrie peuvent donc être proposés, qui ne pouvaient pas être envisagés par les moyens de fabrication classiques. Une géométrie basée sur la théorie des développantes sphériques, combinée à une spirale logarithmique a été développée, puis usinée. De plus, des corrections de bombé ou de profil peuvent être définies afin d'éviter les contacts en bords de denture. Ce type de géométrie analytique offre des possibilités plus simples d'optimisation de l'engrènement. L'optimisation des surfaces peut être réalisée à l'aide du modèle d'engrènement quasi-statique sous charge développé dans le cadre de cette thèse. L'environnement de l'engrenage est pris en compte dans la simulation : déformation des arbres, des dentures et de leurs supports (jantes et voiles) ainsi que les déformations locales de contact. La méthode des coefficients d'influence est utilisée pour résoudre le partage des charges entre toutes les dents instantanément en contact. Une méthode originale, utilisant sur un seul calcul élément finis et la définition de bases de fonctions, permet de calculer rapidement les flexions de denture dans leur environnement. Les déformations de contacts sont, quant à elles, obtenues par une méthode analytique, basée sur les théories de Boussinesq. De plus, des défauts d'assemblage peuvent être intégrés entre le pignon et la roue spiro-conique. Afin de valider les modèles numériques développés, un banc d'essai a été mis en place, permettant la mesure de l'erreur de transmission et la visualisation des portées. Le banc d'essai est intégré dans une fraiseuse numérique 3 axes : le pignon est monté dans la broche de la fraiseuse, le reste du banc étant bridé sur son plateau. Ainsi, des défauts de montage peuvent être appliqués facilement et précisément.

Dans un contexte de vulnérabilité militaire des sous-marins, les ingénieurs et chercheurs doivent être capables de prédire le comportement des structures fissurées. Ainsi, la modélisation de la transition des changements de modes de propagation de fissure (cisaillement-traction et inversement) des métaux sous sollicitations extrêmes devient un outil incontournable ou essentiel. Des critères tridimensionnels de direction de propagation de fissure développés pour une rupture par cisaillement ou par ouverture sont exposés. Des formules de direction de propagation semi-analytiques et analytiques, fonctions des facteurs d'intensité des contraintes et du coefficient de Poisson, sont ainsi proposées. L'interprétation de ces formules laisse envisager la prise en compte des effets tridimensionnels dans de futures simulations 3D de propagation de fissure. Une étude du problème en deux dimensions est également développée, proposant une formule analytique du critère en cisaillement. De plus un algorithme automatique de transition cisaillement-traction a été implémenté dans le code de calcul de dynamique explicite Europlexus, développé par le CEA. Une méthodologie d'identification des paramètres du modèle pour un matériau donné et pour un cas quasi-statique a été proposée. Confronté à l'interprétation de deux expériences connues de propagation dynamique (expériences de Kalthoff et de Ravichandran), le modèle proposé a montré sa pertinence. De plus, afin de mieux connaître le comportement à rupture de l'acier à Haute Limite Élastique Soudable, deux études expérimentales dédiées au suivi de la propagation dynamique d'un front de fissure ont été développées et validées sur des essais de rupture sous chargement quasi-statique et dynamique de type choc. Cette étude expérimentale a permis d'observer que les branchements de fissures, relevés sur les essais sous chargement quasi-statique, n'apparaissent plus sous chargement dynamique et pour des sollicitations en mode I pur. Les méthodes théoriques et numériques développées dans ces travaux de thèse permettent donc de simuler, automatiquement et avec un unique modèle, les changements de modes de rupture au cours d'une propagation dynamique de fissure. De plus, les protocoles expérimentaux exposés dans ce manuscrit permettent d'appréhender les phénomènes de transition cisaillement-traction en soulevant l'importance de la vitesse de sollicitation et du mode de sollicitation de l'essai.

Ce travail concerne l'étude mathématique et l'analyse numérique d'une méthode de résolution du problème de Cauchy basée sur la minimisation d'une fonctionnelle énergétique. Depuis les travaux de J. Hadamard, le problème de Cauchy est connu pour être mal posé et les méthodes de résolution de ce type de problèmes présentent une importante instabilité numérique dans le cas de données bruitées. Dans le premier chapitre, le problème de Cauchy est introduit et des résultats théoriques classiques sont donnés. La méthode énergétique et le problème de minimisation associé sont présentés, la théorie du contrôle optimal est utilisée pour l'étude mathématique de ce problème de minimisation. Le deuxième chapitre est consacré à l'application de la méthode énergétique pour l'équation de la chaleur stationnaire. Une fois le cadre variationnel défini, la discrétisation éléments finis de la méthode et des estimations d'erreur a priori tenant compte des données bruitées sont données. Lorsque les données sont bruitées, l'erreur atteint une valeurs minimale avant d'exploser numériquement tandis que la fonctionnelle atteint assymptotiquement un seuil dépendant du taux de bruit. Une estimation du seuil atteint par la fonctionnelle en fonction du bruit est donnée et aboutit à la proposition d'un critère d'arrêt pour le processus de minimisation permettant de contrôler l'explosion numérique due au bruit. Enfin, les résultats théoriques sont validés numériquement, la robustesse et l'efficacité du critère d'arrêt proposé sont illustrées par différents tests numériques. La méthode énergétique est ensuite appliquée à l'équation de la chaleur en régime transitoire et est analysée en suivant la méthodologie introduite dans le cas stationnaire.