#include <fluidmaterialevaluator.h>

Inheritance diagram for oofem::FluidMaterialEvaluator:

Collaboration diagram for oofem::FluidMaterialEvaluator:

Public Member Functions
	FluidMaterialEvaluator (int i, EngngModel *_master=nullptr)
void	initializeFrom (InputRecord &ir) override
void	solveYourself () override
int	checkConsistency () override
void	doStepOutput (TimeStep *tStep) override
TimeStep *	giveNextStep () override
	Returns next time step (next to current step) of receiver.
const char *	giveClassName () const override
	Returns class name of the receiver.
const char *	giveInputRecordName () const
	EngngModel (int i, EngngModel *_master=NULL)
virtual	~EngngModel ()
	Destructor.
	EngngModel (const EngngModel &)=delete
EngngModel &	operator= (const EngngModel &)=delete
Domain *	giveDomain (int n)
void	setDomain (int i, Domain *ptr, bool iDeallocateOld=true)
int	giveNumberOfDomains ()
	Returns number of domains in problem.
const std::string &	giveDescription () const
const time_t &	giveStartTime ()
bool	giveSuppressOutput () const
virtual ErrorEstimator *	giveDomainErrorEstimator (int n)
virtual MaterialInterface *	giveMaterialInterface (int n)
void	setNumberOfEquations (int id, int neq)
FILE *	giveOutputStream ()
	Returns file descriptor of output file.
std::string	giveOutputBaseFileName ()
std::string	giveReferenceFileName ()
void	letOutputBaseFileNameBe (const std ::string &src)
ContextOutputMode	giveContextOutputMode () const
int	giveContextOutputStep () const
void	setContextOutputMode (ContextOutputMode contextMode)
void	setUDContextOutputMode (int cStep)
double	giveDeltaT ()
	Returns time step size from the time step controlelr.
void	setDeltaT (double dT)
	Returns time step size through the time step controlelr.
void	setProblemMode (problemMode pmode)
void	setParallelMode (bool newParallelFlag)
problemMode	giveProblemMode () const
	Returns domain mode.
void	setProblemScale (problemScale pscale)
problemScale	giveProblemScale () const
	Returns scale in multiscale simulation.
virtual void	setRenumberFlag ()
	Sets the renumber flag to true.
virtual void	resetRenumberFlag ()
	Sets the renumber flag to false.
double	giveSolutionStepTime ()
void	giveAnalysisTime (int &rhrs, int &rmin, int &rsec, int &uhrs, int &umin, int &usec)
void	terminateAnalysis ()
virtual void	restartYourself (TimeStep *tS)
virtual void	solveYourselfAt (TimeStep *tStep)
virtual void	terminate (TimeStep *tStep)
void	saveStepContext (TimeStep *tStep, ContextMode mode)
virtual void	updateYourself (TimeStep *tStep)
virtual void	initializeYourself (TimeStep *tStep)
virtual int	initializeAdaptive (int tStepNumber)
virtual int	giveNumberOfDomainEquations (int di, const UnknownNumberingScheme &num)
virtual double	giveUnknownComponent (ValueModeType, TimeStep , Domain , Dof *)
virtual FieldPtr	giveField (FieldType key, TimeStep *)
virtual FieldPtr	giveField (InternalStateType key, TimeStep *)
EngngModel *	giveMasterEngngModel ()
	Returns the master engnmodel.
virtual double	giveLoadLevel ()
	Returns the current load level.
virtual double	giveEigenValue (int eigNum)
	Only relevant for eigen value analysis. Otherwise returns zero.
virtual void	setActiveVector (int i)
	Only relevant for eigen value analysis. Otherwise does noting.
int	updateSharedDofManagers (FloatArray &answer, const UnknownNumberingScheme &s, int ExchangeTag)
int	exchangeRemoteElementData (int ExchangeTag)
virtual int	giveCurrentNumberOfIterations ()
MPI_Comm	giveParallelComm ()
	Returns the communication object of reciever.
int	packRemoteElementData (ProcessCommunicator &processComm)
int	unpackRemoteElementData (ProcessCommunicator &processComm)
int	packDofManagers (ArrayWithNumbering *src, ProcessCommunicator &processComm)
int	unpackDofManagers (ArrayWithNumbering *dest, ProcessCommunicator &processComm)
ProblemCommunicator *	giveProblemCommunicator (EngngModelCommType t)
void	initializeCommMaps (bool forceInit=false)
virtual int	instanciateYourself (DataReader &dr, InputRecord &ir, const char outFileName, const char desc)
void	Instanciate_init ()
int	instanciateDomains (DataReader &dr)
	Instanciate problem domains by calling their instanciateYourself() service.
int	instanciateMetaSteps (DataReader &dr)
	Instanciate problem meta steps by calling their instanciateYourself() service.
virtual int	instanciateDefaultMetaStep (InputRecord &ir)
	Instanciate default metastep, if nmsteps is zero.
virtual void	updateAttributes (MetaStep *mStep)
void	initMetaStepAttributes (MetaStep *mStep)
virtual void	saveContext (DataStream &stream, ContextMode mode)
virtual void	restoreContext (DataStream &stream, ContextMode mode)
virtual void	updateDomainLinks ()
MetaStep *	giveCurrentMetaStep ()
	Returns current meta step.
virtual TimeStep *	giveCurrentStep (bool force=false)
virtual void	adaptTimeStep (double nIter)
virtual TimeStep *	givePreviousStep (bool force=false)
virtual void	preInitializeNextStep ()
	Does a pre-initialization of the next time step (implement if necessarry).
virtual TimeStep *	giveSolutionStepWhenIcApply (bool force=false)
virtual int	giveNumberOfFirstStep (bool force=false)
int	giveNumberOfMetaSteps ()
	Return number of meta steps.
MetaStep *	giveMetaStep (int i)
	Returns the i-th meta step.
int	giveNumberOfSteps (bool force=false)
virtual double	giveEndOfTimeOfInterest ()
	Returns end of time interest (time corresponding to end of time integration).
int	giveNumberOfTimeStepWhenIcApply ()
	Returns the time step number, when initial conditions should apply.
virtual NumericalMethod *	giveNumericalMethod (MetaStep *mStep)
	Returns reference to receiver's numerical method.
ExportModuleManager *	giveExportModuleManager ()
	Returns receiver's export module manager.
EngngModelTimer *	giveTimer ()
	Returns reference to receiver timer (EngngModelTimer).
virtual double	giveInitialTime ()
	return time at the begining of analysis
virtual double	giveFinalTime ()
virtual int	giveNewEquationNumber (int domain, DofIDItem)
virtual int	giveNewPrescribedEquationNumber (int domain, DofIDItem)
std::string	giveContextFileName (int tStepNumber, int stepVersion) const
std::string	giveDomainFileName (int domainNum, int domainSerNum) const
virtual void	updateComponent (TimeStep tStep, NumericalCmpn cmpn, Domain d)
virtual void	initForNewIteration (Domain d, TimeStep tStep, int iterationNumber, const FloatArray &solution)
virtual void	updateSolution (FloatArray &solutionVector, TimeStep tStep, Domain d)
virtual void	updateInternalRHS (FloatArray &answer, TimeStep tStep, Domain d, FloatArray *eNorm)
virtual void	updateMatrix (SparseMtrx &mat, TimeStep tStep, Domain d)
virtual void	initStepIncrements ()
virtual int	forceEquationNumbering (int i)
virtual int	forceEquationNumbering ()
virtual int	requiresUnknownsDictionaryUpdate ()
virtual bool	requiresEquationRenumbering (TimeStep *tStep)
virtual void	updateDofUnknownsDictionary (DofManager , TimeStep )
virtual int	giveUnknownDictHashIndx (ValueModeType mode, TimeStep *tStep)
virtual bool	newDofHandling ()
virtual ParallelContext *	giveParallelContext (int n)
virtual void	initParallelContexts ()
virtual void	assemble (SparseMtrx &answer, TimeStep tStep, const MatrixAssembler &ma, const UnknownNumberingScheme &s, Domain domain)
virtual void	assemble (SparseMtrx &answer, TimeStep tStep, const MatrixAssembler &ma, const UnknownNumberingScheme &r_s, const UnknownNumberingScheme &c_s, Domain domain)
void	assembleVector (FloatArray &answer, TimeStep tStep, const VectorAssembler &va, ValueModeType mode, const UnknownNumberingScheme &s, Domain domain, FloatArray *eNorms=NULL)
void	assembleVectorFromDofManagers (FloatArray &answer, TimeStep tStep, const VectorAssembler &va, ValueModeType mode, const UnknownNumberingScheme &s, Domain domain, FloatArray *eNorms=NULL)
void	assembleVectorFromElements (FloatArray &answer, TimeStep tStep, const VectorAssembler &va, ValueModeType mode, const UnknownNumberingScheme &s, Domain domain, FloatArray *eNorms=NULL)
void	assembleVectorFromBC (FloatArray &answer, TimeStep tStep, const VectorAssembler &va, ValueModeType mode, const UnknownNumberingScheme &s, Domain domain, FloatArray *eNorms=NULL)
void	assembleExtrapolatedForces (FloatArray &answer, TimeStep tStep, CharType type, Domain domain)
void	assemblePrescribedExtrapolatedForces (FloatArray &answer, TimeStep tStep, CharType type, Domain domain)
virtual int	checkProblemConsistency ()
virtual void	init ()
virtual void	postInitialize ()
virtual void	printOutputAt (FILE file, TimeStep tStep)
virtual void	printOutputAt (FILE file, TimeStep tStep, const IntArray &nodeSets, const IntArray &elementSets)
void	outputNodes (FILE file, Domain &domain, TimeStep tStep, int setNum)
void	outputElements (FILE file, Domain &domain, TimeStep tStep, int setNum)
void	printYourself ()
	Prints state of receiver. Useful for debugging.
virtual void	printDofOutputAt (FILE stream, Dof iDof, TimeStep *tStep)
virtual int	useNonlocalStiffnessOption ()
	Returns nonzero if nonlocal stiffness option activated.
bool	isParallel () const
	Returns true if receiver in parallel mode.
int	giveRank () const
	Returns domain rank in a group of collaborating processes (0..groupSize-1).
int	giveNumberOfProcesses () const
	Returns the number of collaborating processes.
virtual fMode	giveFormulation ()
EngngModelContext *	giveContext ()
	Context requesting service.
virtual int	giveNumberOfSlaveProblems ()
	Returns number of slave problems.
virtual EngngModel *	giveSlaveProblem (int i)
	Returns i-th slave problem.
virtual bool	giveEquationScalingFlag ()
	Returns the Equation scaling flag, which is used to indicate that governing equation(s) are scaled, or non-dimensionalized.
virtual double	giveVariableScale (VarScaleType varId)
	Returns the scale factor for given variable type.
virtual int	estimateMaxPackSize (IntArray &commMap, DataStream &buff, int packUnpackType)
virtual void	balanceLoad (TimeStep *tStep)
virtual LoadBalancer *	giveLoadBalancer ()
virtual LoadBalancerMonitor *	giveLoadBalancerMonitor ()
void	initParallel ()
	Request domain rank and problem size.
EngngModel *	giveEngngModel ()
	Returns reference to itself -> required by communicator.h.
virtual bool	isElementActivated (int elemNum)
virtual bool	isElementActivated (Element *e)
TimeStepController *	giveTimeStepController ()
	Returns the time step controller.
virtual void	drawYourself (oofegGraphicContext &gc)
virtual void	drawElements (oofegGraphicContext &gc)
virtual void	drawNodes (oofegGraphicContext &gc)
virtual void	showSparseMtrxStructure (int type, oofegGraphicContext &gc, TimeStep *tStep)
std::string	errorInfo (const char *func) const
	Returns string for prepending output (used by error reporting macros).

Protected Attributes
double	deltaT = 1.
	Time increment.
int	ndim = 3
IntArray	cmpntFunctions
	Number of spatial dimensions.
int	volFunction = 0
	Time functions controlling each component of the deviatoric part of the stress.
IntArray	sControl
	Time function controlling the volumetric/pressure part.
IntArray	eControl
bool	pressureControl = false
IntArray	vars
std::vector< std ::unique_ptr< GaussPoint > >	gps
std::ofstream	outfile
Protected Attributes inherited from oofem::EngngModel
int	ndomains
	Number of receiver domains.
std ::vector< std ::unique_ptr< Domain > >	domainList
	List of problem domains.
int	numberOfSteps
	Total number of time steps.
int	numberOfEquations
	Total number of equation in current time step.
int	numberOfPrescribedEquations
	Total number or prescribed equations in current time step.
IntArray	domainNeqs
	Number of equations per domain.
IntArray	domainPrescribedNeqs
	Number of prescribed equations per domain.
bool	renumberFlag
	Renumbering flag (renumbers equations after each step, necessary if Dirichlet BCs change).
bool	profileOpt
	Profile optimized numbering flag (using Sloan's algorithm).
int	equationNumberingCompleted
	Equation numbering completed flag.
int	nMetaSteps
	Number of meta steps.
std ::vector< MetaStep >	metaStepList
	List of problem metasteps.
std ::unique_ptr< TimeStep >	stepWhenIcApply
	Solution step when IC (initial conditions) apply.
std ::unique_ptr< TimeStep >	currentStep
	Current time step.
std ::unique_ptr< TimeStep >	previousStep
	Previous time step.
int	number
	Receivers id.
std::string	dataOutputFileName
	Path to output stream.
std::string	coreOutputFileName
	String with core output file name.
FILE *	outputStream
	Output stream.
std::string	referenceFileName
	String with reference file name.
ContextOutputMode	contextOutputMode
	Domain context output mode.
int	contextOutputStep
ExportModuleManager	exportModuleManager
	Export module manager.
InitModuleManager	initModuleManager
	Initialization module manager.
MonitorManager	monitorManager
	Monitor manager.
problemMode	pMode
	Domain mode.
problemScale	pScale
	Multiscale mode.
time_t	startTime
	Solution start time.
EngngModel *	master
	Master e-model; if defined receiver is in maintained (slave) mode.
EngngModelContext *	context
	Context.
EngngModelTimer	timer
	E-model timer.
int	parallelFlag
	Flag indicating that the receiver runs in parallel.
enum fMode	nonLinFormulation
	Type of non linear formulation (total or updated formulation).
std::unique_ptr< ErrorEstimator >	defaultErrEstimator
	Error estimator. Useful for adaptivity, or simply printing errors output.
std::unique_ptr< TimeStepController >	timeStepController
	Time Step controller is responsible for collecting data from analysis, elements, and materials, and select the appropriate timestep size for the next step, or reduce the step in case of convergence problems.
int	rank
	Domain rank in a group of collaborating processes (0..groupSize-1).
int	numProcs
	Total number of collaborating processes.
int	nonlocalExt
	Flag indicating if nonlocal extension active, which will cause data to be sent between shared elements before computing the internal forces.
char	processor_name [PROCESSOR_NAME_LENGTH]
	Processor name.
MPI_Comm	comm
	Communication object for this engineering model.
std::unique_ptr< LoadBalancer >	lb
	Load Balancer.
std::unique_ptr< LoadBalancerMonitor >	lbm
bool	loadBalancingFlag
	If set to true, load balancing is active.
bool	force_load_rebalance_in_first_step
	Debug flag forcing load balancing after first step.
CommunicatorBuff *	commBuff
	Common Communicator buffer.
ProblemCommunicator *	communicator
	Communicator.
ProblemCommunicator *	nonlocCommunicator
	NonLocal Communicator. Necessary when nonlocal constitutive models are used.
std ::vector< ParallelContext >	parallelContextList
	List where parallel contexts are stored.
bool	suppressOutput
	Flag for suppressing output to file.
std::string	simulationDescription

Additional Inherited Members
Public Types inherited from oofem::EngngModel
enum	EngngModel_UpdateMode { EngngModel_Unknown_Mode , EngngModel_SUMM_Mode , EngngModel_SET_Mode }
enum	EngngModelCommType { PC_default , PC_nonlocal }
enum	InitialGuess { IG_None = 0 , IG_Tangent = 1 }
Protected Types inherited from oofem::EngngModel
enum	{ InternalForcesExchangeTag , MassExchangeTag , LoadExchangeTag , ReactionExchangeTag , RemoteElementExchangeTag }
	Message tags. More...
virtual void	packMigratingData (TimeStep *tStep)
virtual void	unpackMigratingData (TimeStep *tStep)

Detailed Description

For testing material behavior, particularly useful for multiscale modeling where one can test a single RVE. The deviatoric and volumetric parts are split. No nodes or elements are used.

Note: The user must take care to ensure that a deviatoric stress is possible, i.e. is all components are controlled for the stress, the time functions must evaluate to something deviatoric (e.g. all zeros).

This model will output data in its own way since it does not contain any actual FE-results. It will not call the output manager at all.

Author: Mikael Öhman

Definition at line 69 of file fluidmaterialevaluator.h.

Constructor & Destructor Documentation

◆ FluidMaterialEvaluator()

oofem::FluidMaterialEvaluator::FluidMaterialEvaluator	(	int	i,
		EngngModel *	_master = nullptr )

Definition at line 49 of file fluidmaterialevaluator.C.

References oofem::EngngModel::EngngModel(), and oofem::EngngModel::ndomains.

Member Function Documentation

◆ checkConsistency()

int oofem::FluidMaterialEvaluator::checkConsistency ( )

overridevirtual

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

Returns: Nonzero if receiver check is o.k.

Reimplemented from oofem::EngngModel.

Definition at line 229 of file fluidmaterialevaluator.C.

References oofem::EngngModel::giveDomain(), oofem::Domain::giveMaterial(), and oofem::Domain::giveNumberOfMaterialModels().

◆ doStepOutput()

void oofem::FluidMaterialEvaluator::doStepOutput ( TimeStep * tStep )

overridevirtual

Prints the ouput of the solution step (using virtual this->printOutputAtservice) to the stream detemined using this->giveOutputStream() method and calls exportModuleManager to do output.

Reimplemented from oofem::EngngModel.

Definition at line 241 of file fluidmaterialevaluator.C.

References oofem::__InternalStateTypeToString(), oofem::EngngModel::giveDomain(), oofem::TimeStep::giveIntrinsicTime(), oofem::FluidDynamicMaterial::giveIPValue(), oofem::Domain::giveMaterial(), oofem::Domain::giveNumberOfMaterialModels(), gps, oofem::TimeStep::isTheFirstStep(), outfile, and vars.

Referenced by solveYourself().

◆ giveClassName()

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

inlineoverridevirtual

Returns class name of the receiver.

Implements oofem::EngngModel.

Definition at line 97 of file fluidmaterialevaluator.h.

◆ giveInputRecordName()

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

inline

Definition at line 98 of file fluidmaterialevaluator.h.

References _IFT_FluidMaterialEvaluator_Name.

◆ giveNextStep()

TimeStep * oofem::FluidMaterialEvaluator::giveNextStep ( )

overridevirtual

Returns next time step (next to current step) of receiver.

Reimplemented from oofem::EngngModel.

Definition at line 267 of file fluidmaterialevaluator.C.

References oofem::EngngModel::currentStep, deltaT, oofem::EngngModel::giveNumberOfTimeStepWhenIcApply(), and oofem::EngngModel::previousStep.

Referenced by solveYourself().

◆ initializeFrom()

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

overridevirtual

Initializes receiver according to object description in input reader. InitString 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.

Reimplemented from oofem::EngngModel.

Definition at line 55 of file fluidmaterialevaluator.C.

◆ solveYourself()

void oofem::FluidMaterialEvaluator::solveYourself ( )

overridevirtual

Starts solution process. Implementation should invoke for each time step solveYourselfAt function with time step as parameter. Time steps are created using giveNextStep function (this will set current time step to newly created, and updates previous step).

Todo: We should ask the material model for this initial guess.

Todo: More flexible control of convergence needed. Something relative?

Reimplemented from oofem::EngngModel.

Definition at line 80 of file fluidmaterialevaluator.C.

References oofem::FloatArray::at(), oofem::FloatMatrix::at(), oofem::FloatMatrix::beSubMatrixOf(), cmpntFunctions, oofem::FluidDynamicMaterial::computeDeviatoricStress2D(), oofem::FluidDynamicMaterial::computeDeviatoricStress3D(), oofem::FloatMatrix::computeFrobeniusNorm(), oofem::FloatArray::computeNorm(), oofem::FluidDynamicMaterial::computeTangents2D(), oofem::FluidDynamicMaterial::computeTangents3D(), doStepOutput(), eControl, oofem::Function::evaluateAtTime(), oofem::FluidDynamicMaterialStatus::giveDeviatoricStrainRateVector(), oofem::EngngModel::giveDomain(), oofem::Domain::giveFunction(), oofem::TimeStep::giveIntrinsicTime(), oofem::Domain::giveMaterial(), giveNextStep(), oofem::Domain::giveNumberOfMaterialModels(), oofem::EngngModel::giveOutputBaseFileName(), oofem::FloatArray::giveSize(), oofem::Material::giveStatus(), gps, oofem::FluidDynamicMaterialStatus::letDeviatoricStrainRateVectorBe(), ndim, oofem::norm(), oofem::EngngModel::numberOfSteps, OOFEM_ERROR, OOFEM_LOG_RELEVANT, OOFEM_WARNING, outfile, pressureControl, sControl, oofem::FloatMatrix::solveForRhs(), oofem::EngngModel::timer, oofem::tolerance, oofem::GaussPoint::updateYourself(), volFunction, and oofem::FloatArray::zero().

Member Data Documentation

◆ cmpntFunctions

IntArray oofem::FluidMaterialEvaluator::cmpntFunctions

protected

Number of spatial dimensions.

Definition at line 75 of file fluidmaterialevaluator.h.

Referenced by initializeFrom(), and solveYourself().

◆ deltaT

double oofem::FluidMaterialEvaluator::deltaT = 1.

protected

Time increment.

Definition at line 72 of file fluidmaterialevaluator.h.

Referenced by giveNextStep(), and initializeFrom().

◆ eControl

IntArray oofem::FluidMaterialEvaluator::eControl

protected

Definition at line 77 of file fluidmaterialevaluator.h.

Referenced by initializeFrom(), and solveYourself().

◆ gps

std::vector< std :: unique_ptr< GaussPoint > > oofem::FluidMaterialEvaluator::gps

protected

Definition at line 82 of file fluidmaterialevaluator.h.

Referenced by doStepOutput(), and solveYourself().

◆ ndim

int oofem::FluidMaterialEvaluator::ndim = 3

protected

Definition at line 74 of file fluidmaterialevaluator.h.

Referenced by initializeFrom(), and solveYourself().

◆ outfile

std::ofstream oofem::FluidMaterialEvaluator::outfile

protected

Definition at line 84 of file fluidmaterialevaluator.h.

Referenced by doStepOutput(), and solveYourself().

◆ pressureControl

bool oofem::FluidMaterialEvaluator::pressureControl = false

protected

Definition at line 78 of file fluidmaterialevaluator.h.

Referenced by initializeFrom(), and solveYourself().

◆ sControl

IntArray oofem::FluidMaterialEvaluator::sControl

protected

Time function controlling the volumetric/pressure part.

Definition at line 77 of file fluidmaterialevaluator.h.

Referenced by initializeFrom(), and solveYourself().

◆ vars

IntArray oofem::FluidMaterialEvaluator::vars

protected

Definition at line 80 of file fluidmaterialevaluator.h.

Referenced by doStepOutput(), and initializeFrom().

◆ volFunction

int oofem::FluidMaterialEvaluator::volFunction = 0

protected

Time functions controlling each component of the deviatoric part of the stress.

Definition at line 76 of file fluidmaterialevaluator.h.

Referenced by initializeFrom(), and solveYourself().

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

Public Member Functions

Protected Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ FluidMaterialEvaluator()

Member Function Documentation

◆ checkConsistency()

◆ doStepOutput()

◆ giveClassName()

◆ giveInputRecordName()

◆ giveNextStep()

◆ initializeFrom()

◆ solveYourself()

Member Data Documentation

◆ cmpntFunctions

◆ deltaT

◆ eControl

◆ gps

◆ ndim

◆ outfile

◆ pressureControl

◆ sControl

◆ vars

◆ volFunction