A modified model for simulating latent hardening during the plastic deformation of rate-dependent FCC polycrystals

A strain hardening model for the plastic deformation of rate-dependent FCC crystals is proposed based on experimental observations previously reported for single crystals. This model, which is an extension of that employed by et al. [1983], includes both the self-hardening and latent hardening of the slip systems. The differential hardening of the latent systems is assumed to arise from the interaction between glide dislocations and forests. With this hardening model and a rate-sensitive crystal plasticity theory, the deformation behavior of FCC polycrystals can be predicted from the deformation response of the constituent single crystals. As examples, the uniaxial tensile behaviour of pure aluminum and copper polycrystals is simulated using the extended model, and the results are compared with published experimental data. The effects of latent hardening on polycrystal deformation, especially on flow stress and the formation of tensile textures, are discussed.     

Study of limit strains for FCC and BCC sheet metal using polycrystal plasticity

In this research, we analyze forming-limit strains of FCC and BCC materials using a viscoplastic self-consistent polycrystal model (VPSC) in conjunction with the Marciniak–Kuczynski (MK) approach. In particular, our work is focused on the theoretical analysis and comparison between FCC and BCC crystal structures made by Inal et al. [Inal, K., Neale, K.W., Aboutajeddine, A., 2005. Forming limit comparison for FCC and BCC sheets, International Journal of Plasticity, 21, 1255–1266]. These authors performed their simulations based on a generalized Taylor-type polycrystal model (MK-FC), finding a remarkably low forming-limit curve for the FCC material and an extremely high forming-limit curve for the BCC material, in the biaxial stretching range. We verified that our predictions are similar to Inal’s results for both FCC and BCC materials when the MK-FC model is used. However, MK-VPSC calculations do not give such extreme values, and we believe that this theory predicts much more reliable results for both FCC and BCC crystallographic assumptions. We also found that localized necking depends on texture evolution in the vicinity of equi-biaxial stretching, through the sharpness of the predicted yield surface. Finally, it is shown that the MK-VPSC’s predictions are in good agreement with experimental data for AA5182-O and a DQ-type steel-sheet metal.     

Non-local crystal plasticity model with intrinsic SSD and GND effects

A strain gradient-dependent crystal plasticity approach is presented to model the constitutive behaviour of polycrystal FCC metals under large plastic deformation. In order to be capable of predicting scale dependence, the heterogeneous deformation-induced evolution and distribution of geometrically necessary dislocations (GNDs) are incorporated into the phenomenological continuum theory of crystal plasticity. Consequently, the resulting boundary value problem accommodates, in addition to the ordinary stress equilibrium condition, a condition which sets the additional nodal degrees of freedom, the edge and screw GND densities, proportional (in a weak sense) to the gradients of crystalline slip. Next to this direct coupling between microstructural dislocation evolutions and macroscopic gradients of plastic slip, another characteristic of the presented crystal plasticity model is the incorporation of the GND-effect, which leads to an essentially different constitutive behaviour than the statistically stored dislocation (SSD) densities. The GNDs, by their geometrical nature of locally similar signs, are expected to influence the plastic flow through a non-local back-stress measure, counteracting the resolved shear stress on the slip systems in the undeformed situation and providing a kinematic hardening contribution. Furthermore, the interactions between both SSD and GND densities are subject to the formation of slip system obstacle densities and accompanying hardening, accountable for slip resistance. As an example problem and without loss of generality, the model is applied to predict the formation of boundary layers and the accompanying size effect of a constrained strip under simple shear deformation, for symmetric double-slip conditions.     

Anisotropic hardening in single crystals and the plasticity of polycrystals

Based on the dislocation structures developed during plastic deformation, an anisotropic hardening law is developed to describe the latent hardening behavior of slip systems under multislip. This theory incorporates the concept of isotropic hardening, kinematic hardening, and the two-parameter representation; it automatically includes the strength differential between the forward and reversed slips and between the acute and obtuse cross slips. The self-hardening modulus of a slip system is found to be “associated” with the latent hardening law involved, and, based on some experimental evidence, two specific sets of self-hardening modulus are suggested. An important feature of this associated modulus is that the slip system with a soft latent hardening (e.g., the reversed system with a Bauschinge effect) will have an enhanced self-hardening modulus. This newly developed hardening law, together with its associated latent hardening moduli, is then applied to examine the strain-hardening behavior of a polycrystal. Although crystals with a stronger latent hardening will, in general, also lead to a stranger strain-hardening for the polycrystal, the stress-strain behavior of the polycrystal using the kinematic hardening law of single crystals is found to be not necessarily softer than that using the isotropic hardening law. Within the range of experimentally measured latent hardening ratio of slip systems, the anisotropic theory is also used to calculate the motion of yield surface of a polycrystal. The general results, employing four selected types of anisotropic hardening, all show the essential features of experimental observations by Phillips and his co-workers. The application is highlighted with a reasonably successful quantitative modeling of initial and subsequent yield surfaces of an aluminum.     

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.     

Micro–macro approach for the rock salt behaviour

Rock salt is considered as a pure aggregate of halite (mineral NaCl) crystals and its behaviour is investigated by a micro–macro approach. The behaviour of the polycrystalline aggregate is deduced from the properties of the constituent halite crystals. A model for the elastoplastic behaviour of halite crystal has been deduced from experimental data available in the literature. The basic equations of the micro–macro model for the polycrystalline medium and the calculation method are then presented and the elastoplastic behaviour of rock salt is investigated by this method. The hardening effects obtained for the polycrystal are found to be very different from those obtained for FCC metal polycrystals. The differences are explained as a consequence of differences of families of glide systems in these crystals. Finally, the internal stresses in the polycrystal are studied in order to elucidate the origin of cracking and damage of the rock salt.     

Predictions of forming limit diagrams using a rate-dependent polycrystal self-consistent plasticity model

In the sheet-metal forming industry, forming-limit strains have been a useful tool for quantifying metals formability. However, the experimental measurement of these strains is a difficult, time consuming and expensive process. It would be useful if strains calculated with a theoretical model could replace many of the experimental measurements. In this research, we analyze forming-limit strains of metals using a rate-dependent plasticity, polycrystal, self-consistent (VPSC) model in conjunction with the Marciniak–Kuczynski (M–K) approach. Previous researchers have studied forming limit diagrams (FLDs) based on the full-constraints Taylor model. This is the first time, to the authors’ knowledge, that the self-consistent approach has been introduced to simulate the polycrystal FLD behavior. Numerous microstructural factors characterizing the material have a strong influence on the FLD, so our model includes the effects of slip hardening, strain-rate sensitivity, anisotropy and initial texture. Finally, the calculation of the FLD with a more realistic scale transition successfully predicts some of the experimental tendencies that the Taylor model cannot reproduce for aluminum alloys AA6116-T4 and AA5182-O.     

PREDICTION OF EARING IN TEXTURED ALUMINIUM SHEETS BASED ON CRYSTAL PLASTICITY

The phenomenon of earing is investigated in the present study based on the theory of crystal plasticity with the dynamic explicit finite element program developed. Firstly texture analysis is carried out of rolled aluminium alloy Al5052 by means of X-ray technique. Then from the texture coefficients an analytical expression for the orientation distribution function (ODF) is derived making use of the computer algebraic language Mathematica4.0, which makes it easier to discretize the ODF into a series of Eulerian angles representing the distribution of lattices and further the preferred orientation (texture) of crystals of the original sheets. For the polycrystal model, the material is described using crystal plasticity where each material point in grains with each grain modelled as an FCC crystal with 12 distinct slip systems. The modified Taylor theory of crystal plasticity is used and only the initial texture is taken into consideration during large plastic deformation. Numerical simulation of earing has been performed for an aluminium sheet with texture and one with crystals exhibiting random distribution to demonstrate the effect of texture of materials on their plastic anisotropy and formability. Project supported by the National Natural Science Foundation of China (No. 59875025).     

Rate-independent crystalline and polycrystalline plasticity,application to FCC materials

This paper deals with the simulation of the mechanical response and texture evolution of cubic crystals and polycrystals for a rate-independent elastic–plastic constitutive law. No viscous effects are considered. An algorithm is introduced to treat the difficult case of multi-surface plasticity. This algorithm allows the computation of the mechanical response of a single crystal. The corresponding yield surface is made of the intersection of several hyper-planes in the stress space. The problem of the multiplicity of the slip systems is solved thanks to a pseudo-inversion method. Self and latent hardening are taken into account. In order to compute the response of a polycrystal, a Taylor homogenization scheme is used. The stress–strain response of single crystals and polycrystals is computed for various loading cases. The texture evolution predicted for compression, plane strain compression and simple shear are compared with the results given by a visco-plastic polycrystalline model.     

Strain hardening in 2D discrete dislocation dynamics simulations: A new ‘2.5D’ algorithm

The two-dimensional discrete dislocation dynamics (2D DD) method, consisting of parallel straight edge dislocations gliding on independent slip systems in a plane strain model of a crystal, is often used to study complicated boundary value problems in crystal plasticity. However, the absence of truly three dimensional mechanisms such as junction formation means that forest hardening cannot be modeled, unless additional so-called ‘2.5D’ constitutive rules are prescribed for short-range dislocation interactions. Here, results from three dimensional dislocation dynamics (3D DD) simulations in an FCC material are used to define new constitutive rules for short-range interactions and junction formation between dislocations on intersecting slip systems in 2D. The mutual strengthening effect of junctions on preexisting obstacles, such as precipitates or grain boundaries, is also accounted for in the model. The new ‘2.5D’ DD model, with no arbitrary adjustable parameters beyond those obtained from lower scale simulation methods, is shown to predict athermal hardening rates, differences in flow behavior for single and multiple slip, and latent hardening ratios. All these phenomena are well-established in the plasticity of crystals and quantitative results predicted by the model are in good agreement with experimental observations.     

Simulation of the behaviour of FCC polycrystals during reversed torsion

Taylor-type polycrystal plasticity models with various single slip hardening laws are evaluated by studying the large strain behaviour of FCC polycrystals during reversed torsion. The hardening laws considered include the model of Asaro and Needleman (“Texture Development and Strain Hardening in Rate Dependent Polycrystals,” Acta Metall. (1985), 34, 1553) as well as a power-law and an exponential version of that, and a more recent model by Bassani and Wu (“Latent Hardening in Single Crystals II. Analytical Characterization and Predictions,” Proc. R. Soc. Lond. (1991), A435, 21). The material parameters for the various hardening laws are fitted to experimental compression data for copper and then used to predict reversed large strain torsion of tubes. Differences under “free-end” (axially unconstrained) or “fixed-end” (axially constrained) conditions between predictions and experimental observations are discussed in detail. In addition to the torque response, the Swift effects upon twist reversal are studied.     

NUMERICAL SIMULATION OF NON-PLANAR PLASTIC ANISOTROPY OF COLD ROLLING POLYCRYSTALLINE MATERIALS

The plastic anisotropy of sheet metal is usually caused by preferred orientation of grains, developed by mechanical deformation and thermal treatment. In the present study, a Taylor-like polycrystal model suggested by Asaro and Needleman is applied to investigate the evolution of the anisotropic behavior of a face centered cubic (FCC) polycrystalline metal, which is considered having {111} (110) slip systems, by stretching it along an arbitrary direction after it has undergonea plane-strata compression that rationally simulates the cold rolling process of FCC polycrystalline pure aluminium. By using the Taylor-like polycrystal model, pole figures are obtained to describe the texture development of polycrystalline aggregate after plane-strain compression, and then the plastic anisotropy of polycrystalline aggregate is evaluated by stretching the polycrystalline aggregate in different direction in term of yield stress. According to the results, the contours of longitudinal flow stress in three-dimensional orientation space are given and analyzed. Experiment results similar to the prediction of planar anisotropy can be found inthe literature written by Takahashi et al. that in directly show the correctness of the prediction of non-planar plastic anisotropy by this analysis.     

Transition theories of elastic-plastic deformation of metallic polycrystals

Summary A general approach to the problem of determination of elastoplastic behavior of metallic polycrystals at finite deformation is presented. The relation between moving dislocation density and global slip rate for grains is developed. Transition to grain response is obtained by introducing the hardening matrix. Field equations for heterogeneous elastoplastic metals are transformed into an integral equation, using Green functions technique. This allows to find the spin of the lattice related to texture formation.Scale transition is achieved by a self-consistent approximation of the integral equation. New results concerning BCC metals (sheet steel) are presented. They apply to tensile test, Lankford coefficient, initial and subsequent yield surfaces, and evolution of the internal state of the polycrystal: second-order residual stress, stored energy and texture evolution.     

Inelastic behaviour of polycrystalline metals under complex loading condition

In this paper, a simple mechanical model of polycrystal is suggested on the basis of the crystal plasticity. This model incorporates interactions among grains in polycrystal and interactions among slip systems in a crystal grain component. Especially, an equation for a slip strain of slip system is based on a theory of thermo-activation motion of dislocation. In this way, the model can take into account rate (creep) properties and temperature effects as well as an induced plastic anisotropy depending on deformation path, under multi-axial condition.Some computations are done by using the model for typical sets of strain paths and compared with corresponding experimental results of aluminum alloy at 200°C. The computational results reproduce accurately typical features observed in the experiments.     

Predictions of plastic stress-strain relations of polycrystals based on the Lin model

A uniaxial stress-strain curve is calculated for a polycrystal consisting of numerous nonhardening FCC crystals based on the Lin model. The coefficient matrix relating the deviation of plastic strain to the one of stress in the Lin model turned out to be singular. The simplex method is used to determine slips of all slip systems in constituent crystals. The calculated curves lie between the KBW model and the Hill model.     

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.     

Model for plasticity effects in metals under nonproportional cyclic loading

A possible physical mechanism for additional hardening is proposed on the basis of an analysis of experiments on nonproportional cyclic loading of metals. A model for an elastoplastic polycrystal with a hardening law taking into account the interaction of slip systems is developed. The effect of additional hardening for elliptic strain paths and the shapes of stress paths and hysteresis loops typical of elliptic strain paths are described qualitatively. A violation of the assumption of the local determinacy and orientations of stress paths is considered for the square strain paths taking place in tests of chromium-nickel austenite stainless steels. Perm' State Technical University, Perm', 614600. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 40, No. 6, pp. 144–151, November–December, 1999.     

Nanocontact between BCC tungsten and FCC nickel using the quasicontinuum method

The nanocontact between BCC tungsten and FCC nickel is investigated by using the multiscale simulation method (quasicontinuum method). The new proposed interatomic potential – Extended Finnis-Sinclair potential is adopted and compared with the conventional Ackland Finnis-Sinclair potential in the nanoindentation model. A series of multiscale simulations of contact formation, indentation, subsequent pulling, and separation between a tungsten indenter and a nickel substrate are performed by means of the quasicontinuum method. Different plastic mechanisms are observed for different crystal directions and indenter shapes. Besides the adherence between the tungsten indenter and the nickel substrate, the geometrical shape, which is formed by the slip planes inside the substrate, plays a vital role in the material transfer between the indenter and the substrate. The formation and evolution of necking and cracking in the nanocontact model are also discussed.     

Wedge indentation into elastic-plastic single crystals, 1: Asymptotic fields for nearly-flat wedge

Asymptotic stress and deformation fields under the contact point singularities of a nearly-flat wedge indenter and of a flat punch are derived for elastic ideally-plastic single crystals with three effective in-plane slip systems that admit a plane strain deformation state. Face-centered cubic (FCC), body-centered cubic (BCC), and hexagonal-close packed (HCP) crystals are considered. The asymptotic fields for the flat punch are analogous to those at the tip of a stationary crack, so a potential solution is that the deformation field consists entirely of angular constant stress plastic sectors separated by rays of plastic deformation across which stresses change discontinuously. The asymptotic fields for a nearly-flat wedge indenter are analogous to those of a quasistatically growing crack tip fields in that stress discontinuities can not exist across sector boundaries. Hence, the asymptotic fields under the contact point singularities of a nearly-flat wedge indenter are significantly different than those under a flat punch. A family of solutions is derived that consists entirely of elastically deforming angular sectors separated by rays of plastic deformation across which the stress state is continuous. Such a solution can be found for FCC and BCC crystals, but it is shown that the asymptotic fields for HCP crystals must include at least one angular constant stress plastic sector. The structure of such fields is important because they play a significant role in the establishment of the overall fields under a wedge indenter in a single crystal. Numerical simulations—discussed in detail in a companion paper—of the stress and deformation fields under the contact point singularity of a wedge indenter for a FCC crystal possess the salient features of the analytical solution.