The interplay of elastic anisotropy and disorder dictates many of the properties of ferroic materials, specifically martensites. We use a phase field model for ferroelastic athermal materials to study their response to an increasing external stress that couples to the strain order parameter. We show that these systems evolve through avalanches, and study the avalanche-size distribution for ferroelastic systems (large anisotropy and/or small disorder) and for the strain glass (small anisotropy and/or large disorder) using various statistical analysis techniques including the maximum likelihood method. The model predicts that in the former case the distribution is subcritical or power-law (in agreement with experimental observations) whereas in the latter case it becomes supercritical. Our results are consistent with experiments on martensitic materials and we predict specific avalanche behavior that can be tested and used as an alternative means to characterize strain glasses.