Adaptive mesh refinement and coarsening for cohesive zone modeling of dynamic fracture

Park, Kyoungsoo; Paulino, Gláucio H.; Celes, Waldemar; Espinha, Rodrigo

doi:10.1002/nme.3163

Cited by 77 publications

(36 citation statements)

References 56 publications

Supporting

Mentioning

Contrasting

Unclassified

Order By: Relevance

“…Other challenges arise from the complexity of implementing computational procedures for fracture surface evolution, such as remeshing of complex and evolving fracture surfaces (Carter et al 2000), phase field methods (Hakim and Karma 2009), cohesive zone models (Park et al 2012), atomistic (Coffman et al 2008), peridynamics (Youn and Bobaru 2010) and other discrete approaches.…”

mentioning

confidence: 99%

Spiral to flat fracture transition for notched rods under torsional loading

Zehnder

Zella

2015

Int J Fract

View full text Add to dashboard Cite

The transition of the fracture surface in notched and unnotched PMMA rods under torsional loading is investigated experimentally. The experiments show that for crack initiation from a notch, resulting in dynamic fracture propagation, as the notch depth increases, the surface transitions from a spiral fracture to a nominally flat, but faceted surface. The facets coarsen and merge as the crack grows toward the center of the rod. The resulting fracture surfaces as well as the multiple internal micro-cracks are visualized using micro CT scanning.

show abstract

mentioning

confidence: 99%

Spiral to flat fracture transition for notched rods under torsional loading

Zehnder

Zella

2015

Int J Fract

View full text Add to dashboard Cite

show abstract

“…In this study, the Lagrange basis shape functions are used in conjunction with an existing finite element mesh. It should be noted that, for an element split of linear triangular elements, the strain energy is conserved when the Lagrange basis shape functions are utilized . This is because two split elements reproduce the displacement field of the original element.…”

Section: Element Splitting Schemementioning

confidence: 99%

“…It should be noted that, for an element split of linear triangular elements, the strain energy is conserved when the Lagrange basis shape functions are utilized. 30 This is because two split elements reproduce the displacement field of the original element. However, strain energy conservation is not guaranteed when node relocation is utilized.…”

Section: Element Splitting Schemementioning

confidence: 99%

Removing mesh bias in mixed‐mode cohesive fracture simulation with stress recovery and domain integral

Choi

Park

2019

Numerical Meth Engineering

Self Cite

View full text Add to dashboard Cite

Summary To remove mesh bias and provide an accurate crack path representation in mixed‐mode investigation, a novel stress recovery technique is proposed in conjunction with a domain integral and element splits. Based on a domain integral and stress recovery technique, a maximum strain energy release rate is estimated to determine a crack path direction. Then, for a given crack path direction, continuum elements are split, and a cohesive surface element is adaptively inserted. One notes that the proposed stress recovery technique provides a more accurate stress field than a standard stress evaluation procedure. The proposed computational framework is verified and validated by solving mode‐I and mixed‐mode examples. Computational results demonstrate that the domain integral with the stress recovery accurately evaluates a crack path, even with a lower‐quality mesh and under a biaxial stress state. Furthermore, the cohesive surface element approach, with the element split in conjunction with the stress recovery and the domain integral, predicts mixed‐mode fracture behaviors while removing mesh bias in the crack path representation. Additionally, the condition numbers of stiffness matrices are within the same order of magnitude during cohesive fracture simulation.

show abstract

“…To avoid this problem, alternative approaches have been developed allowing for a locally bounded refinement zone that follows the crack tip. Common strategies use adaptive hierarchical enrichment , self‐adaptive finite elements , adaptive h ‐refinement , or error‐driven adaptive re‐meshing . All these approaches significantly increase the accuracy of the approximation at the crack tip.…”

Section: Introductionmentioning

confidence: 99%

Multi‐level hp‐adaptivity for cohesive fracture modeling

Zander

Ruess

Bog

et al. 2016

Numerical Meth Engineering

View full text Add to dashboard Cite

Discretization induced oscillations in the load-displacement curve are a well known problem for simulations of cohesive crack growth with finite elements. The problem results from an insufficient resolution of the complex stress state within the cohesive zone ahead of the crack tip. This work demonstrates that the hp-version of the finite element method is ideally suited to resolve this complex and localized solution characteristic with high accuracy and low computational effort. To this end, we formulate a local and hierarchic mesh refinement scheme that follows dynamically the propagating crack tip. In this way, the usually applied static a priori mesh refinement along the complete potential crack path is avoided, which significantly reduces the size of the numerical problem. Studying systematically the influence of h-, p-and hp-refinement, we demonstrate why the suggested hp-formulation allows to capture accurately the complex stress state at the crack front preventing artificial snapthrough and snap-back effects. This allows to decrease significantly the number of degrees of freedom and the simulation runtime. Furthermore, we show that by combining this idea with the finite cell method, the crack propagation within complex domains can be simulated efficiently without resolving the geometry by the mesh.Keywords: cohesive fracture, automatic hp-adaptivity, arbitrary hanging nodes, dynamic meshes, finite cell method Changes to the manuscript resulting from remarks of Reviewer one are marked in red. Changes resulting from remarks of Reviewer two are marked in blue. Changes resulting from remarks of both reviewers are marked in green. Other textual changes are marked in brown.

show abstract

Adaptive mesh refinement and coarsening for cohesive zone modeling of dynamic fracture

Cited by 77 publications

References 56 publications

Spiral to flat fracture transition for notched rods under torsional loading

Spiral to flat fracture transition for notched rods under torsional loading

Removing mesh bias in mixed‐mode cohesive fracture simulation with stress recovery and domain integral

Multi‐level hp‐adaptivity for cohesive fracture modeling

Contact Info

Product

Resources

About