“…The construction of analytical energy functions is the central problem in direct stability analysis of electric power systems [1], [5], [8], [9], [16], [17], [21], [22], [23], [24], [30]. Traditionally, direct methods have been developed for networkreduction models, where all the load are expressed in constant impedances and the entire network representation is reduced to the generator internal buses [22], [24], [30].…”