This class implements a fibered cross section in a finite element problem. More...

#include <fiberedcs.h>

Inheritance diagram for oofem::FiberedCrossSection:

Collaboration diagram for oofem::FiberedCrossSection:

Public Member Functions
	FiberedCrossSection (int n, Domain *d)

virtual	~FiberedCrossSection ()

virtual void	giveRealStress_3d (FloatArray &answer, GaussPoint gp, const FloatArray &reducedStrain, TimeStep tStep)

virtual void	giveRealStress_PlaneStrain (FloatArray &answer, GaussPoint gp, const FloatArray &reducedStrain, TimeStep tStep)

virtual void	giveRealStress_PlaneStress (FloatArray &answer, GaussPoint gp, const FloatArray &reducedStrain, TimeStep tStep)

virtual void	giveRealStress_1d (FloatArray &answer, GaussPoint gp, const FloatArray &reducedStrain, TimeStep tStep)

virtual void	giveRealStress_Warping (FloatArray &answer, GaussPoint gp, const FloatArray &reducedStrain, TimeStep tStep)

virtual void	giveStiffnessMatrix_3d (FloatMatrix &answer, MatResponseMode mode, GaussPoint gp, TimeStep tStep)
	Method for computing the stiffness matrix. More...

virtual void	giveStiffnessMatrix_PlaneStress (FloatMatrix &answer, MatResponseMode mode, GaussPoint gp, TimeStep tStep)

virtual void	giveStiffnessMatrix_PlaneStrain (FloatMatrix &answer, MatResponseMode mode, GaussPoint gp, TimeStep tStep)

virtual void	giveStiffnessMatrix_1d (FloatMatrix &answer, MatResponseMode mode, GaussPoint gp, TimeStep tStep)

virtual void	giveGeneralizedStress_Beam2d (FloatArray &answer, GaussPoint gp, const FloatArray &generalizedStrain, TimeStep tStep)
	Computes the generalized stress vector for given strain and integration point. More...

virtual void	giveGeneralizedStress_Beam3d (FloatArray &answer, GaussPoint gp, const FloatArray &generalizedStrain, TimeStep tStep)

virtual void	giveGeneralizedStress_Plate (FloatArray &answer, GaussPoint gp, const FloatArray &generalizedStrain, TimeStep tStep)

virtual void	giveGeneralizedStress_Shell (FloatArray &answer, GaussPoint gp, const FloatArray &generalizedStrain, TimeStep tStep)

virtual void	giveGeneralizedStress_MembraneRot (FloatArray &answer, GaussPoint gp, const FloatArray &generalizedStrain, TimeStep tStep)

virtual void	giveGeneralizedStress_PlateSubSoil (FloatArray &answer, GaussPoint gp, const FloatArray &generalizedStrain, TimeStep tStep)

virtual void	giveCharMaterialStiffnessMatrix (FloatMatrix &answer, MatResponseMode mode, GaussPoint gp, TimeStep tStep)
	Computes the stiffness matrix of receiver in given integration point, respecting its history. More...

virtual void	give2dBeamStiffMtrx (FloatMatrix &answer, MatResponseMode mode, GaussPoint gp, TimeStep tStep)
	Computes the stiffness matrix for 2d beams. More...

virtual void	give3dBeamStiffMtrx (FloatMatrix &answer, MatResponseMode mode, GaussPoint gp, TimeStep tStep)
	Computes the stiffness matrix for 2d beams. More...

virtual void	give2dPlateStiffMtrx (FloatMatrix &answer, MatResponseMode mode, GaussPoint gp, TimeStep tStep)
	Method for computing 2d plate stiffness matrix. More...

virtual void	give3dShellStiffMtrx (FloatMatrix &answer, MatResponseMode mode, GaussPoint gp, TimeStep tStep)
	Method for computing 3d shell stiffness matrix. More...

virtual void	giveMembraneRotStiffMtrx (FloatMatrix &answer, MatResponseMode mode, GaussPoint gp, TimeStep tStep)
	Method for computing membrane stiffness matrix with added drilling stiffness. More...

virtual void	give2dPlateSubSoilStiffMtrx (FloatMatrix &answer, MatResponseMode mode, GaussPoint gp, TimeStep tStep)
	Method for computing subsoil stiffness matrix for plates. More...

virtual bool	isCharacteristicMtrxSymmetric (MatResponseMode mode)
	Check for symmetry of stiffness matrix. More...

virtual double	give (int aProperty, GaussPoint *gp)
	Returns the value of cross section property. More...

virtual FloatArray *	imposeStressConstrainsOnGradient (GaussPoint gp, FloatArray gradientStressVector3d)
	Returns modified gradient of stress vector, which is used to bring stresses back to yield surface. More...

virtual FloatArray *	imposeStrainConstrainsOnGradient (GaussPoint gp, FloatArray gradientStrainVector3d)
	Returns modified gradient of strain vector, which is used to compute plastic strain increment. More...

virtual int	giveIPValue (FloatArray &answer, GaussPoint gp, InternalStateType type, TimeStep tStep)
	Returns the integration point corresponding value in Reduced form. More...

void	giveFiberMaterialStiffnessMatrix (FloatMatrix &fiberMatrix, MatResponseMode mode, GaussPoint layerGp, TimeStep tStep)
	Method for computing 1d fiber stiffness matrix of receiver. More...

virtual double	give (CrossSectionProperty a, GaussPoint *gp)
	Returns the value of cross section property at given point. More...

virtual const char *	giveInputRecordName () const

virtual const char *	giveClassName () const

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

virtual void	createMaterialStatus (GaussPoint &iGP)

virtual void	printYourself ()
	Prints receiver state on stdout. Useful for debugging. More...

double	computeIntegralThickWidth ()

MaterialMode	giveCorrespondingSlaveMaterialMode (MaterialMode)

GaussPoint *	giveSlaveGaussPoint (GaussPoint *gp, int)

virtual contextIOResultType	saveIPContext (DataStream &stream, ContextMode mode, GaussPoint *gp)
	Stores integration point state to output stream. More...

virtual contextIOResultType	restoreIPContext (DataStream &stream, ContextMode mode, GaussPoint *gp)
	Reads integration point state to output stream. More...

virtual int	checkConsistency ()
	Allows programmer to test some internal data, before computation begins. More...

virtual int	packUnknowns (DataStream &buff, TimeStep tStep, GaussPoint ip)
	Pack all necessary data of integration point (according to element parallel_mode) into given communication buffer. More...

virtual int	unpackAndUpdateUnknowns (DataStream &buff, TimeStep tStep, GaussPoint ip)
	Unpack and updates all necessary data of given integration point (according to element parallel_mode) into given communication buffer. More...

virtual int	estimatePackSize (DataStream &buff, GaussPoint *ip)
	Estimates the necessary pack size to hold all packed data of receiver. More...

virtual void	giveFirstPKStresses (FloatArray &answer, GaussPoint gp, const FloatArray &reducedFIncrement, TimeStep tStep)
	Computes the First Piola-Kirchoff stress vector for a given deformation gradient and integration point. More...

virtual void	giveCauchyStresses (FloatArray &answer, GaussPoint gp, const FloatArray &reducedFIncrement, TimeStep tStep)
	Computes the Cauchy stress vector for a given increment of deformation gradient and given integration point. More...

virtual void	giveStiffnessMatrix_dPdF (FloatMatrix &answer, MatResponseMode rMode, GaussPoint gp, TimeStep tStep)
	Computes the material stiffness matrix dPdF of receiver in a given integration point, respecting its history. More...

virtual void	giveStiffnessMatrix_dCde (FloatMatrix &answer, MatResponseMode rMode, GaussPoint gp, TimeStep tStep)
	Computes the material stiffness matrix dCde of receiver in a given integration point, respecting its history. More...

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

virtual	~StructuralCrossSection ()
	Destructor. More...

virtual void	giveEshelbyStresses (FloatArray &answer, GaussPoint gp, const FloatArray &reducedvF, TimeStep tStep)
	Computes the Eshelby stress vector. More...

virtual void	give3dBeamSubSoilStiffMtrx (FloatMatrix &answer, MatResponseMode mode, GaussPoint gp, TimeStep tStep)
	Method for computing subsoil stiffness matrix for 2d beams. More...

virtual void	give3dDegeneratedShellStiffMtrx (FloatMatrix &answer, MatResponseMode rMode, GaussPoint gp, TimeStep tStep)
	Method for computing 3d shell stiffness matrix on degenerated shell elements. More...

virtual int	testCrossSectionExtension (CrossSectExtension ext)
	Returns nonzero, if receiver implements required extension. More...

virtual Material *	giveMaterial (IntegrationPoint *ip)
	Returns the material associated with the GP. More...

virtual Interface *	giveMaterialInterface (InterfaceType t, IntegrationPoint *ip)

void	giveRealStresses (FloatArray &answer, GaussPoint gp, const FloatArray &reducedStrain, TimeStep tStep)
	Computes the real stress vector for given strain and integration point. More...

virtual void	giveRealStress_3dDegeneratedShell (FloatArray &answer, GaussPoint gp, const FloatArray &reducedStrain, TimeStep tStep)

virtual void	giveGeneralizedStress_3dBeamSubSoil (FloatArray &answer, GaussPoint gp, const FloatArray &generalizedStrain, TimeStep tStep)

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

virtual	~CrossSection ()
	Destructor. More...

int	giveSetNumber () const

virtual bool	hasProperty (CrossSectionProperty a)
	Returns true if the dictionary contains the requested property. More...

virtual double	give (CrossSectionProperty a, const FloatArray &coords, Element *elem, bool local=true)
	Returns the value of cross section property at given point (belonging to given element). More...

virtual int	setupIntegrationPoints (IntegrationRule &irule, int npoints, Element *element)
	Sets up integration rule for the given element. More...

virtual int	setupIntegrationPoints (IntegrationRule &irule, int npointsXY, int npointsZ, Element *element)
	Sets up integration rule for the given element. More...

virtual double	predictRelativeComputationalCost (GaussPoint *ip)
	Returns the weight representing relative computational cost of receiver The reference cross section is integral model in plane stress. More...

virtual double	giveRelativeSelfComputationalCost ()
	Returns the weight representing relative computational cost of receiver The reference element is integral model in plane stress. More...

virtual double	predictRelativeRedistributionCost (GaussPoint *gp)

virtual void	giveInputRecord (DynamicInputRecord &input)
	Setups the input record string of receiver. More...

Public Member Functions inherited from oofem::FEMComponent
	FEMComponent (int n, Domain *d)
	Regular constructor, creates component with given number and belonging to given domain. More...

virtual	~FEMComponent ()
	Virtual destructor. More...

Domain *	giveDomain () const

virtual void	setDomain (Domain *d)
	Sets associated Domain. More...

int	giveNumber () const

void	setNumber (int num)
	Sets number of receiver. More...

virtual void	updateLocalNumbering (EntityRenumberingFunctor &f)
	Local renumbering support. More...

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

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

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

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

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

Protected Member Functions
double	giveArea ()

Protected Attributes
IntArray	fiberMaterials
	Material of each fiber. More...

FloatArray	fiberThicks
	Thickness for each fiber. More...

FloatArray	fiberWidths
	Width for each fiber. More...

int	numberOfFibers
	Number of fibers. More...

double	thick
	Total thickness. More...

double	width
	Total width. More...

double	area
	Total area. More...

FloatArray	fiberYcoords

FloatArray	fiberZcoords

Protected Attributes inherited from oofem::CrossSection
Dictionary	propertyDictionary
	Dictionary for storing cross section parameters (like dimensions). More...

int	setNumber

Protected Attributes inherited from oofem::FEMComponent
int	number
	Component number. More...

Domain *	domain
	Link to domain object, useful for communicating with other FEM components. More...

Detailed Description

This class implements a fibered cross section in a finite element problem.

A cross section is an attribute of a domain. It is usually also attribute of many elements.

The attribute 'propertyDictionary' contains all the properties of a layered cross section, like thickness and width of each layer. The attribute 'layerMaterials' contains an array of Materials corresponding to each layer.

it uses master - slave GaussPoint approach, where master gp has more slaves gp. slave gp represent for each fiber material point. It's coordinate sections contains y,z-coordinates from mid-section. the slaves are manageg completely ( created, saved their context.,,,) from this class. Master gp only deletes slaves in destructor.

Tasks:

Returning standard material stiffness marices (like 3d stress-strain, 2d plane , plate, 3dbeam, 2d beam ..) according to current state determined by parameter StressMode by calling gp->material->GiveMaterialStiffnessMatrix (....) and by possible modifying returned matrix. (for example in layered mode approach each layer is asked for 3dMatrialStiffnes and this is integrated for example over thickness for plate bending problems)
Returning RealStress state in gauss point and for given Stress mode.
Returning a properties of cross section like thickness or area.

Definition at line 89 of file fiberedcs.h.

Constructor & Destructor Documentation

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

inline

Definition at line 102 of file fiberedcs.h.

virtual oofem::FiberedCrossSection::~FiberedCrossSection ( )

inlinevirtual

Member Function Documentation

int oofem::FiberedCrossSection::checkConsistency ( )

virtual

Allows programmer to test some internal data, before computation begins.

For example, one may use this function, to ensure that element has material with required capabilities is assigned to element. This must be done after all mesh components are instanciated.

Returns: Nonzero if receiver is consistent.

Implements oofem::StructuralCrossSection.

Definition at line 634 of file fiberedcs.C.

References oofem::IntArray::at(), fiberMaterials, oofem::FEMComponent::giveClassName(), oofem::FEMComponent::giveDomain(), oofem::Domain::giveMaterial(), oofem::IntArray::giveSize(), and OOFEM_WARNING.

Referenced by giveClassName().

double oofem::FiberedCrossSection::computeIntegralThickWidth ( )

Referenced by giveClassName().

void oofem::FiberedCrossSection::createMaterialStatus ( GaussPoint & iGP )

virtual

Implements oofem::StructuralCrossSection.

Definition at line 445 of file fiberedcs.C.

References oofem::IntArray::at(), oofem::Material::CreateStatus(), oofem::FEMComponent::domain, fiberMaterials, oofem::Domain::giveMaterial(), giveSlaveGaussPoint(), numberOfFibers, and oofem::GaussPoint::setMaterialStatus().

Referenced by giveClassName().

virtual int oofem::FiberedCrossSection::estimatePackSize	(	DataStream &	buff,
		GaussPoint *	ip
	)

inlinevirtual

Estimates the necessary pack size to hold all packed data of receiver.

The corresponding material model service is invoked. The nature of packed data is typically material model dependent.

Parameters

buff	Communication buffer.
ip	Integration point.

Returns: Estimate of pack size.

Implements oofem::CrossSection.

Definition at line 196 of file fiberedcs.h.

References OOFEM_ERROR.

virtual double oofem::FiberedCrossSection::give	(	int	aProperty,
		GaussPoint *	gp
	)

inlinevirtual

Returns the value of cross section property.

Parameters

aProperty	Id of requested property.
gp	Integration point.

Returns: Property value.

Reimplemented from oofem::CrossSection.

Definition at line 141 of file fiberedcs.h.

References giveFiberMaterialStiffnessMatrix(), giveIPValue(), imposeStrainConstrainsOnGradient(), imposeStressConstrainsOnGradient(), and OOFEM_ERROR.

double oofem::FiberedCrossSection::give	(	CrossSectionProperty	a,
		GaussPoint *	gp
	)

virtual

Returns the value of cross section property at given point.

The default implementation assumes constant properties stored in propertyDictionary.

Parameters

a	Id of requested property.
gp	Integration point

Returns: Property value.

Reimplemented from oofem::CrossSection.

Definition at line 583 of file fiberedcs.C.

References oofem::CS_Area, oofem::CS_Thickness, oofem::CS_Width, oofem::CrossSection::give(), giveArea(), thick, and width.

void oofem::FiberedCrossSection::give2dBeamStiffMtrx	(	FloatMatrix &	answer,
		MatResponseMode	mode,
		GaussPoint *	gp,
		TimeStep *	tStep
	)

virtual

Computes the stiffness matrix for 2d beams.

Parameters

answer	The requested matrix.
mode	Material response mode.
gp	Integration point.
tStep	Time step.

Implements oofem::StructuralCrossSection.

Definition at line 224 of file fiberedcs.C.

References OOFEM_ERROR.

Referenced by giveCharMaterialStiffnessMatrix(), and ~FiberedCrossSection().

void oofem::FiberedCrossSection::give2dPlateStiffMtrx	(	FloatMatrix &	answer,
		MatResponseMode	mode,
		GaussPoint *	gp,
		TimeStep *	tStep
	)

virtual

Method for computing 2d plate stiffness matrix.

Parameters

answer	Stiffness matrix.
mode	Material response mode.
gp	Integration point, which load history is used.
tStep	Time step (most models are able to respond only when tStep is current time step).

Implements oofem::StructuralCrossSection.

Definition at line 291 of file fiberedcs.C.

References OOFEM_ERROR.

Referenced by giveCharMaterialStiffnessMatrix(), and ~FiberedCrossSection().

void oofem::FiberedCrossSection::give2dPlateSubSoilStiffMtrx	(	FloatMatrix &	answer,
		MatResponseMode	mode,
		GaussPoint *	gp,
		TimeStep *	tStep
	)

virtual

Method for computing subsoil stiffness matrix for plates.

Parameters

answer	Stiffness matrix.
mode	Material response mode.
gp	Integration point, which load history is used.
tStep	Time step (most models are able to respond only when tStep is current time step).

Implements oofem::StructuralCrossSection.

Definition at line 310 of file fiberedcs.C.

References OOFEM_ERROR.

Referenced by ~FiberedCrossSection().

void oofem::FiberedCrossSection::give3dBeamStiffMtrx	(	FloatMatrix &	answer,
		MatResponseMode	mode,
		GaussPoint *	gp,
		TimeStep *	tStep
	)

virtual

Computes the stiffness matrix for 2d beams.

Parameters

answer	The requested matrix.
mode	Material response mode.
gp	Integration point.
tStep	Time step.

Todo:: This must be wrong, it will use the last evaluated G (from the last fiber), outside the loop. FIXME!

Implements oofem::StructuralCrossSection.

Definition at line 231 of file fiberedcs.C.

References oofem::FloatArray::at(), oofem::FloatMatrix::at(), fiberThicks, fiberWidths, fiberYcoords, fiberZcoords, giveFiberMaterialStiffnessMatrix(), giveSlaveGaussPoint(), numberOfFibers, oofem::FloatMatrix::resize(), and oofem::FloatMatrix::zero().

Referenced by giveCharMaterialStiffnessMatrix(), and ~FiberedCrossSection().

void oofem::FiberedCrossSection::give3dShellStiffMtrx	(	FloatMatrix &	answer,
		MatResponseMode	mode,
		GaussPoint *	gp,
		TimeStep *	tStep
	)

virtual

Method for computing 3d shell stiffness matrix.

Parameters

answer	Stiffness matrix.
mode	Material response mode.
gp	Integration point, which load history is used.
tStep	Time step (most models are able to respond only when tStep is current time step).

Implements oofem::StructuralCrossSection.

Definition at line 298 of file fiberedcs.C.

References OOFEM_ERROR.

Referenced by giveCharMaterialStiffnessMatrix(), and ~FiberedCrossSection().

double oofem::FiberedCrossSection::giveArea ( )

protected

Definition at line 597 of file fiberedcs.C.

References area, oofem::FloatArray::at(), fiberThicks, fiberWidths, and numberOfFibers.

Referenced by give(), and giveStiffnessMatrix_dCde().

virtual void oofem::FiberedCrossSection::giveCauchyStresses	(	FloatArray &	answer,
		GaussPoint *	gp,
		const FloatArray &	reducedFIncrement,
		TimeStep *	tStep
	)

inlinevirtual

Computes the Cauchy stress vector for a given increment of deformation gradient and given integration point.

The service should use previously reached equilibrium history variables. Also it should update temporary history variables in status according to newly reached state. The temporary history variables are moved into equilibrium ones after global structure equilibrium has been reached by iteration process. Elements should always pass their requests to their cross section model, which performs necessary integration over its volume and invokes necessary material services for corresponding material model defined for given integration point.

Parameters

answer	Contains the Cauchy stress.
gp	Integration point.
reducedFIncrement	Increment of the deformation gradient vector in reduced form.

Todo:: should this then be in a multiplicative way? /JB

Parameters

tStep Current time step (most models are able to respond only when tStep is current time step).

Implements oofem::StructuralCrossSection.

Definition at line 205 of file fiberedcs.h.

References OOFEM_ERROR.

void oofem::FiberedCrossSection::giveCharMaterialStiffnessMatrix	(	FloatMatrix &	answer,
		MatResponseMode	mode,
		GaussPoint *	gp,
		TimeStep *	tStep
	)

virtual

Computes the stiffness matrix of receiver in given integration point, respecting its history.

The algorithm should use temporary or equilibrium history variables stored in integration point status to compute and return required result. Elements should always pass their requests to their cross section model, which performs necessary integration over its volume and invokes necessary material services for corresponding material model defined for given integration point.

Parameters

answer	Contains result.
mode	Material response mode.
gp	Integration point.
tStep	Time step (most models are able to respond only when tStep is current time step).

Implements oofem::StructuralCrossSection.

Definition at line 203 of file fiberedcs.C.

References give2dBeamStiffMtrx(), give2dPlateStiffMtrx(), give3dBeamStiffMtrx(), give3dShellStiffMtrx(), oofem::GaussPoint::giveMaterialMode(), and OOFEM_ERROR.

Referenced by ~FiberedCrossSection().

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

inlinevirtual

Returns: Class name of the receiver.

Implements oofem::FEMComponent.

Definition at line 169 of file fiberedcs.h.

References checkConsistency(), computeIntegralThickWidth(), createMaterialStatus(), giveCorrespondingSlaveMaterialMode(), giveSlaveGaussPoint(), initializeFrom(), printYourself(), restoreIPContext(), and saveIPContext().

MaterialMode oofem::FiberedCrossSection::giveCorrespondingSlaveMaterialMode ( MaterialMode masterMode )

Definition at line 567 of file fiberedcs.C.

References OOFEM_ERROR.

Referenced by giveClassName(), and giveSlaveGaussPoint().

void oofem::FiberedCrossSection::giveFiberMaterialStiffnessMatrix	(	FloatMatrix &	fiberMatrix,
		MatResponseMode	mode,
		GaussPoint *	layerGp,
		TimeStep *	tStep
	)

Method for computing 1d fiber stiffness matrix of receiver.

Default implementation computes 3d stiffness matrix using give3dMaterialStiffnessMatrix and reduces it to 1d fiber stiffness using reduce method described above. However, this reduction is quite time consuming and if it is possible, it is recommended to overload this method and provide direct method for computing particular stiffness matrix.

Parameters

fiberMatrix	Stiffness matrix.
mode	Material response mode.
layerGp	Integration point.
tStep	Time step (most models are able to respond only when tStep is current time step).

Definition at line 611 of file fiberedcs.C.

References oofem::IntArray::at(), oofem::FEMComponent::domain, fiberMaterials, oofem::Domain::giveMaterial(), oofem::GaussPoint::giveNumber(), and oofem::StructuralMaterial::giveStiffnessMatrix().

Referenced by give(), and give3dBeamStiffMtrx().

virtual void oofem::FiberedCrossSection::giveFirstPKStresses	(	FloatArray &	answer,
		GaussPoint *	gp,
		const FloatArray &	reducedFIncrement,
		TimeStep *	tStep
	)

inlinevirtual

Computes the First Piola-Kirchoff stress vector for a given deformation gradient and integration point.

The service should use previously reached equilibrium history variables. Also it should update temporary history variables in status according to newly reached state. The temporary history variables are moved into equilibrium ones after global structure equilibrium has been reached by iteration process. Elements should always pass their requests to their cross section model, which performs necessary integration over its volume and invokes necessary material services for corresponding material model defined for given integration point.

Parameters

answer	Contains the First Piola-Kirchoff stresses.
gp	Integration point.
reducedFIncrement	Increment of the deformation gradient vector in reduced form.

Todo:: should this then be in a multiplicative way? /JB

Parameters

tStep Current time step (most models are able to respond only when tStep is current time step).

Implements oofem::StructuralCrossSection.

Definition at line 203 of file fiberedcs.h.

References OOFEM_ERROR.

void oofem::FiberedCrossSection::giveGeneralizedStress_Beam2d	(	FloatArray &	answer,
		GaussPoint *	gp,
		const FloatArray &	generalizedStrain,
		TimeStep *	tStep
	)

virtual

Computes the generalized stress vector for given strain and integration point.

Parameters

answer	Contains result.
gp	Integration point.
generalizedStrain	Strain vector in reduced generalized form.
tStep	Current time step (most models are able to respond only when tStep is current time step).

Implements oofem::StructuralCrossSection.

Definition at line 115 of file fiberedcs.C.

References OOFEM_ERROR.

Referenced by ~FiberedCrossSection().

void oofem::FiberedCrossSection::giveGeneralizedStress_Beam3d	(	FloatArray &	answer,
		GaussPoint *	gp,
		const FloatArray &	generalizedStrain,
		TimeStep *	tStep
	)

virtual

void oofem::FiberedCrossSection::giveGeneralizedStress_MembraneRot	(	FloatArray &	answer,
		GaussPoint *	gp,
		const FloatArray &	generalizedStrain,
		TimeStep *	tStep
	)

virtual

Implements oofem::StructuralCrossSection.

Definition at line 191 of file fiberedcs.C.

References OOFEM_ERROR.

Referenced by ~FiberedCrossSection().

void oofem::FiberedCrossSection::giveGeneralizedStress_Plate	(	FloatArray &	answer,
		GaussPoint *	gp,
		const FloatArray &	generalizedStrain,
		TimeStep *	tStep
	)

virtual

Implements oofem::StructuralCrossSection.

Definition at line 177 of file fiberedcs.C.

References OOFEM_ERROR.

Referenced by ~FiberedCrossSection().

void oofem::FiberedCrossSection::giveGeneralizedStress_PlateSubSoil	(	FloatArray &	answer,
		GaussPoint *	gp,
		const FloatArray &	generalizedStrain,
		TimeStep *	tStep
	)

virtual

Implements oofem::StructuralCrossSection.

Definition at line 197 of file fiberedcs.C.

References OOFEM_ERROR.

Referenced by ~FiberedCrossSection().

void oofem::FiberedCrossSection::giveGeneralizedStress_Shell	(	FloatArray &	answer,
		GaussPoint *	gp,
		const FloatArray &	generalizedStrain,
		TimeStep *	tStep
	)

virtual

Implements oofem::StructuralCrossSection.

Definition at line 184 of file fiberedcs.C.

References OOFEM_ERROR.

Referenced by ~FiberedCrossSection().

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

inlinevirtual

Returns: Input record name of the receiver.

Implements oofem::FEMComponent.

Definition at line 168 of file fiberedcs.h.

References _IFT_FiberedCrossSection_Name.

int oofem::FiberedCrossSection::giveIPValue	(	FloatArray &	answer,
		GaussPoint *	ip,
		InternalStateType	type,
		TimeStep *	tStep
	)

virtual

Returns the integration point corresponding value in Reduced form.

Parameters

answer	contain corresponding ip value, zero sized if not available
ip	Integration point.
type	Determines the type of internal variable.
tStep	Time step.

Returns: Nonzero if o.k, zero otherwise.

Todo:: For now, create material status according to the first fiber material

Reimplemented from oofem::CrossSection.

Definition at line 388 of file fiberedcs.C.

References oofem::IntArray::at(), fiberMaterials, oofem::FEMComponent::giveDomain(), oofem::CrossSection::giveIPValue(), oofem::Domain::giveMaterial(), oofem::Material::giveStatus(), oofem::StructuralMaterialStatus::giveStrainVector(), and oofem::StructuralMaterialStatus::giveStressVector().

Referenced by give().

void oofem::FiberedCrossSection::giveMembraneRotStiffMtrx	(	FloatMatrix &	answer,
		MatResponseMode	mode,
		GaussPoint *	gp,
		TimeStep *	tStep
	)

virtual

Method for computing membrane stiffness matrix with added drilling stiffness.

Parameters

answer	Stiffness matrix.
mode	Material response mode.
gp	Integration point, which load history is used.
tStep	Time step (most models are able to respond only when tStep is current time step).

Implements oofem::StructuralCrossSection.

Definition at line 304 of file fiberedcs.C.

References OOFEM_ERROR.

Referenced by ~FiberedCrossSection().

void oofem::FiberedCrossSection::giveRealStress_1d	(	FloatArray &	answer,
		GaussPoint *	gp,
		const FloatArray &	reducedStrain,
		TimeStep *	tStep
	)

virtual

Implements oofem::StructuralCrossSection.

Definition at line 71 of file fiberedcs.C.

References OOFEM_ERROR.

Referenced by ~FiberedCrossSection().

void oofem::FiberedCrossSection::giveRealStress_3d	(	FloatArray &	answer,
		GaussPoint *	gp,
		const FloatArray &	reducedStrain,
		TimeStep *	tStep
	)

virtual

Implements oofem::StructuralCrossSection.

Definition at line 50 of file fiberedcs.C.

References OOFEM_ERROR.

Referenced by ~FiberedCrossSection().

void oofem::FiberedCrossSection::giveRealStress_PlaneStrain	(	FloatArray &	answer,
		GaussPoint *	gp,
		const FloatArray &	reducedStrain,
		TimeStep *	tStep
	)

virtual

Implements oofem::StructuralCrossSection.

Definition at line 57 of file fiberedcs.C.

References OOFEM_ERROR.

Referenced by ~FiberedCrossSection().

void oofem::FiberedCrossSection::giveRealStress_PlaneStress	(	FloatArray &	answer,
		GaussPoint *	gp,
		const FloatArray &	reducedStrain,
		TimeStep *	tStep
	)

virtual

Implements oofem::StructuralCrossSection.

Definition at line 64 of file fiberedcs.C.

References OOFEM_ERROR.

Referenced by ~FiberedCrossSection().

void oofem::FiberedCrossSection::giveRealStress_Warping	(	FloatArray &	answer,
		GaussPoint *	gp,
		const FloatArray &	reducedStrain,
		TimeStep *	tStep
	)

virtual

Implements oofem::StructuralCrossSection.

Definition at line 78 of file fiberedcs.C.

References OOFEM_ERROR.

Referenced by ~FiberedCrossSection().

GaussPoint * oofem::FiberedCrossSection::giveSlaveGaussPoint	(	GaussPoint *	gp,
		int	i
	)

Definition at line 456 of file fiberedcs.C.

References oofem::FloatArray::at(), fiberYcoords, fiberZcoords, oofem::GaussPoint::gaussPoints, giveCorrespondingSlaveMaterialMode(), oofem::GaussPoint::giveIntegrationRule(), oofem::GaussPoint::giveMaterialMode(), oofem::GaussPoint::giveSlaveGaussPoint(), numberOfFibers, and OOFEM_ERROR.

Referenced by createMaterialStatus(), give3dBeamStiffMtrx(), giveClassName(), giveGeneralizedStress_Beam3d(), restoreIPContext(), and saveIPContext().

void oofem::FiberedCrossSection::giveStiffnessMatrix_1d	(	FloatMatrix &	answer,
		MatResponseMode	mode,
		GaussPoint *	gp,
		TimeStep *	tStep
	)

virtual

Implements oofem::StructuralCrossSection.

Definition at line 106 of file fiberedcs.C.

References OOFEM_ERROR.

Referenced by ~FiberedCrossSection().

void oofem::FiberedCrossSection::giveStiffnessMatrix_3d	(	FloatMatrix &	answer,
		MatResponseMode	mode,
		GaussPoint *	gp,
		TimeStep *	tStep
	)

virtual

Method for computing the stiffness matrix.

Parameters

answer	Stiffness matrix.
mode	Material response mode.
gp	Integration point, which load history is used.
tStep	Time step (most models are able to respond only when tStep is current time step).

Implements oofem::StructuralCrossSection.

Definition at line 85 of file fiberedcs.C.

References OOFEM_ERROR.

Referenced by ~FiberedCrossSection().

virtual void oofem::FiberedCrossSection::giveStiffnessMatrix_dCde	(	FloatMatrix &	answer,
		MatResponseMode	mode,
		GaussPoint *	gp,
		TimeStep *	tStep
	)

inlinevirtual

Computes the material stiffness matrix dCde of receiver in a given integration point, respecting its history.

The algorithm should use temporary or equilibrium history variables stored in integration point status to compute and return required result. Elements should always pass their requests to their cross section model, which performs necessary integration over its volume and invokes necessary material services for corresponding material model defined for given integration point.

Parameters

answer	Contains result.
mode	Material response mode.
gp	Integration point.
tStep	Time step (most models are able to respond only when tStep is current time step).

Implements oofem::StructuralCrossSection.

Definition at line 209 of file fiberedcs.h.

References giveArea(), and OOFEM_ERROR.

virtual void oofem::FiberedCrossSection::giveStiffnessMatrix_dPdF	(	FloatMatrix &	answer,
		MatResponseMode	mode,
		GaussPoint *	gp,
		TimeStep *	tStep
	)

inlinevirtual

Computes the material stiffness matrix dPdF of receiver in a given integration point, respecting its history.

The algorithm should use temporary or equilibrium history variables stored in integration point status to compute and return required result. Elements should always pass their requests to their cross section model, which performs necessary integration over its volume and invokes necessary material services for corresponding material model defined for given integration point.

Parameters

answer	Contains result.
mode	Material response mode.
gp	Integration point.
tStep	Time step (most models are able to respond only when tStep is current time step).

Implements oofem::StructuralCrossSection.

Definition at line 207 of file fiberedcs.h.

References OOFEM_ERROR.

void oofem::FiberedCrossSection::giveStiffnessMatrix_PlaneStrain	(	FloatMatrix &	answer,
		MatResponseMode	mode,
		GaussPoint *	gp,
		TimeStep *	tStep
	)

virtual

Implements oofem::StructuralCrossSection.

Definition at line 99 of file fiberedcs.C.

References OOFEM_ERROR.

Referenced by ~FiberedCrossSection().

void oofem::FiberedCrossSection::giveStiffnessMatrix_PlaneStress	(	FloatMatrix &	answer,
		MatResponseMode	mode,
		GaussPoint *	gp,
		TimeStep *	tStep
	)

virtual

Implements oofem::StructuralCrossSection.

Definition at line 92 of file fiberedcs.C.

References OOFEM_ERROR.

Referenced by ~FiberedCrossSection().

FloatArray * oofem::FiberedCrossSection::imposeStrainConstrainsOnGradient	(	GaussPoint *	gp,
		FloatArray *	gradientStressVector3d
	)

virtual

Returns modified gradient of strain vector, which is used to compute plastic strain increment.

Imposes zeros on places, where zero strain occurs or energetically connected stress is prescribed to be zero.

See also: imposeStressConstrainsOnGradient

Parameters

gp	Integration point.
gradientStressVector3d	General 3d stress gradient.

Reimplemented from oofem::StructuralCrossSection.

Definition at line 356 of file fiberedcs.C.

References oofem::FloatArray::at(), oofem::GaussPoint::giveMaterialMode(), oofem::FloatArray::giveSize(), oofem::StructuralCrossSection::imposeStrainConstrainsOnGradient(), and OOFEM_ERROR.

Referenced by give().

FloatArray * oofem::FiberedCrossSection::imposeStressConstrainsOnGradient	(	GaussPoint *	gp,
		FloatArray *	gradientStressVector3d
	)

virtual

Returns modified gradient of stress vector, which is used to bring stresses back to yield surface.

Method imposes zeros on places, where zero stress occurs. if energetically connected strain is zero, we do not impose zero there, because stress exist and must be taken into account when computing yield function. In such case a problem is assumed to be full 3d with some explicit strain equal to 0. On the other hand, if some stress is imposed to be zero, we understand such case as subspace of 3d case (like a classical plane stress problem, with no tracing of e_z, sigma_z)

Parameters

gp	Integration point.
gradientStressVector3d	General 3d stress gradient.

Reimplemented from oofem::StructuralCrossSection.

Definition at line 318 of file fiberedcs.C.

References oofem::FloatArray::at(), oofem::GaussPoint::giveMaterialMode(), oofem::FloatArray::giveSize(), oofem::StructuralCrossSection::imposeStressConstrainsOnGradient(), and OOFEM_ERROR.

Referenced by give().

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

virtual

Initializes receiver according to object description stored in input record.

This function is called immediately after creating object using constructor. Input record can be imagined as data record in component database belonging to receiver. Receiver may use value-name extracting functions to extract particular field from record.

See also: IR_GIVE_FIELD; IR_GIVE_OPTIONAL_FIELD

Parameters

ir	Input record to initialize from.

Returns: IRResultType

Reimplemented from oofem::CrossSection.

Definition at line 405 of file fiberedcs.C.

Referenced by giveClassName().

bool oofem::FiberedCrossSection::isCharacteristicMtrxSymmetric ( MatResponseMode rMode )

virtual

Check for symmetry of stiffness matrix.

Default implementation returns true. It can be moved to base Cross section class in the future.

Parameters

rMode Response mode of material.

Returns: True if stiffness matrix of receiver is symmetric.

Todo:: As far as I can see, it only uses diagonal components for the 3dbeam, but there is no way to check here.

Implements oofem::StructuralCrossSection.

Definition at line 620 of file fiberedcs.C.

References oofem::FEMComponent::domain, oofem::Domain::giveMaterial(), and numberOfFibers.

Referenced by ~FiberedCrossSection().

virtual int oofem::FiberedCrossSection::packUnknowns	(	DataStream &	buff,
		TimeStep *	tStep,
		GaussPoint *	ip
	)

inlinevirtual

Pack all necessary data of integration point (according to element parallel_mode) into given communication buffer.

The corresponding material model service for particular integration point is invoked. The nature of packed data is material model dependent. Typically, for material of "local" response (response depends only on integration point local state) no data are exchanged. For "nonlocal" constitutive models the send/receive of local values which undergo averaging is performed between local and corresponding remote elements.

Parameters