Fusion cross-sections are computed for the $^{40}$Ca$+^{40}$Ca system over a wide energy range with two microscopic approaches where the only phenomenological input is the Skyrme energy density functional. The first method is based on the coupled-channels formalism, using the bare nucleus-nucleus potential calculated with the frozen Hartree-Fock technique and the deformation parameters of vibrational states computed with the time-dependent Hartree-Fock (TDHF) approach. The second method is based on the density-constrained TDHF method to generate nucleus-nucleus potentials from TDHF evolution. Both approaches incorporate the effect of couplings to internal degrees of freedoms in different ways. The predictions are in relatively good agreement with experimental data.