- index_iThe index i of ij for the strain tensor.
C++ Type:MooseEnum
Unit:(no unit assumed)
Controllable:No
Description:The index i of ij for the strain tensor.
- index_jThe index j of ij for the strain tensor.
C++ Type:MooseEnum
Unit:(no unit assumed)
Controllable:No
Description:The index j of ij for the strain tensor.
- 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
BVStrainComponentAux
buildconstruction:Undocumented Class
The BVStrainComponentAux 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 components of the strain tensor.
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
- (test/tests/viscoelasticity/linear_kelvin.i)
- (test/tests/viscoelasticity/linear_maxwell.i)
- (test/tests/elasticity/shear.i)
- (test/tests/viscoelasticity/burger.i)
- (examples/viscoelasticity/burger/burger.i)
- (test/tests/interfaces/stick_3D.i)
- (test/tests/interfaces/stick_2D.i)
- (examples/viscoelasticity/lubby2/lubby2.i)
- (test/tests/interfaces/constant_friction_2D.i)
- (examples/viscoelasticity/linear_kelvin/linear_kelvin.i)
- (test/tests/interfaces/constant_friction_3D.i)
- (test/tests/viscoelasticity/lubby2.i)
- (test/tests/elasticity/triax.i)
(test/tests/viscoelasticity/linear_kelvin.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]
[strain_yy]
order = CONSTANT
family = MONOMIAL
[]
[stress_yy]
order = CONSTANT
family = MONOMIAL
[]
[]
[AuxKernels]
[strain_yy_aux]
type = BVStrainComponentAux
variable = strain_yy
index_i = y
index_j = y
execute_on = 'TIMESTEP_END'
[]
[stress_yy_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 = '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]
[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/linear_maxwell.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]
[strain_yy]
order = CONSTANT
family = MONOMIAL
[]
[stress_yy]
order = CONSTANT
family = MONOMIAL
[]
[]
[AuxKernels]
[strain_yy_aux]
type = BVStrainComponentAux
variable = strain_yy
index_i = y
index_j = y
execute_on = 'TIMESTEP_END'
[]
[stress_yy_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 = 'front back'
value = 0.0
[]
[v_x]
type = FunctionDirichletBC
variable = disp_x
boundary = 'right'
function = '-0.5*t'
[]
[v_y]
type = FunctionDirichletBC
variable = disp_y
boundary = 'top'
function = '0.5*t'
[]
[]
[Materials]
[elasticity]
type = BVMechanicalMaterial
displacements = 'disp_x disp_y disp_z'
bulk_modulus = 1.0
shear_modulus = 1.0
inelastic_models = 'viscoelastic'
[]
[viscoelastic]
type = BVMaxwellViscoelasticUpdate
viscosity = 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
[]
(test/tests/elasticity/shear.i)
[Mesh]
type = GeneratedMesh
dim = 3
nx = 8
ny = 8
nz = 2
xmin = 0
xmax = 10
ymin = 0
ymax = 10
zmin = 0
zmax = 2.5
[]
[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]
[strain_xy]
order = CONSTANT
family = MONOMIAL
[]
[stress_xy]
order = CONSTANT
family = MONOMIAL
[]
[]
[AuxKernels]
[strain_xy_aux]
type = BVStrainComponentAux
variable = strain_xy
index_i = x
index_j = y
execute_on = 'TIMESTEP_END'
[]
[stress_xy_aux]
type = BVStressComponentAux
variable = stress_xy
index_i = x
index_j = y
execute_on = 'TIMESTEP_END'
[]
[]
[Functions]
[disp_y_func]
type = ParsedFunction
expression = 'm*t*x'
symbol_names = 'm'
symbol_values = '-0.1'
[]
[]
[BCs]
[no_x]
type = DirichletBC
variable = disp_x
boundary = 'left right bottom top front back'
value = 0.0
[]
[no_z]
type = DirichletBC
variable = disp_z
boundary = 'left right bottom top front back'
value = 0.0
[]
[disp_y_plate]
type = FunctionDirichletBC
variable = disp_y
boundary = 'left right bottom top front back'
function = disp_y_func
[]
[]
[Materials]
[elasticity]
type = BVMechanicalMaterial
displacements = 'disp_x disp_y disp_z'
young_modulus = 10.0e+09
poisson_ratio = 0.25
[]
[]
[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 = 2.0
dt = 1.0
[]
[Outputs]
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
[]
(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/interfaces/stick_3D.i)
[Mesh]
[file]
type = FileMeshGenerator
file = mesh-3D.msh
[]
[interface_gen]
type = BVFaultInterfaceGenerator
input = file
sidesets = 'interface'
[]
[]
[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
[]
[]
[InterfaceKernels]
[traction_x]
type = BVMechanicalInterface
boundary = 'interface'
component = x
variable = disp_x
neighbor_var = disp_x
[]
[traction_y]
type = BVMechanicalInterface
boundary = 'interface'
component = y
variable = disp_y
neighbor_var = disp_y
[]
[traction_z]
type = BVMechanicalInterface
boundary = 'interface'
component = Z
variable = disp_z
neighbor_var = disp_z
[]
[]
[AuxVariables]
[strain_xy]
order = CONSTANT
family = MONOMIAL
[]
[stress_xy]
order = CONSTANT
family = MONOMIAL
[]
[]
[AuxKernels]
[strain_xy_aux]
type = BVStrainComponentAux
variable = strain_xy
index_i = x
index_j = y
execute_on = 'TIMESTEP_END'
[]
[stress_xy_aux]
type = BVStressComponentAux
variable = stress_xy
index_i = x
index_j = y
execute_on = 'TIMESTEP_END'
[]
[]
[Functions]
[disp_x_func]
type = ParsedFunction
expression = 'm*t*y'
symbol_names = 'm'
symbol_values = '0.1'
[]
[]
[BCs]
[Periodic]
[all]
variable = 'disp_x disp_y disp_z'
auto_direction = 'x'
[]
[]
[no_y]
type = DirichletBC
variable = disp_y
boundary = 'bottom top'
value = 0.0
[]
[disp_x_plate]
type = FunctionDirichletBC
variable = disp_x
boundary = 'bottom top'
function = disp_x_func
[]
[no_z]
type = DirichletBC
variable = disp_z
boundary = 'bottom top'
value = 0.0
[]
[]
[Materials]
[elasticity]
type = BVMechanicalMaterial
displacements = 'disp_x disp_y disp_z'
young_modulus = 10.0e+09
poisson_ratio = 0.25
[]
[interface]
type = BVMechanicalInterfaceMaterial
boundary = 'interface'
displacements = 'disp_x disp_y disp_z'
normal_stiffness = 1.0e+12
tangent_stiffness = 1.0e+12
[]
[]
[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 = 2.0
dt = 1.0
[]
[Outputs]
exodus = true
[]
(test/tests/interfaces/stick_2D.i)
[Mesh]
[file]
type = FileMeshGenerator
file = mesh-2D.msh
[]
[interface_gen]
type = BVFaultInterfaceGenerator
input = file
sidesets = 'interface'
[]
[]
[Variables]
[disp_x]
order = FIRST
family = LAGRANGE
[]
[disp_y]
order = FIRST
family = LAGRANGE
[]
[]
[Kernels]
[stress_x]
type = BVStressDivergence
component = x
variable = disp_x
[]
[stress_y]
type = BVStressDivergence
component = y
variable = disp_y
[]
[]
[InterfaceKernels]
[traction_x]
type = BVMechanicalInterface
boundary = 'interface'
component = x
variable = disp_x
neighbor_var = disp_x
[]
[traction_y]
type = BVMechanicalInterface
boundary = 'interface'
component = y
variable = disp_y
neighbor_var = disp_y
[]
[]
[AuxVariables]
[strain_xy]
order = CONSTANT
family = MONOMIAL
[]
[stress_xy]
order = CONSTANT
family = MONOMIAL
[]
[]
[AuxKernels]
[strain_xy_aux]
type = BVStrainComponentAux
variable = strain_xy
index_i = x
index_j = y
execute_on = 'TIMESTEP_END'
[]
[stress_xy_aux]
type = BVStressComponentAux
variable = stress_xy
index_i = x
index_j = y
execute_on = 'TIMESTEP_END'
[]
[]
[Functions]
[disp_x_func]
type = ParsedFunction
expression = 'm*t*y'
symbol_names = 'm'
symbol_values = '0.1'
[]
[]
[BCs]
[Periodic]
[all]
variable = 'disp_x disp_y'
auto_direction = 'x'
[]
[]
[no_y]
type = DirichletBC
variable = disp_y
boundary = 'bottom top'
value = 0.0
[]
[disp_x_plate]
type = FunctionDirichletBC
variable = disp_x
boundary = 'bottom top'
function = disp_x_func
[]
[]
[Materials]
[elasticity]
type = BVMechanicalMaterial
displacements = 'disp_x disp_y'
young_modulus = 10.0e+09
poisson_ratio = 0.25
[]
[interface]
type = BVMechanicalInterfaceMaterial
boundary = 'interface'
displacements = 'disp_x disp_y'
normal_stiffness = 1.0e+12
tangent_stiffness = 1.0e+12
[]
[]
[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 = 2.0
dt = 1.0
[]
[Outputs]
exodus = true
[]
(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/interfaces/constant_friction_2D.i)
[Mesh]
[file]
type = FileMeshGenerator
file = mesh-2D.msh
[]
[fault_gen]
type = BVFaultInterfaceGenerator
input = file
sidesets = 'interface'
[]
[]
[Variables]
[disp_x]
order = FIRST
family = LAGRANGE
[]
[disp_y]
order = FIRST
family = LAGRANGE
[]
[]
[Kernels]
[stress_x]
type = BVStressDivergence
component = x
variable = disp_x
[]
[stress_y]
type = BVStressDivergence
component = y
variable = disp_y
[]
[]
[InterfaceKernels]
[traction_x]
type = BVMechanicalInterface
boundary = 'interface'
component = x
variable = disp_x
neighbor_var = disp_x
[]
[traction_y]
type = BVMechanicalInterface
boundary = 'interface'
component = y
variable = disp_y
neighbor_var = disp_y
[]
[]
[AuxVariables]
[strain_xy]
order = CONSTANT
family = MONOMIAL
[]
[stress_xy]
order = CONSTANT
family = MONOMIAL
[]
[stress_yy]
order = CONSTANT
family = MONOMIAL
[]
[normal_stress]
order = CONSTANT
family = MONOMIAL
[]
[shear_stress]
order = CONSTANT
family = MONOMIAL
[]
[]
[AuxKernels]
[strain_xy_aux]
type = BVStrainComponentAux
variable = strain_xy
index_i = x
index_j = y
execute_on = 'TIMESTEP_END'
[]
[stress_xy_aux]
type = BVStressComponentAux
variable = stress_xy
index_i = x
index_j = y
execute_on = 'TIMESTEP_END'
[]
[stress_yy_aux]
type = BVStressComponentAux
variable = stress_yy
index_i = y
index_j = y
execute_on = 'TIMESTEP_END'
[]
[normal_stress_aux]
type = BVFaultNormalStressAux
variable = normal_stress
boundary = 'interface'
execute_on = 'TIMESTEP_END'
[]
[shear_stress_aux]
type = BVFaultShearStressAux
variable = shear_stress
boundary = 'interface'
execute_on = 'TIMESTEP_END'
[]
[]
[Functions]
[disp_x_func]
type = ParsedFunction
expression = 'm*t*y'
symbol_names = 'm'
symbol_values = '0.1'
[]
[]
[BCs]
[Periodic]
[left_right]
auto_direction = x
primary = 'left'
secondary = 'right'
variable = 'disp_x disp_y'
[]
[]
[no_y]
type = DirichletBC
variable = disp_y
boundary = 'bottom top'
value = 0.0
[]
[disp_x_plate]
type = FunctionDirichletBC
variable = disp_x
boundary = 'bottom top'
function = disp_x_func
[]
[]
[Materials]
[elasticity]
type = BVMechanicalMaterial
displacements = 'disp_x disp_y'
young_modulus = 10.0
poisson_ratio = 0.25
initial_stress = '0.0 1.6 0.0'
[]
[interface]
type = BVMechanicalInterfaceMaterial
boundary = 'interface'
displacements = 'disp_x disp_y'
normal_stiffness = 1.0e+03
tangent_stiffness = 1.0e+03
friction_model = 'constant_friction'
[]
[constant_friction]
type = BVConstantFrictionUpdate
friction = 0.5
[]
[]
# [Preconditioning]
# [hypre]
# type = SMP
# full = true
# petsc_options_iname = '-pc_type -pc_hypre_type'
# petsc_options_value = 'hypre boomeramg'
# []
# []
[Preconditioning]
[asm]
type = SMP
petsc_options = '-snes_ksp_ew'
petsc_options_iname = '-ksp_type -ksp_rtol -ksp_max_it
-pc_type
-sub_pc_type
-snes_type -snes_atol -snes_rtol -snes_max_it -snes_linesearch_type
-ksp_gmres_restart'
petsc_options_value = 'fgmres 1e-10 100
asm
ilu
newtonls 1e-12 1e-08 100 basic
201'
[]
[]
[Executioner]
type = Transient
solve_type = 'NEWTON'
start_time = 0.0
end_time = 4.0
dt = 1.0
[]
[Outputs]
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
[]
[]
(test/tests/interfaces/constant_friction_3D.i)
[Mesh]
[file]
type = FileMeshGenerator
file = mesh-3D.msh
[]
[fault_gen]
type = BVFaultInterfaceGenerator
input = file
sidesets = 'interface'
[]
[]
[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
[]
[]
[InterfaceKernels]
[traction_x]
type = BVMechanicalInterface
boundary = 'interface'
component = x
variable = disp_x
neighbor_var = disp_x
[]
[traction_y]
type = BVMechanicalInterface
boundary = 'interface'
component = y
variable = disp_y
neighbor_var = disp_y
[]
[traction_z]
type = BVMechanicalInterface
boundary = 'interface'
component = Z
variable = disp_z
neighbor_var = disp_z
[]
[]
[AuxVariables]
[strain_xy]
order = CONSTANT
family = MONOMIAL
[]
[stress_xy]
order = CONSTANT
family = MONOMIAL
[]
[stress_yy]
order = CONSTANT
family = MONOMIAL
[]
[normal_stress]
order = CONSTANT
family = MONOMIAL
[]
[shear_stress]
order = CONSTANT
family = MONOMIAL
[]
[]
[AuxKernels]
[strain_xy_aux]
type = BVStrainComponentAux
variable = strain_xy
index_i = x
index_j = y
execute_on = 'TIMESTEP_END'
[]
[stress_xy_aux]
type = BVStressComponentAux
variable = stress_xy
index_i = x
index_j = y
execute_on = 'TIMESTEP_END'
[]
[stress_yy_aux]
type = BVStressComponentAux
variable = stress_yy
index_i = y
index_j = y
execute_on = 'TIMESTEP_END'
[]
[normal_stress_aux]
type = BVFaultNormalStressAux
variable = normal_stress
boundary = 'interface'
execute_on = 'TIMESTEP_END'
[]
[shear_stress_aux]
type = BVFaultShearStressAux
variable = shear_stress
boundary = 'interface'
execute_on = 'TIMESTEP_END'
[]
[]
[Functions]
[disp_x_func]
type = ParsedFunction
expression = 'm*t*y'
symbol_names = 'm'
symbol_values = '0.1'
[]
[]
[BCs]
[Periodic]
[all]
variable = 'disp_x disp_y disp_z'
auto_direction = 'x z'
[]
[]
[no_y]
type = DirichletBC
variable = disp_y
boundary = 'bottom top'
value = 0.0
[]
[disp_x_plate]
type = FunctionDirichletBC
variable = disp_x
boundary = 'bottom top'
function = disp_x_func
[]
[no_z]
type = DirichletBC
variable = disp_z
boundary = 'bottom top'
value = 0.0
[]
[]
[Materials]
[elasticity]
type = BVMechanicalMaterial
displacements = 'disp_x disp_y disp_z'
young_modulus = 10.0
poisson_ratio = 0.25
initial_stress = '0.0 1.6 0.0'
[]
[interface]
type = BVMechanicalInterfaceMaterial
boundary = 'interface'
displacements = 'disp_x disp_y disp_z'
normal_stiffness = 1.0e+03
tangent_stiffness = 1.0e+03
friction_model = 'constant_friction'
[]
[constant_friction]
type = BVConstantFrictionUpdate
friction = 0.5
[]
[]
# [Preconditioning]
# [hypre]
# type = SMP
# full = true
# petsc_options_iname = '-pc_type -pc_hypre_type'
# petsc_options_value = 'hypre boomeramg'
# []
# []
[Preconditioning]
[asm]
type = SMP
petsc_options = '-snes_ksp_ew'
petsc_options_iname = '-ksp_type -ksp_rtol -ksp_max_it
-pc_type
-sub_pc_type
-snes_type -snes_atol -snes_rtol -snes_max_it -snes_linesearch_type
-ksp_gmres_restart'
petsc_options_value = 'fgmres 1e-10 100
asm
ilu
newtonls 1e-12 1e-08 100 basic
201'
[]
[]
[Executioner]
type = Transient
solve_type = 'NEWTON'
start_time = 0.0
end_time = 4.0
dt = 1.0
[]
[Outputs]
exodus = true
[]
(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/elasticity/triax.i)
[Mesh]
[gen]
type = GeneratedMeshGenerator
dim = 3
nx = 4
ny = 4
nz = 8
xmin = -0.25
xmax = 0.25
ymin = -0.25
ymax = 0.25
zmin = 0
zmax = 1
[]
[no_ux]
type = ExtraNodesetGenerator
input = gen
coord = '0.0 -0.25 0.0; 0.0 0.25 0.0'
new_boundary = 'no_disp_x'
[]
[no_uy]
type = ExtraNodesetGenerator
input = no_ux
coord = '-0.25 0.0 0.0; 0.25 0.0 0.0'
new_boundary = 'no_disp_y'
[]
[]
[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]
[strain_zz]
order = CONSTANT
family = MONOMIAL
[]
[stress_zz]
order = CONSTANT
family = MONOMIAL
[]
[pressure]
order = CONSTANT
family = MONOMIAL
[]
[dev_stress]
order = CONSTANT
family = MONOMIAL
[]
[]
[AuxKernels]
[strain_zz_aux]
type = BVStrainComponentAux
variable = strain_zz
index_i = z
index_j = z
execute_on = 'TIMESTEP_END'
[]
[stress_zz_aux]
type = BVStressComponentAux
variable = stress_zz
index_i = z
index_j = z
execute_on = 'TIMESTEP_END'
[]
[pressure_aux]
type = BVPressureAux
variable = pressure
execute_on = 'INITIAL TIMESTEP_END'
[]
[dev_stress_aux]
type = BVMisesStressAux
variable = dev_stress
execute_on = 'INITIAL TIMESTEP_END'
[]
[]
[BCs]
[BVPressure]
[pressure_conf]
boundary = 'left right bottom top front'
displacement_vars = 'disp_x disp_y'
value = 10.0e+06
[]
[]
[no_x]
type = DirichletBC
variable = disp_x
boundary = 'no_disp_x'
value = 0.0
[]
[no_y]
type = DirichletBC
variable = disp_y
boundary = 'no_disp_y'
value = 0.0
[]
[no_z]
type = DirichletBC
variable = disp_z
boundary = 'back'
value = 0.0
[]
[axial_load]
type = FunctionDirichletBC
variable = disp_z
boundary = 'front'
function = '-1.0e-02*t'
[]
[]
[Materials]
[elasticity]
type = BVMechanicalMaterial
displacements = 'disp_x disp_y disp_z'
young_modulus = 10.0e+09
poisson_ratio = 0.25
initial_stress = '-10.0e+06 -10.0e+06 -10.0e+06'
[]
[]
[Preconditioning]
[hypre]
type = SMP
full = true
petsc_options_iname = '-pc_type -pc_hypre_type -snes_atol'
petsc_options_value = 'hypre boomeramg 1.0e-12'
[]
[]
[Executioner]
type = Transient
automatic_scaling = true
solve_type = 'NEWTON'
start_time = 0.0
end_time = 2.0
dt = 1.0
[]
[Outputs]
exodus = true
[]