Les composites stratifiés sont aujourd'hui fortement utilisés dans l'industrie. Ils permettent en effet une conception optimisée des structures et offrent ainsi une réponse aux contraintes énergétiques auxquelles font face les industriels. Cependant, la prévision par le calcul de leur comportement sous charge et de leur dégradation progressive soulève une difficulté : la modélisation du matériau doit être réalisée à une échelle très inférieure à celle de la structure. On se concentre ici sur les modèles de zone cohésive qui, depuis les vingt dernières années, ont démontré leur aptitude à la prédiction du délaminage. La taille des problèmes résultant de la simulation de l'évolution de structures industrielles par ce type de modèle est considérable. Leur résolution requiert alors l'utilisation de stratégies de calcul intensif. Dans ces travaux, le problème non linéaire posé à chaque piquet de temps du schéma d'intégration temporel est résolu par une méthode de décomposition de domaine à trois échelles. L'introduction des différents niveaux de résolution permet de traiter les informations de l'échelle fine, de l'échelle du pli et de l'échelle de la structure de manière adaptée au cours des itérations, et ainsi d'assurer l'efficacité de la méthode proposée, en termes de taux de convergence et d'aptitude à la parallélisation. Cette stratégie est ensuite aménagée pour résoudre les difficultés liées aux non-linéarités fortement localisées et aux éventuelles instabilités engendrées par le comportement adoucissant du matériau simulé.

Lors de la tolérance aux dommages en fatigue de structures aéronautiques, l'estimation des vitesses de propagation de fissures existantes permet d'évaluer la tenue résiduelle de la structure et de définir des intervalles d'inspection et de réparations éventuelles. Une analyse complète ne pouvant être menée sur un modèle éléments finis dit « boulon par boulon », des analyses détaillées sont réalisées dans les zones fissurées et incluses dans un modèle général grossier de l'avion. Dans ce travail, une stratégie de calcul multiéchelle répondant à la problématique d'analyse locale globale pour le suivi de fissure et permettant de lever les difficultés de maillage est proposée. Deux techniques sont mises en synergie : une approche multiéchelle basée sur une méthode de décomposition de domaine pour traiter efficacement la localisation des déformations dans la zone fissurée et une technique d'enrichissement local, la X-FEM, pour décrire la géométrie de la fissure indépendamment du maillage.