Two-scale finite element method for piezoelectric problem in periodicstructure

The prediction of the mechanical and electric properties of piezoelectric fibre composites has become an active research area in recent years. By means of introducing a boundary layer problem, some new kinds of two-scale finite element methods for solutions to the electric potential and the displacement for composite material in periodic struc- ture under the coupled piezoelectricity are derived. The coupled two-scale relation of the electric potential and the displacement is set up, and some finite element...     

Accurate modelling of piezoelectric plates: single-layered plate

Summary  The paper presents an efficient two-dimensional approach to piezoelectric plates in the framework of linear theory of piezoelectricity. The approximation of the through-the-thickness variations accounts for the shear effects and a refinement of the electric potential. Using a variational formalism, electromechanically coupled plate equations are obtained for the generalized stress resultants as well as for the generalized electric inductions. The latter are deduced from the conservative electric charge equation, which plays a crucial role in the present model. Emphasis is placed on the boundary conditions at the plate faces. The model is used to examine some problems for cylindrical bending of a single simply supported plate. Number of situations are examined for a piezoelectric plate subject to (i) an applied electric potential, (ii) a surface density of force, and (iii) a surface density of electric charge. The through-thickness distributions of electromechanical quantities (displacements, stresses, electric potential and displacement) are obtained, and compared with results provided by finite element simulations and by a simplified plate model without shear effects. A good agreement is observed between the results coming from the present plate model and finite element computations, which ascertains the effectiveness of the proposed approach to piezoelectric plates. Received 17 July 2000; accepted for publication 26 September 2000     

Multi-scale FE computation for the structures of composite materials with small periodic configuration under condition of coupled thermoelasticity

In this paper, the multi-scale computational method for a structure of composite materials with a small periodic configuration under the coupled thermoelasticity condition is presented. The two-scale asymptotic (TSA) expression of the displacement and the increment of temperature for composite materials with a small periodic configuration under the condition of thermoelasticity are briefly shown at first, then the multi-scale finite element algorithms based on TSA are discussed. Finally the numerical results evaluated by the multi-scale computational method are shown. It demonstrates that the basic configuration and the increment of temperature strongly influence the local strains and local stresses inside a basic cell. The project supported by the National Natural Science Foundation of China (19932030) and Special Funds for Major State Basic Research Projects     

Multi-scale analysis for the structure of composite materials with small periodic configuration under condition of coupled thermoelasticity

The two-scale asymptotic expression of the solution for the increment of temperature in a structure with a small periodic configuration is presented first, and the two-scale asymptotic expression of the displacement for the structure under the coupled thermoelasticity condition is then derived in this paper. In the asymptotic expressions the two-scale coupled relation between the increment of temperature and displacement is included. The approximate solutions and its error estimations are given. The project supported by the National Natural Science Foundation of China (19932030) and Special Funds for Major State Basic Research Projects     

Finite element analysis of the piezoelectric vibrations of quartz plate resonators with higher-order plate theory

A finite element formulation of the piezoelectric vibrations of quartz resonators based on Mindlin plate theory is derived. The higher-order plate theory is employed for the development of a collection of successively higher-order plate elements which can be effective for a broad frequency range including the fundamental and overtone modes of thickness-shear vibrations. The presence of electrodes is also considered for their mechanical effects.The mechanical displacements and electric potential are combined into a generalized displacement field, and the subsequent derivations are carried out with all the generalized equations. Through the standard finite element procedure, the vibration frequency, the vibration mode shapes and the electric potential distribution are obtained. The frequency spectra are compared with some well-known experimental results with good agreement.Our previous experience with finite element analysis of high-frequency quartz plate vibrations leads us to believe that memory and computing time will always remain as key issues despite the advances in computers. Hence, the use of sparse matrix techniques, efficient eigenvalue solvers, and other reduction procedures are explored.     

Characterization of macroscopic tensile strength of polycrystalline metals with two-scale finite element analysis

The objective of this contribution is to develop an elastic-plastic-damage constitutive model for crystal grain and to incorporate it with two-scale finite element analyses based on mathematical homogenization method, in order to characterize the macroscopic tensile strength of polycrystalline metals. More specifically, the constitutive model for single crystal is obtained by combining hyperelasticity, a rate-independent single crystal plasticity and a continuum damage model. The evolution equations, stress update algorithm and consistent tangent are derived within the framework of standard elastoplasticity at finite strain. By employing two-scale finite element analysis, the ductile behaviour of polycrystalline metals and corresponding tensile strength are evaluated. The importance of finite element formulation is examined by comparing performance of several finite elements and their convergence behaviour is assessed with mesh refinement. Finally, the grain size effect on yield and tensile strength is analysed in order to illustrate the versatility of the proposed two-scale model.     

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.     

压电材料中切口/接头端部平面电弹性场奇异性有限元分析

为了对平面载荷作用下压电材料中切口或接头端部附近电弹性场奇异性问题进行分析，首先以应力平衡方程、Maxwell方程和和边界条件为基础，得到一种求解压电材料特征问题的弱式方程；其次，假定楔形切口或接头端部附近单元内位移和电势沿径向分布为指数形式，而周向方向分布则采用泡函数插值，将其代入弱式方程，建立一种只需对楔形切口或接头端部附近周边进行离散的一维简单有限元方法．压电材料的极化轴可以是任意方向．利用该有限元模型讨论了楔形切口角度、极化轴方向和边界条件对奇性场的影响．通过和其它特定情况下的现有解相比，证实了该文有限元数值方法的有效性，而且精度很高．     

Macroscopical non-linear material model for ferroelectric materials inside a hybrid finite element formulation

A new approach for modeling hysteretic non-linear ferroelectric ceramics is presented, based on a fully ferroelectric/ferroelastic coupled macroscopic material model. The material behavior is described by a set of yield functions and the history dependence is stored in internal state variables representing the remanent polarization and the remanent strain. For the solution of the electromechanical coupled boundary value problem, a hybrid finite element formulation is used. Inside this formulation the electric displacement is available as nodal quantity (i.e. degree of freedom) which is used instead of the electric field to determine the evolution of remanent polarization. This involves naturally the electromechanical coupling. A highly efficient integration technique of the constitutive equations, defining a system of ordinary differential equations, is obtained by a customized return mapping algorithm. Due to some simplifications of the algorithm, an analytical solution can be calculated. The automatic differentiation technique is used to obtain the consistent tangent operator. Altogether this has been implemented into the finite element code FEAP via a user element. Extensive verification tests are performed in this work to evaluate the behavior of the material model under pure electrical and mechanical as well as coupled and multi-axial loading conditions.     

压电复合材料层板柱面弯曲的精化理论

本研究旨在建立精确的压电复合材料层板理论。位移场和电势场采用近似表达，其沿板厚的分布通过构造高精度的位移分布函数和电势分布函数来描述。这两个函数由三雏弹性平衡方程和静电平衡方程的特解来导出，从而满足复杂的力电耦合关系和各类连续条件，保证了本文理论的高精度。本文理论仅涉及4个位移和电势变量，且不随层数的增加而增多，较之变量随层数而增多的分层理论简单得多，平衡方程形式简单；也便于发展成有限元等数值模型。通过与三维精确解比较，算例显示了本文理论的高精度和有效性。     

基于高阶位移模式的压电层合板动力有限元模型

建立了含压电片层合板的有限元动力学模型。以位于压电层上下表面处的电场强度和层间电压为未知量,给出了三次函数的电势分布模式,采用Reddy的高阶剪切理论描述板的位移场,假设板厚度方向的正应力为零给出了减缩的本构方程,采用有限元方法,基于Hamilton原理导出结构的动力学方程,然后用静态缩聚的方法压缩掉电场自由度和次要的位移自由度。最后用四边形矩形单元求解了一对称铺层和非对称铺层悬臂板的固有频率,并与ANSYS结果对比验证了本文模型的精确性。     

智能梁的静态形状控制

直接从位移函数出发，用Hamilton原理推导了压电智能梁结构的有限元方程，并采用灵敏度分析方法来确定结构的优化控制电压，在此基础上建立了结构自适应梁。最后采用数值分析来验证其有效性，从仿真结果可以看出文中的方法对梁结构的静态形状控制，有较好的控制精度。     

Topology optimization of piezoelectric nanostructures

We present an extended finite element formulation for piezoelectric nanobeams and nanoplates that is coupled with topology optimization to study the energy harvesting potential of piezoelectric nanostructures. The finite element model for the nanoplates is based on the Kirchoff plate model, with a linear through the thickness distribution of electric potential. Based on the topology optimization, the largest enhancements in energy harvesting are found for closed circuit boundary conditions, though significant gains are also found for open circuit boundary conditions. Most interestingly, our results demonstrate the competition between surface elasticity, which reduces the energy conversion efficiency, and surface piezoelectricity, which enhances the energy conversion efficiency, in governing the energy harvesting potential of piezoelectric nanostructures.     

Static analysis of simply supported functionally graded and layered magneto-electro-elastic plates

In this article, static analysis of functionally graded, anisotropic and linear magneto-electro-elastic plates have been carried out by semi-analytical finite element method. A series solution is assumed in the plane of the plate and finite element procedure is adopted across the thickness of the plate such a way that the three-dimensional character of the solution is preserved. The finite element model is derived based on constitutive equation of piezomagnetic material accounting for coupling between elasticity, electric and magnetic effect. The present finite element is modeled with displacement components, electric potential and magnetic potential as nodal degree of freedom. The other fields are calculated by post-computation through constitutive equation. The functionally graded material is assumed to be exponential in the thickness direction. The numerical results obtained by the present model are in good agreement with available functionally graded three-dimensional exact benchmark solutions given by Pan and Han [Pan, E., Han, F., in press. Green’s function for transversely isotropic piezoelectric functionally graded multilayered half spaces. Int. J. Solids Struct.]. Numerical study includes the influence of the different exponential factor, magneto-electro-elastic properties and effect of mechanical and electric type of loading on induced magneto-electro-elastic fields. In addition further study has been carried out on non-homogeneous transversely isotropic FGM magneto-electro-elastic plate available in the literature [Chen, W.Q., Lee, K.Y., Ding, H.J., 2005. On free vibration of non-homogeneous transversely isotropic magneto-electro-elastic plates].     

A 3D multi-field element for simulating the electromechanical coupling behavior of dielectric elastomers

We propose a multi-field implicit finite element method for analyzing the electromechanical behavior of dielectric elastomers. This method is based on a four-field variational principle, which includes displacement and electric potential for the electromechanical coupling analysis, and additional independent fields to address the incompressible constraint of the hyperelastic material. Linearization of the variational form and finite element discretization are adopted for the numerical implementation. A general FEM program framework is developed using C ++ based on the open-source finite element library deal.II to implement this proposed algorithm. Numerical examples demonstrate the accuracy, convergence properties, mesh-independence properties, and scalability of this method. We also use the method for eigenvalue analysis of a dielectric elastomer actuator subject to electromechanical loadings. Our finite element implementation is available as an online supplementary material.     

2D dynamic analysis of cracks and interface cracks in piezoelectric composites using the SBFEM

The dynamic stress and electric displacement intensity factors of impermeable cracks in homogeneous piezoelectric materials and interface cracks in piezoelectric bimaterials are evaluated by extending the scaled boundary finite element method (SBFEM). In this method, a piezoelectric plate is divided into polygons. Each polygon is treated as a scaled boundary finite element subdomain. Only the boundaries of the subdomains need to be discretized with line elements. The dynamic properties of a subdomain are represented by the high order stiffness and mass matrices obtained from a continued fraction solution, which is able to represent the high frequency response with only 3–4 terms per wavelength. The semi-analytical solutions model singular stress and electric displacement fields in the vicinity of crack tips accurately and efficiently. The dynamic stress and electric displacement intensity factors are evaluated directly from the scaled boundary finite element solutions. No asymptotic solution, local mesh refinement or other special treatments around a crack tip are required. Numerical examples are presented to verify the proposed technique with the analytical solutions and the results from the literature. The present results highlight the accuracy, simplicity and efficiency of the proposed technique.     

Finite element modeling of smart piezo structure: considering dependence of piezoelectric coefficients on electric field

In this work, variation of piezoelectric strain coe?cient and permittivity with change in electric field is included in constitutive equations of piezoelectricity and used for developing finite element model of a plate instrumented with piezoelectric patches. Simulation results show that nodal displacement response and sensor voltage is less controlled if variation of piezoelectric strain coe?cient and permittivity with change in electric field is not included in finite element modelling as compared to the case in which these variations are included in finite element modeling. An experiment is performed so as to validate simulation results.     

均布荷载作用下悬臂磁电弹性梁的解析解

对磁电弹性平面问题进行了研究，给出了用拟调和位移函数表达的通解，进而以试凑法按平面应力问题推导出了均布荷载作用下悬臂磁电弹性粱的解析解，所得解有易于理解、便于校对、形式统一简洁的特点。本文还将计算结果与压电材料和弹性材料相应结果进行了分析、比较，为验证各种数值计算方法提供了参考依据。     

PZT-4紧凑拉伸试样的断裂分析

基于线性压电材料的复势理论,通过解析分析,导出了一种分析有限压电板裂纹问题的解析数值方法. 首先,计算了含中心裂纹有限板的断裂参数,与Woo和Wang的解析数值法(Int J Fract, 1993, 62: 203$\sim$218)相比较,表明该方法具有很高的精度和很好的计算效率. 随后,采用该方法和有限元法计算了PZT-4紧凑拉伸试样在绝缘裂纹面边界条件下断裂时的断裂参数,发现各断裂参数的临界值分散性很大,不能作为压电材料的单参数断裂准则. 进而,针对试样真实的裂隙形状,采用有限元法计算了裂隙尖端的应力、电位移场,比较了裂隙内介质的介电性能对裂隙尖端场的影响,计算了带微裂纹的真实裂隙模型的断裂参数并进行了理论分析.      

直接有限元法求解广义磁热弹二维旋转问题

为了验证直接有限元法求解广义磁热弹耦合旋转问题的有效性及准确性,该文基于Lord和Shulman(L-S)广义热弹性理论,采用直接有限元方法,求解了置于磁场中的旋转半无限大体受热冲击作用的动态响应问题.文中给出了L-S型广义磁热弹耦合旋转问题的控制方程,建立了L-S型广义磁热弹旋转问题的虚位移原理,推导得到了相应的有限...