Characteristic crack-tip distances in fracture criteria: Is crack propagation discontinuous?

The present paper describes a possible mechanism for discontinuous crack advance in which surface separation occurs initially not at the crack-tip itself but within the crack-tip plastic zone of size r_p, at the mid-point of the crack-tip characteristic distance d (identified here with the finite growth step Δa), i.e., at the region of maximum opening tensile stress, spreading towards (and also away from) the crack-tip. The crack extension occurs when the crack-tip is reached and full opening over the distance d is completed.Finite element analyses show that this mechanism causes the formation of a rippled crack face surface in elastic-plastic materials in which irreversible plastic deformations take place during each growth step, in sharp contrast with the smooth surface created in ideal elastic materials in which all deformations are fully reversible. Some pictorial evidence of void formation ahead of the crack tip and of ripples during propagation, found in the literature, is presented.Although the present analysis is from a continuum standpoint it is acknowledged that micro structural features and mechanisms can condition the fracture events taking place in the process zone.The implication to the brittle-ductile transition of the dependence of the energy release rate, G^ΔΞ, on the ratio q (=Δa/r_p) is also discussed.     

Fatigue properties and fatigue fracture mechanisms of SiC whiskers or SiC particulate-reinforced aluminium composites

Fatigue properties and fracture mechanisms were examined for three commercially fabricated aluminium matrix composites containing SiC whiskers (SiC_w) and SiC particles (SiC_p) using a rotating bending test. The fatigue strengths were over 60% higher for SiC_w/A2024 composites than that for the unreinforced rolled material, while for the SiC_p/A357 composites, fatigue strengths were also higher than that for the unreinforced reference material. For the SiC_p/A356 composites at a volume fraction of 20%, the fatigue strength was slightly higher than that of the unreinforced material. Fractography revealed that the Mode I fatigue crack was initiated by the Stage I mechanism for the SiC_w/A2024 and SiC_p/A357 composites, while for the SiC_p/A356 composite, the fatigue crack initiated at the voids situated beneath the specimen surfaces. On the other hand, the fatigue crack propagated to the whisker/matrix interface following the formation of dimple patterns or the formation of striation patterns for SiC_w/A2024 composites, while for the SiC_p/A356 and SiC_p/A357 composites the fatigue crack propagated in the matrix near the crack origin and striation patterns were found. Near final failure, dimple patterns, initiated at silicon carbide particles, were frequently observed. Mode I fatigue crack initiation and propagation models were proposed for discontinuous fibre-reinforced aluminium composites. It is suggested that the silicon carbide whiskers or particles would have a very significant effect on the fatigue crack initiation and crack propagation near the fatigue limit.     

Promotion mechanism of electric current on SiCp/Al composite material deformation

To solve the problem of poor formability of SiC_p/Al composite material, enhance the efficiency and reduce the consumption, current is adopted in hot forming process. Under the action of current, the plastic forming performance has made a huge improvement. It is verified that the detour flow and concentration of electrons may lead to a local liquid layer. In addition, the current flows through the crack tip can relax the stress concentration and prevent the propagation of crack. Otherwise, drifting electrons would collide with dislocation which generates electron wind force to make the dislocations arrange along the current direction. Therefore the electric current can enhance the formability of SiC_p/Al composite material and improve the quality of forming part.     

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

Cohesive modeling of dynamic fracture in functionally graded materials

The dynamic fracture of functionally graded materials (FGMs) is modeled using an explicit cohesive volumetric finite element scheme that incorporates spatially varying constitutive and failure properties. The cohesive element response is described by a rate-independent bilinear cohesive failure model between the cohesive traction acting along the cohesive zone and the associated crack opening displacement. A detailed convergence analysis is conducted to quantify the effect of the material gradient on the ability of the numerical scheme to capture elastodynamic wave propagation. To validate the numerical scheme, we simulate dynamic fracture experiments performed on model FGM compact tension specimens made of a polyester resin with varying amounts of plasticizer. The cohesive finite element scheme is then used in a parametric study of mode I dynamic failure of a Ti/TiB FGM, with special emphasis on the effect of the material gradient on the initiation, propagation and arrest of the crack.     

Modeling end notched flexure tests to establish cohesive element Mode II fracture parameters

The objective of this work was to establish Mode II fracture parameters for cohesive elements that can be further utilized to evaluate Mode II interfacial fracture strength of polyurea/AISI 4340 steel composite structures. To obtain the fracture parameters, end notched flexure (ENF) tests were conducted to validate proposed finite element models. The fracture behavior observed from the tests was highly nonlinear and large plastic deformations were involved during crack formation and propagation. A strain incompatibility model was introduced to describe the nonlinearity prior to fracture. This nonlinear and plastic behavior made Linear Elastic Fracture Mechanics (LEFM) approaches not applicable to approximate the fracture parameters. As a result of these experimental observations, finite element analyses of the ENF tests were performed to develop the necessary fracture parameters for cohesive elements selected to replicate the failure modes. Good agreement between the selected numerical models and experimental data was observed.     

Numerical simulation of fatigue crack growth behavior by crack-tip blunting

In ductile metals one of basic mechanisms for fatigue crack growth is that based on crack-tip blunting under the maximum load and re-sharpening of the crack-tip under minimum load. In this paper, simulations of fatigue crack growth by crack-tip blunting using ANSYS finite element code are presented. This investigation focuses solely on simulation of fatigue crack growth due to crack-tip plasticity only. As such, any material damage and its fracture is not considered. Due to high plastic deformation the present simulations utilize a remeshing technique which allows applying a number of load cycles without terminating the simulation due to the error caused by excessive mesh distortion. The simulations were conducted using a center cracked specimen under various loading conditions including different load ranges and load ratios R = −1, 0 and 0.333. It is shown that fatigue crack growth (FCG) slows down with number of cycles towards a steady state value. The simulated FCG data for constant amplitude loading follow the Paris power law relationship and also indicate a typical R-ratio dependence. It can be noted that for all load cases with load ratios R > 0 no crack closure in the vicinity of the crack-tip wake was observed.     

Analysis of mode I delamination of z-pinned composites using a non-dimensional analytical model

We present a non-dimensional analytical model for crack propagation in a z-pinned double cantilever beam specimen (DCB) under mode I loading. Effect of various design parameters on the crack bridging length and apparent fracture toughness are investigated using this model. The efficacy of the analytical model is evaluated by comparing the results with 3D finite element (FE) simulations of the DCB. In the FE model the z-pins are modeled as discrete nonlinear elements. Bi-linear cohesive elements are used ahead of the crack tip to account for the interlaminar fracture toughness of the composite material. The results for load–deflection and crack length obtained from the analytical model and the FE model are compared and found to be in good agreement. The proposed non-dimensional analytical model will be useful in the design and analysis of translaminar reinforcements for composite structures.     

A full‐discontinuous Galerkin formulation of nonlinear Kirchhoff–Love shells: elasto‐plastic finite deformations,parallel computation,and fracture applications

Because of its ability to take into account discontinuities, the discontinuous Galerkin (DG) method presents some advantages for modeling cracks initiation and propagation. This concept has been recently applied to three‐dimensional simulations and to elastic thin bodies. In this last case, the assumption of small elastic deformations before cracks initiation or propagation reduces drastically the applicability of the framework to a reduced number of materials. To remove this limitation, a full‐DG formulation of nonlinear Kirchhoff–Love shells is presented and is used in combination with an elasto‐plastic finite deformations model. The results obtained by this new formulation are in agreement with other continuum elasto‐plastic shell formulations. Then, this full‐DG formulation of Kirchhoff–Love shells is coupled with the cohesive zone model to perform thin body fracture simulations. As this method considers elasto‐plastic constitutive laws in combination with the cohesive model, accurate results compared with the experiments are found. In particular, the crack path and propagation rate of a blasted cylinder are shown to match experimental results. One of the main advantages of this framework is its ability to run in parallel with a high speed‐up factor, allowing the simulation of ultra fine meshes. Copyright © 2012 John Wiley & Sons, Ltd.     

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.     

Automated dynamic fracture procedure for modelling mixed-mode crack propagation using explicit time integration brick finite elements

Several fracture codes have been developed in recent years to perform analyses of dynamic crack propagation in arbitrary directions. However, general-purpose, commercial finite-element software which have capabilities to do fracture analyses are still limited in their use to stationary cracks and crack propagation along trajectories known a priori . In this paper, we present an automated fracture procedure implemented in the large-scale, nonlinear, explicit, finite-element code DYNA3D which can be used to simulate dynamic crack propagation in arbitrary directions. The model can be used to perform both generation- and application-phase simulations of self-similar as well as non-self-similar dynamic crack propagation in linear elastic structures without user intervention. It is developed based on dynamic fracture mechanics concepts and implemented for three-dimensional solid elements. Energy approach is used in the model to check for crack initiation/propagation. Dynamic energy release rate and stress intensity factors are determined from far-field finite-element field solutions using finite-domain integrals. Fracture toughness is input as a function of crack-tip velocity, and when the criterion for crack growth is satisfied, an element deletion-and-replacement re-meshing procedure is used along with a gradual nodal release technique to update the crack geometry and model the crack propagation. Direction of crack propagation is determined using the maximum circumferential stress criterion. Numerical simulations of experiments involving non-self-similar crack propagation are performed, and results are presented as verification examples.     

Effects of particle size on fatigue crack initiation and small crack growth in SiC particulate-reinforced aluminium alloy composites

Fatigue crack initiation and small crack growth were studied under axial loading using powder metallurgy 2024 aluminum-matrix composites reinforced with SiC particles of three different sizes of 5, 20 and 60 μm. The 5 and 20 μm SiC_p/Al composites exhibited nearly the same fatigue strength as the unreinforced alloy, while the 60 μm SiC_p/Al composite showed a significantly lower fatigue strength due to its inferior crack initiation resistance that could be attributed to interface debonding between particles and the matrix. Small crack growth behaviour was different depending on stress level. At a low applied stress, the addition of SiC particles enhanced the growth resistance, particularly in the composites reinforced with coarser particles, while at a high applied stress, the 60 μm SiC_p/Al composite showed a considerably low growth resistance, which could be attributed to interaction and coalescence of multiple cracks. In the 5 μm SiC_p/Al composite, small cracks grew avoiding particles and thus few particles appearing on the fracture surfaces were seen, particularly in small crack size region. In the 20 and 60 μm SiC_p/Al composites, they grew along interfaces between particles and the matrix and the number of particles appearing on the fracture surfaces increased with increasing crack size or maximum stress intensity factor.     

Impact and dynamic fracture of tough polymers by thermal decohesion in a Dugdale zone

At low temperatures and high loading rates, normally tough crystalline thermoplastics may undergo a transition from ductile tearing to brittle rapid crack propagation (RCP). It is proposed here that RCP — characterised by low toughness, high crack speed (>100 m/s) and a macroscopically smooth fracture surface — occurs by self-sustained melting of a layer, one chain length thick, at each cohesive surface of a crack-tip craze, due to adiabatic heating. Initiation of RCP from a rapidly loaded sharp notch, i.e. impact fracture, requires both the formation of this melt layer, and sufficient crack extension force to propagate it. A schematic linear-elastic analysis based on the Dugdale model accounts both for the measured dynamic fracture resistance, and for the variation of impact fracture resistance with impact speed, in two pipe-grade polyethylenes and in a neat and a rubber-toughened polyacetal. It is concluded that crack initiation resistance, unlike dynamic fracture resistance, cannot be defined as a geometry-dependent material property.     

Fracture analysis of dissimilar Al‐Al friction stir welded joints under tensile/shear loading

In this research, fracture of dissimilar friction stir welded (FSWed) joint made of Al 7075‐T6 and Al 6061‐T6 aluminum alloys is investigated in the cracked semi‐circular bend (CSCB) specimen under mixed mode I/II loading. Due to the elastic‐plastic behavior of the welded material and the existence of significant plastic deformations around the crack tip at the propagation instance, fracture prediction of the FSWed specimens needs some failure criteria in the context of the elastic‐plastic fracture mechanics which are very complicated and time‐consuming. For this purpose, the Equivalent Material Concept (EMC) is used herein by which the tensile behavior of the welded material is equated with that of a virtual brittle material. By combining EMC with the 2 brittle fracture criteria, namely the maximum tangential stress (MTS) and mean stress (MS) criteria, the load‐carrying capacity (LCC) of the FSWed CSCB specimens is predicted. Comparison of the experimental results and theoretical predictions from the 2 criteria showed that both criteria could accurately predict the LCC of the cracked specimens. Moreover, as the contribution of mode II loading increases, the size of the plastic region around the crack tip at failure increases, leading to increasing the LCC.     

Crack growth resistance behavior of a functionally graded material: computational studies

This paper describes crack growth resistance simulation in a ceramic/metal functionally graded material (FGM) using a cohesive zone ahead of the crack front. The plasticity in the background (bulk) material follows J₂ flow theory with the flow properties determined by a volume fraction based, elastic-plastic model (extension of the original Tamura-Tomota-Ozawa model). A phenomenological, cohesive zone model with six material-dependent parameters (the cohesive energy densities and the peak cohesive tractions of the ceramic and metal phases, respectively, and two cohesive gradation parameters) describes the constitutive response of the cohesive zone. Crack growth occurs when the complete separation of the cohesive surfaces takes place. The crack growth resistance of the FGM is characterized by a rising J-integral with crack extension (averaged over the specimen thickness) computed using a domain integral (DI) formulation. The 3-D analyses are performed using WARP3D, a fracture mechanics research finite element code, which incorporates solid elements with graded elastic and plastic properties and interface-cohesive elements coupled with the functionally graded cohesive zone model. The paper describes applications of the cohesive zone model and the DI method to compute the J resistance curves for both single-edge notch bend, SE(B), and single-edge notch tension, SE(T), specimens having properties of a TiB/Ti FGM. The numerical results show that the TiB/Ti FGM exhibits significant crack growth resistance behavior when the crack grows from the ceramic-rich region into the metal-rich region. Under these conditions, the J-integral is generally higher than the cohesive energy density at the crack tip even when the background material response remains linearly elastic, which contrasts with the case for homogeneous materials wherein the J-integral equals the cohesive energy density for a quasi-statically growing crack.     

Experimental and numerical investigations of mixed mode crack growth resistance of a ductile adhesive joint

Polymeric adhesive joints are extensively employed in various industrial and technological applications. It has been observed that in ductile adhesive joints, interface fracture is a common mode of failure which may involve stable crack propagation followed by catastrophic growth. The objectives of this paper are to investigate the effects of bondline thickness and mode mixity on the steady state energy release rate J_ss of such a joint. To this end, a combined experimental and numerical investigation of interfacial crack growth is carried out using a modified compact tension shear specimen involving two aluminium plates bonded by a thin ductile adhesive layer. A cohesive zone model along with a simple traction versus separation law is employed in the finite element simulations of crack growth. It is observed that J_ss increases strongly as mode II loading is approached. Also, it enhances with bondline thickness in the above limit. These trends are rationalized by examining the plastic zones obtained from the numerical simulations. The numerically generated J_ss values are found to agree well with the corresponding experimental results.     

Cohesive fracture modeling of elastic-plastic crack growth in functionally graded materials

This work investigates elastic-plastic crack growth in ceramic/metal functionally graded materials (FGMs). The study employs a phenomenological, cohesive zone model proposed by the authors and simulates crack growth by the gradual degradation of cohesive surfaces ahead of the crack front. The cohesive zone model uses six material-dependent parameters (the cohesive energy densities and the peak cohesive tractions of the ceramic and metal phases, respectively, and two cohesive gradation parameters) to describe the constitutive response of the material in the cohesive zone. A volume fraction based, elastic-plastic model (extension of the original Tamura-Tomota-Ozawa model) describes the elastic-plastic response of the bulk background material. The numerical analyses are performed using WARP3D, a fracture mechanics research finite element code, which incorporates solid elements with graded elastic and plastic properties and interface-cohesive elements coupled with the functionally graded cohesive zone model. Numerical values of volume fractions for the constituents specified at nodes of the finite element model set the spatial gradation of material properties with isoparametric interpolations inside interface elements and background solid elements to define pointwise material property values. The paper describes applications of the cohesive zone model and the computational scheme to analyze crack growth in a single-edge notch bend, SE(B), specimen made of a TiB/Ti FGM. Cohesive parameters are calibrated using the experimentally measured load versus average crack extension (across the thickness) responses of both Ti metal and TiB/Ti FGM SE(B) specimens. The numerical results show that with the calibrated cohesive gradation parameters for the TiB/Ti system, the load to cause crack extension in the FGM is much smaller than that for the metal. However, the crack initiation load for the TiB/Ti FGM with reduced cohesive gradation parameters (which may be achieved under different manufacturing conditions) could compare to that for the metal. Crack growth responses vary strongly with values of the exponent describing the volume fraction profile for the metal. The investigation also shows significant crack tunneling in the Ti metal SE(B) specimen. For the TiB/Ti FGM system, however, crack tunneling is pronounced only for a metal-rich specimen with relatively smaller cohesive gradation parameter for the metal.     

An embedded atom hyperelastic constitutive model and multiscale cohesive finite element method

Based on embedded atom method (EAM), an embedded atom hyperelastic (EAH) constitutive model is developed. The proposed EAH constitutive model provides a multiscale formalism to determine mesoscale or macroscale material behavior by atomistic information. By combining the EAH with cohesive zone model (CZM), a multiscale embedded atom cohesive finite element model (EA-cohesive FEM) is developed for simulating failure of materials at mesoscale and macroscale, e.g. fracture and crack propagation etc. Based on EAH, the EA-cohesive FEM applies the Cauchy-Born rule to calculate mesoscale or macroscale material response for bulk elements. Within the cohesive zone, a generalized Cauchy-Born rule is applied to find the effective normal and tangential traction-separation cohesive laws of EAH material. Since the EAM is a realistic semi-empirical interatomic potential formalism, the EAH constitutive model and the EA-cohesive FEM are physically meaningful when it is compared with experimental data. The proposed EA-cohesive FEM is validated by comparing the simulation results with the results of large scale molecular dynamics simulation. Simulation result of dynamic crack propagation is presented to demonstrate the capacity of EA-cohesive FEM in capturing the dynamic fracture.     

Residual strength prediction of a complex structure using crack extension analyses

The residual strength of a flat panel (thickness 7.6 mm) with five stringers, machined from a monolithic block of Al2024-T351 material, which contained a crack that divided the central stringer, was to be predicted during a Round Robin organised by ASTM. The initial crack tips were right ahead of the stringers #2 and #4, respectively, so that crack branching along the skin and into the stringers occurred after initiation. The prediction has been achieved using finite element simulations including crack extension, for which a cohesive model was utilised. Conventional material properties, yield and ultimate strength as well as experimental results from M(T) specimens in terms of force, COD and Δa, were given. The residual strength prediction was performed in two-steps: First the crack extension parameters for the cohesive model, the cohesive strength, T₀, and the cohesive energy, Γ₀, were determined by numerical reproduction of the results of the M(T) specimen. With the optimised parameters, the five-stringer panel was modelled. These steps were conducted by two different finite element models: by a shell and a 3D finite element mesh. It turned out that it is possible to analyse the structure with both models. In the 3D case, the residual strength prediction was conservative and the deviation of the predicted from the experimental value was below 9%. The results of the shell simulation were even closer to the experiment (deviation approximately 3%), but the simulation was non-conservative.     

Dynamic effects of inclusions and microcracks on a main crack

The shielding and amplification effects of multiple inclusions and microcracks on the tip or the growth path of a main crack under dynamic loading are investigated using numerical simulations. The simulations employ a combined numerical tool based on the time-domain boundary element method together with the sub-region technique for multi-regions, the maximum circumferential stress criterion for crack-growth direction and discrete modelling for the crack propagation. New elements of constant length are added to the moving crack-tip to simulate the growth at an adaptable time-step according to the crack propagation criterion. Of particular interest is the study of the effects of the sizes and positions of inclusions and microcracks and the material combinations on the dynamic stress intensity factors and the crack growth path. The numerical results demonstrate the crack-tip shielding and amplification effects of inclusions and micro-cracks.