Saint-Venant End Effects in Anti-plane Shear for Classes of Linear Piezoelectric Materials

This paper is concerned with further investigation of the effect of mechanical/electrical coupling on the decay of Saint-Venant end effects in linear piezoelectricity. Saint-Venant's principle and related results for elasticity theory have received considerable attention in the literature but relatively little is known about analogous issues in piezoelectricity. The current rapidly developing smart structures technology provides motivation for the investigation of such problems. The decay of Saint-Venant end effects is investigated in the context of anti-plane shear deformations for linear homogeneous piezoelectric solids. For a rather general class of anisotropic piezoelectric materials, the governing partial differential equations of equilibrium are a coupled system of second-order partial differential equations for the mechanical displacement u and electric potential ?. The traction boundary-value problem with prescribed surface charge can be formulated as an oblique derivative boundary-value problem for this elliptic system. Energy-decay estimates using differential inequality methods are used to study the axial decay of solutions on a semi-infinite strip subjected to non-zero boundary conditions only at the near end. This analysis is carried out for a rather general class of materials (the tetragonal ${\bar 4}$ crystal class). The boundary-value problem involves a full coupling of mechanical and electrical effects. There are four independent material constants appearing in the problem. An explicit estimated decay rate (a lower bound for the actual decay rate) is obtained in terms of two dimensionless piezoelectric parameters d ₀,r, the first of which provides a measure of the degree of piezoelectric coupling. The estimated decay rate is shown to be monotone decreasing with increasing values of the coupling parameter d ₀. In the limit as d ₀→0, we recover the exact decay rate for the purely mechanical case. Thus, for the tetragonal ${\bar 4}$ class of materials, piezoelectric end effects are predicted to penetrate further into the strip than their elastic counterparts, confirming recent results obtained in other contexts in linear piezoelectricity.     

A theoretical analysis of piezoelectric/composite laminate with larger-amplitude deflection effect,Part II: Hermite differential quadrature method and application

Incorporating the effects of larger-amplitude deflection and electro-elastical properties of piezoelectric lamina, the Hamilton’s variation principle was used to deduce the fundamental formulations of smart anisotropic composite plate in Part I in terms of Reddy’s simple higher-order theory. In order to solve the five highly coupled nonlinear partial differential equations with complicated overlapping boundary conditions, a novel numerical method-Hermite differential quadrature (HDQ) method was developed to implement the differential equations with complicated overlapping boundary conditions. Based on the presently developed HDQ method, any orders derivatives of the unknown functions or any boundary conditions can be point-collocation-based discretized by a set of point-values along x- and y-direction. Then, a system of complete algebraic nonlinear equations can be constructed to calculate out the final point-values of the mid-plane displacements by using the governing equations and relative boundary conditions with HDQ method. Finally, some detailed numerical examples for the anisotropic piezoelectric/composite laminate with the distributed poling directions of piezoelectric layer and fiber orientations of composite layers were studied to validate the developed theoretical analysis model and HDQ numerical method.     

Analysis of 2D infinite anisotropic elastic strips and semi-strips

Summary  Using Stroh's formalism and the theory of analytic functions, simple and explicit solutions for a line dislocation in an infinite anisotropic elastic strip are obtained. The two boundaries of the strip are free of traction. The problem of a dislocation in an anisotropic elastic semi-infinite strip with traction-free boundaries is also studied. A set of singular integral equations governing the unknown functions is derived. When the medium is orthogonal anisotropic and the coordinate axes x ₁ x ₂ x ₃ are coincident with the material principal axes, all the eigenvalues of the material coefficient matrix are pure imaginary. Explicit expressions of the unknown functions are given for this case. The results obtained are valid not only for plane and anti-plane problems but also for coupled problems between in-plane and out-of-plane deformations. Received 30 October 2000; accepted for publication 28 March 2001     

On the asymptotic behavior of solutions of linear second-order boundary-value problems on a semi-infinite strip

This paper treats the asymptotic behavior of solutions of a linear secondorder elliptic partial differential equation defined on a two-dimensional semiinfinite strip. The equation has divergence form and variable coefficients. Such equations arise in the theory of steady-state heat conduction for inhomogeneous anisotropic materials, as well as in the theory of anti-plane shear deformations for a linearized inhomogeneous anisotropic elastic solid. Solutions of such equations that vanish on the long sides of the strip are shown to satisfy a theorem of Phragmén-Lindelöf type, providing estimates for the rate of growth or decay which are optimal for the case of constant coefficients. The results are illustrated by several examples. The estimates obtained in this paper can be used to assess the influence of inhomogeneity and anisotropy on the decay of end effects arising in connection with Saint-Venant's principle.     

Saint-Venant end effects for plane deformation of sandwich strips

The purpose of this paper is to draw attention to the fact that the routine application of Saint-Venant's principle in the solution of elasticity problems for sandwich type structures is not justified in general. This is illustrated in the context of the plane problem of elasticity for a sandwich strip composed of two dissimilar isotropic materials. The exponential decay of end effects is characterized in terms of a complex eigenvalue. For the case of a sandwich with relatively soft middle core, the characteristic decay length is shown to be much greater than that for an homogeneous isotropic strip. The results are analogous to those obtained previously by the authors for highly anisotropic and composite materials.     

End effects for anti-plane shear deformations of sandwich structures

The purpose of this research is to further investigate the effects of material inhomogeneity and the combined effects of material inhomogeneity and anisotropy on the decay of Saint-Venant end effects. Saint-Venant decay rates for self-equilibrated edge loads in symmetric sandwich structures are examined in the context of anti-plane shear for linear anisotropic elasticity. The problem is governed by a second-order, linear, elliptic, partial differential equation with discontinuous coefficients. The most general anisotropy consistent with a state of anti-plane shear is considered, as well as a variety of boundary conditions. Anti-plane or longitudinal shear deformations are one of the simplest classes of deformations in solid mechanics. The resulting deformations are completely characterized by a single out-of-plane displacement which depends only on the in-plane coordinates. They can be thought of as complementary deformations to those of plane elasticity. While these deformations have received little attention compared with the plane problems of linear elasticity, they have recently been investigated for anisotropic and inhomogeneous linear elasticity. In the context of linear elasticity, Saint-Venant's principle is used to show that self-equilibrated loads generate local stress effects that quickly decay away from the loaded end of a structure. For homogeneous isotropic linear elastic materials this is well-documented. Self-equilibrated loads are a class of load distributions that are statically equivalent to zero, i.e., have zero resultant force and moment. When Saint-Venant's principle is valid, pointwise boundary conditions can be replaced by more tractable resultant conditions. It is shown in the present study that material inhomogeneity significantly affects the practical application of Saint-Venant's principle to sandwich structures.     

DECAY RATE OF SAINT-VENANT END EFFECTS FOR PLANE DEFORMATIONS OF PIEZOELECTRIC-PIEZOMAGNETIC SANDWICH STRUCTURES

This paper is concerned with the decay of Saint-Venant end effects for plane deformations of piezoelectric （PE）-piezomagnetic （PM） sandwich structures, where a PM layer is located between two PE layers with the same material properties or reversely. The end of the sandwich structure is subjected to a set of self-equilibrated magneto-electro-elastic loads. The upper and lower surfaces of the sandwich structure axe mechanically free, electrically open or shorted as well as magnetically open or shorted. Firstly the constitutive equations of PE mate- rials and PM materials for plane strain are given and normalized. Secondly, the simplified state space approach is employed to arrange the constitutive equations into differential equations in a matrix form. Finally, by using the transfer matrix method, the characteristic equations for eigen- values or decay rates axe derived. Based on the obtained characteristic equations, the decay rates for the PE-PM-PE and PM-PE-PM sandwich structures are calculated. The influences of the electromagnetic boundary conditions, material properties of PE layers and volume fraction on the decay rates are discussed in detail.     

THE EIGENVALUE PROBLEM AND SAINT-VENANT DECAY RATE FOR A NONHOMOGENEOUS SEMI-INFINITE STRIP

The eigenvalue problem about a nonhomogeneous semi-infinite strip is investigated using the methodology proposed by Papkovich and Fadle for homogeneous plane problems. Two types of nonhomogeneity are considered： （i） the elastic modulus varying with the thickness coor- dinate x exponentially, （ii） it varying with the length coordinate y exponentially. The eigenvalues for the two cases are obtained numerically in plane strain and plane stress states, respectively. By considering the smallest positive eigenvalue, tile Saint-Venant Decay rates are estimated, which indicates material nonhomogeneity has a signifcant influence on the Saint-Venant end effect.     

经典梁热过屈曲问题的解析解

基于非线性经典梁理论,建立了控制轴向和横向变形的基本方程,将两个非线性方程化简为一个关于横向挠度的四阶非线性积分-微分方程。对于本文所考虑的三类边界条件,该方程与相应的边界条件构成了微分特征值问题;直接求解该问题,得到热过屈曲构形的解析解,该解是外加热载荷的函数。为考察热载荷以及边界条件的影响,根据得到的解析解给出了一些数值算例,讨论了梁过屈曲行为的性质。本文得到的解析解可用于验证或改进各类近似理论和数值方法。     

Saint-Venant Decay Rates for an Isotropic Inhomogeneous Linearly Elastic Solid in Anti-Plane Shear

The purpose of this research is to investigate the effects of material inhomogeneity on the decay of Saint-Venant end effects in linear isotropic elasticity. This question is addressed within the context of anti-plane shear deformations of an inhomogeneous isotropic elastic solid. The mathematical issues involve the effects of spatial inhomogeneity on the decay rates of solutions to Dirichlet or Neumann boundary-value problems for a second-order linear elliptic partial differential equation with variable coefficients on a semi-infinite strip. The elastic coefficients are assumed to be smooth functions of the transverse coordinate. The estimated rate of exponential decay with distance from the loaded end (a lower bound for the exact rate of decay) is characterized in terms of the smallest positive eigenvalue of a Sturm–Liouville problem with variable coefficients. Analytic lower bounds for this eigenvalue are used to obtain the desired estimated decay rates. Numerical techniques are also employed to assess the accuracy of the analytic results. A related eigenvalue optimization question is discussed and its implications for the issue of material tailoring is addressed. The results of this paper are applicable to continuously inhomogeneous materials and, in particular, to functionally graded materials. This revised version was published online in August 2006 with corrections to the Cover Date.     

Self-adjoint differential equations arising in finite elasticity for small superimposed deformations

For a number of problems involving small deformations superimposed upon large non-homogeneous deformations of isotropic incompressible hyperelastic materials the governing fourth order linear ordinary differential equations are shown to be self-adjoint. Moreover it is shown that at least in principle every fourth order linear self-adjoint differential equation can be factorized in terms of a single second order self-adjoint operator. These two results unify the derivation of a number of solutions of these equations which have been derived previously by a variety of ad-hoc procedures. Two problems are posed for the interested reader. Firstly the problem of obtaining these self-adjoint differential equations directly from the underlying variational principle and secondly the problem of obtaining the explicit factorizations for the differential equations given.     

Reflection and refraction of plane quasi-longitudinal waves at an interface of two piezoelectric media under initial stresses

This paper is devoted to study a problem of reflection and refraction of quasi-longitudinal waves under initial stresses at an interface of two anisotropic piezoelectric media with different properties. One of the two media is aluminum nitride, which is considered the down piezoelectric medium and the above medium is chosen as PZT-5H ceramics. The two piezoelectric media welded are assumed to be anisotropic of a type of a transversely isotropic crystals (hexagonal crystal structure, class 6 mm). The equations of motion and constitutive relations for the piezoelectric media have been written. Suitable boundary conditions are used to obtain the reflection and refraction coefficients. For an incidence of quasi-longitudinal plane waves, four independent-type amplitude ratios of elastic displacement components for plane waves, called quasi-longitudinal (qP) and quasi-shear vertical (qSV) waves, are shown to exist. Also, it is observed that there exist four dependent amplitude ratios of electric potential, which are proportional to the previous four types. Finally, it is found that the coefficients of reflection and refraction are functions of angle of incidence, elastic constants, piezoelectric potential parameters and the initial stresses. Numerical computations and the results obtained are depicted graphically. In the end, a particular case has been reduced from the present study. This investigation is considered important because the initial stresses in such practical problems are inevitable and may result in frequency shift, a change in the velocity of surface waves and controlling the selectivity of a filter compensation of the devices.     

Internally Pressurized Radially Polarized Piezoelectric Cylinders

The piezoelectric phenomenon has been exploited in science and engineering for decades. Recent advances in smart structures technology have lead to a resurgence of interest in piezoelectricity, and in particular, in the solution of fundamental boundary-value problems. In this paper, we develop an analytic solution to the axisymmetric problem of an infinitely long, radially polarized, radially orthotropic piezoelectric hollow circular cylinder. The cylinder is subjected to uniform internal pressure, or a constant potential difference between its inner and outer surfaces, or both. An analytic solution to the governing equilibrium equations (a coupled system of second-order ordinary differential equations) is obtained. On application of the boundary conditions, the problem is reduced to solving a system of linear algebraic equations. The stress distributions in the cylinder are obtained numerically for two typical piezoceramics of technological interest, namely PZT-4 and BaTiO₃. It is shown that the hoop stresses in a cylinder composed of these materials can be made virtually uniform throughout the cross-section by applying an appropriate set of boundary conditions. This revised version was published online in June 2006 with corrections to the Cover Date.     

Electrostatic stress analysis of an anisotropic piezoelectric half-plane under surface electromechanical loading

The generalized stress components on an anisotropic piezoelectric half-plane boundary under surface electromechanical loading are investigated. It is found that the behaviors of generalized stress components are related to matrices Γ and Ω, which have the same form as those for the purely elastostatic problem. Matrices Γ and Ω contain all the electro-mechanical coupling phenomena of the generalized stress components. All elements of matrices Γ and Ω are expressed explicitly in terms of generalized elastic stiffness for monoclinic piezoelectric materials with the plane of symmetry at x₃ = 0 and for transversely isotropic piezoelectric materials in which the coupled effects between the mechanical (electrical) deformations induced by electrical (mechanical) loadings are studied analytically. A numerical example of the electro-mechanical coupling behavior for PZT-4 is also given.     

Saint-Venant decay rates for an anisotropic and non-homogeneous mixture of elastic solids in anti-plane shear

The purpose of this research is to investigate the influence of material inhomogeneity and anisotropy on the decay of Saint-Venant end effects in anti-plane shear deformations of linear mixtures of elastic solids. The spatial decay of solutions of a boundary value problem with variable coefficients on a semi-infinite strip is investigated. The results may be interpreted in terms of a Saint-Venant principle for anti-plane shear deformations of linear anisotropic mixtures of elastic solids. As our first results have a very general point of view, we study some examples in detail.     

Decoupled formulation of piezoelectric elasticity under generalized plane deformation and its application to wedge problems

This paper presents the formulation of piezoelectric elasticity under generalized plane deformation derived from the three-dimensional theory. There are four decoupled classes in the generalized plane deformation formulation, i.e. when l₃(μ)=l₂^*(μ)=0, l₃(μ)=l₃^*(μ)=0, l₃^*(μ)=l₂^*(μ)=0 or l₃(μ)=l₃^*(μ)=l₂^*(μ)=0. Only the inplane fields of the first class and the antiplane field of the second class include the piezoelectric effect. Several examples of wedge problem often encountered in smart structures, such as sensors or actuators are studied to examine the stress singularity near the apex of the structure. The bonded materials to the PZT-4 wedge are PZT-5, graphite/epoxy or aluminum (conductor). The influencing factors on the singular behavior of the electro-elastic fields include the wedge angle, material type, poling direction, and the boundary and interface conditions. The numerical results of the first case are compared with Xu's graphs and some comments are made in detail. In addition, some new results regarding the antiplane stress singularity of the second class are obtained via the case study. The coupled singularity solutions under generalized plane deformation are also investigated to seek the conditions of the weakest or vanishing singular stress fields.     

On the one-dimensional modeling of camber bending deformations in active anisotropic slender structures

The paper presents a one-dimensional model for anisotropic active slender structures that captures arbitrary cross-sectional deformations. The 1-D geometrically-nonlinear static problem is derived by an asymptotic reduction process from the equations of 3-D electroelasticity. In addition to the conventional (bending–extension–shear–twist) beam strain measures, it includes a Ritz approximation to account for arbitrary deformation shapes of the finite-size cross-sections. As a particular case, closed-form analytical expressions are derived for the linear static equilibrium of a composite thin strip with surface-mounted piezoelectric actuators. This solution is based on a boundary-layer approximation to the static equilibrium equations in regions where Saint-Venant’s principle for elastic bodies cannot be applied and includes camber bending deformations to account for the local bimoments induced by the distributed actuation in a finite-width strip.     

Mode III fracture problem of an arbitrarily oriented crack in a FGPM strip bonded to a FGPM half plane

This paper studies the mode III electro-elastic field of a cracked functionally graded piezoelectric strip bonded to a functionally graded piezoelectric half plane. The crack is oriented in arbitrary direction. The material properties along x-axis vary in exponential form. By using the Fourier transform, the problem can be formulated into a system of singular integral equations and solved by applying the Gauss–Chebyshev integration formula. The effects come from the edge, crack orientation and the nonhomogeneous material parameters on intensity factors are discussed graphically.     

A piezoelectric material strip with a crack perpendicular to its boundary surfaces

A cracked piezoelectric material strip under combining mechanical and electrical loads is considered. The crack is vertical to the top and bottom edges of the strip. The edges of the strip are parallel to the x-axis and perpendicular to the z-axis. When a piezoelectric ceramic is poled, it exhibits transversely isotropic behavior. Among many possible poled axis orientations, a particular orientation when the poling direction lies parallel to x-axis is examined in this paper. Both impermeable crack and permeable crack assumptions are considered. Numerical results are included for three kinds of fracture mechanics specimens, namely an edge-cracked strip, a double edge-cracked strip, and a center-cracked strip, subjected to uniform tensions and uniform electric displacement loads simultaneously, at the far ends. In addition, an edge-cracked strip under pure bending and uniform electric displacement loads at the far ends is also investigated in this paper.     

Elastodynamic Green's functions for a laminated piezoelectric cylinder

Elastodynamic Green's functions for a piezoelectric structure represent the electro-mechanical response due to a steady-state point source as either a unit force or a unit charge. Herein, Green's functions for a laminated circular piezoelectric cylinder are constructed by means of the superposition of modal data from the spectral decomposition of the operator of the equations governing its dynamic behavior. These governing equations are based on a semi-analytical finite element formulation where the discretization occurs through the cylinder's thickness. Examples of a homogeneous PZT-4 cylinder and a two-layer cylinder composed of a PZT-4 material at crystal orientations of ±30° with the longitudinal axis are presented. Numerical implementation details for these two circular cylinders show the convergence and accuracy of these Green's functions.