#include <macrolspace.h>

Inheritance diagram for oofem::MacroLSpace:

Collaboration diagram for oofem::MacroLSpace:

Public Member Functions
	MacroLSpace (int n, Domain *d)
virtual	~MacroLSpace ()
const char *	giveInputRecordName () const override
const char *	giveClassName () const override
void	initializeFrom (InputRecord &ir, int priority) override
void	postInitialize () override
	Performs post initialization steps.
void	computeStiffnessMatrix (FloatMatrix &answer, MatResponseMode rMode, TimeStep *tStep) override
void	computeField (ValueModeType mode, TimeStep *tStep, const FloatArray &lcoords, FloatArray &answer) override
virtual void	changeMicroBoundaryConditions (TimeStep *tStep)
	Related to setting the boundary conditions of micro problem.
void	giveInternalForcesVector (FloatArray &answer, TimeStep *tStep, int useUpdatedGpRecord=0) override
virtual void	evalInterpolation (FloatArray &answer, const std::vector< FloatArray > &coords, const FloatArray &gcoords)
void	updateYourself (TimeStep *tStep) override
Public Member Functions inherited from oofem::LSpace
	LSpace (int n, Domain *d)
virtual	~LSpace ()
FEInterpolation *	giveInterpolation () const override
Interface *	giveInterface (InterfaceType it) override
int	testElementExtension (ElementExtension ext) override
void	SPRNodalRecoveryMI_giveSPRAssemblyPoints (IntArray &pap) override
void	SPRNodalRecoveryMI_giveDofMansDeterminedByPatch (IntArray &answer, int pap) override
int	SPRNodalRecoveryMI_giveNumberOfIP () override
SPRPatchType	SPRNodalRecoveryMI_givePatchType () override
void	NodalAveragingRecoveryMI_computeNodalValue (FloatArray &answer, int node, InternalStateType type, TimeStep *tStep) override
void	HuertaErrorEstimatorI_setupRefinedElementProblem (RefinedElement refinedElement, int level, int nodeId, IntArray &localNodeIdArray, IntArray &globalNodeIdArray, HuertaErrorEstimatorInterface ::SetupMode sMode, TimeStep tStep, int &localNodeId, int &localElemId, int &localBcId, IntArray &controlNode, IntArray &controlDof, HuertaErrorEstimator ::AnalysisMode aMode) override
void	HuertaErrorEstimatorI_computeNmatrixAt (GaussPoint *gp, FloatMatrix &answer) override
Element_Geometry_Type	giveGeometryType () const override
void	drawRawGeometry (oofegGraphicContext &gc, TimeStep *tStep) override
void	drawDeformedGeometry (oofegGraphicContext &gc, TimeStep *tStep, UnknownType type) override
void	drawScalar (oofegGraphicContext &gc, TimeStep *tStep) override
void	drawSpecial (oofegGraphicContext &gc, TimeStep *tStep) override
void	drawTriad (FloatArray &, int isurf)
Public Member Functions inherited from oofem::Structural3DElement
	Structural3DElement (int n, Domain *d)
virtual	~Structural3DElement ()
	Destructor.
void	initializeFrom (InputRecord &ir, int priority) override
MaterialMode	giveMaterialMode () override
int	computeNumberOfDofs () override
void	giveDofManDofIDMask (int inode, IntArray &answer) const override
double	computeVolumeAround (GaussPoint *gp) override
double	giveCharacteristicLength (const FloatArray &normalToCrackPlane) override
void	giveMaterialOrientationAt (FloatArray &x, FloatArray &y, FloatArray &z, const FloatArray &lcoords)
void	computeStressVector (FloatArray &answer, const FloatArray &strain, GaussPoint gp, TimeStep tStep) override
void	computeConstitutiveMatrixAt (FloatMatrix &answer, MatResponseMode rMode, GaussPoint gp, TimeStep tStep) override
void	computeConstitutiveMatrix_dPdF_At (FloatMatrix &answer, MatResponseMode rMode, GaussPoint gp, TimeStep tStep) override
void	computeInitialStressMatrix (FloatMatrix &answer, TimeStep *tStep) override
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	computeStiffnessMatrix_withIRulesAsSubcells (FloatMatrix &answer, MatResponseMode rMode, TimeStep *tStep)
void	giveInternalForcesVector (FloatArray &answer, TimeStep *tStep, int useUpdatedGpRecord=0) override
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	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	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)
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	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	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 (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 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)
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	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	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.
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 Member Functions inherited from oofem::ZZNodalRecoveryModelInterface
	ZZNodalRecoveryModelInterface (Element *element)
	Constructor.
virtual bool	ZZNodalRecoveryMI_computeNValProduct (FloatMatrix &answer, InternalStateType type, TimeStep *tStep)
virtual void	ZZNodalRecoveryMI_computeNNMatrix (FloatArray &answer, InternalStateType type)
Public Member Functions inherited from oofem::Interface
	Interface ()
	Constructor.
virtual	~Interface ()
Public Member Functions inherited from oofem::SPRNodalRecoveryModelInterface
	SPRNodalRecoveryModelInterface ()
	Constructor.
Public Member Functions inherited from oofem::NodalAveragingRecoveryModelInterface
	NodalAveragingRecoveryModelInterface ()
	Constructor.
Public Member Functions inherited from oofem::SpatialLocalizerInterface
	SpatialLocalizerInterface (Element *element)
virtual int	SpatialLocalizerI_containsPoint (const FloatArray &coords)
int	SpatialLocalizerI_BBoxContainsPoint (const FloatArray &coords)
virtual void	SpatialLocalizerI_giveBBox (FloatArray &bb0, FloatArray &bb1)
virtual double	SpatialLocalizerI_giveClosestPoint (FloatArray &lcoords, FloatArray &closest, const FloatArray &gcoords)
Public Member Functions inherited from oofem::HuertaErrorEstimatorInterface
	HuertaErrorEstimatorInterface ()
	Constructor.

Protected Attributes
IntArray	microMasterNodes
	Array containing the node mapping from microscale (which microMasterNodes corresponds to which macroNode).
IntArray	microBoundaryNodes
IntArray	microDOFs
bool	firstCall
MicroMaterial *	microMaterial
Domain *	microDomain
EngngModel *	microEngngModel
int	iteration
	Information of iteration number.
IntArray	microBoundaryDofManager
	Stores node number on the boundary in the triplets.
FloatMatrix	stiffMatrix
int	stiffMatrxFileNoneReadingWriting
	Process with external file for the storage of stiffness matrix 0-None, 1-read, 2-write.
FloatArray	internalMacroForcesVector
	Array containg the force vector from nodes (if condensation is skipped, use this vector).
TimeStep *	lastStiffMatrixTimeStep
	Last time step when stiffness matrix was assembled.
Protected Attributes inherited from oofem::LSpace
bool	reducedShearIntegration
Protected Attributes inherited from oofem::Structural3DElement
bool	matRotation
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.

Static Private Attributes
static ParamKey	IPK_MacroLSpace_microMasterNodes
static ParamKey	IPK_MacroLSpace_microBoundaryNodes
static ParamKey	IPK_MacroLSpace_stiffMatrixFileName

Additional Inherited Members
Public Types inherited from oofem::HuertaErrorEstimatorInterface
enum	SetupMode { CountMode = 0 , NodeMode = 1 , ElemMode = 2 , BCMode = 3 }
	Mode for problem setup. More...
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::LSpace
int	giveNumberOfIPForMassMtrxIntegration () override
void	computeBmatrixAt (GaussPoint *gp, FloatMatrix &answer, int li, int ui) override
int	computeLoadLSToLRotationMatrix (FloatMatrix &answer, int iSurf, GaussPoint *gp) override
Protected Member Functions inherited from oofem::Structural3DElement
void	computeBmatrixAt (GaussPoint *gp, FloatMatrix &answer, int lowerIndx=1, int upperIndx=ALL_STRAINS) override
void	computeBHmatrixAt (GaussPoint *gp, FloatMatrix &answer) override
void	computeGaussPoints () override
void	giveEdgeDofMapping (IntArray &answer, int iEdge) const override
double	computeEdgeVolumeAround (GaussPoint *gp, int iEdge) override
int	computeLoadLEToLRotationMatrix (FloatMatrix &answer, int iEdge, GaussPoint *gp) override
int	testElementExtension (ElementExtension ext) override
void	computeSurfaceNMatrixAt (FloatMatrix &answer, int iSurf, GaussPoint *gp)
void	giveSurfaceDofMapping (IntArray &answer, int) const override
double	computeSurfaceVolumeAround (GaussPoint *gp, int) override
int	computeLoadLSToLRotationMatrix (FloatMatrix &answer, int, GaussPoint *gp) override
Protected Member Functions inherited from oofem::NLStructuralElement
int	checkConsistency () override
virtual void	computePointLoadVectorAt (FloatArray &answer, Load load, TimeStep tStep, ValueModeType mode, bool global=true)
virtual int	computeLoadGToLRotationMtrx (FloatMatrix &answer)
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)
Protected Member Functions inherited from oofem::HuertaErrorEstimatorInterface
void	setupRefinedElementProblem1D (Element element, RefinedElement refinedElement, int level, int nodeId, IntArray &localNodeIdArray, IntArray &globalNodeIdArray, HuertaErrorEstimatorInterface ::SetupMode mode, TimeStep tStep, int nodes, FloatArray corner, FloatArray &midNode, int &localNodeId, int &localElemId, int &localBcId, IntArray &controlNode, IntArray &controlDof, HuertaErrorEstimator ::AnalysisMode aMode, const char *edgetype)
void	setupRefinedElementProblem2D (Element element, RefinedElement refinedElement, int level, int nodeId, IntArray &localNodeIdArray, IntArray &globalNodeIdArray, HuertaErrorEstimatorInterface ::SetupMode mode, TimeStep tStep, int nodes, FloatArray corner, FloatArray midSide, FloatArray &midNode, int &localNodeId, int &localElemId, int &localBcId, IntArray &controlNode, IntArray &controlDof, HuertaErrorEstimator ::AnalysisMode aMode, const char quadtype)
void	setupRefinedElementProblem3D (Element element, RefinedElement refinedElement, int level, int nodeId, IntArray &localNodeIdArray, IntArray &globalNodeIdArray, HuertaErrorEstimatorInterface ::SetupMode mode, TimeStep tStep, int nodes, FloatArray corner, FloatArray midSide, FloatArray midFace, FloatArray &midNode, int &localNodeId, int &localElemId, int &localBcId, int hexaSideNode[1][3], int hexaFaceNode[1][3], IntArray &controlNode, IntArray &controlDof, HuertaErrorEstimator ::AnalysisMode aMode, const char *hexatype)
Static Protected Attributes inherited from oofem::LSpace
static FEI3dHexaLin	interpolation
static ParamKey	IPK_LSpace_reducedShearIntegration
	[optional] Reduced shear integration flag.
Static Protected Attributes inherited from oofem::Structural3DElement
static ParamKey	IPK_Structural3DElement_materialCoordinateSystem
	[optional] Material coordinate system (local) for the element.
Static Protected Attributes inherited from oofem::NLStructuralElement
static ParamKey	IPK_NLStructuralElement_nlgeoflag

Detailed Description

This class implements a macroelement. It is derived from eight-node brick element. The stiffness matrix is computed from underlying RVE and is condensed to 24 DoFs to corner nodes.

Author: Vit Smilauer

Definition at line 60 of file macrolspace.h.

Constructor & Destructor Documentation

◆ MacroLSpace()

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

Definition at line 55 of file macrolspace.C.

References firstCall, iteration, lastStiffMatrixTimeStep, oofem::LSpace::LSpace(), microDomain, microEngngModel, and microMaterial.

◆ ~MacroLSpace()

oofem::MacroLSpace::~MacroLSpace ( )

virtual

Definition at line 65 of file macrolspace.C.

Member Function Documentation

◆ changeMicroBoundaryConditions()

void oofem::MacroLSpace::changeMicroBoundaryConditions ( TimeStep * tStep )

virtual

Related to setting the boundary conditions of micro problem.

Definition at line 161 of file macrolspace.C.

References _IFT_BoundaryCondition_PrescribedValue, _IFT_ConstantFunction_f, _IFT_GeneralBoundaryCondition_timeFunct, oofem::FloatArray::at(), oofem::classFactory, evalInterpolation(), oofem::Element::giveNode(), oofem::Element::giveNumberOfNodes(), microBoundaryDofManager, microBoundaryNodes, microDomain, microMaterial, OOFEM_ERROR, oofem::DynamicInputRecord::setField(), and oofem::DynamicInputRecord::setRecordKeywordField().

Referenced by giveInternalForcesVector().

◆ computeField()

void oofem::MacroLSpace::computeField	(	ValueModeType	mode,
		TimeStep *	tStep,
		const FloatArray &	lcoords,
		FloatArray &	answer )

inlineoverridevirtual

Computes the unknown vector interpolated at the specified local coordinates. Used for exporting data and mapping fields.

See also: giveElementDofIDMask The unknown vector should match the element field as specified by the element dof IDs.

Parameters

mode	Mode (total, increment, etc) of the output
tStep	Time step to evaluate at
lcoords	Local coordinates to evaluate at
answer	Results

Reimplemented from oofem::Element.

Definition at line 78 of file macrolspace.h.

References OOFEM_ERROR.

◆ computeStiffnessMatrix()

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

overridevirtual

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.

Definition at line 108 of file macrolspace.C.

References firstCall, oofem::FEMComponent::giveDomain(), oofem::TimeStep::giveIntrinsicTime(), oofem::Element::giveMaterial(), oofem::FEMComponent::giveNumber(), oofem::TimeStep::giveTargetTime(), oofem::Element::isActivated(), lastStiffMatrixTimeStep, microBoundaryNodes, microDomain, microEngngModel, microMasterNodes, microMaterial, oofem::microScale, OOFEM_ERROR, OOFEM_LOG_INFO, and stiffMatrix.

◆ evalInterpolation()

void oofem::MacroLSpace::evalInterpolation	(	FloatArray &	answer,
		const std::vector< FloatArray > &	coords,
		const FloatArray &	gcoords )

virtual

Evaluates shape function at a given pointnodal representation of real internal forces obtained from microProblem.

Parameters

answer	Array of shape function values at given node.
coords	Coordinates of nodes defining the interpolation geometry.
gcoords	Global coordinates of point of interest.

Definition at line 285 of file macrolspace.C.

References oofem::LSpace::interpolation.

Referenced by changeMicroBoundaryConditions(), and giveInternalForcesVector().

◆ giveClassName()

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

inlineoverridevirtual

Returns: Class name of the receiver.

Reimplemented from oofem::LSpace.

Definition at line 71 of file macrolspace.h.

◆ giveInputRecordName()

const char * oofem::MacroLSpace::giveInputRecordName ( ) const

inlineoverridevirtual

Returns: Input record name of the receiver.

Reimplemented from oofem::LSpace.

Definition at line 70 of file macrolspace.h.

References _IFT_MacroLSpace_Name.

◆ giveInternalForcesVector()

void oofem::MacroLSpace::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::StructuralElement.

Definition at line 231 of file macrolspace.C.

References oofem::FloatArray::at(), changeMicroBoundaryConditions(), oofem::StructuralEngngModel::computeReaction(), evalInterpolation(), oofem::DofManager::giveCoordinates(), oofem::TimeStep::giveIntrinsicTime(), oofem::FEMComponent::giveNumber(), oofem::DofManager::giveNumberOfDofs(), oofem::TimeStep::giveTargetTime(), internalMacroForcesVector, iteration, microBoundaryDofManager, microDomain, microEngngModel, microMaterial, OOFEM_LOG_INFO, oofem::FloatArray::resize(), and oofem::FloatArray::zero().

◆ initializeFrom()

void oofem::MacroLSpace::initializeFrom	(	InputRecord &	ir,
		int	priority )

overridevirtual

Reimplemented from oofem::LSpace.

Definition at line 68 of file macrolspace.C.

References _IFT_MacroLspace_stiffMatrxFileName, oofem::FEMComponent::domain, oofem::InputRecord::hasField(), IPK_MacroLSpace_microBoundaryNodes, IPK_MacroLSpace_microMasterNodes, IR_GIVE_OPTIONAL_FIELD, microBoundaryNodes, microMasterNodes, oofem::FEMComponent::number, OOFEM_ERROR, PM_UPDATE_PARAMETER, and stiffMatrxFileNoneReadingWriting.

◆ postInitialize()

void oofem::MacroLSpace::postInitialize ( )

overridevirtual

Performs post initialization steps.

Reimplemented from oofem::Element.

Definition at line 92 of file macrolspace.C.

References oofem::ComponentInputException::ctElement, oofem::FEMComponent::domain, IPK_MacroLSpace_microBoundaryNodes, IPK_MacroLSpace_microMasterNodes, microBoundaryDofManager, microBoundaryNodes, microMasterNodes, oofem::FEMComponent::number, and PM_ELEMENT_ERROR_IFNOTSET.

◆ updateYourself()

void oofem::MacroLSpace::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.

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

Reimplemented from oofem::Element.

Definition at line 297 of file macrolspace.C.

References oofem::Element::integrationRulesArray, iteration, microEngngModel, and OOFEM_LOG_INFO.

Member Data Documentation

◆ firstCall

bool oofem::MacroLSpace::firstCall

protected

Definition at line 101 of file macrolspace.h.

Referenced by computeStiffnessMatrix(), and MacroLSpace().

◆ internalMacroForcesVector

FloatArray oofem::MacroLSpace::internalMacroForcesVector

protected

Array containg the force vector from nodes (if condensation is skipped, use this vector).

Definition at line 113 of file macrolspace.h.

Referenced by giveInternalForcesVector().

◆ IPK_MacroLSpace_microBoundaryNodes

ParamKey oofem::MacroLSpace::IPK_MacroLSpace_microBoundaryNodes

staticprivate

Definition at line 64 of file macrolspace.h.

Referenced by initializeFrom(), and postInitialize().

◆ IPK_MacroLSpace_microMasterNodes

ParamKey oofem::MacroLSpace::IPK_MacroLSpace_microMasterNodes

staticprivate

Definition at line 63 of file macrolspace.h.

Referenced by initializeFrom(), and postInitialize().

◆ IPK_MacroLSpace_stiffMatrixFileName

ParamKey oofem::MacroLSpace::IPK_MacroLSpace_stiffMatrixFileName

staticprivate

Definition at line 65 of file macrolspace.h.

◆ iteration

int oofem::MacroLSpace::iteration

protected

Information of iteration number.

Definition at line 106 of file macrolspace.h.

Referenced by giveInternalForcesVector(), MacroLSpace(), and updateYourself().

◆ lastStiffMatrixTimeStep

TimeStep* oofem::MacroLSpace::lastStiffMatrixTimeStep

protected

Last time step when stiffness matrix was assembled.

Definition at line 115 of file macrolspace.h.

Referenced by computeStiffnessMatrix(), and MacroLSpace().

◆ microBoundaryDofManager

IntArray oofem::MacroLSpace::microBoundaryDofManager

protected

Stores node number on the boundary in the triplets.

Definition at line 108 of file macrolspace.h.

Referenced by changeMicroBoundaryConditions(), giveInternalForcesVector(), and postInitialize().

◆ microBoundaryNodes

IntArray oofem::MacroLSpace::microBoundaryNodes

protected

Definition at line 99 of file macrolspace.h.

Referenced by changeMicroBoundaryConditions(), computeStiffnessMatrix(), initializeFrom(), and postInitialize().

◆ microDOFs

IntArray oofem::MacroLSpace::microDOFs

protected

Definition at line 100 of file macrolspace.h.

◆ microDomain

Domain* oofem::MacroLSpace::microDomain

protected

Definition at line 103 of file macrolspace.h.

Referenced by changeMicroBoundaryConditions(), computeStiffnessMatrix(), giveInternalForcesVector(), and MacroLSpace().

◆ microEngngModel

EngngModel* oofem::MacroLSpace::microEngngModel

protected

Definition at line 104 of file macrolspace.h.

Referenced by computeStiffnessMatrix(), giveInternalForcesVector(), MacroLSpace(), and updateYourself().

◆ microMasterNodes

IntArray oofem::MacroLSpace::microMasterNodes

protected

Array containing the node mapping from microscale (which microMasterNodes corresponds to which macroNode).

Definition at line 98 of file macrolspace.h.

Referenced by computeStiffnessMatrix(), initializeFrom(), and postInitialize().

◆ microMaterial

MicroMaterial* oofem::MacroLSpace::microMaterial

protected

Definition at line 102 of file macrolspace.h.

Referenced by changeMicroBoundaryConditions(), computeStiffnessMatrix(), giveInternalForcesVector(), and MacroLSpace().

◆ stiffMatrix

FloatMatrix oofem::MacroLSpace::stiffMatrix

protected

Definition at line 109 of file macrolspace.h.

Referenced by computeStiffnessMatrix().

◆ stiffMatrxFileNoneReadingWriting

int oofem::MacroLSpace::stiffMatrxFileNoneReadingWriting

protected

Process with external file for the storage of stiffness matrix 0-None, 1-read, 2-write.

Definition at line 111 of file macrolspace.h.

Referenced by initializeFrom().

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

Public Member Functions

Protected Attributes

Static Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ MacroLSpace()

◆ ~MacroLSpace()

Member Function Documentation

◆ changeMicroBoundaryConditions()

◆ computeField()

◆ computeStiffnessMatrix()

◆ evalInterpolation()

◆ giveClassName()

◆ giveInputRecordName()

◆ giveInternalForcesVector()

◆ initializeFrom()

◆ postInitialize()

◆ updateYourself()

Member Data Documentation

◆ firstCall

◆ internalMacroForcesVector

◆ IPK_MacroLSpace_microBoundaryNodes

◆ IPK_MacroLSpace_microMasterNodes

◆ IPK_MacroLSpace_stiffMatrixFileName

◆ iteration

◆ lastStiffMatrixTimeStep

◆ microBoundaryDofManager

◆ microBoundaryNodes

◆ microDOFs

◆ microDomain

◆ microEngngModel

◆ microMasterNodes

◆ microMaterial

◆ stiffMatrix

◆ stiffMatrxFileNoneReadingWriting