Calculates the inertial torque and inertia proportional damping for nodal rotational inertia. More...

#include <NodalRotationalInertia.h>

Inheritance diagram for NodalRotationalInertia:

Public Types
enum	ResidualTagType { ResidualTagType::NonReference, ResidualTagType::Reference }

enum	TEST_TYPE

typedef DataFileName	DataFileParameterType

Public Member Functions
	NodalRotationalInertia (const InputParameters &parameters)

virtual void	computeJacobian () override

virtual void	computeOffDiagJacobian (unsigned int jvar) override

const MooseVariable &	variable () const override

void	setSubdomains (const std::set< SubdomainID > &sub_ids)

virtual void	computeResidualAndJacobian ()

virtual void	computeOffDiagJacobianScalar (unsigned int)

virtual void	computeNonlocalJacobian ()

virtual void	computeNonlocalOffDiagJacobian (unsigned int)

const SubProblem &	subProblem () const

virtual void	prepareShapes (unsigned int var_num)

virtual bool	enabled () const

std::shared_ptr< MooseObject >	getSharedPtr ()

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

MooseApp &	getMooseApp () const

const std::string &	type () const

virtual const std::string &	name () const

std::string	typeAndName () const

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

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

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

const InputParameters &	parameters () const

MooseObjectName	uniqueName () const

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

virtual void	initialSetup ()

virtual void	timestepSetup ()

virtual void	jacobianSetup ()

virtual void	residualSetup ()

virtual void	subdomainSetup ()

virtual void	customSetup (const ExecFlagType &)

const ExecFlagEnum &	getExecuteOnEnum () const

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

const Function &	getFunctionByName (const FunctionName &name) const

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

bool	hasFunctionByName (const FunctionName &name) const

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

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

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

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

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

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 VectorPostprocessorValue &	getVectorPostprocessorValue (const std::string &param_name, const std::string &vector_name) const

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

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

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

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

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

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

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

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

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

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

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

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

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

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

bool	hasVectorPostprocessorByName (const VectorPostprocessorName &name) const

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

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

const std::vector< SubdomainName > &	blocks () const

unsigned int	numBlocks () const

virtual const std::set< SubdomainID > &	blockIDs () const

unsigned int	blocksMaxDimension () const

bool	hasBlocks (const SubdomainName &name) const

bool	hasBlocks (const std::vector< SubdomainName > &names) const

bool	hasBlocks (SubdomainID id) const

bool	hasBlocks (const std::vector< SubdomainID > &ids) const

bool	hasBlocks (const std::set< SubdomainID > &ids) const

bool	isBlockSubset (const std::set< SubdomainID > &ids) const

bool	isBlockSubset (const std::vector< SubdomainID > &ids) const

bool	hasBlockMaterialProperty (const std::string &prop_name)

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

virtual bool	blockRestricted () const

virtual void	checkVariable (const MooseVariableFieldBase &variable) const

virtual const std::set< BoundaryID > &	boundaryIDs () const

const std::vector< BoundaryName > &	boundaryNames () const

unsigned int	numBoundaryIDs () const

bool	hasBoundary (const BoundaryName &name) const

bool	hasBoundary (const std::vector< BoundaryName > &names) const

bool	hasBoundary (const BoundaryID &id) const

bool	hasBoundary (const std::vector< BoundaryID > &ids, TEST_TYPE type=ALL) const

bool	hasBoundary (const std::set< BoundaryID > &ids, TEST_TYPE type=ALL) const

bool	isBoundarySubset (const std::set< BoundaryID > &ids) const

bool	isBoundarySubset (const std::vector< BoundaryID > &ids) const

bool	hasBoundaryMaterialProperty (const std::string &prop_name) const

virtual bool	boundaryRestricted () const

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

virtual bool	checkVariableBoundaryIntegrity () const

PenetrationLocator &	getPenetrationLocator (const BoundaryName &primary, const BoundaryName &secondary, Order order)

PenetrationLocator &	getQuadraturePenetrationLocator (const BoundaryName &primary, const BoundaryName &secondary, Order order)

NearestNodeLocator &	getNearestNodeLocator (const BoundaryName &primary, const BoundaryName &secondary)

NearestNodeLocator &	getQuadratureNearestNodeLocator (const BoundaryName &primary, const BoundaryName &secondary)

bool	requiresGeometricSearch () const

virtual const VariableValue &	coupledValueByName (const std::string &var_name)

virtual const ArrayVariableValue &	coupledArrayValueByName (const std::string &var_name)

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

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

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

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

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

void	addFEVariableCoupleableVectorTag (TagID tag)

void	addFEVariableCoupleableMatrixTag (TagID tag)

std::set< TagID > &	getFEVariableCoupleableVectorTags ()

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

std::set< TagID > &	getFEVariableCoupleableMatrixTags ()

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

auto &	getWritableCoupledVariables () const

bool	hasWritableCoupledVariables () const

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

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

const ADVariableGradient &	getADDefaultGradient () const

const ADVectorVariableGradient &	getADDefaultVectorGradient () const

const ADVariableSecond &	getADDefaultSecond () const

const ADVectorVariableCurl &	getADDefaultCurl () const

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)

MooseVariableBase *	mooseVariableBase () 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 PostprocessorValue &	getPostprocessorValue (const std::string &param_name, const unsigned int index=0) const

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

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

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

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

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

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

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

const PostprocessorValue &	getPostprocessorValueOldByName (const PostprocessorName &name) const

const PostprocessorValue &	getPostprocessorValueOldByName (const PostprocessorName &name) const

const PostprocessorValue &	getPostprocessorValueOlderByName (const PostprocessorName &name) const

const PostprocessorValue &	getPostprocessorValueOlderByName (const PostprocessorName &name) const

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

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

bool	isVectorPostprocessorDistributedByName (const VectorPostprocessorName &name) const

bool	isVectorPostprocessorDistributedByName (const VectorPostprocessorName &name) const

const Parallel::Communicator &	comm () const

processor_id_type	n_processors () const

processor_id_type	processor_id () const

Static Public Member Functions
static InputParameters	validParams ()

static bool	restricted (const std::set< BoundaryID > &ids)

Public Attributes
	ALL

	ANY

const ConsoleStream	_console

Protected Member Functions
virtual Real	computeQpResidual () override

virtual Real	computeQpJacobian () override

virtual Real	computeQpOffDiagJacobian (unsigned int jvar) override

virtual void	computeResidual () override

Moose::NodeArg	nodeArg () const

virtual void	precalculateResidual ()

virtual void	precalculateJacobian ()

virtual void	precalculateOffDiagJacobian (unsigned int)

const MooseVariableFieldBase &	getVariable (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	addResidualsAndJacobian (Assembly &assembly, const Residuals &residuals, const Indices &dof_indices, Real scaling_factor)

void	addJacobian (Assembly &assembly, const Residuals &residuals, const Indices &dof_indices, Real scaling_factor)

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 bool	hasBlockMaterialPropertyHelper (const std::string &prop_name)

void	initializeBlockRestrictable (const MooseObject *moose_object)

Moose::CoordinateSystemType	getBlockCoordSystem ()

bool	hasBoundaryMaterialPropertyHelper (const std::string &prop_name) const

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

void	checkWritableVar (MooseWritableVariable *var)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

const ADVariableValue &	adZeroValue () const

const ADVariableGradient &	adZeroGradient () const

const ADVariableSecond &	adZeroSecond () const

const GenericVariableValue< is_ad > &	genericZeroValue ()

const GenericVariableValue< false > &	genericZeroValue ()

const GenericVariableValue< true > &	genericZeroValue ()

const GenericVariableGradient< is_ad > &	genericZeroGradient ()

const GenericVariableGradient< false > &	genericZeroGradient ()

const GenericVariableGradient< true > &	genericZeroGradient ()

const GenericVariableSecond< is_ad > &	genericZeroSecond ()

const GenericVariableSecond< false > &	genericZeroSecond ()

const GenericVariableSecond< true > &	genericZeroSecond ()

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

void	addMooseVariableDependency (MooseVariableFieldBase *var)

void	addMooseVariableDependency (const std::vector< MooseVariableFieldBase * > &vars)

virtual const OutputTools< Real >::VariableValue &	value ()

virtual const OutputTools< Real >::VariableValue &	valueOld ()

virtual const OutputTools< Real >::VariableValue &	valueOlder ()

virtual const OutputTools< Real >::VariableValue &	dot ()

virtual const OutputTools< Real >::VariableValue &	dotDot ()

virtual const OutputTools< Real >::VariableValue &	dotOld ()

virtual const OutputTools< Real >::VariableValue &	dotDotOld ()

virtual const VariableValue &	dotDu ()

virtual const VariableValue &	dotDotDu ()

virtual const OutputTools< Real >::VariableGradient &	gradient ()

virtual const OutputTools< Real >::VariableGradient &	gradientOld ()

virtual const OutputTools< Real >::VariableGradient &	gradientOlder ()

virtual const OutputTools< Real >::VariableSecond &	second ()

virtual const OutputTools< Real >::VariableSecond &	secondOld ()

virtual const OutputTools< Real >::VariableSecond &	secondOlder ()

virtual const OutputTools< Real >::VariableTestSecond &	secondTest ()

virtual const OutputTools< Real >::VariableTestSecond &	secondTestFace ()

virtual const OutputTools< Real >::VariablePhiSecond &	secondPhi ()

virtual const OutputTools< Real >::VariablePhiSecond &	secondPhiFace ()

Protected Attributes
const bool	_has_beta
	Booleans for validity of params. More...

const bool	_has_gamma

const bool	_has_rot_velocities

const bool	_has_rot_accelerations

const bool	_has_x_orientation

const bool	_has_y_orientation

AuxiliarySystem *	_aux_sys
	Auxiliary system object. More...

unsigned int	_nrot
	Number of coupled rotational variables. More...

std::vector< const VariableValue * >	_rot
	Value of rotational displacements. More...

std::vector< const VariableValue * >	_rot_old
	Old value of rotational displacements. More...

std::vector< unsigned int >	_rot_vel_num
	Variable numbers for rotational velocity aux variables. More...

std::vector< unsigned int >	_rot_accel_num
	Variable numbers for rotational acceleration aux variables. More...

std::vector< unsigned int >	_rot_variables
	Variable numbers for rotational variables. More...

std::vector< Real >	_rot_accel
	Current acceleration of the node. More...

std::vector< Real >	_rot_vel
	Current velocity of the node. More...

std::vector< Real >	_rot_vel_old
	Old velocity of the node. More...

const Real	_beta
	Newmark time integration parameter. More...

const Real	_gamma
	Newmark time integration parameter. More...

const Real &	_eta
	Mass proportional Rayliegh damping. More...

const Real &	_alpha
	HHT time integration parameter. More...

const unsigned int	_component
	Component along which torque is applied. More...

RankTwoTensor	_inertia
	Moment of inertia tensor in global coordinate system. More...

std::vector< const VariableValue * >	_rot_dot_residual
	Rotational udot residual. More...

std::vector< const VariableValue * >	_rot_vel_old_value
	Old velocity value. More...

std::vector< const VariableValue * >	_rot_dotdot_residual
	Rotational udotdot residual. More...

const VariableValue *	_du_dot_du
	du_dot_du value More...

const VariableValue *	_du_dotdot_du
	du_dotdot_du value More...

const TimeIntegrator &	_time_integrator
	The TimeIntegrator. More...

const VariableValue &	_u_dot

const VariableValue &	_u

bool	_has_save_in

std::vector< MooseVariableFEBase *>	_save_in

std::vector< AuxVariableName >	_save_in_strings

bool	_has_diag_save_in

std::vector< MooseVariableFEBase *>	_diag_save_in

std::vector< AuxVariableName >	_diag_save_in_strings

FEProblemBase &	_fe_problem

MooseVariable &	_var

const Node *const &	_current_node

unsigned int	_qp

const std::set< SubdomainID > *	_sub_ids

SubProblem &	_subproblem

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

const MaterialData *	_blk_material_data

GeometricSearchData &	_geometric_search_data

bool	_requires_geometric_search

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

Detailed Description

Calculates the inertial torque and inertia proportional damping for nodal rotational inertia.

Definition at line 23 of file NodalRotationalInertia.h.

Constructor & Destructor Documentation

◆ NodalRotationalInertia()

NodalRotationalInertia::NodalRotationalInertia ( const InputParameters & parameters )

Definition at line 64 of file NodalRotationalInertia.C.

   : TimeNodalKernel(parameters),
     _has_beta(isParamValid("beta")),
     _has_gamma(isParamValid("gamma")),
     _has_rot_velocities(isParamValid("rotational_velocities")),
     _has_rot_accelerations(isParamValid("rotational_accelerations")),
     _has_x_orientation(isParamValid("x_orientation")),
     _has_y_orientation(isParamValid("y_orientation")),
     _nrot(coupledComponents("rotations")),
     _rot(_nrot),
     _rot_old(_nrot),
     _rot_vel_num(_nrot),
     _rot_accel_num(_nrot),
     _rot_variables(_nrot),
     _rot_accel(_nrot),
     _rot_vel(_nrot),
     _rot_vel_old(_nrot),
     _beta(_has_beta ? getParam<Real>("beta") : 0.1),
     _gamma(_has_gamma ? getParam<Real>("gamma") : 0.1),
     _eta(getParam<Real>("eta")),
     _alpha(getParam<Real>("alpha")),
     _component(getParam<unsigned int>("component")),
     _rot_dot_residual(_nrot),
     _rot_vel_old_value(_nrot),
     _rot_dotdot_residual(_nrot),
     _time_integrator(_sys.getTimeIntegrator(_var.number()))
 {
   if (_has_beta && _has_gamma && _has_rot_velocities && _has_rot_accelerations)
   {
     _aux_sys = &(_fe_problem.getAuxiliarySystem());
     if (coupledComponents("rotational_velocities") != _nrot ||
         coupledComponents("rotational_accelerations") != _nrot)
       mooseError(
           "NodalRotationalInertia: rotational_velocities and rotational_accelerations should "
           "be same size "
           "as rotations.");
 
     for (unsigned int i = 0; i < _nrot; ++i)
     {
       MooseVariable * rot_var = getVar("rotations", i);
       MooseVariable * rot_vel_var = getVar("rotational_velocities", i);
       MooseVariable * rot_accel_var = getVar("rotational_accelerations", i);
 
       _rot[i] = &rot_var->dofValues();
       _rot_old[i] = &rot_var->dofValuesOld();
 
       _rot_vel_num[i] = rot_vel_var->number();
       _rot_accel_num[i] = rot_accel_var->number();
 
       _rot_variables[i] = coupled("rotations", i);
     }
   }
   else if (!_has_beta && !_has_gamma && !_has_rot_velocities && !_has_rot_accelerations)
   {
     for (unsigned int i = 0; i < _nrot; ++i)
     {
       MooseVariable * rot_var = getVar("rotations", i);
       _rot_vel_old_value[i] = &rot_var->dofValuesDotOld();
       _rot_dot_residual[i] =
           &coupledVectorTagDofValue("rotations", _time_integrator.uDotFactorTag(), i);
       _rot_dotdot_residual[i] =
           &coupledVectorTagDofValue("rotations", _time_integrator.uDotDotFactorTag(), i);
 
       if (i == 0)
       {
         _du_dot_du = &rot_var->dofValuesDuDotDu();
         _du_dotdot_du = &rot_var->dofValuesDuDotDotDu();
       }
     }
   }
   else
     mooseError("NodalRotationalInertia: Either all or none of `beta`, `gamma`, "
                "`rotational_velocities` and `rotational_accelerations` should be provided as "
                "input.");
 
   // Store inertia values in inertia tensor
   _inertia.zero();
   _inertia(0, 0) = getParam<Real>("Ixx");
   _inertia(0, 1) = getParam<Real>("Ixy");
   _inertia(0, 2) = getParam<Real>("Ixz");
   _inertia(1, 0) = _inertia(0, 1);
   _inertia(1, 1) = getParam<Real>("Iyy");
   _inertia(1, 2) = getParam<Real>("Iyz");
   _inertia(2, 0) = _inertia(0, 2);
   _inertia(2, 1) = _inertia(1, 2);
   _inertia(2, 2) = getParam<Real>("Izz");
 
   // Check for positive definiteness of matrix using Sylvester's criterion.
   const Real det1 = _inertia(0, 0);
   const Real det2 = _inertia(0, 0) * _inertia(1, 1) - _inertia(0, 1) * _inertia(1, 0);
   const Real det3 = _inertia.det();
 
   if (det1 < 1e-6 || det2 < 1e-6 || det3 < 1e-6)
     mooseError("NodalRotationalInertia: The moment of inertia tensor should be positive definite.");
   if (_has_x_orientation && _has_y_orientation)
   {
     const RealGradient x_orientation = getParam<RealGradient>("x_orientation");
     const RealGradient y_orientation = getParam<RealGradient>("y_orientation");
 
     if (std::abs(x_orientation.norm() - 1.0) > 1e-6 || std::abs(y_orientation.norm() - 1.0) > 1e-6)
       mooseError("NodalRotationalInertia: x_orientation and y_orientation must be unit vectors.");
 
     Real sum = x_orientation(0) * y_orientation(0) + x_orientation(1) * y_orientation(1) +
                x_orientation(2) * y_orientation(2);
 
     if (std::abs(sum) > 1e-4)
       mooseError("NodalRotationalInertia: x_orientation and y_orientation should be perpendicular "
                  "to each other.");
 
     // Calculate z orientation as a cross product of the x and y orientations
     RealGradient z_orientation;
     z_orientation(0) = (x_orientation(1) * y_orientation(2) - x_orientation(2) * y_orientation(1));
     z_orientation(1) = (x_orientation(2) * y_orientation(0) - x_orientation(0) * y_orientation(2));
     z_orientation(2) = (x_orientation(0) * y_orientation(1) - x_orientation(1) * y_orientation(0));
 
     // Rotation matrix from global to original beam local configuration
     RankTwoTensor orientation;
     orientation(0, 0) = x_orientation(0);
     orientation(0, 1) = x_orientation(1);
     orientation(0, 2) = x_orientation(2);
     orientation(1, 0) = y_orientation(0);
     orientation(1, 1) = y_orientation(1);
     orientation(1, 2) = y_orientation(2);
     orientation(2, 0) = z_orientation(0);
     orientation(2, 1) = z_orientation(1);
     orientation(2, 2) = z_orientation(2);
 
     RankTwoTensor global_inertia = orientation.transpose() * _inertia * orientation;
     _inertia = global_inertia;
   }
   else if ((_has_x_orientation && !_has_y_orientation) ||
            (!_has_x_orientation && _has_y_orientation))
     mooseError("NodalRotationalInertia: Both x_orientation and y_orientation should be provided if "
                "x_orientation or "
                "y_orientation is different from global x or y direction, respectively.");
 
   // Check for Explicit and alpha parameter
   if (_alpha != 0 && _time_integrator.isExplicit())
     mooseError("NodalRotationalInertia: HHT time integration parameter can only be used with "
                "Newmark-Beta time integration.");
 
   // Check for Explicit and beta parameter
   if (_has_beta != 0 && _time_integrator.isExplicit())
     mooseError("NodalRotationalInertia: beta time integration parameter can only be used with "
                "Newmark-Beta time integrator.");
 }

Member Function Documentation

◆ computeQpJacobian()

Real NodalRotationalInertia::computeQpJacobian ( )

overrideprotectedvirtual

Reimplemented from TimeNodalKernel.

Definition at line 258 of file NodalRotationalInertia.C.

 {
   if (_dt == 0)
     return 0.0;
   else
   {
     if (_has_beta)
       return _inertia(_component, _component) / (_beta * _dt * _dt) +
              _eta * (1.0 + _alpha) * _inertia(_component, _component) * _gamma / _beta / _dt;
     else if (_time_integrator.isExplicit())
       // for explicit central difference integration, _eta does not appear in the
       // Jacobian (mass matrix), and alpha is zero
       return _inertia(_component, _component) * (*_du_dotdot_du)[_qp];
     else
       // for NewmarkBeta time integrator
       return _inertia(_component, _component) * (*_du_dotdot_du)[_qp] +
              _eta * (1.0 + _alpha) * _inertia(_component, _component) * (*_du_dot_du)[_qp];
   }
 }

◆ computeQpOffDiagJacobian()

Real NodalRotationalInertia::computeQpOffDiagJacobian ( unsigned int jvar )

overrideprotectedvirtual

Reimplemented from TimeNodalKernel.

Definition at line 279 of file NodalRotationalInertia.C.

 {
   unsigned int coupled_component = 0;
   bool rot_coupled = false;
 
   for (unsigned int i = 0; i < _nrot; ++i)
   {
     if (jvar == _rot_variables[i])
     {
       coupled_component = i;
       rot_coupled = true;
       break;
     }
   }
 
   if (_dt == 0)
     return 0.0;
   else if (rot_coupled)
   {
     if (_has_beta)
       return _inertia(_component, coupled_component) / (_beta * _dt * _dt) +
              _eta * (1.0 + _alpha) * _inertia(_component, coupled_component) * _gamma / _beta / _dt;
     else
       return _inertia(_component, coupled_component) * (*_du_dotdot_du)[_qp] +
              _eta * (1.0 + _alpha) * _inertia(_component, coupled_component) * (*_du_dot_du)[_qp];
   }
   else
     return 0.0;
 }

◆ computeQpResidual()

Real NodalRotationalInertia::computeQpResidual ( )

overrideprotectedvirtual

Implements TimeNodalKernel.

Definition at line 212 of file NodalRotationalInertia.C.

 {
   if (_dt == 0)
     return 0;
   else
   {
     if (_has_beta)
     {
       const NumericVector<Number> & aux_sol_old = _aux_sys->solutionOld();
 
       for (unsigned int i = 0; i < _nrot; ++i)
       {
         _rot_vel_old[i] =
             aux_sol_old(_current_node->dof_number(_aux_sys->number(), _rot_vel_num[i], 0));
         const Real rot_accel_old =
             aux_sol_old(_current_node->dof_number(_aux_sys->number(), _rot_accel_num[i], 0));
 
         _rot_accel[i] = 1.0 / _beta *
                         ((((*_rot[i])[_qp] - (*_rot_old[i])[_qp]) / (_dt * _dt)) -
                          _rot_vel_old[i] / _dt - rot_accel_old * (0.5 - _beta));
         _rot_vel[i] =
             _rot_vel_old[i] + (_dt * (1.0 - _gamma)) * rot_accel_old + _gamma * _dt * _rot_accel[i];
       }
 
       Real res = 0.0;
       for (unsigned int i = 0; i < _nrot; ++i)
         res += _inertia(_component, i) * (_rot_accel[i] + _rot_vel[i] * _eta * (1.0 + _alpha) -
                                           _alpha * _eta * _rot_vel_old[i]);
 
       return res;
     }
     else
     {
       // All cases (Explicit, implicit and implicit with HHT)
       // Note that _alpha is ensured to be zero with explicit integration
       Real res = 0.0;
       for (unsigned int i = 0; i < _nrot; ++i)
         res += _inertia(_component, i) * ((*_rot_dotdot_residual[i])[_qp] +
                                           (*_rot_dot_residual[i])[_qp] * _eta * (1.0 + _alpha) -
                                           _alpha * _eta * (*_rot_vel_old_value[i])[_qp]);
       return res;
     }
   }
 }

◆ validParams()

InputParameters NodalRotationalInertia::validParams ( )

static

Definition at line 20 of file NodalRotationalInertia.C.

 {
   InputParameters params = TimeNodalKernel::validParams();
   params.addClassDescription("Calculates the inertial torques and inertia proportional damping "
                              "corresponding to the nodal rotational inertia.");
   params.addCoupledVar("rotations", "rotational displacement variables");
   params.addCoupledVar("rotational_velocities", "rotational velocity variables");
   params.addCoupledVar("rotational_accelerations", "rotational acceleration variables");
   params.addRangeCheckedParam<Real>(
       "beta", "beta>0.0", "beta parameter for Newmark Time integration");
   params.addRangeCheckedParam<Real>(
       "gamma", "gamma>0.0", "gamma parameter for Newmark Time integration");
   params.addRangeCheckedParam<Real>("eta",
                                     0.0,
                                     "eta>=0.0",
                                     "Constant real number defining the eta parameter for"
                                     "Rayleigh damping.");
   params.addRangeCheckedParam<Real>("alpha",
                                     0.0,
                                     "alpha >= -0.3333 & alpha <= 0.0",
                                     "Alpha parameter for mass dependent numerical damping induced "
                                     "by HHT time integration scheme");
   params.addRequiredRangeCheckedParam<Real>(
       "Ixx", "Ixx>0.0", "Moment of inertia in the x direction.");
   params.addRequiredRangeCheckedParam<Real>(
       "Iyy", "Iyy>0.0", "Moment of inertia in the y direction.");
   params.addRequiredRangeCheckedParam<Real>(
       "Izz", "Izz>0.0", "Moment of inertia in the z direction.");
   params.addParam<Real>("Ixy", 0.0, "Moment of inertia in the xy direction.");
   params.addParam<Real>("Ixz", 0.0, "Moment of inertia in the xz direction.");
   params.addParam<Real>("Iyz", 0.0, "Moment of inertia in the yz direction.");
   params.addParam<RealGradient>(
       "x_orientation", "Unit vector along the x direction if different from global x direction.");
   params.addParam<RealGradient>(
       "y_orientation", "Unit vector along the y direction if different from global y direction.");
   params.addRequiredRangeCheckedParam<unsigned int>(
       "component",
       "component<3",
       "An integer corresponding to the direction "
       "the variable this nodalkernel acts in. (0 for rot_x, "
       "1 for rot_y, and 2 for rot_z).");
   return params;
 }

Member Data Documentation

◆ _alpha

const Real& NodalRotationalInertia::_alpha

protected

HHT time integration parameter.

Definition at line 85 of file NodalRotationalInertia.h.

Referenced by computeQpJacobian(), computeQpOffDiagJacobian(), computeQpResidual(), and NodalRotationalInertia().

◆ _aux_sys

AuxiliarySystem* NodalRotationalInertia::_aux_sys

protected

Auxiliary system object.

Definition at line 46 of file NodalRotationalInertia.h.

Referenced by computeQpResidual(), and NodalRotationalInertia().

◆ _beta

const Real NodalRotationalInertia::_beta

protected

Newmark time integration parameter.

Definition at line 76 of file NodalRotationalInertia.h.

Referenced by computeQpJacobian(), computeQpOffDiagJacobian(), and computeQpResidual().

◆ _component

const unsigned int NodalRotationalInertia::_component

protected

Component along which torque is applied.

Definition at line 88 of file NodalRotationalInertia.h.

Referenced by computeQpJacobian(), computeQpOffDiagJacobian(), and computeQpResidual().

◆ _du_dot_du

const VariableValue* NodalRotationalInertia::_du_dot_du

protected

du_dot_du value

Definition at line 103 of file NodalRotationalInertia.h.

Referenced by NodalRotationalInertia().

◆ _du_dotdot_du

const VariableValue* NodalRotationalInertia::_du_dotdot_du

protected

du_dotdot_du value

Definition at line 106 of file NodalRotationalInertia.h.

Referenced by NodalRotationalInertia().

◆ _eta

const Real& NodalRotationalInertia::_eta

protected

Mass proportional Rayliegh damping.

Definition at line 82 of file NodalRotationalInertia.h.

Referenced by computeQpJacobian(), computeQpOffDiagJacobian(), and computeQpResidual().

◆ _gamma

const Real NodalRotationalInertia::_gamma

protected

Newmark time integration parameter.

Definition at line 79 of file NodalRotationalInertia.h.

Referenced by computeQpJacobian(), computeQpOffDiagJacobian(), and computeQpResidual().

◆ _has_beta

const bool NodalRotationalInertia::_has_beta

protected

Booleans for validity of params.

Definition at line 38 of file NodalRotationalInertia.h.

Referenced by computeQpJacobian(), computeQpOffDiagJacobian(), computeQpResidual(), and NodalRotationalInertia().

◆ _has_gamma

const bool NodalRotationalInertia::_has_gamma

protected

Definition at line 39 of file NodalRotationalInertia.h.

Referenced by NodalRotationalInertia().

◆ _has_rot_accelerations

const bool NodalRotationalInertia::_has_rot_accelerations

protected

Definition at line 41 of file NodalRotationalInertia.h.

Referenced by NodalRotationalInertia().

◆ _has_rot_velocities

const bool NodalRotationalInertia::_has_rot_velocities

protected

Definition at line 40 of file NodalRotationalInertia.h.

Referenced by NodalRotationalInertia().

◆ _has_x_orientation

const bool NodalRotationalInertia::_has_x_orientation

protected

Definition at line 42 of file NodalRotationalInertia.h.

Referenced by NodalRotationalInertia().

◆ _has_y_orientation

const bool NodalRotationalInertia::_has_y_orientation

protected

Definition at line 43 of file NodalRotationalInertia.h.

Referenced by NodalRotationalInertia().

◆ _inertia

RankTwoTensor NodalRotationalInertia::_inertia

protected

Moment of inertia tensor in global coordinate system.

Definition at line 91 of file NodalRotationalInertia.h.

Referenced by computeQpJacobian(), computeQpOffDiagJacobian(), computeQpResidual(), and NodalRotationalInertia().

◆ _nrot

unsigned int NodalRotationalInertia::_nrot

protected

Number of coupled rotational variables.

Definition at line 49 of file NodalRotationalInertia.h.

Referenced by computeQpOffDiagJacobian(), computeQpResidual(), and NodalRotationalInertia().

◆ _rot

std::vector<const VariableValue *> NodalRotationalInertia::_rot

protected

Value of rotational displacements.

Definition at line 52 of file NodalRotationalInertia.h.

Referenced by computeQpResidual(), and NodalRotationalInertia().

◆ _rot_accel

std::vector<Real> NodalRotationalInertia::_rot_accel

protected

Current acceleration of the node.

Definition at line 67 of file NodalRotationalInertia.h.

Referenced by computeQpResidual().

◆ _rot_accel_num

std::vector<unsigned int> NodalRotationalInertia::_rot_accel_num

protected

Variable numbers for rotational acceleration aux variables.

Definition at line 61 of file NodalRotationalInertia.h.

Referenced by computeQpResidual(), and NodalRotationalInertia().

◆ _rot_dot_residual

std::vector<const VariableValue *> NodalRotationalInertia::_rot_dot_residual

protected

Rotational udot residual.

Definition at line 94 of file NodalRotationalInertia.h.

Referenced by computeQpResidual(), and NodalRotationalInertia().

◆ _rot_dotdot_residual

std::vector<const VariableValue *> NodalRotationalInertia::_rot_dotdot_residual

protected

Rotational udotdot residual.

Definition at line 100 of file NodalRotationalInertia.h.

Referenced by computeQpResidual(), and NodalRotationalInertia().

◆ _rot_old

std::vector<const VariableValue *> NodalRotationalInertia::_rot_old

protected

Old value of rotational displacements.

Definition at line 55 of file NodalRotationalInertia.h.

Referenced by computeQpResidual(), and NodalRotationalInertia().

◆ _rot_variables

std::vector<unsigned int> NodalRotationalInertia::_rot_variables

protected

Variable numbers for rotational variables.

Definition at line 64 of file NodalRotationalInertia.h.

Referenced by computeQpOffDiagJacobian(), and NodalRotationalInertia().

◆ _rot_vel

std::vector<Real> NodalRotationalInertia::_rot_vel

protected

Current velocity of the node.

Definition at line 70 of file NodalRotationalInertia.h.

Referenced by computeQpResidual().

◆ _rot_vel_num

std::vector<unsigned int> NodalRotationalInertia::_rot_vel_num

protected

Variable numbers for rotational velocity aux variables.

Definition at line 58 of file NodalRotationalInertia.h.

Referenced by computeQpResidual(), and NodalRotationalInertia().

◆ _rot_vel_old

std::vector<Real> NodalRotationalInertia::_rot_vel_old

protected

Old velocity of the node.

Definition at line 73 of file NodalRotationalInertia.h.

Referenced by computeQpResidual().

◆ _rot_vel_old_value

std::vector<const VariableValue *> NodalRotationalInertia::_rot_vel_old_value

protected

Old velocity value.

Definition at line 97 of file NodalRotationalInertia.h.

Referenced by computeQpResidual(), and NodalRotationalInertia().

◆ _time_integrator

const TimeIntegrator& NodalRotationalInertia::_time_integrator

protected

The TimeIntegrator.

Definition at line 109 of file NodalRotationalInertia.h.

Referenced by computeQpJacobian(), and NodalRotationalInertia().

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

solid_mechanics/include/nodalkernels/NodalRotationalInertia.h
solid_mechanics/src/nodalkernels/NodalRotationalInertia.C

Public Types

Public Member Functions

Static Public Member Functions

Public Attributes

Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ NodalRotationalInertia()

Member Function Documentation

◆ computeQpJacobian()

◆ computeQpOffDiagJacobian()

◆ computeQpResidual()

◆ validParams()

Member Data Documentation

◆ _alpha

◆ _aux_sys

◆ _beta

◆ _component

◆ _du_dot_du

◆ _du_dotdot_du

◆ _eta

◆ _gamma

◆ _has_beta

◆ _has_gamma

◆ _has_rot_accelerations

◆ _has_rot_velocities

◆ _has_x_orientation

◆ _has_y_orientation

◆ _inertia

◆ _nrot

◆ _rot

◆ _rot_accel

◆ _rot_accel_num

◆ _rot_dot_residual

◆ _rot_dotdot_residual

◆ _rot_old

◆ _rot_variables

◆ _rot_vel

◆ _rot_vel_num

◆ _rot_vel_old

◆ _rot_vel_old_value

◆ _time_integrator