Modeling mixed hardening of alkali halides with a modified version of an internal state variables model

In the ductile regime of inelastic flow, alkali halides behave similarly to various metals. Their behavior is fully plastic, shows strong rate and history dependencies, and exhibits anisotropic hardening. In this paper, the authors present the most recent version of SUVIC, a viscoplastic model with internal state variables (ISV). After recalling the main features of alkali halides inelastic response, the modified equations of the model are presented. They are then used to represent the behavior of polycrystalline sodium chloride submitted to conventional triaxial compression (CTC) and reduced triaxial extension (RTE) tests, using results that highlight mixed (kinematic and isotropic) hardening of the material, hence showing a type of Bauschinger effect. A discussion follows.     

Representation of cyclic hardening and softening properties using continuous variables

Our aim is the modeling of cyclic hardening, cyclic softening, cyclic mean stress relaxation, and additional nonproportional cyclic hardening. We do so by means of hardening functionals for back stress and yield stress without employing additional memory surfaces. Rather, we suppose all quantities to evolve simultaneously during elastic-plastic loading in a continuous manner. The basic idea is to formulate evolution equations for the hardening variables, which are of the “hardening/dynamic recovery” format with respect to a transformed arc length. The corresponding transformation is influenced by continuously evolving parameters, measuring strain amplitude and nonproportionality during the recent process history. Although the resulting, model has a very simple structure, it is capable of describing the basic phenomena under quite general loading conditions.     

Finite and symmetric thermomechanical waves in materials with internal state variables

The paper is devoted to the study of wave propagation through materials with internal state variables under two main assumptions: finite speed of propagation and symmetry of acceleration waves. Additional constitutive assumptions are also used: heat flux does not depend on the temperature gradient and the time-derivatives of internal state variables are linear functions of the temperature gradient. Under all these hypotheses, in the neighborhood of a strong equilibrium state, one finds four real and symmetric possible acceleration waves, at least two of them being coupled waves, and heat flux results an internal state variable. All these results are obtained in the general three-dimensional case. As an illustration, the isotropic linear theory is considered, where both acceleration and shock waves are treated.     

Constitutive modeling for anisotropic/asymmetric hardening behavior of magnesium alloy sheets

Magnesium alloy sheets have been extending their field of applications to automotive and electronic industries taking advantage of their excellent light weight property. In addition to well-known lower formability, magnesium alloys have unique mechanical properties which have not been thoroughly studied: high in-plane anisotropy/asymmetry of yield stress and hardening response. The reason of the unusual mechanical behavior of magnesium alloys has been understood by the limited symmetry crystal structure of HCP metals and thus by deformation twinning. In this paper, the phenomenological continuum plasticity models considering the unusual plastic behavior of magnesium alloy sheet were developed for a finite element analysis. A hardening law based on two-surface model was further extended to consider the general stress–strain response of metal sheets including Bauschinger effect, transient behavior and the unusual asymmetry. Three deformation modes observed during the continuous in-plane tension/compression tests were mathematically formulated with simplified relations between the state of deformation and their histories. In terms of the anisotropy and asymmetry of the initial yield stress, the Drucker–Prager’s pressure dependent yield surface was modified to include the anisotropy of magnesium alloy. The numerical formulations and characterization procedures were also presented and finally the correlation of simulation with measurements was performed to validate the proposed theory.     

A critical evaluation and extension of internal state variables consititutive models

This study provides a critical comparison of the internal state variables constitutive model sproposed recently in the literature for metals at elevated temperatures and deformed over a wide range of strain rates. As a result of the observations from this comparative evaluation, a set of modified constitutive relations is also presented. The experimental data available in the literature for 1100 aluminum and a fully dense high purity aluminumm were used along with a nonlinear least squares fitting procedure to estimate the material constants. A nonlinear functional dependence of the effective stress with respect to the internal state variable in the evolution relationship was determined to model more effectively, compared to previous models, the case that a material may exhibit both hardening and softening behaviors, respectively, corresponding to different applied strain rates at the same temperature. The modified relations presented here propose a simple form of this nonlinear function, which accurately models the material response and which allows a more robust implementation in numerical simulations. The deformation of an aluminium gradient specimen was studied for further comparison of the three constitutive models considered and for a demonstration of the validity of this model.     

Behaviour of the acceleration and shock waves in materials with internal state variables

The explicit expressions for the change in the amplitudes of one-dimensional acceleration and shock waves propagating through arbitrary homogeneous materials described by the strain and internal state variables/parameters/are derived. The existence of a critical amplitude β for the acceleration wave and a critical strain gradient λ for the shock wave is established. For an infinitesimal shock wave the general form of the solution of the governing differential equation is furnished. The differential equations for the amplitudes of these two kind of waves are applied to an elastic-viscoplastic material.     

Characterization of the post-necking strain hardening behavior using the virtual fields method

The present study aims at characterizing the post-necking strain hardening behavior of three sheet metals having different hardening behavior. Standard tensile tests were performed on sheet metal specimens up to fracture and heterogeneous logarithmic strain fields were obtained from a digital image correlation technique. Then, an appropriate elasto-plastic constitutive model was chosen. Von Mises yield criterion under plane stress and isotropic hardening law were considered to retrieve the relationship between stress and strain. The virtual fields method (VFM) was adopted as an inverse method to determine the constitutive parameters by calculating the stress fields from the heterogeneous strain fields. The results show that the choice of a hardening law which can describe the hardening behavior accurately is important to derive the true stress–strain curve. Finally, post-necking hardening behavior was successfully characterized up to the initial stage of localized necking using the VFM with Swift and modified Voce laws.     

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.     

Nonlocal gradient-dependent modeling of plasticity with anisotropic hardening

This work addresses the formulation of the thermodynamics of nonlocal plasticity using the gradient theory. The formulation is based on the nonlocality energy residual introduced by Eringen and Edelen (1972). Gradients are introduced for those variables associated with isotropic and kinematic hardening. The formulation applies to small strain gradient plasticity and makes use of the evanescent memory model for kinematic hardening. This is accomplished using the kinematic flux evolution as developed by Zbib and Aifantis (1988). Therefore, the present theory is a four nonlocal parameter-based theory that accounts for the influence of large variations in the plastic strain, accumulated plastic strain, accumulated plastic strain gradients, and the micromechanical evolution of the kinematic flux. Using the principle of virtual power and the laws of thermodynamics, thermodynamically-consistent equations are derived for the nonlocal plasticity yield criterion and associated flow rule. The presence of higher-order gradients in the plastic strain is shown to enhance a corresponding history variable which arises from the accumulation of the plastic strain gradients. Furthermore, anisotropy is introduced by plastic strain gradients in the form of kinematic hardening. Plastic strain gradients can be attributed to the net Burgers vector, while gradients in the accumulation of plastic strain are responsible for the introduction of isotropic hardening. The equilibrium between internal Cauchy stress and the microstresses conjugate to the higher-order gradients frames the yield criterion, which is obtained from the principle of virtual power. Microscopic boundary conditions, associated with plastic flow, are introduced to supplement the macroscopic boundary conditions of classical plasticity. The nonlocal formulation developed here preserves the classical assumption of local plasticity, wherein plastic flow direction is governed by the deviatoric Cauchy stress. The theory is applied to the problems of thin films on both soft and hard substrates. Numerical solutions are presented for bi-axial tension and simple shear loading of thin films on substrates.     

Constitutive modeling of strain range dependent cyclic hardening

In this paper, a new approach for constitutive modeling of strain range dependent cyclic hardening is proposed by extending the kinematic hardening model based on the critical state of dynamic recovery. It is assumed that isotropic, as well as kinematic, hardening consists of several parts, and that each part of isotropic hardening evolves when the corresponding part of kinematic hardening is in the critical state of dynamic recovery. The extended model is capable of simulating the cyclic hardening behavior in which different characteristics of cyclic hardening appear depending on strain range. The model is verified by simulating the relatively large cyclic straining tests of 304 stainless steel at ambient temperature, in which cyclic hardening does not stabilize before rupture if strain range exceeds a certain value. The model is further verified by predicting the history dependence of cyclic hardening under incremental cyclic loading and the maximum plastic strain dependence of strain hardening in cyclic tension.     

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.     

Critical state shear behavior of the soil-structure interface determined by discrete element modeling

The interface between soil and structure can be referred to as a soil-structure system, and its behavior plays an important role in many geotechnical engineering practices. In this study, results are presented from a series of monotonic direct shear tests performed on a sand-structure interface under constant normal stiffness using the discrete element method (DEM). Strain localization and dilatancy behavior of the interface is carefully examined at both macroscopic and microscopic scales. The effects of soil initial relative density and normal stress on the interface shear behavior are also investigated. The results show that a shear band progressively develops along the structural surface as shear displacement increases. At large shear displacement a unique relationship between stress ratio and void ratio is reached in the shear band for a certain normal stress, indicating that a critical state exists in the shear band. It is also found that the thickness and void ratio of the shear band at the critical state decreases with increasing normal stress. Comparison of the DEM simulation results with experimental results provides insight into the shear behavior of a sand-structure interface and offers a means for quantitative modeling of such interfaces based on the critical state soil mechanics.     

Kinematic hardening rules in uncoupled modeling for multiaxial ratcheting simulation

An earlier paper by the authors evaluated the performance of several coupled models in simulating a series of uniaxial and biaxial ratcheting responses. This paper evaluates the performance of various kinematic hardening rules in an uncoupled model for the same set of ratcheting responses. A modified version of the Dafalias–Popov uncoupled model has been demonstrated to perform well for uniaxial ratcheting simulation. However, its performance in multiaxial ratcheting simulation is significantly influenced by the kinematic hardening rules employed in the model. Performances of eight different kinematic hardening rules, when engaged with the modified Dafalias–Popov model, are evaluated against a series of rate-independent multiaxial ratcheting responses of cyclically stabilized carbon steels. The kinematic hardening rules proposed by Armstrong–Frederick, Voyiadjis–Sivakumar, Phillips, Tseng–Lee, Kaneko, Xia–Ellyin, Chaboche and Ohno–Wang are examined. The Armstrong–Frederick rule performs reasonably for one type of the biaxial ratcheting response, but fails in others. The Voyiadjis–Sivakumar rule and its constituents, the Phillips and the Tseng–Lee rules, can not simulate the biaxial ratcheting responses. The Kaneko rule, composed of the Ziegler and the prestress directions, and the Xia–Ellyin rule, composed of the Ziegler and Mroz directions, also fail to simulate the biaxial ratcheting responses. The Chaboche rule, with three decomposed Armstrong–Frederick rules, performs the best for the whole set of ratcheting responses. The Ohno–Wang rule performs well for the data set, except for one biaxial response where it predicts shakedown with subsequent reversal of ratcheting.     

Nonlinear behavior of circular plates with work hardening

Tests were performed on two simply supported plates of aluminum alloy 2024-0, under a central concentrated load, with peak deflections up to 2.6 times the thickness. The load was provided by a small-diameter hard-steel rod. The plates had diameter-to-thickness ratios (D/h) of 20 and 41. Measurements were made of load, deflections and strains; membrane and bending strains were calculated from the test data. The test data are presented in comparison with theoretical predictions generated by the Grumman-developed finiteelement-computer code PLANS, which includes material and geometric nonlinearities. The theoretical prediction was excellent for deflections, and generally good for strains, when the central force was represented by a line load around the loading rod's contact circle. Using a uniform pressure as the central force caused the theory to slightly overpredict the peak deflections and greatly overpredict the peak strains at the plate center. The plates exhibited initial loss of stiffness under the plastic-bending behavior, followed by a restiffening when the large deflections caused a rapidly increasing membrane action that provided much additional resistance to the applied load.     

Visco-plastic constitutive modeling on Ohno-Wang kinematic hardening rule for uniaxial ratcheting behavior of Z2CND18.12N steel

Experimental results of monotonic uniaxial tensile tests at different strain rates and the reversed strain cycling test showed the characteristics of rate-dependence and cyclic hardening of Z2CND18.12N austenitic stainless steel at room temperature, respectively. Based on the Ohno-Wang kinematic hardening rule, a visco-plastic constitutive model incorporated with isotropic hardening was developed to describe the uniaxial ratcheting behavior of Z2CND18.12N steel under various stress-controlled loading conditions. Predicted results of the developed model agreed better with experimental results when the ratcheting strain level became higher, but the developed model overestimated the ratcheting deformation in other cases. A modified model was proposed to improve the prediction accuracy. In the modified model, the parameter m_i of the Ohno-Wang kinematic hardening rule was developed to evolve with the accumulated plastic strain. Simulation results of the modified model proved much better agreement with experiments.     

Suspension-like hardening behavior of HDPE and time-hardening superposition

The rheology of solidifying high-density polyethylene (HDPE) is investigated. Experiments on an HDPE were performed with a novel RheoDSC device. Results agree quantitatively with simulations for a suspension of elastic spheres in a viscoelastic matrix except for very low values of space filling (<5%), indicating that the rheological behavior of the crystallizing melt in the frequency range investigated is purely suspension like. The hardening behavior of the material is characterized in two different ways; a normalized rheological function and a time-hardening superposition (THS) master curve of rheological properties. An improvement is proposed to the procedure for performing THS that was previously used in the literature. Based on this procedure, a novel method for predicting the rheological properties of crystallizing melts is presented.