Quantitative analysis of iron loss taking magnetic hysteresis behavior into account is essential to development of high-efficiency electric machines. In a previous paper, a novel magnetic circuit model incorporating a play model, which is one of the phenomenological models of magnetic hysteresis, was proposed. It was clear that the proposed model can calculate the hysteresis loop of the magnetic reactor with high-speed and high-accuracy. However, this method was applied only for the objects with simple shapes such as a ring core. Thus, it is necessary to extend the applicable range for devices with more complicated shapes such as electric motors. This paper describes that the play model is applied to reluctance network analysis (RNA) by using a permanent magnet (PM) motor as an object of discussion, in order to estimate the iron loss including magnetic hysteresis behavior.