Class for creating an uniaxial 3D two-node link element using Bouc-Wen formulation. More...

#include <UnxBoucWen3DLink.hpp>

Inheritance diagram for UnxBoucWen3DLink:

Collaboration diagram for UnxBoucWen3DLink:

Public Member Functions
	UnxBoucWen3DLink (const std::vector< unsigned int > nodes, std::vector< double > parameters, std::vector< double > variables, const unsigned int dim, const unsigned int dir, double tol=1E-06, unsigned int nmax=50)
	Creates a UnxBoucWen3DLink in a finite element Mesh. More...

	~UnxBoucWen3DLink ()
	Destroys this UnxBoucWen3DLink object. More...

void	CommitState ()
	Save the material states in the element. More...

void	ReverseState ()
	Reverse the material/section states to previous converged state in this element. More...

void	InitialState ()
	Brings the material/section state to its initial state in this element. More...

void	UpdateState ()
	Update the material states in the element. More...

void	SetDomain (std::map< unsigned int, std::shared_ptr< Node > > &nodes)
	Sets the finite element dependance among objects. More...

void	SetDamping (const std::shared_ptr< Damping > &damping)
	Sets the damping model. More...

std::vector< unsigned int >	GetTotalDegreeOfFreedom () const
	Gets the list of total-degree of freedom of this Element. More...

Eigen::MatrixXd	GetStrain () const
	Gets the material/section (generalised) strain. More...

Eigen::MatrixXd	GetStress () const
	Gets the material/section (generalised) stress. More...

Eigen::MatrixXd	GetStrainRate () const
	Gets the material/section (generalised) strain-rate. More...

Eigen::MatrixXd	GetStrainAt (double x3, double x2) const
	Gets the material strain in section at coordinate (x3,x2). More...

Eigen::MatrixXd	GetStressAt (double x3, double x2) const
	Gets the material stress in section at coordinate (x3,x2). More...

Eigen::VectorXd	GetVTKResponse (std::string response) const
	Gets the element internal response in VTK format for Paraview display. More...

double	ComputeEnergy ()
	Computes the element energy for a given deformation. More...

Eigen::MatrixXd	ComputeMassMatrix ()
	Compute the lumped/consistent mass matrix of the element. More...

Eigen::MatrixXd	ComputeStiffnessMatrix ()
	Compute the stiffness matrix of the element using gauss-integration. More...

Eigen::MatrixXd	ComputeDampingMatrix ()
	Compute the damping matrix of the element using gauss-integration. More...

Eigen::MatrixXd	ComputePMLMatrix ()
	Compute the PML history matrix using gauss-integration. More...

Eigen::VectorXd	ComputeInternalForces ()
	Compute the internal (elastic) forces acting on the element. More...

Eigen::VectorXd	ComputeInternalDynamicForces ()
	Compute the elastic, inertial, and viscous forces acting on the element. More...

Eigen::VectorXd	ComputeSurfaceForces (const std::shared_ptr< Load > &surface, unsigned int face)
	Compute the surface forces acting on the element. More...

Eigen::VectorXd	ComputeBodyForces (const std::shared_ptr< Load > &body, unsigned int k=0)
	Compute the body forces acting on the element. More...

Eigen::VectorXd	ComputeDomainReductionForces (const std::shared_ptr< Load > &drm, unsigned int k)
	Compute the domain reduction forces acting on the element. More...

Public Member Functions inherited from Element
	Element (std::string name, const std::vector< unsigned int > nodes, unsigned int ndofs, unsigned int VTKcell, unsigned int SVLcell)
	Creates an Element in a finite element Mesh. More...

virtual	~Element ()=0
	Destroys this Element object. More...

std::string	GetName () const
	Gets the Element Name. More...

unsigned int	GetVTKCellType () const
	Gets the Element VTK cell type. More...

unsigned int	GetVTKGroupType () const
	Gets the Element VTK group type. More...

unsigned int	GetNumberOfNodes () const
	Returns the number of nodes in element. More...

unsigned int	GetNumberOfDegreeOfFreedom () const
	Returns total number of degree of freedom in the element. More...

const std::vector< unsigned int > &	GetNodes () const
	Returns the Node Connectivity Indexes. More...

bool	HasFixedNode (const std::vector< std::shared_ptr< Node > > &nodes) const
	Returns if the element has fixed nodes. More...

Private Member Functions
double	sign (double x) const
	Sign function. More...

Eigen::VectorXd	ComputeRelativeDeformation () const
	Update strain in the element: More...

Eigen::MatrixXd	ComputeLocalAxes () const
	Compute local axes. More...

Eigen::MatrixXd	ComputeRotationMatrix () const
	Compute rotation matrix. More...

Private Attributes
double	Ko
	Initial stiffness of hysteretic component. More...

double	Fy
	Yield force. More...

double	alpha
	Yielding exponent (sharpness of hysteresis loop corners) More...

double	eta
	Yielding exponent (sharpness of hysteresis loop corners) More...

double	beta
	First hysteretic shape parameter. More...

double	gamma
	second hysteretic shape parameter More...

double	Tol
	Newton-Raphson Tolerance. More...

double	z
	Trial Hysteretic evolution parameters. More...

double	zn
	Hysteretic evolution parameters. More...

double	U
	Trial displacements in local coordinates. More...

double	Un
	Displacements in local coordinates. More...

double	qbw
	Non-linear Bouc-Wen internal force. More...

double	qbc
	Non-linear commited Bouc-Wen internal force. More...

double	kbw
	Consistent Bouc-Wen stiffness. More...

double	kbc
	Consistent commited Bouc-Wen stiffness. More...

unsigned int	nMax
	Newton-Raphson maximum number of iterations. More...

unsigned int	Dimension
	The element dimension in global coordinates. More...

unsigned int	Direction
	The element direction in global coordinates. More...

std::vector< std::shared_ptr< Node > >	theNodes
	The Element's Nodes. More...

Detailed Description

Class for creating an uniaxial 3D two-node link element using Bouc-Wen formulation.

See also: Element.hpp Mesh.hpp

Constructor & Destructor Documentation

◆ UnxBoucWen3DLink()

UnxBoucWen3DLink::UnxBoucWen3DLink	(	const std::vector< unsigned int >	nodes,
		std::vector< double >	parameters,
		std::vector< double >	variables,
		const unsigned int	dim,
		const unsigned int	dir,
		double	tol = `1E-06`,
		unsigned int	nmax = `50`
	)

Creates a UnxBoucWen3DLink in a finite element Mesh.

Parameters

nodes	The Node connectivity array of this Element.
parameters	Vector that contains the Bouc-Wen model parameters.
variables	Vector that contains auxiliary model parameters.
dim	The model dimension.
dir	The local direction where this Element is acting.
tol	Tolerance to stop Newton-Raphson iterations.
nmax	Maximum number of iteration for Newton-Raphson.

Note: More details can be found at UnxBoucWen3DLink.

See also: UnxBoucWen3DLink::theNodes, UnxBoucWen3DLink::theDirection, UnxBoucWen3DLink::theDimension.

◆ ~UnxBoucWen3DLink()

UnxBoucWen3DLink::~UnxBoucWen3DLink ( )

Destroys this UnxBoucWen3DLink object.

Member Function Documentation

◆ CommitState()

void UnxBoucWen3DLink::CommitState ( )

virtual

Save the material states in the element.

Note: This function sets the trial states as converged ones in Material.

Implements Element.

◆ ComputeBodyForces()

Eigen::VectorXd UnxBoucWen3DLink::ComputeBodyForces	(	const std::shared_ptr< Load > &	body,
		unsigned int	k = `0`
	)

virtual

Compute the body forces acting on the element.

Parameters

body	Pointer to the Load object that contains this information.
k	The time step at which the body load is evaluated.

Returns: Vector with the Element surface force.

Note: The body force vector can be revisited in UnxBoucWen3DLink.

See also: Assembler::ComputeExternalForceVector(), Integrator::ComputeEffectiveForce().

Implements Element.

◆ ComputeDampingMatrix()

Eigen::MatrixXd UnxBoucWen3DLink::ComputeDampingMatrix ( )

virtual

Compute the damping matrix of the element using gauss-integration.

Returns: Matrix with the Element damping matrix.

Note: The damping matrix can be revisited in UnxBoucWen3DLink.

See also: Assembler::ComputeDampingMatrix(), Integrator::ComputeEffectiveStiffness().

Implements Element.

◆ ComputeDomainReductionForces()

Eigen::VectorXd UnxBoucWen3DLink::ComputeDomainReductionForces	(	const std::shared_ptr< Load > &	drm,
		unsigned int	k
	)

virtual

Compute the domain reduction forces acting on the element.

Parameters

drm	Pointer to the DRM Load object that contains this information.
k	The time step at which the body load is evaluated.

Returns: Vector with the Element domain reduction forces.

Note: The DRM force vector can be revisited in UnxBoucWen3DLink.

See also: Assembler::ComputeExternalForceVector(), Integrator::ComputeEffectiveForce().

Implements Element.

◆ ComputeEnergy()

double UnxBoucWen3DLink::ComputeEnergy ( )

virtual

Computes the element energy for a given deformation.

Returns: Scalar with the element deformation energy.

Implements Element.

◆ ComputeInternalDynamicForces()

Eigen::VectorXd UnxBoucWen3DLink::ComputeInternalDynamicForces ( )

virtual

Compute the elastic, inertial, and viscous forces acting on the element.

Returns: Vector with the Element dynamic internal force.

Note: The internal force vector can be revisited in Element.

See also: Assembler::ComputeDynamicInternalForceVector().

Implements Element.

◆ ComputeInternalForces()

Eigen::VectorXd UnxBoucWen3DLink::ComputeInternalForces ( )

virtual

Compute the internal (elastic) forces acting on the element.

Returns: Vector with the Element internal force.

Note: The internal force vector can be revisited in UnxBoucWen3DLink.

See also: Assembler::ComputeInternalForceVector(), Integrator::ComputeEffectiveForce().

Implements Element.

◆ ComputeLocalAxes()

Eigen::MatrixXd UnxBoucWen3DLink::ComputeLocalAxes ( ) const

private

Compute local axes.

Returns: Matrix with the local axis transformation.

Note: Axes transformation are according to UnxBoucWen3DLink.

◆ ComputeMassMatrix()

Eigen::MatrixXd UnxBoucWen3DLink::ComputeMassMatrix ( )

virtual

Compute the lumped/consistent mass matrix of the element.

Returns: Matrix with the Element mass matrix.

Note: The mass matrix can be revisited in UnxBoucWen3DLink.

See also: Assembler::ComputeMassMatrix(), Integrator::ComputeEffectiveStiffness().

Implements Element.

◆ ComputePMLMatrix()

Eigen::MatrixXd UnxBoucWen3DLink::ComputePMLMatrix ( )

virtual

Compute the PML history matrix using gauss-integration.

Returns: Matrix with the PML Element matrix.

Note: The PML matrix is none existent for this element.

See also: Assembler::ComputePMLHistoryMatrix(), Integrator::ComputeEffectiveStiffness().

Implements Element.

◆ ComputeRelativeDeformation()

Eigen::VectorXd UnxBoucWen3DLink::ComputeRelativeDeformation ( ) const

private

Update strain in the element:

Returns: Vector with the total relative deformation.

◆ ComputeRotationMatrix()

Eigen::MatrixXd UnxBoucWen3DLink::ComputeRotationMatrix ( ) const

private

Compute rotation matrix.

Returns: Matrix with the rotation.

Note: Axes transformation are according to UnxBoucWen3DLink.

◆ ComputeStiffnessMatrix()

Eigen::MatrixXd UnxBoucWen3DLink::ComputeStiffnessMatrix ( )

virtual

Compute the stiffness matrix of the element using gauss-integration.

Returns: Matrix with the Element stiffness matrix.

Note: The stiffness matrix can be revisited in UnxBoucWen3DLink.

See also: Assembler::ComputeStiffnessMatrix(), Integrator::ComputeEffectiveStiffness().

Implements Element.

◆ ComputeSurfaceForces()

Eigen::VectorXd UnxBoucWen3DLink::ComputeSurfaceForces	(	const std::shared_ptr< Load > &	surface,
		unsigned int	face
	)

virtual

Compute the surface forces acting on the element.

Parameters

surface	Pointer to the Load object that contains this information.
k	The time step at which the surface load is evaluated.

Returns: Vector with the Element surface force.

Note: The surface force vector can be revisited in UnxBoucWen3DLink.

See also: Assembler::ComputeExternalForceVector(), Integrator::ComputeEffectiveForce().

Implements Element.

◆ GetStrain()

Eigen::MatrixXd UnxBoucWen3DLink::GetStrain ( ) const

virtual

Gets the material/section (generalised) strain.

Returns: Matrix with the strain at each integration point.

Note: The index (i,j) are the strain and Gauss-point respectively.

Implements Element.

◆ GetStrainAt()

Eigen::MatrixXd UnxBoucWen3DLink::GetStrainAt	(	double	x3,
		double	x2
	)		const

virtual

Gets the material strain in section at coordinate (x3,x2).

Parameters

x3	Local coordinate along the x3-axis.
x2	Local coordinate along the x2-axis.

Returns: Matrix with the strain at coordinate (x3,x2).

Note: The strains are interpolated at this coordinate.

Implements Element.

◆ GetStrainRate()

Eigen::MatrixXd UnxBoucWen3DLink::GetStrainRate ( ) const

virtual

Gets the material/section (generalised) strain-rate.

Returns: Matrix with the strain-rate at each integration point.

Note: The index (i,j) are the strain-rate and Gauss-point respectively.

Implements Element.

◆ GetStress()

Eigen::MatrixXd UnxBoucWen3DLink::GetStress ( ) const

virtual

Gets the material/section (generalised) stress.

Returns: Matrix with the stress at each integration point.

Note: The index (i,j) are the stress and Gauss-point respectively.

Implements Element.

◆ GetStressAt()

Eigen::MatrixXd UnxBoucWen3DLink::GetStressAt	(	double	x3,
		double	x2
	)		const

virtual

Gets the material stress in section at coordinate (x3,x2).

Parameters

x3	Local coordinate along the x3-axis.
x2	Local coordinate along the x2-axis.

Returns: Matrix with the stresses at coordinate (x3,x2).

Note: The stresses are interpolated at this coordinate.

Implements Element.

◆ GetTotalDegreeOfFreedom()

std::vector<unsigned int> UnxBoucWen3DLink::GetTotalDegreeOfFreedom ( ) const

virtual

Gets the list of total-degree of freedom of this Element.

Returns: Vector with the list of degree-of-freedom of this Element.

Implements Element.

◆ GetVTKResponse()

Eigen::VectorXd UnxBoucWen3DLink::GetVTKResponse ( std::string response ) const

virtual

Gets the element internal response in VTK format for Paraview display.

Parameters

response The response to be display in Paraview.

Returns: Vector with the response at the Element center.

Note: The current responses are: "Strain", "Stress".

Implements Element.

◆ InitialState()

void UnxBoucWen3DLink::InitialState ( )

virtual

Brings the material/section state to its initial state in this element.

Note: This function returns the meterial states to the beginning.

Implements Element.

◆ ReverseState()

void UnxBoucWen3DLink::ReverseState ( )

virtual

Reverse the material/section states to previous converged state in this element.

Note: This function returns the trial states to previous converged states at the Material level.

Implements Element.

◆ SetDamping()

void UnxBoucWen3DLink::SetDamping ( const std::shared_ptr< Damping > & damping )

virtual

Sets the damping model.

Parameters

damping Pointer to the damping model.

Note: Several Element objects can share the same damping model.

Implements Element.

◆ SetDomain()

void UnxBoucWen3DLink::SetDomain ( std::map< unsigned int, std::shared_ptr< Node > > & nodes )

virtual

Sets the finite element dependance among objects.

Parameters

nodes The Node list of the Mesh object.

Note: This function sets the relation between Node and Element objects.

See also: lin2DQuad4::theNodes.

Implements Element.

◆ sign()

double UnxBoucWen3DLink::sign ( double x ) const

private

Sign function.

◆ UpdateState()

void UnxBoucWen3DLink::UpdateState ( )

virtual

Update the material states in the element.

Note: This function update the trial states at the Material level.

Implements Element.

Member Data Documentation

◆ alpha

double UnxBoucWen3DLink::alpha

private

Yielding exponent (sharpness of hysteresis loop corners)

◆ beta

double UnxBoucWen3DLink::beta

private

First hysteretic shape parameter.

◆ Dimension

unsigned int UnxBoucWen3DLink::Dimension

private

The element dimension in global coordinates.

◆ Direction

unsigned int UnxBoucWen3DLink::Direction

private

The element direction in global coordinates.

◆ eta

double UnxBoucWen3DLink::eta

private

Yielding exponent (sharpness of hysteresis loop corners)

◆ Fy

double UnxBoucWen3DLink::Fy

private

Yield force.

◆ gamma

double UnxBoucWen3DLink::gamma

private

second hysteretic shape parameter

◆ kbc

double UnxBoucWen3DLink::kbc

private

Consistent commited Bouc-Wen stiffness.

◆ kbw

double UnxBoucWen3DLink::kbw

private

Consistent Bouc-Wen stiffness.

◆ Ko

double UnxBoucWen3DLink::Ko

private

Initial stiffness of hysteretic component.

◆ nMax

unsigned int UnxBoucWen3DLink::nMax

private

Newton-Raphson maximum number of iterations.

◆ qbc

double UnxBoucWen3DLink::qbc

private

Non-linear commited Bouc-Wen internal force.

◆ qbw

double UnxBoucWen3DLink::qbw

private

Non-linear Bouc-Wen internal force.

◆ theNodes

std::vector<std::shared_ptr<Node> > UnxBoucWen3DLink::theNodes

private

The Element's Nodes.

◆ Tol

double UnxBoucWen3DLink::Tol

private

Newton-Raphson Tolerance.

◆ U

double UnxBoucWen3DLink::U

private

Trial displacements in local coordinates.

◆ Un

double UnxBoucWen3DLink::Un

private

Displacements in local coordinates.

◆ z

double UnxBoucWen3DLink::z

private

Trial Hysteretic evolution parameters.

◆ zn

double UnxBoucWen3DLink::zn

private

Hysteretic evolution parameters.

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ UnxBoucWen3DLink()

◆ ~UnxBoucWen3DLink()

Member Function Documentation

◆ CommitState()

◆ ComputeBodyForces()

◆ ComputeDampingMatrix()

◆ ComputeDomainReductionForces()

◆ ComputeEnergy()

◆ ComputeInternalDynamicForces()

◆ ComputeInternalForces()

◆ ComputeLocalAxes()

◆ ComputeMassMatrix()

◆ ComputePMLMatrix()

◆ ComputeRelativeDeformation()

◆ ComputeRotationMatrix()

◆ ComputeStiffnessMatrix()

◆ ComputeSurfaceForces()

◆ GetStrain()

◆ GetStrainAt()

◆ GetStrainRate()

◆ GetStress()

◆ GetStressAt()

◆ GetTotalDegreeOfFreedom()

◆ GetVTKResponse()

◆ InitialState()

◆ ReverseState()

◆ SetDamping()

◆ SetDomain()

◆ sign()

◆ UpdateState()

Member Data Documentation

◆ alpha

◆ beta

◆ Dimension

◆ Direction

◆ eta

◆ Fy

◆ gamma

◆ kbc

◆ kbw

◆ Ko

◆ nMax

◆ qbc

◆ qbw

◆ theNodes

◆ Tol

◆ U

◆ Un

◆ z

◆ zn