https://mooseframework.inl.gov
Public Types | Public Member Functions | Static Public Member Functions | Public Attributes | Static Public Attributes | Protected Member Functions | Protected Attributes | Static Protected Attributes | Private Member Functions | Private Attributes | List of all members
ExplicitDynamicsContactConstraint Class Reference

A ExplicitDynamicsContactConstraint does mechanical contact for explicit dynamics simulations. More...

#include <ExplicitDynamicsContactConstraint.h>

Inheritance diagram for ExplicitDynamicsContactConstraint:
[legend]

Public Types

enum  ResidualTagType { ResidualTagType::NonReference, ResidualTagType::Reference }
 
typedef DataFileName DataFileParameterType
 

Public Member Functions

 ExplicitDynamicsContactConstraint (const InputParameters &parameters)
 
virtual void timestepSetup () override
 
virtual void jacobianSetup () override
 
virtual void residualEnd () override
 
virtual void updateContactStatefulData (bool beginning_of_step)
 
virtual Real computeQpSecondaryValue () override
 
virtual Real computeQpResidual (Moose::ConstraintType type) override
 
virtual void computeJacobian () override
 Computes the jacobian for the current element. More...
 
virtual void computeOffDiagJacobian (unsigned int) override
 Compute off-diagonal Jacobian entries. More...
 
virtual Real computeQpJacobian (Moose::ConstraintJacobianType) override
 
virtual Real computeQpOffDiagJacobian (Moose::ConstraintJacobianType, unsigned int) override
 Compute off-diagonal Jacobian entries. More...
 
bool shouldApply () override
 
void computeContactForce (const Node &node, PenetrationInfo *pinfo, bool update_contact_set)
 
virtual bool isExplicitConstraint () const override
 
virtual bool addCouplingEntriesToJacobian () override
 Return false so that the nonlinear system does not try to add Jacobian entries from the contact forces. More...
 
virtual const std::unordered_set< unsigned int > & getMatPropDependencies () const override
 
virtual void overwriteBoundaryVariables (NumericVector< Number > &, const Node &) const override
 
virtual void computeSecondaryValue (NumericVector< Number > &current_solution)
 
virtual void computeResidual () override
 
virtual void getConnectedDofIndices (unsigned int var_num)
 
virtual bool overwriteSecondaryResidual ()
 
virtual bool overwriteSecondaryJacobian ()
 
virtual MooseVariableprimaryVariable ()
 
BoundaryID primaryBoundary () const
 
BoundaryID secondaryBoundary () const
 
const MooseVariablevariable () const override
 
Real secondaryResidual () const
 
void residualSetup () override
 
std::set< SubdomainIDgetSecondaryConnectedBlocks () const
 
virtual void subdomainSetup () override final
 
void prepareNeighborShapes (unsigned int var_num)
 
virtual void computeResidualAndJacobian ()
 
virtual void computeOffDiagJacobianScalar (unsigned int)
 
virtual void computeNonlocalJacobian ()
 
virtual void computeNonlocalOffDiagJacobian (unsigned int)
 
const SubProblemsubProblem () const
 
virtual void prepareShapes (unsigned int var_num)
 
virtual std::set< std::string > additionalROVariables ()
 
virtual bool enabled () const
 
std::shared_ptr< MooseObjectgetSharedPtr ()
 
std::shared_ptr< const MooseObjectgetSharedPtr () const
 
MooseAppgetMooseApp () const
 
const std::string & type () const
 
virtual const std::string & name () const
 
std::string typeAndName () const
 
std::string errorPrefix (const std::string &error_type) const
 
void callMooseError (std::string msg, const bool with_prefix) const
 
MooseObjectParameterName uniqueParameterName (const std::string &parameter_name) const
 
const InputParametersparameters () const
 
MooseObjectName uniqueName () const
 
const T & getParam (const std::string &name) const
 
std::vector< std::pair< T1, T2 > > getParam (const std::string &param1, const std::string &param2) const
 
const T * queryParam (const std::string &name) const
 
const T & getRenamedParam (const std::string &old_name, const std::string &new_name) const
 
getCheckedPointerParam (const std::string &name, const std::string &error_string="") const
 
bool isParamValid (const std::string &name) const
 
bool isParamSetByUser (const std::string &nm) const
 
void paramError (const std::string &param, Args... args) const
 
void paramWarning (const std::string &param, Args... args) const
 
void paramInfo (const std::string &param, Args... args) const
 
void connectControllableParams (const std::string &parameter, const std::string &object_type, const std::string &object_name, const std::string &object_parameter) const
 
void mooseError (Args &&... args) const
 
void mooseErrorNonPrefixed (Args &&... args) const
 
void mooseDocumentedError (const std::string &repo_name, const unsigned int issue_num, Args &&... args) const
 
void mooseWarning (Args &&... args) const
 
void mooseWarningNonPrefixed (Args &&... args) const
 
void mooseDeprecated (Args &&... args) const
 
void mooseInfo (Args &&... args) const
 
std::string getDataFileName (const std::string &param) const
 
std::string getDataFileNameByName (const std::string &relative_path) const
 
std::string getDataFilePath (const std::string &relative_path) const
 
virtual void initialSetup ()
 
virtual void customSetup (const ExecFlagType &)
 
const ExecFlagEnumgetExecuteOnEnum () const
 
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
 
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
 
bool isImplicit ()
 
Moose::StateArg determineState () 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
 
void setRandomResetFrequency (ExecFlagType exec_flag)
 
unsigned long getRandomLong () const
 
Real getRandomReal () const
 
unsigned int getSeed (std::size_t id)
 
unsigned int getMasterSeed () const
 
bool isNodal () const
 
ExecFlagType getResetOnTime () const
 
void setRandomDataPointer (RandomData *random_data)
 
virtual void meshChanged ()
 
void useVectorTag (const TagName &tag_name, VectorTagsKey)
 
void useVectorTag (TagID tag_id, VectorTagsKey)
 
void useMatrixTag (const TagName &tag_name, MatrixTagsKey)
 
void useMatrixTag (TagID tag_id, MatrixTagsKey)
 
bool isVectorTagged ()
 
bool isMatrixTagged ()
 
bool hasVectorTags () const
 
const std::set< TagID > & getVectorTags (VectorTagsKey) const
 
const std::set< TagID > & getMatrixTags (MatrixTagsKey) const
 
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
 
const std::vector< MooseVariableScalar *> & getCoupledMooseScalarVars ()
 
const std::set< TagID > & getScalarVariableCoupleableVectorTags () const
 
const std::set< TagID > & getScalarVariableCoupleableMatrixTags () const
 
const std::set< MooseVariableFieldBase *> & getMooseVariableDependencies () const
 
std::set< MooseVariableFieldBase *> checkAllVariables (const DofObjectType &dof_object, const std::set< MooseVariableFieldBase * > &vars_to_omit={})
 
std::set< MooseVariableFieldBase *> checkVariables (const DofObjectType &dof_object, const std::set< MooseVariableFieldBase * > &vars_to_check)
 
void addMooseVariableDependency (MooseVariableFieldBase *var)
 
void addMooseVariableDependency (const std::vector< MooseVariableFieldBase * > &vars)
 
MooseVariableBasemooseVariableBase () const
 
MooseVariableField< Real > & mooseVariableField ()
 
MooseVariableFE< Real > * mooseVariable () const
 
MooseVariableFV< Real > * mooseVariableFV () const
 
MooseLinearVariableFV< Real > * mooseLinearVariableFV () 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 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
 
PenetrationLocatorgetPenetrationLocator (const BoundaryName &primary, const BoundaryName &secondary, Order order)
 
PenetrationLocatorgetQuadraturePenetrationLocator (const BoundaryName &primary, const BoundaryName &secondary, Order order)
 
NearestNodeLocatorgetNearestNodeLocator (const BoundaryName &primary, const BoundaryName &secondary)
 
NearestNodeLocatorgetQuadratureNearestNodeLocator (const BoundaryName &primary, const BoundaryName &secondary)
 
bool requiresGeometricSearch () const
 
const Parallel::Communicator & comm () const
 
processor_id_type n_processors () const
 
processor_id_type processor_id () const
 
const GenericMaterialProperty< T, is_ad > & getGenericMaterialProperty (const std::string &name, MaterialData &material_data, const unsigned int state=0)
 
const GenericMaterialProperty< T, is_ad > & getGenericMaterialProperty (const std::string &name, const unsigned int state=0)
 
const GenericMaterialProperty< T, is_ad > & getGenericMaterialProperty (const std::string &name, const unsigned int state=0)
 
const MaterialProperty< T > & getMaterialProperty (const std::string &name, MaterialData &material_data, const unsigned int state=0)
 
const MaterialProperty< T > & getMaterialProperty (const std::string &name, const unsigned int state=0)
 
const MaterialProperty< T > & getMaterialProperty (const std::string &name, const unsigned int state=0)
 
const ADMaterialProperty< T > & getADMaterialProperty (const std::string &name, MaterialData &material_data)
 
const ADMaterialProperty< T > & getADMaterialProperty (const std::string &name)
 
const ADMaterialProperty< T > & getADMaterialProperty (const std::string &name)
 
const MaterialProperty< T > & getMaterialPropertyOld (const std::string &name, MaterialData &material_data)
 
const MaterialProperty< T > & getMaterialPropertyOld (const std::string &name)
 
const MaterialProperty< T > & getMaterialPropertyOld (const std::string &name)
 
const MaterialProperty< T > & getMaterialPropertyOlder (const std::string &name, MaterialData &material_data)
 
const MaterialProperty< T > & getMaterialPropertyOlder (const std::string &name)
 
const MaterialProperty< T > & getMaterialPropertyOlder (const std::string &name)
 
const GenericMaterialProperty< T, is_ad > & getGenericMaterialPropertyByName (const MaterialPropertyName &name, MaterialData &material_data, const unsigned int state)
 
const GenericMaterialProperty< T, is_ad > & getGenericMaterialPropertyByName (const MaterialPropertyName &name, const unsigned int state=0)
 
const GenericMaterialProperty< T, is_ad > & getGenericMaterialPropertyByName (const MaterialPropertyName &name, const unsigned int state=0)
 
const MaterialProperty< T > & getMaterialPropertyByName (const MaterialPropertyName &name, MaterialData &material_data, const unsigned int state=0)
 
const MaterialProperty< T > & getMaterialPropertyByName (const MaterialPropertyName &name, const unsigned int state=0)
 
const MaterialProperty< T > & getMaterialPropertyByName (const MaterialPropertyName &name, const unsigned int state=0)
 
const ADMaterialProperty< T > & getADMaterialPropertyByName (const MaterialPropertyName &name, MaterialData &material_data)
 
const ADMaterialProperty< T > & getADMaterialPropertyByName (const MaterialPropertyName &name)
 
const ADMaterialProperty< T > & getADMaterialPropertyByName (const MaterialPropertyName &name)
 
const MaterialProperty< T > & getMaterialPropertyOldByName (const MaterialPropertyName &name, MaterialData &material_data)
 
const MaterialProperty< T > & getMaterialPropertyOldByName (const MaterialPropertyName &name)
 
const MaterialProperty< T > & getMaterialPropertyOldByName (const MaterialPropertyName &name)
 
const MaterialProperty< T > & getMaterialPropertyOlderByName (const MaterialPropertyName &name, MaterialData &material_data)
 
const MaterialProperty< T > & getMaterialPropertyOlderByName (const MaterialPropertyName &name)
 
const MaterialProperty< T > & getMaterialPropertyOlderByName (const MaterialPropertyName &name)
 
std::pair< const MaterialProperty< T > *, std::set< SubdomainID > > getBlockMaterialProperty (const MaterialPropertyName &name)
 
const GenericMaterialProperty< T, is_ad > & getGenericZeroMaterialProperty (const std::string &name)
 
const GenericMaterialProperty< T, is_ad > & getGenericZeroMaterialProperty ()
 
const GenericMaterialProperty< T, is_ad > & getGenericZeroMaterialPropertyByName (const std::string &prop_name)
 
const MaterialProperty< T > & getZeroMaterialProperty (Ts... args)
 
std::set< 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)
 
bool getMaterialPropertyCalled () const
 
virtual void resolveOptionalProperties ()
 
const GenericMaterialProperty< T, is_ad > & getPossiblyConstantGenericMaterialPropertyByName (const MaterialPropertyName &prop_name, MaterialData &material_data, const unsigned int state)
 
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 hasGenericMaterialProperty (const std::string &name)
 
bool hasGenericMaterialProperty (const std::string &name)
 
bool hasGenericMaterialPropertyByName (const std::string &name)
 
bool hasGenericMaterialPropertyByName (const std::string &name)
 

Static Public Member Functions

static InputParameters validParams ()
 

Public Attributes

const ConsoleStream _console
 

Static Public Attributes

static constexpr PropertyValue::id_type default_property_id
 
static constexpr PropertyValue::id_type zero_property_id
 

Protected Member Functions

void solveImpactEquations (const Node &node, PenetrationInfo *pinfo, const RealVectorValue &distance_gap)
 Determine "Lagrange multipliers" from the iterative solution of the impact problem. More...
 
Real gapOffset (const Node &node)
 
Real getPenalty (const Node &node)
 
virtual const VariableValuecoupledSecondaryValue (const std::string &var_name, unsigned int comp=0)
 
virtual const VariableValuecoupledSecondaryValueOld (const std::string &var_name, unsigned int comp=0)
 
virtual const VariableValuecoupledSecondaryValueOlder (const std::string &var_name, unsigned int comp=0)
 
virtual const VariableGradientcoupledSecondaryGradient (const std::string &var_name, unsigned int comp=0)
 
virtual const VariableGradientcoupledSecondaryGradientOld (const std::string &var_name, unsigned int comp=0)
 
virtual const VariableGradientcoupledSecondaryGradientOlder (const std::string &var_name, unsigned int comp=0)
 
virtual const VariableSecondcoupledSecondarySecond (const std::string &var_name, unsigned int comp=0)
 
virtual const VariableValuecoupledPrimaryValue (const std::string &var_name, unsigned int comp=0)
 
virtual const VariableValuecoupledPrimaryValueOld (const std::string &var_name, unsigned int comp=0)
 
virtual const VariableValuecoupledPrimaryValueOlder (const std::string &var_name, unsigned int comp=0)
 
virtual const VariableGradientcoupledPrimaryGradient (const std::string &var_name, unsigned int comp=0)
 
virtual const VariableGradientcoupledPrimaryGradientOld (const std::string &var_name, unsigned int comp=0)
 
virtual const VariableGradientcoupledPrimaryGradientOlder (const std::string &var_name, unsigned int comp=0)
 
virtual const VariableSecondcoupledPrimarySecond (const std::string &var_name, unsigned int comp=0)
 
const std::set< BoundaryID > & buildBoundaryIDs ()
 
virtual void precalculateResidual ()
 
virtual void precalculateJacobian ()
 
virtual void precalculateOffDiagJacobian (unsigned int)
 
const MooseVariableFieldBasegetVariable (unsigned int jvar_num) const
 
virtual void addUserObjectDependencyHelper (const UserObject &) const
 
virtual void addPostprocessorDependencyHelper (const PostprocessorName &) const
 
virtual void addVectorPostprocessorDependencyHelper (const VectorPostprocessorName &) 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
 
void prepareVectorTag (Assembly &assembly, unsigned int ivar)
 
void prepareVectorTag (Assembly &assembly, unsigned int ivar, ResidualTagType tag_type)
 
void prepareVectorTag (Assembly &assembly, unsigned int ivar, ResidualTagType tag_type)
 
void prepareVectorTag (Assembly &assembly, unsigned int ivar, ResidualTagType tag_type)
 
void prepareVectorTag (Assembly &assembly, unsigned int ivar, ResidualTagType tag_type)
 
void prepareVectorTagNeighbor (Assembly &assembly, unsigned int ivar)
 
void prepareVectorTagLower (Assembly &assembly, unsigned int ivar)
 
void prepareMatrixTag (Assembly &assembly, unsigned int ivar, unsigned int jvar)
 
void prepareMatrixTag (Assembly &assembly, unsigned int ivar, unsigned int jvar, DenseMatrix< Number > &k) const
 
void prepareMatrixTagNonlocal (Assembly &assembly, unsigned int ivar, unsigned int jvar)
 
void prepareMatrixTagNeighbor (Assembly &assembly, unsigned int ivar, unsigned int jvar, Moose::DGJacobianType type)
 
void prepareMatrixTagNeighbor (Assembly &assembly, unsigned int ivar, unsigned int jvar, Moose::DGJacobianType type, DenseMatrix< Number > &k) const
 
void prepareMatrixTagLower (Assembly &assembly, unsigned int ivar, unsigned int jvar, Moose::ConstraintJacobianType type)
 
void accumulateTaggedLocalResidual ()
 
void assignTaggedLocalResidual ()
 
void accumulateTaggedLocalMatrix ()
 
void accumulateTaggedLocalMatrix (Assembly &assembly, unsigned int ivar, unsigned int jvar, const DenseMatrix< Number > &k)
 
void accumulateTaggedLocalMatrix (Assembly &assembly, unsigned int ivar, unsigned int jvar, Moose::DGJacobianType type, const DenseMatrix< Number > &k)
 
void accumulateTaggedNonlocalMatrix ()
 
void assignTaggedLocalMatrix ()
 
void addResiduals (Assembly &assembly, const Residuals &residuals, const Indices &dof_indices, Real scaling_factor)
 
void addResiduals (Assembly &assembly, const DenseVector< T > &residuals, const Indices &dof_indices, Real scaling_factor)
 
void addResiduals (Assembly &assembly, const ADResidualsPacket &packet)
 
void addResidualsAndJacobian (Assembly &assembly, const Residuals &residuals, const Indices &dof_indices, Real scaling_factor)
 
void addResidualsAndJacobian (Assembly &assembly, const ADResidualsPacket &packet)
 
void addJacobian (Assembly &assembly, const Residuals &residuals, const Indices &dof_indices, Real scaling_factor)
 
void addJacobian (Assembly &assembly, const ADResidualsPacket &packet)
 
void addJacobian (Assembly &assembly, DenseMatrix< Real > &local_k, const std::vector< dof_id_type > &row_indices, const std::vector< dof_id_type > &column_indices, Real scaling_factor)
 
void addResidualsWithoutConstraints (Assembly &assembly, const Residuals &residuals, const Indices &dof_indices, Real scaling_factor)
 
void addResidualsAndJacobianWithoutConstraints (Assembly &assembly, const Residuals &residuals, const Indices &dof_indices, Real scaling_factor)
 
void addJacobianWithoutConstraints (Assembly &assembly, const Residuals &residuals, const Indices &dof_indices, Real scaling_factor)
 
void addJacobianElement (Assembly &assembly, Real value, dof_id_type row_index, dof_id_type column_index, Real scaling_factor)
 
void setResidual (SystemBase &sys, const T &residual, MooseVariableFE< T > &var)
 
void setResidual (SystemBase &sys, Real residual, dof_id_type dof_index)
 
void setResidual (SystemBase &sys, SetResidualFunctor set_residual_functor)
 
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
 
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
 
virtual const OutputTools< Real >::VariableValueneighborValue ()
 
const VectorVariableValueneighborValue ()
 
virtual const OutputTools< Real >::VariableValueneighborValueOld ()
 
const VectorVariableValueneighborValueOld ()
 
virtual const OutputTools< Real >::VariableValueneighborValueOlder ()
 
const VectorVariableValueneighborValueOlder ()
 
virtual const OutputTools< Real >::VariableGradientneighborGradient ()
 
virtual const OutputTools< Real >::VariableGradientneighborGradientOld ()
 
virtual const OutputTools< Real >::VariableGradientneighborGradientOlder ()
 
virtual const OutputTools< Real >::VariableSecondneighborSecond ()
 
virtual const OutputTools< Real >::VariableSecondneighborSecondOld ()
 
virtual const OutputTools< Real >::VariableSecondneighborSecondOlder ()
 
virtual const OutputTools< Real >::VariableTestSecondneighborSecondTest ()
 
virtual const OutputTools< Real >::VariablePhiSecondneighborSecondPhi ()
 
virtual const OutputTools< Real >::VariableValuevalue ()
 
virtual const OutputTools< Real >::VariableValuevalueOld ()
 
virtual const OutputTools< Real >::VariableValuevalueOlder ()
 
virtual const OutputTools< Real >::VariableValuedot ()
 
virtual const OutputTools< Real >::VariableValuedotDot ()
 
virtual const OutputTools< Real >::VariableValuedotOld ()
 
virtual const OutputTools< Real >::VariableValuedotDotOld ()
 
virtual const VariableValuedotDu ()
 
virtual const VariableValuedotDotDu ()
 
virtual const OutputTools< Real >::VariableGradientgradient ()
 
virtual const OutputTools< Real >::VariableGradientgradientOld ()
 
virtual const OutputTools< Real >::VariableGradientgradientOlder ()
 
virtual const OutputTools< Real >::VariableSecondsecond ()
 
virtual const OutputTools< Real >::VariableSecondsecondOld ()
 
virtual const OutputTools< Real >::VariableSecondsecondOlder ()
 
virtual const OutputTools< Real >::VariableTestSecondsecondTest ()
 
virtual const OutputTools< Real >::VariableTestSecondsecondTestFace ()
 
virtual const OutputTools< Real >::VariablePhiSecondsecondPhi ()
 
virtual const OutputTools< Real >::VariablePhiSecondsecondPhiFace ()
 
virtual void checkMaterialProperty (const std::string &name, const unsigned int state)
 
void markMatPropRequested (const std::string &)
 
MaterialPropertyName getMaterialPropertyName (const std::string &name) const
 
void checkExecutionStage ()
 
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)
 

Protected Attributes

MooseSharedPointer< DisplacedProblem_displaced_problem
 
const unsigned int _component
 
const ExplicitDynamicsContactModel _model
 
bool _update_stateful_data
 
const unsigned int _mesh_dimension
 
std::vector< unsigned int_vars
 
std::vector< MooseVariable * > _var_objects
 
const bool _has_secondary_gap_offset
 
const MooseVariable *const _secondary_gap_offset_var
 
const bool _has_mapped_primary_gap_offset
 
const MooseVariable *const _mapped_primary_gap_offset_var
 
const Real _penalty
 
const bool _print_contact_nodes
 
NumericVector< Number > & _residual_copy
 
MooseWritableVariable_gap_rate
 Nodal gap rate (output for debugging or analyst perusal) More...
 
const VariableValue_neighbor_vel_x
 X component of velocity at the closest point. More...
 
const VariableValue_neighbor_vel_y
 Y component of velocity at the closest point. More...
 
const VariableValue_neighbor_vel_z
 Z component of velocity at the closest point. More...
 
const SparseMatrix< Number > * _jacobian
 
BoundaryID _secondary
 
BoundaryID _primary
 
MooseVariable_var
 
const MooseArray< Point > & _primary_q_point
 
const QBase *const & _primary_qrule
 
std::set< BoundaryID_boundary_ids
 
PenetrationLocator_penetration_locator
 
const Node *const & _current_node
 
const Elem *const & _current_primary
 
const VariableValue_u_secondary
 
VariablePhiValue _phi_secondary
 
VariableTestValue _test_secondary
 
MooseVariable_primary_var
 
unsigned int _primary_var_num
 
const VariablePhiValue_phi_primary
 
const VariablePhiGradient_grad_phi_primary
 
const VariableTestValue_test_primary
 
const VariableTestGradient_grad_test_primary
 
const VariableValue_u_primary
 
const VariableGradient_grad_u_primary
 
const DofMap & _dof_map
 
const std::map< dof_id_type, std::vector< dof_id_type > > & _node_to_elem_map
 
bool _overwrite_secondary_residual
 
const MooseArray< Real > & _primary_JxW
 
bool _secondary_residual_computed
 
Real _secondary_residual
 
const std::unordered_set< unsigned int_empty_mat_prop_deps
 
std::vector< dof_id_type_connected_dof_indices
 
DenseMatrix< Number_Kne
 
DenseMatrix< Number_Kee
 
DenseMatrix< Number_Ken
 
unsigned int _i
 
unsigned int _j
 
unsigned int _qp
 
SubProblem_subproblem
 
FEProblemBase_fe_problem
 
SystemBase_sys
 
THREAD_ID _tid
 
Assembly_assembly
 
MooseMesh_mesh
 
const bool & _enabled
 
MooseApp_app
 
const std::string _type
 
const std::string _name
 
const InputParameters_pars
 
Factory_factory
 
ActionFactory_action_factory
 
const ExecFlagEnum_execute_enum
 
const ExecFlagType_current_execute_flag
 
const InputParameters_ti_params
 
FEProblemBase_ti_feproblem
 
bool _is_implicit
 
Real_t
 
const Real_t_old
 
int_t_step
 
Real_dt
 
Real_dt_old
 
bool _is_transient
 
MooseApp_restartable_app
 
const std::string _restartable_system_name
 
const THREAD_ID _restartable_tid
 
const bool _restartable_read_only
 
FEProblemBase_mci_feproblem
 
DenseVector< Number_local_re
 
DenseMatrix< Number_local_ke
 
DenseMatrix< Number_nonlocal_ke
 
GeometricSearchData_geometric_search_data
 
bool _requires_geometric_search
 
bool _neighbor_nodal
 
const InputParameters_c_parameters
 
const std::string & _c_name
 
const std::string & _c_type
 
FEProblemBase_c_fe_problem
 
const SystemBase *const _c_sys
 
std::unordered_map< std::string, std::vector< MooseVariableFieldBase *> > _coupled_vars
 
std::vector< MooseVariableFieldBase *> _coupled_moose_vars
 
std::vector< MooseVariable *> _coupled_standard_moose_vars
 
std::vector< VectorMooseVariable *> _coupled_vector_moose_vars
 
std::vector< ArrayMooseVariable *> _coupled_array_moose_vars
 
std::vector< MooseVariableFV< Real > *> _coupled_standard_fv_moose_vars
 
std::vector< MooseLinearVariableFV< Real > *> _coupled_standard_linear_fv_moose_vars
 
const std::unordered_map< std::string, std::string > & _new_to_deprecated_coupled_vars
 
bool _c_nodal
 
bool _c_is_implicit
 
const bool _c_allow_element_to_nodal_coupling
 
THREAD_ID _c_tid
 
std::unordered_map< std::string, std::vector< std::unique_ptr< VariableValue > > > _default_value
 
std::unordered_map< std::string, std::unique_ptr< MooseArray< ADReal > > > _ad_default_value
 
std::unordered_map< std::string, std::unique_ptr< VectorVariableValue > > _default_vector_value
 
std::unordered_map< std::string, std::unique_ptr< ArrayVariableValue > > _default_array_value
 
std::unordered_map< std::string, std::unique_ptr< MooseArray< ADRealVectorValue > > > _ad_default_vector_value
 
VariableValue _default_value_zero
 
VariableGradient _default_gradient
 
MooseArray< ADRealVectorValue_ad_default_gradient
 
MooseArray< ADRealTensorValue_ad_default_vector_gradient
 
VariableSecond _default_second
 
MooseArray< ADRealTensorValue_ad_default_second
 
MooseArray< ADRealVectorValue_ad_default_curl
 
const VariableValue_zero
 
const VariablePhiValue_phi_zero
 
const MooseArray< ADReal > & _ad_zero
 
const VariableGradient_grad_zero
 
const MooseArray< ADRealVectorValue > & _ad_grad_zero
 
const VariablePhiGradient_grad_phi_zero
 
const VariableSecond_second_zero
 
const MooseArray< ADRealTensorValue > & _ad_second_zero
 
const VariablePhiSecond_second_phi_zero
 
const VectorVariableValue_vector_zero
 
const VectorVariableCurl_vector_curl_zero
 
VectorVariableValue _default_vector_value_zero
 
VectorVariableGradient _default_vector_gradient
 
VectorVariableCurl _default_vector_curl
 
VectorVariableDivergence _default_div
 
ArrayVariableValue _default_array_value_zero
 
ArrayVariableGradient _default_array_gradient
 
bool _coupleable_neighbor
 
FEProblemBase_sc_fe_problem
 
const THREAD_ID _sc_tid
 
const Real_real_zero
 
const VariableValue_scalar_zero
 
const Point & _point_zero
 
bool _nodal
 
MooseVariableFE< Real > * _variable
 
MooseVariableFV< Real > * _fv_variable
 
MooseLinearVariableFV< Real > * _linear_fv_variable
 
MooseVariableField< Real > * _field_variable
 
Assembly_mvi_assembly
 
const Parallel::Communicator & _communicator
 
MaterialData_neighbor_material_data
 
const InputParameters_mi_params
 
const std::string _mi_name
 
const MooseObjectName _mi_moose_object_name
 
FEProblemBase_mi_feproblem
 
SubProblem_mi_subproblem
 
const THREAD_ID _mi_tid
 
const Moose::MaterialDataType _material_data_type
 
MaterialData_material_data
 
bool _stateful_allowed
 
bool _get_material_property_called
 
std::vector< std::unique_ptr< PropertyValue > > _default_properties
 
std::unordered_set< unsigned int_material_property_dependencies
 
const MaterialPropertyName _get_suffix
 
const bool _use_interpolated_state
 

Static Protected Attributes

static const unsigned int _no_iterations = 0
 
static const std::string _interpolated_old
 
static const std::string _interpolated_older
 

Private Member Functions

Real computeFaceMass (const NumericVector< Real > &lumped_mass)
 

Private Attributes

std::unordered_map< dof_id_type, Real_dof_to_position
 

Detailed Description

A ExplicitDynamicsContactConstraint does mechanical contact for explicit dynamics simulations.

Definition at line 25 of file ExplicitDynamicsContactConstraint.h.

Constructor & Destructor Documentation

◆ ExplicitDynamicsContactConstraint()

ExplicitDynamicsContactConstraint::ExplicitDynamicsContactConstraint ( const InputParameters parameters)

Definition at line 70 of file ExplicitDynamicsContactConstraint.C.

74  _displaced_problem(parameters.get<FEProblemBase *>("_fe_problem_base")->getDisplacedProblem()),
75  _component(getParam<unsigned int>("component")),
76  _model(getParam<MooseEnum>("model").getEnum<ExplicitDynamicsContactModel>()),
80  _var_objects(3, nullptr),
81  _has_secondary_gap_offset(isCoupled("secondary_gap_offset")),
83  : nullptr),
84  _has_mapped_primary_gap_offset(isCoupled("mapped_primary_gap_offset")),
86  _has_mapped_primary_gap_offset ? getVar("mapped_primary_gap_offset", 0) : nullptr),
87  _penalty(getParam<Real>("penalty")),
88  _print_contact_nodes(getParam<bool>("print_contact_nodes")),
90  _gap_rate(&writableVariable("gap_rate")),
91  _neighbor_vel_x(isCoupled("vel_x") ? coupledNeighborValue("vel_x") : _zero),
92  _neighbor_vel_y(isCoupled("vel_y") ? coupledNeighborValue("vel_y") : _zero),
93  _neighbor_vel_z((_mesh.dimension() == 3 && isCoupled("vel_z")) ? coupledNeighborValue("vel_z")
94  : _zero)
95 {
97 
98  if (isParamValid("displacements"))
99  {
100  // modern parameter scheme for displacements
101  for (unsigned int i = 0; i < coupledComponents("displacements"); ++i)
102  {
103  _vars[i] = coupled("displacements", i);
104  _var_objects[i] = getVar("displacements", i);
105  }
106  }
107 
108  if (parameters.isParamValid("tangential_tolerance"))
109  _penetration_locator.setTangentialTolerance(getParam<Real>("tangential_tolerance"));
110 
111  if (parameters.isParamValid("normal_smoothing_distance"))
112  _penetration_locator.setNormalSmoothingDistance(getParam<Real>("normal_smoothing_distance"));
113 
114  if (parameters.isParamValid("normal_smoothing_method"))
116  parameters.get<std::string>("normal_smoothing_method"));
117 
119  {
120  bool is_correct =
121  (isCoupled("vel_x") && isCoupled("vel_y") && _mesh.dimension() == 2) ||
122  (isCoupled("vel_x") && isCoupled("vel_y") && isCoupled("vel_z") && _mesh.dimension() == 3);
123 
124  if (!is_correct)
125  paramError("vel_x",
126  "Velocities vel_x and vel_y (also vel_z in three dimensions) need to be provided "
127  "for the 'balance' option of solving normal contact in explicit dynamics.");
128  }
129 }
MooseSharedPointer< DisplacedProblem > _displaced_problem
MooseMesh & _mesh
virtual bool isCoupled(const std::string &var_name, unsigned int i=0) const
const VariableValue & _zero
virtual unsigned int coupled(const std::string &var_name, unsigned int comp=0) const
const unsigned int invalid_uint
void setNormalSmoothingDistance(Real normal_smoothing_distance)
std::vector< std::pair< R1, R2 > > get(const std::string &param1, const std::string &param2) const
MooseVariable * getVar(const std::string &var_name, unsigned int comp)
MooseWritableVariable & writableVariable(const std::string &var_name, unsigned int comp=0)
void setTangentialTolerance(Real tangential_tolerance)
bool isParamValid(const std::string &name) const
SystemBase & _sys
TwoMaterialPropertyInterface(const MooseObject *moose_object, const std::set< SubdomainID > &blocks_ids, const std::set< BoundaryID > &boundary_ids)
NodeFaceConstraint(const InputParameters &parameters)
virtual const VariableValue & coupledNeighborValue(const std::string &var_name, unsigned int comp=0) const
virtual unsigned int dimension() const
const ExplicitDynamicsContactModel _model
void paramError(const std::string &param, Args... args) const
const std::set< BoundaryID > & buildBoundaryIDs()
MooseWritableVariable * _gap_rate
Nodal gap rate (output for debugging or analyst perusal)
unsigned int coupledComponents(const std::string &var_name) const
const std::set< SubdomainID > EMPTY_BLOCK_IDS
const VariableValue & _neighbor_vel_z
Z component of velocity at the closest point.
const InputParameters & parameters() const
bool _overwrite_secondary_residual
const VariableValue & _neighbor_vel_y
Y component of velocity at the closest point.
const VariableValue & _neighbor_vel_x
X component of velocity at the closest point.
void setNormalSmoothingMethod(std::string nsmString)
PenetrationLocator & _penetration_locator
virtual NumericVector< Number > & residualGhosted()
bool isParamValid(const std::string &name) const

Member Function Documentation

◆ addCouplingEntriesToJacobian()

virtual bool ExplicitDynamicsContactConstraint::addCouplingEntriesToJacobian ( )
inlineoverridevirtual

Return false so that the nonlinear system does not try to add Jacobian entries from the contact forces.

Returns
bool indicating whether we need to couple Jacobian entries

Reimplemented from NodeFaceConstraint.

Definition at line 73 of file ExplicitDynamicsContactConstraint.h.

73 { return false; }

◆ computeContactForce()

void ExplicitDynamicsContactConstraint::computeContactForce ( const Node &  node,
PenetrationInfo pinfo,
bool  update_contact_set 
)

Definition at line 210 of file ExplicitDynamicsContactConstraint.C.

Referenced by shouldApply().

213 {
214  RealVectorValue distance_vec(node - pinfo->_closest_point);
215 
216  // Adding current velocity to the distance vector to ensure proper contact check
217  dof_id_type dof_x = node.dof_number(_sys.number(), _var_objects[0]->number(), 0);
218  dof_id_type dof_y = node.dof_number(_sys.number(), _var_objects[1]->number(), 0);
219  dof_id_type dof_z =
220  _mesh.dimension() == 3 ? node.dof_number(_sys.number(), _var_objects[2]->number(), 0) : 0;
221  RealVectorValue udotvec = {(*_sys.solutionUDot())(dof_x)*_dt,
222  (*_sys.solutionUDot())(dof_y)*_dt,
223  _mesh.dimension() == 3 ? (*_sys.solutionUDot())(dof_z)*_dt : 0};
224  distance_vec += udotvec;
225 
226  if (distance_vec.norm() != 0)
227  distance_vec += gapOffset(node) * pinfo->_normal;
228 
229  const Real gap_size = -1.0 * pinfo->_normal * distance_vec;
230 
231  // This is for preventing an increment of pinfo->_locked_this_step for nodes that are
232  // captured and released in this function
233  bool newly_captured = false;
234 
235  // Capture nodes that are newly in contact
236  if (update_contact_set && !pinfo->isCaptured() &&
237  MooseUtils::absoluteFuzzyGreaterEqual(gap_size, 0.0, 0.0))
238  {
239  newly_captured = true;
240  pinfo->capture();
241  }
242 
243  if (!pinfo->isCaptured())
244  return;
245 
246  switch (_model)
247  {
249  {
250  const Real penalty = getPenalty(node);
251  RealVectorValue pen_force(penalty * distance_vec);
252  pinfo->_contact_force = pinfo->_normal * (pinfo->_normal * pen_force);
253  break;
254  }
256  solveImpactEquations(node, pinfo, distance_vec);
257  break;
258  default:
259  mooseError("Invalid or unavailable contact model");
260  break;
261  }
262 
263  if (update_contact_set && pinfo->isCaptured() && !newly_captured)
264  {
265  const Real contact_pressure = -(pinfo->_normal * pinfo->_contact_force);
266  if (-contact_pressure >= 0.0)
267  {
268  pinfo->release();
269  pinfo->_contact_force.zero();
270  }
271  }
272 }
MooseMesh & _mesh
RealVectorValue _normal
void solveImpactEquations(const Node &node, PenetrationInfo *pinfo, const RealVectorValue &distance_gap)
Determine "Lagrange multipliers" from the iterative solution of the impact problem.
SystemBase & _sys
virtual unsigned int dimension() const
virtual NumericVector< Number > * solutionUDot()
bool isCaptured() const
const ExplicitDynamicsContactModel _model
unsigned int number() const
RealVectorValue _contact_force
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
bool absoluteFuzzyGreaterEqual(const T &var1, const T2 &var2, const T3 &tol=libMesh::TOLERANCE *libMesh::TOLERANCE)
void mooseError(Args &&... args) const
uint8_t dof_id_type

◆ computeFaceMass()

Real ExplicitDynamicsContactConstraint::computeFaceMass ( const NumericVector< Real > &  lumped_mass)
private
Parameters
lumped_massThe inverted lumped mass vector
Returns
Effictive mass of the face

Definition at line 378 of file ExplicitDynamicsContactConstraint.C.

Referenced by solveImpactEquations().

379 {
380  // Initialize face mass to zero
381  Real mass_face(0.0);
382 
383  // Get the primary side of the current contact
384  const auto primary_side = _current_primary->side_ptr(_assembly.neighborSide());
385 
386  // Get the dofs on the primary (face) side of the contact
387  std::vector<dof_id_type> face_dofs;
388  _primary_var.getDofIndices(primary_side.get(), face_dofs);
389 
390  // Get average mass of face
391  for (const auto dof : face_dofs)
392  mass_face += 1.0 / lumped_mass(dof);
393 
394  mass_face /= face_dofs.size();
395 
396  return mass_face;
397 }
const unsigned int & neighborSide() const
MooseVariable & _primary_var
const Elem *const & _current_primary
virtual void getDofIndices(const Elem *elem, std::vector< dof_id_type > &dof_indices) const override
Assembly & _assembly
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real

◆ computeJacobian()

virtual void ExplicitDynamicsContactConstraint::computeJacobian ( )
inlineoverridevirtual

Computes the jacobian for the current element.

Reimplemented from NodeFaceConstraint.

Definition at line 44 of file ExplicitDynamicsContactConstraint.h.

44 {}

◆ computeOffDiagJacobian()

virtual void ExplicitDynamicsContactConstraint::computeOffDiagJacobian ( unsigned int  )
inlineoverridevirtual

Compute off-diagonal Jacobian entries.

Parameters
jvarThe index of the coupled variable

Reimplemented from NodeFaceConstraint.

Definition at line 50 of file ExplicitDynamicsContactConstraint.h.

50 {}

◆ computeQpJacobian()

virtual Real ExplicitDynamicsContactConstraint::computeQpJacobian ( Moose::ConstraintJacobianType  )
inlineoverridevirtual

Implements NodeFaceConstraint.

Definition at line 52 of file ExplicitDynamicsContactConstraint.h.

52 { return 0.0; }

◆ computeQpOffDiagJacobian()

virtual Real ExplicitDynamicsContactConstraint::computeQpOffDiagJacobian ( Moose::ConstraintJacobianType  ,
unsigned int   
)
inlineoverridevirtual

Compute off-diagonal Jacobian entries.

Parameters
typeThe type of coupling
jvarThe index of the coupled variable

Reimplemented from NodeFaceConstraint.

Definition at line 59 of file ExplicitDynamicsContactConstraint.h.

61  {
62  return 0.0;
63  }

◆ computeQpResidual()

Real ExplicitDynamicsContactConstraint::computeQpResidual ( Moose::ConstraintType  type)
overridevirtual

Implements NodeFaceConstraint.

Definition at line 338 of file ExplicitDynamicsContactConstraint.C.

339 {
340  // We use kinematic contact. But adding the residual helps avoid element inversion.
342  Real resid = pinfo->_contact_force(_component);
343 
344  switch (type)
345  {
346  case Moose::Secondary:
347  return _test_secondary[_i][_qp] * resid;
348 
349  case Moose::Primary:
350  return _test_primary[_i][_qp] * -resid;
351  }
352 
353  return 0.0;
354 }
unsigned int _i
const VariableTestValue & _test_primary
std::map< dof_id_type, PenetrationInfo *> & _penetration_info
const Node *const & _current_node
const std::string & type() const
RealVectorValue _contact_force
VariableTestValue _test_secondary
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
PenetrationLocator & _penetration_locator
unsigned int _qp

◆ computeQpSecondaryValue()

Real ExplicitDynamicsContactConstraint::computeQpSecondaryValue ( )
overridevirtual

Implements NodeFaceConstraint.

Definition at line 331 of file ExplicitDynamicsContactConstraint.C.

332 {
333  // Not used in current implementation.
334  return _u_secondary[_qp];
335 }
const VariableValue & _u_secondary
unsigned int _qp

◆ gapOffset()

Real ExplicitDynamicsContactConstraint::gapOffset ( const Node &  node)
protected

Definition at line 357 of file ExplicitDynamicsContactConstraint.C.

Referenced by computeContactForce().

358 {
359  Real val = 0;
360 
363 
366 
367  return val;
368 }
OutputData getNodalValue(const Node &node) const
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real

◆ getMatPropDependencies()

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

Reimplemented from NodeFaceConstraint.

Definition at line 141 of file ExplicitDynamicsContactConstraint.h.

142 {
144 }
virtual const std::unordered_set< unsigned int > & getMatPropDependencies() const

◆ getPenalty()

Real ExplicitDynamicsContactConstraint::getPenalty ( const Node &  node)
protected

Definition at line 371 of file ExplicitDynamicsContactConstraint.C.

Referenced by computeContactForce().

372 {
373  // TODO: Include normalized penalty values.
374  return _penalty;
375 }

◆ isExplicitConstraint()

virtual bool ExplicitDynamicsContactConstraint::isExplicitConstraint ( ) const
inlineoverridevirtual

Reimplemented from NodeFaceConstraint.

Definition at line 67 of file ExplicitDynamicsContactConstraint.h.

67 { return true; }

◆ jacobianSetup()

virtual void ExplicitDynamicsContactConstraint::jacobianSetup ( )
inlineoverridevirtual

Reimplemented from NodeFaceConstraint.

Definition at line 34 of file ExplicitDynamicsContactConstraint.h.

34 {}

◆ overwriteBoundaryVariables()

virtual void ExplicitDynamicsContactConstraint::overwriteBoundaryVariables ( NumericVector< Number > &  ,
const Node &   
) const
inlineoverridevirtual

Reimplemented from NodeFaceConstraint.

Definition at line 77 of file ExplicitDynamicsContactConstraint.h.

78  {
79  return;
80  };

◆ residualEnd()

virtual void ExplicitDynamicsContactConstraint::residualEnd ( )
inlineoverridevirtual

Reimplemented from NodeFaceConstraint.

Definition at line 35 of file ExplicitDynamicsContactConstraint.h.

35 {}

◆ shouldApply()

bool ExplicitDynamicsContactConstraint::shouldApply ( )
overridevirtual

Reimplemented from NodeFaceConstraint.

Definition at line 182 of file ExplicitDynamicsContactConstraint.C.

183 {
184  if (_current_node->processor_id() != _fe_problem.processor_id())
185  return false;
186 
187  bool in_contact = false;
188 
189  std::map<dof_id_type, PenetrationInfo *>::iterator found =
191  if (found != _penetration_locator._penetration_info.end())
192  {
193  PenetrationInfo * pinfo = found->second;
194  if (pinfo != nullptr)
195  {
196  // This computes the contact force once per constraint, rather than once per quad point
197  // and for both primary and secondary cases.
198  if (_component == 0)
199  computeContactForce(*_current_node, pinfo, true);
200 
201  if (pinfo->isCaptured())
202  in_contact = true;
203  }
204  }
205 
206  return in_contact;
207 }
std::map< dof_id_type, PenetrationInfo *> & _penetration_info
const Node *const & _current_node
void computeContactForce(const Node &node, PenetrationInfo *pinfo, bool update_contact_set)
bool isCaptured() const
FEProblemBase & _fe_problem
PenetrationLocator & _penetration_locator
processor_id_type processor_id() const

◆ solveImpactEquations()

void ExplicitDynamicsContactConstraint::solveImpactEquations ( const Node &  node,
PenetrationInfo pinfo,
const RealVectorValue distance_gap 
)
protected

Determine "Lagrange multipliers" from the iterative solution of the impact problem.

Parameters
nodeThe number of the variable to be checked
pinfoThe component index computed in this routine
distance_gapThe gap distance at the constraint node

Definition at line 275 of file ExplicitDynamicsContactConstraint.C.

Referenced by computeContactForce().

278 {
279  // Momentum balance, uncoupled normal pressure
280  // See Heinstein et al, 2000, Contact-impact modeling in explicit transient dynamics.
281 
282  dof_id_type dof_x = node.dof_number(_sys.number(), _var_objects[0]->number(), 0);
283  dof_id_type dof_y = node.dof_number(_sys.number(), _var_objects[1]->number(), 0);
284  dof_id_type dof_z =
285  _mesh_dimension == 3 ? node.dof_number(_sys.number(), _var_objects[2]->number(), 0) : 0;
286 
287  auto & u_dot = *_sys.solutionUDot();
288 
289  // Get lumped mass value
290  const auto & diag = _sys.getVector("mass_matrix_diag_inverted");
291 
292  Real mass_node = 1.0 / diag(dof_x);
293  Real mass_face = computeFaceMass(diag);
294 
295  Real mass_eff = (mass_face * mass_node) / (mass_face + mass_node);
296 
297  // Include effects of other forces:
298  // Initial guess: v_{n-1/2} + dt * M^{-1} * (F^{ext} - F^{int})
299  Real velocity_x = u_dot(dof_x) + _dt / mass_node * -1 * _residual_copy(dof_x);
300  Real velocity_y = u_dot(dof_y) + _dt / mass_node * -1 * _residual_copy(dof_y);
301  Real velocity_z = u_dot(dof_z) + _dt / mass_node * -1 * _residual_copy(dof_z);
302 
303  Real n_velocity_x = _neighbor_vel_x[0];
304  Real n_velocity_y = _neighbor_vel_y[0];
305  Real n_velocity_z = _neighbor_vel_z[0];
306 
307  RealVectorValue secondary_velocity(
308  velocity_x, velocity_y, _mesh.dimension() == 3 ? velocity_z : 0.0);
309  RealVectorValue closest_point_velocity(
310  n_velocity_x, n_velocity_y, _mesh.dimension() == 3 ? n_velocity_z : 0.0);
311  // Compute initial gap rate
312  Real gap_rate = pinfo->_normal * (secondary_velocity - closest_point_velocity);
313 
314  // Compute the force increment needed to set gap rate to 0
315  RealVectorValue impulse_force = pinfo->_normal * (gap_rate * mass_eff) / _dt;
316  pinfo->_contact_force = impulse_force;
317 
318  // recalculate velocity to determine gap rate
319  velocity_x -= _dt / mass_eff * impulse_force(0);
320  velocity_y -= _dt / mass_eff * impulse_force(1);
321  velocity_z -= _dt / mass_eff * impulse_force(2);
322 
323  // Recalculate gap rate for backwards compatibility
324  secondary_velocity = {velocity_x, velocity_y, _mesh.dimension() == 3 ? velocity_z : 0.0};
325  gap_rate = pinfo->_normal * (secondary_velocity - closest_point_velocity);
326  // gap rate is now always near "0"
327  _gap_rate->setNodalValue(gap_rate);
328 }
MooseMesh & _mesh
RealVectorValue _normal
Real computeFaceMass(const NumericVector< Real > &lumped_mass)
virtual void setNodalValue(const Real &value, unsigned int idx=0)=0
static const std::string velocity_z
Definition: NS.h:48
static const std::string velocity_x
Definition: NS.h:46
SystemBase & _sys
virtual unsigned int dimension() const
virtual NumericVector< Number > * solutionUDot()
static const std::string velocity_y
Definition: NS.h:47
unsigned int number() const
RealVectorValue _contact_force
MooseWritableVariable * _gap_rate
Nodal gap rate (output for debugging or analyst perusal)
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
const VariableValue & _neighbor_vel_z
Z component of velocity at the closest point.
const VariableValue & _neighbor_vel_y
Y component of velocity at the closest point.
const VariableValue & _neighbor_vel_x
X component of velocity at the closest point.
virtual NumericVector< Number > & getVector(const std::string &name)
uint8_t dof_id_type

◆ timestepSetup()

void ExplicitDynamicsContactConstraint::timestepSetup ( )
overridevirtual

Reimplemented from NodeFaceConstraint.

Definition at line 132 of file ExplicitDynamicsContactConstraint.C.

133 {
134  if (_component == 0)
135  {
136  updateContactStatefulData(/* beginning_of_step = */ true);
137  _update_stateful_data = false;
138  _dof_to_position.clear();
139  }
140 }
std::unordered_map< dof_id_type, Real > _dof_to_position
virtual void updateContactStatefulData(bool beginning_of_step)

◆ updateContactStatefulData()

void ExplicitDynamicsContactConstraint::updateContactStatefulData ( bool  beginning_of_step)
virtual

Definition at line 143 of file ExplicitDynamicsContactConstraint.C.

Referenced by timestepSetup().

144 {
145  for (auto & [secondary_node_num, pinfo] : _penetration_locator._penetration_info)
146  {
147  if (!pinfo)
148  continue;
149 
150  const Node & node = _mesh.nodeRef(secondary_node_num);
151  if (node.n_comp(_sys.number(), _vars[_component]) < 1)
152  continue;
153 
154  if (beginning_of_step)
155  {
157  {
158  pinfo->_contact_force_old = pinfo->_contact_force;
159  pinfo->_accumulated_slip_old = pinfo->_accumulated_slip;
160  pinfo->_frictional_energy_old = pinfo->_frictional_energy;
161  pinfo->_mech_status_old = pinfo->_mech_status;
162  }
163  else if (pinfo->_mech_status_old == PenetrationInfo::MS_NO_CONTACT &&
164  pinfo->_mech_status != PenetrationInfo::MS_NO_CONTACT)
165  {
166  mooseWarning("Previous step did not converge. Check results");
167  // The penetration info object could be based on a bad state so delete it
168  delete pinfo;
169  pinfo = nullptr;
170  continue;
171  }
172 
173  pinfo->_starting_elem = pinfo->_elem;
174  pinfo->_starting_side_num = pinfo->_side_num;
175  pinfo->_starting_closest_point_ref = pinfo->_closest_point_ref;
176  }
177  pinfo->_incremental_slip_prev_iter = pinfo->_incremental_slip;
178  }
179 }
MooseMesh & _mesh
std::map< dof_id_type, PenetrationInfo *> & _penetration_info
virtual const Node & nodeRef(const dof_id_type i) const
void mooseWarning(Args &&... args) const
SystemBase & _sys
unsigned int number() const
Executioner * getExecutioner() const
virtual bool lastSolveConverged() const=0
PenetrationLocator & _penetration_locator

◆ validParams()

InputParameters ExplicitDynamicsContactConstraint::validParams ( )
static

Definition at line 35 of file ExplicitDynamicsContactConstraint.C.

36 {
40 
41  params.addRequiredParam<BoundaryName>("boundary", "The primary boundary");
42  params.addParam<BoundaryName>("secondary", "The secondary boundary");
43  params.addRequiredParam<unsigned int>("component",
44  "An integer corresponding to the direction "
45  "the variable this constraint acts on. (0 for x, "
46  "1 for y, 2 for z)");
47  params.addCoupledVar(
48  "displacements",
49  "The displacements appropriate for the simulation geometry and coordinate system");
50  params.addCoupledVar("secondary_gap_offset", "offset to the gap distance from secondary side");
51  params.addCoupledVar("mapped_primary_gap_offset",
52  "offset to the gap distance mapped from primary side");
53  params.set<bool>("use_displaced_mesh") = true;
54  params.addParam<Real>("penalty",
55  1e8,
56  "The penalty to apply. Its optimal value can vary depending on the "
57  "stiffness of the materials");
58  params.addParam<Real>("friction_coefficient", 0, "The friction coefficient");
59  params.addParam<Real>("tangential_tolerance",
60  "Tangential distance to extend edges of contact surfaces");
61  params.addParam<bool>(
62  "print_contact_nodes", false, "Whether to print the number of nodes in contact.");
63 
64  params.addClassDescription(
65  "Apply non-penetration constraints on the mechanical deformation in explicit dynamics "
66  "using a node on face formulation by solving uncoupled momentum-balance equations.");
67  return params;
68 }
void addParam(const std::string &name, const std::initializer_list< typename T::value_type > &value, const std::string &doc_string)
T & set(const std::string &name, bool quiet_mode=false)
static InputParameters validParams()
void addRequiredParam(const std::string &name, const std::string &doc_string)
static InputParameters commonParameters()
Define parameters used by multiple contact objects.
void addCoupledVar(const std::string &name, const std::string &doc_string)
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
static InputParameters validParams()
void addClassDescription(const std::string &doc_string)

Member Data Documentation

◆ _component

const unsigned int ExplicitDynamicsContactConstraint::_component
protected

◆ _displaced_problem

MooseSharedPointer<DisplacedProblem> ExplicitDynamicsContactConstraint::_displaced_problem
protected

Definition at line 93 of file ExplicitDynamicsContactConstraint.h.

◆ _dof_to_position

std::unordered_map<dof_id_type, Real> ExplicitDynamicsContactConstraint::_dof_to_position
private

Definition at line 131 of file ExplicitDynamicsContactConstraint.h.

Referenced by timestepSetup().

◆ _gap_rate

MooseWritableVariable* ExplicitDynamicsContactConstraint::_gap_rate
protected

Nodal gap rate (output for debugging or analyst perusal)

Definition at line 121 of file ExplicitDynamicsContactConstraint.h.

Referenced by solveImpactEquations().

◆ _has_mapped_primary_gap_offset

const bool ExplicitDynamicsContactConstraint::_has_mapped_primary_gap_offset
protected

Definition at line 109 of file ExplicitDynamicsContactConstraint.h.

Referenced by gapOffset().

◆ _has_secondary_gap_offset

const bool ExplicitDynamicsContactConstraint::_has_secondary_gap_offset
protected

Definition at line 107 of file ExplicitDynamicsContactConstraint.h.

Referenced by gapOffset().

◆ _mapped_primary_gap_offset_var

const MooseVariable* const ExplicitDynamicsContactConstraint::_mapped_primary_gap_offset_var
protected

Definition at line 110 of file ExplicitDynamicsContactConstraint.h.

Referenced by gapOffset().

◆ _mesh_dimension

const unsigned int ExplicitDynamicsContactConstraint::_mesh_dimension
protected

Definition at line 102 of file ExplicitDynamicsContactConstraint.h.

Referenced by solveImpactEquations().

◆ _model

const ExplicitDynamicsContactModel ExplicitDynamicsContactConstraint::_model
protected

◆ _neighbor_vel_x

const VariableValue& ExplicitDynamicsContactConstraint::_neighbor_vel_x
protected

X component of velocity at the closest point.

Definition at line 124 of file ExplicitDynamicsContactConstraint.h.

Referenced by solveImpactEquations().

◆ _neighbor_vel_y

const VariableValue& ExplicitDynamicsContactConstraint::_neighbor_vel_y
protected

Y component of velocity at the closest point.

Definition at line 126 of file ExplicitDynamicsContactConstraint.h.

Referenced by solveImpactEquations().

◆ _neighbor_vel_z

const VariableValue& ExplicitDynamicsContactConstraint::_neighbor_vel_z
protected

Z component of velocity at the closest point.

Definition at line 128 of file ExplicitDynamicsContactConstraint.h.

Referenced by solveImpactEquations().

◆ _no_iterations

const unsigned int ExplicitDynamicsContactConstraint::_no_iterations = 0
staticprotected

Definition at line 116 of file ExplicitDynamicsContactConstraint.h.

◆ _penalty

const Real ExplicitDynamicsContactConstraint::_penalty
protected

Definition at line 112 of file ExplicitDynamicsContactConstraint.h.

Referenced by getPenalty().

◆ _print_contact_nodes

const bool ExplicitDynamicsContactConstraint::_print_contact_nodes
protected

Definition at line 114 of file ExplicitDynamicsContactConstraint.h.

◆ _residual_copy

NumericVector<Number>& ExplicitDynamicsContactConstraint::_residual_copy
protected

Definition at line 118 of file ExplicitDynamicsContactConstraint.h.

Referenced by solveImpactEquations().

◆ _secondary_gap_offset_var

const MooseVariable* const ExplicitDynamicsContactConstraint::_secondary_gap_offset_var
protected

Definition at line 108 of file ExplicitDynamicsContactConstraint.h.

Referenced by gapOffset().

◆ _update_stateful_data

bool ExplicitDynamicsContactConstraint::_update_stateful_data
protected

Definition at line 100 of file ExplicitDynamicsContactConstraint.h.

Referenced by timestepSetup().

◆ _var_objects

std::vector<MooseVariable *> ExplicitDynamicsContactConstraint::_var_objects
protected

◆ _vars

std::vector<unsigned int> ExplicitDynamicsContactConstraint::_vars
protected

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