Development Of Higher Order Particle Discretization Scheme For Analysis Of Failure Phenomena

This paper presents the higher order extension of Particle Discretization Scheme (PDS) and its implementation in FEM framework (PDS-FEM) to solve boundary value problems of linear elastic solids, including brittle cracks. Higher order PDS defines an approximation f(d)(x) of a function f(x), defined over domain Omega, as the union of local polynomial approximation of f(x) over each Voronoi tessellation elements of Omega. The support of the local polynomial bases being confined to the domain of each Voronoi element, f(d)(x) consists of discontinuities along each Voronoi boundaries. Considering local polynomial approximations over elements of Delaunay tessellation, PDS define bounded derivatives for this discontinuous f(d)(x). Utilizing the inherent discontinuities in f(d)(x), PDS-FEM proposes a numerically efficient treatment for modeling cracks. This novel use of local polynomial approximations in FEM is verified with a set of linear elastic problems, including mode-I crack tip stress field.