Le procédé " Resin Transfer Molding " consiste à injecter un polymère thermodurcissable à travers des fibres de renfort, qui polymérise ensuite pour former une pièce composite. L'objectif est de modéliser numériquement l'écoulement à différentes échelles, celle des mèches composants le tissu (mésoscopique) et celle des fibres composants les mèches (microscopique). Cette thèse se décompose en deux parties. La première concerne le calcul de perméabilité d'un volume élémentaire représentatif par prise de moyenne des champs de vitesse et de pression, dans le cadre de la méthode d'immersion de domaines. A l'échelle microscopique, le calcul de perméabilité a été validé en utilisant des lois analytiques. A l'échelle mésoscopique, un couplage entre les équations de Stokes et de Darcy, pour les écoulements entre les mèches et à l'intérieur des mèches a été effectué. La seconde partie aborde l'étape d'imprégnation du renfort à l'échelle microscopique, ce qui inclut la modélisation de l'avancée du front de matière entre les fibres avec prise en compte de la tension de surface. Nous avons implémenté des méthodes pour prendre en compte les phénomènes capillaires. Une méthode d'imposition de l'angle de contact statique de l'interface avec une surface solide, et la condition aux limites de Navier permettant d'imposer un glissement à la ligne triple, ont été implémentées et validées. Tous les développements ont été effectués dans une méthode éléments finis mixtes linéaires en vitesse pression, stabilisée par une fonction bulle (MINI-élément), et en utilisant la méthode d'immersion de domaines.

Les Surfaces à Haute Impédance (SHI) ont été largement étudiées pour améliorer toutes sortes de performances des antennes, comme le gain, le facteur de qualité, les formes et dimensions. L'objectif de cette thèse est de modéliser les structures de SHI et de caractériser leurs performances en vue de futures applications aux antennes.Après une brève introduction aux structures SHI et une étude de quelques modèles analytiques fréquemment traités dans la littérature, deux nouvelles méthodes numériques sont proposées pour calculer l'impédance de surface de structures SHI. Ces deux méthodes (dites " méthode du flux de Poynting " et " méthode / ") sont validées sur des structures symétriques, puis mises en service sur des structures de SHI asymétriques. Elles sont également validées par comparaison de résultats analytiques, numériques et expérimentaux.Nous présentons ensuite un modèle équivalent basé sur l'idée de remplacer les structures hétérogènes de SHI par une surface homogène, caractérisée par son impédance surfacique. Ce modèle nous permet d'avoir une prédiction avec un temps de calcul et une occupation de mémoire PC largement réduits.

La plupart des méthodes de calcul utilisées en mécanique des solides, pour résoudre des problèmes non stationnaires sont basées sur un découpage entre le temps et l'espace. Dans la majorité des logiciels existants, les calculs s'effectuent par une discrétisation de l'espace en utilisant des méthodes, comme, par exemple, la Méthode des Éléments Finis et ensuite, une discrétisation en temps par différences finies. On pourrait se demander : que se passerait-il si la discrétisation en espace et en temps s'effectuait sans aucune séparation ? De là, beaucoup de questions se posent, entre autres : serait-il possible de gagner en temps de calcul et en précision ? En résumé, quel est l'intérêt de cette méthode ? La Méthode des Éléments Finis Espace-Temps (STFEM) que nous présenterons dans ce mémoire se distingue des autres méthodes du fait de la possibilité de non séparation des variables espace et temps.

Le comportement mécanique des matériaux granulaires, et en particulier des poudres métalliques ductiles, peut être appréhendé par l'intermédiaire de simulations discrètes dans lesquelles les grains sont modélisés. Cette méthode a l'avantage de permettre l'exploration de chemins de chargement inaccessibles à une approche expérimentale. La méthode développée dans cette thèse consiste à modéliser des assemblages de sphères maillées dans le logiciel d'éléments finis ABAQUS pour effectuer des simulations en appliquant divers chemins de chargement. La première partie de ce travail consiste à utiliser cette méthode pour enrichir la méthode des éléments discrets en l'étendant au domaine de la compression à haute densité relative. La méthode de simulations discrètes par éléments finis donne des informations précieuses sur les déformations non-linéaires des grains aux contacts. Ces informations permettent de proposer un modèle de contact qui gère les interactions complexes entre les diérents contacts à haute densité relative. Par la suite, ce modèle est introduit dans le code d'éléments discrets YADE, et les résultats des deux méthodes peuvent être comparés. Dans la seconde partie de cette thèse, un échantillon numérique est soumis à des sollicitations dans toutes les directions de l'espace des contraintes, pour obtenir les surfaces de charge macroscopiques. L'accent est mis sur l'évolution des surfaces de charge lors de chargements complexes, avec des changements de direction, afin d'obtenir des indications sur la mémoire du matériau. Enfin, l'observation de la microstructure de l'échantillon permet de comprendre les phénomènes micromécaniques qui sont à l'origine de l'évolution macroscopique des surfaces de charge.

Le comportement nonlinéaire des sols observé lors des mouvements sismiques forts est maintenant bien admis et le déploiement des puits accélérométriques a permis des analyses détaillées de la propagation des ondes ainsi qu'une évaluation quantitative des paramètres physiques tels que la vitesse de cisaillement et de compression des ondes et les facteurs d'amortissements en fonction de la déformation. En dépit du nombre grandissant d'études sur ce phénomène, sa connaissance est encore récente et les recherches sur les données de puits accélérométriques restent une étape importante vers la compréhension du comportement complexe in-situ des sédiments soumis à des mouvements sismiques forts.L'objectif de ces travaux est triple. Premièrement, un code d'inversion par algorithme génétique est développé afin d'inverser des données de puits accélérométriques via la théorie des matrices de propagation de Thomson-Haskell. Cette technique nous permet dans un premier temps de valider la structure en une dimension (1D) (e.g., vitesse des ondes de cisaillement, facteurs d' amortissements) d'un puits accélérométrique dans le domaine linéaire et dans un second temps de mettre en évidence de manière quantitative le comportement nonlinéaire des sédiments lors du séisme de Fukuoka, 2005, Japon. Deuxièmement, les résultats de l'inversion sont utilisés pour tester des lois de comportement simples et avancées en utilisant la Méthode des éléments Finis. Les résultats montrent clairement que l'hypothèse bi-linéaire de la loi de comportement simple produit des séries temporelles non réalistes en vitesse et en accélération. L'utilisation d'une loi de comportement avancée mène à de meilleurs résultats, cependant, le nombre de paramètres ajustables pour obtenir des résultats consistants avec l'observation est un obstable inévitable. Troisièmement, afin d'étendre l'étude des effets de site à des dimensions supérieures, des codes 2D et 3D de la Méthode en éléments Spectraux sont développés et validés en comparant leurs résultats dans le domaine linéaire avec ceux obtenus théoriquement ou via d'autres méthodes numériques.

Les Surfaces à Haute Impédance (SHI) ont été largement étudiées pour améliorer toutes sortes de performances des antennes, comme le gain, le facteur de qualité, les formes et dimensions. L'objectif de cette thèse est de modéliser les structures de SHI et de caractériser leurs performances en vue de futures applications aux antennes.Après une brève introduction aux structures SHI et une étude de quelques modèles analytiques fréquemment traités dans la littérature, deux nouvelles méthodes numériques sont proposées pour calculer l'impédance de surface de structures SHI. Ces deux méthodes (dites " méthode du flux de Poynting " et " méthode / ") sont validées sur des structures symétriques, puis mises en service sur des structures de SHI asymétriques. Elles sont également validées par comparaison de résultats analytiques, numériques et expérimentaux.Nous présentons ensuite un modèle équivalent basé sur l'idée de remplacer les structures hétérogènes de SHI par une surface homogène, caractérisée par son impédance surfacique. Ce modèle nous permet d'avoir une prédiction avec un temps de calcul et une occupation de mémoire PC largement réduits.

Dans le cadre de cette thèse, une modélisation par éléments finis continue renforcée discrète est développée. Le modèle est basé sur le concept de la cellule élémentaire, avec utilisation de connecteurs pour prendre en compte la rigidité en tension de la structure, et modéliser l'anisotropie du tissu et un élément coque pour décrire la rigidité en cisaillement et gérer le frottement et le contact avec les outils lors de la simulation du procédé de préformage. Les différents éléments de ce maillage spécifique sont générés automatiquement à l'aide des scripts pythons. Associés à ce choix de modélisation, les caractéristiques du comportement sont déterminées à l'aide d'essais expérimentaux conduits sur le renfort et d'une démarche de méthode inverse. Le comportement de ces différents éléments est choisi d'être linéaire élastique en premier lieu, et par la suite à évoluer pour devenir non linéaire. Ces développements numériques sont réalisés dans le logiciel Abaqus/ Explicit. Dans un second temps le modèle développé est utilisé dans le cadre de simulation de l'étape de préformage des renforts secs. A ce titre un démonstrateur expérimental développé dans le cadre de la thèse sert à obtenir des résultats expérimentaux pour corréler ces travaux de simulations. Des comparaisons avec des simulations mono-couche, hémisphérique, permettre de mettre en évidence l'influence des paramètres de simulations mises en jeu. Dans le cadre de préformages de plusieurs plis, l'influence des caractéristiques du contact / frottement est mise en évidence. Finalement une procédure d'automatisation afin de rendre la simulation plus ergonomique dans un contexte industriel est présentée. Cette procédure consiste à générer le maillage spécifique et mettre la simulation en donné dans Catia. Les simulations peuvent ainsi se réaliser sans le passage par Abaqus/CAE.

Les travaux présentés dans ce mémoire s'inscrivent dans l'étude et l'amélioration des modèles d'endommagement continus régularisés (non locaux), l'objectif étant d'étudier la transition entre un champ d'endommagement continu défini sur l'ensemble d'une structure et un modèle discontinu de fissuration macroscopique.La première étape consiste en l'étude semi-analytique d'un problème unidimensionnel (barre en traction) visant à identifier une famille de lois d'interface permettant de basculer d'une solution non homogène obtenue avec un modèle continu à gradient d'endommagement vers un modèle discontinu de fissuration cohésive. Ce passage continu / discontinu est construit de telle sorte que l'équivalence énergétique entre les deux modèles soit assurée, et reste exacte quelque soit le niveau de dégradation atteint par le matériau au moment où cette transition est déclenchée.Cette stratégie est ensuite étendue au cadre 2D (et 3D) éléments finis dans le cas de la propagation de fissures rectilignes (et planes) en mode I. Une approche explicite basée sur un critère de dépassement d'une valeur " critique " de l'endommagement est proposée afin de coupler les modèles continus et discontinus au sein d'un même calcul quasi-statique par éléments finis. Enfin, plusieurs résultats de simulations menées avec cette approche couplée sont présentés.

Le découpage de tôles minces en alliages de cuivre est un procédé de mise en forme très utilisé pour la fabrication de composants électroniques. Deux aspects permettent la caractérisation du procédé : d'une part, le profil de découpe représente la qualité du produit et d'autre part, l'effort de découpage conduit à la définition de l'outillage (matériaux du poinçon et de la matrice, la vitesse et la puissance des presses, etc.). La modélisation numérique permet d'obtenir des réponses rapides sur les questions concernant ces deux aspects. Au sein de l'équipe modélisation et mise en forme des matériaux du Laboratoire de Mécanique Appliquée R. Chaléat, un code de calcul par éléments finis dédié à la simulation bidimensionnelle du découpage a été développé depuis quelques années (Blankform©).<br />Afin d'obtenir une meilleure prédiction de l'effort de découpage par simulation numérique, notre démarche commence par la caractérisation expérimentale du comportement thermo-élasto-viscoplastique de la nuance de cuivre CuNiP par des essais de traction à températures et à vitesses de déformation imposées. Suite à ces essais, une loi de comportement particulière permettant la prise en compte de la sensibilité du matériau à la vitesse et à la température a été proposée. Après l'implémentation de cette loi dans le code de calcul, l'application du modèle au découpage d'une tôle mince en CuNiP, dans un cadre adiabatique, montre que l'effort maximal de découpage est correctement prédit.<br />Lors d'une opération de découpage, certaines nuances de cuivre, comme le Cua1, présentent un phénomène de recristallisation dynamique. Les essais de traction sur cette nuance nous ont permis de caractériser le comportement mécanique du matériau. Sur cette base spécifique, un modèle de comportement thermo-élasto-viscoplastique avec prise en compte de la recristallisation dynamique a été proposé dans le cadre de la thermodynamique des processus irréversibles. L'application de ce modèle au découpage montre qu'il permet de prédire qualitativement les phénomènes observés expérimentalement pour cette nuance.<br />La modélisation bidimensionnelle du découpage est souvent utilisée pour des empreintes de formes relativement simples. Cependant, pour des formes complexes, il est nécessaire de résoudre le problème par une analyse tridimensionnelle. A partir d'un cas d'étude industriel, une modélisation tridimensionnelle du découpage a été réalisée afin d'étudier l'influence des paramètres technologiques sur la distribution des contraintes et des déformations. Les résultats numériques obtenus corroborent bien les constatations expérimentales en termes de profil de découpe.

Une méthode est proposée pour la modélisation d’unassemblage multimatériaux innovant visant des applicationsdans l’allègement structurel des véhicules. Dans cetassemblage une partie composite thermoplastique, est fixée àune plaque acier texturée par la technologie CMT pins.L’interface est particulièrement complexe et non linéaire : unereprésentation fine du comportement local serait extrêmementpénalisante en temps de calcul. Dans cette optique il a étéchoisi d’orienter la méthodologie vers une modélisation la plussimple possible tout en conservant de bons résultats globaux.Pour ce faire, en s’inspirant de l’état de l’art existant sur lesmultimatériaux, une campagne expérimentale a été menée surdes éprouvettes longitudinales à double recouvrement afin decaractériser cette interface. Une fois les mécanismes locauxinvestigués, deux modèles sont proposés. Le premiernumérique basé sur la méthode des éléments finis etl’introduction d’éléments cohésifs nous permet de valider leshypothèses de modélisation tout en quantifiant la répartitiondes efforts entre les picots. Le deuxième se base sur le calcullocal d’un Volume Élémentaire Représentatif pour établiranalytiquement la loi de comportement de l’interface. Cette loiest ensuite introduite sous la forme d’un ressort non linéaire ausein d’un modèle numérique simplifié de l’éprouvette. Pourfinir ces approches sont appliquées au cas d’étude industriel etles résultats sont validés par une deuxième campagneexpérimentale. / A method is proposed to simulate an innovative multimaterialassembly which has applications in structural lightweight forvehicles. In this assembly, a thermoplastic composite part isfixed on a steel plate, textured by the CMT pins technology.This is an especially complex and nonlinear interface: a finerepresentation of local behaviour would be extremely costlyfor calculation. So, it has been chosen to investigate a model assimple as possible which still demonstrates accurate globalresults.An experimental campaign on double lap shear specimen,inspired by existing state of the art on multimaterial has beenset up to characterize this interface. Once local mechanismshave been understood, two models are proposed and compared.The first is numerical and based on finite elements method andcohesive elements. It allows us to validate the modelhypotheses while describing the effort repartition between thepins. The second one is based on a Representative VolumeElement. It establishes analytically the behaviour law of theinterface. This law is then inserted inside a simplifiednumerical model of the specimen by means of a nonlinearspring. To conclude, these approaches are applied to theindustrial case of study and the result have been validated by asecond experimental campaign.