Back to the Granular-Volcano-Group Home Page Granular-Volcano-Group logo
The ultimate website for understanding granular flows

A Review of Plastic-Frictional Theory
Part. 2
Plastic Potential Theory

You will find the basic facts about Plastic-Frictional Theories (Part. 2) - no details -. If you wanna know more just email me or feel free to ask in the Discussion Forum. I purposely erased all the bibliographical references and detailed equations to keep the text simple and easy to read.

If you need an official reference for the content of this website, please, use:
Dartevelle, S., Numerical and granulometric approaches to geophysical granular flows, Ph.D. thesis, Michigan Technological University, Department of Geological and Mining Engineering, Houghton, Michigan, July 2003.

We have seen on the preceding sections:
I. Introduction
II. Stress space, Slip Planes, Mohr-Coulomb and von Mises stresses
      II.1. Mohr-Coulomb case: a 2D representation of stress (particular case)
      II.2. von Mises case: a 3D representation of stress (general case)

On this page, you will find:
III. Plastic Potential Theory

And on the following pages, you will find the continuation of this frictional course:
IV. Dilatation, consolidation, yield locus, and critical state
      IV.1. Dilatancy
      IV.2. Consolidation
      IV.3. Yield locus and critical state
V. Constitutive Equations for Frictional Granular Flows
      V.1. Beyond the extended von Mises Yield locus
      V.2. Within the extended von Mises Yield locus

Please, don’t forget to sign our Granular Volcano Guest Book.

Done and updated by Sébastien Dartevelle, The WebMaster, Saturday November 22 2003.

III. Plastic Potential Theory

 

The yield condition laws I have reviewed in the preceding section says nothing about the nature of the motion which is initiated at yield. Alone, the Mohr-Coulomb/von Mises relations are helpless for describing a granular media in a continued motion. Constitutive relations are needed for describing the kinematics of motion at yield following the Mohr-Coulomb/von Mises theory. In addition, I must also account for the famous kinematics feature described for the first time by Osborne Reynolds: dilatancy. For a given normal stress, the material will sheared only if a critical shear is attained, once sheared, the granular medium may initially expand. Afterwards, it will settle to a fairly constant volume. This counter-intuitive behavior is of a key importance in soil mechanics.

 

The plastic potential theory will provide us a way to predict the velocity distribution within the granular medium at yield. This theory makes the connection between stress and deformation (or velocity gradient) in using three concepts: a Yield Function (Y), a Plastic Potential Function (G), and a Flow Rule.

 

The yield function (Y) is what we have seen in the previous section (Part. 1). For instance, from the Mohr-Coulomb and the von Mises yield functions, we have:

 

 

      Eq. 1

 

 

where IIdT is the second invariant of the deviator of the stress tensor T. The other symbols have been seen in the 1st Part of this course. IIdT can be expressed as:

 

 

      Eq. 2

 

 

It is clear from all I said in the previous section (Part.1 ), at yield, we must have Y=0 for a plastic material (as indicated by Eq.1). Therefore, at yield, we can see that Eq.1a is the equation of a straight line (which is, in the mathematical terminology, a level curve for Y=0, see Fig.8) in the principal stress plane (), and Eq.1b is the equation of a cone (which is, in the mathematical terminology, a level surface for Y=0, see Fig.9) in the principal stress space (). Both functions have an central axis of symmetry which has the property of being the axis of equal principal stresses, i.e., for Mohr-Coulomb 2D case and for the von Mises 3D case. Hence, this axis of symmetry represents the hydrostatic/isotropic Pressure. Both functions have their vertex on this hydrostatic axis at the origin, when .

 


[Figure 8: Domain of no-deformation (rigid) and domain of plastic deformation for the 2D Mohr-Coulomb case.
The central axis on which the principal stresses are equal represents the hydrostatic Pressure]

[Figure 9: Representation of the 3D von Mises Yield Surface in the principal stress space.
The central axis of symmetry on which all the principal stresses are equal is the hydrostatic/isotropic Pressure.
Plastic deformation occurs at yield on the surface of this cone]

 

The Plastic Potential Function (G) is defined so that the strain rate in any arbitrary directions (Di j) is proportional to the derivative of G with respect to the corresponding stress (Ti j):

 

 

      Eq. 3

 

 

where q is a positive scalar sometimes named "plastic multiplier". This scalar is not a property of the material but rather a property of the flow conditions. Eq.3 can also be written in terms of the principal directions of D and T (using their respective eigenvalues) as:

 

 

      Eq. 4

 

 

Now, we are almost there, since we now specifically play with the rate of deformation (D), hence the velocity gradient. So, the remaining thing to do is to find out what is G exactly? This is given by the flow rule. One of the most widely used flow rule is the associated flow rule, which states that the Plastic Potential Function (G) is equivalent to the Yield function (Y), G=Y, hence, Eq. 4 can be rewritten as:

 

 

      Eq. 5

 

 

It is worth saying that there are other flow rules which do not assume G=Y, those are named non-associated (I won’t comment this any further). It is clear that in Eq. 5, the Right-Hand-Side (RHS) is nothing but the gradient of the Yield Function with respect to the stresses, which is important to notice for defining the properties of the plastic potential flow. But before doing so, let’s indeed derived Eq. 5 using the Yield Function defined by Eq. 1b (extended von Mises) and the definition of the second invariant of the deviator of the stress tensor (IIdT) in Eq. 2c, we have in 3D:

 

 

      Eq. 6

 

 

This equation can be safely generalized as:

 

 

      Eq. 7

 

 

where D, I, and are the rate-of-strain tensor, the unit tensor and the deviator of the stress tensor respectively. The deviator of the stress tensor is:

 

 

      Eq. 8

 

 

where is the spherical part of the stress tensor which is the average in the 3-direction of space of the normal stresses (i.e., one third of the sum of the diagonal element of the stress tensor).

 

From those equations, we can now state the two key properties of the Plastic Potential Theory:

 

1- the co-axiality or alignment condition: the principal axe of rate of deformation are aligned with those of stresses. This to be in agreement with the intuitive idea that the material should respond to unequal stresses by contracting in the direction of greater stress and expand in the direction of lesser stress. In other words, it is equivalent to say that there is no shear strain on planes on which there is no shear stress. This is exactly what shows Eq. 6 and Eq. 7 since whenever the shear stress (Ti j) is zero on a given plane, so is the shear strain (Di j). I will come back to this property latter in this course.

 

2- the normality condition: which is shown by the original equation of the Plastic Potential Theory (Eq. 5). Indeed, if we consider the particular level surface Y=0 (at yield) in Eq. 5, we know from the properties of gradient that Di must be a vector perpendicular to the Yield (level) Surface, in the principal stress space.

 

Now, I’m pretty so you’re gonna ask: Sooooooooooooo what? Why the hell in the world do we want to know all that? Well, those last equations actually give us a connection between the velocity gradient within the granular flow and the state of stress within it. And basically, that’s what we were looking for, because it can be easily implemented into momentum equations in a computer model as we must know such relations.

 

Now, it would be interesting to calculate the divergence of the velocity field, which is nothing but minus the volumetric rate of deformation of the granular body (i.e., first invariant of D):

 

 

          Eq. 9

 

 

Therefore, with both Mohr-Coulomb and von Mises cases, the granular body, at yield, suffers dilatation (divergence of the velocity field is positive, i.e., that is extension). On the one hand, that is good because it is in agreement with Reynolds' Principle of Dilatancy and therefore that is what we expect from those laws. However, on the other hand, it is bad as Eq. 9 leads to a continued expansion of the granular medium without any limits, which is physically impossible.

 

But there is another trouble. Indeed, dilatancy is one phenomenon that can be observed in granular material at yield but there are other processes that may happen as well depending on the physical properties of the granular material and on its previous deformation history. If the granular medium is well-compacted, it will expand when sheared (that is dilatancy), whereas loosely packed, a granular medium at large normal stress may contract (that is consolidation or contractancy). Moreover, the expansion or the contraction takes place only until a critical bulk density is reached, afterwards under continuous stress, deformation will continue without any change in the bulk density. And that is not shown at all in the Eq. 9. In addition, it can be shown that the rate of energy dissipation is always zero if we use a simple Mohr-Coulomb/von Mises law, which is an unacceptable results for frictional processes.

 

So, we got a problem down here … In conclusion, the Mohr-Coulomb/von Mises theory cannot be applied straightforwardly for describing flowing frictional granular material without committing unphysical instabilities and inconsistencies. The simple Mohr-coulomb/von Mises law must be -somehow- adapted to the plastic potential flow theory for describing flowing granular medium accurately. But before getting into the constitutive equations and the possible corrections, it is necessary to review the basic behaviors of granular material at yield. This is what I’m gonna do in this following section:

 

~~~~~~~~~ Part 3: IV. Critical State Theory ~~~~~~~~~

 

If you have enjoyed this Plastic-Frictional (Part 2) webpage, pleaaase, before you leave, sign my Guestbook. It's all I ask!

[Back to Top]

Other Granular Volcano Group Webpages:

- What is a Granular Medium?

A complete review of Plastic-Frictional theories:

   - Part 1. Introduction, Mohr-Coulomb, and von Mises Stresses

   - Part 2. Plastic Potential Theory

   - Part 3. Critical State Theory

   - Part 4. Constitutive Equations for frictional granular flow

- All I wanna know about viscous stress!

- Granular Theory: an Overview

- Compute Your Own Atmospheric Profile

Numerical Results:

   - Go to the Numerical Results Page (Introduction and all Results)

   - Go to the Plinian Cloud simulations Page

Welcome Page

Home Page

You may also enjoy those pages:

| Guestbook | Discussion Forum | Chatroom | Volcano Links | Granular Links | Awards | Logos | Who is the Webmaster? | Ph.D. |

_____________________________________________________

Any problems or comments with this web site can be reported to the Web Master.