An important alteration of the equivalent loads profile has been observed in the electrical energy distribution systems, for the last years. Such fact is due to the significant increment of the electronic processors of electric energy that, in general, behave as nonlinear loads, generating harmonic distortions in the currents and voltages along the electric network. The effects of these nonlinear loads, even if they are concentrated in specific sections of the network, are present along the branch circuits, affecting the behavior of the entire electric network. For the evaluation of this phenomenon it is necessary the analysis of the harmonic currents flow and the understanding of the causes and effects of the consequent voltage harmonic distortions. The usual tools for calculation the harmonic flow consider one-line equivalent networks, balanced and symmetrical systems. Therefore, they are not tools appropriate for analysis of the operation and the influence/interaction of mitigation elements. In this context, this work proposes the development of a computational tool for the analysis of the three-phase harmonic propagation using Norton modified models and considering the real nature of unbalanced electric systems operation. © 2011 IEEE.