Frequency domain modeling of nonlinear end stop behavior in Tuned Mass Damper systems under single- and multi-harmonic excitations