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.     

An efficient higher order zigzag theory for composite and sandwich beams subjected to thermal loading

A new efficient higher order zigzag theory is presented for thermal stress analysis of laminated beams under thermal loads, with modification of the third order zigzag model by inclusion of the explicit contribution of the thermal expansion coefficient α₃ in the approximation of the transverse displacement w. The thermal field is approximated as piecewise linear across the thickness. The displacement field is expressed in terms of the thermal field and only three primary displacement variables by satisfying exactly the conditions of zero transverse shear stress at the top and the bottom and its continuity at the layer interfaces. The governing equations are derived using the principle of virtual work. Fourier series solutions are obtained for simply-supported beams. Comparison with the exact thermo-elasticity solution for thermal stress analysis under two kinds of thermal loads establishes that the present zigzag theory is generally very accurate and superior to the existing zigzag theory for composite and sandwich beams.     

A new simple third-order shear deformation theory of plates

An improved simple third-order shear deformation theory for the analysis of shear flexible plates is presented in this paper. This new plate theory is composed of three parts: the simple third-order kinematics of displacements reduced from the higher-order displacement field derived previously by the author; a system of 10th-order differential equilibrium equations in terms of the three generalized displacements of bending plates; five boundary conditions at each edge of plate boundaries. Although the resulting displacement field is the same as that proposed by Murthy, the variational consistent governing equations and the associated proper boundary conditions are derived and identified in this work for the first time in the literature. The applications and accuracy of the present shear deformation theory of plates are demonstrated by analytically solving the differential governing equations of a twisting plate, a bending beam and two bending plates to which the 3-D elasticity solutions are available, and excellent agreements are achieved even for the torsion of a plate with square cross-section as well the local effects of stresses at plate boundaries can be characterized accurately. These analytical solutions clearly show that the simple third-order shear deformation theory developed in this work indeed gives better results than the first-order shear deformation theories and other simple higher-order shear deformation theories, since the present third-order shear flexible theory is based on a more rigorous kinematics of displacements and consists of not only a system of variational consistent differential equations, but also a group of consistent boundary conditions associated with the differential equations. The present simple third-order shear deformation theory can easily be applied to the static and dynamic finite element analysis of laminated plates just like the applications of other popular shear flexible plate theories, and improved results could be obtained from the present simple third-order shear deformable theories of plates.     

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.     

带剪切型压电激励器的板的高阶剪切变形理论

利用压电材料固有的正，逆压电效应可以对结构变形和振动进行控制。与外加电场与极化方向平行于板厚度的压电材料的拉伸作动机制相比，外加电场与极化方向垂直的压电材料的剪切作动机制可以在作动器内产生较小的应力，从而降低作动器边界产生分层破坏的危险。本文对于压电材料的剪切作动机制进行研究，应用三阶剪切变形理论建立带剪切型压电激励器的智能层合板模型。采用哈密顿原理导出带剪切型压电激励器的层合板的控制方程。采用空间法得到了各种边界条件组合条件下板的解析解。数值算例对一三层板采用高阶和一阶剪切变形理论进行计算，结果表明两种理论所得的变形曲线很相似。但对于厚度剪切型激励器而言，由于激励器是引起板的剪切变形，而高阶剪切变形理论比一阶剪切变形理论能更好地反映结构的剪切应变能，因此高阶剪切变形理论可以提供板变形的更为精确的解。因此，对于厚度剪切型激励器，剪切变形理论的选取对于板变形结果的好坏有重要的作用。     

Multiple mode large amplitude vibration of thick orthotropic circular plates

This paper is analytically concerned with the large amplitude vibration of thick orthotropic circular plates incorporating the effects of transverse shear and rotatory inertia. Von Kármán-type field equations written in terms of the three displacement components of the plate are utilized to obtain solutions to clamped stress-free and immovable plates. By means of Galerkin's technique and a numerical Runge-Kutta procedure a multiple-mode analysis is carried out in both cases. Exact solutions are reported for two of the three governing equations. Effects of transverse shear deformation and modal interaction are found to be significant for orthotropic thick plates. The method given here could be extended to the multiple-mode analysis of circular plates with other boundary conditions.     

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     

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

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

Creep buckling and post-buckling analysis of the laminated piezoelectric viscoelastic functionally graded plates

The creep buckling and post-buckling of the laminated piezoelectric viscoelastic functionally graded material (FGM) plates are studied in this research. Considering the transverse shear deformation and geometric nonlinearity, the Von Karman geometric relation of the laminated piezoelectric viscoelastic FGM plates with initial deflection is established. And then nonlinear creep governing equations of the laminated piezoelectric viscoelastic FGM plates subjected to an in-plane compressive load are derived on the basis of the elastic piezoelectric theory and Boltzmann superposition principle. Applying the finite difference method and the Newmark scheme, the whole problem is solved by the iterative method. In numerical examples, the effects of geometric nonlinearity, transverse shear deformation, the applied electric load, the volume fraction and the geometric parameters on the creep buckling and post-buckling of laminated piezoelectric viscoelastic FGM plates with initial deflection are investigated.     

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.     

压电弹性层合板静力机耦合特性的解析解

应用经典权壳理论和压电理论对层合板结构进行简化。对四边简支压电层合板在不同电学边界条件下,包括四边电学开路和电学短路。上下表面受外电压及无外加电压作用,进行分析,求得了电势和挠度的解析表达式,给出了压电层和基体的挠度、电势分布图。     

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.     

Modeling photo-induced deformation of glassy splay-bend and twist nematic sheets

Azobenzene-containing glassy nematic sheets deform in response to light in a complicated way depending on director distribution. To quantify the large-deflected deformation, a theoretical model is developed for the sheets with typical splay-bend and twist director distributions. A third-order in-plane displacement assumption is adopted to characterize the effect of transverse shear deformation, and the necessity is discussed through two examples for which analytical solutions are obtainable. Though this work is an extension of the third-order shear deformable theory for anisotropic laminates, it involves some new ingredients such as varying spontaneous strains and special material symmetries. The results are expected useful for analysis and design of the glassy nematic sheets in actuation applications.     

Coupled refined layerwise theory for dynamic free and forced response of piezoelectric laminated composite and sandwich beams

This work extends a previously presented coupled refined layerwise theory to dynamic analysis of piezoelectric laminated composite and sandwich beams. Contrary to most of the available theories, all the kinematic and stress boundary conditions are satisfied at the interfaces of the piezoelectric layers with the non-zero longitudinal electric field. Moreover, both electrical transverse normal strains and transverse flexibility are taken into account for the first time in the present theory. In the presented formulation a high-order polynomial, an exponential expression and a layerwise term containing the electric field are included in the describing expression of the in-plane displacement of the beam. For the transverse displacement, the coupled refined model uses a combination of continuous piecewise fourth-order polynomials with a layerwise representation of electrical unknowns. The electric field is also approximated as linear across the thickness direction of piezoelectric layers. One of advantages of the present theory is that the mechanical number of the unknown parameters is very small and is independent of the number of the layers. For validation of the proposed model, various free and forced vibration tests for thin and thick laminated/sandwich piezoelectric beams are carried out. For various electrical and mechanical boundary conditions, excellent correlation has been found between the results obtained from the proposed formulation with those resulted from the three-dimensional theory of piezoelasticity.     

Numerical solution of polarization saturation/dielectric breakdown model in 2D finite piezoelectric media

The polarization saturation (PS) model [Gao, H., Barnett, D.M., 1996. An invariance property of local energy release rates in a strip saturation model of piezoelectric fracture. Int. J. Fract. 79, R25–R29; Gao, H., Zhang, T.Y., Tong, P., 1997. Local and global energy release rates for an electrically yielded crack in a piezoelectric ceramic. J. Mech. Phys. Solids 45, 491–510], and the dielectric breakdown (DB) model [Zhang, T.Y., Zhao, M.H., Cao, C.F., 2005. The strip dielectric breakdown model. Int. J. Fract. 132, 311–327] explain very well some experimental observations of fracture of piezoelectric ceramics. In this paper, the nonlinear hybrid extended displacement discontinuity-fundamental solution method (NLHEDD-FSM) is presented for numerical analysis of both the PS and DB models of two-dimensional (2D) finite piezoelectric media under impermeable and semi-permeable electric boundary conditions. In this NLHEDD-FSM, the solution is expressed approximately by a linear combination of fundamental solutions of the governing equations, which includes the extended point force fundamental solutions with sources placed at chosen points outside the domain of the problem under consideration, and the extended Crouch fundamental solutions with extended displacement discontinuities placed on the crack and the electric yielding zone. The coefficients of the fundamental solutions are determined by letting the approximated solution satisfy certain conditions on the boundary of the domain, on the crack face and the electric yielding zone. The zero electric displacement intensity factor in the PS model or the zero electric field strength intensity factor in the DB model at the outer tips of the electric yielding zone is used as a supplementary condition to determine the size of the electric yielding zone. Iteration approaches are adopted in the NLHEDD-FSM. The electric yielding zone is determined, and the extended intensity factors and the local J-integral are calculated for center cracks in piezoelectric strips. The effects of finite domain size, saturation property and different electric boundary conditions, as well as different models on the electric yielding zone and the local J-integral, are studied.     

A new trigonometric shear deformation theory for isotropic,laminated composite and sandwich plates

A new trigonometric shear deformation theory for isotropic and composite laminated and sandwich plates, is developed. The new displacement field depends on a parameter “m”, whose value is determined so as to give results closest to the 3D elasticity bending solutions. The theory accounts for adequate distribution of the transverse shear strains through the plate thickness and tangential stress-free boundary conditions on the plate boundary surface, thus a shear correction factor is not required. Plate governing equations and boundary conditions are derived by employing the principle of virtual work. The Navier-type exact solutions for static bending analysis are presented for sinusoidally and uniformly distributed loads. The accuracy of the present theory is ascertained by comparing it with various available results in the literature. The results show that the present model performs as good as the Reddy’s and Touratier’s shear deformation theories for analyzing the static behavior of isotropic and composite laminated and sandwich plates.     

A hybrid stress approach for laminated composite plates within the First-order Shear Deformation Theory

This paper presents a hybrid stress approach for the analysis of laminated composite plates. The plate mechanical model is based on the so called First-order Shear Deformation Theory, rationally deduced from the parent three-dimensional theory. Within this framework, a new quadrilateral four-node finite element is developed from a hybrid stress formulation involving, as primary variables, compatible displacements and elementwise equilibrated stress resultants. The element is designed to be simple, stable and locking-free. The displacement interpolation is enhanced by linking the transverse displacement to the nodal rotations and a suitable approximation for stress resultants is selected, ruled by the minimum number of parameters. The transverse stresses through the laminate thickness are reconstructed a posteriori by simply using three-dimensional equilibrium. To improve the results, the stress resultants entering the reconstruction process are first recovered using a superconvergent patch-based procedure called Recovery by Compatibility in Patches, that is properly extended here for laminated plates. This preliminary recovery is very efficient from the computational point of view and generally useful either to accurately evaluate the stress resultants or to estimate the discretization error. Indeed, in the present context, it plays also a key role in effectively predicting the shear stress profiles, since it guarantees the global convergence of the whole reconstruction strategy, that does not need any correction to accommodate equilibrium defects. Actually, this strategy can be adopted together with any plate finite element. Numerical testing demonstrates the excellent performance of both the finite element and the reconstruction strategy.     

An efficient and simple refined theory for nonlinear bending analysis of functionally graded sandwich plates

In this paper, an efficient and simple refined theory is presented for nonlinear bending analysis of functionally graded sandwich plates. The theory presented is variationally consistent, does not require the shear correction factor, and gives rise to transverse shear stress variations such that the transverse shear stresses vary parabolically across the plate thickness, satisfying shear-stress-free surface conditions. The energy concept along with the present theory and the first- and third-order shear deformation theories is used to predict the large deflection and the stress distribution across the thickness of functionally graded sandwich plates.     

Axisymmetric bending analysis of thick functionally graded circular plates using fourth-order shear deformation theory

In the present article, axisymmetric bending and stretching of functionally graded (FG) circular plates subjected to uniform transverse loading based on fourth-order shear deformation plate theory (FOST) have been studied. Using a fourth-order shear deformation theory, the solutions for deflection and rotation functions of FG plates are presented in terms of the corresponding quantities for a homogeneous plate using the classical plate theory (CPT), from which solutions one can easily obtain the FOST solutions for axisymmetric bending of FG circular plates. It is assumed that the effective mechanical properties of the functionally graded plates through the thickness are continuous functions of the volume fractions of the constituent parts which are themselves defined by a power-law function. Numerical results for maximum deflection and shear stress are presented for various percentages of ceramic–metal volume fractions. These results are also compared with those obtained from the first-order shear deformation plate theory of Mindlin (FST), the third-order shear deformation plate theory of Reddy (TST) as well as the exact three-dimensional elasticity solution. It is found that although the maximum deflections obtained using FOST and TST are close to each other, the through-thickness shear stress is predicted more accurately by the FOST formulation than by the TST.     

An advanced higher-order theory for laminated composite plates with general lamination angles

This paper proposes a higher-order shear deformation theory to predict the bending response of the laminated composite and sandwich plates with general lamination configurations.The proposed theory a priori satisfies the continuity conditions of transverse shear stresses at interfaces.Moreover,the number of unknown variables is independent of the number of layers.The first derivatives of transverse displacements have been taken out from the inplane displacement fields,so that the C 0 shape functions are only required during its finite element implementation.Due to C 0 continuity requirements,the proposed model can be conveniently extended for implementation in commercial finite element codes.To verify the proposed theory,the fournode C 0 quadrilateral element is employed for the interpolation of all the displacement parameters defined at each nodal point on the composite plate.Numerical results show that following the proposed theory,simple C 0 finite elements could accurately predict the interlaminar stresses of laminated composite and sandwich plates directly from a constitutive equation,which has caused difficulty for the other global higher order theories.