Les CMC oxyde/oxyde sont de bons candidats pour des applications thermostructurales. Le comportement mécanique et les mécanismes d’endommagement de deux composites alumine/alumine à renforts tissés bi- et tridimensionnels ont été étudiés et comparés. La microstructure de ces CMC à matrice faible a été caractérisée à partir de porosimétrie et de CND, tel que thermographie IR, scan ultrasonore et tomographie X, ce qui a permis de mettre en évidence la présence de défauts initiaux. Le comportement mécanique en traction, ainsi qu’en compression dansle cas du CMC à renfort bidimensionnel, dans la direction des fibres ainsi que dans la direction ±45°, aété étudié à température ambiante. Afin d’exploiter pleinement ces essais, nous avons eu recours à plusieurs méthodes d’extensométrie et de suivi d’endommagement, telles que la thermographie IR et l’émission acoustique. Les propriétés mécaniques à rupture ainsi que le module de Young du CMC à renfort bidimensionnel développé à l’Onera se sont avérées supérieures à celles disponibles dans la littérature. Les mécanismes d’endommagement des matériaux ont été déterminés à partir d’observations post mortem au MEB et d’essais in situ dans un MEB, ce qui a permis d’évaluer la nocivité des défauts initiaux. Enfin, l’étude du comportement mécanique de ces composites a permisde proposer un modèle d’endommagement tridimensionnel qui permettra de poursuivre le développement de ces matériaux grâce à du calcul de structure. A l’issue de cette thèse, des pistes d’amélioration des procédés d’élaboration et de choix d’instrumentation à utiliser pour les futures études, notamment en ce qui concerne le suivi d’endommagement, ont également été proposées. / Oxide/oxide CMCs are good candidates for thermostructural applications. Themechanical behaviour and damage mechanisms of two alumina/alumina composites with two andthree dimensional woven reinforcements were studied and compared. The microstructure of theseweak matrix CMCs was characterized by porosimetry and NDT methods, such as IR thermography,ultrasound scanning and X-ray tomography, which highlighted initial defects. The mechanicalbehaviour was studied through tensile tests, as well as compression tests in the case of the twodimensionalreinforced CMC. These tests were conducted at room temperature, in the fibres directionsand in the ±45° direction. In order to fully exploit these tests, several extensometry and damagemonitoring methods, such as IR thermography and acoustic emission, were used. Young’s moduli andmaximum stresses and strains of the two-dimensional reinforced CMC developed at Onera appearedto be higher than those available in the literature. The damage mechanisms of the materials weredetermined by post mortem SEM observations and in situ testing in a SEM, which made it possible toassess the nocivity of initial defects. Studying the mechanical behaviour of these composites finallyenabled the development of a three-dimensional damage model that will facilitate the furtherdevelopment of such materials, through finite element analysis. Finally, some improvements regardingthe manufacturing processes and the instrumentation for damage monitoring were suggested forfuture studies.

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 aubes de turboréacteurs sont aujourd'hui sujets à de nouvelles technologiques novatrices : les matériaux composites tissées 3D. Ces modifications matériaux permettent des réductions de masse significatives, et conduisent à des performances énergétiques et environnementales accrues. Ces structures tissées possèdent généralement des architectures complexes et requièrent une bonne caractérisation mécanique. Des outils de modélisation textile et d'homogénéisation ont été développés, la géométrie et les propriétés physiques de l'architecture tissée y sont analysées. Toutefois, il est important de prévoir les endommagements engendrés dans ces matériaux au niveau des zones de contact soumises à du fretting. En effet, les pieds d'aubes de turboréacteurs font face à des sollicitations de type fretting. Il en résulte deux types d'endommagements : l'amorçage et propagation de fissures et l'usure des surfaces en contact. Quantifier ces endommagements requiert une analyse fine du contact à partir d'un outil de calcul à la fois robuste et rapide. Pour mener à bien cet objectif, un code semi-analytique permettant aujourd'hui de traiter des problèmes hétérogènes de contacts élasto-plastiques est développé. La méthode d'inclusion équivalente proposée au sens d'Eshelby est utilisée dans le solveur de contact pour décrire l'effet des hétérogénéités (cavités, inclusions, fibres ou torons). L'un des corps en contact peut contenir de multiples hétérogénéités de formes parallélépipédiques, ellipsoïdales, et ses formes dégénérées (ellipsoïde oblate, ellipsoïde prolate, sphère, cylindre, disque plat,...). Cette méthode est modifiée et améliorée afin de prendre en compte les influences mutuelles entre les inclusions hétérogènes voisines, et la décohésion à l'interface hétérogénéité/matrice. Un premier couplage avec un modèle de contact aube/disque Eléments Finis permet de prendre en compte les effets de structures. Cette méthode est analogue à un zoom structural lorsqu'un maillage fin est nécessaire dans des zones de contact à fort gradient de contraintes soumises à un endommagement de type fretting. Un second couplage avec le logiciel WiseTex permet de décrire la géométrie réelle du tissage du composite, les propriétés matériaux des fibres et de la matrice. Une segmentation fine du modèle permet de discrétiser les mèches en de multiples hétérogénéités ellipsoïdales équivalentes et d'appliquer la méthode d'inclusion équivalente modifiée. A partir de cette mise en donnée d'un matériau composite tissé revêtu, des calculs de contact tridimensionnel multiéchelle en situation de fretting sont réalisés avec succès.

Ce travail de thèse est consacré à l'étude des renforts textiles techniques 2D à l'échelle mésoscopique. La méthode des éléments finis est utilisée pour résoudre un problème aux limites, fortement non linéaire, dans le domaine du renfort fibreux. Deux nouveaux modes de déformations pour le comportement des mèches de renforts ont été développés et caractérisés. Le premier mode est la compression longitudinale, qui permet de traduire la faible rigidité des mèches lorsqu'elles subissent une dilatation négative dans le sens des fibres. La relation conflictuelle sur le plan numérique entre la rigidité en tension, très forte, et la rigidité en compression longitudinale, très faible, peut être résolue via trois méthodes : réduction du pas de temps critique, addition de la contribution en tension ou avec une nouvelle stratégie pour l'actualisation du champ de contrainte. Le second mode de déformation est la dilatation transversale des mèches considérée comme conséquence directe de la compression longitudinale. Ce phénomène d'expansion de matière dans les directions transverses peut être observé avec un essai de compression longitudinale in-situ sous tomographie X et est communément ignorer. Un protocole expérimental a été mise en place pour mesurer cette dilatation transversale des mèches et un coefficient de Poisson a pu être caractérisé par identification inverse. Une campagne expérimentale a permis d'identifier les paramètres matériau du modèle mésoscopique et les résultats de simulations sont comparés aux images issues d'essai mécanique in situ sous tomographe. / This thesis is devoted to the mesoscopic study on the performance of textile reinforcements. F.E. simulation is carried out on a mesoscale model for the fibrous material, based on which two kinds of new deformation modes are developed. The first one is a longitudinal compression mode, which is used to reflect the small stiffness when the yarn is compressed longitudinally. The incompatibility problem between the small longitudinal compression stiffness and the large tension stiffness are solved by three different strategies: constraining the critical step time, adding the nonlinear tension part, or using a new strategy to update the stress. The second one is transverse expansion mode that could reflect the influence from longitudinal deformation to transverse deformation. This deformation could be found in tomography view but was ignored by the former researches. An experiment is designed to measure the expansion magnitude, and the geometrical inverse fitting process is applied to measure the value of the longitudinal-transverse Poisson ratio. The parameters of the mesoscale model are measured by a series of mechanical experiments and the simulation results are verified by the tomography methodology.

Les aubes de turboréacteurs sont aujourd’hui sujets à de nouvelles technologiques novatrices : les matériaux composites tissées 3D. Ces modifications matériaux permettent des réductions de masse significatives, et conduisent à des performances énergétiques et environnementales accrues. Ces structures tissées possèdent généralement des architectures complexes et requièrent une bonne caractérisation mécanique. Des outils de modélisation textile et d’homogénéisation ont été développés, la géométrie et les propriétés physiques de l’architecture tissée y sont analysées. Toutefois, il est important de prévoir les endommagements engendrés dans ces matériaux au niveau des zones de contact soumises à du fretting. En effet, les pieds d’aubes de turboréacteurs font face à des sollicitations de type fretting. Il en résulte deux types d’endommagements : l’amorçage et propagation de fissures et l’usure des surfaces en contact. Quantifier ces endommagements requiert une analyse fine du contact à partir d’un outil de calcul à la fois robuste et rapide. Pour mener à bien cet objectif, un code semi-analytique permettant aujourd’hui de traiter des problèmes hétérogènes de contacts élasto-plastiques est développé. La méthode d’inclusion équivalente proposée au sens d’Eshelby est utilisée dans le solveur de contact pour décrire l’effet des hétérogénéités (cavités, inclusions, fibres ou torons). L’un des corps en contact peut contenir de multiples hétérogénéités de formes parallélépipédiques, ellipsoïdales, et ses formes dégénérées (ellipsoïde oblate, ellipsoïde prolate, sphère, cylindre, disque plat,...). Cette méthode est modifiée et améliorée afin de prendre en compte les influences mutuelles entre les inclusions hétérogènes voisines, et la décohésion à l’interface hétérogénéité/matrice. Un premier couplage avec un modèle de contact aube/disque Eléments Finis permet de prendre en compte les effets de structures. Cette méthode est analogue à un zoom structural lorsqu’un maillage fin est nécessaire dans des zones de contact à fort gradient de contraintes soumises à un endommagement de type fretting. Un second couplage avec le logiciel WiseTex permet de décrire la géométrie réelle du tissage du composite, les propriétés matériaux des fibres et de la matrice. Une segmentation fine du modèle permet de discrétiser les mèches en de multiples hétérogénéités ellipsoïdales équivalentes et d’appliquer la méthode d’inclusion équivalente modifiée. A partir de cette mise en donnée d’un matériau composite tissé revêtu, des calculs de contact tridimensionnel multiéchelle en situation de fretting sont réalisés avec succès. / New baseline turbofan engines feature advanced blade technology made of 3D woven composites for a significant reduction in weight and an increase in energy and environmental performance. Woven structures generally have complex architectures wich require high level of mechanical model. Modeling textile and homogenization tools have been developed to provide precisely the geometry and physical properties of the woven architecture. However, it is important to predict the damage mode under contact zone which are subjected to fretting. Indeed, the blade roots of turbofan engines are damaged by fretting phenomenon. Fretting modes generate two kinds of damage : (i) the initiation and propagation cracks, (ii) wear on contact surfaces. Quantify these two kinds of damage requires a detailed contact analysis from a robust and fast contact solver. In order to reach this goal, a semi-analytical solver allows to solve heterogeneous elasto-plastic contact problems. Equivalent inclusion method in the sense of Eshelby allows to describe accurately the effect of inhomogeneities (cavities, inclusions, fibers or strands). Only one of contact bodies contains multiple heterogeneous of cuboïdal and ellipsoidal shapes, and their degenerated forms (oblate spheroid, prolate ellipsoid, sphere, cylinder, flat disk,...). This method is modified and improved in order to take into account the mutual influence between neighboring heterogeneous inclusions and decohesion at the interface heterogeneity/matrix. A first coupling with a finite element model of blade/disk contact allows to take into account the effects of structure. This method is analogous to a structural zoom when structural fine mesh is necessary within areas of high contact stress gradient subjected to fretting. A second coupling with the software WiseTex allows to describe the geometry of the actual weaving of the composite, the material properties of the fibers and the matrix. A fine segmentation of the numerical model allows to discretize reinforcements to multiple equivalent ellipsoidal heterogeneities and to apply the modified equivalent inclusion method. Three-dimensional simulation of frictional contact model are successfully completed from these data of a coating woven composite.

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.

Le procédé de thermo-estampage est une voie intéressante pour la production en grande série de pièces composites. Cependant ce procédé est complexe à maitriser et simuler, en raison des phénomènes multi-physiques mis en jeu (déformation textiles, choc thermique, frottements...) ce qui engendre des campagnes par essai-erreurs qui peuvent être très coûteuses. Cette étude s'intéresse à la mesure et à la caractérisation du procédé d'estampage et des matériaux utilisés afin de simuler le refroidissement de la matrice et la mise en forme du textile. Des défauts récurrents sont évoqués ainsi que des solutions industrielles afin de les résoudre. La simulation thermique permet de déterminer le temps de consolidation nécessaire afin d'optimiser les temps de production en fonction des matériaux et de leur épaisseur. La simulation de la mise en forme textile permet de prédire la faisabilité d'une pièce et l'orientation des fibres afin de définir au mieux les pièces suivant les cas de charges statiques et dynamiques. Les apports de ce travail sont les suivants : la réalisation de mesures thermiques du flan durant un estampage et du choc thermique en surface du stratifié, la réalisation d'une méthodologie efficiente pour analyser la faisabilité d'une pièce complexe dans un contexte industriel grâce à des simulations de mise en forme couplées à des essais expérimentaux. Enfin une méthode d'analyse du comportement en cisaillement plan, pour des renforts dont les fils de chaîne et de trames ne sont pas orthogonaux est proposée. / The thermo-stamping process is a promising way for the mass production of the composite parts. However this process is complex to master and simulate due to the multi-physics background (textile deformation, thermal shock, rubbing...) and trial and error tests campaigns can be expensive. This study focuses on the measurement and assessment of the process and materials behaviour, to simulate the cooling down of the matrix and the forming of the woven. Typical defects are mentioned as well as associated industrial solutions to solve them. The simulation makes it possible to determine the consolidation time necessary in order to optimize the manufacturing time as a function of the material used and of its thickness. Furthermore the forming simulation shows the feasibility of the part and the fibre orientation to design the product for the static and crash cases. The contributions of this work are the following: thermal measurements of the pre-consolidated plate during the stamping phase and the thermal chock at the surface of the composites, an efficient method to analyse the feasibility of a complex shape in an industrial context is proposed through forming simulations compared to the experiment. Finally, a methodology for the analysis of the in-plane shearing behaviour of a woven fabric with non-orthogonal warp and weft yarn is proposed.

Le calcul de la durée de vie des Composites tissés à Matrice Céramique (CMC) nécessite de déterminer l’évolution de la densité de fissures dans le matériau(pouvant atteindre 10 mm-1). Afin de les représenter finement on se propose de travailler à l’échelle mésoscopique. Les méthodes de type Embedded Finite Element (EFEM) nous ont paru être les plus adaptées au problème. Elles permettent une représentation discrète des fissures sans introduire de degrés de liberté additionnels.Notre choix s’est porté sur une EFEM s’affranchissant d’itérations élémentaires et appelée Augmented Finite Element Method (AFEM). Une variante d’AFEM, palliant des lacunes de la méthode originale, a été développée. Nous avons démontré que,sous certaines conditions, AFEM et la méthode des éléments finis classique (FEM) étaient équivalentes. Nous avons ensuite comparé la précision d’AFEM et de FEM pour représenter des discontinuités fortes et faibles. Les travaux de thèse se concluent par des exemples d’application de la méthode aux CMC. / Computing the lifetime of woven Ceramic Matrix Composites (CMC) requires evaluating the crack density in the material (which can reach 10 mm-1). Numerical simulations at the mesoscopic scale are needed to precisely estimate it. Embedded Finite Element Methods (EFEM) seem to be the most appropriate to do so. They allow for a discrete representation of cracks with no additional degrees of freedom.We chose to work with an EFEM free from local iterations named the Augmented Finite Element Method (AFEM). Improvements over the original AFEM have been proposed. We also demonstrated that, under one hypothesis, the AFEM and the classical Finite Element Method (FEM) are fully equivalent. We then compare the accuracy of the AFEM and the classical FEM to represent weak and strong discontinuities. Finally, some examples of application of AFEM to CMC are given.

Ce travail porte sur le développement de modèles micromécaniques pour le calcul de la réponse homogénéisée de matériaux architecturés, en particulier des matériaux se présentant sous forme de treillis répétitifs. Les matériaux architecturés et micro-architecturés couvrent un domaine très large de de propriétés mécaniques, selon la connectivité nodale, la disposition géométrique des éléments structuraux, leurs propriétés mécaniques, et l'existence d'une possible hiérarchie structurale. L'objectif principal de la thèse est la prise en compte des nonlinéarités géométriques résultant des évolutions importantes de la géométrie initiale du treillis, causée par une rigidité de flexion des éléments structuraux faible en regard de leur rigidité en extension. La méthode dite d'homogénéisation discrète est développée pour prendre en compte les non linéarités géométriques pour des treillis quais périodiques; des schémas incrémentaux sont construits qui reposent sur la résolution incrémentale et séquentielle des problèmes de localisation - homogénéisation posés sur une cellule de base identifiée, soumise à un chargement contrôlé en déformation. Le milieu continu effectif obtenu est en général un milieu micropolaire anisotrope, dont les propriétés effectives reflètent la disposition des éléments structuraux et leurs propriétés mécaniques. La réponse non affine des treillis conduit à des effets de taille qui sont pris en compte soit par un enrichissement de la cinématique par des variables de microrotation ou par la prise en compte des seconds gradients du déplacement. La construction de milieux effectifs du second gradient est faite dans un formalisme de petites perturbations. Il est montré que ces deux types de milieu effectif sont complémentaires en raison de l'analogie existant lors de la construction théorique des réponses homogénéisées, et par le fait qu'ils fournissent des longueurs internes en extension, flexion et torsion. Des applications à des structures tissées et des membranes biologiques décrites comme des réseaux de filaments quais-périodiques ont été faites. Les réponses homogénéisées obtenues sont validées par des comparaisons avec des simulations par éléments finis réalisées sur un volume élémentaire représentatif de la structure. Les schémas d'homogénéisation ont été implémentés dans un code de calcul dédié, alimenté par un fichier de données d'entrée de la géométrie du treillis et de ses propriétés mécaniques. Les modèles micromécaniques développés laissent envisager du fait de leur caractère prédictif la conception de nouveaux matériaux architecturés permettant d'élargir les frontières de l'espace 'matériaux-propriétés' / The present thesis deals with the development of micromechanical schemes for the computation of the homogenized response of architectured materials, focusing on periodical lattice materials. Architectured and micro-architectured materials cover a wide range of mechanical properties according to the nodal connectivity, geometrical arrangement of the structural elements, their moduli, and a possible structural hierarchy. The principal objective of the thesis is the consideration of geometrical nonlinearities accounting for the large changes of the initial lattice geometry, due to the small bending stiffness of the structural elements, in comparison to their tensile rigidity. The so-called discrete homogenization method is extended to the geometrically nonlinear setting for periodical lattices; incremental schemes are constructed based on a staggered localization-homogenization computation of the lattice response over a repetitive unit cell submitted to a controlled deformation loading. The obtained effective medium is a micropolar anisotropic continuum, the effective properties of which accounting for the geometrical arrangement of the structural elements within the lattice and their mechanical properties. The non affine response of the lattice leads to possible size effects which can be captured by an enrichment of the classical Cauchy continuum either by adding rotational degrees of freedom as for the micropolar effective continuum, or by considering second order gradients of the displacement field. Both strategies are followed in this work, the construction of second order grade continua by discrete homogenization being done in a small perturbations framework. We show that both strategies for the enrichment of the effective continuum are complementary due to the existing analogy in the construction of the micropolar and second order grade continua by homogenization. The combination of both schemes further delivers tension, bending and torsion internal lengths, which reflect the lattice topology and the mechanical properties of its structural elements. Applications to textiles and biological membranes described as quasi periodical networks of filaments are considered. The computed effective response is validated by comparison with FE simulations performed over a representative unit cell of the lattice. The homogenization schemes have been implemented in a dedicated code written in combined symbolic and numerical language, and using as an input the lattice geometry and microstructural mechanical properties. The developed predictive micromechanical schemes offer a design tool to conceive new architectured materials to expand the boundaries of the 'material-property' space

Au cours des dernières décennies, le développement de nouveaux matériaux a progressé pour les applications liées à la mécanique. De nouvelles générations de composites ont été développées, qui peut offrir des avantages par rapport aux tapis unidirectionnels renforcés de fibres couramment utilisés les matériaux prennent alors le nom de woven fabrics. Le comportement de ce matériau est fortement influencé par la micro-structure du matériau. Dans la thèse, les modèles mécaniques et les schémas numériques capables de modéliser les comportement des tissus et des matériaux de réseau généraux ont été développés. Le modèle prend en compte la micro-structure au moyen d'une technique d'homogénéisation. Les fibres dans le réseau ont été traités comme des micro-poutres, ayant une rigidité à la fois en extension et en flexion, avec différents types de connexions. La procédure développée a été appliquée pour obtenir les modèles mécaniques homogénéisés pour certains types de réseaux de fibres biaxiaux et quadriaxiaux, simulant soit des réseaux de fibres (en ce cas a été supposé parmi les fibres) ou des tissus avec une interaction négligeable entre les faisceaux de fibres et en empêchant tout glissement relatif (dans ce cas, les connexions ont été simulés au moyen de pivots). Différentes géométries ont été analysées, y compris la cas dans lesquels les fibres ne sont pas orthogonales. On obtient généralement un premier milieu à gradient mais, dans certains cas, la procédure d'homogénéisation lui-même indique qu'un continuum d'ordre supérieur est mieux adapté pour représenter la déformation de la micro-structure. Des résultats spéciaux ont été obtenus dans le cas de fibres reliées par pivots. Dans ce cas, un matériau orthotrope à module de cisaillement nul a été obtenu. Un tel matériau a un tenseur constitutif elliptique, il peut donc conduire à des concentrations de contrainte. Cependant, il a été montré que certaines considérations sur le comportement physique de tels réseaux indiqué que les termes d'ordre supérieur inclus dans l'expansion des forces internes et des déformations, de sorte qu'un matériau de gradient de déformation a été obtenu. Les résultats obtenus peuvent être utilisés pour la conception de matériaux spécifiques nécessitant des propriétés. Bien que le modèle de référence soit un matériau de réseau, les résultats obtenus peuvent être appliqué à d'autres types similaires de microstructures, comme des matériaux pantographiques, des micro-dispositifs composé de micro-poutres, etc. Ils étaient limités à la gamme d'élasticité linéaire, qui est petite déformation et comportement élastique linéaire. Ensuite, les simulations numériques ont été axées sur les tests d'extension et les tests de biais. Le obtenu configurations déformées sont conformes aux tests expérimentaux de la littérature, tant pour tissus équilibrés et non équilibrés. De plus, une comparaison entre les premier et deuxième gradients des prédictions numériques ont été effectuées. Il a été observé que les prédictions de deuxième gradient mieux simuler les preuves expérimentales. / In the past decades there has been an impressive progress in the development of new materials for mechanical related applications. New generations of composites have been developed, that can offer advantages over the unidirectional fibre-reinforced mats commonly used then materials take the name of woven fabrics. The behaviour of this material is strongly influenced by the micro-structure of the material. In the thesis mechanical models and a numerical scheme able to model the mechanical behaviour of woven fabrics and general network materials have been developed. The model takes in to account the micro-structure by means of a homogenization technique. The fibres in the network have been treated like microbeams, having both extensional and bending stiffness, with different types of connection, according to the pattern and detail of the network. The developed procedure was applied for obtaining the homogenized mechanical models for some types of biaxial and quadriaxial networks of fibres, simulating either fibre nets (in this case rigid connection were assumed among the fibres) or tissues with negligible interaction between the fibre bundles, and with relative sliding prevented (in this case the connections were simulated by means of pivots). Different geometries were analysed, including the cases in which the fibres are not orthogonal. A first gradient medium is usually obtained but, in some cases, the homogenization procedure itself indicates that a higher order continuum is better fit to represent the deformation of the micro-structure. Special results were obtained for the case of fibres connected by pivots. In this cases an orthotropic material with zero shear modulus was obtained. Such a material has a not elliptic constitutive tensor, thus it can lead to strain concentrations. However, it was shown that some considerations about the physical behaviour of such networks indicated that higher order terms had to be included in the expansion of the internal forces and deformations, so that a strain gradient material was obtained. The results obtained can be used for the design of specific materials requiring ad-hoc properties. Although the reference model is a network material, the results obtained can be applied to other similar kinds of microstructures, like pantographic materials, micro devices composed by microbeams etc. They have been limited at the range of linear elasticity, that is small deformation and linear elastic behaviour. Then, numerical simulations were focused on extension tests and bias tests. The obtained deformed configurations are consistent with the literature experimental tests, both for balanced and unbalanced tissues. Moreover, a comparison between first and second gradient numerical predictions was performed. It was observed that second gradient predictions better simulate the experimental evidences.