The growing demand for electricity increases the complexity of electric power systems which, when combined with limitations in the expansion of transmission systems, results in the operation of electrical systems closer to their limits, making them vulnerable to voltage stability problems. In this context, there is a gap in the market for the development of computational tools that can represent the electrical systems more appropriately, thereby improving the conditions of analysis. The present study formulates three non-classical load flow representations: multiple swing buses, primary frequency control, and frequency dependent transmission network parameters. Once used in the load flow problem solving, these models are also extended to allow the calculation of voltage stability indices of load buses, controlled voltage buses and swing buses. Numerical tests with a 6-bus test system are presented to verify the applicability and adequacy of the proposed models.