On the development of a harmonic balance method for aeroelastic analysis