Implementation of sparse nonlinear solver with indirect control. More...

#include <calmls.h>

Inheritance diagram for oofem::CylindricalALM:

Collaboration diagram for oofem::CylindricalALM:

Public Member Functions
	CylindricalALM (Domain d, EngngModel m)

virtual	~CylindricalALM ()

virtual NM_Status	solve (SparseMtrx &K, FloatArray &R, FloatArray R0, FloatArray &X, FloatArray &dX, FloatArray &F, const FloatArray &internalForcesEBENorm, double &ReachedLambda, referenceLoadInputModeType rlm, int &nite, TimeStep )
	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 IRResultType	initializeFrom (InputRecord *ir)

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

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

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

virtual void	setDomain (Domain *d)

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

virtual const char *	giveClassName () const

virtual const char *	giveInputRecordName () const

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

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

virtual	~SparseNonLinearSystemNM ()
	Destructor. More...

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

IRResultType	initializeFrom (InputRecord *ir)

virtual void	convertPertMap ()

virtual void	applyPerturbation (FloatArray *displacement)

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

Protected Types
enum	calm_ControlType { calm_hpc_off = 0, calm_hpc_on, calml_hpc }
	CALM mode type; determines the calm step length control. More...

enum	calm_NR_ModeType { calm_modifiedNRM, calm_fullNRM, calm_accelNRM }
	Controlling mode of Newton-Raphson Method. More...

typedef std::set< DofIDItem >	__DofIDSet

Protected Member Functions
void	convertHPCMap ()

int	computeDeltaLambda (double &deltaLambda, const FloatArray &dX, const FloatArray &deltaXt, const FloatArray &deltaX_, const FloatArray &R, double RR, double eta, double deltaL, double DeltaLambda0, int neq)

void	search (int istep, FloatArray &prod, FloatArray &eta, double amp, double maxeta, double mineta, int &status)

bool	checkConvergence (const FloatArray &R, const FloatArray *R0, const FloatArray &F, const FloatArray &X, const FloatArray &ddX, double Lambda, double RR0, double RR, double drProduct, const FloatArray &internalForcesEBENorm, int nite, bool &errorOutOfRange)
	Evaluates the convergence criteria. More...

void	do_lineSearch (FloatArray &X, const FloatArray &XInitial, const FloatArray &deltaX_, const FloatArray &deltaXt, const FloatArray &dXm1, FloatArray &dX, FloatArray &ddX, const FloatArray &R, const FloatArray R0, const FloatArray &F, double &DeltaLambda, double &DeltaLambdam1, double &deltaLambda, double &Lambda, double &ReachedLambda, double RR, double &drProduct, TimeStep tStep)
	Perform line search optimization of step length. More...

Protected Attributes
int	nsmax

int	maxRestarts

double	Psi

double	deltaL

double	minStepLength

double	maxStepLength

int	solved

int	numberOfRequiredIterations

calm_NR_ModeType	calm_NR_Mode

calm_NR_ModeType	calm_NR_OldMode

int	calm_NR_ModeTick

int	calm_MANRMSteps

int	minIterations
	Minimum hard number of iteration.s. More...

int	calm_hpc_init
	Variables for HyperPlaneControl. More...

calm_ControlType	calm_Control

FloatArray	calm_HPCWeights

IntArray	calm_HPCIndirectDofMask
	Array containing equation numbers of dofs under indirect control. More...

IntArray	calm_HPCDmanDofSrcArray
	Input array containing dofmanagers and corresponding dof numbers under indirect control. More...

FloatArray	calm_HPCDmanWeightSrcArray
	Input array of dofman weights (for hpcmode 2). More...

std::unique_ptr< SparseLinearSystemNM >	linSolver
	Linear system solver. More...

LinSystSolverType	solverType
	linear system solver ID. More...

int	lsFlag
	Line search flag. More...

double	ls_tolerance
	Line search tolerance. More...

double	amplifFactor
	Line search amplification factor. More...

double	maxEta
	Line search parameters (limits). More...

double	minEta

int	nccdg
	Number of convergence criteria dof groups. More...

std::vector< __DofIDSet >	ccDofGroups
	Convergence criteria dof groups. More...

FloatArray	rtolf
	Relative unbalanced force tolerance for each group. More...

FloatArray	rtold
	Relative iterative displacement change tolerance for each group. More...

ParallelContext *	parallel_context
	Parallel context for computing norms, dot products and such. More...

Protected Attributes inherited from oofem::SparseNonLinearSystemNM
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...

Additional Inherited Members
Public Types inherited from oofem::SparseNonLinearSystemNM
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...

Detailed Description

Implementation of sparse nonlinear solver with indirect control.

It uses Cylindrical Arc-Length Method algorithm.

DESCRIPTION : Perform solution of non-linear problem in the form Kt * deltaR = g where g is defined as g = g(Lambda,X) = Lambda* R - F(X). During iteration process, we want g became zero vector.

This method uses Modified Newton Raphson iteration scheme

If we solve non-linear static we can interpret symbols as follows:

Kt - tangential stiffness dX - increment of displacements g - vector of unbalanced forces (at the end should be zero one) Lambda - Load level Rt - Quasi Total Load vector Rt = R0 + Lambda*R X - total displacement vector F(X) - Nodal representation of (real) internal forces. Psi - control parameter (=0. means displacement control) Bergan_k0 - variable used to compute bergan's parameters of current stiffness. calm_NR_Mode - variable controlling the mode of NRM (ModifiedNR, Full NRM (stiffness update after each iteration), Modified Accelerated NRM (we perform iteration with stiffness matrix updated only after calm_MANRMSteps) calm_NR_OldMode - variable containing the old mode of NRM, which will be restored after calm_NR_ModeTick iterations. calm_NR_ModeTick - see calm_NR_OldMode. calm_MANRMSteps - if calm_NR_Mode == calm_accelNRM, it specifies, that new updated stiffness matrix is assembled after calm_MANRMSteps. calm_Control - variable indicating the ALM control. calm_HPCIndirectDofMask - Mask, telling which dofs are used for HPC. calm_HPCWeights - dofs weights in constrain.

Tasks:

solving problem: solveYourselfAt.
returning results (increment of displacement, reached level of loading and so on)

Variable description :

K(N,N) = Stiffness matrix (assumed positive definite)
dX(N) = Iterative increment of displacement
Rt = "quasi" total load vector
R0 = Fixed (Initial) load vector
deltaL = Step length
deltaLambda = Increment of load level
Psi = Determines loading or displacement control
dX = Current total increment
F = Nodal representation of (real) internal forces
DeltaLambda = Current total increment of load level
Lambda = Total load level
RTOL= Convergence tolerance
NSMAX = Maximum number of sweeps allowed

OUTPUT : (after call solveYourselfAt)

dX = Reached displacement increment
Lambda = Total load level
nite = Number of iterations required to fulfill balance
status = NM_status with flags set to reached state

Definition at line 144 of file calmls.h.

Member Typedef Documentation

typedef std :: set< DofIDItem > oofem::CylindricalALM::__DofIDSet

protected

Definition at line 161 of file calmls.h.

Member Enumeration Documentation

enum oofem::CylindricalALM::calm_ControlType

protected

CALM mode type; determines the calm step length control.

Enumerator
calm_hpc_off	Full ALM with quadratic constrain and all dofs.
calm_hpc_on	Full ALM with quadratic constrain, taking into account only selected dofs.
calml_hpc	Linearized ALM (only displacements), taking into account only selected dofs with given weight.

Definition at line 148 of file calmls.h.

enum oofem::CylindricalALM::calm_NR_ModeType

protected

Controlling mode of Newton-Raphson Method.

Enumerator
calm_modifiedNRM	Keeps the old tangent.
calm_fullNRM	Updates the tangent every iteration.
calm_accelNRM	Updates the tangent after a few steps.

Definition at line 155 of file calmls.h.

Constructor & Destructor Documentation

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

oofem::CylindricalALM::~CylindricalALM ( )

virtual

Definition at line 103 of file calmls.C.

Member Function Documentation

bool oofem::CylindricalALM::checkConvergence	(	const FloatArray &	R,
		const FloatArray *	R0,
		const FloatArray &	F,
		const FloatArray &	X,
		const FloatArray &	ddX,
		double	Lambda,
		double	RR0,
		double	RR,
		double	drProduct,
		const FloatArray &	internalForcesEBENorm,
		int	nite,
		bool &	errorOutOfRange
	)

protected

Evaluates the convergence criteria.

Definition at line 478 of file calmls.C.

References oofem::ParallelContext::accumulate(), oofem::FloatArray::add(), oofem::FloatArray::at(), CALM_MAX_REL_ERROR_BOUND, calm_SMALL_ERROR_NUM, ccDofGroups, oofem::NumericalMethod::domain, oofem::Element_local, oofem::Domain::giveDofManagers(), oofem::Domain::giveElements(), oofem::ParallelContext::localNorm(), nccdg, OOFEM_LOG_INFO, parallel_context, rtold, rtolf, oofem::FloatArray::subtract(), oofem::FloatArray::sum(), oofem::FloatArray::times(), and oofem::FloatArray::zero().

Referenced by solve().

int oofem::CylindricalALM::computeDeltaLambda	(	double &	deltaLambda,
		const FloatArray &	dX,
		const FloatArray &	deltaXt,
		const FloatArray &	deltaX_,
		const FloatArray &	R,
		double	RR,
		double	eta,
		double	deltaL,
		double	DeltaLambda0,
		int	neq
	)

protected

Returns: If 0 then ok, 1 then failure (restart).

Definition at line 987 of file calmls.C.

References oofem::ParallelContext::accumulate(), oofem::IntArray::at(), oofem::FloatArray::at(), calm_Control, calm_hpc_off, calm_hpc_on, calm_HPCIndirectDofMask, calm_HPCWeights, calm_SMALL_NUM, calml_hpc, deltaL, oofem::IntArray::giveSize(), oofem::ParallelContext::localDotProduct(), oofem::ParallelContext::localNorm(), OOFEM_ERROR, parallel_context, and Psi.

Referenced by do_lineSearch(), and solve().

void oofem::CylindricalALM::convertHPCMap ( )

protected

Definition at line 886 of file calmls.C.

References oofem::IntArray::at(), oofem::FloatArray::at(), calm_Control, calm_HPCDmanDofSrcArray, calm_HPCDmanWeightSrcArray, calm_HPCIndirectDofMask, calm_HPCWeights, calml_hpc, oofem::NumericalMethod::domain, oofem::Domain::giveDofManagers(), oofem::IntArray::giveSize(), oofem::ParallelContext::isLocal(), OOFEM_WARNING, parallel_context, oofem::IntArray::resize(), and oofem::FloatArray::resize().

Referenced by solve().

void oofem::CylindricalALM::do_lineSearch	(	FloatArray &	X,
		const FloatArray &	XInitial,
		const FloatArray &	deltaX_,
		const FloatArray &	deltaXt,
		const FloatArray &	dXm1,
		FloatArray &	dX,
		FloatArray &	ddX,
		const FloatArray &	R,
		const FloatArray *	R0,
		const FloatArray &	F,
		double &	DeltaLambda,
		double &	DeltaLambdam1,
		double &	deltaLambda,
		double &	Lambda,
		double &	ReachedLambda,
		double	RR,
		double &	drProduct,
		TimeStep *	tStep
	)

protected

Perform line search optimization of step length.

Definition at line 1204 of file calmls.C.

References oofem::FloatArray::add(), amplifFactor, oofem::FloatArray::at(), oofem::FloatArray::clear(), computeDeltaLambda(), deltaL, oofem::NumericalMethod::domain, oofem::NumericalMethod::engngModel, oofem::FloatArray::giveSize(), oofem::TimeStep::incrementStateCounter(), oofem::InternalRhs, oofem::ParallelContext::localDotProduct(), oofem::ParallelContext::localNorm(), ls_tolerance, maxEta, minEta, OOFEM_LOG_INFO, parallel_context, search(), and oofem::EngngModel::updateComponent().

Referenced by solve().

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

inlinevirtual

Reimplemented from oofem::SparseNonLinearSystemNM.

Definition at line 239 of file calmls.h.

virtual double oofem::CylindricalALM::giveCurrentStepLength ( )

inlinevirtual

Returns step length.

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

Returns: Current step length

Reimplemented from oofem::SparseNonLinearSystemNM.

Definition at line 221 of file calmls.h.

virtual const char* oofem::CylindricalALM::giveInputRecordName ( ) const

inlinevirtual

Definition at line 240 of file calmls.h.

References _IFT_CylindricalALM_Name.

SparseLinearSystemNM * oofem::CylindricalALM::giveLinearSolver ( )

virtual

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 from oofem::SparseNonLinearSystemNM.

Definition at line 969 of file calmls.C.

References oofem::classFactory, oofem::ClassFactory::createSparseLinSolver(), oofem::NumericalMethod::domain, oofem::NumericalMethod::engngModel, linSolver, OOFEM_ERROR, and solverType.

Referenced by initializeFrom(), and solve().

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

virtual

Reimplemented from oofem::NumericalMethod.

Definition at line 675 of file calmls.C.

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

inlinevirtual

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

arc length methods).

Reimplemented from oofem::SparseNonLinearSystemNM.

Definition at line 224 of file calmls.h.

virtual void oofem::CylindricalALM::reinitialize ( )

inlinevirtual

Reinitializes the receiver.

This is used, when topology of problem has changed (for example after adaptive refinement or load transfer in parallel applications). This is necessary for numerical methods, that cache some data between solution steps and that may depend on domain or problem topology. The caching of data by receiver is intended only for speeding up the calculation, but numerical method must be always able to generate this data again. This method clears receiver cached data dependent on topology, when it changes.

Reimplemented from oofem::NumericalMethod.

Definition at line 233 of file calmls.h.

contextIOResultType oofem::CylindricalALM::restoreContext	(	DataStream &	stream,
		ContextMode	mode,
		void *	obj = `NULL`
	)

virtual

Reimplemented from oofem::NumericalMethod.

Definition at line 957 of file calmls.C.

References oofem::CIO_IOERR, oofem::CIO_OK, deltaL, oofem::DataStream::read(), and THROW_CIOERR.

contextIOResultType oofem::CylindricalALM::saveContext	(	DataStream &	stream,
		ContextMode	mode,
		void *	obj = `NULL`
	)

virtual

Reimplemented from oofem::NumericalMethod.

Definition at line 945 of file calmls.C.

References oofem::CIO_IOERR, oofem::CIO_OK, deltaL, THROW_CIOERR, and oofem::DataStream::write().

void oofem::CylindricalALM::search	(	int	istep,
		FloatArray &	prod,
		FloatArray &	eta,
		double	amp,
		double	maxeta,
		double	mineta,
		int &	status
	)

protected

Definition at line 1118 of file calmls.C.

References oofem::FloatArray::at(), and oofem::max().

Referenced by do_lineSearch().

virtual void oofem::CylindricalALM::setDomain ( Domain * d )

inlinevirtual

Reimplemented from oofem::NumericalMethod.

Definition at line 227 of file calmls.h.

virtual void oofem::CylindricalALM::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 from oofem::SparseNonLinearSystemNM.

Definition at line 222 of file calmls.h.

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

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.