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.     

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.     

An energy-based crack growth criterion for modelling elastic–plastic ductile fracture

A crack growth criterion is derived based on the Griffith energy concept and the cohesive zone model for modelling fracture in elastic–plastic ductile materials. The criterion is implemented in the finite element context by a virtual crack extension technique. An automatic modelling of the ductile fracture process is realised by combining a local remeshing procedure and the criterion. The validity of the derived criterion is examined by modelling a compact tension specimen.     

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.     

On the extension of the Gurson-type porous plasticity models for prediction of ductile fracture under shear-dominated conditions

One of the major drawbacks of the Gurson-type of porous plasticity models is the inability of these models to predict material failure under low stress triaxiality, shear dominated conditions. This study addresses this issue by combining the damage mechanics concept with the porous plasticity model that accounts for void nucleation, growth and coalescence. In particular, the widely adopted Gurson–Tvergaard–Needleman (GTN) model is extended by coupling two damage parameters, representing the volumetric damage (void volume fraction) and the shear damage, respectively, into the yield function and flow potential. The effectiveness of the new model is illustrated through a series of numerical tests comparing its performance with existing models. The current model not only is capable of predicting damage and fracture under low (even negative) triaxiality conditions but also suppresses spurious damage that has been shown to develop in earlier modifications of the GTN model for moderate to high triaxiality regimes. Finally the modified GTN model is applied to predict the ductile fracture behavior of a beta-treated Zircaloy-4 by coupling the proposed damage modeling framework with a recently developed J₂–J₃ plasticity model for the matrix material. Model parameters are calibrated using experimental data, and the calibrated model predicts failure initiation and propagation in various specimens experiencing a wide range of triaxiality and Lode parameter combinations.     

Effects of the strain-hardening exponent on two-parameter characterizations of surface-cracks under large-scale yielding

The influence of strain hardening exponent on two-parameter J-Q near tip opening stress field characterization with modified boundary layer formulation and the corresponding validity limits are explored in detail. Finite element simulations of surface cracked plates under uniaxial tension are implemented for loads exceeding net-section yield. The results from this study provide numerical methodology for limit analysis and demonstrate the strong material dependencies of fracture parameterization under large scale yielding. Sufficient strain hardening is shown to be necessary to maintain J-Q predicted fields when plastic flow progresses through the remaining ligament. Lower strain hardening amplifies constraint loss due to stress redistribution in the plastic zone and increases the ratio of tip deformation to J. The onset of plastic collapse is marked by shape change and/or rapid relaxation of tip fields compared to those predicted by MBL solutions and thus defining the limits of J-Q dominance. A radially independent Q-parameter cannot be evaluated for the low strain hardening material at larger deformations within a range where both cleavage and ductile fracture mechanisms are present. The geometric deformation limit of near tip stress field characterization is shown to be directly proportional to the level of stress the material is capable of carrying within the plastic zone. Accounting for the strain hardening of a material provides a more adjusted and less conservative limit methodology compared to those generalized by the yield strength alone. Results from this study are of relevance to establishing testing standards for surface cracked tensile geometries.     

On some topics for the numerical simulation of ductile fracture

In this work, we analyze some aspects of the macroscopic Gurson–Tvergaard–Needleman (GTN) constitutive model when it is addressed to solve ductile fracture problems by means of numerical simulations:     

Void nucleation and growth in mineral-filled PVC – An experimental and numerical study

The nucleation and growth of voids in mineral-filled PVC have been investigated through experimental and numerical studies. Uniaxial tensile specimens were deformed in tension to different elongation levels and then unloaded. The macroscopic strain fields were recorded by use of digital image correlation. After the test, the microstructure of the deformed specimens was investigated in a scanning electron microscope. It was found that the observed volume strain on the macroscale is related to void growth on the microscale. In addition, finite element simulations were performed on unit cell models representing the microstructure of the material in a simplified manner. The numerical simulations demonstrate macroscopic dilation as a result of void growth. Moreover, the numerical simulations indicate that the experimentally observed stress-softening response of the PVC composite material may result from matrix-particle debonding.     

One, no one, and one hundred thousand crack propagation laws: A generalized Barenblatt and Botvina dimensional analysis approach to fatigue crack growth

Barenblatt and Botvina with elegant dimensional analysis arguments have elucidated that Paris’ power-law is a weak form of scaling, so that the Paris’ parameters C and m should not be taken as material constants. On the contrary, they are expected to depend on all the dimensionless parameters of the problem, and are really “constants” only within some specific ranges of all these. In the present paper, the dimensional analysis approach by Barenblatt and Botvina is generalized to explore the functional dependencies of m and C on more dimensionless parameters than the original Barenblatt and Botvina, and experimental results are interpreted for a wider range of materials including both metals and concrete. In particular, we find that the size-scale dependencies of m and C and the resulting correlation between C and m are quite different for metals and for quasi-brittle materials, as it is already suggested from the fact the fatigue crack propagation processes lead to m=2-5 in metals and m=10-50 in quasi-brittle materials. Therefore, according to the concepts of complete and incomplete self-similarities, the experimentally observed breakdowns of the classical Paris’ law are discussed and interpreted within a unified theoretical framework. Finally, we show that most attempts to address the deviations from the Paris’ law or the empirical correlations between the constants can be explained with this approach. We also suggest that “incomplete similarity” corresponds to the difficulties encountered so far by the “damage tolerant” approach which, after nearly 50 years since the introduction of Paris’ law, is still not a reliable calculation of damage, as Paris himself admits in a recent review.     

Modeling the evolution of stress due to differential shrinkage in powder-processed functionally graded metal–ceramic composites during pressureless sintering

Pressureless sintering of powder-processed functionally graded materials is being pursued to economically produce metal–ceramic composites for a variety of high-temperature (e.g., thermal protection) and energy-absorbing (e.g., armor) applications. During sintering, differential shrinkage induces stresses that can compromise the integrity of the components. Because the strength evolves as the component is sintered, it is important to model how the evolution of the differential shrinkage governs the stress distribution in the component in order to determine when the strength will be exceeded and cracking initiated. In this investigation, a model is proposed that describes the processing/microstructure/property/performance relationship in pressurelessly sintered functionally graded plates and rods. This model can be used to determine appropriate shrinkage rates and gradient architectures for a given component geometry that will prevent the component from cracking during pressureless sintering by balancing the evolution of strength, which is assumed to be a power law function of the porosity, with the evolution of stress. To develop this model, the powder mixture is considered as a three-phase material consisting of voids, metal particles, and ceramic particles. A micromechanical thermal elastic–viscoplastic constitutive model is then proposed to describe the thermomechanical behavior of the composite microstructure. The subsequent evolution of the thermomechanical properties of the matrix material during sintering is assumed to obey a power law relationship with the level of porosity, which is directly related to the shrinkage strain, and was refined to account for the evolving interparticle cohesion of the matrix phase due to sintering. These thermomechanical properties are incorporated into a 2-D thermomechanical finite element analysis to predict the stress distributions and distortions that arise from the evolution of differential shrinkage during the pressureless sintering process. Differential shrinkage results were verified quantitatively through comparison with the shape profile for a pressurelessly sintered functionally graded nickel–alumina composite plate with a cylindrical geometry, and the stress distribution results verified from qualitative observations of the absence or presence of cracking as well as the location in specimens with different gradient architectures. The cracking was mitigated using a reverse gradient at one end of the specimen, and the resulting distortions associated with the shape profile were determined to be no more than 15% reduced from the predictions. The effects of geometry were also studied out-of-plane by transforming the plate into a rod through an increase in thickness, while in-plane effects were studied by comparing the results from the cylindrical specimen with a specimen that has a square cross-sectional geometry. By transforming from a plate to a rod geometry, the stress no longer exceeds critical levels and cracks do not form. The results from the in-plane geometric study indicated that critical stresses were reached in the square geometry at temperatures 100 °C less than in the cylindrical geometry. Additionally, the location of primary cracking was shifted towards the metal-rich end of the specimen, while the stress distribution associated with this shift and the lower temperature for the critical stress resulted in secondary cracking.