Elastoplastic Soil Material
The Elastoplastic Soil Material feature is used to model stress-strain relationships that are nonlinear even at infinitesimal strains. It is available in the Solid Mechanics interface. This material model requires a Geomechanics Module license (see https://www.comsol.com/products/specifications/). Elastoplastic Soil Material is available for 3D, 2D, and 2D axisymmetry.
By adding the following subnodes to the Elastoplastic Soil Material node you can incorporate other effects:
Add an External Stress node in case you need to define a pore pressure in a porous soil. The Pore pressure pA is user defined by default. The default value is 1 atm, but you can change it to another value or expression for the pore fluid pressure. If there are other physics interfaces (like Darcy’s Law) in the model that make a pressure variable available, such variables will be available in the list.
Coordinate System Selection
The Global coordinate system is selected by default. The Coordinate system list contains any additional coordinate systems that the model includes (except boundary coordinate systems). The coordinate system is used when stresses or strains are presented in a local system. The coordinate system must have orthonormal coordinate axes, and be defined in the material frame. Many of the possible subnodes inherit the coordinate system settings.
Elastoplastic Soil Material
Select a Material model from the list: Modified Cam-Clay, Modified Structured Cam-Clay, Extended Barcelona Basic, Hardening Soil, or Hardening Soil Small Strain.
Density
All elastoplastic soil models have density as an input. The default Density ρ uses values From material. For User defined enter another value or expression.
If the material has a temperature dependent density or other material property, and From material is selected, the Volume reference temperature list will appear in the Model Input section. As a default, the value of Tref is obtained from a Common model input. You can also select User defined to enter a value or expression for the reference temperature locally.
Default Model Inputs and Model Input in the COMSOL Multiphysics Reference Manual.
Modified Cam-clay
The Modified Cam-Clay options adds the equations and interface for defining the material properties for the modified Cam-clay soil model.
From the Specify list, define the elastic properties either in terms of Poisson’s ratio or Shear modulus.
The defaults for the Poisson’s ratio ν or the Shear modulus G, Density ρ, Slope of critical state line M, Swelling index κ, Compression index λ, Initial void ratio e0 , and Void ratio at reference pressure eref are taken From material. For User defined enter other values or expressions
Enter a value or expression for the Reference pressure pref, and the Initial consolidation pressure pc0.
For the Slope of critical state line you can alternatively select Match to Mohr–Coulomb criterion or Match to Matsuoka–Nakai criterion, which then matches the slope of the virgin consolidation line to the angle of internal friction. Select the Angle of internal friction ϕ as From material or User defined.
For the Initial void ratio, you can alternatively select From void ratio at reference pressure, which then computes the initial void ratio from the reference pressure, the void ratio at reference pressure, the initial consolidation pressure, and the swelling and compression indexes. Then select the Void ratio at reference pressure eref as From material or User defined.
See also The Modified Cam-Clay Soil Model in the Structural Mechanics Theory chapter.
Isotropic Compression with Modified Cam-Clay Material Model: Application Library path Geomechanics_Module/Verification_Examples/isotropic_compression
Modified Structured Cam-Clay
From the Specify list, define the elastic properties either in terms of Poisson’s ratio or Shear modulus.
The defaults for the Poisson’s ratio ν or the Shear modulus G, Density ρ, Slope of critical state line M, Swelling index for structured clay κs, Compression index for destructured clay λd, Angle of internal friction ϕ, Initial structure strength pbi, Destructuring index for volumetric deformation dv, Destructuring index for shear deformation ds, Plastic potential shape parameter ξ, Initial void ratio e0, Void ratio at reference pressure for destructured clay erefd, Additional void ratio at initial yielding Δei, and Critical equivalent deviatoric plastic strain edcp are taken From material. For User defined enter other values or expressions.
Enter a value or expression for the Reference pressure pref, and the Initial consolidation pressure pc0.
For the Slope of critical state line you can alternatively select Match to Mohr–Coulomb criterion or Match to Matsuoka–Nakai criterion, which then matches the slope of the virgin consolidation line to the angle of internal friction. Select the Angle of internal friction ϕ as From material or User defined.
For the Initial void ratio, you can alternatively select From void ratio at reference pressure for destructured clay, which then computes the initial void ratio from the reference pressure, the void ratio at reference pressure for destructured clay, the initial consolidation pressure, and the swelling and compression indexes. Select the Void ratio at reference pressure for destructured clay erefd as From material or User defined.
See also The Modified Structured Cam-Clay Soil Model in the Structural Mechanics Theory chapter.
Extended Barcelona Basic
From the Specify list, define the elastic properties either in terms of Poisson’s ratio or Shear modulus.
The defaults for the Poisson’s ratio ν or the Shear modulus G, Density ρ, Slope of critical state line M, Swelling index κ, Swelling index for changes in suction κs, Compression index at saturation λ0, Compression index for changes in suction λs, Angle of internal friction ϕ, Weight parameter w, Soil stiffness parameter m, Plastic potential shape parameter bs, Tension to suction ratio k, Initial void ratio e0, Void ratio at reference pressure and saturation eref0, and Initial yield value for suction sy0 are taken From material. For User defined enter other values or expressions.
Enter a value or expression for the Initial suction s0, Suction s, the Reference pressure pref, and the Initial consolidation pressure pc0.
For the Slope of critical state line you can alternatively select Match to Mohr–Coulomb criterion or Match to Matsuoka–Nakai criterion, which then matches the slope of the virgin consolidation line to the angle of internal friction. Select the Angle of internal friction ϕ as From material or User defined.
For the Initial void ratio, you can alternatively select From void ratio at reference pressure and saturation, which then computes the initial void ratio from the reference pressure, the void ratio at reference pressure and saturation, the initial consolidation pressure, and the swelling and compression indexes. Then select the Void ratio at reference pressure and saturation eref0 as From material or User defined.
See also The Extended Barcelona Basic Soil Model in the Structural Mechanics Theory chapter.
Hardening Soil
Select the Failure criterionMohr–Coulomb, Matsuoka–Nakai, or Panteghini–Lagioia.
Select the Mobilized dilatancy angle Soreide, Rowe, Modified Rowe, Wehnert, Li–Dafalias, Rowe–Li–Dafalias, or User defined.
The Reference initial stiffness for primary loading Eiref, Reference stiffness for unloading and reloading Eurref, Bulk modulus in compression Kc, Poisson’s ratio ν, Density ρ, Stress exponent m, Cohesion c, Angle of internal friction ϕ, Dilatation angle ψ, Ellipse aspect ratio R, and Initial void ratio e0 are taken From material. For User defined enter other values or expressions.
Enter a value or expression for the Failure ratio Rf, the Reference pressure pref, and the Initial consolidation pressure pc0.
Select the Include dilatancy cutoff check box if needed. The defaults for the Maximum void ratio emax is taken From material. For User defined enter other value or expression. Enter a value or expression for the Initial volumetric strain εvol0.
For the Reference initial stiffness for primary loading Eiref, you can alternatively select From reference failure stiffness, which then calculates Eiref from the Reference failure stiffness E50ref.
For the Bulk modulus in compression Kc, you can alternatively select From swelling to compression ratio.
For the Ellipse aspect ratio R you can alternatively select From coefficient of earth pressure at rest. Select the Coefficient of earth pressure at rest k0nc as From material, From angle of internal friction, or User defined.
See also The Hardening Soil Model in the Structural Mechanics Theory chapter.
Hardening Soil Small Strain
Select the Failure criterionMohr–Coulomb, Matsuoka–Nakai, or Panteghini–Lagioia.
Select the Mobilized dilatancy angle Soreide, Rowe, Modified Rowe, Wehnert, Li–Dafalias, Rowe–Li–Dafalias, or User defined.
The defaults for the Reference initial stiffness for primary loading Eiref, Reference stiffness for unloading and reloading Eurref, Initial shear modulus G0, Bulk modulus in compression Kc, Poisson’s ratio ν, Density ρ, Reference shear strain γ, Stress exponent m, Cohesion c, Angle of internal friction ϕ, Dilatation angle ψ, Ellipse aspect ratio R, and Initial void ratio e0 are taken From material. For User defined enter other values or expressions.
Enter a value or expression for the Failure ratio Rf, the Reference pressure pref, and the Initial consolidation pressure pc0.
Select the Include dilatancy cutoff check box if needed. The defaults for the Maximum void ratio emax is taken From material. For User defined enter other value or expression. Enter a value or expression for the Initial volumetric strain εvol0.
Load Reversal Points
In case of cyclic loading the load reversal points are automatically detected. But in many cases the load reversal points are known a priori. The Load Reversal Points section enables to set the load reversal points manually based on a loading parameter.
For the Reference initial stiffness for primary loading Eiref, you can alternatively select From reference failure stiffness, which then calculates Eiref from the Reference failure stiffness E50ref.
For the Initial shear modulus G0, you can alternatively select From reference initial shear modulus, which then calculates G0 from the Reference initial shear modulus G0ref.
For the Bulk modulus in compression Kc, you can alternatively select From swelling to compression ratio.
For the Ellipse aspect ratio R you can alternatively select From coefficient of earth pressure at rest. Select the Coefficient of earth pressure at rest k0nc as From material, From angle of internal friction, or User defined.
See also The Hardening Soil Small Strain Model in the Structural Mechanics Theory chapter.
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.
The default is None. Select Implicit Gradient to add nonlocal regularization to the equivalent plastic strain. Enter a value for the:
Length scale, lint. The length scale should not exceed the maximum element size of the mesh.
Nonlocal coupling modulus, Hnl. This stiffness is the penalization of the difference between the local and nonlocal variables. A larger value enforces the equivalent plastic strain εpe to be closer to the nonlocal equivalent plastic strain εpe,nl.
See also Nonlocal Plasticity in the Structural Mechanics Theory chapter.
Geometric Nonlinearity
The settings in this section control the overall kinematics, the definition of the strain decomposition, and the behavior of inelastic contributions, for the material.
Select a FormulationFrom study step (default), Total Lagrangian, or Geometrically linear to set the kinematics of the deformation and the definition of strain. When From study step is selected, the study step controls the kinematics and the strain definition.
With the default From study step, a total Lagrangian formulation for large strains is used when the Include geometric nonlinearity check box is selected in the study step. If the check box is not selected, the formulation is geometrically linear, with a small strain formulation.
To have full control of the formulation, select either Total Lagrangian, or Geometrically linear. When Total Lagrangian is selected, the physics will force the Include geometric nonlinearity check box in all study steps.
When inelastic deformations are present, such as for plasticity, the elastic deformation can be obtained in two different ways: using additive decomposition of strains or using multiplicative decomposition of deformation gradients.
Select a Strain decompositionAutomatic (default), Additive, or Multiplicative to decide how the inelastic deformations are treated. This option is not available when the formulation is set to Geometrically linear.
When Automatic is selected, a multiplicative or additive decomposition is used with a total Lagrangian formulation, depending on the Include geometric nonlinearity check box status in the study step.
Select Additive to force an additive decomposition of strains.
Select Multiplicative to force a multiplicative decomposition of deformation gradients. This option is only visible if Formulation is set to Total Lagrangian.
The Strain decomposition input is only visible for material models that support both additive and multiplicative decomposition of deformation gradients.
See Lagrangian Formulation, Deformation Measures, and Inelastic Strain Contributions in the Structural Mechanics Theory chapter.
See Modeling Geometric Nonlinearity in the Structural Mechanics Modeling chapter.
See Study Settings in the COMSOL Multiphysics Reference Manual.
Energy Dissipation
Select the Calculate dissipated energy check box as needed to compute the energy dissipated by Creep, Plasticity, Viscoplasticity, or Viscoelasticity.
To display this section, click the Show More Options button () and select Advanced Physics Options in the Show More Options dialog box.
Discretization
This section is available with the Implicit gradient nonlocal plasticity model. Select the shape function for the Nonlocal equivalent plastic strain εpe,nl Automatic, Linear, Quadratic Lagrange, Quadratic serendipity, Cubic Lagrange, Cubic serendipity, Quartic Lagrange, Quartic serendipity, or Quintic Lagrange. The available options depend on the 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.
Quadrature Settings
Select the Reduced integration check box to reduce the integration points for the weak contribution of the feature. Select a method for Hourglass stabilizationAutomatic, Manual, or None to use in combination with the reduced integration scheme. The default Automatic stabilization technique is based on the shape function and shape order of the displacement field.
Control the hourglass stabilization scheme by using the Manual option. Select Shear stabilization (default) or Volumetric stabilization.
When Shear stabilization is selected, enter a stabilization shear modulus, Gstb. The value should be in the order of magnitude of the equivalent shear modulus.
When Volumetric stabilization is selected, enter a stabilization bulk modulus, Kstb. The value should be in the order of magnitude of the equivalent bulk modulus.
See also Reduced Integration and Hourglass Stabilization in the Structural Mechanics Theory chapter.
Advanced
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.
Location in User Interface
Context Menus
Ribbon
Physics tab with Solid Mechanics selected: