Variational Phase Field Modeling of Fracture Caused by Fluid Transport in Poroelastic Solids

The numerical modeling of failure mechanisms due to fracture based on sharp crack discontinuities is extremely demanding and suffers in situations with complex crack topologies. This drawback can be overcome by recently developed diffusive crack modeling concepts, which are based on the introduction of a crack phase field. Such an approach is conceptually in line with gradient-extended continuum damage models which include internal length scales. In this paper, we extend our recently outlined mechanical framework [1–3] towards the phase field modeling of fracture in the coupled problem of fluid transport in deforming porous media. Here, extremely complex crack patterns may occur due to drying or hydraulic induced fracture, the so called fracking. We develop new variational potentials for Biot-type fluid transport in porous media at finite deformations coupled with phase field fracture. It is shown, that this complex coupled multi-field problem is related to an intrinsic mixed variational principle for the evolution problem. This principle determines the rates of deformation, fracture phase field and fluid content along with the fluid potential. We develop a robust computational implementation of the coupled problem based on the potentials mentioned above and demonstrate its performance by the numerical simulation of complex fracture patterns. (© 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Phase-Field Modeling of Fracture in Anisotropic Media

In the recent years phase-field modeling of fracture has become a promising tool to describe complex crack patterns in all kinds of solid materials. Many of the models assume an isotropic material behavior, which of course is not a meaningful assumption for e.g. biological tissues such as arterial walls. Since the phase-field approach introduces an additional (smeared) phase describing the evolution of the crack, this method is well suited to be extended to anisotropic materials without thinking about an adaption of the discretization technique. Anisotropy can be incorporated in several ways, like by an extension of the surface energy, i.e. by making the energy release rate orientation dependent, as considered in [1]. Our ansatz is based on a pure geometrical approach, namely on an anisotropic formulation of the crack surface itself. Here, we will focus on transversely isotropic and cubically anisotropic solids, where the latter one makes the incorporation of the second gradient of the crack phase field necessary. At the end one numerical example is shown, which conceptually shows the influence of the anisotropy on the crack path. (© 2015 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Ductile failure with gradient plasticity coupled to the phase-field fracture at finite strains

Most metals fail in a ductile fashion, i.e, fracture is preceded by significant plastic deformation. The modeling of failure in ductile metals must account for complex phenomena at micro-scale, such as nucleation, growth and coalescence of micro-voids. In this work, we start with von-Mises plasticity model without considering void generation. The modeling of macroscopic cracks can be achieved in a convenient way by the continuum phase field approaches to fracture, which are based on the regularization of sharp crack discontinuities [1]. This avoids the use of complex discretization methods for crack discontinuities and can account for complex crack patterns. The key aspect of this work is the extension of the energetic and the stress-based phase field driving force function in brittle fracture to account for a coupled elasto-plastic response in line with our recent work [3]. We develop a new theoretical and computational framework for the phase field modeling of ductile fracture in elastic-plastic solids. To account for large strains, the constitutive model is constructed in the logarithmic strain space, which simplify the model equations and results in a formulation similar to small strains. We demonstrate the modeling capabilities and algorithmic performance of the proposed formulation by representative simulations of ductile failure mechanisms in metals. (© 2015 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Towards Phase Field Modeling of Ductile Fracture in Gradient-Extended Elastic-Plastic Solids

The modeling of failure in ductile metals must account for complex phenomena at a micro-scale as well as the final rupture at the macro-scale. Within a top-down viewpoint, this can be achieved by the combination of a micro-structure-informed elastic-plastic model with a concept for the modeling of macroscopic crack discontinuities. In this context, it is important to account for material length scales and thermo-mechanical coupling effects due to dissipative heating. This can be achieved by the construction of non-standard, gradient-enhanced models of plasticity with a full embedding into continuum thermodynamics [1,2]. The modeling of macroscopic cracks can be achieved in a convenient way by recently developed continuum phase field approaches to fracture based on regularized crack discontinuities. This avoids the use of complex discretization methods for crack discontinuities, and can account for complex crack patterns within a pure continuum formulation. Moreover, the phase field modeling of fracture is related to gradient theories of continuum damage mechanics, and fits nicely the structure of constitutive models for gradient plasticity. The main focus of this work is the extensions to gradient thermoplasticity and phase field formulation of ductile fracture, conceptually in line with the work [3]. (© 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A Phase Field Model for Ductile to Brittle Failure Mode Transition

The computational modeling of failure mechanisms in solids due to fracture based on sharp crack discontinuities suffers in dynamic problems with complex crack topologies. This can be overcome by a diffusive crack modeling based on the introduction of a crack phase field. We outline a conceptual framework for phase field models of crack propagation in brittle elastic and ductile elastic-plastic solids under dynamic loading and investigate the ductile to brittle failure mode transition observed in the experiment performed by Kalthoff and Winkeler [3]. We develop incremental variational principles and consider their numerical implementations by multi-field finite element methods. To this end, we define energy storage and dissipation functions for the plastic flow including the fracture phase field. The introduction of local history fields that drive the evolution of the crack phase field inspires the construction of robust operator split schemes. (© 2012 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A Phase-Field Model of Ductile Fracture at Finite Strains

This work outlines a variational-based framework for the phase field modeling of ductile fracture in elastic-plastic solids at large strains. The phase field approach regularizes sharp crack discontinuities within a pure continuum setting by a specific gradient damage model with geometric features rooted in fracture mechanics. Based on the recent works [1, 2], the phase field model of ductile fracture is linked to a formulation of gradient plasticity at finite strains in order to ensure the crack to evolve inside the plastic zones. The thermodynamic formulation is based on the definition of a constitutive work density function including the stored elastic energy and the dissipated work due to plasticity and fracture. The proposed canonical theory is shown to be governed by a rate-type minimization principle, which determines the coupled multi-field evolution problem. Another aspect is the regularization towards a micromorphic gradient plasticity-damage setting which enhances the robustness of the finite element formulation. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Phase Field Modeling of Crack Propagation at Large Strains with Application to Rubbery Polymers

The computational modeling of failure mechanisms in solids due to fracture based on sharp crack discontinuities exhibits drawbacks in situations with complex crack topologies. This drawback is overcome by diffusive crack modeling based on the introduction of a fracture phase field characterizing via an auxiliary variable the crack topology. In the following we extend recent advances in phase-field-type fracture based on operator split techniques, suggested in Miehe et al. [1], to the modeling of crack propagation in geometrically large deforming solids e.g. rubber-like materials. An extremely robust algorithmic treatment based on an operator split scheme is introduced consisting of three steps. Updating i) a local history-field containing the maximum reference energy, ii) the fracture phase field, and iii) the displacement field. We demonstrate the performance of proposed phase field formulation for largely deforming solids by means of a representative numerical example. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Phase Field Modeling of Brittle and Ductile Fracture

The phase field modeling of brittle fracture was a topic of intense research in the last few years and is now well-established. We refer to the work [1-3], where a thermodynamically consistent framework was developed. The main advantage is that the phase-field-type diffusive crack approach is a smooth continuum formulation which avoids the modeling of discontinuities and can be implemented in a straightforward manner by multi-field finite element methods. Therefore complex crack patterns including branching can be resolved easily. In this paper, we extend the recently outlined phase field model of brittle crack propagation [1-3] towards the analysis of ductile fracture in elastic-plastic solids. In particular, we propose a formulation that is able to predict the brittle-to-ductile failure mode transition under dynamic loading that was first observed in experiments by Kalthoff and Winkler [4]. To this end, we outline a new thermodynamically consistent framework for phase field models of crack propagation in ductile elastic-plastic solids under dynamic loading, develop an incremental variational principle and consider its robust numerical implementation by a multi-field finite element method. The performance of the proposed phase field formulation of fracture is demonstrated by means of the numerical simulation of the classical Kalthoff-Winkler experiment that shows the dynamic failure mode transition. (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Configurational-Force-Based Adaptive FE Solver for a Phase Field Model of Fracture

The computational modeling of failure mechanisms in solids due to fracture based on sharp crack discontinuities suffers in situations with complex crack topologies. This can be overcome by diffusive crack modeling, based on the introduction of a crack phase field as outlined in [1, 2]. Following these formulations, we outline a thermodynamically consistent framework for phase field models of crack propagation in elastic solids, develop incremental variational principles and, as an extension to [1, 2], consider their numerical implementations by an efficient h-adaptive finite element method. A key problem of the phase field formulation is the mesh density, which is required for the resolution of the diffusive crack patterns. To this end, we embed the computational framework into an adaptive mesh refinement strategy that resolves the fracture process zones. We construct a configurational-force-based framework for h-adaptive finite element discretizations of the gradient-type diffusive fracture model. We develop a staggered computational scheme for the solution of the coupled balances in physical and material space. The balance in the material space is then used to set up indicators for the quality of the finite element mesh and accounts for a subsequent h-type mesh refinement. The capability of the proposed method is demonstrated by means of a numerical example. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Phase field modeling of complex crack patterns in multi-physics problems

Recently developed continuum phase field models for brittle fracture show excellent modeling capability in situations with complex crack topologies including branching in the small and large strain applications. This work presents a generalization towards fully coupled multi-physics problems at large strains. A modular concept is outlined for the linking of the diffusive crack modeling with complex multi field material response, where the focus is put on the model problem of finite thermo-elasticity. This concerns a generalization of crack driving forces from the energetic definitions towards stress-based criteria, the constitutive modeling of degradation of non-mechanical fluxes on generated crack faces. Particular assumptions are made on the generation of convective heat exchanges approximating surface load integrals of the sharp crack approach by distinct volume integrals. The coupling effect is also shown in generation of cracks due to thermally induced stress states. We finally demonstrate the performance of the phase field formulation of fracture at large strains by means of representative numerical examples. (© 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A Phase Field Model for Three-Dimensional Dynamic Fracture and its Efficient Numerical Implementation

The numerical modeling of dynamic failure mechanisms in solids due to fracture based on sharp crack discontinuities suffers in situations with complex crack topologies and demands the formulation of additional branching criteria. This drawback can be overcome by a diffusive crack modeling, which is based on the introduction of a crack phase field. We focus on the extension of a recently developed phase field model for fracture from the quasi-static setting towards the dynamic setting. It is obtained by taking into account inertial terms and associated dynamic integrators. The introduction of a history field, containing a maximum fracture-driving energy, provides a very transparent representation of the balance equation that governs the diffusive crack topology. In particular, it allows for the construction of an extremely robust operator split technique. In a subsequent step, the proposed model is extended to three dimensional problems. The dynamic treatment opens the door to the analysis of complex fracture phenomena like multiple crack branching and crack arrest. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Phase Field Modeling of Fracture in Plates and Shells

The numerical modeling of failure mechanisms in plates and shells due to fracture based on sharp crack discontinuities is extremely demanding and suffers in situations with complex crack topologies. This drawback can be overcome by a diffusive crack modeling, which is based on the introduction of a crack phase field. In this paper, we extend ideas recently outlined in [1, 2] towards the phase field modeling of fracture in dimension-reduced continua with application to Kirchhoff plates and shells. The introduction of history fields, containing the maximum reference energy obtained in history, provides a very transparent representation of the coupled balance equations and allows the construction of an extremely robust operator split technique. The performance of the proposed models is demonstrated by means of representative numerical examples. (© 2012 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Phase field simulation of thermomechanical fracture

In Francfort and Marigo's variational free-discontinuity formulation of brittle fracture [1] cracking is regarded as an energy minimization process, where the total energy is minimized with respect to any admissible crack set and displacement field. No additional criterion is needed to determine crack paths, branching of cracks and crack initiations. However, a direct discretization of the model is faced with significant technical problems, as it involves minimizations in a set of possibly discontinuous functions. A regularized version of the model has been introduced by Bourdin [2] and based on this, we use the concept of a continuum phase field model to simulate cracking processes. Cracks are indicated by the order parameter of the phase field model and cracking can be regarded as a phase transition problem. Additionally, introducing the heat equation into the model, it is capable to also take account of crack propagation due to thermal stresses. In the numerical implementation, crack parameter as well as temperature are treated as additional degrees of freedom and the coupled field equations are solved using the finite element method together with an implicit time integration scheme. (© 2009 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A novel treatment of crack boundary conditions in phase field models of fracture

In the phase field approach for fracture an additional scalar field is introduced in order to describe the state of the material between intact and fully broken. So far, for the loading dependent degradation of stiffness (damage) either the volumetric-deviatoric split of strain [1, 2] or the spectral decomposition [3, 4] is used. In contrast to such an isotropic degradation of stiffness, the fully broken state represents a crack with a particular orientation. Both aforementioned approaches do not take the crack orientation into account. This may lead to the violation of the crack boundary conditions. In order to satisfy these conditions the phase field approach is modified here by taking the orientation of the crack into account. (© 2015 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A phase field model for fracture

The variational formulation of brittle fracture as formulated for example by Francfort and Marigo in [1], where the total energy is minimized with respect to any admissible crack set and displacement field, allows the identification of crack paths, branching of preexisting cracks and even crack initiation without additional criteria. For its numerical treatment a continuous approximation of the model in the sense of Γ-convergence has been presented by Bourdin in [2]. In the regularized Francfort–Marigo model cracks are represented by an additional field variable (secondary variable) s∈[0,1] which is 0 if the material is cracked and 1 if it is undamaged. In this work, we reinterpret the crack variable as a phase field order parameter and address cracking as a phase transition problem. The crack growth is governed by the evolution equation of the order parameter which resembles the Ginzburg–Landau equation. The numerical treatment is done by finite elements combined with an implicit Euler scheme for the time integration. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Mathematical modelling of bottom evolution by particle deposition and resuspension

A two-dimensional numerical model for the evolution of a bottom due to particle deposition and resuspension by a fluid flow is here presented. A computational fluid dynamic approach is used to calculate the flow field and a Lagrangian particle tracking technique is applied to solve the dispersed phase. The evolution of the lower boundary is simulated taking into account the mass conservation of the solid phase and the geotechnical properties of the granular material. The model is characterized by two important features. First, fluid dynamics are coupled with the bottom evolution due to particle deposition and resuspension. This permits to use the model to simulate complex flow fields as well as complex time-evolving geometries. Second, the dispersed phase is calculated by a Lagrangian approach, which retains the discrete information of the individual particles of the granular bottom which may be of interest for some industrial processes (coating) and environmental flows (sediment stratification). First consistency checks have been performed for some deposition and resuspension test cases with fluid at rest. The model has also been tested by comparison with a physical experiment of deposition inside a cavity. Finally, as an example of possible applications of industrial and environmental interest, the model has been applied to investigate particle deposition in rectangular cavities and the evolution of a sand heap by a fluid flow.