Cet article présente une nouvelle méthode de calcul des états d'équilibre magnétohydrodynamiques 3D (MHD) à l'aide de réseaux de neurones artificiels. Comparés aux solveurs conventionnels, nous paramétrons les modes de Fourier comme des réseaux de neurones artificiels et minimisons le résidu de force non linéaire global dans l'espace réel grâce à une technique d'optimisation du premier ordre. Ainsi, nous obtenons le même niveau de minimum résiduel que celui calculé par les codes conventionnels à un coût de calcul compétitif. Au prix d'un coût de calcul plus élevé, nous obtenons un minimum résiduel plus faible via les réseaux de neurones, suggérant une nouvelle borne inférieure pour le résidu de force. Avec une complexité minimale des réseaux de neurones, nous nous attendons à des améliorations significatives non seulement dans la résolution des états d'équilibre uniques, mais aussi dans le calcul de modèles de réseaux de neurones valables pour les distributions d'équilibre continues.