$$ \frac{d^2l}{dt^2} + c \frac{dl}{dt} +k= F_0\sqrt{1-l^2}$$

Either a closed form or something like given initial conditions, the max value of $l$. For small $l$, this boils down to simple step response. But for larger variations it doesn't work.