Piezoelectricity in a Layered Shell

This tutorial is intended as a simple example showing how to model piezoelectric devices using the layered shell functionality.

Two cases of material orientation are investigated. In the first case, the pole axis is normal to the shell surface, which results in a change in thickness of the deformed shell. In the second case, the pole axis is tangential to the shell, which leads to the shell bending.

This model requires either MEMS or AC/DC Module.

Model Definition

In this tutorial, you model a thin structure with three layers. The top and bottom layers are electric conductors, and the middle one is a piezoelectric material. The geometry of the structure is shown in Figure 1.

Figure 1: Layered shell geometry representation in 3D.

The modeling approach in this tutorial is based on the layered shell technology available in COMSOL Multiphysics. Thus, you represent the geometry by a surface in 3D space, which gives the base of the layered structure; see Figure 2.

Figure 2: Model geometry representation as a surface.

You also use the surface for meshing, which gives the in-plane part of the discretized representation of the model; see Figure 3.

Figure 3: Meshed base selection.

The out-of-plane structure of the shell is modeled using an extra dimension attached to the base selection. The extra dimension together with the base selection span a local 3-dimensional space in which the equations are solved.

The extra dimension is set up using a special Layered Material node, which allows you to change the number of layers as well as their thickness and in-plane rotation, and to change the number of mesh elements used in each layer in the extra dimension. You can also select the materials to be used in each layer; see Figure 4.

Figure 4: Layered material.

MaterialS

The structure in this tutorial has three layers. The top and bottom ones represent 1-mm thick conductors made of aluminum. The layer in the middle has a thickness of 4 mm and is composed of a piezoelectric material called Lead Zirconate Titanate (PZT-5H).

For a piezoelectric material, the inelastic contributions to strain due to applied electric field E can be written in the strain-charge form as

For materials with tetragonal symmetry (class 4 mm) such as PZT-5H, the coupling matrix has the following structure (Ref. 1):

(1)

Here, the pole direction is along the third coordinate axis (the z-axis). Such orientation is assumed for all the piezoelectric material data available in COMSOL Material Library. For PZT-5H, the following values are used: d31 = −2.74·10−10 C/N, d33 = 5.93·10−10 C/N, and d51 = 7.41·10−10 C/N.

When the electric field has only one component, Ez, the only nonzero contributions to the strain are: εxx = εyy = d31Ez and εzz = d33Ez.

In this example, the electric field component Ez is positive. Hence, the piezoelectric layer will tend to stretch in the z direction and shrink in the other two directions.

When the pole direction is changed to the x-axis, the first and last lines in the coupling matrix will be swapped. As a result, the only nonzero contribution to the strain will be εxz = d15Ez so that the piezo layer will be sheared in the xz-plane.

Constraints

At x = 0, the whole cross section of the layered shell (the yz-plane) is mechanically fixed; see Figure 1. At the same cross section, the ground condition is applied to the top conducting layer, and a fixed electric potential of 20 V is applied to the bottom one. All other boundaries are free.

Results and Discussion

Two separate stationary studies are used for computing the electric potential and deformation for two different cases of piezoelectric material orientation. In both cases, the electric potential distribution is similar; see Figure 5.

Figure 5: Electric potential distribution.

In the first study, the material has a default orientation with the z-axis as the pole direction, which is normal to the shell. The structure deformation is shown in Figure 6.

In the second study, the material orientation is changed so that the pole direction is in-plane along the x-axis. The nature of resulting structural deformation is significantly different; see Figure 7.

The results show that using a correct choice of the material orientation is essential when modeling piezoelectric applications.

Figure 6: Vertical displacement for the case where the z-axis is the pole direction in the piezoelectric layer.

Figure 7: Vertical displacement for the case where the x-axis is the pole direction in the piezoelectric layer.

Notes About the COMSOL Implementation

The device is modeled in 3D using the predefined multiphysics interface . Two physics interfaces, structural and , will automatically be added to the model together with a multiphysics coupling feature called .

The interface will contain a node, where you can select geometry boundary for the base surface, and also select certain layers within the layered material. The interface will contain a node , where a similar selection should be made. It is important to have the same selections under both interfaces.

All settings for the material properties and orientation can be found under the node within the interface. This includes the structural, dielectric, and coupling properties.

Reference

1. A.L. Kholkin, N.A. Pertsev, and A.V. Goltsev, “Piezoelectricity and Crystal Symmetry.” In A. Safari and E.K. Akdogan (eds.), Piezoelectric and Acoustic Materials for Transducer Applications, Springer, Boston, MA, 2008.

Composite_Materials_Module/Multiphysics/piezoelectric_layered

Modeling Instructions

From the menu, choose .

New

In the window, click

Model Wizard

In the window, click

In the tree, select .

Click .

Click

In the tree, select .

Click

Geometry 1

In the window, under click .

In the window for , locate the section.

From the list, choose .

Work Plane 1 (wp1)

In the toolbar, click

Work Plane 1 (wp1)>Plane Geometry

In the window, click .

Work Plane 1 (wp1)>Rectangle 1 (r1)

In the toolbar, click

In the window for , locate the section.

In the text field, type 50.

In the text field, type 30.

In the window, right-click and choose .

Click the

button in the toolbar.

This geometry represents the base surface of the layered shell.

Next, add the materials aluminum and PZT-5H from the Material Library.

Add Material

In the toolbar, click

to open the window.

Go to the window.

In the tree, select .

Click

In the tree, select .

Click

In the toolbar, click

to close the window.

Add a layered material node and defined a cross section of the layered shell. Here you can also change the number of mesh elements used in the layered shell extra dimension.

Global Definitions

Layered Material 1 (lmat1)

In the window, under right-click and choose .

In the window for , locate the section.

Click

twice.

In the table, enter the following settings:


Layer	Material	Rotation (deg)	Thickness	Mesh elements
Layer 1	Aluminum (mat1)	0.0	1[mm]	2
Layer 2	Lead Zirconate Titanate (PZT-5H) (mat2)	0.0	4[mm]	4
Layer 3	Aluminum (mat1)	0.0	1[mm]	2

You can preview the shell cross section.

Click to expand the section. In the text field, type 8/30.

Locate the section. Click in the upper-right corner of the section.

Add a layered material link node to make the layered material available in . Using this node, you can also select a coordinate system defining the orientation. Note that only boundary coordinate systems can be selected.

Materials

Layered Material Link 1 (llmat1)

In the window, under right-click and choose .

Next, select only the middle layer in the shell to use the piezoelectric material.

Electric Currents in Layered Shells (ecis)

Piezoelectric Layer 1