Multi-time probability density functions of the dynamic non-Gaussian response of structures