Multiscale constitutive model for subscale flow problems, typically sintering. More...

#include <fe2fluidmaterial.h>

Inheritance diagram for oofem::FE2FluidMaterial:

Collaboration diagram for oofem::FE2FluidMaterial:

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

virtual	~FE2FluidMaterial ()
	Destructor. 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 int	checkConsistency ()
	Allows programmer to test some internal data, before computation begins. More...

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

virtual void	computeDeviatoricStressVector (FloatArray &stress_dev, double &r_vol, GaussPoint gp, const FloatArray &eps, double pressure, TimeStep tStep)
	Computes the deviatoric stress vector and volumetric strain rate from given deviatoric strain and pressure. More...

virtual void	computeDeviatoricStressVector (FloatArray &answer, GaussPoint gp, const FloatArray &eps, TimeStep tStep)
	Computes the deviatoric stress vector from given strain. More...

virtual void	giveDeviatoricStiffnessMatrix (FloatMatrix &answer, MatResponseMode mode, GaussPoint gp, TimeStep tStep)
	Computes the deviatoric stiffness; $\frac{\partial\sigma_{\mathrm{dev}}}{\partial \epsilon_{\mathrm{dev}}}$ . More...

virtual void	giveStiffnessMatrices (FloatMatrix &dsdd, FloatArray &dsdp, FloatArray &dedd, double &dedp, MatResponseMode mode, GaussPoint gp, TimeStep tStep)
	Computes the 4 tangents of the compressible material response in 3D. More...

virtual double	giveEffectiveViscosity (GaussPoint gp, TimeStep tStep)
	Gives the effective viscosity for the given integration point. More...

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

virtual const char *	giveClassName () const

virtual const char *	giveInputRecordName () const

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

virtual	~FluidDynamicMaterial ()
	Destructor. More...

virtual void	updateInternalState (const FloatArray &vec, GaussPoint gp, TimeStep tStep)
	Updates internal state of material according to new state vector. More...

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

Private Attributes
std::string	inputfile

Static Private Attributes
static int	n = 1

Additional Inherited Members
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...

Detailed Description

Multiscale constitutive model for subscale flow problems, typically sintering.

The material uses the MixedGradientPressureBC to perform computational homogenization. The requirement for the supplied subscale flow problem is:

It must have a MixedGradientPressureBC. It should be the only Dirichlet boundary condition.
It must have boundary elements along the outer boundary such that the total volume (including internal pores) can be computed.

Author: Mikael Öhman

Definition at line 132 of file fe2fluidmaterial.h.

Constructor & Destructor Documentation

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

inline

Constructor.

Creates material with given number, belonging to given domain.

Parameters

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

Definition at line 144 of file fe2fluidmaterial.h.

virtual oofem::FE2FluidMaterial::~FE2FluidMaterial ( )

inlinevirtual

Destructor.

Definition at line 146 of file fe2fluidmaterial.h.

References oofem::FEMComponent::checkConsistency(), oofem::FEMComponent::giveInputRecord(), oofem::IntegrationPointStatus::gp, oofem::MaterialStatus::initializeFrom(), and oofem::FE2FluidMaterialStatus::pressure.

Member Function Documentation

int oofem::FE2FluidMaterial::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.

Reimplemented from oofem::Material.

Definition at line 307 of file fe2fluidmaterial.C.

void oofem::FE2FluidMaterial::computeDeviatoricStressVector	(	FloatArray &	stress_dev,
		double &	epsp_vol,
		GaussPoint *	gp,
		const FloatArray &	eps,
		double	pressure,
		TimeStep *	tStep
	)

virtual

Computes the deviatoric stress vector and volumetric strain rate from given deviatoric strain and pressure.

Defaults to incompressible function.

Parameters

stress_dev	Deviatoric stress.
epsp_vol	Volumetric strain-rate (minus the volumetric part of the eps argument).
gp	Integration point.
eps	Strain-rate.
pressure	Pressure.
tStep	Time step.

Reimplemented from oofem::FluidDynamicMaterial.

Definition at line 80 of file fe2fluidmaterial.C.

References oofem::FE2FluidMaterialStatus::bc, oofem::MixedGradientPressureBC::computeFields(), oofem::FE2FluidMaterialStatus::giveBC(), oofem::EngngModel::giveCurrentStep(), oofem::FE2FluidMaterialStatus::giveRVE(), oofem::FluidDynamicMaterialStatus::letDeviatoricStrainRateVectorBe(), oofem::FluidDynamicMaterialStatus::letDeviatoricStressVectorBe(), oofem::FE2FluidMaterialStatus::letPressureBe(), oofem::FE2FluidMaterialStatus::markOldTangents(), oofem::MixedGradientPressureBC::setPrescribedDeviatoricGradientFromVoigt(), oofem::MixedGradientPressureBC::setPrescribedPressure(), oofem::FE2FluidMaterialStatus::setTimeStep(), and oofem::EngngModel::solveYourselfAt().

void oofem::FE2FluidMaterial::computeDeviatoricStressVector	(	FloatArray &	answer,
		GaussPoint *	gp,
		const FloatArray &	eps,
		TimeStep *	tStep
	)

virtual

Computes the deviatoric stress vector from given strain.

Parameters

answer	Deviatoric stress.
gp	Integration point.
eps	Strain-rate.
tStep	Time step.

Implements oofem::FluidDynamicMaterial.

Definition at line 64 of file fe2fluidmaterial.C.

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

MaterialStatus * oofem::FE2FluidMaterial::CreateStatus ( GaussPoint * gp ) const

virtual

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 302 of file fe2fluidmaterial.C.

References oofem::FE2FluidMaterialStatus::FE2FluidMaterialStatus(), and oofem::FEMComponent::giveDomain().

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

inlinevirtual

Returns: Class name of the receiver.

Implements oofem::FEMComponent.

Definition at line 166 of file fe2fluidmaterial.h.

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

virtual

Computes the deviatoric stiffness; $\frac{\partial\sigma_{\mathrm{dev}}}{\partial \epsilon_{\mathrm{dev}}}$ .

Parameters

answer	Stiffness matrix.
mode	Mode of result.
gp	Integration point.
tStep	Time step.

Implements oofem::FluidDynamicMaterial.

Definition at line 105 of file fe2fluidmaterial.C.

References oofem::FE2FluidMaterialStatus::computeTangents(), oofem::FE2FluidMaterialStatus::giveDeviatoricTangent(), and OOFEM_ERROR.

double oofem::FE2FluidMaterial::giveEffectiveViscosity	(	GaussPoint *	gp,
		TimeStep *	tStep
	)

virtual

Gives the effective viscosity for the given integration point.

Parameters

gp	Gauss point of interest.
tStep	Time step.

Returns: The effective viscosity in the point.

Implements oofem::FluidDynamicMaterial.

Definition at line 421 of file fe2fluidmaterial.C.

References oofem::FloatMatrix::at(), oofem::FE2FluidMaterialStatus::computeTangents(), oofem::FE2FluidMaterialStatus::giveDeviatoricTangent(), and oofem::GaussPoint::giveMaterialMode().

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

virtual

Setups the input record string of receiver.

Parameters

input Dynamic input record to be filled by receiver.

Reimplemented from oofem::Material.

Definition at line 255 of file fe2fluidmaterial.C.

References _IFT_FE2FluidMaterial_fileName, oofem::Material::giveInputRecord(), and oofem::DynamicInputRecord::setField().

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

inlinevirtual

Returns: Input record name of the receiver.

Implements oofem::FEMComponent.

Definition at line 167 of file fe2fluidmaterial.h.

References _IFT_FE2FluidMaterial_Name.

int oofem::FE2FluidMaterial::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.

Todo:: What should one call this value? Relation between pressure and volumetric strain-rate.

Reimplemented from oofem::FluidDynamicMaterial.

Definition at line 262 of file fe2fluidmaterial.C.

References DEBUG_ERR, oofem::FluidDynamicMaterial::giveIPValue(), oofem::FE2FluidMaterialStatus::givePressure(), oofem::FE2FluidMaterialStatus::giveVOFFraction(), oofem::FE2FluidMaterialStatus::giveVolumetricPressureTangent(), oofem::norm(), OOFEM_ERROR, oofem::FE2FluidMaterialStatus::pressure, and oofem::FloatArray::zero().

void oofem::FE2FluidMaterial::giveStiffnessMatrices	(	FloatMatrix &	dsdd,
		FloatArray &	dsdp,
		FloatArray &	dedd,
		double &	dedp,
		MatResponseMode	mode,
		GaussPoint *	gp,
		TimeStep *	tStep
	)

virtual

Computes the 4 tangents of the compressible material response in 3D.

Parameters

dsdd	$\frac{\partial\sigma_{\mathrm{dev}}}{\partial \epsilon_{\mathrm{dev}}}$ .
dsdp	$\frac{\partial\sigma_{\mathrm{dev}}}{\partial p}$
dedd	$\frac{\partial\epsilon_{\mathrm{vol}}}{\partial\epsilon_{\mathrm{dev}}}$
dedp	$\frac{\partial\epsilon_{\mathrm{vol}}}{\partial p}$
mode	Mode of result.
gp	Integration point.
tStep	Time step.

Reimplemented from oofem::FluidDynamicMaterial.

Definition at line 116 of file fe2fluidmaterial.C.

References oofem::FloatArray::at(), oofem::FloatArray::beDifferenceOf(), oofem::FloatMatrix::computeFrobeniusNorm(), oofem::FloatArray::computeNorm(), oofem::FE2FluidMaterialStatus::computeTangents(), DEBUG_ERR, oofem::FE2FluidMaterialStatus::giveDeviatoricPressureTangent(), oofem::FE2FluidMaterialStatus::giveDeviatoricTangent(), oofem::FE2FluidMaterialStatus::giveVolumetricDeviatoricTangent(), oofem::FE2FluidMaterialStatus::giveVolumetricPressureTangent(), oofem::norm(), OOFEM_ERROR, oofem::FE2FluidMaterialStatus::pressure, oofem::FloatArray::printYourself(), oofem::FloatMatrix::printYourself(), oofem::FloatMatrix::setColumn(), oofem::FloatArray::subtract(), oofem::FloatMatrix::subtract(), oofem::FloatArray::times(), oofem::FloatMatrix::times(), and oofem::FloatArray::zero().

IRResultType oofem::FE2FluidMaterial::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::Material.

Definition at line 248 of file fe2fluidmaterial.C.

References _IFT_FE2FluidMaterial_fileName, oofem::Material::initializeFrom(), and IR_GIVE_FIELD.

Member Data Documentation

std :: string oofem::FE2FluidMaterial::inputfile

private

Definition at line 135 of file fe2fluidmaterial.h.

int oofem::FE2FluidMaterial::n = 1

staticprivate

Definition at line 136 of file fe2fluidmaterial.h.

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

Public Member Functions

Private Attributes

Static Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation