Base class for mortar-based cohesive zone model. More...

#include <CohesiveZoneModelBase.h>

Inheritance diagram for CohesiveZoneModelBase:

Public Types
typedef DataFileName	DataFileParameterType

Public Member Functions
	CohesiveZoneModelBase (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 ADReal	nodalTangentialPressure (const Node &node, unsigned int direction) const override

virtual bool	usesNodalNormalDerivatives () const override
	Whether this formulation includes displacement derivatives of the secondary nodal normals.

virtual void	timestepSetup () override

virtual const ADVariableValue &	czmGlobalTraction (unsigned int i) const

virtual const ADVariableValue &	contactPressure () const override

virtual ADReal	nodalContactPressure (const Node &node) const override

virtual const VariableTestValue &	tractionBasis () const override

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

virtual Real	getNormalLagrangeMultiplier (const Node *const node) const

virtual Real	getDeltaTangentialLagrangeMultiplier (const Node *const, const unsigned int) const

virtual bool	getActiveSetState (const Node *const node) const

virtual bool	isAugmentedLagrangianConverged () override

virtual void	augmentedLagrangianSetup () override

virtual void	updateAugmentedLagrangianMultipliers () 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.

virtual Real	getNormalGap (const Node *const) const

bool	shouldRecordNodalNormalDerivatives () const
	Whether the current assembly records nodal-normal derivatives for this object's system.

const MooseVariable *	dispVar (unsigned int component) const

const ADRealVectorValue &	contactNormal (const Elem &lower_secondary_elem, unsigned int nodal_index) const
	Return the cached contact normal for the supplied lower-dimensional secondary element node.

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

virtual Real	nodalScale (const DofObject *const) const
	Node-based Lagrange-multiplier scaling factor kappa_j of Popp et al.

virtual bool	usesNodalScaling () const

virtual Real	normalizeCDivisor (const DofObject *const, const Real covered_normalization) const

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

virtual Real	getFrictionalContactPressure (const Node *const, const unsigned int) const

virtual Real	getAccumulatedSlip (const Node *const, const unsigned int) const

virtual Real	getTangentialVelocity (const Node *const, const unsigned int) const

virtual void	threadJoin (const UserObject &) override final

virtual Real	spatialValue (const Point &) const

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

void	setPrimaryThreadCopy (UserObject *primary)

UserObject *	primaryThreadCopy ()

SubProblem &	getSubProblem () const

bool	shouldDuplicateInitialExecution () 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)

std::set< UserObjectName >	getDependObjects () const

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

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

unsigned int	systemNumber () const

virtual bool	needThreadedCopy () const

virtual bool	enabled () const

std::shared_ptr< MooseObject >	getSharedPtr ()

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

bool	isKokkosObject () const

MooseApp &	getMooseApp () const

const std::string &	type () const

const std::string &	name () const

std::string	typeAndName () const

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

MooseObjectName	uniqueName () const

const InputParameters &	parameters () const

const hit::Node *	getHitNode () const

bool	hasBase () const

const std::string &	getBase () 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 &name) const

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

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

void	paramWarning (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

std::string	messagePrefix (const bool hit_prefix=true) const

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

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

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

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

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

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

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

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

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

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

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

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

void	callMooseError (std::string msg, const bool with_prefix, const hit::Node *node=nullptr, const bool show_trace=true) 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 UserObjectBase &	getUserObjectBase (const std::string &param_name, bool is_dependency=true) const

const UserObjectBase &	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)

Moose::Kokkos::MaterialProperty< T, dimension >	getKokkosMaterialPropertyByName (const std::string &prop_name_in)

Moose::Kokkos::MaterialProperty< T, dimension >	getKokkosMaterialPropertyOldByName (const std::string &prop_name)

Moose::Kokkos::MaterialProperty< T, dimension >	getKokkosMaterialPropertyOlderByName (const std::string &prop_name)

Moose::Kokkos::MaterialProperty< T, dimension >	getKokkosMaterialProperty (const std::string &name)

Moose::Kokkos::MaterialProperty< T, dimension >	getKokkosMaterialPropertyOld (const std::string &name)

Moose::Kokkos::MaterialProperty< T, dimension >	getKokkosMaterialPropertyOlder (const std::string &name)

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

std::pair< Moose::Kokkos::MaterialProperty< T, dimension >, std::set< SubdomainID > >	getKokkosBlockMaterialProperty (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)

virtual 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

Moose::Kokkos::VariableValue	kokkosCoupledVectorTagValueByName (const std::string &var_name, const std::string &tag_name, unsigned int comp=0) const

Moose::Kokkos::VariableValue	kokkosCoupledVectorTagValuesByName (const std::string &var_name, const std::string &tag_name) const

Moose::Kokkos::VariableGradient	kokkosCoupledVectorTagGradientByName (const std::string &var_name, const std::string &tag_name, unsigned int comp=0) const

Moose::Kokkos::VariableGradient	kokkosCoupledVectorTagGradientsByName (const std::string &var_name, const std::string &tag_name) const

Moose::Kokkos::VariableValue	kokkosCoupledVectorTagNodalValueByName (const std::string &var_name, const std::string &tag_name, unsigned int comp=0) const

Moose::Kokkos::VariableValue	kokkosCoupledVectorTagNodalValuesByName (const std::string &var_name, const std::string &tag_name) const

Moose::Kokkos::VariableValue	kokkosCoupledVectorTagDofValueByName (const std::string &var_name, const std::string &tag_name, unsigned int comp=0) const

Moose::Kokkos::VariableValue	kokkosCoupledVectorTagDofValuesByName (const std::string &var_name, const std::string &tag_name) const

Moose::Kokkos::VariableValue	kokkosCoupledVectorTagValue (const std::string &var_name, const std::string &tag_param_name, unsigned int comp=0) const

Moose::Kokkos::VariableValue	kokkosCoupledVectorTagValues (const std::string &var_name, const std::string &tag_param_name) const

Moose::Kokkos::VariableGradient	kokkosCoupledVectorTagGradient (const std::string &var_name, const std::string &tag_param_name, unsigned int comp=0) const

Moose::Kokkos::VariableGradient	kokkosCoupledVectorTagGradients (const std::string &var_name, const std::string &tag_param_name) const

Moose::Kokkos::VariableValue	kokkosCoupledVectorTagNodalValue (const std::string &var_name, const std::string &tag_param_name, unsigned int comp=0) const

Moose::Kokkos::VariableValue	kokkosCoupledVectorTagNodalValues (const std::string &var_name, const std::string &tag_param_name) const

Moose::Kokkos::VariableValue	kokkosCoupledVectorTagDofValue (const std::string &var_name, const std::string &tag_param_name, unsigned int comp=0) const

Moose::Kokkos::VariableValue	kokkosCoupledVectorTagDofValues (const std::string &var_name, const std::string &tag_param_name) const

Moose::Kokkos::VectorVariableValue	kokkosCoupledVectorTagVectorValueByName (const std::string &var_name, const std::string &tag_name, unsigned int comp=0) const

Moose::Kokkos::VectorVariableValue	kokkosCoupledVectorTagVectorValuesByName (const std::string &var_name, const std::string &tag_name) const

Moose::Kokkos::VectorVariableGradient	kokkosCoupledVectorTagVectorGradientByName (const std::string &var_name, const std::string &tag_name, unsigned int comp=0) const

Moose::Kokkos::VectorVariableGradient	kokkosCoupledVectorTagVectorGradientsByName (const std::string &var_name, const std::string &tag_name) const

Moose::Kokkos::VectorVariableValue	kokkosCoupledVectorTagVectorNodalValueByName (const std::string &var_name, const std::string &tag_name, unsigned int comp=0) const

Moose::Kokkos::VectorVariableValue	kokkosCoupledVectorTagVectorNodalValuesByName (const std::string &var_name, const std::string &tag_name) const

Moose::Kokkos::VectorVariableValue	kokkosCoupledVectorTagVectorDofValueByName (const std::string &var_name, const std::string &tag_name, unsigned int comp=0) const

Moose::Kokkos::VectorVariableValue	kokkosCoupledVectorTagVectorDofValuesByName (const std::string &var_name, const std::string &tag_name) const

Moose::Kokkos::VectorVariableValue	kokkosCoupledVectorTagVectorValue (const std::string &var_name, const std::string &tag_param_name, unsigned int comp=0) const

Moose::Kokkos::VectorVariableValue	kokkosCoupledVectorTagVectorValues (const std::string &var_name, const std::string &tag_param_name) const

Moose::Kokkos::VectorVariableGradient	kokkosCoupledVectorTagVectorGradient (const std::string &var_name, const std::string &tag_param_name, unsigned int comp=0) const

Moose::Kokkos::VectorVariableGradient	kokkosCoupledVectorTagVectorGradients (const std::string &var_name, const std::string &tag_param_name) const

Moose::Kokkos::VectorVariableValue	kokkosCoupledVectorTagVectorNodalValue (const std::string &var_name, const std::string &tag_param_name, unsigned int comp=0) const

Moose::Kokkos::VectorVariableValue	kokkosCoupledVectorTagVectorNodalValues (const std::string &var_name, const std::string &tag_param_name) const

Moose::Kokkos::VectorVariableValue	kokkosCoupledVectorTagVectorDofValue (const std::string &var_name, const std::string &tag_param_name, unsigned int comp=0) const

Moose::Kokkos::VectorVariableValue	kokkosCoupledVectorTagVectorDofValues (const std::string &var_name, const std::string &tag_param_name) const

Moose::Kokkos::VariableValue	kokkosCoupledValue (const std::string &var_name, unsigned int comp=0) const

Moose::Kokkos::VariableValue	kokkosCoupledValues (const std::string &var_name) const

Moose::Kokkos::VariableGradient	kokkosCoupledGradient (const std::string &var_name, unsigned int comp=0) const

Moose::Kokkos::VariableGradient	kokkosCoupledGradients (const std::string &var_name) const

Moose::Kokkos::VariableValue	kokkosCoupledNodalValue (const std::string &var_name, unsigned int comp=0) const

Moose::Kokkos::VariableValue	kokkosCoupledNodalValues (const std::string &var_name) const

Moose::Kokkos::VariableValue	kokkosCoupledDofValue (const std::string &var_name, unsigned int comp=0) const

Moose::Kokkos::VariableValue	kokkosCoupledDofValues (const std::string &var_name) const

Moose::Kokkos::VariableValue	kokkosCoupledValueOld (const std::string &var_name, unsigned int comp=0) const

Moose::Kokkos::VariableValue	kokkosCoupledValuesOld (const std::string &var_name) const

Moose::Kokkos::VariableGradient	kokkosCoupledGradientOld (const std::string &var_name, unsigned int comp=0) const

Moose::Kokkos::VariableGradient	kokkosCoupledGradientsOld (const std::string &var_name) const

Moose::Kokkos::VariableValue	kokkosCoupledNodalValueOld (const std::string &var_name, unsigned int comp=0) const

Moose::Kokkos::VariableValue	kokkosCoupledNodalValuesOld (const std::string &var_name) const

Moose::Kokkos::VariableValue	kokkosCoupledDofValueOld (const std::string &var_name, unsigned int comp=0) const

Moose::Kokkos::VariableValue	kokkosCoupledDofValuesOld (const std::string &var_name) const

Moose::Kokkos::VariableValue	kokkosCoupledValueOlder (const std::string &var_name, unsigned int comp=0) const

Moose::Kokkos::VariableValue	kokkosCoupledValuesOlder (const std::string &var_name) const

Moose::Kokkos::VariableGradient	kokkosCoupledGradientOlder (const std::string &var_name, unsigned int comp=0) const

Moose::Kokkos::VariableGradient	kokkosCoupledGradientsOlder (const std::string &var_name) const

Moose::Kokkos::VariableValue	kokkosCoupledNodalValueOlder (const std::string &var_name, unsigned int comp=0) const

Moose::Kokkos::VariableValue	kokkosCoupledNodalValuesOlder (const std::string &var_name) const

Moose::Kokkos::VariableValue	kokkosCoupledDofValueOlder (const std::string &var_name, unsigned int comp=0) const

Moose::Kokkos::VariableValue	kokkosCoupledDofValuesOlder (const std::string &var_name) const

Moose::Kokkos::VariableValue	kokkosCoupledDot (const std::string &var_name, unsigned int comp=0) const

Moose::Kokkos::VariableValue	kokkosCoupledDots (const std::string &var_name) const

Moose::Kokkos::VariableValue	kokkosCoupledNodalDot (const std::string &var_name, unsigned int comp=0) const

Moose::Kokkos::VariableValue	kokkosCoupledNodalDots (const std::string &var_name) const

Moose::Kokkos::VectorVariableValue	kokkosCoupledVectorValue (const std::string &var_name, unsigned int comp=0) const

Moose::Kokkos::VectorVariableGradient	kokkosCoupledVectorGradient (const std::string &var_name, unsigned int comp=0) const

Moose::Kokkos::VectorVariableValue	kokkosCoupledVectorNodalValue (const std::string &var_name, unsigned int comp=0) const

Moose::Kokkos::VectorVariableValue	kokkosCoupledVectorDofValue (const std::string &var_name, unsigned int comp=0) const

Moose::Kokkos::VectorVariableValue	kokkosCoupledVectorValues (const std::string &var_name) const

Moose::Kokkos::VectorVariableGradient	kokkosCoupledVectorGradients (const std::string &var_name) const

Moose::Kokkos::VectorVariableValue	kokkosCoupledVectorNodalValues (const std::string &var_name) const

Moose::Kokkos::VectorVariableValue	kokkosCoupledVectorDofValues (const std::string &var_name) const

Moose::Kokkos::VectorVariableValue	kokkosCoupledVectorValueOld (const std::string &var_name, unsigned int comp=0) const

Moose::Kokkos::VectorVariableGradient	kokkosCoupledVectorGradientOld (const std::string &var_name, unsigned int comp=0) const

Moose::Kokkos::VectorVariableValue	kokkosCoupledVectorNodalValueOld (const std::string &var_name, unsigned int comp=0) const

Moose::Kokkos::VectorVariableValue	kokkosCoupledVectorDofValueOld (const std::string &var_name, unsigned int comp=0) const

Moose::Kokkos::VectorVariableValue	kokkosCoupledVectorValuesOld (const std::string &var_name) const

Moose::Kokkos::VectorVariableGradient	kokkosCoupledVectorGradientsOld (const std::string &var_name) const

Moose::Kokkos::VectorVariableValue	kokkosCoupledVectorNodalValuesOld (const std::string &var_name) const

Moose::Kokkos::VectorVariableValue	kokkosCoupledVectorDofValuesOld (const std::string &var_name) const

Moose::Kokkos::VectorVariableValue	kokkosCoupledVectorValueOlder (const std::string &var_name, unsigned int comp=0) const

Moose::Kokkos::VectorVariableGradient	kokkosCoupledVectorGradientOlder (const std::string &var_name, unsigned int comp=0) const

Moose::Kokkos::VectorVariableValue	kokkosCoupledVectorNodalValueOlder (const std::string &var_name, unsigned int comp=0) const

Moose::Kokkos::VectorVariableValue	kokkosCoupledVectorDofValueOlder (const std::string &var_name, unsigned int comp=0) const

Moose::Kokkos::VectorVariableValue	kokkosCoupledVectorValuesOlder (const std::string &var_name) const

Moose::Kokkos::VectorVariableGradient	kokkosCoupledVectorGradientsOlder (const std::string &var_name) const

Moose::Kokkos::VectorVariableValue	kokkosCoupledVectorNodalValuesOlder (const std::string &var_name) const

Moose::Kokkos::VectorVariableValue	kokkosCoupledVectorDofValuesOlder (const std::string &var_name) const

Moose::Kokkos::VectorVariableValue	kokkosCoupledVectorDot (const std::string &var_name, unsigned int comp=0) const

Moose::Kokkos::VectorVariableValue	kokkosCoupledVectorNodalDot (const std::string &var_name, unsigned int comp=0) const

Moose::Kokkos::VectorVariableValue	kokkosCoupledVectorDots (const std::string &var_name) const

Moose::Kokkos::VectorVariableValue	kokkosCoupledVectorNodalDots (const std::string &var_name) const

Moose::Kokkos::ADVariableValue	kokkosADCoupledVectorTagValueByName (const std::string &var_name, const std::string &tag_name, unsigned int comp=0) const

Moose::Kokkos::ADVariableValue	kokkosADCoupledVectorTagValuesByName (const std::string &var_name, const std::string &tag_name) const

Moose::Kokkos::ADVariableGradient	kokkosADCoupledVectorTagGradientByName (const std::string &var_name, const std::string &tag_name, unsigned int comp=0) const

Moose::Kokkos::ADVariableGradient	kokkosADCoupledVectorTagGradientsByName (const std::string &var_name, const std::string &tag_name) const

Moose::Kokkos::ADVariableValue	kokkosADCoupledVectorTagNodalValueByName (const std::string &var_name, const std::string &tag_name, unsigned int comp=0) const

Moose::Kokkos::ADVariableValue	kokkosADCoupledVectorTagNodalValuesByName (const std::string &var_name, const std::string &tag_name) const

Moose::Kokkos::ADVariableValue	kokkosADCoupledVectorTagDofValueByName (const std::string &var_name, const std::string &tag_name, unsigned int comp=0) const

Moose::Kokkos::ADVariableValue	kokkosADCoupledVectorTagDofValuesByName (const std::string &var_name, const std::string &tag_name) const

Moose::Kokkos::ADVariableValue	kokkosADCoupledVectorTagValue (const std::string &var_name, const std::string &tag_param_name, unsigned int comp=0) const

Moose::Kokkos::ADVariableValue	kokkosADCoupledVectorTagValues (const std::string &var_name, const std::string &tag_param_name) const

Moose::Kokkos::ADVariableGradient	kokkosADCoupledVectorTagGradient (const std::string &var_name, const std::string &tag_param_name, unsigned int comp=0) const

Moose::Kokkos::ADVariableGradient	kokkosADCoupledVectorTagGradients (const std::string &var_name, const std::string &tag_param_name) const

Moose::Kokkos::ADVariableValue	kokkosADCoupledVectorTagNodalValue (const std::string &var_name, const std::string &tag_param_name, unsigned int comp=0) const

Moose::Kokkos::ADVariableValue	kokkosADCoupledVectorTagNodalValues (const std::string &var_name, const std::string &tag_param_name) const

Moose::Kokkos::ADVariableValue	kokkosADCoupledVectorTagDofValue (const std::string &var_name, const std::string &tag_param_name, unsigned int comp=0) const

Moose::Kokkos::ADVariableValue	kokkosADCoupledVectorTagDofValues (const std::string &var_name, const std::string &tag_param_name) const

Moose::Kokkos::ADVariableValue	kokkosADCoupledValue (const std::string &var_name, unsigned int comp=0) const

Moose::Kokkos::ADVariableValue	kokkosADCoupledValues (const std::string &var_name) const

Moose::Kokkos::ADVariableGradient	kokkosADCoupledGradient (const std::string &var_name, unsigned int comp=0) const

Moose::Kokkos::ADVariableGradient	kokkosADCoupledGradients (const std::string &var_name) const

Moose::Kokkos::ADVariableValue	kokkosADCoupledNodalValue (const std::string &var_name, unsigned int comp=0) const

Moose::Kokkos::ADVariableValue	kokkosADCoupledNodalValues (const std::string &var_name) const

Moose::Kokkos::ADVariableValue	kokkosADCoupledDofValue (const std::string &var_name, unsigned int comp=0) const

Moose::Kokkos::ADVariableValue	kokkosADCoupledDofValues (const std::string &var_name) const

Moose::Kokkos::ADVariableValue	kokkosADCoupledValueOld (const std::string &var_name, unsigned int comp=0) const

Moose::Kokkos::ADVariableValue	kokkosADCoupledValuesOld (const std::string &var_name) const

Moose::Kokkos::ADVariableGradient	kokkosADCoupledGradientOld (const std::string &var_name, unsigned int comp=0) const

Moose::Kokkos::ADVariableGradient	kokkosADCoupledGradientsOld (const std::string &var_name) const

Moose::Kokkos::ADVariableValue	kokkosADCoupledNodalValueOld (const std::string &var_name, unsigned int comp=0) const

Moose::Kokkos::ADVariableValue	kokkosADCoupledNodalValuesOld (const std::string &var_name) const

Moose::Kokkos::ADVariableValue	kokkosADCoupledDofValueOld (const std::string &var_name, unsigned int comp=0) const

Moose::Kokkos::ADVariableValue	kokkosADCoupledDofValuesOld (const std::string &var_name) const

Moose::Kokkos::ADVariableValue	kokkosADCoupledValueOlder (const std::string &var_name, unsigned int comp=0) const

Moose::Kokkos::ADVariableValue	kokkosADCoupledValuesOlder (const std::string &var_name) const

Moose::Kokkos::ADVariableGradient	kokkosADCoupledGradientOlder (const std::string &var_name, unsigned int comp=0) const

Moose::Kokkos::ADVariableGradient	kokkosADCoupledGradientsOlder (const std::string &var_name) const

Moose::Kokkos::ADVariableValue	kokkosADCoupledNodalValueOlder (const std::string &var_name, unsigned int comp=0) const

Moose::Kokkos::ADVariableValue	kokkosADCoupledNodalValuesOlder (const std::string &var_name) const

Moose::Kokkos::ADVariableValue	kokkosADCoupledDofValueOlder (const std::string &var_name, unsigned int comp=0) const

Moose::Kokkos::ADVariableValue	kokkosADCoupledDofValuesOlder (const std::string &var_name) const

Moose::Kokkos::ADVariableValue	kokkosADCoupledDot (const std::string &var_name, unsigned int comp=0) const

Moose::Kokkos::ADVariableValue	kokkosADCoupledDots (const std::string &var_name) const

Moose::Kokkos::ADVariableValue	kokkosADCoupledNodalDot (const std::string &var_name, unsigned int comp=0) const

Moose::Kokkos::ADVariableValue	kokkosADCoupledNodalDots (const std::string &var_name) const

Moose::Kokkos::Scalar< const Real >	kokkosCoupledDotDu (const std::string &var_name, unsigned int comp=0) const

Moose::Kokkos::VariableValue	kokkosZeroValue () const

Moose::Kokkos::VariableGradient	kokkosZeroGradient () const

Moose::Kokkos::VariableValue	kokkosZeroNodalValue () 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	hasKokkosMaterialProperty (const std::string &name)

bool	hasKokkosMaterialProperty (const std::string &name)

bool	hasKokkosMaterialPropertyByName (const std::string &name)

bool	hasKokkosMaterialPropertyByName (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

Moose::Kokkos::Function	getKokkosFunction (const std::string &name) const

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

Moose::Kokkos::Function	getKokkosFunctionByName (const FunctionName &name) const

const T &	getKokkosFunctionByName (const FunctionName &name) const

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

bool	hasKokkosFunctionByName (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.

const std::array< ADRealVectorValue, 2 > &	contactTangents (const Elem &lower_secondary_elem, unsigned int nodal_index) const
	Return contact tangents for the supplied lower-dimensional secondary element node.

virtual const VariableTestValue &	tangentialTractionBasis (unsigned int) const

Static Public Member Functions
static InputParameters	validParams ()

static void	addNodalNormalCoupling (InputParameters &params)
	Request the secondary face point-neighbor matrix coupling that nodal-normal derivatives need.

static void	callMooseError (MooseApp const app, const InputParameters &params, std::string msg, const bool with_prefix, const hit::Node node, const bool show_trace=true)

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, NameFunc &&name_func)

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

static InputParameters	triangulationParams ()

Public Attributes
	usingCombinedWarningSolutionWarnings

const ConsoleStream	_console

Static Public Attributes
static const std::string	type_param

static const std::string	name_param

static const std::string	unique_name_param

static const std::string	app_param

static const std::string	moose_base_param

static const std::string	kokkos_object_param

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 class	AdaptivityNormalPenalty { SIMPLE , BUSSETTA }
	The adaptivity method for the penalty factor at augmentations. More...

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

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

virtual const VariableTestValue &	test () const override

virtual bool	constrainedByOwner () const override

virtual void	prepareJumpKinematicQuantities ()

virtual void	computeFandR (const Node *const node)

virtual void	computeCZMTraction (const Node *const)=0
	Encapsulate the CZM constitutive behavior.

virtual void	computeDamage (const Node *const)=0

virtual void	computeGlobalTraction (const Node *const node)
	Compute global traction for mortar application.

template<class T >
T	normalizeQuantity (const std::unordered_map< const DofObject , T > &map, const Node const node)
	Normalize mortar quantities (remove mortar integral scaling)

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.

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.

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.

ADPoint	nodalCoordinate (const Node &node) const
	Add displacement derivatives to the node's current coordinate.

bool	usesFaceCoupledValues () const override

virtual void	addPostprocessorDependencyHelper (const PostprocessorName &name) const override

virtual void	addVectorPostprocessorDependencyHelper (const VectorPostprocessorName &name) const override

virtual void	addUserObjectDependencyHelper (const UserObjectBase &uo) const override

void	addReporterDependencyHelper (const ReporterName &reporter_name) override

void	flagInvalidSolutionInternal (const InvalidSolutionID invalid_solution_id) const

InvalidSolutionID	registerInvalidSolutionInternal (const std::string &message, const bool warning) const

const ReporterContextBase &	getReporterContextBaseByName (const ReporterName &reporter_name) const

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, const std::string &tag_param_name, unsigned int comp=0) const

const VariableValue &	coupledMatrixTagScalarValue (const std::string &var_name, const std::string &tag_param_name, 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

Moose::Kokkos::VariableValue	kokkosCoupledVectorTagScalarValueByName (const std::string &var_name, const std::string &tag_name, unsigned int comp=0) const

Moose::Kokkos::VariableValue	kokkosCoupledVectorTagScalarValue (const std::string &var_name, const std::string &tag_param_name, unsigned int comp=0) const

Moose::Kokkos::VariableValue	kokkosCoupledScalarValue (const std::string &var_name, unsigned int comp=0) const

Moose::Kokkos::VariableValue	kokkosCoupledScalarValueOld (const std::string &var_name, unsigned int comp=0) const

Moose::Kokkos::VariableValue	kokkosCoupledScalarValueOlder (const std::string &var_name, unsigned int comp=0) const

Moose::Kokkos::VariableValue	kokkosCoupledScalarDot (const std::string &var_name, unsigned int comp=0) const

Moose::Kokkos::Scalar< const Real >	kokkosCoupledScalarDotDu (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)

virtual void	getKokkosMaterialPropertyHook (const std::string &, const unsigned int)

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)

std::vector< const MooseVariableFieldBase * >	getFieldVars (const std::string &var_name) const

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

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

std::vector< const T * >	getVarsHelper (const std::string &var_name) const

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)

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 unsigned int	_ndisp
	Number of displacement components.

std::vector< const GenericVariableGradient< true > * >	_grad_disp
	Coupled displacement gradients.

std::vector< const GenericVariableGradient< true > * >	_grad_disp_neighbor
	Coupled displacement and neighbor displacement gradient.

std::unordered_map< const DofObject *, ADRankTwoTensor >	_dof_to_rotation_matrix
	* Kinematics/displacement jump quantities * Map from degree of freedom to rotation matrix

std::unordered_map< const DofObject *, ADRealVectorValue >	_dof_to_interface_displacement_jump
	Map from degree of freedom to local displacement jump.

ADRankTwoTensor	_F_interpolation
	Deformation gradient for interpolation.

ADRankTwoTensor	_F_neighbor_interpolation
	Deformation gradient for interpolation of the neighbor projection.

std::unordered_map< const DofObject *, ADRankTwoTensor >	_dof_to_F
	Map from degree of freedom to secondary, interpolated deformation gradient tensor.

std::unordered_map< const DofObject *, ADRankTwoTensor >	_dof_to_F_neighbor
	Map from degree of freedom to neighbor, interpolated deformation gradient tensor.

std::unordered_map< const DofObject *, ADRankTwoTensor >	_dof_to_interface_F
	Map from degree of freedom to interface deformation gradient tensor.

std::unordered_map< const DofObject *, ADRankTwoTensor >	_dof_to_interface_R
	Map from degree of freedom to interface rotation tensor.

std::vector< ADVariableValue >	_czm_interpolated_traction
	The global traction.

const Real	_penalty_friction
	The penalty factor for the frictional constraints.

const Real	_friction_coefficient
	The friction coefficient.

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

std::unordered_map< dof_id_type, std::pair< ADTwoVector, TwoVector > > &	_dof_to_accumulated_slip
	Map from degree of freedom to current and old accumulated slip.

std::unordered_map< dof_id_type, std::pair< ADTwoVector, TwoVector > > &	_dof_to_tangential_traction

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

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

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

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

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

std::unordered_map< const DofObject *, ADRealVectorValue >	_dof_to_czm_traction
	Total Lagrangian stress to be applied on CZM interface.

std::unordered_map< const DofObject *, ADRealVectorValue >	_dof_to_displacement_jump
	Displacement jump on CZM interface.

std::unordered_map< dof_id_type, std::pair< ADReal, Real > > &	_dof_to_damage
	Damage values (pair of current and old) on CZM interface.

const Real	_penalty
	The penalty factor.

const Real	_penalty_multiplier
	penalty growth factor for augmented Lagrange

const Real	_penetration_tolerance
	penetration tolerance for augmented Lagrange contact

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

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

std::unordered_map< dof_id_type, Real > &	_dof_to_lagrange_multiplier
	Map from node id to augmented lagrange multiplier.

std::unordered_map< dof_id_type, Real >	_dof_to_lagrange_multiplier_start
	Augmented lagrange multipliers at the start of the current timestep, restored when the timestep is retried so the retry starts from the accepted state.

int &	_t_step_old_normal
	The timestep index at which _dof_to_lagrange_multiplier_start was last saved.

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

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

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

Real	_dt_old
	previous timestep size

const bool	_use_physical_gap
	Use scaled or physical gap.

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

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

enum PenaltyWeightedGapUserObject::AdaptivityNormalPenalty	_adaptivity_normal

FEProblemBase &	_fe_problem
	The base finite element problem.

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

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

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

const MooseVariable *const	_disp_x_var
	The x displacement variable.

const MooseVariable *const	_disp_y_var
	The y displacement variable.

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

const MooseVariable *const	_disp_z_var
	The z displacement variable.

const ADVariableValue &	_secondary_disp_x
	x-displacement on the secondary face

const ADVariableValue &	_primary_disp_x
	x-displacement on the primary face

const ADVariableValue &	_secondary_disp_y
	y-displacement on the secondary face

const ADVariableValue &	_primary_disp_y
	y-displacement on the primary face

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

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

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

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

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

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

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

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

const Real *	_normalization_ptr = nullptr

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

std::unordered_map< const Node *, ADRealVectorValue >	_ad_nodal_normals
	AD nodal normals, each averaged over the secondary faces containing that node.

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

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

unsigned int	_i = 0
	Test function index.

const Moose::CoordinateSystemType &	_coord_sys

const THREAD_ID	_tid

SubProblem &	_subproblem

SystemBase &	_sys

Assembly &	_assembly

const bool	_duplicate_initial_execution

std::set< std::string >	_depend_uo

const bool &	_enabled

MooseApp &	_app

Factory &	_factory

ActionFactory &	_action_factory

const std::string &	_type

const std::string &	_name

const InputParameters &	_pars

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 MortarInterfaceWarehouse &	_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 bool	_is_kokkos_object

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< MooseVariableField< Real > * >	_coupled_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.

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

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

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

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

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

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

SystemBase &	_sys
	Reference to the EquationSystem object.

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

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

const ADVariableValue &	_secondary_x_dot
	x-velocity on the secondary face

const ADVariableValue &	_primary_x_dot
	x-velocity on the primary face

const ADVariableValue &	_secondary_y_dot
	y-velocity on the secondary face

const ADVariableValue &	_primary_y_dot
	y-velocity on the primary face

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

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

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

std::unordered_map< const Node *, std::array< ADRealVectorValue, 2 > >	_ad_nodal_tangents
	AD contact tangent frames, keyed on the secondary node and reused by every test function and quadrature point that interpolates the frictional traction.

Static Protected Attributes
static const std::string	_interpolated_old

static const std::string	_interpolated_older

Private Types
enum	FuncAge

enum	VarType

Private Member Functions
const ADVariableValue *	getADDefaultValue (const std::string &var_name) const

void	checkVar (const std::string &var_name) const

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

const Moose::FunctionBase *	getKokkosFunctionByNameHelper (const FunctionName &name) const

const UserObjectBase &	getUserObjectFromFEProblem (const UserObjectName &object_name, const THREAD_ID tid=0) const

const T &	castUserObject (const UserObjectBase &uo_base, const std::string &param_name="") const

void	mooseObjectError (const std::string &param_name, std::stringstream &oss) const

const std::string &	userObjectType (const UserObjectBase &uo) const

const std::string &	userObjectName (const UserObjectBase &uo) const

const PostprocessorName &	getPostprocessorNameInternal (const std::string &param_name, const unsigned int index, const bool allow_default_value=true) const

bool	isDefaultPostprocessorValueByName (const PostprocessorName &name) const

PostprocessorValue	getDefaultPostprocessorValueByName (const PostprocessorName &name) const

void	checkParam (const std::string &param_name, const unsigned int index=std::numeric_limits< unsigned int >::max()) const

bool	postprocessorsAdded () const

const VectorPostprocessorValue &	getVectorPostprocessorByNameHelper (const VectorPostprocessorName &name, const std::string &vector_name, bool broadcast, std::size_t t_index) const

const VectorPostprocessorContext< VectorPostprocessorValue > &	getVectorPostprocessorContextByNameHelper (const VectorPostprocessorName &name, const std::string &vector_name) const

bool	vectorPostprocessorsAdded () const

bool	reportersAdded () const

void	possiblyCheckHasReporter (const ReporterName &reporter_name, const std::string &param_name="") const

RestartableDataValue &	registerRestartableDataOnApp (std::unique_ptr< RestartableDataValue > data, THREAD_ID tid) const

void	registerRestartableNameWithFilterOnApp (const std::string &name, Moose::RESTARTABLE_FILTER filter)

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

virtual std::string	meshPropertyPrefix (const std::string &data_name) const

const RestartableDataValue &	getMeshPropertyInternal (const std::string &data_name, const std::string &prefix) const

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

const VariableValue *	getDefaultValue (const std::string &var_name) const

Moose::Kokkos::Variable	kokkosCoupledVectorTagScalarVariable (const std::string &var_name, const std::string &tag_name, unsigned int comp) const

void	setNormals ()

Moose::MaterialDataType	getMaterialDataType (const std::set< BoundaryID > &boundary_ids) const

unsigned int	getMaxQps () const

void	addConsumedPropertyName (const MooseObjectName &obj_name, const std::string &prop_name)

const OutputTools< T >::VariableValue &	vectorTagValueHelper (const std::string &var_names, TagID tag, unsigned int index=0) const

const OutputTools< T >::VariableValue &	vectorTagValueHelper (const std::string &var_names, const std::string &tag_name, unsigned int index=0) const

const OutputTools< T >::VariableValue &	vectorTagDofValueHelper (const std::string &var_name, TagID tag, unsigned int comp=0) const

const OutputTools< T >::VariableValue &	vectorTagDofValueHelper (const std::string &var_name, const std::string &tag_name, unsigned int comp=0) const

void	requestStates (const std::string &var_name, const TagName &tag_name, const unsigned int comp)

void	checkFuncType (const std::string var_name, VarType t, FuncAge age) const

const VectorVariableValue *	getDefaultVectorValue (const std::string &var_name) const

const ArrayVariableValue *	getDefaultArrayValue (const std::string &var_name) const

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

const RealVectorValue &	getDefaultNodalValue (const std::string &var_name, unsigned int) const

const RealEigenVector &	getDefaultNodalValue (const std::string &var_name, unsigned int) const

const Moose::Functor< T > &	getDefaultFunctor (const std::string &var_name) const

Moose::Kokkos::Variable	kokkosCoupledVectorTagVariable (const std::string &var_name, const std::string &tag_name, unsigned int comp) const

Moose::Kokkos::Variable	kokkosCoupledVectorTagVariables (const std::string &var_name, const std::string &tag_name) const

Moose::Kokkos::Variable	kokkosCoupledVectorTagVectorVariable (const std::string &var_name, const std::string &tag_name, unsigned int comp) const

Moose::Kokkos::Variable	kokkosCoupledVectorTagVectorVariables (const std::string &var_name, const std::string &tag_name) const

Moose::Kokkos::Variable	kokkosZeroVariable () const

const PostprocessorValue &	getPostprocessorValueInternal (const std::string &param_name, unsigned int index, std::size_t t_index) const

const PostprocessorValue &	getPostprocessorValueInternal (const std::string &param_name, unsigned int index, std::size_t t_index) const

const PostprocessorValue &	getPostprocessorValueByNameInternal (const PostprocessorName &name, std::size_t t_index) const

const PostprocessorValue &	getPostprocessorValueByNameInternal (const PostprocessorName &name, std::size_t t_index) const

void	possiblyCheckHasVectorPostprocessor (const std::string &param_name, const std::string &vector_name) const

void	possiblyCheckHasVectorPostprocessor (const std::string &param_name, const std::string &vector_name) const

void	possiblyCheckHasVectorPostprocessorByName (const VectorPostprocessorName &name, const std::string &vector_name) const

void	possiblyCheckHasVectorPostprocessorByName (const VectorPostprocessorName &name, const std::string &vector_name) const

Static Private Member Functions
static const hit::Node *	getHitNode (const InputParameters &params)

static std::string	messagePrefix (const InputParameters &params, const bool hit_prefix)

Private Attributes
std::unordered_map< std::string, std::unique_ptr< VariableValue > >	_default_value

UserObject *	_primary_thread_copy

std::set< std::string >	_supplied_uo

const ParallelParamObject &	_parent

const MooseBase &	_si_moose_base

const FEProblemBase *	_si_problem

ExecFlagEnum	_empty_execute_enum

const MooseObject &	_fni_object

const InputParameters &	_fni_params

FEProblemBase &	_fni_feproblem

const THREAD_ID	_fni_tid

const MooseObject &	_uoi_moose_object

const FEProblemBase &	_uoi_feproblem

const THREAD_ID	_uoi_tid

const MooseObject &	_ppi_moose_object

const InputParameters &	_ppi_params

const FEProblemBase &	_ppi_feproblem

std::map< PostprocessorName, std::unique_ptr< PostprocessorValue > >	_default_values

const bool	_broadcast_by_default

const MooseObject &	_vpi_moose_object

const FEProblemBase &	_vpi_feproblem

const THREAD_ID	_vpi_tid

const InputParameters &	_ri_params

FEProblemBase &	_ri_fe_problem_base

const ReporterData &	_ri_reporter_data

const MooseObject &	_ri_moose_object

const InputParameters &	_dni_params

FEProblemBase &	_dni_feproblem

const MooseObject *const	_dni_moose_object_ptr

const InputParameters &	_si_params

FEProblemBase &	_si_feproblem

THREAD_ID	_si_tid

const RestartableDataMapName	_metaname

std::string	_restartable_name

MooseApp &	_meta_data_app

const MooseObject *const	_meta_data_object

const InputParameters &	_sc_parameters

const std::string &	_sc_name

const bool	_sc_is_implicit

std::unordered_map< std::string, std::vector< MooseVariableScalar * > >	_coupled_scalar_vars

std::unordered_map< std::string, std::unique_ptr< ADVariableValue > >	_dual_default_value

std::vector< MooseVariableScalar * >	_coupled_moose_scalar_vars

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

std::set< TagID >	_sc_coupleable_vector_tags

std::set< TagID >	_sc_coupleable_matrix_tags

const AutomaticMortarGeneration *	_amg

const MooseObject &	_mi_moose_object

const bool	_mi_boundary_restricted

const std::set< SubdomainID > &	_mi_block_ids

const std::set< BoundaryID > &	_mi_boundary_ids

std::vector< std::unique_ptr< OptionalMaterialPropertyProxyBase< MaterialPropertyInterface > > >	_optional_property_proxies

unsigned int	_coupleable_max_qps

std::unordered_map< std::string, std::vector< unsigned int > >	_optional_var_index

std::unordered_map< std::string, std::vector< MooseVariableScalar * > >	_c_coupled_scalar_vars

std::set< TagID >	_fe_coupleable_vector_tags

std::set< TagID >	_fe_coupleable_matrix_tags

const bool	_is_fv

const MooseObject *const	_obj

const std::set< std::string >	_older_state_tags

std::vector< std::set< MooseWritableVariable * > >	_writable_coupled_variables

const Real	_kokkos_default_value_zero

AugmentedLagrangianContactProblemInterface *const	_augmented_lagrange_problem
	augmented Lagrange problem and iteration number

const unsigned int &	_lagrangian_iteration_number

static const unsigned int	_no_iterations = 0

Detailed Description

Base class for mortar-based cohesive zone model.

Definition at line 20 of file CohesiveZoneModelBase.h.

Member Enumeration Documentation

◆ AdaptivityNormalPenalty

enum class PenaltyWeightedGapUserObject::AdaptivityNormalPenalty

strongprotectedinherited

The adaptivity method for the penalty factor at augmentations.

Enumerator
SIMPLE
BUSSETTA

Definition at line 134 of file PenaltyWeightedGapUserObject.h.

134{ SIMPLE, BUSSETTA } _adaptivity_normal;

PenaltyWeightedGapUserObject::_adaptivity_normal

enum PenaltyWeightedGapUserObject::AdaptivityNormalPenalty _adaptivity_normal

PenaltyWeightedGapUserObject::AdaptivityNormalPenalty::BUSSETTA

@ BUSSETTA

PenaltyWeightedGapUserObject::AdaptivityNormalPenalty::SIMPLE

@ SIMPLE

Constructor & Destructor Documentation

◆ CohesiveZoneModelBase()

CohesiveZoneModelBase::CohesiveZoneModelBase ( const InputParameters & parameters )

Definition at line 51 of file CohesiveZoneModelBase.C.

  : WeightedGapUserObject(parameters),
    PenaltyWeightedGapUserObject(parameters),
    WeightedVelocitiesUserObject(parameters),
    _ndisp(coupledComponents("displacements")),
    _penalty_friction(isParamValid("penalty_friction") ? getParam<Real>("penalty_friction")
                                                       : getParam<Real>("penalty")),
    _friction_coefficient(getParam<Real>("friction_coefficient")),
    _dof_to_accumulated_slip(
        declareRestartableData<std::unordered_map<dof_id_type, std::pair<ADTwoVector, TwoVector>>>(
            "dof_to_accumulated_slip",
            std::unordered_map<dof_id_type, std::pair<ADTwoVector, TwoVector>>{})),
    _dof_to_tangential_traction(
        declareRestartableData<std::unordered_map<dof_id_type, std::pair<ADTwoVector, TwoVector>>>(
            "dof_to_tangential_traction",
            std::unordered_map<dof_id_type, std::pair<ADTwoVector, TwoVector>>{})),
    _epsilon_tolerance(1.0e-40),
    _dof_to_damage(declareRestartableData<std::unordered_map<dof_id_type, std::pair<ADReal, Real>>>(
        "dof_do_damage", std::unordered_map<dof_id_type, std::pair<ADReal, Real>>{}))
{
  _czm_interpolated_traction.resize(_ndisp);
 
  for (unsigned int i = 0; i < _ndisp; ++i)
  {
    _grad_disp.push_back(&adCoupledGradient("displacements", i));
    _grad_disp_neighbor.push_back(&adCoupledNeighborGradient("displacements", i));
  }
 
  // Set non-intervening components to zero
  for ([[maybe_unused]] const auto i : make_range(_ndisp, 3u))
  {
    _grad_disp.push_back(&_ad_grad_zero);
    _grad_disp_neighbor.push_back(&_ad_grad_zero);
  }
 
  if (_augmented_lagrange_problem)
    mooseError("CohesiveZoneModelBase constraints cannot be enforced with an augmented "
               "Lagrange approach.");
}

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 397 of file PenaltyWeightedGapUserObject.C.

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

Referenced by PenaltyWeightedGapUserObject::bussettaAdaptivePenalty().

◆ addNodalNormalCoupling()

void WeightedGapUserObject::addNodalNormalCoupling ( InputParameters & params )

staticinherited

Request the secondary face point-neighbor matrix coupling that nodal-normal derivatives need.

MortarConsumerInterface supplies the geometric and algebraic relationship manager; this adds the coupling relationship manager, which is attached through a separate system lifecycle.

Definition at line 40 of file WeightedGapUserObject.C.

{
  params.set<bool>("ghost_point_neighbors") = true;
  const auto configure_point_neighbors =
      [](const InputParameters & obj_params, InputParameters & rm_params)
  {
    rm_params.set<bool>("use_displaced_mesh") = obj_params.get<bool>("use_displaced_mesh");
    rm_params.set<BoundaryName>("secondary_boundary") =
        obj_params.get<BoundaryName>("secondary_boundary");
    rm_params.set<BoundaryName>("primary_boundary") =
        obj_params.get<BoundaryName>("primary_boundary");
    rm_params.set<SubdomainName>("secondary_subdomain") =
        obj_params.get<SubdomainName>("secondary_subdomain");
    rm_params.set<SubdomainName>("primary_subdomain") =
        obj_params.get<SubdomainName>("primary_subdomain");
    // Objects that disable ghost_point_neighbors do not include nodal-normal derivatives.
    rm_params.set<bool>("ghost_point_neighbors") = obj_params.get<bool>("ghost_point_neighbors");
  };
  params.addRelationshipManager("AugmentSparsityOnInterface",
                                Moose::RelationshipManagerType::COUPLING,
                                configure_point_neighbors);
}

Referenced by LMWeightedGapUserObject::validParams(), LMWeightedVelocitiesUserObject::validParams(), and PenaltyWeightedGapUserObject::validParams().

◆ adPhysicalGap()

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

inlineinherited

Definition at line 128 of file WeightedGapUserObject.h.

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

Referenced by PenaltyWeightedGapUserObject::selfFinalize().

◆ augmentedLagrangianSetup()

void PenaltyWeightedGapUserObject::augmentedLagrangianSetup ( )

overridevirtualinherited

Implements AugmentedLagrangeInterface.

Definition at line 318 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 374 of file PenaltyWeightedGapUserObject.C.

{
  // 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);
  }
}

Referenced by PenaltyWeightedGapUserObject::updateAugmentedLagrangianMultipliers().

◆ computeCZMTraction()

virtual void CohesiveZoneModelBase::computeCZMTraction ( const Node * const )

protectedpure virtual

Encapsulate the CZM constitutive behavior.

Implemented in BilinearMixedModeCohesiveZoneModel.

Referenced by reinit().

◆ computeDamage()

virtual void CohesiveZoneModelBase::computeDamage ( const Node * const )

protectedpure virtual

Implemented in BilinearMixedModeCohesiveZoneModel.

◆ computeFandR()

void CohesiveZoneModelBase::computeFandR ( const Node *const node )

protectedvirtual

Definition at line 144 of file CohesiveZoneModelBase.C.

{
  // First call does not have maps available
  const bool return_boolean = _dof_to_F.find(node) == _dof_to_F.end();
  if (return_boolean)
    return;
 
  const auto normalized_F = normalizeQuantity(_dof_to_F, node);
  const auto normalized_F_neighbor = normalizeQuantity(_dof_to_F_neighbor, node);
 
  // This 'averaging' assumption below can probably be improved upon.
  _dof_to_interface_F[node] = 0.5 * (normalized_F + normalized_F_neighbor);
 
  for (const auto i : make_range(Moose::dim))
    for (const auto j : make_range(Moose::dim))
    {
      if (!std::isfinite(MetaPhysicL::raw_value(normalized_F(i, j))))
        throw MooseException("The deformation gradient on the secondary surface is not finite in "
                             "CohesiveZoneModelBase. MOOSE needs to cut the time step size.");
      if (!std::isfinite(MetaPhysicL::raw_value(normalized_F_neighbor(i, j))))
        throw MooseException("The deformation gradient on the primary surface is not finite in "
                             "CohesiveZoneModelBase. MOOSE needs to cut the time step size.");
    }
 
  const auto dof_to_interface_F_node = libmesh_map_find(_dof_to_interface_F, node);
 
  const ADFactorizedRankTwoTensor C = dof_to_interface_F_node.transpose() * dof_to_interface_F_node;
  const auto Uinv = MathUtils::sqrt(C).inverse().get();
  _dof_to_interface_R[node] = dof_to_interface_F_node * Uinv;
 
  // Transform interface jump according to two rotation matrices
  const auto global_interface_displacement = _dof_to_interface_displacement_jump[node];
  _dof_to_interface_displacement_jump[node] =
      (_dof_to_interface_R[node] * _dof_to_rotation_matrix[node]).transpose() *
      global_interface_displacement;
}

Referenced by reinit().

◆ computeGlobalTraction()

void CohesiveZoneModelBase::computeGlobalTraction ( const Node *const node )

protectedvirtual

Compute global traction for mortar application.

Definition at line 461 of file CohesiveZoneModelBase.C.

{
  // First call does not have maps available
  const auto it = _dof_to_czm_traction.find(node);
  if (it == _dof_to_czm_traction.end())
    return;
 
  const auto local_traction_vector = it->second;
  const auto rotation_matrix = libmesh_map_find(_dof_to_rotation_matrix, node);
 
  _dof_to_czm_traction[node] = rotation_matrix * local_traction_vector;
}

Referenced by reinit().

◆ computeQpIProperties()

void CohesiveZoneModelBase::computeQpIProperties ( )

overrideprotectedvirtual

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

Reimplemented from WeightedGapUserObject.

Reimplemented in BilinearMixedModeCohesiveZoneModel.

Definition at line 130 of file CohesiveZoneModelBase.C.

{
  WeightedVelocitiesUserObject::computeQpIProperties();
  // Get the _dof_to_weighted_gap map
  const auto * const dof = cast_ptr<const DofObject *>(_lower_secondary_elem->node_ptr(_i));
 
  _dof_to_weighted_displacements[dof] += (*_test)[_i][_qp] * _qp_displacement_nodal;
 
  // TODO: Probably better to interpolate the deformation gradients.
  _dof_to_F[dof] += (*_test)[_i][_qp] * _F_interpolation;
  _dof_to_F_neighbor[dof] += (*_test)[_i][_qp] * _F_neighbor_interpolation;
}

Referenced by BilinearMixedModeCohesiveZoneModel::computeQpIProperties().

◆ computeQpProperties()

void CohesiveZoneModelBase::computeQpProperties ( )

overrideprotectedvirtual

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 111 of file CohesiveZoneModelBase.C.

{
  WeightedVelocitiesUserObject::computeQpProperties();
 
  // Compute F and R.
  const auto F = (ADRankTwoTensor::Identity() +
                  ADRankTwoTensor::initializeFromRows(
                      (*_grad_disp[0])[_qp], (*_grad_disp[1])[_qp], (*_grad_disp[2])[_qp]));
  const auto F_neighbor = (ADRankTwoTensor::Identity() +
                           ADRankTwoTensor::initializeFromRows((*_grad_disp_neighbor[0])[_qp],
                                                               (*_grad_disp_neighbor[1])[_qp],
                                                               (*_grad_disp_neighbor[2])[_qp]));
 
  // TODO in follow-on PRs: Trim interior node variable derivatives
  _F_interpolation = F * (_JxW_msm[_qp] * _coord[_qp]);
  _F_neighbor_interpolation = F_neighbor * (_JxW_msm[_qp] * _coord[_qp]);
}

Referenced by BilinearMixedModeCohesiveZoneModel::computeQpProperties().

◆ constrainedByOwner()

virtual bool CohesiveZoneModelBase::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.

Reimplemented in BilinearMixedModeCohesiveZoneModel.

Definition at line 44 of file CohesiveZoneModelBase.h.

44{ return false; }

Referenced by finalize().

◆ contactNormal()

const ADRealVectorValue & WeightedGapUserObject::contactNormal	(	const Elem &	lower_secondary_elem,
		unsigned int	nodal_index
	)		const

inherited

Return the cached contact normal for the supplied lower-dimensional secondary element node.

The raw value is the stored normalized secondary nodal normal. Coordinate derivatives are included while this object's formulation and the current assembly mode require them.

Definition at line 222 of file WeightedGapUserObject.C.

{
  mooseAssert(nodal_index < lower_secondary_elem.n_nodes(),
              "Nodal normal index must refer to a node on the secondary element.");
  mooseAssert(shouldRecordNodalNormalDerivatives(),
              "AD contact normals should only be requested while recording nodal-normal "
              "derivatives.");
  const Node * const node = lower_secondary_elem.node_ptr(nodal_index);
  const auto normal_it = _ad_nodal_normals.find(node);
  if (normal_it != _ad_nodal_normals.end())
    return normal_it->second;
 
  amg().computeADNodalNormals([this](const Node & coordinate_node)
                              { return nodalCoordinate(coordinate_node); },
                              _ad_nodal_normals);
  return libmesh_map_find(_ad_nodal_normals, node);
}

Referenced by WeightedGapUserObject::computeQpIProperties(), WeightedVelocitiesUserObject::contactTangents(), and NormalMortarMechanicalContact::precomputeQpQuantities().

◆ contactPressure()

const ADVariableValue & PenaltyWeightedGapUserObject::contactPressure ( ) const

overridevirtualinherited

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

Implements WeightedGapUserObject.

Definition at line 126 of file PenaltyWeightedGapUserObject.C.

{
  return _contact_pressure;
}

◆ contactTangentialPressureDirOne()

const ADVariableValue & CohesiveZoneModelBase::contactTangentialPressureDirOne ( ) const

overridevirtual

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

Implements WeightedVelocitiesUserObject.

Definition at line 449 of file CohesiveZoneModelBase.C.

{
  return _frictional_contact_traction_one;
}

◆ contactTangentialPressureDirTwo()

const ADVariableValue & CohesiveZoneModelBase::contactTangentialPressureDirTwo ( ) const

overridevirtual

Implements WeightedVelocitiesUserObject.

Definition at line 455 of file CohesiveZoneModelBase.C.

{
  return _frictional_contact_traction_two;
}

◆ contactTangents()

const std::array< ADRealVectorValue, 2 > & WeightedVelocitiesUserObject::contactTangents	(	const Elem &	lower_secondary_elem,
		unsigned int	nodal_index
	)		const

inherited

Return contact tangents for the supplied lower-dimensional secondary element node.

Raw values are the stored mechanical-contact nodal tangents; derivatives are included when the internal quasistatic normal-derivative path and AD derivative recording are enabled.

Definition at line 67 of file WeightedVelocitiesUserObject.C.

{
  mooseAssert(nodal_index < lower_secondary_elem.n_nodes(),
              "Nodal tangent index must refer to a node on the secondary element.");
  mooseAssert(shouldRecordNodalNormalDerivatives(),
              "AD contact tangents should only be requested while recording nodal-normal "
              "derivatives.");
  const Node * const node = lower_secondary_elem.node_ptr(nodal_index);
  auto it = _ad_nodal_tangents.find(node);
  if (it == _ad_nodal_tangents.end())
    it = _ad_nodal_tangents
             .emplace(node,
                      Moose::Mortar::householderTangents(
                          contactNormal(lower_secondary_elem, nodal_index)))
             .first;
 
#ifndef NDEBUG
  // A cached frame must belong to the current geometry. The tolerance matches the stored versus AD
  // nodal normal check in AutomaticMortarGeneration::computeADNodalNormals.
  const auto & stored_tangents = amg().getNodalTangents(lower_secondary_elem);
  for (const auto direction : make_range(2))
    mooseAssert(
        (MetaPhysicL::raw_value(it->second[direction]) - stored_tangents[direction][nodal_index])
                .norm() < 1e-10,
        "The stored and AD secondary nodal tangents must use the same geometry state.");
#endif
 
  return it->second;
}

Referenced by WeightedVelocitiesUserObject::computeQpIProperties(), and TangentialMortarMechanicalContact::precomputeQpQuantities().

◆ czmGlobalTraction()

const ADVariableValue & CohesiveZoneModelBase::czmGlobalTraction ( unsigned int i ) const

virtual

Definition at line 475 of file CohesiveZoneModelBase.C.

{
  mooseAssert(i <= 3,
              "Internal error in czmGlobalTraction. Index exceeds the traction vector size.");
 
  return _czm_interpolated_traction[i];
}

Referenced by MortarGenericTraction::computeQpResidual().

◆ dispVar()

const MooseVariable * WeightedGapUserObject::dispVar ( unsigned int component ) const

inherited

Parameters

component The displacement component index (0=x, 1=y, 2=z)

Returns: This object's displacement variable for that component, or nullptr if the problem is not three-dimensional and component is 2

Definition at line 103 of file WeightedGapUserObject.C.

{
  switch (component)
  {
    case 0:
      return _disp_x_var;
    case 1:
      return _disp_y_var;
    case 2:
      return _disp_z_var;
    default:
      mooseError("Invalid displacement component '", component, "' requested.");
  }
}

Referenced by ComputeFrictionalForceLMMechanicalContact::ComputeFrictionalForceLMMechanicalContact(), ComputeWeightedGapLMMechanicalContact::ComputeWeightedGapLMMechanicalContact(), NormalMortarMechanicalContact::initialSetup(), and TangentialMortarMechanicalContact::initialSetup().

◆ 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 277 of file WeightedGapUserObject.h.

{
  return _dof_to_weighted_gap;
}

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

◆ 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 146 of file WeightedVelocitiesUserObject.h.

{
  return _dof_to_weighted_tangential_velocity;
}

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

◆ execute()

void WeightedGapUserObject::execute ( )

overridevirtualinherited

Implements MortarUserObject.

Definition at line 293 of file WeightedGapUserObject.C.

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

◆ finalize()

void CohesiveZoneModelBase::finalize ( )

overridevirtual

Reimplemented from PenaltyWeightedGapUserObject.

Reimplemented in BilinearMixedModeCohesiveZoneModel.

Definition at line 431 of file CohesiveZoneModelBase.C.

{
  WeightedVelocitiesUserObject::finalize();
  PenaltyWeightedGapUserObject::selfFinalize();
 
  const bool send_data_back = !constrainedByOwner();
 
  Moose::Mortar::Contact::communicateR2T(
      _dof_to_F, _subproblem.mesh(), _nodal, _communicator, send_data_back);
 
  Moose::Mortar::Contact::communicateR2T(
      _dof_to_F_neighbor, _subproblem.mesh(), _nodal, _communicator, send_data_back);
 
  Moose::Mortar::Contact::communicateRealObject(
      _dof_to_weighted_displacements, _subproblem.mesh(), _nodal, _communicator, send_data_back);
}

Referenced by BilinearMixedModeCohesiveZoneModel::finalize().

◆ 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 191 of file WeightedGapUserObject.h.

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

Referenced by PenaltyWeightedGapUserObject::nodalContactPressure(), LMWeightedGapUserObject::nodalScale(), LMWeightedGapUserObject::normalizeCDivisor(), reinit(), PenaltyFrictionUserObject::reinit(), and PenaltyWeightedGapUserObject::selfFinalize().

◆ getAccumulatedSlip()

virtual Real WeightedGapUserObject::getAccumulatedSlip	(	const Node * const	,
		const unsigned int
	)		const

inlinevirtualinherited

Parameters

node	Node pointer
component	Component of the local slip vector

Returns: The accumulated slip at the node

Reimplemented in PenaltyFrictionUserObject.

Definition at line 146 of file WeightedGapUserObject.h.

  {
    mooseError("Not available in base class.");
  }

◆ getActiveSetState()

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

inlinevirtualinherited

Definition at line 48 of file PenaltyWeightedGapUserObject.h.

  {
    const auto it = _dof_to_normal_pressure.find(_subproblem.mesh().nodePtr(node->id()));
    return it != _dof_to_normal_pressure.end() && it->second > 0;
  }

Referenced by MortarUserObjectAux::MortarUserObjectAux().

◆ getDeltaTangentialLagrangeMultiplier()

virtual Real PenaltyWeightedGapUserObject::getDeltaTangentialLagrangeMultiplier	(	const Node * const	,
		const unsigned int
	)		const

inlinevirtualinherited

Reimplemented in PenaltyFrictionUserObject.

Definition at line 44 of file PenaltyWeightedGapUserObject.h.

  {
    return 0.0;
  };

◆ getFrictionalContactPressure()

virtual Real WeightedGapUserObject::getFrictionalContactPressure	(	const Node * const	,
		const unsigned int
	)		const

inlinevirtualinherited

Parameters

node	Node pointer
component	Component of the frictional pressure vector

Returns: The frictional contact pressure at the node

Reimplemented in PenaltyFrictionUserObject.

Definition at line 135 of file WeightedGapUserObject.h.

  {
    mooseError("Not available in base class.");
  }

◆ getNormalContactPressure()

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

overridevirtualinherited

Parameters

node	Node pointer

Returns: The normal contact pressure at the node

Implements WeightedGapUserObject.

Definition at line 177 of file PenaltyWeightedGapUserObject.C.

{
  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 304 of file WeightedGapUserObject.C.

{
  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;
}

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

◆ getNormalLagrangeMultiplier()

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

virtualinherited

Definition at line 188 of file PenaltyWeightedGapUserObject.C.

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

Referenced by MortarUserObjectAux::MortarUserObjectAux().

◆ getTangentialVelocity()

virtual Real WeightedGapUserObject::getTangentialVelocity	(	const Node * const	,
		const unsigned int
	)		const

inlinevirtualinherited

Parameters

node	Node pointer
component	Component of the local slip vector

Returns: The tangential velocity at the node with local components

Reimplemented in PenaltyFrictionUserObject.

Definition at line 157 of file WeightedGapUserObject.h.

  {
    mooseError("Not available in base class.");
  }

Referenced by MortarUserObjectAux::MortarUserObjectAux().

◆ initialize()

void CohesiveZoneModelBase::initialize ( )

overridevirtual

Reimplemented from PenaltyWeightedGapUserObject.

Reimplemented in BilinearMixedModeCohesiveZoneModel.

Definition at line 231 of file CohesiveZoneModelBase.C.

{
  // these functions do not call WeightedGapUserObject::initialize to avoid double initialization
  WeightedVelocitiesUserObject::selfInitialize();
  PenaltyWeightedGapUserObject::selfInitialize();
 
  // instead we call it explicitly here
  WeightedGapUserObject::initialize();
  _dof_to_F.clear();
  _dof_to_F_neighbor.clear();
  _dof_to_interface_displacement_jump.clear();
  _dof_to_interface_F.clear();
  _dof_to_interface_R.clear();
  _dof_to_czm_traction.clear();
  for (auto & map_pr : _dof_to_rotation_matrix)
    map_pr.second.setToIdentity();
}

Referenced by BilinearMixedModeCohesiveZoneModel::initialize().

◆ initialSetup()

void WeightedGapUserObject::initialSetup ( )

overridevirtualinherited

Reimplemented from MortarUserObject.

Definition at line 119 of file WeightedGapUserObject.C.

{
  MortarUserObject::initialSetup();
  _test = &test();
 
  if (!usesNodalNormalDerivatives())
    return;
 
  if (getParam<bool>("interpolate_normals"))
    paramError("interpolate_normals",
               "Nodal-normal derivatives require normalized secondary nodal normals and cannot be "
               "combined with quadrature-point normal interpolation.");
 
  const std::array<std::pair<const MooseVariable *, const char *>, 3> displacement_variables{
      {{_disp_x_var, "disp_x"}, {_disp_y_var, "disp_y"}, {_disp_z_var, "disp_z"}}};
  for (const auto & [variable, parameter_name] : displacement_variables)
    if (variable)
    {
      if (!variable->isNodal())
        paramError(parameter_name,
                   "Nodal-normal derivatives require a nodal displacement variable.");
      if (&variable->sys() != &_sys)
        paramError(parameter_name,
                   "Nodal-normal derivatives require displacement variables in the nonlinear "
                   "system assembled by this contact object.");
    }
}

◆ isAugmentedLagrangianConverged()

bool PenaltyWeightedGapUserObject::isAugmentedLagrangianConverged ( )

overridevirtualinherited

Implements AugmentedLagrangeInterface.

Reimplemented in PenaltyFrictionUserObject.

Definition at line 257 of file PenaltyWeightedGapUserObject.C.

{
  using std::max, std::min, std::abs;
 
  Real max_positive_gap = 0.0;
  Real min_negative_gap = 0.0;
 
  // Get maximum gap to ascertain whether we are converged.
  for (auto & [dof_object, wgap] : _dof_to_weighted_gap)
  {
    const auto gap = physicalGap(wgap);
    {
      // Check condition for nodes that are active
      if (gap < 0 || _dof_to_lagrange_multiplier[dof_object->id()] < 0.0)
      {
        max_positive_gap = max(max_positive_gap, gap);
        min_negative_gap = min(min_negative_gap, gap);
      }
    }
  }
 
  // Communicate the extreme gap values in parallel.
  std::vector<Real> recv;
  if (this->_communicator.rank() == 0)
    recv.resize(this->_communicator.size());
  this->_communicator.gather(
      0, -min_negative_gap > max_positive_gap ? min_negative_gap : max_positive_gap, recv);
 
  if (this->_communicator.rank() == 0)
  {
    min_negative_gap = *std::min_element(recv.begin(), recv.end());
    max_positive_gap = *std::max_element(recv.begin(), recv.end());
 
    // report the gap value with the largest magnitude
    const Real reported_gap =
        -min_negative_gap > max_positive_gap ? min_negative_gap : max_positive_gap;
    if (abs(reported_gap) > _penetration_tolerance)
    {
      mooseInfoRepeated("Penetration tolerance fail max_gap = ",
                        reported_gap,
                        " (gap_tol=",
                        _penetration_tolerance,
                        "). Iteration number is: ",
                        _lagrangian_iteration_number,
                        ".");
      return false;
    }
    else
      mooseInfoRepeated("Penetration tolerance success max_gap = ",
                        reported_gap,
                        " (gap_tol=",
                        _penetration_tolerance,
                        "). Iteration number is: ",
                        _lagrangian_iteration_number,
                        ".");
  }
 
  return true;
}

Referenced by PenaltyFrictionUserObject::isAugmentedLagrangianConverged().

◆ nodalContactPressure()

ADReal PenaltyWeightedGapUserObject::nodalContactPressure ( const Node & node ) const

overridevirtualinherited

Parameters

node	Secondary node of the current secondary lower-dimensional element

Returns: The nodal normal contact pressure at that node, carrying its derivatives

Consumers interpolate each nodal pressure with that node's normal to build the contact traction vector \(\sum_j \Phi_j z_{n,j} \mathbf{n}_j\). This nodal traction interpolation is the discrete transpose of the weighted gap and balances the two sides of the interface; see Popp et al., Comput Methods Appl Mech Eng 264 (2013) 67, Eq. (11). The default reports that this operation is not available.

Reimplemented from WeightedGapUserObject.

Definition at line 215 of file PenaltyWeightedGapUserObject.C.

{
  // A node only enters this map once it has received a weighted-gap contribution, so a node on a
  // partially integrated element may be absent. Absent means no contact pressure, matching how the
  // interpolated pressure in reinit() treats it.
  return findValue(_dof_to_normal_pressure, cast_ptr<const DofObject *>(&node), ADReal(0));
}

◆ nodalCoordinate()

ADPoint WeightedGapUserObject::nodalCoordinate ( const Node & node ) const

protectedinherited

Add displacement derivatives to the node's current coordinate.

Definition at line 242 of file WeightedGapUserObject.C.

{
  mooseAssert(shouldRecordNodalNormalDerivatives(),
              "AD nodal coordinates should only be requested while recording nodal-normal "
              "derivatives.");
  mooseAssert(Moose::doDerivatives(_subproblem, _sys),
              "AD nodal coordinates require derivative recording.");
  ADPoint point = node;
 
  const std::array<std::pair<const MooseVariable *, unsigned int>, 3> displacement_variables{
      {{_disp_x_var, 0}, {_disp_y_var, 1}, {_disp_z_var, 2}}};
  for (const auto & [variable, component] : displacement_variables)
    if (variable)
    {
      const auto sys_num = variable->sys().number();
      const auto var_num = variable->number();
      // Enriched elements can have mapping nodes without independent displacement degrees of
      // freedom. Their positions enter the geometric map, but they contribute no independent
      // displacement derivative.
      if (node.n_dofs(sys_num, var_num))
        Moose::derivInsert(
            point(component).derivatives(), node.dof_number(sys_num, var_num, 0), 1.0);
    }
 
  return point;
}

Referenced by WeightedGapUserObject::contactNormal().

◆ nodalScale()

virtual Real WeightedGapUserObject::nodalScale ( const DofObject * const ) const

inlinevirtualinherited

Node-based Lagrange-multiplier scaling factor kappa_j of Popp et al.

(2013), implemented by LMWeightedGapUserObject. The unit default recovers the unscaled formulation.

Reimplemented in LMWeightedGapUserObject.

Definition at line 108 of file WeightedGapUserObject.h.

108{ return 1; }

Referenced by ComputeWeightedGapLMMechanicalContact::enforceConstraintOnDof(), ComputeFrictionalForceLMMechanicalContact::enforceConstraintOnDof(), and ComputeFrictionalForceLMMechanicalContact::enforceConstraintOnDof3d().

◆ nodalTangentialPressure()

ADReal CohesiveZoneModelBase::nodalTangentialPressure	(	const Node &	node,
		unsigned int	direction
	)		const

overridevirtual

Parameters

node	Secondary node of the current secondary lower-dimensional element
direction	Tangent direction, 0 or 1

Returns: The nodal frictional pressure at that node, expressed in that node's local tangent frame and carrying its derivatives

Consumers interpolate each nodal pressure with that node's tangent to form the frictional traction vector and preserve the discrete transpose relation; see WeightedGapUserObject::nodalContactPressure. The default reports that this operation is not available.

Reimplemented from WeightedVelocitiesUserObject.

Definition at line 250 of file CohesiveZoneModelBase.C.

{
  mooseAssert(direction < 2, "There are at most two tangent directions.");
  return libmesh_map_find(_dof_to_tangential_traction, node.id()).first(direction);
}

◆ normalizeCDivisor()

virtual Real WeightedGapUserObject::normalizeCDivisor	(	const DofObject * const	,
		const Real	covered_normalization
	)		const

inlinevirtualinherited

Parameters

dof	Node/element dof object
covered_normalization	The covered-region normalization for this dof (the second entry of dofToWeightedGap()'s value)

Returns: The divisor to use for normalize_c. Defaults to the covered-region normalization passed in, so normalize_c alone is unaffected; overridden by LMWeightedGapUserObject to return a coverage-independent (full-element) integral instead when node-based scaling is active, so the two mechanisms do not double-count coverage.

Reimplemented in LMWeightedGapUserObject.

Definition at line 122 of file WeightedGapUserObject.h.

  {
    return covered_normalization;
  }

Referenced by ComputeWeightedGapLMMechanicalContact::enforceConstraintOnDof().

◆ normalizeQuantity()

template<class T >

T CohesiveZoneModelBase::normalizeQuantity	(	const std::unordered_map< const DofObject *, T > &	map,
		const Node *const	node
	)

inlineprotected

Normalize mortar quantities (remove mortar integral scaling)

Definition at line 60 of file CohesiveZoneModelBase.h.

  {
    return libmesh_map_find(map, node) / _dof_to_weighted_gap[node].second;
  }

Referenced by BilinearMixedModeCohesiveZoneModel::computeCriticalDisplacementJump(), BilinearMixedModeCohesiveZoneModel::computeCZMTraction(), BilinearMixedModeCohesiveZoneModel::computeEffectiveDisplacementJump(), computeFandR(), BilinearMixedModeCohesiveZoneModel::computeFinalDisplacementJump(), and BilinearMixedModeCohesiveZoneModel::computeModeMixity().

◆ physicalGap()

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

inlineinherited

Compute physical gap from integration gap quantity.

Definition at line 99 of file WeightedGapUserObject.h.

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

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

◆ prepareJumpKinematicQuantities()

void CohesiveZoneModelBase::prepareJumpKinematicQuantities ( )

protectedvirtual

Definition at line 406 of file CohesiveZoneModelBase.C.

{
  // Compute rotation matrix from secondary surface
  // Rotation matrices and local interface displacement jump.
  for (const auto i : make_range(_test->size()))
  {
    const Node * const node = _lower_secondary_elem->node_ptr(i);
 
    // First call does not have maps available
    const bool skip_node = _dof_to_weighted_gap.find(node) == _dof_to_weighted_gap.end();
    if (skip_node)
      continue;
 
    _dof_to_rotation_matrix[node] = CohesiveZoneModelTools::computeReferenceRotation<true>(
        _normals[i], _subproblem.mesh().dimension());
 
    // Every time we used quantities from a map we "denormalize" it from the mortar integral.
    // See normalizing member functions.
    if (auto it = _dof_to_weighted_displacements.find(node);
        it != _dof_to_weighted_displacements.end())
      _dof_to_interface_displacement_jump[node] = it->second;
  }
}

Referenced by reinit().

◆ reinit()

void CohesiveZoneModelBase::reinit ( )

overridevirtual

Reimplemented from PenaltyWeightedGapUserObject.

Definition at line 258 of file CohesiveZoneModelBase.C.

{
  // Normal contact pressure with penalty
  PenaltyWeightedGapUserObject::reinit();
 
  // Compute all rotations that were created as material properties in CZMComputeDisplacementJump
  prepareJumpKinematicQuantities();
 
  // 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);
 
    // Compute the weighted nodal deformation gradient and rotation tensors.
    computeFandR(node);
 
    // The call below is a 'macro' call. Create a utility function or user object for it.
    computeCZMTraction(node);
 
    // Build final traction vector
    computeGlobalTraction(node);
 
    // Compute mechanical contact until end of method.
    const auto penalty_friction = findValue(
        _dof_to_local_penalty_friction, cast_ptr<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->id()];
    auto & [accumulated_slip, old_accumulated_slip] = _dof_to_accumulated_slip[node->id()];
 
    Real normal_lm = -1;
    if (auto it = _dof_to_lagrange_multiplier.find(node->id());
        it != _dof_to_lagrange_multiplier.end())
      normal_lm = it->second;
 
    // Keep active set fixed from second Uzawa loop
    if (normal_lm < -TOLERANCE && normal_pressure > TOLERANCE)
    {
      using std::abs;
 
      const auto & damage = _dof_to_damage[node->id()].first;
 
      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 = std::abs(MetaPhysicL::raw_value(slip_distance).cwiseAbs()(0)) +
                               std::abs(MetaPhysicL::raw_value(slip_distance).cwiseAbs()(1));
 
      const auto slip_norm = (MetaPhysicL::raw_value(slip_distance)).norm();
      const auto friction_limit = 0.4 * _friction_coefficient * damage * abs(normal_pressure);
 
      if (slip_metric > _epsilon_tolerance && penalty_friction * slip_norm > friction_limit)
      {
        inner_iteration_penalty_friction =
            MetaPhysicL::raw_value(friction_limit / (penalty_friction * slip_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 =
          (MetaPhysicL::raw_value(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 + slip_distance.cwiseAbs();
    }
    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);
    }
  }
 
  for (const auto i : index_range(_czm_interpolated_traction))
  {
    _czm_interpolated_traction[i].resize(_qrule_msm->n_points());
 
    for (const auto qp : make_range(_qrule_msm->n_points()))
      _czm_interpolated_traction[i][qp] = 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);
 
    // End of CZM bilinear computations
    auto it = _dof_to_czm_traction.find(node);
    if (it == _dof_to_czm_traction.end())
      continue;
 
    const auto & test_i = (*_test)[i];
    for (const auto qp : make_range(_qrule_msm->n_points()))
      for (const auto idx : index_range(_czm_interpolated_traction))
        _czm_interpolated_traction[idx][qp] += test_i[qp] * (it->second)(idx);
  }
}

◆ selfFinalize()

void PenaltyWeightedGapUserObject::selfFinalize ( )

protectedinherited

Definition at line 146 of file PenaltyWeightedGapUserObject.C.

{
  using std::min;
 
  // 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->id()] : 0.0;
 
    // keep the negative normal pressure (compressive per convention here)
    auto normal_pressure = 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;
  }
}

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

◆ selfInitialize() [1/2]

void PenaltyWeightedGapUserObject::selfInitialize ( )

protectedinherited

Definition at line 132 of file PenaltyWeightedGapUserObject.C.

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

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

◆ selfInitialize() [2/2]

void WeightedVelocitiesUserObject::selfInitialize ( )

protectedinherited

Definition at line 187 of file WeightedVelocitiesUserObject.C.

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

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

◆ selfTimestepSetup() [1/2]

void PenaltyWeightedGapUserObject::selfTimestepSetup ( )

protectedinherited

Definition at line 224 of file PenaltyWeightedGapUserObject.C.

{
  // Let's not clear the LMs for improved performance in
  // nonlinear problems
 
  // timestepSetup() runs once per time step attempt, so on a retried timestep (e.g. --test-restep
  // or a rejected step) the LMs hold the discarded attempt's values; restore the accepted ones
  const bool retried_timestep = _fe_problem.timeStep() == _t_step_old_normal;
  _t_step_old_normal = _fe_problem.timeStep();
  if (retried_timestep)
    _dof_to_lagrange_multiplier = _dof_to_lagrange_multiplier_start;
  else
    _dof_to_lagrange_multiplier_start = _dof_to_lagrange_multiplier;
 
  // 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;
}

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

◆ selfTimestepSetup() [2/2]

void WeightedVelocitiesUserObject::selfTimestepSetup ( )

inlineprotectedinherited

Definition at line 87 of file WeightedVelocitiesUserObject.h.

87{}

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

◆ shouldRecordNodalNormalDerivatives()

bool WeightedGapUserObject::shouldRecordNodalNormalDerivatives ( ) const

inherited

Whether the current assembly records nodal-normal derivatives for this object's system.

Definition at line 97 of file WeightedGapUserObject.C.

{
  return usesNodalNormalDerivatives() && Moose::doDerivatives(_subproblem, _sys);
}

Referenced by WeightedGapUserObject::computeQpIProperties(), WeightedVelocitiesUserObject::computeQpIProperties(), WeightedGapUserObject::contactNormal(), WeightedVelocitiesUserObject::contactTangents(), WeightedGapUserObject::nodalCoordinate(), NormalMortarMechanicalContact::precomputeQpQuantities(), and TangentialMortarMechanicalContact::precomputeQpQuantities().

◆ tangentialTractionBasis()

virtual const VariableTestValue & WeightedVelocitiesUserObject::tangentialTractionBasis ( unsigned int ) const

inlinevirtualinherited

Parameters

direction Tangent direction, 0 or 1

Returns: The basis belonging to the Lagrange multiplier variable representing the tangential contact pressure in that direction. Consumers use this basis for both the interpolation loop and the coefficient lookup performed by nodalTangentialPressure. The default returns tractionBasis() for formulations that share one basis across normal and tangential directions. Formulations with distinct tangential Lagrange multipliers, such as LMWeightedVelocitiesUserObject, override it.

Reimplemented in LMWeightedVelocitiesUserObject.

Definition at line 69 of file WeightedVelocitiesUserObject.h.

  {
    return tractionBasis();
  }

Referenced by TangentialMortarMechanicalContact::precomputeQpQuantities().

◆ test()

const VariableTestValue & CohesiveZoneModelBase::test ( ) const

overrideprotectedvirtual

Returns: The test function associated with the weighted gap

Reimplemented from PenaltyWeightedGapUserObject.

Definition at line 105 of file CohesiveZoneModelBase.C.

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

◆ timestepSetup()

void CohesiveZoneModelBase::timestepSetup ( )

overridevirtual

Reimplemented from PenaltyWeightedGapUserObject.

Definition at line 182 of file CohesiveZoneModelBase.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 = MetaPhysicL::raw_value(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();
 
  // save off damage from the last timestep
  for (auto & map_pr : _dof_to_damage)
  {
    auto & [damage, old_damage] = map_pr.second;
    old_damage = {MetaPhysicL::raw_value(damage)};
    damage = {0.0};
  }
}

◆ tractionBasis()

const VariableTestValue & PenaltyWeightedGapUserObject::tractionBasis ( ) const

overridevirtualinherited

Returns: The basis that interpolates the contact traction. Under the Petrov-Galerkin approach, this is the dual basis; test() is the auxiliary standard basis used by the weighted gap. The default returns test().

Reimplemented from WeightedGapUserObject.

Definition at line 119 of file PenaltyWeightedGapUserObject.C.

{
  // Penalty contact interpolates traction and weighted gap with the same basis.
  return test();
}

◆ updateAugmentedLagrangianMultipliers()

void PenaltyWeightedGapUserObject::updateAugmentedLagrangianMultipliers ( )

overridevirtualinherited

Implements AugmentedLagrangeInterface.

Reimplemented in PenaltyFrictionUserObject.

Definition at line 330 of file PenaltyWeightedGapUserObject.C.

{
  using std::abs;
 
  for (const auto & [dof_object, wgap] : _dof_to_weighted_gap)
  {
    auto & penalty = _dof_to_local_penalty[dof_object];
    if (penalty == 0.0)
      penalty = _penalty;
 
    const auto gap = getNormalGap(static_cast<const Node *>(dof_object));
    auto & lagrange_multiplier = _dof_to_lagrange_multiplier[dof_object->id()];
 
    const auto possible_normalization =
        (_use_physical_gap ? libmesh_map_find(_dof_to_weighted_gap, dof_object).second : 1.0);
 
    // Update penalty (the factor of 1/4 is suggested in the literature, the limit on AL iteration
    // caps the penalty increase)
    // Before we were updating the LM before adapting the penalty factor
    if (lagrange_multiplier + gap * penalty / possible_normalization <= 0)
      lagrange_multiplier += gap * penalty / possible_normalization;
    else
      lagrange_multiplier = 0.0;
 
    const auto previous_gap = _dof_to_previous_gap[dof_object];
    Real eval_tn = 0;
 
    if (_adaptivity_normal == AdaptivityNormalPenalty::SIMPLE)
    {
      if (abs(gap) > 0.25 * abs(previous_gap) && abs(gap) > _penetration_tolerance)
        penalty *= _penalty_multiplier;
    }
    else if (_adaptivity_normal == AdaptivityNormalPenalty::BUSSETTA)
      bussettaAdaptivePenalty(previous_gap, gap, penalty, eval_tn);
 
    // Heuristics to bound the penalty factor
    if (penalty < _penalty)
      penalty = _penalty;
    else if (penalty > _max_penalty_multiplier * _penalty)
      penalty = _max_penalty_multiplier * _penalty;
  }
}

Referenced by PenaltyFrictionUserObject::updateAugmentedLagrangianMultipliers().

◆ usesNodalNormalDerivatives()

virtual bool CohesiveZoneModelBase::usesNodalNormalDerivatives ( ) const

inlineoverridevirtual

Whether this formulation includes displacement derivatives of the secondary nodal normals.

Reimplemented from PenaltyWeightedGapUserObject.

Definition at line 35 of file CohesiveZoneModelBase.h.

35{ return false; }

◆ usesNodalScaling()

virtual bool WeightedGapUserObject::usesNodalScaling ( ) const

inlinevirtualinherited

Returns: Whether node-based Lagrange-multiplier scaling is active

Reimplemented in LMWeightedGapUserObject.

Definition at line 111 of file WeightedGapUserObject.h.

111{ return false; }

◆ validParams()

InputParameters CohesiveZoneModelBase::validParams ( )

static

Definition at line 27 of file CohesiveZoneModelBase.C.

{
  InputParameters params = WeightedVelocitiesUserObject::validParams();
  params += PenaltyWeightedGapUserObject::validParams();
  params.addRequiredCoupledVar("displacements",
                               "The string of displacements suitable for the problem statement");
  params.addClassDescription(
      "Base class for mortar-based cohesive zone model. To handle frictional cohesive interfaces, "
      "it computes the mortar frictional contact forces using a penalty approach.");
  params.addParam<Real>(
      "penalty_friction",
      "The penalty factor for frictional interaction. If not provided, the normal "
      "penalty factor is also used for the frictional problem.");
  params.addParam<Real>(
      "friction_coefficient", 0.0, "The friction coefficient ruling Coulomb friction equations.");
  // Suppress augmented Lagrange parameters. AL implementation for CZM remains to be done.
  params.suppressParameter<Real>("max_penalty_multiplier");
  params.suppressParameter<Real>("penalty_multiplier");
  params.suppressParameter<Real>("penetration_tolerance");
  params.setParameters("penalty", 0.0);
  params.set<bool>("ghost_point_neighbors") = false;
  return params;
}

Referenced by BilinearMixedModeCohesiveZoneModel::validParams().

Member Data Documentation

◆ _3d

bool WeightedVelocitiesUserObject::_3d

protectedinherited

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

Definition at line 138 of file WeightedVelocitiesUserObject.h.

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

◆ _ad_nodal_normals

std::unordered_map<const Node *, ADRealVectorValue> WeightedGapUserObject::_ad_nodal_normals

mutableprotectedinherited

AD nodal normals, each averaged over the secondary faces containing that node.

Definition at line 262 of file WeightedGapUserObject.h.

Referenced by WeightedGapUserObject::contactNormal(), and WeightedGapUserObject::initialize().

◆ _ad_nodal_tangents

std::unordered_map<const Node *, std::array<ADRealVectorValue, 2> > WeightedVelocitiesUserObject::_ad_nodal_tangents

mutableprotectedinherited

AD contact tangent frames, keyed on the secondary node and reused by every test function and quadrature point that interpolates the frictional traction.

Definition at line 142 of file WeightedVelocitiesUserObject.h.

Referenced by WeightedVelocitiesUserObject::contactTangents(), and WeightedVelocitiesUserObject::selfInitialize().

◆ _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 95 of file PenaltyWeightedGapUserObject.h.

Referenced by PenaltyFrictionUserObject::PenaltyFrictionUserObject(), PenaltyWeightedGapUserObject::PenaltyWeightedGapUserObject(), reinit(), PenaltyFrictionUserObject::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 128 of file PenaltyWeightedGapUserObject.h.

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

◆ _contact_pressure

ADVariableValue PenaltyWeightedGapUserObject::_contact_pressure

protectedinherited

The contact pressure on the mortar segument quadrature points.

Definition at line 89 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 239 of file WeightedGapUserObject.h.

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

◆ _czm_interpolated_traction

std::vector<ADVariableValue> CohesiveZoneModelBase::_czm_interpolated_traction

protected

The global traction.

Definition at line 100 of file CohesiveZoneModelBase.h.

Referenced by czmGlobalTraction(), and reinit().

◆ _disp_x_var

const MooseVariable* const WeightedGapUserObject::_disp_x_var

protectedinherited

The x displacement variable.

Definition at line 217 of file WeightedGapUserObject.h.

Referenced by LMWeightedGapUserObject::computeQpINodalScaling(), WeightedGapUserObject::computeQpProperties(), WeightedVelocitiesUserObject::computeQpProperties(), WeightedGapUserObject::dispVar(), WeightedGapUserObject::initialSetup(), LMWeightedGapUserObject::LMWeightedGapUserObject(), WeightedGapUserObject::nodalCoordinate(), PenaltyWeightedGapUserObject::PenaltyWeightedGapUserObject(), test(), PenaltyFrictionUserObject::test(), and PenaltyWeightedGapUserObject::test().

◆ _disp_y_var

const MooseVariable* const WeightedGapUserObject::_disp_y_var

protectedinherited

The y displacement variable.

Definition at line 219 of file WeightedGapUserObject.h.

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

◆ _disp_z_var

const MooseVariable* const WeightedGapUserObject::_disp_z_var

protectedinherited

The z displacement variable.

Definition at line 223 of file WeightedGapUserObject.h.

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

◆ _dof_to_accumulated_slip

std::unordered_map<dof_id_type, std::pair<ADTwoVector, TwoVector> >& CohesiveZoneModelBase::_dof_to_accumulated_slip

protected

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

Definition at line 112 of file CohesiveZoneModelBase.h.

Referenced by reinit(), and timestepSetup().

◆ _dof_to_czm_traction

std::unordered_map<const DofObject *, ADRealVectorValue> CohesiveZoneModelBase::_dof_to_czm_traction

protected

Total Lagrangian stress to be applied on CZM interface.

Definition at line 133 of file CohesiveZoneModelBase.h.

Referenced by BilinearMixedModeCohesiveZoneModel::computeCZMTraction(), computeGlobalTraction(), initialize(), and reinit().

◆ _dof_to_damage

std::unordered_map<dof_id_type, std::pair<ADReal, Real> >& CohesiveZoneModelBase::_dof_to_damage

protected

Damage values (pair of current and old) on CZM interface.

Definition at line 139 of file CohesiveZoneModelBase.h.

Referenced by BilinearMixedModeCohesiveZoneModel::computeCZMTraction(), BilinearMixedModeCohesiveZoneModel::computeDamage(), BilinearMixedModeCohesiveZoneModel::getCohesiveDamage(), reinit(), and timestepSetup().

◆ _dof_to_displacement_jump

std::unordered_map<const DofObject *, ADRealVectorValue> CohesiveZoneModelBase::_dof_to_displacement_jump

protected

Displacement jump on CZM interface.

Definition at line 136 of file CohesiveZoneModelBase.h.

◆ _dof_to_F

std::unordered_map<const DofObject *, ADRankTwoTensor> CohesiveZoneModelBase::_dof_to_F

protected

Map from degree of freedom to secondary, interpolated deformation gradient tensor.

Definition at line 88 of file CohesiveZoneModelBase.h.

Referenced by computeFandR(), computeQpIProperties(), finalize(), and initialize().

◆ _dof_to_F_neighbor

std::unordered_map<const DofObject *, ADRankTwoTensor> CohesiveZoneModelBase::_dof_to_F_neighbor

protected

Map from degree of freedom to neighbor, interpolated deformation gradient tensor.

Definition at line 91 of file CohesiveZoneModelBase.h.

Referenced by computeFandR(), computeQpIProperties(), finalize(), and initialize().

◆ _dof_to_frictional_lagrange_multipliers

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

protected

Map from degree of freedom to augmented lagrange multiplier.

Definition at line 124 of file CohesiveZoneModelBase.h.

Referenced by reinit(), and timestepSetup().

◆ _dof_to_interface_displacement_jump

std::unordered_map<const DofObject *, ADRealVectorValue> CohesiveZoneModelBase::_dof_to_interface_displacement_jump

protected

Map from degree of freedom to local displacement jump.

Definition at line 79 of file CohesiveZoneModelBase.h.

◆ _dof_to_interface_F

std::unordered_map<const DofObject *, ADRankTwoTensor> CohesiveZoneModelBase::_dof_to_interface_F

protected

Map from degree of freedom to interface deformation gradient tensor.

Definition at line 94 of file CohesiveZoneModelBase.h.

Referenced by computeFandR(), and initialize().

◆ _dof_to_interface_R

std::unordered_map<const DofObject *, ADRankTwoTensor> CohesiveZoneModelBase::_dof_to_interface_R

protected

Map from degree of freedom to interface rotation tensor.

Definition at line 97 of file CohesiveZoneModelBase.h.

Referenced by computeFandR(), and initialize().

◆ _dof_to_lagrange_multiplier

std::unordered_map<dof_id_type, Real>& PenaltyWeightedGapUserObject::_dof_to_lagrange_multiplier

protectedinherited

Map from node id to augmented lagrange multiplier.

Keyed by id so that it can be stored as restartable data, since the multipliers are kept between timesteps.

Definition at line 102 of file PenaltyWeightedGapUserObject.h.

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

◆ _dof_to_lagrange_multiplier_start

std::unordered_map<dof_id_type, Real> PenaltyWeightedGapUserObject::_dof_to_lagrange_multiplier_start

protectedinherited

Augmented lagrange multipliers at the start of the current timestep, restored when the timestep is retried so the retry starts from the accepted state.

Definition at line 106 of file PenaltyWeightedGapUserObject.h.

Referenced by PenaltyWeightedGapUserObject::selfTimestepSetup().

◆ _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 113 of file PenaltyWeightedGapUserObject.h.

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

◆ _dof_to_local_penalty_friction

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

protected

Map from degree of freedom to local friction penalty value.

Definition at line 127 of file CohesiveZoneModelBase.h.

Referenced by reinit(), and timestepSetup().

◆ _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 92 of file PenaltyWeightedGapUserObject.h.

Referenced by PenaltyWeightedGapUserObject::getActiveSetState(), PenaltyWeightedGapUserObject::getNormalContactPressure(), PenaltyFrictionUserObject::isAugmentedLagrangianConverged(), PenaltyWeightedGapUserObject::nodalContactPressure(), reinit(), PenaltyFrictionUserObject::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 116 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 93 of file WeightedVelocitiesUserObject.h.

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

◆ _dof_to_rotation_matrix

std::unordered_map<const DofObject *, ADRankTwoTensor> CohesiveZoneModelBase::_dof_to_rotation_matrix

protected

*** Kinematics/displacement jump quantities *** Map from degree of freedom to rotation matrix

Definition at line 76 of file CohesiveZoneModelBase.h.

Referenced by computeFandR(), computeGlobalTraction(), initialize(), and prepareJumpKinematicQuantities().

◆ _dof_to_step_slip

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

protected

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

Definition at line 109 of file CohesiveZoneModelBase.h.

Referenced by timestepSetup().

◆ _dof_to_tangential_traction

std::unordered_map<dof_id_type, std::pair<ADTwoVector, TwoVector> >& CohesiveZoneModelBase::_dof_to_tangential_traction

protected

Definition at line 115 of file CohesiveZoneModelBase.h.

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

◆ _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 251 of file WeightedGapUserObject.h.

Referenced by computeQpIProperties(), finalize(), WeightedGapUserObject::initialize(), and 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 248 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 90 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 119 of file PenaltyWeightedGapUserObject.h.

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

◆ _dt_old

Real PenaltyWeightedGapUserObject::_dt_old

protectedinherited

previous timestep size

Definition at line 122 of file PenaltyWeightedGapUserObject.h.

Referenced by PenaltyWeightedGapUserObject::selfTimestepSetup().

◆ _epsilon_tolerance

const Real CohesiveZoneModelBase::_epsilon_tolerance

protected

Tolerance to avoid NaN/Inf in automatic differentiation operations.

Definition at line 130 of file CohesiveZoneModelBase.h.

Referenced by BilinearMixedModeCohesiveZoneModel::computeCriticalDisplacementJump(), BilinearMixedModeCohesiveZoneModel::computeEffectiveDisplacementJump(), BilinearMixedModeCohesiveZoneModel::computeFinalDisplacementJump(), BilinearMixedModeCohesiveZoneModel::computeModeMixity(), and reinit().

◆ _F_interpolation

ADRankTwoTensor CohesiveZoneModelBase::_F_interpolation

protected

Deformation gradient for interpolation.

Definition at line 82 of file CohesiveZoneModelBase.h.

Referenced by computeQpIProperties(), and computeQpProperties().

◆ _F_neighbor_interpolation

ADRankTwoTensor CohesiveZoneModelBase::_F_neighbor_interpolation

protected

Deformation gradient for interpolation of the neighbor projection.

Definition at line 85 of file CohesiveZoneModelBase.h.

Referenced by computeQpIProperties(), and computeQpProperties().

◆ _fe_problem

FEProblemBase& WeightedGapUserObject::_fe_problem

protectedinherited

The base finite element problem.

Definition at line 205 of file WeightedGapUserObject.h.

Referenced by PenaltyWeightedGapUserObject::selfTimestepSetup(), and PenaltyFrictionUserObject::timestepSetup().

◆ _friction_coefficient

const Real CohesiveZoneModelBase::_friction_coefficient

protected

The friction coefficient.

Definition at line 106 of file CohesiveZoneModelBase.h.

Referenced by reinit().

◆ _frictional_contact_traction_one

ADVariableValue CohesiveZoneModelBase::_frictional_contact_traction_one

protected

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

Definition at line 118 of file CohesiveZoneModelBase.h.

Referenced by contactTangentialPressureDirOne(), and reinit().

◆ _frictional_contact_traction_two

ADVariableValue CohesiveZoneModelBase::_frictional_contact_traction_two

protected

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

Definition at line 121 of file CohesiveZoneModelBase.h.

Referenced by contactTangentialPressureDirTwo(), and reinit().

◆ _grad_disp

std::vector<const GenericVariableGradient<true> *> CohesiveZoneModelBase::_grad_disp

protected

Coupled displacement gradients.

Definition at line 69 of file CohesiveZoneModelBase.h.

Referenced by computeQpProperties().

◆ _grad_disp_neighbor

std::vector<const GenericVariableGradient<true> *> CohesiveZoneModelBase::_grad_disp_neighbor

protected

Coupled displacement and neighbor displacement gradient.

Definition at line 72 of file CohesiveZoneModelBase.h.

Referenced by computeQpProperties().

◆ _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 221 of file WeightedGapUserObject.h.

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

◆ _i

unsigned int WeightedGapUserObject::_i = 0

protectedinherited

Test function index.

Definition at line 273 of file WeightedGapUserObject.h.

Referenced by LMWeightedGapUserObject::computeQpINodalScaling(), WeightedGapUserObject::computeQpIProperties(), BilinearMixedModeCohesiveZoneModel::computeQpIProperties(), computeQpIProperties(), WeightedVelocitiesUserObject::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 267 of file WeightedGapUserObject.h.

Referenced by WeightedVelocitiesUserObject::computeQpIProperties().

◆ _lagrangian_iteration_number

const unsigned int& PenaltyWeightedGapUserObject::_lagrangian_iteration_number

protectedinherited

Definition at line 97 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 131 of file PenaltyWeightedGapUserObject.h.

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

◆ _ndisp

const unsigned int CohesiveZoneModelBase::_ndisp

protected

Number of displacement components.

Definition at line 66 of file CohesiveZoneModelBase.h.

◆ _no_iterations

const unsigned int PenaltyWeightedGapUserObject::_no_iterations = 0

staticprotectedinherited

Definition at line 96 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 214 of file WeightedGapUserObject.h.

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

◆ _normalization_ptr

const Real* WeightedGapUserObject::_normalization_ptr = nullptr

protectedinherited

Definition at line 255 of file WeightedGapUserObject.h.

◆ _penalty

const Real PenaltyWeightedGapUserObject::_penalty

protectedinherited

The penalty factor.

Definition at line 80 of file PenaltyWeightedGapUserObject.h.

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

◆ _penalty_friction

const Real CohesiveZoneModelBase::_penalty_friction

protected

The penalty factor for the frictional constraints.

Definition at line 103 of file CohesiveZoneModelBase.h.

Referenced by reinit(), and timestepSetup().

◆ _penalty_multiplier

const Real PenaltyWeightedGapUserObject::_penalty_multiplier

protectedinherited

penalty growth factor for augmented Lagrange

Definition at line 83 of file PenaltyWeightedGapUserObject.h.

Referenced by PenaltyWeightedGapUserObject::updateAugmentedLagrangianMultipliers().

◆ _penetration_tolerance

const Real PenaltyWeightedGapUserObject::_penetration_tolerance

protectedinherited

penetration tolerance for augmented Lagrange contact

Definition at line 86 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 228 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 232 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 236 of file WeightedGapUserObject.h.

◆ _primary_var

MooseVariableField<Real>& WeightedVelocitiesUserObject::_primary_var

protectedinherited

Reference to the primary variable.

Definition at line 117 of file WeightedVelocitiesUserObject.h.

◆ _primary_x_dot

const ADVariableValue& WeightedVelocitiesUserObject::_primary_x_dot

protectedinherited

x-velocity on the primary face

Definition at line 123 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 129 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 135 of file WeightedVelocitiesUserObject.h.

◆ _qp

unsigned int WeightedGapUserObject::_qp = 0

protectedinherited

Quadrature point index for the mortar segments.

Definition at line 270 of file WeightedGapUserObject.h.

Referenced by LMWeightedGapUserObject::computeQpINodalScaling(), WeightedGapUserObject::computeQpIProperties(), BilinearMixedModeCohesiveZoneModel::computeQpIProperties(), computeQpIProperties(), WeightedVelocitiesUserObject::computeQpIProperties(), WeightedGapUserObject::computeQpProperties(), BilinearMixedModeCohesiveZoneModel::computeQpProperties(), computeQpProperties(), WeightedVelocitiesUserObject::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 245 of file WeightedGapUserObject.h.

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

◆ _qp_factor

Real WeightedGapUserObject::_qp_factor

protectedinherited

The value of the LM at the current quadrature point.

Definition at line 211 of file WeightedGapUserObject.h.

Referenced by LMWeightedGapUserObject::computeQpINodalScaling(), 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 208 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 242 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 102 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 108 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 99 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 105 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 226 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 230 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 234 of file WeightedGapUserObject.h.

◆ _secondary_var

MooseVariableField<Real>& WeightedVelocitiesUserObject::_secondary_var

protectedinherited

Reference to the secondary variable.

Definition at line 114 of file WeightedVelocitiesUserObject.h.

◆ _secondary_x_dot

const ADVariableValue& WeightedVelocitiesUserObject::_secondary_x_dot

protectedinherited

x-velocity on the secondary face

Definition at line 120 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 126 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 132 of file WeightedVelocitiesUserObject.h.

◆ _sys

SystemBase& WeightedVelocitiesUserObject::_sys

protectedinherited

Reference to the EquationSystem object.

Definition at line 111 of file WeightedVelocitiesUserObject.h.

Referenced by LMWeightedVelocitiesUserObject::nodalTangentialPressure().

◆ _t_step_old_normal

int& PenaltyWeightedGapUserObject::_t_step_old_normal

protectedinherited

The timestep index at which _dof_to_lagrange_multiplier_start was last saved.

Used to detect a retried timestep, since timestepSetup() runs once per time step attempt.

Definition at line 110 of file PenaltyWeightedGapUserObject.h.

Referenced by PenaltyWeightedGapUserObject::selfTimestepSetup().

◆ _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 96 of file WeightedVelocitiesUserObject.h.

96{{nullptr, nullptr}};

◆ _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 259 of file WeightedGapUserObject.h.

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

◆ _use_physical_gap

const bool PenaltyWeightedGapUserObject::_use_physical_gap

protectedinherited

Use scaled or physical gap.

Definition at line 125 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 254 of file WeightedGapUserObject.h.

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

contact/include/userobjects/CohesiveZoneModelBase.h
contact/src/userobjects/CohesiveZoneModelBase.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

Private Types

Private Member Functions

Static Private Member Functions

Private Attributes

Detailed Description

Member Enumeration Documentation

◆ AdaptivityNormalPenalty

Constructor & Destructor Documentation

◆ CohesiveZoneModelBase()

Member Function Documentation

◆ adaptiveNormalPenalty()

◆ addNodalNormalCoupling()

◆ adPhysicalGap()

◆ augmentedLagrangianSetup()

◆ bussettaAdaptivePenalty()

◆ computeCZMTraction()

◆ computeDamage()

◆ computeFandR()

◆ computeGlobalTraction()

◆ computeQpIProperties()

◆ computeQpProperties()

◆ constrainedByOwner()

◆ contactNormal()

◆ contactPressure()

◆ contactTangentialPressureDirOne()

◆ contactTangentialPressureDirTwo()

◆ contactTangents()

◆ czmGlobalTraction()

◆ dispVar()

◆ dofToWeightedGap()

◆ dofToWeightedVelocities()

◆ execute()

◆ finalize()

◆ findValue()

◆ getAccumulatedSlip()

◆ getActiveSetState()

◆ getDeltaTangentialLagrangeMultiplier()

◆ getFrictionalContactPressure()

◆ getNormalContactPressure()

◆ getNormalGap()

◆ getNormalLagrangeMultiplier()

◆ getTangentialVelocity()

◆ initialize()

◆ initialSetup()

◆ isAugmentedLagrangianConverged()

◆ nodalContactPressure()

◆ nodalCoordinate()

◆ nodalScale()

◆ nodalTangentialPressure()

◆ normalizeCDivisor()

◆ normalizeQuantity()

◆ physicalGap()

◆ prepareJumpKinematicQuantities()

◆ reinit()

◆ selfFinalize()

◆ selfInitialize() [1/2]

◆ selfInitialize() [2/2]

◆ selfTimestepSetup() [1/2]

◆ selfTimestepSetup() [2/2]

◆ shouldRecordNodalNormalDerivatives()

◆ tangentialTractionBasis()

◆ test()

◆ timestepSetup()

◆ tractionBasis()

◆ updateAugmentedLagrangianMultipliers()

◆ usesNodalNormalDerivatives()

◆ usesNodalScaling()

◆ validParams()

Member Data Documentation

◆ _3d