#include "b2element_field_transfer.H"

Inheritance diagram for b2000::ElementCFD3DField< SHAPE, EULER >:

[legend]

Collaboration diagram for b2000::ElementCFD3DField< SHAPE, EULER >:

[legend]

Public Member Functions
std::pair< int, Node const >	get_nodes () const override

void	set_nodes (std::pair< int, Node const > nodes_) override

int	face_field_order (const int face_id, const std::string &field_name) override

bool	face_field_linear_on_dof (const int edge_id, const std::string &field_name) override

void	face_geom (const int face_id, const b2linalg::Vector< double, b2linalg::Vdense_constref > &internal_coor, b2linalg::Vector< double > &geom, b2linalg::Matrix< double, b2linalg::Mrectangle > &d_geom_d_icoor) override

void	set_property (ElementProperty *property_) override

const ElementProperty *	get_property () const override

Public Member Functions inherited from b2000::TypedElement< double >
virtual void	get_constraint (Model &model, const bool linear, const double time, const b2linalg::Matrix< double, b2linalg::Mrectangle_constref > &dof, b2linalg::Index &dof_numbering, b2linalg::Vector< double, b2linalg::Vdense > &constraint, b2linalg::Matrix< double, b2linalg::Mcompressed_col > &trans_d_constraint_d_dof, b2linalg::Vector< double, b2linalg::Vdense > &d_constraint_d_time)

virtual void	edge_field_value (const int edge_id, const std::string &field_name, const double internal_coor, const b2linalg::Matrix< double, b2linalg::Mrectangle_constref > &dof, const double time, b2linalg::Vector< double, b2linalg::Vdense > &value, b2linalg::Vector< double, b2linalg::Vdense > &d_value_d_icoor, b2linalg::Index &dof_numbering, b2linalg::Matrix< double, b2linalg::Mrectangle > &d_value_d_dof, b2linalg::Index &d_value_d_dof_dep_col=b2linalg::Index::null)

virtual void	face_field_value (const int face_id, const std::string &field_name, const b2linalg::Vector< double, b2linalg::Vdense_constref > &internal_coor, const b2linalg::Matrix< double, b2linalg::Mrectangle_constref > &dof, const double time, b2linalg::Vector< double, b2linalg::Vdense > &value, b2linalg::Matrix< double, b2linalg::Mrectangle > &d_value_d_icoor, b2linalg::Index &dof_numbering, b2linalg::Matrix< double, b2linalg::Mrectangle > &d_value_d_dof, b2linalg::Index &d_value_d_dof_dep_col)

virtual void	body_field_value (const std::string &field_name, const b2linalg::Vector< double, b2linalg::Vdense_constref > &internal_coor, const b2linalg::Matrix< double, b2linalg::Mrectangle_constref > &dof, const double time, b2linalg::Vector< double, b2linalg::Vdense > &value, b2linalg::Matrix< double, b2linalg::Mrectangle > &d_value_d_icoor, b2linalg::Index &dof_numbering, b2linalg::Matrix< double, b2linalg::Mrectangle > &d_value_d_dof)

virtual void	field_volume_integration (const b2linalg::Vector< double, b2linalg::Vdense_constref > &dof, const b2linalg::Vector< double, b2linalg::Vdense_constref > &dofdot, const b2linalg::Vector< double, b2linalg::Vdense_constref > &dofdotdot, const Field< double > &f, b2linalg::Index &dof_numbering, b2linalg::Vector< double, b2linalg::Vdense > &discretised_field, b2linalg::Matrix< double, b2linalg::Mpacked > &d_discretised_field_d_dof)

virtual void	field_face_integration (const b2linalg::Vector< double, b2linalg::Vdense_constref > &dof, const b2linalg::Vector< double, b2linalg::Vdense_constref > &dofdot, const b2linalg::Vector< double, b2linalg::Vdense_constref > &dofdotdot, const Field< double > &field, int face, b2linalg::Index &dof_numbering, b2linalg::Vector< double, b2linalg::Vdense > &discretised_field, b2linalg::Matrix< double, b2linalg::Mpacked > &d_discretised_field_d_dof)

virtual void	field_edge_integration (const b2linalg::Vector< double, b2linalg::Vdense_constref > &dof, const b2linalg::Vector< double, b2linalg::Vdense_constref > &dofdot, const b2linalg::Vector< double, b2linalg::Vdense_constref > &dofdotdot, const Field< double > &f, int edge, b2linalg::Index &dof_numbering, b2linalg::Vector< double, b2linalg::Vdense > &discretised_field, b2linalg::Matrix< double, b2linalg::Mpacked > &d_discretised_field_d_dof)

void	AssembleElementEffectiveSystem (const solver::TypedSolver< double, MATRIX_FORMAT > &solver, b2linalg::Matrix< double, typename MATRIX_FORMAT::dense > &k_eff, b2linalg::Vector< double, b2linalg::Vdense > &f_eff, b2linalg::Index &index)
	Assembles the element effective matrix and vector.

b2linalg::Matrix< double, typename MATRIX_FORMAT::dense >	AssembleElementEffectiveMatrix (const solver::TypedSolver< double, MATRIX_FORMAT > &solver, b2linalg::Index &index)
	Assembles the element effective matrix.

double	ComputeElementError (const solver::TypedSolver< double, MATRIX_FORMAT > &solver)
	TODO.

b2linalg::Vector< double, b2linalg::Vdense >	AssembleElementEffectiveVector (const solver::TypedSolver< double, MATRIX_FORMAT > &solver, b2linalg::Index index)
	Assembles the element effective vector.

void	ComputeElementGradient (const solver::TypedSolver< double, MATRIX_FORMAT > &solver, GradientContainer *const gradient_container)
	Compute the gradients for post-processing.

virtual void	get_value (Model &model, const bool linear, const EquilibriumSolution equilibrium_solution, const double time, const double delta_time, const b2linalg::Matrix< double, b2linalg::Mrectangle_constref > &dof, GradientContainer gradient_container, SolverHints solver_hints, b2linalg::Index &dof_numbering, b2linalg::Vector< double, b2linalg::Vdense > &value, const std::vector< bool > &d_value_d_dof_flags, std::vector< b2linalg::Matrix< double, b2linalg::Mrectangle > > &d_value_d_dof, b2linalg::Vector< double, b2linalg::Vdense > &d_value_d_time)
	Compute the internal forces and their derivatives of the element.

virtual void	get_value (Model &model, const bool linear, const EquilibriumSolution equilibrium_solution, const double time, const double delta_time, const b2linalg::Matrix< double, b2linalg::Mrectangle_constref > &dof, GradientContainer gradient_container, SolverHints solver_hints, b2linalg::Index &dof_numbering, b2linalg::Vector< double, b2linalg::Vdense > &value, const std::vector< bool > &d_value_d_dof_flags, std::vector< b2linalg::Matrix< double, b2linalg::Mpacked > > &d_value_d_dof, b2linalg::Vector< double, b2linalg::Vdense > &d_value_d_time)
	Compute the internal forces and their derivatives of the element.

Public Member Functions inherited from b2000::Element
virtual const std::vector< VariableInfo >	get_value_info () const

virtual std::pair< int, VariableInfo >	get_constraint_info ()

	Element (const bool symmetric)
	Sets the symmetry property of this element. Call this constructor if your newly defined element is not symmetric.

size_t	get_id () const

void	set_id (size_t id_)

int	get_elem_type_number () const

void	set_elem_type_number (int elno)

const std::string &	get_object_name () const override

void	set_object_name (const std::string &name)

int	get_number_of_dof () const override

size_t	set_global_dof_numbering (size_t index) override

std::pair< size_t, size_t >	get_global_dof_numbering () const override

virtual void	get_dof_numbering (b2linalg::Index &dof_numbering)

virtual void	set_additional_properties (const RTable &rtable)

virtual void	init (Model &model)

virtual std::pair< size_t, Element ** >	get_subelements ()

iterator	begin ()

iterator	end () const

virtual int	edge_field_order (const int edge_id, const std::string &field_name)

virtual bool	edge_field_linear_on_dof (const int edge_id, const std::string &field_name)

virtual int	edge_field_polynomial_sub_edge (const int edge_id, const std::string &field_name, b2linalg::Vector< double, b2linalg::Vdense > &sub_nodes)

virtual void	edge_geom (const int edge_id, const double internal_coor, b2linalg::Vector< double > &geom, b2linalg::Vector< double > &d_geom_d_icoor)

virtual int	face_field_polynomial_sub_face (const int face_id, const std::string &field_name, b2linalg::Matrix< double, b2linalg::Mrectangle > &sub_nodes, std::vector< Triangle > &sub_faces)

virtual int	body_field_order (const std::string &field_name)

virtual bool	body_field_linear_on_dof (const std::string &field_name)

virtual int	body_field_polynomial_sub_volume (const std::string &field_name, b2linalg::Matrix< double, b2linalg::Mrectangle > &sub_nodes, std::vector< Tetrahedral > &sub_volumes)

virtual void	body_geom (const b2linalg::Vector< double, b2linalg::Vdense_constref > &internal_coor, b2linalg::Vector< double > &geom, b2linalg::Matrix< double, b2linalg::Mrectangle > &d_geom_d_icoor)

virtual std::pair< int, int >	get_state_buffer_size () const

virtual std::pair< int , double >	get_state_buffer (std::pair< int , double > buffers) const

virtual std::pair< const int , const double >	set_state_buffer (std::pair< const int , const double > buffers)

Public Member Functions inherited from b2000::Object
virtual	~Object ()

Additional Inherited Members
Public Types inherited from b2000::Element
enum	VariableInfo { zero = 0 , constant = 1 , linear = 2 , nonlinear = 3 }

Static Public Attributes inherited from b2000::TypedElement< double >
static type_t	type
	Element interface typed by the dof.

Static Public Attributes inherited from b2000::Object
static ObjectType	type

Detailed Description

template<typename SHAPE, bool EULER>
class b2000::ElementCFD3DField< SHAPE, EULER >

CFD element for 3D analysis

Member Function Documentation

◆ face_field_linear_on_dof()

template<typename SHAPE , bool EULER>

bool b2000::ElementCFD3DField< SHAPE, EULER >::face_field_linear_on_dof	(	const int	face_id,
		const std::string &	field_name
	)

inlineoverridevirtual

Returns: Whether the integrated field depends linearly on the dof (true), or whether the field depends in nonlinear fashion (false).

Parameters

face_id	The face number, starting at 1, according to the B2000++ conventions.
field_name	The field name.

Solid elements depend linearly on the dof. Shell elements depend in nonlinear fashion on the side faces and on the lower and upper surfaces. Their reference surface depends linearly on the dof.

Reimplemented from b2000::Element.

◆ face_field_order()

template<typename SHAPE , bool EULER>

int b2000::ElementCFD3DField< SHAPE, EULER >::face_field_order	(	const int	face_id,
		const std::string &	field_name
	)

inlineoverridevirtual

Returns: The integration order over the element face. This depends on the shape functions. For instance, linear elements return 1, quadratic elements return 2, etc.

Parameters

face_id	The face number, starting at 1, according to the B2000++ conventions.
field_name	The field name. Usually, the integration order is independent of the field name.

Reimplemented from b2000::Element.

◆ face_geom()

template<typename SHAPE , bool EULER>

void b2000::ElementCFD3DField< SHAPE, EULER >::face_geom	(	const int	face_id,
		const b2linalg::Vector< double, b2linalg::Vdense_constref > &	internal_coor,
		b2linalg::Vector< double > &	geom,
		b2linalg::Matrix< double, b2linalg::Mrectangle > &	d_geom_d_icoor
	)

inlineoverridevirtual

Compute the face geometry at a specified integration point.

Parameters

face_id	The face number, starting at 1, according to the B2000++ conventions (input).
internal_coor	The face-internal coordinates of the integration point (input).
geom	The branch-global coordinates of the integration point (vector, output).
d_geom_d_icoor	The derivatives of the branch-global coordinates w.r.t. the internal coordinates; these are the tangential vectors (matrix, output).

Reimplemented from b2000::Element.

◆ get_nodes()

template<typename SHAPE , bool EULER>

std::pair< int, Node *const * > b2000::ElementCFD3DField< SHAPE, EULER >::get_nodes ( ) const

inlineoverridevirtual

Must be overridden.

Returns: A pair containing the number of nodes and a const pointer to an array of pointers to Node instances.

See also: set_nodes()

Reimplemented from b2000::Element.

◆ get_property()

template<typename SHAPE , bool EULER>

const ElementProperty * b2000::ElementCFD3DField< SHAPE, EULER >::get_property ( ) const

inlineoverridevirtual

Must be overridden. This implementation returns a 0-pointer.

Returns: A pointer to the element property.

See also: set_property()

Reimplemented from b2000::Element.

◆ set_nodes()

template<typename SHAPE , bool EULER>

void b2000::ElementCFD3DField< SHAPE, EULER >::set_nodes ( std::pair< int, Node *const * > nodes )

inlineoverridevirtual

Must be overridden. Set the nodes of the element. This function is called by the Domain object.

todo: not all elements override this function!

Parameters

nodes A pair (first, second) with first being the number of nodes, and second a pointer to a const array containing pointers to the Node instances. The length of the array corresponds to the number of nodes.

Implementations should verify the number of nodes. The checking of the correct node type is optional.

It is recommended to perform a dynamic cast to the exact node type, this allows the use of the get_coor_s() and get_dof_numbering_s() functions without the overhead of the virtual functions Node::get_coor() and DegreesOfFreedom::get_global_dof_numbering().

Example

A quadrilateral 2D stress element named Q4.S.2D.EXAMPLE. It has 3 coordinates (the Z-coordinate set to 0) and displacement dof's in X- and Y-direction. The following using statement gives the alias NodeType for the exact Node type:

using NodeType = HNode<Coor<Trans<X,Y,Z>>,Dof<DTrans<DX,DY>>>;

The class has an array of NodeType pointers.

private:

NodeType* nodes[4];

This array is initialized in the set_nodes() function:

void
set_nodes(std::pair<int, Node* const *> nodes_) {
    if (nodes_.first != 4)
        SizeError() << "Element " << get_object_name() << " of type "
            "Q4.S.2D.EXAMPLE needs 4 nodes." << THROW;
    for (int i = 0; i != 4; ++i) {
        nodes[i] = dynamic_cast<NodeType*>(nodes_.second[i]);
        if (nodes[i] == 0)
            TypeError() << "Element " << get_object_name() <<
                " of type Q4.S.2D.EXAMPLE requires nodes of type " <<
                typeid(NodeType) << "." << THROW;
    }
}

The code to retrieve the coordinates of all 4 nodes of the element, without the overhead of virtual functions could then look like this:

double coor[4][3];
for (int i = 0; i != 4; ++i)
    NodeType::get_coor_s(coor[i], nodes[i]);

The dof_numbering index vector (as given in the get_value(), get_constraint(), etc. functions) can then be initialized without the overhead of virtual functions:

{
    dof_numbering.resize(8);
    size_t* ptr = &dof_numbering[0];
    for (int i = 0; i != 4; ++i)
        ptr = NodeType::get_global_dof_numbering_s(ptr, nodes[i]);
}

Reimplemented from b2000::Element.

◆ set_property()

template<typename SHAPE , bool EULER>

void b2000::ElementCFD3DField< SHAPE, EULER >::set_property ( ElementProperty * property )

inlineoverridevirtual

Must be overridden. Set the element property. This function is called by the Domain object. This implementation performs no operation.

todo: not all elements override this function!

Parameters

property A pointer to the ElementProperty instance.

ElementProperty instances contain material parameters and additional element properties (such as the thickness for shell elements).

Depending on the type of element (stress, heat-transfer, 2D, 3D, etc.), the appropriate property class should be used. The Domain object constructs a "generic" property object for each defined material, and passes a pointer to the Element instance by calling the set_property() function. The implementation of set_property() must perform a dynamic cast to the desired property class and check the result.

Example

An element with 4 nodes and type Q4.S.2D.EXAMPLE. The class has a private variable for storing the pointer to the property. In this case, the element is a 2D stress element, and accordingly, the property class is SolidMechanics2D.

SolidMechanics2D* property;

b2000::SolidMechanics2D

Definition b2solid_mechanics.H:562

The function is implemented as follows:

void
set_property(ElementProperty* property_) {
    property = dynamic_cast<SolidMechanics2D*>(property_);
    if (property == 0)
        TypeError() <<
            "Element " << get_object_name() << " of type Q4.S.2D.EXAMPLE "
            "does not accept materials of type " << typeid(*property_) <<
            " but requires instead materials derived from the class " <<
            typeid(SolidMechanics2D) << "." << THROW;
}

Reimplemented from b2000::Element.

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

elements/smr/b2element_field_transfer.H

	b2api B2000++ API Reference Manual, VERSION 4.6

Public Member Functions

Additional Inherited Members

Detailed Description

Member Function Documentation

◆ face_field_linear_on_dof()

◆ face_field_order()

◆ face_geom()

◆ get_nodes()

◆ get_property()

◆ set_nodes()

Example

◆ set_property()

Example