OOFEM 3.0
Loading...
Searching...
No Matches
oofem::AbaqusUserElement Class Reference

#include <AbaqusUserElement.h>

Inheritance diagram for oofem::AbaqusUserElement:
Collaboration diagram for oofem::AbaqusUserElement:

Public Member Functions

 AbaqusUserElement (int n, Domain *d)
 Constructor.
virtual ~AbaqusUserElement ()
 Destructor.
void computeConsistentMassMatrix (FloatMatrix &answer, TimeStep *tStep, double &mass, const double *ipDensity=NULL) override
void computeStiffnessMatrix (FloatMatrix &answer, MatResponseMode rMode, TimeStep *tStep) override
void giveInternalForcesVector (FloatArray &answer, TimeStep *tStep, int useUpdatedGpRecord=0) override
virtual void giveInternalForcesVector (FloatArray &answer, TimeStep *tStep, FloatArray &U, FloatMatrix &DU, int useUpdatedGpRecord)
int computeNumberOfDofs () override
void giveDofManDofIDMask (int inode, IntArray &answer) const override
void computeField (ValueModeType mode, TimeStep *tStep, const FloatArray &lcoords, FloatArray &answer) override
void updateYourself (TimeStep *tStep) override
void updateInternalState (TimeStep *tStep) override
bool hasTangent () const
virtual void letTempTangentBe (FloatMatrix &src)
virtual FloatMatrixletTempRhsBe (FloatMatrix &src)
virtual FloatArrayletTempSvarsBe (FloatArray &src)
virtual const FloatArraygiveStateVector () const
virtual const FloatArraygiveTempStateVector () const
virtual const FloatMatrixgiveTempTangent ()
void initializeFrom (InputRecord &ir, int priority) override
void giveInputRecord (DynamicInputRecord &input) override
void postInitialize () override
 Performs post initialization steps. Called after all components are created and initialized.
const char * giveClassName () const override
const char * giveInputRecordName () const override
integrationDomain giveIntegrationDomain () const override
Element_Geometry_Type giveGeometryType () const override
Public Member Functions inherited from oofem::StructuralElement
 StructuralElement (int n, Domain *d)
virtual ~StructuralElement ()
 Destructor.
void giveCharacteristicMatrix (FloatMatrix &answer, CharType, TimeStep *tStep) override
void giveCharacteristicVector (FloatArray &answer, CharType type, ValueModeType mode, TimeStep *tStep) override
virtual void computeMassMatrix (FloatMatrix &answer, TimeStep *tStep)
virtual void computeLumpedMassMatrix (FloatMatrix &answer, TimeStep *tStep)
virtual void giveMassMtrxIntegrationgMask (IntArray &answer)
void computeStiffnessMatrix_withIRulesAsSubcells (FloatMatrix &answer, MatResponseMode rMode, TimeStep *tStep)
virtual void computeInitialStressMatrix (FloatMatrix &answer, TimeStep *tStep)
virtual void computeLumpedInitialStressMatrix (FloatMatrix &answer, TimeStep *tStep)
void computeField (ValueModeType mode, TimeStep *tStep, const FloatArray &lcoords, FloatArray &answer) override
virtual void giveInternalForcesVector_withIRulesAsSubcells (FloatArray &answer, TimeStep *tStep, int useUpdatedGpRecord=0)
virtual void computeStrainVector (FloatArray &answer, GaussPoint *gp, TimeStep *tStep)
int giveIPValue (FloatArray &answer, GaussPoint *gp, InternalStateType type, TimeStep *tStep) override
virtual void computeResultingIPTemperatureAt (FloatArray &answer, TimeStep *tStep, GaussPoint *gp, ValueModeType mode)
virtual void computeResultingIPEigenstrainAt (FloatArray &answer, TimeStep *tStep, GaussPoint *gp, ValueModeType mode)
virtual void computeNmatrixAt (const FloatArray &iLocCoord, FloatMatrix &answer)
void updateBeforeNonlocalAverage (TimeStep *tStep) override
virtual void giveNonlocalLocationArray (IntArray &locationArray, const UnknownNumberingScheme &us)
virtual void addNonlocalStiffnessContributions (SparseMtrx &dest, const UnknownNumberingScheme &s, TimeStep *tStep)
int adaptiveUpdate (TimeStep *tStep) override
void updateInternalState (TimeStep *tStep) override
void updateYourself (TimeStep *tStep) override
int checkConsistency () override
void giveInputRecord (DynamicInputRecord &input) override
const char * giveClassName () const override
int giveInternalStateAtNode (FloatArray &answer, InternalStateType type, InternalStateMode mode, int node, TimeStep *tStep) override
void showSparseMtrxStructure (CharType mtrx, oofegGraphicContext &gc, TimeStep *tStep) override
 Shows sparse structure.
void showExtendedSparseMtrxStructure (CharType mtrx, oofegGraphicContext &gc, TimeStep *tStep) override
 Shows extended sparse structure (for example, due to nonlocal interactions for tangent stiffness).
void computeLoadVector (FloatArray &answer, BodyLoad *load, CharType type, ValueModeType mode, TimeStep *tStep) override
void computeBoundarySurfaceLoadVector (FloatArray &answer, BoundaryLoad *load, int boundary, CharType type, ValueModeType mode, TimeStep *tStep, bool global=true) override
void computeBoundaryEdgeLoadVector (FloatArray &answer, BoundaryLoad *load, int boundary, CharType type, ValueModeType mode, TimeStep *tStep, bool global=true) override
virtual void computeEdgeNMatrix (FloatMatrix &answer, int boundaryID, const FloatArray &lcoords)
 computes edge interpolation matrix
virtual void computeSurfaceNMatrix (FloatMatrix &answer, int boundaryID, const FloatArray &lcoords)
StructuralCrossSectiongiveStructuralCrossSection ()
 Helper function which returns the structural cross-section for the element.
virtual void createMaterialStatus ()
Public Member Functions inherited from oofem::Element
 Element (int n, Domain *aDomain)
 Element (const Element &src)=delete
Elementoperator= (const Element &src)=delete
virtual ~Element ()
 Virtual destructor.
void giveLocationArray (IntArray &locationArray, const UnknownNumberingScheme &s, IntArray *dofIds=NULL) const
void giveLocationArray (IntArray &locationArray, const IntArray &dofIDMask, const UnknownNumberingScheme &s, IntArray *dofIds=NULL) const
virtual void giveBoundaryLocationArray (IntArray &locationArray, const IntArray &bNodes, const UnknownNumberingScheme &s, IntArray *dofIds=NULL)
virtual void giveBoundaryLocationArray (IntArray &locationArray, const IntArray &bNodes, const IntArray &dofIDMask, const UnknownNumberingScheme &s, IntArray *dofIds=NULL)
virtual int giveNumberOfDofs ()
virtual int giveNumberOfInternalDofManagers () const
virtual DofManagergiveInternalDofManager (int i) const
virtual void setInternalDofManager (int num, std::unique_ptr< DofManager > dm)
virtual double giveCharacteristicValue (CharType type, TimeStep *tStep)
virtual void computeTangentFromSurfaceLoad (FloatMatrix &answer, BoundaryLoad *load, int boundary, MatResponseMode rmode, TimeStep *tStep)
virtual void computeTangentFromEdgeLoad (FloatMatrix &answer, BoundaryLoad *load, int boundary, MatResponseMode rmode, TimeStep *tStep)
const IntArraygiveBodyLoadList () const
const IntArraygiveBoundaryLoadList () const
void computeVectorOf (ValueModeType u, TimeStep *tStep, FloatArray &answer)
void computeVectorOf (const IntArray &dofIDMask, ValueModeType u, TimeStep *tStep, FloatArray &answer, bool padding=false)
void computeBoundaryVectorOf (const IntArray &bNodes, const IntArray &dofIDMask, ValueModeType u, TimeStep *tStep, FloatArray &answer, bool padding=false)
void computeVectorOf (PrimaryField &field, const IntArray &dofIDMask, ValueModeType u, TimeStep *tStep, FloatArray &answer, bool padding=false)
void computeVectorOfPrescribed (ValueModeType u, TimeStep *tStep, FloatArray &answer)
void computeVectorOfPrescribed (const IntArray &dofIDMask, ValueModeType type, TimeStep *tStep, FloatArray &answer)
virtual int computeNumberOfGlobalDofs ()
int computeNumberOfPrimaryMasterDofs ()
virtual bool computeGtoLRotationMatrix (FloatMatrix &answer)
virtual bool giveRotationMatrix (FloatMatrix &answer)
virtual bool computeDofTransformationMatrix (FloatMatrix &answer, const IntArray &nodes, bool includeInternal)
virtual void giveInternalDofManDofIDMask (int inode, IntArray &answer) const
virtual void giveElementDofIDMask (IntArray &answer) const
virtual double computeVolumeAround (GaussPoint *gp)
virtual double computeVolumeAreaOrLength ()
 Computes the volume, area or length of the element depending on its spatial dimension.
double computeMeanSize ()
virtual double computeVolume ()
virtual double computeArea ()
virtual double computeLength ()
virtual IntArray giveBoundaryEdgeNodes (int boundary, bool includeHierarchical=false) const
virtual IntArray giveBoundarySurfaceNodes (int boundary, bool includeHierarchical=false) const
virtual IntArray giveBoundaryNodes (int boundary) const
virtual std::unique_ptr< IntegrationRulegiveBoundaryEdgeIntegrationRule (int order, int boundary)
virtual std::unique_ptr< IntegrationRulegiveBoundarySurfaceIntegrationRule (int order, int boundary)
int giveDofManagerNumber (int i) const
const IntArraygiveDofManArray () const
void addDofManager (DofManager *dMan)
DofManagergiveDofManager (int i) const
NodegiveNode (int i) const
virtual ElementSidegiveSide (int i) const
virtual FEInterpolationgiveInterpolation () const
virtual FEInterpolationgiveInterpolation (DofIDItem id) const
virtual const FEInterpolationgetGeometryInterpolation () const
virtual MaterialgiveMaterial ()
int giveMaterialNumber () const
CrossSectiongiveCrossSection ()
int getActivityTimeFunctionNumber ()
void setActivityTimeFunctionNumber (int funcIndx)
void setMaterial (int matIndx)
virtual void setCrossSection (int csIndx)
virtual int giveNumberOfDofManagers () const
void setNumberOfDofManagers (int i)
 Sets number of element dof managers.
virtual int giveNumberOfNodes () const
void setDofManagers (const IntArray &dmans)
void setDofManager (int id, int dm)
void setBodyLoads (const IntArray &bodyLoads)
void setIntegrationRules (std ::vector< std ::unique_ptr< IntegrationRule > > irlist)
virtual MaterialMode giveMaterialMode ()
virtual int giveIntegrationRuleLocalCodeNumbers (IntArray &answer, IntegrationRule &ie)
int giveRegionNumber ()
virtual void initializeYourself (TimeStep *timeStepWhenICApply)
virtual bool isActivated (TimeStep *tStep)
virtual bool isCast (TimeStep *tStep)
virtual void initForNewStep ()
virtual Element_Geometry_Type giveEdgeGeometryType (int id) const
 Returns the receiver edge geometry type.
virtual Element_Geometry_Type giveSurfaceGeometryType (int id) const
 Returns the receiver surface geometry type.
virtual int giveSpatialDimension ()
virtual int giveNumberOfBoundarySides ()
 Returns number of boundaries (entities of element_dimension-1: points, edges, surfaces).
virtual int giveNumberOfEdges () const
virtual int giveNumberOfSurfaces () const
virtual int giveDefaultIntegrationRule () const
virtual IntegrationRulegiveDefaultIntegrationRulePtr ()
int giveNumberOfIntegrationRules ()
virtual IntegrationRulegiveIntegrationRule (int i)
std::vector< std ::unique_ptr< IntegrationRule > > & giveIntegrationRulesArray ()
virtual int testElementExtension (ElementExtension ext)
int giveGlobalIPValue (FloatArray &answer, GaussPoint *gp, InternalStateType type, TimeStep *tStep)
virtual double giveLengthInDir (const FloatArray &normalToCrackPlane)
virtual double giveCharacteristicLength (const FloatArray &normalToCrackPlane)
double giveCharacteristicLengthForPlaneElements (const FloatArray &normalToCrackPlane)
double giveCharacteristicLengthForAxisymmElements (const FloatArray &normalToCrackPlane)
virtual double giveCharacteristicSize (GaussPoint *gp, FloatArray &normalToCrackPlane, ElementCharSizeMethod method)
virtual double giveParentElSize () const
virtual int computeGlobalCoordinates (FloatArray &answer, const FloatArray &lcoords)
virtual bool computeLocalCoordinates (FloatArray &answer, const FloatArray &gcoords)
virtual int giveLocalCoordinateSystem (FloatMatrix &answer)
virtual void giveLocalCoordinateSystemVector (InternalStateType isttype, FloatArray &answer)
virtual void computeMidPlaneNormal (FloatArray &answer, const GaussPoint *gp)
virtual int adaptiveMap (Domain *oldd, TimeStep *tStep)
virtual int mapStateVariables (Domain &iOldDom, const TimeStep &iTStep)
virtual int adaptiveFinish (TimeStep *tStep)
void updateLocalNumbering (EntityRenumberingFunctor &f) override
template<class T>
void ipEvaluator (T *src, void(T ::*f)(GaussPoint *gp))
 Integration point evaluator, loops over receiver IP's and calls given function (passed as f parameter) on them. The IP is parameter to function f.
template<class T, class S>
void ipEvaluator (T *src, void(T ::*f)(GaussPoint *, S &), S &_val)
 Integration point evaluator, loops over receiver IP's and calls given function (passed as f parameter) on them. The IP is parameter to function f as well as additional array.
virtual void drawYourself (oofegGraphicContext &gc, TimeStep *tStep)
virtual void drawAnnotation (oofegGraphicContext &gc, TimeStep *tStep)
virtual void 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)
int giveLabel () const
int giveGlobalNumber () const
void setGlobalNumber (int num)
elementParallelMode giveParallelMode () const
void setParallelMode (elementParallelMode _mode)
 Sets parallel mode of element.
virtual elementParallelMode giveKnotSpanParallelMode (int) const
int packUnknowns (DataStream &buff, TimeStep *tStep)
int unpackAndUpdateUnknowns (DataStream &buff, TimeStep *tStep)
int estimatePackSize (DataStream &buff)
const IntArraygivePartitionList () const
void setPartitionList (IntArray &pl)
virtual double predictRelativeComputationalCost ()
virtual double giveRelativeSelfComputationalCost ()
virtual double predictRelativeRedistributionCost ()
IntArraygiveBodyLoadArray ()
 Returns array containing load numbers of loads acting on element.
IntArraygiveBoundaryLoadArray ()
 Returns array containing load numbers of boundary loads acting on element.
void initializeFrom (InputRecord &ir, int priority) override
void initializeFinish () override
void postInitialize () override
 Performs post initialization steps.
void saveContext (DataStream &stream, ContextMode mode) override
void restoreContext (DataStream &stream, ContextMode mode) override
void printOutputAt (FILE *file, TimeStep *tStep) override
virtual const IntArray giveLocation ()
virtual void recalculateCoordinates (int nodeNumber, FloatArray &coords)
void setSharedEdgeID (int iedge, int globalID)
void setSharedSurfaceID (int isurf, int globalID)
const IntArraygiveSharedEdgeIDs () const
const IntArraygiveSharedSurfaceIDs () const
Public Member Functions inherited from oofem::FEMComponent
 FEMComponent (int n, Domain *d)
virtual ~FEMComponent ()=default
 Virtual destructor.
DomaingiveDomain () const
virtual void setDomain (Domain *d)
int giveNumber () const
void setNumber (int num)
virtual void initializeFrom (InputRecord &ir)
virtual void printYourself ()
 Prints receiver state on stdout. Useful for debugging.
virtual InterfacegiveInterface (InterfaceType t)
std::string errorInfo (const char *func) const
 Returns string for prepending output (used by error reporting macros).

Protected Member Functions

void computeStressVector (FloatArray &answer, const FloatArray &strain, GaussPoint *gp, TimeStep *tStep) override
void computeConstitutiveMatrixAt (FloatMatrix &answer, MatResponseMode rMode, GaussPoint *gp, TimeStep *tStep) override
void computeBmatrixAt (GaussPoint *, FloatMatrix &, int=1, int=ALL_STRAINS) override
int giveIPValue (FloatArray &answer, GaussPoint *gp, InternalStateType type, TimeStep *tStep) override
virtual void computePointLoadVectorAt (FloatArray &answer, Load *load, TimeStep *tStep, ValueModeType mode, bool global=true)
virtual void giveEdgeDofMapping (IntArray &answer, int iEdge) const
virtual void giveSurfaceDofMapping (IntArray &answer, int iSurf) const
virtual double computeEdgeVolumeAround (GaussPoint *gp, int iEdge)
virtual double computeSurfaceVolumeAround (GaussPoint *gp, int iSurf)
virtual int computeLoadGToLRotationMtrx (FloatMatrix &answer)
virtual int computeLoadLEToLRotationMatrix (FloatMatrix &answer, int iEdge, GaussPoint *gp)
virtual int computeLoadLSToLRotationMatrix (FloatMatrix &answer, int iSurf, GaussPoint *gp)
virtual int giveNumberOfIPForMassMtrxIntegration ()
void condense (FloatMatrix *stiff, FloatMatrix *mass, FloatArray *load, IntArray *what)
virtual void setupIRForMassMtrxIntegration (IntegrationRule &iRule)
virtual void computeBodyLoadVectorAt (FloatArray &answer, Load *load, TimeStep *tStep, ValueModeType mode)
virtual void computeGaussPoints ()

Private Attributes

void * uelobj
 Dynamically loaded uel.
int nCoords
 Coord transf.
IntArray dofs
FloatMatrix coords
int numSvars
 Number of status variables.
FloatArray props
 Element properties.
FloatArray svars
 Status variables.
FloatArray tempSvars
int jtype
 Element type.
int nrhs = 2
 Element rhs.
FloatMatrix rhs
FloatMatrix tempRHS
FloatMatrix amatrx
 Element amatrx.
FloatMatrix tempAmatrx
int ndofel = 0
 ndofel
int mcrd = 2
 mcrd
int mlvarx = 1
 mlvarx
FloatArray U
 Inputs to element routines. Velocity and Acceleration currently ignored.
FloatArray V
FloatArray A
FloatMatrix DU
IntArray lFlags
 LFlags.
int kinc = 1
 kinc, kstep
int kstep = 1
int npredef = 1
 predef
FloatArray predef
FloatArray energy
 energy
int ndLoad = 0
 loads
int mdLoad = 0
FloatMatrix adlmag
FloatMatrix ddlmag
int * jdltype = nullptr
double params [3]
 params
IntArray jprops
 jprops
bool hasTangentFlag
 Keeps track of whether the tangent has been obtained already.
void(* uel )(double *rhs, double *amatrx, double *svars, double energy[8], int *ndofel, int *nrhs, int *nsvars, double *props, int *nprops, double *coords, int *mcrd, int *nnode, double *u, double *du, double *v, double *a, int *jtype, double time[2], double *dtime, int *kstep, int *kinc, int *jelem, double params[3], int *ndload, int *jdltyp, double *adlmag, double *predef, int *npredef, int *lflags, int *mvarx, double *ddlmag, int *mdload, double *pnewdt, int *jprops, int *njprop, double *period)
 Pointer to the dynamically loaded uel-function (translated to C).
std::string filename
 File containing the uel function.

Static Private Attributes

static ParamKey IPK_AbaqusUserElement_userElement
static ParamKey IPK_AbaqusUserElement_numcoords
static ParamKey IPK_AbaqusUserElement_dofs
static ParamKey IPK_AbaqusUserElement_numsvars
static ParamKey IPK_AbaqusUserElement_properties
static ParamKey IPK_AbaqusUserElement_type
static ParamKey IPK_AbaqusUserElement_name

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
Protected Attributes inherited from oofem::StructuralElement
std::unique_ptr< FloatArrayinitialDisplacements
 Initial displacement vector, describes the initial nodal displacements when element has been casted.
Protected Attributes inherited from oofem::Element
int numberOfDofMans
 Number of dofmanagers.
IntArray dofManArray
 Array containing dofmanager numbers.
int material
 Number of associated material.
int crossSection
 Number of associated cross section.
IntArray bodyLoadArray
IntArray boundaryLoadArray
std::vector< std ::unique_ptr< IntegrationRule > > integrationRulesArray
FloatMatrix elemLocalCS
 Transformation material matrix, used in orthotropic and anisotropic materials, global->local transformation.
int activityTimeFunction
 Element activity time function. If defined, nonzero value indicates active receiver, zero value inactive element.
int globalNumber
int numberOfGaussPoints
elementParallelMode parallel_mode
 Determines the parallel mode of the element.
IntArray partitions
IntArray globalEdgeIDs
IntArray globalSurfaceIDs
Protected Attributes inherited from oofem::FEMComponent
int number
 Component number.
Domaindomain
 Link to domain object, useful for communicating with other FEM components.

Detailed Description

UEL interface from Abaqus user elements.

The function prototype for UEL is: SUBROUTINE UEL(RHS,AMATRX,SVARS,ENERGY,NDOFEL,NRHS,NSVARS, PROPS,NPROPS,COORDS,MCRD,NNODE,U,DU,V,A,JTYPE,TIME, DTIME,KSTEP,KINC,JELEM,PARAMS,NDLOAD,JDLTYP,ADLMAG, PREDEF,NPREDF,LFLAGS,MLVARX,DDLMAG,MDLOAD,PNEWDT, JPROPS,NJPROP,PERIOD)

DIMENSION RHS(MLVARX,*),AMATRX(NDOFEL,NDOFEL), SVARS(NSVARS),ENERGY(8),PROPS(*),COORDS(MCRD,NNODE), U(NDOFEL),DU(MLVARX,*),V(NDOFEL),A(NDOFEL),TIME(2), PARAMS(3),JDLTYP(MDLOAD,*),ADLMAG(MDLOAD,*), DDLMAG(MDLOAD,*),PREDEF(2,NPREDF,NNODE),LFLAGS(*), JPROPS(*)

Author
Giovanni

Definition at line 76 of file AbaqusUserElement.h.

Constructor & Destructor Documentation

◆ AbaqusUserElement()

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

Constructor.

Definition at line 63 of file AbaqusUserElement.C.

References hasTangentFlag, oofem::StructuralElement::StructuralElement(), uel, and uelobj.

◆ ~AbaqusUserElement()

oofem::AbaqusUserElement::~AbaqusUserElement ( )
virtual

Destructor.

Definition at line 67 of file AbaqusUserElement.C.

References uelobj.

Member Function Documentation

◆ computeBmatrixAt()

void oofem::AbaqusUserElement::computeBmatrixAt ( GaussPoint * gp,
FloatMatrix & answer,
int lowerIndx = 1,
int upperIndx = ALL_STRAINS )
inlineoverrideprotectedvirtual

Computes the geometrical matrix of receiver in given integration point. The product of this matrix (assembled at given integration point) and element displacement vector is element strain vector. If lowerIndx and upperIndx parameters are specified, answer is formed only for strains within this interval. This will affects the size of answer.

Parameters
gpIntegration point for which answer is computed.
answerGeometric matrix of receiver.
lowerIndxIf specified, answer is formed only for strain with index equal and greater than lowerIndx. This parameter has default value 1 (answer is formed from first strain).
upperIndxIf specified, answer is formed only for strain with index less and equal than upperIndx. This parameter has default value ALL_STRAINS (answer is formed for all strains).

Implements oofem::StructuralElement.

Definition at line 228 of file AbaqusUserElement.h.

References ALL_STRAINS, and OOFEM_ERROR.

◆ computeConsistentMassMatrix()

void oofem::AbaqusUserElement::computeConsistentMassMatrix ( FloatMatrix & answer,
TimeStep * tStep,
double & mass,
const double * ipDensity = NULL )
overridevirtual

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
answerMass matrix.
tStepTime step.
massTotal mass of receiver.

Reimplemented from oofem::StructuralElement.

Definition at line 306 of file AbaqusUserElement.C.

References ndofel, oofem::FloatMatrix::resize(), and oofem::FloatMatrix::zero().

◆ computeConstitutiveMatrixAt()

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

Computes constitutive matrix of receiver. Default implementation uses element cross section giveCharMaterialStiffnessMatrix service.

Parameters
answerConstitutive matrix.
rModeMaterial response mode of answer.
gpIntegration point for which constitutive matrix is computed.
tStepTime step.

Implements oofem::StructuralElement.

Definition at line 224 of file AbaqusUserElement.h.

References OOFEM_ERROR.

◆ computeField()

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

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

See also
giveElementDofIDMask The unknown vector should match the element field as specified by the element dof IDs.
Parameters
modeMode (total, increment, etc) of the output
tStepTime step to evaluate at
lcoordsLocal coordinates to evaluate at
answerResults

Reimplemented from oofem::Element.

Definition at line 182 of file AbaqusUserElement.h.

References OOFEM_ERROR.

◆ computeNumberOfDofs()

int oofem::AbaqusUserElement::computeNumberOfDofs ( )
inlineoverridevirtual

Computes or simply returns total number of element's local DOFs. Must be defined by particular element.

Returns
Number of local DOFs of element.

Reimplemented from oofem::Element.

Definition at line 180 of file AbaqusUserElement.h.

◆ computeStiffnessMatrix()

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

Computes numerically stiffness matrix of receiver. Default implementation computes element stiffness using \( K=\int_v B^{\mathrm{T}} D B \mathrm{d}V \) formulae, where \( B \) is element geometric matrix and \( D \) is material stiffness matrix. No geometrical nonlinearity is taken into account. NUmerical integration procedure uses integrationRulesArray for numerical integration. Support for reduced or selected integration is implemented. The individual integration rules are assumed to correspond to different terms from which the overall matrix is assembled.

If there is one integration rule, the full integration of all coefficients is performed. Otherwise, integration is performed using following rules. Each integration rule can specify start and end strain index of strain vector components for which is valid. It is necessary to ensure that these start and end indexes, dividing geometrical matrix into blocks, are not overlapping and that each strain component is included.

Then stiffness matrix is obtained as summation of integrals \( I_{ij}=\int_v B^{\mathrm{T}}_i D_{ij} B_j \mathrm{d}V \) where \( B_i \) is i-th block of geometrical matrix and \( D_{ij} \) is corresponding constitutive sub-matrix. The geometrical matrix is obtained using computeBmatrixAt service and the constitutive matrix is obtained using computeConstitutiveMatrixAt service. The \( I_{ij} \) integral is evaluated using such integration rule, which is valid for i-th or j-th block and has smaller number of integration points.

For higher numerical performance, only one half of stiffness matrix is computed and answer is then symmetrized. Therefore, if element matrix will be generally nonsymmetric, one must specialize this method. Finally, the result is transformed into global coordinate system (or nodal coordinate system, if it is defined).

Parameters
answerComputed stiffness matrix (symmetric).
rModeResponse mode.
tStepTime step.

Reimplemented from oofem::StructuralElement.

Definition at line 187 of file AbaqusUserElement.C.

References DU, giveInternalForcesVector(), giveTempTangent(), hasTangent(), and U.

◆ computeStressVector()

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

Computes the stress vector of receiver at given integration point, at time step tStep. The nature of these stresses depends on the element's type.

Parameters
answerStress vector.
strainStrain vector.
gpIntegration point.
tStepTime step.

Implements oofem::StructuralElement.

Definition at line 221 of file AbaqusUserElement.h.

References OOFEM_ERROR.

◆ giveClassName()

const char * oofem::AbaqusUserElement::giveClassName ( ) const
inlineoverridevirtual
Returns
Class name of the receiver.

Implements oofem::FEMComponent.

Definition at line 215 of file AbaqusUserElement.h.

◆ giveDofManDofIDMask()

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

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
inodeMask is computed for local dofmanager with inode number.
answerMask for node.

Reimplemented from oofem::Element.

Definition at line 182 of file AbaqusUserElement.C.

References dofs.

◆ giveGeometryType()

Element_Geometry_Type oofem::AbaqusUserElement::giveGeometryType ( ) const
inlineoverridevirtual

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.

Implements oofem::Element.

Definition at line 218 of file AbaqusUserElement.h.

◆ giveInputRecord()

void oofem::AbaqusUserElement::giveInputRecord ( DynamicInputRecord & input)
overridevirtual

◆ giveInputRecordName()

const char * oofem::AbaqusUserElement::giveInputRecordName ( ) const
inlineoverridevirtual
Returns
Input record name of the receiver.

Implements oofem::FEMComponent.

Definition at line 216 of file AbaqusUserElement.h.

References _IFT_AbaqusUserElement_Name.

◆ giveIntegrationDomain()

integrationDomain oofem::AbaqusUserElement::giveIntegrationDomain ( ) const
inlineoverridevirtual

Returns integration domain for receiver, used to initialize integration point over receiver volume. Default behavior is taken from the default interpolation.

Returns
Integration domain of element.

Reimplemented from oofem::Element.

Definition at line 217 of file AbaqusUserElement.h.

References oofem::_UnknownIntegrationDomain.

◆ giveInternalForcesVector() [1/2]

◆ giveInternalForcesVector() [2/2]

void oofem::AbaqusUserElement::giveInternalForcesVector ( FloatArray & answer,
TimeStep * tStep,
int useUpdatedGpRecord = 0 )
overridevirtual

Returns equivalent nodal forces vectors. Useful for nonlinear analysis. Default implementation computes result as \( F=\int_v B^{\mathrm{T}} \sigma \mathrm{d}V \), where \( \sigma \) is the real element stress vector obtained using computeStressVector service (if useUpdatedGpRecord=0) or (if useUpdatedGpRecord=1) from integration point status. The geometric matrix is obtained using computeBmatrixAt service. Integration is performed using default integration rule, which should produce always valid results, assuming that strains used for computation of stresses are valid.

Parameters
answerInternal nodal forces vector.
tStepTime step.
useUpdatedGpRecordIf equal to zero, the stresses in integration points are computed (slow but safe), else if nonzero the stresses are taken directly from integration point status (should be derived from StructuralMaterialStatus) (fast, but engineering model must ensure valid status data in each integration point).

Reimplemented from oofem::StructuralElement.

Definition at line 214 of file AbaqusUserElement.C.

References oofem::Element::computeVectorOf(), DU, oofem::Domain::giveClassName(), oofem::FEMComponent::giveDomain(), giveInternalForcesVector(), and U.

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

◆ giveIPValue()

int oofem::AbaqusUserElement::giveIPValue ( FloatArray & answer,
GaussPoint * gp,
InternalStateType type,
TimeStep * tStep )
inlineoverrideprotectedvirtual

Returns the integration point corresponding value in full form.

Parameters
answerContain corresponding integration point value, zero sized if not available.
gpIntegration point to check.
typeDetermines the type of internal variable.
tStepTime step.
Returns
Nonzero if o.k, zero otherwise.
Todo
Which "error type" should be used? Why are there several? I don't see the point of this enum when there could be different function calls just as well (and different IST values)

Reimplemented from oofem::Element.

Definition at line 232 of file AbaqusUserElement.h.

References OOFEM_ERROR.

◆ giveStateVector()

virtual const FloatArray & oofem::AbaqusUserElement::giveStateVector ( ) const
inlinevirtual

Definition at line 200 of file AbaqusUserElement.h.

References svars.

Referenced by giveInternalForcesVector().

◆ giveTempStateVector()

virtual const FloatArray & oofem::AbaqusUserElement::giveTempStateVector ( ) const
inlinevirtual

Definition at line 203 of file AbaqusUserElement.h.

References tempSvars.

◆ giveTempTangent()

virtual const FloatMatrix & oofem::AbaqusUserElement::giveTempTangent ( )
inlinevirtual

Definition at line 206 of file AbaqusUserElement.h.

References tempAmatrx.

Referenced by computeStiffnessMatrix().

◆ hasTangent()

bool oofem::AbaqusUserElement::hasTangent ( ) const
inline

Definition at line 187 of file AbaqusUserElement.h.

References hasTangentFlag.

Referenced by computeStiffnessMatrix().

◆ initializeFrom()

◆ letTempRhsBe()

virtual FloatMatrix & oofem::AbaqusUserElement::letTempRhsBe ( FloatMatrix & src)
inlinevirtual

Definition at line 194 of file AbaqusUserElement.h.

References tempRHS.

Referenced by giveInternalForcesVector().

◆ letTempSvarsBe()

virtual FloatArray & oofem::AbaqusUserElement::letTempSvarsBe ( FloatArray & src)
inlinevirtual

Definition at line 197 of file AbaqusUserElement.h.

References tempSvars.

Referenced by giveInternalForcesVector().

◆ letTempTangentBe()

virtual void oofem::AbaqusUserElement::letTempTangentBe ( FloatMatrix & src)
inlinevirtual

Definition at line 190 of file AbaqusUserElement.h.

References hasTangentFlag, and tempAmatrx.

Referenced by giveInternalForcesVector().

◆ postInitialize()

◆ updateInternalState()

void oofem::AbaqusUserElement::updateInternalState ( TimeStep * tStep)
overridevirtual

Updates element state after equilibrium in time step has been reached. Default implementation updates all integration rules defined by integrationRulesArray member variable. Doing this, all integration points and their material statuses are updated also. All temporary history variables, which now describe equilibrium state are copied into equilibrium ones. The existing internal state is used for update.

Parameters
tStepTime step for newly reached state.
See also
Material::updateYourself
IntegrationRule::updateYourself
GaussPoint::updateYourself
Element::updateInternalState

Reimplemented from oofem::Element.

Definition at line 208 of file AbaqusUserElement.C.

References giveInternalForcesVector().

◆ updateYourself()

void oofem::AbaqusUserElement::updateYourself ( TimeStep * tStep)
overridevirtual

Updates element state after equilibrium in time step has been reached. Default implementation updates all integration rules defined by integrationRulesArray member variable. Doing this, all integration points and their material statuses are updated also. All temporary history variables, which now describe equilibrium state are copied into equilibrium ones. The existing internal state is used for update.

Parameters
tStepTime step for newly reached state.
See also
Material::updateYourself
IntegrationRule::updateYourself
GaussPoint::updateYourself
Element::updateInternalState

Reimplemented from oofem::Element.

Definition at line 199 of file AbaqusUserElement.C.

References amatrx, hasTangentFlag, rhs, svars, tempAmatrx, tempRHS, tempSvars, and oofem::StructuralElement::updateYourself().

Member Data Documentation

◆ A

FloatArray oofem::AbaqusUserElement::A
private

Definition at line 116 of file AbaqusUserElement.h.

Referenced by giveInternalForcesVector(), and postInitialize().

◆ adlmag

FloatMatrix oofem::AbaqusUserElement::adlmag
private

Definition at line 136 of file AbaqusUserElement.h.

Referenced by giveInternalForcesVector().

◆ amatrx

FloatMatrix oofem::AbaqusUserElement::amatrx
private

Element amatrx.

Definition at line 104 of file AbaqusUserElement.h.

Referenced by postInitialize(), and updateYourself().

◆ coords

FloatMatrix oofem::AbaqusUserElement::coords
private

Definition at line 85 of file AbaqusUserElement.h.

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

◆ ddlmag

FloatMatrix oofem::AbaqusUserElement::ddlmag
private

Definition at line 136 of file AbaqusUserElement.h.

Referenced by giveInternalForcesVector().

◆ dofs

IntArray oofem::AbaqusUserElement::dofs
private

Definition at line 84 of file AbaqusUserElement.h.

Referenced by giveDofManDofIDMask(), giveInputRecord(), and initializeFrom().

◆ DU

FloatMatrix oofem::AbaqusUserElement::DU
private

◆ energy

FloatArray oofem::AbaqusUserElement::energy
private

energy

Definition at line 130 of file AbaqusUserElement.h.

Referenced by giveInternalForcesVector(), and postInitialize().

◆ filename

std::string oofem::AbaqusUserElement::filename
private

File containing the uel function.

Definition at line 158 of file AbaqusUserElement.h.

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

◆ hasTangentFlag

bool oofem::AbaqusUserElement::hasTangentFlag
private

Keeps track of whether the tangent has been obtained already.

Definition at line 146 of file AbaqusUserElement.h.

Referenced by AbaqusUserElement(), hasTangent(), letTempTangentBe(), and updateYourself().

◆ IPK_AbaqusUserElement_dofs

ParamKey oofem::AbaqusUserElement::IPK_AbaqusUserElement_dofs
staticprivate

Definition at line 162 of file AbaqusUserElement.h.

Referenced by giveInputRecord(), and initializeFrom().

◆ IPK_AbaqusUserElement_name

ParamKey oofem::AbaqusUserElement::IPK_AbaqusUserElement_name
staticprivate

Definition at line 166 of file AbaqusUserElement.h.

◆ IPK_AbaqusUserElement_numcoords

ParamKey oofem::AbaqusUserElement::IPK_AbaqusUserElement_numcoords
staticprivate

Definition at line 161 of file AbaqusUserElement.h.

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

◆ IPK_AbaqusUserElement_numsvars

ParamKey oofem::AbaqusUserElement::IPK_AbaqusUserElement_numsvars
staticprivate

Definition at line 163 of file AbaqusUserElement.h.

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

◆ IPK_AbaqusUserElement_properties

ParamKey oofem::AbaqusUserElement::IPK_AbaqusUserElement_properties
staticprivate

Definition at line 164 of file AbaqusUserElement.h.

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

◆ IPK_AbaqusUserElement_type

ParamKey oofem::AbaqusUserElement::IPK_AbaqusUserElement_type
staticprivate

Definition at line 165 of file AbaqusUserElement.h.

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

◆ IPK_AbaqusUserElement_userElement

ParamKey oofem::AbaqusUserElement::IPK_AbaqusUserElement_userElement
staticprivate

Definition at line 160 of file AbaqusUserElement.h.

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

◆ jdltype

int* oofem::AbaqusUserElement::jdltype = nullptr
private

Definition at line 137 of file AbaqusUserElement.h.

Referenced by giveInternalForcesVector().

◆ jprops

IntArray oofem::AbaqusUserElement::jprops
private

jprops

Definition at line 143 of file AbaqusUserElement.h.

Referenced by giveInternalForcesVector().

◆ jtype

int oofem::AbaqusUserElement::jtype
private

◆ kinc

int oofem::AbaqusUserElement::kinc = 1
private

kinc, kstep

Definition at line 123 of file AbaqusUserElement.h.

Referenced by giveInternalForcesVector().

◆ kstep

int oofem::AbaqusUserElement::kstep = 1
private

Definition at line 123 of file AbaqusUserElement.h.

Referenced by giveInternalForcesVector().

◆ lFlags

IntArray oofem::AbaqusUserElement::lFlags
private

LFlags.

Definition at line 120 of file AbaqusUserElement.h.

Referenced by giveInternalForcesVector(), and postInitialize().

◆ mcrd

int oofem::AbaqusUserElement::mcrd = 2
private

mcrd

Definition at line 110 of file AbaqusUserElement.h.

Referenced by giveInternalForcesVector(), and postInitialize().

◆ mdLoad

int oofem::AbaqusUserElement::mdLoad = 0
private

Definition at line 134 of file AbaqusUserElement.h.

Referenced by giveInternalForcesVector().

◆ mlvarx

int oofem::AbaqusUserElement::mlvarx = 1
private

mlvarx

Definition at line 113 of file AbaqusUserElement.h.

Referenced by giveInternalForcesVector(), and postInitialize().

◆ nCoords

int oofem::AbaqusUserElement::nCoords
private

Coord transf.

Definition at line 83 of file AbaqusUserElement.h.

Referenced by initializeFrom(), and postInitialize().

◆ ndLoad

int oofem::AbaqusUserElement::ndLoad = 0
private

loads

Definition at line 133 of file AbaqusUserElement.h.

Referenced by giveInternalForcesVector().

◆ ndofel

int oofem::AbaqusUserElement::ndofel = 0
private

ndofel

Definition at line 107 of file AbaqusUserElement.h.

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

◆ npredef

int oofem::AbaqusUserElement::npredef = 1
private

predef

Definition at line 126 of file AbaqusUserElement.h.

Referenced by giveInternalForcesVector(), and postInitialize().

◆ nrhs

int oofem::AbaqusUserElement::nrhs = 2
private

Element rhs.

Definition at line 100 of file AbaqusUserElement.h.

Referenced by giveInternalForcesVector(), and postInitialize().

◆ numSvars

int oofem::AbaqusUserElement::numSvars
private

Number of status variables.

Definition at line 88 of file AbaqusUserElement.h.

Referenced by giveInputRecord(), giveInternalForcesVector(), initializeFrom(), and postInitialize().

◆ params

double oofem::AbaqusUserElement::params[3]
private

params

Definition at line 140 of file AbaqusUserElement.h.

Referenced by giveInternalForcesVector().

◆ predef

FloatArray oofem::AbaqusUserElement::predef
private

Definition at line 127 of file AbaqusUserElement.h.

Referenced by giveInternalForcesVector(), and postInitialize().

◆ props

FloatArray oofem::AbaqusUserElement::props
private

Element properties.

Definition at line 91 of file AbaqusUserElement.h.

Referenced by giveInputRecord(), giveInternalForcesVector(), and initializeFrom().

◆ rhs

FloatMatrix oofem::AbaqusUserElement::rhs
private

Definition at line 101 of file AbaqusUserElement.h.

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

◆ svars

FloatArray oofem::AbaqusUserElement::svars
private

Status variables.

Definition at line 94 of file AbaqusUserElement.h.

Referenced by giveStateVector(), postInitialize(), and updateYourself().

◆ tempAmatrx

FloatMatrix oofem::AbaqusUserElement::tempAmatrx
private

Definition at line 104 of file AbaqusUserElement.h.

Referenced by giveTempTangent(), letTempTangentBe(), and updateYourself().

◆ tempRHS

FloatMatrix oofem::AbaqusUserElement::tempRHS
private

Definition at line 101 of file AbaqusUserElement.h.

Referenced by letTempRhsBe(), and updateYourself().

◆ tempSvars

FloatArray oofem::AbaqusUserElement::tempSvars
private

Definition at line 94 of file AbaqusUserElement.h.

Referenced by giveTempStateVector(), letTempSvarsBe(), and updateYourself().

◆ U

FloatArray oofem::AbaqusUserElement::U
private

Inputs to element routines. Velocity and Acceleration currently ignored.

Definition at line 116 of file AbaqusUserElement.h.

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

◆ uel

void(* oofem::AbaqusUserElement::uel) (double *rhs, double *amatrx, double *svars, double energy[8], int *ndofel, int *nrhs, int *nsvars, double *props, int *nprops, double *coords, int *mcrd, int *nnode, double *u, double *du, double *v, double *a, int *jtype, double time[2], double *dtime, int *kstep, int *kinc, int *jelem, double params[3], int *ndload, int *jdltyp, double *adlmag, double *predef, int *npredef, int *lflags, int *mvarx, double *ddlmag, int *mdload, double *pnewdt, int *jprops, int *njprop, double *period)
private

Pointer to the dynamically loaded uel-function (translated to C).

Definition at line 149 of file AbaqusUserElement.h.

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

◆ uelobj

void* oofem::AbaqusUserElement::uelobj
private

Dynamically loaded uel.

Definition at line 80 of file AbaqusUserElement.h.

Referenced by AbaqusUserElement(), postInitialize(), and ~AbaqusUserElement().

◆ V

FloatArray oofem::AbaqusUserElement::V
private

Definition at line 116 of file AbaqusUserElement.h.

Referenced by giveInternalForcesVector(), and postInitialize().


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

This page is part of the OOFEM-3.0 documentation. Copyright Copyright (C) 1994-2025 Borek Patzak Bořek Patzák
Project e-mail: oofem@fsv.cvut.cz
Generated at for OOFEM by doxygen 1.15.0 written by Dimitri van Heesch, © 1997-2011