Simple isotropic damage based model for interface elements. More...

#include <intmatisodamage.h>

Inheritance diagram for oofem::IntMatIsoDamage:

Collaboration diagram for oofem::IntMatIsoDamage:

Public Member Functions
	IntMatIsoDamage (int n, Domain *d)
	Constructor. More...

virtual	~IntMatIsoDamage ()
	Destructor. More...

virtual const char *	giveInputRecordName () const

virtual bool	hasAnalyticalTangentStiffness () const
	Tells if the model has implemented analytical tangent stiffness. More...

virtual void	giveEngTraction_3d (FloatArray &answer, GaussPoint gp, const FloatArray &jump, TimeStep tStep)

virtual void	giveFirstPKTraction_3d (FloatArray &answer, GaussPoint gp, const FloatArray &jump, const FloatMatrix &F, TimeStep tStep)

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

virtual void	computeEquivalentJump (double &kappa, const FloatArray &jump)
	Computes the equivalent jump measure from given jump vector (full form). More...

virtual void	computeDamageParam (double &omega, double kappa)
	computes the value of damage parameter omega, based on given value of equivalent strain. More...

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

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

virtual MaterialStatus *	CreateStatus (GaussPoint *gp) const
	Creates new copy of associated status and inserts it into given integration point. More...

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

virtual	~StructuralInterfaceMaterial ()
	Destructor. More...

virtual void	giveFirstPKTraction_1d (FloatArray &answer, GaussPoint gp, const FloatArray &jump, const FloatMatrix &reducedF, TimeStep tStep)
	Computes the first Piola-Kirchoff traction vector for given total jump/gap and integration point. More...

virtual void	giveFirstPKTraction_2d (FloatArray &answer, GaussPoint gp, const FloatArray &jump, const FloatMatrix &reducedF, TimeStep tStep)

virtual void	giveEngTraction_1d (FloatArray &answer, GaussPoint gp, const FloatArray &jump, TimeStep tStep)

virtual void	giveEngTraction_2d (FloatArray &answer, GaussPoint gp, const FloatArray &jump, TimeStep tStep)

virtual void	give1dStiffnessMatrix_dTdj (FloatMatrix &answer, MatResponseMode rMode, GaussPoint gp, TimeStep tStep)
	Gives the tangent: $\frac{\partial T}{\partial j}$ . More...

virtual void	give2dStiffnessMatrix_dTdj (FloatMatrix &answer, MatResponseMode rMode, GaussPoint gp, TimeStep tStep)

virtual void	give3dStiffnessMatrix_dTdj (FloatMatrix &answer, MatResponseMode rMode, GaussPoint gp, TimeStep tStep)

virtual void	give1dStiffnessMatrix_Eng (FloatMatrix &answer, MatResponseMode rMode, GaussPoint gp, TimeStep tStep)

void	give1dStiffnessMatrix_dTdj_Num (FloatMatrix &answer, GaussPoint gp, TimeStep tStep)

void	give2dStiffnessMatrix_dTdj_Num (FloatMatrix &answer, GaussPoint gp, TimeStep tStep)

void	give3dStiffnessMatrix_dTdj_Num (FloatMatrix &answer, GaussPoint gp, TimeStep tStep)

void	give1dStiffnessMatrix_Eng_Num (FloatMatrix &answer, GaussPoint gp, TimeStep tStep)

void	give2dStiffnessMatrix_Eng_Num (FloatMatrix &answer, GaussPoint gp, TimeStep tStep)

void	give3dStiffnessMatrix_Eng_Num (FloatMatrix &answer, GaussPoint gp, TimeStep tStep)

virtual const char *	giveClassName () const

virtual FloatArray	giveInterfaceStrength ()

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

virtual	~Material ()
	Destructor. More...

virtual bool	isCharacteristicMtrxSymmetric (MatResponseMode rMode)
	Returns true if stiffness matrix of receiver is symmetric Default implementation returns true. More...

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

virtual bool	hasProperty (int aProperty, GaussPoint *gp)
	Returns true if 'aProperty' exists on material. More...

virtual void	modifyProperty (int aProperty, double value, GaussPoint *gp)
	Modify 'aProperty', which already exists on material. More...

double	giveCastingTime ()

virtual bool	isActivated (TimeStep *tStep)

virtual int	hasNonLinearBehaviour ()
	Returns nonzero if receiver is non linear. More...

virtual int	hasMaterialModeCapability (MaterialMode mode)
	Tests if material supports material mode. More...

virtual int	hasCastingTimeSupport ()
	Tests if material supports casting time. More...

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

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	initMaterial (Element *element)
	Optional function to call specific procedures when initializing a material. More...

virtual MaterialStatus *	giveStatus (GaussPoint *gp) const
	Returns material status of receiver in given integration point. 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 double	predictRelativeComputationalCost (GaussPoint *gp)
	Returns the weight representing relative computational cost of receiver The reference material model is linear isotropic material - its weight is set to 1.0 The other material models should compare to this reference model. More...

virtual double	predictRelativeRedistributionCost (GaussPoint *gp)
	Returns the relative redistribution cost of the receiver. More...

virtual void	initTempStatus (GaussPoint *gp)
	Initializes temporary variables stored in integration point status at the beginning of new time step. More...

virtual int	setIPValue (const FloatArray &value, GaussPoint *gp, InternalStateType type)
	Sets the value of a certain variable at a given integration point to the given value. 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
virtual void	give2dStiffnessMatrix_Eng (FloatMatrix &answer, MatResponseMode rMode, GaussPoint gp, TimeStep tStep)

virtual void	give3dStiffnessMatrix_Eng (FloatMatrix &answer, MatResponseMode rMode, GaussPoint gp, TimeStep tStep)

Protected Attributes
double	kn
	Elastic properties (normal moduli). More...

double	ks
	Shear moduli. More...

double	ft
	Tension strength. More...

double	gf
	Fracture energy. More...

double	e0
	Limit elastic deformation. More...

double	maxOmega
	Maximum limit on omega. The purpose is elimination of a too compliant material which may cause convergency problems. Set to something like 0.99 if needed. More...

bool	semiExplicit

Protected Attributes inherited from oofem::Material
Dictionary	propertyDictionary
	Property dictionary. More...

double	castingTime
	Casting time. More...

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...

Additional Inherited Members
Public Attributes inherited from oofem::StructuralInterfaceMaterial
bool	useNumericalTangent

Detailed Description

Simple isotropic damage based model for interface elements.

In 2d, the interface elements are used to model contact layer between element edges. The jump vector contains the relative displacements (in normal and shear direction). The traction vector contains the corresponding tractions in normal and tangent direction.

The behaviour of the model is elastic, described by normal and shear stiffness components. Isotropic damage is initiated when the stress reaches the tensile strength. The damage evolution is governed by the normal jump

Definition at line 113 of file intmatisodamage.h.

Constructor & Destructor Documentation

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

Constructor.

Definition at line 48 of file intmatisodamage.C.

References maxOmega, and semiExplicit.

oofem::IntMatIsoDamage::~IntMatIsoDamage ( )

virtual

Destructor.

Definition at line 55 of file intmatisodamage.C.

Member Function Documentation

void oofem::IntMatIsoDamage::computeDamageParam	(	double &	omega,
		double	kappa
	)

virtual

computes the value of damage parameter omega, based on given value of equivalent strain.

Parameters

[out]	omega	Contains result.
	kappa	Equivalent strain measure.
	strain	Total strain vector in full form. (unnecessary?)
	gp	Integration point.

Reimplemented in oofem::IntMatIsoDamageTable.

Definition at line 296 of file intmatisodamage.C.

References e0, ft, and gf.

Referenced by give2dStiffnessMatrix_Eng(), and giveEngTraction_3d().

void oofem::IntMatIsoDamage::computeEquivalentJump	(	double &	kappa,
		const FloatArray &	jump
	)

virtual

Computes the equivalent jump measure from given jump vector (full form).

Parameters

[out]	kappa	Return parameter containing the corresponding equivalent jump.
	jump	Jump vector in full form.
	gp	Integration point.
	tStep	Time step.

Definition at line 290 of file intmatisodamage.C.

References oofem::FloatArray::at().

Referenced by giveEngTraction_3d().

virtual MaterialStatus* oofem::IntMatIsoDamage::CreateStatus ( GaussPoint * gp ) const

inlinevirtual

Creates new copy of associated status and inserts it into given integration point.

Parameters

gp	Integration point where newly created status will be stored.

Returns: Reference to new status.

Reimplemented from oofem::Material.

Definition at line 170 of file intmatisodamage.h.

References oofem::FEMComponent::domain, and oofem::IntMatIsoDamageStatus::IntMatIsoDamageStatus().

void oofem::IntMatIsoDamage::give2dStiffnessMatrix_Eng	(	FloatMatrix &	answer,
		MatResponseMode	rMode,
		GaussPoint *	gp,
		TimeStep *	tStep
	)

protectedvirtual

Todo:: this is only the secant stiffness - tangent stiffness is broken!

Reimplemented from oofem::StructuralInterfaceMaterial.

Definition at line 119 of file intmatisodamage.C.

References oofem::FloatArray::at(), oofem::FloatMatrix::at(), computeDamageParam(), e0, ft, gf, oofem::IntMatIsoDamageStatus::giveKappa(), oofem::Material::giveStatus(), oofem::IntMatIsoDamageStatus::giveTempDamage(), oofem::StructuralInterfaceMaterialStatus::giveTempJump(), oofem::IntMatIsoDamageStatus::giveTempKappa(), kn, ks, maxOmega, oofem::min(), oofem::FloatMatrix::resize(), and oofem::FloatMatrix::zero().

void oofem::IntMatIsoDamage::give3dStiffnessMatrix_Eng	(	FloatMatrix &	answer,
		MatResponseMode	rMode,
		GaussPoint *	gp,
		TimeStep *	tStep
	)

protectedvirtual

Reimplemented from oofem::StructuralInterfaceMaterial.

Definition at line 177 of file intmatisodamage.C.

References oofem::FloatArray::at(), oofem::FloatMatrix::at(), oofem::Material::giveStatus(), oofem::IntMatIsoDamageStatus::giveTempDamage(), oofem::StructuralInterfaceMaterialStatus::giveTempJump(), kn, ks, maxOmega, oofem::min(), oofem::FloatMatrix::resize(), and oofem::FloatMatrix::zero().

void oofem::IntMatIsoDamage::giveEngTraction_3d	(	FloatArray &	answer,
		GaussPoint *	gp,
		const FloatArray &	jump,
		TimeStep *	tStep
	)

virtual

Reimplemented from oofem::StructuralInterfaceMaterial.

Definition at line 60 of file intmatisodamage.C.

References oofem::FloatArray::at(), computeDamageParam(), computeEquivalentJump(), oofem::IntMatIsoDamageStatus::giveDamage(), oofem::IntMatIsoDamageStatus::giveKappa(), oofem::Material::giveStatus(), kn, ks, oofem::StructuralInterfaceMaterialStatus::letTempJumpBe(), oofem::StructuralInterfaceMaterialStatus::letTempTractionBe(), maxOmega, oofem::min(), semiExplicit, oofem::IntMatIsoDamageStatus::setTempDamage(), and oofem::IntMatIsoDamageStatus::setTempKappa().

Referenced by giveFirstPKTraction_3d().

void oofem::IntMatIsoDamage::giveFirstPKTraction_3d	(	FloatArray &	answer,
		GaussPoint *	gp,
		const FloatArray &	jump,
		const FloatMatrix &	F,
		TimeStep *	tStep
	)

virtual

Reimplemented from oofem::StructuralInterfaceMaterial.

Definition at line 109 of file intmatisodamage.C.

References giveEngTraction_3d(), oofem::Material::giveStatus(), and oofem::StructuralInterfaceMaterialStatus::letTempFirstPKTractionBe().

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

virtual

Setups the input record string of receiver.

Parameters

input Dynamic input record to be filled by receiver.

Reimplemented from oofem::StructuralInterfaceMaterial.

Reimplemented in oofem::IntMatIsoDamageTable.

Definition at line 275 of file intmatisodamage.C.

References _IFT_IntMatIsoDamage_ft, _IFT_IntMatIsoDamage_gf, _IFT_IntMatIsoDamage_kn, _IFT_IntMatIsoDamage_ks, _IFT_IntMatIsoDamage_maxOmega, ft, gf, oofem::StructuralInterfaceMaterial::giveInputRecord(), kn, ks, maxOmega, and oofem::DynamicInputRecord::setField().

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

inlinevirtual

Returns: Input record name of the receiver.

Implements oofem::FEMComponent.

Reimplemented in oofem::IntMatIsoDamageTable.

Definition at line 137 of file intmatisodamage.h.

References _IFT_IntMatIsoDamage_Name.

int oofem::IntMatIsoDamage::giveIPValue	(	FloatArray &	answer,
		GaussPoint *	gp,
		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.
gp	Integration point to which the value refers.
type	Determines the type of internal variable.
tStep	Determines the time step.

Returns: Nonzero if the assignment can be done, zero if this type of variable is not supported.

Reimplemented from oofem::StructuralInterfaceMaterial.

Definition at line 221 of file intmatisodamage.C.

References oofem::FloatArray::at(), oofem::IntMatIsoDamageStatus::giveDamage(), oofem::StructuralInterfaceMaterial::giveIPValue(), oofem::IntMatIsoDamageStatus::giveKappa(), oofem::Material::giveStatus(), oofem::IntMatIsoDamageStatus::giveTempDamage(), oofem::FloatArray::resize(), and oofem::FloatArray::zero().

virtual bool oofem::IntMatIsoDamage::hasAnalyticalTangentStiffness ( ) const

inlinevirtual

Tells if the model has implemented analytical tangent stiffness.

If not, the tangent must be computed numerically.

Implements oofem::StructuralInterfaceMaterial.

Definition at line 138 of file intmatisodamage.h.

References oofem::StructuralInterfaceMaterialStatus::F, oofem::FEMComponent::giveInputRecord(), oofem::IntegrationPointStatus::gp, oofem::MaterialStatus::initializeFrom(), oofem::StructuralInterfaceMaterialStatus::jump, and oofem::IntMatIsoDamageStatus::kappa.

IRResultType oofem::IntMatIsoDamage::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::StructuralInterfaceMaterial.

Reimplemented in oofem::IntMatIsoDamageTable.

Definition at line 253 of file intmatisodamage.C.

References _IFT_IntMatIsoDamage_ft, _IFT_IntMatIsoDamage_gf, _IFT_IntMatIsoDamage_kn, _IFT_IntMatIsoDamage_ks, _IFT_IntMatIsoDamage_maxOmega, e0, ft, gf, oofem::StructuralInterfaceMaterial::initializeFrom(), IR_GIVE_FIELD, IR_GIVE_OPTIONAL_FIELD, kn, ks, oofem::max(), maxOmega, and oofem::min().

Member Data Documentation

double oofem::IntMatIsoDamage::e0

protected

Limit elastic deformation.

Definition at line 125 of file intmatisodamage.h.

Referenced by oofem::IntMatIsoDamageTable::computeDamageParam(), computeDamageParam(), give2dStiffnessMatrix_Eng(), oofem::IntMatIsoDamageTable::initializeFrom(), and initializeFrom().

double oofem::IntMatIsoDamage::ft

protected

Tension strength.

Definition at line 121 of file intmatisodamage.h.

Referenced by computeDamageParam(), give2dStiffnessMatrix_Eng(), oofem::IntMatIsoDamageTable::giveInputRecord(), giveInputRecord(), oofem::IntMatIsoDamageTable::initializeFrom(), and initializeFrom().

double oofem::IntMatIsoDamage::gf

protected

Fracture energy.

Definition at line 123 of file intmatisodamage.h.

Referenced by computeDamageParam(), give2dStiffnessMatrix_Eng(), giveInputRecord(), and initializeFrom().

double oofem::IntMatIsoDamage::kn

protected

Elastic properties (normal moduli).

Definition at line 117 of file intmatisodamage.h.

Referenced by give2dStiffnessMatrix_Eng(), give3dStiffnessMatrix_Eng(), giveEngTraction_3d(), oofem::IntMatIsoDamageTable::giveInputRecord(), giveInputRecord(), oofem::IntMatIsoDamageTable::initializeFrom(), and initializeFrom().

double oofem::IntMatIsoDamage::ks

protected

Shear moduli.

Definition at line 119 of file intmatisodamage.h.

Referenced by give2dStiffnessMatrix_Eng(), give3dStiffnessMatrix_Eng(), giveEngTraction_3d(), oofem::IntMatIsoDamageTable::giveInputRecord(), giveInputRecord(), oofem::IntMatIsoDamageTable::initializeFrom(), and initializeFrom().

double oofem::IntMatIsoDamage::maxOmega

protected

Maximum limit on omega. The purpose is elimination of a too compliant material which may cause convergency problems. Set to something like 0.99 if needed.

Definition at line 127 of file intmatisodamage.h.

Referenced by give2dStiffnessMatrix_Eng(), give3dStiffnessMatrix_Eng(), giveEngTraction_3d(), oofem::IntMatIsoDamageTable::giveInputRecord(), giveInputRecord(), oofem::IntMatIsoDamageTable::initializeFrom(), initializeFrom(), and IntMatIsoDamage().

bool oofem::IntMatIsoDamage::semiExplicit

protected

Definition at line 129 of file intmatisodamage.h.

Referenced by giveEngTraction_3d(), and IntMatIsoDamage().

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

Public Member Functions

Protected Member Functions

Protected Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation