Dans ce travail, une méthode générale est mise en place pour le calcul de bornes garanties et précises de l'erreur de discrétisation commise lors des simulations numériques. Elle est appliquée ici au contrôle de quantités locales dans les problèmes de Mécanique (viscoélasticité linéaire) résolus par la Méthode des Eléments Finis. Plusieurs optimisations sont également abordées, parmi lesquelles la prise en compte des effets d'histoire, le caractère non-intrusif de la technique employée, ou l'estimation de l'erreur de modèle. La pertinence de la méthode est démontrée sur diverses applications numériques.