User object that computes tangential pressures due to friction using a penalty approach, following J.C. More...

#include <PenaltyFrictionUserObject.h>

Inheritance diagram for PenaltyFrictionUserObject:

Public Types
typedef DataFileName	DataFileParameterType

Public Member Functions
	PenaltyFrictionUserObject (const InputParameters &parameters)

virtual const ADVariableValue &	contactTangentialPressureDirOne () const override

virtual const ADVariableValue &	contactTangentialPressureDirTwo () const override

virtual void	initialize () override

virtual void	finalize () override

virtual void	reinit () override

virtual void	timestepSetup () override

virtual Real	getFrictionalContactPressure (const Node *const node, const unsigned int component) const override

virtual Real	getAccumulatedSlip (const Node *const node, const unsigned int component) const override

virtual Real	getTangentialVelocity (const Node *const node, const unsigned int component) const override

virtual Real	getDeltaTangentialLagrangeMultiplier (const Node *const node, const unsigned int component) const override

virtual bool	isAugmentedLagrangianConverged () override

virtual void	updateAugmentedLagrangianMultipliers () override

virtual const ADVariableValue &	contactPressure () const override

virtual Real	getNormalContactPressure (const Node *const node) const override

virtual Real	getNormalLagrangeMultiplier (const Node *const node) const

virtual bool	getActiveSetState (const Node *const node) const

virtual void	augmentedLagrangianSetup () override

virtual void	execute () override

virtual void	initialSetup () override

const std::unordered_map< const DofObject *, std::pair< ADReal, Real > > &	dofToWeightedGap () const
	Get the degree of freedom to weighted gap information. More...

virtual Real	getNormalGap (const Node *const) const

Real	physicalGap (const std::pair< ADReal, Real > &gap) const
	Compute physical gap from integration gap quantity. More...

ADReal	adPhysicalGap (const std::pair< ADReal, Real > &gap) const

virtual void	threadJoin (const UserObject &) override final

SubProblem &	getSubProblem () const

bool	shouldDuplicateInitialExecution () const

virtual Real	spatialValue (const Point &) const

virtual const std::vector< Point >	spatialPoints () const

void	gatherSum (T &value)

void	gatherMax (T &value)

void	gatherMin (T &value)

void	gatherProxyValueMax (T1 &proxy, T2 &value)

void	gatherProxyValueMin (T1 &proxy, T2 &value)

void	setPrimaryThreadCopy (UserObject *primary)

UserObject *	primaryThreadCopy ()

std::set< UserObjectName >	getDependObjects () const

virtual bool	needThreadedCopy () const

const std::set< std::string > &	getRequestedItems () override

const std::set< std::string > &	getSuppliedItems () override

unsigned int	systemNumber () const

virtual bool	enabled () const

std::shared_ptr< MooseObject >	getSharedPtr ()

std::shared_ptr< const MooseObject >	getSharedPtr () const

MooseApp &	getMooseApp () const

const std::string &	type () const

virtual const std::string &	name () const

std::string	typeAndName () const

std::string	errorPrefix (const std::string &error_type) const

void	callMooseError (std::string msg, const bool with_prefix) const

MooseObjectParameterName	uniqueParameterName (const std::string &parameter_name) const

const InputParameters &	parameters () const

MooseObjectName	uniqueName () const

const T &	getParam (const std::string &name) const

std::vector< std::pair< T1, T2 > >	getParam (const std::string &param1, const std::string &param2) const

const T *	queryParam (const std::string &name) const

const T &	getRenamedParam (const std::string &old_name, const std::string &new_name) const

T	getCheckedPointerParam (const std::string &name, const std::string &error_string="") const

bool	isParamValid (const std::string &name) const

bool	isParamSetByUser (const std::string &nm) const

void	paramError (const std::string &param, Args... args) const

void	paramWarning (const std::string &param, Args... args) const

void	paramInfo (const std::string &param, Args... args) const

void	connectControllableParams (const std::string &parameter, const std::string &object_type, const std::string &object_name, const std::string &object_parameter) const

void	mooseError (Args &&... args) const

void	mooseErrorNonPrefixed (Args &&... args) const

void	mooseDocumentedError (const std::string &repo_name, const unsigned int issue_num, Args &&... args) const

void	mooseWarning (Args &&... args) const

void	mooseWarningNonPrefixed (Args &&... args) const

void	mooseDeprecated (Args &&... args) const

void	mooseInfo (Args &&... args) const

std::string	getDataFileName (const std::string &param) const

std::string	getDataFileNameByName (const std::string &relative_path) const

std::string	getDataFilePath (const std::string &relative_path) const

virtual void	jacobianSetup ()

virtual void	residualSetup ()

virtual void	subdomainSetup ()

virtual void	customSetup (const ExecFlagType &)

const ExecFlagEnum &	getExecuteOnEnum () const

UserObjectName	getUserObjectName (const std::string &param_name) const

const T &	getUserObject (const std::string &param_name, bool is_dependency=true) const

const T &	getUserObjectByName (const UserObjectName &object_name, bool is_dependency=true) const

const UserObject &	getUserObjectBase (const std::string &param_name, bool is_dependency=true) const

const UserObject &	getUserObjectBaseByName (const UserObjectName &object_name, bool is_dependency=true) const

const std::vector< MooseVariableScalar *> &	getCoupledMooseScalarVars ()

const std::set< TagID > &	getScalarVariableCoupleableVectorTags () const

const std::set< TagID > &	getScalarVariableCoupleableMatrixTags () const

SubdomainID	primarySubdomain () const

SubdomainID	secondarySubdomain () const

bool	onInterface (BoundaryID primary_boundary_id, BoundaryID secondary_boundary_id) const

const GenericMaterialProperty< T, is_ad > &	getGenericMaterialProperty (const std::string &name, MaterialData &material_data, const unsigned int state=0)

const GenericMaterialProperty< T, is_ad > &	getGenericMaterialProperty (const std::string &name, const unsigned int state=0)

const GenericMaterialProperty< T, is_ad > &	getGenericMaterialProperty (const std::string &name, const unsigned int state=0)

const MaterialProperty< T > &	getMaterialProperty (const std::string &name, MaterialData &material_data, const unsigned int state=0)

const MaterialProperty< T > &	getMaterialProperty (const std::string &name, const unsigned int state=0)

const MaterialProperty< T > &	getMaterialProperty (const std::string &name, const unsigned int state=0)

const ADMaterialProperty< T > &	getADMaterialProperty (const std::string &name, MaterialData &material_data)

const ADMaterialProperty< T > &	getADMaterialProperty (const std::string &name)

const ADMaterialProperty< T > &	getADMaterialProperty (const std::string &name)

const MaterialProperty< T > &	getMaterialPropertyOld (const std::string &name, MaterialData &material_data)

const MaterialProperty< T > &	getMaterialPropertyOld (const std::string &name)

const MaterialProperty< T > &	getMaterialPropertyOld (const std::string &name)

const MaterialProperty< T > &	getMaterialPropertyOlder (const std::string &name, MaterialData &material_data)

const MaterialProperty< T > &	getMaterialPropertyOlder (const std::string &name)

const MaterialProperty< T > &	getMaterialPropertyOlder (const std::string &name)

const GenericMaterialProperty< T, is_ad > &	getGenericMaterialPropertyByName (const MaterialPropertyName &name, MaterialData &material_data, const unsigned int state)

const GenericMaterialProperty< T, is_ad > &	getGenericMaterialPropertyByName (const MaterialPropertyName &name, const unsigned int state=0)

const GenericMaterialProperty< T, is_ad > &	getGenericMaterialPropertyByName (const MaterialPropertyName &name, const unsigned int state=0)

const MaterialProperty< T > &	getMaterialPropertyByName (const MaterialPropertyName &name, MaterialData &material_data, const unsigned int state=0)

const MaterialProperty< T > &	getMaterialPropertyByName (const MaterialPropertyName &name, const unsigned int state=0)

const MaterialProperty< T > &	getMaterialPropertyByName (const MaterialPropertyName &name, const unsigned int state=0)

const ADMaterialProperty< T > &	getADMaterialPropertyByName (const MaterialPropertyName &name, MaterialData &material_data)

const ADMaterialProperty< T > &	getADMaterialPropertyByName (const MaterialPropertyName &name)

const ADMaterialProperty< T > &	getADMaterialPropertyByName (const MaterialPropertyName &name)

const MaterialProperty< T > &	getMaterialPropertyOldByName (const MaterialPropertyName &name, MaterialData &material_data)

const MaterialProperty< T > &	getMaterialPropertyOldByName (const MaterialPropertyName &name)

const MaterialProperty< T > &	getMaterialPropertyOldByName (const MaterialPropertyName &name)

const MaterialProperty< T > &	getMaterialPropertyOlderByName (const MaterialPropertyName &name, MaterialData &material_data)

const MaterialProperty< T > &	getMaterialPropertyOlderByName (const MaterialPropertyName &name)

const MaterialProperty< T > &	getMaterialPropertyOlderByName (const MaterialPropertyName &name)

std::pair< const MaterialProperty< T > *, std::set< SubdomainID > >	getBlockMaterialProperty (const MaterialPropertyName &name)

const GenericMaterialProperty< T, is_ad > &	getGenericZeroMaterialProperty (const std::string &name)

const GenericMaterialProperty< T, is_ad > &	getGenericZeroMaterialProperty ()

const GenericMaterialProperty< T, is_ad > &	getGenericZeroMaterialPropertyByName (const std::string &prop_name)

const MaterialProperty< T > &	getZeroMaterialProperty (Ts... args)

std::set< SubdomainID >	getMaterialPropertyBlocks (const std::string &name)

std::vector< SubdomainName >	getMaterialPropertyBlockNames (const std::string &name)

std::set< BoundaryID >	getMaterialPropertyBoundaryIDs (const std::string &name)

std::vector< BoundaryName >	getMaterialPropertyBoundaryNames (const std::string &name)

void	checkBlockAndBoundaryCompatibility (std::shared_ptr< MaterialBase > discrete)

std::unordered_map< SubdomainID, std::vector< MaterialBase *> >	buildRequiredMaterials (bool allow_stateful=true)

void	statefulPropertiesAllowed (bool)

bool	getMaterialPropertyCalled () const

virtual const std::unordered_set< unsigned int > &	getMatPropDependencies () const

virtual void	resolveOptionalProperties ()

const GenericMaterialProperty< T, is_ad > &	getPossiblyConstantGenericMaterialPropertyByName (const MaterialPropertyName &prop_name, MaterialData &material_data, const unsigned int state)

virtual const VariableValue &	coupledNeighborValue (const std::string &var_name, unsigned int comp=0) const

std::vector< const VariableValue *>	coupledNeighborValues (const std::string &var_name) const

std::vector< const VariableValue *>	coupledNeighborValuesOld (const std::string &var_name) const

std::vector< const VariableValue *>	coupledNeighborValuesOlder (const std::string &var_name) const

virtual const ADVariableValue &	adCoupledNeighborValue (const std::string &var_name, unsigned int comp=0) const

const auto &	coupledGenericNeighborValue (const std::string &var_name, unsigned int comp=0) const

const auto &	coupledGenericNeighborGradient (const std::string &var_name, unsigned int comp=0) const

virtual const ADVariableValue &	adCoupledNeighborValueDot (const std::string &var_name, unsigned int comp=0) const

std::vector< const ADVariableValue *>	adCoupledNeighborValues (const std::string &var_name) const

virtual const ADVectorVariableValue &	adCoupledVectorNeighborValue (const std::string &var_name, unsigned int comp=0) const

virtual const VariableValue &	coupledNeighborValueDot (const std::string &var_name, unsigned int comp=0) const

virtual const VariableValue &	coupledNeighborValueDotDu (const std::string &var_name, unsigned int comp=0) const

virtual const VariableValue &	coupledNeighborValueOld (const std::string &var_name, unsigned int comp=0) const

virtual const VariableValue &	coupledNeighborValueOlder (const std::string &var_name, unsigned int comp=0) const

virtual const VariableGradient &	coupledNeighborGradient (const std::string &var_name, unsigned int comp=0) const

virtual std::vector< const VariableGradient *>	coupledNeighborGradients (const std::string &var_name) const

virtual const VariableGradient &	coupledNeighborGradientOld (const std::string &var_name, unsigned int comp=0) const

virtual const VariableGradient &	coupledNeighborGradientOlder (const std::string &var_name, unsigned int comp=0) const

virtual const ADVariableGradient &	adCoupledNeighborGradient (const std::string &var_name, unsigned int comp=0) const

virtual const VectorVariableGradient &	coupledVectorNeighborGradient (const std::string &var_name, unsigned int comp=0) const

virtual const VectorVariableGradient &	coupledVectorNeighborGradientOld (const std::string &var_name, unsigned int comp=0) const

virtual const VectorVariableGradient &	coupledVectorNeighborGradientOlder (const std::string &var_name, unsigned int comp=0) const

virtual const ArrayVariableValue &	coupledArrayNeighborValue (const std::string &var_name, unsigned int comp=0) const

virtual const ArrayVariableGradient &	coupledArrayNeighborGradient (const std::string &var_name, unsigned int comp=0) const

virtual const ArrayVariableGradient &	coupledArrayNeighborGradientOld (const std::string &var_name, unsigned int comp=0) const

virtual const ArrayVariableGradient &	coupledArrayNeighborGradientOlder (const std::string &var_name, unsigned int comp=0) const

virtual const VariableSecond &	coupledNeighborSecond (const std::string &var_name, unsigned int i=0) const

virtual const VariableValue &	coupledNeighborDofValues (const std::string &var_name, unsigned int comp=0) const

virtual const VariableValue &	coupledNeighborDofValuesOld (const std::string &var_name, unsigned int comp=0) const

virtual const VariableValue &	coupledNeighborDofValuesOlder (const std::string &var_name, unsigned int comp=0) const

const std::unordered_map< std::string, std::vector< MooseVariableFieldBase *> > &	getCoupledVars () const

const std::vector< MooseVariableFieldBase *> &	getCoupledMooseVars () const

const std::vector< MooseVariable *> &	getCoupledStandardMooseVars () const

const std::vector< VectorMooseVariable *> &	getCoupledVectorMooseVars () const

const std::vector< ArrayMooseVariable *> &	getCoupledArrayMooseVars () const

void	addFEVariableCoupleableVectorTag (TagID tag)

void	addFEVariableCoupleableMatrixTag (TagID tag)

std::set< TagID > &	getFEVariableCoupleableVectorTags ()

const std::set< TagID > &	getFEVariableCoupleableVectorTags () const

std::set< TagID > &	getFEVariableCoupleableMatrixTags ()

const std::set< TagID > &	getFEVariableCoupleableMatrixTags () const

auto &	getWritableCoupledVariables () const

bool	hasWritableCoupledVariables () const

const ADVariableValue *	getADDefaultValue (const std::string &var_name) const

const ADVectorVariableValue *	getADDefaultVectorValue (const std::string &var_name) const

const ADVariableGradient &	getADDefaultGradient () const

const ADVectorVariableGradient &	getADDefaultVectorGradient () const

const ADVariableSecond &	getADDefaultSecond () const

const ADVectorVariableCurl &	getADDefaultCurl () const

bool	hasUserObject (const std::string &param_name) const

bool	hasUserObject (const std::string &param_name) const

bool	hasUserObject (const std::string &param_name) const

bool	hasUserObject (const std::string &param_name) const

bool	hasUserObjectByName (const UserObjectName &object_name) const

bool	hasUserObjectByName (const UserObjectName &object_name) const

bool	hasUserObjectByName (const UserObjectName &object_name) const

bool	hasUserObjectByName (const UserObjectName &object_name) const

const GenericMaterialProperty< T, is_ad > &	getGenericNeighborMaterialProperty (const std::string &name, const unsigned int state=0)

const GenericMaterialProperty< T, is_ad > &	getGenericNeighborMaterialProperty (const std::string &name, const unsigned int state=0)

const MaterialProperty< T > &	getNeighborMaterialProperty (const std::string &name, const unsigned int state=0)

const MaterialProperty< T > &	getNeighborMaterialProperty (const std::string &name, const unsigned int state=0)

const ADMaterialProperty< T > &	getNeighborADMaterialProperty (const std::string &name)

const ADMaterialProperty< T > &	getNeighborADMaterialProperty (const std::string &name)

const MaterialProperty< T > &	getNeighborMaterialPropertyOld (const std::string &name)

const MaterialProperty< T > &	getNeighborMaterialPropertyOld (const std::string &name)

const MaterialProperty< T > &	getNeighborMaterialPropertyOlder (const std::string &name)

const MaterialProperty< T > &	getNeighborMaterialPropertyOlder (const std::string &name)

const GenericMaterialProperty< T, is_ad > &	getGenericNeighborMaterialPropertyByName (const std::string &name, const unsigned int state=0)

const GenericMaterialProperty< T, is_ad > &	getGenericNeighborMaterialPropertyByName (const std::string &name, const unsigned int state=0)

const MaterialProperty< T > &	getNeighborMaterialPropertyByName (const std::string &name, const unsigned int state=0)

const MaterialProperty< T > &	getNeighborMaterialPropertyByName (const std::string &name, const unsigned int state=0)

const ADMaterialProperty< T > &	getNeighborADMaterialPropertyByName (const std::string &name)

const ADMaterialProperty< T > &	getNeighborADMaterialPropertyByName (const std::string &name)

const GenericOptionalMaterialProperty< T, is_ad > &	getGenericOptionalMaterialProperty (const std::string &name, const unsigned int state=0)

const GenericOptionalMaterialProperty< T, is_ad > &	getGenericOptionalMaterialProperty (const std::string &name, const unsigned int state=0)

const OptionalMaterialProperty< T > &	getOptionalMaterialProperty (const std::string &name, const unsigned int state=0)

const OptionalMaterialProperty< T > &	getOptionalMaterialProperty (const std::string &name, const unsigned int state=0)

const OptionalADMaterialProperty< T > &	getOptionalADMaterialProperty (const std::string &name)

const OptionalADMaterialProperty< T > &	getOptionalADMaterialProperty (const std::string &name)

const OptionalMaterialProperty< T > &	getOptionalMaterialPropertyOld (const std::string &name)

const OptionalMaterialProperty< T > &	getOptionalMaterialPropertyOld (const std::string &name)

const OptionalMaterialProperty< T > &	getOptionalMaterialPropertyOlder (const std::string &name)

const OptionalMaterialProperty< T > &	getOptionalMaterialPropertyOlder (const std::string &name)

MaterialBase &	getMaterial (const std::string &name)

MaterialBase &	getMaterial (const std::string &name)

MaterialBase &	getMaterialByName (const std::string &name, bool no_warn=false)

MaterialBase &	getMaterialByName (const std::string &name, bool no_warn=false)

bool	hasMaterialProperty (const std::string &name)

bool	hasMaterialProperty (const std::string &name)

bool	hasMaterialPropertyByName (const std::string &name)

bool	hasMaterialPropertyByName (const std::string &name)

bool	hasADMaterialProperty (const std::string &name)

bool	hasADMaterialProperty (const std::string &name)

bool	hasADMaterialPropertyByName (const std::string &name)

bool	hasADMaterialPropertyByName (const std::string &name)

bool	hasGenericMaterialProperty (const std::string &name)

bool	hasGenericMaterialProperty (const std::string &name)

bool	hasGenericMaterialPropertyByName (const std::string &name)

bool	hasGenericMaterialPropertyByName (const std::string &name)

const Function &	getFunction (const std::string &name) const

const Function &	getFunctionByName (const FunctionName &name) const

bool	hasFunction (const std::string &param_name) const

bool	hasFunctionByName (const FunctionName &name) const

bool	isDefaultPostprocessorValue (const std::string &param_name, const unsigned int index=0) const

bool	hasPostprocessor (const std::string &param_name, const unsigned int index=0) const

bool	hasPostprocessorByName (const PostprocessorName &name) const

std::size_t	coupledPostprocessors (const std::string &param_name) const

const PostprocessorName &	getPostprocessorName (const std::string &param_name, const unsigned int index=0) const

const VectorPostprocessorValue &	getVectorPostprocessorValue (const std::string &param_name, const std::string &vector_name) const

const VectorPostprocessorValue &	getVectorPostprocessorValue (const std::string &param_name, const std::string &vector_name, bool needs_broadcast) const

const VectorPostprocessorValue &	getVectorPostprocessorValueByName (const VectorPostprocessorName &name, const std::string &vector_name) const

const VectorPostprocessorValue &	getVectorPostprocessorValueByName (const VectorPostprocessorName &name, const std::string &vector_name, bool needs_broadcast) const

const VectorPostprocessorValue &	getVectorPostprocessorValueOld (const std::string &param_name, const std::string &vector_name) const

const VectorPostprocessorValue &	getVectorPostprocessorValueOld (const std::string &param_name, const std::string &vector_name, bool needs_broadcast) const

const VectorPostprocessorValue &	getVectorPostprocessorValueOldByName (const VectorPostprocessorName &name, const std::string &vector_name) const

const VectorPostprocessorValue &	getVectorPostprocessorValueOldByName (const VectorPostprocessorName &name, const std::string &vector_name, bool needs_broadcast) const

const ScatterVectorPostprocessorValue &	getScatterVectorPostprocessorValue (const std::string &param_name, const std::string &vector_name) const

const ScatterVectorPostprocessorValue &	getScatterVectorPostprocessorValueByName (const VectorPostprocessorName &name, const std::string &vector_name) const

const ScatterVectorPostprocessorValue &	getScatterVectorPostprocessorValueOld (const std::string &param_name, const std::string &vector_name) const

const ScatterVectorPostprocessorValue &	getScatterVectorPostprocessorValueOldByName (const VectorPostprocessorName &name, const std::string &vector_name) const

bool	hasVectorPostprocessor (const std::string &param_name, const std::string &vector_name) const

bool	hasVectorPostprocessor (const std::string &param_name) const

bool	hasVectorPostprocessorByName (const VectorPostprocessorName &name, const std::string &vector_name) const

bool	hasVectorPostprocessorByName (const VectorPostprocessorName &name) const

const VectorPostprocessorName &	getVectorPostprocessorName (const std::string &param_name) const

T &	getSampler (const std::string &name)

Sampler &	getSampler (const std::string &name)

T &	getSamplerByName (const SamplerName &name)

Sampler &	getSamplerByName (const SamplerName &name)

virtual void	meshChanged ()

virtual void	meshDisplaced ()

PerfGraph &	perfGraph ()

const PostprocessorValue &	getPostprocessorValue (const std::string &param_name, const unsigned int index=0) const

const PostprocessorValue &	getPostprocessorValue (const std::string &param_name, const unsigned int index=0) const

const PostprocessorValue &	getPostprocessorValueOld (const std::string &param_name, const unsigned int index=0) const

const PostprocessorValue &	getPostprocessorValueOld (const std::string &param_name, const unsigned int index=0) const

const PostprocessorValue &	getPostprocessorValueOlder (const std::string &param_name, const unsigned int index=0) const

const PostprocessorValue &	getPostprocessorValueOlder (const std::string &param_name, const unsigned int index=0) const

virtual const PostprocessorValue &	getPostprocessorValueByName (const PostprocessorName &name) const

virtual const PostprocessorValue &	getPostprocessorValueByName (const PostprocessorName &name) const

const PostprocessorValue &	getPostprocessorValueOldByName (const PostprocessorName &name) const

const PostprocessorValue &	getPostprocessorValueOldByName (const PostprocessorName &name) const

const PostprocessorValue &	getPostprocessorValueOlderByName (const PostprocessorName &name) const

const PostprocessorValue &	getPostprocessorValueOlderByName (const PostprocessorName &name) const

bool	isVectorPostprocessorDistributed (const std::string &param_name) const

bool	isVectorPostprocessorDistributed (const std::string &param_name) const

bool	isVectorPostprocessorDistributedByName (const VectorPostprocessorName &name) const

bool	isVectorPostprocessorDistributedByName (const VectorPostprocessorName &name) const

const Distribution &	getDistribution (const std::string &name) const

const T &	getDistribution (const std::string &name) const

const Distribution &	getDistribution (const std::string &name) const

const T &	getDistribution (const std::string &name) const

const Distribution &	getDistributionByName (const DistributionName &name) const

const T &	getDistributionByName (const std::string &name) const

const Distribution &	getDistributionByName (const DistributionName &name) const

const T &	getDistributionByName (const std::string &name) const

const Parallel::Communicator &	comm () const

processor_id_type	n_processors () const

processor_id_type	processor_id () const

const std::unordered_map< const DofObject *, std::array< ADReal, 2 > > &	dofToWeightedVelocities () const
	Get the degree of freedom to weighted velocities information. More...

Static Public Member Functions
static InputParameters	validParams ()

static void	sort (typename std::vector< T > &vector)

static void	sortDFS (typename std::vector< T > &vector)

static void	cyclicDependencyError (CyclicDependencyException< T2 > &e, const std::string &header)

Public Attributes
const ConsoleStream	_console

Static Public Attributes
static constexpr PropertyValue::id_type	default_property_id

static constexpr PropertyValue::id_type	zero_property_id

static constexpr auto	SYSTEM

static constexpr auto	NAME

Protected Types
enum	AdaptivityFrictionalPenalty { AdaptivityFrictionalPenalty::SIMPLE, AdaptivityFrictionalPenalty::FRICTION_LIMIT }
	The adaptivity method for the penalty factor at augmentations. More...

enum	AdaptivityNormalPenalty { AdaptivityNormalPenalty::SIMPLE, AdaptivityNormalPenalty::BUSSETTA }
	The adaptivity method for the penalty factor at augmentations. More...

Protected Member Functions
virtual const VariableTestValue &	test () const override

virtual bool	constrainedByOwner () const override

void	selfInitialize ()

void	selfFinalize ()

void	selfTimestepSetup ()

void	bussettaAdaptivePenalty (const Real previous_gap, const Real gap, Real &penalty, Real &eval_tn)
	Adaptive, local penalty for AL. More...

void	adaptiveNormalPenalty (const Real previous_gap, const Real gap, Real &penalty)
	See Algorithm 3 of 'The adapted augmented Lagrangian method: a new method for the resolution of the mechanical frictional contact problem', Comput Mech (2012) 49: 259-275. More...

virtual void	computeQpProperties ()
	Computes properties that are functions only of the current quadrature point (`_qp`), e.g. More...

virtual void	computeQpIProperties ()
	Computes properties that are functions both of `_qp` and `_i`, for example the weighted gap. More...

template<typename K , typename V >
V	findValue (const std::unordered_map< K, V > &map, const K &key, const V &default_value=0) const
	Find a value in a map or return a default if the key doesn't exist. More...

virtual void	addPostprocessorDependencyHelper (const PostprocessorName &name) const override

virtual void	addVectorPostprocessorDependencyHelper (const VectorPostprocessorName &name) const override

virtual void	addUserObjectDependencyHelper (const UserObject &uo) const override

void	addReporterDependencyHelper (const ReporterName &reporter_name) override

const ReporterName &	getReporterName (const std::string &param_name) const

T &	declareRestartableData (const std::string &data_name, Args &&... args)

ManagedValue< T >	declareManagedRestartableDataWithContext (const std::string &data_name, void *context, Args &&... args)

const T &	getRestartableData (const std::string &data_name) const

T &	declareRestartableDataWithContext (const std::string &data_name, void *context, Args &&... args)

T &	declareRecoverableData (const std::string &data_name, Args &&... args)

T &	declareRestartableDataWithObjectName (const std::string &data_name, const std::string &object_name, Args &&... args)

T &	declareRestartableDataWithObjectNameWithContext (const std::string &data_name, const std::string &object_name, void *context, Args &&... args)

std::string	restartableName (const std::string &data_name) const

const T &	getMeshProperty (const std::string &data_name, const std::string &prefix)

const T &	getMeshProperty (const std::string &data_name)

bool	hasMeshProperty (const std::string &data_name, const std::string &prefix) const

bool	hasMeshProperty (const std::string &data_name, const std::string &prefix) const

bool	hasMeshProperty (const std::string &data_name) const

bool	hasMeshProperty (const std::string &data_name) const

std::string	meshPropertyName (const std::string &data_name) const

PerfID	registerTimedSection (const std::string &section_name, const unsigned int level) const

PerfID	registerTimedSection (const std::string &section_name, const unsigned int level, const std::string &live_message, const bool print_dots=true) const

std::string	timedSectionName (const std::string &section_name) const

bool	isCoupledScalar (const std::string &var_name, unsigned int i=0) const

unsigned int	coupledScalarComponents (const std::string &var_name) const

unsigned int	coupledScalar (const std::string &var_name, unsigned int comp=0) const

libMesh::Order	coupledScalarOrder (const std::string &var_name, unsigned int comp=0) const

const VariableValue &	coupledScalarValue (const std::string &var_name, unsigned int comp=0) const

const ADVariableValue &	adCoupledScalarValue (const std::string &var_name, unsigned int comp=0) const

const GenericVariableValue< is_ad > &	coupledGenericScalarValue (const std::string &var_name, unsigned int comp=0) const

const GenericVariableValue< false > &	coupledGenericScalarValue (const std::string &var_name, const unsigned int comp) const

const GenericVariableValue< true > &	coupledGenericScalarValue (const std::string &var_name, const unsigned int comp) const

const VariableValue &	coupledVectorTagScalarValue (const std::string &var_name, TagID tag, unsigned int comp=0) const

const VariableValue &	coupledMatrixTagScalarValue (const std::string &var_name, TagID tag, unsigned int comp=0) const

const VariableValue &	coupledScalarValueOld (const std::string &var_name, unsigned int comp=0) const

const VariableValue &	coupledScalarValueOlder (const std::string &var_name, unsigned int comp=0) const

const VariableValue &	coupledScalarDot (const std::string &var_name, unsigned int comp=0) const

const ADVariableValue &	adCoupledScalarDot (const std::string &var_name, unsigned int comp=0) const

const VariableValue &	coupledScalarDotDot (const std::string &var_name, unsigned int comp=0) const

const VariableValue &	coupledScalarDotOld (const std::string &var_name, unsigned int comp=0) const

const VariableValue &	coupledScalarDotDotOld (const std::string &var_name, unsigned int comp=0) const

const VariableValue &	coupledScalarDotDu (const std::string &var_name, unsigned int comp=0) const

const VariableValue &	coupledScalarDotDotDu (const std::string &var_name, unsigned int comp=0) const

const MooseVariableScalar *	getScalarVar (const std::string &var_name, unsigned int comp) const

const std::set< SubdomainID > &	getHigherDimSubdomainIDs () const

const std::set< BoundaryID > &	getBoundaryIDs () const

const AutomaticMortarGeneration &	amg () const

bool	interpolateNormals () const

virtual void	checkMaterialProperty (const std::string &name, const unsigned int state)

void	markMatPropRequested (const std::string &)

MaterialPropertyName	getMaterialPropertyName (const std::string &name) const

void	checkExecutionStage ()

virtual void	coupledCallback (const std::string &, bool) const

virtual bool	isCoupled (const std::string &var_name, unsigned int i=0) const

virtual bool	isCoupledConstant (const std::string &var_name) const

unsigned int	coupledComponents (const std::string &var_name) const

VariableName	coupledName (const std::string &var_name, unsigned int comp=0) const

std::vector< VariableName >	coupledNames (const std::string &var_name) const

virtual unsigned int	coupled (const std::string &var_name, unsigned int comp=0) const

std::vector< unsigned int >	coupledIndices (const std::string &var_name) const

virtual const VariableValue &	coupledValue (const std::string &var_name, unsigned int comp=0) const

std::vector< const VariableValue *>	coupledValues (const std::string &var_name) const

std::vector< const VectorVariableValue *>	coupledVectorValues (const std::string &var_name) const

const GenericVariableValue< is_ad > &	coupledGenericValue (const std::string &var_name, unsigned int comp=0) const

const GenericVariableValue< false > &	coupledGenericValue (const std::string &var_name, unsigned int comp) const

const GenericVariableValue< true > &	coupledGenericValue (const std::string &var_name, unsigned int comp) const

const GenericVectorVariableValue< is_ad > &	coupledGenericVectorValue (const std::string &var_name, unsigned int comp=0) const

const GenericVectorVariableValue< false > &	coupledGenericVectorValue (const std::string &var_name, unsigned int comp) const

const GenericVectorVariableValue< true > &	coupledGenericVectorValue (const std::string &var_name, unsigned int comp) const

std::vector< const GenericVariableValue< is_ad > *>	coupledGenericValues (const std::string &var_name) const

std::vector< const GenericVariableValue< false > *>	coupledGenericValues (const std::string &var_name) const

std::vector< const GenericVariableValue< true > *>	coupledGenericValues (const std::string &var_name) const

const GenericVariableValue< is_ad > &	coupledGenericDofValue (const std::string &var_name, unsigned int comp=0) const

const GenericVariableValue< false > &	coupledGenericDofValue (const std::string &var_name, unsigned int comp) const

const GenericVariableValue< true > &	coupledGenericDofValue (const std::string &var_name, unsigned int comp) const

const GenericVariableValue< is_ad > &	coupledGenericDot (const std::string &var_name, unsigned int comp=0) const

const GenericVariableValue< false > &	coupledGenericDot (const std::string &var_name, unsigned int comp) const

const GenericVariableValue< true > &	coupledGenericDot (const std::string &var_name, unsigned int comp) const

const GenericVariableValue< is_ad > &	coupledGenericDotDot (const std::string &var_name, unsigned int comp=0) const

const GenericVariableValue< false > &	coupledGenericDotDot (const std::string &var_name, unsigned int comp) const

const GenericVariableValue< true > &	coupledGenericDotDot (const std::string &var_name, unsigned int comp) const

virtual const VariableValue &	coupledValueLower (const std::string &var_name, unsigned int comp=0) const

const ADVariableValue &	adCoupledValue (const std::string &var_name, unsigned int comp=0) const

std::vector< const ADVariableValue *>	adCoupledValues (const std::string &var_name) const

const ADVariableValue &	adCoupledLowerValue (const std::string &var_name, unsigned int comp=0) const

const ADVectorVariableValue &	adCoupledVectorValue (const std::string &var_name, unsigned int comp=0) const

std::vector< const ADVectorVariableValue *>	adCoupledVectorValues (const std::string &var_name) const

virtual const VariableValue &	coupledVectorTagValue (const std::string &var_names, TagID tag, unsigned int index=0) const

virtual const VariableValue &	coupledVectorTagValue (const std::string &var_names, const std::string &tag_name, unsigned int index=0) const

std::vector< const VariableValue *>	coupledVectorTagValues (const std::string &var_names, TagID tag) const

std::vector< const VariableValue *>	coupledVectorTagValues (const std::string &var_names, const std::string &tag_name) const

virtual const ArrayVariableValue &	coupledVectorTagArrayValue (const std::string &var_names, TagID tag, unsigned int index=0) const

virtual const ArrayVariableValue &	coupledVectorTagArrayValue (const std::string &var_names, const std::string &tag_name, unsigned int index=0) const

std::vector< const ArrayVariableValue *>	coupledVectorTagArrayValues (const std::string &var_names, TagID tag) const

std::vector< const ArrayVariableValue *>	coupledVectorTagArrayValues (const std::string &var_names, const std::string &tag_name) const

virtual const VariableGradient &	coupledVectorTagGradient (const std::string &var_names, TagID tag, unsigned int index=0) const

virtual const VariableGradient &	coupledVectorTagGradient (const std::string &var_names, const std::string &tag_name, unsigned int index=0) const

std::vector< const VariableGradient *>	coupledVectorTagGradients (const std::string &var_names, TagID tag) const

std::vector< const VariableGradient *>	coupledVectorTagGradients (const std::string &var_names, const std::string &tag_name) const

virtual const ArrayVariableGradient &	coupledVectorTagArrayGradient (const std::string &var_names, TagID tag, unsigned int index=0) const

virtual const ArrayVariableGradient &	coupledVectorTagArrayGradient (const std::string &var_names, const std::string &tag_name, unsigned int index=0) const

std::vector< const ArrayVariableGradient *>	coupledVectorTagArrayGradients (const std::string &var_names, TagID tag) const

std::vector< const ArrayVariableGradient *>	coupledVectorTagArrayGradients (const std::string &var_names, const std::string &tag_name) const

virtual const VariableValue &	coupledVectorTagDofValue (const std::string &var_name, TagID tag, unsigned int index=0) const

virtual const VariableValue &	coupledVectorTagDofValue (const std::string &var_names, const std::string &tag_name, unsigned int index=0) const

const ArrayVariableValue &	coupledVectorTagArrayDofValue (const std::string &var_name, const std::string &tag_name, unsigned int comp=0) const

std::vector< const VariableValue *>	coupledVectorTagDofValues (const std::string &var_names, TagID tag) const

std::vector< const VariableValue *>	coupledVectorTagDofValues (const std::string &var_names, const std::string &tag_name) const

virtual const VariableValue &	coupledMatrixTagValue (const std::string &var_names, TagID tag, unsigned int index=0) const

virtual const VariableValue &	coupledMatrixTagValue (const std::string &var_names, const std::string &tag_name, unsigned int index=0) const

std::vector< const VariableValue *>	coupledMatrixTagValues (const std::string &var_names, TagID tag) const

std::vector< const VariableValue *>	coupledMatrixTagValues (const std::string &var_names, const std::string &tag_name) const

virtual const VectorVariableValue &	coupledVectorValue (const std::string &var_name, unsigned int comp=0) const

virtual const ArrayVariableValue &	coupledArrayValue (const std::string &var_name, unsigned int comp=0) const

std::vector< const ArrayVariableValue *>	coupledArrayValues (const std::string &var_name) const

MooseWritableVariable &	writableVariable (const std::string &var_name, unsigned int comp=0)

virtual VariableValue &	writableCoupledValue (const std::string &var_name, unsigned int comp=0)

void	checkWritableVar (MooseWritableVariable *var)

virtual const VariableValue &	coupledValueOld (const std::string &var_name, unsigned int comp=0) const

std::vector< const VariableValue *>	coupledValuesOld (const std::string &var_name) const

std::vector< const VectorVariableValue *>	coupledVectorValuesOld (const std::string &var_name) const

virtual const VariableValue &	coupledValueOlder (const std::string &var_name, unsigned int comp=0) const

std::vector< const VariableValue *>	coupledValuesOlder (const std::string &var_name) const

virtual const VariableValue &	coupledValuePreviousNL (const std::string &var_name, unsigned int comp=0) const

virtual const VectorVariableValue &	coupledVectorValueOld (const std::string &var_name, unsigned int comp=0) const

virtual const VectorVariableValue &	coupledVectorValueOlder (const std::string &var_name, unsigned int comp=0) const

virtual const ArrayVariableValue &	coupledArrayValueOld (const std::string &var_name, unsigned int comp=0) const

virtual const ArrayVariableValue &	coupledArrayValueOlder (const std::string &var_name, unsigned int comp=0) const

virtual const VariableGradient &	coupledGradient (const std::string &var_name, unsigned int comp=0) const

std::vector< const VariableGradient *>	coupledGradients (const std::string &var_name) const

const ADVariableGradient &	adCoupledGradient (const std::string &var_name, unsigned int comp=0) const

const ADVariableGradient &	adCoupledGradientDot (const std::string &var_name, unsigned int comp=0) const

std::vector< const ADVariableGradient *>	adCoupledGradients (const std::string &var_name) const

const GenericVariableGradient< is_ad > &	coupledGenericGradient (const std::string &var_name, unsigned int comp=0) const

const GenericVariableGradient< false > &	coupledGenericGradient (const std::string &var_name, unsigned int comp) const

const GenericVariableGradient< true > &	coupledGenericGradient (const std::string &var_name, unsigned int comp) const

std::vector< const GenericVariableGradient< is_ad > *>	coupledGenericGradients (const std::string &var_name) const

std::vector< const GenericVariableGradient< false > *>	coupledGenericGradients (const std::string &var_name) const

std::vector< const GenericVariableGradient< true > *>	coupledGenericGradients (const std::string &var_name) const

const ADVectorVariableGradient &	adCoupledVectorGradient (const std::string &var_name, unsigned int comp=0) const

const ADVariableSecond &	adCoupledSecond (const std::string &var_name, unsigned int comp=0) const

const ADVectorVariableSecond &	adCoupledVectorSecond (const std::string &var_name, unsigned int comp=0) const

virtual const VariableGradient &	coupledGradientOld (const std::string &var_name, unsigned int comp=0) const

std::vector< const VariableGradient *>	coupledGradientsOld (const std::string &var_name) const

virtual const VariableGradient &	coupledGradientOlder (const std::string &var_name, unsigned int comp=0) const

virtual const VariableGradient &	coupledGradientPreviousNL (const std::string &var_name, unsigned int comp=0) const

virtual const VariableGradient &	coupledGradientDot (const std::string &var_name, unsigned int comp=0) const

virtual const VariableGradient &	coupledGradientDotDot (const std::string &var_name, unsigned int comp=0) const

virtual const VectorVariableGradient &	coupledVectorGradient (const std::string &var_name, unsigned int comp=0) const

virtual const VectorVariableGradient &	coupledVectorGradientOld (const std::string &var_name, unsigned int comp=0) const

virtual const VectorVariableGradient &	coupledVectorGradientOlder (const std::string &var_name, unsigned int comp=0) const

virtual const ArrayVariableGradient &	coupledArrayGradient (const std::string &var_name, unsigned int comp=0) const

virtual const ArrayVariableGradient &	coupledArrayGradientOld (const std::string &var_name, unsigned int comp=0) const

virtual const ArrayVariableGradient &	coupledArrayGradientOlder (const std::string &var_name, unsigned int comp=0) const

virtual const ArrayVariableGradient &	coupledArrayGradientDot (const std::string &var_name, unsigned int comp=0) const

virtual const VectorVariableCurl &	coupledCurl (const std::string &var_name, unsigned int comp=0) const

virtual const VectorVariableCurl &	coupledCurlOld (const std::string &var_name, unsigned int comp=0) const

virtual const VectorVariableCurl &	coupledCurlOlder (const std::string &var_name, unsigned int comp=0) const

const ADVectorVariableCurl &	adCoupledCurl (const std::string &var_name, unsigned int comp=0) const

virtual const VectorVariableDivergence &	coupledDiv (const std::string &var_name, unsigned int comp=0) const

virtual const VectorVariableDivergence &	coupledDivOld (const std::string &var_name, unsigned int comp=0) const

virtual const VectorVariableDivergence &	coupledDivOlder (const std::string &var_name, unsigned int comp=0) const

virtual const VariableSecond &	coupledSecond (const std::string &var_name, unsigned int comp=0) const

virtual const VariableSecond &	coupledSecondOld (const std::string &var_name, unsigned int comp=0) const

virtual const VariableSecond &	coupledSecondOlder (const std::string &var_name, unsigned int comp=0) const

virtual const VariableSecond &	coupledSecondPreviousNL (const std::string &var_name, unsigned int comp=0) const

virtual const VariableValue &	coupledDot (const std::string &var_name, unsigned int comp=0) const

std::vector< const VariableValue *>	coupledDots (const std::string &var_name) const

virtual const VariableValue &	coupledDotDot (const std::string &var_name, unsigned int comp=0) const

virtual const VariableValue &	coupledDotOld (const std::string &var_name, unsigned int comp=0) const

virtual const VariableValue &	coupledDotDotOld (const std::string &var_name, unsigned int comp=0) const

const ADVariableValue &	adCoupledDot (const std::string &var_name, unsigned int comp=0) const

std::vector< const ADVariableValue *>	adCoupledDots (const std::string &var_name) const

const ADVariableValue &	adCoupledDotDot (const std::string &var_name, unsigned int comp=0) const

const ADVectorVariableValue &	adCoupledVectorDot (const std::string &var_name, unsigned int comp=0) const

virtual const VectorVariableValue &	coupledVectorDot (const std::string &var_name, unsigned int comp=0) const

virtual const VectorVariableValue &	coupledVectorDotDot (const std::string &var_name, unsigned int comp=0) const

virtual const VectorVariableValue &	coupledVectorDotOld (const std::string &var_name, unsigned int comp=0) const

virtual const VectorVariableValue &	coupledVectorDotDotOld (const std::string &var_name, unsigned int comp=0) const

virtual const VariableValue &	coupledVectorDotDu (const std::string &var_name, unsigned int comp=0) const

virtual const VariableValue &	coupledVectorDotDotDu (const std::string &var_name, unsigned int comp=0) const

virtual const ArrayVariableValue &	coupledArrayDot (const std::string &var_name, unsigned int comp=0) const

virtual const ArrayVariableValue &	coupledArrayDotDot (const std::string &var_name, unsigned int comp=0) const

virtual const ArrayVariableValue &	coupledArrayDotOld (const std::string &var_name, unsigned int comp=0) const

virtual const ArrayVariableValue &	coupledArrayDotDotOld (const std::string &var_name, unsigned int comp=0) const

virtual const VariableValue &	coupledDotDu (const std::string &var_name, unsigned int comp=0) const

virtual const VariableValue &	coupledDotDotDu (const std::string &var_name, unsigned int comp=0) const

const VariableValue &	coupledArrayDotDu (const std::string &var_name, unsigned int comp=0) const

const T &	coupledNodalValue (const std::string &var_name, unsigned int comp=0) const

const Moose::ADType< T >::type &	adCoupledNodalValue (const std::string &var_name, unsigned int comp=0) const

const T &	coupledNodalValueOld (const std::string &var_name, unsigned int comp=0) const

const T &	coupledNodalValueOlder (const std::string &var_name, unsigned int comp=0) const

const T &	coupledNodalValuePreviousNL (const std::string &var_name, unsigned int comp=0) const

const T &	coupledNodalDot (const std::string &var_name, unsigned int comp=0) const

virtual const VariableValue &	coupledNodalDotDot (const std::string &var_name, unsigned int comp=0) const

virtual const VariableValue &	coupledNodalDotOld (const std::string &var_name, unsigned int comp=0) const

virtual const VariableValue &	coupledNodalDotDotOld (const std::string &var_name, unsigned int comp=0) const

virtual const VariableValue &	coupledDofValues (const std::string &var_name, unsigned int comp=0) const

std::vector< const VariableValue *>	coupledAllDofValues (const std::string &var_name) const

virtual const VariableValue &	coupledDofValuesOld (const std::string &var_name, unsigned int comp=0) const

std::vector< const VariableValue *>	coupledAllDofValuesOld (const std::string &var_name) const

virtual const VariableValue &	coupledDofValuesOlder (const std::string &var_name, unsigned int comp=0) const

std::vector< const VariableValue *>	coupledAllDofValuesOlder (const std::string &var_name) const

virtual const ArrayVariableValue &	coupledArrayDofValues (const std::string &var_name, unsigned int comp=0) const

virtual const ADVariableValue &	adCoupledDofValues (const std::string &var_name, unsigned int comp=0) const

const ADVariableValue &	adZeroValue () const

const ADVariableGradient &	adZeroGradient () const

const ADVariableSecond &	adZeroSecond () const

const GenericVariableValue< is_ad > &	genericZeroValue ()

const GenericVariableValue< false > &	genericZeroValue ()

const GenericVariableValue< true > &	genericZeroValue ()

const GenericVariableGradient< is_ad > &	genericZeroGradient ()

const GenericVariableGradient< false > &	genericZeroGradient ()

const GenericVariableGradient< true > &	genericZeroGradient ()

const GenericVariableSecond< is_ad > &	genericZeroSecond ()

const GenericVariableSecond< false > &	genericZeroSecond ()

const GenericVariableSecond< true > &	genericZeroSecond ()

bool	checkVar (const std::string &var_name, unsigned int comp=0, unsigned int comp_bound=0) const

const MooseVariableFieldBase *	getFEVar (const std::string &var_name, unsigned int comp) const

const MooseVariableFieldBase *	getFieldVar (const std::string &var_name, unsigned int comp) const

MooseVariableFieldBase *	getFieldVar (const std::string &var_name, unsigned int comp)

const T *	getVarHelper (const std::string &var_name, unsigned int comp) const

T *	getVarHelper (const std::string &var_name, unsigned int comp)

MooseVariable *	getVar (const std::string &var_name, unsigned int comp)

const MooseVariable *	getVar (const std::string &var_name, unsigned int comp) const

VectorMooseVariable *	getVectorVar (const std::string &var_name, unsigned int comp)

const VectorMooseVariable *	getVectorVar (const std::string &var_name, unsigned int comp) const

ArrayMooseVariable *	getArrayVar (const std::string &var_name, unsigned int comp)

const ArrayMooseVariable *	getArrayVar (const std::string &var_name, unsigned int comp) const

void	validateExecutionerType (const std::string &name, const std::string &fn_name) const

std::vector< T >	coupledVectorHelper (const std::string &var_name, const Func &func) const

const T &	getReporterValue (const std::string &param_name, const std::size_t time_index=0)

const T &	getReporterValue (const std::string &param_name, ReporterMode mode, const std::size_t time_index=0)

const T &	getReporterValue (const std::string &param_name, const std::size_t time_index=0)

const T &	getReporterValue (const std::string &param_name, ReporterMode mode, const std::size_t time_index=0)

const T &	getReporterValueByName (const ReporterName &reporter_name, const std::size_t time_index=0)

const T &	getReporterValueByName (const ReporterName &reporter_name, ReporterMode mode, const std::size_t time_index=0)

const T &	getReporterValueByName (const ReporterName &reporter_name, const std::size_t time_index=0)

const T &	getReporterValueByName (const ReporterName &reporter_name, ReporterMode mode, const std::size_t time_index=0)

bool	hasReporterValue (const std::string &param_name) const

bool	hasReporterValue (const std::string &param_name) const

bool	hasReporterValue (const std::string &param_name) const

bool	hasReporterValue (const std::string &param_name) const

bool	hasReporterValueByName (const ReporterName &reporter_name) const

bool	hasReporterValueByName (const ReporterName &reporter_name) const

bool	hasReporterValueByName (const ReporterName &reporter_name) const

bool	hasReporterValueByName (const ReporterName &reporter_name) const

const GenericMaterialProperty< T, is_ad > *	defaultGenericMaterialProperty (const std::string &name)

const GenericMaterialProperty< T, is_ad > *	defaultGenericMaterialProperty (const std::string &name)

const MaterialProperty< T > *	defaultMaterialProperty (const std::string &name)

const MaterialProperty< T > *	defaultMaterialProperty (const std::string &name)

const ADMaterialProperty< T > *	defaultADMaterialProperty (const std::string &name)

const ADMaterialProperty< T > *	defaultADMaterialProperty (const std::string &name)

virtual void	computeQpProperties () override
	Computes properties that are functions only of the current quadrature point (`_qp`), e.g. More...

virtual void	computeQpIProperties () override
	Computes properties that are functions both of `_qp` and `_i`, for example the weighted gap. More...

void	selfInitialize ()

void	selfTimestepSetup ()

Static Protected Member Functions
static std::string	meshPropertyName (const std::string &data_name, const std::string &prefix)

static void	trimDerivative (dof_id_type remove_derivative_index, ADReal &dual_number)

static void	trimInteriorNodeDerivatives (const std::map< unsigned int, unsigned int > &primary_ip_lowerd_map, const Variables &moose_var, DualNumbers &ad_vars, const bool is_secondary)

Protected Attributes
const Real	_penalty
	The normal penalty factor. More...

const Real	_penalty_friction
	The penalty factor for the frictional constraints. More...

const Real	_slip_tolerance
	Acceptable slip distance for augmented Lagrange convergence. More...

const Real	_friction_coefficient
	The friction coefficient. More...

std::unordered_map< const DofObject *, std::pair< TwoVector, TwoVector > >	_dof_to_step_slip
	Map from degree of freedom to current and old step slip. More...

std::unordered_map< const DofObject *, std::pair< TwoVector, TwoVector > >	_dof_to_accumulated_slip
	Map from degree of freedom to current and old accumulated slip. More...

std::unordered_map< const DofObject *, std::pair< ADTwoVector, TwoVector > >	_dof_to_tangential_traction
	Map from degree of freedom to current and old tangential traction. More...

ADVariableValue	_frictional_contact_traction_one
	The first frictional contact pressure on the mortar segment quadrature points. More...

ADVariableValue	_frictional_contact_traction_two
	The second frictional contact pressure on the mortar segment quadrature points. More...

std::unordered_map< const DofObject *, TwoVector >	_dof_to_frictional_lagrange_multipliers
	Map from degree of freedom to augmented lagrange multiplier. More...

std::unordered_map< const DofObject *, Real >	_dof_to_local_penalty_friction
	Map from degree of freedom to local friction penalty value. More...

const Real	_penalty_multiplier_friction
	Penalty growth factor for augmented Lagrange. More...

enum PenaltyFrictionUserObject::AdaptivityFrictionalPenalty	_adaptivity_friction

const Real	_epsilon_tolerance
	Tolerance to avoid NaN/Inf in automatic differentiation operations. More...

const Real	_penalty_multiplier
	penalty growth factor for augmented Lagrange More...

const Real	_penetration_tolerance
	penetration tolerance for augmented Lagrange contact More...

ADVariableValue	_contact_pressure
	The contact pressure on the mortar segument quadrature points. More...

std::unordered_map< const DofObject *, ADReal >	_dof_to_normal_pressure
	Map from degree of freedom to normal pressure for reporting. More...

std::unordered_map< const DofObject *, Real >	_dof_to_lagrange_multiplier
	Map from degree of freedom to augmented lagrange multiplier. More...

std::unordered_map< const DofObject *, Real >	_dof_to_local_penalty
	Map from degree of freedom to local penalty value. More...

std::unordered_map< const DofObject *, Real >	_dof_to_previous_gap
	Map from degree of freedom to previous AL iteration gap values. More...

const Real &	_dt
	Current delta t... or timestep size. More...

Real	_dt_old
	previous timestep size More...

const bool	_use_physical_gap
	Use scaled or physical gap. More...

const MooseVariable *const	_aux_lm_var
	The auxiliary Lagrange multiplier variable (used together whith the Petrov-Galerkin approach) More...

const Real	_max_penalty_multiplier
	Maximum multiplier applied to the initial penalty factor in AL. More...

enum PenaltyWeightedGapUserObject::AdaptivityNormalPenalty	_adaptivity_normal

FEProblemBase &	_fe_problem
	The base finite element problem. More...

ADReal	_qp_gap
	The value of the gap at the current quadrature point. More...

Real	_qp_factor
	The value of the LM at the current quadrature point. More...

const bool	_nodal
	Whether the dof objects are nodal; if they're not, then they're elemental. More...

const MooseVariable *const	_disp_x_var
	The x displacement variable. More...

const MooseVariable *const	_disp_y_var
	The y displacement variable. More...

const bool	_has_disp_z
	For 2D mortar contact no displacement will be specified, so const pointers used. More...

const MooseVariable *const	_disp_z_var
	The z displacement variable. More...

const ADVariableValue &	_secondary_disp_x
	x-displacement on the secondary face More...

const ADVariableValue &	_primary_disp_x
	x-displacement on the primary face More...

const ADVariableValue &	_secondary_disp_y
	y-displacement on the secondary face More...

const ADVariableValue &	_primary_disp_y
	y-displacement on the primary face More...

const ADVariableValue *const	_secondary_disp_z
	z-displacement on the secondary face More...

const ADVariableValue *const	_primary_disp_z
	z-displacement on the primary face More...

const MooseArray< Real > &	_coord
	Member for handling change of coordinate systems (xyz, rz, spherical) More...

ADRealVectorValue	_qp_gap_nodal
	Vector for computation of weighted gap with nodal normals. More...

ADRealVectorValue	_qp_displacement_nodal
	Vector for computation of relative displacement (determines mixity ratio in interface problems) More...

std::unordered_map< const DofObject *, std::pair< ADReal, Real > >	_dof_to_weighted_gap
	A map from node to weighted gap and normalization (if requested) More...

std::unordered_map< const DofObject *, ADRealVectorValue >	_dof_to_weighted_displacements
	A map from node to weighted displacements. More...

const ADReal *	_weighted_gap_ptr = nullptr
	A pointer members that can be used to help avoid copying ADReals. More...

const Real *	_normalization_ptr = nullptr

const VariableTestValue *	_test = nullptr
	A pointer to the test function associated with the weighted gap. More...

bool	_is_weighted_gap_nodal = true
	Whether the weighted gap is associated with nodes or elements (like for a CONSTANT MONOMIAL Lagrange multiplier). More...

unsigned int	_qp = 0
	Quadrature point index for the mortar segments. More...

unsigned int	_i = 0
	Test function index. More...

SubProblem &	_subproblem

SystemBase &	_sys

const THREAD_ID	_tid

Assembly &	_assembly

const Moose::CoordinateSystemType &	_coord_sys

const bool	_duplicate_initial_execution

std::set< std::string >	_depend_uo

const bool &	_enabled

MooseApp &	_app

const std::string	_type

const std::string	_name

const InputParameters &	_pars

Factory &	_factory

ActionFactory &	_action_factory

const ExecFlagEnum &	_execute_enum

const ExecFlagType &	_current_execute_flag

MooseApp &	_restartable_app

const std::string	_restartable_system_name

const THREAD_ID	_restartable_tid

const bool	_restartable_read_only

FEProblemBase &	_mci_feproblem

FEProblemBase &	_mdi_feproblem

MooseApp &	_pg_moose_app

const std::string	_prefix

FEProblemBase &	_sc_fe_problem

const THREAD_ID	_sc_tid

const Real &	_real_zero

const VariableValue &	_scalar_zero

const Point &	_point_zero

FEProblemBase &	_mci_fe_problem

SubProblem &	_mci_subproblem

const THREAD_ID	_mci_tid

MooseMesh &	_mci_mesh

Assembly &	_mci_assembly

const MortarData &	_mortar_data

const BoundaryID	_secondary_id

const BoundaryID	_primary_id

const SubdomainID	_secondary_subdomain_id

const SubdomainID	_primary_subdomain_id

std::set< BoundaryID >	_secondary_set

std::set< BoundaryID >	_boundary_ids

std::set< SubdomainID >	_higher_dim_subdomain_ids

const bool	_interpolate_normals

const MooseArray< Point > &	_phys_points_secondary

const MooseArray< Point > &	_phys_points_primary

const libMesh::QBase *const &	_qrule_msm

const libMesh::QBase *const &	_qrule_face

Elem const *const &	_lower_secondary_elem

Elem const *const &	_lower_primary_elem

const std::vector< Real > &	_JxW_msm

const Elem *const &	_msm_elem

std::vector< Point >	_normals

MaterialData &	_neighbor_material_data

const InputParameters &	_mi_params

const std::string	_mi_name

const MooseObjectName	_mi_moose_object_name

FEProblemBase &	_mi_feproblem

SubProblem &	_mi_subproblem

const THREAD_ID	_mi_tid

const Moose::MaterialDataType	_material_data_type

MaterialData &	_material_data

bool	_stateful_allowed

bool	_get_material_property_called

std::vector< std::unique_ptr< PropertyValue > >	_default_properties

std::unordered_set< unsigned int >	_material_property_dependencies

const MaterialPropertyName	_get_suffix

const bool	_use_interpolated_state

bool	_neighbor_nodal

const InputParameters &	_c_parameters

const std::string &	_c_name

const std::string &	_c_type

FEProblemBase &	_c_fe_problem

const SystemBase *const	_c_sys

std::unordered_map< std::string, std::vector< MooseVariableFieldBase *> >	_coupled_vars

std::vector< MooseVariableFieldBase *>	_coupled_moose_vars

std::vector< MooseVariable *>	_coupled_standard_moose_vars

std::vector< VectorMooseVariable *>	_coupled_vector_moose_vars

std::vector< ArrayMooseVariable *>	_coupled_array_moose_vars

std::vector< MooseVariableFV< Real > *>	_coupled_standard_fv_moose_vars

std::vector< MooseLinearVariableFV< Real > *>	_coupled_standard_linear_fv_moose_vars

const std::unordered_map< std::string, std::string > &	_new_to_deprecated_coupled_vars

bool	_c_nodal

bool	_c_is_implicit

const bool	_c_allow_element_to_nodal_coupling

THREAD_ID	_c_tid

std::unordered_map< std::string, std::vector< std::unique_ptr< VariableValue > > >	_default_value

std::unordered_map< std::string, std::unique_ptr< MooseArray< ADReal > > >	_ad_default_value

std::unordered_map< std::string, std::unique_ptr< VectorVariableValue > >	_default_vector_value

std::unordered_map< std::string, std::unique_ptr< ArrayVariableValue > >	_default_array_value

std::unordered_map< std::string, std::unique_ptr< MooseArray< ADRealVectorValue > > >	_ad_default_vector_value

VariableValue	_default_value_zero

VariableGradient	_default_gradient

MooseArray< ADRealVectorValue >	_ad_default_gradient

MooseArray< ADRealTensorValue >	_ad_default_vector_gradient

VariableSecond	_default_second

MooseArray< ADRealTensorValue >	_ad_default_second

MooseArray< ADRealVectorValue >	_ad_default_curl

const VariableValue &	_zero

const VariablePhiValue &	_phi_zero

const MooseArray< ADReal > &	_ad_zero

const VariableGradient &	_grad_zero

const MooseArray< ADRealVectorValue > &	_ad_grad_zero

const VariablePhiGradient &	_grad_phi_zero

const VariableSecond &	_second_zero

const MooseArray< ADRealTensorValue > &	_ad_second_zero

const VariablePhiSecond &	_second_phi_zero

const VectorVariableValue &	_vector_zero

const VectorVariableCurl &	_vector_curl_zero

VectorVariableValue	_default_vector_value_zero

VectorVariableGradient	_default_vector_gradient

VectorVariableCurl	_default_vector_curl

VectorVariableDivergence	_default_div

ArrayVariableValue	_default_array_value_zero

ArrayVariableGradient	_default_array_gradient

bool	_coupleable_neighbor

const Parallel::Communicator &	_communicator

std::unordered_map< const DofObject *, std::array< ADReal, 2 > >	_dof_to_weighted_tangential_velocity
	A map from node to two weighted tangential velocities. More...

std::unordered_map< const DofObject *, std::array< ADReal, 2 > >	_dof_to_real_tangential_velocity
	A map from node to two interpolated, physical tangential velocities. More...

std::array< const ADReal *, 2 >	_tangential_vel_ptr = {{nullptr, nullptr}}
	An array of two pointers to avoid copies. More...

std::array< ADReal, 2 >	_qp_tangential_velocity
	The value of the tangential velocity values at the current quadrature point. More...

std::array< ADReal, 2 >	_qp_real_tangential_velocity
	The value of the "real" tangential velocity values at the current quadrature point. More...

ADRealVectorValue	_qp_tangential_velocity_nodal
	The value of the tangential velocity vectors at the current node. More...

ADRealVectorValue	_qp_real_tangential_velocity_nodal
	The value of the real tangential velocity vectors at the current node. More...

SystemBase &	_sys
	Reference to the EquationSystem object. More...

MooseVariableField< Real > &	_secondary_var
	Reference to the secondary variable. More...

MooseVariableField< Real > &	_primary_var
	Reference to the primary variable. More...

const ADVariableValue &	_secondary_x_dot
	x-velocity on the secondary face More...

const ADVariableValue &	_primary_x_dot
	x-velocity on the primary face More...

const ADVariableValue &	_secondary_y_dot
	y-velocity on the secondary face More...

const ADVariableValue &	_primary_y_dot
	y-velocity on the primary face More...

const ADVariableValue *const	_secondary_z_dot
	z-velocity on the secondary face More...

const ADVariableValue *const	_primary_z_dot
	z-velocity on the primary face More...

bool	_3d
	Automatic flag to determine whether we are doing three-dimensional work. More...

Static Protected Attributes
static const std::string	_interpolated_old

static const std::string	_interpolated_older

AugmentedLagrangianContactProblemInterface *const	_augmented_lagrange_problem
	augmented Lagrange problem and iteration number More...

const unsigned int &	_lagrangian_iteration_number

static const unsigned int	_no_iterations = 0

Detailed Description

User object that computes tangential pressures due to friction using a penalty approach, following J.C.

Simo, T.A. Laursen, An augmented lagrangian treatment of contact problems involving friction, https://doi.org/10.1016/0045-7949(92)90540-G.

Definition at line 22 of file PenaltyFrictionUserObject.h.

Member Enumeration Documentation

◆ AdaptivityFrictionalPenalty

enum PenaltyFrictionUserObject::AdaptivityFrictionalPenalty

strongprotected

The adaptivity method for the penalty factor at augmentations.

Enumerator
SIMPLE
FRICTION_LIMIT

Definition at line 92 of file PenaltyFrictionUserObject.h.

92 { SIMPLE, FRICTION_LIMIT } _adaptivity_friction;

SIMPLE

Executioner set up to solve a thermal-hydraulics problem using the SIMPLE algorithm.

Definition: SIMPLE.h:19

PenaltyFrictionUserObject::_adaptivity_friction

enum PenaltyFrictionUserObject::AdaptivityFrictionalPenalty _adaptivity_friction

◆ AdaptivityNormalPenalty

enum PenaltyWeightedGapUserObject::AdaptivityNormalPenalty

strongprotectedinherited

The adaptivity method for the penalty factor at augmentations.

Enumerator
SIMPLE
BUSSETTA

Definition at line 121 of file PenaltyWeightedGapUserObject.h.

121 { SIMPLE, BUSSETTA } _adaptivity_normal;

SIMPLE

Executioner set up to solve a thermal-hydraulics problem using the SIMPLE algorithm.

Definition: SIMPLE.h:19

PenaltyWeightedGapUserObject::_adaptivity_normal

enum PenaltyWeightedGapUserObject::AdaptivityNormalPenalty _adaptivity_normal

Constructor & Destructor Documentation

◆ PenaltyFrictionUserObject()

PenaltyFrictionUserObject::PenaltyFrictionUserObject ( const InputParameters & parameters )

Definition at line 60 of file PenaltyFrictionUserObject.C.

   : WeightedGapUserObject(parameters),
     PenaltyWeightedGapUserObject(parameters),
     WeightedVelocitiesUserObject(parameters),
     _penalty(getParam<Real>("penalty")),
     _penalty_friction(isParamValid("penalty_friction") ? getParam<Real>("penalty_friction")
                                                        : getParam<Real>("penalty")),
     _slip_tolerance(isParamValid("slip_tolerance") ? getParam<Real>("slip_tolerance") : 0.0),
     _friction_coefficient(getParam<Real>("friction_coefficient")),
     _penalty_multiplier_friction(getParam<Real>("penalty_multiplier_friction")),
     _adaptivity_friction(
         getParam<MooseEnum>("adaptivity_penalty_friction").getEnum<AdaptivityFrictionalPenalty>()),
     _epsilon_tolerance(1.0e-40)
 
 {
   if (!_augmented_lagrange_problem == isParamValid("slip_tolerance"))
     paramError("slip_tolerance",
                "This parameter must be supplied if and only if an augmented Lagrange problem "
                "object is used.");
 }

Member Function Documentation

◆ adaptiveNormalPenalty()

void PenaltyWeightedGapUserObject::adaptiveNormalPenalty	(	const Real	previous_gap,
		const Real	gap,
		Real &	penalty
	)

protectedinherited

See Algorithm 3 of 'The adapted augmented Lagrangian method: a new method for the resolution of the mechanical frictional contact problem', Comput Mech (2012) 49: 259-275.

Definition at line 360 of file PenaltyWeightedGapUserObject.C.

Referenced by PenaltyWeightedGapUserObject::bussettaAdaptivePenalty().

 {
   const bool condition_one = std::abs(std::abs(previous_gap / gap) - 1.0) < 0.01;
   const bool condition_two = std::abs(gap) > 1.01 * std::abs(previous_gap);
 
   if (previous_gap * gap < 0)
   {
     if (previous_gap > _penetration_tolerance)
       penalty = std::abs(penalty * previous_gap / gap * (std::abs(gap) + _penetration_tolerance) /
                          (gap - previous_gap));
     else
       penalty = std::abs(penalty * previous_gap / 10.0 / gap);
   }
   else if (std::abs(gap) > _penetration_tolerance)
   {
     if (std::abs(gap - previous_gap) >
         std::max(gap / 10.0, std::max(previous_gap / 10.0, 5.0 * _penetration_tolerance)))
       penalty *= 10.0;
     else if (condition_one && gap < 10.0 * _penetration_tolerance)
       penalty *= MathUtils::pow(std::sqrt(std::abs(gap) / _penetration_tolerance - 1.0) + 1.0, 2);
     else if (condition_two)
       penalty *= 10.0 * previous_gap / gap;
     else
       penalty *= std::sqrt(std::abs(gap) / _penetration_tolerance - 1.0) + 1.0;
   }
 }

◆ adPhysicalGap()

ADReal WeightedGapUserObject::adPhysicalGap ( const std::pair< ADReal, Real > & gap ) const

inlineinherited

Definition at line 62 of file WeightedGapUserObject.h.

Referenced by PenaltyWeightedGapUserObject::selfFinalize().

62 { return gap.first / gap.second; }

◆ augmentedLagrangianSetup()

void PenaltyWeightedGapUserObject::augmentedLagrangianSetup ( )

overridevirtualinherited

Implements AugmentedLagrangeInterface.

Definition at line 283 of file PenaltyWeightedGapUserObject.C.

 {
   // Loop over all nodes for which a gap has been computed
   for (auto & [dof_object, wgap] : _dof_to_weighted_gap)
   {
     const Real gap = physicalGap(wgap);
     // Store previous augmented lagrange iteration gap
     _dof_to_previous_gap[dof_object] = gap;
   }
 }

◆ bussettaAdaptivePenalty()

void PenaltyWeightedGapUserObject::bussettaAdaptivePenalty	(	const Real	previous_gap,
		const Real	gap,
		Real &	penalty,
		Real &	eval_tn
	)

protectedinherited

Adaptive, local penalty for AL.

See 'The adapted augmented Lagrangian method: a new method for the resolution of the mechanical frictional contact problem', Comput Mech (2012) 49: 259-275

Definition at line 337 of file PenaltyWeightedGapUserObject.C.

Referenced by PenaltyWeightedGapUserObject::updateAugmentedLagrangianMultipliers().

 {
   // Positive gaps means no contact
   if (previous_gap > 0.0)
   {
     penalty = _penalty;
     eval_tn = 0.0;
   }
   else
   {
     if (previous_gap * gap < 0)
       eval_tn = 0.0;
     else
       eval_tn = penalty * previous_gap;
 
     adaptiveNormalPenalty(previous_gap, gap, penalty);
   }
 }

◆ computeQpIProperties() [1/2]

void WeightedVelocitiesUserObject::computeQpIProperties ( )

overrideprotectedvirtualinherited

Computes properties that are functions both of _qp and _i, for example the weighted gap.

Reimplemented from WeightedGapUserObject.

Reimplemented in BilinearMixedModeCohesiveZoneModel.

Definition at line 103 of file WeightedVelocitiesUserObject.C.

Referenced by BilinearMixedModeCohesiveZoneModel::computeQpIProperties().

 {
   WeightedGapUserObject::computeQpIProperties();
 
   const auto & nodal_tangents = amg().getNodalTangents(*_lower_secondary_elem);
   // Get the _dof_to_weighted_tangential_velocity map
   const DofObject * const dof =
       _is_weighted_gap_nodal ? static_cast<const DofObject *>(_lower_secondary_elem->node_ptr(_i))
                              : static_cast<const DofObject *>(_lower_secondary_elem);
 
   _dof_to_weighted_tangential_velocity[dof][0] +=
       (*_test)[_i][_qp] * _qp_tangential_velocity_nodal * nodal_tangents[0][_i];
   _dof_to_real_tangential_velocity[dof][0] +=
       (*_test)[_i][_qp] * _qp_real_tangential_velocity_nodal * nodal_tangents[0][_i];
 
   // Get the _dof_to_weighted_tangential_velocity map for a second direction
   if (_3d)
   {
     _dof_to_weighted_tangential_velocity[dof][1] +=
         (*_test)[_i][_qp] * _qp_tangential_velocity_nodal * nodal_tangents[1][_i];
 
     _dof_to_real_tangential_velocity[dof][1] +=
         (*_test)[_i][_qp] * _qp_real_tangential_velocity_nodal * nodal_tangents[1][_i];
   }
 }

◆ computeQpIProperties() [2/2]

void WeightedGapUserObject::computeQpIProperties ( )

protectedvirtualinherited

Computes properties that are functions both of _qp and _i, for example the weighted gap.

Reimplemented in WeightedVelocitiesUserObject, and BilinearMixedModeCohesiveZoneModel.

Definition at line 121 of file WeightedGapUserObject.C.

Referenced by WeightedVelocitiesUserObject::computeQpIProperties(), and WeightedGapUserObject::execute().

 {
   mooseAssert(_normals.size() == _lower_secondary_elem->n_nodes(),
               "Making sure that _normals is the expected size");
 
   // Get the _dof_to_weighted_gap map
   const auto * const dof = static_cast<const DofObject *>(_lower_secondary_elem->node_ptr(_i));
 
   auto & [weighted_gap, normalization] = _dof_to_weighted_gap[dof];
 
   weighted_gap += (*_test)[_i][_qp] * _qp_gap_nodal * _normals[_i];
   normalization += (*_test)[_i][_qp] * _qp_factor;
 
   _dof_to_weighted_displacements[dof] += (*_test)[_i][_qp] * _qp_displacement_nodal;
 }

◆ computeQpProperties() [1/2]

void WeightedVelocitiesUserObject::computeQpProperties ( )

overrideprotectedvirtualinherited

Computes properties that are functions only of the current quadrature point (_qp), e.g.

indepedent of shape functions

Reimplemented from WeightedGapUserObject.

Reimplemented in BilinearMixedModeCohesiveZoneModel.

Definition at line 57 of file WeightedVelocitiesUserObject.C.

Referenced by BilinearMixedModeCohesiveZoneModel::computeQpProperties().

 {
   // Compute the value of _qp_gap
   WeightedGapUserObject::computeQpProperties();
 
   // Trim derivatives
   const auto & primary_ip_lowerd_map = amg().getPrimaryIpToLowerElementMap(
       *_lower_primary_elem, *_lower_primary_elem->interior_parent(), *_lower_secondary_elem);
   const auto & secondary_ip_lowerd_map =
       amg().getSecondaryIpToLowerElementMap(*_lower_secondary_elem);
 
   std::array<const MooseVariable *, 3> var_array{{_disp_x_var, _disp_y_var, _disp_z_var}};
   std::array<ADReal, 3> primary_disp_dot{
       {_primary_x_dot[_qp], _primary_y_dot[_qp], _has_disp_z ? (*_primary_z_dot)[_qp] : 0}};
   std::array<ADReal, 3> secondary_disp_dot{
       {_secondary_x_dot[_qp], _secondary_y_dot[_qp], _has_disp_z ? (*_secondary_z_dot)[_qp] : 0}};
 
   trimInteriorNodeDerivatives(primary_ip_lowerd_map, var_array, primary_disp_dot, false);
   trimInteriorNodeDerivatives(secondary_ip_lowerd_map, var_array, secondary_disp_dot, true);
 
   const ADReal & prim_x_dot = primary_disp_dot[0];
   const ADReal & prim_y_dot = primary_disp_dot[1];
   const ADReal * prim_z_dot = nullptr;
   if (_has_disp_z)
     prim_z_dot = &primary_disp_dot[2];
 
   const ADReal & sec_x_dot = secondary_disp_dot[0];
   const ADReal & sec_y_dot = secondary_disp_dot[1];
   const ADReal * sec_z_dot = nullptr;
   if (_has_disp_z)
     sec_z_dot = &secondary_disp_dot[2];
 
   // Build relative velocity vector
   ADRealVectorValue relative_velocity;
 
   if (_3d)
     relative_velocity = {sec_x_dot - prim_x_dot, sec_y_dot - prim_y_dot, *sec_z_dot - *prim_z_dot};
   else
     relative_velocity = {sec_x_dot - prim_x_dot, sec_y_dot - prim_y_dot, 0.0};
 
   // Geometry is averaged and used at the nodes for constraint enforcement.
   _qp_real_tangential_velocity_nodal = relative_velocity;
   _qp_tangential_velocity_nodal = relative_velocity * (_JxW_msm[_qp] * _coord[_qp]);
 }

◆ computeQpProperties() [2/2]

void WeightedGapUserObject::computeQpProperties ( )

protectedvirtualinherited

Computes properties that are functions only of the current quadrature point (_qp), e.g.

indepedent of shape functions

Reimplemented in WeightedVelocitiesUserObject, and BilinearMixedModeCohesiveZoneModel.

Definition at line 66 of file WeightedGapUserObject.C.

Referenced by WeightedVelocitiesUserObject::computeQpProperties(), and WeightedGapUserObject::execute().

 {
   // Trim interior node variable derivatives
   const auto & primary_ip_lowerd_map = amg().getPrimaryIpToLowerElementMap(
       *_lower_primary_elem, *_lower_primary_elem->interior_parent(), *_lower_secondary_elem);
   const auto & secondary_ip_lowerd_map =
       amg().getSecondaryIpToLowerElementMap(*_lower_secondary_elem);
 
   std::array<const MooseVariable *, 3> var_array{{_disp_x_var, _disp_y_var, _disp_z_var}};
   std::array<ADReal, 3> primary_disp{
       {_primary_disp_x[_qp], _primary_disp_y[_qp], _has_disp_z ? (*_primary_disp_z)[_qp] : 0}};
   std::array<ADReal, 3> secondary_disp{{_secondary_disp_x[_qp],
                                         _secondary_disp_y[_qp],
                                         _has_disp_z ? (*_secondary_disp_z)[_qp] : 0}};
 
   trimInteriorNodeDerivatives(primary_ip_lowerd_map, var_array, primary_disp, false);
   trimInteriorNodeDerivatives(secondary_ip_lowerd_map, var_array, secondary_disp, true);
 
   const ADReal & prim_x = primary_disp[0];
   const ADReal & prim_y = primary_disp[1];
   const ADReal * prim_z = nullptr;
   if (_has_disp_z)
     prim_z = &primary_disp[2];
 
   const ADReal & sec_x = secondary_disp[0];
   const ADReal & sec_y = secondary_disp[1];
   const ADReal * sec_z = nullptr;
   if (_has_disp_z)
     sec_z = &secondary_disp[2];
 
   // Compute gap vector
   ADRealVectorValue gap_vec = _phys_points_primary[_qp] - _phys_points_secondary[_qp];
 
   // Generic displacement for interface problems
   _qp_displacement_nodal(0) = prim_x - sec_x;
   _qp_displacement_nodal(1) = prim_y - sec_y;
   if (_has_disp_z)
     _qp_displacement_nodal(2) = *prim_z - *sec_z;
 
   _qp_displacement_nodal *= _JxW_msm[_qp] * _coord[_qp];
 
   gap_vec(0).derivatives() = prim_x.derivatives() - sec_x.derivatives();
   gap_vec(1).derivatives() = prim_y.derivatives() - sec_y.derivatives();
   if (_has_disp_z)
     gap_vec(2).derivatives() = prim_z->derivatives() - sec_z->derivatives();
 
   // Compute integration point quantities: Normals (geometry) is averaged at the node, but not
   // interpolated within the weak integration.
   _qp_gap_nodal = gap_vec * (_JxW_msm[_qp] * _coord[_qp]);
 
   // To do normalization of constraint coefficient (c_n)
   _qp_factor = _JxW_msm[_qp] * _coord[_qp];
 }

◆ constrainedByOwner()

virtual bool PenaltyFrictionUserObject::constrainedByOwner ( ) const

inlineoverrideprotectedvirtual

Returns: Whether the gap constraint will be enforced solely by the owner of the weighted gap or will be enforced in a distributed way (like in a penalty method)

Reimplemented from PenaltyWeightedGapUserObject.

Definition at line 52 of file PenaltyFrictionUserObject.h.

52 { return false; }

◆ contactPressure()

const ADVariableValue & PenaltyWeightedGapUserObject::contactPressure ( ) const

overridevirtualinherited

Returns: The contact force at quadrature points on the mortar segment

Implements WeightedGapUserObject.

Definition at line 113 of file PenaltyWeightedGapUserObject.C.

 {
   return _contact_pressure;
 }

◆ contactTangentialPressureDirOne()

const ADVariableValue & PenaltyFrictionUserObject::contactTangentialPressureDirOne ( ) const

overridevirtual

Returns: The contact force at quadrature points on the mortar segment

Implements WeightedVelocitiesUserObject.

Definition at line 101 of file PenaltyFrictionUserObject.C.

 {
   return _frictional_contact_traction_one;
 }

◆ contactTangentialPressureDirTwo()

const ADVariableValue & PenaltyFrictionUserObject::contactTangentialPressureDirTwo ( ) const

overridevirtual

Implements WeightedVelocitiesUserObject.

Definition at line 107 of file PenaltyFrictionUserObject.C.

 {
   return _frictional_contact_traction_two;
 }

◆ dofToWeightedGap()

const std::unordered_map< const DofObject *, std::pair< ADReal, Real > > & WeightedGapUserObject::dofToWeightedGap ( ) const

inlineinherited

Get the degree of freedom to weighted gap information.

Definition at line 203 of file WeightedGapUserObject.h.

Referenced by ComputeFrictionalForceLMMechanicalContact::incorrectEdgeDroppingPost(), ComputeWeightedGapLMMechanicalContact::incorrectEdgeDroppingPost(), ComputeFrictionalForceLMMechanicalContact::post(), and ComputeWeightedGapLMMechanicalContact::post().

 {
   return _dof_to_weighted_gap;
 }

◆ dofToWeightedVelocities()

const std::unordered_map< const DofObject *, std::array< ADReal, 2 > > & WeightedVelocitiesUserObject::dofToWeightedVelocities ( ) const

inlineinherited

Get the degree of freedom to weighted velocities information.

Definition at line 107 of file WeightedVelocitiesUserObject.h.

Referenced by ComputeFrictionalForceLMMechanicalContact::incorrectEdgeDroppingPost(), and ComputeFrictionalForceLMMechanicalContact::post().

 {
   return _dof_to_weighted_tangential_velocity;
 }

◆ execute()

void WeightedGapUserObject::execute ( )

overridevirtualinherited

Implements MortarUserObject.

Definition at line 160 of file WeightedGapUserObject.C.

 {
   for (_qp = 0; _qp < _qrule_msm->n_points(); _qp++)
   {
     computeQpProperties();
     for (_i = 0; _i < _test->size(); ++_i)
       computeQpIProperties();
   }
 }

◆ finalize()

void PenaltyFrictionUserObject::finalize ( )

overridevirtual

Reimplemented from PenaltyWeightedGapUserObject.

Definition at line 322 of file PenaltyFrictionUserObject.C.

 {
   WeightedVelocitiesUserObject::finalize();
   PenaltyWeightedGapUserObject::selfFinalize();
 }

◆ findValue()

template<typename K , typename V >

V WeightedGapUserObject::findValue	(	const std::unordered_map< K, V > &	map,
		const K &	key,
		const V &	default_value = `0`
	)		const

inlineprotectedinherited

Find a value in a map or return a default if the key doesn't exist.

Definition at line 125 of file WeightedGapUserObject.h.

Referenced by PenaltySimpleCohesiveZoneModel::reinit(), reinit(), and PenaltyWeightedGapUserObject::selfFinalize().

   {
     const auto it = map.find(key);
     if (it == map.end())
       return default_value;
     return it->second;
   }

◆ getAccumulatedSlip()

Real PenaltyFrictionUserObject::getAccumulatedSlip	(	const Node *const	const,
		const unsigned int
	)		const

overridevirtual

Parameters

node	Node pointer
component	Component of the local slip vector

Returns: The accumulated slip at the node

Reimplemented from WeightedGapUserObject.

Definition at line 177 of file PenaltyFrictionUserObject.C.

Referenced by PenaltyMortarUserObjectAux::PenaltyMortarUserObjectAux().

 {
   const auto it = _dof_to_accumulated_slip.find(_subproblem.mesh().nodePtr(node->id()));
 
   if (it != _dof_to_accumulated_slip.end())
     return MetaPhysicL::raw_value(it->second.first(component));
   else
     return 0.0;
 }

◆ getActiveSetState()

virtual bool PenaltyWeightedGapUserObject::getActiveSetState ( const Node *const node ) const

inlinevirtualinherited

Definition at line 45 of file PenaltyWeightedGapUserObject.h.

Referenced by PenaltyMortarUserObjectAux::PenaltyMortarUserObjectAux().

   {
     return _dof_to_normal_pressure.find(_subproblem.mesh().nodePtr(node->id()))->second > 0;
   }

◆ getDeltaTangentialLagrangeMultiplier()

Real PenaltyFrictionUserObject::getDeltaTangentialLagrangeMultiplier	(	const Node *const	node,
		const unsigned int	component
	)		const

overridevirtual

Reimplemented from PenaltyWeightedGapUserObject.

Definition at line 201 of file PenaltyFrictionUserObject.C.

Referenced by PenaltyMortarUserObjectAux::PenaltyMortarUserObjectAux().

 {
   const auto it =
       _dof_to_frictional_lagrange_multipliers.find(_subproblem.mesh().nodePtr(node->id()));
 
   if (it != _dof_to_frictional_lagrange_multipliers.end())
     return MetaPhysicL::raw_value(it->second[component]);
   else
     return 0.0;
 }

◆ getFrictionalContactPressure()

Real PenaltyFrictionUserObject::getFrictionalContactPressure	(	const Node *const	const,
		const unsigned int
	)		const

overridevirtual

Parameters

node	Node pointer
component	Component of the frictional pressure vector

Returns: The frictional contact pressure at the node

Reimplemented from WeightedGapUserObject.

Definition at line 165 of file PenaltyFrictionUserObject.C.

Referenced by PenaltyMortarUserObjectAux::PenaltyMortarUserObjectAux().

 {
   const auto it = _dof_to_tangential_traction.find(_subproblem.mesh().nodePtr(node->id()));
 
   if (it != _dof_to_tangential_traction.end())
     return MetaPhysicL::raw_value(it->second.first(component));
   else
     return 0.0;
 }

◆ getNormalContactPressure()

Real PenaltyWeightedGapUserObject::getNormalContactPressure ( const Node *const const ) const

overridevirtualinherited

Parameters

node	Node pointer

Returns: The normal contact pressure at the node

Reimplemented from WeightedGapUserObject.

Definition at line 162 of file PenaltyWeightedGapUserObject.C.

Referenced by PenaltyMortarUserObjectAux::PenaltyMortarUserObjectAux().

 {
   const auto it = _dof_to_normal_pressure.find(_subproblem.mesh().nodePtr(node->id()));
 
   if (it != _dof_to_normal_pressure.end())
     return MetaPhysicL::raw_value(it->second);
   else
     return 0.0;
 }

◆ getNormalGap()

Real WeightedGapUserObject::getNormalGap ( const Node * const node ) const

virtualinherited

Parameters

node	Node pointer

Returns: The normal gap at the node

Definition at line 171 of file WeightedGapUserObject.C.

Referenced by PenaltyMortarUserObjectAux::PenaltyMortarUserObjectAux(), and PenaltyWeightedGapUserObject::updateAugmentedLagrangianMultipliers().

 {
   const auto it = _dof_to_weighted_gap.find(_subproblem.mesh().nodePtr(node->id()));
 
   // We are returning the physical weighted gap for analysis purposes
   if (it != _dof_to_weighted_gap.end())
     return physicalGap(it->second);
   else
     return 0.0;
 }

◆ getNormalLagrangeMultiplier()

Real PenaltyWeightedGapUserObject::getNormalLagrangeMultiplier ( const Node *const node ) const

virtualinherited

Definition at line 173 of file PenaltyWeightedGapUserObject.C.

Referenced by PenaltyMortarUserObjectAux::PenaltyMortarUserObjectAux().

 {
   const auto it = _dof_to_lagrange_multiplier.find(_subproblem.mesh().nodePtr(node->id()));
 
   if (it != _dof_to_lagrange_multiplier.end())
     return -MetaPhysicL::raw_value(it->second);
   else
     return 0.0;
 }

◆ getTangentialVelocity()

Real PenaltyFrictionUserObject::getTangentialVelocity	(	const Node *const	const,
		const unsigned int
	)		const

overridevirtual

Parameters

node	Node pointer
component	Component of the local slip vector

Returns: The tangential velocity at the node with local components

Reimplemented from WeightedGapUserObject.

Definition at line 189 of file PenaltyFrictionUserObject.C.

 {
   const auto it = _dof_to_real_tangential_velocity.find(_subproblem.mesh().nodePtr(node->id()));
 
   if (it != _dof_to_real_tangential_velocity.end())
     return MetaPhysicL::raw_value(it->second[component]);
   else
     return 0.0;
 }

◆ initialize()

void PenaltyFrictionUserObject::initialize ( )

overridevirtual

Reimplemented from PenaltyWeightedGapUserObject.

Definition at line 154 of file PenaltyFrictionUserObject.C.

 {
   // these functions do not call WeightedGapUserObject::initialize to avoid double initialization
   WeightedVelocitiesUserObject::selfInitialize();
   PenaltyWeightedGapUserObject::selfInitialize();
 
   // instead we call it explicitly here
   WeightedGapUserObject::initialize();
 }

◆ initialSetup()

void WeightedGapUserObject::initialSetup ( )

overridevirtualinherited

Reimplemented from MortarUserObject.

Definition at line 59 of file WeightedGapUserObject.C.

 {
   MortarUserObject::initialSetup();
   _test = &test();
 }

◆ isAugmentedLagrangianConverged()

bool PenaltyFrictionUserObject::isAugmentedLagrangianConverged ( )

overridevirtual

Reimplemented from PenaltyWeightedGapUserObject.

Definition at line 329 of file PenaltyFrictionUserObject.C.

 {
   // save off step slip
   // This method is called at the beginning of the AL iteration.
   for (auto & map_pr : _dof_to_step_slip)
   {
     auto & [step_slip, old_step_slip] = map_pr.second;
     old_step_slip = step_slip;
   }
 
   std::pair<Real, dof_id_type> max_slip{0.0, 0};
 
   for (const auto & [dof_object, traction_pair] : _dof_to_tangential_traction)
   {
     const auto & tangential_traction = traction_pair.first;
     auto normal_pressure = _dof_to_normal_pressure[dof_object];
 
     // We may not find a node in the map because AL machinery gets called at different system
     // configurations. That is, we may want to find a key from a node that computed a zero traction
     // on the verge of not projecting. Since, when we computed the velocity, the system was at a
     // slightly different configuration, that node may not a computed physical weighted velocity.
     // This doesn't seem an issue at all as non-projecting nodes should be a corner case not
     // affecting the physics.
     TwoVector slip_velocity = {0.0, 0.0};
     if (_dof_to_real_tangential_velocity.find(dof_object) != _dof_to_real_tangential_velocity.end())
     {
       const auto & real_tangential_velocity =
           libmesh_map_find(_dof_to_real_tangential_velocity, dof_object);
       slip_velocity = {MetaPhysicL::raw_value(real_tangential_velocity[0]),
                        MetaPhysicL::raw_value(real_tangential_velocity[1])};
     }
 
     // Check slip/stick
     if (_friction_coefficient * normal_pressure < tangential_traction.norm() * (1 + TOLERANCE))
     {
       // If it's slipping, any slip distance is physical.
     }
     else if (slip_velocity.norm() * _dt > _slip_tolerance && normal_pressure > TOLERANCE)
     {
       const auto new_slip =
           std::make_pair<Real, dof_id_type>(slip_velocity.norm() * _dt, dof_object->id());
       if (new_slip > max_slip)
         max_slip = new_slip;
     }
   }
 
   // Communicate max_slip here where all ranks get to
   this->_communicator.max(max_slip);
 
   // Check normal contact convergence now, to make sure all ranks get here
   if (!PenaltyWeightedGapUserObject::isAugmentedLagrangianConverged())
     return false;
 
   // Did we observe any above tolerance slip anywhere?
   if (max_slip.first > _slip_tolerance)
   {
     if (this->_communicator.rank() == 0)
     {
       mooseInfoRepeated("Stick tolerance fails. Slip distance for sticking node ",
                         max_slip.second,
                         " is: ",
                         max_slip.first,
                         ", but slip tolerance is chosen to be ",
                         _slip_tolerance);
     }
     return false;
   }
 
   return true;
 }

◆ physicalGap()

Real WeightedGapUserObject::physicalGap ( const std::pair< ADReal, Real > & gap ) const

inlineinherited

Compute physical gap from integration gap quantity.

Definition at line 57 of file WeightedGapUserObject.h.

Referenced by PenaltyWeightedGapUserObject::augmentedLagrangianSetup(), WeightedGapUserObject::getNormalGap(), and PenaltyWeightedGapUserObject::isAugmentedLagrangianConverged().

   {
     return MetaPhysicL::raw_value(gap.first) / gap.second;
   }

◆ reinit()

void PenaltyFrictionUserObject::reinit ( )

overridevirtual

Reimplemented from PenaltyWeightedGapUserObject.

Definition at line 214 of file PenaltyFrictionUserObject.C.

 {
   // Normal contact pressure with penalty
   PenaltyWeightedGapUserObject::reinit();
 
   // Reset frictional pressure
   _frictional_contact_traction_one.resize(_qrule_msm->n_points());
   _frictional_contact_traction_two.resize(_qrule_msm->n_points()); // 3D
   for (const auto qp : make_range(_qrule_msm->n_points()))
   {
     _frictional_contact_traction_one[qp] = 0.0;
     _frictional_contact_traction_two[qp] = 0.0;
   }
 
   // zero vector
   const static TwoVector zero{0.0, 0.0};
 
   // iterate over nodes
   for (const auto i : make_range(_test->size()))
   {
     // current node
     const Node * const node = _lower_secondary_elem->node_ptr(i);
 
     const auto penalty_friction = findValue(
         _dof_to_local_penalty_friction, static_cast<const DofObject *>(node), _penalty_friction);
 
     // utilized quantities
     const auto & normal_pressure = _dof_to_normal_pressure[node];
 
     // map the tangential traction and accumulated slip
     auto & [tangential_traction, old_tangential_traction] = _dof_to_tangential_traction[node];
     auto & [accumulated_slip, old_accumulated_slip] = _dof_to_accumulated_slip[node];
 
     Real normal_lm = -1;
     if (_dof_to_lagrange_multiplier.find(node) != _dof_to_lagrange_multiplier.end())
       normal_lm = libmesh_map_find(_dof_to_lagrange_multiplier, node);
 
     // Keep active set fixed from second Uzawa loop
     if (normal_lm < -TOLERANCE && normal_pressure > TOLERANCE)
     {
       using namespace std;
 
       const auto & real_tangential_velocity =
           libmesh_map_find(_dof_to_real_tangential_velocity, node);
       const ADTwoVector slip_distance = {real_tangential_velocity[0] * _dt,
                                          real_tangential_velocity[1] * _dt};
 
       // frictional lagrange multiplier (delta lambda^(k)_T)
       const auto & tangential_lm =
           _augmented_lagrange_problem ? _dof_to_frictional_lagrange_multipliers[node] : zero;
 
       // tangential trial traction (Simo 3.12)
       // Modified for implementation in MOOSE: Avoid pingponging on frictional sign (max. 0.4
       // capacity)
       ADTwoVector inner_iteration_penalty_friction = penalty_friction * slip_distance;
 
       const auto slip_metric = MetaPhysicL::raw_value(slip_distance).cwiseAbs().norm();
 
       if (slip_metric > _epsilon_tolerance &&
           penalty_friction * slip_distance.norm() >
               0.4 * _friction_coefficient * std::abs(normal_pressure))
       {
         inner_iteration_penalty_friction =
             MetaPhysicL::raw_value(0.4 * _friction_coefficient * std::abs(normal_pressure) /
                                    (penalty_friction * slip_distance.norm())) *
             penalty_friction * slip_distance;
       }
 
       ADTwoVector tangential_trial_traction =
           old_tangential_traction + tangential_lm + inner_iteration_penalty_friction;
 
       // Nonlinearity below
       ADReal tangential_trial_traction_norm = tangential_trial_traction.norm();
       ADReal phi_trial = tangential_trial_traction_norm - _friction_coefficient * normal_pressure;
       tangential_traction = tangential_trial_traction;
 
       // Simo considers this a 'return mapping'; we are just capping friction to the Coulomb limit.
       if (phi_trial > 0.0)
         // Simo 3.13/3.14 (the penalty formulation has an error in the paper)
         tangential_traction -=
             phi_trial * tangential_trial_traction / tangential_trial_traction_norm;
 
       // track accumulated slip for output purposes
       accumulated_slip = old_accumulated_slip + MetaPhysicL::raw_value(slip_distance).cwiseAbs();
 
       // Keep track of slip vector for adaptive penalty
       auto & [step_slip, old_step_slip] = _dof_to_step_slip[node];
       step_slip = MetaPhysicL::raw_value(slip_distance);
     }
     else
     {
       // reset slip and clear traction
       accumulated_slip.setZero();
       tangential_traction.setZero();
     }
 
     // Now that we have consistent nodal frictional values, create an interpolated frictional
     // pressure variable.
     const auto & test_i = (*_test)[i];
     for (const auto qp : make_range(_qrule_msm->n_points()))
     {
       _frictional_contact_traction_one[qp] += test_i[qp] * tangential_traction(0);
       _frictional_contact_traction_two[qp] += test_i[qp] * tangential_traction(1);
     }
   }
 }

◆ selfFinalize()

void PenaltyWeightedGapUserObject::selfFinalize ( )

protectedinherited

Definition at line 133 of file PenaltyWeightedGapUserObject.C.

Referenced by PenaltySimpleCohesiveZoneModel::finalize(), finalize(), and PenaltyWeightedGapUserObject::finalize().

 {
   // compute new normal pressure for each node
   for (const auto & [dof_object, wgap] : _dof_to_weighted_gap)
   {
     auto penalty = findValue(_dof_to_local_penalty, dof_object, _penalty);
 
     // If _use_physical_gap is true we "normalize" the penalty parameter with the surrounding area.
     auto gap = _use_physical_gap ? adPhysicalGap(wgap) / wgap.second : wgap.first;
 
     const auto lagrange_multiplier =
         _augmented_lagrange_problem ? _dof_to_lagrange_multiplier[dof_object] : 0.0;
 
     // keep the negative normal pressure (compressive per convention here)
     auto normal_pressure = std::min(0.0, penalty * gap + lagrange_multiplier);
 
     // we switch conventins here and consider positive normal pressure as compressive
     _dof_to_normal_pressure[dof_object] = -normal_pressure;
   }
 }

◆ selfInitialize() [1/2]

void WeightedVelocitiesUserObject::selfInitialize ( )

protectedinherited

Definition at line 130 of file WeightedVelocitiesUserObject.C.

Referenced by WeightedVelocitiesUserObject::initialize(), PenaltySimpleCohesiveZoneModel::initialize(), and initialize().

 {
   _dof_to_weighted_tangential_velocity.clear();
   _dof_to_real_tangential_velocity.clear();
 }

◆ selfInitialize() [2/2]

void PenaltyWeightedGapUserObject::selfInitialize ( )

protectedinherited

Definition at line 119 of file PenaltyWeightedGapUserObject.C.

Referenced by PenaltySimpleCohesiveZoneModel::initialize(), initialize(), and PenaltyWeightedGapUserObject::initialize().

 {
   for (auto & dof_pn : _dof_to_normal_pressure)
     dof_pn.second = 0.0;
 }

◆ selfTimestepSetup() [1/2]

void WeightedVelocitiesUserObject::selfTimestepSetup ( )

inlineprotectedinherited

Definition at line 52 of file WeightedVelocitiesUserObject.h.

Referenced by PenaltySimpleCohesiveZoneModel::timestepSetup(), and timestepSetup().

52 {}

◆ selfTimestepSetup() [2/2]

void PenaltyWeightedGapUserObject::selfTimestepSetup ( )

protectedinherited

Definition at line 200 of file PenaltyWeightedGapUserObject.C.

Referenced by PenaltyWeightedGapUserObject::timestepSetup(), PenaltySimpleCohesiveZoneModel::timestepSetup(), and timestepSetup().

 {
   // Let's not clear the LMs for improved performance in
   // nonlinear problems
 
   // reset penalty
   for (auto & dof_lp : _dof_to_local_penalty)
     dof_lp.second = _penalty;
 
   // clear previous gap
   for (auto & dof_pg : _dof_to_previous_gap)
     dof_pg.second = 0.0;
 
   // save old timestep
   _dt_old = _dt;
 }

◆ test()

const VariableTestValue & PenaltyFrictionUserObject::test ( ) const

overrideprotectedvirtual

Returns: The test function associated with the weighted gap

Reimplemented from PenaltyWeightedGapUserObject.

Definition at line 95 of file PenaltyFrictionUserObject.C.

 {
   return _aux_lm_var ? _aux_lm_var->phiLower() : _disp_x_var->phiLower();
 }

◆ timestepSetup()

void PenaltyFrictionUserObject::timestepSetup ( )

overridevirtual

Reimplemented from PenaltyWeightedGapUserObject.

Definition at line 113 of file PenaltyFrictionUserObject.C.

 {
   // these functions do not call WeightedGapUserObject::timestepSetup to avoid double initialization
   WeightedVelocitiesUserObject::selfTimestepSetup();
   PenaltyWeightedGapUserObject::selfTimestepSetup();
 
   // instead we call it explicitly here
   WeightedGapUserObject::timestepSetup();
 
   // Clear step slip (values used in between AL iterations for penalty adaptivity)
   for (auto & map_pr : _dof_to_step_slip)
   {
     auto & [step_slip, old_step_slip] = map_pr.second;
     old_step_slip = {0.0, 0.0};
     step_slip = {0.0, 0.0};
   }
 
   // save off accumulated slip from the last timestep
   for (auto & map_pr : _dof_to_accumulated_slip)
   {
     auto & [accumulated_slip, old_accumulated_slip] = map_pr.second;
     old_accumulated_slip = accumulated_slip;
   }
 
   for (auto & dof_lp : _dof_to_local_penalty_friction)
     dof_lp.second = _penalty_friction;
 
   // save off tangential traction from the last timestep
   for (auto & map_pr : _dof_to_tangential_traction)
   {
     auto & [tangential_traction, old_tangential_traction] = map_pr.second;
     old_tangential_traction = {MetaPhysicL::raw_value(tangential_traction(0)),
                                MetaPhysicL::raw_value(tangential_traction(1))};
     tangential_traction = {0.0, 0.0};
   }
 
   for (auto & [dof_object, delta_tangential_lm] : _dof_to_frictional_lagrange_multipliers)
     delta_tangential_lm.setZero();
 }

◆ updateAugmentedLagrangianMultipliers()

void PenaltyFrictionUserObject::updateAugmentedLagrangianMultipliers ( )

overridevirtual

Reimplemented from PenaltyWeightedGapUserObject.

Definition at line 401 of file PenaltyFrictionUserObject.C.

 {
   PenaltyWeightedGapUserObject::updateAugmentedLagrangianMultipliers();
 
   for (auto & [dof_object, tangential_lm] : _dof_to_frictional_lagrange_multipliers)
   {
     auto & penalty_friction = _dof_to_local_penalty_friction[dof_object];
     if (penalty_friction == 0.0)
       penalty_friction = _penalty_friction;
 
     // normal quantities
     const auto & normal_lm = libmesh_map_find(_dof_to_lagrange_multiplier, dof_object);
 
     // tangential quantities
     const auto & real_tangential_velocity =
         libmesh_map_find(_dof_to_real_tangential_velocity, dof_object);
     const TwoVector slip_velocity = {MetaPhysicL::raw_value(real_tangential_velocity[0]),
                                      MetaPhysicL::raw_value(real_tangential_velocity[1])};
 
     auto & [tangential_traction, old_tangential_traction] = _dof_to_tangential_traction[dof_object];
 
     const TwoVector tangential_trial_traction =
         old_tangential_traction + tangential_lm + penalty_friction * slip_velocity * _dt;
     const Real tangential_trial_traction_norm = tangential_trial_traction.norm();
 
     // Augment
     if (tangential_trial_traction_norm * (1 + TOLERANCE) <=
         std::abs(_friction_coefficient * normal_lm))
     {
       tangential_lm += penalty_friction * slip_velocity * _dt;
     }
     else
     {
       tangential_lm = -tangential_trial_traction / tangential_trial_traction_norm *
                           penalty_friction * normal_lm -
                       old_tangential_traction;
     }
 
     // Update penalty.
     if (_adaptivity_friction == AdaptivityFrictionalPenalty::SIMPLE)
     {
       if (_slip_tolerance < _dt * slip_velocity.norm())
         penalty_friction *= _penalty_multiplier_friction;
 
       // Provide the user the ability of setting this maximum penalty
       if (penalty_friction > _penalty_friction * _max_penalty_multiplier)
         penalty_friction = _penalty_friction * _max_penalty_multiplier;
     }
     else if (_adaptivity_friction == AdaptivityFrictionalPenalty::FRICTION_LIMIT)
     {
       const auto & step_slip = libmesh_map_find(_dof_to_step_slip, dof_object);
       // No change of direction: Adjust penalty factor for the frictional problem
       if (step_slip.first.dot(step_slip.second) > 0.0 && std::abs(normal_lm) > TOLERANCE &&
           _slip_tolerance < _dt * slip_velocity.norm())
       {
         penalty_friction =
             (_friction_coefficient * std::abs(normal_lm)) / 2 / (_dt * slip_velocity.norm());
         // Alternative: accumulated_slip.norm() - old_accumulated_slip.norm()
       }
       // Change of direction: Reduce penalty factor to avoid lack of convergence
       else if (step_slip.first.dot(step_slip.second) < 0.0)
         penalty_friction /= 2.0;
 
       // Heuristics to bound the penalty factor
       if (penalty_friction < _penalty_friction)
         penalty_friction = _penalty_friction;
       else if (penalty_friction > _penalty_friction * _max_penalty_multiplier)
         penalty_friction = _penalty_friction * _max_penalty_multiplier;
     }
   }
 }

◆ validParams()

InputParameters PenaltyFrictionUserObject::validParams ( )

static

Definition at line 24 of file PenaltyFrictionUserObject.C.

 {
   InputParameters params = WeightedVelocitiesUserObject::validParams();
   params += PenaltyWeightedGapUserObject::validParams();
 
   params.addClassDescription(
       "Computes the mortar frictional contact force via a penalty approach.");
   params.addParam<Real>("penalty_friction",
                         "The penalty factor for frictional interaction. If not provide, the normal "
                         "penalty factor is also used for the frictional problem.");
   params.addRequiredParam<Real>("friction_coefficient",
                                 "The friction coefficient ruling Coulomb friction equations.");
   params.addRangeCheckedParam<Real>(
       "slip_tolerance",
       "slip_tolerance > 0",
       "Acceptable slip distance at which augmented Lagrange iterations can be stopped");
   MooseEnum adaptivity_penalty_friction("SIMPLE FRICTION_LIMIT", "FRICTION_LIMIT");
   adaptivity_penalty_friction.addDocumentation(
       "SIMPLE", "Keep multiplying by the frictional penalty multiplier between AL iterations");
   adaptivity_penalty_friction.addDocumentation(
       "FRICTION_LIMIT",
       "This strategy will be guided by the Coulomb limit and be less reliant on the initial "
       "penalty factor provided by the user.");
   params.addParam<MooseEnum>(
       "adaptivity_penalty_friction",
       adaptivity_penalty_friction,
       "The augmented Lagrange update strategy used on the frictional penalty coefficient.");
   params.addRangeCheckedParam<Real>(
       "penalty_multiplier_friction",
       1.0,
       "penalty_multiplier_friction > 0",
       "The penalty growth factor between augmented Lagrange "
       "iterations for penalizing relative slip distance if the node is under stick conditions.");
   return params;
 }

Member Data Documentation

◆ _3d

bool WeightedVelocitiesUserObject::_3d

protectedinherited

Automatic flag to determine whether we are doing three-dimensional work.

Definition at line 103 of file WeightedVelocitiesUserObject.h.

Referenced by WeightedVelocitiesUserObject::computeQpIProperties(), and WeightedVelocitiesUserObject::computeQpProperties().

◆ _adaptivity_friction

enum PenaltyFrictionUserObject::AdaptivityFrictionalPenalty PenaltyFrictionUserObject::_adaptivity_friction

protected

Referenced by updateAugmentedLagrangianMultipliers().

◆ _adaptivity_normal

enum PenaltyWeightedGapUserObject::AdaptivityNormalPenalty PenaltyWeightedGapUserObject::_adaptivity_normal

protectedinherited

Referenced by PenaltyWeightedGapUserObject::updateAugmentedLagrangianMultipliers().

◆ _augmented_lagrange_problem

AugmentedLagrangianContactProblemInterface* const PenaltyWeightedGapUserObject::_augmented_lagrange_problem

protectedinherited

augmented Lagrange problem and iteration number

Definition at line 91 of file PenaltyWeightedGapUserObject.h.

Referenced by PenaltyFrictionUserObject(), PenaltySimpleCohesiveZoneModel::PenaltySimpleCohesiveZoneModel(), PenaltyWeightedGapUserObject::PenaltyWeightedGapUserObject(), PenaltySimpleCohesiveZoneModel::reinit(), reinit(), and PenaltyWeightedGapUserObject::selfFinalize().

◆ _aux_lm_var

const MooseVariable* const PenaltyWeightedGapUserObject::_aux_lm_var

protectedinherited

The auxiliary Lagrange multiplier variable (used together whith the Petrov-Galerkin approach)

Definition at line 115 of file PenaltyWeightedGapUserObject.h.

Referenced by PenaltySimpleCohesiveZoneModel::test(), test(), and PenaltyWeightedGapUserObject::test().

◆ _contact_pressure

ADVariableValue PenaltyWeightedGapUserObject::_contact_pressure

protectedinherited

The contact pressure on the mortar segument quadrature points.

Definition at line 85 of file PenaltyWeightedGapUserObject.h.

Referenced by PenaltyWeightedGapUserObject::contactPressure(), and PenaltyWeightedGapUserObject::reinit().

◆ _coord

const MooseArray<Real>& WeightedGapUserObject::_coord

protectedinherited

Member for handling change of coordinate systems (xyz, rz, spherical)

Definition at line 168 of file WeightedGapUserObject.h.

Referenced by BilinearMixedModeCohesiveZoneModel::computeQpProperties(), WeightedVelocitiesUserObject::computeQpProperties(), and WeightedGapUserObject::computeQpProperties().

◆ _disp_x_var

const MooseVariable* const WeightedGapUserObject::_disp_x_var

protectedinherited

The x displacement variable.

Definition at line 146 of file WeightedGapUserObject.h.

Referenced by WeightedVelocitiesUserObject::computeQpProperties(), WeightedGapUserObject::computeQpProperties(), PenaltyWeightedGapUserObject::PenaltyWeightedGapUserObject(), PenaltySimpleCohesiveZoneModel::test(), test(), and PenaltyWeightedGapUserObject::test().

◆ _disp_y_var

const MooseVariable* const WeightedGapUserObject::_disp_y_var

protectedinherited

The y displacement variable.

Definition at line 148 of file WeightedGapUserObject.h.

Referenced by WeightedVelocitiesUserObject::computeQpProperties(), WeightedGapUserObject::computeQpProperties(), and PenaltyWeightedGapUserObject::PenaltyWeightedGapUserObject().

◆ _disp_z_var

const MooseVariable* const WeightedGapUserObject::_disp_z_var

protectedinherited

The z displacement variable.

Definition at line 152 of file WeightedGapUserObject.h.

Referenced by WeightedVelocitiesUserObject::computeQpProperties(), WeightedGapUserObject::computeQpProperties(), and PenaltyWeightedGapUserObject::PenaltyWeightedGapUserObject().

◆ _dof_to_accumulated_slip

std::unordered_map<const DofObject *, std::pair<TwoVector, TwoVector> > PenaltyFrictionUserObject::_dof_to_accumulated_slip

protected

Map from degree of freedom to current and old accumulated slip.

Definition at line 70 of file PenaltyFrictionUserObject.h.

Referenced by getAccumulatedSlip(), reinit(), and timestepSetup().

◆ _dof_to_frictional_lagrange_multipliers

std::unordered_map<const DofObject *, TwoVector> PenaltyFrictionUserObject::_dof_to_frictional_lagrange_multipliers

protected

Map from degree of freedom to augmented lagrange multiplier.

Definition at line 83 of file PenaltyFrictionUserObject.h.

Referenced by getDeltaTangentialLagrangeMultiplier(), reinit(), timestepSetup(), and updateAugmentedLagrangianMultipliers().

◆ _dof_to_lagrange_multiplier

std::unordered_map<const DofObject *, Real> PenaltyWeightedGapUserObject::_dof_to_lagrange_multiplier

protectedinherited

Map from degree of freedom to augmented lagrange multiplier.

Definition at line 97 of file PenaltyWeightedGapUserObject.h.

Referenced by PenaltyWeightedGapUserObject::getNormalLagrangeMultiplier(), PenaltyWeightedGapUserObject::isAugmentedLagrangianConverged(), PenaltySimpleCohesiveZoneModel::reinit(), reinit(), PenaltyWeightedGapUserObject::selfFinalize(), updateAugmentedLagrangianMultipliers(), and PenaltyWeightedGapUserObject::updateAugmentedLagrangianMultipliers().

◆ _dof_to_local_penalty

std::unordered_map<const DofObject *, Real> PenaltyWeightedGapUserObject::_dof_to_local_penalty

protectedinherited

Map from degree of freedom to local penalty value.

Definition at line 100 of file PenaltyWeightedGapUserObject.h.

Referenced by PenaltyWeightedGapUserObject::selfFinalize(), PenaltyWeightedGapUserObject::selfTimestepSetup(), and PenaltyWeightedGapUserObject::updateAugmentedLagrangianMultipliers().

◆ _dof_to_local_penalty_friction

std::unordered_map<const DofObject *, Real> PenaltyFrictionUserObject::_dof_to_local_penalty_friction

protected

Map from degree of freedom to local friction penalty value.

Definition at line 86 of file PenaltyFrictionUserObject.h.

Referenced by reinit(), timestepSetup(), and updateAugmentedLagrangianMultipliers().

◆ _dof_to_normal_pressure

std::unordered_map<const DofObject *, ADReal> PenaltyWeightedGapUserObject::_dof_to_normal_pressure

protectedinherited

Map from degree of freedom to normal pressure for reporting.

Definition at line 88 of file PenaltyWeightedGapUserObject.h.

Referenced by PenaltyWeightedGapUserObject::getActiveSetState(), PenaltyWeightedGapUserObject::getNormalContactPressure(), isAugmentedLagrangianConverged(), PenaltySimpleCohesiveZoneModel::reinit(), reinit(), PenaltyWeightedGapUserObject::reinit(), PenaltyWeightedGapUserObject::selfFinalize(), and PenaltyWeightedGapUserObject::selfInitialize().

◆ _dof_to_previous_gap

std::unordered_map<const DofObject *, Real> PenaltyWeightedGapUserObject::_dof_to_previous_gap

protectedinherited

Map from degree of freedom to previous AL iteration gap values.

Definition at line 103 of file PenaltyWeightedGapUserObject.h.

Referenced by PenaltyWeightedGapUserObject::augmentedLagrangianSetup(), PenaltyWeightedGapUserObject::selfTimestepSetup(), and PenaltyWeightedGapUserObject::updateAugmentedLagrangianMultipliers().

◆ _dof_to_real_tangential_velocity

std::unordered_map<const DofObject *, std::array<ADReal, 2> > WeightedVelocitiesUserObject::_dof_to_real_tangential_velocity

protectedinherited

A map from node to two interpolated, physical tangential velocities.

Definition at line 58 of file WeightedVelocitiesUserObject.h.

Referenced by WeightedVelocitiesUserObject::computeQpIProperties(), WeightedVelocitiesUserObject::finalize(), getTangentialVelocity(), isAugmentedLagrangianConverged(), PenaltySimpleCohesiveZoneModel::reinit(), reinit(), WeightedVelocitiesUserObject::selfInitialize(), and updateAugmentedLagrangianMultipliers().

◆ _dof_to_step_slip

std::unordered_map<const DofObject *, std::pair<TwoVector, TwoVector> > PenaltyFrictionUserObject::_dof_to_step_slip

protected

Map from degree of freedom to current and old step slip.

Definition at line 67 of file PenaltyFrictionUserObject.h.

Referenced by isAugmentedLagrangianConverged(), reinit(), timestepSetup(), and updateAugmentedLagrangianMultipliers().

◆ _dof_to_tangential_traction

std::unordered_map<const DofObject *, std::pair<ADTwoVector, TwoVector> > PenaltyFrictionUserObject::_dof_to_tangential_traction

protected

Map from degree of freedom to current and old tangential traction.

Definition at line 74 of file PenaltyFrictionUserObject.h.

Referenced by getFrictionalContactPressure(), isAugmentedLagrangianConverged(), reinit(), timestepSetup(), and updateAugmentedLagrangianMultipliers().

◆ _dof_to_weighted_displacements

std::unordered_map<const DofObject *, ADRealVectorValue> WeightedGapUserObject::_dof_to_weighted_displacements

protectedinherited

A map from node to weighted displacements.

Definition at line 180 of file WeightedGapUserObject.h.

Referenced by WeightedGapUserObject::computeQpIProperties(), BilinearMixedModeCohesiveZoneModel::finalize(), WeightedGapUserObject::initialize(), and BilinearMixedModeCohesiveZoneModel::prepareJumpKinematicQuantities().

◆ _dof_to_weighted_gap

std::unordered_map<const DofObject *, std::pair<ADReal, Real> > WeightedGapUserObject::_dof_to_weighted_gap

protectedinherited

A map from node to weighted gap and normalization (if requested)

Definition at line 177 of file WeightedGapUserObject.h.

◆ _dof_to_weighted_tangential_velocity

std::unordered_map<const DofObject *, std::array<ADReal, 2> > WeightedVelocitiesUserObject::_dof_to_weighted_tangential_velocity

protectedinherited

A map from node to two weighted tangential velocities.

Definition at line 55 of file WeightedVelocitiesUserObject.h.

Referenced by WeightedVelocitiesUserObject::computeQpIProperties(), WeightedVelocitiesUserObject::dofToWeightedVelocities(), WeightedVelocitiesUserObject::finalize(), and WeightedVelocitiesUserObject::selfInitialize().

◆ _dt

const Real& PenaltyWeightedGapUserObject::_dt

protectedinherited

Current delta t... or timestep size.

Definition at line 106 of file PenaltyWeightedGapUserObject.h.

Referenced by BilinearMixedModeCohesiveZoneModel::computeDamage(), isAugmentedLagrangianConverged(), PenaltySimpleCohesiveZoneModel::reinit(), reinit(), PenaltyWeightedGapUserObject::selfTimestepSetup(), and updateAugmentedLagrangianMultipliers().

◆ _dt_old

Real PenaltyWeightedGapUserObject::_dt_old

protectedinherited

previous timestep size

Definition at line 109 of file PenaltyWeightedGapUserObject.h.

Referenced by PenaltyWeightedGapUserObject::selfTimestepSetup().

◆ _epsilon_tolerance

const Real PenaltyFrictionUserObject::_epsilon_tolerance

protected

Tolerance to avoid NaN/Inf in automatic differentiation operations.

Definition at line 95 of file PenaltyFrictionUserObject.h.

Referenced by reinit().

◆ _fe_problem

FEProblemBase& WeightedGapUserObject::_fe_problem

protectedinherited

The base finite element problem.

Definition at line 134 of file WeightedGapUserObject.h.

◆ _friction_coefficient

const Real PenaltyFrictionUserObject::_friction_coefficient

protected

The friction coefficient.

Definition at line 64 of file PenaltyFrictionUserObject.h.

Referenced by isAugmentedLagrangianConverged(), reinit(), and updateAugmentedLagrangianMultipliers().

◆ _frictional_contact_traction_one

ADVariableValue PenaltyFrictionUserObject::_frictional_contact_traction_one

protected

The first frictional contact pressure on the mortar segment quadrature points.

Definition at line 77 of file PenaltyFrictionUserObject.h.

Referenced by contactTangentialPressureDirOne(), and reinit().

◆ _frictional_contact_traction_two

ADVariableValue PenaltyFrictionUserObject::_frictional_contact_traction_two

protected

The second frictional contact pressure on the mortar segment quadrature points.

Definition at line 80 of file PenaltyFrictionUserObject.h.

Referenced by contactTangentialPressureDirTwo(), and reinit().

◆ _has_disp_z

const bool WeightedGapUserObject::_has_disp_z

protectedinherited

For 2D mortar contact no displacement will be specified, so const pointers used.

Definition at line 150 of file WeightedGapUserObject.h.

Referenced by WeightedVelocitiesUserObject::computeQpProperties(), WeightedGapUserObject::computeQpProperties(), and PenaltyWeightedGapUserObject::PenaltyWeightedGapUserObject().

◆ _i

unsigned int WeightedGapUserObject::_i = 0

protectedinherited

Test function index.

Definition at line 199 of file WeightedGapUserObject.h.

Referenced by BilinearMixedModeCohesiveZoneModel::computeQpIProperties(), WeightedVelocitiesUserObject::computeQpIProperties(), WeightedGapUserObject::computeQpIProperties(), and WeightedGapUserObject::execute().

◆ _is_weighted_gap_nodal

bool WeightedGapUserObject::_is_weighted_gap_nodal = true

protectedinherited

Whether the weighted gap is associated with nodes or elements (like for a CONSTANT MONOMIAL Lagrange multiplier).

We have this member so that we don't do virtual calls during inner quadrature-point/test-function loops

Definition at line 193 of file WeightedGapUserObject.h.

Referenced by WeightedVelocitiesUserObject::computeQpIProperties().

◆ _lagrangian_iteration_number

const unsigned int& PenaltyWeightedGapUserObject::_lagrangian_iteration_number

protectedinherited

Definition at line 93 of file PenaltyWeightedGapUserObject.h.

Referenced by PenaltyWeightedGapUserObject::isAugmentedLagrangianConverged().

◆ _max_penalty_multiplier

const Real PenaltyWeightedGapUserObject::_max_penalty_multiplier

protectedinherited

Maximum multiplier applied to the initial penalty factor in AL.

Definition at line 118 of file PenaltyWeightedGapUserObject.h.

Referenced by updateAugmentedLagrangianMultipliers(), and PenaltyWeightedGapUserObject::updateAugmentedLagrangianMultipliers().

◆ _no_iterations

const unsigned int PenaltyWeightedGapUserObject::_no_iterations = 0

staticprotectedinherited

Definition at line 92 of file PenaltyWeightedGapUserObject.h.

◆ _nodal

const bool WeightedGapUserObject::_nodal

protectedinherited

Whether the dof objects are nodal; if they're not, then they're elemental.

Definition at line 143 of file WeightedGapUserObject.h.

Referenced by WeightedVelocitiesUserObject::finalize(), WeightedGapUserObject::finalize(), and BilinearMixedModeCohesiveZoneModel::finalize().

◆ _normalization_ptr

const Real* WeightedGapUserObject::_normalization_ptr = nullptr

protectedinherited

Definition at line 184 of file WeightedGapUserObject.h.

◆ _penalty

const Real PenaltyFrictionUserObject::_penalty

protected

The normal penalty factor.

Definition at line 55 of file PenaltyFrictionUserObject.h.

◆ _penalty_friction

const Real PenaltyFrictionUserObject::_penalty_friction

protected

The penalty factor for the frictional constraints.

Definition at line 58 of file PenaltyFrictionUserObject.h.

Referenced by reinit(), timestepSetup(), and updateAugmentedLagrangianMultipliers().

◆ _penalty_multiplier

const Real PenaltyWeightedGapUserObject::_penalty_multiplier

protectedinherited

penalty growth factor for augmented Lagrange

Definition at line 79 of file PenaltyWeightedGapUserObject.h.

Referenced by PenaltyWeightedGapUserObject::updateAugmentedLagrangianMultipliers().

◆ _penalty_multiplier_friction

const Real PenaltyFrictionUserObject::_penalty_multiplier_friction

protected

Penalty growth factor for augmented Lagrange.

Definition at line 89 of file PenaltyFrictionUserObject.h.

Referenced by updateAugmentedLagrangianMultipliers().

◆ _penetration_tolerance

const Real PenaltyWeightedGapUserObject::_penetration_tolerance

protectedinherited

penetration tolerance for augmented Lagrange contact

Definition at line 82 of file PenaltyWeightedGapUserObject.h.

Referenced by PenaltyWeightedGapUserObject::adaptiveNormalPenalty(), PenaltyWeightedGapUserObject::isAugmentedLagrangianConverged(), and PenaltyWeightedGapUserObject::updateAugmentedLagrangianMultipliers().

◆ _primary_disp_x

const ADVariableValue& WeightedGapUserObject::_primary_disp_x

protectedinherited

x-displacement on the primary face

Definition at line 157 of file WeightedGapUserObject.h.

Referenced by WeightedGapUserObject::computeQpProperties().

◆ _primary_disp_y

const ADVariableValue& WeightedGapUserObject::_primary_disp_y

protectedinherited

y-displacement on the primary face

Definition at line 161 of file WeightedGapUserObject.h.

Referenced by WeightedGapUserObject::computeQpProperties().

◆ _primary_disp_z

const ADVariableValue* const WeightedGapUserObject::_primary_disp_z

protectedinherited

z-displacement on the primary face

Definition at line 165 of file WeightedGapUserObject.h.

◆ _primary_var

MooseVariableField<Real>& WeightedVelocitiesUserObject::_primary_var

protectedinherited

Reference to the primary variable.

Definition at line 82 of file WeightedVelocitiesUserObject.h.

◆ _primary_x_dot

const ADVariableValue& WeightedVelocitiesUserObject::_primary_x_dot

protectedinherited

x-velocity on the primary face

Definition at line 88 of file WeightedVelocitiesUserObject.h.

Referenced by WeightedVelocitiesUserObject::computeQpProperties().

◆ _primary_y_dot

const ADVariableValue& WeightedVelocitiesUserObject::_primary_y_dot

protectedinherited

y-velocity on the primary face

Definition at line 94 of file WeightedVelocitiesUserObject.h.

Referenced by WeightedVelocitiesUserObject::computeQpProperties().

◆ _primary_z_dot

const ADVariableValue* const WeightedVelocitiesUserObject::_primary_z_dot

protectedinherited

z-velocity on the primary face

Definition at line 100 of file WeightedVelocitiesUserObject.h.

◆ _qp

unsigned int WeightedGapUserObject::_qp = 0

protectedinherited

Quadrature point index for the mortar segments.

Definition at line 196 of file WeightedGapUserObject.h.

Referenced by BilinearMixedModeCohesiveZoneModel::computeQpIProperties(), WeightedVelocitiesUserObject::computeQpIProperties(), WeightedGapUserObject::computeQpIProperties(), BilinearMixedModeCohesiveZoneModel::computeQpProperties(), WeightedVelocitiesUserObject::computeQpProperties(), WeightedGapUserObject::computeQpProperties(), and WeightedGapUserObject::execute().

◆ _qp_displacement_nodal

ADRealVectorValue WeightedGapUserObject::_qp_displacement_nodal

protectedinherited

Vector for computation of relative displacement (determines mixity ratio in interface problems)

Definition at line 174 of file WeightedGapUserObject.h.

Referenced by WeightedGapUserObject::computeQpIProperties(), and WeightedGapUserObject::computeQpProperties().

◆ _qp_factor

Real WeightedGapUserObject::_qp_factor

protectedinherited

The value of the LM at the current quadrature point.

Definition at line 140 of file WeightedGapUserObject.h.

Referenced by WeightedGapUserObject::computeQpIProperties(), and WeightedGapUserObject::computeQpProperties().

◆ _qp_gap

ADReal WeightedGapUserObject::_qp_gap

protectedinherited

The value of the gap at the current quadrature point.

Definition at line 137 of file WeightedGapUserObject.h.

◆ _qp_gap_nodal

ADRealVectorValue WeightedGapUserObject::_qp_gap_nodal

protectedinherited

Vector for computation of weighted gap with nodal normals.

Definition at line 171 of file WeightedGapUserObject.h.

Referenced by WeightedGapUserObject::computeQpIProperties(), and WeightedGapUserObject::computeQpProperties().

◆ _qp_real_tangential_velocity

std::array<ADReal, 2> WeightedVelocitiesUserObject::_qp_real_tangential_velocity

protectedinherited

The value of the "real" tangential velocity values at the current quadrature point.

Definition at line 67 of file WeightedVelocitiesUserObject.h.

◆ _qp_real_tangential_velocity_nodal

ADRealVectorValue WeightedVelocitiesUserObject::_qp_real_tangential_velocity_nodal

protectedinherited

The value of the real tangential velocity vectors at the current node.

Definition at line 73 of file WeightedVelocitiesUserObject.h.

Referenced by WeightedVelocitiesUserObject::computeQpIProperties(), and WeightedVelocitiesUserObject::computeQpProperties().

◆ _qp_tangential_velocity

std::array<ADReal, 2> WeightedVelocitiesUserObject::_qp_tangential_velocity

protectedinherited

The value of the tangential velocity values at the current quadrature point.

Definition at line 64 of file WeightedVelocitiesUserObject.h.

◆ _qp_tangential_velocity_nodal

ADRealVectorValue WeightedVelocitiesUserObject::_qp_tangential_velocity_nodal

protectedinherited

The value of the tangential velocity vectors at the current node.

Definition at line 70 of file WeightedVelocitiesUserObject.h.

Referenced by WeightedVelocitiesUserObject::computeQpIProperties(), and WeightedVelocitiesUserObject::computeQpProperties().

◆ _secondary_disp_x

const ADVariableValue& WeightedGapUserObject::_secondary_disp_x

protectedinherited

x-displacement on the secondary face

Definition at line 155 of file WeightedGapUserObject.h.

Referenced by WeightedGapUserObject::computeQpProperties().

◆ _secondary_disp_y

const ADVariableValue& WeightedGapUserObject::_secondary_disp_y

protectedinherited

y-displacement on the secondary face

Definition at line 159 of file WeightedGapUserObject.h.

Referenced by WeightedGapUserObject::computeQpProperties().

◆ _secondary_disp_z

const ADVariableValue* const WeightedGapUserObject::_secondary_disp_z

protectedinherited

z-displacement on the secondary face

Definition at line 163 of file WeightedGapUserObject.h.

◆ _secondary_var

MooseVariableField<Real>& WeightedVelocitiesUserObject::_secondary_var

protectedinherited

Reference to the secondary variable.

Definition at line 79 of file WeightedVelocitiesUserObject.h.

◆ _secondary_x_dot

const ADVariableValue& WeightedVelocitiesUserObject::_secondary_x_dot

protectedinherited

x-velocity on the secondary face

Definition at line 85 of file WeightedVelocitiesUserObject.h.

Referenced by WeightedVelocitiesUserObject::computeQpProperties().

◆ _secondary_y_dot

const ADVariableValue& WeightedVelocitiesUserObject::_secondary_y_dot

protectedinherited

y-velocity on the secondary face

Definition at line 91 of file WeightedVelocitiesUserObject.h.

Referenced by WeightedVelocitiesUserObject::computeQpProperties().

◆ _secondary_z_dot

const ADVariableValue* const WeightedVelocitiesUserObject::_secondary_z_dot

protectedinherited

z-velocity on the secondary face

Definition at line 97 of file WeightedVelocitiesUserObject.h.

◆ _slip_tolerance

const Real PenaltyFrictionUserObject::_slip_tolerance

protected

Acceptable slip distance for augmented Lagrange convergence.

Definition at line 61 of file PenaltyFrictionUserObject.h.

Referenced by isAugmentedLagrangianConverged(), and updateAugmentedLagrangianMultipliers().

◆ _sys

SystemBase& WeightedVelocitiesUserObject::_sys

protectedinherited

Reference to the EquationSystem object.

Definition at line 76 of file WeightedVelocitiesUserObject.h.

◆ _tangential_vel_ptr

std::array<const ADReal *, 2> WeightedVelocitiesUserObject::_tangential_vel_ptr = {{nullptr, nullptr}}

protectedinherited

An array of two pointers to avoid copies.

Definition at line 61 of file WeightedVelocitiesUserObject.h.

◆ _test

const VariableTestValue* WeightedGapUserObject::_test = nullptr

protectedinherited

A pointer to the test function associated with the weighted gap.

We have this member so that we don't do virtual calls during inner quadrature-point/test-function loops

Definition at line 188 of file WeightedGapUserObject.h.

Referenced by WeightedGapUserObject::execute(), WeightedGapUserObject::initialSetup(), BilinearMixedModeCohesiveZoneModel::prepareJumpKinematicQuantities(), BilinearMixedModeCohesiveZoneModel::reinit(), PenaltySimpleCohesiveZoneModel::reinit(), reinit(), and PenaltyWeightedGapUserObject::reinit().

◆ _use_physical_gap

const bool PenaltyWeightedGapUserObject::_use_physical_gap

protectedinherited

Use scaled or physical gap.

Definition at line 112 of file PenaltyWeightedGapUserObject.h.

Referenced by PenaltyWeightedGapUserObject::selfFinalize(), and PenaltyWeightedGapUserObject::updateAugmentedLagrangianMultipliers().

◆ _weighted_gap_ptr

const ADReal* WeightedGapUserObject::_weighted_gap_ptr = nullptr

protectedinherited

A pointer members that can be used to help avoid copying ADReals.

Definition at line 183 of file WeightedGapUserObject.h.

The documentation for this class was generated from the following files:

contact/include/userobjects/PenaltyFrictionUserObject.h
contact/src/userobjects/PenaltyFrictionUserObject.C

Public Types

Public Member Functions

Static Public Member Functions

Public Attributes

Static Public Attributes

Protected Types

Protected Member Functions

Static Protected Member Functions

Protected Attributes

Static Protected Attributes

Detailed Description

Member Enumeration Documentation

◆ AdaptivityFrictionalPenalty

◆ AdaptivityNormalPenalty

Constructor & Destructor Documentation

◆ PenaltyFrictionUserObject()

Member Function Documentation

◆ adaptiveNormalPenalty()

◆ adPhysicalGap()

◆ augmentedLagrangianSetup()

◆ bussettaAdaptivePenalty()

◆ computeQpIProperties() [1/2]

◆ computeQpIProperties() [2/2]

◆ computeQpProperties() [1/2]

◆ computeQpProperties() [2/2]

◆ constrainedByOwner()

◆ contactPressure()

◆ contactTangentialPressureDirOne()

◆ contactTangentialPressureDirTwo()

◆ dofToWeightedGap()

◆ dofToWeightedVelocities()

◆ execute()

◆ finalize()

◆ findValue()

◆ getAccumulatedSlip()

◆ getActiveSetState()

◆ getDeltaTangentialLagrangeMultiplier()

◆ getFrictionalContactPressure()

◆ getNormalContactPressure()

◆ getNormalGap()

◆ getNormalLagrangeMultiplier()

◆ getTangentialVelocity()

◆ initialize()

◆ initialSetup()

◆ isAugmentedLagrangianConverged()

◆ physicalGap()

◆ reinit()

◆ selfFinalize()

◆ selfInitialize() [1/2]

◆ selfInitialize() [2/2]

◆ selfTimestepSetup() [1/2]

◆ selfTimestepSetup() [2/2]

◆ test()

◆ timestepSetup()

◆ updateAugmentedLagrangianMultipliers()

◆ validParams()

Member Data Documentation

◆ _3d

◆ _adaptivity_friction

◆ _adaptivity_normal

◆ _augmented_lagrange_problem

◆ _aux_lm_var

◆ _contact_pressure

◆ _coord

◆ _disp_x_var

◆ _disp_y_var

◆ _disp_z_var

◆ _dof_to_accumulated_slip

◆ _dof_to_frictional_lagrange_multipliers

◆ _dof_to_lagrange_multiplier

◆ _dof_to_local_penalty

◆ _dof_to_local_penalty_friction

◆ _dof_to_normal_pressure

◆ _dof_to_previous_gap

◆ _dof_to_real_tangential_velocity

◆ _dof_to_step_slip

◆ _dof_to_tangential_traction

◆ _dof_to_weighted_displacements

◆ _dof_to_weighted_gap

◆ _dof_to_weighted_tangential_velocity