A framework for fracture modelling based on the material forces concept with XFEM kinematics

A theoretical and computational framework which covers both linear and non‐linear fracture behaviour is presented. As a basis for the formulation, we use the material forces concept due to the close relation between on one hand the Eshelby energy–momentum tensor and on the other hand material defects like cracks and material inhomogeneities. By separating the discontinuous displacement from the continuous counterpart in line with the eXtended finite element method (XFEM), we are able to formulate the weak equilibrium in two coupled problems representing the total deformation. However, in contrast to standard XFEM, where the direct motion discontinuity is used to model the crack, we rather formulate an inverse motion discontinuity to model crack development. The resulting formulation thus couples the continuous direct motion to the inverse discontinuous motion, which may be used to simulate linear as well as non‐linear fracture in one and the same formulation. In fact, the linear fracture formulation can be retrieved from the non‐linear cohesive zone formulation simply by confining the cohesive zone to the crack tip. These features are clarified in the two numerical examples which conclude the paper. Copyright © 2005 John Wiley & Sons, Ltd.     

A computational framework of configurational-force-driven brittle fracture based on incremental energy minimization

A variational formulation of quasi-static brittle fracture in elastic solids at small strains is proposed and an associated finite element implementation is presented. On the theoretical side, a consistent thermodynamic framework for brittle crack propagation is outlined. It is shown that both the elastic equilibrium response as well as the local crack evolution follow in a natural format by exploitation of a global Clausius–Planck inequality. Here, the canonical direction of the crack propagation associated with the classical Griffith criterion is the direction of the material configurational force which maximizes the local dissipation at the crack tip. On the numerical side, we first consider a standard finite element discretization in the two-dimensional space which yields a discrete formulation of the global dissipation in terms of configurational nodal forces. Next, consistent with the node-based setting, the discretization of the evolving crack discontinuity for two-dimensional problems is performed by the doubling of critical nodes and interface segments of the mesh. A crucial step for the success of this procedure is its embedding into a r-adaptive crack-segment re-orientation algorithm governed by configurational-force-based directional indicators. Here, successive crack propagation is performed by a staggered loading-release algorithm of energy minimization at frozen crack state followed by nodal releases at frozen deformation. We compare results obtained by the proposed formulation with other crack propagation criteria. The computational method proposed is extremely robust and shows an excellent performance for representative numerical simulations.     

An extended finite element method with analytical enrichment for cohesive crack modeling

A recent approach to fracture modeling has combined the extended finite element method (XFEM) with cohesive zone models. Most studies have used simplified enrichment functions to represent the strong discontinuity but have lacked an analytical basis to represent the displacement gradients in the vicinity of the cohesive crack. In this study enrichment functions based upon an existing analytical investigation of the cohesive crack problem are proposed. These functions have the potential of representing displacement gradients in the vicinity of the cohesive crack and allow the crack to incrementally advance across each element. Key aspects of the corresponding numerical formulation and enrichment functions are discussed. A parameter study for a simple mode I model problem is presented to evaluate if quasi‐static crack propagation can be accurately followed with the proposed formulation. The effects of mesh refinement and mesh orientation are considered. Propagation of the cohesive zone tip and crack tip, time variation of the cohesive zone length, and crack profiles are examined. The analysis results indicate that the analytically based enrichment functions can accurately track the cohesive crack propagation of a mode I crack independent of mesh orientation. A mixed mode example further demonstrates the potential of the formulation. Copyright © 2008 John Wiley & Sons, Ltd.     

An implicit adaptive node‐splitting algorithm to assess the failure mechanism of inelastic elastomeric continua

This contribution presents a mesh adaptive crack propagation scheme for the evaluation of the viscoelastic fracture response of elastomers at large strains and up to high loading rates. The approach accounts for micromechanical based features of both elastic and viscoelastic bulk responses of idealized polymer networks. To this end, the Bergstörm–Boyce model is considered to introduce hyperelastic and nonlinear finite viscoelastic responses. Moreover, the crack driving force and the crack driving direction are predicted by the material force approach. A consistent thermodynamic framework for the combined configurational motion in viscoelastic continua at finite strain regime is discussed. The fracture toughness of non‐strain‐crystallizing elastomers shows strong rate dependency and the energy release rate versus the rate of tearing to be a fundamental material property. Therefore, in this contribution, a dynamic fracture criterion, which is a function of the rate of crack growth, is shown to be adequate in numerical simulations. The use of the presented method enables to study fracture behaviour of any material nonlinearity within the implicit time integration. Main feature of the proposed algorithm is restructuring the overall discrete system by duplication of crack front DOFs based on minimization of the overall energy via the Griffith criterion. Copyright © 2014 John Wiley & Sons, Ltd.     

Bounds for element size in a variable stiffness cohesive finite element model

The cohesive finite element method (CFEM) allows explicit modelling of fracture processes. One form of CFEM models integrates cohesive surfaces along all finite element boundaries, facilitating the explicit resolution of arbitrary fracture paths and fracture patterns. This framework also permits explicit account of arbitrary microstructures with multiple length scales, allowing the effects of material heterogeneity, phase morphology, phase size and phase distribution to be quantified. However, use of this form of CFEM with cohesive traction–separation laws with finite initial stiffness imposes two competing requirements on the finite element size. On one hand, an upper bound is needed to ensure that fields within crack‐tip cohesive zones are accurately described. On the other hand, a lower bound is also required to ensure that the discrete model closely approximates the physical problem at hand. Both issues are analysed in this paper within the context of fracture in multi‐phase composite microstructures and a variable stiffness bilinear cohesive model. The resulting criterion for solution convergence is given for meshes with uniform, cross‐triangle elements. A series of calculations is carried out to illustrate the issues discussed and to verify the criterion given. These simulations concern dynamic crack growth in an Al₂O₃ ceramic and in an Al₂O₃/TiB₂ ceramic composite whose phases are modelled as being hyperelastic in constitutive behaviour. Copyright © 2004 John Wiley & Sons, Ltd.     

A finite element method for the computational modelling of cohesive cracks

The present contribution is concerned with the computational modelling of cohesive cracks in quasi‐brittle materials, whereby the discontinuity is not limited to interelement boundaries, but is allowed to propagate freely through the elements. In the elements, which are intersected by the discontinuity, additional displacement degrees of freedom are introduced at the existing nodes. Therefore, two independent copies of the standard basis functions are used. One set is put to zero on one side of the discontinuity, while it takes its usual values on the opposite side, and vice versa for the other set. To model inelastic material behaviour, a discrete damage‐type constitutive model is applied, formulated in terms of displacements and tractions at the surface. Some details on the numerical implementation are given, concerning the failure criterion, the determination of the direction of the discontinuity and the integration scheme. Finally, numerical examples show the performance of the method. Copyright © 2005 John Wiley & Sons, Ltd.     

A refined diffuse cohesive approach for the failure analysis in quasibrittle materials—part I: Theoretical formulation and numerical calibration

In this work, a refined interelement diffuse fracture theoretical model, based on a cohesive finite element approach, is proposed for concrete and other quasibrittle materials. This model takes advantage of a novel micromechanics‐based calibration technique for reducing the artificial compliance associated with the adopted intrinsic formulation. By means of this technique, the required values for the elastic stiffness parameters to obtain nearly invisible cohesive interfaces are provided. Furthermore, the mesh‐induced toughening effect, essentially related to the artificial crack tortuosity caused by the different orientations of the interelement cohesive interfaces, is numerically investigated by performing comparisons with an additional fracture model, newly introduced for the purpose of numerical validation. These comparisons are presented to assess the reliability and the numerical accuracy of the proposed fracture approach.     

A material model to reproduce mixed‐mode fracture in concrete

This paper presents a material model to reproduce crack propagation in cement‐based material specimens under mixed‐mode loading. Its numerical formulation is based on the cohesive crack model, proposed by Hillerborg, and extended for the mixed‐mode case. This model is inspired by former works by Gálvez et al but implemented for its use in a finite element code at a material level, that is to say, at an integration point level. Among its main features, the model is able to predict the crack orientation and can reproduce the fracture behaviour under mixed‐mode fracture loading. In addition, several experimental results found in the literature are properly reproduced by the model.     

Cohesive crack propagation in a random elastic medium

The issue of generating non-Gaussian, multivariate and correlated random fields, while preserving the internal auto-correlation structure of each single-parameter field, is discussed with reference to the problem of cohesive crack propagation. Three different fields are introduced to model the spatial variability of the Young modulus, the tensile strength of the material, and the fracture energy, respectively. Within a finite-element context, the crack-propagation phenomenon is analyzed by coupling a Monte Carlo simulation scheme with an iterative solution algorithm based on a truly-mixed variational formulation which is derived from the Hellinger–Reissner principle. The selected approach presents the advantage of exploiting the finite-element technology without the need to introduce additional modes to model the displacement discontinuity along the crack boundaries. Furthermore, the accuracy of the stress estimate pursued by the truly-mixed approach is highly desirable, the direction of crack propagation being determined on the basis of the principal-stress criterion. The numerical example of a plain concrete beam with initial crack under a three-point bending test is considered. The statistics of the response is analyzed in terms of peak load and load–mid-deflection curves, in order to investigate the effects of the uncertainties on both the carrying capacity and the post-peak behaviour. A sensitivity analysis is preliminarily performed and its results emphasize the negative effects of not accounting for the auto-correlation structure of each random field. A probabilistic method is then applied to enforce the auto-correlation without significantly altering the target marginal distributions. The novelty of the proposed approach with respect to other methods found in the literature consists of not requiring the a priori knowledge of the global correlation structure of the multivariate random field.     

Dynamic response of interfacial finite cracks in orthotropic materials subjected to concentrated loads

A cohesive zone model has been proposed to model crack growth with a part-through process zone in a thin solid. With the solid being modeled in Kirchhoff’s plate theory, the crack with a relatively long, inclined front is modeled as a line discontinuity with a finite cohesive zone within the plate. A cohesive force law is adopted to capture the effect of residual strength and residual rigidity of a plate cross-section gradually cracking through the thickness. It is derived by a plane-strain elasticity analysis of a cross section normal to the part-through crack. It is then applied in the plate formulation of a line crack to simulate its propagation within the plate plane. This model essentially resolves the originally three-dimensional crack problem in two hierarchical steps, i.e., in the thickness and in the in-plane directions. In the present study, the bending case is considered. A boundary element method is applied to numerically derive the cohesive force law and simulate the crack growth in a thin titanium-alloy plate. The computational efficiency of the model is demonstrated. The plate is shown to fracture in a nominally brittle or ductile manner depending on its thickness.     

An embedded cohesive crack model for finite element analysis of mixed mode fracture of concrete*

An embedded cohesive crack model is proposed for the analysis of the mixed mode fracture of concrete in the framework of the Finite Element Method. Different models, based on the strong discontinuity approach, have been proposed in the last decade to simulate the fracture of concrete and other quasi‐brittle materials. This paper presents a simple embedded crack model based on the cohesive crack approach. The predominant local mode I crack growth of the cohesive materials is utilized and the cohesive softening curve (stress vs. crack opening) is implemented by means of a central force traction vector. The model only requires the elastic constants and the mode I softening curve. The need for a tracking algorithm is avoided using a consistent procedure for the selection of the separated nodes. Numerical simulations of well‐known experiments are presented to show the ability of the proposed model to simulate the mixed mode fracture of concrete.     

Extended Voronoi cell finite element model for multiple cohesive crack propagation in brittle materials

This paper introduces an extended Voronoi cell finite‐element model (X‐VCFEM) for modelling cohesive crack propagation in brittle materials with multiple cracks. The cracks are modelled by a cohesive zone model and their incremental directions and growth lengths are determined in terms of the cohesive energy near the crack tip. Extension to VCFEM is achieved through enhancements in stress functions in the assumed stress hybrid formulation. In addition to polynomial terms, the stress functions include branch functions in conjunction with level set methods, and multi‐resolution wavelet functions in the vicinity of crack tips. The wavelet basis functions are adaptively enriched to accurately capture crack‐tip stress concentrations. Conditions and methods of stability are enforced in X‐VCFEM for improved convergence with propagating cracks. Two classes of problems are solved and compared with existing solutions in the literature for validation of the X‐VCFEM algorithms. The first set corresponds to results for static cracks, while in the latter set, the propagation of cohesive cracks are considered. Comparison of X‐VCFEM simulation results with results in literature for several fracture mechanics problems validates the effectiveness of X‐VCFEM. Copyright © 2005 John Wiley & Sons, Ltd.     

A cohesive segments method for the simulation of crack growth

 A numerical method for crack growth is described in which the crack is not regarded as a single discontinuity that propagates continuously. Instead, the crack is represented by a set of overlapping cohesive segments. These cohesive segments are inserted into finite elements as discontinuities in the displacement field by exploiting the partition-of-unity property of shape functions. The cohesive segments can be incorporated at arbitrary locations and orientations and are not tied to any particular mesh direction. The evolution of decohesion of the segments is governed by a cohesive law. The independent specification of bulk and cohesive constitutive relations leads to a characteristic length being introduced into the formulation. The formulation permits both crack nucleation and discontinuous crack growth to be modelled. The implementation is outlined and some numerical examples are presented. Dedicated to the memory of Prof. Mike Crisfield, for his cheerfulness and cooperation as a colleague and friend over many years. The authors wish to express their thanks to Erik-Jan Lingen for his help in the implementation of the model in the JIVE finite element toolbox. AN is grateful for support from the Office of Naval Research through grant N00014-97-1-0179.     

黏聚裂纹阻抗的弯曲梁承载力

在混凝土类软化材料断裂研究中,裂纹端部损伤区被简化为具有黏聚应力分布的非线性裂纹,该黏 聚力对裂纹扩展有阻抗作用。裂纹体的应力强度因子是断裂力学标志载荷作用与几何构型因素的量化表达指标; 黏聚力形成的阻抗强度因子数值,与黏聚裂纹长度和材料极值拉伸应力存在数量关系。通过双K断裂判据,以 带切口的三点弯曲梁为断裂力学模型,分析了裂纹黏聚阻力对断裂过程的影响规律,计算该弯曲梁结构断裂试 样的最大承担载荷;其结果比较符合实验数据。     

Asymptotic stress intensity factor density profiles for smeared-tip method for cohesive fracture

The paper presents a computational approach and numerical data which facilitate the use of the smeared-tip method for cohesive fracture in large enough structures. In the recently developed K-version of the smeared tip method, the large-size asymptotic profile of the stress intensity factor density along a cohesive crack is considered as a material characteristic, which is uniquely related to the softening stress-displacement law of the cohesive crack. After reviewing the K-version, an accurate and efficient numerical algorithm for the computation of this asymptotic profile is presented. The algorithm is based on solving a singular Abel's integral equation. The profiles corresponding to various typical softening stress-displacement laws of the cohesive crack model are computed, tabulated and plotted. The profiles for a certain range of other typical softening laws can be approximately obtained by interpolation from the tables. Knowing the profile, one can obtain with the smeared-tip method an analytical expression for the large-size solution to fracture problems, including the first two asymptotic terms of the size effect law. Consequently, numerical solutions of the integral equations of the cohesive crack model as well as finite element simulations of the cohesive crack are made superfluous. However, when the fracture process zone is attached to a notch or to the body surface and the cohesive zone ends with a stress jump, the solution is expected to be accurate only for large-enough structures.     

3D modelling of strong discontinuities in elastoplastic solids: fixed and rotating localization formulations

This paper is concerned with the incorporation of displacement discontinuities into a continuum theory of elastoplasticity for the modelling of localization processes such as cracking in brittle materials. Based on the strong discontinuity approach (SDA) (Computational Mechanics 1993; 12: 277–296) mesh objective 2D and 3D finite element formulations are developed using linear and quadratic 2D elements as well as 8‐noded 3D elements. In the formulation of the finite‐element model proposed in the paper, the analogy with standard formulations is emphasized. The parameter defining the amplitude of the displacement jump within the finite element is condensed out at the material level without employing the standard static condensation technique. This approach results in linearized constitutive equations formally identical to continuum models. Therefore, the standard return mapping algorithm is used to solve the non‐linear equations. In analogy to concepts used in continuum smeared crack models, a rotating formulation of the SDA is proposed in addition to the standard concept of fixed discontinuities. It is shown that the rotating localization approach reduces locking effects observed in analyses based on fixed localization directions. The applicability of the proposed SDA finite‐element model as well as its numerical performance is investigated by means of a three‐dimensional ultimate load analysis of a steel anchor embedded in a concrete block subjected to a shear force. Copyright © 2003 John Wiley & Sons, Ltd.     

Cohesive surface model for fracture based on a two-scale formulation: computational implementation aspects

The paper describes the computational aspects and numerical implementation of a two-scale cohesive surface methodology developed for analyzing fracture in heterogeneous materials with complex micro-structures. This approach can be categorized as a semi-concurrent model using the representative volume element concept. A variational multi-scale formulation of the methodology has been previously presented by the authors. Subsequently, the formulation has been generalized and improved in two aspects: (i) cohesive surfaces have been introduced at both scales of analysis, they are modeled with a strong discontinuity kinematics (new equations describing the insertion of the macro-scale strains, into the micro-scale and the posterior homogenization procedure have been considered); (ii) the computational procedure and numerical implementation have been adapted for this formulation. The first point has been presented elsewhere, and it is summarized here. Instead, the main objective of this paper is to address a rather detailed presentation of the second point. Finite element techniques for modeling cohesive surfaces at both scales of analysis (FE\(^2\) approach) are described: (i) finite elements with embedded strong discontinuities are used for the macro-scale simulation, and (ii) continuum-type finite elements with high aspect ratios, mimicking cohesive surfaces, are adopted for simulating the failure mechanisms at the micro-scale. The methodology is validated through numerical simulation of a quasi-brittle concrete fracture problem. The proposed multi-scale model is capable of unveiling the mechanisms that lead from the material degradation phenomenon at the meso-structural level to the activation and propagation of cohesive surfaces at the structural scale.     

Combination of the material force concept and the extended finite element method for mixed mode crack growth simulations

A novel approach to simulate crack growth within an extended finite element framework is presented. The introduced approach combines the material force concept and the extended finite element method (xFEM) that is not straight forward and faces the major problem that a crack tip node, which is required for the evaluation of the material force, is not available within an xFEM framework. The introduced concept enables an efficient single step evaluation of the crack state and the crack growth direction based on a continuum mechanics approach and represents an alternative to the common procedure of using the stress intensity factor solution within a stress or energy‐based empirical formulation for the determination of the crack growth direction. Two different approaches are introduced that evaluate the crack tip material force within the xFEM based on a domain or contour approach, both providing equivalent results. After an evaluation of the method, a major focus is set on crack growth investigations with increased complexity, including mixed mode loading and crack interaction with other discontinuities. The influence of different evaluation parameters is studied by comparing the results with empirical, experimental and alternative numerical solutions and confirms the applicability and capability of the proposed combination of both concepts. Copyright © 2010 John Wiley & Sons, Ltd.