User login

Navigation

You are here

1D Plasticity - Isotropic hardening, tutorial with examples.

yawlou's picture

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

Comments

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

  • For
    isotropic hardening, hardening accumulates if plastic strain is positive
    or negative. How?

 

  • When
    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?

 

  • For
    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?

 

  • Hence
    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?

 

 

  • It
    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?

 

 

yawlou's picture

Answers are provided by the questions. 

  • For
    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.
    • When
      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.

    • For
      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.

    • Hence
      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.

    • It
      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 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";}

    Isotropic
    hardening:

    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.

    Kinematic
    hardening

    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.

    yawlou's picture

    Kajal,

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

    Isotropic
    hardening:

    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.

    Kinematic
    hardening

    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

    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";}

    Von
    Mises yield surface

    What I would like to understand is;

    • Let us consider Von Mises
      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.
    • If we take the square of
      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.
    • Now- for isotropic hardening
      (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.
    • For
      kinematic 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.

    Questions

    • I
      dot 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?

     

     

     

     

    Weijie Liu's picture

    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)?

     

    yawlou's picture

    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

     

     

    mohamedmoussa's picture

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

    Subscribe to Comments for "1D Plasticity - Isotropic hardening, tutorial with examples."

    More comments

    Syndicate

    Subscribe to Syndicate