EqualValueEmbeddedConstraint

The EqualValueEmbeddedConstraint block is used to construct a set of EqualValueEmbeddedConstraint objects for applying these constraints on multiple sets of blocks and on multiple variables with minimal input. These individual EqualValueEmbeddedConstraint objects are used to define constraints that tie the values of variables at the nodes in a single block to the containing elements in another block. Prescribing this input can be tedious of these constraints are to be defined for multiple blocks and variables, because many input blocks would be required. The EqualValueEmbeddedConstraint block simplifies that process by allowing for multiple blocks and directions to be defined in a single input file block.

Constructed Objects

This block sets up multiple EqualValueEmbeddedConstraint objects.

Example Input Syntax

[Constraints/EqualValueEmbeddedConstraint/Batch]
  [Constraints/EqualValueEmbeddedConstraint/Batch]
    [Constraints/EqualValueEmbeddedConstraint/Batch]
      primary = '1'
      secondary = '2 3'
      variable = 'disp_x disp_y disp_z'
      penalty = 1e12
      formulation = penalty
    []
  []
[]

Input Parameters

penaltyPenalty parameter used in constraint enforcement for kinematic and penalty formulations.
C++ Type:double
Options:
Description:Penalty parameter used in constraint enforcement for kinematic and penalty formulations.
primaryThe list of ids of the primary block (subdomain)
C++ Type:std::vector<SubdomainName>
Options:
Description:The list of ids of the primary block (subdomain)
secondaryThe list of ids of the secondary block (subdomain)
C++ Type:std::vector<SubdomainName>
Options:
Description:The list of ids of the secondary block (subdomain)
variableThe list of variables on the secondary side
C++ Type:std::vector<NonlinearVariableName>
Options:
Description:The list of variables on the secondary side

Required Parameters

active__all__ If specified only the blocks named will be visited and made active
Default:__all__
C++ Type:std::vector<std::string>
Options:
Description:If specified only the blocks named will be visited and made active
formulationkinematicFormulation used to enforce the constraint
Default:kinematic
C++ Type:MooseEnum
Options:kinematic, penalty
Description:Formulation used to enforce the constraint
inactiveIf specified blocks matching these identifiers will be skipped.
C++ Type:std::vector<std::string>
Options:
Description:If specified blocks matching these identifiers will be skipped.
primary_variableThe list of variables on the primary side (defaults to the specified value of 'variable')
C++ Type:std::vector<VariableName>
Options:
Description:The list of variables on the primary side (defaults to the specified value of 'variable')

Optional Parameters

References

(test/tests/actions/EqualValueEmbeddedConstraintAction.i)

[Mesh]
    file = reinforcedBeam.e
[]

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

[Modules/TensorMechanics/Master]
  [./Concrete_block]
    block = 1
    strain = small
    incremental = true
    generate_output = 'stress_xx stress_xy stress_xz stress_yy stress_yz stress_zz
               strain_xx strain_xy strain_xz strain_yy strain_yz strain_zz
               max_principal_stress mid_principal_stress min_principal_stress
               secondinv_stress thirdinv_stress vonmises_stress
               secondinv_strain thirdinv_strain
               elastic_strain_xx elastic_strain_xy elastic_strain_xz elastic_strain_yy elastic_strain_yz elastic_strain_zz'
    save_in = 'resid_x resid_y resid_z'
  [../]
[]

[Modules/TensorMechanics/LineElementMaster]
  [./Reinforcement_block]
    block = '2 3'
    truss = true
    area = area
    displacements = 'disp_x disp_y disp_z'
    save_in = 'resid_x resid_y resid_z'
  [../]
[]

[Constraints/EqualValueEmbeddedConstraint/Batch]
  primary = '1'
  secondary = '2 3'
  variable = 'disp_x disp_y disp_z'
  penalty = 1e12
  formulation = penalty
[]

[Variables]
  [./disp_x]
  [../]
  [./disp_y]
  [../]
  [./disp_z]
  [../]
[]

[AuxVariables]
  [./resid_x]
  [../]
  [./resid_y]
  [../]
  [./resid_z]
  [../]
  [./area]
    order = CONSTANT
    family = MONOMIAL
  [../]
  [./D]
    order = CONSTANT
    family = MONOMIAL
  [../]
  [./intnl0]
    order = CONSTANT
    family = MONOMIAL
  [../]
  [./intnl1]
    order = CONSTANT
    family = MONOMIAL
  [../]
[]

[AuxKernels]
  [./area2]
    type = ConstantAux
    block = 2
    variable = area
    value = 5.09e-4 # 509 mm2
    execute_on = 'initial timestep_begin'
  [../]
  [./area3]
    type = ConstantAux
    block = 3
    variable = area
    value = 2.00e-4 # 200 mm2
    execute_on = 'initial timestep_begin'
  [../]
[]

[BCs]
  [./loading]
    type = FunctionDirichletBC
    variable = disp_x
    boundary = '2'
    function = '1E-5*t'
  [../]
  [./left_support_x]
    type = FunctionDirichletBC
    variable = disp_x
    boundary = '1'
    function = '0.'
  [../]
  [./left_support_y]
    type = FunctionDirichletBC
    variable = disp_y
    boundary = '1'
    function = '0.'
  [../]
  [./left_support_z]
    type = FunctionDirichletBC
    variable = disp_z
    boundary = '1'
    function = '0.'
  [../]
[]

[Materials]
  # Material Properties of concrete
  [./elasticity_tensor_concrete]
    type = ComputeIsotropicElasticityTensor
    block = 1
    youngs_modulus = 3.1E10
    poissons_ratio = 0.18
  [../]
  [./stress_concrete]
    type = ComputeFiniteStrainElasticStress
    block = 1
  [../]

  # Material Properties of steel reinforcement
  [./truss]
    type = LinearElasticTruss
    block = '2 3'
    youngs_modulus = 2e11
  [../]
[]

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

[Executioner]
  solve_type = 'NEWTON'
  nl_max_its = 100
  nl_abs_tol = 1.E-5
  nl_rel_tol = 1E-3

  line_search = none

  petsc_options_iname = '-pc_type'
  petsc_options_value = 'lu'

  petsc_options = '-snes_converged_reason'

  type = Transient
  end_time = 10
  dt = 1
[]

[Outputs]
  file_base = EqualValueEmbeddedConstraint_out
  print_linear_residuals = false
  exodus = true
  csv = false
[]

Constructed Objects
Example Input Syntax
Input Parameters
References