Grain-oriented electrical steels (GOES) are widely used in the manufacturing of high efficiency energy conversion systems thanks to their excellent magnetic properties along the rolling direction. Nevertheless, on the one hand, GOES exhibit strong magnetic anisotropy, resulting in distinct magnetic properties regarding the direction of the applied magnetic field to the rolling direction. On the other hand, GOES undergoes changes in their magnetic properties due to the industrial manufacturing process and mechanical constraints during operation. In general, stress, and particularly compressive stress, deteriorates magnetic properties. This paper deals with the modeling of the magnetic behavior of GOES under compressive stress. The orientation distribution function (ODF) based approach recently applied to describe the B-H first magnetization curves of a conventional GOES without stress as well as under applied uniaxial mechanical tensile stress is extended to account for the effect of compressive stress, with different considerations of modeling. The magnetic properties and the sensitivity of the ODF-based model to the magnetization direction as well as to the applied compressive stress are analyzed. Oscillation issues inherent to the ODF-based model are also discussed.