Shear modulus

Shear modulus
Common symbols	G, S, μ
SI unit	Pa
Derivations from other quantities	G = τ / γ = E / [2(1 + ν)]

Shear strain

In solid mechanics, the shear modulus or modulus of rigidity, denoted by G, or sometimes S or μ, is a measure of the elastic shear stiffness of a material and is defined as the ratio of shear stress to the shear strain:

$G\stackrel{\mathrm{d}\mathrm{e}\mathrm{f}}{=}\frac{{\tau }_{xy}}{{\gamma }_{xy}}=\frac{F/A}{\Delta x/l}=\frac{Fl}{A\Delta x}$

where

${\tau }_{xy}=F/A$ = shear stress

$F$ is the force which acts

$A$ is the area on which the force acts

${\gamma }_{xy}$ = shear strain. In engineering $:=\Delta x/l=\tan \theta$, elsewhere $:=\theta$

$\Delta x$ is the transverse displacement

$l$ is the initial length of the area.

The derived SI unit of shear modulus is the pascal (Pa), although it is usually expressed in gigapascals (GPa) or in thousand pounds per square inch (ksi). Its dimensional form is M¹L⁻¹T⁻², replacing force by mass times acceleration.

Explanation

Material	Typical values for shear modulus (GPa) (at room temperature)
Diamond	478.0
Steel	79.3
Iron	52.5
Copper	44.7
Titanium	41.4
Glass	26.2
Aluminium	25.5
Polyethylene	0.117
Rubber	0.0006
Granite	24
Shale	1.6
Limestone	24
Chalk	3.2
Sandstone	0.4
Wood	4

The shear modulus is one of several quantities for measuring the stiffness of materials. All of them arise in the generalized Hooke's law:

• Young's modulus E describes the material's strain response to uniaxial stress in the direction of this stress (like pulling on the ends of a wire or putting a weight on top of a column, with the wire getting longer and the column losing height),
• the Poisson's ratio ν describes the response in the directions orthogonal to this uniaxial stress (the wire getting thinner and the column thicker),
• the bulk modulus K describes the material's response to (uniform) hydrostatic pressure (like the pressure at the bottom of the ocean or a deep swimming pool),
• the shear modulus G describes the material's response to shear stress (like cutting it with dull scissors).

These moduli are not independent, and for isotropic materials they are connected via the equations

$E=2G(1+\nu )=3K(1-2\nu )$

The shear modulus is concerned with the deformation of a solid when it experiences a force perpendicular to one of its surfaces while its opposite face experiences an opposing force (such as friction). In the case of an object shaped like a rectangular prism, it will deform into a parallelepiped. Anisotropic materials such as wood, paper and also essentially all single crystals exhibit differing material response to stress or strain when tested in different directions. In this case, one may need to use the full tensor-expression of the elastic constants, rather than a single scalar value.

One possible definition of a fluid would be a material with zero shear modulus.

Shear waves

File:SpiderGraph ShearModulus.GIF

Influences of selected glass component additions on the shear modulus of a specific base glass.

In homogeneous and isotropic solids, there are two kinds of waves, pressure waves and shear waves. The velocity of a shear wave, $(v_s)$ is controlled by the shear modulus,

$v_s=\sqrt{\frac{G}{\rho }}$

where

G is the shear modulus

$\rho$ is the solid's density.

Shear modulus of metals

Shear modulus of copper as a function of temperature. The experimental data are shown with colored symbols.

The shear modulus of metals is usually observed to decrease with increasing temperature. At high pressures, the shear modulus also appears to increase with the applied pressure. Correlations between the melting temperature, vacancy formation energy, and the shear modulus have been observed in many metals.

Several models exist that attempt to predict the shear modulus of metals (and possibly that of alloys). Shear modulus models that have been used in plastic flow computations include:

1. the Varshni-Chen-Gray model developed by and used in conjunction with the Mechanical Threshold Stress (MTS) plastic flow stress model.
2. the Steinberg-Cochran-Guinan (SCG) shear modulus model developed by and used in conjunction with the Steinberg-Cochran-Guinan-Lund (SCGL) flow stress model.
3. the Nadal and LePoac (NP) shear modulus model that uses Lindemann theory to determine the temperature dependence and the SCG model for pressure dependence of the shear modulus.

Varshni-Chen-Gray model

The Varshni-Chen-Gray model (sometimes referred to as the Varshni equation) has the form:

$\mu (T)={\mu }_0-\frac{D}{\exp (T_0/T)-1}$

where ${\mu }_0$ is the shear modulus at $T=0K$, and $D$ and $T_0$ are material constants.

SCG model

The Steinberg-Cochran-Guinan (SCG) shear modulus model is pressure dependent and has the form

$\mu (p,T)={\mu }_0+\frac{\partial \mu }{\partial p}\frac{p}{{\eta }^{\frac{1}{3}}}+\frac{\partial \mu }{\partial T}(T-300);\eta :=\frac{\rho }{{\rho }_0}$

where, μ₀ is the shear modulus at the reference state (T = 300 K, p = 0, η = 1), p is the pressure, and T is the temperature.

NP model

The Nadal-Le Poac (NP) shear modulus model is a modified version of the SCG model. The empirical temperature dependence of the shear modulus in the SCG model is replaced with an equation based on Lindemann melting theory. The NP shear modulus model has the form:

$\mu (p,T)=\frac{1}{𝒥\left(\hat{T}\right)}[({\mu }_0+\frac{\partial \mu }{\partial p}\frac{p}{{\eta }^{\frac{1}{3}}})(1-\hat{T})+\frac{\rho }{Cm}T];C:=\frac{{\left(6{\pi }^2\right)}^{\frac{2}{3}}}{3}{f}^2$

where

$𝒥(\hat{T}):=1+\exp [-\frac{1+1/\zeta }{1+\zeta /(1-\hat{T})}]\text{for}\hat{T}:=\frac{T}{T_m}\in [0,6+\zeta ],$

and μ₀ is the shear modulus at absolute zero and ambient pressure, ζ is an area, m is the atomic mass, and f is the Lindemann constant.

Shear relaxation modulus

The shear relaxation modulus $G(t)$ is the time-dependent generalization of the shear modulus $G$:

$G=\underset{t\to \mathrm{\infty }}{\lim }G(t)$.

See also

• Elasticity tensor
• Dynamic modulus
• Impulse excitation technique
• Shear strength
• Seismic moment

References

Homogeneous isotropic linear elastic materials have their elastic properties uniquely determined by any two quantities among these; thus, given any two, any other of the elastic moduli can be calculated according to these formulas, provided both for 3D materials (first part of the table) and for 2D materials (second part).

3D Formulae
Knowns	Bulk modulus (K)	Young's modulus (E)	Lamé's first parameter (λ)	Shear modulus (G)	Poisson's ratio (ν)	P-wave modulus (M)	Notes
(K, E)			⁠3K − E/3 − ⁠E/3K⁠⁠	⁠E/3 − ⁠E/3K⁠⁠	⁠1/2⁠ − ⁠E/6K⁠	⁠3K + E/3 − ⁠E/3K⁠⁠
(K, λ)		⁠9K(K − λ)/3K − λ⁠		⁠3(K − λ)/2⁠	⁠λ/3K − λ⁠	3K − 2λ
(K, G)		⁠9KG/3K + G⁠	K − ⁠2G/3⁠		⁠3K − 2G/6K + 2G⁠	K + ⁠4G/3⁠
(K, ν)		3K(1 − 2ν)	⁠3Kν/1 + ν⁠	⁠3K(1 − 2ν)/2(1 + ν)⁠		⁠3K(1 − ν)/1 + ν⁠
(K, M)		⁠9K(M − K)/3K + M⁠	⁠3K − M/2⁠	⁠3(M − K)/4⁠	⁠3K − M/3K + M⁠
(E, λ)	⁠E + 3λ + R/6⁠			⁠E − 3λ + R/4⁠	− ⁠E + R/4λ⁠ − ⁠1/4⁠	⁠E − λ + R/2⁠	R = ±(E2 + 9λ2 + 2Eλ)⁠1/2⁠
(E, G)	⁠EG/3(3G − E)⁠		⁠G(E − 2G)/3G − E⁠		⁠E/2G⁠ − 1	⁠G(4G − E)/3G − E⁠
(E, ν)	⁠E/3 − 6ν⁠		⁠Eν/(1 + ν)(1 − 2ν)⁠	⁠E/2(1 + ν)⁠		⁠E(1 − ν)/(1 + ν)(1 − 2ν)⁠
(E, M)	⁠3M − E + S/6⁠		⁠M − E + S/4⁠	⁠3M + E − S/8⁠	⁠E + S/4M⁠ − ⁠1/4⁠		S = ±(E2 + 9M2 − 10EM)⁠1/2⁠
(λ, G)	λ + ⁠2G/3⁠	⁠G(3λ + 2G)/λ + G⁠			⁠λ/2(λ + G)⁠	λ + 2G
(λ, ν)	⁠λ/3⁠(1 + ⁠1/v⁠)	λ(⁠1/ν⁠ − 2ν − 1)		λ(⁠1/2ν⁠ − 1)		λ(⁠1/ν⁠ − 1)
(λ, M)	⁠M + 2λ/3⁠	⁠(M − λ)(M+2λ)/M + λ⁠		⁠M − λ/2⁠	⁠λ/M + λ⁠
(G, ν)	⁠2G(1 + ν)/3 − 6ν⁠	2G(1 + ν)	⁠2 G ν/1 − 2ν⁠			⁠2G(1 − ν)/1 − 2ν⁠
(G, M)	M − ⁠4G/3⁠	⁠G(3M − 4G)/M − G⁠	M − 2G		⁠M − 2G/2M − 2G⁠
(ν, M)	⁠M(1 + ν)/3(1 − ν)⁠	⁠M(1 + ν)(1 − 2ν)/1 − ν⁠	⁠M ν/1 − ν⁠	⁠M(1 − 2ν)/2(1 − ν)⁠
2D Formulae
Knowns	(K)	(E)	(λ)	(G)	(ν)	(M)	Notes
(K2D, E2D)			⁠2K2D(2K2D − E2D)/4K2D − E2D⁠	⁠K2DE2D/4K2D − E2D⁠	⁠2K2D − E2D/2K2D⁠	⁠4K2D^2/4K2D − E2D⁠
(K2D, λ2D)		⁠4K2D(K2D − λ2D)/2K2D − λ2D⁠		K2D − λ2D	⁠λ2D/2K2D − λ2D⁠	2K2D − λ2D
(K2D, G2D)		⁠4K2DG2D/K2D + G2D⁠	K2D − G2D		⁠K2D − G2D/K2D + G2D⁠	K2D + G2D
(K2D, ν2D)		2K2D(1 − ν2D)	⁠2K2Dν2D/1 + ν2D⁠	⁠K2D(1 − ν2D)/1 + ν2D⁠		⁠2K2D/1 + ν2D⁠
(E2D, G2D)	⁠E2DG2D/4G2D − E2D⁠		⁠2G2D(E2D − 2G2D)/4G2D − E2D⁠		⁠E2D/2G2D⁠ − 1	⁠4G2D^2/4G2D − E2D⁠
(E2D, ν2D)	⁠E2D/2(1 − ν2D)⁠		⁠E2Dν2D/(1 + ν2D)(1 − ν2D)⁠	⁠E2D/2(1 + ν2D)⁠		⁠E2D/(1 + ν2D)(1 − ν2D)⁠
(λ2D, G2D)	λ2D + G2D	⁠4G2D(λ2D + G2D)/λ2D + 2G2D⁠			⁠λ2D/λ2D + 2G2D⁠	λ2D + 2G2D
(λ2D, ν2D)	⁠λ2D(1 + ν2D)/2ν2D⁠	⁠λ2D(1 + ν2D)(1 − ν2D)/ν2D⁠		⁠λ2D(1 − ν2D)/2ν2D⁠		⁠λ2D/ν2D⁠
(G2D, ν2D)	⁠G2D(1 + ν2D)/1 − ν2D⁠	2G2D(1 + ν2D)	⁠2 G2D ν2D/1 − ν2D⁠			⁠2G2D/1 − ν2D⁠
(G2D, M2D)	M2D − G2D	⁠4G2D(M2D − G2D)/M2D⁠	M2D − 2G2D		⁠M2D − 2G2D/M2D⁠