Le sujet de cette thèse se situe à la frontière entre les mathématiques appliquées et la mécanique. Il s'agit d'étudier, sous un angle mathématique, des problèmes piézoélectriques, c'est à dire des problèmes couplant action mécanique et action électrique, en présence de contact et de frottement. On s'intéresse notamment à l'analyse variationnelle et numérique de ces problèmes. La thèse comporte quatre parties. La première partie contient l'ensemble des outils mathématiques, numériques et mécaniques nécessaires à une bonne compréhension du travail réalisé par la suite. La deuxième partie aborde deux problèmes statiques de contact frottant entre un corps électro-élastique et une fondation. Dans le cas des formulations primales des problèmes, nous prouvons l'existence, l'unicité, et la dépendance continue des solutions faibles, exprimées en termes de déplacements et de potentiel électrique. Une formulation duale équivalente au problème précédent, exprimée en termes de contraintes et de déplacement électrique, est étudiée pour laquelle des résultats d'existence et d'unicité sont établis. La troisième partie aborde deux modèles de contact piézoélectrique dans un cadre évolutif de type quasistatique. Pour chaque modèle, on présente un résultat d'existence. Dans la quatrième partie, le travail porte sur l'analyse numérique et les simulations par différences finies en temps (Euler implicite) et éléments finis en espace. Dans le cas statique, on traite un problème électro-élastique avec contact frottant de type compliance. Le problème discret est posé et les estimations a priori de l'erreur sont obtenues. Le problème est écrit sous forme d'un Lagrangien augmenté couplé à un algorithme de type Newton généralisé. S'ensuit une simulation numérique en dimension deux d'espace et des vérifications numériques de convergence. Des résultats similaires sont obtenus dans le cas d'un problème de contact quasistatique électro-viscoélastique.