Piezoelectric composite lattices with high sensitivity

Similar documents
Cellular solid structures with unbounded thermal expansion. Roderic Lakes. Journal of Materials Science Letters, 15, (1996).

Analysis of high loss viscoelastic composites

ELECTROMECHANICAL RESPONSE OF PIEZOELECTRIC FOAMS

PIEZOELECTRIC TECHNOLOGY PRIMER

Giant enhancement in effective piezoelectric sensitivity by pyroelectric coupling

Reciprocity failure in piezoelectric polymer composite

Three-dimensional stiff cellular structures with negative Poisson's ratio

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

Properties of a chiral honeycomb with a Poisson's ratio -1

Lecture 4 Honeycombs Notes, 3.054

INTRODUCTION TO PIEZO TRANSDUCERS

Chapter 2 Surface Acoustic Wave Motor Modeling and Motion Control

FINITE ELEMENT ANALYSIS OF EFFECTIVE MECHANICAL PROPERTIES, VIBRATION AND ACOUSTIC PERFORMANCE OF AUXETIC CHIRAL CORE SANDWICH STRUCTURES

Finite Element Analysis of Piezoelectric Cantilever

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

Reduced warp in torsion of reticulated foam due to Cosserat elasticity: experiment

Fig. 1. Circular fiber and interphase between the fiber and the matrix.

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

Design of a fastener based on negative Poisson's ratio foam adapted from

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

FOAM TOPOLOGY BENDING VERSUS STRETCHING DOMINATED ARCHITECTURES

Material parameters for electrostriction

A comprehensive numerical homogenisation technique for calculating effective coefficients of uniaxial piezoelectric fibre composites

Mechanics PhD Preliminary Spring 2017

PHASE-FIELD SIMULATION OF DOMAIN STRUCTURE EVOLUTION IN FERROELECTRIC THIN FILMS

Mechanical Engineering Ph.D. Preliminary Qualifying Examination Solid Mechanics February 25, 2002

H Kalathur 1, 2 1, 2,3. , R S Lakes. 3

Extreme Cosserat elastic cube structure with large magnitude of negative Poisson s ratio

PHYSICAL REVIEW B 71,

Piezoelectric Composites as Bender Actuators

Stability of Cosserat solids: size effects, ellipticity and waves

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

Module 7: Micromechanics Lecture 29: Background of Concentric Cylinder Assemblage Model. Introduction. The Lecture Contains

MCE603: Interfacing and Control of Mechatronic Systems

Thickness Optimization of a Piezoelectric Converter for Energy Harvesting

Modeling and analysis of the electromechanical behavior of surface-bonded piezoelectric actuators using finite element method

Computational Modelling of Vibrations Transmission Loss of Auxetic Lattice Structure

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

Piezoelectric Resonators ME 2082

Piezoelectric Vibration Energy Harvesting. Characteristics of Barium Titanate Laminates

Stresses in Curved Beam

3.032 Problem Set 2 Solutions Fall 2007 Due: Start of Lecture,

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

A Note on Suhir s Solution of Thermal Stresses for a Die-Substrate Assembly

ENGN 2340 Final Project Report. Optimization of Mechanical Isotropy of Soft Network Material

16.21 Techniques of Structural Analysis and Design Spring 2003 Unit #5 - Constitutive Equations

Level 7 Postgraduate Diploma in Engineering Computational mechanics using finite element method

Exercise: concepts from chapter 8

DAMPING CONTROL OF A PZT MULTILAYER VIBRATION USING NEGATIVE IMPEDANCE CIRCUIT

Numerical analyses of cement-based piezoelectric smart composites

Transducer design simulation using finite element method

XI. NANOMECHANICS OF GRAPHENE

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

Effect of fibre shape on transverse thermal conductivity of unidirectional composites

PHYSICAL PROPERTIES OF CRYSTALS

STRESS STRAIN AND DEFORMATION OF SOLIDS, STATES OF STRESS

Passive Damping Characteristics of Carbon Epoxy Composite Plates

HEALTH MONITORING OF PLATE STRUCTURE USING PIEZO ELECTRIC PATCHES AND CURVATURE MODE SHAPE

Solid State Physics (condensed matter): FERROELECTRICS

The strain response of silicone dielectric elastomer actuators

A short review of continuum mechanics

Optimal piezo-electro-mechanical coupling to control plate vibrations

Composite materials: mechanical properties

Mechanics of Inflatable Fabric Beams

EE C247B / ME C218 INTRODUCTION TO MEMS DESIGN SPRING 2014 C. Nguyen PROBLEM SET #4

Curvature of a Cantilever Beam Subjected to an Equi-Biaxial Bending Moment. P. Krulevitch G. C. Johnson

EE 5344 Introduction to MEMS CHAPTER 6 Mechanical Sensors. 1. Position Displacement x, θ 2. Velocity, speed Kinematic

Lecture 19. Measurement of Solid-Mechanical Quantities (Chapter 8) Measuring Strain Measuring Displacement Measuring Linear Velocity

Micromechanical analysis of FRP hybrid composite lamina for in-plane transverse loading

THE TOPOLOGICAL DESIGN OF MATERIALS WITH SPECIFIED THERMAL EXPANSION USING A LEVEL SET-BASED PARAMETERIZATION METHOD

Small-Scale Effect on the Static Deflection of a Clamped Graphene Sheet

7.4 The Elementary Beam Theory

6. Bending CHAPTER OBJECTIVES

46th AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics, and Materials Conference April 2005 Austin, Texas

Structural Health Monitoring Using Smart Piezoelectric Material

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

The fracture toughness of planar lattices: imperfection sensitivity

Bending of Simply Supported Isotropic and Composite Laminate Plates

Computational Analysis for Composites

Macroscopic Elastic Constitutive Relationship of Cast-in-Place Hollow-Core Slabs

Analysis of a ceramic bimorph piezoelectric gyroscope

Understand basic stress-strain response of engineering materials.

A Vertical Electrostatic Actuator with Extended Digital Range via Tailored Topology

Testing and analysis of high frequency electroelastic characteristics of piezoelectric transformers

Elastic and piezoelectric fields in substrates GaAs 001 and GaAs 111 due to a buried quantum dot

Module 7: Micromechanics Lecture 34: Self Consistent, Mori -Tanaka and Halpin -Tsai Models. Introduction. The Lecture Contains. Self Consistent Method

Exploring Piezoelectric Properties of Wood and Related Issues in Mathematical Description. Igor Dobovšek

SENSOR DESIGN FOR PIEZOELECTRIC CANTILEVER BEAM ENERGY HARVESTERS

THE BENDING STIFFNESSES OF CORRUGATED BOARD

S.-W. Ricky Lee, M.-L. Zhu & H.L. Wong Department of Mechanical Engineering, Hong Kong University of Science & Technology, Hong Kong

Buckling Analysis of Isotropic Circular Plate with Attached Annular Piezoceramic Plate

Highly piezoelectric, thermal stable ferroelectrets from cyclic olefin copolymer. Yan Li, Hui Wang, Changchun Zeng ANTEC 2015, Orlando

Negative incremental bulk modulus in foams

CHAPTER THREE SYMMETRIC BENDING OF CIRCLE PLATES

CHAPTER -6- BENDING Part -1-

Application of piezoelectric actuators to active control of composite spherical caps

3. BEAMS: STRAIN, STRESS, DEFLECTIONS

TABLE OF CONTENTS. Mechanics of Composite Materials, Second Edition Autar K Kaw University of South Florida, Tampa, USA

Advanced Structural Analysis EGF Section Properties and Bending

EE C245 ME C218 Introduction to MEMS Design

Transcription:

Piezoelectric composite lattices with high sensitivity Roderic Lakes Department of Engineering Physics, Engineering Mechanics Program, Department of Materials Science, Rheology Research Center, University of Wisconsin, 1500 Engineering Drive, Madison, WI 53706-1687 February 5, 2014 Abstract Lattice structures are presented with a piezoelectric sensitivity, d, much larger in magnitude than that of material comprising the lattice ribs. Large sensitivity is achieved by using piezoelectric bimorph elements as ribs; these bend in response to electric input, giving rise to a much larger displacement than for axial elements. Complex electrical connectivity is helpful but not necessary; surface contact can be sufficient. Lattices are amenable to piezoelectric vibration damping. Piezoelectric sensitivity can be arbitrarily large and is unbounded. adapted from Lakes, R. S., Piezoelectric composite lattices with high sensitivity, Philosophical Magazine Letters, 94, (1), 37-44 (2014). 1 1 lakes@engr.wisc.edu 1

1 Introduction 1.1 Piezoelectric materials Piezoelectric materials produce an electric polarization when stressed mechanically and deform in response to an electric field. They are always anisotropic. Specifically [1], the strain ɛ ij, depends upon stress σ kl via the elastic compliance J ijkl, upon electric field E k in piezoelectric materials with modulus tensor d kij at constant temperature, and on temperature change T (via the thermal expansion α ij ). Moreover the electric displacement vector D i depends on electric field via K ij which is the dielectric tensor at constant stress and temperature, and it can depend on temperature change T via the pyroelectric effect; p i is the pyroelectric coefficient at constant stress: ɛ ij = J ijkl σ kl + d kij E k + α ij T (1) D i = d ijk σ jk + K ij E j + p i T (2) If there are phase angles, material properties associated with response to sinusoidal input in time become complex quantities in the context of viscoelasticity and dielectric relaxation [2]. Properties are expressed in the reduced notation in which the two indices associated with stress or strain are reduced to one, so d 333 d 33. The most sensitive piezoelectric materials are ceramics for which the sensitivity of charge to force d 33 is from 100 to 600 pc/n, in comparison with 10 to 20 pc/n, for piezoelectric polymers. The sensitivity of displacement to input voltage has the same value by virtue of reciprocity, provided the material is elastic; the units pm/volt are equivalent to pc/n. Piezoelectric composite materials have been conceptualized with particulate, fibrous, platelet microstructure or they may have a more complex geometry. A set of rule-of-mixtures type equations [3] was provided to approximate the effective piezoelectric coefficients of piezoelectric composites of simple structure. For d 31 it is a Voigt equation; for d 33 the piezoelectric coefficients are weighted by the dielectric constant as well as the volume fraction. Fibrous piezoelectric composites have also been analyzed [4]. Piezoelectric composites have advantages over homogeneous piezoelectric ceramic materials. For example, they may offer increased electromechanical coupling efficiency, also a decrease in specific acoustic impedance, and suppression of unwanted modes of vibration [5]. Too, the voltage coefficient g 33 (open circuit electric field divided by stress) of a laminar composite of piezoelectric and compliant dielectric layers can substantially exceed that of the piezoelectric constituent alone [3]. 1.2 Bounds on properties Bounds on various physical properties have been developed [6]. For example, the upper and lower bounds of stiffness of two phase and multi-phase composites have been obtained by Hashin [7], and by Hashin and Shtrikman [8]. Further analysis on composites is reviewed by Hashin [9] and Milton [6]. Bounds on thermoelastic properties of two phase composites [10] are weighted averages of the expansion values of the constituents. Bounds similar to elastic Voigt and Reuss bounds and to Hashin-Shtrikman bounds have been determined for piezoelectric composites [11] [12]. Arbitrarily high positive or negative thermal expansions can be achieved in lattice type composites with void content [13] [14] or in dense composites with interfaces that allow slip. Composites in which the inclusions are voids or cells, filled with air or liquid, are called cellular solids. Lattices, including two-dimensional honeycombs, are cellular solids in which the cells have a regular repeating structure. Zero expansion can be attained in such lattices [15]. These composites contain 2

rib elements that themselves have composite microstructure, corresponding to a low level of structural hierarchy [16]. Each rib element is a bi-layer made of two bonded layers of differing thermal expansion coefficient. Such materials can exceed the limits imposed by bounds that tacitly assume there is no void space and no slip between constituents. In the present research, we develop piezoelectric lattices based on composite rib elements. We find that piezoelectric sensitivity d can be arbitrarily greater than that of material in the ribs. 2 Lattice structures and properties 2.1 Lattice sensitivity analysis Several triangular piezoelectric lattices are shown in Figure 1. Rib elements that undergo bending in response to an electrical signal are chosen as piezoelectric elements. Each rib of length L and full thickness h in the lattice consists of two layers of piezoelectric material, arranged as a bimorph, a kind of bender element. When a voltage is applied to such a rib, one layer expands, and the other contracts, giving rise to bending. In the figures, dark and light shading of the the ribs denotes different polarization directions of layers in the bimorph ribs. Some electrical connections from +, - electrodes are indicated schematically by thin lines. Commercial bimorph strips in our laboratory have electrodes near the ends, i. e. the shortest sides; the remaining surfaces are covered with an insulating layer; end connection is not shown in the figures to avoid clutter with the end contact. Figure 1(a) shows a fully bend dominated lattice in which each rib is free to bend in response to an electrical signal. A fully triangulated lattice with curved ribs shown in Figure 1(b) approaches stretch dominated behavior if the curvature is small. Specifically the elastic (Young s) modulus E normalized to the modulus E s of the solid comprising the ribs, is E/E s [ h L ]. This goes as the inverse first power of the rib slenderness corresponding to greater stiffness than bend dominated structures, as shown for a bi-layer thermoelastic lattice [15]. Rib curvature is needed to achieve thermoelastic or piezoelectric response: change in curvature results in overall dimensional changes in the lattice. Ribs of sufficient curvature give rise to bend dominated behavior E/E s [ h L ]3. Straight ribs are of interest in that available bimorph elements can be used, therefore such lattices are analyzed. For the edge contact lattice in Figure 1(a), the ribs are free to bend unless the lattice is subject to a constraint. Stiff spacers (small rectangles) transmit the bending motion among adjacent cells. If the spacers are electrically conductive, the outer surfaces, if conductive, of the cells are at a common potential that can be imposed by a voltage V applied at the surface, but wiring is still needed to impose a difference in potential on the inner surfaces. The deflection per volt for a single bimorph cantilever rib of length L and full thickness h containing two anti-parallel piezoelectric layers of sensitivity d 31 and thickness h/2 is [17] (u/v ) rib = 3 2 d 31[ L h ]2. Elastic structural compliance (displacement / force) is u/f = 2 w S 11[ L h ]3 with w as the width and S 11 as the material compliance, inverse of Young s modulus E. A bender in which one layer is piezoelectric and the other layer is not piezoelectric, is called a unimorph; the analysis is more complicated [18] than for a bimorph. The lattice of bimorph ribs contains multiple cells of angle θ between ribs: n cells thick in the z direction, m cells thick in the x direction. Each of the n layers of cells has two horizontal ribs that contribute fully 1 4 (u/v ) rib, because the center of a rib moves 1/4 as much as the end of a cantilever. Each layer also has two oblique ribs that contribute 1 4 (u/v ) ribsin( θ 2 ). So the sensitivity is as follows: d eff 3 33 = d 31 2 2n 4 (1 + sin(θ 2 ))[L h ]2. (3) 3

(a) (b) Figure 1: Triangular piezoelectric lattices. (a), fully bend-dominated lattice. (b), lattice with curved ribs. The sign of d eff 33 depends on the rib orientation; interchange of light and dark layers results in a reversal of sign. This sensitivity is much larger in magnitude than the intrinsic piezoelectric coefficient of the material from which the ribs are made, owing to the enhancement of displacement by bending. For the lattice segment shown, n = 2 and θ = 60. Slender ribs give rise to enhanced piezoelectric sensitivity. The in-plane elastic modulus E of this lattice normalized to the modulus cos( θ 2 ) of the solid rib material is, neglecting the spacer compliance, E/E s = 1 2 1+sin( θ )[ h L ]3. The transverse 2 sensitivity has contributions from m cells in the horizontal direction. In each cell there two ribs that contribute 1 4 (u/v ) ribcos( θ 3 2m 2 ). So deff 31 = d 31 2 4 cos( θ 2 )[ L h ]2. A piezoelectric lattice based on two dimensional chiral honeycomb [19] is shown in Figure 2. This was originally devised for its elastic properties: a negative Poisson s ratio that tends to the lower bound -1, independent of strain [19] as ribs become more slender; slender beam theory is used throughout. For the corresponding piezoelectric lattice, reversal in bimorph orientation is incorporated in each rib to provide for the appropriate deformation shape of the rib, a sigmoid curve. Electrical insulation must therefore be provided at the rib midpoints. The elastic modulus of this type of chiral honeycomb is E/E s = 3[ h L c ] 3 [ Lc r ]2, with r as the node radius, L c as the 4

Figure 2: Chiral piezoelectric lattice. rib length and E s as the solid rib modulus. Deformation is bend-dominated; the circular nodes impose bending moments on the ribs as the lattice is compressed or stretched. Elastic properties are isotropic in-plane provided the ribs are of equal length so the angle between them is θ = 60 as shown. The in-plane Poisson s ratio approaches -1 for slender ribs. The displacement u of the rib-node attachment point of a cell in terms of the bend angle is u = rφ; bend angle is between the deformed rib and its initial slope at its attachment to the node. Bend angle can result from moments from external forces [19] or from electrical input. For the latter case, the angle per volt for a single bimorph rib of length L is φ/v = 3d 31 [ L ] [17]. The strain is ɛ = φr/r with R = ( Lc h 2 r )2 + 4 as the distance between node centers. The lattice sensitivity is d eff 33 = uy V = 1 V nɛl ccos( θ 2 ) because n cells are assumed in the y direction and the projection of length in the y direction is taken. The effective piezoelectric sensitivity of the lattice is as follows with L c = 2L: d eff 3 33 = d 31 2 n cos( θ 2 ) ( Lc r )2 + 4 [ L c h ]2 (4) The piezoelectric sensitivity d increases as the square of the rib slenderness; the elastic modulus decreases as the cube of rib slenderness as with the edge contact lattice in Figure 1(a). Chiral lattices are bend dominated yet they allow an independent control over stiffness and piezoelectric sensitivity via the node radius r. As r becomes small in comparison with rib length, the elastic modulus increases as the inverse square of r; the piezoelectric sensitivity decreases linearly with r. The axial strain in all directions is equal because this structure easily accommodates area change but resists shape change, so if the dimension in each direction is the same, d eff 31 = d eff 33. 2.2 Continuum view of lattices 2.2.1 Equivalent continuum and applications As for the continuum view and scaling, there is no length scale in the classical theory of elasticity or piezoelectricity, so lattices may be of any size, from macroscopic to nano-scale. The above 5

Figure 3: Piezoelectric lattice that allows a voltage to be applied at the surface. Electrically conductive spacers (small rectangles) provide an electrical and mechanical path through the lattice. piezoelectric lattices differ from elastic or thermoelastic lattices in that they should not be viewed as equivalent continua. Elastic lattices or foams deform in response to stress at the surface, and thermoelastic lattices deform in response to overall temperature changes. The above piezoelectric lattices, by contrast, receive electrical input via connectivity to appropriate ribs. That results in a dependence of effective sensitivity on the number of cells in a particular direction, in contrast to elastic and thermoelastic lattices. Similarly, stacks of axial piezoelectric elements have electrical connectivity such that a low voltage applied results in a sum of displacements of each element. If an equivalent continuum amenable to surface electrode contact is desired, a structure such as that in Figure 3 may be considered. The spacers here can be electrically conductive or they may have sufficiently high dielectric constant to provide a much lower impedance than the piezoelectric elements. Suppose there are n layers in the z direction. The overall displacement is the sum of the displacement of each layer; the imposed voltage is shared among the layers so there is no dependence of sensitivity on cell number. The sensitivity is d eff 3 33 = d 31 2 [ L h ]2. Sensitivity is much greater than that of an axial stack, due to bending displacement. All the motion is in the z direction so d eff 31 = 0 in contrast to the prior lattices. The void content of this structure, in contrast to the prior ones, is controllable independently of the rib slenderness, by varying the vertical thickness of the spacers. Void content can be made arbitrarily small. Macroscopic piezoelectric lattices may be assembled from commercial bimorph or unimorph bender elements; for macroscopic or microscopic lattices, prototype methods previously used for composites [20] are applicable. The lattices presented here are distinct from piezoelectric superlattices [21] which are layered structures on a fine scale that exhibit dispersion of waves [22]. Applications are envisaged in the context of actuators, energy harvesting, and structures that change shape in response to stimuli. Airplane wings that change shape have been considered in the context of chiral honeycomb [23]. Vibration damping in structures is enabled by viscoelastic damping that occurs in the piezoelectric material and interface layers, and this may be incorporated into lattices via the concept of Forward [24]. Specifically, appropriately matched electrical resistance connected to the circuit provides mechanical damping via the electromechanical coupling that is inherent in piezoelectric materials. Such damping may be useful in satellites and other aerospace applications in which vibration damping is beneficial. The captured energy may also be harvested. 6

2.2.2 Two-phase bounds on effective properties and how to exceed them The piezoelectric sensitivity is unbounded for these lattices; arbitrarily large values can be attained via sufficiently slender ribs. Thin ribs give rise to reduced stiffness. Too, these lattices are bend dominated so they are much more compliant than stretch-dominated structures that are optimized for stiffness. As in the case of thermal expansion, mathematical bounds for two phase composites do not apply to lattices because the bound analysis tacitly assumes the two phases are perfectly attached together with no slip or void space. Such an assumption, often not stated explicitly, does not apply to the lattice structures considered here. The lattice concept can be readily extended to three dimensions. For example, cubical cells can be envisaged with spacers similar to those used to separate triangles in Figure 1(a). The chiral lattice allows some recovery of stiffness if nodes of small size are used. In the limit of zero node size, the lattice becomes fully triangulated and stretch dominated, with no piezoelectric sensitivity; the bending analysis does not apply to that case. 3 Conclusions To conclude, piezoelectric lattices based on bimorph ribs can have a piezoelectric d sensitivity much larger in magnitude than that of material comprising the lattice ribs. Piezoelectric sensitivity is unbounded. High sensitivity is achieved using ribs that bend in response to an input voltage. Chiral lattices with Poisson s ratio approching -1 in plane provide additional control over the sensitivity and modulus. Acknowledgement Support from the MRSEC program and from the Petroleum Research Fund is gratefully acknowledged. References [1] J. F. Nye, Physical Properties of Crystals, Oxford, Clarendon, 1976. [2] R. S. Lakes, Viscoelastic Materials, Cambridge University Press, 2009. [3] E. Newnham, D. P. Skinner, and L. E. Cross, Mater. Res. Bull. 13(1978) p. 525. [4] K. Schulgasser, J. Mech. Phys. Solids 40 (1992) p. 473. [5] W. A. Smith, IEEE Proc. Ultrasonic Symposium (1989) p. 755. [6] G. Milton, The Theory of Composites, Cambridge University Press, 2002. [7] Z. Hashin, J. Appl. Mech., Trans. ASME 84E (1962) p. 143. [8] Z. Hashin, and S. Shtrikman, J. Mechs. Phys. Solids 11 (1963) p. 127. [9] Z. Hashin, J. Appl. Mech. 50 (1983) p. 481. [10] J. L. Cribb, Nature 220 (1968) p. 576. [11] P. Bisegna, R. Luciano, J. Mechs. Phys. Solids 44(1996) p. 583. 7

[12] J. Y. Li, M. L. Dunn, Phil. Mag. A, 81 (2001) p. 903. [13] R. S. Lakes, J. Materials Sci. Lett. 15 (1996) p. 475. [14] R. S. Lakes, Appl. Phys. Lett. 90 (2007) 221905. [15] J. Lehman and R. S. Lakes, J. Intelligent Material Sys. Struct., 23 (2012) p. 1263. [16] R. S. Lakes, Nature 361 (1993) p. 511. [17] J. G. Smits and S. I. Dalke, IEEE Proc. 1989 Ultrason. Symp. (1989) p. 781. [18] J. G. Smits and W. S. Choi, IEEE Trans. Ultrason., Ferroelectrics and Frequency Control 38 (1991) p. 256. [19] D. Prall, and R. S. Lakes, Int. J. Mechanical Sci. 39 (1996) p. 305. [20] A. Safari, M. Allahverdi, E. K. Akdogan, J. Materials Sci. 41 (2006) p. 177. [21] J. L. Bleustein, Appl. Phys. Lett. 13 (1968) p. 412. [22] L. Fernandez, V.R. Velasco, Surface Sci. 185 (1987) p. 175. [23] J. Martin, J. Heyder-Bruckner, C. Remillat, F. Scarpa, K. Potter, M. Ruzzene, Phys. Stat. Sol. (b) 245 (2008) p. 1521. [24] R. L. Forward, Appl. Opt. 18, (1979) p. 690. 8