Damage
Use the Damage subnode to model damage and cracking in brittle materials according to various criteria. It is available with Linear Elastic Material in the Solid Mechanics and Layered Shell interfaces, and with Linear Elastic Material, Layered in the Shell interface. Phase field damage is also available for Hyperelastic Material in the Solid Mechanics interface.
The Damage subnode is only available with some COMSOL products (see https://www.comsol.com/products/specifications/). Damage is available for 3D, 2D, and 2D axisymmetry.
Shell Properties

This section is only present when Damage is used as a subnode to:
Linear Elastic Material in the Layered Shell interface. See the documentation for the Damage node in the Layered Shell chapter.
Linear Elastic Material, Layered in the Shell interface. See the documentation for the Damage node in the Shell and Plate chapter.
Damage
Select the type of Damage modelScalar damage, Mazars damage for concrete, or Phase field damage. Then follow the instructions below.
Scalar Damage
Select the type of Equivalent strainRankine, stress; Rankine, strain; Smooth Rankine, stress; Smooth Rankine, strain; Norm of elastic strain tensor; or User defined. See Strain-based Damage Models for details.
The Activate damage in compression check box is available for the Rankine, stress; Rankine, strain; Smooth Rankine, stress; Smooth Rankine, strain; or User defined equivalent strain definitions, and it is not selected by default. When selected, the damage evolution law is applied on the total undamaged stress tensor.
Select the type of Damage evolutionLinear strain softening, Exponential strain softening, Polynomial strain softening, Multilinear strain softening, or User defined. See Damage Evolution for details.
For Linear strain softening, Polynomial strain softening, Multilinear strain softening, or Exponential strain softening enter the Tensile strength σts, the default is to take the value From material. Change to User defined to enter other value or expression.
For Linear strain softening or Exponential strain softening select the type of Strain softening input and enter the Fracture energy per area Gf, the Fracture energy per volume gf, or the Strain softening parameter εf accordingly. The default is to take the value From material. Change to User defined to enter other value or expression. The available options depend on the spatial regularization method selected.
For Multilinear strain softening, enter the Shape factor λ. The default value is 0.5.
For User defined, enter an expression for the Damage evolution function d(κ).
Select the type of Spatial regularization methodNone, Crack band, or Implicit gradient. See Spatial Regularization for details.
For the Crack band method select the type of Crack band calculationElement volume/area or Element size. Note that when the Crack band method is selected, only the Fracture energy per area is available as Strain softening input.
For the Implicit gradient method enter the Internal length scale lint. If the Fracture energy per area was selected as Strain softening input, enter also the Characteristic size of the damage dissipation zone hdmg. The Implicit gradient method is available in the Solid Mechanics and Layered Shell interfaces.
Select the type of Viscous regularization methodNone or Delayed Damage.
For the Delayed damage method, enter the Characteristic time τ. The Delayed damage method is intended for time-dependent studies, and adds no contributions for other study types. See Viscous Regularization for details.
Mazars Damage for Concrete
Select the type of Equivalent strainMazars, Modified Mazars, or User defined. See Mazars Damage for Concrete for details.
Enter the Shear exponent β, the default is set to 1.06.
Select the Tensile damage evolutionLinear strain softening, Exponential strain softening, Mazars damage evolution function, or User defined.
For Linear strain softening or Exponential strain softening enter the Tensile strength σts, the default is to take the value From material. Also select the type of Tensile strain softening and enter the Fracture energy per area Gft, the Fracture energy per volume gft, or the Strain softening parameter εft accordingly. The default is to take the value From material. Change to User defined to enter other value or expression.
For Mazars damage evolution function, enter the Tensile strain threshold ε0t, and the Tensile damage evolution parameters At and Bt.
For User defined, enter an expression for the Tensile damage evolution function dt(κ).
Select the Compressive damage evolutionMazars damage evolution function or User defined.
For Mazars damage evolution function, enter the Compressive strain threshold ε0c, and the Compressive damage evolution parameters Ac and Bc.
For User defined, enter an expression for the Compressive damage evolution function dc(κ).
Select the type of Spatial regularization methodNone, Crack band, or Implicit gradient. See Spatial Regularization for details.
For the Crack band method select the type of Crack band calculationElement volume/area or Element size. Note that when the Crack band method is selected, only the Fracture energy per area is available as Strain softening input.
For the Implicit gradient method enter the Length scale lint. If the Fracture energy per area was selected as Strain softening input, enter also the Characteristic size of the damage dissipation zone hdmg. The Implicit gradient method is available in the Solid Mechanics and Layered Shell interfaces.
Select the type of Viscous regularization methodNone or Delayed Damage.
For the Delayed damage method, enter the Characteristic time τ. The Delayed damage method is intended for time-dependent studies, and adds no contributions for other study types. See Viscous Regularization for details.
Phase Field Damage
The Phase field damage model is available with the Linear Elastic Material or the Hyperelastic Material in the Solid Mechanics interface. See Phase Field Damage Models for details.
Select the type of Crack driving forceElastic strain energy density, Total strain energy density, Principal stress criterion, or User defined. See Crack Driving Force for details.
For Elastic strain energy density and Total strain energy density, enter a value for the Critical energy release rate Gc, and the Strain energy threshold Gc0. The default for the Critical energy release rate is to take the value From material. Change to User defined to enter other value or expression.
For Principal stress criterion, enter a value for the Critical fracture stress σc, and the Post-peak slope parameter ξ. The default for the Critical fracture stress is to take the value From material. Change to User defined to enter other value or expression.
For User defined, enter an expression for the dimensionless crack driving force.
For all options, enter a value for the Length scale lint.
Select the Damage evolution function to use — Power law; Cubic, Borden; or User defined. See Damage Evolution for Phase Field Damage for details.
For Power law, enter a value for the Exponent m. The default is to use a quadratic function so that m = 2.
For Cubic, Borden, enter a value for the Model parameter s. This option is available when Crack driving force is set to Elastic strain energy density or Total strain energy density.
For User defined, enter an expression for the Damage evolution function d(ϕ).
Select an option for Exclude compressive energyVolumetric only; Spectral decomposition, stress; Spectral decomposition, strain; or No split. For Hyperelastic Material, only some options are available from the list. See Strain Energy Split for details.
The Viscous regularization check box is not selected by default. When selected, a viscous term is added to the evolution of the crack phase field in time-dependent studies. Enter a value for the Characteristic time τ. See Viscous Regularization for details.
Discretization
This section is available with the Implicit gradient regularization method and the Phase field damage model. Select the shape function for the Nonlocal equivalent strain εeq, or the Crack phase field ϕ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.
Advanced
Enter the Maximum damage. The default value is 0.99999.
When a Mixed Formulation is selected in the parent Linear Elastic Material, the damage model can give spurious results. The Implicit gradient regularization method is more stable in this respect, and it is recommended when using the mixed formulation.
Modeling Damage in the Structural Mechanics Modeling chapter.
Damage Models in the Structural Mechanics Theory chapter.
Location in User Interface
Context Menus
Ribbon
Physics tab with Linear Elastic Material, Hyperelastic Material, or Linear Elastic Material, Layered, node selected in the model tree: