Total number of DOFs of a member is the summation of DOFs of the two nodes. Therefore, for a 3D beam, the stiffness matrix is of size 12 x 12. The elastic stiffness matrix in the local coordinate system with shear deformation is given [Ref. 8] as follows:

where: ; ; ; ;