In this article, a theoretical model in nonlinear elasticity is presented to analyze the influence of the percentage of steel fibers and the effect of the fiber-concrete bond stress on the shear force at the rupture of beams subjected to the combined effect of bending moment, normal force, and shear force. For a given beam section, it is defined by a succession of layers of concrete and longitudinal steel elements. Each layer is defined by its height hi, width bi, and position relative to one end of the section YGi. Each longitudinal steel element is also defined by its cross-sectional area and position relative to one end of the section. The steel fibers are defined by volume percentages of 0.5%, 1%, 1.5%, and 2%, taking into account the mechanical nonlinearity of the materials. This model is based on the multilayer analysis of sections and an iterative solution procedure for each layer and each section, considering a given longitudinal deformation state and shear stress. The global equilibrium of the sections is analyzed under the assumption of flat longitudinal deformations but with, in principle, an interdependence of longitudinal normal stresses and shear stresses. In this study, using the principle of virtual work, equilibrium equations for deformations and stresses, as well as partial compatibility equations between concrete deformations and mean deformations, are derived. Comparative examples between ordinary reinforced concrete beams with variable shapes and reinforcement details, and those reinforced with steel fibers, are presented to demonstrate the accuracy of the proposed model for simulating the nonlinear shear response of beams.