Modeling transients in the mechanical response of copper due to strain path changes

When copper is deformed to large strains its texture and microstructure change drastically, leading to plastic anisotropy and extended transients when it is reloaded along a different strain path. For predicting these transients, we develop a constitutive model for polycrystalline metals that incorporates texture and grain microstructure. The directional anisotropy in the single crystals is considered to be induced by variable latent hardening associated with cross-slip, cut-through of planar dislocation walls, and dislocation-based reversal mechanisms. These effects are introduced in a crystallographic hardening model which is, in turn, implemented into a polycrystal model. This approach successfully explains the flow response of OFHC Cu pre-loaded in tension (compression) and reloaded in tension (compression), and the response of OFHC Cu severely strained in shear by equal channel angular extrusion and subsequently compressed in each of the three orthogonal directions. This new theoretical framework applies to arbitrary strain path changes, and is fully anisotropic.     

A multiscale approach for modeling scale-dependent yield stress in polycrystalline metals

Modeling of scale-dependent characteristics of mechanical properties of metal polycrystals is studied using both discrete dislocation dynamics and continuum crystal plasticity. The initial movements of dislocation arc emitted from a Frank-Read type dislocation source and bounded by surrounding grain boundaries are examined by dislocation dynamics analyses system and we find the minimum resolved shear stress for the FR source to emit at least one closed loop. When the grain size is large enough compared to the size of FR source, the minimum resolved shear stress levels off to a certain value, but when the grain size is close to the size of the FR source, the minimum resolved shear stress shows a sharp increase. These results are modeled into the expression of the critical resolved shear stress of slip systems and continuum mechanics based crystal plasticity analyses of six-grained polycrystal models are made. Results of the crystal plasticity analyses show a distinct increase of macro- and microscopic yield stress for specimens with smaller mean grain diameter. Scale-dependent characteristics of the yield stress and its relation to some control parameters are discussed.     

Elastic/crystalline viscoplastic finite element analyses of single- and poly-crystal sheet deformations and their experimental verification

The elastic/crystalline viscoplastic constitutive equation, based on a newly proposed hardening-softening evolution equation, is introduced into the dynamic-explicit finite element code “Itas-Dynamic.” In the softening evolution equation, the effective distance and the angle between each slip system of a crystal are introduced to elucidate the interaction between the slip systems, which causes a decrease of dislocation density. The polycrystal sheet is modeled by Voronoi polygons, which correspond to the crystal grains; and by the selected orientations, which can relate to the texture, they are assigned to the integration points of the finite elements. We propose a direct crystal orientation assignment method, which means that each integration point of finite element has an assigned orientation, and its orientation can be rotated independently. Therefore, this inhomogeneous polycrystal model can consider the plastic induced texture development and subsequent anisotropy evolution. The parameters of the constitutive equation are identified by uni-axial tension tests carried out on single crystal sheets. Numerical results obtained for sheet tensions are compared with experimental ones to confirm the validity of our finite element code. Further, we investigate the following subjects: (1) how the initial orientation of single crystal affects slip band formation and strain localization; (2) how the grain size and particular orientations of the grain affect the strain localization in case of a polycrystal sheet. It is confirmed that the orientation of a single crystal can be related to the primary slip system and the deformation induced activation of that system, which in turn can be related to the slip band formation of the single crystal sheet. Further, in case of a polycrystal sheet, the larger the grain size, the more the strain localizes at a specific crystal, which has the particular orientation. It is confirmed through comparisons with experiments that our finite element code can predict the localization of strain in sheets and consequently can estimate the formability of sheet metals.     

A crystal plasticity model for strain-path changes in metals

A model is proposed that deals with the transient mechanical anisotropy during strain-path changes in metals. The basic mechanism is assumed to be latent hardening or softening of the slip systems, dependent on if they are active or passive during deformation, reflecting microstructural mechanisms that depend on the deformation mode rather than on the crystallography. The new model captures the experimentally observed behaviour of cross hardening in agreement with experiments for an AA3103 aluminium alloy. Generic results for strain reversals qualitatively agree with two types of behaviour reported in the literature – with or without a plateau on the stress–strain curve. The influence of the model parameters is studied through detailed calculations of the response of three selected parameter combinations, including the evolution of yield surface sections subsequent to 10% pre-strain. The mathematical complexity is kept to a minimum by avoiding explicit predictions related directly to underpinning microstructural changes. The starting point of the model is a combination of conventional texture and work hardening approaches, where an adapted full-constraints Taylor theory and a simple single-crystal work-hardening model for monotonic strain are used. However, the framework of the model is not restricted to these particular models.     

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.     

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.     

Experimental observations of evolving yield loci in biaxially strained AA5754-O

Experimental measurement of the plastic biaxial mechanical response for an aluminum alloy (AA5754-O) sheet metal is presented. Traditional methods of multiaxial sheet metal testing require the use of finite element analysis (FEA) or other assumptions of material response to determine the multiaxial true stress versus true strain behavior of the as-received sheet material. The method used here strives to produce less ambiguous measurements of data for a larger strain range than previously possible, through a combination of the Marciniak flat bottom ram test and an X-ray diffraction technique for stress measurement. The study is performed in conjunction with a study of the microstructural changes that occur during deformation, and these microstructural results are briefly mentioned in this work. Issues of calibration and applicability are discussed, and results are presented for uniaxial (U), plane strain (PS), and balanced biaxial (BB) extension. The results show repeatable behavior (within quantified uncertainties) for U to 20%, PS to almost 15%, and BB to above 20% in-plane strains. The results are first compared with three common yield locus models (von Mises’, Hill’48, and Hosford’79), and show some unexpected results in the shape change of the yield locus at high strain levels (>5% strain). These changes include the rotation of the locus toward the von Mises surface and elongation in the balanced biaxial direction. Comparison with a more complex yield locus model (Yld2000-2d with eight adjustable parameters) showed that the locus elongation in the biaxial direction could be fit well (for a specific level of work), but at the detriment of fit to the plane strain data. Artificially large plastic strain ratios would be needed to match both the biaxial and plane strain behavior even with this more complex model.     

Micromechanical modeling of the elasto-viscoplastic behavior of semi-crystalline polymers

A micromechanically based constitutive model for the elasto-viscoplastic deformation and texture evolution of semi-crystalline polymers is developed. The model idealizes the microstructure to consist of an aggregate of two-phase layered composite inclusions. A new framework for the composite inclusion model is formulated to facilitate the use of finite deformation elasto-viscoplastic constitutive models for each constituent phase. The crystalline lamellae are modeled as anisotropic elastic with plastic flow occurring via crystallographic slip. The amorphous phase is modeled as isotropic elastic with plastic flow being a rate-dependent process with strain hardening resulting from molecular orientation. The volume-averaged deformation and stress within the inclusions are related to the macroscopic fields by a hybrid interaction model. The uniaxial compression of initially isotropic high density polyethylene (HDPE) is taken as a case study. The ability of the model to capture the elasto-plastic stress-strain behavior of HDPE during monotonic and cyclic loading, the evolution of anisotropy, and the effect of crystallinity on initial modulus, yield stress, post-yield behavior and unloading-reloading cycles are presented.     

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.     

Scaling function, anisotropy and the size of RVE in elastic random polycrystals

This article is focused on the identification of the size of the representative volume element (RVE) in linear elastic randomly structured polycrystals made up of cubic single crystals. The RVE is approached by setting up stochastic Dirichlet and Neumann boundary value problems consistent with the Hill(-Mandel) macrohomogeneity condition. Within this framework we introduce a scaling function that relates the single crystal anisotropy to the scale of observation. We derive certain exact characteristics of the scaling function and postulate others based on detailed calculations on copper, lithium, tantalum, magnesium oxide and antimony-yttrium. In deriving the above, we make use of the fact that cubic crystals and polycrystals have a uniquely determined scale-independent bulk modulus. It turns out that the scaling function is exact in the single crystal anisotropy. A methodology to develop a material selection diagram that clearly separates the microscale from the macroscale is proposed. The proposed scaling function not only bridges the length scales but also unifies the treatment of a wide spectrum of cubic crystals. Although the scope of this article is restricted to aggregates made up of cubic-shaped and cubic-symmetry single crystals, the concept of the scaling function can be generalized to other crystal shapes and classes as well as to scaling of other elastic/inelastic properties.     

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.     

Characterization of yielding behavior of polycrystalline metals with single crystal plasticity based on representative characteristic length

Single crystal plasticity based on a representative characteristic length is proposed and introduced into a homogenization approach based on finite element analyses, which are applied to characterization of distinctive yielding behaviors of polycrystalline metals, yield-point elongation, and grain size strengthening. The computational manner for an implicit stress update is derived with the framework of a standard multi-surface plasticity at finite strain, where the evolution of the characteristic lengths are numerically converted from the accumulated slips of all of slip systems by exploiting the mathematical feature of the characteristic length as the intermediate function of the plastic internal variables. Furthermore, a constitutive model for a single crystal reproduces the stress–strain curve divided into three parts. Using two-scale finite element analysis, the macroscopic stress–strain response with yield-point elongation under a situation of low dislocation density is reproduced. Finally, the grain size effect on the yield strength is analyzed with modeling of the grain boundary in the context of the proposed constitutive model and is discussed from both macroscopic and microscopic views.     

Constitutive modeling for anisotropic/asymmetric hardening behavior of magnesium alloy sheets: Application to sheet springback

The constitutive model for the unusual asymmetric hardening behavior of magnesium alloy sheet presented in a companion paper (Lee, M.G., Wagoner, R.H., Lee, J.K., Chung, K., Kim, H.Y., 2008. Constitutive modeling for anisotropic/asymmetric hardening behavior of magnesium alloy sheet, Int. J. Plasticity 24(4), 545–582) was applied to the springback prediction in sheet metal forming. The implicit finite element program ABAQUS was utilized to implement the developed constitutive equations via user material subroutine. For the verification purpose, the springback of AZ31B magnesium alloy sheet was measured using the unconstrained cylindrical bending test of Numisheet (Numisheet ’2002 Benchmark Problem, 2002. In: Yang, D.Y., Oh, S.I., Huh, H., Kim, Y.H. (Eds.), Proceedings of 5th International Conference and Workshop on Numerical Simulation of 3D Sheet Forming Processes, Jeju, Korea) and 2D draw bend test. With the specially designed draw bend test the direct restraining force and long drawn distance were attainable, thus the measurement of the springback could be made with improved accuracy comparable with conventional U channel draw bend test. Besides the developed constitutive models, other models based on isotropic constitutive equations and the Chaboche type kinematic hardening model were also considered. Comparisons were made between simulated results by the finite element analysis and corresponding experiments and the newly proposed model showed enhanced prediction capability, which was also supported by the simple bending analysis adopting asymmetric stress–strain response.     

Gradient-based non-linear microstructure design

Fourier analysis is implemented on the orientation distribution of a polycrystalline microstructure. The linearity and convexity of the Fourier space, with respect to orientation, allows one to consider all possible distributions by considering all linear combinations of single-grain orientations. The limits of the Fourier space are therefore defined by the solutions to a set of linear programming problems. A unique approach to the linear programming, similar to the Krylov subspace methods for obtaining solutions to linear systems, is presented. The method is particularly efficient for this application where a large number of independent variables is often required. These solutions are then used as the constraints in the gradient-based optimization of non-linear functions within the Fourier space. In the example, Taylor yield theory and an anisotropic solution for the stress concentration around a hole in a plate of cubic-orthotropic polycrystalline material are expressed as non-linear functions within the Fourier space. The maximum obtainable ratio of Taylor factor to stress concentration for any polycrystalline orientation distribution in copper is found to be 1.22, more than double the minimum value.     

Homogenization of intergranular fracture towards a transient gradient damage model

This paper focuses on the intergranular fracture of polycrystalline materials, where a detailed model at the meso-scale is translated onto the macro-level through a proposed homogenization theory. The bottom-up strategy involves the introduction of an additional macro-kinematic field to characterize the average displacement jump within the unit cell. Together with the standard macro-strain field, the underlying processes are propagated onto the macro-scale by imposing the equivalence of power and energy at the two scales. The set of macro-governing equations and constitutive relations are next extracted naturally as per standard thermodynamics procedure. The resulting homogenized microforce balance recovers the so-called ‘implicit’ gradient expression with a transient nonlocal interaction. The homogenized gradient damage model is shown to fully regularize the softening behavior, i.e. the structural response is made mesh-independent, with the damage strain correctly localizing into a macroscopic crack, hence resolving the spurious damage growth observed in many conventional gradient damage models. Furthermore, the predictive capability of the homogenized model is demonstrated by benchmarking its solutions against reference meso-solutions, where a good match is obtained with minimal calibrations, for two different grain sizes.     

On a geometrically nonlinear damage model based on a multiplicative decomposition of the deformation gradient and the propagation of microcracks

We aim to derive a damage model for materials damaged by microcracks. The evolution of the cracks shall be governed by the maximum energy release rate, which was recently shown to be a direct consequence of the variational principle of a body with a crack (Arch. Appl. Mech. 69 (5) (1999) 337). From this, we get the path of the growing crack by introducing a series of thermodynamically equivalent straight cracks. The equivalence of the energy dissipated by microcrack growth and the damage dissipation leads to our damage evolution law. This evolution law will be embedded in a finite deformation framework based on a multiplicative decomposition into elastic and damage parts. As a consequence of this, we can present the anisotropic damaged elasticity tensor with the help of push and pull operations. The connection of this approach to other well known damage theories will be shown and the advantages of a finite element framework will be worked out. Numerical examples show the possibilities of the proposed model.     

Modeling aspects of low plastic strain amplitude multiaxial cyclic plasticity in conventional and ultrafine grain nickel

Aspects of the cyclically saturated responses of initially annealed, conventional grain size (average grain diameter of approximately 50 μm) and electrodeposited, ultrafine grain (grains from 20 to 500 nm) nickel to reversed proportional and 90° out-of-phase axial-torsional, strain-controlled cycling at a nominally constant equivalent inelastic strain amplitude of approximately 100 μ strain are reported. An anisotropic, axial-torsional subspace version of Abdel-Karim and Ohno’s kinematic hardening model is presented. Within the framework of conventional small strain, rate-independent plasticity, this approach is used to model the responses. An anhysteretic, phenomenolically based, magnetomechanical model is coupled to the rate-independent plasticity model to include the cyclic magnetostriction response. The kinematic hardening parameter determination scheme, using the proportional path responses, is described. The model correlations achieved are presented and the ability of the resulting models to capture the 90° out-of-phase responses is examined. The model parameter sets, as determined from the proportional responses, require small changes to result in more accurate correlation of the 90° out-of-phase responses and the implications of this are discussed. The relative values of the model parameters between the two materials reflect the initial microstructures. Persistent mean stresses associated with mean total strains imposed are successfully modeled for the proportional strain path responses but not for the 90° out-of-phase responses.