Les logiciels d'éléments finis tels que PLAXIS sont devenus des outils incontournables pour la simulation et la conception des ouvrages géotechniques. Si les lois de comportement classiques sont généralement suffisantes pour les applications courantes, elles atteignent leurs limites lorsqu'il s'agit de problématiques plus complexes nécessitant la prise en compte de l'anisotropie, des déformations différées et du comportement thermoplastique des géomatériaux. L'interface User Defined Soil Model (UDSM) de PLAXIS permet d'implémenter des lois de comportement avancées intégrant différents phénomènes et couplages multiphysiques. Dans ce travail, un modèle UDSM élasto-viscoplastique anisotrope intégrant les couplages thermo-hydro-mécaniques (THM) est développé. Le modèle inclut plusieurs mécanismes d'écrouissage dépendant des déformations plastiques ainsi qu'une fonction d'écrouissage thermique permettant de reproduire les phénomènes d'effondrement thermique observés dans les argiles. Le modèle est validé à partir de plusieurs cas tests représentatifs, puis appliqué à la modélisation numérique d'un système de stockage d'énergie thermique en aquifère à haute température (HT-ATES).