Hi,
unfortunately there are two "bugs" in the harmonic response analysis:
1- it is still in linear domain
2- any load is applied as a sine function, isochronically (though you can shift phases between one another).
At this point, I see no other solution (sigh...) than to run a full transient analysis. By the way, you won't need the static run any more, but every substep will take the same time to solve as the static run, so be careful if your model is huge...
Regards