Block structured adaptive mesh refinement and strong form elasticity approach to phase field fracture with applications to delamination, crack branching and crack deflection
Vinamra Agrawal,Brandon Runnels
DOI: https://doi.org/10.1016/j.cma.2021.114011
IF: 6.588
2021-11-01
Computer Methods in Applied Mechanics and Engineering
Abstract:<p><span>Fracture is a ubiquitous phenomenon in most composite engineering structures, and is often the responsible mechanism for <a class="topic-link" href="/topics/engineering/catastrophic-failure">catastrophic failure</a><span><span>. Over the past several decades, many approaches have emerged to model and predict crack failure. The phase field method for fracture uses a surrogate damage field to model <a class="topic-link" href="/topics/engineering/crack-propagation">crack propagation</a>, eliminating the arduous need for explicit crack meshing. In this work a novel numerical framework is proposed for implementing hybrid phase field fracture in heterogeneous materials. The proposed method is based on the "reflux-free" method for solving, in strong form, the equations of </span><a class="topic-link" href="/topics/engineering/linear-elasticity">linear elasticity</a><span><span> on a block-structured <a class="topic-link" href="/topics/engineering/adaptive-mesh-refinement">adaptive mesh refinement</a><span> (BSAMR) mesh. The use of BSAMR enables highly efficient and scalable regridding, facilitates the use of temporal subcycling for explicit time integration, and allows for ultra-high refinement at crack boundaries with minimal computational cost. The method is applied to a variety of simple heterogeneous structures: laminates, wavy interfaces, and <a class="topic-link" href="/topics/mathematics/circular-inclusion">circular inclusions</a>. In each case a non-dimensionalized parameter study is performed to identify regions of behavior, varying both the geometry of the problem and the relative fracture </span></span><a class="topic-link" href="/topics/engineering/energy-release-rate">energy release rate</a>. In the laminate and wavy interface cases, regions of </span></span></span><a class="topic-link" href="/topics/engineering/delamination">delamination</a> and fracture correspond to simple analytical predictions. For the circular inclusions, the modulus ratio of the inclusion is varied as well as the delamination energy release rate and the problem geometry. In this case, a wide variety of behaviors was observed, including deflection, splitting, delamination, and pure fracture.</p>
mechanics,engineering, multidisciplinary,mathematics, interdisciplinary applications