Computational modeling of mixed-mode fatigue crack growth using extended finite element methods

The extended finite element method (XFEM) combined with a cyclic cohesive zone model (CCZM) is discussed and implemented for analysis of fatigue crack propagation under mixed-mode loading conditions. Fatigue damage in elastic-plastic materials is described by a damage evolution equation in the cohesive zone model. Both the computational implementation and the CCZM are investigated based on the modified boundary layer formulation under mixed-mode loading conditions. Computational results confirm that the maximum principal stress criterion gives accurate predictions of crack direction in comparison with known experiments. Further popular multi-axial fatigue criteria are compared and discussed. Computations show that the Findley criterion agrees with tensile stress dominant failure and deviates from experiments for shear failure. Furthermore, the crack propagation rate under mixed mode loading has been investigated systematically. It is confirmed that the CCZM can agree with experiments.     

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.     

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.     

Interface Cohesive Elements to Model Matrix Crack Evolution in Composite Laminates

In this paper, the transverse matrix (resin) cracking developed in multidirectional composite laminates loaded in tension was numerically investigated by a finite element (FE) model implemented in the commercially available software Abaqus/Explicit 6.10. A theoretical solution using the equivalent constraint model (ECM) of the damaged laminate developed by Soutis et al. was employed to describe matrix cracking evolution and compared to the proposed numerical approach. In the numerical model, interface cohesive elements were inserted between neighbouring finite elements that run parallel to fibre orientation in each lamina to simulate matrix cracking with the assumption of equally spaced cracks (based on experimental measurements and observations). The stress based traction-separation law was introduced to simulate initiation of matrix cracking and propagation under mixed-mode loading. The numerically predicted crack density was found to depend on the mesh size of the model and the material fracture parameters defined for the cohesive elements. Numerical predictions of matrix crack density as a function of applied stress are in a good agreement to experimentally measured and theoretically (ECM) obtained values, but some further refinement will be required in near future work.     

Prediction of 3D small fatigue crack propagation in shot-peened specimens

Shot peening has been widely applied in industrial design to improve fatigue durability of high loaded machine components. The compressive residual stress induced by shot peening is in general assumed to be responsible for the improvement of material fatigue strength. In the present work a cyclic cohesive zone model is extended to analyze three-dimensional fatigue crack growth in shot-peened specimens. Fatigue crack growth behaviors in both unpeened and peened specimens are investigated using 3D finite element analysis. The parameters of the cohesive zone model have been identified in 2D unpeened specimens and are applied to predict peened specimens directly. The results indicate that shot peening strongly affects crack initiation time and crack profiles, but has little effect on crack propagation rate. It implies that the shot peening will hardly change Paris’ law used for the damage tolerant design.     

Mixed-mode quasi-static failure criteria for adhesively-bonded pultruded GFRP joints

The strain energy release rates of adhesively-bonded pultruded GFRP joints were determined experimentally. The crack propagated in the adherend along paths outside the symmetry plane accompanied by fiber bridging. A new method, designated the “extended global method”, was introduced to facilitate mode partitioning in the mixed-mode experiments. Non-linear finite element models were developed in order to quantify the effect of the observed fiber bridging on crack propagation. An exponential traction-separation cohesive law was used to model the fiber bridging zone and calculate the energy release rate due to the fiber bridging, while the virtual crack closure technique was used for calculation of the fracture components at the crack tip. Experimental, analytical and numerical analyses were used to establish quasi-static mixed-mode failure criteria for crack initiation and propagation. The derived mixed-mode failure criteria can be used for simulating progressive crack propagation in other joint configurations comprising the same adhesive and adherends.     

Interface crack initiation at V-notches along adhesive bonding in weakly bonded polymers subjected to mixed-mode loading

In this paper, interface crack initiation at V-notches along adhesive in bonded Polycarbonate (PC) and Poly Methyl Methacrylate (PMMA) subjected to mixed-mode loading conditions was investigated based on a combined experimental, finite element and matched asymptotic analysis. The V-notch specimens with an adhesive interface starting from its tip made at different notch angles were tested under three-point bending conditions. The experimental observations show that the specimens mainly fail by cracks along the interface. Also, the load at the crack initiation increases when the notch angle increases. The computational results are then used to explain and to correlate with the experimental data. A two-fold criterion developed by Leguillon (Eur J Mech A/Solids 21:61?C72 2002) that requires a simultaneous satisfaction of both Griffith energy and stress conditions for the crack initiation at a notch in the specimen made of a homogeneous brittle material is first extended for V-notch specimens under mixed-mode loading conditions and then used to estimate the crack initiation load. The estimated loads appear to agree well with the experimental data. Finally, an inverse method is proposed to estimate the values of fracture toughness at different mode mixity ratios.     

Modeling of cohesive crack growth using an adaptive mesh refinement via the modified-SPR technique

In this paper, an adaptive finite element procedure is presented in modeling of mixed-mode cohesive crack propagation via the modified superconvergent path recovery technique. The adaptive mesh refinement is performed based on the Zienkiewicz–Zhu error estimator. The weighted-SPR recovery technique is employed to improve the accuracy of error estimation. The Espinosa–Zavattieri bilinear cohesive zone model is applied to implement the traction-separation law. It is worth mentioning that no previous information is necessary for the path of crack growth and no region of the domain is necessary to be filled by the cohesive elements. The maximum principal stress criterion is employed for predicting the direction of extension of the cohesive crack in order to implement the cohesive elements. Several numerical examples are analyzed numerically to demonstrate the capability and efficiency of proposed computational algorithm.     

Low‐cycle fatigue delamination initiation and propagation in fibre metal laminates

The aim of this paper is to develop an elastic–plastic‐damage constitutive law and a tool for simulation of delamination initiation and propagation in fibre metal laminates (FMLs) under low‐cycle fatigue loading regime. In the previous studies, the significance of plasticity in delamination growth and modelling of FMLs was not considered. Hence, cohesive zone law that combines the damage evolution with plasticity is developed. The new fatigue damage model is implemented as user‐written subroutines that links with ansys based on the cohesive finite element method. The cohesive zone model constitutive law has been verified by modelling of the delaminated adhesively bonded aluminium joint under normal and shear loadings and compared with the available results in the literature. The developed procedure and tool have been used for the analyses of DCB and ENF specimens under uniform and variable loadings. The obtained results for progressive damage and delamination and stress–strain curves are discussed in this paper.     

Numerical simulation of fatigue crack propagation in compressive residual stress fields of notched round bars

In this paper, the influence of the residual compressive stresses induced by roller burnishing on fatigue crack propagation in the fillet of notched round bar is investigated. A 3D finite element simulation model of rolling has allowed to introduce a residual stress profile as an initial condition. After the rolling process, fatigue loading has been applied to three‐point bending specimens in which an initial crack has been introduced. A numerical predictive method of crack propagation in roller burnished specimens has also been implemented. It is based on a step‐by‐step process of stress intensity factor calculations by elastic finite element analyses. These stress intensity factor results are combined with the Paris law to estimate the fatigue crack growth rate. In the case of roller burnished specimens, a numerical modification concerning experimental crack closure has to be considered. This method is applied to three specimens: without roller burnishing, and with two levels of roller burnishing (type A and type B). In all these cases, the computational finite element predictions of fatigue crack growth rate agree well with the experimental measurements. The developed model can be easily extended to crankshafts in real operating conditions.     

Extended finite element method for cohesive crack growth

The extended finite element method allows one to model displacement discontinuities which do not conform to interelement surfaces. This method is applied to modeling growth of arbitrary cohesive cracks. The growth of the cohesive zone is governed by requiring the stress intensity factors at the tip of the cohesive zone to vanish. This energetic approach avoids the evaluation of stresses at the mathematical tip of the crack. The effectiveness of the proposed approach is demonstrated by simulations of cohesive crack growth in concrete.     

Discrete cohesive zone model for mixed-mode fracture using finite element analysis

The discrete cohesive zone model (DCZM) is implemented using the finite element (FE) method to simulate fracture initiation and subsequent growth when material non-linear effects are significant. Different from the widely used continuum cohesive zone model (CCZM) where the cohesive zone model is implemented within continuum type elements and the cohesive law is applied at each integral point, DCZM uses rod type elements and applies the cohesive law as the rod internal force vs. nodal separation (or rod elongation). These rod elements have the provision of being represented as spring type elements and this is what is considered in the present paper. A series of 1D interface elements was placed between node pairs along the intended fracture path to simulate fracture initiation and growth. Dummy nodes were introduced within the interface element to extract information regarding the mesh size and the crack path orientation. To illustrate the DCZM, three popular fracture test configurations were examined. For pure mode I, the double cantilever beam configuration, using both uniform and biased meshes were analyzed and the results show that the DCZM is not sensitive to the mesh size. Results also show that DCZM is not sensitive to the loading increment, either. Next, the end notched flexure for pure mode II and, the mixed-mode bending were studied to further investigate the approach. No convergence difficulty was encountered during the crack growth analyses. Therefore, the proposed DCZM approach is a simple but promising tool in analyzing very general two-dimensional crack growth problems. This approach has been implemented in the commercial FEA software ABAQUS^® using a user defined subroutine and should be very useful in performing structural integrity analysis of cracked structures by engineers using ABAQUS^®.     

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.     

A cohesive model for fatigue failure of polymers

A cohesive failure model is proposed to simulate fatigue crack propagation in polymeric materials. The model relies on the combination of a bi-linear cohesive failure law used for fracture simulations under monotonic loading and an evolution law relating the cohesive stiffness, the rate of crack opening displacement and the number of cycles since the onset of failure. The fatigue component of the cohesive model involves two parameters that can be readily calibrated based on the classical log-log Paris failure curve between the crack advance per cycle and the range of applied stress intensity factor. The paper also summarizes a semi-implicit implementation of the cohesive model into a cohesive-volumetric finite element framework, allowing for the simulation of a wide range of fatigue fracture problems.     

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

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

Numerical simulation of plasticity induced fatigue crack opening and closure for autofrettaged intersecting holes

The autofrettage of intersecting holes leads to extremely high compressive residual stress fields. These stresses in combination with the plastic deformations decelerate fatigue cracks initiated at the hole intersection notch. Simulations of plasticity induced crack closure of such cracks are presented based on the strip yield and a finite element model. The strip yield model has been extended to allow for an input of residual stresses coming from elsewhere, e.g. from a finite element calculation or measurements. The calculations are applied for constant as well as variable amplitude loading. The numerical expense of the finite element based modelling for variable amplitude loading is still too high if millions of cycles have to be considered. Therefore, a new approximation method is proposed introducing compensatory load sequences. Simulation results are compared to experimentally determined results showing good agreement. However, the accuracy of crack initiation life estimates has turned out to provide a high potential for further improvement.     

Cohesive-zone-model formulation and implementation using the symmetric Galerkin boundary element method for homogeneous solids

A new symmetric boundary integral formulation for cohesive cracks growing in the interior of homogeneous linear elastic isotropic media with a known crack path is developed and implemented in a numerical code. A crack path can be known due to some symmetry implications or the presence of a weak or bonded surface between two solids. The use of a two-dimensional exponential cohesive law and of a special technique for its inclusion in the symmetric Galerkin boundary element method allows us to develop a simple and efficient formulation and implementation of a cohesive zone model. This formulation is dependent on only one variable in the cohesive zone (relative displacement). The corresponding constitutive cohesive equations present a softening branch which induces to the problem a potential instability. The development and implementation of a suitable solution algorithm capable of following the growth of the cohesive zone and subsequent crack growth becomes an important issue. An arc-length control combined with a Newton–Raphson algorithm for iterative solution of nonlinear equations is developed. The boundary element method is very attractive for modeling cohesive crack problems as all nonlinearities are located along the boundaries (including the crack boundaries) of linear elastic domains. A Galerkin approximation scheme, applied to a suitable symmetric integral formulation, ensures an easy treatment of cracks in homogeneous media and excellent convergence behavior of the numerical solution. Numerical results for the wedge split and mixed-mode flexure tests are presented.     

Viscoelastic-damage interface model formulation with friction to simulate the delamination growth in mode II shear

In this paper a formulation of a viscoelastic-damage interface model with friction in mode-II is presented. The cohesive constitutive law contains elastic and damage regimes. It has been assumed that the shear stress in the elastic regime follows the viscoelastic properties of the matrix material. The three element Voigt model has been used for the formulation of relaxation modulus of the material. Damage evolution proceeds according to the bilinear cohesive constitutive law combined with friction stress consideration. Combination of damage and friction is based on the presumption that the damaged area, related to an integration point, can be dismembered into the un-cracked area with the cohesive damage and cracked area with friction. Samples of a one element model have been presented to see the effect of parameters on the cohesive constitutive law. A comparison between the predicted results with available results of end-notched flexure specimens in the literature is also presented to verify the model. Transverse crack tension specimens are also simulated for different applied displacement velocities.     

Failure mode of laser welds in lap‐shear specimens of high strength low alloy (HSLA) steel sheets

In this paper, the failure mode of laser welds in lap‐shear specimens of non‐galvanized SAE J2340 300Y high strength low alloy steel sheets under quasi‐static loading conditions is examined based on experimental observations and finite element analyses. Laser welded lap‐shear specimens with reduced cross sections were made. Optical micrographs of the cross sections of the welds in the specimens before and after tests are examined to understand the microstructure and failure mode of the welds. Micro‐hardness tests were also conducted to provide an assessment of the mechanical properties in the base metal, heat‐affected and fusion zones. The micrographs indicate that the weld failure appears to be initiated from the base metal near the boundary of the base metal and the heat‐affected zone at a distance away from the pre‐existing crack tip, and the specimens fail due to the necking/shear of the lower left load carrying sheets. Finite element analyses based on non‐homogenous multi‐zone material models were conducted to model the ductile necking/shear failure and to obtain the J integral solutions for the pre‐existing cracks. The results of the finite element analyses are used to explain the ductile failure initiation sites and the necking/shear of the lower left load carrying sheets. The J integral solutions obtained from the finite element analyses based on the 3‐zone finite element model indicate that the J integral for the pre‐existing cracks at the failure loads are low compared to the fracture toughness and the specimens should fail in a plastic collapse or necking/shear mode. The effects of the sheet thickness on the failure mode were then investigated for laser welds with a fixed ratio of the weld width to the thickness. For the given non‐homogenous material model, the J integral solutions appear to be scaled by the sheet thickness. With consideration of the plastic collapse failure mode and fracture initiation failure mode, a critical thickness can be obtained for the transition of the plastic collapse or necking/shear failure mode to the fracture initiation failure mode. Finally, the failure load is expressed as a function of the sheet thickness according to the governing equations based on the two failure modes. The results demonstrate that the failure mode of welds of thin sheets depends on the sheet thickness, ductility of the base metal and fracture toughness of the heat‐affected zone. Therefore, failure criteria based on either the plastic collapse failure mode or the fracture initiation failure mode should be used cautiously for welds of thin sheets.