BOUNDARY
CONDITIONS
Polarization/Local
Axes
Data -> Conditions -> Volumes -> Geometry -> Polarization
(Cartesian) -> Define Polarization (P1) (0,0,1) direction
Choose 3PointsXZ -> Center (0,0,0) -> Define X axis by putting
(1,0,0) -> Then define Z axis by putting (0,0,1) -> Escape
Assign Polarization (Local-Axes P1) -> Choose all the PZT volumes
-> Finish
Potential
Data -> Conditions -> Surfaces -> Potential 0.0 (or Ground
in ATILA 5.2.4 or higher version) -> Assign the bottom (Ring
and Dot) 9 surfaces of the PZT
Forced Potential 1.0 -> Assign the top Ring 4 surfaces
Finally, Floating -> Assign the top Dot 5 surfaces
MESHING 
Mapped Mesh
Meshing -> Structured -> Volumes -> Select all
volumes -> Then Cancel the further process
Meshing -> Structured -> Lines -> Enter number to divide
line -> Choose "1" for
the lines -> Generate
Use the higher number such as 2
for the higher accurate calculation.