Phase Field Modeling Of Fracture Research Articles

Mechanical property evaluation of 3D multi-phase cement paste microstructures reconstructed using generative adversarial networks

This study proposes an artificial intelligence based framework for reconstructing the 3D multi-phase cement paste microstructure to evaluate its mechanical properties using simulation. The reconstruction of cement paste microstructures is performed using modified generative adversarial networks (GANs) based on microstructural images from micro-CT. For computational efficiency, 2D microstructures are first reconstructed and then extended to 3D microstructures. The reconstructed microstructures exhibit the same microstructural features as the original microstructures when characterized by probability functions. Mechanical properties such as stiffness and tensile strength are evaluated for the original and reconstructed specimens using a phase-field fracture model, and similar behaviors are observed. The results confirm that the reconstructed virtual microstructures can be used to supplement the real microstructures in evaluating the mechanical properties of 3D multi-phase cement paste. This approach thus provides a critical element of a data-driven approach to correlating its microstructure and properties.

Modelling and simulation of natural hydraulic fracturing applied to experiments on natural sandstone cores

AbstractUnder in-situ conditions, natural hydraulic fractures (NHF) can occur in permeable rock structures as a result of a rapid decrease of pore water accompanied by a local pressure regression. Obviously, these phenomena are of great interest for the geo-engineering community, as for instance in the framework of mining technologies. Compared to induced hydraulic fractures, NHF do not evolve under an increasing pore pressure resulting from pressing a fracking fluid in the underground but occur and evolve under local pore-pressure reductions resulting in tensile stresses in the rock material. The present contribution concerns the question under what quantitative circumstances NHF emerge and evolve. By this means, the novelty of this article results from the combination of numerical investigations based on the Theory of Porous Media with a tailored experimental protocol applied to saturated porous sandstone cylinders. The numerical investigations include both pre-existing and evolving fractures described by use of an embedded phase-field fracture model. Based on this procedure, representative mechanical and hydraulic loading scenarios are simulated that are in line with experimental investigations on low-permeable sandstone cylinders accomplished in the Porous Media Lab of the University of Stuttgart. The values of two parameters, the hydraulic conductivity of the sandstone and the critical energy release rate of the fracture model, have turned out essential for the occurrence of tensile fractures in the sandstone cores, where the latter is quantitatively estimated by a comparison of experimental and numerical results. This parameter can be taken as reference for further studies of in-situ NHF phenomena and experimental results.

A micromagnetic-mechanically coupled phase-field model for fracture and fatigue of magnetostrictive alloys

Magnetostrictive alloys are usually brittle materials with micromagnetic structures. Their structural reliability and durability depend on the complex micromagnetic-mechanical coupling at smaller length scales encompassing the evolution of micromagnetic structures. Herein we propose a micromagnetic-mechanically coupled phase-field model for fracture and fatigue behavior of magnetostrictive alloys with evolution of the micromagnetic structure. The thermodynamically-consistent model is derived from microforce theory, laws of thermodynamics, and Coleman–Noll analysis. The evolution of crack phase-field and magnetization-vector order parameters that are fully coupled is governed by history field dependent Allen–Cahn and Landau–Lifshitz–Gilbert equations, respectively. The model is extended to fatigue by introducing a degradation prefactor for the fracture energy as a function of positive elastic energy. One-dimensional analyses are then presented to anatomize the crack driving forces in terms of fully coupled micromagnetic-mechanical and pure mechanical driving force. We demonstrate the model capabilities by finite-element numerical studies on the micromagnetic domain evolution during the crack propagation and the influence of external magnetic field for type-I, type-II, and three-point bending fracture, as well as for the fracture of a single-edge notched specimen with an elliptical inclusion. The simulation result shows that depending on how micromagnetic domains are switched under micromagnetic-mechanical coupling, the magnetic field can enhance or decrease the critical load. In the presence of inclusion with larger fracture toughness, a crack is found to nucleate in the tri-junction of multi-domain micromagnetic structure owing to the high elastic strain around the tri-junction point. It is further found that a suitable magnetic field promoting magnetization-vector rotation around the crack tip could remarkably improve the fracturing load and fatigue life. The results demonstrate the model promising for the study of micromagnetic-mechanically coupled fracture and fatigue in magnetostrictive alloys.

Phase-field ductile fracture simulations of thermal cracking in additive manufacturing

We present a multiphysics phase-field fracture model for thermo-elasto-plastic solids in the context of finite deformation and apply it to simulate the hot cracking phenomenon during metal additive manufacturing. The model is derived in a thermodynamically consistent manner, with the intercoupling mechanisms among elastoplasticity, phase-field crack and heat transfer comprehensively considered. It involves particularly coupled parameters among these materials physics, e.g. plasticity-dependent degradation function and fracture toughness, damage-dependent yield surface and thermal properties, and temperature-dependent elastoplastic properties and fracture strength. The finite element implementation of the coupled phase-field model is benchmarked with simulations results of a tensile test of an I-shape specimen, encompassing elastoplasticity, hardening, necking, crack initiation and propagation, in contrast to the related experimental results. The validated model is further employed to simulate the multiphysics hot cracking phenomenon in additive manufacturing in the context of both the effective powder-bed model and the powder-resolved model thanks to prior non-isothermal phase-field powder-bed-fusion simulations. Simulation results reveal certain key features of the hot crack and its dependency on process parameters like beam power and scan speed, which are helpful for the fundamental understanding of crack formation mechanisms and process optimization.

Domain decomposition methods and acceleration techniques for the phase field fracture staggered solver

AbstractThe phase field modeling of fracture is able to simulate the nucleation and the propagation of complex crack patterns. However, the relatively small internal lengths that are required usually lead to very fine meshes and high computational costs, especially for three‐dimensional applications. In the present work, additional cost also comes from the implicit dynamics regularization of unstable crack propagations which potentially leads to a large variation of time steps when switching from quasi‐static to dynamic regimes. To reduce the time to solution in this context, this study proposes a domain decomposition framework and acceleration techniques for the phase field fracture staggered solver. The displacement subproblem and the phase field one are solved with parallel domain decomposition solvers. Dual domain decomposition methods provide low cost preconditioner well adapted to the phase field subproblem. For displacement subproblems undergoing unstable crack propagations, primal domain decomposition methods are preferred to be less sensitive to the treatment of floating substructures. Preconditioners performances are assessed and scalability studies over academic test cases, up to 324 subdomains, are presented. Finally, the robustness of the approach is illustrated on two semi‐industrial simulations.

A phase-field fracture model for creep-fatigue behavior

A phase-field fracture model for creep-fatigue behavior is proposed. The fatigue and creep behavior are simulated by modifying the free energy functional through a total degradation function that diminishes fracture toughness with the accumulation of the relevant history variable. Specifically, strain energy density and creep strain serve as the chosen fatigue and creep history variables, respectively. To model the creep-fatigue interaction, we propose several total degradation function formulations comprising fatigue and creep degradation functions. Numerical experiments demonstrate the capability of our model to reproduce the strain-lifetime relationship and crack growth behavior under creep-fatigue loading conditions with varying hold times. Moreover, we simulated examples of plates with holes, demonstrating the capability of our model for complex crack paths.

Phase-field model for 2D cohesive-frictional shear fracture: An energetic formulation

This paper presents a phase-field model for cohesive-frictional shear fracture. The model is derived based on two energetic principles: the energy conservation law and a variational inequality of virtual work that serves as a stability condition. We show that all the governing equations for frictional fracture can be obtained from the above two principles, including the equilibrium condition, the phase-field evolution law and, most importantly, the yield function and flow rule for the plasticity-like frictional slip. The information of crack direction is naturally included in the flow rule. Theoretical and numerical results show that the proposed model is faithfully consistent with the Mohr–Coulomb theory of frictional failure in terms of the crack nucleation stress and direction. In addition, the presented phase-field model has an explicit frictional cohesive law.

A phase-field model for the brittle fracture of Euler–Bernoulli beams coupling stretching and bending

Damage gradient models seek to simulate fracture mechanics by modulating the material stiffness. Within this framework, a singular scalar field representing damage is commonly utilized to globally decrease elastic energy. However, when considering structural models like beams and plates, this approach often fails to adequately capture important aspects due to the interaction between stretching and bending contributions. In this work, we propose a model for planar Euler–Bernoulli beams that incorporates the following features: firstly, we utilize two phase-field damage parameters to describe material damaging, specifically addressing the ‘erosion’ occurring above and below the original, undamaged beam; secondly, we assume a simple linear dependence of the axial stress field on the thickness coordinate, along with linear dependence on the axial force and bending moment. By appropriately identifying the constitutive response, our model effectively considers the coupling between stretching and bending induced during through-the-thickness damage, demonstrating good agreement with 2D observations.

Modeling dynamic crack growth in quasicrystals: Unraveling the role of phonon–phason coupling

In this work, we employed the phase-field fracture (PFF) model for dynamic brittle fracture in quasicrystals (QCs). The dynamic PFF formulation is derived using both elastodynamics and elasto-hydrodynamic theory. The temporal discretization is achieved using an implicit generalized-α method, which is unconditionally stable for a particular choice of constants. The model incorporates two distinct crack driving forces: (a) the contribution of elastic strain energy from both phonon and phason fields and (b) the contribution solely from the elastic energy associated with the phonon field. To validate the implementation of QCs, we have conducted benchmarking against recent results in the literature, particularly in the context of quasi-static loading. We addressed various paradigmatic case studies to demonstrate the dynamic crack nucleation and propagation in planar 2D decagonal QCs, with a primary emphasis on understanding the interplay between phonon and phason fields. Phason walls, representing low-energy crack paths resulting from atomic rearrangements, play a significant role in crack propagation by introducing an additional energy component. Notably, increasing the phonon–phason coupling constant expedites both crack initiation and propagation. Furthermore, variations in the phonon–phason coupling parameter result in different crack trajectories, observed in both uniaxial and biaxial loading scenarios. The developed code can be downloaded from https://github.com/Hirshikesh/Quasicrystal.git.

Voronoi-based 3D phase field model and numerical simulations of granite crack propagation

The distribution of mineral grains within underground rocks and materials is crucial for understanding geological processes and the stability of engineering structures. However, numerical simulation methods that consider 3D mineral particles are still immature. This study presents a 3D phase-field fracture model of mineral grains based on Voronoi polygons, making some beneficial attempts in this regard. The phase-field method does not rely on mesh reconstruction or crack path tracking and can simulate complex crack propagation behavior. Meanwhile, the Voronoi polygon can better describe the microscopic structure of rock with grain boundaries. Rooted in Voronoi polygons, this model efficiently captures the complex geometric shapes of mineral grains and accurately characterizes their positions, shapes, and sizes. Therefore, this study combines Voronoi polygons to establish a 3D phase-field model. Through two numerical simulation cases, the model demonstrates its effectiveness in simulating the impact of three-dimensional mineral grains on crack propagation behavior, showing that cracks tend to propagate at the boundaries or inside softer mineral grains. The modeling approach proposed in this paper has good reference value for numerical simulation studies considering the microscopic structure of rocks.

On the energy decomposition in variational phase-field models for brittle fracture under multi-axial stress states

AbstractPhase-field models of brittle fracture are typically endowed with a decomposition of the elastic strain energy density in order to realistically describe fracture under multi-axial stress states. In this contribution, we identify the essential requirements for this decomposition to correctly describe both nucleation and propagation of cracks. Discussing the evolution of the elastic domains in the strain and stress spaces as damage evolves, we highlight the links between the nucleation and propagation conditions and the modulation of the elastic energy with the phase-field variable. In light of the identified requirements, we review some of the existing energy decompositions, showcasing their merits and limitations, and conclude that none of them is able to fulfil all requirements. As a partial remedy to this outcome, we propose a new energy decomposition, denoted as star-convex model, which involves a minimal modification of the volumetric-deviatoric decomposition. Predictions of the star-convex model are compared with those of the existing models with different numerical tests encompassing both nucleation and propagation.

Optimizing integration point density for exponential finite element shape functions for phase-field modeling of fracture in functionally graded materials

In the realm of fracture mechanics and phase-field modeling, the utilization of exponential finite element (EFE) shape functions has shown promise in accurately predicting fracture responses, particularly in scenarios with intricate crack propagation paths. However, the computational complexity associated with EFE shape functions, including higher integration schemes and orientation requirements, prompts the need for optimization. This study delves into the critical aspect of integration point density and aims to determine the ideal balance between accuracy and computational efficiency in the context of EFE shape functions. The research focuses on functionally graded materials, which inherently exhibit spatial variations in solution quantities due to their graded material system. Through a systematic investigation, the study aims to identify the optimal number of integration points required to maximize the computational advantages offered by EFE shape functions while minimizing the associated computational resources. By conducting a comprehensive analysis across various loading scenarios, such as tension and shear, this study intends to establish a quantitative relationship between integration point density and prediction accuracy for fracture responses.

Phase-field modeling of fracture for ferromagnetic materials through Maxwell’s equation

Electro-active materials are classified as electrostrictive and piezoelectric materials. They deform under the action of an external electric field. Piezoelectric material, as a special class of active materials, can produce an internal electric field when subjected to mechanical stress or strain. In return, there is the converse piezoelectric response, which expresses the induction of the mechanical deformation in the material when it is subjected to the application of the electric field. This work presents a variational-based computational modeling approach for failure prediction of ferromagnetic materials. In order to solve this problem, a coupling between magnetostriction and mechanics is modeled, then the fracture mechanism in ferromagnetic materials is investigated. Furthermore, the failure mechanics of ferromagnetic materials under the magnetostrictive effects is studied based on a variational phase-field model of fracture. Phase-field fracture is numerically challenging since the energy functional may admit several local minima, imposing the global irreversibility of the fracture field and dependency of regularization parameters related discretization size. Here, the failure behavior of a magnetoelastic solid body is formulated based on the Helmholtz free energy function, in which the strain tensor, the magnetic induction vector, and the crack phase-field are introduced as state variables. This coupled formulation leads to a continuity equation for the magnetic vector potential through well-known Maxwell’s equations. Hence, the energetic crack driving force is governed by the coupled magneto-mechanical effects under the magneto-static state. Several numerical results substantiate our developments.

A computational meshfree RPIM approach for phase-field modeling of brittle fracture

A computational meshfree RPIM approach for phase-field modeling of brittle fracture

Damage Mechanics Challenge: Predictions based on the phase field fracture model

In this work, we describe our contribution to the Purdue-SANDIA-LLNL Damage Mechanics Challenge. The phase field fracture model is adopted to blindly estimate the failure characteristics of the challenge test, an unconventional three-point bending experiment on an additively manufactured rock resembling a type of gypsum. The model is formulated in a variationally consistent fashion, incorporating a volumetric–deviatoric strain energy decomposition, and the numerical implementation adopts a monolithic unconditionally stable solution scheme. Our focus is on providing an efficient and simple yet rigorous approach capable of delivering accurate predictions based solely on physical parameters. Model inputs are Young’s modulus E, Poisson’s ratio ν, toughness Gc and strength σc (as determined by the choice of phase field length scale ℓ). We show that a single mode I three-point bending test is sufficient to calibrate the model, and that the calibrated model can then reliably predict the force versus displacement responses, crack paths and surface crack morphologies of more intricate three-point bending experiments that are inherently mixed-mode. Importantly, our peak load, crack trajectory and crack surface morphology predictions for the challenge test, submitted before the experimental data was released, show a remarkable agreement with experiments. The characteristics of the challenge, and how changes in these can impact the predictive abilities of phase field fracture models, are also discussed.

A material energy–momentum flux-driven phase field fracture mechanics model

This work introduces a novel phase-field solution for simulating and predicting fractures in elastic solids. In this phase-field fracture model, the crack growth is driven by the material energy–momentum flux, or the configurational force, inspired by a null divergence conservation law and its variational theory, as opposed to employing the commonly-used strain energy degradation-based crack-driving force. By doing so, the crack growth or material configuration change is solely determined by the physics-based Eshelby energy–momentum tensor density or the Rice J-integral flux, which aligns with the energy release criterion of the celebrated Griffith–Eshelby–Rice (GER) theory. The proposed method not only preserves the integrity of the balance of linear momentum principle at the crack tip region but also provides a correct stress asymptotic field in front of the crack tip, while accurately capturing the crack growth.The key idea of this approach is to use the monotonically increasing and irreversible phase field as the marker of the configurational change and form a phase field variational principle that couples the strain energy density with the configurational strain energy density. By doing so, we only damage the displacement field as a form of diffused crack, as opposed to damaging the stress field in the existing phase-field method. Moreover, from numerical simulations, we discovered two different forms of energy release rates: one drives crack growth by enforcing the crack-tip traction loading condition, and one drives crack growth by enforcing the crack-tip deformation condition, demonstrating the ability of the proposed method to naturally capture two distinct modes of fractures due to two different near-crack tip loading conditions. These advancements and contributions reveal previously unknown insights into fracture mechanics.

An efficient and robust staggered scheme based on adaptive time field for phase field fracture model

The standard staggered scheme in phase field model is computationally slow, so that a large amount of iterations are required. In order to improve the computational efficiency, a novel efficient and robust solution scheme based on adaptive time field is proposed in this paper. This scheme is inspired by the adaptive mesh method for space field, and an error indicator based on phase field and history field variables is established, where the time field is coupled with the phase field and energy field in the space field, so that the time field changes adaptively with the computational iterative process, which improves the convergence rate and reduces the time cost. The numerical results show that the presented novel scheme not only significantly improves the computational efficiency, but also has accurate computational results. Compared with the standard staggered scheme, the time cost of the presented scheme is less than 4% of that of the standard scheme. At the same time, the computational efficiency of the presented solver is faster than the current monolithic BFGS scheme, and it has excellent convergence characteristics.

An elastoplastic phase-field model for quasi-static fracture of nickel-based super-alloys

In the present study, an elastoplastic phase-field model of quasi-static fracture in ductile materials is proposed in the variational framework for J2 plasticity with isotropic hardening, which is suitable to describe the quasi-static behavior of metals as investigated in the performed experiments. These contributions include: (1) the free energy functions for coupling elastic response, plastic yielding and damage evolution are established. (2) The new elastic and plastic energy degradation functions are constructed to quantitatively describe the relationship between energy release and phase-field evolution of elastoplastic materials. (3) Damage evolution and plastic yielding criteria are derived. (4) From a numerical point of view, we derive the governing equations and the corresponding weak forms and the overall solution procedure for the phase-field model is given via the use of a return-mapping algorithm. This phase-field model was validated by a series of tensile experiments on Inconel 718 nickel-based super-alloys standard specimens. In order to compare the simulation results with the experimental results more comprehensively, the digital image correlation (DIC) technique is applied to experimentally investigate the specimen deformation information. In addition, to verify the potential of the model to capture complex cracks, we performed Nooru-Mohamed tests. The numerical simulation results are in good agreements with the results of previous experimental work.

The phase-field fracture model enriched by interpolation cover functions for brittle fracture problems

In this study, interpolation cover functions were introduced into the phase-field fracture model (PFM) to improve its convergence and efficiency. The traditional PFM has proven to be successful in simulating brittle fracture problems. However, one disadvantage is that it requires a relatively dense mesh. The PFM model enriched by interpolation cover functions (EPFM) was used for numerical fracture analysis under a relatively coarse mesh. The rank deficiency problem caused by the introduction of interpolation cover functions was resolved using the rank deficiency repair technique. Higher gradients were captured by introducing interpolation cover functions without changing the mesh structure. The numerical examples showed the improvement in convergence. The high efficiency of the proposed EPFM was demonstrated through a series of numerical examples. The numerical results of the EPFM and PFM under different length-scale parameters and meshes were also compared. The results showed that the proposed EPFM is an excellent alternative for evaluating engineering fracture problems.

Rate- and temperature-dependent ductile-to-brittle fracture transition: Experimental investigation and phase-field analysis for toffee

The mechanical behaviour of many materials, including polymers or natural materials, significantly depends on the rate of deformation. As a consequence, a rate-dependent ductile-to-brittle fracture transition may be observed. For toffee-like caramel, this effect is particularly pronounced. At room temperature, this confectionery may be extensively deformed at low strain rates, while it can behave highly brittle when the rate of deformation is raised. Likewise, the material behaviour does significantly depend on temperature, and even a slight cooling may cause a significant embrittlement.In this work, a thorough experimental investigation of the rate-dependent deformation and fracture behaviour is presented. In addition, the influence of temperature on the material response is studied. The experimental results form the basis for a phase-field modelling of fracture. In order to derive the governing equations of the model, an incremental variational principle is introduced. By means of the validated model, an analysis of the experimentally observed ductile-to-brittle fracture transition is performed. In particular, the coupling between the highly dissipative deformation behaviour of the bulk material and the rate-dependent fracture resistance is discussed.