Micromechanics of coalescence in ductile fracture

Significant progress has been recently made in modelling the onset of void coalescence by internal necking in ductile materials. The aim of this paper is to develop a micro-mechanical framework for the whole coalescence regime, suitable for finite-element implementation. The model is defined by a set of constitutive equations including a closed form of the yield surface along with appropriate evolution laws for void shape and ligament size. Normality is still obeyed during coalescence. The derivation of the evolution laws is carefully guided by coalescence phenomenology inferred from micromechanical unit-cell calculations. The major implication of the model is that the stress carrying capacity of the elementary volume vanishes as a natural outcome of ligament size reduction. Moreover, the drop in the macroscopic stress accompanying coalescence can be quantified for many initial microstructures provided that the microstructure state is known at incipient coalescence. The second part of the paper addresses a more practical issue, that is the prediction of the acceleration rate δ in the Tvergaard-Needleman phenomenological approach to coalescence. For that purpose, a Gurson-like model including void shape effects is used. Results are presented and discussed in the limiting case of a non-hardening material for different initial microstructures and various stress states. Predicted values of δ are extremely sensitive to stress triaxiality and initial spacing ratio. The effect of initial porosity is significant at low triaxiality whereas the effect of initial void shape is emphasized at high triaxiality.     

A three-dimensional crystal plasticity model for duplex Ti–6Al–4V

A rate dependent crystal plasticity model for the α/β Ti–Al alloy Ti–6Al–4V with duplex microstructure is developed and presented herein. Duplex Ti–6Al–4V is a dual-phase alloy consisting of an hcp structured matrix primary α-phase and secondary lamellar α + β domains that are composed of alternating layers of secondary α laths and bcc structured residual β laths. The model accounts for distinct three-dimensional slip geometry for each phase, anisotropic and length scale dependent slip system strengths, the non-planar dislocation core structure of prismatic screw dislocations in the primary α-phase, and crystallographic texture. The model is implemented in the general purpose finite element code (ABAQUS, 2005. Ver 6.5, Hibbitt, Karlsson, and Sorensen, Inc., Pawtucket, RI) via a UMAT subroutine.     

Multiscale modeling of plasticity based on embedding the viscoplastic self-consistent formulation in implicit finite elements

This paper is concerned with the multiscale simulation of plastic deformation of metallic specimens using physically-based models that take into account their polycrystalline microstructure and the directionality of deformation mechanisms acting at single-crystal level. A polycrystal model based on self-consistent homogenization of single-crystal viscoplastic behavior is used to provide a texture-sensitive constitutive response of each material point, within a boundary problem solved with finite elements (FE) at the macroscale. The resulting constitutive behavior is that of an elasto-viscoplastic material, implemented in the implicit FE code ABAQUS. The widely-used viscoplastic selfconsistent (VPSC) formulation for polycrystal deformation has been implemented inside a user-defined material (UMAT) subroutine, providing the relationship between stress and plastic strain-rate response. Each integration point of the FE model is considered as a polycrystal with a given initial texture that evolves with deformation. The viscoplastic compliance tensor computed internally in the polycrystal model is in turn used for the minimization of a suitable-designed residual, as well as in the construction of the elasto-viscoplastic tangent stiffness matrix required by the implicit FE scheme.Uniaxial tension and simple shear of an FCC polycrystal have been used to benchmark the accuracy of the proposed implicit scheme and the correct treatment of rotations for prediction of texture evolution. In addition, two applications are presented to illustrate the potential of the multiscale strategy: a simulation of rolling of an FCC plate, in which the model predicts the development of different textures through the thickness of the plate; and the deformation under 4-point bending of textured HCP bars, in which the model captures the dimensional changes associated with different orientations of the dominant texture component with respect to the bending plane.     

Necking in glassy polymers: Effects of intrinsic anisotropy and structural evolution kinetics in their viscoplastic flow

Two recently proposed developments of the Glass–Rubber constitutive model for glassy polymers treat the viscoplastic deformation as intrinsically anisotropic, and incorporate the kinetics of structural evolution. These features enable the model to capture better the distinctive features of glassy polymers’ constitutive response: post-yield strain-softening and strain-hardening and effects of pre-existing molecular orientation. They have been combined to form a new variant of the model, and the consequences for necking have been explored. Uniaxial extension of prismatic bars was simulated using the finite element method, employing a numerical implementation of the new model, with material parameters of polystyrene. Strain localization predicted with the new model was found to be systematically retarded as compared to predictions with the original (intrinsically isotropic) version of the model, for the same conditions. In particular, the effect of frozen-in molecular orientation was examined. This was found to retard strain localization for stretching parallel to the orientation direction, for both models. But the localization predicted with the new model was always significantly less pronounced than with the original model. Indeed, for sufficiently high pre-orientation (e.g. a uniaxial stretch of 2.2), localization could be effectively prevented with the new model, under conditions when otherwise failure by necking is predicted. Such results can all be explained in terms of a linear stability analysis. They suggest that all previous simulations of necking in glassy polymers made using intrinsically isotropic representations of polymer viscoplasticity may have over-predicted the rate of strain localization.     

Crystal plasticity modeling of slip activity in Ti–6Al–4V under high cycle fatigue loading

Deformation micromechanisms of a Ti–6Al–4V alloy under fatigue loading at room temperature are studied using a three-dimensional crystal plasticity constitutive model. The model employs a minimum set of fitting parameters based on experimental data for Ti–6Al–4V. Single slip is strongly favored through a softening law that affects mainly the driving force for slip on the first activated slip system. Cyclic deformation behavior at the macroscopic scale and at the local scale of grains is analyzed through the simulation of 20 cycles of fatigue on a polycrystalline structure of 900 randomly oriented grains. The progressive activation of slip (basal, prismatic, and pyramidal) is analyzed and compared to experimental observations.     

Size-dependent effective thermoelastic properties of nanocomposites with spherically anisotropic phases

Composites made of semi-crystalline polymers and nanoparticles have a spherulitic microstructure which can be reasonably represented by a spherically anisotropic volume element. Due to the high surface-to-volume ratio of a nanoparticle, the particle-matrix interface stress, usually neglected in determining the effective elastic moduli of particle-reinforced composites, may have a non-negligible effect. To account for the latter in estimating the effective thermoelastic properties of a composite consisting of nanoparticles embedded in a semi-crystalline polymeric matrix, this work adopts a coherent interface model for the nanoparticle-matrix interface and proposes an extended version of the classical generalized-self consistent method. In particular, Eshelby's formulae widely used to calculate the elastic energy change of a homogeneous medium due to the introduction of an inhomogeneity are extended to the thermoelastic case. The nanoparticle size effect on the effective thermoelastic moduli of the composite are theoretically shown and numerically illustrated.     

Effects of spatial grain orientation distribution and initial surface topography on sheet metal necking

The finite element method is used to numerically simulate localized necking in AA6111-T4 under stretching. The measured EBSD data (grain orientations and their spatial distributions) are directly incorporated into the finite element model and the constitutive response at an integration point is described by the single crystal plasticity theory. We assume that localized necking is associated with surface instability, the onset of unstable growth in surface roughening. It is demonstrated that such a surface instability/necking is the natural outcome of the present approach, and the artificial initial imperfection necessitated by the macroscopic M–K approach [Marciniak and Kuczynski (1967). Int. J. Mech. Sci. 9, 609–620] is not relevant in the present analysis. The effects of spatial orientation distribution, material strain rate sensitivity, texture evolution, and initial surface topography on necking are discussed. It is found that localized necking depends strongly on both the initial texture and its spatial orientation distribution. It is also demonstrated that the initial surface topography has only a small influence on necking.     

High temperature creep behaviour of an FeAl intermetallic strengthened by nanoscale oxide particles

The creep behaviour of an FeAl intermetallic strengthened by nanosized oxide particles has been examined at temperatures of 700–825 °C. For all temperatures the strain rate shows a power law dependence on the applied stress. At the lowest temperature and with the highest stresses there is evidence of a threshold stress produced by the difficulty of overcoming the particle barriers, while for higher temperatures as well as at low stresses there is no threshold stress and creep appears to be controlled by general climb. The fine oxide particles produce good strengthening at low temperatures but are more readily overcome at high temperatures due to their very small size and limited attractive relaxation force. Despite such fall in creep strength, this material remains one of the strongest iron aluminides to the temperature range evaluated.     

Effects of texture on shear band formation in plane strain tension/compression and bending

In this study, effects of typical texture components observed in rolled aluminum alloy sheets on shear band formation in plane strain tension/compression and bending are systematically studied. The material response is described by a generalized Taylor-type polycrystal model, in which each grain is characterized in terms of an elastic–viscoplastic continuum slip constitutive relation. First, a simple model analysis in which the shear band is assumed to occur in a weaker thin slice of material is performed. From this simple model analysis, two important quantities regarding shear band formation are obtained: i.e. the critical strain at the onset of shear banding and the corresponding orientation of shear band. Second, the shear band development in plane strain tension/compression is analyzed by the finite element method. Predictability of the finite element analysis is compared to that of the simple model analysis. Third, shear band developments in plane strain pure bending of a sheet specimen with the typical textures are studied. Regions near the surfaces in a bent sheet specimen are approximately subjected to plane strain tension or compression. From this viewpoint, the bendability of a sheet specimen may be evaluated, using the knowledge regarding shear band formation in plane strain tension/compression. To confirm this and to encompass overall deformation of a bent sheet specimen, including shear bands, finite element analyses of plane strain pure bending are carried out, and the predicted shear band formation in bent specimens is compared to that in the tension/compression problem. Finally, the present results are compared to previous related studies, and the efficiency of the present method for materials design in future is discussed.     

Advances in experiments on metal sheets and tubes in support of constitutive modeling and forming simulations

This work is a review of experimental methods for observing and modeling the anisotropic plastic behavior of metal sheets and tubes under a variety of loading paths, such as biaxial compression tests; biaxial tension tests on metal sheets and tubes using closed-loop electrohydraulic testing machines; the abrupt strain path change method for detecting a yield vertex and subsequent yield loci without unloading; in-plane stress reversal tests on metal sheets; and multistage tension tests. Observed material responses are compared with the predictions of phenomenological plasticity models. Special attention is paid to the plastic deformation behavior of materials commonly used in industry, and to verifying the validity of conventional anisotropic yield criteria for those materials and associated flow rules at large plastic strains. The effects of using appropriate anisotropic yield criteria on the accuracy of simulations of forming defects, such as large springback and fracture, are also presented to highlight the importance of accurate material testing and modeling.     

A non-associated constitutive model with mixed iso-kinematic hardening for finite element simulation of sheet metal forming

In this paper an anisotropic material model based on non-associated flow rule and mixed isotropic–kinematic hardening was developed and implemented into a user-defined material (UMAT) subroutine for the commercial finite element code ABAQUS. Both yield function and plastic potential were defined in the form of Hill’s [Hill, R., 1948. A theory of the yielding and plastic flow of anisotropic metals. Proc. R. Soc. Lond. A 193, 281–297] quadratic anisotropic function, where the coefficients for the yield function were determined from the yield stresses in different material orientations, and those of the plastic potential were determined from the r-values in different directions. Isotropic hardening follows a nonlinear behavior, generally in the power law form for most grades of steel and the exponential law form for aluminum alloys. Also, a kinematic hardening law was implemented to account for cyclic loading effects. The evolution of the backstress tensor was modeled based on the nonlinear kinematic hardening theory (Armstrong–Frederick formulation). Computational plasticity equations were then formulated by using a return-mapping algorithm to integrate the stress over each time increment. Either explicit or implicit time integration schemes can be used for this model. Finally, the implemented material model was utilized to simulate two sheet metal forming processes: the cup drawing of AA2090-T3, and the springback of the channel drawing of two sheet materials (DP600 and AA6022-T43). Experimental cyclic shear tests were carried out in order to determine the cyclic stress–strain behavior and the Bauschinger ratio. The in-plane anisotropy (r-value and yield stress directionalities) of these sheet materials was also compared with the results of numerical simulations using the non-associated model. These results showed that this non-associated, mixed hardening model significantly improves the prediction of earing in the cup drawing process and the prediction of springback in the sidewall of drawn channel sections, even when a simple quadratic constitutive model is used.     

A finite element formulation based on non-associated plasticity for sheet metal forming

In the present paper, a finite element formulation based on non-associated plasticity is developed. In the constitutive formulation, isotropic hardening is assumed and an evolution equation for the hardening parameter consistent with the principle of plastic work equivalence is introduced. The yield function and plastic potential function are considered as two different functions with functional form as the yield function of Hill [Hill, R., 1948. Theory of yielding and plastic flow of anisotropic metals. Proc. Roy. Soc. A 193, 281–297] or Karafillis–Boyce associated model [Karafillis, A.P. Boyce, M., 1993. A general anisotropic yield criterion using bounds and a transformation weighting tensor. J. Mech. Phys. Solids 41, 1859–1886]. Algorithmic formulations of constitutive models that utilize associated or non-associated flow rule coupled with Hill or Karafillis–Boyce stress functions are derived by application of implicit return mapping procedure. Capabilities in predicting planar anisotropy of the Hill and Karafillis–Boyce stress functions are investigated considering material data of Al2008-T4 and Al2090-T3 sheet samples. The accuracy of the derived stress integration procedures is investigated by calculating iso-error maps.     

Exactly solvable spherically anisotropic thermoelastic microstructures

Microstructures possessing local spherical anisotropy are considered in this paper. An example is a spherulitic polymer which can be modelled by an assemblage of spheres of all sizes in which a radial direction in every sphere is an axis of local transverse isotropy. Our purpose is to construct effectively isotropic microstructures, with spherically anisotropic and thermoelastic constituents, whose effective bulk modulus, thermal stress and specific heat can be exactly determined. The basic microstructure for which this is achieved in the present paper is the nested composite sphere assemblage of Milgrom and Shtrikman (J. Appl. Phys. 66 (1989) 3429) which was originally formulated for isotropic constituents, in the settings of conductivity and coupled fields with scalar potentials. Here, we allow the phases of this microstructure to be spherically thermoelastic with a symmetry class which can be trigonal, tetragonal, transversely isotropic, cubic or isotropic with respect to a local spherical coordinate system. A rich class of new exact results for two-phase composites and polycrystals are obtained.     

Constitutive framework optimized for myocardium and other high-strain, laminar materials with one fiber family

Central to this analysis is the identification of six rotation invariant scalars α_1-6 that succinctly define the strain in materials that have one family of parallel fibers arranged in laminae. These scalars were chosen so as to minimize covariance amongst the response terms in the hyperelastic limit, and they are termed strain attributes because it is necessary to distinguish them from strain invariants. The Cauchy stress t is expressed as the sum of six response terms, almost all of which are mutually orthogonal for finite strain (i.e. 14 of the 15 inner products vanish). For small deformations, the response terms are entirely orthogonal (i.e. all 15 inner products vanish). A response term is the product of a response function with its associated kinematic tensor. Each response function is a scalar partial derivative of the strain energy W with respect to a strain attribute. Applications for this theory presently include myocardium (heart muscle) which is often modeled as having muscle fibers arranged in sheets. Utility for experimental identification of strain energy functions is demonstrated by showing that common tests on incompressible materials can directly determine terms in W. Since the described set of strain attributes reduces the covariance amongst response terms, this approach may enhance the speed and precision of inverse finite element methods.