Ce travail présente une approche modulaire bien adaptée à l'analyse de situations complexes d'assemblages de structures élastiques. L'approche CONTRAST (CONtact TRidimensionnel dans les Assemblages en STatique) prend facilement en compte les non-linéarités locales dues au contact ou à des éléments technologiques particuliers (joints élastomères, boulons fortement précontraints,...). Elle est basée sur une décomposition de l'assemblage en sous-structures et interfaces, qui sont des entités mécaniques à part entière avec leurs propres équations et leurs propres inconnues. Les interfaces jouent un rôle majeur en modélisant simplement les non-linéarités locales. Le schéma itératif de résolution est issu de la méthode à grand incrément de temps (LATIN method). L'approche CONTRAST est ainsi basée sur une formulation et une stratégie parallèles. Ce parallélisme est avant tout utilisé pour introduire flexibilité et modularité dans la définition du problème et autorise l'utilisation de maillages incompatibles sur les interfaces. Il permet aussi une réduction des coûts de calcul par rapport aux techniques éléments finis classiques. Des comparaisons avec un code de calcul industriel valident le traitement des problèmes de contact et de frottement. Des exemples d'assemblages industriels présentent les possibilités de l'approche ainsi que son comportement sur calculateur parallèle

Caractérisation générale de la superplasticité et du procédé de mise en forme des tôles par pression gazeuse. Présentation d'une modélisation numérique permettant une meilleure maîtrise de ce procédé de fabrication notamment utilisé dans l'industrie aéronautique pour le formage de pièces de structures en tôles d'alliages de titane ou d'aluminium.

Le procédé d'hydroformage repose sur la déformation plastique d'une forme première (tôle, tube) sous l'action combinée d'un fluide sous pression et d'un outil rigide. L'objet de ces travaux est de contribuer à la modélisation de ce procédé. Cette contribution s'articule autour de trois problématiques essentielles: la caractérisation rhéologique, la prédiction des défauts de striction et le développement d'algorithmes adaptés à la simulation de l'hydroformage. - L'identification expérimentale des courbes d'écrouissage en grandes déformations est nécessaire afin d'éviter toute extrapolation intempestive au cours d'une simulation numérique. A cet effet, un essai de gonflage circulaire, analysé selon un modèle homogène, a été mis en place. Des résultats très satisfaisants ont pu ainsi être obtenus grâce notamment à la mesure locale du rayon de courbure de la tôle. - La prédiction de la striction au cours de l'opération d'hydroformage est essentielle pour la maîtrise de ce procédé. Dans cette optique, le Critère de Force Maximum Modifié a fait l'objet d'une étude détaillée. Il a été démontré que la forme locale de la surface de charge constitue le facteur le plus influant sur les prédictions du modèle. Le problème de la dépendance de la CLF au trajet a été également abordé, nous avons pu mettre l'accent sur l'utilité d'une représentation des CLF dans l'espace des déformations équivalentes afin de s'affranchir de cette dépendance. Enfin, nous avons proposé un modèle intégrant l'effet de la sensibilité à la vitesse de déformation. La confrontation avec les résultats expérimentaux de tôles en acier a montré une amélioration des prédictions mais le modèle reste relativement conservatif. - Le troisième volet a concerné la modélisation numérique. D'une part, un schéma d'intégration implicite reposant sur une actualisation itérative du maillage et tenant compte des dérivées du domaine a pu être testé. Les résultats ont montré la faible influence des dérivées du domaine. En revanche, l'actualisation itérative semble améliorer la conservation du volume. D'autre part, il est bien établi que certains matériaux peuvent présenter un potentiel de déformation au-delà du pic de pression. Pour une simulation correcte de ce potentiel de déformation, nous avons développé plusieurs variantes de la méthode de continuation. Reposant sur une technique de prédiction-correction, ces méthodes ont permis de dépasser le pic qui reste un point infranchissables pour la méthode conventionnelle de contrôle de pression.

Travail sur un point important de la construction des ponts : le comportement des bossages d'ancrage de câbles de précontrainte. Le travail théorique prend comme point de départ une recherche très poussée du LCPC du point de vue du calcul élastique, qui comportait une étude expérimentale d'un ancrage par bossage centré dans une plaque appuyée sur la tranche à un about. La thèse se limite à la sollicitation du hourdis par suite de l'application d'une force de précontrainte au moyen d'un bossage (ou de plusieurs). Utilisation de la méthode des coupures modifiée pour une analyse élastique et mise au point d'un programme aux éléments finis non linéaires pour l'analyse non linéaire du comportement à la rupture. Étude expérimentale du groupement des bossages.

A la suite d'une synthèse bibliographique sommaire sur le rôle de la vibration dans le compactage des matériaux granulaires, on examine les bases théoriques relatives à la schématisation du comportement de l'ensemble sol-vibrateur à l'aide des modèles à paramètres concentrés du type masse-ressort-amortisseur ; l'analyse du rôle des caractéristiques mécaniques du matériau dans le comportement vibratoire de l'ensemble sol-rouleau permet d'interpréter plusieurs constatations expérimentales antérieures (essais sur le Vibrex). Cette approche met par ailleurs en évidence l'importance de la variabilité de l'aire de contact sol-rouleau. Les variations de l'aire de contact en fonction des caractéristiques mécaniques du matériau et de la charge statique appliquée fait l'objet d'une étude expérimentale (essais en fosse) et théorique au cours de laquelle la recherche de modèles aux éléments finis élasto-plastiques pouvant expliquer les résultats expérimentaux permet de mieux connaître le rôle des principaux paramètres dont dépend le phénomène réel. La dernière partie est consacrée à l'étude, par la méthode des Eléments Finis, des sollicitations dynamiques induites dans le massif de sol (supposé élastique ou visco-élastique). Afin de simuler les conditions aux limites réelles, on est amené à appliquer la méthode des "bords à amortisseurs visqueux" proposée par Lysmer et Kuhlemeyer.

On développe un modèle viscoélastique non linéaire en grandes déformations, appelé modèle pseudo-linéaire. Ce modèle est basé sur un choix particulier d'une mesure de déformation permettant d'écrire le potentiel d'énergie libre sous forme quadratique et la loi de comportement sous forme d'une convolution simple. Afin de valider ce modèle, on compare le potentiel pseudo-linéaire à un potentiel d'énergie libre écrit sous la forme générale d'un développement de Fréchet d'ordre 4. Les fonctions de relaxation de ce potentiel sont écrites sous forme d'exponentielles décroissantes à un temps caractéristique par décade. Cette comparaison montre une bonne concordance entre les deux descriptions. La formulation proposée permet de calculer simplement la dissipation dans le matériau. Une illustration est donnée dans le cas d'oscillations harmoniques. On développe alors une formulation variationnelle associée au modèle pseudo-linéaire basée sur la minimisation de l'énergie totale. On effectue une approximation linéaire en temps et on propose une procédure numérique de résolution. On construit alors un code de calcul par éléments finis et on le valide dans le cas de comportements simples en comparant les solutions analytiques aux résultats numériques.

On propose une modélisation numérique des plaques multicouches en matériaux composites lorsqu'elles sont soumises à l'impact de charges dynamiques. Cette modélisation permet de déterminer les efforts engendrés par l'impact, le but étant de pouvoir étudier certains endommagements comme le délaminage. La nouvelle modélisation utilisée appelée modélisation multiparticulaire des matériaux multicouches (M4), permet de modéliser l'objet tridimensionnel par une géométrie bidimensionnelle tout en augmentant le nombre de champs cinématiques. Elle fait apparaître de manière naturelle des efforts au niveau des interfaces entre les couches. Ces efforts peuvent être reliés au phénomène de délaminage. L'analyse des résultats de la littérature permet de montrer que cette modélisation est pertinente dans le cas des chocs relativement rapides. Un outil numérique a été réalisé sur la base de la méthode des éléments finis en dynamique. Le code éléments finis obtenu est ensuite validé, puis appliqué à l'analyse des efforts dans une plaque multicouche impactée. Nous montrons sur un exemple, la pertinence de cette modélisation et sa capacité à fournir de manière simple et rapide, des informations utiles pour l'étude des différents aspects du délaminage (initiation, propagation, étendue).

La thèse a pour thème la modélisation numérique et la prédiction de l'endommagement (l'usure ou la fissuration) dû au contact avec frottement avec petits débattements (fretting). Elle comporte deux parties. Dans la première partie, une nouvelle méthode de calcul éléments finis est développe pour les chargements cycliques. Contrairement aux algorithmes classiques incrémentaux, cette méthode est basée sur un algorithme original qui permet d'obtenir directement l'état stabilisé mécanique éventuel. La méthode est itérative et consiste à chercher sur tout le cycle de chargement, en un premier temps des solutions statiquement et cinématiquement admissibles, et puis après des solutions plastiquement admissibles. A la fin de chaque itération, la condition de périodicité est imposée. Et ainsi, lorsque l'algorithme converge, on obtient des champs de contraintes et de déformations périodiques et vérifiant l'équilibre et la loi de comportement sur tout le cycle. La comparaison des temps de calculs par rapport à la méthode incrémentale a démontré la performance de cette méthode. Dans la deuxième partie, cette méthode est mise en oeuvre dans le cas du chargement cyclique mobile de fretting. Ainsi des cartes d'adaptation et d'accommodation sont établies en fonction de divers paramètres. Les résultats de prédiction de l'initiation des fissures de fatigue faite par l'application du critère multiaxial de DANG VAN sont en très bon accord avec les résultats expérimentaux. L'usure est obtenue pour une réponse mécanique de type accommodation et elle est donc le résultat d'un mécanisme de fatigue oligocyclique.

En mécanique des matériaux, la volonté actuelle est de chercher à mieux comprendre certains phénomènes macroscopiques en étudiant la microstructure et les phénomènes physiques à l'échelle microscopique. Cette approche est rendue possible par les nombreux développements dans les techniques d'homogénéisation en mécanique multi-échelles. Dans la présente thèse, on travaille sur les élastomères chargés pour lesquels de nombreuses propriétés mécaniques sont étroitement liées à l'arrangement des particules et agrégats de noir de carbone dans la matrice élastomère, à l'échelle microscopique. La démarche adoptée s'articule autour de deux objectifs principaux. Le premier consiste à modéliser la morphologie de la microstructure du matériau. Pour cela, on met en place un modèle aléatoire de morphologie mathématique décrivant chaque échelle de la microstructure. On propose une méthode d'identification de ce modèle à partir de l'exploitation d'images de microscopie à transmission. Cette méthode trouve son originalité dans le fait d'optimiser le modèle en simulant des images de microscopie à transmission, sur lesquelles on mesure les moments statistiques d'ordre deux et trois que l'on compare aux moments expérimentaux. Cette méthode permet finalement de simuler des microstructures dont la morphologie est très proche de celle du matériau réel. Le second objectif consiste, à partir des microstructures ainsi simulées, à déterminer les propriétés effectives du matériau par le calcul par éléments finis à travers la notion de Volume Elémentaire Représentatif (VER). L'idée est de déterminer la taille du VER par une méthode numérique et statistique en cherchant à estimer les propriétés effectives par une approche de type Monte-Carlo, pour des simulations de microstructures de tailles croissantes. La détermination du VER porte sur les modules élastiques et la conductitivité électrique. De nombreux outils tels que le maillage par éléments finis ou encore le calcul parallèle appliqué aux matériaux présentant un fort contraste sur les propriétés entre les phases, ont été explorés afin de mener à bien ce dernier objectif.

Lors de la fabrication d'une bouteille par injection-soufflage, le conditionnement thermique de la préforme joue un rôle essentiel. Nous proposons une procédure d'optimisation numérique permettant de calculer automatiquement les paramètres de réglage du four infrarouge. Dans un premier temps, l'algorithme d'optimisation de Nelder-Mead est couplé avec des simulations éléments finis de l'étape de soufflage, réalisées avec ABAQUS®. L'objectif est de calculer la distribution de température optimale de la préforme, permettant d'uniformiser l'épaisseur de la bouteille. Dans un second temps, un algorithme de programmation quadratique séquentielle est couplé avec un modèle numérique de chauffage infrarouge 3D développé au laboratoire. Cette méthode permet de calculer les paramètres optimaux pour le réglage du four. Des mesures expérimentales réalisées sur une machine de soufflage semi-industrielle ont permis de valider qualitativement notre approche pour une bouteille de géométrie simple. Les propriétés radiatives du PET sont mesurées par spectrométrie infrarouge. Ces mesures sont exploitées pour calculer l'absorption spectrale du rayonnement. Le modèle de chauffage est validé à l'aide de mesures thermographiques. Un capteur a été développé pour mesurer la résistance thermique de contact entre le polymère et le moule. Le débit d'air injecté dans la préforme a été mesuré, puis appliqué en tant de donnée d'entrée. La pression de soufflage est alors automatiquement calculée à chaque itération par un modèle thermodynamique. Les cinématiques de mise en forme, ainsi que les distributions d'épaisseurs calculées par le modèle sont conformes à celles mesurées sur un pilote de laboratoire.