Tong University, Shanghai , China Published online: 27 May 2014.

Similar documents
Gilles Bourgeois a, Richard A. Cunjak a, Daniel Caissie a & Nassir El-Jabi b a Science Brunch, Department of Fisheries and Oceans, Box

Ankara, Turkey Published online: 20 Sep 2013.

Published online: 05 Oct 2006.

Dissipation Function in Hyperbolic Thermoelasticity

Nacional de La Pampa, Santa Rosa, La Pampa, Argentina b Instituto de Matemática Aplicada San Luis, Consejo Nacional de Investigaciones Científicas

George L. Fischer a, Thomas R. Moore b c & Robert W. Boyd b a Department of Physics and The Institute of Optics,

Use and Abuse of Regression

Open problems. Christian Berg a a Department of Mathematical Sciences, University of. Copenhagen, Copenhagen, Denmark Published online: 07 Nov 2014.

Precise Large Deviations for Sums of Negatively Dependent Random Variables with Common Long-Tailed Distributions

University, Tempe, Arizona, USA b Department of Mathematics and Statistics, University of New. Mexico, Albuquerque, New Mexico, USA

University, Wuhan, China c College of Physical Science and Technology, Central China Normal. University, Wuhan, China Published online: 25 Apr 2014.

Quan Wang a, Yuanxin Xi a, Zhiliang Wan a, Yanqing Lu a, Yongyuan Zhu a, Yanfeng Chen a & Naiben Ming a a National Laboratory of Solid State

Communications in Algebra Publication details, including instructions for authors and subscription information:

The American Statistician Publication details, including instructions for authors and subscription information:

Testing Goodness-of-Fit for Exponential Distribution Based on Cumulative Residual Entropy

Online publication date: 30 March 2011

CCSM: Cross correlogram spectral matching F. Van Der Meer & W. Bakker Published online: 25 Nov 2010.

Published online: 17 May 2012.

To cite this article: Edward E. Roskam & Jules Ellis (1992) Reaction to Other Commentaries, Multivariate Behavioral Research, 27:2,

Online publication date: 01 March 2010 PLEASE SCROLL DOWN FOR ARTICLE

PLEASE SCROLL DOWN FOR ARTICLE

Effects of mass layer dimension on a finite quartz crystal microbalance

Diatom Research Publication details, including instructions for authors and subscription information:

FB 4, University of Osnabrück, Osnabrück

Dresden, Dresden, Germany Published online: 09 Jan 2009.

Guangzhou, P.R. China

OF SCIENCE AND TECHNOLOGY, TAEJON, KOREA

The Homogeneous Markov System (HMS) as an Elastic Medium. The Three-Dimensional Case

Park, Pennsylvania, USA. Full terms and conditions of use:

Derivation of SPDEs for Correlated Random Walk Transport Models in One and Two Dimensions

Full terms and conditions of use:

Characterizations of Student's t-distribution via regressions of order statistics George P. Yanev a ; M. Ahsanullah b a

Erciyes University, Kayseri, Turkey

Full terms and conditions of use:

Tore Henriksen a & Geir Ulfstein b a Faculty of Law, University of Tromsø, Tromsø, Norway. Available online: 18 Feb 2011

Version of record first published: 01 Sep 2006.

Discussion on Change-Points: From Sequential Detection to Biology and Back by David Siegmund

T. Runka a, M. Kozielski a, M. Drozdowski a & L. Szczepańska b a Institute of Physics, Poznan University of

Geometrical optics and blackbody radiation Pablo BenÍTez ab ; Roland Winston a ;Juan C. Miñano b a

Yan-Qing Lu a, Quan Wang a, Yuan-Xin Xi a, Zhi- Liang Wan a, Xue-Jing Zhang a & Nai-Ben Ming a a National Laboratory of Solid State

G. S. Denisov a, G. V. Gusakova b & A. L. Smolyansky b a Institute of Physics, Leningrad State University, Leningrad, B-

Geometric View of Measurement Errors

Analysis of Leakage Current Mechanisms in BiFeO 3. Thin Films P. Pipinys a ; A. Rimeika a ; V. Lapeika a a

University of Thessaloniki, Thessaloniki, Greece

Online publication date: 22 March 2010

Acyclic, Cyclic and Polycyclic P n

Viscosity Models and Vibration Solutions of Piezoelectric Resonators

PLEASE SCROLL DOWN FOR ARTICLE

The Fourier transform of the unit step function B. L. Burrows a ; D. J. Colwell a a

Global Existence of Large BV Solutions in a Model of Granular Flow

The use of the Normalized Difference Water Index (NDWI) in the delineation of open water features

PLEASE SCROLL DOWN FOR ARTICLE. Full terms and conditions of use:

Xiaojun Yang a a Department of Geography, Florida State University, Tallahassee, FL32306, USA Available online: 22 Feb 2007

Effects of middle plane curvature on vibrations of a thickness-shear mode crystal resonator

An Integrated Analysis of Vibrations of Quartz Crystal Plates with the Mindlin Plate Theory for Resonator Design and Optimization

MSE Performance and Minimax Regret Significance Points for a HPT Estimator when each Individual Regression Coefficient is Estimated

Published online: 10 Apr 2012.

PLEASE SCROLL DOWN FOR ARTICLE

PLEASE SCROLL DOWN FOR ARTICLE

PLEASE SCROLL DOWN FOR ARTICLE

Paper V. Acoustic Radiation Losses in Busbars. J. Meltaus, S. S. Hong, and V. P. Plessky J. Meltaus, S. S. Hong, V. P. Plessky.

Melbourne, Victoria, 3010, Australia. To link to this article:

PLEASE SCROLL DOWN FOR ARTICLE

Open Access Response of Quartz Crystal Microbalance Loaded with Single-drop Liquid in Gas Phase

Projective spaces in flag varieties

Astrophysical Observatory, Smithsonian Institution, PLEASE SCROLL DOWN FOR ARTICLE

First published on: 11 November 2010

of Nebraska - Lincoln

A Strongly Convergent Method for Nonsmooth Convex Minimization in Hilbert Spaces

Raman Research Institute, Bangalore, India

Proceedings Effect of Electrode Configuration on High Temperature Thickness Shear Gallium Phosphate Transducer

Fatigue Crack Analysis on the Bracket of Sanding Nozzle of CRH5 EMU Bogie

Sports Technology Publication details, including instructions for authors and subscription information:

P. C. Wason a & P. N. Johnson-Laird a a Psycholinguistics Research Unit and Department of

93. Study on the vibration characteristics of structural of hydrocyclone based on fluid structure interaction

France. Published online: 23 Apr 2007.

Communications in Algebra Publication details, including instructions for authors and subscription information:

Published online: 04 Oct 2006.

COMSOL for Modelling of STW Devices

Online publication date: 29 April 2011

A Simple Approximate Procedure for Constructing Binomial and Poisson Tolerance Intervals

Published online: 04 Jun 2015.

Finite Element Analysis of Piezoelectric Cantilever

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

To link to this article:

PLEASE SCROLL DOWN FOR ARTICLE. Full terms and conditions of use:

C.K. Li a a Department of Mathematics and

Friction Drive Simulation of a SAW Motor with Slider Surface Texture Variation

William Henderson a & J. Scott McIndoe b a Department of Chemistry, The University of Waikato, Hamilton,

35-959, Rzeszów, Poland b Institute of Computer Science, Jagiellonian University,

PLEASE SCROLL DOWN FOR ARTICLE

Online publication date: 12 January 2010

Anisotropic Hall effect in Al 13 TM 4

Orientation dependence of electromechanical properties of relaxor based ferroelectric single crystals

Noise Shielding Using Acoustic Metamaterials

Published online: 27 Aug 2014.

Silvio Franz a, Claudio Donati b c, Giorgio Parisi c & Sharon C. Glotzer b a The Abdus Salam International Centre for Theoretical

A novel type of transverse surface wave propagating in a layered structure consisting of a piezoelectric layer attached to an elastic half-space

András István Fazekas a b & Éva V. Nagy c a Hungarian Power Companies Ltd., Budapest, Hungary. Available online: 29 Jun 2011

Transcription:

This article was downloaded by: [Shanghai Jiaotong University] On: 29 July 2014, At: 01:51 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 1072954 Registered office: Mortimer House, 37-41 Mortimer Street, London W1T 3JH, UK Ferroelectrics Letters Section Publication details, including instructions for authors and subscription information: http://www.tandfonline.com/loi/gfel20 Effects of Electrode Configuration on Vibration Characteristics of Quartz Thickness-Shear Mode Trapped-Energy Resonators Yan Zhang a & Tao Han a a Department of Instrument Science and Engineering, Shanghai Jiao Tong University, Shanghai 200240, China Published online: 27 May 2014. To cite this article: Yan Zhang & Tao Han (2014) Effects of Electrode Configuration on Vibration Characteristics of Quartz Thickness-Shear Mode Trapped-Energy Resonators, Ferroelectrics Letters Section, 41:1-3, 44-50, DOI: 10.1080/07315171.2014.908686 To link to this article: http://dx.doi.org/10.1080/07315171.2014.908686 PLEASE SCROLL DOWN FOR ARTICLE Taylor & Francis makes every effort to ensure the accuracy of all the information (the Content ) contained in the publications on our platform. However, Taylor & Francis, our agents, and our licensors make no representations or warranties whatsoever as to the accuracy, completeness, or suitability for any purpose of the Content. Any opinions and views expressed in this publication are the opinions and views of the authors, and are not the views of or endorsed by Taylor & Francis. The accuracy of the Content should not be relied upon and should be independently verified with primary sources of information. Taylor and Francis shall not be liable for any losses, actions, claims, proceedings, demands, costs, expenses, damages, and other liabilities whatsoever or howsoever caused arising directly or indirectly in connection with, in relation to or arising out of the use of the Content. This article may be used for research, teaching, and private study purposes. Any substantial or systematic reproduction, redistribution, reselling, loan, sub-licensing, systematic supply, or distribution in any form to anyone is expressly forbidden. Terms & Conditions of access and use can be found at http://www.tandfonline.com/page/termsand-conditions

Ferroelectrics Letters, 41:44 50, 2014 Copyright Taylor & Francis Group, LLC ISSN: 0731-5171 print / 1563-5228 online DOI: 10.1080/07315171.2014.908686 Effects of Electrode Configuration on Vibration Characteristics of Quartz Thickness-Shear Mode Trapped-Energy Resonators YAN ZHANG AND TAO HAN Department of Instrument Science and Engineering, Shanghai Jiao Tong University, Shanghai 200240, China 1. Introduction Communicated by Dr. George W. Taylor (Received in final form December 1, 2013) We study the effects of electrode configuration of quartz thickness-shear mode trappedenergy resonators on their basic vibration characteristics. Free and electrically forced vibration analyses are performed on a rectangular trapped-energy resonator by COM- SOL, focusing on the effects of the extended parts of the electrodes leading from the main rectangular parts of the crystal plate. Results show that the effect of the extended parts on resonant frequencies is of the order of 100 ppm. They also displace the mode center which has further implications in resonator frequency stability. These effects are nonnegligible, and should be taken into consideration in resonator design. Keywords Quartz; thickness-shear mode; trapped-energy; electrode configuration; extended parts; COMSOL; FEM Electrodes are necessary parts of crystal resonators. Their most basic function is for electrically exciting mechanical vibrations. In addition, the motional capacitance or admittance of a resonator is an important design consideration and is calculated from the charge (or current) and voltage on the electrodes in an electrically forced vibration. In plate thicknessshear (TSh) mode trapped-energy resonators [1], the electrodes cover the central portion of a crystal plate only. In these resonators, the electrode mass is responsible for what is called energy trapping [2] through which the vibration is confined mainly under the electrodes and decays rapidly outside them, a behavior crucial to resonator mounting. Recently, due to resonator miniaturization, the effects of electrode configuration have become more important in resonator design than before and have been studied from different angles. Electrodes of varying thickness were shown to be effective in producing strong energy trapping [3 7], which provides a possible future alternative to contoured resonators [8] with variable thickness that are known to be difficult to make. Electrodes of unequal thickness on the top and bottom of a plate resonator were studied in [9 12]. Effects of Corresponding author. E-mail: than@sjtu.edu.cn Color versions of one or more of the figures in the article can be found online at www.tandfonline.com/gfel. 44

Electrode on Quartz Trapped-Energy Resonators 45 M N 2w 2b x 1 c 2a x 3 2l Figure 1. A rectangular trapped-energy resonator. misaligned electrodes were treated in [13]. In [14], it was shown that electrodes with corners cause field concentration and should be avoided in general. Electrode shape is another important factor. At present, circular [15, 16] and rectangular electrodes [17, 18] are routinely used in resonator manufacturing. Elliptical electrodes [19 22] are known to be optimal in certain ways but they are relatively more complicated and are not used often. This paper is concerned with the effects of the details of electrode configuration in trapped-energy resonators (see Fig. 1). The electrodes of these resonators are somewhat irregular. The main parts of them are rectangular or circular, covering the central region of a crystal plate only. They are identical on the plate top and bottom. The regular central part of each electrode also has an extension (see Fig. 1) leading from the central part to the edge of the plate, and thus making the shape of the entire electrodes somewhat irregular. In Figure 1, the central part of the top electrode and its extension to the right edge are filled with dark grey. The extended part of the bottom electrode leading to the left edge is bounded by dotted lines and filled with light grey. In theoretical analyses, the extended parts of the electrodes have to be neglected to simplify the problem. While the extended parts of the electrodes can be included in a finite element numerical analysis [23], they did not seem to have been the focus of existing finite analyses. Therefore, in this paper, we perform a finite element analysis to examine the effects of the fine details of real electrodes with both the regular central parts and the extended parts. Since the TSh vibration is strong under the electrodes and near the electrode edges where the extended parts of the electrodes begin, conceivably the extended parts may potentially affect the basic vibration characteristics of the resonator including resonant frequencies and mode shapes which are related to energy trapping, as well as the admittance of the resonator. 2. Finite Element Model Consider an AT-cut quartz resonator whose length, width and thickness are 2l = 25 mm, 2w = 20 mm, and 2h = 1 mm. The electrode is gold and its dimensions are 2a = 15 mm, 2b = 10 mm, and 2h =0.0025 mm. The width of the extended parts c will be varied in the calculation. The mode we are interested in is the trapped fundamental TSh mode which has half a wavelength along the plate thickness and no more the one wavelength in the the x 3 -x 1 plane. It may seem that one does not need many elements to describe such a mode accurately.

46 Y. Zhang and T. Han Figure 2. Finite element mesh (c=0). However, since this is treated as a three-dimensional problem and the plate has a large length/thickness ratio, a considerable amount of elements are in fact needed. Tetrahedron elements are used in the analysis. The finite element mesh is determined by specifying the number of element layers along the plate thickness and the element size. We varied these two numbers to test for convergence and accuracy, and finally chose twelve layers of elements with a maximal size of 0.35 mm. This givesa total number of 123,432 elements when c = 0. When this mesh is varied a little, the change of the frequency of the fundamental TSh mode, the one that we are most interested in, is only a few Hz among about a couple of MHz. The mesh is shown in Fig. 2 in the case of c = 0. In finite element analysis, the ratios among the sides of the elements should be neither too large nor too small. The elements in Fig. 2 look reasonable. The electrodes have only one layer of elements and they are very thin. This is not a concern because during TSh motions the electrodes deform very little and essentially only their inertia matters. 3. Numerical Results and Discussion When c = 0, we denote the frequency of the fundamental TSh mode of interest by ω 0 which is used as a reference for later comparisons. Through a free vibration analysis it is found that ω 0 = 1.622228 MHz. The distribution of the corresponding major displacement component u 1 is shown in Figure 3. When the plate is vibrating at this mode, u 1 at the top

Electrode on Quartz Trapped-Energy Resonators 47 Figure 3. Distribution of u1 of the fundamental TSh mode (c=0). and bottom of the plate at y =±h are out of phase. u 1 is large in the electrode region and decays rapidly outside the electrodes. This is the so-called energy trapping [2]. To examine the effects of the extended parts of the electrodes described by c, weplot all three displacement components along MN in Figure. 1 on the top surface of the plate at y = h for c = 0, 2, 3, and 4 mm, respectively, in Figures 4(a d). Clearly, u 1 is much larger than the other two components, indicating that the mode is the TSh dominated mode of interest. In Fig. 4(a) when the electrodes are without extended parts, the u 1 displacement looks symmetric about x = 0, which agrees with that in theoretical analyses u 1 is usually described by a cosine function of x. When c is not zero, as c increases in Figures 4 (b d), u 1 gradually loses its symmetry. Its amplitude at the left edge with x = 12.5 mm becomes clearly larger than that at the right edge, indicating that energy trapping has been affected and the left edge can no longer be considered as well trapped in Fig. 4(d). It is also visible, although barely, in Fig. 4(d) that the peak of u 1 in fact is no long at x = 0 and has displaced a little to the left. This seemingly small displacement of the peak or mode center has important implications in resonator frequency stability [24]. A perfectly symmetric resonator with a symmetric mode has zero normal acceleration sensitivity. Any displacement of the mode center destroys the symmetry and results in a nonzero normal acceleration sensitivity. This causes concern as the U.S. military is trying to improve the relative acceleration sensitivity of the current technology from 10 10 /g to 10 12 /g in the near future. The resonant frequencies of the fundamental TSh mode for different values of c are summarized in Table 1, along with the maximal values of the real and imaginary parts of the admittance (conductance and susceptance). As c increases, the extended parts of the electrodes become wider. This results in more electrode inertia and lowers the resonant frequencies as expected. This effect if of the order of 100 ppm and is relevant in resonator performance. As the electrode width increases from 0 to 3 mm, the admittance increases too as expected because bigger electrodes collect more charges. However, this trend stops

48 Y. Zhang and T. Han Figure 4. Effect of c on the displacement distributions along MN on the top of the plate (a) c=0, (b) c=2 mm, (c) c=3 mm, (d) c=4 mm. when c goes from 3 to 4 mm, which may be due to that the mode center displacement is in the opposite direction of the electrode extension. The maximal values of the conductance and susceptance occur at slowlychanging frequencies as c varies. This can be seen in Fig. 5. 4. Conclusions COMSOL can calculate quartz resonator frequencies as accurate as a few ppm. This is sufficient for analyzing electrode configuration details. Resonator electrodes are irregular, Table 1 Effect of c on the frequency and admittance of the fundamental TSh mode c (mm) 0 2 3 4 Frequency ω (MHz) ω 0 = 1.622228 1.62213 1.622045 1.621937 Frequency shift ω ω 0 (Hz) 98 198 291 ω ω 0 /ω 0 (ppm) 60 122 179 Max. conductance (S) 0.00637405 0.00677934 0.00688084 0.00681428 Max. susceptance (S) 0.00346394 0.00354615 0.00357256 0.00355788

Electrode on Quartz Trapped-Energy Resonators 49 Figure 5. Effect of c on the conductance and susceptance of the resonator (a) c=0, (b) c=2 mm, (c) c=3 mm, (d) c=4 mm. with extended parts which usually destroy the symmetry of the resonator structures. The extended parts of the electrodes have a frequency effect of the order of 100 ppm, and cause mode center displacement which has implications in resonator frequency stability. These effects indicate that in resonator design attention should also be paid to the extended parts of the electrodes. An electrode configuration that preserves the symmetry of the structure will be ideal. The admittance is not very sensitive to the extended parts of the electrodes. Funding This work was supported by NSFC with Grant No.11174205 and Program for New Century Excellent Talents in University (NCET-12-0357). References 1. H. F. Tiersten, Analysis of trapped-energy resonators operating in overtones of coupled thickness shear and thickness twist. J. Acoust. Soc. Am. 59(4), 879 888 (1976). 2. V. E. Bottom, Introduction to Quartz Crystal Unit Design. New York: Van Nostrand Reinhold, (1982). 3. S. Y. Pao, R. Huang, and C. S. Lam, Energy trapping effect of stepped electrode in miniature ATcut quartz resonator. in Piezoelectricity, Acoustics Waves and Device Applications. Singapore: World Scientific. 152 154 (2007).

50 Y. Zhang and T. Han 4. J. S. Yang, H. Xue, H. Y. Fang, Y. T. Hu, J. Wang, and L. J. Shen, Effects of electrodes with varying thickness on energy trapping in thickness-shear quartz resonators. IEEE Trans. Ultrason. Ferroelect. Freq. Contr. 54(4), 892 895 (2007). 5. J. Wang, L. J. Shen, and J. S. Yang, Effects of electrodes with continuously varying thickness on energy trapping in thickness-shear mode quartz resonators. Ultrasonics. 48(2), 150 154 (2008). 6. J. S. Yang, Z. G. Chen, and H. P. Hu, Electrically forced vibration of a thickness-twist mode piezoelectric resonator with non-uniform electrodes. Acta Mechanica Solida Sinica. 20(3), 266 274 (2007). 7. J. S. Yang, Z. G. Chen, and Y. T. Hu, Vibration of a thickness-twist mode piezoelectric resonator with asymmetric, non-uniform electrodes. IEEE Trans. Ultrason. Ferroelect. Freq. Contr. 55(4), 841 848 (2008). 8. P. Li, F. Jin, and J. S. Yang, Thickness-shear vibration of an AT-cut quartz resonator with a hyperbolic contour. IEEE Trans. Ultrason. Ferroelect. Freq. Contr. 59(5), 1006 1012 (2012). 9. F. Boersma and E. C. van Ballegooyen, Rotated Y-cut quartz crystal with two different electrodes treated as a one-dimensional acoustic composite resonator. J. Acoust. Soc. Am. 62(2), 335 340 (1977). 10. J. A. Kosinski, Thickness vibrations of flat piezoelectric plates with massy electrodes of unequal thickness. in Proceedings of IEEE Ultrasonics Symposium. 70 73 (2003). 11. J. Wang and L. J. Shen, Exact thickness-shear resonance frequency of electroded piezoelectric crystal plates. J. Zhejiang Univ. SCI. 6A(9), 980 985 (2005). 12. J. S. Yang, H. G. Zhou, and W. P. Zhang, Thickness-shear vibration of rotated Y-cut quartz plates with relatively thick electrodes of unequal thickness. IEEE Trans. on Ultrasonics, Ferroelectrics, and Frequency Control. 52(5), 918 922 (2005). 13. H. J. He, J. X. Liu, and J. S. Yang, Effects of mismatched electrodes on an AT-cut quartz resonator. IEEE Trans. Ultrason. Ferroelect. Freq. Contr.59(2), 281 286 (2012). 14. J. S. Yang, H. Xue, and Y. T. Hu, Finite element analysis of stress field concentration near the edge of an electrode. Ferroelectrics Letters Section. 34(3 4), 108 112 (2007). 15. H. F. Tiersten, and Y. S. Zhou, Transversely varying thickness modes in quartz resonators with beveled cylindrical edges. J. Appl. Phys. 76(11), 7201 7208 (1994). 16. H. J. He, J. S. Yang, and Q. Jiang, Thickness-shear and thickness-twist vibrations of circular AT-cut quartz resonators. Acta Mechanica Solida Sinica. accepted. 17. H. F. Tiersten, B. J. Lwo, and B. Dulmet, Transversely varying thickness modes in trapped energy resonators with shallow and beveled contours. J. Appl. Phys. 80(2), 1037 1046 (1996). 18. H. J. He, J. S. Yang, J. A. Kosinski, and J. Wang, Thickness-shear vibration of a rectangular quartz plate with partial electrodes. Acta Mechanica Solida Sinica. 26(2), 121 128 (2013). 19. R. D. Mindlin, Optimal sizes and shapes of electrodes for quartz resonators J. Acoust. Soc. Am.. 43(6), 1329 1331, (1968). 20. Z. T. Yang and J. S. Yang, Optimal electrode shape and size of a few singly rotated quartz and langasite resonators. IEEE Trans. Ultrason. Ferroelect. Freq. Contr. 56(2), 237 238, (2009). 21. Z. T. Yang, J. S. Yang, and Y. T. Hu, Optimal electrode shape and size of doubly rotated quartz plate thickness mode piezoelectric resonators. Appl. Phys. Lett. 92(10), article number 103516 (2008). 22. T. F. Ma, C. Zhang, X. N. Jiang, and G. P. Feng, Thickness-shear mode quartz crystal resonators with optimized elliptical electrodes. Chinese Physics B. 20(4), art. no. 047701 (2011). 23. L. Yang, N. Vichev, and Z. P. Yu, Modal analysis of practical quartz resonators using finite element resonators. IEEE Trans. Ultrason. Ferroelect. Freq. Contr.57(2), 292 298 (2010). 24. Y. S. Zhou and H. F. Tiersten, On the normal acceleration sensitivity of contoured quartzresonators with the mode shape displaced with respect to rectangular supports. J. Appl. Phys. 69(5), 2862 2870 (1991).