Porous Plasticity
Use the Porous Plasticity subnode to define the properties of a plasticity model for a porous material.
Porous Plasticity Model
Use this section to define the plastic properties of the porous material.
Plasticity Model
Select Small plastic strains or Large plastic strains to apply either an additive or multiplicative decomposition between elastic and plastic strains.
Yield Function F
The Yield function F defines the limit of the elastic regime F(σ, σys) ≤ 0.
Select a Yield function F for the porous plasticity criterion — Shima-Oyane, Gurson, Gurson-Tvergaard-Needleman, Fleck-Kuhn-McMeeking, FKM-GTN, or Capped Drucker-Prager.
Shima–Oyane
For Shima-Oyane enter the following data:
Gurson
For Gurson enter the following data:
Gurson–Tvergaard–Needleman
For Gurson-Tvergaard-Needleman enter the following data:
Fleck–Kuhn–McMeeking
For Fleck-Kuhn-McMeeking enter the following data:
FKM–GTN
For FKM-GTN enter the following data:
Capped Drucker–Prager
For Capped Drucker-Prager enter the following data:
The material properties use values From material (default) or User defined.
Void Growth
It is possible to Include void nucleation in tension or Include void growth in shear by selecting the corresponding check box. See the section Void Growth for details.
Isotropic Hardening Model
Select the type of linear or nonlinear isotropic hardening model from the Isotropic hardening model list.
Select Perfectly plastic (ideal plasticity) if the material can undergo plastic deformation without any increase in yield stress. When Capped Drucker-Prager is selected, enter values or expressions to define the semi-axes of the cap under Elliptic cap parameter pa and Elliptic cap parameter pb.
For Linear the default Isotropic tangent modulus ETiso uses values From material (if it exists) or User defined. The flow stress (yield level) σfm is modified as hardening occurs, and it is related to the equivalent plastic strain in the porous matrix εpm as
with
For the linear isotropic hardening model, the flow stress (yield stress) increases proportionally to the equivalent plastic strain in the porous matrix εpm. The Young’s modulus E is taken from the elastic material properties.
Select Ludwik from the list to model nonlinear isotropic hardening. The flow stress (yield level) σfm is modified by the power-law
The Strength coefficient k and the Hardening exponent n use values From material (if it exists) or User defined.
For Power law isotropic hardening, the Hardening exponent n uses the value From material (if it exists) or User defined. The flow stress (yield level) σfm is modified by the power-law
for
The Young’s modulus E is taken from the elastic material properties.
For Hardening function, the isotropic Hardening function σh(εpm) uses values From material or User defined. The flow stress (yield level) σfm is modified as
-
This definition implies that the hardening function σh(εpm) in the Material node must be zero at zero plastic strain. In other words, σfm = σys0 when εpm = 0. With this option it is possible to enter any nonlinear isotropic hardening curve. The hardening function can depend on more variables than the equivalent plastic strain in the porous matrix, for example the temperature. Select User defined to enter any function of the equivalent plastic strain εpm. The variable is named using the scheme <physics>.<elasticTag>.<plasticTag>.epm, for example, solid.lemm1.popl1.epm.
For Exponential hardening, the cap in the Capped Drucker-Prager model evolves with the volumetric strain. Since the volumetric plastic strain εpvol is negative in compression, the limit pressure pb in the cap increases from pb0 as hardening evolves
The Isotropic hardening modulus Kiso, the Maximum plastic volumetric strain εpvol,max and the Ellipse aspect ratio R use values From material (if it exists) or User defined. Enter a value or expression to define the initial semi-axis of the ellipse under the Initial location of the cap pb0.
See also Porous Plasticity, Elliptic Cap, and Elliptic Cap With Hardening in the Structural Mechanics Theory chapter.
To compute the energy dissipation caused by porous compaction, enable the Calculate dissipated energy check box in the Energy Dissipation section of the parent material node (Linear Elastic Material or Nonlinear Elastic Material).
Location in User Interface
Context Menus
Ribbon
Physics tab with Linear Elastic Material or Nonlinear Elastic Material node selected in the model tree: