Piezoelectric materials have widespread applications in modern technical areas such as mechatron- ics, smart structures or microsystem technology, where they serve as sensors or actuators. For the assessment of strength and reliability of piezoelectric structures under combined electrical and mechanical loading, the existence of cracklike defects plays an important role. Meanwhile, piezoelectric fracture mechanics has been established quite well, but its application to realistic crack configurations and loading situations in piezoelec- tric structures requires the use of numerical techniques as finite element methods (FEM) or boundary element methods (BEM). The aim of this paper is to review the state of the art of FEM to compute the coupled electro- mechanical boundary value problem of cracks in 2D and 3D piezoelectric structures under static and dynamic loading. In order to calculate the relevant fracture parameters very precisely and efficiently, the numerical treatment must account for the singularity of the mechanical and electrical fields at crack tips. The follow- ing specialized techniques are presented in detail: (1) special singular crack tip elements, (2) determination of intensity factors KI-KIV from near tip fields, (3) modified crack closure integral, (4) computation of the electromechanical J -integral, and (5) exploitation of interaction integrals. Special emphasis is devoted to a realistic modeling of the dielectric medium inside the crack, leading to specific electric crack face boundary conditions. The accuracy, efficiency, and applicability of these techniques are examined by various example problems and discussed with respect to their advantages and drawbacks for practical applications.