Provides Xfem interface for a structural element. More...

#include <xfemstructuralelementinterface.h>

Inheritance diagram for oofem::XfemStructuralElementInterface:

Collaboration diagram for oofem::XfemStructuralElementInterface:

Public Member Functions
	XfemStructuralElementInterface (Element *e)

virtual	~XfemStructuralElementInterface ()

virtual bool	XfemElementInterface_updateIntegrationRule ()
	Updates integration rule based on the triangulation. More...

MaterialStatus *	giveClosestGP_MatStat (double &oClosestDist, std::vector< std::unique_ptr< IntegrationRule > > &iRules, const FloatArray &iCoord)

double	computeEffectiveSveSize (StructuralFE2MaterialStatus *iFe2Ms)

virtual void	XfemElementInterface_computeConstitutiveMatrixAt (FloatMatrix &answer, MatResponseMode rMode, GaussPoint , TimeStep tStep)

virtual void	XfemElementInterface_computeStressVector (FloatArray &answer, const FloatArray &strain, GaussPoint gp, TimeStep tStep)

virtual bool	hasCohesiveZone () const

virtual void	computeCohesiveForces (FloatArray &answer, TimeStep *tStep)

virtual void	computeGlobalCohesiveTractionVector (FloatArray &oT, const FloatArray &iJump, const FloatArray &iCrackNormal, const FloatMatrix &iNMatrix, GaussPoint &iGP, TimeStep *tStep)

virtual void	computeCohesiveTangent (FloatMatrix &answer, TimeStep *tStep)

virtual void	computeCohesiveTangentAt (FloatMatrix &answer, TimeStep *tStep)

virtual void	XfemElementInterface_computeConsistentMassMatrix (FloatMatrix &answer, TimeStep tStep, double &mass, const double ipDensity=NULL)

virtual IRResultType	initializeCZFrom (InputRecord *ir)

virtual void	giveCZInputRecord (DynamicInputRecord &input)

virtual void	initializeCZMaterial ()

bool	useNonStdCz ()

virtual void	XfemElementInterface_computeDeformationGradientVector (FloatArray &answer, GaussPoint gp, TimeStep tStep)

void	giveIntersectionsTouchingCrack (std::vector< int > &oTouchingEnrItemIndices, const std::vector< int > &iCandidateIndices, int iEnrItemIndex, XfemManager &iXMan)
	Identify enrichment items with an intersection enrichment front that touches a given enrichment item. More...

void	giveSubtriangulationCompositeExportData (std::vector< VTKPiece > &vtkPieces, IntArray &primaryVarsToExport, IntArray &internalVarsToExport, IntArray cellVarsToExport, TimeStep *tStep)
	VTK Interface. More...

void	computeIPAverageInTriangle (FloatArray &answer, IntegrationRule iRule, Element elem, InternalStateType isType, TimeStep *tStep, const Triangle &iTri)
	Help functions for VTK export. More...

Public Member Functions inherited from oofem::XfemElementInterface
virtual const char *	giveClassName () const

std::string	errorInfo (const char *func) const

	XfemElementInterface (Element *e)
	Constructor. More...

virtual	~XfemElementInterface ()

	XfemElementInterface (const XfemElementInterface &src)=delete

XfemElementInterface &	operator= (const XfemElementInterface &src)=delete

void	XfemElementInterface_createEnrBmatrixAt (FloatMatrix &oAnswer, GaussPoint &iGP, Element &iEl)
	Creates enriched B-matrix. More...

void	XfemElementInterface_createEnrBHmatrixAt (FloatMatrix &oAnswer, GaussPoint &iGP, Element &iEl)
	Creates enriched BH-matrix. More...

void	ComputeBOrBHMatrix (FloatMatrix &oAnswer, GaussPoint &iGP, Element &iEl, bool iComputeBH, const FloatArray &iNaturalGpCoord)
	Help function for computation of B and BH. More...

void	XfemElementInterface_createEnrNmatrixAt (FloatMatrix &oAnswer, const FloatArray &iLocCoord, Element &iEl, bool iSetDiscontContribToZero)
	Creates enriched N-matrix. More...

void	XfemElementInterface_createEnrNmatrixAt (FloatMatrix &oAnswer, const FloatArray &iLocCoord, Element &iEl, const std::vector< int > &iLocNodeInd, bool iSetDiscontContribToZero)
	Creates enriched N-matrix for a chosen subset of element nodes. More...

int	XfemElementInterface_giveNumDofManEnrichments (const DofManager &iDMan, XfemManager &iXMan) const
	Computes total number of enrichments in a node. More...

virtual void	XfemElementInterface_partitionElement (std::vector< Triangle > &oTriangles, const std::vector< FloatArray > &iPoints)
	Partitions the element into patches by a triangulation. More...

virtual void	XfemElementInterface_prepareNodesForDelaunay (std::vector< std::vector< FloatArray > > &oPointPartitions, double &oCrackStartXi, double &oCrackEndXi, int iEnrItemIndex, bool &oIntersection)
	Returns an array of array of points. Each array of points defines the points of a subregion of the element. More...

virtual void	XfemElementInterface_prepareNodesForDelaunay (std::vector< std::vector< FloatArray > > &oPointPartitions, double &oCrackStartXi, double &oCrackEndXi, const Triangle &iTri, int iEnrItemIndex, bool &oIntersection)

void	putPointsInCorrectPartition (std::vector< std::vector< FloatArray > > &oPointPartitions, const std::vector< FloatArray > &iIntersecPoints, const std::vector< const FloatArray * > &iNodeCoord) const

void	partitionEdgeSegment (int iBndIndex, std::vector< Line > &oSegments, std::vector< FloatArray > &oIntersectionPoints, const double &iTangDistPadding=0.0)
	Partition a boundary segment to account for cracks cutting the boundary. More...

MaterialMode	giveMaterialMode ()

void	updateYourselfCZ (TimeStep *tStep)

void	computeDisplacementJump (GaussPoint &iGP, FloatArray &oJump, const FloatArray &iSolVec, const FloatMatrix &iNMatrix)

void	computeNCohesive (FloatMatrix &oN, GaussPoint &iGP, int iEnrItemIndex, const std::vector< int > &iTouchingEnrItemIndices)
	Compute N-matrix for cohesive zone. More...

Public Member Functions inherited from oofem::Interface
	Interface ()
	Constructor. More...

virtual	~Interface ()

Public Attributes
Material *	mpCZMat

int	mCZMaterialNum

int	mCSNumGaussPoints

bool	mIncludeBulkJump

bool	mIncludeBulkCorr

std::vector< Triangle >	mSubTri

Public Attributes inherited from oofem::XfemElementInterface
Element *	element

std::vector< int >	mCZEnrItemIndices
	Index of enrichment items associated with cohesive zones. More...

std::vector< std::vector< int > >	mCZTouchingEnrItemIndices
	Indices of enrichment items that give cohesive zone contributions to a given GP, even though the GP is not located on any of these enrichment items. More...

bool	mUsePlaneStrain
	Flag that tells if plane stress or plane strain is assumed. More...

std::vector< std::unique_ptr< IntegrationRule > >	mpCZIntegrationRules

std::vector< std::unique_ptr< IntegrationRule > >	mpCZExtraIntegrationRules

std::vector< std::unique_ptr< IntegrationRule > >	mpCZIntegrationRules_tmp

std::vector< std::unique_ptr< IntegrationRule > >	mpCZExtraIntegrationRules_tmp

std::vector< std::unique_ptr< IntegrationRule > >	mIntRule_tmp

Detailed Description

Provides Xfem interface for a structural element.

Author: Erik Svenning

Date: Feb 14, 2014

Definition at line 50 of file xfemstructuralelementinterface.h.

Constructor & Destructor Documentation

oofem::XfemStructuralElementInterface::XfemStructuralElementInterface ( Element * e )

Definition at line 80 of file xfemstructuralelementinterface.C.

oofem::XfemStructuralElementInterface::~XfemStructuralElementInterface ( )

virtual

Definition at line 92 of file xfemstructuralelementinterface.C.

Member Function Documentation

void oofem::XfemStructuralElementInterface::computeCohesiveForces	(	FloatArray &	answer,
		TimeStep *	tStep
	)

virtual

Definition at line 876 of file xfemstructuralelementinterface.C.

Referenced by oofem::PlaneStress2dXfem::giveInternalForcesVector(), oofem::TrPlaneStress2dXFEM::giveInternalForcesVector(), and oofem::QTrPlaneStress2dXFEM::giveInternalForcesVector().

void oofem::XfemStructuralElementInterface::computeCohesiveTangent	(	FloatMatrix &	answer,
		TimeStep *	tStep
	)

virtual

Definition at line 1136 of file xfemstructuralelementinterface.C.

Referenced by oofem::PlaneStress2dXfem::computeStiffnessMatrix(), oofem::TrPlaneStress2dXFEM::computeStiffnessMatrix(), and oofem::QTrPlaneStress2dXFEM::computeStiffnessMatrix().

void oofem::XfemStructuralElementInterface::computeCohesiveTangentAt	(	FloatMatrix &	answer,
		TimeStep *	tStep
	)

virtual

Definition at line 1493 of file xfemstructuralelementinterface.C.

References hasCohesiveZone().

double oofem::XfemStructuralElementInterface::computeEffectiveSveSize ( StructuralFE2MaterialStatus * iFe2Ms )

Definition at line 765 of file xfemstructuralelementinterface.C.

References oofem::PrescribedGradientHomogenization::domainSize(), oofem::StructuralFE2MaterialStatus::giveBC(), oofem::StructuralFE2MaterialStatus::giveNormal(), and M_PI.

Referenced by computeCohesiveForces(), and computeCohesiveTangent().

void oofem::XfemStructuralElementInterface::computeGlobalCohesiveTractionVector	(	FloatArray &	oT,
		const FloatArray &	iJump,
		const FloatArray &	iCrackNormal,
		const FloatMatrix &	iNMatrix,
		GaussPoint &	iGP,
		TimeStep *	tStep
	)

virtual

Definition at line 1081 of file xfemstructuralelementinterface.C.

References oofem::FloatArray::at(), oofem::FloatArray::beProductOf(), oofem::FloatArray::beTProductOf(), oofem::FloatMatrix::beUnitMatrix(), oofem::FloatArray::beVectorProductOf(), oofem::StructuralInterfaceMaterial::giveFirstPKTraction_3d(), mpCZMat, OOFEM_ERROR, oofem::FloatMatrix::resize(), and oofem::FloatMatrix::setColumn().

Referenced by computeCohesiveForces(), and computeCohesiveTangent().

void oofem::XfemStructuralElementInterface::computeIPAverageInTriangle	(	FloatArray &	answer,
		IntegrationRule *	iRule,
		Element *	elem,
		InternalStateType	isType,
		TimeStep *	tStep,
		const Triangle &	iTri
	)

Help functions for VTK export.

Definition at line 2027 of file xfemstructuralelementinterface.C.

References oofem::FloatArray::add(), oofem::FloatArray::clear(), oofem::Element::giveIPValue(), oofem::Triangle::pointIsInTriangle(), and oofem::FloatArray::times().

Referenced by giveSubtriangulationCompositeExportData().

MaterialStatus * oofem::XfemStructuralElementInterface::giveClosestGP_MatStat	(	double &	oClosestDist,
		std::vector< std::unique_ptr< IntegrationRule > > &	iRules,
		const FloatArray &	iCoord
	)

Definition at line 742 of file xfemstructuralelementinterface.C.

References oofem::FloatArray::distance_square(), and max.

Referenced by XfemElementInterface_updateIntegrationRule().

void oofem::XfemStructuralElementInterface::giveCZInputRecord ( DynamicInputRecord & input )

virtual

Definition at line 1590 of file xfemstructuralelementinterface.C.

References _IFT_XfemElementInterface_CohesiveZoneMaterial, _IFT_XfemElementInterface_NumIntPointsCZ, _IFT_XfemElementInterface_PlaneStrain, mCSNumGaussPoints, mCZMaterialNum, oofem::XfemElementInterface::mUsePlaneStrain, and oofem::DynamicInputRecord::setField().

Referenced by oofem::QTrPlaneStress2dXFEM::giveInputRecord(), oofem::PlaneStress2dXfem::giveInputRecord(), and oofem::TrPlaneStress2dXFEM::giveInputRecord().

void oofem::XfemStructuralElementInterface::giveIntersectionsTouchingCrack	(	std::vector< int > &	oTouchingEnrItemIndices,
		const std::vector< int > &	iCandidateIndices,
		int	iEnrItemIndex,
		XfemManager &	iXMan
	)

Identify enrichment items with an intersection enrichment front that touches a given enrichment item.

Definition at line 1675 of file xfemstructuralelementinterface.C.

References oofem::EnrichmentItem::giveEnrichmentFrontEnd(), oofem::EnrichmentItem::giveEnrichmentFrontStart(), oofem::XfemManager::giveEnrichmentItem(), oofem::EnrichmentFront::giveTipInfo(), and oofem::EnrichmentItem::tipIsTouchingEI().

Referenced by XfemElementInterface_updateIntegrationRule().

void oofem::XfemStructuralElementInterface::giveSubtriangulationCompositeExportData	(	std::vector< VTKPiece > &	vtkPieces,
		IntArray &	primaryVarsToExport,
		IntArray &	internalVarsToExport,
		IntArray	cellVarsToExport,
		TimeStep *	tStep
	)

VTK Interface.

Definition at line 1714 of file xfemstructuralelementinterface.C.

References oofem::FloatArray::add(), oofem::IntArray::at(), oofem::FloatArray::at(), oofem::FloatArray::beDifferenceOf(), oofem::FloatArray::beProductOf(), oofem::VTKXMLExportModule::computeIPAverage(), computeIPAverageInTriangle(), oofem::Element::computeLocalCoordinates(), oofem::BasicGeometry::computeTangentialSignDist(), oofem::Element::computeVectorOf(), oofem::XfemElementInterface::element, oofem::EnrichmentItem::evalLevelSetNormalInNode(), oofem::EnrichmentItem::evalLevelSetTangInNode(), oofem::FEInterpolation::evalN(), oofem::EnrichmentItem::evalNodeEnrMarkerInNode(), oofem::DofManager::giveCoordinates(), oofem::Element::giveDofManager(), oofem::FEMComponent::giveDomain(), oofem::XfemManager::giveEnrichmentItem(), oofem::GeometryBasedEI::giveGeometry(), oofem::DofManager::giveGlobalNumber(), oofem::Element::giveIntegrationRule(), oofem::Element::giveInterpolation(), oofem::Element::giveNumberOfDofManagers(), oofem::XfemManager::giveNumberOfEnrichmentItems(), oofem::XfemTolerances::giveRelLengthTolTight(), oofem::IntArray::giveSize(), oofem::FloatArray::giveSize(), oofem::Domain::giveXfemManager(), oofem::Domain::hasXfemManager(), mSubTri, N, oofem::FloatArray::push_back(), oofem::FloatArray::resize(), oofem::FloatArray::times(), oofem::XfemManager::vtkExportFields, and oofem::XfemElementInterface::XfemElementInterface_createEnrNmatrixAt().

Referenced by oofem::PlaneStress2dXfem::giveCompositeExportData(), oofem::QTrPlaneStress2dXFEM::giveCompositeExportData(), and oofem::TrPlaneStress2dXFEM::giveCompositeExportData().

virtual bool oofem::XfemStructuralElementInterface::hasCohesiveZone ( ) const

inlinevirtual

Definition at line 66 of file xfemstructuralelementinterface.h.

Referenced by computeCohesiveForces(), computeCohesiveTangent(), and computeCohesiveTangentAt().

IRResultType oofem::XfemStructuralElementInterface::initializeCZFrom ( InputRecord * ir )

virtual

Definition at line 1568 of file xfemstructuralelementinterface.C.

References _IFT_XfemElementInterface_CohesiveZoneMaterial, _IFT_XfemElementInterface_NumIntPointsCZ, _IFT_XfemElementInterface_PlaneStrain, IR_GIVE_OPTIONAL_FIELD, oofem::IRRT_OK, mCSNumGaussPoints, mCZMaterialNum, and oofem::XfemElementInterface::mUsePlaneStrain.

Referenced by oofem::QTrPlaneStress2dXFEM::initializeFrom(), oofem::PlaneStress2dXfem::initializeFrom(), and oofem::TrPlaneStress2dXFEM::initializeFrom().

void oofem::XfemStructuralElementInterface::initializeCZMaterial ( )

virtual

Definition at line 1603 of file xfemstructuralelementinterface.C.

References oofem::XfemElementInterface::element, oofem::FEMComponent::giveDomain(), oofem::Domain::giveMaterial(), mCZMaterialNum, mpCZMat, and OOFEM_ERROR.

Referenced by oofem::PlaneStress2dXfem::postInitialize(), oofem::TrPlaneStress2dXFEM::postInitialize(), oofem::QTrPlaneStress2dXFEM::postInitialize(), and XfemElementInterface_updateIntegrationRule().

bool oofem::XfemStructuralElementInterface::useNonStdCz ( )

Definition at line 1615 of file xfemstructuralelementinterface.C.

References oofem::XfemElementInterface::element, oofem::FEMComponent::giveDomain(), oofem::XfemStructureManager::giveUseNonStdCz(), oofem::Domain::giveXfemManager(), and oofem::Domain::hasXfemManager().

Referenced by computeCohesiveForces(), and computeCohesiveTangent().

void oofem::XfemStructuralElementInterface::XfemElementInterface_computeConsistentMassMatrix	(	FloatMatrix &	answer,
		TimeStep *	tStep,
		double &	mass,
		const double *	ipDensity = `NULL`
	)

virtual

Definition at line 1500 of file xfemstructuralelementinterface.C.

References oofem::IntArray::at(), oofem::FloatMatrix::at(), oofem::StructuralElement::computeNmatrixAt(), oofem::Element::computeNumberOfDofs(), oofem::Element::computeVolumeAround(), oofem::XfemElementInterface::element, oofem::CrossSection::give(), oofem::Element::giveIntegrationRule(), oofem::StructuralElement::giveMassMtrxIntegrationgMask(), oofem::FloatMatrix::giveNumberOfRows(), oofem::StructuralElement::giveStructuralCrossSection(), oofem::Element::isActivated(), oofem::IntArray::isEmpty(), OOFEM_ERROR, oofem::FloatMatrix::plusProductSymmUpper(), oofem::FloatMatrix::resize(), oofem::FloatMatrix::symmetrized(), and oofem::FloatMatrix::zero().

Referenced by oofem::PlaneStress2dXfem::computeConsistentMassMatrix(), oofem::QTrPlaneStress2dXFEM::computeConsistentMassMatrix(), and oofem::TrPlaneStress2dXFEM::computeConsistentMassMatrix().

void oofem::XfemStructuralElementInterface::XfemElementInterface_computeConstitutiveMatrixAt	(	FloatMatrix &	answer,
		MatResponseMode	rMode,
		GaussPoint *	gp,
		TimeStep *	tStep
	)

virtual

Definition at line 797 of file xfemstructuralelementinterface.C.

References oofem::XfemElementInterface::element, oofem::Element::giveCrossSection(), oofem::FEMComponent::giveDomain(), oofem::XfemManager::giveEnrichmentItem(), oofem::XfemManager::giveMaterialModifyingEnrItemIndices(), oofem::StructuralCrossSection::giveStiffnessMatrix_PlaneStrain(), oofem::StructuralCrossSection::giveStiffnessMatrix_PlaneStress(), oofem::Domain::giveXfemManager(), oofem::Domain::hasXfemManager(), oofem::EnrichmentItem::isMaterialModified(), oofem::XfemElementInterface::mUsePlaneStrain, and OOFEM_ERROR.

Referenced by oofem::PlaneStress2dXfem::computeConstitutiveMatrixAt(), oofem::TrPlaneStress2dXFEM::computeConstitutiveMatrixAt(), and oofem::QTrPlaneStress2dXFEM::computeConstitutiveMatrixAt().

void oofem::XfemStructuralElementInterface::XfemElementInterface_computeDeformationGradientVector	(	FloatArray &	answer,
		GaussPoint *	gp,
		TimeStep *	tStep
	)

virtual

Definition at line 1639 of file xfemstructuralelementinterface.C.

References oofem::__MaterialModeToString(), oofem::FloatArray::at(), oofem::FloatArray::beProductOf(), oofem::NLStructuralElement::computeBHmatrixAt(), oofem::Element::computeVectorOf(), oofem::XfemElementInterface::element, oofem::GaussPoint::giveMaterialMode(), oofem::StructuralElement::initialDisplacements, OOFEM_ERROR, and oofem::FloatArray::subtract().

Referenced by oofem::PlaneStress2dXfem::computeDeformationGradientVector(), oofem::TrPlaneStress2dXFEM::computeDeformationGradientVector(), and oofem::QTrPlaneStress2dXFEM::computeDeformationGradientVector().

void oofem::XfemStructuralElementInterface::XfemElementInterface_computeStressVector	(	FloatArray &	answer,
		const FloatArray &	strain,
		GaussPoint *	gp,
		TimeStep *	tStep
	)