This vectorpostprocessor computes the Interaction Integral, which is used to compute various fracture mechanics parameters at a crack tip, including KI, KII, KIII, and the T stress. More...

#include <InteractionIntegral.h>

Inheritance diagram for InteractionIntegralTempl< is_ad >:

[legend]

Public Types
typedef DataFileName	DataFileParameterType

Public Member Functions
	InteractionIntegralTempl (const InputParameters &parameters)

virtual void	initialSetup () override

virtual void	initialize () override

virtual void	execute () override

virtual void	finalize () override

virtual void	threadJoin (const UserObject &y) override

SubProblem &	getSubProblem () const

bool	shouldDuplicateInitialExecution () const

virtual Real	spatialValue (const Point &) const

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

void	gatherSum (T &value)

void	gatherMax (T &value)

void	gatherMin (T &value)

void	gatherProxyValueMax (T1 &proxy, T2 &value)

void	gatherProxyValueMin (T1 &proxy, T2 &value)

void	setPrimaryThreadCopy (UserObject *primary)

UserObject *	primaryThreadCopy ()

std::set< UserObjectName >	getDependObjects () const

virtual bool	needThreadedCopy () const

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

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

unsigned int	systemNumber () const

virtual bool	enabled () const

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

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

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

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

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

virtual void	jacobianSetup ()

virtual void	residualSetup ()

virtual void	subdomainSetup ()

virtual void	customSetup (const ExecFlagType &)

const ExecFlagEnum &	getExecuteOnEnum () const

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

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

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

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

const UserObject &	getUserObjectBaseByName (const UserObjectName &object_name, bool is_dependency=true) 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 (const std::set< 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

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

std::vector< SubdomainName >	getMaterialPropertyBlockNames (const std::string &name)

std::set< BoundaryID >	getMaterialPropertyBoundaryIDs (const std::string &name)

std::vector< BoundaryName >	getMaterialPropertyBoundaryNames (const std::string &name)

void	checkBlockAndBoundaryCompatibility (std::shared_ptr< MaterialBase > discrete)

std::unordered_map< SubdomainID, std::vector< MaterialBase *> >	buildRequiredMaterials (bool allow_stateful=true)

void	statefulPropertiesAllowed (bool)

bool	getMaterialPropertyCalled () const

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

virtual void	resolveOptionalProperties ()

const GenericMaterialProperty< T, is_ad > &	getPossiblyConstantGenericMaterialPropertyByName (const MaterialPropertyName &prop_name, MaterialData &material_data, const unsigned int state)

virtual const VariableValue &	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)

void	addMooseVariableDependency (MooseVariableFieldBase *var)

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

bool	isImplicit ()

Moose::StateArg	determineState () 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 unsigned int	getElementIDIndex (const std::string &id_parameter_name, unsigned int comp=0) const

virtual unsigned int	getElementIDIndexByName (const std::string &id_name) const

virtual const dof_id_type &	getElementID (const std::string &id_parameter_name, unsigned int comp=0) const

dof_id_type	getElementID (const Elem *elem, unsigned int elem_id_index) const

virtual const dof_id_type &	getElementIDNeighbor (const std::string &id_parameter_name, unsigned int comp=0) const

virtual const dof_id_type &	getElementIDByName (const std::string &id_name) const

virtual const dof_id_type &	getElementIDNeighborByName (const std::string &id_name) const

bool	hasElementID (const std::string &id_name) const

dof_id_type	maxElementID (unsigned int elem_id_index) const

dof_id_type	minElementID (unsigned int elem_id_index) const

bool	areElemIDsIdentical (const std::string &id_name1, const std::string &id_name2) const

std::unordered_map< dof_id_type, std::set< dof_id_type > >	getElemIDMapping (const std::string &id_name1, const std::string &id_name2) const

std::set< dof_id_type >	getAllElemIDs (unsigned int elem_id_index) const

std::set< dof_id_type >	getElemIDsOnBlocks (unsigned int elem_id_index, const std::set< SubdomainID > &blks) const

std::string	PPName () const

bool	containsCompleteHistory () const

bool	isDistributed () const

const std::set< std::string > &	getVectorNames () const

void	buildOutputHideVariableList (std::set< std::string > variable_names)

const std::set< OutputName > &	getOutputs ()

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 GenericOptionalMaterialProperty< T, is_ad > &	getGenericOptionalMaterialProperty (const std::string &name, const unsigned int state=0)

const GenericOptionalMaterialProperty< T, is_ad > &	getGenericOptionalMaterialProperty (const std::string &name, const unsigned int state=0)

const OptionalMaterialProperty< T > &	getOptionalMaterialProperty (const std::string &name, const unsigned int state=0)

const OptionalMaterialProperty< T > &	getOptionalMaterialProperty (const std::string &name, const unsigned int state=0)

const OptionalADMaterialProperty< T > &	getOptionalADMaterialProperty (const std::string &name)

const OptionalADMaterialProperty< T > &	getOptionalADMaterialProperty (const std::string &name)

const OptionalMaterialProperty< T > &	getOptionalMaterialPropertyOld (const std::string &name)

const OptionalMaterialProperty< T > &	getOptionalMaterialPropertyOld (const std::string &name)

const OptionalMaterialProperty< T > &	getOptionalMaterialPropertyOlder (const std::string &name)

const OptionalMaterialProperty< T > &	getOptionalMaterialPropertyOlder (const std::string &name)

MaterialBase &	getMaterial (const std::string &name)

MaterialBase &	getMaterial (const std::string &name)

MaterialBase &	getMaterialByName (const std::string &name, bool no_warn=false)

MaterialBase &	getMaterialByName (const std::string &name, bool no_warn=false)

bool	hasMaterialProperty (const std::string &name)

bool	hasMaterialProperty (const std::string &name)

bool	hasMaterialPropertyByName (const std::string &name)

bool	hasMaterialPropertyByName (const std::string &name)

bool	hasADMaterialProperty (const std::string &name)

bool	hasADMaterialProperty (const std::string &name)

bool	hasADMaterialPropertyByName (const std::string &name)

bool	hasADMaterialPropertyByName (const std::string &name)

bool	hasGenericMaterialProperty (const std::string &name)

bool	hasGenericMaterialProperty (const std::string &name)

bool	hasGenericMaterialPropertyByName (const std::string &name)

bool	hasGenericMaterialPropertyByName (const std::string &name)

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

const Function &	getFunctionByName (const FunctionName &name) const

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

bool	hasFunctionByName (const FunctionName &name) const

bool	isDefaultPostprocessorValue (const std::string &param_name, const unsigned int index=0) const

bool	hasPostprocessor (const std::string &param_name, const unsigned int index=0) const

bool	hasPostprocessorByName (const PostprocessorName &name) const

std::size_t	coupledPostprocessors (const std::string &param_name) const

const PostprocessorName &	getPostprocessorName (const std::string &param_name, const unsigned int index=0) const

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

bool	hasVectorPostprocessorByName (const VectorPostprocessorName &name) const

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

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

Sampler &	getSampler (const std::string &name)

T &	getSamplerByName (const SamplerName &name)

Sampler &	getSamplerByName (const SamplerName &name)

virtual void	meshChanged ()

virtual void	meshDisplaced ()

PerfGraph &	perfGraph ()

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

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

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

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

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

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

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

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

const PostprocessorValue &	getPostprocessorValueOldByName (const PostprocessorName &name) const

const PostprocessorValue &	getPostprocessorValueOldByName (const PostprocessorName &name) const

const PostprocessorValue &	getPostprocessorValueOlderByName (const PostprocessorName &name) const

const PostprocessorValue &	getPostprocessorValueOlderByName (const PostprocessorName &name) const

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

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

bool	isVectorPostprocessorDistributedByName (const VectorPostprocessorName &name) const

bool	isVectorPostprocessorDistributedByName (const VectorPostprocessorName &name) const

const Distribution &	getDistribution (const std::string &name) const

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

const Distribution &	getDistribution (const std::string &name) const

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

const Distribution &	getDistributionByName (const DistributionName &name) const

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

const Distribution &	getDistributionByName (const DistributionName &name) const

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

const Parallel::Communicator &	comm () const

processor_id_type	n_processors () const

processor_id_type	processor_id () const

Static Public Member Functions
static InputParameters	validParams ()

static MooseEnum	qFunctionType ()
	Get the MooseEnum defining the options for how the q function is defined. More...

static MooseEnum	sifModeType ()
	Get the MooseEnum defining the options for the integral to be computed. More...

static void	sort (typename std::vector< T > &vector)

static void	sortDFS (typename std::vector< T > &vector)

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

Public Attributes
const ConsoleStream	_console

Static Public Attributes
static constexpr PropertyValue::id_type	default_property_id

static constexpr PropertyValue::id_type	zero_property_id

static constexpr auto	SYSTEM

static constexpr auto	NAME

Protected Types
enum	QMethod { QMethod::Geometry, QMethod::Topology }
	Enum used to select the method used to compute the q function used in the fracture integrals. More...

enum	PositionType { PositionType::Angle, PositionType::Distance }
	Enum used to define how the distance along the crack front is measured (angle or distance) More...

enum	SifMethod { SifMethod::KI, SifMethod::KII, SifMethod::KIII, SifMethod::T }
	Enum used to select the type of integral to be performed. More...

Protected Member Functions
Real	computeQpIntegral (const std::size_t crack_front_point_index, const Real scalar_q, const RealVectorValue &grad_of_scalar_q)
	Compute the contribution of a specific quadrature point to a fracture integral for a point on the crack front. More...

void	computeAuxFields (RankTwoTensor &aux_stress, RankTwoTensor &grad_disp, RankTwoTensor &aux_strain, RankTwoTensor &aux_disp)
	Compute the auxiliary fields, including the auxiliary stress and the gradient of the auxiliary displacement for the current point (as defined by the combination of _r and _theta) More...

void	computeTFields (RankTwoTensor &aux_stress, RankTwoTensor &grad_disp)
	Compute the auxiliary fields, including the auxiliary stress and the gradient of the auxiliary displacement for the current point (as defined by the combination of _r and _theta) for the T stress. More...

virtual void	addPostprocessorDependencyHelper (const PostprocessorName &name) const override

virtual void	addVectorPostprocessorDependencyHelper (const VectorPostprocessorName &name) const override

virtual void	addUserObjectDependencyHelper (const UserObject &uo) const override

void	addReporterDependencyHelper (const ReporterName &reporter_name) override

const ReporterName &	getReporterName (const std::string &param_name) const

T &	declareRestartableData (const std::string &data_name, Args &&... args)

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

const T &	getRestartableData (const std::string &data_name) const

T &	declareRestartableDataWithContext (const std::string &data_name, void *context, Args &&... args)

T &	declareRecoverableData (const std::string &data_name, Args &&... args)

T &	declareRestartableDataWithObjectName (const std::string &data_name, const std::string &object_name, Args &&... args)

T &	declareRestartableDataWithObjectNameWithContext (const std::string &data_name, const std::string &object_name, void *context, Args &&... args)

std::string	restartableName (const std::string &data_name) const

const T &	getMeshProperty (const std::string &data_name, const std::string &prefix)

const T &	getMeshProperty (const std::string &data_name)

bool	hasMeshProperty (const std::string &data_name, const std::string &prefix) const

bool	hasMeshProperty (const std::string &data_name, const std::string &prefix) const

bool	hasMeshProperty (const std::string &data_name) const

bool	hasMeshProperty (const std::string &data_name) const

std::string	meshPropertyName (const std::string &data_name) const

PerfID	registerTimedSection (const std::string &section_name, const unsigned int level) const

PerfID	registerTimedSection (const std::string &section_name, const unsigned int level, const std::string &live_message, const bool print_dots=true) const

std::string	timedSectionName (const std::string &section_name) const

virtual bool	hasBlockMaterialPropertyHelper (const std::string &prop_name)

void	initializeBlockRestrictable (const MooseObject *moose_object)

Moose::CoordinateSystemType	getBlockCoordSystem ()

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

const GenericVectorVariableValue< false > &	coupledGenericVectorValue (const std::string &var_name, unsigned int comp) const

const GenericVectorVariableValue< true > &	coupledGenericVectorValue (const std::string &var_name, unsigned int comp) const

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

void	checkWritableVar (MooseWritableVariable *var)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

const ADVariableValue &	adZeroValue () const

const ADVariableGradient &	adZeroGradient () const

const ADVariableSecond &	adZeroSecond () const

const GenericVariableValue< is_ad > &	genericZeroValue ()

const GenericVariableValue< false > &	genericZeroValue ()

const GenericVariableValue< true > &	genericZeroValue ()

const GenericVariableGradient< is_ad > &	genericZeroGradient ()

const GenericVariableGradient< false > &	genericZeroGradient ()

const GenericVariableGradient< true > &	genericZeroGradient ()

const GenericVariableSecond< is_ad > &	genericZeroSecond ()

const GenericVariableSecond< false > &	genericZeroSecond ()

const GenericVariableSecond< true > &	genericZeroSecond ()

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

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

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

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

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

VectorPostprocessorValue &	declareVector (const std::string &vector_name)

const T &	getReporterValue (const std::string &param_name, const std::size_t time_index=0)

const T &	getReporterValue (const std::string &param_name, ReporterMode mode, const std::size_t time_index=0)

const T &	getReporterValue (const std::string &param_name, const std::size_t time_index=0)

const T &	getReporterValue (const std::string &param_name, ReporterMode mode, const std::size_t time_index=0)

const T &	getReporterValueByName (const ReporterName &reporter_name, const std::size_t time_index=0)

const T &	getReporterValueByName (const ReporterName &reporter_name, ReporterMode mode, const std::size_t time_index=0)

const T &	getReporterValueByName (const ReporterName &reporter_name, const std::size_t time_index=0)

const T &	getReporterValueByName (const ReporterName &reporter_name, ReporterMode mode, const std::size_t time_index=0)

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

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

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

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

bool	hasReporterValueByName (const ReporterName &reporter_name) const

bool	hasReporterValueByName (const ReporterName &reporter_name) const

bool	hasReporterValueByName (const ReporterName &reporter_name) const

bool	hasReporterValueByName (const ReporterName &reporter_name) const

const GenericMaterialProperty< T, is_ad > *	defaultGenericMaterialProperty (const std::string &name)

const GenericMaterialProperty< T, is_ad > *	defaultGenericMaterialProperty (const std::string &name)

const MaterialProperty< T > *	defaultMaterialProperty (const std::string &name)

const MaterialProperty< T > *	defaultMaterialProperty (const std::string &name)

const ADMaterialProperty< T > *	defaultADMaterialProperty (const std::string &name)

const ADMaterialProperty< T > *	defaultADMaterialProperty (const std::string &name)

Static Protected Member Functions
static std::string	meshPropertyName (const std::string &data_name, const std::string &prefix)

Protected Attributes
std::size_t	_ndisp
	Number of displacement components. More...

const CrackFrontDefinition *const	_crack_front_definition
	Pointer to the crack front definition object. More...

bool	_treat_as_2d
	Whether to treat a model as 2D for computation of fracture integrals. More...

bool	_using_mesh_cutter
	Is the crack defined by an XFEM cutter mesh. More...

const GenericMaterialProperty< RankTwoTensor, is_ad > &	_stress
	Reference to the stress tensor computed by the material models. More...

const GenericMaterialProperty< RankTwoTensor, is_ad > &	_strain
	Reference to the strain tensor computed by the material models. More...

std::vector< MooseVariableFEBase * >	_fe_vars
	Vector of all coupled variables. More...

const libMesh::FEType &	_fe_type
	FEType object defining order and family of displacement variables. More...

std::vector< const VariableValue * >	_disp
	Displacement variables. More...

std::vector< const VariableGradient * >	_grad_disp
	Gradient of displacements. More...

const bool	_has_temp
	Whether the temperature variable is coupled. More...

const VariableGradient &	_grad_temp
	Gradient of temperature. More...

const MaterialProperty< Real > *	_functionally_graded_youngs_modulus_crack_dir_gradient
	Spatial derivative of the youngs modulus in the crack direction. More...

const MaterialProperty< Real > *	_functionally_graded_youngs_modulus
	Spatial elasticity modulus variable for FGM. More...

Real	_K_factor
	Conversion factor applied to convert interaction integral to stress intensity factor K. More...

bool	_has_symmetry_plane
	Whether the crack plane is also a symmetry plane in the model. More...

Real	_poissons_ratio
	Poisson's ratio of the material. More...

Real	_youngs_modulus
	Young's modulus of the material. More...

const bool	_fgm_crack
	Whether to consider interaction integral and material properties for a crack in functionally graded material. More...

std::size_t	_ring_index
	Index of the ring for the integral computed by this object. More...

const GenericMaterialProperty< RankTwoTensor, is_ad > *const	_total_deigenstrain_dT
	Derivative of the total eigenstrain with respect to temperature. More...

bool	_has_additional_eigenstrain
	Whether the user chooses to add other eigenstrain influence (e.g. irradiation-induced) More...

const VariableGradient *	_additional_eigenstrain_gradient_00
	Gradient used to add contribution to the interaction integral (XX) More...

const VariableGradient *	_additional_eigenstrain_gradient_01
	Gradient used to add contribution to the interaction integral (XY) More...

const VariableGradient *	_additional_eigenstrain_gradient_11
	Gradient used to add contribution to the interaction integral (YY) More...

const VariableGradient *	_additional_eigenstrain_gradient_22
	Gradient used to add contribution to the interaction integral (ZZ) More...

const VariableGradient *	_additional_eigenstrain_gradient_02
	Gradient used to add contribution to the interaction integral (XZ) More...

const VariableGradient *	_additional_eigenstrain_gradient_12
	Gradient used to add contribution to the interaction integral (YZ) More...

std::vector< Real >	_q_curr_elem
	Vector of q function values for the nodes in the current element. More...

const std::vector< std::vector< Real > > *	_phi_curr_elem
	Vector of shape function values for the current element. More...

const std::vector< std::vector< libMesh::RealGradient > > *	_dphi_curr_elem
	Vector of gradients of shape function values for the current element. More...

Real	_kappa
	Kappa Lame constant. More...

Real	_shear_modulus
	Shear modulus. More...

Real	_r
	Radial distance from the current point to the crack front. More...

Real	_theta
	Angle of current point relative to the crack front. More...

unsigned int	_qp
	Current quadrature point index. More...

enum InteractionIntegralTempl::QMethod	_q_function_type

enum InteractionIntegralTempl::PositionType	_position_type

enum InteractionIntegralTempl::SifMethod	_sif_mode

MooseMesh &	_mesh

const Elem *const &	_current_elem

const Real &	_current_elem_volume

const MooseArray< Point > &	_q_point

const QBase *const &	_qrule

const MooseArray< Real > &	_JxW

const MooseArray< Real > &	_coord

SubProblem &	_subproblem

FEProblemBase &	_fe_problem

SystemBase &	_sys

const THREAD_ID	_tid

Assembly &	_assembly

const Moose::CoordinateSystemType &	_coord_sys

const bool	_duplicate_initial_execution

std::set< std::string >	_depend_uo

const bool &	_enabled

MooseApp &	_app

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

MooseApp &	_restartable_app

const std::string	_restartable_system_name

const THREAD_ID	_restartable_tid

const bool	_restartable_read_only

FEProblemBase &	_mci_feproblem

FEProblemBase &	_mdi_feproblem

MooseApp &	_pg_moose_app

const std::string	_prefix

const MaterialData *	_blk_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

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

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

const std::string	_vpp_name

FEProblemBase &	_vpp_fe_problem

const MooseEnum &	_parallel_type

const Parallel::Communicator &	_communicator


VectorPostprocessorValue &	_x
	Vectors computed by this VectorPostprocessor: x,y,z coordinates, position of nodes along crack front, and interaction integral. More...

VectorPostprocessorValue &	_y

VectorPostprocessorValue &	_z

VectorPostprocessorValue &	_position

VectorPostprocessorValue &	_interaction_integral


const GenericMaterialProperty< RankThreeTensor, is_ad > *	_eigenstrain_gradient
	Pointers to optionally-used eigenstrain gradient and body force. More...

const MaterialProperty< RealVectorValue > *	_body_force

Static Protected Attributes
static const std::string	_interpolated_old

static const std::string	_interpolated_older

Detailed Description

template<bool is_ad>
class InteractionIntegralTempl< is_ad >

This vectorpostprocessor computes the Interaction Integral, which is used to compute various fracture mechanics parameters at a crack tip, including KI, KII, KIII, and the T stress.

Definition at line 23 of file InteractionIntegral.h.

Member Enumeration Documentation

◆ PositionType

template<bool is_ad>

enum InteractionIntegralTempl::PositionType

strongprotected

Enum used to define how the distance along the crack front is measured (angle or distance)

Enumerator
Angle
Distance

Definition at line 165 of file InteractionIntegral.h.

165 { Angle, Distance } _position_type;

InteractionIntegralTempl::_position_type

enum InteractionIntegralTempl::PositionType _position_type

◆ QMethod

template<bool is_ad>

enum InteractionIntegralTempl::QMethod

strongprotected

Enum used to select the method used to compute the q function used in the fracture integrals.

Enumerator
Geometry
Topology

Definition at line 161 of file InteractionIntegral.h.

161 { Geometry, Topology } _q_function_type;

InteractionIntegralTempl::_q_function_type

enum InteractionIntegralTempl::QMethod _q_function_type

◆ SifMethod

template<bool is_ad>

enum InteractionIntegralTempl::SifMethod

strongprotected

Enum used to select the type of integral to be performed.

Enumerator
KI
KII
KIII
T

Definition at line 168 of file InteractionIntegral.h.

168 { KI, KII, KIII, T } _sif_mode;

InteractionIntegralTempl::_sif_mode

enum InteractionIntegralTempl::SifMethod _sif_mode

Constructor & Destructor Documentation

◆ InteractionIntegralTempl()

template<bool is_ad>

InteractionIntegralTempl< is_ad >::InteractionIntegralTempl ( const InputParameters & parameters )

Definition at line 116 of file InteractionIntegral.C.

   : ElementVectorPostprocessor(parameters),
     _ndisp(coupledComponents("displacements")),
     _crack_front_definition(&getUserObject<CrackFrontDefinition>("crack_front_definition")),
     _treat_as_2d(false),
     _stress(getGenericMaterialPropertyByName<RankTwoTensor, is_ad>("stress")),
     _strain(getGenericMaterialPropertyByName<RankTwoTensor, is_ad>("elastic_strain")),
     _fe_vars(getCoupledMooseVars()),
     _fe_type(_fe_vars[0]->feType()),
     _disp(coupledValues("displacements")),
     _grad_disp(3),
     _has_temp(isCoupled("temperature")),
     _grad_temp(_has_temp ? coupledGradient("temperature") : _grad_zero),
     _functionally_graded_youngs_modulus_crack_dir_gradient(
         isParamSetByUser("functionally_graded_youngs_modulus_crack_dir_gradient")
             ? &getMaterialProperty<Real>("functionally_graded_youngs_modulus_crack_dir_gradient")
             : nullptr),
     _functionally_graded_youngs_modulus(
         isParamSetByUser("functionally_graded_youngs_modulus")
             ? &getMaterialProperty<Real>("functionally_graded_youngs_modulus")
             : nullptr),
     _K_factor(getParam<Real>("K_factor")),
     _has_symmetry_plane(isParamValid("symmetry_plane")),
     _poissons_ratio(getParam<Real>("poissons_ratio")),
     _youngs_modulus(getParam<Real>("youngs_modulus")),
     _fgm_crack(isParamSetByUser("functionally_graded_youngs_modulus_crack_dir_gradient") &&
                isParamSetByUser("functionally_graded_youngs_modulus")),
     _ring_index(getParam<unsigned int>("ring_index")),
     _total_deigenstrain_dT(
         hasMaterialProperty<RankTwoTensor>("total_deigenstrain_dT")
             ? &getGenericMaterialProperty<RankTwoTensor, is_ad>("total_deigenstrain_dT")
             : nullptr),
     _has_additional_eigenstrain(false),
     _additional_eigenstrain_gradient_00(isCoupled("additional_eigenstrain_00")
                                             ? &coupledGradient("additional_eigenstrain_00")
                                             : nullptr),
     _additional_eigenstrain_gradient_01(isCoupled("additional_eigenstrain_01")
                                             ? &coupledGradient("additional_eigenstrain_01")
                                             : nullptr),
     _additional_eigenstrain_gradient_11(isCoupled("additional_eigenstrain_11")
                                             ? &coupledGradient("additional_eigenstrain_11")
                                             : nullptr),
     _additional_eigenstrain_gradient_22(isCoupled("additional_eigenstrain_22")
                                             ? &coupledGradient("additional_eigenstrain_22")
                                             : nullptr),
     _additional_eigenstrain_gradient_02(isCoupled("additional_eigenstrain_02")
                                             ? &coupledGradient("additional_eigenstrain_02")
                                             : nullptr),
     _additional_eigenstrain_gradient_12(isCoupled("additional_eigenstrain_12")
                                             ? &coupledGradient("additional_eigenstrain_12")
                                             : nullptr),
     _q_function_type(getParam<MooseEnum>("q_function_type").template getEnum<QMethod>()),
     _position_type(getParam<MooseEnum>("position_type").template getEnum<PositionType>()),
     _sif_mode(getParam<MooseEnum>("sif_mode").template getEnum<SifMethod>()),
     _x(declareVector("x")),
     _y(declareVector("y")),
     _z(declareVector("z")),
     _position(declareVector("id")),
     _interaction_integral(declareVector("II_" + Moose::stringify(getParam<MooseEnum>("sif_mode")) +
                                         "_" + Moose::stringify(_ring_index))),
     _eigenstrain_gradient(nullptr),
     _body_force(nullptr)
 {
   if (_has_temp && !_total_deigenstrain_dT)
     mooseError("InteractionIntegral Error: To include thermal strain term in interaction integral, "
                "must both couple temperature in DomainIntegral block and compute "
                "total_deigenstrain_dT using ThermalFractureIntegral material model.");
 
   if ((!_functionally_graded_youngs_modulus &&
        _functionally_graded_youngs_modulus_crack_dir_gradient) ||
       (_functionally_graded_youngs_modulus &&
        !_functionally_graded_youngs_modulus_crack_dir_gradient))
     paramError("functionally_graded_youngs_modulus_crack_dir_gradient",
                "You have selected to compute the interaction integral for a crack in FGM. That "
                "selection requires the user to provide a spatially varying elasticity modulus "
                "that "
                "defines the transition of material properties (i.e. "
                "'functionally_graded_youngs_modulus') and its "
                "spatial derivative in the crack direction (i.e. "
                "'functionally_graded_youngs_modulus_crack_dir_gradient').");
 
   // plane strain
   _kappa = 3.0 - 4.0 * _poissons_ratio;
   _shear_modulus = _youngs_modulus / (2.0 * (1.0 + _poissons_ratio));
 
   // Checking for consistency between mesh size and length of the provided displacements vector
   if (_ndisp != _mesh.dimension())
     paramError("displacements",
                "InteractionIntegral Error: number of variables supplied in 'displacements' must "
                "match the mesh dimension.");
 
   // fetch gradients of coupled variables
   for (std::size_t i = 0; i < _ndisp; ++i)
     _grad_disp[i] = &coupledGradient("displacements", i);
 
   // set unused dimensions to zero
   for (std::size_t i = _ndisp; i < 3; ++i)
     _grad_disp[i] = &_grad_zero;
 
   if (isParamValid("eigenstrain_gradient"))
   {
     _eigenstrain_gradient =
         &getGenericMaterialProperty<RankThreeTensor, is_ad>("eigenstrain_gradient");
     if (_total_deigenstrain_dT)
       paramError("eigenstrain_gradient",
                  "eigenstrain_gradient cannot be specified for materials that provide the "
                  "total_deigenstrain_dT material property");
   }
   if (isParamValid("body_force"))
     _body_force = &getMaterialProperty<RealVectorValue>("body_force");
 
   _has_additional_eigenstrain = _additional_eigenstrain_gradient_00 &&
                                 _additional_eigenstrain_gradient_01 &&
                                 _additional_eigenstrain_gradient_11;
 
   if (_mesh.dimension() == 3 && _has_additional_eigenstrain)
     if (!(_additional_eigenstrain_gradient_22 && _additional_eigenstrain_gradient_12 &&
           _additional_eigenstrain_gradient_02))
       paramError("additional_eigenstrain_gradient_12",
                  "If additional eigenstrains are provided for the computation of the interaction "
                  "integral in three dimensions, make sure the six components are provided.");
 
   if (_has_additional_eigenstrain)
     mooseInfo("A generic eigenstrain provided by the user will be considered in the interaction "
               "integral (via domain integral action).");
 }

Member Function Documentation

◆ computeAuxFields()

template<bool is_ad>

void InteractionIntegralTempl< is_ad >::computeAuxFields	(	RankTwoTensor &	aux_stress,
		RankTwoTensor &	grad_disp,
		RankTwoTensor &	aux_strain,
		RankTwoTensor &	aux_disp
	)

protected

Compute the auxiliary fields, including the auxiliary stress and the gradient of the auxiliary displacement for the current point (as defined by the combination of _r and _theta)

Parameters

aux_stress	Auxiliary stress – computed in this method
grad_disp	Gradient of auxiliary displacement – computed in this method
grad_disp	Auxiliary strain (FGM) – computed in this method
aux_disp	Auxiliary displacement – computed in this method only for axisymmetric (RZ) meshes

Definition at line 611 of file InteractionIntegral.C.

 {
   RealVectorValue k(0.0);
   if (_sif_mode == SifMethod::KI)
     k(0) = 1.0;
   else if (_sif_mode == SifMethod::KII)
     k(1) = 1.0;
   else if (_sif_mode == SifMethod::KIII)
     k(2) = 1.0;
 
   Real t = _theta;
   Real t2 = _theta / 2.0;
   Real tt2 = 3.0 * _theta / 2.0;
   Real st = std::sin(t);
   Real ct = std::cos(t);
   Real st2 = std::sin(t2);
   Real ct2 = std::cos(t2);
   Real stt2 = std::sin(tt2);
   Real ctt2 = std::cos(tt2);
   Real ct2sq = Utility::pow<2>(ct2);
   Real st2sq = Utility::pow<2>(st2);
   Real ct2cu = Utility::pow<3>(ct2);
   Real sqrt2PiR = std::sqrt(2.0 * libMesh::pi * _r);
 
   // Calculate auxiliary stress tensor
   aux_stress.zero();
 
   aux_stress(0, 0) =
       1.0 / sqrt2PiR * (k(0) * ct2 * (1.0 - st2 * stt2) - k(1) * st2 * (2.0 + ct2 * ctt2));
   aux_stress(1, 1) = 1.0 / sqrt2PiR * (k(0) * ct2 * (1.0 + st2 * stt2) + k(1) * st2 * ct2 * ctt2);
   aux_stress(0, 1) = 1.0 / sqrt2PiR * (k(0) * ct2 * st2 * ctt2 + k(1) * ct2 * (1.0 - st2 * stt2));
   aux_stress(0, 2) = -1.0 / sqrt2PiR * k(2) * st2;
   aux_stress(1, 2) = 1.0 / sqrt2PiR * k(2) * ct2;
   // plane stress
   // Real s33 = 0;
   // plane strain
   aux_stress(2, 2) = _poissons_ratio * (aux_stress(0, 0) + aux_stress(1, 1));
 
   aux_stress(1, 0) = aux_stress(0, 1);
   aux_stress(2, 0) = aux_stress(0, 2);
   aux_stress(2, 1) = aux_stress(1, 2);
 
   // Calculate auxiliary strain tensor (only needed to FGM)
   if (_fgm_crack)
   {
     aux_strain.zero();
     RankFourTensor spatial_elasticity_tensor;
     spatial_elasticity_tensor.fillSymmetricIsotropicEandNu(
         (*_functionally_graded_youngs_modulus)[_qp], _poissons_ratio);
     aux_strain = spatial_elasticity_tensor.invSymm() * aux_stress;
   }
 
   // Calculate x1 derivative of auxiliary displacements
   grad_disp.zero();
   grad_disp(0, 0) = k(0) / (4.0 * _shear_modulus * sqrt2PiR) *
                         (ct * ct2 * _kappa + ct * ct2 - 2.0 * ct * ct2cu + st * st2 * _kappa +
                          st * st2 - 6.0 * st * st2 * ct2sq) +
                     k(1) / (4.0 * _shear_modulus * sqrt2PiR) *
                         (ct * st2 * _kappa + ct * st2 + 2.0 * ct * st2 * ct2sq - st * ct2 * _kappa +
                          3.0 * st * ct2 - 6.0 * st * ct2cu);
 
   grad_disp(0, 1) = k(0) / (4.0 * _shear_modulus * sqrt2PiR) *
                         (ct * st2 * _kappa + ct * st2 - 2.0 * ct * st2 * ct2sq - st * ct2 * _kappa -
                          5.0 * st * ct2 + 6.0 * st * ct2cu) +
                     k(1) / (4.0 * _shear_modulus * sqrt2PiR) *
                         (-ct * ct2 * _kappa + 3.0 * ct * ct2 - 2.0 * ct * ct2cu -
                          st * st2 * _kappa + 3.0 * st * st2 - 6.0 * st * st2 * ct2sq);
 
   grad_disp(0, 2) = k(2) / (_shear_modulus * sqrt2PiR) * (st2 * ct - ct2 * st);
 
   // Calculate aux displacement field (used in RZ calculations) - Warp3D manual
   aux_disp.zero();
   if (getBlockCoordSystem() == Moose::COORD_RZ)
   {
     const Real prefactor_rz = std::sqrt(_r / 2.0 / libMesh::pi) / 2.0 / _shear_modulus;
     aux_disp(0, 0) = prefactor_rz * ((k(0) * ct2 * (_kappa - 1.0 + 2.0 * st2sq)) +
                                      (k(1) * st2 * (_kappa + 1.0 + 2.0 * ct2sq)));
     aux_disp(0, 1) = prefactor_rz * ((k(0) * st2 * (_kappa + 1.0 - 2.0 * ct2sq)) -
                                      (k(1) * ct2 * (_kappa - 1.0 - 2.0 * st2sq)));
     aux_disp(0, 2) = prefactor_rz * k(2) * 4.0 * st2;
   }
 }

◆ computeQpIntegral()

template<bool is_ad>

Real InteractionIntegralTempl< is_ad >::computeQpIntegral	(	const std::size_t	crack_front_point_index,
		const Real	scalar_q,
		const RealVectorValue &	grad_of_scalar_q
	)

protected

Compute the contribution of a specific quadrature point to a fracture integral for a point on the crack front.

Parameters

crack_front_point_index	Index to the crack front point
scalar_q	Magnitude of the q vector at the quadrature point
grad_of_scalar_q	Gradient of the magnitude of the q vector at the quadrature point

Returns: Contribution of this quadrature point to the integral

Definition at line 276 of file InteractionIntegral.C.

 {
   // If q is zero, then dq is also zero, so all terms in the interaction integral would
   // return zero. As such, let us avoid unnecessary, frequent computations
   if (scalar_q < TOLERANCE * TOLERANCE * TOLERANCE)
     return 0.0;
 
   // In the crack front coordinate system, the crack direction is (1,0,0)
   RealVectorValue crack_direction(1.0, 0.0, 0.0);
 
   // Calculate (r,theta) position of qp relative to crack front
   _crack_front_definition->calculateRThetaToCrackFront(
       _q_point[_qp], crack_front_point_index, _r, _theta);
 
   RankTwoTensor aux_stress;
   RankTwoTensor aux_du;
   RankTwoTensor aux_strain;
   RankTwoTensor aux_disp;
 
   if (_sif_mode == SifMethod::KI || _sif_mode == SifMethod::KII || _sif_mode == SifMethod::KIII)
     computeAuxFields(aux_stress, aux_du, aux_strain, aux_disp);
   else if (_sif_mode == SifMethod::T)
     computeTFields(aux_stress, aux_du);
 
   auto grad_disp = RankTwoTensor::initializeFromRows(
       (*_grad_disp[0])[_qp], (*_grad_disp[1])[_qp], (*_grad_disp[2])[_qp]);
 
   // Rotate stress, strain, displacement and temperature to crack front coordinate system
   RealVectorValue grad_q_cf =
       _crack_front_definition->rotateToCrackFrontCoords(grad_of_scalar_q, crack_front_point_index);
   RankTwoTensor grad_disp_cf =
       _crack_front_definition->rotateToCrackFrontCoords(grad_disp, crack_front_point_index);
   RankTwoTensor stress_cf = _crack_front_definition->rotateToCrackFrontCoords(
       MetaPhysicL::raw_value((_stress)[_qp]), crack_front_point_index);
   RankTwoTensor strain_cf = _crack_front_definition->rotateToCrackFrontCoords(
       MetaPhysicL::raw_value((_strain)[_qp]), crack_front_point_index);
   RealVectorValue grad_temp_cf =
       _crack_front_definition->rotateToCrackFrontCoords(_grad_temp[_qp], crack_front_point_index);
 
   if (getBlockCoordSystem() == Moose::COORD_RZ)
     strain_cf(2, 2) = 0.0;
 
   RankTwoTensor dq;
   dq(0, 0) = crack_direction(0) * grad_q_cf(0);
   dq(0, 1) = crack_direction(0) * grad_q_cf(1);
   dq(0, 2) = crack_direction(0) * grad_q_cf(2);
 
   // Calculate interaction integral terms
 
   // Term1 = stress * x1-derivative of aux disp * dq
   RankTwoTensor tmp1 = dq * stress_cf;
   Real term1 = aux_du.doubleContraction(tmp1);
 
   // Term2 = aux stress * x1-derivative of disp * dq
   RankTwoTensor tmp2 = dq * aux_stress;
   Real term2 = grad_disp_cf(0, 0) * tmp2(0, 0) + grad_disp_cf(1, 0) * tmp2(0, 1) +
                grad_disp_cf(2, 0) * tmp2(0, 2);
 
   // Term3 = aux stress * strain * dq_x   (= stress * aux strain * dq_x)
   Real term3 = -dq(0, 0) * aux_stress.doubleContraction(strain_cf);
 
   // Term4 (thermal strain term) = q * aux_stress * alpha * dtheta_x
   // - the term including the derivative of alpha is not implemented
   Real term4 = 0.0;
   if (_has_temp)
   {
     Real sigma_alpha =
         aux_stress.doubleContraction(MetaPhysicL::raw_value((*_total_deigenstrain_dT)[_qp]));
     term4 = scalar_q * sigma_alpha * grad_temp_cf(0);
   }
 
   Real term4a = 0.0; // Gradient of arbitrary eigenstrain material contribution
   if (_eigenstrain_gradient)
   {
     // Thermal strain gradient term in Nakamura and Parks, IJSS, 1992:
     // alpha * dT/dx_k*aux_stress*scalar_q
     // Generalization to the gradient of an arbitrary eigenstrain:
     // d_eigenstrain/dx_k*aux_stress*scalar_q
     const RealVectorValue & crack_dir =
         _crack_front_definition->getCrackDirection(crack_front_point_index);
     RankTwoTensor eigenstrain_grad_in_crack_dir =
         crack_dir * MetaPhysicL::raw_value((*_eigenstrain_gradient)[_qp]);
     term4a = scalar_q * aux_stress.doubleContraction(eigenstrain_grad_in_crack_dir);
   }
 
   // Gradient of arbitrary eigenstrain variable contribution.
   // Note that, unlike approaches 4 and 4a, approach 4b can integrate any type of eigenstrain, i.e.
   // it does not have to be thermal or volumetric.
   Real term4b = 0.0;
   // Term4b (eigen strain term) = q * aux_stress * d{eigen_strain}_x
   if (_has_additional_eigenstrain)
   {
     RankTwoTensor I_x(RankTwoTensor::initNone);
     I_x(0, 0) = 1.0;
     RankTwoTensor I_xy(RankTwoTensor::initNone);
     I_xy(0, 1) = 1.0;
     I_xy(1, 0) = 1.0;
     RankTwoTensor I_y(RankTwoTensor::initNone);
     I_y(1, 1) = 1.0;
     RankTwoTensor I_z(RankTwoTensor::initNone);
     I_z(2, 2) = 1.0;
 
     const RealVectorValue & crack_dir =
         _crack_front_definition->getCrackDirection(crack_front_point_index);
 
     const auto eigenstrain_gradient_x =
         I_x.mixedProductJkI((*_additional_eigenstrain_gradient_00)[_qp]);
     const auto eigenstrain_gradient_xy =
         I_xy.mixedProductJkI((*_additional_eigenstrain_gradient_01)[_qp]);
     const auto eigenstrain_gradient_y =
         I_y.mixedProductJkI((*_additional_eigenstrain_gradient_11)[_qp]);
 
     GenericRankThreeTensor<is_ad> eigenstrain_gradient_z;
     if (_additional_eigenstrain_gradient_22)
       eigenstrain_gradient_z = I_z.mixedProductJkI((*_additional_eigenstrain_gradient_22)[_qp]);
 
     RankTwoTensor eigenstrain_grad_in_crack_dir =
         crack_dir * (MetaPhysicL::raw_value(eigenstrain_gradient_x) +
                      MetaPhysicL::raw_value(eigenstrain_gradient_xy) +
                      MetaPhysicL::raw_value(eigenstrain_gradient_y) +
                      MetaPhysicL::raw_value(eigenstrain_gradient_z));
 
     GenericRankThreeTensor<is_ad> eigenstrain_gradient_xz;
     GenericRankThreeTensor<is_ad> eigenstrain_gradient_yz;
     // Add terms for three dimensions
     if (_mesh.dimension() == 3)
     {
       // Add influence of xz and yz.
       RankTwoTensor I_xz(RankTwoTensor::initNone);
       I_xz(0, 2) = I_xz(2, 0) = 1.0;
       eigenstrain_gradient_xz = I_xz.mixedProductJkI((*_additional_eigenstrain_gradient_02)[_qp]);
       RankTwoTensor I_yz(RankTwoTensor::initNone);
       I_yz(1, 2) = I_yz(2, 1) = 1.0;
       eigenstrain_gradient_yz = I_yz.mixedProductJkI((*_additional_eigenstrain_gradient_12)[_qp]);
 
       eigenstrain_grad_in_crack_dir = crack_dir * (MetaPhysicL::raw_value(eigenstrain_gradient_xz) +
                                                    MetaPhysicL::raw_value(eigenstrain_gradient_yz) +
                                                    MetaPhysicL::raw_value(eigenstrain_gradient_x) +
                                                    MetaPhysicL::raw_value(eigenstrain_gradient_xy) +
                                                    MetaPhysicL::raw_value(eigenstrain_gradient_y) +
                                                    MetaPhysicL::raw_value(eigenstrain_gradient_z));
     }
 
     term4b = scalar_q * aux_stress.doubleContraction(eigenstrain_grad_in_crack_dir);
   }
 
   Real term5 = 0.0; // Body force
   if (_body_force)
   {
     // Body force term in Nakamura and Parks, IJSS, 1992:
     // b_i*aux_du*crack_dir*scalar_q
     const RealVectorValue & crack_dir =
         _crack_front_definition->getCrackDirection(crack_front_point_index);
     term5 = -scalar_q * MetaPhysicL::raw_value((*_body_force)[_qp]) * aux_du * crack_dir;
   }
 
   Real term6 = 0.0;
   if (_fgm_crack && scalar_q != 0)
   {
     // See Equation 49 of J.-H. Kim and G. H. Paulino. Consistent formulations of the interaction
     // integral method for fracture of functionally graded materials. Journal of Applied Mechanics,
     // 72(3) 351-364 2004.
 
     // Term 6_1 = Cijkl,j(x) epsilon_{kl}^{aux} disp_{i,1}
     RankTwoTensor cijklj_epsilonkl_aux; // Temporary second order tensor.
     RankFourTensor cijklj;
     cijklj.fillSymmetricIsotropicEandNu(1.0, _poissons_ratio);
     cijklj *= (*_functionally_graded_youngs_modulus_crack_dir_gradient)[_qp];
     cijklj_epsilonkl_aux = cijklj * aux_strain;
 
     const Real term6_a = grad_disp_cf(0, 0) * cijklj_epsilonkl_aux(0, 0) +
                          grad_disp_cf(1, 0) * cijklj_epsilonkl_aux(0, 1) +
                          grad_disp_cf(2, 0) * cijklj_epsilonkl_aux(0, 2);
 
     // Term 6_2 = -Cijkl,1(x) epsilon_{kl} epsilon_{ij}^{aux}
     RankTwoTensor cijkl1_epsilonkl;
     cijkl1_epsilonkl = -cijklj * strain_cf;
     const Real term6_b = cijkl1_epsilonkl.doubleContraction(aux_strain);
 
     term6 = (term6_a + term6_b) * scalar_q;
   }
 
   // Add terms for possibly RZ (axisymmetric problem).
   // See Nahta and Moran, 1993, Domain integrals for axisymmetric interface crack problems, 30(15)
   // Note: Sign problem in Eq. (26). Also Kolosov constant is wrong in Eq. (13)
   Real term7 = 0.0;
   if (getBlockCoordSystem() == Moose::COORD_RZ)
   {
     const Real u_r = (*_disp[0])[_qp];
     const Real radius_qp = _q_point[_qp](0);
 
     const Real term7a =
         (u_r / radius_qp * aux_stress(2, 2) + aux_disp(0, 0) / radius_qp * stress_cf(2, 2) -
          aux_stress.doubleContraction(strain_cf)) /
         radius_qp;
 
     // term7b1: sigma_aux_ralpha * u_alpha,r
     const Real term7b1 = aux_stress(0, 0) * grad_disp_cf(0, 0) +
                          aux_stress(0, 1) * grad_disp_cf(1, 0) +
                          aux_stress(0, 2) * grad_disp_cf(2, 0);
 
     const Real term7b = 1.0 / radius_qp *
                         (term7b1 - aux_stress(2, 2) * grad_disp_cf(0, 0) +
                          stress_cf(2, 2) * (aux_du(0, 0) - aux_disp(0, 0) / radius_qp));
 
     term7 = (term7a + term7b) * scalar_q;
   }
 
   Real q_avg_seg = 1.0;
   if (!_crack_front_definition->treatAs2D())
   {
     q_avg_seg =
         (_crack_front_definition->getCrackFrontForwardSegmentLength(crack_front_point_index) +
          _crack_front_definition->getCrackFrontBackwardSegmentLength(crack_front_point_index)) /
         2.0;
   }
 
   Real eq = term1 + term2 + term3 + term4 + term4a + term4b + term5 + term6 + term7;
 
   if (getBlockCoordSystem() == Moose::COORD_RZ)
   {
     const Real radius_qp = _q_point[_qp](0);
     eq *= radius_qp;
 
     std::size_t num_crack_front_points = _crack_front_definition->getNumCrackFrontPoints();
     if (num_crack_front_points != 1)
       mooseError("Crack front has more than one point, but this is a 2D-RZ problem. Please revise "
                  "your input file.");
 
     const Point * crack_front_point = _crack_front_definition->getCrackFrontPoint(0);
     eq = eq / (*crack_front_point)(0);
   }
 
   return eq / q_avg_seg;
 }

◆ computeTFields()

template<bool is_ad>

void InteractionIntegralTempl< is_ad >::computeTFields	(	RankTwoTensor &	aux_stress,
		RankTwoTensor &	grad_disp
	)

protected

Compute the auxiliary fields, including the auxiliary stress and the gradient of the auxiliary displacement for the current point (as defined by the combination of _r and _theta) for the T stress.

Parameters

aux_stress	Auxiliary stress – computed in this method
grad_disp	Gradient of auxiliary displacement – computed in this method

Definition at line 699 of file InteractionIntegral.C.

 {
 
   Real t = _theta;
   Real st = std::sin(t);
   Real ct = std::cos(t);
   Real stsq = Utility::pow<2>(st);
   Real ctsq = Utility::pow<2>(ct);
   Real ctcu = Utility::pow<3>(ct);
   Real oneOverPiR = 1.0 / (libMesh::pi * _r);
 
   aux_stress.zero();
   aux_stress(0, 0) = -oneOverPiR * ctcu;
   aux_stress(0, 1) = -oneOverPiR * st * ctsq;
   aux_stress(1, 0) = -oneOverPiR * st * ctsq;
   aux_stress(1, 1) = -oneOverPiR * ct * stsq;
   aux_stress(2, 2) = -oneOverPiR * _poissons_ratio * (ctcu + ct * stsq);
 
   grad_disp.zero();
   grad_disp(0, 0) = oneOverPiR / (4.0 * _youngs_modulus) *
                     (ct * (4.0 * Utility::pow<2>(_poissons_ratio) - 3.0 + _poissons_ratio) -
                      std::cos(3.0 * t) * (1.0 + _poissons_ratio));
   grad_disp(0, 1) = -oneOverPiR / (4.0 * _youngs_modulus) *
                     (st * (4.0 * Utility::pow<2>(_poissons_ratio) - 3.0 + _poissons_ratio) +
                      std::sin(3.0 * t) * (1.0 + _poissons_ratio));
 }

◆ execute()

template<bool is_ad>

void InteractionIntegralTempl< is_ad >::execute ( )

overridevirtual

Implements ElementVectorPostprocessor.

Definition at line 516 of file InteractionIntegral.C.

 {
   // calculate phi and dphi for this element
   const std::size_t dim = _current_elem->dim();
   std::unique_ptr<FEBase> fe(FEBase::build(dim, _fe_type));
   fe->attach_quadrature_rule(const_cast<QBase *>(_qrule));
   _phi_curr_elem = &fe->get_phi();
   _dphi_curr_elem = &fe->get_dphi();
   fe->reinit(_current_elem);
 
   // calculate q for all nodes in this element
   std::size_t ring_base = (_q_function_type == QMethod::Topology) ? 0 : 1;
 
   for (auto icfp = beginIndex(_interaction_integral); icfp < _interaction_integral.size(); icfp++)
   {
     _q_curr_elem.clear();
     for (std::size_t i = 0; i < _current_elem->n_nodes(); ++i)
     {
       const Node * this_node = _current_elem->node_ptr(i);
       Real q_this_node;
 
       if (_q_function_type == QMethod::Geometry)
         q_this_node = _crack_front_definition->DomainIntegralQFunction(
             icfp, _ring_index - ring_base, this_node);
       else // QMethod::Topology
         q_this_node = _crack_front_definition->DomainIntegralTopologicalQFunction(
             icfp, _ring_index - ring_base, this_node);
 
       _q_curr_elem.push_back(q_this_node);
     }
 
     for (_qp = 0; _qp < _qrule->n_points(); _qp++)
     {
       Real scalar_q = 0.0;
       RealVectorValue grad_of_scalar_q(0.0, 0.0, 0.0);
 
       for (std::size_t i = 0; i < _current_elem->n_nodes(); ++i)
       {
         scalar_q += (*_phi_curr_elem)[i][_qp] * _q_curr_elem[i];
 
         for (std::size_t j = 0; j < _current_elem->dim(); ++j)
           grad_of_scalar_q(j) += (*_dphi_curr_elem)[i][_qp](j) * _q_curr_elem[i];
       }
 
       if (getBlockCoordSystem() != Moose::COORD_RZ)
         _interaction_integral[icfp] +=
             _JxW[_qp] * _coord[_qp] * computeQpIntegral(icfp, scalar_q, grad_of_scalar_q);
       else
         _interaction_integral[icfp] +=
             _JxW[_qp] * computeQpIntegral(icfp, scalar_q, grad_of_scalar_q);
     }
   }
 }

◆ finalize()

template<bool is_ad>

void InteractionIntegralTempl< is_ad >::finalize ( )

overridevirtual

Reimplemented from ElementVectorPostprocessor.

Definition at line 572 of file InteractionIntegral.C.

 {
   gatherSum(_interaction_integral);
 
   for (auto i = beginIndex(_interaction_integral); i < _interaction_integral.size(); ++i)
   {
     if (_has_symmetry_plane)
       _interaction_integral[i] *= 2.0;
 
     const auto cfp = _crack_front_definition->getCrackFrontPoint(i);
     _x[i] = (*cfp)(0);
     _y[i] = (*cfp)(1);
     _z[i] = (*cfp)(2);
 
     if (_position_type == PositionType::Angle)
       _position[i] = _crack_front_definition->getAngleAlongFront(i);
     else
       _position[i] = _crack_front_definition->getDistanceAlongFront(i);
 
     if (_sif_mode == SifMethod::T && !_treat_as_2d)
       _interaction_integral[i] +=
           _poissons_ratio * _crack_front_definition->getCrackFrontTangentialStrain(i);
 
     _interaction_integral[i] *= _K_factor;
   }
 }

◆ initialize()

template<bool is_ad>

void InteractionIntegralTempl< is_ad >::initialize ( )

overridevirtual

Implements ElementVectorPostprocessor.

Definition at line 253 of file InteractionIntegral.C.

 {
   std::size_t num_pts;
   if (_treat_as_2d && _using_mesh_cutter == false)
     num_pts = 1;
   else
     num_pts = _crack_front_definition->getNumCrackFrontPoints();
 
   _x.assign(num_pts, 0.0);
   _y.assign(num_pts, 0.0);
   _z.assign(num_pts, 0.0);
   _position.assign(num_pts, 0.0);
   _interaction_integral.assign(num_pts, 0.0);
 
   for (const auto * fe_var : _fe_vars)
   {
     if (fe_var->feType() != _fe_type)
       mooseError("All coupled displacements must have the same type");
   }
 }

◆ initialSetup()

template<bool is_ad>

void InteractionIntegralTempl< is_ad >::initialSetup ( )

overridevirtual

Reimplemented from ElementVectorPostprocessor.

Definition at line 245 of file InteractionIntegral.C.

 {
   _treat_as_2d = _crack_front_definition->treatAs2D();
   _using_mesh_cutter = _crack_front_definition->usingMeshCutter();
 }

◆ qFunctionType()

template<bool is_ad>

MooseEnum InteractionIntegralTempl< is_ad >::qFunctionType ( )

static

Get the MooseEnum defining the options for how the q function is defined.

Returns: MooseEnum with options for q function

Definition at line 25 of file InteractionIntegral.C.

Referenced by InteractionIntegralTempl< is_ad >::validParams().

 {
   return MooseEnum("Geometry Topology", "Geometry");
 }

◆ sifModeType()

template<bool is_ad>

MooseEnum InteractionIntegralTempl< is_ad >::sifModeType ( )

static

Get the MooseEnum defining the options for the integral to be computed.

Returns: MooseEnum with options for the integral to be computed

Definition at line 32 of file InteractionIntegral.C.

Referenced by InteractionIntegralTempl< is_ad >::validParams().

 {
   return MooseEnum("KI KII KIII T", "KI");
 }

◆ threadJoin()

template<bool is_ad>

void InteractionIntegralTempl< is_ad >::threadJoin ( const UserObject & y )

overridevirtual

Implements ElementVectorPostprocessor.

Definition at line 601 of file InteractionIntegral.C.

 {
   const auto & uo = static_cast<const InteractionIntegralTempl<is_ad> &>(y);
 
   for (auto i = beginIndex(_interaction_integral); i < _interaction_integral.size(); ++i)
     _interaction_integral[i] += uo._interaction_integral[i];
 }

◆ validParams()

template<bool is_ad>

InputParameters InteractionIntegralTempl< is_ad >::validParams ( )

static

Definition at line 39 of file InteractionIntegral.C.

 {
   InputParameters params = ElementVectorPostprocessor::validParams();
   params.addRequiredCoupledVar(
       "displacements",
       "The displacements appropriate for the simulation geometry and coordinate system");
   params.addCoupledVar("temperature",
                        "The temperature (optional). Must be provided to correctly compute "
                        "stress intensity factors in models with thermal strain gradients.");
   params.addParam<MaterialPropertyName>(
       "functionally_graded_youngs_modulus",
       "Spatially varying elasticity modulus variable. This input is required when "
       "using the functionally graded material capability.");
   params.addParam<MaterialPropertyName>(
       "functionally_graded_youngs_modulus_crack_dir_gradient",
       "Gradient of the spatially varying Young's modulus provided in "
       "'functionally_graded_youngs_modulus' in the direction of crack extension.");
   params.addRequiredParam<UserObjectName>("crack_front_definition",
                                           "The CrackFrontDefinition user object name");
   MooseEnum position_type("Angle Distance", "Distance");
   params.addParam<MooseEnum>(
       "position_type",
       position_type,
       "The method used to calculate position along crack front.  Options are: " +
           position_type.getRawNames());
   params.addRequiredParam<Real>(
       "K_factor", "Conversion factor between interaction integral and stress intensity factor K");
   params.addParam<unsigned int>("symmetry_plane",
                                 "Account for a symmetry plane passing through "
                                 "the plane of the crack, normal to the specified "
                                 "axis (0=x, 1=y, 2=z)");
   params.addRequiredParam<Real>("poissons_ratio", "Poisson's ratio for the material.");
   params.addRequiredParam<Real>("youngs_modulus", "Young's modulus of the material.");
   params.set<bool>("use_displaced_mesh") = false;
   params.addRequiredParam<unsigned int>("ring_index", "Ring ID");
   params.addParam<MooseEnum>("q_function_type",
                              InteractionIntegralTempl<is_ad>::qFunctionType(),
                              "The method used to define the integration domain. Options are: " +
                                  InteractionIntegralTempl<is_ad>::qFunctionType().getRawNames());
   params.addRequiredParam<MooseEnum>(
       "sif_mode",
       InteractionIntegralTempl<is_ad>::sifModeType(),
       "Stress intensity factor to calculate. Choices are: " +
           InteractionIntegralTempl<is_ad>::sifModeType().getRawNames());
   params.addParam<MaterialPropertyName>("eigenstrain_gradient",
                                         "Material defining gradient of eigenstrain tensor");
   params.addParam<MaterialPropertyName>("body_force", "Material defining body force");
   params.addCoupledVar("additional_eigenstrain_00",
                        "Optional additional eigenstrain variable that will be accounted for in the "
                        "interaction integral (component 00 or XX).");
   params.addCoupledVar("additional_eigenstrain_01",
                        "Optional additional eigenstrain variable that will be accounted for in the "
                        "interaction integral (component 01 or XY).");
   params.addCoupledVar("additional_eigenstrain_11",
                        "Optional additional eigenstrain variable that will be accounted for in the "
                        "interaction integral (component 11 or YY).");
   params.addCoupledVar("additional_eigenstrain_22",
                        "Optional additional eigenstrain variable that will be accounted for in the "
                        "interaction integral (component 22 or ZZ).");
   params.addCoupledVar("additional_eigenstrain_02",
                        "Optional additional eigenstrain variable that will be accounted for in the "
                        "interaction integral (component 02 or XZ).");
   params.addCoupledVar("additional_eigenstrain_12",
                        "Optional additional eigenstrain variable that will be accounted for in the "
                        "interaction integral (component 12 or XZ).");
   params.addParamNamesToGroup(
       "additional_eigenstrain_00 additional_eigenstrain_01 additional_eigenstrain_11 "
       "additional_eigenstrain_22 additional_eigenstrain_02 additional_eigenstrain_12",
       "Generic eigenstrains for the computation of the interaction integral.");
   params.addClassDescription(
       "Computes the interaction integral, which is used to compute various "
       "fracture mechanics parameters at a crack tip, including KI, KII, KIII, "
       "and the T stress.");
   return params;
 }

Member Data Documentation

◆ _additional_eigenstrain_gradient_00

template<bool is_ad>

const VariableGradient* InteractionIntegralTempl< is_ad >::_additional_eigenstrain_gradient_00

protected

Gradient used to add contribution to the interaction integral (XX)

Definition at line 131 of file InteractionIntegral.h.

Referenced by InteractionIntegralTempl< is_ad >::InteractionIntegralTempl().

◆ _additional_eigenstrain_gradient_01

template<bool is_ad>

const VariableGradient* InteractionIntegralTempl< is_ad >::_additional_eigenstrain_gradient_01

protected

Gradient used to add contribution to the interaction integral (XY)

Definition at line 133 of file InteractionIntegral.h.

Referenced by InteractionIntegralTempl< is_ad >::InteractionIntegralTempl().

◆ _additional_eigenstrain_gradient_02

template<bool is_ad>

const VariableGradient* InteractionIntegralTempl< is_ad >::_additional_eigenstrain_gradient_02

protected

Gradient used to add contribution to the interaction integral (XZ)

Definition at line 139 of file InteractionIntegral.h.

Referenced by InteractionIntegralTempl< is_ad >::InteractionIntegralTempl().

◆ _additional_eigenstrain_gradient_11

template<bool is_ad>

const VariableGradient* InteractionIntegralTempl< is_ad >::_additional_eigenstrain_gradient_11

protected

Gradient used to add contribution to the interaction integral (YY)

Definition at line 135 of file InteractionIntegral.h.

Referenced by InteractionIntegralTempl< is_ad >::InteractionIntegralTempl().

◆ _additional_eigenstrain_gradient_12

template<bool is_ad>

const VariableGradient* InteractionIntegralTempl< is_ad >::_additional_eigenstrain_gradient_12

protected

Gradient used to add contribution to the interaction integral (YZ)

Definition at line 141 of file InteractionIntegral.h.

Referenced by InteractionIntegralTempl< is_ad >::InteractionIntegralTempl().

◆ _additional_eigenstrain_gradient_22

template<bool is_ad>

const VariableGradient* InteractionIntegralTempl< is_ad >::_additional_eigenstrain_gradient_22

protected

Gradient used to add contribution to the interaction integral (ZZ)

Definition at line 137 of file InteractionIntegral.h.

Referenced by InteractionIntegralTempl< is_ad >::InteractionIntegralTempl().

◆ _body_force

template<bool is_ad>

const MaterialProperty<RealVectorValue>* InteractionIntegralTempl< is_ad >::_body_force

protected

Definition at line 184 of file InteractionIntegral.h.

Referenced by InteractionIntegralTempl< is_ad >::InteractionIntegralTempl().

◆ _crack_front_definition

template<bool is_ad>

const CrackFrontDefinition* const InteractionIntegralTempl< is_ad >::_crack_front_definition

protected

Pointer to the crack front definition object.

Definition at line 89 of file InteractionIntegral.h.

◆ _disp

template<bool is_ad>

std::vector<const VariableValue *> InteractionIntegralTempl< is_ad >::_disp

protected

Displacement variables.

Definition at line 103 of file InteractionIntegral.h.

◆ _dphi_curr_elem

template<bool is_ad>

const std::vector<std::vector<libMesh::RealGradient> >* InteractionIntegralTempl< is_ad >::_dphi_curr_elem

protected

Vector of gradients of shape function values for the current element.

Definition at line 147 of file InteractionIntegral.h.

◆ _eigenstrain_gradient

template<bool is_ad>

const GenericMaterialProperty<RankThreeTensor, is_ad>* InteractionIntegralTempl< is_ad >::_eigenstrain_gradient

protected

Pointers to optionally-used eigenstrain gradient and body force.

Definition at line 183 of file InteractionIntegral.h.

Referenced by InteractionIntegralTempl< is_ad >::InteractionIntegralTempl().

◆ _fe_type

template<bool is_ad>

const libMesh::FEType& InteractionIntegralTempl< is_ad >::_fe_type

protected

FEType object defining order and family of displacement variables.

Definition at line 101 of file InteractionIntegral.h.

◆ _fe_vars

template<bool is_ad>

std::vector<MooseVariableFEBase *> InteractionIntegralTempl< is_ad >::_fe_vars

protected

Vector of all coupled variables.

Definition at line 99 of file InteractionIntegral.h.

◆ _fgm_crack

template<bool is_ad>

const bool InteractionIntegralTempl< is_ad >::_fgm_crack

protected

Whether to consider interaction integral and material properties for a crack in functionally graded material.

Definition at line 123 of file InteractionIntegral.h.

◆ _functionally_graded_youngs_modulus

template<bool is_ad>

const MaterialProperty<Real>* InteractionIntegralTempl< is_ad >::_functionally_graded_youngs_modulus

protected

Spatial elasticity modulus variable for FGM.

Definition at line 113 of file InteractionIntegral.h.

Referenced by InteractionIntegralTempl< is_ad >::InteractionIntegralTempl().

◆ _functionally_graded_youngs_modulus_crack_dir_gradient

template<bool is_ad>

const MaterialProperty<Real>* InteractionIntegralTempl< is_ad >::_functionally_graded_youngs_modulus_crack_dir_gradient

protected

Spatial derivative of the youngs modulus in the crack direction.

Definition at line 111 of file InteractionIntegral.h.

Referenced by InteractionIntegralTempl< is_ad >::InteractionIntegralTempl().

◆ _grad_disp

template<bool is_ad>

std::vector<const VariableGradient *> InteractionIntegralTempl< is_ad >::_grad_disp

protected

Gradient of displacements.

Definition at line 105 of file InteractionIntegral.h.

Referenced by InteractionIntegralTempl< is_ad >::InteractionIntegralTempl().

◆ _grad_temp

template<bool is_ad>

const VariableGradient& InteractionIntegralTempl< is_ad >::_grad_temp

protected

Gradient of temperature.

Definition at line 109 of file InteractionIntegral.h.

◆ _has_additional_eigenstrain

template<bool is_ad>

bool InteractionIntegralTempl< is_ad >::_has_additional_eigenstrain

protected

Whether the user chooses to add other eigenstrain influence (e.g. irradiation-induced)

Definition at line 129 of file InteractionIntegral.h.

Referenced by InteractionIntegralTempl< is_ad >::InteractionIntegralTempl().

◆ _has_symmetry_plane

template<bool is_ad>

bool InteractionIntegralTempl< is_ad >::_has_symmetry_plane

protected

Whether the crack plane is also a symmetry plane in the model.

Definition at line 117 of file InteractionIntegral.h.

◆ _has_temp

template<bool is_ad>

const bool InteractionIntegralTempl< is_ad >::_has_temp

protected

Whether the temperature variable is coupled.

Definition at line 107 of file InteractionIntegral.h.

Referenced by InteractionIntegralTempl< is_ad >::InteractionIntegralTempl().

◆ _interaction_integral

template<bool is_ad>

VectorPostprocessorValue& InteractionIntegralTempl< is_ad >::_interaction_integral

protected

Definition at line 178 of file InteractionIntegral.h.

Referenced by InteractionIntegralTempl< is_ad >::threadJoin().

◆ _K_factor

template<bool is_ad>

Real InteractionIntegralTempl< is_ad >::_K_factor

protected

Conversion factor applied to convert interaction integral to stress intensity factor K.

Definition at line 115 of file InteractionIntegral.h.

◆ _kappa

template<bool is_ad>

Real InteractionIntegralTempl< is_ad >::_kappa

protected

Kappa Lame constant.

Definition at line 149 of file InteractionIntegral.h.

Referenced by InteractionIntegralTempl< is_ad >::InteractionIntegralTempl().

◆ _ndisp

template<bool is_ad>

std::size_t InteractionIntegralTempl< is_ad >::_ndisp

protected

Number of displacement components.

Definition at line 87 of file InteractionIntegral.h.

Referenced by InteractionIntegralTempl< is_ad >::InteractionIntegralTempl().

◆ _phi_curr_elem

template<bool is_ad>

const std::vector<std::vector<Real> >* InteractionIntegralTempl< is_ad >::_phi_curr_elem

protected

Vector of shape function values for the current element.

Definition at line 145 of file InteractionIntegral.h.

◆ _poissons_ratio

template<bool is_ad>

Real InteractionIntegralTempl< is_ad >::_poissons_ratio

protected

Poisson's ratio of the material.

Definition at line 119 of file InteractionIntegral.h.

Referenced by InteractionIntegralTempl< is_ad >::InteractionIntegralTempl().

◆ _position

template<bool is_ad>

VectorPostprocessorValue& InteractionIntegralTempl< is_ad >::_position

protected

Definition at line 177 of file InteractionIntegral.h.

◆ _position_type

template<bool is_ad>

enum InteractionIntegralTempl::PositionType InteractionIntegralTempl< is_ad >::_position_type

protected

◆ _q_curr_elem

template<bool is_ad>

std::vector<Real> InteractionIntegralTempl< is_ad >::_q_curr_elem

protected

Vector of q function values for the nodes in the current element.

Definition at line 143 of file InteractionIntegral.h.

◆ _q_function_type

template<bool is_ad>

enum InteractionIntegralTempl::QMethod InteractionIntegralTempl< is_ad >::_q_function_type

protected

◆ _qp

template<bool is_ad>

unsigned int InteractionIntegralTempl< is_ad >::_qp

protected

Current quadrature point index.

Definition at line 157 of file InteractionIntegral.h.

◆ _r

template<bool is_ad>

Real InteractionIntegralTempl< is_ad >::_r

protected

Radial distance from the current point to the crack front.

Definition at line 153 of file InteractionIntegral.h.

◆ _ring_index

template<bool is_ad>

std::size_t InteractionIntegralTempl< is_ad >::_ring_index

protected

Index of the ring for the integral computed by this object.

Definition at line 125 of file InteractionIntegral.h.

◆ _shear_modulus

template<bool is_ad>

Real InteractionIntegralTempl< is_ad >::_shear_modulus

protected

Shear modulus.

Definition at line 151 of file InteractionIntegral.h.

Referenced by InteractionIntegralTempl< is_ad >::InteractionIntegralTempl().

◆ _sif_mode

template<bool is_ad>

enum InteractionIntegralTempl::SifMethod InteractionIntegralTempl< is_ad >::_sif_mode

protected

◆ _strain

template<bool is_ad>

const GenericMaterialProperty<RankTwoTensor, is_ad>& InteractionIntegralTempl< is_ad >::_strain

protected

Reference to the strain tensor computed by the material models.

Definition at line 97 of file InteractionIntegral.h.

◆ _stress

template<bool is_ad>

const GenericMaterialProperty<RankTwoTensor, is_ad>& InteractionIntegralTempl< is_ad >::_stress

protected

Reference to the stress tensor computed by the material models.

Definition at line 95 of file InteractionIntegral.h.

◆ _theta

template<bool is_ad>

Real InteractionIntegralTempl< is_ad >::_theta

protected

Angle of current point relative to the crack front.

Definition at line 155 of file InteractionIntegral.h.

◆ _total_deigenstrain_dT

template<bool is_ad>

const GenericMaterialProperty<RankTwoTensor, is_ad>* const InteractionIntegralTempl< is_ad >::_total_deigenstrain_dT

protected

Derivative of the total eigenstrain with respect to temperature.

Definition at line 127 of file InteractionIntegral.h.

Referenced by InteractionIntegralTempl< is_ad >::InteractionIntegralTempl().

◆ _treat_as_2d

template<bool is_ad>

bool InteractionIntegralTempl< is_ad >::_treat_as_2d

protected

Whether to treat a model as 2D for computation of fracture integrals.

Definition at line 91 of file InteractionIntegral.h.

◆ _using_mesh_cutter

template<bool is_ad>

bool InteractionIntegralTempl< is_ad >::_using_mesh_cutter

protected

Is the crack defined by an XFEM cutter mesh.

Definition at line 93 of file InteractionIntegral.h.

◆ _x

template<bool is_ad>

VectorPostprocessorValue& InteractionIntegralTempl< is_ad >::_x

protected

Vectors computed by this VectorPostprocessor: x,y,z coordinates, position of nodes along crack front, and interaction integral.

Definition at line 174 of file InteractionIntegral.h.

◆ _y

template<bool is_ad>

VectorPostprocessorValue& InteractionIntegralTempl< is_ad >::_y

protected

Definition at line 175 of file InteractionIntegral.h.

◆ _youngs_modulus

template<bool is_ad>

Real InteractionIntegralTempl< is_ad >::_youngs_modulus

protected

Young's modulus of the material.

Definition at line 121 of file InteractionIntegral.h.

Referenced by InteractionIntegralTempl< is_ad >::InteractionIntegralTempl().

◆ _z

template<bool is_ad>

VectorPostprocessorValue& InteractionIntegralTempl< is_ad >::_z

protected

Definition at line 176 of file InteractionIntegral.h.

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

solid_mechanics/include/vectorpostprocessors/InteractionIntegral.h
solid_mechanics/src/vectorpostprocessors/InteractionIntegral.C

Public Types

Public Member Functions

Static Public Member Functions

Public Attributes

Static Public Attributes

Protected Types

Protected Member Functions

Static Protected Member Functions

Protected Attributes

Static Protected Attributes

Detailed Description

template<bool is_ad> class InteractionIntegralTempl< is_ad >

Member Enumeration Documentation

◆ PositionType

◆ QMethod

◆ SifMethod

Constructor & Destructor Documentation

◆ InteractionIntegralTempl()

Member Function Documentation

◆ computeAuxFields()

◆ computeQpIntegral()

◆ computeTFields()

◆ execute()

◆ finalize()

◆ initialize()

◆ initialSetup()

◆ qFunctionType()

◆ sifModeType()

◆ threadJoin()

◆ validParams()

Member Data Documentation

◆ _additional_eigenstrain_gradient_00

◆ _additional_eigenstrain_gradient_01

◆ _additional_eigenstrain_gradient_02

◆ _additional_eigenstrain_gradient_11

◆ _additional_eigenstrain_gradient_12

◆ _additional_eigenstrain_gradient_22

◆ _body_force

◆ _crack_front_definition

◆ _disp

◆ _dphi_curr_elem

◆ _eigenstrain_gradient

◆ _fe_type

◆ _fe_vars

◆ _fgm_crack

◆ _functionally_graded_youngs_modulus

◆ _functionally_graded_youngs_modulus_crack_dir_gradient

◆ _grad_disp

◆ _grad_temp

◆ _has_additional_eigenstrain

◆ _has_symmetry_plane

◆ _has_temp

◆ _interaction_integral

◆ _K_factor

◆ _kappa

◆ _ndisp

◆ _phi_curr_elem

◆ _poissons_ratio

◆ _position

◆ _position_type

◆ _q_curr_elem

◆ _q_function_type

◆ _qp

◆ _r

◆ _ring_index

◆ _shear_modulus

◆ _sif_mode

◆ _strain

◆ _stress

◆ _theta

◆ _total_deigenstrain_dT

◆ _treat_as_2d

◆ _using_mesh_cutter

◆ _x

◆ _y

◆ _youngs_modulus

◆ _z

template<bool is_ad>
class InteractionIntegralTempl< is_ad >