VM231

VM231
Piezoelectric Rectangular Strip Under Pure Bending Load

Overview

Reference: V. Z. Parton, B. A. Kudryavtsev, N. A, Senik, Applied Mechanics: Soviet Review, Vol. 2: Electromagnetoelasticity, pg. 28.
Analysis Type(s): Static
Element Type(s):
2D Coupled-Field Solid Elements (PLANE13 and PLANE222)
2D 8-Node Coupled-Field Solid Elements (PLANE223)
Input Listing: vm231.dat

Test Case

A piezoceramic (PZT-4) rectangular strip occupies the region |x| l, |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 = ± l. Determine the electro-elastic field distribution in the strip.

This test case is also solved using Ansys Mechanical. See vm-wb-mech-109.

Figure 387: Piezoelectric Strip Problem Sketch

Piezoelectric Strip Problem Sketch

Material PropertiesGeometric PropertiesLoading
l = 1mm
h = 0.5 mm
Pressure slope
σ1 = -20 N/mm3

Constitutive Matrices for PZT4 Polarized Along the Y Axis

PZT4 Dielectric Matrix [εr]

PZT4 Piezoelectric Matrix [e] C/m2

PZT4 Stiffness matrix [c] x 10-10 N/m2

Analysis Assumptions and Modeling Notes

The problem is solved using three different coupled field elements, namely, PLANE13, PLANE222, and PLANE223 elements.

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

For the first case, the finite element model uses a single 4-node, 2D quadrilateral element (PLANE13) with extra displacement and potential modes to produce a correct response to the pure bending load. For the second case, the finite element model uses a single 4-node, 2D quadrilateral element (PLANE222) with enhanced strain formulation. For the third case, the finite element model uses a single, 8-node 2D quadrilateral element (PLANE223).

Results Comparison

Node 3TargetMechanical APDLRatio
PLANE13
UX, μm-0.110-0.1101.000
UY, μm0.1150.1151.000
VOLT, V27.37827.3791.000
Stress x, N/mm2 -10.000-10.0001.000
Electric Field y, V/mm-109.511-109.5151.000
PLANE222
UX, μm-0.110-0.1101.000
UY, μm0.1150.1151.000
VOLT, V27.37827.3791.000
Stress x, N/mm2-10.000-10.0001.000
Electric Field y, V/mm-109.511-109.5151.000
PLANE223
UX, μm-0.110-0.1101.000
UY, μm0.1150.1151.000
VOLT, V27.37827.3791.000
Stress x, N/mm2 -10.000-10.0001.000
Electric Field y, V/mm-109.511-109.5151.000