We review the description of ferroelastic transitions in terms of spin models. We show how one can systematically obtain a pseudo-spin Hamiltonian from the Landau energy describing the first order transition between Austenite/Martensite phases. It is shown that a Local Mean-field approximation predicts the same microstructure as the continuous Landau model in terms of strain variables. This method can be applied to a wide range of two and three dimensional transitions. We then demonstrate how quenched disorder in such pseudo-spin models yields the existence of a glass phase, characterized by the Edwards-Anderson order parameter. Our approach uses Mean-field approximation and Monte-Carlo simulations (using Zero Field Cooling/Field Cooling experiments) to study the influence of the long-range interactions. Although our model captures the salient features of a ferroelastic material in the presence of disorder, the influence of the disorder on the high symmetry austenite phase is not quite consistent with expected behavior. We examine different means of introducing disorder that can improve upon the results.