Cohesive zone laws for fatigue crack growth: Numerical field projection of the micromechanical damage process in an elasto-plastic medium

Cohesive zone failure models are widely used to simulate fatigue crack propagation under cyclic loading, but the model parameters are phenomenological and are not closely tied to the underlying micromechanics of the problem. In this paper, we will inversely extract the cohesive zone laws for fatigue crack growth in an elasto-plastic ductile solid using a field projection method (FPM), which projects the equivalent tractions and separations at the cohesive crack-tip from field information outside the process zone. In our small-scale yielding model, a single row of discrete voids is deployed directly ahead of a crack in an elasto-plastic medium subjected to cyclic mode I K-field loading. Damage accumulation under cyclic loading is captured by the growth of voids within the micro-voiding zone ahead of the crack, while the evolution of the cohesive zone law representing the micro-voiding zone is inversely extracted via the FPM. We show that the field-projected cohesive zone law captures the essential micromechanisms of fatigue crack growth in the ductile medium: from loading and unloading hysteresis caused by void growth and plastic hardening, to the softening damage locus associated with crack propagation via a void by void growth mechanism. The results demonstrate the effectiveness of the FPM in obtaining a micromechanics-based cohesive zone law in-place of phenomenological models, which opens the way for a unified treatment of fatigue crack problems.     

Cohesive-zone laws for void growth — I. Experimental field projection of crack-tip crazing in glassy polymers

A hybrid framework for inverse analysis of crack-tip cohesive-zone model is developed in this two-part paper to measure cohesive-zone laws of void growth in polymers by combining analytical, experimental, and numerical approaches. This paper focuses on experimental measurements of the cohesive-zone laws for two nonlinear fracture processes in glassy polymers, namely multiple crazing in crack-growth toughening of rubber-toughened high-impact polystyrene (HIPS) and crazing of steady-state crack growth in polymethylmethacrylate (PMMA) under a methanol environment. To this end, electronic speckle pattern interferometry (ESPI) is first applied to measure the crack-tip displacement fields surrounding the fracture process zones in these polymers. These fields are subsequently equilibrium smoothed and used in the extraction of the cohesive-zone laws via an analytical solution method of the inverse problem, the planar field projection method (P-FPM) [Hong, S., Kim, K.-S., 2003. Extraction of cohesive-zone laws from elastic far-fields of a cohesive crack tip: a field projection method. Journal of the Mechanics and Physics of Solids 51, 1267-1286]. Results show that the proposed framework of the P-FPM could provide a systematic way of finding the shape of the cohesive-zone laws governed by the different micro-mechanisms in the fracture processes. In HIPS, inter-particle multiple crazing develops and the craze zone broadens ahead of a crack-tip under mechanical loading. The corresponding cohesive-zone relationship of the multiple-craze zone is found to be highly convex, which indicates effectiveness of rubber particle toughening. It is also observed that the effective peak traction, 7 MPa, in the crack-tip cohesive zone of HIPS (30% rubber content) is lower than the uniaxial yield stress of 9 MPa, presumably due to stress multi-axiality effects. In contrast, in PMMA, methanol localizes the crack-tip craze, weakening the craze traction for craze-void initiation to about 9 MPa and the fibril pull-out stress to less than 6 MPa. This reduction in cohesive traction, coupled with a strongly concave traction-separation cohesive-zone relationship, signifies environmental embrittlement of PMMA. These experimentally determined cohesive-zone laws are compared with detailed numerical analyses of effective microscale-void growth ahead of a crack tip in Part II.     

Void growth and coalescence in ductile solids with stage III and stage IV strain hardening

State of the art ductile fracture models often rely on simple power laws to describe the strain hardening of the matrix material. Power laws do not distinguish between the two main stages of hardening observed in polycrystals, referred to as stage III and stage IV hardening, and which emerge from the evolution of the dislocation substructure. The aim of this study is to couple a physics based strain hardening law including these two stages to a micromechanics based ductile damage model. One of the main motivations is that, the stage IV constant hardening rate stage, occurring only at large strain, will be attained in most ductile failure problems if not at the overall level of deformation, at least locally around the growing voids. Furthermore, proper modelling of the stage III involving dislocation storage and recovery terms and the transition to stage IV provides a link with the underlying physical mechanisms of deformation and with the microstructure. First, in order to evaluate the effects of the stage III and stage IV hardening on void growth and coalescence, an extensive parametric study is performed on two-dimensional (2D) axisymmetric finite element (FE) unit cell calculations, using a Kocks-Mecking type hardening law. The cell calculations demonstrate that accounting for the stage IV hardening can have a profound effect on delaying void coalescence and increasing the ductility. The magnitude of the recovery term during stage III has also a significant effect on the void growth rate. Then, the Kocks-Mecking law is incorporated into the Gologanu-Leblond-Devaux (GLD) porous plasticity model supplemented by two different versions of the Thomason void coalescence criterion. The predictions of the damage model are in good agreement with the results of the FE calculations in terms of the stress-strain curves, the evolution of void shape and porosity, as well as the strain value at the onset of void coalescence.     

The size dependence of micro-toughness in ductile fracture

Micro-toughness in ductile fracture is defined as the plastic work dissipated per unit fracture surface area in the material separation processes of void growth and coalescence. A micromechanics model for the estimation of the size dependence of micro-toughness in ductile fracture is presented. Size effects are incorporated in the model using the conventional mechanism-based strain gradient plasticity (CMSG) theory. A finite element model of an axisymmetric representative unit cell with an initial spherical void is used to validate model predictions. Two characteristic length scales emerge from the model. The initial void radius sets the scale for the initial spherical void growth. For the subsequent void coalescence, the scale is set by the width of the intervoid ligament. Energy dissipation in ductile fracture is found to be dominated by the mechanisms of coalescence, and the micro-toughness in ductile fracture is found to be size dependent for dimple sizes approximately one order of magnitude larger than the material length scale.     

The effect of strain softening in the matrix material during void growth

A Finite element analysis has been employed to investigate the growth of an initially spherical void embedded in a cylinder of elastic-plastic material. The boundary displacement of this cylindrical cell is regulated by the value of a parameter α which controls the radial shrinkage of the cell as it elongates. A large strain analysis was used and results for both strain hardening and strain softening (after an appropriate amount of hardening has taken place) have been obtained. The effects of different mean tensile stresses, equivalent strains and initial void volume fractions have also been included. The numerical work shows relationships between the mechanical and geometrical variables that may favour ductile fracture by void coalescence or by shear decohesion.     

伴随微孔洞生长的裂纹弹塑性定常扩展

裂纹在韧性材料中扩展时,将们随着微孔洞的萌生和生长,孔洞的萌生和深化将直接影响着材料的总体断裂韧性和强度,以往的研究主要集中在将裂纹的扩展刻划为微孔洞的萌生、生长和汇合这样一个过程。从传统的断裂过程区模型出发研究微孔洞的萌生和生长对材料总体断裂韧性的影响,通过采用Ｇｕｒｓｏｎ模型,建立塑性增量本构关系,然后针对定常扩展情况直接进行分析,孔洞对材料断裂韧性的影响由本构关系刻划,而在孔洞汇合模型中,上     

Theoretical development and experimental validation of a thermally dissipative cohesive zone model for dynamic fracture of amorphous polymers

A thermally dissipative cohesive zone model is developed for predicting the temperature increase at the tip of a crack propagating dynamically in a nominally brittle material exhibiting a cohesive-type failure such as crazing. The model assumes that fracture energy supplied to the crack tip region that is in excess of that needed for the creation of new free surfaces during crack advance is converted to heat within the cohesive zone. Bulk dissipation mechanisms, such as plasticity, are not accounted for. Several cohesive traction laws are examined, and the model is then used to make predictions of crack tip heating at various crack propagation speeds in the nominally brittle amorphous polymer PMMA, observed to fail by a crazing-type mechanism. The heating predictions are compared to experimental data where the temperature field surrounding a high speed crack in PMMA was measured. Measurements are made in real time using a multi-point high speed HgCdTe infrared radiation detector array. At the same time as temperature, simultaneous measurement of fracture energy is made by a strain gauge technique, and crack tip speed is monitored through a resistance ladder method. Material strength can be estimated through uniaxial tension tests, thus minimizing the need for parameter fitting in the stress-opening traction law. Excellent agreement between experiments and theory is found for two of the cohesive traction law temperature predictions, but only for the case where a single craze is active during the dynamic fracture of PMMA, i.e. crack tip speed up to approximately 0.2c_R. For higher speed fracture where subsurface damage becomes prominent, the line dissipation model of a cohesive zone is inadequate, and a distributed damage model is needed.     

A comparison of cohesive zone modeling and classical fracture mechanics based on near tip stress field

A mode III crack with a cohesive zone in a power-law hardening material is studied under small scale yielding conditions. The cohesive law follows a softening path with the peak traction at the start of separation process. The stress and strain fields in the plastic zone, and the cohesive traction and separation displacement in the cohesive zone are obtained. The results show that for a modest hardening material (with a hardening exponent N = 0.3), the stress distribution in a large portion of the plastic zone is significantly altered with the introduction of the cohesive zone if the peak cohesive traction is less than two times yield stress, which implies the disparity in terms of the fracture prediction between the classical approach of elastic–plastic fracture mechanics and the cohesive zone approach. The stress distributions with and without the cohesive zone converge when the peak cohesive traction becomes infinitely large. A qualitative study on the equivalency between the cohesive zone approach and the classical linear elastic fracture mechanics indicates that smaller cracks require a higher peak cohesive traction than that for longer cracks if similar fracture initiations are to be predicted by the two approaches.     

Vapor pressure and void size effects on failure of a constrained ductile film

To achieve certain properties, semiconductor adhesives and molding compounds are made by blending filler particles with polymer matrix. Moisture collects at filler particle/polymer matrix interfaces and within voids of the composite. At reflow temperatures, the moisture vaporizes. The rapidly expanding vapor creates high internal pressure on pre-existing voids and particle/matrix interfaces. The simultaneous action of thermal stresses and internal vapor pressure drives both pre-existing and newly nucleated voids to grow and coalesce causing material failure. Particularly susceptible are polymeric films and adhesives joining elastic substrates, e.g. Ag filled epoxy. Several competing failure mechanisms are studied including: near-tip void growth and coalescence with the crack; extensive void growth and formation of an extended damaged zone emanating from the crack; and rapid void growth at highly stressed sites at large distances ahead of the crack, leading to multiple damaged zones. This competition is driven by the interplay between stress elevation induced by constrained plastic flow and stress relaxation due to vapor pressure assisted void growth.A model problem of a ductile film bonded between two elastic substrates, with a centerline crack, is studied. The computational study employs a Gurson porous material model incorporating vapor pressure effects. The formation of multiple damaged zones is favored when the film contains small voids or dilute second-phase particle distribution. The presence of large voids or high vapor pressure favor the growth of a self-similar damage zone emanating from the crack. High vapor pressure accelerates film cracking that can cause device failures.     

Void growth to coalescence in a non-local material

The size-effect in metals containing distributed spherical voids is analyzed numerically using a finite strain generalization of a length scale dependent plasticity theory. Results are obtained for stress-triaxialities relevant in front of a crack tip in an elastic-plastic metal. The influence of different material length parameters in a multi-parameter theory is studied, and it is shown that the important length parameter is the same as under purely hydrostatic loading. It is quantified how micron scale voids grow less rapidly than larger voids, and the implications of this in the overall strength of the material is emphasized. The size effect on the onset of coalescence is studied, and results for the void volume fraction and the strain at the onset of coalescence are presented. It is concluded that for cracked specimens not only the void volume fraction, but also the typical void size is of importance to the fracture strength of ductile materials.     

Decaying Turbulence in the Generalised Burgers Equation

This work is concerned with the derivation of optimal scaling laws, in the sense of matching lower and upper bounds on the energy, for a solid undergoing ductile fracture. The specific problem considered concerns a material sample in the form of an infinite slab of finite thickness subjected to prescribed opening displacements on its two surfaces. The solid is assumed to obey deformation-theory of plasticity and, in order to further simplify the analysis, we assume isotropic rigid-plastic deformations with zero plastic spin. When hardening exponents are given values consistent with observation, the energy is found to exhibit sublinear growth. We regularize the energy through the addition of nonlocal energy terms of the strain-gradient plasticity type. This nonlocal regularization has the effect of introducing an intrinsic length scale into the energy. Under these assumptions, ductile fracture emerges as the net result of two competing effects: whereas the sublinear growth of the local energy promotes localization of deformation to failure planes, the nonlocal regularization stabilizes this process, thus resulting in an orderly progression towards failure and a well-defined specific fracture energy. The optimal scaling laws derived here show that ductile fracture results from localization of deformations to void sheets, and that it requires a well-defined energy per unit fracture area. In particular, fractal modes of fracture are ruled out under the assumptions of the analysis. The optimal scaling laws additionally show that ductile fracture is cohesive in nature, that is, it obeys a well-defined relation between tractions and opening displacements. Finally, the scaling laws supply a link between micromechanical properties and macroscopic fracture properties. In particular, they reveal the relative roles that surface energy and microplasticity play as contributors to the specific fracture energy of the material.     

韧性断裂裂纹尖端的物理特征和力学行为

本文综述了金属材料发生韧性断裂的裂尖钝化,孔洞成核和长大,孔洞-孔洞或孔洞-裂尖汇合导致裂纹扩展的物理特征,和描述这些特征的本构理论及物理模型。      

铝合金 2A12韧性断裂机制的实验研究

用轴对称圆棒静拉伸实验研究了铝合金2A12韧性断裂的特点,并探讨了临界空穴扩张比参数-VGC在实验材料上的适用性。实验测试了实验材料的拉伸特性参数及VGC参数。实验证实实验材料在不同应力状态下会发生断裂形式的转变,光滑试样表现为剪切破坏,而缺口试样为拉伸破坏,VGC参数对材料的断裂形式较为敏感。用SEM观测了试样断口显示缺口试样为空穴聚合型断裂机制。VGC参数只适用于以拉伸型断裂为主的缺口试样,对以剪切断裂为主的光滑拉伸试样不适用。     

Mechanical modelling of ductile fracture

Mechanical models of material failure by void growth to coalescence are described to give a brief overview of methods applied in the analysis of ductile fracture. Approximate constitutive relations for porous ductile materials are discussed, modelling both the nucleation and growth of voids. The application of the material models is illustrated by numerical analyses for a tensile test specimen and for dynamic, ductile crack growth. Unstable void growth is a relevant mechanism in ductile materials subject to a high level of triaxial tension. The analysis of such cavitation instabilities in elastic-perfectly plastic materials is discussed for axisymmetric stress states, and the relevance to metal/ceramic components is emphasized.General Lecture presented at the 10th Italian National Congress of Theoretical and Applied Mechanics; AIMETA, Pisa, October 1990.     

Site of ductile fracture initiation in cold forging: A finite element model

Ductile fracture occurs due to micro-void nucleation, growth and finally coalescence into micro-crack. In this study a new ductile fracture condition that based on the microscopic phenomena of void nucleation, growth and coalescence was proposed. Using this condition and combining with finite element model to predict the fracture locations in bulk metal forming, the results show that it is in close accordance with observations of some experimental specimens. However, in order to obtaining the high trustiness many experiments have to be carried out.     

Pressure-sensitive ductile layers – II. 3D models of extensive damage

The mechanisms of void growth and coalescence in ductile polymeric layers, taking into account the effects of pressure-sensitivity, α, and plastic dilatancy, β, are explored in this two-part paper. In Part I, a two-dimensional model containing discrete cylindrical voids was used to simulate void growth and coalescence ahead of a crack. This paper extends the previous work by explicitly modeling initially spherical voids in a three-dimensional configuration. Damage predictions from the present 3D model for low yield strain adhesives are found to be in good agreement with both the 2D model in Part I and the computational cell element model. Significant discrepancies in the damage predictions, however, exist among all three models for high yield strain adhesives (e.g. polymers). The present 3D study also discusses the increasing damage level and its spatial extent with pressure-sensitivity, as well as the exacerbation of these effects arising from the deviation from an associated flow rule. In fact, both high porosity and high pressure-sensitivity promote void interaction. In addition, pressure-sensitivity increases the oblacity of the voids and reduces the intervoid ligament spacing over a wide range of load levels. These effects are compounded as the fracture process zone thickness decreases relative to the adhesive thickness. Results further show that both the adhesive toughness levels and the critical porosity governing the onset of void coalescence are significantly lowered with increasing pressure-sensitivity.     

Modeling of crack growth in ductile solids: a three-dimensional analysis

A population of several spherical voids is included in a three-dimensional, small scale yielding model. Two distinct void growth mechanisms, put forth by [Int. J. Solids Struct. 39 (2002) 3581] for the case of a two-dimensional model containing cylindrical voids, are well contained in the model developed in this study for spherical voids. A material failure criterion, based on the occurrence of void coalescence in the unit cell model, is established. The critical ligament reduction ratio, which varies with stress triaxiality and initial porosity, is used to determine ligament failure between the crack tip and the nearest void. A comparison of crack initiation toughness of the model containing cylindrical voids with the model containing spherical voids reveals that the material having a sizeable fraction of spherical voids is tougher than the material having cylindrical voids. The proposed material failure determination method is then used to establish the fracture resistance curve (J–R curve) of the material. For a ductile material containing a small volume fraction of microscopic voids initially, the void by void growth mechanism prevails, which results in a J–R curve having steep slope. On the other hand, for a ductile material containing a large volume fraction of initial voids, the multiple voids interaction mechanism prevails, which results in a flat J–R curve. Next, the effect of T-stress on fracture resistance is examined. Finally, nucleation and growth of secondary microvoids and their effects on void coalescence are briefly discussed.     

Modelling of plastic flow localisation and damage development in friction stir welded 6005A aluminium alloy using physics based strain hardening law

Plastic flow localisation and ductile failure during tensile testing of friction stir welded aluminium specimens are investigated with a specific focus on modelling the local, finite strain, hardening response. In the experimental part, friction stir welds in a 6005A-T6 aluminium alloy were prepared and analysed using digital image correlation (DIC) during tensile testing as well as scanning electron microscopy (SEM) on polished samples and on fracture surfaces. The locations of the various regions of the weld were determined based on hardness measurements, while the flow behaviour of these zones was extracted from micro-tensile specimens cut parallel to the welding direction. The measured material properties and weld topology were introduced into a 3D finite element model, fully coupled with the damage model. A Voce law hardening model involving a constant stage IV is used within an enhanced Gurson type micro-mechanical damage model, accounting for void nucleation, growth and coalescence, as well as void shape evolution. The stage IV hardening, observed in Simar et al. (2010), was found to increase the stiffness during plastic flow localisation as well as to postpone the onset of fracture as determined by the void coalescence criterion. Furthermore, the presence of a second population of voids was concluded to strongly affect the fracture strain of the high strength regions of the welds. This modelling effort links the microstructure and process parameters to macroscopic parameters relevant to the optimisation of the welds.     

Potential flow model of cavitation-induced interfacial fracture in a confined ductile layer

Fracture of a thin ductile layer sandwiched between stiff substrates often results from growth and coalescence of microscopic cavities ahead of an extending crack. Cavitation induced by plastic flow in a confined, ductile layer is analyzed here to evaluate the interfacial fracture toughness of such sandwich structures. For rigid-plastic materials, a new method is proposed in which the potential flow field of a fluid is used to approximate the plastic deformation. The principle of virtual work rate is applied to determine the equivalent traction-separation law. The method is demonstrated and validated for spherically symmetric cavity growth, for which an exact solution exists. We then study in detail the growth of an initially spherical cavity in a cylindrical bar of finite length subject to uniform traction at its ends. The results show that the stress-separation curves depend strongly on initial cavity size and the strain-hardening exponent, and weakly on the nominal strain. The method has clear advantages over numerical methods, such as finite-element analysis, for parametric study of cavity growth with large plastic deformation.