Generalized Hamilton's Principle for Inelastic Bodies within Non-equilibrium Thermodynamics

Within the thermodynamic framework with internal variables, the classical Hamilton's principle for elastic bodies is extended to inelastic bodies composed of materials whose free energy densities are point functions of internal variables, or the so-termed Green-inelastic bodies, subject to finite deformation and non-conservative external forces. Yet this general result holds true even without the Green-inelasticity presumption under a more general interpretation of the infinitesimal internal rearrangement. Three special cases are discussed following the generalized form: (a) the Green-elastic bodies whose free energy can be identified with the strain energy; (b) the Green-inelastic bodies composed of materials compliant with the additive decomposition of strain; and (c) the Green-inelastic bodies undergoing isothermal relaxation processes where the thermodynamic forces conjugate to internal variables, or the so-termed internal forces prove to be potential forces. This paper can be viewed as an extension of Yang et al. [1].


Introduction
Hamilton's principle, in connection with its generalized forms, is one of the most effective approaches in studying various non-linear systems subject to dynamic actions.It usually takes a brief and neat expression with clear and definite physical meanings pertaining to work and energy.

OPEN ACCESS
The efforts to extend and apply Hamilton's principle to continuum cast back to one and half century ago, as indicated by Truesdell and Toupin [2] and Biot [3], and it continues to attract a lot of interest.More frequently, it is introduced as a starting point in the studies of elastic media or structures, e.g., Toupin [4] on elastic dielectrics, Mindlin [5] on linear elasticity with micro-structure, Batra [6] on thermo-elastic fluids and solids, Simo et al. [7] on the Hamiltonian structure of non-linear elasticity, Grinfeld and Norris [8] on fluid-permeable elastic continua and Ma et al. [9] on microstructuredependent Timoshenko beams.
Additionally, efforts have been made trying to establish generalized forms of Hamilton's principle for the inelastic continuum [1,[10][11][12][13][14].These studies share the common background of non-equilibrium thermodynamics where the vast phenomena of inelasticity are abstracted in a general way known as the internal variables.As remarked by Rice [15], one of the initiators of such concept, this theory is general enough to encompass relevant material behaviors while so simple that the essential content is kept clear.Stolz [10,11] proposed the functional approach in nonlinear dynamics where the notions of Lagrange equations and, for the case of conservative systems, Hamilton's principle, were established for inelastic continuum and shockwaves, as reviewed by Germain [12].But in these earlier studies, the internal variables were not treated as generalized coordinates like the displacements and temperature and some preconditions in thermodynamics were not sufficiently stressed.Recent representative studies include Sievers and Anthony [13] and Anthony [14], mainly on complex fluid and heat conduction and Yang et al. [1] on Green-inelastic solids.
This paper can be viewed as an extension of Yang et al. [1] in four aspects: (a) the most general form of the Hamilton's principle for Green-inelastic bodies is given and the classical form for those undergoing relaxation processes, as focused by Yang [1], is handled as a special case; (b) it is expounded that the general form holds true even for non-Green-inelasticity under a more general interpretation of the infinitesimal internal rearrangement; (c) Green-inelastic bodies compliant with the additive decomposition of strain is investigated as a special case; and (d) the entire problem is discussed under lager deformation situations.
The paper is arranged as follows: Sections 2 and 3 introduce the essential backgrounds of large deformation continuum mechanics and thermodynamics with internal variables, respectively.In Section 4, the generalized Hamilton's principle for inelastic bodies is derived from the principle of virtual work, and in Section 5 three significant special cases are investigated.

Prerequisites in Large Deformation Continuum Mechanics
In this section, we simply formulate the prerequisites of our primary concerns in large (or finite) deformation continuum mechanics following the monograph of Truesdell and Toupin [2], where the integrated background is available.

Work Conjugate Strain and Stress in Lagrange Formulism
Consider a body with volume 0 V and mass density 0 ρ in the reference configuration and volume V and mass density ρ in the displaced or current configuration.To describe finite deformation, let The inverse of ( ) ,t F X exists if the material points in the neighborhood of X are in one-to-one correspondence with their current positions: The Lagrange strain tensor, denoted E, objective and symmetric, measuring finite deformation from an arbitrary reference configuration, can be derived: where the superscript T denotes transposition and I is the second-rank unit tensor.
In the Eulerian description, the deformation power is defined as: where T is the real (or Cauchy) stress and D the deformation-rate tensor given by ( ) , in which L is an Eulerian tensor denoting velocity gradient, = ∂ ∂ L v x , with ( ) being the velocity field in Eulerian coordinates.Two Lagrangian expressions for d P can be shown: where T is known as the Piola stress: ( ) T F (6) and S the Kirchhoff stress: (7) with J denoting the Jacobian determinant, det J = F .In the sense of Equation ( 5), S is termed by Hill [16] the conjugate stress of the Lagrange strain E, or S and E are work conjugates.

The Principle of Virtual Work in Lagrange Formulism
The Lagrangian equations of motion is given by: where a and b denote respectively acceleration and the prescribed body force per unit mass., where p u and p t are the prescribed displacement and surface traction per unit area.)The virtual strain field įE can be shown due to the definition of E: ( ) Multiplying the left-hand side of Equation ( 8) by į i u , integrating over the reference volume 0 V and applying the divergence theorem, we obtain one Lagrangian form of the principle of virtual work: ( ) where p t is the prescribed surface traction per unit reference area with The principle of virtual work Equation (10), which almost epitomizes the integrated boundary value problem for large deformation bodies except for the constitutive relations, i.e., in our case the relation of S and E discussed in the following section, serves as the starting point to find the generalized Hamilton's principle for inelastic bodies.

Thermodynamics with Internal Variables
The internal-variable theory of inelastic continua described by Rice [15] belongs to the domain of thermodynamics of irreversible processes, or non-equilibrium thermodynamics.A collection of discrete scalar internal variables are introduced to represent the extent of micro-structural rearrangement within the material sample with macroscopically homogeneous stress or strain and temperature.Thus, entropy and temperature may be defined at a non-equilibrium state associated with a fictitious "constrained equilibrium state" at which the internal variables are fixed someway on their actual values at the non-equilibrium state.It follows that the existence of a free-energy density can be assumed from which the entropy density and the stress can be derived.Some detailed and overall results on internal variable and constrained equilibrium state have been obtained by Kestin and Rice [17], Germain, Nguyen and Suquet [18], Maugin and Muschik [19], Chaboche [20,21] and Yang et al. [22,23] and others.
Consider a material sample with mass density 0 ρ and volume 0 V measured in a reference state and at a reference temperature 0 ϑ .Its thermodynamic state can be fully described by the state variables ( ) ϑ S ξ .Our problem is discussed in the state space denoted ( ) where E is the Lagrange strain, ϑ temperature and ξ a set of scalar internal variables,

Green-Inelastic Material
Inelastic materials whose free energy densities are point functions of internal variables were termed Green-inelastic materials by Yang et al. [1], which constitute Green-inelastic bodies.For a green-inelastic material sample, the specific free energy density φ in ( ) Neighboring constrained equilibrium states corresponding to different sets of internal variables are related by the following Gibbs equation: (12) where S denotes the Kirchhoff stress: and η the specific entropy density: and įω the variational dissipation work per unit reference volume: where f denotes the set of thermodynamic forces conjugate to internal variables ξ , , which are termed internal forces as compared with external forces: Non-equilibrium thermodynamics extends the thermostatic relation (12) to dynamic processes with į replaced by d dt where t denotes time: ( ) Div 0 where Div denotes the divergence operator with respect to Lagrangian coordinates, and F h is the Lagrangian heat-flux vector and h the Eulerian heat-flux vector.This inequality must be obeyed in any process and at any state, particularly when the temperature gradient vanishes.Thus, the dissipation inequality, or Kelvin inequality should be imposed on the evolution of internal variables:

Non-Green-Inelastic Material
If the presumption of Green-inelastic material is violated in that φ is not a point function of ξ but instead depends on their path history, Equation (11) does not stand.Rice [15] took the transformed version of Equation ( 11), ( ) , where PIR symbolically denotes the pattern of internal rearrangement.However, the pattern of Equation ( 12), packed with Equations ( 13)- (15), remains unchanged under a more general interpretation of į α ξ , that at any given pattern of internal rearrangement of the material sample, a set of discrete scalar infinitesimals For the change in the pattern of internal rearrangement symbolically denoted įPIR resulting from the imposition of įξ , the internal forces f exists but cannot be simply determined by Equation ( 16).More detailed descriptions on this general case are provided by Hill and Rice [24] and Rice [25].
As shown in the next section, the derivation of the generalized Hamilton's principle is based on the unchanged Equations ( 12)-( 15), so the result holds true even without the presumption of Green-inelasticity.But the altered meanings and characteristics of ξ must be recognized: α ξ is no longer a state variable and the infinitesimal į α ξ does not necessarily represent an infinitesimal change in α ξ , although the above interpretation may include this as a special case, as indicated by Rice [15].

Generalized Hamilton's Principle for Green-Inelastic Bodies
To discuss Hamilton's principle, it is essential to clarify the system of generalized coordinates first.With the groundwork laid in Sections 2 and 3, a displaced configuration of an inelastic body can be described by the displacement field u and the internal variables ξ on each material point, which builds up the system of generalized coordinates.u and ξ respectively determine the external and internal configuration of the inelastic body and thus are term external and internal generalized coordinates.The idea of viewing the internal variables as generalized coordinates is inspired by Yang [26] who proposed Hamilton's principle of entropy production for creep and relaxation process of the Green-inelastic materials under macroscopically homogeneous and constant stress or strain and temperature.
Submitting Equation ( 12) into Equation (10) to eliminate : į S E , and then integrating over an isothermal process ( į 0 ϑ = ) with constant time duration 0 1 [ , ] t t leads to: where, Φ is the free energy: K is the kinetic energy: (22) and in Lagrangian description it can be shown that: įA − is the variational work of the non-conservative external forces: (24) įΩ is the variational work done by the internal forces, or the variational dissipation work: If, as is usually assumed, the external forces are conservative, Equation ( 20) can be rewritten as: where L is called the Lagrangian function: and A is the potential energy of the external forces.Equations ( 20) and ( 26) are the generalized Hamilton's principle for green-inelastic bodies subject to respectively non-conservative and conservative external forces.The results are brief and neat as we supposed the Hamilton's principle to be.It is emphasized again that these results hold true for non-Green inelastic bodies described in Section 3.2, if only the different connotation of įΩ being regarded.

Three Significant Special Cases
Now let's evaluate some special cases of certain significance, also to have a better understanding of the generalized form.External forces are assumed conservative in this section just for convenience.

Green-Elastic Bodies
For green-elastic bodies, the internal variables vanish and every local state is an equilibrium state.
The free energy density φ reduces to the Helmholtz free energy density ( ) where T is the absolute temperature.
The Hamilton's principle for Green-elastic bodies subject to conservative external forces in isothermal processes is: where Ψ is the Helmholtz free energy: To obtain the classical Hamilton's principle, the small (or infinitesimal) deformation situation is discussed.For small deformation elastic bodies (of which V and 0 V , ρ and 0 ρ are not distinguished) in an isothermal process with T at a fixed value, Ψ is usually identified with the strain-energy W: where w is the strain-energy volume density such that the stress-strain relation is Then the classical Hamilton's principle is obtained by replacing Ψ with W: It is clarified here, in passing, that in quite a few materials the presumption of linear elasticity, i.e., : = C σ ε , is taken in the demonstration of the classical Hamilton's principle, which is, as indicated by Equation (29), not a necessity.The classical Hamilton's principle stands for any non-linear elasticity if and only if the strain-energy exists.
Comparing Equation (31) with Equation ( 26), we can draw the conclusion that the classical Hamilton's principle for Green-elastic bodies can be extended to Green-inelastic bodies by replacing the strain energy with the specific free energy while considering the dissipation work done by the internal forces.

Green-Inelastic Bodies with Additive Decomposition of Strain
The additive decomposition of strain into elastic and inelastic parts is frequently assumed in studies of inelastic constitutive relations.It actually comes of a more general theory cored by the multiplicative decomposition in the domain of large deformation kinematics.The multiplicative theory, contributed by Green and Naghdi [27] and Casey [28] and others, unfortunately, has too large an inclusion and too many assumptions to be applied here, so we follow the way of Lubliner [29] and only consider small deformation cases.
The infinitesimal strain ε is assumed to be decomposable into elastic part depending on stress and temperature and inelastic part depending on internal variables only: In large deformation situations, a similar additive decomposition of strain with certain approximations, e i = + E E E , is proposed by Green and Naghdi [27] and conditioned by Casey [28], which is adopted in the works of Simo and Ortiz [30] and others, though, not in ours.The decomposition Equation (32), as shown by Lubliner [29], accords with the existence of a free-energy density Likewise, considering isothermal processes, comparing Equation (36) with its general form (12), we obtain the generalized Hamilton's principle for Green-inelastic bodies with additive decomposition of strain: where, the Lagrange function L is: and the transformed variational dissipation work įΩ is: Compared with the general form ( 26), the advantage of the current form (38) is obvious that only the elastic part of the free energy is directly involved, which in most cases lies shallower than the inelastic part.

Green-Inelastic Bodies Undergoing Isothermal Relaxation Processes
In isothermal relaxation processes, with the external state variables, i.e., the strain and temperature fixed at prescribed values, the state space Since the internal forces are potential forces, the dissipation work done by them through the constant time duration 0 1 [ , ] t t can be rewritten as: where P is the internal potential energy: Submitting Equation (44) back to Equation (26), we obtain the generalized Hamilton's principle for Green-inelastic bodies undergoing isothermal relaxation processes: Compared with Equation (31), it takes a similar form to the classical Hamilton's principle for Green-elastic bodies, just by replacing the strain energy with the summation of the corresponding free energy and internal potential energy, as indicated by Yang et al. [1].

Conclusions
Within the thermodynamic framework with internal variables, the generalized Hamilton's principle for Green-inelastic bodies is established, brief and neat as expected.Compared with the classical Hamilton's principle for Green-elastic bodies, the strain energy is replaced by the specific free energy and the dissipation work done by the internal forces must be involved.The general result holds true even without the Green-inelasticity presumption under a more general interpretation of the infinitesimal internal rearrangement.For Green-inelastic bodies with additive decomposition of strain, only the elastic part of the free energy is directly involved in the generalized Hamilton's principle.For Green-inelastic bodies undergoing isothermal relaxation processes, internal forces prove to be potential forces and the Hamilton's principle takes a classical form.

1 {
local rearrangements within the material sample.
(17) that the local Clausius-Duhem inequality becomes: possible infinitesimal variations in internal arrangement.
if and only if ψ can be decomposed as: equation to the elastic part e ψ leads to:

¦ξ
Ǽ ξ and the external configuration of the body changes only through rigid motions.The free energy density becomes ( ) and it is quite visible that a potential energy density, denoted p, exists in the relaxation subspace: serves as the potential energy of the internal forces of the material sample and is thus termed the internal potential energy density.Following the definition (42), the relation between p and