Le but de ce travail est de développer des outils par éléments finis permettant de calculer des structures sandwich sans ajouter aucune autre hypothèse que celle, très générale, de l'élasticité linéaire. Deux voies ont été suivies dans ce travail. Tout d'abord, nous exposons comment développer des éléments finis hybrides, dans lesquels les contraintes aux interfaces entre les différentes couches sont obtenues par l'intermédiaire de multiplicateurs de Lagrange. Des éléments 2D et 3D, basés sur le principe du travail virtuel et sur la fonctionnelle de Pian et Tong sont développés et validés. Ensuite, après avoir examiné les possibilités offertes par le code de calcul Ansys 5.3, nous proposons une méthode de post-traitement, basée sur la fonctionnelle de Reissner, permettant de recouvrer les contraintes d'interfaces à partir des déplacements nodaux obtenus à partir des éléments en déplacements. Finalement ces deux approches sont confrontées entre elles, tant sur le plan de la qualité des résultats numériques que sur la facilité d'utilisation, de programmation et d'implémentation dans des codes de calculs existants.