#include <structuralinterfaceelementphf.h>

Inheritance diagram for oofem::StructuralInterfaceElementPhF:

Collaboration diagram for oofem::StructuralInterfaceElementPhF:

Public Member Functions
	StructuralInterfaceElementPhF (int n, Domain *d)
void	computeStiffnessMatrix (FloatMatrix &answer, MatResponseMode rMode, TimeStep *tStep) override
void	computeStiffnessMatrix_uu (FloatMatrix &answer, MatResponseMode rMode, TimeStep *tStep)
void	computeStiffnessMatrix_ud (FloatMatrix &answer, MatResponseMode rMode, TimeStep *tStep)
void	computeStiffnessMatrix_dd (FloatMatrix &answer, MatResponseMode rMode, TimeStep *tStep)
void	computeStiffnessMatrix_du (FloatMatrix &answer, MatResponseMode rMode, TimeStep *tStep)
void	giveInternalForcesVector_u (FloatArray &answer, TimeStep *tStep, int useUpdatedGpRecord=0)
void	giveInternalForcesVector_d (FloatArray &answer, TimeStep *tStep, int useUpdatedGpRecord=0)
virtual void	giveDofManDofIDMask_u (IntArray &answer)=0
virtual void	giveDofManDofIDMask_d (IntArray &answer)=0
virtual void	computeLocationArrayOfDofIDs (const IntArray &dofIdArray, IntArray &answer)
virtual void	getLocationArray_u (IntArray &answer)=0
virtual void	getLocationArray_d (IntArray &answer)=0
int	computeNumberOfDofs () override
double	giveInternalLength ()
double	computeFreeEnergy (GaussPoint gp, TimeStep tStep)
virtual double	computeDamageAt (GaussPoint gp, ValueModeType valueMode, TimeStep stepN)
void	computeDisplacementUnknowns (FloatArray &answer, ValueModeType valueMode, TimeStep *stepN)
void	computeDamageUnknowns (FloatArray &answer, ValueModeType valueMode, TimeStep *stepN)
void	computeBStress_u (FloatArray &answer, GaussPoint gp, TimeStep tStep, int useUpdatedGpRecord)
void	computeNStress_d (FloatArray &answer, GaussPoint gp, TimeStep tStep, int useUpdatedGpRecord)
virtual void	computeBd_matrixAt (GaussPoint *, FloatMatrix &)
virtual void	computeNd_matrixAt (const FloatArray &lCoords, FloatMatrix &N)
void	giveInternalForcesVector (FloatArray &answer, TimeStep *tStep, int useUpdatedGpRecord=0) override
virtual void	computeTraction (FloatArray &traction, IntegrationPoint ip, FloatArray &jump, TimeStep tStep)
virtual void	computeCovarBaseVectorsAt (GaussPoint *gp, FloatMatrix &G)=0
void	updateInternalState (TimeStep *tStep) override
void	updateYourself (TimeStep *tStep) override
const char *	giveClassName () const override
virtual void	giveEngTraction (FloatArray &answer, GaussPoint gp, const FloatArray &jump, double damage, TimeStep tStep)
Public Member Functions inherited from oofem::StructuralInterfaceElement
	StructuralInterfaceElement (int n, Domain *d)
int	computeGlobalCoordinates (FloatArray &answer, const FloatArray &lcoords) override
void	giveCharacteristicMatrix (FloatMatrix &answer, CharType, TimeStep *tStep) override
void	giveCharacteristicVector (FloatArray &answer, CharType type, ValueModeType mode, TimeStep *tStep) override
virtual void	computeTraction (FloatArray &traction, IntegrationPoint ip, const FloatArray &jump, TimeStep tStep)
virtual void	computeSpatialJump (FloatArray &answer, GaussPoint gp, TimeStep tStep)
int	giveIPValue (FloatArray &answer, GaussPoint gp, InternalStateType type, TimeStep tStep) override
Interface *	giveInterface (InterfaceType) override
int	checkConsistency () override
void	giveInputRecord (DynamicInputRecord &input) override
StructuralInterfaceCrossSection *	giveInterfaceCrossSection ()
virtual double	computeAreaAround (GaussPoint *gp)=0
Element_Geometry_Type	giveGeometryType () const override
virtual void	giveEngTraction (FloatArray &answer, GaussPoint gp, const FloatArray &jump, TimeStep tStep)
virtual void	giveFirstPKTraction (FloatArray &answer, GaussPoint gp, const FloatArray &jump, const FloatMatrix &F, TimeStep tStep)
virtual void	giveStiffnessMatrix_Eng (FloatMatrix &answer, MatResponseMode rMode, IntegrationPoint ip, TimeStep tStep)
virtual void	giveStiffnessMatrix_dTdj (FloatMatrix &answer, MatResponseMode rMode, IntegrationPoint ip, TimeStep tStep)
Public Member Functions inherited from oofem::Element
	Element (int n, Domain *aDomain)
	Element (const Element &src)=delete
Element &	operator= (const Element &src)=delete
virtual	~Element ()
	Virtual destructor.
void	giveLocationArray (IntArray &locationArray, const UnknownNumberingScheme &s, IntArray *dofIds=NULL) const
void	giveLocationArray (IntArray &locationArray, const IntArray &dofIDMask, const UnknownNumberingScheme &s, IntArray *dofIds=NULL) const
virtual void	giveBoundaryLocationArray (IntArray &locationArray, const IntArray &bNodes, const UnknownNumberingScheme &s, IntArray *dofIds=NULL)
virtual void	giveBoundaryLocationArray (IntArray &locationArray, const IntArray &bNodes, const IntArray &dofIDMask, const UnknownNumberingScheme &s, IntArray *dofIds=NULL)
virtual int	giveNumberOfDofs ()
virtual int	giveNumberOfInternalDofManagers () const
virtual DofManager *	giveInternalDofManager (int i) const
virtual void	setInternalDofManager (int num, std::unique_ptr< DofManager > dm)
virtual double	giveCharacteristicValue (CharType type, TimeStep *tStep)
virtual void	computeLoadVector (FloatArray &answer, BodyLoad load, CharType type, ValueModeType mode, TimeStep tStep)
virtual void	computeBoundarySurfaceLoadVector (FloatArray &answer, BoundaryLoad load, int boundary, CharType type, ValueModeType mode, TimeStep tStep, bool global=true)
virtual void	computeTangentFromSurfaceLoad (FloatMatrix &answer, BoundaryLoad load, int boundary, MatResponseMode rmode, TimeStep tStep)
virtual void	computeTangentFromEdgeLoad (FloatMatrix &answer, BoundaryLoad load, int boundary, MatResponseMode rmode, TimeStep tStep)
virtual void	computeBoundaryEdgeLoadVector (FloatArray &answer, BoundaryLoad load, int edge, CharType type, ValueModeType mode, TimeStep tStep, bool global=true)
const IntArray &	giveBodyLoadList () const
const IntArray &	giveBoundaryLoadList () const
void	computeVectorOf (ValueModeType u, TimeStep *tStep, FloatArray &answer)
void	computeVectorOf (const IntArray &dofIDMask, ValueModeType u, TimeStep *tStep, FloatArray &answer, bool padding=false)
void	computeBoundaryVectorOf (const IntArray &bNodes, const IntArray &dofIDMask, ValueModeType u, TimeStep *tStep, FloatArray &answer, bool padding=false)
void	computeVectorOf (PrimaryField &field, const IntArray &dofIDMask, ValueModeType u, TimeStep *tStep, FloatArray &answer, bool padding=false)
void	computeVectorOfPrescribed (ValueModeType u, TimeStep *tStep, FloatArray &answer)
void	computeVectorOfPrescribed (const IntArray &dofIDMask, ValueModeType type, TimeStep *tStep, FloatArray &answer)
virtual int	computeNumberOfGlobalDofs ()
int	computeNumberOfPrimaryMasterDofs ()
virtual bool	computeGtoLRotationMatrix (FloatMatrix &answer)
virtual bool	giveRotationMatrix (FloatMatrix &answer)
virtual bool	computeDofTransformationMatrix (FloatMatrix &answer, const IntArray &nodes, bool includeInternal)
virtual void	giveDofManDofIDMask (int inode, IntArray &answer) const
virtual void	giveInternalDofManDofIDMask (int inode, IntArray &answer) const
virtual void	giveElementDofIDMask (IntArray &answer) const
virtual void	computeField (ValueModeType mode, TimeStep *tStep, const FloatArray &lcoords, FloatArray &answer)
virtual double	computeVolumeAround (GaussPoint *gp)
virtual double	computeVolumeAreaOrLength ()
	Computes the volume, area or length of the element depending on its spatial dimension.
double	computeMeanSize ()
virtual double	computeVolume ()
virtual double	computeArea ()
virtual double	computeLength ()
virtual IntArray	giveBoundaryEdgeNodes (int boundary, bool includeHierarchical=false) const
virtual IntArray	giveBoundarySurfaceNodes (int boundary, bool includeHierarchical=false) const
virtual IntArray	giveBoundaryNodes (int boundary) const
virtual std::unique_ptr< IntegrationRule >	giveBoundaryEdgeIntegrationRule (int order, int boundary)
virtual std::unique_ptr< IntegrationRule >	giveBoundarySurfaceIntegrationRule (int order, int boundary)
int	giveDofManagerNumber (int i) const
const IntArray &	giveDofManArray () const
void	addDofManager (DofManager *dMan)
DofManager *	giveDofManager (int i) const
Node *	giveNode (int i) const
virtual ElementSide *	giveSide (int i) const
virtual FEInterpolation *	giveInterpolation () const
virtual FEInterpolation *	giveInterpolation (DofIDItem id) const
virtual const FEInterpolation *	getGeometryInterpolation () const
virtual Material *	giveMaterial ()
int	giveMaterialNumber () const
CrossSection *	giveCrossSection ()
int	getActivityTimeFunctionNumber ()
void	setActivityTimeFunctionNumber (int funcIndx)
void	setMaterial (int matIndx)
virtual void	setCrossSection (int csIndx)
virtual int	giveNumberOfDofManagers () const
void	setNumberOfDofManagers (int i)
	Sets number of element dof managers.
virtual int	giveNumberOfNodes () const
void	setDofManagers (const IntArray &dmans)
void	setDofManager (int id, int dm)
void	setBodyLoads (const IntArray &bodyLoads)
void	setIntegrationRules (std ::vector< std ::unique_ptr< IntegrationRule > > irlist)
virtual integrationDomain	giveIntegrationDomain () const
virtual MaterialMode	giveMaterialMode ()
virtual int	giveIntegrationRuleLocalCodeNumbers (IntArray &answer, IntegrationRule &ie)
int	giveRegionNumber ()
virtual void	initializeYourself (TimeStep *timeStepWhenICApply)
virtual bool	isActivated (TimeStep *tStep)
virtual bool	isCast (TimeStep *tStep)
virtual void	initForNewStep ()
virtual Element_Geometry_Type	giveEdgeGeometryType (int id) const
	Returns the receiver edge geometry type.
virtual Element_Geometry_Type	giveSurfaceGeometryType (int id) const
	Returns the receiver surface geometry type.
virtual int	giveSpatialDimension ()
virtual int	giveNumberOfBoundarySides ()
	Returns number of boundaries (entities of element_dimension-1: points, edges, surfaces).
virtual int	giveNumberOfEdges () const
virtual int	giveNumberOfSurfaces () const
virtual int	giveDefaultIntegrationRule () const
virtual IntegrationRule *	giveDefaultIntegrationRulePtr ()
int	giveNumberOfIntegrationRules ()
virtual IntegrationRule *	giveIntegrationRule (int i)
std::vector< std ::unique_ptr< IntegrationRule > > &	giveIntegrationRulesArray ()
int	giveGlobalIPValue (FloatArray &answer, GaussPoint gp, InternalStateType type, TimeStep tStep)
virtual double	giveLengthInDir (const FloatArray &normalToCrackPlane)
virtual double	giveCharacteristicLength (const FloatArray &normalToCrackPlane)
double	giveCharacteristicLengthForPlaneElements (const FloatArray &normalToCrackPlane)
double	giveCharacteristicLengthForAxisymmElements (const FloatArray &normalToCrackPlane)
virtual double	giveCharacteristicSize (GaussPoint *gp, FloatArray &normalToCrackPlane, ElementCharSizeMethod method)
virtual double	giveParentElSize () const
virtual void	updateBeforeNonlocalAverage (TimeStep *tStep)
virtual bool	computeLocalCoordinates (FloatArray &answer, const FloatArray &gcoords)
virtual int	giveLocalCoordinateSystem (FloatMatrix &answer)
virtual void	giveLocalCoordinateSystemVector (InternalStateType isttype, FloatArray &answer)
virtual void	computeMidPlaneNormal (FloatArray &answer, const GaussPoint *gp)
virtual int	adaptiveMap (Domain oldd, TimeStep tStep)
virtual int	mapStateVariables (Domain &iOldDom, const TimeStep &iTStep)
virtual int	adaptiveUpdate (TimeStep *tStep)
virtual int	adaptiveFinish (TimeStep *tStep)
void	updateLocalNumbering (EntityRenumberingFunctor &f) override
template<class T>
void	ipEvaluator (T src, void(T ::f)(GaussPoint *gp))
	Integration point evaluator, loops over receiver IP's and calls given function (passed as f parameter) on them. The IP is parameter to function f.
template<class T, class S>
void	ipEvaluator (T src, void(T ::f)(GaussPoint *, S &), S &_val)
	Integration point evaluator, loops over receiver IP's and calls given function (passed as f parameter) on them. The IP is parameter to function f as well as additional array.
virtual void	drawYourself (oofegGraphicContext &gc, TimeStep *tStep)
virtual void	drawAnnotation (oofegGraphicContext &gc, TimeStep *tStep)
virtual void	drawRawGeometry (oofegGraphicContext &gc, TimeStep *tStep)
virtual void	drawDeformedGeometry (oofegGraphicContext &gc, TimeStep *tStep, UnknownType)
virtual void	drawScalar (oofegGraphicContext &gc, TimeStep *tStep)
virtual void	drawSpecial (oofegGraphicContext &gc, TimeStep *tStep)
virtual void	giveLocalIntVarMaxMin (oofegGraphicContext &gc, TimeStep *tStep, double &emin, double &emax)
virtual int	giveInternalStateAtNode (FloatArray &answer, InternalStateType type, InternalStateMode mode, int node, TimeStep *tStep)
virtual int	giveInternalStateAtSide (FloatArray &answer, InternalStateType type, InternalStateMode mode, int side, TimeStep *tStep)
virtual void	showSparseMtrxStructure (CharType mtrx, oofegGraphicContext &gc, TimeStep *tStep)
	Shows sparse structure.
virtual void	showExtendedSparseMtrxStructure (CharType mtrx, oofegGraphicContext &gc, TimeStep *tStep)
	Shows extended sparse structure (for example, due to nonlocal interactions for tangent stiffness).
int	giveLabel () const
int	giveGlobalNumber () const
void	setGlobalNumber (int num)
elementParallelMode	giveParallelMode () const
void	setParallelMode (elementParallelMode _mode)
	Sets parallel mode of element.
virtual elementParallelMode	giveKnotSpanParallelMode (int) const
int	packUnknowns (DataStream &buff, TimeStep *tStep)
int	unpackAndUpdateUnknowns (DataStream &buff, TimeStep *tStep)
int	estimatePackSize (DataStream &buff)
const IntArray *	givePartitionList () const
void	setPartitionList (IntArray &pl)
virtual double	predictRelativeComputationalCost ()
virtual double	giveRelativeSelfComputationalCost ()
virtual double	predictRelativeRedistributionCost ()
IntArray *	giveBodyLoadArray ()
	Returns array containing load numbers of loads acting on element.
IntArray *	giveBoundaryLoadArray ()
	Returns array containing load numbers of boundary loads acting on element.
void	initializeFrom (InputRecord &ir, int priority) override
void	initializeFinish () override
void	postInitialize () override
	Performs post initialization steps.
void	saveContext (DataStream &stream, ContextMode mode) override
void	restoreContext (DataStream &stream, ContextMode mode) override
void	printOutputAt (FILE file, TimeStep tStep) override
virtual const IntArray	giveLocation ()
virtual void	recalculateCoordinates (int nodeNumber, FloatArray &coords)
void	setSharedEdgeID (int iedge, int globalID)
void	setSharedSurfaceID (int isurf, int globalID)
const IntArray *	giveSharedEdgeIDs () const
const IntArray *	giveSharedSurfaceIDs () const
Public Member Functions inherited from oofem::FEMComponent
	FEMComponent (int n, Domain *d)
virtual	~FEMComponent ()=default
	Virtual destructor.
virtual const char *	giveInputRecordName () const =0
Domain *	giveDomain () const
virtual void	setDomain (Domain *d)
int	giveNumber () const
void	setNumber (int num)
virtual void	initializeFrom (InputRecord &ir)
virtual void	printYourself ()
	Prints receiver state on stdout. Useful for debugging.
std::string	errorInfo (const char *func) const
	Returns string for prepending output (used by error reporting macros).

Public Attributes
double	internalLength

Protected Attributes
FEInterpolation *	interpolation = nullptr
IntArray	loc_u
IntArray	loc_d
Protected Attributes inherited from oofem::StructuralInterfaceElement
FloatArray	initialDisplacements
	Initial displacement vector, describes the initial nodal displacements when element has been casted.
bool	nlGeometry = false
	Flag indicating if geometrical nonlinearities apply.
Protected Attributes inherited from oofem::Element
int	numberOfDofMans
	Number of dofmanagers.
IntArray	dofManArray
	Array containing dofmanager numbers.
int	material
	Number of associated material.
int	crossSection
	Number of associated cross section.
IntArray	bodyLoadArray
IntArray	boundaryLoadArray
std::vector< std ::unique_ptr< IntegrationRule > >	integrationRulesArray
FloatMatrix	elemLocalCS
	Transformation material matrix, used in orthotropic and anisotropic materials, global->local transformation.
int	activityTimeFunction
	Element activity time function. If defined, nonzero value indicates active receiver, zero value inactive element.
int	globalNumber
int	numberOfGaussPoints
elementParallelMode	parallel_mode
	Determines the parallel mode of the element.
IntArray	partitions
IntArray	globalEdgeIDs
IntArray	globalSurfaceIDs
Protected Attributes inherited from oofem::FEMComponent
int	number
	Component number.
Domain *	domain
	Link to domain object, useful for communicating with other FEM components.

Additional Inherited Members
Static Public Attributes inherited from oofem::Element
static ParamKey	IPK_Element_mat
static ParamKey	IPK_Element_crosssect
static ParamKey	IPK_Element_nodes
static ParamKey	IPK_Element_bodyload
static ParamKey	IPK_Element_boundaryload
static ParamKey	IPK_Element_lcs
static ParamKey	IPK_Element_partitions
static ParamKey	IPK_Element_remote
static ParamKey	IPK_Element_activityTimeFunction
static ParamKey	IPK_Element_nip
Protected Member Functions inherited from oofem::StructuralInterfaceElement
virtual void	computeNmatrixAt (GaussPoint *gp, FloatMatrix &answer)=0
virtual void	computeTransformationMatrixAt (GaussPoint *gp, FloatMatrix &answer)=0
virtual int	giveNumberOfIPForMassMtrxIntegration ()
virtual int	testCrossSectionExtension (CrossSectExtension ext)
int	testElementExtension (ElementExtension ext) override
virtual void	computeGaussPoints ()

Detailed Description

Interface element class with phase field (PhF) modeling of damage

Definition at line 60 of file structuralinterfaceelementphf.h.

Constructor & Destructor Documentation

◆ StructuralInterfaceElementPhF()

oofem::StructuralInterfaceElementPhF::StructuralInterfaceElementPhF	(	int	n,
		Domain *	d )

Constructor. Creates structural element with given number, belonging to given domain.

Parameters

n	Element number.
d	Domain to which new material will belong.

Definition at line 53 of file structuralinterfaceelementphf.C.

References internalLength, and oofem::StructuralInterfaceElement::StructuralInterfaceElement().

Referenced by oofem::IntElLine1PhF::IntElLine1PhF().

Member Function Documentation

◆ computeBd_matrixAt()

void oofem::StructuralInterfaceElementPhF::computeBd_matrixAt	(	GaussPoint *	gp,
		FloatMatrix &	answer )

virtual

Definition at line 213 of file structuralinterfaceelementphf.C.

References oofem::FloatMatrix::beProductTOf(), computeCovarBaseVectorsAt(), oofem::Element::giveInterpolation(), and oofem::GaussPoint::giveNaturalCoordinates().

Referenced by computeStiffnessMatrix_dd(), and giveInternalForcesVector_d().

◆ computeBStress_u()

void oofem::StructuralInterfaceElementPhF::computeBStress_u	(	FloatArray &	answer,
		GaussPoint *	gp,
		TimeStep *	tStep,
		int	useUpdatedGpRecord )

◆ computeCovarBaseVectorsAt()

virtual void oofem::StructuralInterfaceElementPhF::computeCovarBaseVectorsAt	(	GaussPoint *	gp,
		FloatMatrix &	G )

pure virtual

Implemented in oofem::IntElLine1PhF.

Referenced by computeBd_matrixAt().

◆ computeDamageAt()

double oofem::StructuralInterfaceElementPhF::computeDamageAt	(	GaussPoint *	gp,
		ValueModeType	valueMode,
		TimeStep *	stepN )

virtual

Definition at line 191 of file structuralinterfaceelementphf.C.

References oofem::FloatArray::at(), computeDamageUnknowns(), oofem::FloatArray::dotProduct(), oofem::Element::giveInterpolation(), and oofem::GaussPoint::giveNaturalCoordinates().

Referenced by computeTraction(), and giveInternalForcesVector_d().

◆ computeDamageUnknowns()

void oofem::StructuralInterfaceElementPhF::computeDamageUnknowns	(	FloatArray &	answer,
		ValueModeType	valueMode,
		TimeStep *	stepN )

Definition at line 70 of file structuralinterfaceelementphf.C.

References oofem::Element::computeVectorOf(), and giveDofManDofIDMask_d().

Referenced by computeDamageAt(), and giveInternalForcesVector_d().

◆ computeDisplacementUnknowns()

void oofem::StructuralInterfaceElementPhF::computeDisplacementUnknowns	(	FloatArray &	answer,
		ValueModeType	valueMode,
		TimeStep *	stepN )

Definition at line 61 of file structuralinterfaceelementphf.C.

References oofem::Element::computeVectorOf(), and giveDofManDofIDMask_u().

◆ computeFreeEnergy()

double oofem::StructuralInterfaceElementPhF::computeFreeEnergy	(	GaussPoint *	gp,
		TimeStep *	tStep )

Todo: put in material model

◆ computeLocationArrayOfDofIDs()

void oofem::StructuralInterfaceElementPhF::computeLocationArrayOfDofIDs	(	const IntArray &	dofIdArray,
		IntArray &	answer )

virtual

Definition at line 79 of file structuralinterfaceelementphf.C.

References oofem::IntArray::at(), oofem::IntArray::followedBy(), oofem::Element::giveDofManager(), oofem::Element::giveNumberOfDofManagers(), oofem::DofManager::giveNumberOfDofs(), oofem::IntArray::giveSize(), oofem::DofManager::hasDofID(), and oofem::IntArray::resize().

◆ computeNd_matrixAt()

void oofem::StructuralInterfaceElementPhF::computeNd_matrixAt	(	const FloatArray &	lCoords,
		FloatMatrix &	N )

virtual

Definition at line 204 of file structuralinterfaceelementphf.C.

References oofem::Element::giveInterpolation(), and N.

Referenced by computeStiffnessMatrix_dd(), and computeStiffnessMatrix_ud().

◆ computeNStress_d()

void oofem::StructuralInterfaceElementPhF::computeNStress_d	(	FloatArray &	answer,
		GaussPoint *	gp,
		TimeStep *	tStep,
		int	useUpdatedGpRecord )

References N.

◆ computeNumberOfDofs()

int oofem::StructuralInterfaceElementPhF::computeNumberOfDofs ( )

overridevirtual

Computes or simply returns total number of element's local DOFs. Must be defined by particular element.

Returns: Number of local DOFs of element.

Reimplemented from oofem::Element.

Definition at line 226 of file structuralinterfaceelementphf.C.

References oofem::Element::giveDofManager(), and oofem::Element::giveNumberOfDofManagers().

Referenced by computeStiffnessMatrix(), and giveInternalForcesVector().

◆ computeStiffnessMatrix()

void oofem::StructuralInterfaceElementPhF::computeStiffnessMatrix	(	FloatMatrix &	answer,
		MatResponseMode	rMode,
		TimeStep *	tStep )

overridevirtual

Computes the stiffness/tangent matrix of receiver. Default implementation computes element stiffness using \( K=\int_{\Gamma} N^{\mathrm{T}} D N \mathrm{d}V \) formulae, where \( N \) is the element geometric matrix such that \( j = N a \) and \( D \) is the stiffness matrix of the interface material. Numerical integration procedure uses integrationRulesArray for numerical integration.

The geometrical matrix is obtained using computeNmatrixAt service and the constitutive matrix is obtained using computeConstitutiveMatrixAt service.

For higher numerical performance, only one half of stiffness matrix is computed and answer is then symmetrized. Therefore, if element matrix will be generally nonsymmetric, one must specialize this method. Finally, the result is transformed into global coordinate system (or nodal coordinate system, if it is defined).

Parameters

answer	Computed stiffness matrix (symmetric).
rMode	Response mode.
tStep	Time step.

Todo: this part is enough to do once

Reimplemented from oofem::StructuralInterfaceElement.

Definition at line 238 of file structuralinterfaceelementphf.C.

References oofem::FloatMatrix::assemble(), computeNumberOfDofs(), computeStiffnessMatrix_dd(), computeStiffnessMatrix_uu(), getLocationArray_d(), getLocationArray_u(), loc_d, loc_u, oofem::FloatMatrix::resize(), and oofem::FloatMatrix::zero().

◆ computeStiffnessMatrix_dd()

void oofem::StructuralInterfaceElementPhF::computeStiffnessMatrix_dd	(	FloatMatrix &	answer,
		MatResponseMode	rMode,
		TimeStep *	tStep )

Definition at line 334 of file structuralinterfaceelementphf.C.

References oofem::FloatMatrix::assemble(), oofem::FloatMatrix::clear(), oofem::StructuralInterfaceElement::computeAreaAround(), computeBd_matrixAt(), computeNd_matrixAt(), oofem::Element::giveDefaultIntegrationRulePtr(), oofem::StructuralInterfaceMaterialPhF::giveDrivingForcePrime(), oofem::StructuralInterfaceElement::giveInterfaceCrossSection(), giveInternalLength(), oofem::FloatMatrix::plusProductUnsym(), oofem::FloatMatrix::resize(), and oofem::FloatMatrix::zero().

Referenced by computeStiffnessMatrix().

◆ computeStiffnessMatrix_du()

void oofem::StructuralInterfaceElementPhF::computeStiffnessMatrix_du	(	FloatMatrix &	answer,
		MatResponseMode	rMode,
		TimeStep *	tStep )

◆ computeStiffnessMatrix_ud()

void oofem::StructuralInterfaceElementPhF::computeStiffnessMatrix_ud	(	FloatMatrix &	answer,
		MatResponseMode	rMode,
		TimeStep *	tStep )

Definition at line 307 of file structuralinterfaceelementphf.C.

References oofem::FloatMatrix::beProductOf(), oofem::FloatMatrix::clear(), oofem::StructuralInterfaceElement::computeAreaAround(), computeNd_matrixAt(), oofem::StructuralInterfaceElement::computeNmatrixAt(), oofem::Element::giveDefaultIntegrationRulePtr(), N, oofem::FloatMatrix::plusProductUnsym(), and S.

◆ computeStiffnessMatrix_uu()

void oofem::StructuralInterfaceElementPhF::computeStiffnessMatrix_uu	(	FloatMatrix &	answer,
		MatResponseMode	rMode,
		TimeStep *	tStep )

Definition at line 272 of file structuralinterfaceelementphf.C.

References oofem::FloatMatrix::beProductOf(), oofem::FloatMatrix::clear(), oofem::StructuralInterfaceElement::computeAreaAround(), oofem::StructuralInterfaceElement::computeNmatrixAt(), oofem::StructuralInterfaceElement::computeTransformationMatrixAt(), oofem::Element::giveCrossSection(), oofem::Element::giveDefaultIntegrationRulePtr(), oofem::StructuralInterfaceElement::giveStiffnessMatrix_Eng(), N, oofem::FloatMatrix::plusProductSymmUpper(), oofem::FloatMatrix::plusProductUnsym(), oofem::FloatMatrix::rotatedWith(), and oofem::FloatMatrix::symmetrized().

Referenced by computeStiffnessMatrix().

◆ computeTraction()

void oofem::StructuralInterfaceElementPhF::computeTraction	(	FloatArray &	traction,
		IntegrationPoint *	ip,
		FloatArray &	jump,
		TimeStep *	tStep )

virtual

Definition at line 369 of file structuralinterfaceelementphf.C.

References computeDamageAt(), oofem::StructuralInterfaceElement::computeTransformationMatrixAt(), giveEngTraction(), and oofem::FloatArray::rotatedWith().

Referenced by giveInternalForcesVector_u(), and updateInternalState().

◆ getLocationArray_d()

virtual void oofem::StructuralInterfaceElementPhF::getLocationArray_d ( IntArray & answer )

pure virtual

Implemented in oofem::IntElLine1PhF.

Referenced by computeStiffnessMatrix(), and giveInternalForcesVector().

◆ getLocationArray_u()

virtual void oofem::StructuralInterfaceElementPhF::getLocationArray_u ( IntArray & answer )

pure virtual

Implemented in oofem::IntElLine1PhF.

Referenced by computeStiffnessMatrix(), and giveInternalForcesVector().

◆ giveClassName()

const char * oofem::StructuralInterfaceElementPhF::giveClassName ( ) const

inlineoverridevirtual

Returns: Class name of the receiver.

Reimplemented from oofem::StructuralInterfaceElement.

Definition at line 139 of file structuralinterfaceelementphf.h.

◆ giveDofManDofIDMask_d()

virtual void oofem::StructuralInterfaceElementPhF::giveDofManDofIDMask_d ( IntArray & answer )

pure virtual

Implemented in oofem::IntElLine1PhF.

Referenced by computeDamageUnknowns(), and giveInternalForcesVector().

◆ giveDofManDofIDMask_u()

virtual void oofem::StructuralInterfaceElementPhF::giveDofManDofIDMask_u ( IntArray & answer )

pure virtual

Implemented in oofem::IntElLine1PhF.

Referenced by computeDisplacementUnknowns(), giveInternalForcesVector(), and giveInternalForcesVector_u().

◆ giveEngTraction()

virtual void oofem::StructuralInterfaceElementPhF::giveEngTraction	(	FloatArray &	answer,
		GaussPoint *	gp,
		const FloatArray &	jump,
		double	damage,
		TimeStep *	tStep )

inlinevirtual

Reimplemented in oofem::IntElLine1PhF.

Definition at line 144 of file structuralinterfaceelementphf.h.

References OOFEM_ERROR.

Referenced by computeTraction().

◆ giveInternalForcesVector()

void oofem::StructuralInterfaceElementPhF::giveInternalForcesVector	(	FloatArray &	answer,
		TimeStep *	tStep,
		int	useUpdatedGpRecord = 0 )

overridevirtual

Returns equivalent nodal forces vectors. Useful for nonlinear analysis. Default implementation computes result as \( F=\int_v B^{\mathrm{T}} \sigma \mathrm{d}V \), where \( \sigma \) is the real element stress vector obtained using computeStressVector service (if useUpdatedGpRecord=0) or (if useUpdatedGpRecord=1) from integration point status. The geometric matrix is obtained using computeBmatrixAt service. Integration is performed using default integration rule, which should produce always valid results, assuming that strains used for computation of stresses are valid.

Parameters

answer	Internal nodal forces vector.
tStep	Time step.
useUpdatedGpRecord	If equal to zero, the stresses in integration points are computed (slow but safe), else if nonzero the stresses are taken directly from integration point status (should be derived from StructuralMaterialStatus) (fast, but engineering model must ensure valid status data in each integration point).

Reimplemented from oofem::StructuralInterfaceElement.

Definition at line 99 of file structuralinterfaceelementphf.C.

References oofem::FloatArray::assemble(), computeNumberOfDofs(), getLocationArray_d(), getLocationArray_u(), giveDofManDofIDMask_d(), giveDofManDofIDMask_u(), giveInternalForcesVector_d(), giveInternalForcesVector_u(), loc_d, loc_u, oofem::FloatArray::resize(), and oofem::FloatArray::zero().

◆ giveInternalForcesVector_d()

void oofem::StructuralInterfaceElementPhF::giveInternalForcesVector_d	(	FloatArray &	answer,
		TimeStep *	tStep,
		int	useUpdatedGpRecord = 0 )

◆ giveInternalForcesVector_u()

void oofem::StructuralInterfaceElementPhF::giveInternalForcesVector_u	(	FloatArray &	answer,
		TimeStep *	tStep,
		int	useUpdatedGpRecord = 0 )

Definition at line 120 of file structuralinterfaceelementphf.C.

References oofem::FloatArray::beProductOf(), oofem::FloatArray::clear(), oofem::StructuralInterfaceElement::computeAreaAround(), oofem::StructuralInterfaceElement::computeNmatrixAt(), computeTraction(), oofem::Element::computeVectorOf(), oofem::Element::giveDefaultIntegrationRulePtr(), giveDofManDofIDMask_u(), oofem::StructuralInterfaceElement::initialDisplacements, N, oofem::FloatArray::plusProduct(), and oofem::FloatArray::subtract().

Referenced by giveInternalForcesVector().

◆ giveInternalLength()

double oofem::StructuralInterfaceElementPhF::giveInternalLength ( )

inline

Definition at line 100 of file structuralinterfaceelementphf.h.

References internalLength.

Referenced by computeStiffnessMatrix_dd(), and giveInternalForcesVector_d().

◆ updateInternalState()

void oofem::StructuralInterfaceElementPhF::updateInternalState ( TimeStep * tStep )

overridevirtual

Updates element state after equilibrium in time step has been reached. Default implementation updates all integration rules defined by integrationRulesArray member variable. Doing this, all integration points and their material statuses are updated also. All temporary history variables, which now describe equilibrium state are copied into equilibrium ones. The existing internal state is used for update.

Parameters

tStep Time step for newly reached state.

See also: Material::updateYourself; IntegrationRule::updateYourself; GaussPoint::updateYourself; Element::updateInternalState

Reimplemented from oofem::StructuralInterfaceElement.

Definition at line 397 of file structuralinterfaceelementphf.C.

References oofem::StructuralInterfaceElement::computeSpatialJump(), computeTraction(), and oofem::Element::integrationRulesArray.

◆ updateYourself()

void oofem::StructuralInterfaceElementPhF::updateYourself ( TimeStep * tStep )

overridevirtual

Updates element state after equilibrium in time step has been reached. Default implementation updates all integration rules defined by integrationRulesArray member variable. Doing this, all integration points and their material statuses are updated also. All temporary history variables, which now describe equilibrium state are copied into equilibrium ones. The existing internal state is used for update.

Parameters

tStep Time step for newly reached state.