Response of hysteretic multidegree-of-freedom systems using harmonic balance method