Dans ce travail, nous mettons en place une méthode fournissant des bornes à la fois garanties et précises de l'erreur de discrétisation commise sur des quantités locales lors de simulations numériques menées par la Méthode des Éléments Finis. On se focalise sur les problèmes de viscoélasticité linéaire décrits par variables internes. A partir d'une méthode générale développée au LMT Cachan, plusieurs optimisations sont abordées parmi lesquelles la prise en compte des effets d'histoire, le caractère non intrusif de la méthode ou l'estimation de l'erreur de modèle. La pertinence de la méthode est démontrée sur diverses applications numériques en 1D, 2D et 3D.