Polynomial approach for modeling a piezoelectric disc resonator partially covered with electrodes

The frequency spectrum of a partially metallized piezoelectric disc resonator was studied using Legendre polynomials. The formulation, based on three-dimensional equations of linear elasticity, takes into account the high contrast between the electroded and non-electroded regions. The mechanical displacement components and the electrical potential were expanded in a double series of orthonormal functions and were introduced into the equations governing wave propagation in piezoelectric media. The boundary and continuity conditions were automatically incorporated into the equations of motion by assuming position-dependent physical material constants or delta-functions. The incorporation of electrical sources is illustrated. Structure symmetry was used to reduce the number of unknowns. The vibration characteristics of the piezoelectric discs were analyzed using a three-dimensional modelling approach with modal and harmonic analyses. The numerical results are presented as resonance and anti-resonance frequencies, electric input admittance, electromechanical coupling coefficient and field profiles of fully and partially metallized PIC151 and PZT5A resonator discs. In order to validate our model, the results obtained were compared with those published previously and those obtained using an analytical approach.     

Exact dynamic solutions to piezoelectric smart beams including peel stresses I: Theory and application

This work presents exact dynamic solutions to piezoelectric (PZT) smart beams including peel stresses. The governing equations of partial differential forms are firstly derived for a PZT smart beam made of the identical adherends, and then general solutions of the governing equations are studied. The analytical solutions are applied to a cantilever beam with a partially bonded PZT patch to the fixed end. For the given boundary conditions, exact solutions of the steady state motions are obtained. Based on the exact solutions, frequency spectra, natural frequencies, normal mode shapes, harmonic responses of the shear and peel stresses are discussed for the PZT actuator. The details of the numerical results and sensing electric charges will be presented in Part II of this work. The exact dynamic solutions can be directly applied to a PZT bimorph bender. To compare with the classic shear lag model whose numerical demonstrations will be given in Part II, the related equations are also derived for the shear lag rod model and shear lag beam model.     

Electric admittance analysis of quartz crystal resonator in thickness-shear mode induced by array of surface viscoelastic micro-beams

The electric admittance of a compound system composed of a thicknessshear mode(TSM) quartz crystal resonator(QCR) and an array of surface viscoelastic micro-beams(MBs) is studied.The governing equations of the MBs are derived from the Timoshenko-beam theory in consideration of shear deformation.The electrical admittance is described directly in terms of the physical properties of the surface epoxy resin(SU-8) MBs from an electrically forced vibration analysis.It is found that both the inertia effect and the constraint effect of the MBs produce competitive influence on the resonant frequency and admittance of the compound QCR system.By further comparing the numerical results calculated from the Timoshenko-beam model with those from the Euler-beam model,the shear deformation is found to lead to some deviation of an admittance spectrum.The deviations are revealed to be evident around the admittance peak(s) and reach the maximum when a natural frequency of the MBs is identical to the fundamental frequency of the QCR.Besides,a higher order vibration mode of the MBs corresponds to a larger deviation at the resonance.     

AN ANALYSIS OF A CYLINDRICAL THIN SHELL AS A PIEZOELECTRIC TRANSFORMER

The vibration of a circular cylindrical piezoelectric ceramic shell with tangential poling is investigated for transformer application and an analytical solution is obtained.Numer- ical calculations of the output voltage,input admittance and the efficiency of the transformer are presented.The results indicate that the structure analyzed can be used as a piezoelectric transformer with a certain harmonic driving frequency and matching load impedance.     

Limited permeable crack moving along the interface of a piezoelectric bi-material

A plane problem for a crack moving with a subsonic speed along the interface of two piezoelectric semi-infinite spaces is considered. The crack is assumed to be free from mechanical loading. The limited permeable electric condition with an account of electric traction is adopted at its faces. A uniformly distributed mixed mode mechanical loading and an electric flux are prescribed at infinity. The problem is reduced to the Riemann–Hilbert problem by means of introducing a moving coordinate system and assuming that the electric flux is uniformly distributed along the crack region. An exact solution of this problem is proposed. It permits to find in closed form all necessary electromechanical characteristics at the interface and to formulate the equation for the determination of the electric flux. Analysis of this equation confirms the correctness of the assumption concerning the uniform distribution of the electric flux in the crack region. The values of the electric flux are determined by solving the obtained equation. Thereafter, the stress and electric intensity factors as well as their asymptotic fields at the crack tip are also found. The particular case of a crack moving in a homogeneous piezoelectric material is considered. The values of the electric flux and the fracture parameters are found exactly in a simple form for this case. Also, a numerical analysis is performed for a crack propagating with a subsonic speed between PZT4 and PZT5 materials and for a crack moving in PZT4 material. The electric flux in the crack region, stress and electric intensity factors, crack opening and the energy release rate (ERR) are found as functions of the crack speed, loading and electric permeability of the crack medium. The influence of the electric traction on the crack faces upon the mentioned parameters is demonstrated.     

相似电路耦合模型及其在压电—梁结构分析中的应用

根据线性系统微分方程的相拟性确定机械系统的相似电路,由压电材料的本构和动力方法确定压电材料作为电耦合的变压器,从而用一具耦合的电路系统模拟智能结构的电学和力学行为,通过一个压电－梁智能结构的分析,证明了模型的可用性。     

Modeling and analysis of electrostatic MEMS filters

We present an analytical model and closed-form expressions describing the response of a tunable MEMS filter made of two electrostatic resonators coupled by a weak microbeam. The model accounts for the filter geometric and electric nonlinearities as well as the coupling between them. It is obtained by discretizing the distributed-parameter system to produce a reduced-order model. We predict the filter deflection and static pull-in voltage by solving a boundary-value problem (BVP). We also solve an eigenvalue problem (EVP) to determine the filter poles (the natural frequencies delineating the filter bandwidth). We found a good agreement between the results obtained using our model and published experimental results. We found that, when the input and output resonators are mismatched, the first mode is localized in the softer resonator whereas the second mode is localized in the stiffer resonator. We demonstrated that mismatch between the resonators can be countermanded by applying different DC voltages to the resonators. As the effective nonlinearities of the filter grow, multi-valued responses appear and distort the filter performance. Once again, we found that the filter can be tuned to operate linearly by choosing a DC voltage that makes the effective nonlinearities vanish.     

基于传递函数法的兰姆波解析模拟

本文研究了兰姆波传播的解析模型。薄压电传感器（PZT）被用于兰姆波的激发。为模拟PZT的激励,采用PZT边缘的分布式点源近似等效PZT的激励,建立了解析模型的边界条件。结合傅里叶变换,在频率波数域求解了三维波动方程,得到了在单位周期激励下兰姆波的频域响应。结合传递函数法,可以得到任意激励下兰姆波的频域响应。利用傅里叶逆变换,进一步求解了兰姆波在时间-空间域的全局波场。为了验证解析模型,将解析模拟与实验进行了比对。试验测量值与解析模拟结果相吻合。     

Exact static solutions to piezoelectric smart beams including peel stresses. II. Numerical results,comparison and discussion

This part presents the numerical results, comparisons and discussion for the exact static solutions of smart beams with piezoelectric (PZT) actuators and sensors including peel stresses presented in Part I. (International Journal of Solids and Structures, 39, 4677–4695) The actuated stress distributions in the adhesive and the adhesive edge stresses varying with the thickness ratios are firstly obtained and presented. The actuated internal stress resultants and displacements in the host beam are then calculated and compared with those predicted by using the shear lag model. The stresses in the adhesive caused by an applied axial force, bending moment and shear force are calculated, and then used to compute the sensing electric charges for comparison with those predicted using the shear lag model. The numerical results are given for the smart beam with (a) one bonded PZT and (b) two symmetrically bonded PZTs, with a comparison to those predicted using the shear lag model. Novel, simple and more accurate formulas for the equivalent force and bending moment induced by applied electric field are also derived for the host beam with one PZT or two symmetrically bonded PZTs. The symmetric shear stress and the anti-symmetric peel stress components caused by a shear force are discussed. In addition, in the case of PZT edge debonding, the stress redistribution in the adhesive and the self-arresting mechanism are also investigated.     

基于新型裂尖杂交元的压电材料断裂力学研究

提出了一种裂尖邻域杂交元模型,将其与标准杂交应力元结合来求解压电材料裂纹尖 端的奇性电弹场和断裂参数的数值解．裂纹尖端杂交元的建立步骤为：1) 利用高次内插有限元特征法求解特征问题,得到反映裂尖奇异性电弹场状况的特 征值和特征角分布函数;2) 利用广义Hellinger-Reissner变分泛函以及特征问题的解来建立裂尖邻域杂交元模型．该 方法求解电弹场时,摒弃了传统有限元方法中裂尖奇异性场需要借助解析解的做法,也避免 了单纯有限元方法中需要在裂尖端部进行高密度单元划分．采用PZT5板中心裂纹问题 作为考核例,数值结果显示了良好的精确性．作为进一步应用,求解了含中心界面裂纹 的PZT4-PZT5两相压电材料的应力强度因子和电位移强度因子．所有的算例都考虑 了3种裂纹面电边界条件．     

Propagation and localization of two-dimensional in-plane elastic waves in randomly disordered layered piezoelectric phononic crystals

Two-dimensional in-plane wave propagation and localization in the disordered layered piezoelectric phononic crystals with material 6 mm are investigated taking the electromechanical coupling into account. The electric field is approximated as quasi-static. The analytical solutions of elastic waves are obtained. The 6 × 6 transfer matrix between two consecutive unit cells is obtained by means of the mechanical and electrical continuity conditions. The expressions of the localization factor and localization length in the disordered periodic structures are presented by regarding the variables of the mechanical and electrical fields as the elements of the state vector. The numerical results of the localization factors and localization lengths are presented for two kinds of disordered piezoelectric phononic crystals, i.e. ZnO–PZT–5H and PVDF–PZT–5H piezocomposites. It is seen from the results that the incident angle of elastic waves and the thickness of the piezoelectric ceramics have significant effects on the wave localization characteristics. For different piezoelectric phononic crystals, the effects of the incident angle are very different. Moreover, with the increase of the disorder degree, the localization phenomenon is strengthened.     

Reflection of Pulse Waves and Resonance Characteristics of Arterial Beds

The problem of axisymmetric wave flow of a viscous incompressible fluid in a system consisting of a long thin deformable tube and a terminal element that determines the conditions of wave reflection at the tube end is analyzed. An expression for the input admittance of the system is obtained and the dependence of the admittance on the system parameters is investigated. The resonant frequencies at which the admittance amplitude has extrema are found and it is shown that at these frequencies the admittance variations with variation of the terminal-element parameters are maximal. The dependence of the resonant frequency on the tube length is investigated. Possible applications of the results obtained to the hydromechanical interpretation of a novel method of pulse diagnostics are discussed.     

Simulating electric activities of neurons by using PSPICE

Electric activities in the Morris–Lecar neuron and Josephson junction coupled resonator are investigated in a numerical way, and electric circuits are also designed by using the Personal Simulation Program with Integrated Circuit Emphasis (PSPICE). Within the improved Morris–Lecar circuit, a new integrator for the ion channel of potassium is designed, and the transition of electric activities, quiescent state to spiking to bursting to quiescent state could be observed. In the circuit of the Josephson-junction coupled resonator, an equivalent circuit is designed to reproduce several types of electric activity. The detailed parameter regions are detected to generate spiking and bursting states in the electric circuits for neurons, and these results are consistent with the numerical results. Bifurcation diagrams for interspike interval (ISI) vs. the forcing current are calculated to detect the excitability of the neuron model.     

Development, modeling and deflection analysis of hybrid micro actuator with integrated thermal and piezoelectric actuation

Micro actuators are irreplaceable part of motion control in minimized systems. The current study presents an analytical model for a new Hybrid Thermo Piezoelectric micro actuator based on the combination of piezoelectric and thermal actuation mechanisms. The micro actuator structure is a double PZT cantilever beam consisting of two arms with different lengths. The presented micro actuator uses the structure of electrothermal micro actuator in which polysilicon material is replaced by PZT. Also the voltage and poling directions are considered in the lengthwise of PZT beams. As a result, the piezoelectric actuation mechanism is based on d ₃₃ strain coefficient. The tip deflection of micro actuator is obtained using Timoshenko beam theory. Analytical results are compared with FEM results along with other reported results in the literature. The effects of geometrical parameters and PZT material constants on actuator tip deflection are studied to provide an efficient optimization of HTP micro actuator.     

铁电陶瓷PZT53复杂力电耦合行为的实验研究

通过实验研究了平行和垂直于极化方向的正应力对铁电陶瓷锆钛酸铅(PZT53)的电滞回线(E3-P3)和电致应变曲线(E-ε)的影响. 实验发现平行于极化方向的压应力对PZT53陶瓷的电滞回线、电致应变曲线形状以及矫顽场大小都有明显的影响,但是垂直于极化方向的拉、压应力只对PZT53陶瓷的电致应变曲线形状有明显的影响,但对电滞回线形状和矫顽场大小都没有显著影响. 采用畴翻转的模型详细解释了观察到的实验现象,所得结果为建立铁电陶瓷的多轴力、电耦合本构模型,提供了物理基础.     

The modes of an open resonator with plane circular mirrors

Summary A new method is developed to calculate the frequencies, amplitudes, and energy losses of the modes of an open resonator consisting of two plane circular mirrors.The present theory does not suffer from some of the disadvantages of earlier theories. Numerical results are given for a resonator with a small Fresnel number.     

A two-dimensional closed-form solution for the free-vibrations analysis of piezoelectric sandwich plates

This work presents a two-dimensional (2D) closed-form solution for the free-vibrations analysis of simply-supported piezoelectric sandwich plates. It has the originality to consider all components of the electric field and displacement, thus satisfying exactly the electric equilibrium equation. Besides, the formulation considers full layerwise first-order shear deformation theory and through-thickness quadratic electric potential. Its independent mechanical and electric variables are decomposed using Fourier series expansions, then substituted in the derived mechanical and electric 2D equations of motion. The resulting eigenvalue system is then condensed so that only nine mechanical unknowns are retained. After its validation on single- and three-layer piezoelectric, and hybrid sandwich plates, the present approach was then used to analyze thickness modes of a square sandwich plate with piezoceramic faces and elastic cross-ply composite core. It was found that only the first three thickness modes are global, thus can be modeled by the mixed equivalent single-layer/layerwise approach, often retained in the literature; the remaining higher thickness modes being characteristic of sandwich behavior; i.e., dominated by the deformations of either the core or the faces. These results, together with presented through-thickness variations of the mechanical and electric variables clearly recommend full layerwise modeling. Several numerical results are provided for future reference for validation of 2D approximate analytical or numerical approaches; in particular, of 2D piezoelectric adaptive finite elements.     

Piezoelectric composites with periodic multi-coated inhomogeneities

A new, robust homogenization scheme for determination of the effective properties of a periodic piezoelectric composite with general multi-coated inhomogeneities is developed. In this scheme the coating does not have to be thin, the shape and orientation of the inclusion and coatings do not have to be identical, their centers do not have to coincide, their properties do not have to remain uniform, and the microstructure can be with the 2D elliptic or the 3D ellipsoidal inclusions. The development starts from the local electromechanical equivalent inclusion principle through the introduction of the position-dependent equivalent eigenstrain and electric field. Then with a Fourier series expansion and a superposition procedure, the volume-averaged equivalent eigenstrain and electric field for each phase are obtained. The results in turn are used in an energy equivalent criterion to determine the effective properties of the composite. In this model the interphase interactions in each multi-coated particle and the long-range interactions between the periodically distributed particles are fully accounted for. To demonstrate its wide range of applicability, we applied it to examine the properties of several periodic composites: (i) piezoelectric PZT spherical particles in a polymer matrix, (ii) continuous glassy fibers with thin PZT coating in an epoxy matrix, (iii) spherical PZT particles coated by thick or functionally graded piezoelectric layer, (iv) spheroidal voids coated with a thick non-piezoelectric layer in a PZT matrix, and (v) spherical piezoelectric inhomogeneities with eccentric, non-uniform thickness coating. The calculated results reflect the complex nature of interplay between the properties of core, matrix, and coating, as well as whether the coating is uniform, functionally graded, or eccentric. The accuracy of this new scheme is checked against the double-inclusion and other micromechanics models, and good agreement is observed.     

VALIDATION OF NUMERICAL SHALLOW WATER MODELS FOR STRATIFIED SEICHES

A new analytical solution is presented for the case of a stratified seiche. This solution, especially its energetics, is useful for the validation of numerical shallow water models under stratified conditions. The utility of the analytical solution for validation is shown by using it to validate a simple finite difference numerical model. A comparison of the energetics of the numerical and analytical solutions reveals that the model results converge rapidly to the analytical solution with increasing resolution, such that a grid size of 30×30 would appear adequate for validation. In addition to properly resolving the spatial features, good temporal resolution is also necessary for validation, i.e use of a Courant number (Cr) less than one. For example, owing to the numerical dispersion of the present model, using Cr=5/4 rather than Cr =1/4 for the 50×50 grid resulted in 3·6 times larger RMS errors of model versus analytical barotropic available potential energy. This new analytical solution should be applied to a test suite of such validation tools before using such numerical models to simulate the more realistic geophysical flows encountered in lakes, bays, harbours and semi-enclosed seas under stratified conditions. © 1997 by John Wiley & Sons, Ltd.     

Suppression of chaos in a MEMS resonator by delayed position feedback

A delayed position feedback control is applied on DC voltage source for suppressing chaos of a typical MEMS resonator actuated by electrostatic forces. A theoretical necessary condition for chaotic oscillation of the controlled system is presented. Numerical results and the analytical prediction reveal the evolution of dynamical behavior of the system with AC voltage amplitude and the control effect of delayed feedback on reducing chaos of the system. It shows that the delayed feedback control is effective on suppressing chaos of the micro mechanical resonator.