A structural geometric nonlinear analysis, using the finite element method (FEM), depends on the consideration of five aspects: the bending theory, the kinematic description, the strain-displacement relations, the nonlinear solution scheme and the interpolation (shape) functions. As MEF is a numerical solution, the structure discretization provides great influence on the analysis response. However, applying shape functions calculated from the homogenous solution of the differential equation of the problem, the exact behavior of the structure is obtained for a minimum discretization, as for a linear analysis. Thus, this work aims to integrate the solutions for the formulations of the geometric nonlinearity problem, in order to reduce this influence and allow a minimum discretization of the structure, also considering, large displacements and rotations. Then, using an updated Lagrangian kinematic description, considering a higher-order Green strain tensor, The Euler-Bernoulli and Timoshenko beam theories, the nonlinear solutions schemes and the interpolation functions, that includes the influence of axial force, obtained directly from the solution of the equilibrium differential equation of an deformed infinitesimal element, a spatial bar frame element is developed using a complete formulation. The element was implemented in the Framoop, and their results, for a minimum discretization, were compared with conventional formulations, analytical solutions and with the software Mastan2 v3.5. Results clearly show the efficiency of the developed formulation to predict the critical load of plane and spatial structures using a minimum discretization.