This abstract class represent a general base element class for transport problems. More...

#include <transportelement.h>

Inheritance diagram for oofem::TransportElement:

Collaboration diagram for oofem::TransportElement:

Public Types
enum	ElementMode { HeatTransferEM, HeatMass1TransferEM, Mass1TransferEM }

Public Member Functions
	TransportElement (int n, Domain *d, ElementMode em=HeatTransferEM)

virtual	~TransportElement ()

virtual void	giveCharacteristicMatrix (FloatMatrix &answer, CharType type, TimeStep *tStep)
	Computes characteristic matrix of receiver of requested type in given time step. More...

virtual void	giveCharacteristicVector (FloatArray &answer, CharType type, ValueModeType mode, TimeStep *tStep)
	Computes characteristic vector of receiver of requested type in given time step. More...

virtual void	giveDofManDofIDMask (int inode, IntArray &answer) const
	Returns dofmanager dof mask for node. More...

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

virtual void	computeExternalForcesVector (FloatArray &answer, TimeStep *tStep, ValueModeType mode)

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

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

virtual void	computeLoadVector (FloatArray &answer, BodyLoad load, CharType type, ValueModeType mode, TimeStep tStep)
	Computes the contribution of the given body load (volumetric). More...

virtual void	computeBoundarySurfaceLoadVector (FloatArray &answer, BoundaryLoad load, int boundary, CharType type, ValueModeType mode, TimeStep tStep, bool global=true)
	Computes the contribution of the given load at the given boundary surface in global coordinate system. More...

virtual void	computeTangentFromSurfaceLoad (FloatMatrix &answer, SurfaceLoad load, int boundary, MatResponseMode rmode, TimeStep tStep)
	Computes the tangent contribution of the given load at the given boundary. More...

virtual void	computeTangentFromEdgeLoad (FloatMatrix &answer, EdgeLoad load, int boundary, MatResponseMode rmode, TimeStep tStep)
	Computes the tangent contribution of the given load at the given boundary. More...

virtual void	computeBoundaryEdgeLoadVector (FloatArray &answer, BoundaryLoad load, int edge, CharType type, ValueModeType mode, TimeStep tStep, bool global=true)
	Computes the contribution of the given load at the given boundary edge. More...

virtual void	computeField (ValueModeType mode, TimeStep *tStep, const FloatArray &lcoords, FloatArray &answer)
	Computes the unknown vector interpolated at the specified local coordinates. More...

TransportCrossSection *	giveTransportCrossSection ()

virtual Material *	giveMaterial ()

virtual double	giveThicknessAt (const FloatArray &gcoords)
	Gives the thickness at some global coordinate. More...

virtual void	computeCapacityMatrix (FloatMatrix &answer, TimeStep *tStep)
	Computes the capacity matrix of the receiver. More...

virtual void	computeConductivityMatrix (FloatMatrix &answer, MatResponseMode rMode, TimeStep *tStep)
	Computes the conductivity matrix of the receiver. More...

virtual void	computeBCVectorAt (FloatArray &answer, TimeStep *tStep, ValueModeType mode)
	Computes the RHS contribution to balance equation(s) due to boundary conditions. More...

virtual void	computeBCMtrxAt (FloatMatrix &answer, TimeStep *tStep, ValueModeType mode)
	Computes the LHS contribution to balance equation(s) due to boundary conditions. More...

virtual void	computeInternalSourceRhsVectorAt (FloatArray &answer, TimeStep *tStep, ValueModeType mode)
	Computes the contribution to balance equation(s) due to internal sources. More...

virtual void	computeIntSourceLHSMatrix (FloatMatrix &answer, TimeStep *tStep)
	Computes the LHS contribution to balance equation(s) due to material internal source. More...

virtual void	computeIntSourceLHSSubMatrix (FloatMatrix &answer, MatResponseMode rmode, int iri, TimeStep *tStep)
	Computes the part of internal source LHS contribution corresponding to unknown identified by rmode parameter. More...

virtual void	computeFlow (FloatArray &answer, GaussPoint gp, TimeStep tStep)
	Computes a flow vector in an integration point. More...

virtual void	updateInternalState (TimeStep *tStep)
	Updates element state after equilibrium in time step has been reached. More...

virtual int	checkConsistency ()
	Performs consistency check. More...

virtual int	EIPrimaryFieldI_evaluateFieldVectorAt (FloatArray &answer, PrimaryField &pf, const FloatArray &coords, IntArray &dofId, ValueModeType mode, TimeStep *tStep)
	Evaluates the value of field at given point of interest (should be located inside receiver's volume) using element interpolation. More...

int	giveInternalStateAtNode (FloatArray &answer, InternalStateType type, InternalStateMode mode, int node, TimeStep *tStep)
	Returns internal state variable (like stress,strain) at node of element in Reduced form, the way how is obtained is dependent on InternalValueType. More...

virtual void	computeConstitutiveMatrixAt (FloatMatrix &answer, MatResponseMode rMode, GaussPoint gp, TimeStep tStep)
	Computes constitutive matrix of receiver. More...

Public Member Functions inherited from oofem::Element
	Element (int n, Domain *aDomain)
	Constructor. More...

	Element (const Element &src)=delete

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

virtual	~Element ()
	Virtual destructor. More...

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)
	Returns internal state variable (like stress,strain) at side of element in Reduced form If side is possessing DOFs, otherwise recover techniques will not work due to absence of side-shape functions. More...

virtual void	showSparseMtrxStructure (CharType mtrx, oofegGraphicContext &gc, TimeStep *tStep)
	Shows sparse structure. More...

virtual void	showExtendedSparseMtrxStructure (CharType mtrx, oofegGraphicContext &gc, TimeStep *tStep)
	Shows extended sparse structure (for example, due to nonlocal interactions for tangent stiffness) More...

int	giveLabel () const

int	giveGlobalNumber () const

void	setGlobalNumber (int num)
	Sets receiver globally unique number. More...

elementParallelMode	giveParallelMode () const
	Return elementParallelMode of receiver. More...

void	setParallelMode (elementParallelMode _mode)
	Sets parallel mode of element. More...

virtual elementParallelMode	giveKnotSpanParallelMode (int) const
	Returns the parallel mode for particular knot span of the receiver. More...

int	packUnknowns (DataStream &buff, TimeStep *tStep)
	Pack all necessary data of element (according to its parallel_mode) integration points into given communication buffer. More...

int	unpackAndUpdateUnknowns (DataStream &buff, TimeStep *tStep)
	Unpack and updates all necessary data of element (according to its parallel_mode) integration points into given communication buffer. More...

int	estimatePackSize (DataStream &buff)
	Estimates the necessary pack size to hold all packed data of receiver. More...

const IntArray *	givePartitionList () const
	Returns partition list of receiver. More...

void	setPartitionList (IntArray &pl)
	Sets partition list of receiver. More...

virtual double	predictRelativeComputationalCost ()
	Returns the weight representing relative computational cost of receiver The reference element is triangular plane stress element with linear approximation, single integration point and linear isotropic material. More...

virtual double	giveRelativeSelfComputationalCost ()
	Returns the weight representing relative computational cost of receiver The reference element is triangular plane stress element. More...

virtual double	predictRelativeRedistributionCost ()
	Returns the relative redistribution cost of the receiver. More...

IntArray *	giveBodyLoadArray ()
	Returns array containing load numbers of loads acting on element. More...

IntArray *	giveBoundaryLoadArray ()
	Returns array containing load numbers of boundary loads acting on element. More...

virtual IRResultType	initializeFrom (InputRecord *ir)
	Initializes receiver according to object description stored in input record. More...

virtual void	giveInputRecord (DynamicInputRecord &input)
	Setups the input record string of receiver. More...

virtual contextIOResultType	saveContext (DataStream &stream, ContextMode mode, void *obj=NULL)
	Stores receiver state to output stream. More...

virtual contextIOResultType	restoreContext (DataStream &stream, ContextMode mode, void *obj=NULL)
	Restores the receiver state previously written in stream. More...

virtual void	printOutputAt (FILE file, TimeStep tStep)
	Prints output of receiver to stream, for given time step. More...

virtual const char *	giveClassName () const

void	giveLocationArray (IntArray &locationArray, const UnknownNumberingScheme &s, IntArray *dofIds=NULL) const
	Returns the location array (array of code numbers) of receiver for given numbering scheme. More...

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)
	Returns the location array for the boundary of the element. More...

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
	Returns i-th internal element dof manager of the receiver. More...

virtual double	giveCharacteristicValue (CharType type, TimeStep *tStep)
	Computes characteristic value of receiver of requested type in given time step. More...

const IntArray &	giveBodyLoadList () const
	Returns receiver list of bodyloads. More...

const IntArray &	giveBoundaryLoadList () const
	Returns receiver list of boundary loads. More...

void	computeVectorOf (ValueModeType u, TimeStep *tStep, FloatArray &answer)
	Returns local vector of unknowns. More...

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)
	Boundary version of computeVectorOf. More...

void	computeVectorOf (PrimaryField &field, const IntArray &dofIDMask, ValueModeType u, TimeStep *tStep, FloatArray &answer, bool padding=false)
	Returns local vector of unknowns. More...

void	computeVectorOfPrescribed (ValueModeType u, TimeStep *tStep, FloatArray &answer)
	Returns local vector of prescribed unknowns. More...

void	computeVectorOfPrescribed (const IntArray &dofIDMask, ValueModeType type, TimeStep *tStep, FloatArray &answer)
	Returns local vector of prescribed unknowns. More...

virtual int	computeNumberOfDofs ()
	Computes or simply returns total number of element's local DOFs. More...

virtual int	computeNumberOfGlobalDofs ()
	Computes the total number of element's global dofs. More...

int	computeNumberOfPrimaryMasterDofs ()
	Computes the total number of element's primary master DOFs. More...

virtual bool	computeGtoLRotationMatrix (FloatMatrix &answer)
	Returns transformation matrix from global c.s. More...

virtual bool	giveRotationMatrix (FloatMatrix &answer)
	Transformation matrices updates rotation matrix between element-local and primary DOFs, taking into account nodal c.s. More...

virtual bool	computeDofTransformationMatrix (FloatMatrix &answer, const IntArray &nodes, bool includeInternal)
	Returns transformation matrix for DOFs from global coordinate system to local coordinate system in nodes. More...

virtual void	giveInternalDofManDofIDMask (int inode, IntArray &answer) const
	Returns internal dofmanager dof mask for node. More...

virtual void	giveElementDofIDMask (IntArray &answer) const
	Returns element dof mask for node. More...

virtual double	computeVolumeAround (GaussPoint *gp)
	Returns volume related to given integration point. More...

virtual double	computeVolumeAreaOrLength ()
	Computes the volume, area or length of the element depending on its spatial dimension. More...

double	computeMeanSize ()
	Computes the size of the element defined as its length. More...

virtual double	computeVolume ()
	Computes the volume. More...

virtual double	computeArea ()
	Computes the area (zero for all but 2d geometries). More...

virtual double	computeLength ()
	Computes the length (zero for all but 1D geometries) More...

virtual void	giveBoundaryEdgeNodes (IntArray &bNodes, int boundary)
	Returns list of receiver boundary nodes for given edge. More...

virtual void	giveBoundarySurfaceNodes (IntArray &bNodes, int boundary)
	Returns list of receiver boundary nodes for given surface. More...

virtual IntegrationRule *	giveBoundaryEdgeIntegrationRule (int order, int boundary)
	Returns boundary edge integration rule. More...

virtual IntegrationRule *	giveBoundarySurfaceIntegrationRule (int order, int boundary)
	Returns boundary surface integration rule. More...

int	giveDofManagerNumber (int i) const
	Translates local to global indices for dof managers. More...

const IntArray &	giveDofManArray () const

void	addDofManager (DofManager *dMan)

DofManager *	giveDofManager (int i) const

Node *	giveNode (int i) const
	Returns reference to the i-th node of element. More...

virtual ElementSide *	giveSide (int i) const
	Returns reference to the i-th side of element. More...

virtual FEInterpolation *	giveInterpolation () const

virtual FEInterpolation *	giveInterpolation (DofIDItem id) const
	Returns the interpolation for the specific dof id. More...

int	giveMaterialNumber () const

CrossSection *	giveCrossSection ()

void	setMaterial (int matIndx)
	Sets the material of receiver. More...

virtual void	setCrossSection (int csIndx)
	Sets the cross section model of receiver. More...

virtual int	giveNumberOfDofManagers () const

virtual int	giveNumberOfNodes () const
	Returns number of nodes of receiver. More...

void	setDofManagers (const IntArray &dmans)
	Sets receiver dofManagers. More...

void	setBodyLoads (const IntArray &bodyLoads)
	Sets receiver bodyLoadArray. More...

void	setIntegrationRules (std::vector< std::unique_ptr< IntegrationRule > > irlist)
	Sets integration rules. More...

virtual integrationDomain	giveIntegrationDomain () const
	Returns integration domain for receiver, used to initialize integration point over receiver volume. More...

virtual MaterialMode	giveMaterialMode ()
	Returns material mode for receiver integration points. More...

virtual int	giveIntegrationRuleLocalCodeNumbers (IntArray &answer, IntegrationRule &ie)
	Assembles the code numbers of given integration element (sub-patch) This is done by obtaining list of nonzero shape functions and by collecting the code numbers of nodes corresponding to these shape functions. More...

int	giveRegionNumber ()

virtual void	postInitialize ()
	Performs post initialization steps. More...

virtual void	updateYourself (TimeStep *tStep)
	Updates element state after equilibrium in time step has been reached. More...

virtual void	initializeYourself (TimeStep *timeStepWhenICApply)
	Initialization according to state given by initial conditions. More...

virtual bool	isActivated (TimeStep *tStep)

virtual bool	isCast (TimeStep *tStep)

virtual void	initForNewStep ()
	Initializes receivers state to new time step. More...

virtual Element_Geometry_Type	giveGeometryType () const
	Returns the element geometry type. More...

virtual int	giveSpatialDimension ()
	Returns the element spatial dimension (1, 2, or 3). More...

virtual int	giveNumberOfBoundarySides ()

virtual int	giveDefaultIntegrationRule () const
	Returns id of default integration rule. More...

virtual IntegrationRule *	giveDefaultIntegrationRulePtr ()
	Access method for default integration rule. More...

int	giveNumberOfIntegrationRules ()

virtual IntegrationRule *	giveIntegrationRule (int i)

std::vector< std::unique_ptr< IntegrationRule > > &	giveIntegrationRulesArray ()

virtual int	testElementExtension (ElementExtension ext)
	Tests if the element implements required extension. More...

virtual int	giveIPValue (FloatArray &answer, GaussPoint gp, InternalStateType type, TimeStep tStep)
	Returns the integration point corresponding value in full form. More...

int	giveGlobalIPValue (FloatArray &answer, GaussPoint gp, InternalStateType type, TimeStep tStep)

virtual double	giveLengthInDir (const FloatArray &normalToCrackPlane)
	Default implementation returns length of element projection into specified direction. More...

virtual double	giveCharacteristicLength (const FloatArray &normalToCrackPlane)
	Returns the size of element in the given direction, in some cases adjusted (e.g. More...

double	giveCharacteristicLengthForPlaneElements (const FloatArray &normalToCrackPlane)
	Returns the size of element in the given direction if the direction is in the XY plane, otherwise gives the mean size defined as the square root of the element area. More...

double	giveCharacteristicLengthForAxisymmElements (const FloatArray &normalToCrackPlane)
	Returns the size of an axisymmetric element in the given direction if the direction is in the XY plane, otherwise gives the mean distance vrom the symmetry axis multiplied by pi. More...

virtual double	giveCharacteristicSize (GaussPoint *gp, FloatArray &normalToCrackPlane, ElementCharSizeMethod method)
	Returns characteristic element size for a given integration point and given direction. More...

virtual double	giveParentElSize () const
	Returns the size (length, area or volume depending on element type) of the parent element. More...

virtual void	updateBeforeNonlocalAverage (TimeStep *tStep)
	Updates internal element state (in all integration points of receiver) before nonlocal averaging takes place. More...

virtual int	computeGlobalCoordinates (FloatArray &answer, const FloatArray &lcoords)
	Computes the global coordinates from given element's local coordinates. More...

virtual bool	computeLocalCoordinates (FloatArray &answer, const FloatArray &gcoords)
	Computes the element local coordinates from given global coordinates. More...

virtual int	giveLocalCoordinateSystem (FloatMatrix &answer)
	Returns local coordinate system of receiver Required by material models with ortho- and anisotrophy. More...

virtual void	computeMidPlaneNormal (FloatArray &answer, const GaussPoint *gp)
	Computes mid-plane normal of receiver at integration point. More...

virtual int	adaptiveMap (Domain oldd, TimeStep tStep)
	Initializes the internal state variables stored in all IPs according to state in given domain. More...

virtual int	mapStateVariables (Domain &iOldDom, const TimeStep &iTStep)
	Maps the internal state variables stored in all IPs from the old domain to the new domain. More...

virtual int	adaptiveUpdate (TimeStep *tStep)
	Updates the internal state variables stored in all IPs according to already mapped state. More...

virtual int	adaptiveFinish (TimeStep *tStep)
	Finishes the mapping for given time step. More...

virtual void	updateLocalNumbering (EntityRenumberingFunctor &f)
	Local renumbering support. More...

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. More...

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. More...

Public Member Functions inherited from oofem::FEMComponent
	FEMComponent (int n, Domain *d)
	Regular constructor, creates component with given number and belonging to given domain. More...

virtual	~FEMComponent ()
	Virtual destructor. More...

virtual const char *	giveInputRecordName () const =0

Domain *	giveDomain () const

virtual void	setDomain (Domain *d)
	Sets associated Domain. More...

int	giveNumber () const

void	setNumber (int num)
	Sets number of receiver. More...

virtual void	printYourself ()
	Prints receiver state on stdout. Useful for debugging. More...

virtual Interface *	giveInterface (InterfaceType t)
	Interface requesting service. More...

std::string	errorInfo (const char *func) const
	Returns string for prepending output (used by error reporting macros). More...

Public Member Functions inherited from oofem::EIPrimaryFieldInterface
	EIPrimaryFieldInterface ()

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

virtual	~Interface ()

virtual const char *	giveClassName () const =0

Protected Member Functions
virtual void	computeCapacitySubMatrix (FloatMatrix &answer, MatResponseMode rmode, int iri, TimeStep *tStep)
	Computes the matrix $P=\int_\Omega N^{\mathrm{T}} N\;\mathrm{d}\Omega$ . More...

virtual void	computeConductivitySubMatrix (FloatMatrix &answer, int iri, MatResponseMode rmode, TimeStep *tStep)
	Computes the matrix $C=\int_\Omega \sum_i \left( \frac{\partial N}{\partial x_i}\right)^{\mathrm{T}} \frac{\partial N}{\partial x_i} \;\mathrm{d}\Omega$ . More...

void	computeBCSubVectorAt (FloatArray &answer, TimeStep *tStep, ValueModeType mode, int indx)
	Computes the part of RHS due to applied BCs. More...

void	computeBCSubMtrxAt (FloatMatrix &answer, TimeStep *tStep, ValueModeType mode, int indx)
	Computes the part of LHS due to applied BCs. More...

void	computeBodyBCSubVectorAt (FloatArray &answer, Load load, TimeStep tStep, ValueModeType mode, int indx)

void	computeEdgeBCSubVectorAt (FloatArray &answer, Load load, int iEdge, TimeStep tStep, ValueModeType mode, int indx)
	Computes the part of RHS due to applied BCs on particular edge. More...

void	computeSurfaceBCSubVectorAt (FloatArray &answer, Load load, int iSurf, TimeStep tStep, ValueModeType mode, int indx)
	Computes the part of RHS due to applied BCs on particular surface. More...

virtual void	computeNAt (FloatArray &answer, const FloatArray &lcoord)
	Computes the basis functions. More...

virtual void	computeNmatrixAt (FloatMatrix &answer, const FloatArray &lcoords)
	Computes the interpolation matrix corresponding to all unknowns. More...

virtual void	computeBmatrixAt (FloatMatrix &answer, const FloatArray &lcoords)

virtual void	computeGradientMatrixAt (FloatMatrix &answer, const FloatArray &lcoords)
	Computes the gradient matrix corresponding to one unknown. More...

virtual void	computeInternalSourceRhsSubVectorAt (FloatArray &answer, TimeStep *, ValueModeType mode, int indx)
	Computes the contribution to balance equation(s) due to internal sources. More...

virtual void	computeEgdeNAt (FloatArray &answer, int iEdge, const FloatArray &lcoord)
	Computes the basis functions at the edge for one unknown. More...

virtual void	giveEdgeDofMapping (IntArray &mask, int iEdge)
	Gives the node indexes for given edge. More...

virtual double	computeEdgeVolumeAround (GaussPoint *gp, int iEdge)=0
	Computes the length around a integration point on a edge. More...

virtual IntegrationRule *	GetSurfaceIntegrationRule (int approxOrder)

virtual void	computeSurfaceNAt (FloatArray &answer, int iSurf, const FloatArray &lcoord)

virtual double	computeSurfaceVolumeAround (GaussPoint *gp, int iSurf)

virtual void	giveSurfaceDofMapping (IntArray &mask, int iSurf)

virtual int	giveApproxOrder (int unknownIndx)

void	assembleLocalContribution (FloatMatrix &answer, FloatMatrix &src, int ndofs, int rdof, int cdof)
	Assembles the given source matrix of size (ndofs, ndofs) into target matrix answer. More...

void	assembleLocalContribution (FloatArray &answer, FloatArray &src, int ndofs, int rdof)
	Assembles the given source vector of size (ndofs) into target answer. More...

Protected Member Functions inherited from oofem::Element
virtual void	computeGaussPoints ()
	Initializes the array of integration rules member variable. More...

Protected Attributes
ElementMode	emode

Protected Attributes inherited from oofem::Element
int	numberOfDofMans
	Number of dofmanagers. More...

IntArray	dofManArray
	Array containing dofmanager numbers. More...

int	material
	Number of associated material. More...

int	crossSection
	Number of associated cross section. More...

IntArray	bodyLoadArray
	Array containing indexes of loads (body loads and boundary loads are kept separately), that apply on receiver. More...

IntArray	boundaryLoadArray

std::vector< std::unique_ptr< IntegrationRule > >	integrationRulesArray
	List of integration rules of receiver (each integration rule contains associated integration points also). More...

FloatMatrix	elemLocalCS
	Transformation material matrix, used in orthotropic and anisotropic materials, global->local transformation. More...

int	activityTimeFunction
	Element activity time function. If defined, nonzero value indicates active receiver, zero value inactive element. More...

int	globalNumber
	In parallel mode, globalNumber contains globally unique DoFManager number. More...

int	numberOfGaussPoints
	Number of integration points as specified by nip. More...

elementParallelMode	parallel_mode
	Determines the parallel mode of the element. More...

IntArray	partitions
	List of partition sharing the shared element or remote partition containing remote element counterpart. More...

Protected Attributes inherited from oofem::FEMComponent
int	number
	Component number. More...

Domain *	domain
	Link to domain object, useful for communicating with other FEM components. More...

Static Protected Attributes
static const double	stefanBoltzmann = 5.67e-8
	Stefan–Boltzmann constant W/m2/K4. More...

Detailed Description

This abstract class represent a general base element class for transport problems.

In actual implementation, the same approximation order of all unknowns is assumed, but this can be easily implemented.

Definition at line 51 of file transportelement.h.

Member Enumeration Documentation

enum oofem::TransportElement::ElementMode

Enumerator
HeatTransferEM
HeatMass1TransferEM
Mass1TransferEM

Definition at line 54 of file transportelement.h.

Constructor & Destructor Documentation

oofem::TransportElement::TransportElement	(	int	n,
		Domain *	d,
		ElementMode	em = `HeatTransferEM`
	)

Definition at line 69 of file transportelement.C.

oofem::TransportElement::~TransportElement ( )

virtual

Definition at line 75 of file transportelement.C.

Member Function Documentation

void oofem::TransportElement::assembleLocalContribution	(	FloatMatrix &	answer,
		FloatMatrix &	src,
		int	ndofs,
		int	rdof,
		int	cdof
	)

protected

Assembles the given source matrix of size (ndofs, ndofs) into target matrix answer.

The coefficients in src matrix are assumed to represent local sub-matrix for specific DOF. The answer matrix is element global one, and it typically contains the values for different DOFs. In this routine is assumed, that the number of dofs per node is the same. The DOFs in answer are ordered in a same way as the DOFS at the element level are (all Dofs corresponding to given node subsequently, followed by next node DOFS).

Parameters

answer	Receiver matrix.
src	Source matrix containing local sub matrix corresponding to single DOF.
ndofs	Number of DOFs per node (assumed same for all nodes).
rdof	Rows of src are localized into rows of answer corresponding to rdof-th dof.
cdof	Columns of src are localized into columns of answer corresponding to cdof-th dof.

Definition at line 1228 of file transportelement.C.

References oofem::FloatMatrix::at(), and oofem::Element::giveNumberOfDofManagers().

Referenced by computeBCMtrxAt(), computeBCVectorAt(), computeCapacityMatrix(), computeConductivityMatrix(), computeInternalSourceRhsVectorAt(), computeIntSourceLHSMatrix(), and computeSurfaceVolumeAround().

void oofem::TransportElement::assembleLocalContribution	(	FloatArray &	answer,
		FloatArray &	src,
		int	ndofs,
		int	rdof
	)

protected

Assembles the given source vector of size (ndofs) into target answer.

The coefficients in src vector are assumed to represent local sub-vector for specific DOF. The answer vector is element global one, and it typically contains the values for different DOFs. In this routine is assumed, that the number of dofs per node is the same. The DOFs in answer are ordered in a same way as the DOFS at the element level are (all Dofs corresponding to given node subsequently, followed by next node DOFS).

Parameters

answer	Receiver vector.
src	Source vector containing local subvector corresponding to single DOF.
ndofs	Number of DOFs per node (assumed same for all nodes).
rdof	Rows of src are localized into rows of answer corresponding to rdof-th dof.

Definition at line 1244 of file transportelement.C.

References oofem::FloatArray::at(), and oofem::Element::giveNumberOfDofManagers().

int oofem::TransportElement::checkConsistency ( )

virtual

Performs consistency check.

This method is called at startup for all elements in particular domain. This method is intended to check data compatibility. Particular element types should test if compatible material and crossSection both with required capabilities are specified. Derived classes should provide their own analysis specific tests. Some printed input if incompatibility is found should be provided (error or warning member functions). Method can be also used to initialize some variables, since this is invoked after all domain components are instanciated.

Returns: Zero value if check fail, otherwise nonzero.

Reimplemented from oofem::Element.

Definition at line 150 of file transportelement.C.

References oofem::Element::giveCrossSection(), giveMaterial(), oofem::Element::material, and OOFEM_WARNING.

Referenced by giveThicknessAt().

void oofem::TransportElement::computeBCMtrxAt	(	FloatMatrix &	answer,
		TimeStep *	tStep,
		ValueModeType	mode
	)

virtual

Computes the LHS contribution to balance equation(s) due to boundary conditions.

Definition at line 885 of file transportelement.C.

References assembleLocalContribution(), computeBCSubMtrxAt(), oofem::Element::computeNumberOfDofs(), emode, HeatMass1TransferEM, HeatTransferEM, oofem::FloatMatrix::isNotEmpty(), Mass1TransferEM, OOFEM_ERROR, oofem::FloatMatrix::resize(), and oofem::FloatMatrix::zero().

Referenced by giveCharacteristicMatrix(), and giveThicknessAt().

void oofem::TransportElement::computeBCSubMtrxAt	(	FloatMatrix &	answer,
		TimeStep *	tStep,
		ValueModeType	mode,
		int	indx
	)

protected

Computes the part of LHS due to applied BCs.

the result should be localized.

Parameters

answer	Contains the result.
tStep	Solution step.
mode	Mode of value (incremental, total, ...).
indx	Unknown index.

Todo:: Deal with coupled fields

Todo:: Deal with coupled fields

Definition at line 1115 of file transportelement.C.

References oofem::FloatMatrix::assemble(), oofem::IntArray::at(), oofem::Element::boundaryLoadArray, oofem::FloatMatrix::clear(), oofem::Element::computeBoundaryVectorOf(), computeEdgeVolumeAround(), computeEgdeNAt(), computeSurfaceNAt(), computeSurfaceVolumeAround(), oofem::ConvectionBC, oofem::FEMComponent::domain, oofem::FloatArray::dotProduct(), oofem::EdgeLoadBGT, emode, GetSurfaceIntegrationRule(), giveApproxOrder(), oofem::GeneralBoundaryCondition::giveBCGeoType(), oofem::Element::giveBoundaryLoadArray(), giveEdgeDofMapping(), oofem::Element::giveElementDofIDMask(), oofem::Domain::giveLoad(), oofem::Element::giveNumberOfDofManagers(), oofem::BoundaryLoad::giveProperty(), oofem::IntArray::giveSize(), giveSurfaceDofMapping(), oofem::GeneralBoundaryCondition::giveType(), HeatTransferEM, oofem::ScalarFunction::isDefined(), oofem::Load::isDofExcluded(), oofem::GeneralBoundaryCondition::isImposed(), Mass1TransferEM, OOFEM_ERROR, oofem::FloatMatrix::plusDyadSymmUpper(), oofem::BoundaryLoad::propertyMultExpr, oofem::FloatMatrix::resize(), oofem::GaussIntegrationRule::SetUpPointsOnLine(), oofem::SurfaceLoadBGT, oofem::FloatMatrix::symmetrized(), and oofem::FloatMatrix::zero().

Referenced by computeBCMtrxAt(), and giveThicknessAt().

void oofem::TransportElement::computeBCSubVectorAt	(	FloatArray &	answer,
		TimeStep *	tStep,
		ValueModeType	mode,
		int	indx
	)

protected

Computes the part of RHS due to applied BCs.

The part corresponds to unknown identified by indx param. The result should be localized. The part corresponding to Dirichlet BC is not taken into account.

Parameters

answer	Contains the result.
tStep	Solution step.
mode	Mode of value (incremental, total, ...).
indx	Unknown index.

Definition at line 909 of file transportelement.C.

References oofem::FloatArray::add(), oofem::IntArray::at(), oofem::Element::bodyLoadArray, oofem::Element::boundaryLoadArray, computeBodyBCSubVectorAt(), computeEdgeBCSubVectorAt(), computeSurfaceBCSubVectorAt(), oofem::FEMComponent::domain, oofem::EdgeLoadBGT, oofem::GeneralBoundaryCondition::giveBCGeoType(), oofem::Element::giveBoundaryLoadArray(), oofem::Domain::giveLoad(), oofem::Element::giveNumberOfDofManagers(), oofem::IntArray::giveSize(), oofem::Load::isDofExcluded(), oofem::GeneralBoundaryCondition::isImposed(), OOFEM_ERROR, oofem::FloatArray::resize(), oofem::SurfaceLoadBGT, and oofem::FloatArray::zero().

Referenced by computeBCVectorAt(), and giveThicknessAt().

void oofem::TransportElement::computeBCVectorAt	(	FloatArray &	answer,
		TimeStep *	tStep,
		ValueModeType	mode
	)

virtual

Computes the RHS contribution to balance equation(s) due to boundary conditions.

Definition at line 864 of file transportelement.C.

References assembleLocalContribution(), computeBCSubVectorAt(), oofem::Element::computeNumberOfDofs(), emode, HeatMass1TransferEM, HeatTransferEM, Mass1TransferEM, OOFEM_ERROR, oofem::FloatArray::resize(), and oofem::FloatArray::zero().

Referenced by computeExternalForcesVector(), giveThicknessAt(), and oofem::TransportExternalForceAssembler::vectorFromElement().

void oofem::TransportElement::computeBmatrixAt	(	FloatMatrix &	answer,
		const FloatArray &	lcoords
	)

protectedvirtual

Todo:: We should change the transposition in evaldNdx;

Reimplemented in oofem::Lattice2d_mt.

Definition at line 241 of file transportelement.C.

References oofem::FloatMatrix::beTranspositionOf(), emode, oofem::FEInterpolation::evaldNdx(), oofem::Element::giveInterpolation(), oofem::FloatMatrix::giveNumberOfColumns(), oofem::FloatMatrix::giveNumberOfRows(), HeatMass1TransferEM, HeatTransferEM, Mass1TransferEM, oofem::FloatMatrix::resize(), and oofem::FloatMatrix::zero().

Referenced by computeInternalForcesVector(), and giveThicknessAt().

void oofem::TransportElement::computeBodyBCSubVectorAt	(	FloatArray &	answer,
		Load *	load,
		TimeStep *	tStep,
		ValueModeType	mode,
		int	indx
	)

protected

Definition at line 949 of file transportelement.C.

References oofem::FloatArray::add(), oofem::FloatArray::at(), oofem::Element::computeGlobalCoordinates(), computeNAt(), oofem::Load::computeValueAt(), oofem::Element::computeVolumeAround(), oofem::Load::giveApproxOrder(), oofem::Element::giveIntegrationRule(), oofem::Element::giveInterpolation(), oofem::Element::giveNumberOfDofManagers(), and oofem::FloatArray::resize().

Referenced by computeBCSubVectorAt(), and giveThicknessAt().

void oofem::TransportElement::computeBoundaryEdgeLoadVector	(	FloatArray &	answer,
		BoundaryLoad *	load,
		int	edge,
		CharType	type,
		ValueModeType	mode,
		TimeStep *	tStep,
		bool	global = `true`
	)

virtual

Computes the contribution of the given load at the given boundary edge.

In general, the answer should include only relevant DOFs at the edge. The related is giveBoundaryLocationArray method, which should return corresponding code numbers..

Note: Elements which do not have an contribution should resize the vector to be empty.

Parameters

answer	Requested contribution of load (in Global c.s.).
load	Load to compute contribution from.
edge	Edge number.
type	Type of the contribution.
mode	Determines mode of answer.
tStep	Time step when answer is computed.
global	if true (default) then contribution is in global c.s., when false then contribution is in element local c.s.

Todo:: Very similar to function computeBoundaryLoadVector and computeEdgeBCSubVectorAt. They should be all merged together.

Todo:: Deal with coupled fields (I think they should be another class of problems completely).

Todo:: Change the load type to "BoundaryEdgeLoad" maybe?

Reimplemented from oofem::Element.

Definition at line 778 of file transportelement.C.

void oofem::TransportElement::computeBoundarySurfaceLoadVector	(	FloatArray &	answer,
		BoundaryLoad *	load,
		int	boundary,
		CharType	type,
		ValueModeType	mode,
		TimeStep *	tStep,
		bool	global = `true`
	)

virtual

Computes the contribution of the given load at the given boundary surface in global coordinate system.

In general, the answer should include only relevant DOFs at the edge. The related is giveBoundaryLocationArray method, which should return corresponding code numbers.

Note: Elements which do not have an contribution should resize the vector to be empty.

Parameters

answer	Requested contribution of load.
load	Load to compute contribution from.
boundary	Boundary number.
type	Type of the contribution.
mode	Determines mode of answer.
tStep	Time step when answer is computed.
global	if true (default) then contribution is in global c.s., when false then contribution is in element local c.s.

Todo:: Deal with coupled fields (I think they should be another class of problems completely).

Reimplemented from oofem::Element.

Definition at line 597 of file transportelement.C.

void oofem::TransportElement::computeCapacityMatrix	(	FloatMatrix &	answer,
		TimeStep *	tStep
	)

virtual

Computes the capacity matrix of the receiver.

Reimplemented in oofem::Lattice2d_mt.

Definition at line 174 of file transportelement.C.

References assembleLocalContribution(), computeCapacitySubMatrix(), oofem::Element::computeNumberOfDofs(), emode, HeatMass1TransferEM, HeatTransferEM, Mass1TransferEM, OOFEM_ERROR, oofem::FloatMatrix::resize(), and oofem::FloatMatrix::zero().

Referenced by computeInertiaForcesVector(), computeLumpedCapacityVector(), giveCharacteristicMatrix(), and giveThicknessAt().

void oofem::TransportElement::computeCapacitySubMatrix	(	FloatMatrix &	answer,
		MatResponseMode	rmode,
		int	iri,
		TimeStep *	tStep
	)

protectedvirtual

Computes the matrix $P=\int_\Omega N^{\mathrm{T}} N\;\mathrm{d}\Omega$ .

The result should be multiplied by corresponding coefficient and localized.

Parameters

answer	Contains the result.
rmode	Determines the capacity coefficient to be used.
iri	Index of integration rule to use.
tStep	Time step.

Definition at line 322 of file transportelement.C.

References oofem::FloatMatrix::clear(), computeNAt(), oofem::Element::computeVolumeAround(), oofem::TransportMaterial::giveCharacteristicValue(), giveMaterial(), oofem::Element::integrationRulesArray, oofem::FloatMatrix::plusDyadSymmUpper(), and oofem::FloatMatrix::symmetrized().

Referenced by computeCapacityMatrix(), and giveThicknessAt().

void oofem::TransportElement::computeConductivityMatrix	(	FloatMatrix &	answer,
		MatResponseMode	rMode,
		TimeStep *	tStep
	)

virtual

Computes the conductivity matrix of the receiver.

Reimplemented in oofem::Lattice2d_mt.

Definition at line 198 of file transportelement.C.

References assembleLocalContribution(), computeConductivitySubMatrix(), oofem::Element::computeNumberOfDofs(), emode, HeatMass1TransferEM, HeatTransferEM, Mass1TransferEM, OOFEM_ERROR, oofem::FloatMatrix::resize(), and oofem::FloatMatrix::zero().

Referenced by giveCharacteristicMatrix(), and giveThicknessAt().

void oofem::TransportElement::computeConductivitySubMatrix	(	FloatMatrix &	answer,
		int	iri,
		MatResponseMode	rmode,
		TimeStep *	tStep
	)

protectedvirtual

Computes the matrix $C=\int_\Omega \sum_i \left( \frac{\partial N}{\partial x_i}\right)^{\mathrm{T}} \frac{\partial N}{\partial x_i} \;\mathrm{d}\Omega$ .

The result should be multiplied by corresponding coefficient and localized.

Todo:: Is this expression really correct?

Parameters

answer	Contains the result.
iri	Index of integration rule to use.
rmode	Determines the constitutive sub matrix to be used.
tStep	Solution step.

Definition at line 340 of file transportelement.C.

References oofem::FloatMatrix::beProductOf(), computeConstitutiveMatrixAt(), computeGradientMatrixAt(), oofem::Element::computeVolumeAround(), oofem::Element::giveNumberOfDofManagers(), oofem::Element::integrationRulesArray, oofem::FloatMatrix::plusProductSymmUpper(), oofem::FloatMatrix::resize(), oofem::FloatMatrix::symmetrized(), and oofem::FloatMatrix::zero().

Referenced by computeConductivityMatrix(), and giveThicknessAt().

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

virtual

Computes constitutive matrix of receiver.

Parameters

answer	Computed answer.
rMode	Material response mode of answer.
gp	Integration point for which constitutive matrix is computed.
tStep	Time step.

Definition at line 464 of file transportelement.C.

References giveCharacteristicMatrix(), and giveMaterial().

Referenced by computeConductivitySubMatrix(), computeFlow(), and giveThicknessAt().

void oofem::TransportElement::computeEdgeBCSubVectorAt	(	FloatArray &	answer,
		Load *	load,
		int	iEdge,
		TimeStep *	tStep,
		ValueModeType	mode,
		int	indx
	)

protected

Computes the part of RHS due to applied BCs on particular edge.

The part corresponds to unknown identified by indx param. The result should be localized. The part corresponding to Dirichlet BC is not taken into account.

Parameters

answer	Contains the result.
load	Edge load object.
iEdge	Element edge number subjected to bc.
tStep	Solution step.
mode	Mode of value (incremental, total, ...).
indx	Unknown index.

Todo:: Deal with coupled fields

Definition at line 967 of file transportelement.C.

References oofem::FloatArray::add(), oofem::FloatArray::assemble(), oofem::FloatArray::at(), oofem::FEInterpolation::boundaryEdgeLocal2Global(), oofem::Element::computeBoundaryVectorOf(), computeEdgeVolumeAround(), computeEgdeNAt(), oofem::BoundaryLoad::computeValueAt(), oofem::ConvectionBC, oofem::FEMComponent::domain, oofem::FloatArray::dotProduct(), emode, oofem::Load::FT_Entity, oofem::BoundaryLoad::giveApproxOrder(), giveApproxOrder(), oofem::EngngModel::giveContext(), giveEdgeDofMapping(), oofem::Element::giveElementDofIDMask(), oofem::Domain::giveEngngModel(), oofem::FieldManager::giveField(), oofem::EngngModelContext::giveFieldManager(), oofem::Load::giveFormulationType(), oofem::Element::giveInterpolation(), oofem::Element::giveNumberOfDofManagers(), oofem::BoundaryLoad::giveProperty(), oofem::GeneralBoundaryCondition::giveType(), HeatTransferEM, oofem::ScalarFunction::isDefined(), Mass1TransferEM, OOFEM_ERROR, oofem::BoundaryLoad::propertyMultExpr, oofem::RadiationBC, oofem::FloatArray::resize(), oofem::GaussIntegrationRule::SetUpPointsOnLine(), oofem::TransmissionBC, and oofem::FloatArray::zero().

Referenced by computeBCSubVectorAt(), and giveThicknessAt().

virtual double oofem::TransportElement::computeEdgeVolumeAround	(	GaussPoint *	gp,
		int	iEdge
	)

protectedpure virtual

Computes the length around a integration point on a edge.

Implemented in oofem::Lattice2d_mt, oofem::Brick1_ht, oofem::QWedge_ht, oofem::Wedge_ht, oofem::Tetrah1_ht, oofem::Quad1_ht, oofem::Tr1_ht, oofem::QBrick1_ht, oofem::Tr1Darcy, oofem::QuadAxisym1_ht, and oofem::TrAxisym1_ht.

Referenced by computeBCSubMtrxAt(), computeBoundaryEdgeLoadVector(), computeEdgeBCSubVectorAt(), and giveThicknessAt().

void oofem::TransportElement::computeEgdeNAt	(	FloatArray &	answer,
		int	iEdge,
		const FloatArray &	lcoord
	)

protectedvirtual

Computes the basis functions at the edge for one unknown.

Definition at line 274 of file transportelement.C.

References oofem::Element::giveInterpolation().

Referenced by computeBCSubMtrxAt(), computeEdgeBCSubVectorAt(), and giveThicknessAt().

void oofem::TransportElement::computeExternalForcesVector	(	FloatArray &	answer,
		TimeStep *	tStep,
		ValueModeType	mode
	)

virtual

Reimplemented in oofem::Tr1Darcy.

Definition at line 858 of file transportelement.C.

References computeBCVectorAt().

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

virtual

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 1256 of file transportelement.C.

References oofem::FloatArray::beProductOf(), computeNmatrixAt(), and oofem::Element::computeVectorOf().

void oofem::TransportElement::computeFlow	(	FloatArray &	answer,
		GaussPoint *	gp,
		TimeStep *	tStep
	)

virtual

Computes a flow vector in an integration point.

Parameters

answer	Flow vector.
gp	Integration point.
tStep	Time step.

Definition at line 1266 of file transportelement.C.

References oofem::FloatArray::beProductOf(), computeConstitutiveMatrixAt(), computeGradientMatrixAt(), oofem::Element::computeVectorOf(), emode, oofem::Element::giveElementDofIDMask(), oofem::GaussPoint::giveNaturalCoordinates(), oofem::FloatMatrix::giveNumberOfColumns(), oofem::FloatMatrix::giveNumberOfRows(), oofem::FloatArray::giveSize(), HeatMass1TransferEM, HeatTransferEM, oofem::Element::isActivated(), Mass1TransferEM, oofem::FloatArray::negated(), OOFEM_ERROR, oofem::FloatArray::resize(), oofem::FloatMatrix::resize(), oofem::FloatMatrix::setSubMatrix(), and oofem::FloatArray::zero().

Referenced by oofem::TransportMaterial::giveIPValue(), giveThicknessAt(), and oofem::TransportMaterialStatus::printOutputAt().

void oofem::TransportElement::computeGradientMatrixAt	(	FloatMatrix &	answer,
		const FloatArray &	lcoords
	)

protectedvirtual

Computes the gradient matrix corresponding to one unknown.

Todo:: We should change the transposition in evaldNdx;

Reimplemented in oofem::Lattice2d_mt.

Definition at line 264 of file transportelement.C.

References oofem::FloatMatrix::beTranspositionOf(), oofem::FEInterpolation::evaldNdx(), and oofem::Element::giveInterpolation().

Referenced by computeConductivitySubMatrix(), computeFlow(), giveThicknessAt(), and updateInternalState().

void oofem::TransportElement::computeInertiaForcesVector	(	FloatArray &	answer,
		TimeStep *	tStep
	)

virtual

Definition at line 517 of file transportelement.C.

References oofem::FloatArray::beProductOf(), computeCapacityMatrix(), and oofem::Element::computeVectorOf().

Referenced by giveCharacteristicVector().

void oofem::TransportElement::computeInternalForcesVector	(	FloatArray &	answer,
		TimeStep *	tStep
	)

virtual

Todo:: Can/should this part not be built into the flux itself? Probably not? / Mikael

Reimplemented in oofem::Tr1Darcy.

Definition at line 472 of file transportelement.C.

References oofem::FloatArray::beProductOf(), oofem::FloatArray::clear(), computeBmatrixAt(), oofem::TransportMaterial::computeInternalSourceVector(), computeNmatrixAt(), oofem::Element::computeVectorOf(), oofem::Element::computeVolumeAround(), oofem::TransportMaterial::giveFluxVector(), giveMaterial(), oofem::TransportMaterial::hasInternalSource(), oofem::Element::integrationRulesArray, N, and oofem::FloatArray::plusProduct().

Referenced by giveCharacteristicVector().

void oofem::TransportElement::computeInternalSourceRhsSubVectorAt	(	FloatArray &	answer,
		TimeStep *	tStep,
		ValueModeType	mode,
		int	indx
	)

protectedvirtual

Computes the contribution to balance equation(s) due to internal sources.

Definition at line 359 of file transportelement.C.

References oofem::FloatArray::add(), oofem::FloatArray::at(), oofem::FloatArray::clear(), oofem::TransportMaterial::computeInternalSourceVector(), computeNAt(), oofem::Element::computeVolumeAround(), oofem::Element::giveDefaultIntegrationRulePtr(), giveMaterial(), and oofem::TransportMaterial::hasInternalSource().

Referenced by computeInternalSourceRhsVectorAt(), and giveThicknessAt().

void oofem::TransportElement::computeInternalSourceRhsVectorAt	(	FloatArray &	answer,
		TimeStep *	tStep,
		ValueModeType	mode
	)

virtual

Computes the contribution to balance equation(s) due to internal sources.

Reimplemented in oofem::Lattice2d_mt.

Definition at line 386 of file transportelement.C.

References assembleLocalContribution(), computeInternalSourceRhsSubVectorAt(), emode, oofem::FloatArray::giveSize(), HeatMass1TransferEM, HeatTransferEM, oofem::FloatArray::isEmpty(), oofem::FloatArray::isNotEmpty(), Mass1TransferEM, OOFEM_ERROR, oofem::FloatArray::resize(), and oofem::FloatArray::zero().

Referenced by giveThicknessAt(), and oofem::TransportExternalForceAssembler::vectorFromElement().

void oofem::TransportElement::computeIntSourceLHSMatrix	(	FloatMatrix &	answer,
		TimeStep *	tStep
	)

virtual

Computes the LHS contribution to balance equation(s) due to material internal source.

Definition at line 411 of file transportelement.C.

References assembleLocalContribution(), oofem::FloatMatrix::clear(), computeIntSourceLHSSubMatrix(), oofem::Element::computeNumberOfDofs(), emode, giveMaterial(), oofem::TransportMaterial::hasInternalSource(), HeatMass1TransferEM, HeatTransferEM, Mass1TransferEM, OOFEM_ERROR, oofem::FloatMatrix::resize(), and oofem::FloatMatrix::zero().

Referenced by giveThicknessAt().

void oofem::TransportElement::computeIntSourceLHSSubMatrix	(	FloatMatrix &	answer,
		MatResponseMode	rmode,
		int	iri,
		TimeStep *	tStep
	)

virtual

Computes the part of internal source LHS contribution corresponding to unknown identified by rmode parameter.

Definition at line 440 of file transportelement.C.

References oofem::FloatMatrix::clear(), computeNAt(), oofem::Element::computeVolumeAround(), oofem::TransportMaterial::giveCharacteristicValue(), giveMaterial(), oofem::Element::integrationRulesArray, oofem::FloatMatrix::plusDyadSymmUpper(), and oofem::FloatMatrix::symmetrized().

Referenced by computeIntSourceLHSMatrix(), and giveThicknessAt().

void oofem::TransportElement::computeLoadVector	(	FloatArray &	answer,
		BodyLoad *	load,
		CharType	type,
		ValueModeType	mode,
		TimeStep *	tStep
	)

virtual

Computes the contribution of the given body load (volumetric).

Parameters

answer	Requested contribution of load.
load	Load to compute contribution from.
type	Type of the contribution.
mode	Determines mode of answer.
tStep	Time step when answer is computed.

Todo:: We could technically support other types of volume loads. Even a bulk load is a type of transmission b.c.

Todo:: Deal with coupled fields (I think they should be another class of problems completely).

Reimplemented from oofem::Element.

Definition at line 545 of file transportelement.C.

References oofem::FloatMatrix::beNMatrixOf(), oofem::FloatArray::clear(), oofem::Load::computeValueAt(), oofem::Element::computeVectorOf(), oofem::Element::computeVolumeAround(), oofem::ConvectionBC, emode, oofem::FEInterpolation::evalN(), oofem::Load::FT_Entity, oofem::Load::giveApproxOrder(), oofem::Element::giveElementDofIDMask(), oofem::Load::giveFormulationType(), oofem::FEInterpolation::giveIntegrationRule(), oofem::Element::giveInterpolation(), oofem::FEInterpolation::giveInterpolationOrder(), oofem::GeneralBoundaryCondition::giveType(), HeatMass1TransferEM, oofem::FEInterpolation::local2global(), N, oofem::FloatArray::plusProduct(), oofem::RadiationBC, and oofem::TransmissionBC.

void oofem::TransportElement::computeLumpedCapacityVector	(	FloatArray &	answer,
		TimeStep *	tStep
	)

virtual

Definition at line 528 of file transportelement.C.

References oofem::FloatArray::at(), oofem::FloatMatrix::at(), computeCapacityMatrix(), oofem::FloatMatrix::giveNumberOfRows(), oofem::FloatArray::resize(), and oofem::FloatArray::zero().

Referenced by giveCharacteristicVector().

void oofem::TransportElement::computeNAt	(	FloatArray &	answer,
		const FloatArray &	lcoord
	)

protectedvirtual

Computes the basis functions.

Parameters

answer	The basis functions evaluated at lcoord.
lcoord	The local coordinate.

Definition at line 221 of file transportelement.C.

References oofem::FEInterpolation::evalN(), and oofem::Element::giveInterpolation().

Referenced by computeBodyBCSubVectorAt(), computeCapacitySubMatrix(), computeInternalSourceRhsSubVectorAt(), computeIntSourceLHSSubMatrix(), computeNmatrixAt(), and giveThicknessAt().

void oofem::TransportElement::computeNmatrixAt	(	FloatMatrix &	answer,
		const FloatArray &	lcoords
	)

protectedvirtual

Computes the interpolation matrix corresponding to all unknowns.

In the default implementation the same approximation order is assumed, but it can be extended.

Reimplemented in oofem::Lattice2d_mt.

Definition at line 228 of file transportelement.C.

References oofem::FloatMatrix::beNMatrixOf(), computeNAt(), emode, HeatTransferEM, and Mass1TransferEM.

Referenced by computeField(), computeInternalForcesVector(), EIPrimaryFieldI_evaluateFieldVectorAt(), giveThicknessAt(), and updateInternalState().

void oofem::TransportElement::computeSurfaceBCSubVectorAt	(	FloatArray &	answer,
		Load *	load,
		int	iSurf,
		TimeStep *	tStep,
		ValueModeType	mode,
		int	indx
	)

protected

Computes the part of RHS due to applied BCs on particular surface.

The part corresponds to unknown identified by indx param. The result should be localized. The part corresponding to Dirichlet BC is not taken into account.

Parameters

answer	Contains the result.
load	Surface load object.
iSurf	Element surface number subjected to bc.
tStep	Solution step.
mode	Mode of value (incremental, total, ...).
indx	Unknown index.

Todo:: Deal with coupled fields

Definition at line 1040 of file transportelement.C.

References oofem::FloatArray::add(), oofem::FloatArray::assemble(), oofem::FloatArray::at(), oofem::FEInterpolation::boundarySurfaceLocal2global(), oofem::Element::computeBoundaryVectorOf(), computeSurfaceNAt(), computeSurfaceVolumeAround(), oofem::BoundaryLoad::computeValueAt(), oofem::ConvectionBC, oofem::FEMComponent::domain, oofem::FloatArray::dotProduct(), oofem::Element_SurfaceLoadSupport, emode, oofem::Load::FT_Entity, GetSurfaceIntegrationRule(), oofem::BoundaryLoad::giveApproxOrder(), giveApproxOrder(), oofem::EngngModel::giveContext(), oofem::Element::giveElementDofIDMask(), oofem::Domain::giveEngngModel(), oofem::FieldManager::giveField(), oofem::EngngModelContext::giveFieldManager(), oofem::Load::giveFormulationType(), oofem::Element::giveInterpolation(), oofem::Element::giveNumberOfDofManagers(), oofem::BoundaryLoad::giveProperty(), giveSurfaceDofMapping(), oofem::GeneralBoundaryCondition::giveType(), HeatTransferEM, oofem::ScalarFunction::isDefined(), Mass1TransferEM, OOFEM_ERROR, oofem::BoundaryLoad::propertyMultExpr, oofem::RadiationBC, oofem::FloatArray::resize(), oofem::Element::testElementExtension(), oofem::TransmissionBC, and oofem::FloatArray::zero().

Referenced by computeBCSubVectorAt(), and giveThicknessAt().

void oofem::TransportElement::computeSurfaceNAt	(	FloatArray &	answer,
		int	iSurf,
		const FloatArray &	lcoord
	)

protectedvirtual

Definition at line 298 of file transportelement.C.

References oofem::Element::giveInterpolation().

Referenced by computeBCSubMtrxAt(), computeSurfaceBCSubVectorAt(), and GetSurfaceIntegrationRule().

virtual double oofem::TransportElement::computeSurfaceVolumeAround	(	GaussPoint *	gp,
		int	iSurf
	)

inlineprotectedvirtual

Reimplemented in oofem::Brick1_ht, oofem::Tetrah1_ht, and oofem::QBrick1_ht.

Definition at line 251 of file transportelement.h.

References assembleLocalContribution(), giveApproxOrder(), and giveSurfaceDofMapping().

Referenced by computeBCSubMtrxAt(), and computeSurfaceBCSubVectorAt().

void oofem::TransportElement::computeTangentFromEdgeLoad	(	FloatMatrix &	answer,
		EdgeLoad *	load,
		int	boundary,
		MatResponseMode	rmode,
		TimeStep *	tStep
	)

virtual

Computes the tangent contribution of the given load at the given boundary.

Note: Elements which do not have an contribution should resize the vector to be empty.

Parameters

answer	Requested contribution of load.
load	Load to compute contribution from.
boundary	Surface number.
rmode	Mode of the contribution.
tStep	Time step when answer is computed.

Todo:: Deal with coupled fields (I think they should be another class of problems completely).

Reimplemented from oofem::Element.

Definition at line 724 of file transportelement.C.

References oofem::FloatArray::at(), oofem::FloatMatrix::beNMatrixOf(), oofem::FloatArray::beProductOf(), oofem::FEInterpolation::boundaryEdgeGiveNodes(), oofem::FEInterpolation::boundaryEvalN(), oofem::FEInterpolation::boundaryGiveTransformationJacobian(), oofem::FEInterpolation::boundaryLocal2Global(), oofem::FloatMatrix::clear(), oofem::Element::computeBoundaryVectorOf(), oofem::ConvectionBC, emode, oofem::BoundaryLoad::giveApproxOrder(), oofem::Element::giveBoundaryEdgeIntegrationRule(), oofem::Element::giveElementDofIDMask(), oofem::Element::giveInterpolation(), oofem::FEInterpolation::giveInterpolationOrder(), oofem::BoundaryLoad::giveProperty(), giveThicknessAt(), oofem::BoundaryLoad::giveType(), HeatMass1TransferEM, oofem::ScalarFunction::isDefined(), N, OOFEM_ERROR, oofem::FloatMatrix::plusProductSymmUpper(), oofem::BoundaryLoad::propertyMultExpr, oofem::FloatArray::subtract(), and oofem::FloatMatrix::symmetrized().

void oofem::TransportElement::computeTangentFromSurfaceLoad	(	FloatMatrix &	answer,
		SurfaceLoad *	load,
		int	boundary,
		MatResponseMode	rmode,
		TimeStep *	tStep
	)

virtual

Computes the tangent contribution of the given load at the given boundary.

Note: Elements which do not have an contribution should resize the vector to be empty.

Parameters

answer	Requested contribution of load.
load	Load to compute contribution from.
boundary	Surface number.
rmode	Mode of the contribution.
tStep	Time step when answer is computed.

Todo:: Deal with coupled fields (I think they should be another class of problems completely).

Reimplemented from oofem::Element.

Definition at line 673 of file transportelement.C.

References oofem::FloatMatrix::beNMatrixOf(), oofem::FEInterpolation::boundaryEvalN(), oofem::FEInterpolation::boundaryGiveNodes(), oofem::FEInterpolation::boundaryGiveTransformationJacobian(), oofem::FEInterpolation::boundaryLocal2Global(), oofem::FloatMatrix::clear(), oofem::Element::computeBoundaryVectorOf(), oofem::ConvectionBC, oofem::FloatArray::dotProduct(), emode, oofem::BoundaryLoad::giveApproxOrder(), oofem::Element::giveBoundarySurfaceIntegrationRule(), oofem::Element::giveElementDofIDMask(), oofem::Element::giveInterpolation(), oofem::FEInterpolation::giveInterpolationOrder(), oofem::BoundaryLoad::giveProperty(), giveThicknessAt(), oofem::BoundaryLoad::giveType(), HeatMass1TransferEM, oofem::ScalarFunction::isDefined(), N, OOFEM_ERROR, oofem::FloatMatrix::plusProductSymmUpper(), oofem::BoundaryLoad::propertyMultExpr, and oofem::FloatMatrix::symmetrized().

int oofem::TransportElement::EIPrimaryFieldI_evaluateFieldVectorAt	(	FloatArray &	answer,
		PrimaryField &	pf,
		const FloatArray &	coords,
		IntArray &	dofId,
		ValueModeType	mode,
		TimeStep *	tStep
	)

virtual

Evaluates the value of field at given point of interest (should be located inside receiver's volume) using element interpolation.

Todo:

This should use local coordinates instead of having all elements search for it manually.

Shouldn't this just be part of Element? It's very much the core of functionality for elements.

Parameters

answer	Field evaluated at coordinate.
pf	Field to use for evaluation.
coords	Coordinate.
dofId	IDs of DOFs to evaluate.
mode	Mode of field.
tStep	Time step to evaluate at.

Returns: Zero if ok, nonzero when error encountered.

Implements oofem::EIPrimaryFieldInterface.

Definition at line 1360 of file transportelement.C.

References oofem::IntArray::at(), oofem::FloatArray::at(), oofem::FloatMatrix::at(), oofem::Element::computeLocalCoordinates(), computeNmatrixAt(), oofem::Element::computeVectorOf(), oofem::IntArray::findFirstIndexOf(), oofem::Element::giveElementDofIDMask(), oofem::IntArray::giveSize(), oofem::FloatArray::giveSize(), OOFEM_ERROR, oofem::FloatArray::resize(), and oofem::FloatArray::zero().

Referenced by giveThicknessAt().

virtual IntegrationRule* oofem::TransportElement::GetSurfaceIntegrationRule ( int approxOrder )

inlineprotectedvirtual

Reimplemented in oofem::Brick1_ht, oofem::Tetrah1_ht, and oofem::QBrick1_ht.

Definition at line 249 of file transportelement.h.

References computeSurfaceNAt().

Referenced by computeBCSubMtrxAt(), and computeSurfaceBCSubVectorAt().

int oofem::TransportElement::giveApproxOrder ( int unknownIndx )

protectedvirtual

Reimplemented in oofem::Lattice2d_mt, oofem::QuadAxisym1_ht, and oofem::TrAxisym1_ht.

Definition at line 316 of file transportelement.C.

References oofem::Element::giveInterpolation(), and oofem::FEInterpolation::giveInterpolationOrder().

Referenced by computeBCSubMtrxAt(), computeEdgeBCSubVectorAt(), computeSurfaceBCSubVectorAt(), and computeSurfaceVolumeAround().

void oofem::TransportElement::giveCharacteristicMatrix	(	FloatMatrix &	answer,
		CharType	type,
		TimeStep *	tStep
	)

virtual

Computes characteristic matrix of receiver of requested type in given time step.

Parameters

answer	Requested characteristic matrix (stiffness, tangent, ...). If element has no capability to compute requested type of characteristic matrix error function is invoked.
type	Id of characteristic component requested.
tStep	Time step when answer is computed.

Reimplemented from oofem::Element.

Reimplemented in oofem::Tr1Darcy.

Definition at line 95 of file transportelement.C.

References oofem::__CharTypeToString(), oofem::FloatMatrix::add(), oofem::FloatMatrix::at(), computeBCMtrxAt(), computeCapacityMatrix(), computeConductivityMatrix(), oofem::FloatMatrix::giveNumberOfColumns(), oofem::FloatMatrix::giveNumberOfRows(), and OOFEM_ERROR.

Referenced by computeConstitutiveMatrixAt().

void oofem::TransportElement::giveCharacteristicVector	(	FloatArray &	answer,
		CharType	type,
		ValueModeType	mode,
		TimeStep *	tStep
	)

virtual

Computes characteristic vector of receiver of requested type in given time step.

If element has no capability to compute requested type of characteristic vector error function is invoked.

Parameters

answer	Requested characteristic vector.
type	Id of characteristic component requested.
mode	Determines mode of answer.
tStep	Time step when answer is computed.

Reimplemented from oofem::Element.

Reimplemented in oofem::Tr1Darcy.

Definition at line 127 of file transportelement.C.

References oofem::__CharTypeToString(), computeInertiaForcesVector(), computeInternalForcesVector(), computeLumpedCapacityVector(), OOFEM_ERROR, and oofem::FloatArray::resize().

void oofem::TransportElement::giveDofManDofIDMask	(	int	inode,
		IntArray &	answer
	)		const

virtual

Returns dofmanager dof mask for node.

This mask defines the dofs which are used by element in node. Mask influences the code number ordering for particular node. Code numbers are ordered according to node order and dofs belonging to particular node are ordered according to this mask. If element requests dofs using node mask which are not in node then error is generated. This masking allows node to be shared by different elements with different dofs in same node. Elements local code numbers are extracted from node using this mask. Must be defined by particular element.

Parameters

inode	Mask is computed for local dofmanager with inode number.
answer	Mask for node.

Reimplemented from oofem::Element.

Reimplemented in oofem::Lattice2d_mt, and oofem::Tr1Darcy.

Definition at line 80 of file transportelement.C.

References emode, HeatMass1TransferEM, HeatTransferEM, Mass1TransferEM, and OOFEM_ERROR.

void oofem::TransportElement::giveEdgeDofMapping	(	IntArray &	mask,
		int	iEdge
	)

protectedvirtual

Gives the node indexes for given edge.

Definition at line 286 of file transportelement.C.

References oofem::Element::giveInterpolation().

Referenced by computeBCSubMtrxAt(), computeEdgeBCSubVectorAt(), and giveThicknessAt().

int oofem::TransportElement::giveInternalStateAtNode	(	FloatArray &	answer,
		InternalStateType	type,
		InternalStateMode	mode,
		int	node,
		TimeStep *	tStep
	)

virtual

Returns internal state variable (like stress,strain) at node of element in Reduced form, the way how is obtained is dependent on InternalValueType.

The value may be local, or smoothed using some recovery technique. Returns zero if element is unable to respond to request.

Parameters

answer	Contains result, zero sized if not supported.
type	Determines the internal variable requested (physical meaning).
mode	Determines the mode of variable (recovered, local, ...).
node	Node number, for which variable is required.
tStep	Time step.

Returns: Nonzero if o.k, zero otherwise.

Reimplemented from oofem::Element.

Definition at line 1421 of file transportelement.C.

References oofem::FloatArray::at(), oofem::DofManager::end(), oofem::DofManager::findDofWithDofId(), oofem::Element::giveInternalStateAtNode(), oofem::Element::giveNode(), and oofem::FloatArray::resize().

Referenced by oofem::Quad1_ht::drawScalar(), oofem::Tetrah1_ht::drawScalar(), oofem::Brick1_ht::drawScalar(), and giveThicknessAt().

Material * oofem::TransportElement::giveMaterial ( )

virtual

Returns: Reference to the associated material of element.

Reimplemented from oofem::Element.

Definition at line 1409 of file transportelement.C.

References oofem::FEMComponent::domain, oofem::Element::giveCrossSection(), oofem::Domain::giveMaterial(), and oofem::Element::material.

Referenced by oofem::NonStationaryTransportProblem::applyIC(), checkConsistency(), oofem::Lattice2d_mt::computeCapacityMatrix(), computeCapacitySubMatrix(), oofem::Lattice2d_mt::computeConductivityMatrix(), computeConstitutiveMatrixAt(), computeInternalForcesVector(), oofem::Tr1Darcy::computeInternalForcesVector(), computeInternalSourceRhsSubVectorAt(), oofem::Lattice2d_mt::computeInternalSourceRhsVectorAt(), computeIntSourceLHSMatrix(), computeIntSourceLHSSubMatrix(), oofem::Tr1Darcy::computeStiffnessMatrix(), oofem::Lattice2d_mt::updateInternalState(), updateInternalState(), and oofem::NonStationaryTransportProblem::updateYourself().

void oofem::TransportElement::giveSurfaceDofMapping	(	IntArray &	mask,
		int	iSurf
	)

protectedvirtual

Definition at line 307 of file transportelement.C.

References oofem::Element::giveInterpolation().

Referenced by computeBCSubMtrxAt(), computeSurfaceBCSubVectorAt(), and computeSurfaceVolumeAround().

virtual double oofem::TransportElement::giveThicknessAt ( const FloatArray & gcoords )

inlinevirtual

Gives the thickness at some global coordinate.

For solid elements, the value returned is 1.0 (which is the default implementation).

Parameters

gcoords Global coordinates.

Returns: Thickness of element at given coordinate.

Todo:: Move this into the base element?

Reimplemented in oofem::Tr1Darcy, oofem::Quad1_ht, oofem::Tr1_ht, oofem::QuadAxisym1_ht, and oofem::TrAxisym1_ht.

Definition at line 97 of file transportelement.h.

References checkConsistency(), computeBCMtrxAt(), computeBCSubMtrxAt(), computeBCSubVectorAt(), computeBCVectorAt(), computeBmatrixAt(), computeBodyBCSubVectorAt(), computeCapacityMatrix(), computeCapacitySubMatrix(), computeConductivityMatrix(), computeConductivitySubMatrix(), computeConstitutiveMatrixAt(), computeEdgeBCSubVectorAt(), computeEdgeVolumeAround(), computeEgdeNAt(), computeFlow(), computeGradientMatrixAt(), computeInternalSourceRhsSubVectorAt(), computeInternalSourceRhsVectorAt(), computeIntSourceLHSMatrix(), computeIntSourceLHSSubMatrix(), computeNAt(), computeNmatrixAt(), computeSurfaceBCSubVectorAt(), EIPrimaryFieldI_evaluateFieldVectorAt(), giveEdgeDofMapping(), giveInternalStateAtNode(), and updateInternalState().

Referenced by computeBoundarySurfaceLoadVector(), computeTangentFromEdgeLoad(), and computeTangentFromSurfaceLoad().

TransportCrossSection * oofem::TransportElement::giveTransportCrossSection ( )

Definition at line 1402 of file transportelement.C.

References oofem::Element::giveCrossSection().

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

virtual

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.