Linear elasticity

From Wikipedia, the free encyclopedia

Linear elasticity is the mathematical study of how solid objects deform and become internally stressed due to prescribed loading conditions. Linear elasticity relies upon the continuum hypothesis and is applicable at macroscopic (and sometimes microscopic) length scales. Linear elasticity is a simplification of the more general nonlinear theory of elasticity and is a branch of continuum mechanics. The fundamental "linearizing" assumptions of linear elasticity are: "small" deformations (or strains) and linear relationships between the components of stress and strain. In addition linear elasticity is only valid for stress states that do not produce yielding. These assumptions are reasonable for many engineering materials and engineering design scenarios. Linear elasticity is therefore used extensively in structural analysis and engineering design, often through the aid of finite element analysis.

1 Mathematical formulation
2 Isotropic homogeneous media
- 2.1 Elastostatics
- 2.2 Elastodynamics - The Wave equation
  - 2.2.1 Plane waves
  - 2.2.2 Elastic wave
3 See also
4 References

[edit] Mathematical formulation

Note that the Einstein summation convention of summing on repeated indices is used below.

Equations governing a linear elastic boundary value problem are based on three tensor equations:

Equation of motion, which is an expression of Newton's second law:

$\ \sigma_{ij,j}+ F_i = \rho \partial_{tt} u_i$

Engineering notation
$\frac{\partial \sigma_x}{\partial x} + \frac{\partial \tau_{yx}}{\partial y} + \frac{\partial \tau_{zx}}{\partial z} + F_x = \rho \frac{\partial^2 u_x}{\partial t^2}$ $\frac{\partial \tau_{xy}}{\partial x} + \frac{\partial \sigma_y}{\partial y} + \frac{\partial \tau_{zy}}{\partial z} + F_y = \rho \frac{\partial^2 u_y}{\partial t^2}$ $\frac{\partial \tau_{xz}}{\partial x} + \frac{\partial \tau_{yz}}{\partial y} + \frac{\partial \sigma_z}{\partial z} + F_z = \rho \frac{\partial^2 u_z}{\partial t^2}$

Engineering notation

$\frac{\partial \sigma_x}{\partial x} + \frac{\partial \tau_{yx}}{\partial y} + \frac{\partial \tau_{zx}}{\partial z} + F_x = \rho \frac{\partial^2 u_x}{\partial t^2}$

$\frac{\partial \tau_{xy}}{\partial x} + \frac{\partial \sigma_y}{\partial y} + \frac{\partial \tau_{zy}}{\partial z} + F_y = \rho \frac{\partial^2 u_y}{\partial t^2}$

$\frac{\partial \tau_{xz}}{\partial x} + \frac{\partial \tau_{yz}}{\partial y} + \frac{\partial \sigma_z}{\partial z} + F_z = \rho \frac{\partial^2 u_z}{\partial t^2}$

These are 3 independent equations with 6 independent unknowns (stresses).

Strain-displacement equations:

$\varepsilon_{ij} =\frac{1}{2} (u_{j,i}+u_{i,j})$

Engineering notation
$\epsilon_x=\frac{\partial u_x}{\partial x}$	$\gamma_{xy}=\frac{\partial u_x}{\partial y}+\frac{\partial u_y}{\partial x}$
$\epsilon_y=\frac{\partial u_y}{\partial y}$	$\gamma_{yz}=\frac{\partial u_y}{\partial z}+\frac{\partial u_z}{\partial y}$
$\epsilon_z=\frac{\partial u_z}{\partial z}$	$\gamma_{zx}=\frac{\partial u_z}{\partial x}+\frac{\partial u_x}{\partial z}$

which are 6 independent equations relating strains and displacements with 9 independent unknowns (strains and displacements)

Constitutive equations. For elastic materials the Hooke's law represents the material behavior and relates the unknown stresses and strains. The general equation for Hooke's law is:

$\sigma_{ij} = C_{ijkl} \, \varepsilon_{kl}$

These are 6 independent equations relating stresses and strains. The coefficients of the elasticity tensor can always be specified so that $\ C_{ijkl}=C_{klij}=C_{jikl}=C_{ijlk}$ .

where:

$\sigma_{ij}=\sigma_{ji} \,$ is the Cauchy stress tensor
$F_i \,$ are the body forces
$\rho \,$ is the mass density
$u_i \,$ is the displacement
$C_{ijkl}\,$ is the elasticity tensor
$\ \lambda, \mu, \nu,$ and $\ E$ are material constants
$\varepsilon_{ij}=\varepsilon_{ji} \,$ is the strain
$\partial_t$ is $\partial/\partial t$ .

An elastostatic boundary value problem for an isotropic-homogeneous media is a system of 15 independent equations and equal number of unknowns (3 equilibrium equations, 6 strain-displacement equations, and 6 constitutive equations). Specifying the boundary conditions, the boundary value problem is completely defined. To solve the system two approaches can be taken according to boundary conditions of the boundary value problem: a displacement formulation, and a stress formulation.

[edit] Isotropic homogeneous media

In isotropic media, the elasticity tensor gives the relationship between the stresses (resulting internal stresses) and the strains (resulting deformations). For an isotropic medium, the elasticity tensor has no preferred direction: an applied force will give the same displacements (relative to the direction of the force) no matter the direction in which the force is applied. In the isotropic case, the elasticity tensor may be written:

$C_{ijkl} = K \, \delta_{ij}\, \delta_{kl} +\mu\, (\delta_{ik}\delta_{jl}+\delta_{il}\delta_{jk}-\textstyle{\frac{2}{3}}\, \delta_{ij}\,\delta_{kl})$

where K is the bulk modulus (or incompressibility), and $μ$ is the shear modulus (or rigidity), two elastic moduli. If the medium is homogeneous as well, then the elastic moduli will not be a function of position in the medium. The constitutive equation may now be written as:

$\ \sigma_{ij} =K\delta_{ij}\varepsilon_{kk}+2\mu(\varepsilon_{ij}-\textstyle{\frac{1}{3}}\delta_{ij}\varepsilon_{kk})$

This expression separates the stress into a scalar part on the left which may be associated with a scalar pressure, and a traceless part on the right which may be associated with shear forces. A simpler expression is:

$\ \sigma_{ij} =\lambda \delta_{ij} \varepsilon_{kk}+2\mu\varepsilon_{ij}$

where λ is Lamé's first parameter. Since the constitutive equation is simply a set of linear equations, the strain may be expressed as a function of the stresses as: (Sommerfeld 1964)

$\varepsilon_{ij} = \frac{1}{9K}\delta_{ij}\sigma_{kk} + \frac{1}{2\mu}\left(\sigma_{ij}-\textstyle{\frac{1}{3}}\delta_{ij}\sigma_{kk}\right)$

which is again, a scalar part on the left and a traceless shear part on the right. More simply:

$\varepsilon_{ij} =\frac{1}{2\mu}\sigma_{ij}-\frac{\nu}{E}\delta_{ij}\sigma_{kk}$

where ν is Poisson's ratio and E is Young's modulus.

[edit] Elastostatics

Elastostatics is the study of linear elasticity under the conditions of equilibrium, in which all forces on the elastic body sum to zero, and the displacements are not a function of time. The equilibrium equations are then

$\ \sigma_{ij,j}+ F_i = 0$

Engineering notation
$\frac{\partial \sigma_x}{\partial x} + \frac{\partial \tau_{yx}}{\partial y} + \frac{\partial \tau_{zx}}{\partial z} + F_x = 0$ $\frac{\partial \tau_{xy}}{\partial x} + \frac{\partial \sigma_y}{\partial y} + \frac{\partial \tau_{zy}}{\partial z} + F_y = 0$ $\frac{\partial \tau_{xz}}{\partial x} + \frac{\partial \tau_{yz}}{\partial y} + \frac{\partial \sigma_z}{\partial z} + F_z = 0$

Engineering notation

$\frac{\partial \sigma_x}{\partial x} + \frac{\partial \tau_{yx}}{\partial y} + \frac{\partial \tau_{zx}}{\partial z} + F_x = 0$

$\frac{\partial \tau_{xy}}{\partial x} + \frac{\partial \sigma_y}{\partial y} + \frac{\partial \tau_{zy}}{\partial z} + F_y = 0$

$\frac{\partial \tau_{xz}}{\partial x} + \frac{\partial \tau_{yz}}{\partial y} + \frac{\partial \sigma_z}{\partial z} + F_z = 0$

This section will discuss only the isotropic homogeneous case.

[edit] Displacement formulation

In this case, the displacements are prescribed everywhere in the boundary. In this approach, the strains and stresses are eliminated from the formulation, leaving the displacements as the unknowns to be solved for in the governing equations. First, the strain-displacement equations are substituted into the constitutive equations (Hooke's Law), eliminating the strains as unknowns:

$\ \begin{align} \sigma_{ij} &= \lambda \delta_{ij} \varepsilon_{kk}+2\mu\varepsilon_{ij} \\ &= \lambda\delta_{ij}u_{k,k}+\mu\left(u_{i,j}+u_{j,i}\right) \\ \end{align}$

Differentiating yields:

$\ \sigma_{ij,j} = \lambda u_{k,ki}+\mu\left(u_{i,jj}+u_{j,ij}\right)$

Substituting into the equilibrium equation yields:

$\ \lambda u_{k,ki}+\mu\left(u_{i,jj}+u_{j,ij}\right) +F_i=0$

$\ \mu u_{i,jj}+(\mu+\lambda)u_{j,ij}+F_i=0$

where $\ \lambda$ and $\ \mu$ are Lamé parameters. In this way, the only unknowns left are the displacements, hence the name for this formulation. The governing equations obtained in this manner are called Navier-Cauchy equations or, alternatively, the elastostatic equations.

Derivation of Navier-Cauchy equations in Engineering notation

First, the $\ x$ -direction will be considered. Substituting the strain-displacement equations into the equilibrium equation in the $\ x$ -direction we have

$\sigma_x = 2\mu\varepsilon_x + \lambda(\varepsilon_x + \varepsilon_y +\varepsilon_z)=2\mu\frac{\partial u_x}{\partial x}+\lambda\left(\frac{\partial u_x}{\partial x}+\frac{\partial u_y}{\partial y}+\frac{\partial u_z}{\partial z}\right)$

$\ \tau_{xy} = \mu\gamma_{xy}=\mu\left(\frac{\partial u_x}{\partial y}+\frac{\partial u_y}{\partial x}\right)$

$\ \tau_{xz} = \mu\gamma_{zx}=\mu\left(\frac{\partial u_z}{\partial x}+\frac{\partial u_x}{\partial z}\right)$

Then substituting these equations into the equilibrium equation in the $\ x$ -direction we have

$\frac{\partial \sigma_x}{\partial x} + \frac{\partial \tau_{yx}}{\partial y} + \frac{\partial \tau_{zx}}{\partial z} + F_x = 0$

$\ \frac{\partial}{\partial x}\left( 2\mu\frac{\partial u_x}{\partial x}+\lambda\left(\frac{\partial u_x}{\partial x}+\frac{\partial u_y}{\partial y}+\frac{\partial u_z}{\partial z}\right)\right) + \mu\frac{\partial}{\partial y}\left(\frac{\partial u_x}{\partial y}+\frac{\partial u_y}{\partial x}\right)+ \mu\frac{\partial}{\partial z}\left(\frac{\partial u_z}{\partial x}+\frac{\partial u_x}{\partial z}\right)+F_x=0\,$

Solving and rearranging we have

$\left(\lambda+\mu\right)\frac{\partial}{\partial x}\left(\frac{\partial u_x}{\partial x}+\frac{\partial u_y}{\partial y}+\frac{\partial u_z}{\partial z}\right)+\mu\left(\frac{\partial^2 u_x}{\partial x^2}+\frac{\partial^2 u_x}{\partial y^2}+\frac{\partial^2 u_x}{\partial z^2}\right)+F_x=0$

Following the same procedure for the $\ y$ -direction and $\ z$ -direction we have

$\left(\lambda+\mu\right)\frac{\partial}{\partial y}\left(\frac{\partial u_x}{\partial x}+\frac{\partial u_y}{\partial y}+\frac{\partial u_z}{\partial z}\right)+\mu\left(\frac{\partial^2 u_y}{\partial x^2}+\frac{\partial^2 u_y}{\partial y^2}+\frac{\partial^2 u_y}{\partial z^2}\right)+F_y=0$

$\left(\lambda+\mu\right)\frac{\partial}{\partial z}\left(\frac{\partial u_x}{\partial x}+\frac{\partial u_y}{\partial y}+\frac{\partial u_z}{\partial z}\right)+\mu\left(\frac{\partial^2 u_z}{\partial x^2}+\frac{\partial^2 u_z}{\partial y^2}+\frac{\partial^2 u_z}{\partial z^2}\right)+F_z=0$

These last 3 equations are the Navier-Cauchy equations, which can be also expressed in vector notation as

$(\lambda+\mu)\nabla(\nabla\cdot\mathbf{u})+\mu\nabla^2\mathbf{u}+\mathbf{F}=0$

Once the displacement field has been calculated, the displacements can be replaced into the strain-displacement equations to solve for strains, which later are used in the constitutive equations to solve for stresses.

[edit] The biharmonic equation

The elastostatic equation may be written:

$\ (\alpha^2-\beta^2)u_{j,ij}+ \beta^2u_{i,mm}=-F_i$

Taking the divergence of both sides of the elastostatic equation and assuming the force has zero divergence ( $\ F_{i,i}=0$ ) we have

$\ (\alpha^2-\beta^2)u_{j,iij}+\beta^2u_{i,imm} = 0$

Noting that summed indices need not match, and that the partial derivatives commute, the two differential terms are seen to be the same and we have:

$\ \alpha^2u_{j,iij} = 0$

from which we conclude that:

$\ u_{j,iij} = 0$

Taking the Laplacian of both sides of the elastostatic equation, and assuming in addition $\ F_{i,kk}=0$ , we have

$\ (\alpha^2-\beta^2)u_{j,kkij}+\beta^2u_{i,kkmm}=0$

From the divergence equation, the first term on the right is zero (Note: again, the summed indices need not match) and we have:

$\ \beta^2u_{i,kkmm}=0$

from which we conclude that:

$\ u_{i,kkmm}=0$

or, in coordinate free notation $\nabla^4 \mathbf{u}=0$ which is just the biharmonic equation in $\mathbf{u}$ .

[edit] Stress formulation

In this case, the surface tractions are prescribed everywhere on the surface boundary. In this approach, the strains and displacements are eliminated leaving the stresses as the unknowns to be solved for in the governing equations. Once the stress field is found, the strains are then found using the constitutive equations.

There are six independent components of the stress tensor which need to be determined, yet in the displacement formulation, there are only three components of the displacement vector which need to be determined. This means that there are some constraints which must be placed upon the stress tensor, to reduce the number of degrees of freedom to three. Using the constitutive equations, these constraints are derived directly from corresponding constraints which must hold for the strain tensor, which also has six independent components. The constraints on the strain tensor are derivable directly from the definition of the strain tensor as a function of the displacement vector field, which means that these constraints introduce no new concepts or information. It is the constraints on the strain tensor that are most easily understood. If the elastic medium is visualized as a set of infinitesimal cubes in the unstrained state, after the medium is strained, an arbitrary strain tensor may not yield a situation in which the distorted cubes still fit together without overlapping. The constraints on the strain tensor that are required to assure that this is the case were discovered by Saint Venant, and are called the "Saint Venant compatibility equations". These are 81 equations, 6 of which are independent non-trivial equations, which relate the different strain components. These are expressed in index notation as:

$\ \varepsilon_{ij,km}+\varepsilon_{km,ij}-\varepsilon_{ik,jm}-\varepsilon_{jm,ik}=0$

Engineering notation
$\ \frac{\partial^2 \epsilon_x}{\partial y^2} + \frac{\partial^2 \epsilon_y}{\partial x^2} = 2 \frac{\partial^2 \epsilon_{xy}}{\partial x \partial y}$ $\ \frac{\partial^2 \epsilon_y}{\partial z^2} + \frac{\partial^2 \epsilon_z}{\partial y^2} = 2 \frac{\partial^2 \epsilon_{yz}}{\partial y \partial z}$ $\ \frac{\partial^2 \epsilon_x}{\partial z^2} + \frac{\partial^2 \epsilon_z}{\partial x^2} = 2 \frac{\partial^2 \epsilon_{zx}}{\partial z \partial x}$ $\ \frac{\partial^2 \epsilon_x}{\partial y \partial z} = \frac{\partial}{\partial x} \left ( -\frac{\partial \epsilon_{yz}}{\partial x} + \frac{\partial \epsilon_{zx}}{\partial y} + \frac{\partial \epsilon_{xy}}{\partial z}\right)$ $\ \frac{\partial^2 \epsilon_y}{\partial z \partial x} = \frac{\partial}{\partial y} \left ( \frac{\partial \epsilon_{yz}}{\partial x} - \frac{\partial \epsilon_{zx}}{\partial y} + \frac{\partial \epsilon_{xy}}{\partial z}\right)$ $\ \frac{\partial^2 \epsilon_z}{\partial x \partial y} = \frac{\partial}{\partial z} \left ( \frac{\partial \epsilon_{yz}}{\partial x} + \frac{\partial \epsilon_{zx}}{\partial y} - \frac{\partial \epsilon_{xy}}{\partial z}\right)$

Engineering notation

$\ \frac{\partial^2 \epsilon_x}{\partial y^2} + \frac{\partial^2 \epsilon_y}{\partial x^2} = 2 \frac{\partial^2 \epsilon_{xy}}{\partial x \partial y}$

$\ \frac{\partial^2 \epsilon_y}{\partial z^2} + \frac{\partial^2 \epsilon_z}{\partial y^2} = 2 \frac{\partial^2 \epsilon_{yz}}{\partial y \partial z}$

$\ \frac{\partial^2 \epsilon_x}{\partial z^2} + \frac{\partial^2 \epsilon_z}{\partial x^2} = 2 \frac{\partial^2 \epsilon_{zx}}{\partial z \partial x}$

$\ \frac{\partial^2 \epsilon_x}{\partial y \partial z} = \frac{\partial}{\partial x} \left ( -\frac{\partial \epsilon_{yz}}{\partial x} + \frac{\partial \epsilon_{zx}}{\partial y} + \frac{\partial \epsilon_{xy}}{\partial z}\right)$

$\ \frac{\partial^2 \epsilon_y}{\partial z \partial x} = \frac{\partial}{\partial y} \left ( \frac{\partial \epsilon_{yz}}{\partial x} - \frac{\partial \epsilon_{zx}}{\partial y} + \frac{\partial \epsilon_{xy}}{\partial z}\right)$

$\ \frac{\partial^2 \epsilon_z}{\partial x \partial y} = \frac{\partial}{\partial z} \left ( \frac{\partial \epsilon_{yz}}{\partial x} + \frac{\partial \epsilon_{zx}}{\partial y} - \frac{\partial \epsilon_{xy}}{\partial z}\right)$

The strains in this equation are then expressed in terms of the stresses using the constitutive equations, which yields the corresponding constraints on the stress tensor. These constraints on the stress tensor are known as the Beltrami-Michell equations of compatibility:

$\ \sigma_{ij,kk}+\frac{1}{1+\nu}\sigma_{kk,ij}+F_{i,j}+F_{j,i}+\frac{\nu}{1-\nu}\delta_{i,j}F_{k,k}=0$

These constraints, along with the equilibrium equation (or equation of motion for elastodynamics) allow the calculation of the stress tensor field. Once the stress field has been calculated from these equations, the strains can be obtained from the constitutive equations, and the displacement field from the strain-displacement equations.

An alternative solution technique is to express the stress tensor in terms of stress functions which automatically yield a solution to the equilibrium equation. The stress functions then obey a single differential equation which corresponds to the compatibility equations.

[edit] Solutions for elastostatic cases

Thomson's solution - point force in an infinite isotropic medium

The most important solution of the Navier-Cauchy or elastostatic equation is for that of a force acting at a point in an infinite isotropic medium. This solution was found by William Thomson (later Lord Kelvin) in 1848 (Thomson 1848). This solution is the analog of Coulomb's law in electrostatics. A derivation is given in (Landau & Lifshitz § 8). Defining

$a=1-2\nu\,$

$b=2(1-\nu)=a+1\,$

where $ν$ is Poisson's ratio, the solution may be expressed as

$u_i=G_{ik}F_k\,$

where $F k$ is the force vector being applied at the point, and $G i k$ is a tensor Green's function which may be written in Cartesian coordinates as:

$G_{ik}= \frac{1}{4\pi\mu r}\left[ \left(1-\frac{1}{2b}\right)\delta_{ik}+\frac{1}{2b}\frac{x_i x_k}{r^3} \right]$

It may be also compactly written as:

$G_{ik}= \frac{1}{4\pi\mu}\left[\frac{\delta_{ik}}{r}-\frac{1}{2b}\frac{\partial r^2}{\partial x_i\partial x_k}\right]$

and it may be explicitly written as:

$G_{ik}=\frac{1}{4\pi\mu r}\begin{bmatrix} 1-\frac{1}{2b}+\frac{1}{2b}\frac{x^2}{r^2} & \frac{1}{2b}\frac{xy} {r^2} & \frac{1}{2b}\frac{xz} {r^2} \\ \frac{1}{2b}\frac{yx} {r^2} & 1-\frac{1}{2b}+\frac{1}{2b}\frac{y^2}{r^2} & \frac{1}{2b}\frac{yz} {r^2} \\ \frac{1}{2b}\frac{zx} {r^2} & \frac{1}{2b}\frac{zy} {r^2} & 1-\frac{1}{2b}+\frac{1}{2b}\frac{z^2}{r^2} \end{bmatrix}$

In cylindrical coordinates ( $ρ,φ, z$ ) it may be written as:

$G_{ik}=\frac{1}{4\pi \mu r}\begin{bmatrix} 1-\frac{1}{2b}\frac{z^2}{r^2}&0&\frac{1}{2b}\frac{\rho z}{r^2}\\ 0&1-\frac{1}{2b}&0\\ \frac{1}{2b}\frac{z \rho}{r^2}&0&1-\frac{1}{2b}\frac{\rho^2}{r^2} \end{bmatrix}$

It is particularly helpful to write the displacement in cylindrical coordinates for a point force $F z$ directed along the z-axis. Defining $\hat{\mathbf{\rho}}$ and $\hat{\mathbf{z}}$ as unit vectors in the $ρ$ and $z$ directions respectively yields:

$\mathbf{u}=\frac{F_z}{4\pi\mu r}\left[\frac{1}{4(1-\nu)}\,\frac{\rho z}{r^2}\hat{\mathbf{\rho}} + \left(1-\frac{1}{4(1-\nu)}\,\frac{\rho^2}{r^2}\right)\hat{\mathbf{z}}\right]$

It can be seen that there is a component of the displacement in the direction of the force, which diminishes, as is the case for the potential in electrostatics, as 1/r for large r. There is also an additional ρ-directed component.

Boussinesq-Cerruti solution - point force at the origin of an infinite isotropic half-space

Another useful solution is that of a point force acting on the surface of an infinite half-space. It was derived by Boussinesq(Boussinesq 1885) and a derivation is given in (See Landau & Lifshitz § 8). In this case, the solution is again written as a Green's tensor which goes to zero at infinity, and the component of the stress tensor normal to the surface vanishes. This solution may be written as in Cartesian coordinates as:

$G_{ik}=\frac{1}{4\pi\mu}\begin{bmatrix} \frac{b}{r}+\frac{x^2}{r^3}-\frac{ax^2}{r(r+z)^2}+\frac{az}{r(r+z)} & \frac{xy}{r^3}-\frac{axy}{r(r+z)^2}& \frac{xz}{r^3}-\frac{ax}{r(r+z)}\\ \frac{yx}{r^3} -\frac{ayx}{r(r+z)^2}& \frac{b}{r}+\frac{y^2}{r^3}-\frac{ay^2}{r(r+z)^2}+\frac{az}{r(r+z)} & \frac{yz}{r^3} -\frac{ay}{r(r+z)}\\ \frac{zx}{r^3}+\frac{ax}{r(r+z)}& \frac{zy}{r^3}+\frac{ay}{r(r+z)}& \frac{b}{r}+\frac{z^2}{r^3} \end{bmatrix}$

[edit] Elastodynamics - The Wave equation

Elastodynamics is the study of linear elasticity which include variation in time. The most common case considered in elastodynamics is the wave equation. This section will discuss only the isotropic homogeneous case.

If the material is homogeneous (i.e. the elasticity tensor is constant throughout the material), the three basic equations can be combined to form the elastodynamic equation:

$\mu u_{i,jj}+(\mu+\lambda)u_{j,ij}+F_i=\rho\partial_{tt}u_i \,\,\,\,\,\,\,\,\mathrm{or}\,\,\,\,\,\,\,\, \mu\nabla^2\mathbf{u}+(\mu+\lambda)\nabla(\nabla\cdot\mathbf{u})+\mathbf{F}=\rho\frac{\partial^2\mathbf{u}}{\partial t^2}$

From the elastodynamic equation one gets the wave equation

$(\delta_{kl} \partial_{tt}-A_{kl}[\nabla]) \, u_l = \frac{1}{\rho} f_k$

where

$A_{kl}[\nabla]=\frac{1}{\rho} \, \partial_i \, C_{iklj} \, \partial_j$

is the acoustic differential operator, and $δ k l$ is Kronecker delta.

In isotropic media, the elasticity tensor has the form

$C_{ijkl} = K \, \delta_{ij}\, \delta_{kl} +\mu\, (\delta_{ik}\delta_{jl}+\delta_{il}\delta_{jk}-\frac{2}{3}\, \delta_{ij}\,\delta_{kl})$

where $K$ is the bulk modulus (or incompressibility), and $μ$ is the shear modulus (or rigidity), two elastic moduli. If the material is homogeneous (i.e. the elasticity tensor is constant throughout the material), the acoustic operator becomes:

$A_{ij}[\nabla]=\alpha^2 \partial_i\partial_j+\beta^2(\partial_m\partial_m\delta_{ij}-\partial_i\partial_j)\,$

and the acoustic algebraic operator becomes

$A_{ij}[\mathbf{k}]=\alpha^2 k_ik_j+\beta^2(k_mk_m\delta_{ij}-k_ik_j)\,$

where

$\alpha^2=\left(K+\frac{4}{3}\mu\right)/\rho \qquad \beta^2=\mu/\rho$

are the eigenvalues of $A[\hat{\mathbf{k}}]$ with eigenvectors $\hat{\mathbf{u}}$ parallel and orthogonal to the propagation direction $\hat{\mathbf{k}}$ , respectively. In the seismological literature, the corresponding plane waves are called P-waves and S-waves (see Seismic wave).

[edit] Plane waves

A plane wave has the form

$\mathbf{u}[\mathbf{x}, \, t] = U[\mathbf{k} \cdot \mathbf{x} - \omega \, t] \, \hat{\mathbf{u}}$

with $\hat{\mathbf{u}}$ of unit length. It is a solution of the wave equation with zero forcing, if and only if $ω 2$ and $\hat{\mathbf{u}}$ constitute an eigenvalue/eigenvector pair of the acoustic algebraic operator

$A_{kl}[\mathbf{k}]=\frac{1}{\rho} \, k_i \, C_{iklj} \, k_j$

This propagation condition may be written as

$A[\hat{\mathbf{k}}] \, \hat{\mathbf{u}}=c^2 \, \hat{\mathbf{u}}$

where $\hat{\mathbf{k}} = \mathbf{k} / \sqrt{\mathbf{k}\cdot\mathbf{k}}$ denotes propagation direction and $c=\omega/\sqrt{\mathbf{k}\cdot\mathbf{k}}$ is phase velocity.

[edit] Elastic wave

An elastic wave is a type of mechanical wave that propagates in elastic or viscoelastic materials. The elasticity of the material provides the restoring force of the wave. When they occur in the Earth as the result of an earthquake or other disturbance, elastic waves are usually called seismic waves.

Solids are often assumed to exhibit linear elasticity. That is, Hooke's law is assumed to be valid. The restoring force is directly proportional to the displacement (or strain), but is oppositely directed.

[edit] See also

[edit] References

Gurtin M. E., Introduction to Continuum Mechanics, Academic Press 1981
Elastostatics (Kip Thorne)

Landau, L.D.; Lifshitz, E. M. (1986). Theory of Elasticity, 3rd Edition, Oxford, England: Butterworth Heinemann. ISBN 0-7506-2633-X.

Boussinesq, Joseph (1885). Application des potentiels à l'étude de l'équilibre et du mouvement des solides élastiques. Paris, France: Gauthier-Villars.

Mindlin, R. D. (1936). "Force at a point in the interior of a semi-infinite solid". Physics 7: 195-202.

Hertz, Heinrich (1882). "Contact between solid elastic bodies". Journ. f'tir reine iind angewandte Math. 92.

Sommerfeld, Arnold (1964). Mechanics of Deformable Bodies. New York: Academic Press.

Categories: Elasticity (physics)

Linear elasticity

From Wikipedia, the free encyclopedia

Contents

[edit] Mathematical formulation

[edit] Isotropic homogeneous media

[edit] Elastostatics

[edit] Displacement formulation

[edit] The biharmonic equation

[edit] Stress formulation

[edit] Solutions for elastostatic cases

[edit] Elastodynamics - The Wave equation

[edit] Plane waves

[edit] Elastic wave

[edit] See also

[edit] References

Views

Navigation

Interaction

Search

Languages