Libraries-based explicit fracture and fragmentation framework

1. A computer-implemented method, comprising:

providing a numerical simulation mesh comprising a plurality of nodes; and

performing a numerical simulation of a solid material using the mesh and a user-selectable transition threshold associated with a cohesive zone model traction-separation relation; and

comparing local material states of the mesh to the user-selectable transition threshold during the performing the numerical simulation to determine whether to insert cohesive elements in the mesh;

wherein the cohesive elements are not included at an initiation of the numerical simulation.

2. The computer-implemented method of claim 1 , further comprising dynamically inserting cohesive elements in the mesh based on the comparison; and

balancing loads at the inserted cohesive elements.

3. The computer-implemented method of claim 1 , wherein the comparing comprises applying a discrete fracture criterion at a node of the plurality of nodes of the mesh based on the local material states of the mesh and the transition threshold.

4. The computer-implemented method of claim 3 , wherein the discrete fracture criterion defines when the cohesive elements are to be dynamically inserted and is a generalized criterion that can vary for different materials and that is a function of the Cauchy stress tensor and/or the strain tensor and the transition threshold.

5. The computer-implemented method of claim 3 , further comprising, for a met discrete fracture criterion, dynamically inserting cohesive elements in the mesh; and

applying a vector material dependent field function to determine a direction of a fracture from a fracture initiating node.

6. The computer-implemented method of claim 5 , further comprising splitting the fracture initiating node into at least two nodes and mapping the material states to an adjacent velocity node closest to the direction of the fracture.

7. The computer-implemented method of claim 5 , further comprising remeshing at least a portion of the mesh proximate the fracture initiating node to define a mesh shape oriented along the direction of the fracture; and

splitting the fracture initiating node into at least two nodes and mapping the material states to adjacent velocity nodes.

8. The computer-implemented method of claim 1 , wherein the user-selectable transition threshold is defined in a range between 0 and 1 and corresponds to a ratio of a stress and a fracture strength of the solid material.

9. The computer-implemented method of claim 1 , wherein the plurality of nodes comprises velocity nodes associated with material-specific fracture libraries and defining simplex shapes, and material nodes associated with the simplex shapes and associated with material libraries.

10. The computer-implemented method of claim 9 , wherein the performing the numerical simulation includes:

determining velocities and positions for the velocity nodes using one or more numerical analysis techniques;

determining material states for the material nodes based on the material libraries and the determined velocities and positions; and

determining material states for the velocity nodes based on the determined material states for the material nodes.

11. The computer-implemented method of claim 9 , wherein the material-specific fracture libraries and material libraries are user-selectable libraries.

12. The computer-implemented method of claim 11 , wherein the material-specific fracture libraries include fracture material libraries associated with discrete fracture propagation, and wherein the material libraries include continuum material libraries.

13. The computer-implemented method of claim 12 , wherein the continuum material libraries include libraries describing elasticity, plasticity, damage, strain rate, hardening, shock, or phase change.

14. The computer-implemented method of claim 1 , further comprising processing a solid material based on the performed numerical simulation, including one or more of surface material processing, subsurface material processing, oil drilling, fracking, well drilling, subsurface explosions, or mining operations, earthquake rupture simulations, or structural collapse.

15. The computer-implemented method of claim 1 , wherein the insertion of cohesive elements corresponds to a time-continuous transition between continua mechanics modeling and discontinua mechanics modeling;

wherein the local material states include a stress, a strain, a stress and strain, a stress and an additional other material state, a strain and an additional other material state, or a stress, strain and an additional other material state.

16. The computer-implemented method of claim 1 , further comprising removing cohesive elements from the mesh during the performing the numerical simulation based on a comparison of local material states of the mesh to a selected cohesive element removal threshold.

17. An apparatus, comprising:

at least one processor; and

memory coupled to the at least one processor and configured with instructions that, when executed by the at least one processor, cause the apparatus to:

perform a numerical simulation of a solid material using a numerical simulation mesh comprising a plurality of nodes, wherein the numerical simulation includes a user-selectable transition threshold associated with a cohesive zone model traction-separation relation that determines whether to insert cohesive elements in the mesh; wherein the numerical simulation includes: (1) comparing local material states of the mesh to the user-selectable transition threshold during the performing the numerical simulation to determine whether to insert the cohesive elements in the mesh; (2) providing smooth transition of the material states from continuum to discontinuum to ensure the dynamic equilibrium of the system, wherein the cohesive elements are not included at an initiation of the numerical simulation.

18. A computer-implemented method, comprising:

performing a numerical simulation of a solid material using a numerical simulation mesh comprising a plurality of nodes, wherein the numerical simulation includes a user-selectable transition threshold associated with a cohesive zone model traction-separation relation; and

comparing local material states including a strain state of the mesh to the transition threshold during the performing the numerical simulation to determine whether to dynamically insert cohesive elements in the finite element mesh; and

providing smooth transition of the material states from continuum to discontinuum to ensure the dynamic equilibrium of the system;

wherein the cohesive elements are not included at an initiation of the numerical simulation.

19. The computer-implemented method of claim 18 , further comprising associating a continuum material library with a node of the mesh, wherein the associated continuum material library uses a strain state to define a fracture initiation.