OOFEM  2.4
OOFEM.org - Object Oriented Finite Element Solver
oofem::NonLinearDynamic Class Reference

This class implements nonlinear dynamic engineering problem. More...

#include <nlineardynamic.h>

+ Inheritance diagram for oofem::NonLinearDynamic:
+ Collaboration diagram for oofem::NonLinearDynamic:

Public Member Functions

 NonLinearDynamic (int i, EngngModel *_master=NULL)
 
virtual ~NonLinearDynamic ()
 
virtual void solveYourself ()
 Starts solution process. More...
 
virtual void solveYourselfAt (TimeStep *tStep)
 Solves problem for given 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 printDofOutputAt (FILE *stream, Dof *iDof, TimeStep *tStep)
 DOF printing routine. More...
 
virtual void updateYourself (TimeStep *tStep)
 Updates internal state after finishing 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 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 double giveUnknownComponent (ValueModeType type, TimeStep *tStep, Domain *d, Dof *dof)
 Returns requested unknown. More...
 
virtual IRResultType initializeFrom (InputRecord *ir)
 Initializes receiver according to object description in input reader. More...
 
virtual TimeStepgiveNextStep ()
 Returns next time step (next to current step) of receiver. More...
 
virtual NumericalMethodgiveNumericalMethod (MetaStep *mStep)
 Returns reference to receiver's numerical method. 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 * giveInputRecordName () const
 
virtual const char * giveClassName () const
 Returns class name of the receiver. 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)
 This method is responsible for computing unique dictionary id (ie hash value) from given valueModeType and time step. More...
 
void timesMtrx (FloatArray &answer, FloatArray &vec, CharType type, Domain *domain, TimeStep *tStep)
 
TimeDiscretizationType giveInitialTimeDiscretization ()
 
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...
 
virtual LoadBalancergiveLoadBalancer ()
 Returns reference to receiver's load balancer. More...
 
virtual LoadBalancerMonitorgiveLoadBalancerMonitor ()
 Returns reference to receiver's load balancer monitor. More...
 
- 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...
 
virtual void terminate (TimeStep *tStep)
 Terminates the solution of time step. 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
 
EngngModeloperator= (const EngngModel &)=delete
 
DomaingiveDomain (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 ErrorEstimatorgiveDomainErrorEstimator (int n)
 Service for accessing ErrorEstimator corresponding to particular domain. More...
 
virtual MaterialInterfacegiveMaterialInterface (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 int initializeAdaptive (int tStepNumber)
 Initializes the newly generated discretization state according to previous solution. 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...
 
EngngModelgiveMasterEngngModel ()
 Returns the master engnmodel. More...
 
virtual double giveLoadLevel ()
 Returns the current load level. 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...
 
virtual int giveCurrentNumberOfIterations ()
 Returns number of iterations that was required to reach equilibrium - used for adaptive step length in staggered problem. 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...
 
ProblemCommunicatorgiveProblemCommunicator (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...
 
MetaStepgiveCurrentMetaStep ()
 Returns current meta step. More...
 
virtual TimeStepgiveCurrentStep (bool force=false)
 Returns current time step. More...
 
virtual TimeStepgivePreviousStep (bool force=false)
 Returns previous time step. More...
 
TimeStepgenerateNextStep ()
 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 TimeStepgiveSolutionStepWhenIcApply (bool force=false)
 Returns the solution step when Initial Conditions (IC) apply. 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...
 
MetaStepgiveMetaStep (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...
 
ExportModuleManagergiveExportModuleManager ()
 Returns receiver's export module manager. More...
 
EngngModelTimergiveTimer ()
 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 ParallelContextgiveParallelContext (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...
 
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...
 
EngngModelContextgiveContext ()
 Context requesting service. More...
 
virtual int giveNumberOfSlaveProblems ()
 Returns number of slave problems. More...
 
virtual EngngModelgiveSlaveProblem (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...
 
EngngModelgiveEngngModel ()
 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 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)
 
void determineConstants (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 computeExternalLoadReactionContribution (FloatArray &reactions, TimeStep *tStep, int di)
 Computes the contribution external loading to reaction forces in given domain. 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

std::unique_ptr< SparseMtrxeffectiveStiffnessMatrix
 
std::unique_ptr< SparseMtrxmassMatrix
 
LinSystSolverType solverType
 
SparseMtrxType sparseMtrxType
 
int initFlag
 
TimeDiscretizationType initialTimeDiscretization
 
double gamma
 
double beta
 
double eta
 
double delta
 
double a0
 
double a1
 
double a2
 
double a3
 
double a4
 
double a5
 
double a6
 
double a7
 
FloatArray velocityVector
 
FloatArray accelerationVector
 
FloatArray previousLoadVector
 
FloatArray previousVelocityVector
 
FloatArray previousAccelerationVector
 
FloatArray help
 
FloatArray rhs
 
FloatArray rhs2
 
FloatArray previousInternalForces
 
FloatArray previousIncrementOfDisplacement
 
FloatArray previousTotalDisplacement
 
FloatArray totalDisplacement
 
FloatArray incrementOfDisplacement
 
FloatArray internalForces
 
FloatArray forcesVector
 
int currentIterations
 
int totIterations
 
int MANRMSteps
 
int commInitFlag
 
int nonlocalStiffnessFlag
 
std::unique_ptr< SparseNonLinearSystemNMnMethod
 Numerical method used to solve the problem. More...
 
double deltaT
 Intrinsic time increment. More...
 
- 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< MetaStepmetaStepList
 List of problem metasteps. More...
 
std::unique_ptr< TimeStepstepWhenIcApply
 Solution step when IC (initial conditions) apply. More...
 
std::unique_ptr< TimeStepcurrentStep
 Current time step. More...
 
std::unique_ptr< TimeSteppreviousStep
 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
 
ExportModuleManagerexportModuleManager
 Export module manager. More...
 
InitModuleManagerinitModuleManager
 Initialization module manager. More...
 
problemMode pMode
 Domain mode. More...
 
problemScale pScale
 Multiscale mode. More...
 
time_t startTime
 Solution start time. More...
 
EngngModelmaster
 Master e-model; if defined receiver is in maintained (slave) mode. More...
 
EngngModelContextcontext
 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...
 
ErrorEstimatordefaultErrEstimator
 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...
 
CommunicatorBuffcommBuff
 Common Communicator buffer. More...
 
ProblemCommunicatorcommunicator
 Communicator. More...
 
ProblemCommunicatornonlocCommunicator
 NonLocal Communicator. Necessary when nonlocal constitutive models are used. More...
 
std::vector< ParallelContextparallelContextList
 List where parallel contexts are stored. More...
 
bool suppressOutput
 Flag for suppressing output to file. More...
 
std::string simulationDescription
 
LoadBalancerlb
 Load Balancer. More...
 
LoadBalancerMonitorlbm
 
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...
 

Detailed Description

This class implements nonlinear dynamic 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 loading applied can bo of two types:

  • proportional incremental loading
  • non-proportional fixed loading, reflecting the previous history, but could not be scaled (like dead weight).

Tasks:

  • Creating Numerical method for solving nonlinear problem.
  • Assembling tangent stiffness matrix.
  • Interfacing Numerical method to Elements.
  • Managing time steps.

Solution proceedure described in: A SURVEY OF DIRECT TIME-INTEGRATION METHODS IN COMPUTATIONAL STRUCTURAL DYNAMICS - II. IMPLICIT METHODS K. Subbaraj and M. A. Dokainish Computers & Structures Vol. 32. No. 6. pp. 1387-1401, 1989

Author
Andreas Feymark
Alper Cesur

Definition at line 85 of file nlineardynamic.h.

Constructor & Destructor Documentation

oofem::NonLinearDynamic::NonLinearDynamic ( int  i,
EngngModel _master = NULL 
)
oofem::NonLinearDynamic::~NonLinearDynamic ( )
virtual

Definition at line 80 of file nlineardynamic.C.

Member Function Documentation

void oofem::NonLinearDynamic::assemble ( SparseMtrx answer,
TimeStep tStep,
const MatrixAssembler ma,
const UnknownNumberingScheme s,
Domain domain 
)
protectedvirtual

Assembles characteristic matrix of required type into given sparse matrix.

Parameters
answerAssembled matrix.
tStepTime step, when answer is assembled.
sDetermines the equation numbering scheme.
typeCharacteristic components of type type are requested from elements and assembled.
domainSource domain.
Todo:
This rotation matrix is not flexible enough.. it can only work with full size matrices and doesn't allow for flexibility in the matrixassembler.
Note
: Some active bcs still make changes even when they are not applied We should probably reconsider this approach, so that they e.g. just prescribe their lagrange mult. instead.

Reimplemented from oofem::EngngModel.

Definition at line 739 of file nlineardynamic.C.

References oofem::EngngModel::assemble(), oofem::Timer::getUtime(), oofem::Domain::giveElements(), nonlocalStiffnessFlag, OOFEM_LOG_DEBUG, oofem::SparseMtrx::printStatistics(), oofem::Timer::startTimer(), oofem::Timer::stopTimer(), and oofem::EngngModel::timer.

Referenced by giveInitialTimeDiscretization(), proceedStep(), and updateComponent().

int oofem::NonLinearDynamic::estimateMaxPackSize ( IntArray commMap,
DataStream buff,
int  packUnpackType 
)
virtual

Determines the space necessary for send/receive buffer.

It uses related communication map pattern to determine the maximum size needed.

Parameters
commMapCommunication map used to send/receive messages.
buffCommunication buffer.
packUnpackTypeDetermines the type of packed quantity, used by receiver to estimate the size of pack/unpack buffer accordingly.
Returns
Upper bound of space needed.
Todo:
Fix this old ProblemCommMode__NODE_CUT value

Reimplemented from oofem::EngngModel.

Definition at line 843 of file nlineardynamic.C.

References oofem::Element::estimatePackSize(), oofem::Domain::giveDofManager(), oofem::EngngModel::giveDomain(), oofem::Domain::giveElement(), oofem::DataStream::givePackSizeOfDouble(), and oofem::max().

Referenced by giveInitialTimeDiscretization().

virtual const char* oofem::NonLinearDynamic::giveClassName ( ) const
inlinevirtual

Returns class name of the receiver.

Implements oofem::EngngModel.

Definition at line 142 of file nlineardynamic.h.

virtual fMode oofem::NonLinearDynamic::giveFormulation ( )
inlinevirtual

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 143 of file nlineardynamic.h.

References oofem::EngngModel::nonLinFormulation.

virtual const char* oofem::NonLinearDynamic::giveInputRecordName ( ) const
inlinevirtual

Definition at line 141 of file nlineardynamic.h.

References _IFT_NonLinearDynamic_Name.

LoadBalancer * oofem::NonLinearDynamic::giveLoadBalancer ( )
virtual
LoadBalancerMonitor * oofem::NonLinearDynamic::giveLoadBalancerMonitor ( )
virtual

Returns reference to receiver's load balancer monitor.

Reimplemented from oofem::EngngModel.

Definition at line 891 of file nlineardynamic.C.

References oofem::classFactory, oofem::ClassFactory::createLoadBalancerMonitor(), oofem::EngngModel::lbm, and oofem::EngngModel::loadBalancingFlag.

Referenced by giveInitialTimeDiscretization().

TimeStep * oofem::NonLinearDynamic::giveNextStep ( )
virtual
NumericalMethod * oofem::NonLinearDynamic::giveNumericalMethod ( MetaStep mStep)
virtual

Returns reference to receiver's numerical method.

Reimplemented from oofem::EngngModel.

Definition at line 85 of file nlineardynamic.C.

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

Referenced by proceedStep(), and updateDomainLinks().

double oofem::NonLinearDynamic::giveUnknownComponent ( ValueModeType  ,
TimeStep ,
Domain ,
Dof  
)
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.

See also
Dof::giveUnknown

Reimplemented from oofem::EngngModel.

Definition at line 191 of file nlineardynamic.C.

References oofem::Dof::__giveEquationNumber(), accelerationVector, oofem::FloatArray::at(), oofem::EngngModel::giveCurrentStep(), incrementOfDisplacement, OOFEM_ERROR, totalDisplacement, and velocityVector.

virtual int oofem::NonLinearDynamic::giveUnknownDictHashIndx ( ValueModeType  mode,
TimeStep tStep 
)
inlinevirtual

This method is responsible for computing unique dictionary id (ie hash value) from given valueModeType and time step.

This function is used by particular dofs to access unknown identified by given parameters from its dictionary using computed index. Usually the hash algorithm should produce index that depend on time step relatively to actual one to avoid storage of complete history.

Reimplemented from oofem::EngngModel.

Definition at line 145 of file nlineardynamic.h.

References timesMtrx().

void oofem::NonLinearDynamic::initializeYourself ( TimeStep tStep)
virtual

Provides the opportunity to initialize state variables stored in element integration points according to initial conditions using function initializeYourself() on element level.

Should be called when current time step is time step when IC will apply (see EngngModel::giveNumberOfTimeStepWhenIcApply) somewhere from solveYourselfAt function). Implementation must be provided. Default implementation is empty.

Reimplemented from oofem::EngngModel.

Definition at line 305 of file nlineardynamic.C.

References accelerationVector, oofem::FloatArray::at(), deltaT, oofem::Domain::giveDofManagers(), oofem::EngngModel::giveDomain(), oofem::StructuralEngngModel::giveInternalForces(), oofem::EngngModel::giveNumberOfDomainEquations(), oofem::EngngModel::giveNumberOfTimeStepWhenIcApply(), incrementOfDisplacement, initFlag, internalForces, oofem::TimeStep::isTheFirstStep(), previousAccelerationVector, previousIncrementOfDisplacement, previousInternalForces, previousTotalDisplacement, previousVelocityVector, oofem::FloatArray::resize(), oofem::EngngModel::stepWhenIcApply, totalDisplacement, velocityVector, and oofem::FloatArray::zero().

void oofem::NonLinearDynamic::packMigratingData ( TimeStep tStep)
protectedvirtual

Packs receiver data when rebalancing load.

When rebalancing happens, the local numbering will be lost on majority of processors. Instead of identifying values of solution vectors that have to be send/received and then performing renumbering, all solution vectors are assumed to be stored in dof dictionaries before data migration. Then dofs will take care themselves for packing and unpacking. After data migration and local renumbering, the solution vectors will be restored from dof dictionary data back.

Reimplemented from oofem::EngngModel.

Definition at line 907 of file nlineardynamic.C.

References oofem::FloatArray::at(), oofem::Domain::giveDofManager(), oofem::EngngModel::giveDomain(), oofem::Domain::giveNumberOfDofManagers(), and totalDisplacement.

Referenced by giveInitialTimeDiscretization().

void oofem::NonLinearDynamic::printDofOutputAt ( FILE *  stream,
Dof iDof,
TimeStep tStep 
)
virtual

DOF printing routine.

Called by DofManagers to print Dof specific part. Dof class provides component printing routines, but emodel is responsible for what will be printed at DOF level.

Parameters
streamoutput stream
iDofdof to be processed
tStepsolution step

Reimplemented from oofem::EngngModel.

Definition at line 657 of file nlineardynamic.C.

References oofem::Dof::printMultipleOutputAt().

void oofem::NonLinearDynamic::printOutputAt ( FILE *  file,
TimeStep tStep 
)
virtual

Prints output of receiver to output domain stream, for given time step.

Corresponding function for element gauss points is invoked (gaussPoint::printOutputAt).

Reimplemented from oofem::StructuralEngngModel.

Definition at line 640 of file nlineardynamic.C.

References currentIterations, oofem::OutputManager::doDofManOutput(), oofem::OutputManager::doElementOutput(), oofem::EngngModel::giveDomain(), oofem::TimeStep::giveNumber(), oofem::Domain::giveOutputManager(), oofem::TimeStep::giveTargetTime(), nMethod, and oofem::StructuralEngngModel::printReactionForces().

contextIOResultType oofem::NonLinearDynamic::restoreContext ( DataStream stream,
ContextMode  mode 
)
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.

Parameters
streamContext file.
modeDetermines amount of info in stream.
Returns
contextIOResultType.
Exceptions
ContextIOERRexception if error encountered.

Reimplemented from oofem::EngngModel.

Definition at line 695 of file nlineardynamic.C.

References accelerationVector, oofem::CIO_OK, incrementOfDisplacement, oofem::EngngModel::restoreContext(), oofem::FloatArray::restoreYourself(), THROW_CIOERR, totalDisplacement, and velocityVector.

contextIOResultType oofem::NonLinearDynamic::saveContext ( DataStream stream,
ContextMode  mode 
)
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.

Parameters
streamContext stream.
modeDetermines amount of info in stream.
Returns
contextIOResultType.
Exceptions
ContextIOERRIf error encountered.

Reimplemented from oofem::EngngModel.

Definition at line 667 of file nlineardynamic.C.

References accelerationVector, oofem::CIO_OK, incrementOfDisplacement, oofem::EngngModel::saveContext(), oofem::FloatArray::storeYourself(), THROW_CIOERR, totalDisplacement, and velocityVector.

void oofem::NonLinearDynamic::showSparseMtrxStructure ( int  type,
oofegGraphicContext gc,
TimeStep tStep 
)
virtual

Shows the sparse structure of required matrix, type == 1 stiffness.

Reimplemented from oofem::StructuralEngngModel.

Definition at line 767 of file nlineardynamic.C.

References oofem::EngngModel::giveDomain(), and oofem::Domain::giveElements().

Referenced by giveInitialTimeDiscretization().

void oofem::NonLinearDynamic::solveYourself ( )
virtual

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

Reimplemented from oofem::EngngModel.

Definition at line 262 of file nlineardynamic.C.

References commInitFlag, oofem::EngngModel::giveNumberOfDomainEquations(), oofem::EngngModel::giveRank(), oofem::EngngModel::initializeCommMaps(), oofem::EngngModel::isParallel(), OOFEM_LOG_INFO, and oofem::EngngModel::solveYourself().

void oofem::NonLinearDynamic::solveYourselfAt ( TimeStep tStep)
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::EngngModel.

Definition at line 282 of file nlineardynamic.C.

References commInitFlag, oofem::TimeStep::giveNumber(), oofem::EngngModel::giveNumberOfDomainEquations(), oofem::EngngModel::giveRank(), oofem::TimeStep::giveTargetTime(), oofem::EngngModel::initializeCommMaps(), oofem::EngngModel::isParallel(), OOFEM_LOG_INFO, OOFEM_LOG_RELEVANT, and proceedStep().

void oofem::NonLinearDynamic::unpackMigratingData ( TimeStep tStep)
protectedvirtual

Unpacks receiver data when rebalancing load.

When rebalancing happens, the local numbering will be lost on majority of processors. Instead of identifying values of solution vectors that have to be send/received and then performing renumbering, all solution vectors are assumed to be stored in dof dictionaries before data migration. Then dofs will take care themselves for packing and unpacking. After data migration and local renumbering, the solution vectors will be restored from dof dictionary data back.

Reimplemented from oofem::EngngModel.

Definition at line 926 of file nlineardynamic.C.

References oofem::FloatArray::at(), oofem::DofManager_shared, oofem::Domain::giveDofManager(), oofem::EngngModel::giveDomain(), oofem::EngngModel::giveDomainErrorEstimator(), oofem::Domain::giveNumberOfDofManagers(), oofem::EngngModel::giveNumberOfDomainEquations(), incrementOfDisplacement, initFlag, oofem::EngngModel::initializeCommMaps(), nMethod, oofem::ErrorEstimator::reinitialize(), oofem::FloatArray::resize(), and totalDisplacement.

Referenced by giveInitialTimeDiscretization().

void oofem::NonLinearDynamic::updateAttributes ( MetaStep mStep)
virtual

Update receiver attributes according to step metaStep attributes.

Allows the certain parameters or attributes to be updated for particular metastep. The metastep provides the attributes record, from which the corresponding attributes can be read. The service takes a MetaStep parameter. It is recommended, to implement this service in such way, that multiple calls for steps belonging to same MetaStep does not change response. The default implementation updates the numerical method attributes.

Parameters
mStepMeta step.

Reimplemented from oofem::EngngModel.

Definition at line 102 of file nlineardynamic.C.

References _IFT_NonLinearDynamic_delta, _IFT_NonLinearDynamic_deltat, _IFT_NonLinearDynamic_eta, delta, deltaT, eta, oofem::MetaStep::giveAttributesRecord(), IR_GIVE_OPTIONAL_FIELD, OOFEM_ERROR, and oofem::EngngModel::updateAttributes().

void oofem::NonLinearDynamic::updateComponent ( TimeStep tStep,
NumericalCmpn  cmpn,
Domain d 
)
virtual

Updates components mapped to numerical method if necessary during solution process.

Some numerical methods may require updating mapped components during solution process (e.g., updating of tangent stiffness when using updated Newton-Raphson method).

Parameters
tStepTime when component is updated.
cmpnNumerical component to update.
dDomain.

Reimplemented from oofem::EngngModel.

Definition at line 560 of file nlineardynamic.C.

References a0, a1, oofem::FloatArray::add(), assemble(), oofem::FloatArray::beScaled(), currentIterations, delta, effectiveStiffnessMatrix, eta, forcesVector, oofem::Timer::getUtime(), oofem::EngngModel::giveDomain(), oofem::StructuralEngngModel::giveInternalForces(), oofem::Domain::giveNumber(), help, incrementOfDisplacement, internalForces, oofem::InternalRhs, MANRMSteps, massMatrix, oofem::NonLinearLhs, OOFEM_ERROR, OOFEM_LOG_DEBUG, previousInternalForces, rhs2, oofem::Timer::startTimer(), oofem::Timer::stopTimer(), oofem::FloatArray::subtract(), oofem::EngngModel::timer, timesMtrx(), and totIterations.

void oofem::NonLinearDynamic::updateDomainLinks ( )
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::EngngModel.

Definition at line 724 of file nlineardynamic.C.

References oofem::EngngModel::giveCurrentMetaStep(), oofem::EngngModel::giveDomain(), giveLoadBalancer(), giveNumericalMethod(), oofem::NumericalMethod::setDomain(), oofem::LoadBalancer::setDomain(), and oofem::EngngModel::updateDomainLinks().

void oofem::NonLinearDynamic::updateYourself ( TimeStep tStep)
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 547 of file nlineardynamic.C.

References accelerationVector, incrementOfDisplacement, internalForces, previousAccelerationVector, previousIncrementOfDisplacement, previousInternalForces, previousTotalDisplacement, previousVelocityVector, totalDisplacement, totIterations, oofem::StructuralEngngModel::updateYourself(), and velocityVector.

virtual int oofem::NonLinearDynamic::useNonlocalStiffnessOption ( )
inlinevirtual

Returns nonzero if nonlocal stiffness option activated.

Reimplemented from oofem::EngngModel.

Definition at line 144 of file nlineardynamic.h.

References nonlocalStiffnessFlag.

Member Data Documentation

double oofem::NonLinearDynamic::a0
protected

Definition at line 97 of file nlineardynamic.h.

Referenced by determineConstants(), proceedStep(), and updateComponent().

double oofem::NonLinearDynamic::a1
protected

Definition at line 97 of file nlineardynamic.h.

Referenced by determineConstants(), proceedStep(), and updateComponent().

double oofem::NonLinearDynamic::a2
protected

Definition at line 97 of file nlineardynamic.h.

Referenced by determineConstants(), and proceedStep().

double oofem::NonLinearDynamic::a3
protected

Definition at line 97 of file nlineardynamic.h.

Referenced by determineConstants(), and proceedStep().

double oofem::NonLinearDynamic::a4
protected

Definition at line 97 of file nlineardynamic.h.

Referenced by determineConstants(), and proceedStep().

double oofem::NonLinearDynamic::a5
protected

Definition at line 97 of file nlineardynamic.h.

Referenced by determineConstants(), and proceedStep().

double oofem::NonLinearDynamic::a6
protected

Definition at line 97 of file nlineardynamic.h.

Referenced by determineConstants(), and proceedStep().

double oofem::NonLinearDynamic::a7
protected

Definition at line 97 of file nlineardynamic.h.

FloatArray oofem::NonLinearDynamic::accelerationVector
protected
double oofem::NonLinearDynamic::beta
protected

Definition at line 95 of file nlineardynamic.h.

Referenced by determineConstants(), and initializeFrom().

int oofem::NonLinearDynamic::commInitFlag
protected

Definition at line 107 of file nlineardynamic.h.

Referenced by NonLinearDynamic(), solveYourself(), and solveYourselfAt().

int oofem::NonLinearDynamic::currentIterations
protected

Definition at line 106 of file nlineardynamic.h.

Referenced by printOutputAt(), proceedStep(), and updateComponent().

double oofem::NonLinearDynamic::delta
protected

Definition at line 96 of file nlineardynamic.h.

Referenced by proceedStep(), updateAttributes(), and updateComponent().

double oofem::NonLinearDynamic::deltaT
protected

Intrinsic time increment.

Definition at line 112 of file nlineardynamic.h.

Referenced by determineConstants(), giveNextStep(), initializeFrom(), initializeYourself(), and updateAttributes().

std :: unique_ptr< SparseMtrx > oofem::NonLinearDynamic::effectiveStiffnessMatrix
protected

Definition at line 88 of file nlineardynamic.h.

Referenced by proceedStep(), and updateComponent().

double oofem::NonLinearDynamic::eta
protected

Definition at line 96 of file nlineardynamic.h.

Referenced by proceedStep(), updateAttributes(), and updateComponent().

FloatArray oofem::NonLinearDynamic::forcesVector
protected

Definition at line 104 of file nlineardynamic.h.

Referenced by proceedStep(), and updateComponent().

double oofem::NonLinearDynamic::gamma
protected

Definition at line 95 of file nlineardynamic.h.

Referenced by determineConstants(), and initializeFrom().

FloatArray oofem::NonLinearDynamic::help
protected

Definition at line 101 of file nlineardynamic.h.

Referenced by proceedStep(), and updateComponent().

FloatArray oofem::NonLinearDynamic::incrementOfDisplacement
protected
int oofem::NonLinearDynamic::initFlag
protected
TimeDiscretizationType oofem::NonLinearDynamic::initialTimeDiscretization
protected
FloatArray oofem::NonLinearDynamic::internalForces
protected

Definition at line 104 of file nlineardynamic.h.

Referenced by initializeYourself(), proceedStep(), updateComponent(), and updateYourself().

int oofem::NonLinearDynamic::MANRMSteps
protected

Definition at line 106 of file nlineardynamic.h.

Referenced by initializeFrom(), and updateComponent().

std :: unique_ptr< SparseMtrx > oofem::NonLinearDynamic::massMatrix
protected

Definition at line 88 of file nlineardynamic.h.

Referenced by proceedStep(), and updateComponent().

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

Numerical method used to solve the problem.

Definition at line 110 of file nlineardynamic.h.

Referenced by giveNumericalMethod(), printOutputAt(), proceedStep(), and unpackMigratingData().

int oofem::NonLinearDynamic::nonlocalStiffnessFlag
protected
FloatArray oofem::NonLinearDynamic::previousAccelerationVector
protected

Definition at line 100 of file nlineardynamic.h.

Referenced by initializeYourself(), proceedStep(), and updateYourself().

FloatArray oofem::NonLinearDynamic::previousIncrementOfDisplacement
protected

Definition at line 102 of file nlineardynamic.h.

Referenced by initializeYourself(), proceedStep(), and updateYourself().

FloatArray oofem::NonLinearDynamic::previousInternalForces
protected

Definition at line 101 of file nlineardynamic.h.

Referenced by initializeYourself(), proceedStep(), updateComponent(), and updateYourself().

FloatArray oofem::NonLinearDynamic::previousLoadVector
protected

Definition at line 99 of file nlineardynamic.h.

FloatArray oofem::NonLinearDynamic::previousTotalDisplacement
protected

Definition at line 103 of file nlineardynamic.h.

Referenced by initializeYourself(), proceedStep(), and updateYourself().

FloatArray oofem::NonLinearDynamic::previousVelocityVector
protected

Definition at line 100 of file nlineardynamic.h.

Referenced by initializeYourself(), proceedStep(), and updateYourself().

FloatArray oofem::NonLinearDynamic::rhs
protected

Definition at line 101 of file nlineardynamic.h.

Referenced by proceedStep().

FloatArray oofem::NonLinearDynamic::rhs2
protected

Definition at line 101 of file nlineardynamic.h.

Referenced by proceedStep(), and updateComponent().

LinSystSolverType oofem::NonLinearDynamic::solverType
protected

Definition at line 90 of file nlineardynamic.h.

Referenced by initializeFrom().

SparseMtrxType oofem::NonLinearDynamic::sparseMtrxType
protected

Definition at line 91 of file nlineardynamic.h.

Referenced by initializeFrom(), and proceedStep().

FloatArray oofem::NonLinearDynamic::totalDisplacement
protected
int oofem::NonLinearDynamic::totIterations
protected

Definition at line 106 of file nlineardynamic.h.

Referenced by proceedStep(), updateComponent(), and updateYourself().

FloatArray oofem::NonLinearDynamic::velocityVector
protected

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

This page is part of the OOFEM documentation. Copyright (c) 2011 Borek Patzak
Project e-mail: info@oofem.org
Generated at Tue Jan 2 2018 20:07:39 for OOFEM by doxygen 1.8.11 written by Dimitri van Heesch, © 1997-2011