Validation tests of the full-discontinuous Galerkin/extrinsic cohesive law framework of Kirchhoff-Love shells

Due to its ability to account for discontinuities, the discontinuous Galerkin (DG) method presents two main advantages for modeling crack initiations and propagation. On the one hand, it provides an easy way to insert the cohesive elements during the simulation and therefore avoids the drawbacks inherent to the use of an extrinsic cohesive law. On the other hand, the capture of complex crack path requires very thin meshes, and the recourse to a parallel implementation of DG formulations exhibits a high scalability of the resolution scheme. Recently, the authors developed such a DG-fracture framework for Kirchhoff-Love shells in the linear and non-linear ranges. They proved that this framework dissipates, during the fracture process, an amount of energy equal to the fracture energy of the material and that the model is able to propagate the crack with the right speed. In this paper, novel numerical benchmarks are presented to validate the method in various fracture conditions. The two first ones include an initial notch and study the fracture propagation under two different dynamic loadings (impact and blast). The two other ones focus on the fragmentation of initially unbroken specimens due to uniform expansion, in order to demonstrate the ability of the new framework to model crack initiations. Results are in all cases in agreement with the ones reported in the literature.     

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.     

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 fracture framework for Euler–Bernoulli beams based on a full discontinuous Galerkin formulation/extrinsic cohesive law combination

A new full Discontinuous Galerkin discretization of Euler–Bernoulli beam is presented. The main interest of this framework is its ability to simulate fracture problems by inserting a cohesive zone model in the formulation. With a classical Continuous Galerkin method, the use of the cohesive zone model is difficult because inserting a cohesive element between bulk elements is not straightforward. On one hand if the cohesive element is inserted at the beginning of the simulation, there is a modification of the structure stiffness and on the other hand inserting the cohesive element during the simulation requires modification of the mesh during computation. These drawbacks are avoided with the presented formulation as the structure is discretized in a stable and consistent way with full discontinuous elements and inserting cohesive elements during the simulation becomes straightforward. A new cohesive law based on the resultant stresses (bending moment and membrane) of the thin structure discretization is also presented. This model allows propagating fracture while avoiding through‐the‐thickness integration of the cohesive law. Tests are performed to show that the proposed model releases, during the fracture process, an energy quantity equal to the fracture energy for any combination of tension‐bending loadings. Copyright © 2010 John Wiley & Sons, Ltd.     

A rate-dependent cohesive model for simulating dynamic crack propagation in brittle materials

Numerical investigations are conducted to simulate high-speed crack propagation in pre-strained PMMA plates. In the simulations, the dynamic material separation is explicitly modeled by cohesive elements incorporating an initially rigid, linear-decaying cohesive law. Initial attempts using a rate-independent cohesive law failed to reproduce available experimental results as numerical crack velocities consistently overestimate experimental observations. As proof of concept, a phenomenological rate-dependent cohesive law, which bases itself on the physics of microcracking, is introduced to modulate the cohesive law with the macroscopic crack velocity. We then generalize this phenomenological approach by establishing a rate-dependent cohesive law, which relates the traction to the effective displacement and rate of change of effective displacement. It is shown that this new model produces numerical results in good agreement with experimental data. The analysis demonstrates that the simulation of high-speed crack propagation in brittle structures necessitates the use of rate-dependent cohesive models, which account for the complicated rate-process of dynamic fracture at the propagating crack tip.     

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.     

A cohesive XFEM model for simulating fatigue crack growth under mixed-mode loading and overloading

Structures are subjected to cyclic loads that can vary in direction and magnitude, causing constant amplitude mode I simulations to be too simplistic. This study presents a new approach for fatigue crack propagation in ductile materials that can capture mixed-mode loading and overloading. The extended finite element method is used to deal with arbitrary crack paths. Furthermore, adaptive meshing is applied to minimize computation time. A fracture process zone ahead of the physical crack tip is represented by means of cohesive tractions from which the energy release rate, and thus the stress intensity factor can be extracted for an elastic-plastic material. The approach is therefore compatible with the Paris equation, which is an empirical relation to compute the fatigue crack growth rate. Two different models to compute the cohesive tractions are compared. First, a cohesive zone model with a static cohesive law is used. The second model is based on the interfacial thick level set method in which tractions follow from a given damage profile. Both models show good agreement with a mode I analytical relation and a mixed-mode experiment. Furthermore, it is shown that the presented models can capture crack growth retardation as a result of an overload.     

A refined diffuse cohesive approach for the failure analysis in quasibrittle materials—part II: Application to plain and reinforced concrete structures

This work presents the numerical application of the diffuse cohesive interface model introduced in the Part I paper to the failure analysis of plain and reinforced concrete structures, subjected to complex loading conditions, inducing mixed‐mode fracture initiation and propagation. With the aim of capturing the interaction between concrete and steel reinforcements, the adopted fracture model is incorporated in a novel, more general numerical framework for the nonlinear analysis of reinforced concrete structures. Such a framework includes a newly proposed embedded truss model for the reinforcing bars, allowing them to be crossed by the neighboring propagating cracks. Comparisons with available experimental results are provided, assessing the reliability and the numerical accuracy of the proposed concrete model, with reference to plain specimens subjected to single‐crack propagation as well as to reinforced elements subjected to multiple cracking.     

A micromechanical model for a viscoelastic cohesive zone

A micromechanical model for a viscoelastic cohesive zone is formulated herein. Care has been taken in the construction of a physically-based continuum mechanics model of the damaged region ahead of the crack tip. The homogenization of the cohesive forces encountered in this region results in a damage dependent traction-displacement law which is both single integral and internal variable-type. An incrementalized form of this traction-displacement law has been integrated numerically and placed within an implicit finite element program designed to predict crack propagation in viscoelastic media. This research concludes with several example problems on the response of this model for various displacement boundary conditions.     

A computationally efficient cohesive zone model for fatigue

A cohesive zone model has been developed for the simulation of both high and low cycle fatigue crack growth. The developed model provides an alternative approach that reflects the computational efficiency of the well‐established envelop‐load damage model yet can deliver the accuracy of the equally well‐established loading‐unloading hysteresis damage model. A feature included in the new cohesive zone model is a damage mechanism that accumulates as a result of cyclic plastic separation and material deterioration to capture a finite fatigue life. The accumulation of damage is reflected in the loading‐unloading hysteresis curve, but additionally, the model incorporates a fast‐track feature. This is achieved by “freezing in” a particular damage state for one loading cycle over a predefined number of cycles. The new model is used to simulate mode I fatigue crack growth in austenitic stainless steel 304 at significant reduction in the computational cost.     

Quasi‐static crack propagation in plane and plate structures using set‐valued traction‐separation laws

We introduce a numerical technique to model set‐valued traction‐separation laws in plate bending and also plane crack propagation problems. By using of recent developments in thin (Kirchhoff–Love) shell models and the extended finite element method, a complete and accurate algorithm for the cohesive law is presented and is used to determine the crack path. The cohesive law includes softening and unloading to origin, adhesion and contact. Pure debonding and contact are obtained as particular (degenerate) cases. A smooth root‐finding algorithm (based on the trust‐region method) is adopted. A step‐driven algorithm is described with a smoothed law which can be made arbitrarily close to the exact non‐smooth law. In the examples shown the results were found to be step‐size insensitive and accurate. In addition, the method provides the crack advance law, extracted from the cohesive law and the absence of stress singularity at the tip. Copyright © 2007 John Wiley & Sons, Ltd.     

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.     

Numerical techniques for the analysis of crack propagation in cohesive materials

The aim of the paper is the development, assessment and use of suitable numerical procedures for the analysis of the crack evolution in cohesive materials. In particular, homogeneous as well as heterogeneous materials, obtained by embedding short stiff fibres in a cohesive matrix, are considered. Two‐dimensional Mode I fracture problems are investigated. The cohesive constitutive law is adopted to model the process zone occurring at the crack tip. An elasto‐plastic constitutive relationship, able to take into account the processes of fibre debonding and pull‐out, is introduced to model the mechanical response of the short fibres. Two numerical procedures, based on the stress and on the energy approach, are developed to investigate the crack propagation in cohesive as well as fibre‐reinforced materials, characterized by a periodic crack distribution. The results obtained using the stress and energy approaches are compared in order to evaluate the effectiveness of the procedures. Investigations on the size effect for microcracked periodic cohesive materials, and on the beneficial effects of the fibres in improving the composite material response, are developed. Copyright © 2003 John Wiley & Sons, Ltd.     

Effect of residual stress on cleavage fracture toughness by using cohesive zone model

This study presents the effect of residual stresses on cleavage fracture toughness by using the cohesive zone model under mode I, plane stain conditions. Modified boundary layer simulations were performed with the remote boundary conditions governed by the elastic K‐field and T‐stress. The eigenstrain method was used to introduce residual stresses into the finite element model. A layer of cohesive elements was deployed ahead of the crack tip to simulate the fracture process zone. A bilinear traction–separation‐law was used to characterize the behaviour of the cohesive elements. It was assumed that the initiation of the crack occurs when the opening stress drops to zero at the first integration point of the first cohesive element ahead of the crack tip. Results show that tensile residual stresses can decrease the cleavage fracture toughness significantly. The effect of the weld zone size on cleavage fracture toughness was also investigated, and it has been found that the initiation toughness is the linear function of the size of the geometrically similar weld. Results also show that the effect of the residual stress is stronger for negative T‐stress while its effect is relatively smaller for positive T‐stress. The influence of damage parameters and material hardening was also studied.     

Modelling Mesh Independent Transverse Cracks in Laminated Composites with a Simplified Cohesive Segment Method

A methodology is proposed for modelling transverse matrix cracks in laminated composites in a three-dimensional explicit finite element analysis framework. The method is based on the introduction of extra degrees of freedom to represent the displacement discontinuity and the use of a cohesive zone model to determine damage evolution and crack propagation. The model is designed for the analysis of matrix cracks in laminates made of uni-directional fibre-reinforced plies, allowing several assumptions to be made which greatly simplify the algorithm. This was implemented in the commercial software Abaqus/Explicit as a user-defined element subroutine (VUEL). The methodology was verified via the analysis of open-hole tension tests considering both ±45˚ and quasi-isotropic layups. The results were found to be in qualitative agreement with experimental observations in terms of the nucleation and propagation of matrix cracks, the progressive delamination behaviour and the evident interactions between these damage mechanisms.     

Application of the cohesive zone model for analysing the edge crack propagation of steel sheet in the cold rolling process

In the cold rolling process, the expansion and coalescence of micro‐defects can make steel sheet quality descend and create edge crack in the steel sheet. And the edge crack can cause the strip rupture completely. In this research, the cohesive zone model (CZM) was used to analyse the initiation and propagation of edge crack in the cold rolling process with the non‐reversing two‐high mill. A bi‐linear traction–separation law was utilized which is primarily given by the CZM parameters including the cohesive stress, T, and the cohesive energy, Γ. Compared with other popular models such as the Gurson–Tvergaard–Needleman (GTN) model, the CZM presents certain advantages because it requires a smaller number of parameters to be defined. Comparison results of the experiments and simulation illustrated that the CZM can provide accurate prediction for the propagation of edge crack in the cold rolling process. Parametric analysis was carried out and showed that the extent of the crack propagation increases with the increasing of the reduction ratio.     

A fractional rate‐dependent cohesive‐zone model

This paper presents a novel formulation of a hereditary cohesive zone model able to effectively capture rate‐dependent crack propagation along a defined interface, over a wide range of applied loading rates and with a single set of seven input parameters only, as testified by the remarkable agreement with experimental results in the case of a double cantilever beam made of steel adherends bonded along a rubber interface. The formulation relies on the assumption that the measured fracture energy is the sum of a rate‐independent ‘rupture’ energy, related to the rupture of primary bonds at the atomic or molecular level, and of additional dissipation caused by other rate‐dependent dissipative mechanisms present in the material and occurring simultaneously to rupture. The first contribution is accounted for by introducing a damage‐type internal variable, whose evolution follows a rate‐independent law for consistency with the assumption of rate independence of the rupture energy. To account for the additional dissipation, a fractional‐calculus‐based linear viscoelastic model is used, because for many polymers, it is known to capture the material response within an extremely wide range of strain rates much more effectively than classic models based on an exponential kernel. To the authors' knowledge, this is the first application of fractional viscoelasticity to the simulation of fracture. © 2015 The Authors. International Journal for Numerical Methods in Engineering published by John Wiley & Sons Ltd.     

A partition of unity‐based multiscale approach for modelling fracture in piezoelectric ceramics

The development of models for a priori assessment of the reliability of micro electromechanical systems is of crucial importance for the further development of such devices. In this contribution a partition of unity‐based cohesive zone finite element model is employed to mimic crack nucleation and propagation in a piezoelectric continuum. A multiscale framework to appropriately represent the influence of the microstructure on the response of a miniaturized component is proposed. It is illustrated that using the proposed multiscale method a representative volume element exists. Numerical simulations are performed to demonstrate the constitutive homogenization framework. Copyright © 2009 John Wiley & Sons, Ltd.     

On damage accumulations in the cyclic cohesive zone model for XFEM analysis of mixed-mode fatigue crack growth

Predicting mixed-mode fatigue crack propagation is an important and troublesome issue in structure assessment for decades. In the present paper an extended finite element method (XFEM) combined with a new cyclic cohesive zone model (CCZM) is introduced for simulating fatigue crack propagation under mixed-mode loading conditions, which has been implemented in the commercial general purpose software ABAQUS. The algorithm allows introducing a new crack surface at arbitrary locations and directions in a finite element mesh, without re-meshing. The cyclic cohesive zone model is based on the known S–N curves and Goodman diagram for metallic materials and validated by uniaxial tension results. Furthermore, the sensitivity of the model parameter is investigated for mixed-mode fatigue. The virtual crack closure technique has been extended to the cohesive zone model and proposed to calculate the energy release rate for the generalized Paris’ law. Finally, the crack propagation rate and direction under mixed-mode fatigue loading conditions are studied.