Application of Melnikov´s method for computing heteroclinic orbits in a classical SMIB power system model

The main concern of this paper is the application of the Melnikov method for computing heteroclinic orbits in a classical SMIB power system model and testing the one-parameter transversality condition required in the BCU method for transient stability analysis. First, the Melnikov method is used to give a first approximation for the minimum damping coefficient or, alternatively, for the maximum input mechanical power that should be associated with the system so that the transversality condition is satisfied. This first approximation is obtained without solving any differential equations. After that, this approximation is used as an initial condition for a monotonically convergent process which is very easy to implement and does not require knowledge of fault trajectory at all