Nous présentons dans ce travail un modèle numérique performant pour le calcul de la température, du flux thermique, des déplacements et des contraintes dans un multicouche épais. Ce modèle est basé sur une approche "couche discrète " et utilise des fonctions polynomiales de degré 2 dans l'épaisseur de chaque couche. La spécificité du modèle proposé est qu'il assure la continuité, tant pour tant pour la composante normale du flux thermique que pour le vecteur contrainte suivant la normale aux interfaces. Ce modèle autorise la présence de résistance de contact aux interfaces ; il permet de diminuer le nombre de paramètres du modèle tout en assurant une meilleure approximation. Pour valider ce modèle, nous avons été amené à expliciter les solutions analytiques exactes, pour certaines géométries. Nous avons également comparé ce modèle à des modèles couramment utilisés par les codes éléments finis spécifiques aux multicouches et nous avons montré qu'il apportait une amélioration certaine sur la précision des résultats, pour des problèmes de conduction thermique, d'élasticité et de thermoélasticité.