In this talk, we will first present several structure-preserving finite element methods for some MHD models and then present and analyze some iterative methods for the resulting nonlinear algebraic systems. Energy stability is crucial for the well-posedness and convergence of the schemes.