Applications of Tensor Functions in Solid Mechanics
Applications of Tensor Functions in Solid Mechanics
EDITED BY
J.P. BOEHLER
UNIVERSITY OF GRENOBLE
Jean-Paul Boehler
CONTENTS
1. Introduetion
2. Different domains of meehaniea1 anisotropy =
3. Essential features of the anisotropie meehanieal
behavior of rolled sheet-steel 6
1. Introduetion 13
2. Prineiple of Isotropy of Spaee 14
3. Isotropie materials 16
4. Anisotropie materials 18
5. Orthotropie materials 21
6. Representation of the funetion F 23
7. eonclusions 29
1. Introduction 31
2. Representations for isotropic scalar and tensor
funetions 35
3. Representations for non-polynomial anisotropic
scalar and tensor functions
3.1. Method based on the introduction of
structural tensors 40
3.2. Generalization of the Rivlin-Ericksen
method 41
3.3. General anisotropy 42
3.4. Orthotropy 44
3.5. Transverse isotropy 45
3.6. eomparison with representations for
polynomial anisotropic tensor functions 47
4. Representations for non-polynomial isotropic and
orthotropic tensor functions in a two-dimensional
space 49
1. Introduction 67
2. General theory 68
3. Plastic deformations in uniaxial and triaxial
tests
4. Failure criteria for glass/epoxy composites
under confining pressure
4.1. Introduction 80
4.2. General form of the yield criterion for
triaxial tests 80
4.3. Failure modes and directional strengths 81
4.4. Proposed failure criteria 84
4.5. Comparison with experimental results 87
5. Simplified theory
5.1. Plastic behavior 89
5.2. Yield eriteria 93
1. Introduction 203
2. Integrity basis for a second-order tensor 204
3. Simplified characteristic polynomial 207
4. The Hamilton-eayley theorem 210
5. eonstruction of simultaneaus invariants 212
6. eonstruction of invariants by the polarization
process 214
7. Extended characteristic polynomial 215
8. The Lagrange multiplier method 221
1. Introduction 227
2. The damage state in a continuum 229
3. Stresses in a damaged continuum 236
4. Constitutive equations involving damage and
initial anisotropy 240
1. Introduction 251
2. Tensor function of one variable 252
3. Tensor function of two variables 256
4. Interpolation at coincident points 260
5. Polynomials of second-order and fourth-order
tensors 262
6. Simple examples 263
7. Tensorial generalization of Norton's creep law 265
8. Separation of tensor variables 274
1. Introduction 279
2. Isotropy 280
3. Oriented solids 283
4. Modification of the classical flow rule 288
5. Anisotropy expressed through a fourth-rank
tensor 295
LIST OF AUTHORS
Templergraben 55
D-5100 AACHEN
Federal Republic of Germany
IN SOLID MECHANICS
Chapter 1
PHYSICAL MOTIVATION
J.P. Boehler
University of Grenoble, France
1. INTRODUCTION
E (GPa)
225 R0 (MPa)
142
~ ,001
v
139 \
? ~
215 I 136
,
~~
I / ~ ...
133 L
' "
I
I
I
210 ,;ex peri men tal curv e
! 130
../ i
205 I 127 -
124
200
0 30 60
0 30 60 90
!
[!oo
f R 0 (MPa)
R E (MPa)
10
a) Diatomite
0,2
5
10 5 ,---------t---~----~---t---~-----,
psi
12
fO
r
9 0.6
0.4
0.2
6n t
0 '----~~~~3~0~~4~5--~60~--~~~~goo 0 '----1-;j:5:------;;3~0-----:;4~5---6;:;0;-------:;7::;:5--~90°
ß-- ß-
400
'LIR45
390 17 LI
~
~ -c:::J·
380 R90 1
370
LL V L1
'
///
360
RoD
350 •I
Fig.7 -Anisotropie stress-strain
relations. 340 ~V
330
_.!:_ /ß"
Re ;:l V
'~·
320
5 10 15
anisotropy
Fig.lO - Evolution.of oriented
I
yield stresses with the
:-isotropy
rolling reduetion.
eo
0 15 30 45 60 75 90 R9 (MPa>
420
IL_R4 R45
Fig.8 - Anisotropie yield 410 L.
eriterion in simple
tension. 400 iL VRsjo
/V
Re(MPa) Series
390 l Series C
r=15% V :L
Rgo
V
380 370
370 I
L ll
360
1/
L_ / v ~
l
I
I~ I/V
~ I
Series A I
I r=s·;.
I
360 350 '
_L ~
-
I
........
I
~a
350 / II
II '-...... 340
..... ~ .... B 10%
I/I_
II
I
34 0 I 330
I II
,...-
I
r-- r---. E%
33 0
=~ A
90
320
0 2 4 6 8 10 12 14 16
32 0
0 15 30 45 60 75 90 Fig.ll - Evolution of oriented
yield stresses of rolled
Fig.9 - Evolution of the yield sheet-steel with further
eriterion with the rolling applied irreversible
reduetion. prestrains.
Physical Motivation 9
(prestrainings) ?
Fig.ll represents the evolution of the yield stress in the rolling and
transverse directions and at 45°, with respect to the prestrainings €
58
E =0,020
~-----r--_____",,-------".,-. isotropy
90
0 15 30 45 60 75 90 0 15 30 45 60 75 90
62
necking
I ine 2 necking
line 1
R
a
0 15 30 45 60 75 90
Fig .12 - Evolution of the orientation of the necking lines (after [8]).
REFERENCES
J.P. Boehler
University of Grenoble, France
1. INTRODUCTION
T F(D, 0 ( 1)
Initial configuration
F
...
of the material
Agency D Response T
the material
F depends on ~
the behavior
Constitutive relation
D T
} =<> ara'
Finally we obtain:
3. ISOTROPIC MATERIALS
E1 e1
Qv2~ Ov2\
\
Qv1
__.,.
- --av 1
-- v1
~
v1
===t>T
response T are the same. To see that, we apply the reflection S with res-
pect to the direction ~1 to both the material and the agency. We obtain
obviously:
D => D D ( 3)
T => T (4)
T = T (5)
STSt = T (6)
Isotropie material
Finally, we obtain:
VQ E 0 F(D, ~) (7)
It is easy to see that the condition (7) means that the structural tensor
~ is an isotropic tensor:
~=AI
4. ANISOTROPie MATERIAL
39
-lc
~ ?
E 0 : ! = ( 10)
If only the transformations (I, -I) satisfy (10), then the material is
called "general anisotropic material". If other orthogonal transformations
9 satisfy (10), then they constitute a subgroup S of the orthogonal group
0. The subgroup S is called the material symmetry group.
For examp1e, the material symmetry group for a two-dimensional ortho-
tropic material, with privileged directions (v , v ), is given by:
-1 -2
Introduction to the Invariant Formulation 19
s = ~ ± I, S , S
- -1 -2~
l (11)
v1
Anisotropie material
QE S <=) ~ (12)
It follows that:
Q E S =) F(D, ~) F(D, 0 (13)
Conversely, suppose that the relation (15) holds for a given transforma-
*
tion Q . Employing the Principle of Isotropy of Space (14), we obtain:
Finally, considering the relations (13), (16) and (17), we obtain the
following converse implication:
5. ORTHOTROPIC MATERIALS
M (19)
M
cosa. cosa.
n => d
-a. sina. -a. 0
M -1 +1 v,
Thus, the tensor M transforms the unit circle into the segment [ -1, +1] lying
on the axis. The principal directions of M are the privileged directions
~1
M
cos 2 8 s~n8cos81 ( 2 0)
sin8cos8 s1n 28
(x,y)
It is easy to see that the tensor M verifies the following converse impli-
cation:
(21)
Thus, the invariance group of M is the material symmetry group (11), i.e.
M is the structural tensor for a two-dimensional orthotropic material.
We now apply the fundamental result of Section 4. Consider the cons-
titutive equation:
T F(D, M) (22)
6.1 Introduction
T a I
0-
+ a 1-M + a 2-D (25)
a. a.(trD, trD 2 , trMD)
1 1 - - --
I, M, D (26)
- - -
which constitute the generating set of the representation. More precisely
the generators (26) are known s - invariant symmetric second order ten-
sors, where s is the material symmetry group (11) for two-dimensional
24 J.P. Boehler
E E
11 12
D (28)
trD E
1
+ E
2
E
11
+ E
22
tr(QDQt) , VQ E o
t
(29)
E2
1
+ E2
2
= E 11 + E
22
+ 2E 2
12
= tr(QDQ
---
)2 , VQ E 0
Thus, trD and trD 2 admit the same values for arbitrary orthogonal trans-
formations of D: they are called isotropic invariants of D. As they are
invariant under the full orthogonal group O, they are invariant under the
group S a fortiori ; thus, they constitute also orthotropic invariants
of D.
Introduction to the Invariant Formularion 25
trMD E
l l
e (v ,
-l
E
-l
) (30)
VQ E o: F(D, M) = a I
0-
+ a l-M + a 2-D
Q(a I
- 0-
+ a lM
-
+ a 2-D)Q- t (34)
= F(D, M)
T a I + a M + a. D ( 35)
0- l- 2-
y2 E1
v2 E1
<F
v1
+
~ v1 +
T -- o0 1 + o1 M + o2 D
Introduction to the Invariant F ormulation 27
a
11
a
1 2
1 0 sin8cos8 E:
1
0
Cl
0
+ Cl
l
+ Cl
2
(36)
a1 2 a2 2 0 1 sin8cos8 0 E:
2
where
8 (v , E )
-1 -1
a a 1 sin8cos8 (37)
12
M
-2
= I - M = v-2 e v
-2
(38)
28 J.P. Boehler
v2
+1
M
-2
cosa 0
n => f M2:
-Cl
sina
-Cl
sina v,
-1
Thus, the tensor ~2 transforms the unit circle into the segment [-1, +1]
1ying on the ~2 axis. Consider now the structural tensor:
M·
2'
-1 +1 V
1
v, -1 +1 V
1
-1
~
-
= a(S)M- + b(S)M
-2
b(S)I + [a(S) - b(S)]M (41)
7. CONCLUSIONS
REFERENCES
J.P. Boehler
University of Grenoble, France
1. INTRODUCTION
A
_l
• A
_2
• . .. A
_a • w
_l
•w
_2
• • .• Wb.
_
V •
_l
V •
_2
• •• V
_c (1 )
VQ e: s (3)
l
T = a..G.
1-1
where the a.. are arbitrary polynomial scalar functions of the invariants
1
of the integrity basis. The values of the functions a.. are thus polynomial
1
scalar invariants of the arguments (1) under the group s. The representa-
tion (6) is irreducible if the integrity basis is irreducible and if none
of the generators G. can be expressedas a linear combination of the others,
-1
with coefficients a. ..
1
Explicit representations for isotropic and anisotropic polynomial
scalar-valued and tensor-valued functions are presented in Chapters 8 and
9 by Professor SPENCER.
T a..G.
1-1
but, here, the scalar valued-functions a.. are general functions of the in-
1
variants of the functional basis of the arguments (1). The representation
(7) is irreducible if the functional basis is irreducible and if none of
the generators G. can be expressed as a linear combination of the others,
-1
formed with general functions a. .. As is the case for the bases for scalar
1
functions, irreducible generating sets for general tensor functions con-
tain, in general, fewer elements than do the corresponding generating sets
for polynomial tensor functions, because here reductions are obtained with
the help of linear combinations, the coefficients of which are general
scalar-valued functions and not only polynomials.
Representations for Tensor Functions 35
sublists.
The procedure used by SMITH in order to construct the functional
bases for general isotropic scalar functions f constitutes an extension
of the method developed by RIVLIN and ERICKSEN in [10], applied to a par-
ticular case of the arguments (1). The method consists in establishing a
set of basic invariants such that once the values of these invariants are
known, the components of the arguments (1) can be determined uniquely in
a reference frame fully specified beforehand (orientation and sense uf the
Coordinates are specifjed). Such a set of invariants constitutes a func-
tional basis since any invariant is uniquely expressed in terms of these
components. The required isotropy is assured by the preliminary choice of
the reference system, which must be defined from the arguments (1) them-
selves. In order to establish generating sets for the general isotropic
tensor functions F, SMITH employs WANG 1 s procedure.
The representations obtained by WANG and by SMITH are not identical.
Moreover, the functional bases derived by SMITH contain redundant elements.
This has been analysed in [11], where after suitable amendments of the
previously established bases, an irreducible functional basis is obtained.
Reasonings similar to those employed in [11] allow conclusions to be drawn
for the generating sets as well. The complete and irreducible representa-
tions so obtained are given in the Tables I and II below.
For example, the irreducible representation for an isotropic tensor
function of two symmetric second order tensors :
T F(A, B) (8)
(9)
Variables Invariants
w ' w w
-2 ' -3
trW ww
-1-2-3
- 1
A v•A
- 1 ' ~2'
V
A V
- -1- -2-
A, V ' ·Av , •A 2 v
- - 1
V
- <
V
-1 -- 2
V
-1 - -2
w v•W v, W2 v·W v, w v•W 2 v
-l ' w
w -2 ' -
V
- l - -2- -1- -2- -1- -2-
w, V- 1 ' V-2 V ·Wv V •W 2 v
-
- - 1 .... -2' - l - 2
A, W, Av·Wv, A2v•Wv, AWv•W 2 v
- - - -- --- - -
V
-- --
A ,A , V A V ·A - -A1-V 2 ·A- 2-V 1
V
-1 -2 - 1'
V -2 -1-1 -2-2
A, w, Av •Wv - Av •Wv
- -
V
-l '
V -2 -- 1 --2 --1
--2
w w V w V ·W V - w V ·W V
-1 -2' -1 '
V -2 -1-1 -2-2 -1-2 -2-1
Variables Generators
I
-
A A, A2
- -
w wz
- -
V V
-
V
Qll
- -
A A A A + A A , A2 A + A A2 , A A2+ A2 A
-l ' -2 -l-2 -2-l -1-2 -2-1 -l-2 -2-l
- w
A,
- AW-WA, WAW, A2 W-WA 2 , WAWLW 2AW
---- --- - - -- --- - --
A2 v + A2 v
- -- + -- - -
A, V Av Av v, V
- - -- -- -
V Qll Qll Qll Qll V
w w
-l ' -2
w w+ w
-l-2 w' w
-2-l
W2 + W
-1-2
2W '
-2-1
W2 W+
-l-2
w wz
-2-l
w, Wv Wv, Wv + Wv v, Wv W2v + W2 v Wv
V
-- -- V- -- -- --
Qll Qll Qll Qll Qll
- - - -- - - --
V
-l '
V
-2
V
-1
Qll V
-2
+ V-2 1111 V
-1
A, V (v ® Av + Av Qll V ) - (v ® Av + Av Qll V )
- V
-1 ' -2 -1 --2 --2 -1 -2 --1 --1 -2
\:!. '::\' V
-2
(v
- 1
® Wv
--2
+ Wv
-- 2
Qll V
- 1
) - (v- 2 ® Wv
-- l
+ Wv
-- 1
Qll V
- 2
)
A, W, ... Wb'
-a -1 -
v,
-l
... v)
-c
is given by : F = a.G.,
-1-1
where the a. 's are arbitrary scalar-valued functions of the invariants
1
of the functional bases given in Table I and where the G. 's are the
-1
generators given in Table II, considering all (unordered) combinations
of zero, one, two and three variables.
I, A, A2 , B, B2 , AB+BA,
- ... - - - ( 11)
A , A ••. A (12)
-l -2 -a
40 J.P. Boehler
The integrity basis for the list (12) contains the traces of products of
up to six tensors [7], whereas the functional basis contains the traces
only up to the order three (Table I). Similarly, the generating set in
the polynomial case contains products of up to five tensors [31, whereas
in the general case, it contains products up to only two tensors (Table
II). The complete analysis for the reductions introduced in the represen-
tations when considering general functions instead of polynomial functions
is presented in [12].
T ~ •..• ~ )
• w_2 . . . . w_b' v_l. V-2 • • • • v-c • ~-1 • -2 (13)
F(A
- -1
• A -a • w
-2 •... A -1 -p
VQ e: 0 F(A.,
- -1
W.,
-J
v-k' ~-")
"' FÜ.,
- -1
W.,
-J
V-k' ~- "'n) (14)
where F, ~i' ~j' ~k' ~2.are the transformed of !• ~i' ~j' ~k' h by an orthogo-
nal transformation 9·
The condition (14) means that F is an isotropic func-
Representations for Tensor Functions 41
VQ E s F(A.,
- -1
W.,
- J
'!_k'
-
~0 )
·-Tv
F(A.,
- -1
W., '!_k' ~ 0 )
- J - -Tv
(16)
sentations for tensor valued functions are composed of the functional bases
and the generating sets of tensors. The functional bases constitute the
representations for the corresponding scalar-valued functions.
tion 2.) for isotropic functions. The list of invariants obtained is too
long to be reproduced here.
( 17)
M
-11
V
-1 • V
-1
M
-22
V
-2 • V
-2
M
-33
V
-3 • V
-3
M
-23
+ M
-32
V
-2
• V
-3
+ V
-3
• V
-2
M
-31
+ M
-13
V
-3 • V
-1
+ V
-1 • V
-3
(18)
-2 •
M
-12
+M
-21
V
-l
• V
-2
+ V V
-1
1 0 0 0 0 0 0 0 0
M 0 0 0 0 1 0 M 0 0 0
11 3 3
0 0 0 0 0 0 0 0 1
(19)
0 0 0 0 0 1 0 1 0
M
23
+ M32 0 0 1 M31 + M
13
0 0 0 M12 + M21 1 0 0
0 1 0 1 0 0 0 0 0
It is easy to show that the group (17) is the invariance group of the ten-
sors (18), which constitute a basis in the space of symmetric secend order
tensors. The tensors (18) are tt~ structural tensors for general anisotro-
PY·
Thus, the representation for the isotropic function
The analysis of the explicit expressions show that the invariants (21)
are just the independent components in the (v , v , v ) frame of the agen-
(") -1 -2 -3
cies A. : trMk 0 A.
-1 - JC-1
= Ak!"' . Any scalar invariant of the A. can be expressed
-1
uniquely with respect to the invariants (21). Thus, the invariants (21)
constitute a complete functional basis. As they are functionally indepen-
dent, this functional basis is irreducible a fortiori.
Among the generators for the isotropic function (20) appear the six
tensors (18), which constitute a complete set of tensor generators. As
the tensors (18) are linearly independent, this generatingset is irreduci-
ble.
Finally, an irreducible representation for the general anisotropic
tensor function F of the agencies A. is given by:
-1
T ~ M + u M + u M + u (M + M ) +
1-11 2-22 3-33 4 -23 -32
(22)
~. ~.(trM A., trM A., trM A., trM A., trM A., trM A.)
1 1 -11-1 -22-1 -33-1 -23-1 -31-1 -12-1
i=1,2, ... a.
This result could have been established directly. Indeed, the tensor
T can always be expressed as a linear combination of the basic tensors
(18), with the components Tmn of T as coefficients. As the function F does
not verify any symmetry, the components T are arbitrary functions of the
(') mn
6a independent components Ak~ of the agencies ~i'
44 J.P. Boehler
3.4 Orthotropy
s ?~±I- , S , S , S
-1 -2 -3~
l (23)
M
-1
V
-1
~ V
-1
M
-2
= V ~ V
-2 -2
M
-3
V
-3
~ V
-3
(24)
The expressions of the tensors (24) are given by (19). It is easy to show
that the group ( 23) is the invariance group of the tensors ( 24) :
QE S <==> QM.Qt
--1-
= M.
-1
i 1,2,3 . (25)
Thus, the tensors (24) are the structural tensors for orthotropy.
The representation for the isotropic function
MI?
-1
= M.
-1
trM.
-1
1 i = 1,2,3
M.M. = 0
-1-J
i +j i,j = 1,2,3
(27)
(M-1-A+AM ) + (M-2-A+AM
--2
) + (M-3-A+AM
--3
) 2A
--1
trM A + trM-2-A + trM-3-A trA
-1-
Representations for Tensor Functions 45
+ et 6 (M-3-A+AM
--3
) + a 7 (M-1-B+BM
--1
) + et 8 (M-2-B+BM
--2
) +
(28)
trM AB, trM AB, trM AB, trA 2 B, trAB 2 , trA 3 , trB 3 )
-1-- -2-- -3--
s (29)
M (30)
M
-1 ' M
M M
-2 ' -3 M
-1 ' M
-2 ' -3
A A A A +A A M A A +A A M , M A A +A A M
-1 ' -2 -1-2 -2-1 -1-1-2 -2-1-1 -1-2-1 -1-2-1
M A A +A A M , M A A +A A M
-2-l-2 -2-1-2 -2-2-1 -1-2-2
M A A +A A M , M A A +A A M
-3-1-2 -2-1-3 -3-2-1 -l-2-3
T ~ I + ~ M+ ~ A+ ~ B+ ~ (MA+AM) + ~ (MB+BM) +
l- 2- 3- 4- 5 -- -- 6 -- --
A ,A trA A , trMA A
-1 -2 -1-2 --l-2
A ,A ,A
-1 -2 -3
Generators for
Variables general Generators for
functions polynomial functions
I, M. I, M.
A A, MA+AM, A, MA+AM,
-
A2, MA 2+A 2M.
- -- --
A2 , MA 2+A 2M.
A , A
-1 -2
AA + AA , MAA + A A M, MAA
_2_1_ __ 2_1 + A A M,
_1_2 _
-1-2 -2-1 --l-2
A2A + A A2 ,
_2_1
MA 2A + A
__ 1_2 A2M,
_2_1 _
-1-2
A A2 + A2A. MA A2 + A
__ 1_2
2A M.
_2_1 _
-l-2 -2-1
A , A , A MAA A + A
__ 1_2_3 A A M, MA AA + A
__ 1_3_2 A A M,
-1 -2 -3 _3_2_1_ -2-3_1 _
MA AA + A
__ 2_1_3 A A M,
_3_1_2 _ MA AA + A
__ 2_3_1 A A M,
_1_3_2 _
MA AA + A
__ 3_1_2 A A M,
_2_1_3 _ AMAA
_1 __ 2_3 + A A MA
_3_2 __ 1
T = F(A , A , •. . A ) (34)
- - -1 -2 -a
T = ~0 I
-
+ ~.G.
1-1
(35)
where I is the two-dimensional identity tensor, the G. 's are the genera-
-1
tors:
. ...l
A., A~,
~k~~ + ~~~k' A-k-~
2A + A Ak'
-J -J -~-
(36)
~~ Ak + A A2
-k-~
k < ~ j, k, ~ = 1, 2,
50 J.P. Boehler
0 (38)
From the theorem (39), we obtain the following relations for the gene-
rators (36):
1
A~ = -(trA~- tr 2 A.)I + (trA.)A.,
-J 2 -J -J- -J-J
(tr~k)~~ + (tr~~)~k +
(40)
3
trA~
-J 2 trA.
-J
trA:
-J
(41)
Taking into account the relations (40) and (41), we obtain from (36)
and (37) the following irreducible representation for the isotropic func-
tion (34):
~.
~
= ~.(trA.,
~ -J
trA:, trAkA 0
-J - ~
), (42)
S = ~±I,
( - -2 -3)
S , S l (43)
M (44)
T F(A , A , .. . A , M) (45)
- -l -2 -a -
trM 1 (46)
T + Ql a+z-a
A
REFERENCES
J.P. Boehler
University of Grenoble
1. INTRODUCTION
T a I + a M + a D + a (MD+DM) + a D2 + a (MD 2 +D 2 M)
0- l- 2- 3 -- -- 4- 5 -- - -
(1)
a.
1
= a.(trD,
1 -
trD 2 , trD 3 , trMD, trMD 2 )
M V ® V (2)
-3 -3
T a M + a M + a M + a (M D+DM ) + a (M D+DM ) +
1-1 2-2 3-3 4 -1- --1 5 -2- --2
+ a (M D+DM ) + a D2 • (3)
6 -3- --3 7-
M V ® V M M V ® V (4)
-1 -1 -1 -2 -3 -3 -3
2. TRANSVERSE ISOTROPY
T a I + a M + a n + a (Mn+nM)
0- 1- 2- 3 -- --
a a + b trn + c trMn
0 0 0 0
(6)
a1 a 1 + b 1 trn + c 1 trMn
a
2
a2 a
3
a3
T
-u
a I
o-
+ a M
1-
(7)
D - 0 => T
-U
a I + a M
0- 1-
0 (8)
trT
-u = 3a +a
0 1
0
=> a a1 0 (9)
trMT
--u
= a 0 +a 1 0
0
Thus, when the solid admits a natural state, the transversely isotropic
linear elastic behavior depends on six material constants.
Suppose now that the elastic solid is a hyperelastic material, i.e.
its mechanical response T is the derivative of a scalar-valued strain-
energy function W:
aw
T (10)
an
The representation of the strain-energy function is obtained from Table IV
of Chapter 3:
aw + ~ M aw
T atrD ! atrMD - + 2 atrD 2 ~ +
(12)
Identifying (12) with the linear form (6) of the representation (1), we
obtain:
aw b trD + c trMD aw aw
atrD = o o 0
atrD 2 atrD 3
(13)
aw
a
3
Anisotropie Linear Elasticity 59
a aw => c b . (14)
atrD <atrMD) 0 1
T A D (17)
a f d 0 0 0
f a d 0 0 0
d d c 0 0 0
A (18)
0 0 0 g 0 0
0 0 0 0 g 0
0 0 0 0 0
1
-(a-f)
2
b f b d-f c c+a-2d-4g
0 1 1
(19)
a a-f a3 2g-a+f
2
- = [f trD+(d-f)trMD]I
T
... -- - + [ (d-f)trD+(c+a-2d-4g)trMD]M +
- .........
(20)
+ (a-f)D + (2g-a+f)(~~+~~)
3. ORTHOTROPY
T
-u
a M
1-1
+ a 2-2
M +a M
3-3
(22)
D - 0 => T
-u
M + a 2-2
a 1-1 M + a 3-3
M 0 (23)
Thus, when the solid admits a natural state, the orthotropic linear elas-
tic behavior depends on 12 material constants.
In the case of an orthotropic hyperelastic material, the represen-
tation of the scalar-valued strain-energy function is obtained from
Table III of Chapter 3:
T aw M + aw M + ~- M aw ( )
3trM D -1 3trM D -2 tltrM D -3 + tltrM D2 ~1~+~~1 +
-1- -2- 3- -1-
(26)
Identifying (26) with the linear form (21) of the representation (3), we
obtain:
62 J.P. Boehler
aw aw
b 1 trM D + c 1 trM D + d 1 trM D a4
atrM D -1- -2- -3- atrM D2
-1- -1-
aw aw
= b 2 trM D + c trM D + d 2 trM D a5 (27)
atrM D -1- 2 -2- -3- atrM D2
-2- -2-
aw aw aw
b 3 trM D + c 3 trM D + d 3 trM D a6 atrD 3 = 0.
atrM D -1- -2- -3- atrM D2
-3- -3-
aw aw aw aw
atrM D <atrM-2-D) atrM D <atrM D)
=> c1 b2
- 1- -2- - 1-
aw aw aw aw
atrM D <atrM D) atrM D <atrM D)
=> b3 d1 (28)
-1- - 3- - 3- - 1-
aw aw aw aw
atrM D ( atrM D) atrM D <atrM D) => d
2
c
3
- 2- - 3- - 3- - 2-
T = (b trM D + b trM D + d trM D)M + (b trM D+c trM D+c trM D)M +
1 -1- 2 -2- 1 -3- -1 2 -1- 2 -2- 3 -3- -2
+ a (M D+DM ) + a (M D+DM )
5 -2- --2 6 -3- --3
Comments on the invariant forms (29) and (30) are similar to those given
in Section 2 .1.
a f e 0 0 0
f b d 0 0 0
e d c 0 0 0
T AD A (31)
- 0 0 0 g 0 0
0 0 0 0 h 0
0 0 0 0 0 j
b
1
a+2g-2h-2j c
2
b+2h-2j-2g
a
6
g+h-j c
3
d d1 = e b
2
f
4. ISOTROPY
T a I
0-
+ a D + a D2
1- 2-
(34)
a. a.(trD, trD 2 , trD 3 )
1 1 - - -
T a I
0-
+ a 1-D
(35)
a
0
a
0
+ b0 trD
a
1
a1.
T a I (36)
-U 0-
D - 0 => T
-U
a I
0-
0 => a
0
0 (37)
Thus, if the solid admits a natural state, the isotropic linear elasti-
city depends on two material constants:
T (b trD)I
0
+ a 1-D (38)
w l2 b tr 2 D + 21 a 1 trD 2
0 -
(40)
REFERENCES
J.P. Boehler
University of Grenoble, France
1. INTRODUCTION
The essential features of the anisotropic flow law and their conse-
quences in experimental investigations of the mechanical behavior of aniso-
tropic solids are analyzed in Section 3. Specific forms of failure crite-
ria are proposed for glass-fiber reinforced composites and compared with
experimental data (Section 4).
A simplified theory employing a fourth order tensor of plastic ani-
sotropy and a tensorially linear description of the material response is
discussed in Section 5, in order to arrive at a warkable form of the cons-
titutive relation and at simple forms of the failure criteria for trans-
versely isotropic solids.
2. GENERAL THEORY
T a I + a M + a D + a (MD+DM) + a D2 + a (MD 2 +D 2 M)
0- 1- 2- 3 -- -- 4- 5 -- - -
(1)
a. a. (trD, trD 2 , trD 3 , trMD, trMD 2 )
1 1
M V ® V (2)
-3 -3
aT
o if a~ t= ~ ( 3)
au.
1
E. + nu. 0
aE. 1 1
J
~~ ~0, 1
where: n = when i 2, 3 j 1,2,3,4,5, (4)
4, 5
In (4), the following notation is used for the arguments of the functions
u.:
1
l I 2 l I 3
El trD E2 tr D2 E3 tr D3,
l I 2 (5)
E4 trMD E5 tr MD 2 ,
1
u P(p,q,r,s) ul Q(p,q,r,s) u2 R(p,q,r,s,)
0 E2
( 6)
1 1 1
u3
E2
s(p,q,r,s); u 4 E2 T(p,q,r,s,);u 5 E2 u(p,q,r,s,)
2 2
El E3 E4 E5
p q r = s = (7)
E E E2 E2
2 2
The tensor form (8) is equivalent to the system of five scalar rela-
tions which specify the following five invariants of the functional basis
of T in the transversely isotropic case:
_.... - -
and MT 2 +T 2 M from the tensor relation (8) and then taking the traces of
the obtained expressions. Some reductions are obtained by the application
of the generalization for 3 three-dimensional second order tensors of the
Cayley-Hamilton theorem, due to RIVLIN [4].
The Cayley-Hamilton theorem for a three-dimensional second order
tensor is given by:
0 (10)
(trAtrBtrC-trAtrBC-trCtrAB-trBtrAC+trABC+trCBA)I - 0
0 (13)
Singular cases of yield loci correspond to the case where the dimension-
less kinematic variables (7) arenot independent, i.e. when kinematic
constraints are imposed on the motion of the solid.
In the general form (13) of the yield condition, the independent
variables involve, in addition to the isotropic invariants trT, trT 2 ,
trT 3 , two mixed invariants tr~~· tr~~ 2 ,which account for the anisotropic
character of the material behavior.They specify the orientation of the
principal directions of the stress tensor with respect to the privilegcd
axis of the transversely isotropic material. In an arbitrary privileged
orthonormal frame (·, v ), the mixed invariants are expressed by:
-3
trMT T (14)
33
(15)
D
~ I + ~ M+ ~ T + ~ (MT+TM) + ~ T2 + ~ (MT 2 +T 2 M)
0- 1- 2- 3 -- -- 4 5 -- - -
(16)
~.
1
= ~.(p,q,r,s)
1
T
(18)
a.
1
A similar procedure leads to the same general invariant form (13) for
the yield criterion and to the following invariant quasi-linear form f~r
D
= ~0-I + ~
1-
M+ ~
2-
T + ~
3
(MT+TM)
-- --
(19)
~. ~.(p,q,r,s
1 1
In this case, the explicit calculation necessary to derive (13) and (19)
from the quasi-linear form (1) and the homogeneity condition (3) is given
in [3].
0
p
p p p p
~ p
e:v
If we assume that stresses are prescribed during the test, the stress
tensor T is defined in the axes of the specimens by:
p 0 oI
T 0 p
I
0 0 00 n (v ,e ,e )
(20)
-1 -2 -3
For the stress (20), we obtain the following components for the rate of
deformation tensor, as obtained from the law (16) and expressed in the
principal frame (v ,e ,e ) of the stress tensor:
-1 -2 -3
D
11
D2 2
- - - = '!' + '!' 1 sin 2 8 + ('!' 2 +2'!' 3 sin 2 8)p + ('!' +2'!' sin 2 8)p 2
0 4 5
/trD 2
(21)
D
33
- - - = '!' + '!' cos 2 8 + ('!' +2'!' cos 2 8)o + ('!' +2'!' cos 2 8)o 2
0 1 2 3 n 4 5 n
/trD 2
D
23
['!'
1
+'!' (a +p) + '!' (o 2 +p 2 )]sin8cos8
3 n 5 n
D
3 1
D
l 2
o.
Vielding and Failure of Transversely Isotropie Solids 75
One can conclude that the principal directions of the stress tensor T
and the rate of deformation tensor D coincide if, and only if:
[1 1 + ' 3
(o n +p) + ' s (o n2 +p 2 )]sin8cos8 0 (22)
Thus, for the so-called "in-axis" tests, i.e. when the privileged direct-
ion v ofthe material coincide with one of the principal directions of
-3
the stress tensor (8 = 0 or TI/2), the principal directions of T and D
are the same. For "off-axis" tests (8 f. 0 and TI/2), the principal direct-
ions of T and D cannot coincide in general.
In fact, in a standard compression test under confining pressure,
mixed boundary conditions are prescribed, since D rather than T = o
33 33 n
is enforced in the axial direction e of the specimens. Thus, on the ends
-3
of the specimens, the principal directions of the imposed rate of defor-
mation tensor are the axes of the specimens. On the lateral sides, the
principal stresses T 11 = T 22 =p are prescribed, without shear stress. Thus,
on the lateral sides of the specimens, the principal directions of the im-
posed stress tensor are also the axes of the specimens. Now, if the deve-
loped stress and kinematic fields were homogeneous, the axes (v-1 , e-2 , e-3 )
of the specimens would be the principal diructions of both the stress ten-
sor and the rate of deformation tensor. As this is impossible in general
for off-axis tests, the developed stress and kinematic fields in specimens
of transversely isotropic material for 0 < 8 < TI/2 are heterogeneous; in
fact, a shear stress T 2 3 f. 0 appears on the contact of the specimens with
the platens and the specimens tend to assume a S-shape.
Experimental evidence of the S-shaped deformation of off-axis speci-
mens of transversely isotropic materials subjected to standard compres-
sions under confining pressure is given in ([2], [5], [6]) for a strati-
fied rock (diatomite). Fig.2 shows an undeformed specimen,which was drill-
ed out at the orientation 8 = 60°; Fig.2b shows the deformed shape of
the specimen after axial compression (o 13.2 MPa) under a confining
n
pressure p = 6 MPa.
The same reservations apply to other tests and particularly to the
standard uniaxial tensile tests, where the ends of the specimens are
76 J.P. Boehler
a) b)
Fig.2 - Nonhomogeneity of deformation in compression tests
under confining pressure on a stratified rock (after
[2]); a) specimen before loading ; b) specimen after
unloading.
T pl (23)
Dl l D2 2
'jl + 'jl p + 'jl p2'
0 2 4
ltrD 2 ltrD 2
(24)
D3 3
'jl + 'jl l + ('ji 2 +2'jl 3 )p + ('ji 4 +2'jl 5 )p 2
0
/trD 2
D2 3 D3 l Dl 2
0
trD
3'jl + 'jl + (3'j1 +2'jl )p + (3'j1 +2'jl )p 2 (25)
0 l 2 3 4 5
One can conclude that the principal directions of the rate of deformation
tensor D coincide with the privileged directions (·, ·, v) of the trans-
- 3
versely isotropic material. Since D3 3 'i Dl l D22 , the deviatoric part of
D is different from zero; this constitutes an essential difference from
the behavior under hydrostatic pressure of isotropic materials.
An experimental study has been developed in [12] in order to show
that yieldingand hardening of a soft rock occur under hydrostatic pressure
and that the directional properties of an anisotropic material result in
a non-isotropic deformation under isotropic pressures. The material used
was diatomite, a soft stratified rock. Cylindrical specimens were drilled
out from blocks in seven different directions, oriented at 8 0°, 15°,
30°, 45°, 60°, 75o and 90° with respect to the privileged axis v of the
- 3
78 J.P. Boehler
MPar------.------,------,------,
"10
~*~
L1 L2
_j
30
~
L
•l b)
4.1 Introduction
0 (26)
s T - (.!.
3
trT)I
- -
(27)
2
-(o -p)2 2
-(o -p)3
trT 2p+o
n 3 n 9 n
(28)
1
trMS -(o -p)(3cos 2 8-l)
3 n
1 l I2
-trS 2 + ~ trMS tr S2 E =1 if 0 -p ) Ü
3 - 16 n
Only three invariants remain independent and the general form of the
yield condition for triaxial tests is reduced to:
p= 0
P= 25 MPa
p= 50 MPa
p=75MPa
and 90°), failure occurs in planes crossing the mats; this mode of failure
will be called "across mode". Simple compression parallel to the mats
(8 = 90°, p = 0) is a particular case, where failure is produced by an
extension parallel to the mats; when the confining pressure increases,
this particular failure mode is replaced by the "across mode". These
phenomena are visible in Fig.6, where one deformed specimen, after fai-
lure, is presented for each orientation and each confining pressure; the
visible sides are those with inclined mat planes.
Failure modes in simple tensile tests are quite different from that
observed in compression(Fig.7). For the orientations 8 = 0 to 60°, fai-
lure occurs by an extension parallel to the mats, which is progressively
accompanied by a shear stress parallel to the strata when 8 increases
from 1S to 60°. For 8 = 90°, failure occurs in several planes crossing
the mats' planes. At 8 = 7S 0 , mixed phenomena are observed.
For the compression tests under confining pressures, the obtained
values for the limit over-stresses a -p are presented in Table I. What-
n
ever the confining pressure, a marked drop in the directional strength
is observed at the transitions from 8 = 30° to 8 = 4S 0 and from 8= 90°
to 8 = 7S 0 , which correspond to the transitions from the "across mode"
to the "parallel mode" of failure. Thus, the anisotropy of the triaxial
strengths is of the discontinuous type [1]. In simple tensile tests, the
p(MPa) 8 = 0 8 = 1S 0
8 = 30° 8 = 4S 0
8 = 60° 8 = 7S 0
8 = 90°
0 17 21 20 38 50 84 230 traction
0 480 410 340 160 1SS 208 325 compression
2S SOS 4S8 37S 180 16S 22S 383 II
P-
L -------
lan
--------
'""_j
Fig.8 - Failure in "parallel mode" (after [13]).
io -picos8sin8 = ltrMS 2 - tr 2 MS
n -- --
l
0 (o -p)cos 2 8 + p trMS + -trT (32)
3 3 n 3
3p trT - I .\rs 2
2
From relations (32), one can conclude that the proposed criterion (31)
is a particular form of the general criterion (30). Fina1ly, the proposed
criterion (31) admits the following form:
trMS
0 < on-p /2.trS 2 cos 2 8 l+
3
I ~trS 2
2
2
(35)
Jp trT /2.trS
2
2 •
Vielding and Failure of Transversely Isotropie Solids 87
From relation (35), one can conclude that the proposed criterion (34)
is a particular form of the general criterion (30). Griented simple com-
pression tests allow the values of the material constants b and c and
0 0
the constant k of the criterion to be determined.The six other material
constants account for the influence of the confining pressure.
When applied to simple tensile tests, the criterion (34) admits the
reduced form:
(J (36)
n
This form is more general than the corresponding reduced form of the
criterion (33) and will be used for the prediction of the directional
simple tensile strengths. The three constants k, b and c must be recFti-
o 0
culated, because the failure modes in simple tensile tests are different
from the "across mode".
u.--••25
350r---~--~~~hr~----~~~~~~~
250
100
50
fying the failure strengths for the "across mode" predicts limit-stresses
which are higher than those predicted by the criterion (33) specific for
failure in "parallel mode". Finally, the lower envelope of the two pro-
posed criteria furnishes correct predictions for the values of failure
strengths, as well as for the orientations corresponding to a given
failure mode.
Vielding and Failure of Transversely Isotropie So Iids 89
5. SIMPLILFIED THEORY
T AT
_1 _1
T <P
o-I +
<P D
1-
or T= <P A I +
0::;
<P A D
1:::
(39)
0 (40)
90 J.P. Boehler
E trD E
l 2
1
cj>o = A(p,q) cj>l =E B(p,q) (41)
2
Introducing the solutions (42) into the representation (39), the inva-
riant form for the plastic constitutive law of a transversely isotropic
solid is obtained:
T =Al + JL D
E2 -
(43)
A A(p,q) B B(p,q)
trT 3A + pB
(44)
where M v ® v .
-3 -3
T
(47)
+ l $ D + (l - l)$ (MD+DM),
y 1- ß y 1 -- --
where $ and $ are specified in (41). The flow law is obtained by inver-
o 1
sion of (47); this is equivalent to the inversion of (43) and the intro-
duction of (46):
D (a+y-2ß)trMT
1 I 2
A I + M+ 1 T + b (MT+TM) . (48)
B - B B - B -- --
tr D2
The flow law (48) of the simplified theory camplies with the general
theory in its quasi-linear form (18). Finally, the simplified theory
(39) is a particular case of (1).
In any (•, •, v ) privileged orthonormal frame, the component form
- 3
of the flow law (48) is given by:
-A+yT 11 yT12 ßT 1 3
B
D yT 12 -A+yT 2 2 ßT23 (49)
tr1r2D2-
ßT 13 ßT 2 3 -A+aT
33
It is easy to see that the flow law (49) of the simplified theory pos-
sesses the same essential features as the flow law (16) of the general
theory, i.e.:
92 J.P. Boehler
A.
D11 3(2rT 11 rT 2 2 - aT 33 ) D2 3 A.ßT 23
A (51)
D22 -( -rT 11 + 2rT 22 - aT 3 3 ) D31 AßT 31
3
A -rT
D33 -( 2 2 + 2aT 3 3 ) D1 2 AaT 1 2
3 11 - rT
tr1/ZD2
where: A
B
D D D 0 (52)
11 31 12
T
11
1 (T
2 22
+ ~r T
33
) T
31
T
12
= 0 (53)
A
D33 = -2 (-rT 22 + aT 33 ) D
1 2
AßT 1 2 . (54)
Yielding and Failure of Transversely Isotropie Solids 93
(55)
s T - (l trT)I (56)
3 - -
(57)
R 2 (58)
8
2a 2+ ar - 3ß 2 (59)
2y2 + ay - 3ß2
In the case where y > a, depending upon the values of the material plastic
anisotropy constants a, ß and y, the strength variation is of the type
shown in Figure 11. The strength variation:
1
a) admits a maximum for 8
max < .:!!.4 if ß2 < -(2a
3
2 + ay)
1
b) admits a minimum for 8
min > .:!!.4 if ß2 > -(2y2
3
+ ay)
1 1
c) decreases monot ically if -(2a 2 + ay) < ß2 < -(2y2 + ay)'
3 3
(60)
ß (61)
In the case of plane flow, taking into account relations (53), the
yield condition (55) results in:
(62)
where an and crt are the principal stresses in the plane of deformation.
The analysis of the curves (62), as well as the analysis of the criterion
(55) in plane strain and axisymmetric three-dimensiona!. stresses is given
in [1], [20].
For materials with internal friction, a simple criterion can be
developed starting with an invariant form of the Drucker-Prager crite-
rion for isotropic solids. In the simplified theory, the generalization
for transversely isotropic materials of the Drucker-Prager criterion
takes the form, [17], [18]:
=c (63)
96 J.P. Boehler
c
(65)
where the + sign before the constant a holds for the stren~th in simple
compression and the - sign holds for the strengthin simple traction.
Similar curves to thoseof Figure 11 are obtained in this case, except
that their form is now influenced by the constant a. The respective dis-
cussion, as well as the analysis of the criterion (63) in plane flow,
plane stress and axisymmetric three-dimensional stres~es are given in
[1], [20]. Camparisans with experiments are discussed in [2].
REFERENCES
J.P. Boehler
University of Grenoble, France
1. INTRODUCTION
1
f(trT, /trS 2 , ltrS 3 ) 0 s T-(- trT)I ( 1)
- 3 - -
(2)
where
A(P) > 0
; { A(~) )
1: hardening
0 ( 3)
where the ~. are structural tensors, taking into account the anisotropic
-1
material behavior [4], [5]. In order to derive the invariant form of (3),
the theory of representation of anisotropic tensor functions [4], [21] is
very useful, for it indicates the type and the minimal number of basic
invariants involved. By application of this theory, we obtain the general
invariant form of the initial criterion:
102 J.P. Boehler
-- 3 __
f(trT, /trS 2 , /trS 3 , J 1 , J 2 , • •• , Jn) 0 (4)
where the Jk are mixed invariants of the stress deviator S and the struc-
tural tensors ~. and arehomogeneaus of order 1 with respect toS. The
-1
type of anisotropy of the material is specified by the invariance group
of the Jk [4].
For example, the initial yield criterion of an orthotropic material
is given by:
~ 3~ ( 1 1 r:-;;-:; 1 3r.-;;-;- 1 )
8 3 IJk
f ( trT, vtrs~, vtrs~, Jk
)
= 0 => f Itr!, ~tr~ 4 • ~tr~~, 0 (6)
Such a transformation does not modify the invariance group of the Jk. Thus,
the initial type of anisotropy is preserved in the subsequent criterion.
Finally, for initia1ly isotropic or anisotropic materials, the cla-
ssical concept of isotropic hardening is equivalent to a homothetic trans-
formation of the stress tensor:
For this type of isotropic hardening rule, we suggest the name "homothetic
isotropic hardening" (see also Section 2.5 below).
dening rules were proposed [18], [25]. The classical concept of kinematic
hardening assumes that the yield surface undergoes a rigid translation
C(P) in the stress space, where C(P) is a symmetric second order tensor.
This transformation is equivalent to the translation -C(P) ofthe stress
tensor:
1-
By the isotropic part 3c(~)~, which induces a translation of T parallel
to the isotropic axis of the stress space, the initial criterion (1) is
transformed into:
f(trT, ftrS 2 , ~) = 0 s - 1 -
T -(-trT)I
- 3 - -
(11)
trT = trT
The analysis of Section 2.2 shows that the anisotropic part of kine-
matic hardening is a very restrictive type of orthotropic hardening. In
fact, kinematic hardening cannot take into account most observed phenomena
of anisotropic hardening [9]. Different, more elaborate anisotropic har-
dening rules have been proposed [1,2,22,24], in order to account for change
in shape and size and for rotation and translation of the initial yield
criterion. The aforementioned rules are based on anisotropic linear trans-
formations of the stress tensor. A more general approach, based on the
Formulation of Isotropie and Anisotropie Hardening 105
(16)
106 J.P. Boehler
translational
isotropic
hardening
h ;::::_ J initial criterion and
I homothetic isotropic hardening
tr T
(17)
/
' _jsubsequent
....----1
criterion
homothetic
isotropic hardening
initial criterion
and
translational isotropic
hardening
/
Fig.2 - Von Mises materials subjected to isotropic plastic plane strain
2.5 Conclusions
Thus, there appears a need for a more general concept and a rational
formulation of isotropic and anisotropic hardening.
T a I + a S + a S2
o- 1- 2-
1
(19)
a. a.(trT, /trS 2 , ltrS 3 , E)
1 1 - -
3 __
trT = kE(trT,
1 -
/trS 2 , /trS 3
3 __
/trS2 kE(trT, /trS 2 , ltrS 3 ) (20)
2 -
3__ 3 __
/trS 3 kE(trT, /trS 2 , /trS 3 )
3 -
The last property shows that the proposed concept of isotropic har-
dening is the most general. The demonstration of this property will be
Formulation of Isotropie and Anisotropie Hardening 111
given in ref.6.
0 (22)
1
T => T T(T, ~ .• P); P-(-trP)I E :f. 0 (23)
-1 - - 3 - -
where the K~ are mixed invariants of !• ~i and P, which take into account
the initial anisotropy, the induced anisotropy and the coupling between
them [8].
Furtherdetails and developements will be presented in ref. 6. The
main properties of the general rule (24) are the following:
4. EXAMPLES
4.1 Introduction
(25)
Consider first a material obeying in its initial state the von Mises
criterion. In this case, the coefficients involved in (25) are subjected
to the initial conditions: C(O) = A(O) = 0, A(O) = 1, so thatfor P = 0, we
have T(O) S. Employing (25) and rule (21) of isotropic hardening, the
114 J.P. Boehler
(26)
In (26), the coefficient D(E) and A(E) account for the change of size of
the criterion, C(E) describes the translation parallel to the isotropic
axis and A(E) stands for the change in shape. Thus, the very simple stress
transformation rule (25) and the application of the von Mises criterion
provide a hardening rule involving the essential features of our general
concept of isotropic hardening: the initial criterion can change in size,
shape and position along the isotropic axis.
The functions A(E), C(E), D(E) and A(E) must be further specified
from physical observations or (and) phenomenological considerations. The
purpose of this paper is to show the possibilities of the proposed rule.
Thus, we consider the following simple case:
E ~ 0 /3 c
(27)
/(a+2d)E /l+(a-d)E
1
where a, c, d, ~ arematerial constants. For E 0, we obtain the von Mises
criterion:
(29)
+ 2(- -31 + dE)(o l o 2 +o 2 o 3 +o 3 o l ) = 2k 2 (l+~E)
Employing (25) and rule (21) of isotropic hardening, the subsequent cri-
terion, when developed with respect to the principal stresses, is given
by:
,\ = 1 c - 0 A= D=k(E),
(32)
with 0 < k(E) ~ 1, k(O) = 1
}
T IJI M + IJI M + IJI T
1-1 2-2 3-
(34)
'f'. 'f'. ( trM T, trM T, trT 2 , trM P, trM P)
1 1 -1- -2- -1- -2-
where ~1 = :1® : 1 ' ~2 = :2® :2; (v-1 , v-2 ) being the principal directions
of P in the plane of stress.
When restricted to a first-order polynomial transformation, the in-
variant form (34) can be written in the form:
valued functions of trM P E and trM P E ' E and E2 are the princi-
-1- 1 -2- 2 ' 1
pal values of P.
Formulation of Isotropie and Anisotropie Hardening 117
+ 2[C 1 D1 + C2 (D 2 +F)]trM-2-T +
+ (B 21 +B 22 +2FB 1 )tr 2M T+
-1-
(36)
+ (D 21 +D 22 +2FD 2 )tr 2M T+
-2-
+ Fz trT 2 =2 k2,)_2
p 0 => c1 c 2
B
l
D
l
0 B
2
F 1 (37)
We consider the following very simple case for the stress transfor-
mation (35):
C (P)
1 -
CE
1
c ( p)
2 -
CE
2
=> T CE M +
1- 1
CE M + T
2- 2
} (39)
(40)
where S and E are respectively the stress and plastic strain deviators; c
and A are material constants.
The rule (41) involves a kinematic part, which is taken into account
by the following stress transformation:
s => s S - cE (42)
(43)
M + M - (t: +t: )M
}
E E E
1-1 2-2 1 2 - 3
(44)
M V ® V M V ® V
-2
M V ®V
- 1 - 1 - l - 2 -2 - 3 - 3 - 3
When applied toS, the proposed anisotropic hardening rule (36), with
C1 - C2 = 0, can be written in the form:
5. CONCLUSIONS
REFERENCES
16. MURAKAMI S. and OHNO N., A Continuum Theory of Creep and Creep Damage.
In: Creep in Structures, Proc. IUTAM Symp. (Leicester, Sept. 1980),
(Ed. A.R.S. Panter and D.R. Hayhurst), Springer, Berlin, 1981, 422-
444.
17. MURAKAMI, S. and SAWCZUK A., On Description of Rate-Independent Beha-
vior for Prestrained Solids, Arch. Mech., 32 (1979) : 251-264.
18. FRAGER, W., A New Method of Analyzing StressesandStrains in Work-
hardening Plastic Solids, J. Appl. Mech., 23 (4) (1956) : 493-496.
19. REES, D.W.A., A Survey of Hardening in Metallic Materials. In: Failure
Criteria of Structured Media, Proc. CNRS Int. Coll. 351 (Villard-de-
Lans, June 1983), in press.
20. RIVLIN, R.S., Further Remarks on the Stress-Deformation Relations for
Isotropie Materials, J. Rat. Mech. Anal. 4 (1955) : 681-702.
21. SPENCER, A.J.M., The Formulation of Constitutive Equations for Ani-
sotropie Solids. In: Mechanical Behavior of Anisotropie Solids,
Proc. Euromech. Coll. 115-CNRS Int. Coll. 295 (Villard-de-Lans, June
1979), (Ed. J.P. Boehler), Editions du CNRS (Paris) and Martinus
Nijhoff, The Hague, 1982, 3-26.
22. SVENSSON, J.L., Anisotropy and the Bauschinger Effect in Cold Rolled
Aluminium, J. Mech. Engng Sei., 8 (1966) : 162-172.
23. TOUATI A., Comportement Mecanique des Sols Pulverulents sous Fortes
Contraintes, These de Docteur-Ingenieur, Ecole Nationale des Ponts
et Chaussees, Paris, 1982.
24. YOSHIMURA, Y., Hypothetical Theory of Anisotropy and the Bauschinger
Effet Due to PlasticStrain History, Aero. Res. Inst.,Tokyo Univer-
sity, 349 (1959) : 224-246.
25. ZIEGLER, H., A Modification to Prager's Hardening Rule, Q.J. Appl.
Math., 17 (55), (1959) : 55-65.
Chapter 7
J.P. Boehler
University of Grenoble, France
1. INTRODUCTION
2. CONSTITUTIVE RELATION
v1
e2
e1
01
v2
T D (1)
T F(D, P, M) (2)
where T, D and P are the plane parts of the stress, the kinematic and
the prestrain tensors respectively. The structural tensor M
accounts for the initial orthotropy of the material.
As a consequence of the Principle of Isotropy of Space, F is an
isotropic function with respect to its three arguments D, P and M. The
analysis of form (2) shows that the stress tensor T is an orthotropic
function with respect to the two mechanical arguments D and P, [9] (see
also Chapter 3).
An irreducible representation of the constitutive relation (2) is
given by (see Chapter 3):
T=a.I+a.M
1-
+ a.D
2-
+ a.P
3-
0- (3)
a.. = a..(trD, trMD, trP, trMP, trDP, trD 2 , trP 2 )
1 1 -
126 J.P. Boehler
M E GD E M E GD E
- 1 -1 -1 -2 -2 -2
(4)
M
1 I~ ~I(E -1
'E) -2
M
2
M
- 1
+ M
- 2
= I
P
-
= (2trM P- trP)M
-l- --l
+ (trP - trM P)
- -l--
I (6)
Employing expression (6), the invariants trDP and trP 2 involved in the
representation (3) can be expressed as single-valued functions of the
invariants trP, trD, tr M P and trMD:
-l- -l-
}
-l- - - 1- -1- -
(7)
trP 2 2(trMP
- 1-
- tr P)trM P
- - 1-
+ tr 2 P
T <PI+<PM +<PD
}
0- 1 - 1 2-
(8)
+ a. M +
}
T a. I a. D
0- l- 2-
(9)
a.. a.. ( trD, trMD, trD 2 )
1 1 -
Comparing the representations (8) and (9), we conclude that for the
initially isotropic material the prestrain induced anisotropy is an ortho-
tropy of privileged directions E , E . If the prestrain tensor P evolves
-1 -2 -
without rotation of its principal directions, the induced anisotropy re-
mains an orthotropy of privileged directions ~1 , ~2 ; the degree of the
128 J.P. Boehler
induced orthotropy changes with the values of the invariants trP, trMP.
- 1-
We consider now an initially orthotropic material and we first suppose
that the principal directions of ~ coincide with the initially privileged
directions (v , v ) of the material. In this case, the secend order tensor
- 1 - 2
M is identical to the structural tensor M. Introducing the conditions
-1
M
-1
=M and employing relations (6) and (7), the constitutive law (3)
reduces to:
T ß0 I- + ß 1-M + ß 2-D
I (10)
ßi ß.(trD, trMD, trD 2 , trP, tr~~)
1 - ~
Comparing representation (10) with thatfor a non-prestrained orthotro-
pic material (8), we conclude that the orthotropy with privileged direc-
tions (v , v ) is preserved ; the degree of this orthotropy changes with
-1 -2
the values of the invariants trP and trMP.
Consider now the general case, when the principal directions (E , E )
-1 -2
of the prestrain P do not coincide with the initially privileged direc-
tions (v , v ) of the material. Introducing the expressions (6) for P
- 1 - 2
and (7) for the invariants trDP and trP 2 in the representation (3), we
obtain the following equivalent representation:
a) T y I + y 1-M + y 2~1 + y 3-
D
0-
( 11)
b) yi = y.(trD,
1 -
trD 2 , trMD, trM D,
-1-
trP, trMP, trM P)
-1-
y = k(P)y (12)
2 - 1
m(P) M + k(P)M
- 1
(14)
-
The principal directions of the tensor m makes an angle 6(P) with the ini-
tially privileged directions (v , v-2 ), such that:
-1
k(P)
tan[26(P)] sin 2'i' (15)
1-k(P)
where 'i' is the angle between the principal directions of the prestrain
tensor P and the initially priveleged directions (v-1 , v-2 ). We further
suppose that the scalar-valued functions y.1 satisfy the indentities:
(17)
~i ~.(trD, trmD, trD 2 , trP, trrnP)
1 -
130 J.P. Boehler
3. PLASTIC BEHAVIOR
3.1 General invariant forms of the flow law and the yield criterion
aT
0 if - f. 0 (18)
aD
1J2 'I'I+'I'M+'I'T+'I'p
o- 1- 2- 3-
tr D2
(19)
'1'. = 'l'.(trT, trMT, trP, trMP, trTP, trT 2 , trp 2 )
1 1 -
The general invariant forms (19) and (20) in the plane stress case
are derived by a similar procedure than the corresponding forms in the
three-dimensional case (see Chapter 5).
The flow law (19) specifies the rate of deformation tensor only to
within a scalar multiplier. The relation between D and T is anisotropic.
Anisotropie Hardening of Rolled Sheet-Steel 131
A tr 2 T
1
+ A trT 2 + (A trMT + A trT)trMT + (A trT + A trMT)trTP
2 3 -- 4 - -- 5 - 6 -- --
= 2A 7 k 2 (21)
where the A. are polynomials of order two wit~ respect to the invariants
1/ 2 1
trP, tr P2 and trMP. In criterion (21), the terms A1 tr 2 T and A2 trT 2
account for the isotropic part of the initial criterion and the isotropic
hardening induced by ~· the terms A3 tr 2 MT and A4 trTtrMT account for the
orthotropic part of the initial criterion and the evolution of its degree
with P, the term A5 trTtrTP describes the new anisotropy inrluced by P and
the term A6 trMTtrTP stands for the coupling between thB initial aniso-
tropy and the induced anisotropy.
For the non-prestrainedmaterial (P = 0), criterion (21) reduces to
the yield condition for orthotropic materials submitted to plane stress
[4]. Thus, the proposed criterion constitutes a double generalization of
the isotropic Von Mises criterion generalization for initially orthotro-
pic materials generalization for the description of deformation induced
anisotropy.
132 J.P. Boehler
9=0°
2 a1
v'3k
/ '
,/ ' -2
'1-- ''
Fig. 3. - Anisotropie hardening of an initially isotropic metal
Dimensions in mm
v,
.;16r----,--- -----,-------- .
-.... 0
~ t t
360
~12r----r-v~----T-------~
. .....
4 ~--4~-------m=+----~~
..
.
..
.. u
...
u
C/1 C/1 C/1 r"
0 5 10 15 \ rolling direction
ro11ing reductlon
v,) principal axis of prestrain
.----
.I
I
GI
Cl)
"C
GI
I
~
1:
""
a a
Fig. 7 - New testing device for "off-axis"
anisotropic specimens
Anisotropie Hardening of Rolled Sheet-Steel 135
X X
- U = 500
-- - V: 15
- U = 400
·-V=20
---V: 25
u=Cst u = Cst
v=Cst v = Cst
Step : S ... m
S tep : 5 ... m
10 m m
t Om m
A 8
tal programm, we first analyze the flow law for the non-prestrained mate-
rial (P = 0). In the case of the associated flow law (normality rule),
criterion (21) results in the following expression:
D
2 + 2T
3 trT + d 0 trMT)I +(d 0 trT + c 0 trMT)M
(-- (22)
-- -
(23)
Introducing the expression (23) in the flow law (22), one obtains:
3d - 2
0
E: AE: 1 with (24)
2 2(3d
0
+ 3c 0 + 2)
Relation (24) is only valid for incipient flow. The evolution of this
relation during the prestraining can be obtained from the flow law asso-
ciated with criterion (21) through the normality rule. However, the pre-
strains applied in our tests are small ; thus, in a first approximation,
we can admit that relation (24) is valid during the entire prestraining
process. In the case of an incompressible material, the principal values of
the prestraining tensor P are then given by:
E
1
=E E
2
AE E
3
-(1 + A)E (25)
Making use of the expressions (25), the application of the proposed cri-
terion (21) to oriented simple tensile tests results in:
(26)
R~(E)
(~ + A E + B E + B E 2 ) cos 2 8 + ( c + c E + c E 2 ) cos 4 8
+ A2 E 2 ) + ( B0
3 1 1 2 01 2
where Ra is the simple tensile strength in the direction 8 and the coeffi-
Anisotropie Hardening of Rolled Sheet-Steel 137
:#~
42 hba 42 hba
(:13,25%
40 40
~~"%,.
36
~~2,00%
36 36 r.% : 1 - 0
2 - 0,75
34~ ~ 34 34
3 - 1,35
4 - 2,60
32
~ e' 32 e' 32 e'
L0~~15-7.30~45~67
0~7~5~9~0~ L0~~15~30~45~6~0~7~5-9~0~ L0~~15~3-0-4~5-6~0~75~9L0~
REFERENCES
0 • ;
l.J (R)
tensors A, B, C, ... , or A , by Ao ., B C .. , ... , or AiJ. ; of anti-
l.J
0 • ,
-R l.J l.J
symmetric tensors X, Y, ~, or ~, by X .. , Y .. , zo ., or
With these x(R).
l.J ij l.J l.J
second-order tensors we associate 3 x 3 matrices, also denoted by P, Q,
~R; A, B, C, A ; X, Y, Z, X , where, for example, P = (Po.). The use of
- - - - -R - - - -R - l.J
the same notation to denote tensors and the associated matrices of tensor
components should not generally cause confusion. When the meaning is not
clear from the context we shall state explicitly whether matrix or tensor
notation is being employed.
l.J J
with components P .. Q.k and PQT is the second-order tensor with components
l.J J
Pij"'kj,
,) o and the dyadic product u~v is the second-order tensor with
components uov .. The calculation of the components of such products is
1. J
most easily performed by matrix multiplications. Thus if u, v are
interpreted as 3 x 1 column matrices and P, Q as 3 x 3 matrices of tensor
-
components, then the scalar u.v is the single element of the matrix
T
product u v; the components of the vectors uP and Pu are the elements of
Isotropie Invariants and Tensor Functions 143
the matrix products ~T~ and ~~ respectively; the components of the tensor
P>2 are the elements of the matrix product PQ; and the components of u0v
T
are elements of the matrix product uv
The transformation
or MMT I. (4)
Here eS .. is the Kronecker delta symbol, I is the unit tensor or matrix (so
~J
that (cSij) = I) and MT is the transpose of M. Then also
If
-
u .. u .. (9)
lJ lJ
144 A.J.M. Spencer
U .. = M. M. U
~J ~r JS rs
u. (10}
~
p - - - -
f ( U. , V. , . . . , P .. , Q .. , . . . ) = (det M) f (u. , V. , . . . , P .. , Q .. , . . . )
~ ~ ~J ~J - ~ ~ ~J ~J
u.v. = M. U
p
M. V
q
= 6pq Up V q U V u.v.
~ ~ ~p ~q
PP ~ ~
P .. M. M. P 6 p p P ..
~~ ~p ~q pq pq pq PP ~~
Isotropie Invariants and Tensor Functions 145
Then Peano's theorem states that 1 with the possible exception of the
vth order determinants in which a typical column is comprised of
the elements A~P) 1 every polynomial invariant of the tensors can be
~
The results (c) 1 {d) and (e) (proofs in Gurevich [5] and Spencer
[6]) show that the study of algebraic invariants reduces to that of
homogeneaus polynomial invariants.
M(S)
ij
o~J
.. cos8+e l.J
.. kvk( 1 ) sinS+ (l-cos8)v~ 1 )v~ 1 )
~ J
0 0
cos e sin e
-sin 8 cos e
A reflection in the planes normal to ~1 is characterized by the
( l)
matrix ~1 (Rij ) 1 where
Isotropie Invariantsand Tensor Functions 147
( 1)
R ..
~J
-1 0 0
~I 0 1 0
0 0 1
0 0
-1 0
~I
0 -1
The groups {a), (b), (c) are continuous groups, with an infinite
number of elements. Finite sub-groups of the full orthogonal group are
also of interest. In particular crystal symmetries are represented by
thirty-two groups arranged in seven systems. These groups are known as
the crystal classes. The transformations which form these groups are not
listed here. They can be found in Green and Adkins [7], Spencer [6] or in
texts on crystal physics. We consider only two cases, which may also be
regarded as symmetries appropriate to fibre-reinforced materials
reinforced by two families of continuously distributed fibres.
{d) Prismatic symmetry. The transformation group consists of the
matrices
(Pl (Ql
U, U, u .u (P,Q 1,2, ... ,Ml. (lll
l l -P -Q
It can be shown that this set is complete, and that (lll form an integrity
basis under the full orthogonal group for the set of vectors. The result
is due to Cauchy; a proof based on Peano' s theorem is gi ven in \~ey l [ 4]
and Spencer [6].
(P f. Q f. R f. P) (l2l
and the integrity basis comprises the invariants (lll dnd (12). The
expressions (12) are relative invariants for the full orthogonal group.
We note the syzygies
u .u u .u u .u
-P -S -P -T -P -U
u .u u .u
-R -S -R -T '::R·':u
5. ISOTROPIC TENSORS
M. . M. . . . • 14 . . a. . . (14)
llJl l2J2 lllJll JlJ2 . . ·]lJ
( l) (2) (~)
J = a. , . U, U, ••• U, (15)
l 1 l 2 ..• l~ l 1 l 2 lll
where u~P)
l are componen t s o f ab so 1 u t e vec t ors. Th en
a. . . (16)
l 1 l 2 ••• l 11
(17)
Since
6.l.r 6.l.S
e l.J
.. ke rst 6. 6. (18)
Jr JS
6kr 6ks
(19)
152 A.J.M. Spencer
and since is invariant, it is also equal to (19) with each u. and Pjk
I
-
replaced by the corresponding u. and pjk
~
coefficients
( 1)
(a) 7Ti i = tr 2: 1 , (2la)
(b) (2lb)
(o) (22a)
( 1) ( 7)
where 7T. . , •.. , 7T. . are expressions of the form
~J ~J
7T .. or (23)
~J
p(S) In each term of I at most one expression of type (c) or (d) need
ji .
appear. R 11 R21 R 3 are integers (not necessarily different) chosen from
1121 ••• 1 V 1 and S 11 S 21 ••• are integers (not necessarily different)
chosen from l 1 2 1 ••• 1 A. Expressions of the form (21) areabsolute
invariants for the full and proper orthogonal groups; expressions of the
form (22) are absolute invariants for the proper orthogonal group and
relative invariants for the full orthogonal group.
The expressions (21) and (22) give an integrity basis for u and P
-R -S
but the basis is not minimal or even finite. It is therefore necessary
to eliminate the redundant members.
TI (24)
or 1 in components 1
TI ..
l]
Note that the trace of a matrix product is the trace of its transpese
T T T T
tr TI tr TI , tr PQ ... S = tr S ... Q P (26)
tr 'TT :: 0, (28)
same degree, and tr<:; 1-:; 2 ) = o, we say that !;I and !;2 are equiva~ent.
o. (29)
Hence
0 (30)
- (QR+RQ)
......... ............
tr P- (RP+PR)
.... ..........
tr Q- (PQ+QP)
.... ......... ............
tr R
.....--
- P (tr QR- tr Q tr R) - Q (tr RP- tr. .R tr P) - R(tr PQ- tr P
.......................
tr Q)
......... --
= 0. (31)
..... ..... -
tr P 3 S- tr P 2 S tr P + ~ tr PS{ (tr P)
..... ..... -- 2 - tr P 2 }
......
- tr s det P
.....
0 ,
or,
tr P 3 S ::: 0. (32)
tr (PQR+PRQ+QPR+QRP+RPQ+RQP) S - 0 .
.................. ................ ................. .................. --- --- .....
(33)
- --- -- -
tr (P 2 R+PRP+RP 2 ) S ::: 0 .
-
(34)
(35)
- - - --
- tr P 2 (Q 2 R+RQ 2 ) S + tr (Q 2 P 2 R+P 2 RQ 2 ) S .
- .............. - -- (36)
tr P (QPRQ) S - - tr P (Q 2 PR+PRQ 2 ) S
( 37)
( 38)
tr ( PT+:E~) ~9~~~ - 0 .
( 43)
(44)
Isotropie Invariants and Tensor Functions 157
Thus the trace of any matrix product of degree seven (or more) can be
expressed as a polynomial in traces of matrix products of degree less than
seven.
integrity basis for more than six tensors is obtained by combining the
bases for the tensors taken six at a time, in all possible combinations.
tr A, tr A2 , tr A3 •
TWo matriaes A, B
The following invariants satisfy the theorem:
tr ABC (DEF+DFE+EDF+EFD+FDE+FED) - 0 ,
Of these lSO relations, it can be shown that SO are independent; they may
be used to express SO of the invariants in terms of the remaining ten,
for example the ten listed in Table l.
between absolute and axial vectors. With each vector u we can associate
an anti-symmetric tensor U by the relations
(47)
tr TI ,
The same methods as were used to reduce the integrity basis for
symmetric tensors can now be used for symmetric and anti-symmetric tensors.
In addition to the results used earlier based on the Cayley-Hamilton
theorem for general matrices we have available the following results for
traces of matrix products involving anti-symmetric matrices X, Y, Z and
symmetric matrices A, B, C, ... (which may be the unit matrix I):
tr X = 0, tr XA 0, (48)
TABLE l
INVARIANTS OF VECTORS AND TENSORS - PROPER ORTHOGONAL GROUP
Rivlin and Smith [14] have proved that there are no other distinet matrix
identities exeept ( 31) and (52) for symmetrie and anti-symmetrie 3 x 3
matriees.
of the forms (2la,b) are absolute invariants under the full group, and
those of the forms (22a,b) are only relative invariants under the full
( 1) ( 2) ( 7) .
group, where now TI. . , TI .. , ... , TI. . are componen ts of ma tr1.x products
l.J l.J l.J
in symmetric matrices and in anti-symmetric matrices associated with
anti-symmetric tensors and axial vectors. General tensors can be dealt
with by decomposing them into sums of symmetric and anti-symmetric
tensors.
(54)
Then
(2) ( 2)
U, TI, . V, = tr (~+~) ~ 2 , u.TI .. u. = trATI 2 (55)
l. l.J J - - - l. l.J J --
f. :E', . M. M. F
1. l.J 1.r JS rs
-(p) - (q)
u. p .. M !1 P(q)
1. l.J im jm mn
(58)
.t-1. M. <!> (u ~p) ,P (q))
1.r JS rs J mn
of u~p)
and P(q); the concept of a form-invariant tensor is therefore a
J mn
generalization of that of a scalar invariant.
J 1 =v.f. 1 (59)
~ ~
Evidently
f. F .. (60)
~ ~]
J (2)
J1 I I J ( 1)
Cl. Cl.
J2 I I
Cl. Cl.
Cl. Cl.
f.
~
\' I 3J < 1 l ;av
L. Cl. Cl. i F ..
~J
=I (61)
Cl. Cl.
This gives canonical forms· for fi and Fij although 1 because of the
possible existence of syzygies 1 the expressions may contain redundant
terms even if J(a. 1 ) and J(a. 2 ) are e 1 emen t s o f an ~rre
· d uc~'bl e ~n
· t egr~' ty
F .. ~(F .. +F .. )
l.J l.J Jl.
J'2 V. v.F ..
l. J l.J
(p) p (q)
Then J'2 is an invariant of V.' ~, and quadratic in V.' and
l. mn' l.
F .. (62)
l.J
J"2 K .. F ..
l.J l.J
Methods analogous to (a) and (c), but not to (b), can be exployed
when F .. is anti-symmetric.
l.J
As a simple example let us use method {c) above to find the form-
invariant symmetric tensor function of a single symmetric tensor A. From
the tables, the invariants of A and K are
J~ I 1 tr ~ + I 2 tr ~ + I 3 tr ~ 2 ~ ,
then gives
or
(64)
(65)
K 2X .. W.. X .. (W .. -W .. )
~J ~J ~J ~J J~
+ 18 (~2~2~-~~2~2) • (66)
HEFERENCES
13. SPENCER, A.J.H., Isotropie integrity bases for vectors and second-
order tensors, Part II, A~ch.Rat.Mech.Anal. 18 (1965): 51-82
14. RIVLIN, R.S. and SMITH, G.F., On identities for 3 x 3 matrices, Rend.
di Matematica Se~ie VI, 8 (1975): 345-353
15. SMITH, G.F., On isotropic integrity bases, A~ch.Rat.Mech.Anal. 18
(1965) : 282-292
16. PIPKIN, A.C. and RIVLIN, R.S., The formu1ation oi constitutive
equations in continuum physics I, A~ch.Rat.Mech.Anal. 4 (1969): 129-144
17. SPENCER, A.J.M. and RIVLIN, R.S., The theory of matrix po1ynomia1s
and its app1ication to the mechanics of isotropic continua, A~ch.Rat.
Mech.Anal. 2 (1959): 309-336
18. SPENCER, A.J.l1. and RIVLIN, R.S., Further resu1ts in the theory of
matrix po1ynomia1s, A~ch.Rat.Mech.Anal. 4 (1960): 214~230
19. SPENCER, A.J .H., The invariants of six symmetric 3 x 3 matrices, A~ch.
Rat.Mech.Anal. 7 (1961): 64-77
Chapter 9
A.J.M. Spencer
The University of Nottingham, England
l. TRANSVERSE ISOTROPY
\ ( 1) ( 1) ( 1)
L ß Ö ~ Ö ~ •• • ~ V V ••• V
n n k1 1 k2 2 kq q m1 m2 mr
(a) (b)
(1)
( 1) ( 4) ( 1)
( c) (d) V. 1T .. V,
~ ~] J
where 1T~~) are expressions of the form (23) of Chapter 8. However 1 (1)
~]
2. ORTHOTROPIC SYMMETRY
(4)
and so we may discard one of the set ( 4) , say ;: 3 0;: 3 , in favour of the
unit tensor I and take the basic set of anisotropic tensors to be
trA, tr A 3 , (5)
2 2 2 2
shown that (~ 1 ~~ 2 ) , (~ 1 ~ ~ 2 ) and (~ 1 ~~ 2 ) (~ 1 ~ ~ 2 ) are reducible, and so
(5) constitutes the required integrity basis.
2
J = I 1 tr ~ + I 2 tr ~ + I 3 tr ~ ~ +I 4 ~ 1~~ 1
(6)
3. CRYSTAL SY~~TRIES
A ....
l.JJl.
Clearly, the existence of any index symmetries will reduce the number of
independent invariants.
6. 6. 6. 6,
l.p J.q J.r l.S
6. 6. 6. 6.
JP ]q ]r JS
0 ,
6kp 6kq 6kr 6ks
P .. A .. +e .. kuk+w6 .. ,
l.J l.J l.J l.J
where
A .. (5 independent components) ,
l.J
le P ( 3 componen ts) ,
2 rsk rs
w (l scalar),
Anisotropie Invariantsand Additional Results 179
1
A(;J'k) = 6(A .. k+A.k.+A .. k+A.k.+A. .. +A. .. ) .
~ ~ J ~JJ~ J ~ k~J kJ~
-A. n. - ••• - A . . n
(L .. kx-p)Jq ... s (~J ... kp)x-q ... s
Now add (r-l)A(. . k") to eaeh side, and note that for any pair
... ~J X, pq ... s
of indiees t and u
e e A.
Ai ... t ... u ... s- Ai ... u ... t ... s tuv xyv ~ ... x ... y ... s
We also denote
= e . . k) vq ... s
B.. k
~J ... q ... sw uvwA( U~J...
rA(ij ... kt)pq .•. s rA(.. "'P q •.. s +e.1pwB.] ... knJVq ... sw
1] ... kn)
+e. B. n + ... +e 0 B ..
JPW 1 ... kJVq ... sw "'PW 1] ... q ..• sw
Thus the tensor A symmetrized in its first r-1 indiees is the sum of A
symmetrized in its first r indiees, and tensors formed by eontraeting B
with the third order a1ternating tensor. By repeated app1ieation of this
proeess, A ean be expressed entire1y in terms of symmetrie tensors of
order n or 1ess, and isotropie tensors.
12A. 'k = 12A(. 'k) +4e.k c(.) +4e.k c(.) +6e .. B(k)
1J 1J 1 p JP J p 1p 1JP p
where
e A
prs rsk
, c.1p = e A
rsp (ri)s
0) ,
u e B A -A v = e C A -A •
p rsp rs prr rpr p rsp rs (pr)r rrp
In the most genera1 ease, Aijk has 27 independent eomponents, A(ijk) has
10, B(ij) has 6, C(ij) has 5, ui has 3, vi has 3. The deeomposition is
not unique.
A. 'k =
1J
A~ .k+~(A .6.k+A . k.+A k6 .. )
1J pp1 J PPJ 1 PP 1J
This defines the tensor A' whieh has the property 'k = 0. Then A' has
A~
11
seven independent eomponents. In genera1 a traee1ess symmetrie tensor of
ordern has 2n+1 independent eomponents.
Anisotropie Invariants and Additional Results 181
X = Mx
~
-n
= L
-n-n
s I
1 r21T
j
(m)
P(m) 2 1T tr!::'n (l-cos8) d8,
0
-n+l
= {(1-at-n )(1-at
-1
~(a,t,n) ) .•. (1-at )(1-a)(l-at) ...
n }-1
... (1-atn-1 ) (1-at)
n
II {(1-ats)- 1 } .
s=-n
Hence
00
It follows that
1 r2 n
00
i8
L P(m)am 2 1TJ 0 ~(a,e ,n)(l-cos8)d8
m=O
n
= _!_ J21T II {(l-aeisB)- 1 }(1-cos8) d8,
21T 0 s=-n
1 r2n J.1 i8
= 2TI j0 TI
r=1
{<j)(a
r
1 e 1 n ) } (1- cos 8) d8
r
d8.
00
I P ( m) u
m = -
1 jr2n (1-cos8) d8
-----'-:;::------=-=-=---='-'----==---:-~
2TI i8 -i8
0 ' 0 (1-ue ) (1-u) (1-ue )
00
2 (1 8) d8
\ P(m)am = __!_ J 'TT --~~--~--~c~o~s~~~~------~~
L. 2TI 2i8 i8 -i8 -2i8
0 0 (1-ae ) (1-ae ) (1-a) (1-ae ) (1-ae )
= J.:_
2TI
J2TI --~~----(~1_-~c_o_s_8~)~d~8~-----~
i8 -i8 i8 -i8
0 0 (1-ue ) (1-u) (1-ue ) (1-ve ) (1-v) (1-ve )
184 A.J .M. Spencer
1 r2TT ( 1- cos 8) d8
2 Tl J i8 -i8 2i8 i8 -i8 -2i8
0 (1-ue ) (1-u) (1-ue ) (1-ae ) (1-ae ) (1-a) (1-ae ) (l-ae )
s~ ~ S k-8+n ?k-8+nl .
s~ k-8+n 1
REFERENCES
15 (1957) : 309-314
2. SMITH 1 G.F. 1 On transverse1y isotropic functions of vectors 1 symmetric
186 A.J .M. Spencer
A.J.M. Spencer
The University of Nottingham, England
1. KINEI-1ATIC CONSTRAINTS
X= x(X,t), or X.
l
F.
1R
= ax. ;ax
1 R
.
We employ the finite strain tensors C and B, with components CRS and B ..
l]
respectively, where
T
c F F, B FFT
dX. dX.
l
dX. dX.
-- l
-- ___2:. _ 2
CRS FiRFiS B .. F. F.
axR axs l] lR JR axR axR
ax. Clx
l +_2 -eS
.!. ( __
.J
E E,.
lJ 2 ax.J ax.l ij ·
- l. (-~- ___2
avo av OJ
Do~J 0 - 2 "oxo +"OXo .
J ~
Do 0 tr D 0 .
~~
Eo
~~
0 trE o.
( 0) ( 0)
A.a or, in components, A.ao a 3xo /3X ,
~~0 ' ~
aR FiR R ~ R
( 0) ( 0)
a a C
R S RS
( 0)
a 3xo/3X , 1 .
R ~ R
aoaoD.. 0.
~ J ~J
aoa OEO 0.
J l.J
0
1.
R -p~ Ro
l.J
0 = -po l.J
0 0 •
R Ta®a 1 Ro Taoao
l.J
0
1. J
Ro ODO 0
l.J l.J
0 1
1.
2. LINEAR ELASTICITY
W(E, a€>a)
and only three elastic constants ~T' ~L and ß are required. If the
material is inextensible in the direction a, then aEa = 0, and
W = ~ tr E 2 + 2 (~ -~ ) aE 2 a 1
T ~ L T ~~ ~
and only the shear moduli ~T and ~L are required to specify the material.
Incompressible material
S = 2~ E+ß(aEa)a@a+2(~ -~) (a0aE+Ea0a)
~ T- ~~~ ~ ~ L T - ~~ ~~ ~
S
~
= ,\ItrE+2~ E+2(~ -~)
~ ~ T~ L T
(a0aE+Ea0a).
- -~ -~ ~
In each case the appropriate reaction stress must be added to give the
total stress.
tr E 1 aEa 1
where cos 24; is the angle between the two fibre directions. !t is then
Straightforward to write down the most general function W which is a
quadratic in E formed from the above list. If the material is
Kinematic Constraints and Constitutive Equations 193
3. FINITE ELASTICITY
s .. (1)
lJ
tr C , detC
Incompressib~e materia~
I 3 '
the direction ~ 0 , we replace W by W+tT(I 4 -l) and proceed in the same wayo
s T- .!r tr T .
- 3_ -
tr S Ü I aSa = 0 1
and then
J 2 = aS 2 a 1 J 3 = tr ~ 3
(2)
and
196 A.J.M. Spencer
( 3)
with
tr S 0 1 aSa 0, bSb 0.
Then
J 1 I J 3 I
and
The most general yield function which is quadratic in the stress is (for
the case in which the fibre families are mechanically equivalent)
D D +D
-p -e
D
-p
= "{1 (1 11
~ --yS + --y---zj (a®aS+Sa0a)
kT - kL kT - -- -- -
} (5)
(
0 (S-a 0as-sa 0 a) Jl k2
TI J2 < k2LI
D
-p " { V (a0aS+Sa0a)
Jl k2
I J2
J2
k2
LI
k2
'"" "'"'"' ""'"' '""
T I LI
tr D = 0, aD a = 0,
-p -~p~
it follows that there are three such independent parameters, which may be
taken tobe ET, EL and E 3 , where
(6)
Then
Clk Clk
- -TE• + - -TE• ( 7)
Clt:T T dEL LI
(8)
Therefore, when k
T
and k
L
are functions of wp only,
k. T = 2~ak ;awp ,
T
kL = 2~dk /dW
L p
(10)
k
L
condition and flow rule 1 s is replaced by s-a so that the yield surface
1
tra = o 1 aaa 0.
J 3 I tr T- aTa.
Now Y1 and Y2 are yield stresses in tension in the fibre direction 1 and
equibiaxial tension in planes normal to ~~ respectively. In practice
they are much larger than the shear yield stresses k and k . To admit
T L
constraint hardening f is modified to take the form
f 1.
Kinernarie Constraints and Constitutive Equations 201
REFERENCES
J. Betten
Technical University Aachen, F .R. Germany
I . INTRODUCTION
F(a. a. o ) _ F(o .. ) ( I)
~p Jq pq ~J
sV _ tr o V V I, 2, 3 ( 2)
Invariants of Fourth-Order Tensors 205
J I - o (3a,b,c)
ll
respectively. In (3), i.e., (-1) 3-v Jv = oi 1[i 1] ... oiv[iv]' the operation
Jl SI ' (Sa,b,c)
J I' (6a,b)
(6c)
which are identical to the principal invariants (3a,b,c). Thus, the re-
presentations (4a,b) imply isotropy.
The elementary symmetric functions (6a,b,c) are the coefficients of
the characteristic polynomial P3 (A) of the indeterminate A
(7)
( 8)
F ( 9)
F ( a. a . o ; a . a . ak a 1 A ) =: F ( o .. ; A1. J. k 1 ) · ( I 0)
1p Jq pq 1p Jq r s pqrs lJ
y .. or y ( !Ia, b)
lJ a
where i,j,k,l = 1,2,3 or a,ß = 1,2, ... ,9, respectively. In(!!) two pos-
sible indicial notations are used. The transformation laws of the Car-
tesian tensor components can be written as
Y~'. a. a. Y or Y* (12a,b)
lJ lp Jq pq a 3 aßyß
and
208 ]. Betten
a. a. ak a 1 A or: A* a a A , (13a,b)
1p Jq r s pqrs aß a~ ßn ~Tl
where:
aac;; a. a.
lp Jq
( 15)
A A
0 .. or 0 ( 16a,b)
lJ a
A(O) (0)
where A. is a real scalar to be determined, and Aijkl - 6 ik6 j 1 or
aß -
- 6 Ao , whereas 0 .. or 0 are
aß are the components of the unit tensor lJ a
~
0 or 0 (17a,b)
Invariants of Fourth-Order Tensors 209
n
(0)
p (;\) - det(Aijkl - AAijkl)
n I J (A)A n-v
- ( 18)
v=O \)
~n which, as we can see from ( 17a,b), the first indexpair (ij) =a cha-
racterizes the rows and the second one (kl) =ß the columns of a n x n
matrix ~' generally n = 9.
Analogaus to (3), the principal invariants J ~n ( 18) can be deter-
v
mined by performing the operation of alternation [ 19]:
where (-l)n J 0 =I. The right hand side in (19) is equal to the sum of
a ll ( n}
v = v.'( nn_! V )'• · · 1 m~nors
pr~nc~pa · o f order \! ~ n, wh ere \! = I and
v = n lead to trA and det~, respectively.
Assuming we have the usual symmetry conditions
(20a)
or alternatively written as
then we can express the zero power tensor of fourth-order in ( 18) in the
following way:
A(-I) A ( 2 I)
ijpq pqkl
2
J4 - (54 + 85153
I 65251 + 35 22 - 654)/4!
(22)
where
SV - trA V - ( 23)
..
J2A lj J3 ..
.. _ olJ
olJ (24)
I
v=O
J
V
~ n-v)
b)'
:: 0
ßy
ß,y 1, 2, ... ,n ( 25)
n
pQ A( n-V)
I n-v ßy n '
( 26)
V= I
r>z 6 7
2
8
1 - JI - J2 JI - J3 - J3
+ + 21 11 2 I
2 - J2 - J3 + J 11 2 - J4 + J 11 3 + J2(J2 - J2)
I
3 - J3 - J4 + J 11 3 - J5 + J 11 4 + J/J2 - J2)
I
4 - J
4 - J5 + J 11 4 - J6 + JIJ5 + J4(J2 - J2)
1
5 - J
5
- J
6 + J1J5 JIJ6 + Js(Jz - J2)
I
6 - J6 1 11 6 J6(J2 - J2)
I
212 ]. Betten
From (22), (26) and table I, we see, that all invariants formed by
traces of tensor powers can be expressed by scalar polynomials in the
irreducible basic invariants s 1, ... ,s 6 defined in (23).
The scalar coefficients in (26) can be determined by using the re-
cursion relation [22,23]
p-n
(-l)n-I(J + I p-)lQ J) ( 27)
p-n+V ].1= 1 n-v ].1
which is valid for the arbitrary dimension n. For the symmetric tensor
(20) of rank four (n = 6) we can use (27) to find the scalar polynomials
listed in table I, and for a symmetric second-order tensor (n = 3) we
can calculate polynomials listed in [ 16].
2 n-1
v , "'v•Tv,
,.......,,......,
... ,v•T
,....., ,.....,
v (28)
11~ v] -
(29)
V "' I ,2, ... ,5,
B 0 11[ I J
pq qp - 2
B 0
(2)
pq qp
cpq 0 qp = 11[ I J
3
cpq 0(2)
qp - rri I J.
The isotropic special case can be expressed by the isotropic con-
stitutive tensor
( v)
a V 6 lJ
.. okl + bV (o.ko"l (30)
Aijkl 1 J + o.lo.k)
1 J
of power V"' 1,2, . . . . Th~n the simultaneaus invariants IT~v] in (29) are
equal to the principal invariant J 2 (~) for av "' - 1/2, bv "' 1/4 and
equal to the basic invariant s 2 (o) for a "' 0, b "' 1/2. Similarly, the
[v] - V 3 V
invariants 11 3 in (29) are equal to J 3 (~) - J 1 (~)/6 for av "' - 1/6,
bv "' 1/6 and equal to the cubic basic invariant s3 (o) for a "' 0, b
- V V
. . [v] _ (v) [v] _
1/2. Furthermore, the 1nvar1ants 11 1 = oijAijklokl or ITI* =
(v) (2)
- oijAijklokl are equal to s 1 (~) or s 2 (~), respectively, if 3av + 2bv
I.
It must be mentioned, that the irreducible invariants 11[lv] and 11[v]
I*
arenot contained in (29). Furthermore, a vector v of arbitrary dimen-
sion n has only one invariant, namely v 2 in (28).-However, in cantrast
to ~· a second-order tensor in three dimension has three invariants (3a,
214 J. Betten
[v] B .. A~~) o( 2 ) +
( v)
B .. arr 3 /Clo .. 2B .. o .kAk. o - tr BC + 2tr BOC
1J 1J 1J 1Jpq pq 1J J 1pq pq
2A .... - 2J I(~)
1J1J
2ll[vJ.
2 '
Clll[v]
3
B .. - - - = B.. A<.v.) o(2) +
2B A(v)
.. o J'k k.1pq 0 pq _ tr BC + 2tr BoD
1 J Clo .. 1J 1Jpq pq 1J
1J
where 1 .. _
1J
lnvariants of Fourth-Order Tensors 215
(3)
B .. o .. okk
l J J l
(31)
(0) (0)
instead of the spherical tensor AAijkl' where Aijkl 1s the zeropower
tensor defined in (21). Thus we take into consideration the characte-
ristic polynomial
(33)
y .. AX c5 .• + 2WX .. (34)
lJ rr lJ lJ
racterized by the two elastic constants A. and lJ, which are called the
LAME constants. Using the notation Ya = IaßXß, a,ß = 1,2, ... ,6, we see
that the isotropic tensor I in (34) can be represented in the matrix
form
2).J + A. A. A. 0 0 0
A. 2).J + A. A. 0 0 0
A. A. 2).J + A. 0 0 0
(35)
1 ijkl -
0 0 0 2).J 0 0
0 0 0 0 2).J 0
0 0 0 0 0 2lJ
(0) (38)
{ I, I, I, I, I, I} .
Aijkl
We see, in the special case (,\ = 0, ).J = 1/2) the isotropic tensor (35)
tends to the unit tensor (38). Then we find from (32) the simplified
characteristic polynomial ( 18). The determinant of (35) is
5
det(Iijkl) = 32(2).J + 3A)lJ ,
(0)
which tends to det(Aijkl) in the special case (,\ 0, ).J I I 2).
(39J
,p q
P( A, JJ) c( p,q )"JJ, p 0,1,2,3; q 0, I, ... ,6, ( 40)
p,q p + q ~ 6
( 4 I)
0.
218 ]. Betten
Then, the dyads ~ and Xin (I Ia) are coaxial and, instead of (39), the
following determinant is to expand:
(43)
I = KM
6 - I 3 '
Invariants of Fourth-Order Tensors 219
I I I - - I 6 [ 2 ( K I - L I ) + 3M I ] ,
(44)
(45)
We see, in the system (44) the invariants r 1 and r 15 are the trace and
determinant of the fourth-order tensor, respectively.
In the special case (A = 0) only the six invariants r 1, ... ,1 5 , r 15
are relevant, and they are identical to the six invariants in (22), if
furthermore B1 =
B11 = B111 = 0.
In the two dimensional case (i,j,k,l = 1,2) the linear transfor-
mation (I Ia) can be represented in the following matrix notation
220 J. Betten
ZAIII2)
C") (All II
C)
A1122
y22 A2211 A2222 2A2212 x22 (46a)
yl2 A1211 A1222 2A1212 x12
y ~I X (46b)
C)- C°
y22
yl2
+;>., ,, ;>., +
0
;>.,
Considering (46a,b) and (47), the eigenvalue problern (33) leads to the
characteristic polynomial
(SO)
A~. a. a. A ( 5 I)
lJ lp Jq pq
M .. - a. a. 6 0 .. (52)
lJ 1p JP ij lJ
(54)
( A. . - ,\ ( ) 6 .. ) a( ) . 0 . (55)
1J r lJ r J r1
( I)
with 3 eigenvalues \I):::: >- 1 , ... "\( 3 ) - >- 111 and 3 eigenvectors ni
I (3) III
- ni'' .. ,ni - ni where
(r)
n. (56)
1
1s the r-th row of the orthogonal matrix ~· Note, that there 1s no sum
on the bracketed index r in (54) and (55). The condition for (55) to
have non-trivial solutions (56) yields the characteristic equation
de t( A. . - ,\ 6 .. ) = 0.
1J 1J
We can also start from the quadratic form
F A .. x.x. (57)
1J 1 J
2
M 6 .. x.x. - X (58)
1J l J
and yields
(A .. - \6 .. h. 0. (60)
l J l J J l
( 62a)
(62b)
( 63)
<P ijkl
! !
Cl<P( 1111
.... )/Cla(.) 0. or: Cl<P('''')/Cla(') 0 .. (64)
1 u lU ljlj 1 u lJU
Cno sum on the bracketed indices i or j). Going on 1n this way, we find
more details in [27].
The other possibility to apply the LAGRANGE multiplier method is
illustrated in the following. Similarly to (57), we start from the sca-
lar-valued function
F (65)
L 0 ' (66a)
M (66b)
N = ( 6 J. k6 p l + 6 . 1 6 k) X. X.. Xk l - 2S3
J p lp lJ
(66c)
Thus, the quadratic form (65) of the dyad X 1s modified to the form
F - AL - wM - vN ( 6 7)
CltP/dX ! 0 (68)
rs rs
where Iijkl are the components of the fourth-order isotropic tensor de-
fined in (31). In the special case when V= 0 we find from (69) the el-
genvalue problern (33) which yields the characteristic equation (32).
Some more details concerning the system (69) can be found in [27]. Fur-
thermore, applications of the LAGRANGE multiplier method to tensors of
order six are discussed in [27], too.
lnvariants of Fourth-Order Tensors 225
REFERENCES
FORMULATION OF ANISOTROPie
CONSTITUTIVE EQUATIONS
J. Betten
Technical University Aachen, F.R. Germany
I . INTRODUCTION
d .. f. .(o ,w ,A ) ( I)
lJ lJ pq pq pqrs
d L: cp G (2)
0. ~o.
0.
where the G's are symmetric tensor generators of rank two involving the
argument tensors o, w, A. The coefficients cp 0. in (2) are scalar-valued
~ ~ ~
d .. ( 3)
lJ
0 I 2
where Hijkl' Hijkl' Hijkl are the components of fourth-order tensor-
(4a,b)
d .. = f. . ( o ,w) w..
0
g .. ( o, w) , (Sa,b)
1J 1J - - 1J 1J - -
0 .. (6)
1J
s .. Ob)
1J
V. V.
s .. 1 J ( 8)
1J w. w.
1 J
x,
ol bl cI
d 2 S. - E .. k ( dx 3 ) . (dx 1\!2
1 1J J
(9)
d 3S. - E .. k ( dx 1). (dx 2 \!2
1 1J J
2 3 4
d 1S. + d s.1 + d s.1 + d s.1 0. ( 10)
1 1
where aijk- a Eijk' ßijk =ß Eijk' etc. are total skew-symmetric tensor
of rank three, which have the essential cornponents a 1 ~ 3 = a, B123 = ß,
etc., respectively. We see, the vectors d J~Si, ... , d 4 ~s.1 in Fig. lb dif-
fer frorn the corresponding vectors of the undarnaged configuration (Fig.
Ja) only in lengths. It should also be noted that the vector surn is not
equal to the zero vector:
~ ~
d Js~. 2~ J~
+ d S. + d S. + d S.
4~ 1
~ 0. ( 12)
1 1 1 1 1
where, by analogy to ( 10), the vector sum is equal to the zero vector:
The three area vectors d 1Si' ... , d 3 Si in ( 13) are identical to the cor-
4A
responding vectors of the damaged configuration. The fourth vector d S.,
4 1
4~
having the same magnitude as d S., differs from d S.1 not only in length,
1
4A 4
but also in its direction. Therefore, the vectors d Si and d Si are con-
nected by a linear operator ~ of rank two (second-order tensor):
( 15)
\)!, E • ( 17)
lr qk
( 18)
,~, I a ( 19)
'+'ijk = ei ljk +
w1J
.. k ( 20)
which is similar tothat in (7b). Finally, we can use (21) to deduce the
diagonal form:
Inserting (20) into (22a) and replacing o1 j ::: 1ej, etc. we see, that the
secend rank tensor of continuity can be decomposed in terms of dyadics
formed from the base vectors:
I I 2 2 3 3
\)! .. a( .!: ® e) .. + ß( e ® e) .. + y( ~. 0 e) .. ( 22b)
1J - 1J - lJ - - 1J
Formulation of Anisotropie Constitutive Equations 235
0 y 0/
/
/
0/ 0
/ 0 0 -ß
/ /
/ 0 a 0 0
/
/
0 0/
/
0 /
"I i =31
/
0 0 /0 / 0 ß 0
/ -/
0 / 1i =zl
/
0 /0 a / 0 0 y
/ /
/ /
/0 -a 0 Y'li=ll
Fig. 2 Dual Tensor of Continuity
W·lJ'k = W E.k Ö.
J r 1r
= W E.lJ'k <=>w.
1r
= W Er J.kE lJ
.. k/2 = \),' Ö.
1r
(23)
w.. -
lJ
eS .• -
lJ
wlJ
.. (24a,b)
236 ]. Betten
p. = O .. n. ( 26)
1 J1 J
was derived from the principle of equilibrium, where p., n. are the com-
1 1
ponents of the stress vector p and the unit vector n, respectively.
Pi= Pi 4>1nl
In the same way we can also find the corresponding relation for a
damaged continuum:
( 27)
dP.
1
p.dS - p.ds
1 1
dP.1 (28)
Thus, after comparing (26) to (27) we can finally find the actual net-
stress tensor 2• which is simply a transformation from CAUCHY's tensor:
0 ..
1]
w.1r &rJ. 0 ..
]1
<= > o ..
1]
~
0 ..
]1
(29a,b)
0 ..
1] <w.1po.]q + o.1qw.JP )o pq /2 -
cpijpqopq (29c)
from ( 2 9a, b) . We see, the fourth-order tensors ~ and <P de f ined in ( 29c, d)
are only symmetric with respect to two indices:
<p ••
lJpq
(I(; .
1p Jq
o. + o1q
. I(; JP
. + I(; . o.
1q JP
+ o1p
. wJq
. )I 4 +
+ <w . o.
1p Jq
+ o1q
. I(; JP
. - I(; . o . - o1p
1q JP
. wJq
. ) 14 ( 3 I)
yl ß , o/y . (32)
t .. (&ik 1(;(~1)
kj
(33)
lJ
(-1)
t .. c.. 0 => 0 .. c.. t (34)
lJ ljpq pq lJ lJpq pq
where
(-I)
clJpq
.. (I(;. I(; .
1p Jq
+ I(; . I(; . ) I 2 wi th C ..
1q JP lJpq
( 35)
clJpq
.. cJlpq
.. clJqp
.. c ..
pqlJ
(36)
(37a)
.
d 1ag la2 , ß2 , l, aß/2, ßr/2, ra/2l , (37b)
2
t31 = t 13 ' t32 = t23 ' t33 = 03/1
where dP. is the actual force vector (28). According to ( 15) we can
1
write:
(-I)
dP. o . . \(!. dS
A
(40a)
l J1 Jr r
dP.
1
w.1p t pr ds r (40b)
. h can b e mu l t1p
wh 1c . l.1ed b y '~ki I) 1' ( -
so that we have Wk(~l)dP.
l l
tk r dS r , or
after changing the indices:
dP.
1
t ..
Jl
ds.J (4 I)
240 ]. Betten
One must consider the fact that the symmetric tensor-valued functions
(Sa,b) arevalid for an isotropic material in an anisotropic damaged
state. Furthermore, one must differentiate between the anisotropic da-
mage growth and the initial anisotropy resulting from a forming process,
for instance rolling. Then constitutive equations can be found in ex-
pressions such as those in (I).
In formulating constitutive equations like (I) one has to take the fol-
lowing into account: the undamaged case (w ~ O) immediately leads to
(Sa), while the strain rate tensor ~ approaches infinity as w approaches
6. With respect to polynomial representations of constitutive equations
+)
it is convenient to use the tensor
(-I) (-I)
D .. _ ( 6 ..
1J 1J
w .. )
1J
- wlJ
.. (42)
d .. f .. (o , D , A ) (43)
1J 1J pq pq pqrs
must be taken into considerat ion, i. e. , where f .. = f.. are the compo-
1J Jl
nents of a second-order tensor-valued function of two symmetric tensors
of rank two and one symmetric tensor of rank four. The general represen-
+)
Then, the pseudo-force vector (41) 1s defined as dP. - D.. P., for 1n-
stance. 1 1J J
Formulation of Anisotropie Constitutive Equations 241
tation is given through a linear cornbination (2), where the g's are syrn-
rnetric tensor generatorsofrank two involving the argurnent tensors ~.
or [ 8] :
a.ka. 1 fk 1 (o , ... ,A
1 J pq pqrs
) = f .. (a a o , ... ,a a a a A
lJ pt qu tu pt qu rv sw tuvw
) (46)
l
tr d 0 tr L: ~a Qa 0 - L: ~ a tr Qa 0
a a
(4 7)
tr d D tr ~ A•D tr d A•D - d .. A.. D
lJ lJpq pq
242 J. Betten
3 f . .
.
Alternat1vely, we can form tr Q, tr Q2 , and tr d rom the l1near combl-
nation (2). The sets of invariants so formed contain reducible and equi-
valent elements, i.e. redundant elements. In writing down the possible
irreducible invariants, we can use some results from the theory of in-
variants, for instance the following LEMMAs [ 17]:
- If Q is a reducible generator, then tr XG and tr GX are redu-
cible, where ~ is an arbitrary tensor.
- The trace of a matrix product is unaltered by cyclic permutations
of its factors.
- The trace of a matrix product is equal to the trace of the trans-
pose of the product.
- The transpose of a matrix product TI 1s the matrix product formed
by writing down the transposes of the factors of TI in reverse or-
der, for instance, tr ~Qf = tr f Q ~
1 1 1 •
( 48)
I [ A, v] (M~~,)J,V] + i~,)J,V))/2,
G ..
lJ lJ Jl
(49a,b)
Formulation of Anisotropie Constitutive Equations 243
3[A)JV] - (N~~,JJ,V]
G .. ' ' \,v 0,1,2; )J 0, I, ... ,5, (49c)
lJ lJ
(2)
d .. cpoo·. + cplo .. + cp2o .. (52)
lJ lJ lJ lJ
as a special case.
The representation (51) may be extented by using matrix products of
the types:
244 J. Betten
(2) (2)
f .. o .. , f. .D .. , f. .A.. o , f. .A.. o , f. .A.. D , f. .A.. D ,
lJ Jl lJ Jl lJ Jlpq pq lJ Jlpq pq lJ Jlpq pq lJ Jlpq pq
(55)
T •• T •• (56)
lJ Jl
d .. (57)
lJ
(58)
lJ.(&pq '
d .. f.
lJ
1
pq )
d .. f. .(t , A ) (59)
lJ lJ pq pqrs
d .. (60)
lJ
(A,v = 0, 1,2; ~ = 0, 1, ... ,5), where the ~·s are scalar-valued functions
of the integrity basis.
Finally, instead of (58), we should consider the function
d .. f .. (t , T ) ( 6 I)
lJ lJ pq pq
l
2 3
tr _:::, tr T tr T
( 64)
2 2 2 2
tr ~_:::, tr tT , tr ~ _:::, tr t T
[/.]
are valid for arbitrary values A and v, where the symbols ~ijkl and
[v]
nijkl are A and \! several fourth-order tensors defined in the following
way:
(2)2 2 (2)]
+ [ t .k xk.(T) + x.k(T)t,. 12 , ( 6 7a)
1 J ~ 1 ~ l<J
I 2 [ (A) ). 1 0.)]
d .. -2 L: t.k Xk.(T) + X.k(T)tk. . (67b)
1J \=0 1 J ~ 1 ~ J
2
>x 1J.. ( T) L: cp(\,v) T ••
( v)
1J
(68)
~
v=O
Inserting (68) in (67b), the usual polynomial (62) is obtained.
In a similar way one can represent the symmetric tensor-valued
function (57) of three symmetric argument tensors:
I 2 [ (A) \ \ 0.) J
d .. L:
2 \=0 O.k Xk.(D,T) + X.k(D,T)Ok. +
1J 1 J ~ ~ 1 ~ ~ J
+ l ~ [D(\) Ay (T o) + Ay (T o)D(;\)J +
1'k k'J ~ ' ~ (69)
{_,
2 \=0 1'k ~ ' ~ k'J
I 2 [ 0.) \ (A) J
\ z.k(o,D)Tk.
+ -2 L: T.k zk.(o,D) + ,
\=0 1 J ~ ~ 1 ~ ~ J
where nine isotropic tensor functions each containing two variables are
introduced:
\
X .• (D,T)
2 [ (~) (v) (v) (~)] 0, I, 2. (70)
1J 2 L: Cfl( ' ) D. k Tk . + T. k Dk . ; A
J J
v= 0 "'~'V
~ ~
1 1
~.
REFERENCES
13. SPENCER, A.J.M. and Rivlin, R.S., Finite Integrity Bases for Five or
Fewer Symmetrie 3 x 3 Matrices, Arch. Rational Mech. Anal. 2 ( 1958/59):
435-446.
14. SPENCER, A.J.M., Theory of Invariants, in: Continuum Physics. Vol. I,
(ed. A.C. Eringen), Academic Press, New York 1971, 239-353.
15. BETTEN, J., Integrity Basis for a Second-Order and a Fourth-Order Ten-
sor, International J. Math. & Math. Sei. 5 ( 1982): 87-96.
16. BETTEN, J., Irreducihle Invariants of Fourth-Order Tensors, Fifth
International Conference on Mathematical Modelling, University of Cali-
fornia, Berkeley 1985, to be published in the proceedings: Mathematical
Modelling in Science and Technology (eds. X.J.R. Avula, G. Leitmann,
C.D. Mate, Jr. and E.Y. Rodin), Pergarnon Press, New York 1986.
17. RIVLIN, R.S. and SMITH, G.F., Orthogonal Integrity Basis for N Symme-
trie Matrices, in: Gontributions to Mechanic (ed. D. Abir), Pergarnon
Press, Oxford/ ... /Braunschweig 1969, 121-141.
18. SPENCER, A.J.M. and RIVLIN, R.S., The Theory of Matrix Polynomials and
its Application to the Mechanics of Isotropie Continua, Arch. Rational
Mech. Anal. 2 ( 1958/59): 309-336.
19. SMITH, G.F., On Isotropie Function~ of Symmetrie Tensors, Skew-Symme-
tric Tensors, and Vectors, Int. J. Engng. Sc. 9 ( 1971): 899-916.
20. BETTEN, J., On the Representation of the Plastic Potential of Aniso-
tropie Solids, Colloque International du CNRS n° 319, Comportement
plastique des solides anisotropes, Grenoble 1981, published in the
the proceedings: Plastic Behavior of Anisotropie Solids (ed. J.P.
Boehler), CNRS, Paris 1985, 213-228.
J. Betten
Technical University Aachen, F.R. Germany
I . INTRODUCTION
n
y f(x) = L La(x)ya +Rn ( I)
a =1
n
L (x)
a TI ( x - x 6 ) I ( xa - xß ) , a = I ,2, ... ,n (2)
ß=I
s+a
are introduced. lt 1s clear that La(x 8 ) is equal to one for a ß and
equal to zero for a ~ ß. The remainder in (I) is given by
252 J. Betten
R (3)
n-1
where min (x, x 1, ... , xn) < ~ < max (x, x 1, ... , xn).
Since the above interpolation has a unique solution, all other re-
presentations must coincide with the LAGRANGE polynomial (I). But the
LAGRANGE formula is disadvantageous, if we pass from a space of dimen-
sion n to a space of one higher dimension. Then the "NEWTON formula"
(4)
Ill
Y •.
lJ
f. . (X)
lJ -
I Cl.
f"' L .. (XJ + R .. (X)
.... lJ - lJ -
(6)
Cl.= I
R .. (Ba)
~J
(9a)
(9b)
(9c)
( 10)
<). .. a o
<po ij +
am X
'~'I ij +
am x\~)
'~'2 lJ ( I I)
~J
( 12b)
254 ]. Betten
( 12c)
with
( 13)
(2)
y.. = .. + <ll IX.lJ. + <ll2Xl. J.
<ll 0° lJ ' ( 14)
lJ
III
<ll 'J =L 'J 0, I ,2 . ( I 5)
O:=l
The coefficients ( 15) can also be found by solving the system of linear
equations
( 16)
a. a. Y f. .(a. a. X ) ( 17)
lp Jq pq lJ lp Jq pq
y .. f. . (X)
lJ lJ ~
( I 8)
Further terms are not possible, because the second rank tensor ~ has on-
ly three principal values. Moreover, the last term in ( 18) vanishes be-
cause of the HAMILTON-CAYLEY theorem. Thus we have the formula
y .. ( I 9)
lJ
(20a,b)
( 20c)
(2)
y ..
lJ
= lj!Oo .. + \)!IX .. + I)J2X. .
lJ lJ lJ
, ( 2 I)
a - (22a)
0
(22b,c)
coincide with (14), i.e., the scalar coefficients (15) and (22a,b,c) are
identical:
v=0,1,2, (23)
z = f(x,y) (24)
m n
z = I I L (x)L (y)f(x ,y ) + R(x,y)
a 0 a 0
(26)
a= 1 o= 1
(x - x 1) (x - X ) amt<~,,y)
m
R (28a)
m-1 m! axm
n
(y - y I) (y - y ) a f(x,n 1)
n (28b)
R
n-1 n! ayn
Interpolation Methods for Tensor Functions 257
M(xy)
(28c)
Z .. =
lJ
[f. .(X,Y) + f .. (X,Y)]/2 =
lJ~- Jl~~
Z ..
Jl
(29)
III
f .. (X,Y) =
lJ ~ ~
I aL .k(x) 0 Lk.(_Y)f(a
a ,p=I 1 - J - ,p
) + R .. (X,Y)
lJ - -
, (30)
(X.-,\(
1p a+ I) 6 1p
.)(Xk-,\(
p a+ II) 6 pk) ( 3 I a)
( yk g - )J ( p + I ) 6 k q ) ( y qJ. - )J ( p+ I I ) 6 <.LL
.)
PLk. (Y) - (31b)
J ~
(JJ(p) - JJ(p+I))(JJ(p) - u(p+II))
where ,\(a)' a = I,II,III, and )J(p)' p = I,II,III, are the principal va-
lues of the argument tensors X and I• respectively. The remainder in
(30) is given by
f. .(X,Y)
lJ - -
(33)
I I
A .. - (X .. AI 6 .. ) B .. - (Y .. - lli 6 .. ) ( 34a, b)
lJ lJ lJ ' lJ lJ lJ '
2
A .. - (Xik AI 6ik)(Xkj AII 6kj) (3Sa)
lJ
2
B .. = (\k- lli 6ik)(Ykj- llii 6kj) ( 3Sb)
lJ
are used. The nine coefficients a 0 , a 1, ... , c 4 in formula (33) are de-
termined by the nine interpolation conditions
ao - f(I,I)
III III
a2 - L f(a,I/M(a) ' b2 - L f(I,P/N(P)
a=I P=I
Interpolation Methods for Tensor Functions 259
III
c2 - L (f(I,p) - f(II,p))/[(AI- AII)N(p)J
P=I
III
c3 - L ( f(a, I) - f(a,II))/[(~I- ~II)M(a)J
Ct=l
III
L f(a,p)(A(a+I) - A(a+II))(~(p+I) - ~(p+II))/
a,p=I
are used.
The tensor-valued function (30) or, alternatively, the approach (33) can
be represented by the same minimum polynomial
2
( ) 1 ,1, l~t)Y(v) ( 37)
f. . X' y = L 'I' (1t ,v) ik kJ' '
lJ~~ lt,V=O
t 1•
lJ
• - . 1ax PJ.
ay 1p \)! I o.. +
lJ
21);2X l. J. ' (39)
el J 1• _ at~
lq
1axqj. = 2w 2 o..
lJ
(40)
For example, 1n the case of AI = AIII + All' we find from ( 19) and (39):
(4la,b)
( 4 lc)
ao YI ' aI fl (42a,b)
I '
a2 [fI -
(YI - YIII)I(AI- AIII)]I(AI- AIII) (42c)
I
fl a2 = f~ I 12 , (44a,b,c)
I
But in this special case the argument tensor is a spherical one, X.. =
lJ
A1 6 .. , and therefore the formula ( 19) reduces to the trivial result:
lJ
262 ]. Betten
f .. = fiO ... Note, that the interpolation formulas for scalar functions
lJ lJ
y f(x) approaches the TAYLOR expansion for f(x) at x 0 , if we make xa'
a = 1,2, ... ,n, coincide at x 0 .
In extending the interpolation at coincident points to tensor
functions of two variables we need derivatives of the function (29) with
respect to both the variable X and Y. For example, we have to use the
derivatives
JP
.
z - caz. ;ay + az.JP ;aY pl.)/2 (45b)
lJ lp pj
0 0
az ~ az'. az.
. .
az.
z' -
]
4
(~ + ________lE + ~ +
aY . aY
pi
ax
pj
~)
pl
( 45c)
PJ
and so on.
zlJ 0 0 f. .(X
lJ pq
,Y
pq
,A
pqrs
) (46)
z (4 7)
where the G's are symmetric tensor generators of rank two involving the
argument tensors ~. ~. 6· Such generators, for instance, are given by
matrix products of types
Interpolation Methods for Tensor Functions 263
6. SIMPLE EXAMPLES
X •. (48)
lJ
y .. X~ ~/3). (49)
lJ lJ
we calculate:
a = 41/3
0
(50b)
( 50c)
21/3 + 21/3 -
0
22/3 i13
X~ ~/3) 0 21/3 0 ( 5 I)
lJ
21/3- 21/3 +
0
22/3 22/3
y ln X (52)
(53)
ln ~ ln 2 (3~2 0
0 ) (54)
0 3/2
(55b)
or d Ka n , (56a,b)
266 ]. Betten
where oij = oij - okk6ij/3 are the components of the stress deviator o'.
For the special case of incompressible behaviour (dkk = 0) we find from
(58):
(59)
d .. cp lo!. (60)
lJ lJ + cp2o!lJ'.
. {o, O, O) , ( oij
( oij ) V = d 1ag , ) V= d1ag
. {2 I
)Ü• - )Ü• I ) ,
- )Ü (62a,b)
d ( 63)
( 64a)
n-1
a1 = (I + v)Ko . ( 64b)
The derivative fii at the coincident points XII XIII can be determined
in the following way. From (63) we derive:
or
(65b)
n-2 ( 64c)
a 2 = ( I + v - 3vn)Ko
m
't'Q
= l9 ( I - Sv + 6vn)Kon (66a)
32 ( I + v + 23 vn)Kon-1 (66b)
n-2 · ( 66c)
~2 = ( I + v - 3vn)Ko .
d ( 67)
(68)
( I - 8v + 6vn)Kon - 2( I + v - 3vn)Kon ,
( 70a)
For uni-axial loading (62b) we find with (66a,c) the trivial result:
D
. o . . d .. (71)
lJ Jl
(72)
if the constitutive equation (58) is inserted into (7 I). For the uni-
axial case (62a,b,c) with J 1 : : o, J' : : o2 /3, J) :::: 2o2/27 we find from
2
(72) and (66a,b,c)
Kon+ I =
_ d
o (73)
270 ]. Betten
o . . d .. _ od (74)
lJ J 1
.
D (75)
d .. a 1 o~ . ( 76a)
lJ lJ
3 n-1 ,
d .. -2 Ko o .. (76b)
lJ lJ
n-1
d .. l2 K ( 3J ' ) 2 o ~ .
lJ
(76c)
lJ 2
2)
0.
1 (
0 <. 2.) - 2J 1o.. /3 + J 21u-" .. /9 ( 77)
lJ
0
lJ lJ lJ
into equation (58), one can find a constitutive equation of the form
(57) with the scalar coefficients
2
cp5 cpo - J 1cp/ 3 + J lcpzl 9 ( 78a)
0 (79)
0 (82)
o.lJ. o(·l·o*o
'~-'
.. + •I·*E:
lJ
.. + ·I·*P~
'~-'I lJ '~-'2 lJ
) ,
Wo ( I - 2121)!2)/3 ,
2 2
1/! 2 = [ 1/( I + v) - C]/( I + v) E:
C 19+4v 1 - 8 ( 1 + v) J
-z(l+v)2 (9 + 4v) 2
are used.
Another example discussed in [ 17] is the tensorial generalization
of the RAMBERG-OSGOOD stress-strain relation [ 18]
E = o/E + k(o/E)n .
1 ( 2)
(..
~J
= cpocS ~J
.. +cp1o~. + cp2o. .
~J ~J
'
cpo 3
1 1-n)-
-( 0
+
E
t< 1-8V+6Vn)[o/E+k(u/E)n] '
E.. = cp*;:
0 u.. + cp*o.. + cp*2o (2)
..
~J ~J 1 ~J ~J
Now, the ma1n problern is: to determine the nine scalar coefficients
I
d .. = f. .(o,D) 2(Xikykj + Yikxkj) ( 85)
1J 1J - -
l
X.. X.. (o) = cp86 .. + cp'f 0 . . + cp2*0 (. .2)
1J 1J ~
1J 1J 1J
(86)
,\
ct<Ctro ) ~( 0 I' 0 11' 0 III)
~ V ~
l
(2)
y ..
1J
Y.. (D) =
1J ~
<Po o1J
.. + <P ID .. +
1J
<P2D ..
1J
( 87)
lJ V 0, I, 2, (88)
'
where the scalars cp'~ are determined in (78a,b,c). Similarly to ( 15) the
V
coefficients ~ can be found by solving the system of linear equations:
lJ
m
2 (D ) I
~0 + ~IDI + ~2DI I
2 m
~0 + ~IDII + ~2DII ( D ) II
II (89)
m
~ 2 (D ) III
~0 + ~IDIII + 2DIII III
The exponents mi, ... ,miii 1n (89) are determined by using the creep law
(83) in tests on specimens cut along the mutually perpendicular direc-
tions x 1 , x 2 , x 3 .
(-I)
Because of D. . :=
lJ
w.. )
lJ
( ol..J -
and
w.. = d iag {( I - a) , (I - ß), (I - y)}, the principal values in (89) can
lJ
be expressed through [ 14]:
where a,ß,y are fractions which represent the net cross-section elements
of CAUCHY 1 s tetrahedron perpendicular to the coordinate axes [3].
In the case of two equal parameters, for instance a 1 ß = y, the scalars
~ can be determined by using the interpolation method described 1n
lJ
section 4.
Instead of (86) we can use the isotropic tensor function
.t I 1 ( 2)
X .. X .. (o 1 ) cpou .. + cplo .. + cp2o .. ( 9 I)
lJ lJ - lJ lJ lJ
I D()J) l(v)J
d .. ( 92)
lJ 2 + ik 0 kj '
and experimental data. To show this we can start from (74) and find sl-
milarly to (79) the cubic equation
0 ' (94)
(9Sa)
(95b)
(95c)
D- (95d,e)
We see, that the elements of the integrity basis (93) and experimental
Interpolation Methods for Tensor Functions 277
(96)
In the case of darnage (w .. ) and initial anisotropy (A. "kl) the con-
lJ lJ
stitutive equations discussed in [3,5, 14] are tensor-valued functions
like (46) which can be represented as (47), where the scalar functions ~a
REFERENCES
J. Betten
Technical University Aachen, F.R. Germany
I . INTRODUCTION
(2)
0 .. cp06 .. + cpiE .. + cp2E .. ' ( I)
1] 1J 1J 1J
(2a,b,c)
Eliminating the elastic potential, one can find the following sufficient
and necessary conditions
di\(ClF/Clo .. ) (4)
~J
components of the plastic part of the strain tensor E, and the factor dA
~s LAGRANGE's multiplier.
In the following the classical plastic potential theory is compared
with the tensor function theory. It will be shown, that the former theo-
ry is compatible with the latter, if the material is isotropic, and if
additional conditions are fulfiled.
However, for anisotropic materials the plastic potential theory on-
ly furnishes restricted forms of constitutive equations, even if a gene-
ral plastic potential has been assumed [2]. Consequently, the classical
flow rule (4) must be modified for anisotropic solids. In the following,
appropiate modifications are discussed and resulting conditions of "com-
patibility" are derived.
2. ISOTROPY
variants SV = tr ov,
~
V= 1,2,3, of the stress tensor. Thus, considering
the form F = F(s 1,s 2 ,s 3 ) and using the classical flow rule (4), one im-
Tensor Function Theory and Classical Plastic Potential 281
( 2)
t. . ( o ) = cp 0 cS. . + cp 1o. . + cp 2 o. . , (6)
1J pq lJ lJ 1J
( 8)
Thus, in the isotropic special case the plastic potential theory is com-
patible with the tensor function theory, if the conditions of "compati-
bility" in (8) have been fulfiled.
For example, the requirements 1n (8) can be satisfied, if one assu-
mes:
( 9)
F ( I 0)
( I I)
where
( 12)
( 13)
( 14)
cpl 0
0 - 3cp 0 , ( 15)
0
0 - ]cp
I
one can assume, that the right hand side of the system ( 14) is also
equal to zero. It is possible to find the three conditions of "compa-
tibility" in (8). Note, that these conditions are only sufficient. But
if we use CRAMER's rule, we find the following sufficient and necessary
condition
( 16)
which contains (8) as a special case. As has already been mentioned be-
fore, the conditions in (8) are only sufficient, whereas the similar
conditions (3) in elasticity are both sufficient and necessary.
For more complicated examples than the above mentioned the plastic
potential theory is not compatible with the tensor function theory. For
instance, MURAKAMI and SAWCZUK found in [3], that for the model of
prestrained plastic solids the flow theory of classical plasticity only
furnishes restricted forms of constitutive equations. In [4] LEHMANN
proposed an extended form of the classical flow rule (4) with respect to
experimental results and based upon thermodynamical considerations. In
[5] oriented solids are considered and an appropiate modification of the
classical flow rule (4) is discussed.
3. GRIENTED SOLIDS
2
dpE:. . f. .(o )
I
E
( [A,v] M~~·vJ)
lJ lJ pq' pq
A
2 cp(A , V ) MlJ
.. +
Jl
( 17)
A, v=O
where [A,v] several symmetric tensor generators are formed by matrix
products of the forms
The cp(A;V) 's in ( 17) are nine scalar functions of the integrity basis,
the elements of which can be found by forming all the irreducible traces
of the matrix products ( 18):
SA tr aA T tr Av ; A,V 1,2,3
-
V-
}
~
( 19b)
tr aß Q2 tr OA 2 - tr Ao2 tr o2A2
Ql ~~ ' Q3
-
- ~
Q4 -
From the ten invariants in ( 19) only the seven stress-dependent lnva-
riants are essential for the plastic potential:
Thus, by using the flow rule (4), one can find the constitutive equation
( 21)
:\-I
[ :\]
Qpqij -
ao (:\ 7ao.. ~
{_,
[ 0 (a) (:\-I-a)
. 0 . + 0
( a)
. 0
CA-
.
I -a) J 12
(22)
pq lJ a=O pl qJ PJ ql
One can see, that because of (22) the value A cannot be equal to zero.
Therefore, only the seven stress-dependent invariants of the integrity
basis ( I9) are relevant for the plastic potential (20). Furthermore, by
comparing the constitutive equations ( I7) and (2I), one can see, that
the minimum polynomial ( I7) consists of nine tensor generators, whereas
(2I) has only seven terms, i.e., the plastic potential theory furnishes
only restricted forms of constitutive equations, even if a general
plastic potential is assumed in an anisotropic case.
It may be more useful for practical applications to represent the
constitutive equation ( I7) in a canonical form
(24)
0 (0) [2]
Hijkl - Cflco,o)mijkl + cp( 0, I )mijkl + cp(0,2)mijkl (2Sa)
I (0) [ 2J
Hijkl -
cp( I , 0) mi j k l + cp( I, I)mijkl + cp( I ,2)mijkl (2Sb)
2 (0) [2]
Hijkl -
cp(2, O)mijkl + cp(2, I )mijkl + cp(2,2)mijkl (2Sc)
[v]
having the symmetric tensors m \! = 0, I,2, of rank four:
[v]
mijkl (26)
Especially, for \! = 0 one finds in (26) the zero power tensor of rank
four:
286 J. Betten
[ A] 0) 0) 0) (A)
~ ijkl -
cxik 0 jr + xi 1 6 jk + 6 ikxj 1 + oilxjk )/4
[v] ( v) ( v) 6 Y( v) ( v)
n ijkl -
(Yik 0 jr + y i l 0 jk + 1k jl + 0 iryjk )/ 4
have been used [5], where X and Y are two symmetric second-order ten-
sors. Similar identities for non-symmetric tensors were introduced in
[ sJ.
Similarly to (24) one can represent the constitut!ve equation (21),
derived from the plastic potential theory, 1n a canonical form:
dp€ .. (28)
1J
0
hijkl - (29a)
[6]. Therefore, only v = 0, I in ( 18) has tobe taken into account. The
number of tensor generators in the constitutive equation ( 17) is reduced
from nine to six, i.e., the fourth-order tensor-valued functions (25a,b,
c) can be simplified into:
0 (0) (30a)
Hijkl - <p(O,O)mijkl + <p(O, l)mijkl
I ( 0) (30b)
Hijkl - <p( I ,O)mijkl + <p( I, l)mijkl
2 ( 0) (30c)
Hijkl - <p(2,0)mijkl + <p(2,1)mijkl
0 [ ( 0) (32a)
hijkl - <aF/as 1)mijkl + (aF/a~ 1 )mijkl]dA
I [ ( 0) (32b)
hijkl - 2 (aF;as 2 )mijkl + <aF/a~ 3 )mijkl]dA
2 [ ( 0) 02c)
hijkl - 3 (aF/as 3 )mijkl + 0 ]dA .
[2]
c) one can see, that any terms containing mijkl and mijkl do not appear
1n (29c) or in (32c), i.e., the "second-order effect" in (28), characte-
rized through (29c) or (32c), is not influenced by the anisotropy (26)
of the material. Consequently, the flow rule (4) of the classical flow
theory of plasticity should be modified in the anisotropic case.
The results mentioned above can also be applied to perforated ma-
terials or damaged materials, if the anisotropy tensor A.. in (26) is
lJ
substituted with a perforation tensor [9] or a darnage tensor [8, 10].
( 34a)
( 34b)
( 34c)
a: 3F
- 3d/\ 3F
cp(2, I) - 7 d/1 "'dTIJ - a: d/\ ~~4 ' (34g,h,i)
dS3
where the abbreviations
(3Sa,b,c)
have been used. In this way the second-order effect 1n (28) is in-
fluenced through the material tensor ~ because of a: f 0 found in (34h,
i). This influence is not at all possible, if we use the tensor-valued
function (29c) which results from the classical flow rule (4) which does
not contain the parameter a:. Thus, the plastic potential theory is com-
patible with the tensor function theory, if the modified flow rule (33)
is used instead of (4). But among the scalar functions (34a, ... ,i) we
can find additional relationships ("conditions of compatibility") by
eliminating the plastic potential F. This can be done by closely adhe-
ring to the following procedure. First, one must find the partial deri-
vatives 3F/3S 1, ... ,3F/3~ 4 of the plastic potential F from the system of
linear equations (34a, ... ,i). Second, one must apply the {~0)equalities
to every pair of the obove mentioned ten derivatives. Thus, in this way
we have eliminated the plastic potential F and found 45 conditions of
compatibility.
If the vector v 1s regarded as a unit vector, then it is possible-
to find the following identities by first taking (30a,b,c), (31) into
consideration and by then using the modified flow rule (33):
The results of (37a, ... ,f) can be used to find the partial derivatives
of the plastic potential:
ap;as 1 - <p(0,0) 1 •
ap;as 2 - <p< 1,o)L/2, aF/as 3 - <p(2,0)L/3 '
(38)
aF;an 1 - <p(O, I)L' aF;an 3 - 2<p( 2 ,I)L/a. f
In (38)
(39)
a
cp(2, 1) 1 • I + 0 + 0 + 0 - Jil(0,0) 1 '5
R4 1
0:
0 + 2 cp(2, 1) 1 '2 + 0 + 0
- Jil( I ,0) 1 '5
R71
a
0 + 0 + 3 cp ( 2, I) 1 ' 3 + 0 - Jil(2 ,0) 1 '5
R91
a
0 + 0 + 0 + cp ( 2 ' I ) T~ ' 4 - 2'-Pco, I) 1 •s Rl01
( 4 I)
(42)
and
292 J. Betten
(43)
are introduced.
From the equations found in (41) one takes (~) systems each of
which consists of 3 linear equations. In other words one starts from
(38) and combines (~) identities. Then, by using the (i} equalities
(40), one can find 10 systems each with 3 linear equations, for in-
stance:
(44)
a
=> 2cp ( 2 , I ) 1 ' 2 + 0 -zcr(I,0) 1 •5 R71
F,4 - cp(O,I) 1
a
F,5 - 2cp( 2 , !)1/a 0 + cp(2, 1) 1 '4 -zcr<o. 1) 1 '5 Rl01.
(45)
Note, that the coefficient determinants of the systems (44), (45), and
all other systems, which have been found in the same way, are equal to
zero:
0 2cp ( 0, I ) -cp(I,O) 0
cp(2 ,0) - 3cp(O,O)
a a
0 2cp(2,1) 0 - ~(I ,0) 0.
cp(2, I) - Jil(O,O)
a a
0 3cp(2, I) - Ji>(2 ,0) 0 cp(2, I) -~(0,1)
(46)
(4 7)
F (49)
where g 1, ... ,g 5 are arbitrary functions of only one invariant. The rea-
der may note the analogy between (49) and ( 10).
Note, that the conditions (47) are only sufficient, but not neces-
sary. In order to find sufficient and necessary "conditions of compati-
bility" CRAMER's rule should be used. After applying this rule on (44)
and (45) and after taking (46) into account, one can find the following
sufficient and necessary "conditions of compatibility":
0 (50)
and
( 5 I)
In the isotropic special case the ten conditions (52) are reduced to the
first one of the ten equations in (52), which is identical to ( 16).
Tensor Function Theory and Classical Plastic Potential 295
The results (37a, ... ,f) and (52) point out, that for oriented so-
lids the plastic potential theory can be made compatible with the repre-
sentation theory of tensor functions, if the classical flow rule (4) is
modified according to (33) and if, in addition, the "conditions of com-
patibility" in (52) have been fulfiled.
Finally, the special case of incompressibility should be consi-
dered. In such instances it is more practical to use the stress deviator
~· in the integrity basis ( 19a,b) than to use the stress tensor Q·
Furthermore, the modified flow rule
(53)
(54)
(55)
f. .(o ,A ) (56)
~J pq pqrs
= G[;\,w,vl (58)
'A,w,v cro,w,v) ij
(59)
( 60)
Table I Matrix Products of the Fonns (60) for Isotropy (1-l 0) and
Anisotropy <1-1 ~ 0)
~ 0 z ~ 0 z
~~------~--------+---------~
d.m A_t.l-ll 1!-ll A.l.l-ll diZI
0 d..
IJ IJ
0 IJ r r Aijpq dPQ IJPQ PQ
II-II II-II 11-ll 12)
d.k
1 Ak )rr
. dik Akjpq dPQ dik Akjpq dPQ
diZI A l1-1l
ik kjrr
diZI A11-ll d dm A 1 ~ 1 d 121
1k kJPQ PQ 1k kJPQ PQ
0 z
(IJ.)
0 A..
IIJJ
z symmetric
One can see, that the combinations [A,~, v ] are symmetric with
respect to A,V for all ~. i.e. in eqn (61) or in Table 2 we see 33 In-
.
var1ants . I n ad d.1t1on,
. t h e 1nvar1ant
. . A(. 6. ).. must be cons1dered.
· Thus a
llJJ
set of 34 irreducible invariants is found.
298 J. Betten
F ( 62)
and by using the classical flow rule (4) one obtains the constitutive
equation
(64)
f. .(o ,D ,A ) (66)
lJ pq pq pqrs
REFERENCES