Probabilistic aggregation of harmonic loads