MultiDContactConstraint

under construction:Undocumented Class

The MultiDContactConstraint 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.


# MultiDContactConstraint

!syntax description /Constraints/MultiDContactConstraint

## Overview

!! Replace these lines with information regarding the MultiDContactConstraint object.

## Example Input File Syntax

!! Describe and include an example of how to use the MultiDContactConstraint object.

!syntax parameters /Constraints/MultiDContactConstraint

!syntax inputs /Constraints/MultiDContactConstraint

!syntax children /Constraints/MultiDContactConstraint

!syntax description /Constraints/MultiDContactConstraint

Input Parameters

componentAn integer corresponding to the direction the variable this kernel acts in. (0 for x, 1 for y, 2 for z)
C++ Type:unsigned int
Options:
Description:An integer corresponding to the direction the variable this kernel acts in. (0 for x, 1 for y, 2 for z)
masterThe boundary ID associated with the master side
C++ Type:BoundaryName
Options:
Description:The boundary ID associated with the master side
master_variableThe variable on the master side of the domain
C++ Type:std::vector
Options:
Description:The variable on the master side of the domain
slaveThe boundary ID associated with the slave side
C++ Type:BoundaryName
Options:
Description:The boundary ID associated with the slave side
variableThe name of the variable that this constraint is applied to.
C++ Type:NonlinearVariableName
Options:
Description:The name of the variable that this constraint is applied to.

Required Parameters

disp_xThe x displacement
C++ Type:std::vector
Options:
Description:The x displacement
disp_yThe y displacement
C++ Type:std::vector
Options:
Description:The y displacement
disp_zThe z displacement
C++ Type:std::vector
Options:
Description:The z displacement
displacementsThe displacements appropriate for the simulation geometry and coordinate system
C++ Type:std::vector
Options:
Description:The displacements appropriate for the simulation geometry and coordinate system
execute_onLINEARThe list of flag(s) indicating when this object should be executed, the available options include NONE, INITIAL, LINEAR, NONLINEAR, TIMESTEP_END, TIMESTEP_BEGIN, FINAL, CUSTOM.
Default:LINEAR
C++ Type:ExecFlagEnum
Options:NONE INITIAL LINEAR NONLINEAR TIMESTEP_END TIMESTEP_BEGIN FINAL CUSTOM
Description:The list of flag(s) indicating when this object should be executed, the available options include NONE, INITIAL, LINEAR, NONLINEAR, TIMESTEP_END, TIMESTEP_BEGIN, FINAL, CUSTOM.
jacobian_updateFalseWhether or not to update the 'in contact' list every jacobian evaluation (by default it will happen once per timestep
Default:False
C++ Type:bool
Options:
Description:Whether or not to update the 'in contact' list every jacobian evaluation (by default it will happen once per timestep
modelfrictionlessThe contact model to use
Default:frictionless
C++ Type:MooseEnum
Options:frictionless glued coulomb
Description:The contact model to use
normal_smoothing_distanceDistance from edge in parametric coordinates over which to smooth contact normal
C++ Type:double
Options:
Description:Distance from edge in parametric coordinates over which to smooth contact normal
normal_smoothing_methodMethod to use to smooth normals (edge_based|nodal_normal_based)
C++ Type:std::string
Options:
Description:Method to use to smooth normals (edge_based|nodal_normal_based)
orderFIRSTThe finite element order used for projections
Default:FIRST
C++ Type:MooseEnum
Options:FIRST SECOND THIRD FOURTH
Description:The finite element order used for projections
penalty1e+08The penalty to apply. This can vary depending on the stiffness of your materials
Default:1e+08
C++ Type:double
Options:
Description:The penalty to apply. This can vary depending on the stiffness of your materials
tangential_toleranceTangential distance to extend edges of contact surfaces
C++ Type:double
Options:
Description:Tangential distance to extend edges of contact surfaces

Optional Parameters

control_tagsAdds user-defined labels for accessing object parameters via control logic.
C++ Type:std::vector
Options:
Description:Adds user-defined labels for accessing object parameters via control logic.
enableTrueSet the enabled status of the MooseObject.
Default:True
C++ Type:bool
Options:
Description:Set the enabled status of the MooseObject.
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
Options:
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

extra_matrix_tagsThe extra tags for the matrices this Kernel should fill
C++ Type:std::vector
Options:
Description:The extra tags for the matrices this Kernel should fill
extra_vector_tagsThe extra tags for the vectors this Kernel should fill
C++ Type:std::vector
Options:
Description:The extra tags for the vectors this Kernel should fill
matrix_tagssystemThe tag for the matrices this Kernel should fill
Default:system
C++ Type:MultiMooseEnum
Options:nontime system
Description:The tag for the matrices this Kernel should fill
vector_tagsnontimeThe tag for the vectors this Kernel should fill
Default:nontime
C++ Type:MultiMooseEnum
Options:nontime time
Description:The tag for the vectors this Kernel should fill

Tagging Parameters

Input Files

modules/combined/test/tests/simplest_contact/sm/simplest_contact_test_sm.i
modules/combined/test/tests/simplest_contact/simplest_contact_test.i
modules/combined/test/tests/simplest_contact/simplest_contact_skew_test.i
modules/combined/test/tests/simplest_contact/sm/simplest_contact_skew_test_sm.i

modules/combined/test/tests/simplest_contact/sm/simplest_contact_test_sm.i

[Mesh]
  file = simplest_contact.e
[]

[GlobalParams]
  volumetric_locking_correction = false
  displacements = 'disp_x disp_y'
[]

[Variables]
  [./disp_x]
  [../]
  [./disp_y]
    order = FIRST
    family = LAGRANGE
  [../]
[]

[AuxVariables]
  [./penetration]
  [../]
[]

[SolidMechanics]
  [./solid]
    disp_x = disp_x
    disp_y = disp_y
  [../]
[]

[AuxKernels]
  [./penetration]
    type = PenetrationAux
    variable = penetration
    boundary = 2
    paired_boundary = 3
  [../]
[]

[Constraints]
  [./contact_x]
    type = MultiDContactConstraint
    variable = disp_x
    master_variable = disp_x
    component = 0
    slave = 2
    master = 3
    penalty = 1e6
  [../]

  [./contact_y]
    type = MultiDContactConstraint
    variable = disp_y
    master_variable = disp_y
    component = 1
    slave = 2
    master = 3
    penalty = 1e6
  [../]
[]

[BCs]
  [./left_x]
    type = DirichletBC
    variable = disp_x
    boundary = 1
    value = 0.0
  [../]

  [./left_y]
    type = DirichletBC
    variable = disp_y
    boundary = 1
    value = 0.0
  [../]

  [./right_x]
    type = DirichletBC
    variable = disp_x
    boundary = 4
    value = -0.0001
  [../]

  [./right_y]
    type = DirichletBC
    variable = disp_y
    boundary = 4
    value = 0.0
  [../]
[] # BCs

[Materials]
  [./stiffStuff1]
    type = Elastic
    formulation = NonlinearPlaneStrain
    block = 1

    disp_x = disp_x
    disp_y = disp_y

    youngs_modulus = 1e6
    poissons_ratio = 0.3
  [../]
[] # Materials


[Executioner]
  type = Transient

  #Preconditioned JFNK (default)
  solve_type = 'PJFNK'



  petsc_options_iname = '-pc_type -pc_hypre_type -ksp_gmres_restart'
  petsc_options_value = 'hypre    boomeramg      101'


  line_search = 'none'


#  petsc_options_iname = '-pc_type -pc_hypre_type -ksp_gmres_restart'
#  petsc_options_value = 'hypre    boomeramg      101'

  nl_abs_tol = 1e-8

  l_max_its = 100
  nl_max_its = 10
  dt = 1.0
  num_steps = 1
[] # Executioner

[Outputs]
  [./out]
    type = Exodus
    elemental_as_nodal = true
  [../]
[] # Output

modules/combined/test/tests/simplest_contact/simplest_contact_test.i

[Mesh]
  file = simplest_contact.e
[]

[GlobalParams]
  volumetric_locking_correction = false
  displacements = 'disp_x disp_y'
[]

[AuxVariables]
  [./penetration]
  [../]
[]

[Modules/TensorMechanics/Master]
  [./all]
    add_variables = true
    strain = FINITE
  []
[]

[AuxKernels]
  [./penetration]
    type = PenetrationAux
    variable = penetration
    boundary = 2
    paired_boundary = 3
  [../]
[]

[Constraints]
  [./contact_x]
    type = MultiDContactConstraint
    variable = disp_x
    master_variable = disp_x
    component = 0
    slave = 2
    master = 3
    penalty = 1e6
  [../]

  [./contact_y]
    type = MultiDContactConstraint
    variable = disp_y
    master_variable = disp_y
    component = 1
    slave = 2
    master = 3
    penalty = 1e6
  [../]
[]

[BCs]
  [./left_x]
    type = DirichletBC
    variable = disp_x
    boundary = 1
    value = 0.0
  [../]

  [./left_y]
    type = DirichletBC
    variable = disp_y
    boundary = 1
    value = 0.0
  [../]

  [./right_x]
    type = DirichletBC
    variable = disp_x
    boundary = 4
    value = -0.0001
  [../]

  [./right_y]
    type = DirichletBC
    variable = disp_y
    boundary = 4
    value = 0.0
  [../]
[] # BCs

[Materials]
  [./stiffStuff1]
    type = ComputeIsotropicElasticityTensor
    block = '1'
    youngs_modulus = 1e6
    poissons_ratio = 0.3
  [../]
  [./stiffStuff1_stress]
    type = ComputeFiniteStrainElasticStress
    block = '1'
  [../]
[]  # Materials

[Executioner]
  type = Transient

  #Preconditioned JFNK (default)
  solve_type = 'PJFNK'

  petsc_options_iname = '-pc_type -pc_hypre_type -ksp_gmres_restart'
  petsc_options_value = 'hypre    boomeramg      101'


  line_search = 'none'
  nl_abs_tol = 1e-8

  l_max_its = 100
  nl_max_its = 10
  dt = 1.0
  num_steps = 1
[] # Executioner

[Outputs]
  [./out]
    type = Exodus
    elemental_as_nodal = true
  [../]
[] # Output

modules/combined/test/tests/simplest_contact/simplest_contact_skew_test.i

[Mesh]
  file = simplest_contact_skew.e
[]

[GlobalParams]
  volumetric_locking_correction = false
  displacements = 'disp_x disp_y'
[]

[AuxVariables]
  [./penetration]
    order = FIRST
    family = LAGRANGE
  [../]
[]

[Modules/TensorMechanics/Master]
  [./all]
    add_variables = true
    strain = FINITE
  []
[]

[AuxKernels]
  [./penetration]
    type = PenetrationAux
    variable = penetration
    boundary = 2
    paired_boundary = 3
  [../]
[]

[Constraints]
  [./contact_x]
    type = MultiDContactConstraint
    variable = disp_x
    master_variable = disp_x
    component = 0
    slave = 2
    master = 3
    penalty = 1e6
  [../]

  [./contact_y]
    type = MultiDContactConstraint
    variable = disp_y
    master_variable = disp_y
    component = 1
    slave = 2
    master = 3
    penalty = 1e6
  [../]
[]

[BCs]
  [./left_x]
    type = DirichletBC
    variable = disp_x
    boundary = 1
    value = 0.0
  [../]

  [./left_y]
    type = DirichletBC
    variable = disp_y
    boundary = 1
    value = 0.0
  [../]

  [./right_x]
    type = DirichletBC
    variable = disp_x
    boundary = 4
    value = -0.0001
  [../]

  [./right_y]
    type = DirichletBC
    variable = disp_y
    boundary = 4
    value = 0.0
  [../]
[]

[Materials]
  [./stiffStuff1]
    type = ComputeIsotropicElasticityTensor
    block = '1'
    youngs_modulus = 1e6
    poissons_ratio = 0.3
  [../]
  [./stiffStuff1_stress]
    type = ComputeFiniteStrainElasticStress
    block = '1'
  [../]
[]

[Executioner]
  type = Transient
  solve_type = 'PJFNK'

  #petsc_options_iname = '-pc_type -pc_hypre_type -ksp_gmres_restart'
  #petsc_options_value = 'hypre    boomeramg      101'
  petsc_options_iname = '-pc_type -pc_hypre_type -ksp_gmres_restart'
  petsc_options_value = 'hypre    boomeramg      101'

  line_search = 'none'

  nl_abs_tol = 1e-8

  l_max_its = 100
  nl_max_its = 10
  dt = 1.0
  num_steps = 1
[]

[Outputs]
  [./out]
    type = Exodus
    elemental_as_nodal = true
  [../]
[]

modules/combined/test/tests/simplest_contact/sm/simplest_contact_skew_test_sm.i

[Mesh]
  file = simplest_contact_skew.e
[]

[GlobalParams]
  volumetric_locking_correction = false
  displacements = 'disp_x disp_y'
[]

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

[AuxVariables]
  [./penetration]
    order = FIRST
    family = LAGRANGE
  [../]
[]

[SolidMechanics]
  [./solid]
    disp_x = disp_x
    disp_y = disp_y
  [../]
[]

[AuxKernels]
  [./penetration]
    type = PenetrationAux
    variable = penetration
    boundary = 2
    paired_boundary = 3
  [../]
[]

[Constraints]
  [./contact_x]
    type = MultiDContactConstraint
    variable = disp_x
    master_variable = disp_x
    component = 0
    slave = 2
    master = 3
    penalty = 1e6
  [../]

  [./contact_y]
    type = MultiDContactConstraint
    variable = disp_y
    master_variable = disp_y
    component = 1
    slave = 2
    master = 3
    penalty = 1e6
  [../]
[]

[BCs]
  [./left_x]
    type = DirichletBC
    variable = disp_x
    boundary = 1
    value = 0.0
  [../]

  [./left_y]
    type = DirichletBC
    variable = disp_y
    boundary = 1
    value = 0.0
  [../]

  [./right_x]
    type = DirichletBC
    variable = disp_x
    boundary = 4
    value = -0.0001
  [../]

  [./right_y]
    type = DirichletBC
    variable = disp_y
    boundary = 4
    value = 0.0
  [../]
[]

[Materials]
  [./stiffStuff1]
    type = Elastic
    formulation = NonlinearPlaneStrain
    block = 1

    disp_x = disp_x
    disp_y = disp_y

    youngs_modulus = 1e6
    poissons_ratio = 0.3
  [../]
[]


[Executioner]
  type = Transient
  solve_type = 'PJFNK'

  #petsc_options_iname = '-pc_type -pc_hypre_type -ksp_gmres_restart'
  #petsc_options_value = 'hypre    boomeramg      101'
  petsc_options_iname = '-pc_type -pc_hypre_type -ksp_gmres_restart'
  petsc_options_value = 'hypre    boomeramg      101'

  line_search = 'none'

  nl_abs_tol = 1e-8

  l_max_its = 100
  nl_max_its = 10
  dt = 1.0
  num_steps = 1
[]

[Outputs]
  [./out]
    type = Exodus
    elemental_as_nodal = true
  [../]
[]

Input Parameters
Input Files

Application Usage

Physics and Syntax

Application Development

Framework Development

MOOSEDocs

Infrastructure

Questions

Information and Tools