The elastic and anelastic properties of a single crystal of Co-doped pnictide Ba(Fe0.957Co0.043)2As2 have been determined by resonant ultrasound spectroscopy in the frequency range 10–500 kHz, both as a function of temperature through the normalsuperconducting transition (Tc ≈ 12.5 K) and as a function of applied magnetic field up to 12.5 T. Correlation with thermal expansion, electrical resistivity, heat capacity, DC and AC magnetic data from crystals taken from the same synthetic batch has revealed the permeating influence of strain on coupling between order parameters for the ferroelastic (QE) and superconducting (QSC) transitions and on the freezing/relaxation behaviour of vortices. Elastic softening through Tc in zero field can be understood in terms of classical coupling of the order parameter with the shear strain e6, λe6Q2SC, which means that there must be a common strain mechanism for coupling of the form λQ2SCQE. At fields of ~5 T and above, this softening is masked by Debye-like stiffening and acoustic loss processes due to vortex freezing. The first loss peak may be associated with the establishment of superconductivity ... mehron ferroelastic twin walls ahead of the matrix and the second is due to the vortex liquid–vortex glass transition. Strain contrast between vortex cores and the superconducting matrix will contribute significantly to interactions of vortices both with each other and with the underlying crystal structure. These interactions imply that iron-pnictides represent a class of multiferroic superconductors in which strain-mediated coupling occurs between the multiferroic properties (ferroelasticity, antiferromagnetism) and superconductivity.