Steps to Derive k and Assemble into K

Figure 1 is a flowchart illustrating the sequence for computing the stiffness matrix of a simple problem (the concept is similar to more complex problems). Computing the stiffness matrix involves two main steps:

1. Derive local stiffness matrices (k)
2. Assemble k into K, the global stiffness matrix

Matrix Methods: Direct vs. Iterative

Direct vs. Iterative methods
The two approaches available for solving global stiffness matrix (K) in FEM are:

Implications on FEA solver

• For linear simulation e.g. KU = f , Gaussian elimination can be applied directly
• For nonlinear simulation e.g. K(u)U=f, stiffness is dependent on displacements (u). Therefore an iterative method must be used.

See also
The role of matrix solution in the FEM process