The boundary layer phenomena in two-dimensional transversely isotropic piezoelectric media by exact symplectic expansion

• Published in: International journal for numerical methods in engineering 2007-03, Vol.69 (11), p.2381-2408
• Main Authors: Leung, A. Y. T., Xu, Xinsheng, Gu, Qian, Leung, C. T. O., Zheng, J. J.
• Format: Article
• Language: English
• Subjects: anisotropy; Boundary conditions; Boundary layer; boundary layer phenomena; Computational techniques; Decay; elasticity; Exact sciences and technology; Fundamental areas of phenomenology (including applications); Mathematical analysis; Mathematical methods in physics; Mathematical models; Media; Physics; Piezoelectricity; Solid mechanics; Static elasticity (thermoelasticity...); Structural and continuum mechanics; Two dimensional
• ISSN: 0029-5981; 1097-0207
• DOI: 10.1002/nme.1855

A separable variable method is introduced to find the exact homogeneous solutions of a two‐dimensional transversely isotropic piezoelectric media to handle general boundary conditions. The usual method of separable variables for partial differentiation equations cannot be readily applicable due to the tangling of the unknowns and their derivatives. Introducing dual variables of stresses, we obtain a set of first‐order Hamiltonian equations whose eigensolutions are symplectic spanning over the solution space to cover all possible boundary conditions. The solutions consist of two parts. The first part is the derogative zero‐eigenvalue solutions of the Saint Venant type together with all their Jordan chains. The second part is the decaying non‐zero‐eigenvalue solutions describing the boundary layer effects. The classical solutions are actually the zero‐eigenvalue solutions representing the simple extension, bending, equipotential field, and the uniform electric displacement. On the other hand, the non‐zero‐eigenvalue solutions represent the localized solutions, which are sensitive to the boundary conditions and are decaying rapidly with respect to the distance from the boundaries. Some rate‐of‐decay curves of the newly found non‐zero‐eigenvalue solutions are shown by numerical examples. Finally, the complete boundary layer effects are quantified for the first time. Copyright © 2006 John Wiley & Sons, Ltd.