Assessment of crystal plasticity based calculation of the lattice spin of polycrystalline metals for FE implementation

When texture is incorporated in the finite element simulation of a metal forming process, much computer time can be saved by replacing continuous texture and corresponding yield locus updates by intermittent updates after strain intervals of e.g. 20%. The hypothesis that the evolution of the anisotropic properties of a polycrystalline material during such finite interval of plastic deformation can be modelled by just rotating the initial texture instead of continuously updating it by means of a polycrystal deformation model is tested in this work. Two spins for rotating the frame have been assessed: the classical rigid body spin and a crystal plasticity based “Mandel spin” (calculated from the rotated initial texture) which is the average of the spins of all the crystal lattices of the polycrystal. Each of these methods was used to study the evolution of the yield locus and the r-value distribution during the 20% strain interval. The results were compared to those obtained by simulating the texture evolution continuously using a polycrystal deformation model. When the texture was not updated during deformation, it was found that for most initial textures the Mandel spin does not perform better than the rigid body spin, except for some special initial textures for which the Mandel spin is much better. The latter ones are textures which are almost stable for the corresponding strain mode. When the texture was updated after each strain interval of e.g. 20% the Mandel spin performed much better than the rigid body spin.     

Microstructural effects on yield surface evolution in cubic metals using the viscoplastic ?-model

Polycrystalline yield surfaces of metals are a good way to characterize the anisotropy of plastic deformation. The evolution of these surfaces is impossible to accurately reproduce without taking into account the evolution of the material microstructure such as texture development. In this paper, a numerical computation of yield surfaces using the viscoplastic ?-model is proposed. Results concerning face-centered cubic metals subjected to a plane strain compression test are presented. The influence of several mechanical parameters (strain hardening, strain rate sensitivity coefficient and accumulated deformation) on subsequent yield surfaces evolution is studied. The analysis of the change in the shape and size of the yield surfaces shows that the results depend strongly on the parameter ? which controls the strength of the interactions in the polycrystal. In addition, the predictions are compared to the widely used viscoplastic self-consistent model as well as to experimental yield loci taken from the literature for various aluminum alloy sheets. A fairly good qualitative agreement between our ?-model results and the experimental ones is found. The probable links between the parameter ? and the microstructural features such as the stacking fault energy and the grain size of the polycrystal are also briefly discussed.     

Evaluation of self-consistent polycrystal plasticity models for magnesium alloy AZ31B sheet

Various self-consistent polycrystal plasticity models for hexagonal close packed (HCP) polycrystals are evaluated by studying the deformation behavior of magnesium alloy AZ31B sheet under different uniaxial strain paths. In all employed polycrystal plasticity models both slip and twinning contribute to plastic deformation. The material parameters for the various models are fitted to experimental uniaxial tension and compression along the rolling direction (RD) and then used to predict uniaxial tension and compression along the traverse direction (TD) and uniaxial compression in the normal direction (ND). An assessment of the predictive capability of the polycrystal plasticity models is made based on comparisons of the predicted and experimental stress responses and R values. It is found that, among the models examined, the self-consistent models with grain interaction stiffness halfway between those of the limiting Secant (stiff) and Tangent (compliant) approximations give the best results. Among the available options, the Affine self-consistent scheme results in the best overall performance. Furthermore, it is demonstrated that the R values under uniaxial tension and compression within the sheet plane show a strong dependence on imposed strain. This suggests that developing anisotropic yield functions using measured R values must account for the strain dependence.     

Forming simulation of aluminum sheets using an anisotropic yield function coupled with crystal plasticity theory

This paper describes the application of a coupled crystal plasticity based microstructural model with an anisotropic yield criterion to compute a 3D yield surface of a textured aluminum sheet (continuous cast AA5754 aluminum sheet). Both the in-plane and out-of-plane deformation characteristics of the sheet material have been generated from the measured initial texture and the uniaxial tensile curve along the rolling direction of the sheet by employing a rate-dependent crystal plasticity model. It is shown that the stress–strain curves and R-value distribution in all orientations of the sheet surface can be modeled accurately by crystal plasticity if a “finite element per grain” unit cell model is used that accounts for non-uniform deformation as well as grain interactions. In particular, the polycrystal calculation using the Bassani and Wu (1991) single crystal hardening law and experimental electron backscatter data as input has been shown to be accurate enough to substitute experimental data by crystal plasticity data for calibration of macroscopic yield functions. The macroscopic anisotropic yield criterion CPB06ex2 (Plunkett et al., 2008) has been calibrated using the results of the polycrystal calculations and the experimental data from mechanical tests. The coupled model is validated by comparing its predictions with the anisotropy in the experimental yield stress ratio and strain ratios at 15% tensile deformation. The biaxial section of the 3D yield surface calculated directly by crystal plasticity model and that predicted by the phenomenological model calibrated with experimental and crystal plasticity data are also compared. The good agreement shows the strength of the approach. Although in this paper, the Plunkett et al. (2008) yield function is used, the proposed methodology is general and can be applied to any yield function. The results presented here represent a robust demonstration of implementing microscale crystal plasticity simulation with measured texture data and hardening laws in macroscale yield criterion simulations in an accurate manner.     

Mechanical response and texture evolution of AZ31 alloy at large strains for different strain rates and temperatures

In order to study the behavior of material under finite deformation at various strain rates, the responses of AZ31 Mg sheet are measured under uniaxial (tension and compression) and multiaxial (simple shear) loadings along rolling direction (RD), 45° to rolling direction (DD), 90° to rolling direction (TD), and normal to the sheet (ND) to large strains. The material exhibits positive strain rate sensitivity (SRS) at room and elevated temperatures; the SRS is more pronounced at high temperatures and lower strain rates. The r-value of the material under tensile loading at room temperatures is higher in TD at lower strain rate. Texture measurements on several failed specimens are reported under tension and simple shear after finite plastic deformation of about 20% equivalent strain. The as-received material exhibits a strong fiber with equal fractions of grains having the c-axis slightly tilted away from the sheet normal towards both +RD and −RD. Pole figures obtained after tensile loading along the rolling direction (RD) show that the texture of the material strengthens even at low strains, with c-axis perpendicular to the sheet plane and prism planes lining up in a majority of grains. However, the tensile loading axis along TD does not lead to similar texture strengthening; the c-axis distribution appears to be virtually unchanged from the virgin state. The pole figures obtained after in-plane compression along RD brings the c-axes of the grains parallel to the loading direction. The pole figures after simple shear loading show that the c-axis rotates to lie on the sheet plane consistent with a compression axis 45° away on the sheet plane.     

Texture gradient simulations for extrusion and reversible rolling of FCC metals

A texture simulation method is described for some complex plane strain deformation paths during hot shaping of FCC metals. The method employs both finite element calculations and a polycrystal plasticity model based on the Relaxed-Constraints (RC) Taylor hypothesis and a viscoplastic constitutive law. We have considered the {111}<110> slip systems and the {100}, {110}, {112} <110> non-octahedral slip systems. The finite element codes simulate the strain paths of material flow during a shaping process. The local velocity gradients, expressed in the macroscopic reference coordinates, are rewritten in the local flow line coordinates using a kinematic analysis for steady-state flow. Secondly, for the different deformation paths, the RC polycrystal plasticity model is used to numerically simulate the local deformation texture evolutions as a function of depth. Texture simulations are carried out for two deformation processes combining hot compression and shear: extrusion and reversible rolling. For extrusion, the simulated pole figures and ODFs show the typical texture variations through the thickness of an extruded 6082 aluminium alloys, i.e. (β-fibre in the centre and a TD rotated copper component near the surface. It is shown that hot reversible rolling should develop a strong pure shear texture {001}<110> near the surface.     

Modeling lattice strain evolution at finite strains and experimental verification for copper and stainless steel using in situ neutron diffraction

An elasto-plastic self-consistent (EPSC) polycrystal model is extended to account, in an approximate fashion, for the kinematics of large strains, rigid body rotations, texture evolution and grain shape evolution. In situ neutron diffraction measurements of the flow stress, internal strain, texture and diffraction peak intensity evolutions were performed on polycrystalline copper and stainless steel, up to true tensile strains of ε = 0.3. Suitably adjusted slip system hardening model parameters enable the model to quantitatively describe the flow stress of the polycrystalline aggregate. Quantitative predictions of the texture evolution and the internal strain evolution along the stress axis are good, while predictions of transverse internal strains (perpendicular to the tensile loading direction) are less satisfactory. The latter exhibit a large dispersion from grain to grain around a macroscopic average, and the implications of this finding for the interpretation of in situ neutron diffraction method are explored. Finally, as a demonstration of the applicability of the model to problems involving finite rotation, as well as deformation, simulations of simple shear were conducted which predict a texture evolution in agreement with published experimental data, and other modeling approaches as well.     

Effect of a prestrain on the subsequent yielding of low carbon steel sheets: experiments and simulations

A sequence of uniaxial tension experiments have been performed on low carbon steel sheets along different axes and have measured the evolution of the subsequent yield stress after several moderate prestrain levels and for various angles between the two tensile axes. Simulations have been conducted in order to evaluate the respective influence of the development of deformation textures and of intragranular microstructures. The results of simulations have shown the major influence of microstructural effects and the importance of the initial texture rather than its evolution. In agreement with experiments, one finds an amplitude of 20 percent and a maximum effect around 45 degrees for the relative variation of yield stress with the angles between the two tensile axes at a prestrain of 0.15.     

The effect of large prestrains on fatigue

Evaluation of the fatigue behavior of stamped parts requires a knowledge of the effects of large prestrains on fatigue. This paper reports the results of a test program carried out to investigate the effect of tensile and rolling prestrains on the fatigue strength of SAE 1010 and Maxiform 50 steels. Large rolling-induced prestrains were shown to cause a significant increase, while large tensile prestrains (up to 0.9 times the true fracture strain of the steel) caused a moderate increase in the long-life fatigue properties of these steels. Still larger tensile prestrains (0.95 times the true fracture strain) caused the fatigue strength to reduce to its initial value.     

Elastic–plastic behavior of steel sheets under in-plane cyclic tension–compression at large strain

Elastic–plastic behavior of two types of steel sheets for press-forming (an aluminum-killed mild steel and a dual-phase high strength steel of 590 MPa ultimate tensile strength) under in-plane cyclic tension–compression at large strain (up to 25% strain for mild steel and 13% for high strength steel) have been investigated. From the experiments, it was found that the cyclic hardening is strongly influenced by cyclic strain range and mean strain. Transient softening and workhardening stagnation due to the Bauschinger effect, as well as the decrease in Young's moduli with increasing prestrain, were also observed during stress reversals. Some important points in constitutive modeling for such large-strain cyclic elasto-plasticity are discussed by comparing the stress–strain responses calculated by typical constitutive models of mixed isotropic–kinematic hardening with the corresponding experimental observations.     

Basic equations for the microscopic theory of plasticity of polycrystalline materials with given texture

An expression for the yield stress of anisotropic materials is applied to the anisotropic strength of hard rolled copper foils whose crystallographic texture is known. We assume that this crystallographic texture is the only cause of the anisotropic plastic behaviour of the material. The model used for the yield stress is also used to deduce:

1. Stress-strain relations for isotropic polycrystalline materials;
2. A formula for the fully plastic strain tensor, applied to anisotropic hard rolled copper foils.

For the anisotropic copper foils considered the calculated curves of the yield stress and of the strain tensor as a function of the angle x between rolling and tensile direction agree qualitatively with the measured values. However, the theory is not complete, since the yield stress and the plastic strain tensor are both a function of a parameter Q, the fraction of the number of available crystallographic slip planes on which the maximum shear stress has reached the critical value τ_a. We assume that for “fully” plastic deformation a certain critical fraction Q _e of the total number of slip planes has to be active. The fraction Q _e is called the critical active quantity. With the parameter Q _e we adjust the calculated curves to the measured ones. The dependence of Q _e on the properties of the material (e.g. the crystallographic texture) is discussed in Appendix I.     

Shear banding in strain hardening polycrystals during rolling

Initiation and development of shear band (SB) in f.c.c. strain hardening polycrystals during rolling are modelled in terms of crystallographic texture. The constitutive law of the material is expressed in terms of the texture-dependent normalized yield surface and the critical shear stress which evolves with strain. The normalized yield surface is predicted by the Taylor model as a function of rolling texture. It is shown that a rounded vertex (RV) develops at the loading point as the rolling texture becomes more and more marked. A detailed characterization of the RV is carried out. It is found that the normalized curvature radius of the RV decreases from unity towards zero at very large strain. This allows for a small stress perturbation to induce a shear strain perturbation with a large orientation deviation of deformation. By linearized stability analysis, the condition for initiation of SB from the shear strain perturbation is obtained. Development of SB is analysed by solving the established governing equations of shear banding. It is shown that the conditions for SB initiation and saturation of shear localisation depend strongly on the texture. Based on this model problem, a long discussion is carried out which allows a better understanding of the basic physical origin and saturation of SB in strain-hardening polycrystals, as well as the effects of yield surface curvature and yield surface rotation whose general form is derived.     

Experimental investigation of initial and subsequent yield surfaces for laminated metal matrix composites

The aim of this work is to construct yield surfaces to describe initial yielding and characterize hardening behavior of a highly anisotropic material. A methodology for constructing yield surfaces for isotropic materials using axial–torsion loading is extended to highly anisotropic materials. The technique uses a sensitive definition of yielding based on permanent strain rather than offset strain, and enables multiple yield points and multiple yield surfaces to be conducted on a single specimen. A target value of 20 × 10⁻⁶ is used for Al₂O₃ fiber reinforced aluminum laminates having a fiber volume fraction of 0.55. Sixteen radial probes are used to define the yield locus in the axial–shear stress plane. Initial yield surfaces for [0₄], [90₄], and [0/90]₂ fibrous aluminum laminates are well described by ellipses in the axial–shear stress plane having aspect ratios of 10, 2.5, and 3.3, respectively. For reference, the aspect ratio of the Mises ellipse for an isotropic material is 1.73. Initial yield surfaces do not have a tension–compression asymmetry. Four overload profiles (plus, ex, hourglass, and zee) are applied to characterize hardening of a [0/90]₂ laminate by constructing 30 subsequent yield surfaces. Parameters to describe the center and axes of an ellipse are regressed to the yield points. The results clearly indicate that kinematic hardening dominates so that material state evolution can be described by tracking the center of the yield locus. For a nonproportional overload of (σ, τ) = (500, 70) MPa, the center of the yield locus translated to (σ, τ) = (430, 37) MPa and the ellipse major axis was only 110 MPa.     

Strain-rate history and temperature effects on the torsional-shear behavior of a mild steel

Previous investigations on the effects of strain-rate and temperature histories on the mechanical behavior of steel are briefly reviewed. A study is presented on the influence of strain rate and strain-rate history on the shear behavior of a mild steel, over a wide range of temperature Experiments were performed on thin-walled tubular specimens of short gage length, using a torsional split-Hopkinson-bar apparatus adapted to permit quasi-static as well as dynamic straining at different temperatures. The constant-rate behavior was first measured at nominal strain rates of 10^?3 and 10³ s^?1 for ?150, ?100, ?50, 20, 200 and 400°C. Tests were then carried out, at the same temperatures, in which the strain rate was suddenly increased during deformation from the lower to the higher rate at various large values of plastic strain. The increase in rate occurred in a time of the order of 20 μs so that relatively little change of strain took place during the jump. The low strain-rate results show a well-defined elastic limit but no yield drop, a small yield plateau is found at room temperature. The subsequent strain hardening shows a maximum at 200°C, when serrated flow occurs and the ductility is reduced. The high strain-rate results show a considerable drop of stress at yield. The post-yield flow stress decreases steadily with increasing temperature, throughout the temperature range investigated. At room temperature and below, the strain-hardening rate becomes negative at large strains. The adiabatic temperature rise in the dynamic tests was computed on the assumption that the plastic work is entirely converted to heat. This enabled the isothermal dynamic stress-strain curves to be calculated, and showed that considerable thermal softening took place. The initial response to a strain-rate jump is approximately elastic, and has a magnitude which increases with decrease of testing temperature; it is little affected by the amount of prestrain. At 200 and 400° C, a yield drop occurs after the initial stress increment. The post-jump flow stress is always greater than that for the same strain in a constant-rate dynamic test, the strain-hardening rate becoming negative at large strains or low testing temperature. This observed effect of strain-rate history cannot be explained by the thermal softening accompanying dynamic deformation. These and other results concerning total ductility under various strain-rate and temperature conditions show that strain-rate history strongly affects the mechanical behavior of the mild steel tested and, hence, should be taken into account in the formulation of constitutive equations for that material.     

Mechanics of nonbuckling interconnects with prestrain for stretchable electronics

The performance of the flexibility and stretchability of flexible electronics depends on the mechanical structure design, for which a great progress has been made in past years. The use of prestrain in the substrate, causing the compression of the transferred interconnects, can provide high elastic stretchability. Recently, the nonbuckling interconnects have been designed, where thick bar replaces thin ribbon layout to yield scissor-like in-plane deformation instead of in-or out-of-plane buckling modes. The nonbuckling interconnect design achieves significantly enhanced stretchability. However,combined use of prestrain and nonbuckling interconnects has not been explored. This paper aims to study the mechanical behavior of nonbuckling interconnects bonded to the prestrained substrate analytically and numerically. It is found that larger prestrain,longer straight segment, and smaller arc radius yield smaller strain in the interconnects.On the other hand, larger prestrain can also cause larger strain in the interconnects after releasing the prestrain. Therefore, the optimization of the prestrain needs to be found to achieve favorable stretchability.     

High Strain Rate Torsion Properties of Ultrafine-Grained Aluminum

Mechanical properties of most metallic materials can be improved by reducing their grain size. One of the methods used to reduce the grain size even to the nanometer level is the severe plastic deformation processing. Equal Channel Angular Pressing (ECAP) is one of the most promising severe plastic deformation processes for the nanocrystallization of ductile metals. Nanocrystalline and ultrafine grained metals usually have significantly higher strength properties but lower tensile ductility compared to the coarse grained metals. In this work, the torsion properties of ECAP processed ultrafine grained pure 1070 aluminum were studied in a wide range of strain rates using both servohydraulic materials testing machines and Hopkinson Split Bar techniques. The material exhibits extremely high ductility in torsion and the specimens did not fail even after 300% of strain. Pronounced yield point behavior was observed at strain rates 500 s⁻¹ and higher, whereas at lower strain rates the yielding was continuous. The material showed slight strain softening at the strain rate of 10⁻⁴ s⁻¹, almost ideally plastic behavior at strain rates between 10⁻³ s⁻¹ and 500 s⁻¹, and slight but increasing strain hardening at strain rates higher than that. The tests were monitored using digital cameras, and the strain distributions on the surface of the specimens were calculated using digital image correlation. The strain in the specimen localized very rapidly after yielding at all strain rates, and the localization lead to the development of a shear band. At high strain rates the shear band developed faster than at low strain rates.     

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.     

Anisotropic yield criterion for polycrystalline metals using texture and crystal symmetries

An anisotropic yield criterion for polycrystalline metals which uses texture data and takes advantage of crystal symmetries is presented. A linear transformation is developed to map an anisotropic yield surface for a polycrystal to an appropriate isotropic yield surface. The transformation developed reflects the symmetry of the material being modeled. First, the transformation is determined. Then, information regarding the orientation distribution (texture) of the crystals in a polycrystalline aggregate is used to determine, via averaging, the transformation for the polycrystal. The transformation, along with appropriate isotropic yield surface, provides a phenomenological approach to modeling yield, yet accounts for microstructural texture. The approach reduces to the Hill (1950) anisotropic plasticity theory under certain conditions. The yield surfaces and R-values for various face-centered-cubic ( fcc) polycrystalline textures are computed by this method. Results compare favorably with those given by other theories, and with experiment. The method proves to have the computational efficiency of phenomenological approaches to modeling yield, while effectively incorporating the physics of more complex crystallographic approaches.     

Controlling the yield behavior of fat-oil mixtures using cooling rate

The rheological behavior such as yielding of fat crystal networks are dictated by many variables. Among these variables, the shape of the constituent fat cluster is important yet relatively unexplored. In this work, we describe the rheological investigations of a fat-oil system which can be formulated to either contain bundles of needles or spherical clusters by controlling the cooling rate and fat concentration. Fat-oil mixtures containing high-fat concentrations exhibited weak frequency dependence of storage modulus (G ^′) and loss modulus (G ^″). The yielding behavior of the mixtures were investigated by large amplitude oscillatory shear (LAOS) rheology using strain and stress controlled modes. Lissajous-Bowditch plots and orthogonal set of Chebyshev polynomials were used to analyze the non-linearities associated with the yielded fat-oil mixtures. For a given fat concentration, the yield stress of fat networks obtained at low cooling rates (bundles of needles) were similar to that of networks obtained at high cooling rates (spherical clusters). However, after yielding, Lissajous-Bowditch plots suggested that the mixtures comprising of bundles of needles exhibited viscous-like behavior while the spherical clusters exhibited a plastic-like behavior. This was further supported by microscopy images of yielded fat-oil mixtures. Overall, for a given fat concentration, the two different shapes of fat clusters can give rise to networks of similar yield stress values but different behaviors after yielding.