The dynamics of the F+HD reaction has been studied by means of quasiclassical trajectory calculations on an ab initio potential energy surface (PES) at several collision energies. At the collision energy of 85.9 meV and for the DF+H isotopic channel of the reaction, there is a remarkable agreement between calculated and experimental results, in both the center of mass (c.m.) differential cross sections (DCS) and in the simulation of the laboratory (LAB) time of flight (TOF) and angular distributions (AD). The good agreement also extends to the lower collision energy of 58.6 meV for this channel of the reaction. In contrast, the simulation of the LAB angular distributions for the HF+D channel shows strong discrepancies between theory and experiment at both collision energies, which can be traced back to the absence of a forward peak in the calculated c.m. DCS for HF(v′=3). Simulations made from QCT calculations on other PES with important HF(v′=3) forward scattering contributions also fail to reproduce the overall AD. The theoretical findings and especially the roles of translational energy and initial rotational momentum on the dynamics of this reaction are discussed in terms of the topology of the PES through the analysis of individual trajectories.