Spede, along the lines of what electricpete was doing with the fast decoupled method, I think the numerical analysis you are looking for is this.
The best iterative method we have is newton-raphson. Becuase this involves matrices, I am having trouble writing in the basic equation (in here). You can search for Newton-Raphson on Google and get some good info on it. The fast-decoupled method is a derivation of newton-raphson that only works for transmission line systems where conductors are primarily reactive and the system can be approximated such that real power flow only responds to phase relationship and reactive power only responds to voltage magnitude. Since you are attempting to solve the internal network of the machine, this would not work well for your application. Additionally, NR has been adapted to bus-system solutions, while your application would probably be best suited to a more direct mathematial definition of it, so you might also talk to the numerical methods folks. If you treat the control system/freq. drive as a separate entity, then you can do a series of iterative solutions, where you manually adjust the input from the freq. drive as it would be..make sense? Its a little much for a post, so definitely check the web resources on it.