Hypo-Elasticity Model Based upon the Logarithmic Stress Rate

Recently these authors have proved [46, 47] that a smooth spin tensor Ω^log can be found such that the stretching tensor D can be exactly written as an objective corotational rate of the Eulerian logarithmic strain measure ln V defined by this spin tensor, and furthermore that in all strain tensor measures only ln V enjoys this favourable property. This spin tensor is called the logarithmic spin and the objective corotational rate of an Eulerian tensor defined by it is called the logarithmic tensor-rate. In this paper, we propose and investigate a hypo-elasticity model based upon the objective corotational rate of the Kirchhoff stress defined by the spin Ω^log, i.e. the logarithmic stress rate. By virtue of the proposed model, we show that the simplest relationship between hypo-elasticity and elasticity can be established, and accordingly that Bernstein's integrability theorem relating hypo-elasticity to elasticity can be substantially simplified. In particular, we show that the simplest form of the proposed model, i.e. the hypo-elasticity model of grade zero, turns out to be integrable to deliver a linear isotropic relation between the Kirchhoff stress and the Eulerian logarithmic strain ln V, and moreover that this simplest model predicts the phenomenon of the known hypo-elastic yield at simple shear deformation. This revised version was published online in August 2006 with corrections to the Cover Date.     

On The Objective Corotational Rates of Eulerian Strain Measures

In the present paper, some new basis-free expressions for an arbitrary objective corotational rate of the general Eulerian strain measures are provided which are in compact form. Moreover, a complete discussion on the requirements for the continuity of the objective corotational rates are presented.      

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.     

Objective Tensor Rates and Applications in Formulation of Hyperelastic Relations

Following Ogden, a class of objective (Lagrangian and Eulerian) tensors is identified among the second-rank tensors characterizing continuum deformation, but a more general definition of objectivity than that used by Ogden is introduced. Time rates of tensors are determined using convective rates. Sufficient conditions of objectivity are obtained for convective rates of objective tensors. Objective convective rates of strain tensors are used to introduce pairs of symmetric stress and strain tensors conjugate in a generalized sense. The classical definitions of conjugate Lagrangian (after Hill) and Eulerian (after Xiao et al.) stress and strain tensors are particular cases of the definition of conjugacy of stress and strain tensors in the generalized sense used in the present paper. Pairs of objective stress and strain tensors conjugate in the generalized sense are used to formulate constitutive relations for a hyperelastic medium. A family of objective generalized strain tensors is introduced, which is broader than Hill’s family of strain tensors. The basic forms of the hyperelastic constitutive relations are obtained with the aid of pairs of Lagrangian stress and strain tensors conjugate after Hill (the strain tensors in these pairs belong to the family of generalized strain tensors). A method is presented for generating reduced forms of the constitutive relations with the aid of pairs of Lagrangian and Eulerian stress and strain tensors conjugate in the generalized sense which are obtained from pairs of Lagrangian tensors conjugate after Hill by mapping tensor fields on one configuration of a deformable body to tensor fields on another configuration.      

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.     

Some Basis-Free Formulae for the Time Rate and Conjugate Stress of Logarithmic Strain Tensor

In this paper, two kinds of tensor equations are studied and their solutions are derived in general cases. Then, some compact basis-free representations for the time rate and conjugate stress of logarithmic strain tensors are proposed using six different methods. In addition, relations between the coefficients in these expressions are disclosed. Subsequently, all these basis-free expressions given in this paper are validated for the cases of distinct stretches and double coalescence, respectively.     

Numerical Investigation of Large Elastoplastic Strains of Three-Dimensional Bodies

A method of stress—strain analysis of elastoplastic bodies with large displacements, rotations, and finite strains is developed. The incremental loading technique is used within the framework of the arbitrary Lagrangian—Eulerian formulation. Constitutive equations are derived which relate the Jaumann derivative of the Cauchy—Euler stress tensor and the strain rate. The spatial discretization is based on the FEM and multilinear three-dimensional isoparametric approximation. An algorithm of stress—strain analysis of elastic, hyperelastic, and perfectly plastic bodies is given. Numerical examples demonstrate the capabilities of the method and its software implementation __________ Translated from Prikladnaya Mekhanika, Vol. 41, No. 6, pp. 36–43, June 2005.     

A Natural Generalization of Hypoelasticity and Eulerian Rate Type Formulation of Hyperelasticity

According to the classical hypoelasticity theory, the hypoelasticity tensor, i.e. the fourth order Eulerian constitutive tensor, characterizing the linear relationship between the stretching and an objective stress rate, is dependent on the current stress and must be isotropic. Although the classical hypoelasticity in this sense includes as a particular case the isotropic elasticity, it fails to incorporate any given type of anisotropic elasticity. This implies that one can formulate the isotropic elasticity as an integrable-exactly classical hypoelastic relation, whereas one can in no way do the same for any given type of anisotropic elasticity. A generalization of classical theory is available, which assumes that the material time derivative of the rotated stress is dependent on the rotated Cauchy stress, the rotated stretching and a Lagrangean spin, linear and of the first degree in the latter two. As compared with the original idea of classical hypoelasticity, perhaps the just-mentioned generalization might be somewhat drastic. In this article, we show that, merely replacing the isotropy property of the aforementioned stress-dependent hypoelasticity tensor with the invariance property of the latter under an R-rotating material symmetry group R⋆ G ₀, one may establish a natural generalization of classical theory, which includes all of elasticity. Here R is the rotation tensor in the polar decomposition of the deformation gradient and G ₀ any given initial material symmetry group. In particular, the classical case is recovered whenever the material symmetry is assumed to be isotropic. With the new generalization it is demonstrated that any two non-integrable hypoelastic relations based on any two objective stress rates predict quite different path-dependent responses in nature and hence can in no sense be equivalent. Thus, the non-integrable hypoelastic relations based on any given objective stress rate constitute an independent constitutive class in its own right which is disjoint with and hence distinguishes itself from any class based on another objective stress rate. Only for elasticity, equivalent hypoelastic formulations based on different stress rates may be established. Moreover, universal integrability conditions are derived for all kinds of objective corotational stress rates and for all types of material symmetry. Explicit, simple, integrable-exactly hypoelastic relations based on the newly discovered logarithmic stress rate are presented to characterize hyperelasticity with any given type of material symmetry. It is shown that, to achieve the latter goal, the logarithmic stress rate is the only choice among all infinitely many objective corotational stress rates. This revised version was published online in August 2006 with corrections to the Cover Date.     

Some basis-free expressions for stresses conjugate to Hill’s strains through solving the tensor equation AX + XA = C

In the present paper, some new and compact basis-free expressions are derived for representing stresses conjugate to the Lagrangean Hill’s strain measures. Solving a tensor equation in the form of AX + XA = C, a rather simple general expression is obtained for the conjugate stresses in the case of distinct principal stretches. The result decomposes the conjugate stresses additively into two parts, such that one part is coaxial with the right stretch tensor while another part is orthogonal to it. An expression for the case of double-coalescent principal stretches is also obtained. As an example of the general results, simple expressions are obtained for the stress conjugate to the logarithmic strain.     

Cauchy’s stress theorem for stresses represented by measures

A version of Cauchy’s stress theorem is given in which the stress describing the system of forces in a continuous body is represented by a tensor valued measure with weak divergence a vector valued measure. The system of forces is formalized in the notion of an unbounded Cauchy flux generalizing the bounded Cauchy flux by Gurtin and Martins (Arch Ration Mech Anal 60:305–324, 1976). The main result of the paper says that unbounded Cauchy fluxes are in one-to-one correspondence with tensor valued measures with weak divergence a vector valued measure. Unavoidably, the force transmitted by a surface generally cannot be defined for all surfaces but only for almost every translation of the surface. Also conditions are given guaranteeing that the transmitted force is represented by a measure. These results are proved by using a new homotopy formula for tensor valued measure with weak divergence a vector valued measure.      

Conjugate stresses to two-point deformation tensors

In this paper some expressions for stresses conjugate to two-point deformation tensors are derived. These expressions are offered in both the component and basis-free forms. Although, the material time rate of a two-point deformation tensor is not an objective quantity, the stress tensor conjugate to it may be determined. The component-form expressions are obtained by using the notion of conjugacy together with the objectivity of the stress power. The component-form expressions are then extended to a unified basis-free form, using a theorem established for this purpose. The specific results are provided for all different cases of distinct and coalescent principal stretches in a three-dimensional inner product space.     

A finite strain kinematic hardening constitutive model based on Hencky strain: General framework, solution algorithm and application to shape memory alloys

The logarithmic or Hencky strain measure is a favored measure of strain due to its remarkable properties in large deformation problems. Compared with other strain measures, e.g., the commonly used Green-Lagrange measure, logarithmic strain is a more physical measure of strain. In this paper, we present a Hencky-based phenomenological finite strain kinematic hardening, non-associated constitutive model, developed within the framework of irreversible thermodynamics with internal variables. The derivation is based on the multiplicative decomposition of the deformation gradient into elastic and inelastic parts, and on the use of the isotropic property of the Helmholtz strain energy function. We also use the fact that the corotational rate of the Eulerian Hencky strain associated with the so-called logarithmic spin is equal to the strain rate tensor (symmetric part of the velocity gradient tensor). Satisfying the second law of thermodynamics in the Clausius-Duhem inequality form, we derive a thermodynamically-consistent constitutive model in a Lagrangian form. In comparison with the available finite strain models in which the unsymmetric Mandel stress appears in the equations, the proposed constitutive model includes only symmetric variables. Introducing a logarithmic mapping, we also present an appropriate form of the proposed constitutive equations in the time-discrete frame. We then apply the developed constitutive model to shape memory alloys and propose a well-defined, non-singular definition for model variables. In addition, we present a nucleation-completion condition in constructing the solution algorithm. We finally solve several boundary value problems to demonstrate the proposed model features as well as the numerical counterpart capabilities.     

Time Rates of Hill's Strain Tensors

Based on Hill's principal components formulae, two kinds of basis-free time rate formulae for Hill's strain tensors are proposed in this paper, which are obtained for the cases of distinct stretches, double coalescence and triple coalescence. Furthermore, relations between the coefficients in the representation for the strain tensor and its time derivative are disclosed. This revised version was published online in August 2006 with corrections to the Cover Date.     

Elastoplastic deformation of elements of an isotropic solid along paths of small curvature: Constitutive equations incorporating the stress mode

Constitutive equations relating the components of the stress tensor in a Eulerian coordinate system and the linear components of the finite-strain tensor are derived. These stress and strain measures are energy-consistent. It is assumed that the stress deviator is coaxial with the plastic-strain differential deviator and that the first invariants of the stress and strain tensors are in a nonlinear relationship. In the case of combined elastoplastic deformation of elements of the body, this relationship, as well as the relationship between the second invariants of the stress and strain deviators, is determined from fundamental tests on a tubular specimen subjected to proportional loading at several values of stress mode angle (the third invariant of the stress deviator). Methods to individualize these relationships are proposed. The initial assumptions are experimentally validated. The constitutive equations derived underlie an algorithm for solving boundary-value problems __________ Translated from Prikladnaya Mekhanika, Vol. 43, No. 6, pp. 43–55, June 2007.     

Basic Issues Concerning Finite Strain Measures and Isotropic Stress-Deformation Relations

It is indicated that the commonly-used Rivlin–Ericksen representation formula for isotropic tensor functions exhibits some properties that might be undesirable for its reasonable and effective applications. Towards clarification and improvement, a set of three mutually orthogonal tensor generators is introduced to achieve an alternative representation formula for isotropic symmetric tensor-valued functions of a symmetric tensor. This representation formula enables us to express the unknown representative coefficients in terms of simple, explicit tensorial inner products of the argument tensor and the value tensor without involving their eigenvalues. In particular, the tensorial interpolation expressions thus obtained assume a unified form for the three different cases of coalescence of the eigenvalues of the argument tensor. Moreover, each summand in the alternative representation formula is shown to inherit the continuity and differentiability properties of the represented isotropic tensor function. These results are used to study some basic issues concerning finite strain measures and stress-deformation relations of isotropic materials, such as continuity and differentiability properties of the representation, determination of the representative coefficients in terms of experimental data for stress and deformation tensors, and computations of finite strain measures. This revised version was published online in June 2006 with corrections to the Cover Date.     

A new two-point deformation tensor and its relation to the classical kinematical framework and the stress concept

Starting from the issue of what is the correct form for a Legendre transformation of the strain energy in terms of Eulerian and two-point tensor variables we introduce a new two-point deformation tensor, namely H=(F−F^−T)/2, as a possible deformation measure involving points in two distinct configurations. The Lie derivative of H is work conjugate to the first Piola–Kirchhoff stress tensor P. The deformation measure H leads to straightforward manipulations within a two-point setting such as the derivation of the virtual work equation and its linearization required for finite element implementation. The manipulations are analogous to those used for the Lagrangian and Eulerian frameworks. It is also shown that the Legendre transformation in terms of two-point tensors and spatial tensors require Lie derivatives. As an illustrative example we propose a simple Saint Venant–Kirchhoff type of a strain-energy function in terms of H. The constitutive model leads to physically meaningful results also for the large compressive strain domain, which is not the case for the classical Saint Venant–Kirchhoff material.     

Ellipticity conditions of the static equations of nonlinear elasticity

The relations of the nonlinear model of the theory of elasticity are considered. The Cauchy and the strain gradient tensors are taken to be the characteristics of the stress-strain state of a body. Sufficient conditions under which the static equations of elasticity are of elliptic type are established. These conditions are expressed in the form of constraints imposed on the derivatives of the elastic potential with respect to the strain-measure characteristics. The cases of anisotropic and isotropic bodies are treated, including the case where the Almansi tensor is taken to be the strain measure. The plane strain of a body is investigated using actual-state variables. Novosibirsk State University, Novosibirsk 630090. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 40, No. 2, pp. 196–203, March–April, 1999.     

Limiting Chain Extensibility Constitutive Models of Valanis–Landel Type

A new general constitutive model in terms of the principal stretches is proposed to reflect limiting chain extensibility resulting in severe strain-stiffening for incompressible, isotropic, homogeneous elastic materials. The strain-energy density involves the logarithm function and has the general Valanis–Landel form. For specific functions in the Valanis–Landel representation, we obtain particular strain-energies, some of which have been proposed in the recent literature. The stress–stretch response in some basic homogeneous deformations is described for these particular strain-energy densities. It is shown that the stress response in these deformations is similar to that predicted by the Gent model involving the first invariant of the Cauchy–Green tensor. The models discussed here depend on both the first and second invariants.