This paper presents an alternative solution for the power-flow analysis of power systems with distributed generation provided by heterogeneous sources. The proposed simulation approach relies on a suitable interpretation of the power network in terms of a nonlinear circuit in the phasor domain. The above circuit interpretation can be solved directly in the frequency-domain via the combination of a standard tool for circuit analysis with an iterative numerical scheme, providing directly the steady-state solution of the power-flow of a generic distribution network. At each iteration, the resulting circuit turns out to be composed by two decoupled subnetworks, a large linear part and a set of smaller nonlinear pieces accounting for the load characteristics, with evident benefits in terms of the computational time. The feasibility and strength of the proposed simulation scheme have been verified on a large benchmark consisting of the IEEE 8500-node test feeder. Then it is applied to the statistical simulation of a power network accounting for the variability effects of renewable generators. According to the results, the proposed tool provides an effective alternative to the state-of-the-art approaches for power-flow analysis further highlighting the benefits of the application of well-established tools for circuit analysis to power-flow problems.