19.2: Equilibrium differential equations
- Page ID
- 95345
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)
( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\id}{\mathrm{id}}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\kernel}{\mathrm{null}\,}\)
\( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\)
\( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\)
\( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)
\( \newcommand{\vectorA}[1]{\vec{#1}} % arrow\)
\( \newcommand{\vectorAt}[1]{\vec{\text{#1}}} % arrow\)
\( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vectorC}[1]{\textbf{#1}} \)
\( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)
\( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)
\( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\(\newcommand{\avec}{\mathbf a}\) \(\newcommand{\bvec}{\mathbf b}\) \(\newcommand{\cvec}{\mathbf c}\) \(\newcommand{\dvec}{\mathbf d}\) \(\newcommand{\dtil}{\widetilde{\mathbf d}}\) \(\newcommand{\evec}{\mathbf e}\) \(\newcommand{\fvec}{\mathbf f}\) \(\newcommand{\nvec}{\mathbf n}\) \(\newcommand{\pvec}{\mathbf p}\) \(\newcommand{\qvec}{\mathbf q}\) \(\newcommand{\svec}{\mathbf s}\) \(\newcommand{\tvec}{\mathbf t}\) \(\newcommand{\uvec}{\mathbf u}\) \(\newcommand{\vvec}{\mathbf v}\) \(\newcommand{\wvec}{\mathbf w}\) \(\newcommand{\xvec}{\mathbf x}\) \(\newcommand{\yvec}{\mathbf y}\) \(\newcommand{\zvec}{\mathbf z}\) \(\newcommand{\rvec}{\mathbf r}\) \(\newcommand{\mvec}{\mathbf m}\) \(\newcommand{\zerovec}{\mathbf 0}\) \(\newcommand{\onevec}{\mathbf 1}\) \(\newcommand{\real}{\mathbb R}\) \(\newcommand{\twovec}[2]{\left[\begin{array}{r}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\ctwovec}[2]{\left[\begin{array}{c}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\threevec}[3]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\cthreevec}[3]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\fourvec}[4]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\cfourvec}[4]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\fivevec}[5]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\cfivevec}[5]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\mattwo}[4]{\left[\begin{array}{rr}#1 \amp #2 \\ #3 \amp #4 \\ \end{array}\right]}\) \(\newcommand{\laspan}[1]{\text{Span}\{#1\}}\) \(\newcommand{\bcal}{\cal B}\) \(\newcommand{\ccal}{\cal C}\) \(\newcommand{\scal}{\cal S}\) \(\newcommand{\wcal}{\cal W}\) \(\newcommand{\ecal}{\cal E}\) \(\newcommand{\coords}[2]{\left\{#1\right\}_{#2}}\) \(\newcommand{\gray}[1]{\color{gray}{#1}}\) \(\newcommand{\lgray}[1]{\color{lightgray}{#1}}\) \(\newcommand{\rank}{\operatorname{rank}}\) \(\newcommand{\row}{\text{Row}}\) \(\newcommand{\col}{\text{Col}}\) \(\renewcommand{\row}{\text{Row}}\) \(\newcommand{\nul}{\text{Nul}}\) \(\newcommand{\var}{\text{Var}}\) \(\newcommand{\corr}{\text{corr}}\) \(\newcommand{\len}[1]{\left|#1\right|}\) \(\newcommand{\bbar}{\overline{\bvec}}\) \(\newcommand{\bhat}{\widehat{\bvec}}\) \(\newcommand{\bperp}{\bvec^\perp}\) \(\newcommand{\xhat}{\widehat{\xvec}}\) \(\newcommand{\vhat}{\widehat{\vvec}}\) \(\newcommand{\uhat}{\widehat{\uvec}}\) \(\newcommand{\what}{\widehat{\wvec}}\) \(\newcommand{\Sighat}{\widehat{\Sigma}}\) \(\newcommand{\lt}{<}\) \(\newcommand{\gt}{>}\) \(\newcommand{\amp}{&}\) \(\definecolor{fillinmathshade}{gray}{0.9}\)A.2.1 Equilibrium differential equations
Consider the forces acting on a rectangular parallelepiped at point P. The free body diagram is shown in figure A.7. The vector sum of forces is

Fig. A.7 Surface forces and a body force acting on a rectangular parallelepiped .
For small increments in Δxi, , use the Taylor series representation of surface forces results in the equilibrium equation to get eq. (A.78) below.
Arrange the terms in eq. (A.78) to the form
Divide eq. (A.79) by the volume followed by the limit as Δx1Δx2Δx3→0 to get the vector differential equation of force equilibrium at point P as
Substitute eq. (A.70) for the traction vectors in eq. (A.80) to write the equilibrium differential equations in the xi-coordinate directions. In the order of x1, x2, x3 coordinate directions these equations are
Now consider moment equilibrium about the coordinate axes of the rectangular parallelepiped at point P. For moment equilibrium about the x1-axis refer to the free body diagram in figure A.8. The moment arm from point P to the line of action of the normal force acting on the positive x2-face is denoted by εΔx3, where ε is a small numerical value. Parameter ε is not known, but this will not matter in the end result. The moment arm from point P to the line of action of the shear force
acting on the positive x2 face is Δx2∕2. Including all the forces shown in figure A.8, the sum of moments about the x1-axis through point P, counterclockwise positive, is

Fig. A.8 A free body diagram of the parallelepiped at point P for moment equilibrium about the -axis. The
-axis points normal to the page toward the reader.
Use the Taylor series to expand the forces acting on the positive coordinate faces with respect to the forces acting on the negative coordinate faces to get
Expand eq. (A.83) in powers of Δxi to write it as
where H.O.T. means higher order terms, that is, terms of quartic powers and higher in the increments in the coordinates. Notice the terms multiplied by ε are quartic powers of the increments in the coordinates. Division of eq. (A.84) by Δx1Δx2Δx3, followed by the limit of Δxi→0 leads the condition of moment equilibrium about the x1-axis that . Moment equilibrium about the x2-axis leads to
, and moment equilibrium about the x3-axis leads to
. The equations of moment equilibrium at point P are
Hence, the stress matrix (A.71) is symmetric.
A.2.2 Transformation of stresses between two Cartesian coordinate systems
At point P coordinates are linearly related to coordinates
by eq. (A.41). The stress components σij are functions in the variables
, and the stresses
are functions in the variables
. The stress vectors acting on the xi-faces are denoted by
i, and those acting on the
-faces are denoted by
. These stress vectors are written in their respective coordinate systems by
In eq. (A.86) the stress matrices are
Note that the stress matrix [σ] is symmetric by eq. (A.85). The stress transformation equations between the Cartesian coordinate systems and
are determined by selecting the unit normal in eq. (A.77) to be either
,
, or
. First let
such that
in eq. (A.77). From eq. (A.44) we have
. Hence, eq. (A.77) becomes
Second, let such that
. From eq. (A.44) we have
. Hence, eq. (A.77) becomes
Third, let that
and
. Hence,
The three selections for the unit normal in eq. (A.77) relate the tractions acting on the coordinate faces to the tractions acting on the xi-faces by
Substitute the expressions for the stress vectors from eq. (A.86) into eq. (A.91) to get
The inverse eq. (A.44) is
Substitute eq. (A.93) into the right-hand side of eq. (A.92) and rearrange the result to find
To satisfy eq. (A.94) we find that the stress components in the
-system are related to the stress components σij in the xi-system by
Equation (A.95) in compact form is
Pre-multiply eq. (A.96) by [λ]T, post-multiply it by [λ], and note that to find the inverse transformation
The transpose of eq. (A.96) is [λ][σ]T[λ]T, but [σ]T = [σ], so the stress matrix [σ′] is also symmetric. Comparing the strain transformation eq. (A.63) to the stress transformation eq. (A.96), it is clear that the transformation of strains εij is the same form as the transformation of the stresses σij.
A.2.3 Cartesian tensors
A tensor is a system of numbers or functions, whose components obey a certain law of transformation when the independent variables undergo a linear transformation. If the independent variables are the rectangular cartesian systems xi and transforming by the linear relations given by {x′} = [λ]{x} at point P, then the systems obeying certain laws of transformation are called Cartesian tensors.
Definition. A system of order two may be defined to have nine components εij in xi and nine components in
. If
then the functions εij and are the components in their respective variables of a second order Cartesian tensor. Similarly, functions σij and
are the components in their respective variables of a second order Cartesian tensor.
A.3 Linear elastic material law
To this point in the study of the mechanics of a solid body we have eighteen unknown functions of the Cartesian coordinates x1, x2, and x3. These are the three displacements u1, u2, and u3, the six strains ε11, ε22, ε33, γ23, γ31, and γ12, and nine stresses σ11, σ12, σ13, σ21, σ22, σ23, σ31, σ32, and σ33. There are twelve equations relating these unknowns; the six strain-displacement equations (A.35) and (A.36), and the six equilibrium equations (A.81) and (A.85). Therefore we need six more equations to get the number of unknowns equal to number of equations. The additional six equations come from the relations between the strains and the stresses, which express the material law. This relation between strains and stresses for different materials is established by material characterization tests on standard test specimens.
Solid bodies that can instantly recover their original size and shape when the forces producing the deformation are removed are called perfectly elastic. The elastic limit is defined as the greatest stress that can be applied without resulting in permanent strain on release of the stress. Elasticity is applicable to any body provided the stresses do not exceed the elastic limit. For many solid bodies there is a region where the stress is very nearly proportional to strain. The proportional limit is defined as the greatest stress for which the stress is still proportional to the strain. Both the elastic limit and proportional limit cannot be precisely determined from test data since they are defined by the limiting cases of no permanent deformation and no deviation from linearity. In practice the definition of the yield strength of a material is used to determine the limit of elastic behavior. See article 4.2 on page 77 for a discussion on yield criteria.