https://mooseframework.inl.gov
Loading...
Searching...
No Matches
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 | List of all members
BilinearMixedModeCohesiveZoneModel Class Reference

User object that computes bilinear mixed mode traction separation law. More...

#include <BilinearMixedModeCohesiveZoneModel.h>

Inheritance diagram for BilinearMixedModeCohesiveZoneModel:
[legend]

Public Types

typedef DataFileName DataFileParameterType
 

Public Member Functions

 BilinearMixedModeCohesiveZoneModel (const InputParameters &parameters)
 
virtual void initialize () override
 
virtual void finalize () override
 
Real getModeMixityRatio (const Node *const node) const
 
Real getCohesiveDamage (const Node *const node) const
 
Real getLocalDisplacementNormal (const Node *const node) const
 
Real getLocalDisplacementTangential (const Node *const node) const
 
Real getLocalDisplacementTangentialOne (const Node *const node) const
 
Real getLocalDisplacementTangentialTwo (const Node *const node) const
 
Real getLocalDisplacementTangentialEffective (const Node *const node) const
 
Real getCohesiveTractionNormal (const Node *const node) const
 
Real getCohesiveTractionTangentialMagnitude (const Node *const node) const
 
Real getCohesiveTractionTangentialOne (const Node *const node) const
 
Real getCohesiveTractionTangentialTwo (const Node *const node) const
 
Real getCohesiveTractionEffective (const Node *const node) const
 
virtual const ADVariableValue & contactTangentialPressureDirOne () const override
 
virtual const ADVariableValue & contactTangentialPressureDirTwo () const 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  MixedModeCriterion { POWER_LAW , BK }
 Mixed-mode propagation criterion. More...
 
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 bool constrainedByOwner () const override
 
virtual void computeModeMixity (const Node *const node)
 
virtual void computeCriticalDisplacementJump (const Node *const node)
 
virtual void computeFinalDisplacementJump (const Node *const node)
 
virtual void computeEffectiveDisplacementJump (const Node *const node)
 
virtual void computeDamage (const Node *const node) override
 
virtual void computeCZMTraction (const Node *const node) override
 Encapsulate the CZM constitutive behavior.
 
virtual const VariableTestValue & test () const override
 
virtual void prepareJumpKinematicQuantities ()
 
virtual void computeFandR (const Node *const node)
 
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 selfInitialize ()
 
void selfFinalize ()
 
void selfTimestepSetup ()
 
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)
 

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 bool _set_compressive_traction_to_zero
 Zero compressive traction.
 
std::unordered_map< const DofObject *, ADReal > _dof_to_mode_mixity_ratio
 Map from degree of freedom to mode mixity ratio (AD needed?)
 
std::unordered_map< const DofObject *, RealVectorValue > _dof_to_local_czm_traction
 Map from degree of freedom to the local-frame cohesive traction.
 
const MaterialProperty< Real > & _normal_strength
 The normal strength material property.
 
const MaterialProperty< Real > & _shear_strength
 The shear strength material property.
 
const MaterialProperty< Real > & _GI_c
 Fracture parameter mode I.
 
const MaterialProperty< Real > & _GII_c
 Fracture parameter mode II.
 
ADReal _normal_strength_interpolation
 Interpolated value of normal_strength.
 
ADReal _shear_strength_interpolation
 Interpolated value of shear_strength.
 
ADReal _GI_c_interpolation
 Interpolated value of fracture paramter mode I.
 
ADReal _GII_c_interpolation
 Interpolated value of fracture paramter mode II.
 
const Real _penalty_stiffness_czm
 
enum BilinearMixedModeCohesiveZoneModel::MixedModeCriterion _mix_mode_criterion
 
const Real _power_law_parameter
 Power law parameter for bilinear traction model.
 
const Real _viscosity
 Viscosity for damage model.
 
const Real _regularization_alpha
 Parameter for the regularization of the Macaulay bracket.
 
std::unordered_map< const DofObject *, ADReal > _dof_to_normal_strength
 
std::unordered_map< const DofObject *, ADReal > _dof_to_shear_strength
 
std::unordered_map< const DofObject *, ADReal > _dof_to_GI_c
 
std::unordered_map< const DofObject *, ADReal > _dof_to_GII_c
 
std::unordered_map< const DofObject *, ADReal > _dof_to_delta_initial
 
std::unordered_map< const DofObject *, ADReal > _dof_to_delta_final
 
std::unordered_map< const DofObject *, ADReal > _dof_to_delta_max
 
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< const DofObject *, Real > _dof_to_lagrange_multiplier
 Map from degree of freedom to augmented lagrange multiplier.
 
std::unordered_map< const DofObject *, 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
 
SystemBase & _sys
 Reference to the EquationSystem object.
 
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.
 
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

User object that computes bilinear mixed mode traction separation law.

Definition at line 18 of file BilinearMixedModeCohesiveZoneModel.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 133 of file PenaltyWeightedGapUserObject.h.

◆ MixedModeCriterion

Mixed-mode propagation criterion.

Enumerator
POWER_LAW 
BK 

Definition at line 98 of file BilinearMixedModeCohesiveZoneModel.h.

Constructor & Destructor Documentation

◆ BilinearMixedModeCohesiveZoneModel()

BilinearMixedModeCohesiveZoneModel::BilinearMixedModeCohesiveZoneModel ( const InputParameters &  parameters)

Definition at line 71 of file BilinearMixedModeCohesiveZoneModel.C.

77 _set_compressive_traction_to_zero(getParam<bool>("set_compressive_traction_to_zero")),
78 _normal_strength(getMaterialProperty<Real>("normal_strength")),
79 _shear_strength(getMaterialProperty<Real>("shear_strength")),
80 _GI_c(getMaterialProperty<Real>("GI_c")),
81 _GII_c(getMaterialProperty<Real>("GII_c")),
82 _penalty_stiffness_czm(getParam<Real>("penalty_stiffness")),
83 _mix_mode_criterion(getParam<MooseEnum>("mixed_mode_criterion").getEnum<MixedModeCriterion>()),
84 _power_law_parameter(getParam<Real>("power_law_parameter")),
85 _viscosity(getParam<Real>("viscosity")),
86 _regularization_alpha(getParam<Real>("regularization_alpha"))
87{
88
89 // Checks for bilinear traction input.
90 if (!isParamValid("GI_c") || !isParamValid("GII_c") || !isParamValid("normal_strength") ||
91 !isParamValid("shear_strength") || !isParamValid("power_law_parameter") ||
92 !isParamValid("penalty_stiffness"))
93 paramError("GI_c",
94 "The CZM bilinear mixed mode traction parameters GI_c, GII_c, normal_strength, "
95 "shear_strength, and power_law_parameter are required. Revise your input and add "
96 "those parameters if you want to use the bilinear mixed mode traction model. ");
97}
const bool _set_compressive_traction_to_zero
Zero compressive traction.
const MaterialProperty< Real > & _normal_strength
The normal strength material property.
const Real _power_law_parameter
Power law parameter for bilinear traction model.
const Real _regularization_alpha
Parameter for the regularization of the Macaulay bracket.
const MaterialProperty< Real > & _GI_c
Fracture parameter mode I.
const MaterialProperty< Real > & _GII_c
Fracture parameter mode II.
const Real _viscosity
Viscosity for damage model.
const MaterialProperty< Real > & _shear_strength
The shear strength material property.
Base class for mortar-based cohesive zone model.
const InputParameters & parameters() const
void paramError(const std::string &param, Args... args) const
bool isParamValid(const std::string &name) const
User object for computing weighted gaps and contact pressure for penalty based mortar constraints.
Creates dof object to weighted gap map.
Creates dof object to weighted tangential velocities map.

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

398{
399 using std::abs, std::max, std::sqrt;
400
401 const bool condition_one = abs(abs(previous_gap / gap) - 1.0) < 0.01;
402 const bool condition_two = abs(gap) > 1.01 * abs(previous_gap);
403
404 if (previous_gap * gap < 0)
405 {
406 if (previous_gap > _penetration_tolerance)
407 penalty = abs(penalty * previous_gap / gap * (abs(gap) + _penetration_tolerance) /
408 (gap - previous_gap));
409 else
410 penalty = abs(penalty * previous_gap / 10.0 / gap);
411 }
412 else if (abs(gap) > _penetration_tolerance)
413 {
414 if (abs(gap - previous_gap) >
415 max(gap / 10.0, max(previous_gap / 10.0, 5.0 * _penetration_tolerance)))
416 penalty *= 10.0;
417 else if (condition_one && gap < 10.0 * _penetration_tolerance)
418 penalty *= MathUtils::pow(sqrt(abs(gap) / _penetration_tolerance - 1.0) + 1.0, 2);
419 else if (condition_two)
420 penalty *= 10.0 * previous_gap / gap;
421 else
422 penalty *= sqrt(abs(gap) / _penetration_tolerance - 1.0) + 1.0;
423 }
424}
const Real _penetration_tolerance
penetration tolerance for augmented Lagrange contact
auto max(const L &left, const R &right)
CTSub CT_OPERATOR_BINARY CTMul CTCompareLess CTCompareGreater CTCompareEqual _arg template * sqrt(_arg)) *_arg.template D< dtag >()) CT_SIMPLE_UNARY_FUNCTION(tanh
MetaPhysicL::DualNumber< V, D, asd > abs(const MetaPhysicL::DualNumber< V, D, asd > &a)
T pow(T x, int e)

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.

41{
42 params.set<bool>("ghost_point_neighbors") = true;
43 const auto configure_point_neighbors =
44 [](const InputParameters & obj_params, InputParameters & rm_params)
45 {
46 rm_params.set<bool>("use_displaced_mesh") = obj_params.get<bool>("use_displaced_mesh");
47 rm_params.set<BoundaryName>("secondary_boundary") =
48 obj_params.get<BoundaryName>("secondary_boundary");
49 rm_params.set<BoundaryName>("primary_boundary") =
50 obj_params.get<BoundaryName>("primary_boundary");
51 rm_params.set<SubdomainName>("secondary_subdomain") =
52 obj_params.get<SubdomainName>("secondary_subdomain");
53 rm_params.set<SubdomainName>("primary_subdomain") =
54 obj_params.get<SubdomainName>("primary_subdomain");
55 // Objects that disable ghost_point_neighbors do not include nodal-normal derivatives.
56 rm_params.set<bool>("ghost_point_neighbors") = obj_params.get<bool>("ghost_point_neighbors");
57 };
58 params.addRelationshipManager("AugmentSparsityOnInterface",
59 Moose::RelationshipManagerType::COUPLING,
60 configure_point_neighbors);
61}
std::vector< std::pair< R1, R2 > > get(const std::string &param1, const std::string &param2) const
void addRelationshipManager(const std::string &name, Moose::RelationshipManagerType rm_type, Moose::RelationshipManagerInputParameterCallback input_parameter_callback=nullptr)
T & set(const std::string &name, bool quiet_mode=false)

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

317{
318 // Loop over all nodes for which a gap has been computed
319 for (auto & [dof_object, wgap] : _dof_to_weighted_gap)
320 {
321 const Real gap = physicalGap(wgap);
322 // Store previous augmented lagrange iteration gap
323 _dof_to_previous_gap[dof_object] = gap;
324 }
325}
std::unordered_map< const DofObject *, Real > _dof_to_previous_gap
Map from degree of freedom to previous AL iteration gap values.
std::unordered_map< const DofObject *, std::pair< ADReal, Real > > _dof_to_weighted_gap
A map from node to weighted gap and normalization (if requested)
Real physicalGap(const std::pair< ADReal, Real > &gap) const
Compute physical gap from integration gap quantity.
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real

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

376{
377 // Positive gaps means no contact
378 if (previous_gap > 0.0)
379 {
380 penalty = _penalty;
381 eval_tn = 0.0;
382 }
383 else
384 {
385 if (previous_gap * gap < 0)
386 eval_tn = 0.0;
387 else
388 eval_tn = penalty * previous_gap;
389
390 adaptiveNormalPenalty(previous_gap, gap, penalty);
391 }
392}
const Real _penalty
The penalty factor.
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 m...

Referenced by PenaltyWeightedGapUserObject::updateAugmentedLagrangianMultipliers().

◆ computeCriticalDisplacementJump()

void BilinearMixedModeCohesiveZoneModel::computeCriticalDisplacementJump ( const Node *const  node)
protectedvirtual

Definition at line 201 of file BilinearMixedModeCohesiveZoneModel.C.

202{
203 using std::sqrt;
204
205 const auto interface_displacement_jump =
207
208 const auto mixity_ratio = libmesh_map_find(_dof_to_mode_mixity_ratio, node);
209
210 const auto delta_normal_knot =
212 const auto delta_shear_knot =
214
215 _dof_to_delta_initial[node] = delta_shear_knot;
216
217 if (interface_displacement_jump(0) > _epsilon_tolerance)
218 {
219 const auto delta_mixed = sqrt(delta_shear_knot * delta_shear_knot +
220 Utility::pow<2>(mixity_ratio * delta_normal_knot));
221
222 _dof_to_delta_initial[node] = delta_normal_knot * delta_shear_knot *
223 sqrt(1.0 + mixity_ratio * mixity_ratio) / delta_mixed;
224 }
225}
std::unordered_map< const DofObject *, ADReal > _dof_to_shear_strength
std::unordered_map< const DofObject *, ADReal > _dof_to_normal_strength
std::unordered_map< const DofObject *, ADReal > _dof_to_mode_mixity_ratio
Map from degree of freedom to mode mixity ratio (AD needed?)
std::unordered_map< const DofObject *, ADReal > _dof_to_delta_initial
std::unordered_map< const DofObject *, ADRealVectorValue > _dof_to_interface_displacement_jump
Map from degree of freedom to local displacement jump.
const Real _epsilon_tolerance
Tolerance to avoid NaN/Inf in automatic differentiation operations.
T normalizeQuantity(const std::unordered_map< const DofObject *, T > &map, const Node *const node)
Normalize mortar quantities (remove mortar integral scaling)

Referenced by computeCZMTraction().

◆ computeCZMTraction()

void BilinearMixedModeCohesiveZoneModel::computeCZMTraction ( const Node *const  node)
overrideprotectedvirtual

Encapsulate the CZM constitutive behavior.

Implements CohesiveZoneModelBase.

Definition at line 145 of file BilinearMixedModeCohesiveZoneModel.C.

146{
147 using std::max, std::min;
148
149 // First call does not have maps available
150 const bool return_boolean = _dof_to_weighted_gap.find(node) == _dof_to_weighted_gap.end();
151 if (return_boolean)
152 return;
153
154 computeModeMixity(node);
158 computeDamage(node);
159
160 // Split displacement jump into active and inactive parts
161 const auto interface_displacement_jump =
163
164 const ADRealVectorValue delta_active(max(interface_displacement_jump(0), 0.0),
165 interface_displacement_jump(1),
166 interface_displacement_jump(2));
167 const ADRealVectorValue delta_inactive(min(interface_displacement_jump(0), 0.0), 0.0, 0.0);
168
169 // This traction vector is local at this point.
171 -(1.0 - libmesh_map_find(_dof_to_damage, node->id()).first) * _penalty_stiffness_czm *
172 delta_active -
174
175 // Save local-frame traction before the base class overwrites this map with global traction.
176 // Raw values only: this map feeds auxiliary output, so derivatives are not needed.
178}
virtual void computeDamage(const Node *const node) override
virtual void computeFinalDisplacementJump(const Node *const node)
virtual void computeEffectiveDisplacementJump(const Node *const node)
std::unordered_map< const DofObject *, RealVectorValue > _dof_to_local_czm_traction
Map from degree of freedom to the local-frame cohesive traction.
virtual void computeModeMixity(const Node *const node)
virtual void computeCriticalDisplacementJump(const Node *const node)
std::unordered_map< const DofObject *, ADRealVectorValue > _dof_to_czm_traction
Total Lagrangian stress to be applied 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.
auto min(const L &left, const R &right)
auto raw_value(const Eigen::Map< T > &in)

◆ computeDamage()

void BilinearMixedModeCohesiveZoneModel::computeDamage ( const Node *const  node)
overrideprotectedvirtual

Implements CohesiveZoneModelBase.

Definition at line 285 of file BilinearMixedModeCohesiveZoneModel.C.

286{
287 const auto delta_max = libmesh_map_find(_dof_to_delta_max, node);
288 const auto delta_initial = libmesh_map_find(_dof_to_delta_initial, node);
289 const auto delta_final = libmesh_map_find(_dof_to_delta_final, node);
290
291 auto & [damage, damage_old] = _dof_to_damage[node->id()];
292 if (delta_max < delta_initial)
293 damage = 0;
294 else if (delta_max > delta_final)
295 damage = 1.0;
296 else
297 damage = delta_final * (delta_max - delta_initial) / delta_max / (delta_final - delta_initial);
298
299 if (damage < damage_old)
300 // Irreversibility
301 damage = damage_old;
302
303 // Viscous regularization
304 damage = (damage + _viscosity * damage_old / _dt) / (_viscosity / _dt + 1.0);
305}
std::unordered_map< const DofObject *, ADReal > _dof_to_delta_final
std::unordered_map< const DofObject *, ADReal > _dof_to_delta_max
const Real & _dt
Current delta t... or timestep size.

Referenced by computeCZMTraction().

◆ computeEffectiveDisplacementJump()

void BilinearMixedModeCohesiveZoneModel::computeEffectiveDisplacementJump ( const Node *const  node)
protectedvirtual

Definition at line 268 of file BilinearMixedModeCohesiveZoneModel.C.

269{
270 using std::sqrt;
271
272 const auto interface_displacement_jump =
274
275 const auto delta_normal_pos =
276 MathUtils::regularizedHeavyside(interface_displacement_jump(0), _regularization_alpha) *
277 interface_displacement_jump(0);
278
279 _dof_to_delta_max[node] = sqrt(Utility::pow<2>(interface_displacement_jump(1)) +
280 Utility::pow<2>(interface_displacement_jump(2)) +
281 Utility::pow<2>(delta_normal_pos) + _epsilon_tolerance);
282}
T regularizedHeavyside(const T &x, Real smoothing_length)

Referenced by computeCZMTraction().

◆ computeFandR()

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

Definition at line 144 of file CohesiveZoneModelBase.C.

145{
146 // First call does not have maps available
147 const bool return_boolean = _dof_to_F.find(node) == _dof_to_F.end();
148 if (return_boolean)
149 return;
150
151 const auto normalized_F = normalizeQuantity(_dof_to_F, node);
152 const auto normalized_F_neighbor = normalizeQuantity(_dof_to_F_neighbor, node);
153
154 // This 'averaging' assumption below can probably be improved upon.
155 _dof_to_interface_F[node] = 0.5 * (normalized_F + normalized_F_neighbor);
156
157 for (const auto i : make_range(Moose::dim))
158 for (const auto j : make_range(Moose::dim))
159 {
160 if (!std::isfinite(MetaPhysicL::raw_value(normalized_F(i, j))))
161 throw MooseException("The deformation gradient on the secondary surface is not finite in "
162 "CohesiveZoneModelBase. MOOSE needs to cut the time step size.");
163 if (!std::isfinite(MetaPhysicL::raw_value(normalized_F_neighbor(i, j))))
164 throw MooseException("The deformation gradient on the primary surface is not finite in "
165 "CohesiveZoneModelBase. MOOSE needs to cut the time step size.");
166 }
167
168 const auto dof_to_interface_F_node = libmesh_map_find(_dof_to_interface_F, node);
169
170 const ADFactorizedRankTwoTensor C = dof_to_interface_F_node.transpose() * dof_to_interface_F_node;
171 const auto Uinv = MathUtils::sqrt(C).inverse().get();
172 _dof_to_interface_R[node] = dof_to_interface_F_node * Uinv;
173
174 // Transform interface jump according to two rotation matrices
175 const auto global_interface_displacement = _dof_to_interface_displacement_jump[node];
177 (_dof_to_interface_R[node] * _dof_to_rotation_matrix[node]).transpose() *
178 global_interface_displacement;
179}
for(PetscInt i=0;i< nvars;++i)
unsigned int dim
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_F
Map from degree of freedom to secondary, interpolated deformation gradient tensor.
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 *, ADRankTwoTensor > _dof_to_interface_R
Map from degree of freedom to interface rotation tensor.
std::unordered_map< const DofObject *, ADRankTwoTensor > _dof_to_interface_F
Map from degree of freedom to interface deformation gradient tensor.
FactorizedRankTwoTensorTempl< T > transpose() const
IntRange< T > make_range(T beg, T end)

Referenced by CohesiveZoneModelBase::reinit().

◆ computeFinalDisplacementJump()

void BilinearMixedModeCohesiveZoneModel::computeFinalDisplacementJump ( const Node *const  node)
protectedvirtual

Definition at line 228 of file BilinearMixedModeCohesiveZoneModel.C.

229{
230 using std::sqrt, std::pow;
231
232 const auto interface_displacement_jump =
234
235 const auto mixity_ratio = libmesh_map_find(_dof_to_mode_mixity_ratio, node);
236
237 const auto normalized_GI_c = normalizeQuantity(_dof_to_GI_c, node);
238 const auto normalized_GII_c = normalizeQuantity(_dof_to_GII_c, node);
239
240 _dof_to_delta_final[node] =
241 2.0 * normalized_GII_c / normalizeQuantity(_dof_to_shear_strength, node);
242
243 if (interface_displacement_jump(0) > _epsilon_tolerance)
244 {
246 {
247 _dof_to_delta_final[node] =
248 2.0 / _penalty_stiffness_czm / libmesh_map_find(_dof_to_delta_initial, node) *
249 (normalized_GI_c +
250 (normalized_GII_c - normalized_GI_c) *
251 pow(mixity_ratio * mixity_ratio / (1 + mixity_ratio * mixity_ratio),
253 }
255 {
256 const auto Gc_mixed =
257 pow(1.0 / normalized_GI_c, _power_law_parameter) +
258 pow(mixity_ratio * mixity_ratio / normalized_GII_c, _power_law_parameter);
259 _dof_to_delta_final[node] = (2.0 + 2.0 * mixity_ratio * mixity_ratio) /
261 libmesh_map_find(_dof_to_delta_initial, node) *
262 pow(Gc_mixed, -1.0 / _power_law_parameter);
263 }
264 }
265}
ExpressionBuilder::EBTerm pow(const ExpressionBuilder::EBTerm &left, T exponent)
std::unordered_map< const DofObject *, ADReal > _dof_to_GI_c
std::unordered_map< const DofObject *, ADReal > _dof_to_GII_c

Referenced by computeCZMTraction().

◆ computeGlobalTraction()

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

Compute global traction for mortar application.

Definition at line 460 of file CohesiveZoneModelBase.C.

461{
462 // First call does not have maps available
463 const auto it = _dof_to_czm_traction.find(node);
464 if (it == _dof_to_czm_traction.end())
465 return;
466
467 const auto local_traction_vector = it->second;
468 const auto rotation_matrix = libmesh_map_find(_dof_to_rotation_matrix, node);
469
470 _dof_to_czm_traction[node] = rotation_matrix * local_traction_vector;
471}

Referenced by CohesiveZoneModelBase::reinit().

◆ computeModeMixity()

void BilinearMixedModeCohesiveZoneModel::computeModeMixity ( const Node *const  node)
protectedvirtual

Definition at line 181 of file BilinearMixedModeCohesiveZoneModel.C.

182{
183 using std::sqrt;
184
185 const auto interface_displacement_jump =
187
188 if (interface_displacement_jump(0) > _epsilon_tolerance)
189 {
190 const auto delta_s =
191 sqrt(interface_displacement_jump(1) * interface_displacement_jump(1) +
192 interface_displacement_jump(2) * interface_displacement_jump(2) + _epsilon_tolerance);
193
194 _dof_to_mode_mixity_ratio[node] = delta_s / interface_displacement_jump(0);
195 }
196 else
198}

Referenced by computeCZMTraction().

◆ computeQpIProperties()

void BilinearMixedModeCohesiveZoneModel::computeQpIProperties ( )
overrideprotectedvirtual

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

Reimplemented from CohesiveZoneModelBase.

Definition at line 111 of file BilinearMixedModeCohesiveZoneModel.C.

112{
114
115 // Get the _dof_to_weighted_gap map
116 const auto * const dof = cast_ptr<const DofObject *>(_lower_secondary_elem->node_ptr(_i));
117
118 // TODO: Probably better to interpolate the deformation gradients.
121
122 _dof_to_GI_c[dof] += (*_test)[_i][_qp] * _GI_c_interpolation;
123 _dof_to_GII_c[dof] += (*_test)[_i][_qp] * _GII_c_interpolation;
124}
ADReal _GI_c_interpolation
Interpolated value of fracture paramter mode I.
ADReal _shear_strength_interpolation
Interpolated value of shear_strength.
ADReal _GII_c_interpolation
Interpolated value of fracture paramter mode II.
ADReal _normal_strength_interpolation
Interpolated value of normal_strength.
virtual void computeQpIProperties() override
Computes properties that are functions both of _qp and _i, for example the weighted gap.
Elem const *const & _lower_secondary_elem
unsigned int _i
Test function index.
unsigned int _qp
Quadrature point index for the mortar segments.

◆ computeQpProperties()

void BilinearMixedModeCohesiveZoneModel::computeQpProperties ( )
overrideprotectedvirtual

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

indepedent of shape functions

Reimplemented from CohesiveZoneModelBase.

Definition at line 100 of file BilinearMixedModeCohesiveZoneModel.C.

101{
103
108}
virtual void computeQpProperties() override
Computes properties that are functions only of the current quadrature point (_qp),...
const std::vector< Real > & _JxW_msm
const MooseArray< Real > & _coord
Member for handling change of coordinate systems (xyz, rz, spherical)

◆ constrainedByOwner()

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

Definition at line 46 of file BilinearMixedModeCohesiveZoneModel.h.

46{ 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.

224{
225 mooseAssert(nodal_index < lower_secondary_elem.n_nodes(),
226 "Nodal normal index must refer to a node on the secondary element.");
228 "AD contact normals should only be requested while recording nodal-normal "
229 "derivatives.");
230 const Node * const node = lower_secondary_elem.node_ptr(nodal_index);
231 const auto normal_it = _ad_nodal_normals.find(node);
232 if (normal_it != _ad_nodal_normals.end())
233 return normal_it->second;
234
235 amg().computeADNodalNormals([this](const Node & coordinate_node)
236 { return nodalCoordinate(coordinate_node); },
238 return libmesh_map_find(_ad_nodal_normals, node);
239}
void computeADNodalNormals(const std::function< ADPoint(const Node &)> &coordinate, std::unordered_map< const Node *, ADRealVectorValue > &nodal_normals) const
const AutomaticMortarGeneration & amg() const
bool shouldRecordNodalNormalDerivatives() const
Whether the current assembly records nodal-normal derivatives for this object's system.
ADPoint nodalCoordinate(const Node &node) const
Add displacement derivatives to the node's current coordinate.
std::unordered_map< const Node *, ADRealVectorValue > _ad_nodal_normals
AD nodal normals, each averaged over the secondary faces containing that 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 124 of file PenaltyWeightedGapUserObject.C.

125{
126 return _contact_pressure;
127}
ADVariableValue _contact_pressure
The contact pressure on the mortar segument quadrature points.

◆ contactTangentialPressureDirOne()

const ADVariableValue & CohesiveZoneModelBase::contactTangentialPressureDirOne ( ) const
overridevirtualinherited
Returns
The contact force at quadrature points on the mortar segment

Implements WeightedVelocitiesUserObject.

Definition at line 448 of file CohesiveZoneModelBase.C.

449{
451}
ADVariableValue _frictional_contact_traction_one
The first frictional contact pressure on the mortar segment quadrature points.

◆ contactTangentialPressureDirTwo()

const ADVariableValue & CohesiveZoneModelBase::contactTangentialPressureDirTwo ( ) const
overridevirtualinherited

Implements WeightedVelocitiesUserObject.

Definition at line 454 of file CohesiveZoneModelBase.C.

455{
457}
ADVariableValue _frictional_contact_traction_two
The second frictional contact pressure on the mortar segment quadrature points.

◆ 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.

69{
70 mooseAssert(nodal_index < lower_secondary_elem.n_nodes(),
71 "Nodal tangent index must refer to a node on the secondary element.");
73 "AD contact tangents should only be requested while recording nodal-normal "
74 "derivatives.");
75 const Node * const node = lower_secondary_elem.node_ptr(nodal_index);
76 auto it = _ad_nodal_tangents.find(node);
77 if (it == _ad_nodal_tangents.end())
79 .emplace(node,
81 contactNormal(lower_secondary_elem, nodal_index)))
82 .first;
83
84#ifndef NDEBUG
85 // A cached frame must belong to the current geometry. The tolerance matches the stored versus AD
86 // nodal normal check in AutomaticMortarGeneration::computeADNodalNormals.
87 const auto & stored_tangents = amg().getNodalTangents(lower_secondary_elem);
88 for (const auto direction : make_range(2))
89 mooseAssert(
90 (MetaPhysicL::raw_value(it->second[direction]) - stored_tangents[direction][nodal_index])
91 .norm() < 1e-10,
92 "The stored and AD secondary nodal tangents must use the same geometry state.");
93#endif
94
95 return it->second;
96}
std::array< MooseUtils::SemidynamicVector< Point, 9 >, 2 > getNodalTangents(const Elem &secondary_elem) 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.
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 quadratu...
std::array< Vector, 2 > householderTangents(const Vector &normal)
auto norm(const T &a)

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

◆ czmGlobalTraction()

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

Definition at line 474 of file CohesiveZoneModelBase.C.

475{
476 mooseAssert(i <= 3,
477 "Internal error in czmGlobalTraction. Index exceeds the traction vector size.");
478
480}
std::vector< ADVariableValue > _czm_interpolated_traction
The global traction.

Referenced by MortarGenericTraction::computeQpResidual().

◆ dispVar()

const MooseVariable * WeightedGapUserObject::dispVar ( unsigned int  component) const
inherited
Parameters
componentThe 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.

104{
105 switch (component)
106 {
107 case 0:
108 return _disp_x_var;
109 case 1:
110 return _disp_y_var;
111 case 2:
112 return _disp_z_var;
113 default:
114 mooseError("Invalid displacement component '", component, "' requested.");
115 }
116}
void mooseError(Args &&... args) const
const MooseVariable *const _disp_x_var
The x displacement variable.
const MooseVariable *const _disp_y_var
The y displacement variable.
const MooseVariable *const _disp_z_var
The z displacement variable.

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

◆ 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.

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

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

◆ execute()

void WeightedGapUserObject::execute ( )
overridevirtualinherited

Implements MortarUserObject.

Definition at line 293 of file WeightedGapUserObject.C.

294{
295 for (_qp = 0; _qp < _qrule_msm->n_points(); _qp++)
296 {
298 for (_i = 0; _i < _test->size(); ++_i)
300 }
301}
const libMesh::QBase *const & _qrule_msm
virtual void computeQpProperties()
Computes properties that are functions only of the current quadrature point (_qp),...
virtual void computeQpIProperties()
Computes properties that are functions both of _qp and _i, for example the weighted gap.
const VariableTestValue * _test
A pointer to the test function associated with the weighted gap.
unsigned int n_points() const

◆ finalize()

void BilinearMixedModeCohesiveZoneModel::finalize ( )
overridevirtual

Reimplemented from CohesiveZoneModelBase.

Definition at line 308 of file BilinearMixedModeCohesiveZoneModel.C.

309{
311
312 const bool send_data_back = !constrainedByOwner();
313
316
319
321 _dof_to_GI_c, _subproblem.mesh(), _nodal, _communicator, send_data_back);
322
324 _dof_to_GII_c, _subproblem.mesh(), _nodal, _communicator, send_data_back);
325}
virtual void finalize() override
virtual MooseMesh & mesh()=0
SubProblem & _subproblem
const bool _nodal
Whether the dof objects are nodal; if they're not, then they're elemental.
const Parallel::Communicator & _communicator
void communicateRealObject(std::unordered_map< const DofObject *, T > &dof_to_adreal, const MooseMesh &mesh, const bool nodal, const Parallel::Communicator &communicator, const bool send_data_back)

◆ 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.

192 {
193 const auto it = map.find(key);
194 if (it == map.end())
195 return default_value;
196 return it->second;
197 }

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

◆ getAccumulatedSlip()

virtual Real WeightedGapUserObject::getAccumulatedSlip ( const Node * const  ,
const unsigned int   
) const
inlinevirtualinherited
Parameters
nodeNode pointer
componentComponent of the local slip vector
Returns
The accumulated slip at the node

Reimplemented in PenaltyFrictionUserObject.

Definition at line 146 of file WeightedGapUserObject.h.

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

◆ getActiveSetState()

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

Definition at line 48 of file PenaltyWeightedGapUserObject.h.

49 {
50 const auto it = _dof_to_normal_pressure.find(_subproblem.mesh().nodePtr(node->id()));
51 return it != _dof_to_normal_pressure.end() && it->second > 0;
52 }
virtual const Node * nodePtr(const dof_id_type i) const
std::unordered_map< const DofObject *, ADReal > _dof_to_normal_pressure
Map from degree of freedom to normal pressure for reporting.

Referenced by MortarUserObjectAux::MortarUserObjectAux().

◆ getCohesiveDamage()

Real BilinearMixedModeCohesiveZoneModel::getCohesiveDamage ( const Node *const  node) const

Definition at line 339 of file BilinearMixedModeCohesiveZoneModel.C.

340{
341 const auto it = _dof_to_damage.find(node->id());
342
343 if (it != _dof_to_damage.end())
344 return MetaPhysicL::raw_value(it->second.first);
345 else
346 return 0.0;
347}

Referenced by CohesiveZoneMortarUserObjectAux::CohesiveZoneMortarUserObjectAux().

◆ getCohesiveTractionEffective()

Real BilinearMixedModeCohesiveZoneModel::getCohesiveTractionEffective ( const Node *const  node) const

Definition at line 460 of file BilinearMixedModeCohesiveZoneModel.C.

461{
462 const auto it = _dof_to_local_czm_traction.find(_subproblem.mesh().nodePtr(node->id()));
463
464 if (it != _dof_to_local_czm_traction.end())
465 return it->second.norm();
466 else
467 return 0.0;
468}

Referenced by CohesiveZoneMortarUserObjectAux::CohesiveZoneMortarUserObjectAux().

◆ getCohesiveTractionNormal()

Real BilinearMixedModeCohesiveZoneModel::getCohesiveTractionNormal ( const Node *const  node) const

Definition at line 413 of file BilinearMixedModeCohesiveZoneModel.C.

414{
415 const auto it = _dof_to_local_czm_traction.find(_subproblem.mesh().nodePtr(node->id()));
416
417 if (it != _dof_to_local_czm_traction.end())
418 return it->second(0);
419 else
420 return 0.0;
421}

Referenced by CohesiveZoneMortarUserObjectAux::CohesiveZoneMortarUserObjectAux().

◆ getCohesiveTractionTangentialMagnitude()

Real BilinearMixedModeCohesiveZoneModel::getCohesiveTractionTangentialMagnitude ( const Node *const  node) const

Definition at line 424 of file BilinearMixedModeCohesiveZoneModel.C.

426{
427 using std::sqrt;
428
429 const auto it = _dof_to_local_czm_traction.find(_subproblem.mesh().nodePtr(node->id()));
430
431 if (it != _dof_to_local_czm_traction.end())
432 return sqrt(it->second(1) * it->second(1) + it->second(2) * it->second(2));
433 else
434 return 0.0;
435}

Referenced by CohesiveZoneMortarUserObjectAux::CohesiveZoneMortarUserObjectAux().

◆ getCohesiveTractionTangentialOne()

Real BilinearMixedModeCohesiveZoneModel::getCohesiveTractionTangentialOne ( const Node *const  node) const

Definition at line 438 of file BilinearMixedModeCohesiveZoneModel.C.

439{
440 const auto it = _dof_to_local_czm_traction.find(_subproblem.mesh().nodePtr(node->id()));
441
442 if (it != _dof_to_local_czm_traction.end())
443 return it->second(1);
444 else
445 return 0.0;
446}

Referenced by CohesiveZoneMortarUserObjectAux::CohesiveZoneMortarUserObjectAux().

◆ getCohesiveTractionTangentialTwo()

Real BilinearMixedModeCohesiveZoneModel::getCohesiveTractionTangentialTwo ( const Node *const  node) const

Definition at line 449 of file BilinearMixedModeCohesiveZoneModel.C.

450{
451 const auto it = _dof_to_local_czm_traction.find(_subproblem.mesh().nodePtr(node->id()));
452
453 if (it != _dof_to_local_czm_traction.end())
454 return it->second(2);
455 else
456 return 0.0;
457}

Referenced by CohesiveZoneMortarUserObjectAux::CohesiveZoneMortarUserObjectAux().

◆ getDeltaTangentialLagrangeMultiplier()

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

Reimplemented in PenaltyFrictionUserObject.

Definition at line 44 of file PenaltyWeightedGapUserObject.h.

45 {
46 return 0.0;
47 };

◆ getFrictionalContactPressure()

virtual Real WeightedGapUserObject::getFrictionalContactPressure ( const Node * const  ,
const unsigned int   
) const
inlinevirtualinherited
Parameters
nodeNode pointer
componentComponent of the frictional pressure vector
Returns
The frictional contact pressure at the node

Reimplemented in PenaltyFrictionUserObject.

Definition at line 135 of file WeightedGapUserObject.h.

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

◆ getLocalDisplacementNormal()

Real BilinearMixedModeCohesiveZoneModel::getLocalDisplacementNormal ( const Node *const  node) const

Definition at line 350 of file BilinearMixedModeCohesiveZoneModel.C.

351{
352 const auto it = _dof_to_interface_displacement_jump.find(_subproblem.mesh().nodePtr(node->id()));
353 const auto it2 = _dof_to_weighted_gap.find(_subproblem.mesh().nodePtr(node->id()));
354
355 if (it != _dof_to_interface_displacement_jump.end() && it2 != _dof_to_weighted_gap.end())
356 return MetaPhysicL::raw_value(it->second(0) / it2->second.second);
357 else
358 return 0.0;
359}

Referenced by CohesiveZoneMortarUserObjectAux::CohesiveZoneMortarUserObjectAux().

◆ getLocalDisplacementTangential()

Real BilinearMixedModeCohesiveZoneModel::getLocalDisplacementTangential ( const Node *const  node) const

◆ getLocalDisplacementTangentialEffective()

Real BilinearMixedModeCohesiveZoneModel::getLocalDisplacementTangentialEffective ( const Node *const  node) const

Definition at line 392 of file BilinearMixedModeCohesiveZoneModel.C.

394{
395 using std::sqrt;
396
397 const auto it = _dof_to_interface_displacement_jump.find(_subproblem.mesh().nodePtr(node->id()));
398 const auto it2 = _dof_to_weighted_gap.find(_subproblem.mesh().nodePtr(node->id()));
399
400 if (it != _dof_to_interface_displacement_jump.end() && it2 != _dof_to_weighted_gap.end())
401 {
402 // Take raw values before the sqrt: the AD derivative of sqrt divides by zero
403 // when the tangential jump is exactly zero (e.g. pure mode I or post-debond).
404 const auto tangential_one = MetaPhysicL::raw_value(it->second(1) / it2->second.second);
405 const auto tangential_two = MetaPhysicL::raw_value(it->second(2) / it2->second.second);
406 return sqrt(tangential_one * tangential_one + tangential_two * tangential_two);
407 }
408 else
409 return 0.0;
410}

Referenced by CohesiveZoneMortarUserObjectAux::CohesiveZoneMortarUserObjectAux().

◆ getLocalDisplacementTangentialOne()

Real BilinearMixedModeCohesiveZoneModel::getLocalDisplacementTangentialOne ( const Node *const  node) const

Definition at line 368 of file BilinearMixedModeCohesiveZoneModel.C.

369{
370 const auto it = _dof_to_interface_displacement_jump.find(_subproblem.mesh().nodePtr(node->id()));
371 const auto it2 = _dof_to_weighted_gap.find(_subproblem.mesh().nodePtr(node->id()));
372
373 if (it != _dof_to_interface_displacement_jump.end() && it2 != _dof_to_weighted_gap.end())
374 return MetaPhysicL::raw_value(it->second(1) / it2->second.second);
375 else
376 return 0.0;
377}

Referenced by CohesiveZoneMortarUserObjectAux::CohesiveZoneMortarUserObjectAux(), and getLocalDisplacementTangential().

◆ getLocalDisplacementTangentialTwo()

Real BilinearMixedModeCohesiveZoneModel::getLocalDisplacementTangentialTwo ( const Node *const  node) const

Definition at line 380 of file BilinearMixedModeCohesiveZoneModel.C.

381{
382 const auto it = _dof_to_interface_displacement_jump.find(_subproblem.mesh().nodePtr(node->id()));
383 const auto it2 = _dof_to_weighted_gap.find(_subproblem.mesh().nodePtr(node->id()));
384
385 if (it != _dof_to_interface_displacement_jump.end() && it2 != _dof_to_weighted_gap.end())
386 return MetaPhysicL::raw_value(it->second(2) / it2->second.second);
387 else
388 return 0.0;
389}

Referenced by CohesiveZoneMortarUserObjectAux::CohesiveZoneMortarUserObjectAux().

◆ getModeMixityRatio()

Real BilinearMixedModeCohesiveZoneModel::getModeMixityRatio ( const Node *const  node) const

Definition at line 328 of file BilinearMixedModeCohesiveZoneModel.C.

329{
330 const auto it = _dof_to_mode_mixity_ratio.find(_subproblem.mesh().nodePtr(node->id()));
331
332 if (it != _dof_to_mode_mixity_ratio.end())
333 return MetaPhysicL::raw_value(it->second);
334 else
335 return 0.0;
336}

Referenced by CohesiveZoneMortarUserObjectAux::CohesiveZoneMortarUserObjectAux().

◆ getNormalContactPressure()

Real PenaltyWeightedGapUserObject::getNormalContactPressure ( const Node *const  ) const
overridevirtualinherited
Parameters
nodeNode pointer
Returns
The normal contact pressure at the node

Implements WeightedGapUserObject.

Definition at line 175 of file PenaltyWeightedGapUserObject.C.

176{
177 const auto it = _dof_to_normal_pressure.find(_subproblem.mesh().nodePtr(node->id()));
178
179 if (it != _dof_to_normal_pressure.end())
180 return MetaPhysicL::raw_value(it->second);
181 else
182 return 0.0;
183}

◆ getNormalGap()

Real WeightedGapUserObject::getNormalGap ( const Node * const  node) const
virtualinherited
Parameters
nodeNode pointer
Returns
The normal gap at the node

Definition at line 304 of file WeightedGapUserObject.C.

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

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

◆ getNormalLagrangeMultiplier()

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

Definition at line 186 of file PenaltyWeightedGapUserObject.C.

187{
188 const auto it = _dof_to_lagrange_multiplier.find(_subproblem.mesh().nodePtr(node->id()));
189
190 if (it != _dof_to_lagrange_multiplier.end())
191 return -MetaPhysicL::raw_value(it->second);
192 else
193 return 0.0;
194}
std::unordered_map< const DofObject *, Real > _dof_to_lagrange_multiplier
Map from degree of freedom to augmented lagrange multiplier.

Referenced by MortarUserObjectAux::MortarUserObjectAux().

◆ getTangentialVelocity()

virtual Real WeightedGapUserObject::getTangentialVelocity ( const Node * const  ,
const unsigned int   
) const
inlinevirtualinherited
Parameters
nodeNode pointer
componentComponent 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.

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

Referenced by MortarUserObjectAux::MortarUserObjectAux().

◆ initialize()

void BilinearMixedModeCohesiveZoneModel::initialize ( )
overridevirtual

Reimplemented from CohesiveZoneModelBase.

Definition at line 127 of file BilinearMixedModeCohesiveZoneModel.C.

128{
129 // instead we call it explicitly here
131
132 // Avoid accumulating interpolation over the time step
137 _dof_to_GI_c.clear();
138 _dof_to_GII_c.clear();
139 _dof_to_delta_initial.clear();
140 _dof_to_delta_final.clear();
141 _dof_to_delta_max.clear();
142}
virtual void initialize() override

◆ initialSetup()

void WeightedGapUserObject::initialSetup ( )
overridevirtualinherited

Reimplemented from MortarUserObject.

Definition at line 119 of file WeightedGapUserObject.C.

120{
122 _test = &test();
123
125 return;
126
127 if (getParam<bool>("interpolate_normals"))
128 paramError("interpolate_normals",
129 "Nodal-normal derivatives require normalized secondary nodal normals and cannot be "
130 "combined with quadrature-point normal interpolation.");
131
132 const std::array<std::pair<const MooseVariable *, const char *>, 3> displacement_variables{
133 {{_disp_x_var, "disp_x"}, {_disp_y_var, "disp_y"}, {_disp_z_var, "disp_z"}}};
134 for (const auto & [variable, parameter_name] : displacement_variables)
135 if (variable)
136 {
137 if (!variable->isNodal())
138 paramError(parameter_name,
139 "Nodal-normal derivatives require a nodal displacement variable.");
140 if (&variable->sys() != &_sys)
141 paramError(parameter_name,
142 "Nodal-normal derivatives require displacement variables in the nonlinear "
143 "system assembled by this contact object.");
144 }
145}
virtual void initialSetup()
SystemBase & _sys
virtual bool usesNodalNormalDerivatives() const =0
Whether this formulation includes displacement derivatives of the secondary nodal normals.
virtual const VariableTestValue & test() const =0
if(subdm)

◆ isAugmentedLagrangianConverged()

bool PenaltyWeightedGapUserObject::isAugmentedLagrangianConverged ( )
overridevirtualinherited

Implements AugmentedLagrangeInterface.

Reimplemented in PenaltyFrictionUserObject.

Definition at line 255 of file PenaltyWeightedGapUserObject.C.

256{
257 using std::max, std::min, std::abs;
258
259 Real max_positive_gap = 0.0;
260 Real min_negative_gap = 0.0;
261
262 // Get maximum gap to ascertain whether we are converged.
263 for (auto & [dof_object, wgap] : _dof_to_weighted_gap)
264 {
265 const auto gap = physicalGap(wgap);
266 {
267 // Check condition for nodes that are active
268 if (gap < 0 || _dof_to_lagrange_multiplier[dof_object] < 0.0)
269 {
270 max_positive_gap = max(max_positive_gap, gap);
271 min_negative_gap = min(min_negative_gap, gap);
272 }
273 }
274 }
275
276 // Communicate the extreme gap values in parallel.
277 std::vector<Real> recv;
278 if (this->_communicator.rank() == 0)
279 recv.resize(this->_communicator.size());
280 this->_communicator.gather(
281 0, -min_negative_gap > max_positive_gap ? min_negative_gap : max_positive_gap, recv);
282
283 if (this->_communicator.rank() == 0)
284 {
285 min_negative_gap = *std::min_element(recv.begin(), recv.end());
286 max_positive_gap = *std::max_element(recv.begin(), recv.end());
287
288 // report the gap value with the largest magnitude
289 const Real reported_gap =
290 -min_negative_gap > max_positive_gap ? min_negative_gap : max_positive_gap;
291 if (abs(reported_gap) > _penetration_tolerance)
292 {
293 mooseInfoRepeated("Penetration tolerance fail max_gap = ",
294 reported_gap,
295 " (gap_tol=",
297 "). Iteration number is: ",
299 ".");
300 return false;
301 }
302 else
303 mooseInfoRepeated("Penetration tolerance success max_gap = ",
304 reported_gap,
305 " (gap_tol=",
307 "). Iteration number is: ",
309 ".");
310 }
311
312 return true;
313}
void mooseInfoRepeated(Args &&... args)
processor_id_type size() const
processor_id_type rank() const
void gather(const unsigned int root_id, const T &send_data, std::vector< T, A > &recv) const

Referenced by PenaltyFrictionUserObject::isAugmentedLagrangianConverged().

◆ nodalContactPressure()

ADReal PenaltyWeightedGapUserObject::nodalContactPressure ( const Node &  node) const
overridevirtualinherited
Parameters
nodeSecondary 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 213 of file PenaltyWeightedGapUserObject.C.

214{
215 // A node only enters this map once it has received a weighted-gap contribution, so a node on a
216 // partially integrated element may be absent. Absent means no contact pressure, matching how the
217 // interpolated pressure in reinit() treats it.
218 return findValue(_dof_to_normal_pressure, cast_ptr<const DofObject *>(&node), ADReal(0));
219}
DualNumber< Real, DNDerivativeType, true > ADReal
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.

◆ 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.

243{
245 "AD nodal coordinates should only be requested while recording nodal-normal "
246 "derivatives.");
248 "AD nodal coordinates require derivative recording.");
249 ADPoint point = node;
250
251 const std::array<std::pair<const MooseVariable *, unsigned int>, 3> displacement_variables{
252 {{_disp_x_var, 0}, {_disp_y_var, 1}, {_disp_z_var, 2}}};
253 for (const auto & [variable, component] : displacement_variables)
254 if (variable)
255 {
256 const auto sys_num = variable->sys().number();
257 const auto var_num = variable->number();
258 // Enriched elements can have mapping nodes without independent displacement degrees of
259 // freedom. Their positions enter the geometric map, but they contribute no independent
260 // displacement derivative.
261 if (node.n_dofs(sys_num, var_num))
263 point(component).derivatives(), node.dof_number(sys_num, var_num, 0), 1.0);
264 }
265
266 return point;
267}
bool doDerivatives(const SubProblem &subproblem, const SystemBase &sys)
void derivInsert(SemiDynamicSparseNumberArray< Real, libMesh::dof_id_type, NWrapper< N > > &derivs, libMesh::dof_id_type index, Real value)

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
overridevirtualinherited
Parameters
nodeSecondary node of the current secondary lower-dimensional element
directionTangent 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.

252{
253 mooseAssert(direction < 2, "There are at most two tangent directions.");
254 return libmesh_map_find(_dof_to_tangential_traction, node.id()).first(direction);
255}
std::unordered_map< dof_id_type, std::pair< ADTwoVector, TwoVector > > & _dof_to_tangential_traction

◆ normalizeCDivisor()

virtual Real WeightedGapUserObject::normalizeCDivisor ( const DofObject * const  ,
const Real  covered_normalization 
) const
inlinevirtualinherited
Parameters
dofNode/element dof object
covered_normalizationThe 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.

124 {
125 return covered_normalization;
126 }

Referenced by ComputeWeightedGapLMMechanicalContact::enforceConstraintOnDof().

◆ normalizeQuantity()

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

Normalize mortar quantities (remove mortar integral scaling)

Definition at line 60 of file CohesiveZoneModelBase.h.

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

Referenced by computeCriticalDisplacementJump(), computeCZMTraction(), computeEffectiveDisplacementJump(), CohesiveZoneModelBase::computeFandR(), computeFinalDisplacementJump(), and 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.

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

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

◆ prepareJumpKinematicQuantities()

void CohesiveZoneModelBase::prepareJumpKinematicQuantities ( )
protectedvirtualinherited

Definition at line 405 of file CohesiveZoneModelBase.C.

406{
407 // Compute rotation matrix from secondary surface
408 // Rotation matrices and local interface displacement jump.
409 for (const auto i : make_range(_test->size()))
410 {
411 const Node * const node = _lower_secondary_elem->node_ptr(i);
412
413 // First call does not have maps available
414 const bool skip_node = _dof_to_weighted_gap.find(node) == _dof_to_weighted_gap.end();
415 if (skip_node)
416 continue;
417
418 _dof_to_rotation_matrix[node] = CohesiveZoneModelTools::computeReferenceRotation<true>(
420
421 // Every time we used quantities from a map we "denormalize" it from the mortar integral.
422 // See normalizing member functions.
423 if (auto it = _dof_to_weighted_displacements.find(node);
425 _dof_to_interface_displacement_jump[node] = it->second;
426 }
427}
virtual unsigned int dimension() const
std::vector< Point > _normals
std::unordered_map< const DofObject *, ADRealVectorValue > _dof_to_weighted_displacements
A map from node to weighted displacements.

Referenced by CohesiveZoneModelBase::reinit().

◆ reinit()

void CohesiveZoneModelBase::reinit ( )
overridevirtualinherited

Reimplemented from PenaltyWeightedGapUserObject.

Definition at line 258 of file CohesiveZoneModelBase.C.

259{
260 // Normal contact pressure with penalty
262
263 // Compute all rotations that were created as material properties in CZMComputeDisplacementJump
265
266 // Reset frictional pressure
269
270 for (const auto qp : make_range(_qrule_msm->n_points()))
271 {
274 }
275
276 // zero vector
277 const static TwoVector zero{0.0, 0.0};
278
279 // iterate over nodes
280 for (const auto i : make_range(_test->size()))
281 {
282 // current node
283 const Node * const node = _lower_secondary_elem->node_ptr(i);
284
285 // Compute the weighted nodal deformation gradient and rotation tensors.
286 computeFandR(node);
287
288 // The call below is a 'macro' call. Create a utility function or user object for it.
289 computeCZMTraction(node);
290
291 // Build final traction vector
293
294 // Compute mechanical contact until end of method.
295 const auto penalty_friction = findValue(
296 _dof_to_local_penalty_friction, cast_ptr<const DofObject *>(node), _penalty_friction);
297
298 // utilized quantities
299 const auto & normal_pressure = _dof_to_normal_pressure[node];
300
301 // map the tangential traction and accumulated slip
302 auto & [tangential_traction, old_tangential_traction] = _dof_to_tangential_traction[node->id()];
303 auto & [accumulated_slip, old_accumulated_slip] = _dof_to_accumulated_slip[node->id()];
304
305 Real normal_lm = -1;
306 if (auto it = _dof_to_lagrange_multiplier.find(node); it != _dof_to_lagrange_multiplier.end())
307 normal_lm = it->second;
308
309 // Keep active set fixed from second Uzawa loop
310 if (normal_lm < -TOLERANCE && normal_pressure > TOLERANCE)
311 {
312 using std::abs;
313
314 const auto & damage = _dof_to_damage[node->id()].first;
315
316 const auto & real_tangential_velocity =
317 libmesh_map_find(_dof_to_real_tangential_velocity, node);
318 const ADTwoVector slip_distance = {real_tangential_velocity[0] * _dt,
319 real_tangential_velocity[1] * _dt};
320
321 // frictional lagrange multiplier (delta lambda^(k)_T)
322 const auto & tangential_lm =
324
325 // tangential trial traction (Simo 3.12)
326 // Modified for implementation in MOOSE: Avoid pingponging on frictional sign (max. 0.4
327 // capacity)
328 ADTwoVector inner_iteration_penalty_friction = penalty_friction * slip_distance;
329
330 const auto slip_metric = std::abs(MetaPhysicL::raw_value(slip_distance).cwiseAbs()(0)) +
331 std::abs(MetaPhysicL::raw_value(slip_distance).cwiseAbs()(1));
332
333 const auto slip_norm = (MetaPhysicL::raw_value(slip_distance)).norm();
334 const auto friction_limit = 0.4 * _friction_coefficient * damage * abs(normal_pressure);
335
336 if (slip_metric > _epsilon_tolerance && penalty_friction * slip_norm > friction_limit)
337 {
338 inner_iteration_penalty_friction =
339 MetaPhysicL::raw_value(friction_limit / (penalty_friction * slip_norm)) *
340 penalty_friction * slip_distance;
341 }
342
343 ADTwoVector tangential_trial_traction =
344 old_tangential_traction + tangential_lm + inner_iteration_penalty_friction;
345
346 // Nonlinearity below
347 ADReal tangential_trial_traction_norm =
348 (MetaPhysicL::raw_value(tangential_trial_traction)).norm();
349 ADReal phi_trial = tangential_trial_traction_norm - _friction_coefficient * normal_pressure;
350 tangential_traction = tangential_trial_traction;
351
352 // Simo considers this a 'return mapping'; we are just capping friction to the Coulomb limit.
353 if (phi_trial > 0.0)
354 // Simo 3.13/3.14 (the penalty formulation has an error in the paper)
355 tangential_traction -=
356 phi_trial * tangential_trial_traction / tangential_trial_traction_norm;
357
358 // track accumulated slip for output purposes
359 accumulated_slip = old_accumulated_slip + slip_distance.cwiseAbs();
360 }
361 else
362 {
363 // reset slip and clear traction
364 accumulated_slip.setZero();
365 tangential_traction.setZero();
366 }
367
368 // Now that we have consistent nodal frictional values, create an interpolated frictional
369 // pressure variable.
370 const auto test_i = (*_test)[i];
371 for (const auto qp : make_range(_qrule_msm->n_points()))
372 {
373 _frictional_contact_traction_one[qp] += test_i[qp] * tangential_traction(0);
374 _frictional_contact_traction_two[qp] += test_i[qp] * tangential_traction(1);
375 }
376 }
377
378 for (const auto i : index_range(_czm_interpolated_traction))
379 {
381
382 for (const auto qp : make_range(_qrule_msm->n_points()))
383 _czm_interpolated_traction[i][qp] = 0.0;
384 }
385
386 // iterate over nodes
387 for (const auto i : make_range(_test->size()))
388 {
389 // current node
390 const Node * const node = _lower_secondary_elem->node_ptr(i);
391
392 // End of CZM bilinear computations
393 auto it = _dof_to_czm_traction.find(node);
394 if (it == _dof_to_czm_traction.end())
395 continue;
396
397 const auto & test_i = (*_test)[i];
398 for (const auto qp : make_range(_qrule_msm->n_points()))
400 _czm_interpolated_traction[idx][qp] += test_i[qp] * (it->second)(idx);
401 }
402}
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.
const Real _penalty_friction
The penalty factor for the frictional constraints.
virtual void computeGlobalTraction(const Node *const node)
Compute global traction for mortar application.
std::unordered_map< const DofObject *, Real > _dof_to_local_penalty_friction
Map from degree of freedom to local friction penalty value.
const Real _friction_coefficient
The friction coefficient.
std::unordered_map< const DofObject *, TwoVector > _dof_to_frictional_lagrange_multipliers
Map from degree of freedom to augmented lagrange multiplier.
virtual void computeCZMTraction(const Node *const)=0
Encapsulate the CZM constitutive behavior.
virtual void prepareJumpKinematicQuantities()
virtual void computeFandR(const Node *const node)
AugmentedLagrangianContactProblemInterface *const _augmented_lagrange_problem
augmented Lagrange problem and iteration number
std::unordered_map< const DofObject *, std::array< ADReal, 2 > > _dof_to_real_tangential_velocity
A map from node to two interpolated, physical tangential velocities.
unsigned int idx(const ElemType type, const unsigned int nx, const unsigned int i, const unsigned int j)
auto index_range(const T &sizable)
const Number zero

◆ selfFinalize()

void PenaltyWeightedGapUserObject::selfFinalize ( )
protectedinherited

Definition at line 144 of file PenaltyWeightedGapUserObject.C.

145{
146 using std::min;
147
148 // compute new normal pressure for each node
149 for (const auto & [dof_object, wgap] : _dof_to_weighted_gap)
150 {
151 auto penalty = findValue(_dof_to_local_penalty, dof_object, _penalty);
152
153 // If _use_physical_gap is true we "normalize" the penalty parameter with the surrounding area.
154 auto gap = _use_physical_gap ? adPhysicalGap(wgap) / wgap.second : wgap.first;
155
156 const auto lagrange_multiplier =
158
159 // keep the negative normal pressure (compressive per convention here)
160 auto normal_pressure = min(0.0, penalty * gap + lagrange_multiplier);
161
162 // we switch conventins here and consider positive normal pressure as compressive
163 _dof_to_normal_pressure[dof_object] = -normal_pressure;
164 }
165}
std::unordered_map< const DofObject *, Real > _dof_to_local_penalty
Map from degree of freedom to local penalty value.
const bool _use_physical_gap
Use scaled or physical gap.
ADReal adPhysicalGap(const std::pair< ADReal, Real > &gap) const

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

◆ selfInitialize() [1/2]

void PenaltyWeightedGapUserObject::selfInitialize ( )
protectedinherited

Definition at line 130 of file PenaltyWeightedGapUserObject.C.

131{
132 for (auto & dof_pn : _dof_to_normal_pressure)
133 dof_pn.second = 0.0;
134}

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

◆ selfInitialize() [2/2]

void WeightedVelocitiesUserObject::selfInitialize ( )
protectedinherited

◆ selfTimestepSetup() [1/2]

void PenaltyWeightedGapUserObject::selfTimestepSetup ( )
protectedinherited

Definition at line 222 of file PenaltyWeightedGapUserObject.C.

223{
224 // Let's not clear the LMs for improved performance in
225 // nonlinear problems
226
227 // timestepSetup() runs once per time step attempt, so on a retried timestep (e.g. --test-restep
228 // or a rejected step) the LMs hold the discarded attempt's values; restore the accepted ones
229 const bool retried_timestep = _fe_problem.timeStep() == _t_step_old_normal;
231 if (retried_timestep)
233 else
235
236 // reset penalty
237 for (auto & dof_lp : _dof_to_local_penalty)
238 dof_lp.second = _penalty;
239
240 // clear previous gap
241 for (auto & dof_pg : _dof_to_previous_gap)
242 dof_pg.second = 0.0;
243
244 // save old timestep
245 _dt_old = _dt;
246}
virtual int & timeStep() const
std::unordered_map< const DofObject *, Real > _dof_to_lagrange_multiplier_start
Augmented lagrange multipliers at the start of the current timestep, restored when the timestep is re...
int & _t_step_old_normal
The timestep index at which _dof_to_lagrange_multiplier_start was last saved.
FEProblemBase & _fe_problem
The base finite element problem.

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

◆ selfTimestepSetup() [2/2]

void WeightedVelocitiesUserObject::selfTimestepSetup ( )
inlineprotectedinherited

◆ shouldRecordNodalNormalDerivatives()

bool WeightedGapUserObject::shouldRecordNodalNormalDerivatives ( ) const
inherited

◆ tangentialTractionBasis()

virtual const VariableTestValue & WeightedVelocitiesUserObject::tangentialTractionBasis ( unsigned int  ) const
inlinevirtualinherited
Parameters
directionTangent 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.

70 {
71 return tractionBasis();
72 }
virtual const VariableTestValue & tractionBasis() const

Referenced by TangentialMortarMechanicalContact::precomputeQpQuantities().

◆ test()

const VariableTestValue & CohesiveZoneModelBase::test ( ) const
overrideprotectedvirtualinherited
Returns
The test function associated with the weighted gap

Reimplemented from PenaltyWeightedGapUserObject.

Definition at line 105 of file CohesiveZoneModelBase.C.

106{
108}
virtual const FieldVariablePhiValue & phiLower() const override
const MooseVariable *const _aux_lm_var
The auxiliary Lagrange multiplier variable (used together whith the Petrov-Galerkin approach)

◆ timestepSetup()

void CohesiveZoneModelBase::timestepSetup ( )
overridevirtualinherited

Reimplemented from PenaltyWeightedGapUserObject.

Definition at line 182 of file CohesiveZoneModelBase.C.

183{
184 // these functions do not call WeightedGapUserObject::timestepSetup to avoid double initialization
187
188 // instead we call it explicitly here
190
191 // Clear step slip (values used in between AL iterations for penalty adaptivity)
192 for (auto & map_pr : _dof_to_step_slip)
193 {
194 auto & [step_slip, old_step_slip] = map_pr.second;
195 old_step_slip = {0.0, 0.0};
196 step_slip = {0.0, 0.0};
197 }
198
199 // save off accumulated slip from the last timestep
200 for (auto & map_pr : _dof_to_accumulated_slip)
201 {
202 auto & [accumulated_slip, old_accumulated_slip] = map_pr.second;
203 old_accumulated_slip = MetaPhysicL::raw_value(accumulated_slip);
204 }
205
206 for (auto & dof_lp : _dof_to_local_penalty_friction)
207 dof_lp.second = _penalty_friction;
208
209 // save off tangential traction from the last timestep
210 for (auto & map_pr : _dof_to_tangential_traction)
211 {
212 auto & [tangential_traction, old_tangential_traction] = map_pr.second;
213 old_tangential_traction = {MetaPhysicL::raw_value(tangential_traction(0)),
214 MetaPhysicL::raw_value(tangential_traction(1))};
215 tangential_traction = {0.0, 0.0};
216 }
217
218 for (auto & [dof_object, delta_tangential_lm] : _dof_to_frictional_lagrange_multipliers)
219 delta_tangential_lm.setZero();
220
221 // save off damage from the last timestep
222 for (auto & map_pr : _dof_to_damage)
223 {
224 auto & [damage, old_damage] = map_pr.second;
225 old_damage = {MetaPhysicL::raw_value(damage)};
226 damage = {0.0};
227 }
228}
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.
virtual void timestepSetup()

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

118{
119 // Penalty contact interpolates traction and weighted gap with the same basis.
120 return test();
121}
virtual const VariableTestValue & test() const override

◆ updateAugmentedLagrangianMultipliers()

void PenaltyWeightedGapUserObject::updateAugmentedLagrangianMultipliers ( )
overridevirtualinherited

Implements AugmentedLagrangeInterface.

Reimplemented in PenaltyFrictionUserObject.

Definition at line 328 of file PenaltyWeightedGapUserObject.C.

329{
330 using std::abs;
331
332 for (const auto & [dof_object, wgap] : _dof_to_weighted_gap)
333 {
334 auto & penalty = _dof_to_local_penalty[dof_object];
335 if (penalty == 0.0)
336 penalty = _penalty;
337
338 const auto gap = getNormalGap(static_cast<const Node *>(dof_object));
339 auto & lagrange_multiplier = _dof_to_lagrange_multiplier[dof_object];
340
341 const auto possible_normalization =
342 (_use_physical_gap ? libmesh_map_find(_dof_to_weighted_gap, dof_object).second : 1.0);
343
344 // Update penalty (the factor of 1/4 is suggested in the literature, the limit on AL iteration
345 // caps the penalty increase)
346 // Before we were updating the LM before adapting the penalty factor
347 if (lagrange_multiplier + gap * penalty / possible_normalization <= 0)
348 lagrange_multiplier += gap * penalty / possible_normalization;
349 else
350 lagrange_multiplier = 0.0;
351
352 const auto previous_gap = _dof_to_previous_gap[dof_object];
353 Real eval_tn = 0;
354
356 {
357 if (abs(gap) > 0.25 * abs(previous_gap) && abs(gap) > _penetration_tolerance)
358 penalty *= _penalty_multiplier;
359 }
361 bussettaAdaptivePenalty(previous_gap, gap, penalty, eval_tn);
362
363 // Heuristics to bound the penalty factor
364 if (penalty < _penalty)
365 penalty = _penalty;
366 else if (penalty > _max_penalty_multiplier * _penalty)
368 }
369}
const Real _penalty_multiplier
penalty growth factor for augmented Lagrange
const Real _max_penalty_multiplier
Maximum multiplier applied to the initial penalty factor in AL.
void bussettaAdaptivePenalty(const Real previous_gap, const Real gap, Real &penalty, Real &eval_tn)
Adaptive, local penalty for AL.
virtual Real getNormalGap(const Node *const) const

Referenced by PenaltyFrictionUserObject::updateAugmentedLagrangianMultipliers().

◆ usesNodalNormalDerivatives()

virtual bool CohesiveZoneModelBase::usesNodalNormalDerivatives ( ) const
inlineoverridevirtualinherited

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 BilinearMixedModeCohesiveZoneModel::validParams ( )
static

Definition at line 28 of file BilinearMixedModeCohesiveZoneModel.C.

29{
31
32 params.addClassDescription("Computes the bilinear mixed mode cohesive zone model.");
33
34 // Input parameters for bilinear mixed mode traction.
35 params.addParam<MaterialPropertyName>("GI_c",
36 "Critical energy release rate in normal direction.");
37 params.addParam<MaterialPropertyName>("GII_c",
38 "Critical energy release rate in shear direction.");
39 params.addParam<MaterialPropertyName>("normal_strength", "Tensile strength in normal direction.");
40 params.addParam<MaterialPropertyName>("shear_strength", "Tensile strength in shear direction.");
41 params.addParam<Real>("power_law_parameter", "The power law parameter.");
42 MooseEnum criterion("POWER_LAW BK", "BK");
43 params.addParam<Real>("viscosity", 0.0, "Viscosity for damage model.");
44 params.addParam<MooseEnum>(
45 "mixed_mode_criterion", criterion, "Option for mixed mode propagation criterion.");
46 params.addDeprecatedParam<bool>(
47 "lag_displacement_jump",
48 false,
49 "Whether to use old displacement jumps to compute the effective displacement jump.",
50 "This parameter has no effect in the mortar cohesive zone model and will be removed. Use "
51 "'viscosity' to regularize the damage evolution.");
52 params.addParam<bool>("set_compressive_traction_to_zero",
53 false,
54 "Zero compressive traction (set to true, allowing the use of standard "
55 "zero-penetration mortar contact constraints in "
56 "the normal direction).");
57 params.addParam<Real>(
58 "regularization_alpha", 1e-10, "Regularization parameter for the Macaulay bracket.");
60 "penalty_stiffness", "penalty_stiffness > 0.0", "Penalty stiffness for CZM.");
62 "GI_c GII_c normal_strength shear_strength power_law_parameter viscosity "
63 "mixed_mode_criterion regularization_alpha "
64 "penalty_stiffness",
65 "Bilinear mixed mode traction");
66 // End of input parameters for bilinear mixed mode traction.
67
68 return params;
69}
static InputParameters validParams()
void addParamNamesToGroup(const std::string &space_delim_names, const std::string group_name)
void addDeprecatedParam(const std::string &name, const T &value, const std::string &doc_string, const std::string &deprecation_message)
void addParam(const std::string &name, const std::initializer_list< typename T::value_type > &value, const std::string &doc_string)
void addClassDescription(const std::string &doc_string)
void addRangeCheckedParam(const std::string &name, const T &value, const std::string &parsed_function, const std::string &doc_string)

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

◆ _augmented_lagrange_problem

AugmentedLagrangianContactProblemInterface* const PenaltyWeightedGapUserObject::_augmented_lagrange_problem
protectedinherited

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

Referenced by CohesiveZoneModelBase::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(), computeQpProperties(), CohesiveZoneModelBase::computeQpProperties(), and WeightedVelocitiesUserObject::computeQpProperties().

◆ _czm_interpolated_traction

std::vector<ADVariableValue> CohesiveZoneModelBase::_czm_interpolated_traction
protectedinherited

The global traction.

Definition at line 100 of file CohesiveZoneModelBase.h.

Referenced by CohesiveZoneModelBase::czmGlobalTraction(), and CohesiveZoneModelBase::reinit().

◆ _disp_x_var

const MooseVariable* const WeightedGapUserObject::_disp_x_var
protectedinherited

◆ _disp_y_var

const MooseVariable* const WeightedGapUserObject::_disp_y_var
protectedinherited

◆ _disp_z_var

const MooseVariable* const WeightedGapUserObject::_disp_z_var
protectedinherited

◆ _dof_to_accumulated_slip

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

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

Definition at line 112 of file CohesiveZoneModelBase.h.

Referenced by CohesiveZoneModelBase::reinit(), and CohesiveZoneModelBase::timestepSetup().

◆ _dof_to_czm_traction

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

Total Lagrangian stress to be applied on CZM interface.

Definition at line 133 of file CohesiveZoneModelBase.h.

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

◆ _dof_to_damage

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

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

Definition at line 139 of file CohesiveZoneModelBase.h.

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

◆ _dof_to_delta_final

std::unordered_map<const DofObject *, ADReal> BilinearMixedModeCohesiveZoneModel::_dof_to_delta_final
protected

◆ _dof_to_delta_initial

std::unordered_map<const DofObject *, ADReal> BilinearMixedModeCohesiveZoneModel::_dof_to_delta_initial
protected

◆ _dof_to_delta_max

std::unordered_map<const DofObject *, ADReal> BilinearMixedModeCohesiveZoneModel::_dof_to_delta_max
protected

◆ _dof_to_displacement_jump

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

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
protectedinherited

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

Definition at line 88 of file CohesiveZoneModelBase.h.

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

◆ _dof_to_F_neighbor

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

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

Definition at line 91 of file CohesiveZoneModelBase.h.

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

◆ _dof_to_frictional_lagrange_multipliers

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

Map from degree of freedom to augmented lagrange multiplier.

Definition at line 124 of file CohesiveZoneModelBase.h.

Referenced by CohesiveZoneModelBase::reinit(), and CohesiveZoneModelBase::timestepSetup().

◆ _dof_to_GI_c

std::unordered_map<const DofObject *, ADReal> BilinearMixedModeCohesiveZoneModel::_dof_to_GI_c
protected

◆ _dof_to_GII_c

std::unordered_map<const DofObject *, ADReal> BilinearMixedModeCohesiveZoneModel::_dof_to_GII_c
protected

◆ _dof_to_interface_displacement_jump

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

◆ _dof_to_interface_F

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

Map from degree of freedom to interface deformation gradient tensor.

Definition at line 94 of file CohesiveZoneModelBase.h.

Referenced by CohesiveZoneModelBase::computeFandR(), and CohesiveZoneModelBase::initialize().

◆ _dof_to_interface_R

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

Map from degree of freedom to interface rotation tensor.

Definition at line 97 of file CohesiveZoneModelBase.h.

Referenced by CohesiveZoneModelBase::computeFandR(), and CohesiveZoneModelBase::initialize().

◆ _dof_to_lagrange_multiplier

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

◆ _dof_to_lagrange_multiplier_start

std::unordered_map<const DofObject *, 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 105 of file PenaltyWeightedGapUserObject.h.

Referenced by PenaltyWeightedGapUserObject::selfTimestepSetup().

◆ _dof_to_local_czm_traction

std::unordered_map<const DofObject *, RealVectorValue> BilinearMixedModeCohesiveZoneModel::_dof_to_local_czm_traction
protected

Map from degree of freedom to the local-frame cohesive traction.

Stored as raw values (no AD derivatives) since it only feeds auxiliary output.

Definition at line 68 of file BilinearMixedModeCohesiveZoneModel.h.

Referenced by computeCZMTraction(), getCohesiveTractionEffective(), getCohesiveTractionNormal(), getCohesiveTractionTangentialMagnitude(), getCohesiveTractionTangentialOne(), getCohesiveTractionTangentialTwo(), and initialize().

◆ _dof_to_local_penalty

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

◆ _dof_to_local_penalty_friction

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

Map from degree of freedom to local friction penalty value.

Definition at line 127 of file CohesiveZoneModelBase.h.

Referenced by CohesiveZoneModelBase::reinit(), and CohesiveZoneModelBase::timestepSetup().

◆ _dof_to_mode_mixity_ratio

std::unordered_map<const DofObject *, ADReal> BilinearMixedModeCohesiveZoneModel::_dof_to_mode_mixity_ratio
protected

Map from degree of freedom to mode mixity ratio (AD needed?)

Definition at line 64 of file BilinearMixedModeCohesiveZoneModel.h.

Referenced by computeCriticalDisplacementJump(), computeFinalDisplacementJump(), computeModeMixity(), getModeMixityRatio(), and initialize().

◆ _dof_to_normal_pressure

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

◆ _dof_to_normal_strength

std::unordered_map<const DofObject *, ADReal> BilinearMixedModeCohesiveZoneModel::_dof_to_normal_strength
protected

◆ _dof_to_previous_gap

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

◆ _dof_to_real_tangential_velocity

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

◆ _dof_to_rotation_matrix

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

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

Definition at line 76 of file CohesiveZoneModelBase.h.

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

◆ _dof_to_shear_strength

std::unordered_map<const DofObject *, ADReal> BilinearMixedModeCohesiveZoneModel::_dof_to_shear_strength
protected

◆ _dof_to_step_slip

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

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

Definition at line 109 of file CohesiveZoneModelBase.h.

Referenced by CohesiveZoneModelBase::timestepSetup().

◆ _dof_to_tangential_traction

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

◆ _dof_to_weighted_displacements

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

◆ _dof_to_weighted_gap

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

◆ _dof_to_weighted_tangential_velocity

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

◆ _dt

const Real& PenaltyWeightedGapUserObject::_dt
protectedinherited

◆ _dt_old

Real PenaltyWeightedGapUserObject::_dt_old
protectedinherited

previous timestep size

Definition at line 121 of file PenaltyWeightedGapUserObject.h.

Referenced by PenaltyWeightedGapUserObject::selfTimestepSetup().

◆ _epsilon_tolerance

const Real CohesiveZoneModelBase::_epsilon_tolerance
protectedinherited

Tolerance to avoid NaN/Inf in automatic differentiation operations.

Definition at line 130 of file CohesiveZoneModelBase.h.

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

◆ _F_interpolation

ADRankTwoTensor CohesiveZoneModelBase::_F_interpolation
protectedinherited

Deformation gradient for interpolation.

Definition at line 82 of file CohesiveZoneModelBase.h.

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

◆ _F_neighbor_interpolation

ADRankTwoTensor CohesiveZoneModelBase::_F_neighbor_interpolation
protectedinherited

Deformation gradient for interpolation of the neighbor projection.

Definition at line 85 of file CohesiveZoneModelBase.h.

Referenced by CohesiveZoneModelBase::computeQpIProperties(), and CohesiveZoneModelBase::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
protectedinherited

The friction coefficient.

Definition at line 106 of file CohesiveZoneModelBase.h.

Referenced by CohesiveZoneModelBase::reinit().

◆ _frictional_contact_traction_one

ADVariableValue CohesiveZoneModelBase::_frictional_contact_traction_one
protectedinherited

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

Definition at line 118 of file CohesiveZoneModelBase.h.

Referenced by CohesiveZoneModelBase::contactTangentialPressureDirOne(), and CohesiveZoneModelBase::reinit().

◆ _frictional_contact_traction_two

ADVariableValue CohesiveZoneModelBase::_frictional_contact_traction_two
protectedinherited

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

Definition at line 121 of file CohesiveZoneModelBase.h.

Referenced by CohesiveZoneModelBase::contactTangentialPressureDirTwo(), and CohesiveZoneModelBase::reinit().

◆ _GI_c

const MaterialProperty<Real>& BilinearMixedModeCohesiveZoneModel::_GI_c
protected

Fracture parameter mode I.

Definition at line 77 of file BilinearMixedModeCohesiveZoneModel.h.

Referenced by computeQpProperties().

◆ _GI_c_interpolation

ADReal BilinearMixedModeCohesiveZoneModel::_GI_c_interpolation
protected

Interpolated value of fracture paramter mode I.

Definition at line 89 of file BilinearMixedModeCohesiveZoneModel.h.

Referenced by computeQpIProperties(), and computeQpProperties().

◆ _GII_c

const MaterialProperty<Real>& BilinearMixedModeCohesiveZoneModel::_GII_c
protected

Fracture parameter mode II.

Definition at line 80 of file BilinearMixedModeCohesiveZoneModel.h.

Referenced by computeQpProperties().

◆ _GII_c_interpolation

ADReal BilinearMixedModeCohesiveZoneModel::_GII_c_interpolation
protected

Interpolated value of fracture paramter mode II.

Definition at line 92 of file BilinearMixedModeCohesiveZoneModel.h.

Referenced by computeQpIProperties(), and computeQpProperties().

◆ _grad_disp

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

Coupled displacement gradients.

Definition at line 69 of file CohesiveZoneModelBase.h.

Referenced by CohesiveZoneModelBase::computeQpProperties().

◆ _grad_disp_neighbor

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

Coupled displacement and neighbor displacement gradient.

Definition at line 72 of file CohesiveZoneModelBase.h.

Referenced by CohesiveZoneModelBase::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

◆ _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

◆ _max_penalty_multiplier

const Real PenaltyWeightedGapUserObject::_max_penalty_multiplier
protectedinherited

Maximum multiplier applied to the initial penalty factor in AL.

Definition at line 130 of file PenaltyWeightedGapUserObject.h.

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

◆ _mix_mode_criterion

enum BilinearMixedModeCohesiveZoneModel::MixedModeCriterion BilinearMixedModeCohesiveZoneModel::_mix_mode_criterion
protected

◆ _ndisp

const unsigned int CohesiveZoneModelBase::_ndisp
protectedinherited

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 finalize(), CohesiveZoneModelBase::finalize(), WeightedGapUserObject::finalize(), WeightedVelocitiesUserObject::finalize(), and LMWeightedGapUserObject::finalizeNodalScaling().

◆ _normal_strength

const MaterialProperty<Real>& BilinearMixedModeCohesiveZoneModel::_normal_strength
protected

The normal strength material property.

Definition at line 71 of file BilinearMixedModeCohesiveZoneModel.h.

Referenced by computeQpProperties().

◆ _normal_strength_interpolation

ADReal BilinearMixedModeCohesiveZoneModel::_normal_strength_interpolation
protected

Interpolated value of normal_strength.

Definition at line 83 of file BilinearMixedModeCohesiveZoneModel.h.

Referenced by computeQpIProperties(), and computeQpProperties().

◆ _normalization_ptr

const Real* WeightedGapUserObject::_normalization_ptr = nullptr
protectedinherited

Definition at line 255 of file WeightedGapUserObject.h.

◆ _penalty

const Real PenaltyWeightedGapUserObject::_penalty
protectedinherited

◆ _penalty_friction

const Real CohesiveZoneModelBase::_penalty_friction
protectedinherited

The penalty factor for the frictional constraints.

Definition at line 103 of file CohesiveZoneModelBase.h.

Referenced by CohesiveZoneModelBase::reinit(), and CohesiveZoneModelBase::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().

◆ _penalty_stiffness_czm

const Real BilinearMixedModeCohesiveZoneModel::_penalty_stiffness_czm
protected

◆ _penetration_tolerance

const Real PenaltyWeightedGapUserObject::_penetration_tolerance
protectedinherited

◆ _power_law_parameter

const Real BilinearMixedModeCohesiveZoneModel::_power_law_parameter
protected

Power law parameter for bilinear traction model.

Definition at line 105 of file BilinearMixedModeCohesiveZoneModel.h.

Referenced by computeFinalDisplacementJump().

◆ _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

◆ _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 CohesiveZoneModelBase::computeQpIProperties(), and WeightedGapUserObject::computeQpProperties().

◆ _qp_factor

Real WeightedGapUserObject::_qp_factor
protectedinherited

◆ _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().

◆ _regularization_alpha

const Real BilinearMixedModeCohesiveZoneModel::_regularization_alpha
protected

Parameter for the regularization of the Macaulay bracket.

Definition at line 111 of file BilinearMixedModeCohesiveZoneModel.h.

Referenced by computeEffectiveDisplacementJump().

◆ _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.

◆ _set_compressive_traction_to_zero

const bool BilinearMixedModeCohesiveZoneModel::_set_compressive_traction_to_zero
protected

Zero compressive traction.

Definition at line 61 of file BilinearMixedModeCohesiveZoneModel.h.

Referenced by computeCZMTraction().

◆ _shear_strength

const MaterialProperty<Real>& BilinearMixedModeCohesiveZoneModel::_shear_strength
protected

The shear strength material property.

Definition at line 74 of file BilinearMixedModeCohesiveZoneModel.h.

Referenced by computeQpProperties().

◆ _shear_strength_interpolation

ADReal BilinearMixedModeCohesiveZoneModel::_shear_strength_interpolation
protected

Interpolated value of shear_strength.

Definition at line 86 of file BilinearMixedModeCohesiveZoneModel.h.

Referenced by computeQpIProperties(), and computeQpProperties().

◆ _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 109 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(), CohesiveZoneModelBase::prepareJumpKinematicQuantities(), CohesiveZoneModelBase::reinit(), PenaltyFrictionUserObject::reinit(), and PenaltyWeightedGapUserObject::reinit().

◆ _use_physical_gap

const bool PenaltyWeightedGapUserObject::_use_physical_gap
protectedinherited

◆ _viscosity

const Real BilinearMixedModeCohesiveZoneModel::_viscosity
protected

Viscosity for damage model.

Definition at line 108 of file BilinearMixedModeCohesiveZoneModel.h.

Referenced by computeDamage().

◆ _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: