Dans ce travail, les auteurs considèrent le problème du calcul des contraintes dans un multicouche avec bords chargés et avec couches et interfaces élastoplastiques endommageables. Le modèle multiparticulaire M4-5N originalement développé pour le calcul des efforts d’interface dans des multicouches élastiques est modifié pour le calcul des phénomènes non linéaires du problème. La méthode LATIN est ensuite employée pour la résolution du problème non linéaire. Les auteurs montrent à la fin quelques exemples d’application pour lesquels le modèle peut aider à la prédiction du délaminage.