#include <nlineardynamic.h>

Inheritance diagram for oofem::NonLinearDynamic:

Collaboration diagram for oofem::NonLinearDynamic:

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

Protected Member Functions
void	assemble (SparseMtrx &answer, TimeStep tStep, const MatrixAssembler &ma, const UnknownNumberingScheme &, Domain domain) override
void	proceedStep (int di, TimeStep *tStep)
void	determineConstants (TimeStep *tStep)
void	packMigratingData (TimeStep *tStep) override
void	unpackMigratingData (TimeStep *tStep) override
Protected Member Functions inherited from oofem::StructuralEngngModel
void	printReactionForces (TimeStep tStep, int id, FILE out)
virtual void	computeExternalLoadReactionContribution (FloatArray &reactions, TimeStep *tStep, int di)
void	updateInternalState (TimeStep *tStep)
void	printOutputAt (FILE file, TimeStep tStep) override

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

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

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

◆ NonLinearDynamic()

oofem::NonLinearDynamic::NonLinearDynamic	(	int	i,
		EngngModel *	master = nullptr )

Definition at line 69 of file nlineardynamic.C.

References commInitFlag, forcesVector, incrementOfDisplacement, initFlag, initialTimeDiscretization, internalForces, oofem::StructuralEngngModel::internalVarUpdateStamp, oofem::EngngModel::ndomains, oofem::StructuralEngngModel::StructuralEngngModel(), oofem::TD_ThreePointBackward, and totalDisplacement.

◆ ~NonLinearDynamic()

oofem::NonLinearDynamic::~NonLinearDynamic ( )

virtual

Definition at line 80 of file nlineardynamic.C.

Member Function Documentation

◆ assemble()

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

overrideprotectedvirtual

Assembles characteristic matrix of required type into given sparse matrix.

Parameters

answer	Assembled matrix.
tStep	Time step, when answer is assembled.
s	Determines the equation numbering scheme.
type	Characteristic components of type type are requested from elements and assembled.
domain	Source 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 796 of file nlineardynamic.C.

References oofem::StructuralElement::addNonlocalStiffnessContributions(), oofem::Domain::giveElements(), nonlocalStiffnessFlag, OOFEM_LOG_DEBUG, oofem::SparseMtrx::printStatistics(), and oofem::EngngModel::timer.

◆ determineConstants()

void oofem::NonLinearDynamic::determineConstants ( TimeStep * tStep )

protected

Definition at line 488 of file nlineardynamic.C.

References a0, a1, a2, a3, a4, a5, a6, beta, deltaT, gamma, oofem::TimeStep::giveTimeDiscretization(), oofem::TimeStep::giveTimeIncrement(), oofem::TimeStep::isTheFirstStep(), OOFEM_ERROR, OOFEM_LOG_DEBUG, oofem::TD_Newmark, oofem::TD_ThreePointBackward, and oofem::TD_TwoPointBackward.

Referenced by proceedStep().

◆ estimateMaxPackSize()

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

overridevirtual

Determines the space necessary for send/receive buffer. It uses related communication map pattern to determine the maximum size needed.

Parameters

commMap	Communication map used to send/receive messages.
buff	Communication buffer.
packUnpackType	Determines 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 907 of file nlineardynamic.C.

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

◆ giveClassName()

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

inlineoverridevirtual

Returns class name of the receiver.

Implements oofem::EngngModel.

Definition at line 145 of file nlineardynamic.h.

◆ giveFormulation()

fMode oofem::NonLinearDynamic::giveFormulation ( )

inlineoverridevirtual

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

Reimplemented from oofem::EngngModel.

Definition at line 146 of file nlineardynamic.h.

References oofem::EngngModel::nonLinFormulation.

◆ giveInitialTimeDiscretization()

TimeDiscretizationType oofem::NonLinearDynamic::giveInitialTimeDiscretization ( )

inline

Definition at line 151 of file nlineardynamic.h.

References initialTimeDiscretization.

Referenced by oofem::FluidStructureProblem::initializeYourself().

◆ giveInputRecordName()

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

inline

Definition at line 144 of file nlineardynamic.h.

References _IFT_NonLinearDynamic_Name.

◆ giveLoadBalancer()

LoadBalancer * oofem::NonLinearDynamic::giveLoadBalancer ( )

overridevirtual

Returns reference to receiver's load balancer.

Reimplemented from oofem::EngngModel.

Definition at line 939 of file nlineardynamic.C.

References oofem::classFactory, oofem::EngngModel::giveDomain(), oofem::EngngModel::lb, and oofem::EngngModel::loadBalancingFlag.

Referenced by updateDomainLinks().

◆ giveLoadBalancerMonitor()

LoadBalancerMonitor * oofem::NonLinearDynamic::giveLoadBalancerMonitor ( )

overridevirtual

Returns reference to receiver's load balancer monitor.

Reimplemented from oofem::EngngModel.

Definition at line 955 of file nlineardynamic.C.

References oofem::classFactory, oofem::EngngModel::lbm, and oofem::EngngModel::loadBalancingFlag.

◆ giveNextStep()

TimeStep * oofem::NonLinearDynamic::giveNextStep ( )

overridevirtual

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

Reimplemented from oofem::EngngModel.

Definition at line 213 of file nlineardynamic.C.

References oofem::EngngModel::currentStep, deltaT, oofem::EngngModel::giveMetaStep(), oofem::EngngModel::giveNumberOfFirstStep(), initialTimeDiscretization, oofem::microScale, oofem::EngngModel::nMetaSteps, OOFEM_ERROR, oofem::EngngModel::previousStep, and oofem::EngngModel::pScale.

◆ giveNumericalMethod()

NumericalMethod * oofem::NonLinearDynamic::giveNumericalMethod ( MetaStep * mStep )

overridevirtual

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

◆ giveUnknownComponent()

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

overridevirtual

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 181 of file nlineardynamic.C.

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

◆ giveUnknownDictHashIndx()

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

inlineoverridevirtual

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

◆ initializeFrom()

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

overridevirtual

Initializes receiver according to object description in input reader. InitString can be imagined as data record in component database belonging to receiver. Receiver may use value-name extracting functions to extract particular field from record.

Reimplemented from oofem::EngngModel.

Definition at line 123 of file nlineardynamic.C.

References _IFT_EngngModel_lstype, _IFT_EngngModel_smtype, _IFT_NonLinearDynamic_beta, _IFT_NonLinearDynamic_ddtScheme, _IFT_NonLinearDynamic_deltat, _IFT_NonLinearDynamic_gamma, _IFT_NonLinearDynamic_nonlocalext, _IFT_NonLinearDynamic_nonlocstiff, _IFT_NRSolver_manrmsteps, beta, oofem::CBT_static, oofem::EngngModel::commBuff, oofem::EngngModel::communicator, deltaT, gamma, oofem::EngngModel::giveNumberOfProcesses(), oofem::EngngModel::giveRank(), oofem::InputRecord::hasField(), initialTimeDiscretization, IR_GIVE_OPTIONAL_FIELD, oofem::EngngModel::isParallel(), MANRMSteps, oofem::EngngModel::nonlocalExt, nonlocalStiffnessFlag, oofem::EngngModel::nonlocCommunicator, OOFEM_ERROR, OOFEM_LOG_INFO, solverType, sparseMtrxType, oofem::TD_Newmark, oofem::TD_ThreePointBackward, and oofem::TD_TwoPointBackward.

◆ initializeYourself()

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

overridevirtual

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 292 of file nlineardynamic.C.

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

◆ packMigratingData()

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

overrideprotectedvirtual

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 971 of file nlineardynamic.C.

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

◆ printDofOutputAt()

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

overridevirtual

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

stream	output stream
iDof	dof to be processed
tStep	solution step

Reimplemented from oofem::EngngModel.

Definition at line 722 of file nlineardynamic.C.

References oofem::Dof::printMultipleOutputAt().

◆ printOutputAt()

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

overridevirtual

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

Definition at line 705 of file nlineardynamic.C.

References currentIterations, oofem::EngngModel::giveDomain(), oofem::TimeStep::giveNumber(), oofem::TimeStep::giveTargetTime(), nMethod, and oofem::StructuralEngngModel::printReactionForces().

◆ proceedStep()

void oofem::NonLinearDynamic::proceedStep	(	int	di,
		TimeStep *	tStep )

protected

Definition at line 345 of file nlineardynamic.C.

References a0, a1, a2, a3, a4, a5, a6, accelerationVector, oofem::assemble(), oofem::EngngModel::assembleVector(), oofem::classFactory, oofem::TimeStep::convergedReason, oofem::CR_CONVERGED, currentIterations, delta, determineConstants(), effectiveStiffnessMatrix, eta, forcesVector, oofem::EngngModel::giveCurrentMetaStep(), oofem::EngngModel::giveDomain(), oofem::EngngModel::giveNumberOfDomainEquations(), giveNumericalMethod(), help, incrementOfDisplacement, initFlag, internalForces, oofem::StructuralEngngModel::internalForcesEBENorm, oofem::EngngModel::LoadExchangeTag, massMatrix, nMethod, nonlocalStiffnessFlag, oofem::TimeStep::numberOfIterations, OOFEM_ERROR, OOFEM_LOG_DEBUG, previousAccelerationVector, previousIncrementOfDisplacement, previousInternalForces, previousTotalDisplacement, previousVelocityVector, oofem::FloatArray::resize(), rhs, rhs2, oofem::SMT_CompCol, sparseMtrxType, timesMtrx(), totalDisplacement, totIterations, oofem::EngngModel::updateSharedDofManagers(), velocityVector, and oofem::FloatArray::zero().

Referenced by solveYourselfAt().

◆ restoreContext()

void oofem::NonLinearDynamic::restoreContext	(	DataStream &	stream,
		ContextMode	mode )

overridevirtual

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

stream	Context file.
mode	Determines amount of info in stream.

Exceptions

ContextIOERR exception if error encountered.

Reimplemented from oofem::EngngModel.

Definition at line 756 of file nlineardynamic.C.

References accelerationVector, oofem::CIO_OK, incrementOfDisplacement, THROW_CIOERR, totalDisplacement, and velocityVector.

◆ saveContext()

void oofem::NonLinearDynamic::saveContext	(	DataStream &	stream,
		ContextMode	mode )

overridevirtual

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

stream	Context stream.
mode	Determines amount of info in stream.

Exceptions

ContextIOERR If error encountered.

Reimplemented from oofem::EngngModel.

Definition at line 732 of file nlineardynamic.C.

References accelerationVector, oofem::CIO_OK, incrementOfDisplacement, THROW_CIOERR, totalDisplacement, and velocityVector.

◆ showSparseMtrxStructure()

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

overridevirtual

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

Reimplemented from oofem::EngngModel.

Definition at line 824 of file nlineardynamic.C.

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

◆ solveYourself()

void oofem::NonLinearDynamic::solveYourself ( )

overridevirtual

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

Reimplemented from oofem::EngngModel.

Definition at line 249 of file nlineardynamic.C.

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

◆ solveYourselfAt()

void oofem::NonLinearDynamic::solveYourselfAt ( TimeStep * tStep )

overridevirtual

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

◆ timesMtrx()

void oofem::NonLinearDynamic::timesMtrx	(	FloatArray &	answer,
		FloatArray &	vec,
		CharType	type,
		Domain *	domain,
		TimeStep *	tStep )

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.

Definition at line 847 of file nlineardynamic.C.

References oofem::FloatArray::at(), oofem::FloatMatrix::at(), oofem::IntArray::at(), oofem::Element_remote, oofem::Element::giveCharacteristicMatrix(), oofem::Domain::giveElement(), oofem::Element::giveLocationArray(), oofem::Domain::giveNumber(), oofem::EngngModel::giveNumberOfDomainEquations(), oofem::Domain::giveNumberOfElements(), oofem::FloatMatrix::giveNumberOfRows(), oofem::Element::giveParallelMode(), oofem::Element::giveRotationMatrix(), oofem::IntArray::giveSize(), oofem::FloatMatrix::isNotEmpty(), oofem::EngngModel::MassExchangeTag, OOFEM_ERROR, oofem::FloatArray::resize(), oofem::FloatMatrix::rotatedWith(), oofem::EngngModel::updateSharedDofManagers(), and oofem::FloatArray::zero().

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

◆ unpackMigratingData()

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

overrideprotectedvirtual

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 990 of file nlineardynamic.C.

References oofem::DofManager_shared, oofem::Domain::giveDofManager(), oofem::EngngModel::giveDomain(), oofem::EngngModel::giveDomainErrorEstimator(), oofem::Domain::giveNumberOfDofManagers(), oofem::EngngModel::giveNumberOfDomainEquations(), oofem::DofManager::giveParallelMode(), incrementOfDisplacement, initFlag, oofem::EngngModel::initializeCommMaps(), nMethod, and totalDisplacement.

◆ updateAttributes()

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

overridevirtual

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

mStep Meta 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, and OOFEM_ERROR.

◆ updateComponent()

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

overridevirtual

Reimplemented from oofem::EngngModel.

Definition at line 625 of file nlineardynamic.C.

References a0, a1, oofem::assemble(), currentIterations, delta, effectiveStiffnessMatrix, eta, forcesVector, oofem::EngngModel::giveDomain(), help, incrementOfDisplacement, internalForces, oofem::StructuralEngngModel::internalForcesEBENorm, oofem::InternalRhs, MANRMSteps, massMatrix, oofem::NonLinearLhs, OOFEM_ERROR, OOFEM_LOG_DEBUG, previousInternalForces, rhs2, oofem::EngngModel::timer, timesMtrx(), totIterations, and oofem::StructuralEngngModel::updateInternalRHS().

◆ updateDomainLinks()

void oofem::NonLinearDynamic::updateDomainLinks ( )

overridevirtual

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 781 of file nlineardynamic.C.

References oofem::EngngModel::giveCurrentMetaStep(), oofem::EngngModel::giveDomain(), giveLoadBalancer(), and giveNumericalMethod().

◆ updateInternalRHS()

void oofem::NonLinearDynamic::updateInternalRHS	(	FloatArray &	answer,
		TimeStep *	tStep,
		Domain *	d,
		FloatArray *	eNorm )

overridevirtual

Updates the solution (guess) according to the new values. Callback for nonlinear solvers (e.g. Newton-Raphson).

Parameters

solutionVector	New solution.
tStep	Time when component is updated.
d	Domain.
eNorm	Optional per-element norm (for normalization).

Reimplemented from oofem::EngngModel.

Definition at line 558 of file nlineardynamic.C.

References a0, a1, currentIterations, delta, eta, forcesVector, oofem::EngngModel::giveDomain(), help, incrementOfDisplacement, initFlag, internalForces, oofem::TimeStep::isTheFirstStep(), massMatrix, OOFEM_LOG_DEBUG, previousInternalForces, rhs2, oofem::EngngModel::timer, timesMtrx(), totIterations, and oofem::StructuralEngngModel::updateInternalRHS().

Referenced by initializeYourself().

◆ updateMatrix()

void oofem::NonLinearDynamic::updateMatrix	(	SparseMtrx &	mat,
		TimeStep *	tStep,
		Domain *	d )

overridevirtual

Updates the solution (guess) according to the new values. Callback for nonlinear solvers (e.g. Newton-Raphson).

Note: For performance, the matrix should keep it's non-zero structure between calls, so the caller should make sure not to clear the matrix object before called.

Parameters

solutionVector	New solution.
tStep	Time when component is updated.
d	Domain.

Reimplemented from oofem::EngngModel.

Definition at line 596 of file nlineardynamic.C.

References a0, a1, oofem::SparseMtrx::add(), oofem::assemble(), delta, eta, MANRMSteps, massMatrix, OOFEM_LOG_DEBUG, oofem::EngngModel::timer, oofem::SparseMtrx::times(), totIterations, and oofem::SparseMtrx::zero().

◆ updateSolution()

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

overridevirtual

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

Parameters

solutionVector	New solution.
tStep	Time when component is updated.
d	Domain.

Reimplemented from oofem::EngngModel.

Definition at line 551 of file nlineardynamic.C.

◆ updateYourself()

void oofem::NonLinearDynamic::updateYourself ( TimeStep * tStep )

overridevirtual

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

Definition at line 537 of file nlineardynamic.C.

References accelerationVector, incrementOfDisplacement, internalForces, previousAccelerationVector, previousIncrementOfDisplacement, previousInternalForces, previousTotalDisplacement, previousVelocityVector, totalDisplacement, totIterations, and velocityVector.