Elastic and anelastic properties of ceramic samples of multiferroic perovskites with nominal compositions across the binary join PbZr0.53Ti0.47O3–PbFe0.5Ta0.5O3 (PZT–PFT) have been assembled to create a binary phase diagram and to address the role of strain relaxation associated with their phase transitions. Structural relationships are similar to those observed previously for PbZr0.53Ti0.47O3–PbFe0.5Nb0.5O3 (PZT–PFN), but the magnitude of the tetragonal shear strain associated with the ferroelectric order parameter appears to be much smaller. This leads to relaxor character for the development of ferroelectric properties in the end member PbFe0.5Ta0.5O3. As for PZT–PFN, there appear to be two discrete instabilities rather than simply a reorientation of the electric dipole in the transition sequence cubic–tetragonal–monoclinic, and the second transition has characteristics typical of an improper ferroelastic. At intermediate compositions, the ferroelastic microstructure has strain heterogeneities on a mesoscopic length scale and, probably, also on a microscopic scale. This results in a wide anelastic freezing interval for strain-related defects rather than the freezing of discrete twin walls that would occur in a conventional ferroelastic material. In PFT, however, the acoustic loss behaviour more nearly resembles that due to freezing of conventional ferroelastic twin walls. Precursor softening of the shear modulus in both PFT and PFN does not fit with a Vogel–Fulcher description, but in PFT there is a temperature interval where the softening conforms to a power law suggestive of the role of fluctuations of the order parameter with dispersion along one branch of the Brillouin zone. Magnetic ordering appears to be coupled only weakly with a volume strain and not with shear strain but, as with multiferroic PZT–PFN perovskites, takes place within crystals which have significant strain heterogeneities on different length scales.