oofem::StructuralElementEvaluator Class Reference

#include <structuralelementevaluator.h>

Inheritance diagram for oofem::StructuralElementEvaluator:

Collaboration diagram for oofem::StructuralElementEvaluator:

Protected Member Functions
	StructuralElementEvaluator ()
virtual	~StructuralElementEvaluator ()
virtual void	giveCharacteristicMatrix (FloatMatrix &answer, CharType mtrx, TimeStep *tStep)
virtual void	giveCharacteristicVector (FloatArray &answer, CharType type, ValueModeType mode, TimeStep *tStep)
virtual IntegrationRule *	giveMassMtrxIntegrationRule ()
virtual void	giveMassMtrxIntegrationMask (IntArray &answer)
virtual void	computeLumpedMassMatrix (FloatMatrix &answer, TimeStep *tStep)
virtual void	computeConsistentMassMatrix (FloatMatrix &answer, TimeStep *tStep, double &mass)
virtual Element *	giveElement ()=0
virtual void	computeNMatrixAt (FloatMatrix &answer, GaussPoint *gp)=0
virtual void	computeBMatrixAt (FloatMatrix &answer, GaussPoint *gp)=0
virtual void	computeStiffnessMatrix (FloatMatrix &answer, MatResponseMode rMode, TimeStep *tStep)
virtual double	computeVolumeAround (GaussPoint *gp)
virtual void	giveInternalForcesVector (FloatArray &answer, TimeStep *tStep, bool useUpdatedGpRecord=false)
void	computeVectorOf (ValueModeType u, TimeStep *tStep, FloatArray &answer)
void	computeVectorOf (PrimaryField &field, const IntArray &dofIdMask, ValueModeType u, TimeStep *tStep, FloatArray &answer)
bool	isActivated (TimeStep *tStep)
void	updateInternalState (TimeStep *tStep)
virtual void	computeStressVector (FloatArray &answer, const FloatArray &strain, GaussPoint *gp, TimeStep *tStep)=0
virtual void	computeConstitutiveMatrixAt (FloatMatrix &answer, MatResponseMode rMode, GaussPoint *gp, TimeStep *tStep)=0
void	computeStrainVector (FloatArray &answer, GaussPoint *gp, TimeStep *tStep, FloatArray &u)
virtual int	giveIntegrationElementLocalCodeNumbers (IntArray &answer, Element *elem, IntegrationRule *ie)

Protected Attributes
FloatMatrix	rotationMatrix
int	rotationMatrixDefined
	Flag indicating if transformation matrix has been already computed.

Friends
void	drawIGAPatchDeformedGeometry (Element *elem, StructuralElementEvaluator *se, oofegGraphicContext &gc, TimeStep *tStep, UnknownType)

Detailed Description

This class represent a new concept on how to define elements. Traditionally, Elements are derived from problem specific class (structural element, for example) define their interpolation and implement methods to evaluate shape function matrix, geometrical matrix, etc. This new concept, here represented by StructuralElementEvaluator and derived classes allows to define all problem specific methods (including shape function matrix, geometrical matrix evaluation) to be defined once for all elements of the same type (plane stress elements, space3d elements, etc) just relying on services of finite element interpolation classes. Definition of particular element is then done simply by deriving if from evaluator and providing interpolation.

StructuralElementEvaluator - base class of all structural elements Individual elements supposed to be derived from StructuralElementEvaluator and IGAElement

Todo

Definition at line 56 of file structuralelementevaluator.h.

Constructor & Destructor Documentation

StructuralElementEvaluator()

oofem::StructuralElementEvaluator::StructuralElementEvaluator ( )

protected

Definition at line 51 of file structuralelementevaluator.C.

References rotationMatrix.

Referenced by computeConstitutiveMatrixAt(), drawIGAPatchDeformedGeometry, oofem::PlaneStressStructuralElementEvaluator::PlaneStressStructuralElementEvaluator(), and oofem::Space3dStructuralElementEvaluator::Space3dStructuralElementEvaluator().

~StructuralElementEvaluator()

virtual oofem::StructuralElementEvaluator::~StructuralElementEvaluator ( )

inlineprotectedvirtual

Definition at line 64 of file structuralelementevaluator.h.

Member Function Documentation

computeBMatrixAt()

virtual void oofem::StructuralElementEvaluator::computeBMatrixAt ( FloatMatrix & answer, GaussPoint * gp )

protectedpure virtual

Implemented in oofem::PlaneStressStructuralElementEvaluator, and oofem::Space3dStructuralElementEvaluator.

Referenced by computeStiffnessMatrix(), computeStrainVector(), and giveInternalForcesVector().

computeConsistentMassMatrix()

void oofem::StructuralElementEvaluator::computeConsistentMassMatrix ( FloatMatrix & answer, TimeStep * tStep, double & mass )

protectedvirtual

Computes consistent mass matrix of receiver using numerical integration over element volume. Mass matrix is computed as \( M = \int_V N^{\mathrm{T}} \rho N dV \), where \( N \) is displacement approximation matrix. The number of necessary integration points is determined using this->giveNumberOfIPForMassMtrxIntegration service. Only selected degrees of freedom participate in integration of mass matrix. This is described using dof mass integration mask. This mask is obtained from this->giveMassMtrxIntegrationgMask service. The nonzero mask value at i-th position indicates that i-th element DOF participates in mass matrix computation. The result is in element local coordinate system.

Parameters

answer	Mass matrix.
tStep	Time step.
mass	Total mass of receiver.

Definition at line 206 of file structuralelementevaluator.C.

References oofem::FloatMatrix::at(), oofem::IntArray::at(), computeNMatrixAt(), oofem::Element::computeNumberOfDofs(), computeVolumeAround(), oofem::CrossSection::give(), giveElement(), giveMassMtrxIntegrationMask(), giveMassMtrxIntegrationRule(), oofem::FloatMatrix::giveNumberOfRows(), oofem::StructuralElement::giveStructuralCrossSection(), isActivated(), oofem::IntArray::isEmpty(), OOFEM_ERROR, oofem::FloatMatrix::plusProductSymmUpper(), oofem::FloatMatrix::resize(), oofem::FloatMatrix::symmetrized(), and oofem::FloatMatrix::zero().

Referenced by computeLumpedMassMatrix(), and giveCharacteristicMatrix().

computeConstitutiveMatrixAt()

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

protectedpure virtual

Computes constitutive matrix of receiver.

Parameters

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

Implemented in oofem::PlaneStressStructuralElementEvaluator, and oofem::Space3dStructuralElementEvaluator.

References gc, and StructuralElementEvaluator().

Referenced by computeStiffnessMatrix().

computeLumpedMassMatrix()

void oofem::StructuralElementEvaluator::computeLumpedMassMatrix ( FloatMatrix & answer, TimeStep * tStep )

protectedvirtual

Computes lumped mass matrix of receiver. Default implementation returns lumped consistent mass matrix. Then returns lumped mass transformed into nodal coordinate system. The lumping procedure zeros all off-diagonal members and zeros also all diagonal members corresponding to non-displacement DOFs. Such diagonal matrix is then rescaled, to preserve the element mass. Requires the computeNmatrixAt and giveMassMtrxIntegrationgMask services to be implemented.

Parameters

answer	Lumped mass matrix.
tStep	Time step.

Definition at line 147 of file structuralelementevaluator.C.

References oofem::FloatMatrix::at(), oofem::IntArray::at(), computeConsistentMassMatrix(), oofem::Element::computeNumberOfDofs(), oofem::Element::giveDofManDofIDMask(), giveElement(), oofem::Element::giveNumberOfDofManagers(), oofem::FloatMatrix::giveNumberOfRows(), oofem::IntArray::giveSize(), isActivated(), OOFEM_ERROR, oofem::FloatMatrix::resize(), oofem::FloatMatrix::times(), and oofem::FloatMatrix::zero().

Referenced by giveCharacteristicMatrix().

computeNMatrixAt()

virtual void oofem::StructuralElementEvaluator::computeNMatrixAt ( FloatMatrix & answer, GaussPoint * gp )

protectedpure virtual

Computes the matrix for which the unknown field is obtained, typically [N1, 0, N2, 0, ...; 0, N1, 0, N2, ...].

Implemented in oofem::PlaneStressStructuralElementEvaluator, and oofem::Space3dStructuralElementEvaluator.

Referenced by computeConsistentMassMatrix(), and drawIGAPatchDeformedGeometry.

computeStiffnessMatrix()

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

protectedvirtual

Definition at line 372 of file structuralelementevaluator.C.

References oofem::FloatMatrix::assemble(), oofem::FloatMatrix::beProductOf(), oofem::FloatMatrix::clear(), computeBMatrixAt(), computeConstitutiveMatrixAt(), oofem::Element::computeNumberOfDofs(), computeVolumeAround(), oofem::Element_remote, oofem::Element::giveCrossSection(), giveElement(), giveIntegrationElementLocalCodeNumbers(), oofem::Element::giveIntegrationRule(), oofem::Element::giveInterpolation(), oofem::Element::giveNumberOfIntegrationRules(), oofem::FEInterpolation::hasSubPatchFormulation(), oofem::StructuralCrossSection::isCharacteristicMtrxSymmetric(), oofem::FloatMatrix::plusProductSymmUpper(), oofem::FloatMatrix::plusProductUnsym(), oofem::FloatMatrix::resize(), oofem::FloatMatrix::symmetrized(), and oofem::FloatMatrix::zero().

Referenced by giveCharacteristicMatrix().

computeStrainVector()

void oofem::StructuralElementEvaluator::computeStrainVector ( FloatArray & answer, GaussPoint * gp, TimeStep * tStep, FloatArray & u )

protected

Computes large strain constitutive matrix of receiver.

Parameters

answer	Constitutive matrix.
rMode	Material response mode of answer.
gp	Integration point for which constitutive matrix is computed.
tStep	Time step. Optimized version, allowing to pass element displacements as parameter. Standard version has a huge performance leak; in typical IGA element the element vector is VERY large and its querying for each point take more time than strain evaluation. And this has to be done for each integration point. This optimized version allows to assemble displacement vector only once (for all IP) and pass this vector as parameter

Definition at line 312 of file structuralelementevaluator.C.

References oofem::FloatArray::at(), oofem::IntArray::at(), oofem::FloatArray::beProductOf(), computeBMatrixAt(), giveElement(), giveIntegrationElementLocalCodeNumbers(), oofem::GaussPoint::giveIntegrationRule(), oofem::GaussPoint::giveMaterialMode(), oofem::FloatMatrix::giveNumberOfColumns(), oofem::IntArray::giveSize(), isActivated(), oofem::FloatArray::resize(), and oofem::FloatArray::zero().

Referenced by oofem::BsplinePlaneStressElement::drawScalar(), oofem::NURBSPlaneStressElement::drawScalar(), oofem::NURBSSpace3dElement::drawScalar(), oofem::TSplinePlaneStressElement::drawScalar(), giveInternalForcesVector(), and updateInternalState().

computeStressVector()

virtual void oofem::StructuralElementEvaluator::computeStressVector ( FloatArray & answer, const FloatArray & strain, GaussPoint * gp, TimeStep * tStep )

protectedpure virtual

Computes the stress vector.

Parameters

answer	Stress vector.
strain	Strain vector.
gp	Integration point for which stress is computed.
tStep	Time step.

Implemented in oofem::PlaneStressStructuralElementEvaluator, and oofem::Space3dStructuralElementEvaluator.

Referenced by giveInternalForcesVector(), and updateInternalState().

computeVectorOf() [1/2]

void oofem::StructuralElementEvaluator::computeVectorOf ( PrimaryField & field, const IntArray & dofIdMask, ValueModeType u, TimeStep * tStep, FloatArray & answer )

inlineprotected

Definition at line 121 of file structuralelementevaluator.h.

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

computeVectorOf() [2/2]

void oofem::StructuralElementEvaluator::computeVectorOf ( ValueModeType u, TimeStep * tStep, FloatArray & answer )

inlineprotected

Definition at line 118 of file structuralelementevaluator.h.

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

Referenced by drawIGAPatchDeformedGeometry.

computeVolumeAround()

virtual double oofem::StructuralElementEvaluator::computeVolumeAround ( GaussPoint * gp )

inlineprotectedvirtual

Reimplemented in oofem::PlaneStressStructuralElementEvaluator, and oofem::Space3dStructuralElementEvaluator.

Definition at line 116 of file structuralelementevaluator.h.

Referenced by computeConsistentMassMatrix(), computeStiffnessMatrix(), and giveInternalForcesVector().

giveCharacteristicMatrix()

void oofem::StructuralElementEvaluator::giveCharacteristicMatrix ( FloatMatrix & answer, CharType mtrx, TimeStep * tStep )

protectedvirtual

Reimplemented in oofem::BsplinePlaneStressElement, oofem::NURBSPlaneStressElement, oofem::NURBSSpace3dElement, and oofem::TSplinePlaneStressElement.

Definition at line 128 of file structuralelementevaluator.C.

References oofem::__CharTypeToString(), computeConsistentMassMatrix(), computeLumpedMassMatrix(), computeStiffnessMatrix(), and OOFEM_ERROR.

giveCharacteristicVector()

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

protectedvirtual

Todo

Only for total value mode type (?)

Todo

Only for total value mode type (?)

Reimplemented in oofem::BsplinePlaneStressElement, oofem::NURBSPlaneStressElement, oofem::NURBSSpace3dElement, and oofem::TSplinePlaneStressElement.

Definition at line 116 of file structuralelementevaluator.C.

References oofem::FloatArray::clear(), and giveInternalForcesVector().

giveElement()

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

protectedpure virtual

Implemented in oofem::BsplinePlaneStressElement, oofem::NURBSPlaneStressElement, oofem::NURBSSpace3dElement, and oofem::TSplinePlaneStressElement.

Referenced by oofem::Space3dStructuralElementEvaluator::computeBMatrixAt(), computeConsistentMassMatrix(), oofem::PlaneStressStructuralElementEvaluator::computeConstitutiveMatrixAt(), oofem::Space3dStructuralElementEvaluator::computeConstitutiveMatrixAt(), computeLumpedMassMatrix(), oofem::Space3dStructuralElementEvaluator::computeNMatrixAt(), computeStiffnessMatrix(), computeStrainVector(), oofem::PlaneStressStructuralElementEvaluator::computeStressVector(), oofem::Space3dStructuralElementEvaluator::computeStressVector(), computeVectorOf(), computeVectorOf(), oofem::PlaneStressStructuralElementEvaluator::computeVolumeAround(), oofem::Space3dStructuralElementEvaluator::computeVolumeAround(), giveInternalForcesVector(), and updateInternalState().

giveIntegrationElementLocalCodeNumbers()

int oofem::StructuralElementEvaluator::giveIntegrationElementLocalCodeNumbers ( IntArray & answer, Element * elem, IntegrationRule * ie )

protectedvirtual

Assembles the local element 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

Returns

Nonzero if integration rule code numbers differ from element code numbers

Definition at line 81 of file structuralelementevaluator.C.

References oofem::IntArray::at(), oofem::IntArray::clear(), oofem::IntArray::followedBy(), oofem::Element::giveDofManDofIDMask(), oofem::Element::giveGeometryType(), oofem::Element::giveInterpolation(), oofem::IGAIntegrationElement::giveKnotSpan(), oofem::FEInterpolation::giveKnotSpanBasisFuncMask(), oofem::FEInterpolation::giveNsd(), oofem::IntArray::giveSize(), oofem::FEInterpolation::hasSubPatchFormulation(), and oofem::IntArray::resize().

Referenced by computeStiffnessMatrix(), computeStrainVector(), drawIGAPatchDeformedGeometry, and giveInternalForcesVector().

giveInternalForcesVector()

void oofem::StructuralElementEvaluator::giveInternalForcesVector ( FloatArray & answer, TimeStep * tStep, bool useUpdatedGpRecord = false )

protectedvirtual

Todo

This part computes the B matrix again; Inefficient.

Definition at line 260 of file structuralelementevaluator.C.

References oofem::FloatArray::assemble(), oofem::FloatArray::clear(), computeBMatrixAt(), oofem::Element::computeNumberOfDofs(), computeStrainVector(), computeStressVector(), oofem::Element::computeVectorOf(), computeVolumeAround(), giveElement(), giveIntegrationElementLocalCodeNumbers(), oofem::Element::giveIntegrationRule(), oofem::Element::giveInterpolation(), oofem::Element::giveNumberOfIntegrationRules(), oofem::FloatArray::giveSize(), oofem::StructuralMaterialStatus::giveStressVector(), oofem::FEInterpolation::hasSubPatchFormulation(), isActivated(), oofem::FloatArray::plusProduct(), oofem::FloatArray::resize(), and oofem::FloatArray::zero().

Referenced by giveCharacteristicVector().

giveMassMtrxIntegrationMask()

virtual void oofem::StructuralElementEvaluator::giveMassMtrxIntegrationMask ( IntArray & answer )

inlineprotectedvirtual

Returns mask indicating, which unknowns (their type and ordering is the same as element unknown vector) participate in mass matrix integration. Nonzero value at i-th position indicates that corresponding row in interpolation matrix N will participate in mass matrix integration (typically only displacements are taken into account).

Parameters

answer

Integration mask, if zero sized, all unknowns participate. This is default.

Definition at line 81 of file structuralelementevaluator.h.

References oofem::IntArray::clear().

Referenced by computeConsistentMassMatrix().

giveMassMtrxIntegrationRule()

virtual IntegrationRule * oofem::StructuralElementEvaluator::giveMassMtrxIntegrationRule ( )

inlineprotectedvirtual

Returns the integration rule for mass matrices, if relevant.

Returns

Number of integration points for mass matrix.

Definition at line 72 of file structuralelementevaluator.h.

Referenced by computeConsistentMassMatrix().

isActivated()

bool oofem::StructuralElementEvaluator::isActivated ( TimeStep * tStep )

inlineprotected

Definition at line 124 of file structuralelementevaluator.h.

Referenced by computeConsistentMassMatrix(), computeLumpedMassMatrix(), computeStrainVector(), and giveInternalForcesVector().

updateInternalState()

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

protected

Definition at line 340 of file structuralelementevaluator.C.

References computeStrainVector(), computeStressVector(), oofem::Element::computeVectorOf(), oofem::Element_remote, giveElement(), oofem::Element::giveIntegrationRule(), oofem::Element::giveNumberOfIntegrationRules(), and oofem::FloatArray::subtract().

Friends And Related Symbol Documentation

drawIGAPatchDeformedGeometry

void drawIGAPatchDeformedGeometry ( Element * elem, StructuralElementEvaluator * se, oofegGraphicContext & gc, TimeStep * tStep, UnknownType  )

friend

Definition at line 1001 of file iga.C.

References oofem::FloatArray::at(), oofem::IntArray::at(), oofem::FloatArray::beProductOf(), computeNMatrixAt(), computeVectorOf(), gc, giveIntegrationElementLocalCodeNumbers(), oofem::Element::giveIntegrationRule(), oofem::GaussPoint::giveIntegrationRule(), oofem::Element::giveInterpolation(), oofem::IntegrationRule::giveKnotSpan(), oofem::FEInterpolation::giveKnotVector(), oofem::FEInterpolation::giveNsd(), oofem::Element::giveNumberOfIntegrationRules(), oofem::IntArray::giveSize(), oofem::FEInterpolation::local2global(), N, OOFEG_DEFORMED_GEOMETRY_LAYER, OOFEG_DEFORMED_GEOMETRY_WIDTH, OOFEM_ERROR, oofem::FloatArray::resize(), and StructuralElementEvaluator().

Member Data Documentation

rotationMatrix

FloatMatrix oofem::StructuralElementEvaluator::rotationMatrix

protected

Definition at line 59 of file structuralelementevaluator.h.

Referenced by StructuralElementEvaluator().

rotationMatrixDefined

int oofem::StructuralElementEvaluator::rotationMatrixDefined

protected

Flag indicating if transformation matrix has been already computed.

Definition at line 61 of file structuralelementevaluator.h.

---

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

• structuralelementevaluator.h
• structuralelementevaluator.C