Combined continuum damage‐embedded discontinuity model for explicit dynamic fracture analyses of quasi‐brittle materials

In this paper, a novel constitutive model combining continuum damage with embedded discontinuity is developed for explicit dynamic analyses of quasi‐brittle failure phenomena. The model is capable of describing the rate‐dependent behavior in dynamics and the three phases in failure of quasi‐brittle materials. The first phase is always linear elastic, followed by the second phase corresponding to fracture‐process zone creation, represented with rate‐dependent continuum damage with isotropic hardening formulated by utilizing consistency approach. The third and final phase, involving nonlinear softening, is formulated by using an embedded displacement discontinuity model with constant displacement jumps both in normal and tangential directions. The proposed model is capable of describing the rate‐dependent ductile to brittle transition typical of cohesive materials (e.g., rocks and ice). The model is implemented in the finite element setting by using the CST elements. The displacement jump vector is solved for implicitly at the local (finite element) level along with a viscoplastic return mapping algorithm, whereas the global equations of motion are solved with explicit time‐stepping scheme. The model performance is illustrated by several numerical simulations, including both material point and structural tests. The final validation example concerns the dynamic Brazilian disc test on rock material under plane stress assumption. Copyright © 2014 John Wiley & Sons, Ltd.     

Analysis of three‐dimensional crack initiation and propagation using the extended finite element method

An Erratum has been published for this article in International Journal for Numerical Methods in Engineering 2005, 63(8): 1228. We present a new formulation and a numerical procedure for the quasi‐static analysis of three‐dimensional crack propagation in brittle and quasi‐brittle solids. The extended finite element method (XFEM) is combined with linear tetrahedral elements. A viscosity‐regularized continuum damage constitutive model is used and coupled with the XFEM formulation resulting in a regularized ‘crack‐band’ version of XFEM. The evolving discontinuity surface is discretized through a C⁰ surface formed by the union of the triangles and quadrilaterals that separate each cracked element in two. The element's properties allow a closed form integration and a particularly efficient implementation allowing large‐scale 3D problems to be studied. Several examples of crack propagation are shown, illustrating the good results that can be achieved. Copyright © 2005 John Wiley & Sons, Ltd.     

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.     

Failure analysis of quasi‐brittle materials involving multiple mechanisms on fractured surfaces

We present a variational formulation for the quasi‐static boundary value problem of a structure with quasi‐brittle materials, involving (i) unknown states of contact, (ii) deformation‐dependent frictional forces, (iii) crack opening and closing with cohesive traction, and (iv) configuration change due to the initiation and the evolution of cracks, and propose a new finite cover method (FCM) capable of reflecting those multiple mechanisms in the failure analysis. The cover‐division strategy is also introduced to judge the generation of cracks, and to locate and orient them within the framework of the FCM. A relevant numerical algorithm is designed to be consistent with the mathematical representation of the multiple mechanisms. Several numerical examples are presented to validate the proposed method and to demonstrate its promise and potential for evaluating the ultimate strength of quasi‐brittle materials. Copyright © 2006 John Wiley & Sons, Ltd.     

Complete Tangent Stiffness for eXtended Finite Element Method by including crack growth parameters

The eXtended Finite Element Method (XFEM) is a useful tool for modeling the growth of discrete cracks in structures made of concrete and other quasi‐brittle and brittle materials. However, in a standard application of XFEM, the tangent stiffness is not complete. This is a result of not including the crack geometry parameters, such as the crack length and the crack direction directly in the virtual work formulation. For efficiency, it is essential to obtain a complete tangent stiffness. A new method in this work is presented to include an incremental form the crack growth parameters on equal terms with the degrees of freedom in the FEM‐equations. The complete tangential stiffness matrix is based on the virtual work together with the constitutive conditions at the crack tip. Introducing the crack growth parameters as direct unknowns, both equilibrium equations and the crack tip criterion can be handled within the same standard nonlinear iterations. This new solution strategy is believed to provide the modeling capabilities to deal with simultaneous growth of several cracks. A cohesive crack modeling is used. The method is applied to a partly cracked XFEM element of linear strain triangle type with the crack length as the unknown crack growth parameter. In this paper, two examples are given. The first example verifies the theory and the implementation. The second example is the benchmark test three point bending test, where the efficiency of the complete tangential behavior is shown. Copyright © 2013 John Wiley & Sons, Ltd.     

Numerical study on deformations in a cracked viscoelastic body with the extended finite element method

Deformations such as crack opening and sliding displacements in a cracked viscoelastic body are numerically investigated by the extended finite element method (XFEM). The solution is carried out directly in time domain with a mesh not conforming to the crack geometry. The generalized Heaviside function is used to reflect the displacement discontinuity across a crack surface while the basis functions extracted from the viscoelastic asymptotic fields are used to manifest the gradient singularity at a crack tip. With these treatments, the XFEM formulations are derived in an incremental form. In evaluating the stiffness matrix, a selective integration scheme is suggested for problems with high Poisson ratios often encountered in viscoelastic materials over different element types in the XFEM. Numerical examples show that the crack opening displacement and crack sliding displacement are satisfactory.     

A consistent partly cracked XFEM element for cohesive crack growth

Present extended finite element method (XFEM) elements for cohesive crack growth may often not be able to model equal stresses on both sides of the discontinuity when acting as a crack‐tip element. The authors have developed a new partly cracked XFEM element for cohesive crack growth with extra enrichments to the cracked elements. The extra enrichments are element side local and were developed by superposition of the standard nodal shape functions for the element and standard nodal shape functions for a sub‐triangle of the cracked element. With the extra enrichments, the crack‐tip element becomes capable of modelling variations in the discontinuous displacement field on both sides of the crack and hence also capable of modelling the case where equal stresses are present on each side of the crack. The enrichment was implemented for the 3‐node constant strain triangle (CST) and a standard algorithm was used to solve the non‐linear equations. The performance of the element is illustrated by modelling fracture mechanical benchmark tests. Investigations were carried out on the performance of the element for different crack lengths within one element. The results are compared with previously obtained XFEM results applying fully cracked XFEM elements, with computational results achieved using standard cohesive interface elements in a commercial code, and with experimental results. The suggested element performed well in the tests. Copyright © 2007 John Wiley & Sons, Ltd.     

Anomalous behavior of bilinear quadrilateral finite elements for modeling cohesive cracks with XFEM/GFEM

The performance of partition‐of‐unity based methods such as the generalized finite element method or the extended finite element method is studied for the simulation of cohesive cracking. The focus of investigation is on the performance of bilinear quadrilateral finite elements using these methods. In particular, the approximation of the displacement jump field, representing cohesive cracks, by extended finite element method/generalized finite element method and its effect on the overall behavior at element and structural level is investigated. A single element test is performed with two different integration schemes, namely the Newton‐Cotes/Lobatto and the Gauss integration schemes, for the cracked interface contribution. It was found that cohesive crack segments subjected to a nonuniform opening in unstructured meshes (or an inclined crack in a structured finite element mesh) result in an unrealistic crack opening. The reasons for such behavior and its effect on the response at element level are discussed. Furthermore, a mesh refinement study is performed to analyze the overall response of a cohesively cracked body in a finite element analysis. Copyright © 2013 John Wiley & Sons, Ltd.     

Extended finite element method coupled with face‐based strain smoothing technique for three‐dimensional fracture problems

In this work, the extended finite element method (XFEM) is for the first time coupled with face‐based strain‐smoothing technique to solve three‐dimensional fracture problems. This proposed method, which is called face‐based smoothed XFEM here, is expected to combine both the advantages of XFEM and strain‐smoothing technique. In XFEM, arbitrary crack geometry can be modeled and crack advance can be simulated without remeshing. Strain‐smoothing technique can eliminate the integration of singular term over the volume around the crack front, thanks to the transformation of volume integration into area integration. Special smoothing scheme is implemented in the crack front smoothing domain. Three examples are presented to test the accuracy, efficiency, and convergence rate of the face‐based smoothed XFEM. From the results, it is clear that smoothing technique can improve the performance of XFEM for three‐dimensional fracture problems. Copyright © 2015 John Wiley & Sons, Ltd.     

基于扩展有限元法的钢筋混凝土梁复合断裂过程模拟研究

针对钢筋混凝土梁裂纹扩展问题,基于扩展有限元法,建立了预置裂纹的简支混凝土梁三维模型,用粘聚裂纹模型描述裂纹面间的力学行为,采用线性的软化曲线表示裂纹尖端断裂过程区的应变软化行为,分别对素混凝土梁和钢筋混凝土梁的复合断裂过程进行模拟,分析了纵向钢筋对裂纹扩展路径、荷载-挠度和荷载 -CMOD (裂缝开口处张开位移)曲线的影响,并与文献中的试验结果进行对比,计算结果与试验结果吻合良好,展示了扩展有限元法在结构断裂破坏分析方面的独特优势。     

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.     

Mesh‐independent matrix cracking and delamination modeling in laminated composites

The initiation and evolution of transverse matrix cracks and delaminations are predicted within a mesh‐independent cracking (MIC) framework. MIC is a regularized extended finite element method (x‐FEM) that allows the insertion of cracks in directions that are independent of the mesh orientation. The Heaviside step function that is typically used to introduce a displacement discontinuity across a crack surface is replaced by a continuous function approximated by using the original displacement shape functions. Such regularization allows the preservation of the Gaussian integration schema regardless of the enrichment required to model cracking in an arbitrary direction. The interaction between plies is anchored on the integration point distribution, which remains constant through the entire simulation. Initiation and propagation of delaminations between plies as well as intra‐ply MIC opening is implemented by using a mixed‐mode cohesive formulation in a fully three‐dimensional model that includes residual thermal stresses. The validity of the proposed methodology was tested against a variety of problems ranging from simple evolution of delamination from existing transverse cracks to strength predictions of complex laminates withouttextita priori knowledge of damage location or initiation. Good agreement with conventional numerical solutions and/or experimental data was observed in all the problems considered. Published 2011. This article is a US Government work and is in the public domain in the USA.     

Loading‐rate dependence of mode I crack growth in concrete

Mode I crack propagation process of concrete under relatively low loading rates which cover four orders of magnitude (0.2 μm/s to 2.0 mm/s) is investigated with three‐point bending (TPB) beams. All measured material properties exhibit rate sensitivity and follow a log‐linear relationship with the loading rate. A rate‐sensitive softening curve is established. The complete load‐crack mouth opening displacement (P‐CMOD) curve, crack propagation length, and fracture process zone (FPZ) length are simulated based on crack growth criterion with the fitted material parameters under those loading rates. Results show that the simulated P‐CMOD curves agree well with those of experimental measurements. It is clear that the peak load increases with the loading rate and so is the critical crack mouth opening displacement. Moreover, under the same load level, the length of the FPZ and the cohesive stress at the initial crack tip also increase with the increasing loading rate.     

Two‐scale diffusion–deformation coupling model for material deterioration involving micro‐crack propagation

In this paper, we introduce a two‐scale diffusion–deformation coupled model that represents the aging material deterioration of two‐phase materials involving micro‐crack propagations. The mathematical homogenization method is applied to relate the micro‐ and macroscopic field variables, and a weak coupling solution method is employed to solve the two‐way coupling phenomena between the diffusion of scalar fields and the deformation of quasi‐brittle solids. The macroscopic mechanical behavior represented by the derived two‐scale two‐way coupled model reveals material nonlinearity due to micro‐scale cracking induced by the scalar‐field‐induced deformation, which can be simulated by the finite cover method. After verifying the fundamental validity of the proposed model and the analysis method, we perform a simple numerical example to demonstrate their ability to predict aging material deterioration. Copyright © 2010 John Wiley & Sons, Ltd.     

A continuum damage mechanics method for fracture simulation of quasi‐brittle materials

This paper addresses a novel continuum damage‐based method for simulating failure process of quasi‐brittle materials starting from local damage initiation to final fracture. In the developed method, the preset characteristic length field is used to evaluate damage instead of element, which is used to reduce the spurious sensitivity. In addition, damage is only updated in the most dangerous location at a time for considering stress redistribution due to damage evolution, which is used to simulate competitive fracture process. As cases study, representative numerical simulations of two benchmark tests are given to verify the performance of the developed continuum damage‐based method together with a used damage model. The simulation results of the crack paths for two concrete specimens obtained from the developed method matched well with the corresponding experimental results. The results show that the developed continuum damage‐based method is effective and can be used to simulate damage and fracture process of brittle or quasi‐brittle materials. And the simulation results based on the developed method depend only the preset characteristic length field and not grid mesh.     

不同软化曲线形状对裂缝扩展阻力G_R曲线的影响

最近基于裂缝粘聚力提出了描述裂缝扩展全过程阻力变化的GR曲线,所发展的解析解相关于骨料桥联咬合作用造成的非线性断裂过程区的能量耗散,故与该区域材料使用的拉伸软化本构关系即应力-裂缝张开口位移软化曲线有密切联系。鉴于此,该文采用4种不同的软化曲线研究了软化曲线形状对裂缝扩展阻力GR的影响。结果发现,GR曲线对软化曲线敏感,GR曲线的合理性依赖于软化曲线的准确性。在使用准确软化曲线的前提下,GR曲线的特征点与软化曲线的特征点存在相对应关系。     

A comparative study on finite element methods for dynamic fracture

The performance of finite element methods for dynamic crack propagation in brittle materials is studied. Three methods are considered: the extended finite element method (XFEM), element deletion method and interelement crack method. The extended finite element method is a method for arbitrary crack propagation without remeshing. In element deletion methods, elements that meet a fracture criterion are deleted. In interelement crack methods, the crack is limited to element edges; the separation of these edges is governed by a cohesive law. We show that XFEM and interelement method show similar crack speeds and crack paths. However, both fail to predict a benchmark experiment without adjustment of the energy release rate. The element deletion method performs very poorly for the refinements studied, and is unable to predict crack branching.     

Determination of the kink point in the bilinear softening model for concrete

The characterization of the softening curve from experimental results is essential for predicting the fracture behavior of quasi-brittle materials like concrete. Among various shapes (e.g. linear, exponential) to describe the softening behavior of concrete, the bilinear softening relationship has been extensively used and is the model of choice in this work. Currently, there is no consensus about the location of the kink point in the bilinear softening curve. In this study, the location of the kink point is proposed to be the stress at the critical crack tip opening displacement. Experimentally, the fracture parameters required to describe the bilinear softening curve can be determined with the “two-parameter fracture model” and the total work of fracture method based on a single concrete fracture test. The proposed location of the kink point compares well with the range of kink point locations reported in the literature, and is verified by plotting stress profiles along the expected fracture line obtained from numerical simulations with the cohesive zone model. Finally, prediction of experimental load versus crack mouth opening displacement curves validate the proposed location of the kink point for different concrete mixtures and also for geometrically similar specimens with the same concrete mixture. The experiments were performed on three-point bending specimens with concrete mixtures containing virgin coarse aggregate, recycled concrete coarse aggregate (RCA), and a 50-50 blend of RCA and virgin coarse aggregate. The verification and validation studies support the hypothesis of the kink point occurring at the critical crack tip opening displacement.     

用扩展有限元方法模拟混凝土的复合型开裂过程

用扩展有限元法对混凝土梁复合型开裂过程进行了数值模拟。裂纹面间的力学行为采用粘聚裂纹模型来描述,通过引入切向保留刚度考虑剪力分量的影响。开裂方向的计算采用了一种简化的最大切向应力准则。对Arrea和Ingraffea的混凝土梁复合开裂实验进行了数值模拟。计算给出了裂纹萌生、扩展的过程及破坏形态,并获得了与实验结果对比良好的荷载-裂纹开口滑移曲线。结果表明,扩展有限元法通过附加特定的位移模式,使裂纹两侧不连位移场的表达独立于网格划分,是一种能够模拟准脆性材料复合开裂问题的有效方法。     

Three‐dimensional non‐planar crack growth by a coupled extended finite element and fast marching method

A numerical technique for non‐planar three‐dimensional linear elastic crack growth simulations is proposed. This technique couples the extended finite element method (X‐FEM) and the fast marching method (FMM). In crack modeling using X‐FEM, the framework of partition of unity is used to enrich the standard finite element approximation by a discontinuous function and the two‐dimensional asymptotic crack‐tip displacement fields. The initial crack geometry is represented by two level set functions, and subsequently signed distance functions are used to maintain the location of the crack and to compute the enrichment functions that appear in the displacement approximation. Crack modeling is performed without the need to mesh the crack, and crack propagation is simulated without remeshing. Crack growth is conducted using FMM; unlike a level set formulation for interface capturing, no iterations nor any time step restrictions are imposed in the FMM. Planar and non‐planar quasi‐static crack growth simulations are presented to demonstrate the robustness and versatility of the proposed technique. Copyright © 2008 John Wiley & Sons, Ltd.