## You are here

# 1D Plasticity - Isotropic hardening, tutorial with examples.

For those who may be interested,

I have put together a paper describing 1D plasticity for a variety of cases of isotropic hardening. The material is not new, but hopefully it is written with enough detail that it will help beginners learn some basics of computational plasticity. The notation and material closely follows "Computational Inelasticity" written by Simo and Hughes.

Theory and derivations are carefully provided trying to include all steps. Also, five different isotropic hardening algorithms are summarized (no hardening, linear hardening, quadratic hardening, Ramberg-Osgood hardening, Voce Hardening). A numerical algorithm and examples are provided for 1D plasticity incorporated into a standard 2D truss program.

I hope this is helpful to those striving to learn about plasticity. The paper can be downloaded here:

http://people.wallawalla.edu/~louie.yaw/plasticitypublications/1Dplastic...

regards,

Louie

- yawlou's blog
- Log in or register to post comments
- 33391 reads

## Comments

## Dear Sir,I found this

Dear Sir,

I found this paper very useful and shall be grateful if you can help me with these questions. I have been reading it repeatedly and trying to follow things.The questions are from the first,second and third sections only and address very basic issues

isotropic hardening, hardening accumulates if plastic strain is positive

or negative. How?

hardening occurs, the value of the yield stress changes. Does it mean that

when hardening occurs, elastic strain increases which means that when

hardening occurs value of yield stress increases?

isotropic hardening, hardening accumulates if plastic strain is positive

or negative, hence internal hardening variable keep track of the total

change in plastic strain.

Question:

Ø

When will the plastic strain become negative?

Is it that the internal hardening variable will become negative because yield

stress (hence elastic strain) increases? If so, does it mean that plastic

strain may be negative when the material has just entered the plastic regime?

internal hardening variable "α" keeps a track of the total change in

plastic strain.

Question:

Question:

Yield condition (f) is defined as current

level of stress minus the initial yield stress. When you say-initial yield

stress, do you mean when the yielding began?

Question:Does a negative value of ‘f' always denote

that the material has not started yielding- that is, it is still in the elastic

regime?

Question:

Is ‘f' the stress level at a particular

strain?

has been stated that: if ‘f' is positive it means that current level of

stress is above the yield stress- elastic and plastic strains are

increasing

First

question:

Ø

Elastic strain will always increase-correct?

Second

question:

Ø

Does it mean that when f>0, we need to

calculate the flow stress using the material law we are considering?

## Answers are provided by the

Answers are provided by the questions.

isotropic hardening, hardening accumulates if plastic strain is positive

or negative. How? Material tests show that isotropic hardening happens. If a specimen in tension yields and hardens the same specimen loaded in compression will have a higher yield stress because of the hardening that previously took place when it was loaded in tension. This is due to isotropic hardening. The material has a "memory" (it has been affected by tension yielding) and this affects how it responds when it is loaded in compression.

hardening occurs, the value of the yield stress changes. Does it mean that

when hardening occurs, elastic strain increases which means that when

hardening occurs value of yield stress increases?

yes.

Suppose the initial yeild stress was 36 ksi. After hardening has taken place suppose the material stress reaches 40 ksi. Then if you unload the specimen to 0 stress then loaded it again, it will not yield until reaching 40 ksi. Hence the yield stress has increased from 36 initially to 40 ksi due to hardening. As soon as 40 ksi is passed and hardening continues the yield stress will again increase.

isotropic hardening, hardening accumulates if plastic strain is positive

or negative, hence internal hardening variable keep track of the total

change in plastic strain.

yes

Question:

Ø

When will the plastic strain become negative?

Is it that the internal hardening variable will become negative because yield

stress (hence elastic strain) increases? If so, does it mean that plastic

strain may be negative when the material has just entered the plastic regime?

Based on your choice of coordinate system often tension is positive for plastic strain and compression produces negative plastic strain.

No, the internal hardening variable alpha is always positive, it does not ever go negative.

internal hardening variable "α" keeps a track of the total change in

plastic strain.

yes

Question:

Question:

Yield condition (f) is defined as current

level of stress minus the initial yield stress. When you say-initial yield

stress, do you mean when the yielding began?

f is defined as current level of absolute stress minus (initial yield stress plus any change in yield stress due to hardening).

Initial yield stress is the stress that material will yield at the first time it is loaded up to yield.

Question:Does a negative value of ‘f' always denote

that the material has not started yielding- that is, it is still in the elastic

regime?

Negative f means the material is in the elastic range, but it does not mean that it hasn't yielded previously. It might have yielded previously and then it might be unloaded below yield and be in the elastic range and hence have f<0.

Question:

Is ‘f' the stress level at a particular

strain?

no, it is a function that tries to mathematically describe the boundary between elastic strain causing stresses and plastic strain causing stresses.

has been stated that: if ‘f' is positive it means that current level of

stress is above the yield stress- elastic and plastic strains are

increasing

No, be careful here. f positive means that plastic flow is taking place. It does not mean stress is above the yield stress. This is a subtle point, but important. Mathematically we may calculate that f is positive, we might be inclined to think then that our stress is above yield, but that isn't quite the right way to think of it. Really what it means is that plastic flow is taking place if f is positive and we need to determine what the amount of plastic flow(strain) has happened such that f is at 0. Hence, f allows us to determine the amount of plastic flow and hardening that has taken place. We do this by using the kuhn tucker conditions, the consistency parameter and the consistency condition (all the equations and derivations indicated in the paper).

First

question:

Ø

Elastic strain will always increase-correct?

No, when the specimen is being unloaded elastic strain is decreasing. When zero stress is applied there may be permanent nonzero plastic strain remaining, but the elastic strain will be zero. Hence, clearly elastic strain can increase or decrease.

Second

question:

Ø

Does it mean that when f>0, we need to

calculate the flow stress using the material law we are considering?

When f>0 we need to find the consistency parameter that makes f=0 again. We don't want f>0, f is only a mathematical device to allow us to calculate the value of the consistency parameter. Then using the consistency parameter we follow our material law for plastic flow and hardening and update internal variables appropriately.

## Dear Sir,Thank you for

Dear Sir,

Thank you for the response and I have read your paper.I would be grateful if you can correct me if wrong on the following which I feel I should get it right before going into hard math stuff.This is about fundamentals on isotropic and kinematic hardening and my understanding.

800x600

Normal

0

false

false

false

EN-IN

X-NONE

X-NONE

MicrosoftInternetExplorer4

/* Style Definitions */

table.MsoNormalTable

{mso-style-name:"Table Normal";

mso-tstyle-rowband-size:0;

mso-tstyle-colband-size:0;

mso-style-noshow:yes;

mso-style-priority:99;

mso-style-parent:"";

mso-padding-alt:0cm 5.4pt 0cm 5.4pt;

mso-para-margin:0cm;

mso-para-margin-bottom:.0001pt;

mso-pagination:widow-orphan;

font-size:10.0pt;

font-family:"Calibri","sans-serif";

mso-bidi-font-family:"Times New Roman";}

Isotropichardening:

For isotropic hardening, if you plastically

deform a solid, then unload it, then try to reload it again, you will find that

its yield stress (or elastic limit) would have increased compared to what it

was in the first cycle.

Again, when the solid is unloaded and

reloaded, yield stress (or elastic limit) further increases.

This continues until a stage (or a cycle) is

reached that the solid deforms elastically throughout.

This is isotropic hardening.

Essentially, in isotropic hardening the yield

stress (or elastic limit) increases whereas ductility keeps decreases. Hence,

in the nth cycle the solid deforms elastically.

Is this understanding of isotropic hardening

correct Sir ?

I shall be grateful if corrected.Kinematichardening

Isotropic hardening is not useful in

situations where components are subjected to cyclic loading.

Isotropic hardening does not account for

Bauschinger effect and predicts that after a few cycles, the material (solid)

just hardens until it responds elastically

To fix this, alternative laws i.e. kinematic

hardening laws have been introduced

As per these hardening laws, the material

softens in compression and thus can correctly model cyclic behaviour and

Bauschinger effect.

Is this undestanding of kinematic hardening correct Sir?I shall be grateful if corrected.## Kajal, [Answers in

Kajal,

[Answers in brackets] to your questions are provided by the questions below.

Isotropichardening:

For isotropic hardening, if you plastically

deform a solid, then unload it, then try to reload it again, you will find that

its yield stress (or elastic limit) would have increased compared to what it

was in the first cycle. [YES]

Again, when the solid is unloaded and

reloaded, yield stress (or elastic limit) further increases. [YES, if it is reloaded past its previously reached maximum stress.]

This continues until a stage (or a cycle) is

reached that the solid deforms elastically throughout. [ I think I understand your statement, eventually if the cycles of load are always to the same level, then after just one cycle your specimen on subsequent cycles will just be loading and unloading along the elastic line of the stress strain curve. Note also that, a steel bar with 1 pound of load on it has deformed elastically thoughout.]

This is isotropic hardening. [Isotropic hardening just means if you load something in tension past yield, when you unload it, then load it in compression, it will not yield in compression until it reaches the level past yield that you reached when loading it in tension. In other words if the yield stress in tension increases due to hardening the compression yield stress grows the same amount even though you might not have been loading the speciment in compression. It is a type of hardening used in mathematical models to describe plasticity. It might not be absolutely correct for real materials.]

Essentially, in isotropic hardening the yield

stress (or elastic limit) increases whereas ductility keeps decreases. Hence,

in the nth cycle the solid deforms elastically.[Your last sentence here "Hence, in the nth cycle the solid deforms elastically." i'm not sure what you mean. The previous sentences here are ok. As hardening increases, ductility is being used up, or is decreasing. Although, your sentence might make some sense if your cycles are to the same maximum stress each time, in which case after one cycle the material is just being loaded and unloaded along the elastic line of the stress strain curve.]

Is this understanding of isotropic hardening

correct Sir ? [see comments above, you have it mostly right.]

I shall be grateful if corrected.Kinematichardening

Isotropic hardening is not useful in

situations where components are subjected to cyclic loading.[I believe real metals exhibit some isotropic hardening AND some kinematic hardening(I'm not sure if this is absolutely correct). Hence Isotropic hardening is useful. But, kinematic hardening is apparently more realistic for cyclic loading of metals.]

Isotropic hardening does not account for

Bauschinger effect and predicts that after a few cycles, the material (solid)

just hardens until it responds elastically [Yes, I think what you mean is correct.]

To fix this, alternative laws i.e. kinematic

hardening laws have been introduced[Yes]

As per these hardening laws, the material

softens in compression and thus can correctly model cyclic behaviour and

Bauschinger effect.

Is this undestanding of kinematic hardening correct Sir?[I think you have it.]regards,

Louie

## Thanks to Louie Yaw Sir...

Dear Sir,

Many thanks for the reply.My profound gratitude.Thank you for making things clearer for me.

Just going a bit further- realting isotropic and Kinematic hardening to Von Mises yield surface:

800x600

Normal

0

false

false

false

EN-IN

X-NONE

X-NONE

MicrosoftInternetExplorer4

/* Style Definitions */

table.MsoNormalTable

{mso-style-name:"Table Normal";

mso-tstyle-rowband-size:0;

mso-tstyle-colband-size:0;

mso-style-noshow:yes;

mso-style-priority:99;

mso-style-parent:"";

mso-padding-alt:0cm 5.4pt 0cm 5.4pt;

mso-para-margin:0cm;

mso-para-margin-bottom:.0001pt;

mso-pagination:widow-orphan;

font-size:10.0pt;

font-family:"Calibri","sans-serif";

mso-bidi-font-family:"Times New Roman";}

VonMises yield surface

What I would like to understand is;

yield criteria. Now, Von Mises yield criteria is a yield criteria which states

that you cannot have stress states outside the surface defined by the criteria.

The surface does not have fixed dimensions except for perfect plasticity.

both sides of Von Mises criteria, we notice that we get an equation of infinite

cylinder in 3D (axis being the principal stresses). The cylinder is the yield

surface, the radius of the cylinder is the yield stress.

(not kinematic hardening) the correct material law is the one (or the correct

material modelling) is the one in which the yield surface increases in size but

remains the same shape as a result of plastic straining.

Forkinematic hardening, the correct material law is the one in which allows the

yield surface to translate without changing its shape. So- as one deforms the

material in tension, you drag the yield surface thus increasing the stress thus

modelling strain hardening. This softens the material in compression.

QuestionsIdot follow that for - how is it possible that if the yield surface i.e.

cylinder if dragged (i.e. translated) the material hardens in tension but

softens in compression?

Now few questions on your paper of 1D plasticity which is turning out to be so useful for me undestand things and concepts:

Normal

0

false

false

false

EN-GB

X-NONE

X-NONE

/* Style Definitions */

table.MsoNormalTable

{mso-style-name:"Table Normal";

mso-tstyle-rowband-size:0;

mso-tstyle-colband-size:0;

mso-style-noshow:yes;

mso-style-priority:99;

mso-style-parent:"";

mso-padding-alt:0cm 5.4pt 0cm 5.4pt;

mso-para-margin-top:0cm;

mso-para-margin-right:0cm;

mso-para-margin-bottom:10.0pt;

mso-para-margin-left:0cm;

line-height:115%;

mso-pagination:widow-orphan;

font-size:11.0pt;

font-family:"Calibri","sans-serif";

mso-ascii-font-family:Calibri;

mso-ascii-theme-font:minor-latin;

mso-hansi-font-family:Calibri;

mso-hansi-theme-font:minor-latin;

mso-bidi-font-family:"Times New Roman";

mso-bidi-theme-font:minor-bidi;

mso-fareast-language:EN-US;}

1.

What

does the consistency parameter physically denote? Does it denote amount of

plastic flow? Why is that consistency parameter be always greater than zero?

On page 4, you have mentioned the Kuhn

Tuckers condition

2.

Kuhn

Tucker's conditions are conditions (or rules) used to construct mathematical

algorithms which model the process of plastic flow. One of Kuhn Tuckers

condition is : γ >=0.

Now,

γ is the consistency parameter which allows a means for determining the level

of plastic flow and hardening such that f<=0 is satisfied. What is the

relevance of having : γ >=0. What is the physics/math involved here.

3.

Again,

what is the relevance of Kuhn Tuckers condition γ f(σ) = 0. How is the

evaluation justified mathematically and physically?

In box 4.1, which is 1D plasticity

algorithm for general isotropic hardening as well as box 5.1which is 1D

plasticity algorithm with no hardening (perfect plasticity), it has been stated

in step 3 to compute the elastic trial stress,

the trial value for the yield function and test for plastic loading.

4. What is this trial elastic stress actually? Because, you have mentioned

the expression for trial elastic stress as:

sigma_n+1(trial) = E (epsilon_n+1 - epsilon_n^p)

Also you have marked epsilon_n+1 and epsilon_n^p in Figure 2 of the paper.

Normal

0

false

false

false

EN-GB

X-NONE

X-NONE

/* Style Definitions */

table.MsoNormalTable

{mso-style-name:"Table Normal";

mso-tstyle-rowband-size:0;

mso-tstyle-colband-size:0;

mso-style-noshow:yes;

mso-style-priority:99;

mso-style-parent:"";

mso-padding-alt:0cm 5.4pt 0cm 5.4pt;

mso-para-margin-top:0cm;

mso-para-margin-right:0cm;

mso-para-margin-bottom:10.0pt;

mso-para-margin-left:0cm;

line-height:115%;

mso-pagination:widow-orphan;

font-size:11.0pt;

font-family:"Calibri","sans-serif";

mso-ascii-font-family:Calibri;

mso-ascii-theme-font:minor-latin;

mso-hansi-font-family:Calibri;

mso-hansi-theme-font:minor-latin;

mso-bidi-font-family:"Times New Roman";

mso-bidi-theme-font:minor-bidi;

mso-fareast-language:EN-US;}

What does εn+1- εnp

denote? I’m unable to interprate εn+1- εnp from this graph?

## Dear Sir: This paper is

Dear Sir:

This paper is very useful. And I shall be very grateful if you can help me with these questions.

1. "For the 1D plasticity case the algorithmic tangent modulus is equivalent to the elasto-plastic modulus. In higher dimensions this is not true."

I cann't distinguish the two "algorithmic tangent modulus" and "the elasto-plastic modulus";

2. For 3D plasticity with general linear isotropic harding. Consider the Mises yiled function.

When the SIMO's renturn method converge, the relationship of

dσ=[Dsp] * dε should be satisfied? [Dsp] is the elasto-plastic module.

σn+1=σn + [Dsp] * dε, should be consistent with stress calculated by SIMO's Return Method (f = 0, flow rules)?

## answers to your plasticity questions

Hello Weijie,

I will answer your questions in the order that you asked them above.

1. The elasto-plastic modulus comes about from a purely mathematical(not within an incremental algorithm) derivation. (See Simo and Hughes book page 80, 81 for the general derivation and pages 90 and 91 for particular results for Isotropic hardening and Isotropic with Kinematic hardening, 3dimensionl case). These results are entirely "mathematical" and are NOT derived within the framework of an incremental algorithm that you would implement in a finite element program. The elasto-plastic modulus is basically dsigma/depsilon.

The algorithmic tangent modulus comes about from a derivation of dsigma/depsilon within an incremental plasticity algorithm. Such an algorithm is implemented with discrete increments (steps). The algorithmic tangent modulus in this case (3D case) is NOT the same as the purely mathematical elasto plastic modulus. Compare the results with Simo and Hughes on page 91 with the result in box 3.2 on page 124.

2. See also remark number 4 on page 125 of Simo and Hughes which answers your question number 2. As delta t and the consistency parameter approach zero the elasto plastic modulus and the algorithmic elasto plastic modulus become the same. But, for large time steps however they are different. Why is this important? The reason this is important is as follows. Suppose you write a finite element analysis program and want to implement a nonlinear material model for plasticity. For each increment of load that causes plastic flow you will need to impose the consistency condition to determine the increase in plastic strain and elastic strain. This is can be done by using Newton Raphson iterations at the local element level. However, you may also do Newton-Raphson iterations at the global level for your finite element model in order to achieve global equilibrium. In order for your global iterations to have the desireable quadratic rate of convergence (that Newton Raphson iterations can achieve if done correctly) you need to use the (algorithmic) elasto plastic tangent modulus that is derived consistently within the frame work of the incremental algorithm. If you do not do this you will not have a quadratic rate of convergence with your global Newton Raphson iterations(your iterations will take longer and may not even converge). So, if you are using Newton Raphson iterations at the global level for equilibrium, you definitely want to use the consistent algorithmic elasto plastic modulus at the local element level.

You may also benefit by looking at my paper for the 1D case and seeing how the two are derived (elasto plastic modulus equation 4.11 and the algorithmic elasto plastic modulus equation 4.31). It will help to see how these are derived differently. For the 1D case they end up being the same, but in 3 dimensions they do not turn out the same as I mentioned previously above.

One sometimes says that the algorithmic elastoplastic (tangent) modulus is variationally consistent within the plasticity algorithm when Newton Raphson iterations are used at the global level.

This is not easy to put into words, but that is the best I can do to explain this. It is a somewhat tricky concept to convey. Understanding this comes about from reading books and literature carefully and even better taking a class where the instructor can explain this.

I also want to say that I am still learning these concepts myself and if

some others, here on imechanica, want to add their thoughts on this, I

would be happy for your comments and corrections to what I have said.

I hope that helps,

Louie

## Thanks for posting these

Thanks for posting these notes up. They helped me get my head around plasticity modelling.