Experimental and numerical investigation of combined isotropic-kinematic hardening behavior of sheet metals

To prevent a sheet specimen from buckling subjected to a tension-compression cyclic loading, a new fixture has been developed to use with a regular tensile-compression machine. The novelty of this device lies in 4-block wedge design with pre-loaded springs. This design allows blocks to freely move in the vertical direction while providing the normal support to the entire length of the specimen during the tension-compression cycle. The entire test is easy to setup, which is another advantage of this design. In order to measure the strain accurately, the transmission type laser extensometer was utilized together with the implementation of double-side fins in the specimen. Experimental results of tension-compression tests are presented followed by a review of existing testing methods. In order to describe the accurate cyclic tension-compression behavior, the combined isotropic-kinematic hardening law based on the modified Chaboche model and the practical two-surface model based on Dafalias-Popov and Krieg models have been modified in this work, considering the permanent softening behavior during reverse loading and the non-symmetric behavior during reloading. Through tension-compression tests, the material characterization has been performed for three base materials, BH180, DP600 steels and AA6111-T4 sheets.     

Multiaxial cyclic deformation and non-proportional hardening employing discriminating load paths

Some novel discriminating multiaxial cyclic strain paths with incremental and random sequences were used to investigate cyclic deformation behavior of materials with low and high sensitivity to non-proportional loadings. Tubular specimens made of 1050 QT steel with no non-proportional hardening and 304L stainless steel with significant non-proportional hardening were used. 1050 QT steel was found to exhibit very similar behavior under various multiaxial loading paths, whereas significant effects of loading sequence were observed for 304L stainless steel. In-phase cycles with a random sequence of axial-torsion cycles on an equivalent strain circle were found to cause cyclic hardening levels similar to 90° out-of-phase loading of 304L stainless steel. In contrast, straining with a small increment of axial-torsion on an equivalent strain circle results in higher stress than for in-phase loading of 304L stainless steel, but the level of hardening is lower than for 90° out-of-phase loading. Tanaka’s non-proportionality parameter coupled with a Armstrong–Fredrick incremental plasticity model, and Kanazawa et al.’s empirical formulation as a representative of such empirical models were used to predict the stabilized stress response of the two materials under variable amplitude axial-torsion strain paths. Consistent results between experimental observations and predictions were obtained by employing the Tanaka’s non-proportionality parameter. In contrast, the empirical model resulted in significant over-prediction of stresses for 304L stainless steel.     

On the modeling of hardening in metals during non-proportional loading

The purpose of the current work is the formulation and initial application of a phenomenological model for hardening effects in metals subject to non-proportional loading histories characterized by one or more loading-path changes. This model is closely related to the incremental model of Teodosiu and Hu [Teodosiu, C., Hu, Z., 1995. Evolution of the intragranular microstructure at moderate and large strains: modelling and computational significance. In: Shen, S.F., Dawson, P.R. (Eds.), Simulation of Materials Processing: Theory, Methods and Applications. Balkema, Rotterdam, pp. 173–182; Teodosiu, C., Hu, Z., 1998. Microstructure in the continuum modelling of plastic anisotropy. In: Proceedings of 19th Risø International Symposium on Material’s Science: Modelling of Structure and Mechanics of Materials from Microscale to Product. Risø National Laboratory, Roskilde, Denmark, pp. 149–168]. Like their model, the current model captures in particular hardening stagnation after a load reversal as well as cross-hardening after orthogonal loading-path changes. On the other hand, the two models predict qualitatively different behavior during loading-path changes which take place purely in the inelastic range. Such is the case for example during orthogonal loading-path changes from uniaxial tension to simple shear without release, or during monotonic simple shear, or during deep-drawing. As shown by the experimental results reported on in the current work for the mild steel DC06, significant cross-hardening can occur during continuous orthogonal loading-path changes. Beyond this, the current model accounts in an approximate way for the possible effects of texture development on the material behavior with the help of the plastic spin. After investigating the behavior of the current model for various ideal two-stage loading histories (e.g., tension-shear), the current work ends with a comparison of standard combined hardening and current approaches in the context of the simulation of internal stress development and residual stresses during deep-drawing and the resultant springback after ring-splitting.     

A novel approach for anisotropic hardening modeling. Part II: Anisotropic hardening in proportional and non-proportional loadings,application to initially isotropic material

The modeling of anisotropic hardening, in particular for non-proportional loading paths, is a challenging task for advanced macroscopic models. The complex distortion of the yield locus is related to the activation and cross-hardening of different slip systems, depending on crystallographic orientations. These physical mechanisms can be taken into account in polycrystalline models but the computation times are enormous. The novel approach detailed in Part I (Rousselier et al., 2009) consists in: (i) drastically reducing the number of crystallographic orientations to save the computation cost, (ii) applying a parameter calibration procedure to obtain a good agreement with the experimental database. This methodology is first applied here to the anisotropic hardening in the proportional loadings of the strongly anisotropic aluminum alloy of Part I. Very good modeling is achieved with only eight crystallographic orientations. Different levels of additional hardening in biaxial proportional loading as compared to uniaxial loading can be modeled with the same polycrystalline model. For this, only the parameter calibration has to be performed with different databases. The same methodology has also been applied for the modeling of isotropic behavior. The best compromise between model accuracy and numerical cost is obtained with fourteen orientations. The deviations from isotropy are acceptable in all loading directions. Different levels of hardening in orthogonal loading: simple shear followed by simple tension, are achieved without any modification of the model equations. Only the parameter calibration has to be performed with different hardening levels in the database. FE calculations of a deep drawing test have been performed. The CPU time of the polycrystalline model is only five times larger than that with the simple von Mises model. The CPU time with texture evolution is further increased by a factor of two. The effects of texture evolution in rolling of the initially isotropic fcc material have been investigated. The resulting texture and hardening are qualitatively good.     

On the evolution of isotropic and kinematic hardening with finite plastic deformation Part I: compression/tension loading of OFHC copper cylinders

Compression tests followed by tension tests after re-machining were performed on annealed oxygen-free-high-conductivity copper cylinders. These tests were conducted at nine levels of maximum strain ranging from 5 to 50%. From this data, isotropic and kinematic hardening were calculated using 50, 1000 and 2000 microstrain offset definitions. Both isotropic and kinematic hardening were found to depend on the yield definition. Isotropic hardening, which increased with plastic strain with no signs of saturation, also increased with larger offset definition of yield. Kinematic hardening, which increased to 40% strain and appeared to saturate thereafter, decreased with higher offset definitions of yield.     

Behaviors of three BCC metals during non-proportional multi-axial loadings: experiments and modeling

Non-proportional torsion–tension and biaxial-compressive experimental results are presented on tantalum, tantalum alloy with 2.5% tungsten, and AerMet 100 steel. These test results form a comprehensive set of data to show the material behaviors at finite strain and wide strain-rate range. Using the parameter set determined from uniaxial constant strain-rate compressive and tensile tests, the capability of a new constitutive model (Khan, A.S., Liang, R., 1999. Behaviors of three BCC metal over a wide range of strain rates and temperatures: experiments and modeling. International Journal of Plasticity 15, 1089–1109) is shown to accurately predict complex loading paths of current experimental results. Using von Mises equivalent strain, stress, and strain rate, the constitutive model gives excellent predictions of these non-proportional experimental results.     

An elasto-plastic damage constitutive theory and its prediction of evolution of subsequent yield surfaces and elastic constants

Based on pair functional potentials, Cauchy-Born rule and slip mechanism, a material model assembling with spring-bundle components, a cubage component and slip components is established to describe the elasto-plastic damage constitutive relation under finite deformation. The expansion/shrink, translation and distortion of yield surfaces can be calculated based on the hardening rule and Bauschinger effect defined on the slip component level. Both kinematic and isotropic hardening are included. Numerical simulations and predictions under tension, torsion, and combined tension-torsion proportional/non-proportional loading are performed to obtain the evolution of subsequent yield surfaces and elastic constants and compare with two sets of experimental data in literature, one for a very low work hardening aluminum alloy Al 6061-T6511, and another for a very high work hardening aluminum alloy annealed 1100 Al. The feature of the yield surface in shape change, which presents a sharp front accompanied by a blunt rear under proportional loading, is described by the latent hardening and Bauschinger effect of slip components. Further, the evolution law of subsequent yield surfaces under different proportional loading paths is investigated in terms of their equivalence. The numerical simulations under non-proportional loading conditions for annealed 1100 Al are performed, and the subsequent yield surfaces exhibit mixed cross effect because the kinematic hardening and isotropic hardening follow different evolution tendency when loading path changes. The results of non-proportional loading demonstrate that the present model has the ability to address the issue of complex loading due to the introduction of state variables on slip components. Moreover, as an elasto-plastic damage constitutive model, the present model can also reflect the variation of elastic constants through damage defined on the spring-bundle components.     

A large deformation theory for rate-dependent elastic–plastic materials with combined isotropic and kinematic hardening

We have developed a large deformation viscoplasticity theory with combined isotropic and kinematic hardening based on the dual decompositions F=F^eF^p [Kröner, E., 1960. Allgemeine kontinuumstheorie der versetzungen und eigenspannungen. Archive for Rational Mechanics and Analysis 4, 273–334] and [Lion, A., 2000. Constitutive modelling in finite thermoviscoplasticity: a physical approach based on nonlinear rheological models. International Journal of Plasticity 16, 469–494]. The elastic distortion F^e contributes to a standard elastic free-energy ψ^(e), while , the energetic part of F^p, contributes to a defect energy ψ^(p) – these two additive contributions to the total free energy in turn lead to the standard Cauchy stress and a back-stress. Since F^e=FF^p-1 and , the evolution of the Cauchy stress and the back-stress in a deformation-driven problem is governed by evolution equations for F^p and – the two flow rules of the theory.We have also developed a simple, stable, semi-implicit time-integration procedure for the constitutive theory for implementation in displacement-based finite element programs. The procedure that we develop is “simple” in the sense that it only involves the solution of one non-linear equation, rather than a system of non-linear equations. We show that our time-integration procedure is stable for relatively large time steps, is first-order accurate, and is objective.     

Orthotropic strain rate potential for the description of anisotropy in tension and compression of metals

In this paper, a macroscopic anisotropic strain rate potential, which can describe both the anisotropy and tension-compression asymmetry of the plastic response of textured metals is derived. This strain rate potential is the exact work-conjugate of the anisotropic stress potential CPB06 of Cazacu et al. (2006). Application of the developed strain rate potential to HCP high-purity alpha-titanium is 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.     

An atomistically-informed dislocation dynamics model for the plastic anisotropy and tension-compression asymmetry of BCC metals

Atomistic simulations have shown that a screw dislocation in body-centered cubic (BCC) metals has a complex non-planar atomic core structure. The configuration of this core controls their motion and is affected not only by the usual resolved shear stress on the dislocation, but also by non-driving stress components. Consequences of the latter are referred to as non-Schmid effects. These atomic and micro-scale effects are the reason slip characteristics in deforming single and polycrystalline BCC metals are extremely sensitive to the direction and sense of the applied load. In this paper, we develop a three-dimensional discrete dislocation dynamics (DD) simulation model to understand the relationship between individual dislocation glide behavior and macro-scale plastic slip behavior in single crystal BCC Ta. For the first time, it is shown that non-Schmid effects on screw dislocations of both {110} and {112} slip systems must be implemented into the DD models in order to predict the strong plastic anisotropy and tension-compression asymmetry experimentally observed in the stress-strain curves of single crystal Ta. Incorporation of fundamental atomistic information is critical for developing a physics-based, predictive meso-scale DD simulation tool that can connect length/time scales and investigate the underlying mechanisms governing the deformation of BCC metals.     

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.     

A comparison of boundary element and finite element methods for modeling axisymmetric polymeric drop deformation

A modified boundary element method (BEM) and the DEVSS‐G finite element method (FEM) are applied to model the deformation of a polymeric drop suspended in another fluid subjected to start‐up uniaxial extensional flow. The effects of viscoelasticity, via the Oldroyd‐B differential model, are considered for the drop phase using both FEM and BEM and for both the drop and matrix phases using FEM. Where possible, results are compared with the linear deformation theory. Consistent predictions are obtained among the BEM, FEM, and linear theory for purely Newtonian systems and between FEM and linear theory for fully viscoelastic systems. FEM and BEM predictions for viscoelastic drops in a Newtonian matrix agree very well at short times but differ at longer times, with worst agreement occurring as critical flow strength is approached. This suggests that the dominant computational advantages held by the BEM over the FEM for this and similar problems may diminish or even disappear when the issue of accuracy is appropriately considered. Fully viscoelastic problems, which are only feasible using the FEM formulation, shed new insight on the role of viscoelasticity of the matrix fluid in drop deformation. Copyright © 2001 John Wiley & Sons, Ltd.     

Effect of microstructure on the hardening and softening behavIors of polycrystalline shape memory alloys Part I: Micromechanics constitutive modeling

The effects of microstructure and its evolution on the macroscopic superelastic stress-strain response of polycrystalline Shape Memory Alloy (SMA) are studied by a microstructure-based constitutive model developed in this paper. The model is established on the following basis: (1) the transformation conditions of the unconstrained single crystal SMA microdomain (to be distinguished from the bulk single crystal), which serve as the local criterion for the derivation of overall transformation yield conditions of the polycrystal; (2) the micro- to macro-transition scheme by which the connection between the polycrystal aggregates and the single crystal microdomain is established and the macroscopic transformation conditions of the polycrystal SMA are derived; (3) the quantitative incorporation of three microstructure factors (i.e., nucleation, growth and orientation distribution of martensite) into the modeling. These microstructural factors are intrinsic of specific polycrystal SMA systems and the role of each factor in the macroscopic constitutive response is quantitatively modeled. It is demonstrated that the interplay of these factors will result in different macroscopic transformation kinematics and kinetics which are responsible for the observed macroscopic stress-strain hardening or softening response, the latter will lead to the localization and propagation of transformation bands in TiNi SMA. The project supported by the Research Grant Committee (RGC) of Hong Kong SAR, the National Natural Science Foundation of China and the Provincial Natural Foundation of Jiangxi Province of China     

Asymptotic integration of elasto-plastic constitutive relation for isotropic hardening material and its computational precision

In this paper, an asymptotic expansion solution of the constitutive equation of hardening materials is presented. Its 1st asymptotic integration can give an approximate one with good enough accuracy and the second asymptotic one improves the precision of solutions further. The steps of its algorithms are fairly simple and clear, and its computational workload is considerably reduced. It can be easily incorporated into a general purpose finite element program.The Chinese original of this article was published in the Chinese edition ofActa Mechanica Solida Sinica, No. 1, 1986.     

Experimental study and electromechanical model analysis of the nonlinear deformation behavior of IPMC actuators

This paper develops analytical electromechanical formulas to predict the mechanical deformation of ionic polymer–metal composite(IPMC) cantilever actuators under DC excitation voltages. In this research, IPMC samples with Pt and Ag electrodes were manufactured, and the large nonlinear deformation and the effect of curvature on surface electrode resistance of the IPMC samples were investigated experimentally and theoretically. A distributed electrical model was modified for calculating the distribution of voltage along the bending actuator. Then an irreversible thermodynamic model that could predict the curvature of a unit part of an IPMC actuator is combined with the electrical model so that an analytical electromechanical model is developed. The electromechanical model is then validated against the experimental results obtained from Pt-and Ag-IPMC actuators under various excitation voltages. The good agreement between the electromechanical model and the actuators shows that the analytical electromechanical model can accurately describe the large nonlinear quasi-static deflection behavior of IPMC actuators.     

循环硬化材料本构模型的隐式应力积分和有限元实现

针对新发展的、能够描述循环硬化行为应变幅值依赖性的粘塑性本构模型,讨论了它的数值实现方法。首先,为了能够对材料的循环棘轮行为（Ratcheting）和循环应力松弛现象进行描述,对已有的本构模型进行了改进;然后,在改进模型的基础上,建立了一个新的、全隐式应力积分算法,进而推导了相应的一致切线刚度（Consistent Tangent Modulus）矩阵的表达式;最后,通过ABAQUS用户材料子程序UMAT将上述本构模型进行了有限元实现,并通过一些算例对一些构件的循环变形行为进行了有限元数值模拟,讨论了该类本构模型有限元实现的必要性和合理性。     

Analytical removal of singularity and direct evaluation of singular integrals in 3-D elastoplastic finite deformation analysis by BEM

As a further development of the present authors' research work [1,2], in this paper a method of the so-called quadratic pentahedron polar co-ordinate transformation and analytical removal of singularity of Cauchy principal value singular integrals is proposed to evaluate the strongly singular integrals in the sense of Cauchy principal values and the weakly singular integrals over quadratic internal cells in 3-D elastoplastic finite deformation analysis by BEM. First, a quadratic pentahedron polar co-ordinate transformation technique is used to reduce the order of singularity of the singular integrals. Then, a form of Gauss' theorem is introduced to remove the singularity in the Cauchy principal value singular integrals analytically. Therefore, the evaluation of all those strongly and weakly singular integrals can be carried out by standard Gaussian quadrature accurately and efficiently. Numerical examples of the 3-D elastoplastic problem and 3-D finite deformation problem are given to demonstrate that the method possesses good accuracy and numerical stability, and is convenient to implement. The method in this paper can be applied extensively to evaluating the singular integrals over cubic and higher order elements.