Este artículo presenta un nuevo método para calcular estados de equilibrio magnetohidrodinámicos (MHD) 3D mediante redes neuronales artificiales. En comparación con los solucionadores convencionales, parametrizamos los modos de Fourier como redes neuronales artificiales y minimizamos la fuerza residual no lineal global en el espacio real mediante una técnica de optimización de primer orden. Como resultado, logramos el mismo nivel de mínimo residual que el calculado con códigos convencionales a un coste computacional competitivo y, a costa de un mayor coste computacional, obtenemos un mínimo residual menor mediante redes neuronales, lo que sugiere un nuevo límite inferior para la fuerza residual. Con una complejidad mínima de redes neuronales, esperamos mejoras significativas no solo en la resolución de estados de equilibrio individuales, sino también en el cálculo de modelos de redes neuronales válidos para distribuciones de equilibrio continuas.