A model of large load areas for harmonic studies in distribution networks