In the paper, the method of a numerical simulation concerning diagonal crack propagation in con-crete beams was presented. Two beams reinforced longitudinally but without shear reinforcement were considered during the Finite Element Method analysis. In particular, a nonlinear method was used to simulate the crack evaluation in the beams. The analysis was performed using the commercial program ANSYS. In the numerical simulation, the limit surface for concrete described by Willam and Warnke was applied to model the failure of concrete. To solve the FEM-system of equations, the Newton-Raphson method was used. As the results of FEM calculations, the trajectories of total stains and numerical images of smeared cracks were obtained for two analyzed beams: the slender beam S5 of leff = 1.8 m and the short beam S3k of leff = 1.1 m. The applied method allowed to generate both flexural vertical cracks and diagonal cracks in the shear regions. Some differences in the evaluation of crack patterns in the beams were observed. The greater number of flexural vertical cracks which penetrated deeper in the beam S5 caused the lower stiffness and the greater deformation in the beam S5 compared to the short beam S3k. Numerical results were compared with the experimental data from the early tests performed by Słowik [3]. The numerical simulation yielded very similar results as the experiments and it confirmed that the character of failure process altered according to the effective length of the member. The proposed numerical procedure was successfully verified and it can be suitable for numerical analyses of diagonal crack propagation in concrete beams.