Solving numerical difficulties for element-free Galerkin analyses

S. H. Ju, H. H. Hsu

Research output: Contribution to journalArticlepeer-review

6 Citations (Scopus)


This study solves the numerical problems associated with the element-free Galerkin method (EFGM) to perform analyses efficiently in shared-memory computers. The truncation error is generally large for the moving least-squares approximation, and this can be overcome by using orthogonal basis functions, 16-byte floats, or the local origin. Then, the analysis accuracy is similar to that obtained with the reproducing kernel particle approximation. Determining the index array of the global stiffness matrix requires a large amount of computer memory. We thus propose a scheme to overcome this problem using slightly more computer time but much less computer memory. A binary search is also proposed to find the support domain nodes for Gaussian points, and this method is much more efficient than the linear search one. A Fortran module is developed to establish parallel solutions in the EFGM, and the programmer does not need to handle the global stiffness directly.

Original languageEnglish
Pages (from-to)273-281
Number of pages9
JournalComputational Mechanics
Issue number2
Publication statusPublished - 2014 Feb

All Science Journal Classification (ASJC) codes

  • Computational Mechanics
  • Ocean Engineering
  • Mechanical Engineering
  • Computational Theory and Mathematics
  • Computational Mathematics
  • Applied Mathematics


Dive into the research topics of 'Solving numerical difficulties for element-free Galerkin analyses'. Together they form a unique fingerprint.

Cite this