Equations of motion and closed-form solutions for spatial stability and free vibration analysis of shear flexible thin-walled elastic beams were analytically derived by applying Hellinger-Reissner principle in companion paper. In this paper, elastic and geometric stiffness matrices and mass matrix for finite element analysis are evaluated by using isoparametric and Hermitian interpolation polynomials, where isoparametric interpolation functions with 3 or 4 node per element are used in isoparametric element and third and fifth order Hermitian polynomials including shear deformation effects are used in Hermitian element. In order to verify the validity of this formulation, both analytic and numerical solutions for spatial buckling and free vibration problems including shear effects are presented and compared.