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

Holpzaphel_-_Nonlinear-Solid-Mechanics-a-Contin

.pdf
Скачиваний:
225
Добавлен:
27.03.2016
Размер:
24.13 Mб
Скачать
a = 1, 2, 3. In
func-

6.2 Isotropic Hyperelastic Materials

 

219

eq. (6.1 l) we find that

 

 

 

 

 

 

 

9

.8\Jl(b)

_

aw(v)

_1

 

 

(6.42)

8b

-

av

v

'

 

~

 

 

so that (6.38) reads

 

 

 

 

 

 

 

_ J_1 8\Jl(v)

_. J-J. 8\ll(v)

,

(6.43)

u -

a

 

v--

v

8v

 

··v

 

 

 

 

which is another important stress relation characterizi~g the behavior of isotropic hyperelastic materials at finite strains. Note that, since v is the unique square root of b, aw /8b also com:mutes with v.

Constitutive equations -in terms of principal stretches. If the strain-energy

tion \JI is an invariant, we may regard \JI as a function of the principal stretches ,\a, the place of (6.27), we .may rep.resent Win the form

(6.44)

For the stress-free reference configuration the normalization condition (6.4) takes on

the form '11(1, 1, 1) = 0.

Consider the left stretch tensor v = b 112 describing the -defonned state of an isotropic hyperelastic material. .From the e-igenvalue problem (2. l l 6h we know that

,,\a denote the three principal stretches (the real eigenvalues) of v. Since the principal directions of v coincide with those of b (compare with eqs.. (2..116.h and (2. .l l7h)

they also coincide with the principal directions of the Cauchy stress tensor u (recall

representation (5.88)).

 

 

 

= 1, .2, 3,

Consequently, with respect to (6.43) the principal Cauchy stresses aa, .a

simply result in

 

aw

 

 

aa = J-

1

a = 1, 2, 3 ,

(6.45)

 

Ao a\ _,

"o.

with the volume ratio

(6.46)

according to (.5.9.lh.

In addition to (6.45), we introduce equivalent relations for the three principal PiolaKirchhoff stresses Pn and Sa, namely

aw

1

aw

a = 1, .2, 3

<6.41)

~, = a\ ,

s,l = ,

\ ,

,..,,1

An

''a

 

 

(compare with the fo.llowing Example 6.2)., which may be expressed in terms of the Cauchy stresses (6.45) as

a = 1, 2, 3 .. (6.48)

We start

220

6 Hyperelastic Materials

Constitutive relations (6.45) and (6.47) show that principal stresses in an isotropic hyperelastic material depend only upon the principal stretches. They are simply obtained by differentiating the strain-energy function with res.pect to the corresponding principal stretches.

-~,,~,,-~-·· ----·'""''''"'''•:"''''"'''''"'"'"'"''"'"'J..•''''''''''•·'-"•"''Y•'•','•. -';,,,,,,••.,•"o.. •,,,, . , -- ..•"'~':"V'~'"'"""___ ,.,,,,•.,,,....':Y'''''"'••-,••'"'.''-'"':':''~-·'--'''''"'"''":''''''''...'•-'•'..'•'''-.','-''•"''':•'"''"'•'''•""'•'•''•'•'•'•"'''''•''•'''''••'•'•'''''""''','•:

EXAMPLE 6.2 Consider the strain energy w(C) = w(A1, .1\2 , A3 ).

Obtain the con-

stitutive equations in the spectral fonns

 

 

:1

 

 

O' = L aaDa ®Ila ,

(6.49)

a=.l

 

 

3

:1

 

P == L Pciiia ® Na ,

S = L SaNa ® Na ,

(6.50)

a=l

n=l

 

where Ga and Pc.., Sa, a = 1, 2, 3, are the principal values of the Cauchy stress tensor er and the two Piola-Kirchhoff stress tensors P, .S according to the expressions (6.45) and (6..47)., respectively. The orthonormal vectors Na and Du = RNa, a = 1, 2, 3, denote t~e principal referential and spatial directions (axes qf stretch), respectively.

Thes.e constitutive equations describe isotropic response of byperelastic materials

and hold if and only if ;\1 '# ..\2 # .Ai1 # /\1.

Solution. with constitutive equation (6.50)2 which is expressed in tenns of second Pio.la-Kirchhoff stresses Sa, a = 1, 2, 3. We compute the derivative of the isotropic function w{C) with respect to the symmetric tensor C. By means of the chain rule and kinematic relation (2..123), we obtain for .the general case ~\1 # A2 # .1\~ =f )q,

(6.51.)

In (6.51 )2, .,,\~ are the eigenvalues (the squares of the principal stretches) and N0 the corresponding eigenvectors (principal .referential directions) of C (compare with the eigenvalue problem (2.115)). With (6.5 lh and the chain rule we find from (6.1.3)2 that

(.6.52)

which gives the desired .results (6.47h and (6.50h.

By the use of (6.52)2 , the relation according to ( 1.58.) and eq. (2..132) 1, i.e. FN0 = Aufi,u a = 1, 2, 3, the spectral fonn of .the first Piola-Kirchhoff stress tensor ·pmay be

6.2 Isotropic Hyperelastic Materials

221

found from transformation (3.67) as

(6.53)

Similarly, having in mind the results (6.53)4 and (2.132) 1, trans.formation (3.9) gives, using (3.10), the spectral form "'

<T = J

-·1

T

.J

-1

F

( ~ aw

,.

-- T)

~

-1. aw ,.

..

(6.54)

 

FP =

 

~ f);\

(na 0

Na)

= ~ J

;\a[))..

Ila 0

Ila

 

 

 

 

 

 

a='l a

 

 

lt=l~

 

 

 

of the Cauchy stress tensor, where the property 'F(na ®Na)'r

= (FNa) C8H1a, a= l, 2, 3,

was used (compare with relations which are analogous to (1.85).and (.1.58)).

__,,,...•..-:·-·... -·,.,,.....--·-·-·-,·~·----'~,····"'~ ....,,._,,,, ...;,,,:,,, ...;,,:,,,,,,:.,.,:,,,,:,,,,..,,,....,,.,, ...,,......_,,,.,,,,,...,,.,,..~,,,.,,,,,..,,,_.:.•..,...-,.-....,...-,,...._._,,..,.,,.,,,,., ... ~,,.,,,,,,. ..........,....._.....,,..,._....__, ____,___,,___,..__... _.,..- ..-,--...--,-·-,,... - .... ·.'~··· ..-:: ...--,,,-,,...

EXERCISES

l. By analogy with the procedure which led to (6.11 ), obtajn the eq. (6.42) and

relation

 

 

 

8\Jt(F))T =

2FT8W(b)

·

(6.55)

( aF

ab

 

2.Rewrite .the spectral representations of constitutive equations (6.49) and (6.50) for a given strain-energy function of the particular form w= '11(ln.A 1, ln...\ 2, ln,\3 ).

Consider the general case A1 # /\2 # /\3 # /\l •

3.Take the strain energy \JI as a function of the principal stretches characterizing the behavior of isotropic hyperelastic materials. Let at least one principal stretch be equal to the other.

(a)For the case in which we have two equal principal stretches, namely --\1 =

,\2 f. .A:1, obtain the constitutive equations

,

-1.,

aw(,..

,...

....

...... )

_1 , aw A

....

u = ·

.111

D>..i n1 0 n1

+ n2 0 n2

+ .1

;\3 a>..

03 0 n:i ,

 

 

 

 

 

 

 

 

3

 

aw

P = DAi

~

--

..

,.

aw "

A

,

(n1

0 N1

+ 02 0 N2)

+ a>.. o

3

0 N

 

 

 

 

3

3

 

 

 

 

 

 

 

1 aw ...

"

A

..

1 .aw ..

'"'

S = ---(N1 ·®Ni+

N2 ®

N2) + ---N3 0 Na .

.A1 8A1

 

 

 

/\3 8,\i

 

222

6

Hyperelastic Materials

 

(b)

Using the property (l .65)2 for the second-order unit tensor "I and relation

 

(2. l 22h for the rotation tensor R, show that for A1 =

,.,\2 = ~.\:i == A

 

 

3

 

 

 

 

u = a L Da 0 Da = al l

 

 

 

Cl= 1

 

 

 

:J

 

~i

 

 

 

S = SL N,1

 

 

P =PL·Da 0 Nu = PR ,

0 Na = SI ,

 

a="I

 

a=J

 

with the scalar-valued scalar functions er = J- 1;\8'11 / [)/\ and with P

8,T! IDA, s ·=,\-18\J! /D,\.

6.3 Incompressible Hyperelastic Materials

Numerous polymeric mate.rials can sustain finite strains without noticeable volume changes. Such types of material may be regarded as incompressible so that only .isochoric motions are possible. For many cases, this is a common idealization and accepted assumption often invoked in continuum and computational mechanics. In this sect.ion we present the constitutive foundation of incompressible hyperelastic materials.

1-ncompressible hyperelastidty. Materia]s which keep the volume constant throughout a motion are characterized by the incompressibility .constraint

J = 1 '

(6.56)

or hy some other equivalent expressions according to (2.1.77) (recall the expression (2.5 l) for Lhe volume ratio J). In general, a material which is subJected to an internal constraint, of which incompressibility is the most common, is referred to .as a .con- strained material.

In order to derive general constitutive equations for incompressible hyperelastic materials, we may postulate the strain-energy function

\JI= w(F) -

p(J - l)

,

(6.57)

where the strain energy \II is defined for J

= det"F

= 1.

The scalar JJ introduced .in

(6.57) serves as an indeterminate Lagrange .m.ultiplier, which can be identified .as a hydrostatic pressure. Note that the scalar v may only be determined from the equilibrium equations and the boundary conditions. It repres·ents a workless reaction to the kinematic constraint on the deformatio.n field.

= detC
= detb
Incompressible isotropic .hyperelasticity.

6.3 .Incompressible .Hyperelast.ic Materials

223

Differentiating eq. (6.57) with respect to the defonnation gradient F and using identity (2.174), we .arrive at a general constitutive equation for the first Piola-Kirchho.ff stress tensor P. Hence, eq. (6..1) may be adopted in the form

p = -11F-T + _ol_Jl_(F_)

(6.58)

.

DF .

 

An alternative derivation of (6.58) is obtained by reference to the express.ion (6.3)2 For :incompressible hyperelasticity, F is not arbitrary anymore and the e.xp.ressio.ns .in parentheses of (6.3)2 need not be zero. Howeverl (63}2 must be satisfied for every F which is governed by the incompressibility constraint in the form of J = :F-T : F = 0 (recall (2.177);.1). Consequently, adding the zero term to (6.3h, we find that

.

8\JJ (F)

+ pF

-T)

·

= 0 .

(6.59)

( .P -

DF

·

: F

·with standard arguments, the Coleman-Noll procedure implies physic.al expression

(6.58) .

.Multiplying eq. (6.58) by F- 1 .from the .left-hand side, we conclude from (3.·65h that the second Piola-Kirchhoff stress tensor Stakes on the form

(6.60)

where the inverse of relation (2.63), i.e. c-1 = F- l F-T, and identity (6.11) are to .be used.

However, multiplying eq. (6..58) by FT from the right-hand side, we conclude from (3.9) that the symmetric Cauchy stress tensor <T may be expressed as

(J' = -1JI +

8\J!(F)

n

.

DF

F:I

= -JJI _L F

.

 

r .

(ow(F))T

(6.61)

DF

 

The :fundamental constitutive equations (6.58), (6.·60) and (6.61) are the most general forms used to define .incompressible hyperelastic materials at finite strains. Equations (6.60h and (6..61) are associated with (5.97) and (5.96).. Note that the response functions SJ(C) and g(F) occurring in (5.97) and (.5.96) 1 are identified by jj(C) =

2D\J!(C)/8C and g(F) = (O\Jl(F)/8.F)FT = F(8'1r(.F)/8Ffr,i.e. constitutive equation (6.2h for J = 1.

For the case ·of isotropy we have already pointed out that the dependence of 'Iton the Cauchy-Green tensors C orb may be ex.. pressed by their three strain invariants (see eq. (6.27)). However, for the incompressible case we consider the kine.matic constraint, namely .la = l. Therefore, the two principal invariants / 1 and 12 are the only independent deformation variables.

224 6 Hyperelastic Materials

For a review on the theory of incompressible isotropic hyperelasticity see, for -example,

OGDEN [l.982, 1986].

A suitable strain-energy function for incompressible isotropic hyperelastic materi-

als is, in view of (6.27), given by

'1t = W[I1(C),l2(C)] - ~p(Ia -1) = W[I1(b)J2(b)] - ~p(h -1) , (6.62)

where p/2 serves as an indeterminate Lagrange multiplier.

In order to examine the associated constitutive equation in terms of the two prin-

cipal strain invariants

/ 1., / 2 .,

we proceed by deriving (6.62)i with respect to

tensor

C. Analogous to the procedure which led to (6.32) we

 

find,

using the chain rule,

eqs. (6.30)J, (6.31) and the constraint ! 3 = 1, that

 

 

 

 

 

 

 

s = . aw(I1.,I2) _ a"{p{J3 -1)] = _ c-1

'>(aw

 

 

aw)

 

_9awc

<

 

9

 

 

 

1

 

1

-ar2 '

663

... ac

ac

P

+ - a11 +

 

 

1ar2

 

 

· )

which is basically constitutive equation (6.32), in which the term / 3{ 8\J! / 813 )

is sub-

stituted by -p/2.

 

 

 

 

 

 

 

 

 

 

A push-forward operation of (6.63h and an elimination of b 2 in favor of b- 1 (see relation (5.93)) yields two alternative forms of u_, corresponding to eqs. (6.34) and

(6.35), namely

u = -pl+ 2

a\J!

+ 11

aw )

b - 2

aw "

 

( 01

DI2

012 b~ '

(6.64)

 

1

 

 

 

 

 

.

aw

aw

-1

(6~65)

,,,. - -pl

+ ? - b - ? - b

 

.v - .

~ 811

- 8I2

 

 

Note that the scalars pin eqs. (6.64) and (6.65) differ by the term 212 (8\f! /812 ). By· comparing (6.64) and (6.65) with (5.98) we obtain explicitly the response coefficients

 

aw

aw )

a1 = 2 ( O!i + f 18!2

f3

_ 9 aw

'

1 - ""'8/1

 

aw.

'

a2 = -2 8I2 '

/J _

_ 9 aw

P-l -

*-'8/2 .

(.6.66)

(6..67)

In order to find a constitutive equation for incompressible materials which is associ- ated with (6.38), we recall the transformation (6.55). Then, (6.61) gives the-constitutive equation

 

8\J!(b)

 

8\ll(b)

 

.u =-pl+ 2 Db

b =-pl+ 2b ab ·

 

 

 

.

aw

(6.68)

or

 

 

aab = -pc)ub + 2bac -al

 

Jc(J

6.3 Incompressible Hyperelastic Materials

225

in tenns of the spatial strain variable b. This is only valid for .incompressible .isotropic hyperelastic materials.

If we express \JI as .a function of the .three principal stretches Aa we write \JI = '11("\11 A2 , A3 ) - p(J - 1) in the place of (6.57), with the-indeterminate Lagrange mul- tiplier JJ. Using 8.JIoAa = J)_-;;1_, a = 1, 2, 3, which is relation (2 ..1.74) expressed in principal stretches, eqs. (6.45) and (6.47) are then replaced by

 

 

aw

a -

1

?.

T3

 

(6.69)

 

0-ci

= -p + >.a{),\a '

'

 

 

7 _.,

 

 

 

1

{)\J!

,

 

 

 

 

a = 1, 2, 3 ,

 

Pa = - \JJ + a,

 

 

 

 

(6.70)

l'\a

.1\a

 

 

 

 

 

 

 

with the three principal Cauchy stresses 0:0

and the Piola-Kirchhoff stresses

Pa, Sa.

These stress relations -incorporate the unknown scalar p, which .must be determined from the equilibrium equations and the boundary conditions. The incompressibility constraint J = 1 takes on the form

(6.71)

.leaving two independent stretches as the deformation measures. Expressing the first and second Piola-Kirchhoff stresses in terms of the Cauchy stresses (6.69), we obtain, by analogy with (6.48), ~1 = A;1cru and 80 = /\;2cra, a = 1, 21 3.

EXAMPLE 6.3 Consider a thin sheet of .incompressible hyperelastic material which is embedded in a reference frame of (.right-handed) coordinate axes with a fixed set of orthononnal basis vectors ea, a = 1, .2, 3. Suppose that the .axes are aligned with the major faces of the sheet.

A deformation created by the stretch ratios ,.\ 1, /\2 along the directions -e1, e2 results in a (lw11u~ge11eous) biaxial deformation with the kinematic .relation (2.130). The assodated stress state is assumed to be plane throughout the sheet so that the Cauchy stress

·components a1:~, 0"231 a 33 are equal to zero which is in accordance with (3.59).

Show that the biaxial stress state of the homogeneous problem is of the form

(6.72)

(6.73)

(see RIVLIN [1948, eq. {6.5)]), with the principal .invariants 11 = ~W + ..-\~ + .1\)2,,\22 ,

I2 ·=,,\i A~ + /\12 + /\i··2 and l:i = 1.

226•

6. Hyperelastic J\tlaterials

 

Solution.

Since the tensors u and ·b are coaxial for :isotropic elastic materials (recall

p.. 20 l), the principal stresses follow from (6.65),

 

 

a == 1, 2, 3 ,

(6.74)

where ,\~ are the three eigenvalues of the left Cauchy-Green tensor b (see the eigenvalue problem (2. l l 7h,). This relation was ·first presented by RIVLIN fl 948].

With the condition of incompressibility (6.7.l) in the form of A;1 = ()q.1\2 ) -1 and

the boundary condition a~J

= 0 we may determine p explicitly. For a = 3~

we deduce

from (6.74) that

 

 

.(73 ·=0

--+

(6.75)

This result substituted back into eq. (6.74) leads to the nonzero stress components a 1 and a 2 • •

EXERCISES

l.Consider a thin sheet of incompressible hyperefast.ic materia.l (.Ta = 1) with the same setting as formulated in Example 6.3.

(~)Consider a simple tension for w·hkh Ai = .A. Then, obeying incompress-·

ibility constraint .1\ 1,.\~/\3 = 1, the equal stretch ratios in the transverse di- rections are, by symmetry, A2 = .1\3 = .1\- 1/ 2 Show that for this .mode of deformation the homogeneous stress state reduces to a 1 = a, a 2 = aa = 0,

with

w"here the invariants are / 1 = 2.1\- 1 + A2 , / 2 = /\-2 + 2.A.

As a special case of the biaxial deformation, as discussed in Example 6.3, consider an equibia.xial deformation for which ,,\1 = ;\2 = ,,\, ,\i = ,,\ - 2 and a1 .== -a2 = a, a:i = 0. Show that

W 1t

h I

I = _.)\- + A

, J 2 = A

+ .....;\ -.

.

? .., \ -•1

T \"

? ., -'>

A .material

6.4 Compressible Hypcrc)astic ·rv1.aterials

227

(b)Consider a ho.mogeneous pure shear deformation with the kinematic relation ,,\ 1 = ,\, --\2 = 1, Ai = l/ A. (compare with eq. (2.131 )). Show that the nonzero Cauchy stress components are

(6.76)

(6. 77)

. I I

l =

J.

=

\2

\ _ •.,

+ 1.

W:lt 1

2

.1\

+ .,,, ~

2.Consider a th.in sheet of incompressible hyperelastic material with the same setting as formulated in Example 6.3 but subjected to a homogeneous simple shear deformation which is caused by a motion in the form of (2.3) (compare also with

Exercis-e 2 -on p.. 93).

(a)Show that the associated stress state is completely defined by

') aw

a,T!

~

au = -p + 2(1 + c~) BI

- 2·aI6

1

. 2

 

(6. 78)

where Jt > 0, called the shear .modulus, is a measure of resistance to distortion and p is a scalar to be determined from the boundary conditions.

(b) Consider a plane stress state throughout the sheet in the sense that the face of the body nonnal to the direction e;1 is free of surface tractions, i.e. a 1;i = o2a = a3:~ =O. Show that the nonzero Cauchy .stress components ure

0-.12 = JJ.C

6.4 Compressible Hyperelastic Materials

which can undergo changes of volume is said to be compressible. Foamed elastomers~ for example, .are able to sustain finite strains with volume changes. The only restriction on this class of materials is that the vo]ume ratio J must be positive.

Compressible hyperelasti-dty.

228

6 .Hyperelastic Materials

In this section we introduce suitable constitutive equations .in order to characterize compressible hyperelastic materials, and we discuss isotropy as a special case.

Since some materials behave .quite differently in bulk and shear it is mos.t beneficial to split the deformation locaUy into a so-called volumetric part and an isochoric part, origina11y proposed by FLORY [1961] and successfully applied within the context of isothermal finite strain elasticity by, for example, LUBL.INER [1.985], SIMO and TAYLOR [199.la], OGDEN [1997] and within the context ·of finite strain elastoplasticity by, for ex.ample, SIMO et al. [1985] amon.g many others.

In .particular, we consider the deformation gradient F .and the corresponding strain measure C = FTF. Rather than dealing directly with F and C we perform a multiplicative decomposition ofF into volume-changing (dilational) and volume-preserving (distortiona/) parts, often used in elastopl.asticity (see, for example, LEE [1969]). We write

(6.79.)

The terms J 113 I and J2!=11 are associated with volume-changing defomiations, while F and C = F1~F are associated with volume-preserving deformations of the materia'I,

with

 

 

 

 

and

detC = (detF) 2 = 1 ,

(6.80)

where

 

 

 

\ -

J-1/:J \

.a= 1, 2, 3

(6.8 I)

''a -

A.u ,

characterize the so-called modified principal stretches. We call F and C the ·modified deformation grad.font and the modified right Cauchy-Green tensor, respectively. A material for which dilational changes require a much higher exterior work than vo.lumepreserv.iQg changes is called a nearly .incompressible (or sUghtly compressible) material, for which the compressibility effects are small.

The concept of the multiplicative decomposition of F is supported additionally by the field of computational mechanics. For example, to avoid numerical complications in the fi~ite element analysis of slightly compressible materials it is often advantageous to separat'e numerical treatments of the volumetric and isoc.hor.ic pa.rts of the deformation gradient F; this will be ·discussed in Sections .8.5 and 8.6.

Before proceeding to exa.inine constitutive equations for compressible .hyperelastic materials it is first necessary to stick to kinematics and to compute the derivative of the modified right Cauchy-Green tensor C relative to the symmetric tensor C. By means of (5.91.)2, we obtain from (6.31 }2, 8J2 Iac = J2 c-). Using the cha.in .rule we arrive

at

(6.82)

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