Ce travail vise Ce travail vise à améliorer la performance du code éléments finis FORGE®. En effet, la modélisation du contact dans ce code s'avère insatisfaisante pour simuler les procédés de mise en forme des matériaux à faible zone de contact. On présente dans ce mémoire les développements qui ont été réalisés pour améliorer le modèle de contact existant. Dans la première partie, partant de ce modèle, on s'intéresse à un problème spécifique qui est le procédé de laminage circulaire. Un nouveau modèle de traitement du contact dédié à ce procédé est développé. Il repose sur un schéma de réactualisation temporelle en coordonnées cylindriques, une formulation implicite et une description analytique des obstacles. La condition de contact incrémentale est alors imposée exactement en fin d'incrément. On valide ces développements sur plusieurs cas de procédé de laminage circulaire où l'on montre l'apport du nouvel algorithme en terme de précision de gestion de contact et aussi en terme d'amélioration des résultats de simulation et des temps de calcul. Dans la deuxième partie, on a généralisé cette approche à l'ensemble des procédés à faible zone pour lesquels une représentation analytique des outils n'est pas possible. On a développé alors un modèle implicite général basé sur la réactualisation de la condition de contact au cours des itérations de Newton-Raphson. Cette méthode permet de tenir compte exactement de l'évolution de la position relative de la pièce et de l'outil. Ensuite, on généralise le lissage des outils à des pièces de géométrie plus complexes. On étudie alors l'influence de la discrétisation de l'outil sur la précision de la simulation des procédés de mise en forme et particulièrement pour des procédés à faible zone de contact. Une méthode d'interpolation quadratique proposée par Nagata est alors implémentée pour le lissage de la surface des outils discrétisés par des éléments plans. La nouvelle surface interpolée est déterminée d'une façon locale à partir seulement de la position des nœuds et des normales nodales. Ces normales sont calculées à partir des normales aux facettes entourant le nœud considéré en utilisant la méthode des normales votées. L'efficacité de cette méthode combinée à schéma d'intégration implicite est testée pour plusieurs procédés de mise en forme.