On constitutive equations for fiber-reinforced nonlinearly viscoelastic solids

In this paper, we are interested in developing constitutive equations for fiber-reinforced nonlinearly viscoelastic bodies, in particular for transversely isotropic nonlinearly viscoelastic solids. It follows from results in the theory of algebraic invariants that constitutive equations for such materials can be expressed in terms of functions of 18 independent invariants associated with deformation and fiber orientation. These invariants are analyzed, and we obtain restrictions such as positivity of some of them.     

On thermodynamically consistent constitutive equations for fiber-reinforced nonlinearly viscoelastic solids with application to biomechanics

In this paper, we are interested in developing thermodynamically consistent constitutive equations for fiber-reinforced nonlinearly viscoelastic bodies, in particular for transversely isotropic nonlinearly viscoelastic solids, in isothermal processes. It follows from results in the theory of algebraic invariants that constitutive equations for such materials can be expressed in terms of functions of 18 independent invariants associated with deformation and fiber orientation: 10 of them are isotropic invariants and 8 of them are associated with the deformation and the orientation of the fiber. Among the 8 anisotropic invariants just 6 are related to the viscoelastic response. The terms in the Cauchy stress tensor associated to these 6 invariants are analyzed with respect to thermodynamical consistency, and we obtain restrictions for the corresponding constitutive coefficients. This framework is applied to viscoelastic potentials within the context of biomaterials.     

Constitutive structure in coupled non-linear electro-elasticity: Invariant descriptions and constitutive restrictions

A constitutive framework for electro-sensitive materials in the context of non-linear elasticity is analyzed. Constitutive equations are given in terms of energy functions that depend on several invariants. The study includes both the analysis of the invariants, which are present in the energy functions, and the analysis of constitutive restrictions that have to be obeyed by the constitutive functions. Isotropic as well as non-isotropic electro-sensitive elastomers are studied. The set of invariants that describe each material model is analyzed under two homogeneous deformations: (i) an uniaxial elongation and (ii) a simple shear deformation. These deformations are chosen since they are relevant to specific experiments, from which one may try to fit constitutive equations. The constitutive restrictions developed are based on classical ones used for isotropic non-linear elastic materials, in particular, are based on the Baker–Ericksen inequality and the ellipticity condition.     

各向同性率无关材料本构关系的不变性表示

在内变量理论的框架下,针对各向同性率无关材料,使用张量函数表示理论建立了塑性应变全量及增量本构关系的最一般的张量不变性表示. 它们均由3个完备不可约的基张量组合构成,这3个基张量分别是应力的零次幂、一次幂和二次幂. 因此得出,塑性应变、塑性应变增量与应力三者共主轴. 通过对基张量的正交化,给出了本构关系式在主应力空间中的几何解释. 进一步,全量(或增量)本构关系中3个组合因子被表达为应力、塑性应变(或塑性应变增量)的不变量的函数. 当塑性应变(或塑性应变增量)的3个不变量之间满足一定关系时,所给出的本构关系将退化为经典的形变理论(或塑性势理论).最后,还讨论它与奇异屈服面理论的关系,当满足一定条件时,两者是一致的.      

On the Comparison of Two Strategies to Formulate Orthotropic Hyperelasticity

The main goal of this work is to clarify the relation between two strategies to formulate constitutive equations for orthotropic materials at large strains. On the one hand, the classical approach is based on the incorporation of structural tensors into the free energy function via an enriched set of invariants. On the other hand, a fictitious isotropic configuration is introduced which renders an anisotropic, undeformed reference configuration via an appropriate linear tangent map. This formulation results in a reduced (with respect to the more general setting based on structural tensors) but nevertheless physically motivated set of invariants which are related to the invariants defined by structural tensors. As a main conceptual advantage standard isotropic constitutive equations can be applied and moreover, due to the reduced set of physically motivated invariants, the numerical treatment within a finite element setting becomes manageable.     

Nonlinear Orthotropic Elasticity: Only Six Invariants are Independent

It is often declared in the literature that the seven classical invariants used to characterize the strain energy of a compressible orthotropic elastic solid are independent. In this paper, we show that only six of the seven classical invariants are independent, and a syzygy exists between the classical invariants. Consequently, all other sets of seven invariants, proposed in the literature, that are uniquely related to the set of classical invariants, have only six independent invariants.     

The fundamental role of nonlocal and local balance laws of material forces in finite elastoplasticity and damage mechanics

In this paper the fundamental role of independent balance laws of material forces acting on dislocations and microdefects is shown. They enable a thermodynamically consistent formulation of dissipative deformation processes of continua with dislocation motion and defect evolution in the material space on meso- and microlevel.The balance laws of material forces together with the classical balance laws of physical forces and couples, first and second laws of thermodynamics for physical and material space and general constitutive equations are the basis to develop a thermodynamically consistent framework of nonlocal finite elastoplasticity and brittle and ductile damage.It is shown that a weakly-nonlocal formulation of the balance laws of material forces leads to gradient theories, where local theories are obtained, if all gradient contributions are assumed to be small. In this case the local balance laws of material forces together with the constitutive equations represent evolution laws of the material forces. In the classical approach of internal variables they are assumed from the outset with the result that there is a large number of different propositions in the literature.The well-known splitting test of a circular cylinder of concrete is simulated numerically, where the process of deformation in the physical space and defect and plastic evolution in the material space is represented.     

Anisotropic Mullins stress softening of a deformed silicone holey plate

Rubber like materials parts are designed using finite element code in which more and more precise and robust constitutive equations are implemented. In general, constitutive equations developed in literature to represent the anisotropy induced by the Mullins effect present analytical forms that are not adapted to finite element implementation. The present paper deals with the development of a constitutive equation that represents the anisotropy of the Mullins effect using only strain invariants. The efficiency of the modeling is first compared to classical homogeneous experimental tests on a filled silicone rubber. Second, the model is tested on a complex structure. In this aim, a silicone holey plate is molded and tested in tension, its local strain fields are evaluated by means of digital image correlation. The experimental results are compared to the simulations from the constitutive equation implemented in a finite element code. Global measurements (i.e. force and displacement) and local strain fields are successfully compared to experimental measurements to validate the model.     

Towards the solution of finite plane-strain problems for compressible elastic solids

A new approach to the solution of finite plane-strain problems for compressible Isotropie elastic solids is considered. The general problem is formulated in terms of a pair of deformation invariants different from those normally used, enabling the components of (nominal) stress to be expressed in terms of four functions, two of which are rotations associated with the deformation. Moreover, the inverse constitutive law can be written in a simple form involving the same two rotations, and this allows the problem to be formulated in a dual fashion.For particular choices of strain-energy function of the elastic material solutions are found in which the governing differential equations partially decouple, and the theory is then illustrated by simple examples. It is also shown how this part of the analysis is related to the work of F. John on harmonic materials.Detailed consideration is given to the problem of a circular cylindrical annulus whose inner surface is fixed and whose outer surface is subjected to a circular shear stress. We note, in particular, that material circles concentric with the annulus and near its surface decrease in radius whatever the form of constitutive law within the given class. Whether the volume of the material constituting the annulus increases or decreases depends on the form of law and the magnitude of the applied shear stress.     

The Debye Theory of Rotary Diffusion: History, Derivation, and Generalizations

Motivated by the Debye theory of rotary diffusion in a dipolar fluid, we systematically develop a continuum mechanical theory of rotary diffusion. This theory generalizes classical kinematics to include continuous rotary degrees of freedom and introduces an additional balance law associated with the rotary degrees of freedom. Various constitutive relations are proposed in accordance with standard procedures of nonlinear continuum mechanics. The resulting set of equations provides a properly invariant and thermodynamically consistent theory that allows for constitutive nonlinearities. In particular, the classical Debye theory along with the Nernst-Einstein relations are shown to follow from a special case of linear constitutive relations and an assumption of ideality in which the free energy consists only of a classical entropic contribution. Within our theory, the notion of osmotic pressure arises naturally as a consequence of accounting for forces that act conjugate to the rotary degrees of freedom and serves as the driving force for rotary diffusion. Accepted November 18, 2000?Published online April 23, 2001     

Study on the generalized prandtl-reuss constitutive equation and the corotational rates of stress tensor

I.IntroductionTilepl'ogl'ess11as.toifcertainextent,beenmadeintheelastic-plasticconstitutivetheoryatII[litedefbrlllations.Coil'paredwitllotherconstitutiverelations,thegeneralizedPrandtlReuss(P-R)equatiollsareextensivelystudiedandwidelyapplied.IndevelopingthegeneralizedP-Requation.itisusuallyassumedthatthedeformationrate(thesymmetricpartorvelocitygradiellt)isdecolllposedintotheelasticpartandplasticpart.TheplasticLIcf\'l.llliltlollrittcobeystilenormalfi(,xvrilleasillthecaseofinfinitcsilllnld…     

Invariants of a Remarkable Family of Nonlinear Equations

In classical literature, invariants of families of differentialequations were considered for linear equations only, e.g. the renownedLaplace invariants for linear hyperbolic partial differential equationsand invariants of linear ordinary differential equations with variablecoefficients. The restriction to linear equations was essential inpioneering works of Cockle, Laguerre, Halphen, andForsyth for tackling the problem of invariants of differentialequations. Lie regretted that these authors did not use advantagesprovided by his theory of infinite continuous groups, but he himself didnot undertake further developments in this direction.Recently, the present author considered the possibility hinted byLie's remark and introduced the infinitesimal technique in thetheory of invariants of families of differential equations thatwas lacking in old methods. In consequence, a simple unifiedapproach was developed for calculation of invariants of algebraicand differential equations independent on the assumption oflinearity of equations. It was employed recently for calculationof Laplace type invariants for parabolic equations. Here, themethod is applied to calculation of invariants for the family ofnonlinear equations appearing in the problem on linearization ofnonlinear ordinary differential equations.     

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.     

Approximate description of the inelastic deformation of an isotropic material with allowance for the stress mode

The elastoplastic deformation of an isotropic material is described using constitutive equations and allowing for the stress mode. The equations include two nonlinear functions that relate the first and second invariants of the stress and linear-strain tensors to the stress mode angle. It is proposed to use a linear rather than nonlinear relationship between the first invariants of the tensors. This simplification is validated by comparing calculated and experimental strains under loading with constant and variable stress mode angle     

Constitutive equations for describing the elastoplastic deformation of elements of a body along small-curvature paths in view of the stress mode

Equations relating the components of the stress and strain tensors (constitutive equations) are formulated in terms of Euler coordinates. The equations describe the finite elastoplastic deformation of an isotropic body along paths of small curvature. It is assumed that the stress deviator is coaxial with the plastic-strain differential deviator. The relationships between the first and second invariants of the stress and strain tensors in the case of complex elastoplastic deformation of the body’s elements are determined from base tests on tubular specimens loaded along rectilinear paths for several values of the stress mode angle. Methods for specification of these relationships are proposed. The assumptions adopted to derive the constitutive equations are validated experimentally __________ Translated from Prikladnaya Mekhanika, Vol. 42, No. 4, pp. 62–72, April 2006.     

Constitutive equations of elastoplastic isotropic materials that allow for the stress mode

The constitutive equations describing the elastoplastic deformation of an isotropic material and taking into account the stress mode are validated against available experimental data. We propose a method for the approximate determination of the base functions appearing in the constitutive equations and relating the first and second invariants of the stress tensors and the linear components of finite strains. The strain components obtained by this method are compared to experimental data     

Effects of stress invariants and reverse loading on ductile fracture initiation

Recent research studies on ductile fracture of metals have shown that the ductile fracture initiation is significantly affected by the stress state. In this study, the effects of the stress invariants as well as the effect of the reverse loading on ductile fracture are considered. To estimate the reduction of load carrying capacity and ductile fracture initiation, a scalar damage expression is proposed. This scalar damage is a function of the accumulated plastic strain, the first stress invariant and the Lode angle. To incorporate the effect of the reverse loading, the accumulated plastic strain is divided into the tension and compression components and each component has a different weight coefficient. For evaluating the plastic deformation until fracture initiation, the proposed damage function is coupled with the cyclic plasticity model which is affected by all of the stress invariants and pervious plastic deformation history.For verification and evaluation of this damage-plasticity constitutive equation a series of experimental tests are conducted on high-strength steel, DIN 1.6959. In addition finite element simulations are carried out including the integration of the constitutive equations using the modified return mapping algorithm. The modeling results show good agreement with experimental results.     

Initiation of localized shear bands in plane strain

This study is concerned with the initiation of localized shear bands in plane-strain tension and compression. A theoretical framework which views the initiation of such shear bands as a bifurcation phenomenon from a homogeneous equilibrium field in an elastic-plastic body is first briefly reviewed, and then the predictions of the theory are compared with some experimental observations on an aged maraging steel. The experiments support the physical relevance of the theory within the framework of continuum mechanics.Specifically, the comparison between theory and experiment is concerned with the critical strains to localization and the orientation of the shear bands relative to the load axis. The theoretical predictions are only in qualitative agreement with the experimental observations. A better agreement is obtained with use of the constitutive equations corresponding to the classical deformation theory (a simple vertex model) than with use of the constitutive equations of the classical flow theory. It is concluded that better constitutive equations for elastic-plastic materials are needed before theoretical predictions can be obtained which might be expected to be in closer quantitative accord with experiment.     

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.     

Modelling of stress softening in elastomeric materials: foundations of simple theories

This work is concerned with formulation of constitutive relations for materials exhibiting the stress softening phenomenon (known as the Mullins effect) typical observed in elastomeric and other amorphous materials during loading–reloading cycles. It is assumed that microstructural changes in such materials during the deformation process can be represented by a single scalar-valued softening variable whose evolution is accompanied by microforces satisfying their own law of balance, besides the classical laws of mechanics underlying macroscopic deformation of a material. The constitutive equations are then derived in consistency with thermodynamics of irreversible processes with the restriction to purely mechanical theory. The general form of the derived constitutive equations is subsequently simplified through introduction of additional assumptions leading to various models of the stress softening phenomenon. As an illustration of the general theory, it is shown that the so-called pseudo-elastic model proposed in the literature may be derived without an ad hoc postulate of the variational principle.