- boundaryThe primary boundary
C++ Type:BoundaryName
Controllable:No
Description:The primary boundary
- componentAn integer corresponding to the direction the variable this constraint acts on. (0 for x, 1 for y, 2 for z)
C++ Type:unsigned int
Controllable:No
Description:An integer corresponding to the direction the variable this constraint acts on. (0 for x, 1 for y, 2 for z)
- nodal_areaThe nodal area
C++ Type:std::vector<VariableName>
Controllable:No
Description:The nodal area
- variableThe name of the variable that this residual object operates on
C++ Type:NonlinearVariableName
Controllable:No
Description:The name of the variable that this residual object operates on
MechanicalContactConstraint
Apply non-penetration constraints on the mechanical deformation using a node on face, primary/secondary algorithm, and multiple options for the physical behavior on the interface and the mathematical formulation for constraint enforcement
Description
The MechanicalContactConstraint class enforces mechanical contact of different types. Three contact models are available. These are glued, frictionless, and coulomb. Here we have where denote the total traction, is the normal traction, represents the tangential traction, and are the normal and tangential unit vectors, respectively. The normal contact constraint is the same for these three contact types and are detailed in contact module description. The differences among these contact types lie in the constraint equations for the tangential contact, which is described in the following section.
Contact Models
The glued contact does not allow slip between the nodes once they come into contact, i.e., , where is the tangential relative velocity. Due to the no-slip condition, the glued contact is often used to model bonded connections between adjacent bodies that may have mismatching meshes in solid mechanics.
In contrast to the glued contact, the frictionless contact allows slip between the nodes that are in contact. No tangential traction is assumed, i.e., and . The frictionless contact is often utilized when the contact surface is pretty smooth thus tangential fraction along the interface can be omitted.
Different from the above two contact models, the coulomb model accounts for the friction along the interface caused by the relative motion between nodes. This is modeled using the Coulomb law, which computes a tangential traction that is proportional to the normal traction via a slip condition: where is the friction coefficient.
Contact Formulations
To enforce the contact constraints, three types of formulations are available. These are kinematic, penalty, and mortar. The kinematic and penalty formulations have been developed for some time and turn out to be robust for many mechanical contact problems. Recent efforts have been put in developing a mortar-based contact formulation. In this documentation, however, we focus on kinematic and penalty formulations. For readers who are interested in mortar-based mechanical contact, please go to mortar constraint system for more information.
The overall workflow of imposing a mechanical constraint is as follows: 1) Determine whether penetration occurs. If so, constrain displacement; 2) In the normal direction, eliminate the penetration; 3) In the tangential direction, slip must satisfy friction model, no slip is allowed for the glued model, no constraint is applied for the frictionless model, Kuhn-Tucker conditions must be satisfied in the tangential direction for the coulomb model. The mechanical contact enforcement steps for the glued and frictionless models are summarized below:
For penalty:
| Glued | Frictionless |
|---|---|
For kinematic:
| Glued | Frictionless |
|---|---|
Note that the kinematic method uses penalty parameter only to enforce the condition that the secondary node has the desired position on the primary face, so the errors are very small in the converged solution. Unlike in penalty formulations, penalty compliance is not introduced on the contact surface. Therefore, the converged solution has no error due to penalty compliance.
For the frictional contact (coulomb model), the penalty method calculates normal force and tangential predictor force using penalty stiffness and check for frictional capacity. While using the kinematic method, the deformation of surrounding material serves the role of the penalty parameter.
Contact gap offset
Gap offset can be provided to the current mechanical contact constraint. It can be either secondary_gap_offset (gap offset from secondary side) or mapped_primary_gap_offset (gap offset from primary side but mapped to secondary side). Use of these gap offset parameters treats the surfaces as if they were virtually extended (positive offset value) or narrowed (negative offset value) by the specified amount, so that the surfaces are treated as if they are closer or further away than they actually are. There is no deformation within the material in this gap offset region.
Input Parameters
- al_frictional_force_toleranceThe tolerance of the frictional force for augmented Lagrangian method.
C++ Type:double
Controllable:No
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
Controllable:No
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
Controllable:No
Description:The tolerance of the penetration for augmented Lagrangian method.
- capture_tolerance0Normal distance from surface within which nodes are captured
Default:0
C++ Type:double
Controllable:No
Description:Normal distance from surface within which nodes are captured
- connected_secondary_nodes_jacobianTrueWhether to include jacobian entries coupling nodes connected to secondary nodes.
Default:True
C++ Type:bool
Controllable:No
Description:Whether to include jacobian entries coupling nodes connected to secondary nodes.
- displacementsThe displacements appropriate for the simulation geometry and coordinate system
C++ Type:std::vector<VariableName>
Controllable:No
Description:The displacements appropriate for the simulation geometry and coordinate system
- formulationkinematicThe contact formulation
Default:kinematic
C++ Type:MooseEnum
Controllable:No
Description:The contact formulation
- friction_coefficient0The friction coefficient
Default:0
C++ Type:double
Controllable:No
Description:The friction coefficient
- mapped_primary_gap_offsetoffset to the gap distance mapped from primary side
C++ Type:std::vector<VariableName>
Controllable:No
Description:offset to the gap distance mapped from primary side
- modelfrictionlessThe contact model to use
Default:frictionless
C++ Type:MooseEnum
Controllable:No
Description:The contact model to use
- non_displacement_variables_jacobianTrueWhether to include jacobian entries coupling with variables that are not displacement variables.
Default:True
C++ Type:bool
Controllable:No
Description:Whether to include jacobian entries coupling with variables that are not displacement variables.
- normal_smoothing_distanceDistance from edge in parametric coordinates over which to smooth contact normal
C++ Type:double
Controllable:No
Description:Distance from edge in parametric coordinates over which to smooth contact normal
- normal_smoothing_methodMethod to use to smooth normals
C++ Type:MooseEnum
Controllable:No
Description:Method to use to smooth normals
- normalize_penaltyFalseWhether to normalize the penalty parameter with the nodal area for penalty contact.
Default:False
C++ Type:bool
Controllable:No
Description:Whether to normalize the penalty parameter with the nodal area for penalty contact.
- orderFIRSTThe finite element order used for projections
Default:FIRST
C++ Type:MooseEnum
Controllable:No
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
Controllable:No
Description:The penalty to apply. This can vary depending on the stiffness of your materials
- primaryThe boundary ID associated with the primary side
C++ Type:BoundaryName
Controllable:No
Description:The boundary ID associated with the primary side
- primary_secondary_jacobianTrueWhether to include jacobian entries coupling primary and secondary nodes.
Default:True
C++ Type:bool
Controllable:No
Description:Whether to include jacobian entries coupling primary and secondary nodes.
- primary_variableThe variable on the primary side of the domain
C++ Type:std::vector<VariableName>
Controllable:No
Description:The variable on the primary side of the domain
- print_contact_nodesFalseWhether to print the number of nodes in contact.
Default:False
C++ Type:bool
Controllable:No
Description:Whether to print the number of nodes in contact.
- secondaryThe secondary boundary
C++ Type:BoundaryName
Controllable:No
Description:The secondary boundary
- secondary_gap_offsetoffset to the gap distance from secondary side
C++ Type:std::vector<VariableName>
Controllable:No
Description:offset to the gap distance from secondary side
- stick_lock_iterations4294967295Number of times permitted to switch between sticking and slipping in a solution before locking node in a sticked state.
Default:4294967295
C++ Type:unsigned int
Controllable:No
Description:Number of times permitted to switch between sticking and slipping in a solution before locking node in a sticked state.
- stick_unlock_factor1.5Factor by which frictional capacity must be exceeded to permit stick-locked node to slip again.
Default:1.5
C++ Type:double
Controllable:No
Description:Factor by which frictional capacity must be exceeded to permit stick-locked node to slip again.
- tangential_toleranceTangential distance to extend edges of contact surfaces
C++ Type:double
Controllable:No
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
Controllable:No
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
- 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.
- 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
- 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
- extra_matrix_tagsThe extra tags for the matrices this Kernel should fill
C++ Type:std::vector<TagName>
Controllable:No
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<TagName>
Controllable:No
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
Controllable:No
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
Controllable:No
Description:The tag for the vectors this Kernel should fill