I am not sure that harmonic analysis can handle different phases in different nodes as excitation. I always understood harmonic analysis as like to put the object on a shaker table and excite a set of constrained nodes with the same forcing function. But maybe I’m wrong.
Concerning the frequency, you provide a range, not one single frequency, this is the whole idea of harmonic analysis: you sweep the frequency!. If you refer to APDL, the command is HARFRQ. Or do you use Workbench?
Hope this helps,