Structural beam element is widely used for many applications. To satisfy some requirements, it is not uncommon that the beam has to posses a non-uniform distribution of its cross section along the span. A common approach to analyze such a structure is by using a variational principle to develop the element stiffness matrix. In the present paper, an
