Cette thèse étudie la différenciation des lithologies et fluides à partir de données sismique et de puits dans le réservoir gréseux peu consolidé à huile lourde Morichal (ceinture de l'Orénoque, Venezuela). Dans ces réservoirs, cette discrimination est complexe parce que les lithologies ont des propriétés élastiques similaires et qu'il en est de même pour les sables saturés en eau et ceux saturés en huile lourde (mêmes densités). Dans une première partie, cette problématique de discrimination est analysée au niveau des puits par la détermination d'électrofaciès puis en étudiant la relation entre propriétés pétrophysiques et élastiques. De petites variations de la taille des grains de la roche ou du pourcentage de ciment génèrent des changements de vitesses qui pourraient être utilisés pour discriminer lithologies et fluides. Afin de mieux comprendre pourquoi un premier calage puits-sismique a échoué, une deuxième partie de la thèse a été dédiée à la construction d'un modèle 3D du réservoir, habillé en faciès, propriétés pétrophysiques et élastiques. Ce modèle 3D est ensuite utilisé pour réaliser des tests de modélisations sismiques 1D et 2D, afin d'étudier les effets de l'anisotropie et de l'atténuation sur le signal sismique. En effet, l'anisotropie peut entraîner des changements de phase pour les grands angles d'incidence, et l'atténuation due à la présence d'huile lourde peut causer un retard de phase et un affaiblissement de l'amplitude. Ces deux phénomènes peuvent donc perturber le calage puits-sismique. Après analyse de ces tests, des pistes d'amélioration du calage puits-sismique dans ce contexte géologique sont proposées pour le futur.

L'instabilité des cavités souterraines (mines, carrières, tunnels,...) peut induire les mouvements de terrains d'amplitude suffisante pour endommager les bâtiments et les infrastructures en surface. Les méthodes traditionnelles, utilisées dans les pratiques d'ingénieur pour prévoir les déformations dans les structures, sont basées sur les caractéristiques des mouvements de terrain en condition de terrain vierge sans prendre en compte l'effet de la présence des structures en surface. L'objectif de cette thèse est de prédire les déformations des ouvrages en tenant compte de l'influence de l'interaction sol-structure, d'une part ; et d'évaluer la performance d'une solution de protection (tranchée périphérique), d'autre part. Cela a été achevé par la réalisation d'études paramétriques utilisant deux approches complémentaires : une approche expérimentale à l'aide d'un modèle réduit physique 3D sous gravité normale et une modélisation numérique 3D par la méthode des éléments finis. En particulier l'effet d'un certain nombre de paramètres géométriques et mécaniques a pu être investigué dans l'étude de l'interaction sol-structure : la position de la structure par rapport à la cuvette d'affaissement, le poids de la structure et la raideur relative entre le sol et la structure. Concernant l'étude de l'efficacité de tranchées périphériques, l'effet de la position de la structure, de la position de la tranchée vis-à-vis de la structure et de la rigidité de la tranchée a été analysé. Les résultats obtenus ont abouti à une meilleure compréhension du problème d'interaction sol-structure et ont montré l'importance de cet effet qui doit être pris en compte dans l'évaluation de la vulnérabilité du bâti. Le transfert des mouvements du sol à la structure est faible (moins de 2,5%), dans le cas modélisé : structure rigide et interface glissante. Les différents résultats ont permis par ailleurs de mettre en évidence l'efficacité de la tranchée périphérique pour réduire les sollicitations affectant les structures. La tranchée doit être remplie avec un matériau très déformable et surtout placée à une distance de l'ordre d'un mètre de la structure.

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.

La thèse présente une méthodologie pour construire des éléments finis de type " solide-coque " avec intégration réduite en vue des applications à la simulation de la mise en forme tel que l'emboutissage des tôles où ces éléments finis doivent présenter de bonnes aptitudes à modéliser la flexion mais également les situations de laminage de la tôle. A partir des éléments volumiques à 8 nœuds et 3 degrés de liberté par nœud (les 3 composantes du déplacement), un neuvième nœud est rajouté au centre de l'élément. Ce neuvième nœud n'est pourvu que d'un seul degré de liberté, le déplacement le long de la direction de l'épaisseur. Cette direction privilégiée a un nombre de points d'intégration supérieur ou égal à 3 mais l'intégration est réduite au centre de l'élément diminuant très sensiblement les temps CPU par rapport à une intégration complète. Un soin particulier a été pris pour contrôler tous les modes à énergie nulle dus à l'intégration réduite. Ce nœud supplémentaire permet une distribution linéaire de la déformation normale. Avec les lois de comportement complètement 3D ces nouveaux éléments solide-coque donnent des résultats similaires en flexion à ceux obtenus avec des éléments coques et état plan de contrainte. Le neuvième nœud joue le rôle d'un paramètre supplémentaire pour l'interpolation quadratique du déplacement dans la direction de l'épaisseur. Ce degré de liberté a une signification physique et un effort équivalent à une pression normale peut être prescrit. Dans les situations de pression normale et dans le cas du contact, la contrainte normale obtenue est physique ce qui n'est pas le cas de nombreux éléments solide-coque de la littérature. Le pincement ou le laminage des tôles est correctement modélisé. Pour valider ces éléments, un module d'emboutissage en U avec passage et laminage de la bande de tôle sur des rouleaux a été construit au laboratoire. La comparaison entre les efforts d'emboutissage calculés et mesurés est très bonne ainsi que la géométrie des bandes de tôle obtenue après retour élastique.

Dans le domaine de la mise en forme industrielle, des outils de simulation comme le logiciel Pam-stamp 2G permettent entre autres le prototypage et l’optimisation numérique des produits, réduisant ainsi les coûts expérimentaux de mise au point. Les éléments finis de type coques en hypothèse d’état plan de contrainte demeurent les plus utilisés car ils permettent une prise en compte réaliste des déformations majeures de membrane et de flexion. Cependant, de par leur définition, la contrainte normale pouvant apparaitre en cas de compression du flan dans la direction de l’épaisseur ou encore de flexion extrême sur petit rayon est systématiquement omise. De plus, Il existe de nouveaux procédés de mise en forme de tôle, comme le pliage/sertissage/emboutissage avec laminage et/ou écrasement ainsi que l’hydroformage, qui ne peuvent pas être traités avec ces formulations d’éléments coques classiques. L’utilisation de couches d’éléments volumiques est souvent considérée comme une alternative non convenable aux simulations de ces procédés. Outre le très haut coût CPU, s’ajoutent le rendu parfois non réaliste ainsi que la complexité liée à la découpe du maillage du flan. Récemment, des éléments de type “solid-shell” ont été mis à contribution mais requièrent des améliorations quant à leurs lois de comportements. L’objectif d’ESI group consiste en l’élaboration, l’évaluation, l’implémentation et la validation industrielle d’une nouvelle formulation d’élément fini. Cet élément devra permettre la gestion d’une éventuelle variation d’épaisseur avec prise en compte réaliste de la contrainte normale, tout en assurant des résultats dignes de ceux d’une coque conventionnelle en flexion. Une nouvelle formulation de type coque 3D est ainsi d’abord proposée. Des éléments coques triangulaires et quadrangulaire en théorie de Mindlin et de Kirchhoff sont utilisés. Cette approche est d’abord développée dans un solveur quasi-statique implicite de l’INSA de Lyon pour validation numérique sur cas académiques linéaires et non linéaires de référence. Une validation expérimentale sur opération de pliage dépliage est également réalisée. Aux vues des non linéarités dues au contact avec frottement, grandes déformations et grands déplacements posant des soucis de convergence en implicite au cours de la simulation d’opérations de mise en forme, la méthode est ensuite développée dans le solveur explicite de l’INSA de Lyon. Seuls les éléments en théorie de Mindlin sont considérés. Les particularités liées à cette méthode de résolution dynamique comme la matrice de masse, le pas de temps critique et l’optimisation du CPU sont traitées. Une nouvelle méthode de contact dédiée aux opérations de mise en forme impliquant du laminage et/ou de l’amincissement est également proposée. Elle permet une transition automatique d’éléments standards vers des éléments coque "3D", palliant ainsi le remaillage. Des essais de mise en forme en U avec ou sans laminage apportent une validation expérimentale concernant le retour élastique. Une fois validées, les techniques et formulations les plus abouties sont implémentées dans le code industriel dédié à la mise en forme Pam-stamp 2G v2011. Après une vérification sur tests de référence, des applications sont enfin menées sur cas critiques inspirés de procédés industriels complexes et nécessitant essentiellement une loi de comportement 3D. / In the sheet metal forming industry, shell elements in plane-stress assumption are employed, as they perform quite well in simulating the major membrane and flexural large deformations involved. However, the normal stress, caused by compression along thickness direction of the blank or local high bending over very small radii, is hence systematically omitted. Besides, when it comes to unusual and challenging processes such as hydro-forming, thinning/thickening, forming with ironing, bottoming and so on, makeshift solutions such as layers of 3D solid hexahedrons or even recent “solid-shell” elements are no longer appropriate. An innovative 3D finite element formulation methodology overcoming the overcoming the plane-stress definition of classification shell elements, while keeping their very good bending assets is first proposed in this work. The method basically consists in adding a central node endowed with two degrees of freedom at the element center. These two extra translations normal to the element mid-plane give a new quadratic displacement field along the shell normal direction. A derivative normal strain can hence be expressed and a linear normal stress comes via a full 3D constructive law. A very pioneering contact technique, dedicated to forming processes with ironing, thinning/bottoming operations and allowing a usual-to-enhanced automatic element switch is developed as well. Once widely assessed, most interesting achievements are implemented in the dynamic explicit industrial code Pam-stamp 2 G v2011 and evaluated over critical industrial forming processes that require essentially a full 3D strain-stress behavior.

Le sujet concerne l’influence des effets tridimensionnelles sur les champs de déformation et de contrainte au voisinage d’une éprouvette chargée en mode I. La loi de comportement caractéristique des polymères amorphes avec un seuil d’écoulement viscoplastique suivi d’un adoucissement et d’un durcissement à mesure que la déformation augmente est prise en compte. La loi est implantée dans une UMAT abaqus. Les champs au voisinage de l’entaille sont analysés et les résultats 3D comparés à ceux correspondant à un calcul 2D sous l’hypothèse de déformation plane. L’influence de l’épaisseur de l’échantillon est étudiée et nous montrons qu’au-delà d’un rapport épaisseur t sur rayon d’entaille rt, t/rt>20, les champs de déformation plastique sont qualitativement similaires entre les calculs 3D et 2D. En revanche, nous montrons que la répartition des contraintes et notamment celle de la contrainte moyenne est sur-estimée avec un calcul 2D en comparaison à une simulation 3D. Nous prenons en compte par la suite la rupture par craquelage, modélisée avec un modèle cohésif. Une étude paramétrique est menée afin de définir une procédure d’identification des paramètres caractéristiques du modèle cohésif. Par ailleurs les simulations montrent qu’au-delà d’un rapport t/rt supérieur à 20, une ténacité minimum peut être estimée : ceci constitue un résultat important pour la détermination expérimentale de la ténacité des polymères ductiles. / We investigate 3D effect of crack tip palsticity and the influence of the thickness on 3D glassy polymer fracture. The characteristuc constitutive law with a viscoplastic yield stress followed by softening and progressive hardening is accoutne for and implemented in a UMAt routine, in abaqus. The crack tip fields are investigated and 2D plasne strain versus 3D calculations compared. Qualitatively, the palstic distribution are comparableas soon as the ratio thickness over crack tip radius is larger than 20. However, the 2D calculations over estimate the stress distribution compared to the 3D cases. We have accounted for failure by crazing that is described with a cohesive models. A parametric study sheds light on the methodology to use for the calibration of the cohesive parameters. The simulations show that for a ratio thickness over craci tip radius larger than 20, a minimum tuoghness can be observed. This results has implication on the definition of a thickness larger enough experimentally.

De nos jours, l'intégration de pièces composites dans les produits intéresse de plus en plus les industriels, particulièrement dans le domaine des transports. En effet, ces matériaux présentent de nombreux avantages, notamment celui de permettre une diminution de la masse des pièces lorsqu'ils sont correctement exploités. Pour concevoir ces pièces, plusieurs procédés peuvent être utilisés, parmi lesquels le RTM (Resin Transfer Molding) qui consiste en la mise en forme d'un renfort sec (préformage) avant une étape d'injection de résine. Cette étude concerne la première étape du procédé RTM, celle de préformage. L'objectif est de mettre en œuvre une stratégie efficace conduisant à la simulation par éléments finis de la mise en forme des renforts à l'échelle mésoscopique. A cette échelle, le renfort fibreux est modélisé par un enchevêtrement de mèches supposées homogènes. Plusieurs étapes sont alors nécessaires et donc étudiées ici pour atteindre cet objectif. La première consiste à créer un modèle géométrique 3D le plus réaliste possible des cellules élémentaires des renforts considérés. Elle est réalisée grâce à la mise en œuvre d'une stratégie itérative basée sur deux propriétés. D'une part, la cohérence, qui permet d'assurer une bonne description du contact entre les mèches, c'est-à-dire, que le modèle ne contient ni vides ni interpénétrations au niveau de la zone de contact. D'autre part, la variation de la forme des sections de la mèche le long de sa trajectoire qui permet de coller au mieux à la géométrie évolutive des mèches dans le renfort. Grâce à ce modèle et à une définition libre par l'utilisateur de l'architecture tissée, un modèle représentatif de tout type de renfort (2D, interlock) peut être obtenu. La seconde étape consiste à créer un maillage hexaédrique 3D cohérant de ces cellules élémentaires. Basé sur la géométrie obtenue à la première étape. L'outil de maillage créé permet de mailler automatiquement tout type de mèche, quelle que soit sa trajectoire et la forme de ses sections. La troisième étape à franchir consiste, à partir du comportement mécanique du matériau constitutif des fibres et de la structure de la mèche, à mettre en place une loi de comportement du matériau homogène équivalent à un matériau fibreux. Basé sur les récents développements expérimentaux et numériques en matière de loi de comportement de structures fibreuses, un nouveau modèle de comportement est présenté et implémenté. Enfin, une étude des différents paramètres intervenant dans les calculs en dynamique explicite est réalisée. Ces deux derniers points permettent à la fois de faire converger rapidement les calculs et de se rapprocher de la réalité de la déformation des renforts. L'ensemble de la chaîne de modélisation/simulation des renforts fibreux à l'échelle mésoscopique ainsi créée est validée par comparaison d'essais numériques et expérimentaux de renforts sous sollicitations simples.

Les bassins de retenue des eaux pluviales sont utilisés pour préserver la qualité des eaux réceptrices par sédimentation pendant le temps de pluie. Cependant, les efficacités du bassins n'étaient pas satisfaisants en raison de la mal compréhension du processus de sédimentation. Afin de mieux comprendre ces processus dans des ouvrages in situ, cette thèse porte à la fois sur des expérimentations in situ et sur les modélisations de l'hydrodynamique et du transport particulaire dans les bassins de retenue pilotes et in situ. Cette recherche s'est appuyée en grande partie sur le bassin Django Reinhardt (BDR) à Chassieu (volume: 32000 m3, surface: 11000 m2) dans le cadre de l'OTHU et sur les données expérimentales obtenues par Dufresne (2008) et Vosswinkel et al. (2012). Les échantillons de sédiments ont été prélevés et leurs caractéristiques physiques ont été analysées en laboratoire dans le but de cerner leur distribution spatiale. Concernant la modélisation numérique, les simulations de l'hydrodynamique en régime permanent ont été réalisées à l'aide du logiciel CFD Fluent et ont été évaluées à partir de l'analyse de corrélation entre le comportement hydrodynamique du bassin et la distribution spatiale des caractéristiques physiques des sédiments. Les conditions limites sur le fond couramment utilisées et largement décrites dans la littérature ont été testées dans le but de représenter la distribution spatiale des sédiments et l'efficacité de décantation du BDR. Les conditions testées sont : i) contrainte de cisaillement critique ou bed shear stress - BSS et ii) énergie cinétique turbulente critique ou bed turbulent kinetic energy - BTKE. L'approche Euler-Lagrange dite " particle tracking " a été mise en œuvre. En raison de l'échec de prédiction des zones de dépôt à l'aide des conditions limites disponibles (BSS et BTKE), une nouvelle relation a été proposée pour estimer le seuil BTKE. La condition à la limite obtenue en utilisant cette nouvelle relation a été testée sur un bassin pilote (Dufresne, 2008) et sur le BDR en régime permanent. Les résultats obtenus n'étaient pas très satisfaisants concernant la prédiction des zones de dépôt et l'efficacité de décantation dans le bassin BDR, même en considérant une distribution granulométrique non uniforme. Afin de mieux prédire les zones de dépôt dans le BDR, une nouvelle méthode a été proposée en considérant le transport des particules, leur décantation et leur érosion en régime transitoire. Sur la base de la méthode proposée pour le transport des particules, la décantation et l'érosion en régime transitoire, plusieurs modélisations avec différentes conditions limites ont été réalisées dans un bassin de retenue pilote rectangulaire (Vosswinkel et al., 2012). Les prédictions des efficacités et des zones de dépôt en régime transitoire avec la méthode proposée sont satisfaisantes

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. / The design of spiral bevel gears is still very complex nowadays because the tooth geometry, and thus the kinematic performance, come from the manufacturing process of this type of gear. The cutting is related to two major manufacturers: Gleason and Klingelnberg. Many machine settings drive directly the shape of teeth surfaces, their optimization is therefore not intuitive. Due to the progress made during the last decade by the CNC machines and the CAM (Computer Aided Manufacturing) softwares, it becomes possible to manufacture spiral bevel gears of quite correct quality on a 5-axis milling machine. A numerical model was developed in order to generate a simplified type Gleason geometry. This last was then manufactured with a 5-axis milling machine. A metrological study, comparing the teeth obtained with the CAD models, was then carried out to prove that the manufacturing by 5-axis milling machine can be an alternative to conventional cutting methods. New types of geometry can be then proposed, which could not be considered by the conventional methods of manufacturing. Geometry based on the theory of spherical involutes, combined with a logarithmic spiral was developed and then manufatured. In addition, profile and crowning modifications can be defined to avoid the tooth edge contacts. This type of analytical geometry offers simpler possibilities for optimizing the meshing. The surface optimization can be achieved using the quasi-static meshing model under load developed in the context of this thesis. The surroundings of the gear are taken into account in the simulation: deformation of the shafts, of the gears and their supports (rims for example) as well as the local contact deformations. The influence coefficient method is used to solve the load sharing between all the teeth instantaneously in contact. An original method, using only one finite element computation and the definition of a set of functions, can quickly calculate the teeth bending, taking into account their surroundings. The contact deformations are obtained with an analytical method, based on Boussinesq theories. In addition, meshing defects can be integrated between the spiral beval pinion and gear. To validate the numerical model, a test bench was achieved, allowing the measurement of the loaded transmission error and the visualization of the contact patterns. The test bench is integrated inside a numerical 3-axis milling machine: the pinion is mounted in the spindle of the milling machine, when the base of the bench is clamped on its plate. Thus, assembly errors can be imposed easily and accurately.

De nos jours, l’intégration de pièces composites dans les produits intéresse de plus en plus les industriels, particulièrement dans le domaine des transports. En effet, ces matériaux présentent de nombreux avantages, notamment celui de permettre une diminution de la masse des pièces lorsqu’ils sont correctement exploités. Pour concevoir ces pièces, plusieurs procédés peuvent être utilisés, parmi lesquels le RTM (Resin Transfer Molding) qui consiste en la mise en forme d’un renfort sec (préformage) avant une étape d’injection de résine. Cette étude concerne la première étape du procédé RTM, celle de préformage. L’objectif est de mettre en œuvre une stratégie efficace conduisant à la simulation par éléments finis de la mise en forme des renforts à l’échelle mésoscopique. A cette échelle, le renfort fibreux est modélisé par un enchevêtrement de mèches supposées homogènes. Plusieurs étapes sont alors nécessaires et donc étudiées ici pour atteindre cet objectif. La première consiste à créer un modèle géométrique 3D le plus réaliste possible des cellules élémentaires des renforts considérés. Elle est réalisée grâce à la mise en œuvre d’une stratégie itérative basée sur deux propriétés. D’une part, la cohérence, qui permet d’assurer une bonne description du contact entre les mèches, c'est-à-dire, que le modèle ne contient ni vides ni interpénétrations au niveau de la zone de contact. D’autre part, la variation de la forme des sections de la mèche le long de sa trajectoire qui permet de coller au mieux à la géométrie évolutive des mèches dans le renfort. Grâce à ce modèle et à une définition libre par l’utilisateur de l’architecture tissée, un modèle représentatif de tout type de renfort (2D, interlock) peut être obtenu. La seconde étape consiste à créer un maillage hexaédrique 3D cohérant de ces cellules élémentaires. Basé sur la géométrie obtenue à la première étape. L’outil de maillage créé permet de mailler automatiquement tout type de mèche, quelle que soit sa trajectoire et la forme de ses sections. La troisième étape à franchir consiste, à partir du comportement mécanique du matériau constitutif des fibres et de la structure de la mèche, à mettre en place une loi de comportement du matériau homogène équivalent à un matériau fibreux. Basé sur les récents développements expérimentaux et numériques en matière de loi de comportement de structures fibreuses, un nouveau modèle de comportement est présenté et implémenté. Enfin, une étude des différents paramètres intervenant dans les calculs en dynamique explicite est réalisée. Ces deux derniers points permettent à la fois de faire converger rapidement les calculs et de se rapprocher de la réalité de la déformation des renforts. L’ensemble de la chaîne de modélisation/simulation des renforts fibreux à l’échelle mésoscopique ainsi créée est validée par comparaison d’essais numériques et expérimentaux de renforts sous sollicitations simples. / Nowadays, manufacturers, especially in transport, are increasingly interested in integrating composite parts into their products. These materials have, indeed, many benefits, among which allowing parts mass reduction when properly operated. In order to manufacture these parts, several methods can be used, including the RTM (Resin Transfer Molding) process which consists in forming a dry reinforcement (preform) before a resin being injected. This study deals with the first stage of the RTM process, which is the preforming step. It aims to implement an efficient strategy leading to the finite element simulation of fibrous reinforcements at mesoscopic scale. At this scale, the fibrous reinforcement is modeled by an interlacement of yarns assumed to be homogeneous and continuous. Several steps are then necessary and therefore considered here to achieve this goal. The first consists in creating a 3D geometrical model of unit cells as realistic as possible. It is achieved through the implementation of an iterative strategy based on two main properties. On the one hand, consistency, which ensures a good description of the contact between the yarns, that is to say, the model does not contain spurious spaces or interpenetrations at the contact area. On the other hand, the variation of the yarn section shape along its trajectory that enables to stick as much as possible to the evolutive shape of the yarn inside the reinforcement. Using this tool and a woven architecture freely implementable by the user, a model representative of any type of reinforcement (2D, interlock) can be obtained. The second step consists in creating a 3D consistent hexahedral mesh of these unit cells. Based on the geometrical model obtained in the first step, the meshing tool enables to mesh any type of yarn, whatever its trajectory or section shape. The third step consists in establishing a constitutive equation of the homogeneous material equivalent to a fibrous material from the mechanical behavior of the constituent material of fibers and the structure of the yarn. Based on recent experimental and numerical developments in the mechanical behavior of fibrous structures, a new constitutive law is presented and implemented. Finally, a study of the different parameters involved in the dynamic/explicit scheme is performed. These last two points allow both to a quick convergence of the calculations and approach the reality of the deformation of reinforcements. The entire chain modeling/simulation of fibrous reinforcements at mesoscopic scale created is validated by numerical and experimental comparison tests of reinforcements under simple loadings.