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

Holpzaphel_-_Nonlinear-Solid-Mechanics-a-Contin

.pdf
Скачиваний:
213
Добавлен:
27.03.2016
Размер:
24.13 Mб
Скачать
49
= 8ui/axi.

1.8 Gradients and Related Operators

Solution. We denote the skew part of gradu as W, with components H'ii

With the index relation H'ii = -lVii for skew tensors and (l.124)2, we deduce from (l.272h that

(1.276)

 

Gradient of a vector field. The gradient (or derivative) of a smooth vector.field u(x) is defined to be a (second-order) tensorjield. It is denoted by gradu (or Vu, or

au/ax or sometimes more explicitly by V

® u) and according to (1.264)3 is given as

gradu = V

 

8-ui

ei ® ei ,

(1.277)

® u = -

 

 

 

8

 

 

 

 

 

 

J

 

 

 

with Cartesian components (gradu)ij = 8uif8xi. In matrix notation this reads as

 

 

8u1

a'U1

8u1

 

 

 

 

 

 

 

 

Bx1

ax2

8xa

 

 

[gradu] =

 

8u2

a'U2

8u2

 

( 1.278)

 

-

8x2

-

 

 

 

8:1:1

8:c3

 

 

 

 

8u3

GU3

au3

 

 

 

 

a:r1

8x2

8xa

 

 

 

 

 

 

 

 

 

It is easily seen from eq. ( 1.278) that the sum tr(gradu) of the diagonal elements of the matrix is equivalent to the divergence of u characterized in eq. ( 1.271 ). Thus,

tr{gradu) = gradu : I = divu .

(1.279)

In addition, the transposed gradient of a smooth vectorfield u(x) of position xis denoted by gradTu (or uV, or (8u/8x)T or sometimes more explicitly by u ® V). It

1s given as

au

Bu·

(1.280)

gradTU = U 0 V = ei 0 -

= a.1 ei ® ej ,

8Xi

:ti

 

with Cartesian components (gradTu)ij = 8ui/8xi.

 

Divergence and gradient of a (second-order) tensor field.

The vector operator

V dotted into any smooth (second-order) tensor.field A(x) is a vector.field denoted by divA (or V ·A), which is the divergence of A(x). With (l.264)i, representation (1.59) and properties ( 1.53}, ( 1.21) and ( 1.61) we have

(1.281)

8Aii/8xiei

50 1 Introduction to Vectors and Tensors

Note that, especially in the mathematical community, the alternative definition divA =

of the divergence of a second-order tensor field A is frequently used.

The gradient (or derivative) of a smooth (second-order) tensor field A(x) is defined to be a (third-order) tensor.field. It is denoted by gradA (or VA, or DA/Dx, or

more explicitly by V ®A). Using ( 1.264)3 and representation ( 1.59), we may write

DA··

( 1.282)

gradA = V ®A = a.11

ei ® eJ ® eh .

Xk

 

 

Laplacian and Hessian. The vector operator V dotted into itself gives the Laplacian (or the Laplacian operator) \72 (or V · V, or in the literature sometimes denoted by ~). With the rules ( 1.261 )1, ( 1.264)1 and properties ( 1.2 I), ( l.22h we find that

" .

V ·V(•)

a(.)

a2 (.)

 

 

v-(•) =

== V ·- ei =

ei ·

 

 

 

D:ri

Dxia:c:i

3

 

 

 

 

 

 

 

= a~(·~ oij = a2(~) .

 

( 1.283)

 

 

8xi83-j

8.'ti

 

 

The Laplacian V2 operated upon a scalar field cI> yields another scalar field. One example is Laplace's equation \72 cI> = 0. Another example is the inhomogeneous Laplace's equation often referred to as Poisson's equation, i.e. \72 <I> = w. The differential equations of Laplace and Poisson occur in many fields of physics and engineering, for example, in electrostatics, in elasticity theory, in heat conduction problems of solids and so on.

The operator VV (or more explicitly denoted by V ® V) is called the Hessian.

With the rules ( 1.261 )i and ( 1.264h we obtain

 

 

a(.)

a2 ( . )

 

VV(•)=V0V(•)=V®-a· ei= a a· ei®ei

( 1.284)

 

 

:Ci

Xi Xj

 

for the Hessian. Note the properties

 

 

curlcurlu =grad (divu) - \72 u

,

( 1.285)

\7

2

?

?

( 1.286)

 

(u · v) = v-u ·v + 2gradu · gradv + u · v-v .

If a vector field u is both solenoidal (divu = O} and irrotational (curio =

o), then

\7 2 u = o (see eq. ( 1.285)). For this case the vector field u(x) is said to be harmonic. If a scalar field cI> satisfies Laplace's equation \72 cI> = 0, then cI> is said to be harmonic.

EXERCISES

1. Consider a smooth scalar field <I> = x 1:z:2:r3 - :r: 1Find a vector n of unit length normal to the level surface cl>= canst passing through {2, 0, 3).

 

 

 

 

 

1.8

Gradients and Related Operators

51

2.

Assume a force with magnitude F acting in the direction radially away from the

 

origin at point (2a, 3a, 2v'3c) on the surface of a hyperboloid cI>(x) == cI>(1: 1, 1:2 ,

 

'JI

a

2

 

., I

.,

2 I "

 

..

 

:1:3 )= :r:1

 

+ 1:2 er -

:z::1 c- = 1. Compute the components of the force located

 

in the tangent plane to the surface at the point (2a, 3a, 2v'3c).

3.

Consider a scalar field <I>{:r1 , :v2} = e:ra cos(3:i: 1 -

2:r2 }. Compute the directional

 

derivative of cI>

in

the

direction of the line :1: 1

= :r2 /3 -

1 at the point with

coordinates (0, 1), for increasing values of:1: 1

4.If U = u(x) == :r1:1:2J:ae1 + :r1:r2C2 + :r1 e3, determine divu, curio, graclu, V72u, respectively. Verify that (1.285) is satisfied.

5.Apply the operator V to products of smooth scalar fields cl>, \JI, vector fields u, v and tensor fields A. Establish the important identities

div( cI>u)

=

<I>divu + u · grad <I>

,

 

( 1.287)

div{cpA)

=

cI>divA + Agrad<I>

,

 

( 1.288)

div(ATu)

= divA· u +A : gradu

,

( 1.289)

div(u x v)

=

v ·curio - u · curlv

,

 

 

div(u®v)

= (gradu)v+udivv

,

 

( 1.290)

grad (<I>\JI)

= (gradcp) W + <I> grad \JI

,

 

gTad(<I>u)

=

u ® grad<I> +<I> gradu

,

 

grad{u · v)

=

{gracfru) v + (gradTv) u ,

( 1.291)

curl{ cI>u)

= grad<I> x u +<I> curlu

,

 

curl(uxv) =udivv-vdivu+(gradu)v-(gradv)u.

(1.292)

6. Consider the inverse square law defined by the particular vector field

x

u =a lxFI '

with the positive constant a and x = :ciei. The region of u is the three-dimensio- nal space, excluding the origin <lxl # 0). Note that in fluid mechanics this relation characterizes the velocity field of a point source.

(a) Show that the vector field u(x) is harmonic.

(b} Derive a scalar field cl> whose gradient is u. (<I> == -a/lxl>

52

1 Introduction to Vectors and Tensors

1.9 Integral Theorems

In the following we summarize some results of the integral theorems of Gauss and

Stokes which are of essential importance in the field of continuum mechanics.

Divergence theorem. Suppose u(x) and A(x) are any smooth vector and tensor fields defined on some convex three-dimensional region in physical space with volume v, and on a closed surfaces bounding this volume (see Figure 1.10).

n

Surface element ds

" Closed surface s

Volume v

Figure 1.10 Volume v and bounded closed surface s (with infinitesimal surface element ds and associated unit normal vector n).

Then, for u and A we have - without proof

{'Uinids = J8u·- 1 dv

 

 

/

u · nds = / divudv

or

'

(1.293)

 

v

 

.

v

axi

 

s

 

s

 

 

 

JAnds = JdivAdv

or

I

Ai ·n ·ds =

-- dv

(1.294)

 

1 1

I aA;j

 

 

 

 

 

 

 

8x·

 

 

s

 

 

s

 

J

 

 

ti

 

 

v

 

 

(given here in both symbolic and index notation), where n is the outward uni1' normal field acting along the surface s, dv and ds are infinitesimal volume and surface elements at x, respectively. The important transformation of a surface integral into a volume integral (i.e. (1.293), (1.294)) is known as the divergence theorem (or Gauss' divergence theorem). Similar results hold for higher-order tensors.

The surface integral J'l u · n<ls in the expression (1.293) is called the (total) flux

(the name goes back to Mcu·well) of u (or the outward normal flux) out of the total

boundary surface s enclosing 'V.

 

 

 

 

By setting A= cI>I in theorem (l.294) and knowing that div(cI>I) =

grad<I>, recall

identity ( 1.288), we consequently obtain from ( 1.294)

 

 

J<I>nds = Jgrad<I>dv

or

J<I>nids = J;::dv

,

( 1.295)

s

v

s

v

Stokes' theorem.

1.9 Integral Theorems

53

which is known as the Green-Gauss-Ostrogradskii theorem.

For the sake of completeness the theorem of Stokes is briefly summarized. It relates a surface integral, which is valid over any open surface s, to a line integral around the bounding closed curve c in three-dimensional space.

We now introduce a tangent vector to c denoted by dx (with components chi) and an outward unit vector field n normal to s. Curve c has positive orientation relative to n in the sense shown in Figure 1. I I .

The indicated circuit with the adjacent points 1, 2, 3 (1, 2 on curve c and 3 an interior point of the surface s) induced by the orientation of c is related to the direction of n (i.e. a unit vector normal to sat point 3) by the right-hand screw-rule (see Figure 1. I I).

n

Surface element ds

dx

Surfaces

Curve c

Figure 1.11 Open surface.

For a smooth vector field u defined on some region containing s, we have

f u · dx = Jcurio · ncls

or

 

( 1.296)

 

 

 

c

 

c

s

(the proof is omitted).

Note that if the surface s is closed the integral on the left-hand side of eq. ( 1.296) reduces to zero, and then Is curio · nds = 0. The line integral fc u · dx represents the circulation of u around the closed space curve c.

EXERCISES

1. Let cl>, u, and A be smooth scalar, vector and tensor fields defined in v and on s

54

1 Introduction to Vectors and Tensors

and let n be the outward unit normal field acting along the surface s. Show that

/

cI>u · nds =

/

div(<I>u)dv

,

( 1.297)

s

 

v

 

 

 

 

/

<I>Ands =

/

div(cI>A)d-u

,

( 1.298)

 

 

v

 

 

 

 

/

n x uds = /

curludv

,

 

 

 

 

v

 

 

 

 

/

u 0 nds = /

gradudv

,

 

 

 

 

ti

 

 

 

 

/

u ·Ands= /

div(ATu)dv .

(1.299)

s

 

v

 

 

 

 

2.Let u and A be smooth vector and tensor fields defined in v and on .s and let n be the outward normal to s. Use index notation to show that

/ u x Ands=/ (u x divA+£ : AT)dv ,

(1.300)

v

where e is the third-order permutation tensor as expressed in eq. ( 1.143 ).

3.Let cI> and u be smooth scalar and vector fields defined on s and c and let n act on s. Show that

.f<I>dx =

/ n x gradcI>ds ,

c

s

f u x dx =

/[(divu)n - (gradTu)n]ds

c

s

2 Kinematics

In the real world all physical objects are composed of molecules which are formed by atomic and subatomic particles. A microscopic system is studied by means of magnifying instruments such as a microscope. Microscopic studies are effective at the atomic level and very important in the exploration of a variety of physical phenomena. The atomistic point of view, however, is not a useful and adequate approach for common engineering applications. The microscopic approach is used in this text only briefly within Chapter 7.

We use the method of continuum mechanics as a powerful and effective tool to explain various physical phenomena successfully without detailed knowledge of the complexity of their internal (micro)structures. For example, think about water which is made up of billions of molecules. A good approximation is to treat water as a continuous medium characterized by certain field quantities which are associated with the internal structure, such as density, temperature and velocity. From the physical point of view this is an approximation in which the very large numbers of particles are replaced by a few quantities; we consider a macroscopic system. Hence, our primary interface with nature is through these quantities which represent averages over dimensions that are small enough to capture high gradients and to reflect some microstructural effects. Of course the predictions based on macroscopic studies are not exact but good enough for the design of machine elements in engineering.

The subject of continuum mechanics roughly comprises the following basic ingredients:

(i)the study of motion and defonnation (kinematics),

(ii)the study of stress in a continuum (the concept of stress), and

(iii)the mathematical description of the fundamental laws of physics governing the motion of a continuum (balance principles).

The aim of the fallowing three Chapters 2-4 is to derive the essential equations within the basic fields (i)-(iii). Note that the provided results are applicable to all

55

2 Kinematics

classes of materials, regard.less of their internal physical structure. In -order to explain the macroscopic behavior of physical objects, :first of all we must understand the motion and deformation that cause stresses in a material (or are caused by stress.es) arising from forces and moments. Thus, to study the motion and (finite) d.eformation of a continuum, i.e. kinematics, is essential and mainly the aim -of the following chapter.

The reader who requires additional information on the subject of this -chapte.r i-s referred to the monographs by, for example, MALVERN [ l 969], WANG and TRUESDELL

.(1973], CHADWICK [1976], TRUESDELL [1977], GURTIN [1981.a], CIARLBT :[1988],

TRUESDELL and NOLL [1992], MARSDEN and HUGHES [.1994) or OGDEN [1997].

2.1Configurations, and Motions of Continuum Bodies

Under an electron microscope we see the discontinuous atomic structure of matter. The molecules may be crystalline or .randomly oriented. Between the particles there are large gaps. Theories considering the discrete structure of matter are ·molecular and atomistic theories. They are based on a discrete particle approach. For macroscopic systems such theories tend .to become too complicated to yield the desired results and therefore do not.meet our needs. However, under certain circumstances the microscopic approach is indispensable for the study of physical phenomena. A review of atomistic models is :presented in the work of ORTIZ f1999, and references the.rein_].

Notation of a particle nod a continuum body. Macroscopic systems usually can be described successfully with a continuum a_pproach (macroscopic .approach). Such an approach leads to the continuum theory. The continuum theory has been developed .independently of the molecular and atomistic theory and is meeting our needs. A fundamental assu-mp.tion therein states that a body, denoted by B, may be viewed as having a continuous (or at least a piecewise continuous) distribution of matte.r .in space and time. The body is imagined as being a composition of a (continuous) set of parti.. cles (or continuum particles or material points), represented by P E B as indicated

. p·

'JI

.

m. 1gure -·

lt is important to note that the notion 'particle' (or 'continuum .particle' or 'material point') refers to a part of a body and does not imply any assoc.iation with the point mass -of Newtonian mechanics or the discrete particle of the atomistic theory as .mentioned above. A typical continuum particle is an accumulation of a large number of molecules, yet is small enough to be considered as a particle. The behavior of a continuum partic.le is a consequence of the collective behavior of all the molecules constituting that particle.

Hence, in a macroscopic study we are concerned with the mec.hanics of a body in which both mass and volume are continuous (o.r at least piecewise continuous) .func-

ORTIZ :.[1999].

2..1

Configurations, and Motions of Continuum Bodies

57

Reference

Current

 

configuration

 

configuration

 

 

.

.·.

. . .

.

.

,;:)ii:::::::::

:~i,::>::;; :

:. ·:...

:. :·' ·' ._: .·..<:. _:: _::·'....:....:.::.:·.\~---··::_· /: <:.·.:

 

.......

_:;_:._:·.·.:..·...:.-.,.: ..\·/

 

time t = 0

x

,,----~-------~~~,~,

./'

,,

'

',

time t

l

p

 

 

I

 

 

\

 

'

 

.

 

\

 

\ ', B

 

 

 

,1\'

 

x '~~------------~~---~'

x

r.--11~

x 2' :1:2

e2

Figure .2•.1 Configuration and motion of a continuum body.

tions of continuum particles. Such a body is caIIed a continuum body, or just a continuum. A continuum is determined by macroscopic quantities. It has macroscopic dimensions that are .much "larger than the intermolecular spacings.

For the linkage of atomistic and continuum theories and for mixed atomistic/continuum computational schemes the reader is referred to

Configuration. Consider .a continuum body B with particle P E 8 which is embedded in the three-dimensional Euclidean space at a given .instant of time t, as indicated in Figure 2.1.

We introduce a reference frame of right-handed, rectangular coordinate axes at a fixed origin 0 with orthonormal basis vectors ea, a = 1, 2_, 3, according to properties introduced in ( l .18). As the continuum body B moves in space from one instant of

n which is
that in

58 2 Kinematics

time to another it occupies a continuous sequence of geometrical regions denoted by f201 ••• , n. Hence, every particle P of B corresponds to a so-cal led geometrical point having a position in regions n01 ••. , n. The .regions which are occupied by the continuum body Bat a given time-tare known as the configurations of Bat that time t. The continuum body B may have .infinitely many configurations in space.

The geometrical regions are determined uniquely at any instant of ti.me. Region n0 with the position .of a typical point X corresponds to a fixed reference time. The region is referred to as the fixed -reference (or undeformed) configuration of the body

B. A region at ·inithd time t = 0 is referred to as the initial configuration. We agree subsequently that the .initial configuration coincides with the reference configuration; hence, the reference time is at t = 0. Note dynamics the initial configuration is often not chosen as the referenc-e configuration. Now, the po.int X has the position of a particle occupied by P E B at t = 0, and P may he identified by the position vector (or referential posUion) X of point X relative to the fixed origin 0.

We now assume that the region 0 0 of space moves to a new region

occupied by the continuum body B at a subsequent time t > 0. The configuration of Batt is its so-called current (or deformed) configuration. We relate a typical point

X of the reference con.figuration to a point x of the current configuration occupied by a particle P EB at times t = 0 and t > 0, respectively. The :position vector (or current pos"ition) x serves as .a label for the .associated point x with .respect to the ·fixed origin 0. It is often convenient to call X po.int X, associated with the particle P E ·B at t = 0, and to call x point x (or place x)., associated _with P E B at :t.

The components of vectors X = )(AEA and x ·=:cu.ea are considered as being along the axes introduced. We .Jabel .1X";t, il = 1, .2, ·3, as the material (or referential) coordi~ nates -of point X and :r:a, a = 1, 2, 3, as the spatial (or current) coordinates of point x. ·we have assumed that the origins of the sets {EA} and {Ca} of basis vectors -coincide. ln addition, we have agreed to have the same reference frame for the reference and current configurations, which is why we set the basis {EA} identical to {ea}. We use just {ea} in the following.

To denote scalar, vector and tensor-quantities we use uppercase letters when they are evaluated in the reference configuration, and lowercase letters for corresponding quantities in the current configuration. Sometimes we employ the index zero for quantities acting in the reference configuration (for example, n0 , p0 , a0 ... ). This convention is often used in textbooks of continuum mechanics; however, it occasionally leads to con- "fticts with the ·general aim of Chapter l where we agreed to use lowercase, bold-face letters for vectors and uppercase, bold-face letters for second-order tensors.

For example. we use X for the position vector of a point corresponding to the reference configuration despite the convention of Chapter 1 that uppercase, bold-face Latin letters usually denote second-order tensors. A dear and unique notation .is difficult to establish, considering that certain quantities in continuum mechanics are characterized

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