Holpzaphel_-_Nonlinear-Solid-Mechanics-a-Contin
.pdf
|
4.4 Balance of Mechanical Energy |
159 |
|
manipulations .according to (1.95) that |
|
||
/ |
P : Fc1v = / |
PFT : RRTdv + / aTp : Uc1v |
|
no |
no |
flo |
|
|
/ |
r : RRTdV +/Tu :iicIV |
(4..112) |
|
no |
no |
|
Knowing that r is a symmetric tensor and RR'T' is a skew tensor accordin_g to (2.158), we conclude by analogy with property ( 1.1.17) that r : RRT = 0, consequently the first term in (4.112}~ vanishes. Since U is a symmetric tensor the second term gives the desired result (4.11 Oh. II
In summary, the alternative expressions for the .rate of internal mechanical work
Pint a.re the relations (4.1.07_)i, (4.IOl)i, (4.107)4., (4.108h, (4.l09h and (4.1 lOh, and we have finally the important identities
. |
. |
|
.1 |
-1 . |
w;n1 (t) = Ju : d = P: F = S: E |
= 'E: |
2c |
C |
|
= Juu : Dll = sy1nTn : U . |
|
(4."I 13) |
||
The double contraction of a stress tensor and the associated rate of deformation tensor (:for example, P : F) describes the re'll physical power during a dynamical process, i.e. the .rate of internal mechanical work (or stress power) per un.it reference volume, denoted by w1nt. In this sense the stress fields Ju, P, S, ~, Juu, syn1Tu are said to be work conjugate .to the strain "fields d, F, E, c~1 c/2, D1t, U, respective·ly. Hence, for example, P and F is said to be a work COJ\jugatc pair (for a more comprehensive survey of work conjugate pairs see, for -example, ATLURI [1984.J).
Conservative system. Let the scalar-valued functions Ilext .and flint be called the potential energy of the external loading (or the external potential energy) and the total strain energy (or the internal .potentia~ -energy) of the body, respectively.
The sum of IIl~xt and we cal I the total potential energy (or frequently referred to as the energy functional) I1 of the mechanical system_, i.e.
(4..11.4)
A mechanical system. is known as conservative if the external mechanical power Pffxt expended on a certain reg.ion of the body and the internal mechanical -power P1nt ·are expressible as
Pc::.t(t) = - |
DTI(~Xt("t) |
. |
( . |
(4.U5) |
Dt |
= -Ilc::.t |
t) , |
160 |
4 Halance Principles |
|
|
with |
Ilim(t) =IWdV , (4.116) |
|
|
no |
where \JI represents the s·tra·in-energy function (referred to briefly as the strain -energy or stored energy) defined subsequently per unit reference volum.e rather than per unit mass .and dealt with in Section 6J in more detail. In the literature the strain-energy function is sometimes denoted by
By means of the assumptions (4.115) and (4.116), the balance of mechanical energy (4.85) implies that the sum of the total potential energy IT and the kinetic energy K is
conserved (constant) during a dynamical process (u, x). We write |
|
Ilexi (t) + Ilint ( t) + IC(t) = canst ~ |
(4.117) |
.It is important to note, that, for example, most surface tractions and al.I problems associated with dissipation of energy (due to external or intenrnl friction, viscous or pfastic effects) lead to non-mechanical energies which are non-conservative in Lhe sense that they cannot be derived from a potential.
EXERCISES
1.lf only mechanical energy is considered, show that the energy equation is merely
Cauchy's first equation of motion. Thus, take the dot product of each term of divu + b = pv with the spatial velocity field ·v and integrate the result over the current region .n with volume v.
2. A hydrostatic stress state at a certain point is given by the Cauchy stress tensor in the form of u =
(a) Show that the stress.._power Wint per unit referential volume is given by
·. |
JDp |
· |
|
\ |
\Vinl = JCT : d = ]J- D |
t |
|
. |
µ |
|
|
(b)Formulate the balance!~of mechanical energy for this particular stress state
and conclude that if the considered motion is isochoric and if there exists no external mechanical power, the kinetic energy· is conserved (K = const).
3.A rigid-body is rotating about a fixed point 0 with angular velocity w. Show that the kinetic energy may be expressed as
JC(t) = ~w · Dw , |
D =Ip((x · x)I - x ® xjd'IJ , |
|
n |
where the current position of an arbitrary point in the rig-id-body is characterized
by x. Note that.Dis called the lnertia tensor of the body relative to 0 occupying a region n..
4.5 Balance of Energy in Continuum Thermodynamics |
161 |
4..5 Balance of Energy in Continuum Thermodynamics
In this section we consider both mechanical and thermal energy which are essential in many problems of physics and engineering. The effect of other energies (such as electric, chemical ...) on the behavior of a continuum is a topic beyond the scope of this text.
In a thermodynamic context the conservation of mass and the momentum balance principles are supplemented by the balance of energy and the entropy inequality law.
The development of the concepts of energy and entropy has been one of the most important achievements in the evolution of physics. For an extens·ive account of these conce.pts see, for example, TRUESDELL and TOUPIN [ 1960, Chapter E], the two-volume work of l<.ESTIN [1979], or SILHAVY [.1997] and W:rLMANSKI [1998].
In this section we discuss briefly the balance of energy equation in a thermodynamic context. To begin with, we :introduce some important terminology common in thermodynamics.
General notation. A continuum which possesses both mechanical and thermal en- ergy is called a thermodynamic continuum. We say that the thermodynamic state or the condition of a sys.tern is known if all quantities throughout the entire system are known. All quantities characterizing a system .at a certain state are cal.led thermody- namk state variables. They are macroscopic quantities and, in general, they depend on position and time. For example, the thermodynamic state of a thermoelastic solid can be represented in a seven-dimensional state space with six variables corresponding to the strains and one to temperature. The function that describes a certain state yari- able is called a thermodynamic state function. Any equation that interrelates state variables is called equation of state or constitutive equation.
Any change that a system undergoes from one thermodynamic state to another is called a thermodynamic process. The path that connects the two states is parameter- ized by time t. If a system returns to its initial state at the end of a thermodynamic process this system is said to have undergone a so-called cycle (initial and final states are identical).
A non..equili=brium state is a state of imbalance (there exists a gradient of temperature and velocity), while an equilibrium state is a state of balance (of uniform temperature and zero velocity). To study non-equilibrium or equilibrium states is a central goal of continuum thermodynamics.
A system within equilibrium has no tendency to change when it .is isolated from its surroundings. If there is no change in the values of the state variables at any :particle of the system with time we say that this system is in thermodynamic equilibrium or ~hermal equilibrium. The transition from one state of thermodynamic equilibrium to another is studied in thermostatics. A process in a system that remains close to a state
162 4 Balance Principles
of thermodynamic equilibrium at each time is referred to as a quasi-static process or quasi-equilibrium process. A quasi-static process is a sufficiently slow process whereby enough time remains for the system to adjust itself internally. It is basically a process during which the system is in equilibrium at all times; the contributions due to dynamical quantities are negligible.
Heat is the form of (thermal) energy that is transferred between a system and its surroundings (or between two systems) by virtue of a temperature gradient. The term 'heat' is understood to mean heat transfer in thermodynamics.
Thermal power. Let the thermal (non-mechanical) power or the rate of thermal work be denoted by Q and defined by
Q(t) = f qnds + f rd'U = f QNdS + ./ Rdlf , |
(4.118) |
|||
an |
n |
ano |
no |
|
which is represented in the spatial and material descriptions, respectively.
The time-dependent scalar functions qn and QN denote heat fluxes, determining heat per unit time and per unit current and reference surface area, respectively. The total heat fluxes ./~n qnds and fano QNdS measure the rate at which heat enters (inward normal flux) the body across the current and the reference boundary surfaces 80. and ano, respectively.
The time-dependent scalar fields r = r(x, t) and R = R(X, t) in eq. (4.118) denote heat sources per unit time and per unit current and reference volume, respectively (see Figure 4.3). A heat source is a reservoir that supplies energy in the form of heat. The total heat sources Jn rdv and fno Rdlf measure the rate at which heat is generated (or destroyed) into a certain region of a body.
The counterpart of Cauchy's stress theorem (3.3) in continuum mechanics is the Stokes' heat flux theorem in thermodynamics. It postulates that the scalar functions q11 and QN are linear functions of the outward unit normals so that
qn (x, t, n) == -q(x, t) ·n |
or |
qn = -q0 n" , |
QN(X, t, N) = -Q(X, t) · N |
or |
} (4.119) |
|
The time-dependent spatial vector field q = q(x, t) is the so-called Cauchy heat flux (or true heat flux) defined per unit surface area in n, and n is the outward unit normal to an infinitesimal spatial surface element ds E an at the current position x. The Piola-Kirchhoff heat flux (or nominal heat flux) and the outward unit normal to an infinitesimal material surface element dS E an0 at X are denoted by the vectors Q = Q(X, t) and N, respectively. The time-dependent material vector field Q determines the heat flux per unit surface area in no.
The negative signs in eqs. (4.119) are needed (see Figure 4.3) because the unit vectors n and N are outward normals to an and respectively. However, we claim
in (4.118) that heat enters the body (inward normal flux).
4.5 Balance of Energy in Continuum Thermodynamics |
163 |
n
Figure 4.3 Heal flux vectors q, Q and heat sources r, R.
In order to relate the contravariant vectors q and Q to one another, we equate the total heat fluxes J~no QNdS and Jan qnds. Hence, with the fundamental Stokes' heat flux theorem (4.119), Nanson 's formula (2.54) and manipulations according to ( 1.95)
we obtain |
.!Q · NdS = .!q ·nds = |
.!q ·JF-TN<lS |
|
|
|
|
|||
|
Dno |
an |
ano |
|
|
|
1 |
|
(4.120) |
|
= |
.! .JF- q · |
NdS . |
|
|
|
|||
ano
Comparing the left and right-hand sides of (4.120) we conclude that the Piola-Kirchhoff heat flux Q may be related to the Cauchy heat flux q. We have, by analogy with (3.8),
Q = JF- 1q = Jx; 1(qP) or CJA = JF;\:q(l ' (4.121)
which may be viewed as the pull-back operation on the contravariant vector q~ by the inverse motion x- 1 scaled by the volume ratio J. Such a transformation we called the
Pio/a transformation, here relating the heat flux vectors Q and q (compare with the
definition on p. 84).
4.5 Balance of Energy ·in Continuum Thermodynamics |
165 |
First law of thermodynamics in material description. In order to rewrite balance of energy (4.124) in terms of material coordinates, we recall the equivalent forms of the total energy, the externa] mechanical power 'P(~xt and the thermal power Q, that are given by eqs. (4.95), {4.105), and (4.82), (4.98) and (4. I 18), respectively.
Thus, the first law of thermodynamics in the material description reads
gt/(~poV2 +e)dF= /(T·V+QN)ds+/(B·V+R)dF . (4.126)
flo Oho flu
Following arguments analogous to those which led to (4.125), we find the reduced
global form of balance of ·energy in the material description, i.e. |
|
|
g1·ed1/ := |
;·(P: F- DivQ + R)cnr . |
{4.127) |
t ~ |
. |
|
-!lo |
no |
|
I.n order to achieve the local form we must rewrite the term on the left-hand side of eq. (4.127). Since the reference volume 1/ is independent of time we may write D/Dt J~20e dl/ = J~20 edl..... Note that the volume is arbitrary, leading to the local form of the balance of -energy .in the =material description
e= P : F - DivQ + R or (4.128)
presented in symbolic and index notation, .respectively..
The first-order partial differential ·equation ·is due to Kirclzho.fl and holds at any particle of the body for all times.
EXERCISES
l.Assume a continuum body not subjected to body forces and heat sources. Establish the balance equations
Dt.D !(2l r>v.,- |
+ ec)clv = .!(uv - q)·nds , |
|||
D 1·n l 0 |
.> |
an1· |
"I' |
Q) · NdS . |
Dt. (2'p |
V: |
+ e)dV = . |
(P V - |
|
no |
|
/Jiln |
|
|
These relations -equate the rate of -change of the total energy of a continuum body occupying a certain region with the total flux out of the boundary surfaces of that region. The vectors uv - q and pTv - Q are often referred to as the energy flux vectors with respect to the -current and .reference configurations, respectively.
166 |
4 Balance Principles |
2.Starting from (4..124), show that the local spatial form of the balance -equation has the form
a ( i ,, |
) |
. |
1 .., |
. |
|
Dt |
2pv- + ec |
· = chv[uv - q - |
(2'pv- |
+ ec)v1+ b. v + r . |
|
Hint: Use eqs. (2.50), (2.25), (2.175)n and the rule .according to (I .287) for the. left-hand side of (4.1.24).
3.Derive the "local form of the balance equation in the spatial description -corre- sponding to eq. (4.128) and note that we have used the internal energy per unit volume rather than per unit mass (see eq. (4.91)).
4.6Entropy Inequality Principle
The first law of thermodynamics governs the energy transfer within a thermodynamic process, but is insensitive to the direction of the energy transfer. In the .following we derive a much finer principle, that is the second law of thermodynam.ics, which is respo11sible for the direction of an -energy transfer process. We introduc·e the concept of entropy and discuss some of the consequences of the second law of thermodynamics.
Second -law of thermodynamics. Physical observations show that heat always flows from the warmer to the colder region of a body (free from sources o:f heat), not vice versa; mechanical energy can be transformed into heat by friction, and this can never be converted back into mechanical energy.
We introduce a fundamental state variable, the entropy (c01ning from the Greek words i11 and Tpomx» meaning 'in' and 'turning, direction', respectively). It is an important thermodynamic property first described in the works of Clausius .in the second half of the nineteenth century.. The entropy can be viewed as the quantitative measure of microscopic randomness and disorder (see, for example, CALLEN [1985, Chapter 17]).. A physical interpretation ·is provided by the subject of statistical mechanics. For a more detailed physical interpretation the reader is referred to Section 7. l .of this text. We introduce the notations T/c = T/c(X, t) and "7/ = ·11(X, t) for the entropy per unit cur-
.rent and reference volume, respectively, at a certain point and time t. I.n other texts the entropy is frequently designated by s or S..
The entropy possessed by a continuum body occupying a certain region, denoted by S, is defined to be
S(t) = f ·11c(x, t)c:lv = f 11(X, t)dV , |
(4.129) |
|
n |
no |
|
4..6 Entropy Inequality Principle |
l67 |
with -17(X, t) = J(X, t)rJc(x, t), which should be compared with the analogues of (4.1.03) and (4.104), respectively. This definition .is commonly introduced in continuum mechanics and differs from the definitions used in statistical mechanics..
Let the rate ·of entropy input into a certain reg.ion of a continuum body consist of the value of entropy transferred across its boundary surface and the entropy generated (or destroyed) inside that region. We denote it by Q and write
Q(t} -:-- - j |
h · ncls + ./ hlv = - |
./ H · NdS + ./ R<lV |
. |
(4.1.30) |
|
on |
n |
fmo |
no |
|
|
The time-dependent scalar fields f |
r(x, t) |
and R = .R(X., t) |
denote entropy |
||
sources per unit time and per unit current and reference volume, respectively. The total entropy sources J~ 'fdv and ./~20Rdv· measure the rate at which entropy is generated (or destroyed) into a region of a body. The time-dependent vector fields h == h(x, t) and H = H(X, t) determine the Cauchy entropy flux (or true entropy Hux) defined per unit current surface area .in .n and the Piola-Kirchhoff entropy Hux (or nominal entropy flux) per unit reference surface area in :.n0 , respectively.
As usual, the ·Outward unit normals to the infinitesimal surface elements ds E (JO, at x and dS E DQ0 at X are denoted by n and N, respectively.. Since we have defined the unit vectors n and N to be outward normals to 80 and an0 and since we compute the rate of entropy inplll (entropy entering the body}, we need the negative signs in relation
(4.130).
The .difference between the rate of change of entropy S and the rate of entropy input Q into a body determines the total production of entropy per unit time, which we denote by r. We postulate th.at the total entropy production for all thermodynamic processes is never negative, following the mathematical expression
D . |
- |
(4. l.3l) |
r(t) = n-/(t) - |
Q(t) > o , |
which is known as the second la''T of thermodynamics. Unlike the considerations in previous sections, it is an inequality rather than an equation often referred to as the ·entropy inequality principle.
·More explicitly, with eqs. (4.129) 1 and (4.130) 1 we find from {4.131) the global spatial form of the second law of thermodynamics in the notation of continuum mechanics, Le.
r(t) = -D 1· |
·1Jc(X, t)dv + 1·h · nds - |
1·Pdv > 0 . |
(4.132) |
0 t ' |
' |
' |
|
n |
lJn |
n |
|
These relations assert clearly a trend in time by describing the dire.ction of the energy transfer and postulating irreversibility of various thermodynamic processes.
