Controlling the onset of numerical fracture in parallelized implementations of the material point method (MPM) with convective particle domain interpolation (CPDI) domain scaling
Autor: | James Guilkey, Rebecca M. Brannon, Michael A. Homel |
---|---|
Rok vydání: | 2015 |
Předmět: |
Numerical Analysis
Applied Mathematics General Engineering Geometry 02 engineering and technology 01 natural sciences Domain (software engineering) 010101 applied mathematics Arbitrarily large 020303 mechanical engineering & transports 0203 mechanical engineering Solid mechanics Fracture (geology) Vector field 0101 mathematics Scaling Algorithm Material point method Mathematics Interpolation |
Zdroj: | International Journal for Numerical Methods in Engineering. 107:31-48 |
ISSN: | 0029-5981 |
DOI: | 10.1002/nme.5151 |
Popis: | The material point method is well suited for large-deformation problems in solid mechanics but requires modification to avoid cell-crossing errors as well as extension instabilities that lead to numerical (nonphysical) fracture. A promising solution is convected particle domain interpolation (CPDI), in which the integration domain used to map data between particles and the background grid deforms with the particle, based on the material deformation gradient. While eliminating the extension instability can be a benefit, it is often desirable to allow material separation to avoid nonphysical stretching. Additionally, large stretches in material points can complicate parallel implementation of CPDI if a single particle domain spans multiple computational patches. A straightforward modification to the CPDI algorithm allows a user-specified scaling of the particle integration domain to control the numerical fracture response, which facilitates parallelization. Combined with particle splitting, the method can accommodate materials with arbitrarily large failure strains. Used with a smeared damage/softening model, this approach will prevent nonphysical numerical fracture in situations where the material should remain intact, but the effect of a single velocity field on localization may still produce errors in the post-failure response. Details are given for both 2-D and 3-D implementations of the scaling algorithm. Copyright (c) 2015 John Wiley & Sons, Ltd. |
Databáze: | OpenAIRE |
Externí odkaz: |