This base class is an abstraction for all numerical methods solving sparse nonlinear system of equations. More...

#include <sparsenonlinsystemnm.h>

Inheritance diagram for oofem::SparseNonLinearSystemNM:

Collaboration diagram for oofem::SparseNonLinearSystemNM:

Public Types
enum	referenceLoadInputModeType { rlm_total =0, rlm_incremental =1 }
	The following parameter allows to specify how the reference load vector is obtained from given totalLoadVector and initialLoadVector. More...

Public Member Functions
	SparseNonLinearSystemNM (Domain d, EngngModel m)
	Constructor. More...

virtual	~SparseNonLinearSystemNM ()
	Destructor. More...

virtual NM_Status	solve (SparseMtrx &K, FloatArray &R, FloatArray R0, FloatArray &X, FloatArray &dX, FloatArray &F, const FloatArray &internalForcesEBENorm, double &s, referenceLoadInputModeType rlm, int &nite, TimeStep tStep)=0
	Solves the given sparse linear system of equations . More...

virtual double	giveCurrentStepLength ()
	Returns step length. More...

virtual void	setStepLength (double s)
	Sets the step length. More...

virtual bool	referenceLoad () const
	Returns true if reference loads are used (i.e. More...

virtual void	printState (FILE *outputStream)
	Prints status message of receiver to output stream. More...

virtual SparseLinearSystemNM *	giveLinearSolver ()
	Constructs (if necessary) and returns a linear solver. More...

IRResultType	initializeFrom (InputRecord *ir)

virtual void	convertPertMap ()

virtual void	applyPerturbation (FloatArray *displacement)

virtual const char *	giveClassName () const

std::string	errorInfo (const char *func) const
	Error printing helper. More...

Public Member Functions inherited from oofem::NumericalMethod
	NumericalMethod (Domain d, EngngModel m)
	Constructor. More...

virtual	~NumericalMethod ()
	Destructor. More...

EngngModel *	giveEngngModel ()

virtual void	reinitialize ()
	Reinitializes the receiver. More...

virtual void	setDomain (Domain *d)

virtual contextIOResultType	saveContext (DataStream &stream, ContextMode mode, void *obj=NULL)

virtual contextIOResultType	restoreContext (DataStream &stream, ContextMode mode, void *obj=NULL)

Protected Attributes
double	deltaL
	Load level. More...

double	randPertAmplitude
	Amplitude of a random perturbation applied on the solution before the iteration process. More...

int	randSeed

bool	pert_init_needed

IntArray	igp_PertDmanDofSrcArray

FloatArray	igp_PertWeightArray

IntArray	igp_Map

FloatArray	igp_Weight

Protected Attributes inherited from oofem::NumericalMethod
Domain *	domain
	Pointer to domain. More...

EngngModel *	engngModel
	Pointer to engineering model. More...

Detailed Description

This base class is an abstraction for all numerical methods solving sparse nonlinear system of equations.

The purpose of this class is to declare the general interface to all numerical methods solving this kind of problem. This interface allows to use any suitable instance of the Numerical method class to the solve problem, and leave the whole engineering model code, including mapping, unchanged, because all instances of this class provide the common interface.

Definition at line 74 of file sparsenonlinsystemnm.h.

Member Enumeration Documentation

enum oofem::SparseNonLinearSystemNM::referenceLoadInputModeType

The following parameter allows to specify how the reference load vector is obtained from given totalLoadVector and initialLoadVector.

The initialLoadVector describes the part of loading which does not scale.

Enumerator
rlm_total	The reference incremental load vector is defined as totalLoadVector assembled at given time.
rlm_incremental	The reference load vector is obtained as incremental load vector at given time.

Definition at line 94 of file sparsenonlinsystemnm.h.

Constructor & Destructor Documentation

oofem::SparseNonLinearSystemNM::SparseNonLinearSystemNM	(	Domain *	d,
		EngngModel *	m
	)

inline

Constructor.

Definition at line 100 of file sparsenonlinsystemnm.h.

virtual oofem::SparseNonLinearSystemNM::~SparseNonLinearSystemNM ( )

inlinevirtual

Destructor.

Definition at line 102 of file sparsenonlinsystemnm.h.

Member Function Documentation

void oofem::SparseNonLinearSystemNM::applyPerturbation ( FloatArray * displacement )

virtual

Definition at line 103 of file sparsenonlinsystemnm.C.

References oofem::IntArray::at(), oofem::FloatArray::at(), convertPertMap(), oofem::FloatArray::giveSize(), igp_Map, igp_Weight, pert_init_needed, randPertAmplitude, and randSeed.

Referenced by oofem::CylindricalALM::solve().

void oofem::SparseNonLinearSystemNM::convertPertMap ( )

virtual

Definition at line 80 of file sparsenonlinsystemnm.C.

References oofem::IntArray::at(), oofem::FloatArray::at(), oofem::NumericalMethod::domain, oofem::DofManager::giveDofWithID(), oofem::Dof::giveEquationNumber(), oofem::DofManager::giveLabel(), oofem::Domain::giveNode(), oofem::IntArray::giveSize(), igp_Map, igp_PertDmanDofSrcArray, igp_PertWeightArray, igp_Weight, oofem::IntArray::resize(), and oofem::FloatArray::resize().

Referenced by applyPerturbation().

std :: string oofem::SparseNonLinearSystemNM::errorInfo ( const char * func ) const

inline

Error printing helper.

Definition at line 161 of file sparsenonlinsystemnm.h.

virtual const char* oofem::SparseNonLinearSystemNM::giveClassName ( ) const

inlinevirtual

Reimplemented in oofem::CylindricalALM, oofem::NRSolver, oofem::StaggeredSolver, and oofem::DynamicRelaxationSolver.

Definition at line 159 of file sparsenonlinsystemnm.h.

virtual double oofem::SparseNonLinearSystemNM::giveCurrentStepLength ( )

inlinevirtual

Returns step length.

See also: solve For more details on the step length s.

Returns: Current step length

Reimplemented in oofem::CylindricalALM.

Definition at line 131 of file sparsenonlinsystemnm.h.

Referenced by oofem::AdaptiveNonLinearStatic::adaptiveRemap(), and oofem::AdaptiveNonLinearStatic::initializeAdaptiveFrom().

virtual SparseLinearSystemNM* oofem::SparseNonLinearSystemNM::giveLinearSolver ( )

inlinevirtual

Constructs (if necessary) and returns a linear solver.

Public method because some problems require it for sensitivity analysis, etc. even for nonlinear problems (e.g. tangent relations in multiscale simulations).

Reimplemented in oofem::CylindricalALM, and oofem::NRSolver.

Definition at line 153 of file sparsenonlinsystemnm.h.

Referenced by oofem::NonLinearStatic::proceedStep().

IRResultType oofem::SparseNonLinearSystemNM::initializeFrom ( InputRecord * ir )

virtual

Reimplemented from oofem::NumericalMethod.

Reimplemented in oofem::StaggeredSolver.

Definition at line 44 of file sparsenonlinsystemnm.C.

References _IFT_NonLinearStatic_pert, _IFT_NonLinearStatic_pertw, _IFT_NonLinearStatic_randPertAmplitude, _IFT_NonLinearStatic_randSeed, oofem::IntArray::clear(), oofem::FloatArray::clear(), oofem::IntArray::giveSize(), oofem::FloatArray::giveSize(), igp_PertDmanDofSrcArray, igp_PertWeightArray, IR_GIVE_OPTIONAL_FIELD, oofem::IRRT_BAD_FORMAT, oofem::IRRT_OK, OOFEM_WARNING, pert_init_needed, randPertAmplitude, and randSeed.

Referenced by oofem::NRSolver::initializeFrom(), and oofem::CylindricalALM::initializeFrom().

virtual void oofem::SparseNonLinearSystemNM::printState ( FILE * outputStream )

inlinevirtual

Prints status message of receiver to output stream.

Prints the message corresponding to last solve.

Parameters

outputStream Stream to print state to.

Reimplemented in oofem::NRSolver.

Definition at line 147 of file sparsenonlinsystemnm.h.

Referenced by oofem::NonLinearStatic::printOutputAt().

virtual bool oofem::SparseNonLinearSystemNM::referenceLoad ( ) const

inlinevirtual

Returns true if reference loads are used (i.e.

arc length methods).

Reimplemented in oofem::CylindricalALM.

Definition at line 141 of file sparsenonlinsystemnm.h.

virtual void oofem::SparseNonLinearSystemNM::setStepLength ( double s )

inlinevirtual

Sets the step length.

See also: solve For more details on the step length s.

Parameters

s	New step length.

Reimplemented in oofem::CylindricalALM.

Definition at line 137 of file sparsenonlinsystemnm.h.

Referenced by oofem::AdaptiveNonLinearStatic::adaptiveRemap(), and oofem::AdaptiveNonLinearStatic::initializeAdaptiveFrom().

virtual NM_Status oofem::SparseNonLinearSystemNM::solve	(	SparseMtrx &	K,
		FloatArray &	R,
		FloatArray *	R0,
		FloatArray &	X,
		FloatArray &	dX,
		FloatArray &	F,
		const FloatArray &	internalForcesEBENorm,
		double &	s,
		referenceLoadInputModeType	rlm,
		int &	nite,
		TimeStep *	tStep
	)

pure virtual

Solves the given sparse linear system of equations $ s R + R_0 - F(X) = 0 $ .

Total load vector is not passed, it is defined as $ s R + R_0 $ , where s is scale factor.

See also: EngngModel::updateComponent Used to update the stiffness matrix and load vector.

Parameters

K	Coefficient matrix ( $\displaystyle K = \frac{\partial F}{\partial X}$ ; stiffness matrix).
R	Reference incremental RHS (incremental load).
R0	Initial RHS (initial load).
X	Total solution (total displacement).
dX	Increment of solution (incremental displacements).
F	InternalRhs (real internal forces).
internalForcesEBENorm	Norm of internal nodal forces (evaluated on element by element basis) (split into each DOF id).
s	RHS scale factor (load level).
rlm	Reference load mode.
nite	Number of iterations needed.
tStep	Time step to solve for.