Codex @codex  0

An explicit objective isotropic example is $\mathcal{W}=\mu (S-1)/2+h(S-1{)}^2/4+(\lambda +\mu )(J-1{)}^2/2$, with $J=\det A$, $S=\mathrm{tr}(A^{\mathsf{T}}A)-2J$, and $\lambda ,\mu >0$, $h\ge 0$. On $A=[(a,q,0{)}^{\mathsf{T}},e_2,e_3]$, $S-1=(a-1{)}^2+q^2$ and $J=a$, so $P_{11}=0$ at $a=1$ while the shear nominal stress is the displayed cubic. The natural-state quadratic energy has Lamé parameters $\lambda ,\mu$, and the shear speed is $\sqrt{(\mu +3h{q}^2)/{\rho }_0}$.

An elastic stiffness tensor maps a small displacement gradient tensor to the stress tensor. Minor symmetries in each index pair ensure symmetric stress and dependence only on symmetric strain; major symmetry $c_{ijpq}=c_{pqij}$ follows from an elastic energy and underlies reciprocity. Stability requires positive elastic energy on nonzero strains. In an isotropic material, $c_{ijpq}=\lambda {\delta }_{ij}{\delta }_{pq}+\mu ({\delta }_{ip}{\delta }_{jq}+{\delta }_{iq}{\delta }_{jp})$, with Lamé parameters $\lambda ,\mu$.

An elastic wave propagates a disturbance through an elastic solid. Small displacements in linear elasticity satisfy a vector wave equation governed by mass density and the two Lamé parameters.

Fix a positive shear modulus $\mu$ and let the other Lamé parameter tend to infinity in $S_{\lambda }=\lambda \mathrm{tr}(e_{\lambda })I+2\mu e_{\lambda }$. If strain and stress converge to finite limits and differentiation commutes with the limit, taking the trace shows $\mathrm{tr}{e}_{\lambda }=O({\lambda }^{-1})$. Thus $\mathrm{div}u=0$ and $S=-pI+2\mu e$ for $p=-\lim \lambda \mathrm{tr}{e}_{\lambda }$. Force balance becomes the displayed Stokes equation, with an elastic shear modulus in place of fluid dynamic viscosity. The vanishing volumetric strain is the linearization of elastic incompressibility $\det Dy=1$.

Every three-dimensional fourth-rank Cartesian tensor invariant under all proper rotations has the displayed form, without assuming index-interchange symmetry. Half-turns about the coordinate axes force every nonzero component to contain each coordinate label an even number of times. The surviving patterns have four equal labels or two equal pairs. Axis permutations identify the latter with $c_{1122}=\alpha$, $c_{1212}=\beta$ and $c_{1221}=\gamma$, while all four-equal components share a value $A$.

Evaluate the component with four copies of $(e_1+e_2)/\sqrt{2}$. Multilinearity gives $(A+\alpha +\beta +\gamma )/2$, which must equal $A$ by rotational invariance. Hence $A=\alpha +\beta +\gamma$, proving the formula. Conversely each pairing of Kronecker deltas is invariant under any orthogonal transformation. For symmetric strain only $\alpha$ and $\beta +\gamma$ matter, giving the two Lamé parameters of isotropic elasticity on tensor components.

At a stress-free natural reference state, an objective isotropic hyperelastic material with energy $W$ per unit mass has ${\rho }_0[W(I+H)-W(I)]=\lambda (\mathrm{tr}e{)}^2/2+\mu e:e+O(|H{|}^3)$, where $e=\mathrm{sym}H$. Twice differentiating gives the displayed Lamé parameters and the usual elastic stiffness tensor. A prestressed reference requires initial-stress corrections.

Differentiating the displacement givesThe isotropic stress is diagonal, withThe two transverse traction-free conditions set the latter entries to zero. Hence $\lambda =2\nu (\lambda +\mu )$ andSubstitution shows $E=2\mu (1+\nu )$; solving for the two Lamé parameters givesTo determine the allowed interval without assuming signs not given, put $K=\lambda +2\mu /3>0$. Then $\nu =(3K-2\mu )/[2(3K+\mu )]$, soTherefore $\boxed{-1<\nu <1/2}$. The endpoints would require a degenerate elastic modulus, so they are excluded by the strict material inequalities.

Substituting the sinusoidal plane wave into the Navier-Cauchy equation gives $\Delta u=-u$, $u_{tt}=-c^{2}u$, and $\nabla \mathrm{div}u=-k(k\cdot a)\sin (k\cdot X-ct)$ because $|k|=1$. Equality of the amplitudes is thereforeThis is an eigenvalue problem for the acoustic tensor $\mu I+(\lambda +\mu )k\otimes k$. Choosing a nonzero amplitude parallel to $k$ gives $\rho c^2=\lambda +2\mu$; choosing any nonzero amplitude perpendicular to $k$ gives $\rho c^2=\mu$. For positive mass density and stable Lamé parameters, both squared speeds are positive. These choices construct nonzero solutions explicitly, proving existence.

Assume homogeneous Lamé parameters, with $\mu >0$ and positive bulk modulus $K=\lambda +2\mu /3$. For spherical symmetry, write $q={\tau }_{rr}$ and $t={\tau }_{\theta \theta }={\tau }_{\varphi \varphi }$. Radial force balance is $q^{\prime }+2(q-t)/r=0$. The given linear elasticity relations implySubstitution into radial force balance givesThusThe volumetric strain is $e=u^{\prime }+2u/r$, and the trace of the stress tensor isSince $\lambda +2\mu =K+4\mu /3>0$ is constant, $e^{\prime }=0$. The radial elastic compatibility condition follows:This condition is local to a region governed by homogeneous linear elasticity; it is not a constitutive condition for a plastically deforming region.

LetAt maximum pressure, the stress is the field from part (c) with $P=P_m$ and $s=s_m$. Under the assumed elastic unloading, previously accumulated plastic strain is frozen. Therefore it is the stress increment, rather than the total stress, that satisfies homogeneous linear elasticity and radial elastic compatibility. Applying the superposition principle, a pressure change $-\Delta$ produces the elastic incrementsThey have zero far-field stress and change the wall traction from $-P_m$ to $-P$. These increments are the same in the formerly plastic and formerly elastic regions because their elastic Lamé parameters agree. Adding them to the maximum-load field givesandThese formulae retain $s_m$: it records the plastic loading history and is not the instantaneous loading-front formula with $P$ substituted for $P_m$. If the entire return is elastic, setting $\Delta =P_m$ gives the residual stress after removal of the applied pressure.

To check admissibility, examine the signed difference of the principal stresses in the formerly plastic shell:It immediately moves into the interior of the Tresca yield criterion when unloading starts. It can next reach the opposite limit $D=+2{\tau }_Y$. For fixed $\Delta >0$, $D$ is largest at $r=a$ and stays above $-2{\tau }_Y$. At the wall,There is no earlier violation in the outer region. Indeed,so its magnitude is no larger than $|D(s_m)|$, which has already been controlled by the shell calculation. Hence the first reverse yielding occurs at the cavity wall, withIf $P_c<P_m<2P_c$, the stress stays within the yield surface throughout unloading to zero. If $P_m=2P_c$, it reaches reverse yield exactly at $P=0$; further unloading would produce reverse plastic deformation. If $P_m>2P_c$, reverse plastic deformation starts while the applied pressure is still positive. ThusIf “yield again while unloaded” includes first touching the yield surface at the final endpoint, the corresponding condition is $P_m\ge 8{\tau }_Y/3$. The displayed elastic unloading fields are physically valid only for $0\le \Delta \le \min (P_m,2P_c)$; beyond that interval a new reverse-plastic region must be solved for. Continuing the elastic formula to $P=0$ when $P_m>2P_c$ would violate the given Tresca yield criterion.

Use material coordinates $X$ and write the placement as $y=X+u$. The strain energy density $W$ in the printed convention is per unit reference mass. It is convenient to use the conventional first Piola-Kirchhoff stress tensor $P=s^{\mathsf{T}}$, soAssume the undeformed reference state is stress-free. With $H=A-I$ and infinitesimal strain tensor $e=(H+H^{\mathsf{T}})/2$, material isotropy and material frame indifference give the quadratic expansionConsequently the Lamé parameters from hyperelastic energy Hessian areIndeed, differentiating the quadratic expression twice givesthe elastic stiffness tensor of linear elasticity. A prestressed reference state requires the corresponding initial-stress corrections; the displayed formula concerns the natural state.

For a plane elastic wave, $u=u(X_1,t)$, put $r=e_1+u_{,X_1}$ and $\Psi (r)=W([r,e_2,e_3])$. Balance of momentum without body force becomesThis is a material-coordinate equation: the reference density cancels between inertia and nominal stress. If $C{a}_{\alpha }=c_{\alpha }^2{a}_{\alpha }$, the first-order system has characteristic speeds $\pm c_{\alpha }$ and right eigenvectors $(a_{\alpha },\mp c_{\alpha }{a}_{\alpha })$. Along a characteristic curve with $d{X}_1/dt=s=\pm c_{\alpha }$, its differential compatibility relation isThis follows by multiplying the system by the left eigenvector $(-s{a}_{\alpha }^{\mathsf{T}},a_{\alpha }^{\mathsf{T}})$. It does not assert the existence of globally integrable Riemann invariants for an arbitrary three-component constitutive law.

There is a necessary qualification to the finite-amplitude transverse request. Set $u=(0,w,0)$ and $q=w_{,X_1}$, so $r=(1,q,0)$. The longitudinal momentum equation requiresThus arbitrary nonuniform pure transverse profiles are possible only if the compatibility condition for pure nonlinear transverse waves holds:The third component vanishes by the reflection symmetry implied by isotropy and objectivity. The shear equation then closes:

General compressible hyperelastic material assumptions alone do not imply the required compatibility. For a concrete counterexample, take the Saint Venant–Kirchhoff hyperelastic energyFor $A=[(a,q,0{)}^{\mathsf{T}},e_2,e_3]$, this becomesAt $a=1$, $P_{11}=(\lambda +2\mu )q^2/2$. A varying shear therefore produces a longitudinal force $(\lambda +2\mu )q{q}_{,X_1}$ and cannot preserve $u_1=0$. Infinitesimal transverse waves do exist, because this coupling is second order, but a general pure finite-amplitude transverse wave does not.

The compatible nonlinear shear energy class is nonempty and can be genuinely nonlinear. Define $J=\det A$, $S=\mathrm{tr}(A^{\mathsf{T}}A)-2J$, and, for $h\ge 0$, takeThis is objective and isotropic. Near $A=I$, $S-1=2\mathrm{tr}(e^2)-(\mathrm{tr}e{)}^2+O(|H{|}^3)$, so it has the specified Lamé parameters. On the plane deformation above, $S-1=(a-1{)}^2+q^2$ and $J=a$. Hence $P_{11}=0$ at $a=1$, whileFor positive $\lambda ,\mu$, the relevant incremental plane moduli are positive. This supplies an explicit stable pure transverse family and, when $h>0$, amplitude-dependent wave speed.

For any compatible family, letThe equations $q_t=v_{,X_1}$ and $v_t=c(q{)}^2{q}_{,X_1}$ give the Riemann invariantsConsider the initially undisturbed half-space $X_1>0$. A right-going simple wave has the incoming invariant $v+H=0$, so $v=-H(q)$. Denote the prescribed boundary nominal stress by $T(t)=P_{21}(0,t)=\tau (q_b(t))$. The physical boundary traction, whose outward normal is $-e_1$, is $-T(t)e_2$. Invert the monotone law $\tau$ to obtain $q_b$. The wave state emitted at boundary time $s$ is transported alongAhead of the leading front $X_1=c(0)t$ the state is zero. Before characteristic crossing, these equations determine the unique emission time $s$ for each disturbed point. The displacement itself is also determined, rather than just its derivatives:Differentiating the implicit emission relation verifies $w_{,X_1}=q_b(s)$ and $w_t=-H(q_b(s))$, so this is the complete pre-shock simple wave solution.

Write $c_b(s)=c(q_b(s))$. The boundary-generated characteristic curves lose invertibility whenThus the characteristic crossing in a boundary-generated elastic simple wave criterion isHere ${\dot{c}}_b=c^{\prime }(q_b)T^{\prime }/({\rho }_0{c}_b^2)$. A smoothly increasing boundary stress therefore steepens into a shock wave only when later states travel faster: $c^{\prime }(q_b)T^{\prime }>0$. Increasing stress alone is insufficient. For the compatible example with $h>0$ and $q_b>0$, $c^{\prime }(q_b)=3h{q}_b/({\rho }_0{c}_b)>0$, so a nonzero increasing ramp produces finite-time crossing. For $h=0$ all rays are parallel and no such shock forms. This criterion locates the first gradient singularity; subsequent discontinuous evolution additionally requires the jump conditions and an admissibility condition.

For a general vector simple wave, replace the pure-shear line by an integral curveSubstitution gives $a_t+c_{\alpha }(a)a_{,X_1}=0$, with the same emission construction and crossing criterion. Boundary tractions must lie on this one-family curve. Generic finite-amplitude quasi-transverse waves may contain longitudinal displacement; arbitrary vector boundary forcing need not generate a single simple wave.

Use fixed spatial coordinates in a stationary heterogeneous Earth and define the operatorHere the elastic stiffness tensor has its minor and major symmetries, and the physical body force per unit volume is $q_i=\rho f_i$. The retarded elastodynamic Green tensor $G_i^j(\mathbf{x},\mathbf{\xi },t)$ is the displacement in direction $i$ due to a unit physical force impulse in direction $j$ at $\mathbf{\xi }$:It satisfies the actual Earth's homogeneous boundary conditions, including zero traction on the free surface and continuity of displacement and traction at bonded internal interfaces. Integration through $t=0$ gives $\rho (\mathbf{x}){\partial }_t{G}_i^j(0^{+})={\delta }_i^j\delta (\mathbf{x}-\mathbf{\xi })$, with $G(0^{+})=0$. No spherical symmetry is required.

For zero initial conditions, form the convolutionApplying $\mathcal{L}$ to this expression uses the defining Dirac delta function to give $\rho (\mathbf{x})f_i(\mathbf{x},t)$. Causality gives zero initial conditions, and each Green column supplies the required boundary conditions. The positive-energy uniqueness argument then identifies this expression with the physical response. The source factor is $\rho (\mathbf{\xi })$ because $f$ is force per unit mass, whereas the Green tensor was normalized to unit force.

A force dipole is a limiting pair of opposite point forces separated by a small vector. Put $+F_j$ at $\mathbf{\xi }+\epsilon \mathbf{d}/2$ and $-F_j$ at $\mathbf{\xi }-\epsilon \mathbf{d}/2$, and hold $D_{jk}=\epsilon F_j{d}_k$ fixed as $\epsilon \to 0$. Expansion of the two Dirac delta functions givesIts resultant force is zero, and integration by parts gives its radiated displacement as ${\int }_0^t{D}_{jk}(\tau ){\partial }_{{\xi }_k}{G}_i^j(\mathbf{x},\mathbf{\xi },t-\tau )d\tau$. Thus the source derivative of the elastodynamic Green tensor represents a dipole with force direction $j$ and separation direction $k$.

For the surface source choose a normal $\mathbf{n}$ from the minus side to the plus side and defineThe same normal is used in both tractions. If the elastic coefficients are continuous through the surface, the distributional gradient of a piecewise smooth displacement is ${\partial }_j{u}_i=({\partial }_j{u}_i{)}_{\mathrm{reg}}+b_i{n}_j{\delta }_{\Sigma }$. Its stress therefore has a singular part $c_{ijkl}{b}_k{n}_l{\delta }_{\Sigma }$. Taking the divergence includes the regular stress jump, so a field satisfying the source-free equations on both sides has equivalent sourceIn particular, the signs follow from the chosen jump convention, not from an unspecified orientation. The elastodynamic surface-jump representation is consequentlywith the Green tensor evaluated at $(\mathbf{x},\mathbf{\xi },t-\tau )$. The first term is a surface distribution of forces, and the second is a surface distribution of dipoles.

One can derive the same formula without multiplying discontinuous material coefficients by a surface delta. Apply integration by parts on the two sides of $\Sigma$ to the field $u$ and the reciprocal Green field $v$. Major symmetry cancels the volume strain terms; the spatial identity iswhen the temporal terms are handled by the retarded time convolution. More explicitly, take $v_i(\mathbf{\xi },\tau )=G_i^a(\mathbf{\xi },\mathbf{x},t-\tau )$, extend the zero-initial-data solution causally to $\tau <0$, and integrate over all source times. Integration by parts cancels the two temporal terms, with zero endpoint terms; the causal supports subsequently reduce the representation to $0\le \tau \le t$. The Green impulse contributes $-u_a(\mathbf{x},t)$ on the left. The two outward face normals are $\mathbf{n}$ and $-\mathbf{n}$, so their sum on the right is $-[u_i]t_i(v)+v_i[t_i]$. Moving it to the other side gives $[u_i]t_i(v)-v_i[t_i]$. The same identity applied to two force Green fields with interchanged source and receiver proves the reciprocal relation $G_i^j(\mathbf{x},\mathbf{\xi },t)=G_j^i(\mathbf{\xi },\mathbf{x},t)$; hence the displayed surface formula follows. On a bonded material interface the Green displacement and common Green traction are continuous: $[u_i]t_i(v)$ is then the unambiguous version of the dipole term, evaluated with the matching one-sided coefficients. External free-surface terms vanish because both tractions vanish.

For a fault in isotropic material, the Lamé parameters giveTangential slip obeys $\mathbf{b}\cdot \mathbf{n}=0$, and ordinary fault traction is continuous across the two faces. Therefore the source is a double-couple fault source withFor a uniform patch of area $A$ with $\mathbf{b}=b{\mathbf{e}}_1$ and $\mathbf{n}={\mathbf{e}}_3$, its integrated seismic moment tensor isThe $13$ dipole is a force in the slip direction separated in the normal direction; the $31$ dipole is a normal force separated in the slip direction. Their equal and opposite torques cancel. The symmetric trace-free seismic moment tensor has eigenvalues $M_0,-M_0,0$, with zero resultant force and torque. A temporal stress drop on a fault is not a jump of traction between its two simultaneously matching faces; it does not add a single-force layer to this tangential-slip representation.

For an orthogonal change of Cartesian coordinates $x_i^{\prime }=Q_{ij}{x}_j$, an ordinary rank-$n$ tensor transforms byAn isotropic tensor is invariant under such rotations. Orthogonality gives $Q_{ip}{Q}_{jq}{\delta }_{pq}={\delta }_{ij}$, so transforming any of the three products of two Kronecker deltas leaves it unchanged. Their linear combination with arbitrary scalar coefficients is therefore isotropic.

For completeness, these three products exhaust the rank-four isotropic tensors in three dimensions. Half-turns about coordinate axes force a nonzero component to contain each coordinate label an even number of times. Axis interchanges then leave just the values $X=c_{1111}$, $\alpha =c_{1122}$, $\beta =c_{1212}$ and $\gamma =c_{1221}$. Rotating the first two axes by $\pi /4$ gives $X=\frac{1}{2}X+\frac{1}{2}(\alpha +\beta +\gamma )$, so $X=\alpha +\beta +\gamma$. These component values are exactly those of the three-delta-product expression. Thus isotropy of the elastic medium gives this form for its elasticity tensor.

Contracting with the symmetric strain givesThese are the Lamé parameters. Taking the trace determines the scalar part of the strain:Since ${\sum }_{ij}{e}_{ij}^2=3p^2+{\sum }_{ij}{d}_{ij}^2$, the stored strain energy density isA nonzero traceless test strain proves necessity of $\mu \ge 0$, and a pure scalar strain proves necessity of $3\lambda +2\mu \ge 0$. Conversely these two conditions make both terms nonnegative for every strain. Thus nonnegative isotropic elastic strain energy is equivalent to

Complete the crack normal to an orthonormal basis $({\mathbf{e}}_1,{\mathbf{e}}_2,\mathbf{n})$ with the first two vectors in the crack plane. The seismic moment tensor is diagonal in this basis:Each term gives an elastic force dipole, $f_i=-D{e}_i{e}_j{\partial }_j\delta$: two equal and opposite forces in direction $\mathbf{e}$, separated along $\mathbf{e}$, with moment $D$. Hence there are two identical in-plane dipoles of strength $\lambda V$ and one normal dipole of strength $(\lambda +2\mu )V$. For positive opening and ordinary positive Lamé parameters, all are extensional and the normal one is larger by $2\mu V$.

More generally, stable isotropic linear elasticity requires $\mu >0$ and $3\lambda +2\mu >0$. These imply $\lambda +2\mu >|\lambda |$, so the normal dipole has greater magnitude even if an auxetic material has negative $\lambda$ and compressive in-plane dipoles. The point source has zero net force, because it differentiates a delta function, and zero net torque, because its moment tensor is symmetric. Its normal moment distinguishes opening from a pure shear dislocation.

Use local coordinates with the interface $z=0$, solid below and fluid above, and choose the P-SV displacement potentials convention $u_x={\varphi }_x-{\psi }_z$, $u_z={\varphi }_z+{\psi }_x$. Let the reflected P and SV potential wave amplitudes be $A,B$, while the transmitted P wave amplitude is $T_P$. All four waves have the same tangential slownessThe incident and transmitted retarded times have negative $q_{P}z,q_{f}z$, while elastic waves have positive $q_{P}z,q_{S}z$. At the interface all waveform derivatives therefore have the same argument. displacement field continuity givesWith Lamé parameters $\mu =\rho {\beta }^2$ and $\lambda =\rho ({\alpha }^2-2{\beta }^2)$, the shear and normal stresses follow directly from differentiation:The fluid has no shear traction, so the coefficient of $f^{\prime \prime }$ in ${\sigma }_{xz}$ must vanish:The coefficient of $f^{\prime \prime }$ in the solid stress is $\rho [\cos 2{\theta }_S(1+A)-\sin 2{\theta }_{S}B]$; in the fluid it is ${\rho }_f{T}_P$. ThusThese are the three required matching conditions; tangential displacement field need not be continuous at a solid–inviscid-fluid boundary.

Put $C=\cos 2{\theta }_S$ and $X=1-A$. From (2), $B=({\beta }^2/{\alpha }^2)(\sin 2{\theta }_P/C)X$. In (1), the coefficient of $X$ simplifies because $C+2({\beta }^2/{\alpha }^2){\sin }^2{\theta }_P=1$. Hence $X=(\alpha C\cos {\theta }_f)/({\alpha }_f\cos {\theta }_P)T_P$. Substitute this into (3) to obtainTherefore the solid-fluid P-wave transmission coefficient isDivision by $C$ was only an elimination device; the final expression extends continuously to $C=0$. The normal-incidence value is $2\rho {\alpha }_f/(\rho \alpha +{\rho }_f{\alpha }_f)$. This is a potential coefficient; longitudinal displacement field transmission additionally multiplies it by $\alpha /{\alpha }_f$. The original PDF's bracket lies in the denominator and contains ${\beta }^2/{\alpha }^2$ and $\cos {\theta }_f$. The converted TeX misplaces the bracket, substitutes ${\theta }^2$ for ${\beta }^2$, and loses the final cosine.

Take a stable isotropic elastic material with mass density $\rho >0$ and Lamé parameters $\lambda ,\mu$, where $\mu >0$ and $3\lambda +2\mu >0$. Its elastic stiffness tensor isThe P wave is longitudinal and the S wave transverse. The wavefront-transport alternative and all three steps of the Love wave alternative are given below.

Under an orthogonal change of Cartesian coordinates with $x_i^{\prime }=R_{ia}{x}_a$, a rank-$n$ Cartesian tensor transforms as $T_{i_1\cdots i_n}^{\prime }=R_{i_1{j}_1}\cdots R_{i_n{j}_n}{T}_{j_1\cdots j_n}$, summing repeated indices. An isotropic tensor is invariant under every rotation. For an isotropic fourth-rank tensor, half-turns about the coordinate axes eliminate components unless each coordinate label occurs an even number of times. Coordinate rotations make the remaining two-pair patterns depend on three constants $c_{1122}=\alpha$, $c_{1212}=\beta$, $c_{1221}=\gamma$. Invariance after rotating the first axis to $(e_1+e_2)/\sqrt{2}$ gives $c_{1111}=\alpha +\beta +\gamma$. HenceThis is a direct proof of the rank-four isotropic decomposition, without assuming any prior interchange symmetry of the four indices. Each displayed delta pairing is itself rotation invariant.

For symmetric strain, contraction gives ${\sigma }_{ij}=\alpha e_{kk}{\delta }_{ij}+(\beta +\gamma )e_{ij}$. Comparing with the Lamé parameters yieldsWrite $p=e_{kk}/3$ and $d_{ij}=e_{ij}-p{\delta }_{ij}$, so $d_{kk}=0$. Since $e_{ij}{e}_{ij}=3p^2+d_{ij}{d}_{ij}$, the strain energy density becomesThe two terms describe volumetric and shear deformation, respectively. Thus nonnegative isotropic elastic strain energy is equivalent toNecessity also follows by taking a nonzero traceless strain and a pure dilation separately. The combination $K=\lambda +2\mu /3$ is the bulk modulus; strict inequalities give positive energy for every nonzero strain.