Applications of the Trefftz finite element method to anti-plane electroelastic problems are presented in this paper. A dual variational functional is constructed and used to derive Trefftz finite element formulation. Special trial functions which satisfy boundary conditions are also used to develop a special purpose element with local defects. The performance of the proposed element model is assessed by an example and comparison is made with results obtained by other approaches. The Trefftz finite element approach is demonstrated to be ideally suited for the analysis of the anti-plane problem.
IntroductionDuring the past decades the hybrid-Trefftz finite element (FE) model, originating in 1977 [1], has been considerably improved and has now become a highly efficient computational tool for the solution of complex boundary value problems. Up to now, T-elements have been successfully applied to problems of elasticity [2,3], Kirchhoff plates [4, 5], moderately thick Reissner-Mindlin plates [6-8], thick plates [9], general 3-D solid mechanics [10, 11], axisymmetric solid mechanics [12], potential problems [13, 14], shells [15], elastodynamic problems [16-18], transient heat conduction analysis [19], geometrically nonlinear plates [20-23] and materially nonlinear elasticity [24,25]. Further, the concept of special purpose functions has been found to be of great efficiency in dealing with various geometry or load-dependent singularities and local effects (e.g., obtuse or reentrant corners, cracks, circular or elliptic holes, concentrated or patch loads) [2-4, 26, 27].Planar crack problems have been examined by Sabino et al. [28] using the Trefftz boundary element method and by Freitas and Ji [29] using an equilibrium element model. In the present paper, we confine our attention to the applications of Trefftz FE method to anti-plane electroelastic problems. A pair of dual variational functionals is presented and used to derive the corresponding Trefftz FE formulation of piezoelectric materials. Special functions that satisfy traction-free conditions along crack faces are used as trial functions for those elements containing a crack. Numerical solutions to mode III problems obtained by the Trefftz FE method are compared with those obtained by other approaches.
2Basic formulations for anti-plane problem
Basic equationsIn the case of anti-plane shear deformation involving only out-of-plane displacement u z and in-plane electric fields, we havewhere / is electrical potential. The differential governing equation can be written as c 44 r 2 u z þ e 15 r 2 / ¼ 0; e 15 r 2 u z À j 11 r 2 / ¼ 0 in Xwith the constitutive equations r xz r yz