A MPM explicit scheme. More...

#include <mpm_explicit_scheme.hpp>

Inheritance diagram for Kratos::MPMExplicitScheme< TSparseSpace, TDenseSpace >:

Collaboration diagram for Kratos::MPMExplicitScheme< TSparseSpace, TDenseSpace >:

Protected Attributes
Member Variables
ModelPart &	mr_grid_model_part

Protected Attributes inherited from Kratos::Scheme< TSparseSpace, TDenseSpace >
bool	mSchemeIsInitialized

bool	mElementsAreInitialized
	Flag to be used in controlling if the Scheme has been initialized or not. More...

bool	mConditionsAreInitialized
	Flag taking in account if the elements were initialized correctly or not. More...

Type Definitions
typedef Scheme< TSparseSpace, TDenseSpace >	BaseType

typedef BaseType::TDataType	TDataType

typedef BaseType::DofsArrayType	DofsArrayType

typedef Element::DofsVectorType	DofsVectorType

typedef BaseType::TSystemMatrixType	TSystemMatrixType

typedef BaseType::TSystemVectorType	TSystemVectorType

typedef BaseType::LocalSystemVectorType	LocalSystemVectorType

typedef BaseType::LocalSystemMatrixType	LocalSystemMatrixType

typedef ModelPart::ElementsContainerType	ElementsArrayType

typedef ModelPart::ConditionsContainerType	ConditionsArrayType

typedef BaseType::Pointer	BaseTypePointer

typedef ModelPart::NodesContainerType	NodesArrayType
	The arrays of elements and nodes. More...

typedef ModelPart::NodeIterator	NodeIterator
	Definition for the node iterator. More...

typedef std::size_t	SizeType
	Definition of the size type. More...

typedef std::size_t	IndexType
	Definition of the index type. More...

static constexpr double	numerical_limit = std::numeric_limits<double>::epsilon()
	The definition of the numerical limit. More...

	KRATOS_CLASS_POINTER_DEFINITION (MPMExplicitScheme)

	MPMExplicitScheme (ModelPart &grid_model_part)
	Default constructor. More...

virtual	~MPMExplicitScheme ()
	Destructor. More...

BaseTypePointer	Clone () override

void	Initialize (ModelPart &rModelPart) override
	This is the place to initialize the Scheme. This is intended to be called just once when the strategy is initialized. More...

void	Update (ModelPart &r_model_part, DofsArrayType &rDofSet, TSystemMatrixType &A, TSystemVectorType &Dx, TSystemVectorType &b) override

void	UpdateTranslationalDegreesOfFreedom (const ProcessInfo &r_current_process_info, NodeIterator itCurrentNode, const IndexType DisplacementPosition, const double delta_time, const SizeType DomainSize=3)

void	InitializeSolutionStep (ModelPart &r_model_part, TSystemMatrixType &A, TSystemVectorType &Dx, TSystemVectorType &b) override

void	calculateGridVelocityAndApplyDirichletBC (const ProcessInfo rCurrentProcessInfo, bool calculateVelocityFromMomenta=false)
	Apply Dirichlet BCs to nodal velocity field. More...

void	FinalizeSolutionStep (ModelPart &rModelPart, TSystemMatrixType &A, TSystemVectorType &Dx, TSystemVectorType &b) override

void	CalculateUpdatedGridVelocityField (const ProcessInfo &rCurrentProcessInfo, ModelPart &rModelPart)

void	GetDofList (const Element &rCurrentElement, Element::DofsVectorType &ElementalDofList, const ProcessInfo &CurrentProcessInfo) override

void	GetDofList (const Condition &rCurrentCondition, Element::DofsVectorType &rConditionDofList, const ProcessInfo &rCurrentProcessInfo) override

int	Check (const ModelPart &rModelPart) const override

void	CalculateRHSContribution (Element &rCurrentElement, LocalSystemVectorType &RHS_Contribution, Element::EquationIdVectorType &EquationId, const ProcessInfo &rCurrentProcessInfo) override
	This function is designed to calculate just the RHS contribution. More...

void	CalculateRHSContribution (Condition &rCurrentCondition, LocalSystemVectorType &RHS_Contribution, Element::EquationIdVectorType &EquationId, const ProcessInfo &rCurrentProcessInfo) override
	Functions that calculates the RHS of a "condition" object. More...

Additional Inherited Members
Public Types inherited from Kratos::Scheme< TSparseSpace, TDenseSpace >
using	ClassType = Scheme< TSparseSpace, TDenseSpace >
	The definition of the current class. More...

using	TDataType = typename TSparseSpace::DataType
	Data type definition. More...

using	TSystemMatrixType = typename TSparseSpace::MatrixType
	Matrix type definition. More...

using	TSystemVectorType = typename TSparseSpace::VectorType
	Vector type definition. More...

using	LocalSystemMatrixType = typename TDenseSpace::MatrixType
	Local system matrix type definition. More...

using	LocalSystemVectorType = typename TDenseSpace::VectorType
	Local system vector type definition. More...

using	TDofType = Dof< double >
	DoF type definition. More...

using	DofsArrayType = ModelPart::DofsArrayType
	DoF array type definition. More...

using	ElementsArrayType = ModelPart::ElementsContainerType
	Elements containers definition. More...

using	ConditionsArrayType = ModelPart::ConditionsContainerType
	Conditions containers definition. More...

Public Member Functions inherited from Kratos::Scheme< TSparseSpace, TDenseSpace >
	Scheme ()
	Default Constructor. More...

	Scheme (Parameters ThisParameters)
	Constructor with Parameters. More...

	Scheme (Scheme &rOther)

virtual	~Scheme ()

virtual std::string	Info () const
	Turn back information as a string. More...

virtual void	PrintInfo (std::ostream &rOStream) const
	Print information about this object. More...

virtual void	PrintData (std::ostream &rOStream) const
	Print object's data. More...

	KRATOS_CLASS_POINTER_DEFINITION (Scheme)
	Pointer definition of Scheme. More...

virtual ClassType::Pointer	Create (Parameters ThisParameters) const
	Create method. More...

bool	SchemeIsInitialized ()
	This method returns if the scheme is initialized. More...

void	SetSchemeIsInitialized (bool SchemeIsInitializedFlag=true)
	This method sets if the elements have been initialized or not (true by default) More...

bool	ElementsAreInitialized ()
	This method returns if the elements are initialized. More...

void	SetElementsAreInitialized (bool ElementsAreInitializedFlag=true)
	This method sets if the elements have been initialized or not (true by default) More...

bool	ConditionsAreInitialized ()
	This method returns if the conditions are initialized. More...

void	SetConditionsAreInitialized (bool ConditionsAreInitializedFlag=true)
	This method sets if the conditions have been initialized or not (true by default) More...

virtual void	InitializeElements (ModelPart &rModelPart)
	This is the place to initialize the elements. More...

virtual void	InitializeConditions (ModelPart &rModelPart)
	This is the place to initialize the conditions. More...

virtual void	InitializeNonLinIteration (ModelPart &rModelPart, TSystemMatrixType &A, TSystemVectorType &Dx, TSystemVectorType &b)
	unction to be called when it is needed to initialize an iteration. It is designed to be called at the beginning of each non linear iteration More...

virtual void	FinalizeNonLinIteration (ModelPart &rModelPart, TSystemMatrixType &A, TSystemVectorType &Dx, TSystemVectorType &b)
	Function to be called when it is needed to finalize an iteration. It is designed to be called at the end of each non linear iteration. More...

virtual void	Predict (ModelPart &rModelPart, DofsArrayType &rDofSet, TSystemMatrixType &A, TSystemVectorType &Dx, TSystemVectorType &b)
	Performing the prediction of the solution. More...

virtual void	CalculateOutputData (ModelPart &rModelPart, DofsArrayType &rDofSet, TSystemMatrixType &A, TSystemVectorType &Dx, TSystemVectorType &b)
	Functions to be called to prepare the data needed for the output of results. More...

virtual void	CleanOutputData ()
	Functions that cleans the results data. More...

virtual void	Clean ()
	This function is intended to be called at the end of the solution step to clean up memory storage not needed after the end of the solution step. More...

virtual void	Clear ()
	Liberate internal storage. More...

virtual int	Check (ModelPart &rModelPart)

virtual void	CalculateSystemContributions (Element &rElement, LocalSystemMatrixType &LHS_Contribution, LocalSystemVectorType &RHS_Contribution, Element::EquationIdVectorType &rEquationIdVector, const ProcessInfo &rCurrentProcessInfo)
	This function is designed to be called in the builder and solver to introduce the selected time integration scheme. More...

virtual void	CalculateSystemContributions (Condition &rCondition, LocalSystemMatrixType &LHS_Contribution, LocalSystemVectorType &RHS_Contribution, Element::EquationIdVectorType &rEquationIdVector, const ProcessInfo &rCurrentProcessInfo)
	Functions totally analogous to the precedent but applied to the "condition" objects. More...

virtual void	CalculateLHSContribution (Element &rElement, LocalSystemMatrixType &LHS_Contribution, Element::EquationIdVectorType &rEquationIdVector, const ProcessInfo &rCurrentProcessInfo)
	This function is designed to calculate just the LHS contribution. More...

virtual void	CalculateLHSContribution (Condition &rCondition, LocalSystemMatrixType &LHS_Contribution, Element::EquationIdVectorType &rEquationIdVector, const ProcessInfo &rCurrentProcessInfo)
	Functions totally analogous to the precedent but applied to the "condition" objects. More...

virtual void	EquationId (const Element &rElement, Element::EquationIdVectorType &rEquationId, const ProcessInfo &rCurrentProcessInfo)
	This method gets the eqaution id corresponding to the current element. More...

virtual void	EquationId (const Condition &rCondition, Element::EquationIdVectorType &rEquationId, const ProcessInfo &rCurrentProcessInfo)
	Functions totally analogous to the precedent but applied to the "condition" objects. More...

virtual Parameters	GetDefaultParameters () const
	This method provides the defaults parameters to avoid conflicts between the different constructors. More...

Static Public Member Functions inherited from Kratos::Scheme< TSparseSpace, TDenseSpace >
static std::string	Name ()
	Returns the name of the class as used in the settings (snake_case format) More...

Protected Member Functions inherited from Kratos::Scheme< TSparseSpace, TDenseSpace >
virtual Parameters	ValidateAndAssignParameters (Parameters ThisParameters, const Parameters DefaultParameters) const
	This method validate and assign default parameters. More...

virtual void	AssignSettings (const Parameters ThisParameters)
	This method assigns settings to member variables. More...

Detailed Description

template<class TSparseSpace, class TDenseSpace>
class Kratos::MPMExplicitScheme< TSparseSpace, TDenseSpace >

A MPM explicit scheme.

Scheme options include Forward Euler or Central Difference. Stress update options include Update Stress First (USF), Update Stress Last (USL) and Modified Update Stress Last (MUSL).

Author: Peter

Member Typedef Documentation

◆ BaseType

template<class TSparseSpace , class TDenseSpace >

typedef Scheme<TSparseSpace, TDenseSpace> Kratos::MPMExplicitScheme< TSparseSpace, TDenseSpace >::BaseType

◆ BaseTypePointer

template<class TSparseSpace , class TDenseSpace >

typedef BaseType::Pointer Kratos::MPMExplicitScheme< TSparseSpace, TDenseSpace >::BaseTypePointer

◆ ConditionsArrayType

template<class TSparseSpace , class TDenseSpace >

typedef ModelPart::ConditionsContainerType Kratos::MPMExplicitScheme< TSparseSpace, TDenseSpace >::ConditionsArrayType

◆ DofsArrayType

template<class TSparseSpace , class TDenseSpace >

typedef BaseType::DofsArrayType Kratos::MPMExplicitScheme< TSparseSpace, TDenseSpace >::DofsArrayType

◆ DofsVectorType

template<class TSparseSpace , class TDenseSpace >

typedef Element::DofsVectorType Kratos::MPMExplicitScheme< TSparseSpace, TDenseSpace >::DofsVectorType

◆ ElementsArrayType

template<class TSparseSpace , class TDenseSpace >

typedef ModelPart::ElementsContainerType Kratos::MPMExplicitScheme< TSparseSpace, TDenseSpace >::ElementsArrayType

◆ IndexType

template<class TSparseSpace , class TDenseSpace >

typedef std::size_t Kratos::MPMExplicitScheme< TSparseSpace, TDenseSpace >::IndexType

Definition of the index type.

◆ LocalSystemMatrixType

template<class TSparseSpace , class TDenseSpace >

typedef BaseType::LocalSystemMatrixType Kratos::MPMExplicitScheme< TSparseSpace, TDenseSpace >::LocalSystemMatrixType

◆ LocalSystemVectorType

template<class TSparseSpace , class TDenseSpace >

typedef BaseType::LocalSystemVectorType Kratos::MPMExplicitScheme< TSparseSpace, TDenseSpace >::LocalSystemVectorType

◆ NodeIterator

template<class TSparseSpace , class TDenseSpace >

typedef ModelPart::NodeIterator Kratos::MPMExplicitScheme< TSparseSpace, TDenseSpace >::NodeIterator

Definition for the node iterator.

◆ NodesArrayType

template<class TSparseSpace , class TDenseSpace >

typedef ModelPart::NodesContainerType Kratos::MPMExplicitScheme< TSparseSpace, TDenseSpace >::NodesArrayType

The arrays of elements and nodes.

◆ SizeType

template<class TSparseSpace , class TDenseSpace >

typedef std::size_t Kratos::MPMExplicitScheme< TSparseSpace, TDenseSpace >::SizeType

Definition of the size type.

◆ TDataType

template<class TSparseSpace , class TDenseSpace >

typedef BaseType::TDataType Kratos::MPMExplicitScheme< TSparseSpace, TDenseSpace >::TDataType

◆ TSystemMatrixType

template<class TSparseSpace , class TDenseSpace >

typedef BaseType::TSystemMatrixType Kratos::MPMExplicitScheme< TSparseSpace, TDenseSpace >::TSystemMatrixType

◆ TSystemVectorType

template<class TSparseSpace , class TDenseSpace >

typedef BaseType::TSystemVectorType Kratos::MPMExplicitScheme< TSparseSpace, TDenseSpace >::TSystemVectorType

Constructor & Destructor Documentation

◆ MPMExplicitScheme()

template<class TSparseSpace , class TDenseSpace >

Kratos::MPMExplicitScheme< TSparseSpace, TDenseSpace >::MPMExplicitScheme ( ModelPart & grid_model_part )

inlineexplicit

Default constructor.

The MPMExplicitScheme method

Parameters

MaximumDeltaTime	The maximum delta time to be considered
DeltaTimeFraction	The delta ttime fraction
DeltaTimePredictionLevel	The prediction level

◆ ~MPMExplicitScheme()

template<class TSparseSpace , class TDenseSpace >

virtual Kratos::MPMExplicitScheme< TSparseSpace, TDenseSpace >::~MPMExplicitScheme ( )

inlinevirtual

Destructor.

Member Function Documentation

◆ calculateGridVelocityAndApplyDirichletBC()

template<class TSparseSpace , class TDenseSpace >

void Kratos::MPMExplicitScheme< TSparseSpace, TDenseSpace >::calculateGridVelocityAndApplyDirichletBC	(	const ProcessInfo	rCurrentProcessInfo,
		bool	calculateVelocityFromMomenta = `false`
	)

inline

Apply Dirichlet BCs to nodal velocity field.

◆ CalculateRHSContribution() [1/2]

template<class TSparseSpace , class TDenseSpace >

void Kratos::MPMExplicitScheme< TSparseSpace, TDenseSpace >::CalculateRHSContribution	(	Condition &	rCurrentCondition,
		LocalSystemVectorType &	RHS_Contribution,
		Element::EquationIdVectorType &	EquationId,
		const ProcessInfo &	rCurrentProcessInfo
	)

inlineoverridevirtual

Functions that calculates the RHS of a "condition" object.

Parameters

pCondition	The condition to compute
RHS_Contribution	The RHS vector contribution
EquationId	The ID's of the condition degrees of freedom
rCurrentProcessInfo	The current process info instance

Reimplemented from Kratos::Scheme< TSparseSpace, TDenseSpace >.

◆ CalculateRHSContribution() [2/2]

template<class TSparseSpace , class TDenseSpace >

void Kratos::MPMExplicitScheme< TSparseSpace, TDenseSpace >::CalculateRHSContribution	(	Element &	rElement,
		LocalSystemVectorType &	RHS_Contribution,
		Element::EquationIdVectorType &	rEquationIdVector,
		const ProcessInfo &	rCurrentProcessInfo
	)

inlineoverridevirtual

This function is designed to calculate just the RHS contribution.

Parameters

rElement	The element to compute
RHS_Contribution	The RHS vector contribution
rEquationIdVector	The ID's of the element degrees of freedom
rCurrentProcessInfo	The current process info instance

Reimplemented from Kratos::Scheme< TSparseSpace, TDenseSpace >.

◆ CalculateUpdatedGridVelocityField()

template<class TSparseSpace , class TDenseSpace >

void Kratos::MPMExplicitScheme< TSparseSpace, TDenseSpace >::CalculateUpdatedGridVelocityField	(	const ProcessInfo &	rCurrentProcessInfo,
		ModelPart &	rModelPart
	)

inline

◆ Check()

template<class TSparseSpace , class TDenseSpace >

int Kratos::MPMExplicitScheme< TSparseSpace, TDenseSpace >::Check ( const ModelPart & rModelPart ) const

inlineoverridevirtual

This function is designed to be called once to perform all the checks needed on the input provided. Checks can be "expensive" as the function is designed to catch user's errors.

Parameters

r_model_part

Returns: 0 all ok

Reimplemented from Kratos::Scheme< TSparseSpace, TDenseSpace >.

◆ Clone()

template<class TSparseSpace , class TDenseSpace >

BaseTypePointer Kratos::MPMExplicitScheme< TSparseSpace, TDenseSpace >::Clone ( )

inlineoverridevirtual

Clone

Reimplemented from Kratos::Scheme< TSparseSpace, TDenseSpace >.

◆ FinalizeSolutionStep()

template<class TSparseSpace , class TDenseSpace >

void Kratos::MPMExplicitScheme< TSparseSpace, TDenseSpace >::FinalizeSolutionStep	(	ModelPart &	rModelPart,
		TSystemMatrixType &	A,
		TSystemVectorType &	Dx,
		TSystemVectorType &	b
	)

inlineoverridevirtual

Function called once at the end of a solution step, after convergence is reached if an iterative process is needed

Reimplemented from Kratos::Scheme< TSparseSpace, TDenseSpace >.

◆ GetDofList() [1/2]

template<class TSparseSpace , class TDenseSpace >

void Kratos::MPMExplicitScheme< TSparseSpace, TDenseSpace >::GetDofList	(	const Condition &	rCurrentCondition,
		Element::DofsVectorType &	rConditionDofList,
		const ProcessInfo &	rCurrentProcessInfo
	)

inlineoverridevirtual

Function that returns the list of Degrees of freedom to be assembled in the system for a Given Element

Reimplemented from Kratos::Scheme< TSparseSpace, TDenseSpace >.

◆ GetDofList() [2/2]

template<class TSparseSpace , class TDenseSpace >

void Kratos::MPMExplicitScheme< TSparseSpace, TDenseSpace >::GetDofList	(	const Element &	rCurrentElement,
		Element::DofsVectorType &	ElementalDofList,
		const ProcessInfo &	CurrentProcessInfo
	)

inlineoverridevirtual

Function that returns the list of Degrees of freedom to be assembled in the system for a Given Element

Reimplemented from Kratos::Scheme< TSparseSpace, TDenseSpace >.

◆ Initialize()

template<class TSparseSpace , class TDenseSpace >

void Kratos::MPMExplicitScheme< TSparseSpace, TDenseSpace >::Initialize ( ModelPart & rModelPart )

inlineoverridevirtual

This is the place to initialize the Scheme. This is intended to be called just once when the strategy is initialized.

Parameters

rModelPart The model of the problem to solve

Reimplemented from Kratos::Scheme< TSparseSpace, TDenseSpace >.

◆ InitializeSolutionStep()

template<class TSparseSpace , class TDenseSpace >

void Kratos::MPMExplicitScheme< TSparseSpace, TDenseSpace >::InitializeSolutionStep	(	ModelPart &	r_model_part,
		TSystemMatrixType &	A,
		TSystemVectorType &	Dx,
		TSystemVectorType &	b
	)

inlineoverridevirtual

initializes time step solution only for reasons if the time step solution is restarted

Parameters

r_model_part
A	LHS matrix
Dx	incremental update of primary variables
b	RHS Vector

Reimplemented from Kratos::Scheme< TSparseSpace, TDenseSpace >.

◆ KRATOS_CLASS_POINTER_DEFINITION()

template<class TSparseSpace , class TDenseSpace >

Kratos::MPMExplicitScheme< TSparseSpace, TDenseSpace >::KRATOS_CLASS_POINTER_DEFINITION ( MPMExplicitScheme< TSparseSpace, TDenseSpace > )

◆ Update()

template<class TSparseSpace , class TDenseSpace >

void Kratos::MPMExplicitScheme< TSparseSpace, TDenseSpace >::Update	(	ModelPart &	r_model_part,
		DofsArrayType &	rDofSet,
		TSystemMatrixType &	A,
		TSystemVectorType &	Dx,
		TSystemVectorType &	b
	)

inlineoverridevirtual

Performing the update of the solution Incremental update within newton iteration. It updates the state variables at the end of the time step: u_{n+1}^{k+1}= u_{n+1}^{k}+ \Delta u

Parameters

r_model_part
rDofSet	set of all primary variables
A	LHS matrix
Dx	incremental update of primary variables
b	RHS Vector

Working in 2D/3D (the definition of DOMAIN_SIZE is check in the Check method)

Reimplemented from Kratos::Scheme< TSparseSpace, TDenseSpace >.

◆ UpdateTranslationalDegreesOfFreedom()

template<class TSparseSpace , class TDenseSpace >

void Kratos::MPMExplicitScheme< TSparseSpace, TDenseSpace >::UpdateTranslationalDegreesOfFreedom	(	const ProcessInfo &	r_current_process_info,
		NodeIterator	itCurrentNode,
		const IndexType	DisplacementPosition,
		const double	delta_time,
		const SizeType	DomainSize = `3`
	)

inline

Member Data Documentation

◆ mr_grid_model_part

template<class TSparseSpace , class TDenseSpace >

ModelPart& Kratos::MPMExplicitScheme< TSparseSpace, TDenseSpace >::mr_grid_model_part

protected

◆ numerical_limit

template<class TSparseSpace , class TDenseSpace >

constexpr double Kratos::MPMExplicitScheme< TSparseSpace, TDenseSpace >::numerical_limit = std::numeric_limits<double>::epsilon()

staticconstexpr

The definition of the numerical limit.

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

/home/runner/work/Documentation/Documentation/master/applications/ParticleMechanicsApplication/custom_strategies/schemes/mpm_explicit_scheme.hpp

Protected Attributes

Type Definitions

Additional Inherited Members

Detailed Description

template<class TSparseSpace, class TDenseSpace> class Kratos::MPMExplicitScheme< TSparseSpace, TDenseSpace >

Member Typedef Documentation

◆ BaseType

◆ BaseTypePointer

◆ ConditionsArrayType

◆ DofsArrayType

◆ DofsVectorType

◆ ElementsArrayType

◆ IndexType

◆ LocalSystemMatrixType

◆ LocalSystemVectorType

◆ NodeIterator

◆ NodesArrayType

◆ SizeType

◆ TDataType

◆ TSystemMatrixType

◆ TSystemVectorType

Constructor & Destructor Documentation

◆ MPMExplicitScheme()

◆ ~MPMExplicitScheme()

Member Function Documentation

◆ calculateGridVelocityAndApplyDirichletBC()

◆ CalculateRHSContribution() [1/2]

◆ CalculateRHSContribution() [2/2]

◆ CalculateUpdatedGridVelocityField()

◆ Check()

◆ Clone()

◆ FinalizeSolutionStep()

◆ GetDofList() [1/2]

◆ GetDofList() [2/2]

◆ Initialize()

◆ InitializeSolutionStep()

◆ KRATOS_CLASS_POINTER_DEFINITION()

◆ Update()

◆ UpdateTranslationalDegreesOfFreedom()

Member Data Documentation

◆ mr_grid_model_part

◆ numerical_limit

template<class TSparseSpace, class TDenseSpace>
class Kratos::MPMExplicitScheme< TSparseSpace, TDenseSpace >