Optimization design and robust design are significant measures for improving the performance and reliability of electromagnetic devices (EMDs, specifically refer to relays, contactors in this paper). However, the implementation of the above-mentioned design requires substantial calculation; consequently, on the premise of guaranteeing precision, how to improve the calculation speed is a problem that needs to be solved. This paper proposes a new method for establishing an approximate model for the EMD. It builds a relationship between the input and output of the EMD with different coil voltages and air gaps, by using a user-defined interpolating function. The coefficient of the fitting function is determined based on a quantum particle swarm optimization (QPSO) method. The effectiveness of the method proposed in this paper is verified by the electromagnetic force calculation results of an electromagnetic relay with permanent magnet.