Viscoelasticity
Use the Viscoelasticity subnode to add viscous stress contributions to an elastic material model. This material model is available in the Solid Mechanics, Shell, Layered Shell, and Membrane interfaces, and can be used together with Linear Elastic Material, Nonlinear Elastic Material, Hyperelastic Material, Layered Linear Elastic Material, and Layered Hyperelastic Material.
Note: Some options described below are only available with certain COMSOL products (see https://www.comsol.com/products/specifications/).
See also Linear Viscoelasticity and Large Strain Viscoelasticity in the Structural Mechanics Theory chapter.
Shell Properties

This section is only present when Viscoelasticity is used as a subnode to:
Layered Linear Elastic Material or Layered Hyperelastic Elastic Material in the Shell interface. See the documentation for the Viscoelasticity node in the Shell and Plate chapter.
Layered Linear Elastic Material in the Membrane interface. See the documentation for the Viscoelasticity node in the Membrane chapter.
Sketch
Viscoelastic models are often conceptually viewed as rheological networks consisting of spring and damper elements. The sketch section shows the arrangement of springs and dampers for the chosen viscoelastic material model. If fractional derivatives are enabled, the dampers are replaced by so-called spring-pot elements. Relevant input parameters are shown next to the spring, damper, and spring-pot elements. The elastic stiffness contribution of the parent material model is indicated by a light gray spring symbol representing the bulk and shear moduli, K and G, respectively.
Thermal Effects
Viscoelastic properties have a strong dependence on the temperature. For thermorheologically simple materials, a change in the temperature can be transformed directly into a change in the time scale. Thus, the relaxation time is modified to aT(Tm, where aT(T) is a shift function.
Select a Shift function — None, Williams-Landel-Ferry, Arrhenius, Tool-Narayanaswamy-Moynihan, or User defined.
When the default, None, is kept, the shift function aT(T) is set to unity and the relaxation time is not modified.
For Williams-Landel-Ferry enter values or expressions for these properties:
-
Reference temperature Tref The default is 293.15 K.
-
WLF constant 1 C1WLF. The default is 17.44.
-
WLF constant 2 C2WLF. The default is 51.6 K.
For Arrhenius enter values or expressions for these properties:
-
Reference temperature Tref. The default is 293.15 K.
-
For Tool-Narayanaswamy-Moynihan enter values or expressions for these properties:
-
Reference temperature Tref. The default is 293.15 K.
-
-
-
-
For User defined enter a value or an expression for the shift function aT.
Viscoelasticity Model
Select a Material model Generalized Maxwell, Generalized Kelvin-Voigt, Maxwell, Kelvin-Voigt, Standard linear solid, Burgers, or User defined. Then, enter the settings for each option that follows.
From the Viscoelastic strains list select Volumetric, Deviatoric, or Volumetric and deviatoric. Select Volumetric when the viscoelastic behavior applies only to the volumetric deformation. The Deviatoric option (default) applies the viscoelastic relaxation to the shear deformation only. With Volumetric and deviatoric the viscoelastic strain is full.
For some material models, you can select the stiffness to use when solving a stationary problem. Select the Stiffness used in stationary studies Long-term or Instantaneous. With Long-term all dampers are assumed to be fully relaxed, whereas with Instantaneous all dampers are assumed to be rigid.
Generalized Maxwell
For Generalized Maxwell enter the values for the parameters that describe the viscoelastic behavior as a series of spring-dashpot pairs.
Depending on the selection done in the Viscoelastic strains list, for each Branch row enter the stiffness of the spring Km in the Bulk modulus (Pa) column and/or Gm in the Shear modulus (Pa) column, and the relaxation time constant τm in the Relaxation time (s) column for the spring-dashpot pair in branch m.
When the Use fractional derivatives check box is selected, enter the fractional order βm in the Fractional order (1) column for each spring-spring-pot branch.
For large strain viscoelasticity, in each Branch row enter the energy factor of the branch, βvm, in the Energy factor (1) column and the relaxation time constant τm in the Relaxation time (s) column for the spring-dashpot pair.
Use the Add button () to add a row to the table, the Delete button () to delete a row in the table, or the Clear Table button () to clear the whole table.
Use the Load from file button () and the Save to file button () to load and store data for the branches in a text file with space-separated columns.
When the Prune viscoelastic branches check box is selected, enter the Cutoff frequencies flower and fupper. The relaxation times τm are frozen and cannot be changed when the check box is selected. In order to change these settings again, clear the check box because pruning is only performed at the time when the check box is selected. It is also required that the relaxation times for each branch have constant values when Prune viscoelastic branches is selected.
From the Stiffness used in stationary studies list, select either Long-term or Instantaneous. With Long-term, all dampers are assumed to be relaxed; hence the branches do not contribute to the stress. The material stiffness is therefore given by the stiffness in the parent material model (for example, Linear Elastic Material, Nonlinear Elastic Material or Hyperelastic Material). With Instantaneous, all dampers are assumed to be rigid, and the material stiffness is given by springs arranged in parallel.
Generalized Kelvin–Voigt
For Generalized Kelvin-Voigt enter the values for the parameters that describe the viscoelastic behavior of multiple Kelvin–Voigt elements arranged in series.
Depending on the selection done in the Viscoelastic strains list, for each Branch row enter the stiffness of the spring Km in the column labeled Bulk modulus (Pa) and/or Gm in the column labeled Shear modulus (Pa), and the relaxation time τm in the column labeled Relaxation time (s) for the spring-dashpot pair in the element m.
When the Use fractional derivatives check box is selected, enter the fractional order βm in the Fractional order (1) column for each spring-spring-pot branch.
Use the Add button () to add a row to the table, the Delete button () to delete a row in the table, or the Clear Table button () to clear the whole table.
Use the Load from file button () and the Save to file button () to load and store data for the elements in a text file with space-separated columns.
Select the Stiffness used in stationary studies, either Long-term or Instantaneous. With Long-term, all dampers are assumed to be relaxed. The material stiffness is therefore given by springs arranged in series. With Instantaneous, all dampers are assumed to be rigid; hence the viscoelastic branches do not contribution to the strain, and the instantaneous stiffness is determined by the parent material only (for example, Linear Elastic Material, Nonlinear Elastic Material or Hyperelastic Material).
Maxwell
For Maxwell enter the parameters that describes the viscous behavior of a single dashpot connected in series with a spring.
Depending on the selection done in the Viscoelastic strains list, the relaxation time or viscosity is applied to the volumetric, deviatoric, or both volumetric and deviatoric deformation. Select an option from the Relaxation data list and edit the default as needed:
Relaxation time τv. The default is 3000 s.
Viscosity ηv of the dashpot. The default is 6·1013 Pas.
When the Use fractional derivatives check box is selected, enter the fractional order βv of the spring-pot. The default is 0.5 (dimensionless).
Note that the instantaneous stiffness is given by the parent material model (for example, Linear Elastic Material, Nonlinear Elastic Material or Hyperelastic Material).
Kelvin–Voigt
For Kelvin-Voigt enter the values for the parameter that describes the viscous behavior of the single dashpot in parallel with a spring.
Depending on the selection done in the Viscoelastic strains list, the relaxation time or viscosity is applied to the volumetric, deviatoric, or both volumetric and deviatoric deformation. Select an option from the Relaxation data list and edit the default as needed:
Relaxation time τv of the dashpot. The default is 3000 s.
Viscosity ηv of the dashpot. The default is 6·1013 Pas.
For large strain viscoelasticity, enter the Relaxation time τv. The default is 3000 s.
When the Use fractional derivatives check box is selected, enter the fractional order βv of the spring-pot. The default is 0.5 (dimensionless).
Note that the instantaneous stiffness is given by the parent material model (for example, Linear Elastic Material, Nonlinear Elastic Material or Hyperelastic Material).
Standard Linear Solid
For Standard linear solid enter the values for the parameters that describe the viscoelastic behavior of the single spring-dashpot branch.
Depending on the selection done in the Viscoelastic strains list, enter the Bulk modulus and/or the Shear modulus of the spring in the Kv and Gv fields. The default values are 20 GPa.
For linear viscoelasticity, select an option from the Relaxation data list and edit the default as needed:
Relaxation time τv of the dashpot. The default is 3000 s.
Viscosity ηv of the dashpot. The default is 6·1013 Pas.
For large strain viscoelasticity, enter the Relaxation time τv, which default value is 3000 s, and the Energy factor βv of the dashpot. The default is 0.2.
When the Use fractional derivatives check box is selected, enter the fractional order βv of the spring-pot. The default is 0.5 (dimensionless).
Note that the long-term stiffness is given by the parent material model (for example, Linear Elastic Material, Nonlinear Elastic Material or Hyperelastic Material).
Burgers
For Burgers enter the values for the parameter that describes the viscous behavior of the spring dashpot in series with a second spring-dashpot pair.
Depending on the selection done in the Viscoelastic strains list, enter the Bulk modulus and/or the Shear modulus of the second spring in the Kv2 and Gv2 fields. The default values are 20 GPa.
For linear viscoelasticity, select an option from the Relaxation data list and edit the default as needed:
Relaxation time. The default is 3000 s for both dashpots τv1 and τv2.
Viscosity. Enter the viscosity of the dashpots. The default is 6·1013 Pas for both ηv1 and ηv2.
When the Use fractional derivatives check box is selected, enter the fractional orders, βv1 and βv2, of the spring-pot pairs. The default is 0.5 (dimensionless) for each spring-pot.
Note that the instantaneous stiffness is given by the parent material model (for example, Linear Elastic Material, Nonlinear Elastic Material or Hyperelastic Material).
User Defined
When Volumetric is selected from the Viscoelastic strains list, specify the Storage and loss moduli K' and K'', the Storage and loss compliances Q' and Q'', or the Loss factor ηv that defines the complex–valued bulk modulus.
When Deviatoric is selected from the Viscoelastic strains list, specify the Storage and loss moduli G' and G'', the Storage and loss compliances J' and J'', or the Loss factor ηv that defines the complex–valued shear modulus.
When Volumetric and deviatoric is selected from the Viscoelastic strains list, specify the Storage and loss moduli K', K'', G' and G'', the Storage and loss compliances Q', Q'', J' and J'', or the Loss factor ηv that defines the complex–valued bulk and shear moduli.
These expressions can be entered as functions taken directly from interpolated data, or can be analytical expressions of the frequency variable.
The User defined viscoelastic models are applicable in Frequency Domain and Eigenfrequency study steps only.
The internal variables for the frequency f and angular frequency ω are named phys.freq and phys.omega, respectively. Here, phys is the tag of the parent physics (for instance, solid).
Discretization
To display this section, click the Show More Options button () and select Advanced Physics Options in the Show More Options dialog box.
The Use local time integration check box is selected by default. Clear it in case you want to use the global time integration scheme. The check box is only available for the Generalized Maxwell and Standard Linear Solid models. For all other viscoelasticity models, the global time integration is used.
Clear the Use local time integration check box to select the Shape function type Discontinuous Lagrange (default) or Gauss point data for the components of the auxiliary viscoelastic tensor. When the discontinuous Lagrange discretization is used, the shape function order is set as one order lower than the order used for the displacement field. This results fewer degrees of freedom being added to the model than when using Gauss point data. The accuracy does in general not differ much. If you want to enforce that the constitutive law is fulfilled at the integration points, select Gauss point data.
To compute the energy dissipation caused by viscoelasticity, enable the Calculate dissipated energy check box in the Energy Dissipation section of the parent material node.
Viscoelastic Structural Damper: Application Library path Structural_Mechanics_Module/Dynamics_and_Vibration/viscoelastic_damper_frequency
Eigenmodes of a Viscoelastic Structural Damper: Application Library path Structural_Mechanics_Module/Dynamics_and_Vibration/viscoelastic_damper_eigenmodes
Viscoelastic Structural Damper — Transient Analysis: Application Library path Structural_Mechanics_Module/Dynamics_and_Vibration/viscoelastic_damper_transient
Location in User Interface
Context Menus
Ribbon
Physics tab with Linear Elastic Material, Nonlinear Elastic Material, Hyperelastic Material, Layered Linear Elastic Material, Layered Hyperelastic Material node selected in the model tree: