AbstractThis work presents a hybrid shear‐flexible beam‐element, capable of capturing arbitrarily large inelastic displacements and rotations of planar frame structures with just one element per member. Following Reissner's geometrically exact theory, the finite element problem is herein formulated within nonlinear programming principles, where the total potential energy is treated as the objective function and the exact strain‐displacement relations are imposed as kinematic constraints. The approximation of integral expressions is conducted by an appropriate quadrature, and by introducing Lagrange multipliers, the Lagrangian of the minimization program is formed and solutions are sought based on the satisfaction of necessary optimality conditions. In addition to displacement degrees of freedom at the two element edge nodes, strain measures of the centroid act as unknown variables at the quadrature points, while only the curvature field is interpolated, to enforce compatibility throughout the element. Inelastic calculations are carried out by numerical integration of the material stress‐strain law at the cross‐section level. The locking‐free behavior of the element is presented and discussed, and its overall performance is demonstrated on a set of well‐known numerical examples. Results are compared with analytical solutions, where available, and outcomes based on flexibility‐based beam elements and quadrilateral elements, verifying the efficiency of the formulation.