Dynamic ductile fracture in aluminum round bars: experiments and simulations

The application of rate-dependent cohesive elements is validated in simulation of ductile fracture in aluminum round bars under dynamic loading conditions. Smooth and notched round bars made of AA6060-T6 are tested and simulated under quasi-static and dynamic loadings. The smooth round bar is modeled using finite elements that obey Gurson–Tvergaard–Needleman (GTN) formulation as the constitutive equation. Comparing with experimental results, corresponding GTN parameters and rate-dependent plasticity of the alloy are obtained. A single strain rate-dependent GTN element with the obtained parameters is examined under different values of stress triaxiality and loading rates. The resulting stress-elongation curves represent the traction separation law (TSL) for cohesive elements and the variations of the maximum traction and the energy absorbed are investigated. The notched round bars are modeled by axisymmetric continuum and cohesive elements. The undamaged bulk material is elastic-visco plastic and the cohesive elements obey the TSL defined from the single element calculations. The experiments are simulated by these models in which the cohesive elements are rate sensitive and automatically obtain the values of the total strain rate from their adjacent continuum elements to update the values of the cohesive strength during the analysis. The results of the analysis, including maximum load, time of failure and diameter reduction are validated with the experimental results. The effects of element size, rate-dependent plasticity of the material and stress triaxiality are also discussed.     

Three-dimensional modeling of ductile crack growth: Cohesive zone parameters and crack tip triaxiality

For 10 mm thick smooth-sided compact tension specimens made of a pressure vessel steel 20MnMoNi55, the interrelations between the cohesive zone parameters (the cohesive strength, T_max, and the separation energy, Γ) and the crack tip triaxiality are investigated. The slant shear-lip fracture near the side-surfaces is modeled as a normal fracture along the symmetry plane of the specimen. The cohesive zone parameters are determined by fitting the simulated crack extensions to the experimental data of a multi-specimen test. It is found that for constant cohesive zone parameters, the simulated crack extension curves show a strong tunneling effect. For a good fit between simulated and experimental crack growth, both the cohesive strength and the separation energy near the side-surface should be considerably lower than near the midsection. When the same cohesive zone parameters are applied to the 3D model and a plane strain model, the stress triaxiality in the midsection of the 3D model is much lower, the von-Mises equivalent stress is distinctly higher, and the crack growth rate is significantly lower than in the plane strain model. Therefore, the specimen must be considered as a thin specimen. The stress triaxiality varies dramatically during the initial stages of crack growth, but varies only smoothly during the subsequent stable crack growth. In the midsection region, the decrease of the cohesive strength results in a decrease of the stress triaxiality, while the decrease of the separation energy results in an increase of the triaxiality.     

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.     

Accurate simulation of delamination growth under mixed-mode loading using cohesive elements: Definition of interlaminar strengths and elastic stiffness

A methodology for predicting accurately the propagation of delamination under mixed-mode fracture with cohesive elements is proposed. It is shown that changes in the local mode ratio during the evolution of damage under mixed-mode loading can cause errors in the determination of the energy dissipation and result in inaccurate predictions of the global load–displacement response – even under conditions where, according to Linear Elastic Fracture Mechanics, the global mode ratio is constant. To address this difficulty, relations between the interlaminar strengths and the penalty stiffness are proposed which ensure a correct energy dissipation when delamination propagates. The validity of the proposed methodology is demonstrated for different mode ratios by comparison with the corresponding analytical solutions.     

Effect of strength mis-match and dynamic loading on ductile fracture initiation

It has been well known that ductile fracture of steels is accelerated by triaxial stresses. The characteristics of ductile crack initiation in steels are evaluated quantitatively using a two-parameter criterion based on equivalent plastic strain and stress triaxiality.The present study focuses on the effects of geometrical discontinuity, strength mis-match, which can elevate plastic constraint due to heterogeneous plastic straining, and loading rate on the critical condition for ductile fracture initiation using a two-parameter criterion. Fracture initiation testing has been conducted under static and dynamic loading using circumferentially notched round-bar specimens. In order to evaluate the stress/strain state in the specimens, especially under dynamic loading, a thermal elastic-plastic dynamic finite element (FE) analysis considering the temperature rise due to plastic deformation has been carried out.The tensile tests on specimens with an undermatching interlayer showed that the relationship between the critical equivalent plastic strain to initiate ductile fracture and stress triaxiality was equivalent to that obtained on homogeneous specimens under static loading. Moreover, the two-parameter criterion for ductile fracture initiation is shown to be independent of the loading rate. It was demonstrated that the critical global strain to initiate ductile fracture in specimens with strength mis-match under various loading rate can be estimated based on the local criterion, that is two-parameter criterion obtained on homogeneous specimens under static tension, by mean of FE-analysis taken into account accurately both strength mis-match and dynamic loading effects on stress/strain behaviors.     

Simulation of ductile crack growth in thin aluminum panels using 3-D surface cohesive elements

This work describes the formulation and application of a 3-D, interface-cohesive finite element model to predict quasi-static, ductile crack extension in thin aluminum panels for mode I loading and growth. The fracture model comprises an initially zero thickness, interface element with constitutive response described by a nonlinear traction-separation relationship. Conventional volumetric finite elements model the nonlinear (elastic-plastic) response of background (bulk) material. The interface-cohesive elements undergo gradual decohesion between faces of the volumetric elements to create new traction free crack faces. The paper describes applications of the computational model to simulate crack extension in C(T) and M(T) panels made of a 2.3 mm thick, Al 2024-T3 alloy tested as part of the NASA-Langley Aging Aircraft program. Parameters of the cohesive fracture model (peak opening traction and local work of separation) are calibrated using measured load vs. outside surface crack extensions of high constraint (T-stress > 0) C(T) specimens. Analyses of low constraint M(T) specimens, having widths of 300 and 600 mm and various a/W ratios, demonstrate the capabilities of the calibrated model to predict measured loads and outside surface crack extensions. The models capture accurately the strong 3-D effects leading to various degrees of crack front tunneling in the C(T) and M(T) specimens. The predicted crack growth response shows rapid convergence with through-thickness mesh refinement. Adaptive load increment procedures to control the rate of decohesion in the interface elements leads to stable, rapidly converging iterations in the globally implicit solution procedures.     

An irreversible cohesive zone model for interface fatigue crack growth simulation

Fatigue crack growth (FCG) along an interface is studied. Instead of using the Paris equation, the actual process of material separation during FCG is described by the use of an irreversible constitutive equation for the cyclic interface traction-separation behavior within the cohesive zone model (CZM) approach. In contrast to past development of CZMs, the traction-separation behavior does not follows a predefined path. The model definition, its predicted cyclic material separation behavior and application to a numerical study of interface FCG in double-cantilever beam, end-loaded split and mixed-mode beam specimens are reported.     

A numerical analysis of constraint effects in fatigue crack growth by use of an irreversible cohesive zone model

The analysis of constraint effects in fatigue crack growth in multi-layer structures is discussed. The process of material separation under cyclic loading is described by a cohesive zone model (CZM) with an irreversible constitutive relationship. The traction–separation behavior does not follow a predefined path, but is dependent on the evolution of the damage dependent cohesive zone properties. A modified boundary layer model is used in simulations of fatigue crack growth along the centerline crack of the metal layer sandwiched between two elastic substrates. Fatigue crack growth is computed for a series of values of metal layer thickness under constant and variable amplitude loading conditions. The results of the computations demonstrate that certain combinations of load magnitude, layer thickness and material properties results in significant constrain effects in fatigue crack growth. The influence of these constraint effects on fatigue crack growth rates and on crack closure processes is determined. The evolutions of the traction–separation law, the accumulated and current plastic zones, as well as the stress fields during the crack propagation are discussed.     

Constraint effect on the ductile crack growth resistance of circumferentially cracked pipes

A systematic study has been carried out by using 2D axisymmetric models to understand the ductile fracture behavior of pipes with internal and external circumferential cracks. Crack growth resistance curves have been computed using the complete Gurson model. Pipes with various diameter-to-thickness ratios, internal pressure, crack depths and material properties are analyzed. The results have been compared with those of corresponding SENT and standard SENB specimens. It clearly indicates that the SENT specimen is a good representation of circumferentially flawed pipes and an alternative to the conventional standard SENB specimen for the fracture mechanics testing in engineering critical assessment of pipes.     

Modeling of crack tip high inertia zone in dynamic brittle fracture

A phenomenological modification is proposed to the existing cohesive constitutive law of Roy and Dodds to model the crack tip high inertia region proposed by Gao. The modification involves addition of a term which is attributed to fracture mechanisms that result in high energy dissipation around the crack tip. This term is assumed to be a function of external energy per unit volume input into the system. Finite element analysis is performed on PMMA with constant velocity boundary conditions and mesh discretizations based on the work of Xu and Needleman. The cohesive model with the proposed dissipative term is only applied in the high inertia zone and the traditional Roy and Dodds model is applied on cohesive elements in the rest of the domain. The results show that crack propagates in three phases with a speed of 0.35c_R before branching, confirming experimental observations. The modeling of high inertia zone is one of the key aspects to understanding brittle fracture.     

Effect of residual stresses on ductile crack growth resistance

It is well known that residual stresses influence the ductile fracture behaviour. In this paper, a numerical study was performed to assess the effect of residual stresses on ductile crack growth resistance of a typical pipeline steel. A modified boundary layer model was employed for the analysis under plane strain, Mode I loading condition. The residual stress fields were introduced into the finite element model by the eigenstrain method. A sharp crack was embedded in the center of the weld region. The complete Gurson model has been applied to simulate the ductile fracture by microvoid nucleation, growth and coalescence. Results show that tensile residual stresses can significantly reduce the crack growth resistance when the crack growth is small compared with the length scale of the tensile residual stress field. With the crack growth, the effect of residual stresses on the crack growth resistance tends to diminish. The effect of residual stress on ductile crack growth resistance seems independent of the size of geometrically similar welds. When normalized by the weld zone size, the ductile crack growth resistance collapses into one curve, which can be used to assess the structural integrity and evaluate the effect of residual stresses. It has also been found that the effect of residual stresses on crack growth resistance depends on the initial void volume fraction f₀, hardening exponent n and T-stress.     

Simulation of hydrogen assisted stress corrosion cracking using the cohesive model

This paper investigates the effect of hydrogen diffusion on stable crack propagation by using numerical finite element simulations based on the cohesive model. The model with its two common parameters, cohesive strength, T₀, and critical separation, δ₀, and its two additional parameters for stress corrosion cracking, the effective diffusivity, D_eff, and a material parameter, α, which represents the reduction of the cohesive strength, is described. This model is then employed to predict the stable crack propagation in C(T) specimens made from a high strength structural steel which were tested under hydrogen charging conditions in rising displacement tests using various deformation rates. It is shown that, in general, the prediction of stable crack propagation is promising, but may be further improved by the use of a more sophisticated diffusion equation. Finally, the influence of variations of the effective diffusivity and the cohesive strength reduction on the thus simulated crack growth resistance curves is studied.     

Simulation of ductile crack propagation and determination of CTOAs in pipeline steels using cohesive zone modelling

This paper presents recent results of numerical studies on stable crack extension of high toughness steels typical of those in modern gas pipelines using cohesive zone modelling (CZM). The main focus of the work is on the determination of crack‐tip opening angles (CTOAs) of these steels from CZM. Two sets of materials are modelled. The first material set models a typical structural steel, with variable toughness described by four traction–separation (TS) laws. The second set models an X70 pipe steel, with three different TS laws. For each TS law, there are three defining parameters: the maximum cohesive strength, the final separation and the work of separation. The specimens analysed include a crack in an infinite plate (small‐scale yielding, SSY) and a standard drop‐weight tear test (DWTT). Fracture propagation characteristics and values of CTOA are obtained from these two types of specimens. It is shown that cohesive zone models can be successfully used to simulate ductile crack propagation and to numerically measure CTOAs. The ductile crack propagation characteristics and CTOAs obtained from SSY and DWTT specimens are compared for each set of steels. In addition, the CTOA results from numerical CZM of DWTT specimens of X70 steel are compared with those from laboratory tests.     

Simulation of dynamic crack growth using the generalized interpolation material point (GIMP) method

Dynamic crack growth is simulated by implementing a cohesive zone model in the generalized interpolation material point (GIMP) method. Multiple velocity fields are used in GIMP to enable handling of discrete discontinuity on either side of the interface. Multilevel refinement is adopted in the region around the crack-tip to resolve higher strain gradients. Numerical simulations of crack growth in a homogeneous elastic solid under mode-II plane strain conditions are conducted with the crack propagating along a weak interface. A parametric study is conducted with respect to varying impact speeds ranging from 5 m/s to 60 m/s and cohesive strengths from 4 to 35 MPa. Numerical results are compared qualitatively with the dynamic fracture experiments of Rosakis et al. [(1999) Science 284:1337–1340]. The simulations are capable of handling crack growth with crack-tip velocities in both sub-Rayleigh and intersonic regimes. Crack initiation and propagation are the natural outcome of the simulations incorporating the cohesive zone model. For various impact speeds, the sustained crack-tip velocity falls either in the sub-Rayleigh regime or in the region between (c _S is the shear wave speed) and c _D (c _D is the dilatational wave speed) of the bulk material. The Burridge–Andrews mechanism for transition of the crack-tip velocity from sub-Rayleigh to intersonic speed of the bulk material is observed for impact speeds ranging from 9.5 to 60 m/s (for normal and shear cohesive strengths of 24 MPa). Within the intersonic regime, sustained crack-tip velocities between 1.66 c _S (or 0.82 c _D ) and 1.94 c _S (or 0.95 c _D ) were obtained. For the cases simulated in this work, within the stable intersonic regime, the lowest intersonic crack-tip velocity obtained was 1.66 c _S (or 0.82 c _D ).     

Analysis of a rate-dependent cohesive model for dynamic crack propagation

The effect of including rate-dependence in the cohesive zone modeling of steady-state and transient dynamic crack propagation is analyzed. Spontaneous crack propagation simulations are performed using a spectral form of the elastodynamic boundary integral equations, while the solution to the steady-state problem is obtained by solving the governing Cauchy singular equation on the crack plane. The steady-state analysis shows that the existing techniques for solving the Cauchy singular integral equation are not suitable. A solution technique for the underlying Riemann-Hilbert problem for the chosen rate and damage-dependent cohesive law is presented. Under spontaneous propagation conditions, quasi-steady-state speeds slower than the theoretically predicted shear wave speed are possible. Results also show that, due to the dissipation of energy inside the cohesive zone, the energy required for crack propagation increases with the crack speed.     

Formulation of a cohesive zone model for a Mode III crack

Cohesive zone models have been proven effective in modeling crack initiation and propagation phenomena. In this work, a possible form for a Mode III cohesive zone model is formulated from elastic stress and displacement fields around a crack with a cohesive zone ahead of the crack tip. A traction-separation relation for the model is derived as a direct consequence of the formulation, which establishes some intrinsic connections between properties of the cohesive zone and those of the bulk material. Interestingly, this model states that the von Mises effective stress in the cohesive zone is constant, which may be related to the bulk material’s yield stress and is consistent with the assumption made in conventional strip-yield elastic-plastic solutions.     

Numerical simulation of constraint effects in fatigue crack growth

The computational analysis of constraint effects on fatigue crack growth is discussed. An irreversible cohesive zone model is used in the computations to describe the processes of material separation under cyclic loading. This approach is promising for the investigation of fatigue crack growth under constraint as the energy dissipation due to the formation of new crack surface and cyclic plastic deformation is accounted for independently. Fatigue crack growth in multi-layer structures under consideration of different levels of T-stress are conducted with a modified boundary layer model. Fatigue crack growth is computed as a function of layer thickness and T-stress for constant and variable amplitude loading cases.     

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.     

Simulations of dynamic crack propagation in brittle materials using nodal cohesive forces and continuum damage mechanics in the distinct element code LDEC

Experimental data indicates that the limiting crack speed in brittle materials is less than the Rayleigh wave speed. One reason for this is that dynamic instabilities produce surface roughness and microcracks that branch from the main crack. These processes increase dissipation near the crack tip over a range of crack speeds. When the scale of observation (or mesh resolution) becomes much larger than the typical sizes of these features, effective-medium theories are required to predict the coarse-grained fracture dynamics. Two approaches to modeling these phenomena are described and used in numerical simulations. The first approach is based on cohesive elements that utilize a rate-dependent weakening law for the nodal cohesive forces. The second approach uses a continuum damage model which has a weakening effect that lowers the effective Rayleigh wave speed in the material surrounding the crack tip. Simulations in this paper show that while both models are capable of increasing the energy dissipated during fracture when the mesh size is larger than the process zone size, only the continuum damage model is able to limit the crack speed over a range of applied loads. Numerical simulations of straight-running cracks demonstrate good agreement between the theoretical predictions of the combined models and experimental data on dynamic crack propagation in brittle materials. Simulations that model crack branching are also presented.