All cross sections of neutron-induced reactions, angular distributions, energy spectra and double-differential cross sections are consistently calculated and analyzed for 54,56,57,58,natFe at incident neutron energies below 200 MeV based on the nuclear theoretical models, which are the optical model, preequilibrium and equilibrium reaction theories, and the distorted wave Born approximation theory. Theoretical calculated results are compared with existing experimental data and the evaluated results in ENDF/B-VII and JENDL-3. The optical model potential parameters are obtained according to the experimental data of total and nonelastic-scattering cross sections, and elastic-scattering angular distributions.