This paper studies several decoupled penalty methods to overcome the saddle point system of the steady state 2D/3D incompressible magnetohydronamics (MHD). These approaches combine the Oseen iteration and two-level technique with strong uniqueness condition $$0<\frac{\sqrt{2}C_{0}^{2}\max \{1,\sqrt{2}S_{c}\}\Vert {\mathbf{F }}\Vert _{-1}}{(\min \{R_{e}^{-1},S_{c}C_{1}R_{m}^{-1}\})^2}\le 1-\left( \frac{\Vert \mathbf{F }\Vert |_{-1}}{\Vert |\mathbf{F }\Vert _{0}}\right) ^{\frac{1}{2}}<1$$ satisfied. For the convenience of implementation, we employ two different simple Lagrange finite element pairs $$P_{1}b-P_{1}-P_{1}b$$ and $$P_{1}-P_{0}-P_{1}$$ for velocity field, pressure and magnetic field, respectively. Rigorous analysis of the optimal error estimate and stability are provided. We present comprehensive numerical experiments, which indicate the effectiveness of the proposed methods for both two dimensional and three-dimensional problems.