A double projection Bubnov Galerkin methods is developed and is shown to yield identical results as the Petrov Galerkin method based on operator symmetrization. Using the derivations of this paper it is shown that cubic elements applied to the Timoshenko beam problem yield exact nodal values.