Modeling of piezoelectric multilayer ceramics using finite element analysis

• Published in: IEEE transactions on ultrasonics, ferroelectrics, and frequency control ferroelectrics, and frequency control, 1997-11, Vol.44 (6), p.1204-1214
• Main Authors: Goldberg, R.L., Jurgens, M.J., Mills, D.M., Henriquez, C.S., Vaughan, D., Smith, S.W.
• Format: Article
• Language: English
• Subjects: Acoustics; Bioceramics; Biomedical transducers; Ceramic materials; Ceramics; Computer simulation; Electric lines; Exact sciences and technology; Finite element method; Finite element methods; Fundamental areas of phenomenology (including applications); Mathematical models; Medical imaging; Multilayers; Nonhomogeneous media; Physics; Piezoelectric transducers; Signal to noise ratio; Transduction; acoustical devices for the generation and reproduction of sound; Transmission line measurements; Ultrasonic imaging; Ultrasonic transducer arrays; Ultrasonic transducers
• ISSN: 0885-3010; 1525-8955
• DOI: 10.1109/58.656622

For medical ultrasound imaging, 2-D array transducers have greater versatility than linear arrays. Unfortunately, the tiny array elements in a 2-D array have poor signal-to-noise ratio (SNR). We have previously shown that SNR is increased in 2-D array transducers made from piezoelectric multilayer ceramics. Conventional one-dimensional models provide accurate results when comparing multilayer ceramic performance relative to single layer transducers. However, these models are not accurate when comparing simulations directly to measurements. Because multilayer ceramics have a complex structure, a 3-D model, such as finite element analysis, is needed for accurate simulations. We modeled four arrays that were previously fabricated: a single layer and multilayer 1 MHz, 2-D array element, and a single layer and multilayer 2.25 MHz, 1.5-D array element that can focus and steer in azimuth but only steer in the elevation dimension. We compared the simulated and measured impedance plots for each transducer. The finite element analysis plots accurately predicted the impedance for each vibration mode. On the other hand, the one dimensional KLM transmission line model could simulate only the thickness mode vibrations and the results were inaccurate compared to measurements. We also simulated the transmit output pressure for the 2.25 MHz arrays and compared the results to measurements. The simulated pressure vs. time plots and their spectra were accurate when compared to measurements. Finally, we obtained a series of images that show the impulse response vibrations for the 2.25 MHz, arrays. These animations show the vibration modes in the complex multilayer ceramic structure. Measurements were not available to confirm the animations. Our results show that finite element analysis in three dimensions is a valuable tool to predict the performance of multi-layer transducers.