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.     

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.     

Thermally actuated shape-memory polymers: Experiments, theory, and numerical simulations

With the aim of developing a thermo-mechanically coupled large-deformation constitutive theory and a numerical-simulation capability for modeling the response of thermally actuated shape-memory polymers, we have (i) conducted large strain compression experiments on a representative shape-memory polymer to strains of approximately unity at strain rates of 10⁻³ and 10⁻¹ s⁻¹, and at temperatures ranging from room temperature to approximately 30 °C above the glass transition temperature of the polymer; (ii) formulated a thermo-mechanically coupled large-deformation constitutive theory; (iii) calibrated the material parameters appearing in the theory using the stress-strain data from the compression experiments; (iv) numerically implemented the theory by writing a user-material subroutine for a widely used finite element program; and (v) conducted representative experiments to validate the predictive capability of our theory and its numerical implementation in complex three-dimensional geometries. By comparing the numerically predicted response in these validation simulations against measurements from corresponding experiments, we show that our theory is capable of reasonably accurately reproducing the experimental results. As a demonstration of the robustness of the three-dimensional numerical capability, we also show results from a simulation of the shape-recovery response of a stent made from the polymer when it is inserted in an artery modeled as a compliant elastomeric tube.     

FE-analysis of sheet metal forming processes using continuous contact treatment

Discrete meshes cause stepwise propagation of the contact nodes of a sheet despite the fact that the contact region in the actual forming process is altered very smoothly. This can cause problems of convergence and accuracy in contact-sensitive processes, such as a bending process. In this study, a scheme for a continuous contact treatment is proposed in order to consider the more realistic behavior of the contact phenomena during the forming process. For verification of the proposed method, the contact pressures and forming load are evaluated during the compression forming of a tube. The analysis of a hemispherical dome formed without a blank holder is also presented in order to investigate the effects of the proposed algorithm. The results show that the precise deformation mode is predicted by the utilization of the proposed method.     

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.     

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.     

On strain and damage interactions during tearing: 3D in situ measurements and simulations for a ductile alloy (AA2139-T3)

Strain and damage interactions during tearing of a ductile Al-alloy with high work hardening are assessed in situ and in 3D combining two recently developed experimental techniques, namely, synchrotron laminography and digital volume correlation. Digital volume correlation consists of registering 3D laminography images. Via simultaneous assessments of 3D strain and damage at a distance of 1-mm ahead of a notch root of a thin Compact Tension-like specimen, it is found that parallel crossing slant strained bands are active from the beginning of loading in a region where the crack will be slanted. These bands have an intermittent activity but are stable in space. Even at late stages of deformation strained bands can stop their activity highlighting the importance of plasticity on the failure process rather than damage softening. One void is followed over the loading history and seen to grow and orient along the slant strained band at very late stages of deformation. Void growth and strain are quantified. Gurson–Tvergaard–Needleman-type simulations using damage nucleation for shear, which is based on the Lode parameter, are performed and capture slant fracture but not the initial strain fields and in particular the experimentally found slant bands. The band formation and strain distribution inside and outside the bands are discussed further using plane strain simulations accounting for plastic material heterogeneity in soft zones.     

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.     

Coupled thermo-mechanical modelling of bulk-metallic glasses: Theory, finite-element simulations and experimental verification

A three-dimensional, finite-deformation-based constitutive model to describe the behavior of metallic glasses in the supercooled liquid region has been developed. By formulating the theory using the principles of thermodynamics and the concept of micro-force balance [Gurtin, M., 2000. On the plasticity of single crystals: free energy, microforces, plastic-strain gradients. J. Mech. Phys. Solids 48, 989-1036], a kinetic equation for the free volume concentration is derived by augmenting the Helmholtz free energy used for a conventional metallic alloy with a flow-defect free energy which depends on the free volume concentration and its spatial gradient. The developed constitutive model has also been implemented in the commercially available finite-element program ABAQUS/Explicit (2005) by writing a user-material subroutine. The constitutive parameters/functions in the model were calibrated by fitting the constitutive model to the experimental simple compression stress-strain curves conducted under a variety of strain-rates at a temperature in the supercooled liquid region [Lu, J., Ravichandran, G., Johnson, W., 2003. Deformation behavior of the Zr-Ti-Cu-Ni-Be bulk metallic glass over a wide range of strain-rates and temperatures. Acta Mater. 51, 3429-3443].With the model calibrated, the constitutive model was able to reproduce the simple compression stress-strain curves for jump-in-strain-rate experiments to good accuracy. Furthermore stress-strain responses for simple compression experiments conducted at different ambient temperatures within the supercooled liquid region were also accurately reproduced by the constitutive model. Finally, shear localization studies also show that the constitutive model can reasonably well predict the orientation of shear bands for compression experiments conducted at temperatures within the supercooled liquid region [Wang, G., Shen, J., Sun, J., Lu, Z., Stachurski, Z., Zhou, B., 2005. Compressive fracture characteristics of a Zr-based bulk metallic glass at high test temperatures. Mater. Sci. Eng. A 398, 82-87].     

Impact of woven fabric: Experiments and mesostructure-based continuum-level simulations

Woven fabric is an increasingly important component of many defense and commercial systems, including deployable structures, restraint systems, numerous forms of protective armor, and a variety of structural applications where it serves as the reinforcement phase of composite materials. With the prevalence of these systems and the desire to explore new applications, a comprehensive, computationally efficient model for the deformation of woven fabrics is needed. However, modeling woven fabrics is difficult due, in particular, to the need to simulate the response both at the scale of the entire fabric and at the meso-level, the scale of the yarns that compose the weave. Here, we present finite elements for the simulation of the three-dimensional, high-rate deformation of woven fabric. We employ a continuum-level modeling technique that, through the use of an appropriate unit cell, captures the evolution of the mesostructure of the fabric without explicitly modeling every yarn. Displacement degrees of freedom and degrees of freedom representing the change in crimp amplitude of each yarn family fully determine the deformed geometry of the mesostructure of the fabric, which in turn provides, through the constitutive relations, the internal nodal forces. In order to verify the accuracy of the elements, instrumented ballistic impact experiments with projectile velocities of 22-550 m/s were conducted on single layers of Kevlar^® fabric. Simulations of the experiments demonstrate that the finite elements are capable of efficiently simulating large, complex structures.     

The evolution of microstructure during twinning: Constitutive equations,finite-element simulations and experimental verification

In this work, we develop a rate-dependent, finite-deformation and crystal-mechanics-based constitutive theory which describes the twinning in single-crystal metallic materials. Central to the derivation of the constitutive equations are the use of fundamental thermodynamic laws and the principle of micro-force balance [Fried, E., Gurtin, M., 1994. Dynamic solid–solid transitions with phase characterized by an order parameter. Physica D 72, 287–308]. A robust numerical algorithm based on the constitutive model has also been written and implemented in the ABAQUS/Explicit [Abaqus reference manuals, 2007. SIMULIA, Providence, R.I.] finite-element program.     

On the use of a reduced enhanced solid-shell (RESS) element for sheet forming simulations

A recently proposed reduced enhanced solid-shell (RESS) element [Alves de Sousa, R.J., Cardoso, R.P.R., Fontes Valente, R.A., Yoon, J.W., Grácio, J.J., Natal Jorge, R.M., 2005. A new one-point quadrature enhanced assumed strain (EAS) solid-shell element with multiple integration points along thickness: Part I – Geometrically Linear Applications. International Journal for Numerical Methods in Engineering 62, 952–977; Alves de Sousa, R.J., Cardoso, R.P.R., Fontes Valente, R.A., Yoon, J.W., Grácio, J.J., Natal Jorge, R.M., 2006. A new one-point quadrature enhanced assumed strain (EAS) solid-shell element with multiple integration points along thickness: Part II – Nonlinear Applications. International Journal for Numerical Methods in Engineering, 67, 160–188.] is based on the enhanced assumed strain (EAS) method with a one-point quadrature numerical integration scheme. In this work, the RESS element is applied to large-deformation elasto-plastic thin-shell applications, including contact and plastic anisotropy. One of the main advantages of the RESS is its minimum number of enhancing parameters (only one), which when associated with an in-plane reduced integration scheme, circumvents efficiently well-known locking phenomena, leading to a computationally efficient performance when compared to conventional 3D solid elements. It is also worth noting that the element accounts for an arbitrary number of integration points through thickness direction within a single element layer. This capability has proven to be efficient, for instance, for accurately describing springback phenomenon in sheet forming simulations. A physical stabilization procedure is employed in order to correct the element’s rank deficiency. A general elasto-plastic model is also incorporated for the constitutive modelling of sheet forming operations with plastic anisotropy. Several examples including contact, anisotropic plasticity and springback effects are carried out and the results are compared with experimental data.     

Incorporation of yield surface distortion in finite deformation constitutive modeling of rigid-plastic hardening materials based on the Hencky logarithmic strain

In this paper a constitutive model for rigid-plastic hardening materials based on the Hencky logarithmic strain tensor and its corotational rates is introduced. The distortional hardening is incorporated in the model using a distortional yield function. The flow rule of this model relates the corotational rate of the logarithmic strain to the difference of the Cauchy stress and the back stress tensors employing deformation-induced anisotropy tensor. Based on the Armstrong–Fredrick evolution equation the kinematic hardening constitutive equation of the proposed model expresses the corotational rate of the back stress tensor in terms of the same corotational rate of the logarithmic strain. Using logarithmic, Green–Naghdi and Jaumann corotational rates in the proposed constitutive model, the Cauchy and back stress tensors as well as subsequent yield surfaces are determined for rigid-plastic kinematic, isotropic and distortional hardening materials in the simple shear deformation. The ability of the model to properly represent the sign and magnitude of the normal stress in the simple shear deformation as well as the flattening of yield surface at the loading point and its orientation towards the loading direction are investigated. It is shown that among the different cases of using corotational rates and plastic deformation parameters in the constitutive equations, the results of the model based on the logarithmic rate and accumulated logarithmic strain are in good agreement with anticipated response of the simple shear deformation.     

An energy-based anisotropic yield criterion for cellular solids and validation by biaxial FE simulations

The initial yield surface of 2D lattice materials is investigated under biaxial loading using finite element analyses as well as by analytical means. The sensitivity of initial yield surface to the dominant deformation mode is explored by using both low- and high-connectivity topologies whose dominant deformation mode is either local bending or strut stretching, respectively. The effect of microstructural irregularity on the initial yield surface is also examined for both topologies. A pressure-dependent anisotropic yield criterion, which is based on total elastic strain energy density, is proposed for 2D lattice structures, which can be easily extended for application to 3D cellular solids. Proposed criterion uses elastic constants and yield strengths under uniaxial loading, and does not rely on any arbitrary parameter. The analytical framework developed allows the introduction of new scalar measures of characteristic stresses and strains that are capable of representing the elastic response of anisotropic materials with a single elastic master line under multiaxial loading.