It has been suggested that the frictional interaction of neutron superfluids with normal matter in the inner crust of neutron stars dissipates rotational energy of superfluids and generates heat. Incorporating a general formula of internal heating into the detailed numerical codes of thermal evolution, we examine the effects of the internal heating on thermal evolution of neutron stars. We find that when a very stiff equation of state is used, it takes as long as approximately 2 x 10(4) yr for the interior of a neutron star to reach the isothermal state, even if a strong heat source is placed in a thin layer of the inner crust. This time scale reduces to a few hundred years or less when medium to soft equations of state are used. A neutron star cools by neutrino emissions during the earlier stages referred to as the neutrino cooling era, while during the later photon cooling era it cools primarily by emission of photons from its surface. We show that heating rates expected in the current superfluid-crust interaction model can greatly increase the surface temperature in the photon cooling era, significantly changing the thermal evolution of relatively old neutron stars if a stiff equation of state is adopted. This effect is less important for stars with medium to soft equations of state. Even in the presence of an ''exotic'' fast cooling mechanism, the frictional heating may significantly increase the internal temperature already during the neutrino cooling era. We compare the results of our calculations with the Einstein, EXOSAT, and ROSAT results for the thermal emissions from radio pulsars, and discuss the constraints imposed on the thermal evolution models.