This work develops and benchmarks surrogate models for Dynamic Phasor (DP) simulation of electrical drives. DP simulations of complex systems may be time-consuming due to the increased number of equations. Thus, it is desirable to have a datadriven approach to compute the critical state/control variables and power losses. The surrogate models are intended to be used as a steady-state equivalent of the DP simulation model. We consider the Gaussian Process (GP), Multi Layer Perceptron, and Random Forest as surrogate models. Among other techniques, GPs are found to have good accuracy. Moreover, GPs are dataefficient and have desirable properties, such as built-in uncertainty quantification. The study shows that the GP performs better compared to other techniques in terms of the Mean Squared Error of the prediction, while still being very fast to evaluate. We illustrate the potential of these surrogate models to also predict transient behavior.