VM-WB-MECH-109

VM-WB-MECH-109
Piezoelectric Rectangular Strip Under Pure Bending Load

Overview

Reference: Parton, V. Z., Kudryavtsev, B. A., & Senik, N. A. (1990). Applied Mechanics: Soviet Review, Vol. 2: Electromagnetoelasticity. CRC Press.
Solver

Ansys Mechanical

Analysis Type(s): Coupled Field Static
Element Type(s):

2D Plane Stress

Test Case

A piezoceramic (PZT-4) rectangular strip occupies the region |x| ≤ ℓ, |y| ≤ h. The material is oriented such that its polarization direction is aligned with the Y axis. The strip is subjected to the pure bending load σx = σ1y at x = ± ℓ. Determine the electro-elastic field distribution in the strip.

This test case is also solved using Ansys Mechanical APDL. See VM231.

Figure 157: Problem Sketch: Piezoelectric Strip

Problem Sketch: Piezoelectric Strip

Figure 158: Finite Element Model

Finite Element Model

Material PropertiesLoading
ℓ = 1 mm and h = 0.5 mmPressure Slope = -20 N/mm3

Materials Properties Tables

The piezoelectric material properties of PZT4 defined in Engineering Data are given below:

Figure 159: Anisotropic Elasticity :

Anisotropic Elasticity :

Figure 160: Anisotropic Relative Permittivity (Constant Strain):

Anisotropic Relative Permittivity (Constant Strain):

Figure 161: Piezoelectric Matrix (Stress)

Piezoelectric Matrix (Stress)

Analysis Assumptions and Modeling Notes

The finite element model uses a single 8-node, 2-D quadrilateral element (PLANE223).

Only a one-quarter symmetry sector of the rectangular strip is modeled. Symmetric and antisymmetric boundary conditions are applied as shown below.

Figure 162: Symmetric Boundary Condition

Symmetric Boundary Condition

Figure 163: Anti-Symmetric Boundary Condition

Anti-Symmetric Boundary Condition

Figure 164: Boundary Conditions and Loading

Boundary Conditions and Loading

Results Comparison

The Workbench results are compared with target specified in reference book.

Result at Node 3TargetWorkbench AppError (%)
Directional Deformation in X (µm)-0.10953E-07-0.10953E-070.0
Directional Deformation in X (µm)0.1493E-070.1493E-070.00
Electric Voltage (V)27.37927.3790.000
Equivalent Stress (N/mm2)10.000 x 10610 x 1060
Directional Electric Field Intensity in Y (V/m)-109.515E-03-109.515E-030.00

Figure 165: Directional Deformation in X

Directional Deformation in X

Figure 166: Directional Deformation in Y

Directional Deformation in Y