Les propriétés statiques, dynamiques, rhéologiques et la cinétique de scissions et recombinaisons de micelles linéaires auto-assemblées sont étudiées à l'équilibre et sous-écoulement par simulations sur ordinateur, en utilisant un modèle mésoscopique nouveau. Nous représentons les micelles comme des séquences linéaires de billes browniennes dont l'évolution spatio-temporelle est gouvernée par la dynamique de Langevin. Un algorithme de Monte-Carlo contrôle l'ouverture des liens ou la fusion de deux chaînes par les bouts. Un paramètre cinétique o, qui modélise l'effet d'une barrière le long d'un chemin de réaction, est introduit dans notre modèle. A l'équilibre, nous nous concentrons sur les mécanismes de scission/recombinaison aux temps long et court. Nos résultats montrent que pour les temps plus grands que le temps de vie d'une chaîne moyenne, la cinétique est en accord avec le modèle champ-moyen de Cates. L'étude de fonctions de relaxation macroscopique confirme que nos constantes cinétiques effectives obtenues aux temps longs sont pertinentes pour ces relaxations. Pour la situation hors équilibre, nous étudions les effets du couplage entre un écoulement de cisaillement et la cinétique de scission et recombinaison sur les propriétés structurales et rhéologiques du système micellaire. Nous nous plaçons dans un régime semi-dilué et dynamiquement 'unentangled'. Le paramètre o est choisi de façon à ce que la durée de vie d'une chaîne moyenne soit plus courte que son temps de relaxation de Rouse le plus long. Nos analyses font apparaître une longueur dynamique A, le fiagrnent de chaîne dont la durée de vie TA est égale à son temps de Rouse. Nous trouvons que les propriétés telles que la rhéo-fluidification, l'orientation des chaînes et l'étirement des liens sont des fonctions du taux de cisaillement réduit PA= YT* , alors que la longueur moyenne des micelles est une fonction décroissante du taux de cisaillement, indépendamment de la barrière du processus scission/recombinaison / Statics, Dynamics, and Rheological properties of Micellar solutions by Computer Simulation Statics, dynamics, rheology and scission-recombination kinetics of self-assembling linear micelles are investigated at equlibrium state and under shear flow by computer simulations using a newly proposed mesoscopic model. We model the micelles as linear sequences of Brownian beads whose space-time evolution is governed by Langevin dynamics. A Monte Carlo algorithm controls the opening of a bond or the chain-end fusion. A kinetic parameter o, modelling the effect of a potential barrier along a kinetic path, is introduced in our model. For equilibrium state we focus on the analysis of short and long time behaviors of the scission and recombination mechanisms. Our results show that at time scales larger than the life time of the average chain length, the kinetics is in agreement with the mean-field kinetics model of Cates. By studying macroscopic relaxation phenomena such as the average micelle length evolution after a T-jump, the monomer diffusion, and the zero shear relaxation function, we confirm that the effective kinetic constants found are indeed the relevant parameters when macroscopic relaxation is coupled to the kinetics of micelles. For the non-equilibrium situation, we study the coupled effects of the shear flow and the scissionrecombination kinetics, on the structural and rheological properties of this micellar system. Our study is performed in semi-dilute and dynamically unentangled regime conditions. The explored parameter o range is chosen in order for the life time of the average size chain to remain shorter than its intrinsic (Rouse) longest relaxation time. Central to our analysis is the concept of dynamical unit of size A, the chain fiagrnent for which the life time TA and the Rouse time are equal. Shear thinning, chain orientation and bond stretching are found to depend upon the reduced shear rate P1\=y~A while the average micelle size is found to decrease with increasing shear rate, independently of the height of the barrier of the scission-recombination process

Ce mémoire de thèse s'intéresse au comportement du béton sous chargement triaxial sévère, et plus précisément, aux mécanismes de dégradation associés. Le cadre plus général de cette étude est l'amélioration de la tenue des infrastructures soumises à des chargements dynamiques de type impact, pouvant générer des contraintes moyennes dans le béton de l'ordre du GPa. Des essais quasi-statiques sont nécessaires pour caractériser le comportement d'un matériau dans un état de contraintes homogènes. Les travaux présentés dans ce manuscrit proposent d'investiguer les mécanismes de dégradation du béton lors d'essais triaxiaux. Pour cela, deux machines innovantes sont utilisées : une presse triaxiale de grande capacité (Giga) et un tomographe à rayons X. Le processus de réalisation des essais consiste à scanner des échantillons de béton à l'état vierge et après chaque cycle de chargement (charge-décharge) effectué avec la presse triaxiale. Les images issues de ces examens tomographiques permettent d'évaluer la modification de la mésostructure du béton (échelle où les granulats et les macropores apparaissent) selon différents trajets et niveaux de chargement. Deux échantillons ont été testés triaxialement selon ce processus. Le premier, à relativement faible confinement (50 MPa), montre que les mécanismes menant à la rupture sont fortement localisés et principalement initiés à l'échelle microscopique. Le deuxième, testé à un très fort confinement (650 MPa), présente un endommagement diffus avec également un mécanisme de rupture localisé. Un modèle mésoscopique, utilisant la méthode aux éléments discrets, est également présenté dans ce manuscrit. L'originalité de cette approche réside d'une part dans le caractère prédictif du modèle développé, le comportement du mortier, des granulats ou de leur interaction étant identifié a priori ou à partir d'essais réalisés sur mortier et roche. D'autre part, la méthode de construction du maillage éléments discrets s'appuie sur la segmentation 3D d'images tomographiques ce qui permet une comparaison des faciès de ruptures expérimentaux et numériques. Ce modèle aboutit à une bonne prédiction du comportement triaxial du béton, et semble être un outil prometteur d'étude des mécanismes d'endommagement et de l'influence de paramètres de composition du béton.

Ce mémoire de thèse s'intéresse au comportement du béton sous chargement triaxial sévère. L'étude est réalisée dans le contexte plus général de la compréhension du comportement du béton sous impact, ce qui induit des états de contraintes triaxiales de très haute intensité. Afin de reproduire des niveaux de contraintes très élevés avec des chemins de chargement bien contrôlés, des essais triaxiaux statiques sont réalisés sur des échantillons de béton au moyen d'une presse triaxiale de très grande capacité. Nous nous sommes intéressés à l'influence des granulats, qui occupent environ 40% du volume du béton. Plus précisément, les effets de la forme et de la composition chimique des granulats sont étudiés pour des compression triaxiale entre 0 et 650 MPa de confinement . On considère à la fois la réponse macroscopique et les modes de rupture. La forme des granulats semble peu influencer le comportement du béton, alors que les effets de leur composition sont assez importants, quelque soit le confinement. La deuxième partie de la thèse est consacrée à une modélisation numérique du béton. L'objectif principal de ce travail numérique est de développer un modèle très simple en termes de lois d'interaction et d'introduire l'hétérogénéité du béton au niveau mésoscopique. Un modèle mésoscopique du béton est ainsi établi dans le cadre des éléments discrets. On étudie l'influence des propriétés des granulats et de l'interface graulats/mortier sur les courbes contrainte-déformation et les mécanismes d'endommagement. D'autre part, l'influence des vides est étudiée en utilisant des échantillons mono-phasiques avec différentes caractéristiques des vides. Un comportement complexe non linéaire, résultant de la présence, des vides est observé. Par ailleurs, les résultats numériques complètent, et permettent de mieux comprendre, les résultats expérimentaux.

Ce mémoire de thèse s'intéresse au comportement du béton sous chargement triaxial sévère. L'étude est réalisée dans le contexte plus général de la compréhension du comportement du béton sous impact, ce qui induit des états de contraintes triaxiales de très haute intensité. Afin de reproduire des niveaux de contraintes très élevés avec des chemins de chargement bien contrôlés, des essais triaxiaux statiques sont réalisés sur des échantillons de béton au moyen d'une presse triaxiale de très grande capacité. Nous nous sommes intéressés à l'influence des granulats, qui occupent environ 40% du volume du béton. Plus précisément, les effets de la forme et de la composition chimique des granulats sont étudiés pour des compression triaxiale entre 0 et 650 MPa de confinement . On considère à la fois la réponse macroscopique et les modes de rupture. La forme des granulats semble peu influencer le comportement du béton, alors que les effets de leur composition sont assez importants, quelque soit le confinement. La deuxième partie de la thèse est consacrée à une modélisation numérique du béton. L'objectif principal de ce travail numérique est de développer un modèle très simple en termes de lois d'interaction et d'introduire l'hétérogénéité du béton au niveau mésoscopique. Un modèle mésoscopique du béton est ainsi établi dans le cadre des éléments discrets. On étudie l'influence des propriétés des granulats et de l'interface graulats/mortier sur les courbes contrainte-déformation et les mécanismes d'endommagement. D'autre part, l'influence des vides est étudiée en utilisant des échantillons mono-phasiques avec différentes caractéristiques des vides. Un comportement complexe non linéaire, résultant de la présence, des vides est observé. Par ailleurs, les résultats numériques complètent, et permettent de mieux comprendre, les résultats expérimentaux. / This PhD thesis focuses on identifying concrete behavior under high triaxial loading. The study is carried out within a more general context of understanding the behavior of concrete under impact, which induces very high-intensity triaxial stress states. In order to reproduce very high stress levels with well-controlled loading paths, static triaxial tests are carried out on concrete samples by means of a very high-capacity triaxial press. We are interested in the influence of coarse aggregates, which occupy approximately 40% of the concrete volume. More specifically, the effects of coarse aggregate shape and composition on concrete behavior for triaxial compression tests ranging from 0 to 650 MPa of confinement are investigated. Both macroscopic response and failure patterns are concerned. Coarse aggregate shape seems to have only a slight influence on concrete behavior while the effects of coarse aggregate composition are quite important over whole range of confining pressures. The second part of the thesis is devoted to a numerical modeling of concrete. The main objective of this numerical work is to develop a very simple model in terms of interaction laws and introduce concrete heterogeneity at the mesoscopic level. A mesoscopic model of concrete is thus established within the discrete element framework. The effects of coarse aggregate and aggregate/mortar interface properties on the stress-strain curves and damage mechanisms in a concrete sample are investigated. On the other hand, the influence of voids is studied using one-phase samples having different void properties. A complex, nonlinear behavior of a numerical specimen arising from the void structuring is observed. The numerical results complete and allow better understanding of experimental observations.

Cette thèse s’inscrit à mi-chemin entre l’étude de la déformabilité des structures tissées et la valorisation de la fibre de lin pour des applications dans le renforcement des matériaux composites. Le premier objectif de l’étude est de caractériser expérimentalement le comportement en flexion des mèches de différentes structures constituées de fibres de lin ainsi que des tissus de différentes armures. Les travaux ont abordé aussi des paramètres tels que l’humidité relative et la composition (des mèches comélées ou en lin pure). Le deuxième objectif des travaux est d’étudier le comportement en flexion des tissus en fonction du comportement en flexion des mèches. Cette partie a commencé par la modélisation géométrique des renforts tissés dans le but de suivre l’évolution de la section du tissu qui varie dans la direction de la flexion. La modélisation mésoscopique a permis de calculer analytiquement les propriétés géométriques du tissu en particulier son moment quadratique. Les résultats obtenus ont été utilisés dans la simulation de la flexion du tissu. L’étude a permis de voir jusqu’à quel point le comportement de la mèche et le moment quadratique du tissu pilotent le comportement en flexion du tissu. D’après ces travaux, le comportement en flexion du tissu semble être approché de façon satisfaisante sur toute la gamme de longueurs envisagées à partir de ces deux grandeurs sauf pour les forts taux d’humidité où d’autres phénomènes doivent être considérés. L’étude a souligné que la différence entre deux renforts testés expérimentalement peut être anticipée numériquement. Ainsi, le concepteur de tissus sera capable d’anticiper la rigidité expérimentale du tissu pour faire des tissages adaptés à la mise en forme du renfort. Une étude paramétrique de la flexion a été également réalisée dans le but de déduire les paramètres les plus influents sur lesquels il peut jouer. / This thesis is halfway between the study of the deformability of woven structures and the use of flax fibre as reinforcement of composite materials. The first aim of the study is the experimental characterization of the bending behaviour of tows with different structures made of flax fibres and fabrics with different weaves. Parameters such as relative humidity and the composition (100% flax and commingled tows) were also considered. The second aim of the study is to link the bending behaviour of the fabric to the bending behaviour of its constituent tows. This part starts with the geometric modelling of woven fabrics in order to follow the variation of its section in the bending direction. Mesoscopic modelling allows the analytical calculation of the geometric properties of the fabric in particular its moment of inertia. The results obtained were used in the simulation of the fabrics bending to see how far the behaviour depends on the tows bending behaviour and the moment of inertia. The bending behaviour of the fabric seems to be approached satisfactorily from these two factors. This is verified within the range of lengths considered except for high humidity (in this case, other phenomena must be considered). The study pointed out that the difference between two reinforcements tested experimentally can be predicted numerically. Thus, the fabrics designer will be able to anticipate the experimental bending stiffness of the fabric in order to adapt the weaving to the shape forming. A parametric study of the bending was also achieved in order to deduce the most influential parameters of the fabric for an appropriate weaving.

Le béton est un matériau très utilisé dans l'industrie de construction. Une limite essentielle à un usage de ce matériau est sa dégradation par écaillage lorsqu’il est exposé au feu. Le phénomène d'écaillage consiste en une éjection du béton sous forme d'écailles à la surface du béton exposée à l'incendie. La section de béton s'en trouve progressivement réduite. De plus, ce phénomène expose les armatures et peut conduire à une rupture prématurée de structures telles que les tunnels, les gratte-ciels, les centrales nucléaires etc. De nombreuses recherches ont été consacrées à la mise au point de méthodes de prévention de l'écaillage et à la détermination des paramètres qui ont une influence sur ce phénomène. Cependant, la physique qui contrôle l'écaillage n'est pas encore entièrement comprise. L'objectif principal de la thèse est de fournir une meilleure compréhension des mécanismes impliqués dans l’écaillage du béton en utilisant une approche numérique-expérimentale, i.e., la tomographie neutronique couplée à la modélisation numérique avancée à une échelle adéquate.Dans ce travail, les premières mesures 3D de la teneur en eau du béton (grandeur locale indispensable au suivi du processus de déshydratation potentiellement responsable de l’écaillage) soumis à un chargement thermique sévère ont été réalisées à l'aide de tomographies neutroniques rapides. Le suivi de la déshydratation rapide du béton a été possible en réalisant un scan 3D toutes les minutes grâce à la source neutrons de l'Institut Laue Langevin (leader mondial), à Grenoble, France. Cette vitesse d'acquisition est dix fois plus rapide que toute autre étude tomographique rapportée dans la littérature. Un dispositif, adapté à l'imagerie neutronique et aux essais à haute température, a été développé pour réaliser de telles expériences. L'influence de la taille des agrégats sur la distribution de l'humidité au sein de l'échantillon est présentée. Les résultats quantitatifs sur l'accumulation d'humidité derrière le front de déshydratation, connue sous le nom de 'moisture-clog" et considérée comme un des facteurs principaux engendrant un excès de pression, sont également présentés et discutésEn parallèle, un modèle thermo-hydro-mécanique (THM) entièrement couplé a été mis en œuvre sur le logiciel élément fini Cast3M afin d'étudier et prédire le comportement du béton à haute température. Le code nouvellement implémenté est remarquablement plus rapide (20-30 fois) que le code existant sur lequel il est basé. Une approche mésoscopique a été adaptée au modèle pour prendre en compte l'hétérogénéité du béton. D'abord, le modèle est appliqué à des expériences de la littérature, qui étu-dient les paramètres standards tels que la température, la pression du gaz et la perte de masse. En-suite, des profils d'humidité 1D obtenus à partir d'expériences de radiographie neutronique sont utili-sés pour vérifier et améliorer le modèle en termes de lois de comportement critiques telles que les courbes de déshydratation et de rétention d'eau. Enfin, le modèle est utilisé pour prédire la distribution d'humidité 3D mesurée dans ce travail de doctorat par tomographie neutronique. Entre autres, des simulations THM mésoscopiques sont effectuées pour étudier l'influence d'un agrégat sur le front de séchage. / Concrete has been extensively used in the construction industry as a building material. A major drawback of this material is its instability at high temperature, expressed in the form of violent or non-violent detachment of layers or pieces of concrete from the surface of a structural element. This phenomenon, known as fire spalling, can lead to the failure of concrete structures such as tunnels, high rise buildings, nuclear power-plants, underground parkings etc. because the reinforcement steel is directly exposed to high temperature and the designed cross section of the concrete elements (e.g., columns, beams, slabs) is reduced. A lot of research has been dedicated on developing preventing methods for spalling and also on determining the parameters that have an influence on it. However, the physics behind this phenomenon is not yet fully understood.In this doctoral, the first 3D measurements of moisture content in heated concrete, which is believed to be one of the processes directly related to spalling, have been performed using in-situ neutron tomography. In order to follow the fast dehydration process of concrete, one 3D scan (containing 500 radiographs) per minute was captured thanks to the world leading flux at the Institute Laue Langevin (ILL) in Grenoble France. This acquisition speed, which is ten times faster than any other experiment reported in the literature, was sufficient to follow the dehydration process. A dedicated setup, adapted to neutron imaging and high temperature, has been developed for performing such kind of experiments. Concrete samples with different aggregate size have been tested. Quantitative analysis showing the effect of the aggregate size on the moisture distribution is presented. Results on the moisture accumulation behind the drying front, known as the moisture-clog, are also presented and discussed.In parallel, a numerically-efficient coupled thermo-hydro-mechanical (THM) model has been implemented in the finite element software Cast3M for understanding and predicting the complex behavior of concrete at high temperature in the context of spalling. The newly implemented code is remarkably faster (20-30 times) than an existing one, on which it is based. A mesoscopic approach has been adapted to the model for taking into account the heterogeneity of concrete. First the model is applied to experiments from literature monitoring standard parameters such as temperature, gas pressure and mass loss. Then, 1D moisture profiles obtained from neutron radiography experiments are used for verifying and improving the model in terms of some critical constitutive laws such as dehydration and water retention curves. Finally, the model is employed for predicting the 3D moisture distribution measured in this doctoral work via neutron tomography. Among others, mesoscopic THM simulations are performed for investigating the influence of an aggregate on the drying front.

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.

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.

La simulation à l'échelle mésoscopique de la déformation des renforts composites fournit des informations importantes. En particulier, elle donne la direction et la densité de fibres qui conditionne la perméabilité du renfort textile et les propriétés mécaniques du composite final. Ces analyses mésoscopiques par éléments finis dépendent fortement de la qualité de la géométrie initiale du modèle. Certains logiciels ont été développés pour décrire ces géométries de renforts composites. Mais, les géométries obtenues impliquent une simplification (notamment dans la section transversale de mèche) qui peut perturber le calcul de déformation du renfort. Le présent travail présente une méthode directe utilisant la microtomographie à rayon X pour générer des modèles éléments finis, basée sur la géométrie réelle de l'armure textile. Le modèle EF peut être obtenu pour tout type de renfort, plus ou moins complexe. Les problèmes d’interpénétrations de mèches sont évités. Ces modèles sont utilisés avec deux lois de comportement : une loi hypoélastique et une loi hyperélastique. Les propriétés de chacune d'entre elles, ainsi que les grandeurs caractéristiques nécessaires à leur implémentation dans le logiciel ABAQUS sont développées. Enfin, une identification des paramètres matériau à l'aide d'une méthode inverse est proposée. Les résultats obtenus pour les simulations de mise en forme sont comparés avec les résultats expérimentaux et montrent une bonne correspondance entre les deux. / The simulation at meso-scale of textile composite reinforcement deformation provides important information. In particular, it gives the direction and density of the fibres that condition the permeability of the textile reinforcement and the mechanical properties of the final composite. These meso FE analyses are highly dependent on the quality of the initial geometry of the model. Some software have been developed to describe composite reinforcement geometries. The obtained geometries imply simplification that can disrupt the reinforcement deformation computation. The present work presents a direct method using computed microtomography to determine finite element models based on the real geometry of the textile reinforcement. The FE model is obtained for any specificity or variability of the textile reinforcement, more or less complex. The yarns interpenetration problems are avoided. These models are used with two constitutive laws : a hypoelastic law and a hyperelastic one. An analysis of their properties is presented and their implementation in the software ABAQUS is detailed. Finally, an identification method is presented and the results of forming simulations are compared to experimental tests, which shows a good fit between the both.