Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Holpzaphel_-_Nonlinear-Solid-Mechanics-a-Contin

.pdf
Скачиваний:
225
Добавлен:
27.03.2016
Размер:
24.13 Mб
Скачать

3.4 Alternative Stress Tensors

129

2.The Cauchy stress components at a certain point are given with respe-ct to a x1., :1;rcoordinate system with the associated set {e0 }, a = 1, 2, ·of orthonormal basis vectors as

·[

2

1 ] kN/m2 .

[u] = .

1

-3

 

 

Assume a .rotation of {e0 } into the new set {e0 } of orthonormal bas.is vectors according to

Compute the components of the corotated Cauchy stress tensor u 0

of work or heat,

4 Balance Principles

In this ch~pter we provide the classical balance principles and -discuss some of their important consequences. The fundamentaJ balance principles, i.e. conservation of mass, the momentum balance principles and balance of energy, are valid .in all branches of continuum mechanics. They are applicable to any particular material and must be satisfied for all times.

We also discuss another fundamental set of laws that are expressed as inequalities, such as the second law of thermodynamics. The section devoted to continuum thermo- dynamics specifically addresses balance of energy and the entropy .inequality principle. Finally, the structure -of principles is summarized as the ·master balance (inequality) principle.

4.1Conservation of Mass

Every continuum body B possesses mass, denoted by tn. It .is a fundamental physical

.property com.manly de·fined to be a measure of the amount of a material contained in the body B. In order to perform a macroscopic study we assume .that mass is continuously

(or at .least piecewise continuously) distributed over an arbitrary region n (of physical space) with boundary surface an at time t. The mass .is a scalar measure (a positive number) which is invariant during a motion. We exclude concentrated masses such as those used in class"ical Newtonian mechanics.

~losed and ope.a systems. We define a system as a quantity -of mass or .a particular c~llection of matter in space. The complement of a system, i.e. the mass or region outside the system, we call the surroundings, while the surface that separates the system from its surroundings we can the boundary or wall of the system (se.e Figure 4. l ).

A closed system (or control mass) consists of a fixed amount of mass in a properly selected region n .in space with boundary surface .an which depends on ti.me t (see Figure 4.1). No mass can -cross (enter or leave) its boundary, but energy, in the form can cross the boundary. The volume of a closed system does not have

13.1

132

4 Balance Principles

SURROUNDINGS

Boundary or wall of the system

Control

surface

Fixed amount of mass;

Fixed amount of volume;

.Energy can cross the boundary

Mass and energy can cross

 

the control boundary

Figure 4.1

Closed and open systems~

to be fixed. If even energy does not interact between the system and its surroundings, then we say that the boundary is insulated. Such a system is called mechanically and thermally isolated, which is an idealization for a physical system. There always exist electromagnetic and other types of forces which permeate the space. Note that no physical system .is truly isolated.

An open syste.m (or control volume) consists of a fixed amount of volume of

a properly selected reg.ion flc which .is independent of time t (see .Figure 4.1). The enclosing boundary of a control volume, over which both mass and energy can cross (enter or leave), is called a control surface, which we denote by anc.

Conservation of mass. In non-relativistic physics mass cannot be produced or

destroyed. It is assumed that during a motion there are neither mass sources (reservoirs that supply mass) nor mass sinks (reservoirs that abs-orb mass), so that the mass m of a body is a ·conserved quant-ity. Hence, if a particle has a certain mass in the reference

4.1 Conservation of Mass

133

configuration it must stay the same during a motion. Considering a closed system, obviously that holds for the total mass too. We write

rn..(no) = rn(n) > O ,

(4.1)

which holds for all times :t. Relation (4.1) is a statement of a fundamental mechanical

law known as the conservation of mass. The ·boundary surfaces in the reference and configurations, with volume ll and v, are denoted by n 0 and !1, respectively. Note that the mass -rn is independent of the mot.ion and of the region occupied by the body.

Hence, the material time derivative of the mass 1n gives

 

D

 

D

= o .

(4.2)

 

-n1..(no) =

--rn(n)

 

Dt

·

Dt

 

 

The differential form of eq. (4.1) reads

 

 

 

drn(X) = drn{x, t) > 0 1

(4.3)

·with the infinitesimal mass clement drn.

 

 

 

The mass of n 0 and !1

is characterized by continuous (or at least piecewise con-

tinuous) scalar fields, i.e.

p0 = p0 (X) > 0 and p

= p(x, t)

> 0, respectively. They

denote physic.a·1 properties of the same particle. Property p0 is called the .reference mass density .(or just density) and pis -called the spatial ·mass dens.ity during a motion j( = x(X, -t) .. The spatial muss density, a1so known as the density in the motion, depends on place x E n and ti.met throughout the body. Note that {Jo is time-independent

and intrinsically associated with the reference configuration of the body.

He.nee, p0

.

.

depends only on the position X chosen in configuration 0 0 If the density does not depend ·on X ·E n0 , i.e. if Gradp0 === o, the configuration is said to =be homogeneous.

The mass densities at the points X and x are defined by the limit

Pn(X) =

. .

~7n(no)

..

~rn(.s-2)

(4.4)

hm

~F(n ) ,

p(x, t) = hm

~

(n)

 

.:l F{f?o)-ro / ·"' · ·o.

~u(n)-+O

 

v .

 

\Vhere ~1n denotes a continuous function of incremental mass of an :incremental vol-

:·ume element in the reference and current configurations, which we have denoted by :·Liv·and ~'IJ, respectively.

Note that ~l ..~(.n0 ), il·v (n), actually must not tend :to zero since then the limit of p0 , p, would show a discrete distribution according to the atomistic stmcture of matter. Therefore, to obtain representative averages, .6. ll(n0 ) and ~v(n) must be large in terms of an atomistic scale and s.mall in lenns of a length scale of a certain physical ·problem. Usually the ratio of the length scale of a physical problem and the length sea-le of an incremental volume ele~nent -il ll(!10 ) and ilv(n) is of the order 1oa -or

more.

134

4

Balance Prlndples

 

In the differential form eq. (4.4) reads

 

 

drn(X) = p0 (X)dlt

,

drn(x., t) = p(x, t)dv ,

(4.5)

with the standard infinitesimal

volume elements dl·=" and clv defined in the reference

and current configurations, respectively.

 

 

Substituting eq. (4.5) into (4.3) we o·btain

 

Po(X)dl/ = p(x, t)dv > 0

(4.6)

which means that volume increases when density decreases. By integrating the infinitesimal mass over the entire region, we .find the total mass rn of that region. He.nee, an alternative expression for (4.6) reads

rn = / fJo{X)dV = /

p(x, t)dv =coust > 0

(4.7)

no

n

 

 

for all times t, which impHes the rate form

 

Dni

D .

(4.8)

·1h = -

= -.

p(x, t)dv = 0 .

 

Dt

Dt,

 

 

n

Hence, conservation of mass requires that the material time derivative of m.. is zero for all regions n of a body which change with time (mass remains unchanged during the motion of .n).

An equation which holds at .every point of a continuum and for all times, for ex.ample, eq. {4.6), is referred to us the local (or differential) form of that equation (local means pointwise). An equation in which physical quantities over a certain region of space are integrated is referred to as the global (or integral) form of that equation; see, for example., eq .(4.. 7). Consequently we may say that (4.6) is the local form and (4.7) is the global form of conservation of mass.

.In general, local forms are .ideally suited for approximation techniques such as the finite difference method while global forms are the best to start with when the finite

·element method is employed.

Continuity mass equation. We want to find a relationsh.ip between the reference mass density Po (X) E n0 and the spatial mass density p( x, t) E n.

By recalling eqs. (2.50) .and (2.51), i.e. dv = J(X;t)d l/, J = detF(X, f) > 0, we may change the variable of integration in eq. (4.7) from x = x(X, t) to X and we obtain the identity

/ [po(X} - p(x(X, t), f ).J(X, t)]dY = 0 .

(4.9)

C')

.l .. 0

4.1 Conservation of M·ass

135

By assuming that Vis an arbitrc11J' volume of region n0 , we condude that the integrand in (4.9) must vanish everywhere. Hence,

Po(X) = p(x(X, t_L t).J(X, t) ,

(4.l 0)

holds for .all X E n0 It represents the continuity mass equation (continuity stands for constancy of mass) in the material (or Lagrangian) description which is the appropriate descri pt.ion in solid mechanics.

Since the reference mass dens.ity p0 is independent of time we find simply from

(4.10) that

Dpo(X)

_

. (X) _ O

~

(4. l l)

at

-

{Jo

-

which is the rate form of (4.10) in the material description.

EXAMPLE 4.1 Show how the spatial mass density p = p(x~ t) changes with time. In particular, derive the rate form of continuity mass equation in the spatial (or Eulerian) description, which .is (expressed in terms of the velocity components)

rJ(x, t) + p(x, t)divv{x, t) = 0

 

or

 

 

j1 + p~;~: =

0 ,

(4.12)

or in the two equivalent forms

 

 

 

 

 

 

 

 

 

 

 

 

Dp~~,t) + gradp(x, t) · v(x, t) + p(x, t)divv(x, t) = 0

,

(4.13)

8p(x, t)

 

.

(

 

(

)

(

))

=

o

 

(4.14)

Dt

+div

p

x, t

 

v x, t

.

 

for all x E H and for all times t.

 

 

 

 

 

 

 

 

 

 

 

 

Solution. Since 110 = 0 we obtain from (4.10) that

 

 

 

 

 

D (

pJ

)

 

-· .

 

 

 

 

(4.15)

 

Dt

 

= pJ

= O

 

 

 

 

(for simplicity written without arguments of the scalar quantities). In order to express eq. (4.15) .in terms of the spatial velocity components we find using the product rule

~_nd J = Jdivv, i.e. eq. (2.175){;, that

pJ = i)J + p} = J(p + pdivv) = 0 ,

(4.16)

where the material time derivative of the spatial density function p is, having regard to (225), g.iven by the explicit expression

. Dp

Dp

(4.17)

p= Dt

= Dt +gradp·v

p= 0.
forn1 of (4.8)

136" 4 .Balance Principles

Since J > 0 we deduce from (4.16)2 the desired result (4.12), which is the correspond-

ing local 2. With the material time de.rivative of the spatial density function (4.17) and by means of identity (I ~287) we may obtain from (4.12) the two equivalent

forms (4.13) and (4..14). ·•

.• ,_..,,,,,,,,,, ......... ~, ... , .. , .•.,,... :····:·:-· ..····:,,..,,,,...: ..,,,...................................... ,<C••••."'••''···:,,............. :'"'·'·',..........·:,,....: ............

,,......:'···=·•,:,...••. ,.,.....·:..·,,.,.,,.,..,.....:.-.•,...., ....,.: ..,:.,.,.,: ..-: ..·....,...,..:, .., .....,.:,·:··.. ...:.-=····..·····:,•... ..................

,....,..,,............,....:,.:,:·:,::.,,;o,,,.......• __,,,,.,-; .. ,,.,,

lf the density of a -continuum body is constant at any particle, then from relation (4.12) we find with p = 0 a kine:matical restriction which characterizes un isochoric (volume-preserving) motion, i.e. divv = 0 (compare also with eqs. (2.177)5). A continuum body is inco:mpressible if every motion it undergoes satisfies

The rate forms of continuity mass equations (4.12)-(4. "I 4) show how the spatial mass dens.ity p changes as time changes. They represent the continuity mass equation in the spatial description, which is the appropriate description in.fluid dynamics.

Conservation of ·mass for an open system. Snmetimes we system given by a region nr. and boundary control. surface an(..

At a certain time ta control volume contains the mass ·nz.(t) = the re.gion of integration f2c does not depend on t., integration and

mute and we may write

. .

D /&

(

.)

dv =

/~ 8p(x, t) .

-rn(t)

= -

p x, t

 

.

a .

dv .

 

Dt.

 

 

 

t

 

 

He

 

 

 

n,:

 

 

work with an open

fnc p(x~ t)dv. Since differentiation com-

(4.1.8)

App.lying the divergence theorem for a fixed amount of volume f2c we have using

(l .297)

/

div(p(x, t)v(x, t) )dv = / p(x, t)v(x, t) · nc!H ,

(4.19)

Or

an~

 

where n denotes the outward unit vector field perpendicular to the boundary control surface anc. The tenn .l~nc pv. nds determines the flux of {JV out of nc across Drlc.

Integrating the continuity mass equation in the form of (4.1.4) over a certain region nc and using eqs. (4.18) and (4.19), we obtain the conservation of mass for a control

volume in the global form, .i.e.

Dt .p(x, t)dv = - .p(x, t)v(x, t) . nds .

(420)

flc

one

 

Relation (4.20) .asserts that the material ti.me derivative of the .mass inside a control volume nc is equal to the flux -of pv ente.ring f1e across Df1c:. The global form (4.20) -is

widely used in fluid dynamics.

4.1 Conservation of Mass

137

EXERCISES

1. A velocity field of a plane motion has components of the form

V3 = Q ,

where a and /3 are positive constants. Assume that the spatial mass density p is independent of the curren_t position x so that gradp = o.

(a)Express p so that the continuity mass equation :is satisfied.

(b)Find a condition for which the given motion is isochoric.

2.Tuo motions of a continuum body are given in the form

x = (1 + a:(t)t)X ,

with the scalar function a:(t) and the set {ea}, a = 1, 2_, 3, of orthogonal unit vectors.

Find expressions for the spatial mass density pin terms of p0 so that the continuity mass equation is satisfied.

3.Consider a spatial scalar field <I> = <I>(x, t) and a spatial vector field u = u(x, t). Use the rate form of the continuity -mass -equation and obtain the identities

p - =

· + div{p<I>v)

 

p Du = a~u) + div(pv ® u) , (4.21)

D<I>

D(p<I>)

 

 

 

Dt

8t

'

Dt

.t

where v denotes the spatial velocity field.

4.An irrotational motion of an incompressible continuum is given by the s.patial velocity field v = -grad<I>. Show that for this case the scalar field <I> is harmonic.

5.By means of continuity mass equation (4..12) and identity (2.181.) show that the vorticity vector 2w is .related to the spatial mass density p and to the spatial velocity vector v by

D(2w)

Dv

+ (gradv)2w ,

p -.

-

p

= curl D

t

0 t

 

.

 

which is known as the Beltrami vorticity ·equation.

.138

4 Balance Principles

4.2 Reynolds' Transport Theorem

Suppose we have a spatial scalar field <I> = <I>(.x, t) describing some physical quantity (for example, mass, internal energy, entropy, heat or entropy sources) of a particle in space per unit volume at time t:.

Assume <I> to be smooth, so that it is continuously differentiable. Hence, the present status of a continuum body in some three-dimensional region n with volume v at given time t may be characterized by the scalar-valued function

l(t) = / <l>(x, t)dv .

(422)

n

 

The aim is now to compute the material time derivative of the volume integral I(t). Since the region of integration n depends on time t, integration and time differentiation do not commute. Therefore, as a first step I(t) ·must be transformed to the reference configuration. By .changing variables using the motion x = x(X, -t) and the relation dv = J(X_, t)dll we find the time rate of change of I(t) to be

.

o;·

Dr

.

~

(4.23)

<I>(x, t)d'IJ = Dt.

<I>(x(X, t), t)J(X, t)dl' .

 

l(t) =

Dt.

 

 

 

 

 

n

no

 

 

 

Since the region of integration is now time-independent, integration and differentiation commute. Hence, as a second step, from (4.23h we obtain, using the product rule of differentiation,

gt / <I>(x, t)dv =

/ [<i>(x(X, t), t)J(X, t) + <I>(x(X, t), t).i(X, t)] dV , (4.24)

n

no

where 4~ denotes the material time derivative of the spatial scalar fie.Id <I> .according to relation (2..25). In a last step we undo the change of variables and convert the volume integral back to the current configuration. By means of eq. (2.l 75)fi, dv = J(X, t)d\/ and motion x = x(X, t), we find finally that

D ;·<I>(x, t)clv =

;·(<i>(x(X, t), t) + <I>(x(X, t), t) ·~~~·t~]J(X,t)dF

D.t,

'

 

f

,t

n

nu

 

 

.

 

 

.[.

}(X, t) J

 

 

./

<I>(x, t) + cl>(x, t) J(X, t) dv

 

 

n

 

 

 

 

/

(<i> (x, t) + cl>(x, t}divv(x, t)]clv

(4.25)

n

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]