This paper proposes an analytical subdomain model for predicting magnetic field distributions in a three-phase double-stator permanent magnet synchronous machine (DS-PMSM) during open-circuit and on-load conditions. The geometric structure of DS-PMSM is quite challenging since the stator cores are located in the outer and inner parts of the motor, while the rotor magnets are placed between these two stators. Parameters that influence the motor performance in DS-PMSM include stator outer radius, stator inner radius, magnet thickness, magnet arc, slot opening, outer and inner airgap thickness and the number of winding turns. The analytical subdomain model proposed in this paper, which can accurately predict the performances of DS-PMSM with less computational time, has an excellent advantage as a rapid design tool. The model is initially generated using the separation of variables technique in four subdomains, namely, outer airgap, outer magnet, inner magnet, and inner airgap, based on Laplace’s and Poisson’s equations in polar coordinates. The field solutions in each subdomain are derived by applying the appropriate boundary and interface conditions. Furthermore, finite element analysis (FEA) is used to validate the analytical results in fractional DS-PMSM with a different number of slots between outer and inner stators and a non-overlapping winding configuration. The electromagnetic performances that have been evaluated are the slotted airgap flux density, back-emf and output torque. The results demonstrate that the proposed analytical model is able to predict the magnetic field distributions accurately in DS-PMSM.