Journal of Engineering

Volume 2013 (2013), Article ID 516462, 8 pages

http://dx.doi.org/10.1155/2013/516462

Research Article

Some Variational Principles for Coupled Thermoelasticity

Francesco Marotti de Sciarra

Dipartimento di Ingegneria Strutturale, Università di Napoli Federico II, Via Claudio 21, 80125 Napoli, Italy

Received 22 September 2012; Accepted 18 December 2012

Academic Editor: Oronzio Manca

Copyright © 2013 Francesco Marotti de Sciarra. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The nonlinear thermoelasticity of type II proposed by Green and Naghdi is considered. The thermoelastic structural model is formulated in a quasistatic range, and the related thermoelastic variational formulation in the complete set of state variables is recovered. Hence a consistent framework to derive all the variational formulations with different combinations of the state variables is provided, and a family of mixed variational formulations, with different combinations of state variables, is provided starting from the general variational formulation. A uniqueness condition is provided on the basis of a suitable variational formulation.

1. Introduction

Coupled thermomechanical problem arises in a variety of important fields of application, including casting, metal forming, machining and other manufacturing processes, structural models, and others.

Green and Naghdi (GN) introduced a theory in which heat propagates as thermal waves at finite speed and does not necessarily involve energy dissipation. Another property of the GN theory of type II is the fact that the entropy flux vector is determined by means of the same potential as the mechanical stress tensor. Motivated by the procedure presented in [1, 2], this paper is concerned with the formulation of variational principles characterizing the solutions of the coupled thermomechanical problem for the GN model without dissipation, that is, type II. The variational characterization of the thermoelastic problem means the identification of a functional whose stationary points are solutions of the problem. Once this functional is known, the solutions of the problem can be identified with certain extrema of the functional.

Following the pioneering work of Biot [3], the variational forms of the coupled thermoelastic and thermoviscoelastic problems for the classical Fourier model have been investigated in several papers; see, for example, [2, 4–9] for the case of thermomechanical coupling in dissipative materials. Moreover consistent variational principles for structural problems concerning elastic and elastoplastic models in isothermal conditions are well developed; see, for example, [10–12].

On the contrary variational formulations for the GN theory of type II are addressed in a few papers; see, for example, [1, 13, 14].

In this paper the GN thermoelastic coupled structural model without dissipation is formulated in a suitable form so that we can provide the methodology to build the complete family of all the admissible variational formulations associated with the considered GN model. It is well-known that the GN theory of type II has been developed as a dynamic theory, but an important prerequisite for its use is a thorough understanding of the corresponding problems in which the dynamical effects are disregarded; see, for example, [15].

An advantage of such an approach consists in the fact that the boundary-value problem for the thermoelastic model is formulated in such a way that the model can be cast in terms of a multivalued structural operator defined in terms of all the state variables. This operator encompasses in a unique expression the field equation, the constitutive relations, the constraint relations and the initial conditions which describe the considered thermoelastic structural model.

The related non-smooth potential can then be evaluated by a direct integration along a ray in the operator domain and depends on all the state variables involved in the model. Appraising the generalized gradient of the non-smooth potential [10] and imposing its stationarity, the operator formulation of the problem is recovered.

Hence a general procedure to derive variational formulations within the incremental framework for the considered GN coupled thermoelastic model without dissipation is formulated. It is then shown how a family of mixed variational formulations, associated with the considered GN model, can be obtained following a direct and general procedure by enforcing the fulfilment of field equations and constraint conditions.

The possibility to formulate the coupled GN thermoelastic problem in a variational form has a number of consequences and some beneficial effects. For instance, the variational framework allows one to apply the tools of the calculus of variations to the analysis of the solutions of the problem. In particular, conditions for the existence (see, e.g., [16, 17]) and uniqueness of the solution are based on the variational framework.

Accordingly the condition for the uniqueness of the solution of the considered model is provided by means of a minimization principle.

2. Thermoelastic Structural Problem

Throughout the paper bold-face letters are associated with vectors and tensors. The scalar product between dual quantities (simple or double index saturation operation between vectors or tensors) is denoted by . A superimposed dot means differentiation with respect to time and the symbol denotes the gradient operator.

In small strain analysis the theory of thermoelasticity without energy dissipation, as described in Green and Naghdi [18], is considered. Such a theory is based on the introduction of a scalar thermal displacement defined as

where is a point pertaining to a thermoelastic body defined on a regular bounded domain of an Euclidean space, represents the temperature variation from the uniform reference temperature , and is the initial value of the thermal displacement at the time . Accordingly the time derivative of is the temperature variation, that is, .

The mechanical and thermal parts of the thermoelastic model are hereafter defined.

Let denote the linear space of strain tensors and denote the dual space of stress tensors . The inner product in the dual spaces has the mechanical meaning of the internal virtual work, that is

The linear space of displacements is denoted by . The linear space of forces is and is placed in separating duality with by a nondegenerate bilinear form which has the physical meaning of external virtual work. For avoiding proliferation of symbols, the internal and external virtual works are denoted by the same symbol. Conforming displacement fields satisfy homogeneous boundary conditions and belong to a closed linear subspace .

The kinematic operator is a bounded linear operator from to the space of square integrable strain fields . The subspace of external forces is the dual space of . The equilibrium operator is the continuous operator from to which is dual of . Let be the load functional where and denote the tractions and the body forces [12, 19].

The equilibrium equation and the compatibility condition are given by

where , , and , .

The external relation between reactions and displacements is provided by

being a concave function, and the symbol denotes the sub(super)differential of convex (concave) functions [20]. Accordingly, the inverse relation is expressed as

where the concave function represents the conjugate of and the Fenchel’s relation holds

Different expressions can be given to the functional depending on the type of external constraints such as bilateral, unilateral, elastic, or convex. For future reference the expressions of and are specialized to the case of external frictionless bilateral constraints with homogeneous boundary conditions. Noting that the subspace of the external constraint reactions is the orthogonal complement of the subspace of conforming displacements , that is , the functional turns out to be

and a direct evaluation shows that its conjugate is given by:

The proposed framework has the advantage that the formulation of the thermal model is similar to the mechanical one, and the thermoelastic problem turns out to be suitable to build a general variational formulation as shown in the next section.

The linear space of thermal displacement is . The rate of heat flow into the body by heat sources and the boundary heat fluxes belong to the space , dual of , of square integrable fields on . The external thermal forces are collected in the set .

The kinematic thermal operator is a bounded linear operator, and a thermal gradient is said to be thermally compatible if there exists an admissible thermal displacement field such that . The thermal balance equation is given by , where is the dual operator of and is the entropy flux vector.

Constraint conditions can be fit in field equations by noting that the external relation between reactive thermal forces and thermal displacements is provided by the equivalent relations:

being and conjugate convex functionals. The equality (9)₃ represents the Fenchel’s relation.

Homogeneous constraints on thermal displacement fields are formulated by considering that thermal displacement fields belong to the subspace and are said to be admissible. Reactive thermal forces belong to the orthogonal complement of . Then the functional is the indicator of :

and its conjugate is given by

Accordingly the relations governing the quasistatic thermoelastic structural problem without energy dissipation for given mechanical and thermal loads and in the time interval can be collected in the following form:

The initial thermal conditions are considered in the following form:

where and are, respectively, a prescribed initial thermal displacement and temperature in .

The thermoelastic functionals and provide the thermoelastic constitutive relations of the considered GN model.

In the case of a linear coupled thermoelastic behaviour the expression of the functional is

where the first term at the r.h.s. above is the isentropic elastic strain energy and the second one represents the thermoelastic coupling. The parameter is the specific heat at constant strain at the reference state with temperature , is the symmetric and positive definite isentropic elastic moduli fourth-order tensor. The second-order thermal expansion tensor is self-adjoint, that is, . Moreover the thermal function is given by

where is the tensor of conductivity moduli, symmetric and positive definite.

A solution of the thermoelastic structural problem can be achieved by a finite element approach (see e.g., [21, 22]) and can be obtained starting from a suitable mixed variational formulation. Hence the definition of a general variational formulation which allows one to derive variational principles with different combination of the state variables, without ad hoc procedure, plays a central role.

The time integral of the constitutive relation (12)₆ in the interval is given by

where we have set .

Moreover the time integral of the thermal balance equation (12)₂ in the interval is given by

where we have set and . In fact, the constitutive relation (12)₆, evaluated at the initial time , provides the equality which yields the relation .

Accordingly the thermoelastic structural problem (12) is rewritten in the following form where the initial conditions (13) are incorporated into the field equations (12):

3. General Mixed Variational Principle

Introducing the product space and its dual space , the thermoelastic structural problem without dissipation (18) can be collected in terms of the global multivalued structural operators governing the whole problem

The explicit expression of the structural operators , and of the vectors , , are given by:

The operators and turn out to be integrable by virtue of the duality existing between the operators , , , , , , and , , , , , , , the conservativity of , and the conservativity of the super(sub)differentials and .

The related potential can be evaluated by a direct integration along a straight line in the space starting from its origin to get [23]:

Hence it turns out to be

The potential turns out to be linear in , jointly convex with respect to the state variables and jointly concave with respect to . The following statement then holds.

Proposition 1. The set of state variables is a solution of the saddle problem:

if and only if it is a solution of the thermoelastic structural problem without dissipation (18).

The stationary condition of , enforced at the point , yields back the thermoelastic structural problem. In fact the stationarity of is:

which can be rewritten in the following form

From a mechanical point of view, the relation (25)₁ yields the equilibrium equation:

the relation (25)₂ provides the compatibility condition:

the relations (25)₃₋₄ yield the constitutive relations:

the relation (25)₅ yields the external constraint conditions:

the relation (25)₆ yields the thermal balance condition:

the relation (25)₇ yields the thermal compatibility:

the relation (25)₈ yields the thermal constitutive relation:

and, finally, the relation (25)₉ yields the thermal external relation:

Hence, performing the super(sub)differentials appearing in (25), the structural problem (18) is recovered.

3.1. Mixed Variational Principles

A family of potentials can be recovered from the potential by enforcing the fulfilment of field equations and of constitutive relations. All these functionals assume the same value when they are evaluated at a solution point of the GN structural problem.

Hence a variational principle in which the external reactions do not appear as independent state variables can be obtained by imposing, in the expression (22) of the potential , the external relations (18)5,9 in terms of Fenchel’s equalities (6) and (9)₃. Hence it turns out to be

and the following statement holds.

Proposition 2. The set of state variables is a solution of the saddle problem:

if and only if it is a solution of the thermoelastic structural problem without dissipation (18).

The proof of the mixed GN thermoelastic variational formulation reported in the Proposition 2 is provided in the Appendix. The proofs of the following variational formulations are omitted since they are obtained following a similar reasoning.

Imposing in the expression (34) of the potential the compatibility condition (18)₂ follows as

and the following five-field mixed variational principle holds.

Proposition 3. The set of state variables is a solution of the saddle problem:

if and only if it is a solution of the thermoelastic structural problem without dissipation (18).

Imposing in the expression (36) of the potential the thermal compatibility condition (18)₇ follows:

and the following three-field mixed variational formulation holds.

Proposition 4. The set of state variables is a solution of the saddle problem:

if and only if it is a solution of the thermoelastic structural problem without dissipation (18).

From a mechanical point of view, the specialization of the potential to the Cauchy model for a linear coupled thermoelastic behaviour, external frictionless bilateral constrains, and homogeneous thermal boundary conditions is provided. The functionals and , see (14) and (15), are

Hence the potential becomes

under the condition that the displacement and the thermal displacement are conforming and admissible, that is, and , where

is the total potential energy in classical elasticity in terms of the isentropic elastic stiffness, the functionals:

are the mixed thermal part of the potential, and the functionals:

take into account the thermoelastic coupling.

Minimum principles in structural mechanics are relevant to investigate since solution techniques can be exploited, and existence and uniqueness results can be provided by recourse to functional analysis. Actually, uniqueness of the solution is ensured if the functional to be minimized is strictly convex.

Hence it is important to derive a variational formulation in terms of a minimization problem as hereafter reported.

To this end let us note that the constitutive relation (18)₈ can be equivalently expressed in terms of the conjugate of the thermoelastic functional in the following form:

and the Fenchel’s relation hold:

In the case of the Cauchy model with a linear coupled thermoelastic behaviour, the potential can be evaluated as the conjugate of (15) and is given by

Enforcing in the expression of the functional the thermal constitutive relation (18)₈ in terms of the Fenchel’s equality (46) and the external relation (18)₅ in terms of Fenchel’s equality (6), it turns out to be

and the following minimum principle holds.

Proposition 5. The set of state variables is a solution of the minimum problem:

if and only if it is a solution of the thermoelastic structural problem without dissipation (18).

Accordingly, if the thermoelastic functional pertaining to the GN constitutive model is strictly convex, the functional turns out to be strictly convex, and the GN thermoelastic structural model (18) admits a unique solution (if any).

It is worthnoting that, in the linear coupled thermoelastic case, the expression of the functional is given by (47) so that it turns out to be a strictly convex functional. Therefore the GN thermoelastic structural model (18) admits a unique solution. The question of existence of the solution is still an challenge problem; see, for example, [17].

4. Closure

A variational framework for a class of GN coupled thermomechanical boundary-value problem is presented. The thermoelastic structural model is addressed, and the related general mixed thermoelastic variational formulation in the complete set of state variables is derived starting from the structural model. An advantage of the proposed methodology is that it can be applied to a wide range of structural models, and variational formulations can be obtained following a general reasoning. As a consequence, a family of mixed thermoelastic variational formulations in a reduced number of variables is then contributed. Finally, by appealing to a minimum principle, the connection between uniqueness of the solution and convexity is investigated.

Appendix

The stationary condition of , see (34), enforced at the point , yields back the thermoelastic structural problem. In fact the stationarity of is

which can be rewritten in the following form:

The mechanical interpretation of the relation (A.2)₁ yields the equilibrium equation and the external constraint conditions. In fact there exists a reaction such that

The relation (A.2)₂ provides the compatibility condition:

the relations (A.2)₃₋₄ yield the constitutive relations:

the relation (A.2)₅ yields the thermal balance condition and the thermal external relations. In fact there exists a thermal reaction such that

the relation (A.2)₆ yields the thermal compatibility:

and the relation (A.2)₇  yields the thermal constitutive relation:

Hence the structural problem (18) is thus recovered.

Acknowledgments

Research support by “Ministero dell’ Istruzione, dell’Università e della Ricerca” of Italy is kindly acknowledged.

References

S. Bargmann and P. Steinmann, “An incremental variational formulation of dissipative and non-dissipative coupled thermoelasticity for solids,” Heat and Mass Transfer, vol. 45, no. 1, pp. 107–116, 2008.

Q. Yang, L. Stainier, and M. Ortiz, “A variational formulation of the coupled thermo-mechanical boundary-value problem for general dissipative solids,” Journal of the Mechanics and Physics of Solids, vol. 54, no. 2, pp. 401–424, 2006.

M. A. Biot, “Thermoelasticity and irreversible thermodynamics,” Journal of Applied Physics, vol. 27, no. 3, pp. 240–253, 1956.

M. Ben-Amoz, “On a variational theorem in coupled thermoelasticity,” Journal of Applied Mechanics, vol. 32, no. 4, pp. 943–945, 1965.

G. Batra, “On a principle of virtual work for thermo-elastic bodies,” Journal of Elasticity, vol. 21, no. 2, pp. 131–146, 1989.

G. Herrmann, “On variational principles in thermoelasticity and heat conduction,” Quarterly of Applied Mathematics, vol. 22, pp. 151–155, 1963.

G. Lebon, “Variational principles in thermomechanics,” in Recent Developments in Thermomechanics of Solids, G. Lebon and P. Perzina, Eds., CISM Courses and Lectures, no. 262, pp. 221–396, Springer, Wien, Austria, 1980.

J. T. Oden and J. N. Reddy, Methods in Theoretical Mechanics, Springer, Berlin, Germany, 1976.

F. Armero and J. C. Simo, “A priori stability estimates and unconditionally stable product formula algorithms for nonlinear coupled thermoplasticity,” International Journal of Plasticity, vol. 9, no. 6, pp. 749–782, 1993.

F. Marotti de Sciarra, “Novel variational formulations for nonlocal plasticity,” International Journal of Plasticity, vol. 25, no. 2, pp. 302–331, 2009.

F. Marotti de Sciarra, “A finite element for nonlocal elastic analysis,” in Proceedings of the 4th International Conference on Computational Methods for Coupled Problems in Science and Engineering, Kos, Greece, 2011.

F. Marotti de Sciarra, “Hardening plasticity with nonlocal strain damage,” International Journal of Plasticity, vol. 34, pp. 114–138, 2012.

D. S. Chandrasekharaiah, “Variational and reciprocal principles in thermoelasticity without energy dissipation,” Proceedings of the Indian Academy of Sciences: Mathematical Sciences, vol. 108, no. 2, pp. 209–215, 1998.

S. Chiriţa and M. Ciarletta, “Reciprocal and variational principles in linear thermoelasticity without energy dissipation,” Mechanics Research Communications, vol. 37, no. 3, pp. 271–275, 2010.

R. Quintanilla and J. Sivaloganathan, “Aspects of the nonlinear theory of type II thermoelastostatics,” European Journal of Mechanics A, vol. 32, pp. 109–117, 2012.

G. Dal Maso, G. A. Francfort, and R. Toader, “Quasistatic crack growth in nonlinear elasticity,” Archive for Rational Mechanics and Analysis, vol. 176, no. 2, pp. 165–225, 2005.

R. Quintanilla, “Existence in thermoelasticity without energy dissipation,” Journal of Thermal Stresses, vol. 25, no. 2, pp. 195–202, 2002.

A. E. Green and P. M. Naghdi, “Thermoelasticity without energy dissipation,” Journal of Elasticity, vol. 31, no. 3, pp. 189–208, 1993.

R. E. Showalter, Monotone Operators in Banach Space and Nonlinear Partial Differential Equations, American Mathematical Society, Providence, RI, USA, 1997.

R. T. Rockafellar, Convex Analysis, Princeton University Press, Princeton, NJ, USA, 1970.

E. Ruocco and M. Fraldi, “Critical behavior of flat and stiffened shell structures through different kinematical models: a comparative investigation,” Thin-Walled Structures, vol. 60, pp. 205–215, 2012.

E. Ruocco and M. Fraldi, “An analytical model for the buckling of plates under mixed boundary conditions,” Engineering Structures, vol. 38, pp. 78–88, 2012.

G. Romano, L. Rosati, F. Marotti de Sciarra, and P. Bisegna, “A potential theory for monotone multi-valued operators,” Quarterly of Applied Mathematics, vol. 4, no. 4, pp. 613–631, 1993.