12. Finite Element equations and the Direct Stiffness method
Finite Element Equations & Stiffness Matrices
πŸ“‚ Please log in or enroll to access resources
πŸ“Œ Element stiffness matrix...by a different route

If you found the logic of how I developed the element stiffness matrix in this lecture hard to follow, or it just simply didn't make much sense to you, here's an alternative approach that get's you to the same place.

Let’s say our bar is pinned at node 2, so no displacement at end 2. When we apply a force F1F_1 to end 1, we observe a displacement of Ξ΄x1\delta_{x1}. The force-displacement relationship at end one is therefore,

F1=EALβ€…Ξ΄x1(1)F_1 = \frac{EA}{L}\:\delta_{x1} \tag{1}

Now, since the bar is in equilibrium, we know the reaction force at end 2, F2F_2, is the equal to the negative of F1F_1,

F2=βˆ’F1F_2 = -F_1

Therefore, we can say,

F2=βˆ’EALβ€…Ξ΄x1(2)F_2 = -\frac{EA}{L}\:\delta_{x1} \tag{2}

With equations 1 and 2 put to one side, let’s turn our attention to the other end of the element, node 2. Now imagine node one is pinned and a force is applied at node 2. We can run through the exact same steps to construct equations for F1F_1 and F2F_2 in this scenario too.

The force-displacement relationship at node 2 is,

F2=EALβ€…Ξ΄x2(3)F_2 = \frac{EA}{L}\:\delta_{x2} \tag{3}

And again, since the bar is in equilibrium, we can say F1=βˆ’F2F_1 = -F_2, giving us,

F1=βˆ’F2=βˆ’EALβ€…Ξ΄x2(4)F_1 = -F_2 = -\frac{EA}{L}\:\delta_{x2} \tag{4}

Finally, in the case the neither end is pinned and both ends of the bar undergo axial displacements Ξ΄x1\delta_{x1} and Ξ΄x2\delta_{x2}, we can obtain the corresponding forces at each end using superposition of equations 1 and 4 for F1F_1 and 2 and 3 for F2F_2,

F1=EAL(Ξ΄x1βˆ’Ξ΄x2)(5)F_1 = \frac{EA}{L}(\delta_{x1} - \delta_{x2}) \tag{5}
F2=βˆ’EAL(Ξ΄x1+Ξ΄x2)(6)F_2 = -\frac{EA}{L}(\delta_{x1} + \delta_{x2}) \tag{6}

or, in matrix form,

{F1F2}=EAL[1βˆ’1βˆ’11]{Ξ΄x1Ξ΄x2}(7)\begin{Bmatrix} F_1\\ F_2 \end{Bmatrix} = \frac{EA}{L} \begin{bmatrix} 1 & -1\\ -1 & 1 \end{bmatrix} \begin{Bmatrix} \delta_{x_1}\\ \delta_{x_2}\\ \end{Bmatrix} \tag{7}

This gives us the same force-displacement relationship, with the familiar element stiffness matrix emerging, via a slightly different storyline!

Next Lesson
13. The Principle of Minimum Potential Energy