6.1.1 Eigenvalue analysis of a piezoelectric transducer

Product: ABAQUS/Standard  

This problem performs an eigenspectrum analysis of a cylindrical transducer consisting of a piezoelectric material with brass end caps. Various elements are used in the analysis. The elements range from axisymmetric elements to three-dimensional elements, using both lower- and higher-order elements. The basis of the piezoelectric capability in ABAQUS is described in Piezoelectric analysis, Section 2.10.1 of the ABAQUS Theory Manual.

Geometry and material

This problem is identical to the one discussed in a report by Mercer et al. (1987). The structure is shown in Figure 6.1.1–1 and consists of a piezoelectric material PZT4 with brass end caps. The piezoelectric material is electroded on both the inner and outer surfaces.

The properties for PZT4 in a cylindrical system are:

Elasticity Matrix:

Piezoelectric Coupling Matrix (Stress Coefficients):

Dielectric Matrix:

The 1-direction is radial, the 2-direction is axial, and the 3-direction is tangential. From these matrices it is seen that the poling direction is radially outwards from the axis of symmetry. (The order of the stresses in ABAQUS may differ from those typically used in electrical applications. ABAQUS uses the standard mechanical convention, where the stress components are ordered as . See Piezoelectric behavior, Section 20.6.2 of the ABAQUS Analysis User's Manual.)

The brass is elastic and isotropic with a Young's modulus of 104 GPa and a Poisson's ratio of 0.37.

Models

The transducer is modeled with a variety of elements. It is modeled as an axisymmetric structure utilizing both the planar, axisymmetric elements and the three-dimensional elements. For the axisymmetric elements, five meshes employing 4-node, 6-node, and 8-node elements are used in the finite element discretization. The first two meshes use 4-node elements with two levels of refinement, the third mesh uses 6-node elements, and the last two meshes use the 8-node elements with two levels of refinement. Lumped mass matrices are used for the lower-order elements. Consistent mass matrices are used in the higher-order elements. The meshes used for the 4-node and 6-node axisymmetric elements are shown in Figure 6.1.1–2. The meshes used for the 8-node axisymmetric elements are shown in Figure 6.1.1–3.

The three-dimensional model uses a slice of the structure and applies axisymmetric boundary conditions; 8-node and 20-node brick elements are used. The discretization used for each model is shown in Figure 6.1.1–3. These models use the *ORIENTATION option to maintain the proper definitions of the material properties. Also, in order to prescribe the axisymmetric boundary conditions, the nodal degrees of freedom are transformed into a cylindrical coordinate system.

All the models are considered to be open-circuited. The potentials on the inside surface are restrained to zero. The frequencies correspond to those for anti-resonance.

Results and discussion

The solutions obtained with the various ABAQUS models are shown in Table 6.1.1–1. Even for these coarse models, the results are quite close to the experimental results. In addition, the results from ABAQUS for the lower-order axisymmetric elements with lumped mass and the higher-order axisymmetric elements with consistent mass matrices in the computation of both the resonant and antiresonant frequencies match well with the numerical results reported in Mercer et al. The first four mode shapes for the more refined model with CAX8RE elements are shown in Figure 6.1.1–4.

Similar analyses have been performed considering the problem to be closed-circuited to obtain the resonant frequencies. For this situation, the potentials on both the inner and outer surfaces are set to zero. The results also compare well with those given in Mercer et al.

Input files

eigenpiezotrans_cax4e_coarse.inp

Coarse mesh with 4-node axisymmetric elements.

eigenpiezotrans_cax4e_fine.inp

Refined mesh with 4-node axisymmetric elements.

eigenpiezotrans_cax6e.inp

Mesh with 6-node axisymmetric elements.

eigenpiezotrans_cax8re_coarse.inp

Coarse mesh with 8-node axisymmetric elements.

eigenpiezotrans_cax8re_fine.inp

Refined mesh with 8-node axisymmetric elements.

eigenpiezotrans_c3d8e.inp

Mesh with 8-node three-dimensional elements.

eigenpiezotrans_c3d8e.f

User subroutine ORIENT used in eigenpiezotrans_c3d8e.inp.

eigenpiezotrans_c3d20e.inp

Mesh with 20-node three-dimensional elements.

eigenpiezotrans_c3d20e.f

User subroutine ORIENT used in eigenpiezotrans_c3d20e.f.

eigenpiezotrans_elmatrixout.inp

Data that test the use of *ELEMENT MATRIX OUTPUT with piezoelectric elements.

eigenpiezotrans_usr_element.inp

Data for a job that reads in the matrices output in eigenpiezotrans_elmatrixout.inp and performs an eigenvalue analysis.

Reference

Table

Table 6.1.1–1 Piezoelectric transducer eigenvalue estimates.

ModelFrequencies (kHz) for mode number
Type# of Elements12345
CAX4E1314.139.156.266.179.3
CAX4E32018.640.357.864.288.1
CAX6E1020.043.263.270.498.8
CAX8RE519.642.861.066.996.3
CAX8RE8018.640.357.664.287.6
C3D8E1619.841.862.068.795.2
C3D20E1619.742.960.466.591.7
Experimental(1)18.635.454.263.388.8
(1): Experimental results obtained from Mercer et al. (1987).


Figures

Figure 6.1.1–1 Piezoelectric transducer.

Figure 6.1.1–2 Meshes used with 4-node and 6-node axisymmetric elements.

Figure 6.1.1–3 Meshes used with 8-node axisymmetric and 8-node and 20-node three-dimensional elements.

Figure 6.1.1–4 Mode shapes for 8-node axisymmetric elements.