Obtaining the correct piezoelectric material properties from a manufacturer datasheet and transforming those properties into the correct format for simulation can be a hassle.
Let’s take the example of the datasheet below. We will show you how to convert the soft ceramic material “C6” present on this datasheet into a format which can be inserted into the material database of OnScale.
Before we explain how to do that, first we need to look at some of the equations defining piezoelectricity to understand more closely what those material coefficients mean.
What is the Piezoelectric Effect?
There’s no shame in refreshing our memory before going into the equations, as you will see, it will be extremely useful. The piezoelectric effect is a complex way to say that a pressure on a special material will generate some electric voltage. The piezoelectric effect is a reversible process so if you apply a voltage at the 2 extremities of a piezo material, it will then deform mechanically.
Equations of Piezoelectricity
The piezoelectric effect that we previously described can be translated into 2 important behavior laws which describe the combined effect of electrical and elastic mechanical behavior.
The linear electrical behavior of the material:
 D is the electric charge density displacement (electric displacement)
 ε is permittivity (freebody dielectric constant)
 E is electric field strength
The Hooke’s law for elastic materials:
 S is strain
 s is compliance under shortcircuit conditions
 T is stress
These relations may be combined into socalled coupled equations, of which the straincharge form is:
Of course, we need to think about those equations in terms of their matrix equivalent.
Straincharge Matrix relation for the tetragonal PZT C4V
For example, the straincharge for a material of the 4mm (C4v) crystal class (such as a poled piezoelectric ceramic such as tetragonal PZT or BaTiO3) as well as the 6mm crystal class may also be written as (ANSI IEEE 176):
These are the relations we will use to practically calculate the coefficients afterwards, so those matrices are very important.
The 4 Constitutive Piezoelectric equations
In practice, piezoelectric coefficients can be defined in four ways as follows:
And because of that, we can define the previous 2 equations in 4 different “formats” in function of the coefficients we want to use:
Following is a description of all matrix variables used in the piezoelectric constitutive equations:
Symbol  Object Type  Size  Units  Meaning 
vector

6 x 1

stress components (e.g. s_{1})


vector

6 x 1

strain components (e.g. e_{3})


vector

3 x 1

electric field components


vector

3 x 1

electric charge density displacement components


matrix

6 x 6

compliance coefficients


matrix

6 x 6

stiffness coefficients


matrix

3 x 3

electric permittivity


matrix

3 x 6

piezoelectric coupling coefficients for StrainCharge form

What to do with all that?
In a second, more practical approach, Part 2 will show you:
1 Which coefficients OnScale need
2 A practical algorithm in Python to calculate all those coefficients
3 How to input them in OnScale as a new custom piezoelectric material
Ready for Part 2 of this blog tutorial? We start by answering what parameters are needed in OnScale to define a Piezoelectric material. Click here to start reading!