The present work reports a discrete, stress-dependent dynamic compliance spectra method which may be used to predict the mechanical response of nonlinear viscoelastic polymers during strain-defined processes. The method is based on the observation that the real and complex parts of the discrete dynamic compliance frequency components obtained from creep measurements are smooth, easily fit functions of stress. Comparisons between experimental measurements and model calculations show that the model exhibits excellent quantitative agreement with the basis creep measurements at all experimental stress levels. The model exhibits good quantitative agreement with stress relaxation measurements at moderate levels of applied strain. However, the model underestimates the experimental stress relaxation at an applied strain of 3.26%. The stress relaxation error appears to be a real material effect resulting from the different strain character of creep and stress relaxation tests. The model provides a good quantitative agreement with experimental constant strain rate measurements up to approximately 4% strain, after which the model underestimates the experimental flow stress. This effect is explained by the time dependence of the stress-activated configurational changes necessary for large strains in glassy polymers.