Modeling of periodic non-smooth motion using the harmonic balance method