Search

The phase-field model for fracture, despite its popularity and ease of implementation comes with its set of computational challenges. They are the non-convex energy functional, variational inequality due to fracture irreversibility, the need for extremely fine meshes to resolve the fracture. In this manuscript, the focus is on the numerical treatment of variational inequality. In this context, the popular history-variable approach suffers from variationally inconsistency and non-quantifiable nature of the error introduced. A better alternative, the penalisation approach, has the potential to render the stiffness matrix ill-conditioned. In order to circumvent both aforementioned issues, a micromorphic approach towards phase-field fracture modelling is proposed in this manuscript. Within this approach, a micromorphic extension of the energy functional is carried out. This transforms the phase-field into a local variable, while introducing a micromorphic variable that regularises the fracture problem. This reduction the regularity requirements for the phase-field enables an easier implementation of the fracture irreversibility constraint through simple 'max' operation, with system level precision. Numerical experiments carried out on benchmark brittle and quasi-brittle problems demonstrate the applicability and efficacy of the proposed model for a wide range of fracture problems.

In this manuscript, the phase-field fracture irreversibility constraint is transformed into an equality-based constraint using the slack variable approach. The equality-based fracture irreversibility constraint is then introduced in the phase-field fracture energy functional using the Lagrange Multiplier Method and the Penalty method. Both methods are variationally consistent with the conventional variational inequality phase-field fracture problem, unlike the history-variable approach. Thereafter, numerical experiments are carried out on benchmark problems in brittle and quasi-brittle fracture to demonstrate the efficacy of the proposed method.

The phase-field fracture free-energy functional is non-convex with respect to the displacement and the phase field. This results in a poor performance of the conventional monolithic solvers like the Newton-Raphson method. In order to circumvent this issue, researchers opt for the alternate minimization (staggered) solvers. Staggered solvers are robust for the phase-field based fracture simulations as the displacement and the phase-field sub-problems are convex in nature. Nevertheless, the staggered solver requires very large number of iterations (of the order of thousands) to converge. In this work, a robust monolithic solver is presented for the phase-field fracture problem. The solver adopts a fracture energy-based arc-length method and an adaptive under-relaxation scheme. The arc-length method enables the simulation to overcome critical points (snap-back, snap-through instabilities) during the loading of a specimen. The use of an under-relaxation scheme stabilizes the solver by preventing the divergence due to an ill-behaving stiffness matrix. The efficiency of the proposed solver is further amplified with an adaptive mesh refinement scheme based on PHT-splines within the framework of isogeometric analysis. The numerical examples presented in the manuscript demonstrates the efficacy of the solver. All the codes and data-sets accompanying this work will be made available on GitHub (https://github.com/rbharali/IGAFrac).