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

Holpzaphel_-_Nonlinear-Solid-Mechanics-a-Contin

.pdf
Скачиваний:
189
Добавлен:
27.03.2016
Размер:
24.13 Mб
Скачать
as .an eigenvalue.
f/a..
:-..:· 3..2
..\·:Normal

3.2 Extremal Stress Values

119

Extremal Stress Values

and shear stresses vary in magnitude, direction and location. The designer of :\·:-.astructure must prove that a certain material will not fail as a result of these stresses. \Magnitudes, directions and Jocations of extremal stress values must be identified. We ·_(_:~:p~~sent briefly the standard results on extremal stress values used in practical engineer-

::/_J;?.g~..

?..-,M~ximum .and minimum normal stresses. Consider the Cauchy traction vector

.:f:·t(~;t_, n) = u{x, t)n on an arbitrary oriented surface element throu_gh any po.int x E n

\\Aiftime t. First of all, we wisl1 to find the unit vector n at x, indicating the direction of

.<."·:;11a.ximwn and minimum values of normal stresses a.

/:! In order to obtain these maximum and minimum (or so-called extremal) values of ::~}-.//.~~-may app1y the Lag-range-multiplier method and claim the stationary position of

f~nctional .C to be

~>t.:/·:~·::·,...::

.C (n, ,\*)

= n ·un -

,,\*(InI

 

1)

,

 

 

 

 

 

(3.32)

 

 

 

 

 

2 -

 

 

 

··:...·).:=:":.:·:···~ :· ..

 

 

 

 

 

 

 

 

 

 

..:_,:-(or, .in index notation,

 

 

 

 

 

 

 

 

 

:::.-.::::.-.=:·.·.:: .

 

 

 

 

 

 

 

 

 

 

./·):}:.,-'.:::<·:··:

£(na' ;\*)

= naacJbnb -

,..\* (n,,n(l -

1)

.

(3.33)

;):'~e~~.lnl2 -

1 = 0 (n,.na -

1 = 0) characterizes the constraint (auxiliary) condition

':i(:~rid ,,\* the La.grange multi.plier.

 

 

 

 

 

 

 

:iit}·.::::;:,_ ~or the stationary position of£, the derivatives 8£/8n0

and 8£/8.,\* ·must vanish,

\=/T~~~'~·,"w.ith resp.ect .to (3.33),

 

 

 

 

 

 

 

 

 

":··:·=:..:···..... ·

 

 

 

 

 

 

 

 

 

 

 

i)t:,

= O"a11(c~bcnn+ 1Sacnb) -

A*(2na811c)

 

 

anc

= 2(aca·nci

- ,,\*nc) = 0

.,

 

 

 

(3.34)

 

 

 

 

 

 

= nllna -

1 = 0

 

 

 

 

 

(3.35)

 

a,..\*

 

 

 

 

 

i\W~'haVeused 8na/an,,

= Oat.. property (l .22h and the fact that the stress tensor (T is

(:fs·§frifuetric. In symbolic notation (3.34h and (3.35) read as

 

::,:t;:~{ <

o-n - ,,·n = o ,

lnl = 1 .

(3.36)

fi!:)')\O~ince this is an eigenvalue problem the general results and rules of Section 1.4

;/(hold. The Lagrange multiplier ,,\* .may be identified Adopting the

\:ff~tevant notation of Section 1.4, we must solve homogeneous algebraic equations for

{:(J~~".tinknown eigenvalues,..\:, a = 1, 2, 3, and the unknown eigenvectors ila, a, = 1, 2, 3,

(:)fo'·:the form

;;:~:.:}~:v:.·r./~-;,_-:·_:·:... : ..

:!:l:/:':'.l>:_\·:_::(:':...) .....

: : ..

(a = 1_, 2, 3; no summation)

(3..37)

1, 2, 3.

120 3 The Conce.pt of Stress

(compare with eq. (l.166)). In order to obtain the eigenvalues we must solve the characteristic polynomial-of u, which is cubic in 1\:., i.e.

(3.38)

with the three principal stress invariants lu of tensor u (com.pare with eqs. (1.170)- ( 1. l72)), Le.

f..1

(u) = trcr = .,\~

+ ,,\~ + /\;

,

 

(3.39)

J2

{u) = ~-

[{tru)

2

-

tr(cr

2)]

= /\7A;+ A7/\; + ,\;,\;

,

(3.40)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(3.4.1)

The three eigenvalues ,,\~ are here called principal normal stresses and typically are denoted by a,1' a = The principal values ·'1a include both the maximum and

minimum normal stresses among all planes passing through a given x.

The corresponding three orthonormal e.igenvectors iia, which are then characterized through refation (3.37), are the principal directions of u. The normal stress is stationary along these principal directions. Their related norma] planes are known as principal ,planes. Since the stress tensor u is sym·metric the set of eigenvectors form a mutually orthogonal basis. Shear stress components vanish at normal (principal) planes, since obviously aH eigenvectors are normal to their respective principal planes.

We -may rep.resent u in the spectral form

3

 

 

u(x_, t) = L O"utln ® fiu

(3.42)

<t-=l

 

 

(see Section 1.4), satisfying the eigenvalue problem uita

= alliia, a = 1, 2, 3.

Note

that only for isotropic materials do the introduced three principal directions fia

of u

coincide with the principal directions, as defined in Section 2.6..For a brief explanation the reader is referred to Section 5.4, in particular to ·eq. (5.88), and also to Section 6.2.

Maximum .and minimum shear stresses. The next goal is to find the direction of the unit vector n at x that gives the maximum and minimum values of shear stresses·T.

In the following we choose the eigenvectors fia, a = 1, 21 3, of u as a possible set of basis vectors. Then, according to the spectral decomposition (3..42), all non-diagonal· components of the matrix [u] vanish. The components cr1 ~ a 2 , a3 form the only entries of the stress matrixT and [u] is 'diagonalized'.

The Cauchy traction vector ·tp then appears in the simple form

(3.43).:

where n = n 1n1 + n2 n2 + n:-ifi3 denotes the .unit vector normal to an arbitrary oriented_·:_ surface element. The normal stress .and shear stress on that surface element follow from:·.·,

 

 

3.2 Extremal Stress Values

121

eqs. {3.26)i, (3.29) and have the forms

 

 

-

2

2

2

(3.44)

CJ -

n • t0

= -a1 n 1

+ a2n2

+ a3.n3

T2 = lt1112 -

a2

 

 

 

 

 

 

 

 

(3.45)

:·./ where the relations n = naila and (3AJh with the property nu ·ob= 6ab (according to

:). eq. (1.21)) are to be used.

· With the constraint (auxiliary) condition lnl2 - 1 = 0 (~ana - 1 = 0) we are now

.·able to eliminate n 3 from eq. (3.45h. Since the principal stresses a1 , o-2 , a:i are known,

..:._.r2 is a function of n1 and n2 only. In order to obtain the extremal values of r2 we

>differentiate T 2 with respect to n 1 and n.2 and identify with zero:

:-:·.:·.

a 2

 

;c\

{}:I = 2n1 (a1 - a3) { a1 - O";i - 2((a1 - a3)ni + (a2 - l13)nm = 0 ,

(3.4.6)

..

 

 

 

 

 

 

:·.· ..

.:\·.>_::.:. Yet we must find the solutions for the stationary shear stress directions, i.e. tor the

{:\:complete set of principal directions. The first (obvious) solution of (3.46) is obtained

\'./::,.by taking n 1 = n 2

= 0 and n 3 = ±1. If we consider further the two cases n1.

= n3 = 0

({(~Jlcl n2 =

na =

0 as the associated component pairs, we find

that n2 =

±1 and

r/)1·} = ±1,

respectively. In summary, with n = nafia we obtain the three respective unit

:</:·..normal vectors which are

 

 

 

 

~.\\·.<.':··,-.:..

 

 

 

 

 

 

:·::.· .

 

D = ±03

·D = ±D·>

D = ±01

,

(3..47)

. :.::-. ·.

 

"ifi~§frespondingto the set {00 }

of principal dir:ct:ons of u. Obviously, this first solution

!/\):d~scribes

the principal stress planes upon which T

is minimized (in fact r

= 0), as

r//~g:rt:M1y concluded from <3.45).

/)\<._~,::.\\._Thesecond solution of (3.46) may be found by assuming n 1 = O_, which gives, from :_/;.)(~3_"•.46h, the component n2 = ±l/v'2. The condition lnl = I leads to n:1 = ±I/J2~

::\:~N:_):i.Oing so for n2

= 0 we find n 1

= n3 = ±1/./2, while .for ·na = 0 we obtain the

'.:~{};B_Q"inponents n 1

= n 2 = ± 1/./2.

These results substituted back into (3 .45) give the

&\{Wssociated extremal values r 2In summary,

:·... ·....

::_:(:.:·:·.·:· .... ·.·

·.. ·:·:=.:::.

:.~:/("=..<····. ·.

..

·: :····.·.

.... :: ..

:·.=.-·:.-·:·.·:· .

...........·.:.. ·:·

\W.V\\'.::.._:: .. :

: ..":.=.: -~=. : : .

··.·.::.::.· ..

:.:·::.":::·.·.

·.-·::·.".":::: ..

1

±

1

2

1

')

(3.48)

n = ±-n-,

- 03

T

= 4(£12 -

0"3)~

Vi-·

v'2

 

 

 

 

1

 

1

')

1

')

(3.49)

n = ± rnn1 ±

1n·na

T-

= 4(l1t -

£13)-

v2 ·

 

v2

 

 

 

 

1

 

1

 

 

 

(3.50)

n = ±-01 ±-n2

 

 

 

/2

 

v'2

 

 

 

 

where
·.rz: 1l
:r:2 , :1:3 ,
+ O"minl/2.

122

3 The Concept of Stress

Consequently, the maximum magnitude of the shear stress denoted by Tmax by the largest of the three values of (3.48)2-(3..50)2 . We obtain

1

Tmax = 2lamax - t.Tmin·j '

is given

(3.5 l)

where Gmax and D"min denote the maximum and the minimum magnitudes of principal stresses_, respectively. It is important to note that the maximum shear stress acts on a plane that is shifted about an angle of ±45° to the principal plane in which the maximum and .minimum principal stresses act (compare (3.47) with (3.48}1-(3.50h),.

In addition, it can be shown that the normal stress -a which is associated with Tmax

has.the value.a= lamax

EXERClS.ES

l. A Cauchy stress tensor u, whose components a(j1, depend on the coordinates is given in the form of its matrix representation

0

0

-/3:1:a

.a and {3 are constants. For a certain point x, given by its coordinates

(0, /32 , n), find

(a)the three principal stress invariants of tensor u ..

(b)Compute the principal stress components and the associated principal directions of stress. Verify that the principal directions of stress are mutually orthogonal.

(c)Compute the maximum magnitude of shear stress and the plane on which it acts..

2.The Cauchy stress matrix is given as

[u] =

7

0

14 ·i

 

 

0

8

kN/cm

2

0 ·

[

14

0

-4

 

 

:

 

 

(a) Compute the principal stress components .and the associated principal di-

rections.

(b) Compute the maximum shear stress and the plane on which this maximum shear stress acts.

3.3 Examples or States of Stress

123

3.At a given current position x E Q consider a plane whose unit normal n makes equal angles with each of the three principal directions .Da of stress. Such a plane is typically called an octahedra·I plane characterized by the unit normal vector

n = naila, <L = 1, 2, 3, with components n 1 = n 2 = n:-1 = .1/.../3.

(a)With reference to eqs. (3.26)1 and (3.29) verify that the normal stress and the shear stress on the octahedral plane are

 

 

 

 

 

 

I

 

 

 

+ <l2 + 0"3)

 

 

 

 

 

O'oct

= 3{O"t

,

 

 

Toct

1[

(a11

-

a22

)2

+

(

 

 

 

•)

+ (aa:1 -

')

 

= 3

 

 

an - aa:1)~

u1 i)~

 

 

 

2

 

2

 

+ a:n)

2

1./.2

 

 

 

 

 

 

+6(a12

+ a2a

 

J

 

 

 

 

 

 

1 [

 

0'2)

2

 

+

(

02 -

o;i}

2

+ (a:i - 0"1)

211/2

,

 

= 3 (01 -

 

 

 

 

 

where a 1, a 2 ,

u;1 are the principal stresses. The normal and shear stress on

the octahedral plane are typically denoted by O"oct and T0 ct,, respectively.

(b)Find the traction vector, the normal and shear components on the octahedral plane for .a certain point at which the principal stresses are a.1 = 2 kN/cm2_.,

a2 = a, a3 = 11 kN/cm2, where n is a constant. Detennine a so that Toct is the maximum shear·stress..

/ 3.3 Examples of States of Stress

<·.. Since ~1 Cauchy traction vector at given (x, t) .is defined for each surface element repre-

:/·sented by the outward unit .normal n, we may find infinitely many traction vectors. In

!\ .this section we show some important cases.

::::·.. ·State of stress. The set of pairs { (t, n)} at a given po.int is ca11ed the s.tate of stress, ·:-··· ~ompletely .determined by the stress tensor u at this given po.int and time t according <;-- ~? Cauchy's stress theorem. Since t .is a linear combination of n in three dimensions, :/...... t_~ree linear independent pairs (t, n) form a complete basis. The state of stress is said ·:.··...·_to be homoge:neous if the stress tensor does not depend on space coordinates at each

·time t.

·We now _give some important examples for different states of stress.

/·_··. (i) A pure normal stress state at a certain point is ·given by the stress tensor

 

u = a(n~n)

.or

(3.52)

·(see Figure 3.5(a)).

If a 12
pure shear stress
uniform compression

124

3 The Concept of Stress

Post-multiplying (3.52) with the unit vector n we find using rule (1.53) and (3.28)2 that

un = a(n ® n)n

=a (n · n) n =an= tll .

(3.53)

~

n

 

1

 

 

Evidently, tll is along (or opposite) n (compare with Figure 3.4). This stress state is characterized by a normal stress a, the shear stress is zero. The stress a characterizes either pure tension (if a > 0) or pure compression (if a < 0).

Consider a rectangular Cartesian coordinate system. If a 11 = er = const and all other stress components are identically zero, then we have either uniform tension or in the x 1-direction. In the literature these stress states are often referred to as uniaxial tension and uniaxial compression, respectively. This may

be imagined as the stress in a rod (with uniform cross-section) generated by forces applied to its plane ends in the x 1-direction, which is one principal direction of stress.

A uniform tension or compression of a uniform cylindrical rod leads to a deformation, which we had already called uniform extension or compression (see p. 92).

A three-dimensional stress state in which all shear stress components vanish, i.e. aab = 0, a i= b, is said to be triaxial, while in a biaxial stress state (which is associated with biaxial deformations) we have a pair of non-vanishing normal stresses, i.e. a1, a2 ; a2 , a3 ; a 1 , a3 However, in an equibiaxial stress state (which is associated with equibiaxial deformations) the non-vanishing normal stresses have the same value.

(ii) A pure shear stress state at a certain point is given by the stress tensor

u = r(n ® m + m ® n)

or

(3.54)

(see Figure 3.5(b)). Here, the stress r

is related to the directions of the unit vectors n

and m, with property m · n = 0.

 

 

Post-multiplying (3.54) with the unit vector n we obtain with rule (1.53) and (3.28h

un = r(n®m+m®n)n

 

= r((m · n) n + (n · n) m) =rm= t;. .

(3.55)

~~

 

0

I

 

Evidently, t~ is tangential to the surface and along (or opposite) m, i.e. perpendicular to n (compare with Figure 3.4). This stress state is characterized by a

T (or'" pure tangential stress) and the normal stress is zero.

Consider a rectangular Cartesian coordinate system. = a~n = T = const and

all other stress components are identically zero, then this stress state is characterized by a uniform shear stress. A state of uniform shear stress leads to a uniform (or simple) shear deformation, as illustrated in Figure 2.2 (see also Exercise 2 on p. 93 ).

3.3

Examples of States of Stress

125

(a)

(b)

t~ =Tm

 

lt~I = T

n

n

 

tll =an

(c)

n

tU =-po

Figure 3.5 Examples of stress states.

(iii) A hydrostatic stress state at a certain point is given by the stress tensor

u =-vi or (3.56)

(see Figure 3.5(c)). Within a rectangular Cartesian coordinate system that means that

we have normal stresses a 11

= a 22 = a33 = -p and no shear stresses on any plane

containing this point, i.e. a 12

= a~a = a:11 = 0.

 

Post-multiplying (3.56) with the unit vector n, we obtain, using (3.28)2,

un = -(vl)n = -vn = t~! .

(3.57)

Evidently, tU is opposite to the normal vector n. The relationship un =

-pn holds only

if n is a principal direction of u. Note that any three mutually orthogonal directions may be regarded as principal directions (compare also with eq. (1.179)). This stress state is characterized by a scalar p, known as the hydrostatic pressure. In ge~eral, the hydrostatic pressure is a scalar function of time t, and in the literature is often introduced with the opposite sign.

One example of stress state (3.56) is in an (elastic) fluid (without motion) that is

not able to sustain shear stresses. Taking the trace of (3.56), we obtain

p = -

1

 

1

 

-tru

or

p = - 3a"u '

(3.58)

 

3

 

 

 

126

3

The Concept of Stress

.indicating that the pressure of the (elastic) ·fluid is a mean pressure.

The three stress fields described for a certain point x at .t (with u constant) corre-

spond to the states, as illustrated in Fi_gure 3.5.

(iv) A plane stress state at a certain point is given by the relation

 

 

(359)

and a 11 , a 22 _, a 12

are functions of the coordinates :z;1 .and :r.2 only. Consequently, the

~r~1-direction, as represented by it3,

is a principal direction of stress with a zero car..

responding principal stress a~i:1The other two principal directions acting in a pfane normal .to ila are inclined at an angle fJ with the :r; 1 and x 2 direction, where

 

 

2a12

 

 

(3.60)

 

tan2B = ----

 

 

 

 

a1.1 -

0"22

 

 

 

The corresponding maximum and minimum -stresses are given by

 

1

+ 0"22) ±

[1

?

')

] J j'>

(3.61)

2(0"11

4(0"11 -

0"22)- + <J"j2

- 1

defining a biax.ial stress state.

The maximum shear stress Tmax for a plane stress state will be the largest of the three values of (3.48}~-(3.50)2 s-ince cr:1 = 0. The planes of extremal shear stresses form angles of ±45° with the .planes of the principal stresses.

A plane -stress state occurs at any unloaded surface in a continuum body and is of practical interest.

EXERCISES

I. For the case of plane stress, show that -e-q. (3.38) reduces to eq. (3.61 ).

2~ Assume a plane stress state in a rectangular cube bounded by the planes :r:1

±a, :I:2 = ±b, :1;3 = ±c .. The state of stress at a po.int with coordinates ;1~1 , ~r:2, :1::1 in the cube is given by

where n, /3 are constants. For a point with-coordinates (a/2, -b/2,-0), determ:ine.

(a)the principal normal .stresses and the associated principal directions,

(b)the planes_, characterized by the unit normal n_, that give the maximum and minimum shear stresses and the magnitude of the extremal shear stress.

(c)Find the total Cauchy traction vector on each face of this rectangular cube.

3.4 AHernative Stress Tensors

127

3.4 Alternative Stress Tensors

Numerous definitions and names of stress tensors have been proposed in the literature. Each definition has advantages .and disadvantages. In the following we discuss stress tensors used for practical nonlinear analyses ..Most of their components do not have a direct physical interpretation.

Often it is convenient to work with the so-called Kirchhoff stress tensor r, which differs from the Cauchy stress tensor by the volume ratio J. It is a contravariant spatial tensor field parameterized by spatial coordinates, and is defined by

T =Ju or Tut> = ,Ja ab (3.62)

We introduce further the second Piola-Kirchhoff stress tensor S which does not admit a physical .interpretation in terms of surface tractions. The contravariant material tensor field is symmetric .and parameterized by material coordinates. Therefore, .it often represents a very useful stress measure in computational mechanics and in the formu- lation of constitutive equations, in particular, for solids., as we will see in Chapter 6.

. The second Piola-Kirchhoff stress tensor is obtained by the pull-back operation on the contravariant spatial tensor field T~ by the motion x, which is, according to (2.85h,

or

1 -lp-l

(3.63)

SA H = FAu 1Jb Tub

Hence, the Kirchhoff stress tensor is the push-forward of S, i.e.~ using (2.85}1,

 

or

 

(3.64)

Using eqs. (3.6.3)2, (3.62) and (3.8) we obtain the Piola transformation relating the two stress fields S and u, i.e..

S == .IF-tuF-T = F- 1p =ST

or

S

 

JF-t

- tF - l

 

 

p-lp

=

c

(3.65)

 

 

 

 

 

 

A 0

= ~ . An · ~I.lb

fTub

= ·Aa

n IJ

•.J II A

 

with its inverse,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

or

1-'ub -

'

 

uA .t'.bJJ' A.IJ

(3.66)

 

 

 

 

 

rr

-

 

7- l F.

.'fi'

S .

 

From eq. (3.65) we find a fundamental relationship between the first Piola-Kirchhoff stress tensor P introduced in (3.8) and the symmetric second Piola-Kirchhoff stress tensor S, i.e.

P= FS or

(3.67)

128 3 ·The Concept of Stress

In addition to the four stress tensors -given above, we .introduce another important -quantity. .It ·is a material stress tensor formally defined as

Tn = .RTP or (3.68)

The non-symmetric tensor Tu, which -is not in general positive definite, is known as the Biot .stress tensor. With (3.67) and the polar decomposition F = RU we deduce from (3.68) that Tn = RT (FS) = US. Herein, R and U denote the rotation tensor (with defR = 1) and the (positive definite) symmetric right stretch tensor, respectively.

They are according to the polar decomposition of the deformation gradient F (see Sec-

tion 2~6).

Multiplying eq. (3.68) by R from the left-hand side, we obtain the polar decompo-

s.ition

 

 

P=RTn

or

(3.6.9)

for the first Piola-Kirchhoff stress tensor, which is in analogy with that for F. Since the .Biot stress tensor Tn is not positive definite this decomposition .is not unique, in general.

Other examples of stress tensors are the symmetric so..called corotated Cauchy stress tensor uu, as introduced by Green and Naghdi, and the Mande·J stress tensor

~ which is in general not symmetric. These tensors are defined with respect to an intermediate configuration. In .order to obtain the corotated Cauchy stress tensor take relation (3.66) .and use the symmetric right stretch tensor U instead of F. Then apply relation (3.65}1 and the polar decomposition -F = RU to find

 

(3.70)

The .Mande] stress tensor is defined to be

 

E =CS

(3.71)

which is often used to describe .inelastic (plastic) materials.

EXERCISES

l.Consider an infinitesimal resultant (pseudo) force dfn = (TnN)dS, where TnN denotes the Biot traction vec.tor. Using (3.68) and the linear transformations

(3.3)2 and (3 ..1)2, verify that

= . 'r dfn R df ,

which shows that d.fn -only differs from elf (.introduced in (3 ..l )) by the rotation tensor RT.

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