This class is the implementation of triangular PFEM element with linear (and equal order) interpolation of velocity and pressure fields. More...

#include <tr1_2d_pfem.h>

Inheritance diagram for oofem::TR1_2D_PFEM:

Collaboration diagram for oofem::TR1_2D_PFEM:

Public Member Functions
	TR1_2D_PFEM (int n, Domain *aDomain, int particle1, int particle2, int particle3, int mat, int cs)
	Constructor of TR1_2D_PFEM Element. More...

	~TR1_2D_PFEM ()
	Destructor. More...

virtual void	computeDiagonalMassMtrx (FloatArray &answer, TimeStep *tStep)
	Calculates diagonal mass matrix as vector. More...

virtual void	computeDiagonalMassMtrx (FloatMatrix &answer, TimeStep *tStep)
	Calculates diagonal mass matrix. More...

virtual double	computeCriticalTimeStep (TimeStep *tStep)
	Calculates critical time step. More...

virtual const char *	giveClassName () const

virtual const char *	giveInputRecordName () const

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

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

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

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

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

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

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

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

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

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

virtual Element *	giveElement ()

virtual int	giveInternalStateAtNode (FloatArray &answer, InternalStateType type, InternalStateMode mode, int node, TimeStep *atTime)
	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	drawRawGeometry (oofegGraphicContext &)

virtual void	drawScalar (oofegGraphicContext &context)

virtual FEInterpolation *	giveVelocityInterpolation ()
	Returns the interpolation for velocity. More...

virtual FEInterpolation *	givePressureInterpolation ()
	Returns the interpolation for the pressure. More...

virtual FEInterpolation *	giveInterpolation () const

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

virtual const IntArray &	giveVelocityDofMask () const
	Returns mask of velocity Dofs. More...

virtual const IntArray &	givePressureDofMask () const
	Returns mask of pressure Dofs. More...

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

	~PFEMElement2d ()
	Destructor. More...

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

	~PFEMElement ()
	Destructor. More...

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

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

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

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

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

Public Member Functions inherited from oofem::FMElement
	FMElement (int n, Domain *aDomain)

virtual	~FMElement ()

virtual void	updateStabilizationCoeffs (TimeStep *tStep)
	Updates the stabilization coefficients used for CBS and SUPG algorithms. More...

virtual void	computeVectorOfVelocities (ValueModeType mode, TimeStep *tStep, FloatArray &velocities)

virtual void	computeVectorOfPressures (ValueModeType mode, TimeStep *tStep, FloatArray &pressures)

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 void	giveInputRecord (DynamicInputRecord &input)
	Setups the input record string of receiver. More...

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

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

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	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	computeField (ValueModeType mode, TimeStep *tStep, const FloatArray &lcoords, FloatArray &answer)
	Computes the unknown vector interpolated at the specified local coordinates. 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 Material *	giveMaterial ()

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

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

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

Protected Member Functions
virtual void	computeGaussPoints ()
	Initializes the array of integration rules member variable. More...

virtual void	computeDeviatoricStress (FloatArray &answer, GaussPoint gp, TimeStep tStep)
	Computes deviatoric stress vector in given integration point and solution step from given total strain vector. More...

virtual void	computeDeviatoricStressDivergence (FloatArray &answer, TimeStep *tStep)
	Calculates the divergence of the deviatoric stress. More...

virtual void	computeBodyLoadVectorAt (FloatArray &answer, BodyLoad load, TimeStep tStep, ValueModeType mode)
	Computes the load vector due to body load acting on receiver, at given time step. More...

virtual void	computeEdgeBCSubVectorAt (FloatArray &answer, Load load, int iEdge, TimeStep tStep)

Protected Member Functions inherited from oofem::PFEMElement2d
virtual void	computeBMatrix (FloatMatrix &answer, GaussPoint *gp)
	Calculates the element shape function derivative matrix. More...

virtual void	computeStiffnessMatrix (FloatMatrix &answer, MatResponseMode mode, TimeStep *tStep)
	Calculates the stiffness matrix. More...

virtual void	computePressureLaplacianMatrix (FloatMatrix &answer, TimeStep *tStep)
	Calculates the pressure laplacian matrix. More...

virtual void	computeDivergenceMatrix (FloatMatrix &answerx, TimeStep *tStep)
	Calculates the velocity divergence matrix. More...

virtual void	computeGradientMatrix (FloatMatrix &answer, TimeStep *tStep)
	Calculates the pressure gradient matrix. More...

void	computePrescribedRhsVector (FloatArray &answer, TimeStep *tStep, ValueModeType mode)
	Calculates the prescribed velocity vector for the right hand side of the pressure equation. More...

void	computeEdgeNMatrixAt (FloatMatrix &answer, int iedge, GaussPoint *gp)
	Calculates the shape function matrix on an edge. More...

void	computeEgdeNVectorAt (FloatArray &answer, int iedge, GaussPoint *gp)
	Calculates the shape function vector on an edge. More...

double	computeEdgeVolumeAround (GaussPoint *gp, int iEdge)
	Calculates the volume around an edge. More...

void	giveEdgeDofMapping (IntArray &answer, int iEdge) const
	Gives the mapping for degrees of freedom on an edge. More...

Protected Attributes
double	b [3]

double	c [3]

double	area

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 FEI2dTrLin	velocityInterpolation
	Interpolation for velocity unknowns. More...

static FEI2dTrLin	pressureInterpolation
	Interpolation for pressure unknowns. More...

static IntArray	edge_ordering [3]

static IntArray	velocityDofMask
	Mask of velocity Dofs. More...

static IntArray	pressureDofMask
	Mask of pressure Dofs. More...

Detailed Description

This class is the implementation of triangular PFEM element with linear (and equal order) interpolation of velocity and pressure fields.

Author: David Krybus

Definition at line 68 of file tr1_2d_pfem.h.

Constructor & Destructor Documentation

oofem::TR1_2D_PFEM::TR1_2D_PFEM	(	int	n,
		Domain *	aDomain,
		int	particle1,
		int	particle2,
		int	particle3,
		int	mat,
		int	cs
	)

Constructor of TR1_2D_PFEM Element.

Creates an element with number n belonging to domain aDomain.

Parameters

n	Element's number
aDomain	Pointer to the domain to which element belongs.
particle1	number of first PFEMParticle
particle2	number of second PFEMParticle
particle3	number of third PFEMParticle
mat	number of associated Material
cs	number of associated CrossSection

Definition at line 78 of file tr1_2d_pfem.C.

References oofem::IntArray::at(), oofem::Element::material, oofem::Element::numberOfDofMans, oofem::Element::postInitialize(), oofem::Element::setBodyLoads(), oofem::Element::setCrossSection(), and oofem::Element::setDofManagers().

oofem::TR1_2D_PFEM::~TR1_2D_PFEM ( )

Destructor.

Definition at line 93 of file tr1_2d_pfem.C.

Member Function Documentation

int oofem::TR1_2D_PFEM::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::PFEMElement2d.

Definition at line 276 of file tr1_2d_pfem.C.

References area, b, c, oofem::PFEMElement2d::checkConsistency(), oofem::Node::giveCoordinate(), and oofem::Element::giveNode().

Referenced by giveGeometryType().

void oofem::TR1_2D_PFEM::computeBodyLoadVectorAt	(	FloatArray &	answer,
		BodyLoad *	load,
		TimeStep *	tStep,
		ValueModeType	mode
	)

protectedvirtual

Computes the load vector due to body load acting on receiver, at given time step.

Default implementation computes body load vector numerically as $l=\int_V N^{\mathrm{T}} f \rho\;\mathrm{d}V$ using default integration rule. Result is transformed to global c.s.

Parameters

answer	Computed load vector due to body load
load	Body load which contribution is computed.
tStep	Time step.
mode	determines the response mode

Implements oofem::PFEMElement.

Definition at line 141 of file tr1_2d_pfem.C.

References oofem::Load::computeComponentArrayAt(), oofem::FEI2dTrLin::evalN(), oofem::Element::giveCrossSection(), oofem::FloatArray::giveSize(), oofem::FEI2dTrLin::giveTransformationJacobian(), oofem::Element::integrationRulesArray, N, pressureInterpolation, oofem::FloatArray::resize(), and oofem::FloatArray::zero().

Referenced by givePressureDofMask().

double oofem::TR1_2D_PFEM::computeCriticalTimeStep ( TimeStep * tStep )

virtual

Calculates critical time step.

Implements oofem::PFEMElement2d.

Definition at line 307 of file tr1_2d_pfem.C.

References oofem::FloatArray::add(), oofem::FloatArray::at(), oofem::FloatArray::computeNorm(), oofem::Element::computeVectorOf(), oofem::FEMComponent::domain, oofem::FloatArray::dotProduct(), oofem::Node::giveCoordinate(), oofem::Domain::giveEngngModel(), oofem::Element::giveNode(), oofem::EngngModel::giveSolutionStepWhenIcApply(), oofem::TimeStep::giveTimeIncrement(), oofem::min(), oofem::FloatArray::negated(), and oofem::FloatArray::subtract().

Referenced by oofem::PFEM::giveNextStep().

void oofem::TR1_2D_PFEM::computeDeviatoricStress	(	FloatArray &	answer,
		GaussPoint *	gp,
		TimeStep *
	)

protectedvirtual

Computes deviatoric stress vector in given integration point and solution step from given total strain vector.

Implements oofem::PFEMElement2d.

Definition at line 243 of file tr1_2d_pfem.C.

References oofem::FloatArray::at(), b, c, oofem::FluidDynamicMaterial::computeDeviatoricStressVector(), oofem::Element::computeVectorOf(), oofem::Element::giveCrossSection(), and oofem::FloatArray::resize().

Referenced by computeDeviatoricStressDivergence(), and givePressureDofMask().

void oofem::TR1_2D_PFEM::computeDeviatoricStressDivergence	(	FloatArray &	answer,
		TimeStep *	atTime
	)

protectedvirtual

Calculates the divergence of the deviatoric stress.

Implements oofem::PFEMElement2d.

Definition at line 260 of file tr1_2d_pfem.C.

References area, oofem::FloatArray::at(), b, c, computeDeviatoricStress(), oofem::Element::integrationRulesArray, oofem::FloatArray::resize(), and oofem::FloatArray::zero().

Referenced by givePressureDofMask().

void oofem::TR1_2D_PFEM::computeDiagonalMassMtrx	(	FloatArray &	answer,
		TimeStep *
	)

virtual

Calculates diagonal mass matrix as vector.

Implements oofem::PFEMElement.

Definition at line 212 of file tr1_2d_pfem.C.

References area, oofem::FloatArray::at(), oofem::Element::giveCrossSection(), oofem::Element::integrationRulesArray, oofem::FloatArray::resize(), and oofem::FloatArray::zero().

void oofem::TR1_2D_PFEM::computeDiagonalMassMtrx	(	FloatMatrix &	answer,
		TimeStep *
	)

virtual

Calculates diagonal mass matrix.

Implements oofem::PFEMElement.

Definition at line 225 of file tr1_2d_pfem.C.

References area, oofem::FloatMatrix::at(), oofem::Element::giveCrossSection(), oofem::Element::integrationRulesArray, oofem::FloatMatrix::resize(), and oofem::FloatMatrix::zero().

void oofem::TR1_2D_PFEM::computeEdgeBCSubVectorAt	(	FloatArray &	answer,
		Load *	load,
		int	iEdge,
		TimeStep *	tStep
	)

protectedvirtual

void oofem::TR1_2D_PFEM::computeGaussPoints ( )

protectedvirtual

Initializes the array of integration rules member variable.

Element can have multiple integration rules for different tasks. For example structural element family class uses this feature to implement transparent support for reduced and selective integration of some strain components. Must be defined by terminator classes.

See also: IntegrationRule

Implements oofem::PFEMElement2d.

Definition at line 129 of file tr1_2d_pfem.C.

References oofem::_Triangle, and oofem::Element::integrationRulesArray.

Referenced by givePressureDofMask(), and initializeFrom().

int oofem::TR1_2D_PFEM::computeNumberOfDofs ( )

virtual

Computes or simply returns total number of element's local DOFs.

Must be defined by particular element.

Returns: Number of local DOFs of element.

Implements oofem::PFEMElement2d.

Definition at line 98 of file tr1_2d_pfem.C.

Referenced by giveGeometryType().

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

virtual

Returns volume related to given integration point.

Used typically in subroutines, that perform integration over element volume. Should be implemented by particular elements.

See also: GaussPoint

Parameters

gp	Integration point for which volume is computed.

Returns: Volume for integration point.

Reimplemented from oofem::Element.

Definition at line 434 of file tr1_2d_pfem.C.

References oofem::GaussPoint::giveNaturalCoordinates(), oofem::FEI2dTrLin::giveTransformationJacobian(), oofem::GaussPoint::giveWeight(), and velocityInterpolation.

Referenced by giveGeometryType().

void oofem::TR1_2D_PFEM::drawRawGeometry ( oofegGraphicContext & gc )

virtual

Definition at line 463 of file tr1_2d_pfem.C.

References oofem::oofegGraphicContext::getElementColor(), oofem::oofegGraphicContext::getElementEdgeColor(), oofem::Node::giveCoordinate(), oofem::Element::giveNode(), OOFEG_RAW_GEOMETRY_LAYER, OOFEG_RAW_GEOMETRY_WIDTH, and oofem::oofegGraphicContext::testElementGraphicActivity().

Referenced by giveElement().

void oofem::TR1_2D_PFEM::drawScalar ( oofegGraphicContext & context )

virtual

Definition at line 493 of file tr1_2d_pfem.C.

Referenced by giveElement().

virtual const char* oofem::TR1_2D_PFEM::giveClassName ( ) const

inlinevirtual

Returns: Class name of the receiver.

Reimplemented from oofem::PFEMElement2d.

Definition at line 107 of file tr1_2d_pfem.h.

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

Implements oofem::PFEMElement2d.

Definition at line 104 of file tr1_2d_pfem.C.

Referenced by giveElementDofIDMask(), and giveGeometryType().

virtual Element* oofem::TR1_2D_PFEM::giveElement ( )

inlinevirtual

Reimplemented from oofem::PFEMElement2d.

Definition at line 126 of file tr1_2d_pfem.h.

References drawRawGeometry(), drawScalar(), and giveInternalStateAtNode().

void oofem::TR1_2D_PFEM::giveElementDofIDMask ( IntArray & answer ) const

virtual

Returns element dof mask for node.

This mask defines the dof ordering of the element interpolation. Default implementation for most elements, with noteable exceptions such as XFEM and some types of shell elements.

Parameters

ut	Equation DOFs belong to.
answer	DOF mask for receiver.

Implements oofem::PFEMElement2d.

Definition at line 113 of file tr1_2d_pfem.C.

References giveDofManDofIDMask().

Referenced by giveGeometryType().

virtual Element_Geometry_Type oofem::TR1_2D_PFEM::giveGeometryType ( ) const

inlinevirtual

Returns the element geometry type.

This information is assumed to be of general interest, but it is required only for some specialized tasks.

Returns: Geometry type of element.

Reimplemented from oofem::PFEMElement2d.

Definition at line 110 of file tr1_2d_pfem.h.

References checkConsistency(), computeNumberOfDofs(), computeVolumeAround(), giveDofManDofIDMask(), giveElementDofIDMask(), giveInterface(), initializeFrom(), restoreContext(), and saveContext().

virtual const char* oofem::TR1_2D_PFEM::giveInputRecordName ( ) const

inlinevirtual

Returns: Input record name of the receiver.

Implements oofem::FEMComponent.

Definition at line 108 of file tr1_2d_pfem.h.

References _IFT_TR1_2D_PFEM_Name.

Interface * oofem::TR1_2D_PFEM::giveInterface ( InterfaceType t )

virtual

Interface requesting service.

See also: InterfaceType

Returns: Requested interface if implemented, otherwise NULL.

Implements oofem::PFEMElement2d.

Definition at line 237 of file tr1_2d_pfem.C.

Referenced by giveGeometryType().

int oofem::TR1_2D_PFEM::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.

Implements oofem::PFEMElement2d.

Definition at line 443 of file tr1_2d_pfem.C.

References oofem::FloatArray::at(), oofem::Element::giveCrossSection(), oofem::PFEMElement::giveInternalStateAtNode(), oofem::Element::integrationRulesArray, and oofem::FloatArray::resize().

Referenced by drawScalar(), and giveElement().

virtual FEInterpolation* oofem::TR1_2D_PFEM::giveInterpolation ( ) const

inlinevirtual

Returns: Interpolation of the element geometry, or NULL if none exist.

Reimplemented from oofem::Element.

Definition at line 141 of file tr1_2d_pfem.h.

References velocityInterpolation.

virtual FEInterpolation* oofem::TR1_2D_PFEM::giveInterpolation ( DofIDItem id ) const

inlinevirtual

Returns the interpolation for the specific dof id.

Special elements which uses a mixed interpolation should reimplement this method.

Parameters

id	ID of the dof for the for the requested interpolation.

Returns: Appropriate interpolation, or NULL if none exists.

Reimplemented from oofem::Element.

Definition at line 142 of file tr1_2d_pfem.h.

References velocityInterpolation.

virtual const IntArray& oofem::TR1_2D_PFEM::givePressureDofMask ( ) const

inlinevirtual

Returns mask of pressure Dofs.

Implements oofem::PFEMElement.

Definition at line 148 of file tr1_2d_pfem.h.

References computeBodyLoadVectorAt(), computeDeviatoricStress(), computeDeviatoricStressDivergence(), computeEdgeBCSubVectorAt(), computeGaussPoints(), and pressureDofMask.

virtual FEInterpolation* oofem::TR1_2D_PFEM::givePressureInterpolation ( )

inlinevirtual

Returns the interpolation for the pressure.

Implements oofem::PFEMElement2d.

Definition at line 139 of file tr1_2d_pfem.h.

References pressureInterpolation.

virtual const IntArray& oofem::TR1_2D_PFEM::giveVelocityDofMask ( ) const

inlinevirtual

Returns mask of velocity Dofs.

Implements oofem::PFEMElement.

Definition at line 144 of file tr1_2d_pfem.h.

References velocityDofMask.

virtual FEInterpolation* oofem::TR1_2D_PFEM::giveVelocityInterpolation ( )

inlinevirtual

Returns the interpolation for velocity.

Implements oofem::PFEMElement2d.

Definition at line 138 of file tr1_2d_pfem.h.

References velocityInterpolation.

IRResultType oofem::TR1_2D_PFEM::initializeFrom ( InputRecord * ir )

virtual

Initializes receiver acording to object description stored in input record.

Reimplemented from oofem::PFEMElement2d.

Definition at line 120 of file tr1_2d_pfem.C.

References computeGaussPoints(), and oofem::PFEMElement::initializeFrom().

Referenced by giveGeometryType().

contextIOResultType oofem::TR1_2D_PFEM::restoreContext	(	DataStream &	stream,
		ContextMode	mode,
		void *	obj
	)

virtual

Restores the receiver state previously written in stream.

See also: saveContext

Parameters

stream	Input stream.
mode	Determines amount of info available in stream (state, definition, ...).
obj	Special parameter for sending extra information.

Returns: contextIOResultType.

Exceptions

throws an ContextIOERR exception if error encountered.

Reimplemented from oofem::Element.

Definition at line 420 of file tr1_2d_pfem.C.

References oofem::CIO_OK, oofem::Element::restoreContext(), and THROW_CIOERR.

Referenced by giveGeometryType().

contextIOResultType oofem::TR1_2D_PFEM::saveContext	(	DataStream &	stream,
		ContextMode	mode,
		void *	obj
	)

virtual

Stores receiver state to output stream.

Parameters

stream	Output stream.
mode	Determines amount of info required in stream (state, definition, ...).
obj	Special parameter, used only to send particular integration point to material class version of this method.

Returns: contextIOResultType.

Exceptions

throws an ContextIOERR exception if error encountered.

Reimplemented from oofem::Element.

Definition at line 408 of file tr1_2d_pfem.C.

References oofem::CIO_OK, oofem::Element::saveContext(), and THROW_CIOERR.

Referenced by giveGeometryType().

Member Data Documentation

double oofem::TR1_2D_PFEM::area

protected

Definition at line 73 of file tr1_2d_pfem.h.

Referenced by checkConsistency(), computeDeviatoricStressDivergence(), and computeDiagonalMassMtrx().

double oofem::TR1_2D_PFEM::b[3]

protected

Definition at line 71 of file tr1_2d_pfem.h.

Referenced by checkConsistency(), computeDeviatoricStress(), and computeDeviatoricStressDivergence().

double oofem::TR1_2D_PFEM::c[3]

protected

Definition at line 72 of file tr1_2d_pfem.h.

Referenced by checkConsistency(), computeDeviatoricStress(), and computeDeviatoricStressDivergence().

IntArray oofem::TR1_2D_PFEM::edge_ordering

staticprotected

Initial value:

= {
    IntArray(6), IntArray(6), IntArray(6)
}

Definition at line 80 of file tr1_2d_pfem.h.

Referenced by computeEdgeBCSubVectorAt().

IntArray oofem::TR1_2D_PFEM::pressureDofMask

staticprotected

Initial value:

= {
    3, 6, 9
}

Mask of pressure Dofs.

Definition at line 85 of file tr1_2d_pfem.h.

Referenced by givePressureDofMask().

FEI2dTrLin oofem::TR1_2D_PFEM::pressureInterpolation

staticprotected

Interpolation for pressure unknowns.

Definition at line 78 of file tr1_2d_pfem.h.

Referenced by computeBodyLoadVectorAt(), computeEdgeBCSubVectorAt(), and givePressureInterpolation().

IntArray oofem::TR1_2D_PFEM::velocityDofMask

staticprotected

Initial value:

= {
    1, 2, 4, 5, 7, 8
}

Mask of velocity Dofs.

Definition at line 83 of file tr1_2d_pfem.h.

Referenced by giveVelocityDofMask().

FEI2dTrLin oofem::TR1_2D_PFEM::velocityInterpolation

staticprotected

Interpolation for velocity unknowns.

Definition at line 76 of file tr1_2d_pfem.h.

Referenced by computeVolumeAround(), giveInterpolation(), and giveVelocityInterpolation().

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

Public Member Functions

Protected Member Functions

Protected Attributes

Static Protected Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation