Layered2D System

Description

The Layered2D Action is a convenience object that simplifies part of the mechanics system setup for Layered2D models. It performs the following tasks:

Add StressDivergenceTensors and GeneralizedPlaneStrainOffDiag Kernels
Add GeneralizedPlaneStrain ScalarKernels
Add GeneralizedPlaneStrainUserObject UserObjects
Add Strain calculation material (for the chosen strain model)
Add ScalarStrainAux AuxKernel to assemble the scalar out-of-plane strain(s) into one single given out-of-plane strain AuxVariable
Correctly sets the use_displaced_mesh for all objects set up by this Action
Optional: Setup of in-plane displacement variable (with the correct order for the current mesh)
Optional: Setup of scalar out-of-plane-strain variable(s)
Optional: Add coupling between the temperature variable with the out-of-plane strain scalar variable(s) in the GeneralizedPlaneStrainOffDiag Kernel
Optional: Setup of the GeneralizedPlaneStrainReferenceResidual problem

Table 1: Correspondence Among Action Functionality and MooseObjects for the Layered2D Action

Functionality	Replaced Classes	Associated Parameters
In-plane equilibrium conditions for XYZ coordinate system	Stress Divergence Tensors Kernel	`displacements` : a string of the displacement field variables
Scalar out-of-plane strain coupling with in-plane field variables	Generalized Plane Strain Off-diagonal Kernel	`scalar_out_of_plane_strain`: a list of the scalar variables for the out-of-plane strain direction
\|	`displacements` : a string of the displacement field (in-plane) variable
\|	`temperature`: a string of the temperature field variable
Out-of-plane scalar variable equilibrium condition	Generalized Plane Strain ScalarKernel	`scalar_out_of_plane_strain`: a list of the scalar variables for the out-of-plane strain direction
Residual and diagonal Jacobian calculation for scalar out-of-plane strain variables	Generalized Plane Strain UserObject	`scalar_out_of_plane_strain`: a list of the scalar variables for the out-of-plane strain direction
Add the displacement variable	Variables	`add_variables`: boolean
Add the out-of-plane strain scalar variable(s)	ScalarVariables	`add_scalar_variables`: boolean
Calculation of in-plane strain	Compute Plane Finite Strain	`strain`: MooseEnum to select finite or small strain formulations
\| Compute Plane Small Strain	\|
\| Compute Axisymmetric Plane Incremental Strain \|`incremental` : boolean for using a incremental strain formulation
Add AuxVariable and AuxKernel to gather scalar out-of-plane strains\| AuxVariable and ScalarStrainAux	`out_of_plane_strain_name`: Name provided by user to gather scalar_out_of_plane_strains
Add AuxScalarKernel for reference residual problem	Generalized Plane Strain Reference Residual AuxScalarKernel\| `save_in`: List of names of auxiliary variables to save this Kernel's residual contributions to\|

Example Input File Syntax

[Physics<<<{"href": "../../index.html"}>>>]
  [SolidMechanics<<<{"href": "../index.html"}>>>]
    [Layered2D<<<{"href": "index.html"}>>>]
      [fuel]
        add_scalar_variables<<<{"description": "Add the scalar_out_of_plane_strain variables."}>>> = true
        add_variables<<<{"description": "Add the displacement variables"}>>> = true
        out_of_plane_strain_name<<<{"description": "Name provided by user for aux variable to gather scalar_out_of_plane_strains."}>>> = strain_yy
        out_of_plane_direction<<<{"description": "The direction of the out-of-plane strain."}>>> = y
        temperature<<<{"description": "The temperature"}>>> = temp
        fuel_pin_geometry<<<{"description": "User object name which provides subblock index."}>>> = pin_geometry
        eigenstrain_names<<<{"description": "List of eigenstrains to be applied in this strain calculation"}>>> = 'fuel_thermal_strain'
        strain<<<{"description": "Strain formulation"}>>> = finite
        block<<<{"description": "The list of ids of the blocks (subdomain) that the stress divergence kernels will be applied to"}>>> = fuel
        mesh_generator<<<{"description": "The name of the generator to use as the prefix for mesh meta data properties."}>>> = layered2D_mesh
      []
      [clad]
        add_scalar_variables<<<{"description": "Add the scalar_out_of_plane_strain variables."}>>> = true
        add_variables<<<{"description": "Add the displacement variables"}>>> = true
        out_of_plane_strain_name<<<{"description": "Name provided by user for aux variable to gather scalar_out_of_plane_strains."}>>> = strain_yy
        out_of_plane_direction<<<{"description": "The direction of the out-of-plane strain."}>>> = y
        temperature<<<{"description": "The temperature"}>>> = temp
        fuel_pin_geometry<<<{"description": "User object name which provides subblock index."}>>> = pin_geometry
        eigenstrain_names<<<{"description": "List of eigenstrains to be applied in this strain calculation"}>>> = 'clad_thermal_strain'
        strain<<<{"description": "Strain formulation"}>>> = finite
        block<<<{"description": "The list of ids of the blocks (subdomain) that the stress divergence kernels will be applied to"}>>> = clad
        mesh_generator<<<{"description": "The name of the generator to use as the prefix for mesh meta data properties."}>>> = layered2D_mesh
      []
    []
  []
[]

Subblocks

The subblocks of the Layered2D action are what trigger MOOSE objects for Layered2D models to be built. If a Layered2D model can be applied for the whole simulation domain, i.e. the whole simulation has one scalar out-of-plane strain, a single subblock should be used.

If different a subdomain has a different Layered2D model, multiple subblocks with subdomain restrictions can be used. Generally, the number of subblocks is the same as the number of scalar out-of-plane strains the simulation contains.

Parameters supplied in the [Physics/SolidMechanics/Layered2D] block act as defaults for the QuasiStatic Action subblocks.

Input Parameters

fuel_pin_geometryUser object name which provides subblock index.
C++ Type:UserObjectName
Controllable:No
Description:User object name which provides subblock index.
mesh_generatorThe name of the generator to use as the prefix for mesh meta data properties.
C++ Type:MeshGeneratorName
Controllable:No
Description:The name of the generator to use as the prefix for mesh meta data properties.

Required Parameters

absolute_value_vector_tagsThe tag names for extra vectors that the absolute value of the residual should be accumulated into
C++ Type:std::vector<TagName>
Controllable:No
Description:The tag names for extra vectors that the absolute value of the residual should be accumulated into
active__all__ If specified only the blocks named will be visited and made active
Default:__all__
C++ Type:std::vector<std::string>
Controllable:No
Description:If specified only the blocks named will be visited and made active
add_scalar_variablesFalseAdd the scalar_out_of_plane_strain variables.
Default:False
C++ Type:bool
Controllable:No
Description:Add the scalar_out_of_plane_strain variables.
base_nameMaterial property base name
C++ Type:std::string
Controllable:No
Description:Material property base name
constraint_typesType of each constraint: stress or strain.
C++ Type:MultiMooseEnum
Options:strain, stress, none
Controllable:No
Description:Type of each constraint: stress or strain.
decomposition_methodTaylorExpansionMethods to calculate the finite strain and rotation increments
Default:TaylorExpansion
C++ Type:MooseEnum
Options:TaylorExpansion, EigenSolution, HughesWinget
Controllable:No
Description:Methods to calculate the finite strain and rotation increments
extra_vector_tagsThe tag names for extra vectors that residual data should be saved into
C++ Type:std::vector<TagName>
Controllable:No
Description:The tag names for extra vectors that residual data should be saved into
formulationTOTALSelect between the total Lagrangian (TOTAL) and updated Lagrangian (UPDATED) formulations for the new kernel system.
Default:TOTAL
C++ Type:MooseEnum
Options:TOTAL, UPDATED
Controllable:No
Description:Select between the total Lagrangian (TOTAL) and updated Lagrangian (UPDATED) formulations for the new kernel system.
global_strainName of the global strain material to be applied in this strain calculation. The global strain tensor is constant over the whole domain and allows visualization of the deformed shape with the periodic BC
C++ Type:MaterialPropertyName
Unit:(no unit assumed)
Controllable:No
Description:Name of the global strain material to be applied in this strain calculation. The global strain tensor is constant over the whole domain and allows visualization of the deformed shape with the periodic BC
group_scalar_vars_in_reference_residualFalseGroup all scalar variables in reference residual problem.
Default:False
C++ Type:bool
Controllable:No
Description:Group all scalar variables in reference residual problem.
inactiveIf specified blocks matching these identifiers will be skipped.
C++ Type:std::vector<std::string>
Controllable:No
Description:If specified blocks matching these identifiers will be skipped.
initial_eigenstrain_nameA single eigenstrain name that contains the initial strains applied to a problem.
C++ Type:MaterialPropertyName
Unit:(no unit assumed)
Controllable:No
Description:A single eigenstrain name that contains the initial strains applied to a problem.
new_systemFalseIf true use the new LagrangianStressDiverence kernels.
Default:False
C++ Type:bool
Controllable:No
Description:If true use the new LagrangianStressDiverence kernels.
out_of_plane_strain_namescalar_strain_yyName provided by user for aux variable to gather scalar_out_of_plane_strains.
Default:scalar_strain_yy
C++ Type:AuxVariableName
Unit:(no unit assumed)
Controllable:No
Description:Name provided by user for aux variable to gather scalar_out_of_plane_strains.
scalar_out_of_plane_strain_variablesList of scalar_out_of_plane_strain variables.
C++ Type:std::vector<NonlinearVariableName>
Unit:(no unit assumed)
Controllable:No
Description:List of scalar_out_of_plane_strain variables.
targetsFunctions giving the target values of each constraint.
C++ Type:std::vector<FunctionName>
Unit:(no unit assumed)
Controllable:No
Description:Functions giving the target values of each constraint.
use_automatic_differentiationFalseFlag to use automatic differentiation (AD) objects when possible
Default:False
C++ Type:bool
Controllable:No
Description:Flag to use automatic differentiation (AD) objects when possible
verboseFalseDisplay extra information.
Default:False
C++ Type:bool
Controllable:No
Description:Display extra information.
volumetric_locking_correctionFalseFlag to correct volumetric locking
Default:False
C++ Type:bool
Controllable:No
Description:Flag to correct volumetric locking

Optional Parameters

add_variablesFalseAdd the displacement variables
Default:False
C++ Type:bool
Controllable:No
Description:Add the displacement variables
displacementsThe nonlinear displacement variables for the problem
C++ Type:std::vector<VariableName>
Unit:(no unit assumed)
Controllable:No
Description:The nonlinear displacement variables for the problem
scalingThe scaling to apply to the displacement variables
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:The scaling to apply to the displacement variables
temperatureThe temperature
C++ Type:std::vector<VariableName>
Unit:(no unit assumed)
Controllable:No
Description:The temperature

Variables Parameters

additional_generate_outputAdd scalar quantity output for stress and/or strain (will be appended to the list in `generate_output`)
C++ Type:MultiMooseEnum
Options:cauchy_stress_xx, cauchy_stress_xy, cauchy_stress_xz, cauchy_stress_yx, cauchy_stress_yy, cauchy_stress_yz, cauchy_stress_zx, cauchy_stress_zy, cauchy_stress_zz, creep_strain_xx, creep_strain_xy, creep_strain_xz, creep_strain_yx, creep_strain_yy, creep_strain_yz, creep_strain_zx, creep_strain_zy, creep_strain_zz, creep_stress_xx, creep_stress_xy, creep_stress_xz, creep_stress_yx, creep_stress_yy, creep_stress_yz, creep_stress_zx, creep_stress_zy, creep_stress_zz, deformation_gradient_xx, deformation_gradient_xy, deformation_gradient_xz, deformation_gradient_yx, deformation_gradient_yy, deformation_gradient_yz, deformation_gradient_zx, deformation_gradient_zy, deformation_gradient_zz, elastic_strain_xx, elastic_strain_xy, elastic_strain_xz, elastic_strain_yx, elastic_strain_yy, elastic_strain_yz, elastic_strain_zx, elastic_strain_zy, elastic_strain_zz, mechanical_strain_xx, mechanical_strain_xy, mechanical_strain_xz, mechanical_strain_yx, mechanical_strain_yy, mechanical_strain_yz, mechanical_strain_zx, mechanical_strain_zy, mechanical_strain_zz, pk1_stress_xx, pk1_stress_xy, pk1_stress_xz, pk1_stress_yx, pk1_stress_yy, pk1_stress_yz, pk1_stress_zx, pk1_stress_zy, pk1_stress_zz, pk2_stress_xx, pk2_stress_xy, pk2_stress_xz, pk2_stress_yx, pk2_stress_yy, pk2_stress_yz, pk2_stress_zx, pk2_stress_zy, pk2_stress_zz, plastic_strain_xx, plastic_strain_xy, plastic_strain_xz, plastic_strain_yx, plastic_strain_yy, plastic_strain_yz, plastic_strain_zx, plastic_strain_zy, plastic_strain_zz, small_stress_xx, small_stress_xy, small_stress_xz, small_stress_yx, small_stress_yy, small_stress_yz, small_stress_zx, small_stress_zy, small_stress_zz, strain_xx, strain_xy, strain_xz, strain_yx, strain_yy, strain_yz, strain_zx, strain_zy, strain_zz, stress_xx, stress_xy, stress_xz, stress_yx, stress_yy, stress_yz, stress_zx, stress_zy, stress_zz, effective_plastic_strain, effective_creep_strain, firstinv_stress, firstinv_cauchy_stress, firstinv_pk1_stress, firstinv_pk2_stress, firstinv_small_stress, firstinv_strain, hydrostatic_stress, hydrostatic_cauchy_stress, hydrostatic_pk1_stress, hydrostatic_pk2_stress, hydrostatic_small_stress, intensity_stress, intensity_cauchy_stress, intensity_pk1_stress, intensity_pk2_stress, intensity_small_stress, l2norm_mechanical_strain, l2norm_stress, l2norm_cauchy_stress, l2norm_pk1_stress, l2norm_strain, l2norm_elastic_strain, l2norm_plastic_strain, l2norm_creep_strain, max_principal_mechanical_strain, max_principal_stress, max_principal_cauchy_stress, max_principal_pk1_stress, max_principal_pk2_stress, max_principal_small_stress, max_principal_strain, maxshear_stress, maxshear_cauchy_stress, maxshear_pk1_stress, maxshear_pk2_stress, maxshear_small_stress, mid_principal_mechanical_strain, mid_principal_stress, mid_principal_cauchy_stress, mid_principal_pk1_stress, mid_principal_pk2_stress, mid_principal_small_stress, mid_principal_strain, min_principal_mechanical_strain, min_principal_stress, min_principal_cauchy_stress, min_principal_pk1_stress, min_principal_pk2_stress, min_principal_small_stress, min_principal_strain, secondinv_stress, secondinv_cauchy_stress, secondinv_pk1_stress, secondinv_pk2_stress, secondinv_small_stress, secondinv_strain, thirdinv_stress, thirdinv_cauchy_stress, thirdinv_pk1_stress, thirdinv_pk2_stress, thirdinv_small_stress, thirdinv_strain, triaxiality_stress, triaxiality_cauchy_stress, triaxiality_pk1_stress, triaxiality_pk2_stress, triaxiality_small_stress, volumetric_mechanical_strain, volumetric_strain, vonmises_stress, vonmises_cauchy_stress, vonmises_pk1_stress, vonmises_pk2_stress, directional_stress, directional_strain, axial_stress, axial_strain, axial_plastic_strain, axial_creep_strain, axial_elastic_strain, hoop_stress, hoop_strain, hoop_plastic_strain, hoop_creep_strain, hoop_elastic_strain, radial_stress, radial_strain, spherical_hoop_stress, spherical_hoop_strain, spherical_hoop_plastic_strain, spherical_hoop_creep_strain, spherical_hoop_elastic_strain, spherical_radial_stress, spherical_radial_strain
Controllable:No
Description:Add scalar quantity output for stress and/or strain (will be appended to the list in `generate_output`)
additional_material_output_familySpecifies the family of FE shape functions to use for this variable.
C++ Type:MultiMooseEnum
Options:MONOMIAL, LAGRANGE
Controllable:No
Description:Specifies the family of FE shape functions to use for this variable.
additional_material_output_orderSpecifies the order of the FE shape function to use for this variable.
C++ Type:MultiMooseEnum
Options:CONSTANT, FIRST, SECOND, THIRD, FOURTH, FIFTH, SIXTH, SEVENTH, EIGHTH, NINTH
Controllable:No
Description:Specifies the order of the FE shape function to use for this variable.
generate_outputAdd scalar quantity output for stress and/or strain
C++ Type:MultiMooseEnum
Options:cauchy_stress_xx, cauchy_stress_xy, cauchy_stress_xz, cauchy_stress_yx, cauchy_stress_yy, cauchy_stress_yz, cauchy_stress_zx, cauchy_stress_zy, cauchy_stress_zz, creep_strain_xx, creep_strain_xy, creep_strain_xz, creep_strain_yx, creep_strain_yy, creep_strain_yz, creep_strain_zx, creep_strain_zy, creep_strain_zz, creep_stress_xx, creep_stress_xy, creep_stress_xz, creep_stress_yx, creep_stress_yy, creep_stress_yz, creep_stress_zx, creep_stress_zy, creep_stress_zz, deformation_gradient_xx, deformation_gradient_xy, deformation_gradient_xz, deformation_gradient_yx, deformation_gradient_yy, deformation_gradient_yz, deformation_gradient_zx, deformation_gradient_zy, deformation_gradient_zz, elastic_strain_xx, elastic_strain_xy, elastic_strain_xz, elastic_strain_yx, elastic_strain_yy, elastic_strain_yz, elastic_strain_zx, elastic_strain_zy, elastic_strain_zz, mechanical_strain_xx, mechanical_strain_xy, mechanical_strain_xz, mechanical_strain_yx, mechanical_strain_yy, mechanical_strain_yz, mechanical_strain_zx, mechanical_strain_zy, mechanical_strain_zz, pk1_stress_xx, pk1_stress_xy, pk1_stress_xz, pk1_stress_yx, pk1_stress_yy, pk1_stress_yz, pk1_stress_zx, pk1_stress_zy, pk1_stress_zz, pk2_stress_xx, pk2_stress_xy, pk2_stress_xz, pk2_stress_yx, pk2_stress_yy, pk2_stress_yz, pk2_stress_zx, pk2_stress_zy, pk2_stress_zz, plastic_strain_xx, plastic_strain_xy, plastic_strain_xz, plastic_strain_yx, plastic_strain_yy, plastic_strain_yz, plastic_strain_zx, plastic_strain_zy, plastic_strain_zz, small_stress_xx, small_stress_xy, small_stress_xz, small_stress_yx, small_stress_yy, small_stress_yz, small_stress_zx, small_stress_zy, small_stress_zz, strain_xx, strain_xy, strain_xz, strain_yx, strain_yy, strain_yz, strain_zx, strain_zy, strain_zz, stress_xx, stress_xy, stress_xz, stress_yx, stress_yy, stress_yz, stress_zx, stress_zy, stress_zz, effective_plastic_strain, effective_creep_strain, firstinv_stress, firstinv_cauchy_stress, firstinv_pk1_stress, firstinv_pk2_stress, firstinv_small_stress, firstinv_strain, hydrostatic_stress, hydrostatic_cauchy_stress, hydrostatic_pk1_stress, hydrostatic_pk2_stress, hydrostatic_small_stress, intensity_stress, intensity_cauchy_stress, intensity_pk1_stress, intensity_pk2_stress, intensity_small_stress, l2norm_mechanical_strain, l2norm_stress, l2norm_cauchy_stress, l2norm_pk1_stress, l2norm_strain, l2norm_elastic_strain, l2norm_plastic_strain, l2norm_creep_strain, max_principal_mechanical_strain, max_principal_stress, max_principal_cauchy_stress, max_principal_pk1_stress, max_principal_pk2_stress, max_principal_small_stress, max_principal_strain, maxshear_stress, maxshear_cauchy_stress, maxshear_pk1_stress, maxshear_pk2_stress, maxshear_small_stress, mid_principal_mechanical_strain, mid_principal_stress, mid_principal_cauchy_stress, mid_principal_pk1_stress, mid_principal_pk2_stress, mid_principal_small_stress, mid_principal_strain, min_principal_mechanical_strain, min_principal_stress, min_principal_cauchy_stress, min_principal_pk1_stress, min_principal_pk2_stress, min_principal_small_stress, min_principal_strain, secondinv_stress, secondinv_cauchy_stress, secondinv_pk1_stress, secondinv_pk2_stress, secondinv_small_stress, secondinv_strain, thirdinv_stress, thirdinv_cauchy_stress, thirdinv_pk1_stress, thirdinv_pk2_stress, thirdinv_small_stress, thirdinv_strain, triaxiality_stress, triaxiality_cauchy_stress, triaxiality_pk1_stress, triaxiality_pk2_stress, triaxiality_small_stress, volumetric_mechanical_strain, volumetric_strain, vonmises_stress, vonmises_cauchy_stress, vonmises_pk1_stress, vonmises_pk2_stress, directional_stress, directional_strain, axial_stress, axial_strain, axial_plastic_strain, axial_creep_strain, axial_elastic_strain, hoop_stress, hoop_strain, hoop_plastic_strain, hoop_creep_strain, hoop_elastic_strain, radial_stress, radial_strain, spherical_hoop_stress, spherical_hoop_strain, spherical_hoop_plastic_strain, spherical_hoop_creep_strain, spherical_hoop_elastic_strain, spherical_radial_stress, spherical_radial_strain
Controllable:No
Description:Add scalar quantity output for stress and/or strain
material_output_familySpecifies the family of FE shape functions to use for this variable.
C++ Type:MultiMooseEnum
Options:MONOMIAL, LAGRANGE
Controllable:No
Description:Specifies the family of FE shape functions to use for this variable.
material_output_orderSpecifies the order of the FE shape function to use for this variable.
C++ Type:MultiMooseEnum
Options:CONSTANT, FIRST, SECOND, THIRD, FOURTH, FIFTH, SIXTH, SEVENTH, EIGHTH, NINTH
Controllable:No
Description:Specifies the order of the FE shape function to use for this variable.

Output Parameters

automatic_eigenstrain_namesFalseCollects all material eigenstrains and passes to required strain calculator within TMA internally.
Default:False
C++ Type:bool
Controllable:No
Description:Collects all material eigenstrains and passes to required strain calculator within TMA internally.
eigenstrain_namesList of eigenstrains to be applied in this strain calculation
C++ Type:std::vector<MaterialPropertyName>
Unit:(no unit assumed)
Controllable:No
Description:List of eigenstrains to be applied in this strain calculation
incrementalFalseUse incremental or total strain (if not explicitly specified this defaults to incremental for finite strain and total for small strain)
Default:False
C++ Type:bool
Controllable:No
Description:Use incremental or total strain (if not explicitly specified this defaults to incremental for finite strain and total for small strain)
strainSMALLStrain formulation
Default:SMALL
C++ Type:MooseEnum
Options:SMALL, FINITE
Controllable:No
Description:Strain formulation
strain_base_nameThe base name used for the strain. If not provided, it will be set equal to base_name
C++ Type:std::string
Controllable:No
Description:The base name used for the strain. If not provided, it will be set equal to base_name
use_finite_deform_jacobianFalseJacobian for corrotational finite strain
Default:False
C++ Type:bool
Controllable:No
Description:Jacobian for corrotational finite strain

Strain Parameters

blockThe list of ids of the blocks (subdomain) that the stress divergence kernels will be applied to
C++ Type:std::vector<SubdomainName>
Controllable:No
Description:The list of ids of the blocks (subdomain) that the stress divergence kernels will be applied to
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.
diag_save_inThe displacement diagonal preconditioner terms
C++ Type:std::vector<AuxVariableName>
Unit:(no unit assumed)
Controllable:No
Description:The displacement diagonal preconditioner terms
outputsnone Vector of output names where you would like to restrict the output of variables(s) associated with this object
Default:none
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
save_inThe displacement residuals
C++ Type:std::vector<AuxVariableName>
Unit:(no unit assumed)
Controllable:No
Description:The displacement residuals

Advanced Parameters

cylindrical_axis_point1Starting point for direction of axis of rotation for cylindrical stress/strain.
C++ Type:libMesh::Point
Controllable:No
Description:Starting point for direction of axis of rotation for cylindrical stress/strain.
cylindrical_axis_point2Ending point for direction of axis of rotation for cylindrical stress/strain.
C++ Type:libMesh::Point
Controllable:No
Description:Ending point for direction of axis of rotation for cylindrical stress/strain.
directionDirection stress/strain is calculated in
C++ Type:libMesh::Point
Controllable:No
Description:Direction stress/strain is calculated in
spherical_center_pointCenter point of the spherical coordinate system.
C++ Type:libMesh::Point
Controllable:No
Description:Center point of the spherical coordinate system.

Coordinate System Parameters

out_of_plane_directionzThe direction of the out-of-plane strain.
Default:z
C++ Type:MooseEnum
Options:x, y, z
Controllable:No
Description:The direction of the out-of-plane strain.
out_of_plane_pressure_functionFunction used to prescribe pressure (applied toward the body) in the out-of-plane direction (y for 1D Axisymmetric or z for 2D Cartesian problems)
C++ Type:FunctionName
Unit:(no unit assumed)
Controllable:No
Description:Function used to prescribe pressure (applied toward the body) in the out-of-plane direction (y for 1D Axisymmetric or z for 2D Cartesian problems)
out_of_plane_pressure_material0Material used to prescribe pressure (applied toward the body) in the out-of-plane direction
Default:0
C++ Type:MaterialPropertyName
Unit:(no unit assumed)
Controllable:No
Description:Material used to prescribe pressure (applied toward the body) in the out-of-plane direction
out_of_plane_strainVariable for the out-of-plane strain for plane stress models
C++ Type:VariableName
Unit:(no unit assumed)
Controllable:No
Description:Variable for the out-of-plane strain for plane stress models
planar_formulationNONEOut-of-plane stress/strain formulation
Default:NONE
C++ Type:MooseEnum
Options:NONE, WEAK_PLANE_STRESS, PLANE_STRAIN, GENERALIZED_PLANE_STRAIN
Controllable:No
Description:Out-of-plane stress/strain formulation
pressure_factorScale factor applied to prescribed out-of-plane pressure (both material and function)
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Scale factor applied to prescribed out-of-plane pressure (both material and function)
scalar_out_of_plane_strainScalar variable for the out-of-plane strain (in y direction for 1D Axisymmetric or in z direction for 2D Cartesian problems)
C++ Type:VariableName
Unit:(no unit assumed)
Controllable:No
Description:Scalar variable for the out-of-plane strain (in y direction for 1D Axisymmetric or in z direction for 2D Cartesian problems)

Out-Of-Plane Stress/Strain Parameters

Associated Actions

Available Actions

Bison App
Layered2DActionSets up (Aux)variables, materials and (Aux)kernels for layered two dimensional simulations.

(test/tests/layered2D/multi_block.i)

#
# This test checks that the Layered2D master action properly creates the strain
# calculator, scalar variables and kernels, and aux variable that groups all of
# the scalar variables into a single variable for output purposes.
#
# The simulation domain contains a single layer of fuel and cladding and one
# layer for the plenum. The fuel and cladding layer contains a circle of fuel
# with a ring of cladding and the plenum layer contains only a ring of cladding.
# The initial temperature of both the fuel and cladding is set to 300.0 K and
# is linearly increased to 500.0 K over 2 seconds.  The thermal expansion
# coefficientso are taken as constant at 10e-6 K^-1 and 5e-6 K^-1 for the fuel
# and cladding respectively.  Assuming isotropic linear thermal the calculated
# out-of-plane strain due to thermal expansion in each block is:
#
# epsilon_yy = alpha (T-T_o)
#
# where epsilon_yy is the out-of-plane strain, alpha is the thermal expansion
# coefficient, and T_o is the stress free temperature.
#
# Therefore in the fuel the out-of-plane strain is:
#
# epsilon_yy = 10.0e-6 * (500-300) = 2.0e-3.
#
# and in the cladding:
#
# epsilon_yy = 5.0e-6 * (500-300) = 1.0e-3
#
# These answers are verified by the postprocessors.

[GlobalParams]
  order = SECOND
  family = LAGRANGE
  displacements = 'disp_x disp_y disp_z'
[]

[Mesh]
  [layered2D_mesh]
    type = Layered2DMeshGenerator
    num_sectors = 20
    axial_direction = y
    slices_per_block = 1
    clad_thickness = 1
    pellet_mesh_density = coarse
    clad_mesh_density = coarse
    clad_gap_width = 0.414213562
    fuel_height = 0.318309886 # 1/pi
    plenum_height = 0.318309886 # 1/pi
    pellet_bottom_coor = 0
    pellet_inner_radius = 0
    pellet_outer_radius = 1
  []
[]

[Variables]
  [temp]
    initial_condition = 300.0
  []
[]

[Functions]
  [temp_ramp]
    type = PiecewiseLinear
    xy_data = '0.0 300.0
               2.0 500.0'
  []
[]


[UserObjects]
  [pin_geometry]
    type = Layered2DFuelPinGeometry
    mesh_generator = layered2D_mesh
  []
[]

[Physics]
  [SolidMechanics]
    [Layered2D]
      [fuel]
        add_scalar_variables = true
        add_variables = true
        out_of_plane_strain_name = strain_yy
        out_of_plane_direction = y
        temperature = temp
        fuel_pin_geometry = pin_geometry
        eigenstrain_names = 'fuel_thermal_strain'
        strain = finite
        block = fuel
        mesh_generator = layered2D_mesh
      []
      [clad]
        add_scalar_variables = true
        add_variables = true
        out_of_plane_strain_name = strain_yy
        out_of_plane_direction = y
        temperature = temp
        fuel_pin_geometry = pin_geometry
        eigenstrain_names = 'clad_thermal_strain'
        strain = finite
        block = clad
        mesh_generator = layered2D_mesh
      []
    []
  []
[]

[Kernels]
  [heat]
    type = HeatConduction
    variable = temp
  []
[]

[BCs]
  [no_x]
    type = DirichletBC
    boundary = bottom
    variable = disp_x
    value = 0.0
  []
  [no_z]
    type = DirichletBC
    boundary = left
    variable = disp_z
    value = 0.0
  []

  [temp_ramp]
    type = FunctionDirichletBC
    variable = temp
    boundary = '10 2 5'
    function = temp_ramp
  []
[]

[Materials]
  [elasticity_tensor]
    type = ComputeIsotropicElasticityTensor
    block = 'fuel clad'
    youngs_modulus = 10
    poissons_ratio = 0
  []
  [thermal_material]
    type = HeatConductionMaterial
    thermal_conductivity = 1.0e6
    specific_heat = 1.0
    block = 'fuel clad'
  []
  [thermal_expansion_fuel]
    type = ComputeThermalExpansionEigenstrain
    block = fuel
    temperature = temp
    thermal_expansion_coeff = 10e-6
    stress_free_temperature = 300.0
    eigenstrain_name = 'fuel_thermal_strain'
  []
  [stress_fuel]
    type = ComputeFiniteStrainElasticStress
    block = 'fuel'
  []
  [thermal_expansion_clad]
    type = ComputeThermalExpansionEigenstrain
    block = clad
    temperature = temp
    thermal_expansion_coeff = 5e-6
    stress_free_temperature = 300.0
    eigenstrain_name = 'clad_thermal_strain'
  []
  [stress_clad]
    type = ComputeFiniteStrainElasticStress
    block = 'clad'
  []
[]

[Executioner]
  type = Transient

  solve_type = 'PJFNK'

  petsc_options_iname = '-pc_type -pc_factor_mat_solver_package'
  petsc_options_value = 'lu superlu_dist'

  line_search = 'none'

  l_max_its = 50
  l_tol = 1e-5
  nl_max_its = 15
  nl_abs_tol = 1e-10
  nl_rel_tol = 1e-10
  end_time = 2.0
  dt = 1
[]

[Postprocessors]
  [strain_yy_fuel]
    type = ElementAverageValue
    block = fuel
    variable = strain_yy
  []
  [strain_yy_clad]
    type = ElementAverageValue
    block = clad
    variable = strain_yy
  []
[]

[Outputs]
  exodus = true
  hide = 'disp_x disp_z'
[]

Description
Example Input File Syntax
Input Parameters
Associated Actions
Available Actions