#include <nlstructuralelement.h>

Inheritance diagram for oofem::NLStructuralElement:

Collaboration diagram for oofem::NLStructuralElement:

Public Member Functions
	NLStructuralElement (int n, Domain *d)
virtual	~NLStructuralElement ()
	Destructor.
int	giveGeometryMode ()
void	computeFirstPKStressVector (FloatArray &answer, GaussPoint gp, TimeStep tStep)
void	computeCauchyStressVector (FloatArray &answer, GaussPoint gp, TimeStep tStep)
void	computeStiffnessMatrix (FloatMatrix &answer, MatResponseMode rMode, TimeStep *tStep) override
void	computeInitialStressMatrix (FloatMatrix &answer, TimeStep *tStep) override
void	computeStiffnessMatrix_withIRulesAsSubcells (FloatMatrix &answer, MatResponseMode rMode, TimeStep *tStep)
void	giveInternalForcesVector (FloatArray &answer, TimeStep *tStep, int useUpdatedGpRecord=0) override
virtual void	computeConstitutiveMatrix_dPdF_At (FloatMatrix &answer, MatResponseMode rMode, GaussPoint gp, TimeStep tStep)=0
void	giveInternalForcesVector_withIRulesAsSubcells (FloatArray &answer, TimeStep *tStep, int useUpdatedGpRecord=0) override
virtual void	computeDeformationGradientVector (FloatArray &answer, GaussPoint gp, TimeStep tStep)
double	computeCurrentVolume (TimeStep *tStep)
void	initializeFrom (InputRecord &ir, int priority) override
void	giveInputRecord (DynamicInputRecord &input) override
const char *	giveClassName () 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	computeConsistentMassMatrix (FloatMatrix &answer, TimeStep tStep, double &mass, const double ipDensity=NULL)
virtual void	giveMassMtrxIntegrationgMask (IntArray &answer)
void	computeStiffnessMatrix_withIRulesAsSubcells (FloatMatrix &answer, MatResponseMode rMode, TimeStep *tStep)
virtual void	computeLumpedInitialStressMatrix (FloatMatrix &answer, TimeStep *tStep)
void	computeField (ValueModeType mode, TimeStep *tStep, const FloatArray &lcoords, FloatArray &answer) override
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	computeStressVector (FloatArray &answer, const FloatArray &strain, GaussPoint gp, TimeStep tStep)=0
virtual void	computeBmatrixAt (GaussPoint *gp, FloatMatrix &answer, int lowerIndx=1, int upperIndx=ALL_STRAINS)=0
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)
virtual void	computeConstitutiveMatrixAt (FloatMatrix &answer, MatResponseMode rMode, GaussPoint gp, TimeStep tStep)=0
StructuralCrossSection *	giveStructuralCrossSection ()
	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
Element &	operator= (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 DofManager *	giveInternalDofManager (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 IntArray &	giveBodyLoadList () const
const IntArray &	giveBoundaryLoadList () 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	computeNumberOfDofs ()
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	giveDofManDofIDMask (int inode, IntArray &answer) const
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< IntegrationRule >	giveBoundaryEdgeIntegrationRule (int order, int boundary)
virtual std::unique_ptr< IntegrationRule >	giveBoundarySurfaceIntegrationRule (int order, int boundary)
int	giveDofManagerNumber (int i) const
const IntArray &	giveDofManArray () const
void	addDofManager (DofManager *dMan)
DofManager *	giveDofManager (int i) const
Node *	giveNode (int i) const
virtual ElementSide *	giveSide (int i) const
virtual FEInterpolation *	giveInterpolation () const
virtual FEInterpolation *	giveInterpolation (DofIDItem id) const
virtual const FEInterpolation *	getGeometryInterpolation () const
virtual Material *	giveMaterial ()
int	giveMaterialNumber () const
CrossSection *	giveCrossSection ()
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 integrationDomain	giveIntegrationDomain () const
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	giveGeometryType () const =0
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 IntegrationRule *	giveDefaultIntegrationRulePtr ()
int	giveNumberOfIntegrationRules ()
virtual IntegrationRule *	giveIntegrationRule (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 IntArray *	givePartitionList () const
void	setPartitionList (IntArray &pl)
virtual double	predictRelativeComputationalCost ()
virtual double	giveRelativeSelfComputationalCost ()
virtual double	predictRelativeRedistributionCost ()
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
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 IntArray *	giveSharedEdgeIDs () const
const IntArray *	giveSharedSurfaceIDs () const
Public Member Functions inherited from oofem::FEMComponent
	FEMComponent (int n, Domain *d)
virtual	~FEMComponent ()=default
	Virtual destructor.
virtual const char *	giveInputRecordName () const =0
Domain *	giveDomain () 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 Interface *	giveInterface (InterfaceType t)
std::string	errorInfo (const char *func) const
	Returns string for prepending output (used by error reporting macros).

Protected Member Functions
int	checkConsistency () override
virtual void	computeBHmatrixAt (GaussPoint *gp, FloatMatrix &answer)
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 ()

Protected Attributes
int	nlGeometry =0
	Flag indicating if geometrical nonlinearities apply.
Protected Attributes inherited from oofem::StructuralElement
std::unique_ptr< FloatArray >	initialDisplacements
	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.
Domain *	domain
	Link to domain object, useful for communicating with other FEM components.

Static Protected Attributes
static ParamKey	IPK_NLStructuralElement_nlgeoflag

Friends
class	GradientDamageElement
class	PhaseFieldElement
class	XfemStructuralElementInterface

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

Detailed Description

Abstract base class for "structural" finite elements with geometrical nonlinearities. The base class StructuralElement supports nonlinearities at constitutive level. The new implementation of relevant services is extended to incorporate general geometric nonlinearity. The services for computing nonlinear parts of geometrical equations for particular strain components are provided. The updated services for stiffness, strain and internal forces computations are formulated.

The activation of non-linear effects can be generally controlled on element level using nlGeometry, allowing elements to be used also in linear computations. nlGeometry = 0 - Small strain theory based on the engineering stress and strain nlGeometry = 1 - Finite deformation theory formulated in terms of the deformation gradient F and first Piola-Kirchoff stress P in the virtual work

Note: Methods ComputeStressVector and computeStrainVector as they are implemented here assume Total Lagrange approach.

Tasks :

Defining itself
Calculating its contribution to the problem:
- calculating its mass matrix M, its stiffness matrix K, its load vector f, its location array ;
- calculating its contribution to the LHS and RHS of the linear system, using Static, Newmark,etc, formula. These contributions are usually combinations of M,K,f.
Performing end-of-step operations:
- calculating the strains and stresses at its Gauss points ;
- printing its output in the data file and updating itself ;

Author: Jim Brouzoulis (among others)

Definition at line 87 of file nlstructuralelement.h.

Constructor & Destructor Documentation

◆ NLStructuralElement()

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

Constructor. Creates element with given number, belonging to given domain.

Parameters

n	Element number.
d	Domain to which new material will belong.

Definition at line 57 of file nlstructuralelement.C.

References nlGeometry, and oofem::StructuralElement::StructuralElement().

Referenced by oofem::CoupledFieldsElement::CoupledFieldsElement(), oofem::LIBeam2dNL::LIBeam2dNL(), oofem::LIBeam3d2::LIBeam3d2(), oofem::LIBeam3dNL::LIBeam3dNL(), oofem::LIBeam3dNL2::LIBeam3dNL2(), oofem::QTruss1d::QTruss1d(), oofem::Shell7Base::Shell7Base(), oofem::Structural2DElement::Structural2DElement(), oofem::Structural3DElement::Structural3DElement(), oofem::tet21ghostsolid::tet21ghostsolid(), oofem::Truss1d::Truss1d(), oofem::Truss2d::Truss2d(), and oofem::Truss3d::Truss3d().

◆ ~NLStructuralElement()

virtual oofem::NLStructuralElement::~NLStructuralElement ( )

inlinevirtual

Destructor.

Definition at line 102 of file nlstructuralelement.h.

Member Function Documentation

◆ checkConsistency()

int oofem::NLStructuralElement::checkConsistency ( void )

overrideprotectedvirtual

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

Returns: Zero value if check fail, otherwise nonzero.

Reimplemented from oofem::Element.

Reimplemented in oofem::Shell7Base, oofem::Shell7BaseXFEM, oofem::SolidShell, and oofem::TrPlaneStress2dXFEM.

Definition at line 526 of file nlstructuralelement.C.

References nlGeometry, and OOFEM_ERROR.

◆ computeBHmatrixAt()

virtual void oofem::NLStructuralElement::computeBHmatrixAt	(	GaussPoint *	gp,
		FloatMatrix &	answer )

inlineprotectedvirtual

Computes a matrix which, multiplied by the column matrix of nodal displacements, gives the displacement gradient stored by columns. The components of this matrix are derivatives of the shape functions, but they are arranged in a somewhat different way from the usual B matrix.

Parameters

gp	Integration point.
answer	BF matrix at this point.

Reimplemented in oofem::AxisymElement, oofem::LTRSpaceBoundary, oofem::PlaneStrainElement, oofem::PlaneStress2d, oofem::PlaneStress2dXfem, oofem::PlaneStressElement, oofem::QTrPlaneStress2dXFEM, oofem::QTruss1d, oofem::Quad1PlaneStrain, oofem::SolidShell, oofem::Structural2DElement, oofem::Structural3DElement, oofem::tet21ghostsolid, oofem::TrPlaneStress2dXFEM, oofem::Truss1d, oofem::Truss2d, and oofem::Truss3d.

Definition at line 247 of file nlstructuralelement.h.

References OOFEM_ERROR.

Referenced by oofem::GradDpElement::computeDeformationGradientVector(), oofem::GradientDamageElement::computeDeformationGradientVector(), computeDeformationGradientVector(), computeStiffnessMatrix(), oofem::GradientDamageElement::computeStiffnessMatrix_du(), oofem::GradDpElement::computeStiffnessMatrix_ku(), oofem::GradientDamageElement::computeStiffnessMatrix_ud(), oofem::GradDpElement::computeStiffnessMatrix_uk(), oofem::GradDpElement::computeStiffnessMatrix_uu(), oofem::GradientDamageElement::computeStiffnessMatrix_uu(), computeStiffnessMatrix_withIRulesAsSubcells(), giveInternalForcesVector(), oofem::GradientDamageElement::giveInternalForcesVector_u(), giveInternalForcesVector_withIRulesAsSubcells(), oofem::GradDpElement::giveLocalInternalForcesVector(), and oofem::XfemStructuralElementInterface::XfemElementInterface_computeDeformationGradientVector().

◆ computeCauchyStressVector()

void oofem::NLStructuralElement::computeCauchyStressVector	(	FloatArray &	answer,
		GaussPoint *	gp,
		TimeStep *	tStep )

Computes the Cauchy stress tensor on Voigt format. This method will be called if nlGeo = 1 and mode = UL. This method computes the deformation gradient F and passes it on to the crossection which then asks for the stress from the material.

Parameters

answer	Computed stress vector in Voigt form.
gp	Gauss point at which the stress is evaluated.
tStep	Time step.

Definition at line 138 of file nlstructuralelement.C.

References computeDeformationGradientVector(), oofem::StructuralCrossSection::giveCauchyStresses(), and oofem::StructuralElement::giveStructuralCrossSection().

Referenced by giveInternalForcesVector(), and giveInternalForcesVector_withIRulesAsSubcells().

◆ computeConstitutiveMatrix_dPdF_At()

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

pure virtual

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

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::AxisymElement, oofem::LIBeam2dNL, oofem::LIBeam3d2, oofem::LIBeam3dNL2, oofem::LIBeam3dNL, oofem::PlaneStrainElement, oofem::PlaneStressElement, oofem::QTruss1d, oofem::Shell7Base, oofem::Structural3DElement, oofem::tet21ghostsolid, oofem::Truss1d, oofem::Truss2d, and oofem::Truss3d.

Referenced by computeStiffnessMatrix(), and computeStiffnessMatrix_withIRulesAsSubcells().

◆ computeCurrentVolume()

double oofem::NLStructuralElement::computeCurrentVolume ( TimeStep * tStep )

Computes the current volume of element

Definition at line 66 of file nlstructuralelement.C.

References oofem::FloatMatrix::beMatrixForm(), computeDeformationGradientVector(), oofem::Element::giveDefaultIntegrationRulePtr(), oofem::FloatMatrix::giveDeterminant(), oofem::Element::giveInterpolation(), and oofem::FEInterpolation::giveTransformationJacobian().

◆ computeDeformationGradientVector()

void oofem::NLStructuralElement::computeDeformationGradientVector	(	FloatArray &	answer,
		GaussPoint *	gp,
		TimeStep *	tStep )

virtual

Computes the deformation gradient in Voigt form at integration point ip and at time step tStep. Computes the displacement gradient and adds an identitiy tensor.

Parameters

answer	Deformation gradient vector
gp	Gauss point.
tStep	Time step.

Reimplemented in oofem::LTRSpaceBoundary, oofem::PlaneStress2dXfem, oofem::QTrPlaneStress2dXFEM, oofem::tet21ghostsolid, and oofem::TrPlaneStress2dXFEM.

Definition at line 89 of file nlstructuralelement.C.

References oofem::FloatArray::at(), oofem::FloatArray::beProductOf(), computeBHmatrixAt(), oofem::Element::computeVectorOf(), oofem::GaussPoint::giveMaterialMode(), oofem::StructuralElement::initialDisplacements, OOFEM_ERROR, and oofem::FloatArray::subtract().

Referenced by computeCauchyStressVector(), computeCurrentVolume(), computeFirstPKStressVector(), and oofem::MaterialForceEvaluator::computeMaterialForce().

◆ computeFirstPKStressVector()

void oofem::NLStructuralElement::computeFirstPKStressVector	(	FloatArray &	answer,
		GaussPoint *	gp,
		TimeStep *	tStep )

Computes the first Piola-Kirchhoff stress tensor on Voigt format. This method will be called if nlGeo = 1 and mode = TL. This method computes the deformation gradient F and passes it on to the crossection which then asks for the stress from the material.

Note: P is related to S through F*S.

Parameters

answer	Computed stress vector in Voigt form.
gp	Gauss point at which the stress is evaluated.
tStep	Time step.

Definition at line 125 of file nlstructuralelement.C.

References computeDeformationGradientVector(), oofem::StructuralCrossSection::giveFirstPKStresses(), and oofem::StructuralElement::giveStructuralCrossSection().

Referenced by giveInternalForcesVector(), and giveInternalForcesVector_withIRulesAsSubcells().

◆ computeInitialStressMatrix()

void oofem::NLStructuralElement::computeInitialStressMatrix	(	FloatMatrix &	answer,
		TimeStep *	tStep )

inlineoverridevirtual

Computes the initial stiffness matrix of receiver. This method is used only if mode = UL The response is evaluated using \( \int B ({\mathrm{\sigma}}\otimes \delta )B_{\mathrm{H}} \;\mathrm{d}v \), where \( B \) is the classical B-matrix, but computed wrt updated node position

Parameters

answer	Computed initial stiffness matrix.
tStep	Time step.

Reimplemented from oofem::StructuralElement.

Reimplemented in oofem::Structural3DElement, and oofem::Truss3d.

Definition at line 158 of file nlstructuralelement.h.

References OOFEM_ERROR.

Referenced by computeStiffnessMatrix(), and oofem::GradDpElement::computeStiffnessMatrix_uu().

◆ computeStiffnessMatrix()

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

overridevirtual

Computes the stiffness matrix of receiver. The response is evaluated using \( \int B_{\mathrm{H}}^{\mathrm{T}} D B_{\mathrm{H}} \;\mathrm{d}v \), where \( B_{\mathrm{H}} \) is the B-matrix which produces the displacement gradient vector \( H_{\mathrm{V}} \) when multiplied with the solution vector a. Reduced integration are taken into account.

Parameters

answer	Computed stiffness matrix.
rMode	Response mode.
tStep	Time step.

Todo: We probably need overloaded function (like above) here as well.

Todo: Verify that it works with large deformations

Reimplemented from oofem::StructuralElement.

Reimplemented in oofem::PlaneStress2dXfem, oofem::PlaneStressGradDamage, oofem::PlaneStressPhF2d, oofem::QPlaneStrainGrad, oofem::QPlaneStrainGradDamage, oofem::QPlaneStress2dSlip, oofem::QPlaneStressGradDamage, oofem::QPlaneStressPhF2d, oofem::QTrPlaneStrainGradDamage, oofem::QTrPlaneStrainP1, oofem::QTrPlaneStress2dSlip, oofem::QTrPlaneStress2dXFEM, oofem::QTrPlaneStressGradDamage, oofem::QTRSpaceGradDamage, oofem::QTruss1dGradDamage, oofem::Quad1PlaneStrainGradDamage, oofem::Quad1PlaneStrainP1, oofem::QWedgeGradDamage, oofem::Shell7Base, oofem::Shell7BaseXFEM, oofem::SolidShell, oofem::tet21ghostsolid, oofem::TrPlaneStress2dXFEM, oofem::TrPlaneStressGradDamage, oofem::TrPlanestressRotAllman, oofem::Truss1dGradDamage, oofem::Truss3dnl2, and oofem::Truss3dnl.

Definition at line 332 of file nlstructuralelement.C.

References oofem::FloatMatrix::add(), oofem::AL, oofem::FloatMatrix::beProductOf(), oofem::FloatMatrix::beSubMatrixOf(), oofem::FloatMatrix::clear(), computeBHmatrixAt(), oofem::StructuralElement::computeBmatrixAt(), computeConstitutiveMatrix_dPdF_At(), oofem::StructuralElement::computeConstitutiveMatrixAt(), computeInitialStressMatrix(), oofem::Element::computeVolumeAround(), oofem::FEMComponent::domain, oofem::Element::giveDefaultIntegrationRulePtr(), oofem::StructuralCrossSection::giveStiffnessMatrix_dCde(), oofem::StructuralElement::giveStructuralCrossSection(), oofem::Element::integrationRulesArray, oofem::Element::isActivated(), oofem::StructuralCrossSection::isCharacteristicMtrxSymmetric(), nlGeometry, OOFEM_ERROR, oofem::FloatMatrix::plusProductSymmUpper(), oofem::FloatMatrix::plusProductUnsym(), and oofem::FloatMatrix::symmetrized().

◆ computeStiffnessMatrix_withIRulesAsSubcells()

void oofem::NLStructuralElement::computeStiffnessMatrix_withIRulesAsSubcells	(	FloatMatrix &	answer,
		MatResponseMode	rMode,
		TimeStep *	tStep )

Computes the stiffness matrix of receiver. The response is evaluated using \( \int B_{\mathrm{H}}^{\mathrm{T}} D B_{\mathrm{H}} \;\mathrm{d}v \), where \( B_{\mathrm{H}} \) is the B-matrix which produces the displacement gradient vector \( H_{\mathrm{V}} \) when multiplied with the solution vector a.

Note: Reduced intergration is not taken into account. The integration procedure uses an integrationRulesArray for numerical integration. Each integration rule is considered to represent a separate sub-cell/element. Typically this would be used when integration of the element domain needs special treatment, e.g. when using the XFEM.

Parameters

answer	Computed stiffness matrix.
rMode	Response mode.
tStep	Time step.

Definition at line 450 of file nlstructuralelement.C.

References oofem::AL, oofem::FloatMatrix::assemble(), oofem::FloatMatrix::beProductOf(), oofem::FloatMatrix::clear(), computeBHmatrixAt(), oofem::StructuralElement::computeBmatrixAt(), computeConstitutiveMatrix_dPdF_At(), oofem::StructuralElement::computeConstitutiveMatrixAt(), oofem::Element::computeVolumeAround(), oofem::FEMComponent::domain, oofem::Element::giveIntegrationRuleLocalCodeNumbers(), oofem::Element::giveInterpolation(), oofem::StructuralCrossSection::giveStiffnessMatrix_dCde(), oofem::StructuralElement::giveStructuralCrossSection(), oofem::Element::integrationRulesArray, oofem::Element::isActivated(), oofem::StructuralCrossSection::isCharacteristicMtrxSymmetric(), nlGeometry, oofem::FloatMatrix::plusProductSymmUpper(), oofem::FloatMatrix::plusProductUnsym(), and oofem::FloatMatrix::symmetrized().

◆ giveClassName()

const char * oofem::NLStructuralElement::giveClassName ( ) const

inlineoverridevirtual

Returns: Class name of the receiver.

Reimplemented from oofem::Element.

Reimplemented in oofem::PlaneStress2d, oofem::PlaneStress2dXfem, oofem::PlaneStressGradDamage, oofem::PlaneStressPhF2d, oofem::Q27Space, oofem::Q4Axisymm, oofem::Q9PlaneStress2d, oofem::QPlaneStrain, oofem::QPlaneStrainGrad, oofem::QPlaneStrainGradDamage, oofem::QPlaneStress2d, oofem::QPlaneStress2dSlip, oofem::QPlaneStressGradDamage, oofem::QPlaneStressPhF2d, oofem::QSpace, oofem::QSpaceGradDamage, oofem::QTrPlaneStrain, oofem::QTrPlaneStrainGradDamage, oofem::QTrPlaneStrainP1, oofem::QTrPlaneStress2d, oofem::QTrPlaneStress2dSlip, oofem::QTrPlaneStress2dXFEM, oofem::QTrPlaneStressGradDamage, oofem::QTRSpace, oofem::QTRSpaceGradDamage, oofem::QTruss1d, oofem::QTruss1dGradDamage, oofem::Quad1PlaneStrain, oofem::Quad1PlaneStrainGradDamage, oofem::Quad1PlaneStrainP1, oofem::QWedge, oofem::QWedgeGradDamage, oofem::Shell7Base, oofem::Shell7BaseXFEM, oofem::SolidShell, oofem::tet21ghostsolid, oofem::Tr2Shell7, oofem::Tr2Shell7XFEM, oofem::TrPlaneStrain, oofem::TrPlaneStress2d, oofem::TrPlaneStress2dXFEM, oofem::TrPlaneStressGradDamage, oofem::TrPlanestressRotAllman3d, oofem::TrPlanestressRotAllman, oofem::TrPlaneStrRot3d, oofem::TrPlaneStrRot, oofem::Truss1d, oofem::Truss1dGradDamage, oofem::Truss2d, oofem::Truss3d, oofem::Truss3dnl2, and oofem::Truss3dnl.

Definition at line 235 of file nlstructuralelement.h.

◆ giveGeometryMode()

int oofem::NLStructuralElement::giveGeometryMode ( )

inline

Returns the geometry mode describing the formulation used in the internal work 0 - Engineering (small deformation) stress-strain mode 1 - First Piola-Kirchhoff - Deformation gradient mode, P is defined as FS 2 - Second Piola-Kirchhoff - Green-Lagrange strain mode with deformation gradient as input (deprecated and not supported)

Definition at line 110 of file nlstructuralelement.h.

References nlGeometry.

Referenced by oofem::GradientDamageElement::computeStiffnessMatrix_du(), oofem::GradDpElement::computeStiffnessMatrix_ku(), oofem::GradientDamageElement::computeStiffnessMatrix_ud(), oofem::GradDpElement::computeStiffnessMatrix_uk(), oofem::GradDpElement::computeStiffnessMatrix_uu(), oofem::GradientDamageElement::computeStiffnessMatrix_uu(), oofem::BaseMixedPressureElement::computeStressVector(), oofem::GradientDamageElement::computeStressVector_and_localDamageDrivingVariable(), oofem::GradDpElement::computeStressVectorAndLocalCumulatedStrain(), oofem::GradientDamageElement::giveInternalForcesVector_u(), oofem::GradDpElement::giveLocalInternalForcesVector(), and oofem::StructuralMaterialStatus::printOutputAt().

◆ giveInputRecord()

void oofem::NLStructuralElement::giveInputRecord ( DynamicInputRecord & input )

overridevirtual

Setups the input record string of receiver.

Parameters

input Dynamic input record to be filled by receiver.

Reimplemented from oofem::Element.

Reimplemented in oofem::PlaneStress2dXfem, oofem::QTrPlaneStress2dXFEM, and oofem::TrPlaneStress2dXFEM.

Definition at line 518 of file nlstructuralelement.C.

References oofem::StructuralElement::giveInputRecord(), IPK_NLStructuralElement_nlgeoflag, nlGeometry, and oofem::DynamicInputRecord::setField().

◆ giveInternalForcesVector()

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

overridevirtual

Evaluates nodal representation of real internal forces. Necessary transformations are taken into account.

Todo: what is meant?

Parameters

answer	Equivalent nodal forces vector.
tStep	Time step
useUpdatedGpRecord	If equal to zero, the stresses in integration points are computed (slow but safe).

Todo: Is this really what we should do for inactive elements?

Todo: This is actaully inefficient since it constructs B and twice and collects the nodal unknowns over and over.

Todo: is this check really necessary?

Reimplemented from oofem::StructuralElement.

Reimplemented in oofem::PlaneStress2dXfem, oofem::PlaneStressGradDamage, oofem::PlaneStressPhF2d, oofem::QPlaneStrainGrad, oofem::QPlaneStrainGradDamage, oofem::QPlaneStress2dSlip, oofem::QPlaneStressGradDamage, oofem::QPlaneStressPhF2d, oofem::QTrPlaneStrainGradDamage, oofem::QTrPlaneStrainP1, oofem::QTrPlaneStress2dSlip, oofem::QTrPlaneStress2dXFEM, oofem::QTrPlaneStressGradDamage, oofem::QTRSpaceGradDamage, oofem::QTruss1dGradDamage, oofem::Quad1PlaneStrainGradDamage, oofem::Quad1PlaneStrainP1, oofem::QWedgeGradDamage, oofem::Shell7Base, oofem::Shell7BaseXFEM, oofem::SolidShell, oofem::tet21ghostsolid, oofem::TrPlaneStress2dXFEM, oofem::TrPlaneStressGradDamage, oofem::Truss1dGradDamage, oofem::Truss3dnl2, and oofem::Truss3dnl.

Definition at line 151 of file nlstructuralelement.C.

References oofem::AL, oofem::FloatArray::beProductOf(), oofem::FloatArray::clear(), computeBHmatrixAt(), oofem::StructuralElement::computeBmatrixAt(), computeCauchyStressVector(), computeFirstPKStressVector(), oofem::StructuralElement::computeStressVector(), oofem::Element::computeVectorOf(), oofem::Element::computeVolumeAround(), oofem::FEMComponent::domain, oofem::Element::giveCrossSection(), oofem::StructuralMaterialStatus::giveCVector(), oofem::Element::giveDefaultIntegrationRulePtr(), oofem::CrossSection::giveMaterial(), oofem::StructuralMaterialStatus::givePVector(), oofem::StructuralMaterial::giveReducedSymVectorForm(), oofem::StructuralMaterial::giveReducedVectorForm(), oofem::FloatArray::giveSize(), oofem::StructuralMaterial::giveSizeOfVoigtSymVector(), oofem::Material::giveStatus(), oofem::StructuralMaterialStatus::giveStressVector(), oofem::StructuralElement::initialDisplacements, oofem::Element::isActivated(), nlGeometry, oofem::FloatArray::plusProduct(), oofem::FloatArray::resize(), oofem::FloatArray::subtract(), and oofem::FloatArray::zero().

◆ giveInternalForcesVector_withIRulesAsSubcells()

void oofem::NLStructuralElement::giveInternalForcesVector_withIRulesAsSubcells	(	FloatArray &	answer,
		TimeStep *	tStep,
		int	useUpdatedGpRecord = 0 )

overridevirtual

Evaluates nodal representation of real internal forces.

Numerical integration procedure uses integrationRulesArray for numerical integration. The integration procedure uses an integrationRulesArray for numerical integration. Each integration rule is considered to represent a separate sub-cell/element. Typically this would be used when integration of the element domain needs special treatment, e.g. when using the XFEM.

Parameters

answer	Equivalent nodal forces vector.
tStep	Time step.
useUpdatedGpRecord	If equal to zero, the stresses in the integration points are computed (slow but safe).

Reimplemented from oofem::StructuralElement.

Definition at line 244 of file nlstructuralelement.C.

References oofem::AL, oofem::FloatArray::assemble(), oofem::FloatArray::clear(), computeBHmatrixAt(), oofem::StructuralElement::computeBmatrixAt(), computeCauchyStressVector(), computeFirstPKStressVector(), oofem::StructuralElement::computeStrainVector(), oofem::StructuralElement::computeStressVector(), oofem::Element::computeVolumeAround(), oofem::FEMComponent::domain, oofem::StructuralMaterialStatus::giveCVector(), oofem::Element::giveIntegrationRuleLocalCodeNumbers(), oofem::Element::giveInterpolation(), oofem::StructuralMaterialStatus::givePVector(), oofem::FloatArray::giveSize(), oofem::StructuralMaterialStatus::giveStressVector(), oofem::Element::integrationRulesArray, oofem::Element::isActivated(), nlGeometry, oofem::FloatArray::plusProduct(), and oofem::FloatArray::zero().