Characterization of strain rate effects in sheet laser forming

This work presents numerical simulations and experimental validation of sheet laser forming processes using a single-step straight path with different laser beam powers (four levels ranging from 30 W to 120 W) and scanning speeds (four levels ranging from 5 mm/s to 20 mm/s) in graphite-coated AISI 304 stainless steel 0.6-mm-thick sheets. The numerical simulations of these cases are performed via a coupled thermomechanical finite element formulation accounting for large strains, temperature-dependent material properties and convection–radiation phenomena. Firstly, a rate-independent plastic model is used. Although this model adequately predicts the final bending angle for the cases achieving relatively low maximum temperatures, i.e. cases with low laser beam powers and high scanning speeds, it fails in describing the deformation pattern for the cases with higher maximum temperatures, i.e. cases with high laser beam powers and low scanning speeds. Secondly, in order to overcome this drawback, a rate-dependent viscoplastic model including a stress-dependent viscosity law is proposed to simulate the same cases. The final bending angles provided by this model are found to be in good agreement with the experimental measurements for the whole ranges of laser beam power and scanning speed studied in this work. Therefore, the use of this viscoplastic model in the simulation of sheet laser forming allows us to conclude that the strain rate effects, which mainly play a relevant role at high temperatures, can be adequately characterized.     

Micromechanics, macromechanics and constitutive modeling of the elasto-viscoplastic deformation of rubber-toughened glassy polymers

Under certain conditions, such as sufficiently low temperatures, high loading rates and/or highly triaxial stress states, glassy polymers display an unfavorable characteristic—brittleness. A technique used for reducing the brittleness (increasing the fracture toughness) of these materials is rubber toughening. While there is significant qualitative understanding of the mechanical behavior of rubber-toughened polymers, quantitative modeling tools for the large-strain deformation of rubber-toughened glassy polymers are largely lacking.In this paper, we develop a suite of numerical tools to investigate the mechanical behavior of rubber-toughened glassy polymers, with emphasis on rubber-toughened polycarbonate. The rubber particles are modeled as voids in view of their deformation-induced cavitation early during deformation. A three-dimensional micromechanical model of the heterogeneous microstructure is developed to study the effects of initial rubber particle (void) volume fraction on the underlying elasto-viscoplastic deformation mechanisms in the material, and how these mechanisms influence the macroscopic response of the material. A continuum-level constitutive model is developed for the large-strain elasto-viscoplastic deformation of porous glassy polymers, and it is calibrated against micromechanical modeling results for porous polycarbonate. The constitutive model can be used to study various boundary value problems involving rubber-toughened (porous) glassy polymers. As an example, the case of an axisymmetric notched bar is simulated for the case of polycarbonate with varying levels of initial porosity. The quality of the constitutive model calibration is assessed using a multi-scale modeling approach.     

Model identification and FE simulations: Effect of different yield loci and hardening laws in sheet forming

The bi-axial experimental equipment [Flores, P., Rondia, E., Habraken, A.M., 2005a. Development of an experimental equipment for the identification of constitutive laws (Special Issue). International Journal of Forming Processes] developed by Flores enables to perform Bauschinger shear tests and successive or simultaneous simple shear tests and plane strain tests. Flores investigates the material behavior with the help of classical tensile tests and the ones performed in his bi-axial machine in order to identify the yield locus and the hardening model. With tests performed on one steel grade, the methods applied to identify classical yield surfaces such as [Hill, R., 1948. A theory of the yielding and plastic flow of anisotropic materials. Proceedings of the Royal Society of London A 193, 281–297; Hosford, W.F., 1979. On yield loci of anisotropic cubic metals. In: Proceedings of the 7th North American Metalworking Conf. (NMRC), SME, Dearborn, MI, pp. 191–197] ones as well as isotropic Swift type hardening, kinematic Armstrong–Frederick or Teodosiu and Hu hardening models are explained. Comparison with the Taylor–Bishop–Hill yield locus is also provided. The effect of both yield locus and hardening model choices is presented for two applications: plane strain tensile test and Single Point Incremental Forming (SPIF).     

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.     

Thermo-mechanical large deformation response and constitutive modeling of viscoelastic polymers over a wide range of strain rates and temperatures

A phenomenological one-dimensional constitutive model, characterizing the complex and highly nonlinear finite thermo-mechanical behavior of viscoelastic polymers, is developed in this investigation. This simple differential form model is based on a combination of linear and nonlinear springs with dashpots, incorporating typical polymeric behavior such as shear thinning, thermal softening at higher temperatures and nonlinear dependence on deformation and loading rate. Another model, of integral form, namely the modified superposition principle (MSP), is also modified further and used to show the advantage of the newly developed model over MSP. The material parameters for both models are determined for Adiprene-L100, a polyurethane based rubber. The constants once determined are then utilized to predict the behavior under strain rate jump compression, multiple step stress relaxation loading experiment and free end torsion experiments. The new constitutive model shows very good agreement with the experimental data for Adiprene-L100 for the various finite loading paths considered here and provides a flexible framework for a three-dimensional generalization.     

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.     

Path-dependence of the forming limit stresses in a sheet metal

The effect of changing strain paths on the forming limit stresses of sheet metals is investigated using the Marciniak–Kuczyński model and a phenomenological plasticity model with non-normality effects [Kuroda, M., Tvergaard, V., 2001. A phenomenological plasticity model with non-normality effects representing observations in crystal plasticity. J. Mech. Phys. Solids 49, 1239–1263]. Forming limits are simulated for linear stress paths and two types of combined loading: a combined loading consisting of two linear stress paths in which unloading is included between the first and second loadings (combined loading A), and combined loading in which the strain path is abruptly changed without unloading (combined loading B). The forming limit stresses calculated for combined loading A agree well with those calculated for the linear stress paths, while the forming limit curves in strain space depend strongly on the strain paths. The forming limit stresses calculated for the combined loading B do not, however, coincide with those calculated for the linear stress paths. The strain-path dependence of the forming limit stress is discussed in detail by observing the strain localization process.     

Mathematical modelling of the superplastic forming of a long rectangular sheet

Standard finite element software (ANSYS FEM-code) is used to model the superplastic forming of a superplastic sheet alloy into a rectangular die. As distinct from other known approaches the boundary value problem is stated in the present study in terms of the theory of creep. The results of finite element calculations are found to be in a good agreement with corresponding analytical solutions and experimental data for titanium sheet alloy Ti-6Al-4V. The approach suggested can be used in practice for estimating the current geometry and the thickness of the dome and the time intervals both for the constant gas pressure forming and for the constant strain rate forming.     

A constitutive model for the formation of martensite in austenitic steels under large strain plasticity

A constitutive model for diffusionless phase transitions in elastoplastic materials undergoing large deformations is developed. The model takes basic thermodynamic relations as its starting point and the phase transition is treated through an internal variable (the phase fractions) approach. The usual yield potential is used together with a transformation potential to describe the evolution of the new phase. A numerical implementation of the model is presented, along with the derivation of a consistent algorithmic tangent modulus. Simulations based on the presented model are shown to agree well with experimental findings. The proposed model provides a robust tool suitable for large-scale simulations of phase transformations in austenitic steels undergoing extensive deformations, as is demonstrated in simulations of the necking of a bar under tensile loading and also in simulations of a cup deep-drawing process.     

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.     

Experimental investigation and modeling of the tension behavior of polycarbonate with temperature effects from low to high strain rates

The effects of strain rate and temperature on the tension stress–strain responses of polycarbonate are experimentally investigated over a wide range of strain rates (0.001–1700 s⁻¹) and temperatures (0–120 °C). A modified split Hopkinson tension bar is used for high-rate uniaxial tension tests. Experimental results indicate that the stress–strain responses of polycarbonate at high strain rates exhibit the nonlinear characteristics including the obvious yielding and strain softening. The tension behavior is strongly dependent on the strain rate and temperature. The values of yield stress and strain at yield present a dramatic increase at higher strain rates and decrease with the increase in temperature. Moreover, there exists a significant rate-sensitivity transition in the polycarbonate tension yield behavior. Based on the experimental investigation, a physically based three-dimensional elastoplastic constitutive model for the finite deformation of glassy polymers is used to characterize the rate-temperature dependent yield and post-yield behavior of polycarbonate when subjected to tension loading. The model results are shown close to the experimental data within the investigated strain-rate and temperature ranges.     

Investigation of advanced strain-path dependent material models for sheet metal forming simulations

Sheet metal forming processes often involve complex loading sequences. To improve the prediction of some undesirable phenomena, such as springback, physical behavior models should be considered. This paper investigates springback behavior predicted by advanced elastoplastic hardening models which combine isotropic and kinematic hardening and take strain-path changes into account. A dislocation-based microstructural hardening model formulated from physical observations and the more classical cyclic model of Chaboche have been considered in this work. Numerical implementation was carried out in the ABAQUS software using a return mapping algorithm with a combined backward Euler and semi-analytical integration scheme of the constitutive equations. The capability of each model to reproduce transient hardening phenomena at abrupt strain-path changes has been shown via simulations of sequential rheological tests. A springback analysis of strip drawing tests was performed in order to emphasize the impact of several influential parameters, namely: process, numerical and behavior parameters. The effect of the two hardening models with respect to the process parameters has been specifically highlighted.     

Behaviors of three BCC metal over a wide range of strain rates and temperatures: experiments and modeling

The elastic–plastic behaviors of three body-centered cubic metals, tantalum, tantalum alloy with 2.5% tungsten, and AerMet 100 steel, are presented over a wide range of strains (15%), strain rates (10⁻⁶–10⁴ s⁻¹) and temperatures (77–600°F). Johnson-Cook and Zerilli-Armstrong models were found inadequate to describe the observations. A new viscoplastic model is proposed based on these experimental results. The proposed constitutive model gives good correlations with these experimental results and strain-rate jump experiments. In the next paper (Liang, R., Khan A.S., 2000. Behaviors of three BCC metals during non-proportional multi-axial loadings and predictions using a recently proposed model. International Journal of Plasticity, in press), multi-axial loading results on these materials and comparison with the proposed model will be presented.     

Multiaxial and non-proportional loading responses,anisotropy and modeling of Ti–6Al–4V titanium alloy over wide ranges of strain rates and temperatures

Results from a series of multiaxial loading experiments on the Ti–6Al–4V titanium alloy are presented. Different loading conditions are applied in order to get the comprehensive response of the alloy. The strain rates are varied from the quasi-static to dynamic regimes and the corresponding material responses are obtained. The specimen is deformed to large strains in order to study the material behavior under finite deformation at various strain rates. Torsional Kolsky bar is used to achieve shear strain rates up to 1000 s⁻¹. Experiments are performed under non-proportional loading conditions as well as dynamic torsion followed by dynamic compression at various temperatures. The non-proportional loading experiments comprise of an initial uniaxial loading to a certain level of strain followed by biaxial loading, using a channel-type die at various rates of loadings. All the non-proportional experiments are carried out at room temperature. Experiments are also performed to investigate the anisotropic behavior of the alloy. An orthotropic yield criterion [proposed by Cazacu, O., Plunkett, B., Barlat, F., 2005. Orthotropic yield criterion for hexagonal closed packed metals. International Journal of Plasticity 22, 1171–1194.] for anisotropic hexagonal closed packed materials with strength differential is used to generate the yield surface. Based on the definition of the effective stress of this yield criterion, the observed material response for the different loading conditions under large deformation is modeled using the Khan–Huang–Liang (KHL) equation assuming isotropic hardening. The model constants used in the present study, were pre-determined from the extensive uniaxial experiments presented in the earlier paper [Khan, A.S., Suh, Y.S., Kazmi R., 2004. Quasi-static and dynamic loading responses and constitutive modeling of titanium alloys. International Journal of Plasticity 20, 2233–2248]. The model predictions are found to be extremely close to the observed material response.     

Mechanically alloyed nanocrystalline iron and copper mixture: behavior and constitutive modeling over a wide range of strain rates

A comprehensive study on the response of a nanocrystalline iron and copper mixture (80% Fe and 20% Cu) to quasi-static and dynamic loading is performed. The constitutive model developed earlier by Khan, Huang & Liang (KHL) is extended to include the responses of nanocrystalline metallic materials. The strain rate and grain size dependent behaviors of porous nanocrystalline iron-copper mixture were determined experimentally for both static and dynamic loading. A viscoplastic model is formulated by associating the modified KHL model (representing the fully dense matrix behavior), and Gurson's plastic potential which provides the yield criteria for porous material. Simulations of uniaxial compressive deformations of iron-copper mixture with different initial porosity, grain size and at a wide range of strain rate (10⁻⁴ to 10³ s⁻¹) are made. The numerical results correlate well with the experimental observations.     

Multiscale modeling of the plasticity in an aluminum single crystal

This paper describes a numerical, hierarchical multiscale modeling methodology involving two distinct bridges over three different length scales that predicts the work hardening of face centered cubic crystals in the absence of physical experiments. This methodology builds a clear bridging approach connecting nano-, micro- and meso-scales. In this methodology, molecular dynamics simulations (nanoscale) are performed to generate mobilities for dislocations. A discrete dislocations numerical tool (microscale) then uses the mobility data obtained from the molecular dynamics simulations to determine the work hardening. The second bridge occurs as the material parameters in a slip system hardening law employed in crystal plasticity models (mesoscale) are determined by the dislocation dynamics simulation results. The material parameters are computed using a correlation procedure based on both the functional form of the hardening law and the internal elastic stress/plastic shear strain fields computed from discrete dislocations. This multiscale bridging methodology was validated by using a crystal plasticity model to predict the mechanical response of an aluminum single crystal deformed under uniaxial compressive loading along the [4 2 1] direction. The computed strain-stress response agrees well with the experimental data.     

Modeling and validation of the large deformation inelastic response of amorphous polymers over a wide range of temperatures and strain rates

A robust physically consistent three-dimensional constitutive model is developed to describe the finite mechanical response of amorphous polymers over a wide range of temperatures and strain rates, including the rubbery region and for impact loading rates. This thermomechanical model is based on an elastic–viscoplastic rheological approach, wherein the effects of temperature, strain rate, and hydrostatic pressure are accounted for. Intramolecular, as well as intermolecular, interactions under large elastic–inelastic behavior are considered for the mechanisms of deformation and hardening. For a wide range of temperatures and strain rates, our simulated results for poly(methyl methacrylate) (PMMA) and polycarbonate (PC) are in good agreement with experimental observations.     

A constitutive model of aluminum foam for crash simulations

A new constitutive model for metallic foams is developed to overcome the deficiency of existing models in commercial finite element codes. The proposed constitutive model accounts for volume changes under hydrostatic compression and combines the hydrostatic pressure and von Mises stress into one yield function. The change of the compressibility of the metallic foam is handled in the constitutive model by allowing for shape changes of the yield surface in the hydrostatic pressure-von Mises stress space. The backward Euler method is adopted to integrate the constitutive equations to achieve numerical accuracy and stability. The model is implemented into LS-DYNA as a user-defined subroutine, verified with existing solutions, and validated with foam testing data. The verified and validated model is then utilized in the crushing simulations of foam-filled columns with square and hexagonal cross-sections. Two constitutive models are studied: the first using an exponential function to describe the relationship of plastic Poisson's ratio with respect to true strain and the second using linear interpolation function as an alternative approximation. The new foam model provides satisfactory prediction of crushing forces and deformed shapes compared to experimental results. Additionally, the new foam model was shown to have better numerical stability and accuracy than existing LS-DYNA built-in material models.     

Multi-level modeling of woven glass/epoxy composite for multilayer printed circuit board applications

The objective of this paper is to develop a hybrid homogenization method to predict the elastic properties of a common woven glass/epoxy composite substrate for multilayer circuit board applications. Comprehensive high resolution 3D finite element (FE) models of a quarter of the repeated unit cell (RUC) for the woven glass/epoxy composite were developed based on different micromechanical schemes. . Specifically, four different micromechanics schemes were investigated: self-consistent, Mori–Tanaka, three-phase approach and composite cylinder assemblage (CCA). The element based strain concentration matrices were determined and used to obtain the homogenized woven glass/epoxy composite properties via a specially developed MATLAB code. Attention was further devoted to the predictions of the homogenized elastic moduli of the multilayer printed circuit board (PCB). The results from our simulations, based on Mori–Tanaka and CCA, are in good agreement with existing experimental results, indicating that the newly proposed homogenization scheme can be used as a design tool to predict the overall properties of woven composite materials typically used in multilayer PCB applications.