A coupled electromechanical analysis of a piezoelectric layer bonded to an elastic substrate: Part I,development of governing equations

This two-part contribution presents a novel and efficient method to analyze the two-dimensional (2-D) electromechanical fields of a piezoelectric layer bonded to an elastic substrate, which takes into account the fully coupled electromechanical behavior. In Part I, Hellinger–Reissner variational principle for elasticity is extended to electromechanical problems of the bimaterial, and is utilized to obtain the governing equations for the problems concerned. The 2-D electromechanical field quantities in the piezoelectric layer are expanded in the thickness-coordinate with seven one-dimensional (1-D) unknown functions. Such an expansion satisfies exactly the mechanical equilibrium equations, Gauss law, the constitutive equations, two of the three displacement–strain relations as well as one of the two electric field-electric potential relations. For the substrate the fundamental solutions of a half-plane subjected to a vertical or horizontal concentrated force on the surface are used. Two differential equations and two singular integro-differential equations of four unknown functions, the axial force, N, the moment, M, the average and the first moment of electric displacement, D₀ and D₁, as well as the associated boundary conditions have been derived rigorously from the stationary conditions of Hellinger–Reissner variational functional. In contrast to the thin film/substrate theory that ignores the interfacial normal stress the present one can predict both the interfacial shear and normal stresses, the latter one is believed to control the delamination initiation.     

The effect of bonding layer properties on the dynamic behaviour of surface-bonded piezoelectric sensors

The dynamic behaviour of piezoelectric sensors depends on the bonding condition along the interface between the sensors and the host structure. This paper provides a comprehensive theoretical study of the effect of the bonding layer on the coupled electromechanical characteristics of a piezoelectric sensor bonded to an elastic substrate, which is subjected to a high frequency elastic wave. A sensor model with a viscoelastic bonding layer, which undergoes a shear deformation, is proposed to simulate the two dimensional electromechanical behaviour of the integrated system. Analytical solution of the problem is provided by using Fourier transform and solving the resulting integral equations in terms of the interfacial stress. Numerical simulation is conducted to study the effect of the bonding layer upon the dynamic response of the sensor under different loading frequencies. The results indicate that the modulus and the thickness of the bonding layer have significant effects on sensor response, but the viscosity of the bonding layer is relatively less important.     

Antiplane mechanical and inplane electric time-dependent load applied to two coplanar cracks in piezoelectric ceramic material

This work is concerned with the dynamic response of two coplanar cracks in a piezoelectric ceramic under antiplane mechanical and inplane electric time-dependent load. The cracks are assumed to act either as an insulator or as a conductor. Laplace and Fourier transforms are used to reduce the mixed boundary value problems to Cauchy-type singular integral equations in Laplace transform domain. A numerical Laplace inversion algorithm is used to determine the dynamic stress and electric displacement factors that depend on time and geometry. A normalized equivalent parameter describing the ratio of the equivalent magnitude of electric load to that of mechanical load is introduced in the numerical computation of the dynamic stress intensity factor (DSIF) which has a similar trend as that for the pure elastic material. The results show that the dynamic electric field will impede or enhance crack propagation in a piezoelectric ceramic material at different stages of the dynamic electromechanical load. Moreover, the electromechanical response is greatly affected by the ratio of the crack length to the ligament between the cracks. The stress and electric displacement intensity factor can be combined by the energy density factor or function to address the fracture of piezoelectric materials under the combined influence of electromechanical loading.     

Effect of initial stress on the propagation behavior of Love waves in a layered piezoelectric structure

The propagation behavior of Love waves in a layered piezoelectric structure with an initial stress is investigated in this article. It involves a thin piezoelectric layer bonded perfectly to an elastic substrate. Solutions of the mechanical displacement and electrical potential function are obtained for the piezoelectric layer and elastic substrate by solving the coupled electromechanical field equations. The phase velocity equations of the Love wave propagation and the stress fields in the layered piezoelectric structure are obtained for electrical open and short cases on the free surface, respectively. The effect of the initial stress on the phase velocity, the stress fields and the coupled electromechanical factor are discussed, respectively. Three sets of piezoelectric layer–elastic substrate systems are considered, i.e. BaTiO₃ ceramic layer–borosilicate glass substrate, PZT-5H ceramic layer–borosilicate glass substrate, and PZT-5H ceramic layer–SiO₂ glass substrate. It is seen that the phase velocity of the Love wave propagation decreases with the increase of the magnitude of the initial stress. The coupled electromechanical factor increases remarkably, as the magnitude of the initial the stress is greater than 100 MPa. This is useful for the design of acoustic surface wave devices.     

Multimode vibration control using several piezoelectric transducers shunted with a multiterminal network

In this paper a new approach is presented to reduce vibrations for one- and two-dimensional mechanical structures, as beam or thin plates, by means of several piezoelectric transducers shunted with a proper electric network system. The governing equations of the whole system are coupled to each other through the direct and converse piezoelectric effect. More in detail, the mechanical equations are expressed in accordance with the modal theory considering n vibration modes and the electrical equations reduce to the one-dimensional charge equation of electrostatics for each of n considered piezoelectric transducers. In this electromechanical system, a shunting electric device forms an electric subsystem working as multi degrees of freedom (dof’s) damped vibration absorber for the mechanical subsystem. Herein, it is introduced a proper transformation of the electric coordinates in order to approximate the governing equations for the whole shunted system with n uncoupled, single mode piezoelectric shunting systems that can be readily damped by the methods reported in literature. A further numerical optimisation problem on the spatial distribution of the piezoelectric elements allows to achieve a better performance. Numerical case studies of two relevant systems, a double clamped beam and a fully clamped plate, allow to take into account issues relative to the proposed approach. Laboratory experiments carried out in real time on a beam clamped at both ends consent to validate the proposed technique.     

The complete synchronization condition in a network of piezoelectric micro-beams

This work deals with the dynamics of a network of piezoelectric micro-beams (a stack of disks). The complete synchronization condition for this class of chaotic nonlinear electromechanical system with nearest-neighbor diffusive coupling is studied. The nonlinearities within the devices studied here are in both the electrical and mechanical components. The investigation is made for the case of a large number of coupled discrete piezoelectric disks. The problem of chaos synchronization is described and converted into the analysis of the stability of the system via its differential equations. We show that the complete synchronization of N identical coupled nonlinear chaotic systems having shift invariant coupling schemes can be calculated from the synchronization of two of them. According to analytical, semi-analytical predictions and numerical calculations, the transition boundaries for chaos synchronization state in the coupled system are determined as a function of the increasing number of oscillators.

     

Transient response of piezoelectric composites caused by the sudden formation of localized defects

A continuum model is presented which is capable of generating the transient electroelastic field in piezoelectric composites of periodic microstructure, caused by the sudden appearance of localized defects. These defects are simulated by associating to every one of the ten piezoelectric parameters of the constituents a distinct damage variable. This procedure enables the modeling of localized cracks, soft and stiff inclusions and cavities. As a result, the constitutive equations of the piezoelectric phases appear in a specific form that includes eigen-electromechanical field variables which represent these defects. The method of solution is based on the combination of two distinct approach. In the first one, the representative cell method is employed according to which the periodic composite, which is discretized into several cells, is reduced to a problem of a single cell in the discrete Fourier transform domain. The resulting coupled elastodynamic and electric equations, initial, boundary and interfacial conditions in the transform domain are solved by employing a wave propagation in piezoelectric composite analysis which forms the second approach. The method of solution is verified by comparison with an analytical solution for the transient response of a piezoelectric material with a semi-infinite mode III-crack. Several applications are presented for the sudden formation of cracks in homogeneous and layered piezoelectric materials which are subjected to various types of electromechanical loading, and for the sudden appearance of a cavity. The effect of electromechanical coupling on the dynamic response is discussed.     

精化杂交压电固体单元的研究

基于力、电耦合问题的三类交量广义交分原理，提出了广义杂交压电单元列式。为了进一步改进单元的性能和保证单元能够通过分片检验，通过引入非协调模式、放松电学方程约束条件和单元间的弱连续性条件，建立了新的、修正的广义交分原理，在此基础上成功地引入了应力、应交的正交化插值模式，从而建立了精化杂交压电单元法，它继承了常规精化杂交单元的全部优点。文中所推导的八节点精化杂交压电固体单元列式完全避免了矩阵求逆运算，较广义杂交压电单元和杂交应力压电单元均显著提高了计算效率。数值算例表明，与同类型其他单元相比，该单元明显具有更好的对歪斜网格的适应性。     

压电涂层-功能梯度压电界面层-基底结构的二维接触问题研究

在多层压电元件中,由于界面处材料成分和性质的突变,常常导致界面处应力集中,使得界面处出现开裂或蠕变现象,从而大大缩短了压电元件的使用寿命。功能梯度压电材料作为界面层,可有效的缓解界面材料不匹配导致的破坏。本文主要研究利用功能梯度压电材料界面层连接压电涂层和基底,分析三层结构在圆柱型压头作用下的力电响应。利用傅里叶积分变换技术,本文将压电涂层-功能梯度压电层-基底结构在刚性圆柱压头作用下的二维平面应变接触问题转化为带有柯西核的奇异积分方程。运用高斯-切比雪夫积分公式,将奇异积分方程转化为线性方程组并对其进行数值求解,得到压电涂层-功能梯度压电层-基底结构在圆柱形压头作用下的应力分布和电位移分布。数值结果表明,梯度压电材料参数的变化对结构中的力电响应具有重要的影响。本文研究结果对于利用功能梯度压电界面层消除界面处的应力不连续导致的界面破坏具有重要的理论指导意义,研究结果可为功能梯度压电材料界面层的设计提供帮助。     

Analytical modeling of sandwich beam for piezoelectric bender elements

Piezoelectric bender elements are widely used as electromechanical sensors and actuators.An analytical sandwich beam model for piezoelectric bender elements was developed based on the first-order shear deformation theory(FSDT),which assumes a single rotation angle for the whole cross-section and a quadratic distribution function for coupled electric potential in piezoelectric layers,and corrects the effect of transverse shear strain on the electric displacement integration.Free vibration analysis of simply- supported bender elements was carried out and the numerical results showed that,solu- tions of the present model for various thickness-to-length ratios are compared well with the exact two-dimensional solutions,which presents an efficient and accurate model for analyzing dynamic electromechanical responses of bender elements.     

Electroelastic field of a piezoelectric annular finite cylinder

This paper presents a theoretical study of a piezoelectric annular cylinder under axisymmteric electromechanical loading. The piezoelectric material is assumed to be transversely isotropic and the general solutions of the governing equations are obtained in terms of a Fourier–Bessel series containing Bessel functions of the first and second kind. The boundary-value problems for vertical pressure and an electric charge loading applied to the ends of an annular cylinder are solved by expanding the applied loading in terms of a Fourier–Bessel series. Selected numerical results for the electroelastic field of an annular cylinder are presented for different aspect ratios of a cylinder and material properties.     

曲梁压电俘能器强迫振动的格林函数解

本文运用格林函数法求解了曲梁压电俘能器在强迫振动下的解析解.运用微分法分析了压电层合曲梁结构面内各内力,根据曲梁压电 俘能器的动力学方程组,基于压电本构关系,建立了包含径向阻尼但不考虑俘能器曲梁结构部分的轴向力以及轴向惯性项的Prescott力 电耦合模型. 采用Laplace变换法求得了耦合振动方程的格林函数解.根据叠加原理和格林函数的物理意义,对耦合的系统方程解耦进而 求得强迫振动下曲梁压电俘能器的输出电压. 数值计算中,通过与现有文献的解析解进行对比,验证了本文解析解的有效性,并研究了阻 尼、电阻等重要物理参数对压电函数和谐振频率的影响.通过与有关传统直梁压电俘能器研究成果的对比,体现了曲梁压 电俘能器Prescott模型的高效集能特性. 数值分析研究表明:(1)使得曲梁俘能器达到最大输出电压时连接的最优负载电 阻为1 M$\Omega$;(2)通过更换适当的基底材料,降低材料的弹性模量,可以改变曲梁俘能器的高基频现象,以使结构适应 更复杂的工作环境,但这会导致俘能器的工作效率降低.      

A coupled refined high-order global-local theory and finite element model for static electromechanical response of smart multilayered/sandwich beams

In the present study, a coupled refined high-order global-local theory is developed for predicting fully coupled behavior of smart multilayered/sandwich beams under electromechanical conditions. The proposed theory considers effects of transverse normal stress and transverse flexibility which is important for beams including soft cores or beams with drastic material properties changes through depth. Effects of induced transverse normal strains through the piezoelectric layers are also included in this study. In the presence of non-zero in-plane electric field component, all the kinematic and stress continuity conditions are satisfied at layer interfaces. In addition, for the first time, conditions of non-zero shear and normal tractions are satisfied even while the bottom or the top layer of the beam is piezoelectric. A combination of polynomial and exponential expressions with a layerwise term containing first order differentiation of electrical unknowns is used to introduce the in-plane displacement field. Also, the transverse displacement field is formulated utilizing a combination of continuous piecewise fourth-order polynomial with a layerwise representation of electrical unknowns. Finally, a quadratic electric potential is used across the thickness of each piezoelectric layer. It is worthy to note that in the proposed shear locking-free finite element formulation, the number of mechanical unknowns is independent of the number of layers. Excellent correlation has been found between the results obtained from the proposed formulation for thin and thick piezoelectric beams with those resulted from the three-dimensional theory of piezoelasticity. Moreover, the proposed finite element model is computationally economic.     

Transient response of a piezoelectric ceramic strip with an eccentric crack under electromechanical impacts

The transient response of a piezoelectric strip with an eccentric crack normal to the strip boundaries under applied electromechanical impacts is considered. By using the Laplace transform, the mixed initial-boundary-value problem is reduced to triple series equations, then to a singular integral equation of the first kind by introducing an auxiliary function. The Lobatto–Chebyshev collocation technique is adopted to solve numerically the resulting singular integral equation. Dynamic field intensity factors and energy release rate are obtained for both a permeable crack and an impermeable crack. The effects of the crack position and the material properties on the dynamic stress intensity factor are examined and numerical results are presented graphically.     

A phenomenological and extended continuum approach for modelling non-equilibrium flows

This paper presents a new technique that combines Grad’s 13-moment equations (G13) with a phenomenological approach to rarefied gas flows. This combination and the proposed solution technique capture some important non-equilibrium phenomena that appear in the early continuum-transition flow regime. In contrast to the fully coupled 13-moment equation set, a significant advantage of the present solution technique is that it does not require extra boundary conditions explicitly; Grad’s equations for viscous stress and heat flux are used as constitutive relations for the conservation equations instead of being solved as equations of transport. The relative computational cost of this novel technique is low in comparison to other methods, such as fully coupled solutions involving many moments or discrete methods. In this study, the proposed numerical procedure is tested on a planar Couette flow case, and the results are compared to predictions obtained from the direct simulation Monte Carlo method. This test case highlights the presence of normal viscous stresses and tangential heat fluxes that arise from non-equilibrium phenomena, which cannot be captured by the Navier–Stokes–Fourier constitutive equations or phenomenological modifications.      

Analysis of an interfacial crack in a piezoelectric bi-material via the extended Green’s functions and displacement discontinuity method

Based on the extended Stroh formalism, we first derive the extended Green’s functions for an extended dislocation and displacement discontinuity located at the interface of a piezoelectric bi-material. These include Green’s functions of the extended dislocation, displacement discontinuities within a finite interval and the concentrated displacement discontinuities, all on the interface. The Green’s functions are then applied to obtain the integro-differential equation governing the interfacial crack. To eliminate the oscillating singularities associated with the delta function in the Green’s functions, we represent the delta function in terms of the Gaussian distribution function. In so doing, the integro-differential equation is reduced to a standard integral equation for the interfacial crack problem in piezoelectric bi-material with the extended displacement discontinuities being the unknowns. A simple numerical approach is also proposed to solve the integral equation for the displacement discontinuities, along with the asymptotic expressions of the extended intensity factors and J-integral in terms of the discontinuities near the crack tip. In numerical examples, the effect of the Gaussian parameter on the numerical results is discussed, and the influence of different extended loadings on the interfacial crack behaviors is further investigated.     

Permeable cracks between two dissimilar piezoelectric materials

An interfacial crack with electrically permeable surfaces between two dissimilar piezoelectric ceramics under electromechanical loading is investigated. An exact expression for singular stress and electric fields near the tip of a permeable crack between two dissimilar anisotropic piezoelectric media are obtained. The interfacial crack-tip fields are shown to consist of both an inverse square root singularity and a pair of oscillatory singularities. It is found that the singular fields near the permeable interfacial crack tip are uniquely characterized by the real valued stress intensity factors proposed in this paper. The energy release rate is obtained in terms of the stress intensity factors. The exact solution of stress and electric fields for a finite interfacial crack problem is also derived.     

初应力对压电层状结构声表面波传播性能的影响

研究了压电层状结构中初应力对广义Ｒａｙｌｅｉｇｈ波传播相速度和机电耦合性能的影响,通过求解含初应力的运动微分方程,对自由界面电学开路和短路两种情况得到了相应的相速度方程。给出了具体的数值算例,所得结果对于提高和改善声表面波器件性能有参考意义。     

磁电热弹耦合材料三维多裂纹超奇异积分法

用超奇异积分方程法将多场耦合载荷作用下磁电热弹耦合材料内含任意形状和位置三维多裂纹问题转化为求解一以广义位移间断为未知函数的超奇异积分方程组问题,退化得到内含任意形状平行三维多裂纹问题的超奇异积分方程组;推导出平行三维多裂纹问题的裂纹前沿广义奇异应力场解析表达式、定义了广义(应力、应变能)强度因子和广义能量释放率;应用有限部积分概念及体积力法,为超奇异积分方程组建立了数值求解方法,编制了FORTRAN程序,以平行双裂纹为例,通过典型算例,研究了广义(应力、应变能)强度因子随裂纹位置、裂纹形状及材料参数变化规律,得到裂纹断裂评定准则. 最后,分析了裂纹间干扰、屏蔽作用及其在工程实际中的应用.      

A ferroelectric and ferroelastic 3D hexahedral curvilinear finite element

An isoparametric 3D electromechanical hexahedral finite element integrating a 3D phenomenological ferroelectric and ferroelastic constitutive law for domain switching effects is proposed. The model presents two internal variables which are the ferroelectric polarization (related to the electric field) and the ferroelastic strain (related to the mechanical stress). An implicit integration technique of the constitutive equations based on the return-mapping algorithm is developed. The mechanical strain tensor and the electric field vector are expressed in a curvilinear coordinate system in order to handle the transverse isotropy behavior of ferroelectric ceramics. The hexahedral finite element is implemented into the commercial finite element code Abaqus® via the subroutine user element. Some linear (piezoelectric) and non linear (ferroelectric and ferroelastic) benchmarks are considered as validation tests.