Modeling crack micro-branching using finite elements with embedded strong discontinuities

The emphasis of this work lies in the development of a numerical method which is capable of representing the complex physical phenomena arising in the case of crack branching in brittle materials. In particular, the formation of crack micro-branches needs to be accounted for when it comes to the prediction of the propagation pattern of crack macro-branches which will ultimately lead to the failure of the material. This is achieved by numerically modeling the failure zones within the individual finite elements based on the concept of the embedded finite element method, where all the information with regard to the geometry of the failure zone is stored locally on the element level leading to a very efficient methodology capable of discretely resolving the failure zone. The main feature of the current work is the redundancy of the branching criterion based on crack tip velocity and that both, micro- as well as macro-branches can be modeled. Whether a micro-crack develops into a macro-crack solely depends on the local state of the material as it is outlined based on the application of the proposed numerical scheme on a rectangular block with a pre-existing notch set under tension. A comparison of the oscillatory behavior of the obtained crack tip velocity every time a micro-crack develops with experimental results from the literature is provided. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Modeling quasi-static crack growth with the embedded finite element method on multiple levels

The current work presents the multilevel approach of the embedded finite element method which is obtained by combining features of the method of domain decomposition with those of the standard embedded finite element method. The conventional requirement of fine mesh in a possible failure zone is rendered unnecessary with the new approach thereby reducing the computational expense. In addition, it is also possible to stop a propagating crack-tip in the middle of a finite element. In this approach, the finite elements at the failure-prone zone where cracks or shear bands, referred to as strong discontinuities which represent jumps in the displacement field, can form and propagate based on some failure criterion are treated as separate sub-boundary value problems which are adaptively discretized during the run time into a number of sub-elements and subjected to a kinematic constraint on their boundary. Each sub-element becomes equally capable of developing a strong discontinuity depending upon its state of stress. A linear displacement based constraint is applied initially which is modified accordingly as soon as a strong discontinuity propagates through the boundary of the main finite element. At the local equilibrium, the coupling between the quantities at two different levels of discretization is obtained by matching the virtual energies due to admissible variations of the main finite element and its constituent sub-elements. (© 2012 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

平面弹性裂纹分析的一种有效边界元方法

提出了一种简单而有效的平面弹性裂纹应力强度因子的边界元计算方法．该方法由Crouch与Starfield建立的常位移不连续单元和闫相桥最近提出的裂尖位移不连续单元构成A·D2在该边界元方法的实施过程中,左、右裂尖位移不连续单元分别置于裂纹的左、右裂尖处,而常位移不连续单元则分布于除了裂尖位移不连续单元占据的位置之外的整个裂纹面及其它边界．算例(如单向拉伸无限大板中心裂纹、单向拉伸无限大板中圆孔与裂纹的作用)说明平面弹性裂纹应力强度因子的边界元计算方法是非常有效的．此外,还对双轴载荷作用下有限大板中方孔分支裂纹进行了分析．这一数值结果说明平面弹性裂纹应力强度因子的边界元计算方法对有限体中复杂裂纹的有效性,可以揭示双轴载荷及裂纹体几何对应力强度因子的影响．     

Plastic zone development in dynamic tear-type test specimens

The J-integral is a fracture criterion, which permits measurement of the fracture toughness of a specimen where fracture occurs in the elastic–plastic regime. An understanding of the ratio of plastic zone size (radius) to the crack tip blunting (stretch zone) is required to determine the upper temperature for transition curves, where elastic–plastic fracture becomes invalid and general yielding occurs. This study endeavors to acquire this ratio using finite element techniques. The development of the plastic zone in dynamic tear (DT) specimens and non-standard three-point bending fracture test specimens was the main focus of the study. The ABAQUS finite element software was used to model the elastic–plastic behaviour of the specimens. The cracks in the specimens were induced by pressing the notch followed by fatigue cracking at 30–40% of the limit loads. The shapes of these cracks were adequately modelled in the finite element analysis. The specimens were made of 350WT steel and 304 stainless steel materials and were loaded until fixed amounts of permanent deformation were recorded. Results were obtained in the form of plots, showing the progression of the plastic zone around the crack tip. For each case, mid-point plastic deflection, stretch zone width and plastic zone radius were computed.     

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)     

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.     

A Finite Element Approach for the Simulation of Quasi-Brittle Fracture

In the context of a strong discontinuity approach, we propose a finite element formulation with an embedded displacement discontinuity. The basic assumption of the proposed approach is the additive split of the total displacement field in a continuous and a discontinuous part. An arbitrary crack splits the linear triangular finite element into two parts, namely a triangular and a quadrilateral part. The discontinuous part of the displacement field in the quadrilateral portion is approximated using linear shape functions. For these purposes, the quadrilateral portion is divided into two triangular parts which is in this way similar to the approach proposed in [5]. In contrast, the discretisation is different compared to formulations proposed in [1] and [3], where the discontinuous part of the displacement field is approximated using bilinear shape functions. The basic theory of the underlying finite element formulation and a cohesive interface model to simulate brittle fracture are presented. By means of representative numerical examples differences and similarities of the present formulation and the formulations proposed in [1] and [3] are highlighted. (© 2005 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)     

双I—型裂纹断裂动力学问题的非局部理论解

研究了非局部理论双中I－型裂纹弹性波散射的力学问题,并利用富里叶变换使本问题的求解转换为三重积分方程的求解,进而采用新方法和利用一维非局部积分核代替二维非局部积分核来确定裂纹尖端的应力状态,这种方法就是Schmidt方法,所得结是比艾林根研究断裂静力学问题的结果准确和更加合理,克服了艾林根研究断裂静力学问题时遇到的数学困难,与经典弹性解相比,裂纹尖端不再出现物理意义下不合理的应力奇异性,并能够解释宏观裂纹与微观裂纹的力学问题。     

Evaluation of mixed mode stress intensity factors for interface cracks using EFGM

This paper presents the implementation of element free Galerkin method for the stress analysis of structures having cracks at the interface of two dissimilar materials. The material discontinuity at the interface has been modeled using a jump function with a jump parameter that governs its strength. The jump function enriches the approximation by the addition of special shape function that contains discontinuities in the derivative. The trial and test functions of the weak form are constructed using moving least-square interpolants in each material domain. An intrinsic enrichment criterion with enriched basis has been used to model the crack tip stress fields. The mixed mode (complex) stress intensity factors for bi-material interface cracks are numerically evaluated using the modified domain form of interaction integral. The numerical results are obtained for edge and center cracks lying at the bi-material interface, and are found to be in good agreement with the reference solutions for the interfacial crack problems.     

Multiple crack fatigue growth modeling by displacement discontinuity method with crack-tip elements

This paper presents a numerical approach for modeling multiple crack fatigue growth in a plane elastic infinite plate. It involves a generation of Bueckner’s principle, a displacement discontinuity method with crack-tip elements (a boundary element method) proposed recently by the author and an extension of Paris’ law to a multiple crack problem under mixed-mode loading. Because of an intrinsic feature of the boundary element method, a general multiple crack growth problem can be solved in a single-region formulation. In the numerical simulation, for each increment of crack extension, remeshing of existing boundaries is not necessary. Crack extension is conveniently modeled by adding new boundary elements on the incremental crack extension to the previous crack boundaries. Fatigue growth modeling of an inclined crack in an infinite plate under biaxial cyclic loads is taken into account to illustrate the effectiveness of the present numerical approach. As an example, the present numerical approach is used to study the fatigue growth of three parallel cracks with same length under uniaxial cyclic load. Many numerical results are given.     

基于比例边界有限元法和灰狼优化算法的裂纹尖端位置识别

基于比例边界有限元法(SBFEM)和灰狼优化(GWO)算法,提出了一种裂纹尖端识别方法。首先,借助SBFEM解决断裂力学问题特有的优势,快速准确地计算出反演所需的测点位移,并验证了正问题求解的正确性。其次,建立与裂纹尖端位置有关的目标函数,将求解裂纹尖端位置转换为求解目标函数最小值的优化问题。最后,采用GWO算法对目标函数进行了优化,进而搜索裂纹尖端的最佳位置。数值算例结果表明:利用SBFEM的高精度、半解析的优点,在反演过程中采用其求解正问题是非常有效的;GWO算法具有良好的全局收敛性,且相比经典的粒子群算法,能够更快速准确地搜索出裂纹尖端的位置;GWO算法具有较好的抗噪性。     

A boundary element analysis for stress intensity factors of multiple circular arc cracks in a plane elasticity plate

In this paper, a numerical approach for analyzing interacting multiple cracks in infinite linear elastic media is presented. By extending Bueckner’s principle suited for a crack to a general system containing multiple interacting cracks, the original problem is divided into a homogeneous problem (the one without cracks) subjected to remote loads and a multiple crack problem in an unloaded body with applied tractions on the crack surfaces. Thus, the results in terms of the stress intensity factors (SIFs) can be obtained by considering the latter problem, which is analyzed easily by means of the displacement discontinuity method with crack-tip elements proposed recently by the author. Test examples are given to illustrate that the numerical approach is very accurate for analyzing interacting multiple cracks in an infinite linear elastic media under remote uniform stresses. In addition, the displacement discontinuity method with crack-tip elements is used to analyze a multiple crack problem in a finite plate. It is found that the boundary element method is also very accurate for investigating interacting multiple cracks in a finite plate. Specially, a generalization of Bueckner’s principle and the displacement discontinuity method with crack-tip elements are used to analyze multiple circular arc crack problems in infinite plate in tension (including: Two Collinear Circular Arc Cracks, Three Collinear Circular Arc Cracks, Two Parallel Circular Arc Cracks, Three Parallel Circular Arc Cracks and Two Circular Arc Cracks) in a plane elasticity plate. Many results are given.     

A Phase Field Approach for Dynamic Fracture

In phase field fracture models cracks are indicated by the value of a scalar field variable which interpolates smoothly between broken and undamaged material. The evolution equation for this crack field is coupled to the mechanical field equations in order to model the mutual interaction between the crack evolution and mechanical quantities. In finite element simulations of crack growth at comparatively slow loading velocities, a quasi-static phase field model yields reasonable results. However, the simulation of fast loading or the nucleation of new cracks challenges the limits of such a formulation. Here, the quasi-static phase field model predicts brutal crack extension with an artificially high crack speed. In this work, we analyze to which extend a dynamic formulation of the mechanical part of the phase field model can overcome this paradox created by the quasi-static formulation. In finite element simulations, the impact of the dynamic effects is studied, and differences between the crack propagation behavior of the quasi-static model and the dynamic formulation are highlighted. (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Hexahedral finite element with an open crack

A method of creating the stiffness matrix of a hexahedral eight-node finite element with a single, nonpropagating, transverse, one-edge crack at half of its length is presented in this paper. The crack was modelled by adding an additional flexibility matrix to that of the noncracked element. The terms of the additional matrix have been calculated by use of the laws of fracture mechanics. Employing the elaborated element a numerical test has been worked out, the results of which are compared with the data of analytical solutions accessible in the literature, and a high conformity with them has been obtained. The element presented in the paper may be applied to the static and dynamic analysis of different types of structural elements with material defects in the form of cracks. The described method of creating the stiffness matrix of the element allows to create different kinds of finite elements with cracks provided that the stress intensity factors for a given type of crack are known.     

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)     

Subcritical and critical states of a crack with failure zones

The problem on the stress–strain state near a mode I crack in an infinite plate is solved in the frame of a cohesive zone model. The complex variable method of Muskhelishvili is used to obtain the crack opening displacements caused by the cohesive traction, which models the failure zone at the crack tip, as well as by the external load. The finite stress condition and logarithmic singularity of the derivative of the separation with respect to the coordinate at the tip of a physical crack are taken into account.The cohesive traction distribution is sought in a piecewise linear form, nodal values of which are being numerically chosen to satisfy the traction-separation law. According to this law, the cohesive traction is coupled with the corresponding separation and fracture toughness. The tips of the physical crack and cohesive zone (geometric variables) along with the discrete cohesive traction are used as the problem parameters determining the stress-strain state. If the crack length is included in the set, then the critical crack size can be found for the given loading intensity.The obtained determining system of equations is solved numerically. To find the initial point for a standard numerical algorithm, the asymptotic determining system is derived. In this system, the geometric variables can be easily eliminated, which make it possible to linearize the system.In the numerical examples, the one-parameter traction-separation laws are used. Influence of the shape parameters of the law on the critical crack size and the corresponding cohesive length is studied. The possibility of using asymptotic solutions for determining the critical parameters is analysed. It is established that the critical crack length slightly depends on the shape parameter, while the cohesive length shows a strong dependence on the shape of cohesive laws.     

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)     

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)     

Similarity and dimensional analysis in the plane problem of the propagation of a vertical hydraulic fracture crack in an elastic medium

The problem of the growth of a vertical hydraulic fracture crack in an unbounded elastic medium under the pressure produced by a viscous incompressible fluid is studied qualitatively and by numerical methods. The fluid motion is described in the approximation of lubrication theory. Near the crack tip a fluid-free domain may exist. To find the crack length, Irwin’s fracture criterion is used. The symmetry groups of the equations describing the hydraulic fracture process are studied for all physically meaningful cases of the degeneration of the problem with respect to the control parameters. The condition of symmetry of the system of equations under the group of scaling and time-shift transformations enables the self-similar variables and the form of the time dependence of the quantities involved in the problem to be found. It is established that at non-zero rock pressure the well-known solution of Spence and Sharp is an asymptotic form of the initial-value problem, whereas the solution of Zheltov and Khristianovich is a limiting self-similar solution of the problem. The problem of the formation of a hydraulic fracture crack taking into account initial data is solved using numerical methods, and the problem of arriving at asymptotic mode is investigated. It is shown that the solution has a self-similar asymptotic form for any initial conditions, and the convergence of the exact solutions to the asymptotic forms is non-uniform in space and time.