- variableThe name of the variable that this object applies to
C++ Type:AuxVariableName
Unit:(no unit assumed)
Controllable:No
Description:The name of the variable that this object applies to
BVEqvStrainRateAux
buildconstruction:Undocumented Class
The BVEqvStrainRateAux has not been documented. The content listed below should be used as a starting point for documenting the class, which includes the typical automatic documentation associated with a MooseObject; however, what is contained is ultimately determined by what is necessary to make the documentation clear for users.
Class for outputting the equivalent strain rate.
Overview
Example Input File Syntax
Input Parameters
- blockThe list of blocks (ids or names) that this object will be applied
C++ Type:std::vector<SubdomainName>
Unit:(no unit assumed)
Controllable:No
Description:The list of blocks (ids or names) that this object will be applied
- boundaryThe list of boundaries (ids or names) from the mesh where this object applies
C++ Type:std::vector<BoundaryName>
Unit:(no unit assumed)
Controllable:No
Description:The list of boundaries (ids or names) from the mesh where this object applies
- check_boundary_restrictedTrueWhether to check for multiple element sides on the boundary in the case of a boundary restricted, element aux variable. Setting this to false will allow contribution to a single element's elemental value(s) from multiple boundary sides on the same element (example: when the restricted boundary exists on two or more sides of an element, such as at a corner of a mesh
Default:True
C++ Type:bool
Unit:(no unit assumed)
Controllable:No
Description:Whether to check for multiple element sides on the boundary in the case of a boundary restricted, element aux variable. Setting this to false will allow contribution to a single element's elemental value(s) from multiple boundary sides on the same element (example: when the restricted boundary exists on two or more sides of an element, such as at a corner of a mesh
- execute_onLINEAR TIMESTEP_ENDThe list of flag(s) indicating when this object should be executed. For a description of each flag, see https://mooseframework.inl.gov/source/interfaces/SetupInterface.html.
Default:LINEAR TIMESTEP_END
C++ Type:ExecFlagEnum
Unit:(no unit assumed)
Options:NONE, INITIAL, LINEAR, NONLINEAR_CONVERGENCE, NONLINEAR, POSTCHECK, TIMESTEP_END, TIMESTEP_BEGIN, MULTIAPP_FIXED_POINT_END, MULTIAPP_FIXED_POINT_BEGIN, FINAL, CUSTOM, PRE_DISPLACE
Controllable:No
Description:The list of flag(s) indicating when this object should be executed. For a description of each flag, see https://mooseframework.inl.gov/source/interfaces/SetupInterface.html.
- prop_getter_suffixAn optional suffix parameter that can be appended to any attempt to retrieve/get material properties. The suffix will be prepended with a '_' character.
C++ Type:MaterialPropertyName
Unit:(no unit assumed)
Controllable:No
Description:An optional suffix parameter that can be appended to any attempt to retrieve/get material properties. The suffix will be prepended with a '_' character.
- use_interpolated_stateFalseFor the old and older state use projected material properties interpolated at the quadrature points. To set up projection use the ProjectedStatefulMaterialStorageAction.
Default:False
C++ Type:bool
Unit:(no unit assumed)
Controllable:No
Description:For the old and older state use projected material properties interpolated at the quadrature points. To set up projection use the ProjectedStatefulMaterialStorageAction.
Optional Parameters
- control_tagsAdds user-defined labels for accessing object parameters via control logic.
C++ Type:std::vector<std::string>
Unit:(no unit assumed)
Controllable:No
Description:Adds user-defined labels for accessing object parameters via control logic.
- enableTrueSet the enabled status of the MooseObject.
Default:True
C++ Type:bool
Unit:(no unit assumed)
Controllable:Yes
Description:Set the enabled status of the MooseObject.
- seed0The seed for the master random number generator
Default:0
C++ Type:unsigned int
Unit:(no unit assumed)
Controllable:No
Description:The seed for the master random number generator
- use_displaced_meshFalseWhether or not this object should use the displaced mesh for computation. Note that in the case this is true but no displacements are provided in the Mesh block the undisplaced mesh will still be used.
Default:False
C++ Type:bool
Unit:(no unit assumed)
Controllable:No
Description:Whether or not this object should use the displaced mesh for computation. Note that in the case this is true but no displacements are provided in the Mesh block the undisplaced mesh will still be used.
Advanced Parameters
Input Files
- (examples/viscoelasticity/blanco-martin/blanco-martin-RTL-2024.i)
- (examples/viscoelasticity/lubby2/lubby2.i)
- (test/tests/viscoelasticity/lubby2.i)
- (test/tests/viscoelasticity/blanco-martin-lemaitre.i)
- (examples/viscoelasticity/linear_kelvin/linear_kelvin.i)
- (examples/viscoelasticity/blanco-martin/blanco-martin-lemaitre-2024.i)
- (examples/viscoelasticity/burger/burger.i)
- (test/tests/viscoelasticity/blanco-martin-rtl.i)
- (test/tests/viscoelasticity/burger.i)
(examples/viscoelasticity/blanco-martin/blanco-martin-RTL-2024.i)
# Modified Lemaitre creep model
# See Blanco-Martin et al. (2023)
# Parameters
# Units: stress in MPa, time in days, strain in m / m
E = 12000
nu = 0.3
alpha = 0.575
kr1 = 1.302
beta1 = 3.053
kr2 = 0.091
beta2 = 1.053
A = 100
n = 9
A1 = 0.034
n1 = 1.499
P = 5.0
Q = 5.0
[Mesh]
type = GeneratedMesh
dim = 3
# nx = 5
# ny = 10
# nz = 5
nx = 1
ny = 1
nz = 1
xmin = 0
xmax = 65e-03
ymin = 0
ymax = 130e-03
zmin = 0
zmax = 65e-03
[]
[Variables]
[disp_x]
order = FIRST
family = LAGRANGE
[]
[disp_y]
order = FIRST
family = LAGRANGE
[]
[disp_z]
order = FIRST
family = LAGRANGE
[]
[]
[Kernels]
[stress_x]
type = BVStressDivergence
component = x
variable = disp_x
[]
[stress_y]
type = BVStressDivergence
component = y
variable = disp_y
[]
[stress_z]
type = BVStressDivergence
component = z
variable = disp_z
[]
[]
[AuxVariables]
[eqv_stress]
order = CONSTANT
family = MONOMIAL
[]
[eqv_strain]
order = CONSTANT
family = MONOMIAL
[]
[eqv_strain_rate]
order = CONSTANT
family = MONOMIAL
[]
[eqv_creep_strain_L]
order = CONSTANT
family = MONOMIAL
[]
[strain_yy]
order = CONSTANT
family = MONOMIAL
[]
[stress_yy]
order = CONSTANT
family = MONOMIAL
[]
[]
[AuxKernels]
[eqv_stress_aux]
type = BVMisesStressAux
variable = eqv_stress
execute_on = 'TIMESTEP_END'
[]
[eqv_strain_aux]
type = BVEqvStrainAux
variable = eqv_strain
execute_on = 'TIMESTEP_END'
[]
[eqv_strain_rate_aux]
type = BVEqvStrainRateAux
variable = eqv_strain_rate
execute_on = 'TIMESTEP_END'
[]
[eqv_creep_strain_L_aux]
type = ADMaterialRealAux
variable = eqv_creep_strain_L
property = eqv_creep_strain_L
execute_on = 'TIMESTEP_END'
[]
[strain_zz_aux]
type = BVStrainComponentAux
variable = strain_yy
index_i = y
index_j = y
execute_on = 'TIMESTEP_END'
[]
[stress_zz_aux]
type = BVStressComponentAux
variable = stress_yy
index_i = y
index_j = y
execute_on = 'TIMESTEP_END'
[]
[]
[Functions]
[loading2]
type = ParsedFunction
expression = 'if(t<=20,5.5,if(t<=55,6,if(t<=76,7,if(t<=97,10,if(t<=118,15,if(t<=139,20,25))))))'
[]
[]
[BCs]
[no_x]
type = DirichletBC
variable = disp_x
boundary = 'left'
value = 0.0
[]
[no_y]
type = DirichletBC
variable = disp_y
boundary = 'bottom'
value = 0.0
[]
[no_z]
type = DirichletBC
variable = disp_z
boundary = 'back'
value = 0.0
[]
[BVPressure]
[pressure_right]
boundary = 'right'
displacement_vars = 'disp_x disp_y disp_z'
value = ${P}
[]
[pressure_front]
boundary = 'front'
displacement_vars = 'disp_x disp_y disp_z'
value = ${P}
[]
[pressure_top]
boundary = 'top'
displacement_vars = 'disp_x disp_y disp_z'
function = loading2
[]
[]
[]
[Materials]
[elasticity]
type = BVMechanicalMaterial
displacements = 'disp_x disp_y disp_z'
young_modulus = ${E}
poisson_ratio = ${nu}
initial_stress = '-${P} -${Q} -${P}'
inelastic_models = 'viscoelastic'
[]
[viscoelastic]
type = BVBlancoMartinModelUpdate
alpha = ${alpha}
kr1 = ${kr1}
beta1 = ${beta1}
kr2 = ${kr2}
beta2 = ${beta2}
A1 = ${A1}
n1 = ${n1}
A = ${A}
n = ${n}
B = 0.0
m = ${n}
[]
[]
[Postprocessors]
[eqv_strain_rate]
type = ElementAverageValue
variable = eqv_strain_rate
outputs = csv
[]
[eqv_strain]
type = ElementAverageValue
variable = eqv_strain
outputs = csv
[]
[eqv_strain_L]
type = ElementAverageValue
variable = eqv_creep_strain_L
outputs = csv
[]
[strain_zz]
type = ElementAverageValue
variable = strain_yy
outputs = csv
[]
[stress_zz]
type = ElementAverageValue
variable = stress_yy
outputs = csv
[]
[q]
type = ElementAverageValue
variable = eqv_stress
outputs = csv
[]
[]
[Preconditioning]
active = 'hypre'
[hypre]
type = SMP
full = true
petsc_options = '-snes_ksp_ew'
petsc_options_iname = '-pc_type -pc_hypre_type
-snes_atol -snes_rtol -snes_stol -snes_max_it
-snes_linesearch_type'
petsc_options_value = 'hypre boomeramg
1.0e-10 1.0e-12 0 20
basic'
[]
[superlu]
type = SMP
full = true
petsc_options = '-snes_ksp_ew -snes_converged_reason -ksp_converged_reason -ksp_gmres_modifiedgramschmidt -ksp_diagonal_scale -ksp_diagonal_scale_fix'
petsc_options_iname = '-snes_type
-snes_atol -snes_rtol -snes_max_it
-pc_type -pc_factor_mat_solver_package
-snes_linesearch_type'
petsc_options_value = 'newtonls
1e-10 1e-12 50
lu superlu_dist
l2'
[]
[asm]
type = SMP
petsc_options = '-snes_ksp_ew'
petsc_options_iname = '-ksp_type
-pc_type
-sub_pc_type
-snes_type -snes_atol -snes_rtol -snes_max_it -snes_linesearch_type
-ksp_gmres_restart'
petsc_options_value = 'fgmres
asm
ilu
newtonls 1e-10 1e-10 120 basic
201'
[]
[]
[Executioner]
type = Transient
solve_type = 'NEWTON'
automatic_scaling = true
start_time = 0.0
end_time = 200 # 200 days
dt = 0.02
timestep_tolerance = 1.0e-10
[]
[Outputs]
perf_graph = true
exodus = true
[csv]
type = CSV
execute_on = 'timestep_end'
[]
[](examples/viscoelasticity/lubby2/lubby2.i)
[Mesh]
type = GeneratedMesh
dim = 3
nx = 20
ny = 20
nz = 2
xmin = 0
xmax = 1
ymin = 0
ymax = 1
zmin = 0
zmax = 0.1
[]
[Variables]
[disp_x]
order = FIRST
family = LAGRANGE
[]
[disp_y]
order = FIRST
family = LAGRANGE
[]
[disp_z]
order = FIRST
family = LAGRANGE
[]
[]
[Kernels]
[stress_x]
type = BVStressDivergence
component = x
variable = disp_x
[]
[stress_y]
type = BVStressDivergence
component = y
variable = disp_y
[]
[stress_z]
type = BVStressDivergence
component = z
variable = disp_z
[]
[]
[AuxVariables]
[eqv_stress]
order = CONSTANT
family = MONOMIAL
[]
[eqv_strain]
order = CONSTANT
family = MONOMIAL
[]
[eqv_strain_rate]
order = CONSTANT
family = MONOMIAL
[]
[Kelvin_creep_strain]
order = CONSTANT
family = MONOMIAL
[]
[strain_yy]
order = CONSTANT
family = MONOMIAL
[]
[]
[AuxKernels]
[eqv_stress_aux]
type = BVMisesStressAux
variable = eqv_stress
execute_on = 'TIMESTEP_END'
[]
[eqv_strain_aux]
type = BVEqvStrainAux
variable = eqv_strain
execute_on = 'TIMESTEP_END'
[]
[eqv_strain_rate_aux]
type = BVEqvStrainRateAux
variable = eqv_strain_rate
execute_on = 'TIMESTEP_END'
[]
[Kelvin_creep_strain_aux]
type = ADMaterialRealAux
variable = Kelvin_creep_strain
property = eqv_Kelvin_creep_strain
execute_on = 'TIMESTEP_END'
[]
[strain_yy_aux]
type = BVStrainComponentAux
variable = strain_yy
index_i = y
index_j = y
execute_on = 'TIMESTEP_END'
[]
[]
[BCs]
[no_x]
type = DirichletBC
variable = disp_x
boundary = 'left'
value = 0.0
[]
[no_y]
type = DirichletBC
variable = disp_y
boundary = 'bottom'
value = 0.0
[]
[no_z]
type = DirichletBC
variable = disp_z
boundary = 'front back'
value = 0.0
[]
[BVPressure]
[pressure_right]
boundary = 'right'
displacement_vars = 'disp_x disp_y disp_z'
value = 1.0
[]
[pressure_top]
boundary = 'top'
displacement_vars = 'disp_x disp_y disp_z'
value = -1.0
[]
[]
[]
[Materials]
[elasticity]
type = BVMechanicalMaterial
displacements = 'disp_x disp_y disp_z'
bulk_modulus = 1.0
shear_modulus = 2.0
initial_stress = '-1.0 1.0 0.0'
inelastic_models = 'viscoelastic'
[]
[viscoelastic]
type = BVLubby2ModelUpdate
viscosity_maxwell = 10.0
viscosity_kelvin = 1.0
shear_modulus_kelvin = 1.0
sigma_0 = 1.0
m_1 = 0.327
m_2 = 0.267
m_G = 0.254
[]
[]
[Preconditioning]
active = 'hypre'
[hypre]
type = SMP
full = true
petsc_options = '-snes_ksp_ew -snes_converged_reason -ksp_converged_reason -ksp_diagonal_scale -ksp_diagonal_scale_fix'
petsc_options_iname = '-pc_type -pc_hypre_type
-snes_atol -snes_rtol -snes_stol -snes_max_it -snes_linesearch_type'
petsc_options_value = 'hypre boomeramg
1.0e-10 1.0e-12 0 20
basic'
[]
[]
[Postprocessors]
[e]
type = ElementAverageValue
variable = eqv_strain
outputs = csv
[]
[e_rate]
type = ElementAverageValue
variable = eqv_strain_rate
outputs = csv
[]
[Kelvin_creep]
type = ElementAverageValue
variable = Kelvin_creep_strain
outputs = csv
[]
[]
[Executioner]
type = Transient
solve_type = 'NEWTON'
start_time = 0.0
end_time = 10.0
dt = 0.05
[]
[Outputs]
perf_graph = true
execute_on = 'TIMESTEP_END'
print_linear_residuals = false
exodus = true
[csv]
type = CSV
[]
[]
(test/tests/viscoelasticity/lubby2.i)
[Mesh]
type = GeneratedMesh
dim = 3
nx = 8
ny = 8
nz = 2
xmin = 0
xmax = 1
ymin = 0
ymax = 1
zmin = 0
zmax = 0.1
[]
[Variables]
[disp_x]
order = FIRST
family = LAGRANGE
[]
[disp_y]
order = FIRST
family = LAGRANGE
[]
[disp_z]
order = FIRST
family = LAGRANGE
[]
[]
[Kernels]
[stress_x]
type = BVStressDivergence
component = x
variable = disp_x
[]
[stress_y]
type = BVStressDivergence
component = y
variable = disp_y
[]
[stress_z]
type = BVStressDivergence
component = z
variable = disp_z
[]
[]
[AuxVariables]
[eqv_stress]
order = CONSTANT
family = MONOMIAL
[]
[eqv_strain]
order = CONSTANT
family = MONOMIAL
[]
[eqv_strain_rate]
order = CONSTANT
family = MONOMIAL
[]
[Kelvin_creep_strain]
order = CONSTANT
family = MONOMIAL
[]
[strain_yy]
order = CONSTANT
family = MONOMIAL
[]
[]
[AuxKernels]
[eqv_stress_aux]
type = BVMisesStressAux
variable = eqv_stress
execute_on = 'TIMESTEP_END'
[]
[eqv_strain_aux]
type = BVEqvStrainAux
variable = eqv_strain
execute_on = 'TIMESTEP_END'
[]
[eqv_strain_rate_aux]
type = BVEqvStrainRateAux
variable = eqv_strain_rate
execute_on = 'TIMESTEP_END'
[]
[Kelvin_creep_strain_aux]
type = ADMaterialRealAux
variable = Kelvin_creep_strain
property = eqv_Kelvin_creep_strain
execute_on = 'TIMESTEP_END'
[]
[strain_yy_aux]
type = BVStrainComponentAux
variable = strain_yy
index_i = y
index_j = y
execute_on = 'TIMESTEP_END'
[]
[]
[BCs]
[no_x]
type = DirichletBC
variable = disp_x
boundary = 'left'
value = 0.0
[]
[no_y]
type = DirichletBC
variable = disp_y
boundary = 'bottom'
value = 0.0
[]
[no_z]
type = DirichletBC
variable = disp_z
boundary = 'front back'
value = 0.0
[]
[BVPressure]
[pressure_right]
boundary = 'right'
displacement_vars = 'disp_x disp_y disp_z'
value = 1.0
[]
[pressure_top]
boundary = 'top'
displacement_vars = 'disp_x disp_y disp_z'
value = -1.0
[]
[]
[]
[Materials]
[elasticity]
type = BVMechanicalMaterial
displacements = 'disp_x disp_y disp_z'
bulk_modulus = 1.0
shear_modulus = 1.0
initial_stress = '-1.0 1.0 0.0'
inelastic_models = 'viscoelastic'
[]
[viscoelastic]
type = BVLubby2ModelUpdate
viscosity_maxwell = 10.0
viscosity_kelvin = 1.0
shear_modulus_kelvin = 1.0
m_1 = 0.327
m_2 = 0.267
m_G = 0.254
[]
[]
[Preconditioning]
[hypre]
type = SMP
full = true
petsc_options_iname = '-pc_type -pc_hypre_type'
petsc_options_value = 'hypre boomeramg'
[]
[]
[Executioner]
type = Transient
solve_type = 'NEWTON'
start_time = 0.0
end_time = 1.0
dt = 0.1
[]
[Outputs]
exodus = true
[]
(test/tests/viscoelasticity/blanco-martin-lemaitre.i)
# Modified Lemaitre creep model
# See Blanco-Martin et al. (2024)
# Parameters
# Units: stress in MPa, time in days, strain in m / m
E = 12000
nu = 0.3
alpha = 0.326
kr1 = 0.7
beta1 = 2.922
kr2 = 0.009
beta2 = 0.867
P = 5.0
Q = 5.0
Q1 = 5.5
[Mesh]
type = GeneratedMesh
dim = 3
# nx = 5
# ny = 10
# nz = 5
nx = 1
ny = 1
nz = 1
xmin = 0
xmax = 65e-03
ymin = 0
ymax = 130e-03
zmin = 0
zmax = 65e-03
[]
[Variables]
[disp_x]
order = FIRST
family = LAGRANGE
[]
[disp_y]
order = FIRST
family = LAGRANGE
[]
[disp_z]
order = FIRST
family = LAGRANGE
[]
[]
[Kernels]
[stress_x]
type = BVStressDivergence
component = x
variable = disp_x
[]
[stress_y]
type = BVStressDivergence
component = y
variable = disp_y
[]
[stress_z]
type = BVStressDivergence
component = z
variable = disp_z
[]
[]
[AuxVariables]
[eqv_stress]
order = CONSTANT
family = MONOMIAL
[]
[eqv_strain]
order = CONSTANT
family = MONOMIAL
[]
[eqv_strain_rate]
order = CONSTANT
family = MONOMIAL
[]
[eqv_creep_strain_L]
order = CONSTANT
family = MONOMIAL
[]
[strain_yy]
order = CONSTANT
family = MONOMIAL
[]
[stress_yy]
order = CONSTANT
family = MONOMIAL
[]
[]
[AuxKernels]
[eqv_stress_aux]
type = BVMisesStressAux
variable = eqv_stress
execute_on = 'TIMESTEP_END'
[]
[eqv_strain_aux]
type = BVEqvStrainAux
variable = eqv_strain
execute_on = 'TIMESTEP_END'
[]
[eqv_strain_rate_aux]
type = BVEqvStrainRateAux
variable = eqv_strain_rate
execute_on = 'TIMESTEP_END'
[]
[eqv_creep_strain_L_aux]
type = ADMaterialRealAux
variable = eqv_creep_strain_L
property = eqv_creep_strain
execute_on = 'TIMESTEP_END'
[]
[strain_zz_aux]
type = BVStrainComponentAux
variable = strain_yy
index_i = y
index_j = y
execute_on = 'TIMESTEP_END'
[]
[stress_zz_aux]
type = BVStressComponentAux
variable = stress_yy
index_i = y
index_j = y
execute_on = 'TIMESTEP_END'
[]
[]
[BCs]
[no_x]
type = DirichletBC
variable = disp_x
boundary = 'left'
value = 0.0
[]
[no_y]
type = DirichletBC
variable = disp_y
boundary = 'bottom'
value = 0.0
[]
[no_z]
type = DirichletBC
variable = disp_z
boundary = 'back'
value = 0.0
[]
[BVPressure]
[pressure_right]
boundary = 'right'
displacement_vars = 'disp_x disp_y disp_z'
value = ${P}
[]
[pressure_front]
boundary = 'front'
displacement_vars = 'disp_x disp_y disp_z'
value = ${P}
[]
[pressure_top]
boundary = 'top'
displacement_vars = 'disp_x disp_y disp_z'
function = ${Q1}
[]
[]
[]
[Materials]
[elasticity]
type = BVMechanicalMaterial
displacements = 'disp_x disp_y disp_z'
young_modulus = ${E}
poisson_ratio = ${nu}
initial_stress = '-${P} -${Q} -${P}'
inelastic_models = 'viscoelastic'
[]
# [viscoelastic]
# type = BVBlancoMartinModelUpdate
# alpha = ${alpha}
# kr1 = ${kr1}
# beta1 = ${beta1}
# kr2 = ${kr2}
# beta2 = ${beta2}
# A1 = ${A1}
# n1 = ${n1}
# A = ${A}
# n = ${n}
# B = 0.0
# m = ${n}
# []
[viscoelastic]
type = BVModifiedLemaitreModelUpdate
alpha = ${alpha}
kr1 = ${kr1}
beta1 = ${beta1}
kr2 = ${kr2}
beta2 = ${beta2}
[]
[]
[Preconditioning]
active = 'hypre'
[hypre]
type = SMP
full = true
petsc_options = '-snes_ksp_ew'
petsc_options_iname = '-pc_type -pc_hypre_type
-snes_atol -snes_rtol -snes_stol -snes_max_it
-snes_linesearch_type'
petsc_options_value = 'hypre boomeramg
1.0e-10 1.0e-12 0 20
basic'
[]
[superlu]
type = SMP
full = true
petsc_options = '-snes_ksp_ew -snes_converged_reason -ksp_converged_reason -ksp_gmres_modifiedgramschmidt -ksp_diagonal_scale -ksp_diagonal_scale_fix'
petsc_options_iname = '-snes_type
-snes_atol -snes_rtol -snes_max_it
-pc_type -pc_factor_mat_solver_package
-snes_linesearch_type'
petsc_options_value = 'newtonls
1e-10 1e-12 50
lu superlu_dist
l2'
[]
[asm]
type = SMP
petsc_options = '-snes_ksp_ew'
petsc_options_iname = '-ksp_type
-pc_type
-sub_pc_type
-snes_type -snes_atol -snes_rtol -snes_max_it -snes_linesearch_type
-ksp_gmres_restart'
petsc_options_value = 'fgmres
asm
ilu
newtonls 1e-10 1e-10 120 basic
201'
[]
[]
[Executioner]
type = Transient
solve_type = 'NEWTON'
automatic_scaling = true
start_time = 0.0
end_time = 1
dt = 0.02
timestep_tolerance = 1.0e-10
[]
[Outputs]
perf_graph = true
exodus = true
[](examples/viscoelasticity/linear_kelvin/linear_kelvin.i)
[Mesh]
type = GeneratedMesh
dim = 3
nx = 20
ny = 20
nz = 2
xmin = 0
xmax = 1
ymin = 0
ymax = 1
zmin = 0
zmax = 0.1
[]
[Variables]
[disp_x]
order = FIRST
family = LAGRANGE
[]
[disp_y]
order = FIRST
family = LAGRANGE
[]
[disp_z]
order = FIRST
family = LAGRANGE
[]
[]
[Kernels]
[stress_x]
type = BVStressDivergence
component = x
variable = disp_x
[]
[stress_y]
type = BVStressDivergence
component = y
variable = disp_y
[]
[stress_z]
type = BVStressDivergence
component = z
variable = disp_z
[]
[]
[AuxVariables]
[eqv_stress]
order = CONSTANT
family = MONOMIAL
[]
[eqv_strain]
order = CONSTANT
family = MONOMIAL
[]
[eqv_strain_rate]
order = CONSTANT
family = MONOMIAL
[]
[creep_strain]
order = CONSTANT
family = MONOMIAL
[]
[strain_yy]
order = CONSTANT
family = MONOMIAL
[]
[]
[AuxKernels]
[eqv_stress_aux]
type = BVMisesStressAux
variable = eqv_stress
execute_on = 'TIMESTEP_END'
[]
[eqv_strain_aux]
type = BVEqvStrainAux
variable = eqv_strain
execute_on = 'TIMESTEP_END'
[]
[eqv_strain_rate_aux]
type = BVEqvStrainRateAux
variable = eqv_strain_rate
execute_on = 'TIMESTEP_END'
[]
[creep_strain_aux]
type = ADMaterialRealAux
variable = creep_strain
property = eqv_creep_strain
execute_on = 'TIMESTEP_END'
[]
[strain_yy_aux]
type = BVStrainComponentAux
variable = strain_yy
index_i = y
index_j = y
execute_on = 'TIMESTEP_END'
[]
[]
[BCs]
[no_x]
type = DirichletBC
variable = disp_x
boundary = 'left'
value = 0.0
[]
[no_y]
type = DirichletBC
variable = disp_y
boundary = 'bottom'
value = 0.0
[]
[no_z]
type = DirichletBC
variable = disp_z
boundary = 'front back'
value = 0.0
[]
[BVPressure]
[pressure_right]
boundary = 'right'
displacement_vars = 'disp_x disp_y disp_z'
value = 1.0
[]
[pressure_top]
boundary = 'top'
displacement_vars = 'disp_x disp_y disp_z'
value = -1.0
[]
[]
[]
[Materials]
[elasticity]
type = BVMechanicalMaterial
displacements = 'disp_x disp_y disp_z'
bulk_modulus = 1.0
shear_modulus = 1.0
initial_stress = '-1.0 1.0 0.0'
inelastic_models = 'viscoelastic'
[]
[viscoelastic]
type = BVKelvinViscoelasticUpdate
viscosity = 1.0
shear_modulus = 1.0
[]
[]
[Preconditioning]
active = 'hypre'
[hypre]
type = SMP
full = true
petsc_options = '-snes_ksp_ew -snes_converged_reason -ksp_converged_reason -ksp_diagonal_scale -ksp_diagonal_scale_fix'
petsc_options_iname = '-pc_type -pc_hypre_type
-snes_atol -snes_rtol -snes_stol -snes_max_it -snes_linesearch_type'
petsc_options_value = 'hypre boomeramg
1.0e-10 1.0e-12 0 20
basic'
[]
[]
[Postprocessors]
[e]
type = ElementAverageValue
variable = eqv_strain
outputs = csv
[]
[e_rate]
type = ElementAverageValue
variable = eqv_strain_rate
outputs = csv
[]
[e_creep]
type = ElementAverageValue
variable = creep_strain
outputs = csv
[]
[]
[Executioner]
type = Transient
solve_type = 'NEWTON'
start_time = 0.0
end_time = 10.0
dt = 0.05
[]
[Outputs]
perf_graph = true
execute_on = 'TIMESTEP_END'
print_linear_residuals = false
exodus = true
[csv]
type = CSV
[]
[]
(examples/viscoelasticity/blanco-martin/blanco-martin-lemaitre-2024.i)
# Modified Lemaitre creep model
# See Blanco-Martin et al. (2024)
# Parameters
# Units: stress in MPa, time in days, strain in m / m
E = 12000
nu = 0.3
alpha = 0.326
kr1 = 0.7
beta1 = 2.922
kr2 = 0.009
beta2 = 0.867
P = 5.0
Q = 5.0
[Mesh]
type = GeneratedMesh
dim = 3
# nx = 5
# ny = 10
# nz = 5
nx = 1
ny = 1
nz = 1
xmin = 0
xmax = 65e-03
ymin = 0
ymax = 130e-03
zmin = 0
zmax = 65e-03
[]
[Variables]
[disp_x]
order = FIRST
family = LAGRANGE
[]
[disp_y]
order = FIRST
family = LAGRANGE
[]
[disp_z]
order = FIRST
family = LAGRANGE
[]
[]
[Kernels]
[stress_x]
type = BVStressDivergence
component = x
variable = disp_x
[]
[stress_y]
type = BVStressDivergence
component = y
variable = disp_y
[]
[stress_z]
type = BVStressDivergence
component = z
variable = disp_z
[]
[]
[AuxVariables]
[eqv_stress]
order = CONSTANT
family = MONOMIAL
[]
[eqv_strain]
order = CONSTANT
family = MONOMIAL
[]
[eqv_strain_rate]
order = CONSTANT
family = MONOMIAL
[]
[eqv_creep_strain_L]
order = CONSTANT
family = MONOMIAL
[]
[strain_yy]
order = CONSTANT
family = MONOMIAL
[]
[stress_yy]
order = CONSTANT
family = MONOMIAL
[]
[]
[AuxKernels]
[eqv_stress_aux]
type = BVMisesStressAux
variable = eqv_stress
execute_on = 'TIMESTEP_END'
[]
[eqv_strain_aux]
type = BVEqvStrainAux
variable = eqv_strain
execute_on = 'TIMESTEP_END'
[]
[eqv_strain_rate_aux]
type = BVEqvStrainRateAux
variable = eqv_strain_rate
execute_on = 'TIMESTEP_END'
[]
[eqv_creep_strain_L_aux]
type = ADMaterialRealAux
variable = eqv_creep_strain_L
property = eqv_creep_strain
execute_on = 'TIMESTEP_END'
[]
[strain_zz_aux]
type = BVStrainComponentAux
variable = strain_yy
index_i = y
index_j = y
execute_on = 'TIMESTEP_END'
[]
[stress_zz_aux]
type = BVStressComponentAux
variable = stress_yy
index_i = y
index_j = y
execute_on = 'TIMESTEP_END'
[]
[]
[Functions]
[loading2]
type = ParsedFunction
expression = 'if(t<=20,5.5,if(t<=55,6,if(t<=76,7,if(t<=97,10,if(t<=118,15,if(t<=139,20,25))))))'
[]
[]
[BCs]
[no_x]
type = DirichletBC
variable = disp_x
boundary = 'left'
value = 0.0
[]
[no_y]
type = DirichletBC
variable = disp_y
boundary = 'bottom'
value = 0.0
[]
[no_z]
type = DirichletBC
variable = disp_z
boundary = 'back'
value = 0.0
[]
[BVPressure]
[pressure_right]
boundary = 'right'
displacement_vars = 'disp_x disp_y disp_z'
value = ${P}
[]
[pressure_front]
boundary = 'front'
displacement_vars = 'disp_x disp_y disp_z'
value = ${P}
[]
[pressure_top]
boundary = 'top'
displacement_vars = 'disp_x disp_y disp_z'
function = loading2
[]
[]
[]
[Materials]
[elasticity]
type = BVMechanicalMaterial
displacements = 'disp_x disp_y disp_z'
young_modulus = ${E}
poisson_ratio = ${nu}
initial_stress = '-${P} -${Q} -${P}'
inelastic_models = 'viscoelastic'
[]
# [viscoelastic]
# type = BVBlancoMartinModelUpdate
# alpha = ${alpha}
# kr1 = ${kr1}
# beta1 = ${beta1}
# kr2 = ${kr2}
# beta2 = ${beta2}
# A1 = ${A1}
# n1 = ${n1}
# A = ${A}
# n = ${n}
# B = 0.0
# m = ${n}
# []
[viscoelastic]
type = BVModifiedLemaitreModelUpdate
alpha = ${alpha}
kr1 = ${kr1}
beta1 = ${beta1}
kr2 = ${kr2}
beta2 = ${beta2}
[]
[]
[Postprocessors]
[eqv_strain_rate]
type = ElementAverageValue
variable = eqv_strain_rate
outputs = csv
[]
[eqv_strain]
type = ElementAverageValue
variable = eqv_strain
outputs = csv
[]
[eqv_strain_L]
type = ElementAverageValue
variable = eqv_creep_strain_L
outputs = csv
[]
[strain_zz]
type = ElementAverageValue
variable = strain_yy
outputs = csv
[]
[stress_zz]
type = ElementAverageValue
variable = stress_yy
outputs = csv
[]
[q]
type = ElementAverageValue
variable = eqv_stress
outputs = csv
[]
[]
[Preconditioning]
active = 'hypre'
[hypre]
type = SMP
full = true
petsc_options = '-snes_ksp_ew'
petsc_options_iname = '-pc_type -pc_hypre_type
-snes_atol -snes_rtol -snes_stol -snes_max_it
-snes_linesearch_type'
petsc_options_value = 'hypre boomeramg
1.0e-10 1.0e-12 0 20
basic'
[]
[superlu]
type = SMP
full = true
petsc_options = '-snes_ksp_ew -snes_converged_reason -ksp_converged_reason -ksp_gmres_modifiedgramschmidt -ksp_diagonal_scale -ksp_diagonal_scale_fix'
petsc_options_iname = '-snes_type
-snes_atol -snes_rtol -snes_max_it
-pc_type -pc_factor_mat_solver_package
-snes_linesearch_type'
petsc_options_value = 'newtonls
1e-10 1e-12 50
lu superlu_dist
l2'
[]
[asm]
type = SMP
petsc_options = '-snes_ksp_ew'
petsc_options_iname = '-ksp_type
-pc_type
-sub_pc_type
-snes_type -snes_atol -snes_rtol -snes_max_it -snes_linesearch_type
-ksp_gmres_restart'
petsc_options_value = 'fgmres
asm
ilu
newtonls 1e-10 1e-10 120 basic
201'
[]
[]
[Executioner]
type = Transient
solve_type = 'NEWTON'
automatic_scaling = true
start_time = 0.0
end_time = 200 # 200 days
dt = 0.02
timestep_tolerance = 1.0e-10
[]
[Outputs]
perf_graph = true
exodus = true
[csv]
type = CSV
execute_on = 'timestep_end'
[]
[](examples/viscoelasticity/burger/burger.i)
[Mesh]
type = GeneratedMesh
dim = 3
nx = 20
ny = 20
nz = 2
xmin = 0
xmax = 1
ymin = 0
ymax = 1
zmin = 0
zmax = 0.1
[]
[Variables]
[disp_x]
order = FIRST
family = LAGRANGE
[]
[disp_y]
order = FIRST
family = LAGRANGE
[]
[disp_z]
order = FIRST
family = LAGRANGE
[]
[]
[Kernels]
[stress_x]
type = BVStressDivergence
component = x
variable = disp_x
[]
[stress_y]
type = BVStressDivergence
component = y
variable = disp_y
[]
[stress_z]
type = BVStressDivergence
component = z
variable = disp_z
[]
[]
[AuxVariables]
[eqv_stress]
order = CONSTANT
family = MONOMIAL
[]
[eqv_strain]
order = CONSTANT
family = MONOMIAL
[]
[eqv_strain_rate]
order = CONSTANT
family = MONOMIAL
[]
[Kelvin_creep_strain]
order = CONSTANT
family = MONOMIAL
[]
[strain_yy]
order = CONSTANT
family = MONOMIAL
[]
[]
[AuxKernels]
[eqv_stress_aux]
type = BVMisesStressAux
variable = eqv_stress
execute_on = 'TIMESTEP_END'
[]
[eqv_strain_aux]
type = BVEqvStrainAux
variable = eqv_strain
execute_on = 'TIMESTEP_END'
[]
[eqv_strain_rate_aux]
type = BVEqvStrainRateAux
variable = eqv_strain_rate
execute_on = 'TIMESTEP_END'
[]
[Kelvin_creep_strain_aux]
type = ADMaterialRealAux
variable = Kelvin_creep_strain
property = eqv_Kelvin_creep_strain
execute_on = 'TIMESTEP_END'
[]
[strain_yy_aux]
type = BVStrainComponentAux
variable = strain_yy
index_i = y
index_j = y
execute_on = 'TIMESTEP_END'
[]
[]
[BCs]
[no_x]
type = DirichletBC
variable = disp_x
boundary = 'left'
value = 0.0
[]
[no_y]
type = DirichletBC
variable = disp_y
boundary = 'bottom'
value = 0.0
[]
[no_z]
type = DirichletBC
variable = disp_z
boundary = 'front back'
value = 0.0
[]
[BVPressure]
[pressure_right]
boundary = 'right'
displacement_vars = 'disp_x disp_y disp_z'
value = 1.0
[]
[pressure_top]
boundary = 'top'
displacement_vars = 'disp_x disp_y disp_z'
value = -1.0
[]
[]
[]
[Materials]
[elasticity]
type = BVMechanicalMaterial
displacements = 'disp_x disp_y disp_z'
bulk_modulus = 1.0
shear_modulus = 1.0
initial_stress = '-1.0 1.0 0.0'
inelastic_models = 'viscoelastic'
[]
[viscoelastic]
type = BVBurgerModelUpdate
viscosity_maxwell = 10.0
viscosity_kelvin = 1.0
shear_modulus_kelvin = 1.0
[]
[]
[Preconditioning]
active = 'hypre'
[hypre]
type = SMP
full = true
petsc_options = '-snes_ksp_ew -snes_converged_reason -ksp_converged_reason -ksp_diagonal_scale -ksp_diagonal_scale_fix'
petsc_options_iname = '-pc_type -pc_hypre_type
-snes_atol -snes_rtol -snes_stol -snes_max_it -snes_linesearch_type'
petsc_options_value = 'hypre boomeramg
1.0e-10 1.0e-12 0 20
basic'
[]
[]
[Postprocessors]
[e]
type = ElementAverageValue
variable = eqv_strain
outputs = csv
[]
[e_rate]
type = ElementAverageValue
variable = eqv_strain_rate
outputs = csv
[]
[Kelvin_creep]
type = ElementAverageValue
variable = Kelvin_creep_strain
outputs = csv
[]
[]
[Executioner]
type = Transient
solve_type = 'NEWTON'
start_time = 0.0
end_time = 10.0
dt = 0.05
[]
[Outputs]
perf_graph = true
execute_on = 'TIMESTEP_END'
print_linear_residuals = false
exodus = true
[csv]
type = CSV
[]
[]
(test/tests/viscoelasticity/blanco-martin-rtl.i)
# Modified Lemaitre creep model
# See Blanco-Martin et al. (2023)
# Parameters
# Units: stress in MPa, time in days, strain in m / m
E = 12000
nu = 0.3
alpha = 0.575
kr1 = 1.302
beta1 = 3.053
kr2 = 0.091
beta2 = 1.053
A = 100
n = 9
A1 = 0.034
n1 = 1.499
P = 5.0
Q = 5.0
Q1 = 5.5
[Mesh]
type = GeneratedMesh
dim = 3
# nx = 5
# ny = 10
# nz = 5
nx = 1
ny = 1
nz = 1
xmin = 0
xmax = 65e-03
ymin = 0
ymax = 130e-03
zmin = 0
zmax = 65e-03
[]
[Variables]
[disp_x]
order = FIRST
family = LAGRANGE
[]
[disp_y]
order = FIRST
family = LAGRANGE
[]
[disp_z]
order = FIRST
family = LAGRANGE
[]
[]
[Kernels]
[stress_x]
type = BVStressDivergence
component = x
variable = disp_x
[]
[stress_y]
type = BVStressDivergence
component = y
variable = disp_y
[]
[stress_z]
type = BVStressDivergence
component = z
variable = disp_z
[]
[]
[AuxVariables]
[eqv_stress]
order = CONSTANT
family = MONOMIAL
[]
[eqv_strain]
order = CONSTANT
family = MONOMIAL
[]
[eqv_strain_rate]
order = CONSTANT
family = MONOMIAL
[]
[eqv_creep_strain_L]
order = CONSTANT
family = MONOMIAL
[]
[strain_yy]
order = CONSTANT
family = MONOMIAL
[]
[stress_yy]
order = CONSTANT
family = MONOMIAL
[]
[]
[AuxKernels]
[eqv_stress_aux]
type = BVMisesStressAux
variable = eqv_stress
execute_on = 'TIMESTEP_END'
[]
[eqv_strain_aux]
type = BVEqvStrainAux
variable = eqv_strain
execute_on = 'TIMESTEP_END'
[]
[eqv_strain_rate_aux]
type = BVEqvStrainRateAux
variable = eqv_strain_rate
execute_on = 'TIMESTEP_END'
[]
[eqv_creep_strain_L_aux]
type = ADMaterialRealAux
variable = eqv_creep_strain_L
property = eqv_creep_strain_L
execute_on = 'TIMESTEP_END'
[]
[strain_zz_aux]
type = BVStrainComponentAux
variable = strain_yy
index_i = y
index_j = y
execute_on = 'TIMESTEP_END'
[]
[stress_zz_aux]
type = BVStressComponentAux
variable = stress_yy
index_i = y
index_j = y
execute_on = 'TIMESTEP_END'
[]
[]
[BCs]
[no_x]
type = DirichletBC
variable = disp_x
boundary = 'left'
value = 0.0
[]
[no_y]
type = DirichletBC
variable = disp_y
boundary = 'bottom'
value = 0.0
[]
[no_z]
type = DirichletBC
variable = disp_z
boundary = 'back'
value = 0.0
[]
[BVPressure]
[pressure_right]
boundary = 'right'
displacement_vars = 'disp_x disp_y disp_z'
value = ${P}
[]
[pressure_front]
boundary = 'front'
displacement_vars = 'disp_x disp_y disp_z'
value = ${P}
[]
[pressure_top]
boundary = 'top'
displacement_vars = 'disp_x disp_y disp_z'
function = ${Q1}
[]
[]
[]
[Materials]
[elasticity]
type = BVMechanicalMaterial
displacements = 'disp_x disp_y disp_z'
young_modulus = ${E}
poisson_ratio = ${nu}
initial_stress = '-${P} -${Q} -${P}'
inelastic_models = 'viscoelastic'
[]
[viscoelastic]
type = BVBlancoMartinModelUpdate
alpha = ${alpha}
kr1 = ${kr1}
beta1 = ${beta1}
kr2 = ${kr2}
beta2 = ${beta2}
A1 = ${A1}
n1 = ${n1}
A = ${A}
n = ${n}
B = 0.0
m = ${n}
[]
[]
[Preconditioning]
active = 'hypre'
[hypre]
type = SMP
full = true
petsc_options = '-snes_ksp_ew'
petsc_options_iname = '-pc_type -pc_hypre_type
-snes_atol -snes_rtol -snes_stol -snes_max_it
-snes_linesearch_type'
petsc_options_value = 'hypre boomeramg
1.0e-10 1.0e-12 0 20
basic'
[]
[superlu]
type = SMP
full = true
petsc_options = '-snes_ksp_ew -snes_converged_reason -ksp_converged_reason -ksp_gmres_modifiedgramschmidt -ksp_diagonal_scale -ksp_diagonal_scale_fix'
petsc_options_iname = '-snes_type
-snes_atol -snes_rtol -snes_max_it
-pc_type -pc_factor_mat_solver_package
-snes_linesearch_type'
petsc_options_value = 'newtonls
1e-10 1e-12 50
lu superlu_dist
l2'
[]
[asm]
type = SMP
petsc_options = '-snes_ksp_ew'
petsc_options_iname = '-ksp_type
-pc_type
-sub_pc_type
-snes_type -snes_atol -snes_rtol -snes_max_it -snes_linesearch_type
-ksp_gmres_restart'
petsc_options_value = 'fgmres
asm
ilu
newtonls 1e-10 1e-10 120 basic
201'
[]
[]
[Executioner]
type = Transient
solve_type = 'NEWTON'
automatic_scaling = true
start_time = 0.0
end_time = 1 #
dt = 0.02
timestep_tolerance = 1.0e-10
[]
[Outputs]
perf_graph = true
exodus = true
[](test/tests/viscoelasticity/burger.i)
[Mesh]
type = GeneratedMesh
dim = 3
nx = 8
ny = 8
nz = 2
xmin = 0
xmax = 1
ymin = 0
ymax = 1
zmin = 0
zmax = 0.1
[]
[Variables]
[disp_x]
order = FIRST
family = LAGRANGE
[]
[disp_y]
order = FIRST
family = LAGRANGE
[]
[disp_z]
order = FIRST
family = LAGRANGE
[]
[]
[Kernels]
[stress_x]
type = BVStressDivergence
component = x
variable = disp_x
[]
[stress_y]
type = BVStressDivergence
component = y
variable = disp_y
[]
[stress_z]
type = BVStressDivergence
component = z
variable = disp_z
[]
[]
[AuxVariables]
[eqv_stress]
order = CONSTANT
family = MONOMIAL
[]
[eqv_strain]
order = CONSTANT
family = MONOMIAL
[]
[eqv_strain_rate]
order = CONSTANT
family = MONOMIAL
[]
[Kelvin_creep_strain]
order = CONSTANT
family = MONOMIAL
[]
[strain_yy]
order = CONSTANT
family = MONOMIAL
[]
[]
[AuxKernels]
[eqv_stress_aux]
type = BVMisesStressAux
variable = eqv_stress
execute_on = 'TIMESTEP_END'
[]
[eqv_strain_aux]
type = BVEqvStrainAux
variable = eqv_strain
execute_on = 'TIMESTEP_END'
[]
[eqv_strain_rate_aux]
type = BVEqvStrainRateAux
variable = eqv_strain_rate
execute_on = 'TIMESTEP_END'
[]
[Kelvin_creep_strain_aux]
type = ADMaterialRealAux
variable = Kelvin_creep_strain
property = eqv_Kelvin_creep_strain
execute_on = 'TIMESTEP_END'
[]
[strain_yy_aux]
type = BVStrainComponentAux
variable = strain_yy
index_i = y
index_j = y
execute_on = 'TIMESTEP_END'
[]
[]
[BCs]
[no_x]
type = DirichletBC
variable = disp_x
boundary = 'left'
value = 0.0
[]
[no_y]
type = DirichletBC
variable = disp_y
boundary = 'bottom'
value = 0.0
[]
[no_z]
type = DirichletBC
variable = disp_z
boundary = 'front back'
value = 0.0
[]
[BVPressure]
[pressure_right]
boundary = 'right'
displacement_vars = 'disp_x disp_y disp_z'
value = 1.0
[]
[pressure_top]
boundary = 'top'
displacement_vars = 'disp_x disp_y disp_z'
value = -1.0
[]
[]
[]
[Materials]
[elasticity]
type = BVMechanicalMaterial
displacements = 'disp_x disp_y disp_z'
bulk_modulus = 1.0
shear_modulus = 1.0
initial_stress = '-1.0 1.0 0.0'
inelastic_models = 'viscoelastic'
[]
[viscoelastic]
type = BVBurgerModelUpdate
viscosity_maxwell = 10.0
viscosity_kelvin = 1.0
shear_modulus_kelvin = 1.0
[]
[]
[Preconditioning]
[hypre]
type = SMP
full = true
petsc_options_iname = '-pc_type -pc_hypre_type'
petsc_options_value = 'hypre boomeramg'
[]
[]
[Executioner]
type = Transient
solve_type = 'NEWTON'
start_time = 0.0
end_time = 1.0
dt = 0.1
[]
[Outputs]
exodus = true
[]