This work presents a crack propagation analysis procedure on a computer program based on the consistent boundary element formulation for two-dimensional problems. This method has as one of its main features the exact solution of the singularity problems present in the formulation. In addition, with this methodology it is possible to represent the crack geometry with micrometric openings, similar to the cracks presented in laboratory tests. In this study, the propagation results in three structures with different geometries are analyzed, each structure subjected to different load combinations, in order to reproduce pure loading modes I and II, as well as mixed loading modes. A study is carried out on the size of the increments used in the models and on the propagation angle, making it possible to determine that the ideal size of the elements of new sections should be limited to the same dimension of the neighboring elements, avoiding possible numerical errors, while the propagation angle can be determined using the stress intensity factors (FIT) KI e KII, employing the concept of maximum principal stress. The FIT is obtained through reciprocal displacements close to the crack tip, and a study is carried out with a reference example to measure the reliability of the technique, with differences of at most 7 per cent. The performance observed using the methodology adopted in this study is compared with other results found in literature, showing similar crack propagation paths in all simulations. In the course of the chapters, the concepts of linearly elastic fracture mechanics and the adopted crack geometry are explained, as well as the development of the computational code.