Porous Plasticity
Use the Porous Plasticity subnode to define the properties of a plasticity model for a porous material.
The Porous Plasticity node is only available with some COMSOL products (see https://www.comsol.com/products/specifications/). The material model is available for 3D, 2D, and 2D axisymmetry.
Porous Plasticity Model
Use this section to define the plastic properties of the porous material.
Material Model
Select the Material model 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:
Select an Effective void volume fractionBilinear, Asymptotic, or User defined.
For Bilinear and Asymptotic, enter the following data:
For User defined, enter the effective void volume fraction as a function of for example the void volume fraction, variable <phys>.f.
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.
See also Porous Plasticity, Elliptic Cap, and Elliptic Cap with Hardening in the Structural Mechanics Theory chapter.
Void Nucleation and 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 Nucleation and Growth for details.
When Include void nucleation in tension is selected, enter the Void volume fraction of nucleating voids, the Standard deviation for void nucleation, and the Mean strain for void nucleation. For each property use the value From material or enter a User defined value or expression.
When Include void growth in shear is selected, enter the Void growth rate parameter. Use the value From material or enter a User defined value or expression.
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
Nonlocal Plasticity Model
Nonlocal plasticity can be used to facilitate for example the modeling of material softening. Typical examples that involve material softening are finite strain plasticity and soil plasticity. In these situations, standard (local) plasticity calculations reveal a mesh fineness and topology dependence, where a mesh refinement fails to produce a physically sound solution. Nonlocal plasticity adds regularization to the equivalent plastic strain, thereby stabilizing the solution.
Select None, or Implicit Gradient. By selecting Implicit Gradient, regularization can added the equivalent matrix plastic strain and the void volume fraction. These two hardening variables describe different mechanisms and therefore often require different model parameters.
If Equivalent matrix plastic strain is selected, enter values for the:
Nonlocal coupling modulus, Hnl, which can be seen as a penalization of the difference between the local and nonlocal equivalent matrix plastic strains. A larger value will force the local equivalent matrix plastic strain to be closer to the nonlocal variable.
If Void volume fraction is selected, enter a value for the Length scale, void volume fraction, lint,f.
The two length scales lint,m and lint,f can be chosen independent of each other, but should not exceed the mesh size. It is also not recommended to define them as functions of the mesh size.
See also Nonlocal Plasticity in the Structural Mechanics Theory chapter.
Discretization
This section is available with the Implicit gradient nonlocal plasticity model. Select the shape function for the Nonlocal equivalent matrix plastic strain εpm,nl and the Nonlocal void volume fraction fnl Automatic, Linear, Quadratic Lagrange, Quadratic serendipity, Cubic Lagrange, Cubic serendipity, Quartic Lagrange, Quartic serendipity, or Quintic Lagrange. The options available depends on the chosen order of the displacement field.
To display this section, click the Show More Options button () and select Advanced Physics Options in the Show More Options dialog box.
Advanced
Enter the Maximum damage, which defines the residual stiffness of the model. The default value is 0.995.
Select the Local method to solve the plasticity problem — Automatic or Backward Euler. When Backward Euler is selected, it is possible to specify the maximum number of iterations and the relative tolerance used to solve the local plasticity equations. Enter the following settings:
Maximum number of local iterations. To determine the maximum number of iteration in the Newton loop when solving the local plasticity equations. The default value is 25 iterations.
Relative tolerance. To check the convergence of the local plasticity equations based on the step size in the Newton loop. The final tolerance is computed based on the current solution of the local variable and the entered value. The default value is 1e-6.
To display this section, click the Show More Options button () and select Advanced Physics Options in the Show More Options dialog box.
See also Numerical Solution of the Elastoplastic Conditions 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: