We develop a model that quantifies constitutive nonlinearities and hysteresis inherent to ferroelastic compounds, with emphasis placed on shape memory alloys. We formulate the model in two steps. First, we use the Landau theory of phase transitions to characterize the effective Gibbs free energy for both single-crystal and polycrystalline ferroelastics. The resulting nonlinear equations model equilibrium material behavior in the absence of impurities. Second, we incorporate pinning losses to account for the energy required to move domain walls across material inclusions. The full model is analogous to those developed by Jiles and Atherton for ferromagnetic compounds and Smith and Hom for ferroelectric materials. We illustrate aspects of the model through numerical simulations and comparisons with experimental stress-strain data.Keywords: Ferroelastic hysteresis; shape memory alloy; domain wall theory; ferroelastic domain; superelasticity; Landau theory of phase transitions;anyhysteretic strain; domain wall pinning; pinning sites. Public reporting burden for the collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data sources, gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any other aspect of this collection of information, including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for Information Operations and Reports, 1215 Jefferson Davis Highway, Suite 1204, Arlington VA 22202-4302. Respondents should be aware that notwithstanding any other provision of law, no person shall be subject to a penalty for failing to comply with a collection of information if it does not display a currently valid OMB control number.