[en] The measured neutron energy spectra produced from thick targets by light-mass heavy ion beams were analyzed with the phenomenological hybrid model of equilibrium and pre-equilibrium emissions. The neutron energy spectra are fitted to two (in the energy region of incident projectile energy lower than 100 MeV) and three (in the energy region of incident projectile energy higher than 100 MeV) Maxwellian-type components. The lowest energy component corresponds to the evaporation neutrons from a compound nucleus having a nuclear temperature independent to the neutron emission angle and the higher two components to the pre-equilibrium neutron emission having a nuclear temperature depending on the angle. The total neutron yield was evaluated by surveying several published papers and estimating the neutron yield below a few MeV by fitting some spectra measured above that energy to the Maxwellian distribution. In order to use for the neutron activation cross section measurement and the induced radioactivity study, a pseudo-monoenergetic neutron beam course in the energy range from 15 to 40 MeV was installed by using the Be(p, n) reaction. The peak neutron energy changes with the proton energy and the beam intensity was relatively high. (author)