Piezoelectric Behavior

A piezoelectric material responds to an electric potential gradient by straining, while stress causes an electric potential gradient in the material. This coupling between electric potential gradient and strain is the material's piezoelectric property. The material will also have a dielectric property so that an electrical charge exists when the material has a potential gradient. Piezoelectric material behavior is discussed in Piezoelectric analysis.

The mechanical properties of the material must be modeled by linear elasticity (Linear Elastic Behavior). The mechanical behavior can be defined by

in terms of the piezoelectric stress coefficient matrix, $e_{mij}^{\phi }$, or by

in terms of the piezoelectric strain coefficient matrix, $d_{mkl}^{\phi }$. The electrical behavior is defined by

where

${\sigma }_{ij}$

is the mechanical stress tensor;

${\epsilon }_{ij}$

is the strain tensor;

$q_i$

is the electric “displacement” vector;

$D_{ijkl}^E$

is the material's elastic stiffness matrix defined at zero electrical potential gradient (short circuit condition);

$e_{mij}^{\phi }$

is the material's piezoelectric stress coefficient matrix, defining the stress ${\sigma }_{ij}$ caused by the electrical potential gradient $E_m$ in a fully constrained material (it can also be interpreted as the electrical displacement $q_m$ caused by the applied strain ${\epsilon }_{ij}$ at a zero electrical potential gradient);

$d_{mkl}^{\phi }$

is the material's piezoelectric strain coefficient matrix, defining the strain ${\epsilon }_{kl}$ caused by the electrical potential gradient $E_m$ in an unconstrained material (an alternative interpretation is given later in this section);

$\phi$

is the electrical potential;

$D_{ij}^{\phi \left(\epsilon \right)}$

is the material's dielectric property, defining the relation between the electric displacement $q_i$ and the electric potential gradient $E_j$ for a fully constrained material; and

$E_i$

is the electrical potential gradient vector, $-\partial \phi /\partial x_i$.

The material's electrical and electro-mechanical coupling behaviors are, thus, defined by its dielectric property, $D_{ij}^{\phi \left(\epsilon \right)}$, and its piezoelectric stress coefficient matrix, $e_{mij}^{\phi }$, or its piezoelectric strain coefficient matrix, $d_{mkl}^{\phi }$. These properties are defined as part of the material definition (Material Data Definition).

Alternative forms of the piezoelectric constitutive equations are presented in this section. These forms of the equations involve material properties that cannot be used directly as input for Abaqus/Standard. However, they are related to the Abaqus/Standard input through simple relations that are presented in Piezoelectric analysis. The intent of this section is to draw connections between the Abaqus/Standard terminology and input to that used commonly in the piezoelectricity community. The mechanical behavior can also be defined by

in terms of the piezoelectric coefficient matrix $g_{mij}^{\phi }$, and the stiffness matrix $D_{ijkl}^q$, which defines the mechanical properties at zero electrical displacement (open circuit condition). Likewise, the electrical behavior can also be defined by

in terms of the dielectric matrix $D_{ij}^{\phi \left(\sigma \right)}$ for an unconstrained material or by

where

$D_{ijkl}^q$

is the material's elastic stiffness matrix defined at zero electrical displacement;

$d_{ijk}^{\phi }$

is the material's piezoelectric strain coefficient matrix used earlier, and based on the equations, may alternatively be interpreted as the electrical displacement $q_i$ caused by the stress ${\sigma }_{jk}$ at zero electrical potential gradient;

$g_{mkl}^{\phi }$

is the material's piezoelectric coefficient matrix, which can be interpreted as defining either the strain ${\epsilon }_{kl}$ caused by the electrical displacement $q_m$ in an unconstrained material or the electrical potential gradient $E_m$ caused by the stress ${\sigma }_{kl}$ at zero electrical displacement; and

$D_{ij}^{\phi \left(\sigma \right)}$

is the material's dielectric property, defining the relation between the electric displacement $q_i$ and the electric potential gradient $E_j$ for an unconstrained material.

These are useful relationships that are often seen in the piezoelectric literature. In Piezoelectric analysis the properties $g_{mkl}^{\phi }$, $D_{ij}^{\phi \left(\sigma \right)}$, and $D_{ijkl}^q$ are expressed in terms of the properties $d_{mkl}^{\phi }$, $D_{ij}^{\phi \left(\epsilon \right)}$, and $D_{ijkl}^E$, that are used as input for Abaqus/Standard.

The dielectric matrix can be isotropic, orthotropic, or fully anisotropic. For non-isotropic dielectric materials a local orientation for the material directions must be specified (Orientations). The entries of the dielectric matrix (what are referred to as “dielectric constants” in Abaqus) refer to what is more commonly known in the literature as the permittivity of the material.

The piezoelectric material properties can be defined by giving the stress coefficients, $e_{mij}^{\phi }$ (this is the default), or by giving the strain coefficients, $d_{mkl}^{\phi }$. In either case, 18 components must be given in the following order (substitute d for e for strain coefficients):

The first index on these coefficients refers to the component of electric displacement (sometimes called the electric flux), while the last pair of indices refers to the component of mechanical stress or strain.

Thus, the piezoelectric components causing electrical displacement in the 1-direction are all given first, then those causing electrical displacement in the 2-direction, and then those causing electrical displacement in the 3-direction. (Some references list these coupling terms in a different order.)

PIEZOELECTRIC, TYPE=S

PIEZOELECTRIC, TYPE=E

Property module: material editor: Electrical/MagneticPiezoelectric: Type: Stress or Strain

Abaqus does not account for piezoelectric effects in the total energy balance equation, which can lead to an apparent imbalance of the total energy of the model in some situations. For example, if a piezoelectric truss is fixed at one end point and subjected to a potential difference between its two end points, it deforms due to the piezoelectric effect. Subsequently if the truss is held fixed in this deformed configuration and the potential difference removed, strain energy will be generated due to the constraints. This results in an equivalent increase in the total energy of the model.