A specialized version of DisplacementCubitBcData that includes an additional rotation offset. More...

#include <src/boundary_conditions/EssentialDisplacementCubitBcData.hpp>

Inheritance diagram for MoFEM::DisplacementCubitBcDataWithRotation:

Collaboration diagram for MoFEM::DisplacementCubitBcDataWithRotation:

Public Member Functions
	DisplacementCubitBcDataWithRotation ()

Public Member Functions inherited from MoFEM::DisplacementCubitBcData
struct	__attribute__ ((packed)) _data_
	attributes of DisplacementCubitBcData More...

std::size_t	getSizeOfData () const
	get data structure size More...

const void *	getDataPtr () const
	get pointer to data structure More...

	DisplacementCubitBcData ()

MoFEMErrorCode	fill_data (const std::vector< char > &bc_data)
	get data from structure More...

MoFEMErrorCode	set_data (void *tag_ptr, unsigned int size) const
	set data on structure More...

Public Member Functions inherited from MoFEM::GenericCubitBcData
virtual MoFEMErrorCode	fill_data (const std::vector< char > &bc_data)
	get data from structure More...

virtual MoFEMErrorCode	set_data (void *tag_ptr, unsigned int size) const
	set data on structure More...

virtual std::size_t	getSizeOfData () const =0
	get data structure size More...

virtual const void *	getDataPtr () const =0
	get pointer to data structure More...

virtual const CubitBCType &	getType () const
	get data type More...

	GenericCubitBcData (const CubitBCType type)

virtual	~GenericCubitBcData ()

Static Public Member Functions
static FTensor::Tensor1< double, 3 >	GetRotDisp (const FTensor::Tensor1< double, 3 > &angles, FTensor::Tensor1< double, 3 > coordinates, FTensor::Tensor1< double, 3 > offset=FTensor::Tensor1< double, 3 >{ 0., 0., 0.})
	Calculates the rotated displacement given the rotation angles, coordinates, and an optional offset. More...

Public Attributes
std::array< double, 3 >	rotOffset

Public Attributes inherited from MoFEM::DisplacementCubitBcData
_data_	data

Public Attributes inherited from MoFEM::GenericCubitBcData
const CubitBCType	tYpe
	Type of boundary condition. More...

Detailed Description

A specialized version of DisplacementCubitBcData that includes an additional rotation offset.

DisplacementCubitBcDataWithRotation extends the DisplacementCubitBcData structure with a rotation offset and a method to calculate the rotated displacement given the rotation angles and coordinates.

Definition at line 24 of file EssentialDisplacementCubitBcData.hpp.

Constructor & Destructor Documentation

◆ DisplacementCubitBcDataWithRotation()

MoFEM::DisplacementCubitBcDataWithRotation::DisplacementCubitBcDataWithRotation ( )

inline

Definition at line 27 of file EssentialDisplacementCubitBcData.hpp.

27: rotOffset{0, 0, 0} {}

MoFEM::DisplacementCubitBcDataWithRotation::rotOffset

std::array< double, 3 > rotOffset

Definition: EssentialDisplacementCubitBcData.hpp:26

Member Function Documentation

◆ GetRotDisp()

static FTensor::Tensor1< double, 3 > MoFEM::DisplacementCubitBcDataWithRotation::GetRotDisp	(	const FTensor::Tensor1< double, 3 > &	angles,
		FTensor::Tensor1< double, 3 >	coordinates,
		FTensor::Tensor1< double, 3 >	offset = `FTensor::Tensor1<double, 3>{ 0., 0., 0.}`
	)

inlinestatic

Calculates the rotated displacement given the rotation angles, coordinates, and an optional offset.

Parameters

angles	A FTensor::Tensor1 containing the rotation angles.
coordinates	A FTensor::Tensor1 containing the coordinates.
offset	An optional FTensor::Tensor1 containing the offset (default is {0., 0., 0.}).

Returns: FTensor::Tensor1<double, 3> representing the rotated displacement.

Definition at line 40 of file EssentialDisplacementCubitBcData.hpp.

                                                                         {
                 0., 0., 0.}) {
 
    FTensor::Index<'i', 3> i;
    FTensor::Index<'j', 3> j;
    FTensor::Index<'k', 3> k;
 
    FTensor::Tensor1<double, 3> rotated_displacement;
    rotated_displacement(i) = 0;
 
    auto get_rotation = [&](auto &omega) {
      FTensor::Tensor2<double, 3, 3> rotation_matrix;
 
      constexpr auto kronecker_delta = FTensor::Kronecker_Delta<int>();
      rotation_matrix(i, j) = kronecker_delta(i, j);
 
      const double angle = sqrt(omega(i) * omega(i));
      if (std::abs(angle) < std::numeric_limits<double>::epsilon())
        return rotation_matrix;
 
      FTensor::Tensor2<double, 3, 3> t_omega;
      t_omega(i, j) = FTensor::levi_civita<double>(i, j, k) * omega(k);
      const double a = sin(angle) / angle;
      const double sin_squared = sin(angle / 2.) * sin(angle / 2.);
      const double b = 2. * sin_squared / (angle * angle);
      rotation_matrix(i, j) += a * t_omega(i, j);
      rotation_matrix(i, j) += b * t_omega(i, k) * t_omega(k, j);
 
      return rotation_matrix;
    };
 
    auto rotation = get_rotation(angles);
    FTensor::Tensor1<double, 3> coordinate_distance;
    coordinate_distance(i) = offset(i) - coordinates(i);
    rotated_displacement(i) =
        coordinate_distance(i) - rotation(j, i) * coordinate_distance(j);
 
    return rotated_displacement;
  }

Member Data Documentation

◆ rotOffset

std::array<double, 3> MoFEM::DisplacementCubitBcDataWithRotation::rotOffset

Definition at line 26 of file EssentialDisplacementCubitBcData.hpp.

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

src/boundary_conditions/EssentialDisplacementCubitBcData.hpp

Public Member Functions

Static Public Member Functions