The finite element analysis of high frequency vibrations of quartz crystal plates is a necessary process required in the design of quartz crystal resonators of precision types for applications in filters and sensors. The anisotropic materials and extremely high frequency in radiofrequency range of resonators determine that vibration frequency spectra are complicated with strong couplings of large number of different vibration modes representing deformations which do not appear in usual structural problems. For instance, the higher-order thickness-shear vibrations usually representing the sharp deformation of thin plates in the thickness direction, expecting the analysis is to be done with refined meshing schemes along the relatively small thickness and consequently the large plane area. To be able to represent the precise vibration mode shapes, a very large number of elements are needed in the finite element analysis with either the three-dimensional theory or the higher-order plate theory, although considerable reduction of numbers of degree-of-freedom (DOF) are expected for the two-dimensional analysis without scarifying the accuracy. In this paper, we reviewed the software architecture for the analysis and demonstrated the evaluation and tuning of parameters for the improvement of the analysis with problems of elements with a large number of DOF in each node, or a problem with unusually large bandwidth of the banded stiffness and mass matrices in comparison with conventional finite element formulation. Such a problem can be used as an example for the optimization and tuning of problems from multi-physics analysis which are increasingly important in applications with excessive large number of DOF and bandwidth in engineering.

A nonlinear analysis of high-frequency thickness-shear vibrations of AT-cut quartz crystal plates is presented with the two-dimensional finite element method. We expanded both kinematic and constitutive nonlinear Mindlin plate equations and then truncated them to the first-order equations as an approximation, which is used later for the formulation of nonlinear finite element analysis with all zeroth- and first-order displacements and electric potentials. The matrix equation of motion is established with the first-order harmonic approximation and the generalized nonlinear eigensystem is solved by a direct iterative procedure. A backbone curve and corresponding mode shapes are obtained and analyzed. The nonlinear finite element program is developed based on earlier linear edition and can be utilized to predict nonlinear characteristics of miniaturized quartz crystal resonators in the design process.