A LeeNewmark class defines a solver using Newmark algorithm with Lee damping model. More...

#include <LeeNewmark.h>

Inheritance diagram for LeeNewmark:

Collaboration diagram for LeeNewmark:

Public Member Functions
	LeeNewmark (unsigned, vec &&, vec &&, double, double)

int	initialize () override

int	process_constraint () override

int	process_constraint_resistance () override

void	assemble_resistance () override

void	print () override

Public Member Functions inherited from LeeNewmarkBase
	LeeNewmarkBase (unsigned, double, double, StiffnessType=StiffnessType::CURRENT)

int	initialize () override

int	update_internal (const mat &) final

int	solve (mat &, const mat &) final

int	solve (mat &, const sp_mat &) final

int	solve (mat &, mat &&) final

int	solve (mat &, sp_mat &&) final

vec	get_force_residual () final

vec	get_displacement_residual () final

void	commit_status () final

void	clear_status () final

void	reset_status () final

Public Member Functions inherited from Newmark
	Newmark (unsigned=0, double=.25, double=.5)

void	assemble_resistance () override

void	assemble_matrix () override

int	update_trial_status () override

void	update_parameter (double) override

vec	from_incre_velocity (const vec &, const uvec &) override

vec	from_incre_acceleration (const vec &, const uvec &) override

void	print () override

Public Member Functions inherited from ImplicitIntegrator
constexpr IntegratorType	type () const override

bool	time_independent_matrix () const override

	Integrator (unsigned=0)

	Integrator (const Integrator &)=delete

	Integrator (Integrator &&)=delete

Public Member Functions inherited from Integrator
	Integrator (unsigned=0)

	Integrator (const Integrator &)=delete

	Integrator (Integrator &&)=delete

Integrator &	operator= (const Integrator &)=delete

Integrator &	operator= (Integrator &&)=delete

	~Integrator () override=default

void	set_domain (const weak_ptr< DomainBase > &)

shared_ptr< DomainBase >	get_domain () const

virtual int	initialize ()

virtual constexpr IntegratorType	type () const

void	set_time_step_switch (bool)

bool	allow_to_change_time_step () const

void	set_matrix_assembled_switch (bool)

bool	matrix_is_assembled () const

virtual bool	has_corrector () const

virtual bool	time_independent_matrix () const

virtual int	process_load ()

virtual int	process_constraint ()

virtual int	process_criterion ()

virtual int	process_modifier ()

virtual int	process_load_resistance ()

virtual int	process_constraint_resistance ()

void	record () const

virtual void	assemble_resistance ()

virtual void	assemble_matrix ()

virtual vec	get_force_residual ()

virtual vec	get_displacement_residual ()

virtual vec	get_auxiliary_residual ()

virtual sp_mat	get_reference_load ()

virtual const vec &	get_trial_displacement () const

virtual void	update_load ()

virtual void	update_constraint ()

virtual void	update_trial_load_factor (double)

virtual void	update_trial_load_factor (const vec &)

virtual void	update_from_ninja ()

virtual void	update_trial_time (double)

virtual void	update_incre_time (double)

virtual int	update_trial_status ()

virtual int	correct_trial_status ()

virtual int	sync_status (bool)

virtual int	update_internal (const mat &)

mat	solve (const mat &)

mat	solve (const sp_mat &)

mat	solve (mat &&)

mat	solve (sp_mat &&)

virtual int	solve (mat &, const mat &)

virtual int	solve (mat &, const sp_mat &)

virtual int	solve (mat &, mat &&)

virtual int	solve (mat &, sp_mat &&)

virtual void	erase_machine_error (vec &) const

void	stage_and_commit_status ()

virtual void	stage_status ()

virtual void	commit_status ()

virtual void	clear_status ()

virtual void	reset_status ()

virtual void	update_parameter (double)

virtual vec	from_incre_velocity (const vec &, const uvec &)

virtual vec	from_incre_acceleration (const vec &, const uvec &)

virtual vec	from_total_velocity (const vec &, const uvec &)

virtual vec	from_total_acceleration (const vec &, const uvec &)

vec	from_incre_velocity (double, const uvec &)

vec	from_incre_acceleration (double, const uvec &)

vec	from_total_velocity (double, const uvec &)

vec	from_total_acceleration (double, const uvec &)

Public Member Functions inherited from Tag
	Tag (unsigned=0)

	Tag (const Tag &)=default

	Tag (Tag &&)=default

Tag &	operator= (const Tag &)=delete

Tag &	operator= (Tag &&)=delete

virtual	~Tag ()=default

void	set_tag (unsigned) const

unsigned	get_tag () const

void	enable ()

void	disable ()

void	guard ()

void	unguard ()

bool	is_active () const

bool	is_guarded () const

virtual void	print ()

Protected Attributes
shared_ptr< MetaMat< double > >	current_mass = nullptr

shared_ptr< MetaMat< double > >	current_stiffness = nullptr

shared_ptr< MetaMat< double > >	current_geometry = nullptr

Protected Attributes inherited from LeeNewmarkBase
const uword	n_block

const StiffnessType	stiffness_type

bool	first_iteration = true

bool	if_iterative = false

vec	current_internal

vec	trial_internal

vec	residual

unique_ptr< MetaMat< double > >	stiffness = nullptr

shared_ptr< Factory< double > >	factory = nullptr

Protected Attributes inherited from Newmark
double	C0 = 0.

double	C1 = 0.

double	C2 = 0.

double	C3 = 0.

double	C4 = 0.

double	C5 = 0.

Additional Inherited Members
Public Types inherited from LeeNewmarkBase
enum class	StiffnessType { INITIAL , CURRENT , TRIAL }

Protected Member Functions inherited from LeeNewmarkBase
virtual uword	get_total_size () const =0

virtual void	update_stiffness () const =0

virtual void	update_residual () const =0

int	erase_top_left_block () const

Detailed Description

A LeeNewmark class defines a solver using Newmark algorithm with Lee damping model.

Remarks:

Both dense and sparse storage schemes for original matrices are considered. The matrix–vector product may potentially be faster with dense formulation.

Author: tlc

Date: 25/05/2020

Version: 0.1.1

Constructor & Destructor Documentation

◆ LeeNewmark()

LeeNewmark::LeeNewmark	(	unsigned	T,
		vec &&	X,
		vec &&	F,
		double	A,
		double	B
	)

Member Function Documentation

◆ assemble_resistance()

void LeeNewmark::assemble_resistance ( )

overridevirtual

Reimplemented from Integrator.

Here is the call graph for this function:

◆ initialize()

int LeeNewmark::initialize ( )

overridevirtual

Reimplemented from Integrator.

Here is the call graph for this function:

◆ print()

void LeeNewmark::print ( )

overridevirtual

Reimplemented from Tag.

◆ process_constraint()

int LeeNewmark::process_constraint ( )

overridevirtual

The main task of this method is to apply constraints (of various forms implemented in various methods). Combinations of different types need to be considered: 1) homogeneous, 2) inhomogeneous, 3) linear, 4) nonlinear. Combinations of different methods need to be considered: 1) penalty, 2) multiplier. On exit, the global stiffness matrix should be updated, the global residual vector should be updated.

Reimplemented from Integrator.

Here is the call graph for this function:

◆ process_constraint_resistance()

int LeeNewmark::process_constraint_resistance ( )

overridevirtual

This method is similar to process_constraint(), but it only updates the global residual vector. The global stiffness matrix is not touched as in some solving schemes, the global stiffness matrix is only assembled and factorised once at the beginning. Subsequent iterations do not assemble the global stiffness matrix again and reuse the factorised matrix. In this case, the factorised matrix cannot be modified.

Reimplemented from Integrator.

Here is the call graph for this function:

Member Data Documentation

◆ current_geometry

shared_ptr<MetaMat<double> > LeeNewmark::current_geometry = nullptr

protected

◆ current_mass

shared_ptr<MetaMat<double> > LeeNewmark::current_mass = nullptr

protected

◆ current_stiffness

shared_ptr<MetaMat<double> > LeeNewmark::current_stiffness = nullptr

protected

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

Solver/Integrator/LeeNewmark.h
Solver/Integrator/LeeNewmark.cpp

Public Member Functions

Protected Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ LeeNewmark()

Member Function Documentation

◆ assemble_resistance()

◆ initialize()

◆ print()

◆ process_constraint()

◆ process_constraint_resistance()

Member Data Documentation

◆ current_geometry

◆ current_mass

◆ current_stiffness