Mohr Coulomb and Drucker Prager Plasticity in Optistruct

solon
solon
Altair Employee
edited October 21 in Altair HyperWorks

In metal plasticity, typically only the deviatoric stress is the cause of yielding. Plotting the yield surface in principal stress space, will result in a cylinder whose axis is the hydrostatic pressure line. This indicates that yielding is independent of the hydrostatic stress state, which means that we could have infinite hydrostatic compression with no yielding occurring. On the other hand, rock plasticity theories also use a similar concept except that the requirement of pressure-dependence of the yield surface demands a relaxation of the above assumption.

For materials, like soils and rocks, the frictional and dilatational effects are significant. The plastic behavior depends on the hydrostatic pressure as the internal friction is proportional to the normal force. Two commonly used yield surfaces/failure criteria are the Mohr-Coulomb model and the Drucker-Prager model.

In Optistruct, Mohr Coulomb and Drucker Prager Plasticity is currently supported only for Solid elements in Implicit analysis (NLSTAT/NLTRAN for SMDISP/LGDISP).

 

1. Mohr Coulomb

Coulomb's theory is a mathematical framework that characterizes the behavior of brittle materials, such as concrete, soil, rocks under shear and normal stress conditions. Typically, it is relevant for materials where the compressive strength significantly surpasses the tensile strength.

The Mohr-Coulomb criterion posits that failure takes place when the shear stress at any point within a material attains a value that is linearly related to the normal stress acting on the same plane. This model relies on the graphical representation of Mohr's circle, which illustrates the stress states at failure corresponding to the maximum and minimum principal stresses. The failure line is defined as the optimal straight line that tangentially intersects these Mohr's circles.

Plastic yielding begins when, on a plane in the body, the shearing stress, image, and the normal stress, image, reach the critical combination:

image

Where image is cohesion and image is friction angle.

image

Mohr Coulomb Plasticity in principal stress space:

 

 

image

Yield Criteria:

image

Flow rule:

image

Where, image and image are the maximum and the minimum principal stress of a stress state, image is cohesion, image is friction angle, image is dilatancy angle.

Dilatancy angle image represents the ratio of plastic volume change over plastic shear strain.

When image, it is associative flow, otherwise, it is non-associate flow

When image 0.0, it becomes Tresca Plasticity

 

Optistruct Set Up:

Using MAT1+MATS1 with the Optional continuation line “MC” for Mohr-Coulomb plasticity model:

image

Where:

‘MC’:             FLAG for Mohr Coulomb plasticity

COHE:         Cohesion

FRICA:         Friction Angle (No default (Real >= 0.0 and =< 89.9))

DILA:           Dilatancy Angle (Default = blank = FRICA (Real >= 0.0 and =< FRICA))

 

2. Drucker Prager

Drucker-Prager plasticity model is suitable for confined concrete, rock, soil that has the following characteristics: pressure dependence, path dependence, non-associative flow rule, work or strain hardening, and limited tensile strength. In a plasticity model, the parameters related to friction angle and cohesion affect the hardening and yielding criteria, while the parameter related to plastic dilation determines the flow rule.

Drucker-Prager yield criterion uses a modified von Mises yield criteria to incorporate the effects of pressure for massive structures.

Linear Drucker-Prager law:

For Linear Drucker-Prager, when the yield surface is plotted in principal stress space, it would resemble a cone. This indicates that if there is hydrostatic tension, the yield strength would get smaller and if hydrostatic compression increases, the yield strength would get bigger.

A blue cone with yellow lines and a pointDescription automatically generated

 

 

The yielding function for the DP model is in the form of:

image

Where, q and p are von-Mises stress and mean stress of a stress state, respectively. And β and d are friction angle and cohesive, respectively. (These material constants are provided by the user).

The angle of internal friction β describes the slope of the yield surface. If β=0, this would imply no dependence on hydrostatic pressure – effectively, this would change the behavior to the Mises yield surface.

The flow rule of the DP plasticity follows the form of:

image

Where, φ is the dilatancy angle (The dilatancy angle φ describes the flow potential).

If β=φ, the flow is associative.

If φ=0, no inelastic volumetric strains will be produced.



A line with arrows and lettersDescription automatically generated with medium confidence

Yield Criteria in p-q plane

Optistruct Set Up:

Optional continuation lines for Linear Drucker-Prager plasticity model (only for Implicit analysis)

image

Where:

‘LDP’:          FLAG for Linear Drucker Prager model

TYPE:           type of YIELD

                     = ‘COMP’, compression yield stress

                     = ‘TENS’, tension yield stress

                     = ‘COHE’, cohesion

YIELD:          initial yield stress (or d)

FRICA:         Friction Angle β (No default (Real >= 0.0 and =< 89.9))

DILA:           Dilatancy Angle φ (Default = blank = FRICA (Real >= 0.0 and =< FRICA))

 

When TYPE = ‘COMP’, users need to provide image in YIELD field, and OS converts image to cohesion image by 

image

When TYPE = ‘TENS’, users need to provide image in YIELD field, and OS converts image to cohesion image by 

image

When TYPE = ‘COHE’, users provide image in YIELD field, and OS use it directly without conversion

 

Material constants conversion between MC and DP.

In plane strain condition, the material constants from Mohr-Coulomb model can be converted to Drucker-Prager model by the following formulation:

image

 

 

 

3. Optistruct Example - Bearing Capacity of Soil

 

image

 

The soil model will be 5m wide, 5m long and 0.1m thick.

The footing an area equal to AXZ= 0.5x0.1=0.05 m2.

Enforced Displacement of -0.005 m for the footing nodes. (the nodes at the area of the footing are connected via MPC)

A screen shot of a computerDescription automatically generated

Arc-Length method is used NLPCI (NLCTRL or NLPARM continuation line):

A screenshot of a computerDescription automatically generated

 

Example 1- Mohr Coulomb: MC_model.fem

Material Set -Up:

In MAT1 we define the Young’s Modulus to be equal to 107 kPa, poisson’s ratio as 0.49.

In MATS1 material card we activate the MC continuation line, and we define cohesion equal to 490000 kPa and friction angle as 20o. Since we leave dilatancy angle as default it’s value will be equal to the friction angle (associative flow rule).

A screenshot of a computerDescription automatically generated A screenshot of a computerDescription automatically generated

 

Results:

 

Equivalent plastic strain can be requested via STRAIN(PLASTIC) = ALL

Example 2- Drucker Prager: DP_model.fem

The parameters for Drucker Prager Plasticity are calculated using the following equations, since we enforced plane strain condition:

A math equations with numbers and linesDescription automatically generated with medium confidence

Material Set -Up:

A screenshot of a computerDescription automatically generated

Comparison of the two Placiticty Models:

image

image

 

image

image

 

A screenshot of a computerDescription automatically generated

Load Bearing Curve (SPCF-y vs DISP-y)

 

To plot the above curve in HyperView + Hypergraph:

  • In HV select the nodes of the footing and plot Displacement and SPC forces in y direction.

A screenshot of a computerDescription automatically generated

  • In HG select all curves containing the displacements and right click, Multiple Curves Math, Add, to summate the results in a single curve. Then, repeat the same for the curves with the SPCF results.

A screenshot of a computerDescription automatically generated

  • Go to Define, click on “+” sign to add a new curve, select as source Math and Expression Builder as Curves. Then select the y axis values of the Summated Displacement curve, as X component in the vectors and the y axis values of the Summated SPCF curve, as Y component in the vectors.

A screenshot of a computerDescription automatically generated