Behavior of a viscoelastic and piezoelectric material under electromechanical loading has been systematically investigated in this study. The material has a viscoelastic and piezoelectric anisotropy. In terms of behavior the material responds to the environment prompting it through elastic stress, dissipative stress and electrical polarization fields, whose constitutive equations have been obtained. While elastic stress and electrical polarization fields are derived from the thermodynamic potential, dissipative stress was formed as a tensorial function that depends on certain arguments. After necessary information was obtained on the constitutive functions, power series expansions were made based on the assumption that such functions are analytic. Considering physical application conditions where mechanical interactions are assumed first to be linear and electro-mechanical ones - non-linear, orders of the terms in power series have been determined accordingly. Then assuming that deformation gradients, (or displacement gradients and rates of strain) are very small and electrical interactions are linear too; constitutive equations can be easily linearized. As a result, linear constitutive equations obtained were used as substitutes in balance equations, yielding field equations.