Methodes Numerique

Download as pdf or txt
Download as pdf or txt
You are on page 1of 32

F INITE D EFORMATIONS

David Ryckelynck, Jacques Besson

Centre des Matériaux, Mines ParisTech, UMR CNRS 7633


BP 87 Evry cedex 91003, France
Notations

scalar a a
vector ~a ai
2nd order tensor a aij
4th order tensor a aijkl
matrices a
Voigt notation a→a
V(a.b)

Products
. c = ~a.~b c = ai bi
~c = a.~b ci = aij bj
c = a.b cij = aik bkj
: c=a:b c = aij bij
c=a:b cij = aijkl bkl
⊗ c = ~a ⊗ ~b cij = ai bj
c =a⊗b cijkl = aij bkl
[Einstein convention]

MAGIS Finite Strains 2/1


Outline

Displacements and deformation


Stress measures
Constitutive equations

MAGIS Finite Strains 3/1


Displacement and Deformation

MAGIS Finite Strains 4/1


~
Displacement field around X

~
~ + dX
~u (X ~ ) = ~u (X ~ ) + ∂ u .d X
~ ) + d ~u = ~u (X ~
~
∂X
∂~u
„ «
~ + d ~u =
d ~x = d X 1+ ~
.d X
~
∂X
Transformation gradient
∂~x ∂xi
F = FiI =
~
∂X ∂XI

MAGIS Finite Strains 5/1


Objectivity

Rigid body transformation


~x 0 = Q(t).~x + ~c (t)
Quantities are objective if they are related by the rotation tensor as:

m0 = m
~u 0 = Q(t).~u
0
T = Q(t).T .Q(t)T

Generalization
0
T(n) = ~u1 ⊗ · · · ⊗ ~un objective if T(n) = ~u10 ⊗ · · · ⊗ ~un0 where ~ui0 = Q.~u1

F is not objective
∂~x 0 ∂~x 0 ∂~x
F0 = = . = Q.F
∂X~ ∂~x ∂ X~
B = F .F T is objective
B 0 = F 0 .F 0T = Q.F .F T .Q T = Q.B.Q T

MAGIS Finite Strains 6/1


. . . and invariance

Quantities are invariant if they remain unchanged by the transformation

m0 = m, ~u 0 = ~u , T0 = T

C = F T .F is invariant
C 0 = F 0T .F 0 = F T .Q T .Q.F = C

MAGIS Finite Strains 7/1


Rates

~ be an infinitesimal segment in the current configuration. One gets:


Let δ~x = F .δ X

dδ~x ~
dF .δ X dF ~
= = ~ = Ḟ .F −1 .δ~x = L.δ~x
.δ X = Ḟ .δ X
dt dt dt
L can be separated into symmetric (D) and an skew-symmetric (W ) parts:

L = Ḟ .F −1 = D + W

D characterizes strain rate in the following way:

d “ 1 2” dδ~x 1 dδ~x 2
δ~x .δ~x = .δ~x 2 + δ~x 1 . = (Dij + Wij )δxj1 δxi2 + δxi1 (Dij + Wij )δxj2
dt dt dt

=
h i h
0
Dij δxj1 δxi2 + Dij δxi1 δxj2 + Wij δxj1 δxi2 =+ Wij δxi1 δxj2
i

= 2δ~x 1 .D.δ~x 2

MAGIS Finite Strains 8/1


For a transformation such that: F 0 = Q(t).F

L0 = Ḟ 0 .F 0−1 = (Q̇.F + Q.Ḟ ).F −1 .Q −1


= Q̇.Q T + Q.(D + W ).Q T Q −1 = Q T
= Q.D.Q T + Q̇.Q T + Q.W .Q T
= D0 + W 0

With

D0 = Q.D.Q T
W0 = Q.W .Q T + Q̇.Q T

Note that Q̇.Q T is skew-symmetric as:

Q.Q T = 1 ⇒ Q̇.Q T + Q.Q̇ T = 0

Only D is objective

MAGIS Finite Strains 9/1


R.U — V .U decomposition
The deformation gradient F can be decomposed, using the polar decomposition theorem, into
a product of two second-order tensors
F = R.U = V .R
FiJ = RiK UKJ = Vik RkJ
with (R rotation tensor)
R.R T = 1 U = UT V = VT U = R T .V .R

MAGIS Finite Strains 10/1


Calculation of R and U

C = F T .F ≡ U.U
det C = (det F )2 > 0 and det C 6= 0
Eigen frame for C = P T .C 0 .P
0 1 0√ 1
c1 0 0 c1 0 0

C 0 = @ 0 c2 0 A → U0 = @ 0 c2 0 A

0 0 c3 0 0 c3

U = P T .U 0 .P
Consequently
R = F .U −1 and R T = U −1 .F T
so that
R.R T = F .U −1 .U −1 .F T = F .C −1 .F T = F .(F −1 .F −T ).F T = 1
idem for V

MAGIS Finite Strains 11/1


Volume variation

Jacobian of the transformation


V
J = det F = >0
V0
so that: Z Z
• dΩ = • JdΩ0
Ω Ω0

MAGIS Finite Strains 12/1


Some strain measures

Several rotation-independent symmetric deformation tensors are used in mechanics.


Right Cauchy-Green deformation tensor [Lagrangian tensor]

C = F T .F CIJ = FIkT FkJ = FkI FkJ

Left Cauchy-Green deformation tensor [Eulerian tensor]

B = F .F T T
Bij = FiK FKj = FiK FKj

Some finite strain tensors


Objective or invariant
Must be 0 for F = 1
Must correspond to the small deformation theory for a first order Taylor expansion with respect to F

1 1“ 2 ”
Green-Lagrange E= (C − 1) = U −1
2 2
Biot strain tensor E Biot = U − 1
Logarithmic strain tensor E log = log U

MAGIS Finite Strains 13/1


Principal stretches: λi

U and V have the same eigenvalues λi and can be expressed as:

3
X 3
X
U= ~i ⊗ N
λi N ~i V = λi ~ni ⊗ ~ni
i=1 i=1

so that:
3
X 3
X
C= ~i ⊗ N
λ2i N ~i B= λ2i ~ni ⊗ ~ni
i=1 i=1

One has:
3
X
V = R.U.R T = ~ i ) ⊗ (R.N
λi (R.N ~i)
i=1

Note also that:


∂λi 1 ~ ~i
= Ni ⊗ N
∂C 2λi

MAGIS Finite Strains 14/1


Some strain measures: examples

Tension loading
ux = (∆L/L)x, uy = (∆l/l)y , uz = (∆l/l)z
0 1
1 + ∆L
L
0 0
F =@ 0 1 + ∆ll 0 A
B C
∆l
0 0 1+ l

∆L 2 ∆l 2
„ « „ «
C11 = 1+ C22 = C33 = 1 +
L l
! !
∆L 2
„ « „ «2
1 ∆L 1 ∆l
E11 = 1+ −1 ≈ E22 = E33 = 1+ −1
2 L L 2 l
Biot ∆L Biot Biot ∆L
E11 = , E22 = E33 =
L L

MAGIS Finite Strains 15/1


Simple shear

1 + γ2
0 1 0 1 0 1
1 γ 0 1 γ 0 γ 0
F = @0 1 0A C = @γ γ2 +1 0A B=@ γ 1 0A = R.C.R T
0 0 1 0 0 1 0 0 1

det F = det C = 1
Eigenvalues of C

1 2 1 1 2 1
q q
1 1+ γ + 4γ 2 + γ 4 1+ γ − 4γ 2 + γ 4
2 2 2 2
R has the form 0 1
cos θ − sin θ 0
R = @ sin θ cos θ 0A
0 0 1

Solving B = R.C.R T , yields


θ = − arctan γ/2
Finally
0 2 γ 1
p p 0
B 4 + γ2 4 + γ2 C
U = R T .F = B γ2
B C
Bp γ 2 C
p +p 0C
4 + γ2 4 + γ2 4 + γ2
@ A
0 0 1

MAGIS Finite Strains 16/1


First order Taylor expansion 0 1
1 γ 0
C = @γ 1 0A = U
0 0 1
γ
θ=−
2

MAGIS Finite Strains 17/1


Stress measures and work equivalence

Cauchy stress σ (σij )


Stress measure in the final (t) configuration
Velocity gradient
L = Ḟ .F −1 = D + W
Work Z
w= σ : D dΩ

Kirchhoff stress
Z Z Z
w= σ : D dΩ = Jσ : D dΩ0 = τ : D dΩ0
Ω Ω0 Ω0

τ = Jσ τij = Jσij

MAGIS Finite Strains 18/1


Green-Lagrange strain tensor

1 “ ”
E = F .F T − 1
2
1 “ ”
Ė = Ḟ .F T + F .Ḟ T
2
1 “ ”
or F −1 .Ė.F −T = F −1 .Ḟ + Ḟ T .F −T
2
1“ ”
F −1 .Ė.F −T L + LT = D Ė = F .D.F T
2
Second Piola-Kirchhoff stress tensor
−T
σij Dij dΩ = σij FiK ĖKL FLj−1 JdΩ0
−1 −T
= Jσij FKi FjL ĖKL Ω0
−1
= JFKi σij FjL−T ĖKL Ω0
= SKL ĖKL Ω0

with
1
S = JF −1 .σ.F −T or σ= F .S.F T
J
The Second Piola-Kirchhoff stress, S, is the work conjugate of the Green-Lagrange strain
tensor, E.

MAGIS Finite Strains 19/1


First Piola-Kirchhoff/Boussinesq stress tensor

σ : DdΩ = σ : Ḟ .F −1 JdΩ0 because σ is symmetric


−1 −T
= Jσij ḞiK FKj dΩ0 = Jσij FjK ḞiK dΩ0
= ΠiK ḞiK dΩ0
= Π : Ḟ dΩ0

The First Piola-Kirchhoff stress, Π, is the work conjugate of the transformation gradient F

1
Π = Jσ.F −T or σ= Π.F T
J

Π is not symmetric 

Finally
Z Z Z Z
σ : D dΩ = τ : D dΩ0 = S : Ė dΩ0 = Π : Ḟ dΩ0
Ω Ω0 Ω0 Ω0

MAGIS Finite Strains 20/1


Interpretation of the various stress measures: tensile test

Transformation gradient — Cauchy stress [final configuation]


0 1 0 1
Fk 0 0 σ 0 0
2
F =@0 F⊥ 0 A J = Fk F⊥ σ = @0 0 0A
0 0 F⊥ 0 0 0

Kirchhoff stress
2
τ = Jσ τ = Fk F⊥ σ τ =σ for incompressible materials

Second Piola-Kirchhoff
2
F⊥ 1
S = JF −1 .σ.F −T S= σ for incompressible materials S = σ
Fk Fk2

First Piola-Kirchhoff stress


1
Π = Jσ.F −T Π= σ ≈ engineering stress
Fk

MAGIS Finite Strains 21/1


Stress rates - Jaumann stress rate

Recall of the relation


W 0 = Q.W .Q T + Q̇.Q T
so that:

Q̇ = W 0 .Q − Q.W
Q̇ T = −Q T .W 0 + W .Q T

For an objective displacement vector ~u 0 = Q.~u , one gets:

~u˙ 0 = Q̇.~u + Q.~u˙ = (W 0 .Q − Q.W ).~u + Q.~u˙ = W 0 .~u 0 − Q.W .~u + Q.~u˙

So that:
~u˙ 0 − W 0 .~u 0 = Q.(~u˙ − W .~u )
This allows to define an objective derivative of vectors (Jaumann rate):

~u J = ~u˙ − W .~u

MAGIS Finite Strains 22/1


Following the same methodology for second order tensors:

Ṫ 0 = Q̇.T .Q T + Q.T .Q̇ T + Q.Ṫ .Q T


= (W 0 .Q − Q.W ).T .Q T + Q.T .(−Q T .W 0 + W .Q T ) + Q.Ṫ .Q T
= W 0 .T 0 − Q.W .T .Q T − T 0 .W 0 + Q.T .W .Q T + Q.Ṫ .Q T

which is rewritten as:


“ ”
Ṫ 0 − W 0 .T 0 + T 0 .W 0 = Q. Ṫ − W .T + T .W .Q T

An objective derivative (Jaumann derivative) is then obtained for second order tensors:

T J = Ṫ − W .T + T .W

MAGIS Finite Strains 23/1


Stress rates - Truesdell stress rate
Recall the relation between the Cauchy and the second Piola-Kirchhoff stress:
1
σ= F .S.F T S = JF −1 .σ.F −T
J
As S is invariant the following stress rate will be objective (it corresponds to the transport of
the rate of the second Piola-Kirchhoff stress):
◦ 1
σ= F .Ṡ.F T 6= σ̇
J
Noting that:
Ṡ = J̇F −1 .σ.F −T + J Ḟ −1 .σ.F −T + JF −1 .σ̇.F −T + JF −1 .σ.Ḟ −T
So that:
◦ J̇
σ= σ + F .Ḟ −1 .σ + σ.Ḟ −T .F T + σ̇
J
Note that
d −1 −1 −1
(F .F ) = 0 = Ḟ .F + F .Ḟ
dt
so that
−1 −T T T
F .Ḟ = −L and Ḟ .F = −L
One also has
J̇/J = TrL
Finally, one obtains the Truesdell stress rate

σ = TrL σ − L.σ − σ.LT + σ̇

MAGIS Finite Strains 24/1


Stress rates - Green-Naghdi stress rate

One defines the rotated stress σ R

σ R = R.σ.R T or σ = R T .σ R .R

Following the same methodology as for the Truesdell rate, one gets:


σ = R T .σ̇ R .R

which defines an objective rate


Noting that
σ̇ R = Ṙ.σ.R T + R.σ.Ṙ T + R.σ̇.R T
one gets:

σ = σ̇ + R T .Ṙ.σ + σ.Ṙ T .R = σ̇ − Ω.σ + σ.Ω
with
Ω = Ṙ.R T
Green-Naghdi stress rate

σ = σ̇ − Ω.σ + σ.Ω

MAGIS Finite Strains 25/1


Constitutive equations : hyperelasticity

Hyperelasticity is often used for elastomers


One first defines a strain energy density function W which depends on C
For isotropic materials, W only depends on the invariants of C

I1 = TrC
1“ ”
I2 = (TrC)2 − TrC.C
2
I3 = det C for incompressible materials: I3 = 1
J = det F I3 = J 2

The second Piola-Kirchhoff stress is then given by:

∂W ∂W
S= =2
∂E ∂C

Mooney-Rivlin law
W = C1 (I1 − 3) + C2 (I2 − 3)
Ogden
N
X µp “ αp α α

W (λ1 , λ2 , λ3 ) = λ1 + λ2 p + λ3 p − 3
αp
p=1

MAGIS Finite Strains 26/1


Use of Penn invariant for nearly incompressible materials:

1
F →F = F such that det F = 1
(det F )1/3

then
C = F T .F
and

I1 = TrC
I2 = TrC.C
det C = 1

Modified strain energy density function

W = C1 (I 1 − 3) + C2 (I 2 − 3)

MAGIS Finite Strains 27/1


Constitutive equations : hypo-elasticity

The constitutive equations are written in a rate form relating any objective stress rate to the
deformation rate D:
◦ 
σ J , σ, σ, · · · = Λ : D
These constitutive equations may be path dependent . . . not physical

MAGIS Finite Strains 28/1


Constitutive equations : F e .F p decomposition

One assume an elastic (F e ) / plastic (F p ) transformation decomposition

F = F e .F p

The decomposition defines an intermediate state :

The deformation rate is given by:

L = Ḟ .F −1 = Ḟ e .F e−1 + F e .Ḟ p .F p−1 .F e−1 = Le + F e .Lp .F e−1

Express Lp = D p + W p

MAGIS Finite Strains 29/1


Crystal plasticity
X ` X
Dp = Wp =
´ ` ´
~ s ⊗ ~ns + m
γ̇s m ~ s ⊗ ~ns γ̇s m~ s ⊗ ~ns − m
~ s ⊗ ~ns
s s

~ s ⊗ ~ns ))
γ̇s = γ̇s (T : (m
Isotropic von Mises plasticity
3 T0
Dp = ṗ 0 Wp = 0
2 Teq
T rotated stress (various possibilities)

MAGIS Finite Strains 30/1


Constitutive equations : corotational formulations

The constitutive equation is expressed between the rotated stress

σ R = R.σ.R T

and any stress measure constructed using U


The small strain formalism can be used for the constitutive equation
The corresponding objective stress stress rate is the Green-Naghdi rate.

MAGIS Finite Strains 31/1


Constitutive equations : corotational formulations

The constitutive equations are expressed using:

σ Q = Q.σ.Q T

where Q is obtained so that the instantaneous rotation rate of the medium wih respect to the
frame is zero:
W 0 = Q̇.Q T + Q.W .Q T = 0
so that
Q̇ = −Q.W
The corresponding strain tensor is:
Z
εQ = Q.D.Q T dt
t

The constitutive equations then relate:

σ Q = f (εQ ) small strain formalism

The corresponding objective stress stress rate is the Jaumann rate.

MAGIS Finite Strains 32/1

You might also like

pFad - Phonifier reborn

Pfad - The Proxy pFad of © 2024 Garber Painting. All rights reserved.

Note: This service is not intended for secure transactions such as banking, social media, email, or purchasing. Use at your own risk. We assume no liability whatsoever for broken pages.


Alternative Proxies:

Alternative Proxy

pFad Proxy

pFad v3 Proxy

pFad v4 Proxy