Nous nous sommes intéressés dans ce travail à l’étude du vieillissement physique d’un non-tissé fabriqué à l’aide de fibres textiles issues du recyclage de bouteilles plastiques selon le procédé par voie sèche. L’application de ce non-tissé est l’isolation thermique des bâtiments français. Après avoir mis en évidence la présence de deux réseaux de fibres, nous nous sommes intéressés au phénomène de relaxation structurale. Le vieillissement physique du non-tissé étant régi par l’enrobage d’une des fibres utilisées dans la fabrication de celui-ci, nous avons observé l’évolution de l’enthalpie de relaxation en fonction de différents couples temps de vieillissement, écart de température par rapport à sa température de transition vitreuse. Les résultats ont mis en évidence un état d’avancement du mécanisme de relaxation quasi-nul pour des températures très éloignées de la température de transition vitreuse de l’enrobage. L’évolution de cette même enthalpie de relaxation a été étudiée pour différents non-tissés issus de différents lots de fabrication des cinq dernières années. Il a été mis en évidence un mécanisme de relaxation lent dans les conditions expérimentales nous rassurant quant à la stabilité du produit fini dans le temps. / We focused on physical ageing of a non-woven produced by recycled fibers from plastic bottles. This non-woven is produced by drylaid-textile process and is used for building thermal insulation. Two networks are visible and structural relaxation has been investigated through the study of evolution of one fiber skin part with different couple time, temperature. Physical ageing of polymers can be followed by structural relaxation process. We showed that structural relaxation process is very long for temperatures of use far from glass transition temperature. Same results have been shown on the global non-woven stored in non-insulated box governed by Normandy weather. The combination of the two previous results reinforce our basic idea ; the insulator properties variation will be stable for long time only if temperature of use is far from glass transition temperature.

Face aux défis environnementaux actuels, les industriels ont mis en œuvre de nouveaux matériaux recyclables et permettant une réduction significative de la masse. Le développement de la résine thermoplastique Elium par ARKEMA s’inscrit dans cette problématique. L’utilisation de cette résine pour la fabrication de pièces composites qui peuvent être sujettes à des dommages d’impact, nécessite au préalable des études, dans le but de comprendre leurs mécanismes de ruine sous ce type de sollicitation. Ainsi, la présente thèse propose une contribution à l’analyse multi-échelle de la tenue à l’impact des composites stratifiés à base de la résine Elium. Une étude expérimentale préliminaire a permis de confirmer la meilleure résistance à l’impact des composites à matrice Elium acrylique, comparativement à celles des composites thermodurcissables conventionnels. Ensuite, les performances à l’impact des composites stratifiés ont été améliorées par l’introduction de copolymères à blocs dans la matrice. Ces derniers sont capables de former des micelles de tailles nanométriques et ainsi d’améliorer la ténacité de la matrice acrylique. Les effets de l’énergie d’impact, de la température et de la composition en nanocharges sur la réponse du matériau composite ont été analysés. Afin de proposer un outil d’aide à la prédiction de la réponse à l’impact des matériaux fibres de verre/Acrylique, deux stratégies de modélisation ont été retenues. La première modélisation (macroscopique) considère le pli tissé du stratifié comme un matériau homogène tandis que la seconde (mésoscopique) utilise une description géométrique de l’ondulation et de l’entrecroisement des torons noyés dans la résine Elium. Ces deux modèles considèrent des zones cohésives à l’interface entre les plis adjacents pour simuler le délaminage interlaminaire. Des essais de délaminage (expérimentaux et numériques) ont permis d’alimenter le modèle d’endommagement de l’interface interplis. D’autre part, des essais de caractérisation du comportement mécanique et de l’endommagement du matériau couplés à l’homogénéisation multi-échelle des matériaux par la Mécanique du Génome de Structure ont permis d’identifier les paramètres du modèle macroscopique. A l’échelle mésoscopique, le modèle géométrique a été réalisé grâce au logiciel Texgen. Ce logiciel permet d’obtenir une description approchée mais réaliste de l’ondulation des torons de fibres. La même description a servi à l’homogénéisation numérique multi-échelle des stratifiés étudiés. La simulation numérique de l’impact basse vitesse a été effectuée au moyen du logiciel d’éléments finis ABAQUS/Explicit. Les modèles de comportement du matériau ont été implémentés via la routine utilisateur VUMAT. Les résultats obtenus offrent une bonne corrélation avec les données expérimentales / In the race for light materials able of meeting modern environmental challenges, an acrylic resin (Elium) has been developed. Elium is a thermoplastic resin able to replace thermosetting matrices, which are widespread nowadays in the industrial world. The present study aims to evaluate the impact resistance and to understand the failure mechanisms of composite laminates based on acrylic matrix under impact loading. We provide a contribution to the multiscale analysis of the impact resistance of laminated composite.First, the impact resistance and the damage tolerance of the acrylic resin based composites were compared with those of conventional composites. Then, the impact performance of the laminated composites has been enhanced by adding copolymer blocks to the liquid acrylic resin. These copolymers are able to form micelles of nanometer sizes, which lead to the improvement of both the acrylic matrix fracture toughness and the impact resistance. The effects of the impact energy, temperature, and composition in nano-copolymers have also been investigated.In order to provide a numerical tool for the prediction of the impact response of the glass fiber/Acrylic laminates, two strategies have been analyzed. The first one, performed at the macroscopic scale, considers the woven ply of the laminate as homogeneous material, and the second one (at the mesoscopic scale), deals with a realistic geometrical description of the yarns undulation. Both models use cohesive zones at the interface between the adjacent plies, to simulate the delamination. For this purpose, experimental and numerical delamination tests were performed to feed the inter-ply damage model. Mechanical tests for material characterization were also performed on specimens in order to identify the ply-damage model parameters. The Mechanics of Structure Genome (MSG) and a finite element based micromechanics approaches were then conducted to evaluate the effective thermomechanical properties of the yarns and the plain woven composite laminate. The realistic topological and morphological textures of the composite were accounted through Texgen software. These numerical impact simulations were performed using the finite element software ABAQUS/Explicit. Both models were implemented through a user material subroutine VUMAT. The obtained results appear in a good agreement with the experimental data and confirm the relevance of the proposed approach.

Ce travail de thèse a pour objectif de comprendre la fragilité des composites à matrice polymère (PVC) contenant des particules de bois. Dans un premier temps, il s’agit d’examiner le rôle de chacun des constituants : matrice et charge, sur les propriétés mécaniques du composite. Cette étude a montré que l’influence des paramètres relatifs à la microstructure de la matrice PVC est minoritaire face à celle des paramètres de la farine de bois. Pour mettre en évidence le rôle des particules de bois, une méthode d’analyse de leur dispersion et leur orientation par une technique 3D non destructive, a été employée : la tomographie aux rayons X. Un traitement spécifique des données a permis de quantifier la taille, l’orientation et la connectivité des particules. Ainsi, on a montré que l’existence d’agglomérats et/ou d’un réseau percolant fragilise fortement le composite ce qui explique la chute de la résistance au choc lorsque le taux de farine de bois augmente. Pour optimiser le compromis fragilité/ductilité, différentes solutions ont été envisagées. Dans un deuxième temps, nos travaux s'intéressent à l'évolution de l'endommagement de matériaux composites possédant différents taux de farine de bois. Cette évolution est mesurée à l'aide de deux techniques : l'émission acoustique et la tomographie aux rayons X in-situ. La décohésion à l’interface entre la matrice PVC et les particules de bois est le mécanisme initiant l’endommagement. La transition ductile/fragile des composites s’explique par la coalescence de ces cavités. L’ajout de modifiant choc ne réduit pas suffisamment ce mécanisme. Les agents de couplage sont donc la voie à explorer afin d’améliorer l’interface matrice/particules et remédier au défaut d’adhésion qui pénalise ces composites. / The aim of the present study is to understand the origin of the brittleness of PVC/wood flour composites (WPC). The influence of the matrix and of the filler is investigated: fiber parameters are key factors that significantly affect the mechanical properties of WPC. Fiber dispersion, orientation and particles connectivity can be accurately assessed by using the non-destructive X-ray tomography technique. The qualitative and quantitative analysis of tomography experimental results allows a 3D visualisation of the composite. The brittleness of the composites is caused by agglomerates and/or the presence of a percolated network: it explains why impact resistance decreases when the wood flour content increase. To optimise the rigidity/ductility compromise, several solutions are considered. The other part of this work consists of combining acoustic emission and X-ray tomography to study the evolution of damaging in composite materials with different wood flour content. An original coupling of these two techniques is implemented. The main failure mechanism involved in such composites is the decohesion between the matrix and the wood particles. The ductile/brittle transition occurs when the cavities coalescence becomes the dominant process over the PVC matrix failure. Adding impact modifiers is not sufficient to improve the ductility of the composite. Coupling agents are the route to be explored so as to improve the interfacial adhesion between the matrix and the particles.

Notre étude a porté sur la compréhension microscopique des mécanismes de déformation du polyamide 6 (PA6) à l’état massif. Par des traitements thermiques appropriés, on a obtenu un jeu d’échantillons présentant des microstructures semi-cristallines variées, avec différentes formes cristallographiques (allotropes : α, γ ou β), différents taux de cristallinité (de 24 à 35%), différentes périodes de l’empilement des lamelles cristallines (de 7 à 12nm), ceci pour deux masses moléculaires différentes. Les propriétés mécaniques en traction ont été caractérisées au-dessus et au-dessous de la transition vitreuse de la phase amorphe (Tg). Les différents matériaux présentent des différences notables de comportements. Le taux de cristallinité semble être le facteur prédominant au-dessus de Tg, mais d’autres facteurs sont à prendre en compte en dessous de Tg, comme la forme cristalline et la morphologie lamellaire (épaisseur et facteur de forme). Grâce à un dispositif expérimental fabriqué sur mesure, des essais de traction ont été suivis par diffusion des rayons X aux petits (SAXS) et grands angles (WAXS) sur la ligne D2AM, ESRF, pour caractériser les déformations à l’échelle des empilements lamellaires et à l’échelle de la maille cristalline. Dans les échantillons présentant principalement de la phase cristalline β, les lamelles tendent à s’orienter perpendiculairement à la direction de traction (TD). Ce mécanisme d’orientation local (que nous appelons « modèle de réseau de chaînes ») est induit par la transmission des contraintes par les chaînes amorphes reliant les lamelles cristallines adjacentes. L’allongement local est plus faible que l’allongement macroscopique dans les lamelles perpendiculaire à TD, ce qui implique que les lamelles inclinées doivent être cisaillées. De plus, la phase β se transforme en phase α aux fortes extensions. Dans les échantillons présentant principalement de la phase α (la plus rigide), au-dessus de Tg, dans le régime élastique, les chaînes tendent d’abord à s’orienter perpendiculairement à TD, ce qui implique que les lamelles s’orientent parallèlement à TD (« modèle de lamelles rigides »). Ensuite, dans le régime plastique, une majeure partie des lamelles se réoriente perpendiculairement à TD, comme dans le « modèle de réseau de chaînes », tandis qu’une fraction mineure reste orientée parallèlement à TD. Une morphologie fibrillaire fortement orientée est finalement obtenue pour tous les échantillons quelle que soit la température. / Mechanical properties of bulk polyamide 6 (PA6) have been studied in relation to microscopic deformation mechanisms. By applying various thermal treatments, sets of samples with different semi-crystalline microstructures, namely various crystalline allotropic forms (denoted α, γ and β) and different values of the crystallinity index (from 24 to 35%) and of the long period of the lamellar stacks (from 7 to 12 nm), have been obtained, for two different molecular masses. Mechanical properties have been measured in the linear (viscoelastic) and nonlinear (plastic) regimes below and above the glass transition of the amorphous phase (Tg). Differences of behavior have been observed depending on the microstructure. While the crystallinity index seems to be the predominant factor affecting the mechanical behavior above Tg, other structural parameters such as the crystalline form and the lamellar morphology (thickness and aspect ratio) have to be taken into account below Tg. Deformations at the scales of lamellar stacks and of the crystalline unit cell have been characterized by small and wide angle X-ray scattering (SAXS and WAXS) performed in-situ during tensile tests. In samples with predominantly β phase, lamellae tend to orient perpendicular to the tensile direction (TD). This orientation mechanism (which we denote as ‘Chain Network model’) is driven by the amorphous chains which transmit the stress between adjacent lamellae. The tensile strain in lamellar stacks perpendicular to TD is lower than the macroscopic tensile strain, which must be compensated by increased shear in inclined stacks. Also, at high extension ratios, the β phase transforms into α phase. In samples with predominantly α phase and above Tg, morphology changes are more complex. In a first step, chains orient perpendicular to TD, which implies that lamellar planes tend to orient parallel to TD, possibly due to their high aspect ratio (denoted as ‘Rigid Lamella’ model). In a second step, beyond the yield, a major fraction of crystallites then reorients normal to TD, i.e. chains themselves become parallel to TD, while a minor fraction remains oriented along TD. A highly oriented fibrillar morphology is ultimately obtained in all cases.

Dendrimers are a class of monodisperse polymeric macromolecules with a well defined and highly branched three-dimensional architecture. Their well-defined structure and structural precision makes them outstanding candidates for the development of new types of multifunctional super-molecules and materials with applications in medicine and pharmacy, catalysis, electronics, optoelectronics, etc. Liquid Crystalline Dendrimers (LCDrs) are a relatively new class of super-molecules which are based on the functionalization of common dendrimers with mesogenic (liquid crystalline) units. The combination of the fascinating molecular properties of the common dendrimers with the directionality of the mesogenic units have produced a novel class of liquid crystal forming super-mesogens (LCDRs) with unique molecular properties that allow novel ways of supramolecular self-assembly and self-organisation. This work is mainly concerned with the computational modelling of LCRs. A coarse grain strategy is adopted for the development of computational tractable models which take explicitly into account the specific architecture, the extended flexibility and the shape anisotropy of the mesogenic units of LCDRs. The developed force field applies easily to a variety of dendritic architectures. Utilizing Monte Carlo computer simulations we study the structural and conformational behavior of single LCDrs and of systems of LCDrs either in confined geometries or in the bulk. Special emphasis is given on the modeling of the response of LCDRS on externally applied alignment fields. External fields might be fictitious aligning potentials which mimic electric or magnetic fields or fields induced by the confining substrates. The surface alignment of liquid crystalline dendrimers (LCDrs) is a key factor for many of their potential applications. We present results from Monte Carlo simulations of LCDrs adsorbed on flat, impenetrable aligning substrates. A tractable coarse-grained force field for the inter-dendritic and the dendrimer-substrate interactions is introduced. The developed force field is based on modifications of well-known interaction potentials that can be used either with MC or with molecular dynamics simulations. We investigate the conformational and ordering properties of single, end-functionalized LCDrs under homeotropic, random (or degenerate) planar and nidirectional planar aligning substrates. Depending on the anchoring constrains to the mesogenic units of the LCDr and on temperature, a variety of stable ordered LCDr states, differing in their topology, are observed and analyzed. The influence of the dendritic generation and core functionality on the surface-induced ordering of the LCDrs are examined. The study has been extended to system of LCDrs confined in nano-pores of different shapes and sizes under several anchoring conditions. Two basic confining geometries (pores) considered in this work: slit and cylindrical pores. In each confining geometry, different anchoring conditions are imposed. The Isobaric-Isothermal (NPT) Monte Carlo Simulation is used to investigate the thermodynamic and structural properties of these nano-confined systems. The ransmission of orientational and positional ordering from the surface to the middle region of the pore depends on the size of the pore as well as on temperature and on anchoring strength. In the case of cylindrical pore, alignment propagation is short ranged compared to that of slit-pore. As a benchmark of our coarse-grained modelling strategy, we have extended and tested our coarse grained Force Field for the study of Janus-like dendrimers confined on planar substrates. The obtained results indicate the capability of our model to capture successfully the highly amphipilic nature of these class dendrimers and their self-organisation properties. / Τα δενδριμερή είναι μία κατηγορία μονοδιάσπαρτων πολυμερικών μακρομορίων με δενδρόμορφη τρισδιάστατη αρχιτεκτονική. Η μοριακή αρχιτεκτονική τους και η μονοδιασπορά τους καθιστούν τα δενδριμερή ιδανικά ως πολυλειτουργικά υπερ-μόρια με εφαρμογές στην ιατρική και τη φαρμακολογία, την κατάλυση, την ηλεκτρονική και οπτοηλεκτρονική κλπ. Τα Υγρό-Κρυσταλλικά Δενδριμερή (ΥΚΔ) είναι μια σχετικά νέα κατηγορία υπερ-μορίων που βασίζονται στη χημική τροποποίηση των κοινών δενδριμερών με μεσογόνες (υγρόκρυσταλλικές) μοριακές μονάδες. Ο συνδυασμός των ιδιαίτερων μοριακών ιδιοτήτων των κοινών δενδριμερών με την κατευθυντικότητα των μεσογόνων έχουν οδηγήσει σε μια νέα κατηγορία υπερ-μεσογόνων (ΥΚΔ) με μοναδικές μοριακές ιδιότητες που επιτρέπουν νέους τρόπους (υπερ)μοριακής αυτο-συναρμολόγησης και αυτο-οργάνωσης. Η εργασία αυτή ασχολείται με τη μοντελοποίηση και την υπολογιστική προσομοίωση ΥΚΔ. Εισάγονται οι αρχές για μια αδροποιημένη μοντελοποίηση ΥΚΔ που να λαμβάνει ρητά υπόψη την ειδική αρχιτεκτονική, την εκτεταμένη μοριακή ευκαμψία και την ανισοτροπία σχήματος των μεσογόνων του ΥΚΔ. Τα δυναμικά αλληλεπιδράσεων που εισάγονται επιτρέπουν τη μοντελοποίηση ποικίλων δενδριτικών αρχιτεκτονικών. Με την χρήση υπολογιστικών προσομοιώσεων Monte Carlo μελετώνται οι μοριακές ιδιότητες απλών ΥΚΔ διαφόρων γενεών και αρχιτεκτονικών καθώς και η θερμοδυναμική συμπεριφορά και οι μετατροπές φάσεων συστημάτων ΥΚΔ. Ιδιαίτερη έμφαση δίνεται στην μοντελοποίηση της απόκρισης των LCDRS σε εξωτερικά πεδία που μπορούν να προκαλέσουν ευθυγράμμισης των μεσογόνων ομάδων του ΥΚΔ. Τα εξωτερικά εφαρμοζόμενα πεδία μπορεί να είναι δυναμικά ευθυγράμμισης που μιμούνται τα ηλεκτρικά ή μαγνητικά πεδία ή πεδία που επάγονται από τους γεωμετρικούς περιορισμούς (συνοριακές συνθήκες) που επιβάλλονται στο υλικό όταν βρίσκεται κοντά σε επιφάνειες ή περιορισμένο εντός πόρων. Η δυνατότητα επίτευξης κοινού μοριακού προσανατολισμού στις μεσοφάσεις από ΥΚΔ αποτελεί βασικό παράγοντα για πολλές από τις πιθανές εφαρμογές τους. Παρουσιάζονται αποτελέσματα προσομοιώσεων Monte Carlo Μόντε ΥΚΔ σε επαφή με επίπεδο, αδιαπέραστο υπόστρωμα που έχει τη δυνατότητα προσρόφησης (αγκύρωσης) των μεσογόνων μονάδων του ΥΚΔ υπό επιθυμητό προσανατολισμό. Τα αποτελέσματα βασίζονται σε κατάλληλα αδροποιημένο πεδίο δυνάμεων για την περιγραφή των αλληλεπιδράσεων μεταξύ των δενδριμερών καθώς και του δενδριμερούς με το υπόστρωμα. Ανάλογα με τον τύπο μοριακής αγκύρωσης στην επιφάνεια και τη θερμοκρασία, μια ποικιλία από διαφορετικούς τρόπους οργάνωσης του ΥΚΔ στην επιφάνεια παρατηρούνται και αναλύονται. Η μελέτη έχει επεκταθεί επίσης σε συστήματα ΥΚΔ περιορισμένα σε νανο-πόρους διαφόρων σχήματα και μεγεθών κάτω από διάφορες συνθήκες μοριακής αγκύρωσης. Οι δύο βασικές γεωμετρίες περιορισμού (πόροι) που μελετούνται αναφέρονται σε παραλληλεπίπεδους και κυλινδρικούς πόρους. Σε κάθε γεωμετρία επιβάλλονται διαφορετικές συνθήκες αγκύρωσης. Οι προσομοιώσεις Monte Carlo έγιναν στην ισοβαρή συλλογή (ΝΡΤ) και διερευνήθηκε η θερμοδυναμική συμπεριφορά καθώς και η μοριακή οργάνωση των συστημάτων υπό νανο-εγκλεισμό. Τα συστήματα αυτά παρουσιάζουν πλούσια θερμοδυναμική συμπεριφορά. Η μοριακή οργάνωση καθώς και η μετάδοση του προσανατολισμού και της τάξης θέσεων από την επιφάνεια προς την μεσαία περιοχή των πόρων εξαρτάται το σχήμα και τι μέγεθος του πόρου, από τη θερμοκρασία καθώς και από τις συνθήκες μοριακής αγκύρωσης στις επιφάνειες του πόρου. Για έλεγχο της αποτελεσματικότητας της στρατηγικής αδροποιημένης μοντελοποίησης που αναπτύχθηκε μελετήθηκαν επίσης αμφίφυλα δενδριμερή τύπου Janus περιορισμένα σε επίπεδη επιφάνεια. Τα αποτελέσματα των προσομοιώσεων έδειξαν την ικανότητα του μοντέλου μας να αποτυπώσει με επιτυχία το την αμφίφυλη φύση αυτών των δενδρομερών και να περιγράψει με επιτυχία διαφορετικούς τύπους αυτοργάνωσης που σχετίζονται με την αμφιφυλικότητα αυτών των μορίων και τον συνεπαγόμενο νανο-φασικό διαχωρισμό τους.

Les travaux de la présente thèse portent sur le recalage de modèles dynamiques d’aubes de turbinecomposites à matrice céramique. Ils s’inscrivent dans le cadre de la quantification d’incertitudes pour la validation de modèles et ont pour objectif de fournir des outils d’aide à la décision pour les ingénieurs desbureaux d’études. En effet, la dispersion importante observée lors des campagnes expérimentales invalidel’utilisation des méthodes de recalage déterministe. Après un état de l’art sur la relation entre les incertitudeset la physique, l’approche de Vérification & Validation a été introduite comme approche permettantd’assurer la crédibilité des modèles numériques. Puis, deux méthodes de recalages stochastiques, permettantde déterminer la distribution statistique des paramètres, ont été comparées sur un cas académique. La priseen compte des incertitudes n’élude pas les potentielles compensations entre paramètres. Par conséquent, desindicateurs ont été développés afin de détecter la présence de ces phénomènes perturbateurs. Ensuite, lathéorie info-gap a été employée en tant que moyen de modéliser ces méconnaissances. Une méthode derecalage stochastique robuste a ainsi été proposée, assurant un compromis entre la fidélité du modèle auxessais et la robustesse aux méconnaissances. Ces outils ont par la suite été appliqués sur un modèle éléments / This work is focused on the stochastic updating of ceramic matrix composite turbine blade model. They arepart of the uncertainty quantification framework for model validation. The aim is to enhance the existing toolused by the industrial decision makers. Indeed, consequent dispersion was measured during the experimentalcampaigns preventing the use of deterministic approaches. The first part of this thesis is dedicated to therelationship between mechanical science and uncertainty. Thus, Verification and Validation was introduced asthe processes by which credibility in numerical models is established. Then two stochastic updatingtechniques, able to handle statistic distribution, were compared through an academic example. Nevertheless,taking into account uncertainties doesn’t remove potential compensating effects between parameters.Therefore, criteria were developed in order to detect these disturbing phenomena. Info-gap theory wasemployed as a mean to model these lack of knowledge. Paired with the stochastic updating method, a robuststochasticapproach has been proposed. Results demonstrate a trade-off relationship between the model’sfidelity and robustness. The developed tools were applied on a ceramic matrix composite turbine blade finiteelement model.

Afin de pouvoir prédire le comportement des renforts de composites 3D interlock au cours d'un procédé de mise en forme, il est nécessaire de connaitre la position des mèches dans le renfort durant la phase de préformage du procédé. Les travaux présentés ici traitent de la simulation du préformage de renforts 3D épais à l'aide d'un élément fini hexaédrique semi-discret spécifique. En utilisant le principe des travaux virtuels, on distingue le travail interne virtuel dû à la tension des mèches des autres travaux virtuels. La raideur due aux tensions de mèches, qui constitue la contribution principale de la rigidité du matériau, est prise en compte à l'aide de barres incluses dans les éléments. Les rigidités dues aux autres sollicitations, comme la compression transverse, les cisaillements ou les frottements inter-mèches, sont décrites par un matériau continu additionnel. La combinaison de ce modèle discret du premier ordre et d'un matériau continu hyperélastique anisotrope dit du second ordre, pour formuler le comportement du matériau va permettre la simulation du préformage des renforts tissés épais. Conjointement aux travaux sur la simulation, des travaux expérimentaux pour l'identification des paramètres matériau de la loi de comportement ont été définis et réalisés. Ces paramètres concernent les deux parties de la formulation du comportement. / In order to simulate 3D interlock composite reinforcement behavior during forming process, it is necessary to predict yarns positions in the fabric during the preforming stage of the process. The present work deals with thick 3D interlock fabric forming simulation using specific hexahedral semi-discrete finite elements. Using the virtual work principle, we distinguish the virtual internal work due to tensions in yarns from other internal virtual works. The stiffness relative to yarns tension which is the main part of the rigidity is described by bars within the elements. The other rigidities - like transverse compression, shears or friction between yarns - are depicted by a continuous additional material. A combination of this "first order" discrete model and a continuous orthotropic hyperelastic "second order" material formulation will enable us to simulate the interlock preforming process. Jointly to the simulation work, we also had to specify and perform experimental testing identification of material parameters. These parameters concern both parts of the model.

La simulation du procédé de fabrication de renforts fibreux secs est un enjeu majeur pour l’étude de l’élaboration de matériaux composites, dont l’utilisation dans les industries de pointe s’intensifie rapidement. Ainsi, l’influence du métier à tisser sur la qualité des renforts est primordiale dans la caractérisation de leurs propriétés mécaniques. Une campagne d’essais expérimentaux est tout d’abord réalisée, de manière à identifier les phénomènes physiques mis en jeu. Les différents modes de déformation de la mèche sont ainsi étudiés : élongation, compaction, cisaillement et distorsion. Est étudié également le comportement en flexion et en frottement, afin de mieux appréhender l’effet du procédé de tissage sur les mèches. Deux types de lois de comportement élastiques sont envisagés : une loi hypoélastique et une loi hyperélastique. Sont développées les propriétés de chacune d’entre elles, ainsi que les grandeurs caractéristiques nécessaires à leur implémentation dans le logiciel commercial ABAQUS/Standard. Les algorithmes de deux subroutines sont présentés, correspondant à l’une ou l’autre de ces lois. Le choix est fait de modéliser le comportement mécanique de la mèche à l’aide d’une loi hyperélastique isotrope transverse de type St-Venant, par l’intermédiaire de la subroutine ABAQUS/Standard UANISOHYPER_INV. Enfin, une identification des paramètres matériau à l’aide d’une méthode inverse est proposée. Sont comparés les résultats obtenus par simulation avec les résultats expérimentaux. La loi de comportement alors déterminée permet de mettre en place des simulations de procédé de tissage. / Simulating the manufacturing process of woven preforms is a major stack for understanding the development of composite materials, used in high performance industries. The effect of the weaving loom on the preforms is very important to caracterize their mechanicals properties. Experimental tests are realised to identify the physical phenomenon. Different deformation modes are studied : elongation, compaction, shear and distortion. The bending and friction behavior are also important to understand the effect of weaving process. Two constitutive laws are considered : a hypoelastic law and a hyperelastic law. An analyse of their properties is presented, and their implementation in a commercial software, ABAQUS/Standard, is detailed. In this purpose, two subroutines can be used. The modelisation of the mechanical behavior of the tows is finally realised with a transversely isotropic hyperelastic St-Venant model, with the subroutine ABAQUS/Standard UANISOHYPER_INV. To conclude, an identification method is presented and the simulated results are compared to experimental tests. The obtained consitutive behavior is finally used to simulate the weaving process.

L'industrie aéronautique doit faire face aux nouvelles exigences environnementales, tout particulièrement concernant la réduction de la consommation des énergies fossiles. L'utilisation de matériaux composites plus léger permet de répondre en partie à cette attente. Pour limiter les coûts lors de la fabrication et du développement des composites à renforts tissés 3D, il est nécessaire d'utiliser des outils de simulation performants. En particulier, les outils existants, qui discrétisent à une échelle mésoscopique l'architecture des tissus 3D, ne tiennent pas compte de l'influence du procédé de fabrication sur la constitution de la structure textile. Si des outils numériques dédiés à la modélisation du procédé de tressage et de tricotage sont disponibles, il n'en est rien concernant le tissage. Cette étude avait donc pour but de s'intéresser plus particulièrement à la simulation du prodécé de tissage pour pouvoir obtenir une structure de tissu sèche déformée numériquement. La production de différentes architectures de tissu en verre E dans notre laboratoire nous a permis d'observer les différents éléments en contact avec le fil ou le tissu sur la machine à tisser, par le biais de l'utilisation d'une caméra rapide par exemple. Le développement d'un modèle numérique par éléments finis reproduisant le procédé de tissage a été réalisé. Une loi de comportement isotrope transverse fut utilisée pour modéliser les fils de verre. Des premières simulations numériques encourageantes pour la fabrication d'un tissu d'armure toile et d'un tissu d'armure croisé 2-2 sont présentées et comparées avec les tissus réels produits correspondants. / The aeronautical industry faces new challenges regarding the reduction of fossil fuel consumption. One way to address this issue is to use lighter composite materials. The ability to predict the geometry and the mechanical properties of the unit cell is necessary in order to develop 3D reinforcements in composite materials for these aeronautical applications. There is a difficulty to get realistic geometries for these unit cells due to the complexity of their architecture. Currently, existing tools which model 3D fabrics at a meso scale don't take into account manufacturing process influence on the shape modification of the textile structure. There is already some numerical tools that can model the braiding or knitting process, but none have been developed for weaving so far. Consequently, this study deals with the numerical simulation of the weaving process to obtain a deformed dry fabric structure. During the weaving process of E-glass fabrics, achieved in our laboratory, it has been observed that large deformations led to the modification of transverse section of meshes, or local density changes, that can modify the fabrics mechanical resistance. For this reason, a numerical tool of the weaving process, based on finite element modelling, has been developped to predict these major deformations and their influences on the final textile structure. The correlation between numerical results and fabrics produced with glass fibres has been achieved for plain weave and 2-2 twill.

Le travail présenté dans ce mémoire avait pour objectif d’étudier le processus de la pénétration d'eau dans les composites en carbone/époxyde dans un premier temps, et dans un deuxième temps, d’étudier l’effet de la prise en eau par ces matériaux sur les performances mécaniques des composites et leur joints collés. L'intégration de ces phénomènes physiques dans la modélisation numérique est d'une grande importance dans la prédiction de la durabilité d’une structure en composite subissant un vieillissement hygrothermique. Par conséquent, ce travail consiste non seulement en des observations expérimentales, mais aussi en des simulations numériques. Des corrélations entre les résultats obtenus permettent d’une part de mieux comprendre ce qui se passe dans un système composite avec l’assemblage collé soumis à des charges mécaniques, de l’initiation d’endommagement jusqu’à la rupture finale ; d'autre part, de valider un modèle numérique robuste dans le but de la conception et de l’optimisation. Les originalités de ce travail se situent à différents niveaux en proposant : 1. un nouveau modèle de diffusion à deux-phases permettant de mieux décrire l’effet de l’épaisseur des stratifiés sur la pénétration de l’eau; 2. un nouveau modèle RPM « Residual Property Model » afin de prévoir la dégradation des propriétés mécaniques due à la prise en eau ; 3. une nouvelle loi de traction-séparation linéaire-exponentiel pour décrire la courbe-R observée dans les essais DCB en mode I pur sur les composites stratifiés afin de les intégrer plus facilement dans les modèles numériques / Le travail présenté dans ce mémoire avait pour objectif d’étudier le processus de la pénétration d'eau dans les composites en carbone/époxyde dans un premier temps, et dans un deuxième temps, d’étudier l’effet de la prise en eau par ces matériaux sur les performances mécaniques des composites et leur joints collés. L'intégration de ces phénomènes physiques dans la modélisation numérique est d'une grande importance dans la prédiction de la durabilité d’une structure en composite subissant un vieillissement hygrothermique. Par conséquent, ce travail consiste non seulement en des observations expérimentales, mais aussi en des simulations numériques. Des corrélations entre les résultats obtenus permettent d’une part de mieux comprendre ce qui se passe dans un système composite avec l’assemblage collé soumis à des charges mécaniques, de l’initiation d’endommagement jusqu’à la rupture finale ; d'autre part, de valider un modèle numérique robuste dans le but de la conception et de l’optimisation. Les originalités de ce travail se situent à différents niveaux en proposant : 1. un nouveau modèle de diffusion à deux-phases permettant de mieux décrire l’effet de l’épaisseur des stratifiés sur la pénétration de l’eau; 2. un nouveau modèle RPM « Residual Property Model » afin de prévoir la dégradation des propriétés mécaniques due à la prise en eau ; 3. une nouvelle loi de traction-séparation linéaire-exponentiel pour décrire la courbe-R observée dans les essais DCB en mode I pur sur les composites stratifiés afin de les intégrer plus facilement dans les modèles numériques