This paper proposes a low-order geometrically exact flexible beam formulation based on the utilization of generic beam shape functions to approximate distributed kinematic properties of the deformed structure. The proposed nonlinear beam shapes approach is in contrast to the majority of geometrically nonlinear treatments in the literature in which element-based—and hence high-order—discretizations are adopted. The kinematic quantities approximated specifically pertain to shear and extensional gradients as well as local orientation parameters based on an arbitrary set of globally referenced attitude parameters. In developing the dynamic equations of motion, an Euler angle parametrization is selected as it is found to yield fast computational performance. The resulting dynamic formulation is closed using an example shape function set satisfying the single generic kinematic constraint. The formulation is demonstrated via its application to the modelling of a series of static and dynamic test cases of both simple and non-prismatic structures; the simulated results are verified using MSC Nastran and an element-based intrinsic beam formulation. Through these examples, it is shown that the nonlinear beam shapes approach is able to accurately capture the beam behaviour with a very minimal number of system states.