Holpzaphel_-_Nonlinear-Solid-Mechanics-a-Contin
.pdf
2.5 Strain Tensors |
79 |
··.: ·Consequently, given the nine components Fa.4 , it is easy to compute the six components
·. |
CAB = Cn.-t via (2.63), but given CAH it is impossible to compute the nine components |
|
·: |
~1.AA With definition (2.63) and eqs.. (l. l.01.) and (2.51.) we find that |
|
|
detC = (detF)2 .= J 2 > O . |
(2.65) |
··:·The so-called Piola deformation tensor, denoted by B, is defined by the inverse of the
·..dghtCauchy~Oreen tensor, i.e. "B = c-1, with c-1 = (FTF)- 1 ==
As a further strain measure we define the change in the squared lengths, i.e.. With (2.62)a, the use of the unit tensor I and eq. (2.57)2 we have
E = -(F· F- I) |
or |
E~w = ~(Fa,i.Fi,11 - i5A11) , |
|
1 |
T |
|
- |
2 |
|
|
|
|
|
|
|
(/\de)2-
(2.66)
(2.·67)
.·.· where the introduced normalization factor 1/2 wiH be evident within the linear theory. ·..···: ·This expression describes a strain measure .in the direction of a0 at point X E .n0 . In _:_.· (2.67) we have introduced the commonly used strain tensor E, which is known as the · ~reen-Lagrange strain .tensor. Since I and C ~re sym.metric we deduce from (2.67)
.... that E = ET also.
So far the introduced strain tensors operate so]ely on )he material vectors a0 , X. ·..· Thus, C, .E and their inverse are also referred to as material strain tensors.
:·.·$pati.al strain tensors. In order to relate strain measures to quantities which are
··:_·.:.associated with the current configuration we continue with arguments entirely similar
._:..· ~<? those just used.
. . =·
....
...
The geometry in the .current COJ{/iguration is given by
|
|
y-x· |
|
|
|
||
y = y + (x - x) = x + IY - xi I |
I = x + dx |
|
(2.68) |
||||
|
. |
y-x |
|
|
|
|
|
dx = dfa |
and |
d£ = IYxi _, |
a= |
y-x |
. |
(2.69) |
|
ly-xl |
|||||||
|
|
|
|
|
|
||
::··>··The (spatial) .length of the spatial line -element dx |
== y - x is given by d€ = IY - xi, |
· wi_th d£/lxl ~ 1 (lxl # 0). The unit vector .a, lal |
= 1, acts at the current position x |
:_\..···.·. ~nd points in the direction of the spatial line element, which is the direction of A110 (see :...> figure 2.7). Since I-al = 1., we find using (2.69) 1 that
dx · dx = dea ·dE.a = di~ |
(2.70) |
·.The vector a may be viewed as a spatial element characterizing the direction .of a fiber
.··in the current configuration.
80 |
2 Kinematics |
time t
lal = 1
time l = 0
AH- 1 = F- 1a
IA;l·I=A-I
Figure 2.7 Deformation of a spatial line element with length d£ into n mate1ial line clement with length A- l d€.
In order to perform the relationship between ).no and a we recall eqs. (2.59)i,
(2.61):.i and definition (2.69h to give Fa0 = ...\a. Hence, from (2.60) we deduce that
An0 = A8 . |
(2.71) |
Note that the vector quantities a and dx are naturally associated with the current con- jigu111tion of the body.
Using Taylor's expansion, the associated position vector Y E S20 , which is described by the inverse motion x-1 (y' ·t)' may be express-ed by -means of (2.68), (2.69) and the inverse of the deformation _gradient (2.40). By the chain rule we find, by anal-
·Ogy with eq. (2.58), that |
|
Y = x-1 (y, t) = x-1 (x + dfa, t) |
|
= x-l(x,t)+dEF- 1(x,t)a+o(y-x) . |
(2.72) |
2.5 Strain Tensors |
81 |
In an analogous way to (2.-60), we now define the stretch vector .A0 |
in the direction |
ofthe unit vector a at x E s-2, i.e.
(2.73)
The length of a material line element (originally in the direction of a) is obtained from
(2.72h ·byneglecting terms of order d£2 . With definitions ( 1.15)1 and (2.73), we obtain, by means of Y - X = deF-1a that
(2.74)
w~ere the length of the inverse stretch vector .,x;1 is the inverse stretch ratio ,.\ -1. (or -~imply the (see Figure 2.7).
The square of A- 1 follows with definitions (1.15.h, (2.73) and property (1.81) as
|
|
|
|
(2.75) |
or |
-1 |
= |
p-tp-1 |
(2.76) |
l1ab |
A.a Ab · |
T.he strain tensor b- 1 is the inverse of the left Cauchy-Green tensor b, which is de- fined by
or |
(2.77) |
_(F ..is on the le.ft). In the literature the left Cauchy-Green tensor b is sometimes referred to· as the Finger d.eformation tensor. lt is an important .strain measure in terms of :spatial coordinates. The left Cauchy-Green tensor is symmetric .and positive definite at
each x E 0,
:. |
'l' |
'1' rr |
= b |
'r |
and |
u · bu > 0 for all u # o . |
(2~78) |
·... |
b = FF ·= (F |
F) |
|
||||
\VI"th definition (2.77) and ( 1.. l0 I), (2.5 I) we find consequently that |
|
||||||
|
|
|
|
detb = {det.F)2 = J2 > O . |
(2.79) |
||
:-.._·: \·..· As a last important strain measure we define the change in the squared lengths, i.e.
de2 - (,.\- 1di) 2 . |
With (2.75h, the use of the unit tensor I and eq. (2.70h we find a |
||
relation expressed sole.Jy through quantities in fl. Thus, |
|
||
~(dt2 - |
(,\- 1 d£")2] = |
~[dt2 - (dfa). F-Tr 1(cl€a)] = clx ·edx , |
(2.80) |
2 |
|
2 |
|
e =!{I - rTF- 1) |
or |
(2.81.) |
|
2
82 2 Kinematics
Relation (2..81) describes a strain measure in the direction of .a at place x E 0. We have
.introduced the commonly used symmetric strain tensor e, which is well-known as the EuJ.er·Almansi strain tensor.
Since the strain tensors b, e and their inverse operate on the spatial vectors a, x, we
call them spatial strain tensors.
·we now consider the case -C = FTF = I. From the relation (2.56h we know that the distance between any two neighboring points .X and Y located in the reference
configuration .is de.- Xf. On the -other hand the distance between the associated
neighboring points x and y located in the current configuration is given via (2.6 I )a, i.e. ,\de: = ly--xJ. = I we conclude from (2.62),1 that the line element is unstretched, i.e. ;\ = 1, and consequently de = fy - xi = fY - XI~ Hence, the distance between any two points is unchanged during such .a mot.ion. This means that there is no relative motion of points under x. Since C = I we conclude additionally from relation (2.67) that the strain tensor E vanishes identically, which means that the body does not change its size and shape {no changes in distances and angles).
This particular motion, which preserves the distance between any pair of points of a continuum body, .is called a rigid-body motion and is dealt with in more detail in Section 5.2. Hence, a rigid-body motion induces no strains .and consequently no stresses. A body which is only able to undergo a rigid-body motion .is said to be a rigid body. The idealization that a body is rigid is often cons-idered in engineering dynamics.
Push...forward, pull-hacl\ operation. As already seen_, ve-ctor and tensor-valued
quantities may be resolved along triads of basis vectors belonging to either the ref- erence or the current configuration. Additionally, there are two-paint tensors which are associated with both configurations, one ex.ample being the defonnation gradient (2.39). The transformations between material and spatial quantities are typically called and a pulJ...back operation (familiar .in differential geom...
etry) and are denoted by short-hand x* (•) and x;- 1( •), respectively. .In the literature the pull-back operation is often written as x* (•).
ln particular, a push-forward is an operation which transforms a vector or tensor- valued quantity based on the reference configuration to the current configuration.. Since the Euler-Almansi strain tensor e is defined with respect to spatial coordinates we can compute .it as a push-forward of the Green-Lagrange strain tensor E, which is given in terms of material coordinates. From eq.. (2.81) we conclude, using definition (2.67)~
that
e = ~(I - F-TF- 1) = F-T[~FT(I - F-TF-1)F]F-1
2 |
2 |
= F-T[~{FTF - |
I)]F- 1 = F-TEF-1 |
= x*(E) . |
(2.82) |
|
.2.5 |
Strain Tensors |
83 |
|
A pull-back is an inverse operation, which trunsforms a vector or tensor-valued |
||
|
quantity based on the current configuration to the reference configuration. Similarly to |
||
|
the above, the pull-back of e is |
|
|
|
E = ~(FTF- |
I) = FT[~F-T(FTF- 1)F-1]F |
|
|
2 |
2 |
|
|
= FT[~(I - |
F-TF- 1)]:F .. F'1'eF |
|
..:.=·. |
? |
··· ... |
|
|
= x; 1(e) . |
|
(2.83) |
|
As can._ be seen from eqs. (2.82) and (2.83) the transformations are based on multi- |
||
.· |
plications by one description of the deformation gradient, i.e. F, F- 1, FT, F-T. Which |
||
r~~m of the deformation gradient we have to take depends on the tensor to be trans- |
|||
(::_.:· |
formed. |
|
|
|
Following ·MARSDEN and HUGHES [ 1994] we indicate covariant tensors by (• )·~ |
||
|
and contravaritmt tensors by (• )d |
(for the notions 'covariant' and 'contravariant' the |
|
.. ·. |
reader is referred to Section 1.6). The push-forward and pull-back operations on co" |
||
.....:...:. |
i1aria11t second-order tensors (such as Eb, C~, eb, (b-1)~) are according to |
||
|
|
|
(2.84) |
:.'.-.-:_.·An example was given previously in eqs. (2..82) and {2.83), Le. e |
= x*(E~) and E := |
||
::·.·:..::··· :x;1 ( eb), wh.ich provide the relations.hips between the -material and spatial quantities, _:_:·_:·:.:::.... ~~e. e = F-TEF- 1 and E = FTeF, respectively.
·However, the push-forward and pull-back operations on con.travariant second-order
tensors (suc11 as (c-1) tt, bu and most of the common stress tensors) are according to
. .. ... . ~: :·.
(2.85)
·· :: ·. ~-~ t~e following chapters we use covariant strain tensors in combination with con..
"travariant stress tensors.
/·_:._._.....: · For completeness we write down the push-forward and pull-back operations on
.._covariant vectors~ i.e..
(2.86)
..
... and of contravariant vectors, i.e.
(2.87)
·:..:·.- _.. Finally we provide the so-called P.iol.a transformation of a spatial vector field
if~ u(x, t), with components u 0 , i..e.
(2.88)
2 Kinematics
where U - U(X, t) denotes a material vector field with the components UA. The transformation from u to U involves the pull-back of .u scaled by the volume ratio J,.
The inverse of eq. (2.88), i.e. u = .1- 1x*{U), involves the push-forward operation on
U here scaled by the inverse of the volume ratio, i.e. .1-1•
For a more complete source on the underlying concept the =mathematically oriented reader is referred to the book by MARSDEN and HUGHES (1994, .and references therein].
EXERCISES
1.For a given material point ex.press conditions on the right Cauchy-Green tensor C which ensure that
(a)no stretch occurs in a specified directfon ao of a fiber,
(b)no change in the angle between a pair of specified directions (a0 1, a0 2 ) takes place, and
(c)no change occurs in an infinitesimal surface element els placed in a plane
perpendicular .to a given direction (set ds = dS)..
2. Let a body undergo a homogeneous deformation defined by
,.. |
- |
|
r |
|
, X |
t ·' |
|||
.c 1 |
- |
a... |
|
|
where a, /J, 'Y are constants.
(a)Determine· the components of the material and spatial strain tensors C, E and b? e, respectively.
(b)Take the values /J = -cos() and 'Y = --sinO and show that if a = 1 the
strains are zero. Explain that for a = 1 the map corresponds to a rotation of magnitude fJ about the ){1-axis.
3. If the deformation of a body is defined by
, ............. |
.. |
' |
•1,3 - |
...A3 |
the body is said to be in a state of plane strain. Defonnations and strains occur only .in planes :i:a = const .and do not depend on
x3 are zero,. Determine the matrix representations of tensors F, E, e.
2.6 .Rotation, Stretch Tensors |
85 |
·.·...-·
4.Using relations (2.45h and (2.46:h show that the strain tensors E and =e may be expressed in terms of the displacement gradient tensor according to
··.... |
|
|
E = ~(Gracfru + GradU) + ~GradTUGradU |
(2.89) |
|||||||||
;: - · |
|
|
|
|
|
|
|
|
|
|
|
|
|
i//::...-: |
.:· |
.= |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
T |
u + gradu) - |
] |
|
T |
u gradu |
., |
(2.90) |
||
|
|
|
·e = -{grad |
|
.:..:.grad |
|
|||||||
|
|
|
2 |
|
|
|
2 |
|
|
|
|
|
|
....... |
·.·· |
or in index notation, in terms of material and spatial coordinates, as |
|
||||||||||
|
|
|
|
||||||||||
·.··.. :. ·. |
1 ( 8U11 |
DU.4) |
1 DUc 8Uc |
|
|||||||||
|
|
|
EA11 = 9 |
|
!1-:t·'" |
+a".\·"".. |
+?fl)·.. |
a~·r |
.11 |
|
|||
|
|
|
- |
|
u .."\.A |
.J.. \ JJ |
|
"-'u ...<'\.A |
.il. |
|
|||
.... |
1 ( Duh |
8ua ) |
1 Bue Due |
|
|
||||||||
: . . . |
|
|
|||||||||||
·.··. ·.. |
: . |
|
|
|
+ -,- |
- --- |
|
|
|||||
|
|
|
Cab= - |
|
- |
|
|
||||||
:··(.>·.-::· |
:·. |
2 |
|
8:1:" |
[};1:,, |
2 |
8:1:(J a~rb |
|
|
||||
:···. |
|
|
with implied summations on C and c, respectively. |
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
_:·.=·:_: |
|
|
|
|
|
|
|
|
|
|
|
|
|
,::/..:\.::::·:.· |
|
|
|
|
|
|
|
|
|
|
|
||
{<2.6 |
R.otation, Stretch Tensors |
|
|
|
|
|
|
|
|||||
......::.:-:.-·----------~------------~------~
..:::_/·:.
:.._:._.:{in the following we decompose a local motion, characterized by the nonsingular (in-
'./ vertible) tensor F(X, t), into a pure stretch and a pure rotation. As above, the argu-
::.....\)iients of the tensors are omitted, for convenience. However, they will be employed
;:.::·;:-.:·:::\Vhen additional information is .needed.
·:.· ......
;/ P~tardecomposition. At each point X E 0 0 and each time t, we have the following i: (~nique polar decomposition of the deformation gradient F:
...... :...... |
·.. |
F =RU= vR |
or |
(2.91) |
|
|
|
||
|
|
RTR=I |
V = VT |
(2~92) |
///.: rhis is a fundamental theorem in continuum .mechanics,. |
|
|||
:(::\\-... : |
|
In (2.91.) U and v define unique, positive definite, symmetric tensors, which we call |
||
/_:·/._the right (or material) stretch tensor and the left (or -spatial) stretch tensor, respec-
:·_:.·:.:.:_:_·_:tively. They =measure local stretching (or contraction) along their mutually orthogonal
·.·. :eigenvectors, that .is a change of local shape. The right stretc·h tensor U (w:ith compo- ·:-..'.. ·.~ents U;u1) is defined with respect to the reference configuration while the left stretch ):):.~te_~s~r v (with components 'V.ab) acts on the current configuration. Note that in this
.-//·\·.~ext the symbol U also stands for the material displacement field and v for the spatial
/·'.(._"ye"locityfield, respectively.
:;_}:_···)::}· ;·... |
The positive definite and symmetric tensors U and v are introduced~ so that |
|
||||||
....... |
u |
2 |
|
|
2 |
|
|
|
..:_.:_:.:_ :: ...__ :-:. |
and |
v |
= vv = |
.b ' |
(2.93) |
|||
... |
= uu = c |
|
||||||
.86 2 Kinema.tics
which is based on the square-root theorem; see, for example, GURTIN [198 la, p. 13]. Consider the relation detC = detb = J 2 we deduce with (2.93) and rule (1..101) the important property
detU = detv = J > 0 . |
(2.94) |
The unique R is a proper orthogonal tensor, with detR = 1_, called the rotation tensor. It measures the local rotation that is a change of local orientation. A rigidbody -rotation about a fixed origin is then characterized {f and only {f U = v = I, so that F = R. Hence, each material line element dX is rotated into a unique spatial 1.ine
element dx (and vice versa) according to (2.38), i.e.
dx = R{X, -t)dX or cb;a = RaA cLY_.\ . (2.95)
On the other hand if R = I, the .deformation is called pure stretch, for which ·.F = U =
·v. Using (2.95), the pointwise polar decomposition (2.91.) may be written as
dx = .R(UdX) |
or |
d:1:,1 = RaA (U;1[Jd.X.r:d , |
(2.96) |
dx = v (RdX) |
or |
&r:a. = 'Vab (Rb.•\ d..Y..:l) . |
(2.97) |
.Relation (2.96) describes a pure stretch of dX by the material tensor U, (UAu), followed by a pure rotation performed by the rotation tensor R, (RaA)· The rotation tensor transforms the material vector UdX into the spatial vector dx. However, relation (2.97) describes a pure rotation of dX performed by the same rotation tensor R, (RaA), which transforms .the material vector dX into the spatial vector RdX. Th.is rotation is followed by .a pure stretch of RdX with the spatial tensor v, (vab), which gives the spatial vector dx. In both cases the rotation tensor R maps between the reference and the current configuration and, therefore, R is a two-point tensor like F.
Relation F = RU is also known as the unique right polar decomposition, while
F = vR is referred to as the unique left polar decomposition.
Show that R is proper orthogonal, i.e. RTR = I and dPtR = 1.
Show further that the multiplicative decomposition F = RU is unique.
Solution. Using R = Fu- 1 and eqs. (2.93)1. and (2.9?h we find that
1 |
1 |
2 |
u- |
1 |
R'l'R = (Fu-lfr(Fu- ) = u-TFTFu- |
= u~Tu |
|
||
= (uu-1fruu- 1 = I . |
|
|
|
(2.98) |
In addition, computing d<;~tR = det (Fu-1 ) we deduce with rule ( 1.10 O and relations (2.51) and (2.94) that
1 |
= 1 . |
(2.99) |
detR = detF·u- |
||
c1et- |
|
|
2.6 Rotation, .Stretch Tenso..S |
87 |
.. .. . In order to show uniqueness of the polar decomposition we suppose that there exist
.(i. ·pos.itive definite, second-order tensors R and UT so that
F =RU= RU |
(2..100) |
|
RTR =I |
- -'r |
(2.101.) |
and U=U |
||
.·...·Then, with ·C = .FTF and eqs. (2.100)2 and (2.101) it follows that
(2.102)
·::...· ·\vhich means using (2.93) 1 that U = U, since C has a unique square root, and hence
.:.: R = R. The proof is similar for F = vR. •
.. We discuss further a physical .interpretation of the right stretch tensor U. Look in ::... · itie direction of n0 at point X E n0 , with laoI = 1. According to (2.60) we find using
.( ¥. = RU that the stretch vector A00 may be expressed as
Au0 = RUao . |
(2.103) |
;>·Then, the square of the stretch ratio ,.\ = (A00 • ~ao)112 may be computed with (2.103 ),
.-._:·.identities ( 1.8 .I), ( 1. .84) and relation RTR = I. Thus, we have
A2 = RU.ao · RUa0 = n0 • U'I'(.RT R)Ua0
= Uao · Uao = IUaof . |
(2.104) |
. Note that the length of A.110 acting .at X E Ho along the direction of the unit vector ao depends only on a state of pure stretch, i.e. U; therefore we call U the stretch tensor.
Another form of A2 may be obtained from (2..104:h using (2.92h, i.e. |
|
.A2 = a0 · UT Ua0 = a0 • U2a0 • |
(2.105) |
Consequently, the right Cauchy-Green tensor C (recall eq. (2.93)i) and the Gree.n- Lag.range strain tensor "E = (C - I) /2 do not include information about the rotation,
that is experienced by a particle during motion. The deformation gradient F ·contains more information th.an strains do in general (rotation and strain-like information).
Finaily we find the relation between v and U, and between b and C. By me.ans of
the polar decomposition (2.91) we have
v = FR'l' = RURT |
and |
(2.106) |
88 2 Kinematics
Consequently, with v2 = RU2 RT and eq. (2.93) we find an important relation between the left and right Cauchy-Green tensors, namely
(2.107)
A general formula for strain measures was introduced by SETH [1964] and HILL
The generalized strain measures in the Lagrangian and Eulerian descriptions
are defined as
!(un - I) , |
.!.(vn - |
I) , |
if |
n-#0 |
(2.1.08) |
n. |
n |
. |
|
n = 0 : } |
|
lnU , |
Inv , |
|
if |
|
where n is a rea] number (not necessarily an integer).
For the special cases n = 0 and n = 1 we obtain strain tensors associated with the names .Henck)' and Blot, respectively, while for n = 2 and n = -2 we obtain the
Green-Lagrange strain tensor (U2 - 1)/2 (or also (v2 - 1)/2) and the Euler-Almansi
.strain tensor (I - u- (or also (I - v-2 )/2), which we have introduced in (2.67) and (2.81 ), respectively. For more details see, for example, MAN and Guo [1993].
The Hencky strain tensor in the material and spatial form, Le. lnU and Inv, is of particular interest in .nonlinear constitutive theories. Hecause of the logarithmic func- tions, these tensors are decomposed additively into so-called volumetric and isochoric parts (for a more detailed exposition of the decomposition .procedure of strain measures see, for example, p.
EXAMPLE 2.7 Assume a certain two-dimens.ional ·motion which is given in the form of the deformation gradient F(X, t) = Bx(X, 't)/DX. All tensor quantities are defined with respect to the orthonormal basis {ea}, c\! = l_, 2.
Show that fo.r the case of two dimensions the polar decomposition F = RU may be
given by the closed-form expression |
|
U = [11 (C) + 2J(C)tt./2..[C + J(C)I] |
(2.109) |
for the right stretch tensor, and by the rotation tensor R which follows by means of the inverse of expression (2.109), i.e. R = F0- 1• In eq. (2.109), / 1(C) = C11 + C22 and
J(C) = (C11 C22 )"112 characterize the first scalar invariant of the right Cauchy-Green tensor C and the volume ratio, respectively.
Solution. |
In two dimensions, the characteristic polynomial of U gives a quadratic |
|
equation in ,\, i.e.. |
|
|
|
.a= 1, 2 , |
(2.1.10_) |
