Eulerian constitutive model for finite deformation plasticity with anisotropic hardening

A constitutive model is presented for finite strain plasticity. The model incorporates both isotropic and kinematic hardening of the Ziegler type. The corotational rate used here is in line with the theory suggested by Paulun and Pecherski (1985) but not necessarily confined to the von Mises type yield criterion and the Prager hardening rule. The aspect of integration of the corotational rates is also discussed here. The use of the integration of the material rate of tensors with time as a substitute for the proper integration with time of corotational rates leads to mathematical inconsistencies of the theory of Lie derivatives. The problem of simple shear is investigated and compared with other works.     

Internal-variable constitutive model for rate-independent plasticity with hardening saturation surface

Summary An elastic-plastic material model with internal variables and thermodynamic potential, not admitting hardening states out of a saturation surface, is presented. The existence of such a saturation surface in the internal variables space — a consequence of the boundedness of the energy that can be stored in the material's internal micro-structure — encompasses, in case of general kinematic/isotropic hardening, a one-parameter family of envelope surfaces in the stress space, which in turn is enveloped by a limit surface. In contrast to a multi-surface model, noad hoc rules are required to avoid the intersection between the yield and bounding/envelope surface. The flow laws of the proposed model are studied in case of associative plasticity with the aid of the maximum intrinsic dissipation theorem. It is shown that the material behaves like a standard one as long as its hardening state either is not saturated, or undergoes a desaturation from a saturated hardening state, whereas, for saturated hardening states not followed by desaturation, it conforms to a combined yielding law in which the static internal variable rates obey a nonlinear hardening rule similar to that of analogous models of the literature. Additionally, the material is shown to behave as a perfectly plastic material for a class of (critical) saturated hardening states for which the stress state is on the limit surface. For nonassociative material models, it is shown that, under a special choice of the plastic and saturation potentials and through a suitable parameter identification, the well-known Chaboche model is reproduced. A few numerical examples are presented to illustrate the associative material response under monotonic and cyclic loadings.Dedicated to Prof. Dr. Dr. h. c. Franz Ziegler on the occasion of his 60th birthday     

Effective constitutive equations for fiber-reinforced viscoplastic composites exhibiting anisotropic hardening

Effective elastic-viscoplastic stress-strain relations are derived for fiber-reinforced composites whose constituents are elastic-viscoplastic materials displaying anisotropic hardening. The derivation is based on a recently developed high-order continuum theory with microstructure for the modeling of viscoplastic composites, and is generalized here to incorporate anisotropic hardening effects. A specific reduction of the theory gives the effective rate-dependent elastic-plastic behavior of the composite which exhibits plastic anisotropy. In the special case of perfectly elastic constituents, the approximate overall moduli of the fiber-reinforced composite are obtained. Rate-dependent average stress-strain curves are given for numerous modes of cyclic loading of the composite. The effective behavior of periodically bilaminated viscoplastic composites is determined as a special case.     

On the derivation of the constitutive equations of ideal plasticity

Summary TheLèvy-Mises equations for an ideal plastic material are derived from a general implicit relation between the stress and rate of deformation tensors. This derivation provides an alternative to the usuala priori postulates of a yield criterion and flow rule.

---

Zur Herleitung der Materialgleichungen ideal plastischer Stoffe

Zusammenfassung DieLèvy-Mises-Gleichungen für ein ideal plastisches Material wurden aus einer allgemeinen impliziten Beziehung zwischen dem Spannungstensor und dem Tensor der Deformations-geschwindigkeiten hergeleitet. Diese Herleitung bietet eine Alternative zu den üblichena priori Postulaten einer Fließbedingung und eines Fließgesetzes.

     

Cyclic plastic behavior analysis based on the micromorphic mixed hardening plasticity model

In this paper, a small strain micromorphic elasto-plastic model with isotropic/kinematic hardening is presented for modeling the size effect and Bauschinger effect in material with microstructure. A nonlinear kinematic hardening model is embedded into the micromorphic framework by employing a backstress, a micro-backstress and a micro-couple-backstress in a physical way. The material intrinsic length scale is introduced in the constitutive law, leading to the presence of higher order stress. The present model is further implemented into a 2D plane strain finite element frame with a fully implicit stress integration scheme. The generalized consistent tangent modulus is derived to achieve the parabolic convergence of the global nodal force equilibrium equation. Two numerical examples, including a thin film and a plate with underlying structures subjected to cyclic loading, are analyzed to verify the theoretical developments and numerical formulations. Plastic behaviors in micromorphic continuum, such as size effect, Bauschinger effect, ratcheting effect and plastic shakedown phenomenon, are investigated.     

Combined hardening and softening constitutive model of plasticity: precursor to shear slip line failure

 In this work, we present a finite element model capable of describing both the plastic deformation which accumulates during the hardening phase as the precursor to failure and the failure process leading to softening phenomena induced by shear slip lines. This is achieved by activating subsequently hardening and softening mechanisms with the localization condition which separates them. The chosen model problem of von Mises plasticity is addressed in detail, along with particular combination of mixed and enhanced finite element approximations which are selected to control the locking phenomena and guarantee mesh-invariant computation of plastic dissipation. Several numerical simulations are presented in order to illustrate the ability of the presented model to predict the final orientation of the shear slip lines for the case of non-proportional loading. Dedicated to the memory of Prof. Mike Crisfield, for his cheerfulness and cooperation as a colleague and friend over many years. This work was supported by the French Ministry of Research and ACI research program. This support is gratefully acknowledged.     

On the numerical integration of a class of pressure-dependent plasticity models with mixed hardening

Radial return algorithms, for both three-dimensional and plane stress situations, are developed for a class of pressure-dependent plasticity models (formulated in state variables) with mixed hardening. The consistent tangent matrix has been developed which, among other advantages, does not require numerical inversion. The algorithms, for Gurson's mixed hardening model, are incorporated in a finite element program to solve several simple uniaxial tension problems. When compared with numerically integrated solutions, the results from the finite element analysis show excellent agreement. Finally, results of the axisymmetric necking problem, using Gurson's mixed hardening model, are presented as an example of application.     

Automatic consistency integrations for the stress update of J2 elastoplastic rate constitutive equations with combined hardening

A novel method is introduced to study numerical integrations of J₂ elastoplastic rate constitutive equations with general combined hardening. The basic idea is to transform the usual time rate constitutive equations into those with reference to the equivalent plastic strain. By virtue of tensorial matrix operations, we show that these transformed equations may be converted to a linear differential system governing the shifted stress and the plastic multiplier. From this system, we derive explicit integrations for the shifted stress and then for the back stress and the Cauchy stress. We demonstrate that these results are accurate up to within a third order term of the equivalent plastic strain increment. In particular, for pure kinematic hardening, we show that the integrations obtained can achieve automatic enforcement of both the plastic consistency condition and the loading condition, thus bypassing the numerical treatment of the latter two. Furthermore, we explain that, with the new algorithm for the stress update, the continuum tangent moduli may be used to ensure a quadratic rate of convergency in Newton's iteration scheme for the balance equation. Numerical examples suggest that the new algorithm may be more accurate and efficient than the widely used return algorithm. Copyright © 2014 John Wiley & Sons, Ltd.     

A nonuniform hardening plasticity model for concrete materials

For the most part, the developments of constitutive models for for the concrete by the theory of plasticity have in the past been made to scarch for a suitable failure surface. The initial yield surface is usually assumed to have the same shape as the failure surface but with a reduced size. The subsequent loading surfaces are then obtained by the uniform expansion of the initial one. This approach is found generally inadequate in predicting the deformational behavior of concrete for a wide range of loading conditions.

The present non-uniform hardening plasticity model adopts the most sophisticated failure model of Willam-Warnke or Hsieh-Ting-Chen as the bounding surface; assume an initial yield surface with a shape that is different from the failure surface; proposes a nonuniform hardening rule for the subsequent loading surfaces with a hydrostatic pressure and Lode-angle dependent plasticity modulus; and utilizes a nonassociated flow rule for a general formulation.

The work-hardening stress-strain behaviors of concrete based on the present model are found in good agreement with experimental results involving a wide range of stress states and different types of concrete material. The important features of inelastic behavior of concrete, including brittle failure in tension; ductile behavior in compression; hydrostatic sensitivities; and volumetric dilation under compressive loadings can all be represented by this improved constitutive model.     

Thermodynamic format and heat generation of isotropic hardening plasticity

Summary Heat generation due to plastic deformation of metals and steel is studied. Whereas in many investigations it is assumed that the fraction η of the plastic work transformed into heat is constant throughout the deformation process, the fraction η is here derived from thermodynamic considerations in a large-strain setting. It is shown that for elasto-plasticity the fraction η follows as a result of the choice of free energy, potential function and yield function. Taking the stress-strain response and the dissipative properties of the material as basis for calibration, it is shown that the thermodynamic framework of a thermoplastic material is non-unique for the general situation of non-associated plasticity. In the investigation conducted here, the mechanical response and the portion of the plastic work converted into heat (or into stored energy) during plastic deformations is predicted by means of isotropic hardening von Mises plasticity. It is shown that for a situation in which the internal variable is taken as the effective plastic, close fitting to experimental data requires a non-associated format of the evolution law for the internal variable.     

A new fixed‐point algorithm for hardening plasticity based on non‐linear mixed variational inequalities

Moving from the seminal papers of Han and Reddy, we propose a fixed‐point algorithm for the solution of hardening plasticity two‐dimensional problems. The continuous problem may be classified as a mixed non‐linear non‐differentiable variational inequality of the second type and is discretized by means of a truly mixed finite‐element scheme. One of the main peculiarities of our approach is the use of the composite triangular element of Johnson and Mercier for the approximation of the stress field. The non‐differentiability is coped with via regularization whereas the non‐linearity is approached with a fixed‐point iterative strategy. Numerical results are proposed that investigate the sensitivity of the approach with respect to the mesh size and the regularization parameter ε. The simplicity of the proposed fixed‐point scheme with respect to more classical return mapping approaches seems to represent one of the key features of our algorithm. Copyright © 2003 John Wiley & Sons, Ltd.     

A new strain hardening model for rate-dependent crystal plasticity

This paper presents a crystal plasticity based finite element analysis employing the new microstructure-based strain hardening model recently presented by Saimoto and Van Houtte (2011) [7] to simulate formability and texture evolution in the commercial aluminum alloy 5754. Simulations are performed to compare the predictive capability of the new hardening model against the common work hardening models using a rate-dependent plasticity formulation. The parameters in the numerical models are calibrated using the X-ray data published by Iadicola et al. (2008) [9] for the aluminum sheet alloy 5754. The predictions of the model for balanced biaxial tension and in-plane plane-strain tests are compared against experimental observations presented in Iadicola et al. (2008) [9]. It is concluded that the new model provides the best predictions of the large strain behavior of Aluminum sheet alloy 5754 subjected to various strain paths.     

A variational constitutive model for porous metal plasticity

This paper presents a variational formulation of viscoplastic constitutive updates for porous elastoplastic materials. The material model combines von Mises plasticity with volumetric plastic expansion as induced, e.g., by the growth of voids and defects in metals. The finite deformation theory is based on the multiplicative decomposition of the deformation gradient and an internal variable formulation of continuum thermodynamics. By the use of logarithmic and exponential mappings the stress update algorithms are extended from small strains to finite deformations. Thus the time-discretized version of the porous-viscoplastic constitutive updates is described in a fully variational manner. The range of behavior predicted by the model and the performance of the variational update are demonstrated by its application to the forced expansion and fragmentation of U-6%Nb rings.     

Rate-type constitutive equations and elastic-plastic materials

This paper is concerned with a class of rate-type constitutive equations, which characterize the behavior of elastic-plastic solids under large deformation. The theory is developed without à priori introduction of the concept of plastic strain and involves only the stress and the total strain as basic ingredients, apart from the thermodynamic variables. The relationship between the results obtained here and those of a more familiar form of a theory of elastic-plastic solids is indicated.     

Incremental constitutive equation for large strain elasto plasticity

Starting from the standard theory with intermediate configuration for finite deformations of an isotropic elasto-plastic material with isotropic hardening, rate type constitutive equations are obtained. The small elastic strain approximation is then discussed and it is shown that, in this approximation, these equations reduce to Hill's formalism of large strain elasto-plasticity obtained from the classical Prandtl-Reuss relations of infinitesimal plasticity by substituting for the infinitesimal strain rate, stress and stress rate respectively the rate of deformation tensor, the Cauchy stress tensor and the Jaumann stress rate tensor. The limiting case of perfect plasticity is also investigated.