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.     

A unified residual-based thermodynamic framework for strain gradient theories of plasticity

A unified thermodynamic framework for gradient plasticity theories in small deformations is provided, which is able to accommodate (almost) all existing strain gradient plasticity theories. The concept of energy residual (the long range power density transferred to the generic particle from the surrounding material and locally spent to sustain some extra plastic power) plays a crucial role. An energy balance principle for the extra plastic power leads to a representation formula of the energy residual in terms of a long range stress, typically of the third order, a macroscopic counterpart of the micro-forces acting on the GNDs (Geometrically Necessary Dislocations). The insulation condition (implying that no long range energy interactions are allowed between the body and the exterior environment) is used to derive the higher order boundary conditions, as well as to ascertain a principle of the plastic power redistribution in which the energy residual plays the role of redistributor and guarantees that the actual plastic dissipation satisfies the second thermodynamics principle. The (nonlocal) Clausius-Duhem inequality, into which the long range stress enters aside the Cauchy stress, is used to derive the thermodynamic restrictions on the constitutive equations, which include the state equations and the dissipation inequality. Consistent with the latter inequality, the evolution laws are formulated for rate-independent models. These are shown to exhibit multiple size effects, namely (energetic) size effects on the hardening rate, as well as combined (dissipative) size effects on both the yield strength (intrinsic resistance to the onset of plastic strain) and the flow strength (resistance exhibited during plastic flow). A friction analogy is proposed as an aid for a better understanding of these two kinds of strengthening effects. The relevant boundary-value rate problem is addressed, for which a solution uniqueness theorem and a minimum variational principle are provided. Comparisons with other existing gradient theories are presented. The dissipation redistribution mechanism is illustrated by means of a simple shear model.     

Representation of the glass-transition in mechanical and thermal properties of glass-forming materials: A three-dimensional theory based on thermodynamics with internal state variables

In the vicinity of the glass transition, glass-forming materials exhibit pronounced frequency-dependent changes in the mechanical material properties, the thermal expansion behaviour and the specific heat. The frequency dependence becomes apparent under harmonic stress, strain or temperature excitations. The Prigogine-Defay ratio is a characteristic number which connects the changes in magnitude of these quantities at the glass transition. In order to represent the thermoviscoelastic properties of glass-forming materials in continuum mechanics, a three-dimensional approach which is based on the Gibbs free energy as thermodynamic potential is developed in this article. The Gibbs free energy depends on the stress tensor, the temperature and a set of internal variables which is introduced to take history-dependent phenomena into account. In the vicinity of an equilibrium reference state, the specific Gibbs free energy is approximated up to second order terms. Evaluating the Clausius-Duhem inequality, the constitutive relations for the strain tensor, the entropy and the internal variables are derived. In comparison with other approaches, the entropy, the strain tensor and the internal variables are functionals not only of the stress tensor but also of the temperature. Applying harmonic temperature- or stress-controlled excitations, complex frequency-dependent relations for the specific heat under constant stress, for the thermal expansion coefficients as well as for the dynamic mechanical compliance are obtained. The frequency-dependence of these quantities depicts the experimentally observed behaviour of glass-forming materials as published in literature. Under the assumption of isotropic material behaviour, it is shown that the developed theory is compatible with the Prigogine-Defay inequality for arbitrary values of the material parameters.     

On thermoelastic dielectrics

Constitutive equations for a linear thermoelastic dielectric are derived from the energy balance equation assuming dependence of the stored energy function on the strain tensor, the polarization vector, the polarization gradient tensor and entropy. A method is indicated for constructing a hierarchy of constitutive equations for materials with arbitrary symmetry by introducing various thermodynamic potentials. Maxwell's relations are constructed for the thermodynamic potential W^L. The entropy inequality is used to obtain stability conditions for an elastic dielectric in equilibrium under prescribed boundary constraints. Frequencies are explicitly determined for a plane wave propagating along the x₁-axis in an infinite centro-symmetric isotropic thermoelastic dielectric.     

On Entropy Flux of Transversely Isotropic Elastic Bodies

Recently (Liu in J. Elast. 90:259–270, 2008) thermodynamic theory of elastic (and viscoelastic) material bodies has been analyzed based on the general entropy inequality. It is proved that for isotropic elastic materials, the results are identical to the classical results based on the Clausius-Duhem inequality (Coleman and Noll in Arch. Ration. Mech. Anal. 13:167–178, 1963), for which one of the basic assumptions is that the entropy flux is defined as the heat flux divided by the absolute temperature. For anisotropic elastic materials in general, this classical entropy flux relation has not been proved in the new thermodynamic theory. In this note, as a supplement of the theory presented in (Liu in J. Elast. 90:259–270, 2008), it will be proved that the classical entropy flux relation need not be valid in general, by considering a transversely isotropic elastic material body.      

Evolution Equations for Non-Simple Viscoelastic Solids

Objectivity and compatibility with thermodynamics of evolution equations are examined in connection with the modelling of viscoelastic solids. The purpose of the paper is to show that the evolution equation for the stress is eventually obtained by means of a tensorial internal variable within the framework of the reference configuration. The non-simple character is realized by gradients of the internal variable. The thermodynamic analysis is developed by investigating the entropy inequality in the reference configuration and allowing for a non-zero extra-entropy flux. It follows that the evolution for the Cauchy stress tensor involves the Oldroyd derivative, irrespective of the form of the non-local terms.     

On continuum thermodynamics

Within the scope of classical continuum thermodynamics, we elaborate on the basic concepts and adopt a different approach from usual to the formulation of conservation laws and an entropy production inequality, both for a single phase continuum and for a mixture of any number of constituents. These conservation laws and the entropy inequality can be regarded as applicable to both local and nonlocal problems. In the case of a single phase continuum and for a simple material which is homogeneous in its reference configuration, under fairly mild smoothness assumptions, we prove that all the conservation laws reduce to the usual classical ones and the entropy production inequality reduces to the Clausius-Duhem inequality. Some attention is given to possible redundancies in the basic concepts, as well as to alternative forms of the energy equation and the entropy inequality. The latter is particularly significant in regard to different but equivalent formulations of mixture theory.     

Existence of entropy as a consequence of asymptotic stability

For a class of physical systems whose temporal evolution is governed by ordinary differential equations, the consequences of an assumption of asymptotic stability for equilibrium states in isolation remarkably resemble various forms of the second law of thermodynamics. Here we apply a known converse to Lyapunov's stability theorem to motivate both Gibbs' theory of thermostatics and the use of the Clausius-Duhem inequality for systems which are out of equilibrium and exchanging heat with their surroundings. We also discuss conditions under which the entropy of a system can be expressed as a sum of the entropies of its material points.     

考虑温度效应的弹性损伤一般理论

现有的各种损伤理论基本上都是关于等温问题的 ,且在不同程度上依赖于某些经验假设。本文在严格的不可逆热力学理论基础之上 ,建立了考虑温度效应的弹性损伤一般理论。推导出热弹性各向同性与各向异性损伤材料全部本构方程的一般形式 ,其中包括应力 应变关系、熵密度方程、损伤对偶张量表达式、热 固 损伤耦合的热传导方程和损伤演化方程。它们的特殊形式包含了等温各向同性与各向异性弹性损伤的本构方程     

On the thermodynamics of smooth muscle contraction

Cell function is based on many dynamically complex networks of interacting biochemical reactions. Enzymes may increase the rate of only those reactions that are thermodynamically consistent. In this paper we specifically treat the contraction of smooth muscle cells from the continuum thermodynamics point of view by considering them as an open system where matter passes through the cell membrane. We systematically set up a well-known four-state kinetic model for the cross-bridge interaction of actin and myosin in smooth muscle, where the transition between each state is driven by forward and reverse reactions. Chemical, mechanical and energy balance laws are provided in local forms, while energy balance is also formulated in the more convenient temperature form. We derive the local (non-negative) production of entropy from which we deduce the reduced entropy inequality and the constitutive equations for the first Piola–Kirchhoff stress tensor, the heat flux, the ion and molecular flux and the entropy. One example for smooth muscle contraction is analyzed in more detail in order to provide orientation within the established general thermodynamic framework. In particular the stress evolution, heat generation, muscle shorting rate and a condition for muscle cooling are derived.     

Tresca-type plastic materials in the theory of hypoelasticity II. Optical constitutive equations and birefringence in simple shear

Behavior of a Tresca type plastic dielectric is investigated theoretically from a continuum mechanical point of view. The optical constitutive equations are defined as special cases of a hypo-elastic dielectric of grade two. The singularity condition of the constitutive equations satisfies the Tresca yield criterion. The index deviator tensor is proportional to the stress deviator tensor and, then, the birefringence and the extinction angle are expressed by the stress deviator. Their numerical variations with the angle of shear in simple shear deformation are shown.     

Finite Third-order Gradient Elasticity and Thermoelasticity

A constitutive format for the third-order gradient elasticity is suggested. It includes both isotropic and anisotropic non-linear behavior under finite deformations. Appropriate invariant stress and strain variables are introduced, which allow for reduced forms of the elastic energy law that identically fulfill the objectivity requirement. After working out the transformation behavior under a change of the reference placement, the symmetry transformations for third-order materials can be introduced. After the mechanical third-order theory, an extension to thermoelasticity is given, and necessary and sufficient conditions are derived from the Clausius-Duhem inequality.     

Extended thermodynamics of molecular ideal gases

The objective of extended thermodynamics of molecular ideal gases is the determination of the 17 fields ofmass density, velocity, energy density, pressure deviator, heat flux, intrinsic energy density and intrinsic heat flux. The intrinsic energy represents the rotational or the vibrational energy of the molecules. The necessary field equations are based upon balance laws and the system of equations is closed by constitutive relations which are characteristic for the gas under consideration. The generality of the constitutive relations is restricted by theprinciple of material frame indifference, and by the entropy principle. These principles reduce the constitutive coefficients of all fluxes to the thermal and caloric equation of state of the gas and provide inequalities for the transport coefficients. The transport coefficients can be related to the shear viscosity, the heat conductivity, and the coefficients of self-diffusion and attenuation of sound waves, so that the field equations become quite specific. The theory is in perfect agreement with the kinetic theory of molecular gases. It is shown that in non-equilibrium the temperature is discontinuous at thermometric walls. The dynamic pressure and the volume viscosity, are discussed and it is shown how extended thermodynamics and ordinary thermodynamics are related.     

Original ArticleOn entropy flux-heat flux relation in thermodynamics with Lagrange multipliers

A typical problem often encountered in thermodynamics with Lagrange multipliers is to deduce the relation between the entropy flux and the heat flux from conditions imposed by the entropy principle employing the general entropy inequality. This problem is usually trivial if linear constitutive relations are assumed, and otherwise it can be quite difficult even though such a relation is expected. In this paper, we shall try to formulate such a problem in general and prove as theorems. Some previous thermodynamic theories with Lagrange multipliers are mentioned as typical examples. Received December 13, 1995     

Anisotropic plastic stress fields at a slowly propagating crack tip

Under the condition that any perfectly plastic stress components at a crack tip are nothing but the functions of 0 only making use of equilibrium equations. Hill anisotropic yield condition and unloading stress-strain relations, in this paper, we derive the general analytical expressions of anisotropic plastic stress fields at the slowly steady propagating tips of plane and anti-plane strain. Applying these general analytical expressions to the concrete cracks, the analytical expressions of anisotropic plastic stress fields at the-slowly steady propagating tips of Mode I and Mode III cracks are obtained. For the isotropic plastic material, the anisotropic plastic stress fields at a slowly propagating crack tip become the perfectly plastic stress fields.     

On a finite strain viscoplastic theory based on a new internal dissipation inequality

This work is focused on the theoretical development and numerical implementation of a viscoplastic law. According to the second law of thermodynamics a dissipation inequality described in the rotated material coordinate system is developed. Based on this dissipation inequality and the principle of maximum dissipation a finite strain viscoplastic model described also in the rotated material coordinate system is formulated. The evolution equations are expressed in terms of the material time derivatives of the rotated elastic logarithmic strain, the accumulated plastic strain and the strain-like tensor conjugate to the rotated back stress. The mathematical structure of this theory is concise and similar to that of the infinitesimal viscoplastic theory. These characteristics make the numerical implementation of this theory easy. The stress integration algorithm and the algorithmic tangent moduli for the infinitesimal theory can be applied to the numerical implementation of the present finite strain theory with a little reformulation. The complicated algorithmic formulations for most of other finite plastic laws can be therefore circumvented. In order to check the effectivity of the present finite strain theory a set of numerical examples under strict deformation conditions are presented. These numerical examples prove the excellent performance of the present viscoplastic material law at describing the finite strain elastoplastic and viscoplastic problems.     

散粒体卸载特性的三轴伸长试验研究

根据岩土体的加卸载变形过程和应力状态,卸载曲线分为弹性卸载段和卸载伸长段.采用三种不同粒径的玻璃珠和标准丰浦砂散粒材料进行了三轴伸长卸载试验,对比分析了不同材料、初始偏应力和不同围压下散粒体的卸载特性.研究表明散粒体的单向伸长卸载强度主要取决于与卸载方向垂直的围压和颗粒的内摩擦,且与围压成线性关系,初始偏应力状态及各向异性对于散粒体的伸长强度无影响,但影响其破坏发展形态.此外,材料特性包括表面粗糙度和级配、孔隙率等对内摩擦及伸长卸载特性均有影响.