A fast object-oriented Matlab implementation of the Reproducing Kernel Particle Method

Ettore Barbieri, Michele Meo

Research output: Contribution to journalArticle

23 Citations (Scopus)

Abstract

Novel numerical methods, known as Meshless Methods or Meshfree Methods and, in a wider perspective, Partition of Unity Methods, promise to overcome most of disadvantages of the traditional finite element techniques. The absence of a mesh makes meshfree methods very attractive for those problems involving large deformations, moving boundaries and crack propagation. However, meshfree methods still have significant limitations that prevent their acceptance among researchers and engineers, namely the computational costs. This paper presents an in-depth analysis of computational techniques to speed-up the computation of the shape functions in the Reproducing Kernel Particle Method and Moving Least Squares, with particular focus on their bottlenecks, like the neighbour search, the inversion of the moment matrix and the assembly of the stiffness matrix. The paper presents numerous computational solutions aimed at a considerable reduction of the computational times: the use of kd-trees for the neighbour search, sparse indexing of the nodes-points connectivity and, most importantly, the explicit and vectorized inversion of the moment matrix without using loops and numerical routines. 2011 Springer-Verlag.
Original languageEnglish
Pages (from-to)581-602
JournalComputational Mechanics
Volume49
Issue number5
DOIs
Publication statusPublished - 2012

Fingerprint Dive into the research topics of 'A fast object-oriented Matlab implementation of the Reproducing Kernel Particle Method'. Together they form a unique fingerprint.

  • Cite this