https://mooseframework.inl.gov
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 | List of all members
BilinearMixedModeCohesiveZoneModel Class Referenceabstract

User object that interface pressure resulting from a simple 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 const ADVariableValueczmGlobalTraction (unsigned int i) const
 
virtual void initialize () override
 
virtual void finalize () override
 
virtual void reinit () override
 
virtual void timestepSetup () 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
 
virtual const ADVariableValuecontactPressure () 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. More...
 
virtual Real getNormalGap (const Node *const) const
 
Real physicalGap (const std::pair< ADReal, Real > &gap) const
 Compute physical gap from integration gap quantity. More...
 
ADReal adPhysicalGap (const std::pair< ADReal, Real > &gap) const
 
virtual 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
 
SubProblemgetSubProblem () const
 
bool shouldDuplicateInitialExecution () const
 
virtual Real spatialValue (const Point &) const
 
virtual const std::vector< Point > spatialPoints () const
 
void gatherSum (T &value)
 
void gatherMax (T &value)
 
void gatherMin (T &value)
 
void gatherProxyValueMax (T1 &proxy, T2 &value)
 
void gatherProxyValueMin (T1 &proxy, T2 &value)
 
void setPrimaryThreadCopy (UserObject *primary)
 
UserObjectprimaryThreadCopy ()
 
std::set< UserObjectName > getDependObjects () const
 
virtual bool needThreadedCopy () const
 
const std::set< std::string > & getRequestedItems () override
 
const std::set< std::string > & getSuppliedItems () override
 
unsigned int systemNumber () const
 
virtual bool enabled () const
 
std::shared_ptr< MooseObjectgetSharedPtr ()
 
std::shared_ptr< const MooseObjectgetSharedPtr () const
 
bool isKokkosObject (IsKokkosObjectKey &&) const
 
MooseAppgetMooseApp () 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 InputParametersparameters () 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
 
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 mooseInfo (Args &&... args) const
 
void callMooseError (std::string msg, const bool with_prefix, const hit::Node *node=nullptr) 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 ExecFlagEnumgetExecuteOnEnum () 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 UserObjectgetUserObjectBase (const std::string &param_name, bool is_dependency=true) const
 
const UserObjectgetUserObjectBaseByName (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)
 
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< SubdomainIDgetMaterialPropertyBlocks (const std::string &name)
 
std::vector< SubdomainName > getMaterialPropertyBlockNames (const std::string &name)
 
std::set< BoundaryIDgetMaterialPropertyBoundaryIDs (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 VariableValuecoupledNeighborValue (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 ADVariableValueadCoupledNeighborValue (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 ADVariableValueadCoupledNeighborValueDot (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const ADVariableValue *> adCoupledNeighborValues (const std::string &var_name) const
 
virtual const ADVectorVariableValueadCoupledVectorNeighborValue (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValuecoupledNeighborValueDot (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValuecoupledNeighborValueDotDu (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValuecoupledNeighborValueOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValuecoupledNeighborValueOlder (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableGradientcoupledNeighborGradient (const std::string &var_name, unsigned int comp=0) const
 
virtual std::vector< const VariableGradient *> coupledNeighborGradients (const std::string &var_name) const
 
virtual const VariableGradientcoupledNeighborGradientOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableGradientcoupledNeighborGradientOlder (const std::string &var_name, unsigned int comp=0) const
 
virtual const ADVariableGradientadCoupledNeighborGradient (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableGradientcoupledVectorNeighborGradient (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableGradientcoupledVectorNeighborGradientOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableGradientcoupledVectorNeighborGradientOlder (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableValuecoupledArrayNeighborValue (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableGradientcoupledArrayNeighborGradient (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableGradientcoupledArrayNeighborGradientOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableGradientcoupledArrayNeighborGradientOlder (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableSecondcoupledNeighborSecond (const std::string &var_name, unsigned int i=0) const
 
virtual const VariableValuecoupledNeighborDofValues (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValuecoupledNeighborDofValuesOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValuecoupledNeighborDofValuesOlder (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 ADVariableValuegetADDefaultValue (const std::string &var_name) const
 
const ADVectorVariableValuegetADDefaultVectorValue (const std::string &var_name) const
 
const ADVariableGradientgetADDefaultGradient () const
 
const ADVectorVariableGradientgetADDefaultVectorGradient () const
 
const ADVariableSecondgetADDefaultSecond () const
 
const ADVectorVariableCurlgetADDefaultCurl () 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::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::Scalar< const RealkokkosCoupledDotDu (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)
 
MaterialBasegetMaterial (const std::string &name)
 
MaterialBasegetMaterial (const std::string &name)
 
MaterialBasegetMaterialByName (const std::string &name, bool no_warn=false)
 
MaterialBasegetMaterialByName (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 FunctiongetFunction (const std::string &name) const
 
const FunctiongetFunctionByName (const FunctionName &name) const
 
bool hasFunction (const std::string &param_name) const
 
bool hasFunctionByName (const FunctionName &name) const
 
bool isDefaultPostprocessorValue (const std::string &param_name, const unsigned int index=0) const
 
bool hasPostprocessor (const std::string &param_name, const unsigned int index=0) const
 
bool hasPostprocessorByName (const PostprocessorName &name) const
 
std::size_t coupledPostprocessors (const std::string &param_name) const
 
const PostprocessorName & getPostprocessorName (const std::string &param_name, const unsigned int index=0) const
 
const VectorPostprocessorValuegetVectorPostprocessorValue (const std::string &param_name, const std::string &vector_name) const
 
const VectorPostprocessorValuegetVectorPostprocessorValue (const std::string &param_name, const std::string &vector_name, bool needs_broadcast) const
 
const VectorPostprocessorValuegetVectorPostprocessorValueByName (const VectorPostprocessorName &name, const std::string &vector_name) const
 
const VectorPostprocessorValuegetVectorPostprocessorValueByName (const VectorPostprocessorName &name, const std::string &vector_name, bool needs_broadcast) const
 
const VectorPostprocessorValuegetVectorPostprocessorValueOld (const std::string &param_name, const std::string &vector_name) const
 
const VectorPostprocessorValuegetVectorPostprocessorValueOld (const std::string &param_name, const std::string &vector_name, bool needs_broadcast) const
 
const VectorPostprocessorValuegetVectorPostprocessorValueOldByName (const VectorPostprocessorName &name, const std::string &vector_name) const
 
const VectorPostprocessorValuegetVectorPostprocessorValueOldByName (const VectorPostprocessorName &name, const std::string &vector_name, bool needs_broadcast) const
 
const ScatterVectorPostprocessorValuegetScatterVectorPostprocessorValue (const std::string &param_name, const std::string &vector_name) const
 
const ScatterVectorPostprocessorValuegetScatterVectorPostprocessorValueByName (const VectorPostprocessorName &name, const std::string &vector_name) const
 
const ScatterVectorPostprocessorValuegetScatterVectorPostprocessorValueOld (const std::string &param_name, const std::string &vector_name) const
 
const ScatterVectorPostprocessorValuegetScatterVectorPostprocessorValueOldByName (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)
 
SamplergetSampler (const std::string &name)
 
T & getSamplerByName (const SamplerName &name)
 
SamplergetSamplerByName (const SamplerName &name)
 
virtual void meshChanged ()
 
virtual void meshDisplaced ()
 
PerfGraphperfGraph ()
 
const PostprocessorValuegetPostprocessorValue (const std::string &param_name, const unsigned int index=0) const
 
const PostprocessorValuegetPostprocessorValue (const std::string &param_name, const unsigned int index=0) const
 
const PostprocessorValuegetPostprocessorValueOld (const std::string &param_name, const unsigned int index=0) const
 
const PostprocessorValuegetPostprocessorValueOld (const std::string &param_name, const unsigned int index=0) const
 
const PostprocessorValuegetPostprocessorValueOlder (const std::string &param_name, const unsigned int index=0) const
 
const PostprocessorValuegetPostprocessorValueOlder (const std::string &param_name, const unsigned int index=0) const
 
virtual const PostprocessorValuegetPostprocessorValueByName (const PostprocessorName &name) const
 
virtual const PostprocessorValuegetPostprocessorValueByName (const PostprocessorName &name) const
 
const PostprocessorValuegetPostprocessorValueOldByName (const PostprocessorName &name) const
 
const PostprocessorValuegetPostprocessorValueOldByName (const PostprocessorName &name) const
 
const PostprocessorValuegetPostprocessorValueOlderByName (const PostprocessorName &name) const
 
const PostprocessorValuegetPostprocessorValueOlderByName (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 DistributiongetDistribution (const std::string &name) const
 
const T & getDistribution (const std::string &name) const
 
const DistributiongetDistribution (const std::string &name) const
 
const T & getDistribution (const std::string &name) const
 
const DistributiongetDistributionByName (const DistributionName &name) const
 
const T & getDistributionByName (const std::string &name) const
 
const DistributiongetDistributionByName (const DistributionName &name) const
 
const T & getDistributionByName (const std::string &name) const
 
const Parallel::Communicator & comm () const
 
processor_id_type n_processors () const
 
processor_id_type processor_id () const
 
const std::unordered_map< const DofObject *, std::array< ADReal, 2 > > & dofToWeightedVelocities () const
 Get the degree of freedom to weighted velocities information. More...
 
virtual const ADVariableValuecontactTangentialPressureDirOne () const =0
 
virtual const ADVariableValuecontactTangentialPressureDirTwo () const =0
 
virtual const ADVariableValuecontactTangentialPressureDirOne () const override
 
virtual const ADVariableValuecontactTangentialPressureDirTwo () const override
 

Static Public Member Functions

static InputParameters validParams ()
 
static void callMooseError (MooseApp *const app, const InputParameters &params, std::string msg, const bool with_prefix, const hit::Node *node)
 
static void sort (typename std::vector< T > &vector)
 
static void sortDFS (typename std::vector< T > &vector)
 
static void cyclicDependencyError (CyclicDependencyException< T2 > &e, const std::string &header)
 

Public Attributes

 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  MixedModeCriterion { MixedModeCriterion::POWER_LAW, MixedModeCriterion::BK }
 Mixed-mode propagation criterion. More...
 
enum  AdaptivityNormalPenalty { AdaptivityNormalPenalty::SIMPLE, AdaptivityNormalPenalty::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. More...
 
virtual void computeQpIProperties () override
 Computes properties that are functions both of _qp and _i, for example the weighted gap. More...
 
virtual bool constrainedByOwner () const override
 
virtual void prepareJumpKinematicQuantities () override
 
virtual void computeFandR (const Node *const node) override
 
virtual void computeBilinearMixedModeTraction (const Node *const node) override
 Encapsulate the CZM constitutive behavior. More...
 
virtual void computeGlobalTraction (const Node *const node) override
 Compute global traction for mortar application. More...
 
template<class T >
normalizeQuantity (const std::unordered_map< const DofObject *, T > &map, const Node *const node)
 Normalize mortar quantities (remove mortar integral scaling) More...
 
virtual const VariableTestValuetest () const override
 
void selfInitialize ()
 
void selfFinalize ()
 
void selfTimestepSetup ()
 
void bussettaAdaptivePenalty (const Real previous_gap, const Real gap, Real &penalty, Real &eval_tn)
 Adaptive, local penalty for AL. More...
 
void adaptiveNormalPenalty (const Real previous_gap, const Real gap, Real &penalty)
 See Algorithm 3 of 'The adapted augmented Lagrangian method: a new method for the resolution of the mechanical frictional contact problem', Comput Mech (2012) 49: 259-275. More...
 
template<typename K , typename V >
findValue (const std::unordered_map< K, V > &map, const K &key, const V &default_value=0) const
 Find a value in a map or return a default if the key doesn't exist. More...
 
virtual void addPostprocessorDependencyHelper (const PostprocessorName &name) const override
 
virtual void addVectorPostprocessorDependencyHelper (const VectorPostprocessorName &name) const override
 
virtual void addUserObjectDependencyHelper (const UserObject &uo) const override
 
void addReporterDependencyHelper (const ReporterName &reporter_name) override
 
void flagInvalidSolutionInternal (const InvalidSolutionID invalid_solution_id) const
 
InvalidSolutionID registerInvalidSolutionInternal (const std::string &message, const bool warning) const
 
const ReporterContextBasegetReporterContextBaseByName (const ReporterName &reporter_name) const
 
const ReporterNamegetReporterName (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 VariableValuecoupledScalarValue (const std::string &var_name, unsigned int comp=0) const
 
const ADVariableValueadCoupledScalarValue (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 VariableValuecoupledVectorTagScalarValue (const std::string &var_name, TagID tag, unsigned int comp=0) const
 
const VariableValuecoupledMatrixTagScalarValue (const std::string &var_name, TagID tag, unsigned int comp=0) const
 
const VariableValuecoupledScalarValueOld (const std::string &var_name, unsigned int comp=0) const
 
const VariableValuecoupledScalarValueOlder (const std::string &var_name, unsigned int comp=0) const
 
const VariableValuecoupledScalarDot (const std::string &var_name, unsigned int comp=0) const
 
const ADVariableValueadCoupledScalarDot (const std::string &var_name, unsigned int comp=0) const
 
const VariableValuecoupledScalarDotDot (const std::string &var_name, unsigned int comp=0) const
 
const VariableValuecoupledScalarDotOld (const std::string &var_name, unsigned int comp=0) const
 
const VariableValuecoupledScalarDotDotOld (const std::string &var_name, unsigned int comp=0) const
 
const VariableValuecoupledScalarDotDu (const std::string &var_name, unsigned int comp=0) const
 
const VariableValuecoupledScalarDotDotDu (const std::string &var_name, unsigned int comp=0) const
 
const MooseVariableScalargetScalarVar (const std::string &var_name, unsigned int comp) const
 
const std::set< SubdomainID > & getHigherDimSubdomainIDs () const
 
const std::set< BoundaryID > & getBoundaryIDs () const
 
const AutomaticMortarGenerationamg () 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 intcoupledIndices (const std::string &var_name) const
 
virtual const VariableValuecoupledValue (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 VariableValuecoupledValueLower (const std::string &var_name, unsigned int comp=0) const
 
const ADVariableValueadCoupledValue (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const ADVariableValue *> adCoupledValues (const std::string &var_name) const
 
const ADVariableValueadCoupledLowerValue (const std::string &var_name, unsigned int comp=0) const
 
const ADVectorVariableValueadCoupledVectorValue (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const ADVectorVariableValue *> adCoupledVectorValues (const std::string &var_name) const
 
virtual const VariableValuecoupledVectorTagValue (const std::string &var_names, TagID tag, unsigned int index=0) const
 
virtual const VariableValuecoupledVectorTagValue (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 ArrayVariableValuecoupledVectorTagArrayValue (const std::string &var_names, TagID tag, unsigned int index=0) const
 
virtual const ArrayVariableValuecoupledVectorTagArrayValue (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 VariableGradientcoupledVectorTagGradient (const std::string &var_names, TagID tag, unsigned int index=0) const
 
virtual const VariableGradientcoupledVectorTagGradient (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 ArrayVariableGradientcoupledVectorTagArrayGradient (const std::string &var_names, TagID tag, unsigned int index=0) const
 
virtual const ArrayVariableGradientcoupledVectorTagArrayGradient (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 VariableValuecoupledVectorTagDofValue (const std::string &var_name, TagID tag, unsigned int index=0) const
 
virtual const VariableValuecoupledVectorTagDofValue (const std::string &var_names, const std::string &tag_name, unsigned int index=0) const
 
const ArrayVariableValuecoupledVectorTagArrayDofValue (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 VariableValuecoupledMatrixTagValue (const std::string &var_names, TagID tag, unsigned int index=0) const
 
virtual const VariableValuecoupledMatrixTagValue (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 VectorVariableValuecoupledVectorValue (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableValuecoupledArrayValue (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const ArrayVariableValue *> coupledArrayValues (const std::string &var_name) const
 
MooseWritableVariablewritableVariable (const std::string &var_name, unsigned int comp=0)
 
virtual VariableValuewritableCoupledValue (const std::string &var_name, unsigned int comp=0)
 
void checkWritableVar (MooseWritableVariable *var)
 
virtual const VariableValuecoupledValueOld (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 VariableValuecoupledValueOlder (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const VariableValue *> coupledValuesOlder (const std::string &var_name) const
 
virtual const VariableValuecoupledValuePreviousNL (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableValuecoupledVectorValueOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableValuecoupledVectorValueOlder (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableValuecoupledArrayValueOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableValuecoupledArrayValueOlder (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableGradientcoupledGradient (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const VariableGradient *> coupledGradients (const std::string &var_name) const
 
const ADVariableGradientadCoupledGradient (const std::string &var_name, unsigned int comp=0) const
 
const ADVariableGradientadCoupledGradientDot (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 ADVectorVariableGradientadCoupledVectorGradient (const std::string &var_name, unsigned int comp=0) const
 
const ADVariableSecondadCoupledSecond (const std::string &var_name, unsigned int comp=0) const
 
const ADVectorVariableSecondadCoupledVectorSecond (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableGradientcoupledGradientOld (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const VariableGradient *> coupledGradientsOld (const std::string &var_name) const
 
virtual const VariableGradientcoupledGradientOlder (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableGradientcoupledGradientPreviousNL (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableGradientcoupledGradientDot (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableGradientcoupledGradientDotDot (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableGradientcoupledVectorGradient (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableGradientcoupledVectorGradientOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableGradientcoupledVectorGradientOlder (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableGradientcoupledArrayGradient (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableGradientcoupledArrayGradientOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableGradientcoupledArrayGradientOlder (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableGradientcoupledArrayGradientDot (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableCurlcoupledCurl (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableCurlcoupledCurlOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableCurlcoupledCurlOlder (const std::string &var_name, unsigned int comp=0) const
 
const ADVectorVariableCurladCoupledCurl (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableDivergencecoupledDiv (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableDivergencecoupledDivOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableDivergencecoupledDivOlder (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableSecondcoupledSecond (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableSecondcoupledSecondOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableSecondcoupledSecondOlder (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableSecondcoupledSecondPreviousNL (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValuecoupledDot (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const VariableValue *> coupledDots (const std::string &var_name) const
 
virtual const VariableValuecoupledDotDot (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValuecoupledDotOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValuecoupledDotDotOld (const std::string &var_name, unsigned int comp=0) const
 
const ADVariableValueadCoupledDot (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const ADVariableValue *> adCoupledDots (const std::string &var_name) const
 
const ADVariableValueadCoupledDotDot (const std::string &var_name, unsigned int comp=0) const
 
const ADVectorVariableValueadCoupledVectorDot (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableValuecoupledVectorDot (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableValuecoupledVectorDotDot (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableValuecoupledVectorDotOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableValuecoupledVectorDotDotOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValuecoupledVectorDotDu (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValuecoupledVectorDotDotDu (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableValuecoupledArrayDot (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableValuecoupledArrayDotDot (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableValuecoupledArrayDotOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableValuecoupledArrayDotDotOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValuecoupledDotDu (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValuecoupledDotDotDu (const std::string &var_name, unsigned int comp=0) const
 
const VariableValuecoupledArrayDotDu (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 >::typeadCoupledNodalValue (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 VariableValuecoupledNodalDotDot (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValuecoupledNodalDotOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValuecoupledNodalDotDotOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValuecoupledDofValues (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const VariableValue *> coupledAllDofValues (const std::string &var_name) const
 
virtual const VariableValuecoupledDofValuesOld (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const VariableValue *> coupledAllDofValuesOld (const std::string &var_name) const
 
virtual const VariableValuecoupledDofValuesOlder (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const VariableValue *> coupledAllDofValuesOlder (const std::string &var_name) const
 
virtual const ArrayVariableValuecoupledArrayDofValues (const std::string &var_name, unsigned int comp=0) const
 
virtual const ADVariableValueadCoupledDofValues (const std::string &var_name, unsigned int comp=0) const
 
const ADVariableValueadZeroValue () const
 
const ADVariableGradientadZeroGradient () const
 
const ADVariableSecondadZeroSecond () 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 MooseVariableFieldBasegetFEVar (const std::string &var_name, unsigned int comp) const
 
const MooseVariableFieldBasegetFieldVar (const std::string &var_name, unsigned int comp) const
 
MooseVariableFieldBasegetFieldVar (const std::string &var_name, unsigned int comp)
 
const T * getVarHelper (const std::string &var_name, unsigned int comp) const
 
T * getVarHelper (const std::string &var_name, unsigned int comp)
 
MooseVariablegetVar (const std::string &var_name, unsigned int comp)
 
const MooseVariablegetVar (const std::string &var_name, unsigned int comp) const
 
VectorMooseVariablegetVectorVar (const std::string &var_name, unsigned int comp)
 
const VectorMooseVariablegetVectorVar (const std::string &var_name, unsigned int comp) const
 
ArrayMooseVariablegetArrayVar (const std::string &var_name, unsigned int comp)
 
const ArrayMooseVariablegetArrayVar (const std::string &var_name, unsigned int comp) const
 
void validateExecutionerType (const std::string &name, const std::string &fn_name) const
 
std::vector< T > coupledVectorHelper (const std::string &var_name, const Func &func) const
 
const T & getReporterValue (const std::string &param_name, const std::size_t time_index=0)
 
const T & getReporterValue (const std::string &param_name, ReporterMode mode, const std::size_t time_index=0)
 
const T & getReporterValue (const std::string &param_name, const std::size_t time_index=0)
 
const T & getReporterValue (const std::string &param_name, ReporterMode mode, const std::size_t time_index=0)
 
const T & getReporterValueByName (const ReporterName &reporter_name, const std::size_t time_index=0)
 
const T & getReporterValueByName (const ReporterName &reporter_name, ReporterMode mode, const std::size_t time_index=0)
 
const T & getReporterValueByName (const ReporterName &reporter_name, const std::size_t time_index=0)
 
const T & getReporterValueByName (const ReporterName &reporter_name, ReporterMode mode, const std::size_t time_index=0)
 
bool hasReporterValue (const std::string &param_name) const
 
bool hasReporterValue (const std::string &param_name) const
 
bool hasReporterValue (const std::string &param_name) const
 
bool hasReporterValue (const std::string &param_name) const
 
bool hasReporterValueByName (const ReporterName &reporter_name) const
 
bool hasReporterValueByName (const ReporterName &reporter_name) const
 
bool hasReporterValueByName (const ReporterName &reporter_name) const
 
bool hasReporterValueByName (const ReporterName &reporter_name) const
 
const GenericMaterialProperty< T, is_ad > * defaultGenericMaterialProperty (const std::string &name)
 
const GenericMaterialProperty< T, is_ad > * defaultGenericMaterialProperty (const std::string &name)
 
const MaterialProperty< T > * defaultMaterialProperty (const std::string &name)
 
const MaterialProperty< T > * defaultMaterialProperty (const std::string &name)
 
const ADMaterialProperty< T > * defaultADMaterialProperty (const std::string &name)
 
const ADMaterialProperty< T > * defaultADMaterialProperty (const std::string &name)
 
void selfInitialize ()
 
void selfTimestepSetup ()
 
virtual const VariableTestValuetest () 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)
 

Static Protected Member Functions

static std::string meshPropertyName (const std::string &data_name, const std::string &prefix)
 
static void trimDerivative (dof_id_type remove_derivative_index, ADReal &dual_number)
 
static void trimInteriorNodeDerivatives (const std::map< unsigned int, unsigned int > &primary_ip_lowerd_map, const Variables &moose_var, DualNumbers &ad_vars, const bool is_secondary)
 

Protected Attributes

const unsigned int _ndisp
 Number of displacement components. More...
 
std::vector< const GenericVariableGradient< true > * > _grad_disp
 Coupled displacement gradients. More...
 
std::vector< const GenericVariableGradient< true > * > _grad_disp_neighbor
 Coupled displacement and neighbor displacement gradient. More...
 
std::unordered_map< const DofObject *, ADRankTwoTensor_dof_to_rotation_matrix
 *** Kinematics/displacement jump quantities *** Map from degree of freedom to rotation matrix More...
 
std::unordered_map< const DofObject *, ADRealVectorValue_dof_to_interface_displacement_jump
 Map from degree of freedom to local displacement jump. More...
 
ADRankTwoTensor _F_interpolation
 Deformation gradient for interpolation. More...
 
ADRankTwoTensor _F_neighbor_interpolation
 Deformation gradient for interpolation of the neighbor projection. More...
 
std::unordered_map< const DofObject *, ADRankTwoTensor_dof_to_F
 Map from degree of freedom to secondary, interpolated deformation gradient tensor. More...
 
std::unordered_map< const DofObject *, ADRankTwoTensor_dof_to_F_neighbor
 Map from degree of freedom to neighbor, interpolated deformation gradient tensor. More...
 
std::unordered_map< const DofObject *, ADRankTwoTensor_dof_to_interface_F
 Map from degree of freedom to interface deformation gradient tensor. More...
 
std::unordered_map< const DofObject *, ADRankTwoTensor_dof_to_interface_R
 Map from degree of freedom to interface rotation tensor. More...
 
std::unordered_map< const DofObject *, ADReal_dof_to_mode_mixity_ratio
 Map from degree of freedom to mode mixity ratio (AD needed?) More...
 
const MaterialProperty< Real > & _normal_strength
 The normal strength material property. More...
 
const MaterialProperty< Real > & _shear_strength
 The shear strength material property. More...
 
const MaterialProperty< Real > & _GI_c
 Fracture parameter mode I. More...
 
const MaterialProperty< Real > & _GII_c
 Fracture parameter mode II. More...
 
ADReal _normal_strength_interpolation
 Interpolated value of normal_strength. More...
 
ADReal _shear_strength_interpolation
 Interpolated value of shear_strength. More...
 
ADReal _GI_c_interpolation
 Interpolated value of fracture paramter mode I. More...
 
ADReal _GII_c_interpolation
 Interpolated value of fracture paramter mode II. More...
 
const Real _penalty_stiffness_czm
 
enum BilinearMixedModeCohesiveZoneModel::MixedModeCriterion _mix_mode_criterion
 
const Real _power_law_parameter
 Power law parameter for bilinear traction model. More...
 
const Real _viscosity
 Viscosity for damage model. More...
 
const Real _regularization_alpha
 Parameter for the regularization of the Macaulay bracket. More...
 
std::vector< ADVariableValue_czm_interpolated_traction
 The global traction. More...
 
std::unordered_map< const DofObject *, ADRealVectorValue_dof_to_czm_traction
 Total Lagrangian stress to be applied on CZM interface. More...
 
std::unordered_map< const DofObject *, ADRealVectorValue_dof_to_displacement_jump
 
const Real _penalty
 The penalty factor. More...
 
const Real _penalty_multiplier
 penalty growth factor for augmented Lagrange More...
 
const Real _penetration_tolerance
 penetration tolerance for augmented Lagrange contact More...
 
ADVariableValue _contact_pressure
 The contact pressure on the mortar segument quadrature points. More...
 
std::unordered_map< const DofObject *, ADReal_dof_to_normal_pressure
 Map from degree of freedom to normal pressure for reporting. More...
 
std::unordered_map< const DofObject *, Real_dof_to_lagrange_multiplier
 Map from degree of freedom to augmented lagrange multiplier. More...
 
std::unordered_map< const DofObject *, Real_dof_to_local_penalty
 Map from degree of freedom to local penalty value. More...
 
std::unordered_map< const DofObject *, Real_dof_to_previous_gap
 Map from degree of freedom to previous AL iteration gap values. More...
 
const Real_dt
 Current delta t... or timestep size. More...
 
Real _dt_old
 previous timestep size More...
 
const bool _use_physical_gap
 Use scaled or physical gap. More...
 
const MooseVariable *const _aux_lm_var
 The auxiliary Lagrange multiplier variable (used together whith the Petrov-Galerkin approach) More...
 
const Real _max_penalty_multiplier
 Maximum multiplier applied to the initial penalty factor in AL. More...
 
enum PenaltyWeightedGapUserObject::AdaptivityNormalPenalty _adaptivity_normal
 
FEProblemBase_fe_problem
 The base finite element problem. More...
 
ADReal _qp_gap
 The value of the gap at the current quadrature point. More...
 
Real _qp_factor
 The value of the LM at the current quadrature point. More...
 
const bool _nodal
 Whether the dof objects are nodal; if they're not, then they're elemental. More...
 
const MooseVariable *const _disp_x_var
 The x displacement variable. More...
 
const MooseVariable *const _disp_y_var
 The y displacement variable. More...
 
const bool _has_disp_z
 For 2D mortar contact no displacement will be specified, so const pointers used. More...
 
const MooseVariable *const _disp_z_var
 The z displacement variable. More...
 
const ADVariableValue_secondary_disp_x
 x-displacement on the secondary face More...
 
const ADVariableValue_primary_disp_x
 x-displacement on the primary face More...
 
const ADVariableValue_secondary_disp_y
 y-displacement on the secondary face More...
 
const ADVariableValue_primary_disp_y
 y-displacement on the primary face More...
 
const ADVariableValue *const _secondary_disp_z
 z-displacement on the secondary face More...
 
const ADVariableValue *const _primary_disp_z
 z-displacement on the primary face More...
 
const MooseArray< Real > & _coord
 Member for handling change of coordinate systems (xyz, rz, spherical) More...
 
ADRealVectorValue _qp_gap_nodal
 Vector for computation of weighted gap with nodal normals. More...
 
ADRealVectorValue _qp_displacement_nodal
 Vector for computation of relative displacement (determines mixity ratio in interface problems) More...
 
std::unordered_map< const DofObject *, std::pair< ADReal, Real > > _dof_to_weighted_gap
 A map from node to weighted gap and normalization (if requested) More...
 
std::unordered_map< const DofObject *, ADRealVectorValue_dof_to_weighted_displacements
 A map from node to weighted displacements. More...
 
const ADReal_weighted_gap_ptr = nullptr
 A pointer members that can be used to help avoid copying ADReals. More...
 
const Real_normalization_ptr = nullptr
 
const VariableTestValue_test = nullptr
 A pointer to the test function associated with the weighted gap. More...
 
bool _is_weighted_gap_nodal = true
 Whether the weighted gap is associated with nodes or elements (like for a CONSTANT MONOMIAL Lagrange multiplier). More...
 
unsigned int _qp = 0
 Quadrature point index for the mortar segments. More...
 
unsigned int _i = 0
 Test function index. More...
 
SubProblem_subproblem
 
SystemBase_sys
 
const THREAD_ID _tid
 
Assembly_assembly
 
const Moose::CoordinateSystemType_coord_sys
 
const bool _duplicate_initial_execution
 
std::set< std::string > _depend_uo
 
const bool & _enabled
 
MooseApp_app
 
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. More...
 
std::unordered_map< const DofObject *, std::array< ADReal, 2 > > _dof_to_real_tangential_velocity
 A map from node to two interpolated, physical tangential velocities. More...
 
std::array< const ADReal *, 2 > _tangential_vel_ptr = {{nullptr, nullptr}}
 An array of two pointers to avoid copies. More...
 
std::array< ADReal, 2 > _qp_tangential_velocity
 The value of the tangential velocity values at the current quadrature point. More...
 
std::array< ADReal, 2 > _qp_real_tangential_velocity
 The value of the "real" tangential velocity values at the current quadrature point. More...
 
ADRealVectorValue _qp_tangential_velocity_nodal
 The value of the tangential velocity vectors at the current node. More...
 
ADRealVectorValue _qp_real_tangential_velocity_nodal
 The value of the real tangential velocity vectors at the current node. More...
 
SystemBase_sys
 Reference to the EquationSystem object. More...
 
MooseVariableField< Real > & _secondary_var
 Reference to the secondary variable. More...
 
MooseVariableField< Real > & _primary_var
 Reference to the primary variable. More...
 
const ADVariableValue_secondary_x_dot
 x-velocity on the secondary face More...
 
const ADVariableValue_primary_x_dot
 x-velocity on the primary face More...
 
const ADVariableValue_secondary_y_dot
 y-velocity on the secondary face More...
 
const ADVariableValue_primary_y_dot
 y-velocity on the primary face More...
 
const ADVariableValue *const _secondary_z_dot
 z-velocity on the secondary face More...
 
const ADVariableValue *const _primary_z_dot
 z-velocity on the primary face More...
 
bool _3d
 Automatic flag to determine whether we are doing three-dimensional work. More...
 
const Real _penalty
 The normal penalty factor. More...
 
const Real _penalty_friction
 The penalty factor for the frictional constraints. More...
 
const Real _friction_coefficient
 The friction coefficient. More...
 
std::unordered_map< const DofObject *, std::pair< TwoVector, TwoVector > > _dof_to_step_slip
 Map from degree of freedom to current and old step slip. More...
 
std::unordered_map< const DofObject *, std::pair< TwoVector, TwoVector > > _dof_to_accumulated_slip
 Map from degree of freedom to current and old accumulated slip. More...
 
std::unordered_map< const DofObject *, std::pair< ADTwoVector, TwoVector > > _dof_to_tangential_traction
 Map from degree of freedom to current and old tangential traction. More...
 
std::unordered_map< const DofObject *, ADReal_dof_to_czm_normal_traction
 Map from degree of freedom to czm normal traction. More...
 
ADVariableValue _frictional_contact_traction_one
 The first frictional contact pressure on the mortar segment quadrature points. More...
 
ADVariableValue _frictional_contact_traction_two
 The second frictional contact pressure on the mortar segment quadrature points. More...
 
std::unordered_map< const DofObject *, TwoVector_dof_to_frictional_lagrange_multipliers
 Map from degree of freedom to augmented lagrange multiplier. More...
 
std::unordered_map< const DofObject *, Real_dof_to_local_penalty_friction
 Map from degree of freedom to local friction penalty value. More...
 
const Real _epsilon_tolerance
 Tolerance to avoid NaN/Inf in automatic differentiation operations. More...
 
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
 
std::unordered_map< const DofObject *, std::pair< ADReal, Real > > _dof_to_damage
 

Static Protected Attributes

static const std::string _interpolated_old
 
static const std::string _interpolated_older
 
AugmentedLagrangianContactProblemInterface *const _augmented_lagrange_problem
 augmented Lagrange problem and iteration number More...
 
const unsigned int_lagrangian_iteration_number
 
static const unsigned int _no_iterations = 0
 

Detailed Description

User object that interface pressure resulting from a simple traction separation law.

Definition at line 18 of file BilinearMixedModeCohesiveZoneModel.h.

Member Enumeration Documentation

◆ AdaptivityNormalPenalty

The adaptivity method for the penalty factor at augmentations.

Enumerator
SIMPLE 
BUSSETTA 

Definition at line 121 of file PenaltyWeightedGapUserObject.h.

121 { SIMPLE, BUSSETTA } _adaptivity_normal;
Executioner set up to solve a thermal-hydraulics problem using the SIMPLE algorithm.
Definition: SIMPLE.h:19
enum PenaltyWeightedGapUserObject::AdaptivityNormalPenalty _adaptivity_normal

◆ MixedModeCriterion

Mixed-mode propagation criterion.

Enumerator
POWER_LAW 
BK 

Definition at line 135 of file BilinearMixedModeCohesiveZoneModel.h.

136  {
137  POWER_LAW,
138  BK
enum BilinearMixedModeCohesiveZoneModel::MixedModeCriterion _mix_mode_criterion

Constructor & Destructor Documentation

◆ BilinearMixedModeCohesiveZoneModel()

BilinearMixedModeCohesiveZoneModel::BilinearMixedModeCohesiveZoneModel ( const InputParameters parameters)

Definition at line 69 of file BilinearMixedModeCohesiveZoneModel.C.

75  _ndisp(coupledComponents("displacements")),
76  _normal_strength(getMaterialProperty<Real>("normal_strength")),
77  _shear_strength(getMaterialProperty<Real>("shear_strength")),
78  _GI_c(getMaterialProperty<Real>("GI_c")),
79  _GII_c(getMaterialProperty<Real>("GII_c")),
80  _penalty_stiffness_czm(getParam<Real>("penalty_stiffness")),
81  _mix_mode_criterion(getParam<MooseEnum>("mixed_mode_criterion").getEnum<MixedModeCriterion>()),
82  _power_law_parameter(getParam<Real>("power_law_parameter")),
83  _viscosity(getParam<Real>("viscosity")),
84  _regularization_alpha(getParam<Real>("regularization_alpha"))
85 {
87 
88  for (unsigned int i = 0; i < _ndisp; ++i)
89  {
90  _grad_disp.push_back(&adCoupledGradient("displacements", i));
91  _grad_disp_neighbor.push_back(&adCoupledGradient("displacements", i));
92  }
93 
94  // Set non-intervening components to zero
95  for (unsigned int i = _ndisp; i < 3; i++)
96  {
97  _grad_disp.push_back(&_ad_grad_zero);
99  }
100 
101  // Checks for bilinear traction input.
102  if (!isParamValid("GI_c") || !isParamValid("GII_c") || !isParamValid("normal_strength") ||
103  !isParamValid("shear_strength") || !isParamValid("power_law_parameter") ||
104  !isParamValid("penalty_stiffness"))
105  paramError("GI_c",
106  "The CZM bilinear mixed mode traction parameters GI_c, GII_c, normal_strength, "
107  "shear_strength, and power_law_parameter are required. Revise your input and add "
108  "those parameters if you want to use the bilinear mixed mode traction model. ");
109 }
void paramError(const std::string &param, Args... args) const
const MooseArray< ADRealVectorValue > & _ad_grad_zero
const MaterialProperty< Real > & _GI_c
Fracture parameter mode I.
const Real _regularization_alpha
Parameter for the regularization of the Macaulay bracket.
WeightedGapUserObject(const InputParameters &parameters)
const InputParameters & parameters() const
const ADVariableGradient & adCoupledGradient(const std::string &var_name, unsigned int comp=0) const
std::vector< const GenericVariableGradient< true > * > _grad_disp
Coupled displacement gradients.
PenaltySimpleCohesiveZoneModel(const InputParameters &parameters)
const unsigned int _ndisp
Number of displacement components.
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.
std::vector< ADVariableValue > _czm_interpolated_traction
The global traction.
const MaterialProperty< Real > & _normal_strength
The normal strength material property.
PenaltyWeightedGapUserObject(const InputParameters &parameters)
unsigned int coupledComponents(const std::string &var_name) const
WeightedVelocitiesUserObject(const InputParameters &parameters)
std::vector< const GenericVariableGradient< true > * > _grad_disp_neighbor
Coupled displacement and neighbor displacement gradient.
enum BilinearMixedModeCohesiveZoneModel::MixedModeCriterion _mix_mode_criterion
bool isParamValid(const std::string &name) const
const Real _power_law_parameter
Power law parameter for bilinear traction model.

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

Referenced by PenaltyWeightedGapUserObject::bussettaAdaptivePenalty().

369 {
370  using std::abs, std::max, std::sqrt;
371 
372  const bool condition_one = abs(abs(previous_gap / gap) - 1.0) < 0.01;
373  const bool condition_two = abs(gap) > 1.01 * abs(previous_gap);
374 
375  if (previous_gap * gap < 0)
376  {
377  if (previous_gap > _penetration_tolerance)
378  penalty = abs(penalty * previous_gap / gap * (abs(gap) + _penetration_tolerance) /
379  (gap - previous_gap));
380  else
381  penalty = abs(penalty * previous_gap / 10.0 / gap);
382  }
383  else if (abs(gap) > _penetration_tolerance)
384  {
385  if (abs(gap - previous_gap) >
386  max(gap / 10.0, max(previous_gap / 10.0, 5.0 * _penetration_tolerance)))
387  penalty *= 10.0;
388  else if (condition_one && gap < 10.0 * _penetration_tolerance)
389  penalty *= MathUtils::pow(sqrt(abs(gap) / _penetration_tolerance - 1.0) + 1.0, 2);
390  else if (condition_two)
391  penalty *= 10.0 * previous_gap / gap;
392  else
393  penalty *= sqrt(abs(gap) / _penetration_tolerance - 1.0) + 1.0;
394  }
395 }
MetaPhysicL::DualNumber< V, D, asd > abs(const MetaPhysicL::DualNumber< V, D, asd > &a)
auto max(const L &left, const R &right)
const Real _penetration_tolerance
penetration tolerance for augmented Lagrange contact
CTSub CT_OPERATOR_BINARY CTMul CTCompareLess CTCompareGreater CTCompareEqual _arg template * sqrt(_arg)) *_arg.template D< dtag >()) CT_SIMPLE_UNARY_FUNCTION(tanh
T pow(T x, int e)

◆ adPhysicalGap()

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

Definition at line 59 of file WeightedGapUserObject.h.

Referenced by PenaltyWeightedGapUserObject::selfFinalize().

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

◆ augmentedLagrangianSetup()

void PenaltyWeightedGapUserObject::augmentedLagrangianSetup ( )
overridevirtualinherited

Implements AugmentedLagrangeInterface.

Definition at line 287 of file PenaltyWeightedGapUserObject.C.

288 {
289  // Loop over all nodes for which a gap has been computed
290  for (auto & [dof_object, wgap] : _dof_to_weighted_gap)
291  {
292  const Real gap = physicalGap(wgap);
293  // Store previous augmented lagrange iteration gap
294  _dof_to_previous_gap[dof_object] = gap;
295  }
296 }
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
std::unordered_map< const DofObject *, Real > _dof_to_previous_gap
Map from degree of freedom to previous AL iteration gap values.
Real physicalGap(const std::pair< ADReal, Real > &gap) const
Compute physical gap from integration gap quantity.
std::unordered_map< const DofObject *, std::pair< ADReal, Real > > _dof_to_weighted_gap
A map from node to weighted gap and normalization (if requested)

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

Referenced by PenaltyWeightedGapUserObject::updateAugmentedLagrangianMultipliers().

347 {
348  // Positive gaps means no contact
349  if (previous_gap > 0.0)
350  {
351  penalty = _penalty;
352  eval_tn = 0.0;
353  }
354  else
355  {
356  if (previous_gap * gap < 0)
357  eval_tn = 0.0;
358  else
359  eval_tn = penalty * previous_gap;
360 
361  adaptiveNormalPenalty(previous_gap, gap, penalty);
362  }
363 }
const Real _penalty
The penalty factor.
void adaptiveNormalPenalty(const Real previous_gap, const Real gap, Real &penalty)
See Algorithm 3 of &#39;The adapted augmented Lagrangian method: a new method for the resolution of the m...

◆ computeBilinearMixedModeTraction()

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

Encapsulate the CZM constitutive behavior.

Reimplemented from PenaltySimpleCohesiveZoneModel.

Definition at line 269 of file BilinearMixedModeCohesiveZoneModel.C.

270 {
271  using std::max, std::min;
272 
273  // First call does not have maps available
274  const bool return_boolean = _dof_to_weighted_gap.find(node) == _dof_to_weighted_gap.end();
275  if (return_boolean)
276  return;
277 
278  computeModeMixity(node);
282  computeDamage(node);
283 
284  // Split displacement jump into active and inactive parts
285  const auto interface_displacement_jump =
287 
288  const ADRealVectorValue delta_active(max(interface_displacement_jump(0), 0.0),
289  interface_displacement_jump(1),
290  interface_displacement_jump(2));
291  const ADRealVectorValue delta_inactive(min(interface_displacement_jump(0), 0.0), 0.0, 0.0);
292 
293  // This traction vector is local at this point.
294  _dof_to_czm_traction[node] =
295  -(1.0 - _dof_to_damage[node].first) * _penalty_stiffness_czm * delta_active -
296  _penalty_stiffness_czm * delta_inactive;
297 }
std::unordered_map< const DofObject *, ADRealVectorValue > _dof_to_interface_displacement_jump
Map from degree of freedom to local displacement jump.
std::unordered_map< const DofObject *, std::pair< ADReal, Real > > _dof_to_damage
auto max(const L &left, const R &right)
std::unordered_map< const DofObject *, ADRealVectorValue > _dof_to_czm_traction
Total Lagrangian stress to be applied on CZM interface.
virtual void computeModeMixity(const Node *const node)
T normalizeQuantity(const std::unordered_map< const DofObject *, T > &map, const Node *const node)
Normalize mortar quantities (remove mortar integral scaling)
virtual void computeDamage(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)
std::unordered_map< const DofObject *, std::pair< ADReal, Real > > _dof_to_weighted_gap
A map from node to weighted gap and normalization (if requested)
auto min(const L &left, const R &right)

◆ computeCriticalDisplacementJump()

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

Definition at line 334 of file BilinearMixedModeCohesiveZoneModel.C.

Referenced by computeBilinearMixedModeTraction().

335 {
336  using std::sqrt;
337 
338  const auto interface_displacement_jump =
340 
341  const auto mixity_ratio = libmesh_map_find(_dof_to_mode_mixity_ratio, node);
342 
343  const auto delta_normal_knot =
345  const auto delta_shear_knot =
347 
348  _dof_to_delta_initial[node] = delta_shear_knot;
349 
350  if (interface_displacement_jump(0) > _epsilon_tolerance)
351  {
352  const auto delta_mixed = sqrt(delta_shear_knot * delta_shear_knot +
353  Utility::pow<2>(mixity_ratio * delta_normal_knot));
354 
355  _dof_to_delta_initial[node] = delta_normal_knot * delta_shear_knot *
356  sqrt(1.0 + mixity_ratio * mixity_ratio) / delta_mixed;
357  }
358 }
std::unordered_map< const DofObject *, ADRealVectorValue > _dof_to_interface_displacement_jump
Map from degree of freedom to local displacement jump.
std::unordered_map< const DofObject *, ADReal > _dof_to_normal_strength
std::unordered_map< const DofObject *, ADReal > _dof_to_delta_initial
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)
std::unordered_map< const DofObject *, ADReal > _dof_to_shear_strength
CTSub CT_OPERATOR_BINARY CTMul CTCompareLess CTCompareGreater CTCompareEqual _arg template * sqrt(_arg)) *_arg.template D< dtag >()) CT_SIMPLE_UNARY_FUNCTION(tanh
std::unordered_map< const DofObject *, ADReal > _dof_to_mode_mixity_ratio
Map from degree of freedom to mode mixity ratio (AD needed?)

◆ computeDamage()

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

Definition at line 418 of file BilinearMixedModeCohesiveZoneModel.C.

Referenced by computeBilinearMixedModeTraction().

419 {
420  const auto delta_max = libmesh_map_find(_dof_to_delta_max, node);
421  const auto delta_initial = libmesh_map_find(_dof_to_delta_initial, node);
422  const auto delta_final = libmesh_map_find(_dof_to_delta_final, node);
423 
424  if (delta_max < delta_initial)
425  _dof_to_damage[node].first = 0;
426  else if (delta_max > delta_final)
427  _dof_to_damage[node].first = 1.0;
428  else
429  _dof_to_damage[node].first =
430  delta_final * (delta_max - delta_initial) / delta_max / (delta_final - delta_initial);
431 
432  if (_dof_to_damage[node].first < _dof_to_damage[node].second)
433  // Irreversibility
434  _dof_to_damage[node].first = _dof_to_damage[node].second;
435 
436  // Viscous regularization
437  _dof_to_damage[node].first =
438  (_dof_to_damage[node].first + _viscosity * _dof_to_damage[node].second / _dt) /
439  (_viscosity / _dt + 1.0);
440 }
const Real & _dt
Current delta t... or timestep size.
std::unordered_map< const DofObject *, ADReal > _dof_to_delta_initial
const Real _viscosity
Viscosity for damage model.
std::unordered_map< const DofObject *, std::pair< ADReal, Real > > _dof_to_damage
std::unordered_map< const DofObject *, ADReal > _dof_to_delta_max
std::unordered_map< const DofObject *, ADReal > _dof_to_delta_final

◆ computeEffectiveDisplacementJump()

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

Definition at line 401 of file BilinearMixedModeCohesiveZoneModel.C.

Referenced by computeBilinearMixedModeTraction().

402 {
403  using std::sqrt;
404 
405  const auto interface_displacement_jump =
407 
408  const auto delta_normal_pos =
409  MathUtils::regularizedHeavyside(interface_displacement_jump(0), _regularization_alpha) *
410  interface_displacement_jump(0);
411 
412  _dof_to_delta_max[node] = sqrt(Utility::pow<2>(interface_displacement_jump(1)) +
413  Utility::pow<2>(interface_displacement_jump(2)) +
414  Utility::pow<2>(delta_normal_pos) + _epsilon_tolerance);
415 }
std::unordered_map< const DofObject *, ADRealVectorValue > _dof_to_interface_displacement_jump
Map from degree of freedom to local displacement jump.
const Real _regularization_alpha
Parameter for the regularization of the Macaulay bracket.
std::unordered_map< const DofObject *, ADReal > _dof_to_delta_max
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)
CTSub CT_OPERATOR_BINARY CTMul CTCompareLess CTCompareGreater CTCompareEqual _arg template * sqrt(_arg)) *_arg.template D< dtag >()) CT_SIMPLE_UNARY_FUNCTION(tanh
T regularizedHeavyside(const T &x, Real smoothing_length)

◆ computeFandR()

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

Reimplemented from PenaltySimpleCohesiveZoneModel.

Definition at line 233 of file BilinearMixedModeCohesiveZoneModel.C.

234 {
235  using std::isfinite;
236 
237  // First call does not have maps available
238  const bool return_boolean = _dof_to_F.find(node) == _dof_to_F.end();
239  if (return_boolean)
240  return;
241 
242  const auto normalized_F = normalizeQuantity(_dof_to_F, node);
243  const auto normalized_F_neighbor = normalizeQuantity(_dof_to_F_neighbor, node);
244 
245  // This 'averaging' assumption below can probably be improved upon.
246  _dof_to_interface_F[node] = 0.5 * (normalized_F + normalized_F_neighbor);
247 
248  for (const auto i : make_range(3))
249  for (const auto j : make_range(3))
250  if (!isfinite(MetaPhysicL::raw_value(normalized_F(i, j))))
251  throw MooseException(
252  "The deformation gradient on the secondary surface is not finite in "
253  "PenaltySimpleCohesiveZoneModel. MOOSE needs to cut the time step size.");
254 
255  const auto dof_to_interface_F_node = libmesh_map_find(_dof_to_interface_F, node);
256 
257  const ADFactorizedRankTwoTensor C = dof_to_interface_F_node.transpose() * dof_to_interface_F_node;
258  const auto Uinv = MathUtils::sqrt(C).inverse().get();
259  _dof_to_interface_R[node] = dof_to_interface_F_node * Uinv;
260 
261  // Transform interface jump according to two rotation matrices
262  const auto global_interface_displacement = _dof_to_interface_displacement_jump[node];
264  (_dof_to_interface_R[node] * _dof_to_rotation_matrix[node]).transpose() *
265  global_interface_displacement;
266 }
std::unordered_map< const DofObject *, ADRealVectorValue > _dof_to_interface_displacement_jump
Map from degree of freedom to local displacement jump.
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_rotation_matrix
*** Kinematics/displacement jump quantities *** Map from degree of freedom to rotation matrix ...
auto raw_value(const Eigen::Map< T > &in)
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_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.
T normalizeQuantity(const std::unordered_map< const DofObject *, T > &map, const Node *const node)
Normalize mortar quantities (remove mortar integral scaling)
IntRange< T > make_range(T beg, T end)
static const std::complex< double > j(0, 1)
Complex number "j" (also known as "i")
static const std::string C
Definition: NS.h:168

◆ computeFinalDisplacementJump()

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

Definition at line 361 of file BilinearMixedModeCohesiveZoneModel.C.

Referenced by computeBilinearMixedModeTraction().

362 {
363  using std::sqrt, std::pow;
364 
365  const auto interface_displacement_jump =
367 
368  const auto mixity_ratio = libmesh_map_find(_dof_to_mode_mixity_ratio, node);
369 
370  const auto normalized_GI_c = normalizeQuantity(_dof_to_GI_c, node);
371  const auto normalized_GII_c = normalizeQuantity(_dof_to_GII_c, node);
372 
373  _dof_to_delta_final[node] =
374  sqrt(2.0) * 2.0 * normalized_GII_c / normalizeQuantity(_dof_to_shear_strength, node);
375 
376  if (interface_displacement_jump(0) > _epsilon_tolerance)
377  {
379  {
380  _dof_to_delta_final[node] =
381  2.0 / _penalty_stiffness_czm / libmesh_map_find(_dof_to_delta_initial, node) *
382  (normalized_GI_c +
383  (normalized_GII_c - normalized_GI_c) *
384  pow(mixity_ratio * mixity_ratio / (1 + mixity_ratio * mixity_ratio),
386  }
388  {
389  const auto Gc_mixed =
390  pow(1.0 / normalized_GI_c, _power_law_parameter) +
391  pow(mixity_ratio * mixity_ratio / normalized_GII_c, _power_law_parameter);
392  _dof_to_delta_final[node] = (2.0 + 2.0 * mixity_ratio * mixity_ratio) /
394  libmesh_map_find(_dof_to_delta_initial, node) *
395  pow(Gc_mixed, -1.0 / _power_law_parameter);
396  }
397  }
398 }
std::unordered_map< const DofObject *, ADRealVectorValue > _dof_to_interface_displacement_jump
Map from degree of freedom to local displacement jump.
std::unordered_map< const DofObject *, ADReal > _dof_to_delta_initial
std::unordered_map< const DofObject *, ADReal > _dof_to_GII_c
const Real _epsilon_tolerance
Tolerance to avoid NaN/Inf in automatic differentiation operations.
std::unordered_map< const DofObject *, ADReal > _dof_to_delta_final
T normalizeQuantity(const std::unordered_map< const DofObject *, T > &map, const Node *const node)
Normalize mortar quantities (remove mortar integral scaling)
ExpressionBuilder::EBTerm pow(const ExpressionBuilder::EBTerm &left, T exponent)
std::unordered_map< const DofObject *, ADReal > _dof_to_shear_strength
CTSub CT_OPERATOR_BINARY CTMul CTCompareLess CTCompareGreater CTCompareEqual _arg template * sqrt(_arg)) *_arg.template D< dtag >()) CT_SIMPLE_UNARY_FUNCTION(tanh
std::unordered_map< const DofObject *, ADReal > _dof_to_GI_c
enum BilinearMixedModeCohesiveZoneModel::MixedModeCriterion _mix_mode_criterion
const Real _power_law_parameter
Power law parameter for bilinear traction model.
MooseUnits pow(const MooseUnits &, int)
std::unordered_map< const DofObject *, ADReal > _dof_to_mode_mixity_ratio
Map from degree of freedom to mode mixity ratio (AD needed?)

◆ computeGlobalTraction()

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

Compute global traction for mortar application.

Reimplemented from PenaltySimpleCohesiveZoneModel.

Definition at line 300 of file BilinearMixedModeCohesiveZoneModel.C.

301 {
302  // First call does not have maps available
303  const bool return_boolean = _dof_to_czm_traction.find(node) == _dof_to_czm_traction.end();
304  if (return_boolean)
305  return;
306 
307  const auto local_traction_vector = libmesh_map_find(_dof_to_czm_traction, node);
308  const auto rotation_matrix = libmesh_map_find(_dof_to_rotation_matrix, node);
309 
310  _dof_to_czm_traction[node] = rotation_matrix * local_traction_vector;
311 }
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_czm_traction
Total Lagrangian stress to be applied on CZM interface.

◆ computeModeMixity()

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

Definition at line 314 of file BilinearMixedModeCohesiveZoneModel.C.

Referenced by computeBilinearMixedModeTraction().

315 {
316  using std::sqrt;
317 
318  const auto interface_displacement_jump =
320 
321  if (interface_displacement_jump(0) > _epsilon_tolerance)
322  {
323  const auto delta_s =
324  sqrt(interface_displacement_jump(1) * interface_displacement_jump(1) +
325  interface_displacement_jump(2) * interface_displacement_jump(2) + _epsilon_tolerance);
326 
327  _dof_to_mode_mixity_ratio[node] = delta_s / interface_displacement_jump(0);
328  }
329  else
330  _dof_to_mode_mixity_ratio[node] = 0;
331 }
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)
CTSub CT_OPERATOR_BINARY CTMul CTCompareLess CTCompareGreater CTCompareEqual _arg template * sqrt(_arg)) *_arg.template D< dtag >()) CT_SIMPLE_UNARY_FUNCTION(tanh
std::unordered_map< const DofObject *, ADReal > _dof_to_mode_mixity_ratio
Map from degree of freedom to mode mixity ratio (AD needed?)

◆ computeQpIProperties()

void BilinearMixedModeCohesiveZoneModel::computeQpIProperties ( )
overrideprotectedvirtual

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

Reimplemented from WeightedVelocitiesUserObject.

Definition at line 137 of file BilinearMixedModeCohesiveZoneModel.C.

138 {
140 
141  // Get the _dof_to_weighted_gap map
142  const auto * const dof = static_cast<const DofObject *>(_lower_secondary_elem->node_ptr(_i));
143 
144  // TODO: Probably better to interpolate the deformation gradients.
145  _dof_to_F[dof] += (*_test)[_i][_qp] * _F_interpolation;
146  _dof_to_F_neighbor[dof] += (*_test)[_i][_qp] * _F_neighbor_interpolation;
149 
150  _dof_to_GI_c[dof] += (*_test)[_i][_qp] * _GI_c_interpolation;
151  _dof_to_GII_c[dof] += (*_test)[_i][_qp] * _GII_c_interpolation;
152 }
std::unordered_map< const DofObject *, ADReal > _dof_to_normal_strength
ADReal _normal_strength_interpolation
Interpolated value of normal_strength.
unsigned int _qp
Quadrature point index for the mortar segments.
ADReal _shear_strength_interpolation
Interpolated value of shear_strength.
std::unordered_map< const DofObject *, ADReal > _dof_to_GII_c
Elem const *const & _lower_secondary_elem
ADRankTwoTensor _F_interpolation
Deformation gradient for interpolation.
std::unordered_map< const DofObject *, ADRankTwoTensor > _dof_to_F_neighbor
Map from degree of freedom to neighbor, interpolated deformation gradient tensor. ...
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 *, ADReal > _dof_to_shear_strength
ADReal _GI_c_interpolation
Interpolated value of fracture paramter mode I.
std::unordered_map< const DofObject *, ADReal > _dof_to_GI_c
unsigned int _i
Test function index.
virtual void computeQpIProperties() override
Computes properties that are functions both of _qp and _i, for example the weighted gap...
ADReal _GII_c_interpolation
Interpolated value of fracture paramter mode II.

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

Definition at line 112 of file BilinearMixedModeCohesiveZoneModel.C.

113 {
115 
116  // Called after computeQpProperties() within the same algorithmic step.
117 
118  // Compute F and R.
119  const auto F = (ADRankTwoTensor::Identity() +
121  (*_grad_disp[0])[_qp], (*_grad_disp[1])[_qp], (*_grad_disp[2])[_qp]));
122  const auto F_neighbor = (ADRankTwoTensor::Identity() +
124  (*_grad_disp_neighbor[1])[_qp],
125  (*_grad_disp_neighbor[2])[_qp]));
126 
127  // TODO in follow-on PRs: Trim interior node variable derivatives
129  _F_neighbor_interpolation = F_neighbor * (_JxW_msm[_qp] * _coord[_qp]);
134 }
ADReal _normal_strength_interpolation
Interpolated value of normal_strength.
unsigned int _qp
Quadrature point index for the mortar segments.
const MaterialProperty< Real > & _GI_c
Fracture parameter mode I.
std::vector< const GenericVariableGradient< true > * > _grad_disp
Coupled displacement gradients.
const std::vector< Real > & _JxW_msm
const MaterialProperty< Real > & _GII_c
Fracture parameter mode II.
ADReal _shear_strength_interpolation
Interpolated value of shear_strength.
static RankTwoTensorTempl Identity()
static const std::string F
Definition: NS.h:165
static RankTwoTensorTempl initializeFromRows(const libMesh::TypeVector< ADReal > &row0, const libMesh::TypeVector< ADReal > &row1, const libMesh::TypeVector< ADReal > &row2)
ADRankTwoTensor _F_interpolation
Deformation gradient for interpolation.
const MaterialProperty< Real > & _shear_strength
The shear strength material property.
const MaterialProperty< Real > & _normal_strength
The normal strength material property.
ADRankTwoTensor _F_neighbor_interpolation
Deformation gradient for interpolation of the neighbor projection.
virtual void computeQpProperties() override
Computes properties that are functions only of the current quadrature point (_qp), e.g.
const MooseArray< Real > & _coord
Member for handling change of coordinate systems (xyz, rz, spherical)
ADReal _GI_c_interpolation
Interpolated value of fracture paramter mode I.
std::vector< const GenericVariableGradient< true > * > _grad_disp_neighbor
Coupled displacement and neighbor displacement gradient.
ADReal _GII_c_interpolation
Interpolated value of fracture paramter mode II.

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

Definition at line 44 of file BilinearMixedModeCohesiveZoneModel.h.

Referenced by finalize().

44 { return false; }

◆ contactPressure()

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

Implements WeightedGapUserObject.

Definition at line 113 of file PenaltyWeightedGapUserObject.C.

114 {
115  return _contact_pressure;
116 }
ADVariableValue _contact_pressure
The contact pressure on the mortar segument quadrature points.

◆ contactTangentialPressureDirOne() [1/2]

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

Implements WeightedVelocitiesUserObject.

Definition at line 255 of file PenaltySimpleCohesiveZoneModel.C.

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

◆ contactTangentialPressureDirOne() [2/2]

virtual const ADVariableValue& WeightedVelocitiesUserObject::contactTangentialPressureDirOne ( ) const
pure virtualinherited
Returns
The contact force at quadrature points on the mortar segment

Implemented in PenaltyFrictionUserObject, LMWeightedVelocitiesUserObject, and PenaltySimpleCohesiveZoneModel.

Referenced by TangentialMortarMechanicalContact::computeQpResidual().

◆ contactTangentialPressureDirTwo() [1/2]

const ADVariableValue & PenaltySimpleCohesiveZoneModel::contactTangentialPressureDirTwo ( ) const
overridevirtualinherited

Implements WeightedVelocitiesUserObject.

Definition at line 261 of file PenaltySimpleCohesiveZoneModel.C.

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

◆ contactTangentialPressureDirTwo() [2/2]

virtual const ADVariableValue& WeightedVelocitiesUserObject::contactTangentialPressureDirTwo ( ) const
pure virtualinherited

◆ czmGlobalTraction()

const ADVariableValue & BilinearMixedModeCohesiveZoneModel::czmGlobalTraction ( unsigned int  i) const
virtual

Definition at line 503 of file BilinearMixedModeCohesiveZoneModel.C.

Referenced by MortarGenericTraction::computeQpResidual().

504 {
505  mooseAssert(i <= 3,
506  "Internal error in czmGlobalTraction. Index exceeds the traction vector size.");
507 
508  return _czm_interpolated_traction[i];
509 }
std::vector< ADVariableValue > _czm_interpolated_traction
The global traction.

◆ dofToWeightedGap()

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

Get the degree of freedom to weighted gap information.

Definition at line 200 of file WeightedGapUserObject.h.

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

201 {
202  return _dof_to_weighted_gap;
203 }
std::unordered_map< const DofObject *, std::pair< ADReal, Real > > _dof_to_weighted_gap
A map from node to weighted gap and normalization (if requested)

◆ dofToWeightedVelocities()

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

Get the degree of freedom to weighted velocities information.

Definition at line 107 of file WeightedVelocitiesUserObject.h.

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

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

◆ execute()

void WeightedGapUserObject::execute ( )
overridevirtualinherited

Implements MortarUserObject.

Definition at line 161 of file WeightedGapUserObject.C.

162 {
163  for (_qp = 0; _qp < _qrule_msm->n_points(); _qp++)
164  {
166  for (_i = 0; _i < _test->size(); ++_i)
168  }
169 }
unsigned int _qp
Quadrature point index for the mortar segments.
const libMesh::QBase *const & _qrule_msm
const VariableTestValue * _test
A pointer to the test function associated with the weighted gap.
virtual void computeQpProperties()
Computes properties that are functions only of the current quadrature point (_qp), e.g.
unsigned int n_points() const
unsigned int _i
Test function index.
virtual void computeQpIProperties()
Computes properties that are functions both of _qp and _i, for example the weighted gap...

◆ finalize()

void BilinearMixedModeCohesiveZoneModel::finalize ( )
overridevirtual

Reimplemented from PenaltySimpleCohesiveZoneModel.

Definition at line 474 of file BilinearMixedModeCohesiveZoneModel.C.

475 {
477 
478  const bool send_data_back = !constrainedByOwner();
479 
481  _dof_to_F, _subproblem.mesh(), _nodal, _communicator, send_data_back);
482 
485 
488 
491 
493  _dof_to_GI_c, _subproblem.mesh(), _nodal, _communicator, send_data_back);
494 
496  _dof_to_GII_c, _subproblem.mesh(), _nodal, _communicator, send_data_back);
497 
500 }
virtual MooseMesh & mesh()=0
std::unordered_map< const DofObject *, ADReal > _dof_to_normal_strength
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)
std::unordered_map< const DofObject *, ADRealVectorValue > _dof_to_weighted_displacements
A map from node to weighted displacements.
const bool _nodal
Whether the dof objects are nodal; if they&#39;re not, then they&#39;re elemental.
const Parallel::Communicator & _communicator
std::unordered_map< const DofObject *, ADReal > _dof_to_GII_c
SubProblem & _subproblem
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 *, ADReal > _dof_to_shear_strength
std::unordered_map< const DofObject *, ADReal > _dof_to_GI_c
void communicateR2T(std::unordered_map< const DofObject *, ADRankTwoTensor > &dof_map_adr2t, const MooseMesh &mesh, const bool nodal, const Parallel::Communicator &communicator, const bool send_data_back)
This function is used to communicate velocities across processes.

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

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

123  {
124  const auto it = map.find(key);
125  if (it == map.end())
126  return default_value;
127  return it->second;
128  }

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

79  {
80  mooseError("Not available in base class.");
81  }
void mooseError(Args &&... args) const

◆ getActiveSetState()

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

Definition at line 45 of file PenaltyWeightedGapUserObject.h.

Referenced by MortarUserObjectAux::MortarUserObjectAux().

46  {
47  return _dof_to_normal_pressure.find(_subproblem.mesh().nodePtr(node->id()))->second > 0;
48  }
virtual MooseMesh & mesh()=0
std::unordered_map< const DofObject *, ADReal > _dof_to_normal_pressure
Map from degree of freedom to normal pressure for reporting.
SubProblem & _subproblem
virtual const Node * nodePtr(const dof_id_type i) const

◆ getCohesiveDamage()

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

Definition at line 523 of file BilinearMixedModeCohesiveZoneModel.C.

Referenced by CohesiveZoneMortarUserObjectAux::CohesiveZoneMortarUserObjectAux().

524 {
525  const auto it = _dof_to_damage.find(_subproblem.mesh().nodePtr(node->id()));
526 
527  if (it != _dof_to_damage.end())
528  return MetaPhysicL::raw_value(it->second.first);
529  else
530  return 0.0;
531 }
virtual MooseMesh & mesh()=0
auto raw_value(const Eigen::Map< T > &in)
std::unordered_map< const DofObject *, std::pair< ADReal, Real > > _dof_to_damage
SubProblem & _subproblem
virtual const Node * nodePtr(const dof_id_type i) const

◆ getDeltaTangentialLagrangeMultiplier()

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

Reimplemented in PenaltyFrictionUserObject.

Definition at line 41 of file PenaltyWeightedGapUserObject.h.

42  {
43  return 0.0;
44  };

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

68  {
69  mooseError("Not available in base class.");
70  }
void mooseError(Args &&... args) const

◆ getLocalDisplacementNormal()

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

Definition at line 534 of file BilinearMixedModeCohesiveZoneModel.C.

Referenced by CohesiveZoneMortarUserObjectAux::CohesiveZoneMortarUserObjectAux().

535 {
536  const auto it = _dof_to_interface_displacement_jump.find(_subproblem.mesh().nodePtr(node->id()));
537  const auto it2 = _dof_to_weighted_gap.find(_subproblem.mesh().nodePtr(node->id()));
538 
539  if (it != _dof_to_interface_displacement_jump.end() && it2 != _dof_to_weighted_gap.end())
540  return MetaPhysicL::raw_value(it->second(0) / it2->second.second);
541  else
542  return 0.0;
543 }
std::unordered_map< const DofObject *, ADRealVectorValue > _dof_to_interface_displacement_jump
Map from degree of freedom to local displacement jump.
virtual MooseMesh & mesh()=0
auto raw_value(const Eigen::Map< T > &in)
SubProblem & _subproblem
virtual const Node * nodePtr(const dof_id_type i) const
std::unordered_map< const DofObject *, std::pair< ADReal, Real > > _dof_to_weighted_gap
A map from node to weighted gap and normalization (if requested)

◆ getLocalDisplacementTangential()

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

Definition at line 546 of file BilinearMixedModeCohesiveZoneModel.C.

Referenced by CohesiveZoneMortarUserObjectAux::CohesiveZoneMortarUserObjectAux().

547 {
548  const auto it = _dof_to_interface_displacement_jump.find(_subproblem.mesh().nodePtr(node->id()));
549  const auto it2 = _dof_to_weighted_gap.find(_subproblem.mesh().nodePtr(node->id()));
550 
551  if (it != _dof_to_interface_displacement_jump.end() && it2 != _dof_to_weighted_gap.end())
552  return MetaPhysicL::raw_value(it->second(1) / it2->second.second);
553  else
554  return 0.0;
555 }
std::unordered_map< const DofObject *, ADRealVectorValue > _dof_to_interface_displacement_jump
Map from degree of freedom to local displacement jump.
virtual MooseMesh & mesh()=0
auto raw_value(const Eigen::Map< T > &in)
SubProblem & _subproblem
virtual const Node * nodePtr(const dof_id_type i) const
std::unordered_map< const DofObject *, std::pair< ADReal, Real > > _dof_to_weighted_gap
A map from node to weighted gap and normalization (if requested)

◆ getModeMixityRatio()

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

Definition at line 512 of file BilinearMixedModeCohesiveZoneModel.C.

Referenced by CohesiveZoneMortarUserObjectAux::CohesiveZoneMortarUserObjectAux().

513 {
514  const auto it = _dof_to_mode_mixity_ratio.find(_subproblem.mesh().nodePtr(node->id()));
515 
516  if (it != _dof_to_mode_mixity_ratio.end())
517  return MetaPhysicL::raw_value(it->second);
518  else
519  return 0.0;
520 }
virtual MooseMesh & mesh()=0
auto raw_value(const Eigen::Map< T > &in)
SubProblem & _subproblem
virtual const Node * nodePtr(const dof_id_type i) const
std::unordered_map< const DofObject *, ADReal > _dof_to_mode_mixity_ratio
Map from degree of freedom to mode mixity ratio (AD needed?)

◆ getNormalContactPressure()

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

Implements WeightedGapUserObject.

Definition at line 164 of file PenaltyWeightedGapUserObject.C.

165 {
166  const auto it = _dof_to_normal_pressure.find(_subproblem.mesh().nodePtr(node->id()));
167 
168  if (it != _dof_to_normal_pressure.end())
169  return MetaPhysicL::raw_value(it->second);
170  else
171  return 0.0;
172 }
virtual MooseMesh & mesh()=0
std::unordered_map< const DofObject *, ADReal > _dof_to_normal_pressure
Map from degree of freedom to normal pressure for reporting.
auto raw_value(const Eigen::Map< T > &in)
SubProblem & _subproblem
virtual const Node * nodePtr(const dof_id_type i) const

◆ getNormalGap()

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

Definition at line 172 of file WeightedGapUserObject.C.

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

173 {
174  const auto it = _dof_to_weighted_gap.find(_subproblem.mesh().nodePtr(node->id()));
175 
176  // We are returning the physical weighted gap for analysis purposes
177  if (it != _dof_to_weighted_gap.end())
178  return physicalGap(it->second);
179  else
180  return 0.0;
181 }
virtual MooseMesh & mesh()=0
SubProblem & _subproblem
virtual const Node * nodePtr(const dof_id_type i) const
Real physicalGap(const std::pair< ADReal, Real > &gap) const
Compute physical gap from integration gap quantity.
std::unordered_map< const DofObject *, std::pair< ADReal, Real > > _dof_to_weighted_gap
A map from node to weighted gap and normalization (if requested)

◆ getNormalLagrangeMultiplier()

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

Definition at line 175 of file PenaltyWeightedGapUserObject.C.

Referenced by MortarUserObjectAux::MortarUserObjectAux().

176 {
177  const auto it = _dof_to_lagrange_multiplier.find(_subproblem.mesh().nodePtr(node->id()));
178 
179  if (it != _dof_to_lagrange_multiplier.end())
180  return -MetaPhysicL::raw_value(it->second);
181  else
182  return 0.0;
183 }
virtual MooseMesh & mesh()=0
auto raw_value(const Eigen::Map< T > &in)
SubProblem & _subproblem
virtual const Node * nodePtr(const dof_id_type i) const
std::unordered_map< const DofObject *, Real > _dof_to_lagrange_multiplier
Map from degree of freedom to augmented lagrange multiplier.

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

Referenced by MortarUserObjectAux::MortarUserObjectAux().

90  {
91  mooseError("Not available in base class.");
92  }
void mooseError(Args &&... args) const

◆ initialize()

void BilinearMixedModeCohesiveZoneModel::initialize ( )
overridevirtual

Reimplemented from PenaltySimpleCohesiveZoneModel.

Definition at line 181 of file BilinearMixedModeCohesiveZoneModel.C.

182 {
183  // instead we call it explicitly here
185 
186  // Avoid accumulating interpolation over the time step
187  _dof_to_F.clear();
188  _dof_to_F_neighbor.clear();
190  _dof_to_normal_strength.clear();
191  _dof_to_shear_strength.clear();
192  _dof_to_GI_c.clear();
193  _dof_to_GII_c.clear();
194  _dof_to_delta_initial.clear();
195  _dof_to_delta_final.clear();
196  _dof_to_delta_max.clear();
197  _dof_to_czm_traction.clear();
200  _dof_to_interface_F.clear();
201  _dof_to_interface_R.clear();
202 
203  for (auto & map_pr : _dof_to_rotation_matrix)
204  map_pr.second.setToIdentity();
205 }
std::unordered_map< const DofObject *, ADRealVectorValue > _dof_to_interface_displacement_jump
Map from degree of freedom to local displacement jump.
std::unordered_map< const DofObject *, ADReal > _dof_to_normal_strength
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_rotation_matrix
*** Kinematics/displacement jump quantities *** Map from degree of freedom to rotation matrix ...
std::unordered_map< const DofObject *, ADReal > _dof_to_delta_initial
std::unordered_map< const DofObject *, ADRankTwoTensor > _dof_to_interface_R
Map from degree of freedom to interface rotation tensor.
std::unordered_map< const DofObject *, ADReal > _dof_to_delta_max
std::unordered_map< const DofObject *, ADReal > _dof_to_GII_c
std::unordered_map< const DofObject *, ADRealVectorValue > _dof_to_czm_traction
Total Lagrangian stress to be applied on CZM interface.
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 *, ADReal > _dof_to_delta_final
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 *, ADReal > _dof_to_czm_normal_traction
Map from degree of freedom to czm normal traction.
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_mode_mixity_ratio
Map from degree of freedom to mode mixity ratio (AD needed?)

◆ initialSetup()

void WeightedGapUserObject::initialSetup ( )
overridevirtualinherited

Reimplemented from MortarUserObject.

Definition at line 60 of file WeightedGapUserObject.C.

61 {
63  _test = &test();
64 }
virtual const VariableTestValue & test() const =0
const VariableTestValue * _test
A pointer to the test function associated with the weighted gap.
virtual void initialSetup()

◆ isAugmentedLagrangianConverged()

bool PenaltyWeightedGapUserObject::isAugmentedLagrangianConverged ( )
overridevirtualinherited

Implements AugmentedLagrangeInterface.

Reimplemented in PenaltyFrictionUserObject.

Definition at line 226 of file PenaltyWeightedGapUserObject.C.

Referenced by PenaltyFrictionUserObject::isAugmentedLagrangianConverged().

227 {
228  using std::max, std::min, std::abs;
229 
230  Real max_positive_gap = 0.0;
231  Real min_negative_gap = 0.0;
232 
233  // Get maximum gap to ascertain whether we are converged.
234  for (auto & [dof_object, wgap] : _dof_to_weighted_gap)
235  {
236  const auto gap = physicalGap(wgap);
237  {
238  // Check condition for nodes that are active
239  if (gap < 0 || _dof_to_lagrange_multiplier[dof_object] < 0.0)
240  {
241  max_positive_gap = max(max_positive_gap, gap);
242  min_negative_gap = min(min_negative_gap, gap);
243  }
244  }
245  }
246 
247  // Communicate the extreme gap values in parallel.
248  std::vector<Real> recv;
249  if (this->_communicator.rank() == 0)
250  recv.resize(this->_communicator.size());
251  this->_communicator.gather(
252  0, -min_negative_gap > max_positive_gap ? min_negative_gap : max_positive_gap, recv);
253 
254  if (this->_communicator.rank() == 0)
255  {
256  min_negative_gap = *std::min_element(recv.begin(), recv.end());
257  max_positive_gap = *std::max_element(recv.begin(), recv.end());
258 
259  // report the gap value with the largest magnitude
260  const Real reported_gap =
261  -min_negative_gap > max_positive_gap ? min_negative_gap : max_positive_gap;
262  if (abs(reported_gap) > _penetration_tolerance)
263  {
264  mooseInfoRepeated("Penetration tolerance fail max_gap = ",
265  reported_gap,
266  " (gap_tol=",
268  "). Iteration number is: ",
270  ".");
271  return false;
272  }
273  else
274  mooseInfoRepeated("Penetration tolerance success max_gap = ",
275  reported_gap,
276  " (gap_tol=",
278  "). Iteration number is: ",
280  ".");
281  }
282 
283  return true;
284 }
MetaPhysicL::DualNumber< V, D, asd > abs(const MetaPhysicL::DualNumber< V, D, asd > &a)
void gather(const unsigned int root_id, const T &send_data, std::vector< T, A > &recv) const
void mooseInfoRepeated(Args &&... args)
processor_id_type rank() const
const Parallel::Communicator & _communicator
auto max(const L &left, const R &right)
const Real _penetration_tolerance
penetration tolerance for augmented Lagrange contact
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
std::unordered_map< const DofObject *, Real > _dof_to_lagrange_multiplier
Map from degree of freedom to augmented lagrange multiplier.
Real physicalGap(const std::pair< ADReal, Real > &gap) const
Compute physical gap from integration gap quantity.
const unsigned int & _lagrangian_iteration_number
std::unordered_map< const DofObject *, std::pair< ADReal, Real > > _dof_to_weighted_gap
A map from node to weighted gap and normalization (if requested)
auto min(const L &left, const R &right)

◆ normalizeQuantity()

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

Normalize mortar quantities (remove mortar integral scaling)

Definition at line 156 of file BilinearMixedModeCohesiveZoneModel.C.

Referenced by computeBilinearMixedModeTraction(), computeCriticalDisplacementJump(), computeEffectiveDisplacementJump(), computeFandR(), computeFinalDisplacementJump(), and computeModeMixity().

158 {
159  return libmesh_map_find(map, node) / _dof_to_weighted_gap[node].second;
160 }
std::unordered_map< const DofObject *, std::pair< ADReal, Real > > _dof_to_weighted_gap
A map from node to weighted gap and normalization (if requested)

◆ physicalGap()

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

Compute physical gap from integration gap quantity.

Definition at line 54 of file WeightedGapUserObject.h.

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

55  {
56  return MetaPhysicL::raw_value(gap.first) / gap.second;
57  }
auto raw_value(const Eigen::Map< T > &in)

◆ prepareJumpKinematicQuantities()

void BilinearMixedModeCohesiveZoneModel::prepareJumpKinematicQuantities ( )
overrideprotectedvirtual

Reimplemented from PenaltySimpleCohesiveZoneModel.

Definition at line 208 of file BilinearMixedModeCohesiveZoneModel.C.

209 {
210  // Compute rotation matrix from secondary surface
211  // Rotation matrices and local interface displacement jump.
212  for (const auto i : make_range(_test->size()))
213  {
214  const Node * const node = _lower_secondary_elem->node_ptr(i);
215 
216  // First call does not have maps available
217  const bool return_boolean = _dof_to_weighted_gap.find(node) == _dof_to_weighted_gap.end();
218  if (return_boolean)
219  return;
220 
221  _dof_to_rotation_matrix[node] = CohesiveZoneModelTools::computeReferenceRotation<true>(
223 
224  // Every time we used quantities from a map we "denormalize" it from the mortar integral.
225  // See normalizing member functions.
228  (libmesh_map_find(_dof_to_weighted_displacements, node));
229  }
230 }
std::unordered_map< const DofObject *, ADRealVectorValue > _dof_to_interface_displacement_jump
Map from degree of freedom to local displacement jump.
virtual MooseMesh & mesh()=0
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_weighted_displacements
A map from node to weighted displacements.
SubProblem & _subproblem
Elem const *const & _lower_secondary_elem
virtual unsigned int dimension() const
std::vector< Point > _normals
const VariableTestValue * _test
A pointer to the test function associated with the weighted gap.
IntRange< T > make_range(T beg, T end)
std::unordered_map< const DofObject *, std::pair< ADReal, Real > > _dof_to_weighted_gap
A map from node to weighted gap and normalization (if requested)

◆ reinit()

void BilinearMixedModeCohesiveZoneModel::reinit ( )
overridevirtual

Reimplemented from PenaltySimpleCohesiveZoneModel.

Definition at line 443 of file BilinearMixedModeCohesiveZoneModel.C.

444 {
445  // Normal contact pressure with penalty
447 
448  for (const auto i : index_range(_czm_interpolated_traction))
449  {
451 
452  for (const auto qp : make_range(_qrule_msm->n_points()))
453  _czm_interpolated_traction[i][qp] = 0.0;
454  }
455 
456  // iterate over nodes
457  for (const auto i : make_range(_test->size()))
458  {
459  // current node
460  const Node * const node = _lower_secondary_elem->node_ptr(i);
461 
462  // End of CZM bilinear computations
463  if (_dof_to_czm_traction.find(node) == _dof_to_czm_traction.end())
464  return;
465 
466  const auto & test_i = (*_test)[i];
467  for (const auto qp : make_range(_qrule_msm->n_points()))
468  for (const auto idx : index_range(_czm_interpolated_traction))
469  _czm_interpolated_traction[idx][qp] += test_i[qp] * _dof_to_czm_traction[node](idx);
470  }
471 }
const libMesh::QBase *const & _qrule_msm
Elem const *const & _lower_secondary_elem
std::unordered_map< const DofObject *, ADRealVectorValue > _dof_to_czm_traction
Total Lagrangian stress to be applied on CZM interface.
std::vector< ADVariableValue > _czm_interpolated_traction
The global traction.
const VariableTestValue * _test
A pointer to the test function associated with the weighted gap.
unsigned int n_points() const
IntRange< T > make_range(T beg, T end)
auto index_range(const T &sizable)

◆ selfFinalize()

void PenaltyWeightedGapUserObject::selfFinalize ( )
protectedinherited

Definition at line 133 of file PenaltyWeightedGapUserObject.C.

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

134 {
135  using std::min;
136 
137  // compute new normal pressure for each node
138  for (const auto & [dof_object, wgap] : _dof_to_weighted_gap)
139  {
140  auto penalty = findValue(_dof_to_local_penalty, dof_object, _penalty);
141 
142  // If _use_physical_gap is true we "normalize" the penalty parameter with the surrounding area.
143  auto gap = _use_physical_gap ? adPhysicalGap(wgap) / wgap.second : wgap.first;
144 
145  const auto lagrange_multiplier =
147 
148  // keep the negative normal pressure (compressive per convention here)
149  auto normal_pressure = min(0.0, penalty * gap + lagrange_multiplier);
150 
151  // we switch conventins here and consider positive normal pressure as compressive
152  _dof_to_normal_pressure[dof_object] = -normal_pressure;
153  }
154 }
ADReal adPhysicalGap(const std::pair< ADReal, Real > &gap) const
std::unordered_map< const DofObject *, ADReal > _dof_to_normal_pressure
Map from degree of freedom to normal pressure for reporting.
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&#39;t exist.
const Real _penalty
The penalty factor.
AugmentedLagrangianContactProblemInterface *const _augmented_lagrange_problem
augmented Lagrange problem and iteration number
const bool _use_physical_gap
Use scaled or physical gap.
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_lagrange_multiplier
Map from degree of freedom to augmented lagrange multiplier.
std::unordered_map< const DofObject *, std::pair< ADReal, Real > > _dof_to_weighted_gap
A map from node to weighted gap and normalization (if requested)
auto min(const L &left, const R &right)

◆ selfInitialize() [1/2]

void WeightedVelocitiesUserObject::selfInitialize ( )
protectedinherited

Definition at line 131 of file WeightedVelocitiesUserObject.C.

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

132 {
135 }
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::unordered_map< const DofObject *, std::array< ADReal, 2 > > _dof_to_weighted_tangential_velocity
A map from node to two weighted tangential velocities.

◆ selfInitialize() [2/2]

void PenaltyWeightedGapUserObject::selfInitialize ( )
protectedinherited

Definition at line 119 of file PenaltyWeightedGapUserObject.C.

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

120 {
121  for (auto & dof_pn : _dof_to_normal_pressure)
122  dof_pn.second = 0.0;
123 }
std::unordered_map< const DofObject *, ADReal > _dof_to_normal_pressure
Map from degree of freedom to normal pressure for reporting.

◆ selfTimestepSetup() [1/2]

void WeightedVelocitiesUserObject::selfTimestepSetup ( )
inlineprotectedinherited

◆ selfTimestepSetup() [2/2]

void PenaltyWeightedGapUserObject::selfTimestepSetup ( )
protectedinherited

Definition at line 202 of file PenaltyWeightedGapUserObject.C.

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

203 {
204  // Let's not clear the LMs for improved performance in
205  // nonlinear problems
206 
207  // reset penalty
208  for (auto & dof_lp : _dof_to_local_penalty)
209  dof_lp.second = _penalty;
210 
211  // clear previous gap
212  for (auto & dof_pg : _dof_to_previous_gap)
213  dof_pg.second = 0.0;
214 
215  // save old timestep
216  _dt_old = _dt;
217 }
const Real _penalty
The penalty factor.
const Real & _dt
Current delta t... or timestep size.
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.

◆ test() [1/2]

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

Reimplemented from PenaltyWeightedGapUserObject.

Definition at line 75 of file PenaltySimpleCohesiveZoneModel.C.

76 {
78 }
const MooseVariable *const _aux_lm_var
The auxiliary Lagrange multiplier variable (used together whith the Petrov-Galerkin approach) ...
virtual const FieldVariablePhiValue & phiLower() const override
const MooseVariable *const _disp_x_var
The x displacement variable.

◆ test() [2/2]

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

Implements WeightedGapUserObject.

Reimplemented in PenaltyFrictionUserObject, and PenaltySimpleCohesiveZoneModel.

Definition at line 107 of file PenaltyWeightedGapUserObject.C.

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

◆ timestepSetup()

void BilinearMixedModeCohesiveZoneModel::timestepSetup ( )
overridevirtual

Reimplemented from PenaltySimpleCohesiveZoneModel.

Definition at line 163 of file BilinearMixedModeCohesiveZoneModel.C.

164 {
165  // instead we call it explicitly here
167 
168  // save off tangential traction from the last timestep
169  for (auto & map_pr : _dof_to_damage)
170  {
171  auto & [damage, old_damage] = map_pr.second;
172  old_damage = {MetaPhysicL::raw_value(damage)};
173  damage = {0.0};
174  }
175 
176  for (auto & [dof_object, delta_tangential_lm] : _dof_to_frictional_lagrange_multipliers)
177  delta_tangential_lm.setZero();
178 }
std::unordered_map< const DofObject *, TwoVector > _dof_to_frictional_lagrange_multipliers
Map from degree of freedom to augmented lagrange multiplier.
auto raw_value(const Eigen::Map< T > &in)
std::unordered_map< const DofObject *, std::pair< ADReal, Real > > _dof_to_damage

◆ updateAugmentedLagrangianMultipliers()

void PenaltyWeightedGapUserObject::updateAugmentedLagrangianMultipliers ( )
overridevirtualinherited

Implements AugmentedLagrangeInterface.

Reimplemented in PenaltyFrictionUserObject.

Definition at line 299 of file PenaltyWeightedGapUserObject.C.

Referenced by PenaltyFrictionUserObject::updateAugmentedLagrangianMultipliers().

300 {
301  using std::abs;
302 
303  for (const auto & [dof_object, wgap] : _dof_to_weighted_gap)
304  {
305  auto & penalty = _dof_to_local_penalty[dof_object];
306  if (penalty == 0.0)
307  penalty = _penalty;
308 
309  const auto gap = getNormalGap(static_cast<const Node *>(dof_object));
310  auto & lagrange_multiplier = _dof_to_lagrange_multiplier[dof_object];
311 
312  const auto possible_normalization =
313  (_use_physical_gap ? libmesh_map_find(_dof_to_weighted_gap, dof_object).second : 1.0);
314 
315  // Update penalty (the factor of 1/4 is suggested in the literature, the limit on AL iteration
316  // caps the penalty increase)
317  // Before we were updating the LM before adapting the penalty factor
318  if (lagrange_multiplier + gap * penalty / possible_normalization <= 0)
319  lagrange_multiplier += gap * penalty / possible_normalization;
320  else
321  lagrange_multiplier = 0.0;
322 
323  const auto previous_gap = _dof_to_previous_gap[dof_object];
324  Real eval_tn = 0;
325 
327  {
328  if (abs(gap) > 0.25 * abs(previous_gap) && abs(gap) > _penetration_tolerance)
329  penalty *= _penalty_multiplier;
330  }
332  bussettaAdaptivePenalty(previous_gap, gap, penalty, eval_tn);
333 
334  // Heuristics to bound the penalty factor
335  if (penalty < _penalty)
336  penalty = _penalty;
337  else if (penalty > _max_penalty_multiplier * _penalty)
338  penalty = _max_penalty_multiplier * _penalty;
339  }
340 }
const Real _max_penalty_multiplier
Maximum multiplier applied to the initial penalty factor in AL.
MetaPhysicL::DualNumber< V, D, asd > abs(const MetaPhysicL::DualNumber< V, D, asd > &a)
const Real _penalty
The penalty factor.
void bussettaAdaptivePenalty(const Real previous_gap, const Real gap, Real &penalty, Real &eval_tn)
Adaptive, local penalty for AL.
const bool _use_physical_gap
Use scaled or physical gap.
virtual Real getNormalGap(const Node *const) const
const Real _penetration_tolerance
penetration tolerance for augmented Lagrange contact
std::unordered_map< const DofObject *, Real > _dof_to_local_penalty
Map from degree of freedom to local penalty value.
enum PenaltyWeightedGapUserObject::AdaptivityNormalPenalty _adaptivity_normal
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
std::unordered_map< const DofObject *, Real > _dof_to_previous_gap
Map from degree of freedom to previous AL iteration gap values.
const Real _penalty_multiplier
penalty growth factor for augmented Lagrange
std::unordered_map< const DofObject *, Real > _dof_to_lagrange_multiplier
Map from degree of freedom to augmented lagrange multiplier.
std::unordered_map< const DofObject *, std::pair< ADReal, Real > > _dof_to_weighted_gap
A map from node to weighted gap and normalization (if requested)

◆ validParams()

InputParameters BilinearMixedModeCohesiveZoneModel::validParams ( )
static

Definition at line 27 of file BilinearMixedModeCohesiveZoneModel.C.

28 {
30 
31  params.addClassDescription("Computes the bilinear mixed mode cohesive zone model.");
32  params.addRequiredCoupledVar("displacements",
33  "The string of displacements suitable for the problem statement");
34 
35  // Input parameters for bilinear mixed mode traction.
36  params.addParam<MaterialPropertyName>("GI_c",
37  "Critical energy release rate in normal direction.");
38  params.addParam<MaterialPropertyName>("GII_c",
39  "Critical energy release rate in shear direction.");
40  params.addParam<MaterialPropertyName>("normal_strength", "Tensile strength in normal direction.");
41  params.addParam<MaterialPropertyName>("shear_strength", "Tensile strength in shear direction.");
42  params.addParam<Real>("power_law_parameter", "The power law parameter.");
43  MooseEnum criterion("POWER_LAW BK", "BK");
44  params.addParam<Real>("viscosity", 0.0, "Viscosity for damage model.");
45  params.addParam<MooseEnum>(
46  "mixed_mode_criterion", criterion, "Option for mixed mode propagation criterion.");
47  params.addParam<bool>(
48  "lag_displacement_jump",
49  false,
50  "Whether to use old displacement jumps to compute the effective displacement jump.");
51  params.addParam<Real>(
52  "regularization_alpha", 1e-10, "Regularization parameter for the Macaulay bracket.");
53  params.addRangeCheckedParam<Real>(
54  "penalty_stiffness", "penalty_stiffness > 0.0", "Penalty stiffness for CZM.");
55  params.addParamNamesToGroup(
56  "GI_c GII_c normal_strength shear_strength power_law_parameter viscosity "
57  "mixed_mode_criterion lag_displacement_jump regularization_alpha "
58  "penalty_stiffness",
59  "Bilinear mixed mode traction");
60  // End of input parameters for bilinear mixed mode traction.
61 
62  // Suppress augmented Lagrange parameters. AL implementation for CZM remains to be done.
63  params.suppressParameter<Real>("max_penalty_multiplier");
64  params.suppressParameter<Real>("penalty_multiplier");
65  params.suppressParameter<Real>("penetration_tolerance");
66  return params;
67 }
void addParam(const std::string &name, const std::initializer_list< typename T::value_type > &value, const std::string &doc_string)
void suppressParameter(const std::string &name)
void addRequiredCoupledVar(const std::string &name, const std::string &doc_string)
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
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)
void addParamNamesToGroup(const std::string &space_delim_names, const std::string group_name)

Member Data Documentation

◆ _3d

bool WeightedVelocitiesUserObject::_3d
protectedinherited

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

Definition at line 103 of file WeightedVelocitiesUserObject.h.

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

◆ _adaptivity_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 115 of file PenaltyWeightedGapUserObject.h.

Referenced by PenaltySimpleCohesiveZoneModel::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 85 of file PenaltyWeightedGapUserObject.h.

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

◆ _coord

const MooseArray<Real>& WeightedGapUserObject::_coord
protectedinherited

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

Definition at line 165 of file WeightedGapUserObject.h.

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

◆ _czm_interpolated_traction

std::vector<ADVariableValue> BilinearMixedModeCohesiveZoneModel::_czm_interpolated_traction
protected

The global traction.

Definition at line 151 of file BilinearMixedModeCohesiveZoneModel.h.

Referenced by BilinearMixedModeCohesiveZoneModel(), czmGlobalTraction(), and 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<const DofObject *, std::pair<TwoVector, TwoVector> > PenaltySimpleCohesiveZoneModel::_dof_to_accumulated_slip
protectedinherited

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

Definition at line 63 of file PenaltySimpleCohesiveZoneModel.h.

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

◆ _dof_to_czm_normal_traction

std::unordered_map<const DofObject *, ADReal> PenaltySimpleCohesiveZoneModel::_dof_to_czm_normal_traction
protectedinherited

Map from degree of freedom to czm normal traction.

Definition at line 70 of file PenaltySimpleCohesiveZoneModel.h.

Referenced by initialize().

◆ _dof_to_czm_traction

std::unordered_map<const DofObject *, ADRealVectorValue> BilinearMixedModeCohesiveZoneModel::_dof_to_czm_traction
protected

Total Lagrangian stress to be applied on CZM interface.

Definition at line 174 of file BilinearMixedModeCohesiveZoneModel.h.

Referenced by computeBilinearMixedModeTraction(), computeGlobalTraction(), initialize(), and reinit().

◆ _dof_to_damage

std::unordered_map<const DofObject *, std::pair<ADReal, Real> > BilinearMixedModeCohesiveZoneModel::_dof_to_damage
protected

◆ _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> BilinearMixedModeCohesiveZoneModel::_dof_to_displacement_jump
protected

Definition at line 176 of file BilinearMixedModeCohesiveZoneModel.h.

◆ _dof_to_F

std::unordered_map<const DofObject *, ADRankTwoTensor> BilinearMixedModeCohesiveZoneModel::_dof_to_F
protected

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

Definition at line 93 of file BilinearMixedModeCohesiveZoneModel.h.

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

◆ _dof_to_F_neighbor

std::unordered_map<const DofObject *, ADRankTwoTensor> BilinearMixedModeCohesiveZoneModel::_dof_to_F_neighbor
protected

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

Definition at line 96 of file BilinearMixedModeCohesiveZoneModel.h.

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

◆ _dof_to_frictional_lagrange_multipliers

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

Map from degree of freedom to augmented lagrange multiplier.

Definition at line 79 of file PenaltySimpleCohesiveZoneModel.h.

Referenced by PenaltySimpleCohesiveZoneModel::reinit(), timestepSetup(), and PenaltySimpleCohesiveZoneModel::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> BilinearMixedModeCohesiveZoneModel::_dof_to_interface_displacement_jump
protected

◆ _dof_to_interface_F

std::unordered_map<const DofObject *, ADRankTwoTensor> BilinearMixedModeCohesiveZoneModel::_dof_to_interface_F
protected

Map from degree of freedom to interface deformation gradient tensor.

Definition at line 99 of file BilinearMixedModeCohesiveZoneModel.h.

Referenced by computeFandR(), and initialize().

◆ _dof_to_interface_R

std::unordered_map<const DofObject *, ADRankTwoTensor> BilinearMixedModeCohesiveZoneModel::_dof_to_interface_R
protected

Map from degree of freedom to interface rotation tensor.

Definition at line 102 of file BilinearMixedModeCohesiveZoneModel.h.

Referenced by computeFandR(), and initialize().

◆ _dof_to_lagrange_multiplier

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

◆ _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> PenaltySimpleCohesiveZoneModel::_dof_to_local_penalty_friction
protectedinherited

Map from degree of freedom to local friction penalty value.

Definition at line 82 of file PenaltySimpleCohesiveZoneModel.h.

Referenced by PenaltySimpleCohesiveZoneModel::reinit(), and PenaltySimpleCohesiveZoneModel::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 105 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> BilinearMixedModeCohesiveZoneModel::_dof_to_rotation_matrix
protected

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

Definition at line 81 of file BilinearMixedModeCohesiveZoneModel.h.

Referenced by computeFandR(), computeGlobalTraction(), initialize(), and 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<TwoVector, TwoVector> > PenaltySimpleCohesiveZoneModel::_dof_to_step_slip
protectedinherited

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

Definition at line 60 of file PenaltySimpleCohesiveZoneModel.h.

◆ _dof_to_tangential_traction

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

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

Definition at line 67 of file PenaltySimpleCohesiveZoneModel.h.

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

◆ _dof_to_weighted_displacements

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

A map from node to weighted displacements.

Definition at line 177 of file WeightedGapUserObject.h.

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

◆ _dof_to_weighted_gap

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

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

Referenced by PenaltyWeightedGapUserObject::selfTimestepSetup().

◆ _epsilon_tolerance

const Real PenaltySimpleCohesiveZoneModel::_epsilon_tolerance
protectedinherited

Tolerance to avoid NaN/Inf in automatic differentiation operations.

Definition at line 85 of file PenaltySimpleCohesiveZoneModel.h.

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

◆ _F_interpolation

ADRankTwoTensor BilinearMixedModeCohesiveZoneModel::_F_interpolation
protected

Deformation gradient for interpolation.

Definition at line 87 of file BilinearMixedModeCohesiveZoneModel.h.

Referenced by computeQpIProperties(), and computeQpProperties().

◆ _F_neighbor_interpolation

ADRankTwoTensor BilinearMixedModeCohesiveZoneModel::_F_neighbor_interpolation
protected

Deformation gradient for interpolation of the neighbor projection.

Definition at line 90 of file BilinearMixedModeCohesiveZoneModel.h.

Referenced by computeQpIProperties(), and computeQpProperties().

◆ _fe_problem

FEProblemBase& WeightedGapUserObject::_fe_problem
protectedinherited

The base finite element problem.

Definition at line 131 of file WeightedGapUserObject.h.

◆ _friction_coefficient

const Real PenaltySimpleCohesiveZoneModel::_friction_coefficient
protectedinherited

The friction coefficient.

Definition at line 57 of file PenaltySimpleCohesiveZoneModel.h.

Referenced by PenaltySimpleCohesiveZoneModel::reinit().

◆ _frictional_contact_traction_one

ADVariableValue PenaltySimpleCohesiveZoneModel::_frictional_contact_traction_one
protectedinherited

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

Definition at line 73 of file PenaltySimpleCohesiveZoneModel.h.

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

◆ _frictional_contact_traction_two

ADVariableValue PenaltySimpleCohesiveZoneModel::_frictional_contact_traction_two
protectedinherited

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

Definition at line 76 of file PenaltySimpleCohesiveZoneModel.h.

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

◆ _GI_c

const MaterialProperty<Real>& BilinearMixedModeCohesiveZoneModel::_GI_c
protected

Fracture parameter mode I.

Definition at line 114 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 126 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 117 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 129 of file BilinearMixedModeCohesiveZoneModel.h.

Referenced by computeQpIProperties(), and computeQpProperties().

◆ _grad_disp

std::vector<const GenericVariableGradient<true> *> BilinearMixedModeCohesiveZoneModel::_grad_disp
protected

Coupled displacement gradients.

Definition at line 74 of file BilinearMixedModeCohesiveZoneModel.h.

Referenced by BilinearMixedModeCohesiveZoneModel(), and computeQpProperties().

◆ _grad_disp_neighbor

std::vector<const GenericVariableGradient<true> *> BilinearMixedModeCohesiveZoneModel::_grad_disp_neighbor
protected

Coupled displacement and neighbor displacement gradient.

Definition at line 77 of file BilinearMixedModeCohesiveZoneModel.h.

Referenced by BilinearMixedModeCohesiveZoneModel(), and 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 147 of file WeightedGapUserObject.h.

Referenced by WeightedVelocitiesUserObject::computeQpProperties(), WeightedGapUserObject::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 190 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 118 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 BilinearMixedModeCohesiveZoneModel::_ndisp
protected

Number of displacement components.

Definition at line 71 of file BilinearMixedModeCohesiveZoneModel.h.

Referenced by BilinearMixedModeCohesiveZoneModel().

◆ _no_iterations

const unsigned int PenaltyWeightedGapUserObject::_no_iterations = 0
staticprotectedinherited

Definition at line 92 of file PenaltyWeightedGapUserObject.h.

◆ _nodal

const bool WeightedGapUserObject::_nodal
protectedinherited

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

Definition at line 140 of file WeightedGapUserObject.h.

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

◆ _normal_strength

const MaterialProperty<Real>& BilinearMixedModeCohesiveZoneModel::_normal_strength
protected

The normal strength material property.

Definition at line 108 of file BilinearMixedModeCohesiveZoneModel.h.

Referenced by computeQpProperties().

◆ _normal_strength_interpolation

ADReal BilinearMixedModeCohesiveZoneModel::_normal_strength_interpolation
protected

Interpolated value of normal_strength.

Definition at line 120 of file BilinearMixedModeCohesiveZoneModel.h.

Referenced by computeQpIProperties(), and computeQpProperties().

◆ _normalization_ptr

const Real* WeightedGapUserObject::_normalization_ptr = nullptr
protectedinherited

Definition at line 181 of file WeightedGapUserObject.h.

◆ _penalty [1/2]

const Real PenaltySimpleCohesiveZoneModel::_penalty
protectedinherited

The normal penalty factor.

Definition at line 51 of file PenaltySimpleCohesiveZoneModel.h.

◆ _penalty [2/2]

const Real PenaltyWeightedGapUserObject::_penalty
protectedinherited

◆ _penalty_friction

const Real PenaltySimpleCohesiveZoneModel::_penalty_friction
protectedinherited

The penalty factor for the frictional constraints.

Definition at line 54 of file PenaltySimpleCohesiveZoneModel.h.

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

◆ _penalty_multiplier

const Real PenaltyWeightedGapUserObject::_penalty_multiplier
protectedinherited

penalty growth factor for augmented Lagrange

Definition at line 79 of file PenaltyWeightedGapUserObject.h.

Referenced by PenaltyWeightedGapUserObject::updateAugmentedLagrangianMultipliers().

◆ _penalty_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 142 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 154 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 158 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 162 of file WeightedGapUserObject.h.

◆ _primary_var

MooseVariableField<Real>& WeightedVelocitiesUserObject::_primary_var
protectedinherited

Reference to the primary variable.

Definition at line 82 of file WeightedVelocitiesUserObject.h.

◆ _primary_x_dot

const ADVariableValue& WeightedVelocitiesUserObject::_primary_x_dot
protectedinherited

x-velocity on the primary face

Definition at line 88 of file WeightedVelocitiesUserObject.h.

Referenced by WeightedVelocitiesUserObject::computeQpProperties().

◆ _primary_y_dot

const ADVariableValue& WeightedVelocitiesUserObject::_primary_y_dot
protectedinherited

y-velocity on the primary face

Definition at line 94 of file WeightedVelocitiesUserObject.h.

Referenced by WeightedVelocitiesUserObject::computeQpProperties().

◆ _primary_z_dot

const ADVariableValue* const WeightedVelocitiesUserObject::_primary_z_dot
protectedinherited

z-velocity on the primary face

Definition at line 100 of file WeightedVelocitiesUserObject.h.

◆ _qp

unsigned int WeightedGapUserObject::_qp = 0
protectedinherited

◆ _qp_displacement_nodal

ADRealVectorValue WeightedGapUserObject::_qp_displacement_nodal
protectedinherited

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

Definition at line 171 of file WeightedGapUserObject.h.

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

◆ _qp_factor

Real WeightedGapUserObject::_qp_factor
protectedinherited

The value of the LM at the current quadrature point.

Definition at line 137 of file WeightedGapUserObject.h.

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

◆ _qp_gap

ADReal WeightedGapUserObject::_qp_gap
protectedinherited

The value of the gap at the current quadrature point.

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

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

◆ _qp_real_tangential_velocity

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

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

Definition at line 67 of file WeightedVelocitiesUserObject.h.

◆ _qp_real_tangential_velocity_nodal

ADRealVectorValue WeightedVelocitiesUserObject::_qp_real_tangential_velocity_nodal
protectedinherited

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

Definition at line 73 of file WeightedVelocitiesUserObject.h.

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

◆ _qp_tangential_velocity

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

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

Definition at line 64 of file WeightedVelocitiesUserObject.h.

◆ _qp_tangential_velocity_nodal

ADRealVectorValue WeightedVelocitiesUserObject::_qp_tangential_velocity_nodal
protectedinherited

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

Definition at line 70 of file WeightedVelocitiesUserObject.h.

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

◆ _regularization_alpha

const Real BilinearMixedModeCohesiveZoneModel::_regularization_alpha
protected

Parameter for the regularization of the Macaulay bracket.

Definition at line 148 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 152 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 156 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 160 of file WeightedGapUserObject.h.

◆ _secondary_var

MooseVariableField<Real>& WeightedVelocitiesUserObject::_secondary_var
protectedinherited

Reference to the secondary variable.

Definition at line 79 of file WeightedVelocitiesUserObject.h.

◆ _secondary_x_dot

const ADVariableValue& WeightedVelocitiesUserObject::_secondary_x_dot
protectedinherited

x-velocity on the secondary face

Definition at line 85 of file WeightedVelocitiesUserObject.h.

Referenced by WeightedVelocitiesUserObject::computeQpProperties().

◆ _secondary_y_dot

const ADVariableValue& WeightedVelocitiesUserObject::_secondary_y_dot
protectedinherited

y-velocity on the secondary face

Definition at line 91 of file WeightedVelocitiesUserObject.h.

Referenced by WeightedVelocitiesUserObject::computeQpProperties().

◆ _secondary_z_dot

const ADVariableValue* const WeightedVelocitiesUserObject::_secondary_z_dot
protectedinherited

z-velocity on the secondary face

Definition at line 97 of file WeightedVelocitiesUserObject.h.

◆ _shear_strength

const MaterialProperty<Real>& BilinearMixedModeCohesiveZoneModel::_shear_strength
protected

The shear strength material property.

Definition at line 111 of file BilinearMixedModeCohesiveZoneModel.h.

Referenced by computeQpProperties().

◆ _shear_strength_interpolation

ADReal BilinearMixedModeCohesiveZoneModel::_shear_strength_interpolation
protected

Interpolated value of shear_strength.

Definition at line 123 of file BilinearMixedModeCohesiveZoneModel.h.

Referenced by computeQpIProperties(), and computeQpProperties().

◆ _sys

SystemBase& WeightedVelocitiesUserObject::_sys
protectedinherited

Reference to the EquationSystem object.

Definition at line 76 of file WeightedVelocitiesUserObject.h.

◆ _tangential_vel_ptr

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

An array of two pointers to avoid copies.

Definition at line 61 of file WeightedVelocitiesUserObject.h.

◆ _test

const VariableTestValue* WeightedGapUserObject::_test = nullptr
protectedinherited

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

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

Definition at line 185 of file WeightedGapUserObject.h.

Referenced by WeightedGapUserObject::execute(), WeightedGapUserObject::initialSetup(), prepareJumpKinematicQuantities(), reinit(), PenaltySimpleCohesiveZoneModel::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 145 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 180 of file WeightedGapUserObject.h.


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