An energetic model of statistic domain behavior is applied on the magnetization process of grain oriented FeSi‐steel sheets. The physical constants of the model are derived from anisotropic energy contributions, initial susceptibility, coercivity, and saturation magnetization. Calculated in the [100], [011], and [111] – directions, the hysteretic magnetization curves indicate good agreement with measurements from the demagnetized state up to saturation. The better correspondence between theory and experiments at weak fields is achieved by using realistic initial occupation probabilities for the domain structure. This assumption leads to complex mutual dependencies of the statistic domain’s volume fractions, which have not been considered in previous works. ©1995 American Institute of Physics.