Exact three-dimensional interface stress and electrode-effect analysis of multilayer piezoelectric transducers under torsion

In this paper, exact three-dimensional analysis for torsion of multilayer piezoelectric transducers is presented to highlight the remarkable effects of electrode and size effects on interface stresses and deformation. A generalized formulation is introduced for torsion of an arbitrarily layered cross section. In order to reach an advanced and better device, multilayer piezoelectric materials with different properties are formulated while axis of torsion and the polarization axes of each layer make different orientations with respect to each other. An exact formulation for n-layer piezoelectric device with rectangular cross section is presented to compensate considerable deviation of previous studies’ results from exact solution due to the assumption of linear distribution for the electric potential. In order to improve the performance of the transducers, being utilized in industry, several case studies are presented in order to investigate the influence of different parameters, e.g. thickness and material properties of electrodes over shear and peel stresses. Besides, a correction factor is introduced to completely compensate the effect of ignoring electrodes.     

Exact static analysis of a rotating piezoelectric spherical shell

This paper gives an exact analysis of a rotating piezoelectric spherical shell with arbitrary thickness. At first, three displacement functions are introduced to simplify the basic equations of a spherically isotropic, piezoelectric medium. By expanding the displacement functions as well as the electric potential in terms of spherical harmonics, the basic equations of equilibrium are converted to an uncoupled Euler type, second order ordinary differential equation and a coupled system of three second order ordinary differential equations. A general solution to the homogeneous equations of equilibrium is then derived. The static analysis of a rotating spherical shell is then performed and numerical example is presented. The project supported by the National Natural Science Foundation of China, the Zhejiang Provincial Natural Science Foundation, and the Japanese Committee of Culture, Education and Science.     

Compensation of torsion in rods by piezoelectric actuation

This paper considers the compensation of torsional deformations in rods with the help of thin integrated piezoelectric actuator layers. A laminated orthotropic rod is considered, for which the material properties of each layer are assumed to be homogenous. For the sake of a generalization, the piezoelectric actuation is expressed in terms of eigenstrains. The main scope is the derivation of a distribution of eigenstrains that is able to completely compensate the angle of twist caused by external torsional moments. Saint Venant’s theory of torsion for laminated orthotropic rods is extended for the presence of eigenstrains, which is performed by introducing an additional warping function. It is shown that the actuating torsional moment is a function of the eigenstrains and the additional warping function. For the example of a rectangular cross section, an analytic solution for the actuating moment and the additional warping function is presented. The results are verified by three-dimensional finite-element computations showing a very good accordance with the theoretical results over a large parameter range.     

Exact electroelastic analysis of functionally graded piezoelectric shells

A paper focuses on implementation of the sampling surfaces (SaS) method for the three-dimensional (3D) exact solutions for functionally graded (FG) piezoelectric laminated shells. According to this method, we introduce inside the nth layer I_n not equally spaced SaS parallel to the middle surface of the shell and choose displacements and electric potentials of these surfaces as basic shell variables. Such choice of unknowns yields, first, a very compact form of governing equations of the FG piezoelectric shell formulation and, second, allows the use of strain–displacement equations, which exactly represent rigid-body motions of the shell in any convected curvilinear coordinate system. It is worth noting that the SaS are located inside each layer at Chebyshev polynomial nodes that leads to a uniform convergence of the SaS method. As a result, the SaS method can be applied efficiently to 3D exact solutions of electroelasticity for FG piezoelectric cross-ply and angle-ply shells with a specified accuracy by using a sufficient number of SaS.     

An efficient coupled layerwise theory for static analysis of piezoelectric sandwich beams

Summary An efficient one-dimensional model is developed for the statics of piezoelectric sandwich beams. Third-order zigzag approximation is used for axial displacement, and the potential is approximated as piecewise linear. The displacement field is expressed in terms of three primary displacement variables and the electric potential variables by satisfying the conditions of zero transverse shear stress at the top and bottom and its continuity at layer interfaces. The deflection field accounts for the piezoelectric transverse normal strain. The governing equations are derived using a variational principle. The present results agree very well with the exact solution for thin and thick highly inhomogeneous simply supported hybrid sandwich beams. The developed theory can accurately model open and closed circuit boundary conditions. The first author is grateful to DST, Government of India, for financial support for this work.     

Asymptotic modeling of linearly piezoelectric slender rods

The piezoelectric thin plate modeling already derived by the authors is extended to rod-like structures. Two models corresponding to sensor or actuator behavior are obtained. The conditions of existence of non-local terms in the limit models are discussed. To cite this article: T. Weller, C. Licht, C. R. Mecanique 336 (2008).     

A semi-analytical solution for static and dynamic analysis of plates with piezoelectric patches

A modified mixed variational principle for piezoelectric materials is established and the state-vector equation of piezoelectric plates is deduced directly from the principle. Then the exact solution of the state-vector equation is simply given, and based on the semi-analytical solution of the state-vector equation, a realistic mathematical model is proposed for static analysis of a hybrid laminate and dynamic analysis of a clamped aluminum plate with piezoelectric patches. Both the plate and patches are considered as two three-dimensional piezoelectric bodies, but the same linear quadrilateral element is used to discretize the plate and patches. This method accounts for the compatibility of generalized displacements and generalized stresses on the interface between the plate and patches, and the transverse shear deformation and the rotary inertia of the plate and patches are also considered in the global algebraic equation system. Meanwhile, there is no restriction on the thickness of plate and patches. The model can be also modified to achieve a semi-analytical solution for the transient responses to dynamic loadings and the vibration control of laminated plate with piezoelectric patches or piezoelectric stiffeners.     

Exact closed-form solutions for the static analysis of multi-cracked gradient-elastic beams in bending

Cracks and other forms of concentrated damage can significantly affect the performance of slender beams under static and dynamic loads. The computational model for such defects often consists of a localised reduction in the flexural stiffness, which is macroscopically equivalent to a beam where the undamaged parts are hinged at the position of the crack, with a rotational spring taking into account the residual stiffness (“discrete spring” model). It has been recently demonstrated that this model is equivalent to an inhomogeneous Euler–Bernoulli beam in which a Dirac’s delta is added to the bending flexibility at the position of each damage (“flexibility crack” model). Since these models concentrate the increased curvature at a single abscissa, a jump discontinuity appears in the field of rotations. This study presents an improved representation of cracked slender beams, based on a general class of gradient elasticity with both stress and strain gradient, which allows smoothing the singularities in the flexibility crack model. Exact closed-form solutions are derived for the static response of slender gradient-elastic beams in flexure with multiple cracks, and the numerical examples demonstrate the effects of the nonlocal mechanical parameters (i.e. length scales of the gradient elasticity) in this context.     

压电层合梁静动力学分析的高精度梁单元

给出了一个对复合材料压电层合梁进行数值分析的高精度压电层合梁单元。基于Shi三阶剪切变形板理论的位移场和Layer-wise理论的电势场,用力-电耦合的变分原理及Hamilton原理推导了压电层合梁单元列式。采用拟协调元方法推导了一个可显式给出单元刚度矩阵的两节点压电层合梁单元,并应用于压电层合梁的力-电耦合弯曲和自由振动分析。计算结果表明,该梁单元给出的梁挠度和固有频率与解析解吻合良好,并优于其它梁单元的计算结果,说明了本文所给压电层合梁单元的可靠性和准确性。研究结果可为力-电耦合作用下压电层合梁的力学分析提供一个简单、精确且高效的压电层合梁单元。     

????????????????????????????

介绍了求解复杂截面闭口薄壁杆件扭转问题的网络理论解法,以普朗特应力函 数解法为基础,通过与电路问题比拟,借鉴了电路理论中的网络理论解法. 本文方法对解决复 杂截面闭口薄壁杆件扭转的工程问题有一定的应用价值. 在教学上,可以起到拓宽学生思路 的作用.     

A novel semi-inverse solution method for elastoplastic torsion of heat treated rods

Torsion rods are a primary component of many power transmission and other mechanical systems. The behavior of these rods under elastoplastic torsion is of major concern for designers. Different methods have so far been proposed which deal with the elastoplastic torsion of rods, most of which assume constant yield stress. This assumption produces rough and inaccurate results when the rods are heat treated, since in the process of heat treatment the form of yield stress distribution across the rod cross section changes. We propose a new method for calculating elastoplastic torsion of rods of simply connected cross section which is based on heat treatment observations. In our method the full plastic stress function is obtained by using the semi-inverse method. Elastoplastic stress function is obtained by generalizing the idea of the membrane analogy and using a piecewise continuous stress function. Since the proposed form of yield stress distribution can not be handled by the current Finite Element packages, we produce a computer package with a 3D graphical interface capable of calculating and displaying the 3D elastoplastic stress function, shear stress contours, and torque-angle of rotation per unit length. We show that our method produces excellent agreement for several known cross sections in comparison to methods which assume constant yield stress.     

Shape sensitivity analysis and the energy momentum tensor for the kinematic and static models of torsion

The aim of this paper is to bridge shape sensitivity analysis and configurational mechanics by means of a widespread use of the shape derivative concept. This technique will be applied as a systematic procedure to obtain the Eshelby’s energy momentum tensor associated to the problem under consideration. In order to highlight special features of this procedure and without loss of generality, we focus our attention in the application of shape sensitivity analysis to the problem of twisted straight bars within the framework of linear elasticity.Kinematic and static variational formulations as well as the direct method of sensitivity analysis are used to perform shape derivatives of both models. Integral expressions of first and second order shape derivatives of the total potential energy and the complementary potential energy with respect to an arbitrary transverse cross-section shape change, are achieved. These integral expressions put in evidence the relationship between shape sensitivity analysis and the first and second order Eshelby’s energy momentum tensors. Also, the null divergence property of these tensors is easily proved by comparing, in each case, the domain and boundary integral shape derivative arrived at. Finally, an example with a known exact solution, corresponding to an elastic bar with elliptical transverse cross-section submitted to twist, is presented in order to illustrate the usefulness of these tensors to compute the corresponding shape derivatives.     

Exact solutions for free vibration of transversely isotropic piezoelectric circular plates

By employing the general solution for coupled three-dimensional equations of a transversely isotropic piezoelectric body, this paper investigates the free vibration of a circular plate made of piezoelectric material. Three-dimensional exact solutions are then obtained under two specified boundary conditions, which can be used for both axisymmetric and non-axisymmetric cases. Numerical results are finally presented. The project supported by the National Natural Science Foundation of China (No. 19872060)     

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.     

内含弹性介质功能梯度压电球壳的径向振动调控

针对内含弹性介质的功能梯度压电球壳的径向振动问题,利用变量替换技术,求得了解析解。该解适用于材料参数沿厚度按幂律变化且密度的梯度指标可与其它物理量的梯度指标不同的功能梯度压电球壳,克服了以往分析中所有材料参数梯度指标均假设为相等的局限。数值结果表明,功能梯度压电球壳的径向振动特性可利用材料参数梯度指标和内部弹性介质的刚度进行有效调控。     

Variation in the length of perfectly elastic rods under torsion

On the basis of elastic constitutive relations that reflect geometrically nonlinear second-order effects, we refine the theory of torsion of rectilinear rods of an arbitrary transverse cross-section. In particular, we obtain a universal formula, independent of the material properties, that determines the longitudinal strain arising as the rod undergoes free torsion. According to this formula, the length of a rod made of an isotropic perfectly elastic material can, in contrast to the traditional concepts, either increase or decrease as the rod undergoes torsion. Moreover, the variation in the length depends only on the geometry of the transverse cross-section.