Skip to main content
Engineering LibreTexts

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

    For a rectangular parallelepiped of arbitrary size around interior point cap P, an external body force vector cap B times delta x sub 1, delta x sub 2, and delta x sub 3, is applied. The body force results in each surface having a corresponding surface force vector cap T, with a subscript corresponding to whichever direction it is pointed, multiplied by the delta x values for the corresponding face. Each direction has two faces, one evaluated at initial location x sub 1, x sub 2, or x sub 3, and an equal and opposite force evaluated on the x plus delta x face for each corresponding direction.

    Fig. A.7   Surface forces and a body force acting on a rectangular parallelepiped Δx1Δx2Δx3math-4581.png.

    T1Δx2Δx3x1+Δx1 T1Δx2Δx3x1 + T2Δx1Δx3x2+Δx2 T2Δx1Δx3x2 + T3Δx1Δx2x3+Δx3 T3Δx1Δx2x3 + BΔx1Δx2Δx3. math-4582.png

    For small increments in Δxi, i = 1,2,3math-4583.png, use the Taylor series representation of surface forces results in the equilibrium equation to get eq. (A.78) below.

    x1 T1Δx2Δx3 Δx1 + x2 T2Δx1Δx3 Δx2 + x3 T3Δx1Δx2 Δx3 + BΔx1Δx2Δx3 + O(Δxi)4 = 0. (A.78) math-4584.png

    Arrange the terms in eq. (A.78) to the form

    T1 x1 + T2 x2 + T3 x3 + BΔx1Δx2Δx3 + O(Δxi)4 = 0 (A.79) math-4585.png

    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

    T1 x1 + T2 x2 + T3 x3 + B = 0. (A.80) math-4586.png

    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

    σ11 x1 + σ21 x2 + σ31 x3 + B1 = 0 σ12 x1 + σ22 x2 + σ32 x3 + B2 = 0 σ13 x1 + σ23 x2 + σ33 x3 + B3 = 0 . (A.81) math-4587.png

    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 (σ22Δx1Δx3)x2+Δx2math-4590.png 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 (σ23Δx1Δx3)x2+Δx2math-4591.png 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

    εΔx3(σ22Δx1Δx3)x2+Δx2 + εΔx3(σ22Δx1Δx3)x2 + Δx2 2 (σ23Δx1Δx3)x2+Δx2 + Δx2 2 (σ23Δx1Δx3)x2 + εΔx2(σ33Δx1Δx2)x3+Δx3 εΔ2(σ33Δx1Δx2)x3 Δx3 2 (σ32Δx1Δx2)x3+Δx3 Δx3 2 (σ32Δx1Δx2)x3 = 0. (A.82) math-4592.png

    A slice of a rectangular parallelepiped around point cap P is shown in the x sub 2 x sub 3 plane. Each face has a normal stress sigma sub x x, and shear stress sigma sub x y, with each stress applied on the x face in the y direction. The left and bottom faces correspond to x sub 2 and x sub 3 locations, repsectively, while the right and top faces are x sub 2 plus delta x sub 2 and x sub 3 plus delta x sub 3, respectively. Each stress is multiplied by the width and depth of each face to give the corresponding surface force, such as sigma sub 2 3 on the left face being multiplied by delta x sub 1 and delta x sub 3.

    Fig. A.8   A free body diagram of the parallelepiped at point P for moment equilibrium about the x1math-4588.png-axis. The x1math-4589.png-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

    εΔx3 σ22 x2 Δx2 + O(Δx22)Δx 1Δx3 + Δx2 2 2σ23 + σ23 x2 Δx2 + O(Δx22)Δx 1Δx3 + εΔx2 σ33 x3 Δx3 + O(Δx32)Δx 1Δx2 Δx3 2 2σ32 + σ32 x3 Δx3 + O(Δx32)Δx 1Δx2 = 0. (A.83) math-4593.png

    Expand eq. (A.83) in powers of Δxi to write it as

    σ23 σ32 Δx1Δx2Δx3 + ε σ22 x2 Δx1Δx2Δx32 + σ33 x3 Δx1Δx22Δx 3 + H.O.T. = 0, (A.84) math-4594.png

    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 σ23 σ32 = 0math-4595.png. Moment equilibrium about the x2-axis leads to σ31 σ13 = 0math-4596.png, and moment equilibrium about the x3-axis leads to σ12 σ21 = 0math-4597.png. The equations of moment equilibrium at point P are

    σ12 = σ21σ13 = σ31σ23 = σ32. (A.85) math-4598.png

    Hence, the stress matrix (A.71) is symmetric.

    A.2.2 Transformation of stresses between two Cartesian coordinate systems

    At point P coordinates (x1,x2,x3)math-4599.png are linearly related to coordinates (x1,x2,x3)math-4600.png by eq. (A.41). The stress components σij are functions in the variables (x1,x2,x3)math-4601.png, and the stresses σijmath-4602.png are functions in the variables (x1,x2,x3)math-4603.png. The stress vectors acting on the xi-faces are denoted by Tinline-67.pngi, and those acting on the ximath-4604.png-faces are denoted by Timath-4605.png. These stress vectors are written in their respective coordinate systems by

    T1 T2 T3 = σ11 σ12 σ13 σ21 σ22 σ23 σ31 σ32 σ33 î1 î2 î3 , and T1 T2 T3 = σ11 σ12 σ13 σ21 σ22 σ23 σ31 σ32 σ33 î1 î2 î3 . (A.86) math-4606.png

    In eq. (A.86) the stress matrices are

    [σ] = σ11 σ12 σ13 σ21 σ22 σ23 σ31 σ32 σ33 , and [σ] = σ11 σ12 σ13 σ21 σ22 σ23 σ31 σ32 σ33 . (A.87) math-4607.png

    Note that the stress matrix [σ] is symmetric by eq. (A.85). The stress transformation equations between the Cartesian coordinate systems (x1,x2,x3)math-4608.png and (x1,x2,x3)math-4609.png are determined by selecting the unit normal in eq. (A.77) to be either î1math-4610.png, î2math-4611.png, or î3math-4612.png. First let n̂ = î1math-4613.png such that T(n) = T1math-4614.png in eq. (A.77). From eq. (A.44) we have n̂ = λ11î1 + λ12î2 + λ13î3math-4615.png. Hence, eq. (A.77) becomes

    T1 = λ 11T1 + λ12T2 + λ13T3. (A.88) math-4616.png

    Second, let n̂ = î2math-4617.png such that T(n) = T2math-4618.png. From eq. (A.44) we have n̂ = λ21î1 + λ22î2 + λ23î3math-4619.png. Hence, eq. (A.77) becomes

    T2 = λ 21T1 + λ22T2 + λ23T3. (A.89) math-4620.png

    Third, let n̂ = î3math-4621.png that T(n) = T3math-4622.png and n̂ = λ31î1 + λ32î2 + λ33î3math-4623.png. Hence,

    T3 = λ 31T1 + λ32T2 + λ33T3. (A.90) math-4624.png

    The three selections for the unit normal in eq. (A.77) relate the tractions acting on the ximath-4625.png coordinate faces to the tractions acting on the xi-faces by

    T1 T2 T3 = λ11 λ12 λ13 λ21 λ22 λ23 λ31 λ32 λ33 T1 T2 T3 . (A.91) math-4626.png

    Substitute the expressions for the stress vectors from eq. (A.86) into eq. (A.91) to get

    σ11 σ12 σ13 σ21 σ22 σ23 σ31 σ32 σ33 î1 î2 î3 = λ11 λ12 λ13 λ21 λ22 λ23 λ31 λ32 λ33 σ11 σ12 σ13 σ21 σ22 σ23 σ31 σ32 σ33 î1 î2 î3 . (A.92) math-4627.png

    The inverse eq. (A.44) is

    î1 î2 î3 = λ11 λ21 λ31 λ12 λ22 λ32 λ13 λ23 λ33 î1 î2 î3 . (A.93) math-4628.png

    Substitute eq. (A.93) into the right-hand side of eq. (A.92) and rearrange the result to find

    σ11 σ12 σ13 σ21 σ22 σ23 σ31 σ32 σ33 λ11 λ12 λ13 λ21 λ22 λ23 λ31 λ32 λ33 σ11 σ12 σ13 σ21 σ22 σ23 σ31 σ32 σ33 λ11 λ21 λ31 λ12 λ22 λ32 λ13 λ23 λ33 î1 î2 î3 = 0 0 0 . (A.94) math-4629.png

    To satisfy eq. (A.94) we find that the stress components σijmath-4630.png in the ximath-4631.png-system are related to the stress components σij in the xi-system by

    σ11 σ12 σ13 σ21 σ22 σ23 σ31 σ32 σ33 = λ11 λ12 λ13 λ21 λ22 λ23 λ31 λ32 λ33 σ11 σ12 σ13 σ21 σ22 σ23 σ31 σ32 σ33 λ11 λ21 λ31 λ12 λ22 λ32 λ13 λ23 λ33 . (A.95) math-4632.png

    Equation (A.95) in compact form is

    [σ] = [λ][σ][λ]T. (A.96) math-4633.png

    Pre-multiply eq. (A.96) by [λ]T, post-multiply it by [λ], and note that [λ]T[λ] = [λ][λ]T = [I]math-4634.png to find the inverse transformation

    [σ] = [λ]T[σ][λ]. (A.97) math-4635.png

    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 ximath-4636.png 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 εijmath-4637.png in ximath-4638.png. If

    [εij] = [λ][ε ij][λ]T (A.98) math-4639.png

    then the functions εij and εijmath-4640.png are the components in their respective variables of a second order Cartesian tensor. Similarly, functions σij and σijmath-4641.png 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.


    This page titled 19.2: Equilibrium differential equations is shared under a CC BY-NC-SA 4.0 license and was authored, remixed, and/or curated by Eric Raymond Johnson (Virginia Tech Libraries' Open Education Initiative) via source content that was edited to the style and standards of the LibreTexts platform.