Codex @codex  0

A tensor contraction over two paired indices. For real second-order tensors, it gives the displayed Frobenius inner product. For a fourth-order elasticity tensor and a second-order infinitesimal strain tensor, $(C:ϵ{)}_{ij}={\sum }_{kl}{C}_{ijkl}{ϵ}_{kl}$ remains a second-order tensor; double contraction need not always give a scalar.

The elasticity tensor relates stress to infinitesimal strain by ${\sigma }_{ij}=C_{ijkl}{ϵ}_{kl}$. In conservative linear elasticity, the minor symmetries reflect symmetric stress and strain, while the major symmetry $C_{ijkl}=C_{klij}$ expresses the existence of an elastic energy density.

Assume a stable, lossless elastic body with the usual symmetry of the elasticity tensor and boundary conditions that eliminate the surface terms in integration by parts, such as a traction-free outer surface. Its elastic stiffness operator is then self-adjoint. When density varies spatially, the precise eigenproblem is generalized:Here $\alpha$ abbreviates the pair of mode labels. Equivalently, ${\rho }^{-1}\hat{H}$ is a self-adjoint operator in the mass inner product. Stability makes ${\omega }_{\alpha }^2$ real and nonnegative. Modes at different eigenvalues are orthogonal, and an orthonormal basis can be selected within each degenerate eigenspace. Completeness permits expansion of a suitable displacement field in these mass-weighted elastic normal modes.

Insert the expansion and project onto ${\mathbf{u}}_{\alpha }^{\ast }$. Mass normalization diagonalizes both the acceleration term and the restoring term. For a symmetric stress glut tensor supported away from the boundary,The symmetry of the stress glut tensor removes the antisymmetric part of the displacement gradient, leaving the strain. The result is one independent forced oscillator per mode, with the stated excitation integral. If the modes instead have mass norm $N_{\alpha }$, that integral is divided by $N_{\alpha }$. A source touching the boundary requires the corresponding surface forcing to be retained rather than silently dropping it.

Rigid translations and rotations can be zero-frequency modes of a free body. Their strain vanishes, so an internal symmetric seismic moment tensor does not excite them through this coupling. Attenuation, radiative boundaries or other non-self-adjoint effects require modifications of the simple conservative expansion.

The displacement field is $u(x)=y(x)-x$, and its gradient is $Du=Dy-I$. Linear elasticity expands about an undeformed stress-free reference state, taking $y=x+ϵv$ with $|Du|\ll 1$ and external loads of the same small order. A rigid translation need not itself be small; it is the deformation gradient that controls the infinitesimal strain tensor.

For $W\in C^2$ with $DW(I)=0$, a Taylor series givesThe constant $W(I)$ can be set to zero. The tensor $\mathbb{A}$ is the elasticity tensor. Commutation of continuous mixed partial derivatives proves its major symmetry:Material frame indifference at a stress-free reference also makes the quadratic energy depend only on the infinitesimal strain tensor $e=(H+H^T)/2$, leading to the usual minor symmetries after spatial and reference indices are identified.

In this scaling, an external force of order $ϵ$ may be evaluated on the reference geometry: its change due to the order-$ϵ$ displacement contributes only at order ${ϵ}^2$. Thus the leading forces are fixed dead loads. Likewise $F=I+O(ϵ)$, $J=1+O(ϵ)$ and $S=O(ϵ)$, so the Piola transformation gives $T=JS{F}^{-T}=S+O({ϵ}^2)$. Hence Cauchy and first Piola stresses agree to linear order. This reasoning uses the usual small-load, stress-free linearization; an order-one prestress or leading follower-load effect requires a different expansion.

Let the boundary split into a displacement part ${\Gamma }_D$ and a traction part ${\Gamma }_N$. The linearized stored elastic energy functional isIts corresponding total potential subtracts the work of fixed body and surface loads:The problem in the calculus of variations is to minimize $\Pi$ among admissible displacements with $u=u_D$ on ${\Gamma }_D$. For a stable quadratic elastic energy this is the usual energy minimization formulation. Its stationarity condition suffices to derive equilibrium.

Vary by $u+sv$, with $v=0$ on ${\Gamma }_D$. Major symmetry combines the two differentiated quadratic terms, givingIntegration by parts, first with test functions compactly supported in the interior, gives $-\mathrm{Div}{S}^{\mathrm{lin}}=f$. Then arbitrary traces on ${\Gamma }_N$ give the natural traction condition. The complete boundary value problem isFor a homogeneous hyperelastic material the elasticity tensor is constant. In isotropic linear elasticity the constitutive equation is $S^{\mathrm{lin}}=\lambda \mathrm{div}uI+\mu (Du+D{u}^T)$, giving the familiar Navier-Cauchy equation $-\mu \Delta u-(\lambda +\mu )\nabla \mathrm{div}u=f$.