The paper presents a novel family of strain-based beam finite elements (FE) for analysis of tensile failure of a reinforced concrete bar (RC bar), with results of the analysis being independent of the applied FE mesh. The composite bar consists of a continuous longitudinal ductile reinforcing bar(s) surrounded by brittle concrete cover, which are considered separately in the model. Longitudinal slip at the contact between the concrete cover and reinforcing bars is allowed, while their relative displacements perpendicular to the axis of the RC bar are prevented. Cracks in concrete cover occur when tensile stress in concrete exceeds its tensile strength. Propagation of partially connected crack, that is, softening of the material at the crack, is described through constitutive law in form of nonlinear relationship between stresses in concrete at the crack and the width of the crack. Each separate crack is considered discretely as a discontinuity in geometry of the element. In the analysis of cracking of concrete, it is commonly assumed that the discrete crack can occur at the nodes of FE only. However, this assumption leads to dependence of the analysis results on the employed FE mesh. The presented family of FE enables occurrence of the crack anywhere along the FE. In order to achieve this, the discrete crack is excluded from equations of FE and additional boundary conditions are introduced at the discontinuity. This approach ensures that the location of the cracks, their number and their propagation are independent of the applied FE mesh. Advantages of the novel family of FE are thoroughly presented in a parametric study which investigates influence of number of FE as well as influence of degrees of interpolation and integration on the cracking of RC bar under tensile loading. Experimental results of tensile tests on the analysed bar are available in literature. It can be concluded that the results obtained with the minimal possible number of novel FE and sufficiently high degree of numerical integration scheme, applied for solving integrals in equations of FE, are considerably more accurate than the results of previous analyses with model of discrete crack at the nodes of FE only.