The peridynamic theory is a nonlocal formulation of continuum mechanics based on integro-differential equations, devised to deal with fracture in solid bodies. In particular, the forces acting on each material point are evaluated as the integral of the nonlocal interactions with all the neighboring points within a spherical region, called "neighborhood". Peridynamic bodies are commonly discretized by means of a meshfree method into a uniform grid of cubic cells. The numerical integration of the nonlocal interactions over the neighborhood strongly affects the accuracy and the convergence behavior of the results. However, near the boundary of the neighborhood, some cells are only partially within the sphere. Therefore, the quadrature weights related to those cells are computed as the fraction of cell volume which actually lies inside the neighborhood. This leads to the complex computation of the volume of several cube-sphere intersections for different positions of the cells. We developed an innovative algorithm able to accurately compute the quadrature weights in 3D peridynamics for any value of the grid spacing (when considering fixed the radius of the neighborhood). Several examples have been presented to show the capabilities of the proposed algorithm. With respect to the most common algorithm to date, the new algorithm provides an evident improvement in the accuracy of the results and a smoother convergence behavior as the grid spacing decreases.

Accurate computation of partial volumes in 3D peridynamics

Scabbia, F
;
Zaccariotto, M;Galvanetto, U
2023

Abstract

The peridynamic theory is a nonlocal formulation of continuum mechanics based on integro-differential equations, devised to deal with fracture in solid bodies. In particular, the forces acting on each material point are evaluated as the integral of the nonlocal interactions with all the neighboring points within a spherical region, called "neighborhood". Peridynamic bodies are commonly discretized by means of a meshfree method into a uniform grid of cubic cells. The numerical integration of the nonlocal interactions over the neighborhood strongly affects the accuracy and the convergence behavior of the results. However, near the boundary of the neighborhood, some cells are only partially within the sphere. Therefore, the quadrature weights related to those cells are computed as the fraction of cell volume which actually lies inside the neighborhood. This leads to the complex computation of the volume of several cube-sphere intersections for different positions of the cells. We developed an innovative algorithm able to accurately compute the quadrature weights in 3D peridynamics for any value of the grid spacing (when considering fixed the radius of the neighborhood). Several examples have been presented to show the capabilities of the proposed algorithm. With respect to the most common algorithm to date, the new algorithm provides an evident improvement in the accuracy of the results and a smoother convergence behavior as the grid spacing decreases.
File in questo prodotto:
File Dimensione Formato  
s00366-022-01725-3.pdf

accesso aperto

Tipologia: Published (publisher's version)
Licenza: Creative commons
Dimensione 2.78 MB
Formato Adobe PDF
2.78 MB Adobe PDF Visualizza/Apri
Pubblicazioni consigliate

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/11577/3457572
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 14
  • ???jsp.display-item.citation.isi??? 7
  • OpenAlex ND
social impact