ComputeWeightedGapLMMechanicalContact

The Karush-Kuhn-Tucker conditions of mechanical contact are:

$var element = document.getElementById("moose-equation-b9fe48f0-695b-4dc4-9f9b-4c7a5a3a1b66");katex.render("\\begin{aligned} g &\\ge 0\\\\ \\lambda &\\ge 0\\\\ g\\lambda &= 0 \\end{aligned}", element, {displayMode:true,throwOnError:false});$

where $var element = document.getElementById("moose-equation-bc17f55c-b97e-49cc-9532-afadbf6804c4");katex.render("g", element, {displayMode:false,throwOnError:false});$ is the gap and $var element = document.getElementById("moose-equation-12cd2477-07a6-40c2-a0f1-3abb4fe1ad8e");katex.render("\\lambda", element, {displayMode:false,throwOnError:false});$ is the contact pressure, a Lagrange multipler variable living on the secondary face. Per (Wohlmuth, 2011) and (Popp and Wall, 2014), the variationally consistent, discretized version of the KKT conditions are:

$var element = document.getElementById("moose-equation-b5d8ea21-a0e2-43a6-b32c-e7da12cc575d");katex.render("\\begin{aligned} (\\tilde{g}_n)_j &\\ge 0\\\\ (\\lambda_n)_j &\\ge 0\\\\ (\\tilde{g}_n)_j (\\lambda_n)_j &= 0 \\end{aligned}", element, {displayMode:true,throwOnError:false});$

where $var element = document.getElementById("moose-equation-5e7f2871-f8d7-452d-942a-1eeb46788dc5");katex.render("n", element, {displayMode:false,throwOnError:false});$ indicates the normal direction, $var element = document.getElementById("moose-equation-44bf9668-26b0-492d-8069-8fa7b6d955c8");katex.render("j", element, {displayMode:false,throwOnError:false});$ denotes the j'th secondary contact interface node, and $var element = document.getElementById("moose-equation-5b3b27af-a48d-4ffc-82e1-99a5761c7d88");katex.render("(\\tilde{g}_n)_j", element, {displayMode:false,throwOnError:false});$ is the discrete weighted gap, computed by:

$var element = document.getElementById("moose-equation-355bd4e5-3b8d-444c-b55d-1fa8a4f049a3");katex.render("(\\tilde{g}_n)_j = \\int_{\\gamma_c^{(1)}} \\Phi_j g_{n,h} dA", element, {displayMode:true,throwOnError:false});$

where $var element = document.getElementById("moose-equation-18d15f1c-66d5-4d9a-9522-b485388bac83");katex.render("\\gamma_c^{(1)}", element, {displayMode:false,throwOnError:false});$ denotes the secondary contact interface, $var element = document.getElementById("moose-equation-fa53f76b-32fd-4278-8cb1-662578c3beff");katex.render("\\Phi_j", element, {displayMode:false,throwOnError:false});$ is the j'th lagrange multiplier test function, and $var element = document.getElementById("moose-equation-a10afe29-3ce3-481e-8267-04a43f74f917");katex.render("g_{n,h}", element, {displayMode:false,throwOnError:false});$ is the discretized version of the gap function.

The ComputeWeightedGapLMMechanicalContact object computes the weighted gap. It does not apply** the KKT conditions. That is done with the ApplyPenetrationConstraintLMMechanicalContact object. Consequently, the two objects must always be used in conjunction.

Computes the weighted gap that will later be used to enforce the zero-penetration mechanical contact conditions

Input Parameters

disp_xThe x displacement variable
C++ Type:std::vector<VariableName>
Options:
Description:The x displacement variable
disp_yThe y displacement variable
C++ Type:std::vector<VariableName>
Options:
Description:The y displacement variable
primary_boundaryThe name of the primary boundary sideset.
C++ Type:BoundaryName
Options:
Description:The name of the primary boundary sideset.
primary_subdomainThe name of the primary subdomain.
C++ Type:SubdomainName
Options:
Description:The name of the primary subdomain.
secondary_boundaryThe name of the secondary boundary sideset.
C++ Type:BoundaryName
Options:
Description:The name of the secondary boundary sideset.
secondary_subdomainThe name of the secondary subdomain.
C++ Type:SubdomainName
Options:
Description:The name of the secondary subdomain.

Required Parameters

compute_lm_residualsTrueWhether to compute Lagrange Multiplier residuals
Default:True
C++ Type:bool
Options:
Description:Whether to compute Lagrange Multiplier residuals
compute_primal_residualsTrueWhether to compute residuals for the primal variable.
Default:True
C++ Type:bool
Options:
Description:Whether to compute residuals for the primal variable.
periodicFalseWhether this constraint is going to be used to enforce a periodic condition. This has the effect of changing the normals vector for projection from outward to inward facing
Default:False
C++ Type:bool
Options:
Description:Whether this constraint is going to be used to enforce a periodic condition. This has the effect of changing the normals vector for projection from outward to inward facing
variableThe name of the lagrange multiplier variable that this constraint is applied to. This parameter may not be supplied in the case of using penalty methods for example
C++ Type:NonlinearVariableName
Options:
Description:The name of the lagrange multiplier variable that this constraint is applied to. This parameter may not be supplied in the case of using penalty methods for example

Optional Parameters

control_tagsAdds user-defined labels for accessing object parameters via control logic.
C++ Type:std::vector<std::string>
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.
implicitTrueDetermines whether this object is calculated using an implicit or explicit form
Default:True
C++ Type:bool
Options:
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
Options:
Description:The seed for the master random number generator
use_displaced_meshFalseWhether 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:False
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<TagName>
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<TagName>
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

References

Alexander Popp and WA Wall. Dual mortar methods for computational contact mechanics–overview and recent developments. GAMM-Mitteilungen, 37(1):66–84, 2014.

@article{popp2014dual,
    author = "Popp, Alexander and Wall, WA",
    title = "Dual mortar methods for computational contact mechanics--overview and recent developments",
    journal = "GAMM-Mitteilungen",
    volume = "37",
    number = "1",
    pages = "66--84",
    year = "2014",
    publisher = "Wiley Online Library"
}

Barbara Wohlmuth. Variationally consistent discretization schemes and numerical algorithms for contact problems. Acta Numerica, 20:569–734, 2011.

@article{wohlmuth2011variationally,
    author = "Wohlmuth, Barbara",
    title = "Variationally consistent discretization schemes and numerical algorithms for contact problems",
    journal = "Acta Numerica",
    volume = "20",
    pages = "569--734",
    year = "2011"
}

Input Parameters
References