oofem::CoupledFieldsElement Class Reference

#include <coupledfieldselement.h>

Inheritance diagram for oofem::CoupledFieldsElement:

Collaboration diagram for oofem::CoupledFieldsElement:

Public Member Functions
	CoupledFieldsElement (int i, Domain *aDomain)
virtual	~CoupledFieldsElement ()
Public Member Functions inherited from oofem::NLStructuralElement
	NLStructuralElement (int n, Domain *d)
virtual	~NLStructuralElement ()
	Destructor.
int	giveGeometryMode ()
void	computeFirstPKStressVector (FloatArray &answer, GaussPoint *gp, TimeStep *tStep)
void	computeCauchyStressVector (FloatArray &answer, GaussPoint *gp, TimeStep *tStep)
void	computeStiffnessMatrix (FloatMatrix &answer, MatResponseMode rMode, TimeStep *tStep) override
void	computeInitialStressMatrix (FloatMatrix &answer, TimeStep *tStep) override
void	computeStiffnessMatrix_withIRulesAsSubcells (FloatMatrix &answer, MatResponseMode rMode, TimeStep *tStep)
void	giveInternalForcesVector (FloatArray &answer, TimeStep *tStep, int useUpdatedGpRecord=0) override
virtual void	computeConstitutiveMatrix_dPdF_At (FloatMatrix &answer, MatResponseMode rMode, GaussPoint *gp, TimeStep *tStep)=0
void	giveInternalForcesVector_withIRulesAsSubcells (FloatArray &answer, TimeStep *tStep, int useUpdatedGpRecord=0) override
virtual void	computeDeformationGradientVector (FloatArray &answer, GaussPoint *gp, TimeStep *tStep)
double	computeCurrentVolume (TimeStep *tStep)
void	initializeFrom (InputRecord &ir, int priority) override
void	giveInputRecord (DynamicInputRecord &input) override
const char *	giveClassName () const override
Public Member Functions inherited from oofem::StructuralElement
	StructuralElement (int n, Domain *d)
virtual	~StructuralElement ()
	Destructor.
void	giveCharacteristicMatrix (FloatMatrix &answer, CharType, TimeStep *tStep) override
void	giveCharacteristicVector (FloatArray &answer, CharType type, ValueModeType mode, TimeStep *tStep) override
virtual void	computeMassMatrix (FloatMatrix &answer, TimeStep *tStep)
virtual void	computeLumpedMassMatrix (FloatMatrix &answer, TimeStep *tStep)
virtual void	computeConsistentMassMatrix (FloatMatrix &answer, TimeStep *tStep, double &mass, const double *ipDensity=NULL)
virtual void	giveMassMtrxIntegrationgMask (IntArray &answer)
void	computeStiffnessMatrix_withIRulesAsSubcells (FloatMatrix &answer, MatResponseMode rMode, TimeStep *tStep)
virtual void	computeLumpedInitialStressMatrix (FloatMatrix &answer, TimeStep *tStep)
void	computeField (ValueModeType mode, TimeStep *tStep, const FloatArray &lcoords, FloatArray &answer) override
virtual void	computeStrainVector (FloatArray &answer, GaussPoint *gp, TimeStep *tStep)
int	giveIPValue (FloatArray &answer, GaussPoint *gp, InternalStateType type, TimeStep *tStep) override
virtual void	computeResultingIPTemperatureAt (FloatArray &answer, TimeStep *tStep, GaussPoint *gp, ValueModeType mode)
virtual void	computeResultingIPEigenstrainAt (FloatArray &answer, TimeStep *tStep, GaussPoint *gp, ValueModeType mode)
virtual void	computeStressVector (FloatArray &answer, const FloatArray &strain, GaussPoint *gp, TimeStep *tStep)=0
virtual void	computeBmatrixAt (GaussPoint *gp, FloatMatrix &answer, int lowerIndx=1, int upperIndx=ALL_STRAINS)=0
virtual void	computeNmatrixAt (const FloatArray &iLocCoord, FloatMatrix &answer)
void	updateBeforeNonlocalAverage (TimeStep *tStep) override
virtual void	giveNonlocalLocationArray (IntArray &locationArray, const UnknownNumberingScheme &us)
virtual void	addNonlocalStiffnessContributions (SparseMtrx &dest, const UnknownNumberingScheme &s, TimeStep *tStep)
int	adaptiveUpdate (TimeStep *tStep) override
void	updateInternalState (TimeStep *tStep) override
void	updateYourself (TimeStep *tStep) override
int	checkConsistency () override
void	giveInputRecord (DynamicInputRecord &input) override
const char *	giveClassName () const override
int	giveInternalStateAtNode (FloatArray &answer, InternalStateType type, InternalStateMode mode, int node, TimeStep *tStep) override
void	showSparseMtrxStructure (CharType mtrx, oofegGraphicContext &gc, TimeStep *tStep) override
	Shows sparse structure.
void	showExtendedSparseMtrxStructure (CharType mtrx, oofegGraphicContext &gc, TimeStep *tStep) override
	Shows extended sparse structure (for example, due to nonlocal interactions for tangent stiffness).
void	computeLoadVector (FloatArray &answer, BodyLoad *load, CharType type, ValueModeType mode, TimeStep *tStep) override
void	computeBoundarySurfaceLoadVector (FloatArray &answer, BoundaryLoad *load, int boundary, CharType type, ValueModeType mode, TimeStep *tStep, bool global=true) override
void	computeBoundaryEdgeLoadVector (FloatArray &answer, BoundaryLoad *load, int boundary, CharType type, ValueModeType mode, TimeStep *tStep, bool global=true) override
virtual void	computeEdgeNMatrix (FloatMatrix &answer, int boundaryID, const FloatArray &lcoords)
	computes edge interpolation matrix
virtual void	computeSurfaceNMatrix (FloatMatrix &answer, int boundaryID, const FloatArray &lcoords)
virtual void	computeConstitutiveMatrixAt (FloatMatrix &answer, MatResponseMode rMode, GaussPoint *gp, TimeStep *tStep)=0
StructuralCrossSection *	giveStructuralCrossSection ()
	Helper function which returns the structural cross-section for the element.
virtual void	createMaterialStatus ()
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	computeTangentFromSurfaceLoad (FloatMatrix &answer, BoundaryLoad *load, int boundary, MatResponseMode rmode, TimeStep *tStep)
virtual void	computeTangentFromEdgeLoad (FloatMatrix &answer, BoundaryLoad *load, int boundary, MatResponseMode rmode, TimeStep *tStep)
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	computeNumberOfDofs ()
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 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	giveGeometryType () const =0
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 ()
virtual int	testElementExtension (ElementExtension ext)
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 int	computeGlobalCoordinates (FloatArray &answer, const FloatArray &lcoords)
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	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	giveInternalStateAtSide (FloatArray &answer, InternalStateType type, InternalStateMode mode, int side, TimeStep *tStep)
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	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.
virtual Interface *	giveInterface (InterfaceType t)
std::string	errorInfo (const char *func) const
	Returns string for prepending output (used by error reporting macros).

Protected Member Functions
virtual void	computeNStressAt (GaussPoint *, FloatArray &)=0
virtual void	computeBStressAt (GaussPoint *, FloatArray &)=0
double	computeVolumeAround (GaussPoint *) override=0
void	computeStiffnessMatrix (FloatMatrix &, MatResponseMode, TimeStep *) override=0
void	giveInternalForcesVector (FloatArray &answer, TimeStep *tStep, int useUpdatedGpRecord) override=0
void	computeVectorOfDofIDs (const IntArray &dofIdArray, ValueModeType valueMode, TimeStep *stepN, FloatArray &answer)
void	computeLocationArrayOfDofIDs (const IntArray &dofIdArray, IntArray &answer)
void	computeStiffnessMatrixGen (FloatMatrix &answer, MatResponseMode rMode, TimeStep *tStep, void(*Nfunc)(GaussPoint *, FloatMatrix), void(*Bfunc)(GaussPoint *, FloatMatrix), void(*NStiffness)(FloatMatrix, MatResponseMode, GaussPoint *, TimeStep *), void(*BStiffness)(FloatMatrix, MatResponseMode, GaussPoint *, TimeStep *), double(*volumeAround)(GaussPoint *))
void	giveInternalForcesVectorGen (FloatArray &answer, TimeStep *tStep, int useUpdatedGpRecord, void(*Nfunc)(GaussPoint *, FloatMatrix), void(*Bfunc)(GaussPoint *, FloatMatrix, int, int), void(*NStress)(GaussPoint *, FloatArray), void(*BStress)(GaussPoint *, FloatArray), double(*volumeAround)(GaussPoint *))
Protected Member Functions inherited from oofem::NLStructuralElement
int	checkConsistency () override
virtual void	computeBHmatrixAt (GaussPoint *gp, FloatMatrix &answer)
virtual void	computePointLoadVectorAt (FloatArray &answer, Load *load, TimeStep *tStep, ValueModeType mode, bool global=true)
virtual void	giveEdgeDofMapping (IntArray &answer, int iEdge) const
virtual void	giveSurfaceDofMapping (IntArray &answer, int iSurf) const
virtual double	computeEdgeVolumeAround (GaussPoint *gp, int iEdge)
virtual double	computeSurfaceVolumeAround (GaussPoint *gp, int iSurf)
virtual int	computeLoadGToLRotationMtrx (FloatMatrix &answer)
virtual int	computeLoadLEToLRotationMatrix (FloatMatrix &answer, int iEdge, GaussPoint *gp)
virtual int	computeLoadLSToLRotationMatrix (FloatMatrix &answer, int iSurf, GaussPoint *gp)
virtual int	giveNumberOfIPForMassMtrxIntegration ()
void	condense (FloatMatrix *stiff, FloatMatrix *mass, FloatArray *load, IntArray *what)
virtual void	setupIRForMassMtrxIntegration (IntegrationRule &iRule)
virtual void	computeBodyLoadVectorAt (FloatArray &answer, Load *load, TimeStep *tStep, ValueModeType mode)
virtual void	computeGaussPoints ()

Protected Attributes
int	nlGeo
Protected Attributes inherited from oofem::NLStructuralElement
int	nlGeometry =0
	Flag indicating if geometrical nonlinearities apply.
Protected Attributes inherited from oofem::StructuralElement
std::unique_ptr< FloatArray >	initialDisplacements
	Initial displacement vector, describes the initial nodal displacements when element has been casted.
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
Static Protected Attributes inherited from oofem::NLStructuralElement
static ParamKey	IPK_NLStructuralElement_nlgeoflag

Detailed Description

Abstract class for gradient formulation of coupled damage-plasticity model(GradDp). Yield function is formulated in the effective stress space and damage is driven by the nonlocal(over-nonlocal) cumulated plastic strain. The new nonlocal degrees of freedom (with the meaning of the nonlocal cumulated plastic strain) are introduced with lower order of approximation functions than the displacement field to avoid spurious stress oscillations

Definition at line 47 of file coupledfieldselement.h.

Constructor & Destructor Documentation

CoupledFieldsElement()

oofem::CoupledFieldsElement::CoupledFieldsElement	(	int	i,
		Domain *	aDomain )

Definition at line 55 of file coupledfieldselement.C.

References nlGeo, and oofem::NLStructuralElement::NLStructuralElement().

~CoupledFieldsElement()

virtual oofem::CoupledFieldsElement::~CoupledFieldsElement ( )

inlinevirtual

Definition at line 54 of file coupledfieldselement.h.

Member Function Documentation

computeBStressAt()

virtual void oofem::CoupledFieldsElement::computeBStressAt ( GaussPoint * , FloatArray &  )

protectedpure virtual

computeLocationArrayOfDofIDs()

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

protected

Definition at line 62 of file coupledfieldselement.C.

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

computeNStressAt()

virtual void oofem::CoupledFieldsElement::computeNStressAt ( GaussPoint * , FloatArray &  )

protectedpure virtual

computeStiffnessMatrix()

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

overrideprotectedpure virtual

Computes numerically stiffness matrix of receiver. Default implementation computes element stiffness using \( K=\int_v B^{\mathrm{T}} D B \mathrm{d}V \) formulae, where \( B \) is element geometric matrix and \( D \) is material stiffness matrix. No geometrical nonlinearity is taken into account. NUmerical integration procedure uses integrationRulesArray for numerical integration. Support for reduced or selected integration is implemented. The individual integration rules are assumed to correspond to different terms from which the overall matrix is assembled.

If there is one integration rule, the full integration of all coefficients is performed. Otherwise, integration is performed using following rules. Each integration rule can specify start and end strain index of strain vector components for which is valid. It is necessary to ensure that these start and end indexes, dividing geometrical matrix into blocks, are not overlapping and that each strain component is included.

Then stiffness matrix is obtained as summation of integrals \( I_{ij}=\int_v B^{\mathrm{T}}_i D_{ij} B_j \mathrm{d}V \) where \( B_i \) is i-th block of geometrical matrix and \( D_{ij} \) is corresponding constitutive sub-matrix. The geometrical matrix is obtained using computeBmatrixAt service and the constitutive matrix is obtained using computeConstitutiveMatrixAt service. The \( I_{ij} \) integral is evaluated using such integration rule, which is valid for i-th or j-th block and has smaller number of integration points.

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.

Reimplemented from oofem::StructuralElement.

computeStiffnessMatrixGen()

void oofem::CoupledFieldsElement::computeStiffnessMatrixGen ( FloatMatrix & answer, MatResponseMode rMode, TimeStep * tStep, void(* Nfunc )(GaussPoint *, FloatMatrix), void(* Bfunc )(GaussPoint *, FloatMatrix), void(* NStiffness )(FloatMatrix, MatResponseMode, GaussPoint *, TimeStep *), void(* BStiffness )(FloatMatrix, MatResponseMode, GaussPoint *, TimeStep *), double(* volumeAround )(GaussPoint *) )

protected

Definition at line 140 of file coupledfieldselement.C.

References oofem::FloatMatrix::beProductOf(), computeVolumeAround(), oofem::Element::giveCrossSection(), oofem::Element::giveIntegrationRule(), N, oofem::FloatMatrix::plusProductSymmUpper(), oofem::FloatMatrix::plusProductUnsym(), oofem::FloatMatrix::resize(), and oofem::FloatMatrix::symmetrized().

computeVectorOfDofIDs()

void oofem::CoupledFieldsElement::computeVectorOfDofIDs ( const IntArray & dofIdArray, ValueModeType valueMode, TimeStep * stepN, FloatArray & answer )

protected

Definition at line 83 of file coupledfieldselement.C.

References oofem::FloatArray::at(), oofem::IntArray::at(), oofem::Element::giveDofManager(), oofem::DofManager::giveDofWithID(), oofem::IntArray::giveSize(), oofem::Dof::giveUnknown(), oofem::DofManager::hasDofID(), oofem::Element::numberOfDofMans, oofem::FloatArray::resize(), and oofem::FloatArray::zero().

computeVolumeAround()

double oofem::CoupledFieldsElement::computeVolumeAround ( GaussPoint * gp )

overrideprotectedpure virtual

Returns volume related to given integration point. Used typically in subroutines, that perform integration over element volume. Should be implemented by particular elements.

See also

GaussPoint

Parameters

gp	Integration point for which volume is computed.

Returns

Volume for integration point.

Reimplemented from oofem::Element.

Referenced by computeStiffnessMatrixGen(), and giveInternalForcesVectorGen().

giveInternalForcesVector()

void oofem::CoupledFieldsElement::giveInternalForcesVector ( FloatArray & answer, TimeStep * tStep, int useUpdatedGpRecord )

overrideprotectedpure virtual

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::StructuralElement.

giveInternalForcesVectorGen()

void oofem::CoupledFieldsElement::giveInternalForcesVectorGen ( FloatArray & answer, TimeStep * tStep, int useUpdatedGpRecord, void(* Nfunc )(GaussPoint *, FloatMatrix), void(* Bfunc )(GaussPoint *, FloatMatrix, int, int), void(* NStress )(GaussPoint *, FloatArray), void(* BStress )(GaussPoint *, FloatArray), double(* volumeAround )(GaussPoint *) )

protected

Definition at line 105 of file coupledfieldselement.C.

References oofem::FloatArray::add(), oofem::FloatArray::beTProductOf(), computeVolumeAround(), oofem::Element::giveIntegrationRule(), oofem::Element::giveNumberOfIntegrationRules(), and N.

Member Data Documentation

nlGeo

int oofem::CoupledFieldsElement::nlGeo

protected

Definition at line 50 of file coupledfieldselement.h.

Referenced by CoupledFieldsElement().

---

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

• coupledfieldselement.h
• coupledfieldselement.C