#include <CFiniteElementProblem.h>

Inheritance diagram for openbeam::CFiniteElementProblem:

Classes
struct	node_used_t

struct	TNodeElementConnection

struct	TProblemDOFIndicesForNode

Public Types
using	constraint_list_t = std::map< size_t, num_t >

using	load_list_t = std::map< size_t, num_t >

Public Member Functions
High-level management
virtual void	clear ()

bool	loadFromStream (std::istream &is, const mrpt::optional_ref< vector_string_t > &errMsg=std::nullopt, const mrpt::optional_ref< vector_string_t > &warnMsg=std::nullopt)

bool	loadFromFile (const std::string &file, const mrpt::optional_ref< vector_string_t > &errMsg=std::nullopt, const mrpt::optional_ref< vector_string_t > &warnMsg=std::nullopt)

bool	saveAsImageSVG (const std::string &file, const DrawStructureOptions &options, const StaticSolveProblemInfo solver_info=nullptr, const MeshOutputInfo meshing_info=nullptr, ImageSaveOutputInfo *out_img_info=nullptr) const

bool	saveAsImagePNG (const std::string &file, const DrawStructureOptions &options, const StaticSolveProblemInfo solver_info=nullptr, const MeshOutputInfo meshing_info=nullptr, ImageSaveOutputInfo *out_img_info=nullptr) const

bool	saveAsImage (const std::string &file, const bool is_svg, const DrawStructureOptions &options, const StaticSolveProblemInfo solver_info=nullptr, const MeshOutputInfo meshing_info=nullptr, ImageSaveOutputInfo *out_img_info=nullptr) const

bool	renderToCairoContext (void _cairo_context, const RenderInitData &ri, const DrawStructureOptions &options, const StaticSolveProblemInfo solver_info, const MeshOutputInfo *meshing_info) const

mrpt::opengl::CSetOfObjects::Ptr	getVisualization (const DrawStructureOptions &options, const StaticSolveProblemInfo &solver_info, const MeshOutputInfo meshing_info=nullptr, const StressInfo stressInfo=nullptr) const

Problem elements
size_t	insertElement (CElement::Ptr el)

template<typename ElementClass , typename... _Args>
size_t	createElement (_Args &&... __args)

CElement::ConstPtr	getElement (size_t i) const

const CElement::Ptr &	getElement (size_t i)

size_t	getNumberOfElements () const

Problem constraints
void	insertConstraint (const size_t dof_index, const num_t value=0)

const constraint_list_t &	getAllConstraints () const

void	setLoadAtDOF (const size_t dof_index, const num_t f)

void	addLoadAtDOF (const size_t dof_index, const num_t f)

const load_list_t &	getOverallLoadsOnDOFs () const

Problem nodes and their reference 3D pose (to this, one must add
the displacement obtained as solution to the problem).
void	setNumberOfNodes (size_t N)

size_t	getNumberOfNodes () const

node_index_t	insertNode (const TRotationTrans3D &p)

void	setNodePose (size_t idx, const TRotationTrans3D &p)

void	setNodePose (size_t idx, const num_t x, const num_t y, const num_t z)

TRotationTrans3D &	getNodePose (size_t i)

const TRotationTrans3D &	getNodePose (size_t i) const

void	getNodeDeformedPosition (size_t i, Vector3 &out_final_point, const StaticSolveProblemInfo &solver_info, const num_t exageration_factor=1) const

num_t	getMaximumDeformedDisplacement (const StaticSolveProblemInfo &solver_info) const

void	getBoundingBox (num_t &min_x, num_t &max_x, num_t &min_y, num_t &max_y, bool deformed=false, const StaticSolveProblemInfo *solver_info=nullptr, num_t deformed_scale_factor=1.0) const

Problem solving
virtual void	updateAll ()

void	assembleProblem (BuildProblemInfo &out_info)

void	solveStatic (StaticSolveProblemInfo &out_info, const StaticSolverOptions &opts=StaticSolverOptions())

void	postProcCalcStress (StressInfo &out_stress, const StaticSolveProblemInfo &solver_info)

std::string	getNodeLabel (const size_t idx) const
	"N%i" or custom label

Protected Types
using	TNodeConnections = std::map< element_index_t, TNodeElementConnection >

Protected Member Functions
virtual void	internalComputeStressAndEquivalentLoads ()

bool	internal_loadFromYaml (const mrpt::containers::yaml &f, const mrpt::optional_ref< vector_string_t > &errMsg, const mrpt::optional_ref< vector_string_t > &warnMsg)

void	internal_parser1_Parameters (const mrpt::containers::yaml &f, EvaluationContext &ctx) const

void	internal_parser2_BeamSections (const mrpt::containers::yaml &f, EvaluationContext &ctx) const

void	internal_parser3_nodes (const mrpt::containers::yaml &f, EvaluationContext &ctx)

void	internal_parser4_elements (const mrpt::containers::yaml &f, EvaluationContext &ctx)

void	internal_parser5_constraints (const mrpt::containers::yaml &f, EvaluationContext &ctx)

void	internal_parser6_node_loads (const mrpt::containers::yaml &f, EvaluationContext &ctx)

void	internal_parser7_element_loads (const mrpt::containers::yaml &f, EvaluationContext &ctx)

void	updateListDoFs ()

void	updateNodeConnections ()

void	updateElementsOrientation ()

void	updateNodesMainOrientation ()

void	internal_getVisualization_nodeLoads (mrpt::opengl::CSetOfObjects &gl, const DrawStructureOptions &options, const StaticSolveProblemInfo &solver_info, const MeshOutputInfo *meshing_info, num_t DEFORMED_SCALE_FACTOR) const

void	internal_getVisualization_constraints (mrpt::opengl::CSetOfObjects &gl, const DrawStructureOptions &options, const StaticSolveProblemInfo &solver_info, const MeshOutputInfo *meshing_info, num_t DEFORMED_SCALE_FACTOR) const

void	internal_getVisualization_distributedLoads (const CStructureProblem &str, mrpt::opengl::CSetOfObjects &gl, const DrawStructureOptions &options, const StaticSolveProblemInfo &solver_info, const MeshOutputInfo *meshing_info, num_t DEFORMED_SCALE_FACTOR) const

void	internal_getVisualization_stressDiagrams (mrpt::opengl::CSetOfObjects &gl, const DrawStructureOptions &options, const StaticSolveProblemInfo &solverInfo, const MeshOutputInfo *meshingInfo, num_t DEFORMED_SCALE_FACTOR, const StressInfo &stressInfo) const

Protected Attributes
std::vector< NodeDoF >	m_problem_DoFs

std::vector< TProblemDOFIndicesForNode >	m_problem_DoFs_inverse_list

std::vector< TNodeConnections >	m_node_connections

std::vector< TRotation3D >	m_nodeMainDirection

Main data
std::deque< node_used_t >	m_node_defined

std::deque< TRotationTrans3D >	m_node_poses

std::vector< std::string >	m_node_labels
	node custom label

std::deque< CElement::Ptr >	m_elements

constraint_list_t	m_DoF_constraints

load_list_t	m_loads_at_each_dof

load_list_t	m_loads_at_each_dof_equivs

std::map< size_t, ElementStress >	m_extra_stress_for_each_element

Problem DOFs (abstract list of DoF, doesn't have quantitative
displacement values)
const std::vector< NodeDoF > &	getProblemDoFs ()

std::string	getProblemDoFsDescription ()

size_t	getDOFIndex (const size_t nNode, const DoF_index n) const

static std::vector< size_t >	complementaryDoFs (const std::vector< size_t > &ds, const size_t nTotalDOFs)

Detailed Description

A complete problem of finite elements. It contains:

Elements: A list of finite elements from which the problem is built of.
Nodes: A list of 6 DoF nodes at certain space locations. Elements are always defined between 2 or more nodes.
Constraints: External loads and DoF movement restrictions. They are applied over the nodes.

See also: CFiniteElementProblem

Member Typedef Documentation

◆ TNodeConnections

using openbeam::CFiniteElementProblem::TNodeConnections = std::map<element_index_t, TNodeElementConnection>

protected

For each node, this map has the list of connected elements and the type of the connection.

Member Function Documentation

◆ addLoadAtDOF()

void openbeam::CFiniteElementProblem::addLoadAtDOF	(	const size_t	dof_index,
		const num_t	f
	)

Adds a load (force or moment) to a given DOF, accumulating to previously existing values if any.

See also: setLoadAtDOF

◆ assembleProblem()

void openbeam::CFiniteElementProblem::assembleProblem ( BuildProblemInfo & out_info )

For any partitioning of the DoFs into fixed (restricted, boundary conditions) and the rest (the free variables free_dof_indices), where DoF indices refer to m_problem_DoFs, computes three sparse matrices: K_{FF} (free-free), K_{BB} (boundary-boundary) and K_{BF} (boundary-free) which define the stiffness matrix of the problem.

Note: The free vs constrained DoF's are automatically determined from m_DoF_constraints

See also: solveStatic

◆ clear()

virtual void openbeam::CFiniteElementProblem::clear ( )

virtual

Delete all elements, nodes and constraints in this problem, completely emptying it.

Reimplemented in openbeam::CStructureProblem.

◆ complementaryDoFs()

static std::vector<size_t> openbeam::CFiniteElementProblem::complementaryDoFs	(	const std::vector< size_t > &	ds,
		const size_t	nTotalDOFs
	)

static

Return the complementary list of DOFs, such as "return \CUP ds = 1:N" and "return \CAP ds = \empty".

◆ createElement()

template<typename ElementClass , typename... _Args>

size_t openbeam::CFiniteElementProblem::createElement ( _Args &&... __args )

inline

Insert a new element in the problem

Returns: The element index

◆ getBoundingBox()

void openbeam::CFiniteElementProblem::getBoundingBox	(	num_t &	min_x,
		num_t &	max_x,
		num_t &	min_y,
		num_t &	max_y,
		bool	deformed = `false`,
		const StaticSolveProblemInfo *	solver_info = `nullptr`,
		num_t	deformed_scale_factor = `1.0`
	)		const

Computes the bounding box of all existing nodes

◆ getDOFIndex()

size_t openbeam::CFiniteElementProblem::getDOFIndex	(	const size_t	nNode,
		const DoF_index	n
	)		const

Returns the 0-based index of the n'th DOF (0-5) of node "nNode" in the list of DOFs to consider in the problem. Callable only after buildAll(). Returns string::npos if the DOF is NOT considered in the problem.

◆ getElement() [1/2]

const CElement::Ptr& openbeam::CFiniteElementProblem::getElement ( size_t i )

Return a pointer to the i'th element in the problem.

◆ getElement() [2/2]

CElement::ConstPtr openbeam::CFiniteElementProblem::getElement ( size_t i ) const

Return a constant pointer to the i'th element in the problem.

◆ getMaximumDeformedDisplacement()

num_t openbeam::CFiniteElementProblem::getMaximumDeformedDisplacement ( const StaticSolveProblemInfo & solver_info ) const

Returns the maximum absolute value translation in any X,Y,Z directions for a static problem solution

◆ getNodeDeformedPosition()

void openbeam::CFiniteElementProblem::getNodeDeformedPosition	(	size_t	i,
		Vector3 &	out_final_point,
		const StaticSolveProblemInfo &	solver_info,
		const num_t	exageration_factor = `1`
	)		const

Returns the 3D location of the final deformed state after solving the problem

◆ getNodePose() [1/2]

TRotationTrans3D& openbeam::CFiniteElementProblem::getNodePose ( size_t i )

inline

The 3D location of a node; its orientation (if different than global coordinates) serves as a local frame for constraints.

◆ getNodePose() [2/2]

const TRotationTrans3D& openbeam::CFiniteElementProblem::getNodePose ( size_t i ) const

inline

The 3D location of a node; its orientation (if different than global coordinates) serves as a local frame for constraints.

◆ getNumberOfElements()

size_t openbeam::CFiniteElementProblem::getNumberOfElements ( ) const

inline

Get the number of elements in the problem

◆ getNumberOfNodes()

size_t openbeam::CFiniteElementProblem::getNumberOfNodes ( ) const

inline

Number of DoF's?

◆ getProblemDoFs()

const std::vector<NodeDoF>& openbeam::CFiniteElementProblem::getProblemDoFs ( )

inline

Return the list of all DoF in the problem .

See also: getProblemDoFsDescription

◆ getProblemDoFsDescription()

std::string openbeam::CFiniteElementProblem::getProblemDoFsDescription ( )

Make a list with all DoF in the problem sorted by nodes.

See also: getProblemDoFs

◆ insertConstraint()

void openbeam::CFiniteElementProblem::insertConstraint	(	const size_t	dof_index,
		const num_t	value = `0`
	)

Adds a new displacement constraint to the problem, at the DoF dof_index (indices wrt list returned by getProblemDoFs() ).

◆ insertElement()

size_t openbeam::CFiniteElementProblem::insertElement ( CElement::Ptr el )

Insert a new element in the problem

Returns: The element index

◆ insertNode()

node_index_t openbeam::CFiniteElementProblem::insertNode ( const TRotationTrans3D & p )

Insert a new node to the problem with the given 6D pose, in global coordinates.

Returns: The index of this new node

◆ internal_loadFromYaml()

bool openbeam::CFiniteElementProblem::internal_loadFromYaml	(	const mrpt::containers::yaml &	f,
		const mrpt::optional_ref< vector_string_t > &	errMsg,
		const mrpt::optional_ref< vector_string_t > &	warnMsg
	)

protected

Internal common implementation of loadFrom*() methods()

◆ internalComputeStressAndEquivalentLoads()

virtual void openbeam::CFiniteElementProblem::internalComputeStressAndEquivalentLoads ( )

inlineprotectedvirtual

In base classes, process loads on elements and populate the m_loads_at_each_dof_equivs and m_extra_stress_for_each_element

◆ loadFromFile()

bool openbeam::CFiniteElementProblem::loadFromFile	(	const std::string &	file,
		const mrpt::optional_ref< vector_string_t > &	errMsg = `std::nullopt`,
		const mrpt::optional_ref< vector_string_t > &	warnMsg = `std::nullopt`
	)

Discard the current problem and loads it by parsing the given text file. The expected format of the file is described in XXXX.

Parameters

[in]	is	The input stream.
[out]	errMsg	The error messages found while parsing the text.
[out]	warnMsg	The warning messages found while parsing the text.

Returns: true on Success

◆ loadFromStream()

bool openbeam::CFiniteElementProblem::loadFromStream	(	std::istream &	is,
		const mrpt::optional_ref< vector_string_t > &	errMsg = `std::nullopt`,
		const mrpt::optional_ref< vector_string_t > &	warnMsg = `std::nullopt`
	)

Discard the current problem and loads it by parsing the given text input stream. The expected format of the file is described in XXXX.

Parameters

[in]	is	The input stream.
[out]	errMsg	The error messages found while parsing the text.
[out]	warnMsg	The warning messages found while parsing the text.

Returns: true on Success

◆ postProcCalcStress()

void openbeam::CFiniteElementProblem::postProcCalcStress	(	StressInfo &	out_stress,
		const StaticSolveProblemInfo &	solver_info
	)

After solving the problem with solveStatic, you can optionally call this post-processing method to evaluate the stress of all the elements.

See also: solveStatic

◆ saveAsImageSVG()

bool openbeam::CFiniteElementProblem::saveAsImageSVG	(	const std::string &	file,
		const DrawStructureOptions &	options,
		const StaticSolveProblemInfo *	solver_info = `nullptr`,
		const MeshOutputInfo *	meshing_info = `nullptr`,
		ImageSaveOutputInfo *	out_img_info = `nullptr`
	)		const

Saves a representation of the problem to a SVG file.

Parameters

options Many parameters and switches to control what will be drawn and hiden.

Returns: true on success

◆ setLoadAtDOF()

void openbeam::CFiniteElementProblem::setLoadAtDOF	(	const size_t	dof_index,
		const num_t	f
	)

Sets the load (force or moment) applied to a given DOF

See also: addLoadAtDOF

◆ setNodePose() [1/2]

void openbeam::CFiniteElementProblem::setNodePose	(	size_t	idx,
		const num_t	x,
		const num_t	y,
		const num_t	z
	)

Set the reference 3D pose of a node, in global coordinates

◆ setNodePose() [2/2]

void openbeam::CFiniteElementProblem::setNodePose	(	size_t	idx,
		const TRotationTrans3D &	p
	)

Set the reference 6D pose of a node, in global coordinates

◆ setNumberOfNodes()

void openbeam::CFiniteElementProblem::setNumberOfNodes ( size_t N )

Set the number of nodes in the problem

◆ solveStatic()

void openbeam::CFiniteElementProblem::solveStatic	(	StaticSolveProblemInfo &	out_info,
		const StaticSolverOptions &	opts = `StaticSolverOptions()`
	)

Solve the static FE problem, returning the resulting displacements and reaction forces.

Note: This method internally calls assembleProblem()

See also: assembleProblem, postProcCalcStress

◆ updateAll()

virtual void openbeam::CFiniteElementProblem::updateAll ( )

virtual

Update all internal lists after changing the problem. This actually calls:

Reimplemented in openbeam::CStructureProblem.

◆ updateElementsOrientation()

void openbeam::CFiniteElementProblem::updateElementsOrientation ( )

protected

Call the "updateOrientationFromNodePositions" in each element

◆ updateListDoFs()

void openbeam::CFiniteElementProblem::updateListDoFs ( )

protected

From the list of elements and their properties and connections, build the list of DoFs relevant to the problem.

◆ updateNodeConnections()

void openbeam::CFiniteElementProblem::updateNodeConnections ( )

protected

From the list of elements, update m_node_connections

◆ updateNodesMainOrientation()

void openbeam::CFiniteElementProblem::updateNodesMainOrientation ( )

protected

Fills up m_nodeMainDirection

Member Data Documentation

◆ m_DoF_constraints

constraint_list_t openbeam::CFiniteElementProblem::m_DoF_constraints

protected

List of constrainsts for each "fixed/constrained" DoF. Map key are indices in m_problem_DoFs. Map values are displacement in that DoF wrt the current node pose. Units are SI (meters or radians).

See also: getProblemDoFs

◆ m_extra_stress_for_each_element

std::map<size_t, ElementStress> openbeam::CFiniteElementProblem::m_extra_stress_for_each_element

protected

The extra stress to add up to each element as computed by internalComputeStressAndEquivalentLoads map: element index -> vector for each "port" -> vector of 6 stresses values.

◆ m_loads_at_each_dof

load_list_t openbeam::CFiniteElementProblem::m_loads_at_each_dof

protected

The vector of overall external loads (F_L) - Map keys are indices of m_problem_DoFs

◆ m_loads_at_each_dof_equivs

load_list_t openbeam::CFiniteElementProblem::m_loads_at_each_dof_equivs

protected

The vector of overall loads (F_L') on each DoF due to distributed forces - Map keys are indices of m_problem_DoFs Updated by assembleProblem() via internalComputeStressAndEquivalentLoads()

◆ m_node_connections

std::vector<TNodeConnections> openbeam::CFiniteElementProblem::m_node_connections

protected

A vector with an item for each node (list in m_node_poses) Each item is a map from element_id -> TNodeElementConnection Length=number of nodes.

◆ m_nodeMainDirection

std::vector<TRotation3D> openbeam::CFiniteElementProblem::m_nodeMainDirection

protected

For each node, the +X direction of the first element touching that node. Filled by updateNodesMainOrientation()

◆ m_problem_DoFs

std::vector<NodeDoF> openbeam::CFiniteElementProblem::m_problem_DoFs

protected

The list of DoFs variables of my problem. Built by updateListDoFs()

See also: m_problem_DoFs_inverse_list

◆ m_problem_DoFs_inverse_list

std::vector<TProblemDOFIndicesForNode> openbeam::CFiniteElementProblem::m_problem_DoFs_inverse_list

protected

Built together with m_problem_DoFs by updateListDoFs , holds the mapping from node index -> 6 numbers, with the index of that global DOF in the list of problem DOFs, or -1 if the DOF is not considered.

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

/home/jlblanco/code/openbeam/libopenbeam/include/openbeam/CFiniteElementProblem.h

Classes

Public Types

Public Member Functions

Protected Types

Protected Member Functions

Protected Attributes

Problem DOFs (abstract list of DoF, doesn't have quantitative

Detailed Description

Member Typedef Documentation

◆ TNodeConnections

Member Function Documentation

◆ addLoadAtDOF()

◆ assembleProblem()

◆ clear()

◆ complementaryDoFs()

◆ createElement()

◆ getBoundingBox()

◆ getDOFIndex()

◆ getElement() [1/2]

◆ getElement() [2/2]

◆ getMaximumDeformedDisplacement()

◆ getNodeDeformedPosition()

◆ getNodePose() [1/2]

◆ getNodePose() [2/2]

◆ getNumberOfElements()

◆ getNumberOfNodes()

◆ getProblemDoFs()

◆ getProblemDoFsDescription()

◆ insertConstraint()

◆ insertElement()

◆ insertNode()

◆ internal_loadFromYaml()

◆ internalComputeStressAndEquivalentLoads()

◆ loadFromFile()

◆ loadFromStream()

◆ postProcCalcStress()

◆ saveAsImageSVG()

◆ setLoadAtDOF()

◆ setNodePose() [1/2]

◆ setNodePose() [2/2]

◆ setNumberOfNodes()

◆ solveStatic()

◆ updateAll()

◆ updateElementsOrientation()

◆ updateListDoFs()

◆ updateNodeConnections()

◆ updateNodesMainOrientation()

Member Data Documentation

◆ m_DoF_constraints

◆ m_extra_stress_for_each_element

◆ m_loads_at_each_dof

◆ m_loads_at_each_dof_equivs

◆ m_node_connections

◆ m_nodeMainDirection

◆ m_problem_DoFs

◆ m_problem_DoFs_inverse_list