Table of contents

The Golden Quad

On the right side we have external quantities, such as the load FF and deformation δ\delta (or displacement uu). On the left side we have internal quantities, such as the stress σ\sigma and strain ε\varepsilon.

The golden quad provides the three relations:

Equilibrium, which connects loads to stresses. These relations are formulated by mechanics and creating a free body diagram and posing the equilibrium equations.

Deformation relations, these connect deformations to strains and are purely geometric. These are sometimes called kinematic relations. Typically something is known on the boundary, which is introduced as a kinematic relation.

Material relations, or constitutive relations, these connect the strains to the stresses using knowledge about the material which is acquired experimentally. Typically material parameters are Young's modulus EE, Poisson's ratio ν\nu and the thermal expansion coefficient α\alpha.

The workflow using the golden quad is:

  1. Formulate all relations, 1,2,3, according to the figure. This is mechanics, solid mechanics, modeling and computational thinking.

  2. Solve the resulting system of equations. This is not mechanics or solid mechanics, but rudimentary mathematics which is best left to the computer to solve using your favorite technical programming language.

{σ=F/Aσ=Eεδ=Lεδ=F  LE  A\left\lbrace \begin{array}{ll} \sigma =F/A & \\ \sigma =E\varepsilon & \\ \delta =L\varepsilon & \end{array}\right.\Rightarrow \delta =\dfrac{F\;L}{E\;A}

Bodies in series:

δ=δi=NiLiEiAi\delta =\sum \delta_i =\sum \dfrac{N_i L_i }{E_i A_i }

Example 1

A linear elastic rod is fastened between two rigid walls and exposed to heating. Determine the stress and reaction forces on the walls.

Solution

Cut and free body diagram

{NP=0,σ=N/AEquilibrium  equationsεtot=δ/L,δ=0Deformation,δ=0  boundary  conditionσ=EεelastMaterial,elastic  relationεT=αΔTThermal  strainεtot=εelast+εTTotal  strain\left\lbrace \begin{array}{ll} N-P=0,\sigma =N/A & \textrm{Equilibrium}\;\textrm{equations}\\ \varepsilon_{\textrm{tot} } =\delta /L,\delta =0 & \textrm{Deformation},\delta =0\;\textrm{boundary}\;\textrm{condition}\\ \sigma =E\varepsilon_{\textrm{elast} } & \textrm{Material},\textrm{elastic}\;\textrm{relation}\\ \varepsilon_T =\alpha \Delta T & \textrm{Thermal}\;\textrm{strain}\\ \varepsilon_{\textrm{tot} } ={\varepsilon_{\textrm{elast} } +\varepsilon }_T & \textrm{Total}\;\textrm{strain} \end{array}\right.

From this system of equations we get 7 unknowns of which these are the ones we are interested in:

P=E  A  αΔTP=-E\;A\;\alpha \Delta T, σ=E  αΔT\sigma =-E\;\alpha \Delta T

dimension analysis checks out and we clearly see that we have negative force and stress as it should be.

Example 2

Same problem as above but two different rods and materials.

Solution:

{S1+N1=0,N1+N2=0,N2+S2=0Equilibrium  equationsσ1=N1/A1,σ2=N2/A2Equilibrium  equationsεE1=δ1/L1,εE2=δ2/L2,δ1+δ2=0Deformation  and  boundary  conditionσ1=EεE1,σ2=EεE2Material,elastic  relationεT1=α1ΔT,εT2=α2ΔTThermal  strainεtot1=εE1+εT1,εtot2=εE2+εT2Total  strain\left\lbrace \begin{array}{ll} -S_1 +N_1 =0,-N_1 +N_2 =0,-N_2 +S_2 =0 & \textrm{Equilibrium}\;\textrm{equations}\\ \sigma_1 =N_1 /A_1 ,\sigma_2 =N_2 /A_2 & \textrm{Equilibrium}\;\textrm{equations}\\ \varepsilon_{E1} =\delta_1 /L_1 ,\varepsilon_{E2} =\delta_2 /L_2 ,\delta_1 +\delta_2 =0 & \textrm{Deformation}\;\textrm{and}\;\textrm{boundary}\;\textrm{condition}\\ \sigma_1 =E\varepsilon_{E1} ,\sigma_2 =E\varepsilon_{E2} & \textrm{Material},\textrm{elastic}\;\textrm{relation}\\ \varepsilon_{\textrm{T1} } =\alpha_1 \Delta T,\varepsilon_{\textrm{T2} } =\alpha_2 \Delta T & \textrm{Thermal}\;\textrm{strain}\\ \varepsilon_{\textrm{tot1} } =\varepsilon_{\textrm{E1} } +\varepsilon_{\textrm{T1} } ,\varepsilon_{\textrm{tot2} } =\varepsilon_{\textrm{E2} } +\varepsilon_{\textrm{T2} } & \textrm{Total}\;\textrm{strain} \end{array}\right.

We now have 14 equations and 14 unknowns from which we can determine the stress in the two bodies and the reaction forces.