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.      

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.     

A variational differential quadrature solution to finite deformation problems of hyperelastic shell-type structures: a two-point formulation in Cartesian coordinates

A new numerical approach is presented to compute the large deformations of shell-type structures made of the Saint Venant-Kirchhoff and Neo-Hookean materials based on the seven-parameter shell theory. A work conjugate pair of the first Piola Kirchhoff stress tensor and deformation gradient tensor is considered for the stress and strain measures in the paper. Through introducing the displacement vector, the deformation gradient, and the stress tensor in the Cartesian coordinate system and by means of the chain rule for taking derivative of tensors, the difficulties in using the curvilinear coordinate system are bypassed. The variational differential quadrature (VDQ) method as a pointwise numerical method is also used to discretize the weak form of the governing equations. Being locking-free, the simple implementation, computational efficiency, and fast convergence rate are the main features of the proposed numerical approach. Some well-known benchmark problems are solved to assess the approach. The results indicate that it is capable of addressing the large deformation problems of elastic and hyperelastic shell-type structures efficiently.     

On the Equivalence of Energetic and Geometric Shear Factors Based on Saint Venant Flexure

The equivalence of the shear compliance tensor obtained via an energy approach and the shear compliance tensor obtained via a geometric/kinematic approach is proven in the context of Saint Venant flexure. Specifically, it is shown that such an equivalence holds, as a general result, for sections of arbitrary geometry when both these tensors are computed on the basis of the long wavelength shear warpage, i.e., of the shear warpage independent from the longitudinal abscissa, provided all long-wavelength terms are included. The equivalence of energetic and kinematic shear factors stems as an immediate consequence of the equivalence of shear compliance tensors. A general analytical proof of this result is provided by analyzing in full detail the 3-dimensional elastostatic solution, inclusive of short-wavelength terminal fields, for a tip loaded cantilever. In particular, the developments reported exploit the objective tensor representations of stress and displacement fields solution to Saint Venant problem recently presented in a companion paper (Serpieri and Rosati, J. Elast., 2013). The well posedness of the concepts of energetic principal shear axes and geometric principal shear axes is shown as a further direct consequence. In particular, principal shear axes are shown to deserve the same legitimacy of principal bending axes both under an energetic and a geometric/kinematic viewpoint.     

Noncanonical Poisson brackets for elastic and micromorphic solids

This paper investigates the Lagrangian-to-Eulerian transformation approach to the construction of noncanonical Poisson brackets for the conservative part of elastic solids and micromorphic elastic solids. The Dirac delta function links Lagrangian canonical variables and Eulerian state variables, producing noncanonical Poisson brackets from the corresponding canonical brackets. Specifying the Hamiltonian functionals generates the evolution equations for these state variables from the Poisson brackets. Different elastic strain tensors, such as the Green deformation tensor, the Cauchy deformation tensor, and the higher-order deformation tensor, are appropriate state variables in Poisson bracket formalism since they are quantities composed of the deformation gradient. This paper also considers deformable directors to comprise the three elastic strain density measures for micromorphic solids. Furthermore, the technique of variable transformation is also discussed when a state variable is not conserved along with the motion of the body.     

Duality in constitutive formulation of finite-strain elastoplasticity based on F=FeFp and F=FF decompositions

A constitutive theory for large elastic–plastic deformations is presented by employing F=F^pF^e decomposition of the total deformation gradient. A duality in constitutive formulation based on this and the well-known Lee's decomposition F=F_eF_p is established for isotropic polycrystalline and single crystal plasticity.     

Derivatives and Rates of the Stretch and Rotation Tensors

General expressions for the derivatives and rates of the stretch and rotation tensors with respect to the deformation gradient are derived. They are both specialized to some of the formulas already available in the literature and used to derive some new ones, in three and two dimensions. Essential ingredients of the treatment are basic elements of differential calculus for tensor valued functions of tensors and recently derived results on the solution of the tensor equation A X + XA= H in the unknown X. This revised version was published online in August 2006 with corrections to the Cover Date.     

A weakened form of Ilyushin's postulate and the structure of self-consistent Eulerian finite elastoplasticity

From a general standpoint in terms of internal variables, we formulate a general theory of self-consistent Eulerian finite elastoplasticity based on the additive decomposition of the Eulerian strain rate, i.e., D=D^e+D^p, as well as two consistency criteria. In this theory, the elastic behaviour is characterized by an exactly integrable elastic rate equation for D^e with a general form of complementary elastic potential. It is assumed that the yield function depends in a general manner on the Kirchhoff stress and the internal variables. Moreover, the plastic rate equation for D^p and the evolution equation for each internal variable are allowed to assume general forms relying on the just-mentioned variables and the stress rate. It is indicated that two consistency criteria, i.e., the self-consistency for the elastic rate equation and Prager's yielding stationarity, lead to the unique choice of objective rates, i.e., the logarithmic rate.The structure of the above theory is further studied and examined by virtue of a weakened form of Ilyushin's postulate. In a spinning frame defining the logarithmic rate, we introduce the notion of standard elastoplastic strain cycle, which starts at a point not on but inside a yield surface and incorporates only one infinitesimal plastic subpath. We show that this type of strain cycle is always possible. Then, by ruling out strain cycles starting at points on yield surfaces we propose a weakened form of Ilyushin's postulate, which says that the changing rate of the stress work done along every standard strain cycle should be non-negative, whenever the incorporated plastic subpath tends to vanish. By virtue of simple, rigorous procedures, we demonstrate that this weakened form of Ilyushin's postulate is adequate to ensure direct results concerning the normality rule and the convexity of the yield surface in the context of the foregoing Eulerian finite elastoplasticity theory. Specifically, with an exactly integrable elastic rate equation defining D^e, we prove that, in the space of the Kirchhoff stresses, the difference (D−D^e) is just the gradient of the yield function multiplied by a plastic multiplier, and thus bears the very kinematical and physical feature of plastic strain rate. Furthermore, we prove that, in the space of the Kirchhoff stresses, the elastic domain bounded by each yield surface should be convex. The main results are derived in a self-contained manner within the context of an Eulerian theory of finite elastoplasticity, without involving issues concerning how to define intermediate stress-free states and plastic strains, etc.     

Basis free relations for the conjugate stresses of the strains based on the right stretch tensor

In this paper, general relations between two different stress tensors T^f and T^g, respectively conjugate to strain measure tensors f(U) and g(U) are found. The strain class f(U) is based on the right stretch tensor U which includes the Seth–Hill strain tensors. The method is based on the definition of energy conjugacy and Hill’s principal axis method. The relations are derived for the cases of distinct as well as coalescent principal stretches. As a special case, conjugate stresses of the Seth–Hill strain measures are then more investigated in their general form. The relations are first obtained in the principal axes of the tensor U. Then they are used to obtain basis free tensorial equations between different conjugate stresses. These basis free equations between two conjugate stresses are obtained through the comparison of the relations between their components in the principal axes, with a possible tensor expansion relation between the stresses with unknown coefficients, the unknown coefficients to be obtained. In this regard, some relations are also obtained for T⁽⁰⁾ which is the stress conjugate to the logarithmic strain tensor lnU.     

Strain Rates and Material Spins

The material time rate of Lagrangean strain measures, objective corotational rates of Eulerian strain measures and their defining spin tensors are investigated from a general point of view. First, a direct and rigorous method is used to derive a simple formula for the gradient of the tensor-valued function defining a general class of strain measures. By means of this formula and the chain rule as well as Sylvester's formula for eigenprojections, explicit basis-free expressions for the material time rate of an arbitrary Lagrangean strain measure can be derived in terms of the right Cauchy–Green tensor and the material time rate of any chosen Lagrangean strain measure, e.g. Hencky's logarithmic strain measure. These results provide a new derivation of Carlson–Hoger's general gradient formula for an arbitrary generalized strain measure and supply a new, rigorous proof for Carlson–Hoger's conjecture concerning the n-dimensional case. Next, by virtue of the aforementioned gradient formula, a general fact for objective corotational rates and their defining spin tensors is disclosed: Let Ω = ϒ ( B, D, W) be any spin tensor that is continuous with respect to B, where B, D and B are the left Cauchy–Green tensor, the stretching tensor and the vorticity tensor. Then the corotational rate of an Eulerian strain measure defined by Ω is objective iff Ω = W + υ ( B, D), where Υ is isotropic. By means of this fact and certain necessary or reasonable requirements, it is further found that a single antisymmetric function of two positive real variables can be introduced to characterize a general class of spin tensors defining objective corotational rates. A general basis- free expression for all such spin tensors and accordingly a general basis-free expression for a general class of objective corotational rates of an arbitrary Eulerian strain measure are established in terms of the left Cauchy–Green tensor B and the stretching tensor B as well as the introduced antisymmetric function. By choosing several particular forms of the latter, all commonly-known spin tensors and corresponding corotational rates are shown to be incorporated into these general formulas in a natural way. In particular, with the aid of these general formulae it is proved that an objective corotational rate of the Eulerian logarithmic strain measure ln V is identical with the stretching tensor D and moreover that in all possible strain tensor measures only ln V enjoys this property. This revised version was published online in August 2006 with corrections to the Cover Date.     

On a formulation for anisotropic elastoplasticity at finite strains invariant with respect to the intermediate configuration

The paper is concerned with a formulation of anisotropic finite strain inelasticity based on the multiplicative decomposition of the deformation gradient F=F_eF_p. A major feature of the theory is its invariance with respect to rotations superimposed on the inelastic part of the deformation gradient. The paper motivates and shows how such an invariance can be achieved. At the heart of the formulation is the mixed-variant transformation of the structural tensor, defined as the tensor product of the privileged directions of the material as given in a reference configuration, under the action of F_p. Issues related to the plastic material spin are discussed in detail. It is shown that, in contrast to the isotropic case, any flow function formulated purely in terms of stress quantities, necessarily exhibits a non-vanishing plastic material spin. The possible construction of spin-free rates is discussed as well, where it is shown that the flow rule must then depend not only on the stress but on the strain as well.     

Existence of a Solution for a Nonlinearly Elastic Plane Membrane Subject to Plane Forces

The two-dimensional nonlinear ‘membrane’ equations for a plate made of a Saint Venant–Kirchhoff material have been justified by D. Fox, A. Raoult and J.C. Simo (1993) by means of the method of formal asymptotic expansions applied to the three-dimensional equations of nonlinear elasticity. This model, which retains the material-frame indifference of the original three dimensional problem in the sense that its energy density is invariant under the rotations of R³, is equivalent to finding the critical points of a functional whose nonlinear part depends on the first fundamental form of the unknown deformed surface. We establish here a local existence result for these equations in the case of the membrane subject to forces parallel to its plane and we give qualitative properties of the solutions found in this fashion in terms of injectivity and of minimization. This revised version was published online in August 2006 with corrections to the Cover Date.     

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.     

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.     

n the fourth order tensor valued function of the stress in return map algorithm

针对各向同性材料,基于一组相互正交的基张量,建立了一套有 效的相关运算方法. 基张量中的两个分别是归一化的二阶单位张量和偏应力张量,另一个则 使用应力的各向同性二阶张量值函数经过归一化构造所得,三者共主轴. 根据张量函数表示 定理,本构方程和返回映射算法中所涉及到的应力的二阶、四阶张量值函数及其逆都由这组 基所表示. 推演结果表明：这些张量之间的运算,表现为对应系数矩阵之间的简单 关系. 其中,四阶张量求逆归结为对应的3\times3系数矩阵求逆,它对二阶张量的变换 则表现为该矩阵对3times 1列阵的变换. 最后,对这些变换关系应用于返回映 射算法的迭代格式进行了相关讨论.     

On micromechanical modeling of particulate composites with inclusions of various shapes

The effective elastic properties of statistically homogeneous two-phase particulate composites are considered. Several first-order micromechanical models are re-written in terms of the inclusion compliance contribution tensor (H-tensor). This tensor is a convenient tool to evaluate contribution of arbitrarily shaped inclusions and cavities to the overall composite properties.For any inclusion shape, the procedure starts with calculation of the H-tensor for a single inclusion. The non-interaction approximation is obtained by direct summation. More advanced micromechanical schemes are derived by substituting the non-interaction inclusion compliance contribution tensor into the formulae provided in the paper. The proposed procedure is illustrated by considering several two-dimensional and three-dimensional examples.     

A large deformation theory for rate-dependent elastic–plastic materials with combined isotropic and kinematic hardening

We have developed a large deformation viscoplasticity theory with combined isotropic and kinematic hardening based on the dual decompositions F=F^eF^p [Kröner, E., 1960. Allgemeine kontinuumstheorie der versetzungen und eigenspannungen. Archive for Rational Mechanics and Analysis 4, 273–334] and [Lion, A., 2000. Constitutive modelling in finite thermoviscoplasticity: a physical approach based on nonlinear rheological models. International Journal of Plasticity 16, 469–494]. The elastic distortion F^e contributes to a standard elastic free-energy ψ^(e), while , the energetic part of F^p, contributes to a defect energy ψ^(p) – these two additive contributions to the total free energy in turn lead to the standard Cauchy stress and a back-stress. Since F^e=FF^p-1 and , the evolution of the Cauchy stress and the back-stress in a deformation-driven problem is governed by evolution equations for F^p and – the two flow rules of the theory.We have also developed a simple, stable, semi-implicit time-integration procedure for the constitutive theory for implementation in displacement-based finite element programs. The procedure that we develop is “simple” in the sense that it only involves the solution of one non-linear equation, rather than a system of non-linear equations. We show that our time-integration procedure is stable for relatively large time steps, is first-order accurate, and is objective.