Piezoelectric materials are increasingly employed throughout the field of sensor and actuator applications. The mechanical and electrical fields are coupled by the constitutive relations. This results in incompatibilities when using standard finite element formulations. As a consequence parasitic approximations and incorrect computation results occur. The present work proposes a concept to avoid these incompatible approximation spaces and the thus occurring computation errors.