The paper concerns the semianalytical modeling of parasitic forces in electromagnetic devices. The Fourier analyses technique is extended to incorporate a position dependent relative permeability and conductivity into the solutions of electromagnetic quantities. The spatial distribution of these material properties are described using Fourier series. As a result, the parasitic reluctance and eddy current force can be accurately modeled. To model transient behavior of the device, multiple time harmonics have been incorporated in the solutions. To validate the developed method, it is applied to a coreless linear motor and compared to finite element results.