3D BEM for the general piezoelectric solids

M. Denda, C. Y. Wang

Research output: Contribution to journalArticle

11 Scopus citations


In this paper we present a new scheme of three-dimensional boundary element method for the general anisotropic piezoelectric solids. We use the Radon transform representation of the three-dimensional fundamental solutions of piezoelectricity and integrate them analytically over the triangular boundary element with the linear interpolation. This reduces the computation for the system matrices G and H from the standard singular surface integrations to the simple regular line integrations and enables a drastic reduction of the computation time. The integrand of the line integral consists of the product of a function dependent and another function independent on the location vectors representing the source and observation points. The latter function depends only on the material and element properties and thus calculated only once for each element and saved for a repeated use in the calculation of G and H matrices and in the post-processing. Exploitation of this favorable structure results in the further reduction of the computation time for very large systems. The implementation of the proposed method with numerical examples will be presented.

Original languageEnglish (US)
Pages (from-to)2950-2963
Number of pages14
JournalComputer Methods in Applied Mechanics and Engineering
Issue number37-40
StatePublished - Aug 1 2009


All Science Journal Classification (ASJC) codes

  • Computational Mechanics
  • Mechanics of Materials
  • Mechanical Engineering
  • Physics and Astronomy(all)
  • Computer Science Applications


  • 3D BEM
  • General piezoelectric solid
  • Radon transform

Cite this