Multiple cohesive crack growth in brittle materials by the extended Voronoi cell finite element model

This paper is aimed at modeling the propagation of multiple cohesive cracks by the extended Voronoi cell finite element model or X-VCFEM. In addition to polynomial terms, the stress functions in X-VCFEM 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. Cracks are modeled by an extrinsic cohesive zone model in this paper. The incremental crack propagation direction and length are adaptively determined by a cohesive fracture energy based criterion. Numerical examples are solved and compared with existing solutions in the literature to validate the effectiveness of X-VCFEM. The effect of cohesive zone parameters on crack propagation is studied. Additionally, the effects of morphological distributions such as length, orientation and dispersion on crack propagation are studied.     

Voronoi cell finite element model based on micropolar theory of thermoelasticity for heterogeneous materials

In this paper, a new ‘Voronoi cell finite element model’ is developed for solving steady-state heat conduction and micropolar thermoelastic stress analysis problems in arbitrary heterogeneous materials. The method is based on the natural discretization of a multiple phase domain into basic structural elements by Dirichlet Tessellation. Tessellation process results in a network of polygons called Voronoi polygons. In this paper, formulations are developed for treating these polygons as elements in a finite element mesh. Furthermore, a composite Voronoi cell finite element model is developed to account for the presence of a second phase inclusion within a polygonal element. Various numerical examples are executed for validating the effectiveness of this model in the analysis of the temperature and stress fields for micropolar elastic materials. Effective material properties are derived for microstructures containing different distributions of second phase.     

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.     

Modeling of interfacial debonding crack in particle reinforced composites using Voronoi cell finite element method

In this paper, Voronoi cell finite element method (VCFEM), introduced by Ghosh and coworkers (1993), is applied to describe the matrix-inclusion interfacial debonding for particulate reinforced composites. In proposed VCFEM, the damage initiation is simulated by partly debonding of the interface under the assumption of the critical normal stress law, and gradual matrix-inclusion separations are simulated with an interface remeshing method that a critical interfacial node at the crack tip is replaced by a node pairs along the debonded matrix-inclusion interface and a more pair of nodes are needed to be added on the crack interface near the crack tip in order to better facilitate the free-traction boundary condition and the jumps of solution. The comparison of the results of proposed VCFEM and commercial finite element packages MARC and ABAQUS. Examples have been given for a single inclusion of gradually interfacial debonding and for a complex structure with 20 inclusions to describe the interfacial damage under plane stress conditions. Good agreements are obtained between the VCFEM and the general finite element method. It appears that this method is a more efficient way to deal with the interfacial damage of composite materials. The financial support by the Special Funds for the National Major Fundamental Research Projects G19990650 and the National Natural Science Foundation of China No. 59871022 are gratefully acknowledged.     

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.     

基于参数变分原理的含夹杂Voronoi单元法及非均质材料弹塑性计算

含夹杂Voronoi单元通过在基体单元中引入一任意夹杂, 可以更好地反映非均质材料中微结构特性。基于参数势能和余能原理, 推导了无夹杂和含夹杂Voronoi单元有限元列式, 并在此基础上形成二次规划求解模型。将含夹杂Voronoi单元应用于非均质材料宏观弹塑性性能预测计算中, 分析了非均质材料中夹杂对其宏观等效弹塑性力学性能的影响。数值结果与其它方法所得结果的比较证明了本文中所给出模型的正确性和工程可适用性。      

Computational analysis of mixed-mode fatigue crack growth in quasi-brittle materials using extended finite element methods

In the present paper an extended finite element method (XFEM) containing strong discontinuity within elements is introduced and implemented in the commercial general purpose software ABAQUS. The algorithm allows introducing a new crack surface at arbitrary locations and directions in elements. To consider fatigue crack nucleation and propagation in quasi-brittle materials the XFEM is combined with a cyclic cohesive model. Accumulative material damage is described by separate evolution equations. The crack path is completely independent of the mesh structure but determined by the mixed-mode loading cases. Numerical simulations illustrate the ability of this method to simulate fracture with unstructured meshes. The computational results agree with known fracture experiment data. Known fatigue observations can be predicted using the present model.     

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.     

Centroidal Voronoi tessellation‐based finite element superconvergence

In this article, a finding on finite element superconvergence is reported. The Laplacian operator with Dirichlet boundary condition is considered. The linear finite element solutions have an O(h^2+α)(α≈0.5)‐superconvergence in l₂ norm at nodes on an almost equilateral triangular mesh generated based on centroidal Voronoi tessellation, for an arbitrary 2D bounded domain. Extensive numerical examples are presented to demonstrate the superconvergence property. Copyright © 2008 John Wiley & Sons, Ltd.     

Application of cohesive zone model in crack propagation analysis in multiphase composite materials

A nonlinear cohesive stress distribution function is employed by relating the cohesive stress to the cohesive zone size (CZS) and the distance from the crack tip to investigate the elastic-plastic fracture behaviors. A crack-inclusion interaction problem is taken as an example to explore the fracture process in the cohesive zone area. The CZS and crack surface opening displacement are evaluated numerically. It is found that for different cohesive parameter combinations, the normalized CZS and crack surface opening displacements change drastically. By reducing the current model to the famous Dugdale model, the results obtained match well with the existing ones.     

An enriched cohesive zone model for delamination in brittle interfaces

Application of standard cohesive zone models in a finite element framework to simulate delamination in brittle interfaces may trigger non‐smooth load–displacement responses that lead to the failure of iterative solution procedures. This non‐smoothness is an artifact of the discretization; and hence it can be avoided by sufficiently refining the mesh leading to unacceptably high computational costs and a low efficiency and robustness. In this paper, a process‐driven hierarchical extension is proposed to enrich the separation approximation in the process zone of a cohesive crack. Some numerical examples show that instead of mesh refinement, a more efficient enriched formulation can be used to prevent a non‐smooth solution. Copyright © 2009 John Wiley & Sons, Ltd.     

Extended finite element method and fast marching method for three-dimensional fatigue crack propagation

A numerical technique for planar three-dimensional fatigue crack growth simulations is proposed. The new technique couples the extended finite element method (X-FEM) to the fast marching method (FMM). In the X-FEM, a discontinuous function and the two-dimensional asymptotic crack-tip displacement fields are added to the finite element approximation to account for the crack using the notion of partition of unity. This enables the domain to be modeled by finite elements with no explicit meshing of the crack surfaces. The initial crack geometry is represented by level set functions, and subsequently signed distance functions are used to compute the enrichment functions that appear in the displacement-based finite element approximation. The FMM in conjunction with the Paris crack growth law is used to advance the crack front. Stress intensity factors for planar three-dimensional cracks are computed, and fatigue crack growth simulations for planar cracks are presented. Good agreement between the numerical results and theory is realized.     

An augmented cohesive zone element for arbitrary crack coalescence and bifurcation in heterogeneous materials

We demonstrate that traditional cohesive zone (CZ) elements cannot be accurate when used in conjunction with solid elements with arbitrary intra‐element cracking capability, because they cannot capture the load transfer between cohesive interfaces and the solid elements when crack bifurcation or coalescence occurs. An augmented cohesive zone (ACZ) element based on the augmented finite element method formulation is therefore proposed. The new element allows for arbitrary separation of the cohesive element in accordance with the crack configuration of the abutting solid elements, thus correctly maintaining the non‐linear coupling between merging or bifurcating cracks. Numerical accuracy and effectiveness of the proposed ACZ element are demonstrated through several examples. Copyright © 2011 John Wiley & Sons, Ltd.     

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 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.     

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.     

Extrinsic cohesive modelling of dynamic fracture and microbranching instability in brittle materials

Dynamic crack microbranching processes in brittle materials are investigated by means of a computational fracture mechanics approach using the finite element method with special interface elements and a topological data structure representation. Experiments indicate presence of a limiting crack speed for dynamic crack in brittle materials as well as increasing fracture resistance with crack speed. These phenomena are numerically investigated by means of a cohesive zone model (CZM) to characterize the fracture process. A critical evaluation of intrinsic versus extrinsic CZMs is briefly presented, which highlights the necessity of adopting an extrinsic approach in the current analysis. A novel topology‐based data structure is employed to enable fast and robust manipulation of evolving mesh information when extrinsic cohesive elements are inserted adaptively. Compared to intrinsic CZMs, which include an initial hardening segment in the traction–separation curve, extrinsic CZMs involve additional issues both in implementing the procedure and in interpreting simulation results. These include time discontinuity in stress history, fracture pattern dependence on time step control, and numerical energy balance. These issues are investigated in detail through a ‘quasi‐steady‐state’ crack propagation problem in polymethylmethacrylate. The simulation results compare reasonably well with experimental observations both globally and locally, and demonstrate certain advantageous features of the extrinsic CZM with respect to the intrinsic CZM. Copyright © 2007 John Wiley & Sons, Ltd.     

基于小波内聚力模型的界面裂纹扩展

利用区间B样条小波良好的局部化性能,将内聚力模型（CZM）引入小波有限元法（WFEM）数值分析中,以区间B样条小波尺度函数作为插值函数,构造小波内聚力界面单元,推导了小波内聚力界面单元刚度矩阵,基于虚拟裂纹闭合技术（VCCT）计算界面裂纹应变能释放率（SERR）,采用β-Κ断裂准则,实现界面裂纹扩展准静态分析。将WFEM和传统有限元法（CFEM） 的SERR数值分析结果与理论解进行比较,结果表明:采用WFEM和CFEM计算的SERR分别为96.60 J/m2 和 101.43 J/m2,2种方法的SERR数值解与理论解相对误差分别为1.85%和3.06%,这明确表明WFEM在计算界面裂纹扩展方面能用较少单元和节点数获得较高的计算精度和效率。在此基础上,探讨了界面裂纹初始长度和双材料弹性模量比对界面裂纹扩展的影响,分析结果表明:界面裂纹尖端等效应力随界面裂纹初始长度的增加而增加;双材料弹性模量比相差越大,界面裂纹越易于扩展,且裂纹扩展长度也越大,因此可通过调节双材料弹性模量比来延缓界面裂纹扩展。