Les récentes découvertes et avancées technologiques dans la compréhension des matériaux ainsi que l'essor des outils informatiques d'aide au calcul ont contribué à la prolifération de matériaux intelligents avec un champ d'applications très large. Cette thèse s'inscrit dans le contexte d'utilisation des actionneurs piézoélectriques plutôt qu'une vision purement matériau. Le but est d'enrichir les bibliothèques de modèles de ces types d'actionneurs afin de faciliter leur prise en compte dans les phases de conception des systèmes complexes les intégrant. Le cahier des charges est que ces modèles incluent le plus possible les non-linéarités tout en restant aisés d'utilisation. Pour atteindre ces objectifs, nous proposons de faire un pont entre le domaine des experts des matériaux et celui de l'ingénieur en suivant une méthodologie claire. Dans un premier temps nous passons en revue les approches existantes dans la littérature ainsi que les solutions offertes par certains logiciels commerciaux. Une analyse des équations constitutives de la piézoélectricité associées aux conditions de fonctionnement de l'actionneur nous permet d'en déduire un premier modèle analogique. Ce dernier est ensuite traduit en bond graph pour en déduire des modèles blocs-diagramme. En plus de cet effort de formalisation, ces premiers modèles se distinguent de ceux proposés par les logiciels commerciaux en prenant mieux en compte la dynamique propre à l'actionneur. Nous proposons deux types de modèles. L'un rend uniquement compte du premier mode de résonance alors que le second rend compte de deux modes de résonance. Ensuite nous proposons des modèles prenant en compte les non-linéarités : l'approche de Preisach pour la modélisation de l'hystérésis statique et l'approche de Voigt dans le cas dynamique. Ces deux approches sont ensuite fusionnées dans le but d'avoir un model plus complet.