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)     

On Configurational-Force-Driven Crack Propagation and Adaptive Mesh Refinement in Finite Inelasticity

We introduce a consistent variational framework for inelasticity at finite strains, yielding dual balances in physical and material space as the Euler equations. The formulation is employed for the simultaneous usage of configurational forces as both driving forces for crack propagation as well as h-adaptive mesh refinement. The theoretical basis builds upon a global balance of internal and external power, where the mechanical response is exclusively governed by two scalar functions, the free energy function and a dissipation potential. The resulting variational structure is exploited in the context of fracture mechanics and yields evolution equations for internal variables. In the discrete setting, we present a geometry model fully separated from the finite element mesh structure that represents structural changes of the material configuration due to crack propagation. Advanced meshing algorithms provide an optimal discretization at the crack tip. Local and global criteria are obtained via error estimators based on configurational forces being interpreted as indicators of an energetic misfit due to an insufficient discretization. The numerical handling is decomposed into a staggered algorithm scheme for the dual set of equilibrium equations in material and physical space and efficient mesh generation tools. Exemplary numerical examples are considered to illustrate the method and to underline the effects of inelastic material behaviour in the presented context. (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Crack growth modeling via 3D automatic adaptive mesh refinement based on modified-SPR technique

In this paper, the three-dimensional automatic adaptive mesh refinement is presented in modeling the crack propagation based on the modified superconvergent patch recovery technique. The technique is developed for the mixed mode fracture analysis of different fracture specimens. The stress intensity factors are calculated at the crack tip region and the crack propagation is determined by applying a proper crack growth criterion. An automatic adaptive mesh refinement is employed on the basis of modified superconvergent patch recovery (MSPR) technique to simulate the crack growth by applying the asymptotic crack tip solution and using the collapsed quarter-point singular tetrahedral elements at the crack tip region. A-posteriori error estimator is used based on the Zienkiewicz–Zhu method to estimate the error of fracture parameters and predict the crack path pattern. Finally, the efficiency and accuracy of proposed computational algorithm is demonstrated by several numerical examples.     

Fractal materials, beams, and fracture mechanics

Continuing in the vein of a recently developed generalization of continuum thermomechanics, in this paper we extend fracture mechanics and beam mechanics to materials described by fractional integrals involving D, d and R. By introducing a product measure instead of a Riesz measure, so as to ensure that the mechanical approach to continuum mechanics is consistent with the energetic approach, specific forms of continuum-type equations are derived. On this basis we study the energy aspects of fracture and, as an example, a Timoshenko beam made of a fractal material; the local form of elastodynamic equations of that beam is derived. In particular, we review the crack driving force G stemming from the Griffith fracture criterion in fractal media, considering either dead-load or fixed-grip conditions and the effects of ensemble averaging over random fractal materials.     

A coupled stress and energy failure model for crack initiation in arbitrarily shaped adhesive lap joints

In this work, crack initiation in adhesive lap joints of arbitrary joint configuration is studied by means of a finite fracture mechanics approach. The analysis is based on a general stress solution for adhesive joints combined with a coupled stress and energy criterion. The instantaneous formation of a crack of finite size is predicted if a stress and energy criterion are satisfied simultaneously. The closed-form analytical solution of the stress field allows for an efficient evaluation of the crack initiation load and corresponding finite crack length. A comparison to experimental results from literature and to numerical results obtained with a cohesive zone model approach shows a good agreement. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Numerical modeling of brittle fracture in porous CFC-materials

Both energy and stress criteria are necessary conditions for fracture but neither the one nor the other is sufficient. A combination of these criteria is proposed in [1]. This combined criterion is used for numerical simulation of crack propagation by the 4-point bending test in porous materials. Examples of such materials are carbon-carbon composites (CFC) [2, 3]. Micrographs of the cross-sections of these materials are used for FEM modeling of the crack propagation on the basis of the proposed criterion. Results of the numerical modeling are compared with experimental results. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Analysis of crack formation in bonded lap joints using finite fracture mechanics on the basis of linear elasticity solutions.

In this work, crack formation and the corresponding failure load of bonded lap joints is analyzed. The analysis is based on linear elasticity solutions for bonded lap joints and makes use of the finite fracture mechanics. A hybrid criterion is applied that states the spontaneous formation of a crack of finite size if a stress and an energy criterion are fulfilled simultaneously. The stress distribution of a linear elasticity solution is used for the stress criterion and for the calculation of the incremental energy release rate which is necessary for definition of the energy criterion. The resulting fracture criterion is compared to literature results and shows a good agreement. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

论复合型脆断的周向应变因子准则

本文讨论了线弹性断裂力学的复合型断裂准则.推荐了周向应变因子准则.     

Adaptive mesh generation for curved domains

This paper considers the technologies needed to support the creation of adaptively constructed meshes for general curved three-dimensional domains and outlines one set of solutions for providing them. A brief review of an effective way to integrate mesh generation/adaptation with CAD geometries is given. A set of procedures that support general h-adaptive refinement based on a mesh metric field is given. This is followed by examples that demonstrate the ability of the procedures to adaptively construct anisotropic meshes for flow problems. A procedure for the generation of strongly graded, curved meshes as needed for effective hp-adaptive simulations is also given.     

A computational framework of three dimensional configurational-force-driven crack propagation

A variational formulation of quasi-static brittle fracture is considered and a new finite-element-based computational framework is developed for propagation of cracks in three-dimensional bodies. We outline a consistent thermodynamical framework for crack propagation in elastic solids and show that the crack propagation direction associated with the classical Griffith criterion is identified by the material configurational force which maximizes the local dissipation at the crack front. The evolving crack discontinuity is realized by the doubling of critical nodes and triangular interface facets of the tetrahedral mesh. The crucial step for the success of the procedure is its embedding into an r-adaptive crack-facet reorientation procedure based on configurational-force-based indicators in conjunction with crack front constraints. We further propose a staggered algorithm which minimizes the stored energy at frozen crack state followed by the successive crack releases at frozen deformation. This constitutes a sequence of positive definite subproblems with successively decreasing overall stiffness, providing a very robust algorithmic setting in the postcritical range. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

沿两条内部平行裂纹受压时材料的断裂问题

顺着平面内的平行裂纹受压材料的断裂问题并不能在线性断裂力学的框架内进行描述，Grif-fith－Irvin型或COC型的断裂判据，虽然可以用来处理经典的线性断裂力学，但对本题则完全不适用，这是因为这些压力对应力强度系数没有影响，与裂缝孔隙值也没有关系[1，2]，这一类问题只能采用新的方法，本文的第一作者曾建议过一种新方法，在这一方法中仍然使用了线性关系，但这种线性关系是从变形固体力学中的非线性方程导出的[3，4，5]．这里必须指出，这种方法曾在变形体稳定性问题中广泛地采用过。作为断裂开始的判据，我们采用了裂缝缺陷附近的局部失稳的判据，在这类情况下，我们认为是从失稳过程引发断裂过程的。     

A BE based finite element for crack propagation processes

A boundary element based finite macro element for the simulation of 3D crack propagation in the framework of linear elastic fracture mechanics is presented. While the major part of the numerical model is discretized with finite elements, a small domain containing the crack is meshed with boundary elements. By means of the Symmetric Galerkin BEM a stiffness formulation for the cracked BE domain is obtained which enables a direct FEM/BEM coupling. All necessary operations for the crack propagation are carried out within this boundary element based finite macro element and exploit the potential of the boundary integral formulation. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Formulation and Application of the Initially Rigid Cohesive Finite Element Method

The presented procedure for cohesive crack propagation is based on an adaptive finite element (FE) implementation, which enables the introduction of cohesive surfaces in dependence on the current crack state. In contrast to already existing formulations, the focus of the present model lies on failure processes that can be described at quasi-static conditions within an implicit framework. Furthermore, an extension for mesh independent crack propagation in terms of an additional mesh adaptive formulation is presented. By the evaluation of the failure criterion considering the preferred crack direction, a new crack tip coordinate is computed and the discretization is accordingly adjusted. The remaining mesh is modified for the new boundary representation. The application of the proposed method is shown by the numerical investigation of a concrete fracture specimen from an experimental research project. (© 2009 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Simulation of brittle crack growth in functionally graded materials

We consider a thermodynamic consistent framework for crack propagation by applying a dissipation inequality to a time dependent migrating control volume. The direction of crack growth is obtained in terms of material forces as a result of the principle of maximum dissipation. In the numerical implementation a staggered algorithm – deformation update for fixed geometry followed by geometry update for fixed deformation – is employed within each time increment. The corresponding mesh is generated by combining Delaunay triangulation with local mesh refinement. A numerical example with inhomogeneous material properties illustrates the capability of the resulting algorithm. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

非线性断裂动力学中的路径无关积分和断裂准则

本文考虑非线性断裂动力学中的路径无关积分和断裂准则.在讨论中计入了动力效应和裂纹的传播现象,考虑了裂纹在非线性弹性介质中的传播以及在弹塑性介质中的传播二种情况,作出了一些相应的路径无关积分.作为例子.讨论了裂纹的定常传播情况.最后,给出了这种路径无关积分的力学意义.说明它可用来作为非线性断裂动力学的一种断裂准则.     

Assessment of asymmetric crack initiation in open-hole plates using a coupled stress and energy criterion

An application of the finite fracture mechanics concept to open-hole plates subject to combined tensile and bending loading is presented. In finite fracture mechanics, the simultaneous satisfaction of both, a stress and an energy criterion, is enforced as a condition for crack initiation. Efficient modeling and closed-form expressions for the dependence of the stress and energy quantities on governing structural and material parameters allow for a comprehensive numerical analysis of the onset of asymmetric crack patterns. The obtained failure load predictions are found to agree well with a cohesive zone model and experimental data from literature. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

On an Energetic Interpretation of a Phase Field Model for Fracture

In the pioneering work by Griffith, it is assumed that a crack propagates, if this is energetically favorable. However, this original formulation requires a pre-existing initial crack. In order to bypass this deficiency of classical Griffith theory, Francfort and Marigo advocate a global variational criterion, where the total energy is minimized with respect to any admissible displacement field and crack set. Bourdin's regularized approximation of this variational formulation makes use of a continuous scalar field to indicate cracks. Based on this regularization a phase field fracture model is formulated. The crack field is assumed to follow a Ginzburg-Landau type evolution equation, and cracking is addressed as a phase transition problem. The coupled problem of mechanical balance equations and the evolution equation is solved using the finite element method combined with an implicit time integration scheme. The numerical solution naturally yields the crack evolution including crack propagation, kinking, branching and initiation without any additional criteria. In this work we study the driving mechanisms behind the crack evolution in the phase field fracture model and compare to the purely energetic considerations of the underlying variational formulation. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)