Abstract: Compared to traditional hysteresis models like Preisach and Jiles-Atherton (J-A), the micromagnetic Landau-Lifshitz-Gilbert (LLG) equation offers advantages such as clearer physical interpretation and higher simulation accuracy. However, its application to macroscopic static hysteresis simulation of electrical steel sheets has been limited due to computational intensity and memory constraints. This paper first proposes a geometric model discretization method for the LLG equation based on the Representative Volume Element (RVE) approach, tailored to the structural characteristics of electrical steel sheets. Subsequently, simplified expressions for the energy terms within the total Gibbs free energy are developed for the millimeter scale. A rapid parameter identification method for the relevant parameters is also introduced. This leads to the derivation of a simplified LLG equation specifically for macroscopic static hysteresis simulation in electrical steel sheets. Finally, the proposed simplified LLG equation is used to simulate static hysteresis loops of industry-grade grain-oriented (GO) and non-oriented (NO) silicon steel sheet samples under various operating conditions. The simulations exhibit strong agreement with measurements. Comparisons in simulation accuracy with the widely applied Preisach and J-A hysteresis models further validate the superiority of the proposed model and methodology.