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.     

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)     

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)     

Strategies for the Computation of Configurational Forces in Dissipative Media

Configurational forces can be interpreted as driving forces on material inhomogeneities such as crack tips. In dissipative media the total configurational force on an inhomogeneity consists of an elastic contribution and a contribution due to the dissipative processes in the material. For the computation of discrete configurational forces acting at the nodes of a finite element mesh, the elastic and dissipative contributions must be evaluated at integration point level. While the evaluation of the elastic contribution is straightforward, the evaluation of the dissipative part is faced with certain difficulties. This is because gradients of internal variables are necessary in order to compute the dissipative part of the configurational force. For the sake of efficiency, these internal variables are usually treated as local history data at integration point level in finite element (FE) implementations. Thus, the history data needs to be projected to the nodes of the FE mesh in order to compute the gradients by means of shape function interpolations of nodal data as it is standard practice. However, this is a rather cumbersome method which does not easily integrate into standard finite element frameworks. An alternative approach which facilitates the computation of gradients of local history data is investigated in this work. This approach is based on the definition of subelements within the elements of the FE mesh and allows for a straightforward integration of the configurational force computation into standard finite element software. The suitability and the numerical accuracy of different projection approaches and the subelement technique are discussed and analyzed exemplarily within the context of a crystal plasticity model. (© 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

On the peculiarities of a crack growth path in anisotropic solid

One of the possibilities to increase the resistance of a structure to catastrophic fracture is to force a main line crack to deviate from its path. In this connection the influence of the elastic moduli of an anisotropic material on the possibilities of crack rotation are studied. In particular a linear elastic problem for a straight Mode I crack, located on a symmetry axis of an orthotropic plane is considered. The strength properties of the material are supposed to be isotropic. For studying a direction of a crack growth path several crack models are considered. It is shown that a thin elongated elliptical hole as a crack model leads to more plausible results concerning crack rotation conditions than an ideal cut model. The maximal tensile stresses are taken as a crack growth criterion. It is shown that for some class of orthotropic materials a crack deviates from the straight path just after it starts to grow even in the conditions of uniaxial normal tension. The problem of the stability of a straight crack path under Mode I loading is also considered. This problem is reduced to the problem of the fracture direction determination for thin elongated elliptical cavity slightly inclined to the initial direction. In the frame of the proposed approach the conditions of instability are obtained. It is shown that for some class of orthotropic materials a straight crack path is unstable in the conditions of uniaxial normal tension. This class of materials is wider than one for which a crack deviates from the straight crack path just after its start. (© 2008 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)     

Energy driven crack propagation at finite strains based on the embedded strong discontinuity approach

New advances in three-dimensional finite element modeling of crack propagation at finite strains are presented. The proposed numerical model is based on the Enhanced Assumed Strain concept. The enhanced part of the deformation gradient is associated with a displacement discontinuity. In contrast to previous works, a new, energy based criterion for crack propagation is presented. The necessity for a tracking algorithm for the crack path is avoided by using more than one discontinuity within each finite element. This leads to a strictly local formulation, i.e., no information about the neighboring elements are required. Further advantages of such a formulation are a symmetric tangent stiffness matrix and the reduction of locking effects. (© 2009 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)     

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)     

Modelling of metal forging using SPH

In solid metal forming processes, such as forging, large distortions in the material present challenging problems for numerical simulation using grid based methods. Computations invariably fail after some level of mesh distortion is reached unless suitable re-meshing is implemented to cope with the mesh distortion arising from the material deformation. The issue of mesh distortion and the subsequent re-meshing are topics of much research for grid based methods. These problems can be overcome by using a mesh-less numerical framework. In this paper, the application of a mesh-less method called Smoothed Particle Hydrodynamics (SPH) for modelling three-dimensional complex forging processes is demonstrated. It is shown that SPH is a useful simulation method for obtaining insights into the material deformation and flow pattern during forging of realistic industrial components. The effect of process parameters and material properties on the quality of the forged component is evaluated via SPH simulations. This includes the determination of forging force required for adequate die filling which is an important criterion for die designs. Material hardening, controlled by the degree of heat treatment, is found to have a profound effect on the material deformation pattern and the final product. Forging defects such as incomplete die filling, asymmetry in forged components, flashing and lap formation are shown to be predicted by SPH. SPH can thus potentially be used both for assessment of the quality of forged products and evaluation of prototype forging system designs.     

Crack propagation calculation using trial cracks

We present a method to perform finite element calculations for crack propagation problems with arbitrary crack directions in two dimensions. The crack direction (angle of propagation) is determined by inserting small “trial cracks” at the crack tip. For each trial crack, the domain is remeshed to allow crack propagation between elements. The trial cracks are then opened and the energy release rate is measured. The optimum crack direction (i.e., the crack direction with maximum energy release) is determined by an optimisation procedure. Although the method is computationally expensive due to the need to perform several calculations for each crack increment, it has the advantage that the energy release rate can be calculated even in cases where other methods fail. After explaining the method, it is applied to some test examples. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

圆形杂质对裂纹扩展的影响

在单轴拉伸载荷作用下，运用分布位错方法对无限大平面内含有一个裂纹和一个任意方向的杂质问题进行求解，得到了裂纹尖端的应力强度因子、应力场以及应变能密度.利用最小应变能密度因子准则来判断裂纹扩展方向.结果显示:软杂质对裂纹尖端应力强度因子、应变能密度和应力场有增强作用，而硬杂质则具有屏蔽作用.在 -30°<θ<30°范围内，杂质对裂纹扩展方向的影响较小，而在 -90°<θ<-30°或30°<θ<90°范围内，杂质对裂纹扩展方向的影响较大.软杂质对裂纹扩展有吸引作用，而硬杂质具有排斥作用.     

Configurational-force-based structural updates and adaptive remeshing strategies

The application of configurational forces in h -adaptive strategies for fracture mechanics and inelasticity is investigated. Starting from a global Clausius-Planck inequality, dual equilibrium conditions are derived by means of a Coleman-type exploitation method. The remaining reduced dissipation inequality is used for the derivation of evolution equations for the internal variables. In fracture mechanics, crack loading conditions as well as a normality rule for the crack propagation are obtained. In the discrete setting, the crack propagation is governed by a configurational-force-driven update of the geometry model. The material balance equation is used to set up a h -adaptive refinement indicator. A relative global criterion is defined used for the decision on mesh refinement. In addition, a criterion on the element level is evaluated controlling the local refinement procedure. The capability of the proposed procedures is demonstrated by means of numerical examples. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Comparison of criteria on the direction of crack extension

There are three famous criteria for the prediction on the direction of crack extension: maximum energy release rate criterion, local symmetric criterion and maximum stress criterion. It has already been proved that these criteria are different. In this paper, we show how to measure the difference between maximum energy release rate criterion and maximum stress criterion. The results in this paper and special solutions in Amestoy and Leblond (Internat. J. Solids Struct. 29 (1992) 465) indicate that it is enough to use the maximum stress criterion for the numerical prediction of crack path, if the stress intensity factor of the shearing mode is not too big.     

A Variational–Based Formulation of Regularized Brittle Fracture

The paper performs a comparative study of variational-based brittle fracture with its gradient-type regularization, and outlines aspects of the numerical implementations of both approaches. The latter smoothes out sharp displacement discontinuities of cracks. On the side of discrete crack modeling, we propose a variational framework of configurational-force-driven crack propagation. The latter provides the basis for the computation of material nodal forces and drives the crack propagation in our proposed finite element framework with adaptive nodal doubling. Such a formulation is of limited applicability for the modeling of crack inititation in homogeneous bodies without defects and in situations with complex crack branching. This can be overcome by a regularized crack modeling. Here, an elliptic approximation of the crack surface term yields a regularized two field functional, where an additional scalar field approximates the set of discontinuities. We provide a simple interpretation of such a transition from the sharp crack to the regularized setting. It results in a smooth continuum-damage-type theory with a specific gradient-damage and hardening terms, depending on a length scale that represents the width of a zone that surrounds the crack. Such a variational framework is implemented by a coupled two-field finite element framework in a staightforward manner. We compare representative numerical results obtained by both methods for selected crack patterns and highlight the pro and contra of both meshes. (© 2009 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Crack – interface crack interactions in functionally graded/homogeneous composite bimaterials subjected to a heat flux

A bimaterial containing an interface crack and consisting of a homogeneous material and a functionally graded material (FGM) with a system of small internal cracks is considered. The thermal fracture of the biomaterial under the action of a heat flux applied to it at infinity is investigated. The problem is studied in the case where the interface crack is much larger than the internal ones. It is assumed that the thermal properties of the FGM are continuous functions of the thickness coordinate. Asymptotic analytical formulas for the thermal stress intensity factors (TSIFs) at the tips of the interface crack are obtained as series in a small parameter (the ratio between sizes of the internal and interface cracks). Then, the critical heat fluxes and the angles of propagation direction of the interface crack are calculated using the criterion of maximum circumferential stress. A parametric analysis shows that the propagation direction of the interface crack depends on the location and orientation of the system of internal cracks. The parameters of inhomogeneity of the FGM affect the value of TSIFs and, hence, the deflection angle of the interface crack.     

Thermo-mechanical modeling of crack propagation in dynamically loaded elastomer specimens using a scaled boundary finite element approach

Recently, a scaled boundary finite element (SBFE) formulation for geometrically and physically nonlinear materials has been developed using the scaled boundary finite element method (SBFEM). The SBFE formulation has been employed to describe plane stress problems of notched and unnotched hyperelastic elastomer specimens. In this contribution, the derived SBFE formulation is extended to nonlinear time- and temperature-dependent material behavior. Subsequently, the SBFE formulation is incorporated into a crack propagation scheme to model crack propagation in cyclically loaded elastomer specimens of the so-called tear fatigue analyzer (TFA). (© 2015 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

An accurate and efficient algorithm for the simulation of fatigue crack growth based on XFEM and combined approximations

This study aims to reduce the computational cost in the simulation of fatigue crack growth process. Extended finite element method and combined approximations (CA) are integrated to form an efficient algorithm for such analysis. In the CA approach, binomial series are used as high quality basis vectors for the reduced basis expression and then crack propagation path can be predicted without solving the complete set of system equations. The validity of presented algorithm is fully investigated through several numerical examples. From these results, it is shown that the presented algorithm is very accurate and can save huge amounts of computational effort.     

A generalized cohesive element technique for arbitrary crack motion

A computational method for arbitrary crack motion through a finite element mesh, termed as the generalized cohesive element technique, is presented. In this method, an element with an internal discontinuity is replaced by two superimposed elements with a combination of original and imaginary nodes. Conventional cohesive zone modeling, limited to crack propagation along the edges of the elements, is extended to incorporate the intra-element mixed-mode crack propagation. Proposed numerical technique has been shown to be quite accurate, robust and mesh insensitive provided the cohesive zone ahead of the crack tip is resolved adequately. A series of numerical examples is presented to demonstrate the validity and applicability of the proposed method.