Inversely Prepolarized Piezoceramic Cantilever. D.P. Sedorook and I.V. Ostrovskii

Similar documents
Solid State Physics (condensed matter): FERROELECTRICS

Transduction Based on Changes in the Energy Stored in an Electrical Field

Finite Element Analysis of Piezoelectric Cantilever

VIBRATION CONTROL OF RECTANGULAR CROSS-PLY FRP PLATES USING PZT MATERIALS

Laser Interferometric Displacement Measurements of Multi-Layer Actuators and PZT Ceramics

Piezoelectric Materials and Devices

COURSE OUTLINE. Introduction Signals and Noise Filtering Sensors: Piezoelectric Force Sensors. Sensors, Signals and Noise 1

Testing and analysis of high frequency electroelastic characteristics of piezoelectric transformers

Piezoelectric Vibration Energy Harvesting. Characteristics of Barium Titanate Laminates

ACOUSTIC EMISSION MEASUREMENTS ON PIEZOELECTRIC/ FERROELECTRIC MATERIALS

INTRODUCTION TO PIEZO TRANSDUCERS

Addition 1. Shear Stack Piezoelectric Elements and Shear Effect Basics

CHAPTER 4 DESIGN AND ANALYSIS OF CANTILEVER BEAM ELECTROSTATIC ACTUATORS

Piezoelectric Resonators ME 2082

Thickness Optimization of a Piezoelectric Converter for Energy Harvesting

FINITE ELEMENT MODELLING OF COMPOSITES USING PIEZOELECTRIC MATERIAL

433. New linear piezoelectric actuator based on traveling wave

PIEZOELECTRIC TECHNOLOGY PRIMER

PERFORMANCE OF HYDROTHERMAL PZT FILM ON HIGH INTENSITY OPERATION

Structural Health Monitoring using Shaped Sensors

Piezo materials. Actuators Sensors Generators Transducers. Piezoelectric materials may be used to produce e.g.: Piezo materials Ver1404

1302. Piezoelectric actuator for micro robot used in nanosatellite

Micromechanical modeling and simulation of piezoceramic materials

Piezoelectric materials for MEMS applications Hiroshi Funakubo Tokyo Institute of Technology

Piezoelectric Actuator for Micro Robot Used in Nanosatellite

ANALYSIS AND NUMERICAL MODELLING OF CERAMIC PIEZOELECTRIC BEAM BEHAVIOR UNDER THE EFFECT OF EXTERNAL SOLICITATIONS

SIMULATION AND OPTIMIZATION OF MEMS PIEZOELECTRIC ENERGY HARVESTER WITH A NON-TRADITIONAL GEOMETRY

Development of PZT and PZN-PT Based Unimorph Actuators for Micromechanical Flapping Mechanisms

Online publication date: 23 June 2010 PLEASE SCROLL DOWN FOR ARTICLE

SURFACE ACOUSTIC WAVE FERROELECTRIC PHONONIC CRYSTAL TUNABLE BY ELECTRIC FIELD

COUPLED FIELD ANALYSIS OF PIEZOELECTRIC CANTILEVER BEAM

I. INTRODUCTION II. SAMPLE PREPARATION JOURNAL OF APPLIED PHYSICS VOLUME 92, NUMBER 5 1 SEPTEMBER

The electric field induced strain behavior of single. PZT piezoelectric ceramic fiber

Optimizing the Design of Polymer Based Unimorph Actuator using COMSOL Multiphysics Vineet Tiwari, Rashiya Sharma, R. K. Dwivedi and Geetika Srivastava

Optimization of a circular piezoelectric bimorph for a micropump driver

The Analysis of Aluminium Cantilever Beam with Piezoelectric Material by changing Position of piezo patch over Length of Beam

Available online at ScienceDirect. Physics Procedia 70 (2015 )

Depolarization of a piezoelectric film under an alternating current field

DELAMINATION CONTROL IN COMPOSITE BEAMS USING PIEZOELECTRIC ACTUATORS

1990. Temperature dependence of soft-doped / hard-doped PZT material properties under large signal excitation and impact on the design choice

Effect of Strain Nodes and Electrode Configuration on Piezoelectric Energy Harvesting From Cantilevered Beams

Determination of full set material constants of piezoceramics from phase velocities

19 th INTERNATIONAL CONGRESS ON ACOUSTICS MADRID, 2-7 SEPTEMBER 2007 A METHOD FOR THE AUTOMATIC CHARACTERIZATION OF PIEZOELECTRIC CERAMIC ELEMENTS

A coupled field finite element model to predict actuation properties of piezoelectrically actuated bistable composites.

Simulation of a New PZT Energy Harvester with a Lower Resonance Frequency Using COMSOL Multiphysics. Zewail City Staff Sep, 2014

A flexoelectric microelectromechanical system on silicon

Single-phase driven ultrasonic motor using two orthogonal bending modes of sandwiching. piezo-ceramic plates

440. Simulation and implementation of a piezoelectric sensor for harmonic in-situ strain monitoring

Magneto-Mechanical Modeling and Simulation of MEMS Sensors Based on Electroactive Polymers

Piezo Theory: Chapter 1 - Physics & Design

Analysis of a ceramic bimorph piezoelectric gyroscope

Study of the Tip Deflection in Static State of a Piezoelectric Polymer based Bimorph Actuator with Varying Thickness and Length Ratios.

Epitaxial piezoelectric heterostructures for ultrasound micro-transducers

ELECTROMECHANICAL RESPONSE OF PIEZOELECTRIC FOAMS

Intelligent Materials and their applications. Autors: Jana Pintea, Ionut Balan INCDIE ICPE-CA Bucharest, ROMANIA

DAMPING CONTROL OF A PZT MULTILAYER VIBRATION USING NEGATIVE IMPEDANCE CIRCUIT

A method for matching the eigenfrequencies of longitudinal and torsional vibrations in a hybrid piezoelectric motor

Supplementary Information

Dielectric Properties of Solids

7.Piezoelectric, Accelerometer and Laser Sensors

DEVELOPMENT OF A BAR-SHAPED ULTRASONIC MOTOR FOR MULTI-DEGREES OF FREEDOM MOTION

Principles of Active Vibration Control: Piezoelectric materials

A High Power Piezoelectric Ultrasonic Linear Micromotor Using Slotted Stator

FLUID-STRUCTURE AND ELECTRIC INTERACTION ANALYSIS OF PIEZOELECTRIC FLAP IN A CHANNEL USING A STRONGLY COUPLED FEM SCHEME

ME 515 Mechatronics. Overview of Computer based Control System

Table of Contents. Chapter 1. Dielectricity, Piezoelectricity, Pyroelectricity and Ferroelectricity... 1

Temperature Dependence of the Dielectric, Elastic and Piezoelectric Material Constants of Lead Zirconate Titanate Ceramics

Validation of High Displacement Piezoelectric Actuator Finite Element Models

Smart elastomers a touch of robotics

Use of Polyvinylindene Fluoride (PVDF) and Lead Zirconate Titanate (PZT) In Structural Health Monitoring

Dae Won Ji and Sang-Joo Kim

Nonlinear Considerations in Energy Harvesting

Modeling of Composite Piezoelectric Structures with the Finite Volume Method

SENSOR DESIGN FOR PIEZOELECTRIC CANTILEVER BEAM ENERGY HARVESTERS

The New Boundary Condition Effect on The Free Vibration Analysis of Micro-beams Based on The Modified Couple Stress Theory

Piezoline. Introduction Instruction for Using 62 PIEZOLINE. Contents

SENSORS and TRANSDUCERS

Fabrication and performance of d 33 -mode lead-zirconate-titanate (PZT) MEMS accelerometers

Measurements of Radial In-plane Vibration Characteristics of Piezoelectric Disk Transducers

Electric field dependent sound velocity change in Ba 1 x Ca x TiO 3 ferroelectric perovskites

Piezoelectric Fabrics for Energy Harvesting F06-GT05

Design Optimization of Mems Based Piezoelectric Energy Harvester For Low Frequency Applications

ARTICLE IN PRESS. Journal of Sound and Vibration

Ultra large deflection of thin PZT/aluminium cantilever beam

Technical Report PZT-Silicon Cantilever Benders

Piezoelectric Multilayer Beam Bending Actuators

Electrostriction of lead zirconate titanate/polyurethane composites

Analysis of Stress Distribution in Piezoelectric. MEMS Energy Harvester Using Shaped Cantilever Structure

The Pennsylvania State University. The Graduate School. College of Engineering ON THE IMPACT OF ELECTRODE GEOMETRY ON THE PIEZOELECTRIC

Optimisation of interdigitated electrodes for piezoelectric actuators and active fibre composites

Institute for Electron Microscopy and Nanoanalysis Graz Centre for Electron Microscopy

MODELLING OF RECIPROCAL TRANSDUCER SYSTEM ACCOUNTING FOR NONLINEAR CONSTITUTIVE RELATIONS

Ravichetan Dharenni, Ashok M H, Santoshkumar Malipatil

Modelling of Different MEMS Pressure Sensors using COMSOL Multiphysics

Effect of Crystal Orientation on Piezoelectric Response of Single-Crystal Piezoelectric Layers

A TRANSFER MATRIX APPROACH FOR ACOUSTIC ANALYSIS OF A MULTILAYERED ACTIVE ACOUSTIC COATING

VIBRATION CONTROL SIMULATION OF LAMINATED COMPOSITE PLATES WITH INTEGRATED PIEZOELECTRICS

The dynamic method of determination of the piezoelectric hydrostatic coefficients

Piezoactuators. Jiří Tůma

Thermal deformation compensation of a composite beam using piezoelectric actuators

Transcription:

Inversely Prepolarized Piezoceramic Cantilever. D.P. Sedorook and I.V. Ostrovskii Department of Physics and Astronomy, University of Mississippi, Oxford, MS 38677, USA A nonuniform piezoelectric cantilever with enhanced amplitude of vibration at resonances is proposed. The cantilever made of lead zirconate titanate (type of PZT5H) containing inversely poled piezoelectric domains shows an advantageous increase in vibration amplitude when compared to a single domain device. The amplitude of vibrations is enhanced when the domain boundaries are located at the nodes of vibration displacement. The vibration amplitude of a cantilever with 2 or 3 periodically inverted domains is mostly affected at the 2 nd or 3 rd resonance frequency, respectively. The amplitudes are calculated using the Finite lement method (FM). Piezoelectric cantilevers are the key components of microelectromechanical systems with various applications 1, which include atomic force microscopy 2, detection of chemical compounds 3, and the actuation of wing flapping mechanisms in micromechanical robotic insects 4,5. In micromechanical insect prototype applications, rectangular cantilever dimensions are oblong in shape and transverse vibration is utilized 5. Cantilevers that are less oblong in shape may also be good microactuators when the longitudinal displacement is utilized 6. The piezoelectric material, polyvinylidene fluoride (PVDF), has been used to demonstrate this effect to some extent by designing a system with a three layer PVDFshim metal composite laminate, which lead to a PVDF modal actuator prototype 7. Also implemented were two inversely poled PVDFfilms situated sidebyside on the shim metal along the structure width 8. In reference [8], a bimorph PVDF actuator was proposed and calculated by using the Finite lement Method (FM). The FM was also used to calculate proper electrode patterns 1

of a three layer PVDF laminate composite plate 9. We have to note that in the cited publications [7,8,9], a low frequency region of hundreds of Hertz was considered and used in the experiments. However, the actuation effect that has been established could be much further enhanced if a material with a higher piezoelectric coupling coefficient than PVDF is used, for example PZT. Despite it is a bulk material, nowadays socalled nanoceramic can be found on the market place, which opens new perspectives for modern applications. Different types of PZT ceramic are commonly used in cantilevers due to the high piezoelectric coupling of this material. A cantilever normally utilizes the zeroth order vibration mode A 0 with the dominant displacement in the direction perpendicular to the piezoelectric plate. In a plate that is less oblong, the longitudinal mode S 0 may cause the total displacement to be in the longitudinal direction. In this paper, we show that the longitudinal pulsation can be reduced to the advantage of the transverse vibration by introducing a set of inversely poled piezoelectric domains along the length of the plate. A single crystal cantilever and a 2domain cantilever with 2 inversely poled domains are shown in Figure 1(a) and (b), respectively. The Z axis of the PZT5H material PIC181 is perpendicular to the cantilever plate, which vibrates in the ZX plane. Both cantilevers have dimensions 22.4 0.7 mm along the X and Z directions respectively. The white block arrows in Figure 1 represent the polarization of piezoelectric domains. The length of the larger domain in Figure 1(b) is 16.1 mm. The smaller domain is 4.9 mmlong. The interdomain boundary is at x = 17.5 mm because the vertical component of acoustical displacement u z is equal to zero at this point for the second resonance of the single domain and 2domain cantilevers. Two metal electrodes are deposited on the top (z = 0.7 mm) and the bottom (z = 0) surfaces. Normally, the bottom electrode is grounded, while an AC voltage up to 1 kv is applied to the top electrode. The resulting electric field inside the cantilever is well 2

below the coercive field for piezoelectric ceramic 10. The device is designed to operate at a temperature that is lower than the Curie temperature to prevent depolarization of the domains. The multidomain cantilevers with unipolar electrode design having one, two, and three inversely pooled ferroelectric domains are presented in Fig. 1. (1) Z V 1 Domain PIC181 u z 0.7mm u x Ground Polarization X 22.4mm (2) 0.7m Z PIC181 Cantilever with two inversely poled domains and Unipolar lectrode Configuration u z PIC181 P X 16.1mm 4.9mm 1.4mm 3

(3) 0.7mm Z PIC181 Cantilever with three inversely poled domains and Unipolar lectrode Configuration u P X 8.75mm 8.75 1.4mm 2.1m Figure 1. The multidomain cantilevers with the unipolar electrodes on top and bottom surfaces. The directions of the polarization vector are noted by the block arrows. Length is 22.4mm, thickness is 0.7mm, and width is 3 mm (not shown). Displacements u z and u x occur at the tip. Panel (1) presents a standard single crystal cantilever, panel (2) presents twoinverselypoled domain cantilever, and panel (3) presents threeinverselypoled domain cantilever. All top electrodes marked by plus () are potential ones, and all bottom electrodes marked by minus () are grounded. The acoustical displacements are calculated by a Finite lement code that was developed following the approach previously used to calculate dispersion curves for acoustic modes in multidomain ferroelectric plates 11, 12. The FM routine is based on the Hamilton variational principle, where the solution is found by minimizing the time variation of the total energy that 4

consists of four parts: kin 0.5 u 2 dx dz is the kinetic energy; d 0.5 Ddxdz is the energy associated with electric field and electrical displacement D ; st tr 0.5 S Tdxdz is the elastic energy due to stress T and strain S, where superscript tr denotes a transposed matrix; W q dx is work done by the excitation source when the s electric potential V cos 2ft is applied to the electrodes; and q s is the surface charge 0 density on the electrodes. The components of displacement u m (m = x, y, z) satisfy the equations of motion 13 and piezoelectric relations u T ij i (1) x j T x ij j x j c ijkl u x k l 1 n x, m 1 e 1 mij, (2) xm D x m m x m u 1 x, m n1 i S emij mj x j x j 1, (3) where: c is the stiffness tensor under constant electric field, e is the piezoelectric coupling S tensor, is the dielectric permittivity tensor under constant strain, and x, m is the Kronecker symbol, n is domain number 1, 2, and 3. The sign of piezoelectric coupling constants is inverted in the second domain (n = 2) for the Y and Z directions, but remains the same for the X direction. This represents a 180º rotation of the crystallographic coordinate system about the X axis in the inverted domains. The loss may be introduced by adding an imaginary part c ijkl to the stiffness constants c ijkl c ijkl ic ijkl, where c c, factor is some small ijkl ijkl parameter, and S c ijkl are 11=1.063x10 8 C 2 /Nm 2 S, 33 =1.328x10 8 C 2 /Nm 2, e =10.169 C/m 2 31, 5

e 33 =18.662 C/m 2, = 15.269 C/m 2, =21.130x10 10 N/m 2, c =14.430x10 10 N/m 2, e15 c11 12 c13 c33 c44 66 =15.997x10 10 N/m 2, =21.409 x10 10 N/m 2, =3.215 x10 10 N/m 2, c =1.585 x10 10 N/m 2 used in the calculations. The attenuation parameter is related to mechanical quality factor Q = 1/(2 The equations (1) (3) are solved together with the boundary conditions of electric potential amplitude V 0 at the top electrode, zero electric potential at the bottom electrode, and zero acoustical displacement at the leftside edge where the cantilever is clamped. In the FM routine, a triangular numerical mesh with linear approximating functions is used. The solution is tested for convergence by gradually increasing mesh density until the results differed by less than 1%. The results of FMcomputations for the 1 st, 2 nd, and 3 rd vibrational modes/resonances are given in Fig. 2 for three possible models including single domain (a), 2domain (b) and 3domain (c) cantilever. 6

Fig. 2. Vertical displacement u z versus frequency in the multidomain piezoelectric cantilevers shown in Fig. 1, which are made of piezoceramic PIC181 with Q = 100. xcitation by an electric field = 14.3 kv/mm. Plots 1, 2 and 3 correspond to the singledomain, 2domain, and 3domain cantilever, respectively, as shown in Fig. 1, panels (1), (2) and (3). Panel (a) presents displacement amplitudes for the 1 st vibrational mode at 1.09 khz, panel (b) presents displacement amplitudes for the 2 nd vibrational mode at 6.77 khz, and panel (c) presents displacement amplitude for the 3 rd vibrational mode at 18.72 khz. 7

The second type of cantilever configurations may be designed with the help of the bipolar electrodes deposited on a single crystal prepolarized bar, as shown in Fig. 3. Actually, the first panel (1) is the same as in Fig. 1(1), but the panels (2) and (3) are very different. (1) Z V 1 Domain PIC181 Cantilever u z 0.7mm u x Ground Polarization Vector X 22.4mm (2) 0.7m Z Two pairs of Bipolar lectrodes deposited on a single crystal vibrator u z PIC181 P X 16.1mm 4.9 mm 1.4mm 8

(3) 0.7mm Z Three pairs of Bipolar lectrodes deposited on a single crystal vibrator u z P X 8.75mm 8.75 mm 1.4mm 2.1mm Fig. 3. The single crystal cantilevers with bypolar electrodes on top and bottom surfaces. The polarization vectors are along the zaxis. Length is 22.4mm, thickness is 0.7mm, and width is 3 mm (not shown). Displacements u z and u x occur at the tip. Panel (1) presents a standard single crystal cantilever similar to those in Fig. 1(1), panel (2) presents two pairs of metal electrodes with positive and negative electric potential on both top and bottom surfaces, and panel (3) presents three pairs of metal electrodes with positive and negative electric potential on both top and bottom surfaces. In the case of twoelectrode bipolar excitation scheme, each electrode has a 1.4 mm gap at the domain boundary. The bottom electrode below the longer domain is connected to the top electrode above the shorter domain, and the top electrode above the longer domain is 9

connected to the bottom electrode below the shorter domain. When this design is used, the inverse voltage is applied to the smaller domain when compared to the longer domain. Letter P in Fig.3 indicates the probe point, where the displacement component u z is calculated. Fig. 4. Transversal vertical displacement u z versus frequency in the single crystal piezoelectric cantilevers with bipolar electrodes shown in Fig. 3(2), (3). Piezoceramic PIC181, Q = 100. xcitation by an electric field = 14.3 kv/mm. Plot 1 10

corresponds to a single pair of electrodes as shown in Fig. 1(1). Plots 2 and 3 correspond to the two and three pairs of electrodes, respectively, as shown in Fig. 3 (2) and Fig. 3(3). Panel (b) presents displacement amplitudes for the 2 nd vibrational mode at 6.81 khz, panel (c) presents displacement amplitude for the 3 rd vibrational mode at 18.9 khz. The FM computation with the configuration shown in Fig. 3(1) is the same as panel (a) of Fig. 2. However, the calculations for other configurations, Fig. 3(2) and Fig.3(3) are different from those in Fig. 2(b) and Fig. 2(c). These results are given in Fig. 4, panels (b) and (c), respectively. At the first resonance, the only vibration node in the plate is at x = 0, where the cantilever is clamped. In the single domain plate, the u z component of displacement is coherent along the X direction. Two different settings can be used for excitation of the 2 nd mode. The unipolar, when rf voltage is applied to top electrodes, and bottom contacts are grounded, like in Fig. 1(2), and bipolar when rf voltage is applied to the top electrode with the bottom contact grounded in the 1 st domain, and the electrode configuration is flipped in the 2 nd domain, like in Fig. 3(2). The vibration amplitude u z for Q=100 at 2 nd mode, or second resonance, is enhanced for the 2domain cantilever contradictory to the single domain plate. The corresponding resonance frequency for the 2domain plate is 6.8 khz. The 2domain cantilever shows an increase of the vibration amplitude compared to the single domain plate. The amplitude enhancement can be explained by comparing the displacement components at different positions x along the cantilever length. The vibration component has a node at x = 17.5 mm. In the 2domain plate, the domain boundary is placed at this node deliberately to improve the performance of the 2domain cantilever compared to its singledomain counterpart. The 11

vibration amplitude u z at x = 22.4 mm increases in the 2domain cantilever in comparison to the single domain cantilever. Another new configuration with three inversely poled domains as in Fig. 1(3) also gives amplitude increase at higher, 3 rd, vibrational mode, as clear from comparison of Fig. 2 (c) and Fig. 4(c). It can be observed from these graphs that the amplitude u z for the 3 domain cantilever is greater than u z for the single domain cantilevers at 3 rd mode including the three pairs of bipolar electrode configuration, plot 3 in Fig. 4(c). In conclusion: 1) We propose a multidomain cantilever model that has improved performance at selected resonance frequencies. 2) The vibration amplitude in the multidomain cantilever with inversely poled 2 or 3 domains increases at the 2 nd or 3 rd resonance, respectively, if the domain boundaries are located at the nodes of vibration displacements. 3) Based on our calculations for several Q in the range from 10 to 100, the twodomain cantilever shows an advantageous increase in vibration amplitude at 2 nd harmonic resonance when compared to the single domain plate. 4) There is an advantage to multidomain cantilevers in that the number of domains can be optimally selected at appropriate resonance frequencies without changing the size of the device. 5) The applications for these results include microelectromechanical devices, where size may be an important parameter. This work, in part, is made possible due to the research grant Nonlinear vibrations of piezoelectric resonators, UM, 2011212. We also are grateful to Dr. Lucien Cremaldi and Dr. Andriy Nadtochiy for useful discussions. 12

References. 1 S. Banerjee, N. Gayathri, S. Dash, A. K. Tyagi, and B. Raj, A comparative study of contact resonance imaging using atomic force microscopy, Appl. Phys. Lett. 86, 211913 (2005). 2 them rkan Aktakka, Rebecca L. Peterson, and Khalil Najafi. High Stroke and High Deflection BulkPZT Diaphragm and Cantilever Micro Actuators and ffect of PreStress on Device Performance. I Jour. of Microelectromechanical Systems. 14 pgs. (2013) http://sitemaker.umich.edu/aktakka/files/aktakka_2013_jmems.pdf 3 A. B. Nadtochiy, T. K. Hollis, and I. V. Ostrovskii, Ferroelectric bimorph cantilever with selfassembled silane layer, Appl. Phys. Lett. 93, 263503 (2008). 4 R. G. Polcawich, J. S. Pulskamp, S. Bedair, G. Smith, R. Kaul, C. Kroninger,. Wetzel, H. Chandrahalim, and S. A. Bhave, Integrated PiezoMMS actuators and sensors, I Sensors, 2193 (2010). 5 M. Sitti, D. Campolo', J. Yan, R. S. Fearing, T. Su, D. Taylor, and T. D. Sands, Development of PZT and PZNPT Based Unimorph Actuators for Micromechanical Flapping Mechanisms, Proceedings of the I Int. Conf. on Robotics and Automation, 3839 (2001). 6 Z. Cao, J. Zhang, H., Kuwano, I Int. Conf. on Nano/Micro ngineering and Molecular Systems, 716 (2001). 7 C. K. Lee and F. C. Moon, Modal sensors/actuators. Journal of Applied Mechanics 57(6), p. 434441 (1990). 8 WooSeok Hwang and Hyun Chul Park, Finite element modeling of piezoelectric sensors and actuators, AIAA Journal 31(5), p. 930 (1993). 9 JungKyu Ryou, KeunYoung Park, and SeungJo Kim, lectrode Pattern Design of Piezoelectric Sensors and Actuators Using Genetic Algorithms, AIAA Journal 36(2), p. 227 223 (1998). 13

10 L. Pintilie and M. Lisca, Polarization reversal and capacitancevoltage characteristic of epitaxial Pb(Zr,Ti)O 3 layers, Appl. Phys. Lett., 86, 192902 (2005). 11 I. V. Ostrovskii, V. A. Klymko, and A. B. Nadtochiy, Plate wave stopbands in periodically poled lithium niobate, JASA xpress Lett. 125(4), L129 (2009). 12 I. V. Ostrovskii, A. B. Nadtochiy, and V. A. Klymko, Velocity dispersion of plate acoustic waves in a multidomain phononic superlattice, Phys. Rev. B 82(1), 014302 (2010). 13 B. A. Auld, Acoustic fields and waves in solids, (R.. Krieger Publishing Co., Inc.; Malabar, FL 32950, 1990), Vol. 2. 14