A procedure for the determination of stress intensity factors of 3D flat cracks is proposed, which uses integral identities for displacements and strains on the crack surface. New formulations of boundary problems of the mathematical theory of linear fracture mechanics have been obtained, which are used to develop finite element models. The unknown functions are nodal displacements on the crack surface. The data given are compared with analytical solutions and results obtained by other methods.Keywords: stress intensity factors, 3D cracks, integral representations of displacements.