Application of corotational rates of the logarithmic strain in constitutive modeling of hardening materials at finite deformations

In this paper a finite deformation constitutive model for rigid plastic hardening materials based on the logarithmic strain tensor is introduced. The flow rule of this constitutive model relates the corotational rate of the logarithmic strain tensor to the difference of the deviatoric Cauchy stress and the back stress tensors. The evolution equation for the kinematic hardening of this model relates the corotational rate of the back stress tensor to the corotational rate of the logarithmic strain tensor. Using Jaumann, Green–Naghdi, Eulerian and logarithmic corotational rates in the proposed constitutive model, stress–strain responses and subsequent yield surfaces are determined for rigid plastic kinematic and isotropic hardening materials in the simple shear problem at finite deformations.     

The appropriate corotational rate, exact formula for the plastic spin and constitutive model for finite elastoplasticity

The exact formulae for the plastic and the elastic spin referred to the deformed configuration are derived, where the plastic spin is a function of the plastic strain rate and the elastic spin a function of the elastic strain rate. With these exact formulae we determine the macroscopic substructure spin that allows us to define the appropriate corotational rate for finite elastoplasticity.Plastic, elastic and substructure spin are considered and simplified for various sub-classes of restricted elastic-plastic strains. It is shown that for the special cases of rigid-plasticity and hypoelasticity the proposed corotational rate is identical with the Green-Naghdi rate, while the ZarembaJaumann rate yields a good approximation for moderately large strains.To compare our exact plastic spin formula with the constitutive assumption for the plastic spin introduced by Dafalias and others, we simplify our result for small elastic-moderate plastic strains and introduce a simplest evolution law for kinematic hardening leading to the Dafalias formula and to an exact determination of its unknown coefficient. It is also shown that contrary to statements in the literature the plastic spin is not zero for vanishing kinematic hardening.For isotropic-elastic material with induced plastic flow undergoing isotropic and kinematic hardening constitutive and evolution laws are proposed. Elastic and plastic Lagrangean and Eulerian logarithmic strain measures are introduced and their material time derivatives and corotational rates, respectively, are considered. Finally, the elastic-plastic tangent operator is derived.The presented theory is implemented in a solution algorithm and numerically applied to the simple shear problem for finite elastic-finite plastic strains as well as for sub-classes of restricted strains. The results are compared with those of the literature and with those obtained by using other corotational rates.     

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.     

Self-consistent Eulerian rate type elasto-plasticity models based upon the logarithmic stress rate

The objective of this article is to suggest new Eulerian rate type constitutive models for isotropic finite deformation elastoplasticity with isotropic hardening, kinematic hardening and combined isotropic-kinematic hardening etc. The main novelty of the suggested models is the use of the newly discovered logarithmic stress rate and the incorporation of a simple, natural explicit integrable-exactly rate type formulation of general hyperelasticity. Each new model is thus subjected to no incompatibility of rate type formulation for elastic behaviour with the notion of elasticity, as encountered by any other existing Eulerian rate type model for elastoplasticity or hypoelasticity. As particular cases, new Prandtl-Reuss equations for elastic-perfect plasticity and elastoplasticity with isotropic hardening, kinematic hardening and combined isotropic-kinematic hardening, respectively, are presented for computational and practical purposes. Of them, the equations for kinematic hardening and combined isotropic–kinematic hardening are, respectively, reduced to three uncoupled equations with respect to the spherical stress component, the shifted stress and the back-stress. The effects of finite rotation on the current strain and stress and hardening behaviour are indicated in a clear and direct manner. As illustrations, finite simple shear responses for the proposed models are studied by means of numerical integration. Further, it is proved that, among all possible (infinitely many) objective Eulerian rate type models, the proposed models are not only the first, but unique, self-consistent models of their kinds, in the sense that the rate type equation used to represent elastic behaviour is exactly integrable to really deliver an elastic relation. ©     

Constitutive relations for isotropic or kinematic hardening at finite elastic–plastic deformations

The rate-type constitutive relations of rate-independent metals with isotropic or kinematic hardening at finite elastic–plastic deformations were presented through a phenomenological approach. This approach includes the decomposition of finite deformation into elastic and plastic parts, which is different from both the elastic–plastic additive decomposition of deformation rate and Lee’s elastic–plastic multiplicative decomposition of deformation gradient. The objectivity of the constitutive relations was dealt with in integrating the constitutive equations. A new objective derivative of back stress was proposed for kinematic hardening. In addition, the loading criteria were discussed. Finally, the stress for simple shear elastic–plastic deformation was worked out.     

Large strain responses of elastic-perfect plasticity and kinematic hardening plasticity with the logarithmic rate: Swift effect in torsion

A new Eulerian rate type elastic-perfectly plastic model has recently been established by utilizing the newly discovered logarithmic rate. It has been proved that this model is unique among the objective elastic-perfectly plastic models with all objective corotational stress rates and other known objective stress rates by virtue of the self-consistency criterion: the hypoelastic formulation intended for elastic behaviour must be exactly integrable to deliver a hyperelastic relation. The finite simple shear response of this model has been studied and shown to be reasonable for both shear and normal stress components. On the other hand, a kinematic hardening plasticity model may be formulated by adopting the logarithmic rate. The objective of this work is to further study the large deformation responses of the foregoing two kinds of idealized models, in particular the well-known Swift effect, in torsion of thin-walled cylindrical tubes. A complete, rigorous analysis is made for the orders of magnitude of all stress components. A closed-form solution is obtained for the kinematic hardening plastic case, and an analytical perturbation solution is derived for the elastic-perfectly plastic case. It is shown that the simple idealized kinematic hardening model with the logarithmic rate, which uses only two classical material constants, i.e., the initial (tensile) yield stress and the hardening modulus, may arrive at satisfactory explanation for and reasonable accord with salient features of experimental observation.     

Anisotropic finite elastoplasticity with nonlinear kinematic and isotropic hardening and application to sheet metal forming

The paper discusses the derivation and the numerical implementation of a finite strain material model for plastic anisotropy and nonlinear kinematic and isotropic hardening. The model is derived from a thermodynamic framework and is based on the multiplicative split of the deformation gradient in the context of hyperelasticity. The kinematic hardening component represents a continuum extension of the classical rheological model of Armstrong–Frederick kinematic hardening. Introducing the so-called structure tensors as additional tensor-valued arguments, plastic anisotropy can be modelled by representing the yield surface and the plastic flow rule as functions of the structure tensors. The evolution equations are integrated by a new form of the exponential map that preserves plastic incompressibility and uses the spectral decomposition to evaluate the exponential tensor functions in closed form. Finally, the applicability of the model is demonstrated by means of simulations of several deep drawing processes and comparisons with experiments.     

A coupled temperature and strain rate dependent yield function for dynamic deformations of bcc metals

A coupled temperature and strain rate microstructure physically based yield function is proposed in this work. It is incorporated along with the Clausius–Duhem inequality and an appropriate free energy definition in a general thermodynamic framework for deriving a three-dimensional kinematical model for thermo-viscoplastic deformations of body centered cubic (bcc) metals. The evolution equations are expressed in terms of the material time derivatives of the elastic strain, accumulated plastic strain (isotropic hardening), and the back stress conjugate tensor (kinematic hardening). The viscoplastic multipliers are obtained using both the Consistency and Perzyna viscoplasticity models. The athermal yield function is employed instead of the static yield function in the case of the Perzyna viscoplasticity model. It is found that the static strain rate value, at which the material shows rate-independent behavior, varies with the material deformation temperature. Computational aspects of the proposed model are addressed through the finite element implementation with an implicit stress integration algorithm. Finite element simulations are performed by implementing the proposed viscoplasticity constitutive models in the commercial finite element program ABAQUS/Explicit [ABAQUS, 2003. User Manual, Version 6.3. Habbitt, Karlsson and Sorensen Inc., Providence, RI] via the user material subroutine coded as VUMAT. Numerical implementation for a simple compression problem meshed with one element is used to validate the proposed model implementation with applications to tantalum, niobium, and vanadium at low and high strain rates and temperatures. The analysis of a tensile shear banding is also investigated to show the effectiveness and the performance of the proposed framework in describing the strain localizations at high velocity impact. Results show mesh independency as a result of the viscoplastic regularization used in the proposed formulation.     

On the application of the plastic spin concept for the description of anisotropic hardening in finite deformation plasticity

The aim of the paper is to fill the gap between the general theoretical formulation of the constitutive relations for plastic spin and practical applications for proper prediction of material behavior at finite plastic deformations and anisotropic hardening. An approximation to the representation of the general constitutive equation for plastic spin is considered and the pertinent substructure corotational rate is applied to formulate the relation for rigid-plastic material with kinematic hardening. The simple shear traction problem is analysed and the proposed model is verified with the experimental results of Swift. The merits of the present proposal vis-à-vis the existing theories are discussed.     

Ordinary state-based peridynamics for plastic deformation according to von Mises yield criteria with isotropic hardening

This study presents the ordinary state-based peridynamic constitutive relations for plastic deformation based on von Mises yield criteria with isotropic hardening. The peridynamic force density–stretch relations concerning elastic deformation are augmented with increments of force density and stretch for plastic deformation. The expressions for the yield function and the rule of incremental plastic stretch are derived in terms of the horizon, force density, shear modulus, and hardening parameter of the material. The yield surface is constructed based on the relationship between the effective stress and equivalent plastic stretch. The validity of peridynamic predictions is established by considering benchmark solutions concerning a plate under tension, a plate with a hole and a crack also under tension.     

Kinematic hardening in large strain plasticity

A finite strain hyper elasto-plastic constitutive model capable to describe non-linear kinematic hardening as well as non-linear isotropic hardening is presented. In addition to the intermediate configuration and in order to model kinematic hardening, an additional configuration is introduced – the center configuration; both configurations are chosen to be isoclinic. The yield condition is formulated in terms of the Mandel stress and a back-stress with a structure similar to the Mandel stress.It is shown that the non-dissipative part of the plastic velocity gradient not governed by the thermodynamical framework and the corresponding quantity associated with the kinematic hardening influence the material behaviour to a large extent when kinematic hardening is present. However, for isotropic elasticity and isotropic hardening plasticity it is shown that the non-dissipative quantities have no influence upon the stress–strain relation.As an example, kinematic hardening von Mises plasticity is considered, which fulfils the plastic incompressibility condition and is independent of the hydrostatic pressure. To evaluate the response and to examine the influence of the non-dissipative quantities, simple shear is considered; no stress oscillations occur.     

模型复合材料弹塑性界面应力分析

由纤维增强弹塑性基体所产生的界面具有弹塑性力学行为。考虑到一般材料的塑性变形都遵循幂硬化规律,对模型复合材料的界面进行弹性和应变硬化状态下的变形规律及其应力分析。以纤维拔出试验为研究模型,将界面分成弹性区和塑性区。利用界面应力剪滞理论,分别建立弹性区和塑性区的界面力学基本方程。选择适当的位移函数满足基本方程及埋入纤维的边界条件,再按位移函数求出弹性区和塑性区的界面剪应力。推导出平均界面剪应力与纤维     

Strain-rate effects in rheological models of inelastic response

Stress–strain response under constant and variable strain-rate is studied for selected models of inelastic behavior. The derived closed-form solutions for uniaxial loading enable simple evaluation of the strain-rate effects on the material response. The effect of an abrupt change of strain-rate is also examined. Non-Newtonian viscosity which decreases with an increasing strain-rate is incorporated in the analysis. Parabolic and hyperbolic hardening are used to describe the plastic response in monotonic loading. A three-dimensional generalization of an elastic–viscoplastic model is employed to study the stress relaxation in simple shear. A combined isotropic–kinematic hardening and the concept of overstress are used in the analysis. The unloading nonlinearity of the stress–strain curve is then discussed.     

Directional distortional hardening at large plastic deformations

This paper extents the directional distortional hardening model of Feigenbaum and Dafalias (2007) into the range of large plastic deformations. This model allows the yield surface to deform such that a region of high curvature develops approximately in the direction of loading and a region of flattening develops on the opposite side. To extend this model into large deformations and in order to ensure positive dissipation and objectivity, hardening rules are derived from thermodynamic conditions in terms of corotational rates. Since this model includes a fourth order tensor-valued hardening internal variable, the corotational rates for fourth order tensors are examined in this work employing the concept of plastic spin. Several choices for plastic spins are presented and used for the simulation of the response under simple shear loading up to 1000% strain.     

A critical state sand model with elastic–plastic coupling

Soil elastic moduli are highly pressure-dependent. Experimental findings have indicated that the elastic shear modulus of sands depends on p^χ, where p is mean principal effective stress and χ is a non-dimensional parameter. χ practically remains unchanged for shear strains less than 10⁻⁵ where the mechanical behavior is purely elastic. However, experiments have revealed that the emergence of plasticity for shear strains larger than 10⁻⁵ provokes a gradual increase in χ. Technically, this observation is an elastic–plastic coupling effect in which plasticity causes to change the elastic characteristics. Here, this issue is considered in hyper-elasticity framework in conjunction with a critical state compatible bounding surface plasticity platform for granular soils. To this aim, constitutive equations linking χ to a proper kinematic hardening parameter are presented. Then, using the proposed approach, a hyper-elastic theory is modified to consider the mentioned elastic–plastic coupling effect in the whole domain of the elastoplastic behavior. Adopting the improved hyper-elasticity necessitates the modification of a number of basic plasticity platform elements. In this regard, dilatancy and plastic hardening modulus of the bounding surface platform are modified. Successful performance of the modified constitutive model is presented against experimental data of loading/unloading triaxial tests.     

Dual finite element methods in homogenization for elastic–plastic fibrous composite material

The problem of homogenization for a periodic, elastic–perfectly plastic, fiber reinforced, composite material is considered. The overall mechanical behavior of the material is described using the anisotropic model of elastic–plastic body with kinematic hardening. The appropriate initial–boundary value problem, set for one repeatable cell of the composite, is solved in order to find the effective constitutive relations. The cell problem is solved using the finite element method formulated in two dual forms: in displacements and in stresses. Stress functions are used in the latter formulation.     

Large Strain Response of Kinematic Hardening Elastoplasticity with the Logarithmic Rate: Swift Effect in Torsion

Recently, a new Eulerian rate type kinematic hardening elastoplasticity model has been established by utilizing the newly discovered logarithmic rate. It has been proved that this model is unique among all the objective kinematic hardening elastoplastic models with all possible objective corotational stress rates and other known objective stress rates by virtue of the self-consistency criterion: the hypoelastic formulation intended for elastic behaviour must be exactly integrable to deliver a hyperelastic relation. The finite simple shear response of this model has been studied and shown to be reasonable for both shear and normal stress components. The objective of this work is to further study the large deformation response of this model, in particular the well-known Swift effect, in torsion of thin-walled cylindrical tubes with free ends. Analytical perturbation solution and numerical solution are presented for the case of linear kinematic hardening at large compressible deformation. It is shown that the prediction from the foregoing model is in good accord with experimental data reported in literature.     

A three-invariant cap model with isotropic–kinematic hardening rule and associated plasticity for granular materials

In this paper, a three-invariant cap model is developed for the isotropic–kinematic hardening and associated plasticity of granular materials. The model is based on the concepts of elasticity and plasticity theories together with an associated flow rule and a work hardening law for plastic deformations of granulars. The hardening rule is defined by its decomposition into the isotropic and kinematic material functions. The constitutive elasto-plastic matrix and its components are derived by using the definition of yield surface, material functions and non-linear elastic behavior, as function of hardening parameters. The model assessment and procedure for determination of material parameters are described. Finally, the applicability of proposed plasticity model is demonstrated in numerical simulation of several triaxial and confining pressure tests on different granular materials, including: wheat, rape, synthetic granulate and sand.     

A nonlocal strain gradient plasticity theory for finite deformations

Strain gradient plasticity for finite deformations is addressed within the framework of nonlocal continuum thermodynamics, featured by the concepts of (nonlocality) energy residual and globally simple material. The plastic strain gradient is assumed to be physically meaningful in the domain of particle isoclinic configurations (with the director vector triad constant both in space and time), whereas the objective notion of corotational gradient makes it possible to compute the plastic strain gradient in any domain of particle intermediate configurations. A phenomenological elastic–plastic constitutive model is presented, with mixed kinematic/isotropic hardening laws in the form of PDEs and related higher order boundary conditions (including those associated with the moving elastic/plastic boundary). Two fourth-order projection tensor operators, functions of the elastic and plastic strain states, are shown to relate the skew-symmetric parts of the Mandel stress and back stress to the related symmetric parts. Consistent with the thermodynamic restrictions therein derived, the flow laws for rate-independent associative plasticity are formulated in a six-dimensional tensor space in terms of symmetric parts of Mandel stresses and related work-conjugate generalized plastic strain rates. A simple shear problem application is presented for illustrative purposes.     

Description of plastic anisotropy effects at large deformations. Part II: the case of transverse isotropy

In a previous paper (see Tsakmakis, 1999) a general thermodynamically consistent (visco-) plasticity theory has been developed, which accounts for anisotropy effects. For simplicity, isotropic hardening has not be regarded, while anisotropy arises from kinematic hardening and orientational evolution of the underlying substructure. In the present paper the capabilities of this theory are discussed for the study case of transverse isotropy. Anisotropy effects are elaborated in the free energy and the yield function by means of structural tensors. Characteristic features of the transversely isotropic model are illustrated for the case of homogeneous simple shear.