Many ferroelastic crystals display at sufficiently low measurement frequencies a huge elastic softening below Tc which is caused by domain wall motion. Materials range from perovskites to iron based superconductors and shape memory materials. We present a model -based on Landau-Ginzburg theory including long range elastic interaction between needle shaped ferroelastic domains -to describe the observed superelastic softening. The theory predicts that the domain wall contribution to the elastic susceptibility is different for improper and proper ferroelastic materials. A test of the theory against experimental data on SrTiO3, KMnF3, LaAlO3, La1−xNdxP5O14 and NH4HC2O4· 1 2 H2O yields excellent agreement.The macroscopic response of materials depends often on structures spanning a broad range of length and time scales. Examples of such inhomogeneous structures are precursor clusters near structural phase transitions, domains and domain walls, interfaces at first order displacive or reconstructive phase transitions, etc. They have been reported in a wide range of functional materials. Prominent cases are the enhancement of piezoelectric response [1] of nanotwinned BaTiO 3 for domain thickness below 50 nm, or the giant piezoelectric effect in ferroelectric relaxors [2] due to the high mobility of polar nanoregions [3]. Generally, domains and domain boundaries in ferroic crastals can act as structural elements for creating novel functional devices [4]. The usefulness of domain boundaries in Domain Boundary Engineering depends on the time scale of their response. Static or pinned domain boundaries can be used as functional units since they can host functional properties (ferroelectric, superconducting, etc.) which are absent in the bulk [5]. Mobile domain boundaries can lead to giant macroscopic responses, depending on the frequency of the changing external field.In the present work we will focus on ferroelastic materials, which depending on the type of coupling between the primary order parameter η and the strains ε can be classified as proper, pseudo proper, improper or coelastic ones [6]. Very often ferroelastic crystals consist of a large number of domains which are separated by domain boundaries [7]. Due to mechanical compatibility these domain boundaries should be planar with well defined orientation. However, in real crystals very often needle or dagger shaped domains appear [6]. There are still many open questions, e.g. concerning the stability of ferroelastic domains, their motion under an applied dynamic stress, the observed domain freezing [8] at sufficiently low temperature, etc.Understanding the macroscopic behaviour of multidomain crystals is important for technological applications as well as for the interpretation of seismic signals of our Earth, since domain wall motion influences the low frequency elastic and anelastic behaviour of minerals at seismic frequencies (1)(2)(3)(4)(5)(6)(7)(8)(9)(10)(11)(12)(13)(14)(15)(16)(17)(18)(19)(20). In recent years we have performed quite detailled low frequency (0.1-...