  Empirical 'backward' galaxy evolution models for IR-bright galaxies are
constrained using multiband IR surveys. A new Monte-Carlo algorithm is
developed for this task. It exploits a large library of realistic Spectral
Energy Distributions (SEDs) of 837 local IR galaxies (IRAS 25$\mu m$ selected)
from the UV (1000{\AA}) to the radio (20cm), including ISO-measured 3--13$\mu
m$ unidentified broad features (UIBs). The basic assumption is that the local
correlation between SEDs and Mid-Infrared (MIR) luminosities can be applied to
earlier epochs of the Universe. Three populations of IR sources are considered
in the evolution models. These include (1) starburst galaxies, (2) normal
late-type galaxies, and (3) galaxies with AGN. A set of models so constructed
are compared with data from the literature. Predictions for number counts,
confusion limits, redshift distributions, and color-color diagrams are made for
multiband surveys using the upcoming SIRTF satellite.
