OOFEM
2.4
OOFEM.org - Object Oriented Finite Element Solver
|
This class implements Adaptive Non-LinearStatic Engineering problem. More...
#include <adaptnlinearstatic.h>
Public Member Functions | |
AdaptiveNonLinearStatic (int i, EngngModel *_master=NULL) | |
virtual | ~AdaptiveNonLinearStatic () |
virtual void | solveYourselfAt (TimeStep *tStep) |
Solves problem for given time step. More... | |
virtual void | updateYourself (TimeStep *tStep) |
Updates internal state after finishing time step. More... | |
virtual IRResultType | initializeFrom (InputRecord *ir) |
Initializes receiver according to object description in input reader. More... | |
virtual double | giveUnknownComponent (ValueModeType mode, TimeStep *tStep, Domain *d, Dof *dof) |
Returns requested unknown. More... | |
virtual double | giveTimeStepLoadLevel (int istep) |
Returns the load level corresponding to given solution step number. More... | |
virtual int | initializeAdaptive (int tStepNumber) |
Initializes the newly generated discretization state according to previous solution. More... | |
virtual int | initializeAdaptiveFrom (EngngModel *sourceProblem) |
Initializes the receiver state according to state of given source problem. More... | |
virtual int | adaptiveRemap (Domain *dNew) |
Remaps the solution state to newly given domain. More... | |
virtual contextIOResultType | saveContext (DataStream &stream, ContextMode mode) |
Stores the state of model to output stream. More... | |
virtual contextIOResultType | restoreContext (DataStream &stream, ContextMode mode) |
Restores the state of model from output stream. More... | |
virtual void | updateDomainLinks () |
Updates domain links after the domains of receiver have changed. More... | |
virtual const char * | giveClassName () const |
Returns class name of the receiver. More... | |
virtual const char * | giveInputRecordName () const |
virtual LoadBalancer * | giveLoadBalancer () |
Returns reference to receiver's load balancer. More... | |
virtual LoadBalancerMonitor * | giveLoadBalancerMonitor () |
Returns reference to receiver's load balancer monitor. More... | |
Public Member Functions inherited from oofem::NonLinearStatic | |
NonLinearStatic (int i, EngngModel *_master=NULL) | |
virtual | ~NonLinearStatic () |
virtual void | solveYourself () |
Starts solution process. More... | |
virtual void | terminate (TimeStep *tStep) |
Terminates the solution of time step. More... | |
virtual void | printOutputAt (FILE *file, TimeStep *tStep) |
Prints output of receiver to output domain stream, for given time step. More... | |
virtual void | updateComponent (TimeStep *tStep, NumericalCmpn, Domain *d) |
Updates components mapped to numerical method if necessary during solution process. More... | |
virtual void | updateAttributes (MetaStep *mStep) |
Update receiver attributes according to step metaStep attributes. More... | |
virtual TimeStep * | giveSolutionStepWhenIcApply (bool force=false) |
Returns the solution step when Initial Conditions (IC) apply. More... | |
virtual TimeStep * | giveNextStep () |
Returns next time step (next to current step) of receiver. More... | |
virtual NumericalMethod * | giveNumericalMethod (MetaStep *mStep) |
Returns reference to receiver's numerical method. More... | |
virtual double | giveLoadLevel () |
Returns the current load level. More... | |
virtual fMode | giveFormulation () |
Indicates type of non linear computation (total or updated formulation). More... | |
virtual int | useNonlocalStiffnessOption () |
Returns nonzero if nonlocal stiffness option activated. More... | |
virtual int | giveUnknownDictHashIndx (ValueModeType mode, TimeStep *tStep) |
For load balancing purposes we store all values so hash is computed from mode value only. More... | |
virtual int | giveCurrentNumberOfIterations () |
Returns number of iterations that was required to reach equilibrium - used for adaptive step length in staggered problem. More... | |
virtual void | showSparseMtrxStructure (int type, oofegGraphicContext &gc, TimeStep *tStep) |
Shows the sparse structure of required matrix, type == 1 stiffness. More... | |
virtual int | estimateMaxPackSize (IntArray &commMap, DataStream &buff, int packUnpackType) |
Determines the space necessary for send/receive buffer. More... | |
Public Member Functions inherited from oofem::LinearStatic | |
LinearStatic (int i, EngngModel *_master=NULL) | |
virtual | ~LinearStatic () |
virtual UnknownNumberingScheme * | giveEquationNumbering () |
Public Member Functions inherited from oofem::StructuralEngngModel | |
StructuralEngngModel (int i, EngngModel *_master=NULL) | |
Creates new StructuralEngngModel with number i, associated to domain d. More... | |
virtual | ~StructuralEngngModel () |
Destructor. More... | |
virtual int | checkConsistency () |
Allows programmer to test some receiver's internal data, before computation begins. More... | |
void | computeReaction (FloatArray &answer, TimeStep *tStep, int di) |
Computes reaction forces. More... | |
void | buildReactionTable (IntArray &restrDofMans, IntArray &restrDofs, IntArray &eqn, TimeStep *tStep, int di) |
Builds the reaction force table. More... | |
Public Member Functions inherited from oofem::EngngModel | |
EngngModel (int i, EngngModel *_master=NULL) | |
Constructor. More... | |
virtual | ~EngngModel () |
Destructor. More... | |
EngngModel (const EngngModel &)=delete | |
EngngModel & | operator= (const EngngModel &)=delete |
Domain * | giveDomain (int n) |
Service for accessing particular problem domain. More... | |
void | setDomain (int i, Domain *ptr, bool iDeallocateOld=true) |
Sets i-th domain of receiver. More... | |
int | giveNumberOfDomains () |
Returns number of domains in problem. More... | |
const std::string & | giveDescription () const |
const time_t & | giveStartTime () |
bool | giveSuppressOutput () const |
virtual ErrorEstimator * | giveDomainErrorEstimator (int n) |
Service for accessing ErrorEstimator corresponding to particular domain. More... | |
virtual MaterialInterface * | giveMaterialInterface (int n) |
Returns material interface representation for given domain. More... | |
void | setNumberOfEquations (int id, int neq) |
FILE * | giveOutputStream () |
Returns file descriptor of output file. More... | |
std::string | giveOutputBaseFileName () |
Returns base output file name to which extensions, like .out .vtu .osf should be added. More... | |
std::string | giveReferenceFileName () |
Returns reference file name. More... | |
void | letOutputBaseFileNameBe (const std::string &src) |
Sets the base output file name. More... | |
ContextOutputMode | giveContextOutputMode () |
Returns domain context output mode. More... | |
int | giveContextOutputStep () |
Returns domain context output step. More... | |
void | setContextOutputMode (ContextOutputMode contextMode) |
Sets context output mode of receiver. More... | |
void | setUDContextOutputMode (int cStep) |
Sets user defined context output mode (it sets contextOutputMode to contextOutputMode), setting contextOutputStep to given value. More... | |
void | setProblemMode (problemMode pmode) |
Sets domain mode to given mode. More... | |
void | setParallelMode (bool newParallelFlag) |
Sets the problem to run in parallel (or not). More... | |
problemMode | giveProblemMode () |
Returns domain mode. More... | |
void | setProblemScale (problemScale pscale) |
Sets scale in multiscale simulation. More... | |
problemScale | giveProblemScale () |
Returns scale in multiscale simulation. More... | |
virtual void | setRenumberFlag () |
Sets the renumber flag to true. More... | |
virtual void | resetRenumberFlag () |
Sets the renumber flag to false. More... | |
double | giveSolutionStepTime () |
Returns the user time of the current simulation step in seconds. More... | |
void | giveAnalysisTime (int &rhrs, int &rmin, int &rsec, int &uhrs, int &umin, int &usec) |
Returns the real and user time for the analysis. More... | |
void | terminateAnalysis () |
Performs analysis termination after finishing analysis. More... | |
virtual void | doStepOutput (TimeStep *tStep) |
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. More... | |
void | saveStepContext (TimeStep *tStep, ContextMode mode) |
Saves context of given solution step, if required (determined using this->giveContextOutputMode() method). More... | |
virtual void | initializeYourself (TimeStep *tStep) |
Provides the opportunity to initialize state variables stored in element integration points according to initial conditions using function initializeYourself() on element level. More... | |
virtual int | giveNumberOfDomainEquations (int di, const UnknownNumberingScheme &num) |
Returns number of equations for given domain in active (current time step) time step. More... | |
virtual FieldPtr | giveField (FieldType key, TimeStep *) |
Returns the smart pointer to requested field, Null otherwise. More... | |
EngngModel * | giveMasterEngngModel () |
Returns the master engnmodel. More... | |
virtual double | giveEigenValue (int eigNum) |
Only relevant for eigen value analysis. Otherwise returns zero. More... | |
virtual void | setActiveVector (int i) |
Only relevant for eigen value analysis. Otherwise does noting. More... | |
int | updateSharedDofManagers (FloatArray &answer, const UnknownNumberingScheme &s, int ExchangeTag) |
Exchanges necessary remote DofManagers data. More... | |
int | exchangeRemoteElementData (int ExchangeTag) |
Exchanges necessary remote element data with remote partitions. More... | |
MPI_Comm | giveParallelComm () |
Returns the communication object of reciever. More... | |
int | packRemoteElementData (ProcessCommunicator &processComm) |
Packs data of local element to be received by their remote counterpart on remote partitions. More... | |
int | unpackRemoteElementData (ProcessCommunicator &processComm) |
Unpacks data for remote elements (which are mirrors of remote partition's local elements). More... | |
int | packDofManagers (ArrayWithNumbering *src, ProcessCommunicator &processComm) |
Packing function for vector values of DofManagers. More... | |
int | unpackDofManagers (ArrayWithNumbering *dest, ProcessCommunicator &processComm) |
Unpacking function for vector values of DofManagers . More... | |
ProblemCommunicator * | giveProblemCommunicator (EngngModelCommType t) |
void | initializeCommMaps (bool forceInit=false) |
virtual int | instanciateYourself (DataReader &dr, InputRecord *ir, const char *outFileName, const char *desc) |
Initializes whole problem according to its description stored in inputStream. More... | |
void | Instanciate_init () |
Initialization of the receiver state (opening the default output stream, empty domain creation, initialization of parallel context, etc) before Initialization form DataReader. More... | |
int | instanciateDomains (DataReader &dr) |
Instanciate problem domains by calling their instanciateYourself() service. More... | |
int | instanciateMetaSteps (DataReader &dr) |
Instanciate problem meta steps by calling their instanciateYourself() service. More... | |
virtual int | instanciateDefaultMetaStep (InputRecord *ir) |
Instanciate default metastep, if nmsteps is zero. More... | |
void | initMetaStepAttributes (MetaStep *mStep) |
Update e-model attributes attributes according to step metaStep attributes. More... | |
MetaStep * | giveCurrentMetaStep () |
Returns current meta step. More... | |
virtual TimeStep * | giveCurrentStep (bool force=false) |
Returns current time step. More... | |
virtual TimeStep * | givePreviousStep (bool force=false) |
Returns previous time step. More... | |
TimeStep * | generateNextStep () |
Generate new time step (and associate metastep). More... | |
virtual void | preInitializeNextStep () |
Does a pre-initialization of the next time step (implement if necessarry) More... | |
virtual int | giveNumberOfFirstStep (bool force=false) |
Returns number of first time step used by receiver. More... | |
int | giveNumberOfMetaSteps () |
Return number of meta steps. More... | |
MetaStep * | giveMetaStep (int i) |
Returns the i-th meta step. More... | |
int | giveNumberOfSteps (bool force=false) |
Returns total number of steps. More... | |
virtual double | giveEndOfTimeOfInterest () |
Returns end of time interest (time corresponding to end of time integration). More... | |
int | giveNumberOfTimeStepWhenIcApply () |
Returns the time step number, when initial conditions should apply. More... | |
ExportModuleManager * | giveExportModuleManager () |
Returns receiver's export module manager. More... | |
EngngModelTimer * | giveTimer () |
Returns reference to receiver timer (EngngModelTimer). More... | |
virtual int | giveNewEquationNumber (int domain, DofIDItem) |
Increases number of equations of receiver's domain and returns newly created equation number. More... | |
virtual int | giveNewPrescribedEquationNumber (int domain, DofIDItem) |
Increases number of prescribed equations of receiver's domain and returns newly created equation number. More... | |
std::string | giveContextFileName (int tStepNumber, int stepVersion) const |
Returns the filename for the context file for the given step and version. More... | |
std::string | giveDomainFileName (int domainNum, int domainSerNum) const |
Returns the filename for the given domain (used by adaptivity and restore) More... | |
virtual void | initStepIncrements () |
Initializes solution of new time step. More... | |
virtual int | forceEquationNumbering (int i) |
Forces equation renumbering on given domain. More... | |
virtual int | forceEquationNumbering () |
Forces equation renumbering on all domains associated to engng model. More... | |
virtual int | requiresUnknownsDictionaryUpdate () |
Indicates if EngngModel requires Dofs dictionaries to be updated. More... | |
virtual bool | requiresEquationRenumbering (TimeStep *tStep) |
Returns true if equation renumbering is required for given solution step. More... | |
virtual void | updateDofUnknownsDictionary (DofManager *, TimeStep *) |
Updates necessary values in Dofs unknown dictionaries. More... | |
virtual ParallelContext * | giveParallelContext (int n) |
Returns the parallel context corresponding to given domain (n) and unknown type Default implementation returns i-th context from parallelContextList. More... | |
virtual void | initParallelContexts () |
Creates parallel contexts. More... | |
virtual void | assemble (SparseMtrx &answer, TimeStep *tStep, const MatrixAssembler &ma, const UnknownNumberingScheme &r_s, const UnknownNumberingScheme &c_s, Domain *domain) |
Assembles characteristic matrix of required type into given sparse matrix. More... | |
void | assembleVector (FloatArray &answer, TimeStep *tStep, const VectorAssembler &va, ValueModeType mode, const UnknownNumberingScheme &s, Domain *domain, FloatArray *eNorms=NULL) |
Assembles characteristic vector of required type from dofManagers, element, and active boundary conditions, into given vector. More... | |
void | assembleVectorFromDofManagers (FloatArray &answer, TimeStep *tStep, const VectorAssembler &va, ValueModeType mode, const UnknownNumberingScheme &s, Domain *domain, FloatArray *eNorms=NULL) |
Assembles characteristic vector of required type from dofManagers into given vector. More... | |
void | assembleVectorFromElements (FloatArray &answer, TimeStep *tStep, const VectorAssembler &va, ValueModeType mode, const UnknownNumberingScheme &s, Domain *domain, FloatArray *eNorms=NULL) |
Assembles characteristic vector of required type from elements into given vector. More... | |
void | assembleVectorFromBC (FloatArray &answer, TimeStep *tStep, const VectorAssembler &va, ValueModeType mode, const UnknownNumberingScheme &s, Domain *domain, FloatArray *eNorms=NULL) |
Assembles characteristic vector of required type from boundary conditions. More... | |
void | assembleExtrapolatedForces (FloatArray &answer, TimeStep *tStep, CharType type, Domain *domain) |
Assembles the extrapolated internal forces vector, useful for obtaining a good initial guess in nonlinear analysis with Dirichlet boundary conditions. More... | |
void | assemblePrescribedExtrapolatedForces (FloatArray &answer, TimeStep *tStep, CharType type, Domain *domain) |
void | assembleVectorFromContacts (FloatArray &answer, TimeStep *tStep, CharType type, ValueModeType mode, const UnknownNumberingScheme &s, Domain *domain, FloatArray *eNorms=NULL) |
virtual int | checkProblemConsistency () |
Allows programmer to test problem its internal data, before computation begins. More... | |
virtual void | init () |
Initializes the receiver state. More... | |
virtual void | postInitialize () |
Performs post-initialization for all the problem contents (which is called after initializeFrom). More... | |
virtual void | printOutputAt (FILE *file, TimeStep *tStep, const IntArray &nodeSets, const IntArray &elementSets) |
void | outputNodes (FILE *file, Domain &domain, TimeStep *tStep, int setNum) |
Outputs all nodes in the given set. More... | |
void | outputElements (FILE *file, Domain &domain, TimeStep *tStep, int setNum) |
Outputs all elements in the given set. More... | |
void | printYourself () |
Prints state of receiver. Useful for debugging. More... | |
virtual void | printDofOutputAt (FILE *stream, Dof *iDof, TimeStep *tStep) |
DOF printing routine. More... | |
bool | isParallel () const |
Returns true if receiver in parallel mode. More... | |
int | giveRank () const |
Returns domain rank in a group of collaborating processes (0..groupSize-1) More... | |
int | giveNumberOfProcesses () const |
Returns the number of collaborating processes. More... | |
EngngModelContext * | giveContext () |
Context requesting service. More... | |
virtual int | giveNumberOfSlaveProblems () |
Returns number of slave problems. More... | |
virtual EngngModel * | giveSlaveProblem (int i) |
Returns i-th slave problem. More... | |
virtual bool | giveEquationScalingFlag () |
Returns the Equation scaling flag, which is used to indicate that governing equation(s) are scaled, or non-dimensionalized. More... | |
virtual double | giveVariableScale (VarScaleType varId) |
Returns the scale factor for given variable type. More... | |
virtual void | balanceLoad (TimeStep *tStep) |
Recovers the load balance between processors, if needed. More... | |
void | initParallel () |
Request domain rank and problem size. More... | |
EngngModel * | giveEngngModel () |
Returns reference to itself -> required by communicator.h. More... | |
virtual bool | isElementActivated (int elemNum) |
virtual bool | isElementActivated (Element *e) |
virtual void | drawYourself (oofegGraphicContext &gc) |
virtual void | drawElements (oofegGraphicContext &gc) |
virtual void | drawNodes (oofegGraphicContext &gc) |
std::string | errorInfo (const char *func) const |
Returns string for prepending output (used by error reporting macros). More... | |
Protected Member Functions | |
void | assembleInitialLoadVector (FloatArray &loadVector, FloatArray &loadVectorOfPrescribed, AdaptiveNonLinearStatic *sourceProblem, int domainIndx, TimeStep *tStep) |
Protected Member Functions inherited from oofem::NonLinearStatic | |
virtual void | assemble (SparseMtrx &answer, TimeStep *tStep, const MatrixAssembler &ma, const UnknownNumberingScheme &, Domain *domain) |
Assembles characteristic matrix of required type into given sparse matrix. More... | |
void | proceedStep (int di, TimeStep *tStep) |
virtual void | updateLoadVectors (TimeStep *tStep) |
virtual void | computeExternalLoadReactionContribution (FloatArray &reactions, TimeStep *tStep, int di) |
Computes the contribution external loading to reaction forces in given domain. More... | |
void | assembleIncrementalReferenceLoadVectors (FloatArray &_incrementalLoadVector, FloatArray &_incrementalLoadVectorOfPrescribed, SparseNonLinearSystemNM::referenceLoadInputModeType _refMode, Domain *sourceDomain, TimeStep *tStep) |
virtual void | packMigratingData (TimeStep *tStep) |
Packs receiver data when rebalancing load. More... | |
virtual void | unpackMigratingData (TimeStep *tStep) |
Unpacks receiver data when rebalancing load. More... | |
Protected Member Functions inherited from oofem::StructuralEngngModel | |
void | printReactionForces (TimeStep *tStep, int id, FILE *out) |
Computes and prints reaction forces, computed from nodal internal forces. More... | |
virtual void | giveInternalForces (FloatArray &answer, bool normFlag, int di, TimeStep *tStep) |
Evaluates the nodal representation of internal forces by assembling contributions from individual elements. More... | |
void | updateInternalState (TimeStep *tStep) |
Updates nodal values (calls also this->updateDofUnknownsDictionary for updating dofs unknowns dictionaries if model supports changes of static system). More... | |
Protected Attributes | |
FloatArray | d2_totalDisplacement |
FloatArray | d2_incrementOfDisplacement |
MeshPackageType | meshPackage |
int | equilibrateMappedConfigurationFlag |
Flag indication whether to restore equilibrium after adaptive remapping. More... | |
bool | preMappingLoadBalancingFlag |
Flag to trigger load balancing before adaptive remapping. More... | |
FloatArray | timeStepLoadLevels |
Array storing the load levels reached in corresponding time steps. More... | |
Protected Attributes inherited from oofem::NonLinearStatic | |
double | prevStepLength |
double | currentStepLength |
FloatArray | totalDisplacement |
FloatArray | incrementOfDisplacement |
FloatArray | internalForces |
FloatArray | initialLoadVector |
A load vector already applied, which does not scales. More... | |
FloatArray | incrementalLoadVector |
Incremental load vector which is applied in loading step (as a whole for direct control or proportionally for indirect control). More... | |
FloatArray | initialLoadVectorOfPrescribed |
A load vector which does not scale for prescribed DOFs. More... | |
FloatArray | incrementalLoadVectorOfPrescribed |
Incremental Load Vector for prescribed DOFs. More... | |
double | loadLevel |
double | cumulatedLoadLevel |
bool | mstepCumulateLoadLevelFlag |
int | currentIterations |
NonLinearStatic_stiffnessMode | stiffMode |
int | loadInitFlag |
int | nonlocalStiffnessFlag |
bool | updateElasticStiffnessFlag |
NM_Status | numMetStatus |
SparseNonLinearSystemNM * | nMethod |
Numerical method used to solve the problem. More... | |
NonLinearStatic_controlType | controlMode |
Characterizes the type of control used. More... | |
double | deltaT |
Intrinsic time increment. More... | |
SparseNonLinearSystemNM::referenceLoadInputModeType | refLoadInputMode |
The following parameter allows to specify how the reference load vector is obtained from given totalLoadVector and initialLoadVector. More... | |
InitialGuess | initialGuessType |
The initial guess type to use before starting the nonlinear solver. More... | |
Protected Attributes inherited from oofem::LinearStatic | |
std::unique_ptr< SparseMtrx > | stiffnessMatrix |
FloatArray | loadVector |
FloatArray | displacementVector |
LinSystSolverType | solverType |
SparseMtrxType | sparseMtrxType |
std::unique_ptr< SparseLinearSystemNM > | nMethod |
Numerical method used to solve the problem. More... | |
int | initFlag |
EModelDefaultEquationNumbering * | equationNumbering |
Protected Attributes inherited from oofem::StructuralEngngModel | |
StateCounterType | internalVarUpdateStamp |
Contains last time stamp of internal variable update. More... | |
FloatArray | internalForcesEBENorm |
Norm of nodal internal forces evaluated on element by element basis (squared) More... | |
Protected Attributes inherited from oofem::EngngModel | |
int | ndomains |
Number of receiver domains. More... | |
std::vector< std::unique_ptr< Domain > > | domainList |
List of problem domains. More... | |
int | numberOfSteps |
Total number of time steps. More... | |
int | numberOfEquations |
Total number of equation in current time step. More... | |
int | numberOfPrescribedEquations |
Total number or prescribed equations in current time step. More... | |
IntArray | domainNeqs |
Number of equations per domain. More... | |
IntArray | domainPrescribedNeqs |
Number of prescribed equations per domain. More... | |
bool | renumberFlag |
Renumbering flag (renumbers equations after each step, necessary if Dirichlet BCs change). More... | |
bool | profileOpt |
Profile optimized numbering flag (using Sloan's algorithm). More... | |
int | equationNumberingCompleted |
Equation numbering completed flag. More... | |
int | nMetaSteps |
Number of meta steps. More... | |
std::vector< MetaStep > | metaStepList |
List of problem metasteps. More... | |
std::unique_ptr< TimeStep > | stepWhenIcApply |
Solution step when IC (initial conditions) apply. More... | |
std::unique_ptr< TimeStep > | currentStep |
Current time step. More... | |
std::unique_ptr< TimeStep > | previousStep |
Previous time step. More... | |
int | number |
Receivers id. More... | |
std::string | dataOutputFileName |
Path to output stream. More... | |
std::string | coreOutputFileName |
String with core output file name. More... | |
FILE * | outputStream |
Output stream. More... | |
std::string | referenceFileName |
String with reference file name. More... | |
ContextOutputMode | contextOutputMode |
Domain context output mode. More... | |
int | contextOutputStep |
ExportModuleManager * | exportModuleManager |
Export module manager. More... | |
InitModuleManager * | initModuleManager |
Initialization module manager. More... | |
problemMode | pMode |
Domain mode. More... | |
problemScale | pScale |
Multiscale mode. More... | |
time_t | startTime |
Solution start time. More... | |
EngngModel * | master |
Master e-model; if defined receiver is in maintained (slave) mode. More... | |
EngngModelContext * | context |
Context. More... | |
EngngModelTimer | timer |
E-model timer. More... | |
int | parallelFlag |
Flag indicating that the receiver runs in parallel. More... | |
enum fMode | nonLinFormulation |
Type of non linear formulation (total or updated formulation). More... | |
ErrorEstimator * | defaultErrEstimator |
Error estimator. Useful for adaptivity, or simply printing errors output. More... | |
int | rank |
Domain rank in a group of collaborating processes (0..groupSize-1). More... | |
int | numProcs |
Total number of collaborating processes. More... | |
int | nonlocalExt |
Flag indicating if nonlocal extension active, which will cause data to be sent between shared elements before computing the internal forces. More... | |
char | processor_name [PROCESSOR_NAME_LENGTH] |
Processor name. More... | |
MPI_Comm | comm |
Communication object for this engineering model. More... | |
CommunicatorBuff * | commBuff |
Common Communicator buffer. More... | |
ProblemCommunicator * | communicator |
Communicator. More... | |
ProblemCommunicator * | nonlocCommunicator |
NonLocal Communicator. Necessary when nonlocal constitutive models are used. More... | |
std::vector< ParallelContext > | parallelContextList |
List where parallel contexts are stored. More... | |
bool | suppressOutput |
Flag for suppressing output to file. More... | |
std::string | simulationDescription |
LoadBalancer * | lb |
Load Balancer. More... | |
LoadBalancerMonitor * | lbm |
bool | loadBalancingFlag |
If set to true, load balancing is active. More... | |
bool | force_load_rebalance_in_first_step |
Debug flag forcing load balancing after first step. More... | |
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 } |
Means to choose methods for finding a good initial guess. More... | |
Protected Types inherited from oofem::EngngModel | |
enum | { InternalForcesExchangeTag, MassExchangeTag, LoadExchangeTag, ReactionExchangeTag, RemoteElementExchangeTag } |
Message tags. More... | |
This class implements Adaptive Non-LinearStatic Engineering problem.
Solution of this problem is performed as a series of increments (loading or displacement). At start of Each increment we assemble new tangent stiffness, and iteratively trying to fulfill balance of external and real internal forces at end of load step (see numerical method ).
The error is estimated at the end of each load increment, and based on reached error, the computation continues, or remeshing happens. The solution is then mapped to new mesh. The solution may be taken from current state or if error is too high then is taken from previous equilibrium state.
Tasks:
Definition at line 74 of file adaptnlinearstatic.h.
oofem::AdaptiveNonLinearStatic::AdaptiveNonLinearStatic | ( | int | i, |
EngngModel * | _master = NULL |
||
) |
Definition at line 72 of file adaptnlinearstatic.C.
References equilibrateMappedConfigurationFlag, meshPackage, and preMappingLoadBalancingFlag.
|
virtual |
Definition at line 84 of file adaptnlinearstatic.C.
|
virtual |
Remaps the solution state to newly given domain.
This includes mapping of source solution, internal variable mapping procedures and optionally restoring new global equilibrium. Given domain becomes new domain of receiver.
Definition at line 445 of file adaptnlinearstatic.C.
References oofem::NonLinearStatic::assemble(), oofem::NonLinearStatic::assembleIncrementalReferenceLoadVectors(), assembleInitialLoadVector(), oofem::IntArray::at(), oofem::classFactory, oofem::ClassFactory::createSparseMtrx(), oofem::NonLinearStatic::currentIterations, d2_incrementOfDisplacement, d2_totalDisplacement, oofem::EngngModel::domainList, oofem::EngngModel::domainNeqs, oofem::EngngModel::domainPrescribedNeqs, oofem::Element_remote, oofem::EngngModel::equationNumberingCompleted, equilibrateMappedConfigurationFlag, oofem::EngngModel::exchangeRemoteElementData(), oofem::EngngModel::forceEquationNumbering(), oofem::Timer::getUtime(), oofem::EngngModel::giveCurrentMetaStep(), oofem::EngngModel::giveCurrentStep(), oofem::SparseNonLinearSystemNM::giveCurrentStepLength(), oofem::EngngModel::giveDomain(), oofem::EngngModel::giveNumberOfDomainEquations(), oofem::EngngModel::giveOutputStream(), oofem::NonLinearStatic::incrementalLoadVector, oofem::NonLinearStatic::incrementalLoadVectorOfPrescribed, oofem::NonLinearStatic::incrementOfDisplacement, oofem::TimeStep::incrementVersion(), oofem::LinearStatic::initFlag, oofem::EngngModel::initializeCommMaps(), oofem::NonLinearStatic::initialLoadVector, oofem::NonLinearStatic::initialLoadVectorOfPrescribed, oofem::EngngModel::initStepIncrements(), oofem::NonLinearStatic::internalForces, oofem::StructuralEngngModel::internalForcesEBENorm, oofem::FloatArray::isNotEmpty(), oofem::EngngModel::isParallel(), oofem::NonLinearStatic::loadLevel, oofem::LinearStatic::loadVector, oofem::EIPrimaryUnknownMapper::mapAndUpdate(), oofem::EngngModel::ndomains, oofem::nls_secantStiffness, oofem::NonLinearStatic::nMethod, oofem::NonLinearStatic::nonlocalStiffnessFlag, oofem::EngngModel::numberOfEquations, oofem::EngngModel::numberOfPrescribedEquations, oofem::NonLinearStatic::numMetStatus, OOFEM_ERROR, OOFEM_LOG_DEBUG, OOFEM_LOG_INFO, OOFEM_LOG_RELEVANT, oofem::EngngModel::parallelContextList, oofem::NonLinearStatic::printOutputAt(), oofem::NonLinearStatic::refLoadInputMode, oofem::NumericalMethod::reinitialize(), oofem::EngngModel::RemoteElementExchangeTag, oofem::IntArray::resize(), oofem::FloatArray::resize(), oofem::SparseNonLinearSystemNM::setStepLength(), oofem::SparseNonLinearSystemNM::solve(), oofem::LinearStatic::sparseMtrxType, oofem::Timer::startTimer(), oofem::NonLinearStatic::stiffMode, oofem::LinearStatic::stiffnessMatrix, oofem::Timer::stopTimer(), oofem::NonLinearStatic::terminate(), oofem::EngngModel::timer, oofem::NonLinearStatic::totalDisplacement, oofem::NonLinearStatic::updateAttributes(), updateDomainLinks(), updateYourself(), and oofem::FloatArray::zero().
Referenced by initializeAdaptive(), and solveYourselfAt().
|
protected |
Definition at line 773 of file adaptnlinearstatic.C.
References _IFT_AdaptiveNonLinearStatic_controlmode, _IFT_AdaptiveNonLinearStatic_ddm, _IFT_AdaptiveNonLinearStatic_refloadmode, _IFT_NonLinearStatic_donotfixload, oofem::FloatArray::add(), oofem::NonLinearStatic::assembleIncrementalReferenceLoadVectors(), oofem::NonLinearStatic::deltaT, oofem::MetaStep::giveAttributesRecord(), oofem::EngngModel::giveDomain(), oofem::MetaStep::giveFirstStepNumber(), oofem::MetaStep::giveLastStepNumber(), oofem::EngngModel::giveMetaStep(), oofem::TimeStep::giveMetaStepNumber(), oofem::TimeStep::giveNumber(), oofem::EngngModel::giveNumberOfDomainEquations(), giveTimeStepLoadLevel(), oofem::InputRecord::hasField(), IR_GIVE_OPTIONAL_FIELD, oofem::nls_directControl, oofem::nls_indirectControl, OOFEM_ERROR, oofem::FloatArray::resize(), oofem::FloatArray::times(), and oofem::FloatArray::zero().
Referenced by adaptiveRemap(), giveInputRecordName(), and initializeAdaptiveFrom().
|
inlinevirtual |
Returns class name of the receiver.
Reimplemented from oofem::NonLinearStatic.
Definition at line 136 of file adaptnlinearstatic.h.
|
inlinevirtual |
Reimplemented from oofem::LinearStatic.
Definition at line 137 of file adaptnlinearstatic.h.
References _IFT_AdaptiveNonLinearStatic_Name, assembleInitialLoadVector(), giveLoadBalancer(), giveLoadBalancerMonitor(), and oofem::LinearStatic::loadVector.
|
virtual |
Returns reference to receiver's load balancer.
Reimplemented from oofem::NonLinearStatic.
Definition at line 966 of file adaptnlinearstatic.C.
References oofem::classFactory, oofem::ClassFactory::createLoadBalancer(), oofem::EngngModel::giveDomain(), oofem::EngngModel::lb, oofem::EngngModel::loadBalancingFlag, and preMappingLoadBalancingFlag.
Referenced by giveInputRecordName().
|
virtual |
Returns reference to receiver's load balancer monitor.
Reimplemented from oofem::NonLinearStatic.
Definition at line 980 of file adaptnlinearstatic.C.
References oofem::classFactory, oofem::ClassFactory::createLoadBalancerMonitor(), oofem::EngngModel::lbm, oofem::EngngModel::loadBalancingFlag, and preMappingLoadBalancingFlag.
Referenced by giveInputRecordName().
|
virtual |
Returns the load level corresponding to given solution step number.
Definition at line 952 of file adaptnlinearstatic.C.
References oofem::FloatArray::at(), oofem::EngngModel::giveNumberOfFirstStep(), oofem::EngngModel::giveNumberOfSteps(), OOFEM_ERROR, and timeStepLoadLevels.
Referenced by assembleInitialLoadVector().
|
virtual |
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.
Reimplemented from oofem::NonLinearStatic.
Definition at line 196 of file adaptnlinearstatic.C.
References oofem::Dof::__giveEquationNumber(), oofem::FloatArray::at(), d2_incrementOfDisplacement, d2_totalDisplacement, oofem::EngngModel::giveCurrentStep(), oofem::Domain::giveNumber(), oofem::NonLinearStatic::giveUnknownComponent(), oofem::FloatArray::isNotEmpty(), and OOFEM_ERROR.
|
virtual |
Initializes the newly generated discretization state according to previous solution.
This process should typically include restoring old solution, instanciating newly generated domain(s) and by mapping procedure.
Reimplemented from oofem::EngngModel.
Definition at line 419 of file adaptnlinearstatic.C.
References adaptiveRemap(), CM_State, oofem::EngngModel::giveContextFileName(), oofem::EngngModel::giveDomain(), oofem::EngngModel::giveDomainFileName(), oofem::Domain::giveSerialNumber(), oofem::EngngModel::initStepIncrements(), oofem::Domain::instanciateYourself(), OOFEM_ERROR, OOFEM_LOG_INFO, oofem::ContextIOERR::print(), and restoreContext().
|
virtual |
Initializes the receiver state according to state of given source problem.
This process should typically include mapping of source solution, internal variable mapping procedures and optionally restoring new global equilibrium.
Definition at line 238 of file adaptnlinearstatic.C.
References oofem::NonLinearStatic::assemble(), oofem::NonLinearStatic::assembleIncrementalReferenceLoadVectors(), assembleInitialLoadVector(), oofem::classFactory, oofem::ClassFactory::createSparseMtrx(), oofem::NonLinearStatic::currentIterations, oofem::EngngModel::currentStep, equilibrateMappedConfigurationFlag, oofem::Timer::getUtime(), oofem::EngngModel::giveCurrentMetaStep(), oofem::EngngModel::giveCurrentStep(), oofem::SparseNonLinearSystemNM::giveCurrentStepLength(), oofem::EngngModel::giveDomain(), oofem::Domain::giveElements(), oofem::EngngModel::giveNumberOfDomainEquations(), oofem::EngngModel::giveOutputStream(), oofem::EngngModel::givePreviousStep(), oofem::NonLinearStatic::incrementalLoadVector, oofem::NonLinearStatic::incrementalLoadVectorOfPrescribed, oofem::NonLinearStatic::incrementOfDisplacement, oofem::LinearStatic::initFlag, oofem::NonLinearStatic::initialLoadVector, oofem::NonLinearStatic::initialLoadVectorOfPrescribed, oofem::NonLinearStatic::internalForces, oofem::StructuralEngngModel::internalForcesEBENorm, oofem::FloatArray::isNotEmpty(), oofem::NonLinearStatic::loadLevel, oofem::LinearStatic::loadVector, oofem::EIPrimaryUnknownMapper::mapAndUpdate(), oofem::nls_secantStiffness, oofem::NonLinearStatic::nMethod, oofem::NonLinearStatic::nonlocalStiffnessFlag, oofem::NonLinearStatic::numMetStatus, OOFEM_ERROR, OOFEM_LOG_INFO, OOFEM_LOG_RELEVANT, oofem::EngngModel::previousStep, oofem::NonLinearStatic::printOutputAt(), oofem::NonLinearStatic::refLoadInputMode, oofem::FloatArray::resize(), oofem::SparseNonLinearSystemNM::setStepLength(), oofem::SparseNonLinearSystemNM::solve(), oofem::LinearStatic::sparseMtrxType, oofem::Timer::startTimer(), oofem::NonLinearStatic::stiffMode, oofem::LinearStatic::stiffnessMatrix, oofem::Timer::stopTimer(), oofem::NonLinearStatic::terminate(), oofem::EngngModel::timer, oofem::NonLinearStatic::totalDisplacement, oofem::NonLinearStatic::updateAttributes(), oofem::NonLinearStatic::updateLoadVectors(), updateYourself(), and oofem::FloatArray::zero().
|
virtual |
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::NonLinearStatic.
Definition at line 89 of file adaptnlinearstatic.C.
References _IFT_AdaptiveNonLinearStatic_equilmc, _IFT_AdaptiveNonLinearStatic_meshpackage, _IFT_AdaptiveNonLinearStatic_preMappingLoadBalancingFlag, oofem::EngngModel::defaultErrEstimator, equilibrateMappedConfigurationFlag, oofem::NonLinearStatic::initializeFrom(), IR_GIVE_OPTIONAL_FIELD, meshPackage, OOFEM_ERROR, and preMappingLoadBalancingFlag.
|
virtual |
Restores the state of model from output stream.
Restores not only the receiver state, but also same function is invoked for all DofManagers and Elements in associated domain. Note that by restoring element context also contexts of all associated integration points (and material statuses) are restored. Each context is associated with unique time step. Only one context per time step is allowed. Restore context function will restore such context, which is related (through its step number) to time step number and version given in obj parameter. Restoring context will change current time step in order to correspond to newly restored context.
stream | Context file. |
mode | Determines amount of info in stream. |
ContextIOERR | exception if error encountered. |
Reimplemented from oofem::NonLinearStatic.
Definition at line 748 of file adaptnlinearstatic.C.
References oofem::CIO_OK, oofem::NonLinearStatic::restoreContext(), oofem::FloatArray::restoreYourself(), THROW_CIOERR, and timeStepLoadLevels.
Referenced by initializeAdaptive().
|
virtual |
Stores the state of model to output stream.
Stores not only the receiver state, but also same function is invoked for all DofManagers and Elements in associated domain. Note that by storing element context also contexts of all associated integration points (and material statuses) are stored.
stream | Context stream. |
mode | Determines amount of info in stream. |
ContextIOERR | If error encountered. |
Reimplemented from oofem::NonLinearStatic.
Definition at line 732 of file adaptnlinearstatic.C.
References oofem::CIO_OK, oofem::NonLinearStatic::saveContext(), oofem::FloatArray::storeYourself(), THROW_CIOERR, and timeStepLoadLevels.
|
virtual |
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::NonLinearStatic.
Definition at line 115 of file adaptnlinearstatic.C.
References adaptiveRemap(), oofem::EngngModel::balanceLoad(), oofem::classFactory, oofem::COM_Always, oofem::MesherInterface::createMesh(), oofem::ClassFactory::createMesherInterface(), oofem::EngngModel::defaultErrEstimator, oofem::equilibratedEM, oofem::ErrorEstimator::estimateError(), oofem::RemeshingCriteria::estimateMeshDensities(), oofem::EngngModel::giveCurrentStep(), oofem::EngngModel::giveDomain(), oofem::ErrorEstimator::giveRemeshingCrit(), oofem::RemeshingCriteria::giveRemeshingStrategy(), oofem::LinearStatic::initFlag, meshPackage, oofem::MesherInterface::MI_NEEDS_EXTERNAL_ACTION, oofem::MesherInterface::MI_OK, oofem::NoRemeshing_RS, OOFEM_ERROR, preMappingLoadBalancingFlag, oofem::NonLinearStatic::proceedStep(), oofem::RemeshingFromCurrentState_RS, oofem::RemeshingFromPreviousState_RS, oofem::EngngModel::setContextOutputMode(), oofem::NonLinearStatic::terminate(), oofem::EngngModel::terminateAnalysis(), and updateYourself().
|
virtual |
Updates domain links after the domains of receiver have changed.
Used mainly after restoring context - the domains may change and this service is then used to update domain variables in all components belonging to receiver like error estimators, solvers, etc, having domains as attributes.
Reimplemented from oofem::NonLinearStatic.
Definition at line 765 of file adaptnlinearstatic.C.
References oofem::EngngModel::defaultErrEstimator, oofem::EngngModel::giveDomain(), oofem::ErrorEstimator::setDomain(), and oofem::NonLinearStatic::updateDomainLinks().
Referenced by adaptiveRemap().
|
virtual |
Updates internal state after finishing time step.
(for example total values may be updated according to previously solved increments). Then element values are also updated (together with related integration points and material statuses).
Reimplemented from oofem::StructuralEngngModel.
Definition at line 178 of file adaptnlinearstatic.C.
References oofem::FloatArray::at(), oofem::TimeStep::giveNumber(), oofem::EngngModel::giveNumberOfSteps(), oofem::FloatArray::isEmpty(), oofem::NonLinearStatic::loadLevel, oofem::FloatArray::resize(), timeStepLoadLevels, and oofem::StructuralEngngModel::updateYourself().
Referenced by adaptiveRemap(), initializeAdaptiveFrom(), and solveYourselfAt().
|
protected |
Definition at line 77 of file adaptnlinearstatic.h.
Referenced by adaptiveRemap(), and giveUnknownComponent().
|
protected |
Definition at line 77 of file adaptnlinearstatic.h.
Referenced by adaptiveRemap(), and giveUnknownComponent().
|
protected |
Flag indication whether to restore equilibrium after adaptive remapping.
Definition at line 80 of file adaptnlinearstatic.h.
Referenced by AdaptiveNonLinearStatic(), adaptiveRemap(), initializeAdaptiveFrom(), and initializeFrom().
|
protected |
Definition at line 78 of file adaptnlinearstatic.h.
Referenced by AdaptiveNonLinearStatic(), initializeFrom(), and solveYourselfAt().
|
protected |
Flag to trigger load balancing before adaptive remapping.
The emodel loadBalancingFlag will trigger load balancing after mapping and optional consistency recovery.
Definition at line 86 of file adaptnlinearstatic.h.
Referenced by AdaptiveNonLinearStatic(), giveLoadBalancer(), giveLoadBalancerMonitor(), initializeFrom(), and solveYourselfAt().
|
protected |
Array storing the load levels reached in corresponding time steps.
It is necessary to keep track of this load level history, because after adaptive restart one has to assemble the initial and total load vectors on new mesh -> and for this it is necessary to know the history of loading. The size of this array is equal to numberOfSolutionSteps and should be stored/restored in every context file.
Definition at line 96 of file adaptnlinearstatic.h.
Referenced by giveTimeStepLoadLevel(), restoreContext(), saveContext(), and updateYourself().