All cross sections, energy spectra, and double-differential cross sections of six outgoing particles (neutron, proton, deuteron, triton, 3He, and alpha particle) for the p+90Zr reaction are consistently calculated and analyzed at incident proton energies below 200 MeV by using nuclear theoretical models. The optical model; distorted wave Born approximation theory; the intranuclear cascade model; the exciton model; Hauser-Feshbach theory and the evaporation model; and, especially, the improved Iwamoto-Harada model of light composite particle (deuteron, triton, 3He, and alpha-particle) emission have been included in the calculation, which is in good agreement with the experimental data.