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.      

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.     

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.     

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.     

Thermomechanical material modelling based on a hybrid free energy density depending on pressure,isochoric deformation and temperature

In order to represent temperature-dependent mechanical material properties in a thermomechanical consistent manner it is common practice to start with the definition of a model for the specific Helmholtz free energy. Its canonical independent variables are the Green strain tensor and the temperature. But to represent calorimetric material properties under isobaric conditions, for example the exothermal behaviour of a curing process or the dependence of the specific heat on the temperature history, the temperature and the pressure should be taken as independent variables. Thus, in the field of calorimetry the Gibbs free energy is usually used as thermodynamic potential whereas in continuum mechanics the Helmholtz free energy is normally applied. In order to simplify the representation of calorimetric phenomena in continuum mechanics a hybrid free energy density is introduced. Its canonical independent variables are the isochoric Green strain tensor, the pressure and the temperature. It is related to the Helmholtz free energy density by a Legendre transformation. In combination with the additive split of the stress power into the sum of isochoric and volumetric terms this approach leads to thermomechanical consistent constitutive models for large deformations. The article closes with applications of this approach to finite thermoelasticity, curing adhesives and the glass transition.     

含随机参数的非线性黏弹性有限元计算

进行了粗粒土与结构接触面单调和循环加载试验，基于宏细观测量结果, 扩展了 损伤概念以 描述该类接触面在受载过程中的物态演化, 及由于物态演化导致的力学特性从初始状态到最终 稳定状态的连续变化过程. 揭示了接触面损伤的细观物理基础主要是接触面内土的颗粒破碎 和剪切压密这两种物态演化；指出接触面的剪胀体应变可以划分为可逆性和不可逆性剪胀体 应变两部分，其中不可逆性剪胀体应变可作为接触面损伤发展的宏观量度，因此其归一化 形式可作为一种损伤因子的定义；提出了建立粗粒土与结构接触面一种损伤本构关系的基本思路.     

Constitutive equations for hot-working of metals

Elevated temperature deformation processing - “hot-working,” is an important step during the manufacturing of most metal products. Central to any successful analysis of a hot-working process is the use of appropriate rate and temperature-dependent constitutive equations for large, interrupted inelastic deformations, which can faithfully account for strain-hardening, the restoration processes of recovery and recrystallization and strain rate and temperature history effects. In this paper we develop a set of phenomenological, internal variable type constitutive equations describing the elevated temperature deformation of metals. We use a scalar and a symmetric, traceless, second-order tensor as internal variables which, in an average sense, represent an isotropic and an anisotropic resistance to plastic flow offered by the internal state of the material. In this theory, we consider small elastic stretches but large plastic deformations (within the limits of texturing) of isotropic materials. Special cases (within the constitutive framework developed here) which should be suitable for analyzing hot-working processes are indicated.     

Stochastic damage model for concrete based on energy equivalent strain

Starting with the framework of conventional elastoplastic damage mechanics, a class of stochastic damage constitutive model is derived based on the concept of energy equivalent strain. The stochastic damage model derived from the parallel element model is adopted to develop the uniaxial damage evolution function. Based on the expressions of damage energy release rates (DERRs) conjugated to the damage variables thermodynamically, the concept and its tensor formulations of energy equivalent strain is proposed to bridge the gap between the uniaxial and the multiaxial constitutive models. Furthermore, a simplified coupling model is proposed to consider the evolution of plastic strain. And the analytical expressions of the constitutive model in 2-D are established from the abstract tensor expression. Several numerical simulations are presented against the biaxial loading test results of concrete, demonstrating that the proposed models can reflect the salient features for concrete under uniaxial and biaxial loading conditions.     

Nonlinear Eulerian thermoelasticity for anisotropic crystals

A complete continuum thermoelastic theory for large deformation of crystals of arbitrary symmetry is developed. The theory incorporates as a fundamental state variable in the thermodynamic potentials what is termed an Eulerian strain tensor (in material coordinates) constructed from the inverse of the deformation gradient. Thermodynamic identities and relationships among Eulerian and the usual Lagrangian material coefficients are derived, significantly extending previous literature that focused on materials with cubic or hexagonal symmetry and hydrostatic loading conditions. Analytical solutions for homogeneous deformations of ideal cubic crystals are studied over a prescribed range of elastic coefficients; stress states and intrinsic stability measures are compared. For realistic coefficients, Eulerian theory is shown to predict more physically realistic behavior than Lagrangian theory under large compression and shear. Analytical solutions for shock compression of anisotropic single crystals are derived for internal energy functions quartic in Lagrangian or Eulerian strain and linear in entropy; results are analyzed for quartz, sapphire, and diamond. When elastic constants of up to order four are included, both Lagrangian and Eulerian theories are capable of matching Hugoniot data. When only the second-order elastic constant is known, an alternative theory incorporating a mixed Eulerian–Lagrangian strain tensor provides a reasonable approximation of experimental data.     

The constitutive equations of finite strain poroelasticity in the light of a micro-macro approach

After recalling the constitutive equations of finite strain poroelasticity formulated at the macroscopic level, we adopt a microscopic point of view which consists of describing the fluid-saturated porous medium at a space scale on which the fluid and solid phases are geometrically distinct. The constitutive equations of poroelasticity are recovered from the analysis conducted on a representative elementary volume of porous material open to fluid mass exchange. The procedure relies upon the solution of a boundary value problem defined on the solid domain of the representative volume undergoing large elastic strains. The macroscopic potential, computed as the integral of the free energy density over the solid domain, is shown to depend on the macroscopic deformation gradient and the porous space volume as relevant variables. The corresponding stress-type variables obtained through the differentiation of this potential turn out to be the macroscopic Boussinesq stress tensor and the pore pressure. Furthermore, such a procedure makes it possible to establish the necessary and sufficient conditions to ensure the validity of an ‘effective stress’ formulation of the constitutive equations of finite strain poroelasticity. Such conditions are notably satisfied in the important case of an incompressible solid matrix.     

A 3D finite strain phenomenological constitutive model for shape memory alloys considering martensite reorientation

Most devices based on shape memory alloys experience both finite deformations and non-proportional loading conditions in engineering applications. This motivates the development of constitutive models considering finite strain as well as martensite variant reorientation. To this end, in the present article, based on the principles of continuum thermodynamics with internal variables, a three-dimensional finite strain phenomenological constitutive model is proposed taking its basis from the recent model in the small strain regime proposed by Panico and Brinson (J Mech Phys Solids 55:2491–2511, 2007). In the finite strain constitutive model derivation, a multiplicative decomposition of the deformation gradient into elastic and inelastic parts, together with an additive decomposition of the inelastic strain rate tensor into transformation and reorientation parts is adopted. Moreover, it is shown that, when linearized, the proposed model reduces exactly to the original small strain model.     

An energy-based damage model of geomaterials—I. Formulation and numerical results

In this paper, an anisotropic damage model is established in strain space to describe the behaviour of geomaterials under compression-dominated stress fields. The research work focuses on rate-independent and small-deformation behaviour during isothermal processes. It is emphasized that the damage variables should be defined microstructurally rather than phenomenologically for geomaterials, and a second-order fabric tensor is chosen as the damage variable. Starting from it, a one-parameter damage-dependent elasticity tensor is deduced based on tensorial algebra and thermodynamic requirements ; a fourth-order damage characteristic tensor, which determines anisotropic damaging, is deduced within the framework of Rice, 1971 normality structure in Part II of this paper. An equivalent state is developed to exclude the macroscopic stress⧹strain explicitly from the relevant constitutive equations. Finally, some numerical results are worked out to illustrate the mechanical behaviour of this model.     

Determination of the Constitutive Relation Parameters of a Metallic Material by Measurement of Temperature Increment in Compressive Dynamic Tests

The mechanical behaviour of a material can be established by an analytic expression called the constitutive relation that shows stress as a function of plastic strain, strain rate, temperature, and possibly other thermo-mechanical variables. The constitutive relation usually includes such parameters as coefficients or exponents that must be determined. At a high strain rate, the heat generated during the deformation process is directly related to the plastic deformation energy of the material. This energy can be calculated from the plastic work, resulting in an expression that includes the constitutive relation parameters as variables. The heat generated can also be estimated by measuring the temperature surface of the specimen during compressive tests using the technique of infrared thermography. The objective of this paper is to present a procedure for determining the constitutive relation parameters by measuring the temperature increase associated with plastic strain in compressive Hopkinson tests. The procedure was applied to estimate the parameters of the Johnson–Cook constitutive relation of an aluminium alloy (Al6082).     

A ferroelectric and ferroelastic 3D hexahedral curvilinear finite element

An isoparametric 3D electromechanical hexahedral finite element integrating a 3D phenomenological ferroelectric and ferroelastic constitutive law for domain switching effects is proposed. The model presents two internal variables which are the ferroelectric polarization (related to the electric field) and the ferroelastic strain (related to the mechanical stress). An implicit integration technique of the constitutive equations based on the return-mapping algorithm is developed. The mechanical strain tensor and the electric field vector are expressed in a curvilinear coordinate system in order to handle the transverse isotropy behavior of ferroelectric ceramics. The hexahedral finite element is implemented into the commercial finite element code Abaqus® via the subroutine user element. Some linear (piezoelectric) and non linear (ferroelectric and ferroelastic) benchmarks are considered as validation tests.     

A Generalized Theory of Stress and Strain Measures in the Classical Continuum Mechanics

A generalized theory of stress and strain tensor measures in the classical continuum mechanics is discussed: the main axioms of the theory are proposed, the general formulas for new tensor measures are derived, arid an energy conjugate theorem is formulated to distinguish the complete Lagrangian class of measures. As a subclass, a simple Lagrangian class of energy conjugate measures of stresses and finite strains is constructed in which the families of holonomic and corotational measures are distinguished. The characteristics of holonomic and corotational measures are studied by comparing the tensor measures of the simple Lagrangian class with one another and with logarithmic measures. For the simple Lagrangian class and its families, their completeness and closure are shown with respect to the choice of a generating pair of energetically conjugate measures. The applications of the new tensor measures in modeling the properties of plasticity, viscoelasticity, and shape memory are mentioned.     

General invariant representations of the constitutive equations for isotropic nonlinearly elastic materials

This paper develops general invariant representations of the constitutive equations for isotropic nonlinearly elastic materials. Different sets of mutually orthogonal unit tensor bases are constructed from the strain argument tensor by using the representation theorem and corresponding irreducible invariants are defined. Their relations and geometrical interpretations are established in three dimensional principal space. It is shown that the constitutive law linking the stress and strain tensors is revealed to be a simple relationship between two vectors in the principal space. Relative to two different sets of the basis tensors, the constitutive equations are transformed according to the transformation rule of vectors. When a potential function is assumed to exist, the vector associated with the stress tensor is expressed in terms of its gradient with respect to the vector associated with the strain tensor. The Hill’s stability condition is shown to be that the scalar product of the increment of those two vectors must be positive. When potential function exists, it becomes to be that the 3 × 3 constitutive matrix derived from its second order derivative with respect to the vector associated with the strain must be positive definite. By decomposing the second order symmetric tensor space into the direct sum of a coaxial tensor subspace and another one orthogonal to it, the closed form representations for the fourth order tangent operator and its inversion are derived in an extremely simple way.     

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.