HYDJET++ is a Monte Carlo event generator merging parametrized soft part inspired by hydrodynamics with hard part containing jets. It has been successful to describe particle production in Pb+Pb collisions at $\sqrt{s_{NN}}=2.76$ TeV energies. In this contribution, particle spectra and collective flow for the top LHC energy $\sqrt{s_{NN}}=5.02$ TeV Pb+Pb collisions are presented for the first time. Specifically, the HYDJET++ model version 2.4 has been used to simulate spectra of charged particles, $D^0$ and $J/\psi$ mesons and related $v_2$ and $v_3$ azimuthal flow harmonics. The particle spectra and flow haromincs are studied in different centrality bins ranging from 0-10% up to 30-40% centrality in midrapidity region for charged particles and $D^0$ mesons and in forward rapidity in case of $J/\psi$ mesons. The simulated results have been compared with the LHC data to tune HYDJET++ parameters.