MTEX offers several ways to compute average material tensors from ODFs or EBSD data. We start by importing some EBSD data of Glaucophane and Epidote.

## Data Correction

next, we correct the data by excluding orientations with large MAD value

## Import the elastic stiffness tensors

The elastic stiffness tensor of Glaucophane was reported in Bezacier et al. 2010 (Tectonophysics) with respect to the crystal reference frame

and the density in g/cm^3

by the coefficients \(C_{ij}\) in Voigt matrix notation

The stiffness tensor in MTEX is defined via the command `stiffnessTensor`

.

The elastic stiffness tensor of Epidote was reported in Aleksandrov et al. 1974 'Velocities of elastic waves in minerals at atmospheric pressure and increasing the precision of elastic constants by means of EVM (in Russian)', Izv. Acad. Sci. USSR, Geol. Ser.10, 15-24, with respect to the crystal reference frame

and the density in g/cm^3

by the coefficients \(C_{ij}\) in Voigt matrix notation

## The Average Tensor from EBSD Data

The Voigt, Reuss, and Hill averages for all phases are computed by

The Voigt and Reuss are averaging schemes for obtaining estimates of the effective elastic constants in polycrystalline materials. The Voigt average assumes that the elastic strain field in the aggregate is constant everywhere, so that the strain in every position is equal to the macroscopic strain of the sample. CVoigt is then estimated by a volume average of local stiffness C(gi), where gi is the orientation given by a triplet of Euler angles and with orientation gi, and volume fraction V(i). This is formally described as

\( \left<T\right>^{\text{Voigt}} = \sum_{m=1}^{M} T(\mathtt{ori}_{m})\)

The Reuss average on the other hand assumes that the stress field in the aggregate is constant, so the stress in every point is set equal to the macroscopic stress. CReuss is therefore estimated by the volume average of local compliance S(gi) and can be described as

\( \left<T\right>^{\text{Reuss}} = \left[ \sum_{m=1}^{M} T(\mathtt{ori}_{m})^{-1} \right]^{-1}\)

For weakly anisotropic phases (e.g. garnet), Voigt and Reuss averages are very close to each other, but with increasing elastic anisotropy, the values of the Voigt and Reuss bounds vary considerably

The estimate of the elastic moduli of a given aggregate nevertheless should lie between the Voigt and Reuss average bounds, as the stress and strain distributions should be somewhere between the uniform strain (Voigt bound) and uniform stress.

Hill (1952) showed that the arithmetic mean of the Voigt and Reuss bounds (called Hill or Voigt-Reuss-Hill average) is very often close to the experimental values (although there is no physical justification for this behavior).

## Averaging the elastic stiffness of an aggregate based on EBSD data

for a single phase (e.g. Glaucophane) the syntax is

## ODF Estimation

Next, we estimate an ODF for the Epidote phase

## The Average Tensor from an ODF

The Voigt, Reuss, and Hill averages for the above ODF are computed by

To visualize the polycrystalline Glaucophane wave velocities we can use the command `plotSeismicVelocities`

More details on averaging the seismic properties considering the modal composition of different phases can be found in here