Material Tensor Average

Computes the average of a RankTwoTensor component over a volume.

Description

The MaterialTensorAverage postprocessor computes the volume average of the Rank-2 tensor component specified by the user. $(1)var element = document.getElementById("moose-equation-d64a58d9-41e4-48e5-9f8e-f6ce1c510422");katex.render(" S^{avg} = \\int_V T_{ij} dV / \\int_V dv", element, {displayMode:true,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ where $var element = document.getElementById("moose-equation-cd813d25-d749-4881-b9e1-e758edf8d33f");katex.render("S^{avg}", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ is the computed volume average quantity and $var element = document.getElementById("moose-equation-db4610c8-b849-4b2e-955e-99837d2545c7");katex.render("T_{ij}", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ is the tensor component selected by the user. The tensor component indices, $var element = document.getElementById("moose-equation-499a56bb-d44a-400d-97f7-cd9ba3cb261d");katex.render("i", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ and $var element = document.getElementById("moose-equation-fbcc77fa-3109-4748-8d34-830de5731d77");katex.render("j", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ , range from 0 to 2 as shown in the reference tensor $(2)var element = document.getElementById("moose-equation-92548c01-ff0a-43a9-9a19-b0ae7755dc8a");katex.render(" T_{ij} = \\begin{bmatrix} 0,0 & 0,1 & 0,2 \\\\ 1,0 & 1,1 & 1,2 \\\\ 2,0 & 2,1 & 2,2 \\end{bmatrix}", element, {displayMode:true,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$

The parameter use_displaced_mesh controls the volume utilized to compute the average. If use_displaced_mesh=true the average is compute utilizing the deformed volume, if use_displaced_mesh=false (default) the average is compute utilizing the initial volume.

Example Input File

[Postprocessors]
  [./szz_avg]
    type = MaterialTensorAverage
    rank_two_tensor = stress
    index_i = 2
    index_j = 2
    use_displaced_mesh = true
  []
  [./szz_int]
    type = MaterialTensorIntegral
    rank_two_tensor = stress
    index_i = 2
    index_j = 2
    use_displaced_mesh = true
  []
  [./szz_avg_aux]
    type = ElementAverageValue
    variable = stress_zz
    use_displaced_mesh = true
  []
  [./szz_int_aux]
    type = ElementIntegralVariablePostprocessor
    variable = stress_zz
    use_displaced_mesh = true
  []
[]

Input Parameters

index_iThe index i of ij for the tensor to output (0, 1, 2)
C++ Type:unsigned int
Controllable:No
Description:The index i of ij for the tensor to output (0, 1, 2)
index_jThe index j of ij for the tensor to output (0, 1, 2)
C++ Type:unsigned int
Controllable:No
Description:The index j of ij for the tensor to output (0, 1, 2)
rank_two_tensorThe rank two material property tensor name
C++ Type:MaterialPropertyName
Controllable:No
Description:The rank two material property tensor name

Required Parameters

blockThe list of blocks (ids or names) that this object will be applied
C++ Type:std::vector<SubdomainName>
Controllable:No
Description:The list of blocks (ids or names) that this object will be applied
execute_onTIMESTEP_ENDThe list of flag(s) indicating when this object should be executed, the available options include FORWARD, ADJOINT, HOMOGENEOUS_FORWARD, ADJOINT_TIMESTEP_BEGIN, ADJOINT_TIMESTEP_END, NONE, INITIAL, LINEAR, NONLINEAR, POSTCHECK, TIMESTEP_END, TIMESTEP_BEGIN, MULTIAPP_FIXED_POINT_END, MULTIAPP_FIXED_POINT_BEGIN, FINAL, CUSTOM.
Default:TIMESTEP_END
C++ Type:ExecFlagEnum
Options:FORWARD, ADJOINT, HOMOGENEOUS_FORWARD, ADJOINT_TIMESTEP_BEGIN, ADJOINT_TIMESTEP_END, NONE, INITIAL, LINEAR, NONLINEAR, POSTCHECK, TIMESTEP_END, TIMESTEP_BEGIN, MULTIAPP_FIXED_POINT_END, MULTIAPP_FIXED_POINT_BEGIN, FINAL, CUSTOM, TRANSFER
Controllable:No
Description:The list of flag(s) indicating when this object should be executed, the available options include FORWARD, ADJOINT, HOMOGENEOUS_FORWARD, ADJOINT_TIMESTEP_BEGIN, ADJOINT_TIMESTEP_END, NONE, INITIAL, LINEAR, NONLINEAR, POSTCHECK, TIMESTEP_END, TIMESTEP_BEGIN, MULTIAPP_FIXED_POINT_END, MULTIAPP_FIXED_POINT_BEGIN, FINAL, CUSTOM.
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
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
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

allow_duplicate_execution_on_initialFalseIn the case where this UserObject is depended upon by an initial condition, allow it to be executed twice during the initial setup (once before the IC and again after mesh adaptivity (if applicable).
Default:False
C++ Type:bool
Controllable:No
Description:In the case where this UserObject is depended upon by an initial condition, allow it to be executed twice during the initial setup (once before the IC and again after mesh adaptivity (if applicable).
control_tagsAdds user-defined labels for accessing object parameters via control logic.
C++ Type:std::vector<std::string>
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
Controllable:Yes
Description:Set the enabled status of the MooseObject.
execution_order_group0Execution order groups are executed in increasing order (e.g., the lowest number is executed first). Note that negative group numbers may be used to execute groups before the default (0) group. Please refer to the user object documentation for ordering of user object execution within a group.
Default:0
C++ Type:int
Controllable:No
Description:Execution order groups are executed in increasing order (e.g., the lowest number is executed first). Note that negative group numbers may be used to execute groups before the default (0) group. Please refer to the user object documentation for ordering of user object execution within a group.
force_postauxFalseForces the UserObject to be executed in POSTAUX
Default:False
C++ Type:bool
Controllable:No
Description:Forces the UserObject to be executed in POSTAUX
force_preauxFalseForces the UserObject to be executed in PREAUX
Default:False
C++ Type:bool
Controllable:No
Description:Forces the UserObject to be executed in PREAUX
force_preicFalseForces the UserObject to be executed in PREIC during initial setup
Default:False
C++ Type:bool
Controllable:No
Description:Forces the UserObject to be executed in PREIC during initial setup
implicitTrueDetermines whether this object is calculated using an implicit or explicit form
Default:True
C++ Type:bool
Controllable:No
Description:Determines whether this object is calculated using an implicit or explicit form
outputsVector of output names where you would like to restrict the output of variables(s) associated with this object
C++ Type:std::vector<OutputName>
Controllable:No
Description:Vector of output names where you would like to restrict the output of variables(s) associated with this object
seed0The seed for the master random number generator
Default:0
C++ Type:unsigned int
Controllable:No
Description:The seed for the master random number generator
use_displaced_meshTrueWhether 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:True
C++ Type:bool
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

(modules/solid_mechanics/test/tests/postprocessors/material_tensor_average_test.i)

(modules/solid_mechanics/test/tests/postprocessors/material_tensor_average_test.i)

[Mesh]
  [./msh]
    type = GeneratedMeshGenerator
    dim = 3
    xmax = 2
    ymax = 2
    zmax = 2
  []
[]

[GlobalParams]
  displacements = 'disp_x disp_y disp_z'
[]

[Physics]
  [SolidMechanics]
    [QuasiStatic]
      [./all]
        strain = FINITE
        add_variables = true
        generate_output = 'stress_zz'
      [../]
    [../]
  [../]
[]

[BCs]
  [./back_z]
    type = DirichletBC
    variable = disp_z
    boundary = back
    value = 0.0
  [../]
  [./left_x]
    type = DirichletBC
    variable = disp_x
    boundary = left
    value = 0.0
  [../]
  [./bottom_y]
    type = DirichletBC
    variable = disp_y
    boundary = bottom
    value = 0.0
  [../]
  [./move_front]
    type = FunctionDirichletBC
    variable = disp_z
    boundary = front
    function = 't/10.'
  [../]
[]

[Materials]
  [./stress]
    type = ComputeFiniteStrainElasticStress
  [../]
  [./elasticity_tensor]
    type = ComputeElasticityTensor
    fill_method = symmetric9
    C_ijkl = '1.5e6 0.75e6 0.75e6 1.5e6 0.75e6 1.5e6 0.375e6 0.375e6 0.375e6'
  [../]
[]

[Postprocessors]
  [./szz_avg]
    type =MaterialTensorAverage
    rank_two_tensor = stress
    index_i = 2
    index_j = 2
    use_displaced_mesh = true
  []
  [./szz_int]
    type =MaterialTensorIntegral
    rank_two_tensor = stress
    index_i = 2
    index_j = 2
    use_displaced_mesh = true
  []
  [./szz_avg_aux]
    type =ElementAverageValue
    variable = stress_zz
    use_displaced_mesh = true
  []
  [./szz_int_aux]
    type =ElementIntegralVariablePostprocessor
    variable = stress_zz
    use_displaced_mesh = true
  []
[]


[Preconditioning]
  [./SMP]
    type = SMP
    full = true
  [../]
[]

[Executioner]
  type = Transient
  solve_type = NEWTON
  petsc_options_iname = '-pc_type '
  petsc_options_value = lu
  nl_rel_tol = 1e-10
  nl_abs_tol = 1e-6
  l_max_its = 20
  start_time = 0.0
  dt = 0.2
  end_time = 1.0
[]

[Outputs]
  csv = true
[]

(modules/solid_mechanics/test/tests/postprocessors/material_tensor_average_test.i)

[Mesh]
  [./msh]
    type = GeneratedMeshGenerator
    dim = 3
    xmax = 2
    ymax = 2
    zmax = 2
  []
[]

[GlobalParams]
  displacements = 'disp_x disp_y disp_z'
[]

[Physics]
  [SolidMechanics]
    [QuasiStatic]
      [./all]
        strain = FINITE
        add_variables = true
        generate_output = 'stress_zz'
      [../]
    [../]
  [../]
[]

[BCs]
  [./back_z]
    type = DirichletBC
    variable = disp_z
    boundary = back
    value = 0.0
  [../]
  [./left_x]
    type = DirichletBC
    variable = disp_x
    boundary = left
    value = 0.0
  [../]
  [./bottom_y]
    type = DirichletBC
    variable = disp_y
    boundary = bottom
    value = 0.0
  [../]
  [./move_front]
    type = FunctionDirichletBC
    variable = disp_z
    boundary = front
    function = 't/10.'
  [../]
[]

[Materials]
  [./stress]
    type = ComputeFiniteStrainElasticStress
  [../]
  [./elasticity_tensor]
    type = ComputeElasticityTensor
    fill_method = symmetric9
    C_ijkl = '1.5e6 0.75e6 0.75e6 1.5e6 0.75e6 1.5e6 0.375e6 0.375e6 0.375e6'
  [../]
[]

[Postprocessors]
  [./szz_avg]
    type =MaterialTensorAverage
    rank_two_tensor = stress
    index_i = 2
    index_j = 2
    use_displaced_mesh = true
  []
  [./szz_int]
    type =MaterialTensorIntegral
    rank_two_tensor = stress
    index_i = 2
    index_j = 2
    use_displaced_mesh = true
  []
  [./szz_avg_aux]
    type =ElementAverageValue
    variable = stress_zz
    use_displaced_mesh = true
  []
  [./szz_int_aux]
    type =ElementIntegralVariablePostprocessor
    variable = stress_zz
    use_displaced_mesh = true
  []
[]


[Preconditioning]
  [./SMP]
    type = SMP
    full = true
  [../]
[]

[Executioner]
  type = Transient
  solve_type = NEWTON
  petsc_options_iname = '-pc_type '
  petsc_options_value = lu
  nl_rel_tol = 1e-10
  nl_abs_tol = 1e-6
  l_max_its = 20
  start_time = 0.0
  dt = 0.2
  end_time = 1.0
[]

[Outputs]
  csv = true
[]

Description
Example Input File
Input Parameters
Input Files