Adaptive isogeometric analysis for phase-field modeling of anisotropic brittle fracture

TY - JOUR

T1 - Adaptive isogeometric analysis for phase-field modeling of anisotropic brittle fracture

AU - Chen, Lin

AU - Li, Bin

AU - de Borst, René

N1 - Publisher Copyright: © 2020 John Wiley & Sons, Ltd.

PY - 2020/10/30

Y1 - 2020/10/30

N2 - The surface energy a phase-field approach to brittle fracture in anisotropic materials is also anisotropic and gives rise to second-order gradients in the phase field entering the energy functional. This necessitates C1 continuity of the basis functions which are used to interpolate the phase field. The basis functions which are employed in isogeometric analysis (IGA), such as nonuniform rational B-splines and T-splines naturally possess a higher order continuity and are therefore ideally suited for phase-field models which are equipped with an anisotropic surface energy. Moreover, the high accuracy of spline discretizations, also relative to their computational demand, significantly reduces the fineness of the required discretization. This holds a fortiori if adaptivity is included. Herein, we present two adaptive refinement schemes in IGA, namely, adaptive local refinement and adaptive hierarchical refinement, for phase-field simulations of anisotropic brittle fracture. The refinement is carried out using a subdivision operator and exploits the Bézier extraction operator. Illustrative examples are included, which show that the method can simulate highly complex crack patterns such as zigzag crack propagation. An excellent agreement is obtained between the solutions from global refinement and adaptive refinement, with a reasonable reduction of the computational effort when using adaptivity.

AB - The surface energy a phase-field approach to brittle fracture in anisotropic materials is also anisotropic and gives rise to second-order gradients in the phase field entering the energy functional. This necessitates C1 continuity of the basis functions which are used to interpolate the phase field. The basis functions which are employed in isogeometric analysis (IGA), such as nonuniform rational B-splines and T-splines naturally possess a higher order continuity and are therefore ideally suited for phase-field models which are equipped with an anisotropic surface energy. Moreover, the high accuracy of spline discretizations, also relative to their computational demand, significantly reduces the fineness of the required discretization. This holds a fortiori if adaptivity is included. Herein, we present two adaptive refinement schemes in IGA, namely, adaptive local refinement and adaptive hierarchical refinement, for phase-field simulations of anisotropic brittle fracture. The refinement is carried out using a subdivision operator and exploits the Bézier extraction operator. Illustrative examples are included, which show that the method can simulate highly complex crack patterns such as zigzag crack propagation. An excellent agreement is obtained between the solutions from global refinement and adaptive refinement, with a reasonable reduction of the computational effort when using adaptivity.

KW - T-splines

KW - adaptivity

KW - anisotropic surface energy

KW - isogeometric analysis

KW - phase-field

UR - http://www.scopus.com/inward/record.url?scp=85088263653&partnerID=8YFLogxK

U2 - 10.1002/nme.6457

DO - 10.1002/nme.6457

M3 - 文章

AN - SCOPUS:85088263653

SN - 0029-5981

VL - 121

SP - 4630

EP - 4648

JO - International Journal for Numerical Methods in Engineering

JF - International Journal for Numerical Methods in Engineering

IS - 20

ER -