Fusion excitation functions and angular distributions of evaporation residues (ERs) have been measured for 28 Si + 90,94 Zr systems around the Coulomb barrier using the recoil mass spectrometer, Heavy Ion Reaction Analyzer (HIRA). For both systems, the experimental fusion cross sections are strongly enhanced compared to the predictions of the one-dimensional barrier penetration model (1-d BPM) below the barrier. Coupled channels formalism has been employed to theoretically explain the observed sub-barrier fusion cross section enhancement. The enhancement could be explained by considering the coupling of the low-lying inelastic states of the projectile and target in the 28 Si + 90 Zr system. In the sub-barrier region, the measured fusion cross sections for 28 Si + 94 Zr turned out to be about an order of magnitude higher than the ones for the 28 Si + 90 Zr system, which could not be explained by coupling to inelastic states alone. This observation indicates the importance of multinucleon transfer reaction channels with positive Q values in the sub-barrier fusion cross section enhancement, because 90,94 Zr are believed to have similar collective strengths. This implies that no strong isotopic dependence of fusion cross sections is expected as far as the couplings to collective inelastic states are concerned. In addition, the role of projectile and multiphonon couplings in the enhancement has been explored.