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.     

A simple mixture theory for ν Newtonian and generalized Newtonian constituents

This work presents the development of mathematical models based on conservation laws for a saturated mixture of ν homogeneous, isotropic, and incompressible constituents for isothermal flows. The constituents and the mixture are assumed to be Newtonian or generalized Newtonian fluids. Power law and Carreau–Yasuda models are considered for generalized Newtonian shear thinning fluids. The mathematical model is derived for a ν constituent mixture with volume fractions ${\phi_\alpha}$ using principles of continuum mechanics: conservation of mass, balance of momenta, first and second laws of thermodynamics, and principles of mixture theory yielding continuity equations, momentum equations, energy equation, and constitutive theories for mechanical pressures and deviatoric Cauchy stress tensors in terms of the dependent variables related to the constituents. It is shown that for Newtonian fluids with constant transport properties, the mathematical models for constituents are decoupled. In this case, one could use individual constituent models to obtain constituent deformation fields, and then use mixture theory to obtain the deformation field for the mixture. In the case of generalized Newtonian fluids, the dependence of viscosities on deformation field does not permit decoupling. Numerical studies are also presented to demonstrate this aspect. Using fully developed flow of Newtonian and generalized Newtonian fluids between parallel plates as a model problem, it is shown that partial pressures p _α of the constituents must be expressed in terms of the mixture pressure p. In this work, we propose ${p_\alpha=\phi_\alpha p}$ and ${\sum_\alpha^\nu p_\alpha = p}$ which implies ${\sum_\alpha^\nu \phi_\alpha = 1}$ which obviously holds. This rule for partial pressure is shown to be valid for a mixture of Newtonian and generalized Newtonian constituents yielding Newtonian and generalized Newtonian mixture. Modifications of the currently used constitutive theories for deviatoric Cauchy stress tensor are proposed. These modifications are demonstrated to be essential in order for the mixture theory for ν constituents to yield a valid mathematical model when the constituents are the same. Dimensionless form of the mathematical models is derived and used to present numerical studies for boundary value problems using finite element processes based on a residual functional, that is, least squares finite element processes in which local approximations are considered in ${H^{k,p}\left(\bar{\Omega}^e\right)}$ scalar product spaces. Fully developed flow between parallel plates and 1:2 asymmetric backward facing step is used as model problems for a mixture of two constituents.     

Finite strain––isotropic hyperelasticity

This paper presents a strain energy density for isotropic hyperelastic materials. The strain energy density is decomposed into a compressible and incompressible component. The incompressible component is the same as the generalized Mooney expression while the compressible component is shown to be a function of the volume invariant J only. The strain energy density proposed is used to investigate problems involving incompressible isotropic materials such as rubber under homogeneous strain, compressible isotropic materials under high hydrostatic pressure and volume change under uniaxial tension. Comparison with experimental data is good. The formulation is also used to derive a strain energy density expression for compressible isotropic neo-Hookean materials. The constitutive relationship for the second Piola–Kirchhoff stress tensor and its physical counterpart, involves the contravariant Almansi strain tensor. The stress stretch relationship comprises of a component associated with volume constrained distortion and a hydrostatic pressure which results in volumetric dilation. An important property of this constitutive relationship is that the hydrostatic pressure component of the stress vector which is associated with volumetric dilation will have no shear component on any surface in any configuration. This same property is not true for a neo-Hookean Green’s strain–second Piola–Kirchhoff stress tensor formulation.     

Three dimensional large deformation analysis of phase transformation in shape memory alloys

Shape memory alloys（SMAs）have been explored as smart materials and used as dampers,actuator elements,and smart sensors.An important character of SMAs is its ability to recover all of its large deformations in mechanical loading-unloading cycles without showing permanent deformation.This paper presents a stress-induced phenomenological constitutive equation for SMAs,which can be used to describe the superelastic hysteresis loops and phase transformation between Martensite and Austenite.The Martensite fraction of SMAs is assumed to be dependent on deviatoric stress tensor.Therefore,phase transformation of SMAs is volume preserving during the phase transformation.The model is implemented in large deformation finite element code and cast in the updated Lagrangian scheme.In order to use the Cauchy stress and the linear strain in constitutive laws,a frame indifferent stress objective rate has to be used.In this paper,the Jaumann stress rate is used.Results of the numerical experiments conducted in this study show that the superelastic hysteresis loops arising with the phase transformation can be effectively captured.     

On the partition of rate of deformation in crystal plasticity

Two different partitions of the rate of deformation tensor into its elastic and plastic parts are derived for elastic–plastic crystals in which crystallographic slip is the only cause of plastic deformation. One partition is associated with the Jaumann, and the other with convected rate of the Kirchhoff stress. Different expressions for the plastic part of the rate of deformation are obtained, and corresponding constitutive inequalities discussed. Relationship with the plastic part of the rate of the Lagrangian strain is also given.     

Controllable states of stress for incompressible elastic solids

It is shown that the equilibrium states of Cauchy stress which can exist, in the absence of body force, in every incompressible, homogeneous, isotropic, elastic solid whose deviatoric stress range allows them, must have uniform deviatoric stress invariants. There is at least one such non-uniform stress state. The related problem for incompressible non-Newtonian fluids is also discussed.     

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.     

Stress–strain measures used in the ANSYS package to solve elastoplastic problems at finite strains

It is shown that the Dienes (or the Green–McInnis–Naghdi) derivative is used in the ANSYS package rather than the Jaumann derivative as an objective derivative for the Cauchy stress tensor when solving elastoplastic problems, although the usage of the Jaumann derivative is stated in the ANSYS theory reference manual. In this manual it is also stated that, for these problems, the strain tensor is the Hencky logarithmic strain tensor; however, in reality, this strain tensor is the right nonholonomic strain tensor such that the left nonholonomic strain tensor associated with the right one is generated by the Dienes derivative.     

Irreversible thermodynamics of liquid crystal interfaces

A macroscopic theory for the dynamics of isothermal compressible interfaces between nematic liquid crystalline polymers and isotropic viscous fluids has been formulated using classical irreversible thermodynamics. The theory is based on the derivation of the interfacial rate of entropy production for ordered interfaces, that takes into account interfacial anisotropic viscous dissipation as well as interfacial anisotropic elastic storage. The symmetry breaking of the interface provides a natural decomposition of the forces and fluxes appearing in the entropy production, and singles out the symmetry properties and tensorial dimensionality of the forces and fluxes. Constitutive equations for the surface extra stress tensor and for surface molecular field are derived, and their use in interfacial balance equations for ordered interfaces is identified. It is found that the surface extra stress tensor is asymmetric, since the anisotropic viscoelasticity of the nematic phase is imprinted onto the surface. Consistency of the proposed surface extra stress tensor with the classical Boussinesq constitutive equation appropriate to Newtonian interfaces is demonstrated. The anisotropic viscoelastic nature of the interface between nematic polymers (NPs) and isotropic viscous fluids is demonstrated by deriving and characterizing the dynamic interfacial tension. The theory provides for the necessary theoretical tools needed to describe the interfacial dynamics of NP interfaces, such as capillary instabilities, Marangoni flows, wetting and spreading phenomena.     

弹性大变形问题中的应力状态描述、奇异性问题和余能原理

对弹性大变形理论中的3方面问题进行了综述.首先,对各种应变度量的共轭应力进行综述.大变形问题引起的应力状态描述的复杂性引起了许多学者的兴趣,对这个问题的研究也促进了大变形弹性理论的发展.在各种特定问题中,人们提出了不同的应力张量来描述应力状态,如Caucby应力张量、第一类和第一二类Piola-Kirchhoff应力张...     

Constitutive equation with fractional derivatives for the generalized UCM model

A fundamental problem on the constitutive equation with fractional derivatives for the generalized upper convected Maxwell model (UCM) is studied. The existing investigations on the constitutive equation are reviewed and their limitations or deficiencies are highlighted. By utilizing the convected coordinates approach, a mathematically rigorous constitutive equation with fractional derivatives for the generalized UCM model is proposed, which has an explicit expression for the stress tensor. This model can be reduced to the linear generalized Maxwell model with fractional derivatives, the UCM model and some other existing models. In addition, the rheological properties of this proposed model in the start-up of simple shear and elongation flows are investigated. It is shown that this generalized UCM model can describe the various stress evolution processes and the strain hardening effect of the viscoelastic fluids.     

On the equilibrium response of different modeling approaches for the porosity variation in dry granular matter

A selection of models for the variation in porosity in dry granular flows is investigated and compared on the basis of thermodynamic consistency to illustrate their performance and limitations in equilibrium situations. To this end, the thermodynamic analysis, based on the Müller–Liu entropy principle, is employed to deduce the ultimate constitutive equations at equilibrium. Results show that while all the models deliver appropriate equilibrium expressions of the Cauchy stress tensor for compressible grains, the model in which the variation in porosity is treated kinematically yields a spherical stress tensor for incompressible grains. Only the model in which the variation in porosity is modeled by a dynamic equation can give rise to a non-spherical stress tensor at equilibrium. The present study illuminates the validity and thermodynamic justification of the two modeling approaches for the porosity variation in dry granular matter.     

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.     

Numerical simulations of non-isothermal three-dimensional flows in an extruder by a finite-volume method

This study considers numerical applications of a finite-volume method to steady non-isothermal flows in geometries close to a single-screw extruder. Two geometrical configurations of the channel, with gap and zero gap, are investigated. The simulations concern incompressible fluids obeying different constitutive equations: Newtonian, generalized Newtonian with shear-thinning properties (Carreau–Yasuda law), and two viscoelastic differential models, the upper convected maxwell (UCM) and the Phan–Thien/Tanner (PTT). The temperature dependence is described by a Williams–Landel–Ferry (WLF) equation. For discretizing the equations and unknowns, we use a staggered grid with a QUICK scheme for the convective-type terms and solve the set of governing equations by a decoupled algorithm, stabilized by a pseudo-transient stress term and an elastic viscous stress splitting (EVSS) technique, in the viscoelastic case for the UCM model. The numerical results enable us to state the influence of temperature and rheological properties on the flow characteristics in the geometries investigated and underline the complex behaviour of the materials in such configurations.     

On co-rotational and other rate type constitutive equations

We derive expressions for the dilatational properties of suspensions of gas bubbles in incompressible fluids, using a cell model for the suspension. A cell, consisting of a gas bubble centered in a spherical shell of incompressible fluid, is subjected to a purely dilatational boundary motion and the resulting stress at the cell boundary is obtained. The same dilatational boundary motion is prescribed at the boundary of an “equivalent” cell composed of a one-phase, uniformly compressible fluid with unknown dilatational properties. By specifying that the stress at the boundary of the one-phase cell is equal to the stress at the boundary of the two-phase suspension cell, we obtain expressions for the unknown dilatational properties as a function of observable properties of the suspension. The dilatational viscosity of a suspension with a Newtonian continuous phase and the analogous properties for suspensions with non-Newtonian continuous phases are obtained as functions of the boundary motion, volume fraction of gas, and properties of the incompressible continuous phase. Results are presented for continuous phases which are Newtonian fluids, second-order fluids, and Goddard—Miller model fluids.     

A theory of thermo-visco-elastic-plastic materials: thermomechanical coupling in simple shear

The constitutive relations of a theory of thermo-visco-elastic-plastic continuum have been formulated in Lagrangian form. The Lagrangian strains, strain rates, temperature, temperature rate and temperature gradients are considered as the independent constitutive variables. Three internal state variables (plastic strain tensor, back strain tensor and a scalar hardening parameter) are also incorporated. The axioms of objectivity and equipresence are followed. The Clausius–Duhem inequality is taken as the second law of thermodynamics. Several special theories are deduced based on material symmetries and/or conventionally adopted assumptions. The applications to the formation of shear bands and dynamic crack propagation are discussed.     

Constitutive Relations for Rivlin‐Ericksen Fluids Based on Generalized Rational Approximation

. A well‐known constitutive expression for the stress in an incompressible non‐Newtonian fluid is provided by the representation of the extra stress as a function of the Rivlin‐Ericksen tensors . If this function is ordered in terms of the number of space plus time derivatives and appropriately scaled, one obtains . Truncation at first order yields the usual Newtonian viscous stress while truncation at second order provides the second‐order Rivlin‐Ericksen fluid. Many rheologists believe that in polymeric fluids. However, the requirement causes the rest state of the second‐order fluid to be unstable. This paper shows how the approximation of via generalized rational functions eliminates the instability paradox. (Accepted June 17, 1998)     

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.