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é.