In the traditional boundary element methods, the numerical modelling of cracks is usually carried out by means of a hypersingular fundamental solution. A more natural procedure should make use of fundamental solutions capable of representing the 1/square root of r singularity that arises when one analyses the stress field near the crack tip. This representation has already been made by Dumont and Lopes in 2003, with some refinements achieved by Dumont and Mamani in 2011, in a formulation of the Hybrid Boundary Element Method, where the fundamental solutions are developed from Westergaard-type generalized stress functions for displacement-prescribed crack problems, as proposed by Tada et al in 1993. The present work, which is a continuation of research work done by Dumont and Mamani, carries out a detailed study on the use of these generalized functions for the representation of quantities at the crack tip for problems of elasticity and potential. The results obtained are compared conceptually with the classic developments by Westergaard and Williams. Some results were also analyzed with generalized stress functions for a crack with semi-elliptical opening shapes, besides the use of functions to represent the relative rotation of the crack faces. In addition, the application of the generalized Westergaard stress function as a fundamental solution in the Conventional Boundary Element Method is presented, more specifically for obtaining the matrix G of the system, since the matrix H has already been developed in previous works, with good results. Some numerical examples of application are presented for external boundary, holes and cracks.