Thèse ou mémoire avec insertion d'articles. / L'un des défis majeurs de la simulation numérique de la résistance au crash des structures composites aéronautiques est de pouvoir prédire les endommagements, leur initiation, leur évolution au cours de l'écrasement, et l'énergie absorbée, à partir d'un nombre limité de propriétés matériaux. Le but de la thèse est d'améliorer la compréhension des mécanismes élémentaires impliqués dans l'écrasement des composites à renforts tissés taffetas (plain weave PW) fabriqués en ou hors autoclave à base d'Époxyde et de fibres de carbone et de développer un modèle numérique performant. Ce rapport présente une investigation numérique et expérimentale effectuée au sein du département de Génie Mécanique à l'Université Laval dans le cadre du projet CRIAQ Comp-410 "Impact Modeling of Composite Aircraft Structures" en collaboration avec l'Université de Waterloo, Bombardier Aerospace, Bell Helicopter, le conseil national de recherche du Canada et DRDC-Valcartier. L'objectif principal est de développer des méthodes de simulations numériques prédictives pour différentes vitesses de crash des composites tissés. La démarche consiste en deux grandes étapes, une étude du crash en quasi statique et une autre étude dynamique pour voir l'effet du taux d'endommagement. Des simulations numériques avec des modèles de comportements existants dans les codes commerciaux éléments finis ont été effectuées afin de déceler les avantages et les inconvénients de chaque modèle, ensuite, un nouveau modèle 3D est développé pour tenir compte principalement, du délaminage, de la fragmentation, des bandes de pliage, de l'inélasticité et de l'endommagement. En parallèle, un plan expérimental est mis en œuvre pour la validation des différents modèles numériques et pour faire une étude paramétrique des différents paramètres influents, tels que, la forme des sections des éprouvettes, les initiateurs de crash, l'effet de l'échelle des éprouvettes et les séquences d'empilement. La mise en œuvre d'un modèle mathématique performant pour prédire le comportement mécanique des structures composites soumises au crash est une lourde tâche avec de multiples paramètres de matériaux et mécanismes de déformation à prendre en considération. Cette thèse peut être considérée comme une nouvelle contribution à l'avancement de la compréhension du processus d'endommagement lors du crash des matériaux composites. Le nouveau code numérique éléments finis développé a été validés expérimentalement et une étude de sensibilité des paramètres a été effectuée pour mesurer l'effet et le degré d'influence de chaque paramètre sur la précision de la solution finale. / One of the major challenges of the numerical simulation of the crashworthiness of aeronautical composite structures is to be able to predict damages, their initiation, and evolution during crash, and energy absorbed, from limited material properties. The goal is to improve the understanding of the elementary mechanisms involved during crush of plain weave (PW) fabric CFRP composites made in or out of autoclave to develop a powerful numerical model. This report presents a numerical and experimental investigation carried out within the Department of Mechanical Engineering at Université Laval as part of the CRIAQ Comp-410 project "Impact Modeling of Composite Aircraft Structures" in collaboration with the University of Waterloo, Bombardier Aerospace, Bell Helicopter, National Research Council of Canada and DRDC-Valcartier. The main objective is to develop predictive numerical simulation methods for different crash scenario of woven composites. The approach consists of two main steps, a quasi-static crash study and another dynamic study to see the effect of the strain rate. Numerical simulations with existing material models implemented in finite elements commercial codes were carried out in order to detect the advantages and the disadvantages of each model, then, a new 3D material model is developed to take account, mainly, of delamination, fragmentation, kink-bands, inelasticity and damage. In parallel, an experimental plan is implemented for the validation of the different numerical models and to make a parametric study of the different parameters, such as, the section shapes of the specimens, the triggers, the effect of scale and layup of the specimens. The implementation of a powerful mathematical model to predict the mechanical behavior of composite structures subjected to crash is a heavy task with multiple material parameters and deformation mechanisms to consider. This thesis can be considered as a new contribution to the advancement of the understanding of the process of damage during the crash of composite materials. The new developed material model developed has been validated experimentally and a parameter sensitivity study is performed to measure the effect and degree of influence of each parameter on the accuracy of the final solution.

Ce mémoire présente une étude de la délamination dans des matériaux composites tissés taffetas. Des essais expérimentaux ont été réalisés afin de caractériser la délamination selon les modes purs I et II et pour plusieurs combinaisons des modes mixés I et II, respectivement avec des essais DCB, ENF et MMB de la mécanique de la rupture. Les standards ASTM ont été suivis. Deux stratifiés ont été testés, un [0/90]12 et un [45/-45]12, et leurs résultats sont comparés. Des modèles numériques ont été créés sur Abaqus afin de simuler les essais expérimentaux. La méthode de la zone cohésive a été utilisée dans des analyses statiques en utilisant le solveur Abaqus/Standard. Les résultats des modèles numériques sont comparés aux résultats expérimentaux et leur validité est discutée. / This thesis is about the delamination of plain weave woven composite materials. Some experimental tests have been conducted to characterize the delamination for the pure modes I and II and for several combinations of mixed mode I and II, respectively with the DCB, ENF and MMB fracture mechanics tests. ASTM standards have been followed. Two laminate lay-ups were tested, one [0/90]12 and one [45/-45]12 and their results are compared. Numerical models have been created on Abaqus to simulate the experimental tests. The method of cohesive zone model has been used for static analyses using the Abaqus/Standard solver. Numerical models results are compared to experimental results and their validity is discussed.