### Abstract

A new method for numerical simulation of failure behavior, namely, FEM-β, is proposed. For a continuum model of a deformable body, FEM-β solves a boundary value problem by applying particle discretization to a displacement field; the domain is decomposed into a set of Voronoi blocks and the non-overlapping characteristic functions for the Voronoi blocks are used to discretize the displacement function. By computing average strain and average strain energy, FEM-β obtains a numerical solution of the variational problem that is transformed from the boundary value problem. In a rigorous form, FEM-β is formulated for a variational problem of displacement and stress with different particle discretization, i.e., the non-overlapping characteristic function of the Voronoi blocks and the conjugate Delaunay tessellations, respectively, are used to discretize the displacement and stress functions. While a displacement field is discretized with non-smooth functions, it is shown that a solution of FEM-β has the same accuracy as that of ordinary FEM with triangular elements. The key point of FEM-β is the ease of expressing failure as separation of two adjacent Voronoi blocks owing to the particle discretization that uses non-overlapping characteristic functions. This paper explains these features of FEM-β with results of numerical simulation of several example problems.

Original language | English |
---|---|

Pages (from-to) | 681-703 |

Number of pages | 23 |

Journal | Journal of the Mechanics and Physics of Solids |

Volume | 53 |

Issue number | 3 |

DOIs | |

Publication status | Published - 2005 Mar |

Externally published | Yes |

### Fingerprint

### Keywords

- Discontinuous finite element method
- Failure analysis
- Non-overlapping shape functions
- Particle modeling
- Voronoi and Delaunay tessellation

### ASJC Scopus subject areas

- Mechanical Engineering
- Mechanics of Materials
- Condensed Matter Physics

### Cite this

*Journal of the Mechanics and Physics of Solids*,

*53*(3), 681-703. https://doi.org/10.1016/j.jmps.2004.08.005

**Proposal of FEM implemented with particle discretization for analysis of failure phenomena.** / Hori, Muneo; Oguni, Kenji; Sakaguchi, Hide.

Research output: Contribution to journal › Article

*Journal of the Mechanics and Physics of Solids*, vol. 53, no. 3, pp. 681-703. https://doi.org/10.1016/j.jmps.2004.08.005

}

TY - JOUR

T1 - Proposal of FEM implemented with particle discretization for analysis of failure phenomena

AU - Hori, Muneo

AU - Oguni, Kenji

AU - Sakaguchi, Hide

PY - 2005/3

Y1 - 2005/3

N2 - A new method for numerical simulation of failure behavior, namely, FEM-β, is proposed. For a continuum model of a deformable body, FEM-β solves a boundary value problem by applying particle discretization to a displacement field; the domain is decomposed into a set of Voronoi blocks and the non-overlapping characteristic functions for the Voronoi blocks are used to discretize the displacement function. By computing average strain and average strain energy, FEM-β obtains a numerical solution of the variational problem that is transformed from the boundary value problem. In a rigorous form, FEM-β is formulated for a variational problem of displacement and stress with different particle discretization, i.e., the non-overlapping characteristic function of the Voronoi blocks and the conjugate Delaunay tessellations, respectively, are used to discretize the displacement and stress functions. While a displacement field is discretized with non-smooth functions, it is shown that a solution of FEM-β has the same accuracy as that of ordinary FEM with triangular elements. The key point of FEM-β is the ease of expressing failure as separation of two adjacent Voronoi blocks owing to the particle discretization that uses non-overlapping characteristic functions. This paper explains these features of FEM-β with results of numerical simulation of several example problems.

AB - A new method for numerical simulation of failure behavior, namely, FEM-β, is proposed. For a continuum model of a deformable body, FEM-β solves a boundary value problem by applying particle discretization to a displacement field; the domain is decomposed into a set of Voronoi blocks and the non-overlapping characteristic functions for the Voronoi blocks are used to discretize the displacement function. By computing average strain and average strain energy, FEM-β obtains a numerical solution of the variational problem that is transformed from the boundary value problem. In a rigorous form, FEM-β is formulated for a variational problem of displacement and stress with different particle discretization, i.e., the non-overlapping characteristic function of the Voronoi blocks and the conjugate Delaunay tessellations, respectively, are used to discretize the displacement and stress functions. While a displacement field is discretized with non-smooth functions, it is shown that a solution of FEM-β has the same accuracy as that of ordinary FEM with triangular elements. The key point of FEM-β is the ease of expressing failure as separation of two adjacent Voronoi blocks owing to the particle discretization that uses non-overlapping characteristic functions. This paper explains these features of FEM-β with results of numerical simulation of several example problems.

KW - Discontinuous finite element method

KW - Failure analysis

KW - Non-overlapping shape functions

KW - Particle modeling

KW - Voronoi and Delaunay tessellation

UR - http://www.scopus.com/inward/record.url?scp=12344283945&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=12344283945&partnerID=8YFLogxK

U2 - 10.1016/j.jmps.2004.08.005

DO - 10.1016/j.jmps.2004.08.005

M3 - Article

AN - SCOPUS:12344283945

VL - 53

SP - 681

EP - 703

JO - Journal of the Mechanics and Physics of Solids

JF - Journal of the Mechanics and Physics of Solids

SN - 0022-5096

IS - 3

ER -