|
OOFEM 2.4
OOFEM.org - Object Oriented Finite Element Solver
|
General stabilized SUPG/PSPG element for CFD analysis. More...
#include <supgelement.h>
Inheritance diagram for oofem::SUPGElement:
Collaboration diagram for oofem::SUPGElement:Public Member Functions | |
| SUPGElement (int n, Domain *aDomain) | |
| void | initializeFrom (InputRecord &ir, int priority) override |
| void | giveInputRecord (DynamicInputRecord &input) override |
| Setups the input record string of receiver. | |
| void | giveCharacteristicMatrix (FloatMatrix &answer, CharType type, TimeStep *tStep) override |
| Computes characteristic matrix of receiver of requested type in given time step. | |
| void | giveCharacteristicVector (FloatArray &answer, CharType type, ValueModeType mode, TimeStep *tStep) override |
| Computes characteristic vector of receiver of requested type in given time step. | |
| void | updateStabilizationCoeffs (TimeStep *tStep) override |
| Updates the stabilization coefficients used for CBS and SUPG algorithms. | |
| virtual void | updateElementForNewInterfacePosition (TimeStep *tStep) |
| virtual void | computeAccelerationTerm_MB (FloatMatrix &answer, TimeStep *tStep)=0 |
| Computes acceleration terms (generalized mass matrix with stabilization terms) for momentum balance equations(s). | |
| virtual void | computeAdvectionTerm_MB (FloatArray &answer, TimeStep *tStep)=0 |
| Computes nonlinear advection terms for momentum balance equations(s). | |
| virtual void | computeAdvectionDerivativeTerm_MB (FloatMatrix &answer, TimeStep *tStep)=0 |
| Computes the derivative of advection terms for momentum balance equations(s) with respect to nodal velocities. | |
| virtual void | computeDiffusionTerm_MB (FloatArray &answer, TimeStep *tStep)=0 |
| Computes diffusion terms for momentum balance equations(s). | |
| virtual void | computeDiffusionDerivativeTerm_MB (FloatMatrix &answer, MatResponseMode mode, TimeStep *tStep)=0 |
| Computes the derivative of diffusion terms for momentum balance equations(s) with respect to nodal velocities. | |
| virtual void | computePressureTerm_MB (FloatMatrix &answer, TimeStep *tStep)=0 |
| Computes pressure terms for momentum balance equations(s). | |
| virtual void | computeBCLhsPressureTerm_MC (FloatMatrix &answer, TimeStep *tStep) |
| Computes Lhs terms due to boundary conditions - pressure. | |
| virtual void | computeLSICStabilizationTerm_MB (FloatMatrix &answer, TimeStep *tStep)=0 |
| Computes SLIC stabilization term for momentum balance equation(s). | |
| virtual void | computeLinearAdvectionTerm_MC (FloatMatrix &answer, TimeStep *tStep)=0 |
| Computes the linear advection term for mass conservation equation. | |
| virtual void | computeAdvectionTerm_MC (FloatArray &answer, TimeStep *tStep)=0 |
| Computes advection terms for mass conservation equation. | |
| virtual void | computeAdvectionDerivativeTerm_MC (FloatMatrix &answer, TimeStep *tStep)=0 |
| Computes the derivative of advection terms for mass conservation equation with respect to nodal velocities. | |
| virtual void | computeDiffusionDerivativeTerm_MC (FloatMatrix &answer, TimeStep *tStep)=0 |
| Computes diffusion derivative terms for mass conservation equation. | |
| virtual void | computeDiffusionTerm_MC (FloatArray &answer, TimeStep *tStep)=0 |
| Computes diffusion terms for mass conservation equation. | |
| virtual void | computeAccelerationTerm_MC (FloatMatrix &answer, TimeStep *tStep)=0 |
| Computes acceleration terms for mass conservation equation. | |
| virtual void | computePressureTerm_MC (FloatMatrix &answer, TimeStep *tStep)=0 |
| Computes pressure terms for mass conservation equation. | |
| virtual void | computeBCLhsTerm_MB (FloatMatrix &answer, TimeStep *tStep) |
| Computes Lhs terms due to boundary conditions - velocity. | |
| virtual void | computeBCLhsPressureTerm_MB (FloatMatrix &answer, TimeStep *tStep) |
| Computes Lhs terms due to boundary conditions - pressure. | |
| virtual void | computeBCRhsTerm_MB (FloatArray &answer, TimeStep *tStep)=0 |
| Computes Rhs terms due to boundary conditions. | |
| virtual void | computeBCRhsTerm_MC (FloatArray &answer, TimeStep *tStep)=0 |
| Computes Rhs terms due to boundary conditions. | |
| virtual void | computeSlipWithFrictionBCTerm_MB (FloatMatrix &answer, Load *load, int side, TimeStep *tStep) |
| Computes Lhs term due to applied slip with friction bc. | |
| virtual void | computePenetrationWithResistanceBCTerm_MB (FloatMatrix &answer, Load *load, int side, TimeStep *tStep) |
| Computes Lhs contribution due to applied Penetration bc. | |
| virtual void | computeOutFlowBCTerm_MB (FloatMatrix &answer, int side, TimeStep *tStep) |
| Computes Lhs contribution due to outflow BC. | |
| virtual void | computeHomogenizedReinforceTerm_MB (FloatMatrix &answer, Load *load, TimeStep *tStep) |
| virtual void | computeHomogenizedReinforceTerm_MC (FloatMatrix &answer, Load *load, TimeStep *tStep) |
| virtual double | computeCriticalTimeStep (TimeStep *tStep)=0 |
| Computes the critical time increment. | |
| void | updateInternalState (TimeStep *tStep) override |
| Updates element state after equilibrium in time step has been reached. | |
| int | checkConsistency () override |
| Performs consistency check. | |
| int | giveInternalStateAtNode (FloatArray &answer, InternalStateType type, InternalStateMode mode, int node, TimeStep *tStep) override |
| Returns internal state variable (like stress,strain) at node of element in Reduced form, the way how is obtained is dependent on InternalValueType. | |
| virtual void | giveLocalVelocityDofMap (IntArray &map) |
| virtual void | giveLocalPressureDofMap (IntArray &map) |
Public Member Functions inherited from oofem::FMElement | |
| FMElement (int n, Domain *aDomain) | |
| void | computeVectorOfVelocities (ValueModeType mode, TimeStep *tStep, FloatArray &velocities) |
| void | computeVectorOfPressures (ValueModeType mode, TimeStep *tStep, FloatArray &pressures) |
| FloatArray | computeVectorOfVelocities (ValueModeType mode, TimeStep *tStep) |
| FloatArray | computeVectorOfPressures (ValueModeType mode, TimeStep *tStep) |
Public Member Functions inherited from oofem::Element | |
| Element (int n, Domain *aDomain) | |
| Constructor. | |
| Element (const Element &src)=delete | |
| Element & | operator= (const Element &src)=delete |
| virtual | ~Element () |
| Virtual destructor. | |
| void | giveLocationArray (IntArray &locationArray, const UnknownNumberingScheme &s, IntArray *dofIds=NULL) const |
| Returns the location array (array of code numbers) of receiver for given numbering scheme. | |
| 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. | |
| 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. | |
| virtual void | setInternalDofManager (int num, std::unique_ptr< DofManager > dm) |
| virtual double | giveCharacteristicValue (CharType type, TimeStep *tStep) |
| Computes characteristic value of receiver of requested type in given time step. | |
| virtual void | computeLoadVector (FloatArray &answer, BodyLoad *load, CharType type, ValueModeType mode, TimeStep *tStep) |
| Computes the contribution of the given body load (volumetric). | |
| 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. | |
| virtual void | computeTangentFromSurfaceLoad (FloatMatrix &answer, BoundaryLoad *load, int boundary, MatResponseMode rmode, TimeStep *tStep) |
| Computes the tangent contribution of the given load at the given boundary. | |
| virtual void | computeTangentFromEdgeLoad (FloatMatrix &answer, BoundaryLoad *load, int boundary, MatResponseMode rmode, TimeStep *tStep) |
| Computes the tangent contribution of the given load at the given boundary. | |
| 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. | |
| const IntArray & | giveBodyLoadList () const |
| Returns receiver list of bodyloads. | |
| const IntArray & | giveBoundaryLoadList () const |
| Returns receiver list of boundary loads. | |
| void | computeVectorOf (ValueModeType u, TimeStep *tStep, FloatArray &answer) |
| Returns local vector of unknowns. | |
| 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. | |
| void | computeVectorOf (PrimaryField &field, const IntArray &dofIDMask, ValueModeType u, TimeStep *tStep, FloatArray &answer, bool padding=false) |
| Returns local vector of unknowns. | |
| void | computeVectorOfPrescribed (ValueModeType u, TimeStep *tStep, FloatArray &answer) |
| Returns local vector of prescribed unknowns. | |
| void | computeVectorOfPrescribed (const IntArray &dofIDMask, ValueModeType type, TimeStep *tStep, FloatArray &answer) |
| Returns local vector of prescribed unknowns. | |
| virtual int | computeNumberOfDofs () |
| Computes or simply returns total number of element's local DOFs. | |
| virtual int | computeNumberOfGlobalDofs () |
| Computes the total number of element's global dofs. | |
| int | computeNumberOfPrimaryMasterDofs () |
| Computes the total number of element's primary master DOFs. | |
| virtual bool | computeGtoLRotationMatrix (FloatMatrix &answer) |
| Returns transformation matrix from global c.s. | |
| virtual bool | giveRotationMatrix (FloatMatrix &answer) |
| Transformation matrices updates rotation matrix between element-local and primary DOFs, taking into account nodal c.s. | |
| 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. | |
| virtual void | giveDofManDofIDMask (int inode, IntArray &answer) const |
| Returns dofmanager dof mask for node. | |
| virtual void | giveInternalDofManDofIDMask (int inode, IntArray &answer) const |
| Returns internal dofmanager dof mask for node. | |
| virtual void | giveElementDofIDMask (IntArray &answer) const |
| Returns element dof mask for node. | |
| virtual void | computeField (ValueModeType mode, TimeStep *tStep, const FloatArray &lcoords, FloatArray &answer) |
| Computes the unknown vector interpolated at the specified local coordinates. | |
| virtual double | computeVolumeAround (GaussPoint *gp) |
| Returns volume related to given integration point. | |
| virtual double | computeVolumeAreaOrLength () |
| Computes the volume, area or length of the element depending on its spatial dimension. | |
| double | computeMeanSize () |
| Computes the size of the element defined as its length. | |
| virtual double | computeVolume () |
| Computes the volume. | |
| virtual double | computeArea () |
| Computes the area (zero for all but 2d geometries). | |
| virtual double | computeLength () |
| Computes the length (zero for all but 1D geometries) | |
| virtual IntArray | giveBoundaryEdgeNodes (int boundary, bool includeHierarchical=false) const |
| Returns list of receiver boundary nodes for given edge. | |
| virtual IntArray | giveBoundarySurfaceNodes (int boundary, bool includeHierarchical=false) const |
| Returns list of receiver boundary nodes for given surface. | |
| virtual IntArray | giveBoundaryNodes (int boundary) const |
| Returns list of receiver boundary nodes for given boundary entity. | |
| virtual std::unique_ptr< IntegrationRule > | giveBoundaryEdgeIntegrationRule (int order, int boundary) |
| Returns boundary edge integration rule. | |
| virtual std::unique_ptr< IntegrationRule > | giveBoundarySurfaceIntegrationRule (int order, int boundary) |
| Returns boundary surface integration rule. | |
| int | giveDofManagerNumber (int i) const |
| Translates local to global indices for dof managers. | |
| 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. | |
| virtual ElementSide * | giveSide (int i) const |
| Returns reference to the i-th side of element. | |
| virtual FEInterpolation * | giveInterpolation () const |
| virtual FEInterpolation * | giveInterpolation (DofIDItem id) const |
| Returns the interpolation for the specific dof id. | |
| virtual const FEInterpolation * | getGeometryInterpolation () const |
| Returns the geometry interpolation for the receiver. | |
| virtual Material * | giveMaterial () |
| int | giveMaterialNumber () const |
| CrossSection * | giveCrossSection () |
| int | getActivityTimeFunctionNumber () |
| Gets the activity time function number. | |
| void | setActivityTimeFunctionNumber (int funcIndx) |
| Sets the activity time function of receiver. | |
| void | setMaterial (int matIndx) |
| Sets the material of receiver. | |
| virtual void | setCrossSection (int csIndx) |
| Sets the cross section model of receiver. | |
| virtual int | giveNumberOfDofManagers () const |
| void | setNumberOfDofManagers (int i) |
| Sets number of element dof managers. | |
| virtual int | giveNumberOfNodes () const |
| Returns number of nodes of receiver. | |
| void | setDofManagers (const IntArray &dmans) |
| Sets receiver dofManagers. | |
| void | setDofManager (int id, int dm) |
| Sets receiver specific dofManager. | |
| void | setBodyLoads (const IntArray &bodyLoads) |
| Sets receiver bodyLoadArray. | |
| void | setIntegrationRules (std ::vector< std ::unique_ptr< IntegrationRule > > irlist) |
| Sets integration rules. | |
| virtual integrationDomain | giveIntegrationDomain () const |
| Returns integration domain for receiver, used to initialize integration point over receiver volume. | |
| virtual MaterialMode | giveMaterialMode () |
| Returns material mode for receiver integration points. | |
| 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. | |
| int | giveRegionNumber () |
| virtual void | updateYourself (TimeStep *tStep) |
| Updates element state after equilibrium in time step has been reached. | |
| virtual void | initializeYourself (TimeStep *timeStepWhenICApply) |
| Initialization according to state given by initial conditions. | |
| virtual bool | isActivated (TimeStep *tStep) |
| virtual bool | isCast (TimeStep *tStep) |
| virtual void | initForNewStep () |
| Initializes receivers state to new time step. | |
| virtual Element_Geometry_Type | giveGeometryType () const =0 |
| Returns the element geometry type. | |
| virtual Element_Geometry_Type | giveEdgeGeometryType (int id) const |
| Returns the receiver edge geometry type. | |
| virtual Element_Geometry_Type | giveSurfaceGeometryType (int id) const |
| Returns the receiver surface geometry type. | |
| virtual int | giveSpatialDimension () |
| Returns the element spatial dimension (1, 2, or 3). | |
| virtual int | giveNumberOfBoundarySides () |
| Returns number of boundaries (entities of element_dimension-1: points, edges, surfaces) | |
| virtual int | giveNumberOfEdges () const |
| virtual int | giveNumberOfSurfaces () const |
| virtual int | giveDefaultIntegrationRule () const |
| Returns id of default integration rule. | |
| virtual IntegrationRule * | giveDefaultIntegrationRulePtr () |
| Access method for default integration rule. | |
| 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. | |
| virtual int | giveIPValue (FloatArray &answer, GaussPoint *gp, InternalStateType type, TimeStep *tStep) |
| Returns the integration point corresponding value in full form. | |
| 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. | |
| virtual double | giveCharacteristicLength (const FloatArray &normalToCrackPlane) |
| Returns the size of element in the given direction, in some cases adjusted (e.g. | |
| 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. | |
| 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. | |
| virtual double | giveCharacteristicSize (GaussPoint *gp, FloatArray &normalToCrackPlane, ElementCharSizeMethod method) |
| Returns characteristic element size for a given integration point and given direction. | |
| virtual double | giveParentElSize () const |
| Returns the size (length, area or volume depending on element type) of the parent element. | |
| virtual void | updateBeforeNonlocalAverage (TimeStep *tStep) |
| Updates internal element state (in all integration points of receiver) before nonlocal averaging takes place. | |
| virtual int | computeGlobalCoordinates (FloatArray &answer, const FloatArray &lcoords) |
| Computes the global coordinates from given element's local coordinates. | |
| virtual bool | computeLocalCoordinates (FloatArray &answer, const FloatArray &gcoords) |
| Computes the element local coordinates from given global coordinates. | |
| virtual int | giveLocalCoordinateSystem (FloatMatrix &answer) |
| Returns local coordinate system of receiver Required by material models with ortho- and anisotrophy. | |
| virtual void | giveLocalCoordinateSystemVector (InternalStateType isttype, FloatArray &answer) |
| Returns unit vector of local coordinate system of receiver If local system is equal to global one, returns unit vectors. | |
| virtual void | computeMidPlaneNormal (FloatArray &answer, const GaussPoint *gp) |
| Computes mid-plane normal of receiver at integration point. | |
| virtual int | adaptiveMap (Domain *oldd, TimeStep *tStep) |
| Initializes the internal state variables stored in all IPs according to state in given domain. | |
| 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. | |
| virtual int | adaptiveUpdate (TimeStep *tStep) |
| Updates the internal state variables stored in all IPs according to already mapped state. | |
| virtual int | adaptiveFinish (TimeStep *tStep) |
| Finishes the mapping for given time step. | |
| void | updateLocalNumbering (EntityRenumberingFunctor &f) override |
| Local renumbering support. | |
| template<class T > | |
| void | ipEvaluator (T *src, void(T ::*f)(GaussPoint *gp)) |
| Integration point evaluator, loops over receiver IP's and calls given function (passed as f parameter) on them. The IP is parameter to function f. | |
| template<class T , class S > | |
| void | ipEvaluator (T *src, void(T ::*f)(GaussPoint *, S &), S &_val) |
| Integration point evaluator, loops over receiver IP's and calls given function (passed as f parameter) on them. The IP is parameter to function f as well as additional array. | |
| virtual void | drawYourself (oofegGraphicContext &gc, TimeStep *tStep) |
| virtual void | drawAnnotation (oofegGraphicContext &gc, TimeStep *tStep) |
| virtual void | drawRawGeometry (oofegGraphicContext &gc, TimeStep *tStep) |
| virtual void | drawDeformedGeometry (oofegGraphicContext &gc, TimeStep *tStep, UnknownType) |
| virtual void | drawScalar (oofegGraphicContext &gc, TimeStep *tStep) |
| virtual void | drawSpecial (oofegGraphicContext &gc, TimeStep *tStep) |
| virtual void | giveLocalIntVarMaxMin (oofegGraphicContext &gc, TimeStep *tStep, double &emin, double &emax) |
| virtual int | giveInternalStateAtSide (FloatArray &answer, InternalStateType type, InternalStateMode mode, int side, TimeStep *tStep) |
| 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. | |
| virtual void | showSparseMtrxStructure (CharType mtrx, oofegGraphicContext &gc, TimeStep *tStep) |
| Shows sparse structure. | |
| virtual void | showExtendedSparseMtrxStructure (CharType mtrx, oofegGraphicContext &gc, TimeStep *tStep) |
| Shows extended sparse structure (for example, due to nonlocal interactions for tangent stiffness) | |
| int | giveLabel () const |
| int | giveGlobalNumber () const |
| void | setGlobalNumber (int num) |
| Sets receiver globally unique number. | |
| elementParallelMode | giveParallelMode () const |
| Return elementParallelMode of receiver. | |
| void | setParallelMode (elementParallelMode _mode) |
| Sets parallel mode of element. | |
| virtual elementParallelMode | giveKnotSpanParallelMode (int) const |
| Returns the parallel mode for particular knot span of the receiver. | |
| int | packUnknowns (DataStream &buff, TimeStep *tStep) |
| Pack all necessary data of element (according to its parallel_mode) integration points into given communication buffer. | |
| 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. | |
| int | estimatePackSize (DataStream &buff) |
| Estimates the necessary pack size to hold all packed data of receiver. | |
| const IntArray * | givePartitionList () const |
| Returns partition list of receiver. | |
| void | setPartitionList (IntArray &pl) |
| Sets partition list of receiver. | |
| 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. | |
| virtual double | giveRelativeSelfComputationalCost () |
| Returns the weight representing relative computational cost of receiver The reference element is triangular plane stress element. | |
| virtual double | predictRelativeRedistributionCost () |
| Returns the relative redistribution cost of the receiver. | |
| IntArray * | giveBodyLoadArray () |
| Returns array containing load numbers of loads acting on element. | |
| IntArray * | giveBoundaryLoadArray () |
| Returns array containing load numbers of boundary loads acting on element. | |
| void | initializeFinish () override |
| Finishes the initialization. | |
| void | postInitialize () override |
| Performs post initialization steps. | |
| void | saveContext (DataStream &stream, ContextMode mode) override |
| Stores receiver state to output stream. | |
| void | restoreContext (DataStream &stream, ContextMode mode) override |
| Restores the receiver state previously written in stream. | |
| const char * | giveClassName () const override |
| void | printOutputAt (FILE *file, TimeStep *tStep) override |
| Prints output of receiver to stream, for given time step. | |
| virtual const IntArray | giveLocation () |
| virtual void | recalculateCoordinates (int nodeNumber, FloatArray &coords) |
| void | setSharedEdgeID (int iedge, int globalID) |
| Set global edge ID corresponding to given local edge id. | |
| void | setSharedSurfaceID (int isurf, int globalID) |
| const IntArray * | giveSharedEdgeIDs () const |
| const IntArray * | giveSharedSurfaceIDs () const |
Public Member Functions inherited from oofem::FEMComponent | |
| FEMComponent (int n, Domain *d) | |
| Regular constructor, creates component with given number and belonging to given domain. | |
| virtual | ~FEMComponent ()=default |
| Virtual destructor. | |
| virtual const char * | giveInputRecordName () const =0 |
| Domain * | giveDomain () const |
| virtual void | setDomain (Domain *d) |
| Sets associated Domain. | |
| int | giveNumber () const |
| void | setNumber (int num) |
| Sets number of receiver. | |
| virtual void | initializeFrom (InputRecord &ir) |
| Initializes receiver according to object description stored in input record. | |
| virtual void | printYourself () |
| Prints receiver state on stdout. Useful for debugging. | |
| virtual Interface * | giveInterface (InterfaceType t) |
| Interface requesting service. | |
| std::string | errorInfo (const char *func) const |
| Returns string for prepending output (used by error reporting macros). | |
Protected Member Functions | |
| virtual void | computeDeviatoricStrain (FloatArray &answer, GaussPoint *gp, TimeStep *tStep)=0 |
| virtual void | computeDeviatoricStress (FloatArray &answer, const FloatArray &eps, GaussPoint *gp, TimeStep *tStep)=0 |
| virtual void | computeTangent (FloatMatrix &answer, MatResponseMode mode, GaussPoint *gp, TimeStep *tStep)=0 |
Protected Member Functions inherited from oofem::Element | |
| virtual void | computeGaussPoints () |
| Initializes the array of integration rules member variable. | |
Protected Attributes | |
| IntArray | boundarySides |
| Array of boundary sides. | |
| IntArray | boundaryCodes |
| Boundary sides codes. | |
| double | t_supg = 0. |
| Stabilization coefficients, updated for each solution step in updateStabilizationCoeffs() | |
| double | t_pspg = 0. |
| double | t_lsic = 0. |
Protected Attributes inherited from oofem::Element | |
| int | numberOfDofMans |
| Number of dofmanagers. | |
| IntArray | dofManArray |
| Array containing dofmanager numbers. | |
| int | material |
| Number of associated material. | |
| int | crossSection |
| Number of associated cross section. | |
| IntArray | bodyLoadArray |
| Array containing indexes of loads (body loads and boundary loads are kept separately), that apply on receiver. | |
| IntArray | boundaryLoadArray |
| std::vector< std ::unique_ptr< IntegrationRule > > | integrationRulesArray |
| List of integration rules of receiver (each integration rule contains associated integration points also). | |
| FloatMatrix | elemLocalCS |
| Transformation material matrix, used in orthotropic and anisotropic materials, global->local transformation. | |
| int | activityTimeFunction |
| Element activity time function. If defined, nonzero value indicates active receiver, zero value inactive element. | |
| int | globalNumber |
| In parallel mode, globalNumber contains globally unique DoFManager number. | |
| int | numberOfGaussPoints |
| Number of integration points as specified by nip. | |
| elementParallelMode | parallel_mode |
| Determines the parallel mode of the element. | |
| IntArray | partitions |
| List of partition sharing the shared element or remote partition containing remote element counterpart. | |
| IntArray | globalEdgeIDs |
| List of global element edge IDs. | |
| IntArray | globalSurfaceIDs |
| List of global element surface IDs. | |
Protected Attributes inherited from oofem::FEMComponent | |
| int | number |
| Component number. | |
| Domain * | domain |
| Link to domain object, useful for communicating with other FEM components. | |
Static Protected Attributes | |
| static ParamKey | IPK_SUPGElement_bsides |
| static ParamKey | IPK_SUPGElement_bcodes |
Additional Inherited Members | |
Static Public Attributes inherited from oofem::Element | |
| static ParamKey | IPK_Element_mat |
| static ParamKey | IPK_Element_crosssect |
| static ParamKey | IPK_Element_nodes |
| static ParamKey | IPK_Element_bodyload |
| static ParamKey | IPK_Element_boundaryload |
| static ParamKey | IPK_Element_lcs |
| static ParamKey | IPK_Element_partitions |
| static ParamKey | IPK_Element_remote |
| static ParamKey | IPK_Element_activityTimeFunction |
| static ParamKey | IPK_Element_nip |
General stabilized SUPG/PSPG element for CFD analysis.
Definition at line 53 of file supgelement.h.
| oofem::SUPGElement::SUPGElement | ( | int | n, |
| Domain * | aDomain | ||
| ) |
Definition at line 61 of file supgelement.C.
|
overridevirtual |
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.
Reimplemented from oofem::Element.
Reimplemented in oofem::SUPGElement2, oofem::TR1_2D_SUPG, and oofem::TR21_2D_SUPG.
Definition at line 327 of file supgelement.C.
|
pure virtual |
Computes acceleration terms (generalized mass matrix with stabilization terms) for momentum balance equations(s).
Implemented in oofem::SUPGElement2, oofem::TR1_2D_SUPG, oofem::TR1_2D_SUPG2, oofem::TR1_2D_SUPG2_AXI, and oofem::TR1_2D_SUPG_AXI.
Referenced by oofem::SUPGTangentAssembler::matrixFromElement(), and oofem::SUPGInternalForceAssembler::vectorFromElement().
|
pure virtual |
Computes acceleration terms for mass conservation equation.
Implemented in oofem::SUPGElement2, oofem::TR1_2D_SUPG, oofem::TR1_2D_SUPG2, oofem::TR1_2D_SUPG2_AXI, and oofem::TR1_2D_SUPG_AXI.
Referenced by oofem::SUPGTangentAssembler::matrixFromElement(), and oofem::SUPGInternalForceAssembler::vectorFromElement().
|
pure virtual |
Computes the derivative of advection terms for momentum balance equations(s) with respect to nodal velocities.
Implemented in oofem::SUPGElement2, oofem::TR1_2D_SUPG, oofem::TR1_2D_SUPG2, oofem::TR1_2D_SUPG2_AXI, and oofem::TR1_2D_SUPG_AXI.
Referenced by oofem::SUPGTangentAssembler::matrixFromElement().
|
pure virtual |
Computes the derivative of advection terms for mass conservation equation with respect to nodal velocities.
Implemented in oofem::SUPGElement2, oofem::TR1_2D_SUPG, oofem::TR1_2D_SUPG2, oofem::TR1_2D_SUPG2_AXI, and oofem::TR1_2D_SUPG_AXI.
Referenced by oofem::SUPGTangentAssembler::matrixFromElement().
|
pure virtual |
Computes nonlinear advection terms for momentum balance equations(s).
Implemented in oofem::SUPGElement2, oofem::TR1_2D_SUPG, oofem::TR1_2D_SUPG2, oofem::TR1_2D_SUPG2_AXI, and oofem::TR1_2D_SUPG_AXI.
Referenced by oofem::SUPGInternalForceAssembler::vectorFromElement().
|
pure virtual |
Computes advection terms for mass conservation equation.
Implemented in oofem::SUPGElement2, oofem::TR1_2D_SUPG, oofem::TR1_2D_SUPG2, oofem::TR1_2D_SUPG2_AXI, and oofem::TR1_2D_SUPG_AXI.
Referenced by oofem::SUPGInternalForceAssembler::vectorFromElement().
|
virtual |
Computes Lhs terms due to boundary conditions - pressure.
Reimplemented in oofem::TR1_2D_SUPG_AXI.
Definition at line 276 of file supgelement.C.
References oofem::FloatMatrix::add(), oofem::FloatMatrix::clear(), oofem::GeneralBoundaryCondition::giveType(), and oofem::OutFlowBC.
Referenced by oofem::SUPGTangentAssembler::matrixFromElement(), and oofem::SUPGInternalForceAssembler::vectorFromElement().
|
virtual |
Computes Lhs terms due to boundary conditions - pressure.
Definition at line 304 of file supgelement.C.
References oofem::FloatMatrix::add(), oofem::BodyLoadBGT, oofem::FloatMatrix::clear(), oofem::GeneralBoundaryCondition::giveBCGeoType(), oofem::GeneralBoundaryCondition::giveBCValType(), and oofem::ReinforceBVT.
Referenced by oofem::SUPGTangentAssembler::matrixFromElement(), and oofem::SUPGInternalForceAssembler::vectorFromElement().
|
virtual |
Computes Lhs terms due to boundary conditions - velocity.
Reimplemented in oofem::TR1_2D_SUPG_AXI.
Definition at line 231 of file supgelement.C.
References oofem::FloatMatrix::add(), oofem::BodyLoadBGT, oofem::FloatMatrix::clear(), oofem::GeneralBoundaryCondition::giveBCGeoType(), oofem::GeneralBoundaryCondition::giveBCValType(), oofem::GeneralBoundaryCondition::giveType(), oofem::PenetrationWithResistance, oofem::ReinforceBVT, and oofem::SlipWithFriction.
Referenced by oofem::SUPGTangentAssembler::matrixFromElement(), and oofem::SUPGInternalForceAssembler::vectorFromElement().
|
pure virtual |
Computes Rhs terms due to boundary conditions.
Implemented in oofem::SUPGElement2, oofem::TR1_2D_SUPG, oofem::TR1_2D_SUPG2, oofem::TR1_2D_SUPG2_AXI, and oofem::TR1_2D_SUPG_AXI.
|
pure virtual |
Computes Rhs terms due to boundary conditions.
Implemented in oofem::SUPGElement2, oofem::TR1_2D_SUPG, oofem::TR1_2D_SUPG2, oofem::TR1_2D_SUPG2_AXI, and oofem::TR1_2D_SUPG_AXI.
|
pure virtual |
Computes the critical time increment.
Implemented in oofem::Quad10_2D_SUPG, oofem::Tet1_3D_SUPG, oofem::TR1_2D_SUPG, oofem::TR1_2D_SUPG2, oofem::TR1_2D_SUPG2_AXI, oofem::TR21_2D_SUPG, and oofem::SUPGElement2.
Referenced by oofem::SUPG::giveNextStep().
|
protectedpure virtual |
Implemented in oofem::SUPGElement2, oofem::TR1_2D_SUPG, and oofem::TR1_2D_SUPG_AXI.
|
protectedpure virtual |
|
pure virtual |
Computes the derivative of diffusion terms for momentum balance equations(s) with respect to nodal velocities.
Implemented in oofem::SUPGElement2, oofem::TR1_2D_SUPG, oofem::TR1_2D_SUPG2, oofem::TR1_2D_SUPG2_AXI, and oofem::TR1_2D_SUPG_AXI.
Referenced by oofem::SUPGTangentAssembler::matrixFromElement().
|
pure virtual |
Computes diffusion derivative terms for mass conservation equation.
Implemented in oofem::SUPGElement2, oofem::TR1_2D_SUPG, oofem::TR1_2D_SUPG2, oofem::TR1_2D_SUPG2_AXI, and oofem::TR1_2D_SUPG_AXI.
Referenced by oofem::SUPGTangentAssembler::matrixFromElement().
|
pure virtual |
Computes diffusion terms for momentum balance equations(s).
Implemented in oofem::SUPGElement2, oofem::TR1_2D_SUPG, oofem::TR1_2D_SUPG2, oofem::TR1_2D_SUPG2_AXI, and oofem::TR1_2D_SUPG_AXI.
Referenced by oofem::SUPGInternalForceAssembler::vectorFromElement().
|
pure virtual |
Computes diffusion terms for mass conservation equation.
Implemented in oofem::SUPGElement2, oofem::TR1_2D_SUPG, oofem::TR1_2D_SUPG2, oofem::TR1_2D_SUPG2_AXI, and oofem::TR1_2D_SUPG_AXI.
Referenced by oofem::SUPGInternalForceAssembler::vectorFromElement().
|
inlinevirtual |
Reimplemented in oofem::TR1_2D_SUPG.
Definition at line 183 of file supgelement.h.
References oofem::FloatMatrix::clear(), and OOFEM_WARNING.
|
inlinevirtual |
Reimplemented in oofem::TR1_2D_SUPG.
Definition at line 187 of file supgelement.h.
References oofem::FloatMatrix::clear(), and OOFEM_WARNING.
|
pure virtual |
Computes the linear advection term for mass conservation equation.
Implemented in oofem::SUPGElement2, oofem::TR1_2D_SUPG, oofem::TR1_2D_SUPG2, oofem::TR1_2D_SUPG2_AXI, and oofem::TR1_2D_SUPG_AXI.
Referenced by oofem::SUPGTangentAssembler::matrixFromElement(), and oofem::SUPGInternalForceAssembler::vectorFromElement().
|
pure virtual |
Computes SLIC stabilization term for momentum balance equation(s).
Implemented in oofem::SUPGElement2, oofem::TR1_2D_SUPG, oofem::TR1_2D_SUPG2, oofem::TR1_2D_SUPG2_AXI, and oofem::TR1_2D_SUPG_AXI.
Referenced by oofem::SUPGTangentAssembler::matrixFromElement(), and oofem::SUPGInternalForceAssembler::vectorFromElement().
|
inlinevirtual |
Computes Lhs contribution due to outflow BC.
Reimplemented in oofem::TR1_2D_SUPG, and oofem::TR1_2D_SUPG_AXI.
Definition at line 177 of file supgelement.h.
References oofem::FloatMatrix::clear(), and OOFEM_WARNING.
|
inlinevirtual |
Computes Lhs contribution due to applied Penetration bc.
Reimplemented in oofem::TR1_2D_SUPG, and oofem::TR1_2D_SUPG_AXI.
Definition at line 170 of file supgelement.h.
References oofem::FloatMatrix::clear(), and OOFEM_WARNING.
|
pure virtual |
Computes pressure terms for momentum balance equations(s).
Implemented in oofem::SUPGElement2, oofem::TR1_2D_SUPG, oofem::TR1_2D_SUPG2, oofem::TR1_2D_SUPG2_AXI, and oofem::TR1_2D_SUPG_AXI.
Referenced by oofem::SUPGTangentAssembler::matrixFromElement(), and oofem::SUPGInternalForceAssembler::vectorFromElement().
|
pure virtual |
Computes pressure terms for mass conservation equation.
Implemented in oofem::SUPGElement2, oofem::TR1_2D_SUPG, oofem::TR1_2D_SUPG2, oofem::TR1_2D_SUPG2_AXI, and oofem::TR1_2D_SUPG_AXI.
Referenced by oofem::SUPGTangentAssembler::matrixFromElement(), and oofem::SUPGInternalForceAssembler::vectorFromElement().
|
inlinevirtual |
Computes Lhs term due to applied slip with friction bc.
Reimplemented in oofem::TR1_2D_SUPG, and oofem::TR1_2D_SUPG_AXI.
Definition at line 163 of file supgelement.h.
References oofem::FloatMatrix::clear(), and OOFEM_WARNING.
|
protectedpure virtual |
|
overridevirtual |
Computes characteristic matrix of receiver of requested type in given time step.
| 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::SUPGElement2.
Definition at line 95 of file supgelement.C.
References oofem::FloatMatrix::assemble(), oofem::FloatMatrix::isNotEmpty(), OOFEM_ERROR, oofem::FloatMatrix::resize(), oofem::FloatMatrix::times(), oofem::VST_Density, oofem::VST_Velocity, and oofem::FloatMatrix::zero().
|
overridevirtual |
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.
| 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::SUPGElement2.
Definition at line 155 of file supgelement.C.
References oofem::FloatArray::assemble(), oofem::FloatArray::beProductOf(), oofem::FloatArray::clear(), OOFEM_ERROR, oofem::FloatArray::resize(), and oofem::FloatArray::zero().
|
overridevirtual |
Setups the input record string of receiver.
| input | Dynamic input record to be filled by receiver. |
Reimplemented from oofem::Element.
Reimplemented in oofem::SUPGElement2, oofem::TR1_2D_SUPG, oofem::TR1_2D_SUPG2, and oofem::TR1_2D_SUPG2_AXI.
Definition at line 84 of file supgelement.C.
References oofem::DynamicInputRecord::setField().
|
overridevirtual |
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.
| 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. |
Reimplemented from oofem::Element.
Reimplemented in oofem::SUPGElement2, oofem::TR1_2D_SUPG, oofem::TR1_2D_SUPG2, oofem::TR1_2D_SUPG2_AXI, and oofem::TR21_2D_SUPG.
Definition at line 355 of file supgelement.C.
References oofem::FloatArray::at(), oofem::DofManager::end(), oofem::DofManager::findDofWithDofId(), and oofem::FloatArray::resize().
|
inlinevirtual |
Reimplemented in oofem::Quad10_2D_SUPG, oofem::TR1_2D_SUPG, and oofem::TR21_2D_SUPG.
Definition at line 206 of file supgelement.h.
Referenced by oofem::SUPGTangentAssembler::matrixFromElement(), and oofem::SUPGInternalForceAssembler::vectorFromElement().
|
inlinevirtual |
Reimplemented in oofem::Quad10_2D_SUPG, oofem::TR1_2D_SUPG, and oofem::TR21_2D_SUPG.
Definition at line 205 of file supgelement.h.
Referenced by oofem::SUPGTangentAssembler::matrixFromElement(), and oofem::SUPGInternalForceAssembler::vectorFromElement().
|
overridevirtual |
Reimplemented from oofem::Element.
Reimplemented in oofem::SUPGElement2, oofem::TR1_2D_SUPG, oofem::TR1_2D_SUPG2, and oofem::TR1_2D_SUPG2_AXI.
Definition at line 67 of file supgelement.C.
References OOFEM_ERROR, and PM_UPDATE_PARAMETER_AND_REPORT.
|
inlinevirtual |
Reimplemented in oofem::SUPGElement2, oofem::TR1_2D_SUPG2, and oofem::TR1_2D_SUPG2_AXI.
Definition at line 79 of file supgelement.h.
Referenced by oofem::SUPG::updateElementsForNewInterfacePosition().
|
overridevirtual |
Updates element state after equilibrium in time step has been reached.
Default implementation updates all integration rules defined by integrationRulesArray member variable. Doing this, all integration points and their material statuses are updated also. All temporary history variables, which now describe equilibrium state are copied into equilibrium ones. The existing internal state is used for update.
| tStep | Time step for newly reached state. |
Reimplemented from oofem::Element.
Reimplemented in oofem::SUPGElement2.
Definition at line 339 of file supgelement.C.
Referenced by oofem::SUPG::applyIC().
|
inlineoverridevirtual |
Updates the stabilization coefficients used for CBS and SUPG algorithms.
| tStep | Active time step. |
Reimplemented from oofem::FMElement.
Reimplemented in oofem::Tet1_3D_SUPG, oofem::TR1_2D_SUPG, oofem::TR1_2D_SUPG2, oofem::TR1_2D_SUPG2_AXI, oofem::TR1_2D_SUPG_AXI, and oofem::TR21_2D_SUPG.
Definition at line 78 of file supgelement.h.
|
protected |
Boundary sides codes.
Definition at line 59 of file supgelement.h.
|
protected |
Array of boundary sides.
Definition at line 57 of file supgelement.h.
|
staticprotected |
Definition at line 68 of file supgelement.h.
|
staticprotected |
Definition at line 67 of file supgelement.h.
|
protected |
Definition at line 65 of file supgelement.h.
|
protected |
Definition at line 65 of file supgelement.h.
|
protected |
Stabilization coefficients, updated for each solution step in updateStabilizationCoeffs()
Definition at line 65 of file supgelement.h.