In computational fracture mechanics, great benefits are obtained from the reduced modeling dimension order and the accurate integral formulation of the boundary element method (BEM). However, the direct representation of co-planar surfaces (i.e., cracks) causes a degeneration of the standard displacement BEM formulation which can only be circumvented with special modeling techniques. Aiming to simplify the generalized application of the BEM to fracture mechanics problems, this paper presents a two-dimensional crack modeling approach. The method uses the direct BEM displacement formulation within a single-domain model to efficiently and precisely calculate any mixed mode crack tip stress intensity factor. Details of the application of the method are presented, while its accuracy and reliability are demonstrated through numerous comparisons with benchmark results.