Les sources de couplage thermomécanique dans les matériaux viscoélastiques sont multiples : thermo-élasticité, dissipation visqueuse, évolution des caractéristiques mécaniques avec la température. La simulation numérique de ces couplages en calcul des structures présente encore un certain nombre de défis, spécialement lorsque les effets de couplage sont très marqués (couplage fort). De nombreuses approches algorithmiques ont été proposées dans la littérature pour ce type de problème. Ces méthodes vont des approches monolithiques, traitant simultanément l'équilibre mécanique et l'équilibre thermique, aux approches étagées, traitant alternativement chacun des sous-problèmes mécanique et thermique. La difficulté est d'obtenir un bon compromis entre les aspects de précision, stabilité numérique et coût de calcul. Récemment, une approche variationnelle des problèmes couplés a été proposée par Yang et al. (2006), qui permet d'écrire les équations d'équilibre mécanique et thermique sous la forme d'un problème d'optimisation d'une fonctionnelle scalaire. Cette approche variationnelle présente notamment les avantages de conduire à une formulation numérique à structure symétrique, et de permettre l'utilisation d'algorithmes d'optimisation. Dans ce travail on utilise l'approche variationnelle pour résoudre le problème thermo-visco-élastique fortement couplé, puis on compare plusieurs schémas algorithmiques afin de trouver celui qui présente les meilleures performances.