Material stability analysis of particle methods

Abstract

Material instabilities are precursors to phenomena such as shear bands and fracture. Therefore, numerical methods that are intended for failure simulation need to reproduce the onset of material instabilities with reasonable fidelity. Here the effectiveness of particle discretizations in reproducing of the onset of material instabilities is analyzed in two dimensions. For this purpose, a simplified hyperelastic law and a Blatz-Ko material are used. It is shown that the Eulerian kernels used in smooth particle hydrodynamics severely distort the domain of material stability, so that material instabilities can occur in stress states that should be stable. In particular, for the uniaxial case, material instabilities occur at much lower stresses, which is often called the tensile instability. On the other hand, for Lagrangian kernels, the domain of material stability is reproduced very well. We also show that particle methods without stress points exhibit instabilities due to rank deficiency of the discrete equations.

Meshless methods: A review and computer implementation aspects

“…Libersky et al [71] were the first to employ SPH in solid mechanics (impact). Since the original SPH version suffered from spurious instabilities and inconsistencies [97,9,101], many improvements were incorporated into SPH [12,88,20,21,61,62,35,36,100,105]. While SPH and their corrected versions were based on a strong form, other methods were developed in the 1990s, based on a weak form.…”

A temperature‐related homogenization technique and its implementation in the meshfree particle method for nanoscale simulations

“…However, tensile instability occurs in meshfree particle methods as long as Eulerian kernels are employed even when the stress point integration scheme is utilized. Xiao and Belytschko [23] performed material stability analysis of meshfree particle methods. They pointed out that Eulerian kernels severely distorted the material instability while Lagrangian kernels can exactly reproduce material instability.…”

“…The perturbed first Piola-Kirchhoff stress tensor is (23) where C i jkl is the first tangential stiffness tensor, and it is…”

An extended Galerkin weak form and a point interpolation method with continuous strain field and superconvergence using triangular mesh

“…To solve engineering problems, powerful numerical methods have been developed, such as the Finite Element Methods (FEM) [1][2][3][4], Finite Difference Methods [5,6], Finite Volume Methods (FVM) [7,8], and recently Meshfree Methods [9][10][11][12][13][14][15][16][17][18][19][20][21][22][23][24][25][26][27][28]. These methods and techniques not only provide solution tools for many engineering problems, but also extend our minds in the quest for even more effective and robust computational methods.…”