CrackMeshCut3DUserObject: (1) reads in a mesh describing the crack surface, (2) uses the mesh to do initial cutting of 3D elements, and (3) grows the mesh based on prescribed growth functions. More...

#include <CrackMeshCut3DUserObject.h>

Inheritance diagram for CrackMeshCut3DUserObject:

[legend]

Public Types
typedef DataFileName	DataFileParameterType

Public Member Functions
	CrackMeshCut3DUserObject (const InputParameters &parameters)

virtual void	initialSetup () override

virtual void	initialize () override

virtual const std::vector< Point >	getCrackFrontPoints (unsigned int num_crack_front_points) const override
	get a set of points along a crack front from a XFEM GeometricCutUserObject More...

virtual const std::vector< RealVectorValue >	getCrackPlaneNormals (unsigned int num_crack_front_points) const override
	get a set of normal vectors along a crack front from a XFEM GeometricCutUserObject More...

virtual bool	cutElementByGeometry (const Elem *elem, std::vector< Xfem::CutEdge > &cut_edges, std::vector< Xfem::CutNode > &cut_nodes) const override

virtual bool	cutElementByGeometry (const Elem *elem, std::vector< Xfem::CutFace > &cut_faces) const override

virtual bool	cutFragmentByGeometry (std::vector< std::vector< Point >> &frag_edges, std::vector< Xfem::CutEdge > &cut_edges) const override

virtual bool	cutFragmentByGeometry (std::vector< std::vector< Point >> &frag_faces, std::vector< Xfem::CutFace > &cut_faces) const override

void	findActiveBoundaryNodes ()
	Find all active boundary nodes in the cutter mesh Find boundary nodes that will grow; nodes outside of the structural mesh are inactive. More...

std::vector< int >	getFrontPointsIndex ()
	Get crack front points in the active segment -1 means inactive; positive is the point's index in the Crack Front Definition starting from 0. More...

void	setSubCriticalGrowthSize (std::vector< Real > &growth_size)
	Return growth size at the active boundary to the mesh cutter. More...

unsigned int	getNumberOfCrackFrontPoints () const
	Return the total number of crack front points. More...

MeshBase &	getCutterMesh () const
	Get a reference to the cutter mesh. More...

virtual void	execute () override

virtual void	threadJoin (const UserObject &y) override

virtual void	finalize () override

virtual bool	cutElementByGeometry (const Elem *elem, std::vector< Xfem::CutEdge > &cut_edges, std::vector< Xfem::CutNode > &cut_nodes) const =0
	Check to see whether a specified 2D element should be cut based on geometric conditions. More...

virtual bool	cutElementByGeometry (const Elem *elem, std::vector< Xfem::CutFace > &cut_faces) const =0
	Check to see whether a specified 3D element should be cut based on geometric conditions. More...

virtual bool	cutFragmentByGeometry (std::vector< std::vector< Point >> &frag_edges, std::vector< Xfem::CutEdge > &cut_edges) const =0
	Check to see whether a fragment of a 2D element should be cut based on geometric conditions. More...

virtual bool	cutFragmentByGeometry (std::vector< std::vector< Point >> &frag_faces, std::vector< Xfem::CutFace > &cut_faces) const =0
	Check to see whether a fragment of a 3D element should be cut based on geometric conditions. More...

unsigned int	getInterfaceID () const
	Get the interface ID for this cutting object. More...

void	setInterfaceID (unsigned int interface_id)
	Set the interface ID for this cutting object. More...

bool	shouldHealMesh () const
	Should the elements cut by this cutting object be healed in the current time step? More...

virtual CutSubdomainID	getCutSubdomainID (const Node *) const
	Get CutSubdomainID telling which side the node belongs to relative to the cut. More...

CutSubdomainID	getCutSubdomainID (const Elem *elem) const
	Get the CutSubdomainID for the given element. More...

bool	usesMesh () const
	Getter for if a cutter mesh is used in a derived class. More...

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

bool	isKokkosObject (IsKokkosObjectKey &&) const

MooseApp &	getMooseApp () const

const std::string &	type () const

const std::string &	name () const

std::string	typeAndName () const

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

MooseObjectName	uniqueName () const

const InputParameters &	parameters () const

const hit::Node *	getHitNode () const

bool	hasBase () const

const std::string &	getBase () const

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

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

const T *	queryParam (const std::string &name) const

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

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

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

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

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

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

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

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

std::string	messagePrefix (const bool hit_prefix=true) const

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

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

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

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

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

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

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

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

void	callMooseError (std::string msg, const bool with_prefix, const hit::Node *node=nullptr) const

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

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

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

virtual void	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)

Moose::Kokkos::MaterialProperty< T, dimension >	getKokkosMaterialPropertyByName (const std::string &prop_name)

Moose::Kokkos::MaterialProperty< T, dimension >	getKokkosMaterialPropertyOldByName (const std::string &prop_name)

Moose::Kokkos::MaterialProperty< T, dimension >	getKokkosMaterialPropertyOlderByName (const std::string &prop_name)

Moose::Kokkos::MaterialProperty< T, dimension >	getKokkosMaterialProperty (const std::string &name)

Moose::Kokkos::MaterialProperty< T, dimension >	getKokkosMaterialPropertyOld (const std::string &name)

Moose::Kokkos::MaterialProperty< T, dimension >	getKokkosMaterialPropertyOlder (const std::string &name)

std::pair< const MaterialProperty< T > *, std::set< SubdomainID > >	getBlockMaterialProperty (const MaterialPropertyName &name)

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

const GenericMaterialProperty< T, is_ad > &	getGenericZeroMaterialProperty ()

const GenericMaterialProperty< T, is_ad > &	getGenericZeroMaterialPropertyByName (const std::string &prop_name)

const MaterialProperty< T > &	getZeroMaterialProperty (Ts... args)

std::set< SubdomainID >	getMaterialPropertyBlocks (const std::string &name)

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

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

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

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

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

void	statefulPropertiesAllowed (bool)

virtual bool	getMaterialPropertyCalled () const

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

virtual void	resolveOptionalProperties ()

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

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

Moose::Kokkos::VariableValue	kokkosCoupledVectorTagValueByName (const std::string &var_name, const std::string &tag_name, unsigned int comp=0)

Moose::Kokkos::VariableValue	kokkosCoupledVectorTagValuesByName (const std::string &var_name, const std::string &tag_name)

Moose::Kokkos::VariableGradient	kokkosCoupledVectorTagGradientByName (const std::string &var_name, const std::string &tag_name, unsigned int comp=0)

Moose::Kokkos::VariableGradient	kokkosCoupledVectorTagGradientsByName (const std::string &var_name, const std::string &tag_name)

Moose::Kokkos::VariableValue	kokkosCoupledVectorTagNodalValueByName (const std::string &var_name, const std::string &tag_name, unsigned int comp=0)

Moose::Kokkos::VariableValue	kokkosCoupledVectorTagNodalValuesByName (const std::string &var_name, const std::string &tag_name)

Moose::Kokkos::VariableValue	kokkosCoupledVectorTagValue (const std::string &var_name, const std::string &tag_param_name, unsigned int comp=0)

Moose::Kokkos::VariableValue	kokkosCoupledVectorTagValues (const std::string &var_name, const std::string &tag_param_name)

Moose::Kokkos::VariableGradient	kokkosCoupledVectorTagGradient (const std::string &var_name, const std::string &tag_param_name, unsigned int comp=0)

Moose::Kokkos::VariableGradient	kokkosCoupledVectorTagGradients (const std::string &var_name, const std::string &tag_param_name)

Moose::Kokkos::VariableValue	kokkosCoupledVectorTagNodalValue (const std::string &var_name, const std::string &tag_param_name, unsigned int comp=0)

Moose::Kokkos::VariableValue	kokkosCoupledVectorTagNodalValues (const std::string &var_name, const std::string &tag_param_name)

Moose::Kokkos::VariableValue	kokkosCoupledValue (const std::string &var_name, unsigned int comp=0)

Moose::Kokkos::VariableValue	kokkosCoupledValues (const std::string &var_name)

Moose::Kokkos::VariableGradient	kokkosCoupledGradient (const std::string &var_name, unsigned int comp=0)

Moose::Kokkos::VariableGradient	kokkosCoupledGradients (const std::string &var_name)

Moose::Kokkos::VariableValue	kokkosCoupledNodalValue (const std::string &var_name, unsigned int comp=0)

Moose::Kokkos::VariableValue	kokkosCoupledNodalValues (const std::string &var_name)

Moose::Kokkos::VariableValue	kokkosCoupledValueOld (const std::string &var_name, unsigned int comp=0)

Moose::Kokkos::VariableValue	kokkosCoupledValuesOld (const std::string &var_name)

Moose::Kokkos::VariableGradient	kokkosCoupledGradientOld (const std::string &var_name, unsigned int comp=0)

Moose::Kokkos::VariableGradient	kokkosCoupledGradientsOld (const std::string &var_name)

Moose::Kokkos::VariableValue	kokkosCoupledNodalValueOld (const std::string &var_name, unsigned int comp=0)

Moose::Kokkos::VariableValue	kokkosCoupledNodalValuesOld (const std::string &var_name)

Moose::Kokkos::VariableValue	kokkosCoupledValueOlder (const std::string &var_name, unsigned int comp=0)

Moose::Kokkos::VariableValue	kokkosCoupledValuesOlder (const std::string &var_name)

Moose::Kokkos::VariableGradient	kokkosCoupledGradientOlder (const std::string &var_name, unsigned int comp=0)

Moose::Kokkos::VariableGradient	kokkosCoupledGradientsOlder (const std::string &var_name)

Moose::Kokkos::VariableValue	kokkosCoupledNodalValueOlder (const std::string &var_name, unsigned int comp=0)

Moose::Kokkos::VariableValue	kokkosCoupledNodalValuesOlder (const std::string &var_name)

Moose::Kokkos::VariableValue	kokkosCoupledDot (const std::string &var_name, unsigned int comp=0)

Moose::Kokkos::VariableValue	kokkosCoupledDots (const std::string &var_name)

Moose::Kokkos::VariableValue	kokkosCoupledNodalDot (const std::string &var_name, unsigned int comp=0)

Moose::Kokkos::VariableValue	kokkosCoupledNodalDots (const std::string &var_name)

Moose::Kokkos::Scalar< const Real >	kokkosCoupledDotDu (const std::string &var_name, unsigned int comp=0)

Moose::Kokkos::VariableValue	kokkosZeroValue ()

Moose::Kokkos::VariableGradient	kokkosZeroGradient ()

Moose::Kokkos::VariableValue	kokkosZeroNodalValue ()

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

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

bool	hasKokkosMaterialProperty (const std::string &name)

bool	hasKokkosMaterialPropertyByName (const std::string &name)

bool	hasKokkosMaterialPropertyByName (const std::string &name)

bool	hasGenericMaterialProperty (const std::string &name)

bool	hasGenericMaterialProperty (const std::string &name)

bool	hasGenericMaterialPropertyByName (const std::string &name)

bool	hasGenericMaterialPropertyByName (const std::string &name)

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

const Function &	getFunctionByName (const FunctionName &name) const

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

bool	hasFunctionByName (const FunctionName &name) const

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 void	callMooseError (MooseApp const app, const InputParameters &params, std::string msg, const bool with_prefix, const hit::Node node)

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

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

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

Public Attributes
const ConsoleStream	_console

Static Public Attributes
static const std::string	type_param

static const std::string	name_param

static const std::string	unique_name_param

static const std::string	app_param

static const std::string	moose_base_param

static const std::string	kokkos_object_param

static constexpr PropertyValue::id_type	default_property_id

static constexpr PropertyValue::id_type	zero_property_id

static constexpr auto	SYSTEM

static constexpr auto	NAME

Protected Types
enum	GrowthDirectionEnum { GrowthDirectionEnum::MAX_HOOP_STRESS, GrowthDirectionEnum::FUNCTION }
	Enum to for crack growth direction. More...

enum	GrowthRateEnum { GrowthRateEnum::FATIGUE, GrowthRateEnum::FUNCTION }
	Enum to for crack growth rate. More...

Protected Member Functions
virtual bool	intersectWithEdge (const Point &p1, const Point &p2, const std::vector< Point > &_vertices, Point &point) const
	Check if a line intersects with an element. More...

bool	findIntersection (const Point &p1, const Point &p2, const std::vector< Point > &vertices, Point &point) const
	Find directional intersection along the positive extension of the vector from p1 to p2. More...

bool	isInsideEdge (const Point &p1, const Point &p2, const Point &p) const
	Check if point p is inside the edge p1-p2. More...

Real	getRelativePosition (const Point &p1, const Point &p2, const Point &p) const
	Get the relative position of p from p1. More...

bool	isInsideCutPlane (const std::vector< Point > &_vertices, const Point &p) const
	Check if point p is inside a plane. More...

void	findBoundaryNodes ()
	Find boundary nodes of the cutter mesh This is a simple algorithm simply based on the added angle = 360 degrees Works fine for planar cutting surface for curved cutting surface, need to re-work this subroutine to make it more general. More...

void	findBoundaryEdges ()
	Find boundary edges of the cutter mesh. More...

void	sortBoundaryNodes ()
	Sort boundary nodes to be in the right order along the boundary. More...

Real	findDistance (dof_id_type node1, dof_id_type node2)
	Find distance between two nodes. More...

void	refineBoundary ()
	If boundary nodes are too sparse, add nodes in between. More...

void	findActiveBoundaryDirection ()
	Find growth direction at each active node. More...

void	growFront ()
	Grow the cutter mesh. More...

void	sortFrontNodes ()
	Sort the front nodes. More...

void	findFrontIntersection ()
	Find front-structure intersections. More...

void	refineFront ()
	Refine the mesh at the front. More...

void	triangulation ()
	Create tri3 elements between the new front and the old front. More...

void	joinBoundary ()
	Join active boundaries and inactive boundaries to be the new boundary. 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 ReporterContextBase &	getReporterContextBaseByName (const ReporterName &reporter_name) const

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

void	initializeBlockRestrictable (const MooseObject *moose_object)

void	initializeKokkosBlockRestrictable (const Moose::Kokkos::Mesh *mesh)

Moose::CoordinateSystemType	getBlockCoordSystem ()

KOKKOS_FUNCTION dof_id_type	numKokkosBlockElements () const

KOKKOS_FUNCTION dof_id_type	numKokkosBlockNodes () const

KOKKOS_FUNCTION dof_id_type	numKokkosBlockSides () const

KOKKOS_FUNCTION ContiguousElementID	kokkosBlockElementID (ThreadID tid) const

KOKKOS_FUNCTION ContiguousElementID	kokkosBlockNodeID (ThreadID tid) const

KOKKOS_FUNCTION auto	kokkosBlockElementSideID (ThreadID tid) const

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

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

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

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

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

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

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

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

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

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

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

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

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

bool	hasReporterValueByName (const ReporterName &reporter_name) const

bool	hasReporterValueByName (const ReporterName &reporter_name) const

bool	hasReporterValueByName (const ReporterName &reporter_name) const

bool	hasReporterValueByName (const ReporterName &reporter_name) const

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

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

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

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

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

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


void	serialize (std::string &serialized_buffer)
	Methods to pack/unpack the _marked_elems_2d and _marked_elems_3d data into a structure suitable for parallel communication. More...

void	deserialize (std::vector< std::string > &serialized_buffers)

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

Protected Attributes
const unsigned int	_cut_elem_nnode = 3
	The cutter mesh has triangluar elements only. More...

const unsigned int	_cut_elem_dim = 2

MooseMesh &	_mesh
	The structural mesh. More...

CrackFrontDefinition *	_crack_front_definition
	The crack front definition. More...

std::vector< dof_id_type >	_crack_front_points
	updated crack front definition they are in the same order as defined in the input but the number of nodes may increase its difference from _front is that: _front does not necessarily follow the order of crack front definition therefore, _crack_front_points is generated from _front with the order of crack front definition limitation: this approach does not currently support the growth of one crack front into two More...

const GrowthDirectionEnum	_growth_dir_method
	The direction method for growing mesh at the front. More...

const GrowthRateEnum	_growth_rate_method
	The rate method for growing mesh at the front. More...

const unsigned int	_elem_dim = 3
	The structural mesh must be 3D only. More...

const Real	_const_intersection = 0.01
	Used to define intersection points. More...

Real	_size_control
	Used for cutter mesh refinement and front advancement. More...

unsigned int	_n_step_growth
	Number of steps to grow the mesh. More...

bool	_stop
	Variables to help control the work flow. More...

bool	_grow

std::vector< dof_id_type >	_boundary
	Boundary nodes of the cutter mesh. More...

std::vector< std::vector< dof_id_type > >	_active_boundary
	Active boundary nodes where growth is allowed. More...

std::vector< unsigned int >	_inactive_boundary_pos
	Inactive boundary. More...

std::vector< dof_id_type >	_tracked_crack_front_points
	Front nodes that are grown from the crack front definition defined in the input therefore, they are (1) in the same order as defined in the input and (2) the number of nodes does not change. More...

bool	_cfd

std::set< Xfem::CutEdge >	_boundary_edges
	Edges at the boundary. More...

std::map< dof_id_type, std::vector< dof_id_type > >	_boundary_map
	A map of boundary nodes and their neighbors. More...

std::vector< std::vector< Point > >	_active_direction
	Growth direction for active boundaries. More...

std::vector< Real >	_growth_size
	Growth size for the active boundary in a subcritical simulation. More...

std::vector< unsigned long int >	_dn
	Fatigue life. More...

std::vector< unsigned long int >	_n

std::vector< std::vector< dof_id_type > >	_front
	New boundary after growth. More...

bool	_is_mesh_modified
	Indicator that shows if the cutting mesh is modified or not in this calculation step. More...

unsigned int	_num_crack_front_points
	Total number of crack front points in the mesh cutter. More...

const Function *	_func_x
	Parsed functions of front growth. More...

const Function *	_func_y

const Function *	_func_z

const Function *	_func_v

std::unique_ptr< MeshBase >	_cutter_mesh
	The xfem cutter mesh. More...

std::shared_ptr< XFEM >	_xfem
	Pointer to the XFEM controller object. More...

unsigned int	_interface_id
	Associated interface id. More...

bool	_heal_always
	Heal the mesh. More...

int	_last_step_initialized
	Time step information needed to advance a 3D crack only at the real beginning of a time step. More...

const bool	_uses_mesh
	bool to set if CrackFrontPointsProvider derived objects use a cutter mesh More...

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

Factory &	_factory

ActionFactory &	_action_factory

const std::string &	_type

const std::string &	_name

const InputParameters &	_pars

const ExecFlagEnum &	_execute_enum

const ExecFlagType &	_current_execute_flag

MooseApp &	_restartable_app

const std::string	_restartable_system_name

const THREAD_ID	_restartable_tid

const bool	_restartable_read_only

FEProblemBase &	_mci_feproblem

FEProblemBase &	_mdi_feproblem

MooseApp &	_pg_moose_app

const std::string	_prefix

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 bool	_is_kokkos_object

const Moose::MaterialDataType	_material_data_type

MaterialData &	_material_data

bool	_stateful_allowed

bool	_get_material_property_called

std::vector< std::unique_ptr< PropertyValue > >	_default_properties

std::unordered_set< unsigned int >	_material_property_dependencies

const MaterialPropertyName	_get_suffix

const bool	_use_interpolated_state

const InputParameters &	_c_parameters

const std::string &	_c_name

const std::string &	_c_type

FEProblemBase &	_c_fe_problem

const SystemBase *const	_c_sys

std::unordered_map< std::string, std::vector< MooseVariableFieldBase *> >	_coupled_vars

std::vector< MooseVariableFieldBase *>	_coupled_moose_vars

std::vector< MooseVariable *>	_coupled_standard_moose_vars

std::vector< VectorMooseVariable *>	_coupled_vector_moose_vars

std::vector< ArrayMooseVariable *>	_coupled_array_moose_vars

std::vector< MooseVariableField< Real > *>	_coupled_fv_moose_vars

const std::unordered_map< std::string, std::string > &	_new_to_deprecated_coupled_vars

bool	_c_nodal

bool	_c_is_implicit

const bool	_c_allow_element_to_nodal_coupling

THREAD_ID	_c_tid

std::unordered_map< std::string, std::vector< std::unique_ptr< VariableValue > > >	_default_value

std::unordered_map< std::string, std::unique_ptr< MooseArray< ADReal > > >	_ad_default_value

std::unordered_map< std::string, std::unique_ptr< VectorVariableValue > >	_default_vector_value

std::unordered_map< std::string, std::unique_ptr< ArrayVariableValue > >	_default_array_value

std::unordered_map< std::string, std::unique_ptr< MooseArray< ADRealVectorValue > > >	_ad_default_vector_value

VariableValue	_default_value_zero

VariableGradient	_default_gradient

MooseArray< ADRealVectorValue >	_ad_default_gradient

MooseArray< ADRealTensorValue >	_ad_default_vector_gradient

VariableSecond	_default_second

MooseArray< ADRealTensorValue >	_ad_default_second

MooseArray< ADRealVectorValue >	_ad_default_curl

const VariableValue &	_zero

const VariablePhiValue &	_phi_zero

const MooseArray< ADReal > &	_ad_zero

const VariableGradient &	_grad_zero

const MooseArray< ADRealVectorValue > &	_ad_grad_zero

const VariablePhiGradient &	_grad_phi_zero

const VariableSecond &	_second_zero

const MooseArray< ADRealTensorValue > &	_ad_second_zero

const VariablePhiSecond &	_second_phi_zero

const VectorVariableValue &	_vector_zero

const VectorVariableCurl &	_vector_curl_zero

VectorVariableValue	_default_vector_value_zero

VectorVariableGradient	_default_vector_gradient

VectorVariableCurl	_default_vector_curl

VectorVariableDivergence	_default_div

ArrayVariableValue	_default_array_value_zero

ArrayVariableGradient	_default_array_gradient

bool	_coupleable_neighbor

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 Parallel::Communicator &	_communicator


std::map< unsigned int, std::vector< Xfem::GeomMarkedElemInfo2D > >	_marked_elems_2d
	Containers with information about all 2D and 3D elements marked for cutting by this object. More...

std::map< unsigned int, std::vector< Xfem::GeomMarkedElemInfo3D > >	_marked_elems_3d

Static Protected Attributes
static const std::string	_interpolated_old

static const std::string	_interpolated_older

Detailed Description

CrackMeshCut3DUserObject: (1) reads in a mesh describing the crack surface, (2) uses the mesh to do initial cutting of 3D elements, and (3) grows the mesh based on prescribed growth functions.

Definition at line 25 of file CrackMeshCut3DUserObject.h.

Member Enumeration Documentation

◆ GrowthDirectionEnum

enum CrackMeshCut3DUserObject::GrowthDirectionEnum

strongprotected

Enum to for crack growth direction.

Enumerator
MAX_HOOP_STRESS
FUNCTION

Definition at line 96 of file CrackMeshCut3DUserObject.h.

   {
     MAX_HOOP_STRESS,
     FUNCTION
   };

◆ GrowthRateEnum

enum CrackMeshCut3DUserObject::GrowthRateEnum

strongprotected

Enum to for crack growth rate.

Enumerator
FATIGUE
FUNCTION

Definition at line 105 of file CrackMeshCut3DUserObject.h.

   {
     FATIGUE,
     FUNCTION
   };

Constructor & Destructor Documentation

◆ CrackMeshCut3DUserObject()

CrackMeshCut3DUserObject::CrackMeshCut3DUserObject ( const InputParameters & parameters )

Definition at line 54 of file CrackMeshCut3DUserObject.C.

   : MeshCutUserObjectBase(parameters),
     _mesh(_subproblem.mesh()),
     _growth_dir_method(getParam<MooseEnum>("growth_dir_method").getEnum<GrowthDirectionEnum>()),
     _growth_rate_method(getParam<MooseEnum>("growth_rate_method").getEnum<GrowthRateEnum>()),
     _n_step_growth(getParam<unsigned int>("n_step_growth")),
     _is_mesh_modified(false),
     _func_x(parameters.isParamValid("growth_direction_x") ? &getFunction("growth_direction_x")
                                                           : nullptr),
     _func_y(parameters.isParamValid("growth_direction_y") ? &getFunction("growth_direction_y")
                                                           : nullptr),
     _func_z(parameters.isParamValid("growth_direction_z") ? &getFunction("growth_direction_z")
                                                           : nullptr),
     _func_v(parameters.isParamValid("growth_rate") ? &getFunction("growth_rate") : nullptr)
 {
   _grow = (_n_step_growth == 0 ? 0 : 1);
 
   if (_grow)
   {
     if (!isParamValid("size_control"))
       mooseError("Crack growth needs size control");
 
     _size_control = getParam<Real>("size_control");
 
     if (_growth_dir_method == GrowthDirectionEnum::FUNCTION &&
         (_func_x == nullptr || _func_y == nullptr || _func_z == nullptr))
       mooseError("function is not specified for the function method that defines growth direction");
 
     if (_growth_dir_method == GrowthDirectionEnum::FUNCTION && _func_v == nullptr)
       mooseError("function is not specified for the function method that defines growth rate");
 
     if (_growth_dir_method == GrowthDirectionEnum::FUNCTION && _func_v == nullptr)
       mooseError("function with a variable is not specified for the fatigue method that defines "
                  "growth rate");
 
     if (isParamValid("crack_front_nodes"))
     {
       _tracked_crack_front_points = getParam<std::vector<dof_id_type>>("crack_front_nodes");
       _num_crack_front_points = _tracked_crack_front_points.size();
       _cfd = true;
     }
     else
       _cfd = false;
   }
 
   if ((_growth_dir_method == GrowthDirectionEnum::MAX_HOOP_STRESS ||
        _growth_rate_method == GrowthRateEnum::FATIGUE) &&
       !_cfd)
     mooseError("'crack_front_nodes' is not specified to use crack growth criteria!");
 
   // test element type; only tri3 elements are allowed
   for (const auto & cut_elem : _cutter_mesh->element_ptr_range())
   {
     if (cut_elem->n_nodes() != _cut_elem_nnode)
       mooseError("The input cut mesh should include tri elements only!");
     if (cut_elem->dim() != _cut_elem_dim)
       mooseError("The input cut mesh should have 2D elements only!");
   }
 }

Member Function Documentation

◆ cutElementByGeometry() [1/4]

bool CrackMeshCut3DUserObject::cutElementByGeometry	(	const Elem *	elem,
		std::vector< Xfem::CutEdge > &	cut_edges,
		std::vector< Xfem::CutNode > &	cut_nodes
	)		const

overridevirtual

Definition at line 177 of file CrackMeshCut3DUserObject.C.

 {
   mooseError("invalid method for 3D mesh cutting");
   return false;
 }

◆ cutElementByGeometry() [2/4]

bool CrackMeshCut3DUserObject::cutElementByGeometry	(	const Elem *	elem,
		std::vector< Xfem::CutFace > &	cut_faces
	)		const

overridevirtual

Definition at line 186 of file CrackMeshCut3DUserObject.C.

 {
   bool elem_cut = false;
 
   if (elem->dim() != _elem_dim)
     mooseError("The structural mesh to be cut by a surface mesh must be 3D!");
 
   for (unsigned int i = 0; i < elem->n_sides(); ++i)
   {
     // This returns the lowest-order type of side.
     std::unique_ptr<const Elem> curr_side = elem->side_ptr(i);
     if (curr_side->dim() != 2)
       mooseError("In cutElementByGeometry dimension of side must be 2, but it is ",
                  curr_side->dim());
     unsigned int n_edges = curr_side->n_sides();
 
     std::vector<unsigned int> cut_edges;
     std::vector<Real> cut_pos;
 
     for (unsigned int j = 0; j < n_edges; j++)
     {
       // This returns the lowest-order type of side.
       std::unique_ptr<const Elem> curr_edge = curr_side->side_ptr(j);
       if (curr_edge->type() != EDGE2)
         mooseError("In cutElementByGeometry face edge must be EDGE2, but type is: ",
                    libMesh::Utility::enum_to_string(curr_edge->type()),
                    " base element type is: ",
                    libMesh::Utility::enum_to_string(elem->type()));
       const Node * node1 = curr_edge->node_ptr(0);
       const Node * node2 = curr_edge->node_ptr(1);
 
       for (const auto & cut_elem : _cutter_mesh->element_ptr_range())
       {
         std::vector<Point> vertices;
 
         for (auto & node : cut_elem->node_ref_range())
         {
           Point & this_point = node;
           vertices.push_back(this_point);
         }
 
         Point intersection;
         if (intersectWithEdge(*node1, *node2, vertices, intersection))
         {
           cut_edges.push_back(j);
           cut_pos.emplace_back(getRelativePosition(*node1, *node2, intersection));
         }
       }
     }
 
     // if two edges of an element are cut, it is considered as an element being cut
     if (cut_edges.size() == 2)
     {
       elem_cut = true;
       Xfem::CutFace mycut;
       mycut._face_id = i;
       mycut._face_edge.push_back(cut_edges[0]);
       mycut._face_edge.push_back(cut_edges[1]);
       mycut._position.push_back(cut_pos[0]);
       mycut._position.push_back(cut_pos[1]);
       cut_faces.push_back(mycut);
     }
   }
   return elem_cut;
 }

◆ cutElementByGeometry() [3/4]

virtual bool GeometricCutUserObject::cutElementByGeometry	(	const Elem *	elem,
		std::vector< Xfem::CutEdge > &	cut_edges,
		std::vector< Xfem::CutNode > &	cut_nodes
	)		const

pure virtualinherited

Check to see whether a specified 2D element should be cut based on geometric conditions.

Parameters

elem	Pointer to the libMesh element to be considered for cutting
cut_edges	Data structure filled with information about edges to be cut
cut_nodes	Data structure filled with information about nodes to be cut

Returns: bool true if element is to be cut

Implemented in GeometricCut3DUserObject.

Referenced by GeometricCutUserObject::execute().

◆ cutElementByGeometry() [4/4]

virtual bool GeometricCutUserObject::cutElementByGeometry	(	const Elem *	elem,
		std::vector< Xfem::CutFace > &	cut_faces
	)		const

pure virtualinherited

Check to see whether a specified 3D element should be cut based on geometric conditions.

Parameters

elem	Pointer to the libMesh element to be considered for cutting
cut_faces	Data structure filled with information about edges to be cut

Returns: bool true if element is to be cut

Implemented in GeometricCut3DUserObject.

◆ cutFragmentByGeometry() [1/4]

bool CrackMeshCut3DUserObject::cutFragmentByGeometry	(	std::vector< std::vector< Point >> &	frag_edges,
		std::vector< Xfem::CutEdge > &	cut_edges
	)		const

overridevirtual

Definition at line 257 of file CrackMeshCut3DUserObject.C.

 {
   mooseError("invalid method for 3D mesh cutting");
   return false;
 }

◆ cutFragmentByGeometry() [2/4]

bool CrackMeshCut3DUserObject::cutFragmentByGeometry	(	std::vector< std::vector< Point >> &	frag_faces,
		std::vector< Xfem::CutFace > &	cut_faces
	)		const

overridevirtual

Definition at line 265 of file CrackMeshCut3DUserObject.C.

 {
   // TODO: Need this for branching in 3D
   mooseError("cutFragmentByGeometry not yet implemented for 3D mesh cutting");
   return false;
 }

◆ cutFragmentByGeometry() [3/4]

virtual bool GeometricCutUserObject::cutFragmentByGeometry	(	std::vector< std::vector< Point >> &	frag_edges,
		std::vector< Xfem::CutEdge > &	cut_edges
	)		const

pure virtualinherited

Check to see whether a fragment of a 2D element should be cut based on geometric conditions.

Parameters

frag_edges	Data structure defining the current fragment to be considered
cut_edges	Data structure filled with information about fragment edges to be cut

Returns: bool true if fragment is to be cut

Implemented in GeometricCut3DUserObject.

Referenced by GeometricCutUserObject::execute().

◆ cutFragmentByGeometry() [4/4]

virtual bool GeometricCutUserObject::cutFragmentByGeometry	(	std::vector< std::vector< Point >> &	frag_faces,
		std::vector< Xfem::CutFace > &	cut_faces
	)		const

pure virtualinherited

Check to see whether a fragment of a 3D element should be cut based on geometric conditions.

Parameters

frag_faces	Data structure defining the current fragment to be considered
cut_faces	Data structure filled with information about fragment faces to be cut

Returns: bool true if fragment is to be cut

Implemented in GeometricCut3DUserObject.

◆ deserialize()

void GeometricCutUserObject::deserialize ( std::vector< std::string > & serialized_buffers )

protectedinherited

Definition at line 219 of file GeometricCutUserObject.C.

Referenced by GeometricCutUserObject::finalize().

 {
   mooseAssert(serialized_buffers.size() == _app.n_processors(),
               "Unexpected size of serialized_buffers: " << serialized_buffers.size());
 
   // The input string stream used for deserialization
   std::istringstream iss;
 
   // Loop over all datastructures for all processors to perfrom the gather operation
   for (unsigned int rank = 0; rank < serialized_buffers.size(); ++rank)
   {
     // skip the current processor (its data is already in the structures)
     if (rank == processor_id())
       continue;
 
     // populate the stream with a new buffer and reset stream state
     iss.clear();
     iss.str(serialized_buffers[rank]);
 
     // Load the communicated data into temporary structures
     std::map<unsigned int, std::vector<Xfem::GeomMarkedElemInfo2D>> other_marked_elems_2d;
     std::map<unsigned int, std::vector<Xfem::GeomMarkedElemInfo3D>> other_marked_elems_3d;
     dataLoad(iss, other_marked_elems_2d, this);
     dataLoad(iss, other_marked_elems_3d, this);
 
     // merge the data in with the current processor's data
     _marked_elems_2d.insert(other_marked_elems_2d.begin(), other_marked_elems_2d.end());
     _marked_elems_3d.insert(other_marked_elems_3d.begin(), other_marked_elems_3d.end());
   }
 }

◆ execute()

void GeometricCutUserObject::execute ( )

overridevirtualinherited

Implements ElementUserObject.

Definition at line 64 of file GeometricCutUserObject.C.

 {
   if (_current_elem->dim() == 2)
   {
     std::vector<Xfem::CutEdge> elem_cut_edges;
     std::vector<Xfem::CutNode> elem_cut_nodes;
     std::vector<Xfem::CutEdge> frag_cut_edges;
     std::vector<std::vector<Point>> frag_edges;
 
     EFAElement2D * EFAElem = _xfem->getEFAElem2D(_current_elem);
 
     // Don't cut again if elem has been already cut twice
     if (!EFAElem->isFinalCut())
     {
       // get fragment edges
       _xfem->getFragmentEdges(_current_elem, EFAElem, frag_edges);
 
       // mark cut edges for the element and its fragment
       bool cut = cutElementByGeometry(_current_elem, elem_cut_edges, elem_cut_nodes);
       if (EFAElem->numFragments() > 0)
         cut |= cutFragmentByGeometry(frag_edges, frag_cut_edges);
 
       if (cut)
       {
         Xfem::GeomMarkedElemInfo2D gmei2d;
         gmei2d._elem_cut_edges = elem_cut_edges;
         gmei2d._elem_cut_nodes = elem_cut_nodes;
         gmei2d._frag_cut_edges = frag_cut_edges;
         gmei2d._frag_edges = frag_edges;
         _marked_elems_2d[_current_elem->id()].push_back(gmei2d);
       }
     }
   }
   else if (_current_elem->dim() == 3)
   {
     std::vector<Xfem::CutFace> elem_cut_faces;
     std::vector<Xfem::CutFace> frag_cut_faces;
     std::vector<std::vector<Point>> frag_faces;
 
     EFAElement3D * EFAElem = _xfem->getEFAElem3D(_current_elem);
 
     // Don't cut again if elem has been already cut twice
     if (!EFAElem->isFinalCut())
     {
       // get fragment edges
       _xfem->getFragmentFaces(_current_elem, EFAElem, frag_faces);
 
       // mark cut faces for the element and its fragment
       bool cut = cutElementByGeometry(_current_elem, elem_cut_faces);
       // TODO: This would be done for branching, which is not yet supported in 3D
       // if (EFAElem->numFragments() > 0)
       //  cut |= cutFragmentByGeometry(frag_faces, frag_cut_faces, _t);
 
       if (cut)
       {
         Xfem::GeomMarkedElemInfo3D gmei3d;
         gmei3d._elem_cut_faces = elem_cut_faces;
         gmei3d._frag_cut_faces = frag_cut_faces;
         gmei3d._frag_faces = frag_faces;
         _marked_elems_3d[_current_elem->id()].push_back(gmei3d);
       }
     }
   }
 }

◆ finalize()

void GeometricCutUserObject::finalize ( )

overridevirtualinherited

Implements ElementUserObject.

Definition at line 251 of file GeometricCutUserObject.C.

 {
   // for single processor runs we do not need to do anything here
   if (_app.n_processors() > 1)
   {
     // create send buffer
     std::string send_buffer;
 
     // create byte buffers for the streams received from all processors
     std::vector<std::string> recv_buffers;
 
     // pack the complex datastructures into the string stream
     serialize(send_buffer);
 
     // broadcast serialized data to and receive from all processors
     _communicator.allgather(send_buffer, recv_buffers);
 
     // unpack the received data and merge it into the local data structures
     deserialize(recv_buffers);
   }
 
   for (const auto & it : _marked_elems_2d)
     for (const auto & gmei : it.second)
       _xfem->addGeomMarkedElem2D(it.first, gmei, _interface_id);
 
   for (const auto & it : _marked_elems_3d)
     for (const auto & gmei : it.second)
       _xfem->addGeomMarkedElem3D(it.first, gmei, _interface_id);
 
   _marked_elems_2d.clear();
   _marked_elems_3d.clear();
 }

◆ findActiveBoundaryDirection()

void CrackMeshCut3DUserObject::findActiveBoundaryDirection ( )

protected

Find growth direction at each active node.

Definition at line 700 of file CrackMeshCut3DUserObject.C.

Referenced by initialize().

 {
   mooseAssert(!(_cfd && _active_boundary.size() != 1),
               "crack-front-definition using the cutter mesh only supports one active crack front "
               "segment for now");
 
   _active_direction.clear();
 
   for (unsigned int i = 0; i < _active_boundary.size(); ++i)
   {
     std::vector<Point> temp;
     Point dir;
 
     if (_inactive_boundary_pos.size() != 0)
     {
       for (unsigned int j = 0; j < 3; ++j)
         dir(j) = 0;
       temp.push_back(dir);
     }
 
     unsigned int i1 = 1;
     unsigned int i2 = _active_boundary[i].size() - 1;
     if (_inactive_boundary_pos.size() == 0)
     {
       i1 = 0;
       i2 = _active_boundary[i].size();
     }
 
     if (_growth_dir_method == GrowthDirectionEnum::FUNCTION)
       // loop over active front points
       for (unsigned int j = i1; j < i2; ++j)
       {
         Node * this_node = _cutter_mesh->node_ptr(_active_boundary[i][j]);
         mooseAssert(this_node, "Node is NULL");
         Point & this_point = *this_node;
         dir(0) = _func_x->value(0, this_point);
         dir(1) = _func_y->value(0, this_point);
         dir(2) = _func_z->value(0, this_point);
 
         temp.push_back(dir);
       }
     // determine growth direction based on KI and KII at the crack front
     else if (_growth_dir_method == GrowthDirectionEnum::MAX_HOOP_STRESS)
     {
       const VectorPostprocessorValue & k1 = getVectorPostprocessorValueByName("II_KI_1", "II_KI_1");
       const VectorPostprocessorValue & k2 =
           getVectorPostprocessorValueByName("II_KII_1", "II_KII_1");
       mooseAssert(k1.size() == k2.size(), "KI and KII VPPs should have the same size");
       mooseAssert(k1.size() == _active_boundary[0].size(),
                   "the number of crack front nodes in the self-similar method should equal to the "
                   "size of VPP defined at the crack front");
       mooseAssert(_crack_front_points.size() == _active_boundary[0].size(),
                   "the number of crack front nodes should be the same in _crack_front_points and "
                   "_active_boundary[0]");
 
       // the node order in _active_boundary[0] and _crack_front_points may be the same or opposite,
       // their correspondence is needed
       std::vector<int> index = getFrontPointsIndex();
 
       for (unsigned int j = i1; j < i2; ++j)
       {
         int ind = index[j];
         Real theta = 2 * std::atan((k1[ind] - std::sqrt(k1[ind] * k1[ind] + k2[ind] * k2[ind])) /
                                    (4 * k2[ind]));
 
         // growth direction in crack front coord (cfc) system based on the max hoop stress criterion
         RealVectorValue dir_cfc;
 
         // growth direction in global coord system based on the max hoop stress criterion
         RealVectorValue dir;
 
         dir_cfc(0) = std::cos(theta);
         dir_cfc(1) = std::sin(theta);
         dir_cfc(2) = 0;
         dir = _crack_front_definition->rotateFromCrackFrontCoordsToGlobal(dir_cfc, ind);
 
         temp.push_back(dir);
       }
     }
     else
       mooseError("This growth_dir_method is not pre-defined!");
 
     if (_inactive_boundary_pos.size() != 0)
     {
       for (unsigned int j = 0; j < 3; ++j)
         dir(j) = 0;
       temp.push_back(dir);
     }
 
     _active_direction.push_back(temp);
   }
 
   // normalize the directional vector
   Real maxl = 0;
 
   for (unsigned int i = 0; i < _active_direction.size(); ++i)
     for (unsigned int j = 0; j < _active_direction[i].size(); ++j)
     {
       Point pt = _active_direction[i][j];
       Real length = std::sqrt(pt * pt);
       if (length > maxl)
         maxl = length;
     }
 
   for (unsigned int i = 0; i < _active_direction.size(); ++i)
     for (unsigned int j = 0; j < _active_direction[i].size(); ++j)
       _active_direction[i][j] /= maxl;
 }

◆ findActiveBoundaryNodes()

void CrackMeshCut3DUserObject::findActiveBoundaryNodes ( )

Find all active boundary nodes in the cutter mesh Find boundary nodes that will grow; nodes outside of the structural mesh are inactive.

Definition at line 640 of file CrackMeshCut3DUserObject.C.

Referenced by ParisLaw::execute(), and initialize().

 {
   _active_boundary.clear();
   _inactive_boundary_pos.clear();
 
   std::unique_ptr<PointLocatorBase> pl = _mesh.getPointLocator();
   pl->enable_out_of_mesh_mode();
 
   unsigned int n_boundary = _boundary.size();
 
   // if the node is outside of the structural model, store its position in _boundary to
   // _inactive_boundary_pos
   for (unsigned int j = 0; j < n_boundary; ++j)
   {
     Node * this_node = _cutter_mesh->node_ptr(_boundary[j]);
     mooseAssert(this_node, "Node is NULL");
     Point & this_point = *this_node;
 
     const Elem * elem = (*pl)(this_point);
     if (elem == nullptr)
       _inactive_boundary_pos.push_back(j);
   }
 
   unsigned int n_inactive_boundary = _inactive_boundary_pos.size();
 
   // all nodes are inactive, stop
   if (n_inactive_boundary == n_boundary)
     _stop = 1;
 
   // find and store active boundary segments in "_active_boundary"
   if (n_inactive_boundary == 0)
     _active_boundary.push_back(_boundary);
   else
   {
     for (unsigned int i = 0; i < n_inactive_boundary - 1; ++i)
     {
       if (_inactive_boundary_pos[i + 1] - _inactive_boundary_pos[i] != 1)
       {
         std::vector<dof_id_type> temp;
         for (unsigned int j = _inactive_boundary_pos[i]; j <= _inactive_boundary_pos[i + 1]; ++j)
         {
           temp.push_back(_boundary[j]);
         }
         _active_boundary.push_back(temp);
       }
     }
     if (_inactive_boundary_pos[n_inactive_boundary - 1] - _inactive_boundary_pos[0] <
         n_boundary - 1)
     {
       std::vector<dof_id_type> temp;
       for (unsigned int j = _inactive_boundary_pos[n_inactive_boundary - 1]; j < n_boundary; ++j)
         temp.push_back(_boundary[j]);
       for (unsigned int j = 0; j <= _inactive_boundary_pos[0]; ++j)
         temp.push_back(_boundary[j]);
       _active_boundary.push_back(temp);
     }
   }
 }

◆ findBoundaryEdges()

void CrackMeshCut3DUserObject::findBoundaryEdges ( )

protected

Find boundary edges of the cutter mesh.

Definition at line 395 of file CrackMeshCut3DUserObject.C.

Referenced by initialSetup().

 {
   _boundary_edges.clear();
 
   std::vector<dof_id_type> corner_elem_id;
   unsigned int counter = 0;
 
   std::vector<dof_id_type> node_id(_cut_elem_nnode);
   std::vector<bool> is_node_on_boundary(_cut_elem_nnode);
 
   for (const auto & cut_elem : _cutter_mesh->element_ptr_range())
   {
     for (unsigned int i = 0; i < _cut_elem_nnode; ++i)
     {
       node_id[i] = cut_elem->node_ptr(i)->id();
       is_node_on_boundary[i] = (_boundary_map.find(node_id[i]) != _boundary_map.end());
     }
 
     if (is_node_on_boundary[0] && is_node_on_boundary[1] && is_node_on_boundary[2])
     {
       // this is an element at the corner; all nodes are on the boundary but not all edges are on
       // the boundary
       corner_elem_id.push_back(counter);
     }
     else
     {
       // for other elements, find and store boundary edges
       for (unsigned int i = 0; i < _cut_elem_nnode; ++i)
       {
         // if both nodes on an edge are on the boundary, it is a boundary edge.
         if (is_node_on_boundary[i] && is_node_on_boundary[(i + 1 <= 2) ? i + 1 : 0])
         {
           dof_id_type node1 = node_id[i];
           dof_id_type node2 = node_id[(i + 1 <= 2) ? i + 1 : 0];
           if (node1 > node2)
             std::swap(node1, node2);
 
           Xfem::CutEdge ce;
 
           if (node1 > node2)
             std::swap(node1, node2);
           ce._id1 = node1;
           ce._id2 = node2;
 
           _boundary_edges.insert(ce);
         }
       }
     }
     ++counter;
   }
 
   // loop over edges in corner elements
   // if an edge is shared by two elements, it is not an boundary edge (is_edge_inside = 1)
   for (unsigned int i = 0; i < corner_elem_id.size(); ++i)
   {
     auto elem_it = _cutter_mesh->elements_begin();
 
     for (dof_id_type j = 0; j < corner_elem_id[i]; ++j)
       ++elem_it;
     Elem * cut_elem = *elem_it;
 
     for (unsigned int j = 0; j < _cut_elem_nnode; ++j)
     {
       bool is_edge_inside = 0;
 
       dof_id_type node1 = cut_elem->node_ptr(j)->id();
       dof_id_type node2 = cut_elem->node_ptr((j + 1 <= 2) ? j + 1 : 0)->id();
       if (node1 > node2)
         std::swap(node1, node2);
 
       unsigned int counter = 0;
       for (const auto & cut_elem2 : _cutter_mesh->element_ptr_range())
       {
         if (counter != corner_elem_id[i])
         {
           for (unsigned int k = 0; k < _cut_elem_nnode; ++k)
           {
             dof_id_type node3 = cut_elem2->node_ptr(k)->id();
             dof_id_type node4 = cut_elem2->node_ptr((k + 1 <= 2) ? k + 1 : 0)->id();
             if (node3 > node4)
               std::swap(node3, node4);
 
             if (node1 == node3 && node2 == node4)
             {
               is_edge_inside = 1;
               goto endloop;
             }
           }
         }
         ++counter;
       }
     endloop:
       if (is_edge_inside == 0)
       {
         // store boundary edges
         Xfem::CutEdge ce;
 
         if (node1 > node2)
           std::swap(node1, node2);
         ce._id1 = node1;
         ce._id2 = node2;
 
         _boundary_edges.insert(ce);
       }
       else
       {
         // this is not a boundary edge; remove it from existing edge list
         for (auto it = _boundary_edges.begin(); it != _boundary_edges.end();)
         {
           if ((*it)._id1 == node1 && (*it)._id2 == node2)
             it = _boundary_edges.erase(it);
           else
             ++it;
         }
       }
     }
   }
 }

◆ findBoundaryNodes()

void CrackMeshCut3DUserObject::findBoundaryNodes ( )

protected

Find boundary nodes of the cutter mesh This is a simple algorithm simply based on the added angle = 360 degrees Works fine for planar cutting surface for curved cutting surface, need to re-work this subroutine to make it more general.

Definition at line 383 of file CrackMeshCut3DUserObject.C.

Referenced by initialSetup().

 {
   auto boundary_node_ids = MeshTools::find_boundary_nodes(*_cutter_mesh);
   for (auto it = boundary_node_ids.cbegin(); it != boundary_node_ids.cend(); it++)
   {
     dof_id_type id = *it;
     std::vector<dof_id_type> neighbors;
     _boundary_map[id] = neighbors;
   }
 }

◆ findDistance()

Real CrackMeshCut3DUserObject::findDistance	(	dof_id_type	node1,
		dof_id_type	node2
	)

protected

Find distance between two nodes.

Definition at line 575 of file CrackMeshCut3DUserObject.C.

Referenced by refineBoundary(), refineFront(), and triangulation().

 {
   Node * n1 = _cutter_mesh->node_ptr(node1);
   mooseAssert(n1 != nullptr, "Node is NULL");
   Node * n2 = _cutter_mesh->node_ptr(node2);
   mooseAssert(n2 != nullptr, "Node is NULL");
   Real distance = (*n1 - *n2).norm();
   return distance;
 }

◆ findFrontIntersection()

void CrackMeshCut3DUserObject::findFrontIntersection ( )

protected

Find front-structure intersections.

Definition at line 888 of file CrackMeshCut3DUserObject.C.

Referenced by initialize().

 {
   ConstBndElemRange & range = *_mesh.getBoundaryElementRange();
 
   for (unsigned int i = 0; i < _front.size(); ++i)
   {
     if (_front[i].size() >= 2)
     {
       std::vector<Point> pint1;
       std::vector<Point> pint2;
       std::vector<Real> length1;
       std::vector<Real> length2;
 
       Real node_id = _front[i][0];
       Node * this_node = _cutter_mesh->node_ptr(node_id);
       mooseAssert(this_node, "Node is NULL");
       Point & p2 = *this_node;
 
       if (_front[i].size() >= 4)
         node_id = _front[i][2];
       else
         node_id = _front[i][1];
 
       this_node = _cutter_mesh->node_ptr(node_id);
       mooseAssert(this_node, "Node is NULL");
       Point & p1 = *this_node;
 
       node_id = _front[i].back();
       this_node = _cutter_mesh->node_ptr(node_id);
       mooseAssert(this_node, "Node is NULL");
       Point & p4 = *this_node;
 
       if (_front[i].size() >= 4)
         node_id = _front[i][_front[i].size() - 3];
       else
         node_id = _front[i][_front[i].size() - 2];
 
       this_node = _cutter_mesh->node_ptr(node_id);
       mooseAssert(this_node, "Node is NULL");
       Point & p3 = *this_node;
 
       bool do_inter1 = 1;
       bool do_inter2 = 1;
 
       std::unique_ptr<PointLocatorBase> pl = _mesh.getPointLocator();
       pl->enable_out_of_mesh_mode();
       const Elem * elem = (*pl)(p1);
       if (elem == nullptr)
         do_inter1 = 0;
       elem = (*pl)(p4);
       if (elem == nullptr)
         do_inter2 = 0;
 
       for (const auto & belem : range)
       {
         Point pt;
         std::vector<Point> vertices;
 
         elem = belem->_elem;
         std::unique_ptr<const Elem> curr_side = elem->side_ptr(belem->_side);
         for (unsigned int j = 0; j < curr_side->n_nodes(); ++j)
         {
           const Node * node = curr_side->node_ptr(j);
           const Point & this_point = *node;
           vertices.push_back(this_point);
         }
 
         if (findIntersection(p1, p2, vertices, pt))
         {
           pint1.push_back(pt);
           length1.push_back((pt - p1) * (pt - p1));
         }
         if (findIntersection(p3, p4, vertices, pt))
         {
           pint2.push_back(pt);
           length2.push_back((pt - p3) * (pt - p3));
         }
       }
 
       if (length1.size() != 0 && do_inter1)
       {
         auto it1 = std::min_element(length1.begin(), length1.end());
         Point inter1 = pint1[std::distance(length1.begin(), it1)];
         inter1 += (inter1 - p1) * _const_intersection;
 
         Node * this_node = Node::build(inter1, _cutter_mesh->n_nodes()).release();
         _cutter_mesh->add_node(this_node);
 
         mooseAssert(_cutter_mesh->n_nodes() - 1 > 0,
                     "The cut mesh should have at least one element.");
         unsigned int n = _cutter_mesh->n_nodes() - 1;
 
         auto it = _front[i].begin();
         _front[i].insert(it, n);
 
         if (_cfd)
           _tracked_crack_front_points[_tracked_crack_front_points.size() - 1] = n;
       }
 
       if (length2.size() != 0 && do_inter2)
       {
         auto it2 = std::min_element(length2.begin(), length2.end());
         Point inter2 = pint2[std::distance(length2.begin(), it2)];
         inter2 += (inter2 - p2) * _const_intersection;
 
         Node * this_node = Node::build(inter2, _cutter_mesh->n_nodes()).release();
         _cutter_mesh->add_node(this_node);
 
         dof_id_type n = _cutter_mesh->n_nodes() - 1;
 
         auto it = _front[i].begin();
         unsigned int m = _front[i].size();
         _front[i].insert(it + m, n);
 
         if (_cfd)
           _tracked_crack_front_points[0] = n;
       }
     }
   }
 }

◆ findIntersection()

bool CrackMeshCut3DUserObject::findIntersection	(	const Point &	p1,
		const Point &	p2,
		const std::vector< Point > &	vertices,
		Point &	point
	)		const

protected

Find directional intersection along the positive extension of the vector from p1 to p2.

Definition at line 305 of file CrackMeshCut3DUserObject.C.

Referenced by findFrontIntersection().

 {
   bool has_intersection = false;
 
   Plane elem_plane(vertices[0], vertices[1], vertices[2]);
   Point point = vertices[0];
   Point normal = elem_plane.unit_normal(point);
 
   std::array<Real, 3> plane_point = {{point(0), point(1), point(2)}};
   std::array<Real, 3> planenormal = {{normal(0), normal(1), normal(2)}};
   std::array<Real, 3> p_begin = {{p1(0), p1(1), p1(2)}};
   std::array<Real, 3> p_end = {{p2(0), p2(1), p2(2)}};
   std::array<Real, 3> cut_point = {{0.0, 0.0, 0.0}};
 
   if (Xfem::plane_normal_line_exp_int_3d(
           &plane_point[0], &planenormal[0], &p_begin[0], &p_end[0], &cut_point[0]) == 1)
   {
     Point p(cut_point[0], cut_point[1], cut_point[2]);
     Real dotp = ((p - p1) * (p2 - p1)) / ((p2 - p1) * (p2 - p1));
     if (isInsideCutPlane(vertices, p) && dotp > 1)
     {
       pint = p;
       has_intersection = true;
     }
   }
   return has_intersection;
 }

◆ getCrackFrontPoints()

const std::vector< Point > CrackMeshCut3DUserObject::getCrackFrontPoints ( unsigned int int ) const

overridevirtual

get a set of points along a crack front from a XFEM GeometricCutUserObject

Returns: A vector which contains all crack front points

Implements CrackFrontPointsProvider.

Definition at line 1217 of file CrackMeshCut3DUserObject.C.

 {
   std::vector<Point> crack_front_points(number_crack_front_points);
   // number_crack_front_points is updated via
   // _crack_front_definition->updateNumberOfCrackFrontPoints(_crack_front_points.size())
   if (number_crack_front_points != _crack_front_points.size())
     mooseError("number_points_from_provider does not match the number of nodes given in "
                "crack_front_nodes");
   for (unsigned int i = 0; i < number_crack_front_points; ++i)
   {
     dof_id_type id = _crack_front_points[i];
     Node * this_node = _cutter_mesh->node_ptr(id);
     mooseAssert(this_node, "Node is NULL");
     Point & this_point = *this_node;
     crack_front_points[i] = this_point;
   }
   return crack_front_points;
 }

◆ getCrackPlaneNormals()

const std::vector< RealVectorValue > CrackMeshCut3DUserObject::getCrackPlaneNormals ( unsigned int int ) const

overridevirtual

get a set of normal vectors along a crack front from a XFEM GeometricCutUserObject

Returns: A vector which contains all crack front normals

Implements CrackFrontPointsProvider.

Definition at line 1237 of file CrackMeshCut3DUserObject.C.

 {
   std::vector<RealVectorValue> crack_plane_normals(number_crack_front_points);
 
   // build the node-to-elems map
   std::unordered_map<dof_id_type, std::vector<dof_id_type>> node_to_elems_map;
   node_to_elems_map.clear();
   for (const auto & elem : _cutter_mesh->element_ptr_range())
     for (auto & node : elem->node_ref_range())
       node_to_elems_map[node.id()].push_back(elem->id());
 
   // build the elem-to-normal map
   std::unordered_map<dof_id_type, RealVectorValue> elem_to_normal_map;
   elem_to_normal_map.clear();
   for (const auto & elem : _cutter_mesh->element_ptr_range())
   {
     Point & p1 = *elem->node_ptr(0);
     Point & p2 = *elem->node_ptr(1);
     Point & p3 = *elem->node_ptr(2);
     Plane elem_plane(p3, p2, p1); // to match the current normal of 0,0,-1;
     RealVectorValue normal = elem_plane.unit_normal(p1);
     elem_to_normal_map[elem->id()] = normal;
   }
 
   // for any front node, the normal is averaged based on the normals of all elements sharing this
   // node this code may fail when the front node has no element connected to it, e.g. refinement at
   // step 1 has to be disabled
   for (unsigned int i = 0; i < number_crack_front_points; ++i)
   {
     dof_id_type id = _crack_front_points[i];
     std::vector<dof_id_type> elems = node_to_elems_map[id];
     unsigned int n_elem = elems.size();
 
     RealVectorValue normal_avr = 0;
     for (unsigned int j = 0; j < n_elem; ++j)
       normal_avr += elem_to_normal_map[elems[j]];
     normal_avr = normal_avr / n_elem;
     crack_plane_normals[i] = normal_avr;
   }
   return crack_plane_normals;
 }

◆ getCutSubdomainID() [1/2]

virtual CutSubdomainID GeometricCutUserObject::getCutSubdomainID ( const Node * ) const

inlinevirtualinherited

Get CutSubdomainID telling which side the node belongs to relative to the cut.

The returned ID contains no physical meaning, but should be consistent throughout the simulation.

Parameters

node	Pointer to the node

Returns: An unsigned int indicating the side

Definition at line 184 of file GeometricCutUserObject.h.

Referenced by CutElementSubdomainModifier::computeSubdomainID(), CutSubdomainIDAux::computeValue(), and XFEM::getCutSubdomainID().

   {
     mooseError("Objects that inherit from GeometricCutUserObject should override the "
                "getCutSubdomainID method");
     return 0;
   }

◆ getCutSubdomainID() [2/2]

CutSubdomainID GeometricCutUserObject::getCutSubdomainID ( const Elem * elem ) const

inherited

Get the CutSubdomainID for the given element.

Parameters

node	Pointer to the element

Returns: The CutSubdomainID

Definition at line 285 of file GeometricCutUserObject.C.

 {
   return _xfem->getCutSubdomainID(this, elem);
 }

◆ getCutterMesh()

MeshBase & MeshCutUserObjectBase::getCutterMesh ( ) const

inherited

Get a reference to the cutter mesh.

Returns: reference to the cutter mesh

Definition at line 58 of file MeshCutUserObjectBase.C.

Referenced by XFEMCutMeshOutput::output().

 {
   mooseAssert(_cutter_mesh, "MeshCutUserObjectBase::getCutterMesh _cutter_mesh is nullptr");
   return *_cutter_mesh;
 }

◆ getFrontPointsIndex()

std::vector< int > CrackMeshCut3DUserObject::getFrontPointsIndex ( )

Get crack front points in the active segment -1 means inactive; positive is the point's index in the Crack Front Definition starting from 0.

Definition at line 1280 of file CrackMeshCut3DUserObject.C.

Referenced by ParisLaw::execute(), and findActiveBoundaryDirection().

 {
   // Crack front definition using the cutter mesh currently only supports one active crack front
   // segment
   unsigned int ibnd = 0;
   unsigned int size_this_segment = _active_boundary[ibnd].size();
   unsigned int n_inactive_nodes = _inactive_boundary_pos.size();
 
   std::vector<int> index(size_this_segment, -1);
 
   unsigned int i1 = n_inactive_nodes == 0 ? 0 : 1;
   unsigned int i2 = n_inactive_nodes == 0 ? size_this_segment : size_this_segment - 1;
 
   // loop over active front points
   for (unsigned int j = i1; j < i2; ++j)
   {
     dof_id_type id = _active_boundary[ibnd][j];
     auto it = std::find(_crack_front_points.begin(), _crack_front_points.end(), id);
     index[j] = std::distance(_crack_front_points.begin(), it);
   }
 
   return index;
 }

◆ getInterfaceID()

unsigned int GeometricCutUserObject::getInterfaceID ( ) const

inlineinherited

Get the interface ID for this cutting object.

Returns: the interface ID

Definition at line 162 of file GeometricCutUserObject.h.

162 { return _interface_id; };

GeometricCutUserObject::_interface_id

unsigned int _interface_id

Associated interface id.

Definition: GeometricCutUserObject.h:203

◆ getNumberOfCrackFrontPoints()

unsigned int CrackMeshCut3DUserObject::getNumberOfCrackFrontPoints ( ) const

Return the total number of crack front points.

This function is currently not called anywhere in the code. Ideally, in a future update, the number of crack front points will be managed by CrackFrontPointsProvider instead of CrackFrontDefinition. In that case, getNumberOfCrackFrontPoints() defined here may be used to override a virtual function defined in CrackFrontPointsProvider

Definition at line 1311 of file CrackMeshCut3DUserObject.C.

 {
   return _num_crack_front_points;
 }

◆ getRelativePosition()

Real CrackMeshCut3DUserObject::getRelativePosition	(	const Point &	p1,
		const Point &	p2,
		const Point &	p
	)		const

protected

Get the relative position of p from p1.

Definition at line 345 of file CrackMeshCut3DUserObject.C.

 {
   Real full_len = (p2 - p1).norm();
   Real len_p1_p = (p - p1).norm();
   return len_p1_p / full_len;
 }

◆ growFront()

void CrackMeshCut3DUserObject::growFront ( )

protected

Grow the cutter mesh.

Definition at line 810 of file CrackMeshCut3DUserObject.C.

Referenced by initialize().

 {
   _front.clear();
 
   for (unsigned int i = 0; i < _active_boundary.size(); ++i)
   {
     std::vector<dof_id_type> temp;
 
     unsigned int i1 = 1;
     unsigned int i2 = _active_boundary[i].size() - 1;
     if (_inactive_boundary_pos.size() == 0)
     {
       i1 = 0;
       i2 = _active_boundary[i].size();
     }
 
     for (unsigned int j = i1; j < i2; ++j)
     {
       Node * this_node = _cutter_mesh->node_ptr(_active_boundary[i][j]);
       mooseAssert(this_node, "Node is NULL");
       Point & this_point = *this_node;
       Point dir = _active_direction[i][j];
 
       Point x;
 
       if (_growth_rate_method == GrowthRateEnum::FUNCTION)
         for (unsigned int k = 0; k < 3; ++k)
         {
           Real velo = _func_v->value(0, Point(0, 0, 0));
           x(k) = this_point(k) + dir(k) * velo;
         }
       else if (_growth_rate_method == GrowthRateEnum::FATIGUE)
       {
         // get the number of loading cycles for this growth increament
         if (j == i1)
         {
           unsigned long int dn = (unsigned long int)_func_v->value(0, Point(0, 0, 0));
           _dn.push_back(dn);
           _n.push_back(_n.size() == 0 ? dn : dn + _n[_n.size() - 1]);
         }
 
         Real growth_size = _growth_size[j];
         for (unsigned int k = 0; k < 3; ++k)
           x(k) = this_point(k) + dir(k) * growth_size;
       }
       else
         mooseError("This growth_rate_method is not pre-defined!");
 
       this_node = Node::build(x, _cutter_mesh->n_nodes()).release();
       _cutter_mesh->add_node(this_node);
 
       dof_id_type id = _cutter_mesh->n_nodes() - 1;
       temp.push_back(id);
 
       if (_cfd)
       {
         auto it = std::find(_tracked_crack_front_points.begin(),
                             _tracked_crack_front_points.end(),
                             _active_boundary[0][j]);
         if (it != _tracked_crack_front_points.end())
         {
           unsigned int pos = std::distance(_tracked_crack_front_points.begin(), it);
           _tracked_crack_front_points[pos] = id;
         }
       }
     }
 
     _front.push_back(temp);
   }
 }

◆ initialize()

void CrackMeshCut3DUserObject::initialize ( )

overridevirtual

Reimplemented from GeometricCutUserObject.

Definition at line 139 of file CrackMeshCut3DUserObject.C.

 {
   _is_mesh_modified = false;
 
   if (_grow)
   {
     if (_t_step == 1)
       _last_step_initialized = 1;
 
     _stop = 0;
 
     if (_t_step > 1 && _t_step != _last_step_initialized)
     {
       _last_step_initialized = _t_step;
 
       for (unsigned int i = 0; i < _n_step_growth; ++i)
       {
         if (_stop != 1)
         {
           findActiveBoundaryNodes();
           findActiveBoundaryDirection();
           _is_mesh_modified = true;
           growFront();
           sortFrontNodes();
           if (_inactive_boundary_pos.size() != 0)
             findFrontIntersection();
           refineFront();
           triangulation();
           joinBoundary();
         }
       }
     }
   }
   if (_cfd)
     _crack_front_definition->isCutterModified(_is_mesh_modified);
 }

◆ initialSetup()

void CrackMeshCut3DUserObject::initialSetup ( )

overridevirtual

Reimplemented from ElementUserObject.

Definition at line 115 of file CrackMeshCut3DUserObject.C.

 {
   if (_cfd)
   {
     _crack_front_definition =
         &_fe_problem.getUserObject<CrackFrontDefinition>("crackFrontDefinition");
     _crack_front_points = _tracked_crack_front_points;
   }
 
   if (_grow)
   {
     findBoundaryNodes();
     findBoundaryEdges();
     sortBoundaryNodes();
   }
 
   if (_growth_rate_method == GrowthRateEnum::FATIGUE)
   {
     _dn.clear();
     _n.clear();
   }
 }

◆ intersectWithEdge()

bool CrackMeshCut3DUserObject::intersectWithEdge	(	const Point &	p1,
		const Point &	p2,
		const std::vector< Point > &	_vertices,
		Point &	point
	)		const

protectedvirtual

Check if a line intersects with an element.

Definition at line 274 of file CrackMeshCut3DUserObject.C.

 {
   bool has_intersection = false;
 
   Plane elem_plane(vertices[0], vertices[1], vertices[2]);
   Point point = vertices[0];
   Point normal = elem_plane.unit_normal(point);
 
   std::array<Real, 3> plane_point = {{point(0), point(1), point(2)}};
   std::array<Real, 3> planenormal = {{normal(0), normal(1), normal(2)}};
   std::array<Real, 3> edge_point1 = {{p1(0), p1(1), p1(2)}};
   std::array<Real, 3> edge_point2 = {{p2(0), p2(1), p2(2)}};
   std::array<Real, 3> cut_point = {{0.0, 0.0, 0.0}};
 
   if (Xfem::plane_normal_line_exp_int_3d(
           &plane_point[0], &planenormal[0], &edge_point1[0], &edge_point2[0], &cut_point[0]) == 1)
   {
     Point temp_p(cut_point[0], cut_point[1], cut_point[2]);
     if (isInsideCutPlane(vertices, temp_p) && isInsideEdge(p1, p2, temp_p))
     {
       pint = temp_p;
       has_intersection = true;
     }
   }
   return has_intersection;
 }

◆ isInsideCutPlane()

bool CrackMeshCut3DUserObject::isInsideCutPlane	(	const std::vector< Point > &	_vertices,
		const Point &	p
	)		const

protected

Check if point p is inside a plane.

Definition at line 355 of file CrackMeshCut3DUserObject.C.

Referenced by findIntersection(), and intersectWithEdge().

 {
   unsigned int n_node = vertices.size();
 
   Plane elem_plane(vertices[0], vertices[1], vertices[2]);
   Point normal = elem_plane.unit_normal(vertices[0]);
 
   bool inside = false;
   unsigned int counter = 0;
 
   for (unsigned int i = 0; i < n_node; ++i)
   {
     unsigned int iplus1 = (i < n_node - 1 ? i + 1 : 0);
     Point middle2p = p - 0.5 * (vertices[i] + vertices[iplus1]);
     const Point side_tang = vertices[iplus1] - vertices[i];
     Point side_norm = side_tang.cross(normal);
     Xfem::normalizePoint(middle2p);
     Xfem::normalizePoint(side_norm);
     if (middle2p * side_norm <= 0.0)
       counter += 1;
   }
   if (counter == n_node)
     inside = true;
   return inside;
 }

◆ isInsideEdge()

bool CrackMeshCut3DUserObject::isInsideEdge	(	const Point &	p1,
		const Point &	p2,
		const Point &	p
	)		const

protected

Check if point p is inside the edge p1-p2.

Definition at line 337 of file CrackMeshCut3DUserObject.C.

Referenced by intersectWithEdge().

 {
   Real dotp1 = (p1 - p) * (p2 - p1);
   Real dotp2 = (p2 - p) * (p2 - p1);
   return (dotp1 * dotp2 <= 0.0);
 }

◆ joinBoundary()

void CrackMeshCut3DUserObject::joinBoundary ( )

protected

Join active boundaries and inactive boundaries to be the new boundary.

Definition at line 1177 of file CrackMeshCut3DUserObject.C.

Referenced by initialize().

 {
   if (_inactive_boundary_pos.size() == 0)
   {
     _boundary = _front[0];
     _boundary.pop_back();
     return;
   }
 
   std::vector<dof_id_type> full_front;
 
   unsigned int size1 = _active_boundary.size();
 
   for (unsigned int i = 0; i < size1; ++i)
   {
     unsigned int size2 = _active_boundary[i].size();
 
     dof_id_type bd1 = _active_boundary[i][size2 - 1];
     dof_id_type bd2 = _active_boundary[i + 1 < size1 ? i + 1 : 0][0];
 
     full_front.insert(full_front.end(), _front[i].begin(), _front[i].end());
 
     auto it1 = std::find(_boundary.begin(), _boundary.end(), bd1);
     unsigned int pos1 = std::distance(_boundary.begin(), it1);
     auto it2 = std::find(_boundary.begin(), _boundary.end(), bd2);
     unsigned int pos2 = std::distance(_boundary.begin(), it2);
 
     if (pos1 <= pos2)
       full_front.insert(full_front.end(), _boundary.begin() + pos1, _boundary.begin() + pos2 + 1);
     else
     {
       full_front.insert(full_front.end(), _boundary.begin() + pos1, _boundary.end());
       full_front.insert(full_front.end(), _boundary.begin(), _boundary.begin() + pos2 + 1);
     }
   }
 
   _boundary = full_front;
 }

◆ refineBoundary()

void CrackMeshCut3DUserObject::refineBoundary ( )

protected

If boundary nodes are too sparse, add nodes in between.

Definition at line 586 of file CrackMeshCut3DUserObject.C.

 {
   std::vector<dof_id_type> new_boundary_order(_boundary.begin(), _boundary.end());
 
   mooseAssert(_boundary.size() >= 2, "Boundary should have at least two nodes");
 
   for (unsigned int i = _boundary.size() - 1; i >= 1; --i)
   {
     dof_id_type node1 = _boundary[i - 1];
     dof_id_type node2 = _boundary[i];
 
     Real distance = findDistance(node1, node2);
 
     if (distance > _size_control)
     {
       unsigned int n = static_cast<unsigned int>(distance / _size_control);
       std::array<Real, 3> x1;
       std::array<Real, 3> x2;
 
       Node * n1 = _cutter_mesh->node_ptr(node1);
       mooseAssert(n1 != nullptr, "Node is NULL");
       Point & p1 = *n1;
       Node * n2 = _cutter_mesh->node_ptr(node2);
       mooseAssert(n2 != nullptr, "Node is NULL");
       Point & p2 = *n2;
 
       for (unsigned int j = 0; j < 3; ++j)
       {
         x1[j] = p1(j);
         x2[j] = p2(j);
       }
 
       for (unsigned int j = 0; j < n; ++j)
       {
         Point x;
         for (unsigned int k = 0; k < 3; ++k)
           x(k) = x2[k] - (x2[k] - x1[k]) * (j + 1) / (n + 1);
 
         Node * this_node = Node::build(x, _cutter_mesh->n_nodes()).release();
         _cutter_mesh->add_node(this_node);
 
         dof_id_type id = _cutter_mesh->n_nodes() - 1;
         auto it = new_boundary_order.begin();
         new_boundary_order.insert(it + i, id);
       }
     }
   }
 
   _boundary = new_boundary_order;
   mooseAssert(_boundary.size() > 0, "Boundary should not have zero size");
   _boundary.pop_back();
 }

◆ refineFront()

void CrackMeshCut3DUserObject::refineFront ( )

protected

Refine the mesh at the front.

Definition at line 1010 of file CrackMeshCut3DUserObject.C.

Referenced by initialize().

 {
   std::vector<std::vector<dof_id_type>> new_front(_front.begin(), _front.end());
 
   for (unsigned int ifront = 0; ifront < _front.size(); ++ifront)
   {
     unsigned int i1 = _front[ifront].size() - 1;
     if (_inactive_boundary_pos.size() == 0)
       i1 = _front[ifront].size();
 
     for (unsigned int i = i1; i >= 1; --i)
     {
       unsigned int i2 = i;
       if (_inactive_boundary_pos.size() == 0)
         i2 = (i <= _front[ifront].size() - 1 ? i : 0);
 
       dof_id_type node1 = _front[ifront][i - 1];
       dof_id_type node2 = _front[ifront][i2];
       Real distance = findDistance(node1, node2);
 
       if (distance > _size_control)
       {
         unsigned int n = static_cast<int>(distance / _size_control);
         std::array<Real, 3> x1;
         std::array<Real, 3> x2;
 
         Node * this_node = _cutter_mesh->node_ptr(node1);
         mooseAssert(this_node, "Node is NULL");
         Point & p1 = *this_node;
         this_node = _cutter_mesh->node_ptr(node2);
         mooseAssert(this_node, "Node is NULL");
         Point & p2 = *this_node;
 
         for (unsigned int j = 0; j < 3; ++j)
         {
           x1[j] = p1(j);
           x2[j] = p2(j);
         }
 
         for (unsigned int j = 0; j < n; ++j)
         {
           Point x;
           for (unsigned int k = 0; k < 3; ++k)
             x(k) = x2[k] - (x2[k] - x1[k]) * (j + 1) / (n + 1);
 
           Node * this_node = Node::build(x, _cutter_mesh->n_nodes()).release();
           _cutter_mesh->add_node(this_node);
 
           dof_id_type id = _cutter_mesh->n_nodes() - 1;
 
           auto it = new_front[ifront].begin();
           new_front[ifront].insert(it + i, id);
         }
       }
     }
   }
 
   _front = new_front;
 
   if (_cfd)
   {
     if (_front[0][0] == _tracked_crack_front_points[0] &&
         _front[0].back() == _tracked_crack_front_points.back())
       _crack_front_points = _front[0];
     else if (_front[0][0] == _tracked_crack_front_points.back() &&
              _front[0].back() == _tracked_crack_front_points[0])
     {
       _crack_front_points = _front[0];
       std::reverse(_crack_front_points.begin(), _crack_front_points.end());
     }
     else
       mooseError("the crack front and the tracked crack front definition must match in terms of "
                  "their end nodes");
 
     _num_crack_front_points = _crack_front_points.size();
     _crack_front_definition->updateNumberOfCrackFrontPoints(_num_crack_front_points);
   }
 }

◆ serialize()

void GeometricCutUserObject::serialize ( std::string & serialized_buffer )

protectedinherited

Methods to pack/unpack the _marked_elems_2d and _marked_elems_3d data into a structure suitable for parallel communication.

Definition at line 206 of file GeometricCutUserObject.C.

Referenced by GeometricCutUserObject::finalize().

 {
   // stream for serializing the _marked_elems_2d and _marked_elems_3d data structures to a byte
   // stream
   std::ostringstream oss;
   dataStore(oss, _marked_elems_2d, this);
   dataStore(oss, _marked_elems_3d, this);
 
   // Populate the passed in string pointer with the string stream's buffer contents
   serialized_buffer.assign(oss.str());
 }

◆ setInterfaceID()

void GeometricCutUserObject::setInterfaceID ( unsigned int interface_id )

inlineinherited

Set the interface ID for this cutting object.

Parameters

the	interface ID

Definition at line 168 of file GeometricCutUserObject.h.

Referenced by XFEM::addGeometricCut().

168 { _interface_id = interface_id; };

GeometricCutUserObject::_interface_id

unsigned int _interface_id

Associated interface id.

Definition: GeometricCutUserObject.h:203

◆ setSubCriticalGrowthSize()

void CrackMeshCut3DUserObject::setSubCriticalGrowthSize ( std::vector< Real > & growth_size )

Return growth size at the active boundary to the mesh cutter.

Definition at line 1305 of file CrackMeshCut3DUserObject.C.

Referenced by ParisLaw::execute().

 {
   _growth_size = growth_size;
 }

◆ shouldHealMesh()

bool GeometricCutUserObject::shouldHealMesh ( ) const

inlineinherited

Should the elements cut by this cutting object be healed in the current time step?

Returns: true if the cut element should be healed

Definition at line 175 of file GeometricCutUserObject.h.

Referenced by XFEM::cutMeshWithEFA().

175 { return _heal_always; };

GeometricCutUserObject::_heal_always

bool _heal_always

Heal the mesh.

Definition: GeometricCutUserObject.h:206

◆ sortBoundaryNodes()

void CrackMeshCut3DUserObject::sortBoundaryNodes ( )

protected

Sort boundary nodes to be in the right order along the boundary.

Definition at line 515 of file CrackMeshCut3DUserObject.C.

Referenced by initialSetup().

 {
   _boundary.clear();
 
   for (auto it = _boundary_edges.begin(); it != _boundary_edges.end(); ++it)
   {
     dof_id_type node1 = (*it)._id1;
     dof_id_type node2 = (*it)._id2;
 
     mooseAssert(_boundary_map.find(node1) != _boundary_map.end(),
                 "_boundary_map does not have this key");
     mooseAssert(_boundary_map.find(node2) != _boundary_map.end(),
                 "_boundary_map does not have this key");
 
     _boundary_map.find(node1)->second.push_back(node2);
     _boundary_map.find(node2)->second.push_back(node1);
   }
 
   auto it = _boundary_map.begin();
   while (it != _boundary_map.end())
   {
     if (it->second.size() != 2)
       mooseError(
           "Boundary nodes in the cutter mesh must have exactly two neighbors; this one has: ",
           it->second.size());
     ++it;
   }
 
   auto it2 = _boundary_edges.begin();
   dof_id_type node1 = (*it2)._id1;
   dof_id_type node2 = (*it2)._id2;
   _boundary.push_back(node1);
   _boundary.push_back(node2);
 
   for (unsigned int i = 0; i < _boundary_edges.size() - 1; ++i)
   {
     mooseAssert(_boundary_map.find(node2) != _boundary_map.end(),
                 "_boundary_map does not have this key");
 
     dof_id_type node3 = _boundary_map.find(node2)->second[0];
     dof_id_type node4 = _boundary_map.find(node2)->second[1];
 
     if (node3 == node1)
     {
       _boundary.push_back(node4);
       node1 = node2;
       node2 = node4;
     }
     else if (node4 == node1)
     {
       _boundary.push_back(node3);
       node1 = node2;
       node2 = node3;
     }
     else
       mooseError("Discontinuity in cutter boundary");
   }
 }

◆ sortFrontNodes()

void CrackMeshCut3DUserObject::sortFrontNodes ( )

protected

Sort the front nodes.

Definition at line 882 of file CrackMeshCut3DUserObject.C.

Referenced by initialize().

884 {

885 }

◆ threadJoin()

void GeometricCutUserObject::threadJoin ( const UserObject & y )

overridevirtualinherited

Implements ElementUserObject.

Definition at line 130 of file GeometricCutUserObject.C.

 {
   const auto & gcuo = static_cast<const GeometricCutUserObject &>(y);
 
   for (const auto & it : gcuo._marked_elems_2d)
   {
     mooseAssert(_marked_elems_2d.find(it.first) == _marked_elems_2d.end(),
                 "Element already inserted in map from a different thread");
     _marked_elems_2d[it.first] = it.second;
   }
   for (const auto & it : gcuo._marked_elems_3d)
   {
     mooseAssert(_marked_elems_3d.find(it.first) == _marked_elems_3d.end(),
                 "Element already inserted in map from a different thread");
     _marked_elems_3d[it.first] = it.second;
   }
 }

◆ triangulation()

void CrackMeshCut3DUserObject::triangulation ( )

protected

Create tri3 elements between the new front and the old front.

Definition at line 1090 of file CrackMeshCut3DUserObject.C.

Referenced by initialize().

 {
 
   mooseAssert(_active_boundary.size() == _front.size(),
               "_active_boundary and _front should have the same size!");
 
   if (_inactive_boundary_pos.size() == 0)
   {
     _active_boundary[0].push_back(_active_boundary[0][0]);
     _front[0].push_back(_front[0][0]);
   }
 
   // loop over active segments
   for (unsigned int k = 0; k < _front.size(); ++k)
   {
     unsigned int n1 = _active_boundary[k].size();
     unsigned int n2 = _front[k].size();
 
     unsigned int i1 = 0;
     unsigned int i2 = 0;
 
     // stop when all nodes are associated with an element
     while (!(i1 == n1 - 1 && i2 == n2 - 1))
     {
       std::vector<dof_id_type> elem;
 
       dof_id_type p1 = _active_boundary[k][i1]; // node in the old front
       dof_id_type p2 = _front[k][i2];           // node in the new front
 
       if (i1 != n1 - 1 && i2 != n2 - 1)
       {
         dof_id_type p3 = _active_boundary[k][i1 + 1]; // next node in the old front
         dof_id_type p4 = _front[k][i2 + 1];           // next node in the new front
 
         elem.push_back(p1);
         elem.push_back(p2);
 
         Real d1 = findDistance(p1, p4);
         Real d2 = findDistance(p3, p2);
 
         if (d1 < d2)
         {
           elem.push_back(p4);
           i2++;
         }
 
         else
         {
           elem.push_back(p3);
           i1++;
         }
       }
 
       else if (i1 == n1 - 1)
       {
         dof_id_type p4 = _front[k][i2 + 1]; // next node in the new front
 
         elem.push_back(p1);
         elem.push_back(p2);
         elem.push_back(p4);
         i2++;
       }
 
       else if (i2 == n2 - 1)
       {
         dof_id_type p3 = _active_boundary[k][i1 + 1]; // next node in the old front
 
         elem.push_back(p1);
         elem.push_back(p2);
         elem.push_back(p3);
         i1++;
       }
 
       Elem * new_elem = Elem::build(TRI3).release();
 
       for (unsigned int i = 0; i < _cut_elem_nnode; ++i)
       {
         mooseAssert(_cutter_mesh->node_ptr(elem[i]) != nullptr, "Node is NULL");
         new_elem->set_node(i, _cutter_mesh->node_ptr(elem[i]));
       }
 
       _cutter_mesh->add_elem(new_elem);
     }
   }
 }

◆ usesMesh()

bool CrackFrontPointsProvider::usesMesh ( ) const

inlineinherited

Getter for if a cutter mesh is used in a derived class.

Returns: bool indicating if a cutter mesh is used in the derived class

Definition at line 40 of file CrackFrontPointsProvider.h.

Referenced by CrackFrontDefinition::initialSetup().

40 { return _uses_mesh; }

CrackFrontPointsProvider::_uses_mesh

const bool _uses_mesh

bool to set if CrackFrontPointsProvider derived objects use a cutter mesh

Definition: CrackFrontPointsProvider.h:44

◆ validParams()

InputParameters CrackMeshCut3DUserObject::validParams ( )

static

Definition at line 27 of file CrackMeshCut3DUserObject.C.

 {
   InputParameters params = MeshCutUserObjectBase::validParams();
   MooseEnum growthDirection("MAX_HOOP_STRESS FUNCTION", "FUNCTION");
   params.addParam<MooseEnum>(
       "growth_dir_method", growthDirection, "choose from FUNCTION, MAX_HOOP_STRESS");
   MooseEnum growthRate("FATIGUE FUNCTION", "FUNCTION");
   params.addParam<MooseEnum>("growth_rate_method", growthRate, "choose from FUNCTION, FATIGUE");
   params.addParam<FunctionName>("growth_direction_x",
                                 "Function defining x-component of crack growth direction");
   params.addParam<FunctionName>("growth_direction_y",
                                 "Function defining y-component of crack growth direction");
   params.addParam<FunctionName>("growth_direction_z",
                                 "Function defining z-component of crack growth direction");
 
   params.addParam<FunctionName>("growth_rate", "Function defining crack growth rate");
   params.addParam<Real>(
       "size_control", 0, "Criterion for refining elements while growing the crack");
   params.addParam<unsigned int>("n_step_growth", 0, "Number of steps for crack growth");
   params.addParam<std::vector<dof_id_type>>("crack_front_nodes",
                                             "Set of nodes to define crack front");
   params.addClassDescription("Creates a UserObject for a mesh cutter in 3D problems");
   return params;
 }

Member Data Documentation

◆ _active_boundary

std::vector<std::vector<dof_id_type> > CrackMeshCut3DUserObject::_active_boundary

protected

Active boundary nodes where growth is allowed.

Definition at line 133 of file CrackMeshCut3DUserObject.h.

Referenced by findActiveBoundaryDirection(), findActiveBoundaryNodes(), getFrontPointsIndex(), growFront(), joinBoundary(), and triangulation().

◆ _active_direction

std::vector<std::vector<Point> > CrackMeshCut3DUserObject::_active_direction

protected

Growth direction for active boundaries.

Definition at line 151 of file CrackMeshCut3DUserObject.h.

Referenced by findActiveBoundaryDirection(), and growFront().

◆ _boundary

std::vector<dof_id_type> CrackMeshCut3DUserObject::_boundary

protected

Boundary nodes of the cutter mesh.

Definition at line 130 of file CrackMeshCut3DUserObject.h.

Referenced by findActiveBoundaryNodes(), joinBoundary(), refineBoundary(), and sortBoundaryNodes().

◆ _boundary_edges

std::set<Xfem::CutEdge> CrackMeshCut3DUserObject::_boundary_edges

protected

Edges at the boundary.

Definition at line 145 of file CrackMeshCut3DUserObject.h.

Referenced by findBoundaryEdges(), and sortBoundaryNodes().

◆ _boundary_map

std::map<dof_id_type, std::vector<dof_id_type> > CrackMeshCut3DUserObject::_boundary_map

protected

A map of boundary nodes and their neighbors.

Definition at line 148 of file CrackMeshCut3DUserObject.h.

Referenced by findBoundaryEdges(), findBoundaryNodes(), and sortBoundaryNodes().

◆ _cfd

bool CrackMeshCut3DUserObject::_cfd

protected

Definition at line 142 of file CrackMeshCut3DUserObject.h.

Referenced by CrackMeshCut3DUserObject(), findActiveBoundaryDirection(), findFrontIntersection(), growFront(), initialize(), initialSetup(), and refineFront().

◆ _const_intersection

const Real CrackMeshCut3DUserObject::_const_intersection = 0.01

protected

Used to define intersection points.

Definition at line 117 of file CrackMeshCut3DUserObject.h.

Referenced by findFrontIntersection().

◆ _crack_front_definition

CrackFrontDefinition* CrackMeshCut3DUserObject::_crack_front_definition

protected

The crack front definition.

Definition at line 86 of file CrackMeshCut3DUserObject.h.

Referenced by findActiveBoundaryDirection(), initialize(), initialSetup(), and refineFront().

◆ _crack_front_points

std::vector<dof_id_type> CrackMeshCut3DUserObject::_crack_front_points

protected

updated crack front definition they are in the same order as defined in the input but the number of nodes may increase its difference from _front is that: _front does not necessarily follow the order of crack front definition therefore, _crack_front_points is generated from _front with the order of crack front definition limitation: this approach does not currently support the growth of one crack front into two

Definition at line 93 of file CrackMeshCut3DUserObject.h.

Referenced by findActiveBoundaryDirection(), getCrackFrontPoints(), getCrackPlaneNormals(), getFrontPointsIndex(), initialSetup(), and refineFront().

◆ _cut_elem_dim

const unsigned int CrackMeshCut3DUserObject::_cut_elem_dim = 2

protected

Definition at line 80 of file CrackMeshCut3DUserObject.h.

Referenced by CrackMeshCut3DUserObject().

◆ _cut_elem_nnode

const unsigned int CrackMeshCut3DUserObject::_cut_elem_nnode = 3

protected

The cutter mesh has triangluar elements only.

Definition at line 79 of file CrackMeshCut3DUserObject.h.

Referenced by CrackMeshCut3DUserObject(), findBoundaryEdges(), and triangulation().

◆ _cutter_mesh

std::unique_ptr<MeshBase> MeshCutUserObjectBase::_cutter_mesh

protectedinherited

The xfem cutter mesh.

Definition at line 32 of file MeshCutUserObjectBase.h.

◆ _dn

std::vector<unsigned long int> CrackMeshCut3DUserObject::_dn

protected

Fatigue life.

Definition at line 157 of file CrackMeshCut3DUserObject.h.

Referenced by growFront(), and initialSetup().

◆ _elem_dim

const unsigned int CrackMeshCut3DUserObject::_elem_dim = 3

protected

The structural mesh must be 3D only.

Definition at line 114 of file CrackMeshCut3DUserObject.h.

◆ _front

std::vector<std::vector<dof_id_type> > CrackMeshCut3DUserObject::_front

protected

New boundary after growth.

Definition at line 161 of file CrackMeshCut3DUserObject.h.

Referenced by findFrontIntersection(), growFront(), joinBoundary(), refineFront(), and triangulation().

◆ _func_v

const Function* CrackMeshCut3DUserObject::_func_v

protected

Definition at line 269 of file CrackMeshCut3DUserObject.h.

Referenced by CrackMeshCut3DUserObject(), and growFront().

◆ _func_x

const Function* CrackMeshCut3DUserObject::_func_x

protected

Parsed functions of front growth.

Definition at line 266 of file CrackMeshCut3DUserObject.h.

Referenced by CrackMeshCut3DUserObject(), and findActiveBoundaryDirection().

◆ _func_y

const Function* CrackMeshCut3DUserObject::_func_y

protected

Definition at line 267 of file CrackMeshCut3DUserObject.h.

Referenced by CrackMeshCut3DUserObject(), and findActiveBoundaryDirection().

◆ _func_z

const Function* CrackMeshCut3DUserObject::_func_z

protected

Definition at line 268 of file CrackMeshCut3DUserObject.h.

Referenced by CrackMeshCut3DUserObject(), and findActiveBoundaryDirection().

◆ _grow

bool CrackMeshCut3DUserObject::_grow

protected

Definition at line 127 of file CrackMeshCut3DUserObject.h.

Referenced by CrackMeshCut3DUserObject(), initialize(), and initialSetup().

◆ _growth_dir_method

const GrowthDirectionEnum CrackMeshCut3DUserObject::_growth_dir_method

protected

The direction method for growing mesh at the front.

Definition at line 102 of file CrackMeshCut3DUserObject.h.

Referenced by CrackMeshCut3DUserObject(), and findActiveBoundaryDirection().

◆ _growth_rate_method

const GrowthRateEnum CrackMeshCut3DUserObject::_growth_rate_method

protected

The rate method for growing mesh at the front.

Definition at line 111 of file CrackMeshCut3DUserObject.h.

Referenced by CrackMeshCut3DUserObject(), growFront(), and initialSetup().

◆ _growth_size

std::vector<Real> CrackMeshCut3DUserObject::_growth_size

protected

Growth size for the active boundary in a subcritical simulation.

Definition at line 154 of file CrackMeshCut3DUserObject.h.

Referenced by growFront(), and setSubCriticalGrowthSize().

◆ _heal_always

bool GeometricCutUserObject::_heal_always

protectedinherited

Heal the mesh.

Definition at line 206 of file GeometricCutUserObject.h.

◆ _inactive_boundary_pos

std::vector<unsigned int> CrackMeshCut3DUserObject::_inactive_boundary_pos

protected

Inactive boundary.

Definition at line 136 of file CrackMeshCut3DUserObject.h.

Referenced by findActiveBoundaryDirection(), findActiveBoundaryNodes(), getFrontPointsIndex(), growFront(), initialize(), joinBoundary(), refineFront(), and triangulation().

◆ _interface_id

unsigned int GeometricCutUserObject::_interface_id

protectedinherited

Associated interface id.

Definition at line 203 of file GeometricCutUserObject.h.

Referenced by GeometricCutUserObject::finalize(), and GeometricCutUserObject::GeometricCutUserObject().

◆ _is_mesh_modified

bool CrackMeshCut3DUserObject::_is_mesh_modified

protected

Indicator that shows if the cutting mesh is modified or not in this calculation step.

Definition at line 164 of file CrackMeshCut3DUserObject.h.

Referenced by initialize().

◆ _last_step_initialized

int GeometricCutUserObject::_last_step_initialized

protectedinherited

Time step information needed to advance a 3D crack only at the real beginning of a time step.

Definition at line 209 of file GeometricCutUserObject.h.

Referenced by initialize().

◆ _marked_elems_2d

std::map<unsigned int, std::vector<Xfem::GeomMarkedElemInfo2D> > GeometricCutUserObject::_marked_elems_2d

protectedinherited

Containers with information about all 2D and 3D elements marked for cutting by this object.

Definition at line 212 of file GeometricCutUserObject.h.

Referenced by GeometricCutUserObject::deserialize(), GeometricCutUserObject::execute(), GeometricCutUserObject::finalize(), GeometricCutUserObject::initialize(), GeometricCutUserObject::serialize(), and GeometricCutUserObject::threadJoin().

◆ _marked_elems_3d

std::map<unsigned int, std::vector<Xfem::GeomMarkedElemInfo3D> > GeometricCutUserObject::_marked_elems_3d

protectedinherited

Definition at line 213 of file GeometricCutUserObject.h.

Referenced by GeometricCutUserObject::deserialize(), GeometricCutUserObject::execute(), GeometricCutUserObject::finalize(), GeometricCutUserObject::initialize(), GeometricCutUserObject::serialize(), and GeometricCutUserObject::threadJoin().

◆ _mesh

MooseMesh& CrackMeshCut3DUserObject::_mesh

protected

The structural mesh.

Definition at line 83 of file CrackMeshCut3DUserObject.h.

Referenced by findActiveBoundaryNodes(), and findFrontIntersection().

◆ _n

std::vector<unsigned long int> CrackMeshCut3DUserObject::_n

protected

Definition at line 158 of file CrackMeshCut3DUserObject.h.

Referenced by growFront(), and initialSetup().

◆ _n_step_growth

unsigned int CrackMeshCut3DUserObject::_n_step_growth

protected

Number of steps to grow the mesh.

Definition at line 123 of file CrackMeshCut3DUserObject.h.

Referenced by CrackMeshCut3DUserObject(), and initialize().

◆ _num_crack_front_points

unsigned int CrackMeshCut3DUserObject::_num_crack_front_points

protected

Total number of crack front points in the mesh cutter.

Definition at line 167 of file CrackMeshCut3DUserObject.h.

Referenced by CrackMeshCut3DUserObject(), getNumberOfCrackFrontPoints(), and refineFront().

◆ _size_control

Real CrackMeshCut3DUserObject::_size_control

protected

Used for cutter mesh refinement and front advancement.

Definition at line 120 of file CrackMeshCut3DUserObject.h.

Referenced by CrackMeshCut3DUserObject(), refineBoundary(), and refineFront().

◆ _stop

bool CrackMeshCut3DUserObject::_stop

protected

Variables to help control the work flow.

Definition at line 126 of file CrackMeshCut3DUserObject.h.

Referenced by findActiveBoundaryNodes(), and initialize().

◆ _tracked_crack_front_points

std::vector<dof_id_type> CrackMeshCut3DUserObject::_tracked_crack_front_points

protected

Front nodes that are grown from the crack front definition defined in the input therefore, they are (1) in the same order as defined in the input and (2) the number of nodes does not change.

Definition at line 140 of file CrackMeshCut3DUserObject.h.

Referenced by CrackMeshCut3DUserObject(), findFrontIntersection(), growFront(), initialSetup(), and refineFront().

◆ _uses_mesh

const bool CrackFrontPointsProvider::_uses_mesh

protectedinherited

bool to set if CrackFrontPointsProvider derived objects use a cutter mesh

Definition at line 44 of file CrackFrontPointsProvider.h.

Referenced by CrackFrontPointsProvider::usesMesh().

◆ _xfem

std::shared_ptr<XFEM> GeometricCutUserObject::_xfem

protectedinherited

Pointer to the XFEM controller object.

Definition at line 200 of file GeometricCutUserObject.h.

Referenced by GeometricCutUserObject::execute(), GeometricCutUserObject::finalize(), GeometricCutUserObject::GeometricCutUserObject(), and GeometricCutUserObject::getCutSubdomainID().

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

xfem/include/userobjects/CrackMeshCut3DUserObject.h
xfem/src/userobjects/CrackMeshCut3DUserObject.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

Member Enumeration Documentation

◆ GrowthDirectionEnum

◆ GrowthRateEnum

Constructor & Destructor Documentation

◆ CrackMeshCut3DUserObject()

Member Function Documentation

◆ cutElementByGeometry() [1/4]

◆ cutElementByGeometry() [2/4]

◆ cutElementByGeometry() [3/4]

◆ cutElementByGeometry() [4/4]

◆ cutFragmentByGeometry() [1/4]

◆ cutFragmentByGeometry() [2/4]

◆ cutFragmentByGeometry() [3/4]

◆ cutFragmentByGeometry() [4/4]

◆ deserialize()

◆ execute()

◆ finalize()

◆ findActiveBoundaryDirection()

◆ findActiveBoundaryNodes()

◆ findBoundaryEdges()

◆ findBoundaryNodes()

◆ findDistance()

◆ findFrontIntersection()

◆ findIntersection()

◆ getCrackFrontPoints()

◆ getCrackPlaneNormals()

◆ getCutSubdomainID() [1/2]

◆ getCutSubdomainID() [2/2]

◆ getCutterMesh()

◆ getFrontPointsIndex()

◆ getInterfaceID()

◆ getNumberOfCrackFrontPoints()

◆ getRelativePosition()

◆ growFront()

◆ initialize()

◆ initialSetup()

◆ intersectWithEdge()

◆ isInsideCutPlane()

◆ isInsideEdge()

◆ joinBoundary()

◆ refineBoundary()

◆ refineFront()

◆ serialize()

◆ setInterfaceID()

◆ setSubCriticalGrowthSize()

◆ shouldHealMesh()

◆ sortBoundaryNodes()

◆ sortFrontNodes()

◆ threadJoin()

◆ triangulation()

◆ usesMesh()

◆ validParams()

Member Data Documentation

◆ _active_boundary

◆ _active_direction

◆ _boundary

◆ _boundary_edges

◆ _boundary_map

◆ _cfd

◆ _const_intersection

◆ _crack_front_definition

◆ _crack_front_points

◆ _cut_elem_dim

◆ _cut_elem_nnode

◆ _cutter_mesh

◆ _dn

◆ _elem_dim

◆ _front

◆ _func_v