Magnetostrictive materials are known to have hysteresis with respect to both magnetic field H and mechanical stress sigma. Also, the order in which H and sigma are applied is sometimes significant. This paper proposes one generalization of the classical Preisach model in which these properties are handled by treating a combination of H and sigma as an effective field H(e) (H(z)sigma). The classical Preisach model can thus be generalized to M(f) = integral integral mu(alpha,beta) gamma-alpha-beta H(a) (H(t), sigma(t),alpha,beta) d-alpha d-beta -H(e) less-than-or-equal-to beta less-than-or-equal-to H(g) If H(e) also depends on alpha and beta. In this paper, the function H(g) (H,sigma,alpha,beta) will be assumed to have a certain simple form. It is shown that all necessary parameters can then be determined from a comparatively small number of experiments and that the above expression can be reduced to a single integral. To test the validity of the model, experiments where the two components H and sigma have been varied in many different ways have been performed on Terfenol-D and compared to results computed from the model. Some of these results are presented in the paper. This stress-dependent model is found to have an accuracy comparable to that of the classical Preisach model.