Lyapunov stability results are given for differential/algebraic models of power systems which include the effect of generator damping and nonlinear loads. The global dynamical structure of such a system is studied in terms of multivalued energy functions defined on so-called ‘voltage causal regions’ where voltage behaviour is predicted from angle behaviour. These regions are separated by ‘impasse surfaces’ related to singularity in the load flow equations.
