Tables for
Volume D
Physical properties of crystals
Edited by A. Authier

International Tables for Crystallography (2013). Vol. D, ch. 1.3, pp. 72-100
doi: 10.1107/97809553602060000902

Chapter 1.3. Elastic properties

A. Authiera* and A. Zarembowitchb

aInstitut de Minéralogie et de Physique des Milieux Condensés, 4 Place Jussieu, 75005 Paris, France, and bLaboratoire de Physique des Milieux Condensés, Université P. et M. Curie, 75252 Paris CEDEX 05, France
Correspondence e-mail:

In this chapter, the strain and stress tensors are defined and their main properties are derived. The elastic tensors, elastic stiffnesses and elastic compliances are then introduced. Their variation with orientation, depending on the crystal class, is given in the case of Young's modulus. The next part is devoted to the propagation of waves in continuous media (linear dynamic elasticity and the Christoffel matrix); the relation between the velocity and the elastic constants is given for the cubic, hexagonal and tetragonal classes. The experimental determination of elastic constants and their pressure and temperature dependence are discussed in separate sections. The last two sections of the chapter concerns nonlinear elasticity (second and higher-order elastic constants) and nonlinear dynamical elasticity.

1.3.1. Strain tensor

| top | pdf | Introduction, the notion of strain field

| top | pdf |

Let us consider a medium that undergoes a deformation. This means that the various points of the medium are displaced with respect to one another. Geometrical transformations of the medium that reduce to a translation of the medium as a whole will therefore not be considered. We may then suppose that there is an invariant point, O, whose position one can always return to by a suitable translation. A point P, with position vector [{\bf OP} = {\bf r}], is displaced to the neighbouring point P′ by the deformation defined by [{\bf PP}' = {\bf u} ({\bf r}). ]The displacement vector [{\bf u}({\bf r})] constitutes a vector field. It is not a uniform field, unless the deformation reduces to a translation of the whole body, which is incompatible with the hypothesis that the medium undergoes a deformation. Let Q be a point that is near P before the deformation (Fig.[link]). Then one can write [\hbox{d}{\bf r} = {\bf PQ}\semi \quad {\bf r} + \hbox{d}{\bf r} = {\bf OQ}. ]


Figure | top | pdf |

Displacement vector, [{\bf u}({\bf r})].

After the deformation, Q is displaced to Q′ defined by [{\bf QQ}' = {\bf u} ({\bf r} + \hbox{d}{\bf r}). ]

In a deformation, it is more interesting in general to analyse the local, or relative, deformation than the absolute displacement. The relative displacement is given by comparing the vectors [{\bf P}'{\bf Q}' = \hbox{d}{\bf r}' ] and PQ. Thus, one has [{\bf P}'{\bf Q}' = {\bf P}'{\bf P} + {\bf PQ} + {\bf QQ}'. ]Let us set [\left. \matrix{\hbox{d} {\bf r}' = \hbox{d} {\bf r} + {\bf u} ({\bf r} + \hbox{d} {\bf r}) - {\bf u} ({\bf r}) \hfill \cr \noalign{\vskip5pt} \hbox{d}{\bf u} = {\bf u} ({\bf r} + \hbox{d} {\bf r}) - {\bf u} ({\bf r}) = \hbox{d} {\bf r}' - \hbox{d} {\bf r}.\cr}\right\} \eqno( ]Replacing [{\bf u}({\bf r} + \hbox{d} {\bf r})] by its expansion up to the first term gives [\left. \matrix{\hbox{d}u_{i} = {\displaystyle{\partial u_{i} \over \partial x_{j}}}\ \hbox{d}x_{j} \hfill \cr \noalign{\vskip5pt} \hbox{d}x'_{i} = \hbox{d}x_{i} + {\displaystyle{\partial u_{i} \over \partial x_{j}}}\ \hbox{d}x_{j}.\cr}\right\} \eqno( ]

If we assume the Einstein convention (see Section[link] ), there is summation over j in ([link] and ([link]. We shall further assume orthonormal coordinates throughout Chapter 1.3; variance is therefore not apparent and the positions of the indices have no meaning; the Einstein convention then only assumes repetition of a dummy index. The elements [\hbox{d}x_{i}] and [\hbox{d}x'_{i}] are the components of dr and [\hbox{d}{\bf r}'], respectively. Let us put [M_{ij} = \partial u_{i}/\partial x_{j}\semi \ B_{ij} = M_{ij}+ \delta _{ij}, ]where [\delta _{ij}] represents the Kronecker symbol; the [\delta_{ij}]'s are the components of matrix unity, I. The expressions ([link] can also be written using matrices M and B: [\left. \matrix{\hbox{d}u_{i} = M_{ij} \hbox{d}x_{j}\cr \noalign{\vskip5pt} \hbox{d}x'_{i} = B_{ij}\hbox{d}x_{j}.\cr}\right\} \eqno( ]The components of the tensor [M_{ij}] are nonzero, unless, as mentioned earlier, the deformation reduces to a simple translation. Two cases in particular are of interest and will be discussed in turn:

  • (i) The components [M_{ij}] are constants. In this case, the deformation is homogeneous.

  • (ii) The components [M_{ij}] are variables but are small compared with unity. This is the practical case to which we shall limit ourselves in considering an inhomogeneous deformation. Homogeneous deformation

| top | pdf |

If the components [M_{ij}] are constants, equations ([link] can be integrated directly. They become, to a translation, [\left. \matrix{u_{i} = M_{ij} x_{j}\cr \noalign{\vskip5pt} x'_{i} = B_{ij}x_{j}.\cr}\right\} \eqno( ] Fundamental property of the homogeneous deformation

| top | pdf |

The fundamental property of the homogeneous deformation results from the fact that equations ([link] are linear: a plane before the deformation remains a plane afterwards, a crystal lattice remains a lattice. Thermal expansion is a homogeneous deformation (see Chapter 1.4[link] ). Spontaneous strain

| top | pdf |

Some crystals present a twin microstructure that is seen to change when the crystals are gently squeezed. At rest, the domains can have one of two different possible orientations and the influence of an applied stress is to switch them from one orientation to the other. If one measures the shape of the crystal lattice (the strain of the lattice) as a function of the applied stress, one obtains an elastic hysteresis loop analogous to the magnetic or electric hysteresis loops observed in ferromagnetic or ferroelectric crystals. For this reason, these materials are called ferroelastic (see Chapters 3.1[link] to [link] 3.3[link] and Salje, 1990[link]). The strain associated with one of the two possible shapes of the crystal when no stress is applied is called the macroscopic spontaneous strain. Cubic dilatation

| top | pdf |

Let [{\bf e}_{i}] be the basis vectors before deformation. On account of the deformation, they are transformed into the three vectors [{\bf e}'_{i} = B_{ij}{\bf e}_{j}.]The parallelepiped formed by these three vectors has a volume V′ given by [V' = ({\bf e}'_{1},{\bf e}'_{2},{\bf e}'_{3}) = \Delta (B) ({\bf e}_{1},{\bf e}_{2},{\bf e}_{3}) = \Delta (B)V, ]where [\Delta (B)] is the determinant associated with matrix B, V is the volume before deformation and [({\bf e}_{1},{\bf e}_{2},{\bf e}_{3}) = ({\bf e}_{1} \wedge {\bf e}_{2})\cdot {\bf e}_{3} ]represents a triple scalar product.

The relative variation of the volume is [{V' - V \over V} = \Delta (B) - 1. \eqno(]It is what one calls the cubic dilatation. [\Delta (B)] gives directly the volume of the parallelepiped that is formed from the three vectors obtained in the deformation when starting from vectors forming an orthonormal base. Expression of any homogeneous deformation as the product of a pure rotation and a pure deformation

| top | pdf |

  • (i) Pure rotation: It is isometric. The moduli of the vectors remain unchanged and one direction remains invariant, the axis of rotation. The matrix B is unitary: [BB^{T}= 1.]

  • (ii) Pure deformation: This is a deformation in which three orthogonal directions remain invariant. It can be shown that B is a symmetric matrix: [B = B^{T}.]The three invariant directions are those of the eigenvectors of the matrix; it is known in effect that the eigenvectors of a symmetric matrix are real.

  • (iii) Arbitrary deformation: the matrix B, representing an arbitrary deformation, can always be put into the form of the product of a unitary matrix [B_{1}], representing a pure rotation, and a symmetric matrix [B_{2}], representing a pure deformation. Let us put [B = B_{1}B_{2}]and consider the transpose matrix of B: [B^{T}= B_{2}^T B_{1}^T = B_{2}\left(B_{1}\right)^{-1}. ]The product [B^{T}B] is equal to [B^{T}B = \left(B_{2}\right)^{2}.]This shows that we can determine [B_{2}] and therefore [B_{1}] from B. Quadric of elongations

| top | pdf |

Let us project the displacement vector [{\bf u}({\bf r}) ] on the position vector OP (Fig.[link]), and let [u_{r}] be this projection. The elongation is the quantity defined by [{u_{r} \over r} = {{\bf u \cdot r} \over r^2} = {M_{ij}x_ix_j \over r^2}, ]where [x_{1}], [x_{2}], [x_{3}] are the components of r. The elongation is the relative variation of the length of the vector r in the deformation. Let A and S be the antisymmetric and symmetric parts of M, respectively: [A = {M - M^T \over 2}; \quad S = {M + M^T \over 2}. ]


Figure | top | pdf |

Elongation, [u_{r}/r].

Only the symmetric part of M occurs in the expression of the elongation: [{u_{r} \over r} = {S_{ij}x_ix_j \over r^2}. \eqno( ]

The geometrical study of the elongation as a function of the direction of r is facilitated by introducing the quadric associated with M: [S_{ij} y_iy_j = \varepsilon, \eqno( ]where [\varepsilon] is a constant. This quadric is called the quadric of elongations, Q. S is a symmetric matrix with three real orthogonal eigenvectors and three real eigenvalues, [\lambda_{1}], [\lambda_{2}], [\lambda_{3}]. If it is referred to these axes, equation ([link] is reduced to [\lambda_{1}\left(y_{1}\right)^{2}\lambda_{2}\left(y_{2}\right)^{2}\lambda_{3} \left(y_{3}\right)^{2} = \varepsilon. ]

One can discuss the form of the quadric according to the sign of the eigenvalues [\lambda_{i}]:

  • (i) [\lambda_{1}], [\lambda_{2}], [\lambda_{3}] have the same sign, and the sign of [\varepsilon ]. The quadric is an ellipsoid (Fig.[link]). One chooses [\varepsilon = +1] or [\varepsilon = - 1 ], depending on the sign of the eigenvalues.


    Figure | top | pdf |

    Quadric of elongations. The displacement vector, [{\bf u}({\bf r})], at P in the deformed medium is parallel to the normal to the quadric at the intersection, p, of OP with the quadric. (a) The eigenvalues all have the same sign, the quadric is an ellipsoid. (b) The eigenvalues have mixed signs, the quadric is a hyperboloid with either one sheet (shaded in light grey) or two sheets (shaded in dark grey), depending on the sign of the constant [\varepsilon] [see equation ([link]]; the cone asymptote is represented in medium grey. For a practical application, see Fig.[link] .

  • (ii) [\lambda_{1}], [\lambda_{2} ], [\lambda_{3}] are of mixed signs: one of them is of opposite sign to the other two. One takes [\varepsilon = \pm 1 ]. The corresponding quadric is a hyperboloid whose asymptote is the cone [S_{ij} y_iy_j = 0.]

    According to the sign of [\varepsilon], the hyperboloid will have one sheet outside the cone or two sheets inside the cone (Fig.[link]). If we wish to be able to consider any direction of the position vector r in space, it is necessary to take into account the two quadrics.

In order to follow the variations of the elongation [u_{r}/r ] with the orientation of the position vector, one associates with r a vector y, which is parallel to it and is defined by [{\bf y} = {\bf r}/k; \quad {\bf r} = k {\bf y},]where k is a constant. It can be seen that, in accordance with ([link] and ([link], the expression of the elongation in terms of y is [u_{r}/r = \varepsilon/y^{2}. ]

Thus, the elongation is inversely proportional to the square of the radius vector of the quadric of elongations parallel to OP. In practice, it is necessary to look for the intersection p of the parallel to OP drawn from the centre O of the quadric of elongations (Fig.[link]):

  • (i) The eigenvalues all have the same sign; the quadric Q is an ellipsoid: the elongation has the same sign in all directions in space, positive for [\varepsilon = +1] and negative for [\varepsilon = -1].

  • (ii) The eigenvalues have different signs; two quadrics are to be taken into account: the hyperboloids corresponding, respectively, to [\varepsilon = \pm 1]. The sign of the elongation is different according to whether the direction under consideration is outside or inside the asymptotic cone and intersects one or the other of the two hyperboloids.

Equally, one can connect the displacement vector [{\bf u}({\bf r}) ] directly with the quadric Q. Using the bilinear form [f({\bf y}) = M_{ij}y_{i}y_{j},]the gradient of [f({\bf y})], [\boldnabla (f)], has as components [\partial f/\partial y^{i} = M_{ij}y_{j} = u_{i}.]

One recognizes the components of the displacement vector u, which is therefore parallel to the normal to the quadric Q at the extremity of the radius vector Op parallel to r.

The directions of the principal axes of Q correspond to the extremal values of y, i.e. to the stationary values (maximal or minimal) of the elongation. These values are the principal elongations.

If the deformation is a pure rotation [\eqalign{&B = \pmatrix{\cos \theta & \sin \theta & 0\cr - \sin \theta & \cos \theta & 0\cr 0 & 0 &1\cr}, \cr\noalign{\vskip5.5pt}\qquad &M = \pmatrix{\cos \theta - 1 & \sin \theta & 0\cr - \sin \theta & \cos \theta - 1 &0\cr 0 & 0 &0\cr}.} ]Hence we have [M_{ij}y_{i}y_{j} = (\cos\theta - 1)\left(y_{1}- y_{2}\right) = \varepsilon. ]

The quadric Q is a cylinder of revolution having the axis of rotation as axis. Arbitrary but small deformations

| top | pdf | Definition of the strain tensor

| top | pdf |

If the deformation is small but arbitrary, i.e. if the products of two or more components of [M_{ij}] can be neglected with respect to unity, one can describe the deformation locally as a homogeneous asymptotic deformation. As was shown in Section[link], it can be put in the form of the product of a pure deformation corresponding to the symmetric part of [M_{ij}], [S_{ij} ], and a pure rotation corresponding to the asymmetric part, [A_{ij}]: [\left.\eqalign{ S_{ij} & = S_{ji} = \textstyle{1 \over 2} \left( \displaystyle{\partial u_i \over \partial x_j} + {\partial u_j \over \partial x_i} \right)\cr A_{ij} & = - A_{ji} =\ \textstyle{1 \over 2} \left( \displaystyle{\partial u_i \over \partial x_j} - {\partial u_j \over \partial x_i} \right).\cr }\right\} \eqno{(} ]Matrix B can be written [B = I + A + S, ]where I is the matrix identity. As the coefficients [\partial u_{i}/ \partial x_{j}] of [M_{ij}] are small, one can neglect the product [A\times S] and one has [B = (I + A)(I + S).][(I + S)] is a symmetric matrix that represents a pure deformation. [(I + A)] is an antisymmetric unitary matrix and, since A is small, [(I + A)^{-1}= (I - A).]Thus, [(I + A) ] represents a rotation. The axis of rotation is parallel to the vector with coordinates [\left.\eqalign{\Omega_{1} &= \textstyle{{1 \over 2}} \left( \displaystyle{\partial u_3 \over \partial x_2} + {\partial u_2 \over \partial x_3} \right) = A_{32} \cr \Omega_{2} & =\ \textstyle{1 \over 2} \left( \displaystyle{\partial u_1 \over \partial x_3} + {\partial u_3 \over \partial x_1} \right) = A_{13}\cr \Omega_{3} & = \textstyle{1 \over 2} \left( \displaystyle{\partial u_2 \over \partial x_1} + {\partial u_1 \over \partial x_2} \right) = A_{21},\cr}\right\} ]which is an eigenvector of [(I + A)]. The magnitude of the rotation is equal to the modulus of this vector.

In general, one is only interested in the pure deformation, i.e. in the form of the deformed object. Thus, one only wishes to know the quantities [(I + S)] and the symmetric part of M. It is this symmetric part that is called the deformation tensor or the strain tensor. It is very convenient for applications to use the simplified notation due to Voigt:[\matrix{S_{1} = \displaystyle{\partial u_{1} \over \partial x_{1}}\semi\hfill & S_{2}= \displaystyle{\partial u_{2} \over \partial x_{2}}\semi\hfill & S_{3}= \displaystyle{\partial u_{3} \over \partial x_{3}}\semi \hfill\cr S_{4} = \displaystyle{\partial u_{3} \over \partial x_{2}} + \displaystyle{\partial u_{2} \over \partial x_{3}}\semi \hfill &S_{5}= \displaystyle{\partial u_{3} \over \partial x_{1}} + {\partial u_{1} \over \partial x_{3}}\semi\hfill & S_{6}= \displaystyle{\partial u_{2} \over \partial x_{1}} + {\partial u_{1} \over \partial x_{2}}.\cr} ]One may note that [\matrix{S_{1}= S_{11}\semi\hfill &S_{2}= S_{22}\semi \hfill &S_{3}=S_{33}\semi\hfill\cr S_{4}=S_{23}+ S_{32}\semi \hfill&S_{5}=S_{31}+S_{13}\semi \hfill&S_{6}=S_{12} + S_{21}.\hfill\cr} ]The Voigt strain matrix S is of the form [\pmatrix{S_{1} &S_{6} &S_{5}\cr S_{6} &S_{2} &S_{4}\cr S_{5} &S_{4} &S_{3}\cr}. ] Geometrical interpretation of the coefficients of the strain tensor

| top | pdf |

Let us consider an orthonormal system of axes with centre P. We remove nothing from the generality of the following by limiting ourselves to a planar problem and assuming that point P′ to which P goes in the deformation lies in the plane [x_{1}Px_{2} ] (Fig.[link]). Let us consider two neighbouring points, Q and R, lying on axes [Px_{1}] and [Px_{2}], respectively ([PQ = \hbox{d} x_{1} ], [PR = \hbox{d} x_{2}]). In the deformation, they go to points Q′ and R′ defined by [\eqalign{{\bf QQ'}&: \left\{\matrix {\hbox{d}x'_{1} = \hbox{d}x_{1} + \left(\partial u_{1} / \partial x_{1}\right)\hbox{d}x_{1}\hfill\cr \hbox{d}x'_{2} = \left(\partial u_{2} / \partial x_{1}\right)\hbox{d}x_{1}\hfill\cr \hbox{d}x'_{3} = 0\hfill\cr}\right. \cr {\bf RR'}&: \left\{\matrix {\hbox{d}x'_{1} = \left(\partial u_{1} / \partial x_{2}\right)\hbox{d}x_{2}\hfill\cr \hbox{d}x'_{2} = \hbox{d}x_{2} + \left(\partial u_{2} / \partial x_{2}\right)\hbox{d}x_{2}\hfill\cr \hbox{d}x'_{3} = 0.\hfill\cr}\right.} ]


Figure | top | pdf |

Geometrical interpretation of the components of the strain tensor. [Ox_{1}], [Ox_{2}], [Ox_{3} ]: axes before deformation; [Ox_{1}'], [Ox_{2}'], [Ox_{3}']: axes after deformation.

As the coefficients [\partial u_{i}/\partial x_{j}] are small, the lengths of PQ′ and PR′ are hardly different from PQ and PR, respectively, and the elongations in the directions [Px_{1}] and [Px_{2}] are [\eqalign{{P' Q' - PQ \over PQ} &= { \hbox{d}x'_{1} - \hbox{d}x_{1} \over \hbox{d}x_{1}} = {\partial u_{1} \over \partial x_{1}} = S_{1}\cr {P' R' - PR \over PR} &= { \hbox{d}x'_{2} - \hbox{d}x_{2} \over \hbox{d}x_{2}} = {\partial u_{2} \over \partial x_{2}} = S_{2}.\cr} ]

The components [S_{1}], [S_{2}], [S_{3}] of the principal diagonal of the Voigt matrix can then be interpreted as the elongations in the three directions [Px_{1}], [Px_{2}] and [Px_{3}]. The angles α and β between PQ and PQ′, and PR and PR′, respectively, are given in the same way by [\alpha = \hbox{d}x'_{2}/\hbox{d}x_{1} = \partial u_{2}/\partial x_{1}\semi\ \ \beta =\hbox{d}x'_{1}/dx_{2} = \partial u_{1}/\partial x_{2}. ]One sees that the coefficient [S_{6}] of Voigt's matrix is therefore [S_{6}= {\partial u_{2} \over \partial x_{1}} + {\partial u_{1} \over \partial x_{2}} = \alpha + \beta. ]The angle [\alpha + \beta] is equal to the difference between angles [{\bf PQ}\wedge {\bf PR}] before deformation and [{\bf P'' Q''}\wedge {\bf P' R'}] after deformation. The nondiagonal terms of the Voigt matrix therefore represent the shears in the planes parallel to [Px_{1}], [Px_{2}] and [Px_{3}], respectively.

To summarize, if one considers a small cube before deformation, it becomes after deformation an arbitrary parallelepiped; the relative elongations of the three sides are given by the diagonal terms of the strain tensor and the variation of the angles by its nondiagonal terms.

The cubic dilatation ([link] is [\Delta (B) - 1 = S_{1}+ S_{2}+ S_{3}](taking into account the fact that the coefficients [S_{ij}] are small). Particular components of the deformation

| top | pdf | Simple elongation

| top | pdf |

Matrix M has only one coefficient, [e_{1}], and reduces to (Fig.[link]) [\pmatrix{e_{1} &0 &0\cr 0 &0 &0\cr 0 &0 &0\cr}.]The quadric of elongations is reduced to two parallel planes, perpendicular to [Ox_{1}], with the equation [x_{1} = \pm 1/\sqrt{\vert e \vert} ].


Figure | top | pdf |

Special deformations. The state after deformation is represented by a dashed line. (a) Simple elongation; (b) pure shear; (c) simple shear. Pure shear

| top | pdf |

This is a pure deformation (without rotation) consisting of the superposition of two simple elongations along two perpendicular directions (Fig.[link]) and such that there is no change of volume (the cubic dilatation is zero): [\pmatrix{e_{1} &0 &0\cr 0 &-e_{1} &0\cr 0 &0 &0\cr}. ]The quadric of elongations is a hyperbolic cylinder. Simple shear

| top | pdf |

Matrix [M_{ij}] has one coefficient only, a shear (Fig.[link]): [\pmatrix{0 &s &0\cr 0 &0 &0\cr 0 &0 &0\cr}. ]The matrix is not symmetrical, as it contains a component of rotation. Thus we have [\left.\eqalign{x'_{1} &= x_{1} + sx_{2}\cr x'_{2} &= x_{2}\cr x'_{3} &= x_{3}.\cr}\right\} ]One can show that the deformation is a pure shear associated with a rotation around [Ox_{3}].

1.3.2. Stress tensor

| top | pdf | General conditions of equilibrium of a solid

| top | pdf |

Let us consider a solid C, in movement or not, with a mass distribution defined by a specific mass ρ at each point. There are two types of force that are manifested in the interior of this solid.

  • (i) Body forces (or mass forces), which one can write in the form [{\bf F} \hbox{ d}m = {\bf F} \rho \hbox{ d}\tau, ]where dτ is a volume element and dm a mass element. Gravity forces or inertial forces are examples of body forces. One can also envisage body torques (or volume couples), which can arise, for example, from magnetic or electric actions but which will be seen to be neglected in practice.

  • (ii) Surface forces or stresses. Let us imagine a cut in the solid along a surface element dσ of normal n (Fig.[link]). The two lips of the cut that were in equilibrium are now subjected to equal and opposite forces, R and [{\bf R}' = -{\bf R}], which will tend to separate or draw together these two lips. One admits that, when the area element dσ tends towards zero, the ratio [{\bf R}/\hbox{d}\sigma] tends towards a finite limit, [{\bf T}_{n}], which is called stress. It is a force per unit area of surface, homogeneous to a pressure. It will be considered as positive if it is oriented towards the same side of the surface-area element dσ as the normal n and negative in the other case. The choice of the orientation of n is arbitrary. The pressure in a liquid is defined in a similar way but its magnitude is independent of the orientation of n and its direction is always parallel to n. On the other hand, in a solid the constraint [{\bf T}_{n}] applied to a surface element is not necessarily normal to the latter and the magnitude and the orientation with respect to the normal change when the orientation of n changes. A stress is said to be homogeneous if the force per unit area acting on a surface element of given orientation and given shape is independent of the position of the element in the body. Other stresses are inhomogeneous. Pressure is represented by a scalar, and stress by a rank-two tensor, which will be defined in Section[link].


    Figure | top | pdf |

    Definition of stress: it is the limit of R dσ when the surface element dσ tends towards zero. R and R′ are the forces to which the two lips of the small surface element cut within the medium are subjected.

Now consider a volume V within the solid C and the surface S which surrounds it (Fig.[link]). Among the influences that are exterior to V, we distinguish those that are external to the solid C and those that are internal. The first are translated by the body forces, eventually by volume couples. The second are translated by the local contact forces of the part external to V on the internal part; they are represented by a surface density of forces, i.e. by the stresses [{\bf T}_{n}] that depend only on the point Q of the surface S where they are applied and on the orientation of the normal n of this surface at this point. If two surfaces S and S′ are tangents at the same point Q, the same stress acts at the point of contact between them. The equilibrium of the volume V requires:

  • (i) For the resultant of the applied forces and the inertial forces: [\textstyle\int \int\limits_{S}\displaystyle {\bf T}_{n} \hbox{ d}\sigma + \textstyle\int \int\int\limits_{V}\displaystyle {\bf F}\rho \hbox{ d}\tau = {\hbox{d}\over \hbox{d}t} \left\{\textstyle\int \int \int\limits_{V}\displaystyle {\bf v} \hbox{ d}\tau \right\}. \eqno( ]


    Figure | top | pdf |

    Stress, [{\bf T}_{n}], applied to the surface of an internal volume.

  • (ii) For the resultant moment: [\textstyle\int \int\limits_{S}\displaystyle{\bf OQ} \wedge {\bf T}_{n} \hbox{ d}\sigma + \textstyle\int \int\int\limits_{V}\displaystyle {\bf OP} \wedge {\bf F} \rho \hbox{ d}\tau = {\hbox{ d}\over \hbox{ d}t}\left\{\textstyle\int \int \int\limits_{V}\displaystyle {\bf OP} \wedge {\bf v}\hbox{ d}\tau \right\}, \eqno( ]where Q is a point on the surface S, P a point in the volume V and v the velocity of the volume element dτ.

The equilibrium of the solid C requires that:

  • (i) there are no stresses applied on its surface and

  • (ii) the above conditions are satisfied for any volume V within the solid C. Definition of the stress tensor

| top | pdf |

Using the condition on the resultant of forces, it is possible to show that the components of the stress [{\bf T}_{n}] can be determined from the knowledge of the orientation of the normal n and of the components of a rank-two tensor. Let P be a point situated inside volume V, [Px_{1}], [Px_{2}] and [Px_{3}] three orthonormal axes, and consider a plane of arbitrary orientation that cuts the three axes at Q, R and S, respectively (Fig.[link]). The small volume element PQRS is limited by four surfaces to which stresses are applied. The normals to the surfaces PRS, PSQ and PQR will be assumed to be directed towards the interior of the small volume. By contrast, for reasons that will become apparent later, the normal n applied to the surface QRS will be oriented towards the exterior. The corresponding applied forces are thus given in Table[link]. The volume PQRS is subjected to five forces: the forces applied to each surface and the resultant of the volume forces and the inertial forces. The equilibrium of the small volume requires that the resultant of these forces be equal to zero and one can write [-{\bf T}_{n} \hbox{ d}\sigma + {\bf T}_{1} \hbox{ d}\sigma_{1} + {\bf T}_{2}\hbox{ d}\sigma_{2} + {\bf T}_{3}\hbox{ d}\sigma_{3} + {\bf F}\rho \hbox{ d}\tau = 0 ](including the inertial forces in the volume forces).

Table| top | pdf |
Stresses applied to the faces surrounding a volume element

[\alpha_{1}], [\alpha_{2}] and [\alpha_{3}] are the direction cosines of the normal n to the small surface QRS.

FaceAreaApplied stressApplied force
QRS [-{\bf T}_{n}] [-{\bf T}_{n}\hbox{ d}\sigma ]
PRS [\hbox{d}\sigma_{1} = \alpha_{1} \hbox{ d}\sigma ] [{\bf T}_{1}] [{\bf T}_{n} \hbox{ d}\sigma_{1} ]
PSQ [\hbox{d}\sigma_{2} = \alpha_{2} \hbox{ d}\sigma ] [{\bf T}_{2}] [{\bf T}_{n} \hbox{ d}\sigma_{2} ]
PQR [\hbox{d}\sigma_{3} = \alpha_{3} \hbox{ d}\sigma ] [{\bf T}_{3}] [{\bf T}_{n} \hbox{ d}\sigma_{3} ]

Figure | top | pdf |

Equilibrium of a small volume element.

As long as the surface element dσ is finite, however small, it is possible to divide both terms of the equation by it. If one introduces the direction cosines, [\alpha_{i}], the equation becomes [-{\bf T}_{n} + {\bf T}_{1}\hbox{ d}\alpha_{1} + {\bf T}_{2}\hbox{ d}\alpha_{2} + {\bf T}_{3}\hbox{ d}\alpha_{3} + {\bf F}\rho \hbox{ d}\tau/\hbox{d}\sigma= 0. ]When dσ tends to zero, the ratio [\hbox{d}\sigma/\hbox{d}\tau ] tends towards zero at the same time and may be neglected. The relation then becomes [{\bf T}_{n} = {\bf T}_{i}\alpha^{i}. \eqno( ]This relation is called the Cauchy relation, which allows the stress [{\bf T}_{n}] to be expressed as a function of the stresses [{\bf T}_{1}], [{\bf T}_{2}] and [{\bf T}_{3}] that are applied to the three faces perpendicular to the axes, [Px_{1}], [Px_{2}] and [Px_{3}]. Let us project this relation onto these three axes: [T_{nj} = T_{ij}\alpha_{i}. \eqno(]The nine components [T_{ij}] are, by definition, the components of the stress tensor. In order to check that they are indeed the components of a tensor, it suffices to make the contracted product of each side of ([link] by any vector [x_{i}]: the left-hand side is a scalar product and the right-hand side a bilinear form. The [T_{ij}]'s are therefore the components of a tensor. The index to the far left indicates the face to which the stress is applied (normal to the [x_{1}], [x_{2} ] or [x_{3}] axis), while the second one indicates on which axis the stress is projected. Condition of continuity

| top | pdf |

Let us return to equation ([link] expressing the equilibrium condition for the resultant of the forces. By replacing [{\bf T}_{n}] by the expression ([link], we get, after projection on the three axes, [\textstyle\int \int\limits_{S}\displaystyle T_{ij}\hbox{ d}\sigma_i + \textstyle\int\int\int\limits_{V}\displaystyle F_{j} \rho \hbox{ d} \tau = 0, ]where [\hbox{d}\sigma_{i} = \alpha_{i}\hbox{ d}\sigma] and the inertial forces are included in the volume forces. Applying Green's theorem to the first integral, we have [\textstyle\int\int\limits_{S}\displaystyle T_{ij}\ \hbox{d}\sigma_{i} = \textstyle\int\int\int\limits_{V}\displaystyle \left[\partial T_{ij}/\partial x_{i} \right] \hbox{d}\tau. ]

The equilibrium condition now becomes [\textstyle\int\int\int\limits_{V}\displaystyle \left[\partial T_{ij}/\partial x_{i} + F_{j} \rho \right] \hbox{d} \tau = 0. ]In order that this relation applies to any volume V, the expression under the integral must be equal to zero, [\partial T_{ij}/\partial x_{i} + F_{j} \rho = 0, \eqno{(} ]or, if one includes explicitly the inertial forces, [\partial T_{ij}/\partial x_{i} + F_{j} \rho = - \rho \ \partial^{2} x_{j} / \partial t^{2}. \eqno{(} ]This is the condition of continuity or of conservation. It expresses how constraints propagate throughout the solid. This is how the cohesion of the solid is ensured. The resolution of any elastic problem requires solving this equation in terms of the particular boundary conditions of that problem. Symmetry of the stress tensor

| top | pdf |

Let us now consider the equilibrium condition ([link] relative to the resultant moment. After projection on the three axes, and using the Cartesian expression ( of the vectorial products, we obtain [\textstyle\int\int\limits_{S}\displaystyle \textstyle{1 \over 2}\displaystyle \varepsilon_{ijk} (x_{i} T_{lj} - x_{j} T_{li}) \hbox{d} \sigma_{l} + \textstyle\int\int\int\limits_{V}\displaystyle \left[\textstyle{1\over 2}\displaystyle \varepsilon_{ijk} \rho (x_{i} F_{j} - x_{j} F_{i}) + \Gamma_{k}\right] \hbox{d}\tau = 0. ](including the inertial forces in the volume forces). [\varepsilon_{ijk}] is the permutation tensor. Applying Green's theorem to the first integral and putting the two terms together gives[\int\int\int\limits_{V}\displaystyle \left\{\textstyle{1\over 2}\displaystyle \varepsilon_{ijk} \left[{\partial \over \partial x_{l}} (x_{i} T_{lj} - x_{j} T_{li}) + \rho (x_{i} F_{j} - x_{j} F_{i})\right] + \Gamma_{k} \right\} \hbox{d} \tau = 0. ]

In order that this relation applies to any volume V within the solid C, we must have [\textstyle{1\over 2}\displaystyle \varepsilon_{ijk} \left[{\partial \over \partial x_{l}} (x_{i} T_{lj} - x_{j} T_{li})\right] + \Gamma_{k} = 0 ]or [\textstyle{1\over 2}\displaystyle \varepsilon_{ijk} \left[x_{i} \left({\partial T_{lj}\over \partial x_{l}} + F_{j} \rho \right) - x_{j} \left({\partial T_{li}\over \partial x_{l}} + F_{i} \rho \right) + T_{ij} - T_{ji} \right] + \Gamma_{k} = 0. ]

Taking into account the continuity condition ([link], this equation reduces to [\textstyle{1\over 2}\displaystyle \varepsilon_{ijk} [T_{ij} - T_{ji}]+ \Gamma_{k} = 0. ]

A volume couple can occur for instance in the case of a magnetic or an electric field acting on a body that locally possesses magnetic or electric moments. In general, apart from very rare cases, one can ignore these volume couples. One can then deduce that the stress tensor is symmetrical: [T_{ij} - T_{ji} = 0.]

This result can be recovered by applying the relation ([link] to a small volume in the form of an elementary parallelepiped, thus illustrating the demonstration using Green's theorem but giving insight into the action of the constraints. Consider a rectangular parallelepiped, of sides [2 \Delta x_{1}], [2 \Delta x_{2}] and [2 \Delta x_{3}], with centre P at the origin of an orthonormal system whose axes [Px_{1}], [Px_{2}] and [Px_{3}] are normal to the sides of the parallelepiped (Fig.[link]). In order that the resultant moment with respect to a point be zero, it is necessary that the resultant moments with respect to three axes concurrent in this point are zero. Let us write for instance that the resultant moment with respect to the axis [Px_{3}] is zero. We note that the constraints applied to the faces perpendicular to [Px_{3}] do not give rise to a moment and neither do the components [T_{11}], [T_{13}], [T_{22} ] and [T_{23}] of the constraints applied to the faces normal to [Px_{1}] and [Px_{2}] (Fig.[link]). The components [T_{12}] and [T_{21}] alone have a nonzero moment.


Figure | top | pdf |

Symmetry of the stress tensor: the moments of the couples applied to a parallelepiped compensate each other.

For face 1, the constraint is [T_{12} + (\partial T_{12}/ \partial x_{1})\Delta x_{1} ] if [T_{12}] is the magnitude of the constraint at P. The force applied at face 1 is [\left[T_{12} + {\partial T_{12}\over \partial x_{1}}\Delta x_{1}\right]4 \Delta x_{2} \Delta x_{3} ]and its moment is[\left[T_{12} + {\partial T_{12}\over \partial x_{1}}\Delta x_{1}\right]4 \Delta x_{2} \Delta x_{3} \Delta x_{1}. ]

Similarly, the moments of the force on the other faces are[\displaylines{\hbox{Face 1}':\quad-\left[T_{12} + {\partial T_{12}\over \partial x_{1}}(-\Delta x_{1})\right]4 \Delta x_{2} \Delta x_{3} (-\Delta x_{1})\semi\hfill\cr \hbox{Face 2}\phantom{'}:\quad\phantom{-}\left[T_{21} + {\partial T_{21}\over \partial x_{2}}\Delta x_{2})\right]4 \Delta x_{1} \Delta x_{3} \Delta x_{2};\hfill\cr \hbox{Face 2}':\quad-\left[T_{21} + {\partial T_{21}\over \partial x_{2}}(-\Delta x_{2})\right]4 \Delta x_{1} \Delta x_{3} (-\Delta x_{2}).\hfill} ]

Noting further that the moments applied to the faces 1 and 1′ are of the same sense, and that those applied to faces 2 and 2′ are of the opposite sense, we can state that the resultant moment is [[T_{12} - T_{21}]8 \Delta x_{1} \Delta x_{2} \Delta x_{3} = [T_{12} - T_{21}]\Delta \tau, ]where [8 \Delta x_{1} \Delta x_{2} \Delta x_{3} = \Delta \tau ] is the volume of the small parallelepiped. The resultant moment per unit volume, taking into account the couples in volume, is therefore [T_{12} - T_{21} + \Gamma_{3}.]It must equal zero and the relation given above is thus recovered. Voigt's notation – interpretation of the components of the stress tensor

| top | pdf | Voigt's notation, reduced form of the stress tensor

| top | pdf |

We shall use frequently the notation due to Voigt (1910)[link] in order to express the components of the stress tensor: [\matrix{T_{1} = T_{11}\semi\hfill&T_{2} = T_{22}\semi\hfill&T_{3} = T_{33}\semi\hfill\cr T_{4} = T_{23} = T_{32}\semi\hfill&T_{5} = T_{31} = T_{13}\semi\hfill&T_{6} = T_{12} = T_{21}.\hfill\cr} ]It should be noted that the conventions are different for the Voigt matrices associated with the stress tensor and with the strain tensor (Section[link]).

The Voigt matrix associated with the stress tensor is therefore of the form [\pmatrix{T _{1} &T_{6} &T _{5}\cr T_{6} &T_{2} &T_{4}\cr T_{5} &T_{4} &T_{3}\cr}. ] Interpretation of the components of the stress tensor – special forms of the stress tensor

| top | pdf |

  • (i) Uniaxial stress: let us consider a solid shaped like a parallelepiped whose faces are normal to three orthonormal axes (Fig.[link]). The terms of the main diagonal of the stress tensor correspond to uniaxial stresses on these faces. If there is a single uniaxial stress, the tensor is of the form [\pmatrix{0 &0 &0\cr 0 &0 &0\cr 0 &0 &T_{3}\cr}. ]


    Figure | top | pdf |

    Special forms of the stress tensor. (a) Uniaxial stress: the stress tensor has only one component, [T_{33}]; (b) pure shear stress: [T_{22} = - T_{11}]; (c) simple shear stress: [T_{21} = T_{12}].

    The solid is submitted to two equal and opposite forces, [T_{33} S_{3}] and [- T_{33} S_{3}], where [S_{3}] is the area of the face of the parallelepiped that is normal to the [Ox_{3}] axis (Fig.[link]). The convention used in general is that there is a uniaxial compression if [T_{3}\leq 0] and a uniaxial traction if [T_{3}\geq 0], but the opposite sign convention is sometimes used, for instance in applications such as piezoelectricity or photoelasticity.

  • (ii) Pure shear stress: the tensor reduces to two equal uniaxial constraints of opposite signs (Fig.[link]): [\pmatrix{T_{1} &0 &0\cr 0 &-T_{1} &0\cr 0 &0 &0\cr}. ]

  • (iii) Hydrostatic pressure: the tensor reduces to three equal uniaxial stresses of the same sign (it is spherical): [\pmatrix{-p &0 &0\cr 0 &-p &0\cr 0 &0 &-p\cr},]where p is a positive scalar.

  • (iv) Simple shear stress: the tensor reduces to two equal nondiagonal terms (Fig.[link]c), for instance [T_{12} = T_{21} = T_{6}]. [T_{12}] represents the component parallel to [Ox_{2}] of the stress applied to face 1 and [T_{21}] represents the component parallel to [Ox_{1}] of the stress applied to face 2. These two stresses generate opposite couples that compensate each other. It is important to note that it is impossible to have one nondiagonal term only: its effect would be a couple of rotation of the solid and not a deformation. Boundary conditions

| top | pdf |

If the surface of the solid C is free from all exterior action and is in equilibrium, the stress field [T_{ij}] inside C is zero at the surface. If C is subjected from the outside to a distribution of stresses [T_{n}] (apart from the volume forces mentioned earlier), the stress field inside the solid is such that at each point of the surface [T_{nj} = T_{ij}\alpha_{i}, ]where the [\alpha_{j}]'s are the direction cosines of the normal to the surface at the point under consideration. Local properties of the stress tensor

| top | pdf |

  • (i) Normal stress and shearing stress: let us consider a surface area element dσ within the solid, the normal n to this element and the stress [{\bf T}_{n} ] that is applied to it (Fig.[link]).


    Figure | top | pdf |

    Normal ([\boldnu]) and shearing ([\boldtau]) stress.

    The normal stress, [\boldnu], is, by definition, the component of [{\bf T}_{n}] on n, [\boldnu = {\bf n} ({\bf T}_{n}\cdot {\bf n})]and the shearing stress, [\boldtau], is the projection of [{\bf T}_{n}] on the surface area element,[\boldtau = {\bf n} \wedge ({\bf T}_{n} \wedge {\bf n}) = {\bf T}_{n} - \boldnu. ]

  • (ii) The stress quadric: let us consider the bilinear form attached to the stress tensor: [f({\bf y}) = T_{ij}y_{i}y_{j}. ]The quadric represented by [f({\bf y}) = \varepsilon ]is called the stress quadric, where [\varepsilon = \pm 1 ]. It may be an ellipsoid or a hyperboloid. Referred to the principal axes, and using Voigt's notation, its equation is [y_{i}^2T_{i} = \varepsilon.]To every direction n of the medium, let us associate the radius vector y of the quadric (Fig.[link]) through the relation [{\bf n} = k {\bf y}. ]The stress applied to a small surface element dσ normal to n, [{\bf T}_{n}], is [{\bf T}_{n} = k \boldnabla (f) ]and the normal stress, ν, is [\nu = \alpha_{i}T_{i} = 1/y^{2}, ]where the [\alpha_{i}]'s are the direction cosines of n.


    Figure | top | pdf |

    The stress quadric: application to the determination of the stress applied to a surface element. The surface of the medium is shaded in light grey and a small surface element, dσ, is shaded in medium grey. The stress at P is proportional to [\boldnabla(f)] at the intersection of OP with the stress quadric.

  • (iii) Principal normal stresses: the stress tensor is symmetrical and has therefore real eigenvectors. If we represent the tensor with reference to a system of axes parallel to its eigenvectors, it is put in the form [\pmatrix{T_{1} &0 &0\cr 0 &T_{2} &0\cr 0 &0 &T_{3}\cr}. ][T_{1}], [T_{2}] and [T_{3}] are the principal normal stresses. The mean normal stress, T, is defined by the relation [T = (T_{1} + T_{2} + T_{3})/3 ]and is an invariant of the stress tensor. Energy density in a deformed medium

| top | pdf |

Consider a medium that is subjected to a stress field [T_{ij} ]. It has sustained a deformation indicated by the deformation tensor S. During this deformation, the forces of contact have performed work and the medium has accumulated a certain elastic energy W. The knowledge of the energy density thus acquired is useful for studying the properties of the elastic constants. Let the medium deform from the deformation [S_{ij}] to the deformation [S_{ij} + \delta S_{ij}] under the influence of the stress field and let us evaluate the work of each component of the effort. Consider a small elementary rectangular parallelepiped of sides [2 \Delta x_{1}], [2 \Delta x_{2} ], [2 \Delta x_{3}] (Fig.[link]). We shall limit our calculation to the components [T_{11}] and [T_{12}], which are applied to the faces 1 and 1′, respectively.


Figure | top | pdf |

Determination of the energy density in a deformed medium. PP′ represents the displacement of the small parallelepiped during the deformation. The thick arrows represent the forces applied to the faces 1 and 1′.

In the deformation [\delta S], the point P goes to the point P′, defined by [{\bf PP}' = {\bf u}({\bf r}). ]A neighbouring point Q goes to Q′ such that (Fig.[link]) [{\bf PQ} = \Delta {\bf r} ;\quad {\bf P}'{\bf Q}' = \delta {\bf r}'. ]The coordinates of [\delta {\bf r}'] are given by [\delta x'_{i} = \delta \Delta x_{i} + \delta S_{ij}\delta x_{j}. ]

Of sole importance is the relative displacement of Q with respect to P and the displacement that must be taken into account in calculating the forces applied at Q. The coordinates of the relative displacement are [\delta x'_{i} - \delta \Delta x_{i} = \delta S_{ij}\delta x_{j}. ]

We shall take as the position of Q the point of application of the forces at face 1, i.e. its centre with coordinates [\Delta x_{1}, 0, 0] (Fig.[link]). The area of face 1 is [4 \Delta x_{2} \Delta x_{3}] and the forces arising from the stresses [T_{11}] and [T_{12}] are equal to [4 \Delta x_{2} \Delta x_{3} T_{11}] and [4\Delta x_{2} \Delta x_{3} T_{12}], respectively. The relative displacement of Q parallel to the line of action of [T_{11}] is [\Delta x_{1}\delta S_{11}] and the corresponding displacement along the line of action of [T_{12}] is [\Delta x_{1}\delta S_{21}]. The work of the corresponding forces is therefore [\eqalignno{ \hbox{for} \,\, T_{11}: \quad &4 \Delta x_{1} \Delta x_{2} \Delta x_{3} T_{11} \delta S_{11}\cr \hbox{for} \,\, T_{12}: \quad &4 \Delta x_{1} \Delta x_{2} \Delta x_{3} T_{11} \delta S_{21}.\cr} ]

The work of the forces applied to the face 1′ is the same ([T_{11}], [T_{12}] and [x_{1}] change sign simultaneously). The works corresponding to the faces 1 and 1′ are thus [T_{11} \delta S_{11} \Delta \tau] and [T_{12} \delta S_{21} \Delta \tau] for the two stresses, respectively. One finds an analogous result for each of the other components of the stress tensor and the total work per unit volume is [\delta W = T_{ij} \delta S_{ji}. \eqno(]

1.3.3. Linear elasticity

| top | pdf | Hooke's law

| top | pdf |

Let us consider a metallic bar of length [l_{o}] loaded in pure tension (Fig.[link]). Under the action of the uniaxial stress [T = F/{\cal A}] (F applied force, [{\cal A}] area of the section of the bar), the bar elongates and its length becomes [l = l_{o} + \Delta l]. Fig.[link] relates the variations of Δl and of the applied stress T. The curve representing the traction is very schematic and does not correspond to any real case. The following result, however, is common to all concrete situations:

  • (i) If [0 \,\,\lt\,\, T \,\,\lt\,\, T_{o} ], the deformation curve is reversible, i.e. if one releases the applied stress the bar resumes its original form. To a first approximation, the curve is linear, so that one can write Hooke's law: [{\Delta l \over l} = {1 \over E} T, \eqno(]where E is the elastic stiffness, also called Young's modulus. The physical mechanism at the origin of elasticity is the deformation of the chemical bonds between atoms, ions or molecules in the solid, which act as so many small springs. The reaction of these springs to an applied stress is actually anharmonic and Hooke's law is only an approximation: a Taylor expansion up to the first term. A rigorous treatment of elasticity requires nonlinear phenomena to be taken into account. This is done in Section 1.3.6[link]. The stress below which the strain is recoverable when the stress is removed, [T_{o}], is called the elastic limit.


    Figure | top | pdf |

    Bar loaded in pure tension.


    Figure | top | pdf |

    Schematic stress–strain curve. T: stress; [T_{o}]: elastic limit; [\Delta l /l]: elongation; the asterisk symbolizes the rupture.

  • (ii) If [T \,\,\gt\,\, T_{o}], the deformation curve is no longer reversible. If one releases the applied stress, the bar assumes a permanent deformation. One says that it has undergone a plastic deformation. The region of the deformation is ultimately limited by rupture (symbolized by an asterisk on Fig.[link]). The plastic deformation is due to the formation and to the movement of lattice defects such as dislocations. The material in its initial state, before the application of a stress, is not free in general from defects and it possesses a complicated history of deformations. The proportionality constant between stresses and deformations in the elastic region depends on the interatomic force constants and is an intrinsic property, very little affected by the presence of defects. By contrast, the limit, [T_{o}], of the elastic region depends to a large extent on the defects in the material and on its history. It is an extrinsic property. For example, the introduction of carbon into iron modifies considerably the extent of the elastic region.

    The extents of the elastic and plastic regions vary appreciably from one material to another. Fragile materials, for instance, have a much reduced plastic region, with a clear break. Elastic constants

| top | pdf | Definition

| top | pdf |

Young's modulus is not sufficient to describe the deformation of the bar: its diameter is reduced, in effect, during the elongation. One other coefficient, at least, is therefore necessary. In a general way, let us consider the deformation of a continuous anisotropic medium under the action of a field of applied stresses. We will generalize Hooke's law by writing that at each point there is a linear relation between the components [T_{ij}] of the stress tensor and the components [S_{ij}] of the strain tensor: [\eqalign{S_{ij} &= s_{ijkl}T_{kl}\cr T_{ij} &= c_{ijkl}S_{kl}.\cr} \eqno( ]The quantities [s_{ijkl}] and [c_{ijkl}] are characteristic of the elastic properties of the medium if it is homogeneous and are independent of the point under consideration. Their tensorial nature can be shown using the demonstration illustrated in Section[link] . Let us take the contracted product of the two sides of each of the two equations of ([link] by the components [x_{i}] and [y_{j}] of any two vectors, x and y: [\eqalignno{S_{ij}x_{i}y{_j} &= s_{ijkl}T_{kl}x_{i}y_j\cr T_{ij}x_{i}y{_j} &= c_{ijkl}S_{kl}x_{i}y _{j}.\cr} ]The left-hand sides are bilinear forms since [S_{ij}] and [T_{ij}] are second-rank tensors and the right-hand sides are quadrilinear forms, which shows that [s_{ijkl}] and [c_{ijkl}] are the components of fourth-rank tensors, the tensor of elastic compliances (or moduli) and the tensor of elastic stiffnesses (or coefficients), respectively. The number of their components is equal to 81.

Equations ([link] are Taylor expansions limited to the first term. The higher terms involve sixth-rank tensors, [s_{ijklmn}] and [c_{ijklmn}], with [3^{6} = 729 ] coefficients, called third-order elastic compliances and stiffnesses and eighth-rank tensors with [3^{8} = 6561] coefficients, called fourth-order elastic compliances and stiffnesses. They will be defined in Section[link]. Tables for third-order elastic constants are given in Fumi (1951[link], 1952[link], 1987[link]). The accompanying software to this volume enables these tables to be derived for any point group. Matrix notation – reduction of the number of independent components

| top | pdf |

It is convenient to write the relations ([link] in matrix form by associating with the stress and strain tensors column matrices [1\times 9] and with the tensors of the elastic stiffnesses, c, and of the elastic compliances, s, square matrices [9\times 9] (Section[link] ); these two [9\times 9] matrices are inverse to one another. The number of independent components of the fourth-rank elastic tensors can be reduced by three types of consideration:

  • (i) Intrinsic symmetry: it was shown in Section[link] that tensors representing principal properties are symmetric. This is the case of the elastic tensors and can be shown directly using expression ([link] of the energy stored per unit volume in the medium when we allow it to deform from the state [S_{ij}] to the state [S_{ij} + \delta S_{ij}] under the action of the stress [T_{ij}]: [\delta W = T_{ij} \delta S_{ji}. ]

    Applying relation ([link], we get [\partial W / \partial S_{ij} = c_{ijkl}S_{kl}. \eqno(]Hence, one has by further differentiation [\partial ^{2}W/ (\partial S_{ij}\partial S_{kl}) = c_{ijkl}. ]

    Nothing is changed by interchanging the role of the pairs of dummy indices ij and lk: [\partial ^{2}W/ (\partial S_{kl}\partial S_{ij}) = c_{klij}. ]

    Since the energy is a state function with a perfect differential, one can interchange the order of the differentiations: the members on the left-hand sides of these two equations are therefore equal; one then deduces [c_{ijkl}= c_{klij}. \eqno(]

    The tensor of elastic stiffnesses (and also the tensor of elastic compliances) is thus symmetrical. As shown in Section[link] , the number of their independent components is therefore reduced to 45.

  • (ii) Symmetry of the strain and stress tensors: the strain tensor [S_{ij}] is symmetric by definition (Section[link]) because rotations are not taken into account and the stress tensor [T_{ij}] is symmetric (Section[link]) because body torques are neglected. For this reason, summation ([link], [S_{ij} = s_{ijkl}T_{kl}], can be factorized [equation ([link] ]: [S_{ij} = \textstyle\sum\limits_{l}\displaystyle s_{ijll}T_{ll} + \textstyle\sum\limits_{k \ne l}\displaystyle \left(s_{ijkl} + s_{ijlk}\right)T_{kl}. ]This shows that the number of independent components of tensor [s_{ijkl}] is reduced. This effect of the symmetry of the strain and stress tensors was discussed systematically in Section[link] . It was shown that [\eqalign{s_{ijkl} &= s_{ijlk} = s_{jikl} = s_{jilk}\cr c_{ijkl} &= c_{ijlk} = c_{jikl} = c_{jilk}\cr} \eqno( ]and that the number of independent elastic compliances or stiffnesses is reduced to 21. They are replaced by two-index coefficients constituting [6\times 6] matrices according to Voigt's notation: [\matrix {& \cr i = j\semi\ \ k = l\cr i = j\semi\ \ k \ne l\cr i \ne j\semi\ \ k = l\cr i \ne j\semi\ \ k \ne l\cr }\left|\matrix{{\rm compliances}\hfill\cr s_{\alpha \beta} = s{iil l}\hfill\cr s_{\alpha \beta} = s{iil l} + s{iil k}\hfill\cr s_{\alpha \beta} = s{ijl l} + s{jil l}\hfill\cr s_{\alpha \beta} = s{ijk l} + s{jik l} + s{ijl k} + s{jil k}\hfill\cr}\right\} \matrix{{\rm stiffnesses}\cr\cr c_{\alpha \beta} = cijkl.\cr\cr\cr} \eqno( ]Using these notations and Voigt's notations for the strain and stress tensors, equations ([link] become: [\eqalign{S_{\alpha \beta} &= s_{\alpha \beta}T_\beta\cr T_{\alpha \beta} &= c_{\alpha \beta}S_{\beta}.\cr} \eqno( ]For instance, the first of these two relations is written, once developed, as[\pmatrix{S _{1}\cr S_{2}\cr S_{3}\cr S_{4}\cr S_{5}\cr S_{6}\cr} = \pmatrix{s_{11} &s_{12} &s_{13} &s_{14} &s_{15} &s_{16}\cr s_{12} &s_{22} &s_{23} &s_{24} &s_{25} &s_{26}\cr s_{13} &s_{23} &s_{33} &s_{34} &s_{35} &s_{36}\cr s_{14} &s_{24} &s_{34} &s_{44} &s_{45} &s_{46}\cr s_{15} &s_{25} &s_{35} &s_{45} &s_{55} &s_{56}\cr s_{16} &s_{26} &s_{36} &s_{46} &s_{56} &s_{66}\cr} \pmatrix{T_{1}\cr T_{2}\cr T_{3}\cr T_{4}\cr T_{5}\cr T_{6}\cr}. \eqno( ]

    Matrices [s_{\alpha \beta}] and [c_{\alpha \beta}] are the inverse of each another. It is important to note that they are matrices and not tensors. One cannot apply to them the usual rules of transformation under a change of base since they are only valid for the components of a tensor. In cases where the coordinate system is changed, it is necessary to use the components [s_{ijkl}] and [c_{ijkl}] or to establish the rules of transformation for [s_{\alpha \beta}] and [c_{\alpha \beta}].

  • (iii) Symmetry of the crystal: the reduction of the number of independent components of the matrices [s_{\alpha \beta} ] and [c_{\alpha \beta}] was discussed in Section[link] . As a summary, Table[link] gives the number of independent components for each Laue class.

    Table| top | pdf |
    Number of independent components of the elastic compliances and stiffnesses for each Laue class

    Laue classNo. of independent components
    [\bar{1}, 1] 21
    [2/m, 2, m] 13
    [mmm, 222, 2mm] 9
    [\bar{3}, 3] 7
    [\bar{3}m, 32, 3m] 6
    [4/m, 4,\bar{4}] 7
    [4/mmm, 422, \bar{4}2m, 42m ] 6
    [6/m, 6, \bar{6}] 5
    [6/mmm, 622, \bar{6}2m, 62m ] 5
    [m\bar{3}, 23] 3
    [m\bar{3}m, 432, \bar{4}32 ] 3
    [\infty A_{\infty}/M C, \infty A_{\infty} ] 2
  • (iv) Cauchy relations: in 1828, the French physicist A.-L. Cauchy expressed the strain energy of a discrete-particle structure when its potential energy depends only on attractive and repulsive cohesive forces between pairs of material points (the so-called central or point-to-point force law). These forces act along the line joining the points and depend on the distances between these points. It can be shown that, if in addition every atom is situated at a centre of symmetry, the elastic constants are totally symmetric in their four indices. Thus,[c_{ijkl} = c_{ikjl} = c_{iljk}. ]This additional symmetry to the symmetries required by the previous considerations implies the following relations between the two-index components:[\matrix{c_{23} = c_{44}\semi \quad\hfil& c_{14} = c_{56} \hfil\cr c_{13} = c_{55}\semi\quad\hfil&c_{25} = c_{46}\hfil\cr c_{12} = c_{66}\semi \quad\hfil& c_{36} = c_{45}.\hfil\cr} ]These are known as the Cauchy relations. The further symmetry that they imply reduces the maximum number of independent constants from 21 to 15. In the spherical system (Section[link] ), this number is reduced to 1 ([c_{11} = 3 c_{12}]). Poisson's ratio, [\nu = - s_{12}/s_{11}], defined in Section[link], is then equal to 0.25, taking into account the relations between [s_{11},s_{12}] and [c_{11}, c_{12}] given in Section[link]. In practice, Table[link] shows that in real crystals [\nu] is never equal to 0.25, and that Cauchy relations do not apply.

    Table| top | pdf |
    Elastic compliances of some cubic materials in (GPa)−1 (after Every & McCurdy, 1992[link])

    Material[s_{11}][s_{12}][s_{44}]ν[s_{44}/2(s_{11}-s_{12}) ]
    Ag 22.9 −9.8 22.1 0.428 0.338
    Al 16.0 −5.8 35.3 0.362 0.810
    C (diamond) 1.12 −0.14 1.83 0.125 0.726
    Cu 15.0 −6.3 13.3 0.42 0.312
    Fe 7.67 −2.83 8.57 0.369 0.408
    Ge 9.73 −2.64 14.9 0.271 0.602
    Mo 2.63 −0.74 9.20 0.273 1.304
    Ni 7.67 −2.93 8.23 0.382 0.388
    Pb 93.7 −43.04 68.0 0.459 0.249
    Si 7.73 −2.15 12.20 0.279 0.636
    W 2.45 −0.69 6.24 0.281 0.995
    LiF 11.6 −3.35 15.8 0.289 0.528
    MgO 4.01 −0.96 6.47 0.239 0.651
    NaCl 22.9 −4.8 78.3 0.133 1.452
    GaAs 11.75 −3.66 16.8 0.311 0.545
    ZnS 19.7 −7.6 22.6 0.386 0.414
    BaTiO3 8.33 −2.68 9.24 0.322 0.420
    Adamantane 240 −79 295 0.329 0.462
    Almandine (Fe3Al2Si3O12) 4.036 −1.093 10.77 0.271 1.050
    Spinel (MgAl2O4) 5.80 −2.05 6.49 0.353 0.413

    The extent to which Cauchy relations are fulfilled is often used to assess the validity of the assumption of a central-force law. It is important to apply such a test only to structures having the necessary symmetry properties for Cauchy relations to hold. For instance, f.c.c. and b.c.c. metals have the required symmetry (an inversion centre at each atomic site), while diamond, silicon and germanium do not. Consequently, any apparent fulfilment of the Cauchy relation [c_{12}= c_{66}], or [\nu \simeq 0.25], for a diamond-type structure offers no ground for conclusions about the nature of the force field. Passage from elastic compliances [s_{\alpha \beta}] to elastic stiffnesses [c_{\alpha \beta}]

| top | pdf |

We have noted already that the matrix [c_{\alpha \beta}] is the inverse of the matrix [s_{\alpha \beta}]. These matrices can be written for cubic and isotropic materials as follows: [\eqalignno{s_{\alpha \beta} &= \pmatrix{s_{11} &s_{12} &s_{12} &0 &0 &0\cr s_{12} &s_{11} &s_{12} &0 &0 &0\cr s_{12} &s_{12} &s_{11} &0 &0 &0\cr 0 &0 &0 &s_{44} &0 &0\cr 0 &0 &0 &0 &s_{44} &0\cr 0 &0 &0 &0 &0 &s_{44}\cr} &\cr c_{\alpha \beta} &= \pmatrix{c _{11} &c_{12} &c_{12} &0 &0 &0\cr c_{12} &c_{11} &c_{12} &0 &0 &0\cr c_{12} &c_{12} &c _{11} &0 &0 &0\cr 0 &0 &0 &c_{44} &0 &0\cr 0 &0 &0 &0 &c_{44} &0\cr 0 &0 &0 &0 &0 &c_{44}\cr}, &\cr} ]where we have, for isotropic materials, [\eqalign{s_{44} &= 2(s_{11} - s_{12})\cr c_{44} &= \textstyle{1 \over 2} (c_{11} - c_{12}).\cr} \eqno( ]We easily find that [\eqalign{s_{11} &= {c_{11} + c_{12} \over (c_{11} - c_{12})(c_{11} + 2 c_{12})}\semi \quad c_{11} = {s_{11} + s_{12} \over (s_{11} - s_{12})(s_{11} + 2 s_{12})}\semi \cr s_{12} &= {- c_{12} \over (c_{11} - c_{12})(c_{11} + 2 c_{12})}\semi \quad c_{12} = {- s_{12}\over (s_{11} - s_{12})(s_{11} + 2 s_{12})}\semi \cr s_{44} &= {1 \over c_{44}}\semi \phantom{= {- c_{12} \over (c_{11} - c_{12})(c_{11} c_{12})}\semi} c_{44}= {1 \over s_{44}}. \cr} ]The coefficient [c_{44}] is sometimes called the rigidity modulus. Elastic strain energy

| top | pdf |

Expression ([link] of the strain energy stored per unit volume in a medium for a small deformation can be integrated when the medium is strained under a stress [T_{ij}] according to linear elasticity. Applying relation ([link], one gets for the density of strain energy[W = \textstyle{1 \over 2} T_{ij}S_{ij} = \textstyle{1 \over 2} c_{ijkl}S_{ij}S_{kl}. \eqno( ] Particular elastic constants

| top | pdf | Volume compressibility

| top | pdf |

Let us apply a hydrostatic pressure (Section[link]). The medium undergoes a relative variation of volume [\Delta V/V = ] [S_{1} + S_{2} + S_{3}] (the cubic dilatation, Section[link]). If one replaces in ([link] the stress distribution by a hydrostatic pressure, one obtains for the components of the strain tensor[\eqalignno{S_{1} &= - p (s_{11} + s_{12} + s_{13}) &\cr S_{2} &= - p (s_{12} + s_{22} + s_{23}) &\cr S_{3} &= - p (s_{13} + s_{23} + s_{13}). &\cr} ]From this, we deduce the volume compressibility, [\chi], which is the inverse of the bulk modulus, κ: [\chi =\kappa^{-1} = - {1\over p} {\delta V \over V} = s_{11} + s_{22} + s_{33} + 2(s_{12} + s_{23} + s_{13}). \eqno( ]This expression reduces for a cubic or isotropic medium to [\chi = \kappa^{-1} = 3(s_{11} + 2s_{12}). \eqno(] Linear compressibility

| top | pdf |

Under the action of a hydrostatic pressure, each vector assumes a different elongation. This elongation is given by equation ([link]: [{u_{r} \over r} = {S_{ij} x_{i}x_{j} \over r^{2}} = S_{ij}\alpha_{i}\alpha_{j} = s_{ijkl}T_{kl}\alpha_{i}\alpha_{j}, ]where the [\alpha_{i}]'s are the direction cosines of r. The coefficient of linear compressibility is, by definition, [(-1/p ) ( u_{r}/r)]. Replacing [T_{kl}] by its value [-p \delta_{kl}], we obtain for the coefficient of linear compressibility [-{1 \over p} {u_{r} \over r} = s_{ijkl}\alpha_{i}\alpha_{j}. ]In the case of a cubic or isotropic medium, this expression reduces to [-{1 \over p} {u_{r} \over r} = s_{11} + 2s_{12}.]

The coefficient of linear compressibility is then equal to one third of the coefficient of volume compressibility. We note that the quadric of elongations is a sphere. Young's modulus, Poisson's ratio

| top | pdf |

If the applied stress reduces to a uniaxial stress, [T_{11} = T], the strain tensor is of the form [S_{\alpha} = s_{1\alpha} T. ]In particular, [S_{1} = s_{11}T\semi \quad S_{2} = s_{12}T. ]We deduce from this that Young's modulus (equation[link] is[E = 1/s_{11}. \eqno(]

The elongation of a bar under the action of a uniaxial stress is characterized by [S_{1}] and the diminution of the cross section is characterized by [S_{2}] and [S_{3}]. For a cubic material, the relative diminution of the diameter is [S_{2} = S_{3} = s_{12}T. ]One deduces from this that [s_{12}] is necessarily of opposite sign to [s_{11}] and one calls the ratio [\nu = -s_{12}/s_{11}] Poisson's ratio.

Putting this value into expression ([link] for the coefficient of compressibility in cubic or isotropic materials gives [\kappa^{-1} = 3 s_{11}(1 - 2 \nu). \eqno(]

As the coefficient of compressibility, by definition, is always positive, we have [0 \,\,\lt\,\, \nu \,\,\lt\,\, 0.5.]

In practice, Poisson's ratio is always close to 0.3. It is a dimensionless number. The quantity [s_{44}/2 (s_{11} - s_{12})] represents the departure from isotropy of the material and is the anisotropy factor. It is to be noted that cubic materials are not isotropic for elastic properties. Table[link] gives the values of [s_{11}], [s_{12}], [s_{44}], ν and [s_{44}/2 (s_{11} - s_{12})] for a few cubic materials. Variation of Young's modulus with orientation

| top | pdf |

It is interesting to calculate Young's modulus in any direction. For this it is sufficient to change the axes of the tensor [s_{ijkl} ]. If A is the matrix associated with the change of axes, leading to the direction [x_{1}] changing to the direction [x'_{1}], then Young's modulus in this new direction is [E' = 1/s'_{11}]with [s'_{11}= s'_{1111} = A_{1i}A_{1j}A_{1k}A_{1l}s_{ijkl}. \eqno( ]The matrix coefficients [A_{1i}] are the direction cosines of [Ox'_{1}] with respect to the axes [Ox_{1}], [Ox_{2}] and [Ox_{3}]. In spherical coordinates, they are given by (Fig.[link]) [A_{11}= \cos \theta \sin \varphi, \quad A_{12}= \sin \theta \sin \varphi, \quad A_{13}= \cos \varphi, ]where θ is the angle between [Ox'_{1}] and [Ox_{1}], and ϕ is the angle between [Ox'_{1} ] and [Ox_{3}]. Using the reduction of [s_{ijkl}] for the various crystal classes (Section[link] ), we find, in terms of the reduced two-index components, the following.


Figure | top | pdf |

Spherical coordinates.

  • (i) Triclinic system (groups 1, [\bar{1}]): [\eqalign{s'_{11} &= \left[s_{11}\cos^{4}\theta + s_{22}\sin^{4}\theta + (2s_{12}+ s_{66})\sin^{2}2\theta /4 \right.\cr &\quad + (s_{16}\cos \theta + s_{26}\sin \theta)\sin 2 \theta]\sin^{4}\varphi\cr &\quad + 2 \big\{\left[(s_{25} + s_{46})\sin \theta + (s_{14}+ s_{56})\cos \theta \right]\sin2 \theta /2 \cr &\quad+ s_{15} \cos^3\theta + s_{24}\sin^{3}\theta \big\}\cos \varphi \sin^{3}\varphi \cr &\quad + \left[(2s_{23}+ s_{44})\sin^{2}\theta + (2s_{13}+ s_{55})\cos^{2}\theta \right.\cr &\quad+ (s_{36}+ s_{45})\sin2 \theta \bigr]\sin^{2}2\varphi /4\cr &\quad + 2 (s_{35}\cos \theta +s_{34}\sin \theta) \cos^{3}\varphi \sin \varphi + s_{33}\cos^{4}\varphi.\cr} ]

  • (ii) Monoclinic system (groups 2, m, [2/m]): [\eqalign{s'_{11} &= \left[s_{11}\cos^{4}\theta + s_{22}\sin^{4}\theta + (2s_{12}+ s_{66})\sin^{2}2\theta /4\right]\sin^{4}\varphi \cr &\quad + 2 \left[(s_{25} + s_{46})\sin^{2}\theta + s_{15}\cos^{2}\theta \right] \cos \varphi \sin^{3}\varphi \cos \theta \cr &\quad + \left[(2s_{23}+ s_{44})\sin^{2}\theta + (2s_{13}+ s_{55})\cos^{2}\theta \right]\sin^{2}2\varphi /4\cr &\quad + 2 s_{35}\cos^{3}\varphi \sin \varphi \cos \theta + s_{33}\cos^{4}\varphi.\cr} ]

  • (iii) Orthorhombic system (groups 222, 2mm, mmm): [\eqalign{s'_{11} &= \left[s_{11}\cos^{4}\theta + s_{22}\sin^{4}\theta + (2s_{12} + s_{66})\sin^{2}2\theta /4\right]\sin^{4}\varphi\cr &\quad + \left[(2s_{23}+ s_{44})\sin^{2}\theta + (2s_{13}+ s_{55})\cos^{2}\theta \right] \sin^{2}2\varphi /4 \cr&\quad+ s_{33}\cos^{4}\varphi.\cr} ]

  • (iv) Trigonal system (groups 3, [\bar{3}]): [\eqalign{s'_{11} &= s_{11}\sin^{4}\varphi + s_{33}\cos^{4}\varphi + (2s_{13} + s_{44})\sin^{2}2\varphi /4 \cr&\quad+ (s_{14}\sin 3\theta - s_{25}\cos 3\theta)\sin 2\varphi \sin^{2}\varphi.} ]

  • (v) Trigonal system (groups 32, 3m, [\bar{3}m]): [\eqalign{s'_{11}&= s_{11}\sin^{4}\varphi + s_{33}\cos^{4}\varphi + (2s_{13}+ s_{44})\sin^{2}2\varphi /4 \cr&\quad+ s_{14}\sin3\theta \sin2\varphi \sin^{2}\varphi.} ]

  • (vi) Tetragonal system (groups 4, [\bar{4}, 4/m]): [\eqalign{s'_{11} &= \left\{s_{11}+ \left[s_{66}- 2(s_{11}- s_{12})\right] \sin^{2}2\theta /4\right\}\sin^{4}\varphi + s_{33}\cos^{4}\varphi\cr &\quad +(2s_{13} + s_{44})\sin^{2} 2\varphi /4 + s_{16} \sin 4\theta \sin^{4} \varphi /2.\cr} ]

  • (vii) Tetragonal system (groups 422, [\bar{4}2m], 4mm, [4/mmm]): [\eqalign{s'_{11}&= \left\{s_{11} + \left[s_{66}- 2(s_{11}- s_{12})\right] \sin^{2}2\theta /4\right\}\sin^{4}\varphi + s_{33}\cos^{4}\varphi \cr&\quad+ (2s_{13}+ s_{44})\sin^{2}2\varphi /4.} ]

  • (viii) Hexagonal system: [s'_{11}= s_{11}\sin^{4}\varphi + s_{33}\cos^{4}\varphi + (2s_{13}+ s_{44})\sin^{2}2\varphi /4. ]

  • (ix) Cubic system: [s'_{11}= s_{11}+ \left[s_{44}- 2(s_{11}- s_{12})\right]\sin^{2}\varphi \left[\cos^{2}\varphi + \sin^{2}2\theta \sin^{2}\varphi /4\right]. ]

    This expression reduces to [s_{11}] if [s_{44} - 2(s_{11} - s_{12}) = 0] and we retrieve the relation between elastic compliances in an isotropic material (Sections[link] and[link]).

The representation surface of [s_{11}], the inverse of Young's modulus, is illustrated in Figure[link] for crystals of different symmetries. As predicted by the Neumann principle, the representation surface is invariant with respect to the symmetry elements of the point group of the crystal but, as stated by the Curie laws, its symmetry can be larger. In the examples of Fig.[link], the symmetry of the surface is the same as that of the point group for aluminium (Fig.[link]), tungsten (Fig.[link]) and sodium chloride (Fig.[link]), which have [m\bar{3}m] as point group, for tin (Fig.[link], [4/mmm]) and for calcite (Fig.[link], [\bar{3}m]). But in the case of zinc (Fig.[link], [6/mmm]), the surface is of revolution and has a larger symmetry. It is interesting to compare the differences in shapes of the representation surfaces for the three cubic crystals, depending on the value of the anisotropy factor, which is larger than 1 for aluminium, smaller than 1 for sodium chloride and close to 1 for tungsten (see Table[link]). In this latter case, the crystal is pseudo-isotropic and the surface is practically a sphere.


Figure | top | pdf |

Representation surface of the inverse of Young's modulus. (a) Al, cubic, anisotropy factor [\gt\,\, 1]; (b) W, cubic, anisotropy factor [= 1]; (c) NaCl, cubic, anisotropy factor [\lt\,\, 1]; (d) Zn, hexagonal; (e) Sn, tetragonal; (f) calcite, trigonal. Isotropic materials

| top | pdf |

The isotropy relation between elastic compliances and elastic stiffnesses is given in Section[link]. For reasons of symmetry, the directions of the eigenvectors of the stress and strain tensors are necessarily the same in an isotropic medium. If we take these directions as axes, the two tensors are automatically diagonalized and the second relation ([link] becomes [\eqalign{T_{1} &= c_{11}S_{1}+ c_{12}(S_{2}+ S_{3})\cr T_{2} &= c_{12}S_{1}+ c_{11}S_{2}+ c_{12}S_{3}\cr T_{3} &= c_{12}(S_{1}+ S_{2}) + c_{11}S_{3}.\cr} ]

These relations can equally well be written in the symmetrical form [\eqalign{T_{1}&= (c_{11}- c_{12})S_{1}+ c_{12}(S_{1}+ S_{2}+ S_{3})\cr T_{2}&= (c_{11}- c_{12})S_{2}+ c_{12}(S_{1}+ S_{2}+ S_{3})\cr T_{3}&= (c_{11}- c_{12})S_{3}+ c_{12}(S_{1}+ S_{2}+ S_{3}).\cr} ]

If one introduces the Lamé constants, [\eqalign{\mu & = (1/2) (c_{11}- c_{12}) = c_{44}\cr \lambda &= c_{12},\cr} ]the equations may be written in the form often used in mechanics: [\eqalign{T_{1} &= 2\mu S_{1} + \lambda (S_{1} +S_{2} + S_{3})\cr T_{2} &= 2\mu S_{2} + \lambda (S_{1} +S_{2} + S_{3})\cr T_{3} &= 2\mu S_{3} + \lambda (S_{1} +S_{2} + S_{3}).\cr} \eqno( ]

Two coefficients suffice to define the elastic properties of an isotropic material, [s_{11}] and [s_{12}], [c_{11}] and [c_{12}], μ and λ, μ and ν, etc. Table[link] gives the relations between the more common elastic coefficients.

Table| top | pdf |
Relations between elastic coefficients in isotropic media

CoefficientIn terms of μ and λIn terms of μ and νIn terms of [c_{11} \hbox{ and } c_{12} ]
[c_{11}] [2\mu + \lambda] [2\mu(1-\nu)(1–2\nu)] [c_{11}]
[c_{12}] λ [2\mu \nu(1–2\nu)] [c_{12}]
[c_{44} = 1/s_{44}] μ μ [(c_{11}-c_{12})/2]
[E = 1/s_{11}] [\mu(2\mu+3\lambda)/(\mu+\lambda) ] [2\mu(1+\nu)] See Section[link]
[s_{12}] [-\lambda/\left[2\mu(2\mu+3\lambda)\right] ] [-\nu/\left[2\mu(1+\nu)\right] ] See Section[link]
κ [3/(2\mu+3\lambda)] [3(1–2\nu)/\left[2\mu(1+\nu)\right] ] [3/(c_{11}+2c_{12})]
[\nu = - s_{12}/s_{11}] [\lambda/\left[2(2\mu+3\lambda)\right] ] ν [c_{11}/(c_{11}+c_{12})] Equilibrium conditions of elasticity for isotropic media

| top | pdf |

We saw in Section[link] that the condition of equilibrium is [\partial T_{ij} /\partial x_i + \rho F_{j} =0. ]

If we use the relations of elasticity, equation ([link], this condition can be rewritten as a condition on the components of the strain tensor: [c_{ijkl }{\partial S_{kl} \over \partial x_j} + \rho F_{i}= 0. ]Recalling that [S_{kl }= \textstyle{1 \over 2}\displaystyle \left[ {\partial u_k \over \partial x_l } + {\partial u_l\over \partial x_k }\right], ]the condition becomes a condition on the displacement vector, [{\bf u}({\bf r})]: [c_{ijkl} {\partial^2 \over \partial x_l \partial x_{j}} + \rho F_{i}= 0. ]In an isotropic orthonormal medium, this equation, projected on the axis [Ox_{1}], can be written with the aid of relations ([link] and ([link]: [\eqalign{&c_{11}{\partial^2 u_1 \over (\partial x_1){}^2 } + c_{12}\left[{\partial^2 u_2 \over \partial x_1 \partial x_2 } + {\partial^2 u_3 \over \partial x_1 \partial x_3 }\right] \cr&\quad+ \textstyle{1\over 2}\displaystyle (c_{11}-c_{12})\left[{\partial^2 u_1 \over (\partial x_2){}^2 } + {\partial^2 u_3 \over \partial x_1 \partial x_3 } + {\partial^2 u_1 \over (\partial x_3){}^2 }\right]+ \rho F_i \cr&\quad\quad= 0.} ]This equation can finally be rearranged in one of the three following forms with the aid of Table[link]. [\eqalign{\textstyle{1 \over 2}\displaystyle (c_{11}-c_{12})\Delta {\bf u} + \textstyle{1 \over 2}\displaystyle (c_{11}+c_{12}) \boldnabla(\boldnabla{\bf u}) + \rho {\bf F}&=0\cr \mu \Delta {\bf u} + (\mu + \lambda) \boldnabla(\boldnabla{\bf u}) + \rho {\bf F}&=0\cr \mu \left[\Delta {\bf u} + {1 \over 1 - 2\nu}\boldnabla(\boldnabla{\bf u}) \right]+ \rho {\bf F}&=0.\cr} \eqno( ]

1.3.4. Propagation of elastic waves in continuous media – dynamic elasticity

| top | pdf | Introduction

| top | pdf |

The elastic properties of materials have been considered in the preceding section in the static state and the elastic constants have been defined in terms of the response of the material to particular static forces. It is effectively the way the elastic constants have been measured in the past, although the measurements could not be very precise. A way of proceeding frequently used now is to excite a mechanical wave in the crystal and measure its propagation velocity or the wavelength associated with a particular frequency. One method consists in sending a train of ultrasonic waves through the crystal; one uses a pulse generator and a piezoelectric transducer glued to the crystal. The elapsed time between the emission of the train of waves and its reception after reflection from the rear face of the sample is then measured. Another method involves producing a system of standing waves after reflection at the inner surface of the crystal and determining the set of resonance frequencies. The experimental techniques will be described in Section[link].

The purpose of the next sections is to establish relations between the wavelength – or the velocity of propagation – and the elastic constants. Equation of propagation of a wave in a material

| top | pdf |

Consider the propagation of a wave in a continuous medium. The elongation of each point will be of the form [{\bf u} = {\bf u_0}\exp(2\pi i\nu t) \exp(-2\pi i{\bf q \cdot r}), \eqno( ]where ν is the frequency and q is the wavevector. The velocity of propagation of the wave is [V = \nu / q. \eqno( ]

We saw in Section[link] that the equilibrium condition is [c_{ijkl}{\partial^2 u_k \over \partial x_l \partial x_{j}}+ \rho F_{i} = 0. ]Here the only volume forces that we must consider are the inertial forces: [c_{ijkl}{\partial^2 u_k \over \partial x_l \partial x_{j}} = \rho {\partial^2 x_i\over \partial t^2}. \eqno( ]

The position vector of the point under consideration is of the form [{\bf r} = {\bf r}_0 + {\bf u},]where only [{\bf u}] depends on the time and [{\bf r}_0] defines the mean position. Equation ([link] is written therefore [c_{ijkl }{\partial^2 u_k \over \partial x_l \partial x_{j}} = \rho {\partial^2 u_i\over \partial t^2 }. \eqno( ]

Replacing u by its value in ([link], dividing by [-4 \pi^{2}] and using orthonormal coordinates, we get [c_{ijkl }u_{k}q_{j}q_{l } = \rho \nu ^{2}u_{i}. \eqno( ]

It can be seen that, for a given wavevector, [\rho \nu^{2} ] appears as an eigenvalue of the matrix [c_{ijkl } u_{k}q_{j}q_{l } ] of which the vibration vector u is an eigenvector. This matrix is called the dynamical matrix, or Christoffel matrix. In order that the system ([link] has a solution other than a trivial one, it is necessary that the associated determinant be equal to zero. It is called the Christoffel determinant and it plays a fundamental role in the study of the propagation of elastic waves in crystals.

Let [\alpha_{1}, \alpha_{2}, \alpha_{3}] be the direction cosines of the wavevector q. The components of the wavevector are [q_{i} = q \alpha_{i}.]With this relation and ([link], the system ([link] becomes [c_{ijkl} u_{k}\alpha_j\alpha_l = \rho \nu ^{2}u_{i}. \eqno( ]Putting [\Gamma_{ik} = c_{ijkl}\alpha_{j}\alpha_{l} \eqno( ]in ([link], the condition that the Christoffel determinant is zero can be written [\Delta \left(\Gamma_{ik} - \rho \nu^{2}\delta _{ik}\right) = 0. \eqno( ]

On account of the intrinsic symmetry of the tensor of elastic stiffnesses, the matrix [\Gamma_{ik}] is symmetrical.

If we introduce into expression ([link] the elastic stiffnesses with two indices [equation ([link]], we find, for instance, for [\Gamma_{11}] and [\Gamma_{12}] [\eqalign{\Gamma_{11} &= c_{11}(\alpha_{1}){}^{2}+ c_{66}(\alpha_{2}){}^{2}+ c_{55}(\alpha_{3}){}^{2} + 2c_{16}\alpha_{1}\alpha_{2}\cr&\quad+ 2c_{15}\alpha_{1}\alpha_{3} +2c_{56}\alpha_{2}\alpha_{3} \cr \Gamma_{12}&= c_{16}(\alpha_{1}){}^{2}+ c_{26}(\alpha_{2}){}^{2} + c_{45}(\alpha_{3}){}^{2}+ (c_{12}+c_{66})\alpha_{1}\alpha_{2} \cr&\quad+ (c_{14}+c_{56})\alpha_{1}\alpha_{3}+(c_{46} + c_{25})\alpha_{2}\alpha_{3}.\cr} ]

The expression for the effective value, [c^{e}_{ijkl}], of the `stiffened' elastic stiffness in the case of piezoelectric crystals is given in Section[link] . Dynamic elastic stiffnesses

| top | pdf |

Equation ([link] may be written [\Gamma_{ik} = \textstyle\sum\limits_{j \ne l}\displaystyle \left[c_{ijkl }+ c_{il kj} \right] \alpha_{j}\alpha_{l}. ]

This shows that in a dynamic process only the sums [[c_{ijkl} + c_{ilkj}] ] can be measured and not [c_{ijkl}] and [c_{il kj}] separately. On the contrary, [c_{ijij}] can be measured directly. In the cubic system therefore, for instance, [c_{1122}] is determined from the measurement of [[c_{1122} + c_{1221}] ] on the one hand and from that of [c_{1221}] on the other hand. Polarization of the elastic waves

| top | pdf |

The Christoffel determinant has three roots and the Christoffel matrix, being Hermitian with real coefficients, has three real eigenvalues and three orthogonal eigenvectors. The wavevector q, therefore, encompasses three waves with vibration vectors [{\bf u}_{1} ], [{\bf u}_{2}], [{\bf u}_{3}] which are perpendicular to one another. In the general case, there is no particular angular relationship between the vibration vectors (or polarization vectors). However, if the latter are parallel to certain symmetry directions in the crystal, one of the vibration vectors is along this direction. The corresponding wave is called longitudinal. The two other waves have their polarization direction perpendicular to the wavevector and are thus transverse. If one of the polarization vectors is almost parallel to the wavevector, which often happens, then one speaks of the vibration as being quasi-longitudinal. Relation between velocity of propagation and elastic stiffnesses

| top | pdf |

We shall limit ourselves to cubic, hexagonal and tetragonal crystals and consider particular cases. Cubic crystals

| top | pdf |

  • (i) The wavevector is parallel to [100]. The Christoffel determinant reduces to [\pmatrix{c_{11} - \rho \nu^{2} &0 &0\cr 0 &c_{44} - \rho \nu^{2} &0\cr 0 &0 &c_{44} - \rho \nu^{2}\cr} = 0. ]The three solutions are given in Table[link]. These results are valid for a wave propagating in any direction in an isotropic medium.

    Table| top | pdf |
    Velocity of propagation when the wavevector is parallel to [100] (cubic crystals)

    Velocity of propagationPolarization vectorNature of the wave
    [\nu_\parallel = \sqrt{c_{11}/\rho} ] [100] Longitudinal
    [\nu_{\perp} = \sqrt{c_{44}/\rho} ] [010] Transverse
    [\nu_{\perp} = \sqrt{c_{44}/\rho} ] Any vector normal to [100] Transverse
  • (ii) The wavevector is parallel to [110]. The direction cosines of the wavevector are [1/\sqrt{2}], [1/\sqrt{2}], 0. The Christoffel determinant assumes the form [\pmatrix{{1 \over 2}(c_{11} + c_{44}) - \rho \nu^{2} &{1 \over 2}(c_{12} + c_{44}) &0\cr {1 \over 2}(c_{12} + c_{44}) &{1 \over 2}(c_{11} + c_{44}) - \rho \nu^{2} &0\cr 0 &0 &c_{44} - \rho \nu^{2}\cr} = 0. ]The three solutions are given in Table[link].

    Table| top | pdf |
    Velocity of propagation when the wavevector is parallel to [110] (cubic crystals)

    Velocity of propagationPolarization vectorNature of the wave
    [\nu_{\parallel} = \sqrt{[c_{44} + {1 \over 2}(c_{11} + c_{12})]/\rho} ] [110] Longitudinal
    [\nu_{\perp} = \sqrt{{1 \over 2}(c_{11} - c_{12})/\rho} ] [[1\bar{1}0]] Transverse
    [\nu_{\perp} = \sqrt{c_{44}/\rho} ] [001] Transverse
  • (iii) The wavevector is parallel to [111]. The Christoffel determinant assumes the form [\pmatrix{c_{11} + 2c_{44} - \rho \nu^{2} &c_{12} + c_{44} &c_{12} + c_{44}\cr c_{12} + c_{44} &c_{11} + 2 c_{44}) - \rho \nu^{2} &c_{12} + c_{44}\cr c_{12} + c_{44} &c_{12} + c_{44} &c_{44} - \rho \nu^{2}\cr} = 0. ]The solutions are given in Table[link].

    Table| top | pdf |
    Velocity of propagation when the wavevector is parallel to [111] (cubic crystals)

    Velocity of propagationPolarization vectorNature of the wave
    [\nu_{\parallel} = \sqrt{(c_{11} + 2c_{12} + 4c_{44})/3\rho} ] [111] Longitudinal
    [\nu_{\perp} = \sqrt{(c_{11} - c_{12} + c_{44})/3\rho} ] Any vector normal to [111] Transverse Hexagonal crystals

| top | pdf |

In hexagonal crystals, there are five independent elastic stiffnesses, [c_{11}], [c_{33}], [c_{12} ], [c_{13}], [c_{44}] and [c_{66} = (c_{11} - c_{12})/2] (Section[link] ).

  • (i) The wavevector is parallel to [001]. The Christoffel determinant reduces to [\pmatrix{c_{44} - \rho \nu^{2} &0 &0\cr 0 &c_{44} - \rho \nu^{2} &0\cr 0 &0 &c_{33} - \rho \nu^{2}\cr} = 0. ]The solutions are given in Table[link].

    Table| top | pdf |
    Velocity of propagation when the wavevector is parallel to [001] (hexagonal crystals)

    Velocity of propagationPolarization vectorNature of the wave
    [\nu_{\parallel} = \sqrt{c_{33})/\rho} ] [100] Longitudinal
    [\nu_{\perp} = \sqrt{c_{44}/\rho} ] Any vector normal to [001] Transverse
  • (ii) The wavevector is parallel to [100]. The Christoffel determinant readily reduces to [\pmatrix{c_{11} - \rho \nu^{2} &0 &0\cr 0 &c_{66} - \rho \nu^{2} &0\cr 0 &0 &c_{44} - \rho \nu^{2}&\cr} = 0. ]The three solutions are given in Table[link].

    Table| top | pdf |
    Velocity of propagation when the wavevector is parallel to [100] (hexagonal crystals)

    Velocity of propagationPolarization vectorNature of the wave
    [\nu_{\parallel} = \sqrt{c_{11}/\rho} ] [100] Longitudinal
    [\nu_{\perp} = \sqrt{c_{66}/\rho} ] [010] Transverse
    [\nu_{\perp} = \sqrt{c_{44}/\rho} ] [001] Transverse Tetragonal crystals (classes [4mm], [\bar{4}2m], [4/mmm])

| top | pdf |

In tetragonal crystals, there are six independent elastic stiffnesses, [c_{11}], [c_{33}], [c_{12} ], [c_{13}], [c_{44}] and [c_{66}] (Section[link] ).

  • (i) The wavevector is parallel to [001]. The Christoffel determinant reduces to [\pmatrix{c_{44} - \rho \nu^{2} &0 &0\cr 0 &c_{44} - \rho \nu^{2} &0\cr 0 &0 &c_{33} - \rho \nu^{2}\cr} = 0. ]The three solutions are given in Table[link].

    Table| top | pdf |
    Velocity of propagation when the wavevector is parallel to [001] (tetragonal crystals)

    Velocity of propagationPolarization vectorNature of the wave
    [\nu_{\parallel} = \sqrt{c_{33}/\rho} ] [100] Longitudinal
    [\nu_{\perp} = \sqrt{c_{44}/\rho} ] [010] Transverse
    [\nu_{\perp} = \sqrt{c_{44}/\rho} ] [001] Transverse
  • (ii) The wavevector is parallel to [100]. The Christoffel determinant reduces to [\pmatrix{c_{11} - \rho \nu^{2} &0 &0\cr 0 &c_{66}-^{2}\rho \nu^{2} &0\cr 0 &0 &c_{44}-^{2}\rho \nu^{2}\cr} = 0. ]The three solutions are given in Table[link].

    Table| top | pdf |
    Velocity of propagation when the wavevector is parallel to [100] (tetragonal crystals)

    Velocity of propagationPolarization vectorNature of the wave
    [\nu_{\parallel} = \sqrt{c_{11}/\rho} ] [100] Longitudinal
    [\nu_{\perp} = \sqrt{c_{66}/\rho} ] [010] Transverse
    [\nu_{\perp} = \sqrt{c_{44}/\rho} ] [001] Transverse Experimental determination of elastic constants

| top | pdf | Introduction

| top | pdf |

As mentioned in Section[link], the elastic constants of a material can be obtained by the elastic response of the material to particular static forces; however, such measurements are not precise and the most often used approach nowadays consists of determining the velocity of ultrasonic waves propagating along different directions of the crystal and calculating the elastic constants from the Christoffel determinants ([link]. The experimental values are often accurate enough to justify the distinction between static and dynamic values of the elastic constants and between phase and group velocities, and the careful consideration of the frequency range of the experiments.

  • (i) Static and dynamic elastic constants. When one measures the elastic response of a material to external static forces, work is done and heat is produced. In general, the external forces are applied slowly (quasi-elastic processes) and the solid body remains in thermal equilibrium with its surroundings, which can be considered as a heat reservoir. In this case, the measured elastic constants are isothermal elastic constants as defined in Section[link] . On the contrary, when ultrasonic waves propagate in a solid body, the compressed regions are slightly hotter than the expanded regions; the associated temperature gradients give rise to irreversible processes of thermal conduction. However, in the megahertz or higher-frequency ranges, the processes of heat exchange between different regions or between the solid and its surroundings are slow compared with the period of the ultrasonic wave and hence can be considered as isentropic. In this case, the measured velocities are isentropic and the elastic stiffnesses deduced from them are the adiabatic elastic stiffnesses as defined in Section[link] . The differences between isothermal and adiabatic elastic stiffnesses can be calculated from equation ([link] : [\eqalign{\delta T_{ij} &= (c_{ijkl})^{\Theta} \delta S_{kl} - \lambda_{ij}\delta \Theta\cr \delta \sigma &= \lambda_{kl} \delta S_{kl} + {\rho C^{S} \over \Theta} \delta \Theta. \cr} ]

    By combining these two equations, it is possible to obtain relations between the isothermal elastic stiffnesses, [(c_{ijkl})^{\Theta} ], and the adiabatic elastic stiffnesses, [(c_{ijkl})^{\sigma} ]: [(c_{ijkl})^{\sigma} = (c_{ijkl})^{\Theta} + {\lambda_{ij}\lambda_{kl}\over \rho c^{S}}, ]where [c^{S}] is the specific heat at constant strain.

    This relation shows that adiabatic elastic stiffnesses are larger than isothermal elastic stiffnesses, at least for the terms of the main diagonal of the elastic matrix. In general, the differences between isothermal and adiabatic elastic stiffnesses are less than 1%. An exception to this statement concerns the temperature region near the critical temperature [T_{c}] where a phase transition occurs.

  • (ii) Frequency dependence of the elastic constants. Dynamic measurements of the elastic constants can be performed at a wide range of ultrasonic frequencies. Currently used techniques for investigating a wide range of frequencies are electronic pulse-echo techniques in the megahertz frequency range and Brillouin scattering in the 10 GHz frequency range. To evaluate the possible differences between dynamic elastic constants determined in different frequency ranges, high-accuracy measurements of the absolute value of the elastic constants are required. From optical resonance techniques in the megahertz frequency range and Brillouin scattering in the gigahertz frequency range, it has been shown that the results are the same within 0.1 to 0.2% in the case of simple and stable crystals such as sodium chloride when the ultrasonic waves can be considered as a `passive probe' (Michard et al., 1971[link]). Larger differences are observed when ultrasonic waves interact with the medium, when they favour or impede structural changes of the material or when they are absorbed for a particular frequency.

  • (iii) Phase velocity and group velocity. When the velocity of ultrasonic waves is frequency dependent (dispersive medium), one has to distinguish phase velocity and group velocity. If a resonance technique is used, the resonance frequencies correspond to phase-matching conditions and the calculation of velocity from the resonance frequencies leads to `phase velocity'. If a pulse-echo technique is used, the transit time of an elastic pulse is measured and the interpretation of the measurement is more complex. In this case, Brillouin (1932)[link] has shown that the head and the tail of the pulse travel with the group velocity. Differences of a few per cent can be observed between phase and group velocities in dispersive media. Resonance technique

| top | pdf |

The use of the resonance technique is a well established approach for determining the velocity of sound in a gas by observing nodes and antinodes of a system of standing waves produced in the so-called Kund tube. In the case of transparent solids, optical means allow us to visualize the standing waves and to measure the wavelength directly (Zarembowitch, 1965[link]). An easier procedure can be used: let us consider a transparent crystal in the shape of a parallelepiped (Fig.[link]). A piezoelectric transducer is glued to the crystal and excited at varying frequencies. If the bonding between the transducer and the crystal is loose enough, the crystal can be considered as free from stress and the sequence of its resonance frequencies is given by[\nu = nV/2l,]where n is an integer, V the phase velocity of the wave in the direction orthogonal to the parallel faces and l the distance between these faces.


Figure | top | pdf |

Resonance technique: standing waves excited in a parallelepiped.

The looseness of the bonding can be checked by the regularity of the arithmetic ratio, [V/2l]. On account of the elasto-optic coupling, a phase grating is associated with the elastic standing-wave system and a light beam can be diffracted by this grating. The intensity of the diffraction pattern is maximum when resonance occurs. A large number of resonance frequencies can be detected, usually more than 100, sometimes 1000 for non-attenuating materials. Consequently, in favourable cases the absolute value of the ultrasonic velocity can be determined with an uncertainty less than [10^{-3} ]. Pulse-echo techniques

| top | pdf |

Pulse-echo techniques are valid for transparent and opaque materials. They are currently used for measuring ultrasonic velocities in solids and can be used in very simple as well as in sophisticated versions according to the required precision (McSkimmin, 1964[link]). In the simplest version (Fig.[link]), an electronic pulse generator excites the mechanical vibrations of a piezo-electric transducer glued to one of two plane-parallel faces of a specimen. An ultrasonic pulse whose duration is of the order of a microsecond is generated and transmitted through the specimen. After reflection at the opposite face, it returns and, when it arrives back at the transducer, it gives rise to an electronic signal, or echo. The whole sequence of such echos is displayed on the screen of an oscilloscope and it is possible to measure from them the time interval for transit. Usually, X-cut quartz crystals or ferroelectric ceramics are used to excite longitudinal waves and Y-cut quartz is used to excite transverse waves. In many cases, a circulator, or gate, is used to protect the receiver from saturation following the main `bang'. This method is rough because the beginning and the end of a pulse are not well characterized. Several improvements have therefore been made, mainly based on interferometric techniques (pulse-superposition method, `sing around' method etc.). Nevertheless, if the absolute value of the ultrasonic velocity is not determined with a high accuracy by using pulse-echo techniques, this approach has proved valuable when relative values of ultrasonic velocities are needed, e.g. temperature and pressure dependences of ultrasonic velocities.

  • (i) Pulse-superposition method. A piezoelectric transducer initiates ultrasonic pulses in the specimen. These pulses echo back and forth within the specimen. A continuous-wave oscillator is used to control the pulse repetition rate. When the repetition rate is adjusted so that the initiation of a pulse coincides with the return of the first echo from the preceding pulses, the change in the signal amplitude indicates superposition. The pulse rate is a measure of the travel time within the specimen.


    Figure | top | pdf |

    Block diagram of the pulse-echo technique.

  • (ii) `Sing around method'. The `sing around' method for measuring the velocity of ultrasonic waves involves the use of two piezoelectric transducers, one at each end of the specimen. One transducer receives an impulse from the electronic generator and converts it into an ultrasonic pulse in the specimen. This pulse, after passing through the specimen, is detected by the receiving transducer. The received pulse triggers the electronic generator to initiate a succeeding pulse. The pulse repetition rate is a very sensitive probe for measuring changes of the ultrasonic velocity in the specimen. Relative variations of [10^{-7}] can be measured, such as temperature or stress dependences of the velocity.

1.3.5. Pressure dependence and temperature dependence of the elastic constants

| top | pdf | Introduction

| top | pdf |

In a solid, the elastic constants are temperature and pressure dependent. As examples, the temperature dependence of the elastic stiffnesses of an aluminium single crystal within its stability domain (the melting point is 933 K) and the pressure dependence of the elastic stiffnesses of the ternary compound KZnF3 within its stability domain (the crystal becomes unstable for a hydrostatic pressure of about 20 GPa) are shown in Figs.[link] and[link], respectively.


Figure | top | pdf |

Temperature dependence of the elastic stiffnesses of an aluminium single crystal (after Every & McCurdy, 1992[link]).


Figure | top | pdf |

Pressure dependence of the elastic stiffness [c_{11} ] of a KZnF3 crystal. Reproduced with permission from Ultrasonics Symposium Proc. IEEE (Fischer et al., 1980[link]). Copyright (1980) IEEE.

We can observe the following trends, which are general for stable crystals:

  • (i) From 0 K to about [\Theta_{D}/5], where [\Theta_{D}] is the Debye temperature, the elastic stiffnesses decrease according to a [\Theta^{4}] law. From [\Theta_{D}/5] to the beginning of the instability domain, the dependence is linear with Θ. In addition, [(\partial c_{ii}/\partial \Theta)_{p} = 0] at 0 K as predicted by the third principle of thermodynamics.

  • (ii) For stable crystals, the pressure dependence of the elastic stiffnesses is linear as long as the applied pressure is small compared to the elastic stiffnesses. As an example, a typical order of magnitude for [c_{11}] in crystals is about 100 GPa and, within the experimental uncertainty, the pressure dependence of [c_{11} ] does not depart from a linear behaviour up to at least 0.2 GPa.

These observations can be quantitatively justified on the basis of an equation of state of a solid: [f(T_{ij},S_{ij},X,\Theta) = 0, ]where [T_{ij}] represents the stress tensor, [S_{ij}] the strain tensor, X the position of the elementary elements of the solid and Θ the temperature.

Different equations of state of solids have been proposed. They correspond to different degrees of approximation that can only be discussed and understood in a microscopic theory of lattice dynamics. The different steps in the development of lattice dynamics, the Einstein model, the Debye model and the Grüneisen model, will be presented in Section[link] . Concerning the temperature and the pressure dependences of the elastic constants, we may notice that rather sophisticated models are needed to describe correctly the general trends mentioned above:

  • (a) In the Einstein model, where the N atoms of a crystal are considered as 3N independent harmonic oscillators, the temperature increase affects only the amplitude of the oscillations and not the average positions of the atoms; consequently, this model can explain neither the thermal expansion nor the temperature dependence of the elastic constants. In addition, this theory ignores the difference between isothermal and adiabatic elastic constants. Similarly, if the oscillators are harmonic, the stiffness of the `springs' connecting atoms does not depend on the distances between atoms and the model cannot therefore explain the pressure dependence of the elastic constants, which requires anharmonic `springs' or, more accurately, anharmonic potentials.

  • (b) In the Debye model, the 3N oscillators are not independent but they are still harmonic. The result is that here again the elastic constants are pressure and temperature independent.

  • (c) In the Grüneisen model, the frequencies of the oscillators are volume-dependent; this so-called `quasiharmonic approximation' can justify the standard behaviour shown in Figs.[link] and[link]. Temperature dependence of the elastic constants

| top | pdf |

Table[link] gives typical values of [(\partial \ln c_{ij}/\partial \Theta)_{p}] for some cubic crystals considered within their stability domain. In column 6, the `elastic Debye temperature' of the crystal, [(\Theta_{D})_{\rm el} ], has been calculated according to the formula [(\Theta_{D})_{\rm el}= hv/k_{B}(3n/4\pi)^{1/3},]where h is the Planck constant, [k_{B}] is the Boltzmann constant, v is an average velocity (see for instance De Launay, 1956[link]) and n is the number of atoms per unit volume.

Table| top | pdf |
Temperature dependence of the elastic stiffnesses for some cubic crystals

MaterialTemperature range (K)[(\partial \ln c_{11}/\partial \Theta)_{p} ] (10−4 K−1)[(\partial \ln c_{44}/\partial \Theta)_{p} ] (10−4 K−1)[(\partial \ln c_{12}/\partial \Theta)_{p} ] (10−4 K−1)[\Theta^{\rm el}] (K)
Al 80–300 −3.1 −4.45 −1.3 430
Cu 80–300 −2.01 −3.33 −1.24 344
Ag 50–300 −2.3 −4.0 −1.5 226
Pb 100–300 −4.4 −1.5 −2.8 105
Si 80–300 −0.81 −10.6 −1.10 648
Ge 150–1000 −1.2 −1.15 −1.10 374
ZnS 100–300 −1.2 −0.65 −0.8 347
NaCl 100–300 −7.8 −2.2 −4.7 321
KCl 80–300 −8.3 −2.1 −3.6 236
KBr 80–300 −7.6 −2.1 7 172

It is interesting to compare [(\Theta_{D})_{\rm el}], the `elastic Debye temperature', with [(\Theta_{D})_{\rm cal} ], the `calorimetric Debye temperature'. The definition of [(\Theta_{D})_{\rm cal}] will be given in Section[link] . It results from the attempt at founding a universal description for the thermal properties of solids when the temperature is expressed as a reduced temperature, [\Theta /(\Theta_{D})_{\rm cal}]; [(\Theta_{D})_{\rm cal}] is obtained from calorimetric measurements at low temperature. It is worth noting that accurate values of low-temperature elastic constants and low-temperature calorimetric measurements lead to an excellent agreement between [(\Theta_{D})_{\rm el}] and [(\Theta_{D})_{\rm cal}] [better than 2 or 3 K (De Launay, 1956[link])]. This agreement demonstrates the validity of the Debye model in the vicinity of 0 K. From Table[link], we can observe that for ionic crystals [(\partial \ln c_{11}/\partial \Theta)_{p} ] is, in general, greater than [(\partial \ln c_{44}/\partial \Theta)_{p} ]. This remark is not valid for covalent and metallic crystals. Typical orders of magnitude are given in Table[link]. These statements concern only general trends valid for stable crystals.

Table| top | pdf |
Order of magnitude of the temperature dependence of the elastic stiffnesses for different types of crystals

Type of crystal[(\partial \ln c_{11}/\partial \Theta)_{p} ] (K−1)[(\partial \ln c_{44}/\partial \Theta)_{p} ] (K−1)
Ionic [-10^{-3}] [-3 \times 10^{-4}]
Covalent [-10^{-4}] [-8 \times 10^{-5}]
Metallic [-2 \times 10^{-4}] [-3 \times 10^{-4}]

In the case of temperature-induced phase transitions, some elastic constants are softened in the vicinity and sometimes far from the critical temperature. As an example, Fig.[link] shows the temperature dependence of [c_{44}] in RbCdF3, CsCdF3 and TlCdF3 single crystals. RbCdF3 and TlCdF3 undergo structural phase transitions at 124 and 191 K, respectively, while CsCdF3 remains stable in this temperature range. The softening of [c_{44}] when the temperature decreases starts more than 100 K before the critical temperature, [\Theta_{c}]. In contrast, Fig.[link] shows the temperature dependence of [c_{11}] in KNiF3, a crystal that undergoes a para–antiferromagnetic phase transition at 246 K; the coupling between the elastic and the magnetic energy is weak, consequently [c_{11}] decreases abruptly only a few degrees before the critical temperature. We can generalize this observation and state that the softening of an elastic constant occurs over a large domain of temperature when this constant is the order parameter or is strongly coupled to the order parameter of the transformation; for instance, in the cooperative Jahn–Teller phase transition in DyVO4, [(c_{11} - c_{12})/2] is the soft acoustic phonon mode leading to the phase transition and this parameter anticipates the phase transition 300 K before it occurs (Fig.[link]).


Figure | top | pdf |

Temperature dependence of the elastic constant [c_{44}] in RbCdF3, CsCdF3 and TlCdF3 crystals; the crystals of RbCdF3 and TlCdF3 undergo structural phase transitions (after Rousseau et al., 1975[link]).


Figure | top | pdf |

Temperature dependence of the elastic constant [c_{11}] in KNiF3, which undergoes a para–antiferromagnetic phase transition. Reprinted with permission from Appl. Phys. Lett. (Nouet et al., 1972[link]). Copyright (1972) American Institute of Physics.


Figure | top | pdf |

Temperature dependence of [(c_{11} - c_{12})/2] in DyVO4, which undergoes a cooperative Jahn–Teller phase transition (after Melcher & Scott, 1972[link]). Pressure dependence of the elastic constants

| top | pdf |

As mentioned above, anharmonic potentials are needed to explain the stress dependence of the elastic constants of a crystal. Thus, if the strain-energy density is developed in a polynomial in terms of the strain, only the first and the second elastic constants are used in linear elasticity (harmonic potentials), whereas higher-order elastic constants are also needed for nonlinear elasticity (anharmonic potentials).

Concerning the pressure dependence of the elastic constants (nonlinear elastic effect), considerable attention has been paid to their experimental determination since they are a unique source of significant information in many fields:

  • (i) In geophysics, a large part of the knowledge we have on the interior of the earth comes from the measurement of the transit time of elastic bursts propagating in the mantle and in the core (in the upper mantle, the average pressure is estimated to be about a few hundred GPa, a value which is comparable to that of the elastic stiffnesses of many materials).

  • (ii) In solid-state physics, the pressure dependence of the elastic constants gives significant indications concerning the stability of crystals. For example, Fig.[link] shows the pressure dependence of the elastic constants of KZnF3, a cubic crystal belonging to the perovskite family. As mentioned previously, this crystal is known to be stable over a wide range of temperature and the elastic stiffnesses [c_{ij}] depend linearly on pressure. It may be noted that, consequently, the third-order elastic constants (TOECs) are constant. On the contrary, we observe in Fig.[link] that the pressure dependence of the elastic constants of TlCdF3, a cubic crystal belonging to the same family but which is known to become unstable when the temperature is decreased to 191 K (Fischer, 1982[link]), is nonlinear even at low pressures. In this case, the development of the strain-energy density in terms of strains cannot be stopped after the terms containing the third-order elastic constants; the contributions of the fourth- and fifth-order elastic constants are not negligible.


    Figure | top | pdf |

    Pressure dependence of the elastic constants [(c_{11} - c_{12})/2] in TlCdF3. Reproduced with permission from Ultrasonics Symposium Proc. IEEE (Fischer et al., 1980[link]). Copyright (1980) IEEE.

  • (iii) For practical use in the case of technical materials such as concrete or worked metals, the pressure dependence of the elastic moduli is also required for examining the effect of applied stresses or of an applied hydrostatic pressure, and for studying residual stresses resulting from loading (heating) and unloading (cooling) the materials.

1.3.6. Nonlinear elasticity

| top | pdf | Introduction

| top | pdf |

In a solid body, the relation between the stress tensor T and the strain tensor S is usually described by Hooke's law, which postulates linear relations between the components of T and S (Section[link]). Such relations can be summarized by (see equation[link][T_{ij} = c_{ijkl} S_{kl}, ]where the [c_{ijkl}]'s are the elastic stiffnesses.

For a solid under finite strain conditions, Hooke's law, valid for infinitesimal deformations, does not hold, and the fundamental definitions for stress and strain must be revisited. Lagrangian and Eulerian descriptions

| top | pdf |

Finite elastic strains may be treated from two different viewpoints using either the Lagrangian (material) or the Eulerian (spatial) descriptions.

Let us consider a fixed rectangular Cartesian coordinate system with axes [x_{i}] ([i = 1, 2, 3]). Any particular position vector r of components ([x_{1}], [x_{2}], [x_{3}]) denotes a point in space. A point that always moves with the material is called a particle or material point. Let every particle be identified by its coordinates at some reference time [t_{0}]. These reference coordinates, referred to the same Cartesian system, will be denoted by ([a_{1}], [a_{2}], [a_{3}]) and the corresponding position vector a. A particular vector a can serve as a name for the particle located at that position at the reference time [t_{0} ].

The vectors r and a both specify a position in a fixed Cartesian frame of reference. At any time, we associate each r with an a by the rule that r is the present position vector of the particle initially at a. This connection between r and a is written symbolically as [{\bf r} = {\bf r}(t, a) \quad \hbox{ or } \quad x_{i} = x_{i} (t, a_{1}, a_{2}, a_{3}), \eqno( ]where [{\bf a} = {\bf r}(t_{0}, {\bf a}) \quad \hbox{ or } \quad a_{i} = x_{i} (t_{0}, a_{1}, a_{2}, a_{3}). \eqno( ]

The coordinates [a_{i}] that identify the particles are called material coordinates. A description that, like ([link], uses (t, [a_{1}], [a_{2}], [a_{3}]) as independent variables is called a material or Lagrangian description.

The converse of ([link] and ([link] may be written [{\bf a} = {\bf a} (t, {\bf r}) \quad \hbox{ or } \quad a_{i} = a_{i} (t, x_{1}, x_{2}, x_{3}), \eqno( ]where [{\bf r} = {\bf a}(t_{0}, {\bf r}) \quad \hbox{ or } \quad x_{i} = a_{i} (t_{0}, x_{1}, x_{2}, x_{3}). \eqno( ]

A spatial description or Eulerian description uses the independent variables (t, [x_{1}], [x_{2}], [x_{3}]), the [x_{i}] being called spatial coordinates.

Now, for the sake of simplicity, we shall work with the Lagrangian formulation exclusively. For more details see, for instance, Thurston (1964)[link] and Wallace (1970[link], 1972[link]). Strain and stress tensors

| top | pdf |

The displacement vector from the reference position of a particle to its new position has as components [u_{i} = x_{i} - a_{i}. \eqno( ]

The term strain refers to a change in the relative positions of the material points in a body. Let a final configuration be described in terms of the reference configuration by setting t equal to a constant in ([link]. Then t no longer appears as a variable and ([link] can be written [x_{i} = x_{i} (a_{1}, a_{2}, a_{3}), ]where the [a_{i}] are the independent variables. It follows that[\hbox{d}x_{j} = {\partial x_{j} \over \partial a_{i}} \hbox{d} a_{i} = \left({\partial u_{j} \over \partial a_{i}} + \delta_{ij}\right) \hbox{d}a_{i}. \eqno( ]

Let now the particle initially at ([a_{1}], [a_{2}], [a_{3}]) move to ([x_{1}], [x_{2}], [x_{3}]). The square of the initial distance to a neighbouring particle whose initial coordinates were [a_{j} + \hbox{d}a_{j}] is [\hbox{d}s^{2} = \hbox{d}a_{j} \hbox{d}a_{j}. ]The square of the final distance to the same neighbouring particle is [\hbox{d}s^{2} = \hbox{d}x_{j} \hbox{d}x_{j}.]

In a material description, the strain components [S_{ik} ] are defined by the following equations: [\hbox{d}x_{j} \hbox{d}x_{j} - \hbox{d}a_{j} \hbox{d}a_{j} = 2 S_{ik}\hbox{d}a_{i} \hbox{d}a_{k}. \eqno( ]Substituting ([link] into ([link], it follows that [\left({\partial u_{j} \over \partial a_{i}} + \delta_{ji}\right) \left({\partial u_{j} \over \partial a_{k}} + \delta_{jk}\right) \hbox{d}a_{i} \hbox{d}a_{k} - \hbox{d}a_{j} \hbox{d}a_{j} = 2 S_{ik}\hbox{d}a_{i} \hbox{d}a_{k}. ]Hence [S_{ik} = \textstyle{1\over 2} \left(\displaystyle{\partial u_{k} \over \partial a_{i}} + {\partial u_{i} \over \partial a_{k}} + {\partial u_{j} \over \partial a_{i}} {\partial u_{j} \over \partial a_{k}}\right). ]

If the products and squares of the displacement derivatives are neglected, the strain components reduce to the usual form of `infinitesimal elasticity' [see equation ([link]]: [S_{ik} = \textstyle{1 \over 2} \left(\displaystyle{\partial u_{i} \over \partial a_{k}} + {\partial u_{k} \over \partial a_{i}}\right). ]

It is often useful to introduce the Jacobian matrix associated with the transformation (a, x). The components of this matrix are [J = \pmatrix{\alpha_{11} &\alpha_{12} &\alpha_{13}\cr \alpha_{21} &\alpha_{22} &\alpha_{23}\cr \alpha_{31} &\alpha_{32} &\alpha_{33}\cr}, ]where [\alpha_{ik} = {\partial x_{i} \over \partial a_{k}} = {\partial u_{j} \over \partial a_{k}} + \delta_{jk}. ]

From the definition of matrix J, one has [\hbox{d}{\bf x} = J\hbox{d}{\bf a}]and [\hbox{d}x^{2} - \hbox{d}a^{2} = \left(\hbox{d}{\bf x}\right)^{T} \hbox{d}{\bf x} - \left(\hbox{d}{\bf a}\right)^{T} \hbox{d}{\bf a} = \left(\hbox{d}{\bf a}\right)^{T} \left(J^{T} J - \delta\right)\hbox{d}{\bf a}, ]where [\left(\hbox{d}{\bf a}\right)^{T}], [\left(\hbox{d}{\bf x}\right)^{T}] and [J^{T}] are the transpose matrices of da, dx and J, respectively, and δ is the Kronecker matrix.

The Lagrangian strain matrix S may then be written symbolically: [S = \textstyle{1\over 2} \left(J^{T}J - \delta\right). \eqno( ]

When finite strains are concerned, we have to distinguish three states of the medium: the natural state, the initial state and the final or present state: The natural state is a state free of stress. The initial state is deduced from the natural state by a homogeneous strain. The final state is deduced from the initial state by an arbitrary strain.

Concerning the stress tensor, as pointed out by Thurston (1964)[link], the stress-deformation relation is complicated in nonlinear elasticity because `the strain is often referred to a natural unstressed state, whereas the stress [T_{ij}] is defined per unit area of the deformed body'. For this reason, the differential of work done by the stress is not equal to the stress components times the differentials of the corresponding strain components. So, following Truesdell & Toupin (1960)[link], we shall introduce a thermodynamic tension tensor [t_{ij}] defined as the first derivative of the energy with respect to strain. If the internal energy U per unit mass is considered, the thermodynamic tension refers to an isentropic process. Then [t_{ij}^{\sigma} = \rho_{0} \left({\partial U \over \partial S_{ij}}\right)_{\sigma}, ]where σ is the entropy and [\rho_{0}] the volumic mass in the initial state.

If the Helmholtz free energy F is considered, the thermodynamic tension refers to an isothermal process. Then[t_{ij}^{\Theta} = \rho_{0} \left({\partial U \over \partial S_{ij}}\right)_{\Theta}, ]where Θ is the temperature. It will be shown in Section[link] that [T_{ij} = (1/J)\alpha_{ik}\alpha_{jl}t_{kl}. ] Second-order and higher-order elastic stiffnesses

| top | pdf |

Following Brugger (1964)[link], the strain-energy density, or strain energy per unit volume Φ, is assumed to be a polynomial in the strain: [\Phi = \Phi_{0} + c_{ij}S_{ij} + {1\over 2!}c_{ijkl}S_{ij}S_{kl} + {1\over 3!}c_{ijklmn}S_{ij}S_{kl}S_{mn}, \eqno( ]where [\Phi = \rho_{0} U(X,S_{ij})], [\Phi_{0} = \rho_{0} U(X,0)], X denotes the configuration of the initial state and the [S_{ij}]'s are the Lagrangian finite strain-tensor components.

If the initial energy and the deformation of the body are both zero, the first two terms in ([link] are zero. Note that [c_{ij}] is a stress and not an intrinsic characteristic of the material. In this expression, the elastic stiffnesses [c_{ijkl}] and [c_{ijklmn}] are the second- and third-order stiffnesses, respectively. Since the strain tensor is symmetric, pairs of subscripts can be interchanged [see equation ([link]]: [\eqalign{c_{ijkl} &= c_{jikl} = c_{ijlk} = c_{jilk},\cr c_{ijklmn} &= c_{jiklmn} = c_{ijlkmn} = c_{jilkmn} = c_{ijklnm} \cr&= c_{jiklnm} = c_{ijlknm} = c_{jilknm}.\cr} ]

More accurately, the isentropic and the isothermal elastic stiffnesses are defined as the nth partial derivatives of the internal energy and the Helmholtz free energy, respectively. For example, the third-order isentropic and isothermal stiffnesses are, respectively, [\eqalign{c_{ijkl mn}^{\sigma} &= \rho_{0} {\partial^{3} U \over \partial S_{ij}\partial S_{kl}\partial S_{mn}},\cr c_{ijkl mn}^{\Theta} &= \rho_{0} {\partial^{3} F \over \partial S_{ij}\partial S_{kl}\partial S_{mn}},\cr} ]where the internal energy, U, is a function of X, [S_{ij}] and σ, and the Helmholtz free energy, F, is a function of X, [S_{ij}] and Θ.

From these definitions, it follows that the Brugger stiffness coefficients depend on the initial state. When no additional information is given, the initial state is the natural state.

The third-order stiffnesses form a sixth-rank tensor containing 36 = 729 components, of which 56 are independent for a triclinic crystal and 3 for isotropic materials. The independent components of a sixth-rank tensor can be obtained for any point group using the accompanying software to this volume. They are listed in Table[link] for each Laue class. They can be reduced to three-index coefficients. Table[link] gives the values of the third-order stiffnesses of some materials.

Table| top | pdf |
Number of independent third-order elastic stiffnesses for each Laue class

Laue classNumber of independent components
[{\bar 1}], 1 56
2/m, 2, m 32
mmm, 222, 2mm 20
[{\bar 3}], 3 20
[{\bar 3}], 32, 3m 14
4/m, 4, [{\bar 4}] 16
4/mmm, 422, [{\bar 4}2m], 42m 12
6/m, 6, [{\bar 6}] 12
6/mmm, 622, [{\bar 6}2m], 62m 10
[m{\bar 3}], 23 8
[m{\bar 3}m], 432, [{\bar 4}32] 6
[\infty (A_\infty / M) C], [\infty A_\infty] 3

Table| top | pdf |
Third-order elastic stiffnesses of some materials in (GPa)−1 (after Every & McCurdy, 1992[link])

(a) Trigonal.

Al2O3 −3870 −1090 −963 55 −289 −39 −922 −131 −302 −1160 −4520 −3340 −1090 −19
CaCO3 −579 −147 −193 218 −41 10 −239 82 −69 −139 −675 −498 −195 33

(b) Tetragonal.

KH2PO4 −538 −291 −259 −258 67 −8 −35 −30 −912 −27 −19 8
TeO2 −160 −600 −140 −110 180 −41 36 −640 −2110 −54 −260 −250

(c) Hexagonal.

Cd −2060 −114 −197 −110 −268 227 −332 −2020 −516 −171
CdS −459 −207 −182 −235 −306 −27 9 −355 −327 −69
Co −6710 −1454 −766 −429 −511 133 −1486 −5788 −6347 −210
Mg −663 −178 30 −76 −86 −30 −58 −864 −726 −193
Zn −1760 −440 −270 −210 −350 −10 250 −2410 −720 −440

(d) Cubic.

LiF − 1920 − 330 −40 100 −325 43
MgO − 4900 −95 −69 113 −659 147
KBr − 532 −49 69 22 −28 −35
KCl − 610 −31 9 19 −35 17
NaCl − 866 −37 38 25 −78 18
Diamond − 6260 −2260 112 −674 −2860 −823
GaAs − 620 −392 −62 8 −274 −43
Ge − 714 −388 −34 −9 −303 −48
Si − 795 −445 −75 15 −310 −86
ZnSe − 827 −136 −511 222 −265 −278
Al − 1224 −373 25 −64 −368 −27
Cu − 1271 −814 −50 −3 −780 −95
Au − 1730 −922 −233 −13 −648 −12
Fe − 2705 −626 −575 −836 −531 −721
Ni − 2104 −1345 59 −180 −757 −42

The three independent constants for isotropic materials are often taken as [c_{123}], [c_{144}] and [c_{456}] and denoted respectively by [\nu_{1} ], [\nu_{2}], [\nu_{3}], the `third-order Lamé constants'.

The `third-order Murnaghan constants' (Murnaghan, 1951[link]), denoted by [l, m, n], are given in terms of the Brugger constants by the relations [l = \textstyle{1\over 2} c_{112} ; \ m = c_{155} ; \ n = 4c_{456}. ]

Similarly, the fourth-order stiffnesses form an eighth-rank tensor containing [3^{8} = 6561] components, 126 of which are independent for a triclinic crystal and 11 for isotropic materials (the independent components of a sixth-rank tensor can be obtained for any point group using the accompanying software to this volume).

For a solid under finite strain conditions, the definition of the elastic compliance tensor has to be reconsidered. In linear elasticity, the second-order elastic compliances [s_{ijkl}] were defined through the relations ([link]: [S_{ij} = s_{ijkl} T_{kl} \quad \hbox{or} \quad s_{ijkl} = {\partial S_{ij} \over \partial T_{kl}}, ]while, in nonlinear elasticity, one has [s_{ijkl} = {\partial S_{ij} \over \partial t_{kl}}, ]where [t_{kl} = \rho_0 {\partial U \over \partial S_{kl}}. ] Expansion of elastic constants for small initial stress

| top | pdf |

In most experiments, the initial stress is small compared with the second-order elastic constants (for example, 1 GPa hydrostatic pressure compared with the usual value cijkl = 100 GPa). Consequently, the deformation between the initial (stressed) state and the natural (unstressed) state is small compared with 1. For this reason, it is convenient to expand the elastic constants in the initial state as a power series in the strain about the natural state. To avoid confusion, we introduce new notations: [\bar{X}] now represents the coordinates in the natural or unstressed state; X represents the coordinates in the initial or homogeneously strained state; [u_{i} = x_{i} - X_{i} ] are the components of displacement. All letters with superscript bar refer to the natural state; for example, [\bar{S}_{ij}] denotes the Lagrangian strain in the natural state; [\bar{U} = U(\bar{X}, \bar{S}_{ij}) ].

Now, in order to relate the properties at X to those at [\bar{X}], we need to specify the strain from [\bar{X}] to X. Let [a_{ij} = {\partial X_{i} \over\partial \bar{X}_{j}} = \bar{\alpha}_{ij}. ]Consequently, [{\partial \bar{S}_{ij} \over \partial S_{mn}} = a_{mi}a_{nj}. ]The second-order elastic constants at X can be expressed in terms of the second- and third-order elastic constants at [\bar{X}]: [c_{ijkl} = \rho_{0} {\partial^{2} U \over \partial S_{ij}\partial S_{kl}} = \rho_{0} {\partial^{2}\bar{U} \over \partial \bar{S}_{mn}\partial \bar{S}_{pq}} a_{im}a_{jn}a_{kp}a_{l q} ]or [c_{ijkl} = {\rho_{0} \over \bar{\rho}_0} \left(\bar{c}_{mnpq} + \bar{c}_{mnpqrs}S_{rs} + {1\over 2!} \bar{c}_{mnpqrstu}S_{rs}S_{tu} + \dots\right) a_{im}a_{jn}a_{kp}a_{l q}. ]This expression holds for both isentropic and isothermal elastic constants. Elastic strain-energy density

| top | pdf |

The elastic strain-energy density has appeared in the literature in various forms. Most of the authors use the Murnaghan constants as long as isotropic solids are concerned. However, most of the literature uses Brugger's thermodynamic definition when anisotropic media are under consideration (Brugger, 1964[link]).

The elastic strain-energy density for an isotropic medium, including third-order terms but omitting terms independent of strain, may be expressed in terms of three strain invariants, since an isotropic material is invariant with respect to rotation: [\Phi = {\lambda + 2\mu \over 2} (I_1)^{2} - 2\mu I_{2} + {l + 2m \over 3} (I_{1})^{3} - 2 m I_{1} I_{2} + nI_{3}, ]where λ and μ are the second-order Lamé constants, [l, m, n] are the third-order Murnaghan constants, and [I_{1}], [I_{2}], [I_{3}] are the three invariants of the Lagrangian strain matrix. These invariants may be written in terms of the strain components as [\eqalignno{I_{1} &= S_{11} + S_{22} + S_{33} &\cr I_{2} &= \left| \matrix{S_{11} &S_{12}\cr S_{21} &S_{22}\cr}\right| + \left| \matrix{S_{22} &S_{23}\cr S_{32} &S_{33}\cr}\right| + \left| \matrix{S_{33} &S_{31}\cr S_{13} &S_{11}\cr}\right| &\cr I_{3} &= \left| \matrix{S_{11} &S_{12} &S_{13}\cr S_{21} &S_{22} &S_{23}\cr S_{31} &S_{32} &S_{33}\cr}\right|. &\cr} ]The elastic strain-energy density for an anisotropic medium (for example a medium belonging to the most symmetrical groups of cubic crystals) is (Green, 1973[link]) [\eqalignno{\Phi &= \textstyle{1\over 2} c_{11} \left[(S_{11})^{2} + (S_{22})^{2} + (S_{33})^{2}\right] + c_{12} \left[S_{11}S_{22} + S_{22}S_{33} + S_{33}S_{11}\right] &\cr &\quad + c_{44}[(S_{12})^{2} + (S_{21})^{2} + (S_{23})^{2} +(S_{32})^{2} + (S_{31})^{2} + (S_{13})^{2}] &\cr &\quad + c_{111} \left[(S_{11})^{3} + (S_{22})^{3} + (S_{33})^{3}\right] &\cr &\quad + c_{112} \left[(S_{11})^{2}(S_{22} + S_{33}) + (S_{22})^{2}(S_{33} + S_{11})\right.\cr&\quad + \left.(S_{33})^{2}(S_{11} + S_{22})\right] &\cr &\quad + \textstyle{1\over 2} c_{144} \left\{S_{11}[(S_{23})^{2} + (S_{32})^{2}] + S_{22}[(S_{31})^{2} + (S_{13})^{2}] \right.\cr&\quad + \left. S_{33}[(S_{12})^{2} + (S_{21})^{2}]\right\} &\cr &\quad + \textstyle{1\over 2}c_{166} \left\{\left[(S_{12})^{2} + (S_{21})^{2}\right](S_{11} + S_{22})\right. \cr&\quad+ \left[(S_{23})^{2}+(S_{32})^{2}\right](S_{22} + S_{33}) &\cr &\quad + \left.\left[(S_{13})^{2} + (S_{31})^{2}\right] (S_{11} + S_{33})\right\} &\cr &\quad + c_{123}S_{11}S_{22}S_{33} + c_{456} \left[S_{12}S_{23}S_{31} + S_{21}S_{32}S_{13}\right]. &\cr} ]

1.3.7. Nonlinear dynamic elasticity

| top | pdf | Introduction

| top | pdf |

In recent years, the measurements of ultrasonic wave velocities as functions of stresses applied to the sample and the measurements of the amplitude of harmonics generated by the passage of an ultrasonic wave throughout the sample are in current use. These experiments and others, such as the interaction of two ultrasonic waves, are interpreted from the same theoretical basis, namely nonlinear dynamical elasticity.

A first step in the development of nonlinear dynamical elasticity is the derivation of the general equations of motion for elastic waves propagating in a solid under nonlinear elastic conditions. Then, these equations are restricted to elastic waves propagating either in an isotropic or in a cubic medium. The next step is the examination of two important cases:

  • (i) the generation of harmonics when finite amplitude ultrasonic waves travel throughout an unstressed medium;

  • (ii) the propagation of small amplitude ultrasonic waves when they travel throughout a stressed medium.

Finally, the concept of natural velocity is introduced and the experiments that can be used to determine the third- and higher-order elastic constants are described. Equation of motion for elastic waves

| top | pdf |

For generality, these equations will be derived in the X configuration (initial state). It is convenient to obtain the equations of motion with the aid of Lagrange's equations. In the absence of body forces, these equations are [{{\rm d} \over {\rm d}t}{\partial L \over \partial x^{\prime}_i} + {\partial \over \partial X_i}{\partial L \over \partial \left( \partial x_i/\partial X_j\right)} = 0 \eqno( ]or [{{\rm d} \over {\rm d}t}{\partial L \over \partial x^{\prime}_i} + {\partial \over \partial X_i}{\partial L \over \partial \alpha_{ij}} = 0, \eqno( ]where L is the Lagrangian per unit initial volume and [\alpha_{ij} =] [ \partial x_i/\partial X_j] are the elements of the Jacobian matrix.

For adiabatic motion[L = \textstyle{1\over 2}\displaystyle \rho_0 x_i^{'2} - \rho_0 U, \eqno( ]where U is the internal energy per unit mass.

Combining ([link] and ([link], it follows that [\rho_0 x{''}_i^{2} = {\partial \over \partial X_j} \left( \rho_0 {\partial U \over \partial S_{l m}} {\partial S_{l m}\over \partial \alpha_{ij}}\right), ]which can be written [\rho_0 x{''}_i^{2} = {\partial \over \partial X_j} \left( \alpha_{il}\alpha_{jm}\rho_0 {\partial U \over \partial S_{l m}}\right) ]since [{\partial S_{l m}\over \partial \alpha_{ij}} = \textstyle{1\over 2}\displaystyle \left( \alpha_{im}\delta_{jl} + \alpha_{il} \delta_{jm}\right). ]

Using now the equation of continuity or conservation of mass: [{\rho_0\over \rho}= J = {\rm det} (a_{ij}),]and the identity of Euler, Piola and Jacobi: [{\partial \over \partial x_j} \left( {1\over J} {\partial x_j \over \partial X_i}\right) = 0, ]we get an expression of Newton's law of motion: [\rho x_i'' = {{\rm d}T_{ij}\over {\rm d} X_j}\ {\rm or} \ \rho u_i'' = {{\rm d}T_{ij}\over {\rm d} X_j} \eqno( ]with [T_{ij} = {\rho_0 \over J} \alpha_{ik}\alpha_{jl} {\partial U \over \partial S_{kl}} = \rho \alpha_{ik}\alpha_{jl} {\partial U \over \partial S_{kl}}. ][T_{ij}] becomes [T_{ij} = {1\over J} \alpha_{ik} \alpha_{jl} t_{kl} ]since [t_{kl} = \rho_0 {\partial U \over \partial S_{kl}}. ][t_{kl}], the thermodynamic tensor conjugate to the variable [S_{kl}/\rho_ 0], is generally denoted as the `second Piola–Kirchoff stress tensor'.

Using Φ, the strain energy per unit volume, Newton's law ([link] takes the form [\rho x_i'' = {\partial \over \partial X_j} \left(\alpha_{jk} {\partial \Phi \over \partial S_{ik}}\right) \ \ {\rm or} \quad \rho u_i'' = {\partial \over \partial X_j} \left(\alpha_{jk} {\partial \Phi \over \partial S_{ik}}\right) ]and [T_{ij} = \alpha_{jk}{\partial \Phi \over \partial S_{ik}}. \eqno( ] Wave propagation in a nonlinear elastic medium

| top | pdf |

As an example, let us consider the case of a plane finite amplitude wave propagating along the [x_1] axis. The displacement components in this case become [u_1 = u_1(X_1,t);\quad u_2 = u_2(X_1,t)\semi\quad u_3 = u_3(X_1,t). ]Thus, the Jacobian matrix [\alpha_{ij}] reduces to [J = \pmatrix {\alpha_{11} &0 &0\cr \alpha_{21} &0 &0\cr \alpha_{31} &0 &0\cr}. ]

The Lagrangian strain matrix is [equation ([link]] [S = \textstyle{1\over 2}\displaystyle \left(J^T J - \delta \right). ]The only nonvanishing strain components are, therefore, [\eqalign{S_{11} &= \textstyle{1\over 2}\displaystyle \left(\alpha_{11}^2 + \alpha_{21}^2 + \alpha_{31}^2\right) -1\cr &= {\partial u_1\over \partial X_1} + \textstyle{1\over 2}\displaystyle \left[\left({\partial u_1\over \partial X_1}\right)^2 + \left({\partial u_2\over \partial X_1}\right)^2 + \left({\partial u_3\over \partial X_1}\right)^2\right]\cr S_{12} &= S_{21} = \textstyle{1\over 2}\displaystyle{\partial u_2\over \partial X_1}\cr S_{13} &= S_{31} = \textstyle{1\over 2}\displaystyle{\partial u_3\over \partial X_1}\cr} ]and the strain invariants reduce to [I_1 = S_{11} ; \ \ I_2 = - (S_{12}S_{21} + S_{13}S_{31}); \ \ I_3 = 0. ] Isotropic media

| top | pdf |

In this case, the strain-energy density becomes [\eqalignno{\Phi &= \textstyle{1\over 2}\displaystyle (\lambda+2\mu)(S_{11})^2 + 2 \mu(S_{12}S_{21} + S_{13}S_{31}) + \textstyle{1\over 3}\displaystyle (l+2m)(S_{11}){}^3&\cr &\quad + 2m S_{11}\left(S_{12}S_{21} + S_{13}S_{31}\right). &(\cr} ]Differentiating ([link] with respect to the strains, we get [\openup4.5pt\eqalign{{\partial \Phi \over \partial S_{11}} & = (\lambda+2\mu)S_{11} + (l+2m)(S_{11}){}^2 + 2m \left(S_{12}S_{21} + S_{13}S_{31}\right)\cr {\partial \Phi \over \partial S_{12}} & =2 \mu S_{21} + 2m S_{11}S_{21}\cr {\partial \Phi \over \partial S_{13}} & = 2\mu S_{31} + 2m S_{11}S_{31}\cr {\partial \Phi \over \partial S_{21}} & = 2 \mu S_{12} + 2m S_{11}S_{12}\cr {\partial \Phi \over \partial S_{31}} & =2\mu S_{13} + 2m S_{11}S_{13}.\cr} ]All the other [{\partial \Phi / \partial S_{ij}} = 0].

From ([link], we derive the stress components: [\openup6pt\eqalign{T_{11} &= \alpha_{1k}{\partial \Phi \over \partial S_{1k}}\semi \ \ T_{12} = \alpha_{2k}{\partial \Phi \over \partial S_{1k}}\semi \ \ T_{13} = \alpha_{3k}{\partial \Phi \over \partial S_{1k}}\semi \cr T_{21} &= \alpha_{1k}{\partial \Phi \over \partial S_{2k}}\semi \ \ T_{22} = \alpha_{2k}{\partial \Phi \over \partial S_{2k}}\semi \ \ T_{23} = \alpha_{3k}{\partial \Phi \over \partial S_{2k}}\semi\cr T_{31} &= \alpha_{1k}{\partial \Phi \over \partial S_{3k}}\semi \ \ T_{32} = \alpha_{2k}{\partial \Phi \over \partial S_{3k}}\semi \ \ T_{33} = \alpha_{3k}{\partial \Phi \over \partial S_{3k}}.\cr} ]Note that this tensor is not symmetric.

For the particular problem discussed here, the three components of the equation of motion are [\eqalign{\rho u_1'' &= {\rm d}T_{11}/{\rm d}X_1,\cr \rho u_2'' &= {\rm d}T_{21}/{\rm d}X_1,\cr \rho u_3'' &= {\rm d}T_{31}/{\rm d}X_1.\cr} ]

If we retain only terms up to the quadratic order in the displacement gradients, we obtain the following equations of motion: [\eqalign{\rho u_{1}'' &= (\lambda + 2\mu){\partial^{2} u_{1} \over \partial X_{1}^{2}} + [3(\lambda + 2\mu) + 2(l + 2m)] {\partial u_{1} \over \partial X_{1}}{\partial^{2} u_{1} \over \partial X_{1}^{2}}\cr &\quad + (\lambda + 2\mu + m) \left[{\partial u_{2} \over \partial X_{1}}{\partial^{2} u_{2} \over \partial X_{1}^{2}} + {\partial u_{3} \over \partial X_{1}}{\partial^{2} u_{3} \over \partial X_{1}^{2}}\right]\cr \rho u_{2}'' &= \mu{\partial^{2} u_{2} \over \partial X_{1}^{2}} + (\lambda + 2\mu + m) \left[{\partial u_{1} \over \partial X_{1}}{\partial^{2} u_{2} \over \partial X_{1}^{2}} + {\partial u_{2} \over \partial X_{1}}{\partial^{2} u_{1} \over \partial X_{1}^{2}}\right]\cr \rho u_{3}'' &= \mu{\partial^{2} u_{3} \over \partial X_{1}^{2}} + (\lambda + 2\mu + m)\left[{\partial u_{1} \over \partial X_{1}}{\partial^{2} u_{3} \over \partial X_{1}^{2}} + {\partial u_{3} \over \partial X_{1}}{\partial^{2} u_{1} \over \partial X_{1}^{2}}\right].\cr} \eqno( ] Cubic media (most symmetrical groups)

| top | pdf |

In this case, the strain-energy density becomes [\eqalignno{\Phi &= \textstyle{1\over 2}\displaystyle c_{11} (S_{11})^{2} + c_{44} \left[(S_{12})^{2} + (S_{21})^{2} + (S_{31})^{2} + (S_{13})^{2}\right]&\cr &\quad+ c_{111}(S_{11})^{3} + \textstyle{1\over 2}\displaystyle c_{166} S_{11} \left[(S_{12})^{2} + (S_{21})^{2} + (S_{31})^{2}+ (S_{13})^{2}\right].\cr &&(\cr} ]Differentiating ([link] with respect to the strain, one obtains [\eqalign{{\partial \Phi \over S_{11}} &= c_{11} S_{11} + 3c_{111}(S_{11})^{2} + \textstyle{1\over 2}\displaystyle c_{166} [(S_{12})^{2} + (S_{21})^{2}\cr&\quad + (S_{31})^{2} + (S_{13})^{2})]\cr {\partial \Phi \over S_{21}} &= 2c_{44}S_{21} + c_{166}S_{11}S_{21}\cr \noalign{\vskip6pt}{\partial \Phi \over S_{31}} &= 2c_{44}S_{31} + c_{166}S_{11}S_{31}.\cr} ]All other [\partial \Phi / S_{ij} = 0]. From ([link], we derive the stress components: [\eqalign{T_{11} &= \alpha_{1k}{\partial \Phi \over \partial S_{1k}}\cr \noalign{\vskip5pt} T_{21} &= \alpha_{1k}{\partial \Phi \over \partial S_{2k}}\cr\noalign{\vskip5pt} T_{31} &= \alpha_{1k}{\partial \Phi \over \partial S_{3k}}.\cr} ]In this particular case, the three components of the equation of motion are [\eqalign{\rho u_{1}'' &= \hbox{d}T_{11}/\hbox{d}X_{1}\cr \rho u_{2}'' &= \hbox{d}T_{21}/\hbox{d}X_{1}\cr \rho u_{3}'' &= \hbox{d}T_{31}/\hbox{d}X_{1}.\cr} ]

If we retain only terms up to the quadratic order in the displacement gradients, we obtain the following equations of motion: [\eqalign{\rho u_{1}'' &= c_{11}{\partial^{2} u_{1} \over \partial X_{1}^{2}} + [3c_{11} + c_{111}] {\partial u_{1} \over \partial X_{1}}{\partial^{2} u_{1} \over \partial X_{1}^{2}}\cr &\quad + (c_{11} + c_{166}) \left[{\partial u_{2} \over \partial X_{1}}{\partial^{2} u_{2} \over \partial X_{1}^{2}} + {\partial u_{3} \over \partial X_{1}}{\partial^{2} u_{3} \over \partial X_{1}^{2}}\right]\cr \rho u_{2}'' &= c_{44}{\partial^{2} u_{2} \over \partial X_{1}^{2}} + (c_{11} + c_{166}) \left[{\partial u_{1} \over \partial X_{1}}{\partial^{2} u_{2} \over \partial X_{1}^2} + {\partial u_{2} \over \partial X_{1}}{\partial^{2} u_{1} \over \partial X_{1}^{2}}\right]\cr \rho u_{3}'' &= c_{44}{\partial^{2} u_{3} \over \partial X_{1}^{2}} + (c_{11} + c_{166}) \left[{\partial u_{1} \over \partial X_{1}}{\partial^{2} u_{3} \over \partial X_{1}^{2}} + {\partial u_{3} \over \partial X_{1}}{\partial^{2} u_{1} \over \partial X_{1}^{2}}\right],\cr} \eqno( ]which are identical to ([link] if we put [c_{11} = \lambda + 2\mu\semi \quad c_{44} = \mu\semi \quad c_{111} = 2(l + 2m)\semi \quad c_{166} = m. ] Harmonic generation

| top | pdf |

The coordinates in the medium free of stress are denoted either a or [{\bar X}]. The notation [{\bar X}] is used when we have to discriminate the natural configuration, [{\bar X}], from the initial configuration X. Here, the process that we describe refers to the propagation of an elastic wave in a medium free of stress (natural state) and the coordinates will be denoted [a_{i}].

Let us first examine the case of a pure longitudinal mode, i.e. [u_{1} = u_{1}(a_{1},t)\semi \quad u_{2} = u_{3} = 0. ]

The equations of motion, ([link] and ([link], reduce to [\rho u_{1}'' = (\lambda + 2\mu){\partial^{2} u_{1} \over \partial a_{1}^{2}} + [3(\lambda + 2\mu) + 2(l + 2m)] {\partial u_{1}\over \partial a_{1}}{\partial^{2} u_{1}\over \partial a_{1}^{2}} ]for an isotropic medium or [\rho u_{1}'' = c_{11}{\partial^{2} u_{1}\over \partial a_{1}^{2}} + [3c_{11} + c_{166}] {\partial u_{1}\over \partial a_{1}}{\partial^{2}u_{1}\over \partial a_{1}^{2}} ]for a cubic crystal (most symmetrical groups) when a pure longitudinal mode is propagated along [100].

For both cases, we have a one-dimensional problem; ([link] and ([link] can therefore be written [\rho u_{1}'' = K_{2}{\partial^{2}u_{1}\over \partial a_{1}^{2}} + [3K_{2} + K_{3}] {\partial u_{1}\over \partial a_{1}}{\partial^{2}u_{1}\over \partial a_{1}^{2}}. \eqno( ]

The same equation is also valid when a pure longitudinal mode is propagated along [110] and [111], with the following correspondence: [\eqalign{[100]\ K_{2} &= c_{11}, \quad K_{3} = c_{111}\cr [110]\ K_{2} &= {c_{11} + c_{12} + 2c_{44}\over 2}, \quad \ K_{3} = {c_{111} + 3c_{112} +12 c_{166} \over 4}\cr [111]\ K_2 &= {c_{11} + 2c_{12} + 4c_{44}\over 3},\cr K_{3} &= {c_{111} + 6c_{112} + 12c_{144} + 24c_{166} + 2c_{123} + 16c_{456}\over 9}.\cr} ]Let us assume that [K_{3} \ll K_{2}]; a perturbation solution to ([link] is [u = u^{0} + u^{1},]where [u^{1} \ll u^{0}] with [\eqalignno{u^{0} &= A \sin (ka - \omega t) &(\cr u^{1} &= Ba \sin 2(ka - \omega t) + C a \cos 2(ka - \omega t). &(\cr} ]

If we substitute the trial solutions into ([link], we find after one iteration the following approximate solution: [u = A \sin (ka - \omega t) - {(kA)^{2}(3K_{2} + K_{3})\over 8\rho c^{2}} a \cos 2(ka - \omega t), ]which involves second-harmonic generation.

If additional iterations are performed, higher harmonic terms will be obtained. A well known property of the first-order nonlinear equation ([link] is that its solutions exhibit discontinuous behaviour at some point in space and time. It can be seen that such a discontinuity would appear at a distance from the origin given by (Breazeale, 1984[link]) [L = - 2{(K_2)^{2}\over 3K_{2} + K_{3}}\rho \omega u_{0}', ]where [u_{0}'] is the initial value for the particle velocity. Small-amplitude waves in a strained medium

| top | pdf |

We now consider the propagation of small-amplitude elastic waves in a homogeneously strained medium. As defined previously, [\bar{X}] or a are the coordinates in the natural or unstressed state. X are the coordinates in the initial or homogeneously strained state. [u_{i} = x_{i} - X_{i}] are the components of displacement from the initial state due to the wave.

Starting from ([link], we get [T_{ij} = {\rho_{0} \over J} \alpha_{ik}\alpha_{jl} {\partial U \over \partial S_{kl}}. ]Its partial derivative is [{\partial T_{ij} \over \partial x_{j}} = {1 \over J}{\partial \over \partial X_{k}} \left[\rho_{0} \alpha_{il} {\partial U \over \partial S_{kl}}\right]. ]If we expand the state function about the initial configuration, it follows that [\eqalign{\rho_{0} U(X_{k}, S_{ij}) =&\ \rho_{0} U(X_{k}) + c_{ij}S_{ij} + \textstyle{1\over 2} \displaystyle c_{ijkl}S_{ij}S_{kl} \cr &+ \textstyle{1\over 6} \displaystyle c_{ijklmn}S_{ij}S_{kl}S_{mn} + \ldots.\cr} ]

The linearized stress derivatives become [{\partial T_{ij} \over \partial x_{j}} = \left[c_{jl}\delta_{ik} + c_{ijkl}\right] {\partial^{2} x_{k} \over \partial X_{j} \partial X_{l}}. ]If we let [D_{ijkl} = \left[c_{jl}\delta_{ik} + c_{ijkl}\right] ], the equation of motion in the initial state is [\rho_{0} u_{i}'' = D_{jkli} {\partial^{2} u_{k} \over \partial X_{j} \partial X_{l}}. \eqno( ]The coefficients [D_{ijkl}] do not present the symmetry of the coefficients [c_{ijkl}] except in the natural state where [D_{ijkl}] and [c_{ijkl}] are equal.

The simplest solutions of the equation of motion are plane waves. We now assume plane sinusoidal waves of the form [u_{i} = A_{i} \exp[i(\omega t - {\bf k\cdot X})], \eqno(]where k is the wavevector.

Substitution of ([link] into ([link] results in [\rho_{0} \omega^{2}A_{j} = D_{ijkl} k_{j}k_{l} A_{k} ]or [\rho_{0} \omega^{2}A_{j} = \Delta_{jk} A_{k} ]with [\Delta_{jk} = D_{ijkl} k_{j}k_{l}].

The quantities [\rho_{0} \omega^{2}A_{j}] and A are, respectively, the eigenvalues and eigenvectors of the matrix [\Delta_{jk}]. Since [\Delta_{jk}] is a real symmetric matrix, the eigenvalues are real and the eigenvectors are orthogonal. Experimental determination of third- and higher-order elastic constants

| top | pdf |

The main experimental procedures for determining the third- and higher-order elastic constants are based on the measurement of stress derivatives of ultrasonic velocities and on harmonic generation experiments. Hydrostatic pressure, which can be accurately measured, has been widely used; however, the measurement of ultrasonic velocities in a solid under hydrostatic pressure cannot lead to the whole set of third-order elastic constants, so uniaxial stress measurements or harmonic generation experiments are then necessary.

In order to interpret wave-propagation measurements in stressed crystals, Thurston (1964)[link] and Brugger (1964)[link] introduced the concept of natural velocity with the following comments:

`According to equation of motion, the wave front is a material plane which has unit normal k in the natural state; a wave front moves from the plane [{\bf k \cdot a = 0}] to the plane [{\bf k \cdot a = L}_{\bf 0}] in the time [L_{0}/W]. Thus W, the natural velocity, is the wave speed referred to natural dimensions for propagation normal to a plane of natural normal k.

In a typical ultrasonic experiment, plane waves are reflected between opposite parallel faces of a specimen, the wave fronts being parallel to these faces. One ordinarily measures a repetition frequency F, which is the inverse of the time required for a round trip between the opposite faces.'

Hence[W = 2L_{0} F.]

In most experiments, the third-order elastic constants and higher-order elastic constants are deduced from the stress derivatives of [\bar{\rho}_{0}W^{2}]. For instance, Table[link] gives the expressions for [\left(\bar{\rho}_{0}W^{2}\right)_{0} ] and [\partial \left(\bar{\rho}_{0}W^{2}\right)_{0} /\partial p ] for a cubic crystal. These quantities refer to the natural state free of stress. In this table, p denotes the hydrostatic pressure and the [\Gamma_{ijkl}]'s are the following linear combinations of third-order elastic constants: [\eqalign{\Gamma_{1111} & = \bar{c}_{111} + 2 \bar{c}_{111}\cr \Gamma_{1122} & = 2\bar{c}_{112} + \bar{c}_{123}\cr \Gamma_{2323} & = \bar{c}_{144} + 2 \bar{c}_{166}.\cr} ]

Table| top | pdf |
Relationships between [\rho W^{2}], its pressure derivatives and the second- and third-order elastic constants

PropagationPolarization[\left(\bar{\rho}_{0}W^{2}\right)_{0} ][\partial \left(\bar{\rho}_{0}W^{2}\right)_{0}/\partial p ]
[100] [100] [\bar{c}_{11}] [-1-\left(2\bar{c}_{11}+\Gamma_{1111}\right)/3\bar{\kappa} ]
[100] [010] [\bar{c}_{44}] [-1-\left(2\bar{c}_{44}+\Gamma_{2323}\right)/3\bar{\kappa} ]
[110] [110] [\left(\bar{c}_{11} + \bar{c}_{12} + 2 \bar{c}_{44}\right)/ 2 ] [-1 -\left(\bar{c}_{11}+\bar{c}_{12}+2\bar{c}_{44} + 0.5\left[\Gamma_{1111}+\Gamma_{1122}+\Gamma_{2323}\right]\right)/3\bar{\kappa} ]
[110] [[1\bar{1}0]] [\left(\bar{c}_{11} - \bar{c}_{12} + \bar{c}_{44}\right)/ 3 ] [-1-\left(\bar{c}_{11}-\bar{c}_{11} + 0.5 \left[\Gamma_{1111} -\Gamma_{1122}\right]\right) / 3 \bar{\kappa} ]
[110] [001] [\bar{c}_{44}] [-1-\left(2\bar{c}_{44}+\Gamma_{2323}\right)/3\bar{\kappa} ]
[111] [111] [\left(\bar{c}_{11} + 2 \bar{c}_{12} + 4 \bar{c}_{44}\right)/ 3 ] [-1-\left(2\bar{c}_{11}+4\bar{c}_{12} + 8\bar{c}_{14} +\left[\Gamma_{1111} +2\Gamma_{1122} + 4\Gamma_{2323}\right]\right) / 9 \bar{\kappa} ]
[111] [[1\bar{1}0]] [\left(\bar{c}_{11} - \bar{c}_{12} + \bar{c}_{44}\right)/ 3 ] [-1-\left(2\bar{c}_{11}-2\bar{c}_{12} + 2\bar{c}_{14} +\left[\Gamma_{1111} -\Gamma_{1122} + \Gamma_{2323}\right]\right) / 9 \bar{\kappa} ]

1.3.8. Glossary

| top | pdf |

[{\bf e}_{i}] covariant basis vector
[A^{T}] transpose of matrix A
[u_{i}] components of the displacement vector
[S_{ij}] components of the strain tensor
[S_{\alpha}] components of the strain Voigt matrix
[T_{ij}] components of the stress tensor
[T_{\alpha}] components of the stress Voigt matrix
p pressure
[\boldnu] normal stress
[\boldtau] shear stress
[s_{ijkl}] second-order elastic compliances
[s_{\alpha\beta}] reduced second-order elastic compliances
[(s_{ijkl})^{\sigma}] adiabatic second-order elastic compliances
[s_{ijklmn}] third-order elastic compliances
[c_{ijkl}] second-order elastic stiffnesses
[(c_{ijkl})^{\sigma}] adiabatic second-order elastic stiffnesses
[c_{\alpha\beta}] reduced second-order elastic stiffnesses
[c_{ijklmn}] third-order elastic stiffnesses
ν Poisson's ratio
E Young's modulus
κ bulk modulus (volume compressibility)
λ, μ Lamé constants
Θ temperature
[c^{S}] specific heat at constant strain
ρ volumic mass
[\Theta_{D}] Debye temperature
[k_{B}] Boltzmann constant
U internal energy
F free energy


Breazeale, M. A. (1984). Determination of third-order elastic constants from ultrasonic harmonic generation. Physical Acoustics, Vol. 17, edited by R. N. Thurston, pp. 2–75. New York: Academic Press.
Brillouin, L. (1932). Propagation des ondes électromagnétiques dans les milieux matériels. Congrès International d'Électricité, Vol. 2, Section 1, pp. 739–788. Paris: Gauthier-Villars.
Brugger, K. (1964). Thermodynamic definition of higher-order elastic coefficients. Phys. Rev. 133, 1611–1612.
De Launay, J. (1956). The theory of specific heats and lattice vibrations. Solid State Physics, Vol. 2, edited by F. Seitz & D. Turnbull, pp. 219–303. New York: Academic Press.
Every, A. G. & McCurdy, A. K. (1992). Second and higher order elastic constants. In Landolt-Börnstein, Group III, Condensed Matter, Volume 29, Subvolume A. Berlin/Heidelberg/New York: Springer.
Fischer, M. (1982). Third- and fourth-order elastic constants of fluoperovskites CsCdF3, TlCdF3, RbCdF3, RbCaF3. J. Phys. Chem. Solids, 43, 673–682.
Fischer, M., Zarembowitch, A. & Breazeale, M. A. (1980). Nonlinear elastic behavior and instabilities in crystals. Ultrasonics Symposium Proc. IEEE, pp. 999–1002.
Fumi, F. G. (1951). Third-order elastic coefficients of crystals. Phys. Rev. 83, 1274–1275.
Fumi, F. G. (1952). Third-order elastic coefficients in trigonal and hexagonal crystals. Phys. Rev. 86, 561.
Fumi, F. G. (1987). Tables for the third-order elastic tensors in crystals. Acta Cryst. A43, 587–588.
Green, R. E. (1973). Treatise on Material Science and Technology, Vol. 3. New York: Academic Press.
McSkimmin, H. J. (1964). Ultrasonic methods for measuring the mechanical properties of solids and fluids. Physical Acoustics, Vol. IA, edited by W. P. Mason, pp. 271–334. New York: Academic Press.
Melcher, R. L. & Scott, B. A. (1972). Soft acoustic modes at the co­operative Jahn–Teller transition in DyVO4. Phys. Rev. Lett. 28, 607–610.
Michard, F., Zarembowitch, A., Vacher, R. & Boyer, L. (1971). Premier son et son zéro dans les nitrates de strontium, barium et plomb. Phonons, edited by M. A. Nusimovici, pp. 321–325. Paris: Flammarion.
Murnaghan, F. D. (1951). Finite Deformation in an Elastic Solid. New York: John Wiley and Sons.
Nouet, J., Zarembowitch, A., Pisarev, R. V., Ferré, J. & Lecomte, M. (1972). Determination of TN for KNiF3 through elastic, magneto-optical and heat capacity measurements. Appl. Phys. Lett. 21, 161–162.
Rousseau, M., Gesland, J. Y., Julliard, J., Nouet, J., Zarembowitch, J. & Zarembowitch, A. (1975). Crystallographic, elastic and Raman scattering investigations of structural phase transitions in RbCdF3 and TlCdF3. Phys. Rev. 12, 1579–1590.
Salje, E. K. H. (1990). Phase Transitions in Ferroelastic and Co-elastic Crystals. Cambridge University Press.
Thurston, R. N. (1964). Wave propagation in fluids and normal solids. Physical Acoustics, Vol. 1A, edited by W. P. Mason, pp. 1–109. New York: Academic Press.
Truesdell, C. & Toupin, R. (1960). The classical field theories. Handbuch der Physik, Vol. III/1, edited by S. Flügge. Berlin, Göttingen, Heidelberg: Springer-Verlag.
Voigt, W. (1910). Lehrbuch der Kristallphysik. 2nd ed. (1929). Leipzig: Teubner. Photoreproduction (1966). New York: Johnson Reprint Corp.
Wallace, D. C. (1970). Thermoelastic theory of stressed crystals and higher-order elastic constants. Solid State Physics, Vol. 25. New York: Academic Press.
Wallace, D. C. (1972). Thermodynamics of Crystals. New York: John Wiley and Sons.
Zarembowitch, A. (1965). Etude théorique et détermination optique des constantes élastiques de monocristaux. Bull. Soc. Fr. Minéral. Cristallogr. 28, 17–49.

to end of page
to top of page