Contact Action

Description

ContactAction can be used to specify mechanical normal and tangential contact using several possible models: - frictionless - glued - coulomb (frictional) and formulations: - kinematic - penalty (normal and tangential) - augmented lagrange - mortar (normal with/without tangential)

Example Input syntax

For normal mortar contact:

[Contact]
  [frictionless]
    mesh = simple_mesh
    master = 2
    slave = 1
    formulation = mortar
    system = constraint
  []
[]

/opt/civet/build_0/moose/modules/contact/test/tests/mechanical-small-problem/frictionless-nodal-lm-mortar-disp-action.i

[GlobalParams]
  displacements = 'disp_x disp_y'
[]

[Mesh]
  [./simple_mesh]
    type = FileMeshGenerator
    file = mesh.e
  [../]
[]

[Problem]
  kernel_coverage_check = false
[]

[Variables]
  [./disp_x]
    block = '1 2'
  [../]
  [./disp_y]
    block = '1 2'
  [../]
[]

[BCs]
  [./left_x]
    type = DirichletBC
    variable = disp_x
    boundary = 'outside_left'
    value = 0.0
  [../]
  [./left_y]
    type = DirichletBC
    variable = disp_y
    boundary = 'outside_left'
    value = 0.0
  [../]
  [./right_x]
    type = FunctionDirichletBC
    variable = disp_x
    boundary = 'outside_right'
    function = '-5e-3 * t'
  [../]
  [./right_y]
    type = FunctionDirichletBC
    variable = disp_y
    boundary = 'outside_right'
    function = 0
  [../]
[]

[Kernels]
  [disp_x]
    type = Diffusion
    variable = disp_x
    block = '1 2'
  []
  [disp_y]
    type = Diffusion
    variable = disp_y
    block = '1 2'
  []
[]

[Debug]
  show_var_residual_norms = 1
[]

[Contact]
  [frictionless]
    mesh = simple_mesh
    master = 2
    slave = 1
    formulation = mortar
    system = constraint
  []
[]

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

[Executioner]
  solve_type = PJFNK
  type = Transient
  num_steps = 10
  dt = 1
  dtmin = 1
  petsc_options_iname = '-pc_type -snes_linesearch_type -pc_factor_shift_type -pc_factor_shift_amount'
  petsc_options_value = 'lu       basic                 NONZERO               1e-15'
[]

[Outputs]
  exodus = true
  hide = 'contact_pressure nodal_area_frictionless penetration'
[]

[Postprocessors]
  [contact]
    type = ContactDOFSetSize
    variable = frictionless_normal_lm
    subdomain = '4'
    execute_on = 'nonlinear timestep_end'
  []
[]

For normal and tangential (frictional) mortar contact:

[Contact]
  [frictional]
    mesh = revised_file_mesh
    master = 20
    slave = 10
    formulation = mortar
    system = constraint
    model = coulomb
    friction_coefficient = 0.1
  []
[]

/opt/civet/build_0/moose/modules/contact/test/tests/bouncing-block-contact/frictional-nodal-min-normal-lm-mortar-fb-tangential-lm-mortar-action.i

starting_point = 2e-1

# We offset slightly so we avoid the case where the bottom of the slave block and the top of the
# master block are perfectly vertically aligned which can cause the backtracking line search some
# issues for a coarse mesh (basic line search handles that fine)
offset = 1e-2

[GlobalParams]
  displacements = 'disp_x disp_y'
  diffusivity = 1e0
  scaling = 1e0
[]

[Mesh]
  [./original_file_mesh]
    type = FileMeshGenerator
    file = long-bottom-block-1elem-blocks-coarse.e
  [../]
  # These sidesets need to be deleted because the contact action adds them automatically. For this
  # particular mesh, the new IDs will be identical to the deleted ones and will conflict if we don't
  # remove the original ones.
  [./delete_3]
    type = BlockDeletionGenerator
    input = original_file_mesh
    block_id = 3
  [../]
  [./revised_file_mesh]
    type = BlockDeletionGenerator
    input = delete_3
    block_id = 4
  [../]
[]

[Variables]
  [./disp_x]
    block = '1 2'
    # order = SECOND
  [../]
  [./disp_y]
    block = '1 2'
    # order = SECOND
  [../]
[]

[Contact]
  [frictional]
    mesh = revised_file_mesh
    master = 20
    slave = 10
    formulation = mortar
    system = constraint
    model = coulomb
    friction_coefficient = 0.1
  []
[]

[ICs]
  [./disp_y]
    block = 2
    variable = disp_y
    value = ${fparse starting_point + offset}
    type = ConstantIC
  [../]
[]

[Kernels]
  [./disp_x]
    type = MatDiffusion
    variable = disp_x
  [../]
  [./disp_y]
    type = MatDiffusion
    variable = disp_y
  [../]
[]

[BCs]
  [./botx]
    type = DirichletBC
    variable = disp_x
    boundary = 40
    value = 0.0
  [../]
  [./boty]
    type = DirichletBC
    variable = disp_y
    boundary = 40
    value = 0.0
  [../]
  [./topy]
    type = FunctionDirichletBC
    variable = disp_y
    boundary = 30
    function = '${starting_point} * cos(2 * pi / 40 * t) + ${offset}'
  [../]
  [./leftx]
    type = FunctionDirichletBC
    variable = disp_x
    boundary = 50
    function = '1e-2 * t'
  [../]
[]

[Executioner]
  type = Transient
  end_time = 200
  dt = 5
  dtmin = .1
  solve_type = 'PJFNK'
  petsc_options = '-snes_converged_reason -ksp_converged_reason -pc_svd_monitor -snes_linesearch_monitor -snes_ksp_ew'
  petsc_options_iname = '-pc_type -pc_factor_shift_type -pc_factor_shift_amount -mat_mffd_err'
  petsc_options_value = 'lu       NONZERO               1e-15                   1e-5'
  l_max_its = 30
  nl_max_its = 20
  line_search = 'none'
[]

[Debug]
  show_var_residual_norms = true
[]

[Outputs]
  exodus = true
  hide = 'contact_pressure nodal_area_frictional penetration'
[]

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

[Postprocessors]
  [./num_nl]
    type = NumNonlinearIterations
  [../]
  [./cumulative]
    type = CumulativeValuePostprocessor
    postprocessor = num_nl
  [../]
  [contact]
    type = ContactDOFSetSize
    variable = frictional_normal_lm
    subdomain = frictional_slave_subdomain
    execute_on = 'nonlinear timestep_end'
  []
[]

Input Parameters

masterThe master surface
C++ Type:BoundaryName
Options:
Description:The master surface
slaveThe slave surface
C++ Type:BoundaryName
Options:
Description:The slave surface

Required Parameters

active__all__ If specified only the blocks named will be visited and made active
Default:__all__
C++ Type:std::vector
Options:
Description:If specified only the blocks named will be visited and made active
al_frictional_force_toleranceThe tolerance of the frictional force for augmented Lagrangian method.
C++ Type:double
Options:
Description:The tolerance of the frictional force for augmented Lagrangian method.
al_incremental_slip_toleranceThe tolerance of the incremental slip for augmented Lagrangian method.
C++ Type:double
Options:
Description:The tolerance of the incremental slip for augmented Lagrangian method.
al_penetration_toleranceThe tolerance of the penetration for augmented Lagrangian method.
C++ Type:double
Options:
Description:The tolerance of the penetration for augmented Lagrangian method.
c_normal1Parameter for balancing the size of the gap and contact pressure
Default:1
C++ Type:double
Options:
Description:Parameter for balancing the size of the gap and contact pressure
c_tangential1Parameter for balancing the contact pressure and velocity
Default:1
C++ Type:double
Options:
Description:Parameter for balancing the contact pressure and velocity
capture_tolerance0Normal distance from surface within which nodes are captured
Default:0
C++ Type:double
Options:
Description:Normal distance from surface within which nodes are captured
disp_xThe x displacement
C++ Type:VariableName
Options:
Description:The x displacement
disp_yThe y displacement
C++ Type:VariableName
Options:
Description:The y displacement
disp_zThe z displacement
C++ Type:VariableName
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
formulationkinematicThe contact formulation
Default:kinematic
C++ Type:MooseEnum
Options:kinematic penalty augmented_lagrange tangential_penalty mortar
Description:The contact formulation
friction_coefficient0The friction coefficient
Default:0
C++ Type:double
Options:
Description:The friction coefficient
inactiveIf specified blocks matching these identifiers will be skipped.
C++ Type:std::vector
Options:
Description:If specified blocks matching these identifiers will be skipped.
master_slave_jacobianTrueWhether to include Jacobian entries coupling master and slave nodes.
Default:True
C++ Type:bool
Options:
Description:Whether to include Jacobian entries coupling master and slave nodes.
meshThe mesh generator for mortar method
C++ Type:MeshGeneratorName
Options:
Description:The mesh generator for mortar method
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
C++ Type:MooseEnum
Options:edge_based nodal_normal_based
Description:Method to use to smooth normals
normalize_penaltyFalseWhether to normalize the penalty parameter with the nodal area.
Default:False
C++ Type:bool
Options:
Description:Whether to normalize the penalty parameter with the nodal area.
orderFIRSTThe finite element order: FIRST, SECOND, etc.
Default:FIRST
C++ Type:MooseEnum
Options:CONSTANT FIRST SECOND THIRD FOURTH
Description:The finite element order: FIRST, SECOND, etc.
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
systemDiracKernelSystem to use for constraint enforcement
Default:DiracKernel
C++ Type:MooseEnum
Options:DiracKernel Constraint
Description:System to use for constraint enforcement
tangential_toleranceTangential distance to extend edges of contact surfaces
C++ Type:double
Options:
Description:Tangential distance to extend edges of contact surfaces
tension_release0Tension release threshold. A node in contact will not be released if its tensile load is below this value. No tension release if negative.
Default:0
C++ Type:double
Options:
Description:Tension release threshold. A node in contact will not be released if its tensile load is below this value. No tension release if negative.

Optional Parameters

Description
Example Input syntax
Input Parameters

Application Usage

Physics and Syntax

Application Development

Framework Development

MOOSEDocs

Infrastructure

Questions

Information and Tools