The enhancement of radiative strength function (RSF) in the region of low $gamma$-rays energy ($E_{gamma}leq 12$ MeV), which is caused by the pygmy dipole resonance (PDR), is treated within the phonon damping model (PDM) plus exact thermal pairing (EP) without adding any extra PDR strength function. The numerical calculations performed for $^{161-163}$Dy show that, because of the effect of EP, the EP+PDM can describe reasonably well the total RSF data in both low- and high-energy regions of $gamma$-rays. Consequently, as compared to the conventional description within the phenomenological generalized Lorentzian (GLO) and standard Lorentzian (SLO) models, the EP+PDM calculations can eliminate at least eight free parameters. This indicates the important role of microscopic approaches towards the precise description of the RSF. In particular, temperature is found to have significant contributions to the RSF below the neutron separation energy, questioning again the validity of the Brink-Axel hypothesis in this energy region.