Mech Nonlin Mats 14.5 L01 Adv Plasticity
Mech Nonlin Mats 14.5 L01 Adv Plasticity
Mech Nonlin Mats 14.5 L01 Adv Plasticity
ANSYS Mechanical
Advanced Nonlinear Materials
1
Chapter Overview
In this Lecture, the advanced rate independent nonlinear kinematic hardening
option offered in the Chaboche model is presented.
This plasticity model useful for simulating ratcheting and shakedown.
The elastic region is equal to twice the initial yield stress. This is called the
Bauschinger effect.
1
Subsequent
Yield Surface
'
y
Initial Yield
Surface
2
3
e
3
2y
3
s : s y 0
F
2
where s is the deviatoric stress, y is the uniaxial yield stress, and a is the
back stress (location of the center of the yield surface).
Note in the previous plot that the center of the yield surface has translated
a. Hence, based on the location a, the yield in reversal is still 2sy.
2
C pl
3
4
As explained on the previous slide, the term linear refers to the relationship of
back stress and equivalent plastic strain.
Can be used for cyclic loading since it includes the Bauschinger effect (elastic
region equal to twice the initial yield stress).
However, linear kinematic hardening is recommended for situations where
the strain levels are relatively small (less than 5-10 % true strain).
Because there is only one plastic slope (tangent modulus), this is not representative
of true metals as the hardening is constant.
Multilinear kinematic hardening is different than BKIN, and it uses the Besseling
(a.k.a. sublayer or overlay) model. It characterizes multilinear behavior as a series of
elasto-perfectly plastic subvolumes, each of which yields at different points, so no
back stress a is used.
2
i Ci pl g i i l
3
3
s : s R 0
F
2
Where R is a constant defining the yield stress, similar to linear kinematic
hardening.
7
There is a limiting yield surface, as explained in the next slide. In other words, the
behavior approaches perfectly plastic, unlike linear kinematic hardening, which does
not change slope.
1
Limiting Yield
Surface
C/g
Limiting
value of {a}
e
3
Current Yield
Surface
The nonlinear kinematic hardening term is associated with the translation of yield
surface. A non-zero value of g results in a limiting value of a. This means that,
unlike linear kinematic hardening, the yield surface cannot translate forever in
principal stress space. The translation is limited within a specific region.
The constant R (yield stress), which describes the size of the elastic domain, is added
on the response. If a limiting value of a exists, then a limiting yield surface will also
exist.
Nonlinear kinematic hardening is suitable for large strains and cyclic loading, as it can
simulate the Bauschinger effect. It can model ratchetting and shakedown (discussed
later).
C. Chaboche Model
The Chaboche Kinematic Hardening model is an example of nonlinear
kinematic hardening. As noted previously, the yield function is
3
s : s R 0
F
2
1
Limiting Yield
Surface
C/g
Limiting
value of {a}
3
Current Yield
Surface
10
Where:
2 n
1 dCi
pl
a a i Cie g ia i l
Ta i
3 i 1
Ci dT
i 1
Note that if n=1 and g1=0, Chaboche will reduce to Bilinear Kinematic (no limiting
value for a1).
11
R
a=a1 + a2+ a3
a2
a3
a1
R=160
C1=80000, g1=2000
C2=10000, g2=200
C3=2500, g3=0
If more the five models are necessary, use TB,CHAB and TBDATA,, in a command object
under the geometry branch for the applicable part(s).
14
15
16
2000
Ratchetting
0.06
1500
0.05
1000
Stress (MPa)
Plastic strain
Ratchetting modeled
with Chaboche using:
n=1
R=980,
C1=224000
g1 400
0.04
0.03
0.02
500
-500
0.01
-1000
0.00
-1500
0.00
10
0.01
0.02
Time
0.03
0.04
0.05
Plastic strain
0.008
t
Loading
Controlled Stress
Unsymmetry
1500
0.006
1000
Stress (MPa)
Plastic strain
0.07
2000
Shakedown
0.004
0.002
500
-500
-1000
0.000
0
Time
17
0.06
10
-1500
0.000
0.002
0.004
Plastic strain
0.006
18
On the other hand, a single nonlinear kinematic model (n=1) for the Chaboche model
can capture ratchetting, as shown below.
1500
1000
500
-1.50E-03
-1.00E-03
-5.00E-04
0
0.00E+00
5.00E-04
1.00E-03
-500
-1000
-1500
19
1.50E-03
2.00E-03
2.50E-03
3.00E-03
Unsym
Symm
Ratchetting occurs due to the fact that the initial slope in compression (A-B) is
different from the slope in tension (C-D). Since the loading is unsymmetric, C-D is
nearer to the limiting yield surface, so its slope is more asymptotic.
D
C
A
B
20
21
One kinematic model will have gi0, which will provide a ratchetting effect, as shown
previously.
On the other hand, another model will have gi=0 to provide a stabilization effect.
Recall that gi=0 is equivalent to bilinear kinematic hardening, so there is no
ratchetting.
Together, the two models will provide ratchetting with stabilization after a certain
number of cycles. This is called shakedown.
22
Stress vs Strain
23
KINH Results
24
epl
Determine k approximately
from the elastic domain. k is
usually half the elastic
domain.
Determine the plastic strain
range epl
25
e pl
C1
k
tanh g 1
2
g1
2
References
Lemaitre, J. and Chaboche, J.L., Mechanics of Solid Materials, Cambridge
University Press, 1990.
Chaboche, J.L., Constitutive Equations for Cyclic Plasticity and Cyclic
Viscoplasticity, International Journal of Plasticity, Vol. 5, pp. 247-302, 1989.
27
H. Workshop Exercise
Please refer to your Workshop Supplement:
Workshop 1A: Chaboche
28