OOFEM 3.0
Loading...
Searching...
No Matches
oofem::StationaryMPMSProblem Class Reference

#include <mpmsemodels.h>

Inheritance diagram for oofem::StationaryMPMSProblem:
Collaboration diagram for oofem::StationaryMPMSProblem:

Public Member Functions

 StationaryMPMSProblem (int i, EngngModel *_master)
void initializeFrom (InputRecord &ir) override
void postInitialize () override
void solveYourselfAt (TimeStep *tStep) override
void updateComponent (TimeStep *tStep, NumericalCmpn cmpn, Domain *d) override
void updateSolution (FloatArray &solutionVector, TimeStep *tStep, Domain *d) override
TimeStepgiveNextStep () override
 Returns next time step (next to current step) of receiver.
NumericalMethodgiveNumericalMethod (MetaStep *mStep) override
 Returns reference to receiver's numerical method.
double giveUnknownComponent (ValueModeType mode, TimeStep *tStep, Domain *d, Dof *dof) override
const char * giveInputRecordName () const
const char * giveClassName () const override
 Returns class name of the receiver.
fMode giveFormulation () override
 EngngModel (int i, EngngModel *_master=NULL)
virtual ~EngngModel ()
 Destructor.
 EngngModel (const EngngModel &)=delete
EngngModeloperator= (const EngngModel &)=delete
DomaingiveDomain (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 ErrorEstimatorgiveDomainErrorEstimator (int n)
virtual MaterialInterfacegiveMaterialInterface (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 solveYourself ()
virtual void restartYourself (TimeStep *tS)
virtual void terminate (TimeStep *tStep)
virtual void doStepOutput (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 FieldPtr giveField (FieldType key, TimeStep *)
virtual FieldPtr giveField (InternalStateType key, TimeStep *)
EngngModelgiveMasterEngngModel ()
 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)
ProblemCommunicatorgiveProblemCommunicator (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 ()
MetaStepgiveCurrentMetaStep ()
 Returns current meta step.
virtual TimeStepgiveCurrentStep (bool force=false)
virtual void adaptTimeStep (double nIter)
virtual TimeStepgivePreviousStep (bool force=false)
virtual void preInitializeNextStep ()
 Does a pre-initialization of the next time step (implement if necessarry).
virtual TimeStepgiveSolutionStepWhenIcApply (bool force=false)
virtual int giveNumberOfFirstStep (bool force=false)
int giveNumberOfMetaSteps ()
 Return number of meta steps.
MetaStepgiveMetaStep (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.
ExportModuleManagergiveExportModuleManager ()
 Returns receiver's export module manager.
EngngModelTimergiveTimer ()
 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 initForNewIteration (Domain *d, TimeStep *tStep, int iterationNumber, const FloatArray &solution)
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 ParallelContextgiveParallelContext (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 checkConsistency ()
virtual int checkProblemConsistency ()
virtual void init ()
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.
EngngModelContextgiveContext ()
 Context requesting service.
virtual int giveNumberOfSlaveProblems ()
 Returns number of slave problems.
virtual EngngModelgiveSlaveProblem (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 LoadBalancergiveLoadBalancer ()
virtual LoadBalancerMonitorgiveLoadBalancerMonitor ()
void initParallel ()
 Request domain rank and problem size.
EngngModelgiveEngngModel ()
 Returns reference to itself -> required by communicator.h.
virtual bool isElementActivated (int elemNum)
virtual bool isElementActivated (Element *e)
TimeStepControllergiveTimeStepController ()
 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

SparseMtrxType sparseMtrxType = SMT_Skyline
LinSystSolverType solverType = ST_Direct
std ::unique_ptr< PrimaryFieldunknownsField
 This field stores solution vector. For fixed size of problem, the PrimaryField is used, for growing/decreasing size, DofDistributedPrimaryField applies.
std ::unique_ptr< SparseMtrxeffectiveMatrix
FloatArray solutionVector
FloatArray residualVector
FloatArray eNorm
std::unique_ptr< SparseNonLinearSystemNMnMethod
 Numerical method used to solve the problem.
bool keepTangent = false
IntArray lhsIntegrals
IntArray rhsIntegrals
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< MetaStepmetaStepList
 List of problem metasteps.
std ::unique_ptr< TimeStepstepWhenIcApply
 Solution step when IC (initial conditions) apply.
std ::unique_ptr< TimeStepcurrentStep
 Current time step.
std ::unique_ptr< TimeSteppreviousStep
 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.
EngngModelmaster
 Master e-model; if defined receiver is in maintained (slave) mode.
EngngModelContextcontext
 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< ErrorEstimatordefaultErrEstimator
 Error estimator. Useful for adaptivity, or simply printing errors output.
std::unique_ptr< TimeStepControllertimeStepController
 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< LoadBalancerlb
 Load Balancer.
std::unique_ptr< LoadBalancerMonitorlbm
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.
CommunicatorBuffcommBuff
 Common Communicator buffer.
ProblemCommunicatorcommunicator
 Communicator.
ProblemCommunicatornonlocCommunicator
 NonLocal Communicator. Necessary when nonlocal constitutive models are used.
std ::vector< ParallelContextparallelContextList
 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

Definition at line 55 of file mpmsemodels.h.

Constructor & Destructor Documentation

◆ StationaryMPMSProblem()

oofem::StationaryMPMSProblem::StationaryMPMSProblem ( int i,
EngngModel * _master )
inline

Member Function Documentation

◆ giveClassName()

const char * oofem::StationaryMPMSProblem::giveClassName ( ) const
inlineoverridevirtual

Returns class name of the receiver.

Implements oofem::EngngModel.

Definition at line 230 of file mpmsemodels.h.

◆ giveFormulation()

fMode oofem::StationaryMPMSProblem::giveFormulation ( )
inlineoverridevirtual

Indicates type of non linear computation (total or updated formulation). This is used for example on Nodal level to update coordinates if updated formulation is done, or on element level, when non linear contributions are computed.

Reimplemented from oofem::EngngModel.

Definition at line 231 of file mpmsemodels.h.

References oofem::TL.

◆ giveInputRecordName()

const char * oofem::StationaryMPMSProblem::giveInputRecordName ( ) const
inline

Definition at line 229 of file mpmsemodels.h.

References _IFT_StationaryMPMSProblem_Name.

◆ giveNextStep()

TimeStep * oofem::StationaryMPMSProblem::giveNextStep ( )
inlineoverridevirtual

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

Reimplemented from oofem::EngngModel.

Definition at line 191 of file mpmsemodels.h.

References oofem::EngngModel::currentStep, oofem::EngngModel::giveNumberOfFirstStep(), and oofem::EngngModel::previousStep.

◆ giveNumericalMethod()

NumericalMethod * oofem::StationaryMPMSProblem::giveNumericalMethod ( MetaStep * mStep)
inlineoverridevirtual

Returns reference to receiver's numerical method.

Reimplemented from oofem::EngngModel.

Definition at line 206 of file mpmsemodels.h.

References oofem::EngngModel::giveDomain(), and nMethod.

Referenced by solveYourselfAt().

◆ giveUnknownComponent()

double oofem::StationaryMPMSProblem::giveUnknownComponent ( ValueModeType ,
TimeStep * ,
Domain * ,
Dof *  )
inlineoverridevirtual

Returns requested unknown. Unknown at give time step is characterized by its type and mode and by its equation number. This function is used by Dofs, when they are requested for their associated unknowns.

See also
Dof::giveUnknown

Reimplemented from oofem::EngngModel.

Definition at line 214 of file mpmsemodels.h.

References OOFEM_ERROR, and unknownsField.

◆ initializeFrom()

void oofem::StationaryMPMSProblem::initializeFrom ( InputRecord & ir)
inlineoverridevirtual

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 81 of file mpmsemodels.h.

References oofem::EngngModel::initializeFrom(), IR_GIVE_FIELD, lhsIntegrals, rhsIntegrals, and unknownsField.

◆ postInitialize()

void oofem::StationaryMPMSProblem::postInitialize ( )
inlineoverridevirtual

Performs post-initialization for all the problem contents (which is called after initializeFrom). Currently, it calls Domain::postInitialize for all problem domains.

Reimplemented from oofem::EngngModel.

Definition at line 91 of file mpmsemodels.h.

References oofem::ConnectivityTable::buildSharedBoundaryEntities(), oofem::Domain::giveConnectivityTable(), oofem::EngngModel::giveDomain(), and oofem::EngngModel::postInitialize().

◆ solveYourselfAt()

void oofem::StationaryMPMSProblem::solveYourselfAt ( TimeStep * tStep)
inlineoverridevirtual

Solves problem for given time step. Should assemble characteristic matrices and vectors if necessary and solve problem using appropriate numerical method. After finishing solution, this->updateYourself function for updating solution state and then this->terminate function (for updating nodal and element values) should be called.

Reimplemented from oofem::EngngModel.

Definition at line 99 of file mpmsemodels.h.

References oofem::Integral::assemble_lhs(), oofem::Integral::assemble_rhs(), oofem::classFactory, effectiveMatrix, oofem::EngngModel::forceEquationNumbering(), oofem::EngngModel::giveCurrentMetaStep(), oofem::TimeStep::giveNumber(), oofem::EngngModel::giveNumberOfDomainEquations(), giveNumericalMethod(), lhsIntegrals, oofem::EngngModel::LoadExchangeTag, OOFEM_ERROR, OOFEM_LOG_INFO, residualVector, rhsIntegrals, solutionVector, sparseMtrxType, unknownsField, and oofem::EngngModel::updateSharedDofManagers().

◆ updateComponent()

◆ updateSolution()

void oofem::StationaryMPMSProblem::updateSolution ( FloatArray & solutionVector,
TimeStep * tStep,
Domain * d )
inlineoverridevirtual

Updates the solution (guess) according to the new values. Callback for nonlinear solvers (e.g. Newton-Raphson), and are called before new internal forces are computed.

Parameters
solutionVectorNew solution.
tStepTime when component is updated.
dDomain.

Reimplemented from oofem::EngngModel.

Definition at line 187 of file mpmsemodels.h.

References solutionVector.

Member Data Documentation

◆ effectiveMatrix

std :: unique_ptr< SparseMtrx > oofem::StationaryMPMSProblem::effectiveMatrix
protected

Definition at line 64 of file mpmsemodels.h.

Referenced by solveYourselfAt(), and updateComponent().

◆ eNorm

FloatArray oofem::StationaryMPMSProblem::eNorm
protected

Definition at line 67 of file mpmsemodels.h.

Referenced by updateComponent().

◆ keepTangent

bool oofem::StationaryMPMSProblem::keepTangent = false
protected

Definition at line 72 of file mpmsemodels.h.

◆ lhsIntegrals

IntArray oofem::StationaryMPMSProblem::lhsIntegrals
protected

Definition at line 75 of file mpmsemodels.h.

Referenced by initializeFrom(), solveYourselfAt(), and updateComponent().

◆ nMethod

std::unique_ptr<SparseNonLinearSystemNM> oofem::StationaryMPMSProblem::nMethod
protected

Numerical method used to solve the problem.

Definition at line 70 of file mpmsemodels.h.

Referenced by giveNumericalMethod(), and StationaryMPMSProblem().

◆ residualVector

FloatArray oofem::StationaryMPMSProblem::residualVector
protected

Definition at line 66 of file mpmsemodels.h.

Referenced by solveYourselfAt(), and updateComponent().

◆ rhsIntegrals

IntArray oofem::StationaryMPMSProblem::rhsIntegrals
protected

Definition at line 76 of file mpmsemodels.h.

Referenced by initializeFrom(), and solveYourselfAt().

◆ solutionVector

FloatArray oofem::StationaryMPMSProblem::solutionVector
protected

Definition at line 65 of file mpmsemodels.h.

Referenced by solveYourselfAt(), and updateSolution().

◆ solverType

LinSystSolverType oofem::StationaryMPMSProblem::solverType = ST_Direct
protected

Definition at line 60 of file mpmsemodels.h.

◆ sparseMtrxType

SparseMtrxType oofem::StationaryMPMSProblem::sparseMtrxType = SMT_Skyline
protected

Definition at line 59 of file mpmsemodels.h.

Referenced by solveYourselfAt().

◆ unknownsField

std :: unique_ptr< PrimaryField > oofem::StationaryMPMSProblem::unknownsField
protected

This field stores solution vector. For fixed size of problem, the PrimaryField is used, for growing/decreasing size, DofDistributedPrimaryField applies.

Definition at line 62 of file mpmsemodels.h.

Referenced by giveUnknownComponent(), initializeFrom(), and solveYourselfAt().


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

This page is part of the OOFEM-3.0 documentation. Copyright Copyright (C) 1994-2025 Borek Patzak Bořek Patzák
Project e-mail: oofem@fsv.cvut.cz
Generated at for OOFEM by doxygen 1.15.0 written by Dimitri van Heesch, © 1997-2011