NUMERICAL ALGORITHMS FOR SOLUTIONS OF LARGE EIGENVALUE PROBLEMS IN PIEZOELECTRIC RESONATORS

Two algorithms for eigenvalue problems in piezoelectric finite element analyses are introduced. The first algorithm involves the use of Lanczos method with a new matrix storage scheme, while the second algorithm uses a Rayleigh quotient iteration scheme. In both solution methods, schemes are implemented to reduce storage requirements and solution time. Both solution methods also seek to preserve the sparsity structure of the stiffness matrix to realize major savings in memory. In the Lanczos method with the new storage scheme, the bandwidth of the stiffness matrix is optimized by mixing the electrical degree of freedom with the mechanical degrees of freedom. The unique structural pattern of the consistent mass matrix is exploited to reduce storage requirements. These major reductions in memory requirements for both the stiffness and mass matrices also provided large savings in computational time. In the Rayleigh quotient iteration method, an algorithm for generating good initial eigenpairs is employed to improve its overall convergence rate, and its convergence stability in the regions of closely spaced eigenvalues and repeated eigenvalues. The initial eigenvectors are obtained by interpolation from a coarse mesh. In order for this multi-mesh iterative method to be effective, an eigenvector of interest in the fine mesh must resemble an eigenvector in the coarse mesh. Hence, the method is effective for finding the set of eigenpairs in the low-frequency range, while the Lanczos method with a mixed electromechanical matrix can be used for any frequency range. Results of example problems are presented to show the savings in solution time and storage requirements of the proposed algorithms when compared with the existing algorithms in the literature