Corotational rates in constitutive modeling of elastic-plastic deformation

The principal axes technique is used to develop a new hypoelastic constitutive model for an isotropic elastic solid in finite deformation. The new model is shown to produce solutions that are independent of the choice of objective stress rate. In addition, the new model is found to be equivalent to the isotropic finite elastic model; this is essential if both models describe the same material.

The new hypoelastic model is combined with an isotropic flow rule to form an elastic-plastic rate constitutive equation. Use of the principal axes technique ensures that the stress tensor is coaxial with the elastic stretch tensor and that solutions do not depend on the choice of objective stress rate. The flow rule of von Mises and a parabolic hardening law are used to provide an example of application of the new theory. A solution is obtained for the prescribed deformation of simple rectilinear shear of an isotropic elastic and isotropic elastic-plastic material.     

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.     

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.     

Constitutive modeling of temperature and strain rate dependent elastoplastic hardening materials using a corotational rate associated with the plastic deformation

In this paper, a constitutive model with a temperature and strain rate dependent flow stress (Bergstrom hardening rule) and modified Armstrong-Frederick kinematic evolution equation for elastoplastic hardening materials is introduced. Based on the multiplicative decomposition of the deformation gradient,new kinematic relations for the elastic and plastic left stretch tensors as well as the plastic deformation-dependent spin tensor are proposed. Also, a closed-form solution has been obtained for the elastic and plastic left stretch tensors for the simple shear problem.To evaluate model validity, results are compared with known experimental data for SUS 304 stainless steel, which shows a good agreement with the results of the proposed theoretical model.Finally, the stress-deformation curve, as predicted by the model, is plotted for the simple shear problem at room and elevated temperatures using the same material properties for AA5754-O aluminium alloy.     

Modelling inelastic behaviour of orthotropic metals in a unique alignment of deviatoric plane within the stress space

A finite strain constitutive model to predict the deformation behaviour of orthotropic metals is developed in this paper. The important features of this constitutive model are the multiplicative decomposition of the deformation gradient and a new Mandel stress tensor combined with the new stress tensor decomposition generalized into deviatoric and spherical parts. The elastic free energy function and the yield function are defined within an invariant theory by means of the structural tensors. The Hill’s yield criterion is adopted to characterize plastic orthotropy, and the thermally micromechanical-based model, Mechanical Threshold Model (MTS) is used as a referential curve to control the yield surface expansion using an isotropic plastic hardening assumption. The model complexity is further extended by coupling the formulation with the shock equation of state (EOS). The proposed formulation is integrated in the isoclinic configuration and allows for a unique treatment for elastic and plastic anisotropy. The effects of elastic anisotropy are taken into account through the stress tensor decomposition and plastic anisotropy through yield surface defined in the generalized deviatoric plane perpendicular to the generalized pressure. The proposed formulation of this work is implemented into the Lawrence Livermore National Laboratory-DYNA3D code by the modification of several subroutines in the code. The capability of the new constitutive model to capture strain rate and temperature sensitivity is then validated. The final part of this process is a comparison of the results generated by the proposed constitutive model against the available experimental data from both the Plate Impact test and Taylor Cylinder Impact test. A good agreement between experimental and simulation is obtained in each test.     

On the mechanics of solids with a growing mass

A general constitutive theory of the stress-modulated growth of biomaterials is presented with a particular accent given to pseudo-elastic soft living tissues. The governing equations of the mechanics of solids with a growing mass are revisited within the framework of finite deformation continuum thermodynamics. The multiplicative decomposition of the deformation gradient into its elastic and growth parts is employed to study the growth of isotropic, transversely isotropic, and orthotropic biomaterials. An explicit representation of the growth part of the deformation gradient is given in each case, which leads to an effective incremental formulation in the analysis of the stress-modulated growth process. The rectangular components of the instantaneous elastic moduli tensor are derived corresponding to selected forms of the elastic strain energy function. Physically appealing structures of the stress-dependent evolution equations for the growth induced stretch ratios are proposed.     

A note on nonassociated plastic flow rules

The analytical properties of the constitutive equations in plasticity with a nonassociated flow rule are investigated. Under the assumption of small deformations the directional stiffness (and compliance) rule is considered and the relevant spectral properties of the tangent stiffness tensor are assessed. It is shown that the directional stiffness may be larger than the elastic. It may also be negative in the case of a formally perfectly plastic material and so the nonassociative flow rule represents (spurious) softening in terms of an associated flow rule. The issue of uniqueness at finite strains is briefly addressed, whereby use is made of the tangent stiffness tensor relating the velocity gradient to the first Piola-Kirchhoff stress rate. The relevant spectral properties, which generalise those from the small deformation case, are found explicit. A sufficient condition for uniqueness is given in terms of a critical (upper bound) value of the hardening modulus.     

A class of universal relations in isotropic elasticity theory

A class of universal relations for isotropic elastic materials is described by the tensor equationTB = BT. This simple rule yields at most three component relations which are the generators of many known universal relations for isotropic elasticity theory, including the well-known universal rule for a simple shear. Universal relations for four families of nonhomogeneous deformations known to be controllable in every incompressible, homogeneous and isotropic elastic material are exhibited. These same universal relations may hold for special compressible materials. New universal relations for a homogeneous controllable shear, a nonhomogeneous shear, and a variable extension are derived. The general universal relation for an arbitrary isotropic tensor function of a symmetric tensor also is noted.     

On the Eigenvalues of the Fourth-Order Constitutive Tensor and Loss of Strong Ellipticity in Elastoplasticity

The eigenvalues of the fourth-order constitutive tangent modulus and the corresponding acoustic tensors are analyzed. Explicit expressions of the eigenvalues are made for the nonsymmetric tangent modulus tensor, and in the case of the deviatoric associative rule for the symmetric part of the tangent modulus and its acoustic tensor. In this context, a rate independent infinitesimal elastoplastic model is considered. The expressions of the plastic hardening modulus are summarized for the different local stability criteria (loss of second order work positiveness, loss of ellipticity, and loss of strong ellipticity). The critical hardening modulus and orientation are discussed in detail in the case of loss of ellipticity and loss of strong ellipticity. This analysis is based on the geometric method and linear, isotropic elasticity and deviatoric associative flow rule. In particular, the critical orientation for the loss of strong ellipticity and the classical shear band localization are compared.     

A formulation of anisotropic continuum elastoplasticity at finite strains. Part I: Modelling

A constitutive model for anisotropic elastoplasticity at finite strains is developed together with its numerical implementation. An anisotropic elastic constitutive law is described in an invariant setting by use of structural tensors and the elastic strain measure C_e. The elastic strain tensor as well as the structural tensors are assumed to be invariant in relation to superimposed rigid body rotations. An anisotropic Hill-type yield criterion, described by a non-symmetric Eshelby-like stress tensor and further structural tensors, is developed, where use is made of representation theorems for functions with non-symmetric arguments. The model also considers non-linear isotropic hardening. Explicit results for the specific case of orthotropic anisotropy are given. The associative flow rule is employed and the features of the inelastic flow rule are discussed in full. It is shown that the classical definition of the plastic material spin is meaningless in conjunction with the present formulation. Instead, the study motivates an alternative definition, which is based on the demand that such a quantity must be dissipation-free, as the plastic material spin is in the case of isotropy. Equivalent spatial formulations are presented too. The full numerical treatment is considered in Part II.     

A finite deformation formulation of the 3-parameter viscoelastic fluid

The paper addresses, with reference to the 3-parameter fluid, a method to extend constitutive laws for viscoelastic fluids from small to finite deformations. It relies upon the multiplicative decomposition of the deformation gradient tensor into elastic and inelastic parts. Also care is taken that the second law of thermodynamics is satisfied in every admissible process. To demonstrate the capabilities of the model obtained, simple shear and uniaxial extensional flow are discussed. It turns out that essential effects of material behaviour, which have been observed experimentally, can be predicted by the model.     

An Analysis of Professor J.F. Bell's Research on the Finite Twist and Extension of Thin-walled Polycrystalline Cylindrical Tubes

Professor J.F. Bell's empirical result regarding the rotation factor in the polar decomposition of the deformation gradient for the finite twist–extension of a thin-walled polycrystalline cylindrical metal tube is examined. The correct expression for the rotation is derived and used to show how Bell's result should be interpreted. Some implications for his incremental plasticity equations are also discussed. In particular, they are shown to satisfy appropriate invariance requirements when cast in terms of the variables actually measured by Bell in his experiments. Further consequences of his equations consistent with his data are also derived. Finally, it is shown that his theory furnishes a consistent constitutive statement about the response of isotropic solids provided that the Cauchy stress is constrained to be symmetric.     

Stability of deformation of isotropic hyperelastic bodies

The equations relating stress rates to strain rates are formulated and conditions of stable deformation of isotropic hyperelastic bodies are described. Stress–strain relations are presented for pure shear and uniaxial and axisymmetric loading of a material with a constitutive function obtained by generalization of the constitutive function of Hooke's law. In the case of small strains, the diagrams virtually coincide with the linear diagrams following from Hooke's law. Ramification of solutions and transition to declining diagrams begin at the same time, irrespective of values of the constants of the material, when large stresses of the order of the shear modulus are reached.     

A novel internal dissipation inequality by isotropy and its implication for inelastic constitutive characterization

In the present work a novel inelastic deformation caused internal dissipation inequality by isotropy is revealed. This inequality has the most concise form among a variety of internal dissipation inequalities, including the one widely used in constitutive characterization of isotropic finite strain elastoplasticity and viscoelasticiy. Further, the evolution term describing the difference between the rate of deformation tensor and the “principal rate” of the elastic logarithmic strain tensor is set, according to the standard practice by isotropy, to equal a rank-two isotropic tensor function of the corresponding branch stress, with the tensor function having an eigenspace identical to the eigenspace of the branch stress tensor. Through that a general form of evolution equation for the elastic logarithmic strain is formulated and some interesting and important results are derived. Namely, by isotropy the evolution of the elastic logarithmic strain tensor is embodied separately by the evolutions of its eigenvalues and eigenprojections, with the evolution of the eigenprojections driven by the rate of deformation tensor and the evolution of the eigenvalues connected to specific material behavior. It can be proved that by isotropy the evolution term in the present dissipation inequality stands for the essential form of the evolution term in the extensively applied dissipation inequality.     

A comparative study of kinematic hardening rules at finite deformations

Kinematic hardening models describe a specific kind of plastic anisotropy which evolves with the deformation process. It is well known that the extension of constitutive relations from small to finite deformations is not unique. This applies also to well-established kinematic hardening rules like that of Armstrong-Frederick or Chaboche. However, the second law of thermodynamics offers some possibilities for generalizing constitutive equations so that this ambiguity may, in some extent, be moderated. The present paper is concerned with three possible extensions, from small to finite deformations, of the Armstrong-Frederick rule, which are derived as sufficient conditions for the validity of the second law. All three models rely upon the multiplicative decomposition of the deformation gradient tensor into elastic and plastic parts and make use of a yield function expressed in terms of the so-called Mandel stress tensor. In conformity with this approach, the back-stress tensor is defined to be of Mandel stress type as well. In order to compare the properties of the three models, predicted responses for processes with homogeneous and inhomogeneous deformations are discussed. To this end, the models are implemented in a finite element code (ABAQUS).     

New universal relations for nonlinear isotropic elastic materials

A nonlinear isotropic elastic block is subjected to a homogeneous deformation consisting of simple shear superposed on triaxial extension. Two new relations are established for this deformation which are valid for all nonlinear elastic isotropic materials, and hence are universal relations. The first is a relation between the stretch ratios in the plane of shear and the amount of shear when the deformation is supported only by shear tractions. The second relation is established for a thin-walled cylinder under combined extension, inflation and torsion. Each material element of the cylinder undergoes the same local homogeneous deformation of shear superposed on triaxial extension. The properties of this deformation are used to establish a relation between pressure, twisting moment, angle of twist and current dimensions when no axial force is applied to the cylinder. It is shown that these relations also apply for a mixture of a nonlinear isotropic solid and a fluid.     

Finite Bending of a Rectangular Block of an Elastic Hencky Material

Hencky's elasticity model is an isotropic, finite hyperelastic equation obtained by simply replacing the Cauchy stress tensor and the infinitesimal strain tensor in the classical Hooke's law for isotropic infinitesimal elasticity with the Kirchhoff stress tensor and Hencky's logarithmic strain tensor. A study by Anand in 1979 and 1986 indicates that it is a realistic finite elasticity model that is in good accord with experimental data for a variety of engineering materials for moderate deformations. Most recently, by virtue of well-founded physical grounds and rigorous mathematical procedures it has been demonstrated by these authors that this model may be essential to achieving self-consistent Eulerian rate type theories of finite inelasticity, e.g., the J ₂-flow theory for metal plasticity, etc. Its predictions have been studied for some typical deformation modes, including extension, simple shear and torsion, etc. Here we are concerned with finite bending of a rectangular block. We show that a closed-form solution may be obtained. We present explicit expressions for the bending angle and the bending moment in terms of the maximum or minimum circumferential stretch in a general case of compressible deformations for any assigned stretch normal to the bending plane. In particular, simplified results are derived for the plane strain case and for the case of incompressibility. This revised version was published online in June 2006 with corrections to the Cover Date.     

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.     

The significance of pure measures of distortion in nonlinear elasticity with reference to the Poynting problem

The objective of this paper is to show the significance of expressing the strain energy function in terms of a scalar pure measure of dilatation and a tensor pure measure of distortion, which were essentially introduced by Flory [1]. It is shown that convenient representations for the strain energies of dilatation and distortion, and the pressure and deviatoric Cauchy stress may be recorded in terms of these deformation measures. After specializing to the case of an isotropic material, specific constitutive equations are proposed and the Poynting problem is considered. It is shown that the Poynting effect (extension of a bar in torsion) is significantly influenced by coupling between dilatational and distortional strain energies, which is caused by the dependence of the shear modulus on dilatation.