Analysis of crack growth and crack-tip plasticity in ductile materials using cohesive zone models

Crack initiation and crack growth resistance in elastic plastic materials, dominated by crack-tip plasticity are analyzed with the crack modeled as a cohesive zone. Two different types (exponential and bilinear) of cohesive zone models (CZMs) have been used to represent the mechanical behavior of the cohesive zones. In this work, it is suggested that different forms of CZMs (e.g., exponential, bilinear) are the manifestations of different micromechanisms-based inelastic processes that participate in dissipating energy during the fracture process and each form is specific to each material system. It is postulated that the total energy release rate comprises the plastic dissipation rate in the bounding material and the separation energy rate within the fracture process zone, the latter is determined by CZMs. The total energy release rate then becomes a function of the material properties (e.g., yield strength, strain hardening exponent) and cohesive properties of the fracture process zone (e.g., cohesive strength and cohesive energy), and the form of cohesive zone model (CZM) that determines the rate of energy dissipation in the forward and wake regions of the crack. The effects of material parameters, cohesive zone parameters as well as the form/shape of CZMs in predicting the crack growth resistance and the size of plastic zone (SPZ) surrounding the crack tip are systematically examined. It is found that in addition to the cohesive strength and cohesive energy, the form (shape) of the traction–separation law of CZM plays a very critical role in determining the crack growth resistance (R-curve) of a given material. It is further observed that the shape of the CZM corresponds to inelastic processes active in the forward and wake regions of the crack, and has a profound influence on the R-curve and SPZ.     

A cohesive plastic and damage zone model for dynamic crack growth in rate-dependent materials

Mode I steady-state dynamic crack growth in rate-dependent viscoplastic solids containing damage, under small scale yielding conditions, is analyzed based on a modified cohesive zone model. A multi-scale approach is used to describe the entire non-linear zone consisting of a plastic region and a damage region, each of which has its own constitutive law. Traction in the damage region is characterized by a softening power-law, in terms of the ultimate strength, a softening index and a rate sensitivity factor. In the plastic region, the cohesive law is assumed to be both strain hardening and rate dependent. The critical crack opening displacement at the physical crack-tip controls crack growth. The governing integral equations are derived and solved by a collocation method combined with associated boundary conditions. Numerical results are presented for the traction and opening profiles along the cohesive zone, the fracture energy and lengths of the damage and non-linear zones at different crack speeds and for different material parameters. The importance of factors, such as material softening, plastic deformation, crack speed and viscosity, is identified by parametric studies. In addition, the competition of plastic flow and material damage, and its effect on crack growth, are discussed.     

A material force method for inelastic fracture mechanics

A material force method is proposed for evaluating the energy release rate and work rate of dissipation for fracture in inelastic materials. The inelastic material response is characterized by an internal variable model with an explicitly defined free energy density and dissipation potential. Expressions for the global material and dissipation forces are obtained from a global balance of energy-momentum that incorporates dissipation from inelastic material behavior. It is shown that in the special case of steady-state growth, the global dissipation force equals the work rate of dissipation, and the global material force and J-integral methods are equivalent. For implementation in finite element computations, an equivalent domain expression of the global material force is developed from the weak form of the energy-momentum balance. The method is applied to model problems of cohesive fracture in a remote K-field for viscoelasticity and elastoplasticity. The viscoelastic problem is used to compare various element discretizations in combination with different schemes for computing strain gradients. For the elastoplastic problem, the effects of cohesive and bulk properties on the plastic dissipation are examined using calculations of the global dissipation force.     

Rate sensitivity of mixed mode interface toughness of dissimilar metallic materials: Studied at steady state

Crack propagation in metallic materials produces plastic dissipation when material in front for the crack tip enters the active plastic zone traveling with the tip, and later ends up being part of the residual plastic strain wake. Thus, the macroscopic work required to advance the crack is typically much larger than the work needed in the near tip fracture process. For rate sensitive materials, the amount of plastic dissipation typically depends on the rate at which the material is deformed. A dependency on the crack velocity should therefore be expected. The objective of this paper is to study the macroscopic toughness of crack advance along an interface joining two dissimilar rate dependent materials, characterized by an elastic-viscoplastic material model that approaches the response of a J₂-flow material in the rate independent limit. The emphasis here is on the rate sensitivity of the macroscopic fracture toughness under mixed Mode I/II loading. Moreover, special cases of joined similar rate dependent materials, as well as dissimilar materials where one substrate remains either elastic or approaches the rate independent limit is also included. The numerical analysis is carried out using the SSV model [Suo, Z., Shih, C., Varias, A., 1993. A theory for cleavage cracking in the presence of plastic flow. Acta Metall. Mater. 41, 1551–1557] embedded in a steady state finite element formulation, here assuming plane strain conditions and small-scale yielding. Results are presented for a wide range of material parameters, including noteworthy observations of a characteristic crack velocity at which the macroscopic toughness becomes independent of the material rate sensitivity. The potential of this phenomenon is elaborated on from a modeling point of view.     

Constraint effects in adhesive joint fracture

Constraint effects in adhesive joint fracture are investigated by modelling the adherents as well as a finite thickness adhesive layer in which a single row of cohesive zone elements representing the fracture process is embedded. Both the adhesive and the adherents are elastic-plastic with strain hardening. The bond toughness Γ (work per unit area) is equal to Γ₀+Γ_p, where Γ₀ is the intrinsic work of fracture associated with the embedded cohesive zone response and Γ_p is the extra contribution to the bond toughness arising from plastic dissipation and stored elastic energy within the adhesive layer. The parameters of the model are identified from experiments on two different adhesives exhibiting very different fracture properties. Most of the tests were performed using the wedge-peel test method for a variety of adhesives, adherents and wedge thicknesses. The model captures the constraint effects resulting from the change in Γ_p: (i) the plastic dissipation increases with increasing bond line thickness in the fully plastic regime and then decreases to reach a constant value for very thick adhesive layers; (ii) the plastic dissipation in the fully plastic regime increases drastically as the thickness of the adherent decreases. Finally, this model is used to assess a simpler approach which consists of simulating the full adhesive layer as a single row of cohesive elements.     

基于能量等效原理的应变局部化分析:Ⅰ.一维解析解

基于热力学第一定律和非局部塑性理论,提出了一种求解应变局部化问题的非局部方法.对材料的每一点定义了局部和非局部两种状态空间,局部状态空间的内变量通过非局部权函数映射到非局部空间,成为非局部内变量.在应变软化过程中,局部状态空间中的塑性变形服从正交流动法则,材料的软化律在非局部状态空间中被引入.通过两个状态空间的塑性应变能耗散率的等效,得到了应变软化过程中明确定义的局部化区域以及其中的塑性应变分布.应用本方法导出了一维应变局部化问题的解析解.解析解表明,应变局部化区域的尺寸只与材料内尺度有关;对于高斯型非局部权函数,局部化区域的尺寸大约是材料内尺度的6倍.一维算例表明,局部化区域的塑性应变分布以及载荷-位移曲线仅与材料参数和结构几何尺寸有关,变形局部化区域的尺寸随着材料内尺度的减小而减小,同时塑性应变也随着材料内尺度的减小变得更加集中.当内尺度趋近于零时,应用本文方法得到的解与采用传统的局部塑性理论得到的解相同.     

Intergranular strain evolution near fatigue crack tips in polycrystalline metals

The deformation field near a steady fatigue crack includes a plastic zone in front of the crack tip and a plastic wake behind it, and the magnitude, distribution, and history of the residual strain along the crack path depend on the stress multiaxiality, material properties, and history of stress intensity factor and crack growth rate. An in situ, full-field, non-destructive measurement of lattice strain (which relies on the intergranular interactions of the inhomogeneous deformation fields in neighboring grains) by neutron diffraction techniques has been performed for the fatigue test of a Ni-based superalloy compact tension specimen. These microscopic grain level measurements provided unprecedented information on the fatigue growth mechanisms. A two-scale model is developed to predict the lattice strain evolution near fatigue crack tips in polycrystalline materials. An irreversible, hysteretic cohesive interface model is adopted to simulate a steady fatigue crack, which allows us to generate the stress/strain distribution and history near the fatigue crack tip. The continuum deformation history is used as inputs for the micromechanical analysis of lattice strain evolution using the slip-based crystal plasticity model, thus making a mechanistic connection between macro- and micro-strains. Predictions from perfect grain-boundary simulations exhibit the same lattice strain distributions as in neutron diffraction measurements, except for discrepancies near the crack tip within about one-tenth of the plastic zone size. By considering the intergranular damage, which leads to vanishing intergranular strains as damage proceeds, we find a significantly improved agreement between predicted and measured lattice strains inside the fatigue process zone. Consequently, the intergranular damage near fatigue crack tip is concluded to be responsible for fatigue crack growth.     

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.     

Fatigue damage modeling in solder interconnects using a cohesive zone approach

The objective of this work is to model the fatigue damage process in a solder bump subjected to cyclic loading conditions. Fatigue damage is simulated using the cohesive zone methodology. Damage is assumed to occur at interfaces modeled through cohesive zones in the material, while the bulk material is assumed to be linear elastic. The state of damage at a cohesive zone is incorporated into the cohesive zone constitutive law by a elasticity-based damage variable. The gradual degradation of the solder material and the corresponding damage accumulation throughout the cycling process is accounted for by a damage evolution law which captures the main damage characteristics. The model prediction of the solder bump life-time is shown to be in good agreement with one of the commonly used empirical life-time prediction laws.     

Plastic Strain Energy Model for Rock Salt Under Fatigue Loading

The fatigue test for rock salt is conducted to investigate the effects of stress amplitude, loading frequency and loading rate on the plastic strain energy, from which the evaluation rule of the plastic strain energy is analyzed, which is divided into three stages: cyclic hardening,saturation and cyclic softening. The total accumulated plastic strain energy only depends on the mechanical behavior of rock salt, but is immune to the loading conditions. A novel model for fatigue life prediction is proposed based on the invariance of the total plastic dissipation energy and the stability of the plastic energy per cycle.     

微纳米晶金属的应变率敏感性及应变硬化行为分析

基于亚微米、纳米晶粒组织塑性变形过程中多种变形机制(位错机制、扩散机制及晶界滑动机制)共存,建立了理论模型,用于定量研究亚微米、纳米晶粒组织的塑性变形行为.以铜为模型材料,计算分析了晶粒尺度、应变率以及温度对亚微米、纳米晶粒组织塑性变形行为的影响.结果表明:相比粗晶铜,亚微米晶铜表现出明显的应变率敏感性,并且应变率敏感系数随晶粒尺度及变形速率的减小而增大;同时,增大变形速率或降低变形温度都能提高材料的应变硬化能力,延缓颈缩发生,进而提高材料的延性.计算分析结果与实验报道吻合.     

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.     

Electro-mechanical pre-fracture zones for an electrically permeable interface crack in a piezoelectric bimaterial

A plane strain problem for two piezoelectric half-spaces adhered by a very thin isotropic interlayer with a crack under the action of remote mixed mode mechanical loading and electrical flux is considered. The crack is situated either at an interface or in the interlayer. It is assumed that the substrates are much stiffer than the intermediate layer. Therefore, pre-fracture zones (plastic or damage) arise at the crack continuations. Normal and shear stresses are assumed to be constant in this zones and to satisfy some material equation, which can be taken from theory or derived experimentally. Modeling the pre-fracture zones by the crack continuations with unknown cohesive stresses on their faces reduces the problem to elastic interface crack analysis leading to a Hilbert problem. This problem is solved exactly. The pre-fracture zone lengths and stresses in these zones are found from algebraical and transcendental equations. The latter are derived from the conditions of stress finiteness at the ends of pre-fracture zones and the material equations. The electrical displacement at any point of the pre-fracture zones is found in closed form as well. Particular cases of symmetrical loading and of equivalent properties of the upper and lower bimaterial components are considered. Numerical results corresponding to certain material combinations and interlayer material equations are presented and analysed. In the suggested model, any singularities connected with the crack are eliminated, i.e., all mechanical and electrical characteristics are limited in the near-crack tip region.     

Constitutive modeling for nanocrystalline metals based on cooperative grain boundary mechanisms

In nanocrystalline metals, the plastic deformation is accommodated primarily at the grain boundaries. Yang and Wang (J. Mech. Phys. Solids, 2003) suggested a deformation model based on clusters consisting of nine grains and incorporating both the Ashby-Verrall mechanism and a 30° rotation of closely linked pairs of grains. In the present article, the insertion and rotation processes are considered together as a cooperative deformation mechanisms, and the degree to which each process contributes is determined by the application of the principle of maximum plastic work. Plane strain and three-dimensional constitutive relations based on this concept are derived for which a general stress state drives the orientation evolution of various grain clusters under the Reuss assumption. Detailed calculation shows that the strain rate depends linearly on the stress, with the values of the coefficients in this linear relationship dictated by the microscopic energy dissipation. The deformation contributed to the overall response by the grain boundary mechanism is discussed in the spirit of the Hashin-Shtrikman bounds.     

A unified residual-based thermodynamic framework for strain gradient theories of plasticity

A unified thermodynamic framework for gradient plasticity theories in small deformations is provided, which is able to accommodate (almost) all existing strain gradient plasticity theories. The concept of energy residual (the long range power density transferred to the generic particle from the surrounding material and locally spent to sustain some extra plastic power) plays a crucial role. An energy balance principle for the extra plastic power leads to a representation formula of the energy residual in terms of a long range stress, typically of the third order, a macroscopic counterpart of the micro-forces acting on the GNDs (Geometrically Necessary Dislocations). The insulation condition (implying that no long range energy interactions are allowed between the body and the exterior environment) is used to derive the higher order boundary conditions, as well as to ascertain a principle of the plastic power redistribution in which the energy residual plays the role of redistributor and guarantees that the actual plastic dissipation satisfies the second thermodynamics principle. The (nonlocal) Clausius-Duhem inequality, into which the long range stress enters aside the Cauchy stress, is used to derive the thermodynamic restrictions on the constitutive equations, which include the state equations and the dissipation inequality. Consistent with the latter inequality, the evolution laws are formulated for rate-independent models. These are shown to exhibit multiple size effects, namely (energetic) size effects on the hardening rate, as well as combined (dissipative) size effects on both the yield strength (intrinsic resistance to the onset of plastic strain) and the flow strength (resistance exhibited during plastic flow). A friction analogy is proposed as an aid for a better understanding of these two kinds of strengthening effects. The relevant boundary-value rate problem is addressed, for which a solution uniqueness theorem and a minimum variational principle are provided. Comparisons with other existing gradient theories are presented. The dissipation redistribution mechanism is illustrated by means of a simple shear model.     

Effect of crack-tip shape on the near-tip field in glassy polymer

The physical nature of a crack tip is not absolutely sharp but blunt with finite curvature. In this paper, the effects of crack-tip shape on the stress and deformation fields ahead of blunted cracks in glassy polymers are numerically investigated under Mode I loading and small scale yielding conditions. An elastic–viscoplastic constitutive model accounting for the strain softening upon yield and then the subsequently strain hardening is adopted and two typical glassy polymers, one with strain hardening and the other with strain softening–rehardening are considered in analysis. It is shown that the profile of crack tip has obvious effect on the near-tip plastic field. The size of near-tip plastic zone reduces with the increase of curvature radius of crack tip, while the plastic strain rate and the stresses near crack tip enhance obviously for two typical polymers. Also, the plastic energy dissipation behavior near cracks with different curvatures is discussed for both materials.     

Dynamic plasticity and fracture in high density polycrystals: constitutive modeling and numerical simulation

Presented is a constitutive framework for modeling the dynamic response of polycrystalline microstructures, posed in a thermodynamically consistent manner and accounting for finite deformation, strain rate dependence of flow stress, thermal softening, thermal expansion, heat conduction, and thermoelastic coupling. Assumptions of linear and square-root dependencies, respectively, of the stored energy and flow stresses upon the total dislocation density enable calculation of the time-dependent fraction of plastic work converted to heat energy. Fracture at grain boundary interfaces is represented explicitly by cohesive zone models. Dynamic finite element simulations demonstrate the influences of interfacial separation, random crystallographic orientation, and grain morphology on the high-rate tensile response of a realistic two-phase material system consisting of comparatively brittle pure tungsten (W) grains embedded in a more ductile matrix of tungsten-nickel iron (W-Ni-Fe) alloy. Aspects associated with constitutive modeling of damage and failure in the homogenized material system are discussed in light of the computational results.     

STOCHASTIC ELASTO-PLASTIC FRACTURE ANALYSIS OF ALUMINUM FOAMS

Model I quasi-static nonlinear fracture of aluminum foams is analyzed by considering the effect of microscopic heterogeneity. Firstly, a continuum constitutive model is adopted to account for the plastic compressibility of the metallic foams. The yield strain modeled by a two- parameter Weibull-type function is adopted in the constitutive model. Then, a modified cohesive zone model is established to characterize the fracture behavior of aluminum foams with a cohesive zone ahead of the initial crack. The tensile traction versus local crack opening displacement relation is employed to describe the softening characteristics of the material. And a Weibull statistical model for peak bridging stress within the fracture process zone is used for considering microscopic heterogeneity of aluminum foams. Lastly, the influence of stochastic parameters on the curve of stress-strain is given. Numerical examples are given to illustrate the numerical model presented in this paper and the effects of Weibull parameters and material properties on J-integral are discussed.