Observations are reported on isotactic polypropylene in uniaxial tensile relaxation tests on specimens subjected to tension up to various maximum strains and retraction down to various stresses. Noticeable evolution of shapes of relaxation curves under retraction is revealed with stress at the beginning of relaxation process: with a decrease in this stress, relaxation diagrams characterized by a monotonic decay of stress with time (simple relaxation) become, first, non-monotonic (mixed relaxation), and, finally, monotonically increasing (inverse relaxation). A thorough investigation is performed on the effect of multi-cycle preloading (maximum strain per cycle, minimum stress per cycle, number of cycles, and strain rate) on transition from simple to mixed and to inverse relaxation. It is found that (1) intensity of mixed relaxation increases with maximum strain when this parameter remains below the yield strain and decreases in the post-yield region of deformations, (2) an increase in number of cycles under preloading leads to reduction of intensity of mixed relaxation and transition from mixed to inverse relaxation, (3) inverse relaxation diagrams of specimens subjected to the same preloading program with various strain rates can be superposed to construct a master-curve. A constitutive model is developed in cyclic viscoelastoplasticity of semi crystalline polymers. A polymer is thought of as two-phase continuum composed of amorphous and crystalline regions. Both phases were treated as viscoelastoplastic media whose response was governed by different kinetic equations for evolution of plastic strains and different kinematic equations for changes in relaxation rates and relaxation spectra driven by plastic flow. Good agreement is demonstrated between the experimental data in relaxation tests under retraction and the results of numerical simulation.