The enrichment discussion of the CPDI2 publication* provided the following figure to illustrate the CPDI2 basis functions:
This figure, however, has the shortcoming of not clearly depicting the partition of unity property.
In hindsight, we should have used…
A more detailed and comprehensive version would expand the size of the enriched zone, and it would include more detailed annotation as follows:
As with any MPM method, the algorithm loops over particles. If a particle falls outside of an enriched zone, then its corners affect only regular CPDI grid basis functions. For enriched problems, the CPDI paper identifies three types of particle corners:
- a baseline (regular) corner is located outside the enriched zone. These are beyond the particles marked with triangles.
- an enriched corner is located on the interior of the enriched zone. These are marked with star symbols.
- a transition corner is on the boundary between the enriched and non-enriched zones. These are marked with triangles.
Then three types of particles may be identified:
- A non-enriched particle is located outside the enrichment zone. The basis functions over those particle domains look exactly as they would for an ordinary CPDI implementation, and they are handled the same as in a regular CPDI implementation.
- A fully enriched particle is located entirely inside the enriched zone and has only enriched corners. This particle affects only the enriched DOFs associated with its own corners. The calculation of the average of the shape function and its gradient over the particle domain proceeds the same as in the regular CPDI implementation, but with the sampled shape function values at the corners being the usual Kronecker deltas.
- A transition particle is located inside the enriched zone, but with at least one transition corner. These particles affect both ordinary grid nodes and enriched nodes. The ordinary grid contribution is handled similarly to a regular CPDI implementation except that the sampling of the standard grid shape functions is replaced with zero on the enriched corners. The contribution of the particle to the enriched DOFs is the same as it would be for a fully enriched particle.
*Sadeghirad, A., Brannon, R.M. and Guilkey, J.E. (2013), Second-order convected particle domain interpolation (CPDI2) with enrichment for weak discontinuities at material interfaces. Int. J. Numer. Meth. Engng., 95: 928–952. doi: 10.1002/nme.4526