Simscape 

This example shows how to model a twopole variable reluctance actuator. Initially the rotor is held at 45 degrees by the rotational spring, and then at 0.1 seconds the motor is activated. The rotor is pulled towards zero degrees so as to minimize the magnetic circuit reluctance. The reluctance depends nonlinearly on rotor angle and current, saturating at higher flux values. When saturated, the currentflux relationship is defined by a quadratic function. Torque is then determined by differentiating flux with respect to rotor angle and integrating over the current.
The approach used in this model can be used to incorporate flux as a function of current and rotor angle for any motor as calculated by a finite element magnetic package. Here the underlying equations are approximated by the quadratic function, and a small number of data points are used to calculate the coefficients. An alternative would be to use the 2D lookup table in the Simscape™ Physical Signals sublibrary. Two tables are required, one for the flux and a second derived table of the local derivatives of the flux as a function of rotor angle and current.