#include <fe2fluidmaterial.h>

Inheritance diagram for oofem::FE2FluidMaterial:

Collaboration diagram for oofem::FE2FluidMaterial:

Public Member Functions
	FE2FluidMaterial (int n, Domain *d)
void	initializeFrom (InputRecord &ir) override
void	giveInputRecord (DynamicInputRecord &input) override
int	checkConsistency () override
std::unique_ptr< MaterialStatus >	CreateStatus (GaussPoint *gp) const override
std::pair< FloatArrayF< 6 >, double >	computeDeviatoricStress3D (const FloatArrayF< 6 > &eps, double pressure, GaussPoint gp, TimeStep tStep) const override
FloatArrayF< 6 >	computeDeviatoricStress3D (const FloatArrayF< 6 > &eps, GaussPoint gp, TimeStep tStep) const override
FloatMatrixF< 6, 6 >	computeTangent3D (MatResponseMode mode, GaussPoint gp, TimeStep tStep) const override
Tangents< 6 >	computeTangents3D (MatResponseMode mode, GaussPoint gp, TimeStep tStep) const override
double	giveEffectiveViscosity (GaussPoint gp, TimeStep tStep) const override
int	giveIPValue (FloatArray &answer, GaussPoint gp, InternalStateType type, TimeStep tStep) override
const char *	giveClassName () const override
const char *	giveInputRecordName () const override
Public Member Functions inherited from oofem::FluidDynamicMaterial
	FluidDynamicMaterial (int n, Domain *d)
virtual std::pair< FloatArrayF< 3 >, double >	computeDeviatoricStress2D (const FloatArrayF< 3 > &eps, double pressure, GaussPoint gp, TimeStep tStep) const
virtual FloatArrayF< 3 >	computeDeviatoricStress2D (const FloatArrayF< 3 > &eps, GaussPoint gp, TimeStep tStep) const
virtual FloatArrayF< 4 >	computeDeviatoricStressAxi (const FloatArrayF< 4 > &eps, GaussPoint gp, TimeStep tStep) const
virtual FloatMatrixF< 3, 3 >	computeTangent2D (MatResponseMode mode, GaussPoint gp, TimeStep tStep) const
virtual FloatMatrixF< 4, 4 >	computeTangentAxi (MatResponseMode mode, GaussPoint gp, TimeStep tStep) const
virtual Tangents< 3 >	computeTangents2D (MatResponseMode mode, GaussPoint gp, TimeStep tStep) const
int	giveIPValue (FloatArray &answer, GaussPoint gp, InternalStateType type, TimeStep tStep) override
Public Member Functions inherited from oofem::Material
	Material (int n, Domain *d)
virtual	~Material ()=default
	Destructor.
virtual bool	isCharacteristicMtrxSymmetric (MatResponseMode rMode) const
virtual void	giveCharacteristicMatrix (FloatMatrix &answer, MatResponseMode type, GaussPoint gp, TimeStep tStep) const
	Returns characteristic matrix of the receiver.
virtual void	giveCharacteristicVector (FloatArray &answer, FloatArray &flux, MatResponseMode type, GaussPoint gp, TimeStep tStep) const
	Returns characteristic vector of the receiver.
virtual double	giveCharacteristicValue (MatResponseMode type, GaussPoint gp, TimeStep tStep) const
	Returns characteristic value of the receiver.
virtual double	give (int aProperty, GaussPoint *gp) const
virtual bool	hasProperty (int aProperty, GaussPoint *gp) const
virtual void	modifyProperty (int aProperty, double value, GaussPoint *gp)
double	giveCastingTime () const
virtual bool	isActivated (TimeStep *tStep) const
virtual bool	hasMaterialModeCapability (MaterialMode mode) const
virtual bool	hasCastingTimeSupport () const
virtual int	setIPValue (const FloatArray &value, GaussPoint *gp, InternalStateType type)
void	initializeFrom (InputRecord &ir) override
void	giveInputRecord (DynamicInputRecord &input) override
void	printYourself () override
	Prints receiver state on stdout. Useful for debugging.
virtual void	saveIPContext (DataStream &stream, ContextMode mode, GaussPoint *gp)
virtual void	restoreIPContext (DataStream &stream, ContextMode mode, GaussPoint *gp)
int	checkConsistency () override
virtual void	restoreConsistency (GaussPoint *gp)
virtual int	initMaterial (Element *element)
virtual MaterialStatus *	giveStatus (GaussPoint *gp) const
virtual int	packUnknowns (DataStream &buff, TimeStep tStep, GaussPoint ip)
virtual int	unpackAndUpdateUnknowns (DataStream &buff, TimeStep tStep, GaussPoint ip)
virtual int	estimatePackSize (DataStream &buff, GaussPoint *ip)
virtual double	predictRelativeComputationalCost (GaussPoint *gp)
virtual double	predictRelativeRedistributionCost (GaussPoint *gp)
virtual void	initTempStatus (GaussPoint *gp) const
void	saveContext (DataStream &stream, ContextMode mode) override
void	restoreContext (DataStream &stream, ContextMode mode) override
Public Member Functions inherited from oofem::FEMComponent
	FEMComponent (int n, Domain *d)
virtual	~FEMComponent ()=default
	Virtual destructor.
Domain *	giveDomain () const
virtual void	setDomain (Domain *d)
int	giveNumber () const
void	setNumber (int num)
virtual void	updateLocalNumbering (EntityRenumberingFunctor &f)
virtual void	initializeFrom (InputRecord &ir, int priority)
virtual void	initializeFinish ()
virtual void	postInitialize ()
	Performs post initialization steps. Called after all components are created and initialized.
virtual void	printOutputAt (FILE file, TimeStep tStep)
virtual Interface *	giveInterface (InterfaceType t)
std::string	errorInfo (const char *func) const
	Returns string for prepending output (used by error reporting macros).

Private Attributes
std::string	inputfile

Static Private Attributes
static int	n = 1

Additional Inherited Members
Protected Attributes inherited from oofem::Material
Dictionary	propertyDictionary
double	castingTime
int	preCastingTimeMat
	Material existing before casting time - optional parameter, zero by default.
Protected Attributes inherited from oofem::FEMComponent
int	number
	Component number.
Domain *	domain
	Link to domain object, useful for communicating with other FEM components.

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 130 of file fe2fluidmaterial.h.

Constructor & Destructor Documentation

◆ FE2FluidMaterial()

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 142 of file fe2fluidmaterial.h.

References oofem::FluidDynamicMaterial::FluidDynamicMaterial(), and n.

Member Function Documentation

◆ checkConsistency()

int oofem::FE2FluidMaterial::checkConsistency ( void )

overridevirtual

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

Definition at line 309 of file fe2fluidmaterial.C.

◆ computeDeviatoricStress3D() [1/2]

std::pair< FloatArrayF< 6 >, double > oofem::FE2FluidMaterial::computeDeviatoricStress3D	(	const FloatArrayF< 6 > &	eps,
		double	pressure,
		GaussPoint *	gp,
		TimeStep *	tStep ) const

overridevirtual

Computes the deviatoric stress vector and volumetric strain rate from given deviatoric strain and pressure. Defaults to incompressible function.

Parameters

gp	Integration point.
eps	Strain-rate.
pressure	Pressure.
tStep	Time step.

Returns: Deviatoric stress and volumetric strain-rate (minus the volumetric part of the eps argument).

Reimplemented from oofem::FluidDynamicMaterial.

Definition at line 75 of file fe2fluidmaterial.C.

References oofem::MixedGradientPressureBC::computeFields(), oofem::FE2FluidMaterialStatus::giveBC(), oofem::EngngModel::giveCurrentStep(), oofem::FE2FluidMaterialStatus::giveRVE(), oofem::Material::giveStatus(), 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().

Referenced by computeDeviatoricStress3D().

◆ computeDeviatoricStress3D() [2/2]

FloatArrayF< 6 > oofem::FE2FluidMaterial::computeDeviatoricStress3D	(	const FloatArrayF< 6 > &	eps,
		GaussPoint *	gp,
		TimeStep *	tStep ) const

overridevirtual

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::FloatArrayF< N >::at(), computeDeviatoricStress3D(), and OOFEM_ERROR.

◆ computeTangent3D()

FloatMatrixF< 6, 6 > oofem::FE2FluidMaterial::computeTangent3D	(	MatResponseMode	mode,
		GaussPoint *	gp,
		TimeStep *	tStep ) const

overridevirtual

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

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

◆ computeTangents3D()

FluidDynamicMaterial::Tangents< 6 > oofem::FE2FluidMaterial::computeTangents3D	(	MatResponseMode	mode,
		GaussPoint *	gp,
		TimeStep *	tStep ) const

overridevirtual

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

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}\)
Parameters

mode Mode of result.

gp Integration point.

tStep Time step.

Reimplemented from oofem::FluidDynamicMaterial.

Definition at line 114 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::Material::giveStatus(), oofem::FE2FluidMaterialStatus::giveVolumetricDeviatoricTangent(), oofem::FE2FluidMaterialStatus::giveVolumetricPressureTangent(), oofem::norm(), OOFEM_ERROR, 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().

◆ CreateStatus()

std::unique_ptr< MaterialStatus > oofem::FE2FluidMaterial::CreateStatus ( GaussPoint * gp ) const

overridevirtual

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

References oofem::FEMComponent::domain, inputfile, and n.

◆ giveClassName()

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

inlineoverridevirtual

Returns: Class name of the receiver.

Implements oofem::FEMComponent.

Definition at line 161 of file fe2fluidmaterial.h.

◆ giveEffectiveViscosity()

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

overridevirtual

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

References oofem::FE2FluidMaterialStatus::computeTangents(), oofem::FE2FluidMaterialStatus::giveDeviatoricTangent(), oofem::GaussPoint::giveMaterialMode(), and oofem::Material::giveStatus().

◆ giveInputRecord()

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

overridevirtual

Setups the input record string of receiver.

Parameters

input Dynamic input record to be filled by receiver.

Reimplemented from oofem::FEMComponent.

Definition at line 253 of file fe2fluidmaterial.C.

References _IFT_FE2FluidMaterial_fileName, inputfile, and oofem::DynamicInputRecord::setField().

◆ giveInputRecordName()

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

inlineoverridevirtual

Returns: Input record name of the receiver.

Implements oofem::FEMComponent.

Definition at line 162 of file fe2fluidmaterial.h.

References _IFT_FE2FluidMaterial_Name.

◆ giveIPValue()

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

overridevirtual

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

Definition at line 260 of file fe2fluidmaterial.C.

References DEBUG_ERR, oofem::FE2FluidMaterialStatus::givePressure(), oofem::Material::giveStatus(), oofem::FE2FluidMaterialStatus::giveVOFFraction(), oofem::FE2FluidMaterialStatus::giveVolumetricPressureTangent(), oofem::norm(), OOFEM_ERROR, and oofem::FloatArray::zero().

◆ initializeFrom()

void oofem::FE2FluidMaterial::initializeFrom ( InputRecord & ir )

overridevirtual

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. Note that initializeFrom may be called mutiple times.

See also: IR_GIVE_FIELD; IR_GIVE_OPTIONAL_FIELD

Parameters

ir	Input record to initialize from.
priority	Priority of the input record. This is used to determine the order of initialization

Reimplemented from oofem::FEMComponent.

Definition at line 247 of file fe2fluidmaterial.C.

References _IFT_FE2FluidMaterial_fileName, inputfile, and IR_GIVE_FIELD.

Member Data Documentation

◆ inputfile

std::string oofem::FE2FluidMaterial::inputfile

private

Definition at line 133 of file fe2fluidmaterial.h.

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

◆ n

int oofem::FE2FluidMaterial::n = 1

staticprivate

Definition at line 134 of file fe2fluidmaterial.h.

Referenced by CreateStatus(), and FE2FluidMaterial().

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

◆ FE2FluidMaterial()

Member Function Documentation

◆ checkConsistency()

◆ computeDeviatoricStress3D() [1/2]

◆ computeDeviatoricStress3D() [2/2]

◆ computeTangent3D()

◆ computeTangents3D()

◆ CreateStatus()

◆ giveClassName()

◆ giveEffectiveViscosity()

◆ giveInputRecord()

◆ giveInputRecordName()

◆ giveIPValue()

◆ initializeFrom()

Member Data Documentation

◆ inputfile

◆ n